4. 材料循环本构行为

在单轴疲劳分析中，我们习惯用一条单调应力—应变曲线描述材料刚度与强度。但进入循环载荷后，材料会“记住”自己经历过的塑性变形历史，其响应不再由单调曲线唯一决定。多轴疲劳寿命分析之所以复杂，根本原因就在于：损伤累积发生在循环本构路径上，而不是发生在单调拉伸曲线上。本节围绕循环应力应变曲线、滞回环、循环软化/硬化、记忆效应与多轴循环塑性五个相互关联的概念展开，帮助建立“循环载荷→塑性响应→损伤驱动量”的完整认识。

一、循环应力应变曲线：稳定循环的骨架

对试件施加恒应变幅循环，例如应变幅 Δε/2 = 0.5%，控制应变比 R = −1，在室温下循环若干周次。前几十周应力幅会明显变化，随后逐渐趋于稳定。把多个不同应变幅水平下稳定后的应力幅与应变幅对应起来，就得到循环应力应变曲线。它通常用 Ramberg–Osgood 形式拟合：

Δε/2 = Δσ/(2E) + (Δσ/(2K'))^(1/n')

其中：
Δε/2  应变幅
Δσ/2  应力幅
E     弹性模量
K'    循环强度系数
n'    循环应变硬化指数

以某低合金钢为例（示例参数，非实测）：E = 210 GPa，K' = 1200 MPa，n' = 0.15。当应变幅为 0.5% 时，可反解出应力幅约为 430 MPa 量级。循环曲线与单调曲线往往不重合：若循环曲线位于单调曲线下方，说明材料循环软化；位于上方则循环硬化。工程上判断“用哪条曲线做弹塑性分析”，第一步就是比较这两条曲线。

二、滞回环：一次循环内的能量与路径

把稳定循环中一个完整加卸载周期的应力—应变轨迹画出来，就得到滞回环。它的宽度反映塑性应变幅，面积反映每循环耗散的塑性功。滞回环不是简单闭合曲线，其形状包含丰富信息：

• 环的宽度 Δεp 越大，塑性损伤越显著；
• 环的尖端对应拉、压峰值应力，可判断平均应力与应力比；
• 卸载段斜率接近弹性模量，若明显偏斜，说明存在显著非比例或黏塑性效应；
• 环面积 Wp = ∮σ dε，是能量型疲劳判据的重要输入。

在多轴加载下，滞回环不再是一条平面曲线，而是应力、应变多分量随时间的闭合轨迹。此时常用等效塑性应变幅与等效应力幅来降维描述，但必须注意：等效化会丢失相位信息，对非比例加载可能低估损伤。

三、循环软化与循环硬化：稳定前的过渡

材料在循环初期并不稳定。若在恒定应变幅下应力幅逐周下降，称为循环软化；若逐周上升，称为循环硬化。常见规律是：冷加工硬化材料（如冷轧钢）倾向于循环软化，退火态软材料（如退火铜）倾向于循环硬化。稳定周次可能从几十周到几千周不等，取决于材料、应变幅与温度。

实际分析中，不能直接用第一周或单调曲线代替稳定响应。一个可操作的判断流程是：先做 3～5 个应变幅水平的恒幅循环试验，记录应力幅随周次变化；若在 100～500 周内变化小于 2%，可认为进入稳定；取稳定段拟合 K' 与 n'。若寿命分析关注裂纹萌生早期，则需保留软化/硬化的瞬态段，采用 Chaboche 等随动强化模型逐周积分。

四、记忆效应：材料如何“记住”历史

记忆效应指材料当前屈服行为依赖于此前的加载路径。典型表现是：先拉伸进入塑性再卸载，反向加载时屈服应力低于单调压缩屈服应力，这就是包辛格效应；若中间插入一次反向加载，后续再加载的屈服面位置会随之移动。在随动强化框架下，用背应力张量 α 描述屈服面中心偏移，其演化方程使材料“记住”塑性应变路径。

多轴非比例加载下，记忆效应更明显。例如先拉—扭比例加载，再切换为 90° 相位差非比例加载，等效塑性应变幅可能不变，但应力幅会额外升高，称为非比例附加强化。若忽略这一效应，寿命预测会偏保守或偏危险，取决于所采用的损伤模型。

五、多轴循环塑性：从单轴到多轴的扩展

多轴循环塑性的核心是屈服准则与强化法则。常用 von Mises 屈服面配合随动强化：

f = √(3/2 · (s − α):(s − α)) − σy ≤ 0

s     偏应力张量
α     背应力张量
σy    当前屈服应力

背应力演化可采用 Armstrong–Frederick 形式：dα = (2/3)C dεp − γ α dp，其中 C、γ 控制强化速率与安定性。多轴加载时，塑性流动方向由屈服面法向决定，非比例路径会导致塑性应变增量方向不断旋转，从而产生附加强化与额外的能量耗散。