Double pendulum

In physics and mathematics, in the area of dynamical systems, a double pendulum also known as a chaotic pendulum is a pendulum with another pendulum attached to its end, forming a simple physical system that exhibits rich dynamic behavior with a strong sensitivity to initial conditions.^[1] The motion of a double pendulum is governed by a set of coupled ordinary differential equations and is chaotic.

Analysis and interpretation

Several variants of the double pendulum may be considered; the two limbs may be of equal or unequal lengths and masses, they may be simple pendulums or compound pendulums (also called complex pendulums) and the motion may be in three dimensions or restricted to the vertical plane. In the following analysis, the limbs are taken to be identical compound pendulums of length $ℓ$ and mass $m$ , and the motion is restricted to two dimensions.

Double compound pendulum

Motion of the double compound pendulum (from numerical integration of the equations of motion)

In a compound pendulum, the mass is distributed along its length. If the double pendulum mass is evenly distributed, then the center of mass of each limb is at its midpoint, and the limb has a moment of inertia of $I = .mw-parser-output .sfrac{white-space:nowrap}.mw-parser-output .sfrac.tion,.mw-parser-output .sfrac .tion{display:inline-block;vertical-align:-0.5em;font-size:85%;text-align:center}.mw-parser-output .sfrac .num{display:block;line-height:1em;margin:0.0em 0.1em;border-bottom:1px solid}.mw-parser-output .sfrac .den{display:block;line-height:1em;margin:0.1em 0.1em}.mw-parser-output .sr-only{border:0;clip:rect(0,0,0,0);clip-path:polygon(0px 0px,0px 0px,0px 0px);height:1px;margin:-1px;overflow:hidden;padding:0;position:absolute;width:1px}⁠1/12⁠mℓ2$ about that point.

It is convenient to use the angles between each limb and the vertical as the generalized coordinates defining the configuration of the system. These angles are denoted $θ しーた 1$ and $θ しーた 2$ . The position of the center of mass of each rod may be written in terms of these two coordinates. If the origin of the Cartesian coordinate system is taken to be at the point of suspension of the first pendulum, then the center of mass of this pendulum is at: ${\begin{aligned}x_{1}&={\tfrac {1}{2}}\ell \sin \theta _{1}\\y_{1}&=-{\tfrac {1}{2}}\ell \cos \theta _{1}\end{aligned}}$

and the center of mass of the second pendulum is at ${\begin{aligned}x_{2}&=\ell \left(\sin \theta _{1}+{\tfrac {1}{2}}\sin \theta _{2}\right)\\y_{2}&=-\ell \left(\cos \theta _{1}+{\tfrac {1}{2}}\cos \theta _{2}\right)\end{aligned}}$ This is enough information to write out the Lagrangian.

Lagrangian

The Lagrangian is ${\begin{aligned}L&={\text{kinetic energy}}-{\text{potential energy}}\\&={\tfrac {1}{2}}m\left(v_{1}^{2}+v_{2}^{2}\right)+{\tfrac {1}{2}}I\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-mg\left(y_{1}+y_{2}\right)\\&={\tfrac {1}{2}}m\left({\dot {x}}_{1}^{2}+{\dot {y}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {y}}_{2}^{2}\right)+{\tfrac {1}{2}}I\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-mg\left(y_{1}+y_{2}\right)\end{aligned}}$ The first term is the linear kinetic energy of the center of mass of the bodies and the second term is the rotational kinetic energy around the center of mass of each rod. The last term is the potential energy of the bodies in a uniform gravitational field. The dot-notation indicates the time derivative of the variable in question.

Since (see Chain Rule and List of trigonometric identities) ${\begin{aligned}{\dot {x}}_{1}&={\dot {\theta }}_{1}\left({\tfrac {1}{2}}\ell \cos \theta _{1}\right)&&\rightarrow &{\dot {x}}_{1}^{2}&={\dot {\theta }}_{1}^{2}\left({\tfrac {1}{4}}\ell ^{2}\cos ^{2}\theta _{1}\right)\\[1ex]{\dot {y}}_{1}&={\dot {\theta }}_{1}\left({\tfrac {1}{2}}\ell \sin \theta _{1}\right)&&\rightarrow &{\dot {y}}_{1}^{2}&={\dot {\theta }}_{1}^{2}\left({\tfrac {1}{4}}\ell ^{2}\sin ^{2}\theta _{1}\right)\end{aligned}}$

{\dot {x}}_{1}^{2}+{\dot {y}}_{1}^{2}={\tfrac {1}{4}}{\dot {\theta }}_{1}^{2}\ell ^{2}\left(\cos ^{2}\theta _{1}+\sin ^{2}\theta _{1}\right)={\tfrac {1}{4}}\ell ^{2}{\dot {\theta }}_{1}^{2},

and

