質(zhì)量有缺陷的滾動軸承振動響應(yīng)的循環(huán)方式

1、外文資料譯文外文資料譯文質(zhì)量有缺陷的滾動軸承振動響應(yīng)的循環(huán)方式摘要由于軸承元件間不斷的滑動，在一個旋轉(zhuǎn)軸上的有缺陷的滾動軸承的震動信號是一個不穩(wěn)定狀態(tài)。正是由于這一原因，基于傅立葉分析的傳統(tǒng)的鑒別方式，并不能清晰的發(fā)現(xiàn)缺點的特征和調(diào)整現(xiàn)象，因為包含在第二順序統(tǒng)計方法中的信息被減少了。本文的主要目的是結(jié)合循環(huán)統(tǒng)計和高階統(tǒng)計分析方法，并且考慮到軸承狀態(tài)的鑒定，它能夠帶來更好的結(jié)果。本文證明了循環(huán)分析能引起一個離散譜結(jié)構(gòu)，能夠?qū)S承缺陷的種類和相應(yīng)調(diào)整結(jié)構(gòu)有一個清晰的察覺，從而克服其他現(xiàn)存方法的相應(yīng)問題。符號表x(n) 離散時間系列X(k) 離散傅立葉變換Sxx(k) 一個離散時間系列的能量譜E 統(tǒng)

2、計期望控制k,l 離散頻率變量B(k,l) 光譜fs 抽樣頻率Rx(t,) 自相關(guān)函數(shù) 循環(huán)頻率m(t) 一階時刻T 時間周期 傅立葉級數(shù)定義的自相關(guān)函數(shù)Sx(,f) 譜相關(guān)密度函數(shù)x(t) 時間信號X(f) 一個信號的傅立葉變換 循環(huán)譜 脈沖間隔 沖擊的固定振幅 滾動元件軸承在徑向載荷下的分布函數(shù) 最大載荷密度 載荷分布因素 最大載荷角 滯后時間 沖擊的隨機振幅 品質(zhì)因素 共振頻率 釋放時間 隨意振動頻率B 振幅 軸轉(zhuǎn)速 外滾道球轉(zhuǎn)頻率 內(nèi)滾道球轉(zhuǎn)頻率 中心頻率1.引言滾動軸承的振動狀態(tài)監(jiān)測一直都是一個廣泛研究的課題，因為它是一個最終要的和常用的旋轉(zhuǎn)機械的部件。有缺陷的滾動元件軸承的振動反應(yīng)

3、主要受調(diào)制效應(yīng)支配，并常由缺陷的類型、軸的轉(zhuǎn)速和軸承的幾何形狀來定義。這個概念已經(jīng)被許多基于振動信號譜分析的傳統(tǒng)狀態(tài)監(jiān)測方法所替代。這些方法假定信號統(tǒng)計特征是固定的，而它們的目標(biāo)是使波峰與潛在缺陷的頻率相聯(lián)系。不過由于軸承滾動元件間不斷的滑動，有缺陷的軸承振動反應(yīng)并不是與轉(zhuǎn)速嚴(yán)格鎖定的。因此，機械的振動反應(yīng)與損壞的滾動元件之間不僅表明為周期現(xiàn)象，而且存在一個明顯的非線形、非穩(wěn)定現(xiàn)象。因為功率譜分析存在丟棄所有階段資料的嚴(yán)重缺陷，所以它不能識別兩個階段的系列信號。結(jié)果。HOS（高階統(tǒng)計方法）在有缺陷的軸承狀態(tài)檢測上的應(yīng)用取得了局部成功。同時由于在反應(yīng)上的機械旋轉(zhuǎn)和調(diào)節(jié)作用，較普遍的周期假設(shè)，已廣

