Stability of Extrinsic Cohesive-Zone Model with Penalty-Based Contact in Explicit Dynamic Fragmentation Simulations

本研究指出，在显式动力学破碎模拟中将外在内聚力模型与基于惩罚项的接触结合使用，会导致由于刚度不连续和切换误差而引起的严重非物理能量增长和人工破碎，并最终得出结论：尽管自适应惩罚策略能提供部分缓解，但该方法仍不适用于长期的、能量一致性的模拟。

要点

大局观：模拟玻璃破碎的过程

想象一下，你正试图预测一块石头击中窗户时，玻璃是如何破碎的。你不仅想知道它会碎，还想知道它会碎成多少片、碎片有多大、飞出的速度有多快。为了实现这一目标，科学家们使用计算机模拟。

这篇论文研究了一种专门用于高速爆炸或冲击的特定计算机模拟方法。研究人员发现，他们的模拟程序在向他们“撒谎”。模拟并没有显示出稳定的破碎过程，而是产生了虚假的、无止境的破碎，并凭空创造出了能量。

他们旨在找出：为什么计算机出现了故障？我们该如何修复它？

设定： “胶水”与“弹簧”

为了模拟破碎过程，研究人员在计算机模型中使用了两种主要工具：

1. “胶水”（内聚力区模型/Cohesive Zone Model）： 想象这种材料是由微小的乐高积木组成的。在积木之间，存在着一种看不见的、具有拉伸性的“胶水”。当你拉开积木时，胶水会发生拉伸，并最终断裂。这模拟了裂纹是如何产生和扩展的。
2. “弹簧”（惩罚接触/Penalty Contact）： 一旦胶水断裂且积木分离，它们可能会弹回并撞击彼此。为了防止它们互相穿透（这在物理上是不可能的），计算机使用了一个“弹簧”规则。如果两个积木试图重叠，弹簧就会将它们推开。弹簧越硬，就越难发生重叠。

问题：“蹦蹦床效应”

当他们运行模拟时，计算机的行为就像一个永不停歇的蹦蹦床。

• 症状： 模拟中的总能量不断升高，尽管并没有添加新的能量。
• 结果： 计算机认为材料正在破碎成数百万个微小的、不可能存在的碎片。“碎片计数”（碎片数量）不断无限上升，这在物理学上是不可能的。

研究人员问道：是胶水太弱了？是弹簧太硬了？还是数学本身出了问题？

调查：三个嫌疑对象

团队测试了导致故障的三个可能原因，就像侦探排除嫌疑人一样。

嫌疑人 1：“全新的胶水”（发散的初始刚度）

理论： 当一段胶水刚刚产生时（在它发生拉伸之前），它是极其坚硬的。理论上，它是无限硬的。
测试： 他们检查了这种超强刚度是否导致了计算机计算的不稳定性。
结论： 并非主要元凶。 虽然这种情况可能会导致问题，但在他们的特定测试中，胶水并没有变得足够硬以至于破坏模拟。这是一个误导性的线索。

嫌疑对象 2：“软化”（逐渐减弱）

理论： 随着胶水的拉伸和断裂，它会变得越来越弱（软化）。也许这种强度的变化让计算机产生了困惑。
测试： 他们分析了胶水变弱时的数学过程。
结论： 无罪。 数学证明，当胶水变弱时，损失的能量与用于产生新裂纹表面的能量是完美平衡的。这部分模拟实际上运行正常。

嫌疑对象 3：“切换”（内聚力-接触转换）—— 真凶

理论： 这才是真正的问题所在。想象材料的一个碎片正在振动。它在拉伸（胶水模式），然后又弹回并接触到另一个碎片（接触模式）。

• 在胶水模式下，材料表现得像一种特定类型的弹簧。
• 在接触模式下，材料表现得像另一种类型的弹簧（惩罚弹簧）。

问题在于，每当材料接触时，计算机必须瞬间从一种弹簧规则切换到另一种。这就像开车时，每撞到一个颠簸，车子就会突然从“踩油门”切换到“踩刹车”。
结果： 每当材料在“胶水”和“接触”之间切换时，计算机都会产生一个微小的数学误差。它无意中增加了一点点能量。

• 类比： 想象一个在秋千上荡的孩子。每当孩子荡到最高点时，你都会不小心轻轻推他一下。起初你感觉不到，但荡了1000次后，孩子飞得太高，甚至撞到了天花板。
• 现实情况： 在模拟中，这些微小的能量误差在数百万步中不断累积，导致了“虚假能量”的爆发和无止境的破碎。

提出的“修复方案”及其为何并非真正的解决方案

研究人员尝试了一个聪明的技巧来阻止这种故障。他们让“接触弹簧”的刚度与“胶水弹簧”完全一致。

• 结果： 突然的“切换”消失了。能量不再增长。模拟变得稳定了。
• 代价： 为了让弹簧匹配，当胶水受损时，“接触弹簧”必须变得非常弱。这意味着材料的碎片被允许发生显著的重叠（相互穿透）。
• 结论： 虽然这个方法解决了数学上的故障，但它破坏了物理规律。你不能为了让数字看起来正确，就让固体碎片互相穿透。因此，这个“修复方案”对于诊断问题很有用，但它并不是一个真正的工程解决方案。

最终总结

论文得出结论：在高速破碎模拟中使用“惩罚弹簧”来处理接触问题，在长期准确性方面存在根本性的缺陷。

• 权衡： 你无法鱼与熊掌兼得。如果你让接触弹簧非常硬以防止碎片重叠，你会迫使计算机采取极小、极慢的步长。如果你让它变软以提高速度，你就会得到能量误差和虚假的破碎。
• 未来方向： 作者建议，与其使用“软弹簧”（惩罚法），我们更需要使用“硬规则”（非光滑力学），将接触视为一种严格的法则而非一种弹簧。这将阻止能量泄漏，并实现对物体破碎过程的精确、长期的模拟。

