2. 材料非线性弹塑性拉伸
问题定义与材料参数
单轴拉伸试样采用 100 mm × 10 mm × 10 mm 的实体单元模型，一端固定，另一端施加轴向位移。材料为双线性等向强化弹塑性模型，弹性模量 E = 210 GPa，泊松比 ν = 0.3，初始屈服应力 σy0 = 250 MPa，塑性硬化模量 H = 2 GPa（即屈服后切线模量 Et = E·H/(E+H) ≈ 1.98 GPa）。
分析分两个分析步：Step-1 施加轴向拉伸位移至 2.0 mm（对应名义应变 2%），Step-2 将位移卸载回 0。通过提取加载末态与卸载后的应变分量，可分离弹性应变、塑性应变与残余应变。
双线性等向强化中，屈服面随塑性应变等向扩张，背应力为零。卸载按弹性模量 E 进行，因此残余应变 = 总应变 − 弹性应变，塑性应变在卸载过程中保持不变。
建模关键设置
- 分析类型：Static, General；开启 NLGEOM 以考虑大变形（本例应变 2%，几何非线性影响较小，但保持一致性）。
- 单元类型：C3D8R（八节点线性减缩积分），沿轴向至少划分 20 个单元以捕捉应变分布。
- 边界条件：左端面 ENCASTRE（U1=U2=U3=0），右端面耦合至参考点 RP-1，在 RP-1 上施加 U1 位移。
- 材料定义：Elastic（E=210000 MPa, ν=0.3）+ Plastic（Yield Stress=250, Plastic Strain=0；Yield Stress=254, Plastic Strain=0.002）。
- 场输出：在 Step-1 与 Step-2 中均请求 E（总应变）、PE（塑性应变）、PEEQ（等效塑性应变）、S（应力）、U（位移）。
双线性等向强化应力—应变路径
下图展示单轴拉伸加载—卸载过程中，应力—应变路径在双线性等向强化模型下的演化。加载段沿弹性斜率 E 上升至屈服点 σy0，随后沿硬化斜率 Et 继续上升；卸载段沿弹性斜率 E 返回，最终在零应力处留下残余应变 εres，而塑性应变 εp 在卸载前后保持不变。
关键结果提取与验证
在 Abaqus/CAE 后处理中，通过以下方式提取并验证结果：
- 创建 Field Output 的
E（总应变）、PE（塑性应变）、PEEQ（等效塑性应变）云图，确认加载末态塑性区分布。 - 沿试样轴向定义 Path，绘制
E11与PE11沿路径的分布曲线，检查均匀段应变一致性。 - 提取卸载后同一路径的
E11与PE11，计算残余应变 εres = E11(卸载后) − E11(卸载前) + PE11(卸载前)。 - 验证：卸载后
PEEQ应与加载末态相同（塑性应变不可恢复），E11应等于PE11（弹性应变归零）。
示例数据（教学假设）：加载至名义应变 2% 时，总应变 εtotal = 0.02，屈服应变 εy = 250/210000 ≈ 0.00119，塑性应变 εp = 0.02 − 0.00119 − (σmax−σy0)/E ≈ 0.0188。卸载后弹性应变恢复，残余应变 εres ≈ 0.0188，与塑性应变一致。
常见问题与排查
| 现象 | 可能原因 | 排查方法 |
|---|---|---|
| 卸载后 PEEQ 发生变化 | 材料定义中误用了随动强化或混合强化 | 检查 Plastic 卡片是否仅输入等向强化数据，确认无 Kinematic 关键字 |
| 残余应变与塑性应变不相等 | 卸载未完全回到零应力，或输出变量选择错误 | 确认 Step-2 位移边界已归零，检查 E 与 PE 的输出请求 |
| 屈服点不明显 | 单元积分点应力平滑或网格过粗 | 细化轴向网格，提取积分点应力—应变曲线 |
| 分析不收敛 | 位移步长过大或 NLGEOM 设置不当 | 减小初始增量步，检查材料切线模量是否合理 |