4、泛應(yīng)用于機械狀態(tài)檢測，不過在損壞滾動元件軸承上也有部分成功的案例。高階統(tǒng)計分析與周期分析在機械狀態(tài)檢測上的結(jié)合一直到最近才被考慮，隨后應(yīng)用于齒輪箱的故障檢測。這篇文章的目的是要證明結(jié)合譜，三階統(tǒng)計強大的技術(shù)，能夠識別非線性耦合的頻率，同周期分析一樣能夠提高現(xiàn)有分析技術(shù)的診斷能力，從而導(dǎo)致更明確的診斷結(jié)果。本文第2節(jié)提出了高階統(tǒng)計方法基本概念的簡要概述，周期和循環(huán)譜分析。第三節(jié)提出了一個損壞軸承的動態(tài)響應(yīng)的隨機模型。它表明，即使是像滑動參數(shù)一樣少量的隨機性，頻譜分析卻無法提供有意義的結(jié)果。三節(jié)的模型進一步應(yīng)用的四節(jié)來產(chǎn)生代表軸承內(nèi)外圈缺陷動態(tài)反應(yīng)的信號。這些信號經(jīng)過使用雙頻，周期和循環(huán)譜分析方法

5、的進一步分析，表明了循環(huán)譜分析極大的克服了其他兩種方法的問題，清楚的展示了缺陷的性質(zhì)和明顯的六角形頻率方式。最后，第五節(jié)在工業(yè)設(shè)備安裝測量上驗證結(jié)果。2.高階統(tǒng)計和周期循環(huán)的基本概念在本節(jié)概念的全部理論方法治療能夠發(fā)現(xiàn)在5.11為高階統(tǒng)計方法，在6.7為周期分析和在3為循環(huán)譜分析。2.1高階統(tǒng)計方法高階統(tǒng)計方法延伸到高階的常規(guī)二階功率譜方法。如果信號是高斯概率密度函數(shù)，二階措施就能提供令人滿意的結(jié)果。然而，許多現(xiàn)實生活的信號非高斯。在信號處理中最普遍使用的工具是功率譜。一個離散的時間序列X（n）的功率譜被定義如下，利用離散傅立葉交換（DFT）X（K）的信號：Sxx (k) = EX(k)X*（

6、k） (1)其中E 表示統(tǒng)計期望操作，K是離散變頻，功率譜是在整個頻域中的信號的能力分布。因為這是一個現(xiàn)實的數(shù)量。但并包含任何階段信息。一個可以考慮的頻譜是第二時刻的信號。高階統(tǒng)計方法由延伸的二階方法的定義形成。給高階統(tǒng)計定義的另一個途徑是利用多種譜。因此，第三階譜被稱為譜。像功率譜分解信號的功率一樣，譜在頻域內(nèi)使信號分解成三個時刻。在非線性系統(tǒng)中，頻率元件成對的在一起，由于它們依賴于固定的假設(shè)，功率譜措施不能在其中有效的實施。因此，功率譜措施的有限性導(dǎo)致了高階性能信號評價的需求。像功率譜一樣，譜能夠被DFT信號定義為：B(k, l) = EX(k)X(l)(k + 1). (2)我們看到譜是

7、個復(fù)雜值，由兩個以上的光譜頻率值f1和f2定義，分別由離散頻率指數(shù)k、l和包含階段信息組成。由DFT特性支配的功率譜，不包含在Niquist頻率f2/2的信息中。因此，沒必要為所有頻率對（k，l）計算譜B（K，L）的值，而在平面（K，L）中存在著幾個對稱關(guān)系。非冗余地方被稱為主域并且定義為：k， l: 0 k fs /2,f k,2k + 1 fs . (3)該主域構(gòu)成一個三角形，它又被分為兩個新的三角形，稱為內(nèi)外三角形。由于譜涉及到一個信息的偏斜，所以它被用做追蹤一個對稱的非線性。當(dāng)信號非偏斜時，該譜為0。在頻率譜（K，L）上譜的振幅測定在頻率K、1和K+1間的光譜成分中偶合的數(shù)量。這種偶合

