0% found this document useful (0 votes)
131 views134 pages

Nasa SP 3070

This document provides equations for transferring aerodynamic force and moment data between different reference frames. It includes equations for transferring accelerometer, rate gyro, and angle of attack instrumentation data between vehicle and instrument axes. It also provides transformations to transfer aerodynamic derivatives from one reference center to another, or to a rotated axes system. The equations allow analysis of flight data in various reference frames.

Uploaded by

panbu
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
131 views134 pages

Nasa SP 3070

This document provides equations for transferring aerodynamic force and moment data between different reference frames. It includes equations for transferring accelerometer, rate gyro, and angle of attack instrumentation data between vehicle and instrument axes. It also provides transformations to transfer aerodynamic derivatives from one reference center to another, or to a rotated axes system. The equations allow analysis of flight data in various reference frames.

Uploaded by

panbu
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 134

GAINER and HOFFMAN

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION


f
NASA SP-3070

Thomas G. Gainer and Sherwood Hoffman


Langley Research Center
Hampton, Virginia

Prepared by Langley Research Center

Scientific and Technical In[ormation O[fice 1972

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION


Washington, D.C.
For sale by the National Technical Information Service, Springfield, Virginia 22151 -- Price $3.00
PREFACE

A summary of equations often used in free-flight and wind-tunnel data reduction and
analysis is presented. Included are transfer equations for accelerometer, rate-gyro, and
angle-of-attack instrumentation; axes-system transfers of aerodynamic derivatives; and
methods for measuring moments.of inertia. In general, the equations are in a complete
form; for example, those terms are retained that are missing when planar symmetry is
assumed for airplanes.

iii
TABLE OF CONTENTS

Page
iii
PREFACE .......................................
1
INTRODUCTION ....................................
2
SYMBOLS .......................................

SECTION I - EQUATIONS INVOLVING BASIC FLIGHT MEASUREMENTS ..... 17


GENERAL AXES TRANSFORMATIONS FOR COMPONENTS OF
ACCELERATION, LINEAR VELOCITY, AND ANGULAR VELOCITY ..... 17
Transfer From Vehicle Reference Axes to Instrument Axes .......... 18
Transfer From Instrument Axes to Vehicle Reference Axes .......... 19
CORRECTIONS TO ANGLES OF ATTACK AND SIDESLIP ............ 21
21
Angle of Attack ..................................
Angle of Sideslip ................................ 21
DETERMINATION OF FREE-STREAM MACH NUMBER . . . ......... 21
DETERMINATION OF AERODYNAMIC FORCES AND MOMENTS FROM
ACCELEROMETER AND RATE-GYRO READINGS ............. 22
DIRECT FLIGHT MEASUREMENTS OF LIFT AND DRAG ............ 23

SECTION II - TRANSFER OF AERODYNAMIC FORCE AND MOMENT


COEFFICIENTS AND DERIVATIVES TO ANOTHER REFERENCE
25
CENTER ....................................
25
TRANSFORMATIONS FOR a,/3,V DERIVATIVES ................
X-Axis Force Coefficients and Derivatives ............... . . . 25
Y-Axis Force Coefficients and Derivatives .................. 26
Z-Axis Force Coefficients and Derivatives .............. . . . . 27
X-AXis Moment (Roll) Coefficients and Derivatives .............. 27
Y-Axis Moment (Pitch) Coefficients and Derivatives ............. 28
Z-AXis Moment (Yaw) Coefficients and Derivatives .............. 29
TRANSFORMATIONS FOR u,v,w DERIVATIVES ......... ....... 30
30
X-AXis Force Coefficients and Derivatives .................
Y-AXis Force Coefficients and Derivatives .................. 31
Z-Axis Force Coefficients and Derivatives .................. 31
X-Axis Moment (Roll) Coefficients and Derivatives .............. 32
Y-Axis Moment (Pitch) Coefficients and Derivatives ............. 32
Z-Axis Moment (Yaw)Coefficients and Derivatives .............. 33

v
Page
SECTION Ill - TRANSFER OF AERODYNAMIC FORCE AND MOMENT
COEFFICIENTS AND DERIVATIVES TO A ROTATED AXES SYSTEM 34
DESCRIPTIONS OF AXES SYSTEMS
........................ 34
Body Axes
• " " • • " " • • • • • • • • • • • • • • • • • • • • • • • . . . . . . _. 34
Principal Axes
• " " • " • • • • • • • • • • • • • • • • • • • • • • • • . . . . . . 34
Flight Stability Axes • " " " " • " " • " " • " " " " " " " • " ." " • • • • • • • • • 35

Wind-Tunnel Stability Axes • • • • • " • • • • • • • • • • • • • ' " • • • • • * • 35

Wind Axes
• • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • 35

NOTES ON USE OF TRANSFORMATION EQUATIONS ............... 35


DIRECT TRANSFORMATIONS
............................ 37
Static Force and Moment Coefficients (Direct, Table III) ............ 37
Static-Stability Derivatives (Direct, Table III) 37
Dynamic-Stability Derivatives (Direct, Table III) ................ 39
u,v,w Derivatives (Direct, Table HI) 41
INVERSE TRANSFORMATIONS ........................... 46
Static Force and Moment Coefficients (Inverse, Table IV) ............ 46
Static-Stability Derivatives (Inverse, TableIV) . . . . . . . . . . . . . . . . . 46
Dynamic-Stability Derivatives (Inverse, Table IV) ................ 48
u,v,w Derivatives (Inverse, Table IV) ..........• • • • • • • . . . . . 50

SECTION IV - TRANSFORMATION EQUATIONS FOR MOMENTS


OF INERTIA . ...........
....................... 55
DIRECT TRANSFORMATIONS (TABLE III) 55
INVERSE TRANSFORMATIONS (TABLE IV) . ......... • • • • • • . . . . 56

SECTION V - EQUATIONS OF MOTION FOR SIX DEGREES OF FREEDOM .... 57


FORCE EQUATIONS ........... ...................... 57
X-Axis Forces .................. ................. 57
Y-Axis Forces ....
• • ° ° " • ° ° " • • • • • • • " • ° • • • * • • • • • • • • 58

Z-Axis Forces ........... • • • • • • • • • • • • • • • • • • • • • . . . 59


MOMENT EQUATIONS ................................. 59
X-Axis Moment (Roll) ............................... 59
Y-Axis Moment (Pitch) .............................. 60
Z-Axis Moment (Yaw) ............................... 61
AUXILIARY EQUATIONS . ........ ..................... 61
Relationship Between Euler Angles and Angular Velocities ........... 61
Vehicle Coordinates ................................ 62
Trajectory Parameters (See Fig. 6) ....................... 63
Angle of Attack, Sideslip, and Relative Velocity (See Fig. 7) ........... 63

vi
Page

Wind Corrections (See Fig. 8) .......................... 64

Resolution of Engine Thrust and Torque Into Components Along


Vehicle Reference Axes (See Fig. 9) ..................... 64

Components of Gravitational Acceleration Along X,Y,Z Vehicle


Reference Axes With Earth-Oblateness Effects Included ........... 65

APPENDIX A - SUMMARY OF FREQUENTLY USED FORMS OF AXES


TRANSFORMATIONS AND EQUATIONS OF MOTION .............. 67
EULER ANGLE TRANSFORMATION BETWEEN TWO ORTHOGONAL
AXES SYSTEMS ................................ 67
Direct Transformation .............................. 67
Inverse Transformation ............................ 67
TRANSFORMATIONS FOR ACCELEROMETER AND RATE-GYRO
MEASUREMENTS 68
Case I - Orthogonal Instrument Axes; No Restrictions on
Misalinement Angles ........................ 68

Case II - Nonorthogonal Instrument Axes; Small Misalinement Angles . . . 69


SPECIAL FORMS OF TRANSFORMATIONS FOR AERODYNAMIC FORCE
AND MOMENT COEFFICIENTS AND STABILITY DERIVATIVES ...... 70

Simplified Forms for Transferring Coefficients and Derivatives to


Another Reference Center ........................... 70
Transfer of Coefficients and Derivatives From Body to Wind-Tunnel

Stability Axes ................................. 72

Transfer of Coefficients and Derivatives From Wind-Tunnel


Stability to Body Axes ............................. 73
TRANSFORMATION EQUATIONS FOR MOMENTS AND PRODUCTS
OF INERTIA .................................. 75

Body to Flight Stability Axes ............................ 75


Body to Principal Axes .............................. 75
Flight Stability to Body Axes ........................... 75
Flight Stability to Principal Axes ........................ 75
Principal to Body Axes .............................. 76
Principal to Flight Stability Axes ........................ 76
SPECIAL FORMS OF EQUATIONS OF MOTION .................. 76
Coupled Linear Equations of Motion ....................... 76
Uncoupled Equations of Motion .......................... 78
Wind-Axes Equations for a Point Mass ..................... 79

vii
Page
APPENDIX B - DERIVATION OF EQUATIONS PRESENTED IN SECTION IT
FOR TRANSFER OF COEFFICIENTS TO NEW REFERENCE CENTER 8O

APPENDIX C - METHODS OF MEASURING CENTER-OF-GRAVITY


LOCATIONS AND MOMENTS OF INERTIA OF MODELS AND FLIGHT
VEHICLES .................................... 84
CENTER-OF-GRAVITY LOCATION ........................ 84
Longitudinal c.g. Location • • • 84
Vertical c.g. Location ................. • • • • • • • . . . . . . 84
MOMENTS OF INERTIA .......... . . . . 84
Compound-Pendulum Method (Fig. ii) .............. . . . . . . . 84
Sp ring Method .................................. 85
Spring Method for Full-Scale Vehicles ..................... 85
Torsion-Pendulum Method ............................ 86
Multifilar-Pendulum Methods ........................... 86

APPENDIX D - DETERMINATION OF LONGITUDINAL AND LATERAL

STABILITY DERIVATIVES BY USING SIMPLIFIED LINEAR ANALYSIS .... 87

APPENDIX E - USE OF DIRECTION COSINES AND QUATERNIONS IN


MOTION CALCULATIONS ............................ 90
DIRECTION-COSINE METHOD . • • • " ° " " " " • • • • • • • • • • • • • • • , 90

QUATERNION METHOD .............................. 92


Basic Quaternion Relationships ......................... 92
Relationship Between Euler Parameters and Euler Angles ........... 93

APPENDIXF - SCALING PARAMETERS .......... • • • • • • • . . . . . 94


DYNAMIC-SCALING PARAMETERS ........................ 94
AERODYNAMIC-SCALING PARAMETERS 94
AEROE LASTIC-SCALING PARAMETERS .................... 96
COMMENTS ON SCALING PROCEDURES ......• • • • • • • • • • . . . . . 96
Dynamic Scaling ................................. 96
Aerodynamic Scaling ............................... 96

REFERENCES .................................... 98

BIBLIOGRAPHY ON METHODS OF ANALYZING FLIGHT DATA .......... 99

TABLES ........................................ 102


I - CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS ......... 102
H- RELATIONSHIPS BETWEEN ol,_,V AND u,v,w DERIVATIVES ..... 103
III- ANGLE DESIGNATIONS FOR DIRECT TRANSFORMATIONS ....... 104

viii
Page
IV - ANGLE DESIGNATIONS FOR INVERSE TRANSFORMATIONS ........ 105
V - DESIGNATIONS OF FORCE AND MOMENT COEFFICIENTS FOR
DIFFERENT AXES SYSTEMS ........................ 106

VI - SCALE FACTORS FOR DYNAMIC SCALING ................. IQ7

FIGURES ........................................ 108

1 - Systems of vehicle reference axes and instrument axes ........... 108


2 - Axes systems for transfer from vehicle c.g. to new reference center
by equations of section II .......................... 109
3 - Systems of vehicle reference axes, including body, principal, wind,
flight stability, and wind-tunnel stability .................. 109
4 - Relationship between earth-centered inertial axes, gravity axes, and
vehicle reference axes ........................... 110

5 - Systems of gravity axes and vehicle reference axes ............. 111


6 - Relation of heading angle H and flight-path angle 7 to earth-
centered inertial axes ............................ 112

7 - Resolution of relative velocity into components along vehicle


body axes .................................. 113
8 - Directions of surface and geostrophic winds ............. .... 114
9 - Alinement with respect to vehicle reference axes of thrust and torque
due to rotating mass of engine ....................... 115
10 - Determinatiori of vertical location of center of gravity ............ 115
11 - Measurement of moment of inertia by compound-pendulum method ..... 116
12 - Measurement of moment of inertia by spring method ............ 116
13 - Determination of vertical location of center of gravity and rolling
moments of inertia for full-scale airplanes. (Reproduced
from ref. 3.) ................................. 117
14 - Methods of measuring yawing moments of inertia .............. 118
(a) Torsion pendulum ............................ 118
(b) Bifilar pendulum ............................. 118
15 - Damped angle-of-attack oscillation assumed in analysis of appendix D . . . 119

ix
INTRODUCTION

The equations in this report are the coordinate transformation and motion equations
used in the various tasks associated with free-flight and wind-tunnel data reduction and
analysis. These tasks range from reducing flight data to calculating the motions on a '
digital or analog computer and to applying various techniques for analyzing the data, such
as in references 1 and 2.

While many publications contain a number of these equations, no one contains all
that are usually needed in a complicated aerodynamic analysis; even the more nearly com-
plete reports (refs. 3 and 4, for example) omit the equations for transferring aerodynamic
stability derivatives from one moment reference to another. Moreover, in most cases
the equations are simplified, when they are presented, by assumptions such as small
angles of attack, zero sideslip, and small perturbation motions. Expanded forms of many
of the equations, on the other hand, are needed in special problems that may arise. For
example, parawing vehicles, which have their center of gravity located well below the wing
surface, require the expanded forms of the axes transformations when data measured about

a point on the wing are to be transferred to the center of gravity; reentry motion studies
sometimes involve large-amplitude motions so the complete forms of the transformations,
without the assumptions of small angles of attack or sideslip, are needed. The engineer

working on any of these special problems usually has to derive these equations himself,
and this can be time consuming.

The purpose of this report is to provide the basic equations from which many of the
equations needed in a particular analysis can be generated. A comprehensive summary
of the basic axes transformation and motion equations is included, with most of these
given in their expanded, most general forms. Once these expanded forms are available,
the simpler forms can be written out fairly easily, and yet the general forms are here
when needed for special cases.

The general forms presented include axes transformations that enable transfer back

and forth between any of the five axes systems that are encountered in aerodynamic analy-
sis. Equations of motion are presented that enable calculation of motions anywhere in the
vicinity of the earth. Special problems are also considered; since flight instruments,
such as accelerometers or rate gyros, are not always alined along mutually perpendicular
axes, the procedure for correcting instrument readings for nonorthogonal alinements is
outlined.

In addition to these general forms, many of the simplified forms used frequently in
practical applications are summarized in appendix A.

Other relationships are presented in appendixes B to F. For example, appendix C


summarizes methods used to measure moments of inertia of models and full-scale
vehicles; appendix E discusses the use of the direction-cosine and the quaternion methods,
often used in place of Euler angles in specifying vehicle alinement; appendix F discusses
the scaling parameters used in model testing. However, throughout this paper, the empha-
sis is on providing the basic equations. For discussions of their development and of the
procedures used in their application, the user should turn to general published works on
flight-motion analysis. A comprehensive bibliography of these works is provided _and
includes textbooks and reports dealing with stability, control, and performance as well
as reports discussing various techniques for extracting stability derivatives from flight
data.

SYMBOLS

Throughout this paper, symbols are defined in terms of SI Units with equivalent U.S.
Customary Units given parenthetically. Factors for converting from U.S. Customary to
SI Units are given in table I.

eneralized angle of attack, defined for various axes systems in tables III
A and IV, rad (deg)
ross-sectional area in eq. (F-11) of appendix F, m2 (ft 2)

A _ angle between surface wind vector and plane of local horizontal, measured
perpendicular to plane of local horizontal (fig. 8), rad (deg)

AD acceleration along flight path, g units (lg = 9.807 m/sec 2)

AL acceleration in lift direction, g units (lg = 9.807 m/sec 2)

components of acceleration along X,Y,Z vehicle reference axes


Ax,cg,Ay,cg,Az,cg
at c.g., respectively, g units (lg = 9.807 m/sec2)

components of acceleration indicated by accelerometers along


AX ,i,Ay ,i,A Z ,i
Xi,Yi,Z i instrument axes, respectively, g units
(lg = 9.807 m/sec2)

a speed of sound, m/sec (ft/sec)

a T damping constant defined by eq. (D-4) of appendix D, dimensionless

diameter of circle around which wires or rods are attached in bifilar or

trifilar methods of measuring moments of inertia in eq. (C-9) of


appendix C (see, also, fig. 14(b)), m (in.)
a e
equatorial radius of earth's reference ellipsoid, m (ft)

distance from knife edge to spring (fig. 13), m (ft)


a s

B generalized angle of sideslip, defined for various axes systems in tables III'
and IV, rad (deg)

S t angle between projection of wind vector onto plane of local horizontal and
component of wind velocity tangent to circle of local horizontal (fig. 8),
rad (deg)

b span, m (ft)

C constant in eq. (D-3) of appendix D, rad (deg)

CA axial-force coefficient

crosswind coefficient
CC

CD drag coefficient

CD drag coeff{cient for wind-tunnel stability axes

CL lift coefficient

lift coefficient in steady-state level flight


CL,o

eneralized aerodynamic moment coefficients about X,Y,Z vehicle


.... Moment
reierence axes, respectively,
C/,C m ,Cn
qooD £

kbody-axes moment coefficients

Cm,o
pitching-moment coefficient in steady-state levei flight

CN normal-force coefficient

X,Y,Z vehicle refer-


eneralized aerodynamic force coefficients about
.. . Force
Cx,Cy,Cz
ence axes, respectively, _S
ody-axes force coefficients

3
reference chord, m (ft)

Cp specific heat at constant pressure, J/kg-°C (Btu(thermochemical)/lbm-OF)

Cv specific heat at constant volume, J/kg-OC (Btu(thermochemical)/lbmrOF)

D
erodynamic drag in eqs. (I-25) of Sec. I, N (lbf)
irection-cosine matrix defined in eq. (E-2)of appendix E

dx,dy,d z displacements of centroid of jet interface, measured with respect to X,Y,Z


vehicle reference axes, respectively, m (ft)

E modulus of elasticity, N/m 2 (lbf/ft 2)

e0,e 1 ,e 2 ,e 3 Euler parameters in eq. (E-9) of appendix E

Fr,x,Fr,y,Fr,z forces due to jet control (not main rockets) along X,Y,Z vehicle
reference axes, respectively, N (lbf)

Fx,Fy,F z forces along X,Y,Z vehicle reference axes, respectively, N (lbf)

G transformation matrix composed of Euler parameters

GE geocentric gravitational constant, m3/sec 2 (ft3/sec2)

g acceleration due to gravity, m/sec 2 (ft/sec 2)

gx'gY 'gz components of gravitational acceleration along Xg,Yg,Zg gravity axes,


respectively, for an oblate earth, m/sec 2 (ft/sec2)

H angular direction vehicle is traveling (fig. 6); angle measured positive


clockwise from north, i.e., for H = 90 °, vehicle is traveling toward east,
rad (deg)

h altitude, m (ft)

I moment of inertia, kg-m 2 (slug-ft 2)

I' - Iy
q_oS6

4
Ic
mounting-cradle moment of inertia about vehicle roll axis, kg-m 2
(slug-ft2)

Irm moment of inertia of engine rotating mass, kg-m 2 (slug-ft2)

Ix,Iy,I z moments of inertia about X,Y,Z vehicle reference axes, respectively,


kg-m 2 (slug-ft2)

,!xy,Ixz ,Iy z products of inertia, Ixy= _ xy dm, IXZ = _ xz dm,

Iyz = _ yz dm, kg-m 2 (slug-ft2)

i,j ,k unit vectors having properties given in eqs. (E-10) of appendix E

J2,J3,J4 second, third, fourth zonal harmonics, respectively, of earth's reference


ellipsoid (values given following eq. (V-45) of Sec. V)

j unit vector (see i,j,k)

K spring constant, N/m (lbf/in.)

K1,K2,...,K6 constants defined in eqs. (I-13) of Sec. I

Ko boom misalinement angle at zero g, measured with respect to X vehicle

reference axis; positive with boom inclined above X-axis, rad (deg)

Ks constant for springs shown in fig. 13, N/m (lbf/ft)

values of Ks for springs of different strengths, N/m (lbf/ft)


KS ,I 'Ks ,2

Kt,1 'Kt,2 constants, determined from Ks, 1 and Ks,2, used in measuring moments of
inertia in eqs. (C-6) and (C-7) of appendix C, N/m (lbf/ft)

unit vector (see i,j,k)

torsion-spring constant in eq. (C-8) of appendix C, N-m/rad (ft-lbf/rad)


k
coefficient of heat conduction in eq. (F-9) of appendix F, J/m-sec-K

(Btu (the rmochemical) -in./ft 2 -sec -OF)


L
eocentric
erodynamiclatitude
lift in ofeqs.
vehicle,
(I-25) positive
of Sec. I,in Nnorthern
(ibf)
hemisphere, in Sec. V
and appendix E, rad (deg)

L,M,N aerodynamic moments about X,Y,Z vehicle reference axes, respectively,


N-m (ft-lbf)

Lrm,Mrm,Nrm moments due to gyroscopic action of engine rotating mass about


X,Y,Z vehicle reference axes, respectively, N-m (ft-lbf)

(ft)
ength used in nondimensionalizing moments, m
haracteristic length, m (ft)

T length of wire or rod on which model is suspended in measuring moment of


inertia by multifilar-pendulum method in appendix C, m (ft)

m (ft)
istance between knifefore edgeand toaftsupporting
from weighing scales
spring in eq.
in eq. (C-1)
(C-5) appendix
of appendix C,
C,
m (ft)

_XX,_Xy,...,_ZZ direction cosines defining orientation of rotated (primed) axes

system with respect to initial (unprimed) axes system; _XX is


cosine of angle between X-axes of two systems; _XY, cosine of
angle between X-axis of rotated system and Y-axis of initial sys-

tem; . . .; and _ZZ, cosine of angle between Z-axes of two sys-


tems (eqs. (E-3) of appendix E)

M Mach number

M aerodynamic moment (see L,M,N)

Mrm moment due toengine rotating mass (see Lrm,Mrm,Nrm )

Mr,x,Mr,y,Mr,z moments due to jet controls (not main rockets) about


X,Y,Z vehicle reference axes, respectively, N-m (ft-lbf)

Mx,My,M z moments about X,Y,Z vehicle reference axes, respectively,


N-m (ft-lbf)
m instantaneous mass of vehicle, kg (slugs)

mVoo
m T

q_S

rh time rate of change of vehicle mass (such as due to fuel consumption);


negative when vehicle losing mass, kg/sec (slugs/sec)

mass of weighing cradle used in measuring moments of inertia, kg (slugs)


m e

N aerodynamic moment (see _,M,N)

Nrm moment due to engine rotating mass (see Lrm,Mrm,Nrm)

P period of oscillation, sec

P1 period of oscillation with spring constant Kt,1, sec

P2 period of oscillation with spring constant Kt,2, sec

PM period of oscillation for model alone, sec

PM+S period of oscillation for model plus supporting hardware, sec

PS period of oscillation for supporting hardware, sec

Po,Qo,Ro initial, steady-state values of angular-velocity components along


X,Y,Z vehicle reference axes, respectively, rad/sec

static pressure, N/m 2 (lbf/ft2)

p ,q,r components of angular velocity about X,Y,Z vehicle reference axes,


respectively, rad/sec (deg/sec)

Pi,qi,ri components of angular velocity indicated by rate gyros about

Xi,Yi,Zi instrument axes, respectively, rad/sec

qo angular-velocity component (see Po,Qo,Ro)

7
q
ngular-velocity component (see p,q,r)
uaternion in eq. (E-9) of appendix E

qc impact pressure, N/m 2 (lbf/ft2)

qi angular-velocity component (see pi,qi,ri)

free-stream dynamic pressure, N/m 2 (lbf/ft 2)

conjugate of quaternion in eq. (E-12) of appendix E

distance of c.g. of vehicle from center of earth (fig. 4), m (ft)


R
l/scale factors, with subscripts, as defined in table VI and discussed
__ in appendix F

R e radius of assumed spherical earth, 6378.123 km (20 925 631 ft)

Ro angular-velocity component (see Po,Qo,Ro)

R 1 vertical reaction force (fig. 10) in eq. (C-2) of appendix C, N (lbf)

R 1 ,R2 ,R 3 forces at supports in setups for determining c.g. location in eq. (C-l) of
appendix C, N (lbf)

r angular-velocity component (see p,q,r)

vector expressing distance from c.g. to origin of instrument-axes system,


m (ft)

rg distance from origin of gravity-axes system to vehicle reference center,


m (ft)

ri angular-velocity component (see pi,qi,ri)

vehicle reference area, m 2 (ft2)

eqo

T N (lbf)
Iotal vehicle thrust in eq. (V-36) of Sec. V and in (A-13d) of appendix A,
emperature in eq. (F-9) of appendix F, K (OF)
TD component of thrust along flight path, N (lbf)

To
lrnitial thrust in eq. (V-36) of Sec. V
epresentative temperature of flow in eq. (F-10) of appendix F, K (OF)

Tx,Ty,T z components of thrust along X,Y,Z vehicle reference axes, respectively,


N (lbf)

t time, sea

time to damp to one-half amplitude, sec


tl/2

t1 arbitrary time defined in fig. 15, sec

t2 arbitrary time defined in fig. 15, sec

U o ,V o ,W o initial, steady-state values of linear-velocity components along


X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec)

U,V,W components of vehicle absolute (inertial) velocity along X,Y,Z vehicle


reference axes, respectively, m/sec (ft/sec)

U,V,W components of relative velocity (velocity of vehicle with respect to air),


along X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec)

ui,vi,w i indicated components of relative velocity along Xi,Yi,Zi instrument axes,


respectively, m/sec (ft/sec)

V resultant velocity, m/see (ft/sec)

AV increment in resultant free-stream velocity, m/see (if/see)

Vo velocity component (see Uo,Vo,Wo)

Voo free-stream reference velocity, m/sec (ft/sec)

Vx,Vy,Vz components of geostrophic-wind velocity, due to earth rotation


(atmosphere rotates with earth), relative to X,Y,Z vehicle reference
axes, respectively, m/sec (ft/sec)

V inertial-velocity component (see U ,V ,W)

relative-velocity component (see tiff,W)


relative-velocity component (see ui,vi,wi)

W vehicle weight, N (lbf)

W resultant velocity of earth surface wind, m/sec (ft/sec) :_

We weight of weighing cradle, N (lbf)

WM weight of model alone, N (lbf)

WM+S weight of model plus supporting hardware, N (lbf)

Ws weight of model supporting hardware, N (lbf)

Wo velocity component (see Uo,Vo,Wo)

Wx,Wy,W z components of wind velocity, due to surface winds, relative to


X,Y,Z vehicle reference axes, respectively, m/sec (ft/sec)

W inertial-velocity component Csee u,v,w)

