隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)可靠性評(píng)估

潘　剛， 尚朝軒， 梁玉英， 蔡金燕， 孟亞峰

(軍械工程學(xué)院 電子與光學(xué)工程系，石家莊　050003)

隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)可靠性評(píng)估

潘剛， 尚朝軒， 梁玉英， 蔡金燕， 孟亞峰

(軍械工程學(xué)院 電子與光學(xué)工程系，石家莊050003)

摘要：由于環(huán)境或其他原因，部件可能受到隨機(jī)沖擊，部件性能退化由正常性能退化和隨機(jī)沖擊兩部分構(gòu)成，對(duì)于該類高可靠部件短時(shí)間內(nèi)很難得到足夠的性能數(shù)據(jù)，致使對(duì)部件認(rèn)知存在一定不確定性，無(wú)法準(zhǔn)確估計(jì)系統(tǒng)的可靠性。為實(shí)現(xiàn)對(duì)系統(tǒng)可靠性的準(zhǔn)確估計(jì)，假定沖擊引起的部件性能損傷分布參數(shù)為區(qū)間變量，建立基于區(qū)間變量的部件性能分布模型，給出了部件狀態(tài)性能區(qū)間連續(xù)序列定義和區(qū)間狀態(tài)概率計(jì)算方法，對(duì)傳統(tǒng)的通用生成函數(shù)方法進(jìn)行改進(jìn)，定義了區(qū)間通用生成函數(shù)及其運(yùn)算法則，提出了隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)可靠性評(píng)估方法，并以仿真實(shí)例進(jìn)行了驗(yàn)證說明。該方法不僅克服了性能分布信息缺少，無(wú)法準(zhǔn)確建立狀態(tài)性能分布模型的不足，且具有很強(qiáng)的通用性和工程應(yīng)用價(jià)值。

關(guān)鍵詞：隨機(jī)沖擊；認(rèn)知不確定性；區(qū)間連續(xù)序列；區(qū)間通用生成函數(shù)

20世紀(jì)70年代多態(tài)系統(tǒng)的概念被提出后[5-6]，基于多態(tài)系統(tǒng)模型的系統(tǒng)可靠性研究得到了學(xué)者們的廣泛關(guān)注。在理論方法方面，文獻(xiàn)[4, 11-12]對(duì)多態(tài)系統(tǒng)可靠性的基本概念、評(píng)估方法和優(yōu)化設(shè)計(jì)等方面內(nèi)容進(jìn)行了詳細(xì)的闡述。在工程應(yīng)用方面，多態(tài)系統(tǒng)可靠性的相關(guān)理論已經(jīng)應(yīng)用到電力[7-8]、機(jī)械[9-10]等領(lǐng)域。一些學(xué)者對(duì)上述傳統(tǒng)的多態(tài)系統(tǒng)理論進(jìn)行推廣。Ding等[13-14]在模糊狀態(tài)和模糊狀態(tài)概率給定的前提下，給出了模糊多態(tài)系統(tǒng)的通用定義和可靠性分析方法。鄢民強(qiáng)等[15]提出了一種考慮不完全覆蓋的模糊多狀態(tài)系統(tǒng)的可靠性計(jì)算方法。Li等[16]利用區(qū)間分析理論和通用生成函數(shù)方法分析了多態(tài)系統(tǒng)的區(qū)間可靠性。Li等[17]對(duì)具有隨機(jī)沖擊情況的考慮競(jìng)爭(zhēng)失效的多狀態(tài)系統(tǒng)可靠性模型進(jìn)行分析。對(duì)于具有隨機(jī)沖擊情況的高可靠部件構(gòu)成的系統(tǒng)，通常情況下，很難準(zhǔn)確定量沖擊對(duì)部件或系統(tǒng)性能的影響。通過對(duì)當(dāng)前上述研究分析可知，在具有隨機(jī)沖擊的情況下，無(wú)論是基于 “二態(tài)”假設(shè)，還是采用多態(tài)系統(tǒng)模型，對(duì)于隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)可靠性分析研究相對(duì)較少。

