2. 显式动力学理论基础
显式动力学之所以能高效处理金属成形这类强非线性、大变形、接触状态频繁切换的问题，核心在于它不组装也不求解整体刚度矩阵，而是用中心差分法在时间域上逐步推进。理解这条时间推进链上的四个关键概念——中心差分、条件稳定、临界时间步长、质量缩放与沙漏控制——是后续所有建模决策（网格尺寸、材料密度、时间步设置、单元类型选择）的理论依据。
一、中心差分法：显式推进的数学骨架
显式方法把运动方程写成牛顿第二定律的离散形式。在时刻 tn，节点加速度由内力与外力的合力除以节点质量得到：
a_n = M⁻¹ · (F_ext − F_int)
v_{n+1/2} = v_{n−1/2} + a_n · Δt
u_{n+1} = u_n + v_{n+1/2} · Δt这里 M 采用集中质量矩阵（对角阵），求逆只是逐节点做除法，无需矩阵分解。速度定义在半步 n±1/2，位移定义在整步 n，这种“交错”格式就是中心差分的典型特征。整个推进过程只依赖当前和上一步的状态，不需要迭代收敛判断，因此每个时间步的计算量近似线性于单元数，且天然适合并行。
需要注意：显式方法不检查静力平衡，它求解的是瞬态动力响应。金属成形中惯性效应、应力波传播本身就是物理过程的一部分，这正是显式方法适用的原因。
二、条件稳定与临界时间步长
中心差分法是条件稳定的：只有当时间步长 Δt 小于某个临界值 Δtcrit 时，数值解才不会发散。其物理含义是——一个时间步内，信息（应力波）传播的距离不能超过最小单元的特征尺寸。对一维杆单元，稳定条件为：
Δt ≤ Δt_crit = L_e / c
c = √(E / ρ) （弹性波速）
考虑阻尼与非线性后常取安全系数：
Δt = 0.9 · L_e / c其中 Le 是单元特征长度，c 是材料声速。对三维实体单元，LS-DYNA 内部按单元体积与最大面面积估算特征长度，并取所有单元中最小的 Δtcrit 作为全局时间步。这意味着整个模型的时间步由最小的那个单元决定——一个局部细化的小单元会把全局步长拖到很小，计算时间成倍增加。
三、质量缩放：用密度换时间步
由 Δtcrit = Le/√(E/ρ) 可知，增大密度 ρ 可以增大临界时间步。质量缩放就是人为提高小单元或全部单元的密度，把过小的 Δt 抬升到可接受水平。LS-DYNA 中常用两种方式：
- 固定时间步缩放：指定目标 Δt，程序自动为需要缩放的单元增加密度，使 Δtcrit ≥ 目标值。
- 百分比质量增加：限制总质量增加比例（如 ≤1%～5%），避免惯性效应被过度放大。
质量缩放不是免费的：增加质量会改变惯性力和波速，可能影响薄板回弹、起皱等对惯性敏感的结果。工程上通常只对少数畸变严重或尺寸过小的单元做缩放，并监控“added mass”日志，把总质量增加控制在 1% 以内。若必须大幅缩放，应优先检查网格质量，而不是依赖质量缩放掩盖问题。
四、沙漏控制：单点积分的代价与对策
显式分析为节省时间，常采用单点积分实体单元（如 ELFORM=1）。单点积分只在一个积分点计算应力，计算极快，但会产生一种非物理的零能变形模式——沙漏模式：单元在无应变能增加的情况下发生交替的锯齿状变形，表现为网格出现“棋盘格”式振荡，能量从物理变形中泄漏。
沙漏控制的本质是给这种零能模式施加一个小的抵抗刚度或阻尼。常用方法：
- 粘性沙漏控制（IHQ=1～3）：施加与变形速率相关的阻尼力，适合高速冲击。
- 刚度沙漏控制（IHQ=4～5）：施加与变形量相关的弹性恢复力，适合准静态或低速成形。
- 系数 QH：默认 0.05～0.10，过大将人为刚化模型，过小则无法抑制沙漏。
判断沙漏是否受控，应查看 glstat 中的 hourglass energy：一般要求其小于内能的 5%，若持续增长并超过 10%，说明沙漏严重，需要调整控制方式、细化网格或改用全积分/选择性缩减积分单元。
五、四个概念的联动关系
这四个概念不是孤立的：中心差分法给出推进格式，条件稳定给出 Δt 的上限，临界时间步长把上限与网格尺寸、材料属性绑定，质量缩放提供在网格无法再优化时“抬高上限”的手段，而沙漏控制则保证单点积分单元在稳定推进中不产生虚假能量。实际建模时，建议按以下顺序检查：
- 先看最小单元尺寸和材料声速，估算 Δtcrit，判断是否需要质量缩放。
- 若启用质量缩放，记录 added mass 比例，确保不超过 1%～5%。
- 运行后查看 glstat 的 time step 是否稳定，hourglass energy 是否受控。
- 若沙漏能偏高，优先调整 IHQ 类型与 QH 系数，其次考虑局部网格细化或更换单元公式。
掌握这条“差分—稳定—步长—缩放—沙漏”的因果链，就能在后续建模中把时间步、网格与单元类型的选择从“凭经验试”变成“有依据算”。