在进行显式动力学仿真时,**能量平衡(Energy Balance)**是判断仿真结果可靠性的“金标准”。而在众多能量参数中,**接触能(Contact Energy)**往往是最让工程师头疼的变量之一。
今天我们就来拆解一下 LS-DYNA 中接触能的构成,以及如何通过它来诊断你的模型健康状况。

在LS-DYNA 中,接触能主要记录在 SLEOUT 文件中(需通过 *DATABASE_SLEOUT 开启)。打开这个文件,你会看到四个核心指标:
从动能 (Slave Energy, SE)
主动能 (Master Energy, ME)
摩擦能 (Frictional Energy, FE)
主从能量总和 (Sum of Slave/Master Energy, SoSM):即
SoSM = SE + ME
划重点:
我们在 GLSTAT 文件中看到的 “Sliding Energy”,其实就是 SLEOUT 中所有接触定义的 SoSM 之和。
注意: 摩擦能(FE)已经被包含在 SE 和 ME之中了。如果你想知道单纯为了防止穿透而产生的能量(即法向接触能),你需要手动计算:
E_{normal} = SoSM - FE
在无摩擦的情况下,SE 和 ME 应当大小相近、符号相反,它们的总和 SoSM 应当等于系统存储的接触势能。
当模型中定义了摩擦(如 *CONTACT 中的 FS 和 FD)时,正的接触能是完全正常的。摩擦本质上是做功产生热量的过程,必然会导致能量的增加。
多大算大? 在没有摩擦的情况下,我们希望看到一个很小的净接触能(Net Contact Energy)。工程经验上,如果接触能占到系统峰值内能(Internal Energy)的 10% 以内,通常是可以接受的。
相比正能量,负接触能通常预示着模型存在缺陷。它就像一个不稳定的黑洞,可能会导致计算发散。
初始穿透(Initial Penetration): 这是最常见的原因。如果几何模型没考虑壳单元偏置(Shell Offset),节点一开场就陷在里头,系统会产生巨大的排斥力。
节点“滑脱”: 当一个已穿透的节点从一个主段(Master Segment)滑动到另一个不连续的主段时,系统会瞬间检测到新的穿透,从而产生负功。
观察 GLSTAT 中的曲线。如果 内能曲线 的斜率与 负接触能曲线 正好相反(此消彼长),说明问题高度局部化。
定位工具: 在 LS-PrePost 中使用 Fcomp > Misc > internal energy 查看壳单元。那些能量异常爆表的“热点”区域,通常就是负接触能的源头。
如果你的 SLEOUT 报告了令人不安的负能量,请尝试以下组合拳:
清理几何: 彻底消除初始穿透,仔细检查消息文件(Message File)中的 Warning。
检查冗余: 确保同一对表面没有定义两个重复的接触(Redundant Contact)。
缩小时步: 减小时间步缩放系数(Time step scale factor)。
优化接触参数(推荐配置):
将控制参数恢复默认,但尝试设置 SOFT=1 且 IGNORE=2(可选卡片 C)。
针对尖锐边缘: 在可选卡片 A 中设置 SOFT=2(基于段的接触)。
精细化 SOFT=2 参数:
SBOPT: 若无明显相对滑动,设为 3;若滑动明显,设为 5。
DEPTH: 推荐设为 13 或 23。如果需要考虑壳单元边缘对边缘(Edge-to-Edge)的接触,请设为 25 或 35。
接触能不是一个冷冰冰的数字,它是模型逻辑一致性的体现。当你通过 SLEOUT 锁定问题接触对,并配合 SOFT 选项进行微调时,你离一个完美的仿真结果就不远了。