Abstract The simulation of mathematical models of mechanical systems with closed kinematic chains involves the solution to a system of highly coupled differential-algebraic equations. The numerical stiffness of these systems calls for small time steps in order to insure accuracy. Real-time and interactive forward simulations tend to be difficult to achieve for such systems, especially for large...