A numerical integrator is derived for a class of models that describe the attitude dynamics of a rigid body in the presence of an attitude dependent potential. The numerical integrator is obtained from a discrete variational principle, and exhibits excellent geometric conservation properties. In particular, by performing computations at the level of the Lie algebra, and updating the solution using the matrix exponential, the attitude automatically evolves on the rotation group embedded in the space of matrices. The geometric conservation properties of the numerical integrator imply long time numerical stability. We apply this variational integrator to the uncontrolled 3D pendulum, that is a rigid asymmetric body supported at a frictionless pivot acting under the influence of uniform gravity. Interesting dynamics of the 3D pendulum are exposed
展开▼
