李志農(nóng) 李舒揚(yáng) 柳寶 陶俊勇

摘要 相比較于在完整數(shù)據(jù)下設(shè)備性能退化預(yù)測(cè)，缺失數(shù)據(jù)下的預(yù)測(cè)是更加困難的，也是更有意義的。然而，現(xiàn)有的軸承性能退化預(yù)測(cè)方法都未考慮缺失數(shù)據(jù)下的預(yù)測(cè)，基于此，提出了一種基于無限隱馬爾可夫模型的缺失數(shù)據(jù)下軸承退化預(yù)測(cè)方法。在提出的方法中，通過建立無限隱馬爾可夫預(yù)測(cè)模型，預(yù)測(cè)了滾動(dòng)軸承樣本數(shù)據(jù)在振蕩階段所缺失的數(shù)據(jù)點(diǎn)，形成新的完整數(shù)據(jù)。同時(shí)，再使用建立的預(yù)測(cè)模型對(duì)新的完整數(shù)據(jù)進(jìn)行單步預(yù)測(cè)。實(shí)驗(yàn)結(jié)果表明，與真實(shí)值對(duì)比，得到的預(yù)測(cè)數(shù)據(jù)具有較小的平均誤差值；對(duì)比真實(shí)值、完整數(shù)據(jù)下的預(yù)測(cè)值和新的完整數(shù)據(jù)下的預(yù)測(cè)值，驗(yàn)證了提出方法的有效性，能夠反映滾動(dòng)軸承退化的變化趨勢(shì)。提出的方法可為數(shù)據(jù)缺失下滾動(dòng)軸承的退化趨勢(shì)預(yù)測(cè)提供一種思路，具有重要的理論價(jià)值和工程應(yīng)用價(jià)值。

關(guān)鍵詞 故障診斷; 滾動(dòng)軸承; 無限隱馬爾可夫模型（iHMM）; 性能退化; 趨勢(shì)預(yù)測(cè); 缺失數(shù)據(jù)

引 言

滾動(dòng)軸承作為旋轉(zhuǎn)機(jī)械設(shè)備中回轉(zhuǎn)部件的主體，會(huì)因長時(shí)間處于連續(xù)工作狀態(tài)而引發(fā)嚴(yán)重的故障。因此，對(duì)軸承進(jìn)行性能退化趨勢(shì)預(yù)測(cè)是非常有必要的。目前，針對(duì)旋轉(zhuǎn)機(jī)械中軸承退化趨勢(shì)預(yù)測(cè)的研究取得了一些進(jìn)展，國內(nèi)外學(xué)者提出了一些行之有效的方法。王浩然等［1］使用多目標(biāo)和聲搜索算法優(yōu)化了SVM模型，在對(duì)原始特征集深層次特征提取后，以此為退化指標(biāo)建立了最優(yōu)評(píng)估模型。馬海龍［2］考慮了單一特征量預(yù)測(cè)會(huì)導(dǎo)致全局軸承剩余壽命誤差大，推導(dǎo)出了主元融合多個(gè)指標(biāo)的預(yù)測(cè)方法。陳鵬等［3］在文獻(xiàn)［1?2］的基礎(chǔ)上，針對(duì)軸承退化指標(biāo)提取方法中存在提取不全面的問題，建立了全局和局部領(lǐng)域保持嵌入式的支持向量機(jī)描述模型。周建民等［4］采用模糊C均值和SVDD模型對(duì)滾動(dòng)軸承退化狀態(tài)進(jìn)行識(shí)別。姜萬錄等［5］考慮到傳統(tǒng)的故障診斷方法對(duì)數(shù)據(jù)的處理很復(fù)雜，研究了基于平滑Wigner?Ville分布與深度學(xué)習(xí)相結(jié)合的故障診斷方法。Komshin等［6］提出了一種滾動(dòng)軸承運(yùn)行期間計(jì)量保障和技術(shù)狀態(tài)評(píng)估的新方法。Ali等［7］結(jié)合概率神經(jīng)網(wǎng)絡(luò)和簡(jiǎn)化模糊自適應(yīng)共振理論映射對(duì)軸承退化進(jìn)行識(shí)別。Chegini等［8］在軸承失效實(shí)驗(yàn)中引入了周期強(qiáng)度因子、集成經(jīng)驗(yàn)?zāi)B(tài)分解和小波包分解的混合特征來檢測(cè)軸承開始退化的時(shí)間。Chegini等［9］在Chegini等［8］的研究基礎(chǔ)上使用了基于包絡(luò)諧噪比和快速集成經(jīng)驗(yàn)?zāi)B(tài)分解的特征提取方法。Alia等［10］基于非監(jiān)督感知方法提出了自組織映射分類法來跟蹤軸承缺陷的嚴(yán)重程度。

