基于離散小波變換和變分自編碼器的滾動(dòng)軸承剩余使用壽命預(yù)測(cè)

孟祥龍，丁華，呂彥寶，施瑞

(1.太原理工大學(xué) 機(jī)械與運(yùn)載工程學(xué)院，太原 030024；2.煤礦綜采裝備山西省重點(diǎn)實(shí)驗(yàn)室，太原 030024)

滾動(dòng)軸承廣泛應(yīng)用于眾多機(jī)械設(shè)備中,作為關(guān)鍵的承載部件,長(zhǎng)期處于復(fù)雜、惡劣的工作環(huán)境中,具有故障率高,易損壞的特點(diǎn),是機(jī)械故障的主要源頭,因此軸承剩余使用壽命(Remaining Useful Life, RUL)預(yù)測(cè)對(duì)于機(jī)械設(shè)備的預(yù)測(cè)性維護(hù)有著重要意義[1]。精準(zhǔn)地預(yù)測(cè)軸承剩余使用壽命，可以讓機(jī)器以最小的維護(hù)成本來(lái)獲得最長(zhǎng)的安全運(yùn)行時(shí)間，進(jìn)而提升機(jī)器運(yùn)行過(guò)程中的安全性和經(jīng)濟(jì)性。

人工智能的發(fā)展使數(shù)據(jù)驅(qū)動(dòng)的軸承健康監(jiān)測(cè)維護(hù)[2]相比傳統(tǒng)方法更具優(yōu)越性，基于數(shù)據(jù)驅(qū)動(dòng)的軸承剩余壽命預(yù)測(cè)主要分2步：1)從采集的振動(dòng)加速度信號(hào)中提取軸承性能退化特征；2)構(gòu)建預(yù)測(cè)模型，利用機(jī)器學(xué)習(xí)或深度學(xué)習(xí)卓越的非線性映射能力尋找數(shù)據(jù)特征與剩余使用壽命之間的潛在關(guān)系[3]。特征提取是數(shù)據(jù)驅(qū)動(dòng)方法的重要步驟，其關(guān)鍵在于從監(jiān)測(cè)數(shù)據(jù)中挖掘機(jī)械設(shè)備運(yùn)行狀態(tài)信息，對(duì)剩余使用壽命預(yù)測(cè)精度有很大影響。傳統(tǒng)特征提取方法是從原始振動(dòng)信號(hào)中提取退化特征，如提取時(shí)域的均方根值、頻域的平均幅值以及時(shí)頻域的小波包能量譜等[4]；文獻(xiàn)[5]則利用卷積神經(jīng)網(wǎng)絡(luò)特有的局部卷積、權(quán)值共享等特性直接對(duì)原始信號(hào)進(jìn)行特征自學(xué)習(xí)，省去了專(zhuān)家先驗(yàn)知識(shí)，實(shí)現(xiàn)了端對(duì)端的剩余使用壽命預(yù)測(cè)。然而，原始信號(hào)中含有大量噪聲，軸承的退化特征極易被噪聲淹沒(méi)，影響剩余使用壽命預(yù)測(cè)精度。當(dāng)從不同分析域中提取特征指標(biāo)后，很難判斷其中哪些特征對(duì)軸承缺陷發(fā)展更加敏感[6]；并且這些特征仍是高維，需要通過(guò)降維融合剔除冗余信息，獲得能全面反應(yīng)軸承性能退化趨勢(shì)的綜合指標(biāo)。文獻(xiàn)[7]利用主成分分析從多維數(shù)據(jù)中獲取優(yōu)勢(shì)信息的特性來(lái)實(shí)現(xiàn)特征的融合降維，文獻(xiàn)[8]采用自編碼器構(gòu)建了低維性能退化指標(biāo)；然而，主成分分析處理線性問(wèn)題具有優(yōu)勢(shì)，卻難以充分挖掘軸承變化的有效信息；自編碼器是將高維數(shù)據(jù)映射為低維特定數(shù)值[9]而不是某種分布，得到的低維特征受到約束，退化趨勢(shì)較單一。