${\begin{aligned}{\dot {x}}_{2}&=\ell \left({\dot {\theta }}_{1}\cos \theta _{1}+{\tfrac {1}{2}}{\dot {\theta }}_{2}\cos \theta _{2}\right)&&\rightarrow &{\dot {x}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\cos ^{2}\theta _{1}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \theta _{1}\cos \theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\cos ^{2}\theta _{2}\right)\\{\dot {y}}_{2}&=\ell \left({\dot {\theta }}_{1}\sin \theta _{1}+{\tfrac {1}{2}}{\dot {\theta }}_{2}\sin \theta _{2}\right)&&\rightarrow &{\dot {y}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\sin ^{2}\theta _{1}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\sin \theta _{1}\sin \theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\sin ^{2}\theta _{2}\right)\end{aligned}}$ ${\begin{aligned}{\dot {x}}_{2}^{2}+{\dot {y}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\cos ^{2}\theta _{1}+{\dot {\theta }}_{1}^{2}\sin ^{2}\theta _{1}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\cos ^{2}\theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\sin ^{2}\theta _{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \theta _{1}\cos \theta _{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\sin \theta _{1}\sin \theta _{2}\right)\\[1ex]&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \left(\theta _{1}-\theta _{2}\right)\right),\end{aligned}}$

substituting the coordinates above and rearranging the equation gives ${\begin{aligned}L&={\tfrac {1}{2}}m\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \left(\theta _{1}-\theta _{2}\right)\right)+{\tfrac {1}{24}}m\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-mg\left(y_{1}+y_{2}\right)\\[1ex]&={\tfrac {1}{6}}m\ell ^{2}\left({\dot {\theta }}_{2}^{2}+4{\dot {\theta }}_{1}^{2}+3{{\dot {\theta }}_{1}}{{\dot {\theta }}_{2}}\cos(\theta _{1}-\theta _{2})\right)+{\tfrac {1}{2}}mg\ell \left(3\cos \theta _{1}+\cos \theta _{2}\right).\end{aligned}}$

The Euler-Lagrange equations then give the two following second-order, non-linear differential equations in $(\theta _{1},\theta _{2})$ ^[a]:

${\begin{aligned}m_{2}\left(g\sin(\theta _{1})+\ell _{2}\left({\dot {\theta }}_{2}^{2}\sin(\theta _{1}-\theta _{2})+{\ddot {\theta }}_{2}\cos(\theta _{1}-\theta _{2})\right)+\ell _{1}{\ddot {\theta }}_{1}\right)+m_{1}\left(g\sin(\theta _{1})+\ell _{1}{\ddot {\theta }}_{1}\right)&=0\\g\sin(\theta _{2})+\ell _{1}\left({\ddot {\theta }}_{1}\cos(\theta _{1}-\theta _{2})-{\dot {\theta }}_{1}^{2}\sin(\theta _{1}-\theta _{2})\right)+\ell _{2}{\ddot {\theta }}_{2}&=0\\\end{aligned}}$ No closed form solutions for $θ しーた 1$ and $θ しーた 2$ as functions of time are known, therefore solving the system can only be done numerically, using the Runge Kutta method or similar techniques.

Parametric plot for the time evolution of the angles of a double pendulum. It can be noticed that the graph resembles a Brownian motion.

Chaotic motion

Graph of the time for the pendulum to flip over as a function of initial conditions

Long exposure of double pendulum exhibiting chaotic motion (tracked with an LED)

The double pendulum undergoes chaotic motion, and clearly 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 position when released at rest. Here, the initial value of $θ しーた 1$ ranges along the $x$ -direction from −3.14 to 3.14. The initial value $θ しーた 2$ ranges along the $y$ -direction, from −3.14 to 3.14. The color of each pixel indicates whether either pendulum flips within:

${\sqrt {\frac {\ell }{g}}}$ (black)
$10{\sqrt {\frac {\ell }{g}}}$ (red)
$100{\sqrt {\frac {\ell }{g}}}$ (green)
$1000{\sqrt {\frac {\ell }{g}}}$ (blue) or
$10000{\sqrt {\frac {\ell }{g}}}$ (purple).

Three double pendulums with near identical initial conditions diverge over time, displaying the chaotic nature of the system.

Initial conditions that do not lead to a flip within $10000{\sqrt {\frac {\ell }{g}}}$ 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.$

Within the region defined by this curve, that is if $3\cos \theta _{1}+\cos \theta _{2}>2,$ then it is energetically impossible for either pendulum to flip. Outside this region, the pendulum can flip, but it is a complex question to determine when it will flip. Similar behavior is observed for a double pendulum composed of two point masses rather than two rods with distributed mass.^[2]

The lack of a natural excitation frequency has led to the use of double pendulum systems in seismic resistance designs in buildings, where the building itself is the primary inverted pendulum, and a secondary mass is connected to complete the double pendulum.

Notes

^ The equations were obtained with the following Mathematica code:

Block[
 {m, g, \[Theta], l, L, x, y, v, t},
 x = Accumulate[{Subscript[l, 1], Subscript[l, 2]}* 
    Sin[{Subscript[\[Theta], 1], Subscript[\[Theta], 2]}]];
 y = Accumulate[ {Subscript[l, 1], Subscript[l, 
     2]}*-Cos[{Subscript[\[Theta], 1], Subscript[\[Theta], 2]}]];
 v = D[{x, y} /. Subscript[\[Theta], i_] :> Subscript[\[Theta], i][t],
    t];
 L = Plus @@ ({Subscript[m, 1], Subscript[m, 
      2]}*(1/2 Map[# . # &, Transpose[v]] - ( 
        g y /. Subscript[\[Theta], i_] :> Subscript[\[Theta], i][t])));
 FullSimplify[
   Table[
    D[
      Construct[
         Function,
         L /. Subscript[\[Theta], i]'[t] -> #
         ]'[Subscript[\[Theta], i]'[t]], t
      ] ==
     Construct[
        Function,
        L /. Subscript[\[Theta], i][t] -> #
        ]'[Subscript[\[Theta], i][t]]
    , {i, 2}
    ],
   Assumptions -> {
     Subscript[l, 1] > 0, Subscript[l, 2] > 0,
     Subscript[m, 1] > 0, Subscript[m, 2] > 0
     }
   ] /. h_[t] :> h // Column // TeXForm
 ]

References

^ Levien, R. B.; Tan, S. M. (1993). "Double Pendulum: An experiment in chaos". American Journal of Physics. 61 (11): 1038. Bibcode:1993AmJPh..61.1038L. doi:10.1119/1.17335.
^ Alex Small, Sample Final Project: One Signature of Chaos in the Double Pendulum, (2013). A report produced as an example for students. Includes a derivation of the equations of motion, and a comparison between the double pendulum with 2 point masses and the double pendulum with 2 rods.

Meirovitch, Leonard (1986). Elements of Vibration Analysis (2nd ed.). McGraw-Hill Science/Engineering/Math. ISBN 0-07-041342-8.
Eric W. Weisstein, Double pendulum (2005), ScienceWorld (contains details of the complicated equations involved) and "Double Pendulum" by Rob Morris, Wolfram Demonstrations Project, 2007 (animations of those equations).
Peter Lynch, Double Pendulum, (2001). (Java applet simulation.)
Northwestern University, Double Pendulum Archived 2007-06-03 at the Wayback Machine, (Java applet simulation.)
Theoretical High-Energy Astrophysics Group at UBC, Double pendulum, (2005).

External links

Animations and explanations of a double pendulum and a physical double pendulum (two square plates) by Mike Wheatland (Univ. Sydney)

Interactive Open Source Physics JavaScript simulation with detailed equations double pendulum
Interactive Javascript simulation of a double pendulum
Double pendulum physics simulation from www.myphysicslab.com using open source JavaScript code
Simulation, equations and explanation of Rott's pendulum
Comparison videos of a double pendulum with the same initial starting conditions on YouTube
Double Pendulum Simulator - An open source simulator written in C++ using the Qt toolkit.
Online Java simulator Archived 2022-08-16 at the Wayback Machine of the Imaginary exhibition.

[2] The equations were obtained with the following Mathematica code:
Block[ {m, g, \[Theta], l, L, x, y, v, t}, x = Accumulate[{Subscript[l, 1], Subscript[l, 2]}* Sin[{Subscript[\[Theta], 1], Subscript[\[Theta], 2]}]]; y = Accumulate[ {Subscript[l, 1], Subscript[l, 2]}*-Cos[{Subscript[\[Theta], 1], Subscript[\[Theta], 2]}]]; v = D[{x, y} /. Subscript[\[Theta], i_] :> Subscript[\[Theta], i][t], t]; L = Plus @@ ({Subscript[m, 1], Subscript[m, 2]}*(1/2 Map[# . # &, Transpose[v]] - ( g y /. Subscript[\[Theta], i_] :> Subscript[\[Theta], i][t]))); FullSimplify[ Table[ D[ Construct[ Function, L /. Subscript[\[Theta], i]'[t] -> # ]'[Subscript[\[Theta], i]'[t]], t ] == Construct[ Function, L /. Subscript[\[Theta], i][t] -> # ]'[Subscript[\[Theta], i][t]] , {i, 2} ], Assumptions -> { Subscript[l, 1] > 0, Subscript[l, 2] > 0, Subscript[m, 1] > 0, Subscript[m, 2] > 0 } ] /. h_[t] :> h // Column // TeXForm ]

[1] Levien, R. B.; Tan, S. M. (1993). "Double Pendulum: An experiment in chaos". American Journal of Physics. 61 (11): 1038. Bibcode:1993AmJPh..61.1038L. doi:10.1119/1.17335.

[3] Alex Small, Sample Final Project: One Signature of Chaos in the Double Pendulum, (2013). A report produced as an example for students. Includes a derivation of the equations of motion, and a comparison between the double pendulum with 2 point masses and the double pendulum with 2 rods.

[1]

[a]

[2]