3. 应变-寿命方法:e-N曲线、Coffin-Manson公式与循环应力应变响应

各位,咱们今天聊点实在的。上一章讲了应力-寿命(S-N)方法,说白了就是看应力多大能扛多久。但实际工作中,尤其是那些有缺口的零件、焊缝、或者高温部件,局部早就进入塑性了。这时候再用应力说话,就不太灵了。

我个人习惯,遇到这种局部大应变的情况,直接上应变-寿命(e-N)方法。这个方法的核心思想很简单:疲劳寿命由局部应变决定。你想想看,裂纹总是在应力集中的地方萌生,那个地方的金属早就屈服了,应变才是真正的“元凶”。

3.1 应变-寿命(e-N)曲线

e-N曲线长什么样?跟S-N曲线有点像,但横坐标还是寿命N(循环次数),纵坐标换成了总应变幅Δε/2(或者应变范围Δε)。

这条曲线其实是由两条线叠加而成的:

  • 弹性应变线:高周疲劳区,寿命长,应变小,基本是弹性变形。这条线在双对数坐标下是条直线。
  • 塑性应变线:低周疲劳区,寿命短,应变大,有明显的塑性变形。这条线也是直线,但斜率更陡。

总应变幅就是这两条线的和。嗯,这里要注意:低周疲劳看塑性,高周疲劳看弹性。我在做涡轮盘低循环疲劳分析时,就发现盘心孔边那点地方,每次启动停机都要经历一次大应变循环,塑性应变占主导,用S-N方法算出来寿命偏危险。

核心概念:e-N曲线是总应变幅与失效循环次数的关系曲线,由弹性分量和塑性分量叠加而成。低周疲劳(LCF)以塑性应变为主,高周疲劳(HCF)以弹性应变为主。

3.2 Coffin-Manson公式

说到e-N曲线,就绕不开Coffin-Manson公式。这是整个应变-寿命方法的基石。公式长这样:

Δε/2 = (σ'f / E) * (2Nf)^b + ε'f * (2Nf)^c

别被这一串符号吓到,我来拆开讲:

  • Δε/2:总应变幅,就是你要算的输入量。
  • σ'f:疲劳强度系数,可以理解为材料在单次循环就能拉断的应力水平。单位MPa。
  • E:弹性模量,这个大家都熟。
  • b:疲劳强度指数,弹性线的斜率,一般在-0.05到-0.12之间。
  • ε'f:疲劳延性系数,可以理解为材料在单次循环就能拉断的应变能力。
  • c:疲劳延性指数,塑性线的斜率,一般在-0.5到-0.7之间。
  • 2Nf:反向次数,也就是失效时的循环次数乘以2(一个循环包含两次反向)。

公式左边是总应变幅,右边第一项是弹性应变幅,第二项是塑性应变幅。说白了,就是弹性项和塑性项加起来等于总应变。

我的经验:拿到一个新材料的Coffin-Manson参数,我第一件事就是画两条线——弹性线和塑性线。看它们的交点在哪里。交点对应的寿命,就是材料从低周疲劳过渡到高周疲劳的分界点。这个点很有用,能帮你快速判断工况属于哪个区。

举个例子,某航空铝合金的典型参数:

参数 符号 典型值 单位
疲劳强度系数 σ'f 1100 MPa
疲劳强度指数 b -0.08 -
疲劳延性系数 ε'f 0.35 -
疲劳延性指数 c -0.60 -
弹性模量 E 71000 MPa

代入公式,如果总应变幅是0.01,算出来寿命大概在几千次到几万次之间。如果总应变幅降到0.001,寿命就能到百万次以上。这就是低周和高周的区别。

3.3 循环应力应变响应

好,现在有了e-N曲线和Coffin-Manson公式,但还有一个关键问题没解决:我怎么知道局部应变是多少?

这就引出了循环应力应变响应。材料在循环加载下,应力应变关系不是一成不变的。你想想看,第一次加载和第一百次加载,材料的“脾气”能一样吗?

主要有三种现象:

  • 循环硬化:随着循环次数增加,材料变“硬”了,同样的应变需要更大的应力。我遇到过一些不锈钢材料,刚加载时软绵绵的,循环几百次后明显变硬。
  • 循环软化:跟上面相反,材料变“软”了。高强度钢经常有这个问题。我曾经吃过这个亏——用静载的应力应变曲线去算一个高强钢接头的局部应变,结果寿命预测偏乐观,后来发现材料在循环中软化了。
  • 循环稳定:有些材料比较“老实”,循环几次后就稳定了,应力应变关系不再变化。退火态的铝合金通常是这样。

描述循环应力应变响应,最常用的是Ramberg-Osgood公式

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

其中:

  • ε:总应变
  • σ:应力
  • K':循环强度系数
  • n':循环应变硬化指数

这个公式跟单调拉伸的Ramberg-Osgood公式长得一样,但参数K'和n'是循环加载下测出来的,跟静载参数完全不同。千万别混用!

避坑指南:我曾经犯过一个低级错误——直接用单调拉伸的K和n值去算循环应力应变响应。结果算出来的滞回环形状完全不对,寿命预测偏差了一个数量级。记住:循环参数和单调参数是两码事,必须通过循环试验获得。

有了循环应力应变曲线,再加上Neuber法则应变梯度法,就能从名义应力/应变推算出缺口根部的局部应变。然后代入Coffin-Manson公式,就能算出寿命。

整个流程我画了张图,方便大家理解:

应变-寿命(e-N)方法分析流程 ① 输入条件 几何形状、载荷谱 材料参数(E, K', n'等) ② 名义应力/应变 有限元分析或解析计算 得到名义应力应变历程 ③ 局部应变计算 Neuber法则 / 应变梯度法 结合循环σ-ε曲线 ④ 循环计数 雨流计数法 提取应变循环幅值 ⑤ 损伤计算 Coffin-Manson公式 计算每个循环的损伤 ⑥ 累积损伤 Miner线性累积法则 总损伤 = Σ(Ni/Nfi) 局部应变历程 应变幅值 ⑦ 寿命预测结果 总寿命 Nf = 1 / 总损伤 总损伤值

这张图把整个流程串起来了。从输入条件开始,到最终寿命预测结束。中间每一步都有对应的理论和方法支撑。

我个人觉得,这个流程里最容易出问题的就是第三步——局部应变计算。名义应力好算,但缺口根部的局部应变,受材料循环特性影响很大。如果循环应力应变曲线没选对,后面全白搭。

总结一下本章要点

  • e-N曲线由弹性线和塑性线叠加而成,低周疲劳看塑性,高周疲劳看弹性。
  • Coffin-Manson公式是定量计算的基础,四个参数(σ'f, b, ε'f, c)必须通过试验获得。
  • 循环应力应变响应描述材料在循环加载下的“真实”行为,循环硬化/软化现象不可忽视。
  • Ramberg-Osgood公式配合Neuber法则是工程上最常用的局部应变估算方法。

好了,这一章就到这里。记住,e-N方法不是万能的,但在低周疲劳和局部应变主导的场合,它比S-N方法靠谱得多。下一章咱们聊聊多轴疲劳,那又是另一番天地了。


专注资料整理