Q4.9-3E - Fundamentals Of Differential Equations And Boundary Value Problems

Question

The motion of a mass-spring system with damping is governed by 

$\begin{array}{l}\mathbf{y}\mathbf{"}\left(t\right)\mathbf{+}\mathbf{by}\mathbf{\text{'}}\left(t\right)\mathbf{+}\mathbf{16}\mathbf{y}\left(t\right)\mathbf{=}\mathbf{0}\mathbf{;}\\ \mathbf{y}\left(0\right)\mathbf{=}\mathbf{1}\mathbf{,}\mathbf{y}\mathbf{\text{'}}\left(0\right)\mathbf{=}\mathbf{0}\mathbf{.}\end{array}$

Find the equation of motion and sketch its graph for b = 0,6,8, and 10.

Step-by-Step Solution

Answer

Therefore, the equations of motions are solution for b = 0 is,$\mathbf{y}\mathbf{=}\mathbf{cos}\mathbf{4}\mathbf{t}$ solution for 

b = 6 is $\mathbf{y}\left(t\right)\mathbf{=}{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{+}\frac{3}{\sqrt{7}}{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}$, solution for b = 8 is $\mathbf{y}\mathbf{=}\left(1+4t\right){\mathbf{e}}^{-4t}$ and solution for b = 10 is $\mathbf{y}\mathbf{=}\mathbf{-}\frac{1}{3}{\mathbf{e}}^{-8t}\mathbf{+}\frac{4}{3}{\mathbf{e}}^{-2t}$.

A Sketch of the graph is shown below.

Step 1: General form

The Mass–Spring Oscillator

A damped mass-spring oscillator consists of a mass m attached to a spring fixed at one end, as shown in Figure 4.1. Devise a differential equation that governs the motion of this oscillator, taking into account the forces acting on it due to the spring elasticity, damping friction, and possible external influences.

 

Mass–spring oscillator equation: 

 

 $\begin{array}{c}{\mathbf{F}}_{\mathrm{ext}}\mathbf{=}\left[\mathrm{inertia}\right]\mathbf{y}\mathbf{"}\mathbf{+}\left[\mathrm{damping}\right]\mathbf{y}\mathbf{\text{'}}\mathbf{+}\left[\mathrm{stiffness}\right]\mathbf{y}\\ =\mathrm{my}"+\mathrm{by}\text{'}+\mathrm{ky}\end{array}$  …… (1)

The rule for the bounded equation: Just based on stiffness we can decide whether it is bounded or not if stiffness k > 0 then it is bounded and if k < 0 then it is unbounded.

 

Step 2: Evaluate the equation

Given that, 

$\begin{array}{c}\mathbf{y}\mathbf{"}\left(t\right)\mathbf{+}\mathbf{by}\mathbf{\text{'}}\left(t\right)\mathbf{+}\mathbf{16}\mathbf{y}\left(t\right)\mathbf{=}\mathbf{0}\mathbf{;}\\ \mathbf{y}\left(0\right)\mathbf{=}\mathbf{1}\mathbf{,}\mathbf{y}\mathbf{\text{'}}\left(0\right)\mathbf{=}\mathbf{0}\mathbf{.}\end{array}$ 

 

And the value of.$\mathbf{b}\mathbf{=}\mathbf{0}\mathbf{,}\mathbf{6}\mathbf{,}\mathbf{8}\mathbf{,}\mathbf{10}$

 

Let us take b = 6. Then,

 

   $\mathbf{y}\mathbf{"}\mathbf{+}6\mathbf{y}\mathbf{\text{'}}\mathbf{+}16\mathbf{y}\mathbf{=}\mathbf{0}$…… (2)

 

Since the auxiliary equation is.${\mathbf{r}}^2\mathbf{+}\mathbf{6}\mathbf{r}\mathbf{+}\mathbf{16}\mathbf{=}\mathbf{0}$ And roots are $\mathbf{r}\mathbf{=}\mathbf{-}\mathbf{3}\mathbf{+}\sqrt{7}\mathbf{i}$ and $\mathbf{r}\mathbf{=}\mathbf{-}\mathbf{3}\mathbf{-}\sqrt{7}\mathbf{i}$.

Step3: Find the initial conditions

By the above information find the general solution.

 

The general solution is  $\mathbf{y}\left(t\right)\mathbf{=}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{+}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}$…… (3)

 

Given initial conditions are.$\mathbf{y}\left(0\right)\mathbf{=}1\mathbf{,}\mathbf{y}\mathbf{\text{'}}\left(0\right)\mathbf{=}0$

 

Now, substitute the initial conditions to find the value of c.

$\begin{array}{c}\mathbf{y}\left(t\right)\mathbf{=}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{+}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}\\ \mathbf{y}\left(0\right)\mathbf{=}{\mathbf{c}}_1{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{cos}\left(0\times \sqrt{7}\right)\mathbf{+}{\mathbf{c}}_2{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{sin}\left(\sqrt{7}\times 0\right)\\ 1\mathbf{=}{\mathbf{c}}_1\left(1\right)\mathbf{+}\mathbf{0}\\ {\mathbf{c}}_1\mathbf{=}1\end{array}$

Find the derivative of equation (3). And implement the initial condition.

 

$\mathbf{y}\mathbf{\text{'}}\left(t\right)\mathbf{=}\mathbf{-}\mathbf{3}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{-}\sqrt{7}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}\mathbf{-}\mathbf{3}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}\mathbf{+}\sqrt{7}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}$

 

Then,

 

$\begin{array}{c}\mathbf{y}\mathbf{\text{'}}\left(t\right)\mathbf{=}\mathbf{-}\mathbf{3}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{-}\sqrt{7}{\mathbf{c}}_1{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}\mathbf{-}\mathbf{3}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}\mathbf{+}\sqrt{7}{\mathbf{c}}_2{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\\ \mathbf{y}\mathbf{\text{'}}\left(0\right)\mathbf{=}\mathbf{-}\mathbf{3}{\mathbf{c}}_1{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{cos}\left(0\right)\mathbf{-}\sqrt{7}{\mathbf{c}}_1{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{sin}\left(0\right)\mathbf{-}\mathbf{3}{\mathbf{c}}_2{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{sin}\left(0\right)\mathbf{+}\sqrt{7}{\mathbf{c}}_2{\mathbf{e}}^{\mathbf{-}\mathbf{3}\left(0\right)}\mathbf{cos}\left(0\right)\\ \mathbf{0}\mathbf{=}\mathbf{-}\mathbf{3}{\mathbf{c}}_1\mathbf{+}\sqrt{7}{\mathbf{c}}_2\left(1\right)\\ {\mathbf{c}}_2\mathbf{=}\frac{3}{\sqrt{7}}\\ \mathbf{=}\frac{1}{2}\end{array}$

Now substitute the value of c in equation (3).

 $\mathbf{y}\left(t\right)\mathbf{=}{\mathbf{e}}^{-3t}\mathbf{cos}\sqrt{7}\mathbf{t}\mathbf{+}\frac{3}{\sqrt{7}}{\mathbf{e}}^{-3t}\mathbf{sin}\sqrt{7}\mathbf{t}$

.

Similarly, solution for b = 0 is, $\mathbf{y}\mathbf{=}\mathbf{cos}\mathbf{4}\mathbf{t}$ solution for b = 8 is $\mathbf{y}\mathbf{=}\left(1+4t\right){\mathbf{e}}^{-4t}$and solution for b = 10 is $\mathbf{y}\mathbf{=}\mathbf{-}\frac{1}{3}{\mathbf{e}}^{-8t}\mathbf{+}\frac{4}{3}{\mathbf{e}}^{-2t}$.

Step4: Graph of the motion equation

The graph of the equation is.

 

When b = 0,6,8, and 10 means draw the graph.$\mathbf{y}\mathbf{"}\left(t\right)\mathbf{+}\mathbf{by}\mathbf{\text{'}}\left(t\right)\mathbf{+}\mathbf{16}\mathbf{y}\left(t\right)\mathbf{=}\mathbf{0}$

 

Case (1):

Case (2):

 

Case (3):

Case (4):