Dynamics of Machines - Part II - IFS.pdf

Dynamics of Machines - Part II - IFS.pdf Dynamics of Machines - Part II - IFS.pdf

fam.web.mek.dtu.dk
from fam.web.mek.dtu.dk More from this publisher
28.08.2013 Views

After having created the mechanical model for the physical system, the next step is to derive the equation of motion based on the mechanical model. The mechanical model is composed of a lumped mass m1 (assumption !!!), spring with equivalent stiffness coefficient k1 (calculated using beam theory) and damper with equivalent viscous coefficient d (obtained experimentally). While creating the mechanical model and assuming that the mass is a particle, the equation of motion can be derived, for example, using Newton’s second law: m1¨y1(t) + d1 ˙y1(t) + k1y1(t) = f1(t) (23) ¨y1(t) + d1 ˙y1(t) + m1 k1 y1(t) = m1 f1 (t) (24) m1 ¨y1(t) + 2ξωn ˙y1(t) + ω 2 ny1(t) = f1 (t) (25) m1 1.6.3 Analytical and Numerical Solution of the Equation of Motion After having created the mechanical model (step 1) and derived the mathematical model (step 2) for this mechanical model, the next step is to solve the equation of motion (step 3), aiming at analyzing (step 4) the dynamical behavior of the physical system. Here the analytical and numerical solutions of linear differential equations are presented and compared. Later on you can choose the most convenient way to solve the equations and analyze the dynamical behavior of physical systems. It is important to mention that the numerical procedure is very simple, an integrator of first-order. Other integrators can be used depending on the characteristics of the mechanical models, for example, Runge-Kutta of higher order (third, fourth, etc.) among others. Homogeneous Solution or Transient Solution and Transient Analysis – The homogeneous solution of a linear differential equation is also called transient solution. The homogeneous differential equation is achieved when the right side of the equation is set zero (see eq.(26), or in other words, when no force acts on the system. All analyzes and conclusions obtained from the homogeneous solution are called transient analyzes, and provide information about the dynamical behavior of the system while perturbations of displacement and velocities (in case of second order differential equations of motion) are introduced into the system. ¨y(t) + 2ξωn ˙y(t) + ω 2 ny(t) = 0 (26) yh(t) = Ce λt , (assumption) (27) ˙yh(t) = λCe λt ¨yh(t) = λ 2 Ce λt The assumption (27) and its derivatives are introduced into the differential equation (26), leading to ⇕ λ 2 Ce λt + 2ξωnλCe λt + ω 2 nCe λt = 0 (λ 2 + 2ξωnλ + ω 2 n)Ce λt = 0 (28) 10

Demanding (λ 2 + 2ξωnλ + ω 2 n) = 0, because Ce λt = 0 in eq.(28), one gets two values of λ, or two roots for the equation (λ 2 + 2ξωnλ + ω 2 n) = 0: λ1 = −ξωn − ωn 1 − ξ2 · i λ2 = −ξωn + ωn 1 − ξ2 · i It is important to highlight that all analyzes carried out here will be concentrated in cases where ξ < 1 (sub-critically damped system). Cases where ξ = 1 or ξ > 1 are called critically or super-critically damped systems respectively. In these cases no problem related to amplification of vibration amplitudes can be found while crossing resonances. Using the roots λ1 and λ2 the homogenous solution is: yh(t) = C1e λ1t + C2e λ2t where C1 and C2 are defined as a function of the initial condition of the movement when t = 0: • initial displacement y(0) = yini • initial movement ˙y(0) = vini It is important to point out that these constants have to be calculated by using the general solution, which will be presented later. Permanent Solution and Steady-State Analysis – The permanent solution of a linear differential equation is also called steady-state solution. The complete differential equation is achieved when the right side of the equation is completed with the excitation (see eq.(30)), or in other words, when static or dynamic forces act on the system. All analyzes and conclusions obtained from the permanent solution are called steady-state analyzes, and provide information about the dynamical behavior of the system while excitation forces are introduced into the system. Let us introduce an excitation force which value oscillates in time with frequency ω [rad/s]: ¨y(t) + 2ξωn ˙y(t) + ω 2 ny(t) = f m eiωt (29) (30) yh(t) = Ae iωt , (assumption) (31) ˙yh(t) = jωAe iωt ¨yh(t) = (jω) 2 Ae iωt = −(ω) 2 Ae iωt The assumption (31) and its derivatives are introduced into the differential equation (30), leading to: − (ω) 2 Ae iωt + 2ξωn(jωAe jωt ) + ω 2 nAe iωt = f m eiωt ⇕ (−ω 2 + ω 2 n + i2ξωnω)A = f m 11 (32)

Demanding (λ 2 + 2ξωnλ + ω 2 n) = 0, because Ce λt = 0 in eq.(28), one gets two values <strong>of</strong> λ, or<br />

two roots for the equation (λ 2 + 2ξωnλ + ω 2 n) = 0:<br />

<br />

λ1 = −ξωn − ωn 1 − ξ2 · i<br />

<br />

λ2 = −ξωn + ωn 1 − ξ2 · i<br />

It is important to highlight that all analyzes carried out here will be concentrated<br />

in cases where ξ < 1 (sub-critically damped system). Cases where ξ = 1 or ξ > 1<br />

are called critically or super-critically damped systems respectively. In these cases<br />

no problem related to amplification <strong>of</strong> vibration amplitudes can be found while<br />

crossing resonances. Using the roots λ1 and λ2 the homogenous solution is:<br />

yh(t) = C1e λ1t + C2e λ2t<br />

where C1 and C2 are defined as a function <strong>of</strong> the initial condition <strong>of</strong> the movement when t = 0:<br />

• initial displacement y(0) = yini<br />

• initial movement ˙y(0) = vini<br />

It is important to point out that these constants have to be calculated by using the general<br />

solution, which will be presented later.<br />

Permanent Solution and Steady-State Analysis – The permanent solution <strong>of</strong> a linear<br />

differential equation is also called steady-state solution. The complete differential equation is<br />

achieved when the right side <strong>of</strong> the equation is completed with the excitation (see eq.(30)), or<br />

in other words, when static or dynamic forces act on the system. All analyzes and conclusions<br />

obtained from the permanent solution are called steady-state analyzes, and provide information<br />

about the dynamical behavior <strong>of</strong> the system while excitation forces are introduced into the<br />

system. Let us introduce an excitation force which value oscillates in time with frequency ω<br />

[rad/s]:<br />

¨y(t) + 2ξωn ˙y(t) + ω 2 ny(t) = f<br />

m eiωt<br />

(29)<br />

(30)<br />

yh(t) = Ae iωt , (assumption) (31)<br />

˙yh(t) = jωAe iωt<br />

¨yh(t) = (jω) 2 Ae iωt = −(ω) 2 Ae iωt<br />

The assumption (31) and its derivatives are introduced into the differential equation (30), leading<br />

to:<br />

− (ω) 2 Ae iωt + 2ξωn(jωAe jωt ) + ω 2 nAe iωt = f<br />

m eiωt<br />

⇕<br />

(−ω 2 + ω 2 n + i2ξωnω)A = f<br />

m<br />

11<br />

(32)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!