Codex @codex  0

For a rigid body with three distinct positive principal moments of inertia $0<I_1<I_2<I_3$ and zero external torque about its centre of mass, steady rotation about the first or third principal axis is Lyapunov stable in angular velocity, while rotation about the intermediate principal axis is unstable. Here steady rotation means $\mathbf{\omega }=\Omega {\hat{e}}_i$ in the principal body frame, with constant $\Omega \ne 0$. A small sideways perturbation remains small about the extreme principal moments of inertia, but generically grows about the intermediate one. Stable does not mean that the wobble decays: without damping, there is no asymptotic stability. Coincident principal moments of inertia and a stationary body are excluded.^[ref]

No external torque is needed for a body to tumble. Conservation of angular momentum fixes $\mathbf{L}$ in the fixed space frame, but the body can turn relative to $\mathbf{L}$, so its angular velocity components in the principal body frame change. The torque-free rigid-body invariants constrain this motion. A supported gyroscope with gravitational torque has different governing assumptions.

For $0<I_1<I_2<I_3$, the transverse linear stability rates about steady rotation with signed component $\Omega \ne 0$ areAxes one and three have transverse eigenvalues $\pm i{\nu }_1$ and $\pm i{\nu }_3$, producing bounded oscillations. Axis two has transverse eigenvalues $\pm {\lambda }_2$, giving growing and decaying exponential functions. A neutral longitudinal perturbation corresponds to a small change of $\Omega$. Generic intermediate-axis disturbances contain a growing component; a precisely chosen decaying direction is exceptional. For the normalized intermediate-axis model with $|\Omega |=1$, the values are ${\nu }_1=1/\sqrt{3}$, ${\nu }_3=1$ and ${\lambda }_2=1/\sqrt{3}$.

Proof outline: the linearization of the Euler equations for a torque-free rigid body about the third principal axis gives transverse perturbations obeyingand henceCyclic permutation gives the other rates. For nonlinear Lyapunov stability of the extreme rotations, use the torque-free rigid-body invariants: the conserved quadratic forms $L^2-2I_{1}E$ and $2I_{3}E-L^2$ have strictly positive transverse coefficients about axes one and three respectively. They bound transverse angular velocity, while $E$ bounds its longitudinal change near nonzero rotation. Imaginary transverse eigenvalues alone would not establish this nonlinear conclusion.^[ref]

The practical lesson is to choose a stable principal axis when using rotation to maintain an instrument's pointing direction. Both extreme principal moments of inertia are stable in the ideal torque-free model; if internal damping lowers rotational kinetic energy while approximately preserving angular momentum, the minimum energy $E=L^2/(2I_3)$ instead favors the largest principal moment of inertia. A control system must account for torque, changing mass distributions and flexibility before applying this ideal result.

A principal body frame is an orthonormal basis ${\mathbf{e}}_a(t)$ fixed in the rotating body and aligned with its principal axes. If $R(t)$ is the rotation matrix whose columns are these vectors, differentiating $R^{T}R=I$ shows that $\dot{R}{R}^T$ is skew-symmetric. Every skew-symmetric linear map in three dimensions is cross product with a unique vector $\mathbf{\omega }$, so

For any vector $\mathbf{A}=A_a{\mathbf{e}}_a$,Applying this to the conserved angular momentumgives the Euler equations for a torque-free rigid body:

Multiplying these equations respectively by ${\omega }_i$ and summing proves conservation of kinetic energy; multiplying by $I_i{\omega }_i$ and summing proves conservation of squared angular momentum. Thus the torque-free rigid-body invariants areSolving these two linear equations for the first two squared components givesUsing $I_3{\dot{\omega }}_3=(I_1-I_2){\omega }_1{\omega }_2$ yields the quartic equation for one component of torque-free angular velocity:

Multiply the three Euler equations for a torque-free rigid body by ${\omega }_1,{\omega }_2,{\omega }_3$ and add. The cubic terms cancel, giving conservation of the rotational kinetic energyInstead multiply the equations by $I_1{\omega }_1,I_2{\omega }_2,I_3{\omega }_3$ and add. The cubic terms again cancel, givingthe squared magnitude of the conserved angular momentum. For an asymmetric body these two quadratic forms are not proportional, so the torque-free rigid-body invariants are independent.

The rotational kinetic energy and squared angular momentum are the torque-free rigid-body invariantsandDifferentiating the energy and inserting the Euler equations givesSimilarly,