Add PGS solver
Use DanzigLCPSolver per default
Fix contact constraint for weld joints