Dynamics and Control of Rotorcraft

Size: px
Start display at page:

Download "Dynamics and Control of Rotorcraft"

Transcription

1 Dynamics and Control of Rotorcraft Helicopter Aerodynamics and Dynamics Abhishek Department of Aerospace Engineering Indian Institute of Technology, Kanpur February 3, 2018

2 Overview Flight Dynamics Model 1 Flight Dynamics Model

3 Section 1 Flight Dynamics Model

4 Subsection 1

5 Newton-Euler Equations - 1 F = m a CG + ω V CG + mg ˆk G M = h H dt + ω H M = [I ]{ω} + ω [I ]{ω}

6 Newton Euler Equations - 2 X = m( u + qw rv) + mg sin θ Y = m( v + ru pw) mg sin φ cos θ Z = m(ẇ + pv qu) mg cos φ cos θ L = I xx ṗ I zz (ṙ + pq) (I yy I zz )qr N = I zz ṙ I xz (ṗ qr) (I xx I yy )pq

7 Kinematic Relations Flight Dynamics Model sin θ sin φ φ = p + q + r tan θφ 0 cos θ θ = q cos φ r sin φ ϕ = q sin φ cos φ + r cos φ cos θ

8 Subsection 2

9 Small Perturbation Analysis - 1 Each velocity component is written as some of baseline (trim) value and a time dependent perturbation: u(t) = u 0 + u(t) v(t) = v 0 + v(t) w(t) = w 0 + w(t) p(t) = p 0 + p(t) r(t) = r 0 + r(t)

10 Small Perturbation Analysis - 2 Euler angles and forces and moment are also written in similar way φ(t) = φ 0 + φ(t) θ(t) = θ 0 + φ(t) ϕ(t) = ϕ 0 + ϕ(t) Finally we write forces and moments: X (t) = X 0 + X (t) Y (t) = Y 0 + Y (t)... N(t) = N 0 + N(t) Now we substitute these quantities in Euler equations and kinematic relations.

11 Small Perturbation Analysis - 3 For example: X 0 m + X 0 m = u + q 0w 0 + q 0 w + qw 0 + q w r 0 v 0 r 0 v rv 0 r v + g sin θ 0 cos θ +g cos θ 0 sin θ Small perturbation assumption is made: 2 and HOT are neglected

12 Small Perturbation Analysis - 4 Further the perturbation angle φ, θ are small i.e cos φ 1, sin φ φ etc. Equation becomes: X 0 m + X 0 m = u + q 0w 0 r 0 v 0 + g sin θ 0 + q 0 w r 0 v +w 0 q + g cos θ 0 θ v 0 r Setting perturbations = 0, i.e quantities = 0 corresponds to trim condition. X 0 m = q 0w 0 r 0 v 0 + g sin θ 0 This has to be subtracted from the perturbation equation above.

13 Small Perturbation Analysis - 5 Equation becomes: X m = u + q 0 w r 0 v + w 0 q + g cos θ 0 θ v 0 r Similarly other two force equations are: Y m Z m = v + u 0 r + r 0 u v 0 p p 0 u + g sin θ 0 sin φ 0 θ g cos θ 0 cos φ 0 φ = ẇ + v 0 p + p 0 v u 0 q q 0 u + g sin θ 0 cos φ 0 θ +g cos θ 0 sin φ 0 φ

14 Small Perturbation Analysis - 6 Small perturbation Moment Equilibrium Equations (with dropped): L = I x ṗ I xz ṙ + [(I z I y )r 0 I xz p 0 ]q I xz pq 0 + (I z I y )q 0 M = I y q + [2I xz p 0 + (I x I z )r 0 ]p + [(I x I z )p 0 2I xz r 0 ]r N = I z ṙ I xz ṗ + [(I y I x )p 0 I xz r 0 ]q + (I y I x )q 0 p I xz q 0 r L and M equations are coupled in ṗ and ṙ, therefore can be solved simultaneously to decouple them (as algebraic equation).

