1. 接触碰撞工程问题与LS-DYNA求解定位

接触碰撞是工程中最常见、也最容易被低估的一类力学问题。它看起来只是“两个物体撞在一起”，但真正决定仿真成败的，是碰撞过程中极短时间内的应力波传播、材料非线性、几何大变形以及接触界面的不断变化。汽车正面碰撞、手机跌落、电池包挤压、护栏撞击、包装件运输冲击，这些场景的共同点是：载荷持续时间通常在毫秒量级，结构响应却包含弹性、塑性、断裂、摩擦、分离与再接触等多种机制。用静力学或线性模态方法去处理，往往连基本趋势都难以给出。

本节先建立三个基本认识：接触碰撞问题到底难在哪里；显式动力学为什么适合这类问题；LS-DYNA在碰撞仿真链条中处于什么位置。理解这三点，后续的接触定义、材料模型、沙漏控制、能量检查和结果判读才有落脚点。

1.1 接触碰撞的典型工程场景

从工程对象看，接触碰撞问题大致可以分为四类。第一类是整车与交通碰撞，例如正面100%重叠刚性壁障碰撞、40%偏置可变形壁障碰撞、侧面移动壁障碰撞和行人保护头部冲击。这类问题关注乘员舱完整性、B柱侵入量、加速度峰值和假人损伤指标。第二类是消费电子与电池跌落，例如手机从1.2 m高度跌落到混凝土表面、电池包底部球击或挤压。这类问题关注局部塑性铰、焊点失效、电解液泄漏风险和内部短路可能性。第三类是防护与吸能结构，例如公路护栏、列车吸能元件、头盔缓冲层和物流包装。这类问题关注吸能量、峰值力传递率和回弹行为。第四类是制造与装配接触，例如冲压成形中的板料与模具接触、螺栓预紧、过盈配合压装。这类问题虽然速度较低，但接触状态复杂，同样需要稳定的接触算法。

这些场景的共性可以归纳为：高速度、短历时、强非线性、多体接触、能量转换剧烈。以一辆1.5 t轿车以50 km/h撞击刚性墙为例，初始动能约为144 kJ，碰撞持续时间约80～120 ms，平均减速度可达20g以上。仿真若不能在毫秒级时间步内稳定推进，就无法捕捉峰值力和侵入量的真实波形。

1.2 显式动力学为什么适合碰撞仿真

显式动力学求解的是动力学方程 M a + C v + F_int = F_ext，其中质量矩阵M通常采用集中质量形式，因此不需要组装和求解整体刚度矩阵，也不需要进行迭代收敛判断。每一步的时间步长由Courant条件控制：Δt ≤ L_min / c，其中L_min是最小单元特征长度，c是材料声速。对钢材料，声速约为5000 m/s；若最小单元尺寸为5 mm，则临界时间步长约1 μs。这意味着1 ms的物理过程需要约1000步计算，但每一步计算量小、可高度并行，适合GPU和多核集群加速。

相比之下，隐式方法需要求解大型非线性方程组，在接触状态频繁切换、单元严重畸变时容易不收敛。显式方法虽然时间步长小，但不存在收敛问题，只要时间步稳定，就能持续推进。这正是碰撞、冲击、爆炸、跌落这类“瞬态、强非线性、接触主导”问题的首选求解策略。

1.3 LS-DYNA在碰撞仿真中的角色

LS-DYNA是显式动力学领域的代表性求解器，其核心优势在于：丰富的接触算法（单面接触、面面接触、点面接触、绑定接触、侵蚀接触等）、成熟的材料库（金属塑性、复合材料、泡沫、橡胶、混凝土、生物组织等）、完善的单元技术（壳、实体、梁、离散单元、SPH等）以及强大的并行计算能力。在汽车行业，LS-DYNA长期作为碰撞安全仿真的主力工具；在消费电子和电池领域，它也广泛用于跌落、挤压和冲击分析。

从仿真流程看，LS-DYNA通常处于“前处理—求解—后处理”链条的中段。前处理完成几何清理、网格划分、材料与接触定义、边界与初始条件设置，输出关键字文件；LS-DYNA读取关键字文件进行求解，输出d3plot、binout等结果文件；后处理工具读取结果，提取位移、速度、加速度、应力、接触力、能量等曲线。理解这个链条，有助于在后续章节中明确每一步的输入输出和检查重点。

1.4 课程任务对象与学习路径

本课程以“可复现的接触碰撞建模”为目标，任务对象覆盖三类典型模型：一是刚性墙—薄壁方管轴向压溃，用于理解接触定义、壳单元行为和吸能评价；二是球体—平板低速冲击，用于理解点面接触、局部塑性和回弹；三是简化电池包—刚性柱挤压，用于理解多部件接触、绑定失效和侵入量判读。这三类模型规模适中，既能在普通工作站上完成计算，又能覆盖接触碰撞建模的核心技术点。

学习路径上，建议先建立“关键字文件—求解—结果”的完整闭环，再逐步深入接触算法、材料模型、沙漏控制、时间步与质量缩放、能量检查和网格收敛性。每完成一个模型，都应回答三个问题：接触是否按预期传递力？能量是否守恒？关键响应量是否随网格和参数变化而收敛？只有把这三个问题回答清楚，仿真结果才具备工程参考价值。