For the numerical solution of frictional dynamic contact problems, correct contact points and displacements are determined by iteratively reducing the displacement error vector monotonically toward zero And spurious oscillations are prevented from the solution by enforcing the velocity and acceleration compatibilities of the contact points with the corresponding error vectors. Numerical simulations are conducted to demonstrate the accuracy of the solution and the necessity of the velocity and acceleration compatibilities on the contact surface.