2. 随机振动基础与频域描述

随机振动与正弦定频、扫频振动最大的区别在于：它没有确定的瞬时幅值和相位，任一时刻的加速度值都无法用公式精确预测，只能用统计量描述。工程上处理随机振动时，几乎全部工作都在频域完成，因为时域波形看似杂乱无章，而频域上的功率谱密度（PSD）却稳定、可重复、便于比较。理解 PSD、自谱、互谱、均方根、带宽与中心频率这几个量，是后续把随机载荷转成疲劳损伤的前提。

从时域统计到频域谱密度

随机振动信号通常假设为平稳各态历经过程，即一次足够长的样本就能代表整个随机过程。描述它的基本统计量是均值和方差：均值反映直流分量，方差反映交流功率。对加速度信号 a(t)，其均方值定义为

E[a²] = lim(T→∞) (1/T) ∫₀ᵀ a²(t) dt

均方根 RMS 就是均方值的平方根，单位与加速度相同（如 g 或 m/s²）。在随机振动试验中，总均方根加速度（Grms）是最常用的量级指标，它直接决定振动台推力需求，也粗略反映结构承受的能量水平。

把均方值按频率分解，就得到功率谱密度。PSD 表示单位频率带宽内所含的均方值，单位是 (m/s²)²/Hz 或 g²/Hz。对单边谱，均方根与 PSD 的关系为

Grms = sqrt( ∫[f1→f2] W(f) df )

其中 W(f) 为单边 PSD，积分区间为分析带宽。这条公式是随机振动疲劳仿真的核心桥梁：PSD 曲线下的面积就是均方值，开方即 Grms。若 PSD 为平直谱 W₀，带宽为 Δf，则 Grms = sqrt(W₀·Δf)。例如 W₀ = 0.01 g²/Hz、Δf = 2000 Hz，则 Grms = sqrt(20) ≈ 4.47 g。

自谱、互谱与相干函数

自功率谱（Auto-PSD）描述同一信号自身在各频率上的能量分布，是单点随机载荷的标准表达。互功率谱（Cross-PSD）描述两个信号之间的频域相关关系，它是一个复数，既有幅值又有相位。互谱在疲劳仿真中的典型用途有两个：一是多输入点激励时判断各激励源是否相关，二是通过频响函数与互谱反推载荷。若两个测点信号完全独立，互谱趋近于零；若完全相关，互谱幅值等于两者自谱的几何平均。

相干函数 γ²(f) 由自谱和互谱构造：

γ²(f) = |Sxy(f)|² / (Sxx(f)·Syy(f))

其取值在 0 到 1 之间。γ² 接近 1 说明该频率下两信号线性相关强，互谱可信；γ² 明显偏低则说明存在噪声、非线性或其它未测激励，此时用互谱做载荷识别会失真。在实测数据预处理阶段，通常要求关键频段相干函数大于 0.8 才认为数据可用。

带宽与中心频率

随机振动谱形常用带宽和中心频率来概括。对于窄带随机振动，能量集中在中心频率 f₀ 附近很窄的范围内，其时域波形近似为幅值缓慢变化的“拍”，疲劳损伤与中心频率处的共振密切相关。对于宽带随机振动，能量分布在很宽的频带上，时域波形更杂乱，但高频小幅值循环对损伤的贡献通常有限。

中心频率常取功率谱的重心频率：

f_c = ∫ f·W(f) df / ∫ W(f) df

带宽则可用半功率带宽或有效带宽表示。在疲劳仿真中，带宽决定了应采用哪种损伤计算模型：窄带信号可用 Rayleigh 分布峰值模型，宽带信号则需用 Dirlik、Lalanne 或 Steinberg 三区间法等考虑带宽修正的方法。若误把宽带谱当窄带处理，寿命预测可能偏保守或偏危险，因此带宽判断是建模前必须确认的一步。

PSD 在疲劳仿真中的物理意义

PSD 进入疲劳仿真的路径是：PSD → 频响函数 → 应力 PSD → 谱矩 → 损伤。结构在随机激励下的应力响应功率谱为

Sσσ(f) = |H(f)|² · SFF(f)

其中 H(f) 为从激励到应力的频响函数，SFF 为激励 PSD。得到应力 PSD 后，计算其前四阶谱矩 m₀、m₁、m₂、m₄，即可得到均方根应力、峰值率、不规则因子等参数，再代入相应损伤模型。谱矩定义为

mₙ = ∫ fⁿ · Sσσ(f) df

m₀ 是应力均方值，其平方根为应力 RMS；m₂/m₀ 给出零穿越频率；m₄/m₂ 给出峰值频率。不规则因子 γ = m₂ / sqrt(m₀·m₄) 用于衡量信号接近窄带的程度，γ 越接近 1 越接近窄带。这些量全部来自 PSD 曲线，因此 PSD 的频段截断、谱线分辨率和单位一致性会直接影响寿命结果。

一个可复算的示例

设某安装点随机振动条件为：10～100 Hz 平直谱 0.005 g²/Hz，100～500 Hz 平直谱 0.02 g²/Hz，500～2000 Hz 平直谱 0.003 g²/Hz。分段计算均方值：

第一段：0.005 × (100-10)  = 0.45 g²
第二段：0.02  × (500-100) = 8.00 g²
第三段：0.003 × (2000-500)= 4.50 g²
总均方值 = 12.95 g²
Grms = sqrt(12.95) ≈ 3.60 g

若结构一阶模态为 320 Hz，落在第二段高量级区，则该模态被充分激发，应力 PSD 会在 320 Hz 附近出现明显峰值，疲劳损伤主要由这一段贡献。反之若一阶模态为 80 Hz，处于低量级区，损伤会显著降低。这说明 PSD 与结构频响的匹配关系，比单纯看总 Grms 更能决定寿命。

把上述关系串起来看，随机振动疲劳仿真的频域流程可以概括为：确认 PSD 类型与单位，检查相干函数，计算 Grms 与谱形特征，结合频响函数得到应力 PSD，再提取谱矩并选择匹配带宽的损伤模型。每一步的输入都来自上一步的频域量，任何一处单位或带宽处理错误，都会沿链条放大到寿命结果。因此在实际项目中，建议把 PSD 的频段、分辨率、单位、Grms 和中心频率记录在同一张载荷卡上，作为后续仿真与试验对比的基准。