基于退化量分布時(shí)序分析的產(chǎn)品壽命預(yù)測(cè)

王 立 李曉陽 姜同敏

(北京航空航天大學(xué) 可靠性與系統(tǒng)工程學(xué)院,北京 100191)

基于退化量分布時(shí)序分析的產(chǎn)品壽命預(yù)測(cè)

王 立 李曉陽 姜同敏

(北京航空航天大學(xué) 可靠性與系統(tǒng)工程學(xué)院,北京 100191)

針對(duì)現(xiàn)有性能退化預(yù)測(cè)方法采用單調(diào)回歸函數(shù)對(duì)產(chǎn)品退化隨機(jī)過程描述不夠合理的問題,從產(chǎn)品在各時(shí)刻退化量分布角度出發(fā),分析退化量分布參數(shù)估計(jì)的非平穩(wěn)時(shí)序類型,并對(duì)各參數(shù)估計(jì)分別采用相應(yīng)類型的非平穩(wěn)時(shí)序分析方法建模,進(jìn)而給出一種與實(shí)際退化隨機(jī)過程更為相符的基于退化量分布的產(chǎn)品壽命及可靠度預(yù)測(cè)方法.通過對(duì)某電子產(chǎn)品的退化試驗(yàn),采用所提出的壽命預(yù)測(cè)方法進(jìn)行了該電子產(chǎn)品壽命及可靠度預(yù)測(cè).結(jié)果表明該方法相比傳統(tǒng)方法更符合實(shí)際.

性能退化;壽命預(yù)測(cè);退化量分布;非平穩(wěn)時(shí)序分析

越來越多長(zhǎng)壽命高可靠性產(chǎn)品出現(xiàn),使產(chǎn)品壽命與可靠性評(píng)估更加困難.基于產(chǎn)品性能退化信息預(yù)測(cè)產(chǎn)品壽命及可靠度成為一種有效途徑.

目前性能退化預(yù)測(cè)主要有兩種思路:①將性能退化量隨時(shí)間變化的隨機(jī)過程各樣本函數(shù)稱為退化軌跡,基于退化軌跡預(yù)測(cè).②將性能退化量在不同時(shí)刻所服從分布的參數(shù)看作隨機(jī)變量,基于退化量分布進(jìn)行預(yù)測(cè).

現(xiàn)有性能退化預(yù)測(cè)研究中,文獻(xiàn)[1]假設(shè)退化量為正態(tài)分布進(jìn)行可靠性預(yù)計(jì);文獻(xiàn)[2]給出性能退化量分布參數(shù)的極大似然估計(jì)方法;文獻(xiàn)[3]采用退化軌跡函數(shù)預(yù)測(cè)產(chǎn)品性能退化;文獻(xiàn)[4]采用線性模型來描述退化量分布參數(shù)和退化軌跡.這些文獻(xiàn)均將產(chǎn)品退化軌跡或退化量分布參數(shù)假設(shè)為單調(diào)回歸函數(shù),進(jìn)行性能退化預(yù)測(cè).然而,實(shí)際工程中由于受到環(huán)境等因素影響,性能退化量存在隨機(jī)性及周期性變化,且在這種變化下退化量分布參數(shù)往往具有非平穩(wěn)隨機(jī)性.現(xiàn)有方法對(duì)產(chǎn)品退化隨機(jī)過程描述不夠合理.

非平穩(wěn)時(shí)序包括方差平穩(wěn)時(shí)序和相關(guān)系數(shù)平穩(wěn)時(shí)序等.非平穩(wěn)時(shí)序分析利用現(xiàn)代統(tǒng)計(jì)學(xué)和信息處理技術(shù),能充分挖掘非平穩(wěn)隨機(jī)時(shí)序的自相關(guān)性,刻畫時(shí)序隨機(jī)性波動(dòng)規(guī)律,是一種適于描述退化隨機(jī)過程的方法.文獻(xiàn)[5]給出基于退化軌跡時(shí)序分析的性能退化產(chǎn)品壽命預(yù)測(cè)方法;文獻(xiàn)[6]給出基于退化量分布的性能退化預(yù)測(cè)時(shí)序分析法,但將退化量分布參數(shù)時(shí)序籠統(tǒng)視為方差平穩(wěn)隨機(jī)時(shí)序.由于退化量分布不同參數(shù)屬于不同非平穩(wěn)時(shí)序類型,將所有參數(shù)時(shí)序籠統(tǒng)視為方差平穩(wěn)隨機(jī)時(shí)序的假設(shè)過于簡(jiǎn)單,與實(shí)際情況不完全相符.因此本文從退化量分布角度出發(fā),對(duì)各退化量分布參數(shù),分別采用相應(yīng)非平穩(wěn)時(shí)序分析方法,進(jìn)行性能退化產(chǎn)品壽命預(yù)測(cè)方法研究.

