4. 显式动力学有限元方法原理

跌落冲击是典型的瞬态动力学问题：PCB 组件在几毫秒内经历接触、反弹、应力波传播与局部塑性变形，位移、速度和加速度都随时间剧烈变化。隐式方法需要每步求解大型线性方程组并迭代收敛，在强非线性接触下容易不收敛；显式动力学则不做整体刚度矩阵求逆，而是用中心差分法在时间域上逐步推进，天然适合短时、高非线性、含接触的冲击问题。本节围绕中心差分法、稳定时间步长、质量缩放、沙漏控制四条主线，说明显式方法在跌落分析中的适用边界。

一、中心差分法：显式推进的基本格式

显式动力学以节点加速度为核心变量。在时刻 t_n，由运动方程 M a_n = F_n^ext − F_n^int 直接得到加速度：

a_n = M^(-1) · (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 为质量矩阵。显式方法通常采用集中质量矩阵（lumped mass），即把单元质量按节点集中，使 M 成为对角阵，求逆退化为逐节点除法，无需矩阵分解。内力 F_n^int 由当前构型的单元应力积分得到，接触力也作为外力的一部分加入。速度采用半步长存储（v_n+1/2），位移与加速度在整步长上定义，这种交错格式即中心差分法，具有二阶精度。

关键点在于：显式方法不需要迭代，每一步的计算量近似与单元数成正比，但代价是条件稳定——时间步长不能任意取大。

二、稳定时间步长：Courant 条件

显式积分的稳定性由 Courant–Friedrichs–Lewy 条件约束。对无阻尼线性系统，稳定时间步长满足：

Δt ≤ Δt_crit = 2 / ω_max
Δt_crit ≈ L_e / c,  c = sqrt(E / ρ)

其中 ω_max 为系统最高固有频率，L_e 为最小单元特征长度，c 为材料声速，E 为弹性模量，ρ 为密度。物理含义是：应力波在一个时间步内传播的距离不能超过最小单元尺寸，否则数值解会发散。实际求解器通常取安全系数 0.9 左右，即 Δt = 0.9 L_e/c。

以跌落分析中常见的焊球为例：SnAgCu 焊料 E ≈ 40 GPa，ρ ≈ 7400 kg/m³，则 c ≈ 2325 m/s。若最小单元尺寸为 20 μm，则 Δt_crit ≈ 8.6 ns。若跌落总时长为 5 ms，需要约 58 万步。这个量级说明：显式跌落分析的计算成本主要由最小单元尺寸决定，网格中一个过小的单元会拖慢整个模型。

三、质量缩放：用可控误差换计算效率

质量缩放（mass scaling）通过人为增大密度或质量，降低声速，从而放大稳定时间步长。若希望把时间步长从 Δt 提高到 Δt'，缩放因子为 f = (Δt'/Δt)²，即密度乘以 f。例如把时间步长提高 4 倍，密度需增大 16 倍。

质量缩放不是无代价的：它改变了惯性，会影响冲击响应的高频成分和动能。工程上通常采用两种策略：一是对个别过小单元做局部缩放，只影响局部；二是全局缩放并监控质量增加比例，一般控制在 1%～5% 以内。判断是否可接受，可对比缩放前后关键位置（如焊球、铜柱）的峰值应力和加速度，若差异在工程允许范围内，则认为缩放合理。

四、沙漏控制：一阶减缩积分的副作用

显式分析为降低计算量，常采用一阶减缩积分单元（如 C3D8R）。减缩积分只有一个积分点，会产生零能模式，即沙漏模式：单元在无应变能增加的情况下发生非物理变形，表现为网格出现棋盘状或锯齿状畸变。沙漏能若不受控，会污染结果甚至使计算失稳。

常用控制手段包括：引入沙漏刚度（hourglass stiffness）、采用物理稳定化或人工阻尼。判断沙漏是否受控，可查看沙漏能与内能之比，一般要求小于 5%，严格时小于 1%。若沙漏能持续增长，应加密网格、改用全积分单元或调整沙漏控制系数。

五、跌落分析适用性：为什么显式方法合适

跌落冲击的物理特征与显式方法的优势高度匹配。第一，持续时间短，通常 1～10 ms，显式步数虽多但每步极快，总耗时可控。第二，接触非线性强，焊球与 PCB、外壳与地面之间发生碰撞、滑移和分离，隐式方法需要反复迭代接触状态，显式方法则把接触力直接并入外力，无需收敛判断。第三，材料非线性明显，焊料在冲击下进入塑性甚至蠕变，显式方法对材料本构的更新是逐单元局部计算，不增加全局求解难度。

但显式方法并非万能。若关注的是跌落后的残余应力或长期蠕变，显式积分得到的是瞬态响应，需要再转入隐式或粘塑性分析。若模型以准静态变形为主，显式方法反而效率低。因此，判断是否采用显式方法，可依据三条：载荷持续时间是否远小于结构自振周期、是否含强接触或大变形、是否以瞬态峰值响应为主要目标。

六、参数设置示例与检查要点

以下给出一组用于教学演示的显式跌落分析参数设置示例（数值为合理假设，非实测结果）：

设置完成后，建议按以下顺序检查：先确认最小单元位置与尺寸，估算稳定步长；再决定是否质量缩放，并记录质量增加比例；然后运行短时预分析，观察沙漏能与内能比值是否稳定；最后对比关键节点（焊球顶部、PCB 板边）的加速度和应力峰值，确认结果对步长和缩放不敏感。若沙漏能持续上升或质量增加超过 5%，应回到网格与缩放策略重新调整，而不是直接接受结果。

综上，中心差分法提供了无需迭代的时间推进框架，稳定时间步长划定了计算可行域，质量缩放和沙漏控制分别从效率和精度两端对方法进行调节。理解这四者的相互关系，才能在跌落冲击分析中合理设置参数、正确判读结果。