文 娟， 高宏力

(西南交通大學(xué) 機(jī)械工程學(xué)院 先進(jìn)驅(qū)動(dòng)節(jié)能技術(shù)教育部工程研究中心，成都 610031)

作為一種關(guān)鍵零部件，軸承被廣泛應(yīng)用于各種機(jī)械系統(tǒng)中。通常軸承的工況復(fù)雜，工作環(huán)境惡劣，因而頻繁出現(xiàn)各種故障。一旦軸承出現(xiàn)故障，可能會(huì)造成整個(gè)機(jī)械系統(tǒng)停機(jī)，帶來(lái)災(zāi)難。因此，適時(shí)監(jiān)測(cè)軸承狀態(tài)，并根據(jù)狀態(tài)監(jiān)測(cè)結(jié)果指導(dǎo)維修計(jì)劃對(duì)于提高系統(tǒng)可靠性與降低生產(chǎn)成本具有重要意義。狀態(tài)維修(Condition-Based Maintenance)主要包括兩個(gè)部分，故障診斷和故障預(yù)警[1-3]。故障診斷的主要任務(wù)為判斷系統(tǒng)能否正常工作，一旦不能工作，將會(huì)停機(jī)并更換故障零部件。故障預(yù)警的主要任務(wù)為判斷系統(tǒng)的當(dāng)前狀態(tài)，并根據(jù)當(dāng)前狀態(tài)對(duì)系統(tǒng)的剩余服役壽命進(jìn)行預(yù)測(cè)，并使用預(yù)測(cè)結(jié)果指導(dǎo)維修計(jì)劃的制定。準(zhǔn)確的故障預(yù)警能夠在確保安全的情況下節(jié)約維修費(fèi)用，因而受到了越來(lái)越多研究者的關(guān)注。

通常壽命預(yù)測(cè)可以分為兩種，基于數(shù)據(jù)驅(qū)動(dòng)的方法和基于模型的方法。數(shù)據(jù)驅(qū)動(dòng)的方法試圖通過(guò)機(jī)器學(xué)習(xí)方法建立狀態(tài)監(jiān)測(cè)數(shù)據(jù)與系統(tǒng)退化過(guò)程之間的關(guān)系。Tian等[4]提出了一種利用失效數(shù)據(jù)和未失效懸疑數(shù)據(jù)建立人工神經(jīng)網(wǎng)絡(luò)(Articial Neural Network)壽命預(yù)測(cè)模型的方法。Tran等[5]利用時(shí)間序列技術(shù)和支持向量機(jī)(Support Vector Machines)對(duì)剩余壽命進(jìn)行預(yù)測(cè)。與數(shù)據(jù)驅(qū)動(dòng)方法不同，基于模型的方法主要通過(guò)建立能夠反應(yīng)退化過(guò)程的數(shù)學(xué)模型來(lái)對(duì)剩余壽命進(jìn)行預(yù)測(cè)。Gebraeel等[6-8]在數(shù)學(xué)模型的建立方面完成了許多工作。模型建立后，通常需要利用貝葉斯推斷理論結(jié)合狀態(tài)監(jiān)測(cè)數(shù)據(jù)估計(jì)系統(tǒng)狀態(tài)和模型參數(shù)，并完成壽命預(yù)測(cè)。

