4. 材料非线性基础
金属结构在服役中一旦局部应力超过屈服点，刚度会明显下降，卸载后还会留下残余变形。Abaqus 要正确复现这一行为，必须把材料的弹性、屈服、硬化与塑性流动四件事同时描述清楚。本节围绕这四个概念展开，并给出在 Abaqus 中把单轴拉伸曲线转成塑性参数、再用同一试验做回代验证的完整做法。
一、弹塑性理论：把总应变拆成两部分
小变形弹塑性理论的核心假设是应变加法分解：总应变 ε 等于弹性应变 εᵉ 与塑性应变 εᵖ 之和，即 ε = εᵉ + εᵖ。弹性部分服从胡克定律 σ = E·εᵉ，卸载时按斜率 E 返回；塑性部分不可恢复，且不产生应力增量。判断是否进入塑性的门槛是屈服函数 f(σ, 硬化变量)。当 f < 0 时材料处于弹性状态；f = 0 时材料位于屈服面上，继续加载会产生塑性流动；f > 0 在率无关塑性中不允许出现，求解器会通过返回映射把应力拉回屈服面。这一“弹性预测—塑性修正”的径向返回算法，正是 Abaqus/Standard 每次迭代内部实际执行的计算。
二、屈服准则：金属常用 Mises
金属材料一般对静水压力不敏感，因此 Abaqus 中金属塑性默认采用 Mises 屈服准则。它用等效应力（von Mises 应力）q 与当前屈服应力 σ_y 比较：
q = sqrt( 0.5 * [ (σ1-σ2)^2 + (σ2-σ3)^2 + (σ3-σ1)^2 ] )
屈服条件： q = σ_y( εᵖ )在单轴拉伸中，σ1 = σ，σ2 = σ3 = 0，此时 q 恰好等于轴向应力 σ，所以单轴试验曲线可以直接作为 Mises 屈服面的标定依据。这也是为什么金属塑性参数几乎都从单轴拉伸获得。若材料存在明显拉压不对称（如铸铁、部分高强铝合金），则需改用 Drucker-Prager 或 Hill 等各向异性/压力相关准则，本节不展开。
三、硬化规则：等向硬化与随动硬化
屈服面在塑性变形中会变化，描述这种变化的方式就是硬化规则。Abaqus 中金属塑性最常用的是等向硬化：屈服面中心不动，半径随塑性应变增大而扩大，即 σ_y 随 εᵖ 单调上升。它适合单调加载、比例加载问题，实现简单、参数少。另一类是随动硬化（Abaqus 中为 Combined/kinematic 硬化），屈服面半径基本不变而中心平移，用于描述包辛格效应和循环塑性。若结构经历反复加载、卸载或地震、疲劳工况，仅用等向硬化会高估反向屈服强度，应改用随动或混合硬化。选择哪种规则，取决于载荷路径是否反向，而不是材料牌号本身。
四、流动法则：塑性应变往哪走
流动法则决定塑性应变增量的方向。Abaqus 金属塑性采用关联流动法则，即塑性应变增量垂直于屈服面：dεᵖ = dλ · ∂f/∂σ。对 Mises 屈服面，∂f/∂σ 的方向与偏应力方向一致，因此塑性变形体积不变（不可压），这一结论直接导致两个工程后果：一是塑性区泊松比趋于 0.5，二是必须使用杂交单元（如 C3D8H）或缩减积分单元来避免体积自锁。理解这一点，比记住公式更重要。
五、Abaqus 中金属塑性参数标定实操
假设某低碳钢单轴拉伸试验给出：弹性模量 E = 210000 MPa，泊松比 ν = 0.3，屈服点 σ_y0 = 235 MPa，并在塑性段测得若干数据点。Abaqus 的 Plastic 表要求输入“屈服应力—塑性应变”，因此必须先把总应变换算成塑性应变：εᵖ = ε_total − σ/E。以表中一点 σ = 300 MPa、ε_total = 0.012 为例，弹性应变 = 300/210000 ≈ 0.00143，则 εᵖ ≈ 0.01057。首行必须填 εᵖ = 0、σ = 235，表示初始屈服。若把总应变直接填入，会人为抬高屈服应力，导致结构偏刚。
在 Abaqus/CAE 中的操作路径为：Property 模块 → Material → Mechanical → Plasticity → Plastic，逐行输入上述数据；硬化类型保持 Isotropic。若需随动硬化，则改用 Combined 并补充 Kinematic 参数。材料定义完成后赋给截面，再建立与试验一致的拉伸试样模型：一端固支，另一端施加位移，输出 RF 与 U。
六、单轴拉伸验证与结果判读
验证的关键是“用标定参数回算同一试验”。提取加载端反力 RF 与位移 U，换算成工程应力—应变曲线，与试验曲线叠放比较。判读要点有三：第一，弹性段斜率是否与 E 一致，若不一致说明弹性参数或截面面积有误；第二，初始屈服点是否落在 235 MPa 附近，偏差大通常源于屈服点取值或单位制混乱；第三，塑性段是否贴合，若模拟偏硬，多为 εᵖ 计算错误或漏减弹性应变。一般工程验证要求屈服点误差小于 3%，塑性段整体误差小于 5%。若偏差集中在反向加载或循环工况，则应回到硬化规则，改用随动硬化重新标定。