Three-Dimensional Rigid-Body Dynamics

Introduction and context

Three-dimensional rigid-body dynamics describes the motion of a body whose dimensions and mass distribution remain unchanged while it translates and rotates in space. Unlike planar motion, spatial motion generally involves three coupled components of angular velocity, angular acceleration, torque, and angular momentum. The central difficulty is that angular momentum is not usually parallel to angular velocity. Their relationship depends on the body's mass distribution and is represented by the inertia tensor.

These ideas are fundamental in rotorcraft, spacecraft, vehicle stability, turbine rotors, robotic manipulators, reaction wheels, gyroscopes, and rotating machinery. A successful analysis requires a clear distinction between inertial and body-fixed reference frames, correct use of vector transport, and careful interpretation of principal axes.

Learning objectives

By the end of this chapter, a student should be able to:

1. Define spatial angular velocity and angular acceleration.
2. Construct and interpret a three-dimensional mass moment-of-inertia tensor.
3. Determine principal axes and principal moments of inertia.
4. Apply Euler's equations to a rigid body in a body-fixed frame.
5. Explain gyroscopic torque, free precession, and nutation.
6. Analyze a rotating engineering system and identify modeling limitations.
7. Use computational tools responsibly while preserving physical and dimensional checks.

1. Kinematics of spatial rotation

A rigid body's orientation may be described by three mutually perpendicular unit vectors fixed in the body, or by a rotation matrix [[M:R]] that maps body components into inertial-frame components. A proper rotation matrix satisfies

where [[M:I]] is the identity matrix. The angular velocity vector [[M:ω]] is independent of the particular orientation coordinates used. It is defined so that the velocity of a point fixed in the body is

Here [[M:v_P]] and [[M:v_O]] are inertial velocities, and [[M:r_{P/O}]] is the position of point P relative to reference point O. The equation assumes a rigid body and an inertial observation frame.

Angular acceleration is the inertial derivative of angular velocity:

If angular velocity is resolved in body-fixed axes, the transport theorem gives

for any vector [[M:a]]. Consequently,

when [[M:a = ω]], because [[M:ω × ω = 0]]. Thus the numerical components of angular acceleration in body axes are the time derivatives of the components of angular velocity in those same axes. This special simplification must not be incorrectly applied to an arbitrary vector.

For a point P fixed in the rigid body, its acceleration is

The second term is tangential acceleration and the third is centripetal acceleration. The expression assumes that O is a point fixed in the body and that the body is rigid.

2. Mass moments of inertia and the inertia tensor

For a body of density [[M:ρ]] occupying volume V, the inertia tensor about point O is

where [[M:r]] is the position vector from O, [[M:r² = r·r]], and [[M:I]] is the identity tensor. In Cartesian body axes,

with

and products of inertia

The signs in the matrix arise from the tensor definition; the products themselves are defined here without a minus sign. The tensor has units of kg·m² and is symmetric. Symmetry guarantees real orthogonal eigenvectors and real principal moments.

The angular momentum about O is

when O is fixed in the body and the components of both quantities are expressed in the same axes. In general, [[M:H]] is not parallel to [[M:ω]]. Parallelism occurs when the rotation is about a principal axis or when the inertia tensor is isotropic for the relevant motion.

Parallel-axis theorem

If the center of mass G is known and O is displaced from G by vector [[M:d]], then

where m is total mass. For a planar-looking special case in which the displacement is d along the x-axis, [[M:I_{O,yy} = I_{G,yy} + md²]] and [[M:I_{O,zz} = I_{G,zz} + md²]].

Principal axes

Principal axes are directions for which the products of inertia vanish. They are obtained from

where [[M:e_i]] is a unit eigenvector and [[M:I_i]] is the corresponding principal moment. In principal axes,

