Kernel: SageMath 9.7

C17 More Resonance

Damped Harmonic Oscillator

We found last class that all solutions to

mu(t)+cu(t)+ku(t)=Csin(ωt)mu''(t) + cu'(t) + ku(t) = C\sin(\omega t)

with c>0c > 0 converge to the particular solution

u(t)=C(kmω2)(kmω2)2+(cω)2cos(ωt)+C(cω)2(kmω2)2+(cω)2sin(ωt).u(t) = \frac{C(k - m\omega^2)}{(k - m\omega^2)^2 + (c\omega)^2}\cos(\omega t) + \frac{C(c\omega)^2}{(k - m\omega^2)^2 + (c\omega)^2}\sin(\omega t).

We also saw that this particular solution could be rewritten

u(t)=CG(ω)sin(ωtΦ(ω))u(t) = CG(\omega)\sin(\omega t - \Phi(\omega))

where the gain is

G(ω)=1(kmω2)2+(cω)2G(\omega) = \frac{1}{\sqrt{(k - m\omega^2)^2 + (c\omega)^2}}

and the phase shift is

Φ(ω)=tan1(cωkmω2).\Phi(\omega) = \tan^{-1}\left(\frac{c\omega}{k - m\omega^2} \right).

Finally, we found that the gain GG is maximized at the resonance angular frequency

ωmax=km12(cm)2\omega_{\max} = \sqrt{\frac{k}{m} - \frac{1}{2}\left(\frac{c}{m}\right)^2}

which is close to the unforced angular frequency

ωnat=km(c2m)2.\omega_{\text{nat}} = \sqrt{\frac{k}{m} - \left(\frac{c}{2m}\right)^2}.

Of course, the last two definitions only make sense only if cc is sufficiently small. Otherwise, the gain is maximized at zero and the unforced ODE does not have an oscillatory solution.

Example.

Consider u(t)+u(t)+10u(t)=sin(ωt)u''(t) + u'(t) + 10u(t) = \sin(\omega t). (To get the Greek symbol ω\omega in an input cell, type \omega and the tab key.)

a. Find the particular solution for arbitrary ω\omega. Compare with our earlier work.

b. Plot the gain and phase shift on the interval 0ω50 \leq \omega \leq 5. Compare with our earlier work.

c. Find the resonant radial frequency. Compare with the previous graph.

d. Plot the gain times the forcing term G(ω)sin(ωt)G(\omega)*\sin(\omega t) and the particular solution u(t)u(t) for values of ω\omega less than, equal to, and more than the resonant radial frequency Compare with the earlier parts.

Undamped Harmonic Oscillator

Find the solution to mu(t)+ku(t)=cos(ωt)mu''(t) + ku(t) = \cos(\omega t) with initial conditions u(0)=0u(0) = 0 and u(0)=0u'(0) = 0.. Interpret.

In [0]: