Unit 6, lesson 47

Oscillations, Damping, and Resonance

How Newton's law turns a spring, a damper and an applied force into one second-order equation, what its free, damped and forced solutions look like, and exactly when resonance occurs.

  • Theory part 93
  • Practice part 94
  • Second-Order Linear Equations

A car hits a bump and bounces once before it settles. A wine glass shatters when a singer holds the right note. A radio picks one station out of many. All three come from the same equation with different numbers in it, and solving that equation explains why each one behaves the way it does.

Key terms

Spring-mass system

Definition A mass m>0m > 0 hangs on a spring with spring constant k>0k > 0. Measure the displacement x(t)x(t) from the equilibrium position, with downward positive. Hooke’s law says the spring pulls with force −k-k times its stretch. At equilibrium the spring is stretched by ℓ\ell with kℓ=mgk\ell = mg, so gravity and the spring balance. With displacement xx the stretch is ℓ+x\ell + x, and the net of gravity and spring is mg−k(ℓ+x)=−kxmg - k(\ell + x) = -kx. Adding a damping force −c x′-c\,x' with c≥0c \ge 0 and an applied force F(t)F(t), Newton’s second law gives

mx′′+cx′+kx=F(t).m x'' + c x' + k x = F(t).

In plain words Measuring from the resting position makes gravity disappear from the equation. What is left is the mass times acceleration on one side and the spring, the damper and the outside push on the other.

Example A mass of 0.50.5 kilograms stretches a spring 0.0980.098 meters at rest. With g=9.8g = 9.8 meters per second squared,

k=mgℓ=0.5⋅9.80.098=50 newtons per meter.\begin{aligned} k &= \frac{mg}{\ell} = \frac{0.5 \cdot 9.8}{0.098} \\[0.4em] &= 50 \text{ newtons per meter.} \end{aligned}

Without damping or applied force the model is 0.5 x′′+50 x=00.5\,x'' + 50\,x = 0, that is x′′+100x=0x'' + 100x = 0.

Simple harmonic motion

Definition When c=0c = 0 and F=0F = 0, the equation mx′′+kx=0m x'' + k x = 0 has the general solution

x=c1cos⁡ω0t+c2sin⁡ω0t,ω0=km.\begin{gathered} x = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t, \\[0.4em] \omega_0 = \sqrt{\frac{k}{m}}. \end{gathered}

With R=c12+c22R = \sqrt{c_1^2 + c_2^2} and an angle δ\delta such that c1=Rcos⁡δc_1 = R\cos\delta and c2=Rsin⁡δc_2 = R\sin\delta, this is x=Rcos⁡(ω0t−δ)x = R\cos(\omega_0 t - \delta). The number RR is the amplitude, ω0\omega_0 is the natural angular frequency in radians per second, and the period is 2π/ω02\pi/\omega_0.

In plain words With no friction the mass swings back and forth forever with the same size and the same rhythm. The rhythm depends only on the mass and the stiffness, never on how far you pull it.

Example For the spring above, ω0=10\omega_0 = 10. Suppose x(0)=0.03x(0) = 0.03 meters and x′(0)=−0.4x'(0) = -0.4 meters per second. Then c1=0.03c_1 = 0.03 and 10c2=−0.410c_2 = -0.4, so c2=−0.04c_2 = -0.04. The amplitude and period are

R=0.032+0.042=0.05 meters,T=2π10≈0.63 seconds.\begin{aligned} R &= \sqrt{0.03^2 + 0.04^2} \\[0.4em] &= 0.05 \text{ meters}, \\[0.4em] T &= \frac{2\pi}{10} \approx 0.63 \text{ seconds.} \end{aligned}

The phase satisfies cos⁡δ=0.6\cos\delta = 0.6 and sin⁡δ=−0.8\sin\delta = -0.8, so δ≈−0.927\delta \approx -0.927.

Damping

Definition With damping constant c>0c > 0 and no applied force, the characteristic equation mr2+cr+k=0m r^2 + c r + k = 0 has roots

r=−c±c2−4mk2m.r = \frac{-c \pm \sqrt{c^2 - 4mk}}{2m}.

Since mm, cc and kk are positive, both roots have negative real part, so every solution tends to 00 as t→∞t \to \infty. The sign of c2−4mkc^2 - 4mk sorts the motion into three cases, treated next.

In plain words A damper, such as a shock absorber, always takes energy out. How the motion dies away depends on whether the damping is weak or strong compared with the spring and the mass.

Example For m=1m = 1 and k=4k = 4, the critical damping value comes from c2=4mk=16c^2 = 4mk = 16, so c=4c = 4.

