# Stability derivatives

**Stability derivatives**are a means of linearising the equations of motion of an atmospheric flight vehicle so that conventionalcontrol engineering methods may be applied to assess their stability.The dynamics of atmospheric flight vehicles is potentially very difficult to analyse, because the

forces and moments on the vehicle are seldom simple linear functions of its states. In order to address this problem, and render the analysis of stability and the design of autopilots tractable, it is necessary to deal with linear approximations to the equations of motion. The analysis is then applied to a range of flight conditions.A stability derivative is an incremental change in a force or moment acting on the vehicle corresponding to an incremental change in one of the states.

**Nomenclature**In

aeronautics , the following naming convention is commonly used:The stability axis set differs between aircraft and missiles.Aircraft analysis uses stability axes:

:The x axis points along the velocity vector:The z axis lies in the plane of incidence (i.e. the plane containing the velocity vector and the longitudinal axis), perpendicular to the x axis, such that in level flight, the z axis points downwards.:The y axis completes a right-handed orthogonal set, so that in level flight, it points in the direction of the starboard wing.

Missiles and rockets use principle axes:

:The x axis points along the missile body:The y axis points to starboard:The z axis completes a right-handed orthogonal axis set, i.e. points downwards relative to the xy plane

The following notation is commonly used for the physical quantities involved:::Body Forces: X,Y,Z - for longitudinal, lateral and normal forces.::Moments: L,M,N - for

roll ,pitch and yawing moments respectively::Velocities: u,v,w - for longitudinal, lateral and normal components. :::It is helpful to think of these as projections of the velocity vector on to the three body axes, rather than in terms of the translational motion of the vehicle. As the body rotates relative to thevelocity vector, these components change, even when there is no net change inspeed . In this context the difference between speed and velocity is central to understanding, and not mere pedantry. In general, the body longitudinal axis does not coincide with the velocity vector.::Angular velocities: p,q,r - for roll, pitch and yaw rates**Equations of Motion**The use of stability derivatives is most conveniently demonstrated with missile or rocket configurations, because these exhibit greater symmetry than aeroplanes, and the equations of motion are correspondingly simpler. If it is assumed that the vehicle is roll-controlled, the pitch and yaw motions may be treated in isolation. It is common practice to consider the yaw plane, so that only 2D motion need be considered. Furthermore, it is assumed that thrust equals drag, and the longitudinal equation of motion may be ignored.

The body is oriented at angle $psi$ (psi) with respect to inertial axes. The body is oriented at an angle $eta$ (beta) with respect to the velocity vector, so that the components of velocity in body axes are:::$u=U\; coseta$::$v=U\; sineta$where U is the speed.

The aerodynamic forces are generated with respect to body axes, which is not an inertial frame. In order to calculate the motion, the forces must be referred to inertial axes. This requires the body components of velocity to be resolved through the heading angle $(eta)$ into inertial axes.

Resolving into fixed (inertial) axes:::$u\_f=Ucos(eta)cos(psi)-Usin(eta)sin(psi)=Ucos(eta+psi)$::$v\_f=Usin(eta)cos(psi)+Ucos(eta)sin(psi)=Usin(eta+psi)$The acceleration with respect to inertial axes is found by differentiating these components of velocity with respect to time:::$frac\; \{du\_f\}\{dt\}=frac\; \{dU\}\; \{dt\}\; cos(eta+psi)-Ufrac\; \{d(eta+psi)\}\; \{dt\}\; sin(eta+psi)$::$frac\{dv\_f\}\{dt\}=frac\{dU\}\{dt\}sin(eta+psi)+Ufrac\{d(eta+psi)\}\{dt\}cos(eta+psi)$From

Newton's Second Law , this is equal to the force acting divided by themass . Now forces arise from thepressure distribution over the body, and hence are generated in body axes, and not in inertial axes, so the body forces must be resolved toinertial axes , as Newton's Second Law does not apply in its simplest form to an accelerating frame of reference.Resolving the body forces:::$X\_f=Xcos(psi)-Ysin(psi)$::$Y\_f=Ycos(psi)+Xsin(psi)$Newton's Second Law, assuming constant mass:::$X\_f=mfrac\{du\_f\}\{dt\}$::$Y\_f=mfrac\{dv\_f\}\{dt\}$where "m" is the mass.Equating the inertial values of acceleration and force, and resolving back into body axes, yields the equations of motion:::$X=mfrac\{dU\}\{dt\}cos(eta)-mUfrac\{d(eta+psi)\}\{dt\}sin(eta)$::$Y=mfrac\{dU\}\{dt\}sin(eta)+mUfrac\{d(eta+psi)\}\{dt\}cos(eta)$The sideslip, $eta$, is a small quantity, so the small perturbation equations of motion become:::$X=mfrac\{dU\}\{dt\}$::$Y=mUfrac\{d(eta+psi)\}\{dt\}$The first resembles the usual expression of Newton's Second Law, whilst the second is essentially the

