Math & Physics For Flight Testers, Volume I Axis Transformations
Math & Physics For Flight Testers, Volume I Axis Transformations
Table of Contents
Introduction........................................................................................................................................7.1
Coordinate Systems............................................................................................................................7.1
Inertial Coordinate System.............................................................................................................7.1
Earth Axis System......................................................................................................................7.2
Vehicle Axis Systems.....................................................................................................................7.2
Stability Axes.............................................................................................................................7.3
Wind Axes..................................................................................................................................7.4
Principle Axes............................................................................................................................7.5
Euler Angles.......................................................................................................................................7.7
Force and Velocity Transformations...............................................................................................7.8
Angular Rate Transformations......................................................................................................7.11
Flight Path Angles............................................................................................................................7.15
Transforming Motions Through Axes..............................................................................................7.16
Velocity Transformation...............................................................................................................7.16
Acceleration Transformation........................................................................................................7.17
Specific Angular Momentum.......................................................................................................7.18
Summary..........................................................................................................................................7.20
References........................................................................................................................................7.20
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
Introduction
Various axis and coordinate systems have been developed for specific uses. This chapter first
introduces these systems, then develops the process for transforming rates and accelerations from one
system to another. Understanding this process is required for developing theory for vehicle equations of
motion, simulator and inertial navigation system programming. Because transformations are useful in so
many fields, they are treated as a separate chapter in this text. To follow this chapter thoroughly, a good
understanding of vector analysis and basic matrix algebra including multiplication and inverse operations
is a prerequisite.
Coordinate Systems
There are two primary coordinate systems that are useful in the analysis of vehicle motion; the inertial
and the vehicle coordinate systems. Newton's laws apply only when observed from inertial space, but
practical instrumentation is strapped to the test vehicle. Before attempting to develop transformation
equations, it is first appropriate to define the various inertial and vehicle systems. According to
convention, all coordinate systems used will be right-hand orthogonal.
7.1 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
Earth
Sun
Body Axes
7.2 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
The body axis system is the most straightforward. The axes originate at the aircraft cg and extend
along the body as shown in Figure 7.3. The positive x-axis points forward along an aircraft horizontal
reference line with the positive y-axis out the right wing. The positive z-axis points downward out the
bottom of the aircraft, the xz-plane is usually the vehicle plane of symmetry. Note that lower case xyz are
used. The unit vectors are iˆ , ĵ , and k̂ also have origins at the aircraft cg.
u
iˆ
xis
xa p
ya
xis q
ĵ
z axis
v
r
k̂
w
Figure 7.3 Body Axis Systems
Total aircraft linear velocity can be broken down into the three orthogonal components along the body
axis system:
u = forward velocity, along the positive x-axis
v = side velocity, along the positive y-axis
w = vertical velocity, along the positive z-axis
Similarly, the total aircraft rotational rates can be broken down into the three orthogonal components
about the body axes:
p = roll rate about the x-axis (positive for right roll)
q = pitch rate about the y-axis (positive for pitch up)
r = yaw rate about the z-axis (positive for yaw right)
Application of the body axis system ensures that the moments and products of inertia are constant
(assuming constant mass distribution) and that aerodynamic forces and moments depend only upon the
relative velocity orientation angles and . The body axis system is also the natural frame of reference
for most airframe-mounted instrumentation.
Stability Axes
Assuming zero sideslip, the stability axis system is just like the body system except that the xs axis
points into the relative wind instead of along the nose. This can be easily visualized by rotating the body
system about the y axis by degrees. Obviously then, the difference between xb and xs is , the
difference between zb and zs is also , and the yb and ys axes are coincident. The used in this operation
is the equilibrium flight (trim) value. This initial realignment does not alter the body-fixed nature of the
axis system. In other words, once established, the stability axis system remains fixed to the body for that
7.3 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
application. Of course, each application may have a different value for trim and therefore a different
moment of inertia and product of inertia along the axes.
xb
xs RW
ys, yb zb
zs
Figure 7.4 Stability Axis System
From the geometry of this figure, any force, velocity or acceleration along any stability axis can be
transformed to the body axes as follows:
xb = xs cos zs sin zb = zs cos xs sin yb = ys
Wind Axes
Wind axes are oriented with respect to the flight path of the vehicle, i.e., with respect to the relative
wind VT. If the reference flight condition is symmetric, then sideslip is zero and the wind axes coincide
with the stability axes. Just as the difference between xb and xs illustrates , the difference between xs and
xw illustrates , Figure 7.5a. The equations with the figure show how to transform velocity, force or
acceleration along any wind axis to the stability axes. The relationship of true velocity and its
components to and and the body axis coordinate system is shown in Figure 7.5b.
xb
xs
xw xs = xw cos yw sin xs
ys = yw cos + xw sin RW
zs = zw
u xw
VT
ys
v yb
yw w
(a) (b)
zb
Figure 7.5 Velocity Components and the Aerodynamic Orientation Angles and
The complete transformation [of forces, velocity, or accelerations] from the wind axis system to the
body system is simply the combination of the previous two transforms. This shown below in matrix
form.
7.4 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
7.5 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
7.6 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
Using Figure 7.5, angles and can be expressed in terms of the velocity components:
w
sin
VT cos
w
ASSUMING is small, then cos [see Table 1.1] and
VT
w
ASSUMING is also small, then sin [units in radians, see Table 1.1] and
VT
v
For angle of sideslip: sin
VT
v
If is small, then sin and:
VT
cos sin
(deg) (rad)
0° 1 0 0
5° 0.9962 0.0872 0.0872
10° 0.9848 0.1736 0.1745
15° 0.9659 0.2588 0.2618
20° 0.9396 0.3420 0.3491
25° 0.9063 0.4226 0.4363
30° 0.8660 0.5000 0.5236
Table 7.1 Trigonometric Values for Small Angles
Principle Axes
Principle axes are the natural axes of rotation of the aircraft when only the mass properties are
considered and aerodynamic effects are neglected. The orientation of these axes relative to the aircraft is
a function of the mass distribution and are those axes where all of the products of inertia are reduced to
zero. This can be seen as aligning the axes with "dumbbells" that duplicate the aircraft's mass distribution
along each axis (Figure 7.6).
x y
z z
7.7 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
These are measures of rotational inertia about their respective axes. The integrations performed in
Equation 7.4 can be visualized using Figure 7.7 where the square of the distance for each point of mass is
added up.
d d2 = x2 + z2
x
m1
z
Figure 7.7 Moment of Inertia, Iy
Considering the layout of any aircraft, Iz is always the largest value. Iy > Ix for fuselage-loaded aircraft
and Ix > Iy for wing-loaded aircraft.
The products of inertia are defined as:
Ixy = Iyx xydm (7.6)
Products of inertia are measures of asymmetry. Figure 7.8a shows an object with its mass distributed
about the x and y "body" axes. Integrating this mass according to Equation 7.6 essentially means
concentrating the mass from opposite quadrants into the appropriate size "dumbbells", one positive and
one negative as shown in Figure 7.8b. The principle axis is essentially the average of these dumbbells.
When they have different tilts or sizes due to asymmetric mass distribution, the principle (average) axis
lies along some line different from the body axis. The value for Ixy reflects both the magnitude and the tilt
of this misalignment. When the mass distribution is symmetric about some line (as shown in Figure 7.6a
and 7.6c), the dumbells are symmetrically tilted and massed. In this case, the average weight lies in a line
that is coincident with a body axis. The product of inertia is therefore zero for views having a plane of
symmetry.
7.8 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
x
x Ixy = 0
Principle Axis
y
y
m1
Ixy = (+)
(a) (b)
Figure 7.8 Visualizing the Product of Inertia
The inclination of the xp axis relative to the xb axis has a direct bearing on the inertial moments
experienced about the body axes as reflected by the product of inertia term Ixz in the equations of motion.
The equations of motion would be simplified if these axes were used, but it is difficult to accurately
describe the aircraft motion in this system. Principle axes are not generally used in the basic analysis of
the motion of an aircraft. They are, however, used in more advanced studies such as roll coupling and
spins.
Euler Angles
The orientation of any coordinate system relative to another can be given by three "Euler" angles,
which are consecutive rotations about the z, y, and x axes. They carry one frame into coincidence with
another. In flight dynamics, the Euler angles used are those, which rotate the earth axis system into
coincidence with the relevant vehicle axis system (Figure 7.9).
x
y
z
Earth Plane
7.9 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
- Yaw Angle: The angle between the projection of the vehicle xb-axis onto the horizontal reference
plane and the initial reference position of the earth x-axis. Note that yaw angle is the vehicle
heading only if the initial reference is north.
- Pitch Angle: The angle measured in the vertical plane between the vehicle xb-axis and the
horizontal reference plane.
- Roll Angle: The angle, measured in the yz-plane of the body-axis system, between the y-axis and
the horizontal reference plane. This is the same as bank angle and is a measure of the rotation
(about the x-axis) to return the aircraft to a wings level condition.
The importance of the sequence of the Euler angle rotations cannot be overemphasized. Finite
angular displacements do not behave as vectors. Therefore, if the sequence is performed in a different
order than ,,, the final result will be different. This fact is clearly illustrated by the final aircraft
attitudes shown in Figure 7.10 in which two rotations of equal magnitude are performed about the x and y
axes, but, in opposite order. Addition of a rotation about a third axis does nothing to improve the
outcome.
y
Rotation Sequence 1
Rotate -90° Rotate +90°
about x axis about y axis
x x x
y z z
z
x x
Rotation Sequence 2
Rotate +90° Rotate -90°
about y axis about x axis
x z
y y y
z z
7.10 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
P
z
P1 y
k1
j1
i1
Z r0
R0 x
k0
j0
i0 Y
X
Figure 7 .11 Development of the Transformation Equations
r0 R0 P1
or in component form:
x0iˆ0 y0 ˆj0 z0 kˆ0 X 0iˆ0 Y0 ˆj0 Z 0 kˆ0 x1iˆ1 y1 ˆj1 z1kˆ1 (7.9)
We can determine the components of the above position vector in any direction by forming the dot
product of Equation 7.9 with a unit vector in the desired direction, i.e.:
x0 X 0 x1iˆ1 iˆ0 y1 ˆj1 iˆ0 z1kˆ1 iˆ0
where the dot product of the various unit vectors represents the "direction cosines" between the
coordinates. Direction cosines are discussed thoroughly in Chapter 2, Inertial Navigation Systems.
A point on a rigid body can be defined in terms of body-fixed (xyz) axes and by three independent
Euler angles , and defining the angular orientation of the body axes relative to the inertial axes
(XYZ). To do this, start with the body axes coinciding with the inertial axes at position 0. Then allow the
body axes to rotate about the z0-axis through an angle . The relationship between the two coordinates is
then given by:
x1iˆ1 y1 ˆj z1kˆ1 x0iˆ0 y0 ˆj0 z0 kˆ0
Employing the dot product concept introduced in Equations 7.9, we get the components with respect
to the new axes at position 1 as shown in Figure 7.12:
7.11 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
x0 x1
y0
z0 ,z1 y1
Figure 7 .12 Rotation about z0 through the Angle
x1 x0iˆ0 iˆ1 y0 ˆj0 iˆ1 z0 kˆ0 iˆ1
7.12 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
x1
z1 z2
7.13 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
y2
x2 ,x3
y3
z3 z2
Figure 7 .14 Rotation about x2 through the Angle
Expanding, gives:
x cos cos cos sin sin X
y cos sin sin sin cos cos cos sin sin sin sin cos Y
z sin sin cos sin cos sin cos cos sin sin cos cos Z
(7.10a)
The inverse of the transform matrix for converting from the body axes to the inertial axes coordinate
system can be found with a little effort (by matrix inversion) to be:
7.14 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
X coscos cossin sin sin cos sin sin cossin cos x
Y cossin coscos sin sin sin sin cos cossin sin y
(7.11)
x
y
p1 cos sin 0 0 0
q sin cos 0 0 0
1
r1 0 0 1 Horizontal reference
plane of the moving
earth axis system
Z,z
Figure 7.15 Development of Aircraft Angular Velocities by the Euler Angle Yaw Rate ( rotation)
(Figure 7.16):
The body axis system is then rotated about the y1-axis with an angular velocity of
7.15 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
Figure 7.16 Development of Aircraft Angular Velocities by the Euler Angle Yaw Rate ( rotation)
Finally, the body axis system is rotated about the x-axis with an angular velocity of
(Figure 7.17).
p3 1 0 0 sin sin
q 0
3 cos sin cos cos sin
r3 sin sin cos cos sin cos cos
y
q cos sin j x
Horizo
reference
ntal
plane
cos
P sin i (negative shown)
r cos sin k
c sin
z os
sin
Z
Figure 7.17 Development of Aircraft Angular Velocities by the Euler Angle Yaw Rate ( rotation)
Therefore, the body-axis components of p, q, and r, in terms of inertial-axis rates and Euler angles
are:
p sin (7.12)
7.16 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
cos
q cos sin (7.13)
r sin
cos cos (7.14)
Since , and
, were about the z0, y1, and x2 axes, they are not orthogonal. These equations
can be solved explicitly for , and
, since they are important forms of the equations. This can be
accomplished by first writing the above equations in matrix form:
p 1 0 sin
q 0
cos cos sin (7.15)
r 0 sin cos cos
We can then invert the transformation matrix in the right-hand side of equation 7.15 using the method
described in Vectors and Matrices (Chapter 4 & 5) and premultiplying both sides of the equation by the
inverse to get an equation solved explicitly for , and
, .
1 0 sin
First let A 0 cos
cos sin
0 sin cos cos
Replace each element of [A]T with its cofactor to get the adjoint matrix:
cos sin sin sin cos
adj A 0 cos cos cos sin
0 sin cos
Finally, dividing through by the determinant of [A], gives the inverse matrix:
1 tan sin tan cos
A 1 0 cos sin
0 sin / cos cos / cos
angles are:
p tan ( q sin r cos )
(7.16)
q cos r sin
(7.17)
( q sin r cos ) / cos (7.18)
7.17 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
The following list of angular rate and acceleration transformation is presented without proof. To
avoid confusion, the subscript b denotes body axis, s denotes stability axis, and w denotes the relative
wind axis system.
p b cos
sin
q b cos
sin
cos sin
sin sin
cos cos
rb cos cos sin cos cos sin sin cos
qs qb
q s q b
rs r b cos pb sin
pw ps cos qs sin
qw qs cos ps sin
rw rs
rw rs
7.18 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
The first two parameters above are easily measured using ground-based radar or onboard GPS or inertial
reference systems. If only and the Euler angles are available, then assuming zero winds, the flight
path angles can be calculated as
sin 1 sin cos cos cos sin cos cos sin sin
sin sin cos cos sin cos cos cos sin cos sin cos sin sin sin sin
sin 1
cos
cos sin cos sin sin
sin 1
cos
The rate for each of these parameters can be calculated as follows:
VT cos n x cos sin cos sin sin n y (cos sin ) n z (sin sin cos cos sin )
g
VgT n x cos sin sin sin cos n y cos sin ) n z (sin sin sin cos cos cos
p cos r sin sec sin VgT nx sin n z cos cos cos tan
Actually, these equations describe the velocity vector (angles relative to the air mass). If the air mass
is moving relative to the Earth, as is usually the case, the above equations do not describe the flight path.
m Vcg/XYZ
x
cg
y z
Y Z
Figure 7 .18 True Velocity of Body with Respect to Fixed-earth Axes Coordinate Systems
Velocity Transformation
Considering some rigid body with an element m, vector analysis shows that the velocity of m as seen
from an outside reference (X, Y, Z), is the sum of the body's cg velocity relative to that same reference
plus the body's angular velocity acting through the distance between m and the cg:
VT m / XYZ VT XYZ
r
VT , , and r are vectors defined as:
7.19 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
VT XYZ
uiˆ vˆj wkˆ (7.19)
iˆ ˆj kˆ
V TXYZ uiˆ vˆj wkˆ p q r
xyz
Expanding yields:
V Tm / XYZ uiˆ vˆj wkˆ qz ry iˆ pz rx ˆj py qx kˆ
Rearranging gives: V Tm / XYZ u qz ry iˆ v rx pz ˆj w py qx kˆ (7.22)
The XYZ in VTXYZ is any other reference system, typically the inertial coordinate system. In this case, it
can be considered as the total velocity of the element. Equation 7.22 shows that the total velocity of any
point on a moving body can be described by arranging component velocities according to the body's axes
(a.k.a. unit vectors). Not only can there be linear velocity along each body axis, but an additional
"coupled velocity” can exist. This comes from an angular rate acting at some distance from the cg.
An example of a coupled velocity is seen by an earth-based (inertial) observer who examines the right
wingtip of an aircraft yawing to the right. The positive yaw rate and positive lateral position combine to
give a tip motion towards the rear of the aircraft, a negative velocity. This is shown as the " ry"
term in Equation 7.22.
Acceleration Transformation
From vector analysis, the derivative of the velocity V T in the inertial (fixed-earth) coordinate
system is related to the derivative of V T along the body axis system through the relationship:
dVT dV
dt T VT (7.23)
XYZ dt xyz
Using Equations 7.19 and 7.20, the acceleration equation (7.23) can be written as:
iˆ ˆj kˆ
dVT ˆ
dt u iˆ v ˆj w
k p q r
XYZ u v w
7.20 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
dVT
ˆ
u qw rv iˆ (v ru pw) ˆj ( w pv qu )k
dt XYZ
This total transational acceleration can also be broken down into three components according to the
body's unit vectors
a x u qw rv (7.24)
a y v ru pw (7.25)
a z w pv qu (7.26)
These equations show how motions observed from a moving platform relate to accelerations in
inertial space. Consider Equation 7.22 in the case of an aircraft performing a level acceleration
maneuver: both pitch rate and yaw rate are zero, so the equation simplifies to a x u . This seems
straightforward enough, but if the aircraft is also yawing and translating sideways, then r and v combine
to create an additional acceleration not accounted for by the simple u term. This is another cross-
coupled effect. This effect can be visualized by imagining a top view of an aircraft with its cg moving
right (positive v) while yawing right (positive r). The ensuing acceleration of the cg due to these motions
is opposite to the overall aircraft acceleration u . Equation 7.24 also shows this to be a negative
acceleration along the aircraft's i component. A similar effect occurs with each of the cross-coupled terms
in the equations above.
Acceleration of the aircraft's flight path is used to describe turn capability and other kinematic
information. Flight path acceleration is identical to acceleration (or load factor) along each axis of the
relative wind axis system. The following transformation correlates linear acceleration from the
(measurable) body-axis system to the wind axis system. Note that they are just a specific application of
Equations 7.1 and 7.2.
N xw cos sin 0 cos 0 sin N xb
N sin cos 0 0 1 0 N yb
yw (7.27)
N xw 0 0 1 sin 0 cos N xb
The inverse of this matrix is easily shown to be
N xb cos 0 sin cos 0 sin 0 N xw
N 0 1 0 sin cos 0 N yw
yb
N xb sin 0 cos 0 0 1 N xw
The chapter on inertial navigation systems combines the transformations from Equations 7.10 and
7.27 to convert initially sensed motions & angles to flight path motions.
7.21 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
y m1
r1
r1
H r V
The specific angular momentum of m1 is: m 1 1 (7.30)
First consider the above velocity term. From vector analysis, the rate of change of the radius vector r1
(V1) can be related to the body axis system by:
V1
d r1
dt XYZ
d r1
dt xyz
r1
Assuming the aircraft is a rigid body, then r1 does not change with time and the first term can be
excluded. The above equation then simplifies to:
V1 r1
Applying matrix algebra and Equations 7.20 and 7.21 again gives
iˆ ˆj kˆ
r1 p q r
x y z
7.22 8/10/19
Volume I, Math and Physics for Flight Testers Chapter 7 Axis Transformations
iˆ ˆj kˆ
H
r1 V1 r1 r1 z y z
m
qz ry rx pz py qx
Finally, this matrix can be expanded to show the specific angular momentum of any single element for
each of the three body axis components.
H ˆ
m i y ( py qx) z (rx pz ) p( y z ) q ( xy ) r ( xz)
2 2
(7.31)
H ˆ
m j z ( pz qy ) x( py qx) q( x z ) r ( yz ) p ( xy)
2 2
(7.32)
H ˆ
m k x rx pz y qz ry r x y p xz q yz
2 2
(7.33)
The units for each term above is ft2/sec or [ ft/sec] ft, i.e.; velocity about a moment arm as shown in
Figure 7.19. The total angular momentum for a real object can be calculated by adding up the angular
momentum for each mass element. Another way to determine this is to integrate Equations 7.31 - 7.33
across the density of the object. This will be presented when developing the aircraft Equations of Motion,
Chapter 4.
Summary
This chapter defines the most common axis systems used in aeronautics and establishes the
practices used for transforming motion from one system to another. These procedures can be extended to
other axis systems such as "North, East, Down" as done in Chapter 2, Inertial Navigation Systems. Many
transformations exist which can be used to "simulate" parameters that cannot be practically instrumented
directly. The scope of this text is limited to the most common and useful transformations required for
developing the equations of motion and other basic flight information.
References
Anon., Aircraft Flying Qualities, Chapter 4, Equations of Motion, USAF Test Pilot School notes, AFFTC
Edwards AFB CA, March 1991.
Kalviste, Juri, Flight Dynamics Reference Handbook, Northrop Corp. Aircraft Division, April 1988.
7.23 8/10/19