高可靠性部件構(gòu)成的多態(tài)系統(tǒng)進(jìn)行可靠性分析時(shí)，隨機(jī)沖擊情況下的認(rèn)知不確定性主要表現(xiàn)在兩個(gè)方面：①隨機(jī)沖擊對(duì)部件性能分布的影響具有一定的不確定性，且很難準(zhǔn)確定量；②由于對(duì)系統(tǒng)或部件性能退化機(jī)理認(rèn)知的不完善，而建立的狀態(tài)性能分布模型往往是不準(zhǔn)確的，有時(shí)可能會(huì)出現(xiàn)很大的偏差。鑒于上述不足，本文提出了隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)區(qū)間可靠性分析方法，給出了多態(tài)系統(tǒng)的可靠性分析流程，并以仿真實(shí)例進(jìn)行了驗(yàn)證分析。

1性能退化系統(tǒng)狀態(tài)分析

1.1性能退化部件性能分析

對(duì)隨機(jī)沖擊情況下考慮認(rèn)知不確定的部件性能分布進(jìn)行分析之前，先做如下假設(shè)：

(2) 部件只有一個(gè)性能參數(shù)x，對(duì)應(yīng)一個(gè)性能退化過程，且退化過程不可逆。

(3) 定義N(t)為到t時(shí)刻為止部件受到的沖擊次數(shù)，Sl為第l次沖擊對(duì)部件性能的損傷，S(0)=0，S(t)為到t時(shí)刻為止，隨機(jī)沖擊對(duì)部件性能產(chǎn)生的累計(jì)損傷，則：

假設(shè)沖擊次數(shù)服從參數(shù)為λ的泊松分布，即

(5) 在任意時(shí)刻t，假設(shè)部件性能X(t)服從均值為μ(t)，方差為σ(t)2的正態(tài)分布，其概率密度函數(shù)為f(x;θ(t))。

(6) 同時(shí)考慮正常的性能退化和沖擊所引起的性能退化時(shí)，t時(shí)刻部件的性能可以表示為：

此時(shí)部件的性能分布函數(shù)為：

F(Y)=p{G(t)≤y}=

由于，

最后可得:

(1)

式中，

1.2性能退化部件狀態(tài)概率分析

定義部件的狀態(tài)性能區(qū)間劃分準(zhǔn)則，采用系統(tǒng)最低任務(wù)需求[w]與部件性能分布特性相結(jié)合的方法對(duì)部件的狀態(tài)性能區(qū)間進(jìn)行劃分，具體如下：

圖1　部件性能分布Fig.1 Performance distribution of components

(4) 令部件的最大性能為gmax，為了便于表示假定gε1(t)=μz(t)-3σz(t)，gε2(t)=μz(t)-2σz(t)，gε3(t)=μz(t)-σz(t)，gε4(t)=μz(t)，gε5(t)=μz(t)+σz(t)，gε6(t)=μz(t)+2σz(t)，gε7(t)=μz(t)+3σz(t)。

其中，

其中，

gεx(t)=(gε1(t)+gεl(t))/2，l=1,…,7

其中，

gεx(t)=(gεr+1(t)+gεl(t))/2，1≤r

其中，

其中，

⑥ 當(dāng)[w]>gmax時(shí)，表明部件的性能已不滿足系統(tǒng)分配的最小任務(wù)需求，對(duì)其進(jìn)行狀態(tài)的劃分已經(jīng)沒有實(shí)際的工程意義，故此時(shí)可認(rèn)為部件失效。

(2)

式中，

2運(yùn)算法則

為了降低計(jì)算的難度，在分析時(shí)將狀態(tài)性能區(qū)間連續(xù)序列分為區(qū)間連續(xù)序列上限、區(qū)間連續(xù)序列中值、區(qū)間連續(xù)序列下限三種情況進(jìn)行討論分析，現(xiàn)以區(qū)間連續(xù)序列上限為例進(jìn)行分析說明。