1 退化量分布數(shù)字特征分析

產(chǎn)品退化量分布反映了退化隨機(jī)過程一維分布,其基本數(shù)字特征為均值和方差.若各樣本退化軌跡為單調(diào)回歸函數(shù),退化隨機(jī)過程如圖 1.

可知,退化量分布均值和方差為單調(diào)回歸函數(shù).

然而工程中,由于產(chǎn)品受到各種環(huán)境因素影響,實(shí)際的退化隨機(jī)過程如圖 2.

圖 2 實(shí)際退化隨機(jī)過程

可知,當(dāng)各樣本退化軌跡為隨機(jī)時(shí)序時(shí),退化量分布均值和方差也可能為隨機(jī)時(shí)序.當(dāng)隨機(jī)時(shí)序具有自相關(guān)性時(shí),則時(shí)序取值不僅與當(dāng)前時(shí)刻有關(guān),且與其歷史時(shí)刻取值有關(guān).若仍采用單調(diào)回歸函數(shù)描述退化量分布均值和方差時(shí)序,則反映不出退化過程隨機(jī)波動(dòng)的自相關(guān)性,從而影響退化預(yù)測(cè)結(jié)果.因此,應(yīng)對(duì)退化量分布均值和方差分別采用時(shí)序分析方法描述.

2 退化量分布參數(shù)估計(jì)

本文提出所研究的性能退化過程假設(shè):

假設(shè) 1 產(chǎn)品的性能退化過程總體趨勢(shì)具有單調(diào)性,即性能退化總體趨勢(shì)不可逆.

假設(shè) 2 退化過程中,所有產(chǎn)品采樣時(shí)刻相等.

假設(shè) 3 隨著時(shí)間的變化,退化量分布的類型不變,僅參數(shù)變化.

設(shè)共有 n個(gè)產(chǎn)品樣本,以 yt表示產(chǎn)品在 t時(shí)刻的性能退化量(或性能退化量的某種非線性變換),以 yti表示 yt的第 i個(gè)產(chǎn)品樣本.當(dāng) yt服從某位置-尺度分布時(shí),任意時(shí)刻的位置參數(shù)與尺度參數(shù)即 yt的數(shù)字特征,分別記為 μt和,它們可確定 yt在 t時(shí)刻分布情況.例如:

1)當(dāng)性能退化量服從正態(tài)分布時(shí),yt表示性能退化量,其分布參數(shù)為 μt和

2)當(dāng)性能退化量服從對(duì)數(shù)正態(tài)分布時(shí),yt表示性能退化量的對(duì)數(shù),yt服從正態(tài)分布,其分布參數(shù)為 μt和

3)當(dāng)性能退化量服從形狀參數(shù)為 θt、尺度參數(shù)為 ηt的威布爾分布時(shí),yt表示性能退化量的對(duì)數(shù),yt服從極值分布,其分布參數(shù)為 μt=lnηt和μt和的估計(jì)值可由 yt的樣本均值和樣本方差得到.由于和是按時(shí)間順序變化的隨機(jī)變量,因此是時(shí)間序列,并且是具有一定退化趨勢(shì)的非平穩(wěn)時(shí)間序列.

3 退化量分布參數(shù)估計(jì)平穩(wěn)性分析

