Chapter 4 Rigid Body Dynamics
Total Page:16
File Type:pdf, Size:1020Kb
Chapter 4 Rigid Body Dynamics In Chapter 3 we developed the equations of motion for attitude kinematics. The main results of that chapter involve the description of attitude motion using attitude vari- ables, such as rotation matrices, Euler angles, Euler axis/angle sets, or quaternions. Knowing the values of any of these variables at a particular instant of time allows the attitude analyst to visualize the orientation of one reference frame with respect to another. Examples of mission-oriented analysis include determining the location of a target on the Earth, the direction to the Sun, and the local direction of the Earth’s magnetic field. We also developed in Chapter 3 the differential equations describing how these variables depend on the angular velocity of the reference frame. We did not investigate how the angular velocity, ~ω, is determined, and that is the subject of this chapter. In orbital dynamics, we typically model a spacecraft as a point mass or as a solid, rigid sphere, ignoring the rotational motion of the spacecraft about its mass center. This model is an excellent simplification for studying the motion of a satellite in its orbit — it works well for describing the motion of Europa about Jupiter, of Jupiter about the Sun, and of artificial satellites about the Earth. For artificial satellites, however, we usually need to know the orientation of the satellite with some degree of accuracy, so we need a more detailed model than the point mass model. A real spacecraft, of course, is a complicated piece of machinery with moving parts, flexible appendages, and partially filled fluid containers. However, even with all this complexity, a useful model is the rigid body, and that is the model we develop in this chapter. We begin with a rigorous development of the equations of motion for a rigid body, both translational and rotational. Because the torques applied to spacecraft are typically quite small, we find it useful to consider the special case of torque- free motion. As noted above, all spacecraft include non-rigid elements. The most notable effect of non-rigidity is the effect of energy dissipation, which has important implications for spacecraft design. We also present analysis of two important examples of rigid body dynamics: the Lagrange top, and the unbalanced rotor. 4-1 Copyright C. D. Hall September 20, 2005 4-2 CHAPTER 4. RIGID BODY DYNAMICS 4.1 Newton’s Second Law Newton’s laws are applicable to a particle, or point mass, m. As noted previously, Newton’s Second Law may be expressed as one second-order vector differential equa- tion m~¨r = ~f (4.1) or as two first-order vector differential equations ~p˙ = ~f (4.2) ~r˙ = ~p/m (4.3) Here ~r is the position vector of the particle with respect to an inertially fixed point, ~p is the linear momentum of the particle, and ~f is the total applied force acting on the particle. Either of these equations describe the dynamics of a particle. Integration of these equations requires selection of specific coordinates or position variables. For example, with respect to an inertial reference frame, we might represent ~r and ~f by T ~r = r1ˆi1 + r2ˆi2 + r3ˆi3 = r ˆi (4.4) T n o ~f = f1ˆi1 + f2ˆi2 + f3ˆi3 = f ˆi (4.5) In this case, the scalar second-ordern equationso of motion are mr¨1 = f1 (4.6) m¨r = f mr¨ = f (4.7) ⇔ 2 2 mr¨3 = f3 (4.8) and the scalar first-order equations of motion are p˙1 = f1 (4.9) p˙ = f p˙ = f (4.10) ⇔ 2 2 p˙3 = f3 (4.11) r˙1 = p1/m (4.12) r˙ = p/m r˙ = p /m (4.13) ⇔ 2 2 r˙3 = p3/m (4.14) One approach to developing the equations of motion of a spacecraft is to model the spacecraft as a system of point masses, mi, i = 1, ..., n and then develop the set of 3n second-order equations of motion or the set of 6n first-order equations of motion. This approach involves summations of the mass particles. For example, the n ~ n ~ total mass is m = i=1 mi, and the total applied force is f = i=1 fi. However, if we are to model a complicatedP spacecraft as a system of particles,P the number n will be large and the system of equations will be quite unwieldy. In the following section, we develop the rigid body as an approximate mathematical model for spacecraft attitude motion. Copyright C. D. Hall September 20, 2005 4.2. THE RIGID BODY MODEL 4-3 4.2 The Rigid Body Model Here we introduce the concept of a continuum of differential mass elements, dm.A differential mass element is quite a different concept from a mass particle. Whereas a mass particle has zero volume but finite mass, a differential mass element has infinitesimal volume and infinitesimal mass. The mathematical approach to collecting mass particles is summation, whereas the approach to collecting differential mass elements is integration. A common problem in studying rigid body motion is the problem of integrating over the volume of the body. This operation may have a scalar integrand, a vector integrand, or a tensor integrand, and all three of these types of integrands are impor- tant in rigid body motion. Integrating over the volume of the body is a triple integral operation, but we use the following shorthand notation: r1max r2max(r1) r3max(r1,r2) f(~r) dV = f(~r) dr3 dr2 dr1 (4.15) ZB Zr1min Zr2min(r1) Zr3min(r1,r2) Let us consider two examples, using a rectangular prism shape with a reference frame coinciding with one corner of the prism as shown in Fig. 4.1. The first example is the scalar integral corresponding to the total mass, and the second is the vector integral corresponding to the first moment of inertia. a c bˆ3 bˆ o 2 bˆ1 b Figure 4.1: A rectangular prism for computing mass and first moment of inertia Total mass. To compute the total mass of a rigid body, we need to know the density of the body, which we denote by µ. Units for µ are mass/length3; for example, aluminum has density µ = 2.85 103 kg/m3. In general, the density could vary from point to point within the body, or× it could vary as a function of time. Thus, generally Copyright C. D. Hall September 20, 2005 4-4 CHAPTER 4. RIGID BODY DYNAMICS µ = µ(~r,t). However, we generally assume that µ is constant, both spatially and temporally. Also, we usually use dm = µ dV as the variable of integration, even though m is not an independent variable. The meaning should always be clear. The total mass may then be expressed as m = µ dV = dm (4.16) ZB ZB Carrying out the indicated operations for the prism in Fig. 4.1, we obtain the expected result: a b c m = µ dr3 dr2 dr1 = µabc (4.17) Z0 Z0 Z0 Note that sometimes the total mass is called the zeroth moment of inertia, because it is a generalization of the first and second moments of inertia. First moment of inertia. This vector quantity is defined as ~co = ~ro dm (4.18) ZB Clearly this quantity depends on the origin from which ~r is measured, which motivates use of the superscript o to denote that ~r is measured from o, and the superscript o to denote that ~c is the first moment about o. The matrix version of this equation is co = ro dm (4.19) ZB To compute the first moment, we must choose a reference frame, and the compo- nents depend on that choice of frame. For this example, we use the body frame and corner origin of Fig. 4.1. Thus, ~c = c T bˆ , ~r = r T bˆ , and b { } b { } a b c o T µ 2 2 2 T cb = µ [r1 r2 r3] dr3 dr2 dr1 = a bc ab c abc (4.20) 0 0 0 2 Z Z Z h i Center of mass. We can use this result to determine the center of mass, c, which is defined as the point about which the first moment of inertia is zero; i.e. ~cc = ~0. Mathematically, this definition leads to ~cc = ~rc dm = ~0 (4.21) ZB The position vector from o to any point may be written as ~ro = ~roc +~rc (4.22) where ~roc is the position vector of c with respect to o, which is what we want to find, and ~rc is the position vector from c to a point in the body. The matrix version of this expression is ro = roc + rc (4.23) Copyright C. D. Hall September 20, 2005 4.2. THE RIGID BODY MODEL 4-5 where we omit the subscript b since the reference frame is clear. Substituting this expression into Eq. (4.19), we find co = (roc + rc) dm = roc dm = roc dm = mroc (4.24) ZB ZB ZB Thus, 1 µ roc = co = a2bc ab2c abc2 T (4.25) m 2m h i Using the expression for the mass m developed in Eq. (4.17), we obtain 1 roc = [abc] T (4.26) 2 which is obvious from the figure. However, this same procedure is applicable to more complicated geometric figures and to cases where µ is not constant. We now continue with the development of equations of motion for a rigid body, beginning with the linear momentum equations, and then developing the angular momentum equations. Linear momentum. The linear momentum of a rigid body is defined as ~p = ~v dm (4.27) ZB where ~v is the time derivative of the position vector ~r of a differential mass element dm = µ dV .