翟禹堯， 史賢俊， 韓露， 呂佳朋

(海軍航空大學(xué) 岸防兵學(xué)院， 山東 煙臺(tái) 264001)

0 引言

測(cè)試模型的優(yōu)劣決定了測(cè)試性設(shè)計(jì)水平的高低，當(dāng)前應(yīng)用最廣泛的模型當(dāng)屬多信號(hào)模型，信息流模型緊隨其后[1-2]。信息流模型主要考慮故障與測(cè)試之間的關(guān)系，忽視系統(tǒng)功能，致使其描述能力有限，阻礙了提升的空間。多信號(hào)模型解決了上述模型的一些不足，簡(jiǎn)化了故障模式，采用信號(hào)描述故障和測(cè)試之間關(guān)系有效降低了建模難度。陳春良等[3]建立了某坦克火控系統(tǒng)的多信號(hào)模型，并對(duì)其進(jìn)行測(cè)試性分析，提高了火控系統(tǒng)的測(cè)試性水平?？琢顚挼萚4]采用多信號(hào)模型對(duì)衛(wèi)星電源進(jìn)行了測(cè)試性建模與分析，提出了基于多信號(hào)模型的實(shí)時(shí)故障診斷技術(shù)。孫智等[5]采用分層多信號(hào)流圖對(duì)飛機(jī)空調(diào)系統(tǒng)進(jìn)行了故障診斷，取得了不錯(cuò)的效果。張曄等[6]對(duì)雷達(dá)機(jī)內(nèi)測(cè)試(BIT)的測(cè)試能力進(jìn)行了評(píng)估，并建立了多信號(hào)模型。

多信號(hào)模型雖然解決了信息流模型的不足，但其仍有一些弊端。例如：系統(tǒng)故障與測(cè)試之間的關(guān)系僅用0和1描述，其內(nèi)在聯(lián)系不能完整地表達(dá)出來(lái)，并且面向機(jī)電、控制系統(tǒng)等復(fù)雜系統(tǒng)不再適用；忽略大量有用信息，如漏檢和虛警的影響、延時(shí)性以及故障間的邏輯關(guān)系，導(dǎo)致模型不夠完備，所得測(cè)試性分析結(jié)果不夠全面準(zhǔn)確。不少學(xué)者針對(duì)上述不足做了以下研究，楊鵬[7]、陳希祥等[8]對(duì)模型進(jìn)行了改進(jìn)研究，添加了一定的不確定性知識(shí)。陳希祥等[9]還提出基于本體的描述模型。徐星光等[10]將測(cè)試性結(jié)構(gòu)模型和測(cè)試性貝葉斯網(wǎng)絡(luò)模型相結(jié)合，融合了先驗(yàn)知識(shí)和實(shí)際試驗(yàn)結(jié)果，利用不確定信息提高了裝備測(cè)試性水平。

綜合看來(lái)，不少專家對(duì)當(dāng)前模型進(jìn)行了一定改進(jìn)及調(diào)整，使模型具備了表達(dá)不確定信息等功能，但這些方法沒(méi)有從根本上解決測(cè)試性建模問(wèn)題，如測(cè)試時(shí)延、故障間的邏輯關(guān)系以及故障嚴(yán)酷度等信息。Petri網(wǎng)描述系統(tǒng)狀態(tài)的同時(shí)還表現(xiàn)其行為，由于自身特點(diǎn)使其在很多領(lǐng)域廣泛應(yīng)用[11]。故障Petri網(wǎng)在可靠性、安全性方面應(yīng)用較多[12]，本文不做研究。隨機(jī)Petri網(wǎng)是Petri網(wǎng)的進(jìn)一步擴(kuò)展[13]，可以描述系統(tǒng)狀態(tài)及其動(dòng)態(tài)變化。王瑤[14]對(duì)可靠性、安全性、維修性、測(cè)試性工作開展了深入研究，將系統(tǒng)正常狀態(tài)、故障狀態(tài)、功能失效狀態(tài)引入隨機(jī)Petri網(wǎng)中，起到了良好的效果。文獻(xiàn)[15]采用廣義隨機(jī)有色Petri網(wǎng)(CGSPN)建立系統(tǒng)層次測(cè)試性模型，對(duì)系統(tǒng)層次劃分進(jìn)行了詳細(xì)闡述，添加了編碼方式和嚴(yán)酷度信息，增強(qiáng)了模型的可視化并豐富了模型的內(nèi)容，但是沒(méi)有考慮測(cè)試和故障之間的不確定性問(wèn)題。

