A SYMPLECTIC ALGORITHM FOR DYNAMICS OF RIGID BODY

For the dynamics of a rigid body with a fixed point based on the quaternion and the corresponding generalized momenta,a displacement-based symplectic integration scheme for differential-algebraic equations is proposed and applied to the Lagrange’s equa- tions based on dependent generalized momenta.Numerical experiments show that the algorithm possesses such characters as high precision and preserving system invariants. More importantly,the generalized momenta based Lagrange’s equations show unique ad- vantages over the traditional Lagrange’s equations in symplectic integrations.