Javascript required
Skip to content Skip to sidebar Skip to footer

Find the Steady-state Solution for Each Differential Equation.

17.3: Applications of Second-Order Differential Equations

  1. Last updated
  2. Save as PDF
  • Page ID
    2628
    • Professor (Mathematics) at Massachusetts Institute of Technology (Strang) & University of Wisconsin-Stevens Point (Herman)
    • Publisher: OpenStax CNX

    Learning Objectives

    • Solve a second-order differential equation representing simple harmonic motion.
    • Solve a second-order differential equation representing damped simple harmonic motion.
    • Solve a second-order differential equation representing forced simple harmonic motion.
    • 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 \(\PageIndex{1}\), 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.

    This figure has three images of springs. The first image is a vertical spring in its natural position with length L attached at the top to a fixed point. The second image shows a vertical spring with a mass m attached to the spring, stretching the spring distance s from L. The spring is in equilibrium. The third image is a vertical spring with mass m attached where the spring is in motion, distance x from equilibrium L + s.
    Figure \(\PageIndex{1}\): A spring in its natural position (a), at equilibrium with a mass m attached (b), and in oscillatory motion (c).

    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.mg. 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/sec2 is a pound, so the expression mg can be expressed in pounds. Metric system units are kilograms for mass and m/sec2 for gravitational acceleration. In the metric system, we have \(g=9.8\) m/sec2.

    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 \(−k(s+x).\) 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

    \[\begin{align*}mx″ &=−k(s+x)+mg \\[4pt] &=−ks−kx+mg. \end{align*}\]

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

    \[mx″+kx=0. \nonumber\]

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

    \[x''+ω^2x=0. \nonumber\]

    This differential equation has the general solution

    \[x(t)=c_1 \cos ωt+c_2 \sin ωt, \label{GeneralSol}\]

    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=\dfrac{2π}{ω}\) and the frequency is \(f=\dfrac{1}{T}=\dfrac{ω}{2π}\) (Figure \(\PageIndex{2}\)).

    This figure is the graph of f(t) = sin 2t. It is a periodic, oscillating graph. The period of the graph is represented with a line pointing from one peak to the next. It is labeled with the period  T = 2π/ω.
    Figure \(\PageIndex{2}\): A graph of vertical displacement versus time for simple harmonic motion.

    Example \(\PageIndex{1}\): 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

    \[\begin{align*}mg &=ks\\[4pt] 2 &=k \left(\dfrac{1}{2}\right)\\[4pt] k &=4. \end{align*}\]

    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/sec2.

    \[\begin{align*}W &=mg\\[4pt] 2 &=m(32)\\[4pt] m &=\dfrac{1}{16}\end{align*}\]

    Thus, the differential equation representing this system is

    \[\dfrac{1}{16}x″+4x=0. \nonumber\]

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

    \[x(t)=c_1 \cos (8t)+c_2 \sin (8t). \nonumber\]

    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′(0)=−16\). Applying these initial conditions to solve for \(c_1\) and \(c_2\). gives

    \[x(t)=−2 \sin 8t. \nonumber\]

    The period of this motion is \(\dfrac{2π}{8}=\dfrac{π}{4}\) sec.

    Exercise \(\PageIndex{1}\)

    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?

    Hint

    First find the spring constant.

    Answer

    \(x(t)=0.1 \cos (14t)\) (in meters); frequency is \(\dfrac{14}{2π}\) Hz.

    Writing the general solution in the form \(x(t)=c_1 \cos (ωt)+c_2 \sin(ωt)\) (Equation \ref{GeneralSol}) 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

    \[x(t)=A \sin (ωt+ϕ). \nonumber\]

    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 \(ϕ\) 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 \(ϕ\) such that

    \[c_1 \cos (ωt)+c_2 \sin (ωt)=A \sin (ωt+ϕ). \nonumber\]

    We first apply the trigonometric identity

    \[\sin (α+β)= \sin α \cos β+ \cos α \sin β \nonumber\]

    to get

    \[\begin{align*} c_1 \cos (ωt)+c_2 \sin (ωt) &= A( \sin (ωt) \cos ϕ+ \cos (ωt) \sin ϕ) \\[4pt] &= A \sin ϕ( \cos (ωt))+A \cos ϕ( \sin (ωt)). \end{align*}\]

    Thus,

    \[c1=A \sin ϕ \text{ and } c_2=A \cos ϕ. \nonumber\]

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

    \[\begin{align*}c_1^2+c_2^2 &=A^2 \sin _2 ϕ+A^2 \cos _2 ϕ \\[4pt] &=A^2( \sin ^2 ϕ+ \cos ^2 ϕ) \\[4pt] &=A^2. \end{align*}\]

    Thus,

    \[A=\sqrt{c_1^2+c_2^2}. \nonumber\]

    Now, to find \(ϕ\), go back to the equations for \(c_1\) and \(c_2\), but this time, divide the first equation by the second equation to get

    \[\begin{align*} \dfrac{c_1}{c_2} &=\dfrac{A \sin ϕ}{A \cos ϕ} \\[4pt] &= \tan ϕ. \end{align*} \]

    Then,

    \[\tan ϕ= \dfrac{c_1}{c_2}. \nonumber\]

    We summarize this finding in the following theorem.

    Solution TO THE EQUATION FOR SIMPLE HARMONIC MOTION

    The function \(x(t)=c_1 \cos (ωt)+c_2 \sin (ωt)\) can be written in the form \(x(t)=A \sin (ωt+ϕ)\), where \(A=\sqrt{c_1^2+c_2^2}\) and \( \tan ϕ = \dfrac{c_1}{c_2}\).

    Note that when using the formula \( \tan ϕ=\dfrac{c_1}{c_2}\) to find \(ϕ\), we must take care to ensure \(ϕ\) is in the right quadrant (Figure \(\PageIndex{3}\)).

    This figure is the graph of f(t) = sin 2t. It is a periodic, oscillating graph. The period of the graph is represented with a line pointing from one peak to the next. It is labeled with the period T = 2π/ω. The graph has a phase shift of ϕ/ω so that the sine curve has the value zero to the left of the origin.
    Figure \(\PageIndex{3}\): A graph of vertical displacement versus time for simple harmonic motion with a phase change.

    Example \(\PageIndex{2}\): Expressing the Solution with a Phase Shift

    Express the following functions in the form \(A \sin (ωt+ϕ) \). What is the frequency of motion? The amplitude?

    1. \(x(t)=2 \cos (3t)+ \sin (3t)\)
    2. \(x(t)=3 \cos (2t)−2 \sin (2t)\)

    Solution

    We have

    \[A=\sqrt{c_1^2+c_2^2}=\sqrt{2^2+1^2}=\sqrt{5} \nonumber\]

    and

    \[ \tan ϕ = \dfrac{c_1}{c_2}=\dfrac{2}{1}=2. \nonumber\]

    Note that both \(c_1\) and \(c_2\) are positive, so \(ϕ\) is in the first quadrant. Thus,

    \[ϕ≈1.107 \; \text{rad}, \nonumber\]

    so we have

    \[ x(t) = 2 \cos (3t)+ \sin (3t) =5 \sin (3t+1.107). \nonumber\]

    The frequency is \(\dfrac{ω}{2π}=\dfrac{3}{2π}≈0.477.\) The amplitude is \(\sqrt{5}\).

    1. We have

      \[A=\sqrt{c_1^2+c_2^2}=\sqrt{3^2+2^2}=\sqrt{13} \nonumber\]

      and

      \[ \tan ϕ = \dfrac{c_1}{c_2}= \dfrac{3}{−2}=−\dfrac{3}{2}. \nonumber\]

      Note that \(c_1\) is positive but \(c_2\) is negative, so \(ϕ\) is in the fourth quadrant. Thus,

      \[ϕ≈−0.983 \; \text{rad}, \nonumber\]

      so we have

      \[\begin{align*} x(t) &=3 \cos (2t) −2 \sin (2t) \\ &= \sqrt{13} \sin (2t−0.983). \end{align*}\]

      The frequency is \(\dfrac{ω}{2π}=\dfrac{2}{2π}≈0.318.\) The amplitude is \(\sqrt{13}\).

    Exercise \(\PageIndex{2}\)

    Express the function \(x(t)= \cos (4t) + 4 \sin (4t)\) in the form \(A \sin (ωt+ϕ) \). What is the frequency of motion? The amplitude?

    Hint

    Use the process from the Example \(\PageIndex{2}\).

    Answer

    \(x(t)= \sqrt{17} \sin (4t+0.245), \text{frequency} =\dfrac{4}{2π}≈0.637, A=\sqrt{17}\)

    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; Figure \(\PageIndex{4}\)).

    This figure is a pneumatic cylinder. The cylinder is clear, and the piston can be seen.
    Figure \(\PageIndex{4}\): A dashpot is a pneumatic cylinder that dampens the motion of an oscillating system.

    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 \(−bx′\) for some constant \(b>0\). Again applying Newton's second law, the differential equation becomes

    \[mx″+bx′+kx=0. \nonumber\]

    Then the associated characteristic equation is

    \[mλ^2+bλ+k=0. \nonumber\]

    Applying the quadratic formula, we have

    \[λ=\dfrac{−b±\sqrt{b^2−4mk}}{2m}. \nonumber\]

    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: Overdamped Vibrations

    When \(b^2>4mk\), we say the system is overdamped . The general solution has the form

    \[x(t)=c_1e^{λ_1t}+c_2e^{λ_2t}, \nonumber\]

    where both \(λ_1\) and \(λ_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 \(\PageIndex{5}\) shows what typical critically damped behavior looks like.

    This figure has two graphs labeled (a) and (b). The first graph is a decreasing curve with the horizontal axis as a horizontal asymptote. The second graph initially is a decreasing function but becomes increasing below the horizontal axis. Then, the horizontal axis is also a horizontal asymptote.
    Figure \(\PageIndex{5}\): Behavior of an overdamped spring-mass system, with no change in direction (a) and only one change in direction (b).

    Example \(\PageIndex{3}\): 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 5252 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 \(\dfrac{16}{3}\) ft. Thus, \(16=\left(\dfrac{16}{3}\right)k,\) so \(k=3.\) We also have \(m=\dfrac{16}{32}=\dfrac{1}{2}\), so the differential equation is

    \[\dfrac{5}{2}x′+3x=0. \nonumber\]

    Multiplying through by 2 gives \(x″+5x′+6x=0\), which has the general solution

    \[x(t)=c_1e^{−2t}+c_2e^{−3t}. \nonumber\]

    Applying the initial conditions, \(x(0)=0\) and \(x′(0)=−5\), we get

    \[x(t)=−5e^{−2t}+5e^{−3t}. \nonumber\]

    After 10 sec the mass is at position

    \[x(10)=−5e^{−20}+5e^{−30}≈−1.0305×10^{−8}≈0, \nonumber\]

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

    \[x′(10)=10e^{−20}−15e^{−30}≈2.061×10^{−8}≈0. \nonumber\]

    After only 10 sec, the mass is barely moving.

    Exercise \(\PageIndex{3}\)

    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.

    Hint

    Follow the process from the previous example.

    Answer

    \(x(t)=0.6e^{−2t}−0.2e^{−6t}\)

    Case 2: Critically Damped Vibrations

    When \(b^2=4mk\), we say the system is critically damped. The general solution has the form

    \[x(t)=c_1e^{λ_1t}+c_2te^{λ_1t}, \nonumber\]

    where \(λ_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 \(\PageIndex{6}\) shows what typical critically damped behavior looks like.

    This figure has two graphs labeled (a) and (b). The first graph is in the first quadrant and is a decreasing curve with the horizontal axis as a horizontal asymptote. The second graph initially is a decreasing function but becomes increasing below the horizontal axis. Then, the horizontal axis is also a horizontal asymptote.
    Figure \(\PageIndex{6}\): Behavior of a critically damped spring-mass system. The system graphed in part (a) has more damping than the system graphed in part (b).

    Example \(\PageIndex{4}\): 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(9.8)=0.2k\), so \(k=49.\) Then, the differential equation is

    \[x″+14x′+49x=0, \nonumber\]

    which has general solution

    \[x(t)=c_1e^{−7t}+c_2te^{−7t}. \nonumber\]

    Applying the initial conditions \(x(0)=0\) and \(x′(0)=−3\) gives

    \[x(t)=−3te^{−7t}. \nonumber\]

    Exercise \(\PageIndex{4}\)

    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.

    Hint

    First find the spring constant.

    Answer

    \(x(t)=\dfrac{1}{2}e^{−8t}+4te^{−8t} \)

    Case 3: Undamped Vibrations

    When \(b^2<4mk\), we say the system is underdamped. The general solution has the form

    \[x(t)=e^{αt}(c_1 \cos (βt) + c_2 \sin (βt)), \nonumber\]

    where \(α\) 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 \(\PageIndex{7}\) shows what typical underdamped behavior looks like.

    This figure is an oscillating graph where the amplitude is decreasing. There are red dashed curves at the peaks of the amplitudes showing the pattern of a decreasing amplitude. As t increases, the horizontal axis becomes a horizontal asymptote.
    Figure \(\PageIndex{7}\): Behavior of an underdamped spring-mass system.

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

    Example \(\PageIndex{5}\): 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=\dfrac{16}{3.2}=5\) and \(m=\dfrac{16}{32}=\dfrac{1}{2},\) so the differential equation is

    \[\dfrac{1}{2} x″+x′+5x=0, \; \text{or} \; x″+2x′+10x=0. \nonumber\]

    This equation has the general solution

    \[x(t)=e^{−t} ( c_1 \cos (3t)+c_2 \sin (3t) ) . \nonumber\]

    Applying the initial conditions, \(x(0)=\dfrac{3}{4}\) and \(x′(0)=0,\) we get

    \[x(t)=e^{−t} \bigg( \dfrac{3}{4} \cos (3t)+ \dfrac{1}{4} \sin (3t) \bigg) . \nonumber\]

    Exercise \(\PageIndex{5}\)

    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.

    Hint

    First find the spring constant.

    Answer

    \(x(t)=−0.24e^{−2t} \cos (4t)−0.12e^{−2t} \sin (4t) \)

    Example \(\PageIndex{6}\): 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 picture is a shock absorber on a motorcycle.
    Figure \(\PageIndex{8}\): (credit: modification of work by nSeika, Flickr)

    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 \(\PageIndex{9}\)).

    This figure has two springs attached above at a fixed point. The first spring is labeled,
    Figure \(\PageIndex{9}\): We can use a spring-mass system to model a motorcycle suspension.

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

    \[mx″+bx′+kx=0. \nonumber\]

    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).

    1. Set up the differential equation that models the behavior of the motorcycle suspension system.
    2. 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.
    3. Graph the equation of motion over the first second after the motorcycle hits the ground.

    Solution

    1. We have defined equilibrium to be the point where \(mg=ks\), so we have

      \[\begin{align*} mg &=ks \\ 384 &=k\left(\dfrac{1}{3}\right)\\ k &=1152. \end{align*}\]

      We also have

      \[\begin{align*} W &=mg \\ 384 &=m(32) \\ m &=12. \end{align*}\]

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

      \[12x″+240x′+1152x=0. \nonumber\]

      Dividing through by 12, we get

      \[x''+20x′+96x=0. \nonumber\]

    2. The differential equation found in part a. has the general solution

      \[x(t)=c_1e^{−8t}+c_2e^{−12t}. \nonumber \]

      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(\dfrac{1}{3}\text{ ft}\right)\) below the equilibrium position (with respect to the motorcycle frame), and we have \(x(0)=\dfrac{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′(0)=10.\) Applying these initial conditions, we get \(c_1=\dfrac{7}{2}\) and \(c_2=−\left(\dfrac{19}{6}\right)\),so the equation of motion is

      \[x(t)=\dfrac{7}{2}e^{−8t}−\dfrac{19}{6}e^{−12t}. \nonumber\]

    3. The graph is shown in Figure \(\PageIndex{10}\).
    This figure is the graph of the function f(x) = 7/2e^−8t −19/6e^−12t.  The vertical axis is scaled in increments of 0.05, and the horizontal axis is labeled,
    Figure \(\PageIndex{10}\): Graph of the equation of motion over a time of one second.

    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/sec2, whereas on Mars it is 3.7 m/sec2.

    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 \(\PageIndex{11}\).

    This figure has three images. The first is a picture of the Mars Lander landing on a surface. The second picture is a diagram of the Mars Lander at touchdown, with an uncompressed spring with length L between the Lander and the landing surface. The third image is a diagram of the Lander in equilibrium position after the Lander has landed. The spring is compressed a distance of s.
    Figure \(\PageIndex{11}\): The landing craft suspension can be represented as a damped spring-mass system. (credit "lander": NASA)

    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:

    \[mx''+bx′+kx=0, \nonumber\]

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

    1. 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.
    2. 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.
    3. 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?
    4. 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?
    5. 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?
    6. 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

    \[mx″+bx′+kx=f(t), \nonumber\]

    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

    \[x(t)=c_1x_1(t)+c_2x_2(t)+x_p(t), \nonumber\]

    where \(c_1x_1(t)+c_2x_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, \(\lim \limits_{t \to \infty} c_1x_1(t)+c_2x_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.

    Example \(\PageIndex{7}\): 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 \sin (4t)\) 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(32)=2k,\) so \(k=16\) and the differential equation is

    \[x″+8x′+16x=8 \sin (4t). \nonumber\]

    The general solution to the complementary equation is

    \[c_1e^{−4t}+c_2te^{−4t}. \nonumber\]

    Assuming a particular solution of the form \(x_p(t)=A \cos (4t)+ B \sin (4t)\) and using the method of undetermined coefficients, we find \(x_p (t)=−\dfrac{1}{4} \cos (4t)\), so

    \[x(t)=c_1e^{−4t}+c_2te^{−4t}−\dfrac{1}{4} \cos (4t). \nonumber\]

    At \(t=0,\) the mass is at rest in the equilibrium position, so \(x(0)=x′(0)=0.\) Applying these initial conditions to solve for \(c_1\) and \(c_2,\) we get

    \[x(t)=\dfrac{1}{4}e^{−4t}+te^{−4t}−\dfrac{1}{4} \cos (4t). \nonumber\]

    The transient solution is \(\dfrac{1}{4}e^{−4t}+te^{−4t}\). The steady-state solution is \(−\dfrac{1}{4} \cos (4t).\)

    Exercise \(\PageIndex{6}\)

    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)=68e^{−2}t \cos (4t) \) 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?

    Hint

    Find the particular solution before applying the initial conditions.

    Answer

    \(x(t)=−\dfrac{1}{2} \cos (4t)+ \dfrac{9}{4} \sin (4t)+ \dfrac{1}{2} e^{−2t} \cos (4t)−2e^{−2t} \sin (4t)\)

    \(\text{Transient solution:} \dfrac{1}{2}e^{−2t} \cos (4t)−2e^{−2t} \sin (4t)\)

    \(\text{Steady-state solution:} −\dfrac{1}{2} \cos (4t)+ \dfrac{9}{4} \sin (4t) \)

    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=\dfrac{1}{T}=\dfrac{ω}{2π}\), 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.

    1. Consider the differential equation \(x″+x=0.\) Find the general solution. What is the natural frequency of the system?
    2. Now suppose this system is subjected to an external force given by \(f(t)=5 \cos t.\) Solve the initial-value problem \(x″+x=5 \cos t\), \(x(0)=0\), \(x′(0)=1\).
    3. Graph the solution. What happens to the behavior of the system over time?
    4. 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.

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

      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 of the Tacoma Narrows Bridge "Gallopin' Gertie". https://www.youtube.com/watch?v=j-zczJXSxnw

    1. 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 TV show Mythbusters aired an episode on this phenomenon. Visit this website to learn more about it. Adam Savage also described the experience. Watch this video for his account.

    The RLC Series Circuit

    Consider an electrical circuit containing a resistor, an inductor, and a capacitor, as shown in Figure \(\PageIndex{12}\). 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 \((Ω)\), 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 \(\PageIndex{12}\) 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 and those of you interested in the derivation of these formulas can review the links. 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,

    \[E_L=L\dfrac{dI}{dt}. \nonumber\]

    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,

    \[E_R=RI. \nonumber\]

    Last, the voltage drop across a capacitor is proportional to the charge, \(q,\) on the capacitor, with proportionality constant \(1/C\). Thus,

    \[E_C=\dfrac{1}{C}q. \nonumber\]

    Adding these terms together, we get

    \[L\dfrac{dI}{dt}+RI+\dfrac{1}{C}q=E(t). \nonumber\]

    Noting that \(I=(dq)/(dt)\), this becomes

    \[L\dfrac{d^2q}{dt^2}+R\dfrac{dq}{dt}+\dfrac{1}{C}q=E(t). \nonumber\]

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

    This figure is a diagram of a circuit. It has broken lines at the bottom labeled C. On the left side there is an open circle labeled E. The top has diagonal lines labeled R. The right side has little bumps labeled L.
    Figure \(\PageIndex{12}\): An RLC series circuit can be modeled by the same differential equation as a mass-spring system.

    Series Circuit

    Find the charge on the capacitor in an RLC series circuit where \(L=5/3\) H, \(R=10Ω\), \(C=1/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

    \[\begin{align*} L\dfrac{d^2q}{dt^2}+R\dfrac{dq}{dt}+\dfrac{1}{C}q &=E(t) \\[4pt] \dfrac{5}{3} \dfrac{d^2q}{dt^2}+10\dfrac{dq}{dt}+30q &=300 \\[4pt] \dfrac{d^2q}{dt^2}+6\dfrac{dq}{dt}+18q &=180. \end{align*}\]

    The general solution to the complementary equation is

    \[e^{−3t}(c_1 \cos (3t)+c_2 \sin (3t)). \nonumber\]

    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,

    \[q(t)=e^{−3t}(c_1 \cos (3t)+c_2 \sin (3t))+10. \nonumber\]

    Applying the initial conditions \(q(0)=0\) and \(i(0)=((dq)/(dt))(0)=9,\) we find \(c_1=−10\) and \(c_2=−7.\) So the charge on the capacitor is

    \[q(t)=−10e^{−3t} \cos (3t)−7e^{−3t} \sin (3t)+10. \nonumber\]

    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.

    Exercise \(\PageIndex{7}\)

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

    Hint

    Remember, \(E_L=L((dI)/(dt)).\)

    Answer

    \[q(t)=−25e^{−t} \cos (3t)−7e^{−t} \sin (3t)+25 \nonumber\]

    Key Concepts

    • Second-order constant-coefficient differential equations can be used to model spring-mass systems.
    • An examination of the forces on a spring-mass system results in a differential equation of the form \[mx″+bx′+kx=f(t), \nonumber\] where mm represents the mass, bb is the coefficient of the damping force, \(k\) is the spring constant, and \(f(t)\) represents any net external forces on the system.
    • If \(b=0\), there is no damping force acting on the system, and simple harmonic motion results.
    • If \(b≠0\),the behavior of the system depends on whether \(b^2−4mk>0, b^2−4mk=0,\) or \(b^2−4mk<0.\)
      • If \(b^2−4mk>0,\) the system is overdamped and does not exhibit oscillatory behavior.
      • If \(b^2−4mk=0,\) the system is critically damped. It does not exhibit oscillatory behavior, but any slight reduction in the damping would result in oscillatory behavior.
      • If \(b^2−4mk<0\), the system is underdamped. It exhibits oscillatory behavior, but the amplitude of the oscillations decreases over time.
    • If\(f(t)≠0\), the solution to the differential equation is the sum of a transient solution and a steady-state solution. The steady-state solution governs the long-term behavior of the system.
    • The charge on the capacitor in an RLC series circuit can also be modeled with a second-order constant-coefficient differential equation of the form \[L\dfrac{d^2q}{dt^2}+R\dfrac{dq}{dt}+\dfrac{1}{C}q=E(t), \nonumber\] where \(L\) is the inductance, \(R\) is the resistance, \(C\) is the capacitance, and \(E(t)\) is the voltage source.

    Key Equations

    • Equation of simple harmonic motion \[x″+ω^2x=0 \nonumber\]
    • Solution for simple harmonic motion \[x(t)=c_1 \cos (ωt)+c_2 \sin (ωt) \nonumber\]
    • Alternative form of solution for SHM \[x(t)=A \sin (ωt+ϕ) \nonumber\]
    • Forced harmonic motion \[mx″+bx′+kx=f(t)\nonumber\]
    • Charge in a RLC series circuit \[L\dfrac{d^2q}{dt^2}+R\dfrac{dq}{dt}+\dfrac{1}{C}q=E(t),\nonumber\]

    Glossary

    RLC series circuit
    a complete electrical path consisting of a resistor, an inductor, and a capacitor; a second-order, constant-coefficient differential equation can be used to model the charge on the capacitor in an RLC series circuit
    simple harmonic motion
    motion described by the equation \(x(t)=c_1 \cos (ωt)+c_2 \sin (ωt)\), as exhibited by an undamped spring-mass system in which the mass continues to oscillate indefinitely
    steady-state solution
    a solution to a nonhomogeneous differential equation related to the forcing function; in the long term, the solution approaches the steady-state solution

    Contributors and Attributions

    • Gilbert Strang (MIT) and Edwin "Jed" Herman (Harvey Mudd) with many contributing authors. This content by OpenStax is licensed with a CC-BY-SA-NC 4.0 license. Download for free at http://cnx.org.

    Find the Steady-state Solution for Each Differential Equation.

    Source: https://math.libretexts.org/Bookshelves/Calculus/Book:_Calculus_%28OpenStax%29/17:_Second-Order_Differential_Equations/17.3:_Applications_of_Second-Order_Differential_Equations