## Section 13.6 Damped Harmonic Oscillator

Any frictional force will damp the motion, but viscous drag is a particularly easier damping force to work with analytically.

where \(b\) is the proportionality constant that depends on the viscosity of the medium and the geometry of the moving object. We have a minus sign in Eq. (13.6.1) since drag opposes velocity.

Viscous forces of this type act on objects moving in a fluid, such as air, water, oil, etc., if their speeds are not too great. A viscous damping can also be implemented by pushing and pulling a plunger into a gas cylinder as shown in Figure 13.6.1(b). an arrangement called a dashpot.

### Subsection 13.6.1 Equation of Motion

We place \(x \) along the spring with positive \(x \) axis in the direction in which spring stretches and origin at the position of the block when in equilibrium. With these choices, the \(x\) component of the net force from the two forces on the damped block will be

Therefore, equation of motion of a damped block is

Solving this as a differential equation gives us all possible motions of a damped harmonic oscillator.

### Subsection 13.6.2 Three Types of Damped Oscillators

Solutions of Eq. (13.6.3) tell us about \(x \) at an arbitrary instant \(t\text{,}\) possibly in terms of given \(x_0\) and \(v_0 \text{,}\) the position and velocity at initial instant \(t=0\text{.}\) Once again, this can be done by treating Eq. (13.6.3) as differential equation. We will just borrow the solution found by advanced mathematics.

There are three types of solutions, depending upon the values of the parameters \(m\text{,}\) \(k\text{,}\) and \(b\text{.}\) Actually, only two effective parameters are needed. These are

Sometimes, formulas are written using another damping constant, \(\beta \text{,}\) instead of \(\gamma\text{.}\) For consistency of our treatment here, we will use \(\gamma\) in our formulas.

Depending on the relative values of \(\omega \) and \(\gamma\text{,}\) we get either an oscillating behavior or some other motion. The solutions tell us that there are three types of motion possible.

### Subsection 13.6.3 Underdamped Oscillator Solution

The underdamped solution can be written in many forms, but the following form is very useful in further understanding of what happens to the block as time passes.

where \(\omega^{\prime} \) is the oscilaltion angular frequency of the underdamped oscillator, given by

which shows that the underdamped oscillator oscillates at another frequency than the natural frequency, \(\omega\text{,}\) which was the frequency at which the non-damped oscillator oscillated. The constants \(A \) and \(\phi \) are related to the initial \(x_0\) and \(v_0 \text{.}\) They are very complicated and not very illuminating. Here, I will just list the equations from which you can deduce \(A \) and \(\phi\text{.}\)

Figure 13.6.2 shows a plot of \(x \) versus \(t\) for an underdamped oscillator. Clearly, an underdamped oscillator oscillates with decreasing amplitude.

### Subsection 13.6.4 Overdamped Solution

The overdamped solution shows only the exponential damping and no oscillations as illustrated in Figure 13.6.3. The solution consists of two terms that decay at different rates.

where, \(a \) and \(b\) are

where

Since \(\dfrac{\gamma}{2} \gt \omega \gt \omega^{\prime\prime} \) for overdamping, both \(a \) and \(b\) are positive numbers.

### Subsection 13.6.5 Critically Damped Solution

The critically damped solution shows a combination of exponential damping and a linearly increasing solution that is also damped exponentially. That means there is no oscillation here either.

Figure 13.6.5 illustrates critically damped oscillator. Visually, critically damped and overdamped oscillators appear similar, but when you plot the overdamped and critically damped oscillators under same initial conditions, as in Figure 13.6.5, you would find that critically damped oscillator does a little better job of damping the motion.

A very important application of critically damped oscillation is in car shock absorbers. A shock absorber provides critical damping to the springs of the suspension system and damps a car's motion smoothly and rapidly. Without a shock absorber, a car ride would be very bumpy and dangerous. The picture to the right is a rear shock absorber and spring of a BMW R75/5 motorcycle, where the shock absorber is in the inside space of the spring (Photo credits: uploaded by Jeff Dean for Wikicommons.)

