A Hybrid Multiscale Approach for Rubber Contact

Contact mechanics models based on linearity assumptions, often using the viscoelastic half space theory and numerically implemented with the boundary element method, are known to provide accurate results for small mean square slope of the surface roughness. For large mean square slope, models accounting for finite deformations, often implemented with the non-linear finite element method, are more accurate but lead to a prohibitive computational cost. We propose a new hybrid multiscale approach able to account for the finite deformations arising due to large mean square slope, while keeping a computational cost similar to that associated to linear approaches. The basic strategy is a decomposition of the surface roughness power spectrum into a discrete number of waves, whose spectral range is partitioned into a high mean square slope range and a low mean square slope range. The contact mechanics in the former is accurately solved with the kinematically non-linear model and the results averaged out at the larger wavelength scale in terms of an effective interface interaction law. This law is then applied in the linear simulation involving the scales within the low mean square slope range. The proposed approach is a more accurate alternative to fully linear and a computationally faster alternative to fully non-linear contact mechanics approaches.