A stochastic model for elastoplastic contact of rough surfaces incorporating scale-dependent hardness

This paper presents a novel stochastic model based on compounded Chapman-Kolmogorov equations to analyze elastoplastic contact of rough surfaces by incorporating scale-dependent hardness, thereby predicting contact status evolution across scales and offering new insights for multidisciplinary fields involving multiscale roughness.

The Gist

Imagine two rough surfaces, like two pieces of sandpaper or a tire on a road, being pressed together. Even though they look flat from a distance, if you zoom in, they are actually covered in tiny mountains and valleys. When you push them together, only the very tips of these "mountains" (called asperities) actually touch.

This paper is about understanding exactly what happens at those tiny touch-points, specifically when the materials are soft enough to squish (plastic deformation) but also spring back a little (elasticity).

Here is the breakdown of their discovery using simple analogies:

1. The Problem: The "Zooming" Puzzle

Most old models for how surfaces touch assume that the material's hardness is the same no matter how small you look. But in reality, if you look at a tiny spot, the material often acts harder than if you look at a big spot. This is called the size effect.

Think of it like a crowd of people. If you look at a whole stadium, it's easy to move through. But if you zoom in on just three people standing shoulder-to-shoulder, it's very hard to squeeze through. The "crowd" (the material) feels harder when you look at it up close.

The authors wanted to build a model that accounts for this changing hardness as you zoom in and out of the surface.

2. The Solution: A "Probability Flow" Map

Instead of trying to simulate every single tiny mountain (which would take a supercomputer forever), the authors used a clever mathematical trick called the Chapman-Kolmogorov equation.

The Analogy: Imagine a river flowing downstream.

• The Water: Represents the "pressure" at the contact points.
• The Riverbed: Represents the roughness of the surface.
• The Banks: Represent the limits of the material. If the water gets too high, it spills over the bank (the material yields or squishes).

In the past, scientists thought the water could only flow one way: from a calm pool into a rushing rapid, and once it hit the bank, it stayed there. They assumed that once a spot on the surface got squished (plastic), it would stay squished no matter how much you zoomed in.

The New Discovery: The authors found that when hardness changes with scale, the water can actually flow backward.

• If a spot was squished at a low zoom level, but the material gets harder as you zoom in, that spot might actually "un-squish" and become elastic again.
• They created a map that tracks how the probability of a spot being "touching," "squished," or "not touching" changes as you zoom in and out.

3. The Three Zones of Contact

Using their new math, they identified three distinct states for how two rough surfaces interact, depending on the material properties and the "roughness" of the surface:

1. The "Springy" Zone (Linear Elastic): The surfaces touch, but they act like stiff springs. They squish a little and bounce back. This happens when the material is very hard at small scales.
2. The "Muddy" Zone (Fully Plastic): The surfaces are so soft or the pressure is so high that the mountains just flatten out like wet clay. They don't bounce back.
3. The "Swamp" Zone (Elastoplastic): A mix of both. Some parts are springy, some are squished. This is the messy middle ground that is hardest to predict.

4. The "Traffic Light" Diagram

The most practical part of their work is a new diagram (a chart) that acts like a traffic light for engineers.

• If you know how rough your surface is and how the material's hardness changes with size, you can look at this chart.
• It instantly tells you: "Is this contact mostly springy? Mostly squished? Or a mix?"

They found that previous models were often too pessimistic, thinking surfaces were always "squished" (plastic) when they might actually be "springy" (elastic) if you account for the size effect correctly.

5. Why It Matters (According to the Paper)

The authors state that this new framework helps solve a specific puzzle: Why do some surfaces stay bouncy even under high pressure?

They suggest that the "size effect" (getting harder as you zoom in) is a hidden mechanism that keeps contact points from permanently deforming. This is similar to a phenomenon called "asperity persistence," where tiny contact points can hold more weight than expected because they are effectively "work-hardened" by their own small size.

In Summary:
The paper builds a new, faster, and more accurate mathematical "map" for how rough surfaces touch. It corrects old assumptions by showing that materials can get harder as you look closer, allowing some contact points to stay springy rather than permanently squished. They provide a new chart to help engineers quickly guess if a contact will be springy or squished based on the material and surface roughness.

Technical Summary

Technical Summary: Stochastic Model for Elastoplastic Contact of Rough Surfaces Incorporating Scale-Dependent Hardness

Problem Statement
Natural and manufactured surfaces possess inherent roughness that leads to stress concentrations at contact asperities, often inducing plastic deformation. Modeling this elastoplastic contact is complicated by two competing factors: the self-affine nature of surface topography across multiple scales and the "size effect" of plastic deformation. The latter implies that material hardness is not an intrinsic constant but increases with decreasing indentation size (or increasing magnification) due to the accumulation of Geometrically Necessary Dislocations (GNDs). Existing models, such as those based on Persson's diffusion equation or multi-asperity approaches, often assume constant hardness or struggle to efficiently incorporate scale-dependent hardness while maintaining high resolution. Furthermore, previous formulations often rely on solving complex systems of linear equations (e.g., tri-diagonal matrices in diffusion-based approaches), which can be computationally intensive.

