Variations of Augmented Lagrangian for Robotic Multi-Contact Simulation