机械故障大作业.doc

机械故障诊断课程设计

设计题目：基于小波分析的轴承故障诊断
学 院： 机械工程系
专 业： 机械制造设计及其自动化
班 级：
学 号：
姓 名：
指导老师： 李奕璠
完成日期：

摘 要
利用小波对轴承的振动信号进行小波分解，采用编程快速地在计算机上实现基于小波分析的电机滚动轴承故障诊断，通过基于小波分解系数对含有故障特征频率的第一层细节信号进行小波重构并提取其包络谱，从中检测出故障特征频率，据此判断故障类型。
关键词：滚动轴承故障诊断 小波分析 包络谱
目 录

第一章 绪 论1
第二章 基于小波变换的滚动轴承故障诊断2
2.1 滚动轴承故障的特征频率2
2.2 基于小波分析的轴承故障诊断算法2
第三章 轴承故障诊断的程序分析5
3.1 确定轴承各项参数并计算各部件的故障特征频率5
3.2 故障诊断的结果分析和结论5
3.2.1 第一组数据分析5
3.2.2 第二组数据分析7
3.2.3 第三组数据分析8
3.2.4 第四组数据分析9
参考文献12
附录13

第一章 绪 论
滚动轴承在机械设备中使用非常广泛，其工作状态直接影响整个设备的运行品质,对滚动轴承进行状态监测和故障诊断，能够避免重大事故的发生，获得较大的经济和社会效益。随着生产的需要，对轴承故障的检测方法也越来越多，其中，运用比较广发的集中方法是、功率谱、倒谱、小波分析、人工神经网络、希尔伯特-黄变换、双谱。
小波变换是一种时频分析方法，进行多分辨率分析，即，将信号分解成若干层次的细节信号及概貌信号。对轴承振动信号进行小波变换，提取其中具有故障特征的细节信号进行重构；对重构信号做包络谱分析，从中检测出轴承的故障特征频率，据此判断故障类型。利用软件编程快速地实现了基于小波变换分析的滚动轴承故障判断。
第二章 基于小波变换的滚动轴承故障诊断
2.1 滚动轴承故障的特征频率
滚动轴承由外圈、内圈、滚动体和保持架组成，工作时外圈和轴承座或机壳相连接、固定或相对固定，内圈和机械传动轴相连接，随轴一起转动。当滚动轴承表面发生损伤故障，如内圈、滚动体或外圈出现点蚀、裂纹或剥落等，根据不同的损伤部位，按以下公式分别计算轴承故障的特征频率，如下所示：
外圈故障频率：

内圈故障频率 ：

滚动体故障频率 ：

其中，r为转速，n为滚珠个数；d为滚动体直径，D为轴承节径，α为滚动体接触角。
2.2 基于小波分析的轴承故障诊断算法
一般采用加速度传感器在轴承座上检测滚动轴承的振动信号。若周成表面出现局部损伤，在受载运转时轴承其他零件会周期地撞击损伤点产生低频的冲击信号，其频率即故障频率，但检测该频率主要会遇到2个问题：（a）冲击信号的宽频带性质会激起轴承结构及传感器本身在各自固有频率上发生谐振，故轴承振动信号中还含有故障特征频率的高次谐波分量。（b）由于轴承间隙的存在，冲击信号还要对轴承的高频固有振动信号进行调制。导致固有频率被其它振动所干扰而无法直接通过频谱分析检测出故障特征频率。本文分别采用小波分析和包络谱分析解决上述两个问题。
（1）小波分析提取含故障特征频率的细节信号。小波是一种均值为零，很快衰减的瞬间振荡函数，小波分析是一种时频分析方法，他利用一系列伸缩和平移的小波函数对信号进行展开，该过程等效于用一系列不同频带的高通和低频滤波器将信号分解成若干层次的高频细节信号及低频概貌信号（即信号的主体轮廓），可对信号进行多分辨率分屏，被誉为“数学显微镜”。
小波分析算法的步骤包括分解和重构，为在计算机上实现小波分析，根据二进离散小波变换的快速算法算法进行计算，小波变换公式如下：

(3)
式中

不同类型的小波，如小波、小波、墨西哥草帽小波等，滤波系数和均不相同。序列是的二进离散小波变换。利用小波分解系数重构原信号的公式为：
(4)
根据公式(3)(4)对轴承振动信号进行小波分解和重构可获得其各层概貌信号及细节信号，其中幅度最大的细节信号中包含轴承故障的特征频率。
（2）变换包络谱检测轴承的故障特征频率。含有轴承故障特征频率的细节信号是种调幅信号，它是故障信号对轴承的高频固有振动进行幅度调制形成，设其为式（5）。
（5）
式中

