### Learning Objectives

- 7.3.1. Solve a second-order differential equation representing simple harmonic motion.
- 7.3.2. Solve a second-order differential equation representing damped simple harmonic motion.
- 7.3.3. Solve a second-order differential equation representing forced simple harmonic motion.
- 7.3.4. Solve a second-order differential equation representing charge and current in an RLC series circuit.

We saw in the chapter introduction that second-order linear differential equations are used to model many situations in physics and engineering. In this section, we look at how this works for systems of an object with mass attached to a vertical spring and an electric circuit containing a resistor, an inductor, and a capacitor connected in series. Models such as these can be used to approximate other more complicated situations; for example, bonds between atoms or molecules are often modeled as springs that vibrate, as described by these same differential equations.

### Simple Harmonic Motion

Consider a mass suspended from a spring attached to a rigid support. (This is commonly called a spring-mass system.) Gravity is pulling the mass downward and the restoring force of the spring is pulling the mass upward. As shown in Figure 7.2, when these two forces are equal, the mass is said to be at the equilibrium position. If the mass is displaced from equilibrium, it oscillates up and down. This behavior can be modeled by a second-order constant-coefficient differential equation.

Let $x(t)$ denote the displacement of the mass from equilibrium. Note that for spring-mass systems of this type, it is customary to adopt the convention that down is positive. Thus, a positive displacement indicates the mass is *below* the equilibrium point, whereas a negative displacement indicates the mass is *above* equilibrium. Displacement is usually given in feet in the English system or meters in the metric system.

Consider the forces acting on the mass. The force of gravity is given by $mg\text{.}$ In the English system, mass is in slugs and the acceleration resulting from gravity is in feet per second squared. The acceleration resulting from gravity is constant, so in the English system, $g=32$ ft/sec^{2}. Recall that 1 slug-foot/sec^{2} is a pound, so the expression *mg* can be expressed in pounds. Metric system units are kilograms for mass and m/sec^{2} for gravitational acceleration. In the metric system, we have $g=9.8$ m/sec^{2}.

According to Hooke’s law, the restoring force of the spring is proportional to the displacement and acts in the opposite direction from the displacement, so the restoring force is given by $\text{\u2212}k\left(s+x\right).$ The spring constant is given in pounds per foot in the English system and in newtons per meter in the metric system.

Now, by Newton’s second law, the sum of the forces on the system (gravity plus the restoring force) is equal to mass times acceleration, so we have

However, by the way we have defined our equilibrium position, $mg=ks,$ the differential equation becomes

It is convenient to rearrange this equation and introduce a new variable, called the angular frequency, $\omega \text{.}$ Letting $\omega =\sqrt{k\text{/}m},$ we can write the equation as

This differential equation has the general solution

which gives the position of the mass at any point in time. The motion of the mass is called simple harmonic motion. The period of this motion (the time it takes to complete one oscillation) is $T=\frac{2\pi}{\omega}$ and the frequency is $f=\frac{1}{T}=\frac{\omega}{2\pi}$ (Figure 7.3).

### Example 7.17

#### Simple Harmonic Motion

Assume an object weighing 2 lb stretches a spring 6 in. Find the equation of motion if the spring is released from the equilibrium position with an upward velocity of 16 ft/sec. What is the period of the motion?

#### Solution

We first need to find the spring constant. We have

We also know that weight *W* equals the product of mass *m* and the acceleration due to gravity *g*. In English units, the acceleration due to gravity is 32 ft/sec^{2}.

Thus, the differential equation representing this system is

Multiplying through by 16, we get $x\text{\u2033}+64x=0,$ which can also be written in the form $x\text{\u2033}+({8}^{2})x=0.$ This equation has the general solution

The mass was released from the equilibrium position, so $x(0)=0,$ and it had an initial upward velocity of 16 ft/sec, so ${x}^{\prime}(0)=\mathrm{-16}.$ Applying these initial conditions to solve for ${c}_{1}$ and ${c}_{2}.$ gives

The period of this motion is $\frac{2\pi}{8}=\frac{\pi}{4}$ sec.

A 200-g mass stretches a spring 5 cm. Find the equation of motion of the mass if it is released from rest from a position 10 cm below the equilibrium position. What is the frequency of this motion?

