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

本文提出了一种基于复合查普曼-柯尔莫哥洛夫方程的新型随机模型，通过引入尺度相关硬度来分析粗糙表面的弹塑性接触，从而预测跨尺度的接触状态演变，并为涉及多尺度粗糙度的多学科领域提供新的见解。

要点

想象一下两个粗糙的表面，比如两块砂纸或轮胎与路面接触。尽管从远处看它们很平整，但如果你放大观察，会发现它们其实布满了微小的“山峰”和“山谷”。当你将它们压在一起时，只有这些“山峰”（称为微凸体，asperities）的最顶端才会发生接触。

这篇论文旨在理解在这些微小的接触点上究竟发生了什么，特别是当材料足够柔软、会发生挤压变形（塑性变形）且具有一定的回弹特性（弹性）时。

以下是利用简单类比对他们发现进行的解析：

1. 问题所在：“缩放”谜题

大多数旧有的表面接触模型都假设材料的硬度不随观察尺度的变化而改变。但在现实中，如果你观察一个极小的点，材料表现出的硬度通常比观察大范围区域时要高。这就是所谓的尺寸效应（size effect）。

你可以把这想象成人群。如果你观察整个体育场，移动起来很容易；但如果你放大到只看三个并肩而立的人，挤过去就会变得非常困难。当你观察得越细微，这个“人群”（材料）感觉就越硬。

作者希望建立一个能够解释这种硬度随缩放比例变化的模型的模型。

2. 解决方案：一张“概率流”地图

作者并没有尝试模拟每一个微小的山峰（那需要超级计算机运行很久），而是使用了一个巧妙的数学技巧——查普曼-科尔莫戈罗夫方程（Chapman-Kolmogorov equation）。

类比： 想象一条向河流下游流动的河水。

• 水流： 代表接触点的“压力”。
• 河床： 代表表面的粗糙度。
• 河岸： 代表材料的极限。如果水位过高，水就会漫过河岸（材料发生屈服或挤压变形）。

在过去，科学家认为水只能单向流动：从平静的水池流入湍急的激流，一旦触及河岸，就会停留在那里。他们假设一旦某个点被挤压了（发生塑性变形），无论你如何放大观察，它都会保持被挤压的状态。

新的发现： 作者发现，当硬度随尺度变化时，水流实际上可以向后流动。

• 如果一个点在低缩放级别下被挤压了，但随着缩放倍数的增加，材料变得更硬了，那么那个点实际上可能会“恢复原状”，重新变为弹性状态。
• 他们创建了一张地图，用于追踪一个点处于“接触”、“挤压”或“未接触”状态的概率是如何随缩放过程变化的。

3. 三种接触区域

利用这种新的数学方法，他们根据材料属性和表面的“粗糙度”，识别出了两种表面相互作用的三种不同状态：

1. “弹簧”区（线性弹性）： 表面发生接触，但表现得像坚硬的弹簧。它们会轻微挤压并回弹。这发生在材料在微观尺度下非常硬的情况下。
2. “泥泞”区（完全塑性）： 表面过于柔软或压力过高，导致“山峰”像湿粘土一样被压平。它们不会回弹。
3. “沼泽”区（弹塑性）： 两者的混合体。有些部分是弹性的，有些部分是挤压变形的。这是最复杂、最难以预测的中间地带。

4. “红绿灯”图表

这项工作中最具实用价值的部分是一个新的图表（一种图表），它充当了工程师的“红绿灯”。

• 如果你知道表面的粗糙程度以及材料硬度随尺寸的变化规律，你就可以查看这张表。
• 它能立刻告诉你：“这种接触主要是弹性的？还是主要是挤压变形的？或者是两者的混合？”

他们发现，以往的模型往往过于悲观，认为表面总是处于“挤压”状态（塑性），而如果正确考虑尺寸效应，这些表面实际上可能是“弹簧”状态（弹性）。

5. 为什么这很重要（根据论文所述）

作者指出，这个新框架有助于解决一个特定的谜题：为什么有些表面在承受高压时仍能保持弹性？

他们认为，“尺寸效应”（即随着缩放变小而变硬）是一种隐藏机制，它防止了接触点发生永久性的变形。这类似于一种被称为“微凸体持久性（asperity persistence）”的现象，即微小的接触点由于其自身的微小尺寸，实际上已经产生了“加工硬化”，从而能承载比预期更多的重量。

总结：
这篇论文构建了一个关于粗糙表面如何接触的、更快速且更准确的数学“地图”。通过展示材料在观察尺度变小时会变硬，该研究修正了旧有的假设，使得一些接触点能够保持弹性而非永久变形。他们提供了一张新图表，帮助工程师根据材料和表面粗糙度，快速判断接触属于弹性还是塑性变形。

技术摘要

技术摘要：考虑尺度相关硬度的粗糙表面弹塑性接触随机模型

问题陈述
天然及人工制造的表面具有固有的粗糙度，这会导致接触微凸体处产生应力集中，从而诱发塑性变形。对这种弹塑性接触进行建模非常复杂，原因在于两个相互竞争的因素：表面形貌在多个尺度上的自仿射特性，以及塑性变形的“尺寸效应”。后者意味着材料硬度并非一个内在常数，而是随着观测尺度减小（或放大倍数增加）而增加，这是由于几何必要位错（GNDs）的积累所致。现有的模型（例如基于 Persson 扩散方程或多微凸体方法的研究）通常假设硬度为常数，或者难以在保持高分辨率的同时高效地纳入尺度相关的硬度。此外，以往的公式往往依赖于求解复杂的线性方程组（例如基于扩散方法的三对角矩阵），这在计算上可能非常耗时。

