请教FFT结果意义

skyreason

小弟最近研究滚动轴承状态监测,由于对FFT,小波,HHT刚接触不久
很多问题不太懂!希望各位师兄师姐帮忙!

滚动轴承转动频率23HZ
滚动轴承：选用6208型单排滚珠式轴承，外圈直径：80mm，内圈直径：40mm，滚珠数：10


分析三种情况,分别为轴承正常,内圈故障和外圈故障
x为轴承正常情况时的数据　y为轴承内圈故障时的数据　z为轴承外圈故障时的数据
各１万个点
FFT程序如下:
load x.txt;load y.txt;load z.txt;
x=x(1:1024);y=y(1:1024);z=z(1:1024);
k=fft(x,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=1000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(311);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号正常时频谱图');
k=fft(y,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=1000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(312);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号内圈故障时频谱图');
k=fft(z,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=1000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(313);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号外圈故障时频谱图');


分析图形如下图：




不太明白分析之后说明了什么问题，是不是正常情况下２３ＨＺ附近频率幅度大，其他两个小，所以损坏明显，还是其他什么？
x.txt (101.71 KB, 下载次数: 141)

2007-5-15 19:02 上传

点击文件名下载附件
Ｘ

y.txt (102.67 KB, 下载次数: 84)

2007-5-15 19:02 上传

点击文件名下载附件
Ｙ

z.txt (102.56 KB, 下载次数: 85)

2007-5-15 19:02 上传

点击文件名下载附件
Ｚ

另外做小波分析和ＨＨＴ的时候总提示程序有问题

希望有高手能帮忙分析一下，并解释一下图形
谢谢！

yangzj

测的是加速度还是速度？

skyreason

原帖由 yangzj 于 2007-5-15 19:06 发表
测的是加速度还是速度？

测的加速度!!

yangzj

先把内外环故障特征频率算出来

skyreason

原帖由 yangzj 于 2007-5-15 19:26 发表
先把内外环故障特征频率算出来

内:133HZ
外:96HZ

yangzj

个人认为这个１０００的采样频率太低．轴承帮障一般会激起内外环或滚动体的固有频率，而且为故障特征频率所调制，这个固有频率通常都很高，要采用高的采样频率才能采到

skyreason

原帖由 yangzj 于 2007-5-15 19:45 发表
个人认为这个１０００的采样频率太低．轴承帮障一般会激起内外环或滚动体的固有频率，而且为故障特征频率所调制，这个固有频率通常都很高，要采用高的采样频率才能采到

不好意思!可能没写清楚
采样频率是10000

yangzj

原帖由 skyreason 于 2007-5-15 19:48 发表


不好意思!可能没写清楚
采样频率是10000

呵呵，晕了，你程序里写的是１０００

skyreason

原帖由 yangzj 于 2007-5-15 19:45 发表
个人认为这个１０００的采样频率太低．轴承帮障一般会激起内外环或滚动体的固有频率，而且为故障特征频率所调制，这个固有频率通常都很高，要采用高的采样频率才能采到

师兄有时间的话能否用你的方法做一下
数据这里都有
我想比较一下!
谢谢!:loveliness:

yangzj

轴的转频是２３？
再提供下轴承的滚动体直径，轴承节径和接触角？

yangzj

还有你这个是一个纯粹的轴承结构还是一个齿轮箱或是别的机构

skyreason

不好意思!少写了个0
load x.txt;load y.txt;load z.txt;
x=x(1:1024);y=y(1:1024);z=z(1:1024);
k=fft(x,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=10000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(311);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号正常时频谱图');
k=fft(y,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=10000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(312);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号内圈故障时频谱图');
k=fft(z,1024);%对原始信号进行离散傅立叶变换
p=k.*conj(k)/1024;%计算功率谱密度
ff=10000*(0:511)/1024;%计算变换后不同点所对应的频谱值
subplot(313);plot(ff,p(1:512));
Xlabel('功率谱密度');
Ylabel('频率');
title('信号外圈故障时频谱图');

skyreason

实验说明：
1、本实验以滚动轴承为研究对象。故障轴承置于轴承座上，轴承座可以拆卸，以便更换不同故障类型的轴承。装置中的减速箱是另一课题的研究项目，在此不作介绍。本实验已能确保其处于正常工作状态，不会对轴承信号的测取造成影响。
2、装置加工时确保各元件主轴对中在同一轴线上；同时在底座下方垫上一层厚为10mm的泡沫橡胶板，减少实验台工作时与地板之间的振动，这样使整台装置处于较好的工作状态。
3、为了采集到较好的振动信号，本实验系统中将轴承座顶部铣平，以便传感器能通过底盘的磁铁吸附于其上，测量径向振动，拾取轴承振动信号。这样减少了中间界面，能够尽可能多地获得轴承外圈本身的振动信号。
4、本实验将电机与变频器通过导线连接，变频器通电为电机提供动力，同时控制其转速，使电机带动轴承运转。
5、装置中的磁粉制动器为我们的实验增加一定的负载，防止电机空转。
本实验是对单个滚动轴承进行故障诊断。


是单纯的滚动轴承，采样频率10KHz,主轴转速23Hz
接触角０，节径６０

skyreason

滚动体半径是５直径是１０ＭＭ

skyreason

原帖由 yangzj 于 2007-5-15 19:59 发表
还有你这个是一个纯粹的轴承结构还是一个齿轮箱或是别的机构

师兄在不在?
帮帮忙嘛!