基于信号共振稀疏分解的阶比分析及其在齿轮故障诊断中的应用

本资料包含pdf文件1个，下载需要1积分

Order domain analysis based on resonance-based sparse signaldecomposition and its application to gear fault diagnosisSUN Fun-song，YU De-.fie，CHEN Xiang-min，LI Rong(State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body，Hunan University，Changsha 410082，China)Abstract： An order domain analysis method based on resonance-based sparse signal decomposition was proposed toextracl the fault characteristic frequency from the vibration signal of a fauh gearbox with rotating speed fluctuation. Fheresonance-based spaIsr signal decomposition method decomposes signal into two parts：the high-resonancecomponentand 1ow-resonance component. The low-resonance component is an impact signal consisting of non-t)soillatorytransients with unspecified wave shape and time duration.The chirplet path pursuit algorithm was used to obtain therotating speed signal otthe gearbox.According to the rotating speed signal，the time domain impact signal of the gearboxwas resampled at conslant angle increments.Carrying out the spectral analysis on the resampled impact signal，the orderdomain analysis was accomplished and the final diagnosed results could be obtained accordingly.The proposed approach isot good anti-noise abi]ity ，and is suitable for analyzing the actual vibration signal of a gearbox with rotating speedfluctuation.An practical application example proves the validity and superiority ot the proposed method。

Key words：resonance-based sparse signal decomposition；chirplet；fault diagnosis；gears芮轮是传递动力和旋转运动的关键部件，广泛应Hj于各种机械没备之中。当齿轮出现点蚀、脱落和裂纹等局部故障时，将导致故障轮齿与其他轮齿滑动接触表面之间的润滑油棋破裂，产生冲击现象。该冲包含了齿轮局部故障的重要信息，但由于受齿轮振动信号传递路径(齿轮 -轴 -轴承 -轴承座 -测点)的影响i2，测取的冲击信号比较微弱，往往淹没在齿轮、轴等机械零部件 自身振动的周期谐波信号中，因此，如何从齿轮故障信号中彳效地提取冲击是齿轮局部损伤金项：高等学校博 学科点々项科研基金(20090161110006)；湖南势技计划资助项 日(20l2SK31 84)收稿 蝴 ：2012-06-15 修改稿收到H期：2015-09-19·作嚣 孙 云嵩 女，硕 1987 6月 通讯件 r德介 ，悼 教授，1957年 5月生诊 断的关键 。

近年来，许多学者致力于齿轮局部故障的冲击特征提取研究，取得了较多的研究成果。刘小峰等 。根据瞬态冲击成分和平稳信号在时频空间内分布特征的不同，将代表齿轮故障特征的瞬态冲击从复杂的多分量振动信号中提取出来，该方法的提取效果较好，但计算量大，且阈值不好控制；杨富舂等 根据冲击信号具有显著的非高斯特性来提取齿轮故障信号中冲击成分的频率特征，对恒转速的信号有较好的分析效果，但小适用于变转速情况；李友荣等 利用小波包将信 号分解到不同频段，然后选择有效频段进行故障信 号重构，提取齿轮故障信号中的周期冲击成分，在谐波成分与冲击成分的频率相差较大的情况下具有较好的冲击提取效果，但当两者频带重叠交叉时分析效果则不理想；第 16期 孙云嵩等：基于信号共振稀疏分解的阶比分析及其在齿轮故障诊断中的应用 89袁静等 将相邻系数降噪思想引人到平移不变多小波中，揭示出了齿轮早期裂纹的微弱故障特征信息，降噪效果好，但自适应性不佳，对谐波成分抑制性不足。李辉等 在 MOM机制的基础上，综合利用 TH-MOM阂值更新策略、自适应字典选择等措施对传统 MAC算法进行了改进，较好的提取了轴承多故障特征。

Selesnickl8 最近提出了信号共振稀疏分解方法，该方法根据信号中周期谐波成分与瞬态冲击成分品质因子(定义为中心频率与带宽的比值，用 Q表示)的不同，将信号分解为包含周期谐波成分的高共振分量与包含瞬态冲击的低共振分量。周期谐波成分为窄带信号，具有高的品质因子，而瞬态冲击成分为宽带信号，具有低的品质因子，因而，根据品质因子的差异，可实现周期谐波成分与瞬态冲击成分的分离。当齿轮处于恒定转速时，低共振分量中的冲击将按照-定的时间间隔规律重复出现，但由于受电压波动、负载等影响 ，齿轮往往处于非平稳转速下，导致低共振分量中冲击之问的时间间隔不均匀，并不满足频谱分析的平稳性要求J，因而，需预先对低共振分量进行平稳化。

