$\newcommand{\trace}{\operatorname{tr}} \newcommand{\real}{\operatorname{Re}} \newcommand{\imaginary}{\operatorname{Im}} \newcommand{\lt}{<} \newcommand{\gt}{>} \newcommand{\amp}{&}$

## Section4.4Forcing and Resonance

Resonance happens when the natural frequency and the forcing frequency of an undamped harmonic oscillator are the same. There are many real world examples of resonance. For example, it is possible for a trained singer to shatter a champagne glass by simply singing the correct note (see https://www.youtube.com/watch?v=IZD8ffPwXRo). Shattering a glass by singing the correct note can be modeled by an equation of an undamped harmonic oscillator with sinusoidal forcing. While a certain amount of damping occurs in all physical systems, the damping is often so small that it makes sense to use an undamped harmonic oscillator.

### Subsection4.4.1Resonance

Resonance was responsible for the collapse of the Broughton suspension bridge near Manchester, England in 1831. The collapse occurred when a column of soldiers marched in cadence over the bridge, setting up a periodic force of rather large amplitude. The frequency of the force was approximately equal to the natural frequency of the bridge. Thus, the bridge collapsed when large oscillations occurred. For this reason soldiers are ordered to break cadence whenever they cross a bridge.

The Millennium Bridge, the first new bridge to span the Thames River in London in over 100 years, is a modern example of how resonance can effect a bridge (Figure 4.4.1). This pedestrian bridge, which opened to the public in June 2000, was quickly closed after the bridge experienced high amplitude horizontal oscillations during periods of high traffic. Studies by designers found that the bridge experienced high amplitude horizontal oscillations in response to horizontal forcing at a rate of one cycle per second. Typically, people walk at a rate of two steps per second, so the time between two successive steps of the left foot is about one second. Thus, if people were to walk in cadence, they could set up strong horizontal forcing that would place a destructive load on the bridge. The engineers did not envision this to be a problem since tourists do not generally march in time. However, a video of tourists crossing the bridge revealed the opposite. When the bridge began oscillating, people tended to walk in cadence in order to keep their balance (https://www.youtube.com/watch?v=gQK21572oSU).

Consider the case of an undamped harmonic oscillator with a sinusoidal forcing function. We know that we can model this situation with the second-order linear differential equation

\begin{equation*} x'' + \omega_0^2 x = A \cos \omega t. \end{equation*}

Since the homogeneous solution to this equation is

\begin{equation*} x_h(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t, \end{equation*}

we can see that the natural frequency of the harmonic oscillator is $\omega_0\text{.}$ This is the frequency of the oscillations if there is no forcing term. Depending on whether or not the frequency of the forcing term $\omega$ is equal to the natural frequency, we will have two very different situations. If the two frequencies are equal, then they will re-enforce one another and we will have resonance. If the two frequencies are close but not equal, we will see the phenomenon of beats. We wish to investigate what happens when the driving frequency $\omega$ is equal to the natural frequency and when it is not.

Let us examine the case where the forcing frequency and the natural frequency of the oscillator are the same,

$$x'' + \omega_0^2 x = A \cos \omega_0 t.\label{equation-secondorder04-resonance-5}\tag{4.4.1}$$

Since $A \cos \omega_0 t$ is a solution to the homogeneous equation $x'' + \omega_0^2 x = 0\text{,}$ we cannot assume that a particular solution to equation (4.4.1) has the form $a \cos \omega_0 t + b \sin \omega_0 t\text{.}$ Equivalently, if we use the complex method, we cannot assume that our solution has the form $ae^{i \omega_0t}\text{.}$ Therefore, we will look for a complex solution of the form

\begin{equation*} x_c = a t e^{i \omega_0 t} \end{equation*}

for the complex equation that corresponds to (4.4.1),

$$x'' + \omega_0^2 x = A e^{i \omega_0 t}.\label{equation-secondorder04-resonance-6}\tag{4.4.2}$$

In this case,

\begin{align*} x_c' & = a ( 1 + i \omega_0 t )e^{i \omega_0 t} \\ x_c'' & = a (2 i \omega_0 - \omega_0^2 t) e^{i \omega_0 t}. \end{align*}

Substituting $x_c$ and $x_c''$ into the left-hand side of (4.4.2), we have

\begin{align*} x_c'' + \omega_0^2 x_c & = a (2 i \omega_0 - \omega_0^2 t )e^{i \omega_0 t} + \omega_0^2 a t e^{i \omega_0 t}\\ & = 2ai \omega_0 e^{i \omega t}. \end{align*}

In order for $x_c = a t e^{i \omega_0 t}$ to be a solution to (4.4.2), we must have

\begin{equation*} a = \frac{A}{2 i \omega_0}. \end{equation*}

Thus, our solution to the complex form of the differential equation is

\begin{align*} x_c & = \frac{A}{2i \omega_0} t e^{i \omega_0 t}\\ & = - \frac{Ai}{2 \omega_0} t e^{i \omega_0 t}\\ & = - \frac{Ai}{2 \omega_0} t (\cos \omega_0 t + i \sin \omega_0 t)\\ & = \frac{A}{2 \omega_0} t \sin \omega_0 t - i \frac{A}{2 \omega_0} t \cos \omega_0 t. \end{align*}

The real part of $x_c\text{,}$

\begin{equation*} x_p(t) = \frac{A}{2 \omega_0} t \sin \omega_0 t, \end{equation*}

is a particular solution to (4.4.1). Thus, our general solution is

\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t + \frac{A}{2 \omega_0} t \sin \omega_0 t. \end{equation*}
###### Example4.4.2

Now let us consider the initial-value problem

\begin{align*} x'' + 144 x & = 4 \cos 12t\\ x(0) & = 0\\ x'(0) & = 0. \end{align*}

where $A = 4$ and $\omega_0 = 12\text{.}$ The solution to the homogeneous equation $x'' + 144 x = 0$ is

\begin{equation*} x_h(t) = c_1 \cos 12t + c_2 \sin 12t. \end{equation*}

To find a particular solution to $x'' + 144 x = 4 \cos 12t\text{,}$ we will use the complex method and try to find a particular solution to $x'' + 144 x = 4 e^{12it}\text{.}$ We must assume that the solution has the form $x_c(t) = a t e^{12it}\text{,}$ since $a e^{12it}$ is a solution to the homogeneous equation. As before, we have $x_c'' = a (24i - 144 t) e^{12it}\text{.}$ If we substitute $x_c$ and $x_c''$ into the left-hand side of our differential equation, we have

\begin{equation*} x_c'' + 144 x_c = a (24i - 144 t) e^{12it} + 144 a t e^{12i t} = 24ai e^{12i t}. \end{equation*}

Therefore,

\begin{equation*} a = \frac{4}{24i} = - \frac{1}{6}i, \end{equation*}

and our complex solution is

\begin{equation*} - \frac{1}{6} i e^{12it} = - \frac{1}{6} i (\cos 12t + i \sin 12 t) = \frac{1}{6} \sin 12 t - \frac{1}{6} i \cos 12 t. \end{equation*}

Taking the real part of our complex solution, we have a particular solution

\begin{equation*} x_p(t) = \frac{1}{6} t \sin 12t. \end{equation*}

Thus, the general solution to $x'' + 144 x = 4 \cos 12t$ is

\begin{equation*} x(t) = x_h(t) + x_p(t)= c_1 \cos 12t + c_2 \sin 12t + \frac{1}{6} t \sin 12t. \end{equation*}

Applying the initial conditions, both $c_1 = 0$ and $c_2 = 0\text{.}$ Consequently, the solution to the initial-value problem is

\begin{equation*} x(t) = \frac{1}{6} t \sin 12t. \end{equation*}

The graph of this solution is given in Figure 4.4.3. Notice that our solution grows with time. This growth is due to the fact that the frequency of the forcing term is equal to the natural frequency of the oscillator. Since the force pulls and pushes at a frequency equal to the natural frequency of the oscillator, the amplitude increases with time. This type of behavior is called resonance.

### Subsection4.4.2Beats or the Case $\omega \neq \omega_0$

If the forcing frequency and natural frequency are not equal but close, then we have the phenomenon of beats. For example, a piano tuner uses this phenomenon to make certain that a particular string is correctly tuned. The tuner will strike a tuning fork which vibrates at the correct frequency. Next, the tuner hits the poorly tuned piano string. This modulation gives rise to beats in the tone that are readily audible. When the string is properly tuned, the beats will go away. To see another example of beats, watch the video https://www.youtube.com/watch?v=pRpN9uLiouI.

To understand the situation where the forcing and natural frequencies differ by a small amount, $\omega \neq \omega_0\text{,}$ we will first consider the equation,

$$x'' + \omega_0^2 = A \cos \omega t.\label{equation-secondorder04-beats-1}\tag{4.4.3}$$

The solution to the homogeneous equation $x'' + \omega_0^2 x = 0$ is

\begin{equation*} x_h(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t. \end{equation*}

Since the forcing frequency is not equal to the natural frequency, we can look for a particular solution to the complex of the form of equation (4.4.3),

$$x'' + \omega_0^2 = A e^{i \omega t}.\label{equation-secondorder04-beats-2}\tag{4.4.4}$$

If we assume out solution has the form $x_c = a e^{i\omega t}$ and substitute $x_c$ into the left-hand side of (4.4.4), we find

\begin{align*} x_c'' + \omega_0^2 x_c \amp = -a \omega^2 e^{i \omega t} + \omega_0^2 a e^{i \omega t}\\ \amp = a(\omega_0^2 - \omega^2)e^{i \omega t}. \end{align*}

Therefore, for $x_c = a e^{i\omega t}$ to be a solution for

\begin{equation*} x'' + \omega_0^2 x = Ae^{i\omega t}, \end{equation*}

we must have

\begin{equation*} a = \frac{A}{\omega_0^2 - \omega^2}. \end{equation*}

The real part of $x_c$ is a particular solution to (4.4.3),

\begin{equation*} x_p(t) = \frac{A}{\omega_0^2 - \omega^2} \cos \omega t. \end{equation*}

Thus, the general solution to (4.4.3) is

\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t + \frac{A}{\omega_0^2 - \omega^2} \cos \omega t. \end{equation*}

Now let us examine the case where the motion starts at equilibrium. That is, we will investigate what happens when $x(0) = 0$ and $x'(0) = 0\text{.}$ In this case, we can easily determine that $c_1 = -A/(\omega_0^2 - \omega^2)$ and $c_2 = 0\text{.}$ Hence the solution to our initial value problem is

$$x(t) = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t).\label{equation-secondorder04-resonance-2}\tag{4.4.5}$$

Thus, we have a superposition to two oscillations of different frequencies.

To understand how the superposition of the two frequencies works, we will let

\begin{equation*} \overline{\omega} = \frac{\omega_0 + \omega}{2} \end{equation*}

be the mean frequency and

\begin{equation*} \delta = \frac{\omega_0 - \omega}{2} \end{equation*}

be the half difference. Since $\omega = \overline{\omega} - \delta$ and $\omega_0 = \overline{\omega} + \delta\text{,}$ we know that

\begin{align*} x(t) & = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t)\\ & = \frac{A}{(\overline{\omega} + \delta)^2 - (\overline{\omega} - \delta)^2} (\cos [(\overline{\omega} - \delta) t] - \cos[(\overline{\omega} + \delta) t])\\ & = \frac{A}{4\overline{\omega} \delta} [\cos(\overline{\omega} t - \delta t) - \cos(\overline{\omega}t + \delta t)]\\ & = \frac{A}{4\overline{\omega} \delta} [\cos(\overline{\omega} t) \cos( \delta t) + \sin(\overline{\omega} t) \sin( \delta t) - \cos(\overline{\omega} t) \cos( \delta t) + \sin(\overline{\omega} t) \sin( \delta t)]\\ & = \frac{A}{4\overline{\omega} \delta} [2 \sin(\overline{\omega} t) \sin( \delta t)]\\ & = \frac{A \sin \delta t}{2 \overline{\omega} \delta} \sin \overline{\omega} t. \end{align*}

Thus, we have a sine function, $\sin \overline{\omega} t$ whose amplitude varies according to $|(A/2 \overline{\omega} \delta) \sin \delta t|\text{,}$ and

$$x(t) = \frac{A \sin \delta t}{2 \overline{\omega} \delta} \sin \overline{\omega} t.\label{equation-secondorder04-resonance-3}\tag{4.4.6}$$
###### Example4.4.4

In the initial-value problem

\begin{align} x'' + 144 x \amp = 4 \cos(13t)\label{equation-secondorder04-resonance-7}\tag{4.4.7}\\ x(0) \amp = 0\tag{4.4.8}\\ x'(0) \amp = 0,\tag{4.4.9} \end{align}

we have $A= 4\text{,}$ $\omega_0 = 12\text{,}$ and $\omega = 13\text{.}$ The solution to this equation is

\begin{equation*} x(t) = \frac{4}{25}(\cos 12 t - \cos 13 t). \end{equation*}

The graph of this solution is given in Figure 4.4.5 and illustrate the phenomenon of beats. The two frequencies are almost equal and interfere with each other.

In our example, $\delta = -1/2$ and $\overline{\omega} = 25/2\text{.}$ The factor $\sin \delta t = \sin( - t/2)$ in (4.4.6) oscillates very slowly in comparison to $\sin(\overline{\omega}) = \sin(25t/2)\text{.}$ Thus, we can view the solution (4.4.6) as a fast oscillation with a frequency $\overline{\omega}$ and amplitude

$$\left| \frac{A \sin \delta t}{2 \overline{\omega} \delta} \right| = \left| \frac{8}{25} \sin\left( \frac{t}{2}\right)\right|,\label{equation-secondorder04-resonance-4}\tag{4.4.10}$$

which oscillates much more slowly. If we superimpose the graph of (4.4.10) onto the graph in Figure 4.4.5, we obtain an envelope of the faster oscillation.

### Subsection4.4.3Forced Damped Harmonic Motion

Since perpetual motion does not really exist, we should consider what happens when we add some damping to our harmonic oscillator,

$$x'' + 2c x' + \omega_0^2 x = A \cos \omega t.\label{equation-secondorder04-force-damped-1}\tag{4.4.11}$$
###### Example4.4.7

Let us consider the differential equation

$$x'' + 2x' + 2x = 3 \cos 4t.\label{equation-secondorder04-example-0}\tag{4.4.12}$$

The corresponding homogeneous equation is

\begin{equation*} x'' + 2x' + 2x = 0, \end{equation*}

which has solution

\begin{equation*} x_h(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t. \end{equation*}

To find a particular solution for the equation (4.4.12), we will look for a solution of the form $x_c = a e^{4 i t}$ to the equation

$$x'' + 2x' + 2x = 3e^{4 it}\label{equation-secondorder04-example-3}\tag{4.4.13}$$

and then set $x_p(t)$ to the real part of the solution. Substituting $x_c$ into the lefthand side of (4.4.13), we obtain

\begin{equation*} x_c'' + 2x_c' + 2 x_c = (-14 + 8i) a e^{4it} = P(4i)a e^{4it} = 3e^{4 it}, \end{equation*}

where $P(\lambda) = \lambda^2 + 2 \lambda + 2$ is the characteristic polynomial of $x'' + 2x' + 2x = 0\text{.}$ Therefore,

\begin{equation*} a = \frac{3}{P(4i)} = \frac{3}{-14 + 8i} = - \frac{3}{130}(7 + 4i), \end{equation*}

and the solution to the complex equation is

\begin{align*} x_c \amp = - \frac{3}{130}(7 + 4i)e^{4 i t}\\ \amp = - \frac{3}{130}(7 + 4i)(\cos 4t + i \sin 4t)\\ \amp = \frac{3}{130}[(- 7 \cos 4t + 4 \sin 4t) + i(- 4 \cos 4t - 7 \sin 4t)]. \end{align*}

Thus, we can determine a particular solution by taking the real part of $x_c\text{,}$

\begin{equation*} x_p(t) = \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}

The general solution for (4.4.12) is

\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t + \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}

It is easy to see that all solutions approach $x_p(t)$ as $t \to \infty$ no matter what the initial conditions may be. However, it would be much more useful if we would write the particular solution in the form

\begin{equation*} x_p(t) = K \cos(4 t - \phi ) \end{equation*}

for some constants $K$ and $\phi$ rather than a linear combination of $\cos 4t$ and $\sin 4t\text{.}$

Returning to the general case, the associated homogeneous equation for (4.4.11) is

$$x'' + 2c x' + \omega_0^2 x = 0.\label{equation-secondorder04-force-damped-2}\tag{4.4.14}$$

The characteristic polynomial of (4.4.14),

$$P(\lambda ) = \lambda^2 + 2c \lambda + \omega_0^2,\label{equation-secondorder04-force-damped-3}\tag{4.4.15}$$

has roots

\begin{equation*} \lambda = - c \pm \sqrt{c^2 - \omega_0^2}. \end{equation*}

If our harmonic oscillator is undedamped ($c \lt \omega_0$), then the solution to the homogeneous equation (4.4.14) is

\begin{equation*} x_h = e^{-ct} (c_1 \cos \beta t + c_2 \sin \beta t), \end{equation*}

where $\beta = \sqrt{\omega_0^2 - c^2}\text{.}$

As in Example 4.4.7, we will look for a particular solution to (4.4.11) by examining the equivalent complex differential equation. That is, we shall look for a solution of the form $x_c = a e^{i \omega t}$ to the equation

$$x_c'' + 2c x_c' + \omega_0^2 x_c = Ae^{i \omega t},\label{equation-secondorder04-force-damped-4}\tag{4.4.16}$$

and then set $x_p$ equal to the real part of our solution. Substituting $x_c$ into the right-hand side of (4.4.16), we obtain

\begin{equation*} x_c'' + 2c x_c' + \omega_0^2 x_c = [(i\omega)^2 + 2c(i \omega) + \omega_0^2] a e^{i \omega t} = P(i \omega) a e^{i \omega t}, \end{equation*}

where $P$ is the characteristic polynomial (4.4.15). Thus, equation (4.4.16) becomes

\begin{equation*} P(i \omega) a e^{i \omega t} = Ae^{i \omega t}, \end{equation*}

and

\begin{equation*} x_c(t) = a e ^{i \omega t} = \frac{A}{P(i \omega)} e^{i \omega t} = H(i \omega) A e^{i \omega t}, \end{equation*}

where

\begin{equation*} H(i \omega) = \frac{1}{P(i \omega)}. \end{equation*}

We say that $H(\lambda) = 1 / P(\lambda)$ is the transfer function.

Let us examine the transfer function more closely to see if we can use it to write our particular solution in a more useable form. First, let us write

\begin{equation*} P(i \omega ) = (i\omega)^2 + 2c(i \omega) + \omega_0^2 = (\omega_0^2 - \omega^2) + 2ic \omega. \end{equation*}

in polar form,

\begin{equation*} P(i \omega) = R e^{i \phi} = R(\cos \phi + i \sin \phi), \end{equation*}

where

\begin{equation*} R = \sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2} \end{equation*}