W relat'xve-velocity component (see t_,_,_)

wi relative-velocity component(see ui,vi,wi)

X,Y,Z vehicle reference axes

X,Y,Z aerodynamic forces along X,Y,Z vehicle reference axes, respectively,


N Clbf)

Xe ,Ye, Ze right-handed inertial axes with origin at center of earth (fig. 4)

Xg,Yg,Zg gravity axes with origin at surface of earth (fig. 4)

Xi ,Yi, Z i axes in orthogonal coordinate system with origin at an instrument at point


x,y,z and with system alined at angles _,O,dp with X,Y,Z vehicle
reference axes at c.g., respectively, as shown in fig. 1

x,y,z distances measured from X,Y,Z vehicle reference axes, respectively, m (ft}

10
distances from c.g. to another point on vehicle, m (it)
x,y ,Z

displacements of vehicle with respect to Xg,Yg,Zg gravity axes, respec-


Xg,yg,Zg
tively (figs. 4 and 5), m (It)

coordinates of X-axis accelerometer or rate gyro, measured from


£X,§X,zX
X,Y,Z vehicle reference axes at c.g., respectively, m (it)

coordinates of Y-axis accelerometer or rate gyro, measured from


Xy,yy,Zy
X,Y,Z vehicle reference axes at c.g., respectively, m (it)

coordinates of Z-axis accelerometer or rate gyro, measured from


_Z,§ZfiZ '
X,Y,Z vehicle reference axes at c.g., respectively, m (it)

distance from knife edge to vertical reaction force R1 (fig. I0), m


Yl

distance from knife edge to vehicle c.g. (fig. 13), m (It)


ga

vertical distance to c.g. of weighing cradle (figs. i0 and 13), m (it)


gc

distance from knife edge to c.g. of model, m (it)


zM

distance from knife edge to c.g. of model plus supporting hardware, m


zM+S

distance from knife edge to c.g. of supporting hardware, m (it)


zS

angle of attack, rad (deg)


OL
'_+an-1 _b = tan-1 Wb - VZ _b -
Ub Ub Vx,b

AOl position error, error due to location of angle-of-attack sensor in flow of


vehicle, equal to free-stream angle of attack minus true or corrected local
angle of attack; positive for a free-stream angle of attack greater than
local angle of attack, rad

difference between angle of attack and trim angle of attack at time tl, rad

difference between angle of attack and trim angle of attack at time t2, rad

11
S a
increment in angle of attack caused by boom bending under static air load,
measured with respect to boom center line; positive for upward deflection
of boom, rad

sb change in angle of attack caused by bending of mounting boom under inertial


load, measured with respect to boom center line; positive for upward
deflection of boom, rad

s d change in angle of attack caused by vehicle rotation, measured with respect


to vehicle reference axis, rad

si angle of attack indicated by vane or other sensor, measured with respect to


center line of mounting boom; positive with nose of vane pointed down, rad

s 0 angle between X body axis and X flight stability axis (fig. 3), rad (deg)

s t trim angle of attack, rad (deg)

S u
change in angle of attack caused by upwash from mounting boom, measured
with respect to free-stream velocity; positive for upwash (upward flow
normal to free stream), rad

S V vane floating angle, angle caused by slight mass unbalance or by warp in


vane surface due to imperfections in manufacture, measured with respect
to local velocity vector; positive for nose-up deflections of vane, tad

f
.ngle
in-1 v.._
ofb sideslip,
= sin -1 Vbrad- Vy_
(deg)
b - Wy_ b
V V

sideslip position error, rad

fa
increment in sideslip angle caused by boom bending under static air load,
measured with respect to boom center line; positive for boom deflected
to right, rad

fib change in sideslip angle caused by bending of mounting boom under inertial
loads, measured with respect to boom center line; positive for boom
deflections to right, tad

12
change in sideslip angle caused by vehicle rotation, measured with respect
_d
to vehicle reference axis, rad

sideslip angle indicated by vane or other sensor, measured with respect to


center line of mounting boom; positive with nose of vane pointed to right
of boom, tad

_u
change in sideslip angle caused by sidewash from mounting boom, measured
with respect to free-stream velocity; positive for positive sidewash, rad

sideslip-vane floating angle, angle caused by warp in sideslip vane due to


_V
imperfections in manufacture, measured with respect to local velocity
vector; positive for vane deflected to right, rad

F transformation matrix for orthogonal axes system, defined in eq. (A-2) of

appendix A

transformation matrix for nonorthogonal axes system, defined in


Fnonorthogonal
eq. (I-14) of Sec. I

atio of specific heats in eqs. (I-19) and (I-20) of Sec. I, Cp/C v


light -p at'h angle in eq. (V-17) of Sec. V (fig. 65, rad (deg5

