16.07 Dynamics Fall 2004 Lecture D18-2D Rigid Body Dynamics: Equations of Motion In this lecture, we will particularize the conservation principles presented in the previous lecture to the case in which the system of particles considered is a 2D rigid body. Mass moment of inertia In the previous lecture, we established that the angular momentum of a system of particles relative to the center of mass, G, was >(r G When considering a 2D rigid body, the velocity of any point relative to G consists of a pure rotation and therefore, can be expressed as U r here w is the angular velocity vector perpendicular to the plane of motion. These two equations can be combined to give Here, we have used the vector identity, A x(B C)=(A C)B-(A. B)C, and imposed the fact that r and w are perpendicular for 2d planar bodies If we define the mass moment of inertia of the body, IG, as IG iri
J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 Lecture D18 - 2D Rigid Body Dynamics: Equations of Motion In this lecture, we will particularize the conservation principles presented in the previous lecture to the case in which the system of particles considered is a 2D rigid body. Mass Moment of Inertia In the previous lecture, we established that the angular momentum of a system of particles relative to the center of mass, G, was HG = Xn i=1 (r ′ i × miv ′ i ) When considering a 2D rigid body, the velocity of any point relative to G consists of a pure rotation and, therefore, can be expressed as v ′ i = ω × r ′ i , where ω is the angular velocity vector perpendicular to the plane of motion. These two equations can be combined to give, HG = Xn i=1 (r ′ i × mi(ω × r ′ i )) = Xn i=1 mir ′ i 2ω . Here, we have used the vector identity, A × (B × C) = (A · C)B − (A · B)C, and imposed the fact that r ′ i and ω are perpendicular for 2D planar bodies. If we define the mass moment of inertia of the body, IG, as IG = Xn i=1 mir ′ i 2 , 1
HG=lc The moment of inertia, IG, is a scalar quantity. It is a property of the solid which indicates the way in which the mass of the solid is distributed relative to the center of mass. For example, if most of the mass is far away from the center of mass, ri will be large, resulting in a large moment of inertia. The dimensions of the moment of inertia are [M[L21 When the system of particles is a continuum, the summation over the number of particles is written as an / dm=M where p is the mass density and v is the volume. In this case, the position of the center of mass, rG, is with m=m dm, the total mass of the body When considering three dimensional bodies undergoing two dimensional motion, the moment of inertia needs to be defined with respect to an axis perpendicular to the plane of motion. In this case, we can still use equation 1, with r' replaced by the distance to the axis. It follows from the above definition that the moment of inertia of a composite body about a given point can always be calculated as the sum of the moments of inertia of the different components Exercise Show that the moment of inertia of a homogenous, slender bar of length L and mass M, about an axis perpen dicular to the bar and passing through the center of mass, is IG= ML/12, he moment of inertia of a homogeneous circular ring of mass M and radius R, about an axis perpen- dicular to the plane of the ring and passing through the center of mass, is IG=MR he moment of inertia of a homogeneous circular disc of mass M and radius R, about an axis perpen- dicular to the plane of the disc and passing through the center of mass, is IG=MR/ 2 Parallel Axis(Steiner)Theorem We will often need to find the moment of inertia with respect to a point other than the center of mass. For instance, the moment of inertia with respect to a given point, O, is defined /Pdm=M
then, HG = IGω . The moment of inertia, IG, is a scalar quantity. It is a property of the solid which indicates the way in which the mass of the solid is distributed relative to the center of mass. For example, if most of the mass is far away from the center of mass, r ′ i will be large, resulting in a large moment of inertia. The dimensions of the moment of inertia are [M][L 2 ]. When the system of particles is a continuum, the summation over the number of particles is written as an integral, IG = Z m r ′2 dm = Z V ρr′2 dV , (1) where ρ is the mass density and V is the volume. In this case, the position of the center of mass, rG, is given by rG = R m r dm m , with m = R m dm, the total mass of the body. When considering three dimensional bodies undergoing two dimensional motion, the moment of inertia needs to be defined with respect to an axis perpendicular to the plane of motion. In this case, we can still use equation 1, with r ′ replaced by the distance to the axis. It follows from the above definition that the moment of inertia of a composite body about a given point can always be calculated as the sum of the moments of inertia of the different components. Exercise Show that: - the moment of inertia of a homogenous, slender bar of length L and mass M, about an axis perpendicular to the bar and passing through the center of mass, is IG = ML2/12, - the moment of inertia of a homogeneous circular ring of mass M and radius R, about an axis perpendicular to the plane of the ring and passing through the center of mass, is IG = MR2 , - the moment of inertia of a homogeneous circular disc of mass M and radius R, about an axis perpendicular to the plane of the disc and passing through the center of mass, is IG = MR2/2. Parallel Axis (Steiner) Theorem We will often need to find the moment of inertia with respect to a point other than the center of mass. For instance, the moment of inertia with respect to a given point, O, is defined as IO = Z m r 2 dm = Z V ρr2 dV . 2
Assuming that O is a fixed point, Ho=low. If we know IG, then the moment of inertia with respect to another point, say point O, can be computed easily using the parallel axis theorem. given the relations r2=r.r and r=rG+r, we can then write r2dm=/(v+rar+r2)dm=n/dm+2re·/rdm+/r2dtm=m+, since m rdm=0. From this expression, it follows that the moment of inertia with respect to an arbitrary point is minimum when the point coincides with G. Hence, the minimum value for the moment of inertia is Radius of Gyration t is common to report the moment of inertia of a rigid body in terms of the radius of gyration, k. This is defined as and can be interpreted as the root-mean-square of the mass element distances from the axis of rotation Example Moment of inertia of a rectangular plate with circular cut-out We want to determine the moment of inertia of a rectangular plate of mass M kg about an axis perpendict to the plate and passing through point O ←2m,,1m m lm lm The density of the plate per unit area is Ap-A=0.3498Mkg/ where Ap= 6m2 is the area of the plate without the cut-out, and Ac=3.1416m2 is the area of the circle The moment of inertia of the plate about its center of mass is(Ip)G=pAp(32+22)/12. Using the parallel axis theorem, the moment of inertia of the plate about point O will be(Ip)o=(Ip)G+pAp (1.5+12) Similarly, for the circle, (Io)G=pAc12 and(Ic)o=(I)G+pAc(22+12). Finally, the moment of inertia of the plate with the circle cut-out about point O will b o=(l1)o-(1o=4334-6A=2.5012Mkgm2
Assuming that O is a fixed point, HO = IOω. If we know IG, then the moment of inertia with respect to another point, say point O, can be computed easily using the parallel axis theorem. Given the relations r 2 = r · r and r = rG + r ′ ., we can then write, IO = Z m r 2 dm = Z m (r 2 G + 2rG · r ′ + r ′2 ) dm = r 2 G Z m dm + 2rG · Z m r ′ dm + Z m r ′2 dm = mr2 G + IG , since R m r ′ dm = 0. From this expression, it follows that the moment of inertia with respect to an arbitrary point is minimum when the point coincides with G. Hence, the minimum value for the moment of inertia is IG. Radius of Gyration It is common to report the moment of inertia of a rigid body in terms of the radius of gyration, k. This is defined as k = r I m , and can be interpreted as the root-mean-square of the mass element distances from the axis of rotation. Example Moment of inertia of a rectangular plate with circular cut-out We want to determine the moment of inertia of a rectangular plate of mass M kg about an axis perpendicular to the plate and passing through point O. m m m m m The density of the plate per unit area is ρ = M Ap − Ac = 0.3498M kg/m2 , where Ap = 6m2 is the area of the plate without the cut-out, and Ac = 3.1416m2 is the area of the circle. The moment of inertia of the plate about its center of mass is (Ip)G = ρAp(32 + 22 )/12. Using the parallel axis theorem, the moment of inertia of the plate about point O will be (Ip)O = (Ip)G + ρAp(1.5 2 + 12 ). Similarly, for the circle, (Ic)G = ρAc1 2 and (Ic)O = (Ic)G + ρAc(22 + 12 ). Finally, the moment of inertia of the plate with the circle cut-out about point O will be, IO = (Ip)O − (Ic)O = 4.333ρAp − 6ρAc = 2.5012M kg m2 . 3
The radius of gyration is ko=√To/M=1.581m ee Reference [1(Appendix B)for a more examples on how to calculate moments of inertia. Equations of motion The equations describing the general motion of a rigid body follow from the conservation laws for systems of particles established in the last lecture. Since the general motion of a 2D rigid body can be determined by three parameters(e. g. I and y position of the center of mass, and rotation angle), we will need to supply three equations. Conservation of linear momentum yields one vector equation, or two scalar equations. The dditional condition is conservation of angular momentum We saw in the last lecture that there are several ways to express conservation of angular momentum. In principle, they are all equivalent, but, depending on the problem situation, the use of a particular form may greatly simplify the problem. The conservation of linear momentum for a system of particles yields the vector equation naG F where m is the body mass, aG is the acceleration of the center of mass, and F is the sum of the external forces acting on the bod Conservation of angular momentum about the center of mass, G,is simply IGo=MG here IG is the moment of inertia about the center of mass, a is the angular acceleration(recall that the ngular acceleration(and angular velocity ), is the same for all points in a rigid body), and MG is the total moment of the applied external forces(and moments)about the center of mass Although 3 is a vector equation, a and MG are always perpendicular to the plane of motion, and, therefore equation 3 only yields one scalar equation. It is interesting to note the similarity between equations 2 and 3 The moment of inertia, IG, can be interpreted as a measure of the bodys resistance to changing its angular relocity as a result of applied external moments If the applied forces and moments on the body are known, equations 2 and 3 can be integrated to determine he bodys trajectory. In most practical problems, however, one does not know a priori all the forces but a ombination of some forces and some characteristics of the motion In such cases where forces are unknown the additional equations must come from kinematic conditions(e. g. if the cylinder rolls without sliding, the acceleration of the center of mass is related to the angular velocity acceleration; also, the direction of the acceleration of the center of mass will be known in this case
The radius of gyration is kO = p IO/M = 1.5815m . See Reference [1] (Appendix B) for a more examples on how to calculate moments of inertia. Equations of Motion The equations describing the general motion of a rigid body follow from the conservation laws for systems of particles established in the last lecture. Since the general motion of a 2D rigid body can be determined by three parameters (e.g. x and y position of the center of mass, and rotation angle), we will need to supply three equations. Conservation of linear momentum yields one vector equation, or two scalar equations. The additional condition is conservation of angular momentum. We saw in the last lecture that there are several ways to express conservation of angular momentum. In principle, they are all equivalent, but, depending on the problem situation, the use of a particular form may greatly simplify the problem. The conservation of linear momentum for a system of particles yields the vector equation, maG = F , (2) where m is the body mass, aG is the acceleration of the center of mass, and F is the sum of the external forces acting on the body. Conservation of angular momentum about the center of mass, G, is simply IGα = MG , (3) where IG is the moment of inertia about the center of mass, α is the angular acceleration (recall that the angular acceleration (and angular velocity), is the same for all points in a rigid body), and MG is the total moment of the applied external forces (and moments) about the center of mass. Although 3 is a vector equation, α and MG are always perpendicular to the plane of motion, and, therefore, equation 3 only yields one scalar equation. It is interesting to note the similarity between equations 2 and 3. The moment of inertia, IG, can be interpreted as a measure of the body’s resistance to changing its angular velocity as a result of applied external moments. If the applied forces and moments on the body are known, equations 2 and 3 can be integrated to determine the body’s trajectory. In most practical problems, however, one does not know a priori all the forces but a combination of some forces and some characteristics of the motion. In such cases where forces are unknown, the additional equations must come from kinematic conditions (e.g. if the cylinder rolls without sliding, the acceleration of the center of mass is related to the angular velocity acceleration; also, the direction of the acceleration of the center of mass will be known in this case). 4
Example Cylinder on a ramp Let us consider a uniform cylinder of mass m and radius R rolling without slipping down a ramp of angle o y We consider the fixed reference frame, y, shown in the picture. The equations of motion, 2 and 3, are, in Wsinφ-F nyG N一Wcosφ IG In these equations, W= mg, but the normal force, N, and the friction force, F, are unknown. These two additional unknowns can be determined if we provide two additional kinematic conditions. First, we have that iG=0, from which we can determine N as Second, since the cylinder rolls without sliding, we have that iG=-Ro. Solving for iG, we obtain and F=(IGIG)/R2. For the uniform cylinder, we have that IG= mR2 /2 and iG=(2g sin)/3. If there was no friction, F would be zero, and the acceleration would be simply ig= gsin o. Also, if instead of having a uniform disc, we had a uniform ring with all the mass concentrated at the rim, then IG= mR2 (gsinφ)/2. Rotation about a fixed axis For cases in which there is a fixed point in the body, the motion of the body can be described with a single parameter(e. g. the rotation angle). In principle, we could still consider equations 2 and 3 and use the
Example Cylinder on a ramp Let us consider a uniform cylinder of mass m and radius R rolling without slipping down a ramp of angle φ. We consider the fixed reference frame, xy, shown in the picture. The equations of motion, 2 and 3, are, in this case, mx¨G = W sin φ − F my¨G = N − W cos φ IGα = −F R In these equations, W = mg, but the normal force, N, and the friction force, F, are unknown. These two additional unknowns can be determined if we provide two additional kinematic conditions. First, we have that ¨yG = 0, from which we can determine N as N = W cos φ . Second, since the cylinder rolls without sliding, we have that ¨xG = −Rα. Solving for ¨xG, we obtain x¨G = g sin φ 1 + (IG/mR2) , (4) and F = (IGx¨G)/R2 . For the uniform cylinder, we have that IG = mR2/2 and ¨xG = (2g sin φ)/3. If there was no friction, F would be zero, and the acceleration would be simply ¨xG = g sin φ. Also, if instead of having a uniform disc, we had a uniform ring with all the mass concentrated at the rim, then IG = mR2 and ¨xG = (g sin φ)/2. Rotation about a Fixed Axis For cases in which there is a fixed point in the body, the motion of the body can be described with a single parameter (e.g. the rotation angle). In principle, we could still consider equations 2 and 3 and use the 5