然而，現(xiàn)有的軸承性能退化方法都是建立在完整的數(shù)據(jù)鏈上，忽略了包含缺失數(shù)據(jù)的軸承性能退化情況。在工程項(xiàng)目中，因?qū)嶋H工況環(huán)境多變且惡劣，采集系統(tǒng)會(huì)因此不能得到完整的數(shù)據(jù)集。除此之外，一整套的采集傳輸系統(tǒng)也會(huì)出現(xiàn)各種各樣的意外情況，比如前端傳感器在外部條件變化時(shí)發(fā)生損壞［11］，又或者是數(shù)據(jù)在有線傳輸時(shí)采樣頻率產(chǎn)生波動(dòng)、無線傳輸時(shí)的信號(hào)丟失等，這些情況都會(huì)造成采集數(shù)據(jù)的丟失。對(duì)于不完整的數(shù)據(jù)預(yù)測(cè)是當(dāng)前的一個(gè)難點(diǎn)，也是當(dāng)今迫切需要解決的問題。因此，進(jìn)行缺失數(shù)據(jù)下的退化趨勢(shì)預(yù)測(cè)是非常有必要的。

隱Markov模型（Hidden Markov Model，簡(jiǎn)稱HMM）［12?14］是一種統(tǒng)計(jì)分類模型，在部件退化預(yù)測(cè)領(lǐng)域獲得了一些應(yīng)用。崔昊楊等［15］針對(duì)中國高壓電網(wǎng)中斷路器檢修情況復(fù)雜多變和自身性能衰退的情況，將隱馬爾可夫模型引入到該設(shè)備重要部件的可靠性預(yù)測(cè)中。高山等［16］考慮到如今機(jī)械設(shè)備關(guān)聯(lián)化的大趨勢(shì)，利用全矢譜融合振動(dòng)信號(hào)后，使用HMM進(jìn)行剩余壽命預(yù)測(cè)。Soualhi等［17］為保障軸承的可靠性和安全性，提出基于HMM的診斷模型來預(yù)測(cè)其剩余使用壽命。Liu等［18］提出可描述退化狀態(tài)持續(xù)時(shí)間的改進(jìn)HMM，有效預(yù)測(cè)了刀具磨損的退化狀態(tài)。

然而，基于HMM的預(yù)測(cè)方法在建立模型時(shí)，無法準(zhǔn)確地對(duì)隱狀態(tài)數(shù)目進(jìn)行定義，需要人為設(shè)定隱狀態(tài)數(shù)目，限制了HMM在預(yù)測(cè)中的應(yīng)用。針對(duì)HMM的不足，Beal等［19］提出了無限隱Markov模型（Infinite Hidden Markov Model，iHMM），該模型使用兩級(jí)分層Dirichlet過程來定義非參貝葉斯隱Markov模型。并且通過HDP（Hierarchical Dirichlet Process）集成了隱馬爾可夫模型轉(zhuǎn)移矩陣和發(fā)射矩陣的參數(shù)，使得馬爾可夫模型不用再拘泥于有限多的狀態(tài)和有限的觀察量。因此，無限隱Markov模型所擁有的靈活性和自動(dòng)確定隱狀態(tài)數(shù)的能力要優(yōu)于傳統(tǒng)的隱Markov模型。文獻(xiàn)［20?22］探討了iHMM在機(jī)械故障診斷領(lǐng)域中的應(yīng)用，但是這些方法并未涉及到設(shè)備的性能退化趨勢(shì)預(yù)測(cè)。

基于iHMM的獨(dú)特優(yōu)勢(shì)，本文提出了基于iHMM的缺失數(shù)據(jù)下設(shè)備性能退化預(yù)測(cè)方法，并應(yīng)用到滾動(dòng)軸承性能退化預(yù)測(cè)中。在提出的方法中，以美國智能維護(hù)系統(tǒng)中心（Intelligent Maintenance Systems，IMS）的軸承退化數(shù)據(jù)集為樣本，通過模擬數(shù)據(jù)丟失來展開在缺失數(shù)據(jù)的情況下，滾動(dòng)軸承退化趨勢(shì)的研究。建立起基于缺失數(shù)據(jù)下的無限隱馬爾可夫預(yù)測(cè)模型，通過實(shí)驗(yàn)來進(jìn)一步驗(yàn)證預(yù)測(cè)模型的準(zhǔn)確性。

1 iHMM缺失數(shù)據(jù)退化趨勢(shì)預(yù)測(cè)模型的建立

