ANSYS LS-DYNA常见问题之Sliding Energy问题

3年前浏览4392

编者按

在ANSYS LS-DYNA中，经常遇到Sliding Energy异常，作者在本文中分析了计算模型常出现的2种情况——初始穿透导致的Sliding Energy异常和由于SEGMENT接触导致的Sliding Energy异常，并给出了对应的解决方法。

01问题提出

负的滑移能存在2种情况

■ 第一种情况：GLSTAT中的滑移能在初始状态就出现负值，如图1所示。

图1 第一种能量异常图

■ 第二种情况：CAE工程师在利用ANSYS LS-DYNA进行计算时，很多时候会碰到图2所示的状况，严重的能量不平衡。内能、滑移能在计算前期表现还算正常，但是到计算后期，内能、动能、滑移能严重偏离，导致整个模型能量异常，计算结果不可信。

图2 第二种能量异常图

02解决方案

第一种情况的解决方法（初始穿透）

在仿真模型中出现初始穿透，导致初始的滑移能为负值。仿真穿透模型如图3所示。

产生机理：LS-DYNA在进行接触搜索时，发现存在网格穿透问题，需要把穿透部分的节点移到**透的位置，这可以在LS-DYNA软件中设置。如图4所示，当把初始穿透清除后，其滑移能和整个模型能量表现正常（如图5所示）。

图3 存在初始穿透模型

图4 清除初始穿透模型

图5 第一种能量正常图

初始穿透解决方法有以下几种

手动调整单元节点

在进行模型网格划分和设置零件厚度时必须确认是否存在干涉，即必须考虑壳单元的接触厚度。

如果发生穿透，可以从计算出的message或d3hsp文件中对关键字“initial penetrations”进行搜索，找到相关单元，然后调整单元节点，消除初始穿透。

在一些专门的前处理软件中，例如ANSA，HYPERMESH、SpaceClaim在提交计算前对模型进行穿透检查，可以查出初始穿透的单元，然后进行节点移动，消除穿透。

图6在message文件中查找初始穿透

调整关键字参数

■ 调整*CONTACT中的SFST和SFMT，即调整该接触对从面和主面的接触厚度的比例因子。该方法对微小的初始穿透效果很明显，但是对于初始穿透值很大，效果不好，容易导致计算结果错误。

图7 接触卡片中的SFST和SFMT

■ 设置*CONTACT中的IGNORE参数，该参数默认值为0，即采用*CONTROL_CONTACT中IGNORE的值。可以对每个接触对进行修改，建议修改值为1或者2。

图8 接触卡片中的IGNORE

第二种情况的解决方法

众所周知，ANSYS LS-DYNA存在2种方式设置接触，即有Contact Region和Body Interaction。

■Body Interaction接触即LS-DYNA接触关键字*CONTACT_AUTOMATIC_SINGLE_SURFACE，主要用于整个模型的接触或自接触。

■ Contact Region 用于建立接触对，当几何模型导入到ANSYS以后，Mechanical默认对每个零件建立bond接触。其输出关键字为*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET ，其set类型为segment set。

对于接触面和目标面这两个segment进行接触搜索时，在交界处的节点无法检查出穿透，然而当节点滑到接触厚度中时，就会发现穿透节点的存在，因此系统必定会给该节点施加一个力，把它拉回到接触面上，在这个过程中系统对该节点做功，消耗其接触势能，能量曲线输出表现为负的滑移能。

其解决方法之一是扩大主段的接触面。

由*CONTACT中的关键字MAXPAR进行调整。MAXPAR默认值为0，即表示设置值为1.025，其LS-DYNA推荐值为1.2。

图9 接触卡片中的MAXPAR

其解决方法之二是调换接触算法。

ANSYS LS-DYNA在Mechanical中设定的bond接触是*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET。在LS-DYNA接触选项中带OFFSET，其目的是把基于约束类型的接触算法切换到基于罚函数的接触算法。

因此，当出现滑移能异常，可以把*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET切换为

*CONTACT_TIED_SURFACE_TO_SURFACE。

下文是一个简单的对比，使用bond接触的模型如图10所示，其采用不同的接触算法的能量图如图11所示。

图10 使用Bond接触的模型（黑色部分为接触区域）

（a）采用*CONTACT_TIED_SURFACE_TO_SURFACE_OFFSET

（b）采用*CONTACT_TIED_SURFACE_TO_SURFACE

图11 不同接触算法下的能量平衡图

Bond接触推荐使用类型

■ 实体单元到实体单元的bond接触 也就是说，对于没有节点具有旋转自由度的情况，使用TIED_NODES_TO_SURFACE和TIED_SURFACE_TO_SURFACE类型的接触。这些接触类型可能包括OFFSET或CONSTRAINED_OFFSET选项。

■ 壳单元与壳单元、beam与壳单元的bond接触 也就是说，对于所有节点都具有旋转自由度的情况，可以使用TIED_SHELL_EDGE_TO_SURFACE类型的接触。此接触类型可能包括OFFSET、CONSTRAINED_OFFSET或BEAM_OFFSET选项。

■ TIED_SHELL_EDGE_TO_SOLID 用于将壳单元边与实体单元或beam端与实体单元绑定，即只有从端节点具有旋转自由度的情况。

03总结和建议

此部分，需要用户使用Mechanical中的Contact Property设置接触属性。使用ANSYS LS-DYNA中总是不可避免遇到能量问题，因其涉及判断模型是否正常，计算结果是否正确，所以对能量的判断必须符合要求

■ 对于小的初始穿透，在对模型与试验对标后，没有出现误差，可以采用IGNORE进行处理；但如果影响计算结果，请务必调初始穿透，这要求CAE工程师在网格处理时必须进行网格检查

对于第二种情况，利用GLATAT和SLEOUT查出那一对接触出现能量异常，然后修改参数或者接触算法

参考文献：①LS-DYNA®Theory Manual 07/24/19 (r:11261) [M] 2019

②LS-DYNA_Manual_Volume_I_R11 [M] 2018

③赵海鸥LS-DYNA 动力分析指南[M] 北京：兵器工业出版社，200

仿真体系网格处理 LS-DYNA Mechanical

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2020-10-20