and $\phi$ is the angle defined by the equations

\begin{align*} \cos \phi & = \frac{\omega_0^2 - \omega^2}{\sqrt{(\omega_0^2 - \omega^2)^2 - 4 c^2 \omega^2}},\\ \sin \phi & = \frac{2c\omega}{\sqrt{(\omega_0^2 - \omega^2)^2 - 4 c^2 \omega^2}}. \end{align*}

Since $2 c \omega \gt 0\text{,}$ we know that $\sin \phi \gt 0\text{.}$ Equivalently, $0 \lt \phi \lt \pi\text{.}$ Thus,

\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right). \end{equation*}

Therefore, we can write the transfer function as

\begin{equation*} H(i \omega) = \frac{1}{P(i \omega)} = \frac{1}{R} e^{-i \phi}. \end{equation*}

We define the gain to be

\begin{equation*} G(\omega) = \frac{1}{R} = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}, \end{equation*}

and we will rewrite the transfer function as

\begin{equation*} H(i \omega) = G(\omega) e^{-i \phi(\omega)}. \end{equation*}

Thus, the solution to

\begin{equation*} x_c'' + 2c x_c' + \omega_0^2 x_c = Ae^{i \omega t} \end{equation*}

is

$$x_c(t) = H(i \omega) A e^{i \omega t} = G(\omega) A e^{i(\omega t - \phi)}.\label{equation-secondorder04-force-damped-complex}\tag{4.4.17}$$