傳統(tǒng)的實(shí)現(xiàn)貝葉斯推斷的方法有卡爾曼濾波(Kalman filter, KF)和擴(kuò)展卡爾曼濾波(Extended Kalman filter, EKF)，都已成功應(yīng)用到壽命預(yù)測(cè)中[9]。粒子濾波(Particle filter, PF)是一種利用一群具有特定權(quán)值的粒子近似狀態(tài)估計(jì)值的方法。與KF方法不同，PF算法不受限于高斯假設(shè)，特別適用于非線性非高斯應(yīng)用場(chǎng)景，已經(jīng)在壽命預(yù)測(cè)領(lǐng)域得到了許多應(yīng)用[10-11]。但是，PF算法存在粒子退化的缺陷，在經(jīng)歷了一定的迭代步數(shù)后，許多粒子的權(quán)值變得極小。而粒子的退化會(huì)影響剩余壽命預(yù)測(cè)的準(zhǔn)確度[12]，因此，如何有效降低PF方法中的粒子退化程度，進(jìn)而減小預(yù)測(cè)誤差，是基于PF的壽命預(yù)測(cè)方法中存在的一大問(wèn)題。選擇合適的重要性采樣概率分布函數(shù)是解決這一問(wèn)題的有效方法。無(wú)跡粒子濾波(Unscented Particle filter, UPF)就是一種利用這種方式解決粒子退化問(wèn)題的改進(jìn)粒子濾波算法[13]。UPF結(jié)合了PF和無(wú)跡卡爾曼濾波(Unscented Kalman filter, UKF)的優(yōu)點(diǎn)，利用UKF算法獲得一個(gè)合適的重要性采樣概率分布，從而有效降低粒子退化程度。為了解決基于PF軸承剩余壽命預(yù)測(cè)中的粒子退化問(wèn)題，本文提出一種基于UPF的軸承剩余壽命預(yù)測(cè)方法，即通過(guò)UPF結(jié)合數(shù)學(xué)模型與軸承的狀態(tài)監(jiān)測(cè)信號(hào)估計(jì)軸承的當(dāng)前健康狀態(tài)，并給出其剩余壽命的估計(jì)值。

1 UPF算法

1.1 貝葉斯估計(jì)

大部分的動(dòng)態(tài)系統(tǒng)都可以由兩個(gè)方程描述，狀態(tài)轉(zhuǎn)移方程和觀測(cè)方程為

θt=ft(θt-1,nt-1)

(1)

yt=ht(θt,vt)

(2)

式中：ft—Rnθ×Rnn→Rnθ為上一時(shí)刻到當(dāng)前時(shí)刻的狀態(tài)轉(zhuǎn)移函數(shù)；ht—Rnθ×Rnv→Rny為當(dāng)前狀態(tài)θt的觀測(cè)函數(shù)；nt為獨(dú)立同分布過(guò)程噪聲；vt為獨(dú)立同分布測(cè)量噪聲。

分別使用θ0：t={θ0,…,θt}和y1:t={y1,…,yt}代替當(dāng)前所有狀態(tài)和所有觀測(cè)值。我們的目的是得到后驗(yàn)概率p(θ1:t|y1:t)。根據(jù)馬爾科夫過(guò)程特性，我們無(wú)需追蹤所有已有狀態(tài)，只需遞歸地計(jì)算p(θt|y1：t)便可得到后驗(yàn)概率。具體地，每一時(shí)刻的先驗(yàn)概率密度分布函數(shù)和后驗(yàn)概率密度分布函數(shù)分別如式(3)和式(4)所示

(3)

(4)

1.2 UPF算法

PF算法的本質(zhì)是利用一群具有特定權(quán)值的粒子來(lái)近似系統(tǒng)狀態(tài)，通過(guò)更新粒子權(quán)值和粒子來(lái)實(shí)現(xiàn)最優(yōu)估計(jì)。具體地，系統(tǒng)狀態(tài)后驗(yàn)概率密度分布函數(shù)可以表示為

(5)

標(biāo)準(zhǔn)PF算法中，通常采用系統(tǒng)的先驗(yàn)分布作為重要性采樣分布，這種做法能夠簡(jiǎn)化計(jì)算過(guò)程，但是容易使PF算法陷入粒子退化問(wèn)題中。與標(biāo)準(zhǔn)PF算法不同，UPF利用UKF生成重要性采樣分布，充分利用最新的觀測(cè)值，能夠有效地解決粒子退化問(wèn)題。

UPF的具體實(shí)現(xiàn)步驟如下：

步驟3根據(jù)式(6)利用當(dāng)前時(shí)刻t的測(cè)量值更新各粒子所對(duì)應(yīng)的權(quán)值。

(6)

并根據(jù)式(7)對(duì)權(quán)值進(jìn)行歸一化處理。

(7)

步驟4重采樣。為了增加有效粒子的數(shù)量，提高計(jì)算效率，需要移除權(quán)值極小的粒子，并復(fù)制權(quán)值較大的粒子，這個(gè)過(guò)程稱之為重采樣。一個(gè)簡(jiǎn)單的實(shí)現(xiàn)過(guò)程如下[14]：

