第九讲王慎平显式非线性动态分析北京怡格明思工程技术有限公司北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation1显式动力学方法北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation显式动力学过程概述• 显式动力学求解器与隐式求解器，比如ABAQUS/Standard，是互为补充的。

从用户的角度出发，隐式与显式方法显著的区别为： 显式方法需要小的时间增量。

• 只与模型的最高自然频率相关。

• 与载荷类型和载荷持续时间无关。

• 一般的，增量步的数量级为10,000到1,000,000个增量，但是每个增 量步内的计算费用相对较小。

北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation显式动力学过程概述• 应力波的传播 应力波传播的例子说明了显式 动力学方法的求解过程：没 有迭代，或求解线性方程组。

考虑应力波沿着三个杆单元传 播问题。

在时间增加的过程 中，研究杆的状态。

• 质量被集中到节点。

杆的初始构型，自由端有一个集中力P北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation显式动力学过程概述u1 = &&−u & P & & ⇒ u1 = u1dt ⇒ ε el1 = 1 ⇒ dε el1 = ε el1dt & && l M1∫∫⇒ ε el1 = ε 0 + dε el1 ⇒ σ el1 = Eε el1第一个增量步结束时的构型北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation显式动力学过程概述P − Fel1 && & & && ⇒ u1 = u1old + u1dt u1 = M1∫& ε el1 =&& u2 =Fel1 & && ⇒ u 2 = u 2 dt M2∫u 2 − u1 & & & ⇒ dε el1 = ε el1dt l ⇒ ε el1 = ε1 + dε el1∫⇒ σ el1 = Eε el1第二个增量步开始时杆的构型第三个增量步开始时杆的构型北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation1显式时间积分ABAQUS/Explicit应用中心差分方法对运动方程进行显示的时间积分，应用一个 增量步的动力学条件计算下一个增量步的动力学条件。

在增量步开始时，程序 求解动力学平衡方程，表示为用节点质量矩阵M乘以节点加速度等于节点的合 力（在所施加的外力P与单元内力I之间的差值）：&& Mu = P - I在当前增量步开始时（t时刻），计算加速度为：&& u |(t ) = (M ) −1 ⋅ ( P − I ) |(t )由于显式算法总是采用一个对角的、或者集中的质量矩阵，所以求解加速度并 不复杂，不必同时求解联立方程。

任何节点的加速度是完全取决于节点质量和 作用在节点上的合力，使得节点计算的成本非常低。

北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation1对加速度在时间上进行积分采用中心差分方法，在计算速度的变化时假定加速度为 常数。

