To impose the displacement compatibility condition along the contact area, additional Lagrange multiplier degrees of freedom, , are introduced to the system to minimize the potential function under the contact constraints defined by the function . By adding the contact conditions to the usual variational indicator, the following function is obtained
The friction part of the force transmissibility condition is not enforced yet in Equation (3.5). This condition is handled separately after each iteration and the appropriate boundary conditions are set for the next iteration. This approach causes the status of the boundary conditions applied on the system to change after each iteration. As a result, the quadratic convergence of the Newton-Rapson method may be seriously affected.