工程实际中常采用阶比跟踪方法 来实现信号的平稳化，其基本思想是将时域非平稳振动信号进行等角度重采样，转化为角域平稳振动信号。常用的阶比跟踪方法有硬件阶比跟踪法、计算阶比跟踪法 、基于瞬时频率估计的阶比跟踪法 等。其中，基于瞬时频率估计的阶比跟踪法因无需安装角度编码盘、转速计等硬件，近年来得到了广泛的关注。而在基于瞬时频率估计的阶比跟踪法中，其核心问题是如何从原始振动信号中准确地提取出转速信号。线调频小波路径追踪算法 采用分段拟合的思想，能有效地提取信号中的转速信号。相对于峰值追踪算法，线调频小波路径追踪算法具有拟合精度高、抗噪能力强等优点，因而在机械故障诊断领域得到了应用Ⅲ- 。

由于信号共振稀疏分解方法能有效地实现周期谐波成分与瞬态冲击的分离，凸显齿轮故障特征；而线调频小波路径跟踪算法能从原始振动信号中精确的估计出转速信号，因此，本文将信号共振稀疏分解方法，阶比跟踪方法和线调频小波路径追踪算法相结合，提出- 种基于信号共振稀疏分解的阶比分析方法。该方法首先利用信号共振稀疏分解方法对齿轮振动信号进行分解，得到包含周期谐波成分的高共振分量和包含瞬态冲击成分的低共振分量，然后，采用线调频小波路径跟踪方法从齿轮振动信号中提取出转速信号，根据该转速信号对低共振分量进行等角度重采样，最后对角域重采样信号进行包络谱分析，根据包络谱中故障特征阶比进行齿轮故障诊断。应用实例表明，本文方法能有效地提取出变转速下的齿轮瞬态冲击成分，凸显齿轮故障特征。

1 信号共振稀疏分解方法1.1 信号的共振属性信号的共振属性用品质因子 Q定义。Q越大，信号的频率聚集性越好，具有越高的共振属性；反之，Q越小，信号的时间聚集性越好，具有越低的共振属性。

图 1表明了信号共振属性的概念。图 1左边为信号时域波形图，右边为对应的幅值谱图。图1(a)、(C)为单周期脉冲信号，品质因子 Q较小，定义为低共振信号；图 1(b)、(d)为持续多个周期的脉冲信号，品质因子Q较大，定义为高共振信号。图 1(a)所示信号与图 1(C)所示信号、图1(b)所示信号与图1(d)所示信号之问可通过时间尺度的变化互相转化，时间尺度的变化会引起脉冲信号频率发生变化，但对信号的共振属性没有影响，即具有相同的品质因子。所以高、低共振信号都可能同时包含了低频信号和高频信号。高共振信号可通过具有高Q的基函数来实现稀疏表示，而低共振信号则可通过具有低 Q的基函数来实现稀疏表示。

薹1O0150 2OO o o 4 01譬0曼05 c 00 01 0 2 0 30 4 0 5oO0 01 02 0 30 4 0采样点数N 归-化频率 H图 1 不同品质因子信号的时域波形及频谱Fig.1 Signal with different quality-factor and the specrtum传统的线性滤波方法按频带划分对信号进行分解，但当信号分量的中心频率相近且频带相互重叠时，如图1中(a)与(b)、(c)与(d)所示信号的中心频率重叠在-起，此时线性滤波方法就会失效，而信号共振稀疏分解方法从信号共振属性角度出发，综合考虑了信号中心频率与频率带宽因素，能有效分离中心频率相近且中心频率带相互重叠但具有不同品质因子的信号分量。

