**Przeglądaj wersję html pliku:**

### MultiDOF

CHAPTER 5: Multi degree-of-freedom vibration.

Initially we will just increase the number of degrees-of-freedom to three. Thereafter we will increase the numbers of degrees-of-freedom. It will soon become apparent that the maths will be overwhelming for systems with many degrees-of-freedom. For this reason computers are used extensively in vibration analysis. However it is important to have a good understanding of vibrating systems even when they have many degrees-of-freedom. The results from solving the equations of motion will be given but the mathematical detail will be omitted. In the next chapter more detail will be given of the mathematical solutions of the equations. At this stage it is the object to aim for understanding of the results rather than the means used to obtain them. 5.1 Three degrees-of-freedom There are numerous examples that involve 3 degrees-of-freedom. The one chosen is an extension of the beam problem considered in chapter 3. Consider the beam mounted on springs at each end and an additional spring/mass system mounted on the beam as shown in figure 5.1.

Figure 5.1 System with 3 degrees of freedom We need to introduce another variable x 3 to locate the position of the mass m 3 . The coordinate x 4 is introduced to aid in the deriving of the equations of motion but can be found from x1 and x 2 so does not imply a fourth degree-of-freedom. From geometry, c c c x 4 = x1 + ( x2 − x1 ) = x1 1 − + x2 ................................ (5.1) L L L The position of the centre of mass (x) is given by, x = (x1 + x2 ) / 2 ...................................................... (5.2) and for small angles, θ = (x1 − x2 ) / L ..................................................... (5.3)

Contents

Free body diagrams ignoring gravity forces and working relative to the static equilibrium position are shown in figure 5.2.

Figure 5.2 Free body diagrams. Using Newton's second law for the linear motion of the centre of mass gives mx' ' = − k1x1 − k2x 2 − k 3(x 4 − x3 ) + F ..................................... (5.4) and the angular motion about the centre of mass gives 1 mL2θ' ' = −k1x1L / 2 + k2 x2L / 2 + k 3(x4 − x3 )(c − L / 2) + F( L / 2 − b) .............. (5.5) 12 and for the mass m3 , m 3x3' ' = k 3(x4 − x3 ) .................................................... (5.6) If in equations (5.4) - (5.6) we substitute for x from (5.2), θ from(5.3) and x 4 from (5.1) we obtain 3 equations with just the three degrees-of-freedom x1 , x 2 and x 3 . Thus after much maths, c c m(x1' ' + x2' ' ) = 2 x1 −k1 − k3 1 − + x2 − k2 − k3 + k 3x3 + F .............. (5.7) L L and m(x1' ' − x2' ' ) = k k c 1 1 b c c 1 c c 1 12 x1 − 1 + k 3 1− − + x2 2 + k 3 − + x3 − k3 − + F − L 2 2 L 2 2 L L 2 L L 2 ...... (5.8) Adding (5.7) and (5.8) mx1' ' = c c c c c 6b x1 −4k1 − k3 1− 4 − 6 + x2 2k 2 − k3 4 − 6 + x3k 3 4 − 6 + F 4 − L L L L L L ...... (5.9) Subtracting (5.7) and (5.8)

Contents

mx 2' ' = c c c c c 6b x1 2k1 + k3 1− 2 − 6 + x2 −4k 2 + k3 2 − 6 − x3k 3 2 − 6 − F 2 − L L L L L L ..... (5.10) and c c m 3x3' ' = k 3 x1 1 − + x2 − x3 L L ........................................ (5.11)

Equations (5.9) - (5.11) may be written in matrix form c c 4k1 + k 3 1 − 4 − 6 L L 0 x1' ' c c 0 x2 '' + −2k1 − k3 1− 2 − 6 L L m3 x 3' ' c −k 3 1− L 6b F 4 − L 6b = −F 2 − L 0 .. (5.12) However the stiffness matrix is not symmetrical. If we had used the Lagrange equations the derivation would have been as follows, T= 1 1 1 mx2 + IG θ2 + m3x2 G 3 2 2 2 c c −2k2 + k 3 4 − 6 L L c c 4k 2 − k3 2 − 6 L L c −k 3 L c −k 3 4 − 6 L x1 c −k 3 2 − 6 x2 L x3 k3

