3. LS-DYNA显式求解核心机制

显式动力学之所以能高效处理跌落、碰撞这类强非线性瞬态问题，根本原因在于它不组装整体刚度矩阵、不做迭代求解，而是沿着时间轴一步步“推”出结果。理解这套推进机制，是判断模型能否算得动、算得准、算得稳的前提。本节围绕中心差分法、临界时间步长、质量缩放、沙漏控制与稳定性条件展开，重点说明它们之间的因果链条。

一、中心差分法：显式推进的数学内核

LS-DYNA 的显式时间积分基于中心差分法。它把加速度用位移的二阶中心差分近似，速度用一阶中心差分近似。对节点运动方程 M·a = F_ext − F_int 而言，由于质量矩阵 M 采用集中质量（对角）形式，求逆极其廉价，于是可以直接写出节点加速度：

a_n = M⁻¹ · (F_ext,n − F_int,n)
v_(n+1/2) = v_(n−1/2) + a_n · Δt
u_(n+1)   = u_n + v_(n+1/2) · Δt

注意速度与位移在时间上错开半步，这就是“中心差分”名称的由来。整个流程没有矩阵求逆、没有平衡迭代，每个时间步的计算量近似与单元数成正比，因此可以轻松处理百万级单元模型。代价是条件稳定：时间步必须足够小，否则数值解会发散。

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

显式积分的稳定性由 Courant 条件控制。对无阻尼系统，稳定要求时间步小于临界值：

Δt ≤ Δt_crit = 2 / ω_max

其中 ω_max 是系统最高固有频率。对一维杆单元，可进一步写成 Δt_crit = L / c，L 为单元特征长度，c 为材料声速。声速由弹性模量 E 和密度 ρ 决定：c = √(E/ρ)。这意味着：单元越小、材料越硬、密度越低，临界步长越小。LS-DYNA 会在求解前扫描所有单元，取全局最小值作为实际时间步，并乘以一个安全系数（默认约 0.9）。

三、质量缩放：用“人为增重”换取时间步

当模型中存在少量极小单元时，全局时间步会被它们拖垮。质量缩放的思路是：人为增大这些单元的质量密度，使声速下降、临界步长上升。由 Δt ∝ √(ρ) 可知，若希望把时间步提高 k 倍，密度需增大 k² 倍。LS-DYNA 提供两种方式：

• 固定时间步缩放：通过 *CONTROL_TIMESTEP 设置 DT2MS 为负值，指定目标时间步，程序自动对不满足的单元加质量。
• 百分比质量增加限制：通过 DT2MS 正值或 TSSFAC 控制，避免质量增加过多导致惯性效应失真。

质量缩放不是免费的。增加的质量会改变结构惯性，若总质量增加超过 1%～2%，跌落冲击的加速度峰值和接触力可能明显偏离真实值。因此必须监控 added mass 输出，并优先通过清理网格、合并小单元来解决问题，而不是一味缩放。

四、沙漏控制：单点积分的必要代价

显式分析为节省时间，常采用单点积分（缩减积分）单元。单点积分单元只有一个积分点，无法感知某些变形模式，会产生零能变形——即沙漏模式。沙漏表现为网格出现锯齿状、棋盘状的非物理变形，能量不守恒。LS-DYNA 提供多种沙漏控制方式：

• 粘性沙漏控制：默认方式，通过 *HOURGLASS 设置系数 QH，典型值 0.05～0.10。
• 刚度沙漏控制：适用于大变形、高冲击场景，系数通常取 0.03～0.10。
• 全积分单元：彻底消除沙漏，但计算量显著增加，且可能引入剪切锁死。

判断沙漏是否受控，应查看 glstat 中的沙漏能（hourglass energy）。经验上，沙漏能应小于总内能的 5%，若超过 10% 则结果不可信，需要调整沙漏控制系数或改用全积分单元。

五、稳定性条件的整体关系

上述机制并非孤立，而是形成一条完整的因果链：中心差分法要求条件稳定 → 稳定性由临界时间步长决定 → 临界步长受最小单元和材料声速限制 → 质量缩放可提升步长但引入惯性误差 → 单点积分提升效率但带来沙漏风险 → 沙漏控制消耗额外计算并影响能量。理解这条链条，才能在建模时做出合理取舍。

六、实操中的参数设置与检查

在 LS-DYNA 关键字文件中，与本节直接相关的控制卡片主要有三张。下面给出一个典型跌落分析的示例配置（参数为教学示例，需根据实际模型调整）：

*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
       0.0       0.9         0       0.0    -1.0e-6         0         0         0

*CONTROL_HOURGLASS
$#    ihq       qh
         1       0.10

*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen
         2         2         2         2

其中 DT2MS = -1.0e-6 表示目标时间步为 1 μs，程序会对不满足该步长的单元自动加质量；IHQ = 1 表示采用粘性沙漏控制，QH = 0.10 为控制系数；HGEN = 2 要求输出沙漏能，便于后处理检查。

求解完成后，打开 glstat 文件，重点查看三列数据：总能量（total energy）是否守恒、沙漏能（hourglass energy）占比、附加质量（added mass）百分比。若沙漏能持续增长超过内能的 5%，或附加质量超过 2%，应回到模型检查最小单元尺寸和网格质量，而不是简单调大沙漏系数掩盖问题。

总结来看，显式分析的稳定性是一条环环相扣的链条。中心差分法决定了条件稳定的本质，临界时间步长给出了硬性上限，质量缩放提供了突破上限的手段但附带精度代价，沙漏控制则保障了单点积分单元的可用性。只有同时监控时间步、附加质量与沙漏能三个指标，才能判断一次跌落仿真是否既算得动、又算得准。