Technopedia

The technologist’s encyclopedia

Driven Oscillations

Driven oscillations describe a system subjected to an external periodic force. A distinct but related mechanism is parametric excitation, in which a system parameter (stiffness, length, inductance) is modulated periodically rather than a force being applied. Together they explain resonance, the frequency response of filters, and the operation of parametric amplifiers. They are relevant to mechanical resonance, electrical tuning, bridge dynamics, and MEMS devices.

k c m F₀ cos(ωt)

Authoritative sources

Equation of motion (force-driven)

mx¨+cx˙+kx=F0cos⁡(ωt)m\ddot{x} + c\dot{x} + kx = F_0\cos(\omega t)

where cc is the damping coefficient, kk the stiffness, ω0=k/m\omega_0 = \sqrt{k/m} the natural frequency, and γ=c/(2m)\gamma = c/(2m) the damping rate. The drive has amplitude F0F_0 and angular frequency ω\omega.

Details

Decomposition into transient + steady-state. Because the equation is linear, the full solution is the sum of two independent parts:

x(t)=xtrans(t)+xss(t)x(t) = x_{\text{trans}}(t) + x_{ss}(t)

After a few damping time constants (t≳1/γt \gtrsim 1/\gamma), the transient is gone and only the steady-state remains. This is why a driven system “locks” to the drive frequency: whatever the initial motion, the long-time motion is a pure sinusoid at ω\omega.

t x x(t) total steady-state transient (decays)

Authoritative source for this section: Landau & Lifshitz, Mechanics, §14: Forced Oscillations

Example problem

A driven oscillator has γ=0.5\gamma = 0.5 s⁻¹. How long until the transient amplitude has fallen to 5% of its initial value?

Solution

The transient envelope is e−γte^{-\gamma t}. Set e−0.5 t=0.05e^{-0.5\,t} = 0.05: t=ln⁡(20)/0.5≈2.996/0.5≈6.0t = \ln(20)/0.5 \approx 2.996/0.5 \approx 6.0 s.

Steady-state response

The steady-state solution oscillates at the drive frequency:

xss(t)=X(ω)cos⁡(ωt−δ)x_{ss}(t) = X(\omega)\cos(\omega t - \delta)

where the amplitude and phase lag are:

X(ω)=F0/m(ω02−ω2)2+(2γω)2,tan⁡δ=2γωω02−ω2X(\omega) = \frac{F_0/m}{\sqrt{(\omega_0^2 - \omega^2)^2 + (2\gamma\omega)^2}}, \qquad \tan\delta = \frac{2\gamma\omega}{\omega_0^2 - \omega^2}
Details

Three frequency regimes.

Amplitude response curve. The plot of X(ω)X(\omega) versus ω\omega is the frequency response of the oscillator: a resonance peak centered near ω0\omega_0. The width of the peak at half its maximum height is the bandwidth Δω\Delta\omega.

ω X(ω) Δω peak ω₀

Phase-lag curve. The phase lag δ(ω)\delta(\omega) rises smoothly from 00 to π\pi, passing through π/2\pi/2 exactly at ω=ω0\omega = \omega_0.

ω δ π/2 π 0 ω₀

Authoritative source for this section: OpenStax University Physics Volume 1, §15.6: Forced Oscillations and Resonance

Example problem

For ω0=10\omega_0 = 10 rad/s, γ=0.5\gamma = 0.5 s⁻¹, F0/m=2F_0/m = 2 m/s², compute XX at ω=2ω0=20\omega = 2\omega_0 = 20 rad/s.

Solution

X=2(100−400)2+(2×0.5×20)2=2(−300)2+202=290000+400=290400≈2300.7≈0.00665X = \dfrac{2}{\sqrt{(100 - 400)^2 + (2 \times 0.5 \times 20)^2}} = \dfrac{2}{\sqrt{(-300)^2 + 20^2}} = \dfrac{2}{\sqrt{90000 + 400}} = \dfrac{2}{\sqrt{90400}} \approx \dfrac{2}{300.7} \approx 0.00665 m.

Resonance and quality factor

The peak amplitude (at the displacement resonance ωres\omega_{res}) is, in the weak-damping limit:

Xmax⁡=F0c ω0=F02mγ ω0=Q F0kX_{\max} = \frac{F_0}{c\,\omega_0} = \frac{F_0}{2m\gamma\,\omega_0} = Q\,\frac{F_0}{k}
Details

Displacement vs. velocity resonance. The amplitude X(ω)X(\omega) is maximized at