2.1區(qū)間通用生成函數(shù)定義

部件i的狀態(tài)概率區(qū)間族形式為：

{pi(t)}={[p1(t)],…,[pki(t)],…,[pMi(t)]}

定義部件i的區(qū)間通用生成函數(shù)為：

(3)

式中，i=1,2,…,n，ki=1,2,…,Mi。

2.2區(qū)間通用生成函數(shù)運(yùn)算法則

U(t,z)=Ω(ui(t,z),uj(t,z))=

(4)

根據(jù)系統(tǒng)結(jié)構(gòu)特點(diǎn)定義如下運(yùn)算符：

δ1(ui(t,z),uj(t,z))=

δ2(ui(t,z),uj(t,z))=

3基于區(qū)間通用生成函數(shù)的多態(tài)系統(tǒng)可靠性分析

根據(jù)上述運(yùn)算法則的分析，假定得到的多態(tài)系統(tǒng)區(qū)間通用生成函數(shù)為：

(5)

定義多態(tài)系統(tǒng)的最小區(qū)間任務(wù)性能需求為[w]，則系統(tǒng)的區(qū)間可靠度為：

[R(t)]=P{[G(t)]≥[w]}=

(6)

式中，P{[G(t)]≥[w]}表示[G(t)]-[w]≥0的概率，[G(t)]={[g1(t)],…,[gks(t)],…,[gMs(t)]}。

綜合上述三點(diǎn)，可以將[gks(t)]-[w]的可能度區(qū)間定義為：

p([gks(t)]-[w]≥0)=

因此，系統(tǒng)的區(qū)間可靠度為：

[R]=P{[G(t)]≥[w]}=

4算例分析

圖2　多態(tài)串-并聯(lián)系統(tǒng)Fig.2 Multi-state serial-parallel system

4.1部件和分系統(tǒng)區(qū)間通用生成函數(shù)

4.1.1部件區(qū)間通用生成函數(shù)

表1　部件退化性能、分系統(tǒng)結(jié)構(gòu)函數(shù)和系統(tǒng)結(jié)構(gòu)函數(shù)

部件6可分為2個(gè)狀態(tài)，其狀態(tài)性能定義為g6={[0,70],[70,100]}。

在t=10 000 h時(shí)，根據(jù)2.2節(jié)性能退化部件狀態(tài)概率分析，部件各狀態(tài)性能區(qū)間連續(xù)序列對(duì)應(yīng)的概率區(qū)間為：{p1(t)}={p1,1(t),p1,2(t),p1,3(t)}，其中，p1,1(t)=[0.004 4,0.008 2]，p1,2(t)=[0.064 5,0.107 8]，p1,3(t)=[0.836 2,0.884 0]{p2(t)}={p1(t)}，{p3(t)}={p31(t),p32(t)}，其中，p3,1(t)=[0.037 1,0.049 5]p3,2(t)=[0.935 0,0.950 5]{p3(t)}={p4(t)}={p5(t)}，{p6(t)}={p6,1(t),p6,2(t)}，p6,1(t)=0.051 8，p6,2(t)=0.948 2。

以狀態(tài)性能區(qū)間連續(xù)序列上限方法為例對(duì)部件區(qū)間通用生成函數(shù)和分系統(tǒng)區(qū)間通用生成函數(shù)進(jìn)行分析，其他的方法類似，不再贅述。

根據(jù)狀態(tài)性能區(qū)間連續(xù)序列族定義和性能退化部件狀態(tài)概率求解方法，可得部件1～6的區(qū)間通用生成函數(shù)為：

u1(t,z)=[0.004 4,0.008 2]z[51.97,52.47]+

[0.064 5,0.107 8]z[54.75,55.25]+

[0.836 2,0.884 0]z[69.75,70.25]，

u2(t,z)=u1(t,z),

u3(t,z)=[0.037 1,0.049 5]z[34.77,35.07]+

[0.935 0,0.950 5]z[49.85,50.15]，

u3(t,z)=u4(t,z)=u5(t,z)，