Taking the real part of (4.4.17), our particular solution is

\begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi ). \end{equation*}

It is now clear that particular solution has the same frequency as the forcing term. In addition, $x_p$ is out of phase with the driving force by the amount

\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega}\right). \end{equation*}

The general solution to

\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}

is

\begin{equation*} x(t) = x_h(t) + x_p(t) = e^{-ct} (c_1 \cos(\eta t) + c_2 \sin(\eta t)) + G(\omega) A \cos(\omega t - \phi). \end{equation*}

Since $x_h$ has the factor $e^{-ct}\text{,}$ the homogeneous part of the solution quickly decays to zero as $t \to \infty\text{.}$ For this reason, $x_h$ is called the transient term while $x_p$ is called the steady-state term.

###### Example4.4.8

Let us examine the steady-state solution of the harmonic oscillator in Example 4.4.7. We have already determined that the second-order linear differential equation

\begin{equation*} x'' + 2x' + 2x = 3 \cos 4t \end{equation*}

has solution

\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t + \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}

The natural frequency is $\omega_0 = \sqrt{2} \approx 1.4142$ If we choose the initial conditions so that $c_1 = 0$ and $c_2 = 1\text{,}$ our transient term is

\begin{equation*} x_h(t) = e^{-ct} \sin \beta t = e^{-t} \sin t. \end{equation*}

