2. Abaqus/Explicit求解器基础

Abaqus/Explicit 是汽车碰撞仿真中最常用的求解器，它采用显式时间积分方法推进动力学方程。与 Abaqus/Standard 的隐式求解不同，显式方法不需要在每个增量步组装和求解整体刚度矩阵，而是通过已知的当前状态直接推算下一时刻状态。这一特性使它在处理接触、大变形、材料失效和高度非线性问题时具有明显优势，但也带来了稳定时间增量、质量缩放和沙漏控制等一系列必须理解并正确设置的工程问题。

显式时间积分的基本原理

显式方法求解的运动方程可写为：

M·a = F_ext − F_int
a = M⁻¹·(F_ext − F_int)

其中 M 为质量矩阵，a 为节点加速度，F_ext 为外力，F_int 为内力。由于显式方法通常采用集中质量矩阵，M 为对角矩阵，求逆只需对每个自由度做除法，计算量极小。得到加速度后，通过中心差分法更新速度和位移：

v_(n+1/2) = v_(n−1/2) + a_n · Δt
u_(n+1)   = u_n + v_(n+1/2) · Δt

整个过程不需要迭代，也不存在收敛判断，每个增量步的计算成本基本固定。这正是显式方法能够处理百万级单元整车碰撞模型的原因。但代价是时间增量必须足够小，否则数值解会迅速发散。

稳定时间增量

显式积分的稳定性由 Courant-Friedrichs-Lewy 条件决定。对于无阻尼系统，稳定时间增量近似为：

Δt_stable = 2 / ω_max

其中 ω_max 是系统最高固有频率。在有限元网格中，最高频率通常由最小单元控制，因此工程上常用下式估算：

Δt ≈ L_min / c
c = sqrt(E / ρ)

L_min 为最小单元特征长度，c 为材料波速。以钢为例，E≈210 GPa，ρ≈7.8×10⁻⁹ t/mm³，波速约为 5190 mm/ms。若最小单元尺寸为 5 mm，则稳定时间增量约为 0.96 μs。整车碰撞仿真通常持续 100～150 ms，若全部按此增量计算，需要上亿个增量步，计算量巨大。因此实际工程中必须通过质量缩放或网格控制来调整时间增量。

质量缩放

质量缩放的核心思想是人为增大密度，从而降低波速、增大稳定时间增量。由 Δt ∝ sqrt(ρ/E) 可知，密度增大 n 倍，时间增量增大 sqrt(n) 倍。Abaqus/Explicit 提供固定质量缩放和变质量缩放两种方式。固定质量缩放对整个模型或指定单元集统一放大密度；变质量缩放则根据目标时间增量自动调整每个单元的质量，使所有单元的时间增量接近设定值。

质量缩放必须谨慎使用。一般建议质量增加控制在模型总质量的 1%～5% 以内，且不能改变关键部件的惯性特性。对于汽车碰撞中的薄壁吸能件、B 柱、门槛梁等，过度质量缩放会显著改变变形模式和吸能量，导致结果失真。实际操作中可在 History Output 中监控 DMASS 和 TMASS，若质量增加比例超过 5%，应优先检查网格质量而非继续加大缩放。

沙漏控制

显式方法常采用减缩积分单元以提高效率，但减缩积分单元在弯曲变形下容易出现零能模式，即沙漏。沙漏模式表现为网格呈现锯齿状、无物理意义的变形，且不消耗能量。Abaqus/Explicit 默认对减缩积分单元启用沙漏控制，常用方法包括刚度沙漏控制和粘性沙漏控制。刚度沙漏控制适用于以强度为主的问题，粘性沙漏控制适用于以冲击和波动为主的问题。

判断沙漏是否严重，可查看能量输出中的 ALLAE（人工应变能）与 ALLIE（内能）之比。一般要求 ALLAE 不超过 ALLIE 的 5%，若超过 10% 则结果可信度明显下降。此时应加密网格、改用全积分单元或调整沙漏控制参数。

能量平衡检查

能量平衡是判断显式仿真是否可信的核心依据。Abaqus/Explicit 输出的能量项包括：内能 ALLIE、动能 ALLKE、伪应变能 ALLAE、外部功 ALLWK、阻尼耗散 ALLVD 等。总能量应满足：

ALLKE + ALLIE + ALLVD + ALLAE ≈ ALLWK + ALLKE_0 + ALLIE_0

工程判断中重点关注三点：第一，伪应变能 ALLAE 与内能 ALLIE 的比值应小于 5%；第二，动能 ALLKE 在碰撞过程中应远小于内能 ALLIE，通常峰值不超过内能的 10%～20%；第三，总能量曲线应保持平稳，不应出现突变或持续漂移。若能量不守恒，往往意味着接触穿透、沙漏严重或时间增量设置不当。

碰撞仿真适用性判断

显式方法并非适用于所有工况。对于准静态压溃、低速蠕变或长时间恒载问题，显式方法需要极长的计算时间，此时隐式方法更合适。对于高速碰撞、爆炸冲击、跌落和瞬态大变形问题，显式方法具有明显优势。判断依据可归纳为：载荷持续时间与结构自振周期之比。若载荷持续时间远大于结构自振周期，属于准静态问题，优先选隐式；若载荷持续时间与自振周期相当或更短，属于动态问题，优先选显式。

在实际项目中，建议先完成网格收敛性检查，再确定质量缩放策略，最后通过能量曲线验证结果可信度。若 ALLAE 偏高，优先加密网格或改用全积分单元；若 ALLKE 异常偏大，检查是否存在未约束的刚体运动或接触穿透；若总能量不守恒，检查时间增量是否超过稳定极限。只有能量平衡通过，碰撞仿真的加速度、侵入量和吸能量结果才具备工程参考价值。