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

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

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)

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

Saved successfully!

Ooh no, something went wrong!