基于小波分析的旋转机械振动信号定量特征研究

基于小波分析的旋转机械振动信号定量特征研究1侯敬宏黄树红申弢张燕平（华中科技大学能源与动力工程学院武汉430074）摘要：通过对机械振动信号的连续小波变换，利用小波滤波器良好的时频特性，研究了振动信号经过连续小波变换后的统计特征。在信号的特征提取中，引入“灰度矩”并把一阶矩作为定量指标。对8种典型故障信号的研究表明，这种方法能够简单有效地提取信号的特征，区分振动故障。关键词:小波变换；故障诊断；信号处理；特征提取中图分类号TH165.3文献标识码A现有的旋转机械振动信号的分析方法多只是定性指出机组可能存在的故障，而现实中并不是设备存在隐患就立即停机处理，而是希望根据定量评价故障的危害程度以及发展趋势而作决定，从而有效避免“维修不足”和“过剩维修”所带来的经济损失。小波滤波器是一个具有恒Ｑ(品质因子，定义为滤波器的中心频率与带宽之比)特性的滤波器，因此，其可以随信号频率增高而减小时窗宽度并相应增加频窗宽度，即具有”变焦”分析的特性。故而小波分析在信号的时频分析中具有特殊的地位和性能。从数学上看，对机械振动信号进行小波分析，是将一维信号在二维空间进行分解，使得在一维空间中无法提取的特征在二维空间中得到很好的反映。从物理上看，对振动信号进行连续小波变换和分析，实际上就是将无法从一维时域或者频域中体现的信号特征，在具有良好滤波器特性的小波时频窗中得以体现，以获得信号所包含的特征信息[3]。因此，分析机械振动信号的连续小波系数，从中提取特征，可以更好地反映信号的本质。利用连续小波分析信号时，一般用每个尺度下小波能谱图提取信号的定性特征。这在很大程度上依赖于人的主观判断。如果要实现自主式的状态诊断，则必须要提取能反映信号特征的定量指标。本文引入一种“灰度矩”的统计量，对8种故障信号波形经连续小波变换后的系数矩阵进行处理，提出区分这些故障的定量指标。所选取的8种故障信号为[2]：不平衡，不对中，油膜振荡，轴裂纹，支座松动的轴振动信号，支座松动的支座振动信号，碰摩故障的轴振动信号，碰摩故障的支座振动信号。本文所采用的所有故障信号均为本单位模拟试验台的实验数据。一．连续小波变换一般所讨论的小波，是指一个被称之为母小波或基本小波的函数，经伸缩和平移所产生的函数簇Ψa，b（t），式中a0是尺度因子，b反映位移，其值可正可负。定义满足“容许条件”[1]：d02|)(|……①的函数Ψ（t）为基本小波函数，简称小波函数。若函数f(t)是平方可积函数，即dttf2|)(|，定义：WTx(a，b)=dtabttfa)()(1=f(t)，Ψa，b(t)……②为f(t)的小波变换。上标*代表共轭，x，y代表内积。要实现WTx(a，b)所定义的内积，除少数情况可做解析计算外，大多数情况只能通过计算机做近似数值计算。按数字信号处理的习惯可将式②变成：WTx(a，k)=])([)(aknnfaTn……③1国家自然科学基金项目（50105004）2对每一个固定a值，依次求不同k值下的乘积和，便得到该a值下一组WT系数。f(t)如果是由512个离散点组成的时间序列，每一个给定的a值（尺度），有512个系数值，N个尺度下经小波变换后便组成一个N×512的系数矩阵。二．小波系数矩阵的“灰度矩”信号经连续小波变换后的小波系数往往以灰度图的形式呈现出来。对于旋转机械的振动信号来讲，其故障特征周期出现，可用灰度图表示，如图1所示。显然，从灰度图上可直观分辨不同的故障，但实际故障诊断时，人为判断会存在模糊性和个体差异。那么，能否用一种定量指标来区分灰度图之间的总体差异呢？2004006008001000030609012020040060080010000306090120图1用morl小波分析的不对中与油膜振荡故障的128尺度下的灰度图本文定义一个m×n矩阵的k阶“灰度矩”来定量描述小波灰度图的差异：Mk=minjkijjianm1122)1()1(1……④这里的权值22)1()1(ji表示元素aij与a11之间的“距离”，相当于灰度图中某一个像素点与参考点之间的几何长度。④式不直接采用灰度图中的几何距离，是为了不受灰度图几何尺寸的影响。图2是对应8种振动故障信号morl小波系数矩阵的灰度矩，其中灰度矩的阶数为k=1~50。可以看出，高阶矩与低阶矩对故障的区分度相差不大，因此可以采用一阶灰度矩M1来区分故障的灰度图，以减少计算量，有利于实时诊断。三．故障信号灰度矩M1的处理方法图2同时给出了另外一个结果，即morl小波的灰度矩可以很好地区分8种振动故障。显然，这对于定量故障诊断是极有吸引力的。如果能找到几种故障区分性很好的小波，就有可能从多个角度进行定量故障诊断，减少由于信号畸变而导致某一个小波灰度矩不稳定造成的误诊断。为此，本文进行了以下工作：对采集到的故障信号(假设由512个离散点组成)进行归一化处理，若处理后的序列为{x1，x2，……，05101520253035404550020406080100120M1k图2用morl小波分析8种故障的Mk(k=1~50)分布油膜振荡碰摩-支座松动-支座轴裂纹不对中碰摩-轴不平衡松动-轴3x512}，有151212iix成立；1.选择8种常用的连续小波db4、sym3、coif4、bior6.8、rbio6.8、dmey、meyr、morl分别对每一种故障进行N尺度连续小波变换，得到小波系数矩阵coefs(N×512矩阵)，本文采用MATLAB进行计算；2.对小波系数求绝对值；3.依④式求得coefs矩阵的1阶“灰度矩”M1。之所以对故障信号进行归一化处理，是因为故障的原始波形振幅差异较大，经归一化处理以后，相当于在保证能量守恒的前提下，把不同波形限制在一定的范围内，增强不同波形之间的可比性，使研究的结果具有普遍意义。如果以8种小波为横轴，M1的大小为纵轴，把每一种故障的8M1个用线连起来，便得到对应故障的灰度矩分布线。图3是64尺度的灰度矩分布线。由图3可以看出，8种故障被灰度矩分布线很好的区分开了，只有轴裂纹和碰摩-轴2种故障基本重合在一起。而且db4，sym3，dmey，meyr，morl小波区分度较好。在以下的分析中，将采用dmey、meyr、morl三种小波。如将尺度再增加1倍，改为128尺度，只用dmey，meyr，morl小波进行分析，M1的分布线如图4所示。从图4中可以看出，从上到下8种故障线的相对位置没有发生改变，说明尺度的变化对灰度矩影响不大，这个结果非常有利于故障识别。四.典型故障小波灰度矩M1的分布在实际情况下，同一种典型故障会有不同的时域波形，比如，同一台设备在不同转速下的波形是不同的，二台相同故障的设备的时域波形也可能是不同的。这是因为影响振动的还有其他因素。但是对应同一种故障，不同时刻或位置采集的振动波形应有反映故障的共性特征。那么，在固定尺度的情况下，一定故障是否有一定的灰度矩数值区间呢？如果不同故障的灰度矩区间互不重叠，则可以用信号的小波灰度矩准确识别故障。dmeymeyrmorl05101520255M1图48种故障经3种小波分析后的128尺度M1分布线油膜振荡碰摩-支座松动-支座轴裂纹不对中碰摩-轴不平衡松动-轴Db4Sym3Coif4Bior6.8Rbio6.8dmeymeyrmorl05101520253035M1图38种故障经3种小波分析后的64尺度M1分布线油膜振荡碰摩-支座松动-支座轴裂纹不对中碰摩-轴不平衡松动-轴41.不平衡故障分别提取3种不同转速下的不平衡故障曲线如图5~图7。经dmey、meyr、morl这3种小波分析后，对应的灰度矩分布如图8所示。M1的大小大约在4~6的范围内。2.油膜振荡提取转速分别为9730rpm、9751rpm、9772rpm下的油膜振荡的波形，其中9730rpm的波形及频谱如图9所示。它们经上述3种小波分析后，灰度矩M1的分布如图10。由图10可以看出油膜振荡故障在三种转速下M1的大小集中在18~32之间。3.8种故障信号的灰度矩分布经过相同的计算，8种故障利用上述小变换及频谱后M1的变化区间如表-1所示。将表-1数据绘成灰度矩分布图，参见图11、图12。图11是几种故障的轴振动信号的M1分布区间图。可以看出，几种故障一阶灰度矩的区间并不重合。如果振动信号的灰度矩落入某故障区间，即可初步判断其故障类型。图12上方两条灰度矩分布线分别是碰摩和松动故障的轴承座振动信号灰度矩，它们也有很好的区50100150200250300350400450500-0.02-0.0100.010.0201x00.010.020.030.040.05图-7转速4569下的不平衡故障波形及频谱50100150200250300350400450500-0.0500.050x2100.050.10.150.2图-9转速9730rpm下的油膜振荡波形及频谱50100150200250300350400450500-0.0200.0201x00.020.040.060.08图-5转速4439下的不平衡故障波形及频谱50100150200250300350400450500-0.02-0.0100.010.0201x00.020.040.06图-6转速4496下的不平衡故障波形及频谱M1456944964439图83组不同转速不平衡故障的M1分布线dmeymeyrmorl4546M11020430977297519730图103组不同转速油膜振荡信号的M1分布线dmeymeyrmorl图5转速4439rpm下的不平衡故障波形及频谱图6转速4496rpm下的不平衡故障波形及频谱1X1X1X1/2X图7转速4569rpm下的不平衡故障波形及频谱图9转速9730rpm下的油膜振荡波形及频谱图83组不同转速下不平衡故障的M1分布线图103组不同转速下油膜振荡故障的M1分布线5分度。图12下方的阴影区是不对中故障的轴振动信号灰度矩。可以看出不对中故障的M1区域覆盖了图11中碰摩-轴、不平衡、松动-轴以及轴裂纹的大部分区域，说明这几种故障具有与不对中故障相似的信号特征。可见，还需进一步提取有关征兆和改进分析方法，将这几种故障与不对中故障再区分。这正是作者下一步拟开展的工作。表-18种典型故障在不同情况下的M1值五.结论为了对机械振动信号连续小波系数进行总体定量刻画，并从中提取信号特征，引入了小波系数矩阵的k阶“灰度矩”的概念。研究表明，小波系数矩阵的一阶“灰度矩”能够很好地表征机械振动信号的特征，并且能定量地描述旋转机械的故障，有望成为新的、有效的故障诊断工具，其意义是重大的，在机械故障诊断中有良好的应用前景。参考文献[1]杨福生.小波变换的工程分析与应用.科学出版社，2000[2]申弢.大型旋转机械智能监测诊断中信息融合理论与技术的研究.华中理工大学博士学位论文，1999[3]刘刚，屈梁生.机械信号连续小波系数的统计特性研究.西安交通大学学报.第36卷第3期，2002[4]张志涌.精通MATLAB5.3版.北京航空航天大学出版社，2000典型故障转速统计指标M1值dmeymeyrmorl油膜振荡973032312897512322209772202018碰摩支座346225242235182524223618252422松动支座459615151346521515134683151513轴裂纹3041141413317010.810.21032698.68.27.8不对中2419121211.532744.243.741995.35.24.8碰摩-轴34628.387.335188.387.336188.387.3不平衡443965.95.244965.35.24.84569554.4松动轴45964.34.23.846524.34.23.846834.34.23.8不对中松动-支座摩碰-支座M1dmeymeyrmorlM1322824201612840M1322824201612840图12支座振动