control deflection, rad (deg5

control deflections (aileron, elevator, rudder, respectively) rad (deg)


5a,Se,5r

change in engine rpm


5RPM

angle between Xp principal and X body axes (fig. 35, rad (deg5

strain (elongation per unit length5 measured on model, cm/cm (in./in.)


Em

strain (elongation per unit length) measured on prototype, cm/cm (in./in.5

angle between Xs flight stability and Xwt wind-tunnel stability axes


(fig. 3), rad (deg5

angle between Xp principal and Xs flight stability axes (fig. 3), rad (deg 5

13
X geocentric longitude of vehicle, positive counterclockwise looking in direction
of positive Ze inertial axis (fig. 4), rad

coefficient of viscosity, N-sec/m 2 (slugs/ft-sec)

/2
kinematic viscosity in eq. (F-9) of appendix F, la/p, m2/sec (ft2/sec)

phase angle, rad (deg)

P atmospheric density, kg/m 3 (slugs/ft 3)

(7
ange angle in eq. (V-21) of Sec. V, rad (deg)
urface tension in eq. (F-7) of appendix F, N/m (lbf/ft)

roll angle defined in fig, 10, tad (deg)

Euler angles defining angular alinement of one axes system with respect
to another axes system, rad (deg)

_Pg,Og,_bg Euler angles defining alinement of X,Y,Z vehicle reference axes, respec-
tively, with respect to gravity-axes system (fig. 5), tad (deg)

Euler angles describing alinement of engine-thrust axes (fig. 9), rad (deg)

Euler angles of X-axis accelerometer or rate gyro, measured with respect


to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg)

Euler angles of Y-axis accelerometer or rate gyro, measured with respect


to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg)

_PZ,6Z,q_ z Euler angles of Z-axis accelerometer or rate gyro, measured with respect
to X,Y,Z vehicle reference axes at c.g., respectively, rad (deg)

_e rate of rotation of earth, rad/sec

_rm angular velocity of engine rotating mass, rad/sec

¢o
resultant angular velocity in Sec. I, rad/sec
requency of oscillation in appendix D, rad/sec

14
_X' '_Y' '°_Z' angular-velocity components along X',Y',Z' vehicle reference axes,
respectively, usually measured with respect to gravity-axes or
inertial-axes system, rad/se c

Subs cripts:

b body axes

0 initial conditions

P principal axes

S flight stability axes

w wind axes

wt wind-tunnel stability axes

Notation:

• (Dot) first derivative with respect to time

' (Prime) unless otherwise specified, a primed quantity is one referred to axes sys-
tem located at point x,y,z and/or alined at angles _,0,_ with respect
to initial reference axes system

& perturbation quantity (unless specified otherwise)

Subscripts used with coefficient symbols denote derivatives as follows:

c_ with respect to c_

& with respect to &£/2Voo

with respect to

with respect to _/2Voo

u with respect to u/Voo

15
v with respect to v/V_

W with respect to w/V_o

P with respect to p_/2V_

q with respect to q_/2V_

with respect to r_/2V_

5 with respect to control deflection

V with respect to V/V_o. This derivative, for a force such as the X-force, is
_C x
equal to 2Cx + _(V/Voo)' where the term 2C x represents the change in
8C x
X-force due to changes in free-stream velocity and the term
O(V/V_)
represents the change in CX due to effects such as Mach number or
aeroelastic effects. If these effects are negligible, then CXv = 2Cx,

Cmv = 2Cm, etc.

16
SECTION I

EQUATIONS INVOLVING BASIC FLIGHT MEASUREMENTS

This section summarizes some of the relationships used in working with basic £1ight
measurements. Included are general equations for transferring accelerations, linear
velocities , and angular velocities between any two axes systems on a flight vehicle that
are related by a @,0,q5 Euler rotation sequence; the simplified forms that are frequently
• t
used to correct accelerometer and rate-gyro readings for instrument displacement and
misalinement are given as equations (A-5) to (A-16) in appendix A. Also included are
the correction equations for angles of attack and sideslip, equations for determining flight
Mach number from measurements of static and impact pressures, and equations for deter-
mining vehicle forces and moments from accelerometer and rate-gyro readings.

GENERAL AXES TRANSFORMATIONS FOR COMPONENTS

OF ACCELERATION, LINEAR VELOCITY, AND

ANGULAR VELOCITY

The general equations given here transfer between the two axes systems shown in
figure 1, where the "instrument-axes" system represents any system displaced and/or
misalined with respect to the vehicle reference axes at the c.g. It should be noted that the
Xi,Yi,Z i components of the instrument-axes system need not be mutually perpendicular
(orthogonal) nor referred to the same origin. This would be the case, for example, if
accelerometers were used to measure the separate components of acceleration but the
accelerometers were not orthogonally alined and were located at different points on the
vehicle. Hence, three separate instrument axes have to be considered. For the equa-
tions developed here, one instrument is assumed to be alined along the X-axis of a sys-

tem located at point :_X,_X,_,X and alined at Euler angles _PX,_X,q5 x with the vehicle
reference axes; a second instrument is assumed to be alined along the Y-axis of a sys-

tem located at point _y,:_y,_y and alined at Euler angles _y,0y,_y with the refer-
ence axes, and so on. Nine displacement coordinates and nine Euler angles in all are
needed to define the locations and alinements of the three instruments.

Both axes rotation and axes translation can be performed by using these

equations (in the acceleration equations, axes-translation terms are of the form
-- dV
ainertial = _ + w × V; in the velocity equations, they are of the form

..-" dE
Vinertial = -_- + × r). The transformations from the vehicle reference axes to the
instrument axes are given in equation form; however, those from the instrument axes

17
SECTION I - Continued

to the vehicle reference axes are given only in matrix form, since they involve a matrix
inversion which becomes long and complicated unless certain simplifying assumptions
are made. Two simplified cases, in which (1) instrument axes are orthogonal and (2)
misalinement angles are small, are given as instrument-correction equations in equa-

tions (A-5) to (A-16) of appendix A.

Transfer From Vehicle Reference Axes to Instrument Axes

Acceleration. -

- ( r2 + q2)Xx + (Pq- r):_X + (rp + _l)_.x_COS 0x cos _PX


gAx, i = EgAx,cg

+ (pq + {')9_y- (p2 + r2)_y + (qr-_)_y_COS 0x sin _X

+ <pr- q)XZ + (qr + P)Yz-(q2 + p2)_z]Si n OX (I-1)

+ (pq- I')9x + (rp + q)Zx](COS _y sin Oy sin _y - sin tpy


cos Oy)
gAy,i= EgAx,cg - (r2+ q2)_ x

cos _by)
+ EgAy,cg + _pq+ _,):_y _(p2 + r2)_y+ (qr- _)_y_(Sin d/y sin 0y sin qby+COS _y

(I-2)
+ _Az,cg + (pr -_l)_z + (qr+ _))Yz- (q2+ p2)_z]COS 0y sin _by

gAz, i = _Ax,cg_ (r2 + q2)iX+ (pq- r)YX + (rp + q)ZX_ (cos _Z sin 0z cos _bZ + sin _Z sin _bz)

÷ _Ay,cg + (pq + _)iy _(p2 + r2):_y + (qr-_)_y_ (sin @Z sin 0 z cos qSZ - cos _PZ sin qbZ)

(I-3)
+ _Az,cg + (pr- 51)Xz+ (qr+ _))Yz- (q2+ p2)_z]CO s 0Z cos q5z

Linear velocity.-

_ :(_+_x- r_x)C°S
Oxcos_x+(_+r_-_)cosOX_ _X
(I-4)
-(_ + PYz - qXZ) sin Ox

18
SECTION I - Continued

vi = (6 + qZX-rYx)(sin *Y sin 8y cos _y - sin _y cos *y)

+(_+r_ -p_yl(Sin _y sin ey sin ,y + cos _y cos *YI

(I-5)
+ (W + P:_ z - q_zlsin ,y cos ey

_ =(_+_x _x)(c°_*z cos_z sin 0 Z + sin _Z sin ,Z 1

+ ('_ + r.Xy- pZy)(sin _hz cos *Z sin 0 Z - cos _hZ sin *Z)

(I-6)
+(_+_z _z)c°s_z co_0z
Angular velocity.-

(I-7)
Pi = p cos 0x cos _X + q cos 0X sin _X - r sin 0X

qi = p (COS _y sin 0y sin ,y-Sin _y cos ,y)

r sin ,y cos Oy (I-8)


+ q(sin *y sin _y sin 0y + cos _y cos *y)+

r i = p(cos _Z cos *Z sin 0 z + sin _Z sin *Z)

+ q(sin _Z cos *Z sin 0 Z - COS _Z sin *Z) + r cos 0 z cos *Z (I-9)

Transfer From Instrument Axes to Vehicle Reference Axes

If the axes along which data are measured are not orthogonally alined, the following
procedure must be used to transfer instrument readings to the vehicle reference axes

at the c.g.:

Accelerations are given by

gAx,cg gAx,i

(I-10)
= _nonorthogonal_ -1 gAy ,i + K2
gAy ,cg

gAz,cg gAz ,i
u

19
SECTION I - Continued

angular velocities, by

and linear velocities


= EFnonorthogonal_

by
-1
Er']
qi

i
(I-11)

where
Iil- - _- -1 . (I-12)

K1 = _x(r 2 + q2) _ :_X(pq _ _) _ _x(r p + _) "_

K2 = _:_y(pq + _) + _y(p2 + r 2) _ _y(qr - _)

K3 = -Xz(P r - Cl) - YZ(q r + P) + ZZ(q 2 + p2)


(I-13)
K 4 = ry x - qz x

K 5=p_Y - r:_y

K6 = q:_z - PYz
J

In equations (I-10) to (I-12) _Fnonorthogonal_ -1 is the inverse of a fransforma-


tion matrix made up of the sines and cosines of nine Euler angles, one set _,O,q_ for
each of the three instruments. This transformation matrix is defined as

cos Oxcos _X cos Oxsin _X -sin Ox

cos J/y sin _y sin Oy sin _y sin q_y sin Oy sin qSy cos Oy
(I-14)
_nonorthogonal_ = - sin _y cos qSy + cos _y cos _by

cos _/Z cos _bz sin Oz sin _Z cos _Z sin OZ cos Oz cos _Z

+ sin _Z sin q_Z - cos _Z sin _bz

Except under certain simplifying assumptions, the inverted matrix rFnonorthogonal_ -1


needed to solve these equations is too long and involved to write out; however, the matrix
inversion can be performed fairly easily on a digital computer. Equations for two cases
in which the matrix can be inverted are given in appendix A as equations (A-5) to (A-16).

2O
SECTION I - Continued

CORRECTIONS TO ANGLES OF ATTACK AND SIDESLIP

Angle of Attack
For a vane-type sensor, the true or corrected angle of attack (referred to the vehi-
cle reference axes at the c.g.) is given by
s = s i - s a - s u + s v - K o - s b+ s d +As (I-15)

where AS is the position error, the error due to the location of the sensor in the flow
field of the body, and the other corrections are those due to boom bending, misalinement,
and flow-field effects (see the list of symbols).

Equation (I-15) is for an angle-of-attack vane mounted on an instrument boom extend-


ing from the body or a wing tip but generally applies to any type of sensor. Other types
of sensor and methods of calibration are discussed in reference 5. The importance of the

different types of error is discussed in reference 3.

The angle s d is the correction for vehicle rotation and is given by

Sd = 2_ ql_ 2__ pl_ (I-16)


2V_ _ 2V_

The other errors are determined by calibration as discussed in reference 5.

Angle of Sideslip

The equation for corrected sideslip angle is

/_ = /_i + fia +/_u - _v +/_b + /_d + &_ (i-17)


where the correction for vehicle rotation is given by

2_ r_ 2£ pl_
(I-18)
_d = _ 2Voo + _ 2V_o

DETERMINATION OF FREE-STREAM MACH NUMBER

For flight tests of high-speed aircraft or missiles, the following relationships are
used to determine Mach number from onboard measurements of impact pressure qc
and static pressure 15: For subsonic conditions (M < 1.0),

qc _ - 1 _-1
(I-19)
-_-= 1 +-_M - =(1 + 0.2M) 7/2- 1

For supersonic conditions (M > 1.0), the equation is modified to include the loss in total
pressure behind the shock wave and becomes

21
SECTION I - Continued

qc = _+1 M2 - 2 (. 5.76M 2 h 5/2 -I (I-20)


2 2_ 3/+1
M2 M2
3/- 1)3/-1 i = 1.2M2\5.6M2 _ 0.8)
3/+1 3/+

The indicated Mach number M is determined from tables or plots, based on equa-

tions (I-19) and (I-20), of _c/15 as functions of Mach number.

DETERMINATION OF AERODYNAMIC FORCES AND MOMENTS FROM

ACCELEROMETER AND RATE-GYRO READINGS

The components of the total force acting on the flight vehicle can be determined by
multiplying the corrected accelerometer reading (i.e., corrected for displacement and
misalinement) by vehicle weight. In coefficient form

= WAx,cg-_

CX qooS

WAy,cg_

Cy q_o S _ (I-21)

cz WAz,cg|
qoo S

The total force thus determined includes the static and the dynamic aerodynamic forces,
the engine thrust, the jet damping force, and the reaction and other control forces. It
does not include the components of vehicle weight, however, since gravitational effects
appear as accelerations that are measured directly by the accelerometer; that is, the
products of the weight and the corrected accelerometer readings are

WAx,cg:W({l+wq-vr+g sin Og)=_F x

WAy,cg= WCc+ur-wp-g cos 0g sin _bg)= _Fy (I-22)

W_z,c_=w(_+_, u_ _ cos_ cos,_) =_ _z


where the sine and cosine terms are components of vehicle weight.

The total moments are determined from the rate-gyro readings according to the
equations

22
SECTION I - Continued

_Mx = _)Ix- _tIxy - _'Ixz + q(rI z - PIxz- qIyz ) - r(qIy - rIyz- PIxy )

_My : _Y- _z- _ +r(_Ix-_x_ - r_xz)- _(r_z-_xz- q_Yz) (I-23)

_Mz =_z- _Ixz-_z +_@I_-r_z- _x_)-_(_Ix-q_x_- r_xz)


The moment coefficients are

V" Mx
C/ = qooSb

My
(I-24)
Cm = _ qooS-'--
_ >

Cn = q_Sb

(Although b and c are commonly used to nondimensionalize the moments, these lengths
are arbitrary and any convenient length can be used.)

Methods by which the total forces and moments are broken down into trim, static,

and cross-coupling components so that longitudinal and lateral coefficients can be deter-
mined are quite complicated and beyond the scope of this paper. The simplified method
for determining basic static- and dynamic-stability derivatives from oscillatory flight
motions is developed in appendix D. In this and in most methods, linear systems are used
to represent transient flight data; these linear systems are obtained by assuming small
angles, constant aerodynamic coefficients, constant free-stream conditions, rigid-body
mass and inertial characteristics, and separation of the longitudinal and lateral modes.
For data that do not fit these limitations, exact methods of simulation based on the equa-
tions of motion for six degrees of freedom are used. Other methods that are now becom-
ing popular involve parameter identification in which linear and quasilinear estimation
techniques are used. (See, for example, refs. 2 and 6.)

DIRECT FLIGHT MEASUREMENTS OF LIFT AND DRAG

Lift and drag components for an airplane can be measured directly by using flight-
path accelerometers which differ from body-fixed accelerometers in that they are oper-
ated by a vane or pressure sensor that rotates the sensor into or normal to the stream
direction. Accelerometers measuring parallel and normal to the flight path measure

23
SECTION I - Concluded

vehicle thrust minus drag and lift, respectively, if the component of thrust in the lift
direction can be assumed negligible. In equation form the lift and drag are
TD - D = -WAD_

L = WAL _ _(I-25)

where TD is the component of thrust along the flight path. The accelerations should be

corrected for the usual boom-type errors as described in correction with equation (I-15)
for the angle of attack measured from a boom. Details of this flight-testing technique
are given in refereace 7.

24
SECTION II

TRANSFER OF AERODYNAMIC FORCE AND MOMENT COEFFICIENTS

AND DERIVATIVES TO ANOTHER REFERENCE CENTER


.4

The equations in this section are the general forms for transferring aerodynamic
force and moment coefficients and stability derivatives from a coordinate axes system
with origin at the c.g. to a parallel axes system with origin at a point x,y,z away from
the vehicle c.g. (Simplified forms obtained by assuming zero angles of attack and side-
Slip and neglecting aerodynamic cross derivatives are given as equations (A-17) to (A-43)
in appendix A.) In these general equations, unprimed coefficients are referred to the c.g.;
primed coefficients are referred to the point at x,y,z shown in figure 2.

The transformations are given for both systems of variables used in aerodynamics.
Equations for the _,/3,V,p,q,r system are given by equations (II-l) to (II-40); those for
the u,v,w,p,q,r system are given by equations (II-41) to (11-76). (The relationships
between the two sets of variables are given in table II.) The equations are derived in
appendix B. They are general in that no assumptions are made as to angle of attack or
sideslip and in that all transfer distances, lateral as well as vertical, are included; how-
ever, they are still not complete in that, for the purposes of these transformations,
second-order derivatives with respect to time are assumed negligible and are omitted.
Also, the transformations given for the static forces and moments (derivatives with
respect to _ or _) are the simplified forms that apply only when the body is not under-
going any significant rotation and p, q, and r are essentially zero. The transforma-
tions for static-stability derivatives that apply when there is significant vehicle rotation
can be derived as indicated in appendix B.

TRANSFORMATIONS FOR _,_,V DERIVATIVES

X-Axis Force Coefficients and Derivatives

Cx = Cx (II-1)

C T
Xot Cx_ (II-2)

= (II-3),

Cx& Cx&
(II-4)

T
(n-s)

25
SECTION II - Continued

C T

xp = Cx p _ 2Cx_ _ cos + 2Cx_ (_ cos _ + y sin (2 sin


COS

sin (2 cos _)
(n-6)

COS (2 + + 2C x cos (2 -
sin sin
cos _ CX(2

(II-7)
+ 4C x sin (2 - _ cos (2 cos

C T
Xr = CXr - 2Cx(2 _ sin
cos (2 2Cx_ [ cos _ +

+4C X .cos(2 cos_-_-sin_ (II-8)

Y-Axis Force Coefficients and Derivatives

c_ --
cy
1

Cy(2* %(2 /
' C
(H-9)
Cy_ = y_

c_: c_ |
C_ = Cy_

cos---V _

(II-lO)

'
Cyq = Cyq + c-_s-_\_ cos (2 + _ sin (2 + 2Cyfl cos (2 - _ sin (2 sin

+ 4Cy sin (2 - _- cos cos _ (II-11)

26
SECTION II - Continued

Cyr' =CYr- 2Cy_ _sin_-2Cy_ Cosl3+ cos _ sin

(II-12)
+ 4Cy cos ol cos ¢1 - _- sin

Z-Axis Force Coefficients and Derivatives

Cz = Cz

C' = C
Z_ Z_

C T = (II-13)
z_ Cz_
T

Cz& = Cz&

c_ - CZ_ _

- -- 2Cz Zcos _+ sin_ sin


C' Zp = C Zp 2Cz_ y_
_ cos
cos __ + _

(II-14)
+ 4Cz/_- sin _ - _- sin _ cos _)

C T

Zq = CZq
_cz_c_ _ _)_cz_(_
+ c---_\_ cos el + 7 sin +
£
cos _ - [ sin sin

(II-15)
+ 4C z sin ol - _ cos cos t3

C' Zr = C Zr - 2Czo _ __ sin


cos ol 2Cz/3 cos _ + cos (_ sin

. coscos (iI-16)

X-Axis Moment (Roll) Coefficients and Derivatives

(II-17)
,
C l = C l +_ _ Cy - Cz

C' z § (n-18)
l_ = Cl_ + _ Cy_ - [ Cz_

27
SECTION II - Continued

C'z_ = Cz_ +_cY_


z - _ CZ/_ (]/-19)

C' =
l& Cl& + _" Cy& - _ CZ& (II-20)

:2

(II-2i)
ctk : cl_ + [ cy_- _ Cz_

c':
Ip C/p +_ _'
Cyp - _ CZp __cos
£ cos _(c
/7 lot + _ -_CzD
_ Cyot

+ 2(_ cos /_+_ sin ot sin fl)(C/f_+ _ CYt3 -_cz_)


+,(f
_n_-_s,n
oco_ (II-22)

C'lq = C/q + _-Cyq


z -_ CZq +c-g_ cos ot + _ sin ot + [ Cyot- _ Czot)

s,o cos x_-


s,n -_czD
+ _ Cy_

+4cos/_ sinot -_cos C l+ _c_-_cz) (II-23)

lr Clr + _ CYr - CZr - 2 :_ cos


sin ot
_ lot
+_
_c_-_czD
cos _+_cos ot sin/_)(C//_ + [ Cyfi -_cz_)
cos ot cos /7 - _ sin C1 + _
_- Cy _ .c z (II-24)

Y-Axis Moment (Pitch) Coefficients and Derivatives

, _
C m=C m+_.c z -_-Cx (II-25)

C' x
mo t = Cmo t + _ Czot - _ Cxot (II-26)

28
SECTION II - Continued

, 9` (II-27)
Cm_ = Cm_ + _ Cz_ - _ Cx_

C' = 9` (II-28)
m& Cm& + [ CZ& - [ CX&

, 9` £ (IIL29)
Cm_ = Cm_ +[ CZ_ -[ CX_

C'
rap= Cmp+_CZp
9` _
-_CXp
_ cos a (C mot + _- Czot - _- CX
- 2 _ cos_

cos+ o
+2(} +9`i-Cz¢ _
- _Cx _)

(n-3o)
_
_c z - _c x_ )

C' = Cmq + 9`
[ CZq z CXq + 2 (_ cos ot + [_' sin ot)(c mot + [Czot - [c x

ot_9`
+ 2 sin/3 (} cos _ sin ot) (CI/3 +_Cz_ - 9. Cx_)

+ 4 cos /3(_ sin ot -


cos +_Cz _ )
-_Cx
(ii-31)

C'm r = Cmr + _ CZr - _z CXr - 2 _:_sin


cos ot
/3 (C mot _
+ _ Czot - _- Cxot )

-2/} cos /3 +} cos ot sin _) (Cm/3 + _ c z_ - _- Cx

+ l _cos cos -_sin m +_- CZ -[C


(II-32)

Z-Axis Moment (Yaw) Coefficients and Derivatives

(II-33)
,
Cn = C n + _Cx x
- _ Cy

C' 9, (II-34)
not = Cnot + _ Cxo t - _Cyot

29
lea _D L"- co o
I i I i i I I !
I,iI l,iI I=M I.IM I-..4
lit I,iI
v I..,,.q I-q I.,,..I
v v v v
I N I_,_
!

b4
D

+
C) L)
INI,_ eD 'NI'_
I
I r/l
D D
N _NI'_ _NI'_ D N 'NI,_
C)
I r/1 I ,NI,_ I
_NI'_
I_l,_
co_ I
I '_1,_ qo_ !
+ + N + N
N N c.)
(D
@ L.)
_1'_ '_1'_
°,iI @ + I I:_1_
o +
i _1,_
@ +
'_I_
+ + + c_
O qo_ qo_
r_
D ,N I,_
o
@ "_1 _ V
I .P-I
o
IiI O °_,,,I
r/1 r/1
._.,4 o O ._=4
I r/1 o + o ! r/1 r/1 r_
O O
I,iI "qo_ ,NI,_ r_
I I r/1 ! O
r_ ._=4 D r..) O
I,=,4
,NI,_
'NI'_ 'N I'_ I NI,_4 r/l O
I
,_I_ r_ r/l O
I I I I
O .liI I ©
+ + o
"cO_ I o
N
N
c) N # INI,_ I>4 O
r/1
°r-I

qo..
o r/1 r/1
_1'_ i I>,1,_ 0 O o
i _1'_ ' _,1'_ o r/1 ' _I_ r/l _I_1'_ !
+ N
+ + + _ I,_ IN I,_4 + © +
-qo..
L.) D L_ + + + 0 r..)
÷ I +
II II II II II II II 11
qo_ ._ f-i
.N _.
L.) r..) D

o
SECTION II - Continued

C T (II-43)
Xw CX w

(I1-44)
CXp=CXp +KCXv -_CXw

(ii-_5)
C_q = Cxq - K Cxu + _ Cxw

(II-46)
CXr = CXr + CXu - _ CXv

Y-Axis Force Coefficients and Derivatives

(II-47)
Cy u CY u

(II-48)
Cy v CY v

(II-49)
Cy w CY w

, g _ _ CY w (II-50)
Cyp = Cyp + [ CYv

, _
(II-51)
Cyq= Cyq -[CYu÷[CYw

(n-52)
Cyr = CYr + CYu - _ CY v

Z-Axis Force Coefficients and Derivatives

C ! (II-53)
Zu CZu

C T = (II-54)
ZV CZv

C t = (II-55)
Zw CZw

C' = + _ (II-56)
Zp CZp [ CZv - _CZw

C' = C (II-57)
Zq Zq - _"Cz u +[ CZ w

C' = C + (II-58)
Zr Zr CZ u - K CZv

31
SECTION II - Continued

X-Axis Moment (Roll) Coefficients and Derivatives

C'lu = C/u + z_ CYu - CZu (II-59) ,

C'lv = Clv+[CYv
z - _ CZv , (II-60)

C' (II-61)
lw= C/w+[Cy w-_cz w

C'lp = C/p + _-
_ Cyp -_ C Zp +[t_C lv +_CYv
_ -} CZv )

-_ :_(C lw +z _CYw -_-:_CZw ) (II-62)

C'
lq = C/q + _ Cyq -_
tc Zq-[ 1u +ICY u -£Cz u )

+ x (Clw + _CYw- } CZw) (II-63)

c' +_ (c ,u+_ C u- CZu


)

+_ (II-64)

Y-Axis Moment (Pitch) Coefficients and Derivatives

C T =
mu Cmu + C zu - [ Cxu (II-65)

C' = _i (II-66)
mv Cmv +_-CZv -[CXv

C' = _" (II-67)


mw Cmw+_CZw-_Cx w

C' mp = Crop x
+[CZp z
-_'CXp +_-z (C m v + _ c Zv - _ Cxv

- __. (Caw +9__- C Zw z Cxw )


- _- (II-68)

C' mq = Cmq+_CZq_ -_ _ Cxq - __ (c mu + _- C z u - _ CXu

+_- mw _-CZw-_-CXw (II-69)

32
SECTION II - Concluded

C T
m r = Cmr + _- CZ r - _- Cx r + mu ?" Cz u - _-Cx u )
- (Cmv
+:_ _ _. Cxv ) (n-7o)

Z-Axis Moment (Yaw) Coefficients and Derivatives

C T
= Cnu + _ CXu - _- CY u (II-71)
n u

C' (II-72)
nv = Cnv + _ Cx v - _ CYv

C' _9, (II-73)


n w
= Cn w+_Cxw ICY w

C' +_Cxv
np = Cnp + CXp -[Cyp +[ nv _ _ CYv )

(II-74)
-_(Cnw + _ CX w - _ CYw)

c_ =Cnq+_CXq -_Cyq-_ _(cnu + _- CX u _" CY u

+ _:_"(Ca w + _ Cx w - CYw )
(ii-75)

C'
n r = Cnr+_CXr [CYr+ (Cnu + _ CXu _ x CYu )

_ X(Cnv
_ + _ CXv _ -_
:_ CYv )
(ii-76)

33
SE CTION III

TRANSFER OF" AERODYNAMIC FORCE AND MOMENT COEFFICIENTS

AND DERIVATIVES TO A ROTATED AXES SYSTEM

The equations in this section are the general forms of the axes transformations for
aerodynamic force and moment coefficients and stability derivatives that can be used to
transfer between any two of the five axes systems used in aerodynamic analysis. (Par- ,
ticular forms for transfer between body and wind-tunnel stability axes are given in

appendix A.) Transformations are presented for both _,_,V derivatives and u,v,w
derivatives. (For the relationships between the two sets of derivatives, see table II.)
The equations are developed in terms of a generalized angle of attack A and generalized
sideslip angle B; transformations between any two axes systems can be obtained by sub-
stituting specific angles for A and B in these equations. (See tables III and IV.)

DESCRIPTIONS OF AXES SYSTEMS

The five axes systems considered are those shown in figure 3. While three of
these systems, body, principal, and wind axes, are clearly defined in the literature, there
is some confusion concerning the definition of the stability axes. The stability axes
described in some reports are vehicle or flight stability axes about which the equations

of motion are written; in other reports they are wind-tunnel stability axes about which
aerodynamic data are measured in the wind tunnel. The differences between the two are
pointed out in these brief descriptions of the axes systems.

Body Axes

The orthogonal body-axes system is fixed within the vehicle with the X-axis along
the longitudinal center line of the body, the Y-axis normal to the plane of symmetry, and
the Z-axis in the plane of symmetry. This is the axes system about which aircraft
instruments are usually mounted, Its main advantage in motion calculations is that vehi-
cle moments of inertia about the axes are constant, so that the I terms can be omitted

from the equations of motion. It is the logical system to which to refer velocities, accel-
erations, and stability and control parameters in the study of aircraft handling qualities
because the pilot's orientation with respect to this frame is fixed.

Principal Axes

The principal axes are an orthogonal body-fixed system for which the products of
inertia are zero. The X and Z principal axes lie in the plane of symmetry; the angle

34
SECTION III - Continued

between the X body axis and the X principal axis is usually small so that in many

cases the body axes can be assumed to coincide with the principal axes.

Flight Stability Axes

The flight stability axes (sometimes referred to as vehicle stability axes) are _n
orthogonal body-axes system fixed to the vehicle, the X-axis of which is alined with the
relative wind vector when the vehicle is in a steady-state trim condition but then rotates
with the vehicle after a disturbance as the vehicle changes angle of attack. This system
is preferred in many stability studies because, as with other body-fixed axes, the moments
of inertia about the axes remain constant and also because the motions defined are pri-
marily those about the flight path rather than about body reference lines.

Wind-Tuimel Stability Axes

The wind-tunnel stability axes are the system about which most wind-tunnel data
are obtained. For this system the X-axis is in the same horizontal plane as the relative
wind at all times (fig. 3). The angle ol between the X-axis of this system and the X
body axes is variable. (It is a constant s 0 for the flight stability axes.) This means
that vehicle moments of inertia about the X-axis change. It also means that additional
terms are required in the transformation equations for static-stability derivatives and
for u,v,w derivatives when data are transferred to or from the wind axes or the wind-
tunnel stability axes. These additional terms are designated (1), as in equation (III-7),
for example.

Wind Axes

The wind axes are the system generally used in calculating motions of the vehicle
as a point mass. The X-axis for this system is alined with the relative wind at all times
so that vehicle moments of inertia about this axis change. As with the wind-tunnel sta-
bility axes, additional terms, designated (1), are required in the transformations to or
from the wind axes and either the body, principal, or flight stability axes, since the angle
A between the X wind axis and the X-axis of either of these systems is variable. Also,
since the lateral angle B between the X-axes is variable, there are additional terms,
designated (2), as in equation (III-13), required in the transformations for some of the
lateral derivatives between the wind axes and either of the other axes systems.

NOTES ON USE OF TRANSFORMATION EQUATIONS

! T t

In the transformations that follow, symbols such as Cx,Cy,C z and Cx,Cy,C z


are used in a general sense to designate coefficients and derivatives about corresponding

35
SECTION III- Continued

X',Y',Z' and X,Y,Z axes systems; specific designations to use with each axes system
are given in table V. Two types of transformation are given: "Direct" which include
those defined in table III and "inverse" which include those defined in table IV. Trans-
formations between any two axes systems are obtained by selecting the proper angles for
A and B from table III or IV and then the proper coefficient designations from table V.
For example, in transferring from flight stability axes to body axes (a direct transforma-
tion, according to table III), angle A is replaced by a 0 and angle B equals zero; the
transformation for CX then becomes, by using equation (III-1) and table V,

C x= CX, s cos a 0- CZ, s sina 0

As a second example, in determining the C/p about the wind-tunnel stability axes
from derivatives given about the body axes (an inverse transformation according to
table IV), equation (III-82) and tables IV and V are used to obtain

Clp,w t = Clp cos 2 a + Cnr sin 2 a + (Clr + Cnp)Sin a cos s

In the transformations for static-stability derivatives and also in the transforma-

tions for u,v,w derivatives, the terms designated (1) are included only in transferring
from either wind or wind-tunnel stability axes to either body, principal, or flight stability
axes and not in other transformations. For example, terms (1) are includedwhen the

body-axes derivative CXa (eq. (III-7)) is determined from data given about wind-tunnel
stability axes as in

CXs -- CXs,w t cos s - CZs,w t sin s - Cx,wt sin s - Cz,wt cos s

but are not included when CXs is determined from data about flight stability axes as in

CXs = CXs,s cos s 0 - CZs,s sin s 0

The terms designated (2) are included only in transferring to or from wind axes.

For example, in transferring from wind axes, the body-axes derivative Clfl (eq. (III-16))
is given by

Cl[ 3 = Clfi,w cos s cos fi - Cm/_, w cos s sin fi - Cnfi,w sin s

- Cl,w cos s sin fi - Cm, w cos s cos _fi

but, in transferring from wind-tunnel stability axes, becomes

Clfi = Clfl,w t cos s - Cn_,w t sin s

If the need arises, derivatives with respect to s,[3,V can be converted to deriva-
tives with respect to u,v,w by using the relationships given in table II.

36
SECTION III - Continued
DIRECT TRANSFORMATIONS

Static Force and Moment Coefficients 1 (Direct, Table III)


T

C x= C xcosAcos B - Cy cosAsinB - C Z sinA (III-i)

Cy=C xsinB+CY cos B (III-2)

C z = C xsinAcos B - Cy sinAsinB+C z cosA (III-3)

C l=C 1 cosA cos B - C m cosAsinB- Cn.sinA (III-4)

C m= C 1 sinB+C m cos B (III-5)

C n=C l sinA cos B - C m sinAsinB+C n cosA (III-6)

Static-Stability Derivatives (Direct, Table III)

C T
X_ = CX_ cos A cos B - Cy_ cos A sin B - CZ_ sin A

-_C xsinA cos B+CysinA sinB - C z cosAj (III-7)


V

(1)
where the terms designated (1) are included only in transferring from either wind or wind-
tunnel stability axes to either body, principal, or flight stability axes and not in other
transformations.

Cy_
T ,...
CX_ sin B + Cy(_ cos B (III-8)

C T
Z_- CX_ sin A cos B - CYo _ sin A sin B + Czo l cos A

+C xcos A cos B - Cy cosA sinB - C z sinai (III-9)

(1)

1 Equations for control derivatives and velocity derivatives have this same form;
for example,

C£5 a CX5 a cos A cos B - Cy5 a cos A sin B - CZ6 a sin A

C' = cos A cos B cos A sin B sin A


XV CX V - CYv - C Zv

37
SECTION III- Continued

Cl_ = Cl_ cos A cos B - Cm_ cos A sin B - Cn_ sin A

(m-lo)
- C 1 sinAcos B +C m sinAsinB - C ncos A.2

(1)

'
Cm_
=
Cl_
sin B +
Cm_
cos B (m-11)

C T _-
n_ Cl_ sin A cos B - Cm(_ sin A sin B + Cn_ cos A

+C l cosA cos B -C m cosAsinB- C nsinA (m-12)


k. 2

= Cx/_ cos A cos B - Cy/_ cos A sin B - Cz/_ sin A

(III-13)
-%C xcos A sinB - Cy cos A cos B
_J

(2)

where the terms designated (2) are included only in transferring to or from wind axes.

C_ = CX_ sin B + Cy_ cos B + Cx cos B - Cy sin B_, (III-14)


(2)

C' = sin A cos B sin A sin B + cos A


Z;3 CX;3 - Cy; 3 Cz;3

(III-15)
-_C X sinA sin B - Cy sinA cos BJ
g

(2)

sin A
Cl/_ cos A cos B - Cmfi cos A sin B - Cn/_

(III-16)
-%C1 cosAsinB - C m cosA cos BJ

(2)

(I!I-17)
= Cl_ sin B +Cm/_ cos B+Ckl cos B - C m sin BJ
(2)

38
SECTION III - Continued
t

Cn_ = Cl_ sin A cos B - Cm_ sin A sin B + Cnfi cos A

-C 1 sinA sin B - Cm sinA cos B2 (IiI-18)


y

(21

Dynamic-Stability Derivatives (Direct, Table III)

= _Cx_cos2__ (Cxq+ Cyp).Sin B cos B + Cyq sin 2 B_cos 2 A

cos A (III-19)
+ CZr sin 2 A - E(CXr + CZp) cos B - (CYr + CZq) sin B] sin A

= _CXq C°S2B- Cyp sin2B+(CXp-Cyq)sinBcosB_cosA

- (CZq cos B + CZp sin B)sin A (III-20)

C£r = (CXr cos B- CYr sin I3)cos 2 A- (CZp cos B- CZq sin ]3)sin 2 A

cos A (III-21)
+ECXp c°s2 B + Cyq sin2 B - (CXq + Cyp) sin B cos B - CZr_ sin A

Cyp yp cos 2 B- CXq sin 2 B + (CXp - Cyq)Sin B cos cos A

- (CYr cos B + CXr sin B) sin A (III-22)

C' = cos 2 B + sin 2 B + (III-23)


yq Cyq CXp (CXq + Cyp)sin B cos B

Cyr = (CYr cos B + CXr sin B)cos A

+Ec_cos_
_-c_ s_n_
_+(Cx-c_sin_co_._s*nA (III-24)

C !

Zp (CZp cos B- CZq sin B)cos 2 A -(CXr cos B- CYr sin B)sin 2 A

cos2
o+c, +c, ls ncos
o-CZ s n cos A (III-25)

39
SECTION III- Continued

C t __
Zq (CZq cos B + CZp sin B)cos A

+ ECXq cos2 B - Cyp sin2 B + (CXp - Cyq)Sin B cos B_sin A (III-26)

Zr CZr c°s2 A + Xp c°s2 B + Cyq sin 2 B- (CXq + Cyp)sin B cos sin 2 A

+ _(CXr + CZp) cos B-(CYr + CZq)Sin B_sin A cos A (III-27)

C'lp =_Clp c°s2B-(C/q+Cmp) sinBcosB+Cmqsin2B_cos2A+Cnrsin2A

+ _(C/r + Cnp ) cos B+ (Cmr+ Cnq) sin B_ sin A cos A (III-28)

C'lq =_Clq c°s2B" Cmp sin2B+(Clp-Cmq)sinBcosB_cosA

- (Cnq cos B + Cap sin B) sin A (III-29)

C Y
lr = (C/r cos B- Cmr sin B)cos 2 A -(Cnp cos B- Cnq sin B)sin 2 A

+_Clp cos 2 B+Cmq sin2B-(C/q + Cmp)Sin B cosB-Cnr_Sin A cosA (III-30)

C !

mp Cmp cos 2 B - C/q sin 2 B + (C/p - Cmq )sin B cos B_ cos A

- (Cmr cos B + Clr sin B)sin A (III-31)

Cmq = Cmq cos 2 B + Clp sin 2 B + (Clq + Cmp )sin B cos B (111-32)

C !
m r Cmr cos B + Clr sin B) cos A

+ _Cmp cos 2 B- Clq sin 2 B + (C/p - Cmq)Sin B cos B_ sin A (III-33)

40
SECTION Ill -- Continued

(Cap cos B- Cnq sin B)cos 2 A -(C/r cos B- Cmr sin B)sin 2 A

+ _C/p c°s2 B + Cmq sin2 B- (C/q + Cmp)Sin B cos B- Car]Sin A cos A (111-34)

Cnq cos B + Cnp sin B) cos A

+ [Clq e°s2 B - Cmp sin2 B + (Clp - Cmq)Sin B cos B] sin A (111-35)

C T
n r
= Cn r cos 2 A + _Clp cos 2 B + Cmq sin 2 B - (C/q + Cmp ) sin B cos A
B] sin 2

+ E(Clr + Cnp ) cos B-(Cmr + Cnq)Sin B]sin A cos A (111-36)

u,v,w Derivatives (Direct, Table III)

CX u =_CXuC°S2B sin 2 A
-(CXv+CYu)SinBcosB+Cyvsin2B]cos2A+CZw

+ ¢Cxw
+ +/C +Cz
w s,o
js,o cos.
+l(c_cos----'_ X sin2A cos B - Cy sin 2AsinB+C z sinAcos

(1)

+
_X c°s2 A sin 2 B + Cy cos 2 A sin B cos B (111-37)
J

(2)

C t
Xv = [CXv cos2 B - CYu sin2 B + (CXu - Cyv)Sin B cos B] cos A

- (CZv cos B+CZuSinB)sinA_ (C xcosAsinB+CYcosAcos B_cos B

(2)
(111-38)

41
SECTION III- Continued

CXw cos B - CYw sin B) COS


2 A - (CZu cos B - CZv sin B)sin 2 A

cos A
+ ECXucos2 B + CYv sin2 B- (CXv + CYu)Sin B cos B- CZwJSin A

cos A [,_ sin A


\cos B k_X cos B - Cy sinA sin B+C z cos A))

(1)

t
+sin A sin B(C x cos A sin B + Cy cos A cos B)j
(III-39)
(2)

CY u = Yu cOs2 B- CXv sin 2 B + (CXu- CYv ) sin B cos cos A

-(CYwCOS B+CXwSinB)sinA-\(CxcosB- Cy sin B) cos A sin B,J (III-40)

(2)
t

Cy v = CYv cos 2 B + CXu sin 2 B + (CXv + CYu )sin B cos B

+
k(Cx cos B - Cy sin B)cos Bj (III-41)
v

(2)
T _

Cyw (CYw cos B + CXw sin B)cos A

+ECYuCOS2B-CXvSin2 B+(CXu-CYv)SinBcosB_sinA

- %(Cx cos B - Cy sin B) sin A sin B_2 (III-42)


,

(2)

=(CZu cos B- CZv sin B)cos 2 A -(CXw cos B- CYw sin B)sin 2 A

cos A
+ [CXu cos2 B+CYv sin2 B-(CXv +CYu)SinB cosB- CZw_SinA

-_,(C x cosA cos B - Cy cosA sinB - C z sinA) sinA


cos B2
V

(1)
+%(C xsinAsinB +Cy sinA cos B) cosAsinB (III-43)
J

(2)
42
SECTION III Continued

C _ ___
Zv (CZv cos B + CZ u sin B_/cos A

+ ECXv cos2 B - CYu sin2 B + (CXu - CYv) sin B cos B_ sin A

(III-44)
-k(Cx sin B + Cy cos B) sin A cos B
_2
(2)

C t
= CZw cos2 A + ECXu cos2 B + CYv sin2 B - (CXv + CYu)Sin B c°s B_sin2 A
ZW

+ _(CXw+CZu)COSB-(CYw+CZv)SinB_ sinAc°sA

+%(C XcosAcosB - Cy cosA sinB - C z sinA) c°sA


COS Bj

(i)

+
(in-45)
\(C x sin A sin B + Cy sin A cos B) sin A sin BJ

(2)

Ct
lu = _C/uCOS2 B-(C/v+Cmu)sinBcosB+Cmvsin2B_c°s2A

+ Cnw sin2 A + _(CIN + Cnu)COS B + (Cmw + Cnv)Sin B_ sin A c°s A

+_C l sinAcos B - C m sinAsinB+C n cosA) SinA


COS gj

(1)
(III-46)
+iC 1 sin B + C m cos B)cos 2 A sin B
_2

(2)

C'lv =E Clv c°s2B- Cmu sin2 B+(C/u-Cmv)sinBc°sB_c°sA

(III-47)
-(Cnv cosB+CnusinB)sinA- (ClSinB+Cmc°sB)c°sAc°sB
(2)

43
SECTION HI -- Continued

(C_w
cos
__Cmvs_n
_)cos_
_-<Cnu
cos'-_nv
_n_>s_n_
_
•_C,u
cos_
_*Cmv
sin_
_-<C,v*
Cmu>_n
_cos
O-Cn_S_n
_cos
-_C l sinA cos B - C m sinA sinB +C n cos A) C°S A
cos B}
(1)

+ _(C l sin B +C m cos B) sinA cosAsinB (III-4_)

(2)

Cmu = _Cmuc°S2B-Clvsin2 B+(Clu-Cmv)sinBcosB_cosA

-(CmweosB+C/wsinB)sinA-%(C/ cosB-CmsinB)cosAsinB (III-49)

(2)

C T

m V = Cmv cos 2 B + C/u sin 2 B + (C/v + Cmu )sin B cos B

+L(C/ cos B- Cm sin B) cos BJ (III-50)

(2)

C T

m W = (Cmw cos B + C/w sin B) cos A

+_CmueOS2B -ClvSin2B+(Clu -Cmv)sinBcosB_sinA

C cos B - C m sinB) sinAsin BJ (III-51)


g

(2)

44
SECTION III- Continued

Cnu cos B - Cn v sin B_cos


2/ A - \[C/w cos B - Cm w sin B)sin 2 A
Cnu =

cos A
+ _C/u cos2 B+Cmv sin2 B -(C/v+Cmu)sinBcos B- Cnw_sinA

- L(C l cos A cos B - Cm cosA sinB - C nsin A_/-_'Bj


sin A
.

(1)

(III-52)
+ +C m sinA cos B) cosA sinB_
L(C sinA sinB

(2)

Cnv cos B + Cnu sin B)cos A

+_C/vCOS2 B - Cmu sin2 B+(C/u- Cmv)sinB cos B_sinA

(III-53)
- \(C l sin B + Cm cos B) sin A cos BJ

(2)

C !
nW =CnwCOS2
A+_CluCOS2B+CmvSin2B
-(Clv+Cmu)sinBc°sB_
sin2
A
+_C_,_*_nu)CO_
_-ICm_
+Co_-_n
_3_n_cos
_
+ (C l cosA cos B - Cm cos A sinB - Cn sinA) C°S A
cos

(1)

(III-54)
+
(C/ sin A sin B + Cm sin A cos B) sin A sin BJ

(2)

45
SECTION ]II- Continued

INVERSE TRANSFORMATIONS

Static Force and Moment Coefficients 1 (Inverse, Table IV)

C x= ('C xcosA+C z' sinA ) cos B+ C_ sinB

Cy = - (c£cos A + C 'z sinA )sinB +CyCOS , B (III-56)

T t

C z = -C xsinA+C z cosA (III-57)

C l= (_
C cosA+C ,
nsinA )cos B+C ,
m sinB (III
- 58)

C m= - (c;cosA+C '
nsinA )sinB+C ,
m cos B (III-59)

? T

Cn = -C 1 sinA+C ncosA (III-60)

Static-Stability Derivatives (Inverse, Table IV)

CX_ = (C' x_C°S A+ C' Z_ sin A) cos B + Cyo _ sin B

XSinA- C Z cos cos B (III-61)


)

(1)

Cy_ = -.(c'X_ cosA+ c' Z_ sin A) sin B + Cy_ cos B

+ ,(c'X sinA ' cos


- Cz A) sin B (III-62)
J
V

(i)

CZ_ = -Cx+ sinA+ c' Zo t cosA - c' X cos A- C 'z sin A


/
(III-63)

(1)

1 Equations for control derivatives and velocity derivatives have this same form;
for example,

Cx5 a
=
(Cx0 '' cos A+C
Z5 a
sinA
) cos B +
% sin B

CXv= ( CXv' cosA+C'ZvSinA )cos B+CYvSinB '


46
SECTION III- Continued

Cl_ -- t C' cosA+ C'n_ sin A) cosB+Cm_' sinB

(III-64)
-(C' l sinA- C ncosA
' ) cos B
v J
(1)

cos B
Cmo / = lc _ cosA+ n_ sin sinB+Cmc _

(III-65)
+ C' 1 sinA - C n' cos A) sinB

(I)

I._T T t t (III-66)
Cn_ = -Ul_ sinA+Cnc _ cosA- C l cosA- C nsinA ,)

(1)

Cx_ = x_
C' cos A + CZ/_ sin A 1 cos B + C'Y/3 sinB

(m-67)
l
C xcos
T
A+C T
z sinA /
sinB+Cy !
cos B
2

Cy/3 = - (cx_
' cosA+ c'Z/3 sinA ) sinB + cos B

(m-68)
xc°sA+C z sin cos B - CysinB
• j

(2)

(III-69)
C Z _ = -Cx_ sinA+ C' Z_ cosA

Cl_= t_'l_ c°sA+Cn;3sinA)c°sB+Cm/_sinB

(111-7o)
-(c' l cosA +C' n sin A) sinB ' cos
+C m BJ

(2)

47
SECTION III- Continued

Cm;3 = - Ict l_ cos A + Cn_


t sin A1 sin B + Cm;
t 3 cos B

lt cosA+C nsinA cos B- Cm sinB


(In-71)
- k.IC t 1 t

(2)

t
(III-72)
Cn_ = -Cj /3 sin A + Cn; 3 cos A

Dynamic-Stability Derivatives (Inverse, Table IV)

CXp = [C' Xp cos2A+C' Z r sin2A+(CXr+' C'Zp) sin A cos A] cos 2 B

(III-73)
+ C_q sin 2 B + [(CXq + Cyp) cos A+ <Cyr + C'Zq )sinAlsinBcos B

Cxq = (, CXq cosA+ ., ) zqsinA cos 2 B - ('


Cyp cosA+ ._.,.
) sinA sin 2 B

- C'Xp (III-74)
c°s2 A + C'Zr sin 2 A + C'X r + C'Zp sin A cos A - Cx sin B cos B

CXr __[.,
co..A_c,
('-"rl A]
Xr Zp sin2 A - CXp Z sin A cos cos B

+ cos A - sin sin B (III-75)


(Cyr Cyp A)

Cyp = ICyp cos A+ Cyr sinA ) cos 2 B - ( CXq


' cos A + C'Zq sin A) sin2 B

- C' ' C' ' (III-76)


[C'Xp cos2A+ Zr sin2A+(CXr+ Zp) sinA cosA - Cy q] sin B cos B

Cyq= cos 2 B + [c, c,


Xp c°s2 A+ Zr sin 2 A + (c,c,) Xr + Z sinA cos sin 2
B

- [(%%1 + cos A + (c,%).]


Yr + sin sin B cos B
(III-77)

CYr = (c, Yr cos A - % sin .)cos B

+ E C'
X r
cos 2 A +
c, Zp
sin2 A +
(Cx -C'r)
£ Z
sin A cos sin B (in-78)

48
SECTION III - Continued

B
CZp =_2' Zp cos2A - CXr
' sin2 A - (C Xp
' - C' Z r)Sin A cos A]cos

(III-79)
+ ¢2' Zq cos A-C_q sin A)sin B

C Zq Zq cos A - sin cos B

+ _C' Zp c°s2 A+ (III-80)


CXr "sin 2 A+ (CXp- C' Z r) sinA cos A] sin B

CZr =C' Zr c°s2 A + Ckp sin 2 A - (C£r + C' Zp )sinA cosA (III-81)

B
Clp = C'lp cos 2A+C'nr sin 2A+ ( C_r + C' np) sin A cos A] cos 2

(III-82)
+C' mq sin2B+[(C' lq +Cmp)cosA+I
r' C' m + C _q) sin A] sin B cos B

C/q= lq cosA+Cnq sin cos 2 B - C'mp cosA+Cmr sin sin 2 B

B cos B (III-83)
-_flp cos2A+C'nrSin2A+ IC'lr+Cnp)sinAcosA'
" - Cmq]Sin'

Clr = C'lr cos 2 A - Cnp sin 2 A - C'lp - Cn sinA cos cos B

+ (c'mr cos A - Cmp ' sin A)sin B (III-84)

Cmp = (C'mp cos A + C'mr sinA)c°s2B- ( C'lq cosA+C'nq sinA) sin2B

(III-85)
-_]'lp cos2A+C'n r sin2A+(C'lr +Cnp)SlnAc°sA-Cmq]sinBc°sB'
" '

Cmq = Cmq
' cos 2 B + C'lp cos 2A+Cnr ' sin 2A+ I r +Cnp ,) sinA cos A] sin 2
B

- C' + cosA+ +C sinB cos B (III-86)


[(lqCtmp) (Cmr nq) sinIA ]

49
SECTION III- Continued

Cmr =(c_r cosA - Cmp, sinA) cos B

+ F- C'lr cos 2A+Cnp , sin 2A+ lp - C'r) n sinA cos A] sin B


(III-87)

Cnp =_:' Up cos 2A- C'Ir sin2 A- I C' Ip - C' nr 1 sin A cos A_cos B

+ nq cos A- c, lq sin A) sin B


(III-8_)

= (c, nq cos A _o, A) lq sin cos B

B (III-89)
+ _C' np cos2A+C' lr sin2A+(C' lp - C'nr) sin A cos A_ sin

(III-90)
Cn r = C' nr cos2A+ C/p' sin2A- (C ir
' + Cnp ) sin A cos A

u,v,w Derivatives (Inverse, Table IV)

= C' ' C'


CXu _C' Xu cos2A+ Zw sin2A+ (CXw+ Zu ) sin A cos A_ cos 2 B

+ C'Yv sin2 B+ E(C' Xv+ C' Yu ) cosA+ (C' Yw +Cz v )sinA_sinBcosB

+
xsinA - C z cos A) sinA
,,, _fl

(1)

+ C xT cosA+C zT sinA sinB - Cy cos cosAsinB "(III-91)


_ t t g_

L J
(2)

CXv =(' CXv cosA+ c' Zv sinA ) cos 2 B - (c' Yu cosA+ '
CywSin A) sin 2 B

C' cos 2A+ sin 2A+ + sinAcosA - Cy sinB cos B


Xu Zw Xw Zu

+ cos A + Cz sin sin B + Cy cos cos B (III-92).


)
(2)

5O
SECTION III Continued

CXw C'Xw cos2 A - C'Zu sin 2 A- ( Xu Zw)sinAcos cos B

+ I_' Yw cosA ' u


- CY sinA)sinB- (C' xsinA - C Z cos A) cosA J
V

(I)

+ C cosA+C
'
Z sinA
1sin B - Cy
' cos sinA sinB
J
(III-93)

(2)

"_ C' yuC°SA+Cyw ' sinA)cos 2 B- ( C' Xv cosA+ C' Zv sin A) sin 2 B
CYu=I

- C'Xu cos 2 A + C' Zw sin 2 A+ (C' Xw+ '


Cz_sin A cos A - C'v y sin B cos B

sin Asin B
- k(C' X sin A - C z cos A)
cos B s
(1)

(III-94)
+ cos A+C Z sinA cos B +Cy sin cos A sinB
2
(2)

C' cos 2 A + sin 2 A + + sin A B


CY v = Cyv cos 2 B + Xu Zw Xw Zu cos A] sin 2

- [( C' Xv + C' Yu ) cosA+ - ( CY' w+ Z


C'v) sin A] sin B cos B

(III-95)
- cosA+C_ sinA cos B+Cy sin cos B
J

(2)

CY w = (C' Yw cosA - Cyu


' sinA)cosB

C' C' ' sin B


+ _C' X w cos2A+ Zu sin2A+ ( Xu - CZw)SinA cos

+
_.(C L sinA - C 'Z cos A) cos cosB
A sin Bi

(I)

+ sin B (III-96)
X
T cos
k.[<C A + C z sinA cos B+Cy sin sinA
J

51
SECTION III- Continued

Czu = C'Zu cos 2 A - sin2 A - CXu Z sinA cos cos B

sin A
+
C' cos A - CXv sin sin B + cos A + C z sin (III-97)
Zv "'v- cos
(1)

CZv = (C' Z v cosA '


- CXv sin A) cos B

+ CZu' cos 2A+ sin 2 A+ (III-98)


CLw ( CXu
' -C' Zw ) sin A cos A] sin B

CZ w = C'Zw cos 2 A+ CXu' sin 2 A - \_C' Xw + C' Zu_Sin] A cosA

- IoLcos ,
A + C z sin cos Bs
(III-99)

(1)

B
C/u = _'lu c°s2A+ Cnw' sin2A+ (c,lw + Cnu)
2' sin A cos A] cos

+C' mv sin 2 B + E(C' 1v + C' mu ) cosA+ @' mw+ C' nv )sinA_sinBcosB

,l sinA - C n cosA/sinAI

(i)

+ C cosA+C n sin sin B - C m cos cos A sinB (III-100)


.._]

C/v lv cosA+ nv sinA ) cos 2 B - ( Cmu cosA+Cmw sinA ) sin 2 B

- C' C' ' - '


_,_'I u cos2 A + nw sin2 A + ( I w + Cnu)Sin A cos A Cmv]Sin B cos B

- C cosA+C nsin sinB - C m cos cos B (III-101)


2
V

(2)

52
SECTION III - Continued

C/w __? lw nu sin 2 A _c, ( lu A cos B

+ (,Cmw cosA - Cmu , sinA )sinB - (c'1 sinA '


- C ncosA ) cosA Y

(11

+ C cosA+C nsinA sinB - C m cos sinAsinB (III-102)


J
(2)

Cmu = (Cmu
t cosA+Cmw v sinA) cos 2 B _ (C lv
t cosA+ C nv
t sinA ) sin 2 B

C" cos 2A+Cnwsin 2 A+ + sinA cosA - Cm sinB cos B


lu lw nu

_k( C'1 sinA - C n' cos A) sinAsinB


cos B ,y

(1)

÷ C cosA+C nsin cos B+C m sin cosAsinB (III-103)


J

(2)

Cm v= C'mv cos 2 B + C'lu cos 2 A+Cnwsin 2 A+ lw + Cn smA cos sin 2 B

- [(C'lv + Cmu)COS A+ (C'mw + C_v) sin A] sin B cos B

) o]
T T

- C cos A +C z sinA cos B + C m sin cos B (III-104)


2

(2)

53
SECTION HI -- Concluded

Crow-- t
C!
mw cos A - Cmu ' sin A ) cos B

+ Clw cos 2A+C_usin 2A+ Clu-Cn sinAcos sinB

+ sin ,
A - C n cos c-_s'B .j

(1)

+ C cosA+C n sin cos B+C m sin sinA sinB (III-105)

(2)

Cnu = _:' nu cos2A-C lw


' sin2 A - (2' lu -C'w)n sin A cos A_ cos B

+ (III-106)
QJ nv cosA-C' lv sin A1 sinB + L( C l' cosA+ C n' sinAlSinAcos B2
v
(1)

Cnv =(c' nv cos A - Clv' sin A) cos B

+ '
Cnu cos 2A+ C'lw sin 2A+ C' lu - C n sinAcos sinB (III-107)

Cnw = C' nw cos 2 A+ C' lu sin 2 A- IC' I w + Cn'ul sinA cos A

(C_ cos A + C n' sin A) cos A • (III-108)


_. cos B j
V

(1)

54
SECTION IV

TRANSFORMATION EQUATIONS FOR MOMENTS OF INERTIA

This section gives equations for transferring moments of inertia to a rotated ames
system in a general form, similar to that used in section HI, so that transformations can
be made between any two of the five axes systems used in aerodynamic analysis. As in
section III, these transformations are defined as direct and inverse in tables III and IV,
respectively, and can be performed between any two axes systems by selecting proper
,, angles for A and B from these tables.

These equations are derived from a general transformation of the form


m
t t

I_ -Ixy -Ix_ Ix -Ixy -Ixz

-Ixy Iy -Iyz Iy -Iyz EF] -I (IV-l)


t T T

-Ixz -Iyz Iz -Iyz Iz

where ,r is the transformation matrix given in appendix A and is applied to the aerody-
namic axes by setting _=-B, O=A, and _b=0.

The most commonly used forms of these transformations, those for transfers among
body, principal, and flight stability axes, are given in appendix A.

DIRECT TRANSFORMATIONS (TABLE III)

Ix= IX cos 2 B +2Ixy sinB cos B+Iy sin 2 cos 2A+I z sin 2 A

+ (2Ixz cos B - 2Iyz sin ]3) sin A cos A (IV-2)

Iy=Iycos 2 B +I xsin 2 B - 2Ixy sinB cos B (IV -3)

I z =I Z cos 2 A+ (I xcos 2 B +Iy sin 2 B+2Ixy sinB cos B) sin 2A

- (2Ixz cos B - 2Iyz sin B) sin A cos A (IV-4)

Ixy= s n2 cos A
- (I_z cos_ ÷Ixz _n _)s_nA (IV-5)

55
SECTION IV - Concluded
t

_xz--(_xzcos_ -_,z _n_)co2"A-(_XZco__ -_,Z_n_)s_n_"


A
-(I X cos 2 B +Iy sin 2 B +2Ixy sinB cos B - Iz) sinA cos A (IV-6)

Iyz = (Iyz cos B + IXZ sin B) cos A

+ _xy(eos2 B-sin 2 B)-(I X -Iy)sin B eosB_ sinA

INVERSE TRANSFORMATIONS (TABLE IV)

IX= (I x cos 2 A+ I_, sin2A - 2I_z sin A cos A)cos 2 B+ I_ sin 2 B

- 2 cos A + Iyz ' sin A) sin B cos B (IV-8)

Iy= Iy cos 2 B+ (I_ cos 2 A+ I_. sin 2 A-2I_z sin A cos A)sin 2 B

+2(Ix_ycosA+Iyz sinA) sinBcos B (IV-9)

! T t

I Z =I Z cos 2 A+I Xsin 2 A+2Ixz sinA cosA (iv-lo)

Ixy '
= ( Ixy cos A+Iyz ' sin A) cos 2 B ("Ixy cos A + Iyz
' sin A) sin 2 B

+ I cos 2 A+I ,
z sin 2 A - 2Ixz , sinA cosA - I' sin B cos B (IV-11)

IXZ = Z ( COS A - sin 2 *)+I I xt _

+ (, Iyz cos A - Ixy , sin A) sin B (IV-12)

Iyz = (' Iyz cosA - Ixy


' sin A) cos B + EI'xz(Sin2 A - cos 2 A)

-(I_-I_) sin A cos A_ sin B (iv-13)

56
SECTION V

EQUATIONS OF MOTION FOR SIX DEGREES OF FREEDOM

The equations presented here are the general forms that include the variables ,
likely to be of interest in computing motions in the vicinity of the earth (moon and sun
perturbations are ignored). The terms of each equation are grouped so that various
effects (for example, oblateness of the earth) can be accounted for by adding or omitting
certain terms. Linearized equations and the wind-axes equations for a point mass are
given in appendix A.

The general equations apply to any of the five systems of vehicle reference axes

shown in figure 3. The Euler angles _g,0g,_)g are referred to the gravity-axes sys-
tem, with origin located at the surface of the earth, which rotates with the earth as shown
in figure 4. (The relationship of these gravity axes to the vehicle reference axes is
shown in fig. 5.)

FORCE EQUATIONS

Equations for the forces along the X,Y,Z axes are given in general form as equa-
tions (V-l) to (V-3), respectively. Equations for specialized cases can be obtained from
the general forms as follows:

(1) For an oblate earth, equations (V-43) to (V-45) of auxiliary equations are used
in place of term (2) in X,Y,Z force equations, respectively.

(2) If mass of vehicle is constant (zero thrust is also implied), terms (3) to (5) are
omitted.

(3) For flight outside the atmosphere, terms (6) to (8) are omitted.

X-Axis Forces

• " - + = - r

(1) (2) (3) (4) (5) ..

1 2 + ¥ _. &t • _ilL.
+ _pV_S(Cx,o Cx_a + CXv V'oo
--+_X & _*CXq2v_oj
2Voo

(6)

(v-l)
+ yXl 3#9+ CX_ 2Voo + CXp + CXr 2--_, + CXsaSa + CXse6e + CX6rSr + 0

(7) (s) (9)

57
SECTION V - Continued

where the different terms are defined as follows:

(1) Mass times acceleration

(2) Component of vehicle weight (for a more nearly complete weight component,
one that includes earth-oblateness effects, see auxiliary eqs. (V-43) to (V-45))

(3) Primary rocket thrust (see auxiliary eqs. (V-36) to (V-39))

(4) Jet damping force

(5) Reaction control force

(6) Basic aerodynamic forces

(7) Aerodynamic cross-coupling terms

(8) Aerodynamic control forces

(9) Higher order terms

Expansions of the aerodynamic forces and moments (terms (6) to (9) in force equa-
tions and terms (11) to (14) in moment equations) neglect all aerodynamic partial deriva-
tives with respect to rates of change of velocities and angles except those with respect to
a and /L The forces are expanded in terms of the independent variables ol,/_,V,p,q,r
but could as easily be expanded in terms of the variables u ,v ,w ,p ,q ,r. The X-axis aero-
dynamic force, for example, could also be written as

"1 2 / P_ q_
X = _pV_oS _Cx, o + CXuU + CXvV + CXwW + CXp _ + CXq 2V_o

+ CXr r t + _
2Voo CX65 + Higher order terms]

Y-Axis Forces

_ sin_ _ -kmCrdx-Pdzt+_

(I) (2) (3) (4) (5)

___ _P.!- + rP_


+ _pV:oS Y ,o + Cy_ + CYv + Cy_ 2V_ + Cyp 2V_o CYr 2-'_
v
7

C6)

(V -2)
+ Cyaa + Cy& _ + Cyq
\ -._ ,. - - %M

(7) (8) (9).

where the terms (1) to (9) are defined as in equation (V-I).

58
SECTION V- Continued

Z-Axis Forces

m(w + pv - qu) -rag cos 0g co_ = TZ - +

(1) (2) (3) (4) (5)

v &_ +
+ _pV_S Z,o + CZ_a +CZv_+Cz& 2V_ CZq2v_j

(6)

+ CZ_ + Cz_ -_ + _ + -- + Cz5 5a + CZ5 6e + 0


2v_ czp2v_ CZr2v_ _/
(7) (8) (9)
(v-s)

where the terms (1) to (9) are defined as in equation (V-I).

MOMENT EQUATIONS

Equations of motion involving moments about the X,Y,Z axes are given in general
form as equations (V-4) to (V-6), respectively. Equations for specialized cases can be
obtained from the general forms as follows: _:: .....

(1) If vehicle mass is constant, terms (2) to (4) and (7) to (9) are omitted

(2) If principal £xes are used, terms (3) to (6) are omitted

(3) If vehicle has a plane of symmetry (X-Z plane) but principal axes are not used,
terms (3) and (5) are omitted

X-Axis Moment (Roll)

+ -_ _ k _ )

(1) (2) (3) (4) (5) (6)

= rn_pldy2+dz21+dx(qdy+rdz))_
_ -%(Tydz_+_ + t_
Lrrn

(7), (8) (9) (10)

v _h_ + r_
+ Cla_ + Cl& _
&_
+ C/q
q_
2Vooj
+ CI/3;3+ C/V _ + Cl# 2Voo + CZP 2V_ Clr 2--_j
\ v
v
(i2)
(11)

(V-4) .

+ C/6aSa + C/sr6r) + 0

(13) (14)
59
SECTION V - Continued

where the different terms are defined as follows:

(1), (5), and (6) moments of inertia times angular accelerations

(2) to (4) and (7) jet damping moments

(8) Moments due to main rocket thrust

(9) Moment due to reaction control

(10) Moment due to gyroscopic action of engine rotating mass


(see auxiliary eqs. (V-40) to (V-42))

(11) Basic aerodynamic moments 1

(12) Aerodynamic cross-coupling terms

(13) Moments due to aerodynamic controls

(14) Higher order terms

-Y-Axis Moment (Pitch)

Y +P + qIy - - rIyz + (pq - 9)Iyz - (qr + P)Ix-y +

(1) (2) (3) (4) (5) (6)

= _nI_q (dx2 + dz2)+ dy(p dx + r dz)_ ÷ _(Tx d z _ - TZ dx) ) +_ + Mrm


(7) (8) (9) (lO)

-I- __V + &_ q_


12 pV2S_ Q\
Cm,° + Cm°l_ + Cmv Voo Cm& _ + Cmq 2--_j

(ill

Ak P_ r_ \

+ Cmfi_ + Cm_ 2V_ Cm 6 6a+Cm,. 6e_+ 0 (v-5)


+Cmp2-V-_+Cmr2-_ +_: oe _y
(12) (13) (14)

where terms (1) to (14) are defined as in equation (V-4).

1 See statement following equation (V-l) about terms (11) to (14) in expansion of
aerodynamic forces and moments.

6O
SECTION V - Continued

_Iz +pq(Iv - Ix) + riz - qiyz -pixz


+ (pr
+ +
(i) (2) * (3) (4) (5) (6)

dxl+ _ Nrm
= l:n_r(dx 2 +dy21 +dz(P_ dx + q dY)j1 - v( TX dY LTY +

(7) (st (9) (10)

p_ r_ &_ q_
lpv2sl[ Cn,o + + V + Cnp _ + Cn r 2---_z + Cn_ + Cn_ _ + Cnq 2Vooj
\k Cn/_ Cnv _ + Cn_2V_o k

(117 (12)

(V-6)
vc)

(13) (14)

where terms (1) to (14) are defined as in equation (V-4).

AUXILIARY EQUATIONS

General equations that take into account various relationships and effects in develop-

ing the equations of motion are given as equations (V-7) to (V-45). Equations for special-
ized cases may be obtained as follows:

(1) For flat earth, terms containing and X are omitted

(2) For nonrotating earth, _2e, V x, Vy, and Vz are omitted


are omitted
(3) For no surface winds, terms containing W x, Wy, and WZ

Relationship Between Euler Angles and Angular Velocities

The Euler angles specifying vehicle alinement with the gravity-axes system can be
determined from the angular velocities p,q,r by the equations

r cos qbg q sin qbg


÷ + 1", cos _g tan 0g + (_2 e + _)(sin L
cos Og cos Og

+ cos L sin _hg tan Og)


(v-_)

0g= q COS q_g-r sin _bg-_L sin _g-(_2 e + _)cos L cos


(v-s)

61
SECTION V - Continued

_I_ c°s _g (_e + L


_g=p+qtan agSin _g+rtan 0gCOS _bg+_cos 0g + _)cos
(v-9)
The inverse relationships are

P = _g- _g sin 0g- f_ cos 0g cos _g-(_e + _)(cos L cos 0g sin _Pg

- sin L sin 8g)

q = _g cos qSg + _g sin _g cos 8g- I_(sin Og cos _g sin qSg sin _g cos _g)

- (_e + i)cos L(sin 0g sin _g sin q_g + cos _g cos _g)

- (ae + )sin L sin cos eg

r = -_g sin _g + _g cos q_g cos 0g - f_ (sin 0g cos _g cos _bg + sin _g sin qSg)

-(_2 e + _)cos L(sin Og sin _g cos _g- cos _Pg sin qSg)

- (_e + _)sin L cos qSg cos Og


(V-12)

Vehicle Coordinates

Vehicle coordinates can be computed from the X,Y,Z axes velocity components
and vehicle Euler angles by integrating the equations

_g
: u cos _g cos 0g + v(cos _g sin 0g sin _g- sin _g cos _g)

+ w (cos _g sin 0g cos q_g + sin _g sin _g) (V-13)

_g = u sin _Pg cos Og + v(sin _g sin Og sin qSg + cos _hg cos _g)

+ w(sin _g sin Og cos q_g- cos _g sin _g) (V,14)

_.g = -u sin Og. + v cos Og sin qSg + w cos Og cos qSg (V-15)

rg = _Xg 2 + yg2 + Zg2


(V-16)

62
SECTION V - Continued

Trajectory Parameters (See Fig. 6)

Flight-path angle, longitude, and latitude can be determined from the equations

Y = tan-1 h (V-17)

_(_- Re_2 e cos L) 2 + _2

_g
_2e (V-18)
R e cos L

"1_= -_'g
R--/ (V-19)
The range for a spherical earth is

Range = Rea (V-20)

where, for a total range of less than 50 000 feet,

-= o)cos,oy (V-21)

For a total range greater than 50 000 feet,

or= cos-lEsin Lo sin L + cos L O cos L cos(k ko)_ (V-22)

range for a flat earth is

Range = _Xg 2 + yg2 (V-23)

Angle of Attack, Sideslip, and Relative Velocity (See Fig. 7)

Angle of attack, sideslip, and resultant relative velocity are related to components
of velocity along the vehicle body axes by

= tan -1 Wb
fi"_" (V-24)

fi = sin -1 _7__b
V (V-25)

V = t]b2 + XTb2 + W-'b2 (V-26)

Ub = V cos _ cos /3 = u b - VX, b - WX, b (V-27)

_7b=Vsin/3=v b - Vy, b - Wy, b (V-28)

63
SECTION V - Continued
(V-29)
w--
b= V sin o_ cos _ =w b- VZ, b - WZ, b

where subscripts b; X,b; Y,b; and Z,b denote body-axes components.

Wind Corrections (See Fig. 8)

Geostrophic (due to earth's rotation).- The components of the geostrophic wind

along the X,Y,Z vehicle reference axes are

V x = Re_2 e cos L cos 0g cos _g (V-30)

Vy = -Re_ e cos L (cos _g sin _g - sin Og sin _g cos _g) (V-3i)

V z = Re_2 e cos L(sin _g sin _hg + sin 0g cos qSg cos _hg) (V-32)

Surface winds.- Components of the surface winds along the X,Y,Z vehicle refer-

ence axes are


-- A' B'
W x = W(cos A' sin B' cos 0g cos _g - cos cos cos 0g sin _g

+ sin A' sin 0g) (V-33)

Wy = W_os A' sin B' (sin qSg sin 0g cos _hg- sin _g cos qSg)

- cos A' cos B' (sin _g sin 0g sin ¢>g + cos _g cos @g)

sin A' sin @g cos 0g] (V-34)

WZ= W_os A' sin B' (cos @g cos qSg sin Og + sin _g sin qSg)

- cos A' cos B' (sin @g cos q_g sin 0g- cos @g sin qSg)

sin A' cos _g cos 0g_ (V-35)

Resolution of Engine Thrust and Torque Into Components

Along Vehicle Reference Axes (See Fig. 9)

Vehicle thrust and moment can be resolved into components along the vehicle refer-

ence axes by the equations

°Tu+ 8T (V-36)
T=To+- _- 05Rp M 5RPM

(v-37)
T X = T cos 0rm cos @rm

sin _rm
(v-38)
Ty=T cos 0rm

64
SECTION V - Continued
Tz = -T sin 0rm (V-39)

and, if it is assumed that irm = _rm = 0rm = _rm = 0,

Lrm = Irm_2rm(-q sin erm - r cos 0rm sin _rm) (V-40)

Mrm = Irm_rm(r cos 0rm cos _rm + P sin 0rm ) (V*-41)

Nrm = Irm_2rm(P cos 0rm sin _rm - q cos 0rm cos _rm) (V-42)

Components of Gravitational Acceleration Along X,Y,Z Vehiele Referenee

Axes With Earth-Oblateness Effects Ineluded

If the effect of earth oblateness on the aeeeleration due to gravity is eonsidered

important, terms (2) in the X,Y,Z axes foree equations (V-l) to (V-3) should be replaced,
respectively, by the following X,Y,Z axes weight components:

mg x = m 2 (3 sin L cos L) (5 sin 2 L- 1)cos L

- _5
2
J4 ( sin3
.7 L - 3 sinL) cos L - . . ._cos 0g sin _g

L-1)-2J3(_-_e)3(5 sin3 L- 3 sin L)

- _5
8
J4 (35sin4 L - 30 sin 2 L+ 3)- . (V-43)
• I sin 0g}

mgy 2 2 (3 sin L cos L) +_J3


3 3 (5 sin 2 L - 1) cos
(____e) L

._52 J4l-_ 4
) (7 sin 3 L- 3 sinL) cos L- . . (sin _hg sin 0g sin q_g + cos _Pg cos q_g)

+ _3
j2( )2(3
s,n2 L-i)-2J3(_----e13(5 sin 3 L-3 sin L)

- _58J4 (35 sin 4 (V-44)

65
SECTION V - Concluded

GE ;FT /ae_ 2 3 ae 3 1) cos L


mgz=m _-_2_-)(3 sin L cos L)+ _J3(_)(5 sin 2 L-

-5J4(-_)
cOsL-'2 4(Tsin3L-3sinL) " t (sin_gc°s _gsin0g-cOs _g_in_g)

+ - _" J2 (3 sin 2 L- 1)- 2J 3 (5 sin 3 L- 3 sin L)

sin 4 L - 30 sin 2 L +3) - . (V-45)


• tcos_gC°S 0g_

where

GE = 398 601.2 + 0.4 km3/sec 2 ((14 074 901.1 + 14) X 109 ft3/sec 2)

ae = 6378.160 • 0.005 km (20 925 721.8 + 16 ft)

and J2, J3, and J4 are the second, third, and fourth zonal harmonics. The values
from reference 8 are

J2 = 1082.7(1 + 0.1) × 10 -6

J3 = -2.56(1 + 0.1) × 10 -6

J4 = -1.58(1 ± 0.2) × 10 -6

In equations (V-43) to (V-45), which were derived from the gravitational potential
given in reference 8, oblateness terms through the 4th order are considered and the
earth's longitudinal oblateness is neglected.

Langley Research Center,


National Aeronautics and Space Administration,

Hampton, Va., February 16, 1972.

66
APPENDIX A

SUMMARY OF FREQUENTLY USED FORMS OF AXES TRANSFORMATIONS

AND EQUATIONS OF MOTION

EULER ANGLE TRANSFORMATION BETWEEN

TWO ORTHOGONAL AXES SYSTEMS

For two orthogonal axes systems, an initial X,Y,Z reference system and an

X' ,Y ',Z' system obtained by rotating the initial system through Euler angles _,8,_
(in that order), the transformations between the two systems are given by these
relationships.

Direct Transformation
! t

The transformation from the initial X,Y,Z system to the X' ,Y ,Z system alined

at Euler angles _,O,q_ with the initial system is

Y' = r (A-l)

where

cos e cos cos e sin -sin

sin _ sin O cos _h sin _ sin a sin sin q_ cos a

+ cos _ cos (A-2)


r - sin _ cos

cos _ cos q5 sin O sin _ cos _ sin cos _ cos 0


+ sin _sin q_ - cos _ sin

Inverse Transformation

The transformation from the X ,Y ,_ system back to the X,Y,Z system is

= r-1 y' (A-3)

67
APPENDIX A - Continued
where m

cos e cos sin _) sin 0 cos COS _ cos _ sin 0


- sin _h cos _) + sin _ sin

cos 0 sin sin _ sin 0 sin sin _cos _ sin 0 (A-4)


+ cos cos - cos _ sin

-sin 0 sin q5 cos 0 cos q5 cos 0

TRANSFORMATIONS FOR ACCELEROMETER AND

RATE-GYRO MEASUREMENTS

These equations are simplified forms of the general inverse transformations for
accelerations and angular velocities given in section I; however, they can be used in most
practical applications to correct accelerometer and rate-gyro readings for displacement
and misalinement. For cases that do not fit the assumptions made, the general forms in
section I must be used.

Case I

Orthogonal Instrument Axes; No Restrictions on Misalinement Angles

If X,Y,Z axes accelerometers or rate gyros making flight measurements are


orthogonally alined then, even though the misalinements of the instruments with respect
to the vehicle reference axes are large, their readings can be corrected by using these
equationS.

Acceleration corrections.-

gAx,cg = gAx, i cos O cos _h + gAy,i(sin q5 sin 0 cos _ - sin _ cos _b)

+ gAz,i(cos _ cos _b sin O + sin _ sin _b) + (r 2 + q2)_ x

- (pq - i'):_ x - (pr + _1)_x (A-5)

gAy,cg = gAx, i cos 0 sin _ + gAy,i(sin _ sin 0 sin qS'+ cos @ cos gb)

+ gAz,i(sin _ cos _) sin 0 - cos g_ sin gb) - (pq + i-)_y

+(p2 + r2) y _,(qr - (A-6)

68
APPENDIX A - Continued

gAz,cg = -gAx, i sin 0 + gAy, i sin _ cos 0 + gAz,i cos _b cos 0

- (pr- _t)_ z - (qr + P)Yz + (p2 + q2)_ z (A-7)

Angular -velocity corrections.-

P = Pi cos 0 cos _ + qi(cos _ sin _ sin 0 - sin _ cos _)

+ ri(cos _ cos q5 sin 0 + sin _ sin _b) (A-8)

q = Pi cos 0 sin _ + qi(sin _ sin q5 sin 0 + cos _ cos qS)

+ ri(sin _h cos _b sin 0 - cos _h sin _) (A-9)

(A-10)
r= -Pisin 0+qisin q_ cos 0+r icos 0 cos _b

Case II

Nonorthogonal Instrument Axes; Small Misalinement Angles

If X,Y,Z axes accelerometers or rate gyros are not orthogonally alined but the
misalinements with the vehicle reference axes are small, then these correction equations

can be used. It should be noted that, in equations (A-11) to (A-16), the angles _,0,q5
are in radians.

Acceleration corrections.-

gAx,eg = g(Ax, i - _xAy,i + 0xAz,i ) + (q2 + r2)_ x _ (pq_ r)Yx- (pr + q)ZX (A-11)

gAy,cg = g(,yAx, i + Ay,i_ _byAz,i) _ (pq + _)_y +(p2 + r2)_y _ (qr- f_)_,y (A-12)

(A-13)
gAz,cg = g(-0zAx, i + qSzAy, i + Az,i) - (pr- q):_Z- (qr + P)Yz + (q2 + p2)_ z

Angular-velocity corrections.-

(A-14)
P = Pi - _xqi + 0xri

(A-15)
q = _yPi + qi - _yri

(A-16)
r= -0zp i+ _zqi+ ri

69
APPENDIX A - Continued

SPECIAL FORMS OF TRANSFORMATIONS FOR AERODYNAMIC


FORCE AND MOMENT COEFFICIENTS AND
STABILITY DERIVATIVES

Simplified Forms for Transferring Coefficients and Derivatives


to Another Reference Center

These transformations are the simplified forms obtained from the general equa-
tions of section II by assuming zero angles of attack and sideslip and neglecting aerody _-
namic cross derivatives.
X-axis force coefficients and derivatives.-
!

Cx = C x (A-17)

C_ = CX_ (A-18)

Cx& = Cx& (A-19)

Cx q : Cx q + q-
2_ Cx_ - 4_
7- Cx (A-20)

Y-axis force coefficients and derivatives,-


t

Cy = Cy (A-21)
CT =
yfi Cyfl (A-22)
t

cy_ = Cy_ (A-23)

%-- % *_ % (A-24)

, %
Cyr = _mCyfl+ 49 Cy (A-25)

Z-axis force coefficients and derivatives.-


t

Cz = Cz (A-26)

C_ =
Z_ Czol (A-27)

Cz& = Cz& (A-28)

cz + cz
?cz (A-29)

7O
APPENDIX A - Continued

X-axis moment (roll) coefficients and derivatives.-

(A-30)
C 1 = C 1 + [Cy - Cz

(A-._)
cz_ = ct_ + [ Cy_

' Zc • (A-32)

,
(A-33)
C/p = C/p + _ Cyp + _ C1/_ + C Z _ + _-
2_ 2 Cyfi

, _ 2 _A q + 4__ Cy (A-34)
C/r = C/r + _- CYr - Cl/3 _2 Cy/_ + _ _2

Y-axis moment (pitch) coefficients and derivatives.-

Cm Cm + x7Cz _. (A-35)
TCx

C' = x (A-36)
m_ Cm_ + _- Cz_ - [ Cx_

(A-37)
Cm& = Cm& + [ CZ& - _-CX&

2_ 2
' :_ z CXq + _ Cm_ Z_
Cmq = Cmq + _ CZq- _" -_ C

2_ 4_
4:_ C z + 4g2 Cx (A-38)
Cxo l - -_- Cm -
_2 _2 _2

Z-axis moment (yaw) coefficients and derivatives.-

,
C n = C n +_C z - _-Cy (A-39)

,
(A-40)
Cn/_ = Cni3 - _-CYi3

(A-41)
Cn_ = Cn_ - [ Cy_

Cxa + -_- Cnfi - -- (A-42)


Cnp = Cnp - _- Cyp -

:_ 2_ 2_ 2
C t +4Yc n 4x-_ Cy (A-43)
n r = Cnr - _- CYr - -_- Cn/3 + 7 Cy/_ - _2

71
APPENDIX A - Continued

Transfer of Coefficients and Derivatives From Body


to Wind-Tunnel Stability Axes
These equations are the most frequently used axes transformations. They convert
the coefficients CX, Cy, CZ, etc., measured about body axes in the wind tunnel, into
?
the coefficients about wind-tunnel stability axes CL, Cy, CD, etc. They can be
obtained from the general inverse equations (III-55) to (III-107) of section HI by replac-
ing the angle A by el and letting the angle B equal zero.
)
C D= -C xcos o_ - C z sino_ (A-44)

Cy = Cy (A-45)

CL=C xsino_ -C z cos o_ (A-46)

Cl,wt = C1 cos ol + C n sin ol (A-47)

Cm = Cm (A-48)

Cn,wt = -C l sin _ + C n cos ot (A-49)

C t ____
Dot -CXo l cos ot - CZo t sin _ + CL (A-50)
)
sin
CD& = -Cx& cos _ - Cz& (A-51)

Cyfl = Cy_ (A-52)

= (A-53)
)

CL_ = Cx_ cos _ - CD (A-54)


COS (9/
CL& = CX& (A-55)

Cl_,wt = Cl_ cos ol + Cnt 3 sin ol (A-56)

cl , t : cz cos ol + Cn_ sin ol (A-57)

C/p,wt = Clp 0t (A-58)


cos 2 ot + Cnr sin 2 el + (C/r + Cnp)sin ol cos

C/r,w t = C/r cos 2 ol - Cnp sin o_ cos o_ (A-59)


sin2 _ + (Cnr-C/p)

Cl6,w t = C/6 cos _ + Cn5 sin ol (A-60)

72
APPENDIX A - Continued

Cmol = Cmol (A-61)

Cm& = Cm& (A-62)

Cmq = Cmq (A-63)

Cm6 = Cm 5 (A-64)

(A-65)
Cn/3,wt = -Cl_ sin a + Cnfl cos ol

(A-66)
Cn_,wt = -Cn_ sin ot + Cn_ cos

sin 2 o_ + (Cnr - C/p)sin ot cos (A-67)


Cnp,w t = Cnp cos 2 a - Clr

Cnr,w t = Cnr cos 2 a + C/p (A-68)


sin 2 a - (C/r + Cnp)sin ol cos o_

(A-69)
Cn6,w t = -C/6 sin ot + Cn6 cos

Transfer of Coefficients and Derivatives From

Wind-Tunnel Stability to Body Axes

These transformations convert coefficients about wind-tunnel stability axes CL,


T

Cy, CD, etc., into coefficients about body axes CX, Cy, CZ, etc. They can be
obtained from the general direct transformations in equations (III-1) to (III-36) of sec-
tion III by replacing the angle A by o_ and letting the angle B equal zero.

C x= -C D cos _ +C L sin _ = -C A (A-70)

Cy = Cy (A-71)

C z = -C D sin _ - C L cos cg = -C N (A-72)

C1 = Cl,wt cos ol - Cn,wt sin Ol (A-73)

Cm = Cm (A-74)

C n = Cl,wt sin Ol + Cn,wt cos ot (A-75)

= C' (A-76)
Cxcg - D_ cos a +CLa sina + C N= -CAa
T

CX& = -CD& cos c_ + CL& sin a = -CA& (A-77)

73
APPENDIX A - Continued

= (A-78)

= (A-79)

= C t
Cz_ - D_ sin _ - CL_ cos c_ - C A = -CN_ . (A-80)

___ --C t

Cz& D& sin _ - CL& cos c_ = -CN& (A-81)

sin
Cl_= Cl_,w t cos _ - Cn/_,wt

Cl_ = C 1 /_,wt
. cos _ - Cn_,wt sin (A-83)

C/p = C/p,w t cos 2 _ + Cnr,w t sin2_- (Cnp,wt+C/r,w_Sin


Ol COS Ol (A-84)

C/r = C/r,w t cos2 _ - Cnp,wt sin2 _ + (C/p,wt - Cnr,wt) sin


O/ COS O/ (A-85)

= cos ot sin ol (A-86)


C/5 C/6,wt - Cns,wt

Cm_ = Cm_ (A-87)

Cm&=Cm& (A-88)

Cmq= Cmq (A-89)

Cm5 = Cm 5 (A-90)

Cn/9 = C/_,wt sin o_ + Cnfi,w t cos ol (A-91)

Cn_ = Cl_,w t sin ol + Cn_,w t cos (A-92)

Cnp = Cnp,wt cos 2 _ - C/r,wt sin2 _ + (C/p,w t - Cnr ,wt)sin


O/ COS O_ (A.-93)

Ol COS OZ (A-94)
Cn r = Cnr,wt cos 2 _ + Clp,wt sin 2 _ + (C/r,wt + Cnp ,wt) sin

Cn 6 = C/5,wt (A-95)
sin a + Cns,w t cos

74
APPENDIX A - Continued

TRANSFORMATION EQUATIONS FOR MOMENTS AND PRODUCTS OF INERTIA

Body to Flight Stability Axes

= Ix cos 2 o_0 + I Z sin 2 s 0 - 2Ixz sin s 0 cos .s 0 (A-96)


Ix,s
(A-'97)
Iy,s = Iy

=ixsin 2 s 0+I z cos 2 s 0+2IXZ sin c_0 cos s 0


(A-98)
Iz,s

(A-99)
IXZ,s=(Ix-Iz)sino_0cos °10+Ixz(C°S2 o_0 -sin 2 o_0)

(A-100)
IxY,s = Ixy cos s 0 +Iyz sin s 0

(A-101)
Iyz,s = Iyz cos c_0 - Ixy sin s 0

Body to Principal Axes

(A_I02)
Ix, P =I xcos 2 e +I z sin 2 e - 2Ixz sin e cos e

(A-103)
Iy,p = Iy

(A-104)
iz,p =I z cos 2 e +I xsin 2 e +2Ixz sin e cos e

Flight Stability to Body Axes

(A-105)
IX = iX,s cos 2 c_0 + Iz, s sin 2 s 0 + 2Ixz,s sin c_0 cos c_0

(A-106)
Iy = Iy, s

(A-107)
iz = iX,s sin 2 c_0 + Iz, s cos 2 s 0 - 2Ixz,s sin s 0 cos s 0

IXZ = (Iz,s- Ix,s)sin c_0 c°s _0 + Ixz,s( c°s2 s0-sin2 (_0)


(A-10S)

(A-109)
Iyz = Iyz,s cos s 0 + IxY,s sin c_0

(A-110)
Ixy = IxY,s cos s 0 - Iyz,s sin s 0

Flight Stability to Principal Axes

(A-ill)
Ix,p = Ix, s cos 2 7/+ Iz, s sin 2 _ + 2Ixz,s sin 77 cos 7?

(A-112)
Iy,p = Iy, s

(A-113)
Iz, P = IX, s sin2 _?+ Iz, s cos 2 _?- 2Ixz,s sin _?cos _?

75
APPENDIX A - Continued

Principal to Body Axes

I x=IX, P cos2 e +Iz, P sin 2 e (A-114)

Iy = Iy,p (A-115)

I z=Ix, P sin 2 e +Iz, p cos 2 e (A-116)

IXZ = (Iz, P - Ix,p)sine cos e (A-II7)

Principal to Flight Stability Axes

IX, s = IX, P cos 2 77 + IZ, P sin 2 _? (A-118)

Iy, s = Iy,p (A-119)

IZ, s = IX, P sin 2 _? + IZ, p cos 2 _? (A-120)

Ixz,s =(Ix, P -Iz,p)sin_?cos_/ (A-121)

SPECIAL FORMS OF EQUATIONS OF MOTION

Coupled Linear Equations of Motion

These are the linearized equations that describe the small perturbation motions of
a vehicle about a steady-state flight condition. In steady-state flight, the components of
thrust, aerodynamic force, and other forces in a given direction are balanced by the com-
ponent of vehicle weight and, hence, the initial components do not appear in the equations.
Other specific assumptions made in deriving these equations are

(1) Total-velocity components, angular rates, and angles are equal to a steady-
state value plus a small perturbation value; for example, the velocity com-
ponent along the X-axis is u = U o + Au, where u is the total-velocity
component, Uo is the initial steady-state velodity, and Au is the per-
turbation velocity

(2) Products and squares of perturbations can be neglected

(3) The vehicle has a plane of symmetry (Ixy = Iyz = 0)

(4) Pitch and roll perturbation rates are given in terms of the initial pitch and roll

angles Oo and q_o and not in terms of the total pitch and roll angles
Oo+Ao and 4o+Z_; that is,

76
APPENDIX A - Continued

= Aq COS _bo - Ar sin q5o

and

sin 0o sin 0o
,+ Ar COS qbo
_i = Ap + Aq sin q5o cos OO cos 00

(5) Aerodynamic cross-coupling coefficients are zero

(6) Mass-damping and jet-reaction forces and moments can be neglected

(7) The vehicle is flying at constant thrust (_= 05Rp


Mtcr =0)

These equations are written about the body axes but can be converted to the flight
stability axes by letting Wo = 0 and 0o = Yo" Also, to reduce the number of variables
to six, the incremental velocity components Au, Av, and AW can be expressed in
terms of resultant velocity and incremental angle of attack by using the following equa-

tions (the alternative would be to expand the aerodynamic forces on the right-hand side

in terms of u,v,w derivatives):

For body axis,

AU= AV COS a 0 -V_oAa sin a0

AV = Voo Aa

AW = AV sin s 0 +V_o AS COS S 0

For flight-stability axis,

Au = AV

AV = V_o Aot

AW = Voo A01

The six equations are

A_ +Qo AW + Wo Aq-RoAv- Vo Ar + g(S Aq dt c°s 0°c°s _b°-S Ar dt cos ,0 o sin qSo)

A5
+ _A&_ + CX V -- + CXq CX 5 +...) (A-122)
2m (C X_ As CX& 2Voo AV
Voo A_q_
2V_o, + _

A_r+RoAu+UoAr-PoAW WoAP-g(_Apdtcos 0oC°SqSo+ S Ar dt sin 0o)

2 +
= pV,.___S
2m C Y_ A_ + Cy_ 2Voo + Cyp 2V_ CYr 2V_o
Ar

(A-123)

77
APPENDIX A - Continued

Aw+PoAv+VoAP-QoAu-UoAq+g(_pdtcos Oo sin _o + _ Aq dt sin Oo)

2m
PVo_S(c Za Aa + CZ& 2V_
A&_+ CZv _AV + CZq 2V---_
Aql[ + _ CZ6 A6+ " " .) (A=124)

A_ I x +(Qo z_r + R o hq)(I Z - Iy) - (Po _q +Qo Ap + A_,)Ixz

l_ A_ + q_ 2v_ qp _ qr 2v_ 5 • • (A-125)

Aq iy + (Po Ar + R O Ap)(I x - IZ) + 2(p O Ap - R O Ar)Ixz

= _pVo_S_ Cmo_ _ + Cm V _ + Cm q __.q_!+_2Voo


Cm6 A6+ . . ) (A-126)

A_ Iz + (Po Aq + Qo Ap)(Iy - Ix) + (Qo Ar + R O Aq- A_)Ixz

=lP2 V2St n_ A_ + Cn_ A__+


2Voo Cnp _+
2Voo Cnr _ + Cn6 A5 + " " (A-127)

Uncoupled Equations of Motion

The linearized equations of motion can be uncoupled (lateral motions made inde-
pendent of longitudinal motions) by assuming that the vehicle is in straight and level flight
and that there are no components of initial velocity except Uo, Wo, and Qo in the
initial steady-state condition (i.e., V o = Po = Ro = _o = _bo = 0). Under these assump-
tions, longitudinal equations contain only the variables Au, _0, and Aw (q = 0); and
lateral equations contain only the variables Av, Ar, and Ap.

Longitudinal equations.-

AV
A_ + Qo Aw + W o Aq + gOg cos Og,o + Cxv
As + CX& 2Voo

÷ Cxq 2V----_+
Aq_ _ CX6 A _) (A-128)

78
APPENDIX A - Concluded

2
pv£s / _& _ _v
Avv - Q0 Au - U o hq + g0g sin 0g,o
= _ _Cza aa ÷ Cz& _ + Czv v--_

(A-129)
+ CZq _Aq_ + _ Cz 5 A6 /

A_ Iy P S_ Aa + __ + + + A (A-130)
- ma Cm& 2Voo Cm v_ Cmq _ Cm 5

Lateral equations.-

A;¢ + Uo Ar - W 0 Ap - g(_bg cos 0g,o + _g sin 0g,o )

(A-131)
/

A_ I x + Qo hr(Iz- Iy) - (Qo Ap + A_)Ixz

_ ___._. + Ap
2 (cl/3 ,,/3÷ Cl_ 2Voo C/p + C/r _Ar_ + _ C/5 _61 (A-132)

IZ + Qo Ap(Iy - IX) + (Qo Ar- Ap)Ixz

_ pv2s {c + + Ap
+ Cnr 2V--'-_ Cn5
2 \ n/3 Cn_ 2Voo Cnp 2Vo_ Ar_+_ A6) (A-133)

Wind-Axes Equations for a Point Mass

The wind-axes equations used in ballistic trajectory studies in which the vehicle is
considered to be a point mass are, along the flight path,

m d__V= _pV2S CD _ W sin y + W (A-134)


dt 2
and, normal to the flight path,

mV dy = pV2_____S
CL _ W cos y (A-135)
dt 2

These equations assume a constant thrust acting parallel to the flight path and
neglect jet damping and reaction control forces.

79
APPENDIX B

DERIVATION OF EQUATIONS PRESENTED IN SECTION II

FOR TRANSFER OF COEFFICIENTS TO NEW

REFERENCE CENTER

The aerodynamic forces and moments acting on a flight vehicle can be considered
to be functions of six independent variables: _,/3,V,p,q,r in one system, or u ,v ,w ,p ,q _r
in the alternate system. To evaluate force and moment derivatives at a new reference
center, such as at a new c.g., the derivatives at the new center, which represent changes
in forces or moments with changes in the independent variables, both evaluated at the new
center, are written, for example,

(B-I)

f m f _Y _f m f mT

where the forces X ,Y ,Z and moments L ,M ,N at the new reference center are

related to those at the original reference center by

X =X L =_+_Y -
__, __ __, 9Z ]
(B-2)
Y_, =_Y M M +_Z- _Xl

The variables _,t3,V,p,q,r evaluated at the original reference center can now be
expressed as dependent functions of the independent variables _',/3',V',p',q',r' of the
new system through the equations

8O
APPENDIX B - Continued

= tan. 1 __
w = tan-1 V' sin a' cos _' - py + cl:_
u V' cos c_' cos/3' - q_ + r_

/3 = sin -I _v = sin-I V' sin p' - rx + pz


V V
(B_-3)

V' = Vu2 + v 2 + w 2

:V(V' cos ol' cos fl' -q_. +r_) 2 +(V' sin/3' = rz_+p_)2 +(V' sin or' cos /3' -p:_+q_)2

q q' (B-4)
r r

0 0 0 0 0 0
Since from equations (B-4) --7
8p = ap' 0q' = _q' Or' = 0-'r' and derivatives such as

813 = _ = 0 (they are derivatives of one independent variable with respect to another)
aq, aq,
the derivatives in equations (B-l) can be evaluated from equations (B-3) and (B-4) as

a___.a
=1 _ = 2y_v= ap = __ag_
- = ___r= 0
D_ w D_ t 8_ _ D_ t Da t D_ T

off- 1 as aV _p_ 8cI_ ar


b/3, 7=7=7-7-7 =°
2-Y-
v =i b_2__
= o_E= op = _ = 2s_r
=o
bV' 8V' bV' OV' bV' 8V'

bo_ _ -:_ cos o_ _ = z cos _ + y sin el sin _ o__V = E sin/3 - y sin oz cos /3
ap' V cos/3 ap' V ap'

b__ = 1 _ = b_ff_r
=o (B-5)
ap' _p, op'

- aV = (._ sin oz -.z cos a) cos /3


b.__ =
Oq' :_ cos Volcos/3
+ z sin ot _aq, = (_. cos a xV sin oz) sin fl I_V

8q = I 8p = 8__Er=0
aq' aq, _I'

bet = -_ sin et 0/3 = -x cos _ - y cos ot sin _ b._V =_ cos et cos /3 - z_ sin /3
br' V cos t3 br' V br'

2_.L
=i 2E=___= o
Or' br t br'

81
APPENDIX B - Continued

Substituting equations (B-5) and (B-2) into equations (B-l) and nondimensionalizing yields
the transformation equations given in section IX.
The transformations for the u,v,w derivatives (eqs. (II-41) to (II-77) of Sec. II)
are derived in the same way except that in this system the derivatives are of the form
Jt

aX 0X _u 0X av 0X _ 0X 8p 0X aq_ @X Or
---c=
au 8u au'+ av au'+ _w au'+ _p au'+--_--_u '+ Or au'

o••e•••o•ee•••e•••oe•o•ee•ooo (B-6)

OX oX 8u oX 8v + 0x aw + ox 0p + ox aq_ + OX Or
--F
Or = hi Or' + Ov Or' Ow Or' _p Or' _l Or' Or _r'

The derivatives in equations (B-6) are evaluated from equations expressing velocity com-
ponents at the original reference center in terms of those at the new reference center as

u=u' -q_+r:_
k
v = v' - rx + pz} (B-7)
-I
w=w' -py +qxJ

and are given as

au av _¢¢ ap ___ Or ;
_u' " av' _V' 8p' 8q' Or'

au av 8w
ap, = 0 _p, = _ = -_
(B-8)
au_ _ --_-_=o --_-_=:_
at' aq, at'

_--_--
= :_ --_-_= -_ _----=o
Or' Or' Or'

All other derivatives are zero•

Substituting these derivatives from equations (B-8) along with equations (B-2) into
equations of the form (B-6) and nondimensionalizing then yield the transformations for

u,v,w derivatives given in section II.

In nondimensionalizing, forces are divided by q_S; moments, by qooS_; p, q, r,


_, and _ are multiplied by 1/2V and u,v,w by 1/V. The velocity derivatives, such

as CXv, CYv , and Cmv, are written in terms of static forces and moments by making
use of the following relationship:

82
APPENDIX B - Concluded

= {pV2SCx

therefore,

-- = pVSC
8__XX x + { pV2S aCx
8--_ (B 9)
aV

The derivative _x/aV represents the change in Cx brought about by changes in


velocity V by itself as an independent variable and not by changes in velocity due to
angular-velocity components p,q,r for the rotating body. They are the result of such
'things as aeroelastic and Mach number effects and are usually neglected, so that equa-
tion (B-9) becomes

a_X_X
= pVSC x (B-10)
aV

The relationship in equation (B-10) is also used in nondimensionalizing the derivatives.

The transformations given in section II for the static force and moment coefficients
and the static-stability derivatives (derivatives with respect to _ or _) are simplified
forms that apply only if there is no vehicle rotation and p, q, and r are zero. Although
it probably is not practical to use such forms, more general forms can be derived that
account for significant rotation; for example, it can be shown that the complete form of
a_/8_' from equations (B-3) is

8ol = 1 + iqz - rg) cos _ + (py - q:_) sin


8ol _

and, by similarly evaluating 8_/8_' and 8V/8c_', it can be shown that the more com-

plete transformation for CX_ is

C'=x_ Cx_Ii(q_-ry)
c°s_+(PY-q'q)sin"_+ V cos

- Cx_ _qz - ry)sin o_ V- _ - qg_)cos sin

+ 2CXF(qz - ry)sin ol- (p_- qg_)cos ol]cos ;3 (B-If)


L V J
which reduces to C' = (the form given inSec. II) if p =q= r = 0.
. X_ Cxol

83
APPEND_ C

METHODS OF MEASURING CENTER-OF-GRAVITY LOCATIONS AND

MOMENTS OF INERTIA OF MODELS AND FLIGHT VEHICLES

The methods presented here are a summary of methods used at the Langley
Research Center to measure moments of inertia of rocket-propelled models and the

methods presented in reference 3 for determining c.g. locations and moments of inertia,
of full-scale airplanes. In these methods all oscillations are assumed to be small.

Some other limitations and precautions that must be taken in using these methods, and
which apply to all methods, are discussed in reference 3.

CENTER-OF-GRAVITY LOCATION

Longitudinal c.g. Location

The longitudinal c.g. location can be determined by mounting the aircraft or rocket
on weighing scales. For a typical aircraft installation, two scales measuring the forces

R 1 and R2 are located at wing jack points and a third scale measuring the force R 3
is located at some distance _ forward or aft of the jack points. The longitudinal dis-
tance from the jack points to the c.g. is then given by

= R3_
R1 + R2 + R3 (C-l)

Vertical c.g. Location

In order to determine the vertical c.g. location the vehicle is mounted on a knife

edge as shown in figure 10. The weighing cradle has a weight Wc and a centroid _c,
and the vehicle is supported in various roll and pitch attitudes by a vertical reaction

force R1 acting at a distance Yl" The distance from the knife edge to the airplane c.g.
is obtained from measurements of R1, Yl, and q_ by using the equation

= RlY 1 - Wc_ c sin _b


W sin _ (C-2)

MOMENTS OF INERTIA

Compound-Pendulum Method (Fig. 11)

This method is used mainly to determine pitch and roll moments of inertia on small
models that can easily be mounted with a single attachment point. The moment of inertia
of the model about an axis through the model c.g. is given by

84
APPENDIX C - Continued

.I WM+S_M+S 2 WsZs WMZM 2 (c-3)


= 47r2 PM+S 47r2 PS 2 g

The subscript M+S refers to the model plus the supporting hardware, S to the sup-
from the knife
porting hardware alone, and M to the model alone. The distance zM+S'
edge to the c.g. of the model plus supporting hardware, is obtained from

WS_ s + WM_ M
(C -4)
ZM+ s = WM+S

Spring Method

This method is usually used to measure pitch and roll moments of inertia when the
model is so large that more than one support is needed. The model is mounted on a knife
edge and springs are attached at equal distances on both sides of the knife edge, as shown
in figure 12. The moment of inertia canbe determined from the equation

I = K_-_2 p2 + rngz p2 _ m_2 (c-5)


4_2 4_ 2

Spring Method for Full-Scale Vehicles

This method, which is described in detail in reference 3, uses two sets of springs
with each set having a different spring constant, so that both the moments of inertia and
the vertical c.g. location can be determined. The springs are arranged as shown in fig-
ure 13. The moment of inertia about the point axis through the knife edge is given by

I= (Kt'l - Kt'2)P2 _ roT. 2 _ ic _ male 2 (C-6)

The vertical distance from the knife edge to the c.g. of the model is then determined from

Kt'2 - Kt'l Wjc


= _ (C-7)

\P2)J w
where Wc and _c are the weight and c.g° vertical displacement (measured from the

knife edge) of the weighing cradle that is used. The constants Kt, 1 and Kt, 2 are
determined from the constants for the springs shown in figure 13 by the equations

85
APPENDIX C - Concluded

Kt, 1 = 2Ks,las 2

Kt, 2 = 2Ks,2as 2

Torsion-Pendulum Method

Used for measuring yaw moments of inertia, this method involves mounting the
model on a torsion shaft as shown in figure 14(a). The moment of inertia about the lon-
gitudinal axis of the shaft is given by

I = kp2+S kPS2
(C-8)
4y 2 4_ 2

Multifilar-Pendulum Methods

Essentially the same as the torsion spring method except that the vertical wires or
rods provide the restraining spring moment, an example of a bifilar pendulum is shown in
figure 14(b). The moments of inertia are determined from

p2
I = WM+S M+S 5.2 WS PS 2 5.2 (C-9)
16_2 _' 16_2 _'

where 5. is the diameter of the circle around which the wires or rods are attached and
_' is the length of one of the wires or rods. A system with two wires is called a bifilar

pendulum; a system with three is called a trifilar pendulum. In using these methods, the
wires or rods must be centered about the system mass center.

86
APPENDIX D

DETERMINATION OF LONGITUDINAL AND LATERAL STABILITY

DERIVATIVES BY USING SIMPLIFIED LINEAR ANALYSIS


$

Perhaps the simplest and most direct method of analyzing transient flight data is to
determine static- and dynamic-stability derivatives from the frequency and the time to
damp to half-amplitude of transient motions, such as the damped angle-of-attack oscilla-
tion shown in figure 15. This method has been used effectively to determine the control
response characteristics of rocket-propelled models in reference 9. Even when a more
sophisticated analysis is planned, this method can provide, before the detailed analysis
is begun, a fairly accurate quick-look assessment of the vehicle characteristics from the
basic flight records.

The oscillation in figure 15 represents the transient response to a step-control


input. The method assumes that

(i) The forward velocity and the Mach number of the vehicle are constant

(2) The longitudinal aerodynamic forces vary linearly with angle of attack a, pitch
angle 8, elevator deflection 5e, da/dt, and d_/dt; lateral aerodynamic forces are
assumed to vary linearly with sideslip angle /_, yaw angle _h, rudder deflection 6r,
d/_/dt, and d_/dt

(3) The vehicle is in level flight before the control deflection is applied

(4) Longitudinal and lateral motions are independent of each other (all aerodynamic

cross derivatives such as C/q and Crop are neglected)

Under these assumptions, the longitudinal equations of motion can be simplified to

da (D-I)
mV_° dI_ _-) = CL, o + CLaa + CLSe 5e

Iy d2_ Cm, o + a + Cm_ &6 + Cmq qc + Cm (D-2)


q S_ dt 2 = Cma 2Voo _ 5e 5e

where CL, o and Cm, o are the values of lift and pitching-moment coefficient in
trimmed level flight. The solution to equations (D-l) and (D-2) is of the form

a = Ce a't cos(cot - 4) + a t (D -3)

In equation (D-3), C is a constant that is determined from initial conditions, at is the

trim angle of attack,

87
APPENDIX D - Continued

irCLa - (Cmq CI& _. c .,_ (D -4)


a' = - gLT + ) zvool]

Cma 6 CmqCLa
(D-5)
- I' 2V_ I'm' (a')2

and

m' = --
mVj
qooS
(D-6)
I' = Iy
qooSc

The damping constant a' and the frequency co are determined from measured quanti-
ties; a' can be calculated from the rate of decay of the oscillation (see fig. 15) as

Ao_2

a'= l°ge _ = log eI=._0.693


(D-7)
t2 - t1 tl/2 tl/2

where tl/2 is the time for the oscillation to damp to one-half its initial amplitude. The
frequency co can be calculated from the period of the oscillation P as

2_
(D-8)
P
Equation (D-5) can be rewritten as

CmqCL a
(D-9)
Cma=-I'_a')2+ w2]- _2_oo m t

and equation (D-4), as

4I'V_o/ CLa_ (D-10)


Cmq+Cmh- e ,ka' + 2m'/

The last term of equation (D-9) is usually very small compared to the first term (usually
less than 1 percent) and may be omitted in most cases. (See ref. 9.) Omitting this term
and substituting equations (D-7) and (1)-8) give the equations for the longitudinal static-
and dynamic-stability derivatives in terms of the period and the time to damp to half-
amplitude of the oscillation as

-4_'21y_I + 1-_ _0"693_2_ (D-11)


Cmol = qooS---
_ Lp-"_ 41r2\tl/2 ] ]

88
APPENDIX D - Concluded

(D-12)
(Cm_
+Cm_)--
--_\_ w_c,._
Similar equations for the lateral derivatives can be developed as

(D-13)

-2I z (1.386Voo (D-14)


(Car On_):-__+ _ °_)
For the lateral-derivative equations a small angle of attack is assumed. The values

of Cn_ and (Cnr - Cn_) given by equations (D-13) and (D-14) are those about the stabil-

ity axes.

It should be noted that the mean values of angle of attack or sideslip do not neces-

sarily correspond to the trim value, so that the mean value that is used to determine

tl/2 (eq. (D-7)) must be obtained by selecting positive and negative amplitudes from the
angle-of-attack or sideslip envelopes and determining the mean of these positive and
negative amplitudes.

89
APPENDIX E

USE OF DIRECTION COSINES AND QUATERNIONS

IN MOTION CALCULATIONS

In the equations of motion and throughout this report, vehicle orientation with
respect to given axes is specified in terms of the Euler angles ¢,,0,c>. Two other
methods can be used, however, to specify alinement. One is to give the alinement in ,
terms of direction cosines; the other is to specify alinement in terms of "Euler param_
eters," which are components of a four-parameter quantity called a quaternion. Whereas
there is no particular advantage in one of these methods over the other, both have certain
advantages over the use of the Euler angles. Both eliminate the singularities that occur
when vehicle attitudes approach +90 ° (a condition known as gimbal lock) and, whereas the
equations relating Euler angles to angular rates are nonlinear (eqs. (V-7) to (V-9) of
Sec. V), those relating direction cosines and Euler parameters to angular rates are linear,
so that the computational procedure is simplified, particularly in analog computations.
Either method is preferred to the Euler angle method, therefore, in certain applications.
The basic equations needed in applying these two methods are given here.

DIRECTION-COSINE METHOD

In the direction-cosine method, the axes transformation for any vector has the form

e' =D (E-l)
Z'

where D is the matrix made up of direction cosines and is defined as

_XY

D = _YX _yy fYZ (E -2)

_ZX _ZY -_ZZ

where the elements of D are given in terms of Euler angles as

90
APPENDIX E - Continued

_XX = cos _ cos 0

_XY = sin _ cos 0

_XZ = -sin 0

_YX = -sin g/ cos _ + cos _ sin e sin

= cos _P cos _ + sin _ sin 0 sin (E-3)


_yy

_YZ = cos 0 sin

_ZX = sin _ sin _ + cos _ sin e cos

_ZY = -cos _h sin _ + sin _ sin 0 cos

_ZZ = cos O cos _

It can be shown that D can be determined from the angular velocities by using equa-

tions (E-4) and (E-5) as follows:


B

0 co
Z' -coy,

0 (E -4)
-w Z, coX' D

0
L_wy, - coX,

which is a linear (matrix) differential equation. The angular velocities ¢°X' 'coY' 'coZ'

are those of the pri:ded axes system (see eq. (E-l)) with respect to the unprimed axes
system. In motion calculations where the primed axes are considered as the vehicle
body or other reference axes and the unprimed axes are taken as the gravity-axis sys-

tem (see fig. 5), coX, ,coy, ,W Z, can be determined from the body-axes angular velocities
p,q,r by

-L
coX' !-P-

(E -5)
coy, = q
-D + )cos
L
-(_e + i)sin L
coZ' r
m

or, for a flat nonrotating earth,

coX' = Pl
(E -6)

coY' = _l
coZ'

91
APPENDIX E - Continued

QUATERNION METHOD

Basic Quaternion Relationships


In the quaternion method an axes transformation has the form (ref. 10)

Y' - G
(E-7)

where G is a matrix composed of the Euler parameters e0,el,e2,e 3 and is defined as

(e02 +el2 - e22- e32 ) 2(ele 2 +e0e3) 2(ele 3 - e0e2)

G
2(ele 2 - e0e3) e02 - el2 + e22 - e32 ) 2(e0e I + e2e3) (E -8)

_ 2(ele 3 + e0e2) 2(e2e 3 - e0el) (e02 - el2 _ e22 + e32 )

The Euler parameters are elements of a four-parameter quantity called a quaternion,


defined as

q=e 0 +eli+e2J +e3k (E-9)

where eo,el,e2,e 3 are real numbers and the vectors i,j,k satisfy the following
conditions:
-h
i 2 = j2 = k 2 = -i
jk = -kj = 1_
(E-10)
ij = -ji = k ki -ik

The quantity e0 is the real or scalar part of the quaternion; the terms eli + e2J + e3k
make up the imaginary part. The length or norm of the quaternion is defined as

[ql = _ = +02 + el2 + e22 + e32 (E-11)

The quantity q* is the conjugate of the quaternion and is defined as

q* = e0 - eli - e2J - e3k (E-12)

The transformation in equation (E-7) can also be written in terms of the quaternion and
its Conjugate as

(E-13)

92
APPENDIX E - Concluded

In flight-motion calculations, in which the Euler parameters define the alinement


of the Vehicle body or other reference axes with respect to the gravity axes (see figs. 4
and 5), the differential equation from which the Euler parameters can be determined is

!e
-'0_ --el -e 2
_/.
-e 3
el =i_. (E-14)
e2 2 e3 e0

-e 2
e0 e I
-el_

e0_-
iI r]-G

-
(£e

u__e
+ X)c°s

where G is defined by equation (E-8).

Relationship Between Euler Parameters and Euler Angles

Initial values of e0,el,e2,e 3 for use in performing the calculation indicated by


equation (E-14) can be determined from the initial values of the Euler angles. The Euler
parameters are related to Euler angles by the following equations:

e0 = cos cos _O cos 2_ + sin 2_ sin _O sin 21

O O
eI cos 2_ cos _ sin _ - sin _ sin _ c°s
CE-15)

e 2 = cos sin _ cos + sin cos

e3 sin cos _ cos - cos sin sin

sin 8= -2(ele 3 - e0e2)

2(e0e 1 + e2e3)

tan @=e02 +el 2 _ e22 _ e32 (E-16)

tan _ =
2(o:e2 +e0e3)
e02 - el2 - e22+e

93
APPENDIX F

SCALING PARAMETERS

This appendix summarizes parameters for the three types of scaling used in aero-
dynamic testing: dynamic, aerodynamic, and aeroelastic. These parameters establish
the requirements for a given type of scaling in that each significant parameter must have
the same value for both model and prototype if the simulation is to be valid. In most
cases, it is impossible to satisfy all the scaling requirements at the same time and com_
promise techniques have to be worked out. A knowledge of the scaling parameters is
essential to understanding and applying these techniques.

DYNAMIC-SCALING PARAMETERS

Dynamic-scaiing parameters established the conditions under which the motions


(accelerations, velocities, and angles) and the forces for the model simulate those for
the prototype. These parameters are obtained by requiring that, in the equation describ-
ing the motion of the flight vehicle, the ratio of any one term to another has the same
value for model and prototype. The equation of motion has the general form

Inertial force = Weight + Aerodynamic force (F-l)

Two scaling parameters determined from equation (F-l) are

Froude No. = Vehicle inertial force = mV 2/_ _ V2


Vehicle weight W g_ (F-2)

Mass ratio = Vehicle inertial force = mV2/_ _ W


(F'3)
Aerodynamic force pV2_2 pg_3

In addition, the scaling assumes geometric similarity between model and prototype
as well as similarity in mass distribution. The scale factors given in table VI were
obtained by satisfying all these requirements.

AERODYNAMIC-SCALING PARAMETERS

Aerodynamic-scaling parameters must have the same value for model and proto-
type if the flow field around the model and, hence, the aerodynamic force and moment

coefficients, is to be the same as around the prototype. Aerodynamic-scaling param-


eters determined from the equation of motion of the fluid, in which viscous, pressure,
and gravity forces are considered, are

94
APPENDIX F - Continued

Fluid inertial force PV_ 1_2 pV_


Reynolds No. = (F -4)
Fluid viscous force _Voo_

Much No. = Fluid inertial force = fpV2_2_ 1/2 Voo (F-5)


Fluid pressure force \_/ = _-

Fluid inertial force _ PV 2_2 = _V2 (F-67


Froude No. = Fluid gravity force pg_3 g_

In addition, if surface-tension forces are important, a fourth parameter to be con-


sidered is

Weber No. = Fluid surface-tension force _ _ _ _ (F-7)


Fluid inertial force pV2_2 pV2_

where a is the surface tension per unit length.

The pressure force considered in Much number (eq. (F-5)) is that due to the pres-
sure differential across a shock wave in compressible flow. The pressure-force param-
eter considered in incompressible flow is

Euler No. = Fluid pressure force = _ (F-8)


Fluid inertial force pV2_ 2

The Euler number is usually not important because usually the resultant body forces are

measured; however, it becomes important when body forces are determined from mea-
surements of pressure distribution.

Two scaling parameters determined from the laws of thermodynamics are

Prandtl No. = 1 X Heat added by convection


Reynolds No. Heat added by conduction

v cppV_T/_ = _CPP (F-9)


: V_---_ kT/_2 k

where v = y/p, k is the coefficient of heat conduction (thermal conductivity) of the

fluid, and

Grashof No. - 3g AT (F-10)


2
V_Wo

where AT is the temperature difference between two representative points in the fluid •
and To is a representative temperature.

95
APPENDIX F - Continued
Aeroelastic-Scaling Parameters
The basic aeroelastic-scaling parameters to be satisfied in simulating deforma-
tions caused by aerodynamic loads are
For structural elongation:
EA _ (F-11)
pV2_ 2
For structural bending:
EI
(F-12)
pV2_ 2

in whieh it is assumed that the strains in the model structure are the same as for the

prototype; that is,

ern = 1.0 (F-13)


ep

It is also assumed that neither model nor prototype materials are stressed beyond their
elastic limits and that the ratios of shear-to-tensile stress and of shear-to-elastic

modulus are the same for model and prototype. A more thorough discussion of elastic
scaling is given in reference 11.

COMMENTS ON SCALING PROCEDURES

Dynamic Scaling

In general, tests of dynamically scaled models are limited to subsonic speeds


(below M = 0,6) because of the requirement that the force coefficients must be the same
on model and prototype. For vehicles designed to fly at transonic or supersonic speeds,
the model is usually tested at the correct full-scale Mach number. The motions, loads,
and accelerations of the test models are then reduced to nondimensional coefficients by
utilizing the known mass and inertial characteristics of the models, and these coeffi-
cients are in turn used to calculate the motions, loads, and accelerations of the prototype.
In any case, care must be taken to insure that Reynolds number effects do not introduce
an unacceptable distortion. Reynolds number effects can be neglected in most tests if
the model Reynolds number is above 106; however, in general, the Reynolds number
effects on skin-friction drag must be accounted for.

Aerodynamic Scaling

In wind-tunnel and free-flight model testing, it is impossible to satisfy all the


aerodynamic-scaling parameters at the same time; however, the problem usually reduces

96
APPENDIX E- Concluded

to one in which only two of the parameters, Reynolds number and Mach number, are
important. In most problems, the effect of the fluid gravity on the motion is negligible,
and so the Froude number and the Grashof number of the fluid can be neglected. It can
be shown that the Prandtl number is the same for model and prototype if the ratio of

specific heats Cp//C v is the same for both, and this condition is automatically satisfied
if air is used as the fluid in the model tests. A still further simplification is possible
because Reynolds number effects are limited to the narrow boundary-layer region at the
body surface and the flow outside the boundary layer is practically without the influence

of viscosity. Thus, Reynolds number and Mach number effects can be considered sepa-
rately. Most wind-tunnel tests at transonic and supersonic speeds are made at the full-
scale Mach number. The Reynolds number effects, called scale effects, are accounted

for by making corrections to skin-friction drag coefficient (the coefficient primarily


affected by Reynolds number) or by using boundary-layer trips, wires or a rough grit, to
cause separation at the proper chordwise position and therefore to simulate the full-scale
flow pattern. These techniques are discussed in textbooks on wind-tunnel testing (for
example, ref. 12) and reports dealing with their application, such as reference 13.

97
REFERENCES

1. Breuhaus, Waldemar O.: Resum6 of the Time Vector Method as a Means for Analyz-
ing Aircraft Stability Problems. WADC Tech. Rep. 25-299, U.S. Air Force, Nov.
1952.

2. Taylor, Lawrence W., Jr.; Iliff, Kenneth W.; and Powers, Bruce G.: A Comparison of
Newton-Raphson and Other Methods for Determining Stability Derivatives From
Flight Data. AIAA Paper No. 69-315, Mar. 1969.

3. Wolowicz, Chester H.: Considerations in the Determination of Stability and Control


Derivatives andDynamic Characteristics From Flight Data. AGARD Rep. 549-Pt. I,
1966.

4. Thelander, J.A.: Aircraft Motion Analysis. FDL-TDR-64-70, U.S. Air Force, Mar.
1965.

5. Gracey, William: Summary of Methods of Measuring Angle of Attack on Aircraft.


NACA TN 4351, 1958.

6. Larson, Duane B.: Identification of Parameters by the Method of Quasilinearization.


CAL Rep. No. 164, Cornell Aeronautical Lab., Inc., May 14, 1968.

7. Alien, Willie L.; and Weight, Robert H.: The Use of Flightpath Accelerometers in
Performance Flight Testing. Proceedings of the 13th Annual Air Force Science
& Engineering Symposium, Vol. II, U.S. Air Force, Sept. 1966.

8. Melbourne, William G.; Mulholland, J. Derral; Sjogren, William L.; and Sturms,
Francis M., Jr.: Constants and Related Information for Astrodynamic Calcula-
tions, 1968. Tech. Rep. No. 32-1306 (Contract NAS 7-100), Jet Propulsion Lab.,
California Inst. Technol., July 15, 1968.

9. Gillis, Clarence L.; and Mitchell, Jesse L.: Determination of Longitudinal Stability
and Control Characteristics From Free-Flight Model Tests With Results at Tran-
sonic Speeds for Three Airplane Configurations. NACA Rep. 1337, 1957.

10. Wagner, William E.; and Serold, Arno C.: Formulation on Statistical Trajectory
Estimation Programs. NASA CR-1482, 1970.

11. Murphy, Glenn: Similitude in Engineering. Ronald Press Co., 1950.

12. Pankhurst, R. C.; and Holder, D.W.: Wind-Tunnel Technique. Sir Isaac Pitman
& Sons, Ltd. (London), 1965.

13. Braslow, Albert L.; Hicks, Raymond M.; and Harris, Roy V., Jr.: Use of Grit-Type
Boundary-Layer-Transition Trips on Wind-Tunnel Models. NASA TN D-3579, 1966.

14. Mechtly, E.A.: The International System of Units - Physical Constants and Conver-
sion Factors (Revised). NASA SP-7012, 1969.

98
BIBLIOGRAPHY ON METHODS OF ANALYZING FLIGHT DATA

Breuhaus, Waldemar O.: Resum_ of the Time Vector Method as a Means for Analyzing
Aircraft Stability Problems. WADC Tech. Rep. 25-299, U.S. Air Force, Nov. 1952.
$

Breuhaus, W. O.; and Segel, L.: Notes on Aircraft Dynamic Stability Testing and Analysis
Techniques. Rep. No. CAL-51, Cornell Aeronautical Lab., Inc., Apr. 9, 1953.

Charters, A. C.: The Linearized Equations of Motion Underlying the Dynamic Stability
of Aircraft, Spinning Projectiles, and Symmetrical Missiles. NACA TN 3350, 1955.

Denery, Dallas G.: An Identification Algorithm Which is Insensitive to Initial Parameter


Estimates. AIAA Paper No. 70-34, Jan. 1970.

Doetsch, K.H.: The Time Vector Method for Stability Investigations. R. & M. No. 2945,
Brit. A.R.C., 1957.

Dolbin, Benjamin H., Jr.: A Differential Correction Method for the Identification of Air-
plane Parameters From Flight Test Data. Proceedings of the National Electronics
Conference, Vol. XXV, 1969, pp. 90-94.

Donegan, James J.: Matrix Methods for Determining the Longitudinal-Stability Deriva-
tives of an Airplane From Transient Flight Data. NACA Rep. 1169, 1954.

(Supersedes NACA TN 2902.)

Donegan, James J.; and Pearson, Henry A.: Matrix Method of Determining the
Longitudinal-Stability Coefficients and Frequency Response of an Aircraft From
Transient Flight Data. NACA Rep. 1070, 1952. (Supersedes NACA TN 2370.)

Donegan, James J.; Robinson, Samuel W., Jr.; and Gates, Ordway B., Jr.: Determination
of Lateral-Stability Derivatives and Transfer-Function Coefficients From Frequency-
Response Data for Lateral Motions. NACA Rep. 1225, 1955. (Supersedes NACA
TN 3083.)

Duncan, W.J.: The Principles of the Control and Stability of Aircraft. Cambridge Univ.
Press, 1952.

Eggleston, John M.; and Mathews, Charles W.: Application of Several Methods for Deter-
mining Transfer Functions and Frequency Response of Aircraft From Flight Data.
NACA Rep. 1204, 1954. (Supersedes NACA TN 2997.)

Etkin, Bernard: Dynamics of Flight. John Wiley & Sons, Inc., c.1959.

Gates, S. R.; and Lyon, H. M.: A Continuation of Longitudinal Stability and Control
Analysis. Pt. II. Interpretation of Flight Tests. R. & M. No. 2028, Brit. A.R.C.,
Aug. 1944.

99
Greenberg , Harry." A Survey of Methods for Determining Stability Parameters of an
Airplane From Dynamic Flight Measurements. NASA TN 2340, 1951.

Larson, Duane B.: Identification of Parameters by the Method of Quasilinearization.


CAL Rep. No. 164, Cornell Aeronautical Lab., Inc., May 14, 1968.

Lasdon, L. S.; Mitter, S. K.; and Waren, A.D.: The Conjugate Gradient Method for
Optimal Control Problems. IEEE Trans. Automat. Contr., vol. AC-12, no. 2, Apr.
1967, pp. 132-138.

Milliken, William F., Jr.: Progress in Dynamic Stability and Control Research.
J. Aeronaut. Sci., vol. 14, no. 9, Sept. 1947, pp. 493-519.

Mueller, R.K.: The Graphical Solution of Stability Problems. J. Aeronaut. Sci., vol. 4,
no. 8, June 1937, pp. 324-331.

Neumark, S.: Analysis of Short-Period Longitudinal Oscillations of an Aircraft -


Interpretation of Flight Tests. R. & M. No. 2940, Brit. A.R.C., 1956.

Perkins, Courtland D.; and Hage, Robert E.: Airplane Performance Stability and Control.
John Wiley & Sons, Inc., c.1949.

Rampy, John; and Berry, Donald T.: Determination of Stability Derivatives From Flight
Test Data by Means of High Speed Repetitive Operation Analog Matching.
FTC-TDR-64-8, U.S. Air Force, May 1964. (Available from DDC as AD 440 785.)

Schumacher, Lloyd E.: A Method for Evaluating Aircraft Stability Parameters From
Flight Test Data. USAF Tech. Rep. No. WADC-TR-52-71, U.S. Air Force, June
1952.

Shinbrot, Marvin: A Description and a Comparison of Certain Nonlinear Curve-Fitting


Techniques, With Applications to Analysis of Transient-Response Data. NACA
TN 2622, 1952.

Shinbrot, Marvin: A Least Squares Curve Fitting Method With Applications to the Cal-
culation of Stability Coefficients From Transient-Response Data. NACA TN 2341,
1951.

Shinbrot, Marvin: On the Analysis of Linear and Nonlinear Dynamical Systems From
Transient-Response Data. NACA TN 3288, 1954.

Swanson, Robert S.; and Mastrocola, N.: Survey of Methods for Determining Stability
Derivatives in Free Flight. U.S. Naval Air Missile Test Center, Aug. 2, 1948. •

Wolowicz, Chester H.: Considerations in the Determination of Stability and Control


Derivatives and Dynamic Characteristics From Flight Data. AGARD Rep. 549-Pt. I,
1966.

i00
Wolowicz, Chester H.: Time-Vector Determined Lateral Derivatives of a Swept-Wing
Fighter-Type Airplane With Three Different Vertical Tails at Mach Numbers
Between 0.70 and 1.48. NACA RM H56C20, 1956.

Young, Peter C.: Process Parameter Estimation and Self Adaptive Control. Theory of
Self-Adaptive Control Systems, P. H. Hammond_ ed., Plenum Press, 1966,
pp. 118-140.

Zbrozek, J.K.: On the Extraction of Stability Derivatives From Full-Scale Flight Data.
AGARD Rep. 190, Apr. 1958.

i01
TABLE I.- CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS

Multiply by to obtain
Physical quantity U.S. Customary Unit conversion factor SI Unit
(a) (b)
Angle deg 0.01745329 rad

Length ft 0.3048 m

Velocity ft/sec 0.3048 m/sec


knot (international) 0.5144444 m/sec
knot (U.S. statute) 0.44704 m/sec
mph 0.44704 m/sec

Acceleration ft/sec 2 0.3048 m/sec 2

Pressure lbf/ft2 47.88026 N/m2


atmosphere 1.01325 × 105 N/m 2
ft of H20 (39.2 ° F) 2.98898 × 103 N/m 2
in. of Hg (60 ° F) 3.37685 × 103 N/m 2

Moment of inertia slug-ft 2 1.355818 kg-m 2

Mass slugs 14.59390 kg


Force lbf 4.448222

Moment ft-lbf 1.355818 N-m

Area ft 2 0.09290304 m 2

Mass flow slugs/sec 14.59390 kg/sec


Density slugs/ft3 515.3788 kg/m 3

Spring constant lbf/in. 175.1268 N/m


lbf/ft 14.59390 N/m
Gravitational constant ft3/sec 2 0.02831685 m3/sec 2

Specific heat Btu(thermochemical) 4.184 × 103 joule/kg-°C


Ibm -OF

Kinematic viscosity ft2/sec 0.09290304 m2/sec

Coefficient of Btu(thermochemical) -in. 518.87315 J/m-sec-K


ft2-sec-oF
heat conduction

Temperature OF (c)
oc (c)
OR (c) K

Coefficient of slugs/ft-sec 47.880258 N-sec/m2


viscosity

abased on values in ref 14.


bprefixes to indicate multiples of SI Units are as follows:

Prefix (symbol) Multiple

kilo(k) 103
hecto (h) 102
deka (da) i0
deci (d) i0-I

centi (c) 10-2


milli (m) i0-3
micro (_) i0 -6

c Temperatures related by following formulas:

K=_ 5 (o F + 459.67)
K = oc + 273.15

f_=5°R
9

102
TABLE II.- RELATIONSHIPS BETWEEN o_,_,V AND u,v,w DERIVATIVES

CXa = -CXu sin a cos _ + CXw cos _ cos

sin a sin
CX_ = -CXu cos _ sin fi + CXv cos _ - CXw

CXv = CXu cos a cos _ + CXv sin _ + CXw sin a cos

sin
CXu = -Cxa cos Cx_ cos a sin _ + CXv cos a cos

CX v = Cx_ cos _ + CXv sin

sin a cos
CX w = CX v cos
COS _a Cx_ sin a sin _ + CXv

103
TABLE III.- ANGLE DESIGNATIONS FOR DIRECT a TRANSFORMATIONS

l
Angle designation
Axes conversion I

From To A B
(Unprimed coefficients) (Primed coefficients)
b 0
Principal Body 6

Flight stability Body aO 0

Flight stability Principal Nb 0

Wind-tunnel stability Body 01 0

Wind-tunnel stability Principal (3/ - 6 0

Wind-tunnel stability Flight stability a - a0 0


Wind Body 13/

Wind Principal OZ - 6

Wind Flight stability a - a0


Wind Wind-tunnel stability 0

a Direet transformations (eqs. (IN-I) to (III-54)) represent rotation from ref-


erence (unprimed) axes system through angle -B about Z-axis and then through
angle +A about Y-axis to new (primed) axes system. Inverse transformations
(table IV and eqs. (III-55) to (III-107)) represent reverse of direct.

bAngles e and _/ between the principal axes and the body and flight stabil-

ity axes, respectively, can be determined from moments of inertia by

tan 2e = 2Ixz
I z - IX

2Ixz
tan 2_/= ,s
IX, s - Iz, s

104
TABLE IV.- ANGLE DESIGNATIONS FOR INVERSE a TRANSFORMATIONS

Axes conversion Angle designation

From To A B
(P rimed coefficients) (Unprimed coefficients)

C 0
Body Principal
0
Body Flight stability a0
Wind-tunnel stability 01 0
Body
Body Wind
nb 0
Principal Flight stability
o_ - e b 0
Principal Wind-tunnel stability
Wind 0l - E
Principal
0
Flight stability Wind-tunnel stability a = a0

Flight stability Wind - a0


Wind 0
Wind-tunnel stability

ainverse transformations (eqs. (III-55) to (III-107)) are reverse of direct

(table III and eqs. (III=l) to (III-54)) and represent rotation of primed axes system,
through angle -A about Y-axis and then through angle +B about Z =axis, until
axes coincide with original unprimed system.

b Equations for determining angles e and _? given in footnote in table III.

105
TABLE V.- DESIGNATIONS OF FORCE AND MOMENT COEFFICIENTS

FOR DIFFERENT AXES SYSTEMS

Coefficients for axes system -

Component
Body or Flight Wind -tunnel
Wind
principal stability stability
(a)
t
X-axis force Cx or -C A Cx,s -C D -C D
Y-axis force Cy Cy,s Cy CC
Z-axis force CZ or -C N -C L -C L
CZ,s
X-axis moment (roll) Cl C1 ,s Cl ,wt C1 ,w
Y-axis moment (pitch) Cm Cm,s Cm ,wt Cm ,w
Z-axis moment (yaw) Cn Cn,s Cn,wt Cn,w

a Subscript P sometimes used to denote principal axes.

106
TABLE VI.- SCALE FACTORS a FOR DYNAMIC SCALING

Area b .................................... R/2


Volume ................................... Rl 3
Velocity ............................ Rg1/2 R/1/2
Acceleration ................................. Rg
Mass • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • • •
Rp R/3
Mass-flow ....................... . . . . . . . . . . . RpRg 1/2 R/5/2

Weight ................................... Rp Rg R/3

Force coefficient .............................. 1o0

Moment coefficient ............................. 1.0

Wing loading ................................. Rp Rg R l

Time ..................................... R/1/2 Rg -1/2

Dynamic pressure .............................. Rp Rg R 1

Force .................................... Rp Rg R/3

Angular velocity• . . ............... . o . . . . . . . . . . . RI-1/2 Rgl/2

Angular acceleration ............................ R/-1 Rg

Moment ................................... Rp Rg R/4

Moment of inertia ............... . . . . . . . . . . . . . . . RpR/5

a scale factor is ratio of model quantity to prototype quantity; for example,

Model area = R/2 × Prototype area

b Definitions of symbols:

Rl = Model length
Prototype length

Rg = Acceleration of gravity at model altitude


Acceleration of gravity at prototype altitude

Rp = Air density at model altitude


Air density at prototype altitude

107
Vehicle moment reference center

P Xi
Instrument center

Y
q
r

Yi

Zi

Figure 1.- Systems of vehicle reference axes and instrument axes.

108
Center of gravity

JS- x,
y,_.r _.New reference center to which
data are to be referred

Z_

Figure 2.- Axes systems for transfer from vehicle c.g. to new reference center
by equations of section If.

Yw

Yb ,YP, Ys, Yw t

Figure 3.- Systems of vehicle reference axes, including body, principal, wind,
flight stability, and wind-tunnel stability.

109
X

/ Xg

\
Xe
Z

Z e
Ye

Figure 4.- Relationship between earth-centered inertial axes, gravity axes,


and vehicle reference axes.

110
. , -. • .

X
Vertical velocity of vehicle

Tangent to circle of local latitude Resultant velocity of vehicle V

Plane of local horizontal Horizontal velocity of vehicle

nocrceo
ocoue
y Xe

Ye

Ze

Figure 6.- Relation of heading angle H and flight-path angle y to


earth-centered inertial axes.

112
Xb

Figure 7.- Resolution of relative velocity into components along vehicle body axes.

113
Tangent to circle of local latitude _--_i /

ni!!i to _ \_ S Resultant surface-wind velocity W

_ __ Xe

Ze

Figure 8.- Directions of surface and geostrophic winds.


Y

Figure 9.- Alinement with respect to vehicle reference axes of thrust and torque
dueto rotating mass of engine.

(_ Cradle

RI

Pivot point

Figure 10.- Determination of vertical location of center of gravity.

115
Pivot point

Figure 11.- Measurement of moment of inertia by compound-pendulum method.

k
>
_ Pivot
point
/// /t/// / _'1 II /; ///

..--[--.. ..__ _ _._..

Figure 12.- Measurement of moment of inertia by spring method.

116
Figure 13.- Determination of vertical location of center of gravity and rolling moments of inertia
for full-scale airplanes. (Reproduced from ref. 3.)
Iii111
!
_ Vehicle c.g.

(a) Torsion pendulum.

//i i//l i///


i i

_ _'__

(b) Bifilar pendulum.

Figure 14.- Methods of measuring yawing moments of inertia.

118
a

i ............... Trim angle of


attack

Time ,seconds ._

Figure 15.- Damped angle-of:attack oscillation assumed in analysis of appendix D.

NASA-Langley, 1972 _ 1 L-7201 119


NATIONAL AERONAUTICS AND SPACE ADMISTRATION
WASHINGTON, D.C. 20546
POSTAGE AND FEES PAID
NATIONAL AERONAUTICS AND
OFFICIAL BUSINESS
FIIST OLA{III MAIL SPACE ADMINISTRATION
PENALTY FOR PRIVATE USE S300

If Undeliverable (Section 158


POSTMASTER : Postal Manual ) Do Not Return

"The aeronae/tical and space activities of the Unitect States shall be


conducted so as to contribute . . . to the expansion of human knowl-
edge of phenomena in the atmosphere and space. The Administration
shall provide for the widest practicable and appropriate dissemination
of information concerning its activities and the results thereof."

--NATIONAL AERONAUTICS AND SPACE ACT OF 1958

NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS

TECHNICAL REPORTS: Scientific and TECHNICAL TRANSLATIONS: Information


technical information considered important, published in a foreign language considered
complete, and a lasting contribution to existing to merit NASA distribution in English.
knowledge.
SPECIAL PUBLICATIONS: Information
TECHNICAL NOTES: Information less broad derived from or of value to NASA activities.
in scope but nevertheless of importance as a Publications include conference proceedings,
contribution to existing knowledge. monographs, data compilations, handbooks,
sourcebooks, and special bibliographies.
TECHNICAL MEMORANDUMS:
Information receiving limited distribution TECHNOLOGY UTILIZATION
because of preliminary data, security classifica- PUBLICATIONS: Information on technology
tion, or other reasons. used by NASA that may be of particular
interest in commercial and other non-aerospace
CONTRACTOR REPORTS: Scientific and
applications. Publications include Tech Briefso
technical information generated under a NASA Technology Utilization Reports and
contract or grant and considered an important
Technology Surveys.
contribution to existing knowledge.

Details on the availability of these publications may be obtained from:

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION


W_shinl_on, D.C. 20546

You might also like