u6(t,z)=[0,70]z0.051 8+[70,100]z0.948 2

4.1.2分系統(tǒng)區(qū)間通用生成函數(shù)

由圖2可知分系統(tǒng)1是由部件1和部件2構(gòu)成，且分系統(tǒng)1的性能為部件1和部件2性能的和，根據(jù)第3節(jié)運(yùn)算法則定義的δ1運(yùn)算符，可得分系統(tǒng)1的區(qū)間通用生成函數(shù)為：

Usub1(z,t)=δ1(u1(z,t),u2(z,t))=

[0.000 6,0.001 8]z[106.7,107.7]+

[0.004 2,0.011 6]z[109.5,110.5]+

[0.007 4,0.014 4]z[121.7,122.7]+

[0.107 8,0.190 6]z[124.5,125.5]+

[0.699 2,0.781 6]z[139.5,140.5]

分系統(tǒng)2是由部件3、4和5構(gòu)成，且其性能為部件3、4和5性能的和，同理可得分系統(tǒng)2的區(qū)間通用生成函數(shù)為：

Usub2(z,t)=δ1(u3(z,t),u4(z,t),u5(z,t))=

[0,0.000 1]z[104.3,105.2]+[0.003 9,0.007 0]z[119.4,120.3]+

[0.097 3,0.134 2]z[134.5,135.4]+

[0.817 3,0.858 7]z[149.6,150.5]

4.2系統(tǒng)區(qū)間可靠性分析

4.2.1系統(tǒng)區(qū)間可靠性分析

將分系統(tǒng)1和分系統(tǒng)2采用第2節(jié)運(yùn)算法則定義的δ2運(yùn)算符，可得分系統(tǒng)1、2綜合后的分系統(tǒng)區(qū)間通用生成函數(shù)為：

U1-2=[0.000 4,0.001 5]z[106.7,107.7]+

[0.003 7,0.011 3]z[109.5,110.5]+

[0.003 1,0.006 9]z[119.4,120.3]+

[0.006 8,0.014 3]z[121.7,122.7]+

[0.098 6,0.189 2]z[124.5,125.5]+

[0.068 0,0.104 9]z[134.5,135.4]+

[0.571 4,0.671 2]z[139.5,140.5]

因此，根據(jù)式(7)可得分系統(tǒng)1和2的區(qū)間可靠度為：

R1-2=[0.752 1,0.999 4]

分系統(tǒng)3的可靠度為：

R3=0.948 2

最后可得t=10 000 h時(shí)系統(tǒng)的區(qū)間可靠度為：

R(t)=R1-2(t)R3(t)=[0.713 2,0.947 7]

在t=10 000 h時(shí)采用Monte Carlo(MC)仿真方法進(jìn)行分析，在分析過程中為了充分描述Monte Carlo(MC)仿真方法在每次仿真計(jì)算過程中所得計(jì)算結(jié)果的可能性，取L=1 000次仿真結(jié)果的最大值和最小值，將其作為區(qū)間可靠度的最大值和最小值，可得MC仿真方法的區(qū)間可靠度為[0.894 8,0.956 9]。

如果采用傳統(tǒng)的方法進(jìn)行系統(tǒng)可靠性分析，則一般假設(shè)部件僅有兩個(gè)狀態(tài)，根據(jù)系統(tǒng)總的性能輸出要求，確定部件的失效閾值。對(duì)于部件1,2性能不小于[54.75, 55.25]時(shí)分系統(tǒng)1正常，部件3,4,5性能不小于[35.75,36.25]時(shí)分系統(tǒng)2正常，部件6的性能不小于70時(shí)，分系統(tǒng)3正常，因此，采用傳統(tǒng)方法可得系統(tǒng)的區(qū)間可靠度為[0.483 4,0.606 3]。

采用本文所提方法、傳統(tǒng)可靠性方法、Monte Carlo仿真方法所得系統(tǒng)的區(qū)間可靠性對(duì)比關(guān)系如表2所示。

