AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT

Size: px
Start display at page:

Download "AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT"

Transcription

1 FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING

2 FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be applied to derive the EoM of the flexible aircraft. However, to get used to the equations, we shall start with the already known EoM of a rigid aircraft (supersonic revision using Newtonian Mechanics). In this part, we will formulate the rigid body aircraft dynamics with Lagrangian Mechanics determining potential and kinetic energy expressions in terms of generalized coordinates. Lagrange s equations shall be enforced - and you will rapidly discover that you will need a symbolic mathematics computational tool for, for example, Mathematica.

3 FLÁIO SILESTRE Consider the following aircraft as a rigid-body Application of Lagrangian Mechanics! is the angular velocity relative to the inertial frame is the aircraft velocity relative to the inertial frame dm is a mass element p is the position of the mass element relative to the centre of mass R is the position of the mass element relative to the origin of the inertial frame x b roll pitch CM yaw z b y b longitudinal axis

4 FLÁIO SILESTRE Application of Lagrangian Mechanics ATTENTION! Observe that every vector may be written in any reference frame! x b roll pitch y b I.e. although! is the angular velocity relative to the inertial frame, is does not necessarily BE WRITTEN in the inertial frame. CM z b yaw longitudinal axis Example: in the body frame! b = p q r

5 FLÁIO SILESTRE Lagrange s Equations Consider a dynamical system of n degrees of freedom. Lagrange s EoM read: = Q j, j q j L is the Lagrangian, defined as L = T Q j are the generalised loads, defined as Q j = f R j U d volume integral distributed force, written in the inertial frame element mass position, written in the inertial frame

6 FLÁIO SILESTRE Degrees of freedom three translations: components of R 0 in the body frame R 0 b = x CM y CM x b roll pitch CM y b longitudinal axis z CM yaw three rotations: Euler angles z b,, 6

7 FLÁIO SILESTRE Kinetic Energy Kinetic energy of a mass element dt = 1 dr dt dr dt dm Integrating in the volume of the vehicle T = 1 dr dt dr dt dm d( ) dt is the operator time derivative relative to the inertial frame however, the vector does not need to be written in this frame!!! 7

8 FLÁIO SILESTRE Kinetic Energy The derivative of R relative to the inertial frame is x b roll pitch y b CM longitudinal axis dr dt = dr 0 dt + dp dt z b yaw = + p t +! p ( ) t is the operator time derivative relative to the body frame 8

9 FLÁIO SILESTRE Kinetic Energy Thus + pt +! p + pt +! p T = 1 dm or expanding T = 1 dm + 1 (! p) (! p) dm + 1 p t p t dm + T T 1 p t dm + T (! p) dm T Let s analyse each term. T + p t (! p) dm T 6 9

10 FLÁIO SILESTRE Kinetic Energy Let s analyse each term. T 1= 1 = 1 = 1 m dm dm 10

11 FLÁIO SILESTRE T = 1 = 1 = 1 = 1 = 1!T (! p) (! p) dm (p!) (p!) dm (X p!) (X p!)dm (X p!) T (X p!)dm = 1!T J! X T p X p dm! J = X p = Kinetic Energy antisymmetric matrix 0 z y z 0 x y x 0 inertia matrix J xx J xy J xz J yx J yy J yz J zx J zy J zz 11

12 FLÁIO SILESTRE Kinetic Energy T = 1 =0 p t p t dm x b roll pitch y b T = p t dm CM yaw longitudinal axis =0 z b T 6 = =0 p t (! p) dm rigid-body approximation 1

13 FLÁIO SILESTRE Kinetic Energy T = (! p) dm x b roll pitch CM y b longitudinal axis = apple (! p) dm yaw =! p dm z b =0 1

14 FLÁIO SILESTRE Final kinetic energy expression (assuming rigid-body) is Kinetic Energy T = 1 m + 1!T J! remember our definition of degrees of freedom q j to apply Lagrange s EoM, energy expressions must be written in terms of these variables and their time-derivatives! Let s now look into the potential energy of the system! 1

