1. 课程导入与标定问题定义
在整车碰撞、电池包挤压、行人保护、跌落与冲击防护等工程仿真中,LS-DYNA 的计算结果是否可信,往往不取决于网格有多细、求解器设置有多“高级”,而取决于材料模型参数是否真实反映了所用材料在对应工况下的力学行为。材料模型标定与失效参数识别,正是把“试验室里测到的力—位移、应力—应变、断裂形貌”翻译成“求解器能读懂的卡片参数”的过程。它既是仿真链条的入口,也是误差最容易累积放大的环节。
可以把整个问题抽象成一条清晰的输入—输出链:输入是材料试验数据、试样几何、加载条件与仿真工况;中间是材料模型选择、参数初值估计、数值仿真复现、误差度量与参数迭代;输出是一组可写入 *MAT_* 关键字卡片、并能在目标工况下复现试验响应的参数集。这条链路上任何一环失真,都会让最终仿真“看起来算完了,其实不可信”。
从图中可以看到,标定不是“把试验曲线拟合一下就结束”,而是一个带反馈的闭环:参数迭代后必须回到仿真中复现试验,若误差不满足要求,就要回到模型选择或参数调整环节。同时,误差来源分布在试验、模型、数值和流程四个层面,任何一个层面失控,都会让标定结果失去工程意义。
典型工程场景与标定对象
在汽车与消费电子行业,最常见的标定对象包括:高强钢与铝合金的弹塑性+失效参数、塑料与复合材料的率相关本构、焊点与胶层的连接失效、电池包内部结构的挤压与断裂。以某款 6000 系铝合金板材为例,典型试验矩阵会包含:准静态拉伸(应变率约 0.001 s⁻¹)、中速拉伸(约 1 s⁻¹)、高速拉伸(约 100~1000 s⁻¹)、缺口拉伸(应力三轴度约 0.33~0.6)以及剪切试验(应力三轴度接近 0)。这些试验分别用于识别硬化曲线、应变率敏感性和失效曲面。
需要强调的是,标定目标必须与仿真工况匹配。如果目标工况是低速碰撞中的薄壁梁压溃,那么只做准静态光滑拉伸标定往往不够,因为压溃过程中局部应变率可达 10~100 s⁻¹,且应力状态以剪切和压缩为主。此时应优先补充剪切与缺口试验,并在 *MAT_024 或 *MAT_037 中引入率相关项,必要时配合 *MAT_ADD_EROSION 或 GISSMO 失效模型。
输入输出与参数识别的基本逻辑
标定的输入通常以“力—位移”或“工程应力—工程应变”形式给出。以拉伸试验为例,试验机记录的是力 F 与位移 d,需要先转换为工程应力 σ_e = F/A₀ 和工程应变 ε_e = ΔL/L₀,再转换为真应力 σ_t = σ_e(1+ε_e) 和真应变 ε_t = ln(1+ε_e)。这一步看似基础,但标距选择、颈缩修正和断裂点判定会直接影响后续硬化曲线形状。
对于弹塑性模型,屈服强度 σ_y、硬化模量 E_h 或硬化曲线、应变率参数 C 与 P(Cowper-Symonds 形式)是核心待识别量。失效参数则可能包括等效塑性应变阈值、应力三轴度相关曲线、Lode 角参数以及损伤演化指数。识别过程通常先固定失效参数、拟合硬化段,再固定硬化参数、拟合失效点,最后联合微调。
误差来源与判断依据
判断一组标定参数是否可用,不能只看拟合曲线是否“好看”。更可靠的依据是:在独立工况(未参与标定的试验或仿真)中,力—位移峰值误差、能量吸收误差和断裂位置是否在可接受范围内。工程上常以峰值力误差小于 10%、吸能误差小于 15%、断裂位置与试验一致作为初步通过标准;若用于安全件,还需进一步收紧。
常见的失败模式包括:硬化曲线外推过度导致大应变段应力虚高;失效参数只拟合了单轴拉伸,导致剪切工况下过早或过晚断裂;网格尺寸变化后失效应变未做正则化,结果随网格明显漂移。这些问题在标定阶段若未被识别,进入整车模型后会以“局部穿透”“异常飞溅”“吸能不足”等形式暴露出来,排查成本极高。