Automatic generation and numerical integration of differential-algebraic equations of multibody dynamics (Q2368263)

An algorithm for computing the matrix and vector functions in the numerical integration method of constrained equations of motion in multibody mechanical system dynamics is presented. Cartesian coordinates and Euler parameters have been used as the state variables to construct the mathematical model for a spatial multibody mechanical system. The efficiency and accuracy of computing analytic Jacobian matrix are demonstrated by comparing execution time and numerical results with those obtained by applying a first order linearization method. The proposed computational technique is also suitable for various formulations of the equations of motion, for instance, the Euler angles instead of the Euler parameters can be used for the orientation variables. The final goal is to construct a standard library of joint constraints and force expressions for automatic generation of constrained equations of motion and to enhance the numerical integration methods for multibody dynamics.