8、意味著那兒存在一個二次非線性信號。2.2周期分析大多數(shù)確定的統(tǒng)計信號處理方法考慮到了信號的統(tǒng)計特性，如固定。相反，周期分析認(rèn)為周期性的時間變化是靜止的6.7。如果一個信號的統(tǒng)計時刻直至幾階具有周期性的時間依賴性，它就被稱為幾階循環(huán)。這樣的一個周期性的基本頻率被稱為信號的循環(huán)頻率。信號的第一階循環(huán)與一階時刻有關(guān)：m(t) = Ex(t)= m(t + T). (4)第一階循環(huán)實際應(yīng)用于通過實施轉(zhuǎn)動輪換平均循環(huán)技術(shù)的振動分析中。二階循環(huán)表示了一個周期性的時間變化自相關(guān)函數(shù)，一個信號x(t)的自相關(guān)函數(shù)可以由公式定義為：Rx (t, ) = Ex(t + /2)x(t /2). (5)因為上面的函數(shù)

9、是周期性的，所以它可以擴大到級數(shù)：Rx (t, t ) =R ()ej2(t/2), (6)其中R是確定該循環(huán)自相關(guān)函數(shù)的傅立葉系數(shù)，a=1/t是周期頻率。由于循環(huán)元件的頻率a，循環(huán)自相關(guān)函數(shù)顯示出信號能量的數(shù)量。同樣，一個信號的功率譜密度函數(shù)能夠由它的固定自相關(guān)函數(shù)的傅立葉變換計算出來，譜相關(guān)密度函數(shù)能夠由循環(huán)自相關(guān)函數(shù)的傅立葉變換求出，考慮到時間變化t。譜相關(guān)密度函數(shù)也表示為：Sx (, f) = EX(f /2)X *(f + /2), (7)其中，x(f) 是x(t)的傅立葉變換，根據(jù)公式7，SCDF計算了一個頻率f和被一個頻率位移分解成/2的頻率元件的相關(guān)性。此外，一個信號的周期函數(shù)

10、密度（DCS）可定義為： 周期密度函數(shù)類似于信號外部的功率譜密度的作用和簡單性。2.3 循環(huán)譜這篇文章的目的是集中在損壞軸承信號的某一雙頻特性上，這些缺陷能揭示缺陷的本質(zhì)。循環(huán)譜分析結(jié)合了損壞軸承振動信號的周期和高階屬性。周期分析能識別幾個調(diào)制機械信號的相關(guān)性，通過支配這些信號統(tǒng)計時刻周期的循環(huán)頻率來表示，另一方面，譜能檢測相耦合的組頻率元件。換句話說，算法是建立在小改動譜的計算基礎(chǔ)上，為的是提出循環(huán)頻率的概念。因而，離散信號X(n)的循環(huán)譜可定義為：B (k, l) = E X (k) X (l) X (k + l ). (9)公式9中，循環(huán)譜的介紹需要三個不同的頻率軸：A）.像公式（2）中

11、由指數(shù)K和L分別表示一樣，譜軸線f1和f2由譜分析表示。B）.循環(huán)頻率軸，因為三個不同軸結(jié)果的介紹在實際中幾乎不可能，世界循環(huán)譜分析涉及到公式9的計算僅是為循環(huán)頻率a找一個譜值。因此，等高線圖幾乎產(chǎn)生類似于傳統(tǒng)的譜分析。但是，每一個對應(yīng)一種特定的循環(huán)頻率.滾動軸承振動反應(yīng)的隨機性大量模型已用來描述滾動軸承在不同類型缺陷下的動力行為根據(jù)傳統(tǒng)方式9，由固定軸承缺陷產(chǎn)生的重復(fù)沖擊可以由下面形式的狄克拉三角函數(shù)系列表示為： 其中，d是脈沖周期，被定義為滾動軸承特征缺陷頻率d0是脈沖的振幅,顯示缺陷的嚴(yán)重程度。q(t)是滾動軸承徑向負(fù)荷的分布函數(shù)，近似于著名的stribeck方程9。 但是，一個缺陷軸承