本文在文獻(xiàn)[15]基礎(chǔ)上對(duì)故障間的復(fù)雜性進(jìn)行分析，完善測(cè)試性模型的構(gòu)建；利用專家評(píng)分的先驗(yàn)知識(shí)與小樣本測(cè)試數(shù)據(jù)的后驗(yàn)知識(shí)對(duì)模型中的不確定信息進(jìn)行處理，引入三角模糊數(shù)算法獲取專家知識(shí)，計(jì)算關(guān)于故障檢測(cè)率的先驗(yàn)條件概率分布以歷史測(cè)試數(shù)據(jù)作為后驗(yàn)條件，由貝葉斯公式更新故障檢測(cè)率的概率分布，解決了數(shù)據(jù)容量小的問(wèn)題。

1 基于CGSPN的測(cè)試性模型

1.1 Petri網(wǎng)基本理論

定義1滿足下列條件的4元組PN=(P,T;F,M0)構(gòu)成Petri網(wǎng)，其中：P表示庫(kù)所，T表示變遷，F(xiàn)為弧權(quán)函數(shù)，M0為初始標(biāo)識(shí)。

1)N為基本網(wǎng)。

2)M為系統(tǒng)狀態(tài)標(biāo)識(shí)，M:P→Z，Z為自然數(shù)。

3)變遷發(fā)生規(guī)則：

