A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

当像玻璃或硬化钢这样的脆性材料断裂时，它们并不仅仅是简单地折断成两半。相反，一道裂纹会在固体中疾速蔓延，通常会分支成复杂的破碎网络，将物体粉碎成许多碎片。预测这些裂纹究竟会沿着什么路径行进以及移动速度有多快，是工程学中最困难的挑战之一。如果用于模拟这种破坏的模型存在缺陷，它可能会暗示一座桥梁会向一个安全的方向失效，而实际上它却发生了灾难性的坍塌。几十年来，科学家们一直试图构建能够模拟材料原子结构的计算机模型，将材料视为由微小点通过弹簧连接而成的网格。虽然这些模型功能强大，但它们存在一个奇特的缺陷：裂纹行进的路径往往取决于计算机网格是如何旋转的，而不是取决于材料本身的物理特性。这就像是现实中的一块玻璃，仅仅因为工程师转动了计算机屏幕，其内部的裂纹就会突然改变方向一样。

为了解决这个问题，马蒂亚斯·布劳恩（Matías Braun），西班牙内布里哈大学的一位研究人员，开发了一种模拟动态断裂的新方法。他的方法将传统的点阵网格与“相场”（phase-field）法相结合，这种技术将裂纹视为一个在材料中扩散的平滑、模糊的损伤区域，而不是一个在点与点之间跳跃的尖锐线条。通过融合这两种思路，布劳恩创建了一个允许裂纹向任何方向行进的模型，无论底层的计算机网格如何定向。他利用三种不同的场景测试了这种新方法：一块受到高速撞击的钢板、一块被拉伸直至断裂的有缺口平板，以及一个从内部向外爆炸的厚金属圆柱体。在每种情况下，新模型产生的结果都与现实世界的实验相吻合，并消除了由计算机网格引起的人为偏差，为材料如何真正失效提供了更清晰、更可靠的视角。

布劳恩所解决的核心问题在于计算机如何处理原子的不可见世界。在他早期的工作中，他使用了一种离散模型，其中材料由一个三角形晶格（一种由点和线组成的图案）来表示。当材料受到应力时，点与点之间的连接会断开，从而模拟裂纹。然而，由于这些点的排列是固定的模式，裂纹被迫沿着网格的线条行进。如果网格旋转，裂纹就会改变路径，尽管材料本身并未发生变化。这是一个显著的问题，因为在现实中，裂纹应该遵循由材料属性决定的最小阻力路径，而不是由计算机网格的任意角度决定的。

布劳恩的解决方案是引入相场正则化。他不再强迫特定的点间连接断裂，而是为每个点分配一个代表其受损程度的数值。该数值范围从零（意味着材料完全完好）到一（意味着完全破碎）。裂纹变成了一个数值逐渐增加的区域，创造了一个平滑的过渡而非锯齿状的线条。这使得裂纹可以曲线行进并进行分支，有效地忽略了刚性网格的结构。至关重要的是，布劳恩使用材料的两个基本属性对该系统进行了校准：抗拉强度（即开始产生裂纹所需的力）和断裂能（即维持裂纹生长所需的能量）。通过将模拟锚定在这些真实的物理常数上，该模型确保了结果取决于材料本身，而非数学计算。

为了证明这种新方法有效，布劳恩进行了三项严格测试。第一项是卡尔托夫-温克勒（Kalthoff–Winkler）测试，这是一项经典的实验，其中一个带有预切缺口的钢板受到弹丸的高速撞击。在这种情景下，裂纹会从缺口处以特定角度射出。以往的模型往往会根据网格的旋转情况预测错误的路径角度。在布劳恩的模拟中，裂纹始终以大约70度的角度行进，与实验观察完美匹配。无论网格是与平板对齐还是旋转了15或30度，这一结果都保持不变。裂纹速度也保持稳定，约为材料声速的70%，且不受网格定向的影响。

第二项测试涉及一个在张力作用下的单边缺口平板，这种装置用于研究裂纹在持续拉力下如何缓慢生长。这里的目标是观察裂纹是否会从缺口沿直线运动到对侧边缘。传统的网格模型往往会迫使裂纹根据网格的对齐方式产生锯齿状运动或发生偏离。然而，布劳恩的模型对于每种网格取向都能产生直线裂纹。打破平板所需的力也是一致的，无论网格如何旋转，其变化量均小于2%。这证明了该模型可以在不受计算机网格人为影响的情况下，准确预测材料的路径和强度。

最后也是最复杂的测试是一个受到内部爆炸的厚圆柱体的破碎过程。在这种情景下，没有单一的裂纹受缺口引导；相反，许多裂纹必须自发产生并向外奔袭，将圆柱体粉碎成数十块碎片。这是一个混沌的过程，网格的对称性很容易决定碎片的分布模式。当布劳恩运行模拟并将网格中心对准圆柱轴线时，裂纹呈现出对称模式。然而，当他稍微移动网格使其不再完美居中时，裂纹便呈现出不规则分布，且碎片的数量也发生了变化。这证实了该模型对物理设置具有敏感性，而非仅仅是对网格对称性的反应。模拟预测会有12到19个碎片，这个数字处于其他先进方法的范围内，但略低一些，这可能是因为模型中的平滑损伤区阻止了过于靠近的裂纹形成。

在所有这些测试过程中，研究人员都仔细追踪了系统的能量。他们将材料断裂时损失的能量与产生新裂纹表面所需的能量进行了对比。在某些情况下，损伤模型计算出的能量与实际释放的断裂连接能量之间存在微小的差异。这是混合方法的已知局限性，即驱动裂纹发展的力与断裂键的能量计算方式不同。然而，机械能平衡的准确度保持在0.5%以内，证实了模拟在物理上是合理的。结果表明，虽然该模型并不完美，但它代表了理解动态断裂的一个重要进步。

这项工作的意义不仅在于提升计算机图形效果。通过消除对网格定向的依赖，工程师可以信任他们的碰撞测试、结构失效或材料测试模拟是在反映材料真实的物理特性。该模型成功架起了原子离散性质与固体连续行为之间的桥梁。它表明，通过正确的数学工具，模拟破碎这一混沌且剧烈的过程而不受计算机网格支配是完全可能的。研究结论指出，这种方法为预测材料在极端条件下的失效提供了一种稳健的方式，为设计更安全、更具韧性的结构提供了更清晰的路径。