The amplitude A = F₀/(m·(ω₀² − ω²)) of the previous lesson grows as ω approaches ω₀, and at ω = ω₀ it has a zero in the denominator. That is not a large answer. It is no answer: the form A·cos(ωt) does not solve the equation at all when the push has the natural frequency.
ω = ω₀: m·(−ω₀²·A) + k·A = (k − m·ω₀²)·A = 0
0 = F₀ is impossible for F₀ ≠ 0
The cure is the one used for a repeated root: multiply by t. Try x_p = B·t·sin(ω₀t) and differentiate with the product rule.
x_p = B·t·sin(ω₀t)
x_p′ = B·sin(ω₀t) + B·ω₀·t·cos(ω₀t)
x_p″ = 2B·ω₀·cos(ω₀t) − B·ω₀²·t·sin(ω₀t)
m·x_p″ + k·x_p = 2m·B·ω₀·cos(ω₀t) + B·(k − m·ω₀²)·t·sin(ω₀t)
Theorem · The resonant solutionFor ω = ω₀, the function x_p = (F₀/(2m·ω₀))·t·sin(ω₀t) solves m·x″ + k·x = F₀·cos(ω₀t).
Since k = m·ω₀², the last term of the display above is B·(k − m·ω₀²)·t·sin(ω₀t) = 0, so m·x_p″ + k·x_p = 2m·B·ω₀·cos(ω₀t).
Setting 2m·B·ω₀ = F₀ gives B = F₀/(2m·ω₀), and then the left side is F₀·cos(ω₀t).
□
Example · Resonance with m = 1, k = 4, F₀ = 3Here ω₀ = 2, and the push is 3·cos(2t). Then B = 3/(2·1·2) = 3/4 and x_p = (3/4)·t·sin(2t). The check: x_p″ = 3·cos(2t) − 3t·sin(2t), so x_p″ + 4x_p = 3·cos(2t), the push exactly.
Started at rest at the rest position, the whole motion is x_p, because x_p and its rate are both zero at t = 0. The lab prints it in the forced-response tile.
The amplitude of this motion is 3t/4: it is 3 at t = 4 and 30 at t = 40. It grows without bound, and it is finite at every moment. It does not leap to infinity when the push is switched on; it climbs by the same amount in every unit of time. The model has no damper to stop the climb, and a real structure has one, or breaks.
Beside resonance: beats
Move the push to ω = 3/2 and keep m = 1, k = 4, F₀ = 3. The previous lesson gives A = 12/7, and from rest the motion is x = A·(cos(ωt) − cos(ω₀t)). The identity cos(θ) − cos(ψ) = 2·sin((θ + ψ)/2)·sin((ψ − θ)/2), which follows from the addition and subtraction formulas for the cosine that “Amplitude, Phase and Period” took as given, rewrites it as a product.
x = 2A·sin((ω₀ + ω)t/2)·sin((ω₀ − ω)t/2)
ω = 3/2: (ω₀ + ω)/2 = 7/4 (ω₀ − ω)/2 = 1/4
x = (24/7)·sin(7t/4)·sin(t/4)
Definition · Beat periodIn a motion 2A·sin(p·t)·sin(q·t) with q much smaller than p, the factor 2A·sin(q·t) is a slowly varying envelope for the fast factor. The size of the oscillation swells and shrinks, and the interval between one swell and the next is the beat period, π/q.
Here q = (ω₀ − ω)/2, so the beat period is 2π/|ω₀ − ω|.
For ω = 3/2 that is 2π/(1/2) = 4π, about 12.5664 and rounded, which is what the lab prints in the period tile when the push is on. The envelope (24/7)·sin(t/4) itself repeats every 8π, but the size of the swing, which is its absolute value, repeats every 4π. The largest the swing reaches is 2A = 24/7, twice the steady amplitude.
The two descriptions join. Move ω toward ω₀: the swell height 2A grows without bound and the beat period 2π/|ω₀ − ω| lengthens without bound, so the first swell stretches until it never ends. While t is small, sin((ω₀ − ω)t/2) is close to (ω₀ − ω)t/2, and the envelope is close to F₀·t/(m·(ω₀ + ω)), which tends to the resonant ramp F₀·t/(2m·ω₀) as ω → ω₀. The lab preset for ω = 7/4 shows the swell taller and slower than at 3/2; the preset at ω = 2 is the limit. That the ramp is the limit is what the lab demonstrates, and it is not proved here.