1. 课程导入与问题定义

在工程结构设计中,疲劳失效是最常见也最危险的失效模式之一。桥梁缆索在风振与车流反复作用下的裂纹萌生、航空发动机涡轮盘在离心载荷与热循环下的低周疲劳、风电叶片在随机风载下的寿命衰减——这些场景都有一个共同特征:载荷是随机的、材料性能是分散的、几何与工艺存在偏差,而寿命预测结果往往呈现数量级级别的离散。如果只给出一个"平均寿命"数值,工程决策将失去可靠依据。这正是疲劳寿命不确定性评估要解决的核心问题。

一、为什么必须评估不确定性

传统疲劳寿命预测通常采用确定性方法:给定应力幅 S,通过 S-N 曲线或断裂力学公式计算寿命 N。但实际工程中存在三类不可忽略的随机性来源:

  • 载荷随机性:服役载荷谱并非恒定,应力幅、均值、循环次数均服从某种分布。例如某汽车悬架实测载荷谱中,应力幅在 80~220 MPa 之间波动,循环次数每天约 1.2×10⁴ 次。
  • 材料与工艺分散性:同批次材料的疲劳极限可能相差 15%~30%,缺陷尺寸、表面粗糙度、热处理残余应力都会引入偏差。
  • 模型不确定性:S-N 曲线拟合误差、Miner 线性累积损伤假设的偏差、裂纹扩展 Paris 公式参数 C、m 的统计波动。

当这些随机变量同时作用时,寿命预测结果不再是一个点值,而是一个分布。工程上真正关心的是:在给定可靠度下,寿命下限是多少?失效概率是否低于可接受阈值?

核心目标:建立从"随机输入"到"寿命分布"再到"可靠度指标"的完整量化链路,使设计、运维与检测决策有概率依据,而非依赖单一安全系数。

二、基本术语与符号约定

在进入后续章节前,需要统一以下术语与参照约定。本课程中,应力以拉为正、压为负;应力幅 Sa 定义为循环中最大应力与最小应力差值的一半;应力比 R = Smin/Smax。寿命 N 统一指失效时的循环次数,单位为次。

术语符号含义典型示例
应力幅Sa循环应力半幅值120 MPa
疲劳寿命N失效前循环次数2.5×10⁵ 次
随机变量X输入不确定性量材料强度、载荷幅值
失效概率Pf寿命低于阈值的概率1×10⁻⁴
可靠度指标β安全裕度标准化度量β = 3.72
变异系数COV标准差/均值0.18

其中可靠度指标 β 与失效概率的关系为 Pf = Φ(−β),Φ 为标准正态分布函数。例如 β = 3.72 对应 Pf ≈ 1×10⁻⁴,这是许多工程结构可接受的风险水平。

三、整体学习路径

本课程共 30 章,逻辑上分为五个阶段,形成一条从"认识问题"到"工程落地"的完整链路:

  1. 基础与定义(第1~5章):建立不确定性评估的工程背景、概率统计基础、随机变量建模方法。
  2. 寿命预测模型(第6~12章):S-N 曲线、Miner 损伤累积、断裂力学、局部应力应变法及其不确定性传播。
  3. 不确定性量化方法(第13~20章):蒙特卡洛模拟、一次二阶矩法、响应面法、代理模型与灵敏度分析。
  4. 可靠度与寿命分布(第21~26章):可靠度计算、置信区间、贝叶斯更新、小样本数据处理。
  5. 工程应用与决策(第27~30章):检测周期优化、维修决策、标准规范对照与综合案例。

下面用一张结构图展示本课程的知识拓扑与数据流关系,帮助建立全局认知:

疲劳寿命不确定性评估 · 知识拓扑与数据流 阶段一 · 基础定义 工程背景 / 概率统计 随机变量建模 阶段二 · 寿命模型 S-N / Miner / 断裂力学 局部应力应变法 阶段三 · 不确定性量化 蒙特卡洛 / FORM / 响应面 代理模型 / 灵敏度 阶段四 · 可靠度 可靠度计算 / 置信区间 贝叶斯更新 / 小样本 阶段五 · 工程应用与决策(第27~30章) 检测周期优化 · 维修决策 · 标准规范对照 · 综合案例 随机输入 载荷谱 / 材料参数 / 几何 寿命预测模型 损伤累积 / 裂纹扩展 不确定性传播 抽样 / 代理 / 灵敏度 寿命分布 P_f / β / 置信区间 数据更新与模型修正(贝叶斯反馈) 核心问题:随机输入 → 寿命分布 → 可靠度指标 → 工程决策 本课程目标:掌握从建模、量化到决策的完整方法链,使寿命评估结果具有概率意义

四、一个直观的示例:确定性 vs 不确定性

假设某构件在恒幅载荷 Sa = 150 MPa 下工作,材料 S-N 曲线中位寿命为 N50 = 3.2×10⁵ 次。若只按中位值设计,意味着约 50% 的构件会提前失效——这显然不可接受。若考虑材料分散性(寿命服从对数正态分布,对数标准差 σlnN = 0.22)和载荷波动(Sa 的 COV = 0.12),通过蒙特卡洛模拟 10⁴ 次可得:

示例参数(合理假设,非实测):
  S_a ~ Normal(150, 18²) MPa
  ln(N) ~ Normal(ln(3.2e5), 0.22²)
  Miner 损伤阈值 D_cr = 1.0

模拟结果(示意):
  中位寿命 N_50 ≈ 2.9×10⁵ 次
  5% 分位寿命 N_5 ≈ 1.4×10⁵ 次
  失效概率 P_f(N < 1.5×10⁵) ≈ 6.8×10⁻²
  可靠度指标 β ≈ 1.49

可以看到,考虑不确定性后,5% 分位寿命仅为中位寿命的约 48%。如果工程要求 Pf ≤ 1×10⁻⁴(β ≥ 3.72),则当前设计远不满足,必须降低应力幅、提高材料性能或缩短检测周期。这就是不确定性评估的实际价值:它把"安全"从一个模糊概念转化为可计算、可比较、可决策的概率指标。

学习建议:本课程强调"可计算、可复现"。每章尽量给出明确的输入参数、计算过程与结果判读方法。遇到公式时,先理解其物理含义与适用条件,再关注数值实现。

五、本课程的方法论定位

需要明确的是,不确定性评估不是要否定确定性方法,而是在其基础上增加概率维度。确定性计算给出"寿命是多少",不确定性评估回答"寿命低于某值的可能性有多大"。两者结合,才能支撑完整的工程决策。后续章节将逐步展开:先建立概率统计与随机变量基础,再深入疲劳寿命预测模型,然后系统学习不确定性量化算法,最后落到可靠度计算与工程决策。整条路径的核心逻辑始终是:识别随机源 → 建立概率模型 → 传播不确定性 → 量化可靠度 → 支持决策。

在进入第2章之前,建议先回顾概率论中的随机变量、分布函数、期望与方差等基本概念,并思考一个自己熟悉的疲劳问题:它的随机源有哪些?哪些可以量化?哪些暂时只能靠工程判断?带着这些问题学习,后续方法的适用场景会更加清晰。