Writing the general solution in the form $x(t)={c}_{1}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(\omega t\right)+{c}_{2}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(\omega t\right)$ has some advantages. It is easy to see the link between the differential equation and the solution, and the period and frequency of motion are evident. This form of the function tells us very little about the amplitude of the motion, however. In some situations, we may prefer to write the solution in the form

Although the link to the differential equation is not as explicit in this case, the period and frequency of motion are still evident. Furthermore, the amplitude of the motion, *A*, is obvious in this form of the function. The constant $\varphi $ is called a *phase shift* and has the effect of shifting the graph of the function to the left or right.

To convert the solution to this form, we want to find the values of *A* and $\varphi $ such that

We first apply the trigonometric identity

to get

Thus,

If we square both of these equations and add them together, we get

Thus,

Now, to find $\varphi ,$ go back to the equations for ${c}_{1}$ and ${c}_{2},$ but this time, divide the first equation by the second equation to get

Then,

We summarize this finding in the following theorem.

#### Solution to the Equation for Simple Harmonic Motion

The function $x(t)={c}_{1}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(\omega t\right)+{c}_{2}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(\omega t\right)$ can be written in the form $x(t)=A\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(\omega t+\varphi \right),$ where $A=\sqrt{{c}_{1}^{2}+{c}_{2}^{2}}$ and $\text{tan}\phantom{\rule{0.1em}{0ex}}\varphi =\frac{{c}_{1}}{{c}_{2}}.$

Note that when using the formula $\text{tan}\phantom{\rule{0.1em}{0ex}}\varphi =\frac{{c}_{1}}{{c}_{2}}$ to find $\varphi ,$ we must take care to ensure $\varphi $ is in the right quadrant (Figure 7.4).

### Example 7.18

#### Expressing the Solution with a Phase Shift

Express the following functions in the form $A\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(\omega t+\varphi \right).$ What is the frequency of motion? The amplitude?

- $x(t)=2\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(3t\right)+\text{sin}\phantom{\rule{0.1em}{0ex}}\left(3t\right)$
- $x(t)=3\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(2t\right)-2\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(2t\right)$

#### Solution

- We have

$$A=\sqrt{{c}_{1}^{2}+{c}_{2}^{2}}=\sqrt{{2}^{2}+{1}^{2}}=\sqrt{5}$$

and

$$\text{tan}\phantom{\rule{0.1em}{0ex}}\varphi =\frac{{c}_{1}}{{c}_{2}}=\frac{2}{1}=2.$$

Note that both ${c}_{1}$ and ${c}_{2}$ are positive, so $\varphi $ is in the first quadrant. Thus,

$$\varphi \approx 1.107\phantom{\rule{0.2em}{0ex}}\text{rad,}$$

so we have

$$x(t)=2\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(3t\right)+\text{sin}\phantom{\rule{0.1em}{0ex}}\left(3t\right)=\sqrt{5}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(3t+1.107\right).$$

The frequency is $\frac{\omega}{2\pi}=\frac{3}{2\pi}\approx 0.477.$ The amplitude is $\sqrt{5}.$ - We have

$$A=\sqrt{{c}_{1}^{2}+{c}_{2}^{2}}=\sqrt{{3}^{2}+{2}^{2}}=\sqrt{13}$$

and

$$\text{tan}\phantom{\rule{0.1em}{0ex}}\varphi =\frac{{c}_{1}}{{c}_{2}}=\frac{3}{\mathrm{-2}}=-\frac{3}{2}.$$

Note that ${c}_{1}$ is positive but ${c}_{2}$ is negative, so $\varphi $ is in the fourth quadrant. Thus,

$$\varphi \approx -0.983\phantom{\rule{0.2em}{0ex}}\text{rad,}$$

so we have

$$\begin{array}{cc}\hfill x(t)& =3\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(2t\right)-2\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(2t\right)\hfill \\ & =\sqrt{13}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(2t-0.983\right).\hfill \end{array}$$

The frequency is $\frac{\omega}{2\pi}=\frac{2}{2\pi}\approx 0.318.$ The amplitude is $\sqrt{13}.$

