myphysicslab – simple pendulum
this simulation shows a simple pendulum operating under gravity. for small oscillations the pendulum is linear, but it is non-linear for larger oscillations.
you can change parameters in the simulation such as mass, gravity, and friction (damping). you can drag the pendulum with your mouse to change the starting position. if you don't see the simulation try instructions for enabling java. scroll down to see the math!
try using the graph and changing parameters like mass, length, gravity to answer these questions about the spring simulation:
- what is the relationship between angular acceleration and angle?
- how do mass, length, or gravity affect the relationship between angular acceleration and angle?
- for small oscillations, how do length or gravity affect the period or frequency of the oscillation?
note: leave damping and drive frequency set to zero here (they complicate things). you'll find the answers
below. (hint: try starting the pendulum from an almost vertically "up" position.)
physics - rotational method
the pendulum is modeled as a point mass at the end of a massless rod. we define the following variables:
- θ = angle of pendulum (0=vertical)
- r = length of rod
- t = tension in rod
- m = mass of pendulum
- g = gravitational constant
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 α
- τ = net torque
- i = rotational inertia
- α = θ''= angular acceleration
the rotational inertia about the pivot is i = m r2
. torque can be calculated as the vector cross product of the position vector and the force. the magnitude of the torque due to gravity works out to be τ = −r m g sin θ
. so we have
−r m g sin θ = m r2 α
which simplifies to
this is the equation of motion for the pendulum.
physics - direct method
most students are less familiar with rotational inertia and torque than with the simple mass and acceleration found in newton's second law, f = m a
. to show that there is nothing new in the rotational version of newton's second law, we derive the equation of motion here without the rotational dynamics. as you will see, this method involves more algebra.
we'll need the standard unit vectors, i, j
. we use bold and overline to indicate a vector.
- i = unit vector in horizontal direction
- j = unit vector in vertical direction
the kinematics of the pendulum are then as follows
position = r sin θ i − r cos θ j
velocity = r θ' cos θ i + r θ' sin θ j
acceleration = r(θ'' cos θ i −
θ' 2 sin θ i +
θ'' sin θ j +
θ' 2 cos θ j)
the position is derived by a fairly simple application of trigonometry. the velocity and acceleration are then the first and second derivatives of the position.
next we draw the free body diagram for the pendulum. the forces on the pendulum are the tension in the rod t
and gravity. so we can write the net force as:
f = t cos θ j − t sin θ i − m g j
using newton's law f = m a
and the pendulum acceleration we found earlier, we have
t cos θ j − t sin θ i − m g j =
m r(θ'' cos θ i −
θ' 2 sin θ i +
θ'' sin θ j +
θ' 2 cos θ j)
we can write the vector components of the above equation as separate equations. this gives us two simultaneous equations: the first for the i
component and the second for the j
−t sin θ = m r(θ'' cos θ − θ' 2 sin θ)
t cos θ − m g = m r(θ'' sin θ + θ' 2 cos θ)
next we do some algebraic manipulations to eliminate the unknown t
. multiply the first equation by cos θ
and the second by sin θ
−t sin θ cos θ = m r(θ'' cos2θ − θ' 2 sin θ cos θ)
t cos θ sin θ − m g sin θ =
m r(θ'' sin2θ +
θ' 2 sin θ cos θ)
use the first equation to substitute for t cos θ sin θ
in the second equation and do a little more algebra to get:
−θ'' cos2θ + θ' 2 sin θ cos θ = θ'' sin2θ +
θ' 2 sin θ cos θ +
g⁄r sin θ
with the trig identity cos2θ + sin2θ = 1
this simplifies to equation (1)
θ'' = − g⁄r sin θ
physics - energy method
there is yet a third way to derive the equations of motion for the pendulum. this is to use the "indirect" energy based method associated with the terms "lagrangian", "euler-lagrange equations", "hamiltonian", and others. while this method isn't shown here, you can see an example of it for the pendulum+cart simulation
to solve the equations 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 (1) as two first order equations.
θ' = ω
ω' = − g⁄r sin θ
this is the form needed for using the runge-kutta method.
question: what is the relationship between angular acceleration and angle?
answer: it is a sine wave relationship as given by equation (1):
θ'' = − g⁄r sin θ
question: how do mass, length, or gravity affect the relationship between angular acceleration and angle?
answer: from equation (1) we see that:
- mass doesn't affect the motion at all.
- the amplitude of the sine relationship is proportional to gravity.
- the amplitude of the sine relationship is inversely proportional to length of the pendulum.
question: for small oscillations, how do length or gravity affect the period or frequency of the oscillation?
answer: for small oscillations we can use the approximation that sin θ = θ. then the equation of motion becomes
θ'' = − g⁄r θ
this is a linear relationship. you can see that the graph of acceleration versus angle is a straight line for small oscillations. this is the same form of equation as for the single spring simulation. the analytic solution is
where θ0 is the initial angle and t is time. the period is the time it takes for θ(t) to repeat, so
the frequency of oscillation is the inverse of the period:
so we predict that
- increasing length by 4 times doubles the period and halves the frequency;
- increasing gravity by 4 times halves the period and doubles the frequency;