Strength-degradation phase-field regularization of cohesive fracture: the antiplane case

想象一下，你正试图预测一块玻璃、一团粘土甚至一根金属梁会如何破碎。长期以来，科学家们在看待这个问题时有两种截然不同的观点。一方面是“脆性”视角：材料会突然断裂，就像一根干枯的树枝，此时唯一重要的因素是创造一个新的裂纹表面需要多少能量。另一方面是“延性”视角：材料在破碎前会像太妃糖一样挤压和拉伸，其关键在于材料在压力下发生屈服或变形的程度。

棘手之处在于，现实世界的材料往往兼具两者的特性。一块金属在最终断裂（断裂）之前可能会发生一定的拉伸（塑性）。几十年来，计算机模型一直难以处理这种混合情况。流行的“相场”（phase-field）模型——类似于一种数字模拟，让裂纹能够平滑地生长，而不是在像素之间跳跃——存在一个重大缺陷：它将材料的强度直接与计算机屏幕上裂纹看起来有多“模糊”联系在一起。如果你想模拟一种强韧的材料，你就必须让计算机的“模糊度”变得极小，这需要超级计算机。如果你为了节省时间而增大模糊度，材料就会神奇地变弱。这变成了一场令人沮挫的“三选二”游戏：“你可以拥有强度，或者拥有速度，或者拥有真实的裂纹形状，但你不能同时拥有这三者。”

这篇论文介绍了一种全新的建模方式，终于解开了这个结。作者布莱斯·布尔丁（Blaise Bourdin）和科拉多·毛里尼（Corrado Maurini）提出了一种模型，其中材料的强度是一个独立的、单独的成分，就像它的刚度或韧性一样。他们在一种简化的设置中测试了这个想法，称为“反平面剪切”（想象一下将一副扑克牌向侧面滑动），并发现他们的新模型表现得非常出色。它可以模拟从金属拉伸和挤压（塑性）到干净、突然断裂（脆性断裂）的一切过程，而无需改变游戏规则。最重要的是，他们证明了计算机模拟中的“模糊度”只是一个帮助数学运算工作的数值工具，而不是改变材料强度的物理属性。这意味着科学家终于可以在标准计算机上模拟复杂的破碎行为，而不会因为调整网格大小而改变结果。

新的“强度优先”方法

这篇论文的核心是一种用于模拟物体破碎的新数学配方。作者称之为“强度退化相场正则化”（strength-degradation phase-field regularization）。这个词很拗口，让我们用一个简单的类比来拆解它。

想象一种材料就像一群手拉手的人。在旧模型中，当人群开始感到压力时，人们的手臂（刚度）会变得越来越弱，直到松开。问题在于，他们松开的时机完全取决于你观察他们的紧密程度（“正则化长度”）。如果你观察得很仔细，他们会坚持得更久；如果你离得远，他们就会较早松手。

在这个新模型中，人们并不会削弱自己的手臂。相反，他们有一个最大握力。只要拉力低于这个强度，他们就会保持稳固。一旦拉力达到这个极限，握力就开始退化，他们最终会松手。至关重要的一点是，这个“最大握力”是人们固有的属性，而不是你观察方式带来的副作用。这使得该模型能够处理既强韧又不易拉伸的材料，或者在断裂前会大幅拉伸的材料，且全部在同一个框架内实现。

破碎的三个阶段

作者通过模拟实验观察了他们的模型在不同条件下的表现。他们发现，根据物体的大小和材料属性，他们的单一模型可以自然地过渡到三个不同的失效阶段：