針對(duì)以上問(wèn)題，本文提出一種基于離散小波變換(Discrete Wavelet Transform，DWT)和變分自編碼器(Variational Auto-Encoder，VAE)的滾動(dòng)軸承剩余使用壽命預(yù)測(cè)方法。通過(guò)離散小波變換對(duì)軸承振動(dòng)信號(hào)進(jìn)行逐級(jí)分解以抑制信號(hào)中的環(huán)境噪聲和無(wú)關(guān)成分，從而選取有效信息成分多、噪聲成分少的分量作為后續(xù)算法處理的對(duì)象，實(shí)現(xiàn)降噪預(yù)處理的效果；同時(shí)，利用變分自編碼器強(qiáng)大的無(wú)監(jiān)督特征學(xué)習(xí)能力和抗干擾能力獲得更符合軸承退化趨勢(shì)的綜合特征；最后，通過(guò)分層抽樣劃分?jǐn)?shù)據(jù)集，并引入長(zhǎng)短時(shí)記憶(Long Short-Term Memory，LSTM)網(wǎng)絡(luò)實(shí)現(xiàn)滾動(dòng)軸承剩余使用壽命預(yù)測(cè)。

1 基于DWT和VAE的軸承壽命預(yù)測(cè)

1.1 離散小波變換

(1)

則稱(chēng)為一個(gè)基本小波或母小波，將母函數(shù)ψ(t)經(jīng)伸縮和平移后，可以得到一個(gè)小波序列。

對(duì)于連續(xù)的情況，小波序列為

(2)

連續(xù)小波變換為

(3)

式中：a為伸縮因子；b為平移因子。

對(duì)于離散的情況，小波序列為

(4)

通常取a0=2,b0=1，則 (4) 式變?yōu)?/p>

(5)

3層離散小波變換的原理如圖1所示，其中H(w)為低通濾波器，可以濾除曲線中的高頻部分，得到近似部分(低頻成分)；G(w)為高通濾波器，可以濾除曲線中的低頻部分，得到細(xì)節(jié)成分(高頻成分)；從而將信號(hào)的時(shí)頻結(jié)構(gòu)準(zhǔn)確無(wú)冗余地展現(xiàn)出來(lái)。

圖1 3層離散小波變換原理示意圖

1.2 變分自編碼器

變分自編碼器是一種基于自編碼器的無(wú)監(jiān)督學(xué)習(xí)方法，同時(shí)也是一種生成模型，通過(guò)編碼器獲得潛在空間特征的概率分布，融入隨機(jī)噪聲并以解碼器的形式生成新數(shù)據(jù)，其基本結(jié)構(gòu)如圖2所示。與自編碼器直接將輸入x映射到隱含向量z不同，變分自編碼器的編碼部分是將輸入x映射到隱變量概率分布的參數(shù)(均值、方差)，再加入隨機(jī)噪聲后采樣得到隱變量z，最后由解碼部分生成新的數(shù)據(jù)。

圖2 變分自編碼器基本結(jié)構(gòu)

變分自編碼器包含編碼器、解碼器和損失函數(shù)3部分。編碼器是一個(gè)神經(jīng)網(wǎng)絡(luò)，輸入為x，輸出為隱向量z，參數(shù)為θ，可表示為pθ(z∣x)；解碼器也是一個(gè)神經(jīng)網(wǎng)絡(luò)，輸入為隱向量z，輸出為重構(gòu)的x，參數(shù)為φ，可表示為qφ(x∣z)；損失函數(shù)為

L(θ,φ)=Ez～pθ(z∣x)[lgqφ(x∣z)]+

DKL[pθ(z∣x)p(z)]，

(6)

等式右側(cè)第1項(xiàng)為重構(gòu)損失，與其他自編碼器的損失函數(shù)一樣，目的是讓生成數(shù)據(jù)與原始數(shù)據(jù)盡可能相近；第2項(xiàng)為KL 散度(相對(duì)熵)，是正則項(xiàng)，用來(lái)衡量2種不同概率分布之間的偏離程度。

