A Quasicoordinate Formulation For DynamicSimulation Of Complex Multibody Systems WithConstraints