📚 课程目录第3章 LS-DYNA碰撞建模基础与参数化起点

3. LS-DYNA碰撞建模基础与参数化起点

在参数化碰撞性能优化中，真正决定优化上限的往往不是优化算法本身，而是建模阶段是否把“可调变量”识别清楚。LS-DYNA显式动力学模型由大量关键字卡片构成，其中一部分参数对碰撞响应高度敏感，另一部分则几乎不影响结果。本节的目标是建立一套判断标准：哪些几何尺寸、材料属性、接触参数和沙漏控制参数适合作为参数化变量，哪些应当固定，以及这些变量如何通过关键字卡片被外部脚本驱动。

显式动力学关键设置：先保证模型“算得动、算得准”

碰撞属于高度非线性瞬态问题，LS-DYNA使用中心差分法进行显式时间积分。时间步长由最小单元尺寸和材料波速共同决定，稳定条件可近似写为 Δt ≤ L_min/c，其中 c = √(E/ρ)。这意味着网格细化会直接压缩时间步长，进而增加计算成本。因此在参数化之前，必须确认 *CONTROL_TIMESTEP 的设置方式：是否启用质量缩放（DT2MS）、缩放是否只作用于小于目标步长的单元、缩放质量占比是否控制在合理范围（工程上常关注总质量增加不超过1%～5%）。

同时，*CONTROL_TERMINATION 的终止时间应覆盖整个碰撞过程并留出余量，*DATABASE_BINARY_D3PLOT 的输出间隔要与后续能量、侵入量、加速度曲线的提取需求匹配。若输出过密，D3PLOT文件会迅速膨胀；过疏则无法捕捉峰值。一个可参考的示例配置如下：

*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
       0.0       0.9         0       0.0    -1.0e-6         0         0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
     0.120         0       0.0       0.0       0.0
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid
    0.0020         0         0         0         0

其中 DT2MS 取负值表示仅对时间步长小于该值的单元进行质量缩放，这是碰撞建模中较常见的做法。这些控制参数通常不作为优化变量，但它们是参数化研究能否稳定运行的前提。

材料模型：参数化最活跃的区域

碰撞仿真中金属薄壁结构常用 *MAT_PIECEWISE_LINEAR_PLASTICITY（MAT_024），它通过弹性模量 E、泊松比 ν、密度 ρ、屈服强度 σ_y、切线模量 E_t 以及有效塑性应变失效值 ε_f 描述弹塑性行为。其中 σ_y、E_t、ε_f 对吸能量和峰值力影响显著，是天然的设计变量。例如在材料等级替换或热处理工艺变化的研究中，可将屈服强度设为 250～420 MPa 区间内的连续变量。

对于泡沫、蜂窝等吸能材料，常用 *MAT_CRUSHABLE_FOAM（MAT_063）或 *MAT_HONEYCOMB（MAT_026），其应力—应变曲线通过 *DEFINE_CURVE 定义。此时参数化对象可以是平台应力、密实化应变或曲线缩放系数。需要注意的是，材料曲线若通过曲线ID引用，参数化脚本应同步修改曲线数据，而不是只改材料卡片中的常数。

接触与沙漏控制：影响稳定性的两个关键开关

碰撞模型中接触定义错误会导致穿透或非物理的接触力振荡。常用 *CONTACT_AUTOMATIC_SINGLE_SURFACE 处理自接触，*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE 处理刚性墙与结构、或不同部件之间的接触。关键参数包括静摩擦系数 FS、动摩擦系数 FD、接触厚度 SST/MST 以及 SOFT 选项。摩擦系数在滑移吸能、撕裂和挤出行为中影响明显，可作为参数化变量，但取值范围应受物理约束，例如钢—钢干摩擦通常不超过0.3。

沙漏控制是显式动力学特有的问题。减缩积分单元只有一个积分点，容易产生零能模式。LS-DYNA提供 *HOURGLASS 关键字，常用类型为类型4（Flanagan-Belytschko刚度形式）或类型6（Belytschko-Bindeman）。控制系数 QH 通常取 0.03～0.10。若系数过小，沙漏能占比升高，结果不可信；过大则会使结构人为变刚。判断标准是：沙漏能应低于总内能的5%，理想情况下低于2%。因此 QH 本身也可以作为稳健性参数进行扫描，但它通常不是性能优化变量，而是模型可信度的守门参数。

识别适合参数化的几何与属性变量

几何变量方面，薄壁吸能盒的壁厚是最典型的连续变量，通常可在 1.2～2.5 mm 之间取值；截面形状参数（如多边形边数、圆角半径）属于离散或形状变量；诱导槽的深度、间距和数量则同时影响初始峰值力和变形模式，适合作为多水平离散变量。需要强调的是，几何变量修改往往涉及网格重新划分，若采用参数化建模脚本（如LS-PrePost的脚本接口或Python调用），应保证网格密度随尺寸变化而自适应，否则会出现网格畸变或单元尺寸突变。

属性变量方面，材料屈服强度和失效塑性应变对吸能影响最大，接触摩擦系数对滑移和撕裂影响明显。一个实用的筛选方法是：先做单变量扫描，观察响应变化幅度。若某变量在±10%范围内引起峰值力或吸能量变化超过5%，则值得纳入优化变量集；若变化小于1%，则固定以降低维度。以下是一份可直接参考的变量筛选记录表结构：

变量名	所属卡片	基准值	扫描范围	Fmax变化	Eabs变化	结论
壁厚 t	*SECTION_SHELL	1.8 mm	1.5～2.2 mm	+18%	+22%	纳入优化
屈服强度 σy	*MAT_024	320 MPa	280～400 MPa	+12%	+15%	纳入优化
失效应变 εf	*MAT_024	0.35	0.25～0.45	-6%	+9%	纳入优化
摩擦系数 FS	*CONTACT	0.15	0.05～0.25	+3%	+4%	可纳入
沙漏系数 QH	*HOURGLASS	0.05	0.03～0.10	<1%	<1%	固定，仅做稳健性检查

表中数值为教学示例，用于说明筛选逻辑，并非某具体项目的实测结果。实际使用时，应以自身模型的单变量扫描数据替换。

从建模到参数化的衔接要点

完成变量识别后，需要把变量映射到关键字卡片的具体字段位置。例如壁厚对应 *SECTION_SHELL 的 ELFORM 之后的厚度字段，屈服强度对应 *MAT_024 的 SIGY 字段。参数化脚本应通过字段定位或模板替换的方式修改，而不是简单字符串替换，以免误改注释或其他卡片中的相同数值。建议在脚本中维护一份“变量名—卡片—字段—单位—上下限”的映射表，作为后续优化循环的唯一数据源。

最后，在正式进入优化之前，应至少完成一次基准模型的能量平衡检查：总能量守恒、沙漏能占比、质量增加比例、接触穿透量。只有基准模型可信，参数化扫描和优化结果才有意义。这一步看似基础，却决定了后续所有参数研究的有效性。