next up previous
Next: Driven Circuit Up: Damped and Driven Harmonic Previous: Circuit


Driven Damped Harmonic Oscillation

We saw earlier, in Section 3.1, that when a damped mechanical oscillator is set into motion the oscillations eventually die away due to frictional energy losses. In fact, the only way of maintaining the amplitude of a damped oscillator is to continuously feed energy into the system in such a manner as to offset the frictional losses. A steady (i.e., constant amplitude) oscillation of this type is called driven damped harmonic oscillation. Consider a modified version of the mass-spring system investigated in Section 3.1 in which one end of the spring is attached to the mass, and the other to a moving piston. See Figure 8. Let $x(t)$ be the horizontal displacement of the mass, and $X(t)$ the horizontal displacement of the piston. The extension of the spring is thus $x(t)-X(t)$, assuming that the spring is unstretched when $x=X=0$. Thus, the horizontal force acting on the mass can be written [cf., Equation (55)]
\begin{displaymath}
f = - k\,(x-X)-m\,\nu\,\dot{x}.
\end{displaymath} (92)

The equation of motion of the system then becomes [cf., Equation (56)]
\begin{displaymath}
\ddot{x} + \nu\,\dot{x} + \omega_0^{\,2}\,x = \omega_0^{\,2}\,X,
\end{displaymath} (93)

where $\nu>0$ is the damping constant, and $\omega_0>0$ the undamped oscillation frequency. Suppose, finally, that the piston executes simple harmonic oscillation of angular frequency $\omega>0$ and amplitude $X_0>0$, so that the time evolution equation of the system takes the form
\begin{displaymath}
\ddot{x} + \nu\,\dot{x} + \omega_0^{\,2}\,x = \omega_0^{\,2}\,X_0\,\cos(\omega\,t).
\end{displaymath} (94)

We shall refer to the above as the driven damped harmonic oscillator equation.

Figure 8: A driven oscillatory system
\begin{figure}
\epsfysize =2.in
\centerline{\epsffile{Chapter03/fig03.eps}}
\end{figure}

Now, we would generally expect the periodically driven oscillator shown in Figure 8 to eventually settle down to a steady (i.e., constant amplitude) pattern of oscillation, with the same frequency as the piston, in which the frictional energy loss per cycle is exactly matched by the work done by the piston per cycle (see Exercise 7). This suggests that we should search for a solution to Equation (94) of the form

\begin{displaymath}
x(t)=x_0\,\cos(\omega\,t-\varphi).
\end{displaymath} (95)

Here, $x_0>0$ is the amplitude of the driven oscillation, whereas $\varphi$ is the phase lag of this oscillation (with respect to the phase of the piston oscillation). Now, since
$\displaystyle \dot{x}$ $\textstyle =$ $\displaystyle -\omega\,x_0\,\sin(\omega\,t-\varphi),$ (96)
$\displaystyle \ddot{x}$ $\textstyle =$ $\displaystyle -\omega^2\,x_0\,\cos(\omega\,t-\varphi),$ (97)

Equation (94) becomes
\begin{displaymath}
(\omega_0^{\,2}-\omega^2)\,x_0\,\cos(\omega\,t-\varphi) - \n...
...n(\omega\,t-\varphi)
= \omega_0^{\,2}\,X_0\,\cos(\omega\,t).
\end{displaymath} (98)

However, $\cos(\omega\,t-\varphi)\equiv \cos(\omega\,t)\,\cos\varphi + \sin(\omega\,t)\,\sin\varphi$ and $\sin(\omega\,t-\varphi) \equiv \sin(\omega\,t)\,\cos \varphi - \cos(\omega\,t)\,\sin\varphi$, so we obtain
$\displaystyle \left[x_0\,(\omega_0^{\,2}-\omega^2)\,\cos\varphi + x_0\,\nu\,\omega\,\sin\varphi -\omega_0^{\,2}\,X_0\right]\cos(\omega\,t)$      
$\displaystyle +x_0\left[(\omega_0^{\,2}-\omega^2)\,\sin\varphi-\nu\,\omega\,\cos\varphi\right]\sin(\omega\,t)$ $\textstyle =$ $\displaystyle 0.$ (99)

