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

1、13 ABAQUS/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缓慢和快速加载情况的类比图中显示了两个载满了乘

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

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

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

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

8、力 = 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似乎仍然像是一个高速

9、碰撞速度， 而惯性 力相对于整个结构的刚度已经成为次要的了，如图13-5所示，变形形状显示了很好的准静态响应。图13-5碰撞速度为50 m/s虽然整个结构的响应显示了我们所希望的准静态结果，但通常理想的是将加载时间增加到最低阶模态的周期的10倍以确保解答是真正的准静态。为了更进一步地改进结果，刚环的速度可能会逐渐增大，例如应用一条光滑步骤幅值曲线，从而减缓初始的冲13.2.3 金属成型问题为了获得低成本的求解过程，人为地提高成型问题的速度是必要的，但是，我们能 够把速度提高多少仍可以获得可接受的静态解答呢？如果薄金属板毛坯的变形对应于其最低阶模态的变形形状，可以应用最低阶结构模态的时间周期来指

10、导成型的速度。然 而在成型过程中，刚性的冲模和冲头能够以如此的方式约束冲压，使坯件的变形可能与 结构的模态无关。在这种情况下，一般性的建议是限制冲头的速度小于1%的薄金属板的波速。对于典型的成型过程，冲头速度是在1 m/s的量级上，而钢的波速大约为5000m/s。因此根据这个建议，一个50的因数为冲头提高速度的上限。为了确定一个可接受的冲压速度，建议的方法包括以各种变化的冲压速度运行一系 13-5列的分析，这些速度在 3m/s 至 50m/s 的范围内。由于求解的时间与冲压的速度成反比， 运行分析是以冲压速度从最快到最慢的顺序进行。检查分析的结果，并感受变形形状、 应力和应变是如何随冲压速度而

11、改变的。冲压速度过高的一些表现是与实际不符的、局部化的拉伸与变薄，以及对起皱的抑止。如果你从一个冲压速度开始，例如50 m/s,并从某处减速，在某点上从一个冲压速度到下一个冲压速度的解答将成为相似的，这说明 解答开始收敛于一个准静态的解答。当惯性的影响成为不明显时,在模拟结果之间的区 别也是不明显的。随着人为地增加加载速率,以逐渐和平滑的方式施加载荷成为越来越重要的方式。 例如,最简单的冲压加载方式是在整个成型过程中施加一个定常的速度。在分析开始时,如此加载会对薄金属板坯引起突然的冲击载荷, 在坯件中传播应力波并可能产生不希望 的结果。当加载速率增加时,任何冲击载荷对结果的影响成为更加明显的。

12、应用光滑步 骤幅值曲线,使冲压速度从零逐渐增加可以使这些不利的影响最小化。回弹回弹经常是成型分析的一个重要部分, 因为回弹分析决定了卸载后部件的最终 形状。 尽管 ABAQUS/Explicit 十分适合于成型模拟, 对回弹分析却遇到某些特殊的 困难。在 ABAQUS/Explicit 中进行回弹模拟最主要的问题是需要大量的时间来获得 稳态的结果。特别是必须非常小心地卸载,并且必须引入阻尼以使得求解的时间比 较合理。幸运的是,在 ABAQUS/Explicit 和 ABAQUS/Standard 之间的紧密联系允 许一种更有效的方法。由于回 弹 过 程不涉 及 接 触 , 而 且一 般只包 括

13、 中 度 的 非 线性, 所以 ABAQUS/Standard 可以求解回弹问题, 并且比 ABAQUS/Explicit 求解得更快。 因此, 对于回弹分析更偏爱的方法是将完整的成型模型从ABAQUS/Explicit输入（import）到 ABAQUS/Standard 中进行。在这本指南中不讨论输入功能。13.3 质量放大质量放大（mass scaling）可以在不需要人为提高加载速率的情况下降低运算的成本。 对于含有率相关材料或率相关阻尼（如减震器）的问题,质量放大是惟一能够节省求解 时间的选择。在这种模拟中,不要选择提高加载速度,因为材料的应变率会与加载速率同比例增加。当模型的参数随

