Prikazani su postovi s oznakom Mass. Prikaži sve postove
Prikazani su postovi s oznakom Mass. Prikaži sve postove

Undaped Forced Vibration



Let’s consider a system which consists of block with mass m and spring with stiffness k. On the block acts a force F and this force can be represented by the following equation.
$$F={{F}_{0}}\sin {{\omega }_{0}}t$$
This is a periodic force which has an amplitude of F and forcing frequency ω. The free body diagram for the block when it’s displaced a distance x is shown in next figure.

By applying the second Newton’s law we have:
$$\begin{align} & \sum{{{F}_{x}}=m{{a}_{x}}} \\ & {{F}_{0}}\sin {{\omega }_{0}}t-kx=m\ddot{x}, \\ & m\ddot{x}+kx={{F}_{0}}\sin {{\omega }_{0}}t \\ \end{align}$$
If we divide the previous differential equation with m we will get:
$$\ddot{x}+\frac{k}{m}x=\frac{{{F}_{0}}}{m}\sin {{\omega }_{0}}t$$
This equation is a nonhomogeneous-second order differential equation. The general solution consists of a homogenous and particular solution.
We will get homogenous solution of the differential equation by assuming solution of the following differential equation:
$$\ddot{x}+\omega _{n}^{2}x=0$$
In the previous differential equation the ω is the natural frequency of the system and it can be determined by following expression:
$${{\omega }_{n}}=\sqrt{\frac{k}{m}}$$
The assumed solution for the homogenous part of differential equation is:
$$x=A\sin {{\omega }_{n}}t+B\cos {{\omega }_{n}}t,$$
Now we need to prove that the assumed solution is really solution of the differential equation but before that we need to derive the first and second derivation of the assumed solution:
$$\begin{align} & x=A\sin {{\omega }_{n}}t+B\cos {{\omega }_{n}}t, \\ & \dot{x}=A{{\omega }_{n}}\cos {{\omega }_{n}}t-B{{\omega }_{n}}\sin {{\omega }_{n}}t \\ & \ddot{x}=-A\omega _{n}^{2}\sin {{\omega }_{n}}t-B\omega _{n}^{2}\cos {{\omega }_{n}}t=-\omega _{n}^{2}x \\ \end{align}$$
The next step is to substitute these values of the assumed solution into differential equation. Whit this step we will prove if the assumption is correct or wrong.
$$\begin{align} & -A\omega _{n}^{2}\sin {{\omega }_{n}}t-B\omega _{n}^{2}\cos {{\omega }_{n}}t+\omega _{n}^{2}\left( A\sin {{\omega }_{n}}t+B\cos {{\omega }_{n}}t \right)=0 \\ & -\omega _{n}^{2}x+\omega _{n}^{2}x=0 \\ & 0=0 \\ \end{align}$$
So as you can see we have proved that previous assumed solution is the solution of the homogenous part of the differential equation. All that is left is to assume the particular solution of nonhomogenous differential equation.
The homogenous part of the solution can be written in the following form:
$${{x}_{h}}=Csin\left( {{\omega }_{n}}t+\phi \right)$$
Since the motion is period, the particular solution can be determined by assuming a solution in the following form:
$${{x}_{p}}=X\sin {{\omega }_{0}}t$$
where X is constant. Taking the second time derivate and substituting into the nonhomogenous differential equation:
$$\begin{align} & x=X\sin {{\omega }_{0}}t \\ & \dot{x}=X{{\omega }_{0}}\cos {{\omega }_{0}}t \\ & \ddot{x}=-\omega _{0}^{2}X\sin {{\omega }_{0}}t \\ \end{align}$$
Now we can prove that the assumed solution of the nonhomogenous differential equation:
$$\begin{align} & \ddot{x}+\frac{k}{m}x=\frac{{{F}_{0}}}{m}\sin {{\omega }_{0}}t \\ & -\omega _{0}^{2}X\sin {{\omega }_{0}}t+\frac{k}{m}\left( X\sin {{\omega }_{0}}t \right)=\frac{{{F}_{0}}}{m}\sin {{\omega }_{0}}t \\ & -\omega _{0}^{2}X+\frac{k}{m}X=\frac{{{F}_{0}}}{m} \\ & X\left( \frac{k}{m}-\omega _{0}^{2} \right)=\frac{{{F}_{0}}}{m} \\ & X=\frac{{{F}_{0}}/m}{(k/m)-\omega _{0}^{2}}=\frac{{{F}_{0}}/k}{1-{{\left( {{\omega }_{0}}/{{\omega }_{n}} \right)}^{2}}} \\ \end{align}$$
Substituting expression for X into particular solution we will obtain:
$$x=\frac{{{F}_{0}}/k}{1-{{\left( {{\omega }_{0}}/{{\omega }_{n}} \right)}^{2}}}\sin {{\omega }_{0}}t$$
Now that we obtained homogenous and particular solution of the nonhomogenous differential equation we can write the general solution:
$$\begin{align} & x={{x}_{h}}+{{x}_{p}}, \\ & x=Csin\left( {{\omega }_{n}}t+\phi \right)+\frac{{{F}_{0}}/k}{1-{{\left( {{\omega }_{0}}/{{\omega }_{n}} \right)}^{2}}}\sin {{\omega }_{0}}t \\ \end{align}$$
The homogenous solution defines the FREE VIBRATIONS which depends on the natural frequency of the system and constants C and ϕ. The particular solution describes the forced vibration of the block caused by the applied force
$$F={{F}_{0}}\sin {{\omega }_{0}}t$$
 