12、滾動元件動態(tài)反應(yīng)的重大變化是存在的，當(dāng)上述方程中的某些參數(shù)被假定為隨機的而不是固定的。例如，沖擊影響，排除軸承周向負(fù)荷分布，裝配動剛度的變化，滾動元件和圈體的波動，以及球的配合是否合適等。因此，這一系列的影響應(yīng)視為振幅的隨機調(diào)制。與此同時，滾動軸承經(jīng)歷了大量滑動，結(jié)果，脈動從不在同一周期同一位置發(fā)生。因此，有缺陷的滾動軸承脈動系列的一個更現(xiàn)實的模型可以如下表示： 其中，Ak是沖擊力振幅的隨機變量，具有一個均值d0和假定正常(高斯)的概率密度函數(shù)。k是由于存在滑動的兩個沖擊間的時間差的隨機變量，該滑動被假定為一個0和正常的概率密度函數(shù)。式12的脈沖列產(chǎn)生了脈沖激振力，在軸承和機構(gòu)間產(chǎn)生共振。簡單

13、假定激發(fā)機構(gòu)是一個多自由度線形系統(tǒng)，每個脈沖的結(jié)構(gòu)反映可為： 其中，i=1,M是激起模式，對每個模式i. mi是恢復(fù)時間，fni 是共振頻率，Qi是品質(zhì)因數(shù)。因此，產(chǎn)生于軸承生產(chǎn)缺陷的動態(tài)反應(yīng)可表示為： 其中，符號表示卷積，n(t)是存在的背景噪音。為了展出某些關(guān)鍵參數(shù)隨機性的效果，對振動光譜，如滑動就是一個典型的例子。公式15產(chǎn)生的模擬信號，對應(yīng)的是一個外圈有缺陷軸承的典型反應(yīng)。軸轉(zhuǎn)速是19.35HZ。特性軸承缺陷頻率（球傳頻率外圈）相當(dāng)于3.75倍的主軸轉(zhuǎn)速，導(dǎo)致估計的大約為72.47HZ。該系統(tǒng)的激發(fā)自然頻率假定為3975HZ，信號由16384個樣本組成，抽樣率為40HZ，品質(zhì)因數(shù)=12

14、。最后，假定滑動百分比從0變化到1.44%.圖1顯示了光譜的模擬信號。公式9中定義的模型的預(yù)期結(jié)果是不完善的，一套頻率為72.47HZ的頻率的諧波出現(xiàn)零滑動。在3975HZ的共振頻率附近的頻帶振幅是最大的。此頻率沒有在頻譜中出現(xiàn)，因為它不一定就是一個諧波的頻率。不過，圖1最大重要的特征是即使一個小的滑動，光譜元件特征缺陷就已不明顯，并且激發(fā)頻帶包括冗余信息。原因是二階特性（能量，差額）的所有信息受滑動現(xiàn)象的幅度和發(fā)生率的約束。因而，由軸承滑動產(chǎn)生的大量噪聲和隨機性導(dǎo)致SNR的減少和數(shù)據(jù)屬性的扭曲。由于這一事實，先進統(tǒng)計方法的使用是必要的。4.損壞滾動軸承反應(yīng)的特性為了檢測先進統(tǒng)計分析工具的實效

15、性，譜，二階循環(huán)分析和循環(huán)譜分析都在考慮中，并首先考慮了由隨機模型(15)產(chǎn)生的一系列模擬信號。考慮到的第一信號對應(yīng)于外圈缺陷軸承的典型反應(yīng)。外圈球轉(zhuǎn)動的頻率達到71.4HZ，同時激發(fā)結(jié)構(gòu)自然頻率fn1和fn2分別為2859HZ和3682HZ，選擇的品質(zhì)因數(shù)Q1和Q2為11HZ。模擬信號的采樣頻率為40HZ，它有8192樣本的長度，并且其軸轉(zhuǎn)速fshaft假定為18.3HZ。最終，認(rèn)為軸承的滾動元件在缺陷頻率周圍有從0到1.07%的偏差。圖1 軸承在外滾道有滑動偏移下的信號譜模擬振動反應(yīng)圖2表示了不滑的平穩(wěn)信號連同相應(yīng)的光譜。圖3表示一個變滑百分比高達1.07%的信號。與圖1的大部分光譜相似，

