基于振动信号的齿轮故障诊断方法研究硕士毕业论文..doc

基于EMD的齿轮故障诊断 一 基于EMD（经验模式分解）的振动信号特征提取 1.1 EMD的研究背景 齿轮摩擦是齿轮箱系统中常见的非线性振动故障。当齿轮存在局部摩擦时, 摩擦信号是在动静摩擦的短时间内产生的，具有频带宽的特点。在采集振动信号时混有宽频带随机噪声的噪声信号, 其频带相互叠加。在振动信号中还含有能量较大的、与转速有关的背景信号, 所以利用传统的滤波方法很难仅保留含有故障信息的摩擦信号而将背景信号和噪声信号滤掉。由于用小波分析方法进行故障信号特征提取存在缺点，美籍华人NordenE.Huang等人针对瞬时频率的概念进行深入研究后,创造性地提出了本征模式函数(Intrinsic Mode Function,IMF)的概念以及将任意信号分解为本征模式函数组成的新方法——基于经验模式分解(Empirical Mode Decomposition,EMD)方法，从而赋予了瞬时频率合理的定义、物理意义和求法，初步建立了以瞬时频率为表征信号交变的基本量，以本征模式分量为时域基本信号的新的时频分析方法体系。EMD分解算法的基本思想是：对一给定信号，先获得信号极值点，通过差值获得信号包络，获得均值，与均值的差得到分解的一层信号；如此重复，获得分解结果：，即个IMF和一个残值。由于EMD方法具有自适应特性，适宜于非线性、非平稳信号的分解，该方法应用于齿轮的故障诊断分析中。结果表明，该方法能够突出齿轮振动信号的故障特征，从而提高齿轮故障诊断的准确性[19]。为此，本文将EMD方法应用到齿轮的故障诊断中，利用EMD方法对齿轮信号进行分解，得到若干个IMF，然后对包含故障信息的IMF作功率谱分析，有效提取故障信号特征信息。 1.2 经验模式分解(EMD)方法原理 经验模式分解方法的主要目的是通过对非线性、非平稳信号的分解，获得一系列表征信号特征时间尺度的本征模式函数和一个残余分量，各个IMF是单分量的幅值或者频率调制信号。每个IMF需要满足以下几个条件信号：(1)整个信号中零点数与极点数相等或者相差1；(2)信号上任意一点，由局部极大值点确定的包络线和由局部极小值点确定的包络线的均值均为零，即信号关于时间轴局部对称[20]。 分解过程如下： 设信号为，取其上下包络局部均值组成的系列为，令它们的差值为 (4·1) 对于非线性、非平稳数据来说，一般一次处理不会形成基本模式分量，一些非对称波依然存在。因此要把看作待处理的数据继续重复上述操作。取的均值系列为，得到它们之前的差值为 (4·2) 重复k次操作，就可以得到： (4·3) 当满足基本模式分量的条件后，就会获得了第一个基本模式分量 (4·4) 将基本模式分量从信号中分离出来，可以得到剩余系列 (4·5) 把作为新的数据按上述操作进行处理，依次类推，可得到： …… (4·6) 这个处理过程在满足预先设定的停止标准后就可停止。 将式(4·5)和式(4·6)相加后得到 (4·7) 这样就把一个信号分解为个基本模式函数和一个残余分量[19,20]，其中分解出的个分量分别包含信号从高频到低频的不同频率段成分，而剩余分量是原始信号的中心趋势。 (4·6)由于每个固有模态函数为平稳信号，由式(4·7)可以定义每个固有模态分量的自相关函数为 (4·8) 1.3 基于EMD的振动信号特征提取分析 当齿轮发生局部摩擦故障时，设备振动信号中将出现周期性的冲击，这些周期性的冲击信号由于噪声等干扰，所以很难通过频谱分析等传统方法加以识别，因此本文提出了一种基于EMD的冲击信号提取方法。依照IMF的定义，周期性的冲击信号无法以独立的IMF形式被分解出来，但它会以类似冲击响应信号的形式出现在高频段的IMF中。 假设如下仿真信号： (4·9) 该仿真信号是由一个50Hz的正弦信号和与之同频率的周期性冲击信号构成的，并加入了一定量的均匀白噪声。采样频率为2000Hz单位为微米（）。图4·1表示了该仿真信号的波形及频谱。 图4·1该仿真信号的波形及频谱图 图4·2 表示载入的同频率的冲击信号。 图4·2 50Hz的冲击信号 图4·3表示一定量的白噪声。 图4