An augmented Lagrangian approach to isogeometric analysis of 3D large deformation contact problems