随机疲劳的 “概率核心”：四大损伤模型与概率密度函数解析
1. 概述
随机响应疲劳分析聚焦于随机载荷下结构的疲劳寿命研究。其核心逻辑为：通过随机响应分析获得功率谱密度（PSD）结果，基于该结果计算PSD力矩，进而生成周期数与应力范围的概率密度函数，最终结合材料SN曲线完成疲劳损伤与寿命评估。该方法适用于路面颠簸、机械振动、气流扰动等无固定规律载荷场景，广泛应用于汽车、航空航天、工程机械等领域。
2. 功率谱密度（PSD）力矩
2.1 计算原理
PSD力矩是量化随机应力能量频率分布的核心参数，基于随机响应分析输出的应力PSD结果计算，是后续概率密度函数生成的基础。
2.2 计算公式
2.3 参数说明
2.4 示意图
图 1. PSD 力矩计算示意图
3. 应力范围出现概率计算
3.1 核心定义
针对每个应力范围区间 ，计算该区间内应力范围的出现概率，概率值由用户定义的参数或损伤模型推导。
3.2 计算公式
3.3 参数说明
4. 概率密度函数（周期数与应力范围的概率密度）
4.1 生成逻辑
PSD力矩用于构建应力范围概率密度函数 ，函数形式由指定的损伤模型决定。工程中支持四种损伤模型：DIRLIK、LALANNE、NARROW、Steinberg三波段（THREE）。
4.2 各损伤模型详细说明
4.2.1 DIRLIK模型
适用于一般宽带随机载荷（如多频率叠加的振动），通过三组分分布叠加精准捕捉宽范围应力分布特性，是通用性最强的模型。
（1）概率密度函数公式
（2）衍生参数定义
• 标准化变量： • 不规则因子（载荷带宽表征）： • 统计量： • 权重系数：
• 形状参数：
（3）参数说明
4.2.2 LALANNE模型
适用于宽带随机载荷，公式简化且计算效率高，通过指数分布+误差函数修正适配宽带特性，兼顾精度与速度。
（1）概率密度函数公式
（2）核心参数
• 不规则因子： （与DIRLIK模型定义一致） • ：误差函数，表征正态分布累积概率特性 • ：应力范围（单位：MPa） • ：0阶、2阶、4阶PSD力矩
4.2.3 NARROW模型（窄带模型）
适用于理想窄带随机载荷（如单一频率主导的旋转机械振动），基于瑞利分布，物理意义明确且计算简洁。
（1）概率密度函数公式
（2）关键说明
• 核心参数： 为应力范围（MPa）， 为0阶PSD力矩 • 循环次数计算：默认使用零交叉数 ，而非峰值数 。原因是 的数值计算易受噪声干扰，稳定性较差；理想窄带下零交叉数与峰值数近似相等，不影响精度。
4.2.4 THREE模型（Steinberg三波段模型）
适用于概念设计阶段的快速估算，为离散型简化模型，直接给出三个特征应力范围的概率，无需复杂积分。
（1）概率公式
$$P(S) = \begin{cases}
0.683 & (S = 2\sqrt{m_0}) \
0.271 & (S = 4\sqrt{m_0}) \
0.043 & (S = 6\sqrt{m_0})
\end{cases}$$
（2）核心特性
• 与其他模型的区别： 为概率（而非概率密度），三个概率值之和为1（0.683+0.271+0.043=1.0） • 应力范围依据： 分别对应1倍、2倍、3倍RMS应力（ ），概率取值源于正态分布3σ原则 • 周期数计算：直接通过 总零交叉数得到，无需乘以区间宽度
4.3 示意图
图 2. 概率密度函数（周期数与应力范围的概率密度）示意图
5. 应力范围相关参数设置
概率密度函数的计算精度依赖于应力范围的参数定义，支持四种设置方式，用户可根据工程需求选择：
5.1 上曲面应力范围因子
（1）计算公式
（2）说明
• RMS应力由随机响应子工况直接输出（ ） • 应力范围上限用于限制计算边界，超出该上限的应力在疲劳损伤计算中不予考虑 • 上曲面应力范围因子为工程经验值（通常取3~5），需根据载荷特性调整
5.2 直接指定上曲面应力范围
根据试验数据或工程经验，直接输入应力范围上限值（如200MPa），适用于有明确应力极限要求的场景。
5.3 循环次数
（1）应力范围宽度计算公式
（2）默认规则
• 循环次数默认值为100 • 第一个应力范围区间为 ，后续区间依次递增（如 ），直至覆盖至应力范围上限
5.4 直接定义应力范围宽度
直接输入应力范围区间宽度 （如2MPa），区间宽度越小，概率计算精度越高，但计算量相应增加。
6. 应力范围出现概率的细化计算
6.1 DIRLIK、LALANNE、NARROW模型
• 适用场景：连续型概率密度函数 • 说明： 为第 个区间应力范围 对应的概率密度值， 为区间宽度
6.2 THREE模型
直接采用模型定义的三个概率值（0.683、0.271、0.043），无需额外计算：
• 每个特征应力范围的周期数 = 对应概率 × 总零交叉数 • 仅包含三个应力范围区间，分别为
7. 损伤模型选择指南
7.1 选择依据
损伤模型的核心是匹配随机载荷的应力范围分布特性，即应力范围是集中在特定值附近（窄带）还是覆盖较广范围（宽带）。
7.2 各模型适用场景
7.3 补充说明
• DIRLIK和LALANNE模型能捕捉应力范围分布的广谱特性，适合大多数工程场景 • NARROW模型仅适用于载荷能量集中在单一频率的场景（如电机、泵的稳定振动） • THREE模型基于简化假设，不建议用于关键结构的精确评估
8. 峰值数和零交叉数计算
峰值数和零交叉数是量化随机载荷循环频率的核心参数，用于计算总周期数：
8.1 NARROW、THREE模型（零交叉数）
• 定义：单位时间内应力信号从负变正穿过均值线的次数 • 适用场景：窄带载荷或简化模型，稳定性高于峰值数
8.2 DIRLIK、LALANNE模型（峰值数）
• 定义：单位时间内应力信号出现峰值的次数 • 适用场景：宽带载荷，能更好反映高频成分的循环特性 • 注意事项： 计算易受噪声干扰，需确保PSD数据质量
8.3 参数说明
分别为0阶、2阶、4阶PSD力矩
9. 周期数计算
9.1 总周期数
（1）NARROW、THREE模型
（2）DIRLIK、LALANNE模型
（3）参数说明
9.2 特定应力范围的总周期数
• 定义：第 个应力范围 对应的实际循环次数 • 说明： 为该应力范围的出现概率， 为总周期数
10. 疲劳寿命与损伤计算
10.1 疲劳寿命计算（单应力范围）
根据材料SN曲线，计算单一应力范围 下的疲劳寿命：
参数说明
10.2 总疲劳损伤计算（多应力范围）
采用Miner线性累积损伤理论，叠加各应力范围的损伤贡献：
结果判据
• ：结构未失效，剩余寿命充足 • ：结构恰好达到失效临界状态 • ：结构已失效，存在安全隐患
10.3 平均应力修正
若载荷导致平均应力不为零（如重力、预紧力叠加随机振动），需进行平均应力修正：
1. 定义静态子工况：将产生平均应力的载荷（如重力）单独定义为静态子工况 2. 关联分析事件：在随机疲劳分析事件设置中引用该静态子工况，软件将自动结合静态应力与随机应力，完成平均应力修正
11. 总结
随机响应疲劳分析的核心技术链路为：随机载荷→PSD响应→PSD力矩→概率密度函数→周期数统计→疲劳损伤/寿命。关键要点包括：
1. 选择适配载荷特性的损伤模型（DIRLIK为通用首选） 2. 合理设置应力范围相关参数（上限、区间宽度等），平衡精度与效率 3. 确保PSD数据质量，避免高阶力矩（如 ）计算不稳定 4. 结合材料SN曲线与平均应力修正，提升评估结果的准确性
该方法为随机载荷下结构的疲劳寿命预测提供了科学、系统的解决方案，是工程结构可靠性设计与安全评估的核心技术手段。
