基于時變?yōu)V波經(jīng)驗?zāi)B(tài)分解的轉(zhuǎn)子系統(tǒng)故障診斷

發(fā)布時間：2020-09-15 18:27

　　汽輪機、壓縮機、燃?xì)廨啓C等旋轉(zhuǎn)機械是工業(yè)生產(chǎn)中常見的設(shè)備類型,在工業(yè)領(lǐng)域具有重要的應(yīng)用。研究有效的方法監(jiān)測設(shè)備的健康狀況,診斷故障的類型、位置、嚴(yán)重程度,對于保證機械設(shè)備的正常運行、避免經(jīng)濟(jì)損失和人員傷亡具有重要意義。振動診斷方法是轉(zhuǎn)子系統(tǒng)故障診斷領(lǐng)域普遍采用的方法,主要包括故障機理、故障特征提取和故障模式識別三個方面。因此,本文在時變?yōu)V波經(jīng)驗?zāi)B(tài)分解的基礎(chǔ)上,分別從以上三個方面展開轉(zhuǎn)子系統(tǒng)故障診斷的研究。首先,總結(jié)了轉(zhuǎn)子系統(tǒng)常見的碰磨、裂紋、油膜渦動、不平衡故障的故障特征;建立了典型故障狀態(tài)下的的轉(zhuǎn)子系統(tǒng)動力學(xué)模型,采用計算機數(shù)值分析方法求解各個故障的運動微分方程,得到各個故障的仿真信號,從而為開展轉(zhuǎn)子故障診斷研究奠定理論基礎(chǔ)。為有效提取轉(zhuǎn)子故障信號的特征,將時變?yōu)V波經(jīng)驗?zāi)B(tài)分解(Time Varying Filtering Based Empirical Mode Decomposition,TVFEMD)應(yīng)用于轉(zhuǎn)子系統(tǒng)的故障診斷。首先,針對該方法存在的參數(shù)選取盲目問題,研究了采用粒子群算法進(jìn)行參數(shù)尋優(yōu)的參數(shù)優(yōu)化TVFEMD方法。然后,采用參數(shù)優(yōu)化TVFEMD與Hilbert變換結(jié)合,作為一種時頻分析方法,從轉(zhuǎn)子故障信號中提取故障特征。最后,經(jīng)對轉(zhuǎn)子故障仿真信號和實驗采集的故障信號分析表明,該時頻分析方法能夠從故障信號中有效提取故障特征,且與現(xiàn)有方法相比具有較強的優(yōu)越性。在參數(shù)優(yōu)化TVFEMD與Hilbert變換相結(jié)合的時頻分析方法基礎(chǔ)上,根據(jù)不同故障信號時頻分布特征的區(qū)別,研究了一種基于信息熵和粒子群優(yōu)化支持向量機實現(xiàn)故障類型識別的方法。首先,對轉(zhuǎn)子故障信號進(jìn)行參數(shù)優(yōu)化時變?yōu)V波經(jīng)驗?zāi)B(tài)分解和Hilbert變換,獲取信號的Hilbert時頻圖。然后,計算時頻圖的部分頻帶能量熵。最后,將熵值作為特征向量輸入支持向量機實現(xiàn)對轉(zhuǎn)子狀態(tài)的識別。經(jīng)對仿真和實測轉(zhuǎn)子故障信號進(jìn)行診斷,結(jié)果表明,該方法能夠準(zhǔn)確判別正常、不平衡、碰磨、油膜渦動等轉(zhuǎn)子系統(tǒng)典型狀態(tài),且不受轉(zhuǎn)子轉(zhuǎn)速影響;與EMD結(jié)合時頻熵的方法相比,具有更高的診斷精度。
【學(xué)位單位】：華北電力大學(xué)
【學(xué)位級別】：碩士
【學(xué)位年份】：2019
【中圖分類】：TH133
【部分圖文】：

碰磨,仿真信號

kck1k2czmxPtωo’oa）碰磨轉(zhuǎn)子系統(tǒng)模型 b) 碰磨力模型圖 2-1 轉(zhuǎn)子碰磨故障模型2.1.3 數(shù)值仿真給定系統(tǒng)的參數(shù)為：m 2.8kg， 0.4， 3，q e 1，n 1， 0.4 1.65， f 0.1， 0.3， 2。將參數(shù)帶入公式（2-4）中，使用標(biāo)準(zhǔn)定長的四階 Runge-Kutta 法進(jìn)行數(shù)值求解，以得到轉(zhuǎn)子的振動響應(yīng)。誤差控制量10-3。圖 2-2 為 y 方向的碰磨仿真結(jié)果。從中可以看出，該信號存在豐富的頻成分，由 50 Hz 的基頻及 1/2X、3/2X 等分?jǐn)?shù)諧波及 2X、3X、4X 等高倍諧波成；軸心軌跡呈不規(guī)則的花瓣形；龐加萊圖呈現(xiàn)不規(guī)則曲線。以上特征均符合2-1 所示的轉(zhuǎn)子系統(tǒng)碰磨的故障特征，從而證明使用該仿真信號進(jìn)行碰磨故障究的有效性和準(zhǔn)確性。

轉(zhuǎn)子裂紋,故障仿真

可知系統(tǒng)的一階臨界轉(zhuǎn)速為 radsMk882.5/20 。不同裂紋深度下的剛變量設(shè)置參見文獻(xiàn)[51]。采用四階 Runge-Kutta 法求解公式（2-12）所示的系統(tǒng)的運動微分方程，求解周期 T 2 ，步長為 0. 01 T。初值設(shè) 0.1000000011221122x yxyx y x y 。為消除瞬應(yīng)，舍棄前 500 個周期，計算 100 個周期的動力響應(yīng)。圖 2-5、2-6、2-7、2-8 分別為裂紋深度 A 0. 8 R時不同轉(zhuǎn)速下系統(tǒng)的波形、頻譜、軸心軌跡和龐加萊截面圖。從圖 2-5、2-6 中可以看出，在 1 /3、等分?jǐn)?shù)次臨界轉(zhuǎn)速下，系統(tǒng)響應(yīng)由基頻及其倍頻組成，具有明顯的分?jǐn)?shù)次共象，且隨著轉(zhuǎn)速提高，基頻及倍頻成分增多，幅值加大；軸心軌跡近似呈橢圓且隨轉(zhuǎn)速提高不斷接近圓形；龐加萊映射圖為一個單獨的點。從圖 2-7、2可以看出，隨著轉(zhuǎn)速增加，波形出現(xiàn)調(diào)制現(xiàn)象，不規(guī)則程度不斷加重；頻譜由1 /2等分?jǐn)?shù)次諧波、基頻與倍頻組成，表明系統(tǒng)發(fā)生了油膜振蕩；軸心軌跡環(huán)套小環(huán)的曲線逐步成為圓環(huán)狀；龐加萊映射圖由幾個獨立的點逐步變?yōu)榉馇€。

轉(zhuǎn)子裂紋,故障仿真,龐加萊映射,軸心軌跡