Double pendulum
Encyclopedia : D : DO : DOU : Double pendulum
In horology, a double pendulum is a system of two simple pendulums on a common mounting which move in anti-phase.
In mathematics, in the area of dynamical systems, a double pendulum is a pendulum with another pendulum attached to its end, and is a simple physical system that exhibits rich dynamic behavior. The motion of a double pendulum is governed by a set of coupled ordinary differential equations. Above a certain energy its motion is chaotic.
The double pendulum consists of two thin rods (moment of inertia, [I=\frac M \ell^2]) connected by a pivot and the end of one rod suspended from a pivot. It is natural to define the coordinates to be the angle between each rod and the vertical. These are denoted by θ1 and θ2. The position of the centre of mass of the two rods may be written in terms of these coordinates. These are given by
- [x_1 = \frac \sin \theta_1,]
- [x_2 = \ell \left ( \sin \theta_1 + \frac \sin \theta_2 \right ),]
- [y_1 = -\frac \cos \theta_1]
- [y_2 = -\ell \left ( \cos \theta_1 + \frac \cos \theta_2 \right ).]
Lagrangian
The Lagrangian is given by
- [L = \frac m \left ( v_1^2 + v_2^2 \right ) + \frac I \left ( ^2 + ^2 \right ) - m g \left ( y_1 + y_2 \right ) ]
Plugging in the coordinates above and doing a bit of algebra gives
- [L = \frac m \ell^2 \left [ ^2 + 4 ^2 + 3 cos (theta_1-theta_2) right ] + \frac m g \ell \left ( 3 \cos \theta_1 + \cos \theta_2 \right ).]
- [p_ = \frac} = \frac m \ell^2 \left [ 8 + 3 cos (theta_1-theta_2) right ]]
- [p_ = \frac} = \frac m \ell^2 \left [ 2 + 3 cos (theta_1-theta_2) right ].]
- [ = \frac \frac - 3 \cos(\theta_1-\theta_2) p_}]
- [ = \frac \frac - 3 \cos(\theta_1-\theta_2) p_}.]
- [} = \frac = -\frac m \ell^2 \left [ sin (theta_1-theta_2) + 3 frac sin theta_1 right ]]
- [} = \frac = -\frac m \ell^2 \left [ - sin (theta_1-theta_2) + frac sin theta_2 right ].]
Chaotic motion
The double pendulum undergoes chaotic motion, and shows a sensitive dependence on initial conditions. The image to the right shows the amount of elapsed time before the pendulum "flips over", as a function of initial conditions. Here, the initial value of θ1 ranges along the x-direction, from −3 to 3. The initial value θ2 ranges along the y-direction, from −3 to 3. The colour of each pixel indicates whether either pendulum flips within [10\sqrt] (green), within 100 (red), 1000 (purple) or 10000 (blue). Initial conditions that don't lead to a flip within [10000\sqrt] are plotted white.
The boundary of the central white region is defined in part by energy conservation with the following curve:
- [3 \cos \theta_1 + \cos \theta_2 = 2. \,]
- [3 \cos \theta_1 + \cos \theta_2 > 2, \,]
References
- Eric W. Weisstein, [Double pendulum] (2005), ScienceWorld. (Contains details of the complicated equations involved.)
- Peter Lynch, [Double Pendulum], (2001). (Java applet simulation.)
- Theoretical High-Energy Astrophysics Group at UBC, [Double pendulum], (2005).
From Wikipedia, the Free Encyclopedia. Original article here. Support Wikipedia by contributing or donating.
All text is available under the terms of the GNU Free Documentation License See Wikipedia Copyrights for details.

