1. 课程导入与显式动力学认知

在工程仿真领域，求解器大致分为隐式与显式两大阵营。理解两者的差异，是学习 LS-DYNA 的第一道门槛。隐式方法求解的是形如 K·u = F 的平衡方程，其中 K 为刚度矩阵。它需要在每个时间步组装并求解一个大型线性方程组，因此单步计算成本高，但时间步可以取得很大，适合静力、模态、低频振动等准静态或缓变问题。显式方法则完全不同，它求解的是运动方程 M·a = F_ext − F_int，通过中心差分法在时间上推进：已知当前时刻的位移与速度，直接算出加速度，再积分得到下一时刻的速度与位移。整个过程不需要组装和求逆整体刚度矩阵，单步成本极低，但时间步受稳定性条件限制，必须非常小。

这个稳定性条件就是著名的 Courant–Friedrichs–Lewy（CFL）条件。对于一维杆单元，临界时间步可近似写为 Δt ≤ L / c，其中 L 是最小单元特征长度，c 是材料中的声速。以钢材为例，弹性模量约 210 GPa，密度约 7850 kg/m³，声速 c = √(E/ρ) ≈ 5170 m/s。若最小单元尺寸为 5 mm，则临界时间步约为 0.005 / 5170 ≈ 9.7×10⁻⁷ s，即约 1 微秒。这意味着模拟 10 毫秒的碰撞过程，需要约一万个时间步。显式方法正是用"小步快跑"的方式，换取每一步的廉价计算，从而高效处理短时、强非线性、大变形问题。

LS-DYNA 是显式动力学领域最具代表性的通用求解器。它最初由 Lawrence Livermore 国家实验室开发，专注于爆炸、冲击、碰撞等瞬态动力学问题，后来逐步扩展到多物理场耦合。其核心定位可以概括为：以显式时间积分为主，擅长处理接触、材料非线性、几何非线性高度耦合的短时瞬态事件。与隐式软件（如 Abaqus/Standard、ANSYS Mechanical）相比，LS-DYNA 的优势不在于静力精度，而在于能稳定地"算得动"那些隐式方法难以收敛的极端工况。

典型应用场景包括：汽车整车与零部件碰撞安全分析（正面碰撞、侧面碰撞、行人保护）、金属冲压成形与回弹预测、跌落测试（手机、包装、家电）、爆炸冲击与防护结构设计、鸟撞与叶片包容性分析、以及生物力学中的骨骼与软组织冲击响应。这些场景的共同特征是：作用时间短（毫秒级）、变形大（应变可达数十甚至上百百分比）、接触状态频繁变化、材料可能进入塑性甚至失效断裂。

从学习路径来看，建议按"认知—操作—深入"三层推进。第一层是建立直觉：理解显式积分为什么快、为什么稳、代价是什么。第二层是掌握 LS-DYNA 的基本工作流：前处理（几何清理、网格划分、材料与接触定义）、求解设置（时间步控制、沙漏控制、输出频率）、后处理（应力云图、能量曲线、接触力提取）。第三层是专题深入：材料本构（Johnson-Cook、MAT_24 等）、接触算法（罚函数法、对称罚函数法）、ALE 与 SPH 方法、多物理场耦合。

以一个简单的钢球撞击刚性墙为例，可以快速建立对显式仿真的感性认识。钢球直径 20 mm，密度 7850 kg/m³，弹性模量 210 GPa，泊松比 0.3，初始速度 10 m/s 沿 −Z 方向。刚性墙固定。网格采用六面体单元，最小尺寸约 2 mm。根据 CFL 条件，临界时间步约为 0.002 / 5170 ≈ 3.9×10⁻⁷ s，LS-DYNA 会自动计算并给出建议值，通常乘以一个小于 1 的安全系数（默认 0.9）。总求解时间设为 2 ms，则约需 5000 步。输出接触力与钢球动能曲线，可以观察到动能逐渐转化为内能，接触力先增后减，呈现典型的弹性碰撞特征。若将材料改为塑性模型，则动能不会完全恢复，部分能量被塑性耗散，接触力峰值也会降低。

最后需要明确：显式与隐式不是替代关系，而是互补工具。LS-DYNA 也具备隐式求解能力（如 *CONTROL_IMPLICIT_GENERAL），可用于显式求解前的预加载或显式结束后的回弹计算。工程实践中，常采用"隐式—显式"联合仿真：先用隐式方法完成装配预紧或重力平衡，再将结果导入显式求解器进行碰撞分析。理解两者的边界与衔接方式，是进阶应用的必备能力。