Express the function $x(t)=\text{cos}\phantom{\rule{0.1em}{0ex}}\left(4t\right)+4\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(4t\right)$ in the form $A\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(\omega t+\varphi \right).$ What is the frequency of motion? The amplitude?

### Damped Vibrations

With the model just described, the motion of the mass continues indefinitely. Clearly, this doesn’t happen in the real world. In the real world, there is almost always some friction in the system, which causes the oscillations to die off slowly—an effect called *damping*. So now let’s look at how to incorporate that damping force into our differential equation.

Physical spring-mass systems almost always have some damping as a result of friction, air resistance, or a physical damper, called a *dashpot* (a pneumatic cylinder; see Figure 7.5).

Because damping is primarily a friction force, we assume it is proportional to the velocity of the mass and acts in the opposite direction. So the damping force is given by $\text{\u2212}b{x}^{\prime}$ for some constant $b>0.$ Again applying Newton’s second law, the differential equation becomes

Then the associated characteristic equation is

Applying the quadratic formula, we have

Just as in Second-Order Linear Equations we consider three cases, based on whether the characteristic equation has distinct real roots, a repeated real root, or complex conjugate roots.

#### Case 1: ${b}^{2}>4mk$

In this case, we say the system is *overdamped*. The general solution has the form

where both ${\lambda}_{1}$ and ${\lambda}_{2}$ are less than zero. Because the exponents are negative, the displacement decays to zero over time, usually quite quickly. Overdamped systems do not oscillate (no more than one change of direction), but simply move back toward the equilibrium position. Figure 7.6 shows what typical critically damped behavior looks like.

### Example 7.19

#### Overdamped Spring-Mass System

A 16-lb mass is attached to a 10-ft spring. When the mass comes to rest in the equilibrium position, the spring measures 15 ft 4 in. The system is immersed in a medium that imparts a damping force equal to $\frac{5}{2}$ times the instantaneous velocity of the mass. Find the equation of motion if the mass is pushed upward from the equilibrium position with an initial upward velocity of 5 ft/sec. What is the position of the mass after 10 sec? Its velocity?

#### Solution

The mass stretches the spring 5 ft 4 in., or $\frac{16}{3}$ ft. Thus, $16=\left(\frac{16}{3}\right)k,$ so $k=3.$ We also have $m=\frac{16}{32}=\frac{1}{2},$ so the differential equation is

Multiplying through by 2 gives $x\text{\u2033}+5{x}^{\prime}+6x=0,$ which has the general solution

Applying the initial conditions, $x(0)=0$ and ${x}^{\prime}(0)=\mathrm{-5},$ we get

After 10 sec the mass is at position

so it is, effectively, at the equilibrium position. We have ${x}^{\prime}(t)=10{e}^{\mathrm{-2}t}-15{e}^{\mathrm{-3}t},$ so after 10 sec the mass is moving at a velocity of

After only 10 sec, the mass is barely moving.

### Checkpoint 7.17

A 2-kg mass is attached to a spring with spring constant 24 N/m. The system is then immersed in a medium imparting a damping force equal to 16 times the instantaneous velocity of the mass. Find the equation of motion if it is released from rest at a point 40 cm below equilibrium.

#### Case 2: ${b}^{2}=4mk$

In this case, we say the system is *critically damped*. The general solution has the form

where ${\lambda}_{1}$ is less than zero. The motion of a critically damped system is very similar to that of an overdamped system. It does not oscillate. However, with a critically damped system, if the damping is reduced even a little, oscillatory behavior results. From a practical perspective, physical systems are almost always either overdamped or underdamped (case 3, which we consider next). It is impossible to fine-tune the characteristics of a physical system so that ${b}^{2}$ and $4mk$ are exactly equal. Figure 7.7 shows what typical critically damped behavior looks like.

### Example 7.20

#### Critically Damped Spring-Mass System

A 1-kg mass stretches a spring 20 cm. The system is attached to a dashpot that imparts a damping force equal to 14 times the instantaneous velocity of the mass. Find the equation of motion if the mass is released from equilibrium with an upward velocity of 3 m/sec.

#### Solution

We have $mg=1\left(9.8\right)=0.2k,$ so $k=49.$ Then, the differential equation is