Since $c = 1\text{,}$ $\omega = 4\text{,}$ and $A = 3\text{,}$ the gain is $G(4) \approx 0.0620\text{,}$ the amplitude of the steady state solution is $G(4) \cdot 4 \approx 0.2481\text{.}$ The phase is $\phi = \phi(4) \approx 2.6224.$ The solution and the steady-state solution are given in Figure 4.4.9. As we can see, transient solutions can be quite large when compared to steady-state solutions.

Large transient currents can be destructive in electrical circuits. For example, a light bulb usually burns out when a large transient flows through an already weakened bulb. This, usually occurs at the time the light bulb is turned out or turned off. Large transient currents are particularly harmful to the hard drive in a computer. For this reason, computers should be left on except in the case that they are not being used for a period of several days.

### Subsection4.4.4Important Lessons

• A forced damped harmonic oscillator,
\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
has homogeneous solution
\begin{equation*} x_h = e^{-ct} (c_1 \cos \eta t + c_2 \sin \eta t), \end{equation*}
where $\eta = \sqrt{\omega_0^2 - c^2}\text{.}$ A particular solution is given by
\begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi ), \end{equation*}
where the phase angle $\phi$ is
\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right) \end{equation*}
and the gain $G$ is given by
\begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}
• The amplitude and phase of the steady-state solution,
\begin{equation*} x_p(t) = G(\omega) A \cos(\omega t - \phi), \end{equation*}
are determined by
\begin{equation*} \omega_0^2G = \frac{1}{\sqrt{(1- s^2)^2 + D^2 s^2}}, \end{equation*}
where $\omega = s \omega_0$ and $D = 2c/\omega_0\text{.}$
• In an undamped harmonic oscillators with a sinusoidal forcing term, the natural frequency of the solution interacts with the frequency of the forcing term. As the frequency of the forcing term approaches the natural frequency of the equation, we can observe a phenomenon called resonance. We can use the equation
\begin{equation*} x'' + \omega_0^2 x = A \cos \omega t \end{equation*}
to model an undamped harmonic oscillator with sinusoidal forcing.
• If the driving frequency is not equal to the natural frequency ($\omega \neq \omega_0$), we have a particular solution
\begin{equation*} x_p(t) = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t). \end{equation*}
If the two frequencies are close, we can see a phenomenon called beats.
• If the driving frequency is equal to the natural frequency, we have a particular solution
\begin{equation*} x_p(t) = \frac{A}{2 \omega_0} t \sin \omega_0 t. \end{equation*}
In this case, we can see the resonance phenomenon.
• A forced damped harmonic oscillator,
\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
has homogeneous solution
\begin{equation*} x_h = e^{-ct} (c_1 \cos \eta t + c_2 \sin \eta t), \end{equation*}
where $\eta = \sqrt{\omega_0^2 - c^2}\text{.}$ A particular solution is given by
\begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi), \end{equation*}
where the phase angle $\phi$ is
\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right) \end{equation*}
and the gain $G$ is given by
\begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}