在變分自編碼器中，隱空間先驗(yàn)分布p(z)被設(shè)定為標(biāo)準(zhǔn)正態(tài)分布，正則項(xiàng)的存在是要讓后驗(yàn)分布pθ(z∣x)也接近先驗(yàn)分布。如果沒(méi)有正則項(xiàng)，模型為了減小重構(gòu)損失會(huì)不斷減小隨機(jī)性，也就是編碼器輸出的方差，沒(méi)有隨機(jī)性時(shí)，變分自編碼器無(wú)法生成各種數(shù)據(jù)。若編碼器輸出的z不服從標(biāo)準(zhǔn)正態(tài)分布，將會(huì)在損失函數(shù)中對(duì)編碼器施加懲罰。

變分自編碼器中，pθ(z∣x)確定一個(gè)高斯分布并從中采樣得到隱變量z，然而采樣操作的不可導(dǎo)使模型無(wú)法反向傳播。在實(shí)際模型中，采用重參數(shù)技巧解決后驗(yàn)分布采樣過(guò)程中梯度不可計(jì)算的問(wèn)題，即從均值μ和標(biāo)準(zhǔn)偏差σ的正態(tài)分布中采樣，等價(jià)于先從標(biāo)準(zhǔn)正態(tài)分布中采樣，然后再對(duì)其進(jìn)行變換，即

z=μ+σ⊙，

(7)

1.3 剩余使用壽命標(biāo)簽設(shè)定

在實(shí)際應(yīng)用中，軸承運(yùn)行前期振動(dòng)信號(hào)的振幅十分平穩(wěn)，可視為健康運(yùn)行階段。軸承運(yùn)行一段時(shí)間后，在某一時(shí)刻發(fā)生初始故障，標(biāo)志著軸承開(kāi)始進(jìn)入退化階段，軸承振動(dòng)信號(hào)振幅不再趨于平穩(wěn)而是劇烈振蕩且伴隨著軸承性能的急劇下降，直至軸承失效，如圖3a所示。由于軸承在健康階段并沒(méi)有出現(xiàn)一定的退化演變趨勢(shì)，使用圖3b中線性剩余使用壽命標(biāo)簽進(jìn)行剩余使用壽命預(yù)測(cè)時(shí)效果較差。因此，在軸承運(yùn)行健康階段，將剩余使用壽命設(shè)定為一個(gè)常數(shù)值較為合理[10]，直到出現(xiàn)初始故障，軸承開(kāi)始進(jìn)入退化階段，剩余使用壽命才開(kāi)始逐漸減小，可以采用分段函數(shù)表示[11]，如圖3b中虛線所示。

(a)原始振動(dòng)信號(hào)

1.4 模型介紹

長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)是循環(huán)神經(jīng)網(wǎng)絡(luò)的一個(gè)變體，可以有效地解決循環(huán)神經(jīng)網(wǎng)絡(luò)的梯度爆炸或消失問(wèn)題，其在循環(huán)神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上引入門(mén)控機(jī)制控制信息的累積速度，包括有選擇地加入新信息，有選擇地遺忘之前累積的信息。長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)的具體公式為

it=σ(Wxixt+Whiht-1+bi)，

(8)

ft=σ(Wxfxt+Whfht-1+bf)，

(9)

ot=σ(Wxoxt+Whoht-1+bo)，

(10)

ct=ft⊙ct-1+it⊙tanh(Wxcxt+Whcht-1=bc)，

(11)

ht=ot⊙tanh(ct)，

(12)

式中：it,ft,ot分別為輸入門(mén)、遺忘門(mén)和輸出門(mén)；c，h分別為記憶單元和外部狀態(tài)；σ為 Logistic函數(shù)，輸出區(qū)間為(0,1)；W為各神經(jīng)元的權(quán)重，其與輸入矢量xt和外部狀態(tài)ht-1直接點(diǎn)乘。

