3. 汽车结构碰撞仿真材料建模

在整车碰撞仿真中，材料模型是决定吸能路径、变形模式与乘员伤害结果的核心输入。低碳钢、高强钢与铝合金在车身中承担不同角色：低碳钢用于吸能盒、纵梁等大变形区域，高强钢用于B柱、门槛梁等抗侵入区域，铝合金则用于前防撞梁、减重覆盖件。三类材料的弹塑性行为、应变率敏感性和失效机制差异显著，若统一使用默认弹塑性模型，往往导致碰撞力峰值偏低、断裂位置错误或吸能总量失真。本节围绕弹塑性本构、应变率效应、Johnson-Cook模型、失效准则与GISSMO损伤模型展开，并给出可落地的参数组织方式。

一、弹塑性本构：从屈服到硬化的基本框架

Abaqus/Explicit中金属材料常用*MATERIAL配合*ELASTIC与*PLASTIC定义。弹性段由杨氏模量E和泊松比ν描述；塑性段以“真实应力—真实塑性应变”对输入，第一行塑性应变必须为0，对应初始屈服应力。低碳钢DC04典型值：E=210 GPa，ν=0.3，屈服强度约150~180 MPa，抗拉强度约300 MPa，均匀延伸率约20%。高强钢DP780：屈服约450 MPa，抗拉约780 MPa，硬化段更陡但延伸率较低。铝合金6061-T6：E=69 GPa，ν=0.33，屈服约275 MPa，抗拉约310 MPa，硬化能力弱、断裂应变小。

需要特别注意：Abaqus塑性表使用真实应力与真实塑性应变，而材料单常给工程应力—应变。转换关系为：真实应力σ_t=σ_e(1+ε_e)，真实应变ε_t=ln(1+ε_e)，真实塑性应变ε_p=ε_t−σ_t/E。若直接输入工程值，会导致屈服后应力偏高、吸能偏大。对于颈缩后段，建议用外推或采用损伤模型接管，避免塑性段无限延伸。

二、应变率效应：碰撞速度下的强度提升

碰撞仿真应变率通常在10⁰~10³ s⁻¹。低碳钢对应变率较敏感，高强钢次之，铝合金在室温下敏感性相对较弱但不可忽略。Abaqus中可用*PLASTIC配合*RATE DEPENDENT定义屈服应力随塑性应变率的变化，常用Cowper-Symonds形式：σ_y=σ₀[1+(ε̇/D)^1/p]，其中D、p为材料常数。例如低碳钢可取D=40 s⁻¹、p=5；DP780可取D=800 s⁻¹、p=6；铝合金6061-T6可取D=6500 s⁻¹、p=4。这些数值属于教学示例，实际应以材料卡或试验标定为准。

在Abaqus中，*RATE DEPENDENT需与塑性表配合，且应变率点应覆盖准静态、10、100、1000 s⁻¹等量级。若只给一个高应变率点，中间速率靠插值，可能造成低速工况强度虚高。建议至少给出3~5个速率点，并检查单位制一致：应力MPa、密度t/mm³、时间s时，应变率单位为s⁻¹，无需额外换算。

三、Johnson-Cook模型：兼顾硬化、速率与温度

Johnson-Cook本构将流动应力写成三项乘积：σ=(A+Bε_pⁿ)(1+C ln ε̇^*)(1−T^*m)，其中ε̇^*=ε̇/ε̇₀，T^*=(T−T_r)/(T_m−T_r)。A为屈服应力，B为硬化模量，n为硬化指数，C为速率敏感系数，m为热软化指数。Abaqus中通过*PLASTIC, HARDENING=JOHNSON COOK输入A、B、n、m、T_m、T_r，再配合*RATE DEPENDENT, TYPE=JOHNSON COOK输入C与ε̇₀。以DP780教学示例：A=450 MPa，B=900 MPa，n=0.6，C=0.014，ε̇₀=1 s⁻¹，m=0.8，T_m=1500 K，T_r=293 K。铝合金6061-T6示例：A=275 MPa，B=300 MPa，n=0.35，C=0.002，m=1.0。

JC模型适合高速冲击与绝热温升明显的工况。若碰撞仿真中未开启温度自由度，热软化项不生效，此时应把m设为0或改用纯速率形式，避免参数冗余。若开启*COUPLED TEMPERATURE-DISPLACEMENT，则需同时定义比热、导热系数与Taylor-Quinney系数（常用0.9），否则温升无法反馈到流动应力。

四、失效准则与GISSMO损伤模型

碰撞仿真中，材料失效决定断裂位置与吸能截止点。常用失效准则包括等效塑性应变准则、应力三轴度相关准则和GISSMO。GISSMO（Generalized Incremental Stress State Dependent Damage Model）通过损伤变量D累积：D=∫(ε_p/ε_f(η,ε̇))dε_p，其中ε_f为失效应变，η为应力三轴度。当D达到1时触发损伤，随后进入不稳定段并最终失效删除单元。Abaqus中通过*DAMAGE INITIATION, CRITERION=GISSMO与*DAMAGE EVOLUTION配合，需输入失效曲线、损伤曲线、正则化长度等。

GISSMO的关键输入是失效应变随三轴度的变化曲线。以DP780教学示例：η=−0.33时ε_f=0.9，η=0时ε_f=0.55，η=0.33时ε_f=0.35，η=0.67时ε_f=0.22。铝合金6061-T6对应值整体更低：η=0时ε_f=0.18，η=0.33时ε_f=0.10。低碳钢DC04则较高：η=0时ε_f=0.8，η=0.33时ε_f=0.5。这些曲线应来自材料试验或标定，示例值仅用于说明输入结构。

五、参数组织与检查要点

建议按“材料卡—本构—速率—损伤—单元删除”五层组织INP。以下为DP780教学示例片段，单位制为mm、t、s、MPa：

*MATERIAL, NAME=DP780
*ELASTIC
210000., 0.3
*PLASTIC, HARDENING=JOHNSON COOK
450., 900., 0.6, 0.8, 1500., 293.
*RATE DEPENDENT, TYPE=JOHNSON COOK
0.014, 1.0
*DAMAGE INITIATION, CRITERION=GISSMO
0.55, 0.0, 0.0
*DAMAGE EVOLUTION, TYPE=ENERGY
50.0
*SECTION CONTROLS, NAME=SEC1, ELEMENT DELETION=YES
1, 1., 0.0

检查时重点看四项：一是塑性表首行塑性应变是否为0；二是速率点是否覆盖实际碰撞速率范围；三是GISSMO失效应变曲线是否与三轴度符号约定一致（拉为正、压为负）；四是单元删除后是否出现质量损失过大或能量不守恒。若碰撞力峰值偏高且断裂偏晚，优先检查失效应变是否偏大或正则化长度是否未设；若断裂偏早，则检查三轴度曲线是否过于保守。

实际项目中，建议先以准静态拉伸曲线建立基础塑性段，再用高速拉伸或SHPB数据标定速率项，最后用缺口、剪切、穿孔等试样标定GISSMO失效曲线。三类材料不可混用同一套损伤参数：低碳钢的失效应变明显高于高强钢和铝合金，若把低碳钢的ε_f直接赋给DP780，会导致B柱在碰撞中不断裂、侵入量被低估。反过来，把铝合金的失效曲线用于低碳钢吸能盒，则会在小变形阶段就删除单元，吸能总量严重不足。材料建模的最终目标不是“参数好看”，而是让碰撞力—位移曲线、断裂位置和吸能分配与试验或对标结果一致。