EduKG
← Back

harmonic oscillator kinematics and dynamics

Learning objectives

  • Outcome

    Derive the equation of motion a=−ω2xa = -\omega^2 x from F=−kxF = -kx and Newton's second law, and state ω=k/m\omega = \sqrt{k/m}.

  • Outcome

    Calculate xx, vv and aa at a given time from x(t)=Acos⁡(ωt+φ)x(t) = A\cos(\omega t + \varphi).

  • Outcome

    Apply conservation of mechanical energy 12kx2+12mv2=12kA2\tfrac12kx^2 + \tfrac12mv^2 = \tfrac12kA^2 to find the speed at a given displacement and the maximum speed.

  • Outcome

    Sketch and interpret graphs of x(t)x(t), v(t)v(t), a(t)a(t) over one period, identifying amplitude, period and the phase shifts between them.

Understand the idea

Kinematics and Dynamics of a Harmonic Oscillator

Text

A harmonic oscillator is any system whose restoring force is proportional to the displacement from equilibrium and directed back toward it: F=−kxF = -kx.

Newton's second law then gives ma=−kxma = -kx, so the acceleration obeys a=−ω2xa = -\omega^2 x, where ω=km,ω=2πT;\omega = \sqrt{\frac{k}{m}}, \qquad \omega = \frac{2\pi}{T}; this relation is the dynamics of the oscillator. Solving that equation gives the kinematics of the motion: x(t)=Acos⁡(ωt+φ).x(t) = A\cos(\omega t + \varphi). Differentiating twice yields v(t)=−Aωsin⁡(ωt+φ),a(t)=−Aω2cos⁡(ωt+φ)=−ω2x(t).v(t) = -A\omega\sin(\omega t + \varphi), \qquad a(t) = -A\omega^2\cos(\omega t + \varphi) = -\omega^2 x(t). The amplitude AA and the initial phase φ\varphi are fixed by the initial position and velocity, which relies on the prerequisites concerning potential energy and conservation of mechanical energy. Because the force is conservative with U=12kx2U = \tfrac12 kx^2, the total mechanical energy E=12kx2+12mv2=12kA2=constE = \tfrac12kx^2 + \tfrac12mv^2 = \tfrac12kA^2 = \text{const} is conserved, which lets you find speeds without solving the motion in time.

Worked examples

Direct

Problem

A mass on a spring oscillates as x(t)=Acos⁡(ωt)x(t) = A\cos(\omega t) with A=0.10 mA = 0.10\ \text{m} and period T=2.0 sT = 2.0\ \text{s}. Find xx, vv and aa at t=0.50 st = 0.50\ \text{s}.

Solution

Angular frequency: ω=2π/T=2π/2.0=π≈3.14 rad/s\omega = 2\pi/T = 2\pi/2.0 = \pi \approx 3.14\ \text{rad/s}. Hence ωt=π⋅0.50=π/2≈1.57 rad\omega t = \pi \cdot 0.50 = \pi/2 \approx 1.57\ \text{rad}.

Displacement: x=0.10cos⁡(π/2)=0 mx = 0.10\cos(\pi/2) = 0\ \text{m} — the body is at the equilibrium position.

Velocity: v=−Aωsin⁡(ωt)=−0.10⋅3.14⋅1=−0.31 m/sv = -A\omega\sin(\omega t) = -0.10 \cdot 3.14 \cdot 1 = -0.31\ \text{m/s} (moving toward negative xx).

Acceleration: a=−ω2x=0 m/s2a = -\omega^2 x = 0\ \text{m/s}^2, which agrees with a=−Aω2cos⁡(π/2)=0a = -A\omega^2\cos(\pi/2) = 0. At equilibrium the speed is maximum and the acceleration is zero.

Application

Problem

A block of mass m=0.50 kgm = 0.50\ \text{kg} is attached to a spring of stiffness k=8.0 N/mk = 8.0\ \text{N/m} and released from rest at x0=0.040 mx_0 = 0.040\ \text{m}. Find ω\omega, the period TT, the maximum speed, and the speed at x=0.020 mx = 0.020\ \text{m}.

Solution

Frequency: ω=k/m=8.0/0.50=16=4.0 rad/s\omega = \sqrt{k/m} = \sqrt{8.0/0.50} = \sqrt{16} = 4.0\ \text{rad/s}.

Period: T=2π/ω=2π/4.0=1.57 sT = 2\pi/\omega = 2\pi/4.0 = 1.57\ \text{s}.

Releasing from rest means starting at an endpoint, so the amplitude is A=x0=0.040 mA = x_0 = 0.040\ \text{m}.

Maximum speed (at x=0x = 0): vmax⁡=Aω=0.040⋅4.0=0.16 m/sv_{\max} = A\omega = 0.040 \cdot 4.0 = 0.16\ \text{m/s}.

