Physics · Book 3 · Bachelor Year 1

University Physics — Year 1

University Physics — Year 1 · Bachelor Year 1

14Mechanical Oscillators: Damping and Resonance

Push a child’s swing at the wrong moment and nothing much happens; push it once per period, gently, and within a dozen swings the child is shrieking at the top of the arc. A washing machine shudders violently at one particular speed of its spin cycle and is quiet above it. A wine glass sings when a wet finger circles its rim at just the right rate, and shatters when a singer holds the same note loudly enough. Every one of these is a resonance: a system with a natural frequency, driven near that frequency, accumulates energy cycle after cycle until friction bleeds it off as fast as it comes in. This chapter takes the mechanical oscillator — free, damped, then driven — and derives the resonance curves whose electrical twins appeared in Chapter 8.

14.1 The harmonic oscillator

Definition 14.1 (Harmonic oscillator)

A system whose displacement xx from a stable equilibrium obeys

x¨+ω02x=0\ddot x + \omega_0^2\,x = 0

is a harmonic oscillator of natural angular frequency ω0\omega_0 (period T0=2π/ω0T_0 = 2\pi/\omega_0). The solution x=Acos(ω0t+φ)x = A\cos(\omega_0t + \varphi) has an amplitude and phase fixed by the initial conditions and a period independent of the amplitude (isochronism).

Proposition 14.2 (Where harmonic oscillators come from)

  • Mass mm on a spring kk: ω02=k/m\omega_0^2 = k/m (Hooke’s law).
  • Simple pendulum of length \ell, small angles: ω02=g/\omega_0^2 = g/\ell.
  • Any particle near a stable equilibrium of a potential EpE_p: ω02=Ep(x0)/m\omega_0^2 = E_p''(x_0)/m (Proposition 13.12).
  • LC circuit: ω02=1/LC\omega_0^2 = 1/LC, with qxq \leftrightarrow x, LmL \leftrightarrow m, 1/Ck1/C \leftrightarrow k.

The energy E=12mx˙2+12mω02x2=12mω02A2E = \tfrac12 m\dot x^2 + \tfrac12 m\omega_0^2x^2 = \tfrac12 m\omega_0^2A^2 is constant, sloshing between kinetic and potential form twice per period.

Proof. Newton’s second law with F=kxF = -kx; the pendulum equation θ¨+(g/)sinθ=0\ddot\theta + (g/\ell)\sin\theta = 0 with sinθθ\sin\theta \approx \theta; Taylor’s formula at the minimum of EpE_p; the loop law of the LC circuit. Energy: differentiate and use the equation.

Example 14.3 (A bobbing cylinder)

A vertical cylinder of cross-section SS floats with a depth hh immersed (mg=ρwgShmg = \rho_wgSh, buoyancy from Chapter 21). Pushed down by xx, it receives an extra upward force ρwgSx\rho_wgSx: mx¨=ρwgSxm\ddot x = -\rho_wgSx, so ω02=ρwgS/m=g/h\omega_0^2 = \rho_wgS/m = g/h — the same as a pendulum of length hh. A buoy immersed by 1m1\,\mathrm{m} bobs with a 2s2\,\mathrm{s} period, whatever its mass.

14.2 The damped oscillator

Proposition 14.4 (Damped oscillator)

With a viscous friction αx˙-\alpha\dot x, the equation becomes

x¨+ω0Qx˙+ω02x=0,Q=kmα=mω0α,\ddot x + \frac{\omega_0}{Q}\dot x + \omega_0^2x = 0, \qquad Q = \frac{\sqrt{km}}{\alpha} = \frac{m\omega_0}{\alpha},

the canonical second-order form of Chapter 7: for Q>12Q > \tfrac12 a pseudo-periodic decay x=Aet/τcos(ωt+φ)x = A\eu^{-t/\tau}\cos(\omega t + \varphi) with τ=2Q/ω0\tau = 2Q/\omega_0, ω=ω011/4Q2\omega = \omega_0\sqrt{1 - 1/4Q^2}; the energy decays as e2t/τ\eu^{-2t/\tau} and

Q=2πenergy storedenergy lost per period(Q1).Q = 2\pi\,\frac{\text{energy stored}}{\text{energy lost per period}} \quad (Q \gg 1) .

Proof. Divide mx¨+αx˙+kx=0m\ddot x + \alpha\dot x + kx = 0 by mm. For the energy: EA2e2t/τE \propto A^2\eu^{-2t/\tau} loses the fraction 2T/τ=2π/Q2T/\tau = 2\pi/Q per period when that is small.

Example 14.5 (Ringing times)