表2　三種方法所得系統(tǒng)區(qū)間可靠性對(duì)比分析結(jié)果

通過對(duì)表2分析我們可以得出以下結(jié)論：

1) 采用本文所提方法，分別對(duì)連續(xù)區(qū)間序列下限、中值、上限進(jìn)行分析，所得計(jì)算結(jié)果總體上呈增大的趨勢(shì)，連續(xù)區(qū)間序列的上限方法和下限均未完全包含連續(xù)區(qū)間序列中的狀態(tài)信息，而連續(xù)區(qū)間序列的中值為上限和下限的折中，故所得結(jié)果一定程度上，更加接近客觀值，取連續(xù)區(qū)間序列族下限時(shí)所得結(jié)果隨著時(shí)間的增加與后兩者的偏差越大，在后續(xù)研究中不再分析。

2) 傳統(tǒng)方法、區(qū)間連續(xù)序列中值方法和Monte Carlo仿真方法三者所得結(jié)果進(jìn)行對(duì)比，傳統(tǒng)方法與后兩者相比誤差較大，且采用傳統(tǒng)方法進(jìn)行系統(tǒng)可靠性估計(jì)時(shí)，部件的失效閾值通常是根據(jù)部件自身的統(tǒng)計(jì)規(guī)律確定的，不能根據(jù)具體系統(tǒng)具體分析，有時(shí)可能會(huì)出現(xiàn)一些偏差，進(jìn)一步體現(xiàn)了多態(tài)系統(tǒng)理論在描述系統(tǒng)可靠性時(shí)的準(zhǔn)確性優(yōu)勢(shì)。

4.2.2狀態(tài)性能區(qū)間半徑影響分析

圖3　區(qū)間半徑大小與系統(tǒng)區(qū)間可靠度均值的關(guān)系Fig.3 The relationship between interval radius and the mean of system interval reliability

由圖3可以得出以下結(jié)論：

1) 連續(xù)區(qū)間序列上限方法的最大值點(diǎn)和最小值點(diǎn)分別為(0,0,0.948 1)、(2,2,0.701 4)；連續(xù)區(qū)間序列中值方法的最大值點(diǎn)和最小值點(diǎn)分別為(0,0,0.904 7)、(2,2,0.648 3)。

2) 對(duì)于連續(xù)區(qū)間序列上限方法和連續(xù)區(qū)間序列中值方法，隨著參數(shù)α，β值的增加，系統(tǒng)區(qū)間可靠度均值呈增大趨勢(shì)，其主要原因是當(dāng)α=β=0，區(qū)間連續(xù)序列退化為點(diǎn)區(qū)間，部件狀態(tài)性能區(qū)間中存在的不確定性變?yōu)?，再對(duì)其采用區(qū)間序列上限方法時(shí)，其取值變?yōu)辄c(diǎn)區(qū)間的上邊界，致使[gks(t)]-[w]的可能度增大，從而使得區(qū)間可靠度的均值增大，同理，區(qū)間序列中值法亦是如此。

4.3進(jìn)一步分析

4.3.1部件狀態(tài)點(diǎn)區(qū)間劃分時(shí)系統(tǒng)可靠性分析

為了充分分析本文所提方法的優(yōu)缺點(diǎn)，在t=10 000 h時(shí)，將部件狀態(tài)按傳統(tǒng)的點(diǎn)區(qū)間進(jìn)行劃分，部件的通用生成函數(shù)為：

u1(z,t)=u2(z,t)=0.002 7z[0,51.6]+

0.093 4z[51.6,55]+0.903 8z[55,70]

u3(z,t)=u4(z,t)=u5(z,t)=

0.039 4z[0,34.8]+0.960 6z[34.8,50]

分系統(tǒng)1和分系統(tǒng)2的區(qū)間通用生成函數(shù)分別為：

Usub1(z,t)=δ1(u1(z,t),u2(z,t))=

0.000 5z[51.6,106.6]+0.004 9z[55,106.6]+

0.008 7z[103.2,121.6]+0.168 9z[106.6,125]+0.817 0z[110,140]

