An index 0 differential-algebraic equation formulation for multibody dynamics: Nonholonomic constraints

A method is presented for formulating and numerically integrating index 0 differential-algebraic equations of motion for multibody systems with holonomic and nonholonomic constraints. Tangent space coordinates are defined in configuration and velocity spaces as independent generalized coordinates that serve as state variables in the formulation. Orthogonal dependent coordinates and velocities are used to enforce position, velocity, and acceleration constraints to within specified error tolerances. Explicit and implicit numerical integration algorithms are presented and used in solution of three examples: one planar and two spatial. Numerical results verify that accurate results are obtained, satisfying all three forms of kinematic constraint to within error tolerances embedded in the formulation.