Finite deformation effects are partially taken into account in a more advanced approach in which the contact pressure is com- puted for a fully nonlinear frictionless contact problem, typically using the finite-element method. Subsequently, a (linear) influ- ence coefficient matrix is obtained from off-line finite-element computations, e.g., employing a nodal perturbation technique, and this matrix is used in the EHL solver [17,18]. In other words, in this approach, the nonlinear behaviour of the deforming solid is linearized at the deformed state determined by solving the contact problem. As a drawback, the friction stresses are neglected in this approach, as they are not known a priori and thus cannot be included in the contact analysis.