軸承振動信號瞬時頻率分析診斷方法
本發(fā)明涉及一種對軸承部件的分析診斷方法特別是涉及一種針對軸承振動信號瞬時頻率分析診斷方法。
背景技術(shù)：
軸承是工業(yè)生產(chǎn)設(shè)備中的關(guān)鍵部件由于長期連續(xù)工作在高載荷、高轉(zhuǎn)速下發(fā)生故障的比例極高。同時各類旋轉(zhuǎn)機械設(shè)備不斷向復(fù)雜、高速、高效的方向發(fā)展對于轉(zhuǎn)速的要求日益提高轉(zhuǎn)速波動對軸承故障診斷所帶來的影響已不容忽視。然而從另一個角度看，在轉(zhuǎn)速發(fā)生改變的過程中一些原本微弱的故障信息反而會因此而增強對其診斷也帶來了有利因素。
變速軸承的振動信號在時域中具有時變性、非平穩(wěn)性特點所以采用階比分析使其轉(zhuǎn)化為角域上的平穩(wěn)信號可以避免在信號采集預(yù)處理時丟失故障特征的問題。變速工況下對軸承振動信號進行階比分析的關(guān)鍵是獲取軸承的轉(zhuǎn)速信息。獲取轉(zhuǎn)速信息通常有兩種方法：安裝硬件轉(zhuǎn)速計方法和采用無轉(zhuǎn)速計分析算法。工業(yè)現(xiàn)場眾多設(shè)備的軸承等部件都包含在機殼內(nèi)部這使得轉(zhuǎn)速計的安裝十分困難既費時又浪費資源。所以本發(fā)明采用無轉(zhuǎn)速分析算法來獲取轉(zhuǎn)速該方法無需硬件配合直接從振動信號中獲取速信息可有效避免上述問題。從軸承振動信號中提取含有轉(zhuǎn)速信息的諧波分量對其瞬時頻率進行估計然后根據(jù)瞬時頻率和轉(zhuǎn)速的轉(zhuǎn)換關(guān)系獲得轉(zhuǎn)速信息。變速工況下軸承的振動信號具有非線性、時變性、調(diào)制性和多分量特性等特點各諧波成份十分相近甚至有重疊直接采用短時傅里葉變換、經(jīng)驗?zāi)B(tài)分解等時頻分析方法難以對其進行有效分解產(chǎn)生譜圖模糊瞬時頻率估計不精確等問題這些問題均需要進一步解決。
技術(shù)實現(xiàn)要素：
本發(fā)明的目的在于提供一種針對軸承振動信號瞬時頻率分析診斷方法該方法提出一種基于廣義解調(diào)時頻分析的信號瞬時頻率估計方法直接從軸承的振動信號中提取含有轉(zhuǎn)速信息的諧波分量通過估計其瞬時頻率來獲取轉(zhuǎn)速信息對于在轉(zhuǎn)速波動工況下具有時變性、調(diào)制性及多分量特性特別是諧波成分相近或出現(xiàn)重疊的振動信號的分析十分有效。
本發(fā)明的目的是通過以下技術(shù)方案實現(xiàn)的：
一種針對軸承振動信號瞬時頻率分析診斷方法，油漆測厚儀所述方法為一種基于廣義解調(diào)時頻分析的信號瞬時頻率估計方法該方法包括如下步驟：
步驟1：對變速工況下軸承原始信號進行最大離散小波包變換獲取時頻譜圖然后采用最大能量法分離出時頻譜圖中能量最為顯著的諧波分量用于下一步分析；
步驟2：采用最小二乘法擬合步驟1中得到的諧波分量進而得到其相位函數(shù)；
步驟3：對振動信號進行Hilbert變換得到解析信號其中是的Hilbert變換；
步驟4：根據(jù)步驟2中求取的相位函數(shù)對進行廣義解調(diào)得到信號；
步驟5：為去除負頻率的影響以及便于逆廣義解調(diào)再次對進行Hilbert變換得到新的解析信號其中是的Hilbert變換；
步驟6：構(gòu)建帶通濾波器其中振動信號的時頻譜圖中起始時刻對應(yīng)的頻率作為濾波器的起始頻率與其相鄰的2個分量起始頻率間隔的1/2作為濾波器的頻帶；
步驟7：利用帶通濾波器對解析信號進行帶通濾波分解出平行于時間軸的基頻諧波分量信號；
步驟8：對解析信號進行逆廣義解調(diào)得到解析信號；
步驟9：根據(jù)公式估算解析信號的瞬時頻率值；對離散的瞬時頻率值進行多項式擬合后獲得瞬時頻率估計曲線。本發(fā)明的優(yōu)點與效果是：
本發(fā)明對于軸承在變速工況下的振動信號所表現(xiàn)出的時變性、調(diào)制特性等情況的分析十分有效可解決頻率相近或頻帶重疊情況下的多分量信號的譜圖模糊問題提高分離精準性。本發(fā)明為多分量調(diào)制信號的分析提供了一種途徑利用廣義解調(diào)時頻分析方法并結(jié)合帶通濾波器可有效分離頻率相近或頻帶重疊的諧波分量提高瞬時頻率的估計精度。
技術(shù)特征：
1.一種針對軸承振動信號瞬時頻率分析診斷方法其特征在于所述方法為一種基于廣義解調(diào)時頻分析的信號瞬時頻率估計方法該方法包括如下步驟：
步驟1：對變速工況下軸承原始信號進行最大離散小波包變換獲取時頻譜圖然后采用最大能量法分離出時頻譜圖中能量最為顯著的諧波分量
用于下一步分析；
步驟2：采用最小二乘法擬合步驟1中得到的諧波分量進而得到其相位函數(shù)
；
步驟3：對振動信號進行Hilbert變換得到解析信號
其中
是
的Hilbert變換；
步驟4：根據(jù)步驟2中求取的相位函數(shù)對
進行廣義解調(diào)得到信號
；
步驟5：為去除負頻率的影響以及便于逆廣義解調(diào)再次對進行Hilbert變換得到新的解析信號
其中
是
的Hilbert變換；
步驟6：構(gòu)建帶通濾波器其中振動信號
的時頻譜圖中起始時刻
對應(yīng)的頻率作為濾波器
的起始頻率
與其相鄰的2個分量起始頻率間隔的1/2作為濾波器
的頻帶
；
步驟7：利用帶通濾波器對解析信號
進行帶通濾波分解出平行于時間軸的基頻諧波分量信號
；
步驟8：對解析信號進行逆廣義解調(diào)得到解析信號
；
步驟9：根據(jù)公式估算解析信號
的瞬時頻率值；對離散的瞬時頻率值進行多項式擬合后獲得瞬時頻率估計曲線
。
一種針對軸承振動信號瞬時頻率分析診斷方法涉及一種對軸承部件的分析診斷方法。對軸承振動信號進行最大離散小波包變換采用最大能量法從其取時頻譜圖中分離出時頻譜圖中能量最為顯著的諧波分量并用最小二乘法擬合此諧波分量進而得到其相位函數(shù)；然后根據(jù)得到的相位函數(shù)對振動信號進行廣義解調(diào)分析并構(gòu)建帶通濾波進行濾波；最后，對帶通濾波后的諧波分量進行逆廣義解調(diào)估算逆廣義解調(diào)后信號的瞬時頻率并對其進行曲線擬合最終得到振動信號的瞬時頻率估計。本發(fā)明用于變速工況下對軸承振動信號的瞬時頻率進行估計以實現(xiàn)無轉(zhuǎn)速計階比跟蹤分析變速工況下軸承故障診斷、預(yù)測與健康管理。
技術(shù)研發(fā)人員：齊曉軒;原忠虎;劉英英;杜英魁;都麗
受保護的技術(shù)使用者：沈陽大學
技術(shù)研發(fā)日：2018.06.20
技術(shù)公布日：2018.11.30
相關(guān)推薦：
