4569846498

01.05.2017 Views

Kinematics and dynamics of rigid bodies 71 2.13 Equations of motions If we consider the rigid body, Body 2, shown in Figure 2.36 we can formulate six equations of motion corresponding with the six degrees of freedom resulting from unconstrained motion. From our earlier consideration of linear momentum given in (2.161) we can write d d Σ{ F } { L } m { V G } dt d t 2 1 2 1 2 2 1 Expressing this in the familiar form of Newton’s second law we get (2.224) Σ{ F} m { A G } 2 1 2 2 1 (2.225) The vector equation given in (2.225) will thus yield the three equations associated with the translational motion of the body. It may be noted that for these equations the vectors in (2.225) may be conveniently referred to the fixed ground reference frame O 1 . In the same way that a resultant force acting on the body produces a change in linear momentum, a resultant moment will produce a change in angular momentum. If we consider the expression for angular momentum given in (2.178) we can obtain the equations of motion associated with rotational motion. For the rotational equations it is convenient to refer the vectors to the reference frame O 2 fixed in and rotating with Body 2: d d Σ{ M } { H } I dt d t [ ]{ } 2 G2 1/ 2 2 1/ 2 2 1/ 2 It can also be shown that d d { } d H2 { H2} { 2} { H2} t dt 11 / 12 / 12 / 12 / (2.226) (2.227) { M } / [ I2] / { 2} / { 2} / { H2} / G2 1 2 2 2 1 2 1 2 1 2 (2.228) Body 2 Z 2 {V G2 } 1 Z 1 X 2 O 2 G 2 Y 2 GRF {ω 2 } 1/2 X 1 O 1 Y 1 Fig. 2.36 Rigid body motion

72 Multibody Systems Approach to Vehicle Dynamics Expanding (2.228) gives ⎡M ⎢ ⎢ M ⎣⎢ M x y z ⎤ ⎡I I I ⎥ ⎢ ⎥ ⎢I I I ⎦⎥ ⎢ ⎣ I I I xx xy xz xy yy yz xz yz zz ⎤ ⎡ ⎥ ⎢ ⎥ ⎢ ⎥ ⎦ ⎣⎢ 2x 2y 2z ⎤ ⎡ 0 2z 2y ⎤ ⎡H ⎥ ⎢ ⎥ ⎥ ⎢ ⎢ 2z 0 2x⎥ ⎢ H ⎦⎥ ⎢ ⎥ ⎣ 2y 2x 0 ⎦ ⎣⎢ H (2.229) Substituting in terms for H 2x , H 2y and H 2z now leads to the equations given in (2.230) to (2.231). For convenience we can drop the subscript for Body 2: M I I ( ) I ( ) x xx x xy y x z xz z x y 2 2 ( Izz Iyy ) yz Iyz( y z ) 2x 2y 2z ⎤ ⎥ ⎥ ⎦⎥ (2.230) M I ( ) I I ( ) y xy x y z yy y yz z x y 2 2 xx zz x z xz z x ( I I ) I ( ) (2.231) M I ( ) I ( ) z xz x y z yz y x z 2 2 zz z yy xx x y xy x y I ( I I ) I ( ) (2.232) In summary the rotational equations of motion for Body 2 may be written in vector form as Σ{ M (2.233) G2} 12 / [ I2] 22 / { 2} 12 / [ 2] 12 / [ I2] 22 / { 2} 12 / Hence we can see that in setting up the equations of motion for any rigid body the translational equations for all bodies in a system may conveniently be referred to a single fixed inertial frame O 1 . The rotational equations, however, are better referred to a body centred frame, in this case O 2 . A considerable simplification in these equations will result if frame O 2 is selected such that its axes are the principal axes of the body (I 1 I xx , I 2 I yy , I 3 I zz ) and the products of inertia are zero. The equations that result are known as Euler’s equations of motion: M x I 1 x (I 3 I 2 ) y z (2.234) M y I 2 y (I 1 I 3 ) x z (2.235) M z I 3 z (I 2 I 1 ) x y (2.236) The equations given in (2.234) to (2.236) become even simpler when the motion of a body is constrained so that rotation takes place in one plane only. If, for example, rotation about the x- and y-axes are prevented then equation (2.236) reduces to the more familiar form associated with twodimensional motion: M z I 3 z (2.237) The following example also demonstrates how gyroscopic effects associated with three-dimensional motion may be identified. If we consider the swing arm suspension system shown in Figure 2.37 we can take the suspension arm Body 2 to be constrained by a revolute joint to rotate with a constant angular velocity of 10 rad/s about the axis of the joint as shown. While this motion is in progress the road wheel, Body 3, is also rotating

