Constrained multibody kinematics and dynamics in absolute coordinates: a discussion of three approaches to representing rigid body rotation
share
We compare three approaches to posing the index 3 set of differential algebraic equations (DAEs) associated with the constrained multibody dynamics problem formulated in absolute coordinates. The first approach works directly with the orientation matrix and therefore eschews the need for generalized coordinates used to produce the orientation matrix 𝐀. The approach is informed by the fact that rotation matrices belong to the SO(3) Lie matrix group. The second approach employs Euler parameters, while the third uses Euler angles. In all cases, the index 3 DAE problem is solved via a first order implicit numerical integrator. We note a roughly twofold speedup of 𝐫𝐀 over 𝐫ϵ, and a 1.2 – 1.3 times speedup of 𝐫ϵ over 𝐫𝐩. The tests were carried out in conjunction with four 3D mechanisms. The improvements in simulation speed of the 𝐫𝐀 approach are traced back to a simpler form of the equations of motion and more concise Jacobians that enter the numerical solution. The contributions made herein are twofold. First, we provide first order variations of all the quantities that enter the 𝐫𝐀 formulation when used in the context of implicit integration; i.e., sensitivity of the kinematic constraints for all lower pair joints, as well as the sensitivity of the constraint reaction forces. Second, to the best of our knowledge, there is no other contribution that compares head to head the solution efficiency of 𝐫𝐀, 𝐫𝐩, and 𝐫ϵ in the context of the multibody dynamics problem posed in absolute coordinates.
READ FULL TEXT