第三章 常用仿真材料模型：从线弹性到Johnson-Cook

做插拔力仿真这些年，我最大的体会就是：材料模型选对了，仿真就成功了一半。选错了呢？嗯，结果可能连趋势都对不上。

今天咱们就聊聊几种最常用的材料模型。我会结合自己的项目经验，把每个模型的适用场景、参数含义、校准方法都讲透。

3.1 线弹性模型：最基础的起点

线弹性模型，说白了就是胡克定律的翻版。应力跟应变成正比，比例系数就是弹性模量 E。再加一个泊松比 ν 描述横向变形，齐活。

什么时候用这个模型？我个人习惯是：初步分析、刚度评估、或者材料确实只发生小变形时。比如插拔力仿真中，如果只是评估插针的弹性变形阶段，线弹性完全够用。

但要注意——一旦材料进入塑性，线弹性模型就失效了。我见过不少新手拿线弹性模型算插拔力，结果应力都超过屈服强度了还在那儿算，那数据能对吗？

3.2 双线性弹塑性模型：工程中最常用的选择

双线性弹塑性模型，就是把应力-应变曲线简化为两段直线：弹性段和塑性段。弹性段斜率是 E，塑性段斜率是切线模量 Et。

你想想看，实际材料的应力-应变曲线哪有这么直？但为什么这个模型还这么流行？因为简单、参数少、够用。

我建议，对于大多数金属材料的插拔力仿真，双线性模型是首选。为什么？因为插拔过程中，材料通常只进入轻度塑性，双线性模型能很好地捕捉这个行为。

3.3 多线性弹塑性模型：追求更真实的描述

双线性模型虽然好用，但有些材料的应力-应变曲线拐弯抹角的，两段直线实在拟合不好。这时候就该多线性模型上场了。

多线性模型，就是用多段直线来逼近真实的应力-应变曲线。段数越多，拟合精度越高，但参数也越多。

// ABAQUS 中定义多线性塑性（示例）
*PLASTIC
350, 0.0      // 屈服点：应力350MPa，塑性应变为0
420, 0.02     // 第二点：应力420MPa，塑性应变0.02
480, 0.05     // 第三点：应力480MPa，塑性应变0.05
520, 0.10     // 第四点：应力520MPa，塑性应变0.10
550, 0.15     // 第五点：应力550MPa，塑性应变0.15

这里要注意：输入的应变是塑性应变，不是总应变。总应变 = 弹性应变 + 塑性应变。这个坑我踩过，当时怎么算应力都不对，后来才发现是把总应变当塑性应变输进去了。

3.4 Johnson-Cook模型：考虑应变率和温度

Johnson-Cook模型，简称J-C模型。它跟前面几个模型最大的区别是：考虑了应变率效应和温度软化效应。

公式长这样：

σ = (A + B·εp^n) · (1 + C·ln(ε̇/ε̇0)) · (1 - T*^m)

其中：
A —— 屈服应力（参考应变率下）
B —— 应变硬化系数
n —— 应变硬化指数
C —— 应变率敏感系数
m —— 温度软化指数
T* —— 无量纲温度

说实话，插拔力仿真中用到J-C模型的情况不多。因为插拔速度通常不快，应变率效应不明显。但有一种情况例外——高速插拔或者冲击式插拔。

我记得有个客户做汽车连接器，要求仿真插拔速度达到 100mm/s。普通模型算出来的力偏小，换了J-C模型后，结果跟实测对上了。为什么？因为高速下材料会表现出明显的应变率强化效应。

3.5 如何选择？我的决策流程

讲了这么多，到底该用哪个模型？我给大家梳理一个简单的决策流程：

1. 先问自己：材料会进入塑性吗？ 不会 → 线弹性模型
2. 会进入塑性，但变形不大？ → 双线性弹塑性模型
3. 需要精确描述硬化行为？ → 多线性弹塑性模型
4. 高速加载或温度变化大？ → Johnson-Cook模型

这个流程我用了好多年，基本没出过问题。当然，具体问题还得具体分析。