ωres=ω02−2γ2≈ω0(γ≪ω0)\omega_{res} = \sqrt{\omega_0^2 - 2\gamma^2} \approx \omega_0 \quad (\gamma \ll \omega_0)

not exactly at ω0\omega_0. The displacement peak only coincides with ω0\omega_0 in the weak-damping limit. By contrast, the velocity (and hence the power) response is maximized exactly at ω=ω0\omega = \omega_0, because there the velocity is in phase with the driving force. For high-QQ systems the two frequencies are nearly identical, which is why the distinction is often glossed over.

Quality factor. Q=ω0/(2γ)Q = \omega_0/(2\gamma) quantifies the sharpness of resonance. A high-QQ system (tuning fork, Q∼103Q \sim 10^3; superconducting cavity, Q∼1010Q \sim 10^{10}) has a narrow peak and stores energy for many cycles; a low-QQ system (shock absorber, Q∼1Q \sim 1) has a broad, flat response. Equivalently, QQ is the number of radians of phase over which the free (undriven) oscillation decays by ee; see Damped oscillations for the free-decay analysis.

Bandwidth. The full width at half maximum of the amplitude peak is

Δω≈2γ=ω0Q,i.e.Δωω0≈1Q\Delta\omega \approx 2\gamma = \frac{\omega_0}{Q}, \qquad \text{i.e.} \quad \frac{\Delta\omega}{\omega_0} \approx \frac{1}{Q}

so a higher-QQ system has a narrower resonance. (Note Δω\Delta\omega is a frequency in s⁻¹, whereas 1/Q1/Q is dimensionless; the two are related through ω0\omega_0.)

Practical resonance. In engineering, resonance is exploited (musical instruments, radio tuning, MRI) or avoided (bridge design, building seismic isolation, turbine-blade fatigue). The 1940 Tacoma Narrows Bridge collapse is a classic—though often oversimplified—example of wind-driven resonance.

Authoritative source for this section: Landau & Lifshitz, Mechanics, §14: Forced Oscillations

Example problem

A resonator has ω0=1000\omega_0 = 1000 rad/s and Q=50Q = 50. Estimate its bandwidth Δω\Delta\omega and the fractional bandwidth Δω/ω0\Delta\omega/\omega_0.

Solution

Δω=ω0/Q=1000/50=20\Delta\omega = \omega_0/Q = 1000/50 = 20 rad/s. The fractional bandwidth is Δω/ω0=1/Q=0.02\Delta\omega/\omega_0 = 1/Q = 0.02 (2%).

Energy and power at resonance

The instantaneous power delivered by the drive is P(t)=F0cos⁡(ωt) x˙P(t) = F_0\cos(\omega t)\,\dot{x}. Its cycle average is:

⟨P⟩=12F0ωXsin⁡δ\langle P \rangle = \frac{1}{2} F_0 \omega X \sin\delta
Details

Maximum power transfer at ω=ω0\omega = \omega_0. The average power is proportional to sin⁡δ\sin\delta. Since δ=π/2\delta = \pi/2 exactly at ω=ω0\omega = \omega_0, the power is maximized there—not at the displacement peak ωres\omega_{res}. Physically, at ω=ω0\omega = \omega_0 the velocity x˙\dot{x} is in phase with the force, so the drive does positive work on every half-cycle. This is why velocity resonance (and power resonance) sits at ω0\omega_0 while displacement resonance sits slightly below it.

At resonance the average power is ⟨P⟩ω0=F02/(2c)\langle P\rangle_{\omega_0} = F_0^2/(2c), and the stored energy is E=12kXmax⁡2E = \tfrac{1}{2}k X_{\max}^2. The quality factor can be read as the number of radians of phase over which the stored energy decays by ee: Q≈ω0E/⟨P⟩Q \approx \omega_0 E / \langle P\rangle.

Authoritative source for this section: OpenStax University Physics Volume 1, §15.6: Forced Oscillations and Resonance

Example problem

Using F0=10F_0 = 10 N and c=2c = 2 N·s/m, find the average power delivered at resonance.

Solution

⟨P⟩ω0=F02/(2c)=100/4=25\langle P\rangle_{\omega_0} = F_0^2/(2c) = 100/4 = 25 W.

Electrical analog: the RLC circuit

The same mathematics governs a series RLC circuit driven by a voltage source:

Lq¨+Rq˙+1Cq=V0cos⁡(ωt)L\ddot{q} + R\dot{q} + \frac{1}{C}q = V_0\cos(\omega t)
Details

The mapping between the mechanical and electrical systems is exact:

MechanicalElectricalRole
mass minductance Linertia / energy in field
damping cresistance Rdissipation
stiffness k1/capacitance 1/Crestoring / energy storage
displacement xcharge qstate variable
force F0voltage V0drive

