I’ve been seeing the relationship between torque and angular acceleration stated as
Although this has a nice analogy with the force, mass and acceleration of translational dynamics (), it is actually incorrect. There is an extra term missing from the right-hand side. The correct version is
This extra term is sometimes called the gyroscopic term. It is non-zero when the angular velocity and angular momentum have different directions. Unfortunately there are some lecture notes, slides and articles that make no mention of this extra term, nor do they explain why it has been omitted.
This is particularly the case with material aimed at the development of real-time simulations of rigid-body dynamics, such as game development.
Looking at open source physics libraries, you most likely won’t find this missing term there either, so where has it gone?
The reason for the omission is that the extra term causes instability in numerical simulations with general-purpose solvers. There’s no physical reason to omit it — indeed the extra term has real physical consequences — but where physical accuracy isn’t paramount it is simpler to just remove the term than to employ a special-purpose solver.
This article will derive the torque-angular acceleration relationship in full.
From basic rigid-body dynamics we know that the net torque, , is the time derivative of the angular momentum, ,
The moment of inertia tensor, , is defined as the linear relationship between angular velocity and angular momentum,
Combining these two equations gives
The key thing to note is that the moment of inertia tensor is a time-dependent quantity — it changes as the body rotates. Only when described in a body-fixed rotating frame is it constant.
Now is just the definition of angular acceleration, , so this becomes
The tricky bit is the second term, which involves the time derivative of the inertia tensor. If is the constant inertia tensor in the body-fixed frame and is the current rotation matrix for the body (the body-to-world rotation), then
The rotation matrix evolves due to angular velocity via
where is the anti-symmetric cross-product matrix generated by , giving
The second term is null,
and the remaining term is
Therefore the torque is
and the angular acceleration is
The Gyroscopic Term
In the absence of torque the gyroscopic term means the angular acceleration can still be non-zero. This is the cause of a phenomenum called torque-free precession, where the angular velocity vector gradually rotates around the angular momentum.
The reason for the gyroscopic term’s mysterious absence from numerical simulations is that it can cause instability when used with general-purpose solvers. A full treatment of this instability is beyond the scope of this short article, but suffice it to say that numerical instability is not acceptable. A special solver may be needed to safely include this term, and so as long as a physically accurate result isn’t required, it is simpler to omit it completely. This is a simple and cheap solution to the problem, but it means that the simulation isn’t accurate enough for the scientific level as it cannot exhibit phenomena such as torque-free precession. However, for non-scientific applications, such as games, it can still give a physically-plausible result without risk of instabilities nor cost of special solvers.