4. 主轴轴承故障的振动信号基础
风电主轴轴承的振动诊断，本质上是从加速度传感器采集到的时域波形中，识别出与轴承几何结构和转速相关的周期性冲击成分。要读懂这些成分，必须先理解振动是怎么产生的、特征频率怎么算、调制现象如何形成，以及低速重载工况给信号带来了哪些特殊困难。
一、振动产生机理：从缺陷冲击到结构响应
主轴轴承在正常运转时，滚动体与滚道之间是连续、平稳的接触，振动主要表现为低幅值的随机噪声和转频相关的微弱成分。一旦滚道或滚动体出现剥落、点蚀、压痕等局部缺陷，情况就发生变化：每当滚动体滚过缺陷位置，接触刚度瞬间下降，产生一次短促的冲击力。这个冲击力激励轴承座、主轴和机舱结构，形成一系列宽带脉冲。
关键在于，冲击本身是宽带的，但轴承各部件（内圈、外圈、滚动体、保持架）都有各自的固有频率。冲击能量会激发这些固有频率，使传感器测到的信号呈现为“高频载波被低频冲击周期性调制”的形态。因此，故障诊断通常分两步：先在包络谱中找出冲击的重复频率，再对照特征频率判断缺陷位置。原始频谱中往往看到的是固有频率附近的高频共振峰，而真正的故障信息藏在包络里。
二、故障特征频率计算（BPFO、BPFI、BSF、FTF）
轴承的四个特征频率由几何参数和转速共同决定。设滚动体个数为 n，滚动体直径为 d，节径为 D，接触角为 α，轴转频为 fr，则：
BPFO = (n/2) × fr × (1 − (d/D)·cosα) 外圈故障频率
BPFI = (n/2) × fr × (1 + (d/D)·cosα) 内圈故障频率
BSF = (D/(2d)) × fr × (1 − ((d/D)·cosα)²) 滚动体故障频率
FTF = (1/2) × fr × (1 − (d/D)·cosα) 保持架故障频率以某 2.5 MW 风机主轴轴承为例，假设 n=18，d=45 mm，D=320 mm，α=0°，主轴转速 15 r/min，则 fr=0.25 Hz。代入计算：d/D=0.1406，cosα=1，于是 BPFO≈1.93 Hz，BPFI≈2.57 Hz，BSF≈0.106 Hz，FTF≈0.107 Hz。可以看到，这些频率都极低，BPFO 和 BPFI 仅比转频高几倍，BSF 和 FTF 甚至低于转频。这是风电主轴轴承区别于高速轴承的显著特点。
需要强调：上述数值是教学示例，实际计算必须使用轴承厂家提供的精确几何参数和实测转速。转速波动会直接导致特征频率漂移，因此计算时应以平均转速为基准，并考虑 ±5% 的容差带。
三、转频与调制现象
实际信号中，故障冲击很少以单一频率出现，而是伴随调制。最常见的两种调制是：
- 幅值调制：当内圈存在缺陷时，缺陷随内圈旋转进入、退出承载区，冲击强度周期性变化，调制频率等于转频 fr 或其倍频。包络谱中会在 BPFI 两侧出现 BPFI ± fr 的边带。
- 频率调制：转速波动或载荷变化使冲击间隔轻微变化，表现为特征频率峰附近的展宽或不对称边带。
此外，主轴轴承还常受到叶片通过频率（3P，即三叶片风机的 3 倍转频）和塔影效应的调制。因此，在包络谱中看到 BPFO 附近有 fr 间隔的边带族时，应优先怀疑外圈或内圈故障，而不是简单归因于转速不稳。
四、风电低速重载工况下的信号特点
风电主轴轴承的转速通常只有 10～20 r/min，转频 0.17～0.33 Hz，而故障特征频率往往低于 5 Hz。这带来四个现实困难：
- 冲击能量弱：低速下冲击加速度幅值小，容易被背景噪声淹没，信噪比低。
- 周期长、样本少：一个冲击周期可能长达数秒，短时采样难以覆盖多个周期，频谱分辨率不足。
- 载荷波动大：风速和风向变化导致载荷持续波动，使特征频率漂移、边带模糊。
- 低频段干扰多：塔架摆动、叶片旋转、主轴弯曲等低频成分与故障特征频率重叠，容易误判。
应对策略上，工程中常采用高灵敏度加速度传感器、较长的采样时间（如 60 s 以上）、包络解调分析和阶次跟踪等方法。其中阶次跟踪以转速为参考，把时域信号重采样为角度域，可有效消除转速波动的影响，使特征频率在阶次谱中稳定呈现。
五、从信号到判断：一个可复现的分析流程
把上述知识串起来，实际诊断可按以下顺序推进。第一步，确认转速：从 SCADA 或转速传感器获取主轴平均转速，计算 fr 及四个特征频率，并标注 ±5% 容差带。第二步，采集足够长的振动数据：建议采样率不低于 5 kHz，采样时长不少于 60 s，确保覆盖至少 10 个冲击周期。第三步，做包络解调：对原始信号带通滤波（中心频率选在结构固有频率附近，如 1～3 kHz），再取 Hilbert 包络并做 FFT，得到包络谱。第四步，在包络谱中查找与 BPFO、BPFI、BSF、FTF 对应的峰值，并观察是否存在 fr 间隔的边带。第五步，结合趋势判断：单次峰值可能来自偶然冲击，应对比历史数据，确认该频率成分是否持续增长。