### SubsectionExercises

###### Beats and Resonance

Solve the initial problems in Exercise Group 1–6. Do you find any beats or resonance?

###### 1

$y'' + y = \cos t\text{,}$ $y(0) = 0\text{,}$ $y'(0) = 0$

###### 2

$y'' + y = \cos(11t/10)\text{,}$ $y(0) = 0\text{,}$ $y'(0) = 0$

###### 3

$y'' + 4y = \sin 2 t\text{,}$ $y(0) = 1\text{,}$ $y'(0) = 0$

###### 4

$y'' + 4y = 3 \sin(13t/12)\text{,}$ $y(0) = 1\text{,}$ $y'(0) = -1$

###### 5

$y'' + 16y = 5\cos 4t\text{,}$ $y(0) = 1\text{,}$ $y'(0) = 0$

###### 6

$y'' + 16y = 5\cos(4t +1)\text{,}$ $y(0) = 0\text{,}$ $y'(0) = 0$

###### Plotting the Envelope for Beats

For each of the functions in Exercise Group 7–10 compute the envelope function. Then plot $y(t)$ and the envelope function on the same set of axes.

###### 7

$y = \cos 10 t - \cos 11t\text{.}$

###### 8

$y = \cos 11 t - \cos 12t\text{.}$

###### 9