長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)的循環(huán)單元結(jié)構(gòu)如圖4所示，其計(jì)算過(guò)程為：1)利用上一時(shí)刻的外部狀態(tài)ht-1和當(dāng)前時(shí)刻的輸入xt，計(jì)算出3個(gè)門(mén)it,ft,ot以及候選狀態(tài)ct；2)結(jié)合遺忘門(mén)ft和輸入門(mén)it更新記憶單元ct；3)結(jié)合輸出門(mén)ot將內(nèi)部狀態(tài)的信息傳遞給外部狀態(tài)ht。通過(guò)長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)循環(huán)單元，整個(gè)網(wǎng)絡(luò)可以建立較長(zhǎng)距離的時(shí)序依賴(lài)關(guān)系。

圖4 長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)結(jié)構(gòu)示意圖

2 實(shí)例驗(yàn)證

2.1 數(shù)據(jù)來(lái)源

試驗(yàn)數(shù)據(jù)來(lái)自美國(guó)辛辛那提大學(xué)IMS中心軸承試驗(yàn)臺(tái)[12]，試驗(yàn)軸承為 ZA-2115型滾動(dòng)軸承，軸轉(zhuǎn)速為 2 000 r/min，由彈簧機(jī)構(gòu)施加 26.67 kN 的徑向力加速軸承退化，采用高靈敏度PCB-353B33加速度傳感器收集軸承振動(dòng)信號(hào)，采樣頻率為20.48 kHz，采樣間隔為10 min，采樣時(shí)間1 s，故每次采樣點(diǎn)數(shù)為20 480。選擇第2組試驗(yàn)中軸承1的全壽命周期振動(dòng)數(shù)據(jù)進(jìn)行分析，其在歷時(shí)9 840 min后發(fā)生疲勞失效(拆機(jī)發(fā)現(xiàn)外圈出現(xiàn)裂紋)，其振動(dòng)信號(hào)在后期產(chǎn)生劇烈變化，反映出軸承發(fā)生了急劇退化。

2.2 試驗(yàn)結(jié)果分析

首先需要確定軸承初始故障發(fā)生的時(shí)刻，并以此作為初始退化點(diǎn)來(lái)劃分軸承運(yùn)行的健康階段和退化階段。由于軸承早期故障特征難以檢測(cè)，致使初始退化點(diǎn)難以確定，采用文獻(xiàn)[13]中所提方法進(jìn)行初始退化點(diǎn)檢測(cè)。軸承發(fā)生故障時(shí)會(huì)在其頻譜中有所體現(xiàn)，可以利用頻譜中最大峰值頻率的改變找出初始退化發(fā)生的時(shí)刻。將每次采樣得到的振動(dòng)信號(hào)進(jìn)行快速傅里葉變換得到頻譜，把頻譜中幅值最大的頻率記錄下來(lái)，得到的散點(diǎn)圖如圖5所示：軸承運(yùn)行前期只存在1 000 Hz的頻率，隨著軸承持續(xù)運(yùn)行，第7 000 min以后開(kāi)始持續(xù)出現(xiàn)4 476 Hz的故障頻率，故將第7 000 min確定為軸承的初始退化時(shí)刻。

圖5 初始退化點(diǎn)檢測(cè)

分段函數(shù)將0～7 000 min視為健康狀態(tài)，歸一化后RUL標(biāo)簽為1；第7 000 min之后視為退化狀態(tài)，其標(biāo)簽值遞減直至為0，表征軸承損壞已達(dá)閾值，無(wú)法再正常運(yùn)行，如圖6所示。

圖6 歸一化分段函數(shù)RUL標(biāo)簽

在特征提取方面?？紤]到db小波基函數(shù)具有正交性、緊支撐性以及對(duì)不規(guī)則信號(hào)更敏感等特性，選擇db4小波作為基函數(shù)，將軸承原始振動(dòng)信號(hào)進(jìn)行3層離散小波變換處理，得到近似分量信號(hào)和細(xì)節(jié)分量信號(hào)并從中分別提取標(biāo)準(zhǔn)差、方根幅值、均方根、偏度、峭度以及波形因子特征，結(jié)果如圖7所示。