14、应变率变化时，人为地提高加载速率会人为地改变了分析 的过程。稳定时间增量与材料密度之间的关系如下面的方程所示。如在第923节“稳定极Le限的定义”中所讨论的，模型的稳定极限是所有单元的最小稳定时间增量。它可以表示 成为.：t式中，Le是特征单元长度，Cd是材料的膨胀波速。线弹性材料在泊松比为零时的膨胀波速给出为这里，t是材料密度。根据上面的公式，人为地将材料密度增加因数f2倍，则波速就会降低因数f倍,从而稳定时间增量将提高因数f倍。注意到当全局的稳定极限增加时，进行同样的分析所需要的增量步就会减少，而这正是质量放大的目的。但是，放大质量对惯性效果与人 为地提高加载速率恰好具有相同的影响。因此，

15、过度地质量放大，正像过度地加载速率，可能导致错误的结果。为了确定一个可接受的质量放大因数，所建议的方法类似于确定 一个可接受的加载速率放大因数。两种方法的唯一区别是与质量放大相关的加速因子是质量放大因数的平方根，而与加载速率放大相关的加速因子是与加载速率放大因数成正 比。例如，一个为100倍的质量放大因数恰好对应于10倍的加载速率因数。通过使用固定的或可变的质量放大，可以有多种方法来实现质量放大编程。质量放大的定义也可以随着分析步而改变，允许有很大的灵活性。详细的内容请参阅ABAQUAS 分析用户手册第 7.15.1 节“ Mass scaling ”。14.4 能量平衡评估模拟是否产生了正确

16、的准静态响应，最具有普遍意义的方式是研究模型中的各种能量。下面是在 ABAQUS/Explicit中的能量平衡方程：E| Ev eke efd - Ew - Etotai - constant式中，Ei是内能（包括弹性和塑性应变能），Ev是粘性耗散吸收的能量，Eke是动能,Efd是摩擦耗散吸收的能量，Ew是外力所做的功，Etotai是在系统中的总能量。为了应用一个简单的例子来说明能量平衡，考虑如图13-6所示的一个单轴拉伸实验。准静态实验的能量历史将显示在图13-7中。如果模拟是准静态的，那么外力所做的功是几乎等于系统内部的能量。除非有粘弹性材料、离散的减震器、或者使用了材料 阻尼，否则粘性耗

17、散能量一般地是很小的。由于在模型中材料的速度很小，所以在准静 态过程中，我们已经确定惯性力可以忽略不计。由这两个条件可以推论，动能也是很小 的。作为一般性的规律，在大多数过程中，变形材料的动能将不会超过它的内能的一个小的比例（典型的为 5%到10%）。rLh图13-6单轴拉伸实验图13-7准静态拉伸实验的能量历史FD当比较能量时，请注意ABAQUS/Explicit报告的是整体的能量平衡，它包括了任何含有质量的刚体的动能。由于当评价结果时我们只对变形体感兴趣，当评价能量平衡时 我们应在 Etotal 中扣除刚体的动能。 例如，如果你正在模拟一个采用滚动刚体模具的传输问题，刚体的动能可能占据模

18、型整个动能的很大部分。在这种情况下，你必须扣除与刚体运动有关的动能，然而才可 能做出与内能有意义的比较。13.5 例题： ABAQUS/Explicit 凹槽成型在这个例子中，你将应用 ABAQUS/Explicit 求解第 12 章“接触”中那个凹槽成型 的问题。然后，你将比较分别来自 ABAQUS/Standard 和 ABAQUS/Explicit 分析的结果。 你将修改由 ABAQUS/Standard 分析所创建的模型，这样才能在 ABAQUS/Explicit 中运行它。这些修改包括在材料模型中增加密度，改变单元库，并改变分析步。为了获 得正确的准静态响应， 在运行 ABAQUS/

19、Explicit 分析前， 你将应用在 ABAQUS/Standard 的频率提取过程来确定所需要的计算时间。13.5.1 前处理应用 ABAQUS/Explicit 重新运算模型对于这个模拟，应用 ABAQUS/CAE 修改模型。在本手册的在线文档第 A.12 节 Forming a channel”提供了输入文件。当通过 ABAQUS/CAE运行这个输入文件时， 将创建关于该问题的完整的分析模型。根据下面给出的指导如果你遇到困难，或者如果 你希望检查你的工作，则可以运行这个输入文件。在附录 A“ Example Files ”中，给出 了如何提取和运行输入文件的指导。在开始前，打开关于凹槽

