4.4 HW 4

  4.4.1 Problem 1
  4.4.2 Problem 2
  4.4.3 Problem 3
  4.4.4 Problem 4
  4.4.5 Problem 5
  4.4.6 HW 4 key solution

This HW in one PDF

4.4.1 Problem 1

   4.4.1.1 part(1)
   4.4.1.2 part(2)
   4.4.1.3 Part(3)
   4.4.1.4 Part(4)

pict

SOLUTION:

Note that \(\lambda _{critical}=\omega _{0}\). We are told that \(\lambda =0.1\omega _{0}\) in this problem.

4.4.1.1 part(1)

The resonant frequency (for this case of under-damped) occurs when the steady state amplitude is maximum\[ b=\frac{\frac{f}{m}}{\sqrt{\left ( \omega _{0}^{2}-\omega ^{2}\right ) ^{2}+4\lambda ^{2}\omega ^{2}}}\] This happens when the denominator is minimum. Taking derivative of the denominator w.r.t. \(\omega \) and setting the result to zero gives\begin{align*} \frac{d}{d\omega }\left ( \left ( \omega _{0}^{2}-\omega ^{2}\right ) ^{2}+4\lambda ^{2}\omega ^{2}\right ) & =0\\ 2\left ( \omega _{0}^{2}-\omega ^{2}\right ) \left ( -2\omega \right ) +8\lambda ^{2}\omega & =0\\ 8\lambda ^{2}\omega +4\omega ^{3}-4\omega \omega _{0}^{2} & =0\\ 2\lambda ^{2}+\omega ^{2}-\omega _{0}^{2} & =0\\ \omega ^{2} & =\omega _{0}^{2}-2\lambda ^{2} \end{align*}

Taking the positive root (since \(\omega \) must be positive) gives\[ \omega =\sqrt{\omega _{0}^{2}-2\lambda ^{2}}\] When \(\lambda =0.1\omega _{0}\) the above becomes\begin{align*} \omega & =\sqrt{\omega _{0}^{2}-2\left ( \frac{1}{10}\omega _{0}\right ) ^{2}}\\ & =\sqrt{\frac{98}{100}\omega _{0}^{2}}\\ & =0.98995\omega _{0}\text{ rad/sec} \end{align*}

4.4.1.2 part(2)

Quality factor \(Q\) is defined as\begin{align*} Q & =\frac{\omega _{d}}{2\lambda }\\ & =\frac{\sqrt{\omega _{0}^{2}-\lambda ^{2}}}{2\lambda }\\ & =\frac{\sqrt{\omega _{0}^{2}-\left ( 0.1\omega _{0}\right ) ^{2}}}{2\left ( 0.1\omega _{0}\right ) }\\ & =\frac{\omega _{0}\sqrt{1-0.1^{2}}}{0.2\omega _{0}}\\ & =\frac{\sqrt{1-0.1^{2}}}{0.2} \end{align*}

Therefore\[ \fbox{$Q=4.975$}\]

4.4.1.3 Part(3)

Given \begin{equation} x^{\prime \prime }\left ( t\right ) +2\lambda x^{\prime }+\omega _{0}^{2}x=\frac{f}{m}e^{i\omega t} \tag{1} \end{equation} Assuming the particular solution is \(x_{p}\left ( t\right ) =Be^{i\omega t}\) where \(B=be^{i\phi }\) is the complex amplitude and \(b\) is the amplitude and \(\phi \) is the phase of \(B\). We want to find the phase. Plugging \(x_{p}\left ( t\right ) \) into (1) and simplifying gives\[ B=\frac{\frac{f}{m}}{\omega _{0}^{2}-\omega ^{2}+2\lambda i\omega }\] Hence\begin{align*} \phi & =0-\tan ^{-1}\left ( \frac{2\lambda \omega }{\omega _{0}^{2}-\omega ^{2}}\right ) \\ & =\tan ^{-1}\left ( \frac{-2\lambda \omega }{\omega _{0}^{2}-\omega ^{2}}\right ) \end{align*}

Since \(\lambda =0.1\omega _{0}\) and \(\omega =\frac{\omega _{0}}{2}\) the above becomes\begin{align*} \phi & =\tan ^{-1}\left ( \frac{-2\left ( 0.1\omega _{0}\right ) \frac{\omega _{0}}{2}}{\omega _{0}^{2}-\left ( \frac{\omega _{0}}{2}\right ) ^{2}}\right ) \\ & =\tan ^{-1}\left ( -0.13333\right ) \\ & =-0.13255\text{ rad} \end{align*}

