Accurate simulation of delamination with a resin-rich layer-dependent penalty stiffness based on structural cohesive elements

想象一下，你正在用成千上万块微小的、色彩斑斓的乐高积木搭建一座房子，并将它们粘在一起组成一面超级坚固的墙。在现实世界中，工程师们使用碳纤维和塑料树脂层来建造类似的结构，用于制造飞机机翼或火箭部件。这些材料极其坚固，但它们有一个秘密弱点：层与层之间有时会发生剥离，就像书页被撕开一样。这种剥离被称为“分层”（delamination），如果它在不经意间发生，整个结构可能会失效。为了防止这种情况，科学家们使用计算机程序来模拟这些层是如何分离的。然而，这些程序经常面临一个棘手的难题：为了让数学计算准确，它们需要极其精细（使用数百万个微小的数字积木），这使得计算机运行模拟的时间变得非常漫长。如果数字积木太大，计算机就会产生混乱，预测层会以不可能的方式挤压在一起，从而导致错误的答案。

这篇论文通过发明一种更聪明的方法来告诉计算机这些层在剥离时应该如何表现，解决了这一困惑。作者 Xiaopeng Ai、Christos Kassapoglou 和 Boyang Chen 提出了一种全新的数学“规则书”，用于描述层与层之间“胶水”的行为。他们不再根据整个墙壁的厚度来猜测胶水的硬度，而是观察自然存在于纤维层之间的那层微小且杂乱的额外树脂层。通过一种基于梁弯曲原理的巧妙技巧，他们弄清楚了压力是如何在这一薄层胶水中传播的。他们的新方法允许计算机使用更大型、更粗糙的数字积木，而不会导致数学错误。这意味着工程师可以更快、更准确地模拟复杂的剥离问题，确保使用这些材料制造的飞机和火箭安全可靠。

问题所在：这种“软塌塌”的胶水错误

把复合材料想象成一个由许多层面包（纤维）和每片面包之间的一点点额外黄油（树脂）组成的巨大三明治。当三明治开始撕裂时，黄油层是首先拉伸并断裂的部分。在过去，计算机模拟将这个黄油层视为一个简单的、均匀的弹簧。问题在于，这个“弹簧”往往会根据计算机如何将三明治切割成数字碎片，表现得过硬或过软。

当计算机使用大的数字切片时，它会犯一个特定的错误：它预测撕裂前方的层正受到巨大的挤压力量，就像车祸现场一样，而实际上，它们应该只是在轻轻地剥离。这种“软塌塌”的误差意味着，为了得到正确的答案，工程师必须使用极小的数字切片，这导致模拟需要运行数天时间。这就像是用只有微小、锯齿状的步阶来画一条平滑的曲线；你需要数百万个步阶才能让它看起来正确。

新思路：倾听应力

作者意识到，旧的计算胶水“刚度”的方法过于简单。他们知道，压力（内部力量）并不会均匀地穿过三明治的厚度；它会随着你从顶层移动到底层的变化而变化。为了解决这个问题，他们不仅观察了胶水，还观察了梁弯曲的物理现象。

想象一下弯曲一本厚书。书的外侧页面会被拉伸，内侧页面会被挤压，而中间的页面则保持中性。作者利用同样的逻辑，绘制出了压力随厚度变化的精确图谱。他们推导出了一个全新的公式，充当了一个“智能翻译官”。它能将复杂的、变化的压力模式转化为计算机可以使用的单一、准确的刚度数值。

至关重要的是，即使当计算机将许多同类型的层组合成一个大的“超层”（等效单层模型）时，这个新公式依然有效。以往的方法要求计算机对每一层进行单独建模，这既缓慢又低效。新方法允许计算机将层进行分组，而不会丢失精度，使模拟运行速度大大提升。

证明：测试新的规则书

为了验证他们的想法是否真的奏效，团队使用工程师用来检查数学准确性的标准“基准”问题进行了系列测试。这些测试包括：

1. 双悬臂梁 (DCB)： 一种像打开书本一样将两层拉开的测试。
2. 端部缺口弯曲 (ENF)： 一种让层与层之间相互滑动进行的测试。
3. 混合模式弯曲 (MMB)： 一种结合了拉伸和滑动的测试。
4. 增强型双悬臂梁 (R-DCB)： 一种更复杂的测试，在外部添加了额外的层，模拟真实的加固结构。

在每一次测试中，他们都将这种新的“智能”刚度与旧的“愚笨”刚度进行了对比。结果非常明确：

• 网格收敛性： 当他们使用粗糙（大块）的数字网格时，旧方法失败了。它的答案取决于碎片的大小，并且预测层正在过度挤压。然而，新方法无论数字碎片是大是小，都能给出同样准确的答案。这就像拥有一张无论放大还是缩小看起来都同样完美的地图。
• 准确性： 在 DCB 测试中，新方法显著降低了载荷-位移曲线（即剥离层所需的力）中的误差。对于 7.5 mm 的网格，旧方法的误差超过 50%，而新方法则保持在极低的水平。
• 现实世界对比： 在与真实物理实验进行对比的复杂 R-DCB 测试中，新方法对裂纹路径和破坏所需力量的预测比旧方法要好得多。具体而言，它预测的破坏力第二个峰值比旧方法低 14%，这与现实世界的数据吻合得更加紧密。

这意味着什么

这篇论文并不声称已经解决了复合材料领域的所有问题。它专门针对改进“结构内聚单元”的“惩罚刚度”，这些单元是用于模拟剥离的数字工具。作者通过模拟以及与实验的对比，展示了他们的新公式优于传统公式。

通过使用这种新方法，工程师可以在较粗的网格下模拟复杂 3D 结构中的分层现象。这意味着他们可以在极短的时间内获得准确的结果，而不需要像以前那样耗费大量时间。这就像是从一个缓慢、像素化的视频游戏升级到了流畅的高清游戏，而且不需要超级计算机。作者建议，这种方法未来可以扩展到预测其他类型的损伤，例如层内部的裂纹，但目前，他们已经成功解决了这个“软塌塌的胶水”问题，使复合材料的数字测试变得更快、更可靠。