Spacecraft and Aircraft Dynamics

Size: px
Start display at page:

Download "Spacecraft and Aircraft Dynamics"

Transcription

1 Spacecraft and Aircraft Dynamics Matthew M. Peet Illinois Institute of Technology Lecture 11: Longitudinal Dynamics

2 Aircraft Dynamics Lecture 11 In this Lecture we will cover: Longitudinal Dynamics: Finding dimensional coefficients from non-dimensional coefficients Eigenvalue Analysis Approximate modal behavior short period mode phugoid mode M. Peet Lecture 11: 2 / 25

3 Review: Longitudinal Dynamics Combined Terms u+ θgcosθ 0 = X u u+x w w +X δe δ e +X δt δ T ẇ + θgsinθ 0 u 0 θ = Z u u+z w w +Zẇ ẇ +Z q θ +Z δe δ e +Z δt δ T θ = M u u+m w w +Mẇ ẇ +M q θ +M δe δ e +M δt δ T M. Peet Lecture 11: 3 / 25

4 Force Coefficients Force/Moment Coefficients can be found in Table 3.5 of Nelson M. Peet Lecture 11: 4 / 25

5 Nondimensional Force Coefficients Nondimensional Force/Moment Coefficients can be found in Table 3.3 of Nelson M. Peet Lecture 11: 5 / 25

6 State-Space From the homework, we have a state-space representation of form u ẇ θ = [ A ] u w θ +[ B ][ ] δ e δ T q q Where we get A and B from X u, Z u, etc. Recall: Eigenvalues of A define stability of ẋ = Ax. A is 4 4, so A has 4 eigenvalues. Stable if eigenvalues all have negative real part. M. Peet Lecture 11: 6 / 25

7 Natural Motion We mentioned that A is Stable if eigenvalues all have negative real part. Now we say more: Eigenvalues have the form λ = λ R ±λ I ı If we have a pair of complex eigenvalues, then we have two more concepts: 1. Natural Frequency: ω n = λ 2 R +λ2 I 2. Damping Ratio: d = λ R ω n M. Peet Lecture 11: 7 / 25

8 Natural Frequency Natural frequency is how fast the the motion oscillates. Closely related is the Definition 1. The Period is the time take to complete one oscillation τ = 2π ω n M. Peet Lecture 11: 8 / 25

9 Damping Ratio Damping ratio is how much amplitude decays per oscillation. Even if d is large, may decay slowly is ω n is small Closely related is Definition 2. The Half-Life is the time taken for the amplitude to decay by half. γ =.693 λ R M. Peet Lecture 11: 9 / 25

10 State-Space Example: Uncontrolled Motion C172: V 0 = 132kt, 5000ft. u u ẇ q = w q θ θ M. Peet Lecture 11: 10 / 25

11 State-Space Example: Uncontrolled Motion Using the Matlab command [u,v] = eigs(a), we find the eigenvalues as Phugoid (Long-Period) Mode and Eigenvectors Short-Period Mode and Eigenvectors λ 1,2 =.0209±.18ı ν 1,2 = ± ı λ 3,4 = 4.13±4.39ı 1 0 ν 3,4 = ± ı Notice that this is hard to interpret. Lets scale u and q by equilibrium values. M. Peet Lecture 11: 11 / 25

12 State-Space Example: Uncontrolled Motion After scaling the state by the equilibrium values, we find the eigenvalues unchanged (Why?) as Phugoid (Long-Period) Mode λ 1,2 =.0209±.18ı but clearer Eigenvectors ν 1,2 = ± ı Natural Frequency: ω n =.181rad/s Damping Ratio: d =.115 Period: τ = 34.7s Half-Life: γ = 33.16s Motion dominated by variables u and θ. M. Peet Lecture 11: 12 / 25

13 Modal Illustration M. Peet Lecture 11: 13 / 25

14 State-Space Example: Long Period Mode M. Peet Lecture 11: 14 / 25

