第十七章LS-DYNA的隐式求解

1、第十七章LS-DYNA的隐式求解LS-DYNA作为著名的显式求解程序只能求解瞬态动力问题，对于时间历程较长 的静、动力问题，LS-DYNA的显式中心差分法有它的局限性，而一些与瞬态动 力分析紧密相关的问题要求 LS-DYNA也能够求解，如： 冲压成型后的回弹计 算应力初始化冲击后常时间低频动力响应静力分析特征值分析实际上从950版本开始，LS-DYNA已增加了隐式求解功能。刚开始的应用主要 在冲压成型后的回弹计算，经过 960版到970版的发展，LS-DYNA的隐式求解 功能已大增强，已经能满足以上的求解需要。17.1 显式与隐式的区别：17.1.1 LS-DYNA 显式求解：采用中心差分方法

2、进行显式时间积分Man fn方程非耦合，可以直接求解(显式)但需要常小的时间步保持稳定状态不需要求解刚度矩阵ext fnint适合冲击、穿透等高频非线性动力响应问题17.1.2 LS-DYNA 隐式求解：采用Newmark隐式时间积分M an K un 1f对于线性问题，无条件稳定可以采用大的时间步extn 1 fintn Man对于非线性问题，需要一系列线性逼近(Newton-Raphson)叠代求解需要求解刚度矩阵适合静力问题、低频动力问题及特征值分析。17.2 LS-DYNA中隐式分析的激活及相关关键字在LS-DYNA中，缺省的求解是显式求解，可以通过下面的关键字来激活隐式求解：*CON

3、TROL_IMPLICIT_GENERAL*CONTROL_IMPLICIT_GENERAL$ imflag dt0 iefs nstepsb igso1 0.01 0 0 0其中参数imflag=1激活全隐式求解imflag=0为缺省的显式求解。imflag=2为显式求解后无缝进行隐式回弹求解。该关键字对于所有隐式求解分析来说都是必需的。与隐式求解相关的其它关键字：*CONTROL_IMPLICIT_LINEAR（v960 版本改为 *CONTROL_IMPLICIT_SOLVER）*CONTROL_IMPLICIT_NONLINEAR（在 v960版后改为 *CONTROL_IMPLICI

4、T_SOLUTION）*CONTROL_IMPLICIT_AUTO*CONTROL_IMPLICIT_STABILIZATION*CONTROL_IMPLICIT_DYNAMICS*CONTROL_IMPLICIT_EIGENVALUE*CONTROL_IMPLICIT_BUCKLE （v970）*CONTROL_IMPLICIT_MODES （v970）使用*CONTROL_IMPLICIT_GENERAL 激活隐式求解后，还需要与上面的相关 关键字进行联合，从而对不同的问题设置不同的关键字进行求解。*CONTROL_IMPLICIT_LINEAR（v960 版后改为 *CONTROL_IM

5、PLICIT_SOLVER）*CONTROL_IMPLICIT_SOLVER$ lsolver prntflg negeig0 0 0$该关键字用来设置线性方程求解器的相关参数，主要用来求解下列线性方程： 要注意的是该关键字不是析激活线性分析，该求解器的功能主要是对刚度矩阵进 行转置计算（线性分析还是非线性分析由下面的关键字指定）。*CONTROL_IMPLICIT_NONLINEAR（在 v960版后改为 *CONTROL_IMPLICIT_SOLUTION）*CONTROL_IMPLICIT_SOLUTION$ nlsolvr ilimit maxref dctol ectol rctol

6、 lstol0 0 0 0.0 0.0 0 0$ dnorm divflag inistif nlprint0 0 0 0该关键字是用来指定是线性分析还是非线性分析，若是非线性分析，则指定非线 性叠代求解器类型、控制平衡叠代搜索次数和各种容限。若设置nlsolvr=1 ,则为线性分析，实际上相当于叠代次数为1的非线性分析。若设置nlprint=1 ,则屏幕上输出如下，在求解过程中，也可以使用 CTRL+ nlprint 来输出。K un R 4i'lj LM L I(fe Li IL1. li： IL J. lillbll:! 12： lit J 1JI lLl ICLIF 1IL1i

