CP 9.8 - Damped and forced oscillations

CP 9.8 - Damped and forced oscillations

Resistance changes both the shape and long-term behaviour of an oscillator's motion. You will model resistance proportional to velocity, distinguish light, critical and heavy damping, and solve forced-motion models by combining a transient response with a steady response.

Modelling resistance and classifying damping

Does adding more resistance always bring a displaced object back sooner? It suppresses oscillations, but sufficiently strong resistance can also make the eventual return slower. The roots of a differential equation explain both effects.

You already know

SHM has restoring acceleration proportional to negative displacement. Constant-coefficient homogeneous equations are solved using an auxiliary equation; a forced equation is solved by adding a particular integral to its complementary function.

Measure displacement xx from equilibrium. A linear restoring force is −kx-kx, while resistance proportional to velocity is −cx˙-c\dot x, with k>0k>0, c≥0c\ge0. The resistance changes direction automatically when velocity changes sign. With an applied force F(t)F(t), Newton's second law gives the model.

Linear damped-motion model

mx¨+cx˙+kx=F(t),m>0,k>0,c≥0m\ddot x+c\dot x+kx=F(t),\qquad m>0,\quad k>0,\quad c\ge0

Here mm is mass, cc the damping coefficient and kk the restoring-force coefficient. Free motion has F(t)=0F(t)=0; forced motion has a specified input.

Not in the formula booklet: learn it.

The model assumes constant coefficients, a linear restoring force and resistance proportional to velocity. Any constant force already balanced at equilibrium disappears when displacement is measured from that equilibrium. Large-amplitude real motion may violate the linear restoring-force assumption.

For free motion, use mr2+cr+k=0mr^2+cr+k=0. With c>0c>0, the sign of c2−4mkc^2-4mk gives three cases.

CaseAuxiliary rootsForm and behaviour
Light damping: c2<4mkc^2<4mk−c/(2m)±iωd-c/(2m)\pm\mathrm{i}\omega_{\mathrm d}Decaying sine/cosine motion
Critical damping: c2=4mkc^2=4mkA repeated negative root(A+Bt)e−ct/(2m)(A+Bt)\mathrm{e}^{-ct/(2m)}
Heavy damping: c2>4mkc^2>4mkTwo distinct negative rootsA sum of two decaying exponentials

Here ωd=k/m−c2/(4m2)\omega_{\mathrm d}=\sqrt{k/m-c^2/(4m^2)} is the light-damping angular frequency. These classification and frequency results are not supplied in the formula booklet. The case c=0c=0 is undamped, not lightly damped.

Key idea

Damping can give decaying oscillation or a non-oscillatory return. Inspect the discriminant before choosing a solution form.

Distinguish a repeated root from a complex pair in this check.

Light damping and a decaying envelope

Complex auxiliary roots retain the oscillation, while their negative real part supplies exponential decay. Include the exponential factor when differentiating for initial velocity.

For x¨+2x˙+5x=0\ddot x+2\dot x+5x=0, the roots are −1±2i-1\pm2\mathrm{i}, giving x=e−t(Acos⁡2t+Bsin⁡2t)x=\mathrm{e}^{-t}(A\cos2t+B\sin2t). The oscillatory factor has amplitude A2+B2\sqrt{A^2+B^2}, so the envelope curves x=±A2+B2 e−tx=\pm\sqrt{A^2+B^2}\,\mathrm{e}^{-t} bound the displacement. Actual extrema need not occur where the trigonometric factor alone is maximal.

Worked example 1: a damped release from rest

Question. Solve x¨+2x˙+5x=0\ddot x+2\dot x+5x=0, given x(0)=1x(0)=1, x˙(0)=0\dot x(0)=0.

Plan. Use the complex-root solution, differentiate its product, then apply both conditions.

x=e−t(Acos⁡2t+Bsin⁡2t)x˙=e−t((−A+2B)cos⁡2t+(−B−2A)sin⁡2t)product ruleA=1,−A+2B=0initial data\begin{aligned} x&=\mathrm{e}^{-t}(A\cos2t+B\sin2t)\\ \dot x&=\mathrm{e}^{-t}\bigl((-A+2B)\cos2t+(-B-2A)\sin2t\bigr) && \text{product rule}\\ A&=1,\qquad -A+2B=0 && \text{initial data} \end{aligned}

Thus B=1/2B=1/2.

Why the term −A-A? Differentiating the exponential contributes it even at t=0t=0. Ignoring that term would change the initial velocity.

Answer. x=e−t(cos⁡2t+12sin⁡2t)x=\mathrm{e}^{-t}(\cos2t+\frac12\sin2t).

Check. The roots satisfy r2+2r+5=0r^2+2r+5=0, and the initial derivative is −1+2(1/2)=0-1+2(1/2)=0.

Worked example 2 changes only the initial velocity. The decay rate and oscillatory frequency belong to the equation, so neither changes.

Worked example 2: an inward initial velocity

Question. Solve the same equation with x(0)=1x(0)=1, x˙(0)=−1\dot x(0)=-1.

Plan. Retain the same general form and change the derivative-condition equation.

A=1initial displacement−A+2B=−1initial velocityB=0solve\begin{aligned} A&=1 && \text{initial displacement}\\ -A+2B&=-1 && \text{initial velocity}\\ B&=0 && \text{solve} \end{aligned}

Why can the sine term disappear? The derivative of the decaying cosine already supplies the required initial inward velocity.

Answer. x=e−tcos⁡2tx=\mathrm{e}^{-t}\cos2t.

Check. x˙(0)=−1\dot x(0)=-1, and ∣x∣≤e−t|x|\le\mathrm{e}^{-t}, so displacement approaches zero while changing sign repeatedly.

The oscillation period of the trigonometric factor is π\pi, but the full displacement is not periodic: successive cycles have smaller sizes. Distinguish the time between corresponding phases from a function repeating exactly.

Watch out

Wrong: set 2B=x˙(0)2B=\dot x(0) as in undamped SHM. Right: here x˙(0)=−A+2B\dot x(0)=-A+2B, because the exponential factor also changes with time.

Use another initial velocity with the same damped equation.

Critical and heavy damping

A repeated negative root gives critical damping; two distinct negative roots give heavy damping. Both responses tend to equilibrium without a continuing oscillation. Initial velocity can still affect whether an individual response crosses equilibrium, so do not infer every detail from the classification alone.

[DIAGRAM: asset_name: Light critical and heavy damping from the same initial data; asset_slug: 9fm0-cp98-damping-comparison; caption: With the same mass and restoring coefficient, increasing damping removes oscillation but can lengthen the late return; description: A time-displacement graph compares solutions of x double dot plus c x dot plus 4x equals zero with x(0)=1 and x dot(0)=0. The light curve c=2 crosses the axis, the critical curve c=4 equals (1+2t) exp(-2t), and the heavy curve c=5 equals (4 exp(-t)-exp(-4t))/3. All approach zero, with the heavy curve above the critical curve at later times.]
With the same mass and restoring coefficient, increasing damping removes oscillation but can lengthen the late return

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account