[发明专利]基于自适应时变梳状滤波的变转速齿轮故障信号提取方法有效
申请号: | 202110324209.3 | 申请日: | 2021-03-26 |
公开(公告)号: | CN113074941B | 公开(公告)日: | 2022-08-30 |
发明(设计)人: | 陈向民;段萌;舒国强;张亢;卢绪祥;晋风华;李录平 | 申请(专利权)人: | 长沙理工大学 |
主分类号: | G01M13/045 | 分类号: | G01M13/045 |
代理公司: | 长沙市融智专利事务所(普通合伙) 43114 | 代理人: | 姚瑶 |
地址: | 410114 湖南省*** | 国省代码: | 湖南;43 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 自适应 时变梳状 滤波 转速 齿轮 故障 信号 提取 方法 | ||
1.一种基于自适应时变梳状滤波的变转速齿轮故障诊断方法,其特征在于:包括如下步骤:
步骤S1:采集变转速下的齿轮振动信号x(t)和转速信号fr(t),t=t1,t2,…,tN,N为信号采样长度;
步骤S2:用自适应时变梳状滤波器对齿轮振动信号进行时变梳状滤波,以获取时变梳状滤波信号,所述时变梳状滤波信号为提纯后的齿轮啮合振动成分及其调制边频成分;
步骤S3:根据转速信号对所述时变梳状滤波信号进行等角度重采样得到角域重采样信号;
步骤S4:对所述角域重采样信号进行频谱分析得到阶次谱,并根据所述阶次谱中啮合阶次周围的调制边频阶次诊断变转速下的齿轮故障。
2.根据权利要求1所述的方法,其特征在于:所述自适应时变梳状滤波器是根据转速信号而设计,所述自适应时变梳状滤波器的冲击响应及其频响函数H(ti,ω)分别表示如下:
式中,h(ti,ω)为ti时刻的冲击响应,H(ti,ω)为ti时刻的频响函数,fg(ti)表示ti时刻的啮合频率值;fr(ti)表示ti时刻的转频值;ωi表示ti时刻的滤波器半带宽,ω为角频率,K为滤波阶数,t表示时间。
3.根据权利要求1所述的方法,其特征在于:步骤S2中用自适应时变梳状滤波器对齿轮振动信号进行时变梳状滤波,以获取时变梳状滤波信号的过程如下:
S21:初始化滤波时刻ti=1;
S22:根据滤波时刻ti的自适应时变梳状滤波器的频响函数H(ti,ω)对齿轮振动信号x(t)进行滤波,得到ti时刻的时变梳状滤波信号点s(ti),表示如下:
式中,H[x(t)]表示齿轮振动信号x(t)的频谱信号,IFFT表示FFT逆变换,为y(t)中t=ti时刻的数据;
S23:判断是否存在i<N,若是,则i=i+1,并返回步骤S22,直至i=N,再将滤波得到的所有时刻的时变梳状滤波信号点按时间顺序排列,得到时变梳状滤波信号s(t),t=t1,t2,…,tN。
4.根据权利要求1所述的方法,其特征在于:步骤S3中根据转速信号对所述时变梳状滤波信号进行等角度重采样得到角域重采样信号的过程如下:
S31:采用三次多项式对转速信号fr(t)进行拟合,确定拟合系数a、b、c、d;
fr(t)=at3+bt2+ct+d
S32:基于设定的最大分析阶次Dmax计算出等角度重采样的角度间隔Δθ,如下:
S33:利用所述角度间隔Δθ计算出重采样后的数据长度L,如下:
式中,T为时域采样的总时间,存在T=tN-t1;
S34:根据拟合系数a、b、c、d,用如下公式计算出等角度重采样的键相时标Tn:
式中,T0为时域采样开始时间,即T0=t1;
S35:利用Langrange线性插值公式对时变梳状滤波信号s(t)进行插值,Langrange线性插值公式为:
式中,ω(Tn)为角域重采样信号。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于长沙理工大学,未经长沙理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110324209.3/1.html,转载请声明来源钻瓜专利网。