DOCSLIB.ORG

Course of Lectures in Mechanics Unit-F Dynamics of Rigid Bodies

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Course of Lectures in Mechanics Unit-F Dynamics of Rigid Bodies Course of Lectures in Mechanics Unit-F Dynamics of Rigid Bodies A. K. Kapoor http://0space.org/users/kapoor [email protected]; [email protected] Contents 1 Equations of Motion for a Rigid Body 3 1 Equations of motion in inertial frames . .... 3 § L 1.1 d equation–adifficulty ........................ 3 § dt 1.2 Usingbodyfixedaxes .......................... 4 § 2 Understandingprecession of spinningwheel . ...... 4 § 3 Torque free rotation does not always imply constant ~ω ............ 5 § 4 Principalaxesofarigidbody . 8 § 5 Eulerequations.................................. 10 § 6 EndNotes ..................................... 10 § 2 Examples :: Inertia Tensor and Finding Principal Axes 11 1 LessonOverview ................................. 11 § 1.1 LearningObjectives ........................... 11 § 1.2 Prerequisites ............................... 11 § 2 Computingmomentofinertiatensor . 11 § 2.1 Momentofinertiatensorofacube . 13 § 2.2 MI tensor of a cube relative to a corner . 14 § 2.3 Principal axes of a cube relative to a corner . 15 § 3 L~ is not always parallel to ~ω ........................... 15 § 4 EndNotes ..................................... 17 § 3 Applications of Euler’s Equations 18 1 LessonOverview ................................. 19 § 1.1 LearningGoals .............................. 19 § 1.2 Prerequisites ............................... 19 § 2 Bicyclewheelrevisited ............................. 20 § 3 FreeMotionofaRigidBody .......................... 21 § 3.1 Sphericaltop ............................... 21 § 4 Freesymmetricaltop............................... 22 § 4.1 Rigidrotator ............................... 22 § 1 4.2 Symmetrictop .............................. 22 § 5 Asymmetricaltop................................. 23 § 6 EndNotes ..................................... 24 § 6.1 Whytheearthwobbles ......................... 24 § 6.2 Kowalskeyatop.............................. 25 § Draft Copy for Limited Circulation 2 Lesson 1 Equations of Motion for a Rigid Body §1 Equations of motion in inertial frames The motion of rigid body is described by two sets of equations d2X dL cm = F , = ττ. (1.1) dt2 tot dt Important Recall that the separation of equations of motion of a rigid body takes place only if in the second equation above L is taken to be the angular momentum of the rigid body about the centre of mass of the body. Also, it should be remembered that the forces internal to the system, being of action reaction type do not contribute to the right hand sides of the above equations. The above equations are valid in inertial frames. For several purposes, the above equations are sufficient to get information about the motion of the rigid body. dL §1.1 dt equation – a difficulty Note that the second equation in (1.1) requires writing L w.r.t. the inertial frame. This involves moment of inertia tensor components computed relative to the inertial frame. Since the body is moving the components of the inertia tensor keep changing with time. Computing the moment of inertia tensor requires configuration of the body at time t, which in turn can be determined after equations of motion has been solved . / 3 §1.2 Using body fixed axes The above mentioned difficulty is solved by writing the the rate of change of angular momentum with respect to the frame fixed in the body. For this we use the result dL dL = + ω L (1.2) dt inf dt bf × The time dependence of angular momentum L = Iωin the body frame comes from the time dependence of angular velocity, the inertia tensor relative to the body frame remains constant in time. Thus we write the equation of motion as dL + ω L = ττ. (1.3) dt bf × All reference to the space axes has disappeared and we can drop the suffix bf, and write dL + ω L = ττ. (1.4) dt × In this form the above equation is valid if the body axes are chosen to have the origin at the centre of mass. It turns out that the above equations also hold when one point of the rigid body is fixed and that point is chosen as the origin. §2 Understanding precession of spinning wheel We discuss the dynamics of a spinning wheel. You have already watched this video on the moodle course site. A string is tied to one of the ends, P , of axle of a bicycle wheel and is the other end is supported manually by holding it with hand. The string is kept vertical by, for example tying the other end C of the string to the ceiling. When the end held by hand is released, the wheel falls due to gravity. However, if the wheel is set spinning before being released, it does not fall; it precesses in one direction. The direction of precession gets reversed when if the wheel is set spinning in the opposite direction. How do we understand this? Let the horizontal direction in which the axle points be chosen as X2 axis so that the angular velocity and angular momentum points along the X2. Let the X3 axis be chosen upwards so that the force F = mg due to gravity is along the negative X2 direction. The radius vector ~x of the centre of mass of the wheel is along the positive X axis. So that the torque, ~x F~ is then along the negative X axis. 2 × 1 The change in angular momentum in a short time time ∆t is given by ∆L~ = ~τ∆t and is shown in Fig.1(a). The vector sum of L~ and ∆L~ , drawn in Fig.1(b) is along P D. The implies that after time ∆t the axle has rotated about the X3 axis so that the end C moves to the point D. If the wheel is initially set spinning in the opposite direction, the torque and hence the change in angular momentum ∆L~ remains the same, but the direction of L~ is reversed 4 and the angular momentum after time ∆t will be as shown in Fig.1(c). This means that the end C must move to the point E and hence the axle rotates in the opposite direction. X 3 D X C 2 P E X 1 D L+ ∆L D K ∆ P L L P C C L ∆ L+ L E (a) (b) (c) Fig. 1 Movement of Spinning bicycle wheel [Q/T] := The above analysis has been presented from the point of view of the space axes. Analyze the motion from the point of view of the body frame. §3 Torque free rotation does not always imply constant ~ω )(Short Examples 1 What is the moment of inertia tensor for a rectangular lamina of || sides a and b. Choose the centre of the plate as the origin and X,Y axes parallel to the sides as shown in the figure. 5 Y O α X b a Fig. 2 Figure for the Examples 1 and 2 We divide the plate into small small pieces and find contribution of each piece to the moment of inertia tensor components. So at first divide the lamina into a large number of thin vertical strips parallel to the Y axis. consider one such vertical strip between x and x + dx. We will sum over all such strips by integrating over x from a/2 to a/2. − Next divide the strip at x into small small rectangular pieces by drawing horizontal lines parallel to the Y axis. One such rectangular element is shown as shaded. This element is located at (x,y, 0) has an area of dxdy and its mass will be σdxdy dm where ≡ σ is the mass per unit area. The total mass is M therefore σ = M/ab.The contribution 2 2 to 11 component moment of the inertia tensor is dm(y + z ) and the total value of I11 is obtained by summing (integrating) over all pieces. This is done by integrating over x from a/2 to a/2 over y from b/2 to b/2. So we have − − a/2 b/2 2 Ixx = σ(y + 0)dy dx (∵ z = 0) Z a/2 Z b/2 − n − o a/2 y3 b/2 = σ dx Z 3 a/2 n b/2 o − − a/2 b3 = σ dx Z a/2 12 − b3 a/2 = σ [x] dx 12 a/2 − b3 = σ a 12 Mb2 = (∵ σ = M/(ab)) (1.5) 12 Remember a double integral is computed by evaluated the inner integral over y first Then integrating the result over the outside variable x. In this example the order of integration is unimportant. 6 Ma2 Similarly I = . Also, for this case ( not in general) I = I + I . ( Prove this 22 12 33 xx yy yourself.) Next we compute I12. a/2 b/2 Ixy = σxydy dx − Z a/2 Z b/2 − n − o (1.6) The inner integral, over y, is b/2 y2 b/2 b2 b2 ydy = = = 0. (1.7) Z 2 8 − 8 b/2 b/2 − − Similarly other off diagonal components of the moment of inertia tensor turn out to be zero. This is because of our choice of the coordinate axes is symmetrical w.r.t. the plate. )(Short Examples 2 Find torque acting on a rectangular plate when it rotates with a || constant angular velocity about (a) the Z- axis (b) about its diagonal, so nˆ = 1 (1, 1, 0). √2 The choice of axes is as in the previous example. The moment of inertia tensor is given by b2 0 0 M I = 0 a2 0 (1.8) 12 0 0 (a2 + b2) dω Using dt = 0 in the Euler’s equations dL~ + ~ω L~ = ~τ (1.9) dt × total will give the required torque. Now the components of angular momentum are given by 2 2 Lx b 0 0 ωx b ωx M 2 M 2 Ly = 0 a 0 ωy = a ωy (1.10) 12 12 L 0 0 (a2 + b2) ω (a2 + b2)ω z z z ~ d~ω Note that the components of L are proportional to components of omega~ Hence dt = 0 dL~ means dt = 0. Therefore, The required torque is given by ~τ = ω L.~ (1.11) total × (a) In the first case angular velocity and angular momentum components are given by M ~ω = (0, 0,ω), L~ = (0, 0, (a2 + b2)) (1.12) 12 and required torque vanishes because ~ω and L~ parallel (both being along the Z- axis.) 7 (b) In this case Mω ω = ω(cos α, sin α, 0), L~ (b2 cos α, a2 sin α, 0) (1.13) 12 where α is the angle the diagonal makes with the X- axis.
Load more