15 Small Perturbation Analysis - 7 Small perturbation Kinematic relations: φ = p + tan θ 0 sin φ 0 q + tan θ 0 cos φ 0 r + θ = cos φ 0 q sin φ 0 r ϕ 0 cos θ 0 φ ϕ = sin φ 0 q + cos φ 0 r + ϕ 0 tan θ 0 θ cos θ 0 cos θ 0 ϕ 0 θ cos θ 0 L and M equations are coupled in ṗ and ṙ, therefore can be solved simultaneously to decouple them (as algebraic equation).

16 Introduction to Stability Derivatives - 1 X, Y, Z and L, M, N are functions of a variety of states and controls, for example: X = X(u, v, w, p, q, r, φ, θ, ϕ) Rigid body motion. θ 0, θ 1s, θ 1c, θ t Pilot controls. β, ζ, θ rot Rotor states λ Inflow states Ω Propulsion system states. u FCS Flight control system states. y other states (gust etc.) t time.

17 Introduction to Stability Derivatives - 2 For the time being we will consider only the key parameters, rigid body motion and control inputs. Expand X in a Taylor s series about trimmed flight condition: X = X (u 0, v 0,..., ϕ 0, θ 0,..., (θ 0 ) 0, (θ 1s ) 0, (θ t ) 0 )+ X u (u u 0)+ X v (v v 0)+ Now, X -X (u 0, v 0...) = X X = X u (u u 0) + X v (v v 0) Where, X u is stability derivative, (u u 0) is perturbation u, (v v 0 ) is v

18 Introduction to Stability Derivatives - 3 The partial derivatives are called stability derivatives or control derivative depending on whether denominator is motion quantity or perturbation quantity. These are typically indicated as: X u = X u, Z w = Z M w, θ 1s = M θ1s... etc. Also note: X φ, X θ, X ϕ = 0 Y φ... N φ, N θ... as aerodynamic rotor loads do not change by changing φ or θ or ϕ and keeping everything else constant.

19 Introduction to Stability Derivatives - 4 With above definitions and again dropping for perturbation we can write our equation as: ( X =)X = X u +X v v+x w w+x p p+x q q+x r r+x θ0 θ 0 +X θ1c θ 1c +X θ1s θ 1s = X θt θ t Y= Y u u..... N = N u u + N v v N θt t We can substitute these in small perturbation equation of motion.

20 Introduction to Stability Derivatives - 5 Using u 0 = v cos α cos β, V 0 = v sin β, w 0 = v sin α cos β u X u u (X v + ϕ cos φ 0 cos θ 0 )v (X w + ϕ sin φ 0 cos θ 0 )w X p p (X q v sin α cos β)q X r r v sin βr + g cos θ 0 θ = X θ0 θ 0 + X θ1c θ 1c + X θ1s θ 1s + X θt θ t Final set of 9 equations can be written in a matrix from like this: {ẋ} = [A]{x} + [B]{u} {x} state vector = {u, v, w, p, q, r, φ, θ, ϕ} T {u} control vector = {θ 0, θ 1c, θ 1s, θ t } T A more sophisticated model may have larger number of states involved: {x} = {u, v, w, p, q, r, φ, θ, ϕ, λ 0, λ 1c, λ 1s, β 0, β 1c, β 1s..}

21 Introduction to Stability Derivatives - 6 Depending on the need, the equations can be split in to logitudinal dynamics and lateral-directional dynamics separately or the entire equation can be solved. The set of equations representing longitudinal dynamics for steady level flight case: u ẇ q = θ X u X w X q g cos θ 0 u X θ0 X θ1s ( ) w θ0 Z u Z w Z q g sin θ 0 M u M w M q q θ + Z θ0 M θ0 Z θ1s M θ1s 0 0 θ 1s

22 Introduction to Stability Derivatives - 7 A constant coefficient model is sufficiently accurate for a helicopter, if the motion is slow; i.e if the time required for rotor revolution is small compared with the period of the motion under consideration, i.e perturbation or helicopter body rates ω <1-2 rad/sec [accuracy is acceptable] 1-2< ω <10 rad/s [accuracy becomes questionable] ω >10 rad/s [inaccurate] What is typical rate generated by helicopters during maneuvers: UTTAS pull-up maneuver had highest pitch rate of 30deg/s 0.7 rad/s Choices in FD Constant coefficient stability derivatives. Time dependent stability derivatives Bigger sate vector(w/rotor states, inflow etc) but still constcoefficients.

