Finite-element analysis of contact between elastic self-affine surfaces

S. Hyun,L. Pei,J.-F. Molinari,M. O. Robbins

DOI: https://doi.org/10.1103/PhysRevE.70.026117

2004-07-08

Abstract:Finite element methods are used to study non-adhesive, frictionless contact between elastic solids with self-affine surfaces. We find that the total contact area rises linearly with load at small loads. The mean pressure in the contact regions is independent of load and proportional to the rms slope of the surface. The constant of proportionality is nearly independent of Poisson ratio and roughness exponent and lies between previous analytic predictions. The contact morphology is also analyzed. Connected contact regions have a fractal area and perimeter. The probability of finding a cluster of area $a_c$ drops as $a_c^{-\tau}$ where $\tau $ increases with decreasing roughness exponent. The distribution of pressures shows an exponential tail that is also found in many jammed systems. These results are contrasted to simpler models and experiment.

Materials Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to study the relationship between the contact area and the load, as well as the morphological characteristics of the contact area in the finite - element analysis of non - adhesive and frictionless contact between elastic self - affine surfaces. Specifically, the paper focuses on the following aspects:

1. **Relationship between contact area and load**:
- The paper studies the linear relationship between the contact area \( A \) and the load \( W \) under small loads by the finite - element method. It is found that under small loads, the contact area does increase linearly with the load, and this linear relationship remains stable within a certain range.
- The author also explores that when the proportion of the contact area \( A \) to the total area exceeds 4% or 5%, the relationship between the contact area and the load begins to deviate from linearity.

2. **Average contact pressure**:
- It is found that the average pressure \( \langle p \rangle \equiv \frac{W}{A} \) in the contact area is independent of the load, but is proportional to the root - mean - square slope \( \sqrt{\langle |\nabla h|^2 \rangle} \) of the surface. The proportional constant hardly depends on the Poisson ratio and the roughness index.

3. **Contact morphology**:
- The paper analyzes the morphological characteristics of the connected contact areas and finds that these areas have fractal areas and perimeters. The area distribution \( P(a_c) \) of the contact areas follows a power - law distribution \( P(a_c) \propto a_c^{-\tau} \), where \( \tau \) increases as the roughness index decreases.
- For the roughness index \( H < 0.9 \), the power - law index \( \tau \) is greater than 2, which means that the average contact size is always close to the calculated resolution.

4. **Pressure distribution**:
- The pressure distribution shows an exponential tail, which has also been observed in many pressurized systems. It is found that the pressure distributions under different system sizes, roughness amplitudes and roughness indices can collapse into a universal curve when normalized to their average values.
- The average pressure \( \frac{W}{A} \) can be obtained from the roughness amplitude, the pressure distribution decreases monotonically with the increase of pressure, and decays exponentially at high pressures.

5. **Comparison with theoretical models**:
- The results of the paper are compared with previous theoretical models, especially with the predictions of Bush et al. and Persson. It is found that the proportional constant \( \kappa \) in the actual results is between these two theoretical predictions.
- For different roughness indices \( H \) and Poisson ratios \( \nu \), the change of the proportional constant \( \kappa \) is relatively small, indicating that these parameters have a limited impact on the relationship between the contact area and the load.

In conclusion, through detailed numerical simulations and analyses, this paper provides new insights into the contact behavior between elastic self - affine surfaces, especially in - depth understanding of the relationship between the contact area and the load, contact morphology and pressure distribution. These results are helpful for improving existing theoretical models and providing more accurate predictions for contact mechanics problems in practical applications.