显式求解方法和隐式求解方法对比

《显式求解方法和隐式求解方法对比》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、采用有限元方法开展结构的动力学分析最终归结为求解离散后的常微分方程组。在时域内求解该方程最常用的方法是直接积分法，而又根据求解过程中是否需要迭代求解线性方程组，将直接积分法分为隐式积分方法和显式积分方法两类。隐式积分法认为ｔ＋Δｔ时刻系统的状态不仅与ｔ时刻状态有关，且与ｔ＋Δｔ时刻某些量有关。因此隐式算法是根据tn及tn-1...时刻体系的物理量值建立关于以tn+1时刻物理量为未知量的线性方程组，通过求解方程组确定tn+1时刻的物理量（常用的方法有线性加速度法、常平均加速度法、Newmark方法、Wilson-θ法、Houbolt方法等）

2、。而显式积分法认为ｔ＋Δｔ时刻系统的状态仅与ｔ时刻状态有关可，因此可由tn及tn-1...时刻体系的物理量值直接外推tn+1时刻物理量值（如中心差分法），不需要求解线性方程组，实现了时间离散的解耦。两种算法的比较 ：（1）隐式算法 隐式算法基于虚功原理，要迭代计算。隐式算法在每一增量步内都需要对静态平衡方程进行迭代求解，并且每次迭代都需要求解大型的线性方程组，这一过程需要占用相当数量的计算资源、磁盘空间和内存。理论上在这个算法中的增量步可以很大，但是实际运算中上要受到接触以及摩擦等条件的限制。随着单元数目的增加，计算时间几乎呈平方次增加。

3、由于需要矩阵求逆以及精确积分，对内存要求很高。隐式算法的不利方面就是收敛问题不容易解决，且在开始起皱失稳时，在分叉点处刚度矩阵出现奇异。（2）显式算法显示算法基于动力学方程，无需迭代，包括动态显式和静态显式算法。动态显式算法采用动力学方程的中心差分格式，不用直接求解切线刚度，不需要进行平衡迭代，计算速度快，也不存在收敛控制问题。该算法需要的内存也比隐式算法要少。数值计算过程可以很容易地进行并行计算，程序编制也相对简单。它也有一些不利方面。显式算法要求质量矩阵为对角矩阵，而且只有在单元级计算尽可能少时速度优势才能发挥,因而往往采用减缩积分方

4、法，容易激发沙漏模式，影响应力和应变的计算精度。静态显式法基于率形式的平衡方程组与Euler前插公式，不需要迭代求解。由于平衡方程式仅在率形式上得到满足，所以得出的结果会慢慢偏离正确值。为了减少相关误差，必须每步使用很小的增量，通常一个仿真过程需要多达几千步。由于不需要迭代，所以这种方法稳定性好，但效率低。以下是ＡＢＡＱＵＳ软件中显式和隐式常采用的方法：（摘自帮助文档）结构系统的通用运动学方程为：（1）中心差分法（显式）假定0，，，…，时刻的节点位移，速度与加速度均为已知，现求解时刻的结构响应。中心差分法对加速度，速度的导数采用中心差分代

5、替，即为：(2)将（2）式代入（1）式后整理得到（3）式（3）中分别称为有效质量矩阵，有效载荷矢量。，，，为结构载荷，质量，阻尼，刚度矩阵。求解线性方程组（3），即可获得时刻的节点位移向量，将代回几何方程与物理方程，可得时刻的单元应力和应变。中心差分法在求解瞬时的位移时，只需时刻以前的状态变量和，然后计算出有效质量矩阵，有效载荷矢量，即可求出，故称此解法为显式算法。中心差分法，在开始计算时，需要仔细处理。=0时，要计算，需要知道的值。因此应该有一个起始技术，因而该算法不是自动起步的。由于，，是已知的，由=0时的（2）式可知：中心差分法中时

6、间步长的选择涉及两个方面的约束：数值算法的稳定性和计算时间。中心差分法的实质是用差分代替微分，并且对位移和加速度的导数采用线性外插，这限制了的取值不可过大，否则结果可能失真过大。可以证明：中心差分法是条件稳定的。即当时间步长必须小于由该问题求解方程性质所决定的一个时间步长的临界值。LS-DYNA中，采用“变时间步长法”，即每一时刻的步长由当前结构的稳定性条件来控制。具体算法为：计算每一个单元的极限时间步长，=1,2…，取为下一个时刻的时间步长。各种单元的计算方法如下。（1）1D杆，梁单元其中为时间步长因子，系统默认为0.9。为杆，梁单元的

7、长度。材料声速。（2）2D板，壳单元其中为时间步长因子，系统默认为0.9。为壳单元的最小单元边长度。为材料的声速。（3）3D单元其中和为无量纲常数，默认=1.5，=0.06.为单元等效长度，为单元体积，为单元最大侧面积。为材料声速。时间步因子可由用户设置，减小相当于减少时间步长。设置时间步长因子的关键字为*CONTROL_TIMESTEP，控制参数为TSSFAC。另外，质量缩放可以人为控制时间步长。即调整单元密度来改变时间步长。以壳单元为例说明质量缩放改变时间步长。（假如=1）由得到Newmark法（隐式）Newmark假定在时间间隔内，

8、加速度线性变化，即采用如下的加速度，速度公式：(4)式中，为按积分的精度和稳定性要求可以调整的参数。根据（4）式可给出和用，，表示的表达式，代入（1）式中整理得到（5）其中称之为有效刚度矩阵和