15 FLÁIO SILESTRE Potential Energy Potential (gravitational) energy of a mass element du = R G dm G i = 0 0 g Integrating in the volume of the vehicle U = R G dm = R 0 G dm p G dm U 1 U 1

16 FLÁIO SILESTRE U 1 = = R 0 G = mgh R 0 G dm dm x b roll pitch CM y b longitudinal axis yaw U = p G dm z b = =0 p dm G 16

17 FLÁIO SILESTRE Lagrangian Thus L = T U = 1 m + 1!T J! mgh Let s write L in terms of the chosen variables: x CM, y CM, z CM,,, 17

18 FLÁIO SILESTRE Lagrangian = dr 0 dt = R 0 t +! R 0 R 0 b = x CM y CM z CM R 0 b t = ẋ CM ẏ CM ż CM! b = p q r = 1 0 sin 0 cos sin cos 0 sin cos cos 18

19 FLÁIO SILESTRE Lagrangian Besides R 0 i = x 0 y 0 = x 0 y 0 H = L ib x CM y CM z 0 z CM Applying the last equations, L may be written as we want and the left side of Lagrange s EoM can be determined. j = Q j, j =1,...,n don t get scared!!! 19

20 FLÁIO SILESTRE Generalised Force Let s turn to the right side of Lagrange s EoM: Q j = f R j d volume integral distributed force, written in the inertial frame element mass position, written in the inertial frame Q j = ib f x f y f z j ib x CM y CM z CM + L ib x y z 1 A d aerodynamics thrust 0

21 FLÁIO SILESTRE Q j = ib f x f y f z j ib x CM y CM z CM Generalised Force 1 x + L ib y A d z Note that, for the first three dof's: Q xcm = f T x f y L T ib L ib Q ycm = Q zcm = f z f x f y f z f x f y f z T T L T ib L ib L T ib L ib d = F x d = F y d = F z 1

22 FLÁIO SILESTRE Equations of Motion Applying for all dof s, may simplifications will appear btw. left and right side. You shall get equations in: ẍ CM = ÿ CM = z CM = = = = it shouldn t be too complicate to transform back to u = v = ẇ = ṗ = q = ṙ =

CP1 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 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 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

Chapter 4 The Equations of Motion

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

Spacecraft Dynamics and Control

Spacecraft Dynamics and Control Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 16: Euler s Equations Attitude Dynamics In this Lecture we will cover: The Problem of Attitude Stabilization Actuators Newton

More information

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»

More information

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

Phys 7221 Homework # 8

Phys 7221 Homework # 8 Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with

More information

Flight Dynamics & Control Equations of Motion of 6 dof Rigid Aircraft-Kinematics

Flight Dynamics & Control Equations of Motion of 6 dof Rigid Aircraft-Kinematics Flight Dynamic & Control Equation of Motion of 6 dof Rigid Aircraft-Kinematic Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Outline Rotation Matrix Angular Velocity Euler

More information

SOLUTIONS, PROBLEM SET 11

SOLUTIONS, PROBLEM SET 11 SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that

More information

Rotational & Rigid-Body Mechanics. Lectures 3+4

Rotational & 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 information

Part 8: Rigid Body Dynamics

Part 8: Rigid Body Dynamics Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod

More information

6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant

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

Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient

Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained

More information

General Physics I. Lecture 10: Rolling Motion and Angular Momentum.

General Physics I. Lecture 10: Rolling Motion and Angular Momentum. General Physics I Lecture 10: Rolling Motion and Angular Momentum Prof. WAN, Xin (万歆) 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Outline Rolling motion of a rigid object: center-of-mass motion

More information

Flight Dynamics and Control

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

Lagrange Multipliers

Lagrange Multipliers Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each

More information

In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:

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

Optimal Control, Guidance and Estimation. Lecture 16. Overview of Flight Dynamics II. Prof. Radhakant Padhi. Prof. Radhakant Padhi

Optimal Control, Guidance and Estimation. Lecture 16. Overview of Flight Dynamics II. Prof. Radhakant Padhi. Prof. Radhakant Padhi Optimal Control, Guidance and Estimation Lecture 16 Overview of Flight Dynamics II Prof. Radhakant Padhi Dept. of erospace Engineering Indian Institute of Science - Bangalore Point Mass Dynamics Prof.

