Abaqus中显示动力学分析.docx

ABAQUS中的显式动力学分析

摘要：ABAQUS的分析模块主要由两个求解器完成ABAQUS/Stand和ABAQUS/Explicit,前者主

要用于静态分析和简单的准静态以及接触分析。而后者用于高速碰撞、复杂接触以及材料的

失效和退化仿真分析。ABAQUS/Explicit进行动力学分析最重要的是求解算法、稳定时间极

限和能量平衡问题。它们直接关系到计算的时间和计算所占用的内存以及计算结果的准确性,本文通过理论分析和实际举例验证和分析了ABAQUS/Explicit理论在实际中的应用。

关键词：ABAQUS/Explicit、求解算法、稳定时间极限、能量平衡

张敬东，成焕波

（合肥工业大学机械与汽车工程学院，合肥230009）

0前言

TheExplicitDynamic分析方法起初专门用于对高速动力学的分析，主要原因是隐式求解方

法特别占用内存而且极其耗时。比如在爆炸分析时载荷在瞬间增大而且持续时间非常短，精

确的捕捉结构响应的应力波是极其重要的，而应力波与系统的最高阶固有频率和最小时间增

量步有直接的关系。TheExplicitDynamics对于结构的瞬态响应和随后的震荡接触尤其适

用，而对于结构的刚度在分析的过程中发生显著变化的复杂后屈曲处理以及高度非线性的准静态分析和材料的失效和退化，该方法也能以及其独特的算法和极小的时间增量步完成，并

且不存在隐式求解中经常遇到的收敛问题。

1求解算法

ABAQUS/Explicit在整个运算过程中采用中心差分的算法求解整体平衡方程，运用动态方程

在每个微小的增量步内计算出下一个增量步的动态状况。

假定节点的质量矩阵M，加速度

??

u，外力作用为P，单元的内力I，因此节点的平衡方程是：

M

??

u=PI1

根据方程1可以计算出节点在当前增量步（t）时的节点加速度为：

??

u|(t)=(M)-1(P-I)(t)

由于显示求解器所用的质量矩阵全部为对角矩阵，而无需求解方程组。节点的加速度完全由

质量矩阵和施加在节点上的外力决定，因此求解十分节省资源。

在极短的时间内假定加速度为常值就可以运用中心差分对加速度积分计算出节点的速度方

程：

?

u|（t+

?t

2

?

）=u|（t-

?t

2

(?t|(t??t)??t|(t))??

）+u(t)

2

2

再次运用中心差分对速度积分求解出位移方程：

?t?

u|（t+?t)=u|(t)+?t|(t+?t)

u|(t+)3

2

由1、2、3方程可知模型节点的位移、速度和加速度只与初始增量步的大小有关。因此

ABAQUS/Explicit采用极小的初始时间增量和很多的增量步完成中心差分的运算。

由于对于不同节点的网格类型和网格的形状，后处理时查看的结果可能是节点应力也可能是

单元应力

?

假定单元应变增量为d?，单元的应变率为?，单元应力为?

由材料力学中的应力应变关系可以推导出单元应力应变的关系：

?(t??t)?f(?(t),d?)4

注：上述方程中?t为初始时间增量步。

在ABAQUS/Stand中进行非线性分析时，每一个时间增量步都要求解大量的非对角矩阵平衡方

程组，及其耗时而且在涉及接触时最容易出现不收敛的现象，因而采用ABAQUS/Explicit进行动力学分析具有省时，节约内存和精确性较高的有点等。

2稳定时间极限

?t

stable

ABAQUS/Explicit的算法表明时间增量步?t的大小直接关系到计算的准确性，如果该时间大于最大的允许时间会导致数值不稳定和求解的无限大现象。因此在分析之前必须尽可能精确地估算出稳定时间极限（要求到数量级）。

?）决定，在无阻尼状态下稳定时间极限由系统最高阶频率（

max

2

?t?5

stable

?

max

有阻尼状态下

22

?t?(1????)6

stable

?

max

?—最高阶频率的临界阻尼系数

但如果要确定系统的最高阶频率必须对系统进行线性摄动分析提取多阶模态值，过程较为复

杂。但是系统的最高阶频率与每个单元的最高阶频率紧密相关。并且单元的最高阶频率略高于整体的最高阶频率。为此ABAQUS/Explicit采用更为简便的方法估算稳定极限时间—即计算每个单元的最高阶频率对应的稳定时间极限。

?t=

stable

L

e

C

d

7

Le—单元在各个尺寸方向上的最小值。

Cd—模型材料的波速（这里只适用于弹塑性材料和线弹性材料，超弹性完全不可压缩材料的

波速趋近于无穷）

由7式可以得出稳定时间极限与单元的最小尺寸成正比，与材料的波速成反比。因而网格越

细，所需的稳定时间极限越小。

但是不同材料的波速不同

Cd