將隱馬爾可夫模型擴(kuò)展到無限隱馬爾可夫模型時(shí)，作為一種非參數(shù)貝葉斯模型，需要將傳統(tǒng)的Dirichlet過程（Dirichlet Process，DP）拓展為層次Dirichlet過程（Hierarchical Dirichlet Process，HDP）。它的思想最初源自文檔分類的問題，為的是解決在眾多的文檔中共享無窮多個(gè)類的問題［23?24］。HDP一般通過馬爾可夫鏈和蒙特卡羅的采樣算法來構(gòu)造。圖1為HDP的有向圖模型。具體的HDP構(gòu)造如下式所示：

式中 “～”表示服從分布。式（1）表示第一層Dirichlet過程：G0為HDP中第一層的基分布，H為G0的基分布，r為G0的測(cè)度；式（2）表示第二層Dirichlet過程：Gj為HDP中第二層的基分布，以G0為基分布，將α作為Gj的測(cè)度，因此稱為分層Dirichlet過程。再以Gj為先驗(yàn)，產(chǎn)生θji的值，將θji作為一個(gè)函數(shù)F（?）的參數(shù)，將此函數(shù)作為先驗(yàn)產(chǎn)生xji值。

在HMM中，隱狀態(tài)需要人為確定，并且隱狀態(tài)的數(shù)量是有限的，通過引入HDP，利用其性質(zhì)，實(shí)現(xiàn)了自動(dòng)生成HMM中隱狀態(tài)的數(shù)目。HDP?HMM也被稱為無限隱馬爾可夫模型（iHMM）［24］，iHMM的有向模型圖構(gòu)造如圖2所示。

iHMM的數(shù)學(xué)模型，可以用下式表示：

式中 γ為iHMM中的基測(cè)度；β為iHMM中由γ控制產(chǎn)生的基分布；πk為iHMM中的第一層基分布；θk為iHMM中的第二層基分布；it為iHMM中由θk產(chǎn)生的變量；ot為iHMM中由it產(chǎn)生的變量。

通過建立預(yù)測(cè)模型方法進(jìn)行缺失數(shù)據(jù)的補(bǔ)全有兩種方式：第一種方式是通過前面完整數(shù)據(jù)集進(jìn)行模型的訓(xùn)練，然后預(yù)測(cè)第一個(gè)缺失的數(shù)據(jù)，在預(yù)測(cè)完之后，將這個(gè)預(yù)測(cè)的數(shù)據(jù)加入之前的數(shù)據(jù)集，重新進(jìn)行訓(xùn)練，從而完成對(duì)第二個(gè)缺失數(shù)據(jù)的預(yù)測(cè)。依此類推，第三個(gè)、第四個(gè)、…、第N個(gè)都是如此預(yù)測(cè)。這種方法的優(yōu)點(diǎn)是它將預(yù)測(cè)出來的數(shù)據(jù)運(yùn)用到了極致，適用于原始訓(xùn)練數(shù)據(jù)不是很充足的情況。另一種方式是使用沒有缺失數(shù)據(jù)集的數(shù)據(jù)進(jìn)行訓(xùn)練，得到訓(xùn)練的模型，然后通過這個(gè)模型對(duì)缺失數(shù)據(jù)進(jìn)行一個(gè)一個(gè)的預(yù)測(cè)。這種方式適用于大樣本數(shù)據(jù)訓(xùn)練，由于前期使用大量數(shù)據(jù)訓(xùn)練出的模型進(jìn)行預(yù)測(cè)，所以后續(xù)缺失數(shù)據(jù)之間的預(yù)測(cè)值就不會(huì)互相影響，誤差也不會(huì)越來越大。這里，將采用第一種退化趨勢(shì)預(yù)測(cè)方式。具體的流程如圖3所示。

由圖3可知，基于iHMM的滾動(dòng)軸承退化趨勢(shì)預(yù)測(cè)的步驟如下：

1） 獲得振動(dòng)信號(hào)。

2） 將振動(dòng)信號(hào)進(jìn)行Hilbet變換預(yù)處理。

3） 對(duì)振動(dòng)信號(hào)進(jìn)行特征提取。在這里，我們使用小波熵作為振動(dòng)信號(hào)的退化指標(biāo)，小波熵的計(jì)算過程如下式所示：

式中 Ej為第j尺度下小波分量能量的平均和；E為所有尺度分量的總能量；pi為歸一化后的能量分布；HE為小波熵。在滾動(dòng)軸承預(yù)測(cè)中共計(jì)獲得984個(gè)小波熵值。

4） 進(jìn)行模型訓(xùn)練，第一步是先使用原始樣本進(jìn)行iHMM退化模型訓(xùn)練，然后預(yù)測(cè)出第一個(gè)缺失值。