m 0 0 m 0 0

mL2 IG = (long slender beam) 12 and V = 1 1 1 2 k1x2 + k2 x2 + k3 (x3 − x 4 ) 1 2 2 2 2

We need x1 , x 2 and x 4 in terms of x and θ . From geometry x = (x1 + x2 ) / 2 and θ = (x1 − x2 ) / L so that , x1 = x + and θL θL and x 2 = x − 2 2

L x 4 = x − θ c − 2 so that V = for the q i = x 1 θL 2 1 θL 2 1 L 2 k1 x + + k2 x − + k 3 x3 − x + θ c − 2 2 2 2 2 2

Contents

∂V θL θL L = k1 x + + k 2 x − − k3 x3 − x + θ c − ∂x 2 2 2 ∂T d ∂T and = m1x so that = m1x' ' ∂x dt 2 ∂x θL θL L ∴m 1x'' + k1 x + + k2 x − − k3 x3 − x + θ c − = 0 2 2 2 for the q i = θ ∂V L θL L θL L L = k1 x + − k 2 x − + k 3 c − x3 − x + θ c − ∂θ 2 2 2 2 2 2 ∂T d ∂T and = IG θ so that = IG θ'' ∂θ dt 2 ∂θ I Gθ' ' + k1 for the q i = x 3 L θL L θL L L x + − k 2 x − + k 3 c − x3 − x + θ c − = 0 2 2 2 2 2 2 ∂V L = k3 x3 − x + θ c − ∂x3 2 ∂T = m 3x3 so that d 2 ∂T = m 3x3'' and ∂x 3 dt ∂x3 L m 3x3' ' + k3 x3 − x + θ c − = 0 2 In matrix form m 1 0 0 IG 0 0 0 x'' 0 θ'' m 3 x 3' ' L (k1 − k2 )L / 2 − k 3 c − 2 L2 (k1 + k2 )L2 / 4 + k 3 c − 2 L k3 c − 2

..

k1 + k2 + k 3 L + (k1 − k 2 )L / 2 − k 3 c − 2 −k 3 This may be abbreviated to

x L k 3 c − θ = 0 2 x 3 k3 −k 3

[ M][ X]+ [K][ X] = 0 ........................................................ (5.13)

and in this case both the mass and stiffness matrices are symmetrical. Using the method described in chapter 4 we can find the eigen values and vectors and hence the natural frequencies and mode shapes.

Contents

5.1.1 Natural frequencies If we put λ = ω2 and [ X] = {u} we obtain the characteristic or eigen value equation of the system:

[[K ] − λ[M ]]{u} = 0

For non-zero solutions for {u}

.......................................................... (4.17)

det [[K ] − λ [M]] = 0 .......................................................... (4.18) Making the substitutions, k1 +k 2 +k 3 L det (k1 −k 2 )L/2 −k 3 c− 2 − k3

L 2 L 2 (k1 +k 2 )L2 / 4+k 3 c− 2 L k 3 c− 2 (k1 − k2 )L/ 2−k 3 c−

( )

( )

( ) ( )

m1 0 0 L k 3 c− − λ 0 I G 0 =0 2 0 0 m3 k3

−k3

( )

This is a cubic equation in λ so that there will be three natural frequencies. When each of these values is substituted in (4.17) the associated eigen vector {u} is found and this is the mode shape. For ease of comprehension it is convenient to represent the mode shapes in θL θL terms of the coordinates x1 = x + and x 2 = x − and x 3 as shown in figure 5.3. 2 2