15 State-Space Example: Uncontrolled Motion Using the Matlab command [u,v] = eigs(a), we find the eigenvalues as Short-Period Mode λ 3,4 = 4.13±4.39ı and Eigenvectors ν 3,4 = ± ı Natural Frequency: ω n = 6.03rad/s Damping Ratio: d =.685 Period: τ = 1.04s Half-Life: γ =.167s Motion dominated by variables w and q. M. Peet Lecture 11: 15 / 25

16 Modal Illustration M. Peet Lecture 11: 16 / 25

17 State-Space Example: Short Period Mode M. Peet Lecture 11: 17 / 25

18 State-Space Modal Approximations Now that we know that longitudinal dynamics have two modes: Short Period Mode Phugoid Mode (Long-Period Mode) Short Period Mode: fix u = 0 and w = 0. study variation in θ and q. Similar to Static Longitudinal Stability Long Period Mode: fix q = 0 and w = 0. study variation in θ and u. Now we develop some simplified expressions to study these modes. M. Peet Lecture 11: 18 / 25

19 Short Period Approximation For the short period mode, we have the following dynamics: [ ẇ ] [ ] Z α [ u = u 0 1 w ] [ ] Z δe 0 q M α + M u0 + u 0 u αzα 0 M q +M α q M δe + M α+z δe u 0 [ ] θ = A sp +B q sp δ e δ e To understand stability, we need the eigenvalues of A sp. Eigenvalues are solutions of det(λi A sp ) = 0. Thus we want to solve det(λi A sp ) = λ 2 (M q + Z α u 0 +M α )λ+( Z αm q u 0 M α ) = 0 We use the quadratic formula (Lecture 1): λ 3,4 = 1 2 (M q +M α + Z α )± 1 (M q +M α + Z α ) u 0 2 u 2 Z α 4(m q M α ) 0 u 0 M. Peet Lecture 11: 19 / 25

20 Short Period Approximation Frequency and Damping Ratio λ 3,4 = 1 2 (M q +M α + Z α )± 1 (M q +M α + Z α ) u 0 2 u 2 Z α 4(m q M α ) 0 u 0 This leads to the Approximation Equations: Natural Frequency: Damping Ratio: d sp = 1 2 ω sp = M q Z α u 0 M q +M α + Zα u 0 ω sp M. Peet Lecture 11: 20 / 25

21 Short Period Mode Approximations Example: C127 Approximate Natural Frequency: ω sp = = 6.10rad/s 219 True Natural Frequency: ω sp = 6.03rad/s Approximate Damping Ratio: True Damping Ratio: So, generally good agreement. d sp = d sp =.685rad/s =.683rad/s M. Peet Lecture 11: 21 / 25

22 Long Period Approximation Long period motion considers only motion in u and θ. [ u θ] [ ][ Xu g u = Zu u 0 0 θ] This time, we must solve the simple expression Using the quadratic formula, we get det(λi A) = λ 2 X u λ Z u u 0 g = 0 X u ± λ 1,2 = X 2 u +4Zu u 0 g 2 M. Peet Lecture 11: 22 / 25

23 Long Period Approximation X u ± Xu 2 +4Zu u 0 g λ 1,2 = 2 This leads to the Approximation Equations: Natural Frequency: Damping Ratio: ω lp = Z ug u 0 d lp = X u 2ω lp M. Peet Lecture 11: 23 / 25

24 Long Period Mode Approximations Example: C127 Approximate Natural Frequency: ω lp =.181rad/s True Natural Frequency: ω lp =.208rad/s Approximate Damping Ratio: d lp =.115rad/s True Damping Ratio: d lp =.106rad/s So, generally good, but not as good. Why? M. Peet Lecture 11: 24 / 25

