达索SIMULIA是什么？仿真技术引领工业创新与数字化转型

这是一个在Abaqus显式动力学分析中非常经典且重要的问题。当你在sta文件或消息文件中看到动能（Kinetic Energy）大于总能量（Total Energy）时，这确实是一个需要警惕的信号，但它不一定意味着你的整个模拟都失败了。

下面我将详细解释为什么会发生这种情况、它意味着什么以及如何解决。

为什么动能会大于总能量？

从物理和数值计算的角度来看，这主要有以下几个原因：

1. 能量计算的数值误差（最主要、最常见的原因）

 · 总能量 在Abaqus中通常被计算为：总能量 = 内能 + 动能 + 伪应变能 + 滑移界面能 + ... - 外力做功。

 · 在一个完全守恒的系统中，总能量应该是一个常数。但在有限元计算中，这是不可能的。

 · 问题的核心在于： 当模型中的某些单元发生严重畸变（例如被删除的单元）或出现非常大的塑性变形时，这些单元对内能 的贡献计算会变得不准确，甚至可能被低估。

 · 如果计算出的内能 比实际值小，那么根据能量守恒公式，为了“平衡”，动能 的数值就可能被高估，从而导致动能暂时性地超过总能量。这可以看作是数值计算中的一种“能量泄漏”现象，内能“泄漏”到了动能中。

2. 沙漏效应

 · 在显式分析中，为了计算效率，经常使用减缩积分单元（如C3D8R）。这类单元存在一种非物理的变形模式，称为“沙漏模式”。在这种模式下，单元可以变形但不产生任何应变能（即内能不增加）。

 · 如果沙漏模式失控，节点的速度（动能）会急剧增加，但这些动能并没有转化为有意义的应变能。这会导致动能异常增高，甚至超过总能量。

3. 接触问题

 · 接触定义不当，例如在瞬间碰撞中产生巨大的、不合理的穿透，可能会导致速度（从而动能）的瞬时飙升。

 · 这也可以被视为一种数值上的能量异常。

4. 材料模型的不稳定性

 · 如果材料模型本身存在软化行为（如损伤、失效），并且参数设置不当，可能导致单元在失效时能量计算出现异常。

这说明了什么？有什么后果？

· 这表明你的模拟中可能存在严重的数值不稳定性。

· 结果的可靠性存疑： 当动能异常高时，特别是当它占总能量的比例很大时（例如超过5%-10%），说明模拟的物理过程已经被数值误差所主导。你得到的位移、应力、变形等结果很可能是不准确的。

· 可能预示着模型即将发散： 这种情况通常是模拟完全失败（因为单元过度畸变而导致计算终止）的前兆。

如何诊断和解决？

你需要系统地检查你的模型。以下是详细的排查步骤：

1. 检查沙漏能

· 操作： 在Abaqus/Viewer中绘制ALLIE（内能）、ALLKE（动能） 和ALLAE（伪应变能，即沙漏能+粘性耗散能） 的能量历史曲线。

· 判断标准：

· 沙漏能（ALLAE）应该只占内能（ALLIE）的很小一部分（通常认为1%-5%是可接受的，越小越好）。

· 如果你的沙漏能非常高，甚至和内力能在一个数量级，那么沙漏就是罪魁祸首。

· 解决方案：

· 使用沙漏控制： 在单元属性中，尝试使用增强的沙漏控制（例如 Hourglass = Enhanced 或 Stiffness）。

· 细化网格： 粗糙的网格更容易激发沙漏模式。

· 使用全积分单元： 将C3D8R单元改为C3D8单元。但要注意，这可能会引入剪切自锁问题，且计算成本更高。

· 施加沙漏阻尼。

2. 检查单元变形和删除

· 操作： 在Visualization模块中，打开State->PLY，并打开Show status。播放动画，观察是否有单元因失效准则（如塑性应变）被删除。

· 判断标准： 观察动能激增的时刻，是否正好对应大量单元被删除的时刻。

· 解决方案：

· 调整失效准则： 如果你的材料失效参数过于“脆弱”，可能导致单元过早删除。检查并校准你的材料失效模型。

· 使用更稳定的单元技术： 对于大变形问题，可以考虑使用Abaqus/Explicit中的欧拉域或CEL（耦合的欧拉-拉格朗日）方法，它们能更好地处理极端变形。

· 调整网格： 在可能发生大变形的区域进行局部网格细化。

3. 检查接触定义

· 操作： 检查接触设置，特别是主从面定义、滑移公式和接触刚度。

· 解决方案：

· 确保主从面定义合理（主面选择刚度更大、网格更粗的面）。

· 尝试使用“通用接触”（General Contact），它通常比“面面接触”（Surface-to-Surface Contact）更稳健。

· 如果穿透严重，可以适当增加接触刚度缩放因子（但过大会导致稳定时间增量减小，计算变慢）。

4. 检查加载速率和材料参数

· 问题： 虽然显式动力学可以处理准静态问题，但如果加载速度过快，会引入不必要的动态效应和惯性力，导致动能过高。

· 解决方案：

· 准静态分析： 确保加载时间足够长，使得动能（ALLKE）在整个过程中始终远小于内能（ALLIE），通常建议动能/内能 < 5%。

· 检查材料数据： 确保你的材料参数（尤其是塑性数据）来自可靠的实验，并且单位制一致。

5. 检查质量缩放

· 问题： 为了提高计算效率而使用的质量缩放，如果使用不当，会人为地增加模型的动能。

· 解决方案：

· 避免使用固定时间增量步的质量缩放，因为它会全局地、永久地改变模型的质量。

· 如果使用，尽量将其限制在那些稳定时间增量非常小的特定区域（通过Set来定义）。

· 始终监控由质量缩放添加的总质量，确保它只占总质量的很小一部分。

总结

当遇到“动能大于总能量”的问题时，你的排查流程应该是：

1. 首先绘制能量历史曲线，重点关注沙漏能（ALLAE）的比例。

2. 观察动画，看动能激增的时刻发生了什么（单元删除、严重扭曲、接触穿透）。

3. 根据观察到的现象，有针对性地采取上述解决方案。

记住，在显式动力学分析中，能量平衡 是判断结果可靠性的黄金准则。一个健康的模拟，其总能量应该平滑变化，动能（对于准静态问题）或沙漏能应该保持在很低的水平。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删