The finite element mesh comprises almost 100,000 hexahedral 8-node F-bar elements [28], and the total of 480,181 unknowns including nodal displacements in the body, pressures and void fractions on the lubricated contact surface and the Lagrange multipliers associated with the hanging-node mesh refinement technique.