25 Conclusion In this lecture, we covered: How to find and interpret the eigenvalues and eigenvectors of a state-space matrix Natural Frequency Damping Ratio How to identify Long Period Eigenvlaues/Motion Short Period Eigenvalues/Motion Modal Approximations Phugoid and Short-Period Modes Formulas for natural frequency Formulae for damping ratio M. Peet Lecture 11: 25 / 25

16.333: Lecture # 7. Approximate Longitudinal Dynamics Models

16.333: Lecture # 7. Approximate Longitudinal Dynamics Models 16.333: Lecture # 7 Approxiate Longitudinal Dynaics Models A couple ore stability derivatives Given ode shapes found identify sipler odels that capture the ain responses Fall 24 16.333 6 1 More Stability

More information

16.333: Lecture # 14. Equations of Motion in a Nonuniform Atmosphere. Gusts and Winds

16.333: Lecture # 14. Equations of Motion in a Nonuniform Atmosphere. Gusts and Winds 16.333: Lecture # 14 Equations of Motion in a Nonuniform Atmosphere Gusts and Winds 1 Fall 2004 16.333 12 2 Equations of Motion Analysis to date has assumed that the atmosphere is calm and fixed Rarely

More information

Modal Decomposition and the Time-Domain Response of Linear Systems 1

Modal Decomposition and the Time-Domain Response of Linear Systems 1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout

More information

GENERALIZED STABILITY OF THE TWO-LAYER MODEL

GENERALIZED STABILITY OF THE TWO-LAYER MODEL GENERALIZED STABILITY OF THE TWO-LAYER MODEL The simplest mid-latitude jet model supporting the baroclinic growth mechanism is the two-layer model The equations for the barotropic and baroclinic geostrophic

More information

Chapter 8: Frequency Domain Analysis

Chapter 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 information

ME 563 HOMEWORK # 5 SOLUTIONS Fall 2010

ME 563 HOMEWORK # 5 SOLUTIONS Fall 2010 ME 563 HOMEWORK # 5 SOLUTIONS Fall 2010 PROBLEM 1: You are given the lumped parameter dynamic differential equations of motion for a two degree-offreedom model of an automobile suspension system for small

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real

More information

Chapter 8 Dynamic stability analysis II Longitudinal motion (Lectures 28 to 32)

Chapter 8 Dynamic stability analysis II Longitudinal motion (Lectures 28 to 32) Chapter 8 Dynamic stability analysis II Longitudinal motion (Lectures 28 to 32) Keywords : Stability quartic or characteristic equation for longitudinal motion and its solution ; roots of characteristic

More information

MAE 142 Homework #5 Due Friday, March 13, 2009

MAE 142 Homework #5 Due Friday, March 13, 2009 MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are

More information

Applications Linear Control Design Techniques in Aircraft Control I

Applications Linear Control Design Techniques in Aircraft Control I Lecture 29 Applications Linear Control Design Techniques in Aircraft Control I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics Brief Review

More information

3. Mathematical Properties of MDOF Systems

3. Mathematical Properties of MDOF Systems 3. Mathematical Properties of MDOF Systems 3.1 The Generalized Eigenvalue Problem Recall that the natural frequencies ω and modes a are found from [ - ω 2 M + K ] a = 0 or K a = ω 2 M a Where M and K are

More information

Lecture #AC 3. Aircraft Lateral Dynamics. Spiral, Roll, and Dutch Roll Modes

Lecture #AC 3. Aircraft Lateral Dynamics. Spiral, Roll, and Dutch Roll Modes Lecture #AC 3 Aircraft Lateral Dynamics Spiral, Roll, and Dutch Roll Modes Copy right 2003 by Jon at h an H ow 1 Spring 2003 16.61 AC 3 2 Aircraft Lateral Dynamics Using a procedure similar to the longitudinal

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 2: Drawing Bode Plots, Part 2 Overview In this Lecture, you will learn: Simple Plots Real Zeros Real Poles Complex

More information

11/17/10. Chapter 14. Oscillations. Chapter 14. Oscillations Topics: Simple Harmonic Motion. Simple Harmonic Motion

