On the Numerical Solution of Many-Body Contact Dynamics Problems Formulated as Complementarity Problems