Skip to article frontmatterSkip to article content

Chapter 12 - Simple harmonic motion

12.1Overview

In this chapter, we look at oscillating systems that undergo “simple harmonic motion”, such as the motion of a mass attached to a spring. Many systems in the physical world, such as an oscillating pendulum, can be described by the same mathematical formalism that describes the motion of a mass attached to a spring.

12.2The motion of a spring-mass system

As an example of simple harmonic motion, we first consider the motion of a block of mass mm that can slide without friction along a horizontal surface. The mass is attached to a spring with spring constant kk which is attached to a wall on the other end. We introduce a one-dimensional coordinate system to describe the position of the mass, such that the xx axis is co-linear with the motion, the origin is located where the spring is at rest, and the positive direction corresponds to the spring being extended. This “spring-mass system” is illustrated in Figure 1.

A horizontal spring-mass system oscillating about the origin with an amplitude $A$.

Figure 1:A horizontal spring-mass system oscillating about the origin with an amplitude AA.

We assume that the force exerted by the spring on the mass is given by Hooke’s Law:

F=kxx^\vec F = -kx \hat x

where xx is the position of the mass. The only other forces exerted on the mass are its weight and the normal force from the horizontal surface, which are equal in magnitude and opposite in direction. Therefore, the net force on the mass is the force from the spring.

As we saw in Section Section 8.5, if the spring is compressed (or extended) by a distance AA relative to the rest position, and the mass is then released, the mass will oscillate back and forth between x=±Ax=\pm A[32], which is illustrated in Figure 1. We call AA the “amplitude of the motion”. When the mass is at x=±Ax=\pm A, its speed is zero, as these points correspond to the location where the mass “turns around”.

12.2.1Description using energy

We can describe the motion of the mass using energy, since the mechanical energy of the mass is conserved. At any position, xx, the mechanical energy, EE, of the mass will have a term from the potential energy, UU, associated with the spring force, and kinetic energy, KK:

E=U+K=12kx2+12mv2E = U + K =\frac{1}{2}kx^2 + \frac{1}{2}mv^2

We can find the mechanical energy, EE, by evaluating the energy at one of the turning points. At these points, the kinetic energy of the mass is zero, so E=U(x=A)=1/2kA2E=U(x=A)=1/2kA^2. We can then write the expression for mechanical energy as:

12kx2+12mv2=12kA2\boxed{\frac{1}{2}kx^2 + \frac{1}{2}mv^2 = \frac{1}{2}kA^2}

We can thus always know the speed, vv, of the mass at any position, xx, if we know the amplitude AA:

v(x)=k(A2x2)mv(x) = \sqrt{\frac{k(A^2-x^2)}{m}}

12.2.2Kinematics of simple harmonic motion

We can use Newton’s Second Law to obtain the position, x(t)x(t), velocity, v(t)v(t), and acceleration, a(t)a(t), of the mass as a function of time. The xx component of Newton’s Second Law for the mass attached to the spring can be written:

Fx=kx=ma\sum F_x = -kx = ma

We can write the acceleration in Newton’s Second Law more explicitly as the second derivative of the position, x(t)x(t), with respect to time. If we do this, we can see that Newton’s Second Law for the mass attached to the spring is a differential equation for the function x(t)x(t) (we call it an “equation of motion”):

ma=kxmd2xdt2=kxd2xdt2=kmx\begin{align*} ma &= -kx\nonumber\\ m\frac{d^2x}{dt^2} &= -kx\nonumber\\ \therefore \frac{d^2x}{dt^2} &= -\frac{k}{m}x \end{align*}

We want to find the position function, x(t)x(t). Equation {\ref{eq:simpleharmonicmotion:shmspring}} tells us that the second derivative of x(t)x(t) with respect to time must equal the negative of the x(t)x(t) function multiplied by a constant, k/mk/m. Without having taken a course on differential equations, it might not be obvious what the function x(t)x(t) could be. Several, equivalent functions can satisfy this equation. One possible choice, which we present here as a guess, is\footnote{Other possible guesses that work are Asin(ωt+ϕ)A \sin(\omega t + \phi), and x(t)=Acos(ωt)+Bsin(ωt)x(t) = A\cos(\omega t) + B\sin(\omega t).}:

x(t)=Acos(ωt+ϕ)\boxed{x(t) = A \cos(\omega t + \phi)}

where AA, ω, and ϕ are constants that we need to determine. We can take the second order derivative with respect to time of the function above to verify that it indeed “solves” the differential equation:

x(t)=Acos(ωt+ϕ)ddtx(t)=Aωsin(ωt+ϕ)d2dt2x(t)=ddt(Aωsin(ωt+ϕ))=Aω2cos(ωt+ϕ)d2dt2x(t)=ω2x(t)\begin{align*} x(t) &= A \cos(\omega t + \phi)\\ \frac{d}{dt}x(t) &= -A\omega\sin(\omega t + \phi)\\ \frac{d^2}{dt^2}x(t) &=\frac{d}{dt}\left( -A\omega\sin(\omega t + \phi)\right)= -A\omega^2\cos(\omega t + \phi)\\ \therefore \frac{d^2}{dt^2}x(t) &= - \omega^2 x(t) \end{align*}

The last equation has exactly the same form as (6), which we obtained from Newton’s Second Law, if we define ω as:

ω=km\boxed{\omega = \sqrt{\frac{k}{m}}}

We call ω the “angular frequency” of the spring-mass system. We have found that our guess for x(t)x(t) satisfies the differential equation.

We still need to identify what the constants AA and ϕ have to do with the motion of the mass. The constant AA is the maximal value that x(t)x(t) can take (when the cosine is equal to 1). This corresponds to the amplitude of the motion of the mass, which we already had labelled, AA. The constant, ϕ, is called the “phase” and depends on when we choose t=0t=0 to be. Suppose that we define time t=0t=0 to be when the mass is at x=Ax=A; in that case:

x(t=0)=AAcos(ωt+ϕ)=AAcos(ω(0)+ϕ)=Acos(ϕ)=1ϕ=0\begin{align*} x(t=0) &= A\\ A \cos(\omega t + \phi) &= A\\ A \cos(\omega (0) + \phi) &= A\\ \cos(\phi) &= 1\\ \therefore \phi = 0 \end{align*}

If we define t=0t=0 to be when the mass is at x=Ax=A, then the phase, ϕ, is zero. In general, the value of ϕ can take any value between π-\pi and +π+\pi[33] and, physically, corresponds to our choice of when t=0t=0 (i.e. the position of the mass when we choose t=0t=0).

Since we have determined the position as a function of time for the mass, its velocity and acceleration as a function of time are easily found by taking the corresponding time derivatives:

x(t)=Acos(ωt+ϕ)v(t)=ddtx(t)=Aωsin(ωt+ϕ)a(t)=ddtv(t)=Aω2cos(ωt+ϕ)\begin{align*} x(t) &= A \cos(\omega t + \phi)\\ v(t) &= \frac{d}{dt}x(t) = -A\omega\sin(\omega t + \phi)\\ a(t)&= \frac{d}{dt}v(t) = -A\omega^2\cos(\omega t + \phi) \end{align*}

The position of the mass is described by a sinusoidal function of time; we call this type of motion “simple harmonic motion”. The position and velocity as a function of time for a spring-mass system with m=1kgm=1 {\rm kg}, k=4N/mk=4 {\rm N/m}, A=10mA=10 {\rm m} are shown in Figure 2 for two different choices of the phase, ϕ=0\phi=0 and ϕ=π/2\phi=\pi/2.

Position and velocity as a function of time for a mass-spring system for two different values of the phase, $\phi$.

Figure 2:Position and velocity as a function of time for a mass-spring system for two different values of the phase, ϕ.

We can make a few observations about the position and velocity illustrated in Figure 2:

  • Changing the phase, ϕ, results in an horizontal shift of the functions. A positive phase results in a shift of the functions to the left.
  • The highest speed corresponds to a position of x=0x=0 and the largest position, x=±Ax=\pm A, corresponds to a speed of zero.
  • ϕ=0\phi = 0 corresponds to the “initial condition” at t=0t=0, where the position of the mass is x=Ax=A and its speed is v=0v=0.
  • ϕ=π/2\phi = \pi/2 corresponds to the “initial condition” at t=0t=0, where the position of the mass is x=0x=0 and its velocity is in the negative direction, and with maximal amplitude.
  • The position is always between x=±Ax=\pm A, and the velocity is always between v=±Aωv=\pm A\omega.