A tuning fork at 440Hz440\,\mathrm{Hz} still sounds ten seconds after being struck: its amplitude fell to about 5%5\% in 3τ=10s3\tau = 10\,\mathrm{s}, so τ=3.3s\tau = 3.3\,\mathrm{s} and Q=ω0τ/24600Q = \omega_0\tau/2 \approx 4600. A car’s suspension (Q0.7Q \approx 0.7) stops within one period; a child’s swing (Q20Q \approx 20) loses a third of its energy per swing and needs a push every time.

14.3 Forced oscillations and resonance

Theorem 14.6 (Sinusoidally driven oscillator)

Driven by a force F0cosωtF_0\cos\omega t, the damped oscillator

mx¨+αx˙+kx=F0cosωtm\ddot x + \alpha\dot x + kx = F_0\cos\omega t

settles, after a transient of duration τ\sim\tau, into the forced regime x=Xcos(ωt+φ)x = X\cos(\omega t + \varphi) with, writing xr=ω/ω0x_r = \omega/\omega_0 and X0=F0/kX_0 = F_0/k the static deflection,

X=X0(1xr2)2+xr2/Q2,tanφ=xr/Q1xr2,φ[π,0].X = \frac{X_0}{\sqrt{(1 - x_r^2)^2 + x_r^2/Q^2}}, \qquad \tan\varphi = -\frac{x_r/Q}{1 - x_r^2}, \quad \varphi \in [-\pi, 0] .

The amplitude XX peaks, for Q>1/2Q > 1/\sqrt2, at ωr=ω011/2Q2\omega_r = \omega_0\sqrt{1 - 1/2Q^2} with XmaxQX0X_{\max} \approx QX_0; the velocity amplitude V=ωXV = \omega X peaks exactly at ω0\omega_0, where Vmax=F0/αV_{\max} = F_0/\alpha, and its resonance curve

VVmax=11+Q2(xr1/xr)2\frac{V}{V_{\max}} = \frac{1}{\sqrt{1 + Q^2(x_r - 1/x_r)^2}}

has the bandwidth Δω=ω0/Q\Delta\omega = \omega_0/Q. At ω0\omega_0 the displacement lags the force by π/2\pi/2 and the velocity is in phase with it.

Proof. Complex method (Chapter 8): X(ω2+jωω0/Q+ω02)=F0/m\underline X(-\omega^2 + j\omega\omega_0/Q + \omega_0^2) = F_0/m, so X=(F0/m)/(ω02ω2+jωω0/Q)\underline X = (F_0/m)/(\omega_0^2 - \omega^2 + j\omega\omega_0/Q); divide numerator and denominator by ω02=k/m\omega_0^2 = k/m. Modulus and argument follow; the maximum of XX was located in Proposition 8.10 (same function). Velocity: V=jωX\underline V = j\omega\underline X gives V=ωXV = \omega X; writing (1xr2)/xr=1/xrxr(1 - x_r^2)/x_r = 1/x_r - x_r turns the denominator into 1/Q2+(1/xrxr)2\sqrt{1/Q^2 + (1/x_r - x_r)^2}, i.e. the RLC form. At xr=1x_r = 1, X=jQX0\underline X = -jQX_0: a lag of π/2\pi/2.

Amplitude (left) and phase (right) of the forced oscillator against the driving frequency, for several quality factors. A sharp peak near _0 for large Q, none below Q = 1/√2; the displacement goes from in phase (low ) to opposite (high ), through a π/2 lag at _0 — the sharper the resonance, the more abrupt the switch. Amplitude (left) and phase (right) of the forced oscillator against the driving frequency, for several quality factors. A sharp peak near _0 for large Q, none below Q = 1/√2; the displacement goes from in phase (low ) to opposite (high ), through a π/2 lag at _0 — the sharper the resonance, the more abrupt the switch.
Amplitude (left) and phase (right) of the forced oscillator against the driving frequency, for several quality factors. A sharp peak near ω0\omega_0 for large QQ, none below Q=1/2Q = 1/\sqrt2; the displacement goes from in phase (low ω\omega) to opposite (high ω\omega), through a π/2\pi/2 lag at ω0\omega_0 — the sharper the resonance, the more abrupt the switch.

Remark 14.7 (Two resonances, one circuit)

