MyPhysicsLab – Chaotic Driven Pendulum

PREV HOME NEXT

The Damped Driven Pendulum

A damped driven pendulum is often used as a basic example of a chaotic system. For a chaotic system the future behavior is highly dependent on the exact value of the initial conditions. A tiny change in initial conditions can cause huge changes after a short period of time.

The pendulum is subject to frictional damping, meaning that it will slow down over time if there is no driving force. To see this try setting the drive amplitude to zero, so the only forces are damping and gravity.

The driving force is shown as a curved arrow in the simulation. The length of the arrow corresponds to the amount of torque (twisting force) at a given moment. The direction of the arrow, clockwise or counter-clockwise, shows the direction of this applied force. You can change the amplitude (strength) of the driving force or its frequency (how often it switches directions).

Not all combinations of the parameters (eg. length, gravity, drive amplitude, drive frequency, damping,...) will lead to chaos. Many combinations result in simpler repeating behavior. Click the "single loop" button above to see this. You'll need to wait about a minute for the simulation to settle into the loop.

There is a phenomenon known as period doubling or bifurcation where slowly increasing one of the parameters (such as drive amplitude) causes the number of loops in the repeating behavior to double. Some of the buttons (above) change the drive amplitude gradually so that you see this period doubling occur.

After enough period doublings occur, the system becomes chaotic. Click one of the buttons labelled "chaos" to set the parameters and see this. A curious thing about this chaotic behavior is that it is not entirely random. While you can't predict the exact state of the system at a given time in the future, it is possible to show that the system will follow an elaborate pattern, These patterns are fractals, which are patterns that repeat themselves when you magnify them.

Physics

The pendulum is modeled as a point mass at the end of a massless rod. The damping (friction) is proportional to the angular velocity of the pendulum. There is also an external driving force which provides a periodic torque (twist). We define the following variables:

• θ = angle of pendulum (0=vertical)
• R = length of rod
• m = mass of pendulum
• g = gravitational constant
• b = damping (friction) constant
• A = amplitude of driving force
• k = constant related to frequency of driving force
• t = time in seconds

We will derive the equation of motion for the pendulum using the rotational analog of Newton's second law for motion about a fixed axis, which is ∑ τ = I α where

• ∑ τ = sum of applied torques
• I = rotational inertia
• ω = angular velocity
• α = θ''= angular acceleration

The rotational inertia about the pivot is I = m R². Torque can be calculated as the vector cross product of the position vector and the force.
The torque due to gravity is τ = -R m g sin θ.
The torque due to friction is τ = -b ω.
The torque due to driving force is τ = A cos(k t).
So we have

m R² α = -R m g sin θ - b ω + A cos(k t)

which simplifies to

       -g           -b θ' + A cos(k t)
θ''  = ——— sin θ + ————————————————————
        R                 m R2

This is the equation of motion for the driven damped pendulum.

Numerical Solution

To solve the equation of motion numerically, so that we can drive the simulation, we use the Runge-Kutta method for solving sets of ordinary differential equations. First we define a variable for the angular velocity ω = θ'. Then we can write the second order equation above as two first order equations.

θ' = ω

      -g           -b ω + A cos(k t)
ω'  = ——— sin θ + ———————————————————
       R                m R2

This is the form needed for using the Runge-Kutta method.