4.4.1.4 Part(4)

The steady state amplitude is \(b\) from above, which is found as follows\[ b^{2}=BB^{\ast }\] Where \(B^{\ast }\) is the complex conjugate of \(B=\frac{\frac{f}{m}}{\omega _{0}^{2}-\omega ^{2}+2\lambda i\omega }\). Therefore\begin{align*} b & =\frac{\frac{f}{m}}{\sqrt{\left ( \omega _{0}^{2}-\omega ^{2}\right ) ^{2}+4\lambda ^{2}\omega ^{2}}}\\ & =\frac{f}{m}\frac{1}{\sqrt{\left ( \omega _{0}^{2}-\left ( \frac{\omega _{0}}{2}\right ) ^{2}\right ) ^{2}+4\left ( 0.1\omega _{0}\right ) ^{2}\left ( \frac{\omega _{0}}{2}\right ) ^{2}}}\\ & =\frac{f}{m}\frac{1}{\sqrt{0.572\,5\omega _{0}^{4}}}\\ & =1.3216\frac{f}{m\omega _{0}^{2}} \end{align*}

But \(m\omega _{0}^{2}=k\), the stiffness, hence the above is\[ \fbox{$b=1.3216\frac{f}{k}$}\]

4.4.2 Problem 2

pict

SOLUTION:

The locus the bead describes is an ellipse, since in an ellipse the total distance from any point on it to the points \(A,B\) is always the same

pict

To obtain the potential energy, we move the bead a little from the origin and find how much the bead moved above the origin, as shown in the following diagram

pict

From the above, we see that, by applying pythagoras triangle theorem to the left and to the right triangles, we obtain two equations which we solve for \(h\) in order to obtain the potential energy\begin{align*} s^{2} & =h^{2}+\left ( d+x\right ) ^{2}\\ \left ( 2l-s\right ) ^{2} & =h^{2}+\left ( d-x\right ) ^{2} \end{align*}

Solving for \(h\) gives\[ h=\sqrt{1-\frac{d^{2}}{l^{2}}}\sqrt{l^{2}-x^{2}}\] Therefore \begin{align*} y & =H-h\\ & =H-\sqrt{1-\frac{d^{2}}{l^{2}}}\sqrt{l^{2}-x^{2}} \end{align*}

Hence \begin{align*} U & =mgy\\ & =mg\left ( H-\sqrt{1-\frac{d^{2}}{l^{2}}}\sqrt{l^{2}-x^{2}}\right ) \end{align*}

The kinetic energy is\[ T=\frac{1}{2}m\left ( \dot{x}^{2}+\dot{y}^{2}\right ) \] Therefore the Lagrangian is\begin{align*} L & =T-U\\ & =\frac{1}{2}m\left ( \dot{x}^{2}+\dot{y}^{2}\right ) -mg\left ( H-\sqrt{1-\frac{d^{2}}{l^{2}}}\sqrt{l^{2}-x^{2}}\right ) \end{align*}

The equation of motion in the \(x\) coordinate is now found. From\begin{align*} \frac{\partial L}{\partial x} & =\frac{1}{2}mg\sqrt{1-\frac{d^{2}}{l^{2}}}\frac{\left ( -2x\right ) }{\sqrt{l^{2}-x^{2}}}\\ & =-mg\sqrt{1-\frac{d^{2}}{l^{2}}}\frac{x}{\sqrt{l^{2}-x^{2}}} \end{align*}

And \[ \frac{d}{dt}\frac{\partial L}{\partial x}=m\ddot{x}\] Applying Euler-Lagrangian equation gives\begin{align*} \frac{d}{dt}\frac{\partial L}{\partial x}-\frac{\partial L}{\partial x} & =0\\ \ddot{x}+g\sqrt{1-\frac{d^{2}}{l^{2}}}\frac{x}{\sqrt{l^{2}-x^{2}}} & =0 \end{align*}