which has general solution

Applying the initial conditions $x(0)=0$ and ${x}^{\prime}(0)=\mathrm{-3}$ gives

### Checkpoint 7.18

A 1-lb weight stretches a spring 6 in., and the system is attached to a dashpot that imparts a damping force equal to half the instantaneous velocity of the mass. Find the equation of motion if the mass is released from rest at a point 6 in. below equilibrium.

#### Case 3: ${b}^{2}<4mk$

In this case, we say the system is *underdamped*. The general solution has the form

where $\alpha $ is less than zero. Underdamped systems do oscillate because of the sine and cosine terms in the solution. However, the exponential term dominates eventually, so the amplitude of the oscillations decreases over time. Figure 7.8 shows what typical underdamped behavior looks like.

Note that for all damped systems, $\underset{t\to \infty}{\text{lim}}x(t)=0.$ The system always approaches the equilibrium position over time.

### Example 7.21

#### Underdamped Spring-Mass System

A 16-lb weight stretches a spring 3.2 ft. Assume the damping force on the system is equal to the instantaneous velocity of the mass. Find the equation of motion if the mass is released from rest at a point 9 in. below equilibrium.

#### Solution

We have $k=\frac{16}{3.2}=5$ and $m=\frac{16}{32}=\frac{1}{2},$ so the differential equation is

This equation has the general solution

Applying the initial conditions, $x(0)=\frac{3}{4}$ and ${x}^{\prime}(0)=0,$ we get

### Checkpoint 7.19

A 1-kg mass stretches a spring 49 cm. The system is immersed in a medium that imparts a damping force equal to four times the instantaneous velocity of the mass. Find the equation of motion if the mass is released from rest at a point 24 cm above equilibrium.

### Example 7.22

#### Chapter Opener: Modeling a Motorcycle Suspension System

For motocross riders, the suspension systems on their motorcycles are very important. The off-road courses on which they ride often include jumps, and losing control of the motorcycle when they land could cost them the race.

This suspension system can be modeled as a damped spring-mass system. We define our frame of reference with respect to the frame of the motorcycle. Assume the end of the shock absorber attached to the motorcycle frame is fixed. Then, the “mass” in our spring-mass system is the motorcycle wheel. We measure the position of the wheel with respect to the motorcycle frame. This may seem counterintuitive, since, in many cases, it is actually the motorcycle frame that moves, but this frame of reference preserves the development of the differential equation that was done earlier. As with earlier development, we define the downward direction to be positive.

When the motorcycle is lifted by its frame, the wheel hangs freely and the spring is uncompressed. This is the spring’s natural position. When the motorcycle is placed on the ground and the rider mounts the motorcycle, the spring compresses and the system is in the equilibrium position (Figure 7.10).

This system can be modeled using the same differential equation we used before:

A motocross motorcycle weighs 204 lb, and we assume a rider weight of 180 lb. When the rider mounts the motorcycle, the suspension compresses 4 in., then comes to rest at equilibrium. The suspension system provides damping equal to 240 times the instantaneous vertical velocity of the motorcycle (and rider).

- Set up the differential equation that models the behavior of the motorcycle suspension system.
- We are interested in what happens when the motorcycle lands after taking a jump. Let time $t=0$ denote the time when the motorcycle first contacts the ground. If the motorcycle hits the ground with a velocity of 10 ft/sec downward, find the equation of motion of the motorcycle after the jump.
- Graph the equation of motion over the first second after the motorcycle hits the ground.

#### Solution

- We have defined equilibrium to be the point where $mg=ks,$ so we have

$$\begin{array}{ccc}\hfill mg& =\hfill & ks\hfill \\ \hfill 384& =\hfill & k\left(\frac{1}{3}\right)\hfill \\ \hfill k& =\hfill & 1152.\hfill \end{array}$$

We also have

$$\begin{array}{ccc}\hfill W& =\hfill & mg\hfill \\ \hfill 384& =\hfill & m\left(32\right)\hfill \\ \hfill m& =\hfill & 12.\hfill \end{array}$$

Therefore, the differential equation that models the behavior of the motorcycle suspension is

$$12x\text{\u2033}+240{x}^{\prime}+1152x=0.$$