The correspondence xqx \leftrightarrow q, x˙i\dot x \leftrightarrow i, mLm \leftrightarrow L, αR\alpha \leftrightarrow R, k1/Ck \leftrightarrow 1/C, FeF \leftrightarrow e maps the driven oscillator onto the series RLC of Chapter 8: the velocity resonance is the current resonance (always at ω0\omega_0, bandwidth ω0/Q\omega_0/Q), the amplitude resonance is the charge, i.e. capacitor-voltage, resonance (slightly below ω0\omega_0, absent for low QQ). Power absorbed, Fx˙=12F0Vcosφv\langle F\dot x\rangle = \tfrac12 F_0V\cos\varphi_v, peaks at ω0\omega_0 with the value F02/2αF_0^2/2\alpha, and the half-power points are the velocity bandwidth — exactly as for the RLC.

Example 14.8 (Resonant speeds)

A car’s suspension has f01.3Hzf_0 \approx 1.3\,\mathrm{Hz} (Problem 7.1). On a road with corrugations every λ=10m\lambda = 10\,\mathrm{m} the excitation frequency is v/λv/\lambda: resonance at v=λf0=13m/sv = \lambda f_0 = 13\,\mathrm{m}/\mathrm{s}, about 47km/h47\,\mathrm{km}/\mathrm{h} — the speed at which a washboard track shakes the car hardest, and the reason a damper with Q0.7Q \approx 0.7 is fitted: with Q=5Q = 5 the body would bounce with five times the road’s amplitude. Above resonance (xr1x_r \gg 1) the amplitude falls as X0/xr2X_0/x_r^2: the body “floats” over fast bumps — isolation, the other face of resonance.

14.4 Base excitation: isolation and the seismometer

Proposition 14.9 (Oscillator shaken by its support)

Let the support of a mass–spring–damper move by xg(t)x_g(t) (the road under a car, the ground under a seismometer), and let uu be the displacement of the mass relative to the support and x=xg+ux = x_g + u its absolute displacement (from equilibrium). Then

mu¨+αu˙+ku=mx¨g,m\ddot u + \alpha\dot u + ku = -m\ddot x_g ,

and for xg=Xgcosωtx_g = X_g\cos\omega t, with xr=ω/ω0x_r = \omega/\omega_0:

UXg=xr2(1xr2)2+xr2/Q2,XXg=1+xr2/Q2(1xr2)2+xr2/Q2.\frac{U}{X_g} = \frac{x_r^2}{\sqrt{(1 - x_r^2)^2 + x_r^2/Q^2}}, \qquad \frac{X}{X_g} = \sqrt{\frac{1 + x_r^2/Q^2}{(1 - x_r^2)^2 + x_r^2/Q^2}} .

For xr1x_r \gg 1: UXgU \to X_g (the mass stays still in space — a seismometer reads the ground’s displacement) and X0X \to 0 (the mass is isolated from the shaking); for xr1x_r \ll 1: Ux¨g/ω02U \approx \ddot x_g/\omega_0^2 (an accelerometer) and XXgX \approx X_g (the mass follows the support).

Proof. The spring and damper act on the relative motion: mx¨=kuαu˙m\ddot x = -ku - \alpha\dot u, and x¨=x¨g+u¨\ddot x = \ddot x_g + \ddot u. In complex form U(ω02ω2+jωω0/Q)=ω2Xg\underline U(\omega_0^2 - \omega^2 + j\omega\omega_0/Q) = \omega^2X_g, whence U/XgU/X_g; then X=U+Xg=Xg(1+jxr/Q)/(1xr2+jxr/Q)\underline X = \underline U + X_g = X_g(1 + jx_r/Q)/(1 - x_r^2 + jx_r/Q).

A mass shaken through its support. Left: relative motion U/X_g — an accelerometer below _0 (slope 2 on log scales), a seismometer above. Right: absolute motion X/X_g — the mass follows the support at low frequency, is isolated from it well above _0; a high Q buys nothing but a dangerous peak. A mass shaken through its support. Left: relative motion U/X_g — an accelerometer below _0 (slope 2 on log scales), a seismometer above. Right: absolute motion X/X_g — the mass follows the support at low frequency, is isolated from it well above _0; a high Q buys nothing but a dangerous peak.
A mass shaken through its support. Left: relative motion U/XgU/X_g — an accelerometer below ω0\omega_0 (slope 22 on log scales), a seismometer above. Right: absolute motion X/XgX/X_g — the mass follows the support at low frequency, is isolated from it well above ω0\omega_0; a high QQ buys nothing but a dangerous peak.

Example 14.10 (Why machines sit on soft mounts)

A 50Hz50\,\mathrm{Hz} motor on rubber mounts tuned to f0=10Hzf_0 = 10\,\mathrm{Hz} (xr=5x_r = 5, Q2Q \approx 2): X/Xg=(1+6.25)/(576+6.25)=0.11X/X_g = \sqrt{(1 + 6.25)/(576 + 6.25)} = 0.11 — nine tenths of the vibration stays out of the floor. The mounts must be softer than intuition suggests: stiffening them raises f0f_0 toward 5050 Hz and destroys the isolation. The same formula, read the other way, is the car of Example 14.8 on a washboard road.

