Contact mechanics models based on linearity assumptions, often using the viscoelastic half space theory and numerically implemented with boundary element method, are known to provide accurate results for small mean square slope of surface roughness. For large slope, accounting finite deformations, non-linear more but lead a prohibitive computational cost. We propose new hybrid multiscale approa...