劉 鑫，賈云獻(xiàn)，李孟涵，周 杰

(1.陸軍工程大學(xué)石家莊校區(qū)，河北 石家莊 050003; 2.軍事科學(xué)院系統(tǒng)工程研究院，北京 100082;3.河北工業(yè)大學(xué)，天津 300400)

0 引 言

工程實(shí)際中，系統(tǒng)的性能會(huì)隨時(shí)間不斷退化最終導(dǎo)致系統(tǒng)失效，即發(fā)生退化失效。同時(shí)，由于其外部或者內(nèi)部環(huán)境的變化，系統(tǒng)不可避免地會(huì)受到各種沖擊，如果沖擊量超過(guò)其承受閾值，會(huì)導(dǎo)致系統(tǒng)的突發(fā)失效。大部分的系統(tǒng)都同時(shí)存在著退化失效和突發(fā)失效兩種失效模式，最終系統(tǒng)的失效是兩種失效模式競(jìng)爭(zhēng)的結(jié)果。例如電阻失效可能是高壓擊穿、開(kāi)路或阻值漂移中的一種，前兩種失效模式是突發(fā)失效，后一種則屬于退化失效[1]，最先發(fā)生的失效模式造成了電阻的失效，所以電阻的失效是各種失效模式競(jìng)爭(zhēng)的結(jié)果。

競(jìng)爭(zhēng)失效現(xiàn)象的普遍存在引起了學(xué)者們的廣泛關(guān)注，目前已有許多關(guān)于競(jìng)爭(zhēng)失效模型的研究與應(yīng)用[2-5]。在理論方面，Kristoffer等[6]研究了競(jìng)爭(zhēng)失效模式下的疲勞失效機(jī)理；Zhang等[7]研究了競(jìng)爭(zhēng)失效的加速壽命試驗(yàn)方法；Song等[8]研究了多個(gè)相互影響的部件所組成的系統(tǒng)的可靠性建模。在應(yīng)用方面，王海琨等[9]研究了基于競(jìng)爭(zhēng)失效過(guò)程的剩余壽命預(yù)測(cè)問(wèn)題；Zhang等[10]研究了基于競(jìng)爭(zhēng)失效模型的預(yù)防性維修的維修決策優(yōu)化。同時(shí)，競(jìng)爭(zhēng)失效模型在航空發(fā)動(dòng)機(jī)[11]，金屬化膜脈沖電容器[12]和數(shù)控機(jī)床[13]等裝備或產(chǎn)品的可靠性研究中也得到廣泛應(yīng)用。然而在競(jìng)爭(zhēng)失效模型研究中，對(duì)實(shí)際情況進(jìn)行了大量簡(jiǎn)化，使得計(jì)算結(jié)果與系統(tǒng)實(shí)際可靠性存在一定差距，在工程中可能會(huì)造成可靠度評(píng)估偏高而造成損失。基于此，探索一種貼近系統(tǒng)實(shí)際工況的競(jìng)爭(zhēng)失效模型具有重要研究意義。

針對(duì)此問(wèn)題，黃文平等[14]考慮了退化過(guò)程對(duì)沖擊過(guò)程的影響，即隨著退化量增加，系統(tǒng)抵御沖擊能力會(huì)下降，導(dǎo)致突發(fā)失效閾值產(chǎn)生變化，在此基礎(chǔ)上建立了變失效閾值的競(jìng)爭(zhēng)失效模型。本文在變閾值競(jìng)爭(zhēng)失效模型的基礎(chǔ)上，考慮存在多種不同沖擊類型的情況，對(duì)競(jìng)爭(zhēng)失效模型的適用范圍進(jìn)一步擴(kuò)展。在分析多種沖擊類型同時(shí)存在時(shí)對(duì)系統(tǒng)可靠性影響的基礎(chǔ)上，結(jié)合極值沖擊模型給出競(jìng)爭(zhēng)失效可靠性建模的一般方法。最后，通過(guò)某微型發(fā)動(dòng)機(jī)的可靠性評(píng)估案例，證明該方法的有效性。

