A new approach for the determination of through-thickness and free-edge stresses in composite laminates based on structural elements

想象一下，你正在用成千上万块色彩缤纷的小乐高积木搭建一座摩天大楼。在现实世界中，工程师们并不使用塑料，而是使用“复合材料层合板”——它们就像是用胶水（树脂）浸泡过的织物层叠而成的超强“三明治”。这些材料是轻型飞机和快速汽车背后的秘密武器，因为它们强度极高，但重量却非常轻。然而，这里有一个棘手的难题：当我们弯曲或拉伸这些“三向盘”时，各层之间会想要相互滑动或剥离，尤其是在边缘附近。这被称为“层间剥离”（delamination），它是结构中的“无声杀手”。

为了预测机翼是否能承受住压力，工程师们使用计算机模拟。对于薄结构，他们通常使用一种叫做“壳单元”（shell elements）的捷径。你可以把壳单元想象成计算机程序中的一张平整纸片。它非常擅长计算纸张表面的拉伸或弯曲，但对于研究纸张厚度内部发生的情况却表现得很糟糕。这就像是试图通过只观察三明治顶层面包来理解里面的奶酪和肉一样；你会完全错过中间的填充物。因为这些标准的计算机模型无法“看到”内部的应力（即推动各层分离的力量），它们往往会错过结构失效前的预警信号。本文的研究正是为了填补这一盲点，让工程师能够在不构建每一个胶水颗粒那般庞大、缓慢且复杂的 3D 模型的情况下，看清这些隐藏在内部的力量。

---

本文的核心思想：看见隐形的层级

作者 Xiaopeng Ai、Christos Kassapoglou 和 Boyang Chen 构思了一种巧妙的技巧，可以窥探这些复合材料“三明治”的内部。他们基于一种方法进行开发，该方法将层合板的每一层视为一个“结构单元”（一种高级版的“纸片”），并用“内聚单元”（cohesive elements）将它们连接起来。你可以把这些内聚单元想象成层与层之间看不见的、具有弹性的“胶条”。

在过去，这些“胶条”一直是个谜。工程师们必须猜测它们的刚度应该是多少。如果猜错了，计算机要么会认为层与层之间是用超强环氧树脂粘在一起的（太硬），要么会认为是弱胶带粘在一起的（太软），从而导致关于结构能承受多少应力的错误答案。本文指出，过去猜测这种刚度的方法是有缺陷的，因为它们假设胶水没有厚度，但在现实世界中，层与层之间总会存在一小层富树脂空间。

“富树脂”的秘密配方

该团队的主要发现是一个用于精确计算这些“胶条”应该有多硬的新公式。他们不再靠猜测，而是意识到刚度取决于富树脂层（即织物层之间那个微小的、胶水含量高的区域）的实际属性。

想象一下你在拉伸一叠便利贴。你感受到的力量取决于粘性部分有多厚以及胶水的强度有多大。作者推导出了一个数学规则，该规则指出：“胶条的刚度简单来说就是胶水的强度除以那个微小树脂层的厚度。”

他们在不同类型的复合材料“三明治”上测试了这个想法：

• 十字铺层（Cross-ply laminates）： 像棋盘一样堆叠的层（0° 和 90°）。
• 斜铺层（Angle-ply laminates）： 像锯齿一样呈角度堆叠的层（+45° 和 -45°）。
• 准各向同性层合板（Quasi-isotropic laminates）： 在各个方向表现一致的复杂堆叠。

在这些模拟中，他们对层合板的两端施加了均匀的拉伸，就像拉伸橡皮筋一样。然后，他们将这种新方法与两样东西进行了对比：

1. 解析解（Analytical solutions）： 这是由其他科学家（如 Pagano 和 Mittelstedt）推导出的完美数学答案，作为“金标准”。
2. 实体单元模型（Solid-element models）： 这是重型的、缓慢的计算机模型，它将“三明治”分解成数百万个微小的 3D 块来观察每一个细节。

研究结果

结果令人印象深刻。当作者使用这种新的“富树脂”刚度公式时，他们简单的壳模型几乎完美地预测了内部应力。

• 准确性： 他们的结果与“金标准”数学解以及重型 3D 块模型高度吻合。他们能够准确预测层间法向应力（试图将各层剥离的力量，记作 $\sigma_z$）以及剪切应力（试图让各层相互滑动的力量，记作 $\tau_{xz}$ 和 $\tau_{yz}$）。
• 自由边缘效应（Free-Edge Effects）： 最难预测的事情之一就是材料边缘（即各层暴露在空气中的地方）会发生什么。这也是“应力奇异性”（极端的力峰值）发生的地方。新方法能够准确捕捉到这些峰值，而旧方法往往会漏掉它们，或者需要极其精细且缓慢的网格才能观测到。
• 鲁棒性（稳定性）： 团队还检查了该方法对微小变化的敏感程度。他们问道：“如果树脂层稍微厚一点或薄一点会怎样？如果胶水稍微强一点或弱一点会怎样？”他们发现，该方法非常稳定。即使树脂层厚度在 0.6 × 10⁻³ 英寸 到 0.9 × 10⁻³ 英寸 之间变化，或者胶水刚度在 3 GPa 到 5 GPa 之间变化，结果也不会发生太大变化。这意味着工程师不需要了解胶水的微观细节就能得到准确的答案。

速度提升

或许最令人兴奋的部分是速度。重型的 3D 块模型（“实体单元”模型）虽然准确但很慢。为了获得理想的结果，它们需要为每一层至少切割出 8 个微小的块。然而，新方法在每一层只需要 1 个壳单元。

在测试中，新方法将计算机处理时间缩短了约 45%。对于一个特定的测试案例，旧方法每个计算步骤需要 9.52 秒，而新方法仅需 5.21 秒。当你正在设计一个拥有数千个零件的整个飞机机翼时，这是一个巨大的节省。

这对未来意味着什么

本文并不声称解决了宇宙中的所有问题。作者承认，对于单个层内部非常复杂的应力模式，他们使用了“线性插值”（一条直线）来猜测层与层之间的数值，这可能会错过一些微小的波动细节。然而，对于研究层间如何剥离和滑动的宏观情况，该方法表现得非常出色。

通过证明你可以使用简单、快速的壳单元获得准确的 3D 应力数据，这项工作为工程师更高效地模拟损伤和层间剥离打开了大门。这就像是从一个缓慢、模糊的望远镜升级到了一个清晰、快速的望远镜，让我们能在这些轻量化结构真正起飞之前，就看清其中隐藏的危险。作者暗示，这个框架可以成为未来同时模拟面内损伤（织物裂纹）和层间剥离（层间分离）工具的基础，从而让我们的飞机和汽车更加安全且轻量化。