Now, the only way in which the above equation can be satisfied at all times is if the coefficients of $\cos(\omega\,t)$ and $\sin(\omega\,t)$ separately equate to zero. In other words,
$\displaystyle x_0\,(\omega_0^{\,2}-\omega^2)\,\cos\varphi + x_0\,\nu\,\omega\,\sin\varphi -\omega_0^{\,2}\,X_0$ $\textstyle =$ $\displaystyle 0,$ (100)
$\displaystyle (\omega_0^{\,2}-\omega^2)\,\sin\varphi-\nu\,\omega\,\cos\varphi$ $\textstyle =$ $\displaystyle 0.$ (101)

These two expressions can be combined to give
$\displaystyle x_0$ $\textstyle =$ $\displaystyle \frac{\omega_0^{\,2}\,X_0}{\left[(\omega_0^{\,2}-\omega^2)^2+\nu^2\,\omega^2\right]^{1/2}},$ (102)
$\displaystyle \varphi$ $\textstyle =$ $\displaystyle \tan^{-1}\left(\frac{\nu\,\omega}{\omega_0^{\,2}-\omega^2}\right).$ (103)

This follows because (101) gives
\begin{displaymath}
\tan\varphi = \frac{\nu\,\omega}{\omega_0^{\,2}-\omega^2},
\end{displaymath} (104)

and so
$\displaystyle \cos\varphi$ $\textstyle \equiv$ $\displaystyle \frac{1}{\sqrt{1+\tan^2\varphi}} = \frac{\omega_0^{\,2}-\omega^2}{\left[(\omega_0^{\,2}-\omega^2)^2+\nu^2\,\omega^2\right]^{1/2}},$ (105)
$\displaystyle \sin\varphi$ $\textstyle \equiv$ $\displaystyle \frac{\tan\varphi}{\sqrt{1+\tan^2\varphi}} = \frac{\nu\,\omega}{\left[(\omega_0^{\,2}-\omega^2)^2+\nu^2\,\omega^2\right]^{1/2}}.$ (106)

Hence, substitution into (100) gives (102).

Figure 9: Driven harmonic motion.
\begin{figure}
\epsfysize =2.25in
\centerline{\epsffile{Chapter03/fig04.eps}}
\end{figure}

Let us investigate the dependence of the amplitude, $x_0$, and phase lag, $\varphi$, of the driven oscillation on the driving frequency, $\omega $. This is most easily done graphically. Figure 9 shows $x_0/X_0$ and $\varphi$ plotted as functions of $\omega $ for various different values of $\nu/\omega_0$. In fact, $\nu/\omega_0 = 1/Q_f= 1$, $1/2$, $1/4$, $1/8$, and $1/16$ correspond to the solid, dotted, short-dashed, long-dashed, and dot-dashed curves, respectively. It can be seen that as the amount of damping in the system is decreased the amplitude of the response becomes progressively more peaked at the natural frequency of oscillation of the system, $\omega_0$. This effect is known as resonance, and $\omega_0$ is termed the resonant frequency. Thus, a weakly damped oscillator (i.e., $\nu\ll \omega_0$) can be driven to large amplitude by the application of a relatively small amplitude external driving force which oscillates at a frequency close to the resonant frequency. Note that the response of the oscillator is in phase (i.e., $\varphi\simeq 0$) with the external drive for driving frequencies well below the resonant frequency, is in phase quadrature (i.e., $\varphi=\pi/2$) at the resonant frequency, and is in anti-phase (i.e., $\varphi\simeq \pi$) for frequencies well above the resonant frequency.

According to Equations (82) and (102),

\begin{displaymath}
\frac{x_0(\omega=\omega_0)}{X_0} = \frac{\omega_0}{\nu}
=Q_f.
\end{displaymath} (107)