Figure 5.3 Coordinate definition The amount of maths is significant so a program has been written to calculate the natural frequencies and the associated mode shapes. The results are shown in figure 5.4 for m 1 = 0.5kg; k1 = 10N / m; k 2 = 10N / m; m 3 = 0.2kg; k3 = 10N / m; c / L = 0.8

Contents

ω n1 = 0.4857 Hz

x1 0.2683 {u1} = x2 = 0.6008 x3 1.0

ω n2 = 0.9106 Hz

x1 1.0 {u2} = x2 = 0.0866 x3 −0.4227

It is helpful to see the modes animated.

It is also possible to investigate the effects of changing the parameter values.

Figure 5.4 Natural frequencies and mode shapes

ω n3 = 1.412 Hz

x1 −0.699 {u2} = x2 = 1.0 x3 −0.2247

5.1.2 Forced vibration Consider again the equations of motion mx1' ' = c c c c c 6b x1 −4k1 − k3 1− 4 − 6 + x2 2k 2 − k3 4 − 6 + x3k 3 4 − 6 + F 4 − L L L L L L ...... (5.9) mx 2' ' = c c c c c 6b x1 2k1 + k3 1− 2 − 6 + x2 −4k 2 + k3 2 − 6 − x3k 3 2 − 6 − F 2 − L L L L L L ..... (5.10) and c c m 3x3' ' = k 3 x1 1 − + x2 − x3 L L ........................................ (5.11)

Contents

For steady state vibration the exiting force is Fsinωt and x1=X1sinωt, x2=X2sinωt and x3=X3sinωt. This is because there is no damping. Substituting in (5.9), (5.10) and (5.11) c c c c c X1 4k1 + k3 1− 4 − 6 − mω2 + X2 −2k2 + k 3 4 − 6 + X3k 3 −4 + 6 L L L L L 6b = F4 − L c c c c c X1 −2k1 − k3 1− 2 − 6 + X2 4k 2 − k 3 2 − 6 − mω 2 + X3k 3 2 − 6 L L L L L 6b = −F 2 − L and c c − k 3X1 1− − k 3 X 2 + X3 k 3 − m 3ω 2 = 0 L L

(

)

and writing these in matrix form c c 4k1 + k3 1− 4 − 6 − mω2 L L c c −2k1 − k 3 1− 2 − 6 L L c − k3 1 − L c c −2k2 + k 3 4 − 6 L L c c 4k2 − k3 2 − 6 − mω2 L L c − k3 L 6b F 4 − L 6b = −F 2 − L 0 c k3 −4 + 6 L X1 c k 3 2 − 6 X2 L X3 k 3 − m3ω 2

....................5.14) X1 X2 X , and 3 F F F and hence the response X/F at any position Xpos along the beam may be determined after much maths it is possible to obtain the responses

Figure 5.5 Excitation and response positions.

Contents

Figure 5.6 Response of beam The response shown in figure 5.6 is for the same values as used previously with the excitation (Fpos) and response position (Xpos) as shown, Fpos = 0.3L; Xpos = 0.6L; Dpos = 0.8L; Note that there are infinite resonances at the natural frequencies because there is no damping. To examine the effect of varying the parameters run the program

Note that the resonant frequencies do not change with Fpos or Xpos. To get a better understanding of this system an animation program has been written.

This shows the system vibrating and allows the excitation frequency to be set by dragging the blue line on the graph (see figure 5.7). As shown this is set to a minimum on the graph. At this frequency the beam has no vibration amplitude for the Xpos value. The beam is vibrating as if there was a pivot at Xpos. It is of interest to investigate when the additional spring/mass acts as a detuner.

Contents

Figure 5.7 Excitation frequency variation. The detuned frequency is given by ω = k3 10 = = 4.472rad / s = 0.7118Hz . At this m3 0.5

frequency the beam will have no amplitude of vibration where the detuner is attached. This is because the detuner exhibits infinite dynamic stiffness at this frequency. If Xpos = Dpos then the response will be zero at this frequency. However note that the response graph is drawn at discrete frequencies and the frequency for zero response will be close to 0.7118 Hz but not exactly this value. If Xpos and Dpos are kept the same but varied it will be found that the response is always zero at 0.7118 Hz. 5.1.3 Effects of damping