5） 第二步是將這個(gè)預(yù)測(cè)出的缺失值與原始樣本合并成一個(gè)新的樣本數(shù)據(jù)集，使用這個(gè)樣本數(shù)據(jù)集進(jìn)行第二個(gè)缺失數(shù)據(jù)的預(yù)測(cè)。

6） 依此類推，預(yù)測(cè)出一定數(shù)量的缺失數(shù)據(jù)。

7） 將預(yù)測(cè)出來的缺失數(shù)據(jù)與原數(shù)據(jù)對(duì)比評(píng)價(jià)。

2 實(shí)驗(yàn)研究

本實(shí)驗(yàn)采用來自IMS的數(shù)據(jù)。實(shí)驗(yàn)所需的設(shè)備安裝了4個(gè)ZA?2115雙列滾子軸承，并使用2000 r/min的傳動(dòng)帶驅(qū)動(dòng)轉(zhuǎn)軸恒定轉(zhuǎn)速轉(zhuǎn)動(dòng)。在X軸與Y軸方向分別安裝了加速度傳感器。該數(shù)據(jù)共包含三個(gè)數(shù)據(jù)集，每一個(gè)數(shù)據(jù)集都描述了軸承從測(cè)試到失效的過程。此處采用第二個(gè)數(shù)據(jù)集。在數(shù)據(jù)集中：采樣頻率為20 kHz，采樣間隔為10 min，每隔10 min可采得一個(gè)文件長度為20480個(gè)點(diǎn)的數(shù)據(jù)點(diǎn)集，共計(jì)984組這樣的點(diǎn)集。該實(shí)驗(yàn)共持續(xù)164 h。

該軸承實(shí)驗(yàn)裝置如圖4所示。

圖5為軸承振動(dòng)信號(hào)經(jīng)小波熵特征提取后的984組值。由圖5可知，軸承在0～85.5 h（小波熵提取值為文件第1～512）處于一個(gè)平穩(wěn)過程；在85.5～115 h（小波熵提取值為文件第513～692）處于振蕩過程；在115 h之后（小波熵提取值為文件第693～984）處于故障情況。

同時(shí)加入了峭度作為退化特征信號(hào)用來對(duì)比小波熵，從圖6中可以看出，小波熵在85.5 h后，信號(hào)發(fā)生了明顯的變化，而峭度信號(hào)則在較晚的時(shí)間發(fā)生變化，證明了小波熵在提取軸承退化特征信號(hào)時(shí)的優(yōu)越性。平穩(wěn)時(shí)期的數(shù)據(jù)點(diǎn)不需要對(duì)其進(jìn)行預(yù)測(cè)，同理，故障時(shí)期的點(diǎn)已經(jīng)不再具備預(yù)測(cè)的可能性。而選取振蕩階段經(jīng)小波熵處理后第513～692點(diǎn)的數(shù)據(jù)點(diǎn)開展后續(xù)的缺失數(shù)據(jù)預(yù)測(cè)更有意義。因此，將選擇振蕩過程中第513～642點(diǎn)進(jìn)行后續(xù)缺失數(shù)據(jù)情況下的滾動(dòng)軸承退化趨勢(shì)預(yù)測(cè)實(shí)驗(yàn)。

因?yàn)閕HMM建模需要的觀察值是離散的變量，所以要將小波熵進(jìn)行量化處理，在這里采用的量化處理方法是Lloyds編碼。Lloyds編碼可以實(shí)現(xiàn)小波熵的標(biāo)量化，可以使后續(xù)的結(jié)果經(jīng)處理后在一個(gè)范圍內(nèi)變化。將Lloyds編碼的值設(shè)定為40。將所選取振蕩過程中的小波熵進(jìn)行離散化，得到的結(jié)果如圖7所示。

由圖7可知，經(jīng)Lloyds編碼之后的小波熵特征值在0～40的范圍內(nèi)，在85.5～115 h（編碼后小波熵提取值對(duì)應(yīng)為第513～692點(diǎn)之間）很好地反映了軸承在退化階段的變化情況。

2.1 進(jìn)行iHMM軸承退化數(shù)據(jù)訓(xùn)練

在這個(gè)數(shù)據(jù)段內(nèi)，將第513～613點(diǎn)區(qū)間內(nèi)共計(jì)101個(gè)點(diǎn)作為缺失數(shù)據(jù)情況下iHMM的訓(xùn)練數(shù)據(jù)，缺失數(shù)據(jù)段設(shè)定為第614～623點(diǎn)，這10個(gè)點(diǎn)可以較好反映軸承在產(chǎn)生退化初期時(shí)所表達(dá)的意義。振蕩過程中iHMM訓(xùn)練的相關(guān)參數(shù)如圖8～13所示。