In other words, when the driving frequency matches the resonant frequency the ratio of the amplitude of the driven oscillation to that of the piston oscillation is the quality factor, $Q_f$. Hence, $Q_f$ can be regarded as the resonant amplification factor of the oscillator. Equations (102) and (106) imply that, for a weakly damped oscillator (i.e., $\nu\ll \omega_0$) which is close to resonance [i.e., $\vert\omega-\omega_0\vert\sim \nu\ll \omega_0$],
\begin{displaymath}
\frac{x_0(\omega)}{x_0(\omega=\omega_0)}\simeq \sin\varphi\simeq \frac{\nu}{[4\,(\omega_0-\omega)^2 + \nu^2]^{1/2}}.
\end{displaymath} (108)

This follows because $\omega_0^{\,2}-\omega^2=(\omega_0+\omega)\,(\omega_0-\omega)\simeq 2\,\omega_0\,(\omega_0-\omega)$. Hence, the width of the resonance peak (in angular frequency) is $\Delta\omega = \nu$, where the edges of peak are defined as the points at which the driven amplitude is reduced to $1/\sqrt{2}$ of its maximum value: i.e., $\omega=\omega_0\pm \nu/2$. Note that the phase lag at the low and high frequency edges of the peak are $\pi/4$ and $3\pi/4$, respectively. Furthermore, the fractional width of the peak is
\begin{displaymath}
\frac{\Delta \omega}{\omega_0} = \frac{\nu}{\omega_0} = \frac{1}{Q_f}.
\end{displaymath} (109)

We conclude that the height and width of the resonance peak of a weakly damped ($Q_f\gg 1$) harmonic oscillator scale as $Q_f$ and $Q_f^{-1}$, respectively. Thus, the area under the resonance peak stays approximately constant as $Q_f$ varies.

Now, the force exerted on the system by the piston is

\begin{displaymath}
F(t) = k\,X_0\,\cos(\omega\,t).
\end{displaymath} (110)

Hence, the instantaneous power absorption from the piston becomes
$\displaystyle P(t)$ $\textstyle =$ $\displaystyle F(t)\,\dot{x}(t)$  
  $\textstyle =$ $\displaystyle - k\,X_0\,x_0\,\omega\,\cos(\omega\,t)\,\sin(\omega\,t-\varphi)$  
  $\textstyle =$ $\displaystyle -k\,X_0\,x_0\,\omega\left[\cos(\omega\,t)\,\sin(\omega\,t)\,\cos\varphi - \cos^2(\omega\,t)\,\sin\varphi\right].$ (111)

Thus, the average power absorption during an oscillation cycle is
\begin{displaymath}
\langle P\rangle = \frac{1}{2}\,k\,X_0\,\,x_0\,\omega\,\sin\varphi,
\end{displaymath} (112)

since $\langle \cos(\omega\,t)\,\sin(\omega\,t)\rangle =0$ and $\langle \cos^2(\omega\,t)\rangle =1/2$. Of course, given that the amplitude of the driven oscillation neither grows nor decays, the average power absorption from the piston during an oscillation cycle must be equal to the average power dissipation due to friction (see Exercise 7). Making use of Equations (107) and (108), the mean power absorption when the driving frequency is close to the resonant frequency is
\begin{displaymath}
\langle P\rangle \simeq \frac{1}{2}\,\omega_0\,k\,X_0^{\,2}\,Q_f\left[\frac{\nu^2}{4\,(\omega_0-\omega)^2+\nu^2}\right].
\end{displaymath} (113)

Thus, the maximum power absorption occurs at the resonance (i.e., $\omega=\omega_0$), and the absorption is reduced to half of this maximum value at the edges of the resonance (i.e., $\omega=\omega_0\pm \nu/2$). Furthermore, the peak power absorption is proportional to the quality factor, $Q_f$, which means that it is inversely proportional to the damping constant, $\nu$.

Figure 10: A driven $LCR$ circuit.
\begin{figure}
\epsfysize =2.in
\centerline{\epsffile{Chapter03/fig05.eps}}
\end{figure}


next up previous
Next: Driven Circuit Up: Damped and Driven Harmonic Previous: Circuit
Richard Fitzpatrick 2010-10-11