Figure 5.6 Three degree of freedom system with damping. With damping as shown above we retain the detuner as it has no damping. The transient response involves a decaying vibration and is a major maths exercise to determine. For the steady state solution for sinusoidal excitation the previous undamped solution (equation

Contents

(5.14)) is easily modified. When viscous damping is included in parallel with the springs, wherever there was a kx term we will now have kx + cx' so that equation (5.14) becomes, c c c c c −2( k2 + iωc 2 )+k 3 4−6 k 3 −4+6 4(k1 +iωc1 )+k 3 1− L 4−6 L − mω 2 L L L X1 c c c c 2 k 2−6 c X −2(k1 + iωc1 )−k3 1− 2−6 4(k 2 +iωc 2 )−k3 2−6 − mω 3 L L L L L 2 c c X3 −k3 1− −k 3 k 3 −m 3ω2 L L 6b F 4− L 6b ............................................. (5.15) = −F 2− L 0

( )( ) ( )( ) ( )

( ) ( )

( ) ( )

( ) ( )

after even more maths it is possible to X X X obtain the responses 1 , 2 and 3 F F F (these will be complex, ie having amplitude and phase). The response X/F at any position Xpos along the beam may be determined. The response shown is for m 1 = 1.0kg; k1 = 10N / m; k2 = 10N / m; c1 = c 2 = 0.3Ns / m; Fpos = 0.3L; Xpos = 0.6L; Dpos = 0.8L; m 3 = 0.5kg; k3 = 10N / m;

Figure 5.7 Response of damped system.

It is of interest to note that the highest natural frequency is barely visible. The effect of damping may be examined using the program and varying c1 and c 2 .

There is still a detuner with a detuned frequency ω =

k3 10 = = 4.472rad / s = 0.7118Hz . m3 0.5

At this frequency the beam will have no amplitude of vibration where the detuner is attached. The animation program will allow this to be checked.

Contents

As for the undamped system, if Xpos = Dpos then the response will be zero at this frequency. If Xpos and Dpos are kept the same but varied it will be found that the response is always zero at 0.7118 Hz. 5.2 Five degrees-of-freedom 5.2.1 Natural frequencies To illustrate how the maths compounds and why computers are needed we will consider an axial system with five degrees-of-freedom as shown in figure 5.7.

Figure 5.7 System with five degrees of freedom. Using Newton's second law we obtain five equations of motion. m 1x1' ' = −k1x1 + k 2 (x2 − x1) m 2x 2' ' = −k 2 (x2 − x1) + k 3(x 3 − x2 ) m 3x3' ' = −k 3( x3 − x2 ) + k 4 (x 4 − x3 ) m 4x 4' ' = −k 4 ( x4 − x3 ) + k5( x5 − x 4 ) m 5x5' ' = −k 5( x5 − x 4 ) and if we rearrange and write these in matrix form we have, m1 0 0 0 0 x1' ' (k1 +k2 ) −k 2 0 0 0 x1 0 m x '' −k 0 0 0 2 (k 2 +k 3 ) −k 3 0 0 x2 2 2 0 0 m3 0 0 x3' ' + 0 −k 3 (k3 + k 4 ) −k 4 0 x 3 =0 .....(5.16) 0 0 0 m4 0 x4 '' 0 −k 4 (k 4 +k 5) −k 5 x 4 0 0 0 0 0 m 4 x5' ' 0 0 − k5 k 5 x 5 0 This may be abbreviated to the standard form

[ M][ X]+ [K][ X] = 0 ........................................................ (5.13)

As an example consider

..

Contents

The mass matrix [ M] is then

1.9 0 0 0 0 0 0 0 1.8 0 0 1.7 0 0 kg 0 0 0 0 1.6 0 0 0 0 0 1.5