圖8為經(jīng)iHMM共計(jì)2000次迭代后的隱狀態(tài)K的收斂軌跡，表示iHMM訓(xùn)練過程中軸承退化數(shù)據(jù)隨時(shí)間變化逐漸收斂至固定的值。從圖8中可以看出，在經(jīng)過2000次的迭代之后，隱狀態(tài)K的值逐漸趨于穩(wěn)定，基本穩(wěn)定在2個(gè)或者3個(gè)狀態(tài)。也證實(shí)了經(jīng)訓(xùn)練后的無限隱Markov模型可以自適應(yīng)地確定隱狀態(tài)數(shù)。圖9為iHMM迭代過程中超參數(shù)α的收斂情況，α控制軸承實(shí)際退化狀態(tài)向自身轉(zhuǎn)移的先驗(yàn)趨勢(shì)。從圖9中可以看出，α值在訓(xùn)練中表現(xiàn)得較穩(wěn)定，表明在訓(xùn)練迭代的2000次過程中，自轉(zhuǎn)移先驗(yàn)值變化不大。圖10為iHMM迭代過程中超參數(shù)γ的收斂情況，γ控制觀測(cè)狀態(tài)對(duì)應(yīng)的隱狀態(tài)數(shù)目。從圖10中可以看出，γ值從迭代開始的較大值到迭代200次之后趨于一個(gè)穩(wěn)定的范圍值，間接對(duì)應(yīng)了圖8中隱狀態(tài)數(shù)目的跳變從最開始較大變化到后來較小變化。圖11為iHMM訓(xùn)練中對(duì)振蕩狀態(tài)軸承特征數(shù)據(jù)聯(lián)合似然對(duì)數(shù)值。從圖11中可以看出，在迭代的過程中，聯(lián)合似然對(duì)數(shù)也趨于穩(wěn)定。圖12為iHMM訓(xùn)練結(jié)束后產(chǎn)生的狀態(tài)轉(zhuǎn)移矩陣，狀態(tài)轉(zhuǎn)移矩陣表示軸承退化特征數(shù)據(jù)由一個(gè)狀態(tài)向另一個(gè)狀態(tài)轉(zhuǎn)移所呈現(xiàn)的轉(zhuǎn)移概率矩陣，圖12的轉(zhuǎn)移矩陣表示在軸承退化過程中，選定過程的隱狀態(tài)是以2×2的形式轉(zhuǎn)移變化的。圖13為iHMM訓(xùn)練結(jié)束后產(chǎn)生的發(fā)射矩陣。發(fā)射矩陣表示觀測(cè)到的軸承數(shù)據(jù)與實(shí)際軸承狀態(tài)狀況之間產(chǎn)生的概率，其表示在確定了狀態(tài)轉(zhuǎn)移矩陣為2個(gè)狀態(tài)之間互相轉(zhuǎn)移的前提下，實(shí)際觀測(cè)到的序列在每個(gè)狀態(tài)下都有20個(gè)觀測(cè)序列值，即觀測(cè)序列與隱狀態(tài)序列在每一個(gè)時(shí)刻存在2×20個(gè)概率。

2.2 缺失軸承數(shù)據(jù)預(yù)測(cè)

在這個(gè)數(shù)據(jù)段內(nèi)，將第513～613點(diǎn)區(qū)間內(nèi)共計(jì)101個(gè)點(diǎn)作為缺失數(shù)據(jù)情況下iHMM的訓(xùn)練數(shù)據(jù)，缺失數(shù)據(jù)段設(shè)定為第614～623點(diǎn)。通過訓(xùn)練好的iHMM可以預(yù)測(cè)出第614點(diǎn)。將第614點(diǎn)和第513～613點(diǎn)合并為第513～614點(diǎn)再次進(jìn)行iHMM訓(xùn)練用來預(yù)測(cè)第615點(diǎn)，如此往復(fù)循環(huán)，將缺失數(shù)據(jù)段全部都預(yù)測(cè)出來，預(yù)測(cè)結(jié)果如圖14和15所示，為了更好地與前后數(shù)據(jù)進(jìn)行對(duì)比，在圖中還加入了預(yù)測(cè)缺失數(shù)據(jù)前后的5個(gè)真實(shí)數(shù)據(jù)點(diǎn)作為對(duì)比。

從圖15中可以看出，雖然預(yù)測(cè)出來的十個(gè)數(shù)據(jù)（點(diǎn)614～623）的值與真實(shí)值都有一定的差距，但是可以反映出數(shù)據(jù)的變化趨勢(shì)。為了量化退化過程的趨勢(shì)，引入平均相對(duì)誤差（MRE）、均方根誤差（RMSE）和平均絕對(duì)誤差（MAE）：

