变分模态分解和改进的自适应共振技术在轴承故障特征提取中的应用

李华 伍星 刘韬
摘要: 针对滚动轴承早期故障特征提取困难的问题,提出了基于变分模态分解(Variational Mode Decomposition, VMD)和改进的自适应共振技术的滚动轴承故障特征提取方法。针对轴承故障信号所在频带难以选择的问题,提出了基于改进的自适应共振技术(Improved Adaptive Resonance Technology, IART)的IMF选取方法。首先,确定模态数,提出了峭度最大值的模态数确定方法;然后,对原始振动信号进行VMD分解,获得既定数目的本征模态分量(Intrinsic Mode Function, IMF);其次,利用IART选取包含丰富故障信息的IMF分量;最后,(如有需要)对选取的IMF分量进行基于IART的带通滤波,并进行包络解调分析提取故障特征频率。将该方法应用到轴承仿真数据和实际数据中,能够实现轴承故障特征的精确诊断,证明了该方法的有效性。
关键词: 故障诊断; 滚动轴承; 变分模态分解; 峭度; 改进的自适应共振
中圖分类号:TH165+.3; TN911.7文献标志码: A文章编号: 1004-4523(2018)04-0718-09
DOI:10.16385/j.cnki.issn.1004-4523.2018.04.020
引言
对轴承的状态监测与故障诊断一直是机械设备故障诊断的热点[1]。当滚动轴承发生故障时,其振动信号包含了大量的运行状态信息,表现为非平稳性和多分量性的调制信号,特别在故障早期,由于调制源弱,早期故障信号微弱,并且受周围设备、环境的噪声干扰,导致故障特征频率难以提取、识别[2]。
故障诊断的关键是从原始信号中提取故障特征信号(故障特征频率)。Huang等[3]提出了经验模式分解,此方法具有自适应分解特性,对非平稳和非线性信号的处理具有较高的效率。因此,基于EMD的故障诊断方法层出不穷。但EMD存在模态混叠、端点效应、受采样频率影响较大等不足[4-5]。基于此,Dragomiretskiy等[6]提出了一种自适应信号处理新方法——变分模态分解(VMD),该方法通过迭代搜寻变分模型最优解来确定每个分量的频率中心及带宽,相比于EMD,VMD将信号分解转化为非递归、变分模态分解模态并且具有坚实的理论基础,有效地弥补了EMD和EEMD等的不足。很多学者已将VMD应用于轴承故障诊断研究。Mohanty等[7]将VMD引入轴承故障诊断,并取得了比EMD更好的效果。刘长良等[8]提出了基于变分模态分解的故障特征提取方法,利用观察中心频率的方法确定模态K,并采用模糊C均值聚类进行故障识别。赵洪山等[9]提出了基于变分模态分解和奇异值分解的故障特征提取方法,利用峭度指标选取敏感IMF分量,对选取的IMF进行奇异值分解重构,取得良好效果。刘尚坤等[10]以互信息为判据对原VMD方法进行了迭代停止条件的改进,并结合Teager能量算子的优点,提出了Teager-VMD时频分析新方法。针对故障滚动轴承振动信号的非平稳和非高斯特性,Te Han等[11]提出了一种基于变分模态分解(VMD)和自回归(AR)模型参数的故障特征提取方法。并且在滚动轴承故障诊断领域提出了一种新型随机森林分类器用于模式识别。夏均忠等[12]提出了一种用最大相关峭度解卷积(Maximum Correlated Kurtosis Deconvolution, MCKD)和变分模态分解相结合的方法提取滚动轴承故障特征。Zipeng Li等[13]针对背景噪声严重影响模式个数的设置,提出了一种基于相关性分析的面向独立的VMD方法,用于自适应提取轮对轴承的弱故障和复合故障特征。
当滚动轴承发生故障时,由于振动信号表现出的幅值调制特性,通过包络分析可以得到轴承的故障特征频率[14]。鉴于VMD的模态数需事先确定的问题及其敏感IMF分量的选取问题,本文提出VMD和改进的自适应共振技术相结合的诊断方法。
1基础理论
1.1变分模态分解(VMD)变分模态分解(VMD)[6]是一种基于维纳滤波的自适应信号处理新方法。在VMD算法中,每一个IMF均被称为AM—FM信号,即ukt=Aktcoskt(1)式中ωkt=′kt=dktdt,Akt为信号的瞬时幅值,ωkt为瞬时频率。在t-δ,t+δδ=2π'kt内,可将ukt视为幅值为Akt,频率为ωkt的谐波信号。
第4期李华,等:变分模态分解和改进的自适应共振技术在轴承故障特征提取中的应用振 动 工 程 学 报第31卷VMD的求解过程是一个变分问题的求解过程。假设每个模态是具有中心频率的有限带宽,中心频率和带宽在分解过程中不断更新,则,变分问题可表示为寻求k个模态函数ukt,并使得所有模态函数的估计带宽和最小,模态之和为输入信号f。具体步骤如下:
1.2峭度
峭度(Kurtosis)反映振动信号分布特性的数值统计量,是描述波形尖峰度的无量纲参数[15]。
其数学描述为K=Ex-μ4σ4(6)式中x为所分析的振动信号,μ为信号x的均值,σ为信号x的标准差。
峭度对冲击信号敏感,适合轴承早期故障诊断。滚动轴承无故障时,峭度指标值约等于3。当滚动轴承发生局部故障时,峭度值会增大。
1.3改进的自适应共振技术(IART)
由于传统共振解调技术的缺点,共振频率和带通滤波器的中心频率都是固定的,王平等[16]提出了一种自适应谐振解调技术。在此基础上,本文提出了一种基于功率谱的共振频率选取方法:直接对信号进行功率谱分析,依靠功率谱分析来自动识别原始振动信号和各IMF的共振频率(功率谱反映了能量的集中情况)。并在此基础上确定带宽,设计带通滤波器。由于实际情况下共振频率一般在故障频率的10倍以上,因此设定一个阈值,当振幅峰值出现在故障频率的10倍以下时,将被放弃。
选取共振频率存在如下两种情况:1)当满足阈值条件的功率谱中只有一个幅度峰值时,可以直接选择为中心频率,即共振频率;2)当满足阈值条件的功率谱中存在多条谱线时,选择幅值最大的功率谱线作为中心频率,即共振频率。
然后确定带宽,这里参考文献[12]的设定方法,将带通滤波器带宽设定为3倍故障频率。
2基于VMD和IART的故障特征提取针对滚动轴承早期故障诊断中故障特征微弱难以提取的问题,本文提出VMD和IART相结合的方法。然后,对选取的IMF进行基于IART的带通滤波,进一步提高信噪比,并对滤波后的信号进行包络解调分析,提取故障特征频率。
2.1基于峭度最大值的模态优化
因为峭度指标是故障的一种表征,当存在某个模态的峭度值具有最大值时,可认为此时具有良好的分解效果(在给定的搜索范围内)。基于此,本文提出了基于峭度最大值优化模态数K的方法。具体步骤及流程在2.3节中给出。
2.2基于IART选取敏感IMF分量
轴承振动信号经过VMD分解为一系列IMF分量,若要进行故障特征提取,最好的方法是选取包含故障信息最多的IMF分量进行分析,为此,本文提出了基于IART的敏感IMF选取方法。
对信号进行功率谱分析,选取共振频率,并以此为中心频率设计带通滤波器,将此频带范围作为特征频带区间,利用区间从属关系选取有效IMF分量。用双值特征区间系数代替一般的单值特征参考系数,有效地避免了单值参考系数的不易选取的问题。
首先,对原信号和VMD分解获得的若干IMF分量分别求取功率谱,选取共振频率;然后,以此共振频率作为中心频率设计带通滤波器,并将此滤波频带作为其特征频带区间;其次,比较各个IMF分量的特征频带区间与原信号特征频带区间的从属关系,若满足,则选取为敏感IMF,否则舍弃;如果存在多个敏感IMF,则对其进行叠加、重构。
2.3算法步骤及流程
基于VMD和IART的故障特征提取方法的具体步骤及诊断流程图如图1所示。
(1)获取原始振动信号,初始化模态数K=2,惩罚因子α和带宽τ使用默认值:α=2000,τ=0;
(2)对轴承原始振动信号进行VMD分解,计算各个模态的峭度值,并通过比较获得此模态数下的峭度最大值,然后K=K+1继续以上分析,直到取到K=16为止;比较在各个模态数下所取得的峭度最大值的大小,将最大的峭度值所对应的模态数K判定为最优模态数K。为了确定模态数K的搜索范围和步长,本文借鉴文献[17]给出的模态数K的讨论范围,由于K过大,造成效率低,计算负担重;K过小,易引入噪声,所以搜索范围设为[2,16]。
(3)利用IART选取敏感IMF分量。计算原始振动信号与各个IMF分量的功率谱,获得共振频率,并以此共振频率作为中心频率设计带通滤波器。选取滤波频带从属于原始信号滤波频带的IMF分量为敏感IMF进行信号重构;
(4)对重构信号进行基于IART的带通滤波,此处的带通滤波器由步骤(3)获得,进一步提高SNR;
(5)对步骤(4)中滤波后的信号进行包络解调分析,判断故障发生部位。
3信号仿真分析
对原信号进行11个模态的VMD分解,获得11个IMF分量,由于篇幅限制,仅给出前6个IMF分量的时域波形,如图4所示。由于时域图差别不大,无法直观地选择敏感IMF分量进行分析。因此,需要合适的方法选取敏感IMF分量。
3.2有效IMF选取及带通滤波
利用本文提出的IART选取有效的IMF分量,对原始信号和各IMF分量求取功率谱,按照第1.3节获得共振频率,以此作为中心频率设计带通滤波器,获得频带区间。同样地,由于篇幅限制,只给出了原信号和IMF5-IMF7的功率谱图,如图5所示。
由图5可知,原信号的共振频率选取为3000 Hz。以此作为中心频率设计带通滤波器,根据文献[12],取滤波带宽为3倍的故障特征频率,即Bw=3fin=237 Hz,进一步获得滤波频带为[2882,3119],将其作为原信号的特征频带区间。从各IMF的功率谱可知,IMF6的共振频率在3000 Hz。同样地,设计带通滤波器,Bw=3fin=237 Hz,其特征频带区间为[2882,3119]。可知,IMF6的滤波频带从属于原信号的滤波频带,故选取IMF6为敏感IMF分量进行分析。原信号和IMF5-IMF7的特征频带区间如表1所示。
对IMF6进行包络解调分析,如图6所示。从图中可以清晰地提取的故障特征频率(78.52 Hz)、转频(28.13 Hz)以及二倍频成分,均与理论值很接近,由此可判别故障类型为内圈故障。这充分说明了本文提出的方法的有效性。而且,此时已经几乎沒有噪声的影响,只有少许微弱噪声毛刺,所以最后一步的带通滤波可以舍弃(此处仍给出了带通滤波的包络谱,如图7所示。由图可知,残余的微弱噪声被有效的去除,效果良好)。
4实验验证及分析
为了证明本文提出方法的有效性及优越性,将其用于实际轴承数据分析。本文采用来源于美国西储大学的轴承数据[19],轴承的型号为6205RS JEM SKF,采样频率fs=12000 Hz,试验数据选择转速为1750 r/min,负载为2 hp(1492 W),故障尺寸为0.021″,驱动端轴承在滚动体故障状态下的数据。
4.1VMD分析
如图8(a)和(b)分别是轴承滚动体故障数据的原始信号时域图和包络谱。时域图中由于噪声的严重影响,无法提取冲击特征;包络谱中,虽然能够提取故障特征频率,但受到很严重的宽频带噪声影响,如图中 2000~3000 Hz范围内的波峰尤为突出(尤其红色椭圆框内),无法判别。因此,需要对信号进行预处理,提高信噪比。
为此,对原始振动信号进行VMD分解,初始化模态数K=2,设定K的搜索范围为[2,16](由文献[13]可知,K过大,造成效率低,计算负担重;K过小,易引入噪声),利用本文提出的峭度最大值原则确定最优模态K。经过分析,当K=6时,取得峭度最大值为3.4620,所以取最优模态数为K=6。模态数与峭度最大值关系图如图9所示。
对原信号进行6个模态的VMD分解,获得6个IMF分量,各个IMF分量的时域波形如图9所示。同样地,由于时域图各个IMF分量差别不大,无法直观地选择敏感IMF分量进行分析。因此,利用本文提出的方法选取敏感IMF分量。
4.2有效IMF选取及带通滤波
利用本文提出的基于IART的方法选取有效IMF分量。由于篇幅限制,只给出了原信号和IMF3-IMF6的功率谱图,如图11所示。
如图11(a)所示,根据1.3节所述,可取原始振动信号的共振频率为3345 Hz。以此作为中心频率设计带通滤波器,按文献[12]取滤波带宽为3倍的故障特征频率,即Bw=3fb=412 Hz,进一步获得滤波频带为[3139,3551],并将其作为原信号的特征频带区间。
然后,获取各个IMF分量的共振频率,并设计带通滤波器,获得特征频带区间。经分析,可知IMF5的共振频率为3345 Hz,与原始振动信号的一致。同样地,设计带通滤波器,Bw=3fin=412 Hz,其特征频带区间为[3139,3551]。可知,IMF5的滤波频带从属于原信号的滤波频带,故选取IMF5为敏感分量进行分析。原信号和IMF4-IMF6的特征频带区间如表3所示。
为了证明本文提出的IMF选取方法的优越性,这里利用峭度最大值原则对有效IMF进行选取。对各个IMF分量进行峭度计算,经比较,IMF4具有峭度最大值,因此,选取IMF4为有效分量。对其进行包络解调分析。如图13所示,从图中虽然能够提取转频,但故障特征频率几乎无法提取,而且存在很多未知的波峰干扰。显然,本文图12具有更好的效果。
从图12可看到,虽然已经能够比较清晰地提取故障特征频率,但存在一定的噪声影响,比如在310 Hz附近的波峰。这里,笔者提出利用选取IMF时设计的带通滤波器对选取的IMF进行带通滤波,进一步提高信号的信噪比。如图14所示,经过滤波后,有效地去除了干扰噪声(310 Hz附近的波峰被去除)。证明了基于VMD和IART的故障诊断方法的有效性。
为了验证基于峭度最大值优化模态数K的优势,取K=11,利用本文所述的IART方法选取的敏感IMF分量(经分析,只有IMF7的特征频带与原信号的特征频带满足从属关系,故选取的敏感IMF分量为IMF7)的包络谱如图15所示。容易看出,雖然可以提取出转频,但故障频率几乎淹没在噪声中,由图12所呈现的效果明显的更为理想。由此表明,相比于人为主观决策,基于峭度最大值原则的模态数确定方法提供了理论依据,且更为可靠。
4.3与EEMD对比
采用EEMD算法对上述滚动体故障信号进行分解,得到14个IMF分量,取前3个IMF分量求取频谱,如图16所示。由图可以看到EEMD分解存在频率混叠,不同分量之间的频带差不明显,且互相影响。
4.4基于VMD和IART的弱故障分析
为了能说明本文方法的适用性以及对微弱故障的处理能力,在这里分析了文[19]中轴承的型号为6205RS JEM SKF,故障尺寸为0.007″,转速为1,730 r/min,负载为3 hp(2238 W),采样频率fs=12000 Hz的滚动体故障数据。计算可知,转频fr=28.7 Hz,滚动体故障特征频率fb=134.89 Hz。
4.4.1VMD分析
如图19(a)和(b)分别是轴承滚动体弱故障数据的原始信号时域图和包络谱。时域图中由于噪声的严重影响,无法提取冲击特征;包络谱中由于受到很严重的噪声影响,只能提取到转频,而无法提取故障特征频率。故利用本文方法度其进行处理,提高信噪比。
利用峭度最大值原则确定最优模态K。经分析,当K=11时,具有峭度最大值(3.762),所以最优模态为K=11。对原信号进行11个模态的VMD分解,获得11个IMF分量。
4.4.2有效IMF选取
根据IART选取有效IMF分量,可知,IMF8为有效IMF分量,故对其进行包络解调分析,如图20(a)所示。同样地,利用峭度最大值原则选取IMF3为敏感分量,其包络谱如图20(b)。由图可知,图(a)可以清晰地提取故障特征频率(131.3 Hz)和转频(29.3 Hz),均与理论计算值很接近,而图20(b)只能提取2倍转频,无法提取故障特征频率。由此可知,本文提出的选取方法具有明显优势。证明了本文方法可以有效的应用于微弱特征的提取。由图20(a)可知,此时已经几乎没有噪声影响,所以无需进行二次滤波。
同样地,利用EEMD对信号进行分析,可得14个IMF分量。基于IART选取敏感分量, IMF1和IMF2满足从属关系,将其叠加重构,包络谱如图21所示。由图可知,包络谱仅能提取转频,无法提取故障特征频率。且存在很严重的噪声干扰。故也说明了本文方法具有明显优势。
5结论
1)首先,本文针对模态数K,提出了峭度最大值原则的优化方法,能够有效地确定最优模态数K。并且比人为主观决策更有理论依据更可靠。
2)针对轴承故障信号所在频带难以选择的问题,提出了基于改进的自适应共振技术(IART)的IMF选取方法。在对轴承滚动体故障的两个实例分析中,均取得了比峭度最大值原则选取IMF更好的效果。
3)针对VMD分解后,若几乎不存在噪声影响,则完成诊断过程;若存在的残余噪声对故障特征提取的影响,提出了将VMD与IART相结合的故障特征提取方法。并在对实际轴承故障的分析中证明了其有效性。
4)VMD还存在缺陷,K,α等需事先确定,而且参数的搜索范围建立在前人讨论的基础上,缺乏理论依据,有待进一步完善。而自适应共振的中心频率及带宽的选取也需进一步研究。
参考文献:
[1]Wang H Q, Hou W, Tang G, et al. Fault Detection Enhancement in Rolling Element Bearings via Peak-Based Multiscale Decomposition and Envelope Demodulation[J]. Mathematical Problems in Engineering, 2014, 2014(1):135—142.
[2]丁康, 黄志东, 林慧斌. 一种谱峭度和Morlet 小波的滚动轴承微弱故障诊断方法[J]. 振动工程学报, 2014, 27(1):128—135.
DING Kang, HUANG Zhidong, LIN Huibin. A weak fault diagnosis method for rolling element bearings based on morlet wavelet and spectral kurtosis [J]. Journal of Vibration Engineering, 2014, 27(1):128—135.
[3]Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for non-linear and non-stationary time series analysis[C]. Proceedings of the Royal Society Lond,1998, 454:903—995.
[4]Wu Zhaohua, Huang N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis, 2009, 1(1): 1—41.
[5]Rilling G, Flandrin P. On the influence of sampling on the empirical Mode Decomposition[C]. IEEE International Conference on Acoustics, Speech and Signal Processing. ICASSP 2006 Proceedings,IEEE, 2006:444—447.
[6]Dragomiretskiy K, Zosso D. Variational Mode Decomposition[J]. IEEE Transactions on Signal Processing, 2014, 62(3):531—544.
[7]Mohanty, Gupta K K, Raju K S. Bearing fault analysis using variational mode decomposition[C].International Conference on Industrial and Information Systems. IEEE, 2014:1—6.
[8]劉长良, 武英杰, 甄成刚. 基于变分模态分解和模糊C均值聚类的滚动轴承故障诊断[J]. 中国电机工程学报, 2015, 35(13):3358—3365.
LIU Changliang, WU Yingjie, ZHENG Chenggang. Rolling bearing fault diagnosis based on variational mode decomposition and fuzzy C means clustering [J]. Proceeding of the CSEE, 2015, 35 (13): 3358—3365.
[9]赵洪山, 郭双伟, 高夺. 基于奇异值分解和变分模态分解的轴承故障特征提取[J]. 振动与冲击, 2016, 35(22):183—188.
ZHAO Hongshan, GUO Shuangwei, GAO Duo. Fault feature extraction of bearing faults based on singular value decomposition and variational modal decomposition [J]. Journal of Vibration and Shock, 2016, 35 (22): 183-188.
[10]刘尚坤, 唐贵基, 王晓龙. 基于改进变分模态分解的旋转机械故障时频分析方法[J]. 振动工程学报, 2016, 29(6):1119—1126.
LIUShangkun, TANG Guiji, WANG Xiaolong. Time frequency analysis method for rotary mechanical fault based on improved variational mode decomposition[J]. Journal of Vibration Engineering, 2016, 29(6):1119—1126.
[11]Han T, Jiang D. Rolling bearing fault diagnostic method based on VMD-AR model and random forest classifier[J]. Shock and Vibration, 2016, 2016:1—11.
[12]夏均忠, 趙磊, 白云川,等. 基于MCKD和VMD的滚动轴承微弱故障特征提取[J]. 振动与冲击, 2017, 36(20):78—83.
XIAJunzhong, ZHAO Lei, BAI Yunchuan, et al. Feature extraction for rolling element bearing weak fault based on MCKD and VMD[J]. Journal of Vibration and Shock, 2017, 36(20):78—83.
[13]Li Z, Chen J, Zi Y, et al. Independence-oriented VMD to identify fault feature for wheel set bearing fault diagnosis of high speed locomotive[J]. Mechanical Systems & Signal Processing, 2017, 85:512—529.
[14]Yang Y, Yu D, Cheng J. A fault diagnosis approach for roller bearing based on IMF envelope spectrum and SVM[J]. Measurement, 2007, 40(9-10):943-950.
[15]胡爱军, 马万里, 唐贵基. 基于集成经验模态分解和峭度准则的滚动轴承故障特征提取方法[J]. 中国电机工程学报, 2012, 32(11):106-111.
Hu Aijun, Ma Wanli, Tang Guiji. Rolling bearing fault feature extraction method based on ensemble empirical mode decomposition and kurtosis criterion[J]. Proceedings of the CSEE, 2012, 32(11):106—111.
[16]王平, 廖明夫. 滚动轴承故障诊断的自适应共振解调技术[J]. 航空动力学报, 2005, 20(4):606—612.
WANG Ping, LIAO Ming-fu. Adaptive Demodulated Resonance Technique for the Rolling Bearing Fault Diagnosis[J]. Journal of Aerospace Power, 2005, 20(4): 606—612.
[17]唐贵基, 王晓龙. 变分模态分解方法及其在滚动轴承早期故障诊断中的应用[J]. 振动工程学报, 2016, 29(4):638—648.
TANG Guiji, WANG Xiaolong. Variational mode decomposition method and its application on incipient fault diagnosis of rolling bearing [J]. Journal of Vibration Engineering, 2016, 29 (4): 638—648.
[18]Liu T, Chen J, Dong G, et al. The fault detection and diagnosis in rolling element bearings using frequency band entropy[J]. Proceedings of Institution of Mechanical Engineers Part C Journal of Mechanical Engineering Science, 2013, 227(1):87—99.
[19]http://www.cwru.edu/laboratory/bearing/welcome_over view. htm.
Abstract: According to the difficult problem that the fault features extraction of rolling bearings in early failure duration, an incipient fault diagnosis method for rolling bearings based on the variational mode decomposition (VMD) and the improved adaptive resonance technology (IART) is proposed. According to the problem that the frequency band of the bearing fault signal is difficult to choose, an intrinsic mode function (IMF) selection method based on IART is proposed. Firstly, the mode number determination method based on kurtosis maximum value is proposed here to determine the mode number. Then, the original vibration signal is decomposed by VMD to obtain IMF, and the IMF component with abundant fault information is selected by IART. Finally, the selected IMF component is subjected to band-pass filtering based on IART if necessary, and the fault characteristic frequency is extracted by envelope demodulation analysis. The method can be applied to the bearing simulation data and the actual data, which can realize the accurate diagnosis of bearing fault characteristics and prove the effectiveness of the method.
Key words: fault diagnosis; rolling bearing; variational mode decomposition; kurtosis; improved adaptive resonance technology