$y = \sin 11 t - \sin 10t\text{.}$

###### 10

$y = \sin 12 t - \sin 11t\text{.}$

###### 11

Consider the function $y = \cos 10 t - \cos 11t\text{.}$

1. Find the mean frequency $\overline{\omega} = (\omega_0 + \omega)/2$ and half difference $\delta = (\omega_0 - \omega)/2\text{.}$
2. Use (a) to show that
\begin{equation*} y = \cos 10 t - \cos 11t = 2 \sin\left( \frac{1}{2} t \right) \sin(\left( \frac{21}{2} t \right). \end{equation*}
You will find the following trigonmetric identity very useful,
\begin{equation*} 2 \sin A \sin B = \cos(A - B) - \cos(A + B). \end{equation*}
3. Using (b), plot the graph of $y$ and $2 \sin(t/2)$ on the same axes.
###### 12

For large $t\text{,}$ every solution of

\begin{equation*} y'' + py' + qy = \cos \omega t \end{equation*}

oscillates with angular frequency $\omega$ and amplitude $A$ given by

\begin{equation*} A( \omega, p, q) = \frac{1}{\sqrt{(q - \omega^2)^2 + p^2 \omega^2}}. \end{equation*}

That is, the amplitude of $A$ is a function of the parameters $\omega\text{,}$ $p\text{,}$ and $q\text{.}$

1. Compute $\dfrac{\partial A}{\partial \omega}\text{.}$
2. For fixed $p$ and $q\text{,}$ let $M(p, q)$ denote the maximum value of $A( \omega, p, q)$ as a function of $\omega\text{.}$ Compute an expression for $M(p, q)\text{.}$
3. Set $q = 1$ and plot $M(p, q)$ as a function of $p\text{.}$
4. Explain why $M(p, q)$ is proportional to $1/p$ as $p \to 0\text{.}$
###### 13

Let us examine the amplitude and phase of the steady-state solution,

\begin{equation*} x_p(t) = G(\omega) A \cos(\omega t - \phi), \end{equation*}

where

\begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}

Now let $s = \omega / \omega_0$ and $D = 2c/\omega_0\text{.}$ These new constants, $s$ and $D\text{,}$ measure the ratio of the driving frequency to the natural frequency and the effect of the damping force, respectively. Thus,

\begin{equation*} G = \frac{1}{\omega_0^2\sqrt{(1- s^2)^2 + D^2 s^2}} \end{equation*}

or

\begin{equation*} \omega_0^2G = \frac{1}{\sqrt{(1- s^2)^2 + D^2 s^2}}. \end{equation*}

This expression shows us how the gain varies as $s = \omega/\omega_0$ varies. The natural frequency is fixed in

\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}

and $D = 2c/\omega_0$ is proportional to the damping constant.

### Subsection4.4.5Project—The Tacoma Narrows Bridge

#### SubsubsectionThe History of the Tacoma Narrows Bridge

On July 1, 1940, the Tacoma Narrows Bridge was completed and opened to traffic (Figure 4.4.10). From the first day of operation, the bridge experienced vertical oscillations of several feet. The bridge received the nickname “Galloping Gertie” for its wild behavior. On the morning of November 7, the bridge began undulating up and down with motions as large as three feet. Later, the bridge began to experience torsional oscillations as large as 23 degrees. The bridge finally came crashing down at 11:10 a.m. For a long time, the collapse of the Tacoma Narrows Bridge was attributed to resonance. However, we shall soon see the the story is more complicated.

Construction on first Tacoma Narrows Bridge, A.K.A. “Galloping Gertie”, was completed on July 1, 1940, the bridge was opened to traffic (Figure 4.4.11). Even on the first day of operation, people driving bridge experienced vertical oscillations of several feet. On the morning of November 7, these vertical oscillations were as large as three feet. Later that morning, the bridge began to experience torsional oscillations as great as 23 degrees. The bridge failed castrophically at 11:10 a.m. with the center span falling into Puget Sound below.

There are many interesting stories connected with the first Tacoma Narrows bridge. One of the insurance policies had been written by a local travel agent who pocketed the premium and failed to report the \$800,000 policy to his company. During sentencing at his trial, he pointed out that bridge officials had planned to cancel all policies in one more week after which his crime would have never been discovered. After the collapse of the Tacoma Narrows bridge, the governor of the State of Washington announced, “We are going to build the exact same bridge, exactly as before.” Theodore Von Kármán, director of the Guggenheim Aeronautical Laboratory at the California Institute of Technology and a member of the board of inquiry into the collapse, 18 The Jet Propulsion Laboratory, managed by the California Institute of Technology for the National Aeronautics and Space Administration grew out of the Guggenheim Aeronautical Laboratory. sent a telegram to the governor stating, “If you build the exact same bridge exactly as before, it will fall into the exact same river exactly as before.”

