4.10g - Damped oscillations

4.10g - Damped oscillations

Model a restoring force together with resistance proportional to velocity. You will solve the resulting second order equation and distinguish underdamping, critical damping and overdamping by what their solutions predict.

Building a damping model

Does adding resistance simply make an oscillator swing more slowly, or can it stop repeated swinging altogether? The answer depends on how the resistance compares with the restoring force and the mass.

Suppose displacement xx is measured from equilibrium. A restoring force −kx-kx, with k>0k>0, points towards equilibrium. A resistive force −cx˙-c\dot x, with c>0c>0, opposes velocity. When a particle moves left, x˙<0\dot x<0 and this resistive force points right: resistance changes direction with the motion.

You already know

Newton's second law gives mx¨=Fm\ddot x=F. Without resistance, mx¨=−kxm\ddot x=-kx gives SHM. For a second order homogeneous equation, the auxiliary roots determine whether the solution contains real exponentials, a repeated-root factor, or sine and cosine.

Take a 1 kg1\,\mathrm{kg} particle with restoring force −5x N-5x\,\mathrm{N} and resistance −2x˙ N-2\dot x\,\mathrm{N}, with xx in metres and time in seconds. At x=0.2x=0.2, x˙=−0.5\dot x=-0.5, the restoring force is −1 N-1\,\mathrm{N} and the resistance is +1 N+1\,\mathrm{N}. The acceleration is momentarily zero, although neither force vanishes.

Linear damping model

Damping removes energy from an oscillating system through a resistive effect. For a particle of constant mass m>0m>0, restoring coefficient k>0k>0 and velocity-resistance coefficient c>0c>0, the unforced linear model is

mx¨+cx˙+kx=0m\ddot x+c\dot x+kx=0

It assumes the two force laws remain proportional over the motion and no external driving force acts. Resistance −cx˙-c\dot x is not a restoring force: it depends on velocity, not position.

All coefficients here are constants. The mass has units kg\mathrm{kg}, cc has units kg s−1\mathrm{kg\,s^{-1}}, and kk has units kg s−2\mathrm{kg\,s^{-2}}, so all three terms are forces. This model equation is not in the formula booklet; learn how to construct it from the stated forces.

Worked example 1: forces with different jobs

Question. A 1 kg1\,\mathrm{kg} particle has restoring force −5x N-5x\,\mathrm{N} and resistive force −2x˙ N-2\dot x\,\mathrm{N}. Form its equation of motion and find its acceleration when x=0.2 mx=0.2\,\mathrm{m}, x˙=−0.5 m s−1\dot x=-0.5\,\mathrm{m\,s^{-1}}.

Plan. Add the signed forces, then substitute the position and velocity into the acceleration equation.

x¨=−5x−2x˙F=ma, m=1x¨+2x˙+5x=0standard formx¨=−5(0.2)−2(−0.5)=0substitute both signed values\begin{aligned} \ddot x&=-5x-2\dot x && F=ma,\ m=1\\ \ddot x+2\dot x+5x&=0 && \text{standard form}\\ \ddot x&=-5(0.2)-2(-0.5)=0 && \text{substitute both signed values} \end{aligned}

Why retain the sign of velocity? The resistive force must reverse when the motion reverses.

Answer. x¨+2x˙+5x=0\ddot x+2\dot x+5x=0; the acceleration at that instant is 0 m s−20\,\mathrm{m\,s^{-2}}.

Check. The two forces are −1-1 and +1 N+1\,\mathrm{N}, giving zero resultant.

Change only the resistance coefficient from 22 to 44. The same negative velocity now produces a larger positive resistive force.

Worked example 2: stronger resistance

Question. The same 1 kg1\,\mathrm{kg} particle now has forces −5x N-5x\,\mathrm{N} and −4x˙ N-4\dot x\,\mathrm{N}. Form its equation and find the acceleration at x=0.2x=0.2, x˙=−0.5\dot x=-0.5, in SI units.

Plan. Replace the damping coefficient while keeping the restoring force and mass unchanged.

x¨=−5x−4x˙resultant force divided by massx¨+4x˙+5x=0standard formx¨=−5(0.2)−4(−0.5)=1substitute\begin{aligned} \ddot x&=-5x-4\dot x && \text{resultant force divided by mass}\\ \ddot x+4\dot x+5x&=0 && \text{standard form}\\ \ddot x&=-5(0.2)-4(-0.5)=1 && \text{substitute} \end{aligned}

Why can the acceleration point away from equilibrium? Resistance opposes the present velocity; at this instant it exceeds the restoring force. That does not make resistance a source of energy.

Answer. x¨+4x˙+5x=0\ddot x+4\dot x+5x=0, with acceleration +1 m s−2+1\,\mathrm{m\,s^{-2}}.

