信号分析与检测技术实验报告

1、信号分析与检测技术实验课实验报告专业班级：姓名：学号： 可靠性与系统工程学院2014年6月实验一滚动轴承故障检测与信号分析实验报告一、实验目的与要求1.1 实验目的:1. 了解振动信号采集、分析与处理的整个过程及注意事项；2. 了解并掌握测试仪器的连接、信号的敏感参数选取、测点布置及各注意事项；3. 掌握信号的时域分析、频域分析理论与特点。4. 了解不同形式的故障轴承5. 了解利用振动进行轴承故障诊断原理及方法6. 了解轴承故障的危害性及其表现形式。7. 掌握信号的调制与解调原理与方法1.2 实验要求：1、实验前正确校准系统2、正确布置测点位置3、选取合适的采样参数4、实验室空间比较拥挤，请大

2、家有秩序地进行试验。二、实验原理及结果分析1、试说明什么是采样频率fs，Shannon采样定理是什么？另外，信号采样时，采样频率fs、频率分辨率f、采样数据长度N和采样时间T之间的关系。采样频率fs，也称为采样速度或者采样率，是每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。为保证采样后信号能真实地保留原始模拟信号信息，信号采样频率必须至少为原信号中最高频率成分的2倍，亦称仙农（香农shannon）定理。Xs(Fi)，Fi = i *Fs/N , i = 0，1，2，.，N/2f=Fs/N=1/(t*N)f = 1 /T其中f为频谱分辨率，T为总采样时间，N为采样数据长度

3、。2.试分析一下正常情况和滚动轴承故障时的时域特征的区别？本次实验对滚动轴承正常工作状态下的工况数据进行了采集及分析。其具体时域分析特征指标值如下表所示：序号指标量纲X方向Y方向Z方向EU工程单位EUmVmVmV1最大值(EU)1.087921.592640.555332最小值(EU)-1.05419-1.6246-0.613773平均值(EU)-0.01562-0.00306-0.00314平均幅值(EU)0.16240.228480.09775方根幅值(EU)0.136470.19090.082416有效值(均方根)(EU)0.206840.294240.123497均方值(EU)20.0

4、42780.086580.015258标准差(均方差)(EU)0.206250.294230.123459最大值-平均值(EU)1.103551.595710.5584410最小值-平均值(EU)-1.03857-1.62153-0.6106711方差(EU)20.042540.086570.0152412偏度指标(EU)32.853E-41.536E-4-1.53E-613峭度指标(EU)40.006180.027657.374E-414偏态因数(无)0.032510.00603-8.14E-415峰态因数(无)0.41490.690290.1748716波形因数(无)1.273621.28

5、7831.2638917脉冲因数(无)6.698737.110386.2816718峰值因数(无)5.259585.521214.9700919裕度因数(无)7.971878.510054.97009本次实验对滚动轴承内环故障工作状态下的工况数据进行了采集及分析。其具体时域分析特征指标值如下表所示：序号指标量纲X方向Y方向Z方向EU工程单位EUmVmVmV1最大值(EU)3.76946.381564.236732最小值(EU)-3.44688-5.31059-3.684593平均值(EU)-0.01049-0.02181-0.022924平均幅值(EU)0.348250.476390.2693

6、5方根幅值(EU)0.268870.368790.192956有效值(均方根)(EU)0.520770.72110.455177均方值(EU)20.27120.519990.207188标准差(均方差)(EU)0.520670.720770.454599最大值-平均值(EU)3.779896.403374.2596610最小值-平均值(EU)-3.43639-5.28878-3.6616711方差(EU)20.271090.519510.2066512偏度指标(EU)30.00814-0.007750.0175113峭度指标(EU)40.561942.558320.5160114偏态因数(无)

7、0.05772-0.02070.1864415峰态因数(无)4.646166.478869.0826616波形因数(无)1.495381.513681.6901617脉冲因数(无)10.823613.395615.73218峰值因数(无)7.238048.84979.3079719裕度因数(无)14.01917.303821.9574传感器Z轴方向为轴向，需要对X或Y方向进行分析。在上表中，相关指标分析如下：1）X轴上平均值正常时为-0.01562，内环故障时为-0.01049；平均值反映了信号变化的中心趋势，对比说明在滚动轴承发生内环故障时，振动变化向X轴正方向发生了比较明显的偏移。2）有效

8、值X 轴上正常时为0.20684，内环故障时为0.52077；有效值用以描述振动烈度，对比说明滚动轴承发生内环故障时振动强度明显增大。3）方差X 轴上正常时为0.04254，内环故障时为0.27109；方差反映了信号绕均值的波动程度，对比说明滚动轴承发生内环故障时信号绕均值波动程度明显增大。4）峭度指标X 轴上正常时为0.00618，内环故障时为0.56194；峭度指标可以敏感捕捉信号中的冲击成分，对比说明滚动轴承发生内环故障时收到冲击明显。5）偏态因数X 轴上正常时为0.03251，内环故障时为0.05772；偏态因数用来检验信号偏离正态分布的程度，对比说明滚动轴承发生内环故障时信号较为明显