The motion of the spring is clearly periodic. If the period of the motion is TT, then the position of the mass at time tt will be the same as its position at t+Tt+T. The period of the motion, TT, is easily found:

T=2πω=2πmk\boxed{T=\frac{2\pi}{\omega}=2\pi\sqrt{\frac{m}{k}}}

And the corresponding frequency is given by:

f=1T=ω2π=12πkm\boxed{f = \frac{1}{T}=\frac{\omega}{2\pi}=\frac{1}{2\pi}\sqrt{\frac{k}{m}}}

It should now be clear why ω is called the angular frequency, since it is related to the frequency of the motion.

12.2.3Analogy with uniform circular motion

We can make an analogy between the mathematical description of the motion of a spring-mass system and that of uniform circular motion. Consider a particle that is moving along a circle of radius AA, with constant angular speed ω, as illustrated in Figure 3.

Uniform circular motion of a particle along a circle of radius $A$ with constant angular speed $\omega$.

Figure 3:Uniform circular motion of a particle along a circle of radius AA with constant angular speed ω.

The angular position, θ(t)\theta(t), of the particle is given by:

θ(t)=θ0+ωt\theta(t) = \theta_0 + \omega t

if the particle was located at an angular position θ0\theta_0 at t=0t=0 (θ0=0\theta_0=0 in Figure 3). The xx coordinate of the particle is given by:

x(t)=Acos(θ(t))=Acos(θ0+ωt)x(t) = A\cos(\theta(t)) = A\cos(\theta_0 + \omega t)

We can see that the xx coordinate of the particle has the same functional form as the position for simple harmonic motion. The same is true for the particle’s velocity. The magnitude of the particle’s velocity is given by:

v=ωr=ωAv = \omega r = \omega A

where r=Ar=A is the radius of the circle. The xx component of the particle’s velocity is easily found from the figure and is given by:

vx(t)=vsin(θ(t))=ωAsin(θ0+ωt)v_x(t) = -v\sin(\theta(t)) = -\omega A\sin(\theta_0 + \omega t)

We can visualize simple harmonic motion as if it were the projection onto the xx axis of uniform circular motion with angular speed ω about a circle with radius AA. The phase ϕ corresponds to the angular position of the particle around the circle, θ0\theta_0, at time t=0t=0. When the particle crosses the yy axis (x=0x=0), its velocity is in the xx direction, so the xx component of the velocity is maximal. When the particle crosses the xx axis (x=±Ax=\pm A), the xx component of the velocity is zero.

12.3Vertical spring-mass system

Consider the vertical spring-mass system illustrated in Figure 5.

A vertical spring-mass system.

Figure 5:A vertical spring-mass system.

When no mass is attached to the spring, the spring is at rest (we assume that the spring has no mass). We choose the origin of a one-dimensional vertical coordinate system (yy axis) to be located at the rest length of the spring (left panel of Figure 5). When a mass mm is attached to the spring, the spring will extend and the end of the spring will move to a new equilibrium position, y0y_0, given by the condition that the net force on the mass mm is zero. The only forces exerted on the mass are the force from the spring and its weight. The condition for the equilibrium is thus:

Fy=FgF(y0)=0mgky0=0mg=ky0\begin{align*} \sum F_y = F_g - F(y_0) &=0\\ mg - ky_0 &= 0 \\ \therefore mg &= ky_0 \end{align*}

Now, consider the forces on the mass at some position yy when the spring is extended downwards relative to the equilibrium position (right panel of Figure 5). Newton’s Second Law at that position can be written as:

Fy=mgky=mamd2ydt2=mgky\begin{align*} \sum F_y = mg - ky &= ma\\ \therefore m \frac{d^2y}{dt^2}& = mg - ky \end{align*}

Note that the net force on the mass will always be in the direction so as to “restore” the position of the mass back to the equilibrium position, y0y_0. If the mass had been moved upwards relative to y0y_0, the net force would be downwards.

We can substitute the equilibrium condition, mg=ky0mg = ky_0, into the equation that we obtained from Newton’s Second Law:

