Methods for Computing Periodic Steady-State Jacob White
|
|
- Solomon Manning
- 6 years ago
- Views:
Transcription
1 Introduction to Simulation - Lecture 15 Methods for Computing Periodic Steady-State Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
2 Outline Periodic Steady-state problems Application examples and simple cases Finite-difference methods Formulating large matrices Shooting Methods State transition function Sensitivity matrix Matrix Free Approach
3 Periodic Steady-State Basics Basic Definition () dx t dt = F x() t + u { { () t state input Suppose the system has a periodic input T 2T 3T t Many Systems eventually respond periodically x( t+ T) = x( t) for t >> 0
4 Periodic Steady-State Basics Basic Definition Interesting Property If x satisfies a differential equation which has a unique solution for any initial condition dx t dt ( ()) u() t = F x t + Then if u is periodic with period T and ( + ) = x( t + T ) = x( t) for all t > t0 x t T x t for some t
5 Periodic Steady-State Basics Application Examples Swaying Bridge Periodic Input Wind Response Oscillating Platform Desired Info Oscillation Amplitude
6 Periodic Steady-State Basics Application Examples Communication Integrated Circuit Periodic Input Received Signal at 900Mhz Response filtered demodulated signal Desired Info Distortion
7 Periodic Steady-State Basics Application Examples Automobile Vibration Periodic Input Regularly Spaced Road Bumps Response Car Shakes Desired Info Shake amplitude
8 Periodic Steady-State Basics RLC Circuit Simple Example RLC Filter, Spring+Mass+Dashpot Spring-Mass-Dashpot Force Both Described by Second-Order ODE 2 d x dx M + D + x= u() t dt 2 dt { input
9 Periodic Steady-State Basics Simple Example RLC Filter, Spring+Mass+Dashpot Cont. Both Described by Second-Order ODE 2 d x dx M + D + x= u() t dt 2 dt u(t) = 0 lightly damped (D<<M) Response D t x() t Ke 2M cos φ + M
10 Periodic Steady-State Basics Simple Example RLC Filter, Spring+Mass+Dashpot Cont. Ke D 2M A lightly damped system oscillates many times before settling to a steady-state
11 Periodic Steady-State Basics Computing Steady State Frequency Domain Approach Sinusoidally excited linear time-invariant system dx t dt () e{ i t = Ax t + input Steady-State Solution simple to determine ω () = ( ) 1 i t x t iω A e ω Not useful for nonlinear or time-varying systems
12 Periodic Steady-State Basics Computing Steady State Time Integration Method Time-Integrate Until Steady-State Achieved dx t dt ( ()) l l 1 ˆ ˆ ( l ˆ ) = F x t + u() t x = x + t F x + u( l t) Need many timepoints for lightly damped case!
13 Aside Reviewing Integration Methods Solve with Backward-Euler Nonlinear System dx() t = F x() t + u( t) x { { ( 0 ) = x0 dt state input Initial Condition Backward Euler Equation for timestep l l l 1 l ( ) xˆ xˆ = t F xˆ + u( l t) How do we solve the backward-euler Equation?
14 Aside Reviewing Integration Methods Forward-Euler 1 ˆ 1 ( ) xt x = x(0) + tf x 0, u 0 xt xˆ = xˆ + tf xˆ, u t M L L 1 L 1 ( L) = +, ( ) L 1 xt xˆ xˆ tf xˆ u t Requires just function Evaluations Implicit Methods Backward-Euler Example Backward-Euler ( ˆ ) xt xˆ = x(0) + tf x, u t xt xˆ = xˆ + tf xˆ, u t L L 1 L ˆ ˆ ˆ L xt x = x + tf x, u t L M Nonlinear equation solution at each step Stepwise Nonlinear equation solution needed whenever β0 0
15 Aside Reviewing Integration Methods Implicit Methods Solution with Newton Rewrite the multistep Equation ( t ) k ( l ) l j l j β l + α jx β j tl j l α ˆ ˆ ˆ ˆ 0x t 0f x, u t f x, u = 0 Solve with Newton l, j f xˆ, u( t ) j= 1 j= 1 k b Independent of xˆ l l l, j+ 1 l, j l, j l, j α0i tβ 0 ( xˆ xˆ ) = α ˆ 0x tβ0f ( xˆ, u( tl )) + b x Jacobian ( l, j F x ) Here j is the Newton iteration index
16 Aside Reviewing Integration Methods Implicit Methods Solution with Newton Cont. Newton Iteration: l, j ( ˆ, ) f x u t l 0 0 ( l j + α I tβ x x l j ) = F x l j x, 1,, ˆ ˆ t l k Solution with Newton is very efficient xˆ l t l 3 t l 2 tl 1 t l l, j ( ˆ, ( tl )),0 xˆ l Converged Solution Polynomial Predictor f x u α0i tβ0 α0i as t 0 x Easy to generate a good initial guess using polynomial fitting Jacobian become easy to factor for small timesteps
17 Boundary-Value Problem Basic Formulation Periodicity Constraint Differential Equation Solution d N Differential Equations: x i t = F i x t dt () () N Periodicity Cons traints: x T = x 0 i i
18 Boundary-Value Problem Finite Difference Methods Linear Example Problem dx t dt = Ax t + u t t 0, T x T = x t () () { [ ] input Discretize with Backward-Euler 1 0 ( 1 xˆ = xˆ + t Axˆ + u( t) ) xˆ = xˆ + t Axˆ + u 2 t L L 1 L ˆ ˆ ˆ ( ) ( ) x = x + t Ax + u L t t = Periodicity implies xˆ = xˆ L periodicity constraint T L 0
19 Boundary-Value Problem Finite Difference Methods Linear Example Matrix Form NxL NxL I A 0 0 I t t xˆ u ( 0) xˆ u( t) I I A 0 0 t t M = M 0 0 O O M M 1 1 L 0 0 I I A xˆ u( L t) t t Matrix is almost lower triangular
20 Boundary-Value Problem Finite Difference Methods Nonlinear Problem dx t dt ( ()) () { [ 0, ] = F x t + u t t T x T = x t input periodicity constraint Discretize with Backward-Euler 1 L xˆ xˆ t( F( xˆ 1 ) + u( t) ) 1 xˆ xˆ 2 xˆ 1 t( F( xˆ 2 ) + u( 2 t) ) 2 xˆ H FD = L ˆ L L 1 L x xˆ xˆ t F xˆ + u L t = 0 M ( ) Solve Using Newton s Method
21 Boundary-Value Problem Shooting Method Basic Definitions dx( t) Start with = F ( x() t ) + u() t dt And assume x(t) is unique given x(0). D.E. defines a State-Transition Function Φ ( yt,, t) x( t) where xt is the D.E. solution given x t = y 0
22 Boundary-Value Problem Shooting Method State Transition function Example dx t dt = λx t () Φ (,, ) ( t t ) yt t e 1 0 y 0 1 λ
23 Boundary-Value Problem Solve ( ) Use Newton s method Shooting Method Abstract Formulation H x 0 = Φ x 0,0, T x = 0 xt H ( x,0, T) Φ = J x I H x ( k )( k+ 1 k ) = ( k ) J x x x H x
24 Boundary-Value Problem Shooting Method Computing Newton ( T) To Compute Φ x 0,0, dx( t) Integrate = F( x() t ) + u() t on [0,T] dt (,0, ) Φ x T What is? x x ( 0) + ε x ( 0) x ε ( T) xt Indicates the sensitivity of x(t) to changes in x(0)
25 Boundary-Value Problem Φ x,0, T x Shooting Method Sensitivity Matrix by Perturbation ε1 x T x T x T x T ε N ε 1 ε1 x T x T x T x T ε N N N N N ε 1 L L M L L M M L L M L L ε ε N N
26 Boundary-Value Problem Shooting Method Efficient Sensitivity Evaluation Differentiate the first step of Backward-Euler x ( 0) ( 1 ( 0 ) ( ( 1 ) ) 0 ) xˆ x t F xˆ + u t = ( 0) ( 0) ( ˆ 1) xˆ 1 x 0 F x xˆ 1 t = x x x x I ( ˆ ) ( 0) 1 1 x 0 F x xˆ t = x x( 0) x 0 0 I
27 Boundary-Value Problem Shooting Method Efficient Sensitivity Matrix Cont Applying the same trick on the l-th step I ( ˆ ) l l l 1 F x xˆ xˆ t = x x x ( 0) ( 0) Φ x,0, T x L l= 1 I t ( ˆl ) F x x 1
28 Boundary-Value Problem Shooting Method Observations on Sensitivity Matrix Newton at each timestep uses same matrices Φ x,0, T x L l= 1 ( ˆ l ) F x I t x Timestep Newton Jacobian Formula simplifies in the linear case Φ x,0, T x ( I ta) L 1
29 Shooting Method Matrix-Free Approach Basic Setup Start with dx t dt ( ()) ut () = F xt + Use Newton s method ( ) H x 0 = Φ x 0,0, T x 0 = 0 H ( x,0, T) Φ = J x I H x ( k )( k+ 1 k ) = ( k ) J x x x H x
30 Shooting Method Matrix-Free Approach Matrix-Vector Product Solve Newton equation with Krylov-subspace method ( k x,0, T) Φ I x x = x Φ x T x x b A ( k+ 1 k) k ( k,0, ) Matrix-Vector Product Computation ( x k,0, T) ( x k ε p j,0, T) ( x k,0, T) Φ Φ + Φ x ε j j I p p Krylov method search direction
31 Shooting Method Matrix-Free Approach Convergence for GCR Example dx Ax dt = Shooting-Newton Jacobian 0 eig A real and negative Φ x,0, T I = e AT I x
32 Shooting Method Matrix-Free Approach Convergence for GCR-evals λ1t e 1 AT e I = S O S λnt e 1 1 Many Fast Modes cluster at 1 1 Few Slow Modes larger than 1
33 Summary Periodic Steady-state problems Application examples and simple cases Finite-difference methods Formulating large matrices Shooting Methods State transition function Sensitivity matrix
Introduction to Simulation - Lecture 9. Multidimensional Newton Methods. Jacob White
Introduction to Simulation - Lecture 9 Multidimensional Newton Methods Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Michal Rewienski, and Karen Veroy Outline Quick Review of -D Newton Convergence
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationFast Methods for Integral Equations. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simulation - Lecture 23 Fast Methods for Integral Equations Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Solving Discretized Integral Equations Using
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationAdvanced Computational Methods for VLSI Systems. Lecture 4 RF Circuit Simulation Methods. Zhuo Feng
Advanced Computational Methods for VLSI Systems Lecture 4 RF Circuit Simulation Methods Zhuo Feng 6. Z. Feng MTU EE59 Neither ac analysis nor pole / zero analysis allow nonlinearities Harmonic balance
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationTransduction Based on Changes in the Energy Stored in an Electrical Field
Lecture 6- Transduction Based on Changes in the Energy Stored in an Electrical Field Actuator Examples Microgrippers Normal force driving In-plane force driving» Comb-drive device F = εav d 1 ε oε F rwv
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationVibrations: Second Order Systems with One Degree of Freedom, Free Response
Single Degree of Freedom System 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 5//007 Lecture 0 Vibrations: Second Order Systems with One Degree of Freedom, Free Response Single
More information4.2 Homogeneous Linear Equations
4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this
More informationNumerical Algorithms for ODEs/DAEs (Transient Analysis)
Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no
More informationMath 240: Spring/Mass Systems II
Math 240: Spring/Mass Systems II Ryan Blair University of Pennsylvania Monday, March 26, 2012 Ryan Blair (U Penn) Math 240: Spring/Mass Systems II Monday, March 26, 2012 1 / 12 Outline 1 Today s Goals
More informationComputing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices
Computing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices Alper Demir David Long Jaijeet Roychowdhury Bell Laboratories Murray Hill New Jersey USA Abstract The main effort in oscillator
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationIntroduction to Simulation - Lecture 16. Methods for Computing Periodic Steady-State - Part II Jacob White
Inroducion o Simuaion - ecure 6 Mehods for Compuing Periodic Seady-Sae - Par II Jacob Whie hanks o Deepak Ramaswamy, Micha Rewienski, and Karen Veroy Ouine hree Mehods so far ime inegraion uni seady-sae
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More informationMath Assignment 5
Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationIntroduction to Simulation - Lecture 10. Modified Newton Methods. Jacob White
Introduction to Simulation - Lecture 0 Modified Newton Methods Jacob White Thans to Deepa Ramaswamy, Jaime Peraire, Michal Rewiensi, and Karen Veroy Outline Damped Newton Schemes SMA-HPC 003 MIT Globally
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationChapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs
Chapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs In the previous chapter, we investigated stiffness in ODEs. Recall that an ODE is stiff if it exhibits behavior on widelyvarying timescales.
More informationMath 211. Substitute Lecture. November 20, 2000
1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationMath 240: Spring-mass Systems
Math 240: Spring-mass Systems Ryan Blair University of Pennsylvania Tuesday March 1, 2011 Ryan Blair (U Penn) Math 240: Spring-mass Systems Tuesday March 1, 2011 1 / 15 Outline 1 Review 2 Today s Goals
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More information6.003 Homework #6. Problems. Due at the beginning of recitation on October 19, 2011.
6.003 Homework #6 Due at the beginning of recitation on October 19, 2011. Problems 1. Maximum gain For each of the following systems, find the frequency ω m for which the magnitude of the gain is greatest.
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationLaplace s Equation FEM Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera
Introduction to Simulation - Lecture 19 Laplace s Equation FEM Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera Outline for Poisson Equation
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationLecture 9: Harmonic Loads (Con t)
Lecture 9: Harmonic Loads (Con t) Reading materials: Sections 3.4, 3.5, 3.6 and 3.7 1. Resonance The dynamic load magnification factor (DLF) The peak dynamic magnification occurs near r=1 for small damping
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationMass on a Horizontal Spring
Course- B.Sc. Applied Physical Science (Computer Science) Year- IInd, Sem- IVth Subject Physics Paper- XIVth, Electromagnetic Theory Lecture No. 22, Simple Harmonic Motion Introduction Hello friends in
More informationThe Direct Transcription Method For Optimal Control. Part 2: Optimal Control
The Direct Transcription Method For Optimal Control Part 2: Optimal Control John T Betts Partner, Applied Mathematical Analysis, LLC 1 Fundamental Principle of Transcription Methods Transcription Method
More informationIn this lecture you will learn the following
Module 9 : Forced Vibration with Harmonic Excitation; Undamped Systems and resonance; Viscously Damped Systems; Frequency Response Characteristics and Phase Lag; Systems with Base Excitation; Transmissibility
More informationEven-Numbered Homework Solutions
-6 Even-Numbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complex-valued solution Y
More informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationModeling and Experimentation: Mass-Spring-Damper System Dynamics
Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationEE 230 Lecture 24. Waveform Generators. - Sinusoidal Oscillators
EE 230 Lecture 24 Waveform Generators - Sinusoidal Oscillators Quiz 18 Determine the characteristic equation for the following network without adding an excitation. C R L And the number is? 1 3 8 2? 6
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationQuiescent Steady State (DC) Analysis The Newton-Raphson Method
Quiescent Steady State (DC) Analysis The Newton-Raphson Method J. Roychowdhury, University of California at Berkeley Slide 1 Solving the System's DAEs DAEs: many types of solutions useful DC steady state:
More informationDifferential and Difference LTI systems
Signals and Systems Lecture: 6 Differential and Difference LTI systems Differential and difference linear time-invariant (LTI) systems constitute an extremely important class of systems in engineering.
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationDynamics of structures
Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 One degree of freedom systems in real life 2 1 Reduction of a system to a one dof system Example
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 State Estimator In previous section, we have discussed the state feedback, based on the assumption that all state variables are
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationForced Damped Vibrations
Forced Damped Vibrations Forced Damped Motion Definitions Visualization Cafe door Pet door Damped Free Oscillation Model Tuning a Dampener Bicycle trailer Forced Damped Motion Real systems do not exhibit
More informationCyber-Physical Systems Modeling and Simulation of Continuous Systems
Cyber-Physical Systems Modeling and Simulation of Continuous Systems Matthias Althoff TU München 29. May 2015 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May 2015 1 / 38 Ordinary Differential
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationGeneral Response of Second Order System
General Response of Second Order System Slide 1 Learning Objectives Learn to analyze a general second order system and to obtain the general solution Identify the over-damped, under-damped, and critically
More informationRead textbook CHAPTER 1.4, Apps B&D
Lecture 2 Read textbook CHAPTER 1.4, Apps B&D Today: Derive EOMs & Linearization undamental equation of motion for mass-springdamper system (1DO). Linear and nonlinear system. Examples of derivation of
More information18.03SC Practice Problems 14
1.03SC Practice Problems 1 Frequency response Solution suggestions In this problem session we will work with a second order mass-spring-dashpot system driven by a force F ext acting directly on the mass:
More information1 Oscillations MEI Conference 2009
1 Oscillations MEI Conference 2009 Some Background Information There is a film clip you can get from Youtube of the Tacoma Narrows Bridge called Galloping Gertie. This shows vibrations in the bridge increasing
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationIntroduction to Simulation - Lecture 2. Equation Formulation Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simulation - Lecture Equation Formulation Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Formulating Equations rom Schematics Struts and Joints
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
More informationChapter 8: Frequency Domain Analysis
Chapter 8: Frequency Domain Analysis Samantha Ramirez Preview Questions 1. What is the steady-state response of a linear system excited by a cyclic or oscillatory input? 2. How does one characterize the
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationMath Lab 9 3 November 2016 unid:
Math 2250-004 Lab 9 3 November 2016 unid: Name: Instructions and due date: Due: 10 November 2016 at the start of lab. If extra paper is necessary, please staple it at the end of the packet. For full credit:
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts Signals A signal is a pattern of variation of a physical quantity as a function of time, space, distance, position, temperature, pressure, etc. These quantities are usually
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationMAE 143B - Homework 7
MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationThe Harmonic Balance Method
For Nonlinear Microwave Circuits Hans-Dieter Lang, Xingqi Zhang Thursday, April 25, 2013 ECE 1254 Modeling of Multiphysics Systems Course Project Presentation University of Toronto Contents Balancing the
More information6.003 Homework #6 Solutions
6.3 Homework #6 Solutions Problems. Maximum gain For each of the following systems, find the frequency ω m for which the magnitude of the gain is greatest. a. + s + s ω m = w This system has poles at s
More informationLinear Systems. Chapter Basic Definitions
Chapter 5 Linear Systems Few physical elements display truly linear characteristics. For example the relation between force on a spring and displacement of the spring is always nonlinear to some degree.
More informationIntroduction to Controls
EE 474 Review Exam 1 Name Answer each of the questions. Show your work. Note were essay-type answers are requested. Answer with complete sentences. Incomplete sentences will count heavily against the grade.
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationLecture 3: Hamilton-Jacobi-Bellman Equations. Distributional Macroeconomics. Benjamin Moll. Part II of ECON Harvard University, Spring
Lecture 3: Hamilton-Jacobi-Bellman Equations Distributional Macroeconomics Part II of ECON 2149 Benjamin Moll Harvard University, Spring 2018 1 Outline 1. Hamilton-Jacobi-Bellman equations in deterministic
More information