首页/文章/ 详情

深度解析 LS-DYNA 接触能:从 SLEOUT 到能量平衡的“排雷”指南

7月前浏览543

在进行显式动力学仿真时,**能量平衡(Energy Balance)**是判断仿真结果可靠性的“金标准”。而在众多能量参数中,**接触能(Contact Energy)**往往是最让工程师头疼的变量之一。

今天我们就来拆解一下 LS-DYNA 中接触能的构成,以及如何通过它来诊断你的模型健康状况。


1. 拨开迷雾:接触能的四种面孔


在LS-DYNA 中,接触能主要记录在 SLEOUT 文件中(需通过 *DATABASE_SLEOUT 开启)。打开这个文件,你会看到四个核心指标:

  1. 从动能 (Slave Energy, SE)

  2. 主动能 (Master Energy, ME)

  3. 摩擦能 (Frictional Energy, FE)

  4. 主从能量总和 (Sum of Slave/Master Energy, SoSM):即

     SoSM = SE + ME

划重点:

我们在 GLSTAT 文件中看到的 “Sliding Energy”,其实就是 SLEOUT 中所有接触定义的 SoSM 之和。

注意: 摩擦能(FE)已经被包含在 SE 和 ME之中了。如果你想知道单纯为了防止穿透而产生的能量(即法向接触能),你需要手动计算:

E_{normal} = SoSM - FE

在无摩擦的情况下,SE 和 ME 应当大小相近、符号相反,它们的总和 SoSM 应当等于系统存储的接触势能。



2. 正接触能:摩擦是“罪魁祸首”吗?


当模型中定义了摩擦(如 *CONTACT 中的 FS 和 FD)时,正的接触能是完全正常的。摩擦本质上是做功产生热量的过程,必然会导致能量的增加。

  • 多大算大? 在没有摩擦的情况下,我们希望看到一个很小的净接触能(Net Contact Energy)。工程经验上,如果接触能占到系统峰值内能(Internal Energy)的 10% 以内,通常是可以接受的。



3. 负接触能:仿真中的“红色警报”


相比正能量,负接触能通常预示着模型存在缺陷。它就像一个不稳定的黑洞,可能会导致计算发散。

产生原因:

  1. 初始穿透(Initial Penetration): 这是最常见的原因。如果几何模型没考虑壳单元偏置(Shell Offset),节点一开场就陷在里头,系统会产生巨大的排斥力。

  2. 节点“滑脱”: 当一个已穿透的节点从一个主段(Master Segment)滑动到另一个不连续的主段时,系统会瞬间检测到新的穿透,从而产生负功。

诊断技巧:

观察 GLSTAT 中的曲线。如果 内能曲线 的斜率与 负接触能曲线 正好相反(此消彼长),说明问题高度局部化。

  • 定位工具: 在 LS-PrePost 中使用 Fcomp > Misc > internal energy 查看壳单元。那些能量异常爆表的“热点”区域,通常就是负接触能的源头。



4. 实战方案:如何消除负接触能?


如果你的 SLEOUT 报告了令人不安的负能量,请尝试以下组合拳:

  1. 清理几何: 彻底消除初始穿透,仔细检查消息文件(Message File)中的 Warning。

  2. 检查冗余: 确保同一对表面没有定义两个重复的接触(Redundant Contact)。

  3. 缩小时步: 减小时间步缩放系数(Time step scale factor)。

  4. 优化接触参数(推荐配置):

    • 将控制参数恢复默认,但尝试设置 SOFT=1 且 IGNORE=2(可选卡片 C)。

    • 针对尖锐边缘: 在可选卡片 A 中设置 SOFT=2(基于段的接触)。

  5. 精细化 SOFT=2 参数:

    • SBOPT: 若无明显相对滑动,设为 3;若滑动明显,设为 5。

    • DEPTH: 推荐设为 13 或 23。如果需要考虑壳单元边缘对边缘(Edge-to-Edge)的接触,请设为 25 或 35。



结语

接触能不是一个冷冰冰的数字,它是模型逻辑一致性的体现。当你通过 SLEOUT 锁定问题接触对,并配合 SOFT 选项进行微调时,你离一个完美的仿真结果就不远了。


来源:ABAQUS仿真世界
ACTLS-DYNA显式动力学UM控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-03-12
最近编辑:7月前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 206粉丝 401文章 500课程 0
点赞
收藏
作者推荐

材料失效准则该怎么选?