Since all the vibrating systems are subjected to the friction, the free vibration will dampen out after some time. Because of free vibration dependence on time we can call them transient and forced vibration are steady –state. The reason why we call forced vibrations steady state is because they are the only vibrations that remain. As long as force which causes vibrations of the system is present in the system these vibrations will be active.


Electric Circuit Analogs



The vibrating system and its characteristics can be represented by an electric circuit. Now let’s consider the electrical circuit which is shown in the next figure.
 
Figure 1 - Electric Circuit
The system consists of inductor L, a resistor R and a capacitor C. When the voltage E(t) is applied to the system, it caused a current of magnitude i to flow through the circuit. As the current flows past the inductor the voltage drop is L(di/dt). When the current flows through the resistor the voltage drop is Ri, and when it arrives at the capacitor the drop is 
$$\frac{1}{C}\int{idt}$$
Since the current cannot flow past a capacitor, it’s only possible to measure the charge q acting on the capacitor. The charge can be related to the current by the equation
$$i=\frac{dq}{dt}$$
Thus the voltage drops, which occur across the inductor, resistor, and capacitor becomes
$$L\frac{{{d}^{2}}q}{d{{t}^{2}}},R\frac{dq}{dt},\frac{q}{C}$$
By applying the Kirchhoff’s law which states that the applied voltage balances the sum of the voltage drops around the circuit. Therefore:
$$L\frac{{{d}^{2}}q}{d{{t}^{2}}}+R\frac{dq}{dt}+\frac{1}{C}q=E(t)$$
As you can see the previous differential equation which represents the change of voltage in electrical system is similar to the differential equation which describes the motion of Viscous Dampers Forced Vibration system.
By comparing these to differential equation we can see that these equation have the same form, and hence mathematically the procedure of analyzing an electric circuit is the same as that of analyzing a vibration mechanical system
$$m\frac{{{d}^{2}}x}{d{{t}^{2}}}+c\frac{dx}{dt}+kx=F(t)$$
From this we can derive the analogs between two equations and we will show them in the following table.
Electrical

Mechanical

Electric Charge
q
Displacement
x
Electric Current
i
Velocity
dx/dt
Voltage
E(t)
Applied Force
F(t)
Inductance
L
Mass
m
Resistance
R
Viscous damping coefficient
c
Reciprocal of capacitance
1/C
Spring stiffness
k

Energy Methods



As we all know the simple harmonic motion of a body is due only to gravitational and elastic restoring forces acting on the body. Since these forces are conservative, it is also possible to use the conservation of energy equation in order to obtain the body’s natural frequency or period of vibration.
Now let’s consider a system with the block and spring mode. 
 

Figure 1- Undamped Free Vibrational System


When the block is displaced at distance x from the equilibrium position, the kinetic energy: 
$$T=\frac{1}{2}m{{v}^{2}}=\frac{1}{2}m{{\dot{x}}^{2}}$$

and the potential energy is:
$$V=\frac{1}{2}k{{x}^{2}}.$$

Since the energy is conserved, we can write: 
$$\begin{align} & T+V=\text{constant} \\ & \frac{1}{2}m{{{\dot{x}}}^{2}}+\frac{1}{2}k{{x}^{2}}=\text{constant} \\ \end{align}$$

To get the differential equation that describes the motion of the system we need to derivate the equation and than we get: 
$$\begin{align} & m\dot{x}\ddot{x}+kx\dot{x}=0 \\ & \dot{x}\left( m\ddot{x}+kx \right)=0 \\ \end{align}$$

Analyzing the previous equation we can conclude that the speed or x’ is not always equals zero so the general assumption would be:
$$m\ddot{x}+kx=0$$

As you can see we have got the differential equation which describes the vibrations of the system. If the conservation of energy equation is written for a system of connected bodies, the natural frequency or the equation of motion can also be determined by time differentiation. It is not necessary to dismember the system to account for the internal forces because they do no work.

The Study of Vibrations - Examples



Example 1.1.

A spring is stretched 100 mm by a 10 kg block. If the block is displaced 50 mm downward from its equilibrium position and given a downward velocity of 3 m/s determine the differential equation which describes the motion. Assume that positive displacement is downward. Also, determine the position of the block when t=4 s. 

Figure 1.1 – a) Model of a system, b) Free body diagrams for external and effective forces

Solution: From previous figure we can derive the differential equation.
$$\begin{align} & +\downarrow \sum{{{F}_{y}}}=m{{a}_{y}} \\ & mg-k\left( y+{{y}_{st}} \right)=m\ddot{y} \\ \end{align}$$
where 
 $$k{{y}_{st}}=mg$$
Now we need to substitute the previous equation into the differential equation. With this operation we get: 
$$\begin{align} & mg-ky-mg=m\ddot{y}, \\ & \ddot{y}+\frac{k}{m}y=0, \\ \end{align}$$


Before we can determine  the natural frequency of the system we need to determine the stiffness of the spring. We know that the force the stretch the spring is proportion to the stiffness of the spring and the elongation of the spring. In this case we don’t have a force which acts on spring but we have a block which is attached to the sprig. So we will determine the stiffness of the spring by following expression
$$\begin{align} & F=ky,F=G=mg \\ & k=\frac{G}{y}=\frac{mg}{y}=\frac{10\cdot 9.81}{0.1}=981\text{ N/m}{{\text{m}}^{2}} \\ \end{align}$$
The next step is to determine the natural frequency and we will do that by following expression:
$$\omega =\sqrt{\frac{k}{m}}=\sqrt{\frac{981}{10}}=9.9045\text{ }{{\text{s}}^{-1}}$$
 
