范洪敏

（中國直升機(jī)設(shè)計(jì)研究所，江西 景德鎮(zhèn) 333001）

隨著高可靠性、長壽命的產(chǎn)品的增多，產(chǎn)品壽命數(shù)據(jù)的獲取越來越難。因此，基于性能退化數(shù)據(jù)的可靠性評估方法的研究越來越深入。例如，趙新超等研究了多元性能的可靠性評估問題，周真等使用Bayes方法融合先驗(yàn)信息提高了模型精度。王浩偉等和訾佼佼等分別使用退化量分布方法對電連接器和電主軸進(jìn)行了可靠性評估。

以上研究中，在模型擬合時(shí)，都將初值直接假設(shè)為0或某個(gè)非0常數(shù)，沒有考慮初值的隨機(jī)性和分散性。然而，在工程應(yīng)用中，由于生產(chǎn)工藝和流程不可能完美，不同的產(chǎn)品個(gè)體在初始時(shí)刻的退化性能也不可能相同?；谶@一問題，X. Xiao等和Fan Hongmin等進(jìn)行了基于Wiener過程模型的隨機(jī)初值退化數(shù)據(jù)研究。在這一基礎(chǔ)上，本文以金屬裂紋尺寸退化數(shù)據(jù)為例，使用退化軌跡擬合和退化量分布兩種方法對金屬裂紋尺寸的退化數(shù)據(jù)進(jìn)行可靠性評估。

1 基于退化軌跡擬合的金屬裂紋尺寸退化數(shù)據(jù)可靠性評估方法

退化軌跡擬合法是一種常用的可靠性評估方法。主要步驟為：首先，擬合每個(gè)樣本的退化軌跡；然后，計(jì)算出樣本偽壽命；最后，對偽壽命數(shù)據(jù)進(jìn)行可靠性分析。

文獻(xiàn)[7]和文獻(xiàn)[8]中提供了某金屬產(chǎn)品的裂紋尺寸增長數(shù)據(jù)，每組數(shù)據(jù)的測量間隔為百萬次運(yùn)行周期，失效閾值l=1.3。表1為部分?jǐn)?shù)據(jù)。

表1 部分裂紋尺寸增長數(shù)據(jù)

1．1 模型建立

（1）線性模型

其中，參數(shù)b表示退化速率，A表示退化量初值，將其假定為一個(gè)服從正態(tài)分布N(μ,σ2)的隨機(jī)變量。因此，退化量的密度函數(shù)為：

（2）指數(shù)模型

（3）冪模型

式中，參數(shù)a表示退化速率，B表示初值項(xiàng)，服從對數(shù)正態(tài)分布。因此，退化量的密度函數(shù)為：

1．2 模型選擇

使用赤池信息量準(zhǔn)則（AIC）選擇合適的退化軌跡模型，公式為：

式中，k表示未知參數(shù)個(gè)數(shù)，lnL表示對數(shù)似然函數(shù)值。部分樣品的各模型AIC值見表2。

表2 各模型AIC值

通過對比，對于所有樣品,指數(shù)模型的AIC值都最小，選擇指數(shù)模型擬合退化軌跡。

1．3 計(jì)算偽壽命

顯然，失效時(shí)間T服從正態(tài)分布，其均值即為平均壽命，是標(biāo)志產(chǎn)品平均工作時(shí)長的量。將均值作為偽壽命近似值，得到部分產(chǎn)品偽壽命結(jié)果見表3。

表3 部分偽壽命計(jì)算結(jié)果

1．4 建立壽命模型

分別使用正態(tài)，對數(shù)正態(tài)和威布爾分布擬合偽壽命，并估計(jì)參數(shù)，計(jì)算AIC值。結(jié)果見表4。

表4 壽命模型參數(shù)估計(jì)結(jié)果

根據(jù)表4，選擇對數(shù)正態(tài)分布模型來擬合偽壽命?？煽慷葹椋?/p>

通過計(jì)算，可靠性預(yù)測曲線如圖1所示。

圖1 可靠性預(yù)測曲線

2 基于退化量分布的金屬裂紋尺寸退化數(shù)據(jù)可靠性評估方法

2．1 建立退化量分布模型

分別使用對數(shù)正態(tài)，威布爾和兩參數(shù)指數(shù)分布擬合各時(shí)刻退化數(shù)據(jù)，并且估計(jì)參數(shù)。部分結(jié)果見表5。

根據(jù)表5，選擇對數(shù)正態(tài)分布模型作為退化量分布模型。

2．2 參數(shù)軌跡擬合

對數(shù)正態(tài)模型的分布函數(shù)為：

從表5和分布函數(shù)中可以看出，參數(shù)σ的取值普遍較小。因此，參數(shù)μ即使出現(xiàn)極小的變化，也會對分布函數(shù)產(chǎn)生很大的影響?？紤]到這點(diǎn)，在對參數(shù)σ進(jìn)行擬合時(shí)，增加一項(xiàng)服從正態(tài)分布的隨機(jī)變量作為初值項(xiàng)。

（1）擬合參數(shù)μ

圖2為參數(shù)μ與運(yùn)行周期的關(guān)系圖。從圖中的軌跡趨勢可以看出，使用線性模型進(jìn)行參數(shù)軌跡擬合較為合理。

圖2 參數(shù)μ與運(yùn)行周期關(guān)系圖

模型形式為：

模型中參數(shù)b表示退化速率，參數(shù)a為截距。使用極大似然估計(jì)法擬合模型參數(shù)，可以得到參數(shù)μ的表達(dá)式：。

（2）擬合參數(shù)σ

分別使用冪模型、指數(shù)模型和二次多項(xiàng)式模型擬合參數(shù)σ的軌跡。為了計(jì)算方便，先暫時(shí)不考慮隨機(jī)初值項(xiàng)，三種模型表達(dá)式為：

表5 對數(shù)正態(tài)分布模型參數(shù)估計(jì)結(jié)果

為了選擇合適的模型，將誤差平方和作為反映模型擬合優(yōu)度的量。誤差平方和為各時(shí)刻擬合數(shù)據(jù)和實(shí)際數(shù)據(jù)之間差值的平方和，計(jì)算公式如下：

經(jīng)計(jì)算，指數(shù)模型、冪模型和二次多項(xiàng)式模型的誤差平方和分別為 1.02×10-3、3.41×10-4和3.41×10-5。因此，使用二次多項(xiàng)式模型參數(shù)σ。考慮到初值的隨機(jī)性，在模型中增加一個(gè)隨機(jī)項(xiàng)，即：

2．3 建立可靠度模型

在參數(shù)δ取置信水平為95%時(shí)，可靠性預(yù)測曲線如圖3所示。

圖3 可靠性預(yù)測曲線

圖4為兩種方法結(jié)果比較圖。可以看出，兩種分析方法的可靠性預(yù)測結(jié)果吻合度良好，具有一致性?？梢哉J(rèn)為結(jié)果是可信的。

圖4 可靠度預(yù)測值比較

3 結(jié)語

基于退化數(shù)據(jù)的可靠性評估方法可分為退化軌跡擬合法和退化量分布法。本文中，將隨機(jī)初值引入這兩種方法。具體做法為，在模型擬合時(shí)，增加一個(gè)隨機(jī)變量作為初值項(xiàng)。由于引入了隨機(jī)因素，在擬合偽壽命以及退化量分布參數(shù)時(shí)，大大提高了計(jì)算的精確度和可靠性。可靠度預(yù)測結(jié)果如圖4所示，兩種方法計(jì)算結(jié)果互相對照，可以證明這一方法的合理性。