More information

Lecture (2) Today: Generalized Coordinates Principle of Least Action. For tomorrow 1. read LL 1-5 (only a few pages!) 2. do pset problems 4-6

Lecture (2) Today: Generalized Coordinates Principle of Least Action. For tomorrow 1. read LL 1-5 (only a few pages!) 2. do pset problems 4-6 Lecture (2) Today: Generalized Coordinates Principle of Least Action For tomorrow 1. read LL 1-5 (only a few pages!) 2. do pset problems 4-6 1 Generalized Coordinates The first step in almost any mechanics

More information

From quantum to classical statistical mechanics. Polyatomic ideal gas.

From quantum to classical statistical mechanics. Polyatomic ideal gas. From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,

More information

Lecture 11 Overview of Flight Dynamics I. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 11 Overview of Flight Dynamics I. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 11 Overview of Flight Dynamics I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Point Mass Dynamics Dr. Radhakant Padhi Asst. Professor

More information

Physical Dynamics (SPA5304) Lecture Plan 2018

Physical Dynamics (SPA5304) Lecture Plan 2018 Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle

More information

identify appropriate degrees of freedom and coordinates for a rigid body;

identify appropriate degrees of freedom and coordinates for a rigid body; Chapter 5 Rotational motion A rigid body is defined as a body (or collection of particles) where all mass points stay at the same relative distances at all times. This can be a continuous body, or a collection

More information

UAV Coordinate Frames and Rigid Body Dynamics

UAV Coordinate Frames and Rigid Body Dynamics Brigham Young University BYU ScholarsArchive All Faculty Publications 24-- UAV oordinate Frames and Rigid Body Dynamics Randal Beard beard@byu.edu Follow this and additional works at: https://scholarsarchive.byu.edu/facpub

More information

u (surge) X o p (roll) Body-fixed r o v (sway) w (heave) Z o Earth-fixed X Y Z r (yaw) (pitch)

u (surge) X o p (roll) Body-fixed r o v (sway) w (heave) Z o Earth-fixed X Y Z r (yaw) (pitch) Nonlinear Modelling of Marine Vehicles in Degrees of Freedom Thor I. Fossen and Ola-Erik Fjellstad The Norwegian Institute of Technology Department of Engineering Cybernetics N-0 Trondheim, NORWAY (E-mail:tif@itk.unit.no)

More information

Physical Dynamics (PHY-304)

Physical Dynamics (PHY-304) Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.

More information

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

M2A2 Problem Sheet 3 - Hamiltonian Mechanics MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

More information

COMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.

COMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED. BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section

More information

Applications of Eigenvalues & Eigenvectors

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

2.003 Engineering Dynamics Problem Set 6 with solution

2.003 Engineering Dynamics Problem Set 6 with solution .00 Engineering Dynamics Problem Set 6 with solution Problem : A slender uniform rod of mass m is attached to a cart of mass m at a frictionless pivot located at point A. The cart is connected to a fixed

More information

Video 2.1a Vijay Kumar and Ani Hsieh

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

Dynamics and Control of Rotorcraft

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

Classical Mechanics. Luis Anchordoqui

Classical Mechanics. Luis Anchordoqui 1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

More information

Physics 106b/196b Problem Set 9 Due Jan 19, 2007

Physics 106b/196b Problem Set 9 Due Jan 19, 2007 Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics

More information

Dynamics. 1 Copyright c 2015 Roderic Grupen

Dynamics. 1 Copyright c 2015 Roderic Grupen Dynamics The branch of physics that treats the action of force on bodies in motion or at rest; kinetics, kinematics, and statics, collectively. Websters dictionary Outline Conservation of Momentum Inertia

More information

M E 320 Professor John M. Cimbala Lecture 10

M E 320 Professor John M. Cimbala Lecture 10 M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT

More information

Poincaré (non-holonomic Lagrange) Equations

