A note on stability of motion of a projectile

Size: px
Start display at page:

Download "A note on stability of motion of a projectile"

Transcription

1 Sādhanā, Vol. 6, Part, August 00, pp Printed in India A note on stability of motion of a projectile SDNAI Institute of Armament Technology, Pune 05, India insdnaik@yahoo.com MS received 6 February 000; revised March 00 Abstract. A projectile is stabilised using either gyroscopic or aerodynamic stability. But subcalibre projectiles with sabot have both spin and fins. Separate stability criteria are researched generally for each type of projectile. In this paper a stability criterion which can be used for all such bodies has been developed through the Liapunov second method. eywords. fins. Projectiles; aerodynamic stability; potential function; spin; sabot;. Introduction The basic class of problems in stability of an axis-symmetric body has been discussed by many researchers. Before 930, the condition 0 <σ < dominated every design in flight dynamics. Fowler attempted the stability of high angle of fire through his β-type equations. As awareness of the implications of mathematical modelling to investigate the motion in its full entirety developed during the Second World War, many models were made and linear motion was developed in its completeness. ebby 98 investigated nonlinear motion where the parameter is positive and established that if qs > 0, the ensuing motion is bounded. Rath & Sharma 965 have derived the McShane Murphy stability condition for slowly yawing motion from Fowler s dynamical equations generalised by Rath taking into account the complete aerodynamic force system Nielsen & Synge96. The motion of a nonlinear Lock Fowler missile under the same conditon using the Routh Hurwitz criterion has been discussed by Rath & Namboodiri 980. The dynamical motion of an axi-symmetric projectile in the cross plane moving in the atmosphere Bakshi & Sharma 988 is governed by Ẍ + I + JẊ+ 3 I JX = 0, [ ] 0 v u where J =,X =, and 0 w u = ε C A f g t gx l u, = B ς + ες f p + g p t, A list of symbols is given at the end of the paper 379

2 380 S D Naik 3 = ε B ς f p + g t =εςb C A f + g p t gx l + + g x l u u B gx l ς u ε t ες gx l CA g ε f u, gx l f u p + g p t. For the spin-stabilized projectile 3 < 0 while for aerodynamically stable projectiles 3 > 0. For a subcalibre projectile with sabot, none of these means can be used simply as a stability criterion, as a spin-stabilised body requires high angular spin and finned stabilised systems need large fins and these fins disturb the sabot design. The better way out is obtained by giving some percentage of spin and smaller fins. Obviously for such a body 3 is small with sign indefiniteness. An attempt has been made here to discuss the stability criterion for such a body. It also gives the stability for spin-stabilised as well as fin-stabilised bodies.. Stability criterion The P-method Sharma 977 has been used here to obtain the generating functions V and V. Let N Ẋ + PX, be the generating function. Here P = P +P,P is a symmetric and P a skew-symmetric matrix. The inner product of and gives [ Ẋ t Ẋ + Ẋ t PX+X t P t Ẋ+X t 3 +Q ] Ẋ = Ẋ t P +P t Ẋ+Ẋ + t I + J P + J + Q X + X + X t I + J P + J + Q t Ẋ X t 3 I + J P + P t 3 I J X, 3 which is of the form V. = V. We have introduced an arbitrary matrix Q to develop suitable quadratic forms. The motion is said to be asymptotically stable if V>0 and V <0. Now V and V can be expressed as V = Ẋ t + Xt P t Ẋ+ PX and V = +X t 3 +Q Pt P X, {[ Ẋ t P + P t + X t I + J P + J + Q t P + P t ][ + P +P t + t Ẋ P + P t + I + J P + J + Q ] [ X X t I + J P + J + Q t P + P t + I + J P + J + Q 3 I + J P + P t 3 I + J ]} t X. 5

