2. LS-DYNA显式动力学理论基础

金属成形仿真中，冲压、锻造、旋压等过程往往伴随大变形、接触碰撞和高度非线性，隐式方法在求解这类问题时容易遇到收敛困难。LS-DYNA 采用显式动力学方法，通过中心差分时间积分在时间域上逐步推进，不需要在每个时间步组装和求解大型刚度矩阵，也不进行迭代收敛判断，因此特别适合处理短时、强非线性、以惯性效应为主导的成形问题。理解运动方程、中心差分、条件稳定性、临界时间步长、沙漏控制与质量缩放这六个概念之间的逻辑关系，是正确设置求解参数、判断结果可信度的前提。

运动方程与显式求解的基本思路

显式动力学求解的出发点是离散后的运动方程。对每个节点写出牛顿第二定律，并考虑阻尼与内力，可写成矩阵形式：

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

其中 M 为质量矩阵，C 为阻尼矩阵，a、v 分别为节点加速度与速度向量，F_int 为内力（由单元应力经应变-位移关系集成得到），F_ext 为外力（接触力、压力、重力等）。LS-DYNA 通常采用集中质量矩阵（lumped mass），即 M 为对角阵，这使得加速度可以直接由节点力除以节点质量得到，无需矩阵求逆：

a(t) = M⁻¹ · [F_ext(t) − F_int(t) − C·v(t)]

这一步是显式方法高效的核心：每个时间步的计算量主要花在单元内力计算和接触搜索上，而不是线性方程组求解。内力通过单元在当前构型下的应变计算应力，再积分得到节点力，因此几何非线性、材料非线性都被自然包含在逐步更新中。

中心差分时间积分

LS-DYNA 采用中心差分法对时间进行离散。已知第 n 步的位移 dⁿ、速度 vⁿ⁻¹ᐟ²、加速度 aⁿ，则：

v^(n+1/2) = v^(n−1/2) + a^n · Δt
d^(n+1)   = d^n + v^(n+1/2) · Δt

速度定义在半步时刻，位移定义在整步时刻，这种交错格式使积分具有二阶精度。整个推进过程是显式的：当前步的加速度只依赖当前步的力和已知速度，下一步位移可直接算出，不需要求解未来时刻的耦合方程。对金属成形而言，这意味着每一增量步都能明确给出节点位置、单元应变和应力状态，便于观察板料减薄、回弹趋势和接触压力分布。

条件稳定性与临界时间步长

显式中心差分是有条件稳定的：时间步长必须小于临界值，否则数值解会迅速发散。临界时间步长由系统中最小单元的特征尺寸和材料波速决定，常用估算式为：

Δt_crit = L_e / c
c = sqrt(E / ρ)

其中 L_e 为单元特征长度，c 为材料中的声速（弹性波速），E 为弹性模量，ρ 为密度。对壳单元，LS-DYNA 还会考虑弯曲刚度和剪切修正，实际使用的临界步长通常由程序根据单元几何与材料自动计算，并乘以一个安全系数（默认约 0.9）。如果模型中存在极小单元、严重畸变单元或高刚度材料，Δt_crit 会显著减小，导致计算时间急剧增加。因此，网格质量直接决定求解效率，成形仿真中应尽量避免长宽比过大或过度细化的局部网格。

沙漏控制

显式求解中，壳单元和实体单元常采用单点积分以降低计算量。单点积分只有一个积分点，单元在受力时可能出现一种非物理的、零能量的变形模式，称为沙漏模式。它表现为网格呈现锯齿状或交替扭曲，而单元本身不产生应变能，因此不会被内力约束。如果不加控制，沙漏能会不断累积，污染位移和应力结果。

LS-DYNA 提供两类沙漏控制：粘性沙漏控制和刚度沙漏控制。粘性控制通过引入与变形速率相关的阻尼力抑制沙漏，适用于大多数冲击和成形问题；刚度控制则增加与沙漏变形相关的刚度，适用于需要更严格抑制沙漏的情况。工程上通常关注沙漏能与内能的比值，一般要求控制在 5% 以内。若沙漏能偏高，可通过细化网格、改用全积分单元、调整沙漏系数或优化加载方式改善。在金属成形中，板料弯曲区域和接触区域容易出现沙漏，应重点检查这些部位的沙漏能曲线。

质量缩放

质量缩放是提高显式求解效率的常用手段。由 Δt_crit = L_e / √(E/ρ) 可知，增大密度 ρ 可以增大临界时间步长，从而减少总时间步数。LS-DYNA 允许通过质量缩放将小于目标步长的单元质量人为放大，使全局时间步长提高到设定值。例如，某板料单元特征长度 0.5 mm，材料为钢，E = 210 GPa，ρ = 7.85×10⁻⁹ t/mm³，则声速 c ≈ 5170 mm/s，临界步长约 9.7×10⁻⁵ s。若模型中存在 0.1 mm 的小单元，临界步长降至约 1.9×10⁻⁵ s，计算步数增加约 5 倍。此时可设置目标步长 5×10⁻⁶ s 或对最小单元进行质量缩放，但必须监控质量增加比例。

质量缩放并非没有代价。人为增加质量会改变惯性效应，可能影响高速成形、回弹和薄板起皱的预测精度。一般建议质量增加控制在 1% 以内，最多不超过 5%，并且只对局部小单元进行缩放，避免全局密度改变。对于准静态成形过程，如果加载时间被压缩，质量缩放的影响相对较小；但对于高速冲击或爆炸成形，必须谨慎使用，必要时通过能量比和动量守恒进行校核。

参数设置与结果判断的实践要点

在实际建模中，这六个概念是相互关联的。网格尺寸决定临界步长，临界步长决定计算成本，质量缩放和沙漏控制则是在成本与精度之间做权衡。一个可参考的设置流程是：先检查最小单元尺寸和材料波速，估算自然临界步长；再根据可接受的计算时间设定目标步长，确定是否需要质量缩放；然后选择沙漏控制类型和系数，通常粘性系数取 0.05～0.1；最后在求解过程中监控沙漏能、内能、动能和质量增加曲线。

判断结果是否可信，可以看几个指标：沙漏能与内能之比是否低于 5%；质量增加是否在允许范围内；总能量是否守恒；接触力是否出现异常振荡。如果沙漏能持续上升，说明沙漏控制不足或网格过粗；如果质量增加过大，说明质量缩放过度，惯性效应可能已经影响结果。对这些指标的检查，应贯穿整个求解过程，而不是只看最终变形形状。