4. 信号处理入门

齿轮箱振动信号在时域上表现为复杂的叠加波形:齿轮啮合冲击、轴承滚动体通过、轴系不平衡、结构共振以及背景噪声全部混在一起。若只看时域波形,很难判断哪一频率成分占主导、哪一成分异常。信号处理的任务,就是把时域波形转换到频域,让不同来源的频率成分在频谱上分离开来,从而为故障诊断提供依据。本节围绕傅里叶变换、频谱分析、采样定理、混叠、泄漏与窗函数六个基本概念展开。

一、傅里叶变换:从时域到频域的桥梁

傅里叶变换的核心思想是:任何满足一定条件的周期或非周期信号,都可以分解为一系列不同频率、不同幅值、不同相位的正弦波之和。对于齿轮箱振动这类工程信号,我们实际使用的是离散傅里叶变换(DFT),其快速算法即 FFT。DFT 把长度为 N 的离散时间序列 x[n] 映射为 N 个复数频率点 X[k]:

X[k] = Σ(n=0..N-1) x[n] · e^(-j2πkn/N)
幅值谱 |X[k]|,相位谱 ∠X[k]
频率分辨率 Δf = fs / N = 1 / (N·Ts) = 1 / T
其中 fs 为采样频率,Ts 为采样间隔,T 为采样时长

需要牢记:FFT 得到的是离散谱线,谱线间隔就是频率分辨率 Δf。要分辨两个靠得很近的频率成分(例如齿轮箱中两根轴的转频),采样时长 T 必须足够长,使 Δf 小于两者频率差。

二、频谱分析:齿轮箱诊断的主战场

频谱分析把 FFT 结果以频率为横轴、幅值为纵轴展示。对齿轮箱而言,频谱上通常关注以下特征频率:

  • 轴转频 fr = n / 60(n 为转速 r/min),反映不平衡、轴弯曲等。
  • 齿轮啮合频率 fm = z · fr(z 为齿数),其两侧边带反映齿轮局部损伤。
  • 轴承外圈通过频率 BPFO、内圈通过频率 BPFI、滚动体自转频率 BSF、保持架频率 FTF。
  • 固有频率与共振区,常被冲击激励后形成高频载波。

实际分析中,除了幅值谱,还常用功率谱、包络谱和对数谱。包络谱特别适合轴承早期故障:先对高频共振带做带通滤波,再取包络并做 FFT,故障特征频率会以明显谱线出现。

频谱分析的关键不是“看到峰”,而是把峰与理论特征频率对应起来,并观察其谐波、边带和趋势变化。

三、采样定理与混叠

采样定理(奈奎斯特–香农定理)指出:为了无失真地恢复信号,采样频率 fs 必须大于信号最高频率 fmax 的两倍,即 fs > 2·fmax。fs/2 称为奈奎斯特频率。若信号中含有高于 fs/2 的频率成分,采样后这些高频成分会“折叠”到低频区,形成混叠(Aliasing),在频谱上表现为虚假的低频峰。

混叠对齿轮箱诊断危害很大:一个原本位于 3 kHz 的啮合谐波,若采样率为 4 kHz,可能被折叠到 1 kHz 处,被误判为某个轴承故障频率。防止混叠的根本措施是在 A/D 转换前加抗混叠低通滤波器,把高于 fs/2 的成分滤除,同时根据关注的最高频率合理选择采样率。

常见误区:认为“采样率越高越好”。采样率过高会增大数据量、降低单次分析时长,反而使频率分辨率变差。应按“关注最高频率的 2.5~4 倍”选取采样率,并配合抗混叠滤波。

四、泄漏与窗函数

DFT 隐含假设:所分析的有限长序列是周期延拓的。如果截断处的首尾不连续,周期延拓后就会产生突变,使原本集中的谱线能量向邻近频点扩散,这就是泄漏(Leakage)。泄漏会让弱故障特征被强成分的旁瓣淹没。

解决办法是加窗:用两端渐变为零的窗函数乘以信号,减小截断突变。常用窗函数及其适用场景如下表:

窗函数主瓣宽度旁瓣衰减幅值精度适用场景
矩形窗最窄最差(-13 dB)高整周期采样、瞬态冲击
汉宁窗较宽较好(-31 dB)中一般连续振动、随机信号
汉明窗较宽好(-43 dB)中窄带信号、边带分析
平顶窗最宽好最高精确幅值标定
凯塞窗可调可调可调需要折中主瓣与旁瓣时

选择窗函数的原则是:关注幅值精度选平顶窗;关注频率分辨率和泄漏抑制选汉宁窗;分析冲击成分可用矩形窗或指数窗。加窗会带来幅值损失,需用窗函数的相干增益进行修正。

五、案例:某减速箱轴承外圈故障的频谱判读

某减速箱输入轴转速 1500 r/min,采样率 10 kHz,采样点数 8192,频率分辨率约 1.22 Hz。时域波形出现周期性冲击。加汉宁窗后做 FFT,频谱在 2.8 kHz 附近出现明显共振带,包络谱在 87 Hz 处出现峰值,与理论 BPFO 一致,判定为外圈故障。若不加窗,2.8 kHz 共振带能量泄漏到邻近频点,包络谱中 87 Hz 峰会被抬高底噪掩盖,容易漏判。

六、方法流程与复盘清单

  1. 明确关注最高频率 fmax,按 fs ≥ 2.5·fmax 选择采样率。
  2. 设置抗混叠低通滤波器,截止频率略高于 fmax。
  3. 根据分析目标选择窗函数,记录窗类型与修正系数。
  4. 确定采样时长 T,使 Δf = 1/T 足以分辨目标频率间隔。
  5. 做 FFT,检查频谱是否存在混叠折叠、泄漏旁瓣。
  6. 必要时做包络谱、倒频谱或阶次分析,交叉验证故障频率。
复盘要点:频谱上每一个可疑峰,都应能回答三个问题——它对应哪个理论特征频率?它的谐波和边带是否符合故障机理?改变采样率或窗函数后它是否稳定存在?
齿轮箱振动信号处理链路:从时域到频谱 时域信号 x(t) 齿轮啮合+轴承冲击+噪声 抗混叠低通滤波 截止频率略高于 fmax A/D 采样 fs > 2·fmax N 点,Ts = 1/fs Δf = fs / N 加窗(汉宁/平顶) 抑制泄漏,修正幅值 FFT / DFT X[k] = Σ x[n]·e^(-j2πkn/N) 幅值谱:特征频率定位 fr 2fr fm fm±fr BPFO 混叠:高频折叠为虚假低频 真实 3 kHz 折叠 1 kHz fs/2 边界 折叠 采样率、窗函数与采样时长共同决定频谱的可信度