2. 基于失效物理的应力-强度干涉模型建立
问题定义与建模对象
某机械连接件（M10 8.8级螺栓，承受轴向拉-拉循环载荷）在服役中可能出现两种失效模式：螺栓杆部屈服断裂或螺纹脱扣。本例以杆部疲劳断裂为失效模式，建立应力-强度干涉模型。
已知条件（教学示例数据）：
- 工作载荷 F ~ N(μ_F=28 kN, σ_F=4.2 kN)，服从正态分布
- 螺栓杆部截面积 A = 58 mm²（M10 细牙，d₂≈8.6 mm）
- 材料强度 S ~ N(μ_S=620 MPa, σ_S=45 MPa)，由同批次材料试验统计得到
- 应力 σ = F/A，为载荷的线性函数
应力分布推导
应力是载荷的线性变换，正态分布经线性变换后仍为正态分布：
μ_σ = μ_F / A = 28000 N / 58 mm² = 482.8 MPa
σ_σ = σ_F / A = 4200 N / 58 mm² = 72.4 MPa
σ ~ N(482.8, 72.4²) MPa强度分布 S ~ N(620, 45²) MPa。两者均值差 μ_S − μ_σ = 137.2 MPa，但标准差合成后 σ_combined = √(45² + 72.4²) = 85.2 MPa，干涉区不可忽略。
失效概率计算
应力与强度均服从正态分布且相互独立时，干涉概率（失效概率）为：
z = (μ_S − μ_σ) / √(σ_S² + σ_σ²)
= (620 − 482.8) / √(45² + 72.4²)
= 137.2 / 85.2
= 1.61
P_f = 1 − Φ(1.61) = 1 − 0.9463 = 0.0537即在给定载荷分布下，该连接件的失效概率约为 5.37%，对应可靠度 R ≈ 94.63%。
关键设计变量识别
通过灵敏度分析，识别对失效概率影响最大的设计变量：
| 设计变量 | 当前值 | 变化方向 | 对 P_f 的影响 | 灵敏度等级 |
|---|---|---|---|---|
| 杆部截面积 A | 58 mm² | 增大 | μ_σ 下降，P_f 显著降低 | 高 |
| 材料强度均值 μ_S | 620 MPa | 提高 | P_f 降低 | 高 |
| 载荷均值 μ_F | 28 kN | 降低 | P_f 降低 | 中 |
| 强度标准差 σ_S | 45 MPa | 减小 | P_f 降低 | 中 |
| 载荷标准差 σ_F | 4.2 kN | 减小 | P_f 降低 | 低 |
设计改进方向：将截面积从 58 mm² 增至 65 mm²（如改用 M12 螺栓），μ_σ 降至 430.8 MPa，z 值升至 2.19，P_f 降至约 1.43%。若同时将材料强度均值提升至 650 MPa，P_f 可进一步降至 0.5% 以下。
验证与判读
- 模型验证：用蒙特卡洛模拟 10⁶ 次抽样，统计 S < σ 的频率，应与解析解 5.37% 一致（允许 ±0.3% 抽样误差）。
- 判读准则：若目标可靠度 R_target ≥ 99%，则当前设计不满足，需调整截面积或材料等级。
- 注意事项：应力-强度干涉模型假设 S 与 σ 独立。若强度随载荷循环退化（疲劳），需引入强度退化模型，本例仅考虑单次加载下的静态干涉。