1.2 品质因子可调小波变换二进制小波变换作为-种恒 Q变换(其 Q值由所样函数确定)，在对分段光滑信号的稀疏表示中显示了其有效性，但由于其品质因子相对较低，频率分辨率不高，因而在对频率分辨率要求较高的信号分析中，二进制小波不适用。有理膨胀小波本质上为-种过完备的二进制小波变换，相对于二进制小波具有更高的品质因子和更高的频率分辨率。品质因子可调小波变换 与有理膨胀小波变换类似，具有完全离散、完美重构、适度完备、依赖于两通道滤波器组，并利用离散傅振 动 与 冲 击 2013年第 32卷里叶变换计算等特点。但相对于有理膨胀小波变换，品质因子可调小波变换概念简单；利用基为 2的快速傅里叶算法，计算更加高效；品质因子和冗余度更容易量化。

信号共振稀疏分解方法利用品质因子可调小波变换分别获取高 Q变换与低 Q变换的基函数库，并计算其相应的变换系数∩调品质因子小波变换通过带通滤波器组实现，其两通道滤波器组如图2所示。

(a)分解滤波器组(b)合成滤波器组2 两通道滤波器组Fig.2 Two-channel filter banks图2中 为高通尺度因子，卢 ； 为低通尺度因子， 1-旦，r表示冗余度。子带信号 。(n)的，采样频率为 ， (n)的采样频率为 ， 为原信号(n)的采样频率。

品质因子可调小波变换利用图2(a)所示两通道分解滤波器组以迭代的方式实现信号的分解， 层品质因子可调小波变换如图5所示，图3中 表示信号经过第 层变换得到的高频系数， 表示经过第 层变换得到的低频系数， l，，L。

V ，翟 F r图 3 品质凶子可调小波变换图Fig.3 Plots with adjustable quality factor1.3 高共振分量和低共振分量的分离信号共振稀疏分解方法利用形态分量分析 。。

将信号中各成分按振荡特性进行非线性分离，建立起高共振分量和低共振分量各自的最佳稀疏表示形式。

假定观测信号 可表示为两个信号 与 之和： 1 2 ， l， 2∈R (1)形态分量分析的目的即是从观测信号 中分别估计出源信号 和 。假定信号 和 可分别用基函数库(或框架)S 和 S (S。、S。具有低的相关性，本文中s 、5 分别表示为高、低品质因子可调小波的滤波器组)表示，形态分量分析的-种目标函数可表示为：，(W1，W2) -sI 1-s2 l A l 。I1 A 2 l l。 (2)式中：W ， 分别表示信号 ， 。在框架 S 、S 下的变换系数；A，、A。为正则化参数，A。、A 的取值对分解出的高共振分量与低共振分量的能量分配有影响，给定 A ，增大 A 会使A 所对应分量的能量减少；同时增大A 、A 的值，则会使残余信号能量增大。

在式(2)中，由于 f 范数不可微、且参数较多，使得式(2)的求解变得困难。信号共振稀疏分解方法利用分裂增广拉格朗日搜索算法 l，通过迭代更新变换系数 。， ，使目标函数l，最携，最终实现高共振分量和低共振分量的有效分离。

假设 目标函数 最小时，对应的高共振和低共振变换系数分别为 ， ，则求取的高共振分量和低共振分量的估计值分别表示为： S W1 ， S， (3)2 线调频小波路径追踪算法任意信号.厂(t)可展开为-组基函数的线性组合l厂(t)∑-。 h ，如果该组基函数为正交基，则可用内积计算它们的展开系数，即：a ：<-厂(t)，h >/l h， l (4)展开系数 a 的大小反映了.厂(t)与基函数的相似程度 。