For very small \(x\), we drop the \(x^{2}\) term and the above reduces to\[ \ddot{x}+g\sqrt{1-\frac{d^{2}}{l^{2}}}\frac{x}{l}=0 \] Hence the undamped natural frequency is\[ \omega _{0}^{2}=\frac{g}{l}\sqrt{1-\frac{d^{2}}{l^{2}}}\] or\[ \omega _{0}=\sqrt{\frac{g}{l}\sqrt{1-\frac{d^{2}}{l^{2}}}}\] The period of small oscillation is therefore\begin{align*} T & =\frac{2\pi }{\omega _{0}}\\ & =2\pi \frac{1}{\sqrt{\frac{g}{l}\sqrt{1-\frac{d^{2}}{l^{2}}}}} \end{align*}

4.4.3 Problem 3

   4.4.3.1 SOLUTION method one
   4.4.3.2 another solution

pict

4.4.3.1 SOLUTION method one

The velocity of the particle is as shown in the following diagram

pict

There is no potential energy, and the Lagrangian only comes from kinetic energy.

\begin{align*} v^{2} & =V_{x}^{2}+V_{y}^{2}\\ & =\left ( \dot{r}\cos \theta -r\omega \sin \theta \right ) ^{2}+\left ( \dot{r}\sin \theta +r\omega \cos \theta \right ) ^{2} \end{align*}

Exapnding and simplifying gives\[ v^{2}=\dot{r}^{2}+r^{2}\omega ^{2}\] Hence\[ L=\frac{1}{2}m\left ( \dot{r}^{2}+r^{2}\omega ^{2}\right ) \] And the equation of motion in the radial \(r\) direction is\begin{align*} \frac{d}{dt}\frac{\partial L}{\partial \dot{r}}-\frac{\partial L}{\partial r} & =0\\ \frac{d}{dt}m\dot{r}-mr\omega ^{2} & =0 \end{align*}

Hence the equation of motion is\begin{equation} \fbox{$\ddot{r}-r\omega ^2=0$}\tag{1} \end{equation} The roots of the characteristic equation are \(\pm \omega \), hence the solution is\[ r\left ( t\right ) =c_{1}e^{\omega t}+c_{2}e^{-\omega t}\] At \(t=0\), \(r\left ( 0\right ) =0\) and \(\dot{r}\left ( t\right ) =L\omega \). Using these we can find \(c_{1},c_{2}\). \begin{equation} 0=c_{1}+c_{2}\tag{2} \end{equation} But \(\dot{r}\left ( t\right ) =\omega c_{1}e^{\omega t}-\omega c_{2}e^{-\omega t}\) and at \(t=0\) this becomes\begin{equation} L\omega =\omega c_{1}-\omega c_{2}\tag{3} \end{equation} From (2,3) we solve for \(c_{1},c_{2}\). From (2), \(c_{1}=-c_{2}\) and (3) becomes\begin{align*} L\omega & =-\omega c_{2}-\omega c_{2}\\ c_{2} & =\frac{L\omega }{-2\omega }=\frac{-1}{2}L \end{align*}

Hence \(c_{1}=\frac{1}{2}L\) and the solution is\begin{align*} r\left ( t\right ) & =c_{1}e^{\omega t}+c_{2}e^{-\omega t}\\ & =\frac{1}{2}Le^{\omega t}-\frac{1}{2}Le^{-\omega t}\\ & =L\left ( \frac{e^{\omega t}-e^{-\omega t}}{2}\right ) \end{align*}

Or\[ \fbox{$r\left ( t\right ) =L\left ( \sinh \omega t\right ) $}\] To find the time it takes to reach end of rod, we solve for \(t_{p}\) from\begin{align*} L & =L\left ( \sinh \omega t_{p}\right ) \\ 1 & =\sinh \omega t_{p} \end{align*}

Hence\begin{align*} \omega t_{p} & =\sinh ^{-1}\left ( 1\right ) \\ & =0.88137 \end{align*}

Therefore\[ \fbox{$t_p=\frac{0.88137}{\omega }$ sec}\]

4.4.3.2 another solution

Let the local coordinate frame rotate with the bar, where the bar is oriented along the \(x\) axis of the local body coordinate frame as shown below.

pict