從圖7中剔除對(duì)軸承退化趨勢(shì)不敏感的偏度、峭度以及波形因子，選擇單調(diào)趨勢(shì)性較好的標(biāo)準(zhǔn)差、方根幅值和均方根作為敏感特征構(gòu)建高維特征集，如圖8所示：從近似分量中提取的均方根在退化階段的單調(diào)上升趨勢(shì)更符合指數(shù)型增長(zhǎng)，其明顯優(yōu)于細(xì)節(jié)分量中提取的方根幅值；而從局部放大圖中可以看出，近似分量的均方根在健康階段非常平穩(wěn)，而細(xì)節(jié)分量的方根幅值在健康階段已有輕微的起伏上升，表明細(xì)節(jié)分量特征對(duì)于軸承的損壞更加敏感，能提取到軸承早期微弱的退化信息；因此，單單使用均方根或方根幅值無(wú)法完整表達(dá)軸承整個(gè)壽命周期的退化趨勢(shì)，需要通過(guò)特征融合獲得反映軸承退化的綜合性能指標(biāo)。

圖7 離散小波變換處理后軸承振動(dòng)信號(hào)的特征指標(biāo)Fig.7 Feature indexes of bearing vibration signal after DWT processing

圖8 離散小波變換處理后軸承振動(dòng)信號(hào)的部分敏感特征

利用歸一化操作對(duì)敏感特征構(gòu)成的高維特征集進(jìn)行縮放并輸入變分自編碼器中進(jìn)行無(wú)監(jiān)督特征融合降維，將高維特征集降至一維后所得軸承的深層本質(zhì)特征如圖9所示，其中x軸代表時(shí)間點(diǎn)，y軸代表將高維特征集輸入到變分自編碼器后融合降維得到的一維深層特征值。由圖9可知：由

圖9 變分自編碼器降維得到的一維特征

離散小波變換和變分自編碼器提取出的深層特征具有良好的單調(diào)性，既能體現(xiàn)健康階段的輕微損傷(局部放大圖)，又能在退化階段表現(xiàn)出良好的指數(shù)型單調(diào)上升，能夠綜合表征軸承整個(gè)生命周期的衰退過(guò)程，可以將其作為基本特征輸入長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型進(jìn)行軸承剩余使用壽命預(yù)測(cè)。

試驗(yàn)軸承全壽命周期數(shù)據(jù)共有984個(gè)時(shí)間點(diǎn)，最后8個(gè)時(shí)刻的監(jiān)控?cái)?shù)據(jù)與實(shí)際趨勢(shì)不一致，誤差較大。為避免這一誤差的影響，將這部分?jǐn)?shù)據(jù)剔除。采用分層抽樣的方法將全壽命周期數(shù)據(jù)劃分為訓(xùn)練集和測(cè)試集[14]，以10個(gè)時(shí)間點(diǎn)為一層，以7:3的比例對(duì)每一層數(shù)據(jù)進(jìn)行隨機(jī)抽樣構(gòu)建訓(xùn)練集和測(cè)試集，訓(xùn)練集和測(cè)試集分別占全壽命周期數(shù)據(jù)的70%和30%，最后得到訓(xùn)練集與測(cè)試集樣本數(shù)的比為683：293。

長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型重要的超參數(shù)有時(shí)間步長(zhǎng)、隱藏層狀態(tài)數(shù)、批量大小和學(xué)習(xí)率，超參數(shù)的選擇會(huì)影響預(yù)測(cè)結(jié)果的精確性，因此采用粒子群優(yōu)化算法(Particle Swarm Optimization，PSO) 對(duì)超參數(shù)進(jìn)行設(shè)定，設(shè)置粒子群優(yōu)化算法搜索空間的上、下界分別為[20,128,100,0.01]，[2,8,10,0.000 1]，種群數(shù)量為20，使模型預(yù)測(cè)效果達(dá)到最優(yōu)的超參數(shù)為[10,70,32,0.005 3]。另外，由于Adam作為自適應(yīng)優(yōu)化算法能夠動(dòng)態(tài)更新學(xué)習(xí)率，因此選擇Adam優(yōu)化器進(jìn)行梯度優(yōu)化。