11/17/10. Chapter 14. Oscillations. Chapter 14. Oscillations Topics: Simple Harmonic Motion. Simple Harmonic Motion 11/17/10 Chapter 14. Oscillations This striking computergenerated image demonstrates an important type of motion: oscillatory motion. Examples of oscillatory motion include a car bouncing up and down,

More information

MATH 251 Examination II July 28, Name: Student Number: Section:

MATH 251 Examination II July 28, Name: Student Number: Section: MATH 251 Examination II July 28, 2008 Name: Student Number: Section: This exam has 9 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must be shown.

More information

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution

More information

Autonomous system = system without inputs

Autonomous 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 information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Identification 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 information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real Poles

More information

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky

ELEC 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 information

DESIGN PROJECT REPORT: Longitudinal and lateral-directional stability augmentation of Boeing 747 for cruise flight condition.

DESIGN PROJECT REPORT: Longitudinal and lateral-directional stability augmentation of Boeing 747 for cruise flight condition. DESIGN PROJECT REPORT: Longitudinal and lateral-directional stability augmentation of Boeing 747 for cruise flight condition. Prepared By: Kushal Shah Advisor: Professor John Hodgkinson Graduate Advisor:

More information

MAE 143B - Homework 8 Solutions

MAE 143B - Homework 8 Solutions MAE 43B - Homework 8 Solutions P6.4 b) With this system, the root locus simply starts at the pole and ends at the zero. Sketches by hand and matlab are in Figure. In matlab, use zpk to build the system

More information

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,

More information

April 15, 2011 Sample Quiz and Exam Questions D. A. Caughey Page 1 of 9

April 15, 2011 Sample Quiz and Exam Questions D. A. Caughey Page 1 of 9 April 15, 2011 Sample Quiz Exam Questions D. A. Caughey Page 1 of 9 These pages include virtually all Quiz, Midterm, Final Examination questions I have used in M&AE 5070 over the years. Note that some

More information

Chemistry 431. Lecture 23

Chemistry 431. Lecture 23 Chemistry 431 Lecture 23 Introduction The Larmor Frequency The Bloch Equations Measuring T 1 : Inversion Recovery Measuring T 2 : the Spin Echo NC State University NMR spectroscopy The Nuclear Magnetic

More information

Linear Algebra and ODEs review

Linear Algebra and ODEs review Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A

More information

Laboratory handout 5 Mode shapes and resonance

Laboratory handout 5 Mode shapes and resonance laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful

More information

EE221A Linear System Theory Final Exam

EE221A Linear System Theory Final Exam EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,

More information

Time 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 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 information

Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli

Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)

More information

Modeling II Linear Stability Analysis and Wave Equa9ons

Modeling II Linear Stability Analysis and Wave Equa9ons Modeling II Linear Stability Analysis and Wave Equa9ons Nondimensional Equa9ons From previous lecture, we have a system of nondimensional PDEs: (21.1) (21.2) (21.3) where here the * sign has been dropped

More information

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 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 information

( ) Chapter 3: Free Vibration of the Breakwater. 3.1 Introduction

( ) Chapter 3: Free Vibration of the Breakwater. 3.1 Introduction Chapter : Free Vibration of the Breakwater. Introduction Before conducting forcing analyses on the breakwater, a free vibration study is necessary. This chapter will describe the analysis conducted to

More information

Math 1302, Week 8: Oscillations

Math 1302, Week 8: Oscillations Math 302, Week 8: Oscillations T y eq Y y = y eq + Y mg Figure : Simple harmonic motion. At equilibrium the string is of total length y eq. During the motion we let Y be the extension beyond equilibrium,

More information

MATH 23 Exam 2 Review Solutions

MATH 23 Exam 2 Review Solutions MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution

More information

21.55 Worksheet 7 - preparation problems - question 1:

21.55 Worksheet 7 - preparation problems - question 1: Dynamics 76. Worksheet 7 - preparation problems - question : A coupled oscillator with two masses m and positions x (t) and x (t) is described by the following equations of motion: ẍ x + 8x ẍ x +x A. Write

More information

AA 242B / ME 242B: Mechanical Vibrations (Spring 2016)

AA 242B / ME 242B: Mechanical Vibrations (Spring 2016) AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally

More information

Physics 105, Spring 2015, Reinsch. Homework Assignment 1 (This is a five page document. More hints added Feb. 2, 9pm. Problem statements unchanged.

Physics 105, Spring 2015, Reinsch. Homework Assignment 1 (This is a five page document. More hints added Feb. 2, 9pm. Problem statements unchanged. Physics 105, Spring 2015, Reinsch Homework Assignment 1 (This is a five page document. More hints added Feb. 2, 9pm. Problem statements unchanged.) Due Wednesday, February 4, 5:00 pm Problem 1 A particle

More information

Waves & Oscillations

Waves & Oscillations Physics 42200 Waves & Oscillations Lecture 21 Review Spring 2013 Semester Matthew Jones Midterm Exam: Date: Wednesday, March 6 th Time: 8:00 10:00 pm Room: PHYS 203 Material: French, chapters 1-8 Review

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 12: Overview In this Lecture, you will learn: Review of Feedback Closing the Loop Pole Locations Changing the Gain

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

Problem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:

Problem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations: EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting

More information

Math 308 Exam II Practice Problems

Math 308 Exam II Practice Problems Math 38 Exam II Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..

More information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t

More information

MECH 6091 Flight Control Systems Final Course Project

MECH 6091 Flight Control Systems Final Course Project MECH 6091 Flight Control Systems Final Course Project F-16 Autopilot Design Lizeth Buendia Rodrigo Lezama Daniel Delgado December 16, 2011 1 AGENDA Theoretical Background F-16 Model and Linearization Controller

More information

CDS 101 Precourse Phase Plane Analysis and Stability

CDS 101 Precourse Phase Plane Analysis and Stability CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/

More information

EE364 Review Session 4

EE364 Review Session 4 EE364 Review Session 4 EE364 Review Outline: convex optimization examples solving quasiconvex problems by bisection exercise 4.47 1 Convex optimization problems we have seen linear programming (LP) quadratic

More information

MATH 53H - Solutions to Problem Set III

MATH 53H - Solutions to Problem Set III MATH 53H - Solutions to Problem Set III. Fix λ not an eigenvalue of L. Then det(λi L) 0 λi L is invertible. We have then p L+N (λ) = det(λi L N) = det(λi L) det(i (λi L) N) = = p L (λ) det(i (λi L) N)

More information

kx m x B N 1 C L, M Mg θ

kx m x B N 1 C L, M Mg θ .004 MODELING DYNAMICS AND CONTROL II Spring 00 Solutions to Problem Set No. 7 Problem 1. Pendulum mounted on elastic support. This problem is an execise in the application of momentum principles. Two

More information

18. Linearization: the phugoid equation as example

18. Linearization: the phugoid equation as example 79 18. Linearization: the phugoid equation as example Linearization is one of the most important and widely used mathematical terms in applications to Science and Engineering. In the context of Differential

More information

Introduction to Modern Control MT 2016

Introduction 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 information

Aircraft Flight Dynamics & Vortex Lattice Codes

Aircraft Flight Dynamics & Vortex Lattice Codes Aircraft Flight Dynamics Vortex Lattice Codes AA241X April 14 2014 Stanford University Overview 1. Equations of motion 2. Non-dimensional EOM Aerodynamics 3. Trim Analysis Longitudinal Lateral 4. Linearized

More information

MATH 24 EXAM 3 SOLUTIONS

MATH 24 EXAM 3 SOLUTIONS MATH 4 EXAM 3 S Consider the equation y + ω y = cosω t (a) Find the general solution of the homogeneous equation (b) Find the particular solution of the non-homogeneous equation using the method of Undetermined

More information

22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes

22.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 information

Aeroelasticity. Lecture 4: Flight Flutter Testing. G. Dimitriadis. AERO0032-1, Aeroelasticity and Experimental Aerodynamics, Lecture 4

Aeroelasticity. Lecture 4: Flight Flutter Testing. G. Dimitriadis. AERO0032-1, Aeroelasticity and Experimental Aerodynamics, Lecture 4 Aeroelasticity Lecture 4: Flight Flutter Testing G. Dimitriadis AERO0032-1, Aeroelasticity and Experimental Aerodynamics, Lecture 4 1 Main reference A Historical Overview of Flight Flutter Testing, Michael

More information

Practice Problems For Test 3

Practice Problems For Test 3 Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)

More information

THE subject of the analysis is system composed by

THE subject of the analysis is system composed by MECHANICAL VIBRATION ASSIGNEMENT 1 On 3 DOF system identification Diego Zenari, 182160, M.Sc Mechatronics engineering Abstract The present investigation carries out several analyses on a 3-DOF system.

More information

Flying Qualities Criteria Robert Stengel, Aircraft Flight Dynamics MAE 331, 2018

Flying Qualities Criteria Robert Stengel, Aircraft Flight Dynamics MAE 331, 2018 Flying Qualities Criteria Robert Stengel, Aircraft Flight Dynamics MAE 331, 2018 Learning Objectives MIL-F-8785C criteria CAP, C*, and other longitudinal criteria ϕ/β, ω ϕ /ω d, and other lateral-directional

More information

Practice Problems For Test 3

Practice Problems For Test 3 Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)