(8)

步驟5狀態(tài)估計(jì)。利用更新后的粒子及其權(quán)值估計(jì)當(dāng)前時(shí)刻t的狀態(tài)

(9)

步驟6然后返回至步驟2，重復(fù)之后的步驟，直到t=T(T為測(cè)量值的數(shù)量)。

2 基于UPF的軸承剩余壽命預(yù)測(cè)方法

軸承的整個(gè)工作周期可以分為三個(gè)階段，正常階段，退化階段以及失效階段。當(dāng)軸承處于失效階段時(shí)，意味著軸承已經(jīng)不能正常工作，需要進(jìn)行更換。因此，本文只考慮前兩個(gè)階段的剩余壽命預(yù)測(cè)。軸承的狀態(tài)監(jiān)測(cè)數(shù)據(jù)有振動(dòng)信號(hào)、溫度和聲發(fā)射等。其中，由于振動(dòng)信號(hào)采集方便，而且對(duì)軸承的性能退化十分敏感，因而被廣泛應(yīng)用于軸承的故障診斷和剩余壽命預(yù)測(cè)中[15]。因此，本文從振動(dòng)信號(hào)中提取特征，將其作為軸承性能退化的指標(biāo)。

本文提出的方法如圖1所示。根據(jù)已知的歷史失效軸承數(shù)據(jù)利用隨機(jī)過(guò)程模型建立軸承的退化模型。在針對(duì)具體軸承進(jìn)行壽命預(yù)測(cè)時(shí)，首先判斷其是否開(kāi)始退化，當(dāng)檢測(cè)到其開(kāi)始退化時(shí)，利用UPF方法結(jié)合數(shù)學(xué)模型和測(cè)量數(shù)據(jù)，更新模型參數(shù)，估計(jì)軸承的退化狀態(tài)，并實(shí)現(xiàn)軸承的剩余壽命預(yù)測(cè)。通常，信號(hào)的不同特征描述的是信號(hào)的不同方面，因而能從不同角度反映軸承的退化狀態(tài)。在軸承的故障診斷與壽命預(yù)測(cè)中，峭度對(duì)于早期故障比較敏感，但對(duì)于故障的發(fā)展惡化并不敏感，而均方根值(Root Mean Square, RMS)反映的是振動(dòng)能量隨故障發(fā)展的增長(zhǎng)，因而對(duì)于軸承故障發(fā)展程度較為敏感。因此，本文利用峭度監(jiān)測(cè)故障的出現(xiàn)，采用RMS評(píng)估軸承的故障發(fā)展程度，并預(yù)測(cè)其剩余壽命。

圖1 基于UPF的軸承壽命預(yù)測(cè)方法流程圖Fig.1 Flowchart of the proposed method

2.1 軸承退化檢測(cè)

在軸承的整個(gè)工作周期中，軸承有很長(zhǎng)一段時(shí)間都處于正常工作狀態(tài)。在這個(gè)階段，我們沒(méi)有必要預(yù)測(cè)軸承的剩余壽命。因此，檢測(cè)軸承的退化起點(diǎn)，僅在軸承的性能退化階段預(yù)測(cè)其剩余壽命能夠有效節(jié)約計(jì)算資源。此外，正確區(qū)分正常狀態(tài)與退化狀態(tài)的監(jiān)測(cè)數(shù)據(jù)，能夠有效避免正常狀態(tài)數(shù)據(jù)的干擾，提高壽命預(yù)測(cè)的準(zhǔn)確度。因此，本文在針對(duì)具體軸承進(jìn)行壽命預(yù)測(cè)時(shí)，先判斷其工作狀態(tài)，當(dāng)檢測(cè)到其開(kāi)始退化時(shí)，才開(kāi)始預(yù)測(cè)其剩余壽命。

