Differential Equations · Damping

Damping

Damping removes energy from an oscillation. Adjust ζ to compare oscillating returns with returns that do not cross equilibrium. Time is measured in T₀, the period without damping.

01

Damping in this visualization

Damping is the loss of energy that makes an oscillation die out. For a mass on a spring with a damper, mx′′+cx′+kx=0m x'' + c x' + k x = 0, written x′′+2ζω0x′+ω02x=0x'' + 2\zeta\omega_0 x' + \omega_0^2 x = 0 with the natural frequency ω0=k/m\omega_0 = \sqrt{k/m} and the damping ratio ζ=c/(2mk)\zeta = c/(2\sqrt{mk}). Below 1 the motion is underdamped (it swings), at 1 critically damped, above 1 overdamped (it creeps). For a mass hanging vertically, gravity only moves the rest position down by mg/k; x is measured from there. The page releases the mass from rest at x = 1 and asks whether more damping stops it sooner.

The roots determine the motion; a repeated root also gives a t·e^(λt) term
λ2+2ζω0λ+ω02=0,λ=ω0(−ζ±ζ2−1),σ=−max⁡Re⁡λ\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0, \qquad \lambda = \omega_0\left(-\zeta \pm \sqrt{\zeta^2 - 1}\right), \qquad \sigma = -\max \operatorname{Re}\lambda

The right of the picture is the root plane. The two roots multiply to ω02\omega_0^2, so while the motion swings they lie on the circle of radius ω₀: the real part −ζω0-\zeta\omega_0 sets the decay, the imaginary part ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2} the frequency, so damping also slows each swing. At ζ = 1 the roots meet at −ω₀. Past it one root runs left and the other creeps back towards 0, about −ω0/(2ζ)=−k/c-\omega_0/(2\zeta) = -k/c for large ζ. The slowest part of the motion decays like e−σte^{-\sigma t}, and this decay rate σ, the distance of the rightmost root from the imaginary axis, is largest exactly at ζ = 1. On x(t) the dashed blue curve shows it: the envelope ±e−ζω0t/1−ζ2\pm e^{-\zeta\omega_0 t}/\sqrt{1-\zeta^2} when underdamped (the curve touches it between its peaks), the slow root's term when overdamped.

“Stops” needs a rule. The page uses T1%T_{1\%}, the earliest time after which |x| stays within 1 % of the release for good. Its strip under x(t) magnifies the band. The extrema sit at tn=nπ/ωdt_n = n\pi/\omega_d with ∣x(tn)∣=e−ζω0tn|x(t_n)| = e^{-\zeta\omega_0 t_n}, which decides exactly which peaks leave the band; the last return is then found by bisection on the closed-form solution, so the readout is numeric. A smaller ζ gives a slightly later settling time whenever one more peak leaves the band, so the curve of T1%T_{1\%} against ζ jumps. Its minimum is not at critical damping but at ζ∗=ln⁡100/π2+ln⁡2100≈0.826\zeta^* = \ln 100/\sqrt{\pi^2 + \ln^2 100} \approx 0.826, where the one overshoot touches −1 % exactly: 0.667 T₀ against 1.057 T₀. Critical damping is the fastest only among motions that never cross x = 0, and only from rest; another band moves ζ* (a 2 % band gives about 0.780).

02

Things to notice

  • Two different “fastest”: the decay rate σ peaks at ζ = 1 exactly; the settling time is shortest at ζ*, a little below. A literal 0.826 is not ζ*: it overshoots by 1.0015 %, and its settling time jumps to 0.896 T₀. The page keeps ζ* exact.
  • Doubling a large damping ratio roughly doubles the time to creep back: the slow root's time constant is about c/k.
  • x(t) is the exact solution in all three regimes; the regime of a typed ratio is decided on the number as typed, so 0.99999999999999999 is underdamped. The fluid in the cylinder is only a sketch.
03

Controls

  • Push the ζ fader (it catches at 0.1, 1 and 10), type a fraction or decimal, or press ζ = 1 or ζ*. A root can be dragged along its path; at −ω₀ it catches, and only there does it change branch. Drag the red time cursor to move the mass.
  • Presentation mode opens a full-screen view: ← → change ζ by a fader step, H hides the decay rate and settling times.
04

Related

On the trace–determinant map of the phase portrait page, the damped oscillator is the horizontal line Δ = ω₀².

Further reading: Wikipedia: Damping; OpenStax, University Physics Volume 1, 15.5 Damped Oscillations.