將這10個(gè)預(yù)測(cè)點(diǎn)的數(shù)據(jù)代入到式（13）～（15）中，得出MRE=0.1596，RMSE=10.1193，MAE=3.8000。結(jié)果表明，使用iHMM進(jìn)行滾動(dòng)軸承的退化趨勢(shì)預(yù)測(cè)可以很好反映特征信號(hào)的變化。

在將這十個(gè)缺失數(shù)據(jù)補(bǔ)全之后，歸類到一個(gè)完整的數(shù)據(jù)集中，重新進(jìn)行滾動(dòng)軸承的退化趨勢(shì)預(yù)測(cè)。本次預(yù)測(cè)目的是使用最初完整樣本數(shù)據(jù)集和通過預(yù)測(cè)缺失數(shù)據(jù)補(bǔ)全的數(shù)據(jù)集分別進(jìn)行退化趨勢(shì)預(yù)測(cè)實(shí)驗(yàn)對(duì)比。在后續(xù)的預(yù)測(cè)過程中也使用本文提出的方法，預(yù)測(cè)出一個(gè)值后，將這個(gè)值代入預(yù)測(cè)模型，重新訓(xùn)練，再預(yù)測(cè)下一個(gè)值。進(jìn)行20次的預(yù)測(cè)后，如圖16和17所示。

將缺失數(shù)據(jù)下的預(yù)測(cè)值與真實(shí)值比較：MRE=0.0073，RMSE=1.5811，MAE=1.7000；將完整數(shù)據(jù)下的預(yù)測(cè)值與真實(shí)值比較：MRE=0.0086，RMSE=2.2136，MAE=0.9000。

由圖17可知，通過三條折線的比較，得出使用iHMM在缺失數(shù)據(jù)補(bǔ)全情況下和原始數(shù)據(jù)集都能較好地反映滾動(dòng)軸承的退化趨勢(shì)。不同之處在于，缺失數(shù)據(jù)的iHMM預(yù)測(cè)的效果與完整數(shù)據(jù)iHMM還是存在一定的差距。并且在預(yù)測(cè)步數(shù)的末尾處與真實(shí)值的趨勢(shì)存在差異，然而，這種差異并沒有影響到軸承性能退化的預(yù)測(cè)趨勢(shì)。

3 總 結(jié)

缺失數(shù)據(jù)下的軸承性能退化的預(yù)測(cè)是一直被忽略的一個(gè)關(guān)鍵性問題。本文基于無限隱Marcov模型的獨(dú)特優(yōu)勢(shì)，即可以自適應(yīng)地確定HMM預(yù)測(cè)模型的隱狀態(tài)數(shù)目，提出了基于無限隱Marcov模型的軸承性能退化預(yù)測(cè)方法，并應(yīng)用到缺失數(shù)據(jù)下的軸承性能趨勢(shì)預(yù)測(cè)中。在提出的方法中，通過對(duì)滾動(dòng)軸承的振動(dòng)信號(hào)提取特征，結(jié)合無限隱Marcov模型對(duì)其進(jìn)行在缺失數(shù)據(jù)狀態(tài)下的數(shù)據(jù)預(yù)測(cè)和補(bǔ)全缺失數(shù)據(jù)情況下的退化趨勢(shì)研究。將獲取到的軸承振動(dòng)信號(hào)預(yù)處理后再進(jìn)行退化特征提取，輸入到建立好的iHMM退化趨勢(shì)預(yù)測(cè)模型中訓(xùn)練。結(jié)果表明：針對(duì)軸承故障數(shù)據(jù)在傳輸中可能丟失的情況，采用無限隱Markov模型對(duì)退化趨勢(shì)進(jìn)行單步預(yù)測(cè)，得到的數(shù)據(jù)結(jié)果與原始數(shù)據(jù)相比具有較高的精度。將補(bǔ)全后的數(shù)據(jù)整合至完整數(shù)據(jù)集后，通過iHMM再訓(xùn)練與重新預(yù)測(cè)，很好地反映了數(shù)據(jù)集的退化趨勢(shì)。提出的方法同時(shí)也為解決滾動(dòng)軸承缺失數(shù)據(jù)的預(yù)測(cè)開拓了一種思路，在理論研究和工程應(yīng)用中具有重要的理論意義和工程應(yīng)用價(jià)值。

參考文獻(xiàn)

1王浩然， 袁小芳， 田爭(zhēng)鳴， 等. 基于多目標(biāo)和聲搜索優(yōu)化SVM的伺服電機(jī)滾動(dòng)軸承性能退化評(píng)估［J］. 機(jī)床與液壓， 2020， 48（16）： 176-181.

