基于振動信號的齒輪故障診斷方法研究畢業論文-資料下載頁

　　

【正文】 0 2 4 6 0time一s一signalvalue一一D3一一一圖 38 大偏心齒輪振動信號三層小波分解時域圖0 2 4 6 0time一s一signalvalue一一一A1一一一0 2 4 6 0time一s一signalvalue一一D1一一一0 2 4 6 0time一s一signalvalue一一一A2一一一0 2 4 6 0time一s一signalvalue一一D2一一一0 2 4 6 0time一s一signalvalue一一一A3一一一0 2 4 6 0time一s一signalvalue一一D3一一一圖 39 小偏心齒輪振動信號三層小波分解時域圖小波變換后輸出為信號的時域波形圖，時域波形圖可以看出四種齒輪狀態信號北京化工大學畢業設計(論文)27的明顯不同。從上圖中可以明顯看出無論是低頻還是高頻部分，信號的幅值都比前幾幅圖的幅值增大，對于齒輪局部損壞和斷齒的狀態，時域圖中幅值有明顯的增加。小波變換后，我們可以看到高頻信號的幅值也有所增加。所以通過小波分解后的時域圖對比可以粗略判別信號的正常與故障。第 節 頻域分析在齒輪故障診斷中的應用基于一維離散小波變換，使原始波形細化為三個低頻波和三個高頻波，時域波形是在時空坐標下的表示，表示的是時間和振幅的關系，但是要想從頻率角度對齒輪故障進行分析，就要從頻率和振幅角度分析，因此，在上一章節的基礎上在對信號進行FFT變換（快速傅里葉變換）。 MATLAB 中的 FFT 變換快速傅里葉變換（FFT）， x表示經離散小波變換后分離出來的低頻信號， y表示經離散小波變換后分離出來的高頻信號。程序如下： k=2^nextpow2(length(x))。 f=fs*(0:((k/2)1))/k。 x=x(ones(1,length(x)))39。*mean(x)。 y=y(ones(1,length(y)))39。*mean(y)。 X=fft(x,k)。 Y=fft(y,k)。 Pxx=X.*conj(X)/k。 Pyy=Y.*conj(Y)/k。 kk=1024*20。 f1=f(1:kk)。程序運行后，得到頻譜圖，下面是經小波變換及FFT變換后四種不同狀態的齒輪的頻譜圖。圖3 11為正常狀態齒輪的頻譜圖，圖312為局部損壞狀態齒輪的頻譜圖，圖3 13大偏心狀態齒輪的頻譜圖，圖314小偏心狀態齒輪的頻譜圖。從這四幅圖中可以看出，高頻信號的頻域圖的特征不是很明顯，所以本文以介紹低頻信號的頻域圖北京化工大學畢業設計(論文)28為主。 頻譜分析通過齒輪故障診斷試驗提供的參數我們計算齒輪的嚙合頻率及每根軸的轉頻。本實驗的齒輪箱為一級齒輪傳動，所以只有兩根軸。見圖310。 齒輪箱簡圖圖310 齒輪箱簡圖參數計算： 的轉動頻率: 軸 Ⅰ 1056zrpmfHZ? 的轉動頻率:軸 Ⅱ ?齒輪的嚙合頻率： 125mzzffZ??0 500 1000 1500 202202040frequence(Hz)signalvalue一一A1一一一0 500 1000 1500 2022024x 1013frequence(Hz)signalvalue一一D4一一一0 500 1000 1500 202202040frequence(Hz)signalvalue一一A2一一一0 500 1000 1500 202201x 1010frequence(Hz)signalvalue一一D5一一一0 500 1000 1500 202202040frequence(Hz)signalvalue一一A3一一一0 500 1000 1500 2022024x 107frequence(Hz)signalvalue一一D6一一一圖311 正常狀態齒輪的頻譜圖??1=55??=300????????1=75軸 Ⅰ 軸 Ⅱ北京化工大學畢業設計(論文)290 500 1000 1500 202202040frequence(Hz)signalvalue一一A1一一一0 500 1000 1500 2022024x 1013frequence(Hz)signalvalue一一D4一一一0 500 1000 1500 202202040frequence(Hz)signalvalue一一A2一一一0 500 1000 1500 202201x 1010frequence(Hz)signalvalue一一D5一一一0 500 1000 1500 202202040frequence(Hz)signalvalue一一A3一一一0 500 1000 1500 202205x 107frequence(Hz)signalvalue一一D6一一一圖312 局部損壞狀態齒輪的頻譜圖0 500 1000 1500 2022050100frequence(Hz)signalvalue一一 A1一一一0 500 1000 1500 202201x 1012frequence(Hz)signalvalue一一 D4一一一0 500 1000 1500 2022050100frequence(Hz)signalvalue一一 A2一一一0 500 1000 1500 2022012x 1010frequence(Hz)signalvalue一一 D5一一一0 500 1000 1500 2022050100frequence(Hz)signalvalue一一 A3一一一0 500 1000 1500 2022024x 107frequence(Hz)signalvalue一一 D6一一一圖313 大偏心狀態齒輪的頻譜圖北京化工大學畢業設計(論文)300 500 1000 1500 202201020frequence(Hz)signalvalue一一 A1一一一0 500 1000 1500 2022024x 1013frequence(Hz)signalvalue一一 D4一一一0 500 1000 1500 202201020frequence(Hz)signalvalue一一 A2一一一0 500 1000 1500 2022012x 1010frequence(Hz)signalvalue一一 D5一一一0 500 1000 1500 202201020frequence(Hz)signalvalue一一 A3一一一0 500 1000 1500 202201x 107frequence(Hz)signalvalue一一 D6一一一圖314 小偏心狀態齒輪的頻譜圖第 節 本章結論本文介紹時頻分析中的小波變換理論，并基于一維離散小波變換的基礎上，采取三層小波分解，得到了四種信號的不同時域圖，并對信號進行快速傅里葉變換，得到頻譜圖，并通過頻譜圖分析，進行故障識別及診斷。由于齒輪振動信號非常復雜，內部包含著很多元素，并且誤差較大，所以通過頻譜圖的直接診斷是有一定難度的。本文也只能通過幅值大小來粗略的辨別出故障信號的不同特征。但是可以證明小波變換在齒輪故障診斷的應用中是可行的。北京化工大學畢業設計(論文)31第 4 章 基于 EMD 的齒輪故障診斷第 節 基于 EMD（經驗模式分解）的振動信號特征提取 的研究背景齒輪摩擦是齒輪箱系統中常見的非線性振動故障。當齒輪存在局部摩擦時, 摩擦信號是在動靜摩擦的短時間內產生的，具有頻帶寬的特點。在采集振動信號時混有寬頻帶隨機噪聲的噪聲信號, 其頻帶相互疊加。在振動信號中還含有能量較大的、與轉速有關的背景信號, 所以利用傳統的濾波方法很難僅保留含有故障信息的摩擦信號而將背景信號和噪聲信號濾掉。由于用小波分析方法進行故障信號特征提取存在缺點， 等人針對瞬時頻率的概念進行深入研究后,創造性地提出了本征模式函數(Intrinsic Mode Function,IMF)的概念以及將任意信號分解為本征模式函數組成的新方法——基于經驗模式分解(Empirical Mode Deposition,EMD)方法，從而賦予了瞬時頻率合理的定義、物理意義和求法，初步建立了以瞬時頻率為表征信號交變的基本量，以本征模式分量為時域基本信號的新的時頻分析方法體系。EMD分解算法的基本思想是：對一給定信號，先獲得信號極值點，通過差值獲得信號包絡，獲得均值，與均值的差得到分解的一層信號；如此重復，獲得分解結果：北京化工大學畢業設計(論文)32，即 個IMF和一個殘值 。由于EMD方法具有自適應特性，適宜于非線性、非平rlr穩信號的分解，該方法應用于齒輪的故障診斷分析中。結果表明，該方法能夠突出齒輪振動信號的故障特征，從而提高齒輪故障診斷的準確性 [19]。為此，本文將EMD方法應用到齒輪的故障診斷中，利用EMD方法對齒輪信號進行分解，得到若干個IMF，然后對包含故障信息的IMF作功率譜分析，有效提取故障信號特征信息。 經驗模式分解(EMD)方法原理經驗模式分解方法的主要目的是通過對非線性、非平穩信號 的分解，獲得一()xt系列表征信號特征時間尺度的本征模式函數 和一個殘余分量 ，()0,?()nrt各個IMF是單分量的幅值或者頻率調制信號。每個IMF需要滿足以下幾個條件信號：(1)整個信號中零點數與極點數相等或者相差1；(2)信號上任意一點，由局部極大值點確定的包絡線和由局部極小值點確定的包絡線的均值均為零，即信號關于時間軸局部對稱 [20]。分解過程如下：（1（ 設信號為 ，取其上下包絡局部均值組成的系列為 ，令它們的差值為()xt 1()nt()ht (41)11()()htxnt??對于非線性、非平穩數據來說，一般一次處理不會形成基本模式分量，一些非對稱波依然存在。因此要把 看作待處理的數據繼續重復上述操作。取 的均1()t 1()ht值系列為 ，得到它們之前的差值為2()nt 2()ht (42)12()nt??重復 k 次操作，就可以得到： (43)(1)()kkkhtt?當 滿足基本模式分量的條件后，就會獲得了第一個基本模式分量 ()xt 1()ct (44)1kct?（2（ 將基本模式分量從信號中分離出來，可以得到剩余系列 1()rt (45)11()()rtxct?（3（ 把 作為新的數據按上述操作進行處理，依次類推，可得到：1()rt 212()()tt?北京化工大學畢業設計(論文)33323()rtct??…… (46)1()()nntt?這個處理過程在滿足預先設定的停止標準后就可停止。（4（ 將式(45)和式(46)相加后得到 (47)1()()ninxtctr???這樣就把一個信號 分解為 個基本模式函數 和一個殘余分()xt (0,123.)i n?量 [19,20]，其中分解出的 個分量 分別包含信號從高頻到低頻的不同頻率段成分，()nrt ni而剩余分量 是原始信號的中心趨勢。n(46)由于每個固有模態函數 為平穩信號，由式(47)可以定義每個固有模態分nr量的自相關函數為 (48)()[()]ciiREk?? 基于 EMD 的振動信號特征提取分析當齒輪發生局部摩擦故障時，設備振動信號中將出現周期性的沖擊，這些周期性的沖擊信號由于噪聲等干擾，所以很難通過頻譜分析等傳統方法加以識別，因此本文提出了一種基于 EMD 的沖擊信號提取方法。依照 IMF 的定義，周期性的沖擊信號無法以獨立的 IMF 形式被分解出來，但它會以類似沖擊響應信號的形式出現在高頻段的 IMF 中。假設如下仿真信號： (49)()sin2*50)()xttpnt???該仿真信號是由一個