abaqus系列教程-13ABAQUSExplicit准静态分析.doc

13 ABAQUS/Explicit 准静态分析显式求解方法是一种真正的动态求解过程，它的最初发展是为了模拟高速冲击问 题，在这类问题的求解中惯性发挥了主导性作用当求解动力平衡的状态时，非平衡力 以应力波的形式在相邻的单元之间传播由于最小稳定时间增量一般地是非常小的值， 所以大多少问题需要大量的时间增量步在求解准静态问题上，显式求解方法已经证明是有价值的，另外 ABAQUS/Explicit 在求解某些类型的静态问题方面比 ABAQUS/Standard 更容易在求解复杂的接触问题 时，显式过程相对于隐式过程的一个优势是更加容易此外，当模型成为很大时，显式 过程比隐式过程需要较少的系统资源关于隐式与显式过程的详细比较请参见第 2.4 节“隐式和显式过程的比较”将显式动态过程应用于准静态问题需要一些特殊的考虑根据定义，由于一个静态 求解是一个长时间的求解过程， 所以在其固有的时间尺度上分析模拟常常在计算上是不 切合实际的，它将需要大量的小的时间增量因此，为了获得较经济的解答，必须采取 一些方式来加速问题的模拟但是带来的问题是随着问题的加速，静态平衡的状态卷入 了动态平衡的状态，在这里惯性力成为更加起主导作用的力。

