4. 随机振动疲劳损伤机理

随机振动疲劳寿命仿真的核心任务，是把一段“看起来杂乱无章”的宽带或窄带随机载荷，折算成可比较、可累积的损伤量。要完成这个折算，需要四个关键环节：用 S-N 曲线描述材料在不同应力幅下的寿命能力；用雨流计数把连续随机应力历程拆成一个个闭合的应力循环；用 Miner 线性累积损伤把不同幅值循环的损伤叠加起来；当载荷用功率谱密度（PSD）给出时，则改用频域疲劳损伤理论直接由谱参数估算损伤。理解这四个环节各自的假设与边界，比记住公式更重要。

S-N 曲线：材料寿命能力的基线

S-N 曲线描述恒幅对称循环下应力幅 S_a 与失效循环数 N 的关系。工程上最常用的是幂律形式：

S_a = S_f' · (2N)^b        （Basquin 形式，应变-寿命中常用）
N · S_a^m = C               （应力-寿命幂律形式，m、C 为材料常数）

在双对数坐标下，幂律形式是一条直线，斜率由 m 决定。以某铝合金焊接接头为例（示例参数，非实测）：m = 3.5，当 S_a = 100 MPa 时 N = 2×10⁵ 次，则 C = 2×10⁵ × 100^3.5 ≈ 2.24×10¹²。若应力幅降到 50 MPa，可算得 N ≈ 2.53×10⁶ 次，寿命提高约 12.6 倍——这正是随机振动中“小幅值循环数量巨大、但单次损伤极小”的根源。

S-N 曲线的适用边界必须明确：它来自恒幅、对称、等幅循环试验，平均应力、加载频率、表面状态、残余应力、环境温度都会使其偏移。工程中常用 Goodman 或 Gerber 修正把平均应力折算到等效对称幅值，再查 S-N 曲线。此外，钢材常存在疲劳极限，低于该幅值的循环理论上不产生损伤；铝合金等则通常没有明显疲劳极限，需人为设定截止幅值或截止寿命（如 10⁷ 次），否则小循环会被无限累积。

雨流计数：从随机历程到应力循环

随机应力时程 σ(t) 是一段连续波动的曲线，不能直接代入 S-N 曲线。雨流计数的做法是：把时程旋转 90°，让“雨水”沿峰谷依次流下，凡遇到来自上方的更大峰值（或来自下方的更深谷值）即停止，从而把复杂波形拆成一系列闭合的滞回环，每个环对应一个应力幅 S_a 和均值 S_m。

实际操作中，先做峰谷提取（只保留极值点），再做雨流计数，得到“幅值-均值-循环数”的二维计数矩阵。以一段 60 s、采样率 2000 Hz 的加速度响应为例（示例）：经峰谷压缩后得到约 1.2 万个极值点，雨流计数后归并为 64 个幅值区间 × 16 个均值区间，总循环数约 5800 次。这个矩阵就是后续 Miner 累积的输入。

雨流计数的边界在于：它假设材料对循环的记忆是“封闭的”，即每个闭合环独立计一次损伤，不考虑加载顺序和过载效应。对于有显著过载迟滞、裂纹闭合效应或腐蚀环境的场景，雨流计数给出的损伤会偏保守或偏乐观，需要配合顺序效应修正。

Miner 线性累积损伤：把循环叠加成寿命

Miner 法则假设：每个应力幅 S_a,i 对应寿命 N_i，该幅值下实际循环 n_i 次，则损伤为 D_i = n_i/N_i；总损伤 D = ΣD_i，当 D = 1 时认为失效。若总时长为 T，则寿命 T_f = T/D。

继续上面的示例：若雨流矩阵中 80 MPa 幅值有 1200 次循环，查 S-N 曲线得 N = 4.2×10⁵，则 D₈₀ = 1200/4.2×10⁵ ≈ 2.86×10^-3；50 MPa 有 2600 次，N = 2.53×10⁶，D₅₀ ≈ 1.03×10^-3；其余小幅值循环贡献约 4×10^-4。合计 D ≈ 4.3×10^-3，对应寿命约 60 s / 4.3×10^-3 ≈ 1.4×10⁴ s ≈ 3.9 小时。这个结果说明：在随机振动中，少数大幅值循环往往贡献了主要损伤，而大量小幅值循环虽然数量多，但单次损伤极小。

Miner 法则的边界同样清晰：它忽略加载顺序、忽略损伤非线性累积、忽略残余应力与裂纹扩展阶段，且对 S-N 曲线和截止幅值非常敏感。工程上常把 D = 1 改为 D = 0.3～0.5 作为设计判据，以覆盖分散性。

频域疲劳损伤理论：从 PSD 直接估算

当随机载荷以 PSD G(f) 给出时，不必先做时域雨流计数，可直接由谱矩 m_n = ∫ fⁿ G(f) df 计算损伤。常用方法有窄带近似、Wirsching-Light 修正、Dirlik 模型和 Lalanne 方法。以窄带近似为例：

m0 = ∫ G(f) df          （0 阶谱矩，对应均方值）
m2 = ∫ f² G(f) df       （2 阶谱矩）
m4 = ∫ f⁴ G(f) df       （4 阶谱矩）
ε = sqrt(1 - m2²/(m0·m4))   （带宽参数）
E[P] = sqrt(m4/m2)          （峰值率）
D = (E[P]·T / C) · (√(2m0))^m · Γ(1 + m/2)

其中 Γ 为 Gamma 函数，C、m 来自 S-N 曲线。窄带近似假设应力幅服从 Rayleigh 分布，适用于 ε 接近 0 的窄带响应；当 ε 增大（宽带），窄带近似会高估损伤，需用 Dirlik 等修正模型。Dirlik 模型通过一个指数加两个 Rayleigh 分布的混合概率密度拟合幅值分布，在 ε = 0～1 范围内精度较好，是目前工程软件中应用最广的频域方法之一。

频域方法的适用边界：它要求载荷近似平稳高斯随机过程；若响应存在明显非高斯、非平稳、冲击或间歇性特征，谱矩和 Rayleigh 假设不再成立，应回到时域雨流计数或采用非高斯修正。此外，频域方法对 PSD 的频率分辨率和截止频率敏感，低频段缺失会低估 m0，高频段截断会低估 m4，从而影响峰值率和带宽参数。

四个环节的假设与边界对照