Double pendulum

Last updated October 31, 2024

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 Double-compound-pendulum-dimensioned.svg — Double compound pendulum

Motion of the double compound pendulum (from numerical integration of the equations of motion) Double-compound-pendulum.gif — 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:

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})$ ^{[lower-alpha 1]}:

${\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. Double-pendulum.png — 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

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

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/2Map[#.#&,Transpose[v]]-(gy/.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]

Related Research Articles

In mathematics, the polar coordinate system is a two-dimensional coordinate system in which each point on a plane is determined by a distance from a reference point and an angle from a reference direction. The reference point is called the pole, and the ray from the pole in the reference direction is the polar axis. The distance from the pole is called the radial coordinate, radial distance or simply radius, and the angle is called the angular coordinate, polar angle, or azimuth. Angles in polar notation are generally expressed in either degrees or radians.

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a given point in space is specified by three real numbers: the radial distance $r$ along the radial line connecting the point to the fixed point of origin; the polar angle $θ$ between the radial line and a given polar axis; and the azimuthal angle $φ$ as the angle of rotation of the radial line around the polar axis. (See graphic re the "physics convention".) Once the radius is fixed, the three coordinates (r, θ, φ), known as a 3-tuple, provide a coordinate system on a sphere, typically called the spherical polar coordinates. The plane passing through the origin and perpendicular to the polar axis (where the polar angle is a right angle) is called the reference plane (sometimes fundamental plane).

In physics, equations of motion are equations that describe the behavior of a physical system in terms of its motion as a function of time. More specifically, the equations of motion describe the behavior of a physical system as a set of mathematical functions in terms of dynamic variables. These variables are usually spatial coordinates and time, but may include momentum components. The most general choice are generalized coordinates which can be any convenient variables characteristic of the physical system. The functions are defined in a Euclidean space in classical mechanics, but are replaced by curved spaces in relativity. If the dynamics of a system is known, the equations are the solutions for the differential equations describing the motion of the dynamics.

<span class="mw-page-title-main">Legendre polynomials</span> System of complete and orthogonal polynomials

In mathematics, Legendre polynomials, named after Adrien-Marie Legendre (1782), are a system of complete and orthogonal polynomials with a wide number of mathematical properties and numerous applications. They can be defined in many ways, and the various definitions highlight different aspects as well as suggest generalizations and connections to different mathematical structures and physical and numerical applications.

<span class="mw-page-title-main">Tautochrone curve</span> Curve for which the time to roll to the end is equal for all starting points

A tautochrone curve or isochrone curve is the curve for which the time taken by an object sliding without friction in uniform gravity to its lowest point is independent of its starting point on the curve. The curve is a cycloid, and the time is equal to π times the square root of the radius over the acceleration of gravity. The tautochrone curve is related to the brachistochrone curve, which is also a cycloid.

In mathematics and physical science, spherical harmonics are special functions defined on the surface of a sphere. They are often employed in solving partial differential equations in many scientific fields. The table of spherical harmonics contains a list of common spherical harmonics.

An inverted pendulum is a pendulum that has its center of mass above its pivot point. It is unstable and falls over without additional help. It can be suspended stably in this inverted position by using a control system to monitor the angle of the pole and move the pivot point horizontally back under the center of mass when it starts to fall over, keeping it balanced. The inverted pendulum is a classic problem in dynamics and control theory and is used as a benchmark for testing control strategies. It is often implemented with the pivot point mounted on a cart that can move horizontally under control of an electronic servo system as shown in the photo; this is called a cart and pole apparatus. Most applications limit the pendulum to 1 degree of freedom by affixing the pole to an axis of rotation. Whereas a normal pendulum is stable when hanging downward, an inverted pendulum is inherently unstable, and must be actively balanced in order to remain upright; this can be done either by applying a torque at the pivot point, by moving the pivot point horizontally as part of a feedback system, changing the rate of rotation of a mass mounted on the pendulum on an axis parallel to the pivot axis and thereby generating a net torque on the pendulum, or by oscillating the pivot point vertically. A simple demonstration of moving the pivot point in a feedback system is achieved by balancing an upturned broomstick on the end of one's finger.

In analytical mechanics, generalized coordinates are a set of parameters used to represent the state of a system in a configuration space. These parameters must uniquely define the configuration of the system relative to a reference state. The generalized velocities are the time derivatives of the generalized coordinates of the system. The adjective "generalized" distinguishes these parameters from the traditional use of the term "coordinate" to refer to Cartesian coordinates.

<span class="mw-page-title-main">Cardioid</span> Type of plane curve

In geometry, a cardioid is a plane curve traced by a point on the perimeter of a circle that is rolling around a fixed circle of the same radius. It can also be defined as an epicycloid having a single cusp. It is also a type of sinusoidal spiral, and an inverse curve of the parabola with the focus as the center of inversion. A cardioid can also be defined as the set of points of reflections of a fixed point on a circle through all tangents to the circle.

In physics, a spherical pendulum is a higher dimensional analogue of the pendulum. It consists of a mass $m$ moving without friction on the surface of a sphere. The only forces acting on the mass are the reaction from the sphere and gravity.

In mathematics, the associated Legendre polynomials are the canonical solutions of the general Legendre equation

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

In mathematics and physics, the Christoffel symbols are an array of numbers describing a metric connection. The metric connection is a specialization of the affine connection to surfaces or other manifolds endowed with a metric, allowing distances to be measured on that surface. In differential geometry, an affine connection can be defined without reference to a metric, and many additional concepts follow: parallel transport, covariant derivatives, geodesics, etc. also do not require the concept of a metric. However, when a metric is available, these concepts can be directly tied to the "shape" of the manifold itself; that shape is determined by how the tangent space is attached to the cotangent space by the metric tensor. Abstractly, one would say that the manifold has an associated (orthonormal) frame bundle, with each "frame" being a possible choice of a coordinate frame. An invariant metric implies that the structure group of the frame bundle is the orthogonal group $O(p, q)$ . As a result, such a manifold is necessarily a (pseudo-)Riemannian manifold. The Christoffel symbols provide a concrete representation of the connection of (pseudo-)Riemannian geometry in terms of coordinates on the manifold. Additional concepts, such as parallel transport, geodesics, etc. can then be expressed in terms of Christoffel symbols.

<span class="mw-page-title-main">Routhian mechanics</span> Formulation of classical mechanics

In classical mechanics, Routh's procedure or Routhian mechanics is a hybrid formulation of Lagrangian mechanics and Hamiltonian mechanics developed by Edward John Routh. Correspondingly, the Routhian is the function which replaces both the Lagrangian and Hamiltonian functions. Although Routhian mechanics is equivalent to Lagrangian mechanics and Hamiltonian mechanics, and introduces no new physics, it offers an alternative way to solve mechanical problems.

A pendulum is a body suspended from a fixed support such that it freely swings back and forth under the influence of gravity. When a pendulum is displaced sideways from its resting, equilibrium position, it is subject to a restoring force due to gravity that will accelerate it back towards the equilibrium position. When released, the restoring force acting on the pendulum's mass causes it to oscillate about the equilibrium position, swinging it back and forth. The mathematics of pendulums are in general quite complicated. Simplifying assumptions can be made, which in the case of a simple pendulum allow the equations of motion to be solved analytically for small-angle oscillations.

The swinging Atwood's machine (SAM) is a mechanism that resembles a simple Atwood's machine except that one of the masses is allowed to swing in a two-dimensional plane, producing a dynamical system that is chaotic for some system parameters and initial conditions.

In physics and mathematics, in the area of dynamical systems, an elastic pendulum is a physical system where a piece of mass is connected to a spring so that the resulting motion contains elements of both a simple pendulum and a one-dimensional spring-mass system. For specific energy values, the system demonstrates all the hallmarks of chaotic behavior and is sensitive to initial conditions.At very low and very high energy, there also appears to be regular motion. The motion of an elastic pendulum is governed by a set of coupled ordinary differential equations.This behavior suggests a complex interplay between energy states and system dynamics.

In celestial mechanics, a Kepler orbit is the motion of one body relative to another, as an ellipse, parabola, or hyperbola, which forms a two-dimensional orbital plane in three-dimensional space. A Kepler orbit can also form a straight line. It considers only the point-like gravitational attraction of two bodies, neglecting perturbations due to gravitational interactions with other objects, atmospheric drag, solar radiation pressure, a non-spherical central body, and so on. It is thus said to be a solution of a special case of the two-body problem, known as the Kepler problem. As a theory in classical mechanics, it also does not take into account the effects of general relativity. Keplerian orbits can be parametrized into six orbital elements in various ways.

<span class="mw-page-title-main">Furuta pendulum</span>

The Furuta pendulum, or rotational inverted pendulum, consists of a driven arm which rotates in the horizontal plane and a pendulum attached to that arm which is free to rotate in the vertical plane. It was invented in 1992 at Tokyo Institute of Technology by Katsuhisa Furuta and his colleagues. It is an example of a complex nonlinear oscillator of interest in control system theory. The pendulum is underactuated and extremely non-linear due to the gravitational forces and the coupling arising from the Coriolis and centripetal forces. Since then, dozens, possibly hundreds of papers and theses have used the system to demonstrate linear and non-linear control laws. The system has also been the subject of two texts.

<span class="mw-page-title-main">Lagrangian mechanics</span> Formulation of classical mechanics

In physics, Lagrangian mechanics is a formulation of classical mechanics founded on the stationary-action principle. It was introduced by the Italian-French mathematician and astronomer Joseph-Louis Lagrange in his presentation to the Turin Academy of Science in 1760 culminating in his 1788 grand opus, Mécanique analytique.

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.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[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/2Map[#.#&,Transpose[v]]-(gy/.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]

[lower-alpha 1]

[2]