Damping constantc2−4mkc^2 - 4mkCase
c=2c = 2−12-12underdamped
c=4c = 400critically damped
c=5c = 599overdamped

Underdamped motion

Definition If 0<c2<4mk0 < c^2 < 4mk, the roots are −c2m±iμ-\frac{c}{2m} \pm i\mu with μ=4mk−c22m>0\mu = \frac{\sqrt{4mk - c^2}}{2m} > 0, and

x=e−ct/(2m)(c1cos⁡μt+c2sin⁡μt).x = e^{-ct/(2m)}\big(c_1 \cos \mu t + c_2 \sin \mu t\big).

The mass crosses equilibrium infinitely often, at time gaps of π/μ\pi/\mu, and the oscillation is squeezed by the envelope ±R e−ct/(2m)\pm R\,e^{-ct/(2m)}. The number 2π/μ2\pi/\mu is called the quasi period.

In plain words It still oscillates, a bit more slowly than without damping, and each swing is smaller than the one before.

Example Take x′′+2x′+10x=0x'' + 2x' + 10x = 0 with x(0)=1x(0) = 1 and x′(0)=0x'(0) = 0. The roots are −1±3i-1 \pm 3i, so x=e−t(c1cos⁡3t+c2sin⁡3t)x = e^{-t}(c_1 \cos 3t + c_2 \sin 3t). The data give c1=1c_1 = 1 and −c1+3c2=0-c_1 + 3c_2 = 0, so

x=e−t(cos⁡3t+13sin⁡3t).x = e^{-t}\Big(\cos 3t + \tfrac{1}{3}\sin 3t\Big).

The quasi period is 2π/3≈2.092\pi/3 \approx 2.09 seconds.

Critically damped motion

Definition If c2=4mkc^2 = 4mk, the characteristic equation has the double root −c2m-\frac{c}{2m}, and

x=(c1+c2t) e−ct/(2m).x = (c_1 + c_2 t)\,e^{-ct/(2m)}.

Unless xx is identically zero, the factor c1+c2tc_1 + c_2 t has at most one zero, so the mass passes through equilibrium at most once.

In plain words This is the smallest damping that stops all oscillation. A well tuned door closer aims for it, returning fast without swinging past.

Example Take x′′+4x′+4x=0x'' + 4x' + 4x = 0 with x(0)=1x(0) = 1 and x′(0)=0x'(0) = 0. The double root is −2-2. The data give c1=1c_1 = 1 and −2c1+c2=0-2c_1 + c_2 = 0, so

x=(1+2t) e−2t,x = (1 + 2t)\,e^{-2t},

which stays positive for t≥0t \ge 0 and decreases to 00.

Overdamped motion

Definition If c2>4mkc^2 > 4mk, the roots r1r_1 and r2r_2 are real, distinct and negative, and

x=c1er1t+c2er2t.x = c_1 e^{r_1 t} + c_2 e^{r_2 t}.

As with critical damping, a solution that is not identically zero passes through equilibrium at most once, since c1er1t=−c2er2tc_1 e^{r_1 t} = -c_2 e^{r_2 t} has at most one solution tt. Order the roots so that r2<r1<0r_2 < r_1 < 0. For fixed mm and kk, the slower root satisfies

−ω0<r1<0,ω0=km.-\omega_0 < r_1 < 0, \qquad \omega_0 = \sqrt{\frac{k}{m}}.

So every solution with c1≠0c_1 \ne 0 decays like er1te^{r_1 t}, more slowly in the long run than every nonzero critically damped solution, which decays like e−ω0te^{-\omega_0 t} or t e−ω0tt\,e^{-\omega_0 t}. Initial data that make c1=0c_1 = 0 excite only the fast root and decay faster.

In plain words The damping is so heavy that the mass creeps back to rest. For typical starting data it takes longer to settle than with critical damping, because the heavy damper also resists the spring pulling it home.

Example Take x′′+5x′+4x=0x'' + 5x' + 4x = 0 with x(0)=1x(0) = 1 and x′(0)=0x'(0) = 0. The roots are −1-1 and −4-4. The equations c1+c2=1c_1 + c_2 = 1 and −c1−4c2=0-c_1 - 4c_2 = 0 give c1=43c_1 = \tfrac{4}{3} and c2=−13c_2 = -\tfrac{1}{3}, so

x=43e−t−13e−4t.x = \tfrac{4}{3}e^{-t} - \tfrac{1}{3}e^{-4t}.

This has the same m=1m = 1 and k=4k = 4 as the critical example, where ω0=2\omega_0 = 2. Here c1≠0c_1 \ne 0 and the slow rate is 1<21 < 2. At t=3t = 3 the overdamped mass is still at about 0.0660.066, while the critically damped one is at 7e−6≈0.0177e^{-6} \approx 0.017.