軸承振動(dòng)信號(hào)的峭度是一個(gè)無(wú)量綱參數(shù)，與軸承的負(fù)載及參數(shù)無(wú)關(guān)。通常，峭度指數(shù)對(duì)于軸承的早期故障比較敏感，但對(duì)于故障的發(fā)展與狀態(tài)惡化并不敏感。因此，本文采用振動(dòng)信號(hào)的峭度值來(lái)判斷故障的出現(xiàn)，即確定軸承剩余壽命預(yù)測(cè)的起始點(diǎn)。首先，計(jì)算軸承正常工作狀態(tài)下振動(dòng)信號(hào)峭度值的均值μ與標(biāo)準(zhǔn)差σm，定義軸承狀態(tài)正常的峭度值區(qū)間為[μ-3σm,μ+3σm]。狀態(tài)監(jiān)測(cè)過(guò)程中，當(dāng)tf時(shí)刻的峭度值mf超出該區(qū)間時(shí)，則認(rèn)為軸承振動(dòng)信號(hào)出現(xiàn)異常。由于振動(dòng)信號(hào)的采集以及特征提取中存在一定的隨機(jī)誤差，會(huì)造成誤判。為了消除這一影響，引入Li等研究中的觸發(fā)機(jī)制來(lái)判定軸承的退化起始點(diǎn)，具體步驟如下：

步驟1首先，定義l= 0，當(dāng)峭度值第一次超出3σm時(shí)，定義該時(shí)刻點(diǎn)為FPT0；

步驟2令l=l+1，當(dāng)時(shí)間點(diǎn)tf滿足連續(xù)l+1個(gè)時(shí)間點(diǎn)的峭度值{mf+k}k=0:l滿足{|mf+k-μ|>3σm}k=0:l，則定義時(shí)刻tf為FPTl；

步驟3使l由1逐漸增大直到l滿足FPTl=FPTl-1，則認(rèn)為軸承此時(shí)出現(xiàn)故障，并定義FPTl為軸承的剩余壽命預(yù)測(cè)起始點(diǎn)。

通常，由隨機(jī)噪聲造成的異常狀態(tài)不可能連續(xù)出現(xiàn)l+1次，因此上述觸發(fā)機(jī)制能夠很好地消除隨機(jī)噪聲對(duì)軸承退化檢測(cè)的影響。采用以上方法檢測(cè)到軸承退化后，在后續(xù)的狀態(tài)監(jiān)測(cè)過(guò)程中，根據(jù)軸承振動(dòng)信號(hào)預(yù)測(cè)軸承的剩余壽命，從而為軸承的維護(hù)維修提供參考。

2.2 退化模型

軸承的退化模型是用來(lái)描述軸承健康狀態(tài)與其工作時(shí)間關(guān)系的數(shù)學(xué)模型。由于RMS描述的是信號(hào)強(qiáng)度，能夠反應(yīng)軸承整體損傷，對(duì)于軸承的故障發(fā)展程度較為敏感。因此，當(dāng)使用2.1節(jié)中方法檢測(cè)到軸承退化后，采用RMS來(lái)表征軸承的故障發(fā)展程度。RMS是一個(gè)有量綱參數(shù)，會(huì)受到負(fù)載、安裝等工況的影響，因而容易引發(fā)錯(cuò)誤的預(yù)測(cè)結(jié)果。為了降低這一影響，利用相對(duì)RMS值(Relative Root Mean Square，RRMS)作為健康指標(biāo)，用于軸承剩余壽命預(yù)測(cè)。RRMS的具體定義為

(10)

式中：XRRMS(t)為信號(hào)采集時(shí)刻t的RRMS值；XRMS(t)為該時(shí)刻的RMS絕對(duì)數(shù)值；XRMS(FPTl)為軸承剩余壽命預(yù)測(cè)起始點(diǎn)的RMS絕對(duì)數(shù)值。

通常，機(jī)械設(shè)備的退化過(guò)程是一個(gè)隨機(jī)過(guò)程。機(jī)械設(shè)備隨機(jī)退化過(guò)程的不確定性主要來(lái)源于四個(gè)方面：①臨時(shí)不確定度；②單元-單元不確定度；③非線性不確定度；④測(cè)量不確定度。Si的研究中提出了一種包含前三種不確定度的自適應(yīng)退化模型。結(jié)合該模型與軸承的退化特點(diǎn)，建立軸承的退化模型為

X(t)=a+λtb+σB(t)

(11)