①變遷t∈T，若?P∈t*，M(P)≥1，則t可被觸發(fā)，記作M[t>，其中t*為變遷t的后置集，M(P)為庫(kù)所標(biāo)識(shí)。

②M′為觸發(fā)后標(biāo)識(shí)，則

記作M[t>M′. 其中p為某種故障模式的故障檢測(cè)率；*t為變遷t的前置集；標(biāo)識(shí)M可以用1個(gè)非負(fù)整數(shù)的m維向量表示，向量元素滿足M(i)=M(Pi)，Pi為一個(gè)庫(kù)所，i=1,2,…,m.

可達(dá)性(見(jiàn)定義2)是Petri網(wǎng)建模過(guò)程中分析的重要方式。

定義2Petri網(wǎng)PN=(P,T；F,M0)，系統(tǒng)標(biāo)識(shí)?M1,M2,…,Mi,…,Mk，使得?1≤i≤k，對(duì)于任意變遷?ti∈T:[ti>Mi+1，則稱變遷序列σ=t1,t2,…,ti,…,tk在M1下是使能的，Mk+1從M1是可達(dá)的，記作M1[σ>Mk+1標(biāo)識(shí)。

廣義隨機(jī)Petri網(wǎng)(GSPN)由一個(gè)7元組構(gòu)成[15-16]，記為GSPN=(P,T;F,M0,I,O,λ)，其中：I和O為輸入矩陣、輸出矩陣；λ表示變遷速率或概率。

表1所示為GSPN元素的基本含義以及其符號(hào)表示方法。

1.2 故障模式在Petri網(wǎng)表達(dá)方式

1.2.1 故障模式編碼、嚴(yán)酷度等級(jí)劃分

在實(shí)際工作中，嚴(yán)酷度較高的故障一旦發(fā)生，會(huì)帶來(lái)嚴(yán)重后果，在模型中根據(jù)顏色可以清楚地看到每個(gè)庫(kù)所(故障模式)的嚴(yán)酷度，密切關(guān)注嚴(yán)酷度較高的故障并設(shè)置相應(yīng)的測(cè)試，會(huì)為后續(xù)的故障診斷提供方便，有效保護(hù)裝備系統(tǒng)。根據(jù)國(guó)家軍用標(biāo)準(zhǔn)GJB/Z 1391—2006 故障模式影響及危害性分析指南[17]，故障模式的嚴(yán)酷度分為：災(zāi)難級(jí)、嚴(yán)重級(jí)、輕度級(jí)和輕微級(jí)4類[18]，分別對(duì)應(yīng)紅色(R)、橙色(O)、黃色(Y)和藍(lán)色(B)。應(yīng)當(dāng)說(shuō)明的是，除了上述4類故障模式外，本文還定義了一種功能失效狀態(tài)，對(duì)于具有n種功能模式的元件，其功能失效狀態(tài)亦有n種。不同于上述故障模式，功能失效狀態(tài)稱為失效級(jí)，用白色(W)表示。

表1 GSPN模型的圖形化表示Tab.1 Graphical representation of GSPN

本文定義的顏色屬性用于區(qū)別模型中各故障模式的嚴(yán)酷度，為后續(xù)測(cè)試優(yōu)化選擇的研究提供方案。值得注意的是該顏色屬性與傳統(tǒng)有色Petri網(wǎng)的概念不同[19]，模型結(jié)構(gòu)不會(huì)因?yàn)橛猩珜傩缘南Ф淖儭?/p>

對(duì)復(fù)雜系統(tǒng)建模，需要考慮某個(gè)故障可能存在多個(gè)故障模式，為了降低建模難度并提高系統(tǒng)測(cè)試性，本文采用以下編碼方案[15]：編碼方案：XXX-XXX-XXX-OO.x. 其中，XXX代表系統(tǒng)、現(xiàn)場(chǎng)可更換單元(LRU)和車間可更換單元(SRU)等，可以繼續(xù)劃分到元器件；OO.x中OO為不同嚴(yán)酷度的顏色，x∈n為元件的故障模式種類。

1.2.2 故障傳播分析

系統(tǒng)故障傳播具有橫向性、縱向性、復(fù)雜性、時(shí)間性以及邏輯復(fù)雜性等。故障會(huì)逐級(jí)傳遞，從底層元件傳遞至系統(tǒng)層[14]。對(duì)裝備進(jìn)行分層處理就完成了故障傳播的橫向性、縱向性處理；考慮到導(dǎo)彈長(zhǎng)期貯存、一次性使用的特點(diǎn)，其故障模式特點(diǎn)與飛機(jī)等常規(guī)武器裝備不同，時(shí)間性特點(diǎn)并不明顯?？紤]到篇幅，邏輯復(fù)雜性在本文暫不作研究，只對(duì)故障模式的復(fù)雜性展開論述。

一般來(lái)說(shuō)，導(dǎo)彈故障模式存在相容和相斥兩種形式。相容模式在模型中表達(dá)如下：1個(gè)庫(kù)所表示一種故障模式，令牌流入某個(gè)庫(kù)所，則代表該故障發(fā)生。根據(jù)1.2.1節(jié)所述為庫(kù)所編號(hào)，并根據(jù)故障模式、影響和危害性分析(FMECA)中故障嚴(yán)酷度等級(jí)為庫(kù)所賦予顏色。此外庫(kù)所與故障模式共用1個(gè)編號(hào)，以確保庫(kù)所在模型中的唯一性。圖1(a)所示為相容模式的簡(jiǎn)單實(shí)例，第2個(gè)庫(kù)所編號(hào)為001-001-001-O.1，表示系統(tǒng)的元件編號(hào)為001-001-001，其嚴(yán)酷度等級(jí)為嚴(yán)重級(jí)。

圖1 故障模式在Petri網(wǎng)中的表現(xiàn)形式Fig.1 Failure modes in Petri net

對(duì)于相斥模式，其表示方法可分為以下兩個(gè)步驟：

1)表達(dá)方式與相容模式相同，首先賦予庫(kù)所顏色屬性以及編碼方案。

2)與相容模式不同之處在于使用了禁止弧。假設(shè)系統(tǒng)存在故障模式A和B，若A與B相斥，則在系統(tǒng)中找出B庫(kù)所的前置變遷集合T，在A庫(kù)所與變遷ti(ti∈T)之間添加禁止弧，表示A發(fā)生后B禁止發(fā)生。若B發(fā)生后A禁止發(fā)生，則在庫(kù)所B以及A庫(kù)所的前置變遷間添加禁止弧即可。圖1(b)所示為3個(gè)故障模式相斥時(shí)的Petri網(wǎng)。