Usub2(z,t)=δ1(u3(z,t),u4(z,t),u5(z,t))=

0.000 1z[0,104.5]+0.004 5z[34.8,119.7]+

0.109 1z[69.7,134.8]+0.886 3z[104.5,150]

將分系統(tǒng)1、2綜合后的分系統(tǒng)區(qū)間通用生成函數(shù)為：

U1-2=0.000 1z[0,104.5]+0.000 5z[51.6,106.6]+

0.004 9z[55,110]+0.004 5z[34.8,119.6]+

0.007 7z[103.2,121.6]+0.149 7z[104.5,125]+

0.089 2z[69.7,134.8]+0.724 1z[104.5,140]

當(dāng)系統(tǒng)的最小性能需求[w]=[90,110]，則分系統(tǒng)12的可靠度為：R1-2=0.837 7。最后可得系統(tǒng)的可靠度為：

R(t)=R1-2(t)R3(t)=0.794 4

4.3.2進(jìn)一步對(duì)比分析

為了進(jìn)一步對(duì)比分析，將區(qū)間連續(xù)序列上限方法(狀態(tài)3-2)和區(qū)間連續(xù)序列中值方法(狀態(tài)3-2)兩種方法結(jié)果的均值，點(diǎn)區(qū)間方法、傳統(tǒng)方法和Monte Carlo仿真方法的結(jié)果進(jìn)行比較。

圖4　幾種方法所得可靠性結(jié)果對(duì)比Fig.4 Comparison of reliabilities respectively obtained by the afore-said methods

由圖4可以得出以下結(jié)論：

1) 在傳統(tǒng)方法中，部件失效閾值的選取對(duì)系統(tǒng)可靠性的影響較為敏感，因此其成為影響系統(tǒng)可靠度評(píng)估精度的薄弱環(huán)節(jié)，與其他方法相比誤差較大。

2) 區(qū)間連續(xù)序列中值方法(狀態(tài)3-2)和點(diǎn)區(qū)間方法兩種方法所得結(jié)果相對(duì)比較保守，此外狀態(tài)數(shù)相同的情況下，因大部分系統(tǒng)的失效率隨著時(shí)間的增加逐漸增大，即可靠性降低的速率要增大，故采用區(qū)間連續(xù)序列中值方法要比點(diǎn)區(qū)間方法更加符合系統(tǒng)可靠性的退化規(guī)律，精度更高，一定程度上克服了隨著狀態(tài)數(shù)的增加，計(jì)算過程過于繁瑣的不足。

3) 區(qū)間連續(xù)序列上限方法(狀態(tài)3-2)與Monte Carlo仿真方法結(jié)果相對(duì)存在一定的誤差，隨著時(shí)間的增加偏差增大，原因主要在于仿真的過程中部件的狀態(tài)數(shù)N=20相對(duì)較少，若部件的狀態(tài)數(shù)過多的話，在仿真計(jì)算系統(tǒng)的狀態(tài)數(shù)過于龐大，難以計(jì)算。

眾所周知，采用通用生成函數(shù)方法對(duì)多態(tài)系統(tǒng)可靠性進(jìn)行分析時(shí)，隨著狀態(tài)數(shù)的增加，精度增大，但計(jì)算的復(fù)雜度增加，下面將分析分系統(tǒng)1中部件狀態(tài)為4，分系統(tǒng)2中部件狀態(tài)為3時(shí)，系統(tǒng)區(qū)間可靠度和區(qū)間可靠度均值隨時(shí)間的變化關(guān)系，見圖5。

圖5　部件狀態(tài)取不同值時(shí)系統(tǒng)可靠性的變化Fig.5 Changes of MSS reliability when component’ state number is changed

由表3和圖5可以得出以下結(jié)論：

1) 狀態(tài)數(shù)(4-3)與狀態(tài)數(shù)(3-2)相比，在相同時(shí)刻，狀態(tài)數(shù)(4-3)中區(qū)間連續(xù)序列下限方法和區(qū)間連續(xù)序列中值方法的分析結(jié)果呈增大的趨勢(shì)。