25 −13 0 0 0 0 0 −13 27 −14 −14 29 −15 0 N/m and the stiffness matrix [ K] is then 0 0 0 −15 31 −16 0 0 0 −16 16 Using the eigen value approach and a computer program we find the five modes as shown in Figure 5.8. Mode 1 λ 1 = 0.0652 Mode 2 λ 2 = 0.5619 Mode 3 λ 3 = 1.392 Mode 4 λ 4 = 2.295 Mode 5 λ 5 = 3.212

rad 2 sec 2

rad 2 sec 2

rad 2 sec2

rad 2 sec 2 rad sec

rad 2 sec 2 rad sec

ω n1 = 0.2553

rad rad rad ω n2 = 0.7496 ω n3 = 1.18 sec sec sec 0.9075 1.0 {u}2 = 0.3635 −0.462 −0.9764

ω n4 = 1.515

ω n5 = 1.792

0.3309 0.6048 {u}1 = 0.8085 0.9389 1.0

1.0 −0.6988 −0.1603 −0.1116 1.0 0.4442 {u}3 = −0.944 {u}4 = −0.373 {u}5 = −0.8289 −0.2315 −0.6843 1.0 0.7588 0.5943 −0.4972

Figure 5.8 Natural frequencies and mode shapes of five degree of freedom system. The program used to calculate these values is available and the effect of varying parameter values may be investigated.

Contents

5.2.2 Forced vibration Consider forces applied to each of the masses. This is a common approach as then if we wish to have a force applied to only one mass we may set the others to zero.

Using Newton's second law we obtain five equations of motion. m 1x1' ' = −k1x1 + k 2 (x2 − x1) + F1 m 2x 2' ' = −k 2 (x2 − x1) + k 3(x 3 − x2 ) + F2 m 3x3' ' = −k 3( x3 − x2 ) + k 4 (x 4 − x3 ) + F3 m 4x 4' ' = −k 4 ( x4 − x3 ) + k5( x5 − x 4 ) + F4 m 5x5' ' = −k 5( x5 − x 4 ) + +F 5 and if we rearrange and write these in matrix form we have, m1 0 0 0 0 x1' ' (k1 +k2 ) −k 2 0 0 0 x1 F 1 0 m x '' −k 0 0 0 2 (k 2 +k 3 ) −k 3 0 0 x2 F 2 2 2 0 0 m3 0 0 x3' ' + 0 −k 3 (k3 + k 4 ) −k 4 0 x 3 = F3 0 0 0 m4 0 x4 '' 0 −k 4 (k 4 +k 5) −k 5 x 4 F 4 0 0 0 0 0 m 4 x5' ' 0 0 − k5 k 5 x 5 F5 0 As an example consider steady state sinusoidal excitation with F1(t) = F1 sin ωt and x1=X1sinωt, x2=X2sinωt, x3=X3sinωt, x4=X4sinωt and x5=X5sinωt. We then obtain, (k +k )−m ω 2 X F −k 2 0 0 0 1 1 2 1 1 −k2 (k2 +k3 )−m 2ω 2 − k3 0 0 X2 0 X = 0 2 0 −k 3 (k3 +k 4 )− m3ω − k4 0 3 2 0 0 −k 4 (k 4 + k5 )− m4ω − k5 X 4 0 0 0 0 −k 5 k 5 −m 5ω2 X5 0 ..........5.16) After much maths it is possible to obtain the responses X1 X2 X3 X 4 X , , , and 5 . F1 F1 F1 F1 F1

A computer program has been written to obtain the responses. For the example considered previously and for which the natural frequencies have been found.

Contents

The five responses are shown in figure 5.9.

Figure 5.9 Responses of each mass caused by a force on mass 1.

