1. 课程导入与高速碰撞问题界定

在工程仿真领域，“高速碰撞”并不是一个仅由速度数值定义的绝对概念，它更本质地取决于惯性效应与材料强度、结构刚度之间的相对关系。当冲击过程中惯性力、应变率效应和波传播效应开始主导结构响应时，准静态假设失效，问题就进入了显式动力学擅长的范畴。本课程讨论的对象，正是这类以短时、强非线性、大变形、接触与失效为特征的高速碰撞与冲击问题。

一、高速碰撞问题的典型对象与工况

从工程对象看，本课程覆盖的典型场景包括：汽车碰撞中的车身纵梁压溃、B柱与门槛梁弯折、保险杠吸能盒的渐进折叠；电池包在底部球击与侧柱碰中的箱体侵入与模组挤压；航空领域的叶片包容、鸟撞风挡与机翼前缘；防护工程中的弹体侵彻金属靶板、混凝土靶与复合装甲；以及消费电子跌落、包装运输冲击、头盔吸能等中小尺度问题。它们的共同点是：作用时间通常在毫秒到数十微秒量级，结构经历塑性大变形、断裂、碎片飞散与自接触，且往往伴随应变率从10²到10⁴ s⁻¹量级的跃升。

这些工况可以按“冲击物—被冲击物—响应模式”三个维度归类。例如，刚性弹体撞击薄板，主导响应是局部侵彻与花瓣开裂；刚性墙撞击薄壁管，主导响应是轴向渐进折叠；而两个柔性体对撞，则更多表现为整体弯曲与回弹。明确工况归类，直接决定了后续单元类型、材料模型与接触算法的选择。

二、LS-DYNA的求解器定位

LS-DYNA的核心是显式时间积分。它求解运动方程 M a = F_ext − F_int，采用中心差分法在时间步上推进，无需组装整体刚度矩阵，也不进行迭代求解。这使它天然适合处理接触状态频繁变化、材料非线性强烈、结构可能发生大变形甚至破碎的问题。与之相对，隐式方法在强非线性和接触突变时容易不收敛，而显式方法只要时间步足够小，就能稳定推进。

但显式方法并非没有代价。它的稳定性受Courant条件约束：临界时间步大致为最小单元特征长度除以材料波速。网格越细、材料波速越高，时间步越小，计算成本越大。因此，LS-DYNA的使用者必须同时理解物理问题和数值约束，而不是把模型丢给求解器就等待结果。

三、显式动力学基本假设与求解流程

一个完整的高速碰撞显式分析，通常包含以下关键环节：几何清理与简化、单元类型选择（壳、实体、梁或SPH等无网格方法）、材料本构与失效准则定义、接触定义、边界与初始速度施加、时间步控制、沙漏控制、输出设置以及结果判读。每一步都会影响最终结果的可信度。

以汽车纵梁碰撞为例，若采用壳单元，需要关注厚度方向积分点数量；若材料使用*MAT_PLASTIC_KINEMATIC，则需给出屈服强度、切线模量和应变率参数C、P；接触常用*CONTACT_AUTOMATIC_SINGLE_SURFACE处理自接触；初始速度通过*INITIAL_VELOCITY施加。时间步由程序自动计算，但可通过*CONTROL_TIMESTEP设置质量缩放或时间步缩放系数。这些关键字不是孤立存在的，它们共同构成一个可复现的分析流程。

四、课程任务路线与学习目标

本课程的任务路线可以概括为“建模—求解—判读—优化”四个阶段。第一阶段，学会把工程问题抽象为可计算的几何、单元、材料与接触定义；第二阶段，掌握显式求解的关键控制参数，理解时间步、沙漏、质量缩放对结果的影响；第三阶段，能够从能量曲线、接触力、变形模式中判断结果是否可信；第四阶段，基于判读结论调整模型，形成可复现的分析闭环。

具体到本课程，你将逐步接触以下内容：薄壁管轴向压溃的壳单元建模与材料参数标定；刚性弹侵彻金属靶的SPH与侵蚀接触设置；电池包底部球击的接触与失效定义；以及如何用能量平衡和网格收敛性验证结果。每一章都会给出可操作的参数示例和判断依据，而不是停留在概念介绍。

五、一个可参考的入门示例参数

为了让你对后续内容有直观预期，这里给出一个教学示例：一根长度300 mm、截面40×40 mm、壁厚1.5 mm的铝合金薄壁方管，以10 m/s轴向撞击刚性墙。材料采用*MAT_PLASTIC_KINEMATIC，密度2.7e-9 t/mm³，弹性模量70 GPa，泊松比0.33，屈服强度200 MPa，切线模量800 MPa，应变率参数C=6500、P=4。壳单元厚度方向5个积分点，接触采用自动单面接触，初始速度通过*INITIAL_VELOCITY施加。计算时间设为20 ms，时间步由程序自动控制。预期结果是方管发生渐进折叠，接触力峰值在30–50 kN量级，总吸能量约300–500 J。这些数值属于合理假设的示例，实际结果会因网格尺寸、材料参数和接触设置而变化。

通过这个示例，你可以看到：高速碰撞分析不是简单地“跑一个模型”，而是需要把物理判断、数值约束和工程目标结合起来。接下来的章节，将围绕这条主线逐步展开。