23 Introduction to Stability Derivatives - 8 How to compute stability derivatives? Analytical approach Take derivatives about trim conditions. Tedious process and dependent on model. Numerical approach Numerical perturbation based calculation. Faster, easier and independent of model. Rotor model can be non-linear and comprehensive. Selecting numerical perturbation is not simple.

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

Robot Dynamics - Rotary Wing UAS: Control of a Quadrotor

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

Topic # Feedback Control

Topic # 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 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

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

Equations of Motion for Micro Air Vehicles

Equations of Motion for Micro Air Vehicles Equations of Motion for Micro Air Vehicles September, 005 Randal W. Beard, Associate Professor Department of Electrical and Computer Engineering Brigham Young University Provo, Utah 84604 USA voice: (80

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

Supplementary Section D: Additional Material Relating to Helicopter Flight Mechanics Models for the Case Study of Chapter 10.

Supplementary Section D: Additional Material Relating to Helicopter Flight Mechanics Models for the Case Study of Chapter 10. Supplementary Section D: Additional Material Relating to Helicopter Flight Mechanics Models for the Case Study of Chapter 1. D1 Nonlinear Flight-Mechanics Models and their Linearisation D1.1 Introduction

More information

Problem 1: Ship Path-Following Control System (35%)

Problem 1: Ship Path-Following Control System (35%) Problem 1: Ship Path-Following Control System (35%) Consider the kinematic equations: Figure 1: NTNU s research vessel, R/V Gunnerus, and Nomoto model: T ṙ + r = Kδ (1) with T = 22.0 s and K = 0.1 s 1.

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

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

AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT

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

Sailing Performance and Maneuverability of a Traditional Ryukyuan Tribute Ship

Sailing Performance and Maneuverability of a Traditional Ryukyuan Tribute Ship sia Navigation Conference 009 ailing Performance and Maneuverability of a Traditional Ryukyuan Tribute hip by Yutaka MUYM (Kanazawa Institute of Technology) Hikaru YGI (Tokai University ) Yutaka TERO (Tokai

More information

CHAPTER 1. Introduction

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

Mechanics of Flight. Warren F. Phillips. John Wiley & Sons, Inc. Professor Mechanical and Aerospace Engineering Utah State University WILEY

Mechanics of Flight. Warren F. Phillips. John Wiley & Sons, Inc. Professor Mechanical and Aerospace Engineering Utah State University WILEY Mechanics of Flight Warren F. Phillips Professor Mechanical and Aerospace Engineering Utah State University WILEY John Wiley & Sons, Inc. CONTENTS Preface Acknowledgments xi xiii 1. Overview of Aerodynamics

More information

Time Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors

Time Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors 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

More information

Pitch Control of Flight System using Dynamic Inversion and PID Controller

Pitch Control of Flight System using Dynamic Inversion and PID Controller Pitch Control of Flight System using Dynamic Inversion and PID Controller Jisha Shaji Dept. of Electrical &Electronics Engineering Mar Baselios College of Engineering & Technology Thiruvananthapuram, India

More information

Rough Plane Analysis. Contents

Rough Plane Analysis. Contents Rough Plane Analysis David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy

More information

EECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators.

EECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators. Name: SID: EECS C128/ ME C134 Final Wed. Dec. 14, 211 81-11 am Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth 1 points total. Problem Points Score 1 16 2 12

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

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

Some effects of large blade deflections on aeroelastic stability

Some effects of large blade deflections on aeroelastic stability 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-839 Some effects of large blade deflections on aeroelastic stability

More information

Linear Flight Control Techniques for Unmanned Aerial Vehicles

Linear Flight Control Techniques for Unmanned Aerial Vehicles Chapter 1 Linear Flight Control Techniques for Unmanned Aerial Vehicles Jonathan P. How, Emilio Frazzoli, and Girish Chowdhary August 2, 2012 1 Abstract This chapter presents an overview of linear flight

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

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 AE89 MS Special Problems Report Space Systems Design Lab (SSDL) Guggenheim School of Aerospace Engineering Georgia Institute of Technology

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

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

Flight Dynamics and Control. Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege

Flight Dynamics and Control. Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege Flight Dynamics and Control Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege Previously on AERO0003-1 We developed linearized equations of motion Longitudinal direction

More information

P A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.

P A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,. ? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>

More information

Investigating the Performance of Adaptive Methods in application to Autopilot of General aviation Aircraft

Investigating the Performance of Adaptive Methods in application to Autopilot of General aviation Aircraft I J C T A, 8(5), 2015, pp 2423-2431 International Science Press Investigating the Performance of Adaptive Methods in application to Autopilot of General aviation Aircraft V Rajesari 1 and L Padma Suresh

More information

Feasibility study of a novel method for real-time aerodynamic coefficient estimation

Feasibility study of a novel method for real-time aerodynamic coefficient estimation Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections 11-4-2010 Feasibility study of a novel method for real-time aerodynamic coefficient estimation Phillip Gurbacki

More information

Robot Dynamics Rotorcrafts: Dynamic Modeling of Rotorcraft & Control

Robot Dynamics Rotorcrafts: Dynamic Modeling of Rotorcraft & Control Robot Dynamics Rotorcrafts: Dynamic Modeling of Rotorcraft & Control 5-85- V Marco Hutter, Roland Siegwart and Thomas Stastny Robot Dynamics - Rotary Wing UAS: Propeller Analysis and Dynamic Modeling 7..5

More information

DYNAMIC MODELING AND STABILITY ANALYSIS OF A HIGH ALTITUDE AIRSHIP

DYNAMIC MODELING AND STABILITY ANALYSIS OF A HIGH ALTITUDE AIRSHIP DYNAMIC MODELING AND STABILITY ANALYSIS OF A HIGH ALTITUDE AIRSHIP Ramesh Hun 1, Nandan K Sinha 2, Debasish Chakraborti 3 and P.Balasubramanian 4 1 Research Scholar and Corresponding author, 2 Professor,

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

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15 Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Kinematics Semester 1, 2014-15 1 / 15 Introduction The kinematic quantities used are: position r, linear velocity

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

Kinematics. Basilio Bona. October DAUIN - Politecnico di Torino. Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October / 15

Kinematics. Basilio Bona. October DAUIN - Politecnico di Torino. Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October / 15 Kinematics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October 2013 1 / 15 Introduction The kinematic quantities used are: position r,

More information

Fundamentals of Airplane Flight Mechanics

Fundamentals of Airplane Flight Mechanics David G. Hull Fundamentals of Airplane Flight Mechanics With 125 Figures and 25 Tables y Springer Introduction to Airplane Flight Mechanics 1 1.1 Airframe Anatomy 2 1.2 Engine Anatomy 5 1.3 Equations of

More information

Visual Servoing for a Quadrotor UAV in Target Tracking Applications. Marinela Georgieva Popova

Visual Servoing for a Quadrotor UAV in Target Tracking Applications. Marinela Georgieva Popova Visual Servoing for a Quadrotor UAV in Target Tracking Applications by Marinela Georgieva Popova A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate

More information

DIFFERENTIAL KINEMATICS. Geometric Jacobian. Analytical Jacobian. Kinematic singularities. Kinematic redundancy. Inverse differential kinematics

DIFFERENTIAL KINEMATICS. Geometric Jacobian. Analytical Jacobian. Kinematic singularities. Kinematic redundancy. Inverse differential kinematics DIFFERENTIAL KINEMATICS relationship between joint velocities and end-effector velocities Geometric Jacobian Analytical Jacobian Kinematic singularities Kinematic redundancy Inverse differential kinematics

More information

The Role of Zero Dynamics in Aerospace Systems

The Role of Zero Dynamics in Aerospace Systems The Role of Zero Dynamics in Aerospace Systems A Case Study in Control of Hypersonic Vehicles Andrea Serrani Department of Electrical and Computer Engineering The Ohio State University Outline q Issues

More information

UNCLASSIFIED. Control of Spinning Symmetric Airframes. Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734

UNCLASSIFIED. Control of Spinning Symmetric Airframes. Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734 Control of Spinning Symmetric Airframes Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734 Several missiles and projectiles have periods during their flight where the

More information

Robot Dynamics Fixed-wing UAVs: Dynamic Modeling and Control

Robot Dynamics Fixed-wing UAVs: Dynamic Modeling and Control Robot Dynamics Fixed-wing UAVs: Dynamic Modeling and Control 151-0851-00 V Marco Hutter, Roland Siegwart, and Thomas Stastny 05.12.2017 1 Contents Fixed-wing UAVs 1. ntroduction 2. Aerodynamic Basics 3.

More information

The spectacle of a rolling aircraft close to the tarmac invariably draws a standing ovation during

The spectacle of a rolling aircraft close to the tarmac invariably draws a standing ovation during Approximating the rolling aircraft J G Manathara & S Pradeep Department of Aerospace Engineering, Indian Institute of Science, Bangalore A R T I C L E I N F O Unpublished manuscript. A B S T R A C T This

More information

Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective

Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective Andrea Serrani Department of Electrical and Computer Engineering Collaborative Center for Control Sciences The Ohio State University

More information

144 Chapter 3. Second Order Linear Equations

144 Chapter 3. Second Order Linear Equations 144 Chapter 3. Second Order Linear Equations PROBLEMS In each of Problems 1 through 8 find the general solution of the given differential equation. 1. y + 2y 3y = 0 2. y + 3y + 2y = 0 3. 6y y y = 0 4.

More information

Lecture No. # 09. (Refer Slide Time: 01:00)

Lecture No. # 09. (Refer Slide Time: 01:00) Introduction to Helicopter Aerodynamics and Dynamics Prof. Dr. C. Venkatesan Department of Aerospace Engineering Indian Institute of Technology, Kanpur Lecture No. # 09 Now, I just want to mention because

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state

More information

DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS

DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS DEVIL PHYSICS BADDEST CLASS ON CAMPUS IB PHYSICS OPTION B-1A: ROTATIONAL DYNAMICS Essential Idea: The basic laws of mechanics have an extension when equivalent principles are applied to rotation. Actual

More information

Frequency Domain System Identification for a Small, Low-Cost, Fixed-Wing UAV

Frequency Domain System Identification for a Small, Low-Cost, Fixed-Wing UAV Frequency Domain System Identification for a Small, Low-Cost, Fixed-Wing UAV Andrei Dorobantu, Austin M. Murch, Bernie Mettler, and Gary J. Balas, Department of Aerospace Engineering & Mechanics University

More information

GUST RESPONSE ANALYSIS IN THE FLIGHT FORMATION OF UNMANNED AIR VEHICLES

GUST RESPONSE ANALYSIS IN THE FLIGHT FORMATION OF UNMANNED AIR VEHICLES 28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES GUST RESPONSE ANALYSIS IN THE FLIGHT FORMATION OF UNMANNED AIR VEHICLES Yoshinobu Inada*, Kyohei Ohta*, and Hideaki Takanobu** *Tokai University,

More information

MANY BILLS OF CONCERN TO PUBLIC

MANY BILLS OF CONCERN TO PUBLIC - 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -

More information

Simultaneous-Perturbation-Stochastic-Approximation Algorithm for Parachute Parameter Estimation

Simultaneous-Perturbation-Stochastic-Approximation Algorithm for Parachute Parameter Estimation JOURNAL OF AIRCRAFT Vol. 4, No. 5, September October 005 Simultaneous-Perturbation-Stochastic-Approximation Algorithm for Parachute Parameter Estimation Govindarajan Kothandaraman and Mario A. Rotea Purdue

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

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

Modelling the Dynamic Response of a Morphing Wing with Active Winglets

Modelling the Dynamic Response of a Morphing Wing with Active Winglets AIAA Atmospheric Flight Mechanics Conference and Exhibit 20-23 August 2007, Hilton Head, South Carolina AIAA 2007-6500 Modelling the Dynamic Response of a Morphing Wing with Active Winglets N. Ameri, M.

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

Robust Nonlinear Design of Three Axes Missile Autopilot via Feedback Linearization

Robust Nonlinear Design of Three Axes Missile Autopilot via Feedback Linearization Robust Nonlinear Design of Three Axes Missile Autopilot via Feedback Linearization Abhijit Das, Ranajit Das and Siddhartha Mukhopadhyay, Amit Patra 1 1 Abstract The nonlinearity and coupling of the missile

More information

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total

More information

Dynamic modeling and control system design for tri-rotor UAV

Dynamic modeling and control system design for tri-rotor UAV Loughborough University Institutional Repository Dynamic modeling and control system design for tri-rotor UAV This item was submitted to Loughborough University's Institutional Repository by the/an author.

More information

Aggressive Maneuvering Flight Tests of a Miniature Robotic Helicopter

Aggressive Maneuvering Flight Tests of a Miniature Robotic Helicopter Aggressive Maneuvering Flight Tests of a Miniature Robotic Helicopter Vladislav Gavrilets, Ioannis Martinos, Bernard Mettler, and Eric Feron Massachusetts Institute of Technology, Cambridge MA 2139, USA

More information

Mathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016

Mathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016 Mathematics 36 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 9 and 2, 206 Every rational function (quotient of polynomials) can be written as a polynomial

More information

QUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. x 2. 1 (1 + 2x 2 /3!) ] x 2 1 2x 2 /3

QUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. x 2. 1 (1 + 2x 2 /3!) ] x 2 1 2x 2 /3 QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) Keeping only the lowest power of x needed, we find 1 x 1 sin x = 1 x 1 (x x 3 /6...) = 1 ) 1 (1 x 1 x /3 = 1 [ 1 (1 + x /3!)

More information

CIS 4930/6930: Principles of Cyber-Physical Systems

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

FINITE-STATE DYNAMIC WAKE INFLOW MODELLING FOR COAXIAL ROTORS

FINITE-STATE DYNAMIC WAKE INFLOW MODELLING FOR COAXIAL ROTORS FINITE-STATE DYNAMIC WAKE INFLOW MODELLING FOR COAXIAL ROTORS Felice Cardito, Riccardo Gori, Giovanni Bernardini, Jacopo Serafini, Massimo Gennaretti Department of Engineering, Roma Tre University, Rome,

More information

Direct spatial motion simulation of aircraft subjected to engine failure

Direct spatial motion simulation of aircraft subjected to engine failure Direct spatial motion simulation of aircraft subjected to engine failure Yassir ABBAS*,1, Mohammed MADBOULI 2, Gamal EL-BAYOUMI 2 *Corresponding author *,1 Aeronautical Engineering Department, Engineering

More information

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67

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

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

PARAMETER ESTIMATION AND FLIGHT PATH RECONSTRUCTION USING OUTPUT-ERROR METHOD

PARAMETER ESTIMATION AND FLIGHT PATH RECONSTRUCTION USING OUTPUT-ERROR METHOD 24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES PARAMETER ESTIMATION AND FLIGHT PATH RECONSTRUCTION USING OUTPUT-ERROR METHOD Luiz C. S. Góes*, Elder M. Hemerly**, Benedito C. O. Maciel* *Dept.

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

Power system modelling under the phasor approximation

Power system modelling under the phasor approximation ELEC0047 - Power system dynamics, control and stability Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct October 2018 1 / 16 Electromagnetic transient vs. phasor-mode simulations

More information

review To find the coefficient of all the terms in 15ab + 60bc 17ca: Coefficient of ab = 15 Coefficient of bc = 60 Coefficient of ca = -17

review To find the coefficient of all the terms in 15ab + 60bc 17ca: Coefficient of ab = 15 Coefficient of bc = 60 Coefficient of ca = -17 1. Revision Recall basic terms of algebraic expressions like Variable, Constant, Term, Coefficient, Polynomial etc. The coefficients of the terms in 4x 2 5xy + 6y 2 are Coefficient of 4x 2 is 4 Coefficient

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

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

A Free Wake Linear Inflow Model Extraction Procedure for Rotorcraft Analysis

A Free Wake Linear Inflow Model Extraction Procedure for Rotorcraft Analysis A Free Wake Linear Inflow Model Extraction Procedure for Rotorcraft Analysis Jeffrey D. Keller Senior Associate Robert M. McKillip, Jr. Senior Associate Continuum Dynamics, Inc. Ewing, NJ, USA Daniel A.

More information

Physics 170 Week 9 Lecture 2

Physics 170 Week 9 Lecture 2 Physics 170 Week 9 Lecture 2 http://www.phas.ubc.ca/ gordonws/170 Physics 170 Week 9 Lecture 2 1 Textbook Chapter 1: Section 1.6 Physics 170 Week 9 Lecture 2 2 Learning Goals: We will solve an example

More information

Hamiltonian. March 30, 2013

Hamiltonian. March 30, 2013 Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................

More information

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

A GLOBAL SLIDING MODE CONTROL WITH PRE-DETERMINED CONVERGENCE TIME DESIGN FOR REUSABLE LAUNCH VEHICLES IN REENTRY PHASE

A GLOBAL SLIDING MODE CONTROL WITH PRE-DETERMINED CONVERGENCE TIME DESIGN FOR REUSABLE LAUNCH VEHICLES IN REENTRY PHASE IAA-AAS-DyCoSS-4 -- A GLOBAL SLIDING MODE CONTROL WITH PRE-DETERMINED CONVERGENCE TIME DESIGN FOR REUSABLE LAUNCH VEHICLES IN REENTRY PHASE L. Wang, Y. Z. Sheng, X. D. Liu, and P. L. Lu INTRODUCTION This

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

A. H. Hall, 33, 35 &37, Lendoi

A. H. Hall, 33, 35 &37, Lendoi 7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9

More information

SCHEME OF BE 100 ENGINEERING MECHANICS DEC 2015

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

Rotor reference axis

Rotor reference axis Rotor reference axis So far we have used the same reference axis: Z aligned with the rotor shaft Y perpendicular to Z and along the blade (in the rotor plane). X in the rotor plane and perpendicular do

More information

Satellite attitude control system simulator

Satellite attitude control system simulator Shock and Vibration 5 (28) 395 42 395 IOS Press Satellite attitude control system simulator G.T. Conti and L.C.G. Souza National Institute for Space Research, INPE, Av. dos Astronautas 758, 2227-, S ão

More information

Attitude Dynamics Aims and Objectives Euler Equations of Rotational Motion

Attitude Dynamics Aims and Objectives Euler Equations of Rotational Motion 13 Attitude Dynamics 13.1 Aims and Objectives To present the universal rotational dynamics model applicable to all aerospace vehicles, emphasizing the commonality between the stability and control characteristics

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

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

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

Alternative Expressions for the Velocity Vector Velocity restricted to the vertical plane. Longitudinal Equations of Motion

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

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response .. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........

More information

Engineering Mechanics

Engineering Mechanics 2019 MPROVEMENT Mechanical Engineering Engineering Mechanics Answer Key of Objective & Conventional Questions 1 System of forces, Centoriod, MOI 1. (c) 2. (b) 3. (a) 4. (c) 5. (b) 6. (c) 7. (b) 8. (b)

More information

CHAPTER 72 AREAS UNDER AND BETWEEN CURVES

CHAPTER 72 AREAS UNDER AND BETWEEN CURVES CHAPTER 7 AREAS UNDER AND BETWEEN CURVES EXERCISE 8 Page 77. Show by integration that the area of the triangle formed by the line y, the ordinates and and the -ais is 6 square units. A sketch of y is shown

More information

Aircraft Performance, Stability and control with experiments in Flight. Questions

Aircraft Performance, Stability and control with experiments in Flight. Questions Aircraft Performance, Stability and control with experiments in Flight Questions Q. If only the elevator size of a given aircraft is decreased; keeping horizontal tail area unchanged; then the aircraft

More information

Torque/Rotational Energy Mock Exam. Instructions: (105 points) Answer the following questions. SHOW ALL OF YOUR WORK.

Torque/Rotational Energy Mock Exam. Instructions: (105 points) Answer the following questions. SHOW ALL OF YOUR WORK. AP Physics C Spring, 2017 Torque/Rotational Energy Mock Exam Name: Answer Key Mr. Leonard Instructions: (105 points) Answer the following questions. SHOW ALL OF YOUR WORK. (22 pts ) 1. Two masses are attached

More information

16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity

16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity 6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =

More information