A note on stability of motion of a projectile

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

PERFORMANCE ANALYSIS OF ISL S GUIDED SUPERSONIC PROJECTILE. Pierre Wey

Stability and Control

Theoretical and experimental research of supersonic missile ballistics

Controls Problems for Qualifying Exam - Spring 2014

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

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

AB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT

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

Flight Dynamics and Control

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

Chapter 2: Unconstrained Extrema

Aerodynamic jump: A short range view for long rod projectiles

Investigation of the Dynamics and Modeling of a Triangular Quadrotor Configuration

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

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

ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT

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

Modeling & Control of Hybrid Systems Chapter 4 Stability

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

Kinetic Energy of Rolling

MCE693/793: Analysis and Control of Nonlinear Systems

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

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

Dynamical Systems & Lyapunov Stability

Sample Solutions of Assignment 10 for MAT3270B

Vibrations of stretched damped beams under non-ideal boundary conditions

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

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

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

Course syllabus Engineering Mechanics - Dynamics

A Study of Firing Sidewise from an Airplane

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

The Simulation of Wraparound Fins Aerodynamic Characteristics

Designing Information Devices and Systems I Spring 2015 Note 3

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

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

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

Robot Control Basics CS 685

Parametrically Excited Vibration in Rolling Element Bearings

Computation of CPA Lyapunov functions

Flight Dynamics, Simulation, and Control

CONTROL * ~ SYSTEMS ENGINEERING

Robot Dynamics - Rotary Wing UAS: Control of a Quadrotor

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

CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017

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

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

Conjugate Gradient (CG) Method

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

8.012 Physics I: Classical Mechanics Fall 2008

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

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

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

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

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

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

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

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS

DEVELOPMENT of accurate cost-efficient guided projectiles

Old Dominion University Physics 811 Fall 2008

Stability of Nonlinear Systems An Introduction

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

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

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

On matrix equations X ± A X 2 A = I

3 rd Tutorial on EG4321/EG7040 Nonlinear Control

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

QUATERNION FEEDBACK ATTITUDE CONTROL DESIGN: A NONLINEAR H APPROACH

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

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

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

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

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

CONTROL SYSTEMS ENGINEERING Sixth Edition International Student Version

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

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

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

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

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

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

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

Modeling Realistic Projectile Motion

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

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

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

Predictive Optimal Pulse jet Control for Symmetric Projectiles

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

Table of Laplacetransform

ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS

Introduction to Flight Dynamics

Robust Control of Robot Manipulator by Model Based Disturbance Attenuation

Unit 1: Mechanical Equilibrium

Vibrations in Mechanical Systems

Principles of Rocketry

ME8230 Nonlinear Dynamics

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

Translational and Rotational Dynamics!

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

(Continued on next page)

Transcription:

Sādhanā, Vol. 6, Part, August 00, pp. 379 385. Printed in India A note on stability of motion of a projectile SDNAI Institute of Armament Technology, Pune 05, India e-mail: 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

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

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, 3 + 3 > 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 α + β + 3 + ] β + α X, V = β [Ẋ t X t Ẋ + J α β+ β J α β + β X α β + Asymptotic stability is guarranteed for V>0 and α + β + 3 + β + 3 β + α β ], V <0 and is obtained if > 0, 3a β + α β < 0, 3b

38 S D Naik α β + + β 3 β+α >0. 3c α and β can be arbitrary. In particular, α =,β =0 reduce 3 to < 0 a + 3 + > 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, + β + 3 + β β + β > 0, 6a 6b 3 β+ β >0. 6c β + The value of β can be selected in such a way that V>0. Let β = + β + + β + + 3 +. 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 + + 3 [ + ] > 0. 8

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 + + 3 0. 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 + 3 0. 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, 9 + 3 + + 3 0. 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 = 59.06. 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 = 3.7 0 5

38 S D Naik = 0.0003796 3 = 0.5 0 0 = 0.00735 0 5. 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

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: 37 50 Sharma C 977 On the dynamic stability of a missile. Des. Sci. J. 7: 9 3