Dividing through by 12, we get

$$x\text{\u2033}+20{x}^{\prime}+96x=0.$$ - The differential equation found in part a. has the general solution

$$x(t)={c}_{1}{e}^{\mathrm{-8}t}+{c}_{2}{e}^{\mathrm{-12}t}.$$

Now, to determine our initial conditions, we consider the position and velocity of the motorcycle wheel when the wheel first contacts the ground. Since the motorcycle was in the air prior to contacting the ground, the wheel was hanging freely and the spring was uncompressed. Therefore the wheel is 4 in. $\left(\frac{1}{3}\phantom{\rule{0.2em}{0ex}}\text{ft}\right)$ below the equilibrium position (with respect to the motorcycle frame), and we have $x(0)=\frac{1}{3}.$ According to the problem statement, the motorcycle has a velocity of 10 ft/sec downward when the motorcycle contacts the ground, so ${x}^{\prime}(0)=10.$ Applying these initial conditions, we get ${c}_{1}=\frac{7}{2}$ and ${c}_{2}=\text{\u2212}\left(\frac{19}{6}\right),$ so the equation of motion is

$$x(t)=\frac{7}{2}{e}^{\mathrm{-8}t}-\frac{19}{6}{e}^{\mathrm{-12}t}.$$ - The graph is shown in Figure 7.11.

### Student Project

#### Landing Vehicle

NASA is planning a mission to Mars. To save money, engineers have decided to adapt one of the moon landing vehicles for the new mission. However, they are concerned about how the different gravitational forces will affect the suspension system that cushions the craft when it touches down. The acceleration resulting from gravity on the moon is 1.6 m/sec^{2}, whereas on Mars it is 3.7 m/sec^{2}.

The suspension system on the craft can be modeled as a damped spring-mass system. In this case, the spring is below the moon lander, so the spring is slightly compressed at equilibrium, as shown in Figure 7.12.

We retain the convention that down is positive. Despite the new orientation, an examination of the forces affecting the lander shows that the same differential equation can be used to model the position of the landing craft relative to equilibrium:

where *m* is the mass of the lander, *b* is the damping coefficient, and *k* is the spring constant.

- The lander has a mass of 15,000 kg and the spring is 2 m long when uncompressed. The lander is designed to compress the spring 0.5 m to reach the equilibrium position under lunar gravity. The dashpot imparts a damping force equal to 48,000 times the instantaneous velocity of the lander. Set up the differential equation that models the motion of the lander when the craft lands on the moon.
- Let time $t=0$ denote the instant the lander touches down. The rate of descent of the lander can be controlled by the crew, so that it is descending at a rate of 2 m/sec when it touches down. Find the equation of motion of the lander on the moon.
- If the lander is traveling too fast when it touches down, it could fully compress the spring and “bottom out.” Bottoming out could damage the landing craft and must be avoided at all costs. Graph the equation of motion found in part 2. If the spring is 0.5 m long when fully compressed, will the lander be in danger of bottoming out?
- Assuming NASA engineers make no adjustments to the spring or the damper, how far does the lander compress the spring to reach the equilibrium position under Martian gravity?
- If the lander crew uses the same procedures on Mars as on the moon, and keeps the rate of descent to 2 m/sec, will the lander bottom out when it lands on Mars?
- What adjustments, if any, should the NASA engineers make to use the lander safely on Mars?

### Forced Vibrations

The last case we consider is when an external force acts on the system. In the case of the motorcycle suspension system, for example, the bumps in the road act as an external force acting on the system. Another example is a spring hanging from a support; if the support is set in motion, that motion would be considered an external force on the system. We model these forced systems with the nonhomogeneous differential equation

where the external force is represented by the $f(t)$ term. As we saw in Nonhomogenous Linear Equations, differential equations such as this have solutions of the form

where ${c}_{1}{x}_{1}(t)+{c}_{2}{x}_{2}(t)$ is the general solution to the complementary equation and ${x}_{p}(t)$ is a particular solution to the nonhomogeneous equation. If the system is damped, $\underset{t\to \infty}{\text{lim}}{c}_{1}{x}_{1}(t)+{c}_{2}{x}_{2}(t)=0.$ Since these terms do not affect the long-term behavior of the system, we call this part of the solution the *transient solution***.** The long-term behavior of the system is determined by ${x}_{p}(t),$ so we call this part of the solution the steady-state solution.