本文在模型中設(shè)置正常狀態(tài)庫(kù)所，并為延時(shí)變遷賦予相應(yīng)的隨機(jī)分布，目的是使CGSPN具有表達(dá)元件由正常到故障這一隨機(jī)過(guò)程。實(shí)際工作中，用綠色表示正常狀態(tài)。若某個(gè)元件具有K種故障模式，則需要在模型中設(shè)置K個(gè)正常狀態(tài)庫(kù)所，使得正常庫(kù)所與故障模式的庫(kù)所一一對(duì)應(yīng)。圖2所示為元件4種故障模式的簡(jiǎn)單實(shí)例。從圖2中可以看出：變遷點(diǎn)火變化前，令牌存在于001-001-001.G.x(x=1,2,3,4)中呈現(xiàn)綠色，系統(tǒng)正常；滿足點(diǎn)火條件后，庫(kù)所001-001-001.R.1、001-001-001.O.1、001-001-001.Y.1、001-001-001.B.1存在令牌，系統(tǒng)出現(xiàn)故障。

圖2 變遷在CGSPN中的表達(dá)形式Fig.2 Expression of transition in CGSPN

1.3 構(gòu)建思路

在測(cè)試性模型中，某個(gè)測(cè)試點(diǎn)與測(cè)試項(xiàng)沒(méi)有明確對(duì)應(yīng)關(guān)系，本文對(duì)其不做研究。假設(shè)測(cè)試點(diǎn)與測(cè)試項(xiàng)一一對(duì)應(yīng)，故障檢測(cè)率、故障隔離率的計(jì)算與該假設(shè)無(wú)關(guān)。以該假設(shè)為前提，將測(cè)試項(xiàng)信息添加到CGSPN模型中：

1)測(cè)試項(xiàng)庫(kù)所的建立：為每個(gè)測(cè)試項(xiàng)建立1個(gè)庫(kù)所。

2)故障與測(cè)試之間的關(guān)系用有向流和瞬時(shí)變遷表示。定義測(cè)試庫(kù)所的顏色為黑色(B)，目的是與模型中的其他庫(kù)所區(qū)分。

其基本構(gòu)建思路如下：

1)建立庫(kù)所，元件對(duì)應(yīng)的庫(kù)所包括故障模式、功能失效和正常狀態(tài)3種狀態(tài)；子系統(tǒng)級(jí)別以上包括故障和正常2種狀態(tài)。

2)建立延時(shí)變遷，設(shè)置在正常庫(kù)所和故障庫(kù)所之間，完成系統(tǒng)中元件由正常狀態(tài)到故障發(fā)生的描述。

3)建立兩類瞬時(shí)變遷，第1類設(shè)置在故障庫(kù)所與故障庫(kù)所之間；第2類設(shè)置在故障庫(kù)所與測(cè)試庫(kù)所之間。

4)為起始庫(kù)所設(shè)置令牌。

1.4 基于CGSPN的測(cè)試性模型

CGSPN定義為一個(gè)9元組Σ=(P,T;F,K,M0,α,β,η,Fs)，其中：

1)P={PG,PS,PF,Pt}為庫(kù)所集合，PG、PS、PF分別對(duì)應(yīng)正常、故障和功能失效的庫(kù)所集合，Pt為測(cè)試庫(kù)所集合；

2)T={Td,Ti,Tf},Td為正常庫(kù)所到故障庫(kù)所間的延時(shí)變遷集合,Ti為故障庫(kù)所到故障庫(kù)所間的瞬時(shí)變遷集合，Tf為故障庫(kù)所到測(cè)試庫(kù)所間的瞬時(shí)變遷集合；

3)F為弧權(quán)函數(shù)；

4)K:P映射為{0,1}，是Σ中庫(kù)所容量函數(shù)；

5)M0與前文含義相同；

6)α為一個(gè)映射，表示顏色屬性，α:PG映射為{G},PS映射為{R,O,Y,B},PF映射為{W}，Pt映射為{B}；

7)β為一個(gè)映射，β:Td為延時(shí)變遷服從的分布類型及函數(shù)，表明Σ存在隨機(jī)性，如Weibull(α，β)、Exp(λ)等；

8)η為一個(gè)映射，η:Ti,Tf為瞬時(shí)變遷；

9)Fs為已觸發(fā)的變遷序列集合，系統(tǒng)未運(yùn)行時(shí)Fs=?；Fs可以避免系統(tǒng)陷入循環(huán)，并記錄令牌的傳播路徑[19]，對(duì)于模型分析至關(guān)重要。

需要注意的是，本文主要針對(duì)測(cè)試性模型開展研究，在建模中只需要考慮故障和測(cè)試之間的關(guān)系，得到相關(guān)性矩陣即可，元件由正常向故障衍化的隨機(jī)過(guò)程不做考慮，故在仿真前將PG類庫(kù)所及后置延時(shí)變遷刪除。