Methodology
The authors develop a novel stochastic framework based on the compounded Chapman-Kolmogorov (C-K) equation to solve elastoplastic contact problems involving scale-dependent hardness. The core assumptions and steps include:

1. Markov Process Assumption: The variation of random contact pressure with respect to magnification ($\zeta$) is treated as a Markov process.
2. Integral Equation Formulation: Instead of using a diffusion equation with moving boundary conditions, the authors derive a system of three integral equations describing the evolution of:
   • The probability density function (PDF) of elastic contact pressure, $P_0(p, \zeta)$.
   • The relative plastic contact area, $A^*_{pl}(\zeta)$.
   • The relative non-contact area, $A^*_{nc}(\zeta)$.
3. Scale-Dependent Hardness: The model incorporates a hardness function $H(\zeta) = H_0 \zeta^n$, where $n$ is the hardness exponent. This allows the upper limit of elastic pressure to increase as the observation scale decreases (magnification increases).
4. Backflow Mechanism: A critical theoretical advancement is the relaxation of the "no re-entry" assumption found in constant-hardness models. When hardness increases with scale, surface points previously in plastic contact ($p = H(\zeta')$) can transition back into the elastic contact area or even the non-contact area at higher magnifications. This "backflow" of probability is mathematically captured in the integral equations, distinguishing this model from previous constant-hardness formulations.
5. Numerical Implementation: The magnification axis is discretized using a logarithmic scale. The solution proceeds sequentially through magnification steps, utilizing a transition probability kernel derived from a diffusion approximation with absorbing and non-absorbing boundary conditions.

Key Results
The study presents numerical solutions and convergence tests demonstrating the model's behavior under various conditions:

• Convergence: The model shows excellent convergence with increasing discretization points ($n_\zeta$), yielding results consistent with Persson's theory for linear elastic contact and existing closed-form solutions for constant hardness.
• Contact Status Evolution:
  • Hardness Exponent ($n$): Increasing $n$ (stronger size effect) shifts the contact behavior from fully plastic toward linear elastic. For sufficiently large $n$, the contact remains elastic until full nominal contact is achieved.
  • Hurst Exponent ($H$): Higher $H$ values (smoother surfaces) also promote elastic behavior by reducing the variance of surface height and the rate of pressure redistribution.
  • Symmetry Breaking: Unlike constant-hardness models where elastic pressure PDFs are symmetric about $H_0/2$, the scale-dependent model breaks this symmetry. Probability accumulates near the current hardness limit $H(\zeta)$ due to the backflow mechanism.
• Area-to-Load Relationship: The model predicts a smooth transition from linear elasticity to elastoplastic behavior and finally to full plasticity. The slope of the area-load relationship at low loads ($\kappa$) decreases as the hardness exponent $n$ increases, transitioning from the fully plastic limit to the linear elastic limit.
• Topographic Yield Parameter: The authors propose a new topographic yield parameter, $w(\zeta)$, which quantifies the competition between the redistribution of elastic pressure and the increase in hardness. Unlike previous parameters, this new formulation explicitly incorporates Young's modulus, cut-off wavenumbers, and PSD constants.
• Contact Status Diagram: Using misfit parameters ($M_{le}$ and $M_{fp}$) to quantify the deviation from purely elastic and purely plastic limits, the authors generate a contact status diagram in the $(H, n)$ space. This diagram delineates linear elastic, elastoplastic, and fully plastic regimes. The results suggest that for common rough surfaces ($H \approx 0.8$), even a small size effect ($n > 0$) can lead to predominantly elastic behavior, contrasting with predictions from models assuming constant hardness.

Significance and Claims
The paper claims several significant contributions to the field of contact mechanics:

1. Novel Framework: It pioneers the application of the compounded Chapman-Kolmogorov equation to rough surface contact analysis, offering a stochastic process framework that is distinct from diffusion-based or multi-asperity summation methods.
2. Computational Efficiency: The integral equation approach avoids the need to solve large systems of linear equations (tri-diagonal matrices) required by diffusion-based models, potentially offering a more computationally efficient solution process.
3. Physical Insight into Size Effects: The model provides a mechanism for understanding how scale-dependent hardness alters contact mechanics, specifically the "backflow" of probability from plastic to elastic states, which leads to more elastic behavior at larger scales (lower magnifications) compared to constant-hardness predictions.
4. Practical Utility: The derived contact status diagram offers a tool for the rapid identification of contact regimes (elastic vs. plastic) based on a finite set of mechanical and geometrical properties, aiding in the design and analysis of tribological systems.
5. Broader Applicability: The authors suggest that the integral equations characterizing the evolution of interfacial properties with scale could provide valuable insights for other multidisciplinary fields involving multiscale roughness, such as earthquake mechanics, electrical contact, and contact electrification.

The study concludes that while the model shares a foundational "protuberance-on-protuberance" concept with earlier works (Archard, Gao-Bower, Jackson-Streator), the specific geometric representation (Gaussian rough surfaces with scale-dependent variance) and the inclusion of scale-dependent hardness via the C-K equation lead to distinct governing equations and physical predictions.