72 Multibody Systems Approach to Vehicle Dynamics<br />

Expanding (2.228) gives<br />

⎡M<br />

⎢<br />

⎢<br />

M<br />

⎣⎢<br />

M<br />

x<br />

y<br />

z<br />

⎤ ⎡I I I<br />

⎥ ⎢<br />

⎥<br />

⎢I I I<br />

⎦⎥<br />

⎢<br />

⎣<br />

I I I<br />

xx xy xz<br />

xy yy yz<br />

xz yz zz<br />

⎤ ⎡<br />

⎥ ⎢<br />

⎥ ⎢<br />

<br />

⎥<br />

⎦ ⎣⎢<br />

<br />

2x<br />

2y<br />

2z<br />

⎤ ⎡ 0 2z<br />

2y<br />

⎤ ⎡H<br />

⎥ ⎢<br />

⎥<br />

⎥<br />

⎢ <br />

⎢<br />

2z<br />

0 2x⎥<br />

⎢<br />

H<br />

⎦⎥<br />

⎢<br />

⎥<br />

⎣<br />

2y<br />

2x<br />

0<br />

⎦ ⎣⎢<br />

H<br />

(2.229)<br />

Substituting in terms for H 2x , H 2y and H 2z now leads to the equations given<br />

in (2.230) to (2.231). For convenience we can drop the subscript for Body 2:<br />

M I I ( ) I<br />

( )<br />

x xx x xy y x z xz z x y<br />

2 2<br />

(<br />

Izz Iyy ) yz Iyz( y z<br />

)<br />

2x<br />

2y<br />

2z<br />

⎤<br />

⎥<br />

⎥<br />

⎦⎥<br />

(2.230)<br />

M I ( ) I I<br />

( )<br />

y xy x y z yy y yz z x y<br />

2 2<br />

xx zz x z xz z x<br />

(<br />

I I ) I<br />

( )<br />

(2.231)<br />

M I ( ) I<br />

( )<br />

z xz x y z yz y x z<br />

2 2<br />

zz z yy xx x y xy x y<br />

I (<br />

I I ) I<br />

( )<br />

(2.232)<br />

In summary the rotational equations of motion for Body 2 may be written<br />

in vector form as<br />

Σ{ M (2.233)<br />

G2} 12 / [ I2] 22 / { 2} 12 / [ 2] 12 / [ I2] 22 / { 2}<br />

12 /<br />

Hence we can see that in setting up the equations of motion for any rigid<br />

body the translational equations for all bodies in a system may conveniently<br />

be referred to a single fixed inertial frame O 1 . The rotational equations,<br />

however, are better referred to a body centred frame, in this case O 2 .<br />

A considerable simplification in these equations will result if frame O 2 is<br />

selected such that its axes are the principal axes of the body (I 1 I xx ,<br />

I 2 I yy , I 3 I zz ) and the products of inertia are zero. The equations that<br />

result are known as Euler’s equations of motion:<br />

M x I 1 x (I 3 I 2 ) y z (2.234)<br />

M y I 2 y (I 1 I 3 ) x z (2.235)<br />

M z I 3 z (I 2 I 1 ) x y (2.236)<br />

The equations given in (2.234) to (2.236) become even simpler when the<br />

motion of a body is constrained so that rotation takes place in one plane<br />

only. If, for example, rotation about the x- and y-axes are prevented then<br />

equation (2.236) reduces to the more familiar form associated with twodimensional<br />

motion:<br />

M z I 3 z (2.237)<br />

The following example also demonstrates how gyroscopic effects associated<br />

with three-dimensional motion may be identified. If we consider the<br />

swing arm suspension system shown in Figure 2.37 we can take the suspension<br />

arm Body 2 to be constrained by a revolute joint to rotate with a<br />

constant angular velocity of 10 rad/s about the axis of the joint as shown.<br />

While this motion is in progress the road wheel, Body 3, is also rotating

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

Saved successfully!

Ooh no, something went wrong!