1.5 可達(dá)性算法

故障- 測(cè)試相關(guān)性矩陣(DM)簡(jiǎn)稱D矩陣，D的獲取是測(cè)試性設(shè)計(jì)至關(guān)重要的一步。根據(jù)定義2對(duì)故障傳播過(guò)程進(jìn)行可達(dá)性分析，獲得D矩陣的算法，簡(jiǎn)稱可達(dá)性算法，其基本原理如下：

首先將系統(tǒng)庫(kù)所的初始狀態(tài)記為0，然后根據(jù)仿真需求設(shè)置令牌于某個(gè)初始庫(kù)所中，滿足變遷條件后該令牌會(huì)進(jìn)行流動(dòng)，最后將其經(jīng)過(guò)的庫(kù)所做出相應(yīng)的標(biāo)記如下：將令牌到達(dá)的故障庫(kù)所或測(cè)試庫(kù)所標(biāo)記為1，并將該庫(kù)所作為起始庫(kù)所，讓令牌繼續(xù)傳播。其具體步驟在文獻(xiàn)[15]種給出具體描述，本文不做贅述。

2 不確定信息處理方法

在實(shí)際工程中，測(cè)試與故障存在不確定性問(wèn)題，為了描述故障與測(cè)試之間不確定性問(wèn)題，用三角模糊函數(shù)確定先驗(yàn)專家知識(shí)，與實(shí)驗(yàn)獲得的測(cè)試數(shù)據(jù)作相結(jié)合，完成故障- 測(cè)試不確定相關(guān)性矩陣的確定。

2.1 貝塔分布

貝塔分布是[0,1]區(qū)間上的連續(xù)概率分布，其與二項(xiàng)分布共軛先驗(yàn)，利用該性質(zhì)可以將專家的先驗(yàn)信息和實(shí)際的測(cè)試數(shù)據(jù)相結(jié)合。

通常認(rèn)為故障檢測(cè)率的計(jì)算模型服從二項(xiàng)分布，則故障檢測(cè)率p的先驗(yàn)分布可用Beta(a,b)表示，p的概率密度函數(shù)記為f(p;a,b)。Beta(a,b)為標(biāo)準(zhǔn)化函數(shù)，其概率密度積分等于1，a、b為β分布的組成參數(shù)。

p(t|f)為故障f由測(cè)試t檢測(cè)的先驗(yàn)概率，則有

p(t|f)=f(p;a,b)，

(1)

Beta(a,b)分布表達(dá)式為

(2)

其期望和方差為

(3)

(4)

根據(jù)(3)式和(4)式可得

(5)

2.2 基于三角模糊數(shù)的專家知識(shí)處理方法

在裝備測(cè)試性設(shè)計(jì)初期，相關(guān)測(cè)試性信息不夠充分，一些信息需要依賴專家知識(shí)進(jìn)行確定。本文采用三角模糊函數(shù)對(duì)專家知識(shí)進(jìn)行處理，以確保專家知識(shí)的權(quán)威性和可靠性。l位專家對(duì)p的評(píng)分結(jié)果為1(s1,q1,u1;θ1),2(s2,q2,u2;θ2),…，j(sj,qj,uj;θj)，…,l(sl,ql,ul;θl)。其中，sj(j=1，2，…，l)為p的下限，qj為p的可能最大值，uj為p的上限，θj為第j位專家的評(píng)分權(quán)重。令ωj為歸一化處理后的權(quán)重，

(6)

專家的平均評(píng)分模糊數(shù)為

(7)

式中：s為p的下限期望；q為p的可能最大值的期望；u為p的上限期望。

基于上述處理后，p的先驗(yàn)期望和方差為

(8)

(9)

式中：n=1,2，….

可通過(guò)(8)式、(9)式和(10)式計(jì)算β分布的組成參數(shù)a、b：

(10)

式中：V(p)為分布方差。

最后根據(jù)先驗(yàn)信息得到Dpt|f，

(11)

式中：元素ptl|fk為故障由測(cè)試檢測(cè)到的概率，tl|fk表示故障fk由測(cè)試tl檢測(cè)，k為總的故障個(gè)數(shù)。

2.3 融合測(cè)試數(shù)據(jù)的后驗(yàn)分布確定方法