More information

ME 563 HOMEWORK # 7 SOLUTIONS Fall 2010

ME 563 HOMEWORK # 7 SOLUTIONS Fall 2010 ME 563 HOMEWORK # 7 SOLUTIONS Fall 2010 PROBLEM 1: Given the mass matrix and two undamped natural frequencies for a general two degree-of-freedom system with a symmetric stiffness matrix, find the stiffness

More information

Design of Structures for Earthquake Resistance

Design of Structures for Earthquake Resistance NATIONAL TECHNICAL UNIVERSITY OF ATHENS Design of Structures for Earthquake Resistance Basic principles Ioannis N. Psycharis Lecture 3 MDOF systems Equation of motion M u + C u + K u = M r x g(t) where:

More information

1 Simple Harmonic Oscillator

1 Simple Harmonic Oscillator Physics 1a Waves Lecture 3 Caltech, 10/09/18 1 Simple Harmonic Oscillator 1.4 General properties of Simple Harmonic Oscillator 1.4.4 Superposition of two independent SHO Suppose we have two SHOs described

More information

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions

More information

Lecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow

Lecture 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 information

Math 504 (Fall 2011) 1. (*) Consider the matrices

Math 504 (Fall 2011) 1. (*) Consider the matrices Math 504 (Fall 2011) Instructor: Emre Mengi Study Guide for Weeks 11-14 This homework concerns the following topics. Basic definitions and facts about eigenvalues and eigenvectors (Trefethen&Bau, Lecture

More information

The (Fast) Fourier Transform

The (Fast) Fourier Transform The (Fast) Fourier Transform The Fourier transform (FT) is the analog, for non-periodic functions, of the Fourier series for periodic functions can be considered as a Fourier series in the limit that the

More information

University of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.

University of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA. University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is

More information

ECE 240a - Notes on Spontaneous Emission within a Cavity

ECE 240a - Notes on Spontaneous Emission within a Cavity ECE 0a - Notes on Spontaneous Emission within a Cavity Introduction Many treatments of lasers treat the rate of spontaneous emission as specified by the time constant τ sp as a constant that is independent

More information

I ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l

I ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l Advanced Control Theory Homework #3 Student No.597 Name : Jinseong Kim. The simple pendulum dynamics is given by g sin L A. Derive above dynamic equation from the free body diagram. B. Find the equilibrium

More information

Vibrations Qualifying Exam Study Material

Vibrations Qualifying Exam Study Material Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors

More information

Dynamics and Control Preliminary Examination Topics

Dynamics and Control Preliminary Examination Topics Dynamics and Control Preliminary Examination Topics 1. Particle and Rigid Body Dynamics Meirovitch, Leonard; Methods of Analytical Dynamics, McGraw-Hill, Inc New York, NY, 1970 Chapters 1-5 2. Atmospheric

More information

Topic # /31 Feedback Control Systems

Topic # /31 Feedback Control Systems Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31

More information

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and

More information

Solution of Matrix Eigenvalue Problem

Solution of Matrix Eigenvalue Problem Outlines October 12, 2004 Outlines Part I: Review of Previous Lecture Part II: Review of Previous Lecture Outlines Part I: Review of Previous Lecture Part II: Standard Matrix Eigenvalue Problem Other Forms

More information

PHY217: Vibrations and Waves

PHY217: Vibrations and Waves Assessed Problem set 1 Issued: 5 November 01 PHY17: Vibrations and Waves Deadline for submission: 5 pm Thursday 15th November, to the V&W pigeon hole in the Physics reception on the 1st floor of the GO

More information

San José State University Aerospace Engineering Department AE168: Aerospace Vehicle Dynamics & Control Fall 2016

San José State University Aerospace Engineering Department AE168: Aerospace Vehicle Dynamics & Control Fall 2016 San José State University Aerospace Engineering Department AE168: Aerospace Vehicle Dynamics & Control Fall 216 Course and Contact Information Instructor: Office Location: Email: Office Hours: Class Days/Time:

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect

More information

CDS 101/110: Lecture 3.1 Linear Systems

CDS 101/110: Lecture 3.1 Linear Systems CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,

More information

Properties of Linear Transformations from R n to R m

Properties of Linear Transformations from R n to R m Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation

More information

Physics Lecture 13. P. Gutierrez. Department of Physics & Astronomy University of Oklahoma

Physics Lecture 13. P. Gutierrez. Department of Physics & Astronomy University of Oklahoma Physics 2514 Lecture 13 P. Gutierrez Department of Physics & Astronomy University of Oklahoma P. Gutierrez (University of Oklahoma) Physics 2514 February 23, 2011 1 / 14 Goal Goals for today s lecture:

More information

1-DOF Forced Harmonic Vibration. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 8 Fall 2011

1-DOF Forced Harmonic Vibration. MCE371: Vibrations. Prof. Richter. Department of Mechanical Engineering. Handout 8 Fall 2011 MCE371: Vibrations Prof. Richter Department of Mechanical Engineering Handout 8 Fall 2011 Harmonic Forcing Functions Transient vs. Steady Vibration Follow Palm, Sect. 4.1, 4.9 and 4.10 Harmonic forcing

More information

RAD229: Midterm Exam 2015/2016 October 19, Minutes. Please do not proceed to the next page until the exam begins.

RAD229: Midterm Exam 2015/2016 October 19, Minutes. Please do not proceed to the next page until the exam begins. RAD229: Midterm Exam 2015/2016 October 19, 2015 ---- 75 Minutes Name: Student ID: General Instructions: 1. Write your name legibly on this page. 2. You may use notes including lectures, homework, solutions

More information

CDS 101/110: Lecture 3.1 Linear Systems

CDS 101/110: Lecture 3.1 Linear Systems CDS /: Lecture 3. Linear Systems Goals for Today: Revist and motivate linear time-invariant system models: Summarize properties, examples, and tools Convolution equation describing solution in response

More information

Systems Analysis and Control

Systems Analysis and Control Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using

More information

TTK4150 Nonlinear Control Systems Solution 6 Part 2

TTK4150 Nonlinear Control Systems Solution 6 Part 2 TTK4150 Nonlinear Control Systems Solution 6 Part 2 Department of Engineering Cybernetics Norwegian University of Science and Technology Fall 2003 Solution 1 Thesystemisgivenby φ = R (φ) ω and J 1 ω 1

More information

Multiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with.

Multiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with. Multiple Degree of Freedom Systems The Millennium bridge required many degrees of freedom to model and design with. The first step in analyzing multiple degrees of freedom (DOF) is to look at DOF DOF:

More information

III.H Zeroth Order Hydrodynamics

III.H Zeroth Order Hydrodynamics III.H Zeroth Order Hydrodynamics As a first approximation, we shall assume that in local equilibrium, the density f 1 at each point in space can be represented as in eq.(iii.56), i.e. [ ] p m q, t)) f