Wang Haoran， Yuan Xiaofang， Tian Zhengming， et al. Performance degradation evaluation of servo motor rolling bearing based on multi-objective harmony search optimization SVM［J］. Machine Tool & Hydraulics， 2020， 48（16）： 176-181.

2馬海龍. 基于主元特征融合和SVM的軸承剩余壽命預(yù)測(cè)［J］. 工礦自動(dòng)化， 2019， 45（8）： 74-78.

Ma Hailong. Prediction of bearing remaining life based on principal component feature fusion and SVM［J］. Industry and Mine Automation， 2019， 45（8）： 74-78.

3陳鵬， 趙小強(qiáng). 基于GLNPE-SVDD的滾動(dòng)軸承性能退化評(píng)估方法［J］. 華中科技大學(xué)學(xué)報(bào) （自然科學(xué)版）， 2021， 49（1）： 12-16.

Chen Peng， Zhao Xiaoqiang. Evaluation method of rolling bearing performance degradation based on GLNPE-SVDD［J］. Journal of Huazhong University of Science and Technology （Natural Science Edition）， 2021， 49（1）： 12-16.

4周建民， 游濤， 尹文豪， 等. 基于融合FCM-SVDD模型的滾動(dòng)軸承退化狀態(tài)識(shí)別［J］. 機(jī)械設(shè)計(jì)與研究， 2020， 36（1）： 124-129.

Zhou Jianming， You Tao， Yin Wenhao， et al. Recognition of degradation state of rolling bearing based on fusion FCM-SVDD model［J］. Machine Design and Research， 2020， 36（1）： 124-129.

5姜萬錄， 李振寶， 雷亞飛， 等. 基于深度學(xué)習(xí)的滾動(dòng)軸承故障診斷與性能退化程度識(shí)別方法［J］. 燕山大學(xué)學(xué)報(bào)， 2020， 44（6）： 526-536.

Jiang Wanlu， Li Zhenbao， Lei Yafei， et al. Deep learning-based fault diagnosis and performance degradation identification method of rolling bearing［J］. Journal of Yanshan University， 2020， 44（6）： 526-536.

6Komshin A S， Potapov K G， Pronyakin V I， et al. Operational effectiveness of phase-chronometric and neurodiagnostic methods for controlling rolling-element bearing degradation［J］. Measurement Techniques， 2020， 63（7）： 559-566.

7Ali J B， Saidi L， Mouelhi A， et al. Linear feature selection and classification using PNN and SFAM neural networks for a nearly online diagnosis of bearing naturally progressing degradations［J］. Engineering Applications of Artificial Intelligence， 2015， 42： 67-81.

8Chegini S N， Manjili M J H， Bagheri A， et al. New bearing slight degradation detection approach based on the periodicity intensity factor and signal processing methods［J］. Measurement， 2021， 170（1）： 108696.

9Chegini S N， Manjili M J H， Bagheri A. New fault diagnosis approaches for detecting the bearing slight degradation［J］. Meccanica， 2020， 55（1）： 261-286.

10Alia S， Nasri R， Meddour I， et al. Comparison between sound perception and self-organizing maps in the monitoring of the bearing degradation［J］. The International Journal of Advanced Manufacturing Technology， 2020， 110（7）： 2003-2013.

11鄭明亮. 構(gòu)皮灘電廠#1機(jī)氣隙傳感器損壞原因分析與處理［J］. 裝備維修技術(shù)， 2020， 176（2）：133.

Zheng Mingliang. Goupitan power plant #1. analysis and treatment of damage causes of air gap sensor of engine［J］. Equipment Technology， 2020， 176（2）： 133.

12Baum L E， Eagon J A. An inequality with applications to statistical estimation for probabilistic functions of Markov processes and to a model for ecology［J］. Bulletin of the American Mathematical Society， 1967， 73（3）： 360-363.

13Baum L E， Petrie T. Statistical inference for probabilistic functions of finite state Markov Chains［J］. The Annals of Mathematical Statistics， 1966， 37（6）： 1554-1563.

14Baum L E， Petrie T， Soules G， et al. A maximization technique occurring in the statistical analysis of probabilistic functions of Markov Chains［J］. The Annals of Mathematical Statistics， 1970， 41（1）： 164-171.

15崔昊楊， 夏晟， 周坤， 等. 基于Moffat休息型HMM退化過程的220 kV斷路器可靠性預(yù)測(cè)［J］. 高電壓技術(shù)， 2021， 47（6）： 2108-2116.

Cui Haoyang， Xia Shen， Zhou Kun， et al. Reliability prediction of 220 kV circuit breaker based on the degradation process of Moffat resting HMM［J］. High Voltage Technology， 2021， 47（6）： 2108-2116.