20、成型例题的模型数据库文件，它被创建在第12.6 节“ABAQUS/Standard example ： forming a channel ”。确定一个合适的分析步时间对于一个准静态过程， 在第 13.2 节“加载速率” 中讨论了确定合适的分析步时 间的过程。如果我们知道了坯件的最低阶固有频率，即基( fundamental )频，我们就可 以确定分 析步时间的一 个大致的 下限。一种获 得这个信 息的方法是 在 ABAQUS/Standard 中运行频率分析。 在这个成型分析中， 冲压对坯件产生的变形类 似于它的最低阶模态。因此，如果你想模拟整个结构而并非局部的变形，选择第一 个成型阶段的时

21、间是大于或等于坯件最低阶模态的周期是十分重要的。运行一个固有频率提取过程:1. 将已存在的模型复制成为一个新的模型，命名为Frequency ,并对Fequency 模型进行如下全面的修改：在频率提取分析中，你将用一个单独的频率提取分 析步取代现在所有的分析步。此外，你将删除所有的刚性工具和接触相互作用；它们与确定毛坯的基频无关。2. 在Property模块中，为Steel材料模型增加一个 7800的密度。3在Assembly模块中，删除冲模、冲头和夹具部件的实体。对于频率分析并不需 要这些刚体部件。提示：你可以从工具箱中采用Delete 工具删除这些部件。4进入Step模块，用一个单独的频率

22、提取分析步替代现存的所有分析步。a. 在Step Manager （分析步管理器）中，删除分析步Remove RightCon strai nt、Holder Force 、Establish Con tact II和 MoveRunch 。b. 选择分析步 Establish Con tact I，并点击 Replace 。c. 在Repalce Step（替换分析步）对话框中，从Lin ear Parturbation过程歹U表中选择 Freque ncy ，键入分析步描述为Freque ncy modes ; 选择Lanczos特征值选项，并要求五个特征值。重新命名分析步为Extract

23、Freque ncies。d. 取消DOF Monitor（自由度监视器）选项。注意：由于频率提取分析步是一个线性扰动过程，将忽略材料的非线性性质。 在这个分析中，坯件的左端约束沿x-方向的位移和绕法线的转动；但是，没有约束沿y-方向的位移。因此，提取的第一阶模态将是刚体模态。对于在 ABAQUS/Explicit中的准静态分析，第二阶模态的频率将确定合适的时间段。5.在Interaction模块，删除所有的接触相互作用。6 .进入 Load模块，在 BC Manager （边界条件管理器）中检查在ExtractFreque ncies分析步中的边界条件。除了边界条件名称Cen terBC 以

24、外，删除所有的边界条件。将这个留下的采用了对称边界条件的毛坯约束施加到左端。7. 在创建和提交作业前，如果有必要则重新剖分网格。8. 进入Job模块，创建一个作业，命名为Forming-Frequency，采用如下的作业描述： Channel forming freque ncy an alysis。提交作业进行分析，并监控求解过程。9. 当分析完成时，进入Visualization模块，并打开由这个作业创建的输出数据库文件。从主菜单栏中，选择 Plot -Deformed Shape;或者应用在工具箱中的工具。绘制出一阶屈曲模态的模型变形形状。进一步绘出毛坯的二阶模态，将未变形 的模型形状叠

25、加在模型变形图上。频率分析表明坯件有一个140 Hz的基频，对应的周期为0.00714 s。图13-8显示了第二阶模态的位移形状。对于成型分析，我们现在知道最短的分析步时间为 0.00714 s。2图13-8 由ABAQUS/Standard频率分析的毛坯二阶模态创建ABAQUS/Explicit 成型分析成型过程的目标是采用0.03 m的冲头位移准静态地成型一个凹槽。在选择准静态分析的加载速率时，建议你在开始时用较快的加载速率，并根据需要减小加载速 率，更快地收敛到一个准静态解答。然而，如果你希望在你的第一次分析尝试中就 增加能够得到准静态结果的可能性，你应当考虑分析步时间是比相应的基频缓慢

