[发明专利]一种滚动轴承故障诊断方法及系统有效
申请号: | 201910292988.6 | 申请日: | 2019-04-12 |
公开(公告)号: | CN110006652B | 公开(公告)日: | 2020-04-03 |
发明(设计)人: | 李继猛;王慧;于青文;王向东;张金凤;张云刚 | 申请(专利权)人: | 燕山大学 |
主分类号: | G01M13/045 | 分类号: | G01M13/045 |
代理公司: | 北京高沃律师事务所 11569 | 代理人: | 程华 |
地址: | 066000 河北*** | 国省代码: | 河北;13 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 滚动轴承 故障诊断 方法 系统 | ||
1.一种滚动轴承故障诊断方法,其特征在于,所述诊断方法包括如下步骤:
采集滚动轴承振动信号;
对所述滚动轴承振动信号进行离散傅里叶变换,得到滚动轴承振动信号的Fourier谱;
采用短时能量边界检测方法对所述Fourier谱进行自适应划分,得到多个连续区间;
获取每个连续区间的尺度函数、小波函数和经验小波变换的细节系数和逼近系数;
根据所述尺度函数、小波函数和经验小波变换的细节系数和逼近系数,对每个连续区间内的滚动轴承振动信号进行经验小波分解,得到多个调幅调频模态分量;
计算每个所述调幅调频模态分量的包络谱显著性指数;具体包括:获取第m个调幅调频模态分量的包络谱;m=1,2,…,M,M表示调幅调频模态分量的个数;获取所述包络谱的所有幅值极大值;将所有的幅值极大值按照从高到低的顺序排序,获取前K个幅值极大值;根据K个所述幅值极大值利用公式ESSIm=Fm1(ω)*Fm2(ω)*…*FmK(ω)计算m个调幅调频模态分量的显著性指数ESSIm,Fm1(ω)、Fm2(ω)和FmK(ω)分别表示第m个模态分量的第1个、第2个和第K个幅值极大值;
采用Teager能量算子解调包络谱显著性指数最大的调幅调频模态分量,获得最优解调包络谱;
利用所述最优解调包络谱进行滚动轴承的故障诊断。
2.根据权利要求1所述的滚动轴承故障诊断方法,其特征在于,所述采用短时能量边界检测方法对所述Fourier谱进行自适应划分,得到多个连续区间,具体包括:
将所述Fourier谱的频率范围归一化到[0,π],得到支撑区间;
将所述支撑区间均分为预设个数的子区间;
计算每个子区间的Fourier谱的能量值,获得能量值序列;
获取所述能量值序列的所有极小值点;
将每个所述极小值点所对应的子区间的边界频率的平均值作为新的频域分割区间边界;
根据所述新的频域分割区间边界对所述支撑区间进行区间划分,得到多个连续区间。
3.根据权利要求1所述的滚动轴承故障诊断方法,其特征在于,所述获取每个连续区间的尺度函数、小波函数和经验小波变换的细节系数和逼近系数,具体包括:
根据Little-Paley和Meyer原理,获取每个连续区间的经验尺度函数
其中,β(x)=x4(35-84x+70x2-20x3),ω为所述Fourier谱的频率,τm=γmωm,0<γm<1,ωm-1和ωm分别为第m个连续区间的左边界和右边界;
根据Little-Paley和Meyer原理,获取每个连续区间的小波函数:
其中,τm+1=γm+1ωm+1,0<γm+1<1,ωm和ωm+1分别为第m+1个连续区间的左边界和右边界;
根据所述经验尺度函数和小波函数,获取经验小波变换的细节系数和逼近系数:
其中,*表示取共轭,和分别为第1个连续区间的经验小波变换的逼近系数和第m个连续区间的经验小波变换的细节系数;ψm(τm-t)和分别是ψm(ω)和的Fourier逆变换,ψm(ω)和分别表示第m个连续区间的小波函数和第1个连续区间的经验尺度函数,F-1[·]表示傅里叶逆变换,f(t)表示滚动轴承振动信号,f(ω)表示滚动轴承振动信号的Fourier谱。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于燕山大学,未经燕山大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910292988.6/1.html,转载请声明来源钻瓜专利网。