4. 弹塑性基础参数标定

弹塑性材料模型的标定，本质上是把试验曲线翻译成求解器能读懂的参数。LS-DYNA 中常用的 *MAT_PLASTIC_KINEMATIC（MAT_003）、*MAT_ISOTROPIC_ELASTIC_PLASTIC（MAT_012）以及 *MAT_JOHNSON_COOK（MAT_015）等模型，都共享同一组基础参数：弹性模量 E、泊松比 ν、初始屈服强度 σy、硬化曲线（或硬化模量 Et/切线模量）、以及应变率相关参数（Cowper-Symonds 的 C、P 或 Johnson-Cook 的 A、B、n、C）。本节按“先弹性、再屈服、后硬化、最后应变率”的顺序，说明每一步的识别对象、数据来源与判断依据。

一、弹性模量与泊松比：从弹性段斜率提取

弹性模量来自应力—应变曲线初始直线段的斜率。工程做法是取 0.05%～0.25% 应变区间做线性拟合，避开试验机柔度与引伸计初始接触带来的非线性。以低碳钢为例，若在 0.001～0.003 应变区间拟合得到斜率约 205 GPa，则 E 取 205000 MPa。需要注意单位一致性：LS-DYNA 常用 mm-ton-s 单位制，此时应力单位为 MPa，密度为 t/mm³，E 直接填 205000。

泊松比通常由横向与轴向引伸计同步测量得到，ν = −ε_lat/ε_ax。金属材料一般在 0.28～0.34 之间，若无横向测量数据，可按材料手册取值，但应在标定报告中注明来源。泊松比主要影响弹性阶段的体积响应与颈缩后的横向变形，对屈服前的小变形仿真影响有限，但对薄壁结构或大变形问题不可忽略。

二、屈服面与初始屈服强度：0.2% 偏移法

金属材料没有明显的屈服平台时，采用 0.2% 塑性应变偏移法确定屈服点：从应变轴 0.002 处作一条与弹性段平行的直线，与试验曲线的交点对应的应力即为屈服强度 σy。对于有明显上下屈服点的材料（如部分退火钢），可取上屈服点或下屈服点，但必须在标定说明中明确选择，因为两者差异可达 10% 以上。

屈服面在 LS-DYNA 中由屈服准则描述。MAT_003 使用 von Mises 屈服面，只需 σy 一个参数；MAT_015 使用 Johnson-Cook 屈服面，需要 A（准静态屈服强度）、B（硬化模量）、n（硬化指数）。若试验只有单向拉伸曲线，通常假设各向同性屈服，即拉伸与压缩屈服强度相同。若材料存在拉压不对称（如某些镁合金、聚合物），则需要选用能描述 SD 效应的模型，本节不展开。

三、硬化曲线：从真应力—真应变到塑性段

硬化曲线是标定中最容易出错的部分。试验机输出的是工程应力—工程应变，必须先转换为真应力—真应变：

ε_true = ln(1 + ε_eng)
σ_true = σ_eng × (1 + ε_eng)

转换后，从总应变中减去弹性应变 ε_e = σ_true/E，得到塑性应变 ε_p。以塑性应变为横轴、真应力为纵轴，即得硬化曲线。LS-DYNA 中 MAT_003 可用两种方式输入硬化：一是直接给切线模量 Et（双线性），二是给多段塑性应力—塑性应变曲线（*DEFINE_CURVE）。对于 304 不锈钢，若试验得到 σy=250 MPa、Et=1500 MPa，则双线性近似可用；但若曲线在 5% 塑性应变后明显变缓，建议用多段曲线，至少取 8～12 个点，覆盖到仿真中可能出现的最大塑性应变。

MAT_015 则用 A+B·ε_p^n 拟合硬化段。拟合时先固定 A=σy，再对 B 和 n 做最小二乘。若拟合残差在 3% 以内，可认为硬化曲线可用；若残差偏大，应检查是否在颈缩后仍强行拟合，因为颈缩后真应力已不能由工程公式准确转换，需要借助 Bridgman 修正或数字图像相关（DIC）数据。

四、应变率相关性：Cowper-Symonds 与 Johnson-Cook

应变率敏感性通过不同应变率下的拉伸试验识别。典型做法是在 10⁻³、10⁻¹、1、10² s⁻¹ 四个量级下各做 3 次重复试验，取屈服强度平均值。以 6061-T6 铝合金为例，准静态屈服约 275 MPa，在 1000 s⁻¹ 下屈服约 320 MPa，则可用 Cowper-Symonds 模型拟合：

σ_y(ε̇) = σ_y0 × [1 + (ε̇/C)^(1/P)]

取 σy0=275 MPa，用 1000 s⁻¹ 下 320 MPa 的数据点反算，若 C=6500、P=4，则预测值为 275×[1+(1000/6500)^0.25]≈275×1.16≈319 MPa，与试验吻合。若使用 MAT_015，则应变率项为 1+C·ln(ε̇/ε̇0)，C 由不同应变率屈服强度的斜率确定。注意：应变率参数对仿真结果影响很大，但试验数据往往离散，建议至少用 3 个应变率点拟合，并检查拟合曲线是否单调。

五、标定结果的验证与常见问题

参数识别完成后，必须做一次单轴拉伸仿真反算：用 1 个单元或标准拉伸试样，施加与试验相同的位移边界，提取力—位移曲线，与试验曲线对比。判断标准是：弹性段斜率误差 ≤2%，屈服点误差 ≤5%，硬化段在 5% 塑性应变内误差 ≤8%。若弹性段偏软，通常是试验机柔度未修正；若屈服点偏高，可能是偏移法取值点选错；若硬化段后期偏差大，多半是颈缩后数据未修正。

应变率参数标定中，常见问题是只用一个高应变率点拟合 C、P，导致中间应变率预测失真。建议至少用 3 个应变率点，并检查拟合曲线在 10⁻³～10³ s⁻¹ 范围内是否单调递增。若材料对应变率不敏感（如 304 不锈钢在室温下），可令 C 取极大值或直接关闭应变率效应，避免引入虚假的率相关性。

最后需要强调：弹塑性基础参数不是一次标定就永久有效的。同一牌号材料不同批次、不同热处理状态，屈服强度和硬化曲线都可能变化。若仿真对象涉及焊接热影响区、冷作硬化区或高温工况，应分区标定，而不是用一组母材参数覆盖整个模型。