###### Checkpoint 13.6.6. Deciding Type of Damped Harmonic Oscillator.

Decide what type of damping is in the following oscillators.

(a)-(d) First work out the effective parameters \(\omega\) and \(\gamma\text{.}\)

(a) Overdamped, (b) Overdamped, (c) Underdamped, (d) Critically damped.

(a) Let us first find the two effective parameters.

Since \(\omega \lt \gamma/2\text{,}\) the oscillator is overdamped.

(b) Let us first find the two effective parameters.

Since \(\omega \lt \gamma/2\text{,}\) the oscillator is overdamped.

(c) Let us first find the two effective parameters.

Since \(\omega \gt \gamma/2\text{,}\) the oscillator is underdamped.

(d) Let us first find the two effective parameters.

Since \(\omega = \gamma/2\text{,}\) the oscillator is critically damped.

###### Checkpoint 13.6.7. Lightly Damped Oscillator Solution and Envelop.

The position of a 250-gram lightly damped oscillator is given by the following function of time, \(x(t) = 2\exp{(-0.01\,t)}\cos(4\pi t)\text{,}\) where \(t\) is in seconds and \(x\) in meters.

(a) What is the frequency of oscillations?

(b) Plot \(x\) vs \(t\) for at least 4 cycles.

(c) How long does it take for the envelop of oscillations to drop by \(\dfrac{1}{e}\text{?}\)

(d) How long does it take the envelop to drop by a factor \(\dfrac{1}{e^2}\text{?}\)

(a) Look at the argument of the cosine. (b) plot \(t=0 \) to \(t=2\text{.}\) (c) Set \(0.01\, t = 1\text{.}\) (d) Set \(0.01\, t = 2\text{.}\)

(a) \(2\text{ Hz}\text{,}\) (b) See solution, (c) \(100\text{ sec}\text{,}\) (d) \(200\text{ sec}\text{.}\)

(a) From the argument of the cosine we can find the period of oscillations. The inverse of this period is the frequency. The period of \(\cos(4\pi t)\) is \(T= 0.5\text{ sec}\text{.}\) Hence, the frequency of the underdamped oscillatory is

(b) Plot the given function for the range \(t=[0,\ 2]\text{.}\)

(c) The envelop is dropping in time as \(\exp{(-0.01\,t)}\text{.}\) We need to find the value of \(t\) at which the value becomes \(\exp{(-1)}\text{.}\) This means

(d) The envelop is dropping in time as \(\exp{(-0.01\,t)}\text{.}\) We need to find the value of \(t\) at which the value becomes \(\exp{(-2)}\text{.}\) This means

###### Checkpoint 13.6.8. Displacement of a Harmonic Oscillator Given Initial Displacement and Velocity.

A block of mass \(2\text{ kg}\) is attached to a spring of spring constant \(32\text{ N/m}\) is oscillating with damping constant \(8\text{ kg/s}\text{.}\) The block started at \(x = 2\text{ cm}\) from rest. Where is the block at an arbitrary instant \(t\text{?}\)

First find out what type of damping it is and then use the general solution for \(x\) and \(v\) to fit the initial condition given.

\(3.1\text{ cm}\,\exp(-4 t)\,\cos(\sqrt{12} t -0.86)\text{.}\)

We have three formulas for \(x\) to choose from, depending upon the type of damping. The effective parameters are

Since \(\omega \gt \gamma/2\text{,}\) the damping is an underdamping with oscillation angular frequency

and \(x\) and \(v\) as:

Now, we use the given initial conditions to get

The second equation gives \(\phi=-0.857\text{ rad}\) which is used in the first to get \(A=3.1\text{ cm}\text{.}\)

Therefore, we have the following \(x\) at arbitrary instant.

where \(t\) is in \(\text{sec}\text{.}\)