4、疲劳寿命预测方法:名义应力法(S-N曲线法)、局部应力应变法(Neuber法)、断裂力学法(Paris公式)

各位工程师朋友,大家好。今天我们来聊聊疲劳寿命预测的三大主流方法。说实话,我在钛合金项目上摸爬滚打这么多年,这三种方法几乎每个项目都会用到。它们各有各的脾气,也各有各的适用场景。你想想看,如果选错了方法,那预测结果可能差出好几个数量级,这可不是闹着玩的。

我个人习惯把这三种方法比作「三把尺子」:一把量宏观,一把量局部,一把量裂纹。咱们一个一个来看。

4.1 名义应力法(S-N曲线法)

名义应力法,说白了就是最经典、最直观的方法。它基于材料的S-N曲线(应力-寿命曲线),直接通过名义应力幅来预测寿命。

核心思想: 构件在循环载荷下,只要名义应力幅不超过某个阈值,就能承受无限次循环。超过这个阈值,寿命就由S-N曲线决定。

适用场景:

  • 高周疲劳(N > 10⁴ ~ 10⁵ 次)
  • 结构简单、应力集中不严重的情况
  • 初步设计阶段的快速评估

具体步骤:

  1. 确定载荷谱,计算名义应力幅 Sₐ
  2. 查材料的S-N曲线(或通过试验获得)
  3. 考虑各种修正系数(尺寸系数、表面系数、应力集中系数Kt等)
  4. 用修正后的应力幅查寿命N

S-N曲线的数学表达(Basquin公式):

Sₐ = σ'f · (2N)^b

其中:

  • Sₐ —— 应力幅
  • σ'f —— 疲劳强度系数(≈ 真实断裂强度)
  • b —— 疲劳强度指数(通常在 -0.05 ~ -0.12 之间)
  • N —— 循环次数

重要提示: 钛合金的S-N曲线通常没有明显的水平段(即没有真正的疲劳极限)。这意味着,理论上只要应力不为零,总会发生疲劳破坏。只不过在低应力下,寿命会非常长。

我的经验: 我在做某型飞机起落架钛合金部件时,一开始直接用光滑试样的S-N曲线,结果预测寿命比实际短了3倍。后来发现,表面加工状态的影响太大了。钛合金对表面缺陷特别敏感,抛光与车削的寿命能差一个数量级。所以,我建议一定要用与实际构件表面状态一致的S-N曲线。

4.2 局部应力应变法(Neuber法)

名义应力法有个硬伤——它不考虑局部塑性变形。你想想看,在缺口根部,应力早就超过屈服点了,这时候再用名义应力法,误差就大了。

局部应力应变法就是来解决这个问题的。它关注的是缺口根部的真实应力应变状态。

Neuber法的核心公式:

Kt² = Kσ · Kε

其中:

  • Kt —— 理论应力集中系数(弹性)
  • Kσ —— 真实应力集中系数(σ_local / σ_nominal)
  • Kε —— 真实应变集中系数(ε_local / ε_nominal)

具体步骤:

  1. 计算名义应力 S 和理论应力集中系数 Kt
  2. 用Neuber公式迭代求解局部应力 σ 和局部应变 ε
  3. 结合材料的循环应力-应变曲线(σ-ε 滞回环)
  4. 用Manson-Coffin公式计算裂纹形成寿命

Manson-Coffin公式:

εₐ = εe + εp = (σ'f / E) · (2N)^b + ε'f · (2N)^c

其中:

  • εₐ —— 总应变幅
  • εe —— 弹性应变幅
  • εp —— 塑性应变幅
  • ε'f —— 疲劳延性系数
  • c —— 疲劳延性指数(通常在 -0.5 ~ -0.7 之间)

注意: Neuber法适用于低周疲劳(N < 10⁴ ~ 10⁵ 次),也就是塑性应变占主导的情况。对于高周疲劳,塑性应变很小,Neuber法的优势就不明显了。

我曾经踩过的坑: 有一次做钛合金叶片的缺口分析,我直接用Neuber法算局部应力,结果发现算出来的应力比实际屈服强度还高出一大截。后来才意识到,Neuber法在严重塑性变形时会有偏差。我现在的做法是:先用Neuber法估算,再用有限元验证一下,两者取平均。嗯,这样稳妥多了。

4.3 断裂力学法(Paris公式)

前面两种方法,预测的是「裂纹形成寿命」。但实际构件中,往往存在初始缺陷或微裂纹。这时候,断裂力学法就派上用场了——它预测的是「裂纹扩展寿命」。

Paris公式:

da/dN = C · (ΔK)^m

其中:

  • da/dN —— 裂纹扩展速率(每循环扩展的长度)
  • ΔK —— 应力强度因子幅(Kmax - Kmin)
  • C, m —— 材料常数(由试验确定)

具体步骤:

  1. 确定初始裂纹尺寸 a₀(通常取无损检测的检出极限,如 0.1~0.5 mm)
  2. 计算临界裂纹尺寸 a_c(由断裂韧性 K_IC 决定)
  3. 对Paris公式积分,得到从 a₀ 扩展到 a_c 所需的循环次数

积分公式:

N = ∫_{a₀}^{a_c} da / [C · (ΔK)^m]

关键点: 钛合金的Paris指数 m 通常在 3~4 之间。m值越大,裂纹扩展越快,对初始缺陷越敏感。我见过有些钛合金的m值高达4.5,那真是「一裂就断」,设计时必须留足安全裕度。

三种方法的对比:

方法 适用阶段 输入参数 输出结果 典型应用
名义应力法 裂纹形成 名义应力幅、S-N曲线 总寿命 高周疲劳、初步设计
局部应力应变法 裂纹形成 局部应力应变、循环σ-ε曲线 裂纹形成寿命 低周疲劳、缺口分析
断裂力学法 裂纹扩展 初始裂纹尺寸、Paris参数 裂纹扩展寿命 损伤容限设计、检修周期

我的建议: 在实际工程中,我通常把三种方法结合起来用。比如,用名义应力法做初步筛选,用局部应力应变法评估关键缺口,再用断裂力学法确定检修周期。这样「三管齐下」,心里才踏实。

知识体系框架图

下面这张图,是我自己总结的三种方法的逻辑关系。你看一眼就能明白它们各自的位置和衔接点。

疲劳寿命预测方法体系 名义应力法 局部应力应变法 断裂力学法 核心参数 S-N曲线、应力幅 核心参数 Neuber公式、滞回环 核心参数 Paris公式、ΔK 适用:高周疲劳 初步设计阶段 适用:低周疲劳 缺口/局部塑性区 适用:裂纹扩展 损伤容限设计 工程实践:三法结合,互为补充 名义应力法 → 局部应力应变法 → 断裂力学法

好了,这一章的内容就到这里。三种方法各有千秋,关键是要根据实际情况灵活选用。记住,没有万能的方法,只有最合适的组合。


公众号:蓝海资料掘金营,微信deep3321