Poincaré (non-holonomic Lagrange) Equations Department of Theoretical Physics Comenius University Bratislava fecko@fmph.uniba.sk Student Colloqium and School on Mathematical Physics, Stará Lesná, Slovakia, August 23-29, 2010 We will learn: In which

More information

Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013

Rigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013 Rigid body dynamics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Rigid body dynamics October 2013 1 / 16 Multiple point-mass bodies Each mass is

More information

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect

More information

Classical Mechanics Comprehensive Exam Solution

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

LOWELL WEEKLY JOURNAL

LOWELL WEEKLY JOURNAL Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q

More information

Adaptive Trim and Trajectory Following for a Tilt-Rotor Tricopter Ahmad Ansari, Anna Prach, and Dennis S. Bernstein

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

A Newtonian Development of the Mean-Axis Dynamics with Example and Simulation

A Newtonian Development of the Mean-Axis Dynamics with Example and Simulation A Newtonian Development of the Mean-Axis Dynamics with Example and Simulation A THESIS SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Sally Ann Keyes IN PARTIAL FULFILLMENT

More information

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA

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

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

Exercise 1: Inertia moment of a simple pendulum

Exercise 1: Inertia moment of a simple pendulum Exercise : Inertia moment of a simple pendulum A simple pendulum is represented in Figure. Three reference frames are introduced: R is the fixed/inertial RF, with origin in the rotation center and i along

More information

If the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.

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

Problem Goldstein 2-12

Problem Goldstein 2-12 Problem Goldstein -1 The Rolling Constraint: A small circular hoop of radius r and mass m hoop rolls without slipping on a stationary cylinder of radius R. The only external force is that of gravity. Let

More information

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1

AE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1 AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01

More information

Parachute Dynamic Stability and the Effects of Apparent Inertia

Parachute Dynamic Stability and the Effects of Apparent Inertia Parachute Dynamic Stability and the Effects of Apparent Inertia Jason M. Ginn Georgia Institute of Technology, Atlanta, GA, 333 Ian G. Clark Jet Propulsion Laboratory, California Institute of Technology,

More information

PHYS 705: Classical Mechanics. Euler s Equations

PHYS 705: Classical Mechanics. Euler s Equations 1 PHYS 705: Classical Mechanics Euler s Equations 2 Euler s Equations (set up) We have seen how to describe the kinematic properties of a rigid body. Now, we would like to get equations of motion for it.

More information

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY)

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Rigid Body Rotation Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content Introduction Angular Velocity Angular Momentum

More information

Translational and Rotational Dynamics!

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

Quadrotor Modeling and Control

Quadrotor Modeling and Control 16-311 Introduction to Robotics Guest Lecture on Aerial Robotics Quadrotor Modeling and Control Nathan Michael February 05, 2014 Lecture Outline Modeling: Dynamic model from first principles Propeller

More information