线调频小波路径追踪算法采用的多尺度线调频基元函数如下D( 。 ， ，) ，，，K bg ，，exp-i(。 b t )l，(t) (5)式中：D为基元函数库；h ，，，(t)为多尺度线调频的基元函数；，为动态分析时间段，，[kN2- ～( 1)·N2- ，其中Ⅳ为分析信号的采样长度， 为分析尺度系数， 0，1，2--1；K . ，为归-化 系数，使得l h， ，，I1，n 为频率偏置系数，b 为调频率。根据采样定理， 2b t

式(5)定义的多尺度线性调频基函数在动态分析时间段内的瞬时频率为 2b t。多尺度线性调频基函数通过对信号进行逐段投影分析，计算获得每个时间分析段 内的最大投影系数和对应的线调频基元函数，该基元函数即为在时间分析段中与分析信号最为第 16期 孙云嵩等：基于信号共振稀疏分解的阶比分析及其在齿轮故障诊断中的应用 91相似的频率成分，适合在小的动态分析时间段内逐段拟合频率呈曲线变化的频率。线调频基函数的多尺度特性使得它具有了动态匹配分析信号的特性，而基函数中包含的调频率信息则使得其适合分析频率呈曲线变化的非平稳信号。

当信号与多尺度线性调频基函数越相似时，其投影系数也越大，基元函数的能量也越大，因此要求找到- 种动态分析时间段连接方法，在该连接方法下满足在整个分析时间内使连接的所有基元函数信号的总能量最大，即：max(∑ 1I K ，exp(- ( b ))1 )H ， ，2，∈， (6)兀覆盖整个分析时间段，不能重叠，其对应的最大投影系数和基元函数分别为：卢 p，'，卢如， f7、日 ， ， ， ， ，2，兀的连接方法应保证连接的基函数在整个分析时间段内的总能量最大。线调频小波路径追踪算法提出的连接算法如下：(1)初始化。以i为时间支持区序号，d( )为第 i个时间支持区之前分解信号的总能量，pre( )为连接到第 i个时间支持区的前置时间支持区序号，e(i)为第 i个时间支撑区最大投影系数对应的分解信号的能量 ，初始化时，置 d(i)0，pre(i)0。

(2)对于动态分析时间段集合 t， ∈Z中的每- 个元素， ，查找出与其相邻的所有下-个动态分析时间段集合引 ，即， 中所有元素的起始时间与 ， 相邻。如果：d(i)e( )>d(i) (8)有：d(i)e(i)d(i) pre(j)i (9)兀的连接方法可以保证在整个分析时间段内基元函数组合形成的信号与分析信号最为相似，而基元函数在动态分析时间支撑区兀，。，，2，内的瞬时频率为 2b t，t ∈Ii，对应的线性直线连接形成的频率曲线则是对分解所得信号的瞬时频率估计。

3 基于信号共振稀疏分解的阶比分析当齿轮发生故障时，振动信号主要由瞬态冲击成分和齿轮的啮合频率、转频等谐波成分组成。信号共振稀疏分解方法能够将瞬态冲击成分与谐波干扰成分进行分离，从而去除谐波干扰成分对故障诊断结果的影响。在进行非平稳齿轮振动信号分析时，线调频小波路径追踪算法可以有效的提取转速信号，将非平稳振动信号转化为角域平稳信号，从而满足信号分析对信号平稳性的要求，同样其提取的转速信号也能够将时域非平稳的瞬态冲击成分转化为角域平稳的瞬态冲击成分，以便利于故障诊断。因此本文首先用信号共振稀疏分解方法提取瞬态冲击成分，然后利用线调频小波路径追踪算法获得齿轮箱转速信号，根据转速信号对提取的冲击信号进行等角度重采样，再对角域冲击信号进行阶比分析，完成基于信号共振稀疏分解的阶比分析过程。主要分析步骤如下：(1)对齿轮故障振动信号进行共振稀疏分解，将信号分解为高共振部分和低共振部分；(2)利用线调频小波路径追踪算法对原信号进行瞬时频率估计 ，得到瞬时啮合频率，根据瞬时啮合频率求出瞬时转频；(3)用三阶多项式对瞬时转频进行拟合，有：l厂( )at bt ctd (10)a，b，C，d为三阶多项式系数，由最小二乘法拟合得到 ；(4)确定阶次跟踪的最大分析阶次 D ；(5)计算重采样后数据的长度 Ⅳ和等角度重采样的键相时标 ，根据求出的键相时标 ，利用 Lagrange线性插值公式对低共振分量进行插值，求出低共振分量在角域里对应于键相时标 的幅值 ( )。其中Ⅳ根据式(11)得出，键相时标 根据式(12)得到，幅值 (T )根据式(13)得到。式(11)～(13)的表达式如下 ：r2 f.厂(t)dtNAO (11)式中： 为时域采样的总时间 t)为频率拟合函数。

导 ÷ d 詈 十 十号焉d71o n 1，2，式中： 为时域采样开. ～始时间(tf)tilt ≤ To≤ t (6)对 ( )(n1，2，3，，Ⅳ)进行包络阶次分析，根据得到的低共振分量阶比谱进行故障诊断。