### Media

This website shows a simulation of forced vibrations.

### Example 7.23

#### Forced Vibrations

A mass of 1 slug stretches a spring 2 ft and comes to rest at equilibrium. The system is attached to a dashpot that imparts a damping force equal to eight times the instantaneous velocity of the mass. Find the equation of motion if an external force equal to $f(t)=8\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(4t\right)$ is applied to the system beginning at time $t=0.$ What is the transient solution? What is the steady-state solution?

#### Solution

We have $mg=1\left(32\right)=2k,$ so $k=16$ and the differential equation is

The general solution to the complementary equation is

Assuming a particular solution of the form ${x}_{p}(t)=A\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(4t\right)+B\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(4t\right)$ and using the method of undetermined coefficients, we find ${x}_{p}(t)=-\frac{1}{4}\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(4t\right),$ so

At $t=0,$ the mass is at rest in the equilibrium position, so $x(0)={x}^{\prime}(0)=0.$ Applying these initial conditions to solve for ${c}_{1}$ and ${c}_{2},$ we get

The transient solution is $\frac{1}{4}{e}^{\mathrm{-4}t}+t{e}^{\mathrm{-4}t}.$ The steady-state solution is $-\frac{1}{4}\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(4t\right).$

A mass of 2 kg is attached to a spring with constant 32 N/m and comes to rest in the equilibrium position. Beginning at time $t=0,$ an external force equal to$f(t)=68{e}^{\mathrm{-2}t}\text{cos}\phantom{\rule{0.1em}{0ex}}\left(4t\right)$ is applied to the system. Find the equation of motion if there is no damping. What is the transient solution? What is the steady-state solution?

### Student Project

#### Resonance

Consider an undamped system exhibiting simple harmonic motion. In the real world, we never truly have an undamped system; –some damping always occurs. For theoretical purposes, however, we could imagine a spring-mass system contained in a vacuum chamber. With no air resistance, the mass would continue to move up and down indefinitely.

The frequency of the resulting motion, given by $f=\frac{1}{T}=\frac{\omega}{2\pi},$ is called the *natural frequency of the system*. If an external force acting on the system has a frequency close to the natural frequency of the system, a phenomenon called *resonance* results. The external force reinforces and amplifies the natural motion of the system.

- Consider the differential equation $x\text{\u2033}+x=0.$ Find the general solution. What is the natural frequency of the system?
- Now suppose this system is subjected to an external force given by $f(t)=5\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}t\text{.}$ Solve the initial-value problem $x\text{\u2033}+x=5\phantom{\rule{0.1em}{0ex}}\text{cos}\phantom{\rule{0.1em}{0ex}}t,$ $x(0)=0,$ ${x}^{\prime}(0)=1.$
- Graph the solution. What happens to the behavior of the system over time?
- In the real world, there is always some damping. However, if the damping force is weak, and the external force is strong enough, real-world systems can still exhibit resonance. One of the most famous examples of resonance is the collapse of the Tacoma Narrows Bridge on November 7, 1940. The bridge had exhibited strange behavior ever since it was built. The roadway had a strange “bounce” to it. On the day it collapsed, a strong windstorm caused the roadway to twist and ripple violently. The bridge was unable to withstand these forces and it ultimately collapsed. Experts believe the windstorm exerted forces on the bridge that were very close to its natural frequency, and the resulting resonance ultimately shook the bridge apart.

### Media

This website contains more information about the collapse of the Tacoma Narrows Bridge.

### Media

During the short time the Tacoma Narrows Bridge stood, it became quite a tourist attraction. Several people were on site the day the bridge collapsed, and one of them caught the collapse on film. Watch the video to see the collapse.

- Another real-world example of resonance is a singer shattering a crystal wineglass when she sings just the right note. When someone taps a crystal wineglass or wets a finger and runs it around the rim, a tone can be heard. That note is created by the wineglass vibrating at its natural frequency. If a singer then sings that same note at a high enough volume, the glass shatters as a result of resonance.

### The *RLC* Series Circuit