md2ydt2=mgkymd2ydt2=ky0kymd2ydt2=k(yy0)d2ydt2=km(yy0)\begin{align*} m \frac{d^2y}{dt^2}& = mg - ky \\ m \frac{d^2y}{dt^2}&= ky_0 - ky\\ m \frac{d^2y}{dt^2}&=-k(y-y_0) \\ \therefore \frac{d^2y}{dt^2} &= -\frac{k}{m}(y-y_0) \end{align*}

Consider a new variable, y=yy0y'=y-y_0. This is the same as defining a new yy' axis that is shifted downwards by y0y_0; in other words, this the same as defining a new yy' axis whose origin is at y0y_0 (the equilibrium position) rather than at the position where the spring is at rest. Noting that the second time derivative of y(t)y'(t) is the same as that for y(t)y(t):

d2ydt2=d2dt2(y+y0)=d2ydt2\begin{align*} \frac{d^2y}{dt^2} &= \frac{d^2}{dt^2} (y' + y_0) = \frac{d^2y'}{dt^2}\\ \end{align*}

we can write the equation of motion for the mass, but using y(t)y'(t) to describe its position:

d2ydt2=kmy\begin{align*} \frac{d^2y'}{dt^2} &= -\frac{k}{m}y' \end{align*}

This is the same equation as that for the simple harmonic motion of a horizontal spring-mass system (6), but with the origin located at the equilibrium position instead of at the rest length of the spring. In other words, a vertical spring-mass system will undergo simple harmonic motion in the vertical direction about the equilibrium position. In general, a spring-mass system will undergo simple harmonic motion if a constant force that is co-linear with the spring force is exerted on the mass (in this case, gravity). That motion will be centred about a point of equilibrium where the net force on the mass is zero rather than where the spring is at its rest position.

12.3.1Two-spring-mass system

Consider a horizontal spring-mass system composed of a single mass, mm, attached to two different springs with spring constants k1k_1 and k2k_2, as shown in Figure 6.

A mass attached to two different springs.

Figure 6:A mass attached to two different springs.

We introduce a horizontal coordinate system, such that the end of the spring with spring constant k1k_1 is at position x1x_1 when it is at rest, and the end of the k2k_2 spring is at x2x_2 when it is as rest, as shown in the top panel. A mass mm is then attached to the two springs, and x0x_0 corresponds to the equilibrium position of the mass when the net force from the two springs is zero. We will assume that the length of the mass is negligible, so that the ends of both springs are also at position x0x_0 at equilibrium. You can see in the middle panel of Figure 6 that both springs are in extension when in the equilibrium position. It is possible to have an equilibrium where both springs are in compression, if both springs are long enough to extend past x0x_0 when they are at rest.

If we assume that both springs are in extension at equilibrium, as shown in the figure, then the condition for equilibrium is given by requiring that the sum of the forces on the mass is zero when the mass is located at x0x_0. The extension of the spring on the left is x0x1x_0 - x_1, and the extension of the spring on the right is x2x0x_2-x_0:

Fx=k1(x0x1)+k2(x2x0)=0k1x0+k1x1+k2x2k2x0=0(k1+k2)x0+k1x1+k2x2=0k1x1+k2x2=(k1+k2)x0\begin{align*} \sum F_x = -k_1(x_0-x_1) + k_2 (x_2 - x_0) &= 0\\ -k_1x_0+k_1x_1+k_2x_2-k_2x_0 &=0\\ -(k_1+k_2)x_0 +k_1x_1+k_2x_2 &=0\\ \therefore k_1x_1+k_2x_2 &=(k_1+k_2)x_0 \end{align*}

Note that if the mass is displaced from x0x_0 in any direction, the net force on the mass will be in the direction of the equilibrium position, and will act to “restore” the position of the mass back to x0x_0.

When the mass is at some position xx, as shown in the bottom panel (for the k1k_1 spring in compression and the k2k_2 spring in extension), Newton’s Second Law for the mass is:

k1(xx1)+k2(x2x)=mak1x+k1x1+k2x2k2x=md2xdt2(k1+k2)x+k1x1+k2x2=md2xdt2\begin{align*} -k_1(x-x_1) + k_2 (x_2 - x) &= m a \\ -k_1x +k_1x_1 + k_2 x_2 - k_2 x &= m \frac{d^2x}{dt^2}\\ -(k_1+k_2)x + k_1x_1 + k_2 x_2&= m \frac{d^2x}{dt^2} \end{align*}

Note that, mathematically, this equation is of the form kx+C=ma-kx + C =ma, which is the same form of the equation that we had for the vertical spring-mass system (with C=mgC=mg), so we expect that this will also lead to simple harmonic motion. We can use the equilibrium condition (k1x1+k2x2=(k1+k2)x0k_1x_1+k_2x_2 =(k_1+k_2)x_0) to re-write this equation:

(k1+k2)x+k1x1+k2x2=md2xdt2(k1+k2)x+(k1+k2)x0=md2xdt2(k1+k2)(xx0)=md2xdt2\begin{align*} -(k_1+k_2)x + k_1x_1 + k_2 x_2&= m \frac{d^2x}{dt^2}\\ -(k_1+k_2)x + (k_1+k_2)x_0&= m \frac{d^2x}{dt^2}\\ \therefore -(k_1+k_2) (x-x_0) &= m \frac{d^2x}{dt^2} \end{align*}

Let us define k=k1+k2k=k_1+k_2 as the “effective” spring constant from the two springs combined. We can also define a new coordinate, x=xx0x' = x-x_0, which simply corresponds to a new xx axis whose origin is located at the equilibrium position (in a way that is exactly analogous to what we did in the vertical spring-mass system). We can thus write Newton’s Second Law as:

(k1+k2)(xx0)=md2xdt2kx=md2xdt2d2xdt2=kmx\begin{align*} -(k_1+k_2) (x-x_0) &= m \frac{d^2x}{dt^2}\\ -kx' &= m \frac{d^2x'}{dt^2}\\ \therefore \frac{d^2x'}{dt^2} &= -\frac{k}{m}x' \end{align*}

and we find that the motion of the mass attached to two springs is described by the same equation of motion for simple harmonic motion as that of a mass attached to a single spring. In this case, the mass will oscillate about the equilibrium position, x0x_0, with a an effective spring constant k=k1+k2k=k_1+k_2. Combining the two springs in this way is thus equivalent to having a single spring, but with spring constant k=k1+k2k=k_1+k_2. The angular frequency of the oscillations is given by:

ω=km=k1+k2m\omega = \sqrt{\frac{k}{m}}=\sqrt{\frac{k_1+k_2}{m}}

12.4Simple harmonic motion

In the previous sections, we modelled the motion of a mass attached to a spring and found that its position, x(t)x(t), was described by the following differential equation:

d2xdt=ω2x\boxed{\frac{d^2x}{dt} = -\omega^2x}

A possible solution to that equation was given by:

x(t)=Acos(ωt+ϕ)\boxed{x(t) = A\cos(\omega t+ \phi)}

We then saw that the motion of a vertical spring-mass system, as well as that of a mass attached to two springs, could also be described by (28). Any physical system that can described by (28) is said to undergo “simple harmonic motion”, or to be a “simple harmonic oscillator”. If we find that the physical model of a system leads to (28), then we immediately know that the position of system can be described by (29).

The key physical characteristic of a simple harmonic oscillator is that there is a “restoring force” whose magnitude is proportional to the displacement from the equilibrium position. A restoring force is a force that acts to place the system back in equilibrium, and is thus always in the direction that is opposite of the displacement relative to an equilibrium position. In the three systems that we considered so far, the net force on the mass was always such that it would restore the mass back to the equilibrium position, where the net force on the mass is zero.

Many systems in nature are well modelled as simple harmonic oscillators. Some examples are: the motion of a pendulum as it oscillates, the motion of a buoy bobbing up and down in the sea, the motion of electrons in a shorted capacitor, and the vibrations of atoms in a molecule.

12.5The motion of a pendulum

In this section, we show how and when the motion of a pendulum can be described as simple harmonic motion. Consider the simple pendulum that is constructed from a mass-less string of length, LL, attached to a fixed point on one end and to a point mass mm on the other, as illustrated in Figure 7.

A simple pendulum which oscillates in a vertical plane.

Figure 7:A simple pendulum which oscillates in a vertical plane.

The pendulum can swing in the vertical plane, and we have shown our choice of coordinate system (the zz axis, not shown, is out of the page). The only two forces on the mass are the tension from the string and its weight. We can describe the position of the mass by the angle, θ(t)\theta(t), that the string makes with the vertical. We can model the dynamics of the simple pendulum by considering the net torque and angular acceleration about the axis of rotation that is perpendicular to the plane of the page and that goes through the point on the string that is fixed.

The force of tension cannot create a torque on the mass about the axis of rotation, as it is anti-parallel to the vector from the point of rotation to the mass. The net torque is thus the torque from the force of gravity:

τnet=τg=r×Fg=(Lsinθx^Lcosθy^)×(mgy^)=mgLsinθz^\begin{align*} \vec\tau^{net} &=\vec \tau_g \\ &=\vec r \times \vec F_g = (L\sin\theta \hat x - L\cos\theta \hat y) \times (-mg\hat y)\\ &=-mgL\sin\theta \hat z \end{align*}

where LL is the magnitude of the vector, r\vec r, from the axis of rotation to where the force of gravity is exerted. The net torque is equal to the angular acceleration, α, multiplied by the moment of inertia, II, of the mass:

τnet=IαmgLsinθz^=mL2αgsinθz^=Lα\begin{align*} \vec\tau^{net} &= I\vec\alpha\\ -mgL\sin\theta \hat z&= mL^2 \vec\alpha\\ -g\sin\theta \hat z&= L \vec\alpha \end{align*}

where I=ML2I=ML^2 is the moment of inertia for a point mass a distance LL away from the axis of rotation. For the position illustrated in Figure 7, the angular acceleration of the pendulum is in the negative zz direction (into the page) and corresponds to a clockwise motion for the pendulum, as we would expect. The angular acceleration is the second time derivative of the angle, θ:

α=d2θdt2\alpha = \frac{d^2\theta}{dt^2}

We can thus re-write the equation that we obtained from the rotational dynamics version of Newton’s Second Law as:

gsinθz^=Lαd2θdt2=gLsinθ\begin{align*} -g\sin\theta \hat z&= L \vec\alpha\\ \frac{d^2\theta}{dt^2} &= -\frac{g}{L}\sin\theta \end{align*}

where we only used the magnitudes in the second equation, since all of the angular quantities are in the zz direction. This equation of motion for θ(t)\theta(t) almost looks like the equation for simple harmonic oscillation for the angle θ (except that we have sinθ\sin\theta instead of θ). However, consider the “the small angle approximation”\footnote{Look up the Maclaurin/Taylor series for the sine function!} for the sine function:

sinθθ\sin\theta \approx \theta

If the oscillations of the pendulum are “small”, such that the small angle approximation is valid, then the equation of motion for the pendulum is:

d2θdt2=gLsinθgLθd2θdt2=gLθ(for small θ)\begin{align*} \frac{d^2\theta}{dt^2} &= -\frac{g}{L}\sin\theta \approx -\frac{g}{L}\theta \\ \therefore \frac{d^2\theta}{dt^2} &=-\frac{g}{L}\theta \quad (\text{for small }\theta) \end{align*}

and the angle that the pendulum makes with the vertical is described by the equation for simple harmonic oscillation with angular frequency:

ω=gL\omega = \sqrt{\frac{g}{L}}

The angle, θ, as a function of time is thus described by the function:

θ(t)=θmaxcos(ωt+ϕ)\theta(t) = \theta_{max}\cos(\omega t +\phi)

where θmax\theta_{max} is the maximal amplitude of the oscillations and ϕ is a phase that depends on when we choose to define t=0t=0.

12.5.1The physical pendulum

A physical pendulum is defined as any object that is allowed to rotate in the vertical plane about some axis that goes through the object, as illustrated in Figure 8.

A physical pendulum which oscillates in a vertical plane about an axis through the object.

Figure 8:A physical pendulum which oscillates in a vertical plane about an axis through the object.

The only forces exerted on the pendulum are its weight (exerted at its centre of mass) and a contact force exerted at the axis of rotation. The physical pendulum can be modelled in exactly the same way as the simple pendulum, except that we use the moment of inertia of the object about the axis of rotation. Only the weight results in a torque about the rotation axis, since the contact force is exerted at the rotation axis:

τnet=τg=Iαmghsinθ=Iα=Id2θdt2\begin{align*} \tau^{net} = \tau_g &= I\alpha\\ -mgh\sin\theta &= I\alpha = I \frac{d^2\theta}{dt^2} \end{align*}

where hh is the distance from the axis of rotation to the centre of mass. In the small angle approximation, this becomes:

d2θdt2=mghIθ(for small θ)\begin{align*} \frac{d^2\theta}{dt^2} &=-\frac{mgh}{I}\theta \quad (\text{for small }\theta) \end{align*}

and we find that the physical pendulum oscillates with an angular frequency:

ω=mghI\omega = \sqrt{\frac{mgh}{I}}

12.6Summary

The equation of motion for the position, x(t)x(t), of the mass in a one-dimensional spring-mass system with no friction can be written:

d2xdt2=kmx=ω2x\frac{d^2x}{dt^2}=-\frac{k}{m}x = -\omega^2 x

and has a solution:

x(t)=Acos(ωt+ϕ)x(t) = A\cos(\omega t + \phi)

where AA is the amplitude of the motion, ϕ is the phase, which depends on our choice of initial conditions (when we choose time t=0t=0), and ω:

ω=km\omega = \sqrt{\frac{k}{m}}

is the angular frequency of the motion. The mass will oscillate about an equilibrium position with a period, TT, and frequency, ff, given by:

T=2πω=2πmkf=1T=ω2π=12πkm\begin{align*} T&=\frac{2\pi}{\omega}=2\pi\sqrt{\frac{m}{k}}\\ f&=\frac{1}{T}=\frac{\omega}{2\pi}=\frac{1}{2\pi}\sqrt{\frac{k}{m}} \end{align*}

The velocity and acceleration of the mass are found by taking the time derivatives of the position x(t)x(t):

x(t)=Acos(ωt+ϕ)v(t)=ddtx(t)=Aωsin(ωt+ϕ)a(t)=d2dt2x(t)=ddt(Aωsin(ωt+ϕ))=Aω2cos(ωt+ϕ)\begin{align*} x(t)&= A \cos(\omega t + \phi)\\ v(t)&=\frac{d}{dt}x(t) = -A\omega\sin(\omega t + \phi)\\ a(t)&= \frac{d^2}{dt^2}x(t) =\frac{d}{dt}\left( -A\omega\sin(\omega t + \phi)\right)= -A\omega^2\cos(\omega t + \phi) \end{align*}

The total mechanical energy of the mass, at some position xx, is given by:

E=U+K=12kx2+12mv2=12kA2E =U+K=\frac{1}{2}kx^2+\frac{1}{2}mv^2= \frac{1}{2}kA^2

and is conserved.

Any system that can be described by the equation of motion:

d2xdt2=ω2x\frac{d^2x}{dt^2}= -\omega^2 x

is said to be a simple harmonic oscillator, and its position will be described by:

x(t)=Acos(ωt+ϕ)x(t) = A\cos(\omega t + \phi)

A simple harmonic oscillator will always oscillate about an equilibrium position, where the net force on the oscillator is zero. The net force on a simple harmonic oscillator is always directed towards the equilibrium position, and has a magnitude proportional to the distance of the oscillator from its equilibrium position. The force is called a restoring force. A vertical spring-mass system, and a mass attached to two springs will both undergo simple harmonic motion about their respective equilibrium position.

A simple pendulum will undergo simple harmonic oscillations, if the amplitude of the oscillations is small. The angular frequency for the oscillations of a simple pendulum only depends on the length of the pendulum:

ω=gL\omega = \sqrt{\frac{g}{L}}

This is valid in the small angle approximation, where:

sinθθ\sin\theta \approx \theta

A physical pendulum of mass mm which oscillates about an axis through the object will also undergo simple harmonic oscillation in the small angle approximation. The angular frequency of the oscillations for a physical pendulum is given by:

ω=mghI\omega = \sqrt{\frac{mgh}{I}}

where hh is the distance between the centre of mass and the axis of rotation, and II is the moment of inertia of the object about the rotation axis.

12.7Thinking about the material

12.8Sample problems and solutions

12.8.1Problems

12.8.2Solutions