基于eemd的滚动轴承振动故障特征提取与诊断研究

基于eemd的滚动轴承振动故障特征提取与诊断研究

ID:34870206

大小:12.73 MB

页数:75页

时间:2019-03-12

上传者:U-56225
基于eemd的滚动轴承振动故障特征提取与诊断研究_第1页
基于eemd的滚动轴承振动故障特征提取与诊断研究_第2页
基于eemd的滚动轴承振动故障特征提取与诊断研究_第3页
基于eemd的滚动轴承振动故障特征提取与诊断研究_第4页
基于eemd的滚动轴承振动故障特征提取与诊断研究_第5页
资源描述:

《基于eemd的滚动轴承振动故障特征提取与诊断研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

’,,,一、一_t巧撕'…?W飞>-'‘々一:;"/%;>V.4。.:‘,与//^斗餐''’‘;"、6屯.'';.‘■‘^?穿..巧;冉:-,-、^^巧一.声.;瑞妄巧卓,;'‘户:沪,巧^^詩.;:皆羊據旁證管邁读據紫:..鱗5編蚕巧蔡^觀無菊穂義;編黨璋載巧窠:品9—―’'‘鮮^-如或/A3臾:'"'一-’心’;二V、硕±学位论文%々典、基于EEMD的滚动轴承振动故障特征提取与V;诊断研究;.‘‘‘.y’*一.U:^少ResearchonExtractionandDiagnosisofFauUFeatureof气占:VVibrationofR幻HinBeari打BasedonEEMDVgg'点卢辦',■■--^■*4"^■A一,六w一,?;|;r女'、,.;^产含、产户f.巧,巧片於早、‘’—‘''一命、‘、-‘養"’’心t.,rr:叩二;将严';>■"-^-->...;;;;;::.■.'由苗吗化班诚:W进批戊串玻V姆續巧抑?’种,兴心'成心巧減^在乃户:敵兴於^/巧,:?'議画瞧i^i樣备售:?li^?|iiii‘’-'..^心■.;w嘴?韻鑽類P:辩游疑奎爲郑誘賊難送'‘:海:’.'解鑽落蕉讀群护轉專黄孩羯键吊矣片V.、.:八..>,^201531議繫顯變議^^:年月斯r\:*、‘?*记1:,:;巧巧;满钱,>式淵、.'"‘-i’W.'一.、告乂L*r去.'V、\一.\,yJ;一身, 国内图书分类号:TH17学校代码:10079;国际图书分类号:621密级公开硕i学位论文基于EEMD的滚动轴承振动故障特征提取与诊断硏究硕±研究生;刘觉晓导师;杜冬梅教授申请学位;工学硕±学科:机械工程专业;机械电子工程所在学院;能源动为与机械工程学院答辩日期:2015年3月授予学位单位:华北电力大学 ClassifiedIndex:TH17U:.D.C621ThesisfortheMasterDegreeResearchonExtractionandDianosisofFaultFeatureofgVibrationofRollinBearinBasedonEEMDggCandidate:LiuJuexiaoSuervisor:ProfDuDonmeip.gAcademicDegreeAppliedfor:MasterofEngineeringSpeciality:MechanicalEngineeringSchoolofEnerPowerandMechanical。L.gy,School:r.EnmeennggDateofDefence:March2015,Deree-Conferrin-Instggitu杜on:NorthChinaElectricPowerUnivershy 1华北电力大学硕±学位论文原创性声明本人郑重声明:此处所提交的硕±学位论文《基于EEMD的滚动轴承振动故障特征提取与诊断研究》,是本人在导师指导下,在华北电力大学攻读硕±学位期间独立进行研究工作所取得的成果。据本人所知,论文中除已注明部分外不包含他人已发表或撰写过的研究成果。对本文的研究工作做出重要贡献的个人和集体,均已在文中W明确方式注明。本声明的法律结果将完全由本人承担。■、作者签名:曰期;>年曰iJ喉縣^15^月作华北电力大学硕±学位论文使用授权书《基于EEMD的滚动轴承振动故障特征提取与诊断研究》系本人在华北电力大学攻读硕±学位期间在导师指导下完成的硕±学位论文。本论文的研究成果归华北电力大学所有,本论文的研究内容不得W其它单位的名义发表。本人完全了解华北电力大学关于、保存使用学位论文的规定,同意学校保留并向有关部口送交论文的复印件和电子版本,允许论文被查阅和借阅,学校可W为存在馆际合作关系的兄弟高校用户提供文献传递服务和交换服务。本人授权华北电力大学,可W采用影印、缩印或其他复制手段保存论文,:可W公布论文的全部或部分内容。""本学位论文属于觸在上相应方框内打V:)保密□,在年解密后适用本授权书不保密^‘作者签名:表喉^巧以曰期:年期/7曰I■导师签名:;日期年案月峰日 华北电力大学硕±学位论文摘要一滚动轴承是各种旋转机械中应用最广泛的种通用部件,其运行状态往往直接影响整台机器的性能。,对滚动轴承的状态监测和故障诊断具有重要的现实意义和经济价值一一些包含故障信在滚动轴承发生故障时,种常用的诊断方法就是,从振动信号中提取息的特征参数,从而揭示其故障类型、故障位置、故障程度。因此,如何有效的提取故障特征显得尤为重要。故障特征参数有时域参数和频域参数。常见的提取方法有小波分析和经验模态分解(EmpiricalModeDecomposition,EMD),而集合经验模态分解(EnsembleEmpiricalModeDecomposition,EEMD)作为经验模态分解的优化算法,近年来得到越来越广泛的应用。一、本文采用种方法,基于EEMD分解,提取W峭度均方根值为代表的时域特征参数和W重也频率、均方根频率为代表的频域特征参数,并运用支持向量机(SupportVectorMachine,SVM)进行故障识别。通过与EMD分解的结果进行对比,验证该方法。在诊断故障中有较高的正确率同时,在时域和频域中,多特征参数结合的诊断正确率一要高于单参数。;时域特征参数与频域特征参数相结合的诊断正确率可W达到更高关键词:滚动轴承;振动;故障特征;集合经验模态分解1 华北电力大学硕±学位论文Abstract民olli打gbearingis0打eofthemostwidelyusedpartsinvarityofrotatingmachinery,whoserunningstateoftendirectlafects化eerformanceofthewholemadii打esothefaulty,pdiagnosisofrollinbearinhasimortantracticalsinificanceandeconomicvalue.Whenggppg化ero化ngbearingfault,extractsomecharacteristicparametersincludingthefaultinformatio打:fromvibrationsignals,soas化revealUsfaulttypes,f加Itlocation,faultdegree,whichisacommonlyuseddianosticmethod.Thereforehow化efectivelextractthefaultg,yfeaturesisveryimportant.Thefaultfeatureparametersareclassifiedintotimedomainparametersandfrequencydomainparamet^s.TheWaveletAnalsisandEmiricalModeypDecompositionEMDaretwokindsofmethodforextractincommon.EnsembleEmirical()gpModeDecomposkionEEMDastheotimizatio打alorithmofEMDhasbeenaliedmore()pg,ppre猫dmowidelyinrecentears.yThispaperadopts江setofmethods,whichbasedo打theEEMD,化衍:tractkurtosis,rootmemsqu孤eV加eas化erepresentativeofthetimedomain,and化extractgravityfrequency,rootmeansquarefrequency,frequencystandarddeviationastherepresentativeofthefrequencydomain,andusingofSupportVectorMachine(SVM)forfaultidentification.Verifythattheme1:hodhashighcorrectrateinthe(diagnosisoffaults,throughcomparedwiththeresultsofEMDdecomposition.At化esametime,inthetimedomai打andfrequencydomaincombininmultifeaturearameter化anosiscorrectrateishiherhanthatofsinlet,gpgggparameter.Diagnosisofcharacteristicparametersoftimedomainandfrequencydomainfeatureparameterscombinationcanachievehi^ercorrectrate.Kewords:RollinBearinVibrationFaultFeaturesEEMDygg,,,II 华北电力大学硕±学位论文目录摘要IAbstractIIWiMmt1111.研究背景及意义2国外研1.究现状21.2.1解调方法研究现状21.2.2特征参数提取方法研究现状41.3本文研究目的和内容5第2章滚动轴承的故障特征72.1引言72.2滚动轴承基本结构72.3滚动轴承的故障形式72.4滚动轴承振动机理92.4.1轴承结构特点及加工装配误差引起的振动92.4.2轴承运行故障引起的振动92.5滚动轴承故障诊断常用统计特征102.5.1时域特征参数10■2.5.2频域特征参数112.6本章小结12第3章EEMD提取与SVM诊断算法13311.引言33.2经验模态分解133.2.1EMD原理133.2.2EMD分解算法142.2.3EMD算法存在的问题163.3集合经验模态分解163.3.1EEMD原理1640【]3.3.2EEMD分解算法173.4支持向量机17III 华北电力大学硕±学位论文319.5本章小结第4章基于EEMD时域特征参数的漆动轴承故障诊断204.1引言204.2数据说明204.2.1数据来源204.2.2数据分类214.3信号时域特征224.3.1原始信号时域特征224.3.2EEMD分解的信号时域特征244.4故障诊断方法设计274.4.1数据分组274.4.2诊断流程274.5基于峭度指标的滚动轴承故障诊断274.5.1基于EMD峭度的滚动轴承故障诊断27452EEMD29..基于峭度的滚动轴承故障诊断43.5.3分析与讨论0431.6基于均方根值指标的滚动轴承故障橡断4.6.1基于EMD均方根值的滚动轴承故障诊断314.6.2基于EEMD均方根值的滚动轴承故障诊断324.6.3分析与讨论i334.7基于时域多参数指标的滚动轴承诊断3541基于EMD时域多参数的滚动体故障珍断35.7.436.7.2基于EEMD时域多参数的滚动体故障诊断4.7.3分析与讨论374;38.8本章小结第5章基于EEMD频域特征参数的滚动轴承故障诊405.1引胃405.2信号频域特征405.2.1原始信号频域特征405.2.2EEMD分解的信号频域特征42IV 华北电力大学硕±学位论文545.3基于重也频率指标的滚动轴承故障诊断5.3.1基于EMD重也频率的滚动轴承故障诊断455.3.2基于EEMD重必频率的滚动轴承故障诊断46475.3.3分析与讨论548.4基于均方根频率指标的滚动轴承故障诊断548.4.1基于EMD均方根频率的滚动轴承故障诊断5.4.2基于EEMD均方根频率的滚动轴承故障诊断5051.4.3分析与讨论55.5基于频域多参数指标的滚动轴承故障诊断525.5.1基于EMD频域多参数的滚动轴承故障诊断525.5.2基于EEMD频域多参数的滚动轴承故障诊断535.5.3分析与讨论545.6基于时域频域多参数指标的滚动轴承故障诊断555.6.1基于EMD时域频域多参数的滚动轴承故障诊断565.6.2基于EEMD时域频域多参数的滚动轴承故障诊断引5.6.3分析与讨论5857本章.小结巧第6章结论与展望606.1结论606.2展望60参考^献62;攻读硕±学位期间发表的论文及其它成果65致谢66V 巧北电力人学硕」:学位论文第1章绪论滚动轴承的故障诊断在当今社会很多领域占据着十分重要的地位,不论是国外还是国内都发现了其内在的经济价值。如今国外在这项技术的研究已经发展的比较成熟,并正朝着智能化,集成化的方向努力发展;尽管我国目前还处于初级阶段,但是各高校纷纷结合国外的先进技术提出了自己独特的思想,并应用在实际生产之中,产生了巨大的经济效益。由于大量的旋转机械中都存在着滚动轴承,因此对设备进行实时监测与诊断,W就能够避免很多不必耍的经济损失甚至重大事故。1.1硏究背景及意义随着生产系统和设备的自动化、高速化和精密化的不断提升及其设备自身不断的大型化,在送个科技生产力持续发展的21世纪,各个部口不但对化产效率和成本提出了一更高的要求,而且对设计、使用、维修也制定了更严格的标准,个微小的故障,有时会对整个系统造成稳定性和安全性等方面的影响,甚至会导致整个系统的雍疾。在电力、机械一、航天等领域中,滚动轴承是使用最为广泛的种机械零件,同时也是很易受到损W伤的零件。因此,对滚动轴承这个至关重要的旋转零件进行故障诊断W及定期维修极大的引起了人们的重视。滚动轴承在机械设备中是一一种通用零件,也是机械设备中故障频出的部分之,据统计在旋转机械发生的故障中,大约30%的机械故障是由滚动轴承引起的;尤其是在感应电机故障中,甚至占有40%的比例在齿轮箱故障中,它也仅次于齿轮故障占比达到;了19%么多。查阅相关的资料,在我国,40%必,,机车用的滚动轴承每年有须下车检验而且大约33%的滚动轴承必须被更换,所W,机车轴承的故障诊断也具有很重要的意义和实用价值,。除此W外应用了故障诊断技术W后,机械设备出现故障的概率可W降低75 ̄,2550%%用在维修方面的费用减少可法%。实际的例子还有很多,1985,,。比如年大同电厂断轴毁机事故在短短不到40砂的时间内,发电机组全部损坏,000多方元。1988(直接损失达1年,由于轴系断裂汽P1)轮机发电机姐,位于陕西省的秦岭发电厂发生造成经济损失约3000方元的事故。其实,送样的例子屡见不鲜,,各国都开始关注故障诊断的问题正因为如此。英国根据调一查表统计,每3年单单在维修费用上就要花费亿英镑的巨资,如果多加研究故障诊断去避免故障就能大大缩减费用,因为研究这项技术的成本不高(5千万英镑)。法国的某电厂为了验证故障诊断技术的重要性一,经过两年的状态监测的统计结果显示,电厂共节约了95乂欧元,不过上述统计结果有些因素没有考虑(比如停工造成的损失)。国内方面,也有很多知名的高校做了相似的研究并发表了报告,西安交通大学就是很好的例子,应用自主研发的回转机械故障诊断与监测系统,在沪州天然气化学工业公司作1 华北电力大学硕±学位论文89%57,为试点,投入使用仅仅几年,就使生产加工中机组在大修中的开缸率从降到%W变向带来的经济效益十分显著(高达数千万元人民币)。由此可见,滚动轴承故障诊断技术可W掌握了解轴承的性能状态并且尽早的发现潜在的故障,除此之外,,对机械设备的替理水平和维修效能也会有显著的提升其经济效WW益也会十分的显著。1.2国内外研究现状,当滚动轴承局部存在损伤或缺陷时,在承载运转过程中将产生衰减冲击脉冲力从而激起轴承的高频固有振动。这种高频固有振动作为轴承振动的载波,其幅值将受到这些缺陷引起的脉冲激振力的调制,从而使轴承最终的振动波形表现为复杂的幅值调制波。一因此,解调分析成为机械故障诊断,特别是正确诊断故障部位的种常用的信号处理方法,在实践中解决了很多故障诊断方面的问题然而在机械故障诊断的发展过程中,人们发现最重要最关键而且也是最困难的问题之一就是故障特征信息的提取。对原始信号解调,是为了更好的进行特征提取。1.2.1解调方法研究现状(1)共振解调法。目前園内外普遍采用的行之有效的方法是共振解调法,又称高频共振法或包络分析法,即解调,。巧振解调巧术是对低频冲击所激起的高频共振波进行包络检波和低通滤波获得一个对应于低频冲击的、而又放大并展宽了的共振解调波。通过对此共振解调波进行频谱分析,判定故障的量值和类型。常用的共振解调方法有基于希尔伯特变换的包络解调方法、广义检波滤波解调方法和能量算子解调方法。希尔伯特变换是一种将时域实信号变为时域解析信号的方法。变换所得解析信号(复数)的实部是实信号本身,虚部是实信号的希尔伯特变换,而解析信号的幅值便是实信号的包络。由于对原始信号经过了包络处理,使原始信号变得光滑和平缓。在滚动轴承故障诊断中,由于可W采用希尔伯特变换提取机械故障信号包络的方法获得特征信息,,,同时希尔伯特变换又具有快速算法因此希尔伯特变换在机械故障振动信号的解W调分析中得到了广义的应用。广义检波滤波。高通绝对值分析、检波滤波和平方解调H种解调方法的解调原理都一是对信号取绝对值或绝对值的平方,然后进行低通滤波,可,原理基本致W统称为广W一一一义检波始波解调分析,,国内的些学者对其中的些算法进行了研究但未作为种整体方法对其解调原理和局限性进行详细系统的研究。文献详细地推导研究了广义检波滤波解调分析的解调原理和局限性。一能量算子解调方法。能量算子解调方法是由Teager提出的种简单的信号分析算法,用于分析和跟踪窄带信号的能量,对故障轴承的高频共。人们将其引入滚动轴承诊断中振信号进行包络解调,提取故障特征,试验和应用结果表明,能量算子解调法比传统的2 华北电力大学硕±学位论文ilfl希尔伯特解调法精度更高,速度更快。W上H种方法在应用中虽然取得一定效果一,但也存在些不足。希尔伯特变换法和广义检波滤波法除了存在抗干扰噪声信号的能力较弱1^^外,还存在几个局限性①将)不包括调制信息(故障信息的两时域相加信号;,也W其频率之差作为解调信号而解出②在广义检波滤波解调分析中,取绝对值,需、检波过程、平方过程可能产生混额效应1W,要选择合适的采样频率;③在几种细化解调分析新算法中对选抽后的信号进行频谱分析前不能进行数字低通滤波,有可能会出现调制频率高次谐波成分发生频率混叠而反折到低频部分的现象。能量算子解调方法虽然运算速度较快,但成功地运用该解调方法的前提同样是正确选取载波频段和调制频段,对原始信号进行带通滤波。所选择合适的共振频带进行带通滤波是术确进行解调分析的关键。(2)基于小波变换的解调方法。一1W小波变换是为了克服传统的傅立叶分析的局限性而提出的种时频分析方法,在时,,、频域都具有表征信号局部特征的能力具有多分辨率分析的特点可W探测正常信号中隐含的瞬态突变成分,很适合用来分析滚动轴承故障信号。目前,很多学者将其应usi利巧正交用于滚动轴承包络解调分析中。张中民等小波基将轴承故障振动信号变换到-尺度域时间,对高频段的小波系数进行包络细化谱分析,不仅能检测到轴承故障的存在,而且能有效地识别轴承的故障模态oNG。Niklaou等使用小波包变换作为分析系统振动信号的工具来诊断轴承的局部缺陷,结果表明,小波包变换能有效地提取振动信号的微弱瞬态特征。上文献采用二进离散小波变换和小波包变换,虽然有快速算法,但对于滚动轴承故障诊断而言,存在W下H点不足①时间和尺度二进网格的划分过于粗趟。二进小波的时窗长度按尺度降低方向(频率增加方向)呈二倍速减,其尺度划分由于基于二进划分而跳跃,离散间隙太大而过于粗檔;②二进小波变换要求小波基必""须是正交的,送样它才存在对偶,从而才能有重构公式。基小波的选择往往要与待""""检成分相匹配,而这种匹配性与正交性通常很难兼顾③二进小波变换不具有;""时(移)不变,在检测具有不确定时刻出现的信号的特征,检测方法是必须具有时一不变性的。人们投入了大量的精力來研究如何使.进小波满足时不变性,但这些改进虽,且要求基小波是正交的然可W实现时不变的变换,但是计算量显著提高。(3)基于经验模态分解的解调方法经验模态分解一(EmpiricalModeDecompositionEMD)方法是由Huan提出的种,g新的信号分解方法tW它根据被分析信号的内在结构特征,将复杂信号分解为若干个基,一本模态分量(Intrinsicmodefunction,IMF)之和,而每个IMF分量都可是个调幅和/或调频的单分量信号。EMD方法具有自适应特点,主要体现在该方法的分辨率是根据数据自身的特点,通过曲线拟合确定,突破了依靠先验经验的传统方式。轴承振动信号多为含有多个瞬时频率的多分量信号,借助于EMD恰能实现多分量振动信号的简单化处理。此外,轴承振动信号中的强信号成分往往将弱信号成分淹没,使得弱信号特征无一法体现,些方,而用EMD分解得到的某些IMF分量恰能对应着弱信号成分故避免了3 华北电为大学硕七学位论文法难分析弱信号成分的不足。目前EMD已在轴承故障诊断中得到了广泛的应用,如杨宇提出了EMD与神经网络顯合、EMD与奇异值分解融合等轴承故障诊断方法,在应用中有较高识别;王太P勇叫等瑞利炯、最小二乘支持向量机与EMD巧合,可对轴承故障进行智能诊断:李辉PI1MD和提出基于ETeager能量算子的诊断方法,可用于齿轮箱轴承的故障诊断;雷亚圍提出了EMD、自适应神经模糊推理系统和遗传算法擺合智能模型故障诊断方法,P31并在电力机车轮对轴承的故障诊断中得到应用;张琳将切片双谱理论融入到EMD中,该方法在实现了故障诊断的同时,还有较好的噪声抑制能力。然而EMD有典型的缺陷,原始信号间歇性W及分解过程中有误差,引起分解成的一IMF分量有多个频率成分或者单频率成分分布到几个分量中,送种和EMD分解初衷相违背的现象称为模态混叠。针对EMD的模态混叠,Huang又提出了集合经验模态分PW解(EnsembleEmpiricalModeDecomposition,EEMD)。该方法根据不相关量的统计一特性均值为零的数学原理,将白噪声信号加入到待分解信号中,方面白噪声信号可一弥补原始信号的缺失,另方面白噪声为不相关信号,其统计均值为零,这样既对模态混叠现象有很好的抑制作用生干扰。EEMD是EMD的改进,它,又不会对原始信号产同样具有完备性和近乎正交性,既能达到凸显局部特性的分解目的,还能有效抑制模态混叠。1.22特征参数提取.方法研究现状所谓特征提取,是指通过变换(或映射),把高维的原始特征空间的模态向量用低维的特征空间的新的模态向量來表达,从而找出最具代表性的、最有效的特征的方法。为了从根本上解决故障特征信息提取这个关键问题,人们主要借助于信号处理,特别是现代信号处理的理论和技术手段,从对信号的深度分析中获取更多的信息。前面总结了目前常用的解调和降噪方法,都是为了更好地进行特征提取。(1)基于时域和频域特征参数的提取方法滚动轴承发生故障时,其振动信号中的许多统计特征参数都会随故障的性质及大小发生变化,可作为故障诊断的依据。在轴承故障诊断中应用比较广泛的幅值参数有:有,及无量纲指标中的波形因子量纲指标的均方根值和峰值、峰值因子、脉冲因子峭度因子和裕度因子等。其中,对于元件表面损伤类故障,用峰值指标判断比较敏感,对磨损类故障,用均方根指标有效,而峰值因子既考虑了峰值又考虑了均方根值,所对两类故障都可W判断。但利用特征参数进行故障诊断只对故障的早期阶段较为敏感,当故障较为严重时,其抗干扰性较差,其值与正常状态值接近而产生误判问题。如均方根,从值可W反映轴承总体的劣化状况,但是,而使均方根值增大时,,当有其他非轴承振源、则易产生误判,但是,;峰值峰值因子、峭度等参数虽对冲击故障敏感当有其他冲击源或随机冲击干扰时,它们也会增大且,峰值因子、,当故障进入严重发展阶段时;况峭度等参数处于饱和状态,,,与正常轴承相同失去诊断能力。总之如果原始信号不加4 华北电力大学硕上-学位论文任何预处理,,W上参数的诊断能力较差,随机性较大。为此,为了提高诊断的可靠性需对原始信号进行必要的预处理,其采用多参数诊断法。(2)基于能量的特征参数提取方法在信号的分析过程中,时间尺度和随时间尺度分布的能量是信号的两个最主要参数。事实上,当滚动轴承发生故障时,振动信号各频带的能量会发生变化。因此,在各频带信号的能量中包含了丰富的故障信息,可根据滚动轴承振动信号的各频带能量的变化來进行滚动轴承故障分析。EMD可将信号在不同的尺度上展开,并可提取信号在不同频带上的特征,并从分解后的时间序列中提取轴承信号的能量特征信息作为神经网络的输slii入,对滚动轴承进行自动分类与诊断。杨宇等WEMD作为振动信号的预处理器,WBP(BackPropagation)神经网络作为分类器,提取EMD分解后的平稳IMF分量的能量指标作为网络的特征参数,对故障轴承和正常轴承进行分类。(3)基于滴度量的特征参数提取方法一滴是表明系统无序状态的某种测度,是系统复杂度的种定量描述工具,系统的状态发生变化时,其摘值也会改变。目前,显示其,滴的概念越来越广地应用到众多领域一PSj具有普遍意义,它已近成为研究数字信号处理的重要基础理论之。近年来,很多学一者将烦作为特征参数提取的方法引入到机械故障诊断领域,取得了定的应用成果。杨P世锡等叫尋K摘引入故障诊断领域,利用K赌的计算数值作为故障的特征值來识别机械的故障类型及严重程度。李兵等采用EMD方法对轴承故障振动信号进行分解,得到IMF矩阵,提取IMF矩降的奇异值和各IMF分量的能量值,计算香农滴和瑞利炯两种测度下信号的奇异值和能量分布赌,同时提化原信号的屯个统计特征参数共同构建轴承故障诊断的原始特征集。1.3本文研究目的和内容目前,滚动轴承故障诊断方法都先经过解调等预处理,然后提取特征参数。大多数特征参数需要某种模态识别工具进行辅助诊断,才能得到最终的诊断结论。通常需要大,使用中不是很方便量的数据样本进行训练,如果能找到某些参数直接表征滚动轴承的一一故障状态,将是件非常有意义的工作。本文采用种基于EEMD分解,提取时域、频域特征参数,然后通过支持向量机(SuortVectorMachSVM)ppine进行分类识别的,诊断方法,通过识别正确率证明该方法的可行性和有效性。论文结构安排如下:第2章首先从滚动轴承的基本结构、失效形式、振动机理等方面,从理论层面介绍了滚动轴承故障发生的机理和形式。然后介绍了基于时域和频域的常见特征参数及其提取方法。第3章主要介绍论文中用到的几种算法:首先介绍EMD的产生和原理,在故障诊断中的优势及缺点,;然后介绍了基于EMD发展来的EEMD的产生和改进原理及相较于EMD的优势。最后介绍SVM分类的基本原理及工作流程。5 华北电力大学硕±学位论文第4章主要介绍了W时域特征参数为指标,基于EEMD滚动轴承振动信号故障恃征参数提取方法理论,将原始信号分别进行EMD和EEMD分解,从得到的分量中提取时域特征参数,然后利用SVM进巧训练和分类,从而对故障种类进巧识别。同时,将EMD和EEMD的分类结果进行比较。第5章主要介绍了W频域特征参数为指标,基于EEMD滚动轴承振动信号故障特征参数提取方法理论,将原始信号分别进行EMD和EEMD分解,从得到的分量中提取频域特征参数,然后利用SVM进行训练和分类,从而对故障种类进行识别。同时,将EMD和EEMD的分类结果进行比较。第6章主要总结论文的研究成果,对滚动轴承故障诊断技术亟待解决的关键问题和将来的研究方向进行预测和展望。6 华北电力人学硕上学位论文第2章滚动轴承的故障特征2.1引言滚动轴承是旋转机械重要的支撑部件,要深入了解其故障特征,就要分析其结构特点一,进步总结故障形式及产生原因。由于滚动轴承故障诊断依靠其振动信号的采集、分化因此,要从内部、外部原因明确滚动轴承振动产生机理。然后归纳旌结常见的振动信号统计特征。2.2滚动轴承基本结构滚动轴承由内圈2-、外圈、滚动体和保持架四种元件组成,如图1所示。通常,轴承内圈与机械传动轴的轴颈过盈配合联接一,工作时随轴起转动外圈安装在轴承座、;而箱体或其他支撑物上工作时一,般固定或相对固定。但也有少数外圈回转、内圈不动或内、外圈分别按不同转速回转的使用情况。滚动体是滚动轴承的核也元件,它使相对运动表面间的滑动摩擦变为滚动摩擦。滚动体的形式有球形、圆柱形、圆锥形、鼓形等。在滚动轴承内、外圈上都设有凹槽滚道,它们起着降低接触应力和限制滚动轴承轴向移动的作用。保持架使滚动体等距离分布并减少滚动体间的摩擦和磨损。講\fi图2-1滚动轴承的基本结构2.3滚动轴承的故障形式滚动轴承在运转过程中可能会由于各种原因引起损坏,如装配不当、涧滑不良、水分和异物侵入、腐蚀和过载等都可能会导致轴承过早损伤。即使在安装、涧滑和使用维护都正常的惰况下一,经过段时间运转,轴承也会出现疲劳剥落和磨损而不能正常工作。7 华北电力大学硕±学位论文PS滚动轴承的失效形式主要有1:(1)疲劳剥落滚动轴承的内外滚道和滚动体表面既承受载荷又相对滚动,由于交变载荷的作用,首先在表面下一定深度处最大剪应力处形成裂纹,继而扩展到接触表面使表层发生剥()落坑,最后发展到大片剥落,这种现象就是疲劳剥落。疲劳剥落使机械在运转时产生冲击一、振动和噪声。在正常工作条件下,疲劳剥落往往是滚动轴承失效的主要原因,般所说的轴承寿命就是指轴承的疲劳寿命,轴承的寿命试验就是疲劳试验。试验规程规定,在滚道或滚动体上出现面积为20.5mm的疲劳剥落坑就认为轴承寿命终结。滚动轴承的疲劳寿命分散性很大一批轴承中,同,其最高寿命与最低寿命可W相差几十倍乃至几百一倍,这从另角度说明了滚动轴承故障监测的重要性。(2)磨损滚道与滚动体的相对运动和尘埃异物的侵入引起轴承表面磨损,润婿不良也会加剧磨损。磨损的结果使轴承游隙増大,表面粗髓度增加,降低了轴承运转精度,因而也降低了机器的运动精度,振动和噪声增大。对于精密机械轴承,往往是磨损量限制了轴承的寿命。(3)塑性变形当轴承受到过大的冲击载荷或静载荷时,或因热变形引起额外的载荷,或有硬度很高的异物侵入时都会在滚道表面上形成凹痕、划痕或压痕。送导致轴承在运转时产生剧烈的振动和噪声一。而且,压痕引起的冲击载荷会进步引起附近表面的剥落。(4)绣蚀与电蚀镑蚀是滚动轴承严重的问题之一。高精度轴承可能会由于表面镑蚀导致精度丧失而不能继续工作。水分或酸、碱性物质直接侵入会引起轴承诱蚀。当轴承停止工作后,轴承温度下降到露点、,空气中水分凝结成水滴附在轴承表面上会引起诱蚀。此外,当轴承内部有电流通过化电流有可能通过滚道和滚动体上的接触点化很薄的油膜引起电火花而产生电蚀,在表面上形成磋板状的凹凸不平。(5)裂纹与断裂轴承元件的裂纹和断裂是最危险的一种损坏形式,这主要是由于轴承超负荷运行、金属材料有缺陷和热处理不良所引起的。转速过高、润滑不良、轴承在轴上压配过盈量太大W及过大的热应力都会引起裂纹和断裂。.(6)胶合在润滑不良、高速重载情况下工作时,由于摩擦发热,轴承零件可在极短时间内达到很高的温度,导致表面烧伤及胶合。(7)保持架损坏装配和使用不当可引起保持架发生变形,增加了它与滚动体之间的摩擦,甚至使某些攘动体卡死不能滚动,也有可能造成保持架与内外圈发生摩擦等。送些都将增加轴承的振动、噪声与发热加剧。,导致轴承损坏8 华北电力大学硕I:学位论文2.4滚动轴承振动机理在机械运转时,由于滚动轴承本身的结构特点、加工装配误差和运行过程中出现的故障等内部因素,W及传动轴上其他零部件的运动和力的作用等外部因素,当传动轴W一一定的速度并在定的载荷下运转时,,对轴承和轴承座或外壳组成的振动系统产生激励使该系统振动2-2,其振动产生的机理如图所示,在实际的诊断中,如果不考虑轴承加工和装配误差一,则主要为运行故障这内部因素所引起的振动。结构特点IIH1W1内i轴加工装配 ̄?轴承系统?振动j護I^!I素八ii运行故障H—L—i:外部因素iIII图2-2滚动轴承振动产生原理实际诊断中,通过布置在轴承座或外壳上适当位置的传感器拾取的振动信号是上述各种内部和外部激励施加于滚动轴承系统的综合振动。如何从综合振动中把轴承故障引起的振动信号提取出来,从而有效地诊断出轴承故障,是滚动轴承故障诊断技术的关键所在。为了能够有效地把故障特征从综合振动中提取出来,需要研究滚动轴承内部两种因素引起的振动信号的各自特征2.4.1轴承结构特点及加工装配误差引起的振动当轴承旋转时,,滚动体便在轴承内、外圏滚道上滚动即使对加工装配毫无误差的轴承来说,由于滚动体在不同位置时所受的力大小不同,同时承载的滚动体数目也不同,一这些轴承本身的结构特点造成承载刚度的变化,引起轴承振动,,当轴的转速定时这一振动具有确定性质;轴承元件加工时留下的表面波纹度、粗精度及形位误差和装配误差等因素产生交变激振力使轴承系统振动,虽然送些加工装配因素造成的激振力大都具有周期性特点,但由于实际构成因素十分复杂,各因素之间也不存在特定的关系,所W,总体來说这些激振为随机性较强,含有多种频率成分,那么轴承系统在这些激振力的作用下所产生的振动当然也具有多种频率成分并具有较强的随机性。2.4.2轴承运行故障引起的振动滚动轴承在运行过程中出现的故障按其振动信号的特征不同可分为两大类一:类是一磨损类故障,,;另类是表面损失类故障包括点蚀、剥落、擦伤等下面分别具体介绍。一(1)磨损类故障:般来说,在正常使用情况下,滚动轴承工作表面磨损故障经一历的时间较长,是个渐变性故障。轴承表面磨损后产生的振动同正常轴承的振动具有相同的性质,即两者的波形都是无规则的,随机性强。但磨损后振动水平幅值)明显高(9 华北电力大学硕±学位论文于正常轴承。运就是磨损类故障引起的振动信号的基本特点。由于磨损不会马上引起轴承破坏,其危害程度远小于表面损伤类故障,故通常人们最为重视的还是表面损伤类故障。(2)表面损伤类故障对于表面损伤类故障,当损伤点滚过轴承元件表面时,产生突变的冲击脉冲力,该一脉冲力是宽带信号,所W必然覆盖轴承系统的富频固有振动频率而引起谐振,从而产生周期性冲击振动。送就是损伤类故障引起的振动信号的基本特点。该类信号的处理是本文重点要研究的内容。损伤故障产生的周期性冲击振动成分从性质上可分成两类:其一,是由于轴承元件的工作表面损伤点在运行中反复撞击与之相接触的其它元件""表面而产生低频振动成分,称之为通过振动。其发生周期7;是有规律的,可獻从转=速和轴承的几何尺寸求得。并且,损伤发生在内、外圈或滚动体上时,该频率人1/打)("一"不同。这轴承通过振动发生的频率称为故障特征频率。利用频谱分析诊断轴承故障时,其基本原理就是查看轴承振动信号中有无这些故障特征频率成分。滚动轴承故障一一特征频率般在IkHzW下,是滚动轴承故障的重要特征信息之。""其二,是由于损伤冲击作用而诱发的轴承系统的高频固有振动成分。这里所说的高是相对故障特征频率而言。轴承系统的高频固有振动很复杂,如轴承内、外圈的径向弯曲固有振动、滚动体的固有振动、甚至测振传感器的固有振动等都可由于损伤冲击而产生并反映在轴承的振动信号中。2.5滚动轴承故障诊断常用统计特征,在渡动轴承状态监测和故障诊断的振动特征分析中,特征参数法占有重要位置。特征参数法的优点在于仅用少数指标来解释轴承的状态,结果分析简单和方便。此外,采用特征参数法评定轴承状态时不需要轴桌W前的历史记录一,这是此技术另项具有吸引力的特征。在滚动轴承诊断中常用的特征参数包括有效值、峰值等各种时域恃征参数和重也频率等各种频域参数。下面结合实测数据对各种特征参数进行分析2.5.1时域特征参数,时域特征参数有均方根值、峰值、峰值因子、峭度等,送里仅选择有代表性的两个。特征参数进行介绍,并运用于后文的特征参数提取与诊断中(1)均方根值(RootMeanSuareRMSq,)均方根值又称为有效值。滚动轴承振动的瞬时值随着时间而不断地变化,作为表现这种振动变化大小的方法,广泛地使用有效值。有效值是振动振幅的均方根值,定义如下:。=2-1屑肉品()10 华北电力大学硕±学位论文其中,X的为时域信号。。=品腐口巧’=1...其中,2]\为采集的信号,^为采样点数。切〇}〇,,,〇公式2-1为连续时域信号求均方根值的公式-),22为离散信号求均方根值的公式。(()由于有效值是对时间平均的,,所W对表面磨损那样的无规则振动波形的异常其测定值的变动小,并可W给出恰当的评价。但是,对表面剥落或压痕那样的具有瞬变冲击振动的异常是不适用的。其原因是由于冲击波峰的振幅大,但持续时间短,,如对时间平均那么有无峰值的差异几乎表现不出来。当滚动轴承无故障时一旦轴承出现了损伤,有效值比较稳定;,则产生冲击信号,振动峰值明显増大,但此时有效值尚无明显的増大;当故障不断扩展时,峰值逐步达到极限值后,有效值则开始增大,直至恢复到无故障时的大小。峰值因子对表面磨损故障几乎没有检测能力。(2)峭度(Kurtosis)对于零均值的振动时域信号x(f);,峭度定义为庐(〇如)叫(其中P的为蛛)的概率密度分布函数。离散化的公式为;4:衣(如()*K二2-4F7f()?!>!(()=11__对于其振幅满足正态分布规律的无故障轴承,其峭度指标值约为3,随着故障的出现和发展,峭度指标值具有与峰值因子类似的变化趋势。其优点也在于与轴承的转速、尺寸和负荷无关。但主要适用于点蚀类故障的诊断。2.5.2频域特征参数\时域参数指标判断滚动轴承故障基本上不受轴承型号、转速、载荷、信号绝对水平和测点等的影响一,但只能判断轴承是否有故障。而要进步确定故障的类型,发生故障的元件,就必须根据频谱图中的频率成分处的幅值的大小进行诊断。选取两个主要的频域参数分别为;(1)重也频率:11 华北电力大学硕±学位论文=C2.S/培户()s/d/()J7(2)均方根频率:斯r八(耸三f|各式中的S的均为信号的功率谱。重也频率、均方根频率都是描述功率谱主频带位置变化的。利用傅立叶变换的性质,’=可得信号扣闲〇12...iVJV为采样点数的频域参数重必频率、均方根频率的时域快},,,;)速算法为;、.2-7()请2、4尤巧=(〇,^,户!(>(0C,口广姆訓>x-x-l(〇()^2.9()么其中,/!为采样间隔,iV为采样点数。采用频域参数可W对滚动轴承的故障进行粗略判断。当轴承无故障时,频率成分主,重也频率较小。当轴承出现局部损伤类故障时要集中在低频,由于冲击引起的共振,主频区右移,重也频率增加。2.6本章小结本章系统的介绍了滚动轴承结构,常见故障形式及其原因。并从振动机理入手,阐述了轴承振动包含的信息,为滚动轴承振动信号分析提供T理论依据。从时域、频域角度,详细描述了滚动轴承信号统计特征,为特征参数计算提供了计算依据。12 华北电力大学硕±学位论文第3章EEMD提取与SVM诊断算法3.1引言在本文提出的滚动轴承敌障诊断方法中,,主要运用了H种算法对原信号进行预处理的算法有经验模态分解EMD和改进后的总体经验模态分解EEMD,提取了特征参数—M后,还用到对故障类型进行分类识别的算法支持向量机SV。本章将着重介绍这兰种算法的原理。3.2经验模态分解EMD经验模态分解法是美国国家宇航局美籍华人黄嚮(化E.Huang)等人于1998年创造性地提出的一种新型自适应信号时频处理方法。EMD一被认为是2000年来W傅立叶变换为基础的线性和稳态频谱分析的个重大突破,该方法是依据数据自身的时间尺度特征来进行信号分解,无须预先设定任何基函数。这一点与建立在先验性的谐波基函数和小波基函数上的傅里叶分解与小波分解方法具有本质性的差别。正是由于这样的特点,EMD方法在理论上可W应用于任何类型的信号的分解,因而在处理非平稳及非线性数据上,具有非常明显的优势,适合于分析非线性一经提出就在不同的工程、非平稳信号序列,具有很高的信噪比。所W,EMD方法PU领域得到了迅速有效的应用。3.2.1EMD原理EMD的核也由两部分组成,即瞬时频率和本征模函数,要透彻地理解EMD的原理,首先要理解瞬时频率和本征模函数的概念。(1)瞬时频率假设一信号;C似可写成指数信号的iV项相加,即.'"(=.)。6-刪12]:*口)1其中化为常数。则瞬时频率频率表示)==-2kl…N口)从,,)鸣心直观上,瞬时频率为相位的微分。对于自然界中的实数讯号,如何定义讯号的相位。G油or提出解析讯号法,将实数讯号表示为对应的复数讯号,即可定义复数讯号的大小与相位,将实数讯号的瞬时频率求出。(2)本征模函数若使一个函数或者信号的瞬时频率有意义,那么这个信号或函数必须是对称的并且13 华北电为大学硕±学位论文局部均值等于零,同时其过零点数目和极值点数目相等。基于此,Huang等人提出了本一一征模函数一(IMF)这概念,他们认为信号是由系列本征模函数和个余量组成的,即;=c+-r3玄L,(0口)非),(0C?为其中,IMF分量。3一PHu必须满足W下两个条件3ang指出个本征模函数(MF):一(1)在整个信号长度上,极值点数目和过零点数目必须相等或者至多相差个;(2)在任意数据点处信号上下包络的平均必须等于零。一u然而,在实际中遇到的信号般是不满足上述的理想条件的,于是,Hang又做了如下假设:一一^一(1)任何信号都是由系列的本征模函数和个余量組合而成;(2)在局部范围内,组成信号的各个本征模画数可W是线性的也可W是非线性的,并且各本征模国数的零点和极值点数目相等,同时上下包络线也是关于时间轴对栋的;一一(3)不论在什么时候,个信号都可W包含有系列的本征模函数,如果这些本征模函数发生混叠,就会组成复合信号。3.2.2EMD分解算法fWEMD是通过筛选过程來获取MF:,其具体分解过程为①确定信号;C的的所有局部极值点,然后用H次样条插值分别将所有的局部极大值点和局部极小值点连接起来3-,形成上包络线和下包络线,如图1,包络线包络了所有-的信号数据。3-2②如图,计算上、下包络线的平均值,记为m,:i求出=-hxtm3-4^),{)(一一一③如果如是个IMF,那么它就是信号;c〇)的第个IMF。如果不是个IMF,则将知作为原始数据,重复①、③步骤,得到上、下包络线的平均值,记为mi,,计算==IM/i,>,并判断是否满足F的条件,如果不满足,则重复循环,计算>1it一一=F直到Aw是个IMF。记C,则C为信号JC的的第。l个IMi将Cl从x(句中分离出来,得到;=-rxtc3-5,^{)() ̄二c。重复循环n次;C将n作为原始数据,重复步骤①⑨得到第个MF2,到信号的的n个IMF,于是有:-…r3旬。(一当成为个单调函数不能再从中提取满足MF条件的分量时。这样,循环结束EMD分解将信号表示为:14 华北电力大学硕±学位论文=-。7非),^口)茲|-其中,称为残余函数,代表信号的平均趋势。分解结果如图33。2IIIIIIIIIII1020304050607080901001101203-簡1求出信号的上、下包络线2tIIIIIIIIII1020304050607080901001101203-2计算出信号的上團、下包络线的平均值1尺A/\/WVA八羣-■——i[■IIIIIIIIIII—^?SIII?IIIIIII01020304050目07080如100110120-3图3EMD分解结果15 华北电力大学硕±学位论文33EMD算.2.法存在的问题(1)EMD分解存在边界效应,尽管有很多减小边界效应的方法,但仍不能从根本上消除。(2)现有的EMD分解算法还不算很完善,信号在分解过程中会产生这样或那样的一方面会产生较多的分量误差,,从而带来附加的低频信息。(3)由于信号的复杂性,可能会产生各分量之间存在模态混叠的现象,不利于信一号的进步分析。3.3集合经验模态分解一模态混叠是指IMF分量中包含了尺度差异较大的信号,或者个相似尺度的信号出现在不同的IMF分量中。当出现模态混叠时,会使单独的IMF分量缺乏物理意义,而且会使EMD算法不稳定。一为了减轻模态混叠现象an一一IHug等提出种新的白噪声辅助数据分析方法PyEEMD集合经验模态分解。3.3.1EEMD原理模态混叠之所W会出现,往往与信号的间歇性有关,EMD是按照由島到低的频率次序进行,,而信号的间歇性会导致当分解高频信号不完善时,为了满足算法需要就会一频率中调取一部分成分来弥补该部分的缺失从次高频或者此次分解频率的下,W至于PWEMD将高频成分变成了混合成分。除此W外,如果两个成分的频率和幅值极为相近,基于算法自身的原因可能无法将两种频率完全分辨出来而出现奈舌L,也将影响到分量的一频率成分还有点就是噪声信号往往具有窩频特性,,其频率时常比真实信号还有高在分解过程中,它会取代真实信号高频成分的位置,使得分解的分量混杂着噪声成分,一PW也是模态混叠出现的个原因。在这些原因中,信号的间断性的影响占到了非常大的比重。一主要影响因素一为了消除模态混叠,特别是针对信号的间歇性这,种主要思想是加入特定的成分来改善信号间歇处的不足,使信号变连续,同时还要考虑到加入的特定成分不能对原始信号产生影响,W便保持信号原有特性。基于这些考虑,总体平均经验一模态分解(EEMD)应运而生,种基于噪声辅助的信号处理方法依据统计学理论一中不相关序列的统计均值为零这原理,向待分解信号中加入噪声,由于噪声为随机序列,频率均匀分布,故它在使原始信号中的间断变连续的同时,又因为统计均值为零,可相互抵消其对信号的影响。EEMD是EMD算法的改进,它除了保留了EMD的优点之外,又可有效抑制模态混叠,提高信号分析的精确度。16 华北电力大学硕±学位论文W33.3.2EEMD分解算法一般来说,所有信号的测量数据都是由信号和噪声组成的,即:=-;cfsf+w38()()的()送里,X的是记录的数据,5阳和U0是真实的信号和噪声。为了提高测量的正确性,类似于物理实验中对单个物理量的测量通过多次测量取平均值,在数据球仲引入白噪声,附加的白噪声被认为是在测量过程中可能遇到的随机""噪声:,贝帕第次人为的观测值为'=xf+-W39式的()f()一尽管附加的噪声可能导致较小的信噪比一,但附加的白噪声将提供个均化的相关尺度分布使EMD更加容易。因此,低信噪比并不影咱分解方法,实际上可W帮助避免模态混叠。基于这个观点,加入白噪声这个附加的步骤可W帮助在数据中提取真实的信号,这种方法就是EEMD。其实现过程如下;①为待分解信号蛛)加入白噪声序列,得到信号是;(0②将信号;W按照EMD步驟分解成不同的IMF分量;③将①、②两步重复多次,但是每次都加入不同的白噪声;④求取相对应的IMF分量的平均值作为对应的最终的IMF分量目可表述为;;二f-刪茲5+(310),3.4支持向量机V一支持向量机SVM是apnik等人根据统计学习理论提出的种学习方法,它在解决模态识别的小样本问题中表现出了独特的优势和良好的应用前景,己用于模态识别和特征提取等领域。SVM可W采用少量的时域故障数据样本训练故障分类器,不必进行信号WU预处适来提取特征量,便可实现多故障的识别和诊断。SVM是从线性可分情况下的最优分类面发展而来的-4的两,基本思想可用图3维情况说明。图中,五角星和圆圈分别代表两类样本,H为分类渐H、分别为过各类中i离分类线最近的样本且平行于分类线的直线,它们之间的距离叫做分类间隔。所谓最优分类线就是要求分类线不但能将两类正确分开(训练错误率为0),而且使分类间隔最大。一分类线方程为=我们可W对它进行归化,使得对线性巧分的样本集.3=一乂.£!1..丹丈£灰+11,满足,乂,,,,,少,(;){}.-.-仍f'+61含=...:0!1只.,311义,,[(,)]()22仿3-此时分类间隔等于2/,使间隔最大等价于使H最小。满足条件(1。且使仿最小||III^wai的分类面就叫做最优分类面,H、H2上的训练样本点就称作支持向量。i17 华北电力大学硕±学位论文支持向量/分类间隔紙A相支持向量II玲*最忧分冀线3*4图SVM分类示意图Vapnik给出了求该最优分类面问题的解,即分类函数为:巧=nx+6-1/s32aWg2)w(v)()-M_其中乘子a,为目标函数;nn1=ra"VX3-()2>,2>,乂A/)(巧(^==1111,7在约束>=0i-l2...n314,,,,()■=0-巧狂讯(3=il?对应的样本点下的极大值点,非零a;c<为支持向量,阔值:,b=-.mX3-1yii(巧n。=-化317)弘(=/1上面讨论的是线性问题,对非线性问题,可W通过非线性变换转化为某个高维空间的线性问题一般来讲比较复杂我们可注,在变换空间求最优分类面这种变换,但是意到--,在上面的问题中,不论是寻优函数(313)还是分类函数312都只涉及训综样本之()间的内积运算>1,这样,在高维空间实际上只需进行内积运算,而这种内积运算是可1^用原空间中的函数实现的,我们甚至没有必要知道变换的形式。根据泛函的有关理论,只一一要种核函数满足正定矩阵,它就对应某变换空间中的内积。因化在最优分类面中采用适当的内积函数就可W实现某一非线性变换后的线性分类-,而计算复杂度却没有增加,此时目标函数313变为()18 华北电力大学硕±学位论文nn2--=x?WaaUaKX-{),y,,31,巧JY^jjj)(^(==il/,yi相应的分类函数(3-1巧也变为巧=n+&-xsg及和义19乂()公片|()(3)=/1_-tW常用的核函数有线性核函数、多项式核函数、径向基核函数、Sigmoid核西数。3.5本章小结本章详细介绍了在后文会用到的两种算法的原理及计算公式。在阐明EMD算法的原理后,同时总结了EMD算法存在的问题。从其中横态混叠问题解决办法入手,引入了EEMD算法,详细阐述了EEMD相对于EMD改进的原理。最后介绍了SVM算法的原理与基本步骤,显示了机器学习的优势。19 华北电力大学硕±学位论文第4章基于EEMD时域特征参数的滚动轴承故障诊断4.1引言在滚动轴承状态监测和故障诊断的振动特征分析中,特征参数法占有重要位置。特征参数法的优点在于仅用少数指标来解释轴承的状态。此外,采,结果分析简单和方便用特征参数法评定轴承状态时不需要轴承W前的历史记录一,这是此技术另项具有吸引力的特征。然而,在工程中测取的滚动轴承故障振动信号常常受到随机噪声的干扰,使得信号的信噪比很低,特别是发生早期故障时,其微弱的故障信息完全淹没在噪声中。一而EEMD是种自适应的信号分解方法,利用该方法不仅能分离部分噪声,而且不会对信号造成扭曲,能较好地保留原始信号特征。若将两者结合起来就能取得较好的效果。而在传统的故障诊断中,,,对故障的识别主要依靠时、频谱分析从中找到故障特征并依据经验或者故障信息库进行对比识别,效率和正确率均得不到保证。本章仅用时域特征参数进行滚动轴承故障识别。4.2数据说明4.2.1数据来源(1)实验装置本文实验数据取自美国凯斯西储大学电气工程实验室的滚动轴承故障模拟实验台获取的滚动轴承故障数据。该数据可在美国凯斯西储大学轴承数据中必网站下载。该一一一-实验装置如图41所示,该装置由台2马力的电动机,个扭矩传感器/编码器,个6205-2RSEM功率测试计,还有电子控制器姐成。驱动端轴承型号为SKFJ,风扇端轴承SKF6203-2民SEM4-1型号为J,表列出了驱动端和风扇端轴承尺寸信息。电动化八功率测试计扭矩编巧器:?.三!::!!!!!!^贩.勺图4-1实验测试平台20 华北电为大学硕l:学位论文表4-1驱动端轴承尺寸信息单位nchi ̄11内圈直径外圈直径ii滚动体直径轴承节径-SKF62052RS0.98432..5.04720.590603126137SKF6203-2民S0.66931.57480.47241.122(2)信号采集在风扇端和驱动端轴承座上2一1点钟方向各放置个加速度传感器来采集振动信息。612振动加速度信号由1通道数据记乗仪采集,风扇端轴承故障采样频率为kHz,驱动端轴承故障采样频率有两种:12kHz和48kHz。实验根据电机负载不同,共有四种实验转速:空载时候转速为1797r/min,负载为1马力时候转速为1772r/min,负载为2马力时候转速为1750r/min,负载为3马力时候转速为1730r/min。(3)点蚀故障设置实验通过电火花加工的方式分别为驱动端和风扇端轴承内圈、外圈和滚动体制造点蚀损伤,损伤直径分别为0.007inch、0.021inch、0.028inch和0.040in沈。山于外圈损伤的位置是相对固定的,因此损伤点相对于轴承负荷区的不同位置对电动机/轴承系统的振动响应有直接的影响。为了量化这种影响,驱动端和M扇端轴承外圈的损伤点分别放置在3点钟6、12、点钟点钟王个不同位貴。4.2.2数据分类--本文采用驱动端采样频率为12kHz的故障信号,为了减少变量因素,只用种负载下的数据,确定为负载为1马力转速为772rmn的数据。其中外圈故障采用外圈损伤1/i点在6点钟方向的信号。采集的信号依据不同的参数都保存为单独的.mat文件。il根据点蚀直径的不同,可认为是滚动轴承的故障程度不同,因此7m代表故障早期或者轻微故障,28mil代表故障晚期或者严重故障,这就有了不同故障程度的数据。本文设计的诊断方法,诊断结果不仅指示故障类型,同时可W显示故障程度。假设故障特征参数不仅与故障类型有关,还与故障程度有关,即同种故障不同程度,其特征参数规律不同。W此为依据,,可[^1将信号数据分为12类对应不同故障类型的不同程-度,42。这样划分的优势就是可W判断故障程度,见表,有利于制定不同的维修方案有实用价值。4-2表滚动轴承测试信巧分类单位mli轴承状态内岡故降外岡故降滚动体故巧‘甲 ̄ ̄ ̄ ̄ ̄ ̄11,J点蚀直径71421瓦1142^1141\百M12345678910U一考虑到存在误差的情况,送种分类具有优势SVM分类时,某个2,如果在組信号被识别为3组,这个结果在12组分类中是错误的,但是从故障类型角度来看,送个"",212组分类,错误是可W接受的因为、3组同属于内圈故障。因此,采用在只考21 华北电力大学硕±学位论文虑故障类型而不涉及故障程度时,可W化为4组分类,且识别正确率较高。4.3信号时域特征4.3.1原始信号时域特征(1)不同故障类型时域波形本章基于滚动轴承振动信号的时域特征参数为依据进行故障诊断,在时域信号进行正常与故障的时域图是不同的-上到下分别是正常、内圈EEMD分解前。如图42,从,故障、外圈故障、滚动体故障的时域波形图。正常信号时域波形II1iI1III一CM2-JOElaU微-畏-2-111!:Ii]1.100..060.070.080.090.010.020.030.040050内圈故障信号时域波形IIIII!I1—M2-C0v^|l沒-畏-2--1I1.II1III.100...040.050.060.070.080.090.010020030外圈故障信号时域波形4i■iIII[1I」I!1■J1J_?I—I.090.100.0.040.050.060.070.080.010.02030滚巧体故障信号时域波形 ̄*!IIII1I[IjCM2-巧-导-2-IIII0.1.060.070.080.090001.020030.040.05.00.时间s图4-2滚动轴承正常与故障的时域波形对比22 华北电力乂学硕:|:学位论文,通过对比可看出,正常信号相比于故障信号,加速度幅值较小没有明显冲击,且没有表现出明显的高频周期特征,这是因为正常信号中包含的是基本轴转频率及其谐一。波,总体保持在低频段通过对比时域图,内圈故障与外圈故障有定相似性,都有较强的冲击,且都表现出高频的周期性,但内圈故障显得有些杂乱,说明其中包含的频率成分比外圈故障要多。滚动体故障的时域波形图与内圈、外圈故障的时域波形图有很大区别,,较为杂乱看不出明显的周期性,故障信号、噪声信号等将正常的轴转信号完全掩盖,这是因为滚动体在轴承旋转时,不仅随着内圈转动,还有自身的转动,因此时域波形较为复杂。点蚀直径为7mli■1IIIIII1--CM2齊导-2-IIIII1III00.010020.030.040.050.060.070.080.090.1.点蚀直径为14mil]IIIIIII]CM2-‘脅li畏-2-IJII11III00.010.020.030.040.050.060.070.080.0901.点蚀直径为21milI]I[IIIIiS--I2C导-2-'-IIIIII!II00010020.030......040050.060070080.090.1点蚀直径为28milIIIIIIIIiI'j,巧-2-!-Ij:][jIjI00.010.020.030.040.050.060.070.080.090.1时间s圏4-3滚动轴承内閣故障不同点蚀直径的时域波形对比23 华北电力大学硕±学位论文同时,,对比可W看出正常信号振动加速度幅值很小,而内圈故障、外圈故障的振动幅度放大了10倍左右,而滚动体故障振动幅度则小很多只比正常信号略大。,(2)不同故障程度时域波形-、28,43,7m14i21mmil,W内圈故障为例图数据选择点蚀直径依次为il、nil、il代表故障从早期到严重等四个阶段。通过对比可W看出,故障程度不同,时域波形也表一定差异。故障早期振动的周期性较为明显,且振动幅值也比较大现出;之后随着故障发展,频率成分发生了变化,周期性冲击减少,但振动幅值依然比较大;之后在振动幅值不变的情况下,,,巧击次数增加高频成分增加到了故障晚期周期性表现不明显,;显得杂乱,但振动幅值有所增大。图中的对比可W看出点蚀直径为Mmil时的时域波形一i反而比7ml的波形振幅要小,且频率成分单,应该是滚动体经过直径为14mil的点蚀产生的振动与固有振动形成了消减。,导致频率成分减少4.3.2EEMD分解的信号时域特征W内圈故障为例,分别进行EMD和EEMD分解。(1)EMD分解后的信号时域特征-,得到8个,4。按照EMD进行EMD分解后IMF分量各IMF分量时域波形如图4,的原理分解得到的,,IMFIMF1分量最接近原始信号随着分解不断进行分量成分越一一来越简单,直到最后个残余分量呈现单调减,而IMF也按照不同尺度,从高频到低頻尺度分A。理想的目MD分脾是,8个IMF,还可分豈和残余分望袖互涵加til还原原始信号,。从分解到重构的过程看;其实就是个减法到加法的过程,减法求异剥离出频率(周期)大致相同的IMF,而加法求同,回到原始波形。残余分量其实是个趋势线,即频率极低(周期很长)的波,可W看成是个基底,其它IMF都构筑在它之上。(2)EEMD分解后的信号时域特征-,得到9个IMF分量5。进行EEMD分解后,各IMF分量时域波形如图4对比可EEMDEMDEMD一IMFW看出,与分解过程基本类似,但E相比多个分量。从图中发现EEMD的第6和第7个分量波形图比较接近,而EMD分解的各个分量之间差异较大;""""且EMD的分量的波形比较光滑但不够均匀,,而EEMD的波形却更为均匀均匀一分量中的频率成分较为集中就代表同。而且实际中的振动信号往往包含山化频到高频的复杂频率,甚至存在调制现象。EEMD针对时间尺度的划分更细致,也就从分解的角度,尽可能避免了模态混叠现象,使频率成分的分离更为清晰。一EEMO与EMD的差异在一后章节中会从频域角度进步分析,而且在本章的故障诊断方法中,通过适量的信号训练,经EEMD分解的分量比经EMD分解的分量所包含的故障特征更为正确,在故障诊断识别中,指向性更为明确。24 华北电力火学硕±学位论文21112〇0.0200400600801....0'.5''己_IIIIQ500.020040060080.1...05.Ii!I^'-^vwwAwwv0l^J1lVwvvwvwv\r/Wwvvl/v/)J^A\|V\/lVA/\\y/YJ/\//養y|^_|jJ一'^i^'-0.500.020040.060.080.1.0111..,0^"^mwa^WVWV/wVWvW\麵WVVV\/-''i'01.00.020.040.060.080.101.05111[0-A?/\A^WV\/\/VIkAA/WV\/-'i'005'.00.020.040.060.080.10.021.11—0QgII1I'日-30.020.040.060.080.1X10丰IIII日-0.0203.040.060.080.1X10"""''■■'"弓IIIII。画dIIII|0-0.020.040.060.080.13X10广'5III'___0-^ ̄當:_gIIII00.020.040060.080.1.图44内圈故障EMD分解図25 华北电力大学硕-上学位论文—21-.」I00020.040..060.080.1*05- ̄r1.I2-1jI105_00.020.040.060.080.1-■0III.5^''^''^'"^^wVlrWW'^0^/WwWWtaAWi/Av%^\\//VW\/(yV\/AWLrVw\A^yyy|yy\/y|^|L-0——^-—J—I^.500.020.040.060.080.1-0I^一^一^L1.00.020.040.060.080.1—下一0-,rr.05^0L/A/一??ILij-0■.05000200400600801.....0.05'111。:,'、'、'■'、'-’—?? ̄- ̄ ̄—■—'、、———、■—.'^—-^_一^S0、-_-iI-j0■0500.020.040060080.1..0.02....^I’^0一—I''-002''.00.020.040.060.080.1、^"-.?一0^?I\一■'''^-0.0100020...040.060.08010i.0111i1—— ̄'-■'1'0.0100.020.040.060.080.10051111.'= ̄■—' ̄ ̄w0i^——''-0'.05400.020.040.060.080.1图4-5内罔故障EEMD分解图26 华北电力大学硕±学位论文4.4故障诊断方法设计4.4.1数据分组进行口组分类故障识别,其中正常取80姐数据,内圈故障依4种不同故障程度分别取40组数据,共160组数据:外圈故障依3种不同故障程度分别取40组数据,共120组数据;滚动体故障依4种不同故障程度分别取40组数据,共160组数据,总计520一半输入組数据。不同状态的数据均等分为两部分,SVM支持向量机作为样本进行学习一,另半作为测试样本用来识别。4.4.2诊断流程具体诊断方法流程图如图4-MDEMD12所示,它采用E(E)方法对滚动轴承振动信号进行预处理,然后求出所有MF分量的时域特征参数参数,组成故障特征向量SVM作为分类器来识别滚动轴承的工作状态和故障类型。一(.1)按定的采样频率分别对滚动轴承的几种状态(正常、内圈故障、外圈故障、滚动体故障)进行采样;(2)对每个样本信号进行EMD(EEMD)分解,得到各信号的各IMF分量;(3)求各信号和各信号的各MF分量的时域特征参数,组成特征向量;一(4)将滚动轴承的几种状态(正常、内圈故障、外圈故障、滚动体故障)的部分特征向量输入SVM,对SVM进行训练;(5)将剩下的特征向量作为测试样本输入SVM进行分类计算,并将分类函数的输出结果与测试样本的真实故障类型(或工作状态)进行比较,计算正确率。J^JIJIIII原巧动信EMD/EEMD分提取特征SVM故障模式—^ ̄^ ̄^^号解向量分类器识别图4-6故障诊断流程图4.5基于峭度指标的滚动轴承故障诊断4.5.1基于EMD峭度的滚动轴承故障诊断按照故障诊斯流程,先对原始信号分别进行EMD分解,然后分别计算各IMF分量的峭度值-此作为特征向量进行训练与故障识别。分类识别结果如图4。,W7具体分类正确率如表4-3。从表中可W看出,内圈故障中,对早期故障识别正确率达到100%,对直径14mil点蚀直径的故障识别正确率只有55%,相对于严重故障的80%而言比较低。这可能与之前时域波形分析时发现的波形特征—振动频率成分减少相吻合。,早期与严重故障的诊断正确率均达到100%外圈故障中,而处于14mil的中度故障27 华北电力大学硕±学位论文正确率为65%,总体来说,对外圈故障识别效果较好。滚动体故障中11组21mil,只有严重故障的正确率达到90%,其余均较低,尤其是的正确率仅为20%,严重影响了滚动体故障整体的识别率。7。正常状态正确率为7.5%,整体正确率为76.15%,识别效果还有待提高表4-3基于EMD峭度的分类正确率S类别号正确数量总数量正确率(%)13140而2M20m内圈3112055故障4182090516208062020m外關7132065一故,障820201009102050滚动体10172085故障114202020901218整哥m7^mmm12心T711斗1|■III,;I絕!^I^I1IIII)1---i-■1>10修牛旬曲醒;IIIiij+++固-料*+I;;i1.IIt占IIIIIQILL.i—iGiSHBBEIIIjIIIIIIII---半msm夺新坤II1鲁聲^!!Itll1i——6^主1【和J霖絲;;;;!;命Immmi普种*’IIIIII4-^:I0实,I;际类别;I*;诊断类别1--2msm111JIIIIIIIIIII>II夺I抑术丰IiIII11II!]ItIIQ10酣100150200250300诊断样本图4-7基于EMD峭度的诊断分类图28 华北电力大学硕止学位论文4.5.2基于EEMD峭度的滚动轴承故障诊断EMD分解效果不尽人意,现在按照故障诊断流程,先对原始信号分别进行EEMD分解,然后计算各IMF分量的峭度,W此作为特征向量进行故障识别。分类识别结果如4-8-图。具体分类正确率如表44。表4^4基于EEMD峭度的分类正确率^类别号正确数量总数量正确率(%)S1M401OT240Yq100内圈3202050故障41020905182日1006182090外圈742020故障82020100 ̄ ̄916^io滚动体01520751故障H42015112202000整体^m78.4612wmmTT,1,广IIIIIIIIW来术i夺翊?陵IIIII>II1-■--■10^<tIIII来翊戲WIJI;;■..ItI<.t■8■-…■*f寺!:\IIItI■II牛韦:补I! ̄1----H-61齋msm:J獄|;;!:;;Iflm-IIiB化株特夺IIIIIII(1t4——恥—--靜帝::1 ̄j;;Iii际类别::*诊类别;!:III2--Fii1III11II'i《画田I奇啦I林I本IIIIIJIIIIillIII0I1[II050100150200250300诊断样本图4-8基于EEMD峭度的诊断分类图29 华北电力大学硕±学位论文从表中可W看出,内圈故障中,对早期故障与严重故障识别正确率均达到100%,对直径14mil点蚀直径的故障识别正确率只有50%,而21mil严重故障的正确率为90%。外圈故障中,严重故障的诊断正确率达到100%,早期故障正确率也在90%,而处于14ma的中度故障正确率仅为20%,严重影响了外圈故障整体的识别率。1,只有严重故障的正确率达到100%mil滚动体故障中,其余均较低,尤其是11組2的正确率仅为15%。正常状态正确率为100%,整体正确率为78.46%,识别效果还有待提高。4.5.3分析与讨论如表4-5。,基于EEMD分解的峭度指标故障诊断整体正确率略高于EMD分解针对正常信号的诊断,EEMD为100%,EMD为77.5%,相比之下,W峭度为指标,EMD分解存在将正常信号错误识别为故障信号的风险。在内圈故障中,对第3类的识别正确率均在50%左右,说明内圈故障14niil时,振动的不显著影响了故障的诊断效果。,早期与严重故障的正确率相近,在外圈故障中,但EEMD对14mil的正确率仅为20%低于EMD的诊断正确率,说明EEMD在外圈故障诊断中没有优势。在滚动体故障中,严重故障正确率相近,EEMD对早期故障的正确率较高,但是EMD和EEMD对21mil的正确率均很低,说明该方法对滚动体故障21mil不敏感。总体来说,EEMD在整体正确率有优势,但个别故障诊断效果并不好。表4-5比EMD与EEMD分类对—S类别号EM6正确率(%)EEMD正确率(%)E#I而100 ̄r"m内圏35550故障49090580100 ̄6%外圈76520故障8001001-——滚动体108575故障1120151290100-.整体'KlS7^考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的总正确率,如30 华北电力大学硕±学位论文-6表4。表4-6重新统汁后的分类对比SEMD正确率(%)EEMD正确率(%)3Ws1W内圈82.587.5外圈88.3370滚动体拍.591.25整体77.6986.54EMD与EEMD整体识别正确率均有提高,EEMD高于EMD,但是EEMD在外圈故障识别中,正确率低于EMD,仅有70%,而在内圈故障与滚动体故障识别正确率相化于EMD有明显优势,且不存在将正常信号误判的可能。。峭度指标作为单特征参数在进行故障诊断时,诊断效果有待于提离4.6基于均方根值指标的滾动轴承故障诊断对基于峭度指标的故障诊断分析后,进行基于均方根值指标的故障诊断分析。4.6.1基于EMD均方根值的滚动轴承故障诊断按照故障诊断流程,先对原始信号分别进行EMD分解,然后计算各IMF分量的均方根值-,W此作为特征向量进行故障识别。分类识别结果如图49。具体分类正确率如表-47。表4*7基于EMD均方根值的分类正确率?类别号正确数量总数量正确率(%)1m^ioo2MM100内圈3182090故障4202010052020100 ̄619M外圈7192095_故障820201009142^而滚动体10172085故障11142070122020100整体241m91^31 华北电力大学硕上学位论文————iS12BSBMB1iII了j--■---jeSaWIKBIIi()III)#-4-11马<410'I''I普扣阵8^——f7|jaartwfe^I林III1---1^!*霞0mam’爭111,,IIII-4J:j!aawaffe牛*III拂III!2111i!:o实际类别:i1=+诊断类别IIIjI1''〇l11i0如1001甜200250300诊断样本图4-9基于EMD均方根值的诊断分类图-,从表47中可W看出,内圈故障中其,对不同故障程度的诊断巧别正确率均很高中对直径14mil点蚀直径的故障识别正确率为90%,相对其他略低。外圈故障中100%,故障早期和中期的诊断正确率,严重故障的诊断正确率均达到为95%,整体对外圈故障诊断效果较好。滚动体故障中,只有严重故障的正确率达到100%,其相对较低,尤其是9组和11一;组的正确率为70%,定程度影响了滚动体故障整体的识别率。正常状态正确率为100%,整体正确率为%.76%,识别效果相比于峭度指标有了非一定实用价值常大的提高,说明基于均方根值的故障诊断方法具有。4.6.2基于EEMD均方根值的滚动轴承故障诊断EMD分解的结果己经有很大提高,接下來进行EEMD分解。按照故障诊断流程,先对原始信号分别进行EEMD分解,然后计算各IMF分量的均方根值,W此作为特征-向量进行故障识别。分类结果如图4-10。具体分类正确率如表48。-8,,其从表4中可W看出,内圈故障中对不同故障程度的诊断识别正确率均很高中对直径14m85%,相对其他略低。il点蚀直径的故障识别正确率为1,外圈故障中,早期故障的诊断正确率均达到00%故障中期和严重故障的诊断正确率略低。4m28m1,7mil85%21mil的正滚动体故障中1il与il的正确率为00%的正确率为,确率仅有55%,诊断效果有待提高。3外圈故障效正常诊断全部正确,整体正确率为9.08%,说明该方法诊断内圏故障与32 华北电力大学硕±学位论文果较好,滚动体故障诊断效果相比其他两类故障正确率较低。表4-8基于EEMD均方根值的分类正确率ii类别号正确数量总数量正确率(%)14040iw2M20ITO内圈3172085故障4202010052020100 ̄6^20外圈7192095心。*故賠8182090917WS5滚动体102020100120放障11155122020100整体^^93.08 ̄--_12mmm-r了半1,,IIitI■I半本杏翔田Bi吸辛tIIIIIIII:.如+.、iim110irfIIIIII种誦篇错牛III跨II—一-8一f|jIIIIII>msmIII米I!1!--6nmmm^號靜111I!:I:I■mwIIIL一4—-一一__;tIIIIImmmII丰;料寺I;I2msssm-11111_?;I;II〇实际类别IIs?*?*诊断类别IIIiIl1111I〇1050100巧0200250300诊断样本團4-EEMD均10基于方根值的诊断分类图4.6.3分析与讨论-9EEMD分解的峭度指标故EMD分解。如表4,基于障诊断整体正确率略高于针对正常信号的诊断,EMD与EEMD均为100%正确判断,不存在将市常轴承诊断为存在故障的可能。33 华北电力大学硕±学位论文mil诊断正在内圈故障中,两者趋势相同,巧为14确率略低,其原因与故障信号自身特点有关。在外圈故障中,EMD与EEMD呈现了不同的特点,EMD对严重故障诊断正确率为100%,其他略低,对,。说明;EEMD正相反早期故障的诊断正确率为100%其他略低EEMD对滚动轴承外圈早期故障诊断更有优势,在故障程度判断中,二者可W实现互补。在滚动体故障中,二者趋势接近,但EEMD对21mil故障诊断正确率低于EMD的。效果,说明对该故障程度信号不敏感,但其他程度均有所提高总体来说,EEMD在整体正确率有优势,但个别故障诊断效果并不好。表4-9EMD与EEMD分类对比0S类别号EMD正确率(/。)EEMD正确率(%) ̄S1002内圈39085故障41001005100100695100外圈79595故障8100909TOS5滚动体1085100故障11705512100100技体92^93]^考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的正确率,如表-410。表4-01重新统计后的分类对比HiEMD正确率(%)EEMD正确率(%)M100m内圈97.596.25外圈98.巧95滚动体9095957796整体..1534 华北电力大学硕±学位论文-EMD正从表410中可看出,EMD与E确率均有提島,EMD在内圈故障,其中与外圈故障诊断中更有优势,EEMD在渡动体故障诊断中正确率较髙。在重新统计过程中发现,11组正确率较低是因为部分被识别为9组,因此从故障种类诊断角度来看,EEMD更有优势。综上所述,基于EEMD均方根值的故障诊断方法是可行且有效的,EEMD在整体故障诊断中具有优势,但仍有提鳥空间,尤其是滚动体故障。4.7基于时域多参数指标的滚动轴承诊断一按照单时域参数指标进行分析对比之后,,考虑将W上两种时域参数指标进行结合用综合指标进行故障诊断,并将EMD与EEMD的诊断结果进行对比分析。4^.1基于EMD时域多参数的滚动体故障诊断按照故障诊断流程,将之前EMD分解后,各IMF分量求得的峭度与均方根值前后--排列作为特征向量进行故障识别。分类识别结果如图411。具体分类正确率如表411。表4-11基于EMD时域多参数的分类正确率类别号正确数量总数量正确率(%)Si40^220^100内圈32020100故障4202010052020100 ̄6^^1TO外圈7162080山。立故障82020100917^S5滚动体10162080故障11132065122020100整体242m9^一-从表411中可W看出,内圈故障中,诊断正确率均为100%,与単参数的诊断结果相比,有了很大提离,可见時域多参数法针对内圈故障诊断是十分有效的。外圈故障中,早期和严重故障的诊断正确率均达到100%,14mil故障正确率为80%,相比较低一,与单参数结果相比,高于峭度指标而低于均方根值指标。滚动体故障中,严重故障的正确率为100%,11组21mil故障正确率为65%,但相一比于单参数的诊断结果,有了很大提高。7mil与14mil的正确率同样介于峭度指标与均方根值指标之间。35 华北电力大学硕-上学位论文12-—T,]T'|IIIIt-■4-:II*IIIIIIIIItII—--14+1<10IIIIII)IIII本縫圓料卑I^;II;IIIIII8—f1I1IIIIIIIi'IIws^卡:味!IIIIIIIi-+^msm:;憂61W\IIIII1till4--711}[mmmI;!II!I2了1111;〇实际类另ii;iI^;;:*S品类另ij11|I0—*—■^'050100巧0200巧0300诊断样本阁4-11基于EMD时域多参数的诊断分类凶整体正确率为93.08%,诊断效果与均方根值接近,部分故障诊断有提升,也有部分故障诊断有所下降,原因就在于峭度指标的诊断效果不好,从而降低了多参数的整体诊断正碗率。4.7.2基于EEMD时域多参数的滚动体故障诊断分析完EMD诊断结果之后,按照故障诊断流程,将之前EEMD分解后,各IMF分量求得的峭度与均方根值前后排列作为特征向量进行故障识别。分类识别结果如图4--12。具体分类正确率如表412。-1295%可W看到,内圈故障中,从表4,14mil故障正确率为相比峭度指标有了很一大提高比均方根值有,。,相定提高总体针对内圈故障的诊断是有效的表4-12基于EEMD时域多参数的分类止确率Si类别巧化确数M总数结正确率(%)1^401002^^内圈3192095故障42020100520201006WM外圈7182090故,障A818209036 华北电力大学硕壬学位讫文表4-12(续表)S类别号正确数量总数量正确率(%)9142070滚动体10182090故障11112055122020100 ̄整体2^9154.-12mm)「rT,1s11IIII1^11I#詞warn:I1"10味+'-MHKtII1It——*--—-—8'Wsss^IjI\!mmm1\I1--愚6爭11J粉IIIIIIIill)IIIIIj*I?4]--MWB2)71111'''〇实际类别nagguaamI|+诊巧类别I;;;;L一.1101'0知1朋甜2朋2胡3001诊断样本图4-EEMD时域多参数的诊断分类12基T幽外圈故障中,其诊断正确率趋势与均方根值指标相同,针对早期故障较为有效,总体圧确率介于峭度指标与均方根值指标之间。滚动体故障中,严重故障正确率为100%,11组21mil故障最低,与均方根值指标持平,高于峭度指标诊断结果。9154整体准确率为.%,介于均方根值指标与峭度指标之间。4.7.3分析与讨论4-MD,EMD如表13基于分解的时域多参数指标故障诊断整体正确率略高于EE分解。针对正常信号的诊断,EMD与EEMD均为100%正确判断,不存在将正常轴承诊断为存在故障的可能。在内圈故障中,EMD诊断正确率全部达到100%,而EEMD的14mil诊断正确率略低。在外圈故障中,EMD与EEMD呈现了不同的特点,EMD对14mil故障诊断正确率为80%,EEMD则对14mil、21mil故障诊断准确率均为90%。在滚动体故37 华北电力大学硕±学位论文障中,二者趋势接近,但EEMD对7miK21mil故障诊断正确率低于EMD的效果。表4-13EMD与EEMD分类对比S类别号EMD正确率(%)EEMD正确率(%)止常I2100内圈310095故障410010051001006外圈78090故障8100909EW滚动体108090故障11655512100100整体91089L54考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,如表,得到故障类型分类的正确率-144。-:表414重新统计后的分类对比 ̄' ̄■SEMD正确率(%)EEMD正确率(%)! ̄'M100内圏10098.75外圈93..339667滚动体97..59625整体97.的97.69从表4-14可W看出,EMD与EEMD的诊断正确率都有了较大提高,EMD对内圈故障与滚动体故障的诊断正确率均高于EEMD,而EEMD表现出的是各种故障诊断的。均衡,没有很低的数据一综上,基于时域多参数的滾动轴承故障诊断的结果与单特征参数故障诊断相比,总体有提升,但由于哨度指标的诊断效果并不理想,所W也影响了多参数指标的诊断效果。在EMD于EEMD诊断的比较中,EMD在部分故障种类诊断中优于EEMD诊断效果,但整体上基于EEMD的诊断结果更为均衡。4.8本章小结一本章提出了种基于EEMD分解—,W时域特征中的两个参数峭度、均方根值为38 华北电力大学硕±学位论文特征量进行滚动轴承故障诊断的方法。并对EMD和EEMD分解的诊断结果进行对比,EEMD的诊断正确率要高于EMD,且将均方根值作为特征参数的诊断正确率要远高于一将峭度作为特征参数的正确率,证明均方根值是个有效的时域特征参数。而将两个时域参数组合成多参数诊断方法,其诊断正确率整体有所提高,部分故障诊断正确率介于哨度指标与均方根值指标之间。39 华北电力大学硕±学位论文第5章基于EEMD频域特征参数的滾动轴承故障诊断5.1引言一上章中选取时域特征中具有代表性的两个特征参数作为指标进行故障诊断,证明了时域参数的有效性。特征参数法不仅包括时域恃征,还包括频域特征,本章将重也频率、巧方根频率作为频域指标进行故障诊断。5.2信号频域特征。在滚动轴承故障诊断中,频域特征往往会包含丰富的故障信息相对于时域特征,在信号成分较为简单故障中,可W直接从频谱中找到故障频率。而相对复杂的信号,如一调制信号,或者信号中包含了大量的噪声些处理,,将故障信号掩盖,则通过可W还原故障信号,或者暴露出典型的故障频率特征。5.2.1原始信号频域特征(1)不同故障类型频谱图本章基于滚动轴承振动信号的频域特征参数为依据进行故障诊断,原始信号中,正常与故障信号的频谱图是不同的。选取与第四章中正常与故障时域波形所对应的频谱图进行分析-、。如图51中分别是正常内圈故障、外圈故障、滚动体故障的频谱图。从图5-。正常信号的频1中可W看出,正常信号和故障信号的频谱是有很大区别的率集中在2100HzW下,且幅值较小,低频处幅值较大的是轴转频率,其中2100Hz是lOOOHz处的谐波。0-4000HZ均有内圈故障信号相比于正常信号,増加了高频成分,其频率从分布,4-5,在2000Hz到4000Hz中较为集中表现出冲击性而且幅值比正常信号大了倍,峰值出现在3500Hz附近。外圈故障同样表现出较强的冲击分布在2000Hz到4200Hz之间,峰值出现在3000Hz附近,峰值大于内圈故障峰值,同时,低频处(2000HzW下)则没"有较强烈的频率分布,这点与内圈故障不同。滚动体故障信号则从04000Hz均有分布,-4000HZ且主要集中在2000,峰值出现在%OOHz附化但是幅值明显小于内圈故障与外圈故障,仅仅略大于正常信号,低频处在正常信号中出现频率被抑制了。通过简单的计算,可M得到滚动轴承常见故障的故障频率,但是送些都属于低频,!由于噪声信号。、信号调试等原因,无法直接从频谱图中获得较为明确的故障频率信息(2)不同故障程度的频谱图与第四章相同,W内圈故障为例,数据选择点蚀直径依次为7mil、14mil、21mil、28m-il,通过傅里叶变换,画出对应的频谱,如图52。40 华北电力大学硕上学位论文正常信号頻谱0.1IIIII0--.05醒Qak<l^---.i-■J01000200030004000加006000内圏故障信号频谱Q——-■*■-3IIII0.2迴阻0-—.1—…00100020003000400050006000外圈故障信号频谱0■ ̄■ ̄- ̄.41i1rj。-'2謹y|0100020003000400050006000滾动体故障信号频谱01■ ̄ ̄ ̄■ ̄ ̄ ̄1I.11[化。5--I1||01000200030004000加006000频率Hz圈5-1滚动轴承正常与故障的频谱对比通过对比可W看出一7m,故障程度不同,频谱图也表现出定差异。点蚀直径il与14mil的频谱成分比较接近,在低频处有频率分布,不过幅值较小,在2700Hz与3500此附近有较大能量,但I4mil的幅值远小于7mil,这与第四章中时域特征也是相吻合的,14mil的时域技形幅值比7mil小。21mil的频率集中在3000Hz附近,分布范围比7mil与14mil要小。28mil的频谱中,从500Hz到3500Hz均有较为强烈的频率分布,且幅值^大了很多,其中峰值出现在600Hz附近。41 华北电力大学硕dr学位论文点蚀直径为7mil0.2III1I-0—碧.10100020003000400050006000点蚀直径为14mil0.06riiiii0-’04^11ill010002000300040加加006000点蚀直径为21mil'^0.3iIIIrj0-.2!*"UJii.」—01000200030004000加006000点蚀直径为23mU0.4i111l0.■^政tlj心0化0020003000400050006000频率Hz.图5-2滚动轴承内圈故障不同点烛直径的频谱形对比5.2.2EEMD分解的信号频域特征与第四章相同,,W内圈故障为例,分别进行EMD和EEMD分解并对各分量进行-,画出频谱。EMD分解后IMF分量频谱如图53EEMD分解后IMF分量傅里叶变换,5-4如图。()EMD1分解后的信号频域特点-3,5,进行EMD分解后从高如图,得到8个MF分量,IMF按照不同频率尺度频到低频依次分解。IMF1分量主要包含2500Hz到3500Hz的频率,而IMF2则包含0到2000Hz的频率,频带太宽,IMF3主要频率为400化到80胤Z,显眼与IMF2有重叠,42 华北电力大学硕J:学化论义也就表明400Hz到800Hz频带被分在了两个IMF分量中,送就是模态漉叠现象。04i.IITi'含‘IJJJAAk諭…■■I,Q01000200030004000500060000.02..,,S0.01I含0,,,0100020003000400050006000— ̄0.1i11—0.05L函'-一一4如?-I■QII0100020003000400050006000— ̄—0.02111〇-.圓.1.北版..01000200030004000500060000.02'.0—臟.01101000200030004000500060000'111.01-0.005IH含ft、IIIIIQ01000200030004000加006000-3X1〇^11111'--麵1II.>\\^IQII.0100020003000400050006000-3X1〇——^■■——4\11I00jPh勺-§\V—QIIII0100020003000400050006000-3X1〇——4!r[I苗2--么\、瓜…j—j^_0100020003000400050006000阁5-3内圏故障EMD分解频谱43 华北电*大学硕±学位论文-■— ̄ ̄ ̄ ̄ ̄0II.402--II0i0100020003000400050006000———^11r10.024^——麻的U-'K汰—.,__〇i0100020003000400050006000-——0rn11.05j20IrII01000200030004000500060000州1ui01000200030004000如006000———1r0rr1I.010--I.005一''WIKk-.1]1—■0I0100020003000400050006000-3X1〇21广i1§fLI^4^.二=_=.三」Q^0100020003000400050006000-3X1〇■41r111-麵^—IOVI0100020003000400050006000-3X1〇j ̄ ̄II4II〇〇-I\V-IIIIQ010002000300040005朋06000.。0.01i10-g.005心iII1Q_J01000200030004000日00日6000-^^0.02、Iwig0.01Ljt--I1.-QI1I0100020003000400050006000图5*4内閣故陈EEMD分解频谱44 华北电力大学硕±学位论文(2)EEMD分解后的信号时域特点由图5-4可看到I,EEMD分解后,得到9个分量,在MF2中,lOOOHzl^下的频率明显少了,与IM巧的频带划分更明显,当然IMF3与MF4依然存在部分模态混叠现象,不过由于主要的故障信息包含在前面的IMF分量中,模态混叠现象的影响要比EMD小一些。同时EEMD分解的分量比EMD更多,也就代表在频带划分上,EEMD更细致。5.3基于重瓜频率指标的滚动轴承故障诊断5.3.1基于EMD重心频率的滾动轴承故障诊断按照故障诊断方法流程,先对原始信号分别进行EMD分解,然后计算各IMF分量的重也频率。分类识别结果如图5-5。具体分类正确,W此作为特征向量进行故障识别5-率如表1。-70%从表51中可看出,内圈故障中,对直径14mil点蚀直径的故障识别正确率为,对其他故障程度的诊断识别正确率均很高。.外圈故障中,14mil故障的诊断正确率均达到%%,早期故障的诊断正确率为85%,严重故障的诊断正确率较低,仅有55%。滚动体故障中,只有严重故障的正确率达到%%,14mil故障诊断正确率为90%,其余故障诊断正确率相对较低,尤其是早期故障,其诊断正确率为70%。正常状态正确率为100%,说明不存在将正常信号错误诊断为故障的可能。整体正确率为87.69%,仍有提高空间。5-EMD表1基于重也频率的分类正确率类别号正确数量总数量正确率(%) ̄M14040100220201TO.内圈3152075故障419209552020100 ̄616M80外圈1192095A化故障8112055914M^滚动体10182090故障1117208512192095整体m2m87.6945 华北电力大学硕±学位论文mm—^—?12rT1,ss),IIIIII丰IIII半牛利坦IIIIIII1IIIII??-)+1<1041IIiIIII'(I!^^I;;!:■—一-^争_^寺爭87斗f!:\IIIIIILijgmlMtfo1IIIIIIIIII-1-^1------!巧6^I1卡11狱;;I1;1iI1(IIIIIIIII1i—4-寺本-1t[IIIImmm半琳?普牛牛|||I|牛;1IIIII2了1!11 ̄ill!〇实际类别I__1;;;I1+诊断类别1111I111I〇1110如100巧02002如300诊聯样本圏5-5基于EMD重也频率的诊断分类图5.3.2基于EEMD重心频率的滚动轴承故障诊断按照故障诊断方法流程,先对原始信号分别进行EEMD分解,然后计算各IMF分-、56。量的重屯频率,W此作为特征向量进行故障识别。分类识别结果如图具体分类正5-2确率如表。表5-2基于EEMD重也频率的分类正确率S类則号正确数量总数量正确率(%)M140402^M内圈3172085故障42020100520201006^20外罔72020100故障8152075918M^滚动体10182090故障11172085122020100整体245m94^46 华北电力乂学硕上学位论文■.‘10___…'么____巧_I了II]^-eI牛MIII普I^一;+-4--mmssm10Ia'''>夺irmmm辛一一一-1-18却—本i1m/t^^I1II1--6i!1J蘇1撕;IIII印肺酬:II11‘4WSSSS^:jIjmmm林〇I:;;;121111〇实际类别ji11111!+诊断类别I;''〇i^'0501001凯200250300诊断样本图-6、5基于EEMD重屯频率的诊断分类图5-2可W看出14,对直径mil点蚀直径的故障识别正确率为80%,通过表,内圈故障中对其他故障程度的诊断识别正确率均达到100%。外圈故障中,714100%mil和mil故障的诊断正确率均为,严重故障的诊断正确率较低,为75%,有提高空间。滚动体故障中,21mil故障的正确率为85%,7mil和14mil故障诊断正确率为90%,严重故障诊断全部正确。正常状态正确率为100%。整体正确率为94.23%,基于EEMD重也频率的故障诊断一方法是有效的,故障诊断正确率比较均衡,外圈严重故障诊断需要进步提高。5.3.3分析与讨论如表5-3EEMDEMD,基于分解的重也频率故障诊断正确率高于分解在内圈故障中,EEMD故障诊断术确率全面高于EMD,其中14mil故障正确率略低,其原因与信号本身有关。在外圈故障中,EEMD故障诊断正确率均高于EMD,但在外圈严重故障诊断中,仍有提局空间。在滚动体故障中,EEMD故障诊断正确率较为均衡,均在85%及W上。其中21mil故障诊断正确率较低,这与时域特征参数诊断的结果类似。综上,基于EEMD重也频率的滚动轴承诊断方法针对故障严重程度的诊断是有效的,但其中外圈严重故障诊断有必要进一步提高。47 华北电力大学硕±学位论文5-3表EMD与EEMD分类对比〇类别号EMD正确率(%)EEMD正确率(/〇)^\2内圈37585故障4951005100100680m外圈795100澈膺855759TO^滚动体109090放障11858512%100整体8^9423考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的正确率,如表5 ̄4。表5^重新统计后的分类对比〇EMD正确率(/〇)EEMD正确率(%)'100.内圈92.596.25外圈78.巧91.67滚动体9097.5整体89.姑96.15由表5-4可^看出,EMD与EEMD故障诊断整体正确率均有提高,EEMD高于EMD,!达到了96.15%。外圈故障诊断正确率略低与其他故障类型。EEMD在W重也频率为指祿的故障诊断中相比于EMD有明显化势。重也频率指标作为单特征参数在进行故障诊断时,不仅,诊断效果明显可诊断故障类型,同时可W诊断故障严重程度。5.4基于均方根频率指标的滚动轴承故障诊断5.4.1基于EMD均方根频率的滚动轴承故障诊断按照故障诊断方法流程,,先对原始信号分别进行EMD分解然后计算各IMF分量-的巧方根频率。57。具体分类正,W此作为特征向量进行故障识别分类识别结果如图48 华北电山大学硕i-学位论文确率如表5-5。表5-5基于EMD均方根频率的分类正确率类别号正确数量总数量正确率(%)140402^M内圈3132065故障419209552020100 ̄ ̄618^W外圃7192095故障-8112055913^^滚动体101920%故障1115207512192095整体^^86.92一—121J1IItIIIII11I夺IiI鲁i翔閣巧,I.1IIIIIIII1:1I+?+->——10msmm,<?IIIIIIII[IIIjBgTWr,II:^B?林IItjII'8--"^——f本Trmms^,4^{IIIIIIIIIIIIIIIII—-—-驚6术;11II獄I?II?IIIII),Ifl—M:II!IIIIII(L_-一4细,11jIIIIII'mm制^h^夺m中半III味制IIIT21半110;’iwe1I;1I+诊断类别;I;;;1〇1111110加100150200250300.诊断样本图5-7基于EMD均方根频率的诊断分类图从表5-5中可W看出,在内圈故障诊断中,14mil故障诊断正确率仅为65%,其他故障程度较高。在外圈故障中,严重故障诊断正确率仅有55%。在滚动体故障诊断中,14mil与28mil故障诊断正确率可W达到%%21ml,7mil与i故障诊断正确率较低。整〇一体正确率为86.92/〇,但基于EMD均方根频率的故障诊断正确率较低,需要进步提高。49 华北电力大学硕±学位论文5.4.2基于EEMD均方根频率的滚动轴承故障诊断按照故障诊断方法流程,先对原始信号分别进行EEMD分解,然后计算各IMF分-8量的均方根频率,W此作为特征向量进行故障识别。分类识别结果如图5。具体分类正确率如表5-6。--表5-6基TEEMD均方根频率的分类止确率^类别号正确数豈总数量正确率(%) ̄正常1404021820%内圈3192095放障419209552020100619^^外闡心72020100故P早8182090917W85滚动体10182090故障11152075122020100' ̄' ̄盤体mm93.46■_占--12rT,1|IIII!IIII1I特幸命■II1ii—*一一V?+i101(ggBHfe4III1}I夺爭I利悼I:IIIttII8——f;!II1IIIIII--—-1---——-6I1胃钱J11i::;1-IttI^SSSSSI'^—一4— ̄i-h——一一山—HJI*II??1aiuMXiL\f||y市IIII)-11■211i'〇实际卖I;;;I4诊断巧别!1111IQI1Ii0加100150200250300诊断样本图5-8基于EEMD均方根频率的诊断分类圏50 华北电力大学硕±学位论文由表5-6可W看出,在内圈故障中,基于EEMD均方根频率的故障诊断正确率均达到90%及上,说明针对内圈故障可W实现较好的诊断。在外圈故障中,诊断正确率也达到了90%及W上。在滚动体故障中,21mil故障诊断正确率为75%,是同种故障中最低的。5.4.3分祈与讨论-7EEMD分解的均方根值指标故EMD分解如表5,达到,基于障诊断正确率高于了93.46%。相对于EMD在14mil内圈故障和21mil外圈故障的较低诊断正确率,EEMD的正确率均提升到一90%及W上。当然,EEMD在滚动体故障的诊断方面还需要进步提高。表5-7MD与EEMD分类对比ES类别号EMD正确率(%)EEMD正确率(%)2%内圈36595故障4959551001006%^外圈795100巧障85590厂^^.滚动体109590故障1175751295100''整体8^9146考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的正确率,如表5-8。表5-8重新统计后的故障类型分类SEMD正确率(%)EEMD正确率(%)100內圏9096.25外圈88.33%滚动体化75100整体撕.7797.6951 华北电力大学硕±学位论文EMD与EEMD整体识别正确率均有所提离,EEMD在各类故障类型诊断中正确率均高于EMD,整体故障诊断率达到了97.69%。可见,均方根频率指标作为单持征参数在进行故障诊断时,诊断效果非常明显。在故障类型诊断中,EEMD诊断效果比较好,但是在故障程度识别中,滚动体故障程度识别还有侍提高。因此,W均方根频率为指标的滚动轴承诊断方法是可行的。5.5基于频域多参数指标的滚动轴承故障诊断按照单一频域参数指标进行分析对比之后,,考虑将W上两种时域参数指标进行结合用综合指标进行故障诊断,并将EMD与EEMD的诊断结果进行对比分析。5.5.1基于EMD频域多参数的滚动轴承故障诊断按照故障诊断方法流程,将之前EMD分解后,各IMF分量求得的重也频率与均方-根频率前后排列作为特征向量进行故障识别。分类识别结果如图59。具体分类正确率如表5-9。表5-9基于EMD频域多参数的分类正确率类别号正确数量总数量正确率)(%Mi404022020内圈3142070故障418209052020100 ̄16^80外圈7192095故障8112055,913W65滚动体10192095故障1114207012182090整体瓦im8538一由表5-9,在内圈故障中,频域多参数相比单参数断正确率基本持可W看出,诊平一,其中,21mil故障诊断正确率低于单参数,14mil故障诊断正确率介于重也频率与均方根频率之间。一一参数的故障诊断结果在外圈故障中,趋势与单致,严重故降诊断正确率最低,仅为55%。其次7mil故障诊断正确率较低,为80%,一在滚动体故障中,频域多参数的故障诊断正确率均低于单参数的诊断正确率。7mil与21mil故障诊断正确率分别为65%和70%。52 *华北电如火学硕I;学位论文虽然整体诊断正确率为85.38,但是基于EMD频域多参数故障诊断方法仍有待提高。12厂TMitmi1If|IIIIII料I辛1II牛^IIIII10私41^sssm-1i]IItIIIIItI半夺^辛命利株:I;I;IIIIII8——I1IIdTOaLmfe*IIIIIIIIIIII—+----——6;1!mamI普親靜IIIIIIIIIWSKBIIIIIIItIIIItIL4—wmm林—;1;I1tII,^mm特■牛如林?本蛛■牛!:I\:;IIIIIIL11T11;;IIo实际类别III.+诊断类别I;;;*'LoJ___L」1J050100150200巧0300诊断样本图5-9基于EMD频域多参数的诊断分类凶5.5.2基于EEMD频域多参数的滚动轴承故障诊断按照故障诊断方法流程,将之前EMD分解后,各IMF分量求得的重也频率与均方5-根频率前后排列作为特征向量进行故障识别。分类识别结果如图10。具体分类正确率5-如表10。表5-10基于EEMD频域多参数的分类正确率H类别号正确数量总数量正确率(%)140401TO2M20巧圈3182090故障419209552020100 ̄6M20外圈72020100丄故L障口廿8192095 ̄ ̄91620m滚动体10182090故障11182090122020100整体2^m953853 华北电力大学硕±学位论文5-1410可,mil故障诊断正确率介于重也频率与均方根由表W看出,在内圈故障中。频率之间,其他均有所提高,整体来说针对内圈故障诊断正确率较高一在外圈故障中致,7mil与14mil故障诊断,故障诊断正确率分布趋势与重也频率〇正确率为1〇〇/。,严重故障诊断正确率则高于重也频率的诊断正确率。、,但是其他故障在滚动体故障中,早期故障诊断正确率低于重屯频率与均方根频率程度诊断正确率都有所提高。基于EEMD频域多参数故障诊断效果较好,整体故障诊断正确率达到了95.38%,滚动体早期故障诊断正确率有待进一步提高。■12rT1,1I;II半II夺+10+1mmm1t夺1I鍾瞧料;^I]一——一一—一一一-,;,。8夺7!;f!11111--^---^1i^ssmmh叫v^6+1IfI1^Igjggijgj^L.411||J[|i言ijI-I1!211 ̄1!i1^实际类别!??*11*诊断类别fII;J'JL'0'0501001凯2002前孤0诊断样本图5-EEMD频域多参数的诊断分类图10基于5.5.3分析与讨论5-。如表11,基于EEMD频域多参数故障诊断正确率高于EMD针对正常信号的诊断,EMD与EEMD均为100%正确判断,不存在将正常轴承诊断为存在故障的可能。90%在内圈故障中,EEMD诊断正确率全部在及W上,可见对内圈故障诊断比较有EMD14。效,而的mil诊断正确率较低在外圈故障中,EEMD的故障诊断正确率远高于EMD,而且通过多参数结合,克服了严重外圈故障诊断正确率低的问题。在滚动体故障中,EEMD也高于EMD的诊断正确率,不过在滚动体早期故障诊断一还有待于进。中,步提髙总体来说,EEMD比EMD在故障类型与故障程度诊断中更有优势,诊断正确率也比较均衡。54 华北电力大学硕±学位论文表5-11EMD与EEMD分类对比S类别号EMD正确率(%)EEMD正确率(%)Mim2iTOiro内圈37090故障490955100100 ̄6^外圈195100化n立故障855959?SO滚动体10%90故障1170902190100整体85.3895.38考虑到SVM在识别过程中存在误差,因此,在仅考虑进行故障类型识别的情况下,可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的正确率,如表5-12。5-表12重新统计后的分类对比WMEMD正确率(%)EEMD正确率(%)内圈9096.2583.33%外圈.33滚动体91.25100整体90.38W.46EMD与EEMD整体识别正确率均有所提高,EEMD在各类故障类型诊断中正确率均高于EMD。可见,基于EEMD重也频率和均方根频率的频域多参数指标的故障诊断方法在诊断识别故障类型时,具有很高的诊断率,而且在进行故瞳程度识别时,也能达到较高的正确率一,尤其是内圈故障与外圈故障,滚动体故障诊断还有定的提升空间。5.6基于时域频域多参数指标的滚动轴承故障诊断一参数及多参数故障诊断分析后经过时域与频域的单,其诊断结果比较满意,但是仍有提升空间,尤其是故障程度诊断方面。现考虑将时域参数与频域参数相结合,组成时域频域多参数特征向量,进行故障诊断,并巧EMD分解和EEMD分解诊断结果进行化较。55 华北屯为乂学硕:I学位论文5.6.1基于EMD时域频域多参数的演动轴承故障诊断建立时域频域多参数特征向量,考虑到运算效率,从时域特征参数与频巧特征参数一中只能各选取个特征参指标,频域中重必频率效果较,时域中均义报值指枯效果较好好,将这两个特扯参数W结台。按照故蜡诊断流粒,将之前EMD分解后,各IMF分M求得的均方根值与重必频率前后排列作为特化向量进行故障识别。分类识别结果如图--511。具体分类正确率如表513。表5-EMD13基于时域频域多参数的分类正确率S类别号止确数结总数量化确率(%)140^ ̄2WM100内圈3182090故障4202010052020100 ̄6^M外圈72020100故、,障一8202010091320?滚动体10172085’故巧1;182090122020100整体mm94^-12f了\11II.II1I--.IiMMfc11辛I卡味和本iB1IIIIIIIIIIII10i41(smmm1<I1IIIII,msm:市ItIIi1:L——1.riMMMB!1fwt'It)IIIIIIIIIIIIIvpWBmBItIII万IIIIII---1,1卞?1絲tT;6讯怖mh;;;:::I>iIIIItII-I1iIIiI1一J—____—____—________——__J?..tI1IIJI>iIIKI山I?山Ij卡牛十II1,IIIIIII(2T巧HHK1111〇实际类男|J:iiiI|*诊断类别;:;I;L_i!0J!!10如100巧0200250300诊断样本5-圍11基于EMD频域多参数的诊断分类圏56 华北电力大学硕±学位论文-3m从表51中可W看出,内圈故障中,14il故障诊断正确率为90%,相比之下略低,一其余故障程度诊断正确率均为100%。与单参数的诊断结果相比,与均方根值的故障诊断正确率相同,比重也频率的故障诊断正确率有了很大提离。外圈故障中,所有故障的诊断正确率均达到100%,经均方根的配合,克服了重必频率在严重故障诊断中正确率太低的问题,时域频域多参数对外圈故障诊断达到了最佳效果。滚动体故障中,早期故障的诊断正确率为65%,低于均方根值指标和重也频率指标的诊断正确率,这应该是时域频域参数在组合时,排序方式导致的。14mil与21mil故障诊断正确率介于均方根频率和重也频率之间,总体来说,滚动体故障诊断效果相比于一单参数并不具有优势。整体诊断正确率为94.62%、,可见,时域频域多参数对内圈外圈故障诊断非常准确。5.6.2基于邸MD时域频域多秦数的滾动轴承故障诊断按照故障诊断流程,将之前EMD分解后,各IMF分量求得的均方根值与重也频率前后排列作为特征向量进行故障识别-。。分类识别结果如图512具体分类正确率如表5-14。表5-MD时域频域多参数的分类正确率14基于EEW类别号正确数量总数量正确率(%)140402^^内圏32020100巧障4202010052020100 ̄6^Mm外圈72020100故障82020100 ̄^91520了5療动体102020100巧障11182090122020100整体^m9^5-从表14中可W看出,正常状态、内圈故障、外圈故障,都可W实现100%正确诊断,这是非常理想的诊断方法。滚动体故障中,早期故障是诊断准确率最低的,为75%,m一其次是21il故障诊断正确率为90%。与单参数故障诊断方法相比,不仅提高了均方根值为指标的故障诊断中严重故障的诊断正确率,也提高了重必频率指标的故障诊断中57 华北电力大学硕+学位论文14mil和21mil故障的诊断正确率,但也拉低了早期故障的诊断故障率,所W多参数法不仅仅是之前各项单一参数的优势叠加一,也会出现比单参数诊断正确率更的情况。————_ILI11I]IIIII'I辛SWBB1IIIIIiIItIII+———1wmm1<10iIIIIII111I'IImm牛牛'IIIIIIIJIIIQ———————一一―—―一————_———.IO-———―———————————————r了IIIIIIIIIIIIIIitIIIIIiI6——_--叫——————糕,十i1I;:;;;[IIWBSSBIIIIItII(4__—--USB-j;;1II1III1nagMX4IIII>jBMg|IIIIIIIIIIII2Ti-WHSKB^111i:::〇实1;脉类别___Ii'.止i;1+诊断类别:lIIIIIIjQ!I—丄I050100巧0200250300诊断样本阁5-12基于EMD频域多参数的诊断分类簡5.63.分析与讨论5-EEMD频域多参数故障诊断吊确率高于EMD如表15,基于。表5-15EMD与EEMD分类对比S类别号EMD正确率(%)EEMD正确率(%)^1Ito2内圃390100故除410010051001006外圈7100100故障*81001009?干5滚动体1085100故障11909012100100整体94^9^在内圈故障中,EEMD诊断全部正确,优于EMD的诊断结果。在外圈故障中,EMD58 华北电力大学硕女学位论文—与EEMD样全部诊断正确。在滚动体故障中,EEMD的诊断正确率也高于EMD的诊断正确率。,考虑到SVM在识别过程中存在误差,因此,在仅考虑进斤故障类型识别的情况下可W将错误的分类按照故障类型的大类别重新统计,得到故障类型分类的.巿确率,如表5-。165-表16重新统计后的分类对比miEMD正确率(%)EEMD正确率(%)^m100内圈97.5100外圏100100滚动体97.5100>8整体^.02100EMD整体识别正确率有所提高,EEMD在故障类型诊断中实现了没有错误。可见,基于EEMD时域频域多参数故障诊断方法可100%诊断故障类型,针对故障程度诊断时,内圈故障、外圈故障可W做到100%诊断,而滚动体故障程度诊断正确率也可达到很高的正确率。基于EEMD结合均方根值和重也频率的诊断方法是有效地,针对工业生产具有实用价值。5.7本章小结一EEMD分解一一重也频率本章提出了种基于,W频域特征中的两个参数、均方根频率为特征量进行滚动轴承故障诊断的方法。并对EMD和EEMD分解的诊断结果进行对比。其中,W重也频率为特征参数的诊断效果最化从诊断故障程度角度,识别內圈故障一、外圈故障的故障程度均在%%及^^上,滚动体故障程度识别还存在定误差。将重必频率与均方根频率结合的多参数诊断方法相比于单一参数的故障诊断方法均有巧高的诊断正确率。将时域特征参数均方根值与频域特征参数重也频率相结合的故障诊断方法可W取得同类型中最好的诊断效果。59 华电力大学硕±学位论文第6章结论与展望6.1结论滚动轴承是旋转机械中不可或缺的支撑部件,但是它的结构决定了其易损坏的特点,因此如何对滚动轴承的运行状态进行识别一,尤其是故障类型进行诊断,进步对故障程度进行诊断一,直是研究的重要方向。本文总结了滚动轴承的基本结构、常见故障类型、振动产生机理,为滚动轴承故障识别提供了理论基础。同时提出了基于EEMD时域特征参数和频域特征参数两个方面及将二者结合的滚动轴承特征参数提取方法,并通过实验数据进行了故障识别验证,主要结论如下:(1)时域、频域特征参数包含了滚动轴承故障信息,W此为特征参数进行提取可1^1有效地识别故障类型,同时可^判断故障严重程度;(2)通过频谱分析,可W观察到EMD分解存在的模态混叠现象,同时从原理及频EEMD—定程度克服模态海叠现象谱分析都验证了可W。在滚动轴承故障识别时,基于EEMD分解的故障诊断正确率高于基于EMD的故障诊断正确率;(3)采用基于EEMD时域特征参数提取与故障诊断方法,尤其W均方根值为特征参数,可W有效的对滚动轴承故障类型与故障严重程度进行诊断与识别。将两个时域特征参数结合的多参数诊断方法,可W达到较高的诊断正确率;(4)采用基于EEMD频域特征参数提取与故障诊断方法,尤其W重也频率为特征参数,可W有效的对滚动轴承故障类型与故障严重程度进行诊断与识别。将两个频域特征参数结合的多参数诊断方法,可W达到较高的诊断正确率;(5)采用基于EEMD时域频域特征参数提取与故障诊断方法,将均方根值与重也频率相结合的时域频域多参数诊断方法,可W达到更高的诊断正确率;一(6)在多参数结合进行故障诊断的过程中,可能会出现两种单参数的诊断优势叠加现象一,也可能导致多参数诊断法在部分故障类型诊断中诊断正确率低于单参数的诊断准确率。6.2展望滚动轴承的故障诊断技术发展对于保证旋转机械正常运行有着及其重要的意义。虽一定的研究成果间和精力有限一然本文取得了,但由于时,还有许多工作需要进步研究和发展:(1)局限于数据来源,本文仅论证了基于EEMD特征参数提取方法在诊断滚动轴承点蚀故障中的有效性,而滚动轴承故障类型并不局限于点蚀故障,还有疲劳剥落、磨60 华化电力大学硕±学位论文损等故障,因此,后续还需要通过其他故障类型数据来验证该方法的有效性及通用性;(2)本文中,利用多特征参数进行故障诊断,虽然诊断结果比较好,但仍有提高空间,考虑到不同的持征参数对不同故障类型更有效,今后还应该考虑多种(大于两种)一参数共同诊断识别,从而增加该方法的适用范围,进步提高诊断祀确率。(3)本文中,将原始信号进行EEMD分解后得到若干IMF分量,但众多分量中有一部分是无用的,甚至会干扰其他分量的特征提取,后续应当考虑通过多种手段,将无用的分量滤除,从而提高运算效率,并提高诊断正确率。61 华北电力大学硕±学位论文参考文献1.付新欣.基于小波分析的滚动轴承故障诊断技术研究[D].沈阳航空航天大学2013[],口]周福昌.基于循环平稳信号处理的滚动轴承故障诊断方法研究[D].上海:上海交通大学2006.,-敏:口黄邵毅.设备故障诊断手册机械设备状态检测和故障诊断[M],西安西安]徐,-1998530.交通大学出版社:,[4]离毅龙.数据挖掘技术及其在工程诊断中的应用[D].西安:西安交通大学,2000:-110.口]ZhaoMY,ZhuHLLocomotiveBearingFaultDiagnosisandPredictionMethod口].Locomo-tiveElectricTransmission20055:5860.,()工程机械状态检测与故障诊断技术化汁量与测试技20024-4术:403.阀谢宝义()[7]张家凡?振动信号的包络解调分析方法研究及应用[巧?武汉:武汉理工大学,2008.[8]Randall艮B.Hilberttransformtechniquesinmachinediagnostics[C].Internationalerence-ConfonRotordynamicsToko2006409420.,,y,9了康谢明张彼德等.复解析带通滤波器及其在解调分析中的应用m.振动工程[],,,20005-学报133:38392.,,()10ZhanF,DinK.AnalsisofGeneralizedDemodulatio打of化eThreeAlorithmsand[]ggygheL-timitationsoftheStudy口?JoiraialofVibrationEngineering2002152:243248.],,()朱善安..农业[u]武巧雷,,林瑞仲,等基于能量算子解调法的滚动轴承故障诊断阴200-化械学报334:118120.,,(1)[12]了康,李巍华,朱小勇.齿轮及齿轮箱故障诊断实用技术[M].北京:机械工业出版粗2005.,[13]柱秋华,杨署年.细化包络分析在滚动轴承缺陷诊断中的应用饥.轴承20043:,,()-3431.M-28锦泰小波分析导论[M].西安西安交通大学出版化1995:17.[]崔15张中民..,卢文祥,杨巧子等基于小波系数包络谱的滚动轴承故障诊断饥振动[],工程学报-1998111:6569.,,()16LiuHTianYFWuT.TurbineGovernorModelIncludinUnitCoordinatedControl[],,gSystemanditsApplication[J].AutomationofElecfricPowerSystems,2008,32(22):-103107.[17]林京屈梁生.基于连续小波变换的信号检测技术与故障诊断分析[J].机械工程学,200036-报12:95100.,,()1Sun了anY.Sinularityanalsisusincontiimouswavelettransformforbearinfault[巧Q,ggygg-化anosis口?Me浊anicalSstemsandSinalProcessin2002166:10251041.g]ygg,,()62 华北电力大学硕±学位论文[1却杨宇,于德化程军圣基于EMD的奇异值分解技术在滚动轴承故障诊断中的应-用饥.052425.振动与冲击,20:121,()[20王太勇,何慧龙王国锋.基于经验模态分解和最小二乘法支持向量机的滚动轴承],-%故障诊断机.机械工程学报.2007434:88.,()PULiHZhenHYanSP.TheFaultDianosisofBearinBasedonEMDandTeaer,gQ,ggggEnerOerator-J.VibrationandShock20082710:1518.gyp[],,()[22]雷亚国,何正嘉,紫艳阳.基于混合智能新模型的故障诊断[J].机械工程学报,2008,7-44:112117.()[23]张琳,黄敏.基于EMD与切片双谱的轴承故障诊断方法[J].北京航空航天大学学2010363287-290报.:,,()[24]WuZH,HuangNE.Ensembleempiricalmodedecomposition:anoiseassisteddataana-lsismethodJ.dvancesinAdaveDatanalsis200911141.y[]AptiAy,,():2李友荣肖涵等..基于加权相空间重构降噪及样本摘的齿轮故障分类J[巧呂勇,,,[]25-振动工程学报:4624.,2009,2()662..[巧杨世锡,汪慰军柯尔莫哥洛夫滴及其在故障诊断中的应用阴拆械科学与技术,2000-;191:68.,()27李兵张培林米巧山GA-SVM的轴.基于EMD与.机械强[],等承故障诊断[J]度,,,2010-:58:3233362.()口巧杨国安.机械设备故障诊断实用技术M.北京:中国石化出版社2007.],[29苏文胜.滚动轴承振动信号处理及将征提取方法研究D工大学.大连:大连理[][],2010.口0黄建鸿.基于小波包分析的滚动轴承故障智能诊断D].南昌大学2005.][,1郭奇刘h瑜史立波二ner-e口.基乎Vill]等次EEMD的Wig分布旋转机械故障,,,20-信号分析及试验研究[J].振动与冲击,123113129133.,():口3李揖郑海起唐力佑基于EMD和功率谱的齿轮故障诊断研巧J.振动与冲击]20061-25:133135.,()3张超军郭迅.基于EMD能量俯和支持向量机的齿轮故障诊断方法阴.振口],陈建,20-动与冲击,10290:216220.,(1)[34]鞠萍华,秦树人,秦毅,等.多分辨EMD方法与频域平均在齿轮早期故障诊断中的研究J-.振动与冲击20092859101.[]:7,,()口5]时世晨,单佩韦.基于EEMD的信号处理方法分析和实现饥.现代电子技术,20-11341;8890.,()3.MD.[巧曹冲锋,杨世锡,杨将新大型旋转机械非平稳振动信号的EE降噪方法机振动与冲击2009-38289:巧.,,()[37]田海雷,李洪儒,许谋华.基于EEMD和平滑能量算子解调的轴向柱塞累故障特征提取机-.海军工程大学学报20131:4347.,()口巧陈%李野陈满.EEMD分解在电力系统故障信号检测中的应用化计算机仿真,,63 华北电力大学硕±学位论文:%3福20103.()3-.基于改进证化ertHuan换的机械故障诊断町机械工程学化2011[刊雷亚国g变,475-(:7177.)40林近山.基于EEMD和脚bert变换的齿轮箱故障诊断町化械传动2010345:[],,()62-64.["]程军圣,于德介,杨宇.基于SVM巧EMD包络谱的滚动轴承故障诊断方法阴,2006259-系统工程理论与实践:131136.,,(),4^JoachimsT.MakinlarescaleSVMleaminracticalJ.1999.[gggp[]43LeslieCSElkinENobleWS.Thesectrumkernel:AstrinkernelforSVMrotein[],,pgpdasPa2-sificatio打C//ci行csmoshimonbiocomuti打.2007:566575.[]yppg,4hen-YZhouXSHuanS.OneclassSVMfbrlearnininimaeretrievalC//Imae[句C,,g了gg[]gProcess-in2001.roceedinnernationalferenceon.EE20011:3437.Ps,2001ItConIEg,g,,64 华北电力大学硕±学位论文攻读硕±学位期间发表的论文及其它成果[1]刘觉晓,杜冬梅.EMD与EEMD在滚动轴承故障诊断中的比较.华北电力大学第十二届研究生学术交流年会,2014化 华北电力大学硕止学位论文致谢,岁月如梭,白鞠过隙,转眼间,两年半的硕±生涯即将结束,期间有苦也有乐更!多的是收获!在此向所有关也帮助过我的老师、同学、朋友、家人致W由衷的感谢感谢我的导师,华北电力大学能动学院杜冬梅教授,在论文工作期间,杜冬梅老师严谨的治学态度、渊博的专业科学知识,给了我深深的感动,从她那里除了学到了很多相关专业知识外,我也从她身上学到了很多做人做事的优良作风。为了我的毕业论文的!完成,,杜老师倾注了大量的也血在此谨向我的导师致W深深的谢意感谢何青教授,,在我攻读硕±的两年半的时间里给予我学习上,W及生活上的巨大帮助。感谢何青教授在我科研课题上的帮助,给我指出了研究方向W及在遇到困难时为我指点迷津!,在此谨向我的指导老师致W深深的感谢感谢华北电力大学电站设备状态监测与故障诊断实验室的同学,褚东亮博±、张陳傅±、张昭、李展,W及范世岩、刘猜、徐威、赵龙等人,在硕±期间给予我无私的帮助与鼓励!,带给了我们这个实验室很多难忘的回忆,再次感谢感谢华北电力大学6号宿舍楼113宿舍的室友,吴松畔、高征、曹雷,感谢你们陪!伴我度过两年时间,在我迷茫的时候给我鼓励,在我孤单的时候给我关怀,谢谢一最后,特别感谢养育我的父母,是你们后援与支持,我才路上后顾无忧,只希望在不远的将来能回报你们的培养。66

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
关闭