1 考慮多類型沖擊的競(jìng)爭(zhēng)失效過(guò)程研究

在可靠性研究中，常用沖擊模型來(lái)表示系統(tǒng)受到的沖擊過(guò)程[15]，本文以極值沖擊模型為例來(lái)表示沖擊過(guò)程。在極值沖擊模型中，當(dāng)某一次沖擊量超過(guò)預(yù)先設(shè)定的突發(fā)失效閾值時(shí)，則認(rèn)為系統(tǒng)失效。在基于極值沖擊模型的競(jìng)爭(zhēng)失效中，系統(tǒng)退化量Xs(t)由兩部分組成，一是系統(tǒng)連續(xù)性退化過(guò)程，二是系統(tǒng)由于沖擊過(guò)程所造成的退化量，當(dāng)二者累積之和超過(guò)系統(tǒng)的退化失效閾值H，則系統(tǒng)失效，稱為退化失效；或者某一次的沖擊量或沖擊強(qiáng)度超過(guò)突發(fā)失效閾值D，系統(tǒng)同樣失效，稱為突發(fā)失效?；诖?，競(jìng)爭(zhēng)失效過(guò)程如圖1所示，圖1（a）中，在系統(tǒng)的退化過(guò)程中總共經(jīng)歷了3次沖擊過(guò)程（沖擊量分別為W1、W2和W3，對(duì)應(yīng)退化量Y1、Y2和Y3），且第3次沖擊量超過(guò)系統(tǒng)閾值，此時(shí)，雖然退化量并未達(dá)到系統(tǒng)退化失效閾值H，但是系統(tǒng)依然失效，此即為由沖擊過(guò)程造成的突發(fā)失效；圖1（b）中，幾次沖擊均未超過(guò)系統(tǒng)閾值，只是造成退化量的增長(zhǎng)，最后系統(tǒng)總退化量超過(guò)系統(tǒng)退化失效閾值H導(dǎo)致退化失效。

系統(tǒng)在退化過(guò)程中，隨著退化量增加，其抵御沖擊能力會(huì)相應(yīng)降低，所以在系統(tǒng)退化到一定程度且未失效的情況下，其沖擊失效閾值應(yīng)該相應(yīng)降低。同時(shí)，系統(tǒng)在運(yùn)行過(guò)程中會(huì)受到多種類型的沖擊影響，比如筆記本電腦的電池，其連續(xù)性退化過(guò)程主要是由于內(nèi)部化學(xué)反應(yīng)導(dǎo)致充放電能力變?nèi)?，而影響電池退化量的沖擊過(guò)程則包括過(guò)熱、電壓過(guò)高等因素。因此，本文在變閾值競(jìng)爭(zhēng)失效模型基礎(chǔ)上研究同時(shí)存在兩種或者以上不同沖擊過(guò)程競(jìng)爭(zhēng)失效模型。假設(shè)系統(tǒng)在退化過(guò)程中，當(dāng)總退化量達(dá)到某一設(shè)定的閾值L（L＜H）時(shí)，由于系統(tǒng)抵御外界沖擊的能力減弱，突發(fā)失效的閾值相應(yīng)由D1降為D2，此時(shí)模型主要有3種失效形式，如圖2所示。圖2（a）中表示在退化過(guò)程中，系統(tǒng)受到兩種沖擊影響，第一種沖擊出現(xiàn)了兩次（W1,1和W1,2），第二種沖擊出現(xiàn)一次（W2,1），每次沖擊的沖擊量均未超出最大失效閾值（分別為D1,1和D2,1），最終系統(tǒng)由于總的退化量達(dá)到了失效閾值H而失效，并未出現(xiàn)突發(fā)失效。圖2（b）中，受到同樣的沖擊過(guò)程，雖然系統(tǒng)總退化量遠(yuǎn)未達(dá)到其失效閾值，但是第三次沖擊過(guò)程超過(guò)了失效閾值D2,1而導(dǎo)致系統(tǒng)突發(fā)失效。圖2（c）中，系統(tǒng)在t0時(shí)刻達(dá)到了設(shè)定的閾值L，此時(shí)系統(tǒng)突發(fā)失效的閾值變?yōu)镈1,2和D2,2，接著在t4時(shí)刻，由于沖擊量超過(guò)了閾值D1,2而導(dǎo)致系統(tǒng)突發(fā)失效。

