Double Pendulum
Introduction
A double pendulum is one of the simplest physical systems that exhibits chaos. It consists of a first pendulum arm hanging from a fixed pivot, with a second pendulum arm hanging from the end of the first. Under gravity, this arrangement produces motion that is extraordinarily sensitive to its starting conditions — two double pendulums released from positions that differ by even a fraction of a degree will, within seconds, trace out completely different paths. This behaviour is not the result of randomness or noise; it is a deterministic consequence of the nonlinear equations that govern the system. The double pendulum sits at the crossroads of classical mechanics, oscillation theory, and chaos theory, making it one of the most instructive and visually striking systems in all of physics.
The Physics Explained
A simple single pendulum swings back and forth in a regular, repeating pattern. For small angles the motion is nearly sinusoidal, and the period depends only on the length of the arm and the strength of gravity — not on the mass or the amplitude. This is why pendulums were used in clocks for centuries. The single pendulum is a well-behaved, predictable oscillator.
The double pendulum breaks that predictability. Because the second arm is attached to a moving point rather than a fixed pivot, the two arms interact with each other through their mutual forces. The motion of the upper arm changes the effective gravity experienced by the lower arm, and the swinging of the lower arm pulls back on the upper arm through tension in the connecting rod. These interactions are nonlinear — the forces depend on the angles in a way that cannot be simplified into a clean proportional relationship — and it is this nonlinearity that produces chaos.
The correct framework for analysing the double pendulum is Lagrangian mechanics, a reformulation of Newtonian mechanics that works in terms of energy rather than forces. Instead of resolving all the tensions and reactions at the pivot points, we write down the total kinetic energy and the total potential energy of the system, then apply a mathematical procedure called the Euler-Lagrange equations to derive the equations of motion. The result is a pair of coupled, nonlinear, second-order differential equations — one for each angle. These equations have no closed-form analytical solution; to find how the angles change over time, they must be integrated numerically.
Chaos in the double pendulum means that tiny differences in initial conditions grow exponentially over time, a property quantified by the Lyapunov exponent. If two double pendulums start with angles that differ by a tiny amount, the difference between their trajectories grows roughly as e raised to the power of the Lyapunov exponent multiplied by time. This exponential divergence is the mathematical signature of chaos, and it places an absolute limit on how far into the future the motion can be predicted in practice. Despite this, the system conserves energy perfectly — no chaos here means randomness in energy, only in trajectory.
For small oscillations about the downward equilibrium position, the nonlinear terms are negligible and the system behaves like two coupled linear oscillators. In this regime there are two normal modes: one in which both arms swing in the same direction at the same frequency, and one in which the arms swing in opposite directions at a higher frequency. As the amplitude grows, the nonlinear terms kick in and the motion transitions from regular oscillation to chaos.
Key Equations
Key Variables
| Symbol | Unit | Description |
|---|---|---|
| theta1 | rad | Angle of the upper arm measured from the downward vertical |
| theta2 | rad | Angle of the lower arm measured from the downward vertical |
| theta1_dot | rad/s | Angular velocity of the upper arm |
| theta2_dot | rad/s | Angular velocity of the lower arm |
| theta1_ddot | rad/s^2 | Angular acceleration of the upper arm |
| theta2_ddot | rad/s^2 | Angular acceleration of the lower arm |
| m1 | kg | Mass of the bob at the end of the upper arm |
| m2 | kg | Mass of the bob at the end of the lower arm |
| L1 | m | Length of the upper arm |
| L2 | m | Length of the lower arm |
| g | m/s^2 | Acceleration due to gravity (9.81 m/s^2 at Earth's surface) |
| KE | J | Total kinetic energy of both arms and bobs |
| PE | J | Total gravitational potential energy of both bobs |
| E | J | Total mechanical energy; conserved throughout the motion |
| lambda | 1/s | Lyapunov exponent; measures the rate of divergence of nearby trajectories |
Real World Examples
- Human arm and forearm: The shoulder and elbow joints form a biological double pendulum. When you throw a ball or swing a tennis racket, your nervous system must continuously compensate for the chaotic tendency of the coupled arm segments, which is part of why fine motor control requires so much neural processing power.
- Robotics and prosthetics: Engineers designing robotic arms face the double pendulum problem directly. The coupled nonlinear dynamics of multi-joint limbs must be handled by sophisticated control algorithms that constantly correct for the system's sensitivity to disturbances.
- Acrobatics and gymnastics: A gymnast swinging on the high bar and releasing into a somersault behaves approximately like a double pendulum. Coaches and athletes develop intuition for the nonlinear coupling between the torso and legs, though the mathematics underlying their movements is the same chaotic system studied in physics labs.
- Molecular dynamics: Chains of atoms or polymer molecules can exhibit double-pendulum-like dynamics at the molecular scale. Understanding how coupled oscillators behave chaotically informs models of molecular flexibility and protein folding.
- Chaos demonstrations in education: The double pendulum is one of the most widely used physical demonstrations in university physics courses precisely because its chaos is visible to the naked eye. Two identical pendulums released from very slightly different angles diverge within moments, making the abstract concept of sensitive dependence on initial conditions immediately tangible.
How the Simulation Works
The simulation lets you set the length and mass of each arm using sliders, and choose the starting angles of both the upper and lower arms. When you press the play button, the simulation integrates the two coupled equations of motion forward in time using a fourth-order Runge-Kutta numerical integration scheme. This method estimates the rate of change of all four state variables — theta1, theta2, theta1_dot, and theta2_dot — at several intermediate points within each time step, then combines those estimates to produce a highly accurate update. The step size is kept small enough that the total energy of the system remains very nearly constant throughout the run, confirming that the integration is accurate.
The positions of the two bobs are computed from the angles using standard trigonometry: the upper bob sits at (L1*sin(theta1), -L1*cos(theta1)) relative to the pivot, and the lower bob sits at that position plus (L2*sin(theta2), -L2*cos(theta2)). These coordinates are drawn to the canvas on every frame, with an optional trail showing the path traced by the lower bob. The trail makes the chaotic, space-filling nature of the trajectory immediately visible.
To demonstrate chaos directly, a second ghost pendulum can be enabled. This second pendulum starts with an initial angle that differs from the first by a very small amount — far less than anything visible on screen at the start. You can watch as the two pendulums quickly diverge and begin tracing entirely different paths, illustrating sensitive dependence on initial conditions without any hand-waving. The energy readout confirms that both pendulums conserve their total mechanical energy throughout, showing that chaos and energy conservation coexist in this deterministic system.
Further Reading
- Lagrangian mechanics — the energy-based approach to deriving equations of motion for complex systems with constraints
- Chaos theory and the Lyapunov exponent — the mathematical tools used to quantify and characterise chaotic behaviour
- Simple harmonic motion and the simple pendulum — the linear oscillator that the double pendulum reduces to for small angles
- Runge-Kutta numerical integration — the algorithm used to solve differential equations that have no analytical solution
- Normal modes and coupled oscillators — the regular, non-chaotic behaviour that emerges when oscillation amplitudes are kept small
- Phase space and Poincare sections — graphical tools for visualising and understanding the structure of chaotic trajectories