Solved Problem on Harmonic Oscillations
Português
English
advertisement
A block of mass m is attached to a spring, the spring constant is equal to k and a shock
absorber of damping coefficient b. The block is moved from its equilibrium position O to
a point x0 and released from rest. Determine:
a) The differential equation of motion;
b) The equation solution for the system in the case of underdamped oscillations and the angular frequency of the oscillations.
a) The differential equation of motion;
b) The equation solution for the system in the case of underdamped oscillations and the angular frequency of the oscillations.
Problem data:
- Mass of body: m;
- Spring constant: k;
- Damping coefficient: b;
- Initial position (t = 0): x0;
- Initial speed (t = 0): v0 = 0.
Problem diagram:
We choose a reference frame with a positive direction to the right. The block is moved to position
x0 and released with an initial velocity equal to zero. When released, the spring
force will act to restore the equilibrium position (Figure 1). With this, we write the
Initial Conditions of the problem
\[
\begin{align}
& x(0)=x_{0}\\[10pt]
& v_{0}=\frac{dx(0)}{dt}=0
\end{align}
\]
Solution
a) Applying Newton's Second Law (Figure 1)
\[
\begin{gather}
\bbox[#99CCFF,10px]
{F=m\frac{d^{2}x}{dt^{2}}} \tag{I}
\end{gather}
\]
the forces acting on the block are the spring force
\( {\vec{F}}_{S} \)
and the damping force
\( {\vec{F}}_{D} \)
given, in magnitude, by
\[
\begin{gather}
\bbox[#99CCFF,10px]
{F_{S}=-kx} \tag{II-a}
\end{gather}
\]
\[
\begin{gather}
\bbox[#99CCFF,10px]
{F_{D}=-bv=-b\frac{dx}{dt}} \tag{II-b}
\end{gather}
\]
the minus sign in the spring force means that it acts in the opoosite direction of the displacement
of the block (acts in the direction of restoring equilibrium), in the damping force means that it acts
in the opposite direction of the velocity (acts in the direction of braking the movement).
Substituting the expressions (II-a) and (II-b) into the expression (I)
\[
\begin{gather}
-kx-b\frac{dx}{dt}=m\frac{d^{2}x}{dt^{2}}\\[5pt]
m\frac{d^{2}x}{dt^{2}}+b\frac{dx}{dt}+kx=0
\end{gather}
\]
this is a Second Order Homogeneous Differential Equation. Dividing the equation by the mass m
\[
\begin{gather}
\bbox[#FFCCCC,10px]
{\frac{d^{2}x}{dt^{2}}+\frac{b}{m}\frac{dx}{dt}+\frac{k}{m}x=0}
\end{gather}
\]
b) In the equation of the previous item we will make the following definitions
\[
\begin{gather}
2\gamma \equiv \frac{b}{m} \tag{III-a}\\[10pt]
\omega_{0}^{2}\equiv\frac{k}{m} \tag{III-b}
\end{gather}
\]
\[
\begin{gather}
\frac{d^{2}x}{dt^{2}}+2\gamma \frac{dx}{dt}+\omega_{0}^{2}x=0
\end{gather}
\]
Solution of \( \displaystyle \frac{d^{2}x}{dt^{2}}+2\gamma \frac{dx}{dt}+\omega_{0}^{2}x=0 \)
The solution to this type of equation is found substituting
The angular frequency ω is given by
The solution to this type of equation is found substituting
\[
\begin{align}
& x=\operatorname{e}^{\lambda t}\\[10pt]
& \frac{dx}{dt}=\lambda \operatorname{e}^{\lambda t}\\[10pt]
& \frac{d^{2}x}{dt^{2}}=\lambda ^{2}\operatorname{e}^{\lambda t}
\end{align}
\]
substituting these values into the differential equation
\[
\begin{gather}
\lambda ^{2}\operatorname{e}^{\lambda t}+2\gamma \lambda\operatorname{e}^{\lambda t}+\omega_{0}^{2}\operatorname{e}^{\lambda t}=0\\[5pt]
\operatorname{e}^{\lambda t}\left(\lambda ^{2}+2\gamma \lambda+\omega_{0}^{2}\right)=0\\[5pt]
\lambda ^{2}+2\gamma \lambda +\omega_{0}^{2}=\frac{0}{\operatorname{e}^{\lambda t}}\\[5pt]
\lambda ^{2}+2\gamma\lambda +\omega_{0}^{2}=0
\end{gather}
\]
this is the Characteristic Equation that has a solution
\[
\begin{gather}
\Delta =b^{2}-4ac=\left(2\gamma \right)^{2}-4\times 1\times \omega_{0}^{2}=4\gamma^{2}-4\omega_{0}^{2}=4\left(\gamma^{2}-\omega_{0}^{2}\right)\\[10pt]
\lambda_{1}=\frac{-b+\sqrt{\Delta\;}}{2a}=\frac{-2\gamma +\sqrt{4\left(\gamma^{2}-\omega_{0}^{2}\right)\;}}{2\times 1}=-{\frac{2\gamma}{2}}+\frac{2\sqrt{\gamma^{2}-\omega_{0}^{2}\;}}{2}=-\gamma +\sqrt{\gamma^{2}-\omega_{0}^{2}\;}\\[5pt]
\lambda_{2}=\frac{-b-\sqrt{\Delta\;}}{2a}=\frac{-2\gamma -\sqrt{4\left(\gamma^{2}-\omega_{0}^{2}\right)\;}}{2\times 1}=-{\frac{2\gamma}{2}}-\frac{2\sqrt{\gamma^{2}-\omega_{0}^{2}\;}}{2}=-\gamma -\sqrt{\gamma^{2}-\omega_{0}^{2}\;}
\end{gather}
\]
For the system to oscillate with underdamped oscillations we must have
ω02>γ2, the term in the square root will be
\[
\begin{gather}
\sqrt{-1\times \left(\omega_{0}^{2}-\gamma^{2}\right)\;}=\sqrt{-1}\times \sqrt{\left(\omega_{0}^{2}-\gamma^{2}\right)\;}=\mathsf{i}\sqrt{\left(\omega_{0}^{2}-\gamma^{2}\right)\;}
\end{gather}
\]
where \( \mathsf{i}=\sqrt{-1\;} \).The angular frequency ω is given by
\[
\begin{gather}
\omega =\sqrt{\omega_{0}^{2}-\gamma^{2}\;}
\end{gather}
\]
Using the definitions made in (III-a) and (III-b) for ω02 and γ
\[
\begin{gather}
\bbox[#FFCCCC,10px]
{\omega =\sqrt{\frac{k}{m}-\left(\frac{b}{2m}\right)^{2}\;}}
\end{gather}
\]
The solution of the differential equation will be
Differentiation of the expression (IV) with respect to time
\[
\begin{gather}
x=C_{1}\operatorname{e}^{\lambda_{1}t}+C_{2}\operatorname{e}^{\lambda_{2}t}\\[5pt]
x=C_{1}\operatorname{e}^{\left(-\gamma +\mathsf{i}\sqrt{\omega_{0}^{2}-\gamma^{2}\;}\right)t}+C_{2}\operatorname{e}^{\left(-\gamma-\mathsf{i}\sqrt{\omega_{0}^{2}-\gamma^{2}\;}\right)t}\\[5pt]
x=C_{1}\operatorname{e}^{\left(-\gamma t+\mathsf{i}\sqrt{\omega_{0}^{2}-\gamma^{2}\;}t\right)}+C_{2}\operatorname{e}^{\left(-\gamma t-\mathsf{i}\sqrt{\omega_{0}^{2}-\gamma^{2}\;}t\right)}\\[5pt]
x=C_{1}\operatorname{e}^{-\gamma t}\operatorname{e}^{\mathsf{i}\omega \;t}+C_{2}\operatorname{e}^{-\gamma t}\operatorname{e}^{-\mathsf{i}\omega t}\\[5pt]
x=\operatorname{e}^{-\gamma t}\left(C_{1}\operatorname{e}^{\mathsf{i}\omega t}+C_{2}\operatorname{e}^{-\mathsf{i}\omega t}\right)
\end{gather}
\]
where C1 and C2 are constants of integration, using
Euler's Formula
\( \operatorname{e}^{\mathsf{i}\theta }=\cos \theta +\mathsf{i}\sin \theta \)
\[
\begin{gather}
x=\operatorname{e}^{-\gamma t}\left[C_{1}\left(\cos \omega t+\mathsf{i}\sin \omega t\right)+C_{2}\left(\cos \omega t-\mathsf{i}\sin \omega t\right)\right]\\[5pt]
x=\operatorname{e}^{-\gamma t}\left(C_{1}\cos \omega t+\mathsf{i}C_{1}\sin \omega t+C_{2}\cos \omega t-\mathsf{i}C_{2}\sin \omega t\right)\\[5pt]
x=\operatorname{e}^{-\gamma t}\left[\left(C_{1}+C_{2}\right)\cos \omega t+\mathsf{i}\left(C_{1}-C_{2}\right)\sin \omega t\right]
\end{gather}
\]
defining two new constants α and β in terms of C1 and
C2
\[
\begin{gather}
\alpha \equiv C_{1}+C_{2}\\[5pt]
\text{e}\\[5pt]
\beta \equiv \mathsf{i}(C_{1}-C_{2})
\end{gather}
\]
\[
\begin{gather}
x=\operatorname{e}^{-\gamma t}\left(\alpha \cos \omega t+\beta\sin \omega t\right) \tag{IV}
\end{gather}
\]
where α and β are constants of integration determined by the
Initial Conditions.
Differentiation of the expression (IV) with respect to time
\[
\begin{gather}
x=\underbrace{\operatorname{e}^{-\gamma t}}_{u}\underbrace{\left(\alpha\cos \omega t+\beta \sin \omega t\right)}_{v}
\end{gather}
\]
using the Product Rule for the differentiation of functions
\[
\begin{gather}
(uv)'=u'v+uv'
\end{gather}
\]
where
\( u=\operatorname{e}^{-\gamma t} \)
and
\( v=\left(\alpha \cos \omega t+\beta \sin \omega t\right) \), the term in
parentheses is a sum of functions, the derivative is given by the sum of the derivatives
\[
\begin{gather}
(f+g)'=f'+g'
\end{gather}
\]
and the functions sine and cosine are composite functions, using the Chain Rule
\[
\begin{gather}
\frac{dv[w(t)]}{dt}=\frac{dv}{dw}\frac{dw}{dt}
\end{gather}
\]
with
\( v_{1}=\alpha\cos w \),
\( v_{2}=\beta\sin w \)
and
\( w=\omega t \)
\[
\begin{gather}
\frac{dx}{dt}=\frac{du}{dt}v+u\frac{dv}{dt}\\[5pt]
\frac{dx}{dt}=\frac{du}{dt}v+u\left(\frac{df}{dt}+\frac{dg}{dt}\right)\\[5pt]
\frac{dx}{dt}=\frac{du}{dt}v+u\left(\frac{dv_{1}}{dw}\frac{dw}{dt}+\frac{dv_{2}}{dw}\frac{dw}{dt}\right)\\[5pt]
\frac{dx}{dt}=\frac{d\left(\operatorname{e}^{-\gamma t} \right)}{dt}\left(\alpha \cos \omega t+\beta \sin \omega t\right)+\left(\operatorname{e}^{-\gamma t} \right)\left[\frac{d(\alpha\cos w)}{dw}\frac{d(\omega t)}{dt}+\frac{d(\beta\sin w)}{dw}\frac{d(\omega t)}{dt}\right]\\[5pt]
\frac{dx}{dt}=-\gamma \operatorname{e}^{-\gamma t}\left(\alpha \cos \omega t+\beta \sin \omega t\right)+\operatorname{e}^{-\gamma t}\left[(-\alpha\sin w)(\omega)+(\beta\cos w)(\omega)\right]\\[5pt]
\frac{dx}{dt}=-\gamma \operatorname{e}^{-\gamma t}\left(\alpha \cos \omega t+\beta \sin \omega t\right)+\operatorname{e}^{-\gamma t}\left(-\omega\alpha\sin \omega t+\omega\beta\cos \omega t\right)\\[5pt]
\frac{dx}{dt}=\operatorname{e}^{-\gamma t}\left(-\gamma\alpha \cos \omega t-\gamma\beta \sin \omega t-\omega\alpha\sin \omega t+\omega\beta\cos \omega t\right)\\[5pt]
\frac{dx}{dt}=\operatorname{e}^{-\gamma t}\left[-\alpha \left(\gamma\cos \omega t+\omega\sin \omega t\right)-\beta \left(\gamma\sin \omega t-\omega\cos \omega t\right)\right] \tag{V}
\end{gather}
\]
Substituting Initial Conditions into expressions (IV) amd (V)
\[
\begin{gather}
x(0)=x_{0}=\operatorname{e}^{-\gamma \times 0}\left(\alpha \cancelto{1}{\cos \omega \times 0}+\beta \cancelto{0}{\sin \omega \times 0}\right)\\[5pt]
\alpha=x_{0} \tag{VI}
\end{gather}
\]
\[
\begin{gather}
\frac{dx(0)}{dt}=0=\operatorname{e}^{-\gamma \times 0}\left[-\alpha \left(\gamma \cancelto{1}{\cos \omega \times 0}+\omega \cancelto{0}{\sin \omega \times 0}\right)- \right.\\
\qquad\qquad\qquad\qquad\qquad\qquad-\left. \beta \left(\gamma \cancelto{0}{\sin \omega \times 0}-\omega \cancelto{1}{\cos \omega \times 0}\right)\right]\\[5pt]
0=1\times \left[-x_{0} \gamma+\beta \omega \right]\\[5pt]
-x_{0}\gamma +\beta \omega =0\\[5pt]
\beta=\frac{x_{0}\gamma}{\omega} \tag{VII}
\end{gather}
\]
substituting the constants (VI) and (VII) into expression (IV)
\[
\begin{gather}
x=\operatorname{e}^{-\gamma t}\left(x_{0}\cos \omega t+\frac{x_{0}\gamma}{\omega}\sin \omega t\right)
\end{gather}
\]
\[
\begin{gather}
\bbox[#FFCCCC,10px]
{x(t)=x_{0}\operatorname{e}^{-\gamma t}\left(\cos \omega t+\frac{\gamma}{\omega}\sin \omega t\right)}
\end{gather}
\]
advertisement
Fisicaexe - Physics Solved Problems by Elcio Brandani Mondadori is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License .