14.5 Energy: how a resonance builds up

Proposition 14.11 (Energy balance in the forced regime)

In the steady forced regime the average power supplied by the driving force equals the average power dissipated by friction, 12αV2\tfrac12\alpha V^2. At resonance the stored energy 12mVmax2\tfrac12 mV_{\max}^2 is Q/2πQ/2\pi times the energy supplied per period: the oscillator accumulates the drive’s work over about QQ cycles before losses match the input, which is also the number of cycles (τ/T\sim\tau/T) the transient lasts.

Proof. Average the energy theorem over a period:  ⁣dE/ ⁣dt=0\langle\dd E/\dd t\rangle = 0 in steady state, so Fx˙=αx˙2=12αV2\langle F\dot x\rangle = \langle\alpha\dot x^2\rangle = \tfrac12\alpha V^2. At resonance Vmax=F0/αV_{\max} = F_0/\alpha, stored energy 12m(F0/α)2\tfrac12 m(F_0/\alpha)^2, supplied per period 12α(F0/α)2T\tfrac12\alpha(F_0/\alpha)^2T; ratio m/(αT)=Q/2πm/(\alpha T) = Q/2\pi.

Example 14.12 (The swing)

A swing of length 2.5m2.5\,\mathrm{m} (T0=3.2sT_0 = 3.2\,\mathrm{s}, Q20Q \approx 20) gets a push of 50N50\,\mathrm{N} over 0.5m0.5\,\mathrm{m} each period: 25J25\,\mathrm{J} per cycle. It grows until the loss per period, 2πE/Q2\pi E/Q, equals 25J25\,\mathrm{J}: E=80JE = 80\,\mathrm{J}; with a 30kg30\,\mathrm{kg} child, mg(1cosθmax)=80Jmg\ell(1 - \cos\theta_{\max}) = 80\,\mathrm{J} gives θmax27\theta_{\max} \approx 27^\circ. Pushing at any other rhythm, the work of the pushes would alternate in sign and average to nothing.

Start-up of a resonance (Q = 5, drive switched on at t = 0 at _0): the amplitude grows toward QX_0 with the time constant = 2Q/ _0, i.e. over about Q/π periods — the transient of  superposed on the forced regime.
Start-up of a resonance (Q=5Q = 5, drive switched on at t=0t = 0 at ω0\omega_0): the amplitude grows toward QX0QX_0 with the time constant τ=2Q/ω0\tau = 2Q/\omega_0, i.e. over about Q/πQ/\pi periods — the transient of Chapter 7 superposed on the forced regime.

14.6 Exercises

Exercise 14.1

A 0.20kg0.20\,\mathrm{kg} mass on a spring of stiffness 80N/m80\,\mathrm{N}/\mathrm{m} oscillates with amplitude 5.0cm5.0\,\mathrm{cm}. Angular frequency, frequency, period, maximal speed, total energy.

Solution

Solution of Exercise 14.1.

ω0=80/0.20=20rad/s\omega_0 = \sqrt{80/0.20} = 20\,\mathrm{rad}/\mathrm{s}, f0=3.2Hzf_0 = 3.2\,\mathrm{Hz}, T0=0.31sT_0 = 0.31\,\mathrm{s}; vmax=Aω0=1.0m/sv_{\max} = A\omega_0 = 1.0\,\mathrm{m}/\mathrm{s}; E=12kA2=0.10JE = \tfrac12 kA^2 = 0.10\,\mathrm{J}.

Exercise 14.2

Period of a 1.0m1.0\,\mathrm{m} pendulum on Earth and on the Moon (g=1.62m/s2g = 1.62\,\mathrm{m}/\mathrm{s}^{2}). What happens to the mass–spring period on the Moon?

Solution

Solution of Exercise 14.2.

T=2π/gT = 2\pi\sqrt{\ell/g}: 2.0s2.0\,\mathrm{s} on Earth, 4.9s4.9\,\mathrm{s} on the Moon. The mass–spring period, 2πm/k2\pi\sqrt{m/k}, does not involve gg: unchanged.

Exercise 14.3

An oscillator has f0=2.0Hzf_0 = 2.0\,\mathrm{Hz} and Q=50Q = 50. Decay time τ\tau of the amplitude; time for the energy to halve; number of oscillations before the amplitude is down to 5%5\%.

Solution