Consider an electrical circuit containing a resistor, an inductor, and a capacitor, as shown in Figure 7.10. Such a circuit is called an *RLC* series circuit. *RLC* circuits are used in many electronic systems, most notably as tuners in AM/FM radios. The tuning knob varies the capacitance of the capacitor, which in turn tunes the radio. Such circuits can be modeled by second-order, constant-coefficient differential equations.

Let $I(t)$ denote the current in the *RLC* circuit and $q(t)$ denote the charge on the capacitor. Furthermore, let *L* denote inductance in henrys (H), *R* denote resistance in ohms $\left(\text{\Omega}\right),$ and *C* denote capacitance in farads (F). Last, let $E(t)$ denote electric potential in volts (V).

Kirchhoff’s voltage rule states that the sum of the voltage drops around any closed loop must be zero. So, we need to consider the voltage drops across the inductor (denoted ${E}_{L}$), the resistor (denoted ${E}_{R}$), and the capacitor (denoted ${E}_{C}$). Because the *RLC* circuit shown in Figure 7.10 includes a voltage source, $E(t),$ which adds voltage to the circuit, we have ${E}_{L}+{E}_{R}+{E}_{C}=E(t).$

We present the formulas below without further development. Those of you interested in the derivation of these formulas should consult a physics text. Using Faraday’s law and Lenz’s law, the voltage drop across an inductor can be shown to be proportional to the instantaneous rate of change of current, with proportionality constant *L*. Thus,

Next, according to Ohm’s law, the voltage drop across a resistor is proportional to the current passing through the resistor, with proportionality constant *R*. Therefore,

Last, the voltage drop across a capacitor is proportional to the charge, *q*, on the capacitor, with proportionality constant $1\text{/}C.$ Thus,

Adding these terms together, we get

Noting that $I=\left(dq\right)\text{/}\left(dt\right),$ this becomes

Mathematically, this system is analogous to the spring-mass systems we have been examining in this section.

### Example 7.24

#### The *RLC* Series Circuit

Find the charge on the capacitor in an *RLC* series circuit where $L=5\text{/}3$ H, $R=10\text{\Omega},$ $C=1\text{/}30$ F, and $E(t)=300$ V. Assume the initial charge on the capacitor is 0 C and the initial current is 9 A. What happens to the charge on the capacitor over time?

#### Solution

We have

The general solution to the complementary equation is

Assume a particular solution of the form ${q}_{p}=A,$ where $A$ is a constant. Using the method of undetermined coefficients, we find $A=10.$ So,

Applying the initial conditions $q(0)=0$ and $i(0)=\left(\left(dq\right)\text{/}\left(dt\right)\right)(0)=9,$ we find ${c}_{1}=\mathrm{-10}$ and ${c}_{2}=\mathrm{-7}.$ So the charge on the capacitor is

Looking closely at this function, we see the first two terms will decay over time (as a result of the negative exponent in the exponential function). Therefore, the capacitor eventually approaches a steady-state charge of 10 C.

Find the charge on the capacitor in an *RLC* series circuit where $L=1\text{/}5$ H, $R=2\text{/}5\text{\Omega},$ $C=1\text{/}2$ F, and $E(t)=50$ V. Assume the initial charge on the capacitor is 0 C and the initial current is 4 A.

### Section 7.3 Exercises

A mass weighing 4 lb stretches a spring 8 in. Find the equation of motion if the spring is released from the equilibrium position with a downward velocity of 12 ft/sec. What is the period and frequency of the motion?

A mass weighing 2 lb stretches a spring 2 ft. Find the equation of motion if the spring is released from 2 in. below the equilibrium position with an upward velocity of 8 ft/sec. What is the period and frequency of the motion?

A 100-g mass stretches a spring 0.1 m. Find the equation of motion of the mass if it is released from rest from a position 20 cm below the equilibrium position. What is the frequency of this motion?

A 400-g mass stretches a spring 5 cm. Find the equation of motion of the mass if it is released from rest from a position 15 cm below the equilibrium position. What is the frequency of this motion?

A block has a mass of 9 kg and is attached to a vertical spring with a spring constant of 0.25 N/m. The block is stretched 0.75 m below its equilibrium position and released.

