A constrained variable metric
method for solving canonical equations in constrained multibody systems dynamics is presented.
The framwork of which is based on rewriting the equations as the optimization problem and using
symplectic integral. The algorithm keep calculation stable for time simulation. Numerical
examples show that algorithm is numerically efficient.