式中：X(t)為時(shí)刻t的RRMS值；a為系統(tǒng)的初始狀態(tài)，通常等于0；λ為時(shí)變參數(shù)，代表系統(tǒng)的退化速度；b為系統(tǒng)的非線性度；σB(t)為服從正太分布N(0，σ2)的布朗運(yùn)動(dòng)。

在軸承的狀態(tài)監(jiān)測(cè)過(guò)程中，很難準(zhǔn)確地對(duì)軸承的運(yùn)行狀態(tài)進(jìn)行測(cè)量，因而采用振動(dòng)信號(hào)這種間接測(cè)量方式反應(yīng)其健康狀態(tài)，但間接測(cè)量結(jié)果與其真實(shí)狀態(tài)間存在一定的誤差。同時(shí)，信號(hào)處理與特征提取過(guò)程中都會(huì)產(chǎn)生一定的誤差。因此，得到的特征與軸承的真實(shí)狀態(tài)之間存在誤差，這些誤差統(tǒng)一作為測(cè)量噪聲。因而建立軸承的測(cè)量方程為

Y(t)=X(t)+ν

(12)

式中：ν為軸承的測(cè)量噪聲，服從正太分布N(0,γ2)。

由于工作環(huán)境和材料參數(shù)等各方面的原因，每個(gè)軸承的參數(shù)λ，b取值均不相同。為了準(zhǔn)確地預(yù)測(cè)軸承剩余壽命，必須準(zhǔn)確地估計(jì)各個(gè)參數(shù)的值。根據(jù)軸承的退化模型，可以得到其狀態(tài)轉(zhuǎn)移方程和測(cè)量方程

(13)

式中：η=σ(B(tt)-B(tt-1))服從正太分布N(0,σ2Δt)，Δt=tt-tt-1。

2.3 剩余壽命預(yù)測(cè)

軸承退化模型確定后，便可以通過(guò)1.2節(jié)所述的UPF算法步驟利用測(cè)量值對(duì)模型參數(shù)進(jìn)行更新，得到測(cè)量時(shí)刻t的模型參數(shù)，并估計(jì)此時(shí)的退化狀態(tài)。利用退化模型傳遞現(xiàn)有分布，可以預(yù)測(cè)t+p時(shí)刻的狀態(tài)。在PF算法中，一種比較簡(jiǎn)單的預(yù)測(cè)方法就是計(jì)算每個(gè)粒子在t+p時(shí)刻的狀態(tài)，從而得到系統(tǒng)的狀態(tài)。

(14)

假設(shè)t時(shí)刻得到的粒子能夠準(zhǔn)確地表示t時(shí)刻的系統(tǒng)狀態(tài)，則可以通過(guò)遞推的方式，最終得到t+p時(shí)刻的狀態(tài)。遞推公式為

(15)

3 案例分析

為了闡述提出方法的有效性，采用實(shí)際試驗(yàn)中得到的全壽命周期振動(dòng)信號(hào)對(duì)提出方法進(jìn)行驗(yàn)證。并將提出方法與傳統(tǒng)PF方法進(jìn)行對(duì)比，說(shuō)明提出方法的優(yōu)越性。

3.1 試驗(yàn)臺(tái)與試驗(yàn)數(shù)據(jù)

本文所采用的數(shù)據(jù)由FEMTO-ST機(jī)構(gòu)提供，該數(shù)據(jù)曾用于2012年壽命預(yù)測(cè)與健康管理數(shù)據(jù)挑戰(zhàn)[16]。采集該數(shù)據(jù)的試驗(yàn)平臺(tái)如圖2所示。通過(guò)加速性能退化技術(shù)，該試驗(yàn)系統(tǒng)能夠在較短時(shí)間內(nèi)收集軸承的全壽命周期數(shù)據(jù)。試驗(yàn)中，軸承振動(dòng)信號(hào)的采樣頻率為25.6 kHz，數(shù)據(jù)長(zhǎng)度為2 560，每隔10 s采集一次。試驗(yàn)中，分別采集了橫向和垂向的振動(dòng)信號(hào)。但試驗(yàn)過(guò)程中施加的力為徑向力，因此，橫向振動(dòng)信號(hào)更能反映軸承的工作狀態(tài)。因此，本文選用橫向振動(dòng)信號(hào)進(jìn)行分析，從橫向振動(dòng)信號(hào)中提取特征表征軸承的健康狀況。