結(jié)合成敗型數(shù)據(jù)(x,N)，根據(jù)貝葉斯理論及(1)式，得到p的后驗(yàn)分布為

pb(t|f)=f(p;a+x,b+N-x)，

(12)

式中：N表示測(cè)試次數(shù)；x表示通過(guò)測(cè)試的次數(shù)。推理過(guò)程可參考文獻(xiàn)[19]，則經(jīng)實(shí)際測(cè)試數(shù)據(jù)修正后的D矩陣為

(13)

矩陣中的元素可由(12)式計(jì)算得到，含義為故障可由測(cè)試檢測(cè)的后驗(yàn)概率。根據(jù)Dpbt|f開展測(cè)試性評(píng)價(jià)工作，用FDRg表示fg的故障率，用Ag表示測(cè)試檢測(cè)故障的事件，假設(shè)各個(gè)測(cè)試檢測(cè)同一故障的事件是相互獨(dú)立的，則有

FDRg=pb(t1+t2+…+tl|fg)=
pb[(t1|fg)+(t2|fg)+…+(tl|fg)]=
pb(A1+A2+…+Ag)，

(14)

式中：g=1,2,…,k，從而有

(15)

3 實(shí)例分析

本文選擇文獻(xiàn)[10]的例子進(jìn)行建模分析，導(dǎo)彈測(cè)試性指標(biāo)一般要求故障檢測(cè)率FDR≥0.95，故障隔離率FIR≥0.92. 圖3所示為某型飛航導(dǎo)彈系統(tǒng)結(jié)構(gòu)框圖，系統(tǒng)分為系統(tǒng)級(jí)、分系統(tǒng)級(jí)和設(shè)備器件3個(gè)層級(jí)。表2所示為故障模式及其編碼，包含8個(gè)故障模式，故障率單位為10-6. 測(cè)試有6個(gè)，分別為電氣狀態(tài)字檢查t1、電子控制單元(ECU)狀態(tài)字檢查t2、繼電器動(dòng)作遙測(cè)指示t3、發(fā)動(dòng)機(jī)IO輸出線路電壓測(cè)量t4、發(fā)動(dòng)機(jī)點(diǎn)火信號(hào)的輸出頻率采集t5、發(fā)動(dòng)機(jī)轉(zhuǎn)速表狀態(tài)指示t6. 文獻(xiàn)[10]沒(méi)有提到故障模式相斥，經(jīng)查閱資料和相關(guān)FMECA，確定發(fā)動(dòng)機(jī)點(diǎn)火繼電器故障與發(fā)動(dòng)機(jī)點(diǎn)火裝置故障相斥。故障之間關(guān)系有利于測(cè)試的設(shè)置，有效減少測(cè)試個(gè)數(shù)。例如兩種關(guān)鍵故障需要同時(shí)設(shè)置2個(gè)測(cè)試進(jìn)行檢測(cè)，但是兩種故障模式相斥，不會(huì)同時(shí)發(fā)生，因此只需設(shè)置1個(gè)測(cè)試即可，從而為多故障診斷奠定了基礎(chǔ)。

圖3 某型飛航導(dǎo)彈的系統(tǒng)結(jié)構(gòu)框圖Fig.3 Block diagram of a winged missile system

表2 故障模式及其編碼Tab.2 Failure modes and codes

考慮到故障傳播縱向性，高層次系統(tǒng)的故障一般源于低層次子系統(tǒng)故障，而文獻(xiàn)[10]中實(shí)例模型的傳播方向相反，與實(shí)際不符。采用本文所提模型重新建立該系統(tǒng)的測(cè)試性模型，如圖4所示。圖4中，Pt表示導(dǎo)彈由其子系統(tǒng)的故障模式而引起的功能失效。由圖4可見(jiàn)，相比于文獻(xiàn)[10]以及多信號(hào)模型，本文模型不僅可以描述故障傳播過(guò)程，還可以描述故障的嚴(yán)酷度以及故障模式的復(fù)雜性等屬性。

圖4 某飛航導(dǎo)彈分層CGSPN模型Fig.4 Layered CGSPN model of a winged missile