Speed at x=0.020 mx = 0.020\ \text{m} follows from conservation of energy:

12kx2+12mv2=12kA2  ⇒  v=ωA2−x2.\tfrac12kx^2 + \tfrac12mv^2 = \tfrac12kA^2 \;\Rightarrow\; v = \omega\sqrt{A^2 - x^2}.

Substituting: v=4.00.0402−0.0202=4.01.2×10−3≈0.14 m/sv = 4.0\sqrt{0.040^2 - 0.020^2} = 4.0\sqrt{1.2\times10^{-3}} \approx 0.14\ \text{m/s}.

Energy check: E=12kA2=6.4×10−3 JE = \tfrac12kA^2 = 6.4\times10^{-3}\ \text{J}, U=12kx2=1.6×10−3 JU = \tfrac12kx^2 = 1.6\times10^{-3}\ \text{J}, so K=4.8×10−3 JK = 4.8\times10^{-3}\ \text{J} and v=2K/m≈0.14 m/sv = \sqrt{2K/m} \approx 0.14\ \text{m/s}.

Transfer

Problem

The same oscillator (m=0.50 kgm = 0.50\ \text{kg}, k=8.0 N/mk = 8.0\ \text{N/m}, hence ω=4.0 rad/s\omega = 4.0\ \text{rad/s}) is started with x0=0.030 mx_0 = 0.030\ \text{m} and v0=−0.40 m/sv_0 = -0.40\ \text{m/s}. Write x(t)x(t) and find the first time after t=0t = 0 at which the block passes through equilibrium.

Solution

Find the amplitude from energy (equivalently, from the link between xx and vv):

A=x02+(v0ω)2=0.0302+(−0.404.0)2=9.0×10−4+1.0×10−2≈0.10 m.A = \sqrt{x_0^2 + \left(\frac{v_0}{\omega}\right)^2} = \sqrt{0.030^2 + \left(\frac{-0.40}{4.0}\right)^2} = \sqrt{9.0\times10^{-4} + 1.0\times10^{-2}} \approx 0.10\ \text{m}.

Initial phase from x0=Acos⁡φx_0 = A\cos\varphi: cos⁡φ=0.030/0.104≈0.29\cos\varphi = 0.030/0.104 \approx 0.29, so φ≈±1.28 rad\varphi \approx \pm 1.28\ \text{rad}. Because v0<0v_0 < 0 and v0=−Aωsin⁡φv_0 = -A\omega\sin\varphi, we need sin⁡φ>0\sin\varphi > 0, hence φ=+1.28 rad\varphi = +1.28\ \text{rad}.

Thus x(t)=0.10cos⁡(4.0t+1.28) m.x(t) = 0.10\cos(4.0t + 1.28)\ \text{m}.

Equilibrium means x=0x = 0, i.e. cos⁡(4.0t+1.28)=0\cos(4.0t + 1.28) = 0, so 4.0t+1.28=π/24.0t + 1.28 = \pi/2. Therefore

t=π/2−1.284.0≈0.073 s.t = \frac{\pi/2 - 1.28}{4.0} \approx 0.073\ \text{s}.

At that instant sin⁡(π/2)=1>0\sin(\pi/2) = 1 > 0, so v=−Aω<0v = -A\omega < 0: the block passes equilibrium moving toward negative xx.

Practice

Text

Use the PhET Masses and Springs simulation or an equivalent mass–spring model; do the observing and measuring in the simulation rather than with a real hanging mass.

  1. Set k=8.0 N/mk = 8.0\ \text{N/m} and m=0.50 kgm = 0.50\ \text{kg}. Before running, predict T=2πm/kT = 2\pi\sqrt{m/k} and write the prediction down.
  2. Pull the mass to A=0.040 mA = 0.040\ \text{m}, release it, and record the x(t)x(t) graph. Measure the period from successive crossings of equilibrium and compare it with your prediction.
  3. Repeat with A=0.020 mA = 0.020\ \text{m} and A=0.060 mA = 0.060\ \text{m}. Note whether the period changes with amplitude.
  4. Record the speed at x=0x = 0 and at x=Ax = A; compare with AωA\omega and with zero.
  5. Plot xx, vv and aa over the same time interval and check the phase shifts: vv leads xx by T/4T/4, and aa is opposite to xx.
  6. Change the mass to 2.0 kg2.0\ \text{kg}; predict the new period first, then measure it and compare with 2πm/k2\pi\sqrt{m/k}.

Visual reference

Purpose

Show how the position, velocity and acceleration of a harmonic oscillator are related in time: vv leads xx by a quarter period, and aa is opposite in sign to xx and obeys a=−ω2xa = -\omega^2 x.

Representation

