The flat-Earth aircraft equations of motion are used for short-distance flight only and for observation of the aircraft responses to control inputs. Here the aircraft equations of motion for spherical, rotating Earth are derived. These equations are widely used for aircraft simulation.
Derivation
Nomenclature
then the direction cosine matrix from the ECEF to the body-fixed coordinate system is:
Simplifying by gathering terms gives:
Some of the derivatives can be rewritten as:
If the above definitions are substituted in the acceleration equation the following equation is obtained:
Assuming the Earth's rotational velocity is constant:
This equation gives the derivative of the velocity vector in the body-fixed coordinate system with respect to the EFEC coordinate system. However, since the body-fixed coordinate system rotates with respect to the EFEC coordinate system the following equation is used:
Substitution of the above equation in the previous gives:
The above equation gives the derivative of the velocity vector in the body-fixed coordinate system with respect to the body-fixed coordinate system. The flat-Earth aircraft equations of motion can be identified. However, there are additional terms for the Coriolis and Centripetal acceleration due to the Earth's rotational velocity. Since angular velocities are additive, the angular velocity of the body-fixed coordinate system relative to the ECI the coordinate system in body coordinates is given by:
Further simplifying the above equation gives:
The above are the translational aircraft equations of motion. The rotational aircraft equations of motion are given by:
The angular velocity of the NED coordinate system with respect to the ECEF coordinate system is given by:
where
Then the angular velocity of the body-fixed relative to the NED coordinate system is given by:
Non-linear state-space form
[C]b/a- Direction cosine matrix transformation from coordinate system a to b()a- Vector in coordinate system ab()˙a- Derivative of vector in coordinate system awith respect to coordinate system b()b/a- b with respect to aμ- latitude
l- longitude
h- altitude
Re- Earth radius
In order to derive the aircraft equations for spherical, rotating Earth four coordinate systems are defined. The Earth-centred inertial (ECI) coordinate system has its origin at the centre of the Earth with its xb-axis pointing towards the vernal equinox and its zb-axis pointing towards the Earth's north pole. The Earth-fixed centred coordinate system (EFEC) has its origin at the centre of the Earth with its xb-axis pointing towards the prime meridian and zb-axis towards the Earth's north pole. The north-east-down (NED) coordinate system has its origin at a particular longitude, latitude and altitude with respect to the EFEC coordinate system, its xb-axis is pointing north and its zb-axis is pointing toward the centre of the Earth. The last coordinate system is the body-fixed coordinate system. Its origin coincides with the NED coordinate system origin.
The direction cosine matrices between the EFEC, NED and body-fixed coordinate systems must be defined. The direction cosine matrix from the EFEC to NED coordinate system in terms of longitude (μ) and latitude (l) is:
A position vector r in the ECEF coordinate system gives the position of the NED and the body-fixed coordinate systems (the NED and body-fixed coordinate systems coincide, as mentioned earlier). The time derivative of rwith respect to the ECI coordinate system is given by:
ir˙e=er˙e+ωe/ie×re,
where ωe/iis the earth rotational velocity. Differentiating again with respect to time gives:
and are similar to the flat-Earth rotational equations of motion with the difference that ωb/ib is the aircraft angular velocity vector with respect to the ECI coordinate system. The components of this angular velocity vector are not directly equal to p, q, and r which is the case for the flat-Earth equations of motion. Here they are given by ωb/fb. In order to determine p, q and r, the rotation of the NED with respect to the EFEC coordinate system must be taken into account. Why is rotation present between the NED and EFEC coordinate systems? The NED coordinate system is always tangent to the Earth. If the aircraft starts its journey at l=0 deg and ends its journey at l=90 deg, the NED coordinate system would have rotated in order to stay tangent to the Earth. So the angular velocity of the body-fixed to the ECEF coordinate system can be represented as:
ωb/eb=ωb/fb+[C]b/fωf/ef,
where ωb/fbis the angular velocity of the body-fixed with respect to the NED coordinate system and ωf/efis the angular velocity of the NED with respect to the ECEF coordinate system.