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) 服从均值为零的正态分布,其标准差 σ 等于 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 Factor | 1σ 为默认;疲劳校核常取 3σ |
| 概率分布 | 结果后处理判读 | 高斯假设下按 σ 倍数换算超越概率 |
综上,PSD 是随机振动分析的“载荷语言”,自谱描述单点功率分布,互谱描述多点相关性,RMS 与概率分布则把频域结果翻译成可用于强度与疲劳判断的工程指标。掌握这四者的定义、单位与相互换算关系,才能在 Ansys Mechanical 中正确设置激励、合理解读 1σ/3σ 结果,并对异常量级做出快速判断。