16高山， 周玉平， 陳宏， 等. 全矢HMM在軸承剩余壽命預(yù)測(cè)中的應(yīng)用［J］. 機(jī)械設(shè)計(jì)與制造， 2020，12（12）： 64-67.

Gao Shan， Zhou Yuping， Chen Hong， et al. The application of full vector HMM in the prediction of bearing remaining life［J］. Machinery Design & Manufacture， 2020，12（12）： 64-67.

17Soualhi A， Clerc G， Razik H， et al. Hidden Markov models for the prediction of impending fault[J]. IEEE Transactions on Industrial Electronics， 2016， 63（5）： 3271?3281.

18Liu T S， Zhu K P， Zeng L C， et al. Diagnosis and prognosis of degradation process via hidden semi?Markov model[J]. IEEE?ASME Transactions on Mechatronics， 2018， 23（3）：1456?1466.

19Beal M J， Ghahramani Z， Rasmussen C E. The infinite hidden Markov model［C］. Advances in Neural Information Processing Systems. Vancouver， British Columbia， Canada， 2001： 14.

20李志農(nóng)， 柳寶. 無限隱Markov模型理論及仿真研究［J］. 南昌航空大學(xué)學(xué)報(bào) （自然科學(xué)版）， 2016， 30（2）： 37-43.

Li Zhinong， Liu Bao. Infinite hidden Markov model theory and simulation research［J］. Journal of Nanchang Hangkong University （Natural Sciences）， 2016， 30（2）： 37-43.

21李志農(nóng)， 柳寶， 侯娟. 基于無限隱Markov模型的旋轉(zhuǎn)機(jī)械故障診斷方法研究［J］. 儀器儀表學(xué)報(bào)， 2016， 37（10）： 2185-2192.

Li Zhinong， Liu Bao， Hou Juan. Research on rotating machinery fault diagnosis method based on infinite hidden Markov model［J］. Chinese Journal of Scientific Instrument， 2016， 37（10）： 2185-2192.

22李志農(nóng)， 熊俊偉. 基于無限因子隱Markov模型的旋轉(zhuǎn)機(jī)械故障識(shí)別方法［J］. 失效分析與預(yù)防， 2016， 11（3）： 133-138.

Li Zhinong， Xiong Junwei. Rotating machinery fault identification method based on infinite factor hidden Markov model［J］. Failure Analysis and Prevention， 2016， 11（3）： 133-138.

23周志敏， 高申勇. 分層Dirichlet過程原理及應(yīng)用綜述［J］. 計(jì)算機(jī)應(yīng)用與軟件， 2014， 31（8）： 1-5.

Zhou Zhimin， Gao Shenyong. Summary of the principle and application of hierarchical Dirichlet process［J］. Computer Applications and Software， 2014， 31（8）： 1-5.

24周建英， 王飛躍， 曾大軍. 分層Dirichlet過程及其應(yīng)用綜述［J］. 自動(dòng)化學(xué)報(bào)， 2011， 37（4）： 389-407.

Zhou Jianying， Wang Feiyue， Zeng Dajun. Summary of hierar-chical dirichlet process and its application［J］. Acta Automatica Sinica， 2011， 37（4）： 389-407.

Infinite hidden Markov model and its application in the prediction of bearing degradation trends with missing data

LI Zhi-nong 1 ?LI Shu-yang 1LIU Bao 1TAO Jun-yong 2

1. Key Laboratory of Non-destructive Testing Technology of Ministry of Education， Nanchang Hangkong University， Nanchang 330063， China；

2. Key Laboratory of Equipment Comprehensive Support Technology， National University of Defense Technology， Changsha 410073， China

Abstract Compared with the equipment performance degradation prediction under the complete data， the prediction under the missing data is more difficult and more meaningful. However， the existing prediction methods of bearing performance degradation do not consider the prediction under missing data. Based on the above problem， a bearing degradation prediction method based on infinite hidden Markov model （iHMM） is proposed under the missing data. In the proposed method， an iHMM prediction model with wavelet entropy as the degradation feature is established to predict the missing data points of rolling bearing sample data and form new complete data. Then the proposed prediction model is used to make single-step predictions on the new complete data. The experiment results show that compared with the real value， the obtained prediction data has a smaller average error. Compare the real value， the predicted value under the complete data， and the predicted value under the new complete data， the prediction data obtained by the iHMM prediction model can also well reflect the degradation trend of rolling bearing. The proposed method can provide a new idea for predicting the degradation trend of rolling bearings under the missing data. Therefore， the proposed method has important theoretical value and engineering application value.

Keywords fault diagnosis; rolling bearing; infinite hidden Markov model （iHMM）; performance degradation; trend forecast; missing data