显式和隐式求解方法是有限元中最为关键的知识,对于初学者,可能对其内在的含义还是理解不够,只是记着“大变形用显式,线性小变形用隐式”这样的一般性结论,若是能对显式和隐式有更深层的理解,对于有限元内在的求解方式将会掌握更好。
ABAQUS和LS-DYNA都可以进行显式和隐式求解,不同的是,ABAQUS更擅长隐式求解,而LS-DYNA显式求解更强,至于强在何处,就是另外的话了,这里主要是介绍下显式和隐式的含义。
1、含义
隐式求解,即implicit method,在ABAQUS中,Standard模块主要进行隐式求解的计算,在分析步中进行设置;LS-DYNA则使用关键字*CONTROL_IMPLICIT_GENERAL进行显隐式的设置。
图 1:ABAQUS隐式设置
图 2:LS-DYNA隐式设置
隐式求解的特点是利用迭代的方法求解下个增量步的未知量,即对于一个问题,隐式是将其看作一个整体,进行矩阵的计算,迭代方法一般为Newton-Rapson法,这种可以比作“鲸吞”,因此,隐式求解没有条件稳定,任何大小的时间增量皆可让结果在一定范围内,但由于采用的是迭代的方法,因此有计算收敛性问题。
不同于隐式,显式求解(Explicit method),ABAQUS中由Explicit模块求解,LS-DYNA中默认采用的即为显式求解方法。显式求解利用中央差分法,借助多个时间增量完成模拟。显式只关注前一时刻的状态,它每一步的求解都是基于前一步的结果,通过预先设置的时间增量来递推后面的结果,因此说,显式相当于将一个问题分成很多块,然后一步步去计算,类似“蚕食”。你会发现,上面叙述的时间增量就是很关键的东西了,如果时间增量很大,显式求解偏离真实结果就会很大,一般来说增量必须小于一临界值,才会提供可靠的结果,显式有条件稳定,但无计算收敛问题。后面会针对显式的时间增量进行进一步讨论,因为它在求解过程中还涉及到一个质量缩放的概念。
图 3:显式和隐式求解增量及求解过程 [江*]
2、区别
总结一下显式和隐式的区别,如下图所示。
图 4:显式和隐式求解区别
3、显式质量缩放:
上面显隐式各自的特点会发现:显式求解的高精度往往依赖于极小的时间增量,这类似于数学中的积分,自变量划分的越小,积分结果越精确,这自然产生了一个问题,想得到可用的结果,什么样的时间增量是可行的呢?
其实,ABAQUS和LS-DYNA中都能根据求解模型选择自动时间增量,软件已经帮你计算时间增量,直接求解即可,可是当模型很大时,或者模型某些位置的网格很小,软件得到的自动时间增量会非常小,显式求解的速度会非常慢,那有什么办法解决呢?这里就涉及到显式求解时间增量的由来了。
在显式求解中,求解稳定时间增量和系统的最高自然频率有关,系统没有阻尼时,稳定极限定义为:
有阻尼时稳定极限定义为:
是最高频率模态的临界阻尼部分。可以证明,以逐个单元为基础确定的最高单元频率总是高于有限元组合模型的最高频率。因此,基于单个单元的估算,稳定极限可以用单元长度L和材料波速c重新定义:
因此,单元长度L越小,材料刚度E越大,密度越小,稳定极限越小。通常软件估算得到的稳定极限就是根据上面的公式而来。
知道了自动时间增量原理,如果模型太大或者某些网格很小时,就可以人为的增大时间增量来加快求解速度,这就叫做质量缩放(质量加速)。
4、质量缩放方法
LS-DYNA中的质量缩放通过关键字*CONTROL_TIMESTEP里参数DT2MS实现,设置DT2MS为正值和负值来控制是对整个模型进行质量缩放还是对某些特定的单元进行缩放
LS DYNA质量缩放 Xue Hao,公众号:逐梦之xueLS DYNA质量缩放
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删