计算阶次分析中的阶次混叠现象.docxVIP

计算阶次分析中的阶次混叠现象 旋转机械的振动信号分析是机械故障诊断的重要部分。因为它通常是不稳定的。如果直接分析旋转机械的这些不稳定信号，就会出现严重的频率失真。为了解决这一问题,阶次分析应运而生,其核心思想是对振动信号进行等角度间隔采样。传统的阶次分析方法是基于模拟电路的硬件锁相原理实现对模拟形式的振动信号的等角度间隔采样,其最大缺点是所需设备较多,现场试验过程繁琐。随着计算机运算速度的提高,产生了计算阶次分析方法(Computed Order Tracking,COT)。 COT中存在的问题是如何确定角域重采样阶次,确保阶次谱分析结果无阶次混叠。文献提及了角域采样可能出现阶次混叠的问题,但是对于如何解决这一问题没有给出有效的方法。目前常用的角域重采样阶次确定方法,都是以信号所包含的最大阶次为已知前提,然后根据采样定理确定重采样阶次,而实际工程应用中这样默认则会出现阶次混叠的问题。文献推导了以角域重采样为目标的时域采样频率确定方法,并利用时域采样频率确定后续角域重采样阶次,在一定程度上解决了混叠的问题。但是因角域重采样是插值运算,即便所用的采样频率能满足限制信号阶次的需求,但重采样的精度会受采样频率大小的严重影响,甚至导致阶次谱图无法准确反映信号所包含的阶次。 为解决COT中容易出现的阶次混叠问题,本文深入分析了转速变化时各阶次分量对应的频域范围,提出了先确定信号的滤波范围,后确定阶次的方法,给出了滤波器截止频率选择的原则和滤波后信号角域重采 样时采样阶次的确定原则。仿真信号分析和实际实验测试结果均表明,提出的方法可成功解决阶次混叠问题,并降低角域重采样阶次。该方法无需事先知道被分析信号的最高阶次,只需给定要分析的阶次带宽即可对信号进行分析,分析结果无阶次混叠和虚假阶次。 1 滤波与阶次的关系 COT的过程如图1所示:基于等时间间隔采样得到的转速信号是脉冲信号,用脉冲时标(脉冲的上升或下降沿对应的时间点)确定参考轴转过的角度;利用脉冲时标及其转过的角度按照一定的公式进行插值计算,得到等角度间隔的重采样时标数组;基于该数组对振动信号进行插值,得到等角度间隔的振动信号。对角域重采样后的信号进行傅里叶变换,就可以得到信号的阶次谱。 COT方法的关键在于信号重采样,而造成阶次混叠的原因也出自这个环节。如果重采样阶次不够高,根据采样定理,会出现阶次混叠。本文提出采用滤波器对信号进行滤波的方法,限定信号的最高阶次。然而滤波器只能滤除信号的频率成份,要确定滤波与阶次的关系,需要关注信号不同阶次分量的频率重叠问题。 图2描述了转速变化时,信号各阶次分量对应频率的变化情况。随着转速的升高,参考轴(1阶)的转频升高,各阶次对应的频率也升高,而阶次越高频率升高得越快。当转速从ω1升高到ω3时,2阶信号分量对应的频率从f1升高到f3,3阶信号分量对应的频率从f2升高到f4。如图2,f2~f3这段频率成份包含了高转速段的2阶分量和低转速段的3阶分量,即发生了频率的重叠。此时使用滤波器对信号进行滤波,在保留信号低阶分量完整频域信息的前提下,无法把高阶分量的低转速段信息滤除。 2 滤波后信号的重采样 滤波定阶方法的基本过程是:首先根据信号需要分析的阶次确定带宽和最高转速,分析滤波频率,对信号进行滤波,滤除信号中的高频分量(含部分高阶分量);然后计算滤波后信号在给定转速范围内的最高阶次;最后依据角域重采样的采样定理确定重采样阶次。其流程如图3所示。采用该重采样阶次对信号进行重采样,通过傅里叶变换即得到阶次谱。 设某信号为随转速变化的信号,时域采样频率为fs,参考轴角速度ω变化范围已知,maxω为最高转速,minω为最低转速,需要分析的阶次带宽为Om。下面推导“滤波+定阶”方法确定重采样阶次的具体步骤。 (1) 阶次带宽om的频率分布 在转速变化过程中,需要分析的阶次带宽Om内的阶次分量在转速变化过程中的频率变化范围是: 所以,设定滤波器的截止频率为: 则可以保留需要分析的阶次带宽Om范围内的完整信息。 (2) 频率范围 时域信号经过滤波后,f≥fc的频率分量都被滤掉了。由于k阶信号对应的频率范围是Okminω/2π≤fk≤Ommaxω/2π,当阶次高到一定范围,对应的fk将始终大于fc,则该阶信号将被完全滤除。因此滤波后信号所包含的最高阶次满足: 所以滤波后信号包含的最高阶次为: (3) 确定角域中的重叠采样序列 经过滤波后,通过计算得到了信号隐含的最高阶次,此时依据采样定理确定角域重采样阶次,则可避免出现混叠。 角域重采样阶次为: 3 滤波定阶算法 为了检验上述角域重采样阶次确定方法是否合理,设定了一组仿真信号,并分析其阶次谱。 设参考轴转速信号为:y0=sin2πt2,信号的时间跨度为t=,所以转速信号的转速变化范围是:。若欲分析的信号用y表示: 如