Under this map ω0=1/LC\omega_0 = 1/\sqrt{LC}, γ=R/(2L)\gamma = R/(2L), and Q=ω0L/R=L/C/RQ = \omega_0 L/R = \sqrt{L/C}/R. The resonance peak becomes the tuned circuit: at ω=ω0\omega = \omega_0 the inductive and capacitive reactances cancel, the impedance is purely resistive (Z=RZ = R), and the current is maximal. This is the basis of radio tuning and band-pass filters, and it is the electrical realization of everything in the resonance analysis above.

Authoritative source for this section: OpenStax University Physics Volume 2, §16.5: RLC Series Circuits

Example problem

A series RLC circuit has L=10L = 10 mH, C=25C = 25 nF, R=10R = 10 Ω. Find ω0\omega_0 and QQ.

Solution

ω0=1/LC=1/10−2×25×10−9=1/2.5×10−10≈6.3×104\omega_0 = 1/\sqrt{LC} = 1/\sqrt{10^{-2} \times 25 \times 10^{-9}} = 1/\sqrt{2.5 \times 10^{-10}} \approx 6.3 \times 10^4 rad/s.

Q=ω0L/R=6.3×104×10−2/10≈63Q = \omega_0 L/R = 6.3 \times 10^4 \times 10^{-2}/10 \approx 63.

Parametric excitation

Instead of applying an external force, a system parameter (e.g., length, stiffness, inductance) is modulated periodically:

mx¨+cx˙+k0[1+εcos⁡(Ωt)]x=0m\ddot{x} + c\dot{x} + k_0\bigl[1 + \varepsilon\cos(\Omega t)\bigr]x = 0
Details

This is a Mathieu equation (when ε\varepsilon is small and c=0c = 0). The key feature: energy is pumped into the system by modulating a parameter, not by applying a direct force.

Parametric resonance occurs when the modulation frequency is near twice the natural frequency, Ω≈2ω0\Omega \approx 2\omega_0. Within the instability band the amplitude grows exponentially, provided the modulation depth ε\varepsilon exceeds a threshold set by damping. The instability region in the (Ω,ε)(\Omega, \varepsilon) plane is a “tongue” opening out from Ω=2ω0\Omega = 2\omega_0:

ε Ω/ω₀ 2 ε_th (damping) instability (parametric resonance)

Examples:

Parametric excitation is a subcase of driven oscillations: the “drive” enters through a coefficient rather than a source term. The distinction matters because parametric resonance occurs at Ω=2ω0\Omega = 2\omega_0 (not ω0\omega_0) and can amplify signals without adding a direct force.

Authoritative source for this section: Landau & Lifshitz, Mechanics, §14: Forced Oscillations

Example problem

A pendulum of natural frequency ω0=2\omega_0 = 2 rad/s is driven parametrically. At what modulation frequency does the primary parametric resonance occur?

Solution

The primary parametric resonance is at Ω=2ω0=4\Omega = 2\omega_0 = 4 rad/s.

Simulation

Quizz

Quizz

Problem: A driven damped oscillator has ω0=10\omega_0 = 10 rad/s, γ=0.5\gamma = 0.5 s⁻¹, and F0/m=2F_0/m = 2 m/s². Find the steady-state amplitude at (a) ω=ω0\omega = \omega_0, (b) ω=2ω0\omega = 2\omega_0, and (c) the ratio X(ω0)/X(2ω0)X(\omega_0)/X(2\omega_0).

Solution

(a) At resonance: X=(F0/m)/(2γ ω0)=2/(2×0.5×10)=2/10=0.2X = (F_0/m)/(2\gamma\,\omega_0) = 2/(2 \times 0.5 \times 10) = 2/10 = 0.2 m.

(b) At ω=20\omega = 20 rad/s:

X=2/(100−400)2+(2×0.5×20)2=2/90000+400=2/90400≈2/300.7≈0.00665X = 2/\sqrt{(100-400)^2 + (2 \times 0.5 \times 20)^2} = 2/\sqrt{90000 + 400} = 2/\sqrt{90400} \approx 2/300.7 \approx 0.00665 m.

(c) Ratio: X(ω0)/X(2ω0)=0.2/0.00665≈30X(\omega_0)/X(2\omega_0) = 0.2/0.00665 \approx 30. The amplitude at 2ω02\omega_0 is about 30× smaller than at resonance, illustrating the sharpness of the peak for Q=ω0/(2γ)=10Q = \omega_0/(2\gamma) = 10.