4. 应变疲劳与局部应力应变法
当构件在缺口根部、焊缝趾部或应力集中区进入塑性后，名义应力法（S-N曲线）就失效了——因为S-N曲线描述的是光滑试样的名义应力与寿命关系，而缺口根部的真实应力已经超过屈服，局部应变成为控制损伤的主导参量。应变疲劳（低周疲劳）正是处理这一问题的核心方法：以局部应变幅为控制变量，用Manson-Coffin方程描述应变-寿命关系，再通过Neuber修正把名义载荷换算到缺口根部的局部应力应变，最终估算缺口件寿命。
一、循环应力应变响应：材料在反复加载下的真实行为
单调拉伸曲线只描述一次加载。反复加载时，材料会表现出循环硬化、循环软化或循环稳定三种响应。判断依据是循环应力幅随周次的变化：若在恒定应变幅控制下，应力幅逐周升高，为循环硬化；逐周降低，为循环软化；基本不变，则为循环稳定。例如退火态低碳钢通常循环硬化，冷加工态钢通常循环软化，而许多调质钢在中等应变幅下接近循环稳定。
稳定后的循环应力应变曲线可用Ramberg-Osgood关系描述：
ε_a = σ_a/E + (σ_a/K')^(1/n')
其中：
ε_a —— 总应变幅
σ_a —— 应力幅
E —— 弹性模量
K' —— 循环强度系数
n' —— 循环应变硬化指数该式把总应变幅分解为弹性部分 ε_ae = σ_a/E 和塑性部分 ε_ap = (σ_a/K')^(1/n')。当 ε_ap 可忽略时，材料处于弹性主导；当 ε_ap 与 ε_ae 相当或更大时，塑性损伤不可忽略，必须用应变疲劳方法。工程上常用“过渡寿命”判断：总应变幅中塑性分量占比超过约5%时，建议转入应变疲劳分析。
二、Manson-Coffin方程：应变幅与寿命的定量关系
Manson-Coffin方程把总应变幅写成弹性项与塑性项之和：
Δε/2 = (σ'_f/E)(2N_f)^b + ε'_f(2N_f)^c
其中：
Δε/2 —— 总应变幅
σ'_f —— 疲劳强度系数
b —— 疲劳强度指数（通常 -0.05 ~ -0.12）
ε'_f —— 疲劳延性系数
c —— 疲劳延性指数（通常 -0.5 ~ -0.7）
2N_f —— 失效反向数（N_f 为循环周次）
E —— 弹性模量弹性项 (σ'_f/E)(2N_f)^b 在双对数坐标上是一条斜率为 b 的直线，塑性项 ε'_f(2N_f)^c 是斜率为 c 的直线。两条直线的交点对应过渡寿命 2N_t：低于该寿命时塑性项主导，高于该寿命时弹性项主导。以某低合金钢为例（示例参数，非实测）：σ'_f=1200 MPa，b=-0.08，ε'_f=0.6，c=-0.6，E=210 GPa。当 Δε/2=0.01 时，弹性项约 0.0035，塑性项约 0.0065，塑性占比约65%，说明该应变水平下塑性损伤显著，必须用应变法而非名义应力法。
三、Neuber修正：从名义应力到缺口根部局部应力应变
Neuber修正解决的是“已知名义应力幅 S_a 和应力集中系数 K_t，如何求缺口根部局部应力幅 σ_a 和局部应变幅 ε_a”的问题。其核心关系为：
K_t² · S_a² / E = σ_a · ε_a
即：σ_a · ε_a = (K_t · S_a)² / E该式右侧是已知量，左侧是待求的局部应力与局部应变的乘积。要解出 σ_a 和 ε_a，还需联立循环应力应变曲线：
ε_a = σ_a/E + (σ_a/K')^(1/n')两式联立后，通常用迭代法求解：先假设 σ_a，代入Ramberg-Osgood式算出 ε_a，再检查 σ_a·ε_a 是否等于 (K_t·S_a)²/E，若不相等则调整 σ_a 重算，直到收敛。工程上也可用图解法：在循环应力应变曲线上找一点，使其与原点构成的矩形面积等于 (K_t·S_a)²/E。
以示例说明：某缺口件 K_t=3，名义应力幅 S_a=200 MPa，E=210 GPa，材料 K'=1200 MPa，n'=0.15。则 (K_t·S_a)²/E = (600)²/210000 = 1.714 MPa。假设 σ_a=450 MPa，代入 ε_a=450/210000+(450/1200)^(1/0.15)=0.00214+0.375^6.67≈0.00214+0.0015=0.00364，乘积 σ_a·ε_a=1.638 MPa，略小于1.714，需增大 σ_a。再取 σ_a=460 MPa，ε_a≈0.00219+0.0018=0.00399，乘积=1.835 MPa，略大。取 σ_a≈455 MPa，ε_a≈0.0038，乘积≈1.729 MPa，接近目标。因此局部应力幅约455 MPa，局部应变幅约0.0038。该结果说明缺口根部已进入塑性（455 MPa超过屈服），必须用应变疲劳方法。
四、缺口件寿命估算：完整闭环
得到局部应变幅 Δε/2=0.0038 后，代入Manson-Coffin方程求解 2N_f。仍用前述示例参数：σ'_f=1200 MPa，b=-0.08，ε'_f=0.6，c=-0.6，E=210 GPa。方程变为：
0.0038 = (1200/210000)(2N_f)^(-0.08) + 0.6(2N_f)^(-0.6)
= 0.00571(2N_f)^(-0.08) + 0.6(2N_f)^(-0.6)试算：令 2N_f=10000，弹性项=0.00571×10000^(-0.08)=0.00571×0.479=0.00274；塑性项=0.6×10000^(-0.6)=0.6×0.00398=0.00239；总和=0.00513，大于0.0038。令 2N_f=100000，弹性项=0.00571×0.398=0.00227；塑性项=0.6×0.000631=0.000379；总和=0.00265，小于0.0038。说明寿命在10000~100000反向之间。再试 2N_f=30000：弹性项=0.00571×0.438=0.00250；塑性项=0.6×0.00199=0.00119；总和=0.00369，接近0.0038。取 2N_f≈28000，则 N_f≈14000 周次。该结果为示例估算值，实际工程中需用材料实测参数并考虑表面加工、尺寸效应和平均应力修正。
实际工程中，缺口件寿命估算的完整链条是：名义载荷谱 → 名义应力幅 → Neuber修正 → 局部应力应变 → Manson-Coffin方程 → 寿命 N_f。每一步的输入参数都应来自材料实测或可靠手册，示例中的数值仅用于说明计算过程。对于焊接接头、螺栓孔、键槽等典型缺口，还应叠加残余应力、表面粗糙度和尺寸系数的修正，才能得到可用于工程判断的寿命估算值。