More information

Appendix A Equations of Motion in the Configuration and State Spaces

Appendix A Equations of Motion in the Configuration and State Spaces Appendix A Equations of Motion in the Configuration and State Spaces A.1 Discrete Linear Systems A.1.1 Configuration Space Consider a system with a single degree of freedom and assume that the equation

More information

A Level. A Level Physics. Oscillations (Answers) AQA, Edexcel. Name: Total Marks: /30

A Level. A Level Physics. Oscillations (Answers) AQA, Edexcel. Name: Total Marks: /30 Visit http://www.mathsmadeeasy.co.uk/ for more fantastic resources. AQA, Edexcel A Level A Level Physics Oscillations (Answers) Name: Total Marks: /30 Maths Made Easy Complete Tuition Ltd 2017 1. The graph

More information

Spacecraft and Aircraft Dynamics

Spacecraft and Aircraft Dynamics Spacecraft and Aircraft Dynamics Matthew M. Peet Illinois Institute of Technology Lecture 4: Contributions to Longitudinal Stability Aircraft Dynamics Lecture 4 In this lecture, we will discuss Airfoils:

More information

PHYSICS. Chapter 15 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.

PHYSICS. Chapter 15 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc. PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 15 Lecture RANDALL D. KNIGHT Chapter 15 Oscillations IN THIS CHAPTER, you will learn about systems that oscillate in simple harmonic

