2. 材料非线性弹塑性拉伸模拟
问题定义与材料模型
本案例对一根标准金属拉伸试样进行弹塑性分析。试样标距段长度 L₀ = 50 mm，截面为矩形 12 mm × 3 mm，一端固定，另一端施加轴向位移。材料采用双线性等向强化模型，需要输入弹性模量、泊松比、屈服应力和塑性硬化模量四个关键参数。
双线性等向强化在 Abaqus 中通过 *PLASTIC 关键字定义，每行给出（屈服应力，塑性应变）数据对。第一行塑性应变为 0，对应初始屈服；第二行对应硬化后的应力-塑性应变点。Abaqus 会自动在两点之间线性插值，形成双线性硬化段。
材料参数（教学示例值，模拟 Q235 类低碳钢）：
- 弹性模量 E = 210000 MPa
- 泊松比 ν = 0.30
- 初始屈服应力 σ_y = 235 MPa（对应塑性应变 0）
- 硬化点：σ = 420 MPa（对应塑性应变 0.10）
- 由两点可推得塑性硬化模量 H' = (420−235)/0.10 = 1850 MPa
需要特别注意：Abaqus 的 *PLASTIC 数据中第二列是塑性应变，不是总应变。若误填总应变，硬化斜率会被严重高估。正确换算关系为 ε_plastic = ε_total − σ/E。
建模关键设置
在 Abaqus/Standard 中建立二维平面应力模型即可高效完成本案例。关键设置如下：
- 分析步：Static, General，打开
NLGEOM=ON（虽然本例应变水平不高，但保持与后续大变形案例一致的习惯）。初始增量步 0.05，最小 1e-5，最大 0.1。 - 单元类型：CPE4R（四节点双线性平面应力四边形，减缩积分），网格尺寸约 1 mm，标距段至少 50 个单元。
- 边界条件：左端面完全固定（U1=U2=0），右端面耦合到参考点 RP-1，在 RP-1 上施加 U1 = 5 mm 的位移。
- 输出请求：在 History Output 中请求 RP-1 的
RF1（反力）和U1（位移），用于绘制载荷-位移曲线。
网格收敛性提示：弹塑性问题中，塑性区对应力集中敏感。建议先用 1 mm 网格试算，再将标距段加密到 0.5 mm 对比 RF1 峰值差异。若差异小于 2%，可认为网格收敛。
载荷-位移曲线提取与实验对比
提交作业后，在 Visualization 模块中创建 XY Data：选择 RP-1 的 RF1 作为 Y 轴、U1 作为 X 轴，即可得到模拟载荷-位移曲线。将实验曲线以 XY Data 形式导入同一图表进行对比。
对比时重点关注三个特征量：
| 特征量 | 模拟值（示例） | 实验值（示例） | 偏差判断 |
|---|---|---|---|
| 弹性段斜率 | 由 E 和几何直接决定 | 由拉伸机刚度修正后得到 | 偏差 < 5% 为正常 |
| 屈服载荷 | F_y = σ_y × A = 235 × 36 = 8460 N | 从曲线拐点读取 | 偏差 < 8% 可接受 |
| 塑性段斜率 | 由 H' 和几何决定 | 由实验曲线拟合 | 偏差 < 15% 可接受 |
常见错误：若模拟曲线在屈服后立即出现水平平台（无硬化），说明 *PLASTIC 只填了一行数据，Abaqus 将其视为理想塑性。若塑性段斜率远大于预期，检查第二列是否误填了总应变。
结果判读与验证
完成对比后，按以下标准判断模型是否可用：
- 弹性段斜率与实验值偏差在 5% 以内——说明 E 和截面尺寸正确。
- 屈服点位置与实验拐点偏差在 8% 以内——说明 σ_y 合理。
- 塑性段趋势与实验一致，且最大载荷偏差在 10% 以内——说明 H' 取值合理。
- 若三项均满足，可认为双线性等向强化模型在本案例的应变范围内（ε ≤ 0.10）具有足够精度。
若塑性段偏差较大，可考虑改用多段 *PLASTIC 数据或 *CYCLIC 硬化模型，但这已超出本例范围。