Forcing

Definition A periodic applied force F(t)=F0cos⁡ωtF(t) = F_0 \cos \omega t, with F0≠0F_0 \ne 0 and ω>0\omega > 0, gives mx′′+cx′+kx=F0cos⁡ωtm x'' + c x' + k x = F_0 \cos\omega t. In the undamped case c=0c = 0 with ω≠ω0\omega \ne \omega_0, undetermined coefficients gives the particular solution

xp=F0m(ω02−ω2)cos⁡ωt,x_p = \frac{F_0}{m(\omega_0^2 - \omega^2)}\cos \omega t,

and the general solution is xpx_p plus simple harmonic motion at frequency ω0\omega_0.

In plain words The outside push drives a motion at the frequency of the push, and that forced motion sits on top of the natural swinging at ω0\omega_0.

Example For x′′+9x=16cos⁡5tx'' + 9x = 16\cos 5t we have ω0=3\omega_0 = 3, ω=5\omega = 5, and xp=169−25cos⁡5t=−cos⁡5tx_p = \frac{16}{9 - 25}\cos 5t = -\cos 5t. Starting from rest at equilibrium, x(0)=0x(0) = 0 and x′(0)=0x'(0) = 0 give

x=cos⁡3t−cos⁡5t=2sin⁡t sin⁡4t.x = \cos 3t - \cos 5t = 2\sin t \,\sin 4t.

The fast factor sin⁡4t\sin 4t is wrapped in the slow envelope 2sin⁡t2\sin t. This pattern is called beats.

Resonance

Definition Pure resonance happens when c=0c = 0 and ω=ω0\omega = \omega_0. Then cos⁡ω0t\cos\omega_0 t solves the homogeneous equation, the trial needs a factor tt, and

xp=F02mω0 tsin⁡ω0t,x_p = \frac{F_0}{2m\omega_0}\,t \sin \omega_0 t,

which is unbounded as t→∞t \to \infty. When c>0c > 0 every solution stays bounded, and for F0>0F_0 > 0 the steady periodic solution has amplitude

A(ω)=F0(k−mω2)2+c2ω2.A(\omega) = \frac{F_0}{\sqrt{(k - m\omega^2)^2 + c^2\omega^2}}.

If c2<2mkc^2 < 2mk, then on ω>0\omega > 0 the amplitude A(ω)A(\omega) has its largest value at

ωr=km−c22m2.\omega_r = \sqrt{\frac{k}{m} - \frac{c^2}{2m^2}}.

This peak is called practical resonance, and its frequency ωr\omega_r is less than ω0\omega_0. If c2≥2mkc^2 \ge 2mk, then A(ω)A(\omega) strictly decreases for ω>0\omega > 0 and there is no peak. Both facts follow by writing u=ω2u = \omega^2 and differentiating (k−mu)2+c2u(k - mu)^2 + c^2 u, which gives 2m2u−2mk+c22m^2 u - 2mk + c^2.

In plain words Pushing an undamped system at exactly its natural rhythm adds energy every cycle without limit. Real systems have some damping, so the response peaks but stays finite.

Example For x′′+4x=3cos⁡2tx'' + 4x = 3\cos 2t, here m=1m = 1 and ω0=2=ω\omega_0 = 2 = \omega, so

xp=32⋅1⋅2 tsin⁡2t=34 tsin⁡2t.x_p = \frac{3}{2 \cdot 1 \cdot 2}\,t\sin 2t = \tfrac{3}{4}\,t \sin 2t.

Add damping c=1c = 1. Since c2=1<8=2mkc^2 = 1 < 8 = 2mk, the peak is at ω2=4−12=3.5\omega^2 = 4 - \tfrac{1}{2} = 3.5, where A≈1.549A \approx 1.549. At ω=2\omega = 2 the amplitude is exactly 1.51.5.

RLC circuits

Definition A series circuit with inductance L>0L > 0 henries, resistance R≥0R \ge 0 ohms, capacitance C>0C > 0 farads and applied voltage E(t)E(t) volts satisfies, for the charge q(t)q(t) in coulombs on the capacitor,

Lq′′+Rq′+1C q=E(t),L q'' + R q' + \frac{1}{C}\,q = E(t),

and the current is I=q′I = q' amperes. It is the spring-mass equation with LL for mm, RR for cc and 1/C1/C for kk, so every result above carries over.

In plain words Inductance resists changes in current like mass resists changes in velocity, resistance burns energy like a damper, and the capacitor stores it like a spring.