不同類型的非平穩(wěn)時(shí)序建模方法迥異,因此在對(duì)退化量分布參數(shù)估計(jì)采用時(shí)序分析方法建模前,必須先進(jìn)行平穩(wěn)性分析,以確定退化量分布參數(shù)估計(jì)的非平穩(wěn)時(shí)序類型.

1)產(chǎn)品樣本退化軌跡時(shí)序非平穩(wěn)性分析.

實(shí)際工程中,通常認(rèn)為產(chǎn)品各樣本在 m時(shí)間段內(nèi)退化軌跡 yti(t=1,2,…,m)為相互獨(dú)立的方差平穩(wěn)時(shí)序.由時(shí)序分解原理,yti可分為確定性時(shí)序 dti和平穩(wěn)隨機(jī)時(shí)序 rti的疊加.

2)樣本均值時(shí)序非平穩(wěn)性分析.

3)樣本方差時(shí)序非平穩(wěn)性分析.s2

t時(shí)序表達(dá)式為

式(3)中第一項(xiàng)為確定性時(shí)序.

其中,g(t)為關(guān)于 t的非零函數(shù);bi為常數(shù).式(3)第二項(xiàng)可表示為

根據(jù)文獻(xiàn)[7]對(duì)相關(guān)系數(shù)平穩(wěn)隨機(jī)過程定義:

定義 1 設(shè){xt,t∈ T}是隨機(jī)過程,如果

1){xt,t∈T}是二階矩過程;

2)對(duì)任意的 t,t+τ∈T,有

則{xt,t∈T}是相關(guān)系數(shù)平穩(wěn)隨機(jī)過程.

另外再根據(jù)文獻(xiàn)[8]中的定理:

定理 1 設(shè) xt為相關(guān)系數(shù)平穩(wěn)隨機(jī)過程,則)為平穩(wěn)隨機(jī)過程,其中D(xt)=σ2I2(t),σ為常數(shù).

由此,本文給出推論如下:

推論 1 任一非零確定性時(shí)序與平穩(wěn)時(shí)序之積為相關(guān)系數(shù)平穩(wěn)時(shí)序.

證明 設(shè) zt為任一平穩(wěn)時(shí)序,I(t)為任一關(guān)于 t的非零函數(shù),即非零確定性時(shí)序,則

由此可知 I(t)zt為相關(guān)系數(shù)平穩(wěn)時(shí)序. 證畢

因此,當(dāng) (dti-d-t)滿足式 (4)時(shí) ,式 (3)第二項(xiàng)為相關(guān)系數(shù)平穩(wěn)時(shí)序.

式(3)第三項(xiàng)中,當(dāng) rti為正態(tài)過程時(shí),(rtir-t)為平穩(wěn)時(shí)序且為正態(tài)過程.為常數(shù),并有

因此,本文僅考慮式(3)前兩項(xiàng),將 s2t時(shí)序視為確定性時(shí)序與相關(guān)系數(shù)平穩(wěn)時(shí)序之和.

4 退化量分布參數(shù)估計(jì)時(shí)序建模

產(chǎn)品退化量分布參數(shù)估計(jì)時(shí)序模型是根據(jù)產(chǎn)品樣本的退化軌跡時(shí)序模型推導(dǎo)而得,因此本文首先給出退化軌跡時(shí)序的建模方法.

1)產(chǎn)品樣本退化軌跡時(shí)序建模.

在產(chǎn)品性能退化過程中,產(chǎn)品除自身退化特性外,往往還受到環(huán)境及設(shè)備因素影響,包括周期性和隨機(jī)性影響,因此每個(gè)樣本退化量時(shí)序 yti的確定性部分 dti還可進(jìn)一步分解為單調(diào)趨勢(shì)項(xiàng) fti和周期項(xiàng) cti的疊加

fti可采用線性回歸函數(shù)或可轉(zhuǎn)化為線性的單調(diào)非線性回歸函數(shù)描述為

其中,bi為待定系數(shù);f0i為 fti初值;g(t)為單調(diào)非線性函數(shù),當(dāng) g(t)=t時(shí),式(7)為線性回歸函數(shù).

cti可采用適用于挖掘數(shù)據(jù)潛在周期性規(guī)律的潛周期模型描述,由于同一批受試產(chǎn)品樣本所受到的環(huán)境及設(shè)備影響相同,因此不同樣本的周期項(xiàng) cti理論上也相同,則有