Solution of Exercise 14.3.

τ=2Q/ω0=100/(4π)=8.0s\tau = 2Q/\omega_0 = 100/(4\pi) = 8.0\,\mathrm{s}; energy e2t/τ\propto \eu^{-2t/\tau} halves at t=τln2/2=2.8st = \tau\ln2/2 = 2.8\,\mathrm{s}; 5%5\% after 3τ=24s3\tau = 24\,\mathrm{s}, about 48Q48 \approx Q oscillations.

Exercise 14.4

A mass on a spring (k=100N/mk = 100\,\mathrm{N}/\mathrm{m}, Q=10Q = 10) is driven by a force of amplitude 1.0N1.0\,\mathrm{N}. Static deflection; amplitude at resonance; amplitude at ω=3ω0\omega = 3\omega_0.

Solution

Solution of Exercise 14.4.

X0=F0/k=1.0cmX_0 = F_0/k = 1.0\,\mathrm{cm}; at resonance QX0=10cm\approx QX_0 = 10\,\mathrm{cm}; at 3ω03\omega_0: X0/64+0.09=0.13cmX_0/\sqrt{64 + 0.09} = 0.13\,\mathrm{cm}.

Exercise 14.5 ★★

A car (f0=1.3Hzf_0 = 1.3\,\mathrm{Hz}, Q=0.7Q = 0.7) crosses speed bumps 8.0m8.0\,\mathrm{m} apart. Resonant speed? At 60km/h60\,\mathrm{km}/\mathrm{h}, what fraction of the bump amplitude does the body feel? Comment.

Solution

Solution of Exercise 14.5.

v=Lf0=8.0×1.3=10m/sv = Lf_0 = 8.0 \times 1.3 = 10\,\mathrm{m}/\mathrm{s} (37km/h37\,\mathrm{km}/\mathrm{h}). At 60km/h60\,\mathrm{km}/\mathrm{h}: f=16.7/8.0=2.1Hzf = 16.7/8.0 = 2.1\,\mathrm{Hz}, xr=1.6x_r = 1.6: X/Xg=(1+5.2)/(2.4+5.2)=0.90X/X_g = \sqrt{(1 + 5.2)/(2.4 + 5.2)} = 0.90 — barely isolated; isolation needs xr1x_r \gg 1 (at 120km/h120\,\mathrm{km}/\mathrm{h}, xr=3.2x_r = 3.2: 0.450.45).

Exercise 14.6 ★★

Explain with the energy balance why a swing must be pushed at its own period; for Q=20Q = 20 and a push adding 25J25\,\mathrm{J} per cycle, find the steady-state energy, and the amplitude for a 30kg30\,\mathrm{kg} child on a 2.5m2.5\,\mathrm{m} swing.

Solution

Solution of Exercise 14.6.

The work of a push is positive only when it acts along the motion; pushed every period, all pushes add; at any other rhythm they alternate. Steady state: 2πE/Q=25J2\pi E/Q = 25\,\mathrm{J}, E=80JE = 80\,\mathrm{J}; mg(1cosθmax)=80mg\ell(1 - \cos\theta_{\max}) = 80: 1cosθmax=0.1091 - \cos\theta_{\max} = 0.109, θmax=27\theta_{\max} = 27^\circ.

Exercise 14.7 ★★

Successive maxima of a free oscillation decrease by 10%10\% each period. Find QQ and the bandwidth (in hertz) of its velocity resonance if f0=2.0Hzf_0 = 2.0\,\mathrm{Hz}.

Solution

Solution of Exercise 14.7.

δ=ln(1/0.9)=0.105\delta = \ln(1/0.9) = 0.105, Q=π/δ=30Q = \pi/\delta = 30; Δf=f0/Q=0.067Hz\Delta f = f_0/Q = 0.067\,\mathrm{Hz}.

Exercise 14.8 ★★

Give the phase of the displacement relative to the force well below resonance, at resonance and well above, and the phase of the velocity at resonance. Why is the velocity resonance the one that matters for power?

Solution

Solution of Exercise 14.8.

Displacement: in phase (φ=0\varphi = 0) well below, π/2-\pi/2 at resonance, π-\pi (opposite) well above. Velocity at resonance: in phase with the force. Power is Fv\vect F\cdot\vect v: it is the velocity’s phase, not the displacement’s, that decides how much work the drive does.

Exercise 14.9 ★★

Show that the velocity amplitude V=ωXV = \omega X is maximal exactly at ω0\omega_0, with Vmax=F0/αV_{\max} = F_0/\alpha, and that V/VmaxV/V_{\max} has the form 1/1+Q2(xr1/xr)21/\sqrt{1 + Q^2(x_r - 1/x_r)^2}.

