2. 跌落冲击力学基础
跌落冲击仿真与静力分析最大的区别在于：载荷不是缓慢施加的，而是在毫秒甚至微秒量级内突然作用到结构上。理解冲击动力学的基本概念，是后续正确设置 ANSYS 显式动力学（Explicit Dynamics / LS-DYNA）模型、选择材料本构、判断结果合理性的前提。本节围绕冲击动力学基本概念、应力波传播、惯性效应、应变率效应、材料在冲击下的行为特点，以及跌落高度与冲击速度的换算关系展开。
一、冲击动力学的基本特征
冲击问题有三个显著特征：载荷时间极短、惯性力不可忽略、应力以波的形式传播。在静力分析中，我们假设结构内部各点同时达到平衡，惯性项被忽略；而在冲击中，结构各点的运动状态并不同步，必须求解包含惯性项和阻尼项的动力学方程：
M·ü + C·u̇ + K·u = F(t)其中 M 为质量矩阵，C 为阻尼矩阵，K 为刚度矩阵，F(t) 为随时间快速变化的冲击载荷。跌落冲击中，F(t) 往往不是直接给定的，而是通过接触碰撞产生的，因此还需要引入接触算法和罚函数或拉格朗日乘子法来处理。
冲击问题通常用两个时间尺度描述：载荷持续时间和结构自振周期。当载荷持续时间远小于结构自振周期时，结构来不及整体响应，局部区域先发生变形，这就是所谓的“局部冲击响应”；反之则接近准静态响应。跌落仿真中，手机、平板等消费电子产品的跌落脉宽常在 0.5～3 ms，而整机一阶模态周期可能在 5～20 ms，因此局部响应占主导。
二、应力波传播
冲击载荷作用下，应力和应变不会瞬间传遍整个结构，而是以波的形式从加载点向外传播。在一维杆中，弹性纵波波速为：
c = √(E/ρ)其中 E 为弹性模量，ρ 为密度。以钢材为例，E ≈ 210 GPa，ρ ≈ 7850 kg/m³，则 c ≈ 5170 m/s；铝合金 E ≈ 70 GPa，ρ ≈ 2700 kg/m³，c ≈ 5090 m/s；而 ABS 塑料 E ≈ 2.3 GPa，ρ ≈ 1050 kg/m³，c ≈ 1480 m/s。可见材料越硬越轻，波速越高。
应力波在结构中传播时会遇到界面、自由表面和截面变化，产生反射与透射。当压缩波从构件内部传到自由表面时，会反射为拉伸波，这种拉伸波是导致塑料外壳、玻璃盖板背面开裂的重要原因。理解这一点，对判断跌落仿真中裂纹萌生位置非常关键。
三、惯性效应
冲击过程中，结构各部分的加速度不同，惯性力会显著改变应力分布。以一块均质平板为例，当冲击点位于中心时，中心区域先获得速度，边缘区域由于惯性仍保持静止，于是板内产生弯曲波和膜应力。惯性效应在以下情况尤为突出：
- 结构质量较大、刚度相对较低时，整体惯性响应明显；
- 冲击速度高、脉宽短时，局部惯性力远大于外载；
- 存在集中质量（如电池、摄像头模组）时，局部惯性力会放大连接处的应力。
在 ANSYS 显式动力学中，惯性效应是自动包含的，因为求解的是完整运动方程。但需要注意：如果为了节省计算时间而人为放大质量（质量缩放），会改变惯性效应，必须控制缩放比例，一般建议不超过 1%～5%，并监控动能与内能之比。
四、应变率效应
材料在冲击下的力学行为与准静态下明显不同，最典型的是应变率强化：应变率越高，屈服强度和流动应力越高。常用 Cowper-Symonds 模型描述：
σ_d / σ_s = 1 + (ε̇ / D)^(1/p)其中 σ_d 为动态屈服应力，σ_s 为准静态屈服应力，ε̇ 为应变率，D 和 p 为材料常数。以低碳钢为例，D ≈ 40 s⁻¹，p ≈ 5；当应变率达到 1000 s⁻¹ 时，动态屈服强度可比准静态提高约 50%～80%。塑料的应变率敏感性更强，ABS 在 1000 s⁻¹ 下模量可提高 2～3 倍，但断裂应变下降。
跌落仿真中，冲击点附近的应变率通常在 10²～10⁴ s⁻¹ 量级，因此必须使用率相关材料模型。若仍采用准静态应力-应变曲线，会低估承载能力，导致仿真结果偏保守或与试验不符。
五、材料在冲击下的行为特点
不同材料在冲击下的表现差异很大：
| 材料类型 | 典型行为 | 仿真要点 |
|---|---|---|
| 金属（钢、铝） | 应变率强化、塑性硬化、可能绝热剪切 | Johnson-Cook 或 Cowper-Symonds + 塑性 |
| 塑料（ABS、PC） | 强应变率敏感、大变形、韧性/脆性转变 | 率相关弹性-塑性，必要时加失效准则 |
| 玻璃/陶瓷 | 脆性断裂、裂纹快速扩展、碎片飞溅 | JH-2 或最大主应力失效 + 单元删除 |
| 复合材料 | 分层、纤维断裂、基体开裂 | Hashin 准则、 cohesive 单元 |
此外，冲击下材料还可能发生绝热升温，在高应变率塑性变形中，塑性功大部分转化为热，来不及散失，导致局部软化。这一效应在高速冲击和金属切削仿真中不可忽略，但在普通消费电子跌落中影响相对有限。
六、跌落高度与冲击速度的换算
自由跌落时，忽略空气阻力，冲击速度由能量守恒得到：
v = √(2gh)其中 g = 9.81 m/s²，h 为跌落高度。常见换算如下：
| 跌落高度 h | 冲击速度 v | 典型场景 |
|---|---|---|
| 0.5 m | ≈ 3.13 m/s | 桌面跌落 |
| 1.0 m | ≈ 4.43 m/s | 手持跌落 |
| 1.5 m | ≈ 5.42 m/s | 口袋高度跌落 |
| 2.0 m | ≈ 6.26 m/s | 举高跌落 |
若已知冲击速度，也可反推等效跌落高度：h = v²/(2g)。例如某标准要求 1.2 m 跌落，则冲击速度约为 4.85 m/s。在 ANSYS 中，通常将速度作为初始条件施加到整个模型或冲击物上，而不是直接施加力。若考虑冲击物质量，还需根据动量守恒或接触刚度估算接触力峰值，但显式动力学中一般直接给定初速度即可。
掌握以上六个方面后，就可以在 ANSYS 中合理定义材料模型、施加初始速度、设置接触和求解控制。下一阶段将进入具体操作，但在此之前，建议先用手算校核冲击速度和波速量级，确保仿真设置与物理直觉一致。