基于经验小波变换的机械故障诊断方法

基于经验小波变换的机械故障诊断方法农机12013120306027靳伟摘要：经验小波变换(EWT)是一种新的自适应信号分解方法，该方法继承了EMD和小波分析方法的各自优点，通过提取频域极大值点自适应地分割傅里叶频谱以分离不同的模态，然后在频域自适应地构造带通滤波器组从而构造正交小波函数，以提取具有紧支撑傅立叶频谱的调幅一调频(AM．FM)成分。本文将该方法引用到机械故障诊断中，提出了一种基于经验小波变换的机械故障诊断方法，并与EMD方法进行了对比分析。仿真结果表明，经验小波变换法明显优于EMD方法，能有效地分解出信号的固有模态。与EMD相比较，该方法具有分解的模态少，不存在虚假的模态，计算量小，且在理论上具有易理解性等特点。最后将该方法应用到转子碰磨故障诊断中，实验结果进一步验证了该方法的有效性，能够有效地揭示出碰磨故障数据的频率结构，区分碰磨故障的严重程度。关键词：经验小波变换；固有模态；故障诊断：自适应信号分解；经验模态分解MechanicalfaultdiagnosismethodbasedonempiricalwaveletTransformAbstract:Abstract：Empiricalwavelettransform(EWT)isanewselfadaptivesignaldecompositionmmethod．ThismethodinheritstheadvantagesofEMDandwavelettransform，adaptivelysegmentstheFourierspectrumbyextractingthemaximapointinthefrequencydomaintoseparatethedifferentmodes，andthenconstructsadaptiveband—passfiltersinthefrequencydomainsoastoconstructorthogonalwaveletfunctionsandextractAM—FMcomponentsthathaveacompactsupportFourierspectrum．Here，theEWTisintroducedinthemechanicalfaultdiagnosis．andanewmechanicalfaultdiagnosismethodbasedonEWTisproposed．TheEWTmethodiscomparedwiththetraditionalEMDmethod．ThesimulationresultsshowthattheEWTmethodisobviouslysuperiortotheEMDmethod．Theproposedmethodcaneffectivelydecomposetheintrinsicmodesofthesigna1．ComparedwiththeEMDmethod，thismethodhassomedistinctadvantages，suchaslessdecomposedmodes，novirtualmodes，lesscalculation，easytobeunderstoodintheoryandetc．Finally,theproposedmethodwassuccessfullyappliedtotherub—impactfaultdiagnosisofarotorsystem．Theexperimentresultsshowthattheproposedmethodisefective，andcaneffectivelyrevealthefrequencystructureinrubbingfaultanddiscerntheseverityofrub—impactfault．Keywords：EWT；intrinsicmode；faultdiagnosis；adaptivesignaldecomposition；EMD1引言经验模态分解(empiricalmodedecomposition，EMD[1][是由美国国家宇航局美籍华人黄锷(N．E．Huang)等于1998年创造性地提出的一种新型自适应信号时频分析方法。该方法克服了传统方法中用无意义的谐波分量来表示非平稳、非线性信号的缺陷，并具有良好的时频聚焦性，特别适用于非线性、非平稳信号的分析处理。目前，基于EMD的机械故障诊断方法取得了很大的进展[2]。然而，随着对EMD方法研究的深入，其也存在着一些问题，主要有：1)EMD是一种经验性的方法，它缺乏完备的理论基础，比如，到目前为止，EMD方法分解所得到IMF分量正交性尚未得到理论上的证明。2)EMD由于收敛条件不合理、过包络和欠包络等问题易出现模态混叠。3)EMD要分解出一个IMF分量，需要进行多次的迭代，因此，要得到一个实际信号的所有IMF分量需要较长的计算时间。针对EMD存在的不足，Gillesll[3]结合EMD的白适应性和小波分析的理论框架，提出了一种新的自适应信号处理方法即经验小波变换(empiricalwavelettransform．EWT1。该方法的核心思想是通过对信号的频谱进行自适应划分，构造合适的正交小波滤波器组以提取具有紧支撑傅里叶频谱的AM—FM成分，然后，对提取出的AM．FM模态进行Hilbert变换，得到有意义的瞬时频率和瞬时幅值，进而可以得Hilbert谱。由于EWT是在小波框架下建立的方法，故其理论充分且其计算量远小于传统的EMD方法。本文将EWT引入到机械故障诊断中，提出了一种基于EWT的机械故障诊断方法，并与EMD方法进行了对比分析，仿真结果验证了EWT方法的有效性。最后，将该方法成功地应用到转子不同碰磨严重程度的故障数据分析中，实验结果进一步验证了该方法的有效性。2经验小波变换经验模态分解方法的实质是把一个信号)(tf分解成N个本征模态函数(IMF)kc和残差nr之和：nkNkrctf1)(（1）而文献[13]与原始文献[1]稍微有所不同，经验小波变换的目的是把信号)(tf分解成N+1个本征模态函数)(tfk之和：）tftfkNk()(0(2)这里，一个固有模态函数定义为调幅．调频信号(AM—FM)：(3)在EMD方法中最感兴趣的是如何自适应地从原始信号中筛选出调幅-调频成分）tfk(。传统的EMD方法由于终止条件不合理、欠包络和过包络等问题而造成模态混叠，并且EMD的筛分过程也缺乏可靠的理论依据。为克服此问题，GillesE在EMD的基础上，结合小波分析，提出了经验小波变换方法，该方法是建立在小波分析的理论框架上，根据信号的Fourier谱特性自适应地选择一组小波滤波器组来提取信号的不同AM-FM成分。为了选择合适的小波滤波器组，需要对Fourier谱进行自适应地分割。假设将Fourier支撑[0，л]分割成N个连续的部分，用表示各片段之间的边界，n选择为信号Fourier谱两个相邻极大值点之间的中点(0=0，n=兀)，如图1所示，则每段可以表示为：以每个n为中心，宽度为nnT2定义了一个过渡段，如图1所示的阴影区。确定分割区间n后，经验小波定义为每个n上的带通滤波器，Gilles根据Meyer小波[4]的构造方法构造经验小波，经验小波函数)(n和经验尺度函数)(n定义如下：在EWT方法中，如何分割傅里叶频谱是至关重要的，因为它直接关系到自适应分解的结果。不同的频谱部分对应于以不同特定紧支撑的频率为中心的模态。除去0和¶还需要找到N-1个边界。要找到这样的边界，就要确定N，文献[12]给出了一种确定N的阈值法。为频率域范围内检测到的极大值点的幅值，然后按照递减的规律排列并且归一化到[0，1]，取M+(M1一Mm)为阈值，其中为相对振幅比，∈(0，1)，对于确定的可以令大于阈值的极大值点个数为Ⅳ，并取前Ⅳ个最大的极大值点求边界。为此，可以采用类似于传统的小波变换方法来定义经验小波变换。假设傅里叶变换和逆变换分别记为F[x]和F[x+1]，则细节系数由经验小波函数与信号内积产生：近似系数通过尺度函数与信号内积产生：式中：)(n和)(n分别为经验小波函数和尺度函数。原信号重建如下：式中：符号表示卷积，),0(n和),(nn分别表示),0(tWn和),(tnWn的Fourier变换。经验模态）tfk(定义如下：通过经验小波变换，得到一个信号的经验模态后，就可以对每个经验模态函数进行Hilbert变换，从而得到有意义的瞬时频率和瞬时幅值，进而得到Hilbert谱。3仿真研究考查如下两个仿真信号X1(f)和2(f)，仿真信号X1(t)由调幅调频信号x11(f)和调幅信号12(f)组成，其时域波形如图2(a)所示。这里，采样频率Fs=l024Hz。采样点数是512。仿真信号x2(t)由单调趋势信号X21(f)，调频信号X2(f)和两个不同时间段频率单一的和信号X3(f)信号组成，其时域波形如图2(b)所示。经验小波分解的输出是由一个尺度函数和Ⅳ个小波函数分别滤波的结果。对于信号X1(f)和X2(f)，分别取N=2和N=3。图3显示出了两个仿真信号的频谱和每个滤波器所支撑的已检测到的边界。图4为经验小波变换结果，为了对比分析，在此也给出了EMD的分解结果，如图5所示。由图4(a)和5(a)可以观察到仿真信号Xl(f)的2个模态被有效地分解出来，但是相对于经验小波变换，EMD分解出了过多的模态，且分离出一些本来是相同成份的部分信息，除了高频部分，很难解释这些原来信号中不存在的虚假模态。由图4(b)可知，经验小波变换将最初属于同一个部分的两个模态分解出来，因为这两个模态明显有各自的能量，并且可以被视为两个独立的模态，而由图5(b)可知，EMD则分解不出来，且分解出了过多的虚假模态。因此，可以得出，经验小波变换能够有效地检测出频谱中存在的模态，并能够将模态分解出来。理论上，经验小波变换有坚实的理论基础，可以通过小波变换的方法去理解该方法。然而，EMD缺乏完备的理论基础，例如，到目前为止，EMD方法分解所得到IMF分量正交性尚未得到理论上的证明，存在模态混叠问题，而模态的混叠，影响了瞬时频率的表达。EMD的模态混叠产生的主要原因有以下3点：①从模态分解的终止条件来看，如果所选用的终止条件不合理，那么模态分解就会不充分。②筛选过程中采用了三次样条插值的方法，由于三次样条的欠冲和过冲的缺陷，同样也会造成模态之间的混叠。③从信号特性的角度上看，当信号的时间尺度存在跳跃性变化时，会产生模态的混叠现象。另外，传统的EMD方法要分解出一个IMF分量，需要多次的迭代，计算量大，而EwT是在小波框架下建立的方法，故其计算量远小于传统的EMD方法。由此可知，在处理非线性、非平稳信号方面，EWT方法明显优于EMD方法。4实验研究为进一步验证经验小波变换的有效性，将经验小波变换应用到双盘转子的碰磨数据的分析中，碰磨转子实验台如图6所示[5]，转子由电动机驱动，轴承为滑动轴承，用非接触式电涡流传感器测量垂直与水平方向的振动。转子径向碰摩故障通过以下装置来模拟，通过置换不同内径的定子，可以模拟不同程度的碰摩故障，如图7所示[5]。在转速为3000r／min工况下，采样频率为1．6kHz，采样点数为1024。利用传感器获得的两种不同碰磨程度下的振动信号及其频谱图.从图9中可以看出，时频图得到了很好的分割，主要的倍频都处于支撑边界的中间，且严重碰磨的分割段数Ⅳ比轻微碰磨大。从图10(a)可以看出，轻微碰磨时，高阶频率分量F3表现为周期性的冲击信号，但相对于低阶频率分量明显非常微弱，从图11(a)可知，轻微碰磨的Hilbert谱在1X和2X分量反映明显。幅值基本上比较稳定，且持续存在，而在高阶部分主要集中在6X和8X之间，十分微弱，都是周期性地被激发。通过以上分析可知，EwT变换能够按照频率特征有效的从低频到高频自适应地分解碰磨故障信号。Hilbert[6]谱能够很好的反映出碰