The Damped, Driven Oscillator
Picture a child on a swing. Left alone after one big push, the swing arcs back and forth but each
arc is a little smaller than the last — air resistance and friction in the chains steadily bleed
away the motion until the swing hangs still. That is damping. Now stand behind the
child and give a small shove once every time the swing comes back to you. If you time your shoves to
the swing's own rhythm, those little pushes add up and up until the child is flying — you have
driven the oscillator on resonance. Mistime them, pushing too fast
or too slow, and you fight the swing as often as you help it; it barely moves.
That single playground scene contains the whole physics of this page. A mass on a spring, an
electron in a radio's tuning circuit, a molecule soaking up light, a suspension bridge in a
crosswind — every one of them is a damped, driven oscillator, and every one obeys
the same equation. We will write that equation down, split its solution into a part that dies away
and a part that lives forever, and discover why a tiny periodic push at just the right frequency can
build a catastrophically large response.
Three forces, one equation
Take a block of mass m on a spring, free to slide along a line. Let
x be its displacement from the resting position. Three forces act on it:
-
the spring's restoring force, -k x — always pulling
back towards x=0, with stiffness k;
-
a damping (friction/drag) force, -b\dot{x} — opposing
the velocity \dot{x}, with damping constant b;
-
an external driving force, F_0\cos(\omega t) — a push
that oscillates at our chosen drive frequency \omega.
Newton's second law, m\ddot{x} = \sum F, gathers them into one line. Move
everything but the drive to the left and you have the equation this whole page is about:
-
The equation of motion.
m\,\ddot{x} + b\,\dot{x} + k\,x = F_0\cos(\omega t).
-
Natural frequency. With no damping and no drive the block oscillates at its own
natural angular frequency
\omega_0 = \sqrt{\dfrac{k}{m}}.
-
Damping rate. Writing \gamma = \dfrac{b}{2m}, the
free motion decays with the envelope e^{-\gamma t}.
It is a second-order, linear, non-homogeneous differential equation: second order
because of the \ddot{x}, linear because x and
its derivatives appear only to the first power, and non-homogeneous because of that
F_0\cos(\omega t) sitting on the right-hand side refusing to be zero.
First, switch the drive off: free damped motion
Before we push, let us just let go. Set F_0 = 0 and the equation becomes
homogeneous:
m\,\ddot{x} + b\,\dot{x} + k\,x = 0.
How the block returns to rest depends entirely on how the damping b
compares with the critical value
b_c = 2\sqrt{mk} = 2m\omega_0. There are three regimes:
-
Underdamped (b < 2\sqrt{mk}): the block overshoots and
oscillates, but inside a shrinking envelope. It rings like a struck bell fading away:
x(t) = A\,e^{-\gamma t}\cos(\omega_d t + \phi), wobbling at a slightly
reduced frequency \omega_d = \sqrt{\omega_0^2 - \gamma^2}.
-
Critically damped (b = 2\sqrt{mk}): the fastest
possible return to rest with no overshoot. This is the sweet spot engineers aim a car's
shock absorbers at — you hit a bump and settle at once, without bouncing.
-
Overdamped (b > 2\sqrt{mk}): so much drag that the
block creeps sluggishly back to zero without ever crossing it — like a door on a stiff hydraulic
closer. It never oscillates.
In every case, though, the key fact is this: the free motion decays. Whatever
wobble the block started with fades under that e^{-\gamma t} envelope and,
after a while, is gone. Hold on to that — it is the reason the driven problem eventually simplifies.
Transient + steady state
Now switch the drive back on. Because the equation is linear and non-homogeneous, its general
solution is a sum of two pieces — exactly the recipe from
solving a
second-order non-homogeneous equation:
x(t) = \underbrace{x_h(t)}_{\text{transient}} + \underbrace{x_p(t)}_{\text{steady state}}.
-
x_h(t) is the homogeneous solution — the free damped
motion we just met. It carries the memory of how the system started (the initial position and
velocity). But it wears the e^{-\gamma t} envelope, so it
dies out. We call it the transient.
-
x_p(t) is a particular solution — a single steady
oscillation that keeps in step with the driving force forever. We call it the
steady state.
Wait long enough — a few decay times — and the transient has vanished. All that survives is the
steady state. And here is the crucial, non-obvious fact:
-
It oscillates at the drive frequency, not the natural one. Once the transient
has died, the block moves as
x_p(t) = A(\omega)\,\cos\!\big(\omega t - \delta\big),
oscillating at \omega — the frequency you impose — not at its
own \omega_0.
-
The amplitude depends on how close you drive to resonance.
A(\omega) = \frac{F_0}{\sqrt{(k - m\omega^2)^2 + (b\omega)^2}}.
-
The response lags the drive by a phase angle
\delta, with
\tan\delta = \dfrac{b\omega}{k - m\omega^2}.
You can check the amplitude formula by substituting x_p = A\cos(\omega t - \delta)
back into the equation of motion; the algebra is a lovely exercise in matching
\cos and \sin terms. What matters for us is the
shape of that function A(\omega) — because that shape is
resonance.
The resonance curve
Let us plot the steady-state amplitude A(\omega) against the drive
frequency \omega. To keep the picture clean we set
m = k = F_0 = 1, so the natural frequency sits at
\omega_0 = \sqrt{k/m} = 1. Drag the slider to change the damping
b and watch the peak.
Everything about resonance is in this one picture. When you drive near
\omega_0, the two terms under the square root nearly cancel — the
(k - m\omega^2) term vanishes exactly at
\omega = \omega_0 — leaving only the small
b\omega term in the denominator, so the amplitude soars. The lighter the
damping, the taller and sharper the spike. This is why timing your pushes to the swing's own rhythm
works so spectacularly, why a wine glass shatters at exactly its ringing note, and why soldiers
break step crossing a bridge.
Almost — but not quite, and it is worth being honest about. It is tempting to say the amplitude is
biggest when you drive precisely at the natural frequency \omega_0. In
fact the peak of A(\omega) sits at a slightly lower frequency,
\omega_{\text{peak}} = \sqrt{\omega_0^2 - 2\gamma^2} = \sqrt{\omega_0^2 - \frac{b^2}{2m^2}},
which is pulled below \omega_0 by the damping. For light
damping (small b) the shift is tiny and "resonance is at
\omega_0" is a fine everyday statement. But crank the damping up and the
peak visibly slides to the left — you can see it happen on the chart above. And a second subtlety:
the free ringing frequency \omega_d = \sqrt{\omega_0^2 - \gamma^2}, the
peak-amplitude frequency \omega_{\text{peak}}, and the natural frequency
\omega_0 are three different numbers once damping is present.
They only collapse to one in the ideal, frictionless limit.
The quality factor Q is a single number that captures how
sharp the resonance is:
Q = \frac{\omega_0}{2\gamma} = \frac{m\omega_0}{b} = \frac{\sqrt{mk}}{b}.
A high Q means light damping: a tall, needle-thin
resonance peak and a system that rings for a long time before its free oscillation dies away
(roughly Q oscillations, in fact). A low
Q means heavy damping: a broad, gentle hump and motion that fades in barely
a swing or two. A struck tuning fork has Q in the thousands; a well-tuned
radio circuit picks one station out of the crowd precisely because its Q is
high enough to respond sharply to one frequency and ignore its neighbours. It is the same knob as the
damping b, just expressed as "how many good oscillations do I get".
Worked examples
Example 1 — the natural frequency. A block of mass
m = 0.5\ \text{kg} hangs on a spring of stiffness
k = 200\ \text{N/m}. Its natural angular frequency is
\omega_0 = \sqrt{\frac{k}{m}} = \sqrt{\frac{200}{0.5}} = \sqrt{400} = 20\ \text{rad/s}.
To make it resonate you would drive it at about 20\ \text{rad/s} — roughly
3.2 pushes per second (since f = \omega_0/2\pi).
Example 2 — which damping regime? The same spring and mass are fitted with a damper
of constant b = 15\ \text{kg/s}. Is the free motion underdamped, critical,
or overdamped? Compare b with the critical value:
b_c = 2\sqrt{mk} = 2\sqrt{0.5 \times 200} = 2\sqrt{100} = 20\ \text{kg/s}.
Since b = 15 < 20 = b_c, the system is underdamped: it
will oscillate a few times inside a decaying envelope before settling. Swap in a stronger damper with
b = 25\ \text{kg/s} > 20 and it would be overdamped, creeping back without
a single wobble.
Example 3 — amplitude at resonance. Keep the underdamped case
(m = 0.5, k = 200,
b = 15) and drive it right at resonance,
\omega \approx \omega_0 = 20\ \text{rad/s}, with a force amplitude
F_0 = 6\ \text{N}. At \omega = \omega_0 the
(k - m\omega^2) term is zero, so the amplitude formula collapses to
A(\omega_0) = \frac{F_0}{b\,\omega_0} = \frac{6}{15 \times 20} = \frac{6}{300} = 0.02\ \text{m} = 2\ \text{cm}.
Halve the damping to b = 7.5 and the resonant amplitude
doubles to 4\ \text{cm} — lighter damping, taller peak, exactly
as the resonance curve promised.