基于DMD降噪的滚动轴承故障诊断方法.pdf
《基于DMD降噪的滚动轴承故障诊断方法.pdf》由会员分享,可在线阅读,更多相关《基于DMD降噪的滚动轴承故障诊断方法.pdf(8页珍藏版)》请在咨信网上搜索。
1、第4 6卷第5期武汉科技大学学报V o l.4 6,N o.52 0 2 3年1 0月J o u r n a l o fW u h a nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yO c t.2 0 2 3 收稿日期:2 0 2 2-1 0-0 3 基金项目:国家自然科学基金资助项目(5 1 9 7 5 4 2 7).作者简介:涂福泉(1 9 7 0-),男,武汉科技大学教授,博士.E-m a i l:t u f u q u a n w u s t.e d u.c nD O I:1 0.3 9 6 9/j.i s s
2、n.1 6 7 4-3 6 4 4.2 0 2 3.0 5.0 0 8基于DMD降噪的滚动轴承故障诊断方法涂福泉,杨家瑜,陈 超,罗迎九,吴维崧(武汉科技大学冶金装备及其控制教育部重点实验室,湖北 武汉,4 3 0 0 8 1)摘要:针对轴承振动信号难以剔除噪声的问题,提出一种将H i l b e r t变换、动力学模态分解(DMD)和奇异值差分谱相结合的滚动轴承故障诊断方法。首先将原始信号进行H i l b e r t变换得到包络信号,由包络信号构造H a n k e l矩阵进行动力学模态分解,利用奇异值差分谱确定合适的截断秩后进行信号重构,最后通过频谱分析来提取故障特征。采用该方法对滚动轴
3、承故障仿真信号和实验数据进行分析,结果表明其降噪效果显著,能有效获取轴承故障特征频率。关键词:故障诊断;滚动轴承;动力学模态分解;H i l b e r t变换;奇异值差分谱;信号降噪中图分类号:TH 1 3 3.3 3 文献标志码:A 文章编号:1 6 7 4-3 6 4 4(2 0 2 3)0 5-0 3 7 6-0 8 滚动轴承是机械传动中的重要组成部件,广泛应用于各种机械设备,如风电机组、高速列车、航空发动机等。由于工作环境的复杂性,滚动轴承容易出现点蚀、断裂、磨损等故障,从而影响设备运行的稳定性和可靠性,因此对滚动轴承进行故障诊断具有重要意义。在设备故障诊断时,通常是利用传感器收集的
4、振动信号进行分析。这些振动信号中有大量的噪声,导致其故障特征不明显,难以直接应用,需要将原始信号进行降噪处理,再进行下一步的分析1。经典的信号降噪和特征提取方法包括小波变换(WT)、经验模态分解(EMD)、变分模态分解(VMD)、奇异值分解(S V D)等,然而这些方法都有一些不足之处。例如,WT需要选择合适的母小波及分解层数;EMD公式缺乏严格的理论推导,并且存在模态混叠和端点效应等问题;VMD的参数选择对其分解效果影响较大;S V D对强噪声背景下的信号特征提取效果欠佳。因此,研究人员以上述方法为基础提出了一系列的改进措施。Y a n g等2针对传统的采用单一准则优化的小波参数无法完全描述
5、故障特征的问题,设计了一种多目标优化小波滤波器用于提取轴承故障特征,在强噪声背景下的应用效果较好。滚动轴承早期故障信号存在大量噪声使得故障特征提取十分困难,为此曹玲玲等3提出一种将互补集合经验模态分解(C E EMD)和改进小波阈值算法相结合的方法,可以提高故障信号的信噪比,降噪效果明显。Y i n等4提出一种滚动轴承弱故障特征提取方法,能选择出具有故障特征的内模函数,并且通过新构造的自适应阈值进行信号去噪。赵玮5采用基于VMD和奇异值差分谱的故障特征提取方法,可 以 有 效 降 低 噪 声 的 影 响。鉴 于 传 统VMD方法对高速列车轮对轴承故障特征提取不够准确,李翠省等6设计一种结合集成
6、经验模态分解(E EMD)和改进参数自适应VMD的方法,其在强噪声干扰下仍具有优良的故障特征提取效果。李华等7为了提升S V D方法的降噪性能,提出一种基于相关奇异值比的S V D方法并成功应用于轴承故障诊断。C h e n g等8提出自适应周期模态分解方法,以解决现有滚动轴承故障诊断方法对周期性冲击故障识别能力差的问题。动力学模态分解(d y n a m i cm o d ed e c o m p o s i-t i o n,DMD)是基于K o o p m a n算子的无方程数据驱动分析方法,它不仅有S V D算法的优点,还考虑了机械系统演化过程中的时空特性,具有时空耦合建模的独特优势,在
7、信号降噪和特征提取方面已得到有效应用9。DMD本质上是一种融合奇异值分解和模态分解理论的空间降维方法,因此它也不可避免地继承了两者的缺点,例如在截断秩的选取和含调制成分信号的处理方面存在困难。针对上述问题,本文提出一种H i l b e r t变换和DMD相结合并通过奇异值差分谱确定截断秩的信号降噪和特征提取方法,用于滚动轴承故障诊断。首先使用H i l b e r t变换对原始信号进行包2 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法络解调,消除原信号中调制成分对后续分析的影响,然后进行动力学模态分解,通过奇异值差分谱选择合适的截断秩进行信号重构,从而减少噪声的干扰,以
8、达到较好的特征提取效果。本文最后利用仿真信号和美国凯斯西储大学公开的轴承故障数据集进行方法验证。1 基本理论1.1 H i l b e r t变换对一个连续信号x(t)进行H i l b e r t变换可表示为1 0 x(t)=Hx(t)=1n-x()t-d=x(t)1t(1)式中:Hx(t)表示时域中的H i l b e r t变换,为卷积算子。原信号的解析信号为y(t)=x(t)+jx(t)(2)从解析信号中可以得到原信号x(t)的包络信号为B(t)=x2(t)+x2(t)(3)1.2 DMD算法假设两个连续序列X、X(用矩阵形式表示)满足与最优线性算子的映射关系,即X=A X,DMD算法
9、利用矩阵A的低秩表达式来捕捉连续序列潜在的动态特性1 1-1 2。由于DMD最初是应用于流体力学领域,所收集的信号通常是二维矩阵形式,而轴承振动信号是一维的,因此需要将一维振动信号转换为矩阵形 式,这 里 构 造H a n k e l矩 阵 以 满 足 上 述 要求1 3。对于振动信号X(t)=x(1),x(2),x(3),x(N),其H a n k e l矩阵形式为X=x(1)x(n-1)x(m)x(N-1),X=x(2)x(n)x(m+1)x(N)(4)式中:m=N/2,1nN且n=N-m+1。DMD算法步骤为:1)对X进行奇异值分解X=U V*(5)式中:U、V为酉矩阵,为奇异值对角矩阵
10、。U的每个列 向 量 相 互 正 交,是X的 本 征 正 交 分 解(P O D)模态。取U的前r阶来截断数据,通过保留主要模态来降维。2)利用X的伪逆来求AA=X XT=X V-1U*(6)通常情况下,A包含的数据量很大,故通过合适的投影矩阵将其降维为r阶的矩阵A-:A-=U*rA Ur=U*rX V-1(7)3)计算矩阵A-的特征值和特征向量A-W=W(8)式中:W是矩阵A-的特征向量矩阵,是由矩阵A-的特征值i组成的对角矩阵。4)计算重构矩阵定义矩阵:=X V-1W(9)中每一个列向量i称为DMD的模态,并且每个模态对应着一个特征值i。最后得到由若干个模态组成的DMD重构矩阵:XDMD=
11、ri=1ie x p(wit)bi(1 0)式中:wi=l n(i)/t,bi是每个模态的幅值。1.3 奇异值差分谱奇异值差分谱可以表征奇异值序列的变化,差分谱计算公式如下1 4:di=i-i+1ni=1(i-i+1),i=1,2,n(1 1)式中:i为奇异值。所有的di构成奇异值差分谱,差分谱的峰值与信号特征相关,峰值越大,特征越明显,当峰值降低到一定程度便可截断。2 本文方法本文方法的基本流程如图1所示。首先将原始信号进行H i l b e r t变换得到包络信号,再根据包络信号构造H a n k e l矩阵进行动力学模态分解,然后通过奇异值差分谱确定DMD的截断秩,使用该截断秩进行重构得
12、到重构矩阵,其第一行图1 本文方法的基本流程F i g.1F l o w c h a r t o f t h ep r o p o s e dm e t h o d773武汉科技大学学报2 0 2 3年第5期为重构信号,根据重构信号进行频谱分析,提取出准确的信号特征,从而实现滚动轴承故障的精确诊断。3 仿真信号分析下面采用滚动轴承内圈故障仿真信号来验证本文方法的有效性。仿真信号表达式如下1 5:x(t)=s(t)+n(t)=iA(t)h(t-i T)+n(t)h(t)=e x p(-C t)c o s(2 fnt)A(t)=1+A0c o s(2 frt)(1 2)式中:fr为转频,fn为轴承
13、的共振频率,T为故障周期,C为衰减系数,A0为幅值,n(t)为高斯白噪声。设定采样点数N=40 9 6,高斯白噪声的信噪比为-5d B,仿真信号其余参数的具体取值如表1所示,其中fs为系统的采样频率。表1 滚动轴承内圈故障仿真信号的参数T a b l e 1P a r a m e t e r s o f t h e s i m u l a t e d s i g n a l o f r o l l i n gb e a r i n gw i t h i n n e r-r a c e f a u l t参数A0fr/H zCT/sfs/k H zfn/k H z取值0.3 52 57 5 01
14、/1 2 01 23 图2为仿真信号的时域波形图。由图2(b)可见,时域信号的故障特征已经被噪声淹没,从中无法获得轴承故障频率。为了分析信号的故障特征,对加噪声后的仿真信号进行包络谱分析,结果如图3所示。在图3中可以找到信号的特征频率和部分多倍频,但是包含大量的噪声信息。(a)未加噪声的仿真信号(b)加噪声的仿真信号图2 仿真信号的时域波形图F i g.2T i m e-d o m a i nw a v e f o r m so f t h e s i m u l a t e ds i g n a l s图3 加噪声仿真信号的包络谱F i g.3E n v e l o p e s p e c
15、t r u mo f t h e s i m u l a t e ds i g n a lw i t hn o i s e 为了更加有效地提取轴承的故障特征,使用本文方法进行信号分析。首先将原始信号通过H i l b e r t变换构造解析信号,对解析信号取模得到包络信号,如图4所示。8732 0 2 3年第5期涂福泉,等:基于DMD降噪的滚动轴承故障诊断方法图4 本文方法得到的包络信号F i g.4E n v e l o p e s i g n a l b y t h ep r o p o s e dm e t h o d 根据包络信号构造H a n k e l矩阵,进行动力学模态分解,通过
16、奇异值差分谱确定截断秩的阶数,如图5所示。需要在奇异值分布曲线上选取一个合适的截断秩阶数,但是仅根据图5(a)难以确定最优点,因此可以根据奇异值差分谱进行分析。在理想状态下,差分谱峰值的大小与信号特征相关,峰值越大,特征越明显,当差分谱峰值降低到一定程度而难以区分特征信息时便可截断。如图5(b)所示,通过奇异值差分谱分析,设定奇异值差分阈值为0.0 0 1,取最后一个峰值大于阈值的差分阶数作为截断秩阶数,最终确定的截断秩阶数为3 1。(a)奇异值分布曲线(b)奇异值差分谱图5 仿真信号的奇异值分布及奇异值差分谱F i g.5S i n g u l a rv a l u ed i s t r i
17、 b u t i o na n ds i n g u l a rv a l u ed i f f e r e n c e s p e c t r u mo f t h e s i m u l a t e ds i g n a l 根据设定的截断秩阶数,可以得到3 1个精确DMD模态,然后构造DMD重构矩阵,选取其中第一行为重构信号。通过DMD重构后的信号如图6(a)所示,从中可以观察到故障特征。利用快速傅里叶变换(F F T)对重构信号进行频谱分析,结果如图6(b)所示,可以观察到特征频率及其多倍频,谱线突出,特征频率附近存在一对边频,4倍频附近存在一个上边频。上述结果表明,本文方法能有效降低
18、轴承振动信号中的噪声干扰,可提取出轴承故障特征频率。(a)时域波形(b)频率谱图6 本文方法对仿真信号的分析结果F i g.6A n a l y s i s r e s u l t so f t h e s i m u l a t e ds i g n a l b y t h ep r o p o s e dm e t h o d 图7所示为对比方法(包括S V D、自适应DMD、VMD、E EMD)对仿真 信 号 的 分 析 结 果。S V D重构信号的包络谱中特征频率和部分多倍频谱线明显,但是存在部分噪声;自适应DMD重973武汉科技大学学报2 0 2 3年第5期构信号的包络谱中特征频率和
19、部分多倍频谱线突出,但高于5倍频的特征不明显,部分特征信息丢失;选 取VMD重 构 信 号 中 合 适 的 模 态 分 量(I MF 3)进行分析,其包络谱中故障频率、2倍频和3倍频的谱线较明显,但伴随有部分噪声;选取E EMD重构信号中合适的模态分量(I MF 1)进行分析,其包络谱中故障频率、部分多倍频的谱线较明显,但噪声较多。(a)S V D重构信号的包络谱(b)自适应DMD重构信号的包络谱(c)VMD重构信号I MF 3分量的包络谱(d)E EMD重构信号I MF 1分量的包络谱图7 其他方法对仿真信号的分析结果F i g.7A n a l y s i s r e s u l t so
20、 f t h e s i m u l a t e ds i g n a l b yo t h e rm e t h o d s4 实验验证为了进一步验证本文方法在实际应用中的有效性,以美国凯斯西储大学公开的轴承故障数据集中的内圈故障数据为案例进行分析。轴承故障直径为0.0 0 7英寸(0.0 1 7 7 8c m),转速为1 7 9 7r/m i n,采样频率为1 2k H z。通过理论计算,轴承内圈故障特征频率为1 6 2.1 9H z。由轴承内圈故障信号的时域波形(图8)可知,故障特征被噪声淹没,无法提取出有效的特征信息。该信号的包络谱分析结果见图9,从包络谱中可以提取轴承故障特征,但是还
21、存在噪声干扰。为了更好地提取故障特征信息,使用本文方法对原始信号进行降噪处理。首先将原始信号通过H i l b e r t变换得到包络信号,如图1 0所示。根据包络信号构造H a n k e l矩阵后进行动力学模态分解,并通过奇异值差分谱确定截断秩的阶数。轴承内圈故障信号的奇异值分布及奇异值差分谱如图1 1所示,据此选取合适的截断秩阶数。从图1 1(a)来看,应选取曲线变化趋于稳定的一个点。设定奇异值差分阈值为0.0 0 1,结合图1 1(b),确定截断秩阶数为1 0 5。图8 轴承内圈故障信号的时域波形图F i g.8T i m e-d o m a i nw a v e f o r mo f
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 DMD 滚动轴承 故障诊断 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。