Elastic Response of Rough Surfaces in Partial Contact

We model numerically the partial normal contact of two elastic rough surfaces with highly correlated asperities. Facing surfaces are unmated and described as self-affine with a Hurst exponent H. The numerical algorithm is based on Fourier acceleration and allows efficient simulation of very large systems. The force, F, versus contact area, A, characteristics follows the law F = c A^[(1+H)/2] in accordance with the suggestion of Roux et al. (Europhys. Lett. 23, 277 (1993)). However finite size corrections are very large even for 512X512 systems where the effective exponent is still 20% larger than its asymptotic value.