16、圖3b的光譜未能為檢測和表征缺陷的類型提供清晰而有價值的工具。觀察圖3a波形圖，與圖2a的相反，它指出了脈沖由于滑動效果而“溜走”。結(jié)果，能量的隨機部分大大高于固定部分。那么，周期分析完成。圖4用例子說明為模擬信號制作的SCDF的一部分等高線。SCDF的所有點在斜線上劃分。不同點形成了一系列被認(rèn)為與缺陷軸承的缺陷類型和調(diào)節(jié)機制有關(guān)的復(fù)雜菱形結(jié)構(gòu)。這些點的間距必須與BPFO頻率、它在水平方向的諧波以及它在垂直方向間距的2倍相等。根據(jù)這些說明，這些特征模式將提供自然頻率與故障頻率明顯的相關(guān)性。因而，如果這些條件成立，這種關(guān)聯(lián)就驗證了存在于信號中的具體的調(diào)節(jié)效應(yīng)。(a)(b)圖2 軸承在外滾道為零滑

17、動下的模擬反應(yīng) (a)圖3(b) 軸承在外滾道滑移變化達到1.07%下的模擬反應(yīng)圖4 在圖3a中信號的二階周期分析的輪廓描述圖5 在圖3a中信號譜的輪廓描述從檢測方面得來的有用信息集中于由于滾動軸承在不同類型缺陷下的動態(tài)行為所激發(fā)的地區(qū)。所以，我們的檢測研究集中在水平方向頻寬在2到4KHZ的頻率，因為它是系統(tǒng)自然頻率附近的激起頻率。在圖4中觀察到優(yōu)勢明顯的點形成了特別的菱形結(jié)構(gòu)，與我們期望的相反，縱向與橫向坐標(biāo)點和菱形對角線軸的尺寸不對應(yīng)于軸承的故障頻率或諧波。因而，圖4中出現(xiàn)的模型是不確定的，并且由軸承故障頻率BPFO產(chǎn)生的調(diào)節(jié)效應(yīng)不能輕易的被發(fā)現(xiàn)。從圖4的SCDF圖上可以清楚的看出在與故障

18、頻率相聯(lián)系的 光譜元素和調(diào)制元素間有內(nèi)在相關(guān)度，由于隨機成分信息的缺乏。這一缺乏是由于這一現(xiàn)實：滾動元件的滑移產(chǎn)生了一個隨機性，它擴大了信號中背景噪音的作用。因而，相同循環(huán)頻率的信號元素的連接是引用錯誤的，導(dǎo)致了所有相關(guān)循環(huán)的衰減。譜載是另一種可以特別使用的方法。該譜載作為一個與信號傾斜有關(guān)的非線形轉(zhuǎn)變的指示器出現(xiàn)，盡管不足二階特征，但有用的三階信息能夠在受非線形現(xiàn)象支配的信號中保留。在滑動下，信號的光譜分析可由圖5解釋。不幸的是，這個方法也不能產(chǎn)生有用的結(jié)果。圖5中的等高圖并不能提出一個可以呈現(xiàn)缺陷類型的有用模式。原因是三階統(tǒng)計（譜）的成功依賴于收到的數(shù)據(jù)中呈現(xiàn)的非線形與非高斯的程度。相反，

19、一個被認(rèn)為的高斯信號消失在第三累計階中。在觀察中的信號滑動產(chǎn)生的隨即性， 擴大了信號的高斯性，同時降低了階段關(guān)系和三階統(tǒng)計行為。因而，信號的非高斯性在增長程度和發(fā)生率的下滑上往往是微乎其微的。此外，事實是譜為整個頻域估計的結(jié)果是方法有效性的降低。最后一步的分析涉及到循環(huán)譜方法的應(yīng)用。根據(jù)2.3節(jié)的分析，實際的循環(huán)譜分析涉及到為循環(huán)頻率a的一系列具體值計算的公式9。根據(jù)公式15，最適合的一套循環(huán)頻率包括故障頻率BPFO和激發(fā)自然頻率。盡管BPFO頻率能通過類型軸承和轉(zhuǎn)速的測定，但激發(fā)自然頻率不能被實際估計，因為它們不能出現(xiàn)在譜中。因而，循環(huán)譜在實際應(yīng)用中，循環(huán)頻率a有兩個值可供選擇：A）、故障頻