Recommended publications
  • Mechanics 84 6.1 2D Particle in Cartesian Coordinates
    A Hostetler Handbook v 0.82 Contents Preface v 1 Newton's Laws1 1.1 Many Particles ................................. 2 1.2 Cartesian Coordinates............................. 3 1.3 Polar Coordinates ............................... 4 1.4 Cylindrical Coordinates ............................ 7 1.5 Air Resistance, Friction, and Buoyancy.................... 8 1.6 Charged Particle in a Magnetic Field..................... 17 1.7 Summary: Newton's Laws........................... 20 2 Momentum and Center of Mass 23 2.1 Rocket with no External Force ........................ 24 2.2 Multistage Rockets............................... 25 2.3 Rocket with External Force.......................... 25 2.4 Center of Mass................................. 26 2.5 Angular Momentum of a Single Particle................... 28 2.6 Angular Momentum of a System of Particles ................ 29 2.7 Angular Momentum of a Continuous Mass Distribution .......... 30 2.8 Summary: Momentum and Center of Mass ................. 34 3 Energy 36 3.1 General One-dimensional Systems ...................... 40 3.2 Atwood Machine................................ 43 3.3 Spherically Symmetric Central Forces .................... 44 3.4 The Energy of a System of Particles ..................... 45 3.5 Summary: Energy ............................... 47 4 Oscillations 49 4.1 Simple Harmonic Oscillators.......................... 49 4.2 Damped Harmonic Oscillators......................... 52 4.3 Driven Oscillators ............................... 57 4.4 Summary: Oscillations............................
    [Show full text]
  • Chapter 04 Rotational Motion
    Chapter 04 Rotational Motion P. J. Grandinetti Chem. 4300 P. J. Grandinetti Chapter 04: Rotational Motion Angular Momentum Angular momentum of particle with respect to origin, O, is given by l⃗ = ⃗r × p⃗ Rate of change of angular momentum is given z by cross product of ⃗r with applied force. p m dl⃗ dp⃗ = ⃗r × = ⃗r × F⃗ = ⃗휏 r dt dt O y Cross product is defined as applied torque, ⃗휏. x Unlike linear momentum, angular momentum depends on origin choice. P. J. Grandinetti Chapter 04: Rotational Motion Conservation of Angular Momentum Consider system of N Particles z m5 m 2 Rate of change of angular momentum is m3 ⃗ ∑N l⃗ ∑N ⃗ m1 dL d 훼 dp훼 = = ⃗r훼 × dt dt dt 훼=1 훼=1 y which becomes m4 x ⃗ ∑N dL ⃗ net = ⃗r훼 × F dt 훼 Total angular momentum is 훼=1 ∑N ∑N ⃗ ⃗ L = l훼 = ⃗r훼 × p⃗훼 훼=1 훼=1 P. J. Grandinetti Chapter 04: Rotational Motion Conservation of Angular Momentum ⃗ ∑N dL ⃗ net = ⃗r훼 × F dt 훼 훼=1 Taking an earlier expression for a system of particles from chapter 1 ∑N ⃗ net ⃗ ext ⃗ F훼 = F훼 + f훼훽 훽=1 훽≠훼 we obtain ⃗ ∑N ∑N ∑N dL ⃗ ext ⃗ = ⃗r훼 × F + ⃗r훼 × f훼훽 dt 훼 훼=1 훼=1 훽=1 훽≠훼 and then obtain 0 > ⃗ ∑N ∑N ∑N dL ⃗ ext ⃗ rd ⃗ ⃗ = ⃗r훼 × F + ⃗r훼 × f훼훽 double sum disappears from Newton’s 3 law (f = *f ) dt 훼 12 21 훼=1 훼=1 훽=1 훽≠훼 P.
    [Show full text]
  • Orthogonality
    Part VII Orthogonality 547 Section 31 Orthogonal Diagonalization Focus Questions By the end of this section, you should be able to give precise and thorough answers to the questions listed below. You may want to keep these questions in mind to focus your thoughts as you complete the section. What does it mean for a matrix to be orthogonally diagonalizable and why • is this concept important? What is a symmetric matrix and what important property related to diago- • nalization does a symmetric matrix have? What is the spectrum of a matrix? • Application: The Multivariable Second Derivative Test In single variable calculus, we learn that the second derivative can be used to classify a critical point of the type where the derivative of a function is 0 as a local maximum or minimum. Theorem 31.1 (The Second Derivative Test for Single-Variable Functions). If a is a critical number of a function f so that f 0(a) = 0 and if f 00(a) exists, then if f (a) < 0, then f(a) is a local maximum value of f, • 00 if f (a) > 0, then f(a) is a local minimum value of f, and • 00 if f (a) = 0, this test yields no information. • 00 In the two-variable case we have an analogous test, which is usually seen in a multivariable calculus course. Theorem 31.2 (The Second Derivative Test for Functions of Two Variables). Suppose (a; b) is a critical point of the function f for which fx(a; b) = 0 and fy(a; b) = 0.
    [Show full text]
  • Classical Mechanics and Symplectic Geometry
    CLASSICAL MECHANICS AND SYMPLECTIC GEOMETRY Maxim Jeffs Version: August 4, 2021 CONTENTS CLASSICAL MECHANICS 3 1 Newton’s Laws of Motion 3 2 Conservation Laws 7 3 Lagrangian Mechanics 10 4 Legendre Transform, Hamiltonian Mechanics 15 5 Problems 19 DIFFERENTIAL GEOMETRY 21 6 Constrained Mechanics, Smooth Manifolds 21 7 The Tangent Space 26 8 Differential Forms 30 9 Problems 35 SYMPLECTIC GEOMETRY 38 10 Symplectic Manifolds 38 11 Symplectic Mechanics 42 12 Lagrangian Submanifolds 47 13 Problems 51 SYMMETRIES IN MECHANICS 54 1 14 Lie Groups 54 15 Hamiltonian Group Actions 58 16 Marsden-Weinstein Theorem 64 17 Arnol’d-Liouville Theorem 70 18 The Hamilton-Jacobi Equation 74 19 Problems 75 FURTHER TOPICS 78 20 Pseudoholomorphic Curves 78 21 Symplectic Topology 82 22 Celestial Mechanics 87 23 Rigid Body Mechanics 88 24 Geometric Quantization 89 25 Kolmogorov-Arnol’d-Moser Theory 90 26 Symplectic Toric Manifolds 91 27 Contact Manifolds 92 28 Kähler Manifolds 93 Index 94 Bibliography 97 2 1 NEWTON’S LAWS OF MOTION Use the Force, Luke. Star Wars: A New Hope Firstly, what is classical mechanics? Classical mechanics is that part of physics that describes the motion of large-scale bodies (much larger than the Planck length) moving slowly (much slower than the speed of light). The paradigm example is the motion of the celestial bodies; this is not only the oldest preoccupation of science, but also has very important practical applications to navigation which continue to the present day. We shall have much to say about this example in this class. More than this, classical mechanics is a complete mathematical theory of the universe as it appears at such scales, with a rich mathematical structure; it has been studied by almost every great mathematician from Euler to Hilbert.
    [Show full text]
  • Space Exploration & Space Colonization Handbook
    Space Exploration & Space Colonization Handbook By USI, United Space Industries © Space Exploration & Space Colonization Handbook, 2018. This document lists the work that is involved in space exploration & space colonization. To find out more about a topic search ACOS, Australian Computer Operating System, scan the internet, go to a university library, go to a state library or look for some encyclopedia's & we recommend Encyclopedia Britannica. Eventually PAA, Pan Aryan Associations will be established for each field of space work listed below & these Pan Aryan Associations will research, develop, collaborate, innovate & network. 5 Worlds Elon Musk Could Colonize In The Solar System 10 Exoplanets That Could Host Alien Life 10 Major Players in the Private Sector Space Race 10 Radical Ideas To Colonize Our Solar System - Listverse 10 Space Myths We Need to Stop Believing | IFLScience 100 Year Starship 2001 Mars Odyssey 25 years in orbit: A celebration of the Hubble Space Telescope Recent Patents on Space Technology A Brief History of Time A Brown Dwarf Closer than Centauri A Critical History of Electric Propulsion: The First 50 Years A Few Inferences from the General Theory of Relativity. Einstein, Albert. 1920. Relativity: A Generation Ship - How big would it be? A More Efficient Spacecraft Engine | MIT Technology Review A New Thruster Pushes Against Virtual Particles! A Space Habitat Design A Superluminal Subway: The Krasnikov Tube A Survey of Nuclear Propulsion Technologies for Space Application A Survival Imperative for Space Colonization - The New York Times A `warp drive with more reasonable total energy requirements A new Mechanical Antigravity concept DeanSpaceDrive.Org A-type main-sequence star A.M.
    [Show full text]
  • Arxiv:1807.03867V2 [Physics.Class-Ph] 21 Feb 2019 Eetrve)Adohrtertcladcalculational Paper and Theoretical Other for Advantages
    Analytic formula for the Geometric Phase of an Asymmetric Top Nicholas A. Mecholsky∗ Vitreous State Laboratory The Catholic University of America Washington, DC 20064 (Dated: 1/21/2019) The motion of a handle spinning in space has an odd behavior. It seems to unexpectedly flip back and forth in a periodic manner as seen in a popular YouTube video.1 As an asymmetrical top, its motion is completely described by the Euler equations and the equations of motion have been known for more than a century. However, recent concepts of the geometric phase have allowed a new perspective on this classical problem. Here we explicitly use the equations of motion to find a closed form expression for total phase and hence the geometric phase of the force-free asymmetric top and explore some consequences of this formula with the particular example of the spinning handle for demonstration purposes. As one of the simplest dynamical systems, the asymmetric top should be a canonical example to explore the classical analog of Berry phase. PACS numbers: Keywords: geometric phase, Hannay angle, Berry phase, holonomy, Euler top, rigid body I. INTRODUCTION this is the asymmetric top in force-free motion in space. Some of the dynamic variables are periodic (like the an- gular momentum and angular velocity vectors). However The motion of spinning objects is ubiquitous. From a other variables (like the Euler angles) are not periodic. In seemingly mundane and classical motion, new myster- one cycle of the periodic variables, the system almost re- ies are still being debated and surprising features are turns to its original state except for the variables that are revealed.2–4 One important observation is the identifi- 5–7 not periodic.
    [Show full text]