其中,q為角頻率個(gè)數(shù);λj為角頻率;aj為幅值.

rti可采用傳統(tǒng)平穩(wěn)時(shí)序分析方法中工程應(yīng)用廣泛、建模簡(jiǎn)單且適于預(yù)測(cè)的自回歸模型描述為

其中,pi為自回歸模型階數(shù);ηji為自回歸系數(shù);εti為白噪聲.

2)樣本均值時(shí)序模型.

y-t具有與樣本退化軌跡時(shí)序相同表達(dá)形式:

即式(4)滿足,式(3)第二項(xiàng)為相關(guān)系數(shù)平穩(wěn)時(shí)序.

由式(11)可知,式(3)第一項(xiàng)為單調(diào)趨勢(shì)項(xiàng),以 fst表示趨勢(shì)項(xiàng),xt表示相關(guān)系數(shù)平穩(wěn)項(xiàng),有

由式(3)和式(11),fst模型為

其中,bs為待定系數(shù).

相關(guān)系數(shù)平穩(wěn)項(xiàng) xt的模型為

其中,fxt表示相關(guān)系數(shù)平穩(wěn)項(xiàng) xt的確定性部分,稱為相關(guān)系數(shù)平穩(wěn)趨勢(shì)項(xiàng),是關(guān)于 t的非零函數(shù);rxt表示相關(guān)系數(shù)平穩(wěn)項(xiàng) xt的隨機(jī)性部分,稱為相關(guān)系數(shù)平穩(wěn)隨機(jī)項(xiàng),是平穩(wěn)時(shí)序;σr為 rxt的標(biāo)準(zhǔn)差;σxt為 xt的標(biāo)準(zhǔn)差.

由文獻(xiàn)[8]可知,σxt可由 xt的絕對(duì)值序列|xt|的均值模型得到,即

可知 E|xt|即 fxt,由式(3)和式(11),fxt可采用單調(diào)回歸函數(shù)模型描述為

其中,bx為待定系數(shù).rxt采用自回歸模型描述為

其中,px,ηxj,εxt為 rxt的自回歸 模型階數(shù) 、自回歸系數(shù)、白噪聲.

5 基于退化量分布的壽命預(yù)測(cè)

1)退化量分布參數(shù)估計(jì)時(shí)序預(yù)測(cè).

根據(jù)時(shí)序模型最佳(最小均方誤差)預(yù)測(cè)原理,由式(10),y-t的向前 l步最佳預(yù)測(cè)值為

2)壽命預(yù)測(cè).

通過給定某失效閾值 D,可得到 yt在任意時(shí)刻到達(dá) D的概率,即產(chǎn)品可靠度 Rt.由第 2節(jié)有:

①當(dāng)性能退化量服從正態(tài)分布或?qū)?shù)正態(tài)分布時(shí),yt服從正態(tài)分布,若 yt隨 t單調(diào)上升,產(chǎn)品可靠度為

若 yt隨 t單調(diào)下降,產(chǎn)品可靠度為

②當(dāng)性能退化量服從威布爾分布時(shí),yt服從極值分布,若 yt隨 t單調(diào)上升,產(chǎn)品可靠度為

若 yt隨 t單調(diào)下降,產(chǎn)品可靠度為

此時(shí),產(chǎn)品壽命定義為產(chǎn)品性能穿越失效閾值的概率為 Rt時(shí),所對(duì)應(yīng)的時(shí)間 t,而并非產(chǎn)品性能第一次穿越失效閾值的時(shí)間.

6 實(shí) 例

對(duì)同一批次 13個(gè)某電子產(chǎn)品進(jìn)行退化試驗(yàn),如圖 3.對(duì)產(chǎn)品施加額定電壓,以串聯(lián)電阻兩端電壓值作為觀測(cè)數(shù)據(jù).采樣間距為 1min.試驗(yàn)共進(jìn)行了17000min.

圖 3 試驗(yàn)設(shè)備示意圖

