1. 课程导入与吸能仿真概览
薄壁结构吸能仿真，研究的不是“零件会不会坏”，而是“撞击发生时，结构如何按设计意图逐步变形、把动能转化为可控的耗散能量”。汽车前纵梁、保险杠横梁、电池包边框、轨道车辆端部吸能管、直升机座椅下方的波纹梁，都是典型对象。它们的共同特征是：用较薄的板或管，在轴向或侧向载荷下产生稳定的塑性折叠、撕裂或卷曲，从而在有限行程内吸收大量能量，同时把峰值力控制在乘员或设备可承受的范围内。
以一根常见的方形前纵梁为例：材料为DP590高强钢，截面约80 mm×80 mm，壁厚1.8 mm，长度约420 mm，前端布置诱导槽。整车正面碰撞中，它需要在约250 mm的压溃行程内吸收约35～55 kJ能量，同时把初始峰值力压到约120 kN以下。这个“能量—行程—峰值力”三角关系，就是吸能仿真的核心矛盾：吸能总量要够，压溃要稳定，峰值不能过高。LS-DYNA之所以成为该领域主流工具，正是因为它擅长处理大变形、接触、材料失效和显式时间积分下的瞬态冲击问题。
1.1 为什么薄壁结构适合吸能
薄壁结构吸能的本质，是利用金属薄壁在压缩过程中形成塑性铰和折叠波，把撞击动能转化为塑性功、摩擦功和少量断裂能。与实心结构相比，薄壁管在相同质量下能提供更长的压溃行程和更稳定的力—位移平台。理想轴向压溃曲线通常表现为：初始弹性上升、峰值力、然后进入一段相对平稳的波动平台，最后因密实化而力值快速上升。设计者希望这段平台尽量长、波动尽量小、峰值尽量低。
但实际结构很少是理想轴向压溃。前纵梁在整车中往往同时承受轴向压缩、弯曲和剪切，诱导槽位置、焊点间距、材料应变率敏感性、板厚公差都会改变折叠模式。若诱导不足，可能从端部直接失稳成欧拉弯曲，吸能量骤降；若诱导过强，又可能提前撕裂，行程未用尽就失效。因此，仿真必须能复现“起皱—折叠—压实”的完整过程，而不是只算出一个最大应力。
1.2 LS-DYNA在碰撞与吸能分析中的定位
LS-DYNA采用显式中心差分时间积分，天然适合求解接触、大变形和材料失效。对于薄壁吸能问题，常用单元为壳单元，配合Belytschko-Tsay或全积分壳公式；材料模型常用*MAT_PIECEWISE_LINEAR_PLASTICITY（024）、*MAT_JOHNSON_COOK（015）或*MAT_MODIFIED_JOHNSON_COOK，失效可用等效塑性应变或*MAT_ADD_EROSION控制。接触方面，自接触常用*AUTOMATIC_SINGLE_SURFACE，刚性墙与结构之间用*RIGIDWALL_PLANAR或*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE。
与隐式方法相比，显式方法不需要组装整体刚度矩阵，单步计算便宜，但时间步受最小单元尺寸和声速限制。壳单元稳定时间步约为Δ t ≈ Lc/c，其中Lc为特征长度，c为材料声速。对钢材料，c约5000 m/s；若最小单元尺寸为2 mm，则稳定时间步约0.4 μs。要模拟100 ms的碰撞过程，需要约25万步。这也是为什么吸能仿真中网格尺寸、质量缩放和子循环设置会直接影响计算成本与结果可信度。
1.3 课程目标与学习路径
本课程的目标不是让读者“会点按钮”，而是能独立完成一个薄壁吸能结构的建模、求解、判读和优化闭环。学习路径按六个层次推进：
- 第一层：对象与指标。明确吸能结构类型、压溃行程、吸能量、峰值力、比吸能（SEA）等指标定义。例如，比吸能可写为SEA = Eabsorbed/m，单位为kJ/kg。
- 第二层：几何与网格。掌握壳单元建模、诱导槽处理、网格尺寸对折叠模式的影响，理解沙漏控制与全积分壳的取舍。
- 第三层：材料与失效。掌握*MAT_024、*MAT_015等模型参数含义，理解应变率效应、失效应变和断裂能对压溃曲线的影响。
- 第四层：接触与加载。掌握自接触、刚性墙、初速度、附加质量、约束和边界条件的设置方法，避免初始穿透和接触失效。
- 第五层：求解与稳定性。理解显式时间步、质量缩放、子循环、能量平衡检查，能判断结果是否可信。
- 第六层：结果判读与优化。能提取力—位移曲线、内能曲线、变形模式，识别峰值力、平台力、吸能量和失效位置，并据此调整壁厚、诱导槽和材料分布。
1.4 一个可复现的入门示例
为便于后续章节衔接，这里给出一个入门级示例设置。对象为单根方形薄壁管，前端固定刚性墙，后端施加初速度。示例参数如下：
| 项目 | 设置 | 说明 |
|---|---|---|
| 单元类型 | 壳单元，Belytschko-Tsay | 厚度方向5个积分点 |
| 网格尺寸 | 2 mm | 诱导槽附近局部加密至1 mm |
| 材料 | *MAT_PIECEWISE_LINEAR_PLASTICITY | 密度7.85e-9 t/mm³，弹性模量210 GPa，泊松比0.3 |
| 塑性曲线 | 屈服强度590 MPa，切线模量1.2 GPa | 示例值，实际需按材料卡输入 |
| 失效 | 等效塑性应变0.35 | 示例阈值，用于观察撕裂 |
| 接触 | *AUTOMATIC_SINGLE_SURFACE | 静摩擦0.15，动摩擦0.1 |
| 刚性墙 | *RIGIDWALL_PLANAR | 法向沿轴向，固定 |
| 初速度 | 15 m/s | 沿轴向指向刚性墙 |
| 附加质量 | 300 kg | 通过*ELEMENT_MASS施加于后端节点 |
| 终止时间 | 30 ms | 确保压溃行程充分 |
| 输出 | 力—位移、内能、动能、沙漏能 | 通过*DATABASE_BINARY_D3PLOT和*DATABASE_GLSTAT |
提交计算后，首先检查能量平衡：总能量应基本守恒，沙漏能应低于总内能的5%～10%。若沙漏能过高，说明壳单元沙漏控制不足或网格过粗。然后提取刚性墙接触力与位移，绘制力—位移曲线。若曲线出现异常尖峰，通常来自初始穿透或接触刚度设置不当；若平台力过低且变形模式为整体弯曲，说明诱导不足或边界条件不合理。
1.5 常见误区与判断依据
吸能仿真中最常见的误区有三个。第一，把“最大应力”当作吸能指标。吸能看的是力—位移曲线下面积，不是某一点应力峰值。第二，忽略网格尺寸对折叠波长的影响。折叠波长与壁厚和材料硬化相关，网格过粗会人为抑制折叠，导致平台力偏高或模式失真。第三，不做能量平衡检查。若沙漏能、界面滑移能或质量缩放附加能量占比过大，结果即使“看起来合理”也不可信。
判断结果是否可信，可以按以下顺序检查：总能量是否守恒；沙漏能是否低于内能的10%；质量缩放增加的质量是否低于总质量的1%～2%；接触力是否出现非物理尖峰；变形模式是否与试验或经验一致；吸能量和峰值力是否落在设计目标区间。只有这些检查通过，后续的参数优化才有意义。