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

1、13ABAQUS/Explicit准静态分析 显式求解方法是一种真正的动态求解过程，它的最初发展是为了模拟高速冲击问 题，在这类问题的求解中惯性发挥了主导性作用。当求解动力平衡的状态时，非平衡力 以应力波的形式在相邻的单元之间传播。由于最小稳定时间增量一般地是非常小的值， 所以大多少问题需要大量的时间增量步。 在求解准静态问题上，显式求解方法已经证明是有价值的，另外ABAQUS/Explicit 在求解某些类型的静态问题方面比ABAQUS/Standard更容易。在求解复杂的接触问题 时，显式过程相对于隐式过程的一个优势是更加容易。此外，当模型成为很大时，显式 过程比隐式过程需要较少的系统资源

2、。关于隐式与显式过程的详细比较请参见第2.4 节 “隐式和显式过程的比较”。 将显式动态过程应用于准静态问题需要一些特殊的考虑。根据定义，由于一个静态 求解是一个长时间的求解过程，所以在其固有的时间尺度上分析模拟常常在计算上是不 切合实际的，它将需要大量的小的时间增量。因此，为了获得较经济的解答，必须采取 一些方式来加速问题的模拟。但是带来的问题是随着问题的加速，静态平衡的状态卷入 了动态平衡的状态，在这里惯性力成为更加起主导作用的力。目标是在保持惯性力的影 响不显著的前提下用最短的时间进行模拟。 准静态 （Quasi-static）分析也可以在ABAQUS/Standard 中进行。 当惯性

3、力可以忽略 时，在ABAQUS/Standard中的准静态应力分析用来模拟含时间相关材料响应（蠕变、 膨胀、粘弹性和双层粘塑性）的线性或非线性问题。关于在ABAQUS/Standard中准静 态分析的更多信息，请参阅 ABAQUS 分析用户手册 （ABAQUS Analysis User s Manual） 的第 6.2.5 节“Quasi-static analysis”。 13.1显式动态问题类比 为了使你能够更直观地理解在缓慢、准静态加载情况和快速加载情况之间的区别， 我们应用图13-1 来类比说明。 13-1 缓慢情况快速情况 图 13-1缓慢和快速加载情况的类比 图中显示了两个载满了

4、乘客的电梯。在缓慢的情况下，门打开后你步入电梯。为了 腾出空间，邻近门口的人慢慢地推他身边的人，这些被推的人再去推他身边的人，如此 继续下去。这种扰动在电梯中传播，直到靠近墙边的人表示他们无法移动为止。一系列 的波在电梯中传播，直到每个人都到达了一个新的平衡位置。如果你稍稍加快速度，你 会比前面更用力地推动你身边的人，但是最终每个人都会停留在与缓慢的情况下相同的 位置。 在快速情况下，门打开后你以很高的速度冲入电梯，电梯里的人没有时间挪动位置 来重新安排他们自己以便容纳你。你将会直接地撞伤在门口的两个人，而其他人则没有 受到影响。 对于准静态分析，实际的道理是同样的。分析的速度经常可以提高许多

5、而不会严重 地降低准静态求解的质量；缓慢情况下和有一些加速情况下的的最终结果几乎是一致 的。但是，如果分析的速度增加到一个点，使得惯性影响占主导地位时，解答就会趋向 于局部化，而且结果与准静态的结果是有一定区别的。 13.2加载速率 一个物理过程所占用的实际时间称其为它的固有时间（nature time）。对于一个准 静态过程在固有时间中进行分析，我们一般地有把握假设将得到准确的静态结果。毕竟， 如果实际事件真实地发生在其固有时间尺度内，并在结束时其速度为零，那么动态分析 应该能够得到这样的事实，即分析实际上已经达到了稳态。你可以提高加载速率使相同 的物理事件在较短的时间内发生，只要解答保持与

6、真实的静态解答几乎相同，而且动态 效果保持是不明显的。 13-2 13.2.1光滑幅值曲线 对于准确和高效的准静态分析，要求施加的载荷尽可能地光滑。突然、急促的运动 会产生应力波，它将导致振荡或不准确的结果。以可能最光滑的方式施加载荷要求加速 度从一个增量步到下一个增量步只能改变一个小量。如果加速度是光滑的，随其变化的 速度和位移也是光滑的。 ABAQUS 有一条简单、固定的光滑步骤（smooth step）幅值曲线，它自动地创建一 条光滑的载荷幅值。当你定义一个光滑步骤幅值曲线时，ABAQUS自动地用曲线连接 每一组数据对，该曲线的一阶和二阶导数是光滑的，在每一组数据点上，它的斜率都为 零。