#### SubsubsectionThe Design of the Narrows Bridge

The Tacoma Narrows bridge is a suspension bridge like the Golden Gate bridge. The road bed of the suspension bridge was hung from vertical cables attached to other cables strung between the two towers. If we think of the cables as long springs, it is tempting to model the oscillations of the road bed with an harmonic oscillator, and one's first guess as to the reason for the collapse of the bridge would be resonance.

However, the answer may not be so simple. If the collapse of the bridge was due to resonance, the forcing frequency of a forced harmonic oscillator must be close to its natural frequency, but the wind (the forcing) simply does not behave like this. One explanation might be the fact that cables are not true springs. They act more like rubber bands. Imagine a mass suspended by both a spring and a rubber band. The rubber band acts like a spring when it is stretched, but there is no restoring force if the oscillator is in a compressed position (Figure 4.4.12). 19 For the arguments for and against resonance, consult the references at the end of this section.

#### SubsubsectionThe Spring-Rubber Band Model

Suppose that $y(t)$ is vertical displacement of the mass the spring-rubber band model (Figure 4.4.12), where the displacement is positive in the downward direction and negative in the upward direction. We will also assume that we have a damped system. If this is the case, we can model the motion of the mass with the differential equation

$$my'' = - b y' - k_1 y - k_2 y^+ + F(t),\label{equation-secondorder06-spring-rubber-band-1}\tag{4.4.18}$$

where $y^+ = \max(y, 0)\text{.}$ Here, $m$ is the mass, $b$ is the damping coefficient, $k_1$ is the spring constant, and $F(t)$ is the forcing term. The term $k_2 y^+$ corresponds to the rubber band in our system and can be explained as follows:

\begin{equation*} y^+(t) = \begin{cases} y(t), & \text{the rubber band is stretched}, \\ 0, & \text{the rubber band is slack or tight but not stretched}. \end{cases} \end{equation*}

Thus, the rubber band acts as a spring with spring coefficient $k_2\text{,}$ if the rubber band is stretched, but has no effect on the system if the rubber band is slack. Of course, the term $my''$ is just an application of Newton's Second law of Motion.

We can rewrite equation (4.4.18) in the more familiar form

$$my'' + b y' + k_1 y + k_2 y^+ = F(t),\label{equation-secondorder06-spring-rubber-band-2}\tag{4.4.19}$$

Our first observation is that this equation is a nonlinear due to the $y^+$ term. This may effect how sensitive the system is to initial conditions. A linear harmonic oscillator is not sensitive to initial conditions. For example, suppose we consider the harmonic oscillator

\begin{equation*} y'' + 2 y' + 2y = \sin t \end{equation*}

and examine the behavior of the solutions for different initial conditions. In Figure 4.4.13, we give the solutions for $y(0) = 0$ with $y'(0) = -4\text{,}$ $y'(0) = 1\text{,}$ and $y'(0) = 3\text{.}$ Even though the transient part of the solution can differ a great deal depending on the initial conditions, all of the solutions eventually approach the steady-state solution.

Let us return to our spring-rubber band model and let us assume that $m = 1\text{.}$ We can rewrite equation (4.4.19) as a system of first-order equations,

\begin{align*} y' & = v\\ v' & = - k_1 y - k_2 y^+ - b v + F(t) \end{align*}

or

\begin{align*} x' & = v\\ v' & = - h(y) - b v + F(t), \end{align*}

where

\begin{equation*} h(y) = \begin{cases} k_1 y & \text{if } y \lt 0 \\ (k_1 + k_2) y & \text{if } y \geq 0. \end{cases} \end{equation*}

Let us assume that the forcing term is $F(t) = 10 + \lambda \sin \mu t\text{.}$ If we keep $\lambda$ small and assume that both $y(0)$ and $y'(0)$ are small, there will always be tension on the rubber band. Or in the case of the Tacoma Narrows Bridge, there will always be tension on the cables of the suspension bridge. However, a large initial velocity such as a gust of wind on the bridge, might cause the cables to go slack.

To see how this works in practice, we choose constants $\lambda = 0.1\text{,}$ $\mu = 4\text{,}$ $k_1 = 13\text{,}$ $k_2 = 4\text{,}$ and $b = 0.01\text{,}$ the system becomes

\begin{align*} y' & = v\\ v' & = - h(y) - 0.01v + 10 + 4 \sin 0.1 t \end{align*}

where

\begin{equation*} h(y) = \begin{cases} 13 x & \text{if } x \lt 0 \\ 17 x & \text{if } x \geq 0. \end{cases} \end{equation*}

This is a nonautonomous (and nonlinear) system and can be very sensitive to initial conditions. If we find a numerical solution to this system with $y(0) = 0$ and $v(0) = 4.764\text{,}$ we obtain low amplitude solutions (Figure 4.4.14). However, we only have to change the initial velocity of the system slightly to $v(0) = 4.765$ in order to obtain large amplitude solutions that do not die out over time.

#### SubsubsectionThe Bridge Model

We can now model the bridge by