9、的偏离正太分布。6）峰值因数X 轴上正常时为5.25958，内环故障时为7.23804；峰值因数表示波形是否有冲击的指标，对比说明滚动轴承发生内环故障时发生冲击。时域指标诊断不易直接看出所包含信息与故障的联系，经常用于做定性分析。信号的时域统计分析（包括均值、均方根值、方差、概率密度函数等），在简单情况下可识别系统的状态。根据以往经验以及对滚动轴承正常时和内环故障时获得的数据对比，当滚动轴承发生故障时，时域指标中，峰值和峰值因数会变大，峭度因子也会相应变化。3、振动信号的频谱分析主有哪些？各有何特点？试画出实验时正常和滚动轴承故障时的振动时域波形及其幅值谱图，并给出matlab程序。 时域分析

10、法：1）简单、方便；2）当振源较为复杂时，只能定性诊断，即指出设备有无异常；而无法指出具体的故障部位 频域分析法：各个部件在运转时，都以一定的特征频率进行振动，特征频率所对应的幅值即代表了该部件的振动大小。通过特征频率，可将各个部件的振动分开，实现定位诊断。 倒频域分析法：倒频谱运算，即将复杂的功率谱为一系列卷积或乘积的形式，通过取对数转化为简单的和式形式，便于识别信号的主要分量。优点为：1）将检测信号的功率谱与系统传递函数分离开来；2）便于将主要的信息功率谱从复杂的边频中识别出来。 幅值域分析法：根据概率密度函数曲线的形状进行诊断。当设备发生故障时，振动加剧，大幅值的数目增多，其比例值增大，

11、曲线形状由“尖、瘦”变为“矮、胖”。 相关域分析法：，主要是对信号的自相关函数R（）进行求解，并通过R（）对信号中的周期成分进行分析，从而找出故障振源。主要对周期性的故障信号进行检测。利用滚动轴承正常工况下的实验数据，得到振动时X、Y、Z三通道时域波形图如下： 图1 X通道时域波形图 图2 Y通道时域波形图 图3 Z通道时域波形图利用滚动轴承正常实验中获取的频率峰值数据，得到X、Y、Z幅值谱图如下：图4 X通道幅值谱图 图5 Y通道幅值谱图图6 Z通道幅值谱图利用滚动轴承内环故障工况下的实验数据，得到振动时X、Y、Z三通道域波形图如下： 图7 X通道时域波形图 图8 Y通道时域波形图图9 Z通

12、道时域波形图利用滚动轴承内环故障实验中获取的频率峰值数据，得到X、Y、Z幅值谱图如下： 图10 X通道幅值谱图 图11 Y通道幅值谱图图12 Z通道幅值谱图相关matlab程序如下：plot(a0,a1);xlabel(t/s);ylabel(幅值);title(正常滚动轴承振动时域波形);grid;plot(a0,a2);xlabel(t/s);ylabel(幅值);title(正常滚动轴承振动时域波形);grid;plot(a0,a3)xlabel(t/s);ylabel(幅值);title(正常滚动轴承振动时域波形);grid;fs=10240;%采样频率ts=1/fs;%采样时间间隔N

13、=1024;%采样点数n=0:N-1;y1=fft(a1,1024);%对X 轴信号进行fft 变换mag1=abs(y1)/512;%求幅值乘上后面的2/N 得到幅值f=(0:length(y1)-1)*fs/length(y1);%进行对应的频率转换plot(f(1:N/2),mag1(1:N/2);%做频谱图xlabel(频率/Hz);ylabel(幅值);title(正常滚动轴承振动幅频谱图);grid;y2=fft(a2,1024);%对X 轴信号进行fft 变换mag2=abs(y2)/512;%求幅值乘上后面的2/N 得到幅值f=(0:length(y2)-1)*fs/lengt

14、h(y2);%进行对应的频率转换plot(f(1:N/2),mag2(1:N/2);%做频谱图xlabel(频率/Hz);ylabel(幅值);title(正常滚动轴承振动幅频谱图);grid;y3=fft(a3,1024);%对X 轴信号进行fft 变换mag3=abs(y3)/512;%求幅值乘上后面的2/N 得到幅值f=(0:length(y3)-1)*fs/length(y3);%进行对应的频率转换plot(f(1:N/2),mag3(1:N/2);%做频谱图xlabel(频率/Hz);ylabel(幅值);title(正常滚动轴承振动幅频谱图);grid;4、计算轴承的故障特征频率，

15、并填入下表。轴承外圈固定，内圈（即轴）的旋转频率为fs，n为轴的转速，轴承节径为D，滚动体直径为d，接触角为，滚动体个数为z。假定滚动体与内外圈之间纯滚动接触，则滚动轴承的特征频率计算公式如下：内圈（即轴）的旋转频率： 外环故障频率： 内环故障频率： 滚珠故障频率： 保持架碰外环： 保持架碰内环： 内径（mm）外径(mm)节圆直径(mm)滚动体直径(mm)接触角(deg)滚动体数目(个)255238.56.753012根据以上数据及相应公式，可得到下表1中数据：表1转速n(r/min)外环故障频率内环故障频率滚珠故障频率保持架碰外环故障频率保持架碰内环故障频率2782.8255.7347.24258.4419.6726.715、根据实验数据，将分析结果填入表2-3，并仔细观察振动信号的时域波形和频谱特征