The position vector of the particle is \(\mathbf{r}=\mathbf{i}r\) where \(\mathbf{i}\) is unit vector along the \(x\) axis. Taking time derivative, and using the rotating vector time derivative rule which says that \(\frac{d\mathbf{A}}{dt}=\left ( \frac{d\mathbf{A}}{dt}\right ) _{relative}+\boldsymbol{\omega }\times \mathbf{A}\) where \(\boldsymbol{\omega }\) is the angular velocity of the rotating frame then\begin{equation} \mathbf{\dot{r}}=\mathbf{\dot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{r}\tag{1} \end{equation} To find the acceleration of the particle, we take time derivative one more time\[ \frac{d}{dt}\mathbf{\dot{r}}=\frac{d}{dt}\left ( \mathbf{\dot{r}}_{rel}\right ) +\boldsymbol{\dot{\omega }}\times \mathbf{r}+\boldsymbol{\omega }\times \mathbf{\dot{r}}\] But \(\frac{d}{dt}\left ( \mathbf{\dot{r}}_{rel}\right ) =\mathbf{\ddot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}\) by applying the rule of time derivative of rotating vector again. Therefore the above equation becomes\[ \frac{d}{dt}\mathbf{\dot{r}}=\mathbf{\ddot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}+\boldsymbol{\dot{\omega }}\times \mathbf{r}+\boldsymbol{\omega }\times \mathbf{\dot{r}}\] Replacing \(\mathbf{\dot{r}}\) in the above from its value in (1) gives \begin{align*} \mathbf{\ddot{r}} & =\mathbf{\ddot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}+\boldsymbol{\dot{\omega }}\times \mathbf{r}+\boldsymbol{\omega }\times \left ( \mathbf{\dot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{r}\right ) \\ & =\mathbf{\ddot{r}}_{rel}+\boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}+\boldsymbol{\dot{\omega }}\times \mathbf{r}+\boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}+\mathbf{\omega \times }\left ( \mathbf{\omega }\times \mathbf{r}\right ) \\ & =\mathbf{\ddot{r}}_{rel}+2\left ( \boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}\right ) +\boldsymbol{\dot{\omega }}\times \mathbf{r}+\boldsymbol{\omega }\times \left ( \boldsymbol{\omega }\times \mathbf{r}\right ) \end{align*}

But \(\boldsymbol{\omega }\) is constant (bar rotate with constant angular speed), hence the term \(\boldsymbol{\dot{\omega }}\) above is zero, and the above reduces to\begin{equation} \mathbf{\ddot{r}}=\mathbf{\ddot{r}}_{rel}+2\left ( \boldsymbol{\omega }\times \mathbf{\dot{r}}_{rel}\right ) +\boldsymbol{\omega }\times \left ( \boldsymbol{\omega }\times \mathbf{r}\right ) \tag{2} \end{equation} The above is the acceleration of the particle as seen in the inertial frame. Now we calculate this acceleration by preforming the vector operations above, noting that \(\mathbf{r}=\mathbf{i}r,\boldsymbol{\omega }=\mathbf{k}\omega \), hence (2) becomes\begin{align*} \mathbf{\ddot{r}} & =\mathbf{i}\ddot{r}_{rel}+2\left ( \mathbf{k}\omega \times \mathbf{i}\dot{r}_{rel}\right ) +\mathbf{k}\omega \mathbf{\times }\left ( \mathbf{k}\omega \times \mathbf{i}r\right ) \\ & =\mathbf{i}\ddot{r}_{rel}+2\left ( \mathbf{j}\omega \dot{r}_{rel}\right ) +\mathbf{k}\omega \mathbf{\times }\left ( \mathbf{j}\omega r\right ) \\ & =\mathbf{i}\ddot{r}_{rel}+2\left ( \mathbf{j}\omega \dot{r}_{rel}\right ) -\mathbf{i}\omega ^{2}r\\ & =\mathbf{i}\left ( \ddot{r}_{rel}-\omega ^{2}r\right ) +\mathbf{j}\left ( 2\omega \dot{r}_{rel}\right ) \end{align*}

The particle has an acceleration along \(x\) axis and an acceleration along \(y\) axis. We are interested in the acceleration along \(x\) since this is where the rod is oriented along. The scalar version of the acceleration in the \(x\) direction is\[ a_{x}=\ddot{r}_{rel}-\omega ^{2}r \] Using \(F_{x}=ma_{x}\) and since \(F_{x}=0\) (there is no force on the particle) then the equation of motion along the bar (\(x\) axis) is\[ \ddot{r}_{rel}-\omega ^{2}r=0 \] The roots of the characteristic equation is \(\pm \omega \), hence the solution is\[ r\left ( t\right ) =c_{1}e^{\omega t}+c_{2}e^{-\omega t}\] At \(t=0\), \(r\left ( 0\right ) =0\) and \(\dot{r}\left ( t\right ) =L\omega \). Using these we can find \(c_{1},c_{2}\). \begin{equation} 0=c_{1}+c_{2}\tag{3} \end{equation} But \(\dot{r}\left ( t\right ) =\omega c_{1}e^{\omega t}-\omega c_{2}e^{-\omega t}\) and at \(t=0\) this becomes\begin{equation} L\omega =\omega c_{1}-\omega c_{2}\tag{4} \end{equation} From (3,4) we solve for \(c_{1},c_{2}\). From (3), \(c_{1}=-c_{2}\) and (4) becomes\begin{align*} L\omega & =-\omega c_{2}-\omega c_{2}\\ c_{2} & =\frac{L\omega }{-2\omega }=\frac{-1}{2}L \end{align*}