Solution

Solution of Exercise 14.9.

V=ωX=ω(F0/m)/(ω02ω2)2+(ωω0/Q)2V = \omega X = \omega(F_0/m)/\sqrt{(\omega_0^2 - \omega^2)^2 + (\omega\omega_0/Q)^2}; divide numerator and denominator by ωω0/Q\omega\omega_0/Q: V=(F0Q/mω0)/1+Q2(ω0/ωω/ω0)2V = (F_0Q/m\omega_0)/ \sqrt{1 + Q^2(\omega_0/\omega - \omega/\omega_0)^2}, and F0Q/mω0=F0/αF_0Q/m\omega_0 = F_0/\alpha. The square root is 1\geq 1, equal to 11 exactly at ω=ω0\omega = \omega_0.

Exercise 14.10 ★★★

Compute the average power absorbed at resonance for F0=1.0NF_0 = 1.0\,\mathrm{N}, α=0.050Ns/m\alpha = 0.050\,\mathrm{N}\,\mathrm{s}/\mathrm{m}, and the stored energy if m=0.20kgm = 0.20\,\mathrm{kg}; check that their ratio is Q/ω0Q/\omega_0 with ω0=20rad/s\omega_0 = 20\,\mathrm{rad}/\mathrm{s}.

Solution

Solution of Exercise 14.10.

P=F02/2α=10WP = F_0^2/2\alpha = 10\,\mathrm{W}; Vmax=F0/α=20m/sV_{\max} = F_0/\alpha = 20\,\mathrm{m}/\mathrm{s}, E=12mVmax2=40JE = \tfrac12 mV_{\max}^2 = 40\,\mathrm{J}; E/P=4.0sE/P = 4.0\,\mathrm{s} and Q/ω0=m/α=4.0sQ/\omega_0 = m/\alpha = 4.0\,\mathrm{s}.

Exercise 14.11 ★★★

A tuning fork at 440Hz440\,\mathrm{Hz} is heard for 10s10\,\mathrm{s} (amplitude down to 5%5\%). Estimate τ\tau and QQ, the bandwidth of its resonance, and comment on why it is a good frequency standard.

Solution

Solution of Exercise 14.11.

3τ=10s3\tau = 10\,\mathrm{s}: τ=3.3s\tau = 3.3\,\mathrm{s}, Q=ω0τ/2=2π×440×3.3/24600Q = \omega_0\tau/2 = 2\pi \times 440 \times 3.3/2 \approx 4600; Δf=f0/Q=0.1Hz\Delta f = f_0/Q = 0.1\,\mathrm{Hz}. It responds to, and emits, essentially one frequency: a standard good to a few parts in 10410^4.

Exercise 14.12 ★★★

Derive the bobbing frequency of a floating cylinder (Example 14.3) from Newton’s law and Archimedes’ force ρwgS(h+x)\rho_wgS(h + x); show ω02=g/h\omega_0^2 = g/h; compute the period for h=10cmh = 10\,\mathrm{cm}. Why is the result independent of the mass and of the liquid?

Solution

Solution of Exercise 14.12.

mx¨=mgρwgS(h+x)=ρwgSxm\ddot x = mg - \rho_wgS(h + x) = -\rho_wgSx (since mg=ρwgShmg = \rho_wgSh): ω02=ρwgS/m=g/h\omega_0^2 = \rho_wgS/m = g/h; T=2π0.10/9.81=0.63sT = 2\pi\sqrt{0.10/9.81} = 0.63\,\mathrm{s}. Both the mass and the liquid’s density enter only through the equilibrium depth hh.

A seismograph: a heavy mass on a soft suspension stays still while the ground moves, and the pen writes the difference — the base-excited oscillator of the weekend problem.
A seismograph: a heavy mass on a soft suspension stays still while the ground moves, and the pen writes the difference — the base-excited oscillator of the weekend problem.

14.7 Problem: The seismograph

Problem 14.1

Weekend problem — a mass hangs from a spring inside a box bolted to the ground; the ground shakes: what does the pen record, why the same instrument measures displacement at high frequency and acceleration at low frequency, and why a long-period seismometer cannot simply be a long spring

A mass m=1.0kgm = 1.0\,\mathrm{kg} hangs from a spring of stiffness kk inside a rigid frame fixed to the ground; a damper exerts αu˙-\alpha\dot u, where uu is the displacement of the mass relative to the frame, measured from its equilibrium position. The ground, hence the frame, moves vertically by xg(t)x_g(t). The frame of the distant stars is inertial.

