2. 显式动力学理论基础

显式动力学求解冲击问题的核心，是把时间推进交给一个显式的中心差分格式，而不是像隐式那样在每个增量步里反复迭代求解整体刚度方程。理解这一点，是后续设置稳定时间增量、判断质量缩放是否合理、以及用能量曲线判断结果可信度的前提。本节围绕四个相互咬合的概念展开：中心差分法如何推进、稳定时间增量由什么决定、条件稳定性意味着什么、质量缩放与能量平衡判据如何配合使用。

中心差分法：用当前状态直接算出下一步

在显式动力学中，位移、速度、加速度按时间离散。设第 n 步的位移为 u_n，中心差分法用相邻两步位移之差近似速度，用相邻两步速度之差近似加速度：

u̇_(n+1/2) = (u_(n+1) − u_n) / Δt
ü_n       = (u̇_(n+1/2) − u̇_(n−1/2)) / Δt

把运动方程 M·ü + C·u̇ + F_int = F_ext 代入，并利用质量矩阵 M 通常为对角（集中质量）矩阵这一特点，可以不解线性方程组，直接得到下一步位移：

u_(n+1) = f(M, C, Δt, u_n, u_(n−1), F_int, F_ext)

这就是“显式”的含义：每一步的计算量主要花在单元内力 F_int 的组装上，而不需要矩阵求逆或迭代。代价是时间增量 Δt 必须足够小，否则数值解会发散——这就引出了稳定时间增量。

稳定时间增量：由最小单元与材料波速决定

显式求解的稳定时间增量由系统的最高固有频率决定。工程上常用的是基于单元尺度的估计式：

Δt_stable = L_e / c_d
c_d = sqrt(E / ρ)          （线弹性、一维杆波速）
c_d = sqrt((λ + 2μ) / ρ)   （三维实体膨胀波速）

其中 L_e 是单元特征长度，c_d 是材料中的膨胀波速，E 为弹性模量，ρ 为密度，λ、μ 为拉梅常数。Abaqus/Explicit 会在每个增量步自动估算全局稳定时间增量，并乘以一个安全系数（默认 0.9 左右）作为实际使用的 Δt。也就是说，真正决定计算步数的是模型里“最小、最硬”的那个单元：单元越小、材料波速越高，Δt 越小，步数越多。

举一个教学示例：一块钢材，E ≈ 210 GPa，ρ ≈ 7850 kg/m³，则 c_d ≈ 5170 m/s。若最小单元特征长度 L_e = 1 mm，则 Δt_stable ≈ 1×10⁻³ / 5170 ≈ 1.9×10⁻⁷ s，即约 0.19 μs。若冲击过程持续 2 ms，则至少需要约 1 万步以上。这个数量级关系，是判断模型规模是否可行的第一把尺子。

条件稳定性：为什么显式必须“小步快跑”

中心差分法是条件稳定的：只有当 Δt 小于某个临界值时，数值解才不会随步数指数放大。这个临界值就是上面估算的 Δt_stable。一旦超过，位移、速度会迅速振荡并溢出，Abaqus/Explicit 会报出 “time increment is too large” 或直接终止。因此显式分析不是“想走多大步就走多大步”，而是被最小单元和材料波速牢牢约束。

这也解释了显式与隐式的分工：冲击、碰撞、爆炸这类过程本身持续时间短（毫秒到微秒级），显式虽然步数多，但每步极快，总成本可接受；而准静态问题若用显式，往往需要人为加速加载或质量缩放，否则步数会大到不现实。

质量缩放：用“人为增重”换取更大时间增量

由 Δt ≈ L_e/√(E/ρ) 可知，增大密度 ρ 可以降低波速、增大 Δt。质量缩放（mass scaling）就是有控制地放大局部或全局质量，把过小的 Δt 抬到可接受水平。常见做法有两种：

• 固定时间增量缩放：设定目标 Δt，让 Abaqus 自动对需要缩放的单元加质量，使全局 Δt 不低于该值。
• 直接指定质量缩放因子：对特定单元集或整体乘以一个系数，适合已知问题出在少数小单元上的情形。

质量缩放不是免费的。它改变了系统的惯性，会人为降低波速、改变冲击响应的高频成分。判断是否可接受，不能只看“算得快不快”，而要看能量平衡。

能量平衡判据：用能量曲线判断结果是否可信

Abaqus/Explicit 在历史输出中提供一组能量：内能 ALLIE、动能 ALLKE、伪应变能 ALLAE、外力功 ALLWK、总能量 ETOTAL 等。判断质量缩放是否过度的核心指标是动能与内能的比值，以及伪应变能占比。

实操中，建议先做一次无质量缩放的试算，记录最小稳定增量与能量曲线；再逐步加入质量缩放，观察 ALLKE/ALLIE 是否仍处于可接受区间。若缩放后动能占比明显抬升，说明缩放量已经影响到物理响应，应改为局部缩放或细化网格，而不是继续加大缩放因子。