根據(jù)圖4得到故障相關(guān)性矩陣如表3所示，選定總體、電氣和控制專業(yè)領(lǐng)域?qū)＜以u(píng)分，采用三角模糊數(shù)算法統(tǒng)計(jì)評(píng)分結(jié)果，得出各個(gè)故障模式可被檢測(cè)的先驗(yàn)概率為f(p;a,b)，融合先驗(yàn)信息的概率矩陣為Dpt|f：

限于篇幅，飛行試驗(yàn)后積累的歷史測(cè)試數(shù)據(jù)不一一列舉，根據(jù)表3以及(12)式將測(cè)試歷史數(shù)據(jù)與先驗(yàn)條件概率融合,得到融合后驗(yàn)信息矩陣Dpbt|f:

表3 相關(guān)性矩陣Tab.3 Correlation matrix

選取矩陣元素期望，得到用于測(cè)試性評(píng)價(jià)的D矩陣：

在MATLAB軟件中編程，以FDR3為例進(jìn)行計(jì)算，以說(shuō)明(14)式的計(jì)算方式，矩陣中各個(gè)測(cè)試對(duì)相同故障檢測(cè)的事件相互獨(dú)立，根據(jù)p(AB)=p(A)p(B)和(14)式，有

FDRi=pb(t1+t2+…+tl|fi)
FDR3=pb[(t2|f3)+(t4|f3)+…+(t6|f3)]=
pb(A2+A4+A6)=pb(A2)+pb(A4)+pb(A6)-
pb(A2A4)-pb(A2A6)-pb(A4A6)+pb(A2A4A6)=
pb(A2)+pb(A4)+pb(A6)-pb(A2)pb(A4)-
pb(A2)pb(A6)-pb(A4)pb(A6)+
pb(A2)pb(A4)pb(A6)=
0.75+0.71+0.81-0.75×0.71-0.71×
0.81-0.75×0.81+0.71×0.81×0.75=0.986.

經(jīng)過(guò)計(jì)算，得到每種故障模式的故障檢測(cè)率分別為FDR1=0.995，F(xiàn)DR2=0.992，F(xiàn)DR3=0.986，F(xiàn)DR4=0.996，F(xiàn)DR5=0.972，F(xiàn)DR6=0.912，F(xiàn)DR7=0.924，F(xiàn)DR8=0.970故障模式的故障率，并根據(jù)(15)式可以得到裝備總體的故障檢測(cè)率FDR=0.968. 矩陣中沒(méi)有相同行，故障檢測(cè)即可隔離FIR=1，滿足指標(biāo)要求。

4 結(jié)論

本文針對(duì)傳統(tǒng)測(cè)試性建模定性描述故障和測(cè)試之間關(guān)系的問(wèn)題，提出了一種新的測(cè)試性建模方法。首先在模型中引入嚴(yán)酷度等級(jí)，完成故障模式的分類、編碼和著色；然后對(duì)系統(tǒng)、子系統(tǒng)以及元件在模型中的表達(dá)進(jìn)行了詳細(xì)論述，考慮裝備的復(fù)雜性，將系統(tǒng)進(jìn)行層次劃分，完成模型的構(gòu)建。得到以下主要結(jié)論:

1)相比于文獻(xiàn)[10]以及多信號(hào)模型，CGSPN在完成測(cè)試性建模的基礎(chǔ)上聚焦于故障模式的嚴(yán)酷度、故障模式的復(fù)雜關(guān)系以及故障發(fā)生的隨機(jī)過(guò)程。在模型中根據(jù)顏色可以清楚地看到每個(gè)庫(kù)所的嚴(yán)酷度，密切關(guān)注嚴(yán)酷度較高的故障并設(shè)置相應(yīng)的測(cè)試，可以有效避免危害性故障的發(fā)生，也可以為后續(xù)的測(cè)試優(yōu)化選擇做好鋪墊。

2) 引入三角模糊數(shù)算法獲取專家知識(shí)，將專家數(shù)據(jù)作為先驗(yàn)信息，與實(shí)驗(yàn)獲得的測(cè)試數(shù)據(jù)作為后驗(yàn)數(shù)據(jù)相結(jié)合，用來(lái)解決數(shù)據(jù)量少和不可靠問(wèn)題。

3)通過(guò)實(shí)例分析，得到96.8%的故障檢測(cè)率和100%的故障隔離率，滿足系統(tǒng)測(cè)試行指標(biāo)要求，驗(yàn)證了模型的有效性。

4)虛警率的研究和確定將是后續(xù)的研究方向。