centrifugal acceleration.The equation of motion governing the rotation of the body is derived from the time derivative ofangular momentum :::$N=Cfrac\{d^2psi\}\{dt^2\}$where C is themoment of inertia about the yaw axis.Assuming constant speed, there are only two state variables; $eta$ and $frac\{dpsi\}\{dt\}$, which will be written more compactly as the yaw rate r.There is one force and one moment, which for a given flight condition will each be functions of $eta$, r and their time derivatives. For typical missile configurations the forces and moments depend, in the short term, on $eta$ and r. The forces may be expressed in the form:::$Y=Y\_0\; +\; frac\; \{partial\; Y\}\{partial\; eta\}\; eta\; +frac\; \{partial\; Y\}\{partial\; r\}r$where $Y\_0$ is the force corresponding to the equilibrium condition (usually called thetrim ) whose stability is being investigated.It is common practice to employ a shorthand:::$frac\{partial\; Y\}\{partial\; eta\}=Y\_eta$Thepartial derivative $frac\{partial\; Y\}\{partial\; eta\}$ and all similar terms characterising the increments in forces and moments due to increments in the state variables are called stability derivatives. Typically, $frac\{partial\; Y\}\{partial\; r\}$ is insignificant for missile configurations, so the equations of motion reduce to:::$frac\{deta\}\{dt\}=frac\{Y\_eta\}\{mU\}eta-r$ ::$frac\{dr\}\{dt\}=frac\{N\_eta\}\{C\}eta+frac\{N\_r\}\{C\}r$**tability Derivative Contributions**Each stability derivative is determined by the position, size, shape and orientation of the missile components. In aircraft, the directional stability determines such features as

dihedral of the main planes, size of fin and area oftailplane , but the large number of important stability derivatives involved precludes a detailed discussion within this article. The missile is characterised by only three stability derivatives, and hence provides a useful introduction to the more complex aeroplane dynamics.Consider first $Y\_eta$, a body at an

angle of attack $eta$ generates a lift force in the opposite direction to the motion of the body. For this reason $Y\_eta$ is always negative.At low angles of attack, the lift is generated primarily by the wings, fins and the nose region of the body. The total lift acts at a distance $x\_cp$ ahead of the

centre of gravity (it has a negative value in the figure), this, in missile parlance, is the centre of pressure . If the lift acts ahead of the centre of gravity, the yawing moment will be negative, and will tend to increase the angle of attack, increasing both the lift and the moment further. It follows that the centre of pressure must lie aft of the centre of gravity for static stability. $x\_cp$ is thestatic margin and must be negative forlongitudinal static stability . Alternatively, positive angle of attack must generate positive yawing moment on a statically stable missile, i.e. $N\_eta$ must be positive. It is common practice to design manoeuvrable missiles with near zero static margin (i.e. neutral static stability).The need for positive $N\_eta$ explains why arrows and darts have flights and unguided rockets have fins.

The effect of angular velocity is mainly to decrease the nose lift and increase the tail lift, both of which act in a sense to oppose the rotation. $N\_r$ is therefore always negative. There is a contribution from the wing, but since missiles tend to have small static margins (typically less than a

calibre ), this is usually small. Also the fin contribution is greater than that of the nose, so there is a net force $Y\_r$, but this is usually insignificant compared with $Y\_eta$ and is usually ignored.**Response**Manipulation of the equations of motion yields a second order homogeneous linear differential equation in the angle of attack $eta$::$frac\{d^2eta\}\{dt^2\}-left(frac\{Y\_eta\}\{mU\}+frac\{N\_r\}\{C\}\; ight)frac\{deta\}\{dt\}+left(frac\{N\_eta\}\{C\}+frac\{Y\_eta\}\{mU\}frac\{N\_r\}\{C\}\; ight)eta=0$

The qualitative behavior of this equation is considered in the article on directional stability. Since $Y\_eta$ and $N\_r$ are both negative, the

damping is positive. The stiffness does not only depend on the static stability term $N\_eta$, it also contains a term which effectively determines the angle of attack due to the body rotation. The distance of the center of lift, including this term, ahead of the centre of gravity is called themaneuver margin . It must be negative for stability.This damped oscillation in angle of attack and yaw rate, following a disturbance, is called the 'weathercock' mode, after the tendency of a