7、由于这些一阶和二阶导数都是光滑的，你可以采用位移加载，应用一条光滑步骤幅 值曲线，只用初始的和最终的数据点，而且中间的运动将是光滑的。使用这种载荷幅值 允许你进行准静态分析而不会产生由于加载速率不连续引起的波动。一条光滑步骤幅值 曲线的例子，如图13-2所示。 值 幅 时间 图 13-2采用光滑步骤幅值曲线的幅值定义 13.2.2结构问题 在静态分析中，结构的最低模态通常控制着结构的响应。如果已知最低模态的频率 和相应的周期，你可以估计出得到适当的静态响应所需要的时间。为了说明如何确定适 13-3 当的加载速率， 考虑在汽车门上的一根梁被一个刚性圆环从侧面侵入的变形，如图 13-3 所示。实际

8、的实验是准静态的。 杨氏模量= 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 13-4 最低阶模态的频率大约为250 Hz，它对应于4 ms 的周期。应用在 ABAQUS/Standard 中的特征频率

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

10、成型问题的速度是必要的，但是，我们能 够把速度提高多少仍可以获得可接受的静态解答呢？如果薄金属板毛坯的变形对应于 其最低阶模态的变形形状，可以应用最低阶结构模态的时间周期来指导成型的速度。然 而在成型过程中，刚性的冲模和冲头能够以如此的方式约束冲压，使坯件的变形可能与 结构的模态无关。在这种情况下，一般性的建议是限制冲头的速度小于1%的薄金属板 的波速。对于典型的成型过程，冲头速度是在1m/s 的量级上，而钢的波速大约为5000 m/s。因此根据这个建议，一个50 的因数为冲头提高速度的上限。 为了确定一个可接受的冲压速度，建议的方法包括以各种变化的冲压速度运行一系 13-5 列的分析，这些速

11、度在3m/s至 50m/s的范围内。由于求解的时间与冲压的速度成反比， 运行分析是以冲压速度从最快到最慢的顺序进行。检查分析的结果，并感受变形形状、 应力和应变是如何随冲压速度而改变的。冲压速度过高的一些表现是与实际不符的、局 部化的拉伸与变薄，以及对起皱的抑止。如果你从一个冲压速度开始，例如50m/s，并 从某处减速，在某点上从一个冲压速度到下一个冲压速度的解答将成为相似的，这说明 解答开始收敛于一个准静态的解答。当惯性的影响成为不明显时，在模拟结果之间的区 别也是不明显的。 随着人为地增加加载速率，以逐渐和平滑的方式施加载荷成为越来越重要的方式。 例如，最简单的冲压加载方式是在整个成型过程

12、中施加一个定常的速度。在分析开始时， 如此加载会对薄金属板坯引起突然的冲击载荷，在坯件中传播应力波并可能产生不希望 的结果。当加载速率增加时，任何冲击载荷对结果的影响成为更加明显的。应用光滑步 骤幅值曲线，使冲压速度从零逐渐增加可以使这些不利的影响最小化。 回弹 回弹经常是成型分析的一个重要部分，因为回弹分析决定了卸载后部件的最终 形状。 尽管 ABAQUS/Explicit十分适合于成型模拟，对回弹分析却遇到某些特殊的 困难。在 ABAQUS/Explicit中进行回弹模拟最主要的问题是需要大量的时间来获得 稳态的结果。特别是必须非常小心地卸载，并且必须引入阻尼以使得求解的时间比 较合理。幸

13、运的是，在ABAQUS/Explicit和 ABAQUS/Standard 之间的紧密联系允 许一种更有效的方法。 由 于 回 弹 过 程 不 涉 及 接 触 ， 而 且 一 般 只 包 括 中 度 的 非 线 性 ， 所 以 ABAQUS/Standard 可以求解回弹问题， 并且比 ABAQUS/Explicit求解得更快。 因此， 对于回弹分析更偏爱的方法是将完整的成型模型从ABAQUS/Explicit输入（import） 到 ABAQUS/Standard 中进行。 在这本指南中不讨论输入功能。 13.3质量放大 质量放大（ massscaling）可以在不需要人为提高加载速率的情况