Three vertically stacked graphs sharing one horizontal time axis for one full period TT: x(t)=Acos⁡(ωt)x(t) = A\cos(\omega t), v(t)=−Aωsin⁡(ωt)v(t) = -A\omega\sin(\omega t) and a(t)=−Aω2cos⁡(ωt)a(t) = -A\omega^2\cos(\omega t). A dashed vertical line at t=T/8t = T/8 connects the values on all three graphs; beside them is a small sketch of a mass on a spring with the xx-axis and the equilibrium position marked.

Real-world examples

Car suspension over a speed bump

Scenario

The front axle of a 1400 kg1400\ \text{kg} car carries 700 kg700\ \text{kg}, so each of its two springs supports m=350 kgm = 350\ \text{kg} with stiffness k=2.8×104 N/mk = 2.8\times10^{4}\ \text{N/m}. A speed bump pushes the body down by A=5.0 cmA = 5.0\ \text{cm} and lets go. The angular frequency is ω=k/m=8.9 rad/s\omega=\sqrt{k/m}=8.9\ \text{rad/s}, so the body bounces with period T=2πm/k=0.70 sT=2\pi\sqrt{m/k}=0.70\ \text{s}, while the largest speed and acceleration reach vmax⁡=ωA=0.45 m/sv_{\max}=\omega A=0.45\ \text{m/s} and amax⁡=ω2A=4.0 m/s2a_{\max}=\omega^{2}A=4.0\ \text{m/s}^{2} (about 0.4g0.4g).

Connection

This is a textbook harmonic oscillator: the restoring force is linear, F=−kxF=-kx, and Newton's second law gives a=−kmxa=-\frac{k}{m}x, whose solution is x=Acos⁡(ωt)x=A\cos(\omega t) with ω=k/m\omega=\sqrt{k/m}. The scenario exercises the kinematic pair v=x˙v=\dot x and a=x¨a=\ddot x, which lag a quarter period behind the displacement, and the fact that TT does not depend on the 5.0 cm5.0\ \text{cm} amplitude.

A guitar string

Scenario

The A2 string of a guitar sounds at f=110 Hzf = 110\ \text{Hz}. Its vibrating length is 0.65 m0.65\ \text{m} with linear density 1.5 g/m1.5\ \text{g/m}, giving m=9.8×10−4 kgm = 9.8\times10^{-4}\ \text{kg}; plucked at the middle it swings with amplitude A=1.0 mmA = 1.0\ \text{mm}. Then ω=2πf=691 rad/s\omega = 2\pi f = 691\ \text{rad/s}, the effective stiffness is k=mω2=4.7×102 N/mk = m\omega^{2} = 4.7\times10^{2}\ \text{N/m}, the middle of the string moves at up to vmax⁡=ωA=0.69 m/sv_{\max}=\omega A = 0.69\ \text{m/s}, and its acceleration reaches amax⁡=ω2A=4.8×102 m/s2≈49ga_{\max}=\omega^{2}A = 4.8\times10^{2}\ \text{m/s}^{2}\approx 49g. The stored mechanical energy is E=12kA2=2.3×10−4 JE=\tfrac12 kA^{2}=2.3\times10^{-4}\ \text{J}.

Connection

The plucked string obeys the same kinematics as a mass on a spring: x=Asin⁡(ωt)x = A\sin(\omega t), v=ωAcos⁡(ωt)v = \omega A\cos(\omega t), a=−ω2xa = -\omega^{2}x. The numbers make the phase relations visible — the speed peaks at x=0x=0, while the acceleration peaks at the extreme deflection of 1.0 mm1.0\ \text{mm}, where the speed is zero.

Step counter in a phone

Scenario

A phone in a jogger's pocket records the vertical motion of the torso: the body rises and falls twice per stride, f=2.0 Hzf = 2.0\ \text{Hz}, with amplitude A=3.0 cmA = 3.0\ \text{cm}. The built-in accelerometer samples at 100 Hz100\ \text{Hz} and sees a sinusoidal signal a(t)a(t): ω=2πf=12.6 rad/s\omega = 2\pi f = 12.6\ \text{rad/s} and a peak value amax⁡=ω2A=4.7 m/s2≈0.48ga_{\max}=\omega^{2}A = 4.7\ \text{m/s}^{2}\approx 0.48g. The torso's speed peaks at vmax⁡=ωA=0.38 m/sv_{\max}=\omega A = 0.38\ \text{m/s}.

Connection

The torso is treated as a harmonic oscillator, so its acceleration a=−ω2xa=-\omega^{2}x is proportional to the displacement and always points back towards the middle of the step — that is why the accelerometer trace is sinusoidal and why its amplitude grows as the square of the step frequency. Counting zero crossings of that trace gives the step rate, because one full cycle of x=Acos⁡(ωt)x=A\cos(\omega t) corresponds to one stride.