《故障诊断方法与应用综述.docx》由会员分享，可在线阅读，更多相关《故障诊断方法与应用综述.docx（16页珍藏版）》请在三一文库上搜索。
1、课程名称：故障诊断方法与应用报告题目： 内圈故障诊断实验报告学生班级；研152学生姓名：任课教师：学位类别：评分标准 及分值选题与参阅资料（分值10）报告内容（分值 60 ）报告表述（分值20 ）创新性（分值10 ）评分报告评语：总评分评阅教师：评阅时间年 月日故障诊断实用技术摘要设备故障诊断技术是一种了解和掌握设备在使用过程中的状态，确定其整体或局部是正常或异常，早期发现故障及其原因，并能预报故障发展趋势的技术。安装合适的传感器可以获得故障的特征信号， 通过信号反映故障产生原因。滚动 轴承是机械中的易损元件，据统计旋转机械的故障有 30溢由轴承引起的，它的 好坏对机器的工作状态影响极大。轴承
2、的缺陷会导致机器剧烈振动和产生噪声， 甚至会引起设备的损坏。滚动轴承的振动可由于外部的振源引起， 也可由于轴承 本身的结构特点及缺陷引起。而随着科学技术不断发展和工业化程度的不断提 高，机械设备精密程度、复杂程度及自动化程度不断提高， 凭个人的感观经验对 机械设备进行诊断己经远远不够，因此轴承的状态检测和故障诊断是十分必要 的，已经成为机械设备故障诊断技术的重要内容。 滚动轴承故障监测诊断方法有 很多种，它们各具特点，其中振动信号法应用最广泛。本次实验就是采用振动信 号法对滚动轴承故障实验平台的滚动轴承的故障信号进行分析。目录1绪论1.2轴承内圈故障特征频率 2.3时域无量纲参数分析2.3.1
3、 时域波形2.3.2 傅里叶变换运算分析故障 34通过自相关、互相关、功率谱运算分析故障 .44.1 自相关分析4.4.2 互相关运算分析故障.54.3 功率谱密度6.5 Haar小波分析 7.5.1 小波分解7.5.2 小波降噪9.故障诊断实用技术1绪论随着对滚动轴承的运动学、动力学的深入研究，对于轴承振动信号中的频率成分和轴承零件的几何尺寸及缺陷类型的关系有了比较清楚的了解，加之快速傅里叶变换技术的发展。开创了用频域分析方法来检测和诊断轴承故障的新领域。其中最具代表性的有对钢球共振频率的研究， 对轴承圈自由共振频率的研究。本文主要着重于对滚动轴承内圈磨损的故障研究，主要研究方法为傅里叶变换
4、，功率谱，自相关以及互相关，小波理论。滚动轴承在运行过程中可能会因为各种原因出现故障，如安装不当、异物入 侵、润滑不良、腐蚀和剥落等都会导致轴承出现故障。 安装不当会导致轴承不对 中，使得轴承在运行中，产生一种附加弯矩，给轴承增加附加载荷，形成附加激 励，引起几组强烈振动，严重时会导致转子严重磨损、轴弯曲、联轴器和轴承断 裂等严重后果。即使轴承安装正确，在长期的运行中，由于异物的入侵或则负荷 的作用下，接触面会出现不同程度的金属剥落、 裂痕等现象，进而导致旋转部件 与故障区域接触时产生强烈振动。本次实验主要针对潜在危害很大的裂痕故障信 号进行分析研究。滚动轴承在出现裂痕故障后，随着轴承的旋转，
5、由于旋转部件 与裂痕周期性的碰撞会产生周期性的冲击信号，且周期可以通过轴承结构计算得 出。图1.1所示为滚动轴承基本结构。图1.1滚动轴承基本结构d :滚动体直径D:轴承节径（滚动体所在圆的直径）R:内圈直径Ro :外圈直径a :接触角（滚动体受力方向与轴承径向平面的夹角）Z :滚动体个数w :内圈旋转角速度，正表示顺时针，负表示逆时针假设滚道与滚动体之间无相对滑动，且承受径向，轴向载荷时各部分无变形， 设当滚动轴承外圈固定，fi为内圈旋转频率，则可以推出 Z个滚动体在轴承不 同部件出现故障时，可分别计算出故障特征频率，即内圈故障特征频率：1Z dfbi= -(1 D cos: ) fi外圈故
6、障特征频率：,z dfbo= c(1 - c cos : ) fi2 D滚动体故障特征频率：- cosj d Dfi以上故障特征频率值均为理论计算值， 在实际应用中，由于滚动体可能出现 的相对滑动和受力不均导致的摇摆， 都会导致实际测量值与理论值有偏差， 在实 际实验过程中，我们需要将实际测量值与理论值进行对比分析。2轴承内圈故障特征频率内圈故障的特征频率： 1d .fi = n|Ni _N0|1 口 cos 上式中n为轴承滚动体的个数，Ni为轴承内圈的转动速度，No为轴承外圈的转动速度，d为滚动体的直径，D为滚动体中心所在圆的直径，邛为滚动体受 力方向与内外滚道垂直直线的夹角。进过计算后，将
7、以下各个数据代入上式，n =12, Ni =10, No =0, d=7.5mm,D =39.5mm, = =0,所计算的内圈故障特征频率 fi=76Hz。3时域无量纲参数分析3.1 时域波形内圈故障信号与正常信号采集，轴承外圈固定，转速为 600转/分，采样频率为10000HZ,采样时间为1.6384S共采集16384个点的数据。所采集到的数据 分别保存为故障信号 600 16384 10K.dat,和正常信号 ZC600 16384 10K.dat,下文中称正常信号，故障信号。信号通过以下程序导入MATLAB中，显示出时域由上图观察比较可以发现，正常时的振动信号要平稳很多，振幅较小，而故障
8、信号则振 动剧烈，振幅较大，呈现出周期性的特点。3.2 傅里叶变换运算分析故障傅里叶变换是数字信号处理领域一种很重要的算法。它表明：任何连续测量 的时序或者信号，都可以表示为不用频率的正弦波信号的无限叠加。 根据该原理, 傅里叶变换算法利用直接测量的原始信号，一类价方式来计算信号中不同正弦波信号的频率、振幅和相位。对故障信号和正常信号分别进行傅里叶变换，得出经过FFT处理后的数据,并根据香农采样定理，绘出5000Hz内的频率变化的振幅，如下图3.2和图3.3图3.2故障时频域波形图3.3正常时频域波形从上述图3.3中可以看出，正常状态下，轴承的低频段幅值突出。比较发现，在发生故障的情形下，轴承
9、的振动信号，具有很多噪音，不仅在低频段表现突出，而且在高频段谱峰突出比较明显，而且有多个谱峰，这主要是因为内圈有损伤时，在转动一周内，有时损伤点位于载荷区内，有时位于损伤区外；当损伤点位于载荷区内时，它与滚动体接触时产生脉冲力；当损伤点位于载荷区外时，如果不考虑运动时惯性力的作用，则不产生脉冲力；由此分析可以看出，内圈损伤引起的脉冲力的大小和方向受损点位置的影响，当产生的脉冲力较大时，能够引起高频共振，当产生的脉冲力较小时，引起中低频共振。4通过自相关、互相关、功率谱运算分析故障4.1 自相关分析自相关函数表达了同一过程不同时刻的相互依赖关系，而互相关函数表示不同过程的某一时刻的相互依赖关系。
10、这个在求解某些概率的时候是很有用的。自相关函数的定义和性质定义：描述信号x(t)在一个时刻的取值和另一个时刻取值之间的相似关系1 T. 1 ?.Rx( ) = - o x(t)x(t )dt Rx(k)为为 k式中：T、N 信号观测时间；t、k-时间问隔性质1)自相关函数Rx(T)是懒函数，即RxG)= Rx)；2)当 30时，Rx(0广金;当石加时，Rx。) < Rx(0);3)白噪声Rx(0)=max ,当7为时，氏国二04)周期信号的Rx，仍是周期信号，两者周期相同，但不反映相位信息对故障信号和正常信号分别进彳T自相关函数运算，得出图4.1和图4.2,横坐标为0至16384, 16
11、38532768个采样点数。纵坐标为相关程度图4.1正常信号自相关运算图4.2故障信号自相关运算正常状态的机器振动噪声是大量的、无序的、大小接近的随机冲击的结果， 其频较宽而均匀。机器运行状态不正常时，在随机噪声中将出现有规则的、 周期 性的脉冲，其大小要比随机冲击大的多。采用自相关分析方法：在振动噪声中查出隐藏的周期分量，特别是在故障发生初期，周期信号不明显、直观难以发现 的时候，依靠图3.4和图3.5中Rx(i)的幅值和波动的频率的比较可以看出正常信 号的自相关运算后其图形非常平滑，而故障信号自相关运算后，其中一部分幅值 以及频率被凸现出来，使我们更容易的发现异常，需要进一步的分析。4.2
12、 互相关运算分析故障互相关函数是描述随机信号 X(t),Y(t)在任意两个不同时刻t1, t2,的取值之 间的相关程度。自相关函数是描述随机信号 X(t)在任意两个不同时刻t1, t2,的 取值之间的相关程度。应用定义：互相关函数是描述两个信号之间的相似关系，可为_ , .1 T 1 ?Rxy( )0 x(t)y(t )dtRxy(k)XiykIN id性质：1) Rxy(9的峰值不一定在T = 0,峰值点偏离原点的距离表小两信号取得最大 相关程度的时移02)互相关函数是一非奇非偶的实函数，具有反对称性，RxyG)= Ryx (-T) o3)周期信号的Rxy(9也是同频率的周期信号，且保留了原
13、两信号的相位差信故而对正常信号和故障信号进行互相关运算，得出图4.3,横坐标为0至16384, 1638532768个采样点数。纵坐标为相关程度图4.3互相关运算故可以看出正常信号与故障信号的相关性越靠近 0点的时候越大，随着采样 点数的增加其相关性逐渐减小，因而可以判断出产生滚动轴承产生故障需要进一 步的分析。4.3 功率谱密度功率谱具有单位频率的平均功率量纲。 所以标准叫法是功率谱密度。通过功 率谱密度函数，可以看出随机信号的能量随着频率的分布情况。对于平稳随机过 程的样本函数一般不是绝对可积的，因此不能直接对它进行傅立叶分析。这里采 用利用相关函数的傅里叶变化来定义功率谱密度。所得时的功
14、率谱图形如图4.4, 正常时的功率谱图形故障如图4.5。图4.4故障时的功率谱图形正常信号自相关功率律频率HzJ图4.5正常时的功率谱图形通过比较可以看出，在产生故障时功率谱在低频，中频阶段都有更突出的幅 值产生，说明有更强的能量产生，而且整个频率段中能量比正常信号都要大得多。 故而可以知道故障滚动轴承的功率谱特点， 不仅在低频段表现突出，而且在高频 段谱峰突出比较明显，而且有多个谱峰，这主要是因为内圈有损伤时，在转动一 周内，有时损伤点位于载荷区内，有时位于损伤区外；当损伤点位于载荷区内时， 它与滚动体接触时产生脉冲力；当损伤点位于载荷区外时，如果不考虑运动时惯 性力的作用，则不产生脉冲力；
15、由此分析可以看出，内圈损伤引起的脉冲力的大 小和方向受损点位置的影响，当产生的脉冲力较大时，能够引起高频共振，当产 生的脉冲力较小时，引起中低频共振。5 Haar小波分析与Fourier变换相比，Haar小波变换是空间（时间）和频率的局部变换，因 而能有效地从信号中提取信息。通过伸缩和平移等运算功能可对函数或信号进行 多尺度的细化分析，解决了 Fourier变换不能解决的许多困难问题。小波分析是 时间一尺度分析和多分辨分析的一种新技术，它在信号分析方面的研究取得了有 科学意义和应用价值的成果。信号分析的主要目的是寻找一种简单有效的信号变换方法，使信号所包含的重要信息能显现出来。小波分析属于信号
16、时频分析的一种，它继承和发展了短时 傅立叶变换局部化的思想，同时又克服了窗口大小不随频率变化等缺点， 能够提 供一个随频率改变的时间一频率窗口，是进行信号时频分析和处理的理想工具。 它的主要特点是通过变换能够充分突出问题某些方面的特。5.1 小波分解对故障信号进行Hara小波分解ktz分解恭1)= % + WjT + + WD+ fa现在取j=11,分别计算出V11 Vl0 V 9 V 8 V 5 V 3绘制出图形如下图5.1绘制出 W10 W9 W8 W7 W6 W5 W4 W3 W2 Wi W0 图形如图 5.2 与图 5.31I °ij |/1 卜卜+小胖*护 j 则“什林 3
17、-T. 打 % g 1叫4电噂 十也用册一卜毗 十串 ,4+5L 枷*+*Ji iiIkiiII II知。1WM1500刘 002sg知。3sM好时罟D1111. 儿L JI.r.二4 L j LJu- - - -j. . , f下- prh> rrreui+1I|50C100016002000改 0雯Q111J111Ytl1_in *_" 4i1 嘀 i4d *,7-riL "： .* *-4 * *, ,'01-7 -11-飞川 m d一 itj, 1L . 12004fl0辆。陋1WW备口«» “第中科5息 JUl 丁一 11|I|1
18、00刘。涮400融硒* 口11111111115010015020025D300¥ QIIIIII*-1 "Iiu/XXR11 l_S_UIW_J1_,-_1r产、J , r1mi一(削40EOSO10012DUD图 5.2 系数图形 W10 W9 W8 W7 W6 W5图5.3系数图形 W4 W3 W2 Wi Wo5.2小波降噪信号降噪准则，第一，光滑性：在大部分情况下，降噪后的信号应该至少和 原信号具有同等的光滑性。第二，相似性：降噪后的信号和原信号的方差估计应 该是最坏情况下的最小值。小波分析用于降噪的过程分解过程：选用HARA小波，对信号进行3层小波分解。(2)作用
19、阈值过程：对分解得到的各层系数选择一个阈值，并对细节系数作 用软阈值处理。(3)重建过程：降噪处理后的系数通过小波重建恢复原始信号。小波变换中，对各层系数降噪所需的阈值一般是根据原信号的信号噪声比来 取的，在嘉定噪声为白噪声的情况下，一般应用原信号的小波分解的各层系数的 标准差来衡量。Matlab提供了 wnoisest来实现这个功能。在得到信号的噪声强度以后，我们就可以根据噪声强度来确定各层的阈值， 小波变换中的阈值c,l=wavedec(x1',5, 'haar'); %对信号进行层数为 5的尺度分解 sigma=wnoisest(c,l,1); %通过第一层细节系
20、数估算信号的噪声强度 thr1=wbmpen(c,l,sigma,2) %通过penalty策略确定降噪的阈值为 2 x11=wdencmp( 'gbl' ,c,l,'haar',5,thr1, 's',1);thr2,nkeep=wdcbm(c,l,2); %x12,cxd,lxd,perf0,perf12=wdencmp( 'lvd' ,c,l, 'haar',5,thr2, 's');通过以上2种降噪方式比较，得出图5.4。进彳TFFT变换得到图5.5。13图5.4 1原始故障信号时域图2通过
21、penalty策略确定降噪的阈值降噪3使用Birge-Massart策略确定降噪的阈值降噪图5.5对降噪后的故障信号进行 FFT运用相同的方式，对正常的信号进行降噪，并做FFT变换得图5.6,图5.7图5.6 1原始正常信号时域图2通过penalty策略确定降噪的阈值降噪3使用Birge-Massart策略确定降噪的阈值降噪图5.7对降噪后正常的信号进行 FFTXt FFT变换后的故障降噪信号和正常降噪信号取其幅值最大的10个点以及相应的频率。组成如下表格1和表格2幅值28.11826.55224.51924.00 021.53321.23519.90219.66 018.92118.9171
22、6516687611173214917812697表1故障降噪信号FFT处理后的频率值幅值23.77514.21710.8878.5157.6337.5617.3167.2186.9486.73916521645251021008934926364表2正常降噪信号FFT处理后的频率值由表1和表2中可发现频率87HZ与频率61HZ凸现出来，与特征频率76HZ 相符，故而可以推断出轴承的内圈出现了故障。6.结束语本文主要是对轴承内圈的故障信号进行处理， 并与正常信号进行比较。主要 用到的方法有：快速傅里叶变换（FFT）、功率谱、自相关及互相关运算、小波分 析，无论是哪一种方法，通过作图发现，故障信号与正常信号的差别是很明显的。 其中采用小波分析的方式，有着比较明显的优势。信号的小波分解由于其基函数 的时频局部化特性，可以很好地表示信号的非线性和冲击特征，而对于存在的高 频低幅噪声有极大的略化作用，能有效地区分信号中的突变部分和噪声， 从而实 现信号的降噪。结合小波降噪的优点，对信号进行降噪处理，成功实现了故障信 号特征频率提取，体现了小波在轴承故障诊断中的工程实用价值。 对上述实验结 果的分析，我们可以清楚的观察到滚动轴承的正常信号与故障信号在用不同分析 方法时所得出不同波形以及相应的结论。