Dynamics of Machines - Part II - IFS.pdf
Dynamics of Machines - Part II - IFS.pdf Dynamics of Machines - Part II - IFS.pdf
In order to calculate the stiffness matrix (k11, k12, k21, k22) of the 2 d.o.f. system, Ky = f (16) where K = k11 k12 k21 k22 y = { y1 y2 } T f = { f1 f2 } T one has to solve the following system of equations for case (I) and case (II): k11 k12 k21 k22 y1 · y2 = Case (I): f = { F 0 } T k11 k12 k21 k22 · Case (II): f = { 0 F } T k11 k12 k21 k22 f1 f2 F ·L 3 1 3·EI F 6EI · (2L3 1 + 3L2 1 (L2 − L1)) F 6EI · · (2L31 + 3L21 (L2 − L1)) F ·L3 2 3·EI = = F 0 0 F It leads to a system of 4 equations, which can be written in a matrix form: ⎡ ⎢ ⎣ F ·L 3 1 3·EI F 6EI · (2L31 + 3L 2 0 1(L2 − L1)) 0 0 F ·L 0 3 1 3·EI F 6EI · (2L31 + 3L 2 F 6EI 1(L2 − L1)) · (2L31 + 3L 2 1(L2 − L1)) F ·L 3 0 2 3·EI 0 0 F 6EI 0 · (2L31 + 3L 2 1(L2 − L1)) F ·L 3 ⎤ ⎥ ⎦ 2 3·EI · ⎧ ⎪⎨ · ⎪⎩ ⎫ ⎧ ⎪⎬ ⎪⎨ F 0 = ⎪⎭ ⎪⎩ 0 F ⎫ ⎪⎬ ⎪⎭ (21) k11 k12 k21 k22 Solving this matrix system of order 4 by using the software MATHEMATICA, or by using Cramer’s rule, one achieves the stiffness matrix: EI K = (4L2 − L1)(L1 − L2) 2 ⎡ 12(L2/L1) ⎣ 3 ⎤ 6(L1 − 3L2)/L1 ⎦ (22) 6(L1 − 3L2)/L1 12 Similar procedure can be made in order to get the stiffness matrix of the 3 d.o.f system. This is the motivation of an exercise later on. The results (stiffness matrix with 9 stiffness coefficient) will be presented in the section describing mechanical systems with 3 d.o.f. 8 (17) (18) (19) (20)
1.6 Mechanical Systems with 1 D.O.F. 1.6.1 Physical System and Mechanical Model (a) (b) (c) Figure 4: (a) Real mechanical system composed of a turbine attached to an airplane flexible wing; (b) Laboratory prototype built by a lumped mass attached to a flexible beam); (c) Equivalent mechanical model with 1 D.O.F. for a lumped mass attached to a flexible beam. 1.6.2 Mathematical Model It is important to point out, that the equations of motion in Dynamics of Machinery will frequently have the form of second order differential equations: ¨y(t) = F(y(t), ˙y(t)). Such equations can generally be linearized around an operational position of a physical system, leading to second order linear differential equations. It means that the coefficients which are multiplying the variables ¨y(t) , ˙y(t) , y(t) (co-ordinate chosen for describing the motion of the physical system) do not depend on the variables themselves. In the mechanical model presented in figure 4 these coefficients are constants: m1, d1 and k1. One of the aims of the course Dynamics of Machinery is to help the students to properly find these coefficients so that the equations of motion can really describe the movement of the physical system. The coefficients can be predicted using theoretical or experimental approaches. 9
- Page 1 and 2: DYNAMICS OF MACHINES 41614 PART I -
- Page 3 and 4: 1 Introduction to Dynamical Modelli
- Page 5 and 6: 1.3 Data of the Mechanical System
- Page 7: 1.5 Calculating Stiffness Matrices
- Page 11 and 12: Demanding (λ 2 + 2ξωnλ + ω 2 n
- Page 13 and 14: 1 yini − A det λ1 vini − A C2
- Page 15 and 16: 1.6.4 Analytical and Numerical Solu
- Page 17 and 18: (a) y(t) [m] (b) y(t) [m] (c) y(t)
- Page 19 and 20: 1.6.6 Homogeneous Solution or Free-
- Page 21 and 22: (a) Amplitude [m/s 2 ] x Signal 10
- Page 23 and 24: (a) Amplitude [m/s 2 ] (b) Amplitud
- Page 25 and 26: Imag(A(ω)) [m/N] 0 −1 −2 −3
- Page 27 and 28: 1.6.11 Superposition of Transient a
- Page 29 and 30: (a) Amplitude [m/s 2 ] (b) Amplitud
- Page 31 and 32: 1.7 Mechanical Systems with 2 D.O.F
- Page 33 and 34: ⎧ ⎫ ⎪⎨ ˙y1(t) ⎪⎬ ˙y
- Page 35 and 36: zini = U c + A ⇒ c = U −1 {(zin
- Page 37 and 38: 1.7.4 Modal Analysis using Matlab e
- Page 39 and 40: 0.7 0.6 0.5 0.4 0.3 0.2 0.1 First M
- Page 41 and 42: %__________________________________
- Page 43 and 44: ||y 1 (ω)|| [m/N] Phase [ o] Excit
- Page 45 and 46: ||y i (ω)|| [m/N] Phase [ o] 0.8 0
- Page 47 and 48: Imag(y i (ω)/f 1 (ω)) (i=1,2) [m/
- Page 49 and 50: 1.8 Mechanical Systems with 3 D.O.F
- Page 51 and 52: which could be verified using Modal
- Page 53 and 54: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
- Page 55 and 56: 1.8.4 Theoretical Frequency Respons
- Page 57 and 58: (a) Amplitude [m/s 2 ] (b) Amplitud
In order to calculate the stiffness matrix (k11, k12, k21, k22) <strong>of</strong> the 2 d.o.f. system,<br />
Ky = f (16)<br />
where<br />
K =<br />
k11 k12<br />
k21 k22<br />
<br />
y = { y1 y2 } T<br />
f = { f1 f2 } T<br />
one has to solve the following system <strong>of</strong> equations for case (I) and case (<strong>II</strong>):<br />
k11 k12<br />
k21 k22<br />
<br />
y1<br />
·<br />
y2<br />
<br />
=<br />
Case (I): f = { F 0 } T<br />
k11 k12<br />
k21 k22<br />
<br />
·<br />
Case (<strong>II</strong>): f = { 0 F } T<br />
k11 k12<br />
k21 k22<br />
f1<br />
f2<br />
<br />
F ·L 3 1<br />
3·EI<br />
F<br />
6EI · (2L3 1 + 3L2 1 (L2 − L1))<br />
<br />
F<br />
6EI<br />
·<br />
· (2L31 + 3L21 (L2 − L1))<br />
F ·L3 2<br />
3·EI<br />
<br />
<br />
=<br />
=<br />
F<br />
0<br />
0<br />
F<br />
It leads to a system <strong>of</strong> 4 equations, which can be written in a matrix form:<br />
⎡<br />
⎢<br />
⎣<br />
F ·L 3 1<br />
3·EI<br />
F<br />
6EI · (2L31 + 3L 2 0<br />
1(L2 − L1))<br />
0<br />
0<br />
F ·L<br />
0<br />
3 1<br />
3·EI<br />
F<br />
6EI · (2L31 + 3L 2 F<br />
6EI<br />
1(L2 − L1))<br />
· (2L31 + 3L 2 1(L2 − L1))<br />
F ·L 3 0<br />
2<br />
3·EI<br />
0<br />
0<br />
F<br />
6EI<br />
0<br />
· (2L31 + 3L 2 1(L2 − L1))<br />
F ·L 3 ⎤<br />
⎥<br />
⎦<br />
2<br />
3·EI<br />
·<br />
⎧<br />
⎪⎨<br />
·<br />
⎪⎩<br />
⎫ ⎧<br />
⎪⎬ ⎪⎨<br />
F<br />
0<br />
=<br />
⎪⎭ ⎪⎩<br />
0<br />
F<br />
⎫<br />
⎪⎬<br />
⎪⎭<br />
(21)<br />
k11<br />
k12<br />
k21<br />
k22<br />
Solving this matrix system <strong>of</strong> order 4 by using the s<strong>of</strong>tware MATHEMATICA, or by using<br />
Cramer’s rule, one achieves the stiffness matrix:<br />
EI<br />
K =<br />
(4L2 − L1)(L1 − L2) 2<br />
⎡<br />
12(L2/L1)<br />
⎣<br />
3 ⎤<br />
6(L1 − 3L2)/L1<br />
⎦ (22)<br />
6(L1 − 3L2)/L1 12<br />
Similar procedure can be made in order to get the stiffness matrix <strong>of</strong> the 3 d.o.f system. This is<br />
the motivation <strong>of</strong> an exercise later on. The results (stiffness matrix with 9 stiffness coefficient)<br />
will be presented in the section describing mechanical systems with 3 d.o.f.<br />
8<br />
<br />
<br />
(17)<br />
(18)<br />
(19)<br />
(20)