Double Pendulum Explorer

explore deterministic chaos and the butterfly effect in coupled motion

Loading interactive simulation...

why double pendulums become unpredictable 🖖

This is deterministic chaos, not randomness. Start two runs just 0.1Β° apart: after a while they no longer look related because separation grows roughly exponentially (a positive Lyapunov exponent). Energy keeps shuttling between the two arms through nonlinear coupling, so tiny phase differences get amplified into completely different future states. The equations are exact; what collapses is the prediction horizon.
This pendulum does not just go left-right; it keeps changing its mind about which arm leads the dance.

just two rods and gravity 🖖

There is no motor, no noise generator, no trick β€” only two arms, two pivots, and gravity. The wild motion comes entirely from the second joint, which lets each arm push the other around. Try starting with a small nudge near the bottom: it swings almost like an ordinary pendulum. Chaos only takes over once you give it enough energy for the arms to swing over the top.

the tip can fall faster than gravity 🖖

Release the pendulum from rest with the arms held out horizontally and watch the outer bob: for an instant it can accelerate downward faster than a freely dropped ball β€” more than g. The rigid arms act as a lever and whip, so the inner arm flings the outer one down harder than gravity alone ever could. It is the same effect that lets a hinged stick beat a falling coin to the floor.
Gravity starts the fall, but geometry decides who lands first.

Problem solved in full

  1. The period of two clean oscillations from release at rest 7 steps

    Released from rest at θ₁ = 20Β° and ΞΈβ‚‚ = 10Β°, with m₁ = mβ‚‚ = 1, l₁ = lβ‚‚ = 1 and g = 9.81, and damping c = 0. Angles this small are not chaotic β€” the motion is a sum of two clean oscillations. Find the period of each, then decide which of them the cycle marks on the panel are counting.

    1. Cut the equations of motion above down to size. Small angles mean sin ΞΈ becomes ΞΈ and cos 2Ξ” becomes 1, and the squared-velocity terms are second order in a small amplitude, so they go. Both denominators collapse to 2m₁ + mβ‚‚ βˆ’ mβ‚‚ = 2, and what survives is linear and still coupled: neither bob moves without pushing the other.

    2. Set a yardstick before solving anything. A single 1 m pendulum swings at Ο‰β‚€ = √(g/l), giving the 2.01 s the panel prints as its period hint and stamps on every cycle mark. Hold on to it, because nothing below is going to equal it.

    3. A normal mode is a motion in which both angles oscillate at one shared frequency and keep a fixed ratio. Substitute ΞΈα΅’ = aα΅’ cos Ο‰t and each of the two equations hands you that ratio aβ‚‚/a₁, arrived at from opposite ends.

    4. The two expressions for one ratio must agree, and setting them equal eliminates both amplitudes. What is left is a quadratic in ω², and a quadratic has 2 roots. That is where the second period comes from.

    5. So this pendulum has 2 periods, 2.62 s and 1.09 s, and the 2.01 s on the cycle marks is neither. It lies between them and belongs to neither one.

    6. Each root carries its own shape back. The slow root gives aβ‚‚/a₁ = +1.414: both bobs swing the same way, the lower one 1.414 times wider. The fast root gives βˆ’1.414: they swing in opposition and the joint between them scissors.

    7. A start of 20Β° and 10Β° is neither shape, so both modes are running at once. Split it into a slow amplitude A and a fast amplitude B β€” the top bob needs A + B = 20Β°, the bottom needs √2(A βˆ’ B) = 10Β° β€” and the slow mode takes 13.54Β° of the start while the fast one takes 6.46Β°.

    Answer

    2.62 s and 1.09 s, and their ratio is exactly 1 + √2. An irrational ratio means the 2 modes never come back into step, so the figure the lower bob traces never quite repeats β€” and this is the tame small-angle case, with no chaos in it anywhere. It also settles what the 2.01 s cycle marks are: a single-pendulum ruler laid against a 2-frequency motion, which is why the angles logged at each mark keep drifting instead of returning to 20.0Β° and 10.0Β°. Push the release to 145Β° and even the 2 frequencies go, because sin 145Β° = 0.5736 while the angle itself is 2.531 rad β€” wrong by a factor of 4.4, and the linearisation this whole derivation stands on has nothing left to say.

References (1)

Example problems

  • small angles - Small angle regime is more regular and less sensitive.
  • chaotic kick - Large initial offset produces fast divergence in phase.
  • heavy second bob - Longer second arm with slight damping highlights energy decay.