Robot Control Basics CS 685

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

Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles.

Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles. » ~ $ ) 7 x X ) / ( 8 2 X 39 ««x» ««! «! / x? \» «({? «» q «(? (?? x! «? 8? ( z x x q? ) «q q q ) x z x 69 7( X X ( 3»«! ( ~«x ««x ) (» «8 4 X «4 «4 «8 X «x «(» X) ()»» «X «97 X X X 4 ( 86) x) ( ) z z

More information

Dynamics of Open Chains

Dynamics of Open Chains Chapter 9 Dynamics of Open Chains According to Newton s second law of motion, any change in the velocity of a rigid body is caused by external forces and torques In this chapter we study once again the

More information

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

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top

Physics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem

More information

Video 3.1 Vijay Kumar and Ani Hsieh

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

AN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016

AN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016 AN INTRODUCTION TO LAGRANGE EQUATIONS Professor J. Kim Vandiver October 28, 2016 kimv@mit.edu 1.0 INTRODUCTION This paper is intended as a minimal introduction to the application of Lagrange equations

More information

Lecture Notes Multibody Dynamics B, wb1413

Lecture Notes Multibody Dynamics B, wb1413 Lecture Notes Multibody Dynamics B, wb1413 A. L. Schwab & Guido M.J. Delhaes Laboratory for Engineering Mechanics Mechanical Engineering Delft University of Technolgy The Netherlands June 9, 29 Contents

More information

Lesson Rigid Body Dynamics

Lesson Rigid Body Dynamics Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body

More information

Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System

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

MODELING AND DYNAMIC ANALYSIS OF A MULTI-JOINT MORPHING AIRCRAFT

MODELING AND DYNAMIC ANALYSIS OF A MULTI-JOINT MORPHING AIRCRAFT MODELING AND DYNAMIC ANALYSIS OF A MULTI-JOINT MORPHING AIRCRAFT By DANIEL T. GRANT A THESIS PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR

More information

Leonhard Euler ( September 1783)

Leonhard Euler ( September 1783) LEONHARD EULER (1707 - SEPTEMBER 1783) Leonhard Euler (1707 - September 1783) BEYOND EQUATIONS Leonhard Euler was born in Basle, Switzerland; he was in fact a born mathematician, who went on to become

More information

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS PHYSICS 44 {CLASSICAL MECHANICS 3 SUMMER FINAL EXAM { SOLUTIONS Problem 1 Consider the following Atwood machine comosed of the two masses m 1 and m and attached to a stand with the hel of a sring of strength

More information

Two Posts to Fill On School Board

Two Posts to Fill On School Board Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83

More information

HW3 Physics 311 Mechanics

HW3 Physics 311 Mechanics HW3 Physics 311 Mechanics FA L L 2 0 1 5 P H Y S I C S D E PA R T M E N T U N I V E R S I T Y O F W I S C O N S I N, M A D I S O N I N S T R U C T O R : P R O F E S S O R S T E F A N W E S T E R H O F

More information

Physics 351, Spring 2018, Homework #9. Due at start of class, Friday, March 30, 2018

Physics 351, Spring 2018, Homework #9. Due at start of class, Friday, March 30, 2018 Physics 351, Spring 218, Homework #9. Due at start of class, Friday, March 3, 218 Please write your name on the LAST PAGE of your homework submission, so that we don t notice whose paper we re grading

More information

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort - 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [

More information

Chapter 1 Lecture 2. Introduction 2. Topics. Chapter-1

Chapter 1 Lecture 2. Introduction 2. Topics. Chapter-1 Chapter 1 Lecture 2 Introduction 2 Topics 1.4 Equilibrium of airplane 1.5 Number of equations of motion for airplane in flight 1.5.1 Degrees of freedom 1.5.2 Degrees of freedom for a rigid airplane 1.6

More information

Lecture 38: Equations of Rigid-Body Motion

Lecture 38: Equations of Rigid-Body Motion Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can

More information

Linearized Equations of Motion!

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 information

MatFys. Week 5, Dec , 2005

MatFys. Week 5, Dec , 2005 MatFys Week 5, Dec. 1-18, 005 Lectures During this week s lectures we shall finish our treatment of classical mechanics and begin the quantum theory part of the course with an introduction to Hilbert spaces.

More information

f x f y or else = m a y

f x f y or else = m a y Rigid Body Dynamics If the aim of kinematics is to describe the body motion, the aim of dynamics is to explain it; the history of mechanics shows that the passage from description to explanation requires

More information

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a Stephen Martin PHYS 10 Homework #1 Question 1: A particle starts at rest and moves along a cycloid whose equation is [ ( ) a y x = ± a cos 1 + ] ay y a There is a gravitational field of strength g in the

More information

21.60 Worksheet 8 - preparation problems - question 1:

21.60 Worksheet 8 - preparation problems - question 1: Dynamics 190 1.60 Worksheet 8 - preparation problems - question 1: A particle of mass m moes under the influence of a conseratie central force F (r) =g(r)r where r = xˆx + yŷ + zẑ and r = x + y + z. A.

More information

2.003 Quiz #1 Review

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

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

AE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept. AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Single-track models of an A-double heavy vehicle combination

Single-track models of an A-double heavy vehicle combination Single-track models of an A-double heavy vehicle combination PETER NILSSON KRISTOFFER TAGESSON Department of Applied Mechanics Division of Vehicle Engineering and Autonomous Systems Vehicle Dynamics Group

More information

Stability and Control

Stability and Control Stability and Control Introduction An important concept that must be considered when designing an aircraft, missile, or other type of vehicle, is that of stability and control. The study of stability is

More information

Minimization of Static! Cost Functions!

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

Introduction to Haptic Systems

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

Physics 5153 Classical Mechanics. Canonical Transformations-1

Physics 5153 Classical Mechanics. Canonical Transformations-1 1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant

More information

Stability of the Lagrange Points, L 4 and L 5

Stability of the Lagrange Points, L 4 and L 5 Stability of the Lagrange Points, L 4 and L 5 Thomas Greenspan January 7, 014 Abstract A proof of the stability of the non collinear Lagrange Points, L 4 and L 5. We will start by covering the basics of

More information

P235 Midterm Examination Prof. Cline

P235 Midterm Examination Prof. Cline P235 Mier Exaination Prof. Cline THIS IS A CLOSED BOOK EXAMINATION. Do all parts of all four questions. Show all steps to get full credit. 7:00-10.00p, 30 October 2009 1:(20pts) Consider a rocket fired

More information

PHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010

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

Chapter 8 Part 1. Attitude Dynamics: Disturbance Torques AERO-423

Chapter 8 Part 1. Attitude Dynamics: Disturbance Torques AERO-423 Chapter 8 Part 1 Attitude Dynamics: Disturbance Torques AEO-43 Types of Disturbance Torques Solar Pressure Dominant torque for geosynchronous satellites Gravity Gradient Can be disturbance or control torque

More information

Physics 6010, Fall Relevant Sections in Text: Introduction

Physics 6010, Fall Relevant Sections in Text: Introduction Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational

More information

Constraints. Noninertial coordinate systems

Constraints. Noninertial coordinate systems Chapter 8 Constraints. Noninertial codinate systems 8.1 Constraints Oftentimes we encounter problems with constraints. F example, f a ball rolling on a flo without slipping, there is a constraint linking

More information

' Liberty and Umou Ono and Inseparablo "

' Liberty and Umou Ono and Inseparablo 3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <

More information

Torque and Rotation Lecture 7

Torque and Rotation Lecture 7 Torque and Rotation Lecture 7 ˆ In this lecture we finally move beyond a simple particle in our mechanical analysis of motion. ˆ Now we consider the so-called rigid body. Essentially, a particle with extension

More information

Lecture Notes for PHY 405 Classical Mechanics

Lecture Notes for PHY 405 Classical Mechanics Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter

More information

ESM 3124 Intermediate Dynamics 2012, HW6 Solutions. (1 + f (x) 2 ) We can first write the constraint y = f(x) in the form of a constraint

ESM 3124 Intermediate Dynamics 2012, HW6 Solutions. (1 + f (x) 2 ) We can first write the constraint y = f(x) in the form of a constraint ESM 314 Intermediate Dynamics 01, HW6 Solutions Roller coaster. A bead of mass m can slide without friction, under the action of gravity, on a smooth rigid wire which has the form y = f(x). (a) Find the

More information

9 Kinetics of 3D rigid bodies - rotating frames

9 Kinetics of 3D rigid bodies - rotating frames 9 Kinetics of 3D rigid bodies - rotating frames 9. Consider the two gears depicted in the figure. The gear B of radius R B is fixed to the ground, while the gear A of mass m A and radius R A turns freely

More information

Lagrangian Dynamics: Derivations of Lagrange s Equations

Lagrangian Dynamics: Derivations of Lagrange s Equations Constraints and Degrees of Freedom 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 4/9/007 Lecture 15 Lagrangian Dynamics: Derivations of Lagrange s Equations Constraints and

More information

Properties of surfaces II: Second moment of area

Properties of surfaces II: Second moment of area Properties of surfaces II: Second moment of area Just as we have discussing first moment of an area and its relation with problems in mechanics, we will now describe second moment and product of area of

More information

Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction

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