4. 材料模型标定

材料模型标定是LS-DYNA仿真与试验对标中最基础、也最容易出问题的一环。它的本质不是“把试验数据填进卡片”，而是完成一次参数反演：以试验可观测的宏观响应（力—位移、应力—应变、加速度、吸能量等）为目标，反向确定材料本构模型中无法直接测量的内部参数，再用独立工况验证这组参数是否具备预测能力。若把标定理解为曲线拟合，往往会得到“拟合工况很准、外推工况失效”的卡片。

一、先明确标定的输入、输出与映射关系

标定的输入通常包括三类：一是标准材料试验数据，如准静态单轴拉伸的工程应力—应变曲线、压缩曲线、剪切曲线；二是试样几何与试验条件，如标距段尺寸、应变率、夹持方式；三是仿真侧的建模约定，如单元类型、壳厚积分点、是否使用*MAT_024、*MAT_123等。输出则是材料卡片中的一组参数，例如弹性模量E、泊松比ν、屈服面参数、硬化曲线、失效应变或失效应力。

关键难点在于：试验曲线是“结构响应”，材料卡片参数是“本构参数”，两者之间隔着试样几何、应力状态和数值离散。因此标定必须建立一条清晰的映射链：试验原始信号 → 工程/真应力应变 → 本构参数初值 → 单单元或试样级仿真 → 与试验曲线对比 → 参数修正。下面这张图给出了这一反演与验证链路的整体结构。

二、从试验曲线提取本构参数的实操方法

以常用的金属塑性模型 *MAT_PLASTIC_KINEMATIC（*MAT_003）或 *MAT_PIECEWISE_LINEAR_PLASTICITY（*MAT_024）为例，标定流程可以落到具体数值上。假设某低碳钢准静态拉伸试验得到：标距段初始截面积 A₀=50 mm²，标距 L₀=50 mm，屈服载荷 F_y=12.5 kN，最大载荷 F_max=18.0 kN，断裂时标距伸长 ΔL=12 mm。

第一步，计算弹性模量。在弹性段取两点（F₁=5.0 kN，ΔL₁=0.025 mm）和（F₂=10.0 kN，ΔL₂=0.050 mm），则：

σ₁ = F₁/A₀ = 5000 N / 50 mm² = 100 MPa
ε₁ = ΔL₁/L₀ = 0.025 / 50 = 5.0e-4
σ₂ = 10000 / 50 = 200 MPa
ε₂ = 0.050 / 50 = 1.0e-3
E = (σ₂ - σ₁)/(ε₂ - ε₁) = 100 MPa / 5.0e-4 = 200 GPa

第二步，确定屈服强度。σ_y = F_y/A₀ = 12500/50 = 250 MPa。若采用0.2%偏移法，则需在曲线上找到塑性应变0.002对应的应力，本例假设为260 MPa。

第三步，将工程应力—应变转换为真应力—真应变。在均匀塑性段，真应力 σ_t = σ_e(1+ε_e)，真应变 ε_t = ln(1+ε_e)。例如工程应变 ε_e=0.05 时，工程应力 σ_e=320 MPa，则 σ_t=320×1.05=336 MPa，ε_t=ln(1.05)=0.0488。塑性应变 ε_p = ε_t - σ_t/E = 0.0488 - 336/200000 ≈ 0.0471。

第四步，构造硬化曲线。将均匀塑性段离散为若干点，填入 *DEFINE_CURVE，横坐标为塑性应变，纵坐标为真应力。注意：颈缩之后的数据不能直接作为材料曲线使用，因为此时已不是均匀变形，应通过反演或外推处理。

三、参数反演中的关键判断

参数反演不是一次完成的。常见做法是先用单单元模型（如一个实体单元或壳单元）施加与试验相同的位移边界，快速对比力—位移曲线。若仿真曲线在屈服后明显高于试验，说明硬化模量偏大；若失效应变设置过小，单元会在试验断裂前被删除，表现为力突然掉零。

对于应变率相关材料，如 *MAT_024 配合 *DEFINE_CURVE 的应变率缩放，或 *MAT_003 中的 C、P 参数，需要至少两组不同应变率的试验。Cowper-Symonds 模型形式为：σ_y^d = σ_y[1+(ε̇/D)^(1/P)]。假设准静态屈服250 MPa，应变率100 s⁻¹时屈服320 MPa，则可反推 D 和 P。实际标定中通常固定 P=5 左右，再拟合 D，或通过最小二乘同时优化。

四、验证：用独立工况检验预测能力

标定完成后必须做验证，且验证工况不能与标定工况相同。推荐至少包含三类：一是不同应变率下的拉伸，检验率相关参数；二是不同几何试样，如缺口拉伸，检验应力三轴度影响；三是结构件级试验，如薄壁管轴向压溃，检验吸能量与变形模式。

验证指标建议量化：峰值力误差控制在±5%以内，吸能量误差±10%以内，失效位移误差±15%以内。若结构件仿真出现非物理的单元畸变或删除，应回到卡片检查失效应变、单元尺寸与沙漏控制，而不是继续调材料参数。

最后需要强调：材料卡片是仿真与试验之间的“翻译层”，标定的质量取决于对试验状态和本构假设的理解深度。同一批试验数据，若单元类型、积分方式或失效准则选择不同，最优参数也会变化。因此标定报告应记录试验条件、数据处理方法、仿真设置和验证误差，使参数具备可追溯性和可复现性。