由于噪聲干擾對(duì)非平穩(wěn)時(shí)序分析影響很大,在對(duì)試驗(yàn)觀測(cè)數(shù)據(jù)進(jìn)行處理前,應(yīng)先對(duì)其進(jìn)行濾波處理,以去掉噪聲干擾的影響.

經(jīng)過濾波后的試驗(yàn)觀測(cè)數(shù)據(jù)再經(jīng)初值化預(yù)處理后如圖 4.

圖 4 初值化后的試驗(yàn)觀測(cè)數(shù)據(jù)

采用皮爾遜 χ2擬合優(yōu)度檢驗(yàn)方法對(duì)各時(shí)刻對(duì)應(yīng)的預(yù)處理后數(shù)據(jù)分別進(jìn)行退化量分布假設(shè)檢驗(yàn).根據(jù)產(chǎn)品個(gè)數(shù),本文取 6個(gè)數(shù)據(jù)量值寬度相等的區(qū)間,以保證每個(gè)區(qū)間包含 2個(gè)以上產(chǎn)品樣本,則 χ2理論值為 7.815.若 χ2檢驗(yàn)值小于理論值,則接受數(shù)據(jù)分布假設(shè).對(duì)該數(shù)據(jù)所有時(shí)刻退化量分布類型的 χ2檢驗(yàn)平均值如表 1.

可知,威布爾分布檢驗(yàn)不通過,取 χ2檢驗(yàn)平均值最小的對(duì)數(shù)正態(tài)分布為退化量分布.

取經(jīng)過預(yù)處理后數(shù)據(jù)的對(duì)數(shù)作為退化數(shù)據(jù)進(jìn)行本方法的分析,如圖 5,黑色平滑曲線為采用冪函數(shù)對(duì)各產(chǎn)品樣本退化軌跡擬合的趨勢(shì)項(xiàng).

表 1 退化量分布 χ2檢驗(yàn)平均值

圖 5 退化軌跡時(shí)序及其趨勢(shì)項(xiàng)

對(duì)退化軌跡時(shí)序進(jìn)行趨勢(shì)項(xiàng)和周期項(xiàng)建模,減去這兩項(xiàng),得到退化軌跡時(shí)序隨機(jī)項(xiàng)如圖 6.

圖 6 退化軌跡時(shí)序隨機(jī)項(xiàng)

對(duì)退化軌跡時(shí)序隨機(jī)項(xiàng)采用輪次檢驗(yàn)法進(jìn)行平穩(wěn)性檢驗(yàn).本文將隨機(jī)項(xiàng)按時(shí)間分為 20等份,對(duì)各份的均方值進(jìn)行輪次計(jì)算,在 0.05顯著性水平下,若輪次數(shù)在 64至 125間時(shí),則該數(shù)據(jù)平穩(wěn),各產(chǎn)品結(jié)果如表 2.

表 2 退化軌跡時(shí)序隨機(jī)項(xiàng)平穩(wěn)性檢驗(yàn)

可知全部產(chǎn)品退化軌跡時(shí)序隨機(jī)項(xiàng)平穩(wěn).則該批產(chǎn)品退化軌跡為方差平穩(wěn)時(shí)序.對(duì)隨機(jī)項(xiàng)采用皮爾遜 χ2擬合優(yōu)度正態(tài)分布檢驗(yàn)結(jié)果如表 3.

表 3 退化軌跡時(shí)序隨機(jī)項(xiàng)正態(tài)分布檢驗(yàn)

可知全部產(chǎn)品退化軌跡時(shí)序隨機(jī)項(xiàng)符合正態(tài)分布,則式(3)第三項(xiàng)為平穩(wěn)時(shí)序.

分別對(duì)退化數(shù)據(jù)的樣本均值時(shí)序和樣本方差時(shí)序進(jìn)行建模.樣本均值時(shí)序的趨勢(shì)項(xiàng)模型為

樣本方差時(shí)序的趨勢(shì)項(xiàng)模型為

對(duì)減去趨勢(shì)項(xiàng)后的樣本方差殘差數(shù)據(jù)進(jìn)行平穩(wěn)性檢驗(yàn)結(jié)果為 57,檢驗(yàn)未通過,說明樣本方差時(shí)序并非方差平穩(wěn)時(shí)序.如圖 7.

