Premium
Rigid body dynamics in terms of quaternions: Hamiltonian formulation and conserving numerical integration
Author(s) -
Betsch Peter,
Siebert Ralf
Publication year - 2009
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.2586
Subject(s) - quaternion , discretization , hamiltonian (control theory) , mathematics , variational integrator , rigid body dynamics , rigid body , sylvester's law of inertia , equations of motion , classical mechanics , mathematical analysis , integrator , physics , geometry , mathematical optimization , symmetric matrix , quantum mechanics , eigenvalues and eigenvectors , voltage
In the present paper unit quaternions are used to describe the rotational motion of a rigid body. The unit‐length constraint is enforced explicitly by means of an algebraic constraint. Correspondingly, the equations of motion assume the form of differential‐algebraic equations (DAEs). A new route to the derivation of the mass matrix associated with the quaternion formulation is presented. In contrast to previous works, the newly proposed approach yields a non‐singular mass matrix. Consequently, the passage to the Hamiltonian framework is made possible without the need to introduce undetermined inertia terms. The Hamiltonian form of the DAEs along with the notion of a discrete derivative make possible the design of a new quaternion‐based energy–momentum scheme. Two numerical examples demonstrate the performance of the newly developed method. In this connection, comparison is made with a quaternion‐based variational integrator, a director‐based energy–momentum scheme, and a momentum conserving scheme relying on the discretization of the classical Euler's equations. Copyright © 2009 John Wiley & Sons, Ltd.