崔 策， 石新發(fā)， 賀石中， 李秋秋， 康鑫碩

(廣州機(jī)械科學(xué)研究院有限公司 設(shè)備潤(rùn)滑與檢測(cè)研究所, 廣東 廣州 510700)

引言

近年來(lái)，隨著風(fēng)力發(fā)電行業(yè)的快速發(fā)展，風(fēng)能已成為世界主要可再生能源之一[1]。但由于風(fēng)機(jī)安裝、運(yùn)行、維護(hù)和零部件的生產(chǎn)制造過(guò)程中的管理經(jīng)驗(yàn)不足，導(dǎo)致我國(guó)風(fēng)能利用率較低，風(fēng)機(jī)故障事件頻繁發(fā)生。為提升風(fēng)能轉(zhuǎn)換效率及降低風(fēng)機(jī)故障率， 科研人員在槳葉控制、功率轉(zhuǎn)換及維護(hù)方案上做了大量研究[2]。液力變槳系統(tǒng)以其響應(yīng)快、扭矩大的特點(diǎn)廣泛應(yīng)用于風(fēng)機(jī)葉片的槳距控制機(jī)構(gòu)，因此正確開(kāi)展液力變槳系統(tǒng)可靠性分析對(duì)于風(fēng)機(jī)運(yùn)行安全維護(hù)有著十分積極的作用[3]。

目前，在機(jī)電設(shè)備的可靠性分析過(guò)程中，通常采用基于物理模型或數(shù)據(jù)驅(qū)動(dòng)的分析方法進(jìn)行評(píng)估。物理模型分析法首先充分分析對(duì)象的結(jié)構(gòu)特點(diǎn)及物化特性，采集各項(xiàng)實(shí)驗(yàn)參數(shù)，在對(duì)失效機(jī)理充分分析后以評(píng)估對(duì)象可靠性或剩余壽命，但由于其成本較高，該類(lèi)方法通常應(yīng)用于高精尖領(lǐng)域設(shè)備如航空發(fā)動(dòng)機(jī)、飛機(jī)起落架等。近年來(lái)，數(shù)據(jù)驅(qū)動(dòng)方法由于其成本低、不依賴(lài)設(shè)備失效機(jī)理等特點(diǎn)得到廣泛應(yīng)用，包括基于統(tǒng)計(jì)模型的方法、基于可靠性函數(shù)的方法及基于機(jī)器學(xué)習(xí)的方法[4]。其中統(tǒng)計(jì)模型Weibull分布是常見(jiàn)的設(shè)備可靠性評(píng)估方法，由于三參數(shù)Weibull分布在擬合數(shù)據(jù)方面具有較強(qiáng)適應(yīng)性，幾種常見(jiàn)的分布如指數(shù)分布、正態(tài)分布等都可看成特定參數(shù)的Weibull分布，目前已有大量學(xué)者針對(duì)此類(lèi)方法開(kāi)展可靠性研究。文獻(xiàn)[5]使用最小二乘法對(duì)二參數(shù)的Weibull分布進(jìn)行三參數(shù)估計(jì)，并用航空設(shè)備實(shí)例數(shù)據(jù)驗(yàn)證方法的可行性。但對(duì)于常用的三參數(shù)的 Weibull分布沒(méi)有進(jìn)行相應(yīng)的研究。文獻(xiàn)[6]等對(duì)幾種常用的Weibull分布參數(shù)估計(jì)方法進(jìn)行了說(shuō)明和比較，相比之下，極大似然估計(jì)方法具有比較優(yōu)越的估計(jì)精度。文獻(xiàn)[7]提出一種針對(duì)競(jìng)爭(zhēng)失效模式下的三參數(shù)威布爾分布實(shí)驗(yàn)壽命統(tǒng)計(jì)分析，并采用極大似然模型進(jìn)行實(shí)例驗(yàn)證，結(jié)果表明了分析方法和模型的可行性，但計(jì)算過(guò)程較為復(fù)雜。本研究采用三參數(shù)Weibull分布擬合測(cè)試數(shù)據(jù)，對(duì)于三參數(shù)的Weibull分布建立極大似然方程求解形狀參數(shù)、尺度參數(shù)及位置參數(shù)，為解決計(jì)算過(guò)程復(fù)雜的問(wèn)題，采用遺傳算法迭代求解方程組。遺傳算法是一種模擬自然界選擇與遺傳的機(jī)理尋求最優(yōu)解的優(yōu)化方法。遺傳算法不要求函數(shù)連續(xù)，具有可擴(kuò)展性和潛在的并行性等優(yōu)點(diǎn)，己得到廣泛運(yùn)用。文獻(xiàn)[8]采用遺傳算法和參數(shù)優(yōu)化的支持向量機(jī)建模，同時(shí)進(jìn)行故障特征提取和模型訓(xùn)練，用該模型研究7種典型的制冷機(jī)組故障，取得良好的識(shí)別結(jié)果。