14、下降低运算的成本。 对于含有率相关材料或率相关阻尼（如减震器）的问题，质量放大是惟一能够节省求解 时间的选择。在这种模拟中，不要选择提高加载速度，因为材料的应变率会与加载速率 13-6 同比例增加。当模型的参数随应变率变化时，人为地提高加载速率会人为地改变了分析 的过程。 稳定时间增量与材料密度之间的关系如下面的方程所示。如在第9.2.3 节“稳定极 限的定义”中所讨论的，模型的稳定极限是所有单元的最小稳定时间增量。它可以表示 成为 t e L c d 式中，L e是特征单元 长度， e是特征单元长度， c 是材料的膨胀波速。线弹性材料在泊松比为零时的膨胀 d 波速给出为 c d E 这里，是

15、材料密度。 根据上面的公式，人为地将材料密度增加因数 2 f 倍，则波速就会降低因数f 倍， 从而稳定时间增量将提高因数f倍。注意到当全局的稳定极限增加时，进行同样的分析 所需要的增量步就会减少，而这正是质量放大的目的。但是，放大质量对惯性效果与人 为地提高加载速率恰好具有相同的影响。因此，过度地质量放大，正像过度地加载速率， 可能导致错误的结果。为了确定一个可接受的质量放大因数，所建议的方法类似于确定 一个可接受的加载速率放大因数。两种方法的唯一区别是与质量放大相关的加速因子是 质量放大因数的平方根，而与加载速率放大相关的加速因子是与加载速率放大因数成正 比。例如，一个为100倍的质量放大因

16、数恰好对应于10 倍的加载速率因数。 通过使用固定的或可变的质量放大，可以有多种方法来实现质量放大编程。质量放 大的定义也可以随着分析步而改变，允许有很大的灵活性。详细的内容请参阅 ABAQUAS分析用户手册第7.15.1 节“Mass scaling”。 13-7 14.4能量平衡 评估模拟是否产生了正确的准静态响应，最具有普遍意义的方式是研究模型中的各 种能量。下面是在ABAQUS/Explicit中的能量平衡方程： EEEEEE IVKEFDWtotal constant 式中， EI是内能（包括弹性和塑性应变能），EV是粘性耗散吸收的能量，EKE是动能， EFD是摩擦耗散吸收的能量，E

17、W是外力所做的功，Etotal是在系统中的总能量。 为了应用一个简单的例子来说明能量平衡，考虑如图13-6 所示的一个单轴拉伸实 验。 准静态实验的能量历史将显示在图13-7 中。如果模拟是准静态的，那么外力所做 的功是几乎等于系统内部的能量。除非有粘弹性材料、离散的减震器、或者使用了材料 阻尼，否则粘性耗散能量一般地是很小的。由于在模型中材料的速度很小，所以在准静 态过程中，我们已经确定惯性力可以忽略不计。由这两个条件可以推论，动能也是很小 的。作为一般性的规律，在大多数过程中，变形材料的动能将不会超过它的内能的一个 小的比例（典型的为5%到 10%）。 量 能 时间 图 13-6单轴拉伸实

18、验图 13-7准静态拉伸实验的能量历史 当比较能量时， 请注意 ABAQUS/Explicit报告的是整体的能量平衡，它包括了任何 含有质量的刚体的动能。由于当评价结果时我们只对变形体感兴趣，当评价能量平衡时 13-8 我们应在Etotal中扣除刚体的动能。 例如，如果你正在模拟一个采用滚动刚体模具的传输问题，刚体的动能可能占据模 型整个动能的很大部分。在这种情况下，你必须扣除与刚体运动有关的动能，然而才可 能做出与内能有意义的比较。 13.5例题：ABAQUS/Explicit凹槽成型 在这个例子中，你将应用ABAQUS/Explicit求解第12 章“接触”中那个凹槽成型 的问题。然后，你

19、将比较分别来自ABAQUS/Standard 和 ABAQUS/Explicit分析的结果。 你将修改由ABAQUS/Standard分析所创建的模型，这样才能在ABAQUS/Explicit 中运行它。这些修改包括在材料模型中增加密度，改变单元库，并改变分析步。为了获 得正确的准静态响应，在运行 ABAQUS/Explicit分析前， 你将应用在ABAQUS/Standard 的频率提取过程来确定所需要的计算时间。 13.5.1前处理应用 ABAQUS/Explicit重新运算模型 对于这个模拟，应用ABAQUS/CAE修改模型。在本手册的在线文档第A.12 节 “Forming a cha

