Augmented lagrangian and mass-orthogonal projection methods for constrained multibody dynamics