干旋转机械故障诊断这些年我遇到最多的尴尬场景是设备明明只是台普通离心风机转速却总是稳不住。变频器按工艺曲线加速、减速你用常规FFT频谱去采振动信号屏幕上原本应该立起来的峰值全摊成了一片鼓包再换更长的采样窗口情况没有变好反而更像一段杂乱无章的地震记录。不是传感器坏了也不是采集仪不行而是分析方法用的坐标系不对。这种时候我通常会跟旁边的人说别着急我们换个坐标系看。这个坐标系就是“转轴转过的角度”对应的方法叫阶比分析也叫阶次分析、阶次追踪Order Analysis / Order Tracking。它非常适合旋转机械在转速变化过程中的振动与噪声诊断——风机、泵、压缩机、齿轮箱、发动机、电动机、风电传动链凡是主动轴会变速的设备基本都绕不开它。这篇内容适合正在做设备状态监测、故障诊断、台架试验的工程师和技术人员也适合刚接触振动分析、但已经被“转速一变频谱就糊”困扰的同学。我尽量把原理讲透把实现方法讲清楚最后把实测中容易出问题的细节一起摆出来尽量帮大家少走弯路。1. 转速一波动为什么频谱图就“糊”了1.1 FFT谱分析的前提假设FFT频谱分析的有效性是建立在一个隐含前提上的截取进来的这段信号在时间和频率上基本平稳。也就是说窗口内各振动分量的幅值和频率不能剧烈变化。对旋转机械而言很多振动分量天生就跟着转频走转频一变它们的频率也跟着漂移。举例来说一台风机从1000 r/min匀加速到1100 r/min1倍转频成分在这1秒的采集窗口内频率从16.7 Hz漂移到18.3 Hz。FFT窗口越长频率漂移量就越大能量被摊得越开。最终谱图上的1倍频不再是一根窄峰值而是一团变宽的“鼓包”。谐波阶次越高频率漂移越明显比如齿轮啮合频率如果是23倍频它在同一窗口内漂移幅度就达到了1.6 Hz的23倍谱峰直接找不着。这个问题的本质不是谱图分辨率不够高而是信号的物理特征在做非平稳变化用平稳假设的分析工具去处理它结果必然是“对不齐”。对不齐的谱峰幅值、边带、谐波结构全部失真故障诊断自然无从谈起。1.2 哪些分量“跟转速走”哪些不跟旋转机械里的振动源大致可以分成两类。一类是“阶次分量”频率与转速严格成倍数关系转子不平衡产生1倍频联轴器不对中产生2倍频齿轮啮合产生齒数倍频叶轮叶片通过位置产生叶片数倍频轴承损伤会不同程度地反复冲击并产生与转速相关的调制。这类分量是阶比分析的重点对象因为它们的阶次即频率与转频的比值在物理上具有确定性。另一类是“非阶次分量”结构固有频率、管道声学共振、电网50/60 Hz电磁干扰、变频器载波干扰、流体激振等。这些频率基本不随转速变化或者变化规律与转速无关。在变速工况下把这两类分量混在一起看频率谱阶次分量在“漂”非阶次分量在“原地不动”谱图信息相当混乱。阶比分析的价值就是先把坐标系变换到“以转速为基准的阶次域”让阶次分量稳定下来再去做后续的提取和判断。1.3 时频图为什么不能完全替代阶比分析有人会问既然转速在变那直接用短时傅里叶变换STFT或者时频图把频率随时间变化的曲线画出来不也能看吗时频图确实能提供很多线索但它有几个天然的不足第一时频图横轴是时间诊断要回答的往往是“这个振动是在多少转速、多少阶次下出现的”时间轴需要再经过转速换算属于间接信息。第二STFT的频率分辨率与窗长相互制约转速变化剧烈时频率估计误差大不同阶次的曲线会交叉糊在一起。第三时频图难以直接做幅值平均干扰噪声大故障特征细微变化不容易定量对比。阶比分析相当于做了一次“坐标变换”将时间轴按转轴的角位移重新映射转速变化带来的频率漂移被统一校正随后既可以用FFT角度域谱提取稳定阶次也可以继续做同步平均、包络分析等工具。它和时频图不是替代关系而是更精确的目标定位工具。2. 阶比分析的核心以转轴为时钟重新采样2.1 阶次这个“无量纲”单位先建立最基础的概念。阶次Order定义为振动频率与转轴转频的比值O f / f_r。其中 f 是实际频率单位Hzf_r 是转轴频率也是单位Hzr/min除以60可得。如果一根轴以1800 r/min转动f_r 30 Hz那么1阶就是30 Hz2阶是60 Hz23阶是690 Hz。当转速变化时同样的1阶成分实际频率会从16 Hz变到18 Hz甚至更多但在阶次轴上始终是1阶。这就是为什么阶次域非常适合描述与转速相关的振动成分。旋转机械故障诊断里常见的阶次往往带有明确的物理指向1阶对应不平衡、2阶常与不对中相关分数阶0.5、1.5、2.5往往提示松动或非线性接触齿轮啮合频率的阶次等于齿数叶轮叶片通过频率的阶次等于叶片数。从不抽象每根谱线都能翻译回机械结构本身。2.2 等角度采样把变速信号“掰直”我把阶比分析的核心理解为四个字等角采样。常规数据采集是等时间间隔采样间隔为固定Δt每秒固定多少个点。阶比分析要求的是等角度间隔采样轴每转过固定角度Δθ就采一个点。这样如果某个振动是每转一圈触发一次那么它在重采样后的角域信号里就是一个严格周期信号齿轮啮合每过一个齿触发一次无论转速怎么变一个齿对应角度内的波形特征都稳定重现。也就是说等角度重采样把原本在时间域里被转速伸缩的波形按轴的角度“掰直”了。之后无论是做FFT、同步平均还是包络分析都像是在恒定转速下做一样。这不是复杂技巧而是阶比分析能生效的最重要原因。打个比方你在旋转木马上用固定时间间隔连续拍照背景里的所有景物都是甩动的但如果你把快门触发和旋转的栅格同步让每次拍照都精确发生在同一个角位置照片里相对静止的图案就能被稳定捕捉。阶比分析做的就是“同步快门”。2.3 阶比谱的横坐标与分辨率完成等角度重采样后角域信号在数学上和普通时域信号没有区别只是横坐标从时间变成了“转数”或“圈数角”。对它做FFT结果就是阶比谱横轴是阶次纵轴是幅值。这里有两个关键参数一定要记牢每转采样点数S决定最大可分析阶数。角域“采样率”就是每转S个点按奈奎斯特理论最大分析阶数等于S/2。如果你要看的齿轮啮合阶次是64阶S至少取128实际建议取256以上。每段分析的转数N决定阶次分辨率。阶次分辨率等于1/N。分析10转阶次分辨率为0.1阶可以分辨相距0.1阶的相邻谱峰如果只采了2转分辨率只有0.5阶很多细节就丢了。举个例子轴转速在600~1500 r/min之间变化齿轮齿数为64想同时看到前10阶谐波和齿轮啮合阶次64阶取S1024则最大阶次512阶足够每段分析10转阶次分辨率0.1阶在1000 r/min时对应的频率分辨率约为0.1×16.7≈1.67 Hz。这样的参数配置作为起步完全够用。2.4 转速-阶比彩图Campbell图的使用习惯阶比分析能给出一维的阶比谱也能组合成二维的转速-阶比彩图。横轴是转速纵轴是阶次颜色代表幅值大小。这类图在发动机台架、风电传动链监测中几乎是标配因为一张图就能看全整个升速过程中的振动演化。我拿到变速数据后的习惯是先生成转速-阶比彩图扫一眼哪些彩色条纹竖直稳定说明是固定阶次成分可以把它们沿转速方向做切片得到“阶次跟踪曲线”哪些条纹斜着穿过并和其他条纹交叉出亮点就要警惕共振。比如某阶次的斜线与固有频率水平线交叉的位置往往就是共振转速点。后续再用阶比谱精确定位诊断思路会清楚很多。3. 实现阶次追踪的三条主流技术路线3.1 硬件同步采样最传统的确定性方法阶次分析早期最直接的办法是让采集硬件跟着转速走。具体做法是在轴上安装编码器或键相传感器把转速脉冲接入一个锁相环PLL倍频电路由它控制数据采集卡的采样时钟这样ADC采样率实时跟随转速的倍数变化物理上就实现了每转固定点数。这种硬同步采样方案优点是延时极低、实时性好几乎不占用额外计算缺点也很明显专用硬件贵采样通道受倍频器限制而且锁相环响应有滞后转速突变时容易跟丢。现在大多数便携式振动分析仪和在线监测系统已经不再用纯硬件方案而是走计算路线。3.2 计算阶次跟踪COT工程中最常用计算阶次跟踪Computed Order TrackingCOT是目前现场最主流的方式振动信号和转速脉冲同步采集到同一台采集仪之后在软件里完成角度重构和重采样。完整的流程可以拆成四步同步采集振动通道和转速通道使用同一采样时钟保证时间轴严格对齐。转速曲线重建从键相脉冲或编码器脉冲中提取每转的时间标记再通过插值得到连续转速曲线r(t)。等角度重采样计算角度随时间变化的曲线再反插值求出每个等角度间隔对应的时间点从原始时域信号中取值得到角域信号。阶比谱分析对截取好的角域信号段做FFT得到阶比谱或继续做阶次跟踪、同步平均。COT对硬件要求不高普通多通道采集卡即可算法复杂度也在可控范围。但它有一个隐含假设在短时间段内转速变化是平滑的。如果转速曲线剧烈跳动或者键相脉冲缺失、抖动重采样结果会明显失真这一点到第5部分再展开。3.3 Vold-Kalman阶次跟踪更高级的去噪与分离当目的不只是“看一根阶次谱线”而是要把某个阶次成分的时域波形单独提取出来比如去掉齿轮啮合大分量之后再做轴承包络分析就可以考虑Vold-Kalman阶次跟踪VKF。VKF的核心思路和COT不同它不依赖等角度重采样的“几何补偿”而是建立状态空间模型把目标阶次的幅值和相位设定为待估状态通过卡尔曼滤波和平滑从原始时域信号中分离出该阶次分量。它对转速估计误差的敏感度低于COT适合转速变化很快、转速曲线很难精确测量的场景同时可以同时跟踪多个阶次也可以任意频率间隔分离。代价是算法参数调节比较复杂计算量也更大。滤波带宽参数决定了阶次信号允许的带宽设小了会漏掉真实调制信息设大了会把邻近阶次的能量勾进来。我个人的建议是常规故障诊断先用COT拿到阶比谱性价比最高确实需要做时域重构、分离干扰阶次时再上VKF不要一上来就在参数里挣扎。3.4 一个最小的计算阶次跟踪示例很多人对原理知道了但一上来还是不知道怎么动手。我写一个便于理解的最小示例用Python和NumPy演示等角度重采样的核心步骤。这里假设我们已经通过键相脉冲算出了和振动信号时间轴对齐的瞬时转速数组rpm单位r/min采样率是fs。import numpy as np from scipy.interpolate import CubicSpline fs 25600.0 # 采样率单位Hz spp 1024 # 每转采样点数决定最高阶数 t np.arange(len(vib)) / fs # 1. 由瞬时转速积分得到角度随时间的变化 theta np.cumsum(2.0 * np.pi * rpm / 60.0 / fs) # 2. 生成目标等角度序列 total_turns (theta[-1] - theta[0]) / (2.0 * np.pi) turn_idx np.arange(int(total_turns * spp)) / spp theta_target theta[0] turn_idx * 2.0 * np.pi # 3. 用插值把原始振动信号映射到等角度时刻 interp CubicSpline(theta, vib, extrapolateFalse) vib_angle interp(theta_target) # 4. 对指定圈数的角域信号做FFT得到阶比谱 segment_len 10 * spp # 10转一个分段 seg vib_angle[:segment_len] * np.hanning(segment_len) order_axis np.arange(0, segment_len // 2) / 10 # 0~0.5*S阶 amp_spectrum np.abs(np.fft.rfft(seg)) * 2 / (spp * 10)这段代码把最核心的操作表达得很直观先通过转速积分得到角度曲线再从角度反查时间在原始信号上插值得到角域信号。工程实现里还要处理边界效应、非法插值点、窗函数幅值修正、直流分量剔除等细节但原理就是这三行积分角度、等角插值、FFT。必须强调这个示例假设转速曲线已知且单调连续。实际现场如果转速突变、反转或者键相脉冲丢失需要先对转速数据做滤波和异常剔除否则重采样后的阶比谱全是假峰。3.5 三条技术路线怎么选技术路线核心原理优点缺点适用场景硬件同步采样锁相环倍频控制ADC采样时钟实时性好无需复杂计算硬件专用成本高转速突变响应慢老式专用测试系统、生产监控计算阶次跟踪COT同步采集后用转速积分做等角度插值硬件通用灵活算法可调对转速曲线平滑度有要求绝大多数现场测试、在线监测Vold-Kalman状态空间模型估计各阶次幅值相位转速误差容忍度高可分离单阶参数调节复杂计算量大单阶提取、共振隔离、时域重构4. 典型故障在阶比谱上的“指纹”4.1 基础类故障不平衡、不对中、松动先看最常见的低阶故障特征。转子不平衡阶比谱上以1阶为主幅值通常随转速升高而增大径向测点最明显。不对中常呈现1阶和2阶同时增大且2阶比例突出有时伴随3阶和4阶谐波。结构松动或间隙配合不良会在1阶的整数分频处出现峰值比如0.5阶、1.5阶、2.5阶严重时出现连续分数阶谱线。这些基础特征在恒转速频谱里也能看但在变速工况下阶比谱的最大优势是能稳定给出“几阶”这个不变信息。以前我判断一台加速过程显著振动的水泵用频谱看只能看到一大堆变化峰用阶比谱一看主峰稳定在2阶2阶高过1阶再结合轴向测点振动大直接锁定对中不良省了很多弯弯绕。4.2 齿轮故障啮合阶次与边带结构齿轮箱的故障诊断阶比分析几乎是必修课。齿轮正常啮合时阶比谱上辊子因为啮合冲击而产生的峰值出现在啮合频率对应的阶次也就是齿数Z还有2Z、3Z等高次谐波。齿轮齿面有局部损伤时每次轮齿进入啮合都会受到冲击振动的幅值和相位会被所在轴的转速调制于是啮合阶次两侧会出现大量边带边带间隔恰好是故障轴转频对应的阶次多数情况下就是1阶。当断齿或严重磨损时啮合频率谐波数量增多边带数量增加、幅值升高。更典型的是边带不再对称啮合阶次周围出现“成簇”的谱峰严重时甚至会掩盖啮合峰本身。我遇到过一次挤塑机齿轮箱故障恒转速频谱只能看到一大片难以辨认的高频鼓包但阶比谱上啮合阶次两侧出现了间距为1阶的密集边带最高边带超过啮合主峰再结合时域里明显的周期性冲击间隔停机开箱果然是齿面剥落。另外一个非常实用的工具是同步平均在角域信号里按整数转数截取并平均能有效压制与轴转不同步的随机冲击和干扰齿轮啮合波形会变得非常清晰。很多齿轮故障特征是从时域同步平均波形里先发现再去阶比谱验证的。4.3 滚动轴承变速条件下的包络阶次追踪滚动轴承的故障特征频率比如外圈BPFO、内圈BPFI、滚动体BSF、保持架FTF本身就不是严格的整数阶次而是由滚道直径、滚动体直径、接触角等几何参数计算出的比例值多数落在3~10阶之间。轴承滚动体存在随机滑移频率也不是完全稳定的所以直接在变速轴上进行高阶次谱分析效果往往不如低转速恒定工况理想。工程上更推荐的做法是“包络阶次分析”先对原始振动信号做带通滤波把耦合到轴承固有共振频带的冲击能量提取出来再对包络信号做等角度重采样最后对重采样后的包络信号做FFT得到包络阶比谱。此时轴承故障阶次就非常明显了。我曾经在一台转速按斜坡变化的风机驱动电机上用常规加速度阶比谱没能看到清晰的轴承内圈故障特征角域包络谱却明确在6.6阶附近出现幅值很高的峰值。事后拆检发现内圈滚道确实有一条贯穿性剥落带。变速工况下轴承诊断请一定要把“包络”和“阶比”结合用别单看原始阶比谱。4.4 叶片与叶轮叶片通过频率的阶次叶轮机械里叶片通过频率BPF的阶次等于叶片数。一台9叶片的离心风机BPF的阶次就是9阶一台有17个叶片的轴流风机BPF直接出现在17阶。这个特征在阶比谱上很好识别几乎不会认错。叶片安装不均匀、叶片变形或局部破损时除了整数阶叶片通过分量还会出现分数阶分量比如每转产生一次的非平衡激励在1阶、2阶处增强。气蚀、旋转失速等流体激励则往往会让BPF附近的宽频噪声成分上升表现为阶次谱线上方的“媳峰”或“噪声着地”。叶片故障的阶比分析重点就是看BPF主峰值和高阶谐波的增长趋势以及特有的边带宽度变化。我发现只要结合不同负荷工况下的多组阶比谱对比叶片早期损伤识别并不难。4.5 阶次切片找共振的捷径阶比谱是静态切片而转速-阶比彩图可以沿某一个确定阶次做“切片”得到该阶次随转速变化的幅值曲线。比如风机升速过程中2阶幅值在1400 r/min左右突然飙升随后又回落这大概率是2阶激励与某个结构固有频率发生共振。沿阶次切片能快速划定共振转速范围再配合模态实验或运行变形分析做确认。共振问题在变频驱动设备上非常普遍。转速一旦跨过共振区间哪怕设备平衡没有问题振动也可能超标到报警。我处理过不少“换了联轴器还是抖”的现场案例最后都是通过阶次切片发现是机组底座固有频率被阶次激励激发而非平衡或对中问题。5. 现场实测与数据处理中容易踩的坑5.1 键相信号质量决定一切阶比分析依赖转速通道。键相传感器安装不好后续算法再高明也无济于事。常见的键相测量方式有反光贴纸加光电传感器、电涡流测轴键槽、接近开关测齿轮齿、编码器输出脉冲。无论哪种都要确保脉冲计数完好、边沿清晰、不丢失脉冲。轴面有油污、反光贴翘边、传感器支架松动都会让键相脉冲时间标志带干扰重采样时刻随之抖动阶比谱上便会产生一批伪谱线看起来像“多加了好几阶”。实测中还有一个容易被忽略的问题键相传感器和振动传感器必须保持同一采集时钟同步。异步采集时哪怕时间标签只差几个毫秒转速和振动之间的相位关系也会错乱阶次谱的峰值位置可能偏移边带结构失真。用数据采集软件时请先确认“同一板卡同一时钟”或做了时钟同步。5.2 抗混叠滤波必须在重采样之前完成等角度重采样后角域“采样率”等于每转采样点数S最高可分析阶数是S/2。如果原始振动信号里存在高于S/2对应的阶次的高频分量重采样后就会混叠到低阶区形成假峰且很难分辨。我见过一位工程师测一台高速电主轴S取256仪器没有自动滤波几千赫兹的固有共振信号混叠后在阶比谱上2~3阶附近多出来一串假峰值差点把故障误判为严重不对中。后来加了截止频率按最大分析阶数设计的低通滤波器假峰消失真正的问题才浮现出来。项目里建议把抗混叠滤波放在转速重采样之前且截止频率要结合最大转速和最高关注阶次来设别偷懒直接用固定低通那样可能把真实高频故障信息也滤掉了。5.3 转速曲线的插值精度不能敷衍COT的转速曲线从键相脉冲时间戳重建。如果每转只有一个键相脉冲转速在两个脉冲之间默认按某种方式插值通常用线性或样条插值。转速变化加速度较大时线性插值会造成角位置估计偏差阶比谱峰形变宽、幅值变低严重时出现虚假边带。我的做法是转速变化平缓时线性插值够用转速变化快或测试精度要求高时至少使用三次样条插值如果键相脉冲频率很密编码器几百脉冲每转可以把脉冲角度累加再插值到较高精度。重采样边界处也就是开头和结尾不完整转数建议丢弃一半转左右不要强行补零否则会引入边缘频谱泄漏。5.4 窗函数和每帧转数需要按目标折中角域信号和时域信号一样有频谱泄漏问题。一般建议对每个分析分段加汉宁窗或汉明窗但分段转数N要结合目标阶次分辨率来选。分析转数越多分辨率越细但消耗时间越长变速过程里转速变化也越大可能使一个分段内的信号本身“非平稳”。我的习惯是诊断齿轮和轴承分段转数取8~16转阶次分辨率0.06~0.125阶只看基础故障和共振判断4~8转就够。如果转速变化非常快宁可通过多次分段加平均来压低噪声也别一味延长分段。角域FFT的幅值修正要用窗函数的相干增益系数做归一化否则幅值会偏低这一点和时域FFT完全相同。5.5 留意那些“不跟转速走”的分量阶比谱擅长呈现“跟转速走”的阶次成分但机械系统里那些“不跟转速走”的分量同样重要。结构固有频率、固定频率的电磁干扰、传感器谐振、管道声模态在转速-阶比彩图上表现为斜向或其他形态的条带在阶比谱里则可能以随机峰值甚至宽噪声的形式出现。如果把所有谱峰都按阶次去解释很容易误判。正确的做法是交叉验证把同样的数据画成频率谱和阶比谱并排对比看谱峰在两个域里的形状。跟着转速走、在阶次域里稳定、在频率域里漂移的是真正的阶次成分反之则更可能是固定频率干扰或共振。诊断报告里最好同时给出转速-阶比图和频谱图让结果经得起复核。5.6 低速设备要优先考虑积分和包络转速很低比如几十r/min的大型回转设备一倍频加速度幅值非常小阶比谱上的基础故障特征很容易淹没在底噪里。这种情况下我建议在重采样之前或之后对信号做一次积分把加速度信号转成速度或位移提升低频段的信噪比。积分的代价是会放大低频扰动和直流漂移必须配合高通滤波抑制趋势项。轴承低速工况下的故障诊断也一样原始加速度阶比谱基本看不出东西而包络阶比谱通常效果明显得多。这套“低速必用积分/包络”的经验帮我处理过不少回转窑、大型磨机轴承的疑难问题。5.7 把阶比分析和传统手段叠加使用最后一条经验也是最想强调的阶比分析不是孤立的诊断工具最好与传统手段叠加使用。先做“三件套”起步——时域波形看冲击和调制频谱和包络谱看主要频率带转速-阶比图定位特征阶次再做重点——对可疑阶次做切片对故障轮齿做同步平均必要时再上Vold-Kalman分离单阶信号重构冲击波形。这种“从粗到细”的流程比直接怼着一张阶比谱下结论要可靠得多。阶比谱上的一根高幅值谱线也许只是转速估计误差叠加窗函数泄漏造成的伪峰只有多域交叉验证才能确定拿到的确实是故障“指纹”。最后再分享一个个人习惯我拿到变速振动数据第一件事永远是先看转速曲线确认脉冲干净、转速单调、区间稳定再做后续分析。项目里踩过太多没看转速就开始算阶比的亏比如启动和停机阶段转速反转、键相贴片脱落导致脉冲丢失这些脏数据算出来的阶比谱再漂亮也是垃圾进垃圾出。阶比分析的本质是“把转轴当作主时钟”那么首先要保证的就是这个时钟走得准。做阶比分析的时间越久越会觉得数据质量、键相信号、重采样参数这些“不起眼的细节”往往才是项目成败的关键。