Part I — The equation.

  1. Write the absolute position of the mass as x=xg+u+constx = x_g + u + \text{const}, list the forces on the mass (in the inertial frame), and show that

    mu¨+αu˙+ku=mx¨g.m\ddot u + \alpha\dot u + ku = -m\ddot x_g .
  2. Put it in canonical form; identify ω0\omega_0 and QQ.
  3. The instrument is built with f0=0.50Hzf_0 = 0.50\,\mathrm{Hz} and Q=0.70Q = 0.70: compute kk and α\alpha.
  4. A constant ground acceleration aga_g (the frame tilting, say) gives what steady uu? Interpret.
  5. What does the pen attached to the mass, writing on a drum fixed to the frame, record?

Part II — Harmonic ground motion. The ground moves as xg=Xgcosωtx_g = X_g\cos\omega t; write xr=ω/ω0x_r = \omega/\omega_0.

  1. Using complex amplitudes, show that U=Xgxr2/(1xr2+jxr/Q)\underline U = X_g\,x_r^2/ (1 - x_r^2 + jx_r/Q).
  2. Deduce U/XgU/X_g as a function of xrx_r and QQ.
  3. Limit xr1x_r \gg 1: show UXgU \to X_g and explain physically what the mass does (what is the instrument measuring?).
  4. Limit xr1x_r \ll 1: show Uxr2Xg=ag/ω02U \approx x_r^2X_g = a_g/\omega_0^2 with aga_g the ground’s acceleration amplitude: what is the instrument now?
  5. Value at xr=1x_r = 1; why is Q0.7Q \approx 0.7 a sensible choice?
  6. Sketch U/XgU/X_g against xrx_r (log scales) for Q=0.7Q = 0.7 and Q=5Q = 5.
  7. Phase of uu relative to xgx_g at high frequency; interpret.

Part III — Designing for earthquakes. Distant earthquakes produce surface waves at 0.050.05\, to 0.1Hz0.1\,\mathrm{Hz} with millimeter amplitudes; local tremors reach 55\, to 10Hz10\,\mathrm{Hz}.

  1. To record 0.05Hz0.05\,\mathrm{Hz} waves as displacements one wants f00.05Hzf_0 \ll 0.05\,\mathrm{Hz}, say 0.02Hz0.02\,\mathrm{Hz}. What static stretch would the spring have under the mass’s weight? Comment.
  2. A simple pendulum of the same frequency: what length?
  3. Real long-period instruments use a nearly horizontal “garden-gate” pendulum, or electronic force feedback. Explain in one sentence what each achieves.
  4. For the 0.50Hz0.50\,\mathrm{Hz} instrument, compute the static stretch: feasible?
  5. With it, a 0.10Hz0.10\,\mathrm{Hz} wave of amplitude 1.0mm1.0\,\mathrm{mm} gives what UU? In which regime is the instrument?
  6. A 5.0Hz5.0\,\mathrm{Hz} local tremor of amplitude 1.0mm1.0\,\mathrm{mm}: UU? Regime?

Part IV — The accelerometer in a phone. A micro-machined mass on silicon springs has f0=1.0kHzf_0 = 1.0\,\mathrm{kHz} and Q=0.7Q = 0.7.

  1. A 10Hz10\,\mathrm{Hz} shake of acceleration amplitude gg: compute UU. How is such a displacement read?
  2. Up to what frequency is the reading within 10%10\% of ag/ω02a_g/\omega_0^2?
  3. What time constant τ=2Q/ω0\tau = 2Q/\omega_0 governs its response to a jolt, and why does that matter for detecting a fall?
  4. The seismometer’s τ\tau: value and consequence after a sudden ground step.
  5. Why is critical (or near-critical) damping chosen for both instruments, rather than a high QQ?
  6. The phone is shaken at exactly 1.0kHz1.0\,\mathrm{kHz}: by what factor is the reading off?
  7. Summarize: one device, two regimes, decided by which single comparison.
Solution

Solution of Problem 14.1.

1. Forces: weight (balanced by the static stretch), spring ku-ku, damper αu˙-\alpha\dot u; mx¨=kuαu˙m\ddot x = -ku - \alpha\dot u with x¨=x¨g+u¨\ddot x = \ddot x_g + \ddot u.

2. u¨+(ω0/Q)u˙+ω02u=x¨g\ddot u + (\omega_0/Q)\dot u + \omega_0^2u = -\ddot x_g, ω0=k/m\omega_0 = \sqrt{k/m}, Q=km/αQ = \sqrt{km}/\alpha.