1 Weibull分布建模

在工業(yè)設(shè)備的可靠性及壽命研究中，Weibull分布尤其適用于具有累計(jì)失效分布形式的設(shè)備，因此被廣泛應(yīng)用于各種設(shè)備或元件的壽命或可靠性評(píng)估[9]。通過(guò)形狀參數(shù)、尺度參數(shù)和位置參數(shù)的變化，Weibull分布可以近似描述各種故障率分布規(guī)律，具有較好的泛化性[10]。建立Weibull分布?jí)勖翱煽啃栽u(píng)估模型的難點(diǎn)在于精確的估計(jì)模型參數(shù)。目前，常用的參數(shù)估計(jì)法有圖估計(jì)法、最小二乘法、極大似然估計(jì)法以及回歸分析法等[11]。

使用三參數(shù)的Weibull分布對(duì)液力變槳系統(tǒng)進(jìn)行可靠性評(píng)估，可靠度函數(shù)為：

(1)

失效分布函數(shù)為：

(2)

故障率密度函數(shù)為：

(3)

故障率函數(shù)為：

(4)

其中，β，η，γ分別表示形狀參數(shù)、尺度參數(shù)和位置參數(shù)，當(dāng)形狀參數(shù)β<1時(shí)，λ(t)是t的遞減函數(shù)；當(dāng)β=1時(shí)，故障率λ(t)為固定值；當(dāng)β>1時(shí)，λ(t)是遞增函數(shù)，且此時(shí)故障率曲線拐點(diǎn)為：(e1/β-1)/e1/β。

經(jīng)初步統(tǒng)計(jì)發(fā)現(xiàn)，在大型風(fēng)機(jī)液力變槳系統(tǒng)的運(yùn)行周期內(nèi)，風(fēng)機(jī)液力變槳系統(tǒng)的運(yùn)行初期故障率較高，其故障形式和系統(tǒng)配置、維護(hù)等人為因素相關(guān)，故障率函數(shù)隨時(shí)間衰減；而設(shè)備運(yùn)行中期失效率通常恒定且明顯低于初期水平，失效事件主要以隨機(jī)離散事件為主，如系統(tǒng)密封失效導(dǎo)致進(jìn)水、污染等。而運(yùn)行末期具有典型的損耗特性，如潤(rùn)滑油氧化及磨損累積導(dǎo)致潤(rùn)滑表面失效等。綜上所述，風(fēng)機(jī)液力變槳系統(tǒng)的故障率變化特征與不同參數(shù)Weibull的故障率函數(shù)變化一致，因此適用于液力變槳系統(tǒng)的可靠性研究過(guò)程。

2 參數(shù)估計(jì)

二參數(shù)的Weibull分布建模忽略位置參數(shù)γ的影響，默認(rèn)失效從設(shè)備運(yùn)行開(kāi)始之前便一直存在。當(dāng)采用二參數(shù)Weibull分布進(jìn)行可靠性評(píng)估時(shí)，最小二乘法是常見(jiàn)的參數(shù)估計(jì)方法。最小二乘法主要是通過(guò)最小化理論值與觀測(cè)值之差的平方和，使得擬合對(duì)象無(wú)限接近目標(biāo)。采用最小二乘法進(jìn)行參數(shù)估計(jì)，首先對(duì)可靠性函數(shù)R(t)即式(1)等式兩邊分別取兩次對(duì)數(shù)，可得：

ln[-lnR(t)]=βlnt-βlnη

(5)

因此，將上式簡(jiǎn)化成y=kx+b的形式，則y=ln[-lnR(t)],k=β,b=-βlnη,x=lnt,通過(guò)計(jì)算，得出最小二乘估計(jì)為：

(6)

(7)

三參數(shù)Weibull分布由于加入了故障開(kāi)始時(shí)間(位置參數(shù)γ)，因此相比二參數(shù)Weibull可以更好地模擬系統(tǒng)失效特性。三參數(shù)Weibull分布的參數(shù)估計(jì)通常采用極大似然法求解，極大似然原理是利用已知樣本的結(jié)果信息，反推最有可能(最大概率)導(dǎo)致樣本結(jié)果的模型參數(shù)值，即“模型已知，但參數(shù)未知”， 極大似然估計(jì)是一種有效且精度較高的算法,但其求解效率卻一直受限于復(fù)雜的超越方程及龐大的數(shù)據(jù)本身[12]。似然函數(shù)如式(8)，其中θ為待估算參數(shù)，x表示隨機(jī)變量取值。