3 Stability of motion of a projectile 38 The asymptotic stability is achieved if V<0 and and V <0. That is 3 + Q P t P 0, 6 3 I + J P + P t 3 I + J t I + J P + J + Q t P + P t + I + J P + J + Q > 0. 7 The values of P and Q are arbitrary. One choice of these values Naik 000 is P = I+ J, Q = 3 + P t P I, reduce these conditions to a single condition, > 0. 8 Another choice for P and Q, one can get by generalizing P and Q as P = α J βi, 9 and Q = β + α I. 0 These values are selected as P is expressed as P = p I + p J and Q is diagonal. α, β are scalar quantities. These values of P and Q reduce V and V to [ V = Ẋ t + X t α J βi Ẋ + α J βi +X t α + β ] β + α X, V = β [Ẋ t X t Ẋ + J α β+ β J α β + β X α β + Asymptotic stability is guarranteed for V>0 and α + β β + 3 β + α β ], V <0 and is obtained if > 0, 3a β + α β < 0, 3b

4 38 S D Naik α β + + β 3 β+α >0. 3c α and β can be arbitrary. In particular, α =,β =0 reduce 3 to < 0 a > 0, b 3 + > 0, c which are conditions for spin-stability. Similarly α = 0,β =, reduce the above conditions 3 to < 0, < 0, 3 > 0. 5a 5b 5c These are the conditions for fin-stability. For a body which is stabilised with small spin and fins, one can select α and β that are arbitrary but real and nonzero. A linear relation between α and β has been attempted here. On examining various values, the one that is applicable is α β =. This reduces the stability conditions to β β < 0, + β β β + β > 0, 6a 6b 3 β+ β >0. 6c β + The value of β can be selected in such a way that V>0. Let β = + β + + β The condition for V>0is β > 0. β is a quadratic in β and has two roots. Its maximum value should provide a practical design situation and is obtained when the two roots are equal. β max = [ + ], 7 β > 0, implies [ + ] > 0. 8

5 Stability of motion of a projectile 383 It is same as 6a. The last term in the above equation appears through the interaction of spin and fins combination stability. This term being positive decreases the spin stability situation. This arises due to the coupling in magnus and damping coefficients. We need a strongly damping force due to spin appearing in. For a subcalibre projectile which is stabilized with spin as well as fins, condition 8 can be called the dynamic stability condition. The corresponding condition for spinning bodies is If + 3 0, the above condition is always satisfied. This is static stability for a spinning body. In particular, if we ignore damping parameter, the dynamic stability 8 becomes This contains a parameter µ = which is very important for the design. For this, β max =, gives the fin stability. With this β max defined as 7, the other two conditions for asymptotic stability coming from V <0 are + + < 0, The condition 0 for a marginally stable projectile, 3 = 0, and a damped system < 0 implies that, the magnus coefficient should be negative, a situation which is a generally a design situation. Hence this result is likely to be met. This negative magnus coefficient also strongly satisfies 9. Thus a negative magnus coefficient with a marginal value as sign indefinite is necessary for stability of a subcalibre projectile. 3. Verification Let us consider a projectile with the following data. ρ =.5, m = 3.6, d = 0.08, p = 5, u = 7, l = 0.86, C D = 0.35, f =.9, g =.5, g = 57.6, g p = 5.86, I xx = 3.66, I yy = To verify the stability criteria, we have to calculate the parameters,, 3 and first. Here ε = ρsd m, t = I yy md, B = I xx Iyy and ς = pl u. Generally ε and x are small and ignored. The values of i s are simplified to =