20、nnel”提供了输入文件。当通过ABAQUS/CAE运行这个输入文件时， 将创建关于该问题的完整的分析模型。根据下面给出的指导如果你遇到困难，或者如果 你希望检查你的工作，则可以运行这个输入文件。在附录A“Example Files ”中，给出 了如何提取和运行输入文件的指导。 在开始前，打开关于凹槽成型例题的模型数据库文件，它被创建在第12.6 节 “ABAQUS/Standardexample ：forming a channel”。 确定一个合适的分析步时间 对于一个准静态过程，在第 13.2 节“加载速率” 中讨论了确定合适的分析步时 间的过程。如果我们知道了坯件的最低阶固有频率，即基

21、（fundamental）频，我们 就可 以确定分析步时间的一个大致的下限。一种获得这个信息的方法是在 ABAQUS/Standard 中运行频率分析。在这个成型分析中，冲压对坯件产生的变形类 似于它的最低阶模态。因此，如果你想模拟整个结构而并非局部的变形，选择第一 个成型阶段的时间是大于或等于坯件最低阶模态的周期是十分重要的。 13-9 运行一个固有频率提取过程： 1将已存在的模型复制成为一个新的模型，命名为 Frequency，并对 Fequency 模型进行如下全面的修改：在频率提取分析中，你将用一个单独的频率提取分 析步取代现在所有的分析步。此外，你将删除所有的刚性工具和接触相互作用；

22、 它们与确定毛坯的基频无关。 2在 Property 模块中，为Steel材料模型增加一个7800 的密度。 3在 Assembly 模块中，删除冲模、冲头和夹具部件的实体。对于频率分析并不需 要这些刚体部件。 提示 ：你可以从工具箱中采用Delete工具删除这些部件。 4进入 Step模块，用一个单独的频率提取分析步替代现存的所有分析步。 a. 在StepManager （分析步管理器）中，删除分析步RemoveRight Constraint、HolderForce、 EstablishContact II和 Move Runch 。 b. 选择分析步Establish Contact I

23、，并点击Replace 。 c. 在 Repalce Step（替换分析步）对话框中，从LinearParturbation过程 列表中选择Frequency，键入分析步描述为Frequency modes；选择 Lanczos 特征值选项，并要求五个特征值。重新命名分析步为Extract Frequencies。 d. 取消 DOFMonitor（自由度监视器）选项。 注意 ：由于频率提取分析步是一个线性扰动过程，将忽略材料的非线性性质。 在这个分析中，坯件的左端约束沿x-方向的位移和绕法线的转动；但是，没有 约束沿y-方向的位移。因此，提取的第一阶模态将是刚体模态。对于在 ABAQUS/E

24、xplicit中的准静态分析，第二阶模态的频率将确定合适的时间段。 5在 Interaction 模块，删除所有的接触相互作用。 6进入Load 模块，在BCManager （边界条件管理器）中检查在Extract Frequencies分析步中的边界条件。除了边界条件名称CenterBC以外，删 除所有的边界条件。将这个留下的采用了对称边界条件的毛坯约束施加到左端。 7. 在创建和提交作业前，如果有必要则重新剖分网格。 8. 进入 Job 模块，创建一个作业，命名为Forming-Frequency，采用如下的作 13-10 业描述： Channel forming - frequency

25、analysis。提交作业进行分 析，并监控求解过程。 9当分析完成时，进入Visualization模块，并打开由这个作业创建的输出数据库 文件。从主菜单栏中，选择Plot -Deformed Shape；或者应用在工具箱中的 工具。 绘制出一阶屈曲模态的模型变形形状。进一步绘出毛坯的二阶模态，将未变形 的模型形状叠加在模型变形图上。 频率分析表明坯件有一个140 Hz 的基频，对应的周期为0.00714 s。图 13-8 显 示了第二阶模态的位移形状。对于成型分析，我们现在知道最短的分析步时间 为 0.00714 s。 图 13-8由 ABAQUS/Standard 频率分析的毛坯二阶模态

26、 创建 ABAQUS/Explicit成型分析 成型过程的目标是采用0.03m 的冲头位移准静态地成型一个凹槽。在选择准静 态分析的加载速率时，建议你在开始时用较快的加载速率，并根据需要减小加载速 率，更快地收敛到一个准静态解答。然而，如果你希望在你的第一次分析尝试中就 增加能够得到准静态结果的可能性，你应当考虑分析步时间是比相应的基频缓慢10 到 50 倍的因数。在这个分析中，对于成型分析步，你将从0.007 s的时间开始。这 是基于在ABAQUS/Standard 中进行的频率分析，它显示出毛坯具有140Hz 的基频， 对应于 0.00714s 的时间周期。这个时间周期对应于4.3 m/s