(8)

(9)

因此結(jié)合式(9)，三參數(shù)Weibull分布的概率密度函數(shù)極大似然函數(shù)為：

(10)

對(duì)上式取對(duì)數(shù)得：

l(t;β,η,γ)=lnL(t;β,η,γ)=nlnβ-nβlnη+

(11)

由于參數(shù)極大似然估計(jì)結(jié)果是似然函數(shù)的最大值，當(dāng)似然函數(shù)對(duì)3個(gè)參數(shù)(β,η,γ)連續(xù)可微，且最大值在定義域區(qū)間時(shí)，即對(duì)上式的3個(gè)參數(shù)分別求偏導(dǎo)并令其等于0，求解聯(lián)立方程組即為極大似然估計(jì)量的結(jié)果：

(12)

(13)

(14)

通過(guò)式(12)～式(14)聯(lián)立似然方程組，發(fā)現(xiàn)采用常規(guī)解析方法無(wú)法逐一求解。由于遺傳算法不要求函數(shù)連續(xù)，且具有拓展性強(qiáng)、魯棒性好等優(yōu)點(diǎn)，在參數(shù)求解、組合優(yōu)化等方面已得到廣泛應(yīng)用。遺傳算法借鑒自然界生物種群遺傳進(jìn)化機(jī)理，通過(guò)個(gè)體適應(yīng)能力的判定，及變異、交叉的遺傳方式，在迭代過(guò)程中盡量擴(kuò)大優(yōu)秀個(gè)體及其后代占比，以尋找適應(yīng)能力最高的后代，即函數(shù)最優(yōu)解[13]。本研究將遺傳算法與三參數(shù)Weibull分布的極大似然函數(shù)方程組結(jié)合起來(lái)，逐步求解Weibull的參數(shù)。

針對(duì)傳統(tǒng)遺傳算法的基因交叉及變異率固定不變，可能導(dǎo)致迭代過(guò)程陷入局部最優(yōu)解的現(xiàn)象，提出動(dòng)態(tài)修正基因的交叉概率及變異概率Pm，Pc。對(duì)于適應(yīng)度較大的個(gè)體，降低其變異的概率，同時(shí)增加交叉的概率，以增加優(yōu)秀基因的遺傳后代占比；若個(gè)體的適應(yīng)度低于整體平均適應(yīng)度時(shí)，則增大其變異概率，降低交叉概率。當(dāng)個(gè)體的適應(yīng)值最大，其變異概率為0，并在迭代過(guò)程中直接復(fù)制保留，個(gè)體的交叉及變異如下：

令Pci=1-Pmi，其中fmax，fmin表示個(gè)體最大及最小適應(yīng)度。由于遺傳算法不能直接處理問(wèn)題空間的參數(shù)，因此必須通過(guò)編碼將求解的目標(biāo)轉(zhuǎn)化成具有遺傳特性的染色體或者個(gè)體。在編譯過(guò)程中，個(gè)體可采用二進(jìn)制、十進(jìn)制以及浮點(diǎn)數(shù)等方式編碼。浮點(diǎn)數(shù)在編碼方式上需要遵循一定的標(biāo)準(zhǔn)，計(jì)算過(guò)程相對(duì)繁瑣；十進(jìn)制編碼同樣位數(shù)上能存儲(chǔ)更多的信息，但是變異過(guò)程中誤差變化較大；二進(jìn)制編碼在后續(xù)的交叉和變異等操作中的誤差精度控制相比十進(jìn)制更好，因此本研究采用二進(jìn)制編碼。

通過(guò)觀察式(12)～式(14)發(fā)現(xiàn)， 由于位置參數(shù)對(duì)失效概率密度函數(shù)只起到調(diào)節(jié)位置的作用，求解時(shí)首先固定位置參數(shù)， 則聯(lián)立方程可簡(jiǎn)化為尺度參數(shù)及形狀參數(shù)的二元方程，經(jīng)進(jìn)一步簡(jiǎn)化后可得：

(15)

(16)

3 案例分析

表1記錄了某大型風(fēng)機(jī)液力變槳系統(tǒng)的15組時(shí)序失效數(shù)據(jù)，測(cè)試結(jié)果表示從t=0到監(jiān)測(cè)結(jié)果中首次出現(xiàn)異常故障的時(shí)間。