6 38 S D Naik = = = It can be further observed that a All the conditions 8, 8,9 and 0 are satisfied. b The spin required for stability of this projectile is 0 revs. Hence the projectile is over-stabilized as actually the spin applied is 5 revs. It is thus concluded that sufficient stability criterion meets the technology requirements. c To stabilize this projectile only with the help of spin,the requirement is 387 revs. The spin can be reduced by introducing small fins in the design. In this design the energy saved because of the fins is 85% of that of the spin-stabilised projectile, which can be utilised in any other form.. Conclusions A simple analysis has provided very important design critera for a subcalibre projectile. These are as below. The spin-stability is determined by fins through the magnus parameter. This will be smaller if fins are smaller which is the case for subcalibre projectile. The magnus moment has to be negative. 3 The negativeness of the magnus moment damps the oscillations faster as seen from 9. List of symbols C D drag force coefficient; d diameter of the projectile; f normal force coefficient; g aerodynamic moment coefficient; g damping moment coefficient; g p magnus moment coefficient; l total length of the projectile; I xx moment of inertia of the projectile about X-axis; I yy moment of inertia of the projectile about Y -axis; damping parameter; spin parameter; 3 conservative force parameter; magnus parameter; m mass of the projectile ; p constant spin of projectile; ρ density of air; u avarage speed of the projectile. References Bakshi S, Sharma C 988 A fresh look at the free flight dynamics of spinning projectile. Sadhana : 3 337

7 Stability of motion of a projectile 385 ebby C H 98 The stability of the undamped motion of a projectile subject to a yawing moment of the Fowler-Lock type. Armament Research Establishment Report No. 388 Naik S D 000 Stability criterion for a finned spinning projectile. Int. J. Differential Eqnations Dyn. Syst. 8: 5 63 Namboodiri A V 980 On the dynamic response of a rolling missile due to a nonlinear aerodynamic moment of the Lock Fowler type. Ph D thesis, Department of Mathematics, University of Poona, Pune Nielsen L, Synge J L 96 On the motion of a spinning shell. Q. J. Appl. Math. : 0 06 Rath P C, Sharma J P 965 On the stability of a spinning projectile. Des. Sci. J. 5: Sharma C 977 On the dynamic stability of a missile. Des. Sci. J. 7: 9 3

Research of Flight Characteristics of Rod-Type Projectile with Triangular Cross-section. Dr. Wenjun Yi, Prof. Xiaobing Zhang, Prof.

Research of Flight Characteristics of Rod-Type Projectile with Triangular Cross-section. Dr. Wenjun Yi, Prof. Xiaobing Zhang, Prof. Research of Flight Characteristics of Rod-Type Projectile with Triangular Cross-section Dr. Wenjun Yi, Prof. Xiaobing Zhang, Prof. Jianping Qian Ballistic Research Laboratory of China Outline INTRODUCTION

More information

PERFORMANCE ANALYSIS OF ISL S GUIDED SUPERSONIC PROJECTILE. Pierre Wey

PERFORMANCE ANALYSIS OF ISL S GUIDED SUPERSONIC PROJECTILE. Pierre Wey 3 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16- APRIL 7 PERFORMANCE ANALYSIS OF ISL S GUIDED SUPERSONIC PROJECTILE Pierre Wey French-German Research Institute of Saint-Louis (ISL) P.O.

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

Theoretical and experimental research of supersonic missile ballistics

Theoretical and experimental research of supersonic missile ballistics BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 63, No. 1, 015 DOI: 10.1515/bpasts-015-007 Theoretical and experimental research of supersonic missile ballistics B. ZYGMUNT 1, K. MOTYL

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls Problems for Qualifying Exam - Spring 2014 Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function

More information

Control Systems. Internal Stability - LTI systems. L. Lanari

Control Systems. Internal Stability - LTI systems. L. Lanari Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable

More information

IMPROVEMENT OF ARTILLERY PROJECTILE ACCURACY. C.Grignon, R.Cayzac, S.Heddadj

IMPROVEMENT OF ARTILLERY PROJECTILE ACCURACY. C.Grignon, R.Cayzac, S.Heddadj RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-0 APRIL 007 IMPROVEMENT OF ARTILLERY PROJECTILE ACCURACY C.Grignon, R.Cayzac, S.Heddadj Nexter Munitions Ballistics and Flight Control Department

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

MAGNUS INSTABILITIES AND MODELING FOR A 12.7 mm PROJECTILE WITH ROLL DAMPING FINLETS FROM FREE-FLIGHT TESTS AT MACH 1.7

