超弹性接触的新思路：第三介质法在三维问题中的推广与应用
论文题目：Three-dimensional third medium contact model for hyperelastic contact and pneumatically actuated systems
期刊：Journal of the Mechanics and Physics of Solids
DOI：10.1016/j.jmps.2026.106617
在计算力学中，“接触”一直是一个既经典又让人头疼的问题。
无论是大变形超弹性体、自接触、软体驱动器，还是充气结构中的折叠与闭合，接触处理往往都绕不开几个关键词：接触搜索、主从面、约束施加、非线性收敛、算法鲁棒性。对于做有限元的同学来说，这些内容并不陌生，但真正上手后会发现：理论上清楚，不代表数值实现轻松；模型看似简单，不代表算例就容易收敛。
最近看到一篇很有意思的工作，它提供了一个相当“反常规”的思路：不再显式离散接触界面，而是在潜在接触区域中人为引入一个“第三介质”。这样一来，原本带不等式约束的接触问题，被改写成了一个统一的有限变形边值问题。对做计算力学、尤其是接触算法、软材料、软体机器人方向的同学来说，这篇论文很值得读一读。
一、这篇论文到底在做什么？
作者提出了一个三维第三介质接触（third-medium contact, TMC）框架，用于处理：
超弹性体之间的复杂接触； 大变形下的自接触； 气动驱动系统中的充气/吸气及其诱发的接触问题。
它的核心思想并不复杂：在两个可能发生接触的物体之间，预先填充一层“第三介质”。这层材料在未接触时非常柔软，对主体结构几乎没有影响；但当两个物体逐渐靠近时，第三介质会因为体积压缩而迅速“变硬”，从而像一道屏障一样阻止穿透。
这件事的妙处在于：接触不再依赖传统的接触面搜索与接触约束，而是被吸收到统一的连续体能量泛函中。
二、从能量形式看，它为什么成立？
论文中，主体材料采用超弹性模型，基本应变能密度写成
这里的结构并不陌生： 描述体积变化， 是右 Cauchy-Green 张量， 和 分别对应体积模量和剪切模量。
第三介质的能量则被写成
其中 是一个很小的参数。它的作用非常关键：
在远离接触时，第三介质非常软； 在趋于压扁时， ，于是 ； 对应的能量会急剧增大，从而形成“不可穿透”的屏障效应。
换句话说，这篇论文是把“接触”这件事，转化为了一个通过体积压缩奇异性来实现的障碍机制。
三、真正的难点：网格畸变
第三介质思路虽然漂亮，但它并不是没有代价。
因为一旦接触发生，第三介质区域会被极度压缩，此时最直接的问题不是“穿透”，而是单元畸变。如果畸变失控，计算照样会崩。
因此，作者在这篇论文里真正的重要贡献之一，是提出了一个新的正则化项，用来控制第三介质中的单元质量。
文中新的正则化形式基于增量转角梯度：
其中增量转角向量由
给出。
这个构造很有意思。以往三维中如果直接去处理旋转张量及其梯度，推导和实现都会很重；而这里作者改用“相邻收敛构形之间的增量转角”来近似描述旋转场的空间不均匀性，从而避开了直接求三维旋转张量梯度的复杂性。对做非线性有限元的人来说，这种处理方式是很有启发性的：它不是简单换个符号，而是在数值实现层面主动降低复杂度。
四、接触+气动加载
这篇论文另一个很吸引人的点，是把气动驱动系统也统一纳入了第三介质框架。
作者额外引入了一个气动能量项：
通过这个项，可以在第三介质中产生等效的静水压力应力。于是，充气、吸气、腔体膨胀、腔体闭合，以及由此引发的自接触，都被统一放在同一个能量框架里处理了。
这对软体机器人、柔性机构、可充气结构模拟来说，意义是很直接的：加载与接触不再是两个割裂的问题，而是同属一个统一的连续体模型。
五、数值算例说明了什么？
论文安排了多组很有代表性的算例。
首先是三维盒体内部自接触。作者比较了不同正则化项与不同参数的影响，结果表明：新提出的正则化项在相同参数下可以得到更小的接触间隙，说明它更能抑制单元重叠、提升接触精度。更重要的是，在计算效率上，新正则化项也更有优势：文中给出的时间开销对比中，旧形式约为 4449.11 s，而新形式约为 2469.23 s，改进是比较明显的。
接着是旋转诱导的三维自接触。这个算例更接近复杂自由表面的接触情形。作者发现，旧正则化在旋转角超过 478° 后失效，而新正则化可以支撑到 529°，这说明新方法在大转动、自接触和严重畸变条件下更稳定。
随后是经典的压头压入问题（punch problem）。这类问题的价值在于，它不是封闭腔体内部自接触，而是两个独立实体之间的接触。论文结果表明，该框架同样可以稳定处理这类三维接触问题。
再往后，是带腔体的气动驱动盒体。作者分别模拟了吸气和充气两类工况，并展示了在吸力较大时，结构会向内收缩并发生自接触。这一结果很重要，因为它说明第三介质方法不仅能施加气动作用，还能自动处理气动载荷诱发的接触行为。
最后一个算例是软体机器人气动驱动器。随着内压增加，各个气室逐渐膨胀、闭合并彼此接触，驱动器整体产生弯曲变形。作者还将结果与 ABAQUS 进行了对比，二者在变形模式上具有较好一致性。这说明该方法并不是“只会做概念演示”，而是具备面向工程应用的潜力。
六、这篇论文对我们有什么启发？
我觉得，这篇工作的价值不只是“又做了一个接触算法”，而是它提供了一种很值得重视的建模观念：
把原本依赖界面搜索与约束施加的问题，转化为区域型、能量型、统一连续体的问题。
对于传统接触有限元来说，困难往往集中在界面层面；而第三介质方法则相当于绕开了这个正面战场，从体积分布能量的角度重新组织问题。这种思路对于复杂接触、自接触、多腔体驱动、软材料大变形等场景，确实很有吸引力。
当然，论文也没有回避它当前的限制。比如文中明确提到：现阶段正则化项仍要求二次单元，这会影响计算效率；未来还需要进一步发展一阶格式，并扩展到摩擦接触。但从研究路径来看，这已经是一个非常清晰、非常有潜力的方向了。
如果你正在做接触算法、超弹性有限元、软体机器人，或者对“如何把复杂非线性问题改写成更统一的变分框架”感兴趣，那么这篇文章值得认真读一遍。
有时候，真正好的工作，不只是把已有方法做得更复杂，而是把原本复杂的问题，换一种更统一的方式重新表达。
这篇关于 third-medium contact 的论文，我觉得就属于这种类型。
文献信息
Xu, B.-B., Xue, T., & Wriggers, P. Three-dimensional third medium contact model for hyperelastic contact and pneumatically actuated systems. Journal of the Mechanics and Physics of Solids, 213 (2026), 106617.
DOI: 10.1016/j.jmps.2026.106617