- Find the position function $x(t)$ of the block.
- Find the period and frequency of the vibration.
- Sketch a graph of $x(t).$
- At what time does the block first pass through the equilibrium position?

A block has a mass of 5 kg and is attached to a vertical spring with a spring constant of 20 N/m. The block is released from the equilibrium position with a downward velocity of 10 m/sec.

- Find the position function $x(t)$ of the block.
- Find the period and frequency of the vibration.
- Sketch a graph of $x(t).$
- At what time does the block first pass through the equilibrium position?

A 1-kg mass is attached to a vertical spring with a spring constant of 21 N/m. The resistance in the spring-mass system is equal to 10 times the instantaneous velocity of the mass.

- Find the equation of motion if the mass is released from a position 2 m below its equilibrium position with a downward velocity of 2 m/sec.
- Graph the solution and determine whether the motion is overdamped, critically damped, or underdamped.

An 800-lb weight (25 slugs) is attached to a vertical spring with a spring constant of 226 lb/ft. The system is immersed in a medium that imparts a damping force equal to 10 times the instantaneous velocity of the mass.

- Find the equation of motion if it is released from a position 20 ft below its equilibrium position with a downward velocity of 41 ft/sec.
- Graph the solution and determine whether the motion is overdamped, critically damped, or underdamped.

A 9-kg mass is attached to a vertical spring with a spring constant of 16 N/m. The system is immersed in a medium that imparts a damping force equal to 24 times the instantaneous velocity of the mass.

- Find the equation of motion if it is released from its equilibrium position with an upward velocity of 4 m/sec.
- Graph the solution and determine whether the motion is overdamped, critically damped, or underdamped.

A 1-kg mass stretches a spring 6.25 cm. The resistance in the spring-mass system is equal to eight times the instantaneous velocity of the mass.

- Find the equation of motion if the mass is released from a position 5 m below its equilibrium position with an upward velocity of 10 m/sec.
- Determine whether the motion is overdamped, critically damped, or underdamped.

A 32-lb weight (1 slug) stretches a vertical spring 4 in. The resistance in the spring-mass system is equal to four times the instantaneous velocity of the mass.

- Find the equation of motion if it is released from its equilibrium position with a downward velocity of 12 ft/sec.
- Determine whether the motion is overdamped, critically damped, or underdamped.

A 64-lb weight is attached to a vertical spring with a spring constant of 4.625 lb/ft. The resistance in the spring-mass system is equal to the instantaneous velocity. The weight is set in motion from a position 1 ft below its equilibrium position with an upward velocity of 2 ft/sec. Is the mass above or below the equation position at the end of $\pi $ sec? By what distance?

A mass that weighs 8 lb stretches a spring 6 inches. The system is acted on by an external force of $8\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}8t$ lb. If the mass is pulled down 3 inches and then released, determine the position of the mass at any time.

A mass that weighs 6 lb stretches a spring 3 in. The system is acted on by an external force of $8\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}\left(4t\right)$ lb. If the mass is pulled down 1 inch and then released, determine the position of the mass at any time.

Find the charge on the capacitor in an *RLC* series circuit where $L=40$ H, $R=30\text{\Omega},$ $C=1\text{/}200$ F, and $E(t)=200$ V. Assume the initial charge on the capacitor is 7 C and the initial current is 0 A.

Find the charge on the capacitor in an *RLC* series circuit where $L=2$ H, $R=24\text{\Omega},$ $C=0.005$ F, and $E(t)=12\phantom{\rule{0.1em}{0ex}}\text{sin}\phantom{\rule{0.1em}{0ex}}10t$ V. Assume the initial charge on the capacitor is 0.001 C and the initial current is 0 A.

A series circuit consists of a device where $L=1$ H, $R=20\text{\Omega},$ $C=0.002$ F, and $E(t)=12$ V. If the initial charge and current are both zero, find the charge and current at time *t*.

A series circuit consists of a device where $L=\frac{1}{2}$ H, $R=10\text{\Omega},$ $C=\frac{1}{50}$ F, and $E(t)=250$ V. If the initial charge on the capacitor is 0 C and the initial current is 18 A, find the charge and current at time *t*.