長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型在測(cè)試集上的預(yù)測(cè)結(jié)果如圖10所示，預(yù)測(cè)值與RUL標(biāo)簽有很高的重合度。為評(píng)估剩余使用壽命預(yù)測(cè)的不確定性,采用文獻(xiàn)[15]方法對(duì)性能退化階段進(jìn)行區(qū)間估計(jì),在預(yù)測(cè)值附近設(shè)置 95%置信水平的置信區(qū)間。雖然有些預(yù)測(cè)值偏離RUL標(biāo)簽曲線(局部放大圖)，但絕大部分預(yù)測(cè)值均落在置信區(qū)間內(nèi)，說(shuō)明預(yù)測(cè)值趨勢(shì)與RUL標(biāo)簽曲線趨勢(shì)一致。

圖10 試驗(yàn)軸承在測(cè)試集上的RUL預(yù)測(cè)結(jié)果

2.3 對(duì)比試驗(yàn)

為驗(yàn)證方法的有效性，采用3種特征提取方法：1)直接從原始振動(dòng)信號(hào)中提取時(shí)頻特征；2)對(duì)原始信號(hào)進(jìn)行小波軟硬閾值折衷法去噪后[16]提取時(shí)頻特征；3)對(duì)原始信號(hào)進(jìn)行離散小波變換，對(duì)信號(hào)子分量提取敏感特征。融合降維則使用2種方法，即堆疊自編碼器(Stacked Auto-Encoder，SAE)和變分自編碼器。將3種特征提取方法與2種融合降維方法相互組合，共形成6種方法進(jìn)行對(duì)比，其得到的一維特征如圖11所示。

圖11 不同特征提取和降維融合方法提取的一維深層特征Fig.11 One-dimensional deep features extracted by different eature extraction and dimension reduction fusion methods

由圖11可知：1)原始信號(hào)深層特征在健康階段起伏上升現(xiàn)象明顯，而經(jīng)過(guò)小波軟硬閾值折衷法去噪和離散小波變換處理后所提取深層特征在健康階段十分平穩(wěn)，起伏上升得到大幅抑制，且離散小波變換所提取特征在退化階段呈指數(shù)型上升，更符合軸承退化趨勢(shì)，對(duì)軸承振動(dòng)信號(hào)特征提取的效果更好；2)以小波去噪為例，堆疊自編碼器提取的深層特征在退化階段表現(xiàn)出先輕微下降再緩慢上升，即“愈合”現(xiàn)象，而變分自編碼器提取的深層特征在退化階段先趨于穩(wěn)定再快速上升，退化趨勢(shì)明顯更優(yōu)，說(shuō)明相對(duì)于堆疊自編碼器，變分自編碼器的融合降維效果更佳，可以有效抑制退化階段的“愈合”現(xiàn)象[17]，獲得單調(diào)性更優(yōu)的特征指標(biāo)。

從定量分析的角度，進(jìn)一步對(duì)一維深層特征的趨勢(shì)性、單調(diào)性和魯棒性進(jìn)行分析[18]，其計(jì)算結(jié)果見(jiàn)表1，表明離散小波變換和變分自編碼器算法所提取深層特征的綜合性能更優(yōu)。

表1 不同方法所提取深層特征的定量評(píng)價(jià)指標(biāo)

采用分層抽樣對(duì)每種方法得到的深層特征劃分訓(xùn)練集與測(cè)試集，通過(guò)訓(xùn)練集來(lái)訓(xùn)練長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)模型，并使用測(cè)試集驗(yàn)證預(yù)測(cè)效果，結(jié)果如圖12所示。

