跳转到内容
新建笔记

旋转机械阶次跟踪:测量链与参数选择

阶次跟踪用于区分“随转速一起变化的振动分量”和“固定频率附近的响应”。这篇笔记侧重测量链、分辨率和工程设置;基础波形、FFT 幅值约定与扫频推导见阶次、频率与赫兹。

1. 先选参考轴,再定义阶次

跳转到“1. 先选参考轴,再定义阶次”

记参考轴转速为 R(t)R(t),单位 rpm;转频为 r(t)=R(t)/60r(t)=R(t)/60,单位转/秒。某振动分量的瞬时频率为 f(t)f(t),它相对于该轴的阶次为

O(t)=f(t)r(t)=60f(t)R(t).O(t)=\frac{f(t)}{r(t)}=\frac{60f(t)}{R(t)}.

一阶表示每转振动一次,二阶表示每转振动两次。阶次可以是整数、分数或其他实数;非同步成分的阶次还可能随转速变化。同一信号相对于不同轴会有不同阶次,所以报告必须注明参考轴和传动比。

结果横轴回答的问题
频谱频率 Hz某段数据在哪些频率有能量?
阶次谱阶次每转重复多少次的成分比较强?
阶次跟踪曲线转速或时间选定阶次的幅值、相位如何随工况变化?
阶次瀑布图/色谱阶次与转速升降速过程中,哪些成分与转动同步?

一阶峰值本身不能直接证明不平衡;安装偏心、结构响应以及其他激励也可能贡献一阶成分。诊断还需要方向、相位、工况和趋势证据。

2. 从固定时间采样到等角度采样

跳转到“2. 从固定时间采样到等角度采样”

累计转数定义为

q(t)=q(t0)+∫t0tr(τ) dτ,θ(t)=2πq(t).q(t)=q(t_0)+\int_{t_0}^{t}r(\tau)\,d\tau, \qquad \theta(t)=2\pi q(t).

若选定区间内 q(t)q(t) 单调,就能求出目标转角对应的时间 t(q)t(q),再构造角域信号 xq(q)=x(t(q))x_q(q)=x(t(q))。设每转重采样 MM 点:

qj=q0+jM,tj=t(qj),xq[j]=x^(tj).q_j=q_0+\frac{j}{M},\qquad t_j=t(q_j),\qquad x_q[j]=\widehat{x}(t_j).

其中 x^\widehat{x} 表示由原始采样数据重建/插值得到的信号。在升速时,相邻 tjt_j 越来越近;在降速时越来越远。物理 ADC 可以保持固定采样率,变化的是数字重采样的输出时刻。

转速脉冲确定等转角输出时刻,固定采样ADC数据经数字重采样生成角域信号

Crystal Instruments 原图:下方是振动采集链,上方是转速脉冲、周期估计和等角度采样时刻计算。图中的固定采样率 ADC 与变频输出采样命令属于两个不同环节。图片版权归 Crystal Instruments。

图中以二阶曲线拟合估计脉冲间角度,这是特定实现方法,不是所有阶次算法的共同要求。仅有一次/转脉冲时,脉冲之间的转速变化必须靠模型估计;快速扭振和转速突变可能需要更密集的角度测量。

3. 时间分辨率与阶次分辨率分别计算

跳转到“3. 时间分辨率与阶次分辨率分别计算”

固定采样率为 fsf_s,采集 NN 点时,DFT 的观察时长与频率栅格间隔为

Δt=1fs,T=Nfs,Δf=fsN=1T.\Delta t=\frac{1}{f_s},\qquad T=\frac{N}{f_s},\qquad \Delta f=\frac{f_s}{N}=\frac{1}{T}.

这里的 TT 是 DFT 周期长度;首末样本的实际时间差是 (N−1)/fs(N-1)/f_s。频率栅格间隔也不等于分开两个相近峰所需的全部条件,窗函数主瓣宽度、信噪比和工况稳定性仍会影响可分辨程度。

每转 MM 点,分析块仍为 NN 点,则

Δq=1M,Q=NM,ΔO=1Q=MN.\Delta q=\frac{1}{M},\qquad Q=\frac{N}{M},\qquad \Delta O=\frac{1}{Q}=\frac{M}{N}.

QQ 是分析块覆盖的转数,ΔO\Delta O 是阶次栅格间隔。角域奈奎斯特上限为 M/2M/2 阶;目标阶次应严格低于这个边界,并为实际抗混叠滤波器留过渡带。

例如,每转 M=128M=128 点,取 N=4096N=4096 点:覆盖 32 转,阶次间隔为 1/32=0.031251/32=0.03125,理论奈奎斯特阶次为 64。若分析到 40 阶,剩余范围可以用于滤波过渡,但是否足够仍取决于滤波器设计。