26、10到50倍的因数。在这个分析中，对于成型分析步，你将从 0.007 s的时间开始。这 是基于在 ABAQUS/Standard中进行的频率分析， 它显示出毛坯具有 140Hz的基频， 对应于0.00714 s的时间周期。这个时间周期对应于4.3 m/s的常数冲头速度。你将仔细地检查动能和内能的结果，以检验结果中并没有包含显著的动态影响。将 Standard 模型复制成一个新模型，命名为 Explicit 。如果必要，通过 从位于工具栏下方的 Model （模型）列表中选择 Explicit 模型作为当前的模型。 使所有接下来的模型改变成为 Explicit 模型。在 ABAQUS/Stand

27、ard 分析中，在冲头和坯件之间模拟一个初始的缝隙以便于 接触计算。在 ABAQUS/Explicit 分析中则不需要采取这种预防措施。因此，在 Assembly模块中，沿U2方向平移冲头-0.001 m。在警告对话框中出现的关于相对 和绝对约束中，点击 Yes。在毛坯夹具上施加一个集中力，为了计算夹具的动态反应，必须在刚性体的参 考点上赋予一个点质量。夹具的实际质量是不重要的；而重要的是它的质量必须与 毛坯的质量（0.78 kg ）具有同一个数量级，以使在接触计算中的振荡最小化。选择 数值为 0.1 kg 的点质量。在 Property 模块中，创建一个点的截面定义，命名为 Pointmas

28、s 。在 Edit Section 对话框的 Inertial Properties域中，键入0.1 点质量的值。在参考点 RigidRefHolder 应用这个截面定义。此外，编辑 Steel 材料定义来包括 7800 kg/m 3 的质量密度。进入 Step 模块。你需要为 ABAQUS/Explicit 分析创建两个分析步。在第一个 分析步中施加夹具力；在第二个分析步中施加冲头压下力。除了命名为EstablishContact I 的分析步之外，删除所有其他的分析步，并用一个单一的显式动态分析步 替换这个分析步。键入分析步描述为 Apply holder force ，并指定 0.000

29、1 s 的分析步时间。这个时间对于施加夹具载荷是适合的，因为它是足够长以避免了动 态效果，而且又足够短以防止了对整个作业运行时间的明显冲击。将分析步重新命 名为 Holder force 。创建第二个显式动态分析步，命名为 Displace punch ， 分析步的时间为 0.007 s，键入Apply punch stroke作为分析步的描述。为了帮助确定分析是如何接近于准静态假设， 研究各种能量的历史是非常有用 的。特别有用的是比较动能和内部应变能。 能量历史默认地写入了输出数据库文件。在这个金属成型分析的第一次尝试中，对于施加的夹具力和冲头压力，你将应 用具有默认的光滑参数的表格形式的幅

30、值曲线。进入 Load 模块，为施加的夹具力 创建一个名为 Ramp1 的表格形式的幅值曲线。在表 13-1 中输入幅值数据。为冲头 压力定义第二个表格形式的幅值曲线，命名为Ramp2。在表13-2中输入幅值数据。表13-1 Rampl和Smoothl的递增幅值数据时间（s）幅值0.00.00.00011.0表13-2 Ramp2和Sooth2 的递增幅值数据时间（s）幅值0.00.00.0071.0在Load Manager （载荷管理器）中，在命名为 Holder force 的分析步中 创建一个集中力，命名为 RefHolderForce，在施加的点上指定 RefHolder 和一个沿着

31、CF2方向大小为-440000 的力。对于这个载荷，改变幅值定义为 Ramp1。在BoudaryConditionManager （边界条件 管理 器）中，删除命名为MidLeftBC 和MidRightBC 的边界条件。编辑 RefDieBC边界条件，这样在Holder force 分析步中沿着U2方向的约束为零，不改变其他方向的约束。对于 RefHolderBC边界条件，解除沿着U2方向的约束，而其他方向的约束保持不变。在Displace Punch分析步中，改变位移边界条件RefPunchBC ，使沿着 U2方向的位移为-0.03 m。对于这个边界条件，应用幅值曲线Ramp2。监视自由度