27、的常数冲头速度。你将 仔细地检查动能和内能的结果，以检验结果中并没有包含显著的动态影响。 13-11 将 Standard模型复制成一个新模型，命名为Explicit。如果必要，通过 从位于工具栏下方的Model （模型）列表中选择Explicit模型作为当前的模型。 使所有接下来的模型改变成为Explicit模型。 在 ABAQUS/Standard分析中，在冲头和坯件之间模拟一个初始的缝隙以便于 接触计算。在ABAQUS/Explicit分析中则不需要采取这种预防措施。因此，在 Assembly 模块中，沿U2 方向平移冲头-0.001 m。在警告对话框中出现的关于相对 和绝对约束中，点击

28、Yes。 在毛坯夹具上施加一个集中力，为了计算夹具的动态反应，必须在刚性体的参 考点上赋予一个点质量。夹具的实际质量是不重要的；而重要的是它的质量必须与 毛坯的质量（ 0.78kg）具有同一个数量级，以使在接触计算中的振荡最小化。选择 数值为 0.1 kg 的点质量。在Property模块中，创建一个点的截面定义，命名为 Pointmass。在 Edit Section对话框的 Inertial Properties域中，键入 0.1点质量的值。在参考点RigidRefHolder应用这个截面定义。此外，编辑 Steel材料定义来包括7800 kg/m 3的质量密度。 进入 Step 模块。你

29、需要为ABAQUS/Explicit分析创建两个分析步。在第一个 分析步中施加夹具力；在第二个分析步中施加冲头压下力。除了命名为Establish Contact I 的分析步之外，删除所有其他的分析步，并用一个单一的显式动态分析步 替换这个分析步。键入分析步描述为Applyholder force，并指定0.0001s 的分析步时间。这个时间对于施加夹具载荷是适合的，因为它是足够长以避免了动 态效果，而且又足够短以防止了对整个作业运行时间的明显冲击。将分析步重新命 名为 Holderforce。创建第二个显式动态分析步，命名为Displace punch， 分析步的时间为0.007s，键入

30、Apply punch stroke作为分析步的描述。 为了帮助确定分析是如何接近于准静态假设，研究各种能量的历史是非常有用 的。特别有用的是比较动能和内部应变能。能量历史默认地写入了输出数据库文件。 在这个金属成型分析的第一次尝试中，对于施加的夹具力和冲头压力，你将应 用具有默认的光滑参数的表格形式的幅值曲线。进入Load 模块，为施加的夹具力 创建一个名为Ramp1 的表格形式的幅值曲线。在表13-1 中输入幅值数据。为冲头 压力定义第二个表格形式的幅值曲线，命名为 Ramp2。在表 13-2 中输入幅值数据。 13-12 表 13-1Ramp1 和 Smooth1的递增幅值数据 时间（

31、s）幅值 0.00.0 0.00011.0 表 13-2 Ramp2和 Sooth2的递增幅值数据 时间（ s）幅值 0.00.0 0.0071.0 在 LoadManager（载荷管理器）中，在命名为Holderforce的分析步中 创建一个集中力，命名为RefHolderForce，在施加的点上指定RefHolder和 一个沿着CF2 方向大小为 -440000的力。对于这个载荷， 改变幅值定义为Ramp1。 在BoudaryConditionManager（ 边 界条件 管理 器 ） 中 ， 删 除命名为 MidLeftBC和 MidRightBC的边界条件。编辑RefDieBC边界条件

32、，这样在 Holderforce分析步中沿着U2 方向的约束为零，不改变其他方向的约束。对于 RefHolderBC边界条件， 解除沿着U2 方向的约束， 而其他方向的约束保持不变。 在 DisplacePunch 分析步中，改变位移边界条件RefPunchBC，使沿着U2 方 向的位移为 -0.03m。对于这个边界条件，应用幅值曲线Ramp2。 监视自由度的值 在这个模型中，你将在整个分析步中监视冲头的参考节点的竖向位移（自由度 2）。在ABAQUS/Standard成型分析中，由于已经设置了DOFMonitor监视 RefPunch的竖向位移，所以你无需做出任何改变。 创建网格和定义作业

