Double Pendulum

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 Equation s. Above a certain Energy its motion is Chaotic .

The double pendulum consists of two thin rods ( Moment Of Inertia ,

I=rac{1}{12} M l^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 θ₁ and θ₂. 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 = rac{l}{2} \sin 	heta_1,

x_2 = l \left (  \sin 	heta_1 + rac{1}{2} \sin 	heta_2 
ight ),

y_1 = -rac{l}{2} \cos 	heta_1

and
:

y_2 = -l \left (  \cos 	heta_1 + rac{1}{2} \cos 	heta_2 
ight ).

This is enough information to write out the Lagrangian.

Lagrangian

The Lagrangian is given by
:

L = rac{1}{2} m \left ( v_1^2 + v_2^2 
ight ) + rac{1}{2} I \left ( {\dot 	heta_1}^2 + 
{\dot 	heta_2}^2 
ight ) - m g \left ( y_1 + y_2 
ight )

where the Kinetic Energy is the sum of the kinetic energy of the Center Of Mass of each rod and the kinetic energy about the centres of mass of the rods. The Potential Energy of a body in a uniform gravitational field is given by the potential energy at the center of mass.

Plugging in the coordinates above and doing a bit of algebra gives
:

+ rac{1}{2} m g l \left ( 3 \cos heta_1 + \cos heta_2 ight )

There is only one conserved quantity (the energy), and no conserved Momenta . The two momenta may be written as

:

p_{	heta_1} = rac{\partial L}{\partial {\dot 	heta_1}} = rac{1}{6} m l^2 \left 8 {\dot 	heta_1}  + 3 {\dot 	heta_2} \cos (	heta_1-	heta_2) 
ight

and
:

p_{	heta_2} = rac{\partial L}{\partial {\dot 	heta_2}} = rac{1}{6} m l^2 \left 2 {\dot 	heta_2} + 3 {\dot 	heta_1} \cos (	heta_1-	heta_2)  
ight

These expressions may be inverted to get

:

{\dot 	heta_1} = rac{6}{ml^2} rac{ 2 p_{	heta_1} - 3 \cos(	heta_1-	heta_2) p_{	heta_2}}{16 - 9 \cos^2(	heta_1-	heta_2)}

and
:

{\dot 	heta_2} = rac{6}{ml^2} rac{ 8 p_{	heta_2} - 3 \cos(	heta_1-	heta_2) p_{	heta_1}}{16 - 9 \cos^2(	heta_1-	heta_2)}

The remaining equations of motion are written as

:

{\dot p_{	heta_1}} = rac{\partial L}{\partial 	heta_1} = -rac{1}{2} m l^2 \left {\dot 	heta_1} {\dot 	heta_2} \sin (	heta_1-	heta_2) + 3 rac{g}{l} \sin 	heta_1 
ight

and

:

{\dot p_{	heta_2}} = rac{\partial L}{\partial 	heta_2}
 = -rac{1}{2} m l^2 \left -{\dot 	heta_1} {\dot 	heta_2} \sin (	heta_1-	heta_2) +  rac{g}{l} \sin 	heta_2 
ight

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 θ₁ ranges along the x-direction, from -3 to 3. The initial value θ₂ ranges along the y-direction, from -3 to 3. The colour of each pixel indicates whether either pendulum flips within 10

\sqrt{g/l}

(green), within 100 (red), 1000 (purple) or 10000 (blue). Initial conditions that don't lead to a flip within 10000

\sqrt{g/l}

are plotted white.

The boundary of the central white region is defined in part by energy conservation with the following curve

3 \cos 	heta_1 + \cos 	heta_2  = 2 \,

Within the region defined by this curve, that is if

3 \cos 	heta_1 + \cos 	heta_2  > 2 \,

it is energetically impossible for either pendulum to flip. Outside this region, the pendulums can flip but this is different from determining when they will flip.

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).

Information About

Double Pendulum