Hence \(c_{1}=\frac{1}{2}L\) and the solution is\begin{align*} r\left ( t\right ) & =c_{1}e^{\omega t}+c_{2}e^{-\omega t}\\ & =\frac{1}{2}Le^{\omega t}-\frac{1}{2}Le^{-\omega t}\\ & =L\left ( \frac{e^{\omega t}-e^{-\omega t}}{2}\right ) \\ & =L\left ( \sinh \omega t\right ) \end{align*}

To find the time it takes to reach end of rod, we solve for \(t_{p}\) from\begin{align*} L & =L\left ( \sinh \omega t_{p}\right ) \\ 1 & =\sinh \omega t_{p} \end{align*}

Hence\begin{align*} \omega t_{p} & =\sinh ^{-1}\left ( 1\right ) \\ & =0.88137 \end{align*}

Therefore\[ \fbox{$t_p=\frac{0.88137}{\omega }$ sec}\]

4.4.4 Problem 4

   4.4.4.1 case (1)
   4.4.4.2 case (2)
   4.4.4.3 case(3)

pict

SOLUTION:

Starting with the equation of motion for damped oscillator\[ x^{\prime \prime }+2\lambda x^{\prime }+\omega _{0}^{2}x=0 \] The solution for cases 1,2 (both are underdamped) is\begin{equation} x=e^{-\lambda t}\left ( A\cos \omega _{d}t+B\sin \omega _{d}t\right ) \tag{1} \end{equation} Where \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}\). While the solution for case (3), the critical damped case is\begin{equation} x=\left ( A+tB\right ) e^{-\lambda t} \tag{2} \end{equation} For (1) above, at \(t=0\) we obtain\[ 1=A \] Hence (1) becomes \(x=e^{-\lambda t}\left ( \cos \omega _{d}t+B\sin \omega _{d}t\right ) \), and taking derivative gives\[ \dot{x}=-\lambda e^{-\lambda t}\left ( \cos \omega _{d}t+B\sin \omega _{d}t\right ) +e^{-\lambda t}\left ( -\omega _{d}\sin \omega _{d}t+B\omega _{d}\cos \omega _{d}t\right ) \] At \(t=0\) we have\begin{align*} 0 & =-\lambda +B\omega _{d}\\ B & =\frac{\lambda }{\omega _{d}} \end{align*}

Hence the complete solution for (1) is\begin{align} x & =e^{-\lambda t}\left ( \cos \omega _{d}t+\frac{\lambda }{\omega _{d}}\sin \omega _{d}t\right ) \tag{3}\\ \dot{x} & =-\lambda x+e^{-\lambda t}\left ( -\omega _{d}\sin \omega _{d}t+\lambda \cos \omega _{d}t\right ) \tag{4} \end{align}

Now we find the solution for (2), the critical damped case. At \(t=0\)\[ 1=A \] Hence (2) becomes \(x=\left ( 1+tB\right ) e^{-\lambda t}\), and taking derivative gives\[ \dot{x}=Be^{-\lambda t}-\lambda \left ( 1+tB\right ) e^{-\lambda t}\] At \(t=0\)\begin{align*} 0 & =B-\lambda \\ B & =\lambda \end{align*}

Hence the solution to (2) becomes\begin{align} x & =\left ( 1+\lambda t\right ) e^{-\lambda t}\tag{5}\\ \dot{x} & =\lambda e^{-\lambda t}-\lambda \left ( 1+\lambda t\right ) e^{-\lambda t} \tag{6} \end{align}

Now that the solutions are found, we plot the phase space using the computer, using parametric plot command

4.4.4.1 case (1)