MAGNUS INSTABILITIES AND MODELING FOR A 12.7 mm PROJECTILE WITH ROLL DAMPING FINLETS FROM FREE-FLIGHT TESTS AT MACH 1.7 EB12 19th International Symposium of Ballistics, 7 11 May 2001, Interlaken, Switzerland MAGNUS INSTABILITIES AND MODELING FOR A 12.7 mm PROJECTILE WITH ROLL DAMPING FINLETS FROM FREE-FLIGHT TESTS AT MACH

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

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

Chapter 2: Unconstrained Extrema

Chapter 2: Unconstrained Extrema Chapter 2: Unconstrained Extrema Math 368 c Copyright 2012, 2013 R Clark Robinson May 22, 2013 Chapter 2: Unconstrained Extrema 1 Types of Sets Definition For p R n and r > 0, the open ball about p of

More information

Aerodynamic jump: A short range view for long rod projectiles

Aerodynamic jump: A short range view for long rod projectiles 239 Aerodynamic jump: A short range view for long rod projectiles Mark Bundy US Army Research Laboratory, AMSRL-WM-BC, Aberdeen Proving Ground, MD 21005-5066, USA It is shown that aerodynamic jump for

More information

Investigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration

Investigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration Investigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration TONI AXELSSON Master s Thesis at Aerospace Engineering Supervisor: Arne Karlsson Examiner: Arne Karlsson ISSN 1651-7660

More information

A Simple Design Approach In Yaw Plane For Two Loop Lateral Autopilots

A Simple Design Approach In Yaw Plane For Two Loop Lateral Autopilots A Simple Design Approach In Yaw Plane For Two Loop Lateral Autopilots Jyoti Prasad Singha Thakur 1, Amit Mukhopadhyay Asst. Prof., AEIE, Bankura Unnayani Institute of Engineering, Bankura, West Bengal,

More information

Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis

Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis J. Alexandre dit Sandretto, A. Chapoutot and O. Mullier U2IS, ENSTA ParisTech SYNCOP April 11, 2015 Closed-loop

More information

ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT

ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1599 1604 c World Scientific Publishing Company ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT KEVIN BARONE and SAHJENDRA

More information

Computational angular oscillatory motion for a.50 caliber bullet type firing from a helicopter machine gun

Computational angular oscillatory motion for a.50 caliber bullet type firing from a helicopter machine gun Journal of Computations & Modelling, vol.2, no.4, 2012, 93-108 ISSN: 1792-7625 (print), 1792-8850 (online) Scienpress Ltd, 2012 Computational angular oscillatory motion for a.50 caliber bullet type firing

More information

Modeling & Control of Hybrid Systems Chapter 4 Stability

Modeling & Control of Hybrid Systems Chapter 4 Stability Modeling & Control of Hybrid Systems Chapter 4 Stability Overview 1. Switched systems 2. Lyapunov theory for smooth and linear systems 3. Stability for any switching signal 4. Stability for given switching

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

Kinetic Energy of Rolling

Kinetic Energy of Rolling Kinetic Energy of Rolling A solid disk and a hoop (with the same mass and radius) are released from rest and roll down a ramp from a height h. Which one is moving faster at the bottom of the ramp? A. they

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov

More information

Contribution of Airplane design parameters on Roll Coupling اي داءالبارامترات التصميميه للطائره على ازدواج الحركي

Contribution of Airplane design parameters on Roll Coupling اي داءالبارامترات التصميميه للطائره على ازدواج الحركي International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:06 7 Contribution of Airplane design parameters on Roll Coupling اي داءالبارامترات التصميميه للطائره على ازدواج الحركي

More information

A Six Degree of Freedom Trajectory Analysis of Spin- Stabilized Projectiles

A Six Degree of Freedom Trajectory Analysis of Spin- Stabilized Projectiles nternational onference of omputational Methods in Sciences and Engineering 007,MSE 007, orfu, Greece, 5-30 Sep. 007 A Six egree of Freedom Trajectory Analysis of Spin- Stabilized Projectiles imitrios N.

More information

Dynamical Systems & Lyapunov Stability

Dynamical Systems & Lyapunov Stability Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence

More information

