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

基于小波分析的旋转机械振动信号定量特征研究 侯敬宏 黄树红 申 弢 张燕平 （华中科技大学能源与动力工程学院 武汉 430074） 摘要： 通过对机械振动信号的连续小波变换，利用小波滤波器良好的时频特性，研究了振动信号经过连续小波变换后的统计特征。在信号的特征提取中， 引入“灰度矩”并把一阶矩作为定量指标。对8种典型故障信号的研究表明，这种方法能够简单有效地提取信号的特征，区分振动故障。 关键词: 小波变换；故障诊断；信号处理；特征提取 中图分类号 TH165.3 文献标识码 A 现有的旋转机械振动信号的分析方法多只是定性指出机组可能存在的故障，而现实中并不是设备存在隐患就立即停机处理，而是希望根据定量评价故障的危害程度以及发展趋势而作决定，从而有效避免“维修不足”和“过剩维修”所带来的经济损失。 小波滤波器是一个具有恒Ｑ(品质因子，定义为滤波器的中心频率与带宽之比)特性的滤波器，因此，其可以随信号频率增高而减小时窗宽度并相应增加频窗宽度，即具有”变焦”分析的特性。故而小波分析在信号的时频分析中具有特殊的地位和性能。从数学上看，对机械振动信号进行小波分析，是将一维信号在二维空间进行分解，使得在一维空间中无法提取的特征在二维空间中得到很好的反映。从物理上看，对振动信号进行连续小波变换和分析，实际上就是将无法从一维时域或者频域中体现的信号特征，在具有良好滤波器特性的小波时频窗中得以体现，以获得信号所包含的特征信息[3]。因此，分析机械振动信号的连续小波系数，从中提取特征，可以更好地反映信号的本质。利用连续小波分析信号时，一般用每个尺度下小波能谱图提取信号的定性特征。这在很大程度上依赖于人的主观判断。如果要实现自主式的状态诊断，则必须要提取能反映信号特征的定量指标。 本文引入一种“灰度矩”的统计量，对8种故障信号波形经连续小波变换后的系数矩阵进行处理，提出区分这些故障的定量指标。所选取的8种故障信号为[2]：不平衡，不对中，油膜振荡，轴裂纹，支座松动的轴振动信号，支座松动的支座振动信号，碰摩故障的轴振动信号，碰摩故障的支座振动信号。本文所采用的所有故障信号均为本单位模拟试验台的实验数据。 连续小波变换 一般所讨论的小波，是指一个被称之为母小波或基本小波的函数，经伸缩和平移所产生的函数簇Ψa，b（t），式中a0是尺度因子，b反映位移，其值可正可负。定义满足“容许条件”[1]： ……① 的函数Ψ（t）为基本小波函数，简称小波函数。 若函数f(t)是平方可积函数，即，定义： WTx(a，b)==f(t)，Ψa，b(t) …… ② 为f(t)的小波变换。上标*代表共轭，x，y代表内积。 要实现WTx(a，b)所定义的内积，除少数情况可做解析计算外，大多数情况只能通过计算机做近似数值计算。按数字信号处理的习惯可将式②变成： WTx(a，k)= …… ③ 对每一个固定a值，依次求不同k值下的乘积和，便得到该a值下一组WT系数。f(t)如果是由512个离散点组成的时间序列，每一个给定的a值（尺度），有512个系数值，N个尺度下经小波变换后便组成一个N×512的系数矩阵。 小波系数矩阵的“灰度矩” 信号经连续小波变换后的小波系数往往以灰度图的形式呈现出来。对于旋转机械的振动信号来讲，其故障特征周期出现，可用灰度图表示，如图1所示。显然，从灰度图上可直观分辨不同的故障，但实际故障诊断时，人为判断会存在模糊性和个体差异。那么，能否用一种定量指标来区分灰度图之间的总体差异呢？ 本文定 义 一 个 m×n 矩阵的k阶“灰度矩”来定量描述小波灰度图的差异： Mk= …… ④ 这里的权值表示元素aij与a11之间的“距离”，相当于灰度图中某一个像素点与参考点之间的几何长度。④式不直接采用灰度图中的几何距离，是为了不受灰度图几何尺寸的影响。 图2是对应8种振动故障信号morl小波系数矩阵的灰度矩，其中灰度矩的阶数为k=1~50。可以看出，高阶矩与低阶矩对故障的区分度相差不大，因此可以采用一阶灰度矩M1来区分故障的灰度图，以减少计算量，有利于实时诊断。 故障信号灰度矩M1的处理方法 图2同时给出了另外一个结果，即morl 小波的灰度矩可以很好地区分8种振动故障。显然，这对于定量故障诊断是极有吸引力的。如果能找到几种故障区分性很好的小波，就有可能从多个角度进行定量故障诊断，减少由于信号畸变而导致某一个小波灰度矩不稳定造成的误诊断。为此，本文进行了以下工作： 对采集到的故障信号(假设由512个离散点组成)进行归