Time Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors
|
|
- Arron Ryan
- 5 years ago
- Views:
Transcription
1 Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization of nonlinear models LTI System Response Phase-plane plots Copyright 217 by Robert Stengel. All rights reserved. For educational use only. 1 Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors x = x y ; v = v x v y ; a = a x a y ; j = j x j y ; s = s x s y ( t) ( t) x x t v x v x a x a x 1 x 1 t t 2 2 t 3 6 t 4 24 t 5 12 v x ( ) 1 a x ( ) = y1( ) t t 2 2 t 3 6 t 4 24 j x ( ) 1 y ( t) 1 s x ( ) 1 1 t t 2 2 t 3 6 t 4 24 t 5 12 t t 2 2 t 3 6 v y ( ) 1 c x ( ) v y ( t = ) 1 t t 2 2 t 3 6 t 4 24 a y ( ) 1 1 t t a y ( t) 2 2 t 3 6 y v y a y j y s y c y 2
2 Two-Dimensional Trajectory Solve for Cartesian components separately x Component j x s x c x 6 t 3 6 t 3 36 t 2 24 t 2 9 t 3 t = 36 t 4 36 t t t 3 36 t 2 24 t 2 72 t 5 72 t 5 36 t 4 36 t 4 6 t 3 6 t 3 ( t) ( t) x x t v x v x a x a x y Component j y s y c y 6 t 3 6 t 3 36 t 2 24 t 2 9 t 3 t = 36 t 4 36 t t t 3 36 t 2 24 t 2 72 t 5 72 t 5 36 t 4 36 t 4 6 t 3 6 t 3 ( t) ( t) y y t v y v y a y a y 3 Two-Dimensional Example Required acceleration vector is specified by a( t) = a( ) + j( )t + s( )t ct 3 6 = a control ( t) + a gravity ( t) + a disturbance ( t) + f gravity ( t) + f disturbance ( t) = f control t m t 4
3 Six-Degree-of-Freedom (Rigid Body) Equations of Motion r I = H B I v B v B = 1 m f B B v B = L I B B B = I 1 B ( m B B I B B ) Translational position and velocity r I = x y z ; v B = u v w Rotational position and velocity ( p * = *; + B = q * ) r ( * * * ) 5 nonlinear, and complex Rate of change of Translational Position x I = ( cos cos )u + ( cos sin + sin sin cos )v + ( sin sin + cos sin cos )w y I = ( cos sin )u + ( cos cos + sin sin sin )v + (sin cos + cos sin sin )w z I = (sin )u + ( sin cos )v + ( cos cos )w Rate of change of Translational Velocity u = X / m gsin + rv qw v = Y / m + gsin cos ru + pw w = Z / m + gcos cos + qu pv Rate of change of Angular Position = p + ( qsin + r cos )tan = qcos rsin = ( qsin + r cos )sec Rate of change of Angular Velocity 2 I xx I zz I xz p = I zz L + I xz N { I xz ( I yy I xx I zz ) p + I 2 xz + I zz ( I zz I yy ) r }q pr I xz ( p 2 r 2 ) q = M I xx I zz I yy r = I xz L + I xx N { I xz ( I yy I xx I zz )r + I 2 xz + I xx ( I xx I yy ) p }q 2 I xx I zz I xz 6
4 Multiple Rigid Links Lead to Multiple Constraints Each link is subject to the same 6-DOF rigidbody dynamic equations but each link is constrained to have a single degree of freedom w.r.t. proximal link 7 Newton-Euler Link Dynamics Link dynamics are coupled Proximal-link loads affected by distal-link positions and velocities Distal-link accelerations affected by proximallink motions Joints produce constraints on link motions Net forces and torques at each joint related to velocities and accelerations of the centroids of the links Equations of motion derived directly for each link, with constraints v B = 1 m f B B v B B = I 1 B ( m B B I B B ) x = f [ x,u,w,p,t ] 8
5 Lagrangian Link Dynamics Lagrange s equation derives from Newton s Laws Principle of virtual work D Alembert s principle Dynamic behavior described by work done and energy stored in the system Equations of motion derived from Lagrangian function and Lagrange s equation Generalized coordinate, q n, e.g., ( x, y,z,,,) Generalized force, F n, e.g., f x, f y, f z,m x,m y,m z d dt L( q n, q n ) KE PE dq n dl q n, q n ( L q, q n n = F n ( q n 9 x = f [ x,u,w,p,t ] Hamiltonian Link Dynamics Hamilton s Principle: Lagrange s equation is a necessary condition for an extremum (i.e., maximum or minimum) extremum I = t 2 t 1 L( q n, q n )dt Hamiltonian function H ( p,q) = q n p n Equations of motion derived from an optimization problem L q n, q n Hamilton s equations Generalized momentum, p n, q n = H ( p,q) p n p n = H p,q q n e.g., mv x,mv y,mv z, I xx p, I yy q, I zz r x = f [ x,u,w,p,t ] 1
6 Example: Two-Link Robot Equations of Motion State vector x 1 x 2 x 3 x 4 = Angleof1 st link, 1,rad Angular rateof1 st link, 1,rad/sec Angleof 2 nd link, 2,rad Angular rateof 2 nd link, 2,rad/sec Control vector u 1 u 2 st = 1, torque at 1 joint 2, torque at 2 nd joint McKerrow, When Possible, Simplify the Equations Two-link robot equations of motion Mass, m, located at end of Link 2 Inertias of Links 1 and 2 neglected x 1 = x 2 x 2 = x 3 = x 4 1 x cos 2 2 x 4 sin 2x 3 + u 1 x 3 ml 2 2 x 4 = g l 2 cos x 3 x 2 2 sin 2x 3 + u 2 ml 2 2 from Example 7.1, Lagrangian derivation, McKerrow, 1991, pp
7 Differential Equations Integrated to Produce Time Response x = f [ x,u,w,p,t ], x() given x = x() + t f [ x,u,w,p,t ]dt Numerical integration is an approximation 13 Rectangular and Trapezoidal Integration of Differential Equations Rectangular (Euler) Integration x(t k ) = x(t k1 ) + x(t k1,t k ) [ ] t, t = t k t k1 x(t k1 ) + f x(t k1 ),u(t k1 ),w(t k1 ) Trapezoidal (modified Euler) Integration (ode23) x(t k ) x(t k1 ) + 1 [ 2 x + x 1 2 ] where [ ] t { } t x 1 = f x(t k1 ),u(t k1 ),w(t k1 ) x 2 = f [ x(t k1 ) + x 1 ],u(t k ),w(t k ) ode23 varies step size, t, to reduce numerical error 14
8 Numerical Integration Examples Cosine Wave 1-sec Step How can approximation accuracy be improved? 15 More Complicated Algorithms (e.g., MATLAB) 16
9 1-D Example a net ( t) = ( ) +.6t.36t t Comparison of Exact and Numerically Integrated Trajectories x 1 t Calculate trajectory, given constants for t f = 1 x t v t a t 1 t t 2 2 t 3 6 t 4 24 t 5 12 = 1 t t 2 2 t 3 6 t t t 2 2 t Calculate trajectory by numerical integration = x 2 ( t) = v ( t) = a( t) =.6t.36t t 3 6 x 2 t 18
10 Comparison of Exact and Numerically Integrated Trajectories Polynomial, ode23, and ode45 19 Comparison of Exact and Numerically Integrated Trajectories (Zoom) Polynomial, ode23, and ode45 2
11 Generic Dynamic System Dynamic Process: Current state may depend on prior state x : state dim = (n x 1) u : input dim = (m x 1) w : disturbance dim = (s x 1) p : parameter dim = (l x 1) t : time (independent variable, 1 x 1) Observation Process: Measurement may contain error or be incomplete y : output (error-free) dim = (r x 1) n : measurement error dim = (r x 1) z : measurement dim = (r x 1) 21 Equations of the System x = dx dt Dynamic Equation = f[x,u,w,p,t] Output Equation y = h[x,u] Measurement Equation z = y + n 22
12 Dynamic System Example: Automotive Vehicle Dynamic Process x : dynamic state Position, velocity, angle, angular rate u : input Steering, throttle, brakes w : disturbance Road surface, wind p : parameter Weight, moments of inertia, drag coefficient, spring constants t : time (independent variable) Observation Process y : error-free output Speed, front-wheel angle, engine rpm, acceleration, yaw rate, throttle, brakes, GPS location n : measurement error Perturbations to y z : measurement Sum of y and n 23 Nonlinearity and Time Variation in Dynamic Systems Nonlinear, time-varying dynamics x = f x,u,w,p,t [ ] Nonlinear, time-invariant dynamics x = f x,u,w [ ] Linear, time-varying dynamics x = Fx + Gu + Lw Linear, time-invariant dynamics x = Fx + Gu + Lw NTV NTI LTV LTI 24
13 Nonlinearity and Time Variation in Dynamic Systems Nonlinear, time-varying dynamics Nonlinear, time-invariant dynamics Linear, time-varying dynamics Linear, time-invariant dynamics m, k, B 25 Solutions of Ordinary Differential Dynamic Equations Time-Domain Model (ODE) Solution by Numerical Integration Principle of Superposition Frequency-Domain Model. Nonlinear, timevarying (NTV) Yes No No. Nonlinear, timeinvariant (NTI) Yes No Yes (amplitudedependent, harmonics). Linear, timevarying (LTV) Yes Yes Approximate. Linear, timeinvariant (LTI) Yes Yes Yes 26
14 Comparison of Damped Linear and Nonlinear Systems x 1 = x 2 x 2 = 1x 1 x 2 Linear plus Stiffening Cubic Spring x 1 = x 2 x 2 = 1x 1 1x 1 3 x 2 Linear plus Weakening Cubic Spring x 1 = x 2 Linear Spring Spring Damper x 2 = 1x 1 +.8x 1 3 x 2 Spring Force vs. Displacement NTV or NTI? 27 MATLAB Simulation of Linear and Nonlinear Dynamic Systems MATLAB Main Program Nonlinear and Linear Examples clear tspan = [ 1]; xo = [, 1]; [t1,x1 = ode23(nonlin,tspan,xo); xo = [, 1]; [t2,x2] = ode23(nonlin,tspan,xo); xo = [, 1]; [t3,x3] = ode23(lin,tspan,xo); xo = [, 1]; [t4,x4] = ode23(lin,tspan,xo); subplot(2,1,1) plot(t1,x1(:,1),k,t2,x2(:,1),b,t3,x3(:,1),r,t4,x4(:,1),g) ylabel(position), grid subplot(2,1,2) plot(t1,x1(:,2),k,t2,x2(:,2),b,t3,x3(:,2),r,t4,x4(:,2),g) xlabel(time), ylabel(rate), grid Linear Spring x 1 = x 2 x 2 = 1x 1 x 2 function xdot = Lin(t,x) Linear Ordinary Differential Equation x(1) = Position x(2) = Rate xdot = [x(2) -1*x(1) - x(2)]; Weakening Spring x 1 = x 2 x 2 = 1x 1 +.8x 3 1 x 2 function xdot = NonLin(t,x) Nonlinear Ordinary Differential Equation x(1) = Position x(2) = Rate xdot = [x(2) -1*x(1) +.8*x(1)^3 - x(2)]; 28
15 Linear and Stiffening Cubic Springs: Small and Large Initial Rates x 1 () x 2 () = 1 and 1 Linear and nonlinear responses are indistinguishable with small initial condition 29 Linear and Weakening Cubic Springs: Small and Large Initial Rates x 1 () x 2 () = 1 and 1 3
16 Linearization of Nonlinear Equations Given Nominal (or reference) robot trajectory, control, and disturbance histories x N, u N, w N for t in [t o,t f ] Actual path and inputs x, u, w for t in [t o,t f ] = n 1 = m 1 = s 1 dim x dim u dim w perturbed by Initial condition variation Control variation Disturbance variation 31 State, Control, and Disturbance Perturbations Difference between nominal and actual paths: x(t o ) = x(t o ) x N (t o ) x = x x N Difference between nominal and actual inputs: u = u u N [ Control perturbation] w = w w N Disturbance perturbation [ ] 32
17 Linear Approximation of Perturbation Effects Nominal and actual paths both satisfy same dynamic equations Actual dynamics expressed as sum of nominal terms plus perturbations x = x N + x x N = f[x N,u N,w N,t] x = f[x,u,w,t] { } = f [x N + x],[u N + u],[w N + w,t] Exact f[x N,u N,w N,t]+ f ( x t)x + f ( u t)u + f ( w t)w Approx. Partial-derivative (Jacobian) matrices are evaluated along the nominal path 33 Linearized Equation Approximates Perturbation Dynamics Solve the nominal and perturbation parts separately Nominal (nonlinear) equation x N = f[x N,u N,w N,t], x N () given Perturbation (linear) equation x Fx + Gu + Lw, x() given Approximate total solution x x N + x 34
18 Jacobian Matrices Express Sensitivity to Small Perturbations Along Nominal Trajectory F = f x x=x N (t ) u=u N (t ) w=w N (t ) ; G = f u x=x N (t ) u=u N (t ) w=w N (t ) ; L = f w x=x N (t ) u=u N (t ) w=w N (t ) Sensitivity to state perturbations: stability matrix F = f = x x=x N (t ) u=u N (t ) w=w N (t ) f 1 x 1 f 1 x 2... f 1 x n f 2 f 2 f... 2 x 1 x 2 x n f n x 1 f n x 2... f n x n x=xn (t ) u=u N (t ) w=w N (t ) 35 Sensitivity to Small Control and Disturbance Perturbations Along Nominal Trajectory F = f x x=x N (t ) u=u N (t ) w=w N (t ) ; G = f u x=x N (t ) u=u N (t ) w=w N (t ) ; L = f w x=x N (t ) u=u N (t ) w=w N (t ) Control-effect matrix Disturbance-effect matrix G = f u x=x N (t ) u=u N (t ) w=w N (t ) = f 1 u m f n f... n u 1 u m f 1 u 1... x=xn (t ) u=u N (t ) w=w N (t ) L = f w x=x N (t ) u=u N (t ) w=w N (t ) = f 1 w 1... f 1 w s f n w 1... f n w s x=xn (t ) u=u N (t ) w=w N (t ) 36
19 Stiffening Cubic Spring Example Nonlinear equation, no inputs x 1 = f 1 = x 2 x 2 = f 2 = 1x 1 1x 3 1 x 2 Integrate nonlinear equation to produce nominal path x 1 () x 2 () t f ( f 1N f 2N dt x 1N x 2N in,t f 37 Stiffening Cubic Spring Example Evaluate partial derivatives along the path F( t) = f 1 x 1 = f 2 = 1 3x 2 1N x 1 f 1 x 2 = 1 f 2 x 2 = 1 ( ( ( ( ( G( t) = f 1 u = f 2 u = L( t) = f 1 w = f 2 w = 38
20 Nominal and Perturbation Dynamic Equations Nominal system, NTI x N = f[x N ], x N () given x 1N = x 2N x 2N = 1x 1N 1x 3 1N x 2N x 1N () x 2N () = 9 x 1 x 2 New initial condition: x = Fx, Perturbation system, LTV x() given = 1 ( 1 + 3x 2 1N (1 x 1 x 2 x 1 () x 2 () = 1 x 1 () x 2 () = 1 x 1N generated by nominal solution 39 Nominal and Perturbation Dynamic Solutions for Cubic Spring Example with x N () = If nominal solution remains at equilibrium x N = f[x N ], x N () =, x N = in [, ] Linearization is time-invariant x 1 x 2 = 1 x 1 (1 (1 x 2 = Linear, Time-Invariant (LTI) System 4
21 Comparison of Approximate and Exact Solutions Rate Initial Condition: x 2N () = 9 x 2 () = 1 x 2 = 1 x 2N + x 2 = 1 x N : Nominal x : Perturbation x : Actual x N + x : Approximation 41 Initial-Condition Response of a Linear, Time-Invariant (LTI) Model Doubling the initial condition doubles the output Stability, speed of response, and damping are independent of the initial condition Linear Model - Initial Condition F = [ ;.7814 ]; G = [1-1; 2]; Hx = [1 ; 1]; sys = ss(f, G, Hx,); xo = [1;]; [y1,t1,x1] = initial(sys, xo); xo = [2;]; [y2,t2,x2] = initial(sys, xo); plot(t1,y1,t2,y2) figure xo = [;1]; initial(sys, xo) F = G = 2 H x = I 2 ; H u = 42
22 Step Response of a Linear, Time-Invariant Model Doubling the step input doubles the output Stability, speed of response, and damping are independent of the input Linear Model - Step F = [ ;.7814 ]; G = [1-1; 2]; Hx = [1 ; 1]; sys = ss(f, G, Hx,); sys2 = ss(f, 2*G, Hx,); Step response step(sys, sys2) 43 Response to Combined Initial Condition and Step Input Linear system responses are additive Linear Model - Superposition F = [ ;.7814 ]; G = [1;]; Hx = [1 ; 1]; sys = ss(f, G, Hx,); xo = [1; ]; t = [:.2:2]; u = ones(1,length); [y1,t1,x1] = lsim(sys,u,t,xo); [y2,t2,x2] = lsim(sys,u,t); u = zeros(1,length); [y3,t3,x3] = lsim(sys,u,t,xo); plot(t1,y1,t2,y2,t3,y3) 44
23 Initial-Condition Responses of 1 st - Order LTI Systems are Exponentials x = Fx, x() given State vector is a scalar F = [ a] x = ax x() given 1 st -order system Unstable, a > Stable, a < t x = xdt = axdt = e at x() t 45 Initial-Condition Responses of 2 nd -Order LTI Systems Exponentials and Sinusoids x = Fx x() given F = f 11 f 12 f 21 f 22 x 1 = f 11 x 1 + f 12 x 2 + g 1 u x 2 = f 21 x 1 + f 22 x 2 + g 2 u Stable, Lightly Damped 2 nd -order system Sinusoid with exponential envelope x 1 = A 1 e nt cos n 1 2 t + 1 ( ) x 2 = A 2 e nt cos n 1 2 t + 2 ( ) 46
24 Equilibrium Response of Linear, Time-Invariant Models General equation x = Fx + Gu + Lw At equilibrium, Derivative goes to zero State is unchanging = Fx + Gu + Lw State at equilibrium ( Gu*+Lw *) x* = F 1 Gu*+Lw * = Adj F det F 47 State ( Phase )-Plane Plots Cross-plot of one component against another Time or frequency not shown explicitly in phase plane 2nd-Order Model - Initial Condition Response clear z =.1; Damping ratio wn = 6.28; Natural frequency, rad/s F = [ 1;-wn^2-2*z*wn]; G = [1-1; 2]; Hx = [1 ; 1]; sys = ss(f, G, Hx,); t = [:.1:1]; xo = [1;]; [y1,t1,x1] = initial(sys, xo, t); plot(t1,y1) grid on figure plot(y1(:,1),y1(:,2)) grid on 48
25 Effects of Damping Ratio on State-Plane Plots Damping ratio =.1 Damping ratio =.3 Damping ratio =.1 49 Response of Higher-Order LTI Systems Superposition of Sub-System Responses F System = F System 1 Effect of 1 on 2 Effect of 2 on 1 F System 2 Third-order system with uncoupled 1 st - and 2 nd - order sub-systems Coupling in first row and first column F = f 11 F = 1 f 32 f 33 f 11 f 12 f 13 f 21 1 f 31 f 32 f 33 x 1 x = x 2 x 3 5
26 Examples of Coupled and Uncoupled Third-Order Systems Third-order system with uncoupled 1 st - and 2 nd -order sub-systems F = x 1 x = x 2 x 3 Coupling in first row and first column Position Coupling, x 3 Rate Coupling,x 2 F = F = rd -Order LTI Systems with Coupled Response Initial Condition on x 1 Initial Condition on x 3 52
27 Next Time: Dynamic Effects of Feedback Control 53 Supplemental Material 54
28 Perturbed Initial Conditions Produce Perturbed Path Given Initial condition, control, and disturbance histories x(t ),u, w for t in [t o,t f ] Path (or trajectory) is approximated by executing a numerical algorithm Perturbing the initial condition produces a new path Both paths satisfy the same equations of motion x N = f[x N,u,w,t], x N (t ) given x N : Nominal path x: Perturbed path x = f[x,u,w,t], x(t ) given 55 Equilibrium Response of Second-Order LTI System System description F = f 11 f 12 ; G = g 1 ; L = l 1 f 21 f 22 g 2 l 2 State equilibrium depends on constant input values x 1 * x 2 * f 22 ( f 12 = ( ( f 21 f 11 ) ( f 11 f 22 ( f 12 f 21 ) * + g 1 g 2 / F = f 11 f 22 ( f 12 f 21, ) -. u *+ l 1 * + l 2, -. w * 56
Linearized Equations of Motion!
Linearized Equations of Motion Robert Stengel, Aircraft Flight Dynamics MAE 331, 216 Learning Objectives Develop linear equations to describe small perturbational motions Apply to aircraft dynamic equations
More informationTranslational and Rotational Dynamics!
Translational and Rotational Dynamics Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Copyright 217 by Robert Stengel. All rights reserved. For educational use only.
More informationTime Response of Linear, Time-Invariant (LTI) Systems Robert Stengel, Aircraft Flight Dynamics MAE 331, 2018
Time Response of Linear, Time-Invariant LTI Systems Robert Stengel, Aircraft Flight Dynamics MAE 331, 2018 Learning Objectives Methods of time-domain analysis Continuous- and discrete-time models Transient
More informationControl Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only.
Control Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 2017 Analog vs. digital systems Continuous- and Discretetime Dynamic Models Frequency Response Transfer Functions
More informationAlternative Expressions for the Velocity Vector Velocity restricted to the vertical plane. Longitudinal Equations of Motion
Linearized Longitudinal Equations of Motion Robert Stengel, Aircraft Flig Dynamics MAE 33, 008 Separate solutions for nominal and perturbation flig paths Assume that nominal path is steady and in the vertical
More informationOptimal Control and Estimation MAE 546, Princeton University Robert Stengel, Preliminaries!
Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, 2018 Copyright 2018 by Robert Stengel. All rights reserved. For educational use only. http://www.princeton.edu/~stengel/mae546.html
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationOptimal Control and Estimation MAE 546, Princeton University Robert Stengel, Preliminaries!
Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, 2017 Copyright 2017 by Robert Stengel. All rights reserved. For educational use only. http://www.princeton.edu/~stengel/mae546.html
More informationChapter 4 The Equations of Motion
Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations
More informationMinimization of Static! Cost Functions!
Minimization of Static Cost Functions Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2017 J = Static cost function with constant control parameter vector, u Conditions for
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationManipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA
Manipulator Dynamics 2 Forward Dynamics Problem Given: Joint torques and links geometry, mass, inertia, friction Compute: Angular acceleration of the links (solve differential equations) Solution Dynamic
More informationSTATE VARIABLE (SV) SYSTEMS
Copyright F.L. Lewis 999 All rights reserved Updated:Tuesday, August 05, 008 STATE VARIABLE (SV) SYSTEMS A natural description for dynamical systems is the nonlinear state-space or state variable (SV)
More informationLecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow
Lecture AC-1 Aircraft Dynamics Copy right 23 by Jon at h an H ow 1 Spring 23 16.61 AC 1 2 Aircraft Dynamics First note that it is possible to develop a very good approximation of a key motion of an aircraft
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationRobot Dynamics - Rotary Wing UAS: Control of a Quadrotor
Robot Dynamics Rotary Wing AS: Control of a Quadrotor 5-85- V Marco Hutter, Roland Siegwart and Thomas Stastny Robot Dynamics - Rotary Wing AS: Control of a Quadrotor 7..6 Contents Rotary Wing AS. Introduction
More informationTopic # Feedback Control
Topic #7 16.31 Feedback Control State-Space Systems What are state-space models? Why should we use them? How are they related to the transfer functions used in classical control design and how do we develop
More informationDynamics and Control of Rotorcraft
Dynamics and Control of Rotorcraft Helicopter Aerodynamics and Dynamics Abhishek Department of Aerospace Engineering Indian Institute of Technology, Kanpur February 3, 2018 Overview Flight Dynamics Model
More informationDynamic Optimal Control!
Dynamic Optimal Control! Robert Stengel! Robotics and Intelligent Systems MAE 345, Princeton University, 2017 Learning Objectives Examples of cost functions Necessary conditions for optimality Calculation
More informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationMinimization with Equality Constraints
Path Constraints and Numerical Optimization Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2015 State and Control Equality Constraints Pseudoinverse State and Control Inequality
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationFlight Dynamics and Control
Flight Dynamics and Control Lecture 1: Introduction G. Dimitriadis University of Liege Reference material Lecture Notes Flight Dynamics Principles, M.V. Cook, Arnold, 1997 Fundamentals of Airplane Flight
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationDynamics and control of mechanical systems
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day 2 (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
More informationVideo 2.1a Vijay Kumar and Ani Hsieh
Video 2.1a Vijay Kumar and Ani Hsieh Robo3x-1.3 1 Introduction to Lagrangian Mechanics Vijay Kumar and Ani Hsieh University of Pennsylvania Robo3x-1.3 2 Analytical Mechanics Aristotle Galileo Bernoulli
More informationFirst-Order Low-Pass Filter!
Filters, Cost Functions, and Controller Structures! Robert Stengel! Optimal Control and Estimation MAE 546! Princeton University, 217!! Dynamic systems as low-pass filters!! Frequency response of dynamic
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationVideo 3.1 Vijay Kumar and Ani Hsieh
Video 3.1 Vijay Kumar and Ani Hsieh Robo3x-1.3 1 Dynamics of Robot Arms Vijay Kumar and Ani Hsieh University of Pennsylvania Robo3x-1.3 2 Lagrange s Equation of Motion Lagrangian Kinetic Energy Potential
More informationIf the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.
Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate
More informationStochastic Optimal Control!
Stochastic Control! Robert Stengel! Robotics and Intelligent Systems, MAE 345, Princeton University, 2015 Learning Objectives Overview of the Linear-Quadratic-Gaussian (LQG) Regulator Introduction to Stochastic
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationNonlinear State Estimation! Extended Kalman Filters!
Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal
More informationSTRUCTURAL DYNAMICS BASICS:
BASICS: STRUCTURAL DYNAMICS Real-life structures are subjected to loads which vary with time Except self weight of the structure, all other loads vary with time In many cases, this variation of the load
More informationStochastic and Adaptive Optimal Control
Stochastic and Adaptive Optimal Control Robert Stengel Optimal Control and Estimation, MAE 546 Princeton University, 2018! Nonlinear systems with random inputs and perfect measurements! Stochastic neighboring-optimal
More informationMinimization with Equality Constraints!
Path Constraints and Numerical Optimization Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2017 State and Control Equality Constraints Pseudoinverse State and Control Inequality
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationNumerical Methods for Rigid Multibody Dynamics
Numerical Methods for Rigid Multibody Dynamics Claus Führer Centre for Mathematical Sciences Lund University Lappenranta 2012 Unit 0: Preface These notes serve as a skeleton for the compact course. They
More informationBoth Paths Satisfy the Dynamic Equations
Liearized Equatios ad Modes of Motio Robert Stegel, Aircraft Flight Dyamics MAE 331, 008 Liearizatio of oliear dyamic models Nomial flight path Perturbatios about the omial flight path Modes of motio Nomial
More informationFundamentals Physics. Chapter 15 Oscillations
Fundamentals Physics Tenth Edition Halliday Chapter 15 Oscillations 15-1 Simple Harmonic Motion (1 of 20) Learning Objectives 15.01 Distinguish simple harmonic motion from other types of periodic motion.
More informationState Space Representation
ME Homework #6 State Space Representation Last Updated September 6 6. From the homework problems on the following pages 5. 5. 5.6 5.7. 5.6 Chapter 5 Homework Problems 5.6. Simulation of Linear and Nonlinear
More informationCHAPTER 1. Introduction
CHAPTER 1 Introduction Linear geometric control theory was initiated in the beginning of the 1970 s, see for example, [1, 7]. A good summary of the subject is the book by Wonham [17]. The term geometric
More informationIntroduction to Optimization!
Introduction to Optimization! Robert Stengel! Robotics and Intelligent Systems, MAE 345, Princeton University, 2017 Optimization problems and criteria Cost functions Static optimality conditions Examples
More informationLAGRANGIAN AND HAMILTONIAN
LAGRANGIAN AND HAMILTONIAN A. Constraints and Degrees of Freedom. A constraint is a restriction on the freedom of motion of a system of particles in the form of a condition. The number of independent ways
More informationFree Vibration of Single-Degree-of-Freedom (SDOF) Systems
Free Vibration of Single-Degree-of-Freedom (SDOF) Systems Procedure in solving structural dynamics problems 1. Abstraction/modeling Idealize the actual structure to a simplified version, depending on the
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationTranslational Mechanical Systems
Translational Mechanical Systems Basic (Idealized) Modeling Elements Interconnection Relationships -Physical Laws Derive Equation of Motion (EOM) - SDOF Energy Transfer Series and Parallel Connections
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationAircraft Flight Dynamics Robert Stengel MAE 331, Princeton University, 2018
Aircraft Flight Dynamics Robert Stengel MAE 331, Princeton University, 2018 Course Overview Introduction to Flight Dynamics Math Preliminaries Copyright 2018 by Robert Stengel. All rights reserved. For
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 informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationSystem simulation using Matlab, state plane plots
System simulation using Matlab, state plane plots his lab is mainly concerned with making state plane (also referred to as phase plane ) plots for various linear and nonlinear systems with two states he
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More information6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant
DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional
More informationFundamental principles
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
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 informationAnalytical Dynamics: Lagrange s Equation and its Application A Brief Introduction
Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction D. S. Stutts, Ph.D. Associate Professor of Mechanical Engineering Missouri University of Science and Technology Rolla,
More informationENGG 5402 Course Project: Simulation of PUMA 560 Manipulator
ENGG 542 Course Project: Simulation of PUMA 56 Manipulator ZHENG Fan, 115551778 mrzhengfan@gmail.com April 5, 215. Preface This project is to derive programs for simulation of inverse dynamics and control
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationRobotics. Dynamics. University of Stuttgart Winter 2018/19
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler, joint space control, reference trajectory following, optimal operational
More informationSimple Harmonic Motion
Physics 7B-1 (A/B) Professor Cebra Winter 010 Lecture 10 Simple Harmonic Motion Slide 1 of 0 Announcements Final exam will be next Wednesday 3:30-5:30 A Formula sheet will be provided Closed-notes & closed-books
More informationAircraft Flight Dynamics!
Aircraft Flight Dynamics Robert Stengel MAE 331, Princeton University, 2016 Course Overview Introduction to Flight Dynamics Math Preliminaries Copyright 2016 by Robert Stengel. All rights reserved. For
More informationPeriodic Motion. Periodic motion is motion of an object that. regularly repeats
Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A special kind of periodic motion occurs in mechanical systems
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More informationCIS 4930/6930: Principles of Cyber-Physical Systems
CIS 4930/6930: Principles of Cyber-Physical Systems Chapter 2: Continuous Dynamics Hao Zheng Department of Computer Science and Engineering University of South Florida H. Zheng (CSE USF) CIS 4930/6930:
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More information41514 Dynamics of Machinery
41514 Dynamics of Machinery Theory, Experiment, Phenomenology and Industrial Applications Ilmar Ferreira Santos 1. Recapitulation Mathematical Modeling & Steps 2. Example System of Particle 3. Example
More informationIn this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:
Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this
More informationSimplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System
NASA Technical Memorandum 112868 Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System Nelson J. Groom Langley Research Center, Hampton, Virginia June 1997 National
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More informationIntroduction to Vibration. Professor Mike Brennan
Introduction to Vibration Professor Mie Brennan Introduction to Vibration Nature of vibration of mechanical systems Free and forced vibrations Frequency response functions Fundamentals For free vibration
More information1.053J/2.003J Dynamics and Control I Fall Final Exam 18 th December, 2007
1.053J/2.003J Dynamics and Control I Fall 2007 Final Exam 18 th December, 2007 Important Notes: 1. You are allowed to use three letter-size sheets (two-sides each) of notes. 2. There are five (5) problems
More information4.1 Important Notes on Notation
Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationEE Homework 3 Due Date: 03 / 30 / Spring 2015
EE 476 - Homework 3 Due Date: 03 / 30 / 2015 Spring 2015 Exercise 1 (10 points). Consider the problem of two pulleys and a mass discussed in class. We solved a version of the problem where the mass was
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationSCHEME OF BE 100 ENGINEERING MECHANICS DEC 2015
Part A Qn. No SCHEME OF BE 100 ENGINEERING MECHANICS DEC 201 Module No BE100 ENGINEERING MECHANICS Answer ALL Questions 1 1 Theorem of three forces states that three non-parallel forces can be in equilibrium
More information2.003 Quiz #1 Review
2.003J Spring 2011: Dynamics and Control I Quiz #1 Review Massachusetts Institute of Technology March 5th, 2011 Department of Mechanical Engineering March 6th, 2011 1 Reference Frames 2.003 Quiz #1 Review
More informationOrbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel
Orbital Mechanics Space System Design, MAE 342, Princeton University Robert Stengel Conic section orbits Equations of motion Momentum and energy Kepler s Equation Position and velocity in orbit Copyright
More informationNonholonomic Constraints Examples
Nonholonomic Constraints Examples Basilio Bona DAUIN Politecnico di Torino July 2009 B. Bona (DAUIN) Examples July 2009 1 / 34 Example 1 Given q T = [ x y ] T check that the constraint φ(q) = (2x + siny
More informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Euler-Lagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered Euler-Lagrange formulation manipulator inertia matrix
More informationDYNAMICS OF SERIAL ROBOTIC MANIPULATORS
DYNAMICS OF SERIAL ROBOTIC MANIPULATORS NOMENCLATURE AND BASIC DEFINITION We consider here a mechanical system composed of r rigid bodies and denote: M i 6x6 inertia dyads of the ith body. Wi 6 x 6 angular-velocity
More informationAircraft Equations of Motion: Translation and Rotation Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018
Aircraft Equations of Motion: Translation and Rotation Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018 Learning Objectives What use are the equations of motion? How is the angular orientation of
More informationFirst-Order Low-Pass Filter
Filters, Cost Functions, and Controller Structures Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 218! Dynamic systems as low-pass filters! Frequency response of dynamic systems!
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationAB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT
FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be
More informationIntroduction to Haptic Systems
Introduction to Haptic Systems Félix Monasterio-Huelin & Álvaro Gutiérrez & Blanca Larraga October 8, 2018 Contents Contents 1 List of Figures 1 1 Introduction 2 2 DC Motor 3 3 1 DOF DC motor model with
More informationAdaptive Trim and Trajectory Following for a Tilt-Rotor Tricopter Ahmad Ansari, Anna Prach, and Dennis S. Bernstein
7 American Control Conference Sheraton Seattle Hotel May 4 6, 7, Seattle, USA Adaptive Trim and Trajectory Following for a Tilt-Rotor Tricopter Ahmad Ansari, Anna Prach, and Dennis S. Bernstein Abstract
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Rigid body physics Particle system Most simple instance of a physics system Each object (body) is a particle Each particle
More informationInverse differential kinematics Statics and force transformations
Robotics 1 Inverse differential kinematics Statics and force transformations Prof Alessandro De Luca Robotics 1 1 Inversion of differential kinematics! find the joint velocity vector that realizes a desired
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
More informationAdvanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich
Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion
More informationChapter 6: Momentum Analysis of Flow Systems
Chapter 6: Momentum Analysis of Flow Systems Introduction Fluid flow problems can be analyzed using one of three basic approaches: differential, experimental, and integral (or control volume). In Chap.
More informationAutomatic Control Systems. -Lecture Note 15-
-Lecture Note 15- Modeling of Physical Systems 5 1/52 AC Motors AC Motors Classification i) Induction Motor (Asynchronous Motor) ii) Synchronous Motor 2/52 Advantages of AC Motors i) Cost-effective ii)
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More information