圖2 加速壽命試驗(yàn)臺(tái)Fig.2 Overview of the experimental platform

圖3為一個(gè)典型的軸承全壽命周期振動(dòng)信號(hào)。由圖可以看出，軸承經(jīng)歷了一段較長(zhǎng)的正常工作時(shí)間，然后開(kāi)始退化，一旦開(kāi)始退化后，軸承的性能迅速退化直至失效。

圖3 軸承振動(dòng)信號(hào)Fig.3 The vibration signal of a bearing

3.2 軸承剩余壽命預(yù)測(cè)

本文提出的方法需要先判斷其是否進(jìn)入退化階段，即找出第一個(gè)剩余壽命預(yù)測(cè)點(diǎn)，然后再對(duì)軸承進(jìn)行剩余壽命預(yù)測(cè)。根據(jù)第2部分所述方法，提取軸承振動(dòng)信號(hào)的峭度值與RMS，如圖4與圖5所示。首先，根據(jù)2.1節(jié)所述方法，利用峭度值判斷軸承退化起始點(diǎn)。由圖4可知，當(dāng)軸承正常工作時(shí)，峭度值比較穩(wěn)定，同時(shí)由于隨機(jī)噪聲的存在，這個(gè)過(guò)程中也會(huì)出現(xiàn)異常值。采用2.1節(jié)中所述觸發(fā)機(jī)制消除隨機(jī)噪聲的影響，結(jié)果表明，當(dāng)t=10 910 s時(shí)，軸承開(kāi)始退化，如圖4中的豎線所示。因此，選擇t=10 910 s為軸承的第一個(gè)剩余壽命預(yù)測(cè)點(diǎn)。相應(yīng)地，圖5中的豎線表示軸承的退化起始點(diǎn)，此時(shí)軸承的RMS值變化較大，此后RMS值一直隨著退化程度的加深而增大。

圖4 軸承振動(dòng)信號(hào)的峭度值Fig.4 Kurtosis results

圖5 軸承振動(dòng)信號(hào)的RMS值Fig.5 RMS results

檢測(cè)到軸承開(kāi)始退化后，計(jì)算每一個(gè)信號(hào)采集時(shí)刻的RRMS值，利用2.2節(jié)與2.3節(jié)所述方法完成軸承的剩余壽命預(yù)測(cè)。圖6為退化模型參數(shù)的更新過(guò)程，圖中橫坐標(biāo)0 s表示第一個(gè)剩余壽命預(yù)測(cè)點(diǎn)，后續(xù)橫坐標(biāo)值均表示當(dāng)前時(shí)刻距離第一個(gè)剩余壽命預(yù)測(cè)點(diǎn)的時(shí)間，圖7、圖8、圖9以及圖10中橫坐標(biāo)的意義均與圖6一樣。由圖可知，剛開(kāi)始進(jìn)行參數(shù)估計(jì)的時(shí)候，由于涉及的測(cè)量信息較少，參數(shù)估計(jì)結(jié)果存在較大誤差。隨著測(cè)量數(shù)據(jù)的增多，參數(shù)λ與參數(shù)b均逐漸收斂于真實(shí)值。因此，采用本文提出方法能夠有效地利用測(cè)量數(shù)據(jù)估計(jì)退化模型參數(shù)。

為了說(shuō)明本文提出方法的有效性，將PF方法作為對(duì)比。分別利用兩種方法對(duì)軸承退化階段的狀態(tài)進(jìn)行追蹤，得到的結(jié)果如圖7所示。由圖可知，PF方法和UPF方法都能夠準(zhǔn)確地追蹤RRMS值，在少數(shù)點(diǎn)處，UPF的估計(jì)結(jié)果比PF的估計(jì)結(jié)果更為準(zhǔn)確。

圖6 退化模型參數(shù)更新過(guò)程Fig.6 Update process of model parameters

圖7 軸承退化過(guò)程中RRMS估計(jì)值Fig.7 RRMS estimations of bearing degradation process