Sample Solutions of Assignment 10 for MAT3270B

Sample Solutions of Assignment 10 for MAT3270B Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of

More information

Vibrations of stretched damped beams under non-ideal boundary conditions

Vibrations of stretched damped beams under non-ideal boundary conditions Sādhanā Vol. 31, Part 1, February 26, pp. 1 8. Printed in India Vibrations of stretched damped beams under non-ideal boundary conditions 1. Introduction HAKAN BOYACI Celal Bayar University, Department

More information

ECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008

ECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008 ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are

More information

EE363 Automatic Control: Midterm Exam (4 problems, 90 minutes)

EE363 Automatic Control: Midterm Exam (4 problems, 90 minutes) EE363 Automatic Control: Midterm Exam (4 problems, 90 minutes) ) Block diagram simplification (0 points). Simplify the following block diagram.,i.e., find the transfer function from u to y. Your answer

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

Course syllabus Engineering Mechanics - Dynamics

Course syllabus Engineering Mechanics - Dynamics Course syllabus Engineering Mechanics - Dynamics COURSE DETAILS Type of study programme Study programme Course title Course code ECTS (Number of credits allocated) Course status Year of study Course Web

More information

A Study of Firing Sidewise from an Airplane

A Study of Firing Sidewise from an Airplane Applied Mathematical Sciences, Vol. 4, 2010, no. 1, 1-12 A Study of Firing Sidewise from an Airplane Dimitrios N. Gkritzapis, Dionissios P. Margaris Fluid Mechanics Laboratory (FML) Mechanical Engineering

More information

Accurate Computation of the Trajectory of the Spin and Fin-Stabilized Projectiles

Accurate Computation of the Trajectory of the Spin and Fin-Stabilized Projectiles INTERNATIONAL JOURNAL OF MATHEMATIS AN OMPUTERS IN SIMULATION Accurate omputation of the Trajectory of the Spin and Fin-Stabilized Projectiles imitrios N. Gkritzapis, Elias E. Panagiotopoulos Abstract

More information

The Simulation of Wraparound Fins Aerodynamic Characteristics

The Simulation of Wraparound Fins Aerodynamic Characteristics The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:

More information

Designing Information Devices and Systems I Spring 2015 Note 3

Designing Information Devices and Systems I Spring 2015 Note 3 EECS 16A Designing Information Devices and Systems I Spring 2015 Note 3 Lecture notes by Christine Wang (01/27/2015) Introduction to Vectors Remark. Often, vectors are represented as letters in boldface(x),

More information

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation

More information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

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

Parametrically Excited Vibration in Rolling Element Bearings

Parametrically Excited Vibration in Rolling Element Bearings Parametrically Ecited Vibration in Rolling Element Bearings R. Srinath ; A. Sarkar ; A. S. Sekhar 3,,3 Indian Institute of Technology Madras, India, 636 ABSTRACT A defect-free rolling element bearing has

More information

Computation of CPA Lyapunov functions

Computation of CPA Lyapunov functions Computation of CPA Lyapunov functions (CPA = continuous and piecewise affine) Sigurdur Hafstein Reykjavik University, Iceland 17. July 2013 Workshop on Algorithms for Dynamical Systems and Lyapunov Functions

More information

Flight Dynamics, Simulation, and Control

Flight Dynamics, Simulation, and Control Flight Dynamics, Simulation, and Control For Rigid and Flexible Aircraft Ranjan Vepa CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group, an

More information

CONTROL * ~ SYSTEMS ENGINEERING

CONTROL * ~ SYSTEMS ENGINEERING CONTROL * ~ SYSTEMS ENGINEERING H Fourth Edition NormanS. Nise California State Polytechnic University, Pomona JOHN WILEY& SONS, INC. Contents 1. Introduction 1 1.1 Introduction, 2 1.2 A History of Control

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

5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...

5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References... 1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...

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

/ m U) β - r dr/dt=(n β / C) β+ (N r /C) r [8+8] (c) Effective angle of attack. [4+6+6]

