2. 随机振动分析的数学基础

随机振动与正弦振动最本质的区别在于:正弦振动在任意时刻的瞬时值可由幅值、频率和相位完全确定,而随机振动的瞬时幅值无法用确定函数描述,只能用统计量来刻画。Ansys Mechanical 的随机振动分析(Random Vibration)正是建立在频域统计描述之上的:它不求解某一条确定的时间历程,而是以功率谱密度(PSD)作为输入载荷,通过频响函数把输入的统计特性传递到响应,最终给出响应的统计量。因此,理解 PSD、自谱、互谱、均方根值与概率分布,是正确设置载荷、读懂结果的前提。

一、从时域信号到频域统计描述

一个平稳随机过程 x(t) 在时域上表现为幅值、相位都不可预测的波形,但它具有稳定的统计特性:均值、方差、自相关函数不随时间原点变化。对自相关函数 Rxx(τ) 做傅里叶变换,就得到自功率谱密度 Sxx(f),即自谱。它描述的是信号功率在频率轴上的分布密度,单位是“物理量平方每赫兹”,例如加速度 PSD 的单位为 g²/Hz 或 (m/s²)²/Hz。

PSD 与均方值的关系是随机振动分析中最核心的一条公式:

σ² = ∫ S_xx(f) df   (积分区间为分析频带)
RMS = √σ² = √( ∫ S_xx(f) df )

也就是说,PSD 曲线下的面积等于该信号的均方值,开方后即为均方根值 RMS。这条关系决定了随机振动载荷的定义方式:工程上给出的试验谱通常就是一条 PSD 折线(如 10~2000 Hz 内若干控制点),而不是一条时间波形。

二、自谱、互谱与相关矩阵

自谱 Sxx(f) 是实值、非负的偶函数,只反映单个通道的功率分布,丢失了相位信息。当存在多个激励点、且激励之间存在相关性时,必须引入互谱密度 Sxy(f)。互谱是复值函数,其模反映两通道在频率 f 上的共同功率大小,其相位反映两通道之间的相位差。归一化后得到相干函数:

γ²_xy(f) = |S_xy(f)|² / ( S_xx(f) · S_yy(f) ),  0 ≤ γ² ≤ 1

γ²=1 表示两通道完全相干(如同一振动台同时驱动两点),γ²=0 表示完全不相干(如两点各自独立的随机激励)。在 Ansys 中,多点多轴随机振动(PSD 分析中的多点激励)就需要通过互谱矩阵来描述各激励点之间的相关关系;若忽略互谱、只给自谱,则默认各激励相互独立。这是工程中容易出错的地方:把本应相关的两点当作独立激励,会低估或高估响应。

随机振动频域统计描述与传递关系 时域随机信号 x(t) 幅值/相位不可预测,统计特性平稳 自相关 R_xx(τ) 描述信号自身在不同时刻的相关程度 相关 FFT 自谱 S_xx(f) 功率在频率轴上的分布密度 曲线下面积 = 均方值 σ² RMS = √( ∫ S_xx(f) df ) 例:10~2000 Hz 谱面积 0.36 g² → RMS = 0.6 g 互谱 S_xy(f) 与相干函数 γ²_xy(f) S_xy(f) 为复值:模→共同功率,相位→通道间相位差 γ²_xy = |S_xy|² / (S_xx · S_yy),取值 0 ~ 1 γ²=1 完全相干(同一振动台驱动);γ²=0 完全独立 频响传递:S_yy(f) = |H(f)|² · S_xx(f) H(f) 由模态叠加或直接法求得;响应 PSD 再积分开方得到响应 RMS 1σ 结果 ≈ 68.3% 时间不被超越;3σ 结果 ≈ 99.73% 时间不被超越

三、均方根值与概率分布

对平稳高斯随机过程,瞬时值 x(t) 服从均值为零的正态分布,其标准差 σ 等于 RMS。由此得到工程上最常用的判据:

  • 1σ 结果:瞬时值落在 ±1σ 内的概率约 68.3%,即约 31.7% 的时间会被超越;
  • 2σ 结果:约 95.4% 的时间不被超越;
  • 3σ 结果:约 99.73% 的时间不被超越,是疲劳与强度校核中最常用的保守取值。

Ansys Mechanical 随机振动分析默认输出 1σ 位移、应力、加速度结果,用户可在结果详情中把 Scale Factor 设为 2 或 3 得到 2σ、3σ 结果。需要注意:应力 σ 结果与位移 σ 结果都按同一比例缩放,但应力是二阶量,其 3σ 值并非位移 3σ 值的简单线性对应,必须分别提取。

四、教学示例:单自由度系统的 PSD 响应

设一单自由度系统固有频率 fn=100 Hz,阻尼比 ζ=0.02,受平直加速度 PSD 激励 S0=0.01 g²/Hz(10~500 Hz)。其频响函数幅值平方为:

|H(f)|² = 1 / [ (1 - r²)² + (2ζr)² ],  r = f / f_n

在共振点 r = 1 处:|H|² = 1 / (2ζ)² = 1 / 0.04² = 625

响应均方值近似:σ² ≈ (π/2) · f_n · ζ · S_0 · |H|²
             = 1.5708 × 100 × 0.02 × 0.01 × 625
             ≈ 19.6 g²
RMS ≈ 4.43 g

该估算说明:即使输入 PSD 只有 0.01 g²/Hz,在低阻尼共振放大下响应 RMS 也可达数 g。这正是随机振动分析必须关注共振区 PSD 幅值与阻尼设置的原因。在 Workbench 中,若把阻尼比误设为 0.002,|H|² 将放大 100 倍,RMS 估算值会升至约 14 g,结果将严重失真。

五、在 Ansys Mechanical 中的对应设置

随机振动分析(Random Vibration)通常接在模态分析之后,其载荷分支下需要定义 PSD 激励。关键操作与数学量的对应关系如下:

数学量Ansys 中的位置设置要点
自谱 S_xx(f)Random Vibration → PSD G Acceleration输入频率-幅值折线,单位 g²/Hz;注意频率范围覆盖所有关注模态
互谱 S_xy(f)多点 PSD 激励的相关设置默认各点独立;需相关时通过相关矩阵或相干系数定义
均方根 RMS结果 → Detail → Scale Factor1σ 为默认;疲劳校核常取 3σ
概率分布结果后处理判读高斯假设下按 σ 倍数换算超越概率
PSD 输入的单位必须与模型单位制一致。若模型使用 mm-t-s 单位制,加速度 PSD 应换算为 mm²/s⁴/Hz;若沿用 g²/Hz 而未换算,响应量级会相差约 9800² 倍。
随机振动分析只对线性系统成立。若模型中含接触非线性、大变形或材料非线性,频响函数不再适用,需改用瞬态动力学加随机载荷样本,或采用等效线性化处理。

综上,PSD 是随机振动分析的“载荷语言”,自谱描述单点功率分布,互谱描述多点相关性,RMS 与概率分布则把频域结果翻译成可用于强度与疲劳判断的工程指标。掌握这四者的定义、单位与相互换算关系,才能在 Ansys Mechanical 中正确设置激励、合理解读 1σ/3σ 结果,并对异常量级做出快速判断。