3. 局部应力应变法估算缺口件裂纹萌生寿命
问题定位：缺口根部的真实应力应变
缺口件在名义应力远低于屈服强度时，缺口根部已可能进入塑性。局部应力应变法的核心思路是：用Neuber修正把名义弹性应力/应变换算为缺口根部的真实弹塑性应力/应变，再沿循环应力应变曲线确定稳定滞回环，最后用Manson-Coffin方程把应变幅换算为裂纹萌生寿命。
本节以一块中心圆孔板为例：材料为某低合金钢，弹性模量 E=206000 MPa，循环强度系数 K'=1200 MPa，循环应变硬化指数 n'=0.15，疲劳强度系数 σ'f=950 MPa，疲劳延性系数 ε'f=0.45，疲劳强度指数 b=-0.09，疲劳延性指数 c=-0.55。承受名义应力幅 ΔS/2=180 MPa 的恒幅循环载荷，应力比 R=-1，缺口弹性应力集中系数 Kt=3.0。
Neuber修正：从名义量到局部量
Neuber规则给出缺口根部局部应力幅 Δσ 与局部应变幅 Δε 的乘积约束：
Δσ·Δε = (Kt·ΔS)² / E
其中 ΔS 为名义应力幅。同时局部点必须落在循环应力应变曲线上：
Δε/2 = Δσ/(2E) + (Δσ/(2K'))1/n'
两式联立即可解出 Δσ 和 Δε。下面用迭代法求解：先假设纯弹性 Δσ0=Kt·ΔS=540 MPa，代入Neuber式得 Δε0=(KtΔS)²/(E·Δσ0)=0.00262，再代入循环曲线校核并修正，迭代收敛后得到局部应力幅约 Δσ≈412 MPa，局部应变幅约 Δε≈0.00485（即 εa≈0.00243）。
迭代计算要点：Neuber双曲线 Δσ·Δε=C 与循环应力应变曲线 Δε=f(Δσ) 的交点即为局部响应。若交点处 Δσ/2 < σy'（循环屈服强度），说明缺口根部仍为弹性，可直接用 Kt 放大；否则必须用上述弹塑性解。
寿命估算与结果判读
将 εa=0.00243 代入 Manson-Coffin 方程，弹性项贡献约 0.00100，塑性项贡献约 0.00143。塑性项占比约 59%，说明缺口根部已明显进入塑性，寿命处于低周疲劳区。求解方程得到 2Nf≈1.2×105，即裂纹萌生寿命 Nf≈6.0×104 次循环。
| 参数 | 数值 | 说明 |
|---|---|---|
| 名义应力幅 ΔS/2 | 180 MPa | 远低于屈服，但缺口处已塑性 |
| 弹性应力集中系数 Kt | 3.0 | 纯弹性峰值 540 MPa |
| Neuber常数 C | 1.415 mm·MPa | (KtΔS)²/E |
| 局部应力幅 Δσ | ≈412 MPa | Neuber与循环曲线交点 |
| 局部应变幅 εa | ≈0.00243 | Δε/2 |
| 弹性应变分量 | ≈0.00100 | σ'f/E·(2Nf)b |
| 塑性应变分量 | ≈0.00143 | ε'f·(2Nf)c |
| 裂纹萌生寿命 Nf | ≈6.0×104 | 2Nf≈1.2×105 |
工程判断：若局部应力幅 Δσ/2 小于材料循环屈服强度，Neuber解退化为弹性解，可直接用 Kt 放大名义应力；若塑性项占比超过 50%，说明寿命对塑性参数 ε'f 和 c 更敏感，材料数据应优先保证这两个参数的准确性。
注意：Neuber修正假设缺口根部处于平面应力状态且无残余应力。若缺口件经过喷丸、滚压等表面强化处理，或存在焊接残余应力，需在局部应力中叠加残余应力后再代入Manson-Coffin方程，否则会高估萌生寿命。