1. 课程导入与跌落冲击仿真概览

跌落冲击是消费电子、家电、汽车电子、包装运输等行业最常见的可靠性工况之一。手机从口袋滑落到地面、路由器在物流分拣线上翻滚、电池包在整车碰撞中受到挤压，本质上都属于短时、高速、强非线性的瞬态动力学问题。用静力学或准静态方法去评估这类问题，往往无法捕捉应力波传播、局部塑性铰、接触分离与再接触等关键现象。LS-DYNA 正是为这类问题而生：它以显式时间积分求解瞬态动力学方程，能够在一次求解中同时处理大变形、材料非线性、接触碰撞与失效，是跌落冲击仿真领域事实上的工业标准工具。

为什么跌落冲击要用显式动力学

显式动力学求解的运动方程为 M·a = F_ext − F_int，其中 M 为质量矩阵，a 为节点加速度，F_ext 为外力，F_int 为内力。LS-DYNA 采用集中质量矩阵（lumped mass），使 M 成为对角阵，无需组装与求逆整体刚度矩阵，因此每个时间步只需做向量运算，单步计算量极小。配合中心差分法，位移、速度、加速度按时间递推：

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

显式方法的稳定条件是时间步必须小于临界步长，近似为 Δt_crit ≈ L_min / c，其中 L_min 是最小单元特征尺寸，c 是材料声速。对钢材，c 约为 5000 m/s；若最小单元尺寸为 0.5 mm，则临界步长约 0.1 μs。跌落过程通常持续 5～20 ms，因此需要数万到数十万步。这解释了显式动力学的两个典型特征：单步便宜、总步数巨大，以及网格最小尺寸直接决定计算成本。与之相对，隐式方法每步都要迭代求解大型方程组，在强接触、大变形下容易不收敛，因此跌落冲击几乎不用隐式求解。

典型工程应用场景

LS-DYNA 跌落冲击分析的工程场景大致可分为四类。第一类是消费电子整机跌落，例如手机、平板、TWS 耳机，关注屏幕玻璃破裂、中框塑性变形、内部焊点与连接器失效，通常按 1.0～1.5 m 高度、六面或多角度跌落评估。第二类是包装与物流运输，例如服务器机箱、家电整机带包装跌落，关注缓冲泡沫的吸能与产品内部加速度峰值，常用 ISTA 或 GB/T 4857 系列标准规定跌落高度与姿态。第三类是汽车电子与电池包，例如 ECU、BMS、动力电池在碰撞或跌落工况下的结构完整性与短路风险。第四类是零部件级落锤与冲击试验，例如玻璃盖板、连接器、焊球的落球冲击，用于标定材料失效参数。

这些场景的共同点是：关心的是峰值应力、峰值加速度、塑性应变、能量吸收与失效判据，而不是长期疲劳或稳态响应。因此仿真目标通常设定为“在给定跌落高度与姿态下，判断关键部位是否超过许用应力或应变阈值”。

完整仿真流程

一个可复现的跌落冲击仿真，通常按以下链路推进：几何清理与简化 → 网格划分与质量检查 → 材料本构与失效参数设定 → 接触定义 → 初始速度与重力场施加 → 求解控制与输出设置 → 求解 → 结果判读与验证。其中几何简化要删去对刚度贡献很小的圆角、螺纹、小孔，但必须保留可能成为应力集中源的倒角与卡扣；网格方面，冲击接触区应保证至少 3 层单元，最小单元尺寸与临界步长直接相关；接触常用 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE 或 *CONTACT_AUTOMATIC_SINGLE_SURFACE，并设置合理的软约束与穿透检查；初始速度通过 *INITIAL_VELOCITY 施加，方向与跌落姿态一致；求解时间一般取 5～20 ms，输出 d3plot 用于后处理。

结果判读与验证要点

求解完成后，后处理通常关注三类曲线：关键节点的加速度时程、关键单元的等效应力与塑性应变时程、系统总能量与内能/动能/沙漏能的分配。一个健康的跌落仿真，沙漏能应控制在总能量的 5% 以内，否则说明网格或接触设置存在问题。加速度峰值出现时刻往往对应第一次接触冲击，随后可能出现二次冲击，这与结构回弹再接触有关。若仿真得到的失效位置与试验照片一致、加速度峰值误差在 15% 以内，通常认为模型具备工程可信度；否则需要回到材料失效参数、接触刚度或网格尺寸进行修正。

本课程后续将沿着上述链路逐项展开：从几何与网格、材料本构、接触设置，到求解控制、结果判读与试验对标，最终形成一套可复现的跌落冲击仿真工作流。掌握显式动力学的基本假设与临界步长概念，是理解后续所有操作的前提。