圖12 不同方法所構(gòu)建模型對(duì)試驗(yàn)軸承在測(cè)試集上的剩余使用壽命預(yù)測(cè)結(jié)果Fig.12 RUL prediction results of test bearing on test set by models constructed by different methods

選擇能反映測(cè)量精度的均方根誤差(Root Mean Square Error，RMSE)和能反映預(yù)測(cè)誤差實(shí)際情況的平均絕對(duì)誤差(Mean Absolute Error，MAE)作為模型的評(píng)價(jià)指標(biāo)[19](表2)，評(píng)價(jià)指標(biāo)對(duì)應(yīng)的計(jì)算公式為

表2 不同方法所構(gòu)建預(yù)測(cè)模型的評(píng)價(jià)指標(biāo)

(13)

(14)

分析圖12和表2可以看出，基于離散小波變換和變分自編碼器的滾動(dòng)軸承剩余使用壽命預(yù)測(cè)方法的誤差更小，精度更高。

2.4 不同初始退化點(diǎn)的預(yù)測(cè)結(jié)果

上述分析中，軸承初始退化點(diǎn)被確定為第7 000 min并以此劃分健康階段和退化階段，然而文獻(xiàn)[20]指出軸承在第5 330 min已發(fā)生輕微退化，到第7 000 min已是中度退化狀態(tài)。因此，本文嘗試分析不同初始退化點(diǎn)的選擇是否會(huì)對(duì)所提方法產(chǎn)生影響，故將第5 330 min設(shè)定為初始退化點(diǎn)，則健康階段為0～5 330 min，退化階段為5 330～9 840 min。相應(yīng)的分段函數(shù)如圖13所示，其余處理與上述方法一致，得到的預(yù)測(cè)結(jié)果如圖14所示。

圖13 不同初始退化點(diǎn)時(shí)的剩余使用壽命分段函數(shù)

由圖14可知：對(duì)于從原始振動(dòng)信號(hào)中提取時(shí)頻特征再融合降維得到的一維深層特征，其在5 330～7 000 min區(qū)間起伏上升明顯，故在選取不同初始退化點(diǎn)時(shí)，測(cè)試集(1 600～2 100 min，即5 330×30%～ 7 000×30%)的預(yù)測(cè)效果略優(yōu)于小波去噪和離散小波變換方法，但是第2 100 min之后的預(yù)測(cè)效果卻遠(yuǎn)不如離散小波變換和變分自編碼器相結(jié)合的方法，從整體上看，離散小波變換與變分自編碼器相結(jié)合的方法依舊能達(dá)到最優(yōu)的預(yù)測(cè)效果。

圖14 不同初始退化點(diǎn)時(shí)各種模型在測(cè)試集上的RUL預(yù)測(cè)結(jié)果Fig.14 RUL prediction results of various models on test set at different initial degradation points

3 結(jié)束語(yǔ)

針對(duì)滾動(dòng)軸承工作環(huán)境中噪聲干擾使有效特征難提取的問(wèn)題，提出一種基于離散小波變換和變分自編碼器的軸承剩余使用壽命預(yù)測(cè)方法，試驗(yàn)結(jié)果表明：

1)利用離散小波變換分解信號(hào)后對(duì)其子分量提取特征，能夠有效抑制軸承退化中的“愈合”現(xiàn)象，使提取的特征更具單調(diào)趨勢(shì)性,準(zhǔn)確地反映軸承生命周期的退化過(guò)程。

2)利用變分自編碼器強(qiáng)大的無(wú)監(jiān)督學(xué)習(xí)能力和抗干擾能力，能夠獲得綜合表征軸承退化趨勢(shì)的性能指標(biāo)，從而提高軸承剩余使用壽命預(yù)測(cè)精度。

3)選擇不同初始退化點(diǎn)時(shí)，離散小波變換和變分自編碼器相結(jié)合的特征提取方法依然能精準(zhǔn)的預(yù)測(cè)剩余使用壽命。