2) 區(qū)間連續(xù)序列中值方法與Monte Carlo仿真方法相比結(jié)果雖存在一定誤差，但在一定程度上可以表征系統(tǒng)的可靠性水平，因此在損失一定的精度情況下，本文區(qū)間連續(xù)序列族中值方法不僅解決了傳統(tǒng)方法誤差大的問題，而且克服了Monte Carlo仿真方法仿真規(guī)模大、耗時(shí)長(zhǎng)的不足，方法簡(jiǎn)單，便于實(shí)現(xiàn)。

3) 采用區(qū)間通用生成函數(shù)方法隨著狀態(tài)數(shù)的增加，降低了系統(tǒng)的可靠性評(píng)估結(jié)果與MC仿真方法所得評(píng)估結(jié)果之間的誤差，因此在一定程度上驗(yàn)證隨著部件狀態(tài)數(shù)的增加，系統(tǒng)可靠性評(píng)估結(jié)果精度增大的結(jié)論。

表3　狀態(tài)數(shù)(4-3)時(shí)所得系統(tǒng)區(qū)間可靠性對(duì)比分析結(jié)果

5結(jié)論

本文針對(duì)考慮認(rèn)知不確定性與隨機(jī)沖擊情況的多態(tài)系統(tǒng)，提出了隨機(jī)沖擊情況下考慮認(rèn)知不確定的多態(tài)系統(tǒng)區(qū)間可靠性分析方法，給出了多態(tài)串-并聯(lián)系統(tǒng)可靠性的分析流程，通過仿真實(shí)例分析可得出以下結(jié)論：

(1) 該方法降低了傳統(tǒng)方法中部件失效閾值選取不當(dāng)對(duì)系統(tǒng)可靠度評(píng)估影響，提高了系統(tǒng)可靠性評(píng)估精度；

(2) 該方法在狀態(tài)數(shù)相同的情況下，與傳統(tǒng)點(diǎn)區(qū)間通用生成函數(shù)法相比，不僅精度較高，且一定程度克服了其隨著狀態(tài)數(shù)的增加，計(jì)算過程過于繁瑣的不足；

(3) 該方法與Monte Carlo仿真方法相比，克服了其仿真規(guī)模大、耗時(shí)長(zhǎng)的不足。

參 考 文 獻(xiàn)

[1] Nakagawa T. Shock and damage models in reliability theory[M]. Berlin:Springer, 2006.

[2] Liu Y, Huang H, Pham H. Reliability evaluation of systems with degradation and random shocks, 2008[C]//Reliability and Maintainability Symposium.Las Vegas,NV:RAMS,2008.

[3] Wang Ya-ping. Multi-objective impact for dependent competing system with multiple degradation progress and random shock[D]. New Brunswick,NJ:The State University of New Jersey, 2012.

[4] Lisnianski A, Levitin G. Multi-state system reliability: assessment, optimization and applications[M]. Singapore: World Scientific, 2003.

[5] Barton R M, Damon W. Reliability in a multi-state system[C]//Procedings of the 6th Annual Southeastern Symposium on Systems Theory. Baton Rouge,LA，1974.

[6] Barton R E, Wu A S. Coherent systems with multi-state components[J]. Mathematica of Operations Research, 1978(3):275-281.

[7] Li Y, Zio E. A multi-state model for the reliability assessment of a distributed generation system via universal generating function[J]. Reliability Engineering and System Safety, 2012(106):28-36.

[8] Billinton R Y G. Adequacy assessment of composite power generation model for reliability evaluation.[J]. International Journal of Reliability and Safety, 2008,2(1):79-98.

[9] 史新紅, 齊先軍, 王治國(guó). 基于UGF的發(fā)電系統(tǒng)區(qū)間可靠性評(píng)估及其仿射算法改進(jìn)[J]. 合肥工業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版, 2014,37(3):286-291.