20、率BPFO ，大致為71.4HZ。B)、.對應(yīng)最大峰譜（2891HZ）的頻率fc假定是接近共振頻率。 圖6 循環(huán)頻率等于BPFO頻率(71.4HZ)的信號循環(huán)譜的輪廓描述圖7 循環(huán)頻率等于中心頻率(2891HZ)的信號循環(huán)譜的輪廓描述圖6表示了循環(huán)頻率的首選模式?？煽闯?，顯示的模式由一些不同點的出現(xiàn)所支配。但是，這些點仍不能呈現(xiàn)一個規(guī)則結(jié)構(gòu)，因此，它們不能提供缺陷的指標(biāo)。原因是由于滑動的程度，軸承缺陷頻率BPFO并不是和軸頻率有固定聯(lián)系和嚴(yán)格鎖相的。因而，由于周期性，BPFO并不與相對有聯(lián)系，并且與其他譜元件也無相互聯(lián)系。圖7表示了循環(huán)頻率的第二選擇模式。等高點是受許多清晰存在的不同點支配，它

21、形成了一個大跨度和一個中心頻率為2891HZ的六邊形。圖8是這一地區(qū)進一步的放大顯示。在光譜頻率f1軸與f2軸的頻率2400HZ和3400HZ間的區(qū)域間，兩個特征六邊形清晰的出現(xiàn)在中心頻率對（f1,f2）附近。圖8 放大后的圖7中在2400到3400HZ的循環(huán)譜輪廓點應(yīng)該指出，六邊形的尺寸與軸承損壞頻率和3倍的軸承損壞頻率3×BPFO的和差密切相關(guān)。這個六邊形結(jié)構(gòu)形成了一個頻率模式，它清晰的顯示了軸承損壞的類型。這個明確的結(jié)構(gòu)是由于這樣的現(xiàn)實：循環(huán)譜方法能更好的為接近系統(tǒng)固有頻率的循環(huán)頻率提取周期性時間變化的三階統(tǒng)計。這個頻率與其它譜線的任何一個都相關(guān)。因而，在顯示信號的非線性的能力

22、上，該譜能更好地利用這一特定頻域，而不是在整個頻域，由于循環(huán)固有中心頻率fc的強度。因而，循環(huán)譜分析能夠再次提出一個明確的和離散的頻率模式，從中可以推斷出缺陷的類型和調(diào)控機制。5.工業(yè)應(yīng)用這個測量在Lavrion Attica的一個工業(yè)裝置的火電站的風(fēng)機上進行。測量期間，軸轉(zhuǎn)速在1236rpm(20.6HZ)附近。被測量的軸承是SKF的22244C3型，它的BPFO頻率為軸轉(zhuǎn)速的6.49倍，導(dǎo)致了理論估計的預(yù)期BPFO頻率在134HZ附近。兩個信號被記錄了幾乎9個月左右。每個信號有16384個樣本的長度，并且被記錄的樣本率為20KHZ。信號在圖9中顯示。一系列的沖擊呈現(xiàn)在該波形中，提出一個可能

23、缺陷的初步跡象。第一信號間的脈沖的時間間隔幾乎不變，為0.007秒（圖9a）.這些時間間隔與預(yù)期的BPFO頻率的理論估算相一致。在這個第一信號的光譜上顯示存在兩個頻段，每一個包括特別的，個別的和等間距譜線的頻率間隔大約為142HZ，與自然頻率BPFO相對應(yīng)。(a)(b)圖9 能源發(fā)電廠風(fēng)機軸承上測量的信號波形相反，只有一個激發(fā)頻段出現(xiàn)在第二信號的光譜上，與第一信號的第二頻段相一致(圖10b)。但是，在這種情形下，該段不包含任何主峰，可以代表缺陷的類型。這是由于滑動的程度和頻繁度，它們可以在第二信號的波形圖中觀察到(圖9b)。圖10 在圖9中的信號譜圖11 在圖9中信號的概率密度也有意思的是，第二信號的kurtosis減少，并且它的