变换可对调幅信号进行包络解调，就是从中提取。信号的变换是和的卷积（符号为“*”），公式见式（6）。
（6）
对做傅里叶变换得式（7）
（7）
故的变换可看成是通过一个幅度为1的全通滤波器输出，其频率成分做相移，负频率成分做相移。则调幅信号的变换为式（8）。
（8）
设的解析信号为：,则有式（9）。
（9）
因此可利用变换提取的包络，即故障信号，再用傅里叶变换对其进行功率分析，功率谱中幅度最大处的频率即故障特征频率。

第三章 轴承故障诊断的程序分析
3.1 确定轴承各项参数并计算各部件的故障特征频率
由轴承型号为 6205-2 ，转速1750 可知：
滚珠个数9；滚动体直径7.938；
轴承节径39；滚动体接触角α=0；；
由以上数据计算滚动轴承不同部件故障的特征频率为：
外圈故障频率：

内圈故障频率 ：

滚动体故障频率 ：

3.2 故障诊断的结果分析和结论
3.2.1 第一组数据分析
首先对轴承振动信号文件中数组y第一列运行程序，结果如下图所示：

图1 振动加速度图和功率谱

图2 小波分解和重构图

图3 包络谱
图1为该信号的功率谱，显然从功率谱中难以检测故障频率，图2为对振动信号做3级小波分解和重构所得第1-3层细节信号d13和第三层概貌信号c3，对整体幅度较大的细节信号d1做包络谱（见图3），其幅度最大处的频率27.83为轴承内外圈相对旋转频率，第二大处频率为80.37同滚动体故障特征频率相接近，因此轴承故障类型为滚动体故障；
3.2.2 第二组数据分析
再用同样的方法对振动信号文件中数组y第二列进行处理，得的包络谱如图6所示：

图4 振动加速度图和功率谱

图5 小波分解和重构图

图6 包络谱
图4为该信号的功率谱，图5为对振动信号做3级小波分解和重构所得第1-3层细节信号d13和第三层概貌信号c3，对整体幅度较大的细节信号d1做包络谱（见图6），没有和轴承故障特征频率相同的频率，其幅度最大处的频率1796，第二大处频率为1036，第三大频率为29.3,同内外圈相对特征频率相接近，因此该列轴承数据正常。
3.2.3 第三组数据分析
对振动信号文件中数组y第三列进行处理，得的包络谱如图9所示：

图7 振动加速度图和功率谱

图8 小波分解和重构图

图9 包络谱
图7为该信号的功率谱，图8为对振动信号做3级小波分解和重构所得第1-3层细节信号d13和第三层概貌信号c3，对整体幅度较大的细节信号d1做包络谱（见图9），其幅度最大处的频率104.7，同轴承外圈理论故障特征频率相接近，因此该列轴承数据表示外圈故障。
3.2.4 第四组数据分析
对振动信号文件中数组y第四列进行处理，得的包络谱如图12所示：

图10 振动加速度图和功率谱

图11 小波分解和重构图

图12 包络谱
图10为该信号的功率谱，图11为对振动信号做3级小波分解和重构所得第1-3层细节信号d13和第三层概貌信号c3，对整体幅度较大的细节信号d1做包络谱（见图12），其幅度最大处的频率157.5，同内圈理论故障特征频率相接近，说明此轴承的故障发生在轴承的内圈。

参考文献
[1] 褚福磊,彭志科,冯志鹏,李志农.机械故障诊断中的现代信号处理方法[M].北京：科学出版社.2009:159-160
[2] 葛哲学,陈仲生时频分析技术及其应用[M].北京：人民邮电出版社.2006:164
[3] 杨国安.机械设备故障诊断实用技术[M].北京：中国石化出版社.2007

附录
程序：以第一列数据为例
%清除内存
%调入故障轴承振动信号的数据文件，存入
;
12000设置取样频率，取样点数N
65536;
1;
(1,1);
0:(1)*; %时间轴
(())(,1); %对归一化
(211); %绘制波形
();
('时间');
('振动加速度');
16384;
(); %对做功率谱
(212); %绘制功率谱
((02-1)*(12));
('频率');
('功率谱');
[](,4,'2'); %利用2对进行3级小波分解
c3('a','2',3)重构第1-3层细节d13和第3层盖帽c3
d3('d','2',3);
d2('d','2',2);
d1('d','2',1);
;
(414)(3)('c3')绘制c3
(413)(3)('d3')绘制d3
(412)(2)('d2')绘制d2
(411)(1)('d1')绘制d1
(d1)对d1进行变换，得y
(y)
();对去均值（目的是去除幅度较大的直流分量）
()的功率谱为P
;
((02-1)*(12))绘出d1的包络谱
('频率');
(12)故障频率f1为包络谱中幅度最大处的频率
[1](P);
f11*1;