weathercock to point into wind.**Comments**We have chosen the state variables to be the angle of attack $eta$ and the yaw rate r, and have omitted the speed perturbation u, together with the associated derivatives e.g. $Y\_u$. This may appear arbitrary. However, we know that the timescale of the speed variation is much greater than that of the variation in angle of attack, so that its effects are negligible as far as the directional stability of the vehicle is concerned. Similarly, the effect of roll on yawing motion is also ignored, because missiles generally have low aspect ratio configurations and the roll inertia is much less than the yaw inertia, consequently the roll loop is expected to be much faster than the yaw response, and is consequently ignored. These simplifications of the problem based on "a priori" knowledge, represent an engineer's approach. Mathematicians prefer to keep the problem as general as possible and only simplify it at the end of the analysis, if at all.

Aircraft dynamics is much more complex than missile dynamics, mainly because the simplifications, such as separation of fast and slow modes, and the similarity between pitch and yaw motions, are not obvious from the equations of motion, and are consequently deferred until a late stage of the analysis. Subsonic transport aircraft have high aspect ratio configurations, so that yaw and roll cannot be treated as decoupled. However, this is merely a matter of degree; the basic ideas needed to understand aircraft dynamics are covered in this simpler analysis of missile motion.

**Control Derivatives**Deflection of control surfaces modifies the pressure distribution over the vehicle, and these are dealt with by including perturbations in forces and moments due to control deflection. The fin deflection is normally denoted $zeta$ (zeta). Including these terms, the equations of motion become:::$frac\{deta\}\{dt\}=frac\{Y\_eta\}\{mU\}eta-r+frac\{Y\_zeta\}\{mU\}zeta$ ::$frac\{dr\}\{dt\}=frac\{N\_eta\}\{C\}eta+frac\{N\_r\}\{C\}r+frac\{N\_zeta\}\{C\}zeta$Including the control derivatives enables the response of the vehicle to be studied, and the equations of motion used to design the autopilot.

**ee also***

Longitudinal static stability

* Neutral point

*Aerodynamic center

*Flight dynamics

*Directional Stability **References*** Babister A W: Aircraft Dynamic Stability and Response. Elsever 1980, ISBN 0-08-024768-7

* Friedland B: Control System Design. McGraw-Hill Book Company 1987. ISBN 0-07-100420-3

*Wikimedia Foundation.
2010.*

### Look at other dictionaries:

**USAF Stability and Control DATCOM**— The United States Air Force Stability and Control DATCOM is a collection, correlation, codification, and recording of best knowledge, opinion, and judgment in the area of aerodynamic stability and control prediction methods. It presents… … Wikipedia**Directional stability**— is stability of a moving body or vehicle about an axis which is perpendicular to its direction of motion. Stability of a vehicle concerns itself with the tendency of a vehicle to return to its original direction in relation to the oncoming medium … Wikipedia**Concern over Derivatives**— ▪ 1996 by Patricia Tehan Derivatives had been acquiring a bad name even before the collapse of the London based merchant bank Barings PLC in February 1995 rocked the world s banking community. But the failure of Barings brought derivatives… … Universalium**Flight dynamics**— is the science of air and space vehicle orientation and control in three dimensions. The three critical flight dynamics parameters are the angles of rotation in three dimensions about the vehicle s center of mass, known as pitch , roll and yaw… … Wikipedia**USAF Digital DATCOM**— The United States Air Force Stability and Control Digital DATCOM is a computer program that implements the methods contained in the USAF Stability and Control DATCOM to calculate the static stability, control and dynamic derivative… … Wikipedia**Linearization**— In mathematics and its applications, linearization refers to finding the linear approximation to a function at a given point. In the study of dynamical systems, linearization is a method for assessing the local stability of an equilibrium point… … Wikipedia**Tangent stiffness matrix**— In computational mechanics, a tangent stiffness matrix is a matrix that describes the stiffness of a system in response to small changes in configuration. It represents tangent in that the energy of the system can be thought of as a high… … Wikipedia**MANGLER**— GERMANY (see also List of Individuals) 26.9.1910 Sangerhausen/D 8.2.1980 Celle/D Kurt Werner Mangler graduated as a physicist from University of Halle/Saale in 1933 and in 1934 there submitted a PhD Thesis in mathematics. He stayed from 1936 to… … Hydraulicians in Europe 1800-2000**fin**— ► NOUN 1) an external organ on the body of a fish or other aquatic animal, used for propelling, steering, and balancing. 2) an underwater swimmer s flipper. 3) a projection on an aircraft, rocket, or car, for providing aerodynamic stability.… … English terms dictionary**finned**— fin ► NOUN 1) an external organ on the body of a fish or other aquatic animal, used for propelling, steering, and balancing. 2) an underwater swimmer s flipper. 3) a projection on an aircraft, rocket, or car, for providing aerodynamic stability.… … English terms dictionary