/ m U) β - r dr/dt=(n β / C) β+ (N r /C) r [8+8] (c) Effective angle of attack. [4+6+6] Code No: R05322101 Set No. 1 1. (a) Explain the following terms with examples i. Stability ii. Equilibrium. (b) Comment upon the requirements of stability of a i. Military fighter aircraft ii. Commercial

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

Conjugate Gradient (CG) Method

Conjugate Gradient (CG) Method Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous

More information

New Parametric Affine Modeling and Control for Skid-to-Turn Missiles

New Parametric Affine Modeling and Control for Skid-to-Turn Missiles IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 9, NO. 2, MARCH 2001 335 New Parametric Affine Modeling and Control for Skid-to-Turn Missiles DongKyoung Chwa and Jin Young Choi, Member, IEEE Abstract

More information

8.012 Physics I: Classical Mechanics Fall 2008

8.012 Physics I: Classical Mechanics Fall 2008 MIT OpenCourseWare http://ocw.mit.edu 8.012 Physics I: Classical Mechanics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS INSTITUTE

More information

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A

More information

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background

More information

Weak convergence. Amsterdam, 13 November Leiden University. Limit theorems. Shota Gugushvili. Generalities. Criteria

Weak convergence. Amsterdam, 13 November Leiden University. Limit theorems. Shota Gugushvili. Generalities. Criteria Weak Leiden University Amsterdam, 13 November 2013 Outline 1 2 3 4 5 6 7 Definition Definition Let µ, µ 1, µ 2,... be probability measures on (R, B). It is said that µ n converges weakly to µ, and we then

More information

Computational atmospheric trajectory simulation analysis of spin-stabilised projectiles and small bullets

Computational atmospheric trajectory simulation analysis of spin-stabilised projectiles and small bullets Int. J. Computing Science and Mathematics, Vol. x, No. x, xxxx 1 Computational atmospheric trajectory simulation analysis of spin-stabilised projectiles and small bullets.n. Gkritzapis* Laboratory of Firearms

More information

7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system

7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system 7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically

More information

Robust Stability. Robust stability against time-invariant and time-varying uncertainties. Parameter dependent Lyapunov functions

Robust Stability. Robust stability against time-invariant and time-varying uncertainties. Parameter dependent Lyapunov functions Robust Stability Robust stability against time-invariant and time-varying uncertainties Parameter dependent Lyapunov functions Semi-infinite LMI problems From nominal to robust performance 1/24 Time-Invariant

More information

Susceptibility of F/A-18 Flight Control Laws to the Falling Leaf Mode Part II: Nonlinear Analysis

Susceptibility of F/A-18 Flight Control Laws to the Falling Leaf Mode Part II: Nonlinear Analysis Susceptibility of F/A-18 Flight Control Laws to the Falling Leaf Mode Part II: Nonlinear Analysis Abhijit Chakraborty, Peter Seiler and Gary J. Balas Department of Aerospace Engineering & Mechanics University

More information

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS dx dt DIFFERENTIAL EQUATIONS AND CONTROL PROCESSES N 1, 1997 Electronic Journal, reg. N P23275 at 07.03.97 http://www.neva.ru/journal e-mail: diff@osipenko.stu.neva.ru Control problems in nonlinear systems

More information

DEVELOPMENT of accurate cost-efficient guided projectiles

DEVELOPMENT of accurate cost-efficient guided projectiles JOURNAL OF GUIDANCE,CONTROL, AND DYNAMICS Vol. 33, No. 4, July August 21 Design of a Mortar Projectile with Reciprocating Canards Jonathan Rogers and Mark Costello Georgia Institute of Technology, Atlanta,

More information

Old Dominion University Physics 811 Fall 2008

Old Dominion University Physics 811 Fall 2008 Old Dominion University Physics 811 Fall 008 Projects Structure of Project Reports: 1 Introduction. Briefly summarize the nature of the physical system. Theory. Describe equations selected for the project.

More information

Stability of Nonlinear Systems An Introduction

Stability of Nonlinear Systems An Introduction Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear

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

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR

LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu

More information

On matrix equations X ± A X 2 A = I

On matrix equations X ± A X 2 A = I Linear Algebra and its Applications 326 21 27 44 www.elsevier.com/locate/laa On matrix equations X ± A X 2 A = I I.G. Ivanov,V.I.Hasanov,B.V.Minchev Faculty of Mathematics and Informatics, Shoumen University,

More information

3 rd Tutorial on EG4321/EG7040 Nonlinear Control

3 rd Tutorial on EG4321/EG7040 Nonlinear Control 3 rd Tutorial on EG4321/EG7040 Nonlinear Control Lyapunov Stability Dr Angeliki Lekka 1 1 Control Systems Research Group Department of Engineering, University of Leicester arch 9, 2017 Dr Angeliki Lekka

More information

IBRAHIM EREN. Yildiz Technical University, Mechanical Engineering Department, 34349, Besiktas, Istanbul, Turkey

IBRAHIM EREN. Yildiz Technical University, Mechanical Engineering Department, 34349, Besiktas, Istanbul, Turkey Sādhanā Vol. 33, Part 1, February 008, pp. 45 55. Printed in India Determining large deflections in rectangular combined loaded cantilever beams made of non-linear Ludwick type material by means of different

More information

QUATERNION FEEDBACK ATTITUDE CONTROL DESIGN: A NONLINEAR H APPROACH

QUATERNION FEEDBACK ATTITUDE CONTROL DESIGN: A NONLINEAR H APPROACH Asian Journal of Control, Vol. 5, No. 3, pp. 406-4, September 003 406 Brief Paper QUAERNION FEEDBACK AIUDE CONROL DESIGN: A NONLINEAR H APPROACH Long-Life Show, Jyh-Ching Juang, Ying-Wen Jan, and Chen-zung

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

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

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

Module 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions

Module 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions Cambridge University Engineering Dept. Third Year Module 3F: Systems and Control EXAMPLES PAPER ROOT-LOCUS Solutions. (a) For the system L(s) = (s + a)(s + b) (a, b both real) show that the root-locus

More information

CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version

CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version Norman S. Nise California State Polytechnic University, Pomona John Wiley fir Sons, Inc. Contents PREFACE, vii 1. INTRODUCTION, 1

More information

Computational non-linear structural dynamics and energy-momentum integration schemes

Computational non-linear structural dynamics and energy-momentum integration schemes icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum

More information

Modelling of Opposed Lateral and Longitudinal Tilting Dual-Fan Unmanned Aerial Vehicle

Modelling of Opposed Lateral and Longitudinal Tilting Dual-Fan Unmanned Aerial Vehicle Modelling of Opposed Lateral and Longitudinal Tilting Dual-Fan Unmanned Aerial Vehicle N. Amiri A. Ramirez-Serrano R. Davies Electrical Engineering Department, University of Calgary, Canada (e-mail: namiri@ucalgary.ca).

More information

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) 8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)

More information

Outline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system

Outline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system Outline Control systems Lecture-4 Stability V. Sankaranarayanan Outline Outline 1 Outline Outline 1 2 Concept of Stability Zero State Response: The zero-state response is due to the input only; all the

More information

TRANSONIC AERODYNAMIC AND SCALING ISSUES FOR LATTICE FIN PROJECTILES TESTED IN A BALLISTICS RANGE

TRANSONIC AERODYNAMIC AND SCALING ISSUES FOR LATTICE FIN PROJECTILES TESTED IN A BALLISTICS RANGE EB01 19th International Symposium of Ballistics, 7 11 May 2001, Interlaken, Switzerland TRANSONIC AERODYNAMIC AND SCALING ISSUES FOR LATTICE FIN PROJECTILES TESTED IN A BALLISTICS RANGE Gregg Abate1, Gerald

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

Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015

Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015 Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting

More information

Modeling Realistic Projectile Motion

Modeling Realistic Projectile Motion Modeling Realistic Projectile Motion Matthew Geske 1. Introduction The calculation of simple projectile motion is one of the most basic exercises in physics. Finding, for example, the range of a thrown