简而言之： 计算机之所以产生了一场永无止境的爆炸幻觉，是因为每当破碎的碎片与其他碎片碰撞时，它都会感到困惑。用于防止碎片重叠的“弹簧”方法正是这种困惑的根源，而真正解决问题的办法是彻底改变计算机处理碰撞的方式。

技术摘要

技术摘要：外在内聚力模型结合基于惩罚项接触在显式动力学破碎模拟中的稳定性

问题陈述
动力学破碎模拟对于预测高应变率下的材料响应（如轨道碎片场景中的超高速碰撞）至关重要。虽然结合了外在内聚力模型（CZM）与基于惩罚项接触的显式动力学模拟因其可扩展性而广受欢迎，但它们经常表现出严重的数值不稳定性。在基准测试中，这些不稳定性表现为指数级的能量增长和人为的、非物理性的破碎，从而使碎片质量和速度分布等统计输出失效。尽管这一问题普遍存在，但在显式时间积分中，外在 CZM 与惩罚项接触耦合所导致的不稳定性的具体根源尚未得到系统的隔离或量化。

方法论
作者利用在开源代码 Akantu 中实现的数值框架，采用具有线性牵引-分离律（TSL）的外在 CZM，并结合基于惩罚项的接触公式。模拟使用显式中心差分（CD）时间积分方案。

为了诊断不稳定性，作者：

1. 建立了基准： 模拟了一个二维正方形平板（类似于爆炸的空心球体）的动力学破碎过程，使用 AD-995 氧化铝材料属性，其内聚强度按 Weibull 分布进行采样。
2. 隔离不稳定性机制： 使用简化模型系统地分析了三个潜在的不稳定来源：
   • 发散的初始刚度： 从理论上分析新插入的内聚单元在损伤接近于零（无穷大刚度）时对稳定时间步长的影响。
   • 内聚–接触切换： 使用单自由度（SDOF）系统对内聚刚度（$k^+$）与接触惩罚刚度（$k^-$）之间的突变进行建模，以量化每步产生的能量误差。
   • 软化： 分析内聚刚度逐渐降低（损伤演化）过程中的能量平衡。
3. 量化能量漂移： 使用解析误差估计、相空间诊断和能量增长指标，确定微小的单步误差如何累积成长期的能量漂移。
4. 测试缓解措施： 评估了一种自适应惩罚策略，其中接触刚度（$k^-$）与演化的内聚刚度（$k^+$）动态绑定，以消除刚度不连续性。

主要贡献与结果

• 识别三种机制： 研究隔离并量化了三个不同的不稳定来源：

  1. 发散的初始内聚刚度： 当损伤趋于零时，内聚刚度趋于无穷大，从而限制了稳定时间步长。然而，在所使用的特定高应变率基准测试中，该机制并非观察到的全局不稳定的主要原因，因为内聚单元会迅速越过不稳定区域。
  2. 不连续的刚度跳跃（主要原因）： 在界面处，内聚刚度（$k^+$）与固定接触惩罚（$k^-$）之间的突变造成了牵引-分离响应的不连续。作者证明，在这些状态之间的反复切换会产生交替的能量注入与耗散。在特定的参数范围（刚度比和时间步长）下，这些微小的单步误差会累积，导致无界的能量增长和人为的破碎，即使在满足标准库朗–库朗（CFL）稳定性条件的情况下也是如此。
  3. 由软化引起的连续性缺失： 虽然软化引入了非光滑路径，但研究发现由此产生的算法能量损失与物理断裂耗散完全平衡。因此，软化事件本身是数值保守的，不会触发不稳定性。

• 量化“内聚–接触”不稳定性： 通过使用 SDOF 模型，作者表明，针对具有切换刚度的分段线性系统，经典的稳定性界限（针对线性时不变系统推导）是不充分的。作者绘制了刚度比（$\alpha_{ref} = k^+/k^-$）和归一化时间步长参数空间中的稳定与不稳定区域。结果表明，对于具有显著刚度对比的频繁切换，时间步长必须显著降低至低于标准极限（例如 $\Delta t \lesssim 0.2 \Delta t_c$）才能保持稳定。

• 评估自适应惩罚： 作者测试了一种诊断性修改，即将接触惩罚刚度与内聚刚度绑定（$k^- = k^+(d)$）。这种修改成功消除了刚度不连续性，恢复了能量守恒，并稳定了碎片计数。然而，随着损伤增加，它会引入显著的非物理穿透，使其不适合作为最终的物理补救措施，但作为确认根因的诊断工具非常有效。

意义与主张
论文得出结论，在外在 CZM 与标准惩罚项接触的耦合中，对于需要具有物理意义的碎片统计数据的长时程、能量一致的破碎模拟，该方法并非可行方案。惩罚法的内在局限在于刚度权衡：高刚度防止穿透但限制了时间步长并引入了通过不连续性产生的能量伪影；低刚度允许更大的时间步长但允许穿透。

作者断言，观察到的非物理能量漂移的根源是重复的内聚–接触切换导致的误差累积。他们认为启发式修复（如局部时间步长缩减、质量重分布）仅是缓解症状。相反，本文提倡采用非光滑时间步进方案，在速度层面上强制执行接触（使用冲量和互补约束）。此类方法将消除人工刚度跳跃，使稳定性限制回归到整体动力学和内聚动力学，并允许进行稳健、长时间的模拟，而不产生与基于惩罚项接触相关的能量漂移。