32、的值在这个模型中，你将在整个分析步中监视冲头的参考节点的竖向位移（自由度2）。在 ABAQUS/Sta ndard 成型分析中，由于已经设置了DOF Mon itor 监视RefPu nch 的竖向位移，所以你无需做出任何改变。创建网格和定义作业在网格Mesh模块中，将用于剖分坯件网格的单元族改变为Explicit，并指定 增强沙漏控制，并剖分坯件网格。因为已经将工具模拟成了解析刚性表面，因此无 需将它们剖分网格。在 Job 模块中创建一个作业，命名为 Forming-1 ，给予作业如下的描述： Channel forming - attempt 1。在运行成型分析前，你可能希望知道该分析将需

33、要多少个增量步，进而了解该分析需要多少计算机时间。你可以通过运行数据检查 （data check）分析来获得关于初始稳定时间增量的近似值， 或者你可以应用在第 13.3 节“质量放大” 中的关系式 进行估计。在这个例题中，从一个增量步到下一个增量步的稳定时间增量不会有太 大的变化，因此知道了稳定时间增量，你可以确定完成成型阶段的分析需要多少个 增量步。一旦分析开始，你就能够知道每一个增量步需要多少CPU 时间，进而知道整个分析需要多少 CPU 时间。应用在 13.3 节“质量放大” 中表述的关系式， 关于这个分析的稳定时间增量近似为1X 10-7 s。因此，对于0.007 s的分析步时间，成型

34、阶段需要大约185,000个增量步。将模型保存到模型数据库文件中，并提交作业进行分析。监视求解过程；改正 任何检测到的模拟错误，并调查任何警告信息的原因。完成整个分析可能需要运行 1 0分钟或更长的时间。一旦分析开始运行，在另一个视图窗中会显示出你选择来监视（冲头的竖向位 移）的自由度值的 X-Y 曲线图。从主菜单栏中，选择 Viewport -Job Monitor: Forming-1 ，在分析运行的整个时间中跟踪沿着2-方向冲头位移的发展进程。评价结果的策略在查看我们最关心的结果之前，诸如应力和变形形状，我们需要确定结果是否 是准静态的。一个好的方法是比较动能与内能的历史。在金属成型分析

35、中，大部分 的内能是由于塑性变形产生的。在这个模型中，坯件是动能的主要因素（忽略夹具 的运动，没有与冲头和模具相关的质量）。为了确定是否已经获得了一个可接受的 准静态解答，坯件的动能应该小于其内能的几个百分点。对于更高的精确度，特别 地是对回弹应力感兴趣时，动能应该是更低的。这个方法是非常有用的，因为它应 用于所有类型的金属成型过程，而且不需要任何直观地理解在模型中的应力；许多 成型过程可能是过于复杂，以至于不允许对结果有一个直观的判断。虽然是衡量准静态分析的良好和重要的证明， 仅凭动能与内能的比值还不足以确任解的质量。 你还必须对这两种能量进行独立地评估， 以确定它们是否是合理的。当需要准确

36、的回弹应力结果时，这一部分的评估是更增加了重要性，因为一个高度 精确的回弹应力解答是高度地依赖于准确的塑性结果。即使动能是非常小的量，如 果它包含了高度的振荡，则模型也会经历显著的塑性。一般说来，我们希望光滑加 载以产生光滑的结果；如果加载是光滑的，但是能量的结果是振荡的，则结果可能 是不合适的。由于一个能量的比值无法显示这种行为，所以你也必须研究动能本身 的历史以观察是否是光滑的还是振荡的。如果动能不能显示出准静态的行为， 在某些节点上观察速度的历史可能是有用 的，以帮助理解在各个区域中模型的行为。这种速度历史可以表明在模型的哪些区 域是振荡的，并产生大量的动能。评估结果进入 Visuliz

37、ation 模块，并打开由这个作业 ( Forming-1.odb)创建的输出数据库。绘制动能和内能。创建能量历史的曲线：1从主菜单栏中，选择 Plot - History Output 。显示出整个模型的伪应变能历史曲线。2从主菜单栏中，选择 Result -History Output 。显示出 History Output 对话框。3从变量的列表中，选择 Kinetic energy: ALLKE for Whole Model。4点击 Plot 创建一条 ALLKE 的历史曲线。显示出整个模型的动能历史曲线(见图13-9)。13-1512.502&0IA ；ooDODalx)- 0 8