应用这个速度的变化值加上前一个增量步中点的速度来确定当前增量步中点 的速度： ( ∆t | +∆t | )& u|速度对时间的积分并加上在增量步开始时的位移以确定增量步结束时的位移：∆t (t + ) 2& =u|∆t (t − ) 2+( t +∆t )(t )2∆t ) 2&& u |( t )& u |( t +∆t ) = u |(t ) +∆t |( t +∆t ) u |(t +这样，在增量步开始时提供了满足动力学平衡条件的加速度。

得到了加速度，在时 间上“显式地”前推速度和位移。

所谓“显式”是指在增量步结束时的状态仅依赖于该 增量步开始时的位移、速度和加速度。

这种方法精确地积分常值的加速度。

为了使 该方法产生精确的结果，时间增量必须相当小，这样在增量步中加速度几乎为常数。

由于时间增量步必须很小，一个典型的分析需要成千上万个增量步。

幸运的是，因 为不必同时求解联立方程组，所以每一个增量步的计算成本很低。

大部分的计算成 本消耗在单元的计算上，以此确定作用在节点上的单元内力。

单元的计算包括确定 单元应变和应用材料本构关系（单元刚度）确定单元应力，从而进一步地计算内力。

北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation1显式动力学求解过程总结： 1. 节点计算 a. 动力学平衡方程&& u (t ) = (M ) −1 ⋅ ( P(t ) − I ( t ) )b. 对时间显式积分& u(t +∆t ) 2& =u(t −∆t ) 2+(∆t(t +∆t ) + ∆t(t ) ) 2(t + ∆t ) 2&& ut& u (t +∆t ) = u (t ) + ∆t(t +∆t ) u2. 单元计算 & a. 根据应变速率 ε ，计算单元应变增量 dε b. 根据本构关系计算应力 σσ (t +∆t ) = f (σ ( t ) , d ε )c. 集成节点内力 I( t +∆t ) 3. 设置时间 t 为 t + ∆t ，返回到步骤1。

北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation88显式时间积分方法的优越性1. 显式方法特别适用于求解需要分成许多小的时间增量来达到高精度的 高速动力学时间，诸如冲击，碰撞，爆破问题等； 2. 接触问题和其他一些极度非连续事件在显式方法中很容易表达清楚并 且能够一个节点一个节点地求解而不需要迭代。

节点加速度能够用来 调整外力和内力在接触中的平衡； 3. 显式方法罪显著的特性是没有整体切线刚度矩阵，而这是隐式方法所 必须的。

因为模型的状态为显式求解，所以不需要迭代和收敛准则；北京怡格明思工程技术有限公司PDF 文件使用 "pdfFactory Pro" 试用版本创建 Innovating through simulation显式方法的条件稳定性基于时间段开始时刻t的模型状态，应用显式方法求解，模型的状态通过时间的增量发生变化。

状态能够发生变化而且要保留对问题的精确描述，一般的时间增量非常短。

如果时间增量比最大的时间增量长，此时间增量就是所谓超出了稳定极限。

超过稳定极限的可能后果就是数值不稳定，会导致解答不收敛。

一般不可能精确地确定稳定极限，而是采用保守的估计值。

稳定极限对可靠性和精确性有很大的影响，所以必须一致和保守地确定。

为了计算的有效性，ABAQUS/Explicit 选择尽可能接近而且不超过稳定极限的时间增量。

稳定极限的定义稳定极限是依据系统的最高频率（）来定义的。

无阻尼时稳定极限由下式定义：maxωmax2ω=∆stablet有阻尼时由下式定义：)1(22maxξξω−+=∆stabletξ是具有最高频率的模型的临界阻尼比(为了控制高频振动）系统的实际最高频率是基于复杂的一组相互作用的因素，要计算出确切的值是不大可行的。

代替的办法是应用一个有效、保守的简单估算。

不考虑模型整体，我们估算模型中每个单独构件的最高频率，它常常与扩展的模态有关。

可以观察到由一个个单元为基础确定的最高频率常常比有限元组合模型的最高频率要高。

基于一个个单元的估算，稳定极限可以用单元长度和材料波速重新定义：eLdCdestable CLt=∆因为没有明确怎么确定单元的长度，对于多数单元类型，例如，一个扭曲的四边形单元，上述方程只是实际一个个单元的稳定极限的估算。

作为近似值，可以来用最短的单元尺寸，但是结果估算并不一定总是保守的。

单元长度越短，稳定极限越小。

波速是材料的一个特性：ρECd=其中E是杨氏模量，ρ是密度。

材料的刚度越大，波速越高，结果是稳定极限越小。

密度越高，波速越低，结果是稳定极限越大。

我们对简单的稳定极限定义提供了一些直觉的理解。

稳定极限是扩展波通过由单元特征长度定义的距离的短暂时间。

如果我们知道最小的单元尺寸和材料的波速，我们就能估算稳定极限。

例如，如果最小单元尺寸是5mm，扩展波速是5000m／s，稳定时间增量就在1×10－6s量级。

ABAQUS／Explicit的自动时间增量ABAQUS／Explicit 在分析的过程中自动调整时间增量的大小，以使基于当前模型状态的稳定极限永不越界。

时间增量是自动的，不需用户干涉，甚至不需要建议初始的时间增量。

稳定极限是从数值模型得来的数学概念。

因为有限元程序包含所有的相关细节，所以能够确定出一个有效并且守的稳定极限。

ABAQUS／Explicit 容许用户在必要时可不考虑自动时间增量。

质量缩放控制时间增量因为质量密度影响稳定极限，在某些环境下，缩放质量密度能够潜在地提高分析的效率。

例如，因为许多模型的复杂的离散性，有些区域常常包含着控制稳定极限的非常小或者形状极差的单元。

这些控制单元常常数量很少并且可能存在于局部区域。

通过只增加这些控制单元的质量，稳定极限可以显著的增加，同时对模型整体动力学行为的影响是可以忽略的。

ABAQUS/Explicit 的自动质量缩放功能可以保证这些单元不影响稳定极限。

质量缩放采用了两种基本方法：直接定义一个缩放因子或者给那些质量需要缩放的单元逐个定义稳定的时间增量。