The Neumann-Neumann domain decomposition algorithm is extended to a system which consists of two plane elastic bodies in mutual contact taking into account friction obeying the Coulomb law. The main difficulty is due to highly nonlinear unilateral and friction conditions which are prescribed on contact interfaces.
A fixed point approach is used to ensure continuity of contact stresses. Numerical experiments illustrate that there exists an optimal relaxation parameter which is nearly independent of the value of friction coefficient and the mesh norm.