所有工程材料(如金属、塑料、玻璃和混凝土)的使用寿命都是有限的。疲劳、腐蚀、侵蚀、应力腐蚀开裂或异物损坏等各种因素会影响材料的强度、稳定性和化学成分,从而导致其失效。这篇文章旨在总结一些更常见的失效理论,并提供根据材料类型使用哪些理论的指导......一、介绍在所有外部因素中,腐蚀和机械故障是材料故障的主要原因。机械故障通常分为以下两种类型:延展性失效——这种失效发生在断裂前经历大幅变形的金属中。材料在断裂前存在巨大的塑性变形和明显的颈缩。脆性破坏——这种类型的破坏发生在塑料、岩石、混凝土、玻璃和复合材料等材料中。材料突然断裂而没有任何明显的变形。这些材料在断裂前塑性变形非常小或没有。预测单轴载荷下样本的失效是很简单的。当单轴应力达到延性材料的屈服强度和脆性材料的极限强度时,它就会失效。但当结构承受复杂载荷时,预测失效是很困难的。不存在适用于所有材料的通用失效理论,失效理论的预测必须与实验观察结果一致。本文讨论了延性材料和脆性材料的不同失效理论。比较不同失效理论之间平面应力条件下的屈服面。二、延性材料的失效理论屈服不会由静水应力引起,偏应力通过剪切材料导致形状扭曲。这种偏应力是导致延性材料屈服的原因。良好的延性材料失效理论必须仅考虑偏应力来确定是否断裂。Tresca准则(最大剪切应力准则)该理论指出,当最大剪切应力的绝对值达到简单单轴拉伸试验中屈服时的剪切应力时,材料就会发生失效。该理论表明断裂与静水应力无关。对于平面应力,Tresca准则的屈服面在主应力空间中表示,如下图所示。VonMises屈服准则(最大畸变能准则)该理论指出,当最大应变能等于单轴拉伸试验中屈服时的变形能时,材料就会发生失效。变形能是由材料中存在的剪切应力引起的。使用主应力计算的公式如下:在单轴拉伸试验中,屈服变形能计算如下:根据上述方程,可得失效准则在数学上表示为:该屈服面在平面应力条件下是椭圆形的。VonMises理论与实验数据更加吻合。Tresca屈服面被VonMises屈服面包围。因此,Tresca是更为保守的失效准则。Hosford屈服准则这是用于确定材料在三轴应力状态下塑性变形开始的失效准则。它在数学上表达为,当a=1时,上式代表Tresca准则,a=2时代表VonMises屈服准则。这被视为Tresca和VonMises屈服面之间的插值,如下图所示。Gurson失效准则该屈服准则适用于具有一般椭圆形空隙的延性塑性多孔固体材料。从数学上讲,该失效理论表示为:VonMises理论与Gurson的不同之处在于压力相关函数phi(p,α,f),该函数展开为:当上式中f=0时,它即变成vonMises准则。该材料的屈服应力以一维初始屈服应力y0、相关应变α0、杨氏模量E、硬化指数H和等效塑性应变α表示,如下所示:不同空隙体积分数f值的屈服面与f=0时vonMises屈服面的关系如下图所示。三、脆性材料的失效理论脆性材料的失效是由于断裂而不是屈服。在这些材料中,压缩强度高于拉伸强度。因此,失效理论必须能表达拉伸和压缩的这种差异。Rankine失效理论(最大主应力)该理论假设当最大或最小主应力达到材料的失效强度(σu)时,材料就会发生失效。这是一个简单的理论,但它并没有考虑屈服仅取决于偏应力。因此,这仅适用于脆性材料,不适用于延性材料。Coulomb-Mohr失效准则该理论考虑了失效对静水应力的敏感性,并解释了压缩断裂强度和拉伸断裂强度之间的差异。这假设最大和最小主应力平面包含失效平面的法向。在数学上,该理论表示为:在平面应力条件下,屈服面如下图所示。ModifiedMohr失效准则库仑-莫尔理论与实验数据不太吻合。因此,建议进行修改,需要3个材料常数。这是不同类型土壤的首选理论。它捕获脆性材料中与压力相关的剪切强度、不可恢复的压实和非线性弹性卸载。从数学上讲,这与莫尔方程相同,修改后表示为:在平面应力条件下,屈服面如下图所示。四、最后的想法了解物体失效所需的应力大小非常重要。希望本文能够很好地洞察延性破坏理论和脆性破坏理论之间的区别。它们有助于预测承受负载的不同机械和物理部件的寿命和功能。来源:ABAQUS仿真世界

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