Example Take L=0.5L = 0.5, R=10R = 10, C=0.01C = 0.01 and a constant E=12E = 12 volts, with q(0)=0q(0) = 0 and I(0)=0I(0) = 0. Multiplying by 22 gives q′′+20q′+200q=24q'' + 20q' + 200q = 24, with roots −10±10i-10 \pm 10i and constant particular solution qp=0.12q_p = 0.12. Fitting the data gives

q=0.12−0.12 e−10t⋅(cos⁡10t+sin⁡10t).\begin{aligned} q = 0.12 &- 0.12\,e^{-10t} \\[0.4em] &\cdot (\cos 10t + \sin 10t). \end{aligned}

As a unit check, the final charge satisfies q/C=0.12/0.01=12q/C = 0.12/0.01 = 12 volts, matching the battery.

Transient and steady response

Definition Suppose m,c,k>0m, c, k > 0 and F(t)=F0cos⁡ωtF(t) = F_0\cos\omega t. Every solution is x=xc+xpx = x_c + x_p, where xcx_c solves the homogeneous equation and tends to 00 as t→∞t \to \infty. The part xcx_c is the transient. The periodic particular solution

xp=A(ω)cos⁡(ωt−γ)x_p = A(\omega)\cos(\omega t - \gamma)

is the steady periodic response. Here F0>0F_0 > 0, A(ω)A(\omega) is the amplitude from the resonance section, and the phase lag γ\gamma is the angle in (0,π)(0, \pi) with

cos⁡γ=k−mω2(k−mω2)2+c2ω2,sin⁡γ=c ω(k−mω2)2+c2ω2.\begin{aligned} \cos\gamma &= \frac{k - m\omega^2}{\sqrt{(k - m\omega^2)^2 + c^2\omega^2}}, \\[0.4em] \sin\gamma &= \frac{c\,\omega}{\sqrt{(k - m\omega^2)^2 + c^2\omega^2}}. \end{aligned}

So the response lags the force by less than a quarter period when ω<ω0\omega < \omega_0, by exactly π/2\pi/2 when ω=ω0\omega = \omega_0, and by more when ω>ω0\omega > \omega_0. The steady response does not depend on the initial conditions.

In plain words Whatever the start, the memory of it fades, and the motion settles into the rhythm of the push.

Example For x′′+2x′+2x=5cos⁡tx'' + 2x' + 2x = 5\cos t the trial acos⁡t+bsin⁡ta\cos t + b\sin t gives a+2b=5a + 2b = 5 and b−2a=0b - 2a = 0, so a=1a = 1, b=2b = 2 and

xp=cos⁡t+2sin⁡t,x_p = \cos t + 2\sin t,

with amplitude 5\sqrt{5}, matching A(1)=5/1+4A(1) = 5/\sqrt{1 + 4}. The phase lag has tan⁡γ=2\tan\gamma = 2, so γ≈1.107\gamma \approx 1.107. The transient is e−t(c1cos⁡t+c2sin⁡t)e^{-t}(c_1\cos t + c_2\sin t), and after a few seconds it is too small to see.

Check yourself

  1. A system has m=2m = 2 and k=8k = 8. Which damping constant makes it critically damped, and is c=6c = 6 under or over that?
Show the answer

Critical damping needs c2=4mk=64c^2 = 4mk = 64, so c=8c = 8. With c=6c = 6 we get c2−4mk=36−64<0c^2 - 4mk = 36 - 64 < 0, so the system is underdamped.

  1. For x′′+25x=10cos⁡ωtx'' + 25x = 10\cos\omega t, which ω\omega gives pure resonance, and what is the particular solution then?
Show the answer

The natural frequency is ω0=5\omega_0 = 5, so resonance occurs at ω=5\omega = 5. The trial t(Acos⁡5t+Bsin⁡5t)t(A\cos 5t + B\sin 5t) gives A=0A = 0 and B=1B = 1, so xp=tsin⁡5tx_p = t\sin 5t, which grows without bound.

  1. An RLC circuit has L=1L = 1, R=2R = 2, 1/C=51/C = 5, E=0E = 0, q(0)=1q(0) = 1 and q′(0)=−1q'(0) = -1. Find q(t)q(t).
Show the answer

The equation q′′+2q′+5q=0q'' + 2q' + 5q = 0 has roots −1±2i-1 \pm 2i, so q=e−t(c1cos⁡2t+c2sin⁡2t)q = e^{-t}(c_1\cos 2t + c_2\sin 2t). The data give c1=1c_1 = 1 and −1+2c2=−1-1 + 2c_2 = -1, so c2=0c_2 = 0 and q=e−tcos⁡2tq = e^{-t}\cos 2t.

How this lesson connects

It also uses Second-order constant-coefficient linear equations, Trigonometry and initial conditions and Force, displacement and elementary circuit quantities introduced in the lesson.