显式非线性动态分析.doc

Note ABAQUS 2016有限元分析从入门到精通 Note ABAQUS 2016有限元分析从入门到精通 Note第7章 显式非线性动态分析 Note 第7章 显式非线性动态分析 ·PAGE 180· ·PAGE 179· 第 章 显式非线性动态分析 本章将重点介绍使用ABAQUS进行显示分析的步骤和方法，使读者熟练掌握应用ABAQUS进行显示分析的解题过程。 显式动态程序对于求解各种类型的非线性固体和结构力学问题是一种非常实用有效的工具，这是对隐式求解器，如ABAQUS/Standard的一个补充。从用户的角度来看，显式与隐式方法的区别在于，显式方法需要很小的时间增量步，它仅取决于模型的最高固有频率，而与载荷的类型和持续的时间无关，而隐式方法对时间增量步的大小没有内在的限制，增量的大小通常依赖于精度和收敛情况。 了解显式分析的基础理论 熟悉钢球与钢板撞击过程的分析 任务驱动项目案例 7.1 ABAQUS/Explicit适用的问题类型 在研究探讨显式动态程序怎样工作之前，需要了解ABAQUS/Explicit适合于求解哪类问题。 1．高速动力学事件 最初发展显式动力学方法的目的是分析那些用隐式方法（如ABAQUS/Standard）分析可能极端费时的高速动力学事件。作为这种类型模拟的例子，在第7章中分析了一块钢板在瞬态冲击载荷下的响应。因为短时间迅速施加的巨大载荷，结构的响应在这段时间内变化得非常快。对于捕获动力响应，精确地跟踪板内的应力波是非常重要的。因为应力波与系统的最高阶频率相互关联，所以为了得到精确解答需要许多足够小的时间增量。 2．复杂的接触问题 相对于应用隐式方法建立接触条件的公式，应用显式动力学方法要容易得多。结论是ABAQUS/ Explicit能够比较容易地解决包括许多独立物体相互作用的复杂接触问题。ABAQUS/Explicit特别适合应用于分析受冲击载荷并随后在结构内部发生复杂相互接触作用的结构的瞬间动态响应问题。在接触中所展示的电路板跌落试验就是这类问题的一个例子。在这个例子中，一块插入在泡沫封装中的电路板从1m的高度跌落到地板上。这个问题包括封装与地板之间的冲击，以及在电路板和封装之间的接触条件的迅速变化。 3．复杂的后屈曲问题 ABAQUS/Explicit能够比较容易地解决不稳定的后屈曲（Post?buckling）问题。在这类问题中，随着载荷的施加，结构的刚度会发生剧烈的变化。在后屈曲响应中常常涉及接触相互作用的影响。 4．高度非线性的准静态问题 由于各种不同的原因，ABAQUS/Explicit往往能够有效地解决某些在本质上是静态的问题。准静态（quasi-static）过程模拟问题（包括复杂的接触，如锻造、滚压和薄板成形等过程）一般属于这一类型的问题。薄板成形问题通常包含非常大的膜变形、褶皱和复杂的摩擦接触条件。块体成形问题的特征包括大扭曲、瞬间变形以及与模具之间的相互接触。 5．材料退化和失效 在隐式分析程序中，材料的退化（Degradation）和失效（Failure）常常导致严重的收敛困难，但是ABAQUS/Explicit能够很好地模拟这类材料。材料退化中的一个例子是混凝土开裂的模型，其拉伸裂缝导致了材料的刚度变为负值。金属的延性失效模型是一个材料失效的例子，其材料刚度能够退化并且一直降低到零，在这段时间中，单元从模型中被完全除掉。 这些类型分析的每一个问题都有可能包含温度和热传导的影响。 7.2 动力学显式有限元方法 本节讲述了ABAQUS/Explicit求解器的算法，对隐式和显式时间积分做了比较，并探讨了显式方法的优越性。 7.2.1 显式时间积分 ABAQUS/Explicit应用中心差分方法对运动方程进行显式的时间积分，应用一个增量步的动力学条件计算下一个增量步的动力学条件。在增量步开始时，程序求解动力学平衡方程，表示为用节点质量矩阵M乘以节点加速度ü等于节点的合（在所施加的外力P与单元内力I之间的差值）： （7-1） 在当前增量步开始时（t时刻），计算加速度为 （7-2） 由于显式算法总是采用一个对角的或者集中的质量矩阵，所以求解加速度并不复杂，不必同时求解联立方程。任何节点的加速度完全取决于节点质量和作用在节点上的合力，使得节点计算的成本非常低。 对加速度在时间上进行积分采用中心差分方法，在计算速度的变化时假定加速度为常数。应用这个速度的变化值加上前一个增量步中点的速度来确定当前增量步中点的速度： （7-3） 速度对时间的积分加上在增量步开始时的位