[发明专利]基于奇异值分量频域谱的共振带选择方法有效
申请号: | 201810067275.5 | 申请日: | 2018-01-24 |
公开(公告)号: | CN108287073B | 公开(公告)日: | 2019-12-27 |
发明(设计)人: | 马增强;张安;杨绍普;阮婉莹;陈明义;刘俊君;秦松岩;校美玲;张光跃;李响 | 申请(专利权)人: | 石家庄铁道大学 |
主分类号: | G01M13/045 | 分类号: | G01M13/045 |
代理公司: | 13128 石家庄轻拓知识产权代理事务所(普通合伙) | 代理人: | 王占华 |
地址: | 050000 河*** | 国省代码: | 河北;13 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 共振 矩阵 故障轴承 频域谱 采集 加速度传感器 奇异值分解 自适应选择 包络解调 分量中心 个数确定 故障诊断 后续轴承 频域变换 信号构造 信号频谱 频谱 包络 峭度 主频 轴承 带宽 测试 诊断 分析 | ||
1.一种基于奇异值分量频域谱的共振带选择方法,其特征在于包括如下步骤:
1)使用加速度传感器采集故障轴承的振动信息,将加速度传感器采集到的故障轴承信号进行FFT频域变换;
2)根据信号频谱的主频个数即峰值个数,将频谱划分成n个区域;
3)将采集到的故障轴承信号构造成行数为n+1的Hankel矩阵,对矩阵进行奇异值分解,得到n+1个奇异值分量;
确定Hankel矩阵行数方法如下:
根据信号频谱的主频个数即峰值的个数,将频谱划分成n个区域,Hankel矩阵行数即为n+1;对测得含有故障信息的振动信号x(i),i=1,2...N;基于相空间重构理论,构造出行数为n+1的Hankel矩阵,矩阵形式如下:
式中:N为故障轴承信号长度,n+1为矩阵的行数;
将由故障轴承信号构成的Hankel矩阵进行奇异值分解,分解公式如下:
其中:A是由信号构成的Hankel矩阵,U是m阶正交矩阵,V是n阶正交矩阵,∑=diag(σ1,σ2,…σr)是r阶对角阵,而σi是矩阵A的非零奇异值,0是零矩阵;
将式(1)中的零奇异值去除,将A的奇异值分解写成精简的分量形式,如下:
式中:ui、vi分别是U、V的第i个列向量;
即为分解后的得到的n+1奇异值分量,与Hankel矩阵的行数一致;每一个奇异值分量的频谱只反映了一部分的频带并且顺序靠前的分量频谱反映的是较低频率的频带,顺序靠后的分量反映的是频率较高的频带,各个奇异值分量频谱的频带按照奇异值分量顺序由低到高依次排列;
4)选择第n个分量进行包络解调分析,由包络谱峭度最大值来确定分量中心频率和带宽;
5)根据中心频率以及带宽计算具体的共振带数值。
2.如权利要求1所述的基于奇异值分量频域谱的共振带选择方法,其特征在于根据确定Hankel矩阵行数方法选择共振带所对应的奇异值分量的方法如下:
Hankel矩阵行数与奇异值分解后得到的奇异值分量个数是一致的,在n+1个奇异值分量中,第n个分量即是共振带所对应的分量;各奇异值分量频谱的频带按照奇异值分量顺序由低到高依次排列,且第n+1个分量是高频震动不明显的部分,不能反映共振的情况,所以选择第n个分量作为共振带所对应的分量。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于石家庄铁道大学,未经石家庄铁道大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810067275.5/1.html,转载请声明来源钻瓜专利网。