When we have determined the natural frequency we can insert it into differential equation. Whit this operation we get:
$$\begin{align} & \ddot{y}+{{(9.9045)}^{2}}y=0, \\ & \ddot{y}+98.1y=0, \\ \end{align}$$
All that is left is to apply the boundary conditions. At
$$t=0,y=0,05$$
By applying this boundary conditions we get:
$$\begin{align} & 0.05=A\sin 0+b\cos 0 \\ & b=0.05 \\ \end{align}$$
At
$$t=0,{{v}_{0}}=3m/s$$ $$\begin{align} & {{v}_{0}}=A\omega \cos 0-0=A \\ & A=\frac{{{v}_{0}}}{\omega }=\frac{3}{9.9045}=0.30289 \\ \end{align}$$
Example 1.2.
A 3 kg block is suspended from a spring having a stiffness of k=200N/m. If the block is pushed 50 mm upward from its equilibrium position and then released from rest, determine the equation that describes the motion. What are the amplitude and the frequency of the vibration? Assume that the positive displacement is downward.
From the given data we can determine the frequency of the vibrations.
$$\begin{align} & \omega =\sqrt{\frac{k}{m}}=\sqrt{\frac{200}{3}}=8.165 \\ & f=\frac{\omega }{2\pi }=\frac{8.165}{2\pi }=1.299=1.3Hz \\ & x=A\sin \omega t+B\cos \omega t, \\ \end{align}$$
Boundary conditions are:
$$\begin{align} & t=0,x=-0.05, \\ & t=0,v=0 \\ \end{align}$$
By applying the boundary condition to the solution of the differential equation we get:
$$\begin{align} & 0.05=A\sin 0+B\cos 0, \\ & b=-0.05. \\ & v=A\omega \cos \omega t-B\omega \sin \omega t, \\ & 0=A\omega \cos 0-B\omega \sin 0 \\ & A=0 \\ \end{align}$$
In order to determine the amplitude of the system we will apply the following equation.
$$C=\sqrt{{{A}^{2}}+{{B}^{2}}}=\sqrt{2.5\cdot {{10}^{-3}}}=0.05m=50mm$$

Example 1.3
A spring has a stiffness of 800 N/m. If a 2 kg block is attached to the spring, pushed 50 mm above its equilibrium position, and released from rest, determine the equation that describes the block’s motion. Assume that positive displacement is downward.
$$\begin{align} & \omega =\sqrt{\frac{k}{m}}=\sqrt{\frac{800}{2}}=20{{s}^{-1}}, \\ & y=A\sin \omega t+B\cos \omega t, \\ \end{align}$$
For t=0, y=-0.05 m
$$\begin{align} & -0.05=A\sin 0+B\cos 0, \\ & B=-0.05 \\ \end{align}$$ $$\begin{align} & y=A\sin \omega t+B\cos \omega t, \\ & v=A\omega \cos \omega t-B\omega \sin \omega t, \\ \end{align}$$
For t=0, v=0
$$\begin{align} & v=A\omega \cos \omega t-B\omega \sin \omega t, \\ & 0=A\omega \cos 0-B\omega \sin 0, \\ & A=0 \\ \end{align}$$
Thus we get:
$$y=-0.05\cos \left( 20t \right)$$
Example 1.4.
A spring is stretched 200 mm by a 15 kg block. If the bloc is displaced 100 mm downward from its equilibrium position and given downward velocity of 0.75 m/s determine the equation which describes the motion. What is a phase angle? Assume that positive displacement is downward
$$\begin{align} & F=ky\Rightarrow k=\frac{F}{y}=\frac{15\cdot 9.81}{0.2}=735.75N/m \\ & \omega =\sqrt{\frac{k}{m}}=\sqrt{\frac{735.75}{15}}=7{{s}^{-1}} \\ & y=A\sin \omega t+B\cos \omega t \\ \end{align}$$
When t=0 then y=0.1 m and
When t=0 then v=0.75 m/s
$$\begin{align} & y=A\sin \omega t+B\cos \omega t \\ & 0.1=A\sin 0+B\cos 0 \\ & B=0.1m \\ \end{align}$$ $$\begin{align} & v=A\omega \cos \omega t-B\omega \sin \omega t \\ & 0.75=A\omega \cos 0-B\omega \sin 0 \\ & 0.75=A\omega \\ & A=\frac{0.75}{7}=0.107 \\ \end{align}$$
Now we can insert the values of A and B into the solution.
$$\begin{align} & y=0.107\sin (7t)+0.100\cos (7t) \\ & \phi =arctan\left( \frac{B}{A} \right)=\arctan \left( \frac{0.100}{0.107} \right)={{43.0}^{\circ }} \\ \end{align}$$