7.3 疲劳、断裂与工程建模

疲劳指构件在循环应力远低于强度极限时，经数十万次循环后萌生裂纹并扩展至断裂的现象，设计依据是 S-N 曲线与损伤累积法则。断裂力学则假定裂纹已经存在，用应力强度因子裁决它是否扩展。本节完成疲劳寿命估算、裂纹临界尺寸计算，并以工程建模的完整清单为全册收尾。

结案流程"测、算、判"的最后一步是判。它与第 5 章的接口是应力计算（循环应力的幅值与均值）、与第 6 章的接口是动载荷（应力循环的来源）、与 7.2 节的接口是应力集中的数值评估。

一、疲劳：低于强度极限的慢性病

静强度合格的结构在服役几年后断裂，断口找不到明显塑性变形——这是疲劳的签名。机理分三幕：微裂纹在应力集中处萌生（滑移带挤出挤入）、裂纹每次循环推进一小步（扩展区留下海滩纹）、剩余截面不足以承载时瞬间断掉（瞬断区）。应力幅值越低，能撑的循环次数越多，S-N 曲线记录这层关系；钢在循环一千万次后存在水平段（疲劳极限），低于它理论寿命无限。


import numpy as np

# 某合金钢 S-N 曲线简化：Sa = Sf' * (2N)^b，取 Sf'=900 MPa, b=-0.09
Sf, b = 900.0, -0.09
for Sa in [500, 400, 300, 250]:
    N = 0.5 * (Sa/Sf)**(1/b)
    print(f"应力幅 {Sa} MPa：寿命约 {N:.1e} 次循环")
# 若疲劳极限 240 MPa
print("应力幅压到 240 MPa 以下（疲劳极限），寿命趋于无限")

应力幅从 500 降到 300，寿命从几十万次涨到几十亿次——疲劳对应力水平极度敏感，这也是为什么应力集中系数控制、表面抛光、喷丸强化（引入残余压应力）这些"细节工艺"能决定疲劳寿命的量级。

变幅载荷用 Miner 法则累积损伤：每种应力水平的循环数除以各自寿命，加起来等于 1 时失效判定。


import numpy as np

# 一天载荷谱：三种应力幅各循环若干次
spectrum = [(350, 200), (300, 2000), (250, 20000)]   # (MPa, 次/天)
Sf, b = 900.0, -0.09
D_daily = 0.0
for Sa, n in spectrum:
    N = 0.5 * (Sa/Sf)**(1/b)
    D_daily += n / N
    print(f"应力幅 {Sa} MPa：寿命 {N:.2e} 次，日损伤 {n/N:.3e}")
life_days = 1 / D_daily
life_years = life_days / 365
print(f"日累积损伤 {D_daily:.4f}，估算寿命 {life_days:.0f} 天 = {life_years:.1f} 年")
# 敏感性：最大应力幅升 10%
spectrum2 = [(385, 200), (300, 2000), (250, 20000)]
D2 = sum(n / (0.5*(Sa/Sf)**(1/b)) for Sa, n in spectrum2)
print(f"最大档应力升 10%：寿命降到 {1/D2/365:.1f} 年（降 {(1-D2/ (1/life_days*1) if False else D_daily/D2):.0%} 量级）")

只把最大一档应力抬高 10%，总寿命就明显缩水——Miner 法则展示的"短板效应"是结构健康监测只盯高应力事件的原因。

二、断裂力学：带着裂纹生活

疲劳裂纹萌生后，问题换成断裂力学的语言：裂纹前端附近的应力场强度由应力强度因子 K 描述（等于名义应力乘几何修正因子乘根号派乘裂纹长度）。K 达到材料的断裂韧度时裂纹失稳扩展；低于门槛值时基本不扩展；两者之间按 Paris 公式缓慢扩展。无损检测查出裂纹后，工程关心的正是"还剩多少寿命"。


import numpy as np
from scipy.integrate import quad