\begin{equation*} m y'' + \alpha y' + \beta y + \gamma y^+ = -gm + F(t), \end{equation*}

where

• $y$ is the height of the road bed at the center of the bridge.
• $m$ is the mass of the bridge.
• $\alpha y'$ is the damping term. We assume that $\alpha$ is small since suspension bridges are relatively flexible structures.
• $\beta y$ is the force provided by the material of the bridge pulling the bridge back to $y = 0\text{.}$
• $\gamma y^+$ is the force provided by the pull of the cables.
• $-gm$ is the force provided by gravity.
• $F(t)$ is the periodic force provided by the wind. Yes, wind does provide a periodic force. Just think of how a flag flutters in a wind.

#### SubsubsectionThe Model for Torsion

We have still not accounted for the large torsional oscillations that occurred on November 7, 1940, the day that the Tacoma Narrows bridge collapsed. If we look down the bridge (Figure 4.4.11), we see the road bed suspended from two cables—one on each side of the road. Thus, we can view a cross-section of the bridge as a rod suspended by cables or springs on each end. The bar is free to move vertically and rotate about its center of mass.

If a spring with spring constant $k$ is extended by a distance $y\text{,}$ the potential energy is $ky^2/2\text{.}$ Let $\theta$ be the angle of the rod from the horizontal. If a rod of mass $m$ and length $2l$ rotates about its center of gravity with angular velocity $d\theta/dt\text{,}$ then its kinetic energy is given by 20 Consult your local physicist.

\begin{equation*} \frac{1}{6} m l^2 \left( \frac{d \theta}{dt} \right)^2. \end{equation*}

The potential due to gravity is $-mgy\text{.}$

The extension of one spring is given by $(y - l \sin \theta)^+$ and $(y + l \sin \theta)^+$ in the other, where $y^+ = \max(y,0)\text{.}$ We are making the assumption that the cables remain vertical. Since the cables are relatively long compared to the length of the roadbed, this is a reasonable assumption.

The total potential energy is given by

\begin{equation*} V = \left(\frac{k}{2}\right) \left( \left( (y - l \sin \theta)^+\right)^2 + (y + l \sin \theta)^+ \right)^2 - mgy, \end{equation*}

and the total kinetic energy is given by

\begin{equation*} T = \frac{m}{2} \left( \frac{dy}{dt}\right)^2 + \frac{1}{6} m l^2 \left( \frac{d \theta}{dt} \right)^2. \end{equation*}

Let $L = T - V$ and put

\begin{equation*} \frac{d}{dt} \left( \frac{\delta L}{\delta \theta'} \right) = \frac{\delta L}{\delta \theta} \end{equation*}

and

\begin{equation*} \frac{d}{dt} \left( \frac{\delta L}{\delta y'} \right) = \frac{\delta L}{\delta y} \end{equation*}

where $\delta \theta'$ and $\delta y'$ are small damping terms. Thus, we obtain the equations

\begin{align*} \frac{1}{3} m l^2 \theta'' & = k l \cos \theta [(y - l \sin \theta)^+ - (y + l \sin \theta)^+]\\ y'' & = -k [(y - l \sin \theta)^+ + (y + l \sin \theta)^+] + mg \end{align*}

Since the force exerted by the springs (cables) is not perpendicular to the rod but at an angle $\theta\text{,}$ $\cos \theta$ must appear in the first equation. Adding forcing and damping terms, we have

\begin{align*} \theta'' & = - \delta \theta' + \frac{3k}{ml}\cos \theta [(y - l \sin \theta)^+ - (y + l \sin \theta)^+] + F(t)\\ y'' & = - \delta y' - \frac{k}{m}[(y - l \sin \theta)^+ + (y + l \sin \theta)^+] + g. \end{align*}

If we assume that the cables never lose tension, then

\begin{align*} (y - l \sin \theta)^+ & = y - l \sin \theta\\ (y + l \sin \theta)^+ & = y + l \sin \theta \end{align*}

and we end up with the uncoupled equations

\begin{align*} \theta'' & = - \delta \theta' - \frac{6k}{m}\cos \theta \sin \theta + F(t)\\ y'' & = - \delta y' - \frac{2k}{m} y + g. \end{align*}

The first equation models the torsional motion while the second models the vertical motion.

If the torsional oscillations are very small, then we can assume that

\begin{equation*} \sin \theta \approx \theta \qquad \text{and} \qquad \cos \theta \approx 1, \end{equation*}

or

\begin{align*} \theta'' & = - \delta \theta' - \frac{6k}{m}\theta + F(t)\\ y'' & = - \delta y' - \frac{2k}{m} y + g, \end{align*}

each of which is which is linear. However, the torsional oscillations for the Narrow bridge were as large as 23 degrees. The linear model

\begin{equation*} \theta'' = - 0.01 \theta' - (2.4) \theta + (1.2) \sin(0.06 t) \end{equation*}

and the nonlinear model

\begin{equation*} \theta'' = - 0.01 \theta' - (2.4) \cos \theta \sin \theta + (1.2) \sin(0.06 t) \end{equation*}

behave very differently under the initial conditions $\theta(0) = 1.2$ and $\theta'(0) = 0\text{.}$ See Figure 4.4.15 and Figure 4.4.16