33、在网格 Mesh 模块中，将用于剖分坯件网格的单元族改变为Explicit ，并指定 13-13 增强沙漏控制，并剖分坯件网格。因为已经将工具模拟成了解析刚性表面，因此无 需将它们剖分网格。 在 Job 模块中创建一个作业，命名为Forming-1，给予作业如下的描述： Channel forming -attempt 1。 在运行成型分析前，你可能希望知道该分析将需要多少个增量步，进而了解该 分析需要多少计算机时间。你可以通过运行数据检查（ data check）分析来获得关于 初始稳定时间增量的近似值，或者你可以应用在第13.3节“质量放大” 中的关系式 进行估计。在这个例题中，从一个增量

34、步到下一个增量步的稳定时间增量不会有太 大的变化，因此知道了稳定时间增量，你可以确定完成成型阶段的分析需要多少个 增量步。一旦分析开始，你就能够知道每一个增量步需要多少CPU 时间，进而知 道整个分析需要多少CPU 时间。 应用在 13.3 节“质量放大” 中表述的关系式，关于这个分析的稳定时间增量近 -7s。因此，对于 0.007 s的分析步时间，成型阶段需要大约185,000 个增 似为 1 10 量步。 将模型保存到模型数据库文件中，并提交作业进行分析。监视求解过程；改正 任何检测到的模拟错误，并调查任何警告信息的原因。完成整个分析可能需要运行 10 分钟或更长的时间。 一旦分析开始运行

35、，在另一个视图窗中会显示出你选择来监视（冲头的竖向位 移）的自由度值的X-Y 曲线图。从主菜单栏中，选择Viewport -JobMonitor: Forming-1，在分析运行的整个时间中跟踪沿着2-方向冲头位移的发展进程。 评价结果的策略 在查看我们最关心的结果之前，诸如应力和变形形状，我们需要确定结果是否 是准静态的。一个好的方法是比较动能与内能的历史。在金属成型分析中，大部分 的内能是由于塑性变形产生的。在这个模型中，坯件是动能的主要因素（忽略夹具 的运动，没有与冲头和模具相关的质量）。为了确定是否已经获得了一个可接受的 准静态解答，坯件的动能应该小于其内能的几个百分点。对于更高的精确

36、度，特别 地是对回弹应力感兴趣时，动能应该是更低的。这个方法是非常有用的，因为它应 用于所有类型的金属成型过程，而且不需要任何直观地理解在模型中的应力；许多 成型过程可能是过于复杂，以至于不允许对结果有一个直观的判断。 虽然是衡量准静态分析的良好和重要的证明，仅凭动能与内能的比值还不足以 13-14 确任解的质量。 你还必须对这两种能量进行独立地评估，以确定它们是否是合理的。 当需要准确的回弹应力结果时，这一部分的评估是更增加了重要性，因为一个高度 精确的回弹应力解答是高度地依赖于准确的塑性结果。即使动能是非常小的量，如 果它包含了高度的振荡，则模型也会经历显著的塑性。一般说来，我们希望光滑加

37、 载以产生光滑的结果；如果加载是光滑的，但是能量的结果是振荡的，则结果可能 是不合适的。由于一个能量的比值无法显示这种行为，所以你也必须研究动能本身 的历史以观察是否是光滑的还是振荡的。 如果动能不能显示出准静态的行为，在某些节点上观察速度的历史可能是有用 的，以帮助理解在各个区域中模型的行为。这种速度历史可以表明在模型的哪些区 域是振荡的，并产生大量的动能。 评估结果 进入 Visulization 模块，并打开由这个作业 （Forming-1.odb）创建的输出数据库。 绘制动能和内能。 创建能量历史的曲线： 1从主菜单栏中，选择Plot -History Output。 显示出整个模型的

38、伪应变能历史曲线。 2从主菜单栏中，选择Result -History Output。 显示出 History Output对话框。 3从变量的列表中，选择Kinetic energy:ALLKE for Whole Model。 4点击 Plot 创建一条ALLKE 的历史曲线。 显示出整个模型的动能历史曲线（见图13-9）。 13-15 图 13-9成型分析的动能历史，尝试1 5类似地，创建模型内能的历史曲线，ALLIE （见图 13-10）。 图 13-10成型分析的内能历史，尝试1 13-16 显示在图13-9 中的动能历史发生显著地振荡。另外，动能的历史与坯件的成型没 有明确的关系，

