3. 显式时间积分与稳定性控制

显式动力学之所以能高效处理碰撞、冲击、爆炸等强非线性问题，核心在于它采用中心差分法对时间进行离散，避开了隐式方法中每步求解大型线性方程组的开销。但代价是：显式积分是条件稳定的，时间步长必须小于一个临界值，否则数值解会迅速发散。本节围绕中心差分原理、临界时间步长、质量缩放、沙漏控制与能量平衡检查五个相互关联的问题展开。

一、中心差分法：显式积分的核心机制

在显式动力学中，位移、速度、加速度按交错的时间点存储，速度定义在半步点，位移和加速度定义在整步点。给定 t_n 时刻的位移 u_n，加速度由运动方程直接得到：

a_n = M⁻¹ · (F_ext − F_int)
其中 M 为集中质量矩阵（对角阵），F_int 为单元内力

由于 M 是对角阵，求逆只需逐分量相除，无需矩阵分解，这是显式方法每步计算量极小的根本原因。随后用中心差分更新速度与位移：

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

整个流程无需迭代、无需组装整体刚度矩阵，每步只做向量运算，因此单步成本与单元数近似线性。但这也意味着：一旦 Δt 超过稳定极限，误差会在每一步被放大，最终导致位移、能量爆炸式增长。

二、临界时间步长：稳定性的硬约束

显式积分的稳定条件由 Courant 条件给出。对线弹性材料，临界步长近似为：

Δt_cr = 2 / ω_max = L_e / c
c = √(E / ρ)   （声速，E 为弹性模量，ρ 为密度）
L_e 为单元特征长度

实际计算中，LS-DYNA 会扫描所有单元，取最小特征长度与最大声速的组合，再乘以一个安全系数（默认 0.9 左右）：

Δt = 0.9 · min( L_e / c )

这意味着最小、最硬的单元决定全局时间步长。例如钢材 E≈210 GPa、ρ≈7850 kg/m³，声速约 5170 m/s；若某个单元特征长度仅 1 mm，则 Δt_cr≈0.19 μs。若模型中有 10 万个这样的单元，模拟 10 ms 就需要约 5 万步，计算量巨大。因此，控制最小单元尺寸、避免过度细化，是显式建模的第一原则。

三、质量缩放：用精度换效率

当模型中存在少量极小单元时，全局步长会被严重拖累。质量缩放的基本思路是：人为增大这些单元（或整体）的密度，使声速下降、临界步长上升。LS-DYNA 提供两种方式：

• 全局质量缩放：通过 *CONTROL_TIMESTEP 的 DT2MS 参数指定目标步长，程序自动对不满足步长要求的单元增加质量。若 DT2MS 为正值，质量只加到小单元上；为负值则所有单元都加质量。
• 局部质量缩放：用 *PART 或 *ELEMENT_MASS 对特定部件单独调整，适合已知某区域网格过密的场景。

质量缩放的代价是引入附加惯性，可能改变结构的动力学响应。工程上通常控制总质量增加不超过 1%～5%，并在 d3hsp 文件中检查实际添加的质量比例。若某次缩放导致总质量增加 20%，即使步长变大，结果也可能已不可信。

四、沙漏控制：单点积分的必要补充

显式分析普遍采用单点积分单元以降低计算量，但单点积分会产生沙漏模式——单元在无外力做功的情况下发生交替变形，能量不守恒。LS-DYNA 通过 *HOURGLASS 关键字施加人工刚度或黏性阻尼来抑制。常用类型包括：

• 类型 1（标准刚度）：默认选项，适用于大多数实体单元。
• 类型 4（刚度系数）：对刚度较大的问题更稳定。
• 类型 5/6（黏性）：适合高速冲击，但可能引入额外阻尼。

沙漏控制系数 QM 默认约 0.05～0.1。系数过小无法抑制沙漏，过大则使单元过刚、结果偏硬。判断标准是：沙漏能应小于内能的 5%，超过 10% 通常说明网格或控制参数需要调整。

五、能量平衡检查：结果可信度的最终判据

显式分析结束后，必须检查能量历史。LS-DYNA 通过 glstat（全局）和 matsum（材料）文件输出能量。核心关系为：

总能量 = 动能 + 内能 + 沙漏能 + 滑移界面能 + 阻尼能 + 外部功
理想情况：总能量在无外力时保持常数

检查要点：

• 沙漏能 / 内能 < 5%：否则沙漏控制不足或网格过粗。
• 动能 / 内能：冲击问题中动能应逐渐转化为内能，若动能异常增大，可能发生数值不稳定。
• 总能量曲线：应平滑，若出现突跳或持续上升，说明时间步长或接触设置有问题。
• 质量增加比例：在 d3hsp 中查看，超过 5% 需重新评估质量缩放策略。

中心差分法、临界步长、质量缩放、沙漏控制与能量检查构成一个闭环：步长决定效率，质量缩放调节步长，沙漏控制保证单点积分可用，能量检查则验证整个链条是否可信。任何一环失控，都会让仿真结果失去工程意义。