Singularity-free augmented Lagrangian algorithms for constrained multibody dynamics