显式和隐式动力分析比较

显式分析:用上一步的结果和当前步的结果计算下一步的计算结果。有条件收敛,要求时间步较小。通常做动力分析用这种方法。

隐式分析:用当前步结果和下一步未知结果反复迭代下一步结果,必须通过迭代得到。无条件收敛。是一种能量平衡的结果。通常做静力分析用这种方法。

显式算法:

显式算法最大优点是有较好的稳定性。

动态显式算法采用动力学方程的一些差分格式(如广泛使用的中心差分法、线性加速度法等),不用直接求解切线刚度,不需要进行平衡迭代,计算速度快,时间步长只要取的足够小,一般不存在收敛性问题。 显式算法不需要迭代,也不需要组集总刚,因此需要的内存也比隐式算法要少。并且数值计算过程可以很容易地进行并行计算,程序编制也相对简单。但显式算法要求质量矩阵为对角矩阵,而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法,容易激发沙漏模式,影响应力和应变的计算精度。"

静态显式法基于率形式的平衡方程组与Euler向前差分法,不需要迭代求解。由于平衡方程式仅在率形式上得到满足,所以得出的结果会慢慢偏离正确值。为了减少相关误差,必须每步使用很小的增量。这个方法目前应用比较少。

隐式算法 :

在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这以过程需要占用相当数量的计算资源、磁盘空间和内存。该算法中的增量步可以比较大,至少可以比显式算法大得多,但是实际运算中上要受到迭代次数及非线性程度的限制,需要取一个合理值。通常在一个时间增量下,隐式算法在弹塑性后期此时刚度比较奇异了,难以收敛,也就需要采用二分法把当前的分析步折半继续解方程,并且不断二分下去直到方程求解完成,所以隐式算法时在前面计算速度还不错,但越往后越慢,时间也是不可估量的,可能在地震波20s就算不下去了等等。但是隐式算法有个好处就是只要算了步数其结果一般一定是正确的,不会像显式那样虽然算完了,其结果可能是发散的。

隐式算法在求解运动方程时是当前步结果除了更上一步有关外还跟当前步有关,因此不能解耦大型方程组,需要组集刚度,用得多是newmark法和wilson-theta法,其中的beta、alpha和theta值的确定可改变方法的稳定性、计算效率等。隐式求解法的最大优点是它具有无条件稳定性,即时间步长可以任意大。

计算时间:

上面大致说了下,不过一般经验表明:

通常认为两者有此区别:

从上面大致可以出,有些认识只是表象而已,其中(1)是明显的错误,因为两种方法都可以用于动力学和静力学问题,只是一种选择而已,(2)也是错误的,阻尼的选择不是算法本身,而是求解需求。(03)是算法的本质,也决定了(04)(05)(06)(07)(08),这是算法的特点所在,04~08并非算法本身的特点而是每一步求解方法的特征。同理(09)是算法的特点,决定了(10)(11)。

过(03)(09)可以得到两种方法的计算特点,显式算法是每一步求解为矩阵乘法,时间步选择为条件稳定;隐式算法是每一步求解为线性方程组求解,时间步选择为无条件稳定。

附加说明:

{1)求解线性静力学问题,虽然求解线性方程组,但是没有时步的关系,所以不应将其看作隐式算法。而是静力学隐式方程求解问题。

(2)求解非线性静力学问题,虽然求解过程需要迭代,或者是增量法,但是没有明显的时步问题,所以不应将其看作隐式算法。仍归结为隐式方程求解问题

3)静态松弛法,可以认为是将动力学问题看作静力学问题来解决,每一步达到静力平衡,需要数值阻尼。

4)动态松弛法,可以认为是将静力学问题或者动力学问题,分为时步动力学问题,采用向后时步迭代的思想计算。对于解决静力学问题时,需要人工阻尼。