1. “挤压”阶段（微观尺度屈服）： 当一个小区域开始形成裂纹时，裂纹尖端的材料并不会直接断裂。它会拉伸和变形，产生一个“过程区”，在该区域内材料发生屈服。该模型完美地捕捉到了这一点，展示了一个符合经典塑性理论的变形区。
2. “内聚”阶段： 随着裂纹张开，材料并不会立即断裂。相反，它会保持一种“内聚”力，就像正在缓慢剥离的胶带一样。拉开它所需的力随着间隙变宽而减小。这就是“内聚裂纹”行为，对于理解现实材料如何失效至关重要。
3. “脆性”阶段： 最后，如果间隙变得足够大，材料就会彻底放弃并断裂。裂纹变成了一个干净、锐利的断口，就像最简单的脆性断裂模型那样。

这项工作的精妙之处在于，它不需要在描述这三个阶段时切换不同的数学公式。这是一个连续的过程。模型仅仅根据物理规律，而非人类程序员的选择，从一个状态演变为另一个状态。

“冲浪”实验与“塑性尾迹”

论文中最引人入胜的部分之一是作者称之为“冲浪”的模拟。想象一个已经存在的裂纹，你正以恒定的速度拖动它穿过材料，就像冲浪者乘着浪行进一样。作者以两种方式运行了这次模拟：一次是材料的变形是可逆的（像橡皮筋），另一次是不可逆的（像会保持变形形状的粘土）。

在可逆的情况下，破坏材料所需的能量与理论预测完全吻切。但在不可逆的情况下，有趣的事情发生了。随着裂纹的移动，它留下了一个“塑性尾迹”——一段永久变形的材料痕迹，就像船只留下的航迹。因为材料必须被挤压才能让裂纹通过，所以破坏材料所需的能量比简单脆性裂纹理论所暗示的要高。作者测量了这种“有效韧性”，发现它比材料的基础韧性高出约 30%。这解释了为什么现实中的某些材料看起来比其基本公式预测的更具韧性：裂纹尖端前方的塑性变形正在做额外的功。

通过 V 型缺口架起桥梁

为了证明他们的模型适用于复杂形状，作者模拟了一个“V 型缺口”——即材料中一个尖锐的切口，类似于 Pac-Man（吃豆人）形状的角。他们对这个缺口施加力，并观察发生了什么。

他们发现，随着载荷增加，材料经历了相同的三个阶段：首先，缺口处形成了一个小的塑性区（“挤压”阶段）；然后，一个内聚裂纹开始生长（“胶带”阶段）；最后，一个脆性裂纹突然产生并贯穿整个材料。

这里最令人印象深刻的结果是，启动裂纹所需的力遵循一条“通用法则”。无论缺口是非常尖锐（如裂纹）还是非常钝（如直边），模型的预测都汇聚在同一条曲线上。这条曲线将材料的强度（破坏它需要多大的拉力）与韧性（破坏它需要多少能量）平滑地连接在一起。这表明该模型可以预测任何形状的失效，而无需猜测使用哪种公式。

这为什么重要

作者谨慎地指出，这些结果来自模拟和数学证明，而非对真实材料的物理实验。然而，模拟结果非常稳健，并且完美符合已知的理论极限。他们已经证明，通过改变模型退化强度而非刚度的方式，他们创造了一个既灵活又准确的工具。以前，如果你想模拟一种强韧但脆性的材料，你会陷入困境：你要么使用计算量巨大（需要极其精细的网格，会导致电脑崩溃）的模型，要么使用给出错误答案的模型。而这种新方法允许“正则化长度”（计算机的模糊度）仅仅作为一个数值工具。你可以将其调小以获得锐利的裂纹，或者调大以节省计算能力，而材料的强度和韧性将保持完全不变。

简而言之，这篇论文为破碎提供了一种统一的语言。它将塑性（挤压）、内聚断裂（剥离）和脆性断裂（断裂）的世界整合进了一个单一且一致的框架中。它表明我们不需要在这些不同的失效类型之间做出选择；它们只是材料破碎这一故事中不同的章节。对于工程师和科学家来说，这意味着我们可能很快就能以以往难以企及的细节和精度，来模拟桥梁、飞机或生物组织等复杂结构的失效过程。