2. Abaqus材料非线性分析环境与流程
在Abaqus中做金属弹塑性仿真，第一步不是急着画网格，而是先弄清楚求解器、材料模块、单位制和分析步这四件事如何协同。Abaqus/Standard与Abaqus/Explicit虽然共享同一套材料库，但对非线性问题的处理路径完全不同：前者用隐式Newton-Raphson迭代求解静力平衡，适合准静态拉伸、装配过盈等缓变过程；后者用显式中心差分推进，适合冲击、碰撞、成形等高速瞬态问题。选择哪一个，直接决定了材料参数、分析步类型和输出请求的写法。
求解器与材料模块的对应关系
Abaqus的材料定义集中在Property模块，但不同求解器对同一材料模型的支持程度不同。以最常用的金属塑性为例，Standard支持*PLASTIC（Mises屈服+各向同性硬化）、*KINEMATIC HARDENING（随动硬化，用于循环加载）以及*HYPERELASTIC等；Explicit同样支持这些模型，但更强调密度和稳定时间增量。对于本节要完成的简单金属拉伸试样，采用Standard求解器、Mises各向同性硬化即可，材料参数只需弹性模量、泊松比、屈服应力和塑性应变-应力数据对。
需要特别注意的是，Abaqus没有内置单位制，所有输入必须自洽。工程中最常用的是mm-N-MPa-s体系：长度mm、力N、应力MPa、密度t/mm³。例如钢的密度写作7.85e-9 t/mm³，弹性模量210000 MPa，屈服强度235 MPa。如果误用kg-m-Pa体系，密度会变成7850 kg/m³，弹性模量写成2.1e11 Pa，虽然物理上等价，但一旦混用就会导致质量矩阵和刚度矩阵量级错乱，显式分析甚至无法启动。
分析步设置的关键参数
在Standard中做弹塑性拉伸，分析步类型选Static, General，并打开NLGEOM=ON以考虑几何非线性。增量步控制是收敛的关键：初始增量不宜过大，建议取总时间的5%左右；最小增量步可设为1e-5，防止因局部塑性应变集中导致迭代失败；最大增量步可设为0.1，避免跳过屈服拐点。如果模型在屈服附近出现不收敛，通常不是材料参数错误，而是增量步太大导致Newton迭代在弹塑性交界处反复振荡，此时把初始增量降到0.01往往能直接解决。
边界条件方面，左端面施加ENCASTRE全约束，右端面施加U1=5 mm的位移载荷。这里有一个容易忽略的细节：如果右端只约束U1而放开U2、U3，试样在拉伸过程中会因泊松效应产生横向收缩，这是正确的；但如果右端同时约束了U2和U3，就会人为抑制颈缩，导致载荷-位移曲线偏高。对于简单拉伸试样，右端应只保留U1方向的位移加载，其余方向自由。
输出请求与结果提取
场输出（Field Output）默认每增量步写入，建议至少勾选U（位移）、RF（反力）、S（应力）、PEEQ（等效塑性应变）。其中PEEQ是判断材料是否进入塑性的直接指标：PEEQ=0表示弹性，PEEQ>0表示已屈服。历史输出（History Output）用于提取右端参考点的RF1和U1，在Visualization模块中可以直接绘制载荷-位移曲线。如果曲线在屈服后出现明显下降，说明试样发生了颈缩或局部失稳，这属于正常物理现象，不是计算错误。
一个可复现的示例参数集
以下参数可直接用于建立本节拉伸试样模型，材料为Q235钢，单位制mm-N-MPa-s：
| 项目 | 参数值 | 说明 |
|---|---|---|
| 弹性模量 E | 210000 MPa | 线弹性段斜率 |
| 泊松比 ν | 0.3 | 横向收缩比 |
| 屈服应力 σ_y | 235 MPa | 塑性起始点 |
| 塑性数据对 | (0,235), (0.02,310), (0.05,360), (0.10,400) | ε_p–σ 表格 |
| 密度 ρ | 7.85e-9 t/mm³ | Standard静力分析可忽略，Explicit必需 |
| 分析步总时间 | 1.0 | 无量纲，对应5 mm位移 |
| 初始增量 | 0.05 | 屈服前可适当放大 |
提交计算后，在Monitor中观察增量步收敛情况。如果某增量步出现“Attempt 1–5”反复迭代，说明该步长下塑性流动导致刚度突变，可手动将该步长减半。计算完成后，在Visualization中提取右端RF1–U1曲线，与上表塑性数据对的理论趋势对比：弹性段斜率应接近210000 MPa×截面积20 mm²=4200 N/mm，屈服点对应RF1≈235×20=4700 N。若仿真曲线在4700 N附近出现明显拐点，说明材料定义和分析步设置正确。