圖1 基于極值沖擊模型的競(jìng)爭(zhēng)失效過(guò)程

2 系統(tǒng)競(jìng)爭(zhēng)失效過(guò)程的可靠性建模

為建立競(jìng)爭(zhēng)失效過(guò)程的可靠性模型，根據(jù)上文分析，需要同時(shí)研究系統(tǒng)的退化失效過(guò)程和突發(fā)失效過(guò)程。本文首先建立系統(tǒng)總退化量的退化失效模型，在此基礎(chǔ)上，基于變閾值模型，根據(jù)突發(fā)失效閾值的變化，分兩類情況討論系統(tǒng)競(jìng)爭(zhēng)失效過(guò)程的可靠性建模。

2.1 退化失效模型

對(duì)于任何系統(tǒng)的退化失效過(guò)程，當(dāng)總退化量達(dá)到失效閾值H時(shí)，則認(rèn)為系統(tǒng)失效。為建立系統(tǒng)退化失效模型，首先需要計(jì)算系統(tǒng)總的退化量：

式中：Xs(t)——系統(tǒng)總的退化量；

X(t)——連續(xù)退化過(guò)程的退化量；

S(t)——沖擊過(guò)程所造成的退化量。

系統(tǒng)連續(xù)性退化過(guò)程在工程實(shí)際中可能服從任意規(guī)律，例如服從Wiener過(guò)程或者Gamma過(guò)程等，為簡(jiǎn)化建模過(guò)程，假設(shè)系統(tǒng)連續(xù)性退化過(guò)程為線性退化過(guò)程，即：

其中，φ表示退化過(guò)程的初始退化量，β為退化速率，假設(shè)其服從正態(tài)分布，即

假設(shè)該系統(tǒng)共有l(wèi)種不同的沖擊類型，對(duì)于同一類型沖擊，所有沖擊過(guò)程都為到達(dá)率λj的泊松過(guò)程。每次沖擊過(guò)程造成的退化量為Yj,k（第j種沖擊類型第k次出現(xiàn)），且假設(shè)沖擊量為服從正態(tài)分布的獨(dú)立同分布非負(fù)隨機(jī)變量，即此時(shí)，外界沖擊所造成的退化量為：

其中Nj(t)為第j種沖擊類型出現(xiàn)的次數(shù)。

綜合以上分析，系統(tǒng)不發(fā)生退化失效的概率為：

根據(jù)條件概率的性質(zhì)，可以將其表示為：

由于已經(jīng)假設(shè)沖擊的退化量服從正態(tài)分布，沖擊過(guò)程服從泊松分布，則可以計(jì)算系統(tǒng)不發(fā)生退化失效的概率為：

2.2 競(jìng)爭(zhēng)失效可靠性建模

在競(jìng)爭(zhēng)失效模型中，需要在連續(xù)性退化過(guò)程中同時(shí)考慮突發(fā)失效，在本文中即沖擊過(guò)程。假設(shè)Wj,k（第j種沖擊類型第k次的沖擊量）為沖擊量的大小，且服從獨(dú)立同分布的正態(tài)分布則系統(tǒng)不發(fā)生突發(fā)失效概率為：

在變失效閾值的情況下，同時(shí)考慮多種沖擊類型的存在，系統(tǒng)的可靠性建模十分復(fù)雜，因此，根據(jù)失效閾值的變化將競(jìng)爭(zhēng)失效過(guò)程的可靠性計(jì)算分為兩種情況討論。