方法论
作者开发了一种基于复合 Chapman-Kolmogorov (C-K) 方程的新型随机框架，用于求解涉及尺度相关硬度的弹塑性接触问题。其核心假设和步骤包括：

1. 马尔可夫过程假设： 将随机接触压力随放大倍数 ($\zeta$) 的变化视为一个马尔可夫过程。
2. 积分方程表述： 作者没有使用带有移动边界条件的扩散方程，而是推导了一个由三个积分方程组成的系统，用于描述以下各项的演化：
   • 弹性接触压力的概率密度函数 (PDF) $P_0(p, \zeta)$。
   • 相对塑性接触面积 $A^*_{pl}(\zeta)$。
   • 相对非接触面积 $A^*_{nc}(\zeta)$。
3. 尺度相关硬度： 模型引入了硬度函数 $H(\zeta) = H_0 \zeta^n$，其中 $n$ 是硬度指数。这使得观测尺度减小时（放大倍数增加），弹性压力的上限会随之提高。
4. 回流机制（Backflow Mechanism）： 一个关键的理论进展是放宽了常硬度模型中的“无重入”假设。当硬度随尺度增加时，在较高放大倍数下，先前处于塑性接触状态（$p = H(\zeta')$）的表面点可能会转回到弹性接触区域甚至非接触区域。这种概率的“回流”通过积分方程在数学上得以捕捉，使该模型区别于以往的常硬度公式。
5. 数值实现： 使用对数尺度对放大倍数轴进行离散化。解的过程通过放大倍数步骤顺序进行，利用了由具有吸收和非吸收边界条件的扩散近似所导出的转移概率核。

关键结果
研究通过数值解和收敛性测试，展示了模型在各种条件下的行为：

• 收敛性： 模型表现出极佳的收敛性，随着离散点数 ($n_\zeta$) 的增加，其结果与 Persson 的线性弹性接触理论以及现有关于常硬度的闭式解保持一致。
• 接触状态演化：
  • 硬度指数 ($n$)： 增加 $n$（增强尺寸效应）会使接触行为从完全塑性向线性弹性转变。对于足够大的 $n$，接触在达到完全标称接触之前将保持弹性。
  • 赫斯特指数 ($H$)： 更高的 $H$ 值（更光滑的表面）也通过减少表面高度的方差和压力重新分布的速率来促进弹性行为。
  • 对称性破缺： 与弹性压力 PDF 关于 $H_0/2$ 对称的常硬度模型不同，尺度相关模型打破了这种对称性。由于回流机制的存在，概率会在当前硬度极限 $H(\zeta)$ 附近累积。
• 面积-载荷关系： 模型预测了从线性弹性到弹塑性，再到完全塑性的平滑过渡。在低载荷下，面积-载荷关系的斜率 ($\kappa$) 随着硬度指数 $n$ 的增加而减小，实现了从完全塑性极限向线性弹性极限的过渡。
• 形貌屈服参数： 作者提出了一个新的形貌屈力参数 $w(\zeta)$，用于量化弹性压力重新分布与硬度增加之间的竞争。与以往参数不同，该新公式明确纳入了杨氏模量、截止波数和功率谱密度 (PSD) 常数。
• 接触状态图： 利用失配参数（$M_{le}$ 和 $M_{fp}$）来量化对纯弹性和纯塑性极限的偏离，作者在 $(H, n)$ 空间中生成了接触状态图。该图划定了线性弹性、弹塑性和完全塑性区域。结果表明，对于常见的粗糙表面（$H \approx 0.8$），即使是很小的尺寸效应（$n > 0$）也会导致以弹性行为为主，这与假设常硬度的模型预测形成对比。

意义与主张
本文声称在接触力学领域做出了以下显著贡献：

1. 新颖框架： 它开创了将复合 Chapman-Kolmogorov 方程应用于粗糙表面接触分析的先河，提供了一个有别于基于扩散或多微凸体求和方法的随机过程框架。
2. 计算效率： 积分方程方法避免了求解基于扩散模型的庞大线性方程组（三对角矩阵）的需求，从而提供了更高效的求解过程。
3. 对尺寸效应的物理洞察： 模型提供了一种理解尺度相关硬度如何改变接触力学的机制，特别是从塑性向弹性状态转换的概率“回流”，这导致在较大尺度（较低放大倍数）下表现出比常硬度预测更多的弹性行为。
4. 实用价值： 推导出的接触状态图提供了一个工具，能够基于有限的机械和几何属性快速识别接触状态（弹性 vs. 塑性），有助于摩擦学系统的设计与分析。
5. 广泛适用性： 作者指出，表征界面性质随尺度演化的积分方程可以为涉及多尺度粗糙度的其他多学科领域（如地震力学、电接触和接触电化）提供有价值的见解。

研究结论认为，尽管该模型与早期工作（Archard, Gao-Bower, Jackson-Streator）共享基本的“突起对突起”概念，但其特定的几何表示（具有尺度相关方差的高斯粗糙表面）以及通过 C-K 方程引入的尺度相关硬度，导致了截然不同的控制方程和物理预测。