Abaqus-原来显式计算也可以这么快
本文详细介绍了 Abaqus/explicit 中的准静态分析及其提高求解速度的方法。
我如何知道我的模拟是否是准静态的？
准静态问题是通常可以使用 Abaqus/Standard 解决的问题之一，但由于触或材料的复杂性，可能难以收敛，从而导致大量迭代。具有挑战性的非线性准静态问题通常涉及：
接触条件非常复杂，Abaqus/Standard 可能会因接触问题而无法收敛；
非常大的变形可能导致严重的网格变形；
例如，通常在金属成形分析中，我们面临这样的困难：用 Abaqus/Standard 来模拟这样的问题确实很困难。
示例：模拟拉深过程中的撕裂
Abaqus/Explicit显式准静态分析问题
Abaqus/Explicit 对于高度非线性静态（准静态分析）问题的建模更加有效。对于涉及接触和金属成型等非常大变形的三维问题尤其如此。
应用 Abaqus/Explicit 模拟准静态事件需要特别考虑。在自然时间段内对过程进行建模在计算上是不切实际的。从字面上看，需要数百万次的时间增量。因此，我们在模拟中人为地提高加载过程的速度以获得经济的解决方案。
使用 Abaqus/Explicit 获得经济的准静态分析解决方案的两种方法是：
1、人为提高加载速率
我们可以通过提高加载率来人为地减少该过程的时间尺度。增加的加载速率会减少模拟的时间。将加载速率提高 f 倍，分析速度提高 f 倍。
2、采用质量缩放
它增加了稳定时间增量的大小，因此完成作业所需的增量更少。人为地将材料密度（质量缩放）增加 f*f 倍，可以将分析速度提高 f 倍。
在本文中，我们的重点是提高加载速率的方法。
为了减少 Abaqus/Explicit 分析中所需的增量数量，我们可以比实际过程的时间加快模拟速度，也就是说，我们可以人为地缩短事件的时间周期，或者同等地提高事件发生的速率。正在加载。这会引入可能的错误。如果加载速率增加太多，增加的惯性力将改变预测的响应。在极端情况下，问题将表现出波传播响应。避免此错误的唯一方法是选择不太大的加载速率。
如何判断加载速率是否合适？
以不同的加载速率运行多个模拟：
按照从最快加载速率到最慢加载速率的顺序运行一系列模拟。如您所知，加载速率越慢，分析时间就越长。检查结果（变形的形状、应力、应变和能量），以了解在更改 Abaqus 载荷率时改变模型的影响。金属板材成型过程中刀具速度过高往往会导致不切实际的局部拉伸；成型模拟中工具速度过高会导致喷射（流体动力型响应）；过高的加载速率会导致施加负载附近的高度局部变形；由于初始变形的（非结构性）阻力增加，准静态倒塌分析中的过高加载速率可能会导致载荷。
使用固有频率检查加载速率：准静态分析的主要响应将是第一结构模式。因此，我们使用该模式的频率来估计合适的 Abaqus 加载速率。
1、估计模型的第一固有频率 (f)。在简单模型中，我们可以通过可用的分析关系找到该频率。对于更复杂的模型，首先在 Abaqus 中运行频率分析。
2、使用模型的第一固有频率计算相应的时间段 (T)：T=1/f
3、运行显式分析（步长时间 = T）并估计在此时间 (T) 期间模型冲击方向的全局偏转 (D)。
4、计算冲击速度（V）：V=D/T
5、一般建议是将冲击速度限制在材料波速的 1% 以下。金属中的典型波速为 5000 。
示例（门梁入侵测试）
为了说明确定适当加载速率的问题，请考虑车门中侧防盗梁的变形。实际测试是准静态的。
我们将测试建模为圆形梁（长度为 l、直径为 d、厚度为 t）固定在两端，刚性圆柱体（直径为 D）使梁变形。
在这里，我们检查圆柱体的 20 m/s 和 400 m/s 速度，看看哪一个可以适用于我们的问题。
1、第一模态的频率约为250 Hz：f=250
2、该速率对应于 4 毫秒的周期：T=1/250=0.004 s
3、使用 20 m/sec 的速度，分析显示圆柱体将在 4 ms内被推入梁 0.1
4、冲击速度为：V=D/T= 0.08/0.004= 20 m/s
5、回想一下，金属中的波速约为 5000 m/sec，因此 25 m/sec 的冲击速度约为波速的 0.5%（小于 1%）。
如果我们检查 400 m/s 的速度，将导致大约 4% 的波速（不可接受）。
小结
随着过程速度的增加，静态平衡状态演变成动态平衡状态，惯性力变得更加占主导地位。我们应该尝试在惯性力仍然微不足道的最短时间段（Abaqus 最大加载速率）内对过程进行建模；
除了惯性力之外，问题的某些方面（例如材料行为）也可能与速率相关。在这种情况下，无法更改正在建模的事件的实际时间段。质量缩放方法在此类问题中变得有吸引力。