The program that calculates the response may be run and the parameter values changed and also which mass is excited. It will be found that the number of anti-resonances also changes with the excitation position. It should be noted that all responses have five resonances. These resonances are at the natural frequencies found previously and hence are at the same frequencies for each response. The 'height' of the resonances varies because the computer draws the graph at discrete frequencies. These discrete frequencies are not exactly at the natural frequencies (when the response would be infinite) and the responses are not the same off resonance. The responses have anti-resonances, when the response tends to zero. The 'depth' of these antiresonances again would be zero if the discrete frequency used by the program was at the X 'correct' value. It is of interest to note from figure 5.9 that for 1 there are 4 anti-resonances, F1 Contents

X2 X X there are 3 anti-resonances, for 3 there are 2 anti-resonances, for 4 there is 1 antiF1 F1 F1 X5 resonance, and for there are no anti-resonances, It is informative to view an animation of F1 for the motion at a particular excitation frequency. A program has been written to do this.

A typical output from the program is shown with an additional drawing of the system added at X the bottom. A typical response is shown - 4 . This illustrates that this program also allows F1 the steady state motion of the system to be observed at a particular frequency. This frequency may be changed by 'dragging' the blue line across the graph. In the diagram below the animation is frozen at an extreme of its motion. By comparison with the fixed drawing of the system the amplitudes can be deduced. As the frequency of the animation is set at 0.75 rad/s (the second natural frequency) the animation shows the second mode shape.

Figure 5.10 Output from animation program 'frozen' at extreme of amplitudes.

Contents

We may begin to understand the anti-resonance effect when the blue line is dragged to the appropriate frequency (1.035 rad/s). At this frequency the mass m 4 is not moving. However the mass m 5 to the right is moving with no excitation on it. Thus the mass m 5 and the stiffness k5 act as a detuner. The detuned frequency is k5 / m5 = 1.6 / 1.5 = 1.033rad / s . The discrepancy with 1.035 rad/s is because the program uses discrete frequencies.

5.2.3 Effects of damping

The transient response involves a decaying vibration and is a major maths exercise to determine. For the steady state solution for sinusoidal excitation of mass m 1 the previous undamped solution (equation (5.16)) is easily modified. When viscous damping is included in parallel with the springs, wherever there was a kx term we will now have kx + cx' so that equation (5.16) becomes, (k * +k * )−m ω2 −k* 0 0 0 1 2 1 2 X1 F1 * * + k* )−m ω 2 * X 2 0 −k 2 (k 2 3 −k 3 0 0 2 0 −k * (k* + k* )−m 3ω2 − k* 0 X3 = 0 3 3 4 4 * * + k* )−m ω2 * 0 0 −k 4 (k 4 5 −k 5 X 4 0 4 X 0 0 0 0 −k * k* −m 5ω2 5 5 5 where k* = k1 + iωc1 ; k* = k2 + iωc 2 ; k* = k3 + iωc3 ; k* = k 4 + iωc 4 and k* = k5 + iωc5 . 1 2 3 4 5 The maths is now even more prohibitive. A computer program similar to that for the undamped case allows the various responses to be checked.

A typical response is shown in figure 5.11. The effect of low levels of damping is to reduce the 'height' of the resonance peaks and to raise the anti-resonance 'depth'. If the animation is

Contents

observed for high damping it will be found that the masses do not just vibrate in-phase and out of phase, as is the case with undamped systems.

Figure 5.11 Response of damped five degree of freedom system. 5.4 Higher numbers of degrees of freedom The axial vibration program that has an animation at a selected frequency can also be used to increase and decrease the number of masses. The system may be made free/free or clamped/clamped. The damping may be either viscous or hysteretic. The instructions for running the program are available from it.

5.5 Conclusions As the number of degrees-of-freedom increases the maths to solve the equations increases significantly. In this chapter much of that maths has been omitted and the major trends have been highlighted. It should be clear that an 'N' degree of freedom system will have 'N' natural frequencies and associated mode shapes. When excited, resonance normally occurs at (undamped) or close to (damped) the natural frequencies. The response may have antiresonances when a detuner comes into play. _____________________________________________________ Next chapter

Contents

.

Contents