Damped Oscillatory Motion

According to Eq. (167), a one-dimensional conservative system which is slightly perturbed from a stable equilibrium point (and then left alone) oscillates about this point with a fixed frequency and a constant amplitude. In other words, the oscillations never die away. This is not very realistic, since we know that, in practice, if we slightly perturb a dynamical system (such as a pendulum) from a stable equilibrium point then it will indeed oscillate about this point, but these oscillations will eventually die away due to frictional effects, which are present in all real dynamical systems. In order to model frictional effects, we need to include some sort of frictional drag force in our perturbed equation of motion, (166).

The most common model for a frictional drag force is one which is always directed in the opposite direction to the instantaneous velocity of the object upon which it acts, and is directly proportional to the magnitude of this velocity. Let us adopt this model. So, our drag force can be written

$$f_{drag} = - 2\,m\,\nu\,\frac{dx}{dt},$$ (171)

where $\nu$ is a positive constant. Including such a force in our perturbed equation of motion, (166), we obtain

$$\frac{d^2 x}{dt^2} + 2\,\nu\,\frac{dx}{dt} + \omega_0^{\,2}\,x = 0.$$ (172)

Thus, the positive constant $\nu$ parameterizes the strength of the frictional damping in our dynamical system.

Equation (172) is a linear, second-order, ordinary differential equation, which we suspect possesses oscillatory solutions. There is a standard trick for solving such an equation. We search for complex oscillatory solutions of the form

$$x = a\,{\rm e}^{-{\rm i}\,\omega\,t},$$ (173)

where both $\omega$ and $a$ are, in general, complex. Of course, the physical solution is the real part of the above expression. Note that this method of solution is only appropriate for linear differential equations. Incidentally, the method works because

$${\rm Re}[{\cal L}(x)]\equiv {\cal L}({\rm Re}[x]),$$ (174)

where $x$ is a complex variable, and ${\cal L}$ some real linear differential operator which acts on this variable.

Substituting Eq. (173) into Eq. (172), we obtain

$$a \left[-\omega^2-{\rm i}\,2\,\nu\,\omega + \omega_0^{\,2}\right] {\rm e}^{-{\rm i}\,\omega\,t} = 0,$$ (175)

which reduces to the following quadratic equation for $\omega$:

$$\omega^2 + {\rm i}\,2\,\nu\,\omega - \omega_0^{\,2} = 0.$$ (176)

The solution to this equation is

$$\omega_\pm = - {\rm i}\,\nu \pm \sqrt{\omega_0^{\,2}-\nu^2}.$$ (177)

Thus, the most general physical solution to Eq. (172) is

$$x(t) = {\rm Re}\left[a_+\,{\rm e}^{-{\rm i}\,\omega_+\,t} + a_-\,{\rm e}^{-{\rm i}\,\omega_-\,t}\right],$$ (178)

where $a_\pm$ are two arbitrary complex constants.

We can distinguish three different cases. In the first case, $\nu < \omega_0$, and the motion is said to be underdamped. The most general solution is written

$$x(t) = x_0\,{\rm e}^{-\nu\,t}\,\cos(\omega_r\,t) + \left(\fr... ...nu\,x_0}{\omega_r}\right){\rm e}^{-\nu\,t}\,\sin(\omega_r\,t),$$ (179)

where $\omega_r = \sqrt{\omega_0^{\,2}-\nu^2}$, $x_0 = x(0)$, and $v_0 = dx(0)/dt$. It can be seen that the solution oscillates at some real frequency, $\omega_r$, which is somewhat less than the natural frequency of oscillation of the undamped system, $\omega_0$, but also decays exponentially in time at a rate proportional to the damping coefficient, $\nu$.

In the second case, $\nu=\omega_0$, and the motion is said to be critically damped. The most general solution is written

$$x(t) = \left[x_0 \,(1+\omega_0\,t) + v_0\,t\right] {\rm e}^{-\omega_0\,t}.$$ (180)

It can be seen that the solution now decays without oscillating.

In the third case, $\nu>\omega_0$, and the motion is said to be overdamped. The most general solution is written

$$x(t) = -\left(\frac{v_0 + \nu_-\,x_0}{\nu_+-\nu_-}\right){\r... ...(\frac{v_0+\nu_+\,x_0}{\nu_+-\nu_-}\right){\rm e}^{-\nu_-\,t},$$ (181)

where $\nu_\pm = \nu\pm \sqrt{\nu^2-\omega_0^{\,2}}$. It can be seen that the solution again decays without oscillating, except there are now two independent decay rates. The largest, $\nu_+$, is always greater than the critically damped decay rate, $\omega_0$, whereas the smaller, $\nu_-$, is always less than this decay rate. This means that, in general, the critically damped solution is more rapidly damped than either the underdamped or overdamped solutions.

Figure 20: Damped oscillatory motion.

Figure 20 shows typical examples of underdamped (i.e., $\nu=\omega_0/4$) , critically damped (i.e., $\nu=\omega_0$), and overdamped (i.e,, $\nu=4\,\omega_0$) solutions, calculated with the initial conditions $x_0=1$ and $v_0=0$. Here, $T_0 = 2\pi/\omega_0$. The three solutions correspond to the solid, short-dashed, and long-dashed curves, respectively.