Double Pendulum - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


The Double Pendulum

Main Concept

In this Math App we explore the motion of the double pendulum in a constant gravitational field. It is a rather simple physical model, but nevertheless has a rich dynamics as it undergoes chaotic motion.


Fig. 1 Setup of the double pendulum
 
The double pendulum consists of two single pendula, with one being attached at the end of the other. We consider idealized simple pendula, i.e. all their mass is concentrated at the end point of each cord and they move without any friction. The lengths of the pendula are denoted by  and , their masses by  and , respectively. Furthermore, we measure the displacement of each mass via the two angles  and  as indicated in the plot.

The dynamics of the double pendulum can be studied via the Euler-Lagrange equations. We implement and solve these in the next section using Maple procedures. In the subsequent section we discuss some properties of this solution.

Implementation in Maple

Below, we implement the problem of the double pendulum into Maple and solve its dynamics. We start by defining the positions of the masses using the notation in Fig. 1.

(1.1)

The - and -coordinates of the two masses depend on the time,  via the time dependence of the angles. Now we can define the kinetic and potential energy as follows:

(1.2)

(1.3)

Thus the Lagrange function is given via (1.2) and (1.3) as:  

(1.4)

(1.5)

where in the second step we inserted relations (1.1), used trigonometric relations and simplified via combine and simplify. The equations of motion of the two masses can then be obtained via the Euler-Lagrange equations:

for  as:

(1.6)

(1.7)

(1.8)

(1.9)

The resulting equations  and  govern the dynamics of the two masses. They cannot be solved exactly, but we may try solving them numerically using dsolve. To do so, we set the constants to the following values in SI-units ( is given in  , the lengths  in  and the masses:  in )

Furthermore, we set the following initial conditions for the angles and their time derivatives:

with the derivatives given in . Now we can solve the equations of motion (1.7) and (1.9) via

The procedures for the two angles  and  obtained from this solution:

(1.10)

determine the dynamics of the system. In the following section we explore this solution a bit.

Chaotic Motion

The following plot visualizes the procedures obtained above in the time interval .



Fig. 2 Visualization of the solution for  and 

Note that the solution for the angles in this interval evolves chaotically, i.e. it does not exhibit any periodicity. Furthermore, note that the values of  are not constrained to the interval . Whenever the line crosses a gridline, the pendulum flips over.

The dynamics of the double pendulum is chaotic, i.e. it is highly sensitive to the initial conditions . The following plot contains the trajectory and the angles of the mass  for slightly different initial conditions.  

 

Fig. 3 Chaotic motion

 

The three plots show the motion of the double pendulum as parameterized in the previous section. They only differ in the initial condition for the first mass's initial position:

 

Plot 1:

Plot 2:
Plot 3:

 

While under the first condition,  flips over five times, in the second plot it only does so twice and in the third plot four times within the first 25 seconds. Thus, a small change in the initial state can result in large differences in the later behavior of the system.

 

In the following you can explore the motion of the double pendulum. Using the radio buttons, you can choose the strength of the gravitational field. Furthermore, you can adjust the lengths and masses of the two individual pendula using the sliders. The higher the mass, the bigger the circle representing the mass. The initial position of the pendula can be changed by dragging and dropping the masses directly in the plot. Their initial velocities can be adjusted using the sliders; by selecting Show arrows the animation shows the velocities of the masses relative to their pivot points. Furthermore, you may choose the animation duration. You can also plot the animation alongside the trajectory of  by selecting Show trajectory. After you have adjusted the settings, click Animate and wait for the Play button to appear.

Try to answer the following questions:
1) How does the gravitational field influence the motion of the pendulum? What happens as you turn off gravity ()? Is the motion still chaotic?

2) How do the lengths of the individual pendula influence the motion of the double pendulum? What do you observe for  and ?
3) How do the masses influence the motion of the pendulum? What do you observe in the case ? What happens when you change the masses of the pendula in the case  and ?  
4) Investigate the motion of the pendulum for different initial states. Confirm that slight changes in the initial state can lead to very different behaviors at later times.   

 

gravitational field:

 

 

 

 

 

  

 

 

lengths:


masses:              
 
 

initial velocities:  

   
 

 

Drag and drop the two masses into position directly in the plot.

Note: The calculation of the animation delays the appearance

of the Play button. Depending on the chosen animation

duration and whether "Show arrows" or "Show trajectory" is
checked, you may have to wait a few seconds.

animation duration    


Download Help Document