# 拉板含中心裂纹：K = sigma * sqrt(pi*a)，断裂韧度 Kc=80 MPa√m
Kc, sigma, Y = 80e6, 120e6, 1.12    # 边缘裂纹修正系数
a_cr = (Kc / (Y*sigma))**2 / np.pi
print(f"名义应力 120 MPa 下临界裂纹尺寸 {a_cr*1000:.1f} mm")
# Paris 扩展律：da/dN = C*(dK)^m，C=6.9e-12（MPa√m, m/循环），m=3
C, m = 6.9e-12, 3.0
dK = lambda a: Y * sigma * np.sqrt(np.pi * a) / 1e6   # MPa√m
dNda = lambda a: 1 / (C * dK(a)**m)
N, _ = quad(dNda, 0.5e-3, a_cr)     # 从 0.5 mm 扩到临界
print(f"从 0.5 mm 扩展至断裂需 {N:.3e} 次循环")
# 检测周期：剩余寿命的 1/2 作为检查间隔
print(f"按剩余寿命一半设检测周期：{N/2:.2e} 次循环一查")

从半毫米裂纹到断裂若需要两百万次循环，检测周期定一百万次——"损伤容限"设计就是把断裂力学与无损检测编成一张时间表，承认裂纹存在、管理裂纹扩展。

图疲劳三阶段与断裂判据

三、工程建模：全册的总复盘

把 6 章的武器按流程串起来，任何构件的分析都可以归纳为一张清单：


import numpy as np

# 综合演练：起重机吊臂连接销轴的完整校核流水线（伪代码式清单演示）
# 输入：拉力 150 kN，销轴 d=40 mm 双剪，材料调质钢
F, d = 150e3, 40e-3
A_shear = 2 * np.pi*d**2/4          # 双剪两个剪切面
tau = F / A_shear
print(f"第1步 静强度：切应力 {tau/1e6:.1f} MPa < 许用 190 MPa -> 通过")
# 弯曲（视为短梁近似略）——工程判断：销轴以剪切与挤压为主
sig_b = F / (d * 25e-3) / 1e6       # 挤压（孔径投影）
print(f"第2步 挤压：{sig_b:.1f} MPa < 许用 300 MPa -> 通过")
# 动载系数与疲劳
Kd = 1.3                             # 起升冲击
tau_max = tau * Kd
print(f"第3步 动载放大 {Kd}：切应力峰值 {tau_max/1e6:.1f} MPa")
Sa = (tau_max - tau)/1e6 / 2        # 循环幅值近似
print(f"第4步 疲劳：应力幅约 {Sa:.1f} MPa，对照切应力 S-N 许用值评估寿命")
print("第5步 结论：静强度挤压通过，疲劳寿命给出更换周期建议")

五步走完，静强度、挤压、动载、疲劳各判一次——这份"多失效模式逐一排查"的纪律，比任何一个单独公式都更接近工程力学的本质。

⚠️ 常见坑：把静强度校核当成全部。疲劳、稳定（第 5 章）、共振（第 6 章）都不在静强度视野内；一份只算了静强度的报告，等于只排查了失效模式清单的第一行。

本节要点回顾

疲劳是低应力慢性病：断口三幕（源、海滩纹、瞬断区）是诊断指纹，寿命对应力幅极度敏感。
Miner 损伤线性累积：各档循环数除以寿命求和，高应力事件主导总损伤。
断裂力学承认裂纹存在：应力强度因子对比断裂韧度定临界尺寸，Paris 律算扩展，检测周期取剩余寿命一半。
两条设计哲学并存：无限寿命（压应力幅到疲劳极限以下）与损伤容限（管理裂纹扩展时间表）。
多失效模式排查清单：静强度、挤压、稳定、共振、疲劳逐一过堂，缺一项就留一类风险在门外。

全册至此结案：从第 1 章的第一张自由体图，到本节的寿命时间表，载荷路径侦探完成了从"力在哪"到"还能用多久"的完整推理链。

作者与出处

原作者：灏天文库

来源：灏天文库

整理：灏天文库整理

由灏天文库平台收录，内容或由平台用户上传，仅供学习交流

7.3 疲劳、断裂与工程建模