接触有限元仿真之初,就常听人提起,ANSYS隐式求解更好,ABAQUS显式求解更好,理了理两种方式到底有什么区别。
1.动力学问题
在讨论隐式与显式动力学之前,先讨论一下动力学问题和静力学问题。在求解静力学问题时没有时间的概念,即使在载荷步控制力有Time这个选项,但是这个Time的含义更多的是载荷步。Time前后的求解过程相互没有影响。动力学问题的特点是施加到结构上的外载荷的大小和方向可能随着时间的变化而发生变化,使结构产生速度和加速度。在用有限元求解静力学问题时主要是求矩阵方程组的问题,如下所示,只需要考虑结构的刚度矩阵。当考虑动力响应之后,除了考虑刚度矩阵还需要考虑质量矩阵和阻尼矩阵,分别计入速度和加速度的影响,也即阻尼力和惯性力的作用。
2.动力学问题的求解
隐私求解和显式求解这两个概念相信大家都听的很多,ANSYS的显示动力学LS-DYNA和Abaqus Explict 均是采用显式动力学求解动力问题。
所谓显式和隐式,是指求解方法的不同,即数学上的出发点不一样。并不是说显式只能求动力学问题,隐式只能求静力学问题,只是求解策略不通。隐式求解和时间无关,采用的是牛顿迭代法(线性问题就直接求解线性代数方程组),因此存在一个迭代收敛问题,不收敛就的不到结果。由于两者解题的出发点,所以一般来说显式用于求解和时间相关的动力学问题。隐式用来求解和时间无关的静力学问题。
显式求解方法在ANSYS显式动力学采用的是中心差分法,这个方法在求解温度场时用的较多。基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
对于中心差分法,如果知道最初的两个位移,即可一步一步求出之后的位移速度和加速度,从而得到系统的位移速度和加速度时程曲线。
求解动力学问题的另一种方法是Newmark法,采用这种方法求解属于隐式求解。