為了驗(yàn)證相對(duì)于PF方法，UPF方法能夠有效解決粒子退化問(wèn)題，采用有效粒子數(shù)來(lái)衡量狀態(tài)更新過(guò)程中粒子的退化程度[17]。狀態(tài)估計(jì)過(guò)程中，兩種方法所選擇的粒子數(shù)均為1 000，整個(gè)過(guò)程中，兩種方法的有效粒子數(shù)如圖8所示。由圖可以看出，迭代開(kāi)始后，兩種方法的粒子均迅速退化，但整個(gè)過(guò)程中，UPF方法的有效粒子數(shù)始終大于PF方法的有效粒子數(shù)。

為了驗(yàn)證UPF方法在壽命預(yù)測(cè)中的作用，選用第10個(gè)時(shí)間點(diǎn)和第35個(gè)時(shí)間點(diǎn)的預(yù)測(cè)結(jié)果來(lái)說(shuō)明UPF方法相對(duì)于PF方法的優(yōu)越性。圖9為已知10組數(shù)據(jù)時(shí)的預(yù)測(cè)結(jié)果，由圖可以看出與PF方法相比，UPF的預(yù)測(cè)誤差較小。已知10組數(shù)據(jù)時(shí)，使用PF方法預(yù)測(cè)到RRMS值超過(guò)失效閾值的時(shí)間為580 s，而真實(shí)的失效時(shí)間為490 s，因此，其剩余壽命預(yù)測(cè)誤差為90 s。而當(dāng)使用UPF方法時(shí)，預(yù)測(cè)到RRMS值超過(guò)失效閾值的時(shí)間為520 s，其預(yù)測(cè)誤差為30 s。圖10為已知35組數(shù)據(jù)時(shí)的預(yù)測(cè)結(jié)果，由圖10可知，隨著測(cè)量值的增加，兩種預(yù)測(cè)方法的準(zhǔn)確度都隨之增加。但UPF方法的預(yù)測(cè)結(jié)果依然比PF方法好，此時(shí)，PF預(yù)測(cè)方法預(yù)測(cè)軸承的失效時(shí)間為520 s，預(yù)測(cè)誤差為30 s，而UPF方法預(yù)測(cè)軸承失效時(shí)間為480 s，預(yù)測(cè)誤差為10 s。UPF方法的預(yù)測(cè)誤差低于PF方法的原因可以解釋如下：在UPF方法中，每個(gè)粒子均采用UKF方法進(jìn)行了更新，因而不會(huì)與實(shí)際狀態(tài)產(chǎn)生很大的偏差。在后續(xù)權(quán)值更新過(guò)程中，由于其與測(cè)量值相差不會(huì)過(guò)大，因而避免了大部分粒子權(quán)值過(guò)小發(fā)生退化，能夠更好地估計(jì)系統(tǒng)狀態(tài)。因此，UPF方法能夠有效解決PF方法中的粒子退化問(wèn)題，提高軸承剩余壽命預(yù)測(cè)的準(zhǔn)確度。

圖8 有效粒子數(shù)對(duì)比Fig.8 Comparison of the effective sample sizes

圖9 10組數(shù)據(jù)時(shí)的預(yù)測(cè)結(jié)果Fig.9 Prediction results with 10 measurements

圖10 35組數(shù)據(jù)時(shí)的預(yù)測(cè)結(jié)果Fig.10 Prediction results with 35 measurements

4 結(jié) 論

PF算法是一種廣泛用于解決目標(biāo)追蹤問(wèn)題的方法。近年來(lái)，PF技術(shù)被引入剩余壽命預(yù)測(cè)領(lǐng)域。但隨著迭代次數(shù)的增加，PF算法會(huì)出現(xiàn)粒子退化問(wèn)題。UPF算法結(jié)合PF和UKF，能夠有效解決這一問(wèn)題。本文提出了一種基于UPF的軸承剩余壽命預(yù)測(cè)方法，結(jié)合隨機(jī)過(guò)程模型與UPF方法，對(duì)退化過(guò)程中的軸承狀態(tài)進(jìn)行追蹤，并完成其剩余壽命預(yù)測(cè)。結(jié)果表明，UPF方法能夠有效彌補(bǔ)PF方法粒子退化這一缺陷，與傳統(tǒng)PF方法相比，本文提出方法能夠更加準(zhǔn)確地預(yù)測(cè)軸承的剩余壽命。