39、这表明这个分析是不适合的。在这个分析中，冲头的速度保持为常数， 而主要地依赖于坯件运动的动能却远非是恒定值。 在除了开始阶段以外的整个分析步中，比较图 13-9和图 13-10 表明动能是内能的一 个很小的百分数（小于1%）。即使对于这种严重的加载情况，还是满足了动能必须相 对地小于内能地准则。 尽管模型的动能只是内能的一个小的分数，它还是有一定的振荡。所以，我们应该 以某种方式改变模拟以获得更平滑的解答。 14.5.2成型分析尝试2 即使实际上冲头是以几乎接近于常值的速度运动，第一次模拟尝试的结果表明理想 的方式是采用不同的幅值曲线以允许坯件更光滑地加速。当考虑应用什么类型的加载幅 值时，记

40、住在准静态分析的所有方面，光滑性是重要的。最偏爱的方法是尽可能光滑地 移动冲头，在理想的时间内移动理想的距离。 应用一种光滑地施加的冲头力和一段光滑地施加的冲头距离，我们现在将分析成型 阶段；我们将与前面获得的结果进行比较。关于光滑步骤幅值曲线的解释，请阅读第 13.2.1 节“光滑幅值曲线”。 在 Load 模块中，定义一条光滑步骤幅值曲线，命名为Smooth1 。输入在表13-1 中给出的幅值数据。创建第二条光滑步骤幅值曲线，命名为Smooth2 ，应用在表13-2 中给出的幅值数据。在Holderforce分析步中，修改RefHolderForce载荷，使 它采用Smooth1的幅值。在

41、Displace punch分析步中，修改位移边界条件 RefPunchBC ，使它采用 Smooth2的幅值。 通过设置在分析步开始时的幅值为0.0和在 分析步结束时的幅值为1.0，ABAQUS/Explicit创建了一个幅值定义，它的一阶和二阶 导数都是光滑的。因此，应用一条光滑步骤幅值曲线对位移进行控制，也使我们确信了 其速度和加速度是光滑的。 在 Job模块中，创建一个作业， 命名为 Forming-2，给予作业如下的描述： Channel forming -attempt 2。 将模型保存到模型数据库文件中，并提交作业进行分析。监视求解过程；改正任何 检测到的模拟错误，并调查任何警告

42、信息的原因。完成整个分析可能需要运行10 分钟 或更长的时间。 13-17 评估第二次尝试的结果 动能的结果如图13-11 所示。动能的响应是明显地与坯件的成型相关：在第二 个分析步的中间阶段出现了动能的峰值，它对应于冲头速度最大的时刻。因此，动 能是适当的和合理的。 图 13-11成型分析的动能历史，尝试2 关于第二次尝试的内能历史如图13-12 所示，显示了从零上升到最终值的光滑 增长。再次看出，动能与内能的比值是相当小的，并显示出是可接受的。图13-13 比较了在两次成型尝试中的内能。 13-18 图 13-12成型分析的内能历史，尝试2 图 13-13关于成型分析的两次尝试的内能比较

43、13-19 13.5.3两次成型尝试的讨论 我们评价结果可接受性的初始原则是动能与内能相比必须是小量。我们发现即使对 于最严重的情况，尝试1，这个条件似乎是仍然得到了满足。增加光滑步骤幅值曲线帮 助减小了在动能中的振荡，得到了令人满意的准静态响应。 附加的要求动能和内能的历史必须是适当的和合理的是非常有用的和必 要的，但是它们也增加了评价结果的主观性。在一般更为复杂的成型过程中，强调这些 要求可能是很困难的，因为这些要求的提出需要对成型过程的行为的一些直观考虑。 成型分析的结果 我们现在已经满意了关于成型分析的准静态解答是合适的，我们可以研究感兴 趣的某些其它结果。图13-14 显示了应用AB

44、AQUS/Standard和 ABAQUS/Explicit 得到的在坯件中Mises 应力的比较。 图 13-14在 ABAQUS/Standard （左）和 ABAQUS/Explicit （右） 凹槽成型分析中Mises应力的等值线图 从图中显示在ABAQUS/Standard 和 ABAQUS/Explicit分析中的应力峰值的差 13-20 别在 1%以内，并且在坯件中整个应力的等值线图是非常类似的。为了进一步检验 准静态分析结果的有效性，你应该从两个分析中比较等效塑性应变的结果和最终变 形的形状。 图 13-15 显示了在坯件中等效塑性应变的等值线图，而图13-16 显示了由两个