7、： I ph iI d ：> b, 'ji / l亚阴 址！IM热*潮 EMM I： O” 小 “一 >步 ii, X aUi»ll»>l：liGIIEJ-nil. v. 4 ,即CU：E£It U旧 uJIHI IRli：； II Ml i Hriiijp "I W tfy il. l: K: E tiv.np/ rj"：1p:小dI ii I, ii，|» | inii I; | I >L i l“Hi :；n I 嘈 I U1 I" M；：卜 l| 温闾曲7升学瞄UJI；泡id, 4；

8、<：il « tin hir>jL . u: 3 .狗 i和“灯 E m hjl i. i j 苫胖 H ” 的璘 n-j if* im . 咱 1*CONTROL_IMPLICIT_AUTO*CONTROL_IMPLICIT_AUTO$ iauto iteopt itewin dtmin dtmax1 0 0 0.0 0.0该关键字控制激活自动步长控制，缺省是固定时间步长，当平衡叠代失效后，可 以试着改为自动时间步长。*CONTROL_IMPLICIT_STABILIZATION该关键字仅用于多步回弹分析，施加人工粘度，提高回弹分析的收敛性，具体应 用见冲压成型分析章节

9、。*CONTROL_IMPLICIT_DYNAMICS*CONTROL_IMPLICIT_DYNAMICS$ imass gamma beta0 0.0 0.0该关键字激活是静力还是动力分析的选项。Imass=0为静力分析Imass=1低频动力分析（采用Newmark时间积分）Imass=2低频动力分析（采用模态综合法）*CONTROL_IMPLICIT_EIGENVALUE*CONTROL_IMPLICIT_EIGENVALUE$ neig15该关键字激活特征值分析，可输出各阶模态，neig=15表示输出前15阶模态。不同的应用需要的关键字不一样，如做特征值分析，需要3个关键字：*CONTR

10、OL_IMPLICIT_GENERAL*CONTROL_IMPLICIT_EIGENVALUE*CONTROL_TERMINATION建议使用双精度的求解器，单元公式用18号壳单元公式(下面有具体说明)，求解完后自动生成两个文件：D3EIGV :与D3POLT文件一样的二进制输出文件，包含每一模态下的振型，可 以用LS-PREPOST显示。EIGOUT: ASCII文件，是所有输出模态的文档总结。后面有专门的示例来说明这种分析并怎样进行后处理17.3 LS-DYNA隐式单元公式：一般情况下，不建议使用LS-DYNA缺省的显式单点积分单元来求解隐式问题， LS-DYNA开发有专门的各种隐式单元应

11、用于隐式分析中，下面具体介绍各种单 元：使用*SECTION_OPTION 来定义。17.3.1 隐式壳单元公式：1. S/R Hughes-Liu单元：6号单元公式2 x2选择性缩减积分 每个节点6个自由度(dx, dy, dz, rx, ry, rz)在*CONTROL_SHELL2 .快速壳：16号单元公式每个节点6个自由度(dx, dy, dz, rx, ry, rz)个积分点不能退化为三角形利用沙漏控制8增加翘曲刚度，改善收敛性3 .线性壳：18号单元公式每个节点6个自由度(dx, dy, dz, rx, ry, rz)由于是线弹性，厚度方向2个积分点足够面内2 x2积分点，厚方向建

12、议5非常适合用于回弹计算面内2 x2积分点包含三角形壳公式仅应用于线性问题和特征值分析必须使双精度求解器17.3.2隐式体单元公式：1 .积分S/R体单元：2号单元公式每个节点3个自由度(dx, dy, dz)2 x2 x 2选择性缩减积分CPU开销大2,线性体单元：18号单元公式每个节点3个自由度(dx, dy, dz)3 x2 x 2积分点仅应用于线性问题和特征值分析必须使用双精度求解器1.1.3 3隐式梁单元公式：每个节点6个自由度(dx, dy, dz, rx, ry, rz)长度方向一积分点 截面积分 横向剪切修正因子 可移动的中面方式s, t方向由向节点n3控制1.1.4 4隐式平

13、面应变单元：13号单元公式：1.1.5 5隐式轴对称单元：15号单元公式特殊条件下的4节点壳单元X-Y平面建模Y轴为15号单元的对称轴2个自由 度(dx,dy)建议约束自由度DOF (dz, rx, ry, rz)建议设置NIP=4非常快，很好的收 敛特性17.4 LS-DYNA隐式求解材料：LS-DYNA开发的材料主要应用于显式分析，并不全部支持隐式功能(但在 970 中大多数材料都支持隐式功能)不同的单元公式材料支持的程度也不一样，下面为960的材料支持程度：3DSolid单元1,2,3,4,5,6,7,9,12,13,18,20,24,26,30,49,57,63,91,92,103,1

