ROTATION: a Review of Useful Theorems Involving Proper Orthogonal Matrices Referenced to Three- Dimensional Physical Space

Total Page:16

File Type:pdf, Size:1020Kb

ROTATION: a Review of Useful Theorems Involving Proper Orthogonal Matrices Referenced to Three- Dimensional Physical Space Unlimited Release Printed May 9, 2002 ROTATION: A review of useful theorems involving proper orthogonal matrices referenced to three- dimensional physical space. Rebecca M. Brannon† and coauthors to be determined †Computational Physics and Mechanics T Sandia National Laboratories Albuquerque, NM 87185-0820 Abstract Useful and/or little-known theorems involving33× proper orthogonal matrices are reviewed. Orthogonal matrices appear in the transformation of tensor compo- nents from one orthogonal basis to another. The distinction between an orthogonal direction cosine matrix and a rotation operation is discussed. Among the theorems and techniques presented are (1) various ways to characterize a rotation including proper orthogonal tensors, dyadics, Euler angles, axis/angle representations, series expansions, and quaternions; (2) the Euler-Rodrigues formula for converting axis and angle to a rotation tensor; (3) the distinction between rotations and reflections, along with implications for “handedness” of coordinate systems; (4) non-commu- tivity of sequential rotations, (5) eigenvalues and eigenvectors of a rotation; (6) the polar decomposition theorem for expressing a general deformation as a se- quence of shape and volume changes in combination with pure rotations; (7) mix- ing rotations in Eulerian hydrocodes or interpolating rotations in discrete field approximations; (8) Rates of rotation and the difference between spin and vortici- ty, (9) Random rotations for simulating crystal distributions; (10) The principle of material frame indifference (PMFI); and (11) a tensor-analysis presentation of classical rigid body mechanics, including direct notation expressions for momen- tum and energy and the extremely compact direct notation formulation of Euler’s equations (i.e., Newton’s law for rigid bodies). Computer source code is provided for several rotation-related algorithms. A draft version of this document is available at http://www.me.unm.edu/~rmbrann/gobag.htmlDRAF Intentionally Left Blank DRAFTii Introduction ............................................................................................................ 1 Orthogonal basis & coordinate transformations.......................................... 4 Motivational discussion: Principal basis ...................................................... 4 Orthonormal basis transformations............................................................... 5 Alternative direction cosine matrix................................................................ 10 Coordinate transformations........................................................................... 10 Rotation operations.............................................................................................. 12 Example: rotation of an angle α about the X-axis ........................................ 14 Example: rotation of an angle α about the Y-axis. ....................................... 15 Example: rotation of an angle α about the Z-axis. ....................................... 16 Where does that negative sign go? ................................................................ 17 Faster way to write down the matrix of a rotation tensor ............................. 17 Specific example: A 90˚ rotation about the Z-axis......................................... 18 Axis and angle of rotation.................................................................................. 19 T Euler-Rodrigues formula ..................................................................................... 20 Computing the rotation tensor given axis and angle............................................ 22 Example.......................................................................................................... 23 Another example ............................................................................................ 24 Another similar example................................................................................ 24 Numerical example (for testing computer codes) .......................................... 26 Alternative way to construct the rotation tensor............................................ 27 Some properties of the axial tensor ............................................................... 28 Argyris’s form of the Euler-Rodrigues formula............................................. 30 Corollary to the Euler-Rodrigues formula: Existence of a preferred basis ........................................................................ 31 Computing axis and angle given the rotation tensor............................................ 32 Finding the principal rotation angle ............................................................. 32 Finding the principal rotation axis................................................................ 33 Method 1 algorithm for axis and angle of rotation ....................................... 36 Example.......................................................................................................... 37 Another example ............................................................................................ 38 Numerical Example (for testing codes).......................................................... 39 Method 2 algorithm for computing axis and angle........................................ 40 Rotations contrasted with reflections .............................................................. 41 Quaternion representation of a rotation ........................................................ 43 Shoemake’s form [3]...................................................................................... 43 A more structural direct form ........................................................................ 43 Relationship between quaternion and axis/angle forms....................................... 44 Dyad form of an invertible linear operator.................................................... 45 SPECIAL CASE: lab basis............................................................................. 45 SPECIAL CASE: rotation .............................................................................. 46 Sequential Rotations ............................................................................................ 47 Sequential rotations about fixed (laboratory) axes. ............................................. 47 EULER ANGLES: Sequential rotations about “follower” axes.......................... 49 Converting Euler angles to direction cosines................................................ 49 Example: ....................................................................................................... 49 DRAFiii Converting Euler angles to axis and angle.................................................... 50 Series expression for a rotation......................................................................... 51 Spectrum of a rotation......................................................................................... 53 Sanity check ................................................................................................... 54 Polar decomposition............................................................................................ 55 Difficult definition of the deformation gradient .................................................. 55 Intuitive definition of the deformation gradient................................................... 59 Converse problem: interpreting a deformation gradient matrix ................... 61 The Jacobian of the deformation.......................................................................... 62 Invertibility of a deformation............................................................................... 63 Sequential deformations....................................................................................... 63 Matrix analysis version of the polar decomposition theorem.............................. 64 The polar decomposition theorem — a hindsight intuitive introduction............. 65 A more rigorous (classical) presentation of the polar decomposition theorem ... 68 THEOREM:.................................................................................................... 69 PROOF: ......................................................................................................... 69 Working with the inverse gradient................................................................. 72 The *FAST* way to do a polar decomposition in two dimensions..................... 73 “Mixing” or interpolating rotations................................................................... 75 proposal #1: Map and re-compute the polar decomposition................................ 75 proposal #2: Discard the “stretch” part of a mixed rotation. ............................... 76 Proof: ............................................................................................................. 77 proposal #3: mix the pseudo-rotation-vectors...................................................... 78 proposal #4: mix the quaternions......................................................................... 78 Rates of rotation .................................................................................................... 79 The “spin” tensor ................................................................................................. 79 The angular velocity vector ................................................................................. 80 Angular velocity in terms of axis and angle of rotation......................................
Recommended publications
  • Rotational Motion (The Dynamics of a Rigid Body)
    University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Robert Katz Publications Research Papers in Physics and Astronomy 1-1958 Physics, Chapter 11: Rotational Motion (The Dynamics of a Rigid Body) Henry Semat City College of New York Robert Katz University of Nebraska-Lincoln, [email protected] Follow this and additional works at: https://digitalcommons.unl.edu/physicskatz Part of the Physics Commons Semat, Henry and Katz, Robert, "Physics, Chapter 11: Rotational Motion (The Dynamics of a Rigid Body)" (1958). Robert Katz Publications. 141. https://digitalcommons.unl.edu/physicskatz/141 This Article is brought to you for free and open access by the Research Papers in Physics and Astronomy at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Robert Katz Publications by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. 11 Rotational Motion (The Dynamics of a Rigid Body) 11-1 Motion about a Fixed Axis The motion of the flywheel of an engine and of a pulley on its axle are examples of an important type of motion of a rigid body, that of the motion of rotation about a fixed axis. Consider the motion of a uniform disk rotat­ ing about a fixed axis passing through its center of gravity C perpendicular to the face of the disk, as shown in Figure 11-1. The motion of this disk may be de­ scribed in terms of the motions of each of its individual particles, but a better way to describe the motion is in terms of the angle through which the disk rotates.
    [Show full text]
  • The Rotation Group in Plate Tectonics and the Representation of Uncertainties of Plate Reconstructions
    Geophys. J. In?. (1990) 101, 649-661 The rotation group in plate tectonics and the representation of uncertainties of plate reconstructions T. Chang’ J. Stock2 and P. Molnar3 ‘Department of Mathematics, University of Virginia, Charlottesville, Virginia 22903-3199, USA ‘Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138, USA ’Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA Accepted 1989 December 27. Received 1989 December 27; in original form 1989 August 11 Downloaded from SUMMARY The calculation of the uncertainty in an estimated rotation requires a parametriza- tion of the rotation group; that is, a unique mapping of the rotation group to a point in 3-D Euclidean space, R3. Numerous parametrizations of a rotation exist, http://gji.oxfordjournals.org/ including: (1) the latitude and longitude of the axis of rotation and the angle of rotation; (2) a representation as a Cartesian vector with length equal to the rotation angle and direction parallel to the rotation axis; (3) Euler angles; or (4) unit length quaternions (or, equivalently, Cayley-Klein parameters). The uncertainty in a rotation is determined by the effect of nearby rotations on the rotated data. The uncertainty in a rotation is small, if rotations close to the best fitting rotation degrade the fit of the data by a large amount, and it is large, if only rotations differing by a large amount cause such a degradation. Ideally, we would at California Institute of Technology on August 29, 2014 like to parametrize the rotations in such a way so that their representation as points in R3 would have the property that the distance between two points in R3 reflects the effects of the corresponding rotations on the fit of the data.
    [Show full text]
  • PHYSICS of ARTIFICIAL GRAVITY Angie Bukley1, William Paloski,2 and Gilles Clément1,3
    Chapter 2 PHYSICS OF ARTIFICIAL GRAVITY Angie Bukley1, William Paloski,2 and Gilles Clément1,3 1 Ohio University, Athens, Ohio, USA 2 NASA Johnson Space Center, Houston, Texas, USA 3 Centre National de la Recherche Scientifique, Toulouse, France This chapter discusses potential technologies for achieving artificial gravity in a space vehicle. We begin with a series of definitions and a general description of the rotational dynamics behind the forces ultimately exerted on the human body during centrifugation, such as gravity level, gravity gradient, and Coriolis force. Human factors considerations and comfort limits associated with a rotating environment are then discussed. Finally, engineering options for designing space vehicles with artificial gravity are presented. Figure 2-01. Artist's concept of one of NASA early (1962) concepts for a manned space station with artificial gravity: a self- inflating 22-m-diameter rotating hexagon. Photo courtesy of NASA. 1 ARTIFICIAL GRAVITY: WHAT IS IT? 1.1 Definition Artificial gravity is defined in this book as the simulation of gravitational forces aboard a space vehicle in free fall (in orbit) or in transit to another planet. Throughout this book, the term artificial gravity is reserved for a spinning spacecraft or a centrifuge within the spacecraft such that a gravity-like force results. One should understand that artificial gravity is not gravity at all. Rather, it is an inertial force that is indistinguishable from normal gravity experience on Earth in terms of its action on any mass. A centrifugal force proportional to the mass that is being accelerated centripetally in a rotating device is experienced rather than a gravitational pull.
    [Show full text]
  • Chapter 5 ANGULAR MOMENTUM and ROTATIONS
    Chapter 5 ANGULAR MOMENTUM AND ROTATIONS In classical mechanics the total angular momentum L~ of an isolated system about any …xed point is conserved. The existence of a conserved vector L~ associated with such a system is itself a consequence of the fact that the associated Hamiltonian (or Lagrangian) is invariant under rotations, i.e., if the coordinates and momenta of the entire system are rotated “rigidly” about some point, the energy of the system is unchanged and, more importantly, is the same function of the dynamical variables as it was before the rotation. Such a circumstance would not apply, e.g., to a system lying in an externally imposed gravitational …eld pointing in some speci…c direction. Thus, the invariance of an isolated system under rotations ultimately arises from the fact that, in the absence of external …elds of this sort, space is isotropic; it behaves the same way in all directions. Not surprisingly, therefore, in quantum mechanics the individual Cartesian com- ponents Li of the total angular momentum operator L~ of an isolated system are also constants of the motion. The di¤erent components of L~ are not, however, compatible quantum observables. Indeed, as we will see the operators representing the components of angular momentum along di¤erent directions do not generally commute with one an- other. Thus, the vector operator L~ is not, strictly speaking, an observable, since it does not have a complete basis of eigenstates (which would have to be simultaneous eigenstates of all of its non-commuting components). This lack of commutivity often seems, at …rst encounter, as somewhat of a nuisance but, in fact, it intimately re‡ects the underlying structure of the three dimensional space in which we are immersed, and has its source in the fact that rotations in three dimensions about di¤erent axes do not commute with one another.
    [Show full text]
  • The Spin, the Nutation and the Precession of the Earth's Axis Revisited
    The spin, the nutation and the precession of the Earth’s axis revisited The spin, the nutation and the precession of the Earth’s axis revisited: A (numerical) mechanics perspective W. H. M¨uller [email protected] Abstract Mechanical models describing the motion of the Earth’s axis, i.e., its spin, nutation and its precession, have been presented for more than 400 years. Newton himself treated the problem of the precession of the Earth, a.k.a. the precession of the equinoxes, in Liber III, Propositio XXXIX of his Principia [1]. He decomposes the duration of the full precession into a part due to the Sun and another part due to the Moon, to predict a total duration of 26,918 years. This agrees fairly well with the experimentally observed value. However, Newton does not really provide a concise rational derivation of his result. This task was left to Chandrasekhar in Chapter 26 of his annotations to Newton’s book [2] starting from Euler’s equations for the gyroscope and calculating the torques due to the Sun and to the Moon on a tilted spheroidal Earth. These differential equations can be solved approximately in an analytic fashion, yielding Newton’s result. However, they can also be treated numerically by using a Runge-Kutta approach allowing for a study of their general non-linear behavior. This paper will show how and explore the intricacies of the numerical solution. When solving the Euler equations for the aforementioned case numerically it shows that besides the precessional movement of the Earth’s axis there is also a nu- tation present.
    [Show full text]
  • 4-2 Degrees and Radians
    4-2 Degrees and Radians Write each degree measure in radians as a multiple of π and each radian measure in degrees. 10. 30° SOLUTION: To convert a degree measure to radians, multiply by ANSWER: 14. SOLUTION: To convert a radian measure to degrees, multiply by ANSWER: 120° Identify all angles that are coterminal with the given angle. Then find and draw one positive and one negative angle coterminal with the given angle. 20. 225° SOLUTION: All angles measuring are coterminal with a angle. Sample answer: Let n = 1 and −1. eSolutions Manual - Powered by Cognero Page 1 ANSWER: 225° + 360n°; Sample answer: 585°, −135° 24. SOLUTION: All angles measuring are coterminal with a angle. Sample answer: Let n = 1 and −1. ANSWER: Find the length of the intercepted arc with the given central angle measure in a circle with the given radius. Round to the nearest tenth. 29. , r = 4 yd SOLUTION: ANSWER: 5.2 yd 30. 105°, r = 18.2 cm SOLUTION: Method 1 Convert 105° to radian measure, and then use s = rθ to find the arc length. Substitute r = 18.2 and θ = . Method 2 Use s = to find the arc length. ANSWER: 33.4 cm Find the rotation in revolutions per minute given the angular speed and the radius given the linear speed and the rate of rotation. 36. = 104π rad/min SOLUTION: The angular speed is 104π radians per minute. Each revolution measures 2π radians 104π ÷ 2π = 52 The angle of rotation is 52 revolutions per minute. ANSWER: 52 rev/min 37. v = 82.3 m/s, 131 rev/min SOLUTION: 131 × 2π = 262π The linear speed is 82.3 meters per second with an angle of rotation of 262π radians per minute.
    [Show full text]
  • Hydraulics Manual Glossary G - 3
    Glossary G - 1 GLOSSARY OF HIGHWAY-RELATED DRAINAGE TERMS (Reprinted from the 1999 edition of the American Association of State Highway and Transportation Officials Model Drainage Manual) G.1 Introduction This Glossary is divided into three parts: · Introduction, · Glossary, and · References. It is not intended that all the terms in this Glossary be rigorously accurate or complete. Realistically, this is impossible. Depending on the circumstance, a particular term may have several meanings; this can never change. The primary purpose of this Glossary is to define the terms found in the Highway Drainage Guidelines and Model Drainage Manual in a manner that makes them easier to interpret and understand. A lesser purpose is to provide a compendium of terms that will be useful for both the novice as well as the more experienced hydraulics engineer. This Glossary may also help those who are unfamiliar with highway drainage design to become more understanding and appreciative of this complex science as well as facilitate communication between the highway hydraulics engineer and others. Where readily available, the source of a definition has been referenced. For clarity or format purposes, cited definitions may have some additional verbiage contained in double brackets [ ]. Conversely, three “dots” (...) are used to indicate where some parts of a cited definition were eliminated. Also, as might be expected, different sources were found to use different hyphenation and terminology practices for the same words. Insignificant changes in this regard were made to some cited references and elsewhere to gain uniformity for the terms contained in this Glossary: as an example, “groundwater” vice “ground-water” or “ground water,” and “cross section area” vice “cross-sectional area.” Cited definitions were taken primarily from two sources: W.B.
    [Show full text]
  • Rotation Matrix - Wikipedia, the Free Encyclopedia Page 1 of 22
    Rotation matrix - Wikipedia, the free encyclopedia Page 1 of 22 Rotation matrix From Wikipedia, the free encyclopedia In linear algebra, a rotation matrix is a matrix that is used to perform a rotation in Euclidean space. For example the matrix rotates points in the xy -Cartesian plane counterclockwise through an angle θ about the origin of the Cartesian coordinate system. To perform the rotation, the position of each point must be represented by a column vector v, containing the coordinates of the point. A rotated vector is obtained by using the matrix multiplication Rv (see below for details). In two and three dimensions, rotation matrices are among the simplest algebraic descriptions of rotations, and are used extensively for computations in geometry, physics, and computer graphics. Though most applications involve rotations in two or three dimensions, rotation matrices can be defined for n-dimensional space. Rotation matrices are always square, with real entries. Algebraically, a rotation matrix in n-dimensions is a n × n special orthogonal matrix, i.e. an orthogonal matrix whose determinant is 1: . The set of all rotation matrices forms a group, known as the rotation group or the special orthogonal group. It is a subset of the orthogonal group, which includes reflections and consists of all orthogonal matrices with determinant 1 or -1, and of the special linear group, which includes all volume-preserving transformations and consists of matrices with determinant 1. Contents 1 Rotations in two dimensions 1.1 Non-standard orientation
    [Show full text]
  • Exploring Artificial Gravity 59
    Exploring Artificial gravity 59 Most science fiction stories require some form of artificial gravity to keep spaceship passengers operating in a normal earth-like environment. As it turns out, weightlessness is a very bad condition for astronauts to work in on long-term flights. It causes bones to lose about 1% of their mass every month. A 30 year old traveler to Mars will come back with the bones of a 60 year old! The only known way to create artificial gravity it to supply a force on an astronaut that produces the same acceleration as on the surface of earth: 9.8 meters/sec2 or 32 feet/sec2. This can be done with bungee chords, body restraints or by spinning the spacecraft fast enough to create enough centrifugal acceleration. Centrifugal acceleration is what you feel when your car ‘takes a curve’ and you are shoved sideways into the car door, or what you feel on a roller coaster as it travels a sharp curve in the tracks. Mathematically we can calculate centrifugal acceleration using the formula: V2 A = ------- R Gravitron ride at an amusement park where V is in meters/sec, R is the radius of the turn in meters, and A is the acceleration in meters/sec2. Let’s see how this works for some common situations! Problem 1 - The Gravitron is a popular amusement park ride. The radius of the wall from the center is 7 meters, and at its maximum speed, it rotates at 24 rotations per minute. What is the speed of rotation in meters/sec, and what is the acceleration that you feel? Problem 2 - On a journey to Mars, one design is to have a section of the spacecraft rotate to simulate gravity.
    [Show full text]
  • Rotational Motion CONTENTS CHAPTER-OPENING QUESTION—Guess Now! 8–1 Angular Quantities a Solid Ball and a Solid Cylinder Roll Down a Ramp
    You too can experience rapid rotation—if your stomach can take the high angular velocity and centripetal acceleration of some of the faster amusement park rides. If not, try the slower merry-go-round or Ferris wheel. Rotating carnival rides have rotational kinetic energy as well as angular momentum. Angular acceleration is produced by a net torque, and rotating objects have rotational kinetic energy. P T A E H R C 8 Rotational Motion CONTENTS CHAPTER-OPENING QUESTION—Guess now! 8–1 Angular Quantities A solid ball and a solid cylinder roll down a ramp. They both start from rest at the 8–2 Constant Angular Acceleration same time and place. Which gets to the bottom first? 8–3 Rolling Motion (a) They get there at the same time. (Without Slipping) (b) They get there at almost exactly the same time except for frictional differences. 8–4 Torque (c) The ball gets there first. 8–5 Rotational Dynamics; (d) The cylinder gets there first. Torque and Rotational Inertia (e) Can’t tell without knowing the mass and radius of each. 8–6 Solving Problems in Rotational Dynamics ntil now, we have been concerned mainly with translational motion. We 8–7 Rotational Kinetic Energy discussed the kinematics and dynamics of translational motion (the role 8–8 Angular Momentum and of force). We also discussed the energy and momentum for translational Its Conservation U *8–9 Vector Nature of motion. In this Chapter we will deal with rotational motion. We will discuss the Angular Quantities kinematics of rotational motion and then its dynamics (involving torque), as well as rotational kinetic energy and angular momentum (the rotational analog of linear momentum).
    [Show full text]
  • NASA Tit It-462 &I -- .7, F
    N ,! S’ A TECHNICAL, --NASA Tit it-462 &I REP0R.T .7,f ; .A. SPACE-TIME. TENsOR FQRMULATION _ ,~: ., : _. Fq. CONTINUUM MECHANICS IN ’ -. :- -. GENERAL. CURVILINEAR, : MOVING, : ;-_. \. AND DEFQRMING COORDIN@E. SYSTEMS _. Langley Research Center Ehnpton, - Va. 23665 NATIONAL AERONAUTICSAND SPACE ADhiINISTRATI0.N l WASHINGTON, D. C. DECEMBER1976 _. TECHLIBRARY KAFB, NM Illllll lllllIII lllll luII Ill Ill 1111 006857~ ~~_~- - ;. Report No. 2. Government Accession No. 3. Recipient’s Catalog No. NASA TR R-462 5. Report Date A SPACE-TIME TENSOR FORMULATION FOR CONTINUUM December 1976 MECHANICS IN GENERAL CURVILINEAR, MOVING, AND 6. Performing Organization Code DEFORMING COORDINATE SYSTEMS 7. Author(s) 9. Performing Orgamzation Report No. Lee M. Avis L-10552 10. Work Unit No. 9. Performing Organization Name and Address 1’76-30-31-01 NASA Langley Research Center 11. Contract or Grant No. Hampton, VA 23665 13. Type of Report and Period Covered 2. Sponsoring Agency Name and Address Technical Report National Aeronautics and Space Administration 14. Sponsoring Agency Code Washington, DC 20546 %. ~u~~l&~tary Notes 6. Abstract Tensor methods are used to express the continuum equations of motion in general curvilinear, moving, and deforming coordinate systems. The space-time tensor formula- tion is applicable to situations in which, for example, the boundaries move and deform. Placing a coordinate surface on such a boundary simplifies the boundary-condition treat- ment. The space-time tensor formulation is also applicable to coordinate systems with coordinate surfaces defined as surfaces of constant pressure, density, temperature, or any other scalar continuum field function. The vanishing of the function gradient com- ponents along the coordinate surfaces may simplify the set of governing equations.
    [Show full text]
  • Chapter 16 Two Dimensional Rotational Kinematics
    Chapter 16 Two Dimensional Rotational Kinematics 16.1 Introduction ........................................................................................................... 1 16.2 Fixed Axis Rotation: Rotational Kinematics ...................................................... 1 16.2.1 Fixed Axis Rotation ........................................................................................ 1 16.2.2 Angular Velocity and Angular Acceleration ............................................... 2 16.2.3 Sign Convention: Angular Velocity and Angular Acceleration ................ 4 16.2.4 Tangential Velocity and Tangential Acceleration ....................................... 4 Example 16.1 Turntable ........................................................................................... 5 16.3 Rotational Kinetic Energy and Moment of Inertia ............................................ 6 16.3.1 Rotational Kinetic Energy and Moment of Inertia ..................................... 6 16.3.2 Moment of Inertia of a Rod of Uniform Mass Density ............................... 7 Example 16.2 Moment of Inertia of a Uniform Disc .............................................. 8 16.3.3 Parallel Axis Theorem ................................................................................. 10 16.3.4 Parallel Axis Theorem Applied to a Uniform Rod ................................... 11 Example 16.3 Rotational Kinetic Energy of Disk ................................................ 11 16.4 Conservation of Energy for Fixed Axis Rotation ............................................
    [Show full text]