4 应用实例4.1 断齿及正常齿轮分析将齿轮箱故障试验台上的某-主动齿轮的-个齿进行破坏，即将某-主动齿轮切割-个齿，模拟齿轮断齿故障，在变转速下采集-组断齿齿轮振动信号。输入轴齿轮齿数为 55，输出轴齿轮齿数为 75，采样频率92 振 动 与 冲 击 2013年第 32卷(信号 振稀疏分解提取冲击 ) I --- 二二 二二 - J l瞬念冲， 成分l I矧4 基于信 共振稀疏分解的阶比分析原理框图Fig.4 Schematic of the order domain analysis methodbased on resonance-based sparse signal decomposition为4 096 Hz，采样时长为2 s，正常齿轮和断齿齿轮参数相同。 5为变转速条件下采集到的断齿齿轮振动信弓删 域波形，从中呵以看出断齿信号存在问隔不均匀的冲击。

利用共振稀疏分解方法对断齿振动信号进行分解，取高品质因子 Q 为5，低品质因子 Q 为 1，结果如图6所示，图6(a)为高共振分量，即信号谐波分量，图6(b)为低共振分量，即提取的瞬态冲击，冲击的时序间隔并不均匀，具有非平稳性，很难直接看出有效信息。

f 图 5 断齿齿轮实测信号Fig.5 Measured broken gear vibration sigual利用线调频小波路径追踪方法进行转频估计，得到转速信号如图7中所示。利用拟合转频将提取到的故障冲击信号角域平稳化，对平稳化后的冲击信号部分进行包络阶次分析，得到阶次谱结果如图 8所示。

当齿轮发生故障时，齿轮箱振动信号主要表现为以故障齿轮所在轴转频及其倍频为调制频率的啮合频率调制，且振动能量和包络能量较大。图 8仅可见-个明显的峰值谱线，峰值 1点坐标为(2，0.018 47)，说明调制频率为转频的2倍频，可见该齿轮发生故障，且无干扰分量。分析结果验证了本文所提出的基于信号共振稀疏分解的阶比分析方法的有效性。

为了进行对 比，用线调频小波路径追踪算法对该断齿齿轮振动信号进行阶比分析，得到结果如图9。从图中可以看到有两个明显峰值谱线，峰值点 1、2坐标分别是(0.757 9，0.061 05)、(2.021，0.131 2)，据此分析得齿轮调制频率为0.757 9倍频和2.021倍频，二倍频为齿轮故障引起的调制频率，但0.757 9倍频不是转频及其整数倍频。由于齿轮故障振动信号中存在-些谐波干扰信号，在解调时会作为干扰分量出现在阶比谱中，可以诊断0.757 9倍频为信号中的干扰成分造成的。在基于信号共振稀疏分解的阶比分析方法中，由于将瞬态冲击信号与谐波干扰信号分离，避免了干扰成分的影响，因此更有利于齿轮故障分析。

同样可以得到正常齿轮振动信号经过信号共振稀疏分解后的阶比谱如图10所示。峰值 1、2、3、4点坐标分别为 (1.022，0.002 485)、(2.022，0.002 737)、(4.044，0.002 439)、(5.044，0.002 057)，说明其调制频率分别为-倍、二倍、三倍 、四倍等转频。但与图 8相比，正常齿轮阶比谱中峰值明显相对较校由于齿轮加工误差会引起的齿形误差等问题，导致齿轮在运转过程会引起-定程度的冲击，在振动信号分析过程亦会出现转频调制现象，但其能量小，不会对故障齿轮信号分析造成影响。

I- 山山。I .山 4 j-。jIlfL IL山-I ”0 30 0 35 0 40 0 45 0 50 0 55 0 60t/s(a)高共振分量塞。 州卅 馨jJ l l。 。