Check. The resultant is −1+2=1 N-1+2=1\,\mathrm{N}, agreeing with mx¨m\ddot x.

The sign of each force follows its own physical role, rather than a memorised direction on the page.

Watch out

Wrong: write resistance as −c∣x˙∣-c|\dot x|, always pointing left. Right: the signed component is −cx˙-c\dot x, so it always opposes motion.

Use a new resistance coefficient in the same force model.

Underdamping

With modest resistance, the particle still crosses equilibrium repeatedly, but its oscillations shrink.

For x¨+2x˙+5x=0\ddot x+2\dot x+5x=0, the roots are −1±2i-1\pm2i, so the solution is a sine–cosine oscillation multiplied by e−te^{-t}. Its oscillation has angular frequency 22, whereas the undamped model x¨+5x=0\ddot x+5x=0 has angular frequency 5\sqrt5: damping affects timing as well as size.

Underdamped solution

For mx¨+cx˙+kx=0m\ddot x+c\dot x+kx=0, with m,k,c>0m,k,c>0 and c2<4mkc^2<4mk, the auxiliary roots are

−c2m±iβ,β=4mk−c22m-\frac{c}{2m}\pm i\beta,\qquad \beta=\frac{\sqrt{4mk-c^2}}{2m}

The solution has the form

x=e−ct/(2m)(Acos⁡βt+Bsin⁡βt)x=e^{-ct/(2m)}(A\cos\beta t+B\sin\beta t)

This is underdamping: a non-zero solution oscillates with a decreasing envelope. It repeats its phase every 2π/β2\pi/\beta, but its displacement is not periodic because its size decreases.
Not in the formula booklet: learn it.

You can find the roots directly instead of memorising their expanded formula. The sine–cosine bracket can be written Rsin⁡(βt+ϕ)R\sin(\beta t+\phi); hence ∣x∣≤Re−ct/(2m)|x|\le Re^{-ct/(2m)}. The bounding curves need not meet the actual local maxima, because the exponential is changing while the particle moves.

Worked example 3: release from rest with damping

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. Describe the motion as tt increases.

Plan. Use the complex auxiliary roots and differentiate the full product before applying the velocity condition.

m2+2m+5=0⇒m=−1±2iauxiliary equationx=e−t(Acos⁡2t+Bsin⁡2t)complex rootsx(0)=A=1displacement conditionx˙=e−t[(−A+2B)cos⁡2t+(−B−2A)sin⁡2t]product rulex˙(0)=−1+2B=0⇒B=12velocity condition\begin{aligned} m^2+2m+5=0&\Rightarrow m=-1\pm2i && \text{auxiliary equation}\\ x&=e^{-t}(A\cos2t+B\sin2t) && \text{complex roots}\\ x(0)=A&=1 && \text{displacement condition}\\ \dot x&=e^{-t}[(-A+2B)\cos2t+(-B-2A)\sin2t] && \text{product rule}\\ \dot x(0)=-1+2B=0&\Rightarrow B=\frac12 && \text{velocity condition} \end{aligned}

Why differentiate the exponential too? It contributes −A-A to the initial velocity, so omitting it gives the wrong sine coefficient.

Answer. x=e−t(cos⁡2t+12sin⁡2t)x=e^{-t}(\cos2t+\tfrac12\sin2t). The motion is underdamped: it crosses equilibrium repeatedly and tends to equilibrium with shrinking oscillations.

Check. x(0)=1x(0)=1, x˙(0)=−1+1=0\dot x(0)=-1+1=0. The bound ∣x∣≤52e−t|x|\le\tfrac{\sqrt5}{2}e^{-t} tends to zero.

Change only the initial velocity to −1-1. The equation and decay rate stay fixed, but the initial phase changes.

Worked example 4: a moving start

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

Plan. Keep the same general solution and use the changed velocity condition to find BB.

x=e−t(Acos⁡2t+Bsin⁡2t)m=−1±2iA=1x(0)=1−A+2B=−1x˙(0)=−1B=0solve\begin{aligned} x&=e^{-t}(A\cos2t+B\sin2t) && m=-1\pm2i\\ A&=1 && x(0)=1\\ -A+2B&=-1 && \dot x(0)=-1\\ B&=0 && \text{solve} \end{aligned}

Why is a cosine-only bracket compatible with non-zero initial velocity? The decaying multiplier itself is changing at t=0t=0.

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

Check. x˙=e−t(−cos⁡2t−2sin⁡2t)\dot x=e^{-t}(-\cos2t-2\sin2t), so x˙(0)=−1\dot x(0)=-1; the graph lies between ±e−t\pm e^{-t}.

Continue this lesson.

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

Create free account