## Introduction

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.

## Derivation

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

and

The rotation matrix evolves due to angular velocity via

where is the anti-symmetric cross-product matrix generated by , giving

and

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.