如圖 8,樣本方差殘差數(shù)據(jù)絕對(duì)值的均值模型為

圖 7 樣本方差時(shí)序殘差數(shù)據(jù)

圖 8 樣本方差殘差數(shù)據(jù)絕對(duì)值及其均值

該殘差數(shù)據(jù)除以其絕對(duì)值均值后,再進(jìn)行平穩(wěn)性檢驗(yàn)結(jié)果為 105,檢驗(yàn)通過,說明樣本方差可視為相關(guān)系數(shù)平穩(wěn)時(shí)序與確定性時(shí)序之和.如圖9.

采用第 5節(jié)方法,將退化量分布樣本均值及樣本方差時(shí)序預(yù)測(cè)至 30000min,如圖 10、圖 11.

根據(jù)該電子產(chǎn)品的實(shí)際失效情況的經(jīng)驗(yàn),取產(chǎn)品失效閾值為初始值的 93%,得到可靠度與壽命關(guān)系預(yù)測(cè)結(jié)果.如圖 12.

圖 9 樣本方差相關(guān)系數(shù)平穩(wěn)隨機(jī)項(xiàng)

圖 10 樣本均值時(shí)序及其預(yù)測(cè)

圖 11 樣本方差時(shí)序及其預(yù)測(cè)

圖 12 可靠度與壽命關(guān)系預(yù)測(cè)

可見,采用本文基于退化量分布非平穩(wěn)時(shí)序分析得到的可靠度與壽命關(guān)系曲線相比單調(diào)回歸函數(shù)分析方法能反映出退化隨機(jī)過程的波動(dòng)性規(guī)律.該批電子產(chǎn)品中位壽命預(yù)測(cè)為 24580min.

7 結(jié) 論

本文的研究分析和實(shí)例應(yīng)用表明,對(duì)于實(shí)際工程中大多數(shù)產(chǎn)品,其各樣本退化軌跡為非平穩(wěn)隨機(jī)時(shí)序,則其退化量分布的數(shù)字特征往往也為非平穩(wěn)隨機(jī)時(shí)序.因而本文從退化量分布角度出發(fā),提出樣本均值和樣本方差時(shí)序的非平穩(wěn)時(shí)序分析方法,從而得到反映產(chǎn)品退化隨機(jī)過程波動(dòng)性規(guī)律的產(chǎn)品可靠度與壽命關(guān)系預(yù)測(cè).

實(shí)際上,這種可靠度隨壽命時(shí)間的波動(dòng)性是由于產(chǎn)品可靠度被定義為產(chǎn)品在各時(shí)刻退化量穿越失效閾值的概率,在退化量分布數(shù)字特征隨時(shí)間并非單調(diào)性變化的實(shí)際情況下,產(chǎn)品可靠度也因此不隨時(shí)間單調(diào)下降,而是在單調(diào)下降的總體趨勢(shì)下,還具有一定波動(dòng)性.這種波動(dòng)性既包含與當(dāng)前時(shí)刻相關(guān)的周期性波動(dòng),也包含與歷史時(shí)刻相關(guān)的隨機(jī)性波動(dòng).

由于考慮了產(chǎn)品退化過程中由環(huán)境因素引起的周期性和隨機(jī)性波動(dòng),與現(xiàn)有僅考慮產(chǎn)品各樣本退化軌跡單調(diào)性退化趨勢(shì)的方法相比,本文方法描述的產(chǎn)品退化隨機(jī)過程與實(shí)際情況更相符.

需指出,本文假設(shè)產(chǎn)品的失效閾值為一固定值,然而實(shí)際工程中,產(chǎn)品不同樣本的失效閾值往往是隨機(jī)性的.由于產(chǎn)品的失效閾值對(duì)產(chǎn)品壽命預(yù)測(cè)結(jié)果影響很大,考慮隨機(jī)失效閾值的基于退化量分布產(chǎn)品壽命預(yù)測(cè)方法值得進(jìn)一步研究.