1）系統(tǒng)的性能退化總量小于L，即突發(fā)失效的閾值為Dj,1，此時(shí)系統(tǒng)的可靠度函數(shù)即為總退化量小于L，且每種類型的沖擊每次出現(xiàn)時(shí)的沖擊量都小于突發(fā)失效閾值Dj,1的概率。根據(jù)前面的計(jì)算，可以得到：

最后得到可靠度為：

2）系統(tǒng)的退化量大于或者等于L且小于H，此時(shí)系統(tǒng)的可靠性函數(shù)為：

其中Rshock(t)為在沖擊影響下系統(tǒng)不發(fā)生沖擊失效的可靠度函數(shù)。假設(shè)系統(tǒng)共經(jīng)歷各種沖擊過(guò)程n=次，且u(u＜t)為退化量到達(dá)L的時(shí)刻。令pj為[0,u]上第j類沖擊發(fā)生的次數(shù)，qj為(u,t)上第j類沖擊發(fā)生的次數(shù)，且根據(jù)泊松過(guò)程的性質(zhì)，滿足上面情況的概率為：

令FX(N(u),L,u)為系統(tǒng)退化量小于L且外界沖擊總數(shù)為時(shí)不發(fā)生失效的分布函數(shù)：

對(duì)t求導(dǎo)數(shù)，可以得到其密度函數(shù)為：

綜合以上分析，可以得到其可靠性為：

其中pj+qj=mj,?j。

因?yàn)榇藘煞N狀態(tài)為互斥事件，因此可靠度為二者之和，即：

為了對(duì)所建立的可靠性模型進(jìn)行驗(yàn)證，選取沖擊次數(shù)為0這一特例進(jìn)行檢驗(yàn)。此時(shí)m1=m2=···=mj=0，即不存在任何沖擊過(guò)程，此時(shí)的模型應(yīng)該為變失效閾值的競(jìng)爭(zhēng)失效模型。將其帶入建立的可靠性模型中，可以得到：

此時(shí)，可以得到系統(tǒng)的可靠度為：

此結(jié)果與文獻(xiàn)[14]的模型在形式上是一致的，而區(qū)別在于不發(fā)生外界沖擊的概率是不同的，造成這一結(jié)果的原因是本文的模型考慮了多種類型的沖擊同時(shí)存在的情況。

3 案例分析

本文選取桑迪亞國(guó)家實(shí)驗(yàn)室某微型發(fā)動(dòng)機(jī)為研究對(duì)象，對(duì)所建立的模型進(jìn)行驗(yàn)證。該微型發(fā)動(dòng)機(jī)在運(yùn)行過(guò)程中同時(shí)存在著磨損退化與沖擊過(guò)程，由摩擦造成的磨損退化是微型發(fā)動(dòng)機(jī)失效的主要模式，同時(shí)沖擊會(huì)產(chǎn)生一定量的磨損碎屑，當(dāng)沖擊量大于一定閾值時(shí)也會(huì)導(dǎo)致其失效。該微型發(fā)動(dòng)機(jī)的失效過(guò)程同時(shí)包括連續(xù)退化過(guò)程和沖擊過(guò)程，為典型的競(jìng)爭(zhēng)失效過(guò)程。微型發(fā)動(dòng)機(jī)的數(shù)據(jù)主要來(lái)自文獻(xiàn)[16]，為了對(duì)本文的模型進(jìn)行驗(yàn)證，在原有數(shù)據(jù)基礎(chǔ)上，假設(shè)第二種沖擊類型的存在，為了簡(jiǎn)化計(jì)算，假設(shè)只存在兩種類型的沖擊過(guò)程，且兩種沖擊類型服從到達(dá)率相同的泊松分布，同時(shí)對(duì)微型發(fā)動(dòng)機(jī)的變失效閾值進(jìn)行假設(shè)，其主要參數(shù)如表1所示。

