The numerical techniques are presented to solve the static and dynamic contact problems of deformable bodies having large rotations of the contact surfaces. The contact conditions on the possible contact surfaces are enforced by using the contact error vector, and an iterative scheme similar to augmented Lagrange multiplier method is employed to reduce the contact error vector monotonically. For dynamic contact problems using implicit time integration, a contact error vector is also defined by combining the displacement, velocity, and acceleration on the contact surface. The suggested iterative technique is implemented to ABAQUS by using the UEL subroutine UEL. In this work, after the computing procedures to solve the frictional contact problems are explained, the numerical examples are presented to compare the present solutions with those obtained by ABAQUS.