This preview shows page 1. Sign up to view the full content.
Unformatted text preview: J. Peraire, S. Widnall 16.07 Dynamics Fall 2008 Version 2.0 Lecture L30  3D Rigid Body Dynamics: Tops and Gyroscopes 3D Rigid Body Dynamics: Euler Equations in Euler Angles
In lecture 29, we introduced the Euler angles as a framework for formulating and solving the equations for conservation of angular momentum. We applied this framework to the freebody motion of a symmetrical body whose angular momentum vector was not aligned with a principal axis. The angular moment was however constant. We now apply Euler angles and Euler's equations to a slightly more general case, a top or gyroscope in the presence of gravity. We consider a top rotating about a fixed point O on a flat plane in the presence of gravity. Unlike our previous example of freebody motion, the angular momentum vector is not aligned with the Z axis, but precesses about the Z axis due to the applied moment. Whether we take the origin at the center of mass G or the fixed point O, the applied moment about the x axis is Mx = M gzG sin, where zG is the distance to the center of mass.. Initially, we shall not assume steady motion, but will develop Euler's equations in the Euler angle variables (spin), (precession) and (nutation). 1 Referring to the figure showing the Euler angles, and referring to our study of freebody motion, we have the following relationships between the angular velocities along the x, y, z axes and the time rate of change of the Euler angles. The angular velocity vectors for , and are shown in the figure. Note that these three angular velocity vectors are not orthogonal, giving rise to some cross products when the angular velocities i are calculated about the three principal axes. x = sin sin + cos y = sin cos  sin z = cos + (1) (2) (3) 3D Rigid Body Dynamics: Euler's Equations
We consider a symmetric body, appropriate for a top, for which the moments of inertia Ixx = Iyy = I0 and Izz = I. The angular momentum is then Hx = I0 x Hy = Hz = I0 y Iz (4) (5) (6) For the general motion of a threedimensional body, we have Euler's equations in bodyfixed axes which rotate with the body so that the moment of inertia is constant in time. In this bodyfixed coordinate system, the conservation of angular momentum is d H= ([I]{}) = AppliedM oments dt 2 (7) Since we have chosen to work in a rotating coordinate system so that d dt I = 0, we must pay the price, applying Coriolis theorem to obtain the time derivative of the angular velocity vector in the rotating coordinate system d H = H + H, dt (8) resulting in the "Euler" equations expressed in the x, y, z coordinate system moving with the body. In general, we must rotate with the total angular velocity of the body, so that the governing equation for the conservation of angular momentum become, with = . Mx My Mz = Hx  Hy z + Hz y (9) (10) (11) = Hy  Hz x + Hx z = Hz  Hx y + Hy x where is the rotational angular velocity of our axis system. In this case because of the symmetry of the body, we are free to choose to allow the spin velocity to rotate relative to our bodyfixed axes system and fasten our axis system to and . That is, we follow the motion of the top in and but allow it to spin with angular velocity relative to our coordinate system rotating with angular velocity . So that x = y = sin (12) (13) (14) z = cos. For this choice of coordinate system, we have essentially performed a rotation in and only, leading to the geometry shown. 3 In this coordinate system, since the rotation did not occur, the angular velocity of the body is x = y = sin (15) (16) (17) z = cos + . For this choice of coordinate system, we have Hx = Hy = I0 I0 sin (18) (19) (20) Hz = I( cos + ) resulting in Hx = Hy = I0 I0 (sin + cos (21) (22) (23) Hz = I( cos  sin + ) and H in components is (I  I0 ) cos sin 2 + I sin (I0  I) cos  I 0. (24) (25) (26) Note that the z component is zero. Since the axes chosen are principal axes, the final form of Euler's equations becomes, Mx My Mz = I0 (  2 sin cos ) + I sin ( cos + ) = I0 ( sin + 2 cos )  I( cos + ) = I( + cos  sin ) (27) (28) (29) For a top, Mx = M gzG sin, My = 0 and Mz = 0. These equations are unsteady and nonlinear. We can gain insight by examining the character of some special solutions and constants of the motion. Steady Precession: Gyroscopic Motion
We now consider the steady precession of a top about the Z axis. In terms of the variables we have defined, the top rotates with spin velocity about its principal axis, and precesses with angular velocity while 4 maintaining a constant angle with the vertical axis Z. We then have = constant = 0 = constant = 0 =0 =0 ==0 =0 (30) (31) (32) For a steady precession of a top, Euler's equations reduce to sin(I(cos + )  I0 cos) = Mx = M gsinzG Note that sin cancels, resulting in I  (I0  I)2 cos = M gzG (34) (33) In the usual case for tops or gyroscopes, we have >> so that 2 may be ignored. Therefore, for steady precession, the relationship between the precession angular velocity and the spin angular velocity is = M gzG /(I ) (35) where I = Izz , the moment of inertia of the gyroscope about its spin axis. This result is also true if = 0. For this motion, the angular momentum vector is not aligned with the Z axis as for freebody motion, but is in the plane of z, Z, and rotates around the Z axis according to the applied external moment which is constant and in the x direction. In the limit as >> the angular momentum vector is essentially along the axis of rotation, z, with a slight component due to which can be calculated after has been determined. Some attempt to sketch this is shown in the figure, but even this small correction is an approximation for >> . For >> , this same result can be obtained by less complex approach, but it is important to realize that this is an approximation. For the spin velocity much greater than the precessional velocity, we may take 5 the angular momentum vector as directed along the spin or z axis, H = Ik. Then the applied moment direction), this give a result for in agreement with the previous result = M gzG /(I) will precess the vector with rotation rate H. Since = K = K, (where K is a unit vector in the Z (36) Since dropped out of the previous equation, for = 0, we again have = M gzG /(I). The limit >> is the "gyroscopic" limit where the device behaves as a gyroscope rather than as the more general case of a top. The difference is that, for a gyroscope, is larger than any other rotation rate in the system, such as the angular velocity of an aircraft or spacecraft. This makes the gyroscope a useful basis for many instruments. We shall return to this issue. Unsteady Precession of a Top: Integrals of Motion
For the general case of a top in a gravitational field, we have Euler's Equations. In the most general case, we will have spin , precession and nutation , all varying with time. Mx My Mz = I0 (  2 sin cos ) + I sin (cos + ) = M gzG sin = I0 ( sin + 2 cos )  I( cos + ) = 0 = I( + cos  sin ) = 0 (37) (38) (39) The consequences of My = 0 and Mz = 0 provide two quantities that must be constant/conserved during the motion. (We will later show that these results can also be obtained by applying Lagrange's equation to this system.) d/dt(I0 sin2 + I cos ( + cos)) = d/dt(I( + cos )) = dp dt dp dt =0 =0 (40) (41) It can easily be seen by inspection that differentiating p and setting it equal to 0 yields the equation Mz = 0. The equation from
dp dt is more complex; to obtain My requires combinations of 6 dp dt and dp dt . But the conclusion is that p and p are constants of the motion; we will show later that these can be considered generalized momenta, derived by application of Lagrange's equation. Therefore, the spin angular velocity and the precession angular velocity are instantaneously related to the nutation angle through p  p cos = I0 sin2 p p  p cos  = Isin2 I0 We may therefore write the remaining governing equation as +(
p I0 (42) (43)  p p I cos)( I sin3  p cos ) I0  gM zG =0 I0 (44) This reduces the problem to the motion of a single variable ; once (t) has been determined, the constancy of p /I0 and p /I during the motion, give the relations between the nutation angle (t), the spin angle (t), and precession angle (t). For this more complex motion, all three angles change with time, and the tip of the top traces out a motion, inscribed on the surface of a sphere for visualization, as shown in the figure. These various motions are referred to as unidirectional precession, looping precession and cuspidal motion. For more discussion of these solutions see J.B Marion and S.T. Thornton, Classical Dynamics, Chapter 11. Spinning Top by Lagrange's Equation
The constancy of two momenta obtained by application of Euler's equation can be found perhaps more directly by application of Lagrange's equation. We write the kinetic energy of a spinning top as
2 2 2 T = 1/2I0 (x + y ) + 1/2Iz (45) We use the full form of the angular velocities in the Euler system moving with the spinning top and obtain
2 x = ( sin sin + cos )2 2 y = (sin cos  sin )2 2 z = (cos + )2 (46) (47) (48) 7 resulting in
2 2 x + y = 2 sin2 + 2 2 z = ( cos + )2 . (49) (50) Therefore the kinetic energy is T = 1/2I0 (2 sin2 + 2 ) + 1/2I( cos + )2 For a top, the potential energy is V = M gzG cos so that the Lagrangian is L = 1/2I0 (2 sin2 + 2 ) + 1/2I( cos + )2  M gzG cos . (53) (52) (51) We consider , and to be our generalized coordinates, qi . And we move forward applying Lagrange's that equation, d/dt(L/qi )  V /qi = 0. However, since V /qi = 0 for both qi = and qi = , we conclude L = (I0 sin2 + I cos2 ) + I cos = constant = p L = I( + cos ) = constant = p (54) (55) We define these groupings which must be constant in time as generalized momenta, consistent with our earlier exposition of Lagrange's equations. These are of course the same constants identified using Euler's equations in equations (3839). Lagrange's equations win! The final governing equation is obtained from d L L  =0 (56) dt resulting in +( in agreement with equation (44). We also observe that total energy is a constant of the motion. E = 1/2 I0 (2 sin2 + theta2 ) + 1/2I0 cos + )2 + M HzG cos (58)
p I0  p p I cos)( I sin3  p cos ) I0  gM zG = 0, I0 (57) For more discussion of these solutions see J.B Marion and S.T. Thornton, Classical Dynamics, Chapter 11. Extreme Tops: Gyroscopes
Considerable simplification and practical application results if the angular velocity is an order of magnitude larger than angular velocity associated with the motions of a system of interest, such as an aircraft or a 8 spacecraft. In this case, the gyroscope can be used as an instrument to measure some quantity of interest in a vehicle. In order to use a gyroscope as an instrument, we first have to consider how it is mounted and affixed to the vehicle. The classic mount is called a gimbal. It permits free rotation of the gyroscope about its own axis, and depending on how complex the gimbal system is, free rotation about other axis as well. Consider the gyroscope shown in the figure. It is free to spin about its axis with angular velocity . It is attached to a gimbal which is free to rotate about the support as sketched. The entire apparatus is mounted on a turntable rotating with angular velocity , where >> . If we can measure the orientation of the axis of rotation of the disk relative to the turntable mount, we could possibly use this as an instrument. What does it measure? Consider that the initial angular momentum vector of the disk is in the horizontal plane. The axis of the disk is free to rotate about the horizontal axis through an angle . As the turntable rotates with angular velocity , the axis of the disk will rotate in its only free direction to keep its rate of change of angular momentum zero, since it can experience no moments/torques. There are two source of changes in angular momentum, the unsteady component if changes with time; and the component of the angular momentum of the disk in the horizontal plane that is precessed by the turntable rotation . The unsteady term is I0 , where I0 is the moment of inertia of the disk about the axis perpendicular to its spin axis; the vector is in the horizontal direction. The horizontal component of the angular momentum of the disk is Isin, where I is the moment of inertia of the disk about its axis of rotation. The rate of 9 change is also in the horizonal direction, perpendicular to and of magnitude Isin. The sum of these two terms, which must be zero, is I0 + Isin = 0 (59) This is the familiar equation for the oscillation of a pendulum. If there is a small amount of friction in the support, the system will settle at = 0 and will indicate the direction of rotation of the turntable axis. This is perhaps not too interesting, since we already know the direction of rotation of the turntable. But if we have this device on a rotating body, such as the earth, it will again indicate the direction of the axis of rotation of the body, in this case the earth. We call this direction "North". Thus, this application of the gyroscope could be used as a compass, often called a gyrocompass. ADDITIONAL READING J.L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition 7/9 J.B Marion and S.T. Thornton, Classical Dynamics, Chapter 11 W.T. Thompson, Introduction to Space Dynamics, Chapter 5 J. H. Ginsberg, Advanced Engineering Dynamics, Second Edition, Chapter 8 D. Kleppner, R.J. Kolenkow, An introduction to Mechanics,Chapter 7 10 MIT OpenCourseWare http://ocw.mit.edu 16.07 Dynamics
Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ...
View
Full
Document
This note was uploaded on 05/02/2011 for the course DYNAM 101 taught by Professor Matuka during the Spring '11 term at MIT.
 Spring '11
 MAtuka

Click to edit the document details