同一个 32 转分析块,在 600 rpm 稳态下持续 3.2 s,在 6000 rpm 下持续 0.32 s。因此,更长转数会改善阶次栅格分辨率,也会把更长过程平均到一起。

角域间隔 Δq=1/M\Delta q=1/M 的单位是转,时间间隔 1/fs1/f_s 的单位是秒,两者不能直接写成相等。局部稳态近似下,等角度输出对应的有效时间采样率为

fs,angle(t)=M r(t),Δt(t)≈1M r(t).f_{s,\mathrm{angle}}(t)=M\,r(t),\qquad \Delta t(t)\approx\frac{1}{M\,r(t)}.

原资料中“采样点/最大阶次”为 2.56~4 的数字属于分析仪设计举例。统一记为 M=κOmax⁡M=\kappa O_{\max} 时,κ>2\kappa>2 仅说明超出理论最低采样密度,不能自动保证任何输入都无混叠。

4. 两次采样过程都要考虑混叠

跳转到“4. 两次采样过程都要考虑混叠”

固定时间采集首先需要模拟抗混叠滤波。若最高转频为 rmax⁡r_{\max}、关注到 Omax⁡O_{\max} 阶,目标频带至少要覆盖 Omax⁡rmax⁡O_{\max}r_{\max};ADC 采样率还要为模拟滤波的过渡带留余量。

随后进行角域重采样时,原始数据中高于目标角域奈奎斯特阶次的成分也必须受到抑制。原始 ADC 没有混叠,并不代表降采样后的角域数据自动没有混叠。

例如 M=128M=128、转频 10 Hz 时,角域输出相当于局部 1280 样本/秒,640 Hz 对应 64 阶;原始采集中若存在 2 kHz 的强振动,它在该工况下是 200 阶,不能未经适当滤波就投到这个角域网格。简单线性插值只是插值方法,不等价于完整的带限重采样滤波器。

转速接近零、反转、脉冲丢失或角度不单调时,应分段处理并标注无效区间,不能直接对整段数据求逆映射。

5. 从采集到阶次曲线的设置顺序

跳转到“5. 从采集到阶次曲线的设置顺序”
  1. 定义测量对象。 选择参考轴、测点、方向和物理量;相位测量还要明确键相信号的边沿和机械零位。
  2. 设置固定采集。 依据最高转速、最高关注阶次和其他高频成分选择量程、频带及采样率;检查传感器/采集链是否饱和。
  3. 验证转速。 确认每转脉冲数、缺齿处理、阈值和迟滞,先查看 RPM—时间曲线有无跳变。
  4. 选择角域参数。 确定每转点数 MM、分析转数 QQ、窗函数和重叠率;明确幅值输出为峰值、RMS、功率还是 PSD。
  5. 计算并提取。 对角域数据做 DFT;跟踪指定阶次的幅值/相位,或按有明确归一化的阶次带积分。
  6. 回看原始信号。 阶次异常应能对应到波形、转速与工况变化;保留原始时间轴和 t(q)t(q) 关系便于复核。

6. 启停机、稳态与 Bode 展示

跳转到“6. 启停机、稳态与 Bode 展示”

启停机测试让转速穿过工作范围,可观察选定阶次激励与结构共振相遇时的幅值和相位变化。运行很慢的升速与快速扫速不能简单使用相同块长度;块内转速变化也会影响固定频谱的展宽。

稳态测量常用于趋势、相位比较和动平衡。相位比较须保持通道时序、参考脉冲定义和传感器方向一致。

旋转机械领域常把“选定阶次幅值及相位随转速变化”的成对曲线称为 Bode 图。它与控制理论的频率响应图形式相似,但若没有输入/输出比值和对应假设,不能仅凭这个图名把测得的振动幅值当成传递函数。

7. 原资料的设备能力与使用边界

跳转到“7. 原资料的设备能力与使用边界”

Crystal Instruments 原文描述 Spider/CI Order Tracking 的一种系统配置:可测最多两路转速脉冲、单模块 1~8 路动态输入,并通过多模块扩展到最多 128 路频谱或阶次分析;可记录多阶次、分数阶次、相位及固定频带随转速的变化。原理示例中的 ADC 采样率为 102.4 kHz。

这些是该来源中的产品示例,受硬件模块、软件许可和版本影响,不是阶次分析的通用要求,也不证明当前任意设备具有同样能力。

本文由既有 Crystal Instruments 中文整理稿重新组织。保留其测量链原图和设备示例,补足单位、采样域、参考轴和混叠条件,修正原稿将 Δq\Delta q 与 1/fs1/f_s 直接等同的问题。原文及图片版权归 Crystal Instruments 所有。