For a body with reflection symmetry, some products of inertia are zero immediately. For example, a uniform body symmetric about the x-y plane has [[M:I_{xz = I_{yz = 0}]] when the origin is at the symmetry center. Symmetry is a physical shortcut, whereas diagonalization is the general method.

3. Angular momentum and kinetic energy

The rotational kinetic energy about the center of mass is

This expression is always nonnegative for a physical mass distribution. If the body also translates, total kinetic energy is

The angular momentum about an arbitrary point O is related to that about G by

provided the velocities and position vectors are measured in the same inertial frame. For dynamics about a fixed point, the angular-momentum equation is

For a freely translating rigid body, use the center of mass:

External forces determine translation through [[M:∑F = ma_G]], while external moments determine rotation.

4. Euler's equations

Let body-fixed axes coincide with principal axes at the center of mass. Write [[M:ω = (ω_1,ω_2,ω_3)]] and [[M:M = (M_1,M_2,M_3)]]. Applying the transport theorem to [[M:H = Iω]] produces Euler's equations:

Here [[M:dot(ω_i)]] means the derivative of the body component with respect to time, [[M:I_i]] are constant principal moments, and [[M:M_i]] are applied external-moment components. These equations assume a rigid body, body-fixed principal axes, and an inertial frame for the applied moment balance.

The cross-coupling terms are gyroscopic terms. They do not represent additional forces; they result from the changing orientation of the body-fixed axes. If the body rotates exactly about a principal axis and no disturbance exists, the cross terms may vanish. For a spherical body, [[M:I_1 = I_2 = I_3]], all gyroscopic coupling terms vanish.

For a body-fixed frame that is not principal, the compact vector equation is

where I is constant in the body frame. This is often the safest computational form.

5. Gyroscopic effects and precession

A rapidly spinning rotor possesses angular momentum approximately along its spin axis. If an external moment is applied perpendicular to that angular momentum, the angular-momentum vector changes direction rather than merely changing magnitude. The governing relation is

for steady precession, where [[M:Ω_p]] is the precession angular velocity. If [[M:Ω_p]] is perpendicular to [[M:H]], then

in magnitude. For a rotor with spin inertia [[M:I_s]] and spin speed [[M:ω_s]], [[M:H ≈ I_sω_s]], so

This approximation requires a dominant spin angular momentum, small wobble, and a moment perpendicular to the spin axis. It is not a universal formula for arbitrary three-dimensional rotation.

A familiar example is a top under gravity. For a symmetric top with mass m, center-of-mass distance l from the support, spin angular momentum H, and slow steady precession, the gravitational moment magnitude is approximately mgl. Hence

The precession direction follows the vector equation, not a memorized scalar sign.

Free motion may include nutation, an oscillation of the body axis about a mean direction. For an axisymmetric body, rapid spin can stabilize orientation, but disturbances may produce coning or wobble. In engineering systems, these motions can excite bearings, control actuators, structural modes, and sensors.

6. Stability of torque-free rotation

For torque-free motion, [[M:M = 0]]. If the body rotates exactly about a principal axis, it is an equilibrium rotation. Linearization shows that rotation about the axis associated with the smallest or largest principal moment is stable in the ideal torque-free model, while rotation about the intermediate principal moment is unstable. The latter is known as the intermediate-axis or tennis-racket instability.

This result does not mean every real machine will become unstable about the intermediate axis immediately. Damping, flexible deformation, control systems, bearing forces, and energy dissipation modify the behavior. It does mean that small disturbances can grow in the ideal rigid-body equations.

7. Worked example: angular momentum and gyroscopic torque

A rotor has principal moments [[M:I_1 = 2.0]], [[M:I_2 = 3.0]], and [[M:I_3 = 5.0]] kg·m². Its body-frame angular velocity is

[[M:ω = (4, 2, 10)]] rad/s, and its applied body-frame moment is zero.

The angular momentum is

Using Euler's equations with zero applied moment:

Thus, even without an applied moment, the body components of angular velocity change because the body-fixed axes rotate. Angular momentum remains constant in inertial space, but its components in the body frame do not. The calculation illustrates why setting [[M:M = 0]] does not imply constant body components of angular velocity.

8. Applied case: spacecraft reaction-wheel attitude control

Background

An illustrative spacecraft carries a reaction wheel whose spin axis is aligned approximately with the spacecraft roll axis. The spacecraft is in orbit, where external aerodynamic moments are small, but solar radiation pressure, magnetic interaction, and gravity-gradient effects produce disturbances.

Problem

A control system must rotate the spacecraft without expelling propellant. The wheel changes its angular momentum, and the spacecraft responds with an equal and opposite angular-momentum change.

Analysis and approach

For a wheel-spacecraft system isolated from external moment over a short interval,

If the wheel angular momentum changes by [[M:∆H_w]], the spacecraft change is approximately

The actuator torque is generated by changing wheel speed. In practice, the controller measures attitude and angular rate with gyroscopes, estimates the state, and commands motor torque. A simplified single-axis model is

where the signs depend on the selected positive directions, [[M:I_s]] and [[M:I_w]] are spacecraft and wheel inertias, and [[M:M_dist]] represents disturbance moment. A complete model uses the full inertia tensor and three coupled Euler equations.

Outcome and engineering interpretation

Reaction wheels permit precise pointing but accumulate momentum when persistent disturbance moments act. Once wheel speed approaches a limit, the spacecraft must unload momentum using magnetic torquers or thrusters. A design must therefore consider wheel saturation, bearing life, structural vibration, sensor bias, flexible appendages, and fault tolerance.

Learning points

1. Internal motor torque redistributes angular momentum but cannot change total angular momentum without an external moment.
2. A diagonal inertia model is valid only when axes are principal or when off-diagonal terms are demonstrably negligible.
3. Control performance depends on both rigid-body dynamics and actuator, sensor, and structural limitations.

9. Applications in engineering practice

Rotors and turbines. Imbalance produces periodic forces and moments. Misalignment can create products of inertia and coupling between bending and rotation. Rotor models often combine rigid-body gyroscopic terms with flexible-shaft finite-element models.

Vehicles. Vehicle body roll, pitch, and yaw are coupled when mass distribution is asymmetric or when the velocity and body axes differ. Tire forces, suspension compliance, and aerodynamic moments must be added to the rigid-body equations.

Aerospace systems. Aircraft, spacecraft, and launch vehicles use attitude dynamics for guidance and control. Inertial sensors measure angular rate, while estimators combine sensor data with a dynamic model.

Robotics. A manipulator's links have configuration-dependent inertia tensors. The same principles extend to joint-space equations, but the inertia matrix changes with configuration and interactions among links are significant.

10. Modeling, computation, and limitations

A reliable workflow is:

1. Select inertial and body-fixed frames.
2. Locate the center of mass and specify the reference point.
3. Construct the inertia tensor using geometry, CAD data, measurement, or integration.
4. Check symmetry, units, positive definiteness, and plausible magnitudes.
5. Identify principal axes when useful.
6. Apply force and moment balances.
7. Integrate the equations while monitoring energy, angular momentum, and orientation errors.
8. Compare predictions with experimental data.

Euler angles are intuitive but can suffer from singularities, including gimbal lock. Rotation matrices avoid coordinate singularities but require numerical re-orthogonalization. Quaternions provide an efficient nonsingular orientation representation, subject to the unit constraint

Numerical software and AI-assisted tools can calculate eigenvalues, integrate Euler equations, and fit inertia parameters from test data. They do not replace frame definitions, dimensional analysis, validation, or engineering judgment. A physically impossible negative principal moment, loss of matrix orthogonality, or unexplained energy growth should be treated as a modeling or numerical error.

11. Design and safety considerations

Rotating equipment stores substantial kinetic energy, approximately [[M:T_rot = ½Iω²]] for single-axis rotation. Failure containment, overspeed protection, balancing, guarding, and bearing-load assessment are essential. In aerospace and automotive applications, incorrect sign conventions or neglected cross-coupling can produce unstable control commands. Traceable assumptions, independent checks, and conservative verification are therefore part of professional practice.

The rigid-body model is limited when deformation, fluid-structure interaction, internal damping, bearing flexibility, mass variation, or contact impacts are important. In such cases, rigid-body dynamics should be coupled to flexible-body, multibody, computational-fluid-dynamics, or experimentally identified models.

Chapter Summary

Three-dimensional rigid-body dynamics requires simultaneous treatment of translation, orientation, angular velocity, angular acceleration, angular momentum, and applied moments. The inertia tensor captures mass distribution and connects angular velocity to angular momentum. Its eigenvectors define principal axes, where the tensor becomes diagonal and Euler's equations take their clearest form. The body-frame equation [[M:M = I dot(ω)_B + ω × (Iω)]] explains gyroscopic coupling and why zero applied moment does not imply constant body components of angular velocity. Gyroscopic torque produces precession, while torque-free motion can exhibit nutation and intermediate-axis instability. These principles support the analysis and design of rotors, vehicles, spacecraft, and robots, but reliable engineering practice also requires correct frame selection, numerical validation, safety analysis, and recognition of flexible-body limitations.

Key Terms

Angular velocity: Vector describing the instantaneous rate and axis of rotation.

Angular acceleration: Inertial time derivative of angular velocity.

Inertia tensor: Second-order tensor relating angular velocity to angular momentum.

Mass moment of inertia: Measure of resistance to angular acceleration about an axis.

Product of inertia: Off-diagonal mass-distribution integral coupling two axes.

Principal axis: Eigenvector direction of the inertia tensor.

Principal moment: Eigenvalue of the inertia tensor associated with a principal axis.

Angular momentum: Rotational momentum given by [[M:H = Iω]] about a specified point.

Transport theorem: Relation between inertial and rotating-frame derivatives.

Euler's equations: Body-frame rotational equations of motion.

Gyroscopic torque: Moment associated with changing the direction of angular momentum.

Precession: Rotation of an angular-momentum or body axis about another direction.

Nutation: Oscillatory variation of a spinning body's orientation.

Reaction wheel: Internal wheel used to control spacecraft attitude.

Intermediate-axis instability: Instability of ideal torque-free rotation about the middle principal moment.

Rotation matrix: Orthogonal matrix mapping components between orientations or frames.

Body-fixed frame: Coordinate frame attached to the rigid body.

Inertial frame: Frame in which Newton's laws take their standard form.

Rigid-body assumption: Idealization that distances between all body points remain constant.

Review Questions

5 Short Answer Questions

1. Why are angular momentum and angular velocity not generally parallel in three-dimensional motion?
2. What physical information is represented by an inertia tensor?
3. State the transport theorem for a vector observed from a rotating frame.
4. What is meant by a principal axis?
5. Why can a torque-free rigid body have changing angular-velocity components in body-fixed axes?

5 Long Answer Questions

1. Derive the body-frame vector equation [[M:M = I dot(ω)_B + ω × (Iω)]] from the angular-momentum balance and transport theorem.
2. Explain the construction, symmetry, and physical interpretation of the inertia tensor, including the parallel-axis theorem.
3. Discuss gyroscopic precession and derive the approximate relation [[M:Ω_p ≈ M/(I_sω_s)]] under clearly stated assumptions.
4. Compare rotation matrices, Euler angles, and quaternions for representing rigid-body orientation.
5. Explain the intermediate-axis instability and discuss why real engineering systems may behave differently from the ideal rigid-body model.

3 Analytical / Application Questions

1. A rotor has principal moments 1.5, 2.0, and 4.0 kg·m² and body angular velocity components 3, −2, and 8 rad/s. With zero applied moment, calculate the three initial body-frame angular accelerations using Euler's equations.
2. Design a conceptual reaction-wheel attitude-control strategy for a small spacecraft subject to a constant disturbance moment. Identify the angular-momentum flow, wheel-saturation problem, sensors, and momentum-unloading method.
3. A vehicle model shows large off-diagonal inertia terms after its center of mass is moved. Propose a diagnostic procedure to determine whether the terms arise from an incorrect frame, an incorrect center-of-mass estimate, or a genuine asymmetric mass distribution.

Multiple Choice Questions with Answers

Question 1

Which relation gives angular momentum about the center of mass for a rigid body expressed in the same axes as angular velocity?

A. [[M:H_G = I_Gω]]

B. [[M:H_G = m v_G]]

C. [[M:H_G = I_Gα]]

D. [[M:H_G = ω × I_G]]

Correct Answer: A. [[M:H_G = I_Gω]]

Question 2

What is the defining property of principal axes?

A. They are always parallel to the inertial axes.

B. Products of inertia vanish in those axes.

C. They always pass through an arbitrary point O.

D. They make angular momentum zero.

Correct Answer: B. Products of inertia vanish in those axes.

Question 3

A rotor has inertia 4 kg·m² about its spin axis and spins at 50 rad/s. Its approximate angular momentum magnitude is:

A. 0.08 kg·m²/s

B. 12.5 kg·m²/s

C. 200 kg·m²/s

D. 400 kg·m²/s

Correct Answer: C. 200 kg·m²/s

Question 4

For a body-fixed principal-axis frame, which term represents gyroscopic coupling in the first Euler equation?

A. [[M:I_1 dot(ω_1)]]

B. [[M:(I_3−I_2)ω_2ω_3]]

C. [[M:I_1ω_1²]]

D. [[M:m a_G]]

Correct Answer: B. [[M:(I_3−I_2)ω_2ω_3]]

Question 5

In steady precession with perpendicular applied moment and angular momentum, the magnitude relation is:

A. [[M:M = H/Ω_p]]

B. [[M:M = Ω_pH]]

C. [[M:M = Ω_p/H]]

D. [[M:M = H²Ω_p]]

Correct Answer: B. [[M:M = Ω_pH]]

Question 6

Which ideal torque-free rotation is generally unstable to small disturbances?

A. Rotation about the smallest principal moment.

B. Rotation about the largest principal moment.

C. Rotation about the intermediate principal moment.

D. Rotation about every principal axis equally.

Correct Answer: C. Rotation about the intermediate principal moment.

Question 7

The units of a mass moment of inertia are:

A. kg·m

B. kg/s²

C. kg·m²

D. N·m/s

Correct Answer: C. kg·m²

Question 8

For a freely translating rigid body, the safest standard point for rotational angular-momentum balance is usually:

A. The center of mass.

B. Any point on the body without correction.

C. The farthest geometric corner.

D. The instantaneous point of zero velocity only.

Correct Answer: A. The center of mass.

Question 9

What is a principal limitation of Euler angles?

A. They cannot represent any finite rotation.

B. They violate conservation of angular momentum.

C. They can encounter coordinate singularities such as gimbal lock.

D. They cannot be used in numerical integration.

Correct Answer: C. They can encounter coordinate singularities such as gimbal lock.

Question 10

A reaction wheel changes its angular momentum by ∆H_w in an isolated spacecraft system. Approximately what spacecraft angular-momentum change results?

A. ∆H_s ≈ ∆H_w

B. ∆H_s ≈ 0 for all conditions

C. ∆H_s ≈ −∆H_w

D. ∆H_s ≈ (∆H_w)²

Correct Answer: C. ∆H_s ≈ −∆H_w