At a glance
This sim spins a free box, with no gravity and no torque, about its length, width or thickness, adds a small disturbance, and follows what the spin does next. The box's shape decides the outcome, which matters for a thrown book, a tennis racket and a spacecraft.
What you'll learn
- What the three principal axes and principal moments of a box are.
- How the moments of inertia decide what a small disturbance does about each axis.
- How the time before the first turn depends on the size of the disturbance.
- Which quantities stay constant while the box moves, and which ones change.
What's on screen
- The box, with the face colours above and a line through each axis (the chosen spin axis is thicker).
- L (amber arrow): the angular momentum.
- ω (rose arrow): the angular velocity.
- Body-frame inset (bottom right; in the Guided Inquiry it appears at step 4): the direction of L seen from the box, with a trail.
- Readouts: time since release, kinetic energy T and |L|, plus (in Free Exploration) the measured first turn and an estimate of it.
The controls
Press ⌖ on a row to see where that control sits in the sim. The disturbance and the axis are the two controls that set up a run, and both restart it. The rest change only what you see or when it plays.
| Control | What it changes | Range | |
|---|---|---|---|
| Disturbance ε | The size of the starting disturbance: a small extra spin component εω₀ added to the main spin ω₀. Changing it restarts the run. | 0.002–0.1 opens at 0.020 | |
| Axis of spin | Which principal axis carries the main spin: length e₁ (smallest moment), width e₂ (middle), thickness e₃ (largest). Choosing one restarts the run. | e₁, e₂ or e₃ opens at e₂ in Free Exploration | |
| Readouts | Time since release, T and |L|. Free Exploration also shows the first turn (measured) and the estimate λ⁻¹ln(1/ε) (section 05). | panel | |
| Body-frame inset | The tip of L̂, the unit vector along L, drawn on a sphere fixed to the box, with a blue trail. | stage, bottom right | |
| The stage | Drag it to orbit the view. This changes the viewpoint only, not the motion. | drag | |
| ▶ Play / ⏸ Pause | Starts or freezes the run. | top bar | |
| Speed | Plays the motion faster or slower on screen. The physics stays the same. | 0.25×–4× opens at 1× | |
| ↻ Reset | Restarts from t = 0 with ε back at 0.020 and the axis back at thickness e₃. | top bar |
Fixed: m = 1 kg; box 0.30 × 0.20 × 0.10 m; main spin ω₀ = 2π rad/s, one turn per second of the t readout. At 1× the sim plays 0.4 s of motion per real second.
Three axes, three moments
Every rigid body has three principal axes, at right angles through its centre of mass. Spin the body exactly about one of them and ω and L both lie along it, so the spin does not change. For a uniform box they run through the centres of the faces. Their moments of inertia are:
With no torque, the angular velocity components obey Euler's equations:
They keep two quantities constant, whatever the starting spin: the kinetic energy T and the size |L| of the angular momentum.
One sign decides the stability
Spin about axis k and add a small disturbance. To first order, a small tilt q of the spin away from axis k obeys the equation below, where i and j are the other two axes:
| Spin axis | C (s⁻²) | A small disturbance |
|---|---|---|
| e₁, length (smallest I) | −12.1 | oscillates, period 1.80 s |
| e₂, width (middle I) | +9.11 | grows exponentially |
| e₃, thickness (largest I) | −18.9 | oscillates, period 1.44 s |
- The axes with the smallest and the largest moment are stable. The one in between is not.
- The sign of C comes from the ordering of the moments alone. The size of ε and the spin rate set only the sizes and the timing.
- The growth lasts only while the tilt is small. Sections 05 and 06 follow it further.
How long before the first turn
About the width, the growth rate comes from the box's own moments:
A tilt that starts at a fraction ε of the spin grows until it is of order one, and the estimate for that time is:
This calculator starts at the sim's ε. Moving it does not change the sim.
Making ε ten times smaller adds a fixed ln 10/λ = 0.76 s to the wait, rather than multiplying it by ten. The sim's measured first turn comes about 0.4 s after the estimate at every ε, because the estimate is built from the linearised equations and the sim integrates the full ones. At ε = 0.02 the estimate is 1.30 s and the measured turn is 1.71 s.
Why the turn repeats, and nothing is lost
Because T and |L| are fixed, the tip of L̂ in the box's frame can only move along a closed curve. For a width spin with a small ε, that curve passes close to +e₂ and to −e₂. L̂ lingers near one, swings across to the other (a turn of the box), lingers there, and swings back.
- The turns are regular. At ε = 0.02 they come at 1.71 s, 5.12 s and 8.54 s: the first after half a gap, then every 3.42 s.
- ε sets the timing. A smaller ε makes L̂ linger longer near each of +e₂ and −e₂, but the pattern of turns is the same. Only ε = 0 never turns.
- L does not move. No torque acts, so L is fixed in space. The box moves around it, and what changes is which face points along L.
- Nothing is lost. T and |L| keep the same six decimals through every turn.
Where it matters
- The Dzhanibekov effect. In 1985, aboard Salyut 7, cosmonaut Vladimir Dzhanibekov released a wing nut from its bolt. It spun steadily, then turned over, again and again at regular intervals. The nut was spinning about its middle-moment axis, the unstable one.
- Thrown objects. A tennis racket tossed with a spin about the axis across the head, in the plane of the strings, tends to turn over. The same goes for a phone or a book tossed about its width. Spins about the other two axes stay steady.
- Spacecraft. Spin-stabilised satellites depend on it. Explorer 1 (1958) was spun about its long, smallest-moment axis, and soon after launch its spin shifted to end-over-end rotation.
Go deeper: why Explorer 1 changed its spin
- The sim's body is perfectly rigid. A real body flexes, and the flexing turns rotational energy into heat. No torque acts, so |L| cannot change, but T slowly falls.
- For a fixed |L|, T = |L|²/2I along a principal axis, which is lowest for the largest I. Energy loss therefore drives the spin towards the largest-moment axis.
- So with energy loss the smallest-moment axis is not permanently stable. The sim does not model this.
Summary
Key points
- The sign of C = (Ij − Ik)(Ik − Ii)ω₀²/IiIj decides: negative gives an oscillation, positive gives exponential growth.
- The first turn comes after about ln(1/ε)/λ (the sim measures about 0.4 s more), and each tenfold cut in ε adds a fixed ln10/λ = 0.76 s.
- T and |L| are constant, which confines L̂ to a closed curve in the box, so the turns repeat.
Common mix-ups
- "Something must be pushing the box over." No torque acts. The turn comes from Euler's equations alone.
- "The smallest-moment axis must be the least stable." For a rigid body it is as stable as the largest. Only the middle axis is unstable.
- "A tiny enough disturbance means no turn." Any ε above zero turns. Smaller ε only delays it.
Things to try on your own
- Run the width spin at ε = 0.1 and at ε = 0.002, and compare each measured first turn with its estimate.
- Run the length and the thickness at ε = 0.1 and compare the loops the inset's trail makes around each axis.