SHI Xin-hong，QI Xian-jun，WANG Zhi-guo. Interval reliability estimation of power generating system based on UGF method and its modification by using affine arithmetic[J]. Journal of Hefei University of Technology:Natural Science, 2014,37(3):286-291.

[10] 尚彥龍, 蔡琦, 趙新文,等. 基于UGF和Semi-Markov方法的反應(yīng)堆泵機(jī)組多狀態(tài)可靠性分析[J]. 核動(dòng)力工程, 2012(1): 117-123.

SHANG Yan-long, CAI Qi, ZHAO Xin-wen, et al.Multi-state reliability for pump group in nuclear power system based on UGF and semi-markov process[J]. Nuclear Power Engineering,2012(1):117-123.

[11] Natvig B. Multistate systems reliability theory with applications[M]. Hoboken, NJ, USA: Wiley,2010.

[12] Levitin G. The universal generating function in reliability analysis and optimization[M]. London: Springer, 2005.

[13] Ding Y, Zuo M J, Lisnianski A,et al.Fuzzy multi-state systems: general definitions, and performance assessment[J]. IEEE Transactions on Reliability, 2008,57(4):589-594.

[14] Ding Y, Lisnianski A. Fuzzy universal generating functions for multi-state system reliability assessment[J]. Fuzzy Sets and Systems, 2008，159:307-324.

[15] 鄢民強(qiáng), 楊波, 王展. 不完全覆蓋的模糊多狀態(tài)系統(tǒng)可靠性計(jì)算方法[J]. 西安交通大學(xué)學(xué)報(bào), 2011,45(10):109-114.

YAN Min-qiang, YANG Bo, WANG Zhan. Reliability assessment for multi-state system subject to imperfect fault coverage[J].Journal of Xi’an Jiaotong University, 2011,45(10): 109-114.

[16] Li Chun-yang,Chen Xun, Yi Xiao-shan,et al. Interval-valued reliability analysis of multi-state systems[J]. IEEE Transactions on Reliability, 2011,60(1):323-330.

[17] Li W, Pham H. Reliability modeling of multi-state degraded systems with multi-competing failures and random shocks[J]. IEEE Transactions on Reliability, 2005,54(2):297-303.

[18] 邱志平. 非概率集合理論凸方法及其應(yīng)用[M]. 北京: 國(guó)防工業(yè)出版社, 2005.

Estimation of interval-valued reliability of multi-state system in consideration of epistemic uncertainty under random shock

PAN Gang, SHANG Chao-xuan, LIANG Yu-ying, CAI Jin-yan, MENG Ya-feng

(Department of Electronic and Optic Engineering, Ordnance Engineering College, Shijiazhuang 050003, China)

Abstract:Components may suffer from random shocks in some environments or other conditions. The performance degradation of components consists of normal performance degradation and random shocks. Since it is hard to obtain adequate performance data of high-reliability components within a short time, there are epistemic uncertainties on components and the system reliability cannot be accurately estimated. For the purpose of accurate estimation of system reliability, the distribution parameters of components’ performance deterioration caused by random shocks were assumed as interval variables, the components’ performance distribution model based on the interval variables was built, the definition of interval-continuous sequences of components’ state performances was presented and a method to calculate the interval-valued state probability were provided. The traditionally universal generating function method was improved, the interval-valued universal generating function and its algorithm were defined, a method to assess multi-state system reliability in consideration of epistemic uncertainty under random shocks was proposed, and the verification and illustration were conducted by using simulations. The method not only overcomes the inaccuracy of the reliability analysis model, but also has strong versatility and engineering application value.

Key words:random shock; epistemic uncertainty; interval-continuous sequence; interval-valued universal generating function

基金項(xiàng)目：國(guó)家自然科學(xué)基金(61372039;61271153)

收稿日期：2015-04-02修改稿收到日期：2015-05-31

通信作者尚朝軒 男，博士，教授，博士生導(dǎo)師，1964年4月生

中圖分類號(hào):N945;TB114.3

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2016.10.005

第一作者 潘剛 男，博士生，1987年7月生