1. 随机振动疲劳问题导入
在轨道车辆、航空发动机、车载电子和风电装备等领域，大量结构件并非在单一频率的正弦载荷下工作，而是长期承受宽带随机振动。这类载荷的幅值和频率成分随时间不断变化，用传统的定频疲劳分析方法往往无法解释实际出现的早期开裂。随机振动疲劳寿命仿真要解决的，正是"在给定随机载荷谱下，结构能安全服役多久"这一核心问题。
一个典型的工程失效场景
某车载电子控制盒通过支架安装在车身地板上。整车道路试验中，控制盒在累计行驶约 4 万公里后，支架根部出现贯穿裂纹，而按早期正弦扫频试验评估的寿命远大于此。事后分析发现：正弦扫频只考核了单一共振点的定频响应，而实际路面激励是宽带随机谱，支架在 80～400 Hz 范围内同时被多个模态激发，应力响应呈现随机波动，其峰值分布和循环计数方式与定频完全不同。这个案例说明，随机振动疲劳的关键难点不在"算应力"，而在于如何把随机应力时程转化为可用于疲劳损伤累积的循环统计量。
由此可以明确本节要建立的整体认识：随机振动疲劳寿命仿真的目标，是在给定的功率谱密度（PSD）激励下，预测结构危险点的疲劳寿命或损伤，并判断是否满足设计寿命要求。它的输入通常包括激励 PSD 谱、结构有限元模型、材料 S-N 曲线、阻尼与边界条件；输出则是危险点位置、应力响应 PSD、损伤值 D 与预测寿命。
从激励到寿命的分析链路
随机振动疲劳仿真并不是一步完成的，它是一条从载荷到损伤的传递链。理解这条链上每一环的输入输出，是避免"跑完仿真却不知道结果是否可信"的前提。下图给出该链路的核心结构与关键参数传递关系。
关键参数与判断依据
链路上最容易被忽视的是谱矩与不规则因子。应力响应 PSD 的零阶谱矩 m0 对应应力均方值，二阶谱矩 m2 与峰值率相关，四阶谱矩 m4 用于判断带宽。不规则因子 γ = m2 / √(m0·m4)，当 γ 接近 1 时响应接近窄带，Dirlik 等频域方法精度较高；当 γ 明显小于 1 时，响应呈宽带特征，需谨慎选择损伤模型或改用时域雨流计数。以车载控制盒为例，若仿真得到危险点 m0 = (12 MPa)²、峰值率 νp ≈ 95 Hz、γ ≈ 0.72，则说明响应并非理想窄带，直接套用窄带近似会高估或低估损伤，必须用宽带修正方法。
另一个关键判断是阻尼比。支架类结构常用 ζ = 0.02～0.05，若误取 0.01，共振峰放大倍数近似翻倍，损伤可能被高估数倍。因此仿真报告中应明确写出阻尼来源——是模态试验测得，还是经验取值。
分析流程与结果判读
完整的随机振动疲劳仿真通常按以下顺序推进：先建立并验证有限元模型，通过模态分析确认前若干阶频率与试验或经验值吻合；再施加 PSD 激励，用模态叠加法或直接频响法计算危险点应力响应 PSD；随后提取谱矩，选择频域或时域损伤模型；最后用 Miner 准则累积损伤，得到寿命 T = 1/D。判读结果时，不能只看 D 是否小于 1，还要检查危险点位置是否与工程经验一致、网格是否收敛、模态是否遗漏，以及寿命是否对阻尼和 S-N 曲线斜率敏感。