此外,對(duì)于某些產(chǎn)品由于制造費(fèi)用昂貴,只能得到小樣本退化數(shù)據(jù)的情況,應(yīng)采用本文第 4節(jié)中產(chǎn)品樣本退化軌跡時(shí)序建模方法,即基于至少一個(gè)產(chǎn)品樣本的退化軌跡進(jìn)行產(chǎn)品退化過程描述和壽命預(yù)測(cè).對(duì)于某些產(chǎn)品由于退化過程過緩,只能在一定時(shí)間內(nèi)能到短時(shí)退化數(shù)據(jù)的情況,可采用加速退化試驗(yàn)的方法,通過提高產(chǎn)品所受的應(yīng)力水平,加快產(chǎn)品的退化失效過程,從而在較短時(shí)間內(nèi)即可得到更多退化信息.

References)

[1]Jayaram JSR,Girish T.Reliability prediction through degradation data modeling using a quasi-likelihood approach[C]//Proceedings Annual Reliability＆Maintainability Symposium.New York:IEEE,2005:193-199

[2]Huang W,Duane L.An alternative degradation reliability modeling approach maximum likelihood estimation[J].IEEE Transactions on Reliability,2005,54(2):310-317

[3]Wilson SP,Taylor D.Reliability assessment from fatigue micro crack data[J].IEEE Transactions on Reliability,1997,46(2):165-171

[4]鄧愛民,陳循,張春華,等.基于性能退化數(shù)據(jù)的可靠性評(píng)估[J].宇航學(xué)報(bào),2006,27(3):546-552 Deng Aimin,Chen Xun,Zhang Chunhua,Reliability assessment based on performance degradation data[J].Journal of Astronautics,2006,27(3):546-552(in Chinese)

[5]Wang Li,Li Xiaoyang,Jiang Tongm in,et al.SLD constant-stress ADT data analysis based on time series method[C]//Proceedings 8th International Conference on Reliability Maintainability and Safety.Chengdu,China:IEEE Computer Society,2009:1313-1317

[6]尤琦,趙宇,馬小兵.產(chǎn)品性能可靠性評(píng)估的時(shí)序分析方法[J].北京航空航天大學(xué)學(xué)報(bào),2009,35(5):644-648 You Qi,Zhao Yu,Ma Xiaobing.Performance reliability assessment for products based on time series analysis[J].Journal of Beijing University of Aeronautics and Astronautics,2009,35(5):644-648(in Chinese)

[7]傅惠民.相關(guān)系數(shù)平穩(wěn)過程方法[J].機(jī)械強(qiáng)度,2002,24(3):400-404 Fu Huim in.Method for correlation coefficient stationary process[J].Journal of Mechanical Strength,2002,24(3):400-404(in Chinese)

[8]傅惠民,劉成瑞,馬小兵.時(shí)間序列均值和方差函數(shù)的確定方法[J].機(jī)械強(qiáng)度,2004,26(2):164-169 Fu Huimin,Liu Chengrui,Ma Xiaobing.Method for determining functions of time seriesmean and variance[J].Journal of Mechanical Strength,2004,26(2):164-169(in Chinese)

(編 輯 :婁 嘉)

Life prediction of product based on degradation amount distribution time series analysis

Wang Li Li Xiaoyang Jiang Tongmin

(School of Reliability and Systems Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)

For the problem that traditional life prediction methods described product degradation stochastic process inadequately using monotonous regression function,the non-stationary time series types of different degradation amount distribution parameters estimation was analyzed based on product degradation amount distribution at each time.Different parameters estimation was modeled respectively usingnon-stationary time series analysis method.Hence,a realistic product life and reliability prediction method based on degradation amount distribution time series analysis was proposed.A degradation test of a certain electronic product was conducted,and product life and reliability were obtained by the suggested method.The results show that the suggested life prediction method is more reasonable than traditional life prediction method.

performance degradation;life prediction;degradation amount distribution;non-stationary time series analysis

TB 114.3

A

1001-5965(2011)04-0492-07

2010-08-10

國(guó)家安全重大基礎(chǔ)研究項(xiàng)目(61312803)

王 立(1983-),女,北京人,博士生,xiaolizi1983@hotmail.com.