3. k=mω02=π2=9.9N/mk = m\omega_0^2 = \pi^2 = 9.9\,\mathrm{N}/\mathrm{m}; α=mω0/Q=3.14/0.70=4.5Ns/m\alpha = m\omega_0/Q = 3.14/0.70 = 4.5\,\mathrm{N}\,\mathrm{s}/\mathrm{m}.

4. u=ag/ω02u = -a_g/\omega_0^2: the mass sags by a fixed amount — the instrument is an accelerometer at zero frequency (and a tilt meter).

5. The relative displacement u(t)u(t).

6. ω2U+jω(ω0/Q)U+ω02U=ω2Xg-\omega^2\underline U + j\omega(\omega_0/Q)\underline U + \omega_0^2\underline U = \omega^2X_g; divide by ω02\omega_0^2.

7. U/Xg=xr2/(1xr2)2+xr2/Q2U/X_g = x_r^2/\sqrt{(1 - x_r^2)^2 + x_r^2/Q^2}.

8. UXgU \to X_g: the mass does not move in the inertial frame (its inertia keeps it still while the spring is too soft to drag it); the frame moves around it, and the pen records the ground displacement — a seismometer.

9. Uxr2Xg=ω2Xg/ω02=ag/ω02U \approx x_r^2X_g = \omega^2X_g/\omega_0^2 = a_g/\omega_0^2: an accelerometer.

10. U=QXg=0.7XgU = QX_g = 0.7X_g: no peak; the two regimes join smoothly and the instrument rings not.

11. For Q=0.7Q = 0.7: a line of slope 22 rising to 11 near xr=1x_r = 1, then flat at 11; for Q=5Q = 5: the same with a peak of height 55 at xr=1x_r = 1.

12. UXg\underline U \to -X_g: u=xgu = -x_g, the mass still while the frame moves.

13. Δ=g/ω02=9.81/(2π×0.02)2=620m\Delta\ell = g/\omega_0^2 = 9.81/(2\pi \times 0.02)^2 = 620\,\mathrm{m}: absurd.

14. =g/ω02=620m\ell = g/\omega_0^2 = 620\,\mathrm{m}.

15. Garden gate: the pendulum swings about a nearly vertical axis, so only a tiny fraction of gg restores it and a short arm gives a long period. Force feedback: a coil holds the mass fixed, and the current needed is proportional to the ground acceleration — the stiffness is electronic and adjustable.

16. Δ=9.81/π2=0.99m\Delta\ell = 9.81/\pi^2 = 0.99\,\mathrm{m}: feasible.

17. xr=0.2x_r = 0.2: U/Xg=0.04/0.92+0.08=0.040U/X_g = 0.04/\sqrt{0.92 + 0.08} = 0.040: 40µm40\,\text{µ}\mathrm{m} — accelerometer regime, needs amplification (optical lever or electronics).

18. xr=10x_r = 10: 100/9801+2041.00100/\sqrt{9801 + 204} \approx 1.00: 1.0mm1.0\,\mathrm{mm}, faithful displacement record.

19. ω0=6.3×103rad/s\omega_0 = 6.3 \times 10^{3}\,\mathrm{rad}/\mathrm{s}: U=g/ω02=2.5×107mU = g/\omega_0^2 = 2.5 \times 10^{-7}\,\mathrm{m} — read as the change of a capacitor gap.

20. U/(ag/ω02)=1/(1xr2)2+xr2/Q2U/(a_g/\omega_0^2) = 1/\sqrt{(1 - x_r^2)^2 + x_r^2/Q^2}; with Q=0.7Q = 0.7 it decreases monotonically and stays above 0.90.9 up to xr0.7x_r \approx 0.7: about 700Hz700\,\mathrm{Hz}.

21. τ=2Q/ω0=0.22ms\tau = 2Q/\omega_0 = 0.22\,\mathrm{ms}: it tracks a fall (tens of milliseconds) without lag.

22. τ=1.4/3.14=0.45s\tau = 1.4/3.14 = 0.45\,\mathrm{s}: after a ground step the record settles in half a second with no ringing.

23. A high QQ would make the instrument ring at f0f_0 after every jolt (recording itself rather than the ground) and peak near f0f_0; Q0.7Q \approx 0.7 is flat and quick.

24. U=QXg=0.7XgU = QX_g = 0.7X_g against the ideal XgX_g (or ag/ω02=Xga_g/\omega_0^2 = X_g): 30%30\% low.

25. Compare the signal frequency with f0f_0: well above, the mass is a fixed reference and uu is the displacement; well below, uu is the acceleration over ω02\omega_0^2.

Terms defined in this chapter

See all 393 terms in the glossary