目标是在保持惯性力的影 响不显著的前提下用最短的时间进行模拟准静态 ( Quasi-static )分析也可以在 ABAQUS/Standard 中进行 当惯性力可以忽略 时，在 ABAQUS/Standard 中的准静态应力分析用来模拟含时间相关材料响应(蠕变、 膨胀、粘弹性和双层粘塑性)的线性或非线性问题关于在 ABAQUS/Standard 中准静 态分析的更多信息， 请参阅 ABAQUS 分析用户手册 (ABAQUS Analysis User 's Manual ) 的第 6.2.5 节“ Quasi-static analysis ”13.1 显式动态问题类比为了使你能够更直观地理解在缓慢、准静态加载情况和快速加载情况之间的区别， 我们应用图 13-1 来类比说明快速情况图13-1缓慢和快速加载情况的类比图中显示了两个载满了乘客的电梯在缓慢的情况下，门打开后你步入电梯为了 腾出空间，邻近门口的人慢慢地推他身边的人，这些被推的人再去推他身边的人，如此 继续下去这种扰动在电梯中传播，直到靠近墙边的人表示他们无法移动为止一系列 的波在电梯中传播，直到每个人都到达了一个新的平衡位置。

如果你稍稍加快速度，你 会比前面更用力地推动你身边的人， 但是最终每个人都会停留在与缓慢的情况下相同的在快速情况下，门打开后你以很高的速度冲入电梯，电梯里的人没有时间挪动位置 来重新安排他们自己以便容纳你你将会直接地撞伤在门口的两个人，而其他人则没有 受到影响对于准静态分析，实际的道理是同样的分析的速度经常可以提高许多而不会严重 地降低准静态求解的质量；缓慢情况下和有一些加速情况下的的最终结果几乎是一致 的但是，如果分析的速度增加到一个点，使得惯性影响占主导地位时，解答就会趋向 于局部化，而且结果与准静态的结果是有一定区别的13.2 加载速率一个物理过程所占用的实际时间称其为它的固有时间( n ature time )对于一个准静态过程在固有时间中进行分析， 我们一般地有把握假设将得到准确的静态结果 毕竟,如果实际事件真实地发生在其固有时间尺度内，并在结束时其速度为零，那么动态分析 应该能够得到这样的事实，即分析实际上已经达到了稳态你可以提高加载速率使相同 的物理事件在较短的时间内发生，只要解答保持与真实的静态解答几乎相同，而且动态 效果保持是不明显的13.2.1 光滑幅值曲线对于准确和高效的准静态分析，要求施加的载荷尽可能地光滑。

突然、急促的运动 会产生应力波，它将导致振荡或不准确的结果以可能最光滑的方式施加载荷要求加速 度从一个增量步到下一个增量步只能改变一个小量如果加速度是光滑的，随其变化的 速度和位移也是光滑的ABAQUS有一条简单、固定的光滑步骤(smooth step)幅值曲线，它自动地创建一 条光滑的载荷幅值当你定义一个光滑步骤幅值曲线时， ABAQUS自动地用曲线连接每一组数据对，该曲线的一阶和二阶导数是光滑的，在每一组数据点上，它的斜率都为 零由于这些一阶和二阶导数都是光滑的，你可以采用位移加载，应用一条光滑步骤幅 值曲线，只用初始的和最终的数据点，而且中间的运动将是光滑的使用这种载荷幅值 允许你进行准静态分析而不会产生由于加载速率不连续引起的波动 一条光滑步骤幅值曲线的例子，如图13-2所示图13-2采用光滑步骤幅值曲线的幅值定义13.2.2 结构问题在静态分析中，结构的最低模态通常控制着结构的响应如果已知最低模态的频率 和相应的周期，你可以估计出得到适当的静态响应所需要的时间为了说明如何确定适当的加载速率，考虑在汽车门上的一根梁被一个刚性圆环从侧面侵入的变形, 所示实际的实验是准静态的如图13-3杨氏模量= 200 GPa固定边界泊松比=0.3屈服应力 = 250 MPa 硬化模量 =20 MPa 密度=7800 Kg/m 3 壳厚=3 mm固定边界圆柱梁长=1 m图13-3刚性圆环与梁的碰撞采用不同的加载速率，梁的响应变化很大。

以一个极高的碰撞速度为 400 m/s，在梁中的变形是高度局部化的，如图 13-4所示为了得到一个更好的准静态解答，考虑图13-4碰撞速度为400 m/s最低阶的模态最低阶模态的频率大约为 250 Hz,它对应于4 ms的周期应用在ABAQUS/Standard中的特征频率提取过程可以容易地计算自然频率 为了使梁在4 ms内发生所希望的0.2m的变形，圆环的速度为50 m/s虽然50 m/s似乎仍然像是一个高速碰撞速度， 而惯性 力相对于整个结构的刚度已经成为次要的了，如图 13-5所示，变形形状显示了很好的准静态响应图13-5碰撞速度为50 m/s虽然整个结构的响应显示了我们所希望的准静态结果， 但通常理想的是将加载时间增加到最低阶模态的周期的 10倍以确保解答是真正的准静态为了更进一步地改进结果，刚环的速度可能会逐渐增大，例如应用一条光滑步骤幅值曲线，从而减缓初始的冲13.2.3 金属成型问题为了获得低成本的求解过程，人为地提高成型问题的速度是必要的，但是，我们能 够把速度提高多少仍可以获得可接受的静态解答呢？如果薄金属板毛坯的变形对应于其最低阶模态的变形形状，可以应用最低阶结构模态的时间周期来指导成型的速度。

然 而在成型过程中，刚性的冲模和冲头能够以如此的方式约束冲压，使坯件的变形可能与 结构的模态无关在这种情况下，一般性的建议是限制冲头的速度小于 1%的薄金属板的波速对于典型的成型过程，冲头速度是在 1 m/s的量级上，而钢的波速大约为 5000m/s因此根据这个建议，一个 50的因数为冲头提高速度的上限为了确定一个可接受的冲压速度， 建议的方法包括以各种变化的冲压速度运行一系 13-#列的分析，这些速度在 3m/s 至 50m/s 的范围内由于求解的时间与冲压的速度成反比， 运行分析是以冲压速度从最快到最慢的顺序进行检查分析的结果，并感受变形形状、 应力和应变是如何随冲压速度而改变的冲压速度过高的一些表现是与实际不符的、局部化的拉伸与变薄，以及对起皱的抑止如果你从一个冲压速度开始，例如 50 m/s,并从某处减速，在某点上从一个冲压速度到下一个冲压速度的解答将成为相似的，这说明 解答开始收敛于一个准静态的解答当惯性的影响成为不明显时,在模拟结果之间的区 别也是不明显的随着人为地增加加载速率,以逐渐和平滑的方式施加载荷成为越来越重要的方式 例如,最简单的冲压加载方式是在整个成型过程中施加一个定常的速度。

在分析开始时,如此加载会对薄金属板坯引起突然的冲击载荷, 在坯件中传播应力波并可能产生不希望 的结果当加载速率增加时,任何冲击载荷对结果的影响成为更加明显的应用光滑步 骤幅值曲线,使冲压速度从零逐渐增加可以使这些不利的影响最小化回弹回弹经常是成型分析的一个重要部分, 因为回弹分析决定了卸载后部件的最终 形状 尽管 ABAQUS/Explicit 十分适合于成型模拟, 对回弹分析却遇到某些特殊的 困难在 ABAQUS/Explicit 中进行回弹模拟最主要的问题是需要大量的时间来获得 稳态的结果特别是必须非常小心地卸载,并且必须引入阻尼以使得求解的时间比 较合理幸运的是,在 ABAQUS/Explicit 和 ABAQUS/Standard 之间的紧密联系允 许一种更有效的方法由于回 弹 过 程不涉 及 接 触 , 而 且一 般只包 括 中 度 的 非 线性, 所以 ABAQUS/Standard 可以求解回弹问题, 并且比 ABAQUS/Explicit 求解得更快 因此, 对于回弹分析更偏爱的方法是将完整的成型模型从 ABAQUS/Explicit输入（import）到 ABAQUS/Standard 中进行。

在这本指南中不讨论输入功能13.3 质量放大质量放大（mass scaling）可以在不需要人为提高加载速率的情况下降低运算的成本 对于含有率相关材料或率相关阻尼（如减震器）的问题,质量放大是惟一能够节省求解 时间的选择在这种模拟中,不要选择提高加载速度,因为材料的应变率会与加载速率同比例增加当模型的参数随应变率变化时，人为地提高加载速率会人为地改变了分析 的过程稳定时间增量与材料密度之间的关系如下面的方程所示如在第 923节“稳定极Le限的定义”中所讨论的，模型的稳定极限是所有单元的最小稳定时间增量它可以表示 成为.：t式中，Le是特征单元长度，Cd是材料的膨胀波速线弹性材料在泊松比为零时的膨胀波速给出为这里，t是材料密度根据上面的公式，人为地将材料密度 增加因数f2倍，则波速就会降低因数 f倍,从而稳定时间增量将提高因数 f倍注意到当全局的稳定极限增加时，进行同样的分析所需要的增量步就会减少，而这正是质量放大的目的但是，放大质量对惯性效果与人 为地提高加载速率恰好具有相同的影响 因此，过度地质量放大，正像过度地加载速率，可能导致错误的结果为了确定一个可接受的质量放大因数，所建议的方法类似于确定 一个可接受的加载速率放大因数。

两种方法的唯一区别是与质量放大相关的加速因子是质量放大因数的平方根，而与加载速率放大相关的加速因子是与加载速率放大因数成正 比例如，一个为100倍的质量放大因数恰好对应于 10倍的加载速率因数通过使用固定的或可变的质量放大，可以有多种方法来实现质量放大编程质量放大的定义也可以随着分析步而改变，允许有很大的灵活性详细的内容请参阅ABAQUAS 分析用户手册第 7.15.1 节“ Mass scaling ”14.4 能量平衡评估模拟是否产生了正确的准静态响应， 最具有普遍意义的方式是研究模型中的各种能量下面是在 ABAQUS/Explicit中的能量平衡方程：E| ■ Ev ' eke ' efd - Ew - Etotai - constant式中，Ei是内能（包括弹性和塑性应变能）， Ev是粘性耗散吸收的能量， Eke是动能,Efd是摩擦耗散吸收的能量， Ew是外力所做的功，Etotai是在系统中的总能量为了应用一个简单的例子来说明能量平衡，考虑如图 13-6所示的一个单轴拉伸实验准静态实验的能量历史将显示在图 13-7中如果模拟是准静态的，那么外力所做的功是几乎等于系统内部的能量。