根據(jù)本文提出的考慮多種沖擊類型的競(jìng)爭(zhēng)失效可靠性模型，計(jì)算微型發(fā)動(dòng)機(jī)的可靠度函數(shù)如圖3所示?？梢钥吹剑?0 000轉(zhuǎn)之前，發(fā)動(dòng)機(jī)的可靠度下降相對(duì)緩慢，仍然保持較高的可靠度水平，在50 000轉(zhuǎn)左右，發(fā)動(dòng)機(jī)的可靠度急劇下降，這主要是由于系統(tǒng)退化量達(dá)到變閾值水平，發(fā)動(dòng)機(jī)抵御外界沖擊的能力變?nèi)?，失效閾值降低，所以可靠度迅速下降，?5 000轉(zhuǎn)之后，發(fā)動(dòng)機(jī)性能又相對(duì)平穩(wěn)，可靠度下降又相對(duì)緩慢。文獻(xiàn)[14]中的模型只考慮了變閾值，沒(méi)有考慮存在多種沖擊類型的情況，且由于本案例中第二類沖擊造成的退化量比第一類要大，導(dǎo)致其各個(gè)時(shí)刻的可靠度評(píng)估結(jié)果都比本文模型的結(jié)果高。

為更清楚地進(jìn)行對(duì)比，二者相同時(shí)刻可靠度差值如圖4所示，可以看到，可靠度差值隨時(shí)間迅速增大，然后又逐漸減小，差值在45 000~55 000轉(zhuǎn)達(dá)到最大，主要原因和上文中分析相同，為發(fā)動(dòng)機(jī)退化量達(dá)到規(guī)定閾值，導(dǎo)致沖擊閾值下降，由于本文模型假設(shè)存在兩種沖擊類型，且第二類沖擊造成的沖擊量均值比第一種要大，所以退化量會(huì)首先到達(dá)變閾值臨界點(diǎn)，導(dǎo)致發(fā)動(dòng)機(jī)可靠度迅速下降。本文模型由于綜合考慮了多種因素的影響，與系統(tǒng)真實(shí)的工作環(huán)境更加接近，得到的可靠度結(jié)果也更為全面，可以有效避免對(duì)可靠性的過(guò)高估計(jì)而帶來(lái)的損失。

表1 某微型發(fā)動(dòng)機(jī)競(jìng)爭(zhēng)失效過(guò)程參數(shù)

圖3 微型發(fā)動(dòng)機(jī)壽命可靠度函數(shù)

圖4 基于不同模型的可靠度差值

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

本文針對(duì)競(jìng)爭(zhēng)失效模型連續(xù)性退化過(guò)程中同時(shí)存在多種類型沖擊問(wèn)題，在變閾值競(jìng)爭(zhēng)失效模型基礎(chǔ)上，研究了系統(tǒng)可靠性建模。首先，研究了基于極值沖擊模型的競(jìng)爭(zhēng)失效，在此基礎(chǔ)上，分析了多種沖擊類型同時(shí)存在時(shí)對(duì)系統(tǒng)可靠性的影響，給出了競(jìng)爭(zhēng)失效模型主要失效模式；最后，建立了該模型的解析表達(dá)式。主要結(jié)論如下：1）本文建立的可靠性模型，減少了理想化條件，提高了模型的適用性和合理性；2）考慮存在多種類型沖擊存在的情況，可以得到更為全面的評(píng)估結(jié)果，提高了可靠性評(píng)估精度，避免了對(duì)系統(tǒng)可靠性的樂(lè)觀估計(jì)；3）以某微型發(fā)動(dòng)機(jī)為研究對(duì)象，驗(yàn)證了該模型的有效性，具有一定的工程應(yīng)用和推廣價(jià)值。

本文在競(jìng)爭(zhēng)失效的可靠性建模中，因?yàn)榻＿^(guò)程十分復(fù)雜，為簡(jiǎn)化計(jì)算，均假設(shè)裝備壽命服從正態(tài)分布。因此，如何建立適用范圍廣，適合各種系統(tǒng)的競(jìng)爭(zhēng)失效可靠性模型，值得進(jìn)一步探索研究。