表1 某大型風(fēng)機(jī)液力變槳系統(tǒng)失效數(shù)據(jù)

傳統(tǒng)遺傳算法的停止條件中包含遺傳代數(shù)達(dá)到最大或適應(yīng)值達(dá)到迭代次數(shù)內(nèi)的最佳，設(shè)置初始交叉概率Pc=0.7，變異概率Pm=0.3，群體規(guī)模popsize=500，最大迭代次數(shù)為500。初始種群中個(gè)體以隨機(jī)方式產(chǎn)生，個(gè)體的求解區(qū)間為[0,5]，且當(dāng)?shù)螖?shù)不小于50，且最優(yōu)個(gè)體適應(yīng)度恒定不變時(shí)停止迭代。

由圖1可知，在迭代過(guò)程中，種群個(gè)體平均適應(yīng)度逐步上升，表明在動(dòng)態(tài)修正交叉及變異概率的策略下，群體的平均適應(yīng)度增大，而曲線跌宕現(xiàn)象可能由于變異的隨機(jī)性導(dǎo)致。當(dāng)?shù)恋?0代，出現(xiàn)適應(yīng)度函數(shù)值最大的最優(yōu)個(gè)體，且在50代次內(nèi)，并未出現(xiàn)更優(yōu)個(gè)體，因此求得適應(yīng)度最高個(gè)體即方程最優(yōu)解β=2.02，帶入式(16)和式(14)，同時(shí)計(jì)算出η=845.67，γ=14.52，擬合度檢驗(yàn)AD統(tǒng)計(jì)量為0.447。采用最小二乘法計(jì)算得出的二參Weibull分布的數(shù)據(jù)β=1.67，η=738.63，擬合度檢驗(yàn)AD統(tǒng)計(jì)量為0.494，圖2分別是二、三參數(shù)Weibull分布下的概率圖線性擬合，圖中數(shù)據(jù)點(diǎn)均收斂于95%置信區(qū)間內(nèi)，并分布在直線周?chē)叶叩臄M合優(yōu)度P值均大于0.05，對(duì)比發(fā)現(xiàn)三參數(shù)的AD統(tǒng)計(jì)量更小，表明三參數(shù)Weibull分布函數(shù)擬合程度更優(yōu)，更能表現(xiàn)該風(fēng)機(jī)液力變槳系統(tǒng)的失效特性。

圖1 遺傳算法求解目標(biāo)函數(shù)

圖2 Weibull分布下概率圖線性擬合

圖3 失效率函數(shù)擬合對(duì)比

因此得出大型風(fēng)力發(fā)電機(jī)液力變槳系統(tǒng)的可靠性函數(shù)為：

失效率密度函數(shù)為：

圖4分別是失效密度曲線和可靠性及失效概率曲線，同時(shí)計(jì)算出故障率曲線的拐點(diǎn)為：

(e1/2.02-1)/e1/2.02=0.39

4 結(jié)論

本研究簡(jiǎn)要分析了大型風(fēng)機(jī)液力變槳系統(tǒng)的失效特性，結(jié)合三參數(shù)Weibull分布函數(shù)建立時(shí)序失效模型，其中函數(shù)參數(shù)采用極大似然法進(jìn)行求解，建立了三參數(shù)極大似然方程組。針對(duì)似然方程求解困難的問(wèn)題，提出一種改進(jìn)遺傳算法進(jìn)行求解，在遺傳算法的個(gè)體迭代過(guò)程中動(dòng)態(tài)設(shè)置了交叉及變異概率，同時(shí)對(duì)最優(yōu)的個(gè)體進(jìn)行復(fù)制保留，加快了遺傳算法的搜尋過(guò)程，得到似然函數(shù)最優(yōu)解。

結(jié)合某大型風(fēng)機(jī)液力變槳系統(tǒng)的實(shí)際失效數(shù)據(jù)，分別采用二參數(shù)及三參數(shù)Weibull分布函數(shù)對(duì)其失效分布特征進(jìn)行擬合。 結(jié)果表明， 與二參數(shù)的Weibull分布函數(shù)擬合結(jié)果對(duì)比，三參數(shù)Weibull分布函數(shù)擬合優(yōu)度檢驗(yàn)結(jié)果更優(yōu)，且由于增加了位置參數(shù)的考量，三參數(shù)Weibull分布函數(shù)的擬合結(jié)果更符合該液力變槳系統(tǒng)的失效特征，其失效率曲線的初始階段更接近實(shí)際故障發(fā)生規(guī)律，能夠較好地評(píng)估其可靠性過(guò)程。