45、分析预测的最终变形形状的覆盖图。对于ABAQUS/Standard和 ABAQUS/Explicit 的分析，等效塑性应变的结果彼此相差在5%以内。另外，最终变形形状的比较显 示出显式准静态分析的结果与ABAQUS/Standard 静态分析的结果吻合得极好。 图 13-15在 ABAQUS/Standard （左）和 ABAQUS/Explicit （右） 凹槽成型分析的PEEQ的等值线图 13-21 图 13-16在 ABAQUS/Standard 和 ABAQUS/Explicit成型分析的凹槽最终变形形状 你也应该比较由ABAQUS/Standard和 ABAQUS/Explicit分

46、析预测的稳态冲头 压 力 。 如 图13-17 可 见 ， 由ABAQUS/Explicit预 测 的 稳 态 冲 头 压 力 值 比 由 ABAQUS/Standard 预测的值大约高12%。 在 ABAQUS/Standard 和 ABAQUS/Explicit 结果之间的这个差别主要是源于两个因素。首先，ABAQUS/Explicit规则化了材料 数据。 其次， 在两个分析软件中摩擦效果的处理稍有区别；ABAQUS/Standard 使用 罚函数摩擦，而ABAQUS/Explicit使用动力学摩擦。 13-22 图 13-17对于 ABAQUS/Standard 和 ABAQUS/Exp

47、licit的稳态冲头压力比较 从这些比较中，可以明显看出ABAQUS/Standard 和 ABAQUS/Explicit都有能 力处理诸如本例问题的困难接触分析。然而，在 ABAQUS/Explicit中运行这类分析 有某些优势： 与 ABAQUS/Standard 相比，ABAQUS/Explicit能够更容易地处理复杂 的接触条件和采用较少的分析步和边界条件进行计算。特别地，ABAQUS/Standard 分析需要五个分析步和附加的边界条件以确保正确的边界条件和防止刚体运动。在 ABAQUS/Explicit中完成同样的分析只需要两个分析步和无需附加边界条件。然而， 当选择 ABAQUS

48、/Explicit进行准静态分析时，你必须明确在一个合适的加载速率下 你可能需要进行迭代。在确定加载速率时，建议你开始时采用较快的加载速率，并 根据需要减小加载速率。这可以帮助优化对分析进行求解的时间。 13-23 13.5.4加速分析的方法 现在我们已经获得了一个可接受的成型分析的解答，我们可以尝试采用更短的计算 机时间来获得类似的可接受的结果。因为采用显式动态标准的成型问题的实际时间是过 大的，所以大部分成型分析都需要过多的计算机时间以至于无法按照它们自己的物理时 间尺度进行运算；若使分析在一个可接受的计算机时间范围内运行，常常需要对分析做 出改变以减少计算机成本。有两种节省分析成本的方法

49、： 1人为地增加冲头的速度，从而在一个更短的分析步时间内发生同样的成型过程。 这种方法称为加载速率放大（load rate scaling ）。 2人为地增加单元的质量密度，从而增大稳定时间极限，允许分析采用较少的增 量步。这种方法称为质量放大（mass scaling）。 这两种方法等效地做相同的事情，除非模型具有率相关材料或者阻尼。 确定可接受的质量放大 第 13.2节“加载速率”和第13.2.3 节“金属成型问题”讨论了如何确定可接受 的加载速率或质量的放大因子以加速准静态分析的时间尺度。目标是在保持惯性力 不显著的前提下以最短的时间模拟过程。求解的时间加快多少是有界的，而且还要 能够得

50、到一个有意义的准静态解答。 如在第 13.2 节“加载速率” 中讨论的那样， 我们可以应用同样的方法以确定一 个合适的质量放大因子，如我们已经应用以确定一个合适的加载速率放大因子的方 法。在两种方法之间的区别是加载速率放大因子f与质量放大因子f 2的效果相同。 最初，我们假设分析步的时间为坯件的基频周期的阶数时会产生适当的准静态结 果。通过研究模型的能量和其他的结果，我们相信这些结果是可以接受的。这项技 术产生了大约4.3 m/s 的冲头速度。我们现在将接受采用质量放大的求解时间，并 将结果与我们没有质量放大求解的结果进行比较，以确定由质量放大得到的结果是 否可以接受。我们假设这种放大仅可能降