24 524b4基0 0 2 0 4 0 6 0 8t/s图7 线调频小波路径追踪算法拟合的断齿齿轮转频曲线Fig.7 The rotation frequency of the broken gearobtained by the chirplet path pursuit algorithm，22 l1粤阶比图8 断齿齿轮低共振分量阶比谱Fig.8 The order spectrum of the low-resonancecomponent of the broken gear害第 16期 孙云嵩等：基于信号共振稀疏分解的阶比分析及其在齿轮故障诊断中的应用 93阶比n图9 断齿齿轮振动信号阶比谱Fig.9 The order spectrum of the broken gear2莹 1量罂阶比图 lO 正常齿轮振动信号低共振分量阶比谱Fig.10 The order spectrum of the low-resonancecomponent of the normal gear4.2 轮齿裂纹故障分析利用线切割技术在轮齿根处加工出宽 0.1 mm、深3 mm的槽，并在无负载条件下模拟齿轮齿根裂纹。输入轴齿轮齿数为37，输出轴齿数为 37，在变转速条件下利用加速度信号传感器拾取振动信号，采样频率为2 048 Hz，采样点数为4 096，在变转速情况下采集到的齿轮裂纹振动信号如图 1 1所示。

利用共振稀疏分解方法对裂纹信号进行分解 ，取高品质因子 Q 为5，低品质因子 Q 为 1，得到高共振分量如图 12(a)所示，低共振分量如图 12(b)所示。低共振分量冲击的时序间隔不均匀，具有非平稳性，难以直接看出有效信息。

图 l3为利用线调频小波路径追踪方法得到的转频曲线。利用拟合转频曲线对低共振分量等角度重采样，再对角域平稳信号进行包络阶次分析，得到阶次谱结果如图 14所示。图 l4仅可见-个明显的峰值谱线 ，峰值 1点坐标为(2，0.040 53)，说明调制频率为转频的2倍频，可见该齿轮发生故障，且无干扰分量。分析结果再次验证了本文方法的有效性。

为了进行对 比，利用拟合转频曲线直接对裂纹齿轮振动信号等角度重采样，再对角域平稳信号进行包络阶次分析，得到结果如图 15所示。从图中可以看到有两 个 明显 峰值 谱 线，峰值 点 1、2坐标 分别 是(0.537 9，0.050 76)、(2.005，0.097 75)，据此分析得齿轮调制频率为0.537 9倍频和2.005倍频，二倍频为齿轮故障引起的调制频率，但0。537 9倍频不是转频及其整数倍频，同样可以诊断0.757 9倍频为信号中的干扰成分造成的。从图 14与图 15的对比可以看出，本文方法能够有效分离瞬态冲击信号与谐波干扰信号，去除干扰成分的影响，更有利于齿轮故障分析。

篓 -1g罂t/s(a)高共振分量U U 肪 q9 n髓 1 1 ∞ t/s(b)低共振分量图l2 裂纹齿轮振动信号共振稀疏分解结果Fig.12 The resonance-based sparse signal decompositionresults of the cracked gear vibration signal图 13 线调频小波路径追踪算法拟合的裂纹齿轮转频曲线Fig.13 The rotation frequency of the cracked gearobtained by the chirplet path pursuit algorithm苔- ×j型骥图 l4 裂纹齿轮振动信号低共振分量阶次谱Fig.14 The order spectrum of the low-resonancecomponent of the cracked gear毫耋图 15 裂纹齿轮振动信号阶比谱Fig.15 The order spectrum of the cracked gear5 0 5 O ∞.暑--×≈ 罂。.弓魈 -～振 动 冲 击 2013年第 32卷5 结 论本文将信号共振稀疏分解方法与线调频小波路径追踪算法相结合，提出基于信号共振稀疏分解的阶比分析方法，并将其应用于非平稳转速状态下的齿轮箱故障诊断♂合应用实例分析，得到如下结论(1)信号共振稀疏分解方法能有效分离中心频率相近且中心频率带相互重叠但具有不同品质因子的信号分量～共振分解方法应用于冲击故障特征提取，能够将瞬态冲击与谐波分量分离，去除干扰分量，提高瞬态冲击成分提取能力。

(2)将基于信号共振稀疏分解的阶比分析方法应用于实际齿轮振动信号分析，结果表明该方法能够将谐波分量与冲击分离，从而有效克服了干扰信号对齿轮故障诊断的影响，与线调频小波路径追踪算法结合，保证了阶比分析结果的准确性，为非平稳转速下齿轮箱故障的特征频率的提娶故障识别和诊断提供了-种新的分析方法。