More information

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop

Perturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael

More information

Control Systems. Internal Stability - LTI systems. L. Lanari

Control Systems. Internal Stability - LTI systems. L. Lanari Control Systems Internal Stability - LTI systems L. Lanari definitions (AS) - A system S is said to be asymptotically stable if its state zeroinput response converges to the origin for any initial condition

More information

Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, Preliminaries!

Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, Preliminaries! Optimal Control and Estimation MAE 546, Princeton University Robert Stengel, 2018 Copyright 2018 by Robert Stengel. All rights reserved. For educational use only. http://www.princeton.edu/~stengel/mae546.html

More information

Predictive Optimal Pulse jet Control for Symmetric Projectiles

Predictive Optimal Pulse jet Control for Symmetric Projectiles Predictive Optimal Pulse jet Control for Symmetric Projectiles Bradley T. Burchett Rose-Hulman Institute of Technology, Terre Haute, IN, 4783 The linear theory model of a symmetric projectile is well suited

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

Table of Laplacetransform

Table of Laplacetransform Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun

More information

ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS

ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS Letters International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1579 1597 c World Scientific Publishing Company ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS A. S. HEGAZI,H.N.AGIZA

More information

Introduction to Flight Dynamics

Introduction to Flight Dynamics Chapter 1 Introduction to Flight Dynamics Flight dynamics deals principally with the response of aerospace vehicles to perturbations in their flight environments and to control inputs. In order to understand

More information

Robust Control of Robot Manipulator by Model Based Disturbance Attenuation

Robust Control of Robot Manipulator by Model Based Disturbance Attenuation IEEE/ASME Trans. Mechatronics, vol. 8, no. 4, pp. 511-513, Nov./Dec. 2003 obust Control of obot Manipulator by Model Based Disturbance Attenuation Keywords : obot manipulators, MBDA, position control,

More information

Unit 1: Mechanical Equilibrium

Unit 1: Mechanical Equilibrium Unit 1: Mechanical Equilibrium Chapter: Two Mechanical Equilibrium Big Idea / Key Concepts Student Outcomes 2.1: Force 2.2: Mechanical Equilibrium 2.3: Support Force 2.4: Equilibrium for Moving Objects

More information

Vibrations in Mechanical Systems

Vibrations in Mechanical Systems Maurice Roseau Vibrations in Mechanical Systems Analytical Methods and Applications With 112 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Contents Chapter I. Forced Vibrations

More information

Principles of Rocketry

Principles of Rocketry 1-1 Principles of Rocketry 1-2 Water Rockets BASIC CONCEPTS 1-3 What is a Rocket? A chamber enclosing a gas under pressure. A balloon is a simple example of a rocket. Rubber walls compress the air inside.

More information

ME8230 Nonlinear Dynamics

ME8230 Nonlinear Dynamics ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass

More information

Physics 2210 Fall 2011 David Ailion FINAL EXAM. December 14, 2011

Physics 2210 Fall 2011 David Ailion FINAL EXAM. December 14, 2011 Dd Physics 2210 Fall 2011 David Ailion FINAL EXAM December 14, 2011 PLEASE FILL IN THE INFORMATION BELOW: Name (printed): Name (signed): Student ID Number (unid): u Discussion Instructor: Marc Lindley

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

Time-delay feedback control in a delayed dynamical chaos system and its applications

Time-delay feedback control in a delayed dynamical chaos system and its applications Time-delay feedback control in a delayed dynamical chaos system and its applications Ye Zhi-Yong( ), Yang Guang( ), and Deng Cun-Bing( ) School of Mathematics and Physics, Chongqing University of Technology,

More information

(Continued on next page)

(Continued on next page) (Continued on next page) 18.2 Roots of Stability Nyquist Criterion 87 e(s) 1 S(s) = =, r(s) 1 + P (s)c(s) where P (s) represents the plant transfer function, and C(s) the compensator. The closedloop characteristic

More information