1. Introduction
Modeling flight dynamics with tensors enables building complex simulations by taking the two-step approach from tensor modeling to matrix programming. The new paradigm is called tensor flight dynamics. It generalizes the hitherto traditional vector flight dynamics as taught by Kane [1] , Likins [2] , Hughes [3] , and others. Vector analysis is based on Gibbs [4] , who simplified Hamilton’s [5] quaternions to three dimensional space, while tensor flight dynamics uses the Ricci and Levi-Civita [6] tensor formalism.
A tensor of rank two is in traditional dynamics a dyad, which is expressed in a particular set of orthogonal unit vectors. Special rules apply for the multiplication of a dyad with a vector [2] . Hughes [3] goes one step further and uses the concept of vectrix to express a dyad in basis vectors.
While vector flight dynamics uses only one dyad, i.e., the moment of inertia dyad, tensor flight dynamics uses three tensors of rank two, namely the moment of inertia, rotation, and angular velocity tensor. Simple tensor multiplication is used to combine them with tensors of rank one.
The moment of inertia tensor
of body B (superscript) referred to the center of mass B (subscript) is a real, symmetrical tensor. (Tensors are expressed in bold face, upper case for rank two and lower case for rank one). When expressed in
body coordinates
, the nine elements of the 3 × 3 matrix
provide the numerical values of the axial and products of inertia.
The rotation tensor
describes the orientation of frame B with respect to (wrt) frame A and is an orthogonal tensor. Let B be an airplane and A the airport. Then
models the attitude of the airplane wrt the airport without reference to any coordinate system. If we want to calculate the attitude angles of the aircraft relative to the airport, we introduce two coordinates systems
and
associated with the airplane and the airport, respectively. Their transformation
matrix
, containing the attitude angles, is obtained from the rotation tensor expressed in either coordinate system
, where the overbar indicates the transpose of the matrix.
The angular velocity tensor
models the angular rate of frame B wrt frame A. Because it is a skew-symmetric tensor of rank two it can be contracted to a tensor of rank one,
. (For brevity, tensors of rank two are called tensors, while tensors of rank one are called vectors). In vector analysis the angular velocity vector is called an axial vector, which has to abide by the right-hand rule, because, in actuality, it is contracted from a tensor of rank two.
Tensor flight dynamics elevates the modeling of flight dynamics to a coordinate-invariant form. It is based on Einstein’s Covariance Principle [7] , Page 107, which states, “All physical laws are invariant under all coordinate transformations”. While Einstein based his special and general theory of relativity on the covariance principle, tensor flight dynamics, a branch of classical mechanics, applies this principle to Newton’s laws of motion.
However, to achieve this invariance, the ordinary time derivative, must be replaced by a new time operator called the rotational time derivative ( [8] , Section 4.2.1). The ordinary time derivative applied to a vector destroys its tensor characteristics, if the vector is subjected to a time-dependent coordinate transformation. By replacing the ordinary time operator in Newton’s Second Law by the rotational time derivative, the tensor characteristic of the law is maintained. With p the linear momentum, f the externally applied forces, and
the rotational time derivative wrt the inertial frame I, Newton’s Second Law is expressed in tensor form
which is valid in all Cartesian coordinate systems even those that are related by time dependent coordinate transformations.
Besides the rotational time derivative, the Euler transformation is instrumental to enable tensor flight dynamics. It governs the shift of reference frames. Given frames I and B and the angular velocity tensor of frame I wrt B,
, the rotational time derivative of tensors of rank one x transforms from frame I to frame B like
(1)
This is a tensor relationship valid in all coordinate systems. Its equivalent in vector mechanics is for vector x
The proof of the Euler transformation for tensors of rank one can be found in [8] , Appendix D. Here we present and prove the Euler transformation for tensors of rank two. It is the foundation for deriving the relationship between the rotation tensor and the angular velocity tensor, which hitherto was based on a definition, both in vector mechanics [2] and in tensor flight dynamics [8] , rather than on a strict mathematical derivation.
The Euler transformation of tensors of rank two has some precursors in vector mechanics. In [2] , Equation 8.93, we find for the dyad D and the coordinate systems a and b
However, because of the ordinary time derivative, this relationship is only valid for the two coordinate systems. Similarly, [3] , Appendix B.4, Equation 17 provides the relationship for the dyad D
where
is the vectrix of coordinate system b. Again, the ordinary time derivative (over-dot) limits this expression to the two coordinate systems a and b.
As we will see, the Euler transformation of tensors of rank two has some similarities with these precursors from vector mechanics, but because it is valid for all coordinate systems it is a true tensor concept.
2. Euler Transformation for Tensors of Rank Two
Though the Euler transformation for tensors of rank one is sufficient for almost all flight dynamics situations, this extension to tensors of rank two makes possible a proof of an important kinematic relationship, which is given in Section 3. Additional examples are then provided in Sections 4 and 5.
We start with the formal statement of the new transformation, followed by the proof.
Let A and B be two arbitrary frames related by the angular velocity tensor
. Then, for any tensor X of rank two, the following transformation of the rotational time derivative holds:
(2)
The overbar indicates the transpose.
For the proof, we start with the Euler transformation for tensors of rank one, Equation (1)
and substitute
, where y is an arbitrary tensor of rank one
Apply the chain rule
Now the rotational time derivative of the second term on the left-hand side is transformed from frame A to frame B
and introduced into the previous equation
The second terms on the left- and right-hand sides cancel; and since y is arbitrary, we have
Bringing the second term on the left-hand side to the right-hand side and using the fact that the angular velocity tensor is skew-symmetric,
, we have proven the transformation law of tensors of rank two
3. Angular Velocity Tensor
The angular velocity tensor of rank two
of frame B wrt to frame A is a skew-symmetric tensor and can therefore be contracted to a tensor of rank one
. Whether first or second rank, their derivations in the past were based on geometric consideration, without mathematical substantiation. Now, with the Euler transformation for tensors of rank two, we can provide such a mathematical derivation using the rotation tensor.
The rotation tensor
of rank two establishes the orientation of any frame B wrt frame A. The angular velocity tensor
is then postulated ( [8] , Equation 4.47) to be obtained by
(3)
To prove this relationship, we apply Equation (2) to the rotation tensor
and post-multiply the right and left-hand sides by
The second term on the right-hand side is
, because
is orthogonal and thus
(unit tensor). Together with the left-hand term we get our relationship, Equation (3).
What remains to be shown is that
is satisfied identically. We pre-multiply with
and post-multiply with
. With
we get
(4)
which is almost again in the form of Equation (3). To show that, we resort to a trick using the fact that the rotational time derivative of the unit tensor E is zero
We apply the chain rule
and modify the last term, using
( [8] , Section 4.1.1)
Substituting into Equation (4) we get
But because
( [8] Section 4.2.4.1), we obtain again our Equation (3) with the frames A and B reversed
Because this part of the proof is satisfied identically, we have succeeded in deriving the mathematical relationship of the angular velocity tensor
with the rotation tensor
as shown in Equation (3)
4. Derivation of the Strap-Down INS Equations
Any kind of tactical missile today has a strap-down INS. Its sensors consist of gyros and accelerometers mounted on the missile body frame B. Integrating the accelerometer measurements twice yields the missile position, provided they are first converted to the inertial frame I with the help of the gyro measurements. We will derive these equations that need be programmed for the navigation computer. The gyros measure the angular velocity
of the missile body B wrt to the inertial frame I. This measurement is related to the attitude of the missile by the rotation tensor
, as expressed by Equation (3)
Post-multiplying by
yields the time differential equations that govern the attitude of the missile, given the gyro measurements
(5)
To implement this relationship in the navigation processor we must convert the tensors to matrices by introducing coordinate systems. Though
is measured in body coordinates, we start by applying inertial coordinates to Equation (5)
The rotational derivative wrt to frame I becomes the ordinary time derivative, as the coordinate system
is associated (fixed) with frame I
Because the gyro measurements are made in missile body coordinates,
, we need to introduce the coordinate transformation matrix
of body coordinates
wrt inertial coordinates
,
Now, because the coordinate systems
and
are associated with frames B and I, respectively, there exists a relationship between the rotation tensor
and the transformation matrix
, namely
( [8] , Section 4.1.1), which we now introduce
with
and
( [8] , Section 3.2.1.3)
(6)
here we have the nine linear differential equations to be programmed for the navigation processor. The transfer alignment will provide the initial conditions.
Applying
to the measured accelerations
yields the acceleration in inertial coordinates, ready to be integrated for the missile’s velocity and position after a gravity compensation.
5. Effect of Variable Moment of Inertia on Missile Dynamics
It is common practice in flight dynamics to neglect the time-rate-of-change of mass except in rocketry, where the expulsion of fuel is called thrust, but considered an external force. The same simplification is made in the attitude equations of motion, where the time-rate-of-change of the moment of inertia tensor is neglected.
Thomson [9] discusses this effect as early as 1966. He derives the equations of motion with vector notation, and uses a fixed body point as reference. Quadrelli [10] addresses the variable mass effects of a supersonic inflatable decelerator using a coordinate-free approach based on Reynolds’ transport theorem from continuum mechanics. Here we use our tensor formulation to investigate the effect of the varying moment of inertia on a short-range air-to-air missile and refer the equations of motion to the center of mass.
The attitude equations of motion are governed by Euler’s law: The inertial time-rate-of-change
of the angular momentum
of the body B wrt the inertial frame I, referred to the c.m. B equals the eternal moments
applied at the c.m. B
where the angular momentum is the product of the moment of inertia tensor
of body B referred to the c.m. B and the angular velocity
of body B wrt the inertial frame I
Substitution into the previous equation yields
here the moment of inertia
is a tensor of rank two,
and
are tensors of rank one. We apply the Euler transformation for tensors of rank one, Equation (1), in order to shift the rotational time derivative to frame B
and apply the chain rule to the argument of the rotational time derivative
(7)
In flight dynamics it is usually assumed that
is negligibly small, so
that the attitude equations of motion assume the simpler form
This assumption may be acceptable for aircraft. But is it also justified for missiles, particularly for air-to-air missiles, whose fuel may be 40% of the launch mass and expended in less than 3 seconds? To pursue this investigation, we use a typical short range air-to-air missile, called AIM6 [11] , and use the complete attitude equations of motion, Equation (7). Introducing body coordinates
and multiplying out the matrices to obtain the scalar equations in the state variables
, yields the attitude equations of motion
where the over-dot indicates the ordinary time derivative, and where the moment of inertia for a missile with tetragonal symmetry (
) is given in body coordinates
with the external moments
.
It appears that the changing moment of inertia terms introduce a lag into the three channels. Let’s assume the missile executes a pure yaw maneuver, then
and what remains is
This is a simple first order lag system with the time constant
. However, because
is negative, the time constant is negative and we are faced with destabilization. How serious is this effect?
I ran a horizontal engagement of the prototype missile AIM6 using its 6 DoF simulation [11] and recorded
at 1.345 sec, Mach=2, and lateral acceleration 33 g’s. The time constant at that point is
sec. Since the thrust phase lasts only 2.5 sec, the destabilizing effect should not be too severe.
To see the actual evidence, I used the 6 DoF simulation without autopilot, which is possible because the missile is aerodynamically stable. Figure 1 displays the results of two runs with and without the unsteady moment of inertia, taken around the time of our test point at 1.345 sec.
The unsteady effect lowers the pitch angle and increases the amplitude of the pitch rate. The discrepancies are small and can be neglected. With the autopilot engaged they are not noticeable.
6. Conclusion
With the proof of the Euler transformation for tensors of rank two, all elements of tensor flight dynamics are now in place. This transformation provides us with a mathematically sound derivation of the relationship between angular velocity
Figure 1. Effect of the time-rate-of-change of moment of inertia.
and rotation tensor, which we used for the straight-forward derivation of the navigation equations of a strap-down INS. By applying the transformation to the moment of inertia tensor, we could discuss how its rapidly changing values affect the dynamics of an air-to-air missile. Because the effect is negligible, the common practice of neglecting the time-rate-of-change of moment of inertia is justified.