2. 跌落冲击问题力学基础

跌落冲击与准静态压缩最本质的区别在于：载荷在毫秒甚至微秒量级内完成施加，惯性效应不可忽略，应力以波的形式在结构中传播，材料的屈服强度还会随加载速率升高而改变。如果仍用静力学思路去设置跌落仿真，往往会得到偏乐观的变形量或错误的失效位置。本节围绕冲击动力学基本方程、应力波传播、应变率效应以及能量与动量守恒四条主线，说明它们在 LS-DYNA 跌落模型中的具体体现。

一、冲击动力学基本方程

连续介质在冲击下的运动由三组方程共同描述。第一是动量方程（运动方程）：

ρ · a = ∇·σ + ρ · b

其中 ρ 为密度，a 为加速度，σ 为柯西应力张量，b 为单位质量体力。与静力平衡方程 ∇·σ + ρb = 0 相比，多出的 ρ·a 项正是惯性力的来源。跌落中物体以初速度撞击地面，接触瞬间速度在极短时间内被改变，加速度可达 10³～10⁵ g，惯性项与应力梯度项同量级，因此不能略去。

第二是质量守恒（连续性方程），在拉格朗日描述下体现为密度随体积变形而变化：ρ = ρ₀ · J⁻¹，J 为变形梯度行列式。第三是能量守恒方程，即外力功等于动能、内能与耗散能的变化之和。LS-DYNA 采用显式中心差分法求解动量方程，时间步长受 Courant 条件限制：Δt ≤ L_min / c，其中 L_min 是最小单元特征尺寸，c 是材料声速。这解释了为什么跌落模型中过小的单元会显著拉长计算时间。

二、应力波传播

撞击瞬间，接触面产生的扰动以应力波形式向内部传播。一维弹性杆中纵波波速为：

c = sqrt(E / ρ)

以钢材为例，E ≈ 210 GPa，ρ ≈ 7850 kg/m³，可得 c ≈ 5170 m/s。若构件厚度为 5 mm，波穿过一次仅需约 1 μs。这意味着在跌落持续的数毫秒内，应力波会在构件内部反复反射数十甚至上百次，形成复杂的波叠加。波在自由表面反射时，压缩波变为拉伸波，这种反射拉伸往往是脆性材料或薄壁结构背面开裂的直接原因。

在 LS-DYNA 中，应力波传播的准确性依赖网格密度：要分辨一个波长，通常需要沿传播方向布置 10 个以上单元。若单元尺寸过大，波前会被数值抹平，峰值应力被低估。因此对薄壁、细长结构，网格尺寸常需控制在 0.5～1 mm 量级。

三、应变率效应

材料在高速加载下的屈服强度高于准静态值，这一现象称为应变率效应。工程中常用 Cowper-Symonds 模型描述：

σ_d / σ_s = 1 + (ε̇ / D)^(1/p)

其中 σ_d 为动态屈服应力，σ_s 为静态屈服应力，ε̇ 为应变率，D 与 p 为材料常数。对低碳钢，典型值 D ≈ 40 s⁻¹，p ≈ 5。当 ε̇ = 1000 s⁻¹ 时，动态屈服强度约为静态的 1.5～2 倍。跌落冲击中，接触区应变率常达 10²～10⁴ s⁻¹，若忽略应变率效应，会高估塑性变形、低估结构承载能力。

在 LS-DYNA 中，可通过 *MAT_PLASTIC_KINEMATIC 的 C 和 P 参数（对应 Cowper-Symonds 的 D 和 p）引入应变率强化，或使用 *MAT_JOHNSON_COOK 配合 *DEFINE_CURVE 定义应变率相关屈服。需要注意的是，应变率是单元级变量，网格尺寸会影响计算得到的应变率幅值，因此材料参数应与网格尺度匹配标定。

四、能量守恒与动量守恒的体现

跌落全过程可用能量守恒描述：初始动能 E_k 在撞击后转化为弹性应变能、塑性耗散能、接触摩擦耗散能以及数值沙漏能。LS-DYNA 的 glstat 文件给出各能量分量随时间的变化，是判断模型是否可信的关键依据。经验上，沙漏能与内能之比应控制在 5% 以内，总能量相对初始值的偏差应小于 1%。

动量守恒则体现在接触碰撞中。以质量为 m、初速度为 v₀ 的物体撞击刚性地面为例，若完全弹性反弹，反弹速度约为 v₀；若完全塑性，反弹速度接近零。实际跌落中，恢复系数 e 介于 0 与 1 之间，可通过接触定义中的 SOFT、PENALTY 等参数间接影响。对于自由跌落，初始动能由跌落高度决定：

E_k = m · g · h

例如 0.5 kg 的手机从 1.5 m 高度跌落，初始动能约 7.35 J，撞击速度约 5.42 m/s。这些数值可直接作为 *INITIAL_VELOCITY 的输入，或通过设置重力场与初始高度让求解器自行计算。

五、四条主线在建模中的对应关系

把上述力学基础落到 LS-DYNA 建模动作上，可以形成如下对应：动量方程对应显式时间积分与质量缩放控制；应力波传播对应网格尺寸与声速约束；应变率效应对应材料本构中的率相关参数；能量守恒对应 glstat 能量审计与沙漏控制。理解这层对应关系，才能在出现异常结果时快速定位是物理设置问题还是数值问题。