38、 6 4 146山 01 0.300.002.004.00Tme6:00対“图13-9成型分析的动能历史，尝试 15.类似地，创建模型内能的历史曲线，ALLIE （见图13-10 ）。制12-00.6O.2G.AO1 1 OGJOU 山fflEse-30.00 2M 400r00 DdCTime图13-10 成型分析的内能历史，尝试113-17显示在图 13-9 中的动能历史发生显著地振荡。另外，动能的历史与坯件的成型没 有明确的关系，这表明这个分析是不适合的。在这个分析中，冲头的速度保持为常数， 而主要地依赖于坯件运动的动能却远非是恒定值。在除了开始阶段以外的整个分析步中， 比较图 13-9

39、和图 13-10 表明动能是内能的一 个很小的百分数（小于 1% ）。即使对于这种严重的加载情况，还是满足了动能必须相 对地小于内能地准则。尽管模型的动能只是内能的一个小的分数，它还是有一定的振荡。所以，我们应该 以某种方式改变模拟以获得更平滑的解答。14.5.2 成型分析尝试 2即使实际上冲头是以几乎接近于常值的速度运动， 第一次模拟尝试的结果表明理想 的方式是采用不同的幅值曲线以允许坯件更光滑地加速。 当考虑应用什么类型的加载幅 值时，记住在准静态分析的所有方面，光滑性是重要的。最偏爱的方法是尽可能光滑地 移动冲头，在理想的时间内移动理想的距离。应用一种光滑地施加的冲头力和一段光滑地施加的

40、冲头距离， 我们现在将分析成型 阶段；我们将与前面获得的结果进行比较。关于光滑步骤幅值曲线的解释，请阅读第 13.2.1 节“光滑幅值曲线”。在 Load 模块中，定义一条光滑步骤幅值曲线，命名为Smooth1 。输入在表 13-1中给出的幅值数据。创建第二条光滑步骤幅值曲线，命名为Smooth2 ，应用在表 13-2中给出的幅值数据。在 Holder force分析步中，修改 RefHolderForce 载荷，使它采用 Smooth1 的幅值。在 Displace punch 分析步中，修改位移边界条件 RefPunchBC ，使它采用 Smooth2 的幅值。 通过设置在分析步开始时的幅

41、值为 0.0 和在 分析步结束时的幅值为 1.0，ABAQUS/Explicit 创建了一个幅值定义，它的一阶和二阶 导数都是光滑的。因此，应用一条光滑步骤幅值曲线对位移进行控制，也使我们确信了 其速度和加速度是光滑的。在 Job 模块中，创建一个作业， 命名为 Forming-2 ，给予作业如下的描述： Channel forming - attempt 2。将模型保存到模型数据库文件中，并提交作业进行分析。监视求解过程；改正任何 检测到的模拟错误，并调查任何警告信息的原因。完成整个分析可能需要运行 10 分钟 或更长的时间。评估第二次尝试的结果动能的结果如图13-11所示。动能的响应是明显

42、地与坯件的成型相关：在第二 个分析步的中间阶段出现了动能的峰值，它对应于冲头速度最大的时刻。因此，动 能是适当的和合理的。图13-11 成型分析的动能历史，尝试 2CJ口COUS关于第二次尝试的内能历史如图13-12所示，显示了从零上升到最终值的光滑增长。再次看出，动能与内能的比值是相当小的，并显示出是可接受的。图13-13比较了在两次成型尝试中的内能。13-19aw2.004 006jOOTime0.00图13-12成型分析的内能历史，尝试2IE Actwinpt 1IE Attempt 2A? u 山raEec-图13-13关于成型分析的两次尝试的内能比较13-2113.5.3 两次成型尝

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

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

45、和ABAQUS/Explicit的分析，等效塑性应变的结果彼此相差在5%以内。另外，最终变形形状的比较显示出显式准静态分析的结果与ABAQUS/Standard静态分析的结果吻合得极好。-EX-KFI.? I*.I eeefl亡w V sil-ccc cri-TTC-i- flni- ret Ct r图 13-15 在 ABAQUS/Standard （左）和 ABAQUS/Explicit （右）?EES i斗叶.13-2313-#凹槽成型分析的PEEQ的等值线图13-#图13-16 在ABAQUS/Standard和ABAQUS/Explicit成型分析的凹槽最终变形形状你也应该比较由 A

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

