Rotations of Moment of Inertia Tensor
|
|
- Augusta Daniel
- 5 years ago
- Views:
Transcription
1 Rotations of Moment of Inertia Tensor using Quaternions Mikica B Kocic, , v0.3 There is a wonderful connection between complex numbers (and quaternions, as their extension) and geometry in which, translations correspond to additions, rotations and scaling to multiplications, and reflections to conjugations. However, the beauty of describing spatial rotations with quaternions is shadowed by an expression like: - 2 Iy 2 + z 2 M 2 xy- 2 wz 2 wy+ 2 xz 2 xy+ 2 wz - 2 Ix 2 + z 2 M 2 yz- 2 wx ÿ 2 xz- 2 wy 2 wx+ 2 yz - 2 Ix 2 + y 2 M - 2 Iy 2 + z 2 M 2 xy- 2 wz 2 wy+ 2 xz T ÿ ÿ 2 xy+ 2 wz - 2 Ix 2 + z 2 M 2 yz- 2 wx 2 xz- 2 wy 2 wx+ 2 yz - 2 Ix 2 + y 2 M which actually represents such simple transformation as the rotation of moment of inertia tensor using quaternions [2]. The purpose of this document is to express spatial rotations of a moment of inertia tensor in the pure quaternionic form and to highlight the algebraical meaning behind such rotations by proving some interesting theorems. ü Conventions and Definitions Load quaternions package (see Quat.nb for the annotated source code) Get@ "Quat.m", Path > 8 < D; ô Conventions In this document, matrices, vectors and matrix representation of quaternions are shown in uppercase font (roman-bold font in text) and quaternions and other numbers are shown in lowercase font (italic font in text). Since the symbol I in Mathematica is used for imaginary unit, the selected symbol for moment of inertia is the uppercase script J i.e.. The identity matrix is suffixed with a number of dimensions, i.e. I4 represents the identity matrix in 4. Matrices 4ä4 originated from quaternions may be geometrically partitioned into scalar (real), vector (pure imaginary) and cross-product parts [3]. Such partitions are visualized using row and column divider lines throughout this document: scalar É vector É É É É É vector É cross product É É É É É ô General symbols and variables used in theorems (definitions) Hamilton's quaternions q, q and q 2 with real number components (but not with complex number components as biquaternions; see definition of default $Assumptions bellow): q w, x, y, z D; q w, x, y, z D; q2 w2, x2, y2, z2 D; Identity matrix I4 in 4 : I4 = IdentityMatrix@4D;
2 2 Rotations of Tensors using Quaternions v0.3.nb Common 4 matrices A and B: A = Aww Awx Awy Awz Axw Axx Axy Axz Ayw Ayx Ayy Ayz Azw Azx Azy Azz ; B = Bww Bwx Bwy Bwz Bxw Bxx Bxy Bxz Byw Byx Byy Byz Bzw Bzx Bzy Bzz ; Symmetrical 4 matrix representing moment of inertia tensor : = Ixx Ixy Ixz ; 0 Ixy Iyy Iyz 0 Ixz Iyz Izz Default assumption used for all asserted statements and expressions in this notebook is that previously defined entities q, q, q 2, A, B and are on field of, i.e. that their components are elements in : $Assumptions = latten@8 ToList ê@ 8 q, q, q2 <, A,B, < Reals; ü Different Types of Multiplications The ordirary notation for multiplication in Mathematica is A. B a ä b p ** q multiplication between matrices ("dot" multiplication) cross-product between 3D vectors quaternion multiplication (both in Quat.m and Mathematica Quaternions package) Beside the ordinary notation for multiplication, this document also introduces the following notation based on the center dot ( ) operator. A B p q q A := multiplication between matrices := quaternion multiplication := column-wise quaternion-matrix multiplication Note that the center dot operator in Mathematica is without built-in meaning. ô Matrix multiplication Define center dot operator A B as matrix multiplication (otherwise, Mathematica uses plain dot "." for matrix multiplications). CenterDot@ m? MatrixQ D := Dot@ m D ô Quaternion multiplication Mathematica uses non-commutative multiply (**) operator as symbol for quaternion multiplication. The Quat.m package and this document follow this convention. Here we define also the center dot operator as symbol for quaternion multiplication, just for convinience and esthetics CenterDot@ q? Q D := NonCommutativeMultiply@ q D
3 Rotations of Tensors using Quaternions v0.3.nb 3 ô Column-wise (or horizontal) quaternion-matrix multiplication Let us now define multiplication between quaternion q and matrix M, where it is assumed that matrix is a row vector of quaternions in columns, i.e. M = H q q 2... q n L, so that quaternion multiplication is done between q and q j in every column j. q q j q n qq q j = qq j qq n q q 2 q n q q 2 q n q ÿ q 2 q 22 q 2n q 3 q q 3n := q 2 q 22 q 2n q 3 q q 3n q 4 q 42 q 4n q 4 q 42 q 4n Since the quaternion multiplication is not commutative, there are several cases to distinguish: a) Left multiplication q ÿ M, where a quaternion in each column in M is multiplied by the quaternion q from the left: CenterDot@ q_? Q, mat_? MatrixQ D := Transpose@ TableA ToListA q mat PAll,jT, 9j,, Dimensions@matD P2T = b) Right multiplication M ÿ q, where a quaternion in each column in M is multiplied by the quaternion q from the right: CenterDot@ mat_? MatrixQ, q_? Q D := Transpose@ TableA ToListA mat PAll,jT q, 9j,, Dimensions@matD P2T = c) Compound multiplication from both sides p ÿ M ÿ q, where multiplication of each column in M by the quaternion p from the left and q from the right: CenterDot@ q_? Q, mat_? MatrixQ, q2_? Q D := Transpose@ TableA ToListA q mat PAll,jT q2, 9j,, Dimensions@matD P2T = The multiplication between a quaternion and a matrix is associative, i.e. it holds: q ÿ M ÿ q 2 ã q ÿ HM ÿ q 2 L ã Hq ÿ ML ÿ q 2 which is proved as a theorem in the following sections. Note that for a single column matrix M, the multiplication between a quaternion and a matrix falls back to an ordinary quaternion multiplication. ô Row-wise (or vertical) quaternion-matrix multiplication Row-wise quaternion-matrix multiplication can be defined by transposing matrices and using column-wise multiplication: I q ÿ M ÿ q 2 M ô Conj[] matrix operator Conj[] operator conjugates only the vector part of the matrix, but not the scalar and the cross product parts: Conj@ q_? MatrixQ D := q P,T q P,2T q P,3T q P,4T q P2,T q P2,2T q P2,3T q P2,4T q P3,T q P3,2T q P3,3T q P3,4T q P4,T q P4,2T q P4,3T q P4,4T Note that Conj[] is not a conjugate of a quaternion! The conjugate of a quaternion corresponds to the plain transpose of the matrix.
4 4 Rotations of Tensors using Quaternions v0.3.nb Idea behind quaternion-matrix multiplication Basic idea behind the quaternion-matrix multiplication comes from an expression: 0, 0, 0 D 0, 0, 0 D q, for which it is assumed to be also valid when transforming identity matrix I4 that is equvalent representation of the quaternion [, 0, 0, 0 ] in matrix form, i.e. that there should be kind of multiplication where: H L q H L q, ü Theorems ô Quaternion-matrix multiplication Quaternion multiplication associativity (sanity-check): q q q2 q H q q2l Hq ql q2 Quaternion-matrix multiplication associativity theorems: q A q2 q H A q2 L H q A L q2 q A q H q A L q q A q q H A q L H q I4 L q I4 H I4 q L I4 q Hq I4L A q A HI4 ql A A q Hq I4L HI4 q2l q I4 q2
5 Rotations of Tensors using Quaternions v0.3.nb 5 HI4 q2l HHq I4L A L Hq AL q2 q A q2 HI4 q2l Hq I4L A q A q2 Ï HI4 q2l Hq I4L A q A q2 Ï Hq I4L HI4 q2l A q A q2 Ï Hq I4 q2l A q A q2 Conjugations of quaternion-matrix multiplications: Conj@ q A D Conj@AD q Conj@ A q D q Conj@AD Conj@ q A q2d Conj@ q HA q2ld Conj@A q2d q q2 Conj@AD q Conj@ q A q D q Conj@AD q ô Left and right rotation matrices Because the quaternion multiplication is bilinear, it can be expressed in a matrix form, and in two different ways ([]). Multiplication on the left qp gives L q ÿp where p is now treated as a 4-dimensional column vector p and multiplication on the right pq gives p ÿr q. Let us define L q and R q as quaternion-multiplication with identity matrix: ClearAll@ L, R, Q D L q_ := q I4 R q_ := I4 q Q q_ := L q R q Proof that the quaternion multiplication can be decomposed into L&R matrix parts and that L q R q * p corresponds to qpq * : AssertA L q R q ToVector@qD ToVector@ q q q D, ô Properties of L, R and Q matrices L q IL q M R q IR q M I q R q M A q I R q A M IL q A M q
6 6 Rotations of Tensors using Quaternions v0.3.nb Q q L q R q IR q L q M IR q L q M R q L q Q q R q L q == L q R q Q q Q q q 4 I4 ConjAL q A q ConjA q IL q AM êêassert ConjAL q A q ConjAL q A q ô L & R matrices components 9 "L q ", "L q ", "R q ", "R q " =, 9 L q êê Qorm, L q êê Qorm, R q êê Qorm, R q êê Qorm = = L q L q R q R q w x y z w x y z w x y z w x y z x w z y y z w x z y x w x w z y y z w x z y x w x w z y y z w x z y x w x w z y y z w x z y x w ô Rotation matrix Q components 9 "Q q L q R q R q L q ", "Q q L q R q R q L q " =, 9 SimplifyA Q q, êê Qorm, SimplifyA Q q, êê Qorm = = Q q L q R q R q L q Q q L q R q R q L q y 2 2z 2 2xy 2wz 2Hw y+ xzl 0 2y 2 2z 2 2 Hx y+ wzl 2 wy+ 2xz 0 2 Hx y+ wzl 2x 2 2z 2 2wx+ 2yz 0 2 x y 2wz 2x 2 2z 2 2 Hw x+ yzl 0 2 wy+ 2xz 2Hw x+ yzl 2x 2 2y Hw y+ xzl 2 wx+ 2yz 2x 2 2y 2 ô Shoemake's rotation matrix from quaternion Classical Shoemake's form of the rotation matrix from a quaternion is given as (a transformation matrix based on the Cayley- Klein parameters; see [2]): RotationMatrix4@ q D êê Qorm Iy 2 + z 2 M 2xy 2wz 2wy+ 2xz 0 2xy+ 2wz 2 Ix 2 + z 2 M 2 wx+ 2yz 0 2 wy+ 2xz 2wx+ 2yz 2 Ix 2 + y 2 M
7 Rotations of Tensors using Quaternions v0.3.nb 7 ô Scalar, vector and cross-product parts of L&R matrices Scalar part: 2 IL q + L q M 2 IR q + R q M Vector part: w w w w 2 IL q R qm == 2 IR q L qm 0 x y z x y z IL q R q M == 2 IR q L q M 0 x y z x y z Cross product part: 2 IL q R q M z y 0 z 0 x 0 y x 0 2 IL q R q M z y 0 z 0 x 0 y x 0 Q Q + 2 Q (vector part) form. Note that the vector part is a pure imaginary part of a quaternion. q q D q 2 Q q == L q L q + 2 L q 0 x y z x y z L q L q + 2 L q 2 IR q L qm L q R q 2 L q IL q L q Q q M == 2 IL q R q M
8 8 Rotations of Tensors using Quaternions v0.3.nb quality of rotation matrix from quaternion and uler-rodrigues' formula uler-rodrigues' formula for the rotation matrix corresponding to a rotation by an angle q about a fixed axis specified by the unit vector u. Reference: ulerrodrigues@ θ_, 8 x_, y_, z_ <D:= Cos@θD H Cos@θDL xx yx zx xy yy zy xz yz zz + Sin@θD 0 z y z 0 x y x 0 Proof that the rotation matrix from a quaternion (as defined in Quat.m package) is equivalent to uler-rodrigues' formula. BlockA 8 x, y, z, θ, q<, q = To $AngleAxis@ θ, 8 x, y, z <D; AssertA RotationMatrix@ q D ulerrodrigues@ θ, 8 x, y, z <D, Assumptions π θ πïx 2 + y 2 + z 2 Ï8x, y, z, θ < Reals ô Rotation matrix theorems RM = RotationMatrix4@ q D; RM êê Qorm Iy 2 + z 2 M 2xy 2wz 2wy+ 2xz 0 2xy+ 2wz 2 Ix 2 + z 2 M 2 wx+ 2yz 0 2 wy+ 2xz 2wx+ 2yz 2 Ix 2 + y 2 M RM Q q L q R q Ì RM Q q R q L q, AssertA RM ToVector@qD L q R q ToVector@qD ToVector@ q q q D,
9 Rotations of Tensors using Quaternions v0.3.nb 9 RM A RM Q q A Q q L q R q A L q R q J q H q A q L q N, I4 Hq I4 q L I4 Hq I4 q L IL q R q M I4 IL q R q M L q R q I4 R q L q L q IR q I4 R q M L q, RM A RM Hq I4 q L A Hq I4 q L IL q R q M A IL q R q M L q R q A R q L q L q IR q A R q M L q, Proof that both Q q and R M are orthogonal matrices: Q q Q q IL q R q M IL q R q M L q R q IR q M IL q M L q R q R q L q I4,
10 0 Rotations of Tensors using Quaternions v0.3.nb AssertA Q q Q q Q q Q q RM RM I4, ô Quaternionic meaning of RHqL A RHqL T The hint how quaternion multplication affects a (moment of inertia or any other) tensor A can be seen from the earlier proven theorem: RM A RM J q H q A q L q N, which shows that the tensor A is rotated, at first, by rotating quaternions across each column A c = q ÿ A ÿ q *, and at second, by rotating quaternions across each row A rc = Iq ÿ HA c L ÿ q * M. q ÿ HL q ÿ HL q q 2 q 3 q 4 q 2 q 22 q 23 q 24 q 3 q 32 q 33 q 34 q 4 q 42 q 43 q 44 HL ÿ q * HL ÿ q * The meaning of the compound multiplications across columns and rows is extracting and affecting individual components of quaternions embedded from the previous rotations of the tensor. The moment of inertia tensor can always be reduced to its principal axes in diagonal form. In a matrix representation of a quaternion, the diagonal of the matrix contains the repeated scalar (real) part of the quaternion. (or any quaternion q there exist quaternion p such that qp is scalar.) The principal moment of inertia, on the other hand, has on its diagonal three different scalars in general case, which gives origin to three different quaternions that are embedded and independently shuffled across-columns/rows during rotations. ü References [] berly, David H. with a contrib. by Shoemake, Ken - Game Physics, Morgan Kaufmann, 2004, Chapter 0 "Quaternions", pp , [2] Shoemake, Ken - Uniform random rotations, in: D. Kirk (d.), Graphics Gems III, Academic Press, London, 992, pp , [3] Shoemake, Ken - Quaternions, Department of Computer and Information Science, University of Pennsylvania, , downloaded from C+MS course CS7 at Caltech:
Determination of Locally Varying Directions through Mass Moment of Inertia Tensor
Determination of Locally Varying Directions through Mass Moment of Inertia Tensor R. M. Hassanpour and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering
More informationCourse 2BA1: Hilary Term 2007 Section 8: Quaternions and Rotations
Course BA1: Hilary Term 007 Section 8: Quaternions and Rotations David R. Wilkins Copyright c David R. Wilkins 005 Contents 8 Quaternions and Rotations 1 8.1 Quaternions............................ 1 8.
More informationProgramming Project 2: Harmonic Vibrational Frequencies
Programming Project 2: Harmonic Vibrational Frequencies Center for Computational Chemistry University of Georgia Athens, Georgia 30602 Summer 2012 1 Introduction This is the second programming project
More informationManipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA
Manipulator Dynamics 2 Forward Dynamics Problem Given: Joint torques and links geometry, mass, inertia, friction Compute: Angular acceleration of the links (solve differential equations) Solution Dynamic
More informationCOMP 175 COMPUTER GRAPHICS. Lecture 04: Transform 1. COMP 175: Computer Graphics February 9, Erik Anderson 04 Transform 1
Lecture 04: Transform COMP 75: Computer Graphics February 9, 206 /59 Admin Sign up via email/piazza for your in-person grading Anderson@cs.tufts.edu 2/59 Geometric Transform Apply transforms to a hierarchy
More informationWhy Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms
Why Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms ITCS 3050:Game Engine Programming 1 Geometric Transformations Implementing
More informationProblem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1
Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate
More informationCourse MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions
Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions David R. Wilkins Copyright c David R. Wilkins 2000 2010 Contents 4 Vectors and Quaternions 47 4.1 Vectors...............................
More informationOrientation of Piezoelectric Crystals and Acoustic Wave Propagation
Orientation of Piezoelectric Crystals and Acoustic Wave Propagation Guigen Zhang Department of Bioengineering Department of Electrical and Computer Engineering Institute for Biological Interfaces of Engineering
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationRigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient
Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained
More informationClosed-Form Solution Of Absolute Orientation Using Unit Quaternions
Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding
More informationGroup, Rings, and Fields Rahul Pandharipande. I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S,
Group, Rings, and Fields Rahul Pandharipande I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S, A binary operation φ is a function, S S = {(x, y) x, y S}. φ
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
More informationMultipole moments. Dipole moment. The second moment µ is more commonly called the dipole moment, of the charge. distribution and is a vector
Dipole moment Multipole moments The second moment µ is more commonly called the dipole moment, of the charge distribution and is a vector µ = µ x ˆx + µ y ŷ + µ z ẑ where the α component is given by µ
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS. William Rowan Hamilton
ON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS By William Rowan Hamilton (Proceedings of the Royal Irish Academy, 2 (1844), 424 434.) Edited by David R. Wilkins 1999 On
More informationProblem Set 2 Due Thursday, October 1, & & & & # % (b) Construct a representation using five d orbitals that sit on the origin as a basis:
Problem Set 2 Due Thursday, October 1, 29 Problems from Cotton: Chapter 4: 4.6, 4.7; Chapter 6: 6.2, 6.4, 6.5 Additional problems: (1) Consider the D 3h point group and use a coordinate system wherein
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 16: Euler s Equations Attitude Dynamics In this Lecture we will cover: The Problem of Attitude Stabilization Actuators Newton
More informationHomework 1/Solutions. Graded Exercises
MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both
More informationExercise 1: Inertia moment of a simple pendulum
Exercise : Inertia moment of a simple pendulum A simple pendulum is represented in Figure. Three reference frames are introduced: R is the fixed/inertial RF, with origin in the rotation center and i along
More informationQuaternions. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Quaternions Semester 1, / 40
Quaternions Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-2017 B. Bona (DAUIN) Quaternions Semester 1, 2016-2017 1 / 40 Introduction Complex numbers with unit norm can be used as rotation operators
More informationDeep Learning Book Notes Chapter 2: Linear Algebra
Deep Learning Book Notes Chapter 2: Linear Algebra Compiled By: Abhinaba Bala, Dakshit Agrawal, Mohit Jain Section 2.1: Scalars, Vectors, Matrices and Tensors Scalar Single Number Lowercase names in italic
More informationMatrix Operations. Linear Combination Vector Algebra Angle Between Vectors Projections and Reflections Equality of matrices, Augmented Matrix
Linear Combination Vector Algebra Angle Between Vectors Projections and Reflections Equality of matrices, Augmented Matrix Matrix Operations Matrix Addition and Matrix Scalar Multiply Matrix Multiply Matrix
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More information1. Vectors.
1. Vectors 1.1 Vectors and Matrices Linear algebra is concerned with two basic kinds of quantities: vectors and matrices. 1.1 Vectors and Matrices Scalars and Vectors - Scalar: a numerical value denoted
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationCSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer
CSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2016 Announcements Monday October 3: Discussion Assignment
More informationAppendix Composite Point Rotation Sequences
Appendix Composite Point Rotation Sequences A. Euler Rotations In Chap. 6 we considered composite Euler rotations comprising individual rotations about the x, y and z axes such as R γ,x R β,y R α,z and
More informationLinear Algebra. Chapter 8: Eigenvalues: Further Applications and Computations Section 8.2. Applications to Geometry Proofs of Theorems.
Linear Algebra Chapter 8: Eigenvalues: Further Applications and Computations Section 8.2. Applications to Geometry Proofs of Theorems May 1, 2018 () Linear Algebra May 1, 2018 1 / 8 Table of contents 1
More informationPoint cloud to point cloud rigid transformations. Minimizing Rigid Registration Errors
Pont cloud to pont cloud rgd transformatons Russell Taylor 600.445 1 600.445 Fall 000-015 Mnmzng Rgd Regstraton Errors Typcally, gven a set of ponts {a } n one coordnate system and another set of ponts
More information7.4 Rigid body simulation*
7.4 Rigid body simulation* Quaternions vs. Rotation Matrices - unit quaternions, a better way to represent the orientation of a rigid body. Why? (1) More compact: 4 numbers vs. 9 numbers (2) Smooth transition
More informationCSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer
CSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer Jürgen P. Schulze, Ph.D. University of California, San Diego Spring Quarter 2016 Announcements Project 1 due next Friday at
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More informationAn Introduction to Matrix Algebra
An Introduction to Matrix Algebra EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #8 EPSY 905: Matrix Algebra In This Lecture An introduction to matrix algebra Ø Scalars, vectors, and matrices
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More information6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant
DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional
More informationHilbert s Metric and Gromov Hyperbolicity
Hilbert s Metric and Gromov Hyperbolicity Andrew Altman May 13, 2014 1 1 HILBERT METRIC 2 1 Hilbert Metric The Hilbert metric is a distance function defined on a convex bounded subset of the n-dimensional
More informationEEL6667: Homework #1 Solutions
EEL6667: Homework #1 Solutions Problem 1: Note: homework1.nb is a Mathematica notebook that solves many of the problems in this homework. (a) See homework1.nb. (b) See homework1.nb. Problem :[raig, Exercise.14
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationOctonions. Robert A. Wilson. 24/11/08, QMUL, Pure Mathematics Seminar
Octonions Robert A. Wilson 4//08, QMUL, Pure Mathematics Seminar Introduction This is the second talk in a projected series of five. I shall try to make them as independent as possible, so that it will
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More informationUseful Formulae ( )
Appendix A Useful Formulae (985-989-993-) 34 Jeremić et al. A.. CHAPTER SUMMARY AND HIGHLIGHTS page: 35 of 536 A. Chapter Summary and Highlights A. Stress and Strain This section reviews small deformation
More informationHong Kong Institute of Vocational Education (Tsing Yi) Higher Diploma in Civil Engineering Structural Mechanics. Chapter 2 SECTION PROPERTIES
Section Properties Centroid The centroid of an area is the point about which the area could be balanced if it was supported from that point. The word is derived from the word center, and it can be though
More informationCHAPTER 3 BOOLEAN ALGEBRA
CHAPTER 3 BOOLEAN ALGEBRA (continued) This chapter in the book includes: Objectives Study Guide 3.1 Multiplying Out and Factoring Expressions 3.2 Exclusive-OR and Equivalence Operations 3.3 The Consensus
More informationInverse of a Square Matrix. For an N N square matrix A, the inverse of A, 1
Inverse of a Square Matrix For an N N square matrix A, the inverse of A, 1 A, exists if and only if A is of full rank, i.e., if and only if no column of A is a linear combination 1 of the others. A is
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationChapter 5: Application of Filters to Potential Field Gradient Tensor Data
Chapter 5: Filters 98 Chapter 5: Application of Filters to Potential Field Gradient Tensor Data 5.1 Introduction One of the objectives of this research is to investigate the use of spatial filters on potential
More informationLagrange Multipliers
Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationSubalgebras of the Split Octonions
Subalgebras of the Split Octonions Lida Bentz July 27, 2017 CONTENTS CONTENTS Contents 1 Introduction 3 1.1 Complex numbers......................... 3 1.2 Quaternions............................ 3 1.3
More informationCGAlgebra: a Mathematica package for conformal geometric algebra
arxiv:1711.02513v2 [cs.ms] 23 Aug 2018 CGAlgebra: a Mathematica package for conformal geometric algebra José L. Aragón Centro de Física Aplicada y Tecnología Avanzada, Universidad Nacional Autónoma de
More informationEigen decomposition for 3D stress tensor
Eigen decomposition for 3D stress tensor by Vince Cronin Begun September 13, 2010; last revised September 23, 2015 Unfinished Dra ; Subject to Revision Introduction If we have a cubic free body that is
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 00 Dec 27, 2014.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim100887@aol.com rev 00 Dec 27, 2014 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationAB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT
FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be
More informationComputer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design
CVEN 302-501 Computer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design Date distributed : 12/2/2015 Date due : 12/9/2015 at
More informationReview of Coordinate Systems
Vector in 2 R and 3 R Review of Coordinate Systems Used to describe the position of a point in space Common coordinate systems are: Cartesian Polar Cartesian Coordinate System Also called rectangular coordinate
More informationChapter 4 The Equations of Motion
Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations
More informationRigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013
Rigid body dynamics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Rigid body dynamics October 2013 1 / 16 Multiple point-mass bodies Each mass is
More informationSignal Processing s Importance in Manufacturing of a Special Device
Signal Processing s Importance in Manufacturing of a Special Device MIHAIELA ILIESCU 1, BRINDUŞ COMĂNESCU, EMIL NUŢU 1 1 Manufacturing Department; POLITEHNICA University of Bucharest Splaiul Independenţei
More informationIntroduction. How to use this book. Linear algebra. Mathematica. Mathematica cells
Introduction How to use this book This guide is meant as a standard reference to definitions, examples, and Mathematica techniques for linear algebra. Complementary material can be found in the Help sections
More informationPolyhedral Mass Properties (Revisited)
Polyhedral Mass Properties (Revisited) David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationSimplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System
NASA Technical Memorandum 112868 Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System Nelson J. Groom Langley Research Center, Hampton, Virginia June 1997 National
More informationreview To find the coefficient of all the terms in 15ab + 60bc 17ca: Coefficient of ab = 15 Coefficient of bc = 60 Coefficient of ca = -17
1. Revision Recall basic terms of algebraic expressions like Variable, Constant, Term, Coefficient, Polynomial etc. The coefficients of the terms in 4x 2 5xy + 6y 2 are Coefficient of 4x 2 is 4 Coefficient
More information. ffflffluary 7, 1855.
x B B - Y 8 B > ) - ( vv B ( v v v (B/ x< / Y 8 8 > [ x v 6 ) > ( - ) - x ( < v x { > v v q < 8 - - - 4 B ( v - / v x [ - - B v B --------- v v ( v < v v v q B v B B v?8 Y X $ v x B ( B B B B ) ( - v -
More informationProperties of the stress tensor
Appendix C Properties of the stress tensor Some of the basic properties of the stress tensor and traction vector are reviewed in the following. C.1 The traction vector Let us assume that the state of stress
More informationMatrix Algebra: Vectors
A Matrix Algebra: Vectors A Appendix A: MATRIX ALGEBRA: VECTORS A 2 A MOTIVATION Matrix notation was invented primarily to express linear algebra relations in compact form Compactness enhances visualization
More informationMatrices. Background mathematics review
Matrices Background mathematics review David Miller Matrices Matrix notation Background mathematics review David Miller Matrix notation A matrix is, first of all, a rectangular array of numbers An M N
More informationProperties of surfaces II: Second moment of area
Properties of surfaces II: Second moment of area Just as we have discussing first moment of an area and its relation with problems in mechanics, we will now describe second moment and product of area of
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationContents. Appendix D (Inner Product Spaces) W-51. Index W-63
Contents Appendix D (Inner Product Spaces W-5 Index W-63 Inner city space W-49 W-5 Chapter : Appendix D Inner Product Spaces The inner product, taken of any two vectors in an arbitrary vector space, generalizes
More informationQuantum Information & Quantum Computing
Math 478, Phys 478, CS4803, February 9, 006 1 Georgia Tech Math, Physics & Computing Math 478, Phys 478, CS4803 Quantum Information & Quantum Computing Problems Set 1 Due February 9, 006 Part I : 1. Read
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More information2.003J Spring 2011: Dynamics and Control I On Notation Massachusetts Institute of Technology Department of Mechanical Engineering Feb.
2.J Spring 211: Dynamics and Control I On Notation Department of Mechanical Engineering Feb. 8, 211 On Notation Much of 2.j involves defining then manipulating points, frames of reference, and vectors.
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More information- 1 - Items related to expected use of technology appear in bold italics.
- 1 - Items related to expected use of technology appear in bold italics. Operating with Geometric and Cartesian Vectors Determining Intersections of Lines and Planes in Three- Space Similar content as
More informationMoments of inertia of a cross section
BUDAPEST UNVERSTY OF TECHNOLOGY AND ECONOMCS Moments of inertia of a cross section Made by: David Lehotzky Budapest, March 3, 206. Figure : The cross section under study a b h v d 0 8 4 0.5 5 Table : Parameters
More informationGeometric Fundamentals in Robotics Quaternions
Geometric Fundamentals in Robotics Quaternions Basilio Bona DAUIN-Politecnico di Torino July 2009 Basilio Bona (DAUIN-Politecnico di Torino) Quaternions July 2009 1 / 47 Introduction Quaternions were discovered
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationComplex Variables. Chapter 1. Complex Numbers Section 1.2. Basic Algebraic Properties Proofs of Theorems. December 16, 2016
Complex Variables Chapter 1. Complex Numbers Section 1.2. Basic Algebraic Properties Proofs of Theorems December 16, 2016 () Complex Variables December 16, 2016 1 / 12 Table of contents 1 Theorem 1.2.1
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationGG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS
GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS I Main Topics A Why deal with tensors? B Order of scalars, vectors, and tensors C Linear transformation of scalars and vectors (and tensors) II Why
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
More informationVectors and Matrices
Vectors and Matrices Scalars We often employ a single number to represent quantities that we use in our daily lives such as weight, height etc. The magnitude of this number depends on our age and whether
More informationInverses of Square Matrices
Inverses of Square Matrices A. Havens Department of Mathematics University of Massachusetts, Amherst February 23-26, 2018 Outline 1 Defining Inverses Inverses for Products and Functions Defining Inverse
More informationVector Calculus - GATE Study Material in PDF
Vector Calculus - GATE Study Material in PDF In previous articles, we have already seen the basics of Calculus Differentiation and Integration and applications. In GATE 2018 Study Notes, we will be introduced
More informationGame Engineering: 2D
Game Engineering: 2D CS420-2010F-06 2D Math David Galles Department of Computer Science University of San Francisco 06-0: Back to Basics A Vector is a displacement Vector has both direction and length
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationA Learning Progression for Complex Numbers
A Learning Progression for Complex Numbers In mathematics curriculum development around the world, the opportunity for students to study complex numbers in secondary schools is decreasing. Given that the
More informationAlgebraic Expressions and Identities
ALGEBRAIC EXPRESSIONS AND IDENTITIES 137 Algebraic Expressions and Identities CHAPTER 9 9.1 What are Expressions? In earlier classes, we have already become familiar with what algebraic expressions (or
More informationELEMENTARY MATRIX ALGEBRA
ELEMENTARY MATRIX ALGEBRA Third Edition FRANZ E. HOHN DOVER PUBLICATIONS, INC. Mineola, New York CONTENTS CHAPTER \ Introduction to Matrix Algebra 1.1 Matrices 1 1.2 Equality of Matrices 2 13 Addition
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More information