14、04,106,115 Shell单元1,2,3,4,6,9,18,20,24,36,37,49,91,92,103,104,106,116,117,118Bearnl 元1,3,4,6,9,18,20,24,1002D Solid 单元1,2,3,4,5,6,7,9,12,13,18,20,24,26,57,6317.5 LS-DYNA隐式接触界面：LS-DYNA能用于隐式求解的接触界面如下：SUIU ACLT(),SURI ACT；NODES。一SURF ACEON E_W AY_S U R卜 AC】JO_SU R.FACL |()RMING_SLIR1 AClL_TO_SURFACIi b

15、OR M1 NG_NODES_TO_SU KbAClt b()RMJN<jJJNk_VVAY_XLIRIACT£7(LSURFACk A UTOMAT tC_S URFACE_TO_SURFACEAllTOMATIC_NODI£S_TOURI ACEAUTOMATIC_ONL_WAY_SURI'ACb,TO_SURI ACLA UTOMAT lC_*NGLE_gURF ACETlhL)_SUR1 4Ch_TO_$URhACL_OH SLTT1 ED_NUDLS_TO_§ U R 卜 A 3二_。卜 F SET21J_A ITO MAT I C_SU

16、RIACL _TU_、U RI ACE所有的隐式接触除了 TIED方式外都使用罚函数方式节点发射逻辑在隐式分析中自动失效(可选卡片B中参数SNLOG=1)建议使用重定向法向矢量(ORIEN=1或2)自动接触类型在隐式求解中常失效，不建议使用其他边界条件、约束条件、载荷、输出文件等与显式分析大同小异，这里不再详述。17.6隐式求解非线性收敛问题：应用LS-DYNA进行隐式求解常遇到的问题就是收敛问题，如：到达平衡叠代步数一不满足位移和能量容限，放弃该子步发散一不平衡力增加(而不是减小)，重新构造刚度矩阵，继续叠代 能量爆 炸一能量相比初始值增加百万倍，放弃该子步 负的特征值一这是在使用线性方 程

17、求解器在转置刚度矩阵时出现的错 误，缺省状态下忽略，试着继续运算线性搜索子步尺寸为零一新的位移不降低不平衡力；在没有到达最后叠代容限情况下重新构造刚度矩阵，否则放弃该子步17.6.1 几种诊断工具：使用 NLPRINT 输出开关(NLPRINT=1 )或CTRL+CNLPRINT 得至诊 断信息。负的特征值开关可以使程序在遇到负的特征值时立即用更小的时间步去 计算。CTRL+CCONV可以采用力收敛容限 CTRL+CITER创建每次叠代的 D3ITER数据库文件17.6.2 解决收敛问题的步骤：决定终止的原因(检查错误文件)激活输出开关得到更多的信息(NLPRINT=1)在叠代过程中使用D3I

18、TER检查几何体变形情况仔细坚持关键字输入文件上述收敛问题的解决方法：17.6.3 到达平衡叠代步数错误信息：观察位移和能量规范的收敛过程（查看 D3HSP文件，或激活NLPRINT ）如 果发现收敛容限几乎要满足：一增加更多的叠代步数一稍微放宽收敛容限一使用CTRL+CCONV采用力收敛容限如果收敛过程非常缓慢一改为全牛顿叠代方法（ILIMIT=1 ）一检查是否有不正确的材料特性一增加时间步长如果收敛过程不稳定一降低接触刚度（10倍量级缩放因子）一降低时间步尺寸17.6.4 发散和能量爆炸错误信息：一些PART模型可能变形太快 使用CTRL+CITER创建每次叠代的 D3ITER数据库文件在叠代过程中观察模型的变形情况检查是否有过大的载荷 检查是否有不正确屈服应力或硬化特性 检查接触定义，是否有过大的渗透 降低 时间步尺寸17.6.5 负的特征值错误信息：出现负的特征值有四种可能性：1 .不正确的材料特性：负杨氏模量，负的应力应变斜率等