滚动轴承故障诊断(附MATLAB程序)VIP免费

第二组实验 轴承故障数据： Test2.mat 数据打开后应采用X105_DE_time 作为分析数据，其他可作为参考，转速1797rpm 轴承型号：6205-2RS JEM SKF, 深沟球轴承 采样频率：12k Hz 1 、确定轴承各项参数并计算各部件的故障特征频率 通过以上原始数据可知次轴承的参数为： 轴承转速r=1797r/min；滚珠个数n=9；滚动体直径d=7.938mm； 轴承节径D=39mm；：滚动体接触角α=0 由以上数据计算滚动轴承不同部件故障的特征频率为： 外圈故障频率f1=r/60 * 1/2 * n(1-d/D *cosα)=107.34Hz 内圈故障频率f2=r/60 * 1/2 * n(1+d/D *cosα)=162.21Hz 滚动体故障频率f3=r/60*1/2*D/d*[1-(d/D)^2* cos^2(α)]=70.53Hz 保持架外圈故障频率f4=r/60 * 1/2 * (1-d/D *cosα)=11.92Hz 2 .对轴承故障数据进行时域波形分析 将轴承数据Test2.mat 导入 MATLAB 中直接做 FFT 分析得到时域图如下： 并求得时域信号的各项特征： （1 ）有效值：0 .2 9 0 9 ； （2 ）峰值：1 .5 2 5 6 ； （3 ）峰值因子：5 .2 4 4 1 ；（4 ）峭度：5 .2 7 9 3 ； （5 ）脉冲因子：7 .2 8 8 4 ；（6 ）裕度因子：9 .1 0 8 3 ： 3 .包络谱分析 对信号做EMD 模态分解，分解得到的每一个IMF 信号分别和原信号做相关分析，找出相关系数较大的IMF 分量并对此IMF 分量进行Hilbert 变换。 signalEmpirical Mode Decompositionimf1imf2imf3imf4imf5imf6imf7imf8res. 由图中可以看出经过 EMD 分解后得到的９个IMF 分量和一个残余量。IMF 分量分别和原信号做相关分析后得出相关系数如下： EMD 分量 IMF1 IMF2 IMF3 IMF4 IMF5 相关系数 0.9596 0.1990 0.1096 0.0062 0.0230 EMD 分量 IMF6 IMF7 IMF8 IMF9 IMF10 相关系数 0.0032 0.0045 0.0055 0.0060 0.0062 由上表得：IMF1 的相关系数明显最大，所以选用 IMF1 做Hilbert包络谱分析。所得Hilbert 包络谱图如下： 对包络谱图中幅值较大区域局部放大得到下图 由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显，三个峰值频率由小到大排列分别为58.59Hz、105.5Hz、164.1Hz。把这三个频率数值和前文计算所得的理论值进行比较可知：频率值最大为164.1Hz 和内圈的故障理论计算特征频率f2=162.21Hz相近，说明此轴承的故障发生在轴承的内圈。 clc 程序1：原始信号时域分析及小波去噪处理 clear all z=i...