For \(\lambda =0.05\), and \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}=\sqrt{0.5^{2}-0.05^{2}}=\allowbreak 0.4975\), then equations (3,4) become\begin{align} x & =e^{-0.05t}\left ( \cos 0.4975t+0.1005\sin 0.4975t\right ) \tag{3A}\\ \dot{x} & =-0.05x+e^{-0.05t}\left ( -0.4975\sin 0.4975t+0.05\cos 0.4975t\right ) \tag{4A} \end{align}

Here is the plot generated, showing starting point \(\left ( 1,0\right ) \) with the code used

pict

am = 0.05;  
wn = 0.5;  
wd = Sqrt[wn^2 - lam^2];  
x = Exp[-lam t] (Cos[wd t] + lam/wd Sin[wd t]);  
y = -lam x + Exp[-lam t] (-wd Sin[wd t] + lam Cos[lam t]);  
ParametricPlot[{x, y}, {t, 0, 50}, Frame -> True,  
 GridLines -> Automatic, GridLinesStyle -> LightGray,  
 FrameLabel -> {{"v(t)", None}, {"x(t)",  
    "Phase plot, 50 seconds, case(1)"}}, Epilog -> Disk[{1, 0}, .02],  
 ImageSize -> 400]

4.4.4.2 case (2)

For \(\lambda =0.25\), and \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}=\sqrt{0.5^{2}-0.25^{2}}=0.433\), equations (3,4) become\begin{align} x & =e^{-0.25t}\left ( \cos 0.433t+0.5774\sin 0.433t\right ) \tag{3A}\\ \dot{x} & =-0.05x+e^{-0.25t}\left ( -0.433\sin 0.433t+0.05\cos 0.433t\right ) \tag{4A} \end{align}

Here is the plot generated where the starting point was \(\left ( 1,0\right ) \)

pict

This below is a zoomed in version of the above close to the origin

pict

4.4.4.3 case(3)

For this case, equations (5,6) are used. For \(\lambda =0.5\), equations (5,6) become\begin{align} x & =\left ( 1+0.5t\right ) e^{-0.5t}\tag{5A}\\ \dot{x} & =0.5e^{-0.5t}-0.5\left ( 1+0.5t\right ) e^{-0.5t} \tag{6A} \end{align}

Here is the plot generated, showing starting point \(\left ( 1,0\right ) \) with the code used

pict

lam = 0.5;  
x = (1 + lam*t) Exp[-lam t];  
y = lam*Exp[-lam t] - lam*(1 + lam t) Exp[- lam t]  
ParametricPlot[{x, y}, {t, 0, 30}, Frame -> True,  
 GridLines -> Automatic, GridLinesStyle -> LightGray,  
 FrameLabel -> {{"v(t)", None}, {"x(t)",  
    "Phase plot, 50 seconds, case(3)"}}, Epilog -> Disk[{1, 0}, .02],  
 ImageSize -> 500, PlotRange -> {{-.3, 1.2}, {-.3, .2}},  
 PlotTheme -> "Classic"]

4.4.5 Problem 5

   4.4.5.1 Part (1)
   4.4.5.2 Part(2)
   4.4.5.3 Part(3)

pict

SOLUTION:

Since \(\omega _{0}=\frac{2\pi }{T_{0}}\), then \(\omega _{0}=\frac{2\pi }{1}=2\pi \).

4.4.5.1 Part (1)

For critical damping \(\lambda =\omega _{0}\) and the solution is\begin{align} x\left ( t\right ) & =\left ( A+Bt\right ) e^{-\lambda t}\tag{1}\\ \dot{x}\left ( t\right ) & =Be^{-\lambda t}-\lambda \left ( A+Bt\right ) e^{-\lambda t} \tag{2} \end{align}

Initial conditions are now used to find \(A,B\). At \(t=0\), \(x\left ( 0\right ) =x_{0}=0.1\). From (1) we obtain\[ x_{0}=A \] And since \(\dot{x}\left ( 0\right ) =0\), then from (2)\begin{align*} 0 & =B-\lambda A\\ B & =\lambda A\\ & =\lambda x_{0} \end{align*}