47、00.000.00 4.00 B.OQ 12.00 16 00 20.00 24.00 28.00Punch displacement图13-17 对于ABAQUS/Standard和ABAQUS/Explicit 的稳态冲头压力比较从这些比较中，可以明显看出ABAQUS/Standard和ABAQUS/Explicit 都有能力处理诸如本例问题的困难接触分析。然而，在ABAQUS/Explicit中运行这类分析有某些优势：与ABAQUS/Standard相比，ABAQUS/Explicit 能够更容易地处理复杂 的接触条件和采用较少的分析步和边界条件进行计算。特别地，ABAQUS/Sta n

48、dard分析需要五个分析步和附加的边界条件以确保正确的边界条件和防止刚体运动。在ABAQUS/Explicit中完成同样的分析只需要两个分析步和无需附加边界条件。然而，当选择ABAQUS/Explicit进行准静态分析时，你必须明确在一个合适的加载速率下 你可能需要进行迭代。在确定加载速率时，建议你开始时采用较快的加载速率，并 根据需要减小加载速率。这可以帮助优化对分析进行求解的时间。13.5.4 加速分析的方法现在我们已经获得了一个可接受的成型分析的解答，我们可以尝试采用更短的计算机时间来获得类似的可接受的结果。因为采用显式动态标准的成型问题的实际时间是过大的，所以大部分成型分析都需要过多的

49、计算机时间以至于无法按照它们自己的物理时 间尺度进行运算；若使分析在一个可接受的计算机时间范围内运行，常常需要对分析做 出改变以减少计算机成本。有两种节省分析成本的方法：1 人为地增加冲头的速度，从而在一个更短的分析步时间内发生同样的成型过程。这种方法称为 加载速率放大(load rate scaling )。2 人为地增加单元的质量密度，从而增大稳定时间极限，允许分析采用较少的增量步。这种方法称为 质量放大(mass scaling)。这两种方法等效地做相同的事情，除非模型具有率相关材料或者阻尼。确定可接受的质量放大第13.2节“加载速率”和第13.2.3节“金属成型问题”讨论了如何确定可接

50、受 的加载速率或质量的放大因子以加速准静态分析的时间尺度。目标是在保持惯性力不显著的前提下以最短的时间模拟过程。求解的时间加快多少是有界的，而且还要 能够得到一个有意义的准静态解答。如在第13.2节“加载速率”中讨论的那样，我们可以应用同样的方法以确定一 个合适的质量放大因子，如我们已经应用以确定一个合适的加载速率放大因子的方 法。在两种方法之间的区别是加载速率放大因子f与质量放大因子f2的效果相同。最初，我们假设分析步的时间为坯件的基频周期的阶数时会产生适当的准静态结 果。通过研究模型的能量和其他的结果，我们相信这些结果是可以接受的。这项技 术产生了大约4.3 m/s的冲头速度。我们现在将接

51、受采用质量放大的求解时间，并 将结果与我们没有质量放大求解的结果进行比较，以确定由质量放大得到的结果是否可以接受。我们假设这种放大仅可能降低结果的质量，而不会使其得到改进。目 的是应用质量放大以减少计算机时间，并仍能产生可接受的结果。我们的目标是确定放大因子的值为多少时仍能产生可接受的结果，以及在哪一点上质量放大产生的结果成为不可接受的。为了观察可接受的和不可接受的放大因子的影响，在稳定时间增量尺度上，我们研究放大因子的一个范围从别的，我们选择了.5、.10和5。这些加速因子分别换算成质量放大因子为5、10和25。应用质量放大因子：1. 进入Step模块，并创建一个包含坯件的集合，命名为Biank。2. 编辑分析步 Holder force 。3. 在Edit Step（编辑分析步）对话框中，点击 Mass Scali ng（质量放大）页并选中Use scaling definitions below（使用如下放大定义）。4. 点击Create。接受半自动质量放大的默认选择。选择集合Blank作为施加的区 域，并输入一个5的值作为放大因子。在作业模块，创建一个作业，命名