More information

Solution Set 3. Hand out : i d dt. Ψ(t) = Ĥ Ψ(t) + and

Solution Set 3. Hand out : i d dt. Ψ(t) = Ĥ Ψ(t) + and Physikalische Chemie IV Magnetische Resonanz HS Solution Set 3 Hand out : 5.. Repetition. The Schrödinger equation describes the time evolution of a closed quantum system: i d dt Ψt Ĥ Ψt Here the state

More information

Advanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One

Advanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to

More information

MIT (Spring 2014)

MIT (Spring 2014) 18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your

More information

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017

More information

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems R. M. Murray Fall 2004 CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 101/110 Homework Set #2 Issued: 4 Oct 04 Due: 11 Oct 04 Note: In the upper left hand corner of the first page

More information

Differential Equations and Linear Algebra Exercises. Department of Mathematics, Heriot-Watt University, Edinburgh EH14 4AS

Differential Equations and Linear Algebra Exercises. Department of Mathematics, Heriot-Watt University, Edinburgh EH14 4AS Differential Equations and Linear Algebra Exercises Department of Mathematics, Heriot-Watt University, Edinburgh EH14 4AS CHAPTER 1 Linear second order ODEs Exercises 1.1. (*) 1 The following differential

More information

Math 331 Homework Assignment Chapter 7 Page 1 of 9

Math 331 Homework Assignment Chapter 7 Page 1 of 9 Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a

More information

Relaxation. Ravinder Reddy

Relaxation. Ravinder Reddy Relaxation Ravinder Reddy Relaxation What is nuclear spin relaxation? What causes it? Effect on spectral line width Field dependence Mechanisms Thermal equilibrium ~10-6 spins leads to NMR signal! T1 Spin-lattice

More information