Putting values found for \(A,B\), back into (1) gives\[ x\left ( t\right ) =\left ( x_{0}+\lambda x_{0}t\right ) e^{-\lambda t}\] Since this is critical damping, then \(\lambda =\omega _{0}=2\pi \), hence\[ x\left ( t\right ) =\left ( x_{0}+2\pi x_{0}t\right ) e^{-2\pi t}\] Finally, since \(x_{0}=0.1\) meter, then\[ x\left ( t\right ) =\left ( \frac{1}{10}+\frac{2\pi }{10}t\right ) e^{-2\pi t}\] A plot of the above for \(0\leq t\leq 2s\) is given below

pict

4.4.5.2 Part(2)

For overdamped, \(\lambda >\omega _{0}\) the two roots of the characteristic polynomial are real, hence no oscillation occur. The solution is given by\begin{equation} x\left ( t\right ) =Ae^{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t}+Be^{\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t} \tag{1} \end{equation} \(A,B\,\)are found from initial conditions. When \(t=0\) the above becomes\begin{equation} x_{0}=A+B \tag{2} \end{equation} Taking derivative of (1) gives\[ \dot{x}\left ( t\right ) =A\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) e^{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t}+B\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) e^{\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t}\] At \(t=0\) the above becomes\begin{equation} 0=\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) A+\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) B \tag{3} \end{equation} We have two equations (2,3) which we solve for \(A,B\). From (2), \(A=x_{0}-B\), and (3) becomes\begin{align} 0 & =\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) \left ( x_{0}-B\right ) +\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) B\nonumber \\ 0 & =\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) x_{0}-B\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) +\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) B\nonumber \\ 0 & =\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) x_{0}-2B\sqrt{\lambda ^{2}-\omega _{0}^{2}}\nonumber \\ B & =\frac{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) x_{0}}{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}} \tag{4} \end{align}

Using \(B\) found in (4) then (3) now gives \(A\) as\begin{align*} A & =x_{0}-B\\ & =x_{0}-\frac{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) x_{0}}{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}}\\ & =x_{0}\left ( 1-\frac{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) }{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}}\right ) \\ & =x_{0}\left ( \frac{\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}}{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}}\right ) \end{align*}

Hence the complete solution from (1) becomes\begin{equation} x\left ( t\right ) =x_{0}\left ( \frac{\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}}{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}}\right ) e^{\left ( -\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t}+x_{0}\left ( \frac{-\lambda +\sqrt{\lambda ^{2}-\omega _{0}^{2}}}{2\sqrt{\lambda ^{2}-\omega _{0}^{2}}}\right ) e^{\left ( -\lambda -\sqrt{\lambda ^{2}-\omega _{0}^{2}}\right ) t} \tag{5} \end{equation} The above is now used for each case below to plot the solution..

case (i) \(\lambda =2.2\pi ,\omega _{0}=2\pi ,x_{0}=0.1\), hence (5) becomes\begin{align*} x\left ( t\right ) & =0.1\left ( \frac{2.2\pi +\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -2.2\pi +\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}+0.1\left ( \frac{-2.2\pi +\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -2.2\pi -\sqrt{\left ( 2.2\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}\\ & =0.17e^{-4.032\,2t}-0.07\,e^{-9.791t} \end{align*}

A plot of the above for \(0\leq t\leq 2s\) is given below

pict

case (ii) \(\lambda =4\pi ,\omega _{0}=2\pi ,x_{0}=0.1\), hence (5) becomes\begin{align*} x\left ( t\right ) & =0.1\left ( \frac{4\pi +\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -4\pi +\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}+0.1\left ( \frac{-4\pi +\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -4\pi -\sqrt{\left ( 4\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}\\ & =0.1077e^{-1.6836t}-0.00774e^{-23.449t} \end{align*}

A plot of the above for \(0\leq t\leq 2s\) is given below

pict

case (iii) \(\lambda =10\pi ,\omega _{0}=2\pi ,x_{0}=0.1\), hence (5) becomes\begin{align*} x\left ( t\right ) & =0.1\left ( \frac{10\pi +\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -10\pi +\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}+0.1\left ( \frac{-10\pi +\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}{2\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}}\right ) e^{\left ( -10\pi -\sqrt{\left ( 10\pi \right ) ^{2}-\left ( 2\pi \right ) ^{2}}\right ) t}\\ & =0.101\,e^{-0.634\,73t}-0.001034e^{-62.197t} \end{align*}

A plot of the above for \(0\leq t\leq 2s\) is given below

pict

To compare to the critical damped case, the above three plots are plotted on the same figure against the critical damped case in order to get a better picture and be able to compare the results

pict

From the above we see that critical damping has the fastest decay of the response \(x\left ( t\right ) \). As the damping increases, it takes longer for the response to decay.

4.4.5.3 Part(3)

For the underdamped case, the solution is given by\begin{equation} x\left ( t\right ) =e^{-\lambda t}\left ( A\cos \omega _{d}t+B\sin \omega _{d}t\right ) \tag{1} \end{equation} Where \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}\) and \(A,B\,\) are constant of integration that can be found from initial conditions. And\begin{equation} \dot{x}\left ( t\right ) =-\lambda e^{-\lambda t}\left ( A\cos \omega _{d}t+B\sin \omega _{d}t\right ) +e^{-\lambda t}\left ( -A\omega _{d}\sin \omega _{d}t+B\omega _{d}\cos \omega _{d}t\right ) \tag{2} \end{equation} Applying initial conditions \(x\left ( 0\right ) =x_{0}\) then (1) becomes\[ x_{0}=A \] Applying initial conditions \(\dot{x}\left ( 0\right ) =0\) then (2) becomes\begin{align*} 0 & =-\lambda x_{0}+B\omega _{d}\\ B & =\frac{\lambda x_{0}}{\omega _{d}} \end{align*}

Replacing \(A,B\) back into the solution (1) gives the solution\begin{equation} x\left ( t\right ) =e^{-\lambda t}\left ( x_{0}\cos \omega _{d}t+\frac{\lambda x_{0}}{\omega _{d}}\sin \omega _{d}t\right ) \tag{3} \end{equation} We now use the above solution for the rest of the problem

case(i) \(\lambda =5s^{-1},\omega _{0}=2\pi ,x_{0}=0.1\), hence \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}=\sqrt{\left ( 2\pi \right ) ^{2}-5^{2}}=3.8051\) and (3) becomes\begin{align*} x\left ( t\right ) & =e^{-5t}\left ( 0.1\cos \left ( 3.8051t\right ) +\frac{\left ( 5\right ) \left ( 0.1\right ) }{3.8051}\sin \left ( 3.8051t\right ) \right ) \\ & =e^{-5t}\left ( 0.1\cos \left ( 3.8051t\right ) +0.1314\sin \left ( 3.8051t\right ) \right ) \end{align*}

A plot of the above solution \(x\left ( t\right ) \) for \(0\leq t\leq 2s\) is given below

pict

case(ii) \(\lambda =1s^{-1},\omega _{0}=2\pi ,x_{0}=0.1\), hence \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}=\sqrt{\left ( 2\pi \right ) ^{2}-1^{2}}=6.2031\) and (3) becomes\begin{align*} x\left ( t\right ) & =e^{-t}\left ( 0.1\cos \left ( 6.2031t\right ) +\frac{\left ( 1\right ) \left ( 0.1\right ) }{6.2031}\sin \left ( 6.2031t\right ) \right ) \\ & =e^{-t}\left ( 0.1\cos \left ( 6.2031t\right ) +0.016\sin \left ( 6.2031t\right ) \right ) \end{align*}

A plot of the above solution \(x\left ( t\right ) \) for \(0\leq t\leq 2s\) is given below

pict

case(iii) \(\lambda =0.1s^{-1},\omega _{0}=2\pi ,x_{0}=0.1\), hence \(\omega _{d}=\sqrt{\omega _{0}^{2}-\lambda ^{2}}=\sqrt{\left ( 2\pi \right ) ^{2}-0.1^{2}}=6.2824\) and (3) becomes\begin{align*} x\left ( t\right ) & =e^{-0.1t}\left ( 0.1\cos \left ( 6.2824t\right ) +\frac{\left ( 0.1\right ) \left ( 0.1\right ) }{6.2824}\sin \left ( 6.2824t\right ) \right ) \\ & =e^{-0.1t}\left ( 0.1\cos \left ( 6.2824t\right ) +0.001592\sin \left ( 6.2824t\right ) \right ) \end{align*}

A plot of the above solution \(x\left ( t\right ) \) for \(0\leq t\leq 2s\) is given below

pict

To compare to the critical damped case, the above 3 plots are now plotted on the same figure against the critical damped case in order to get a better picture and be able to compare the results

pict

As the damping becomes smaller, more oscillation occur. The case for \(\lambda =5s^{-1}\) had the smallest oscillation.

4.4.6 HW 4 key solution

pict

pict

pict

pict

pict

pict

pict

pict

pict

pict

pict