2. 显式动力学理论基础

显式动力学求解的核心，是把结构动力学问题写成随时间推进的离散运动方程，再用中心差分法逐步“推”出每个时刻的位移、速度和加速度。理解这条主线，才能明白为什么显式分析适合冲击、碰撞这类短时强非线性问题，也才能判断计算是否稳定、结果是否可信。

一、运动方程：从连续体到离散节点

在有限元离散后，结构动力学方程可写为：

M·a(t) + C·v(t) + F_int(u,t) = F_ext(t)

其中 M 为质量矩阵，C 为阻尼矩阵，F_int 为内力（由单元应力积分得到），F_ext 为外力，a、v、u 分别为节点加速度、速度、位移。显式分析通常采用集中质量矩阵（lumped mass），使 M 为对角阵，从而无需组装和求解整体刚度矩阵，也不需要进行矩阵求逆。这正是显式方法每步计算量小、适合大规模并行计算的根本原因。

在冲击问题中，惯性项和接触力往往占主导，阻尼影响相对次要。因此实际建模时，材料本构、接触定义和网格质量比阻尼设置更值得优先关注。

二、中心差分法：时间推进的“显式”含义

中心差分法用相邻时刻的位移差近似速度和加速度：

v(t) ≈ [u(t+Δt) - u(t-Δt)] / (2Δt)
a(t) ≈ [u(t+Δt) - 2u(t) + u(t-Δt)] / Δt²

代入运动方程并整理，可得位移递推公式：

u(t+Δt) = f(u(t), u(t-Δt), M, F_ext, F_int)

也就是说，下一时刻的位移只依赖当前和上一时刻的已知量，不需要迭代求解方程组。这种“显式”推进方式每步计算量小，但代价是时间增量必须足够小，否则数值解会发散。

三、稳定时间增量与 Courant 条件

显式分析的稳定性由稳定时间增量 Δt_stable 控制。对于无阻尼系统，常用估计为：

Δt_stable = 2 / ω_max

其中 ω_max 是系统最高固有频率。对一维杆单元，可进一步写成 Courant 条件形式：

Δt ≤ L_e / c

这里 L_e 是单元特征长度，c 是材料波速。对线弹性材料，c = √(E/ρ)。这意味着：网格越小、材料刚度越大、密度越小，稳定时间增量就越小，计算步数越多。冲击仿真中，为了捕捉应力波传播，网格往往较密，因此显式分析的时间步长常在 1e-7～1e-9 s 量级。

四、沙漏控制：为什么需要抑制零能模式

显式分析常用减缩积分单元以提高效率，但减缩积分会引入沙漏模式——单元在无应变能增加的情况下发生非物理变形。若不控制，沙漏会迅速累积并污染结果。常用控制方式包括：

• 沙漏刚度：在单元中引入人工刚度，抑制零能变形；
• 沙漏阻尼：通过人工阻尼耗散沙漏能量；
• 增强沙漏控制：适用于大变形和冲击问题，精度更高但计算量略增。

判断沙漏是否受控，不能只看总能量，还要看沙漏能（ALLAE）与内能（ALLIE）的比值。工程上通常要求沙漏能不超过内能的 5%～10%。若超标，应优先细化网格、改用增强沙漏控制或调整单元类型。

五、能量平衡：结果可信度的“体检指标”

显式分析中，能量平衡是判断结果是否可信的重要依据。Abaqus/Explicit 输出的能量项包括：

总能量应近似守恒：外力功 ≈ 内能 + 动能 + 耗散能。若总能量出现明显漂移，或沙漏能、粘性耗散能异常增大，说明模型可能存在网格畸变、接触穿透或沙漏失控等问题。

实际建模时，建议先做一次短时“试算”：只计算几百步，观察稳定时间增量、沙漏能和总能量曲线。若稳定增量过小，优先检查最小单元尺寸；若沙漏能偏高，优先调整单元类型或沙漏控制；若总能量漂移明显，再排查接触、边界条件和材料参数。把这三项作为冲击仿真的常规体检项，能显著降低“算得出来但结果不可信”的风险。