王 志,李星新,王 成,于永利
(陸軍工程大學(xué)石家莊校區(qū),石家莊 050003)
在裝備試驗過程中,為了節(jié)省試驗經(jīng)費和避免重復(fù)工作、拖延進度,裝備測試性試驗一般與維修性試驗同時進行,其樣本的分配和模擬故障策略與維修性試驗中維修作業(yè)樣本的分配與故障模擬相同[1]。這樣的試驗方法雖然確保了測試性驗證試驗故障樣本的真實性,但對于裝備的部分組成單元來說,在有限的試驗期間內(nèi)出現(xiàn)的自然故障數(shù)量較少,甚至無故障?;谠摌颖痉峙浞桨福@些組成單元分配的故障樣本量也常為0,這樣使得測試性驗證不充分、不全面,導(dǎo)致試驗評估結(jié)果存在偏差。而且當(dāng)驗證樣本量能對所有故障模式多次覆蓋時,對維修性驗證來說,即使是模擬(注入)同一故障,其故障修復(fù)時間也不一定相同,而對測試性驗證來說,同一故障模式多次注入意義不大,且容易出現(xiàn)故障率高的故障樣本被重復(fù)選取而故障率低但其他指標高的樣本沒有選中,最終導(dǎo)致樣本結(jié)構(gòu)不合理。
針對上述問題,文獻[2-8]分別對故障樣本分配所考慮的影響因素進行了擴充,包括故障率、危害度、嚴酷度、擴散度、被檢測難度、故障影響、平均修復(fù)時間、故障隔離時間等,并對提高部分因素的計算精度進行了研究。故障樣本的分配大多采用基于數(shù)理統(tǒng)計中分層抽樣的方法[9],主要考慮故障的某一個屬性或多個屬性綜合加權(quán),然后在已知故障樣本總量的情況下,依據(jù)量化屬性值按比例進行抽樣。目前,故障樣本分配問題主要側(cè)重于研究如何更加全面地考慮影響分配比例計算的因子、多因子的加權(quán)策略以及如何提高這些因子的計算精度,故障樣本分配的合理性有了明顯提高。例如,文獻[10]提出了一種故障樣本綜合加權(quán)分配方法,提高了分配結(jié)果的合理性。但是,該方法采用基于專家經(jīng)驗加權(quán)的方法計算指標權(quán)重,權(quán)值分配的主觀性較強,故障樣本分配結(jié)果的置信度不高;文獻[11]提出了一種基于故障屬性的分配方案,比較全面地考慮了各種故障屬性。但是,該方法關(guān)于嚴酷度的計算方法考慮情況比較簡單,計算結(jié)果精度不高,并且忽略了指標權(quán)重的計算;這些方法大都以故障檢測率指標的驗證為例,進行了案例說明。但其適用性還有待進一步證明,文獻驗證案例中涉及的單元故障模式種類數(shù)都比較少,且比例抽樣產(chǎn)生的小數(shù)取整規(guī)則沒有合理的解決。本文主要考慮故障檢測率(FDR,fault detection rate)的驗證,通過對標準中驗證方案進行梳理,對上述相關(guān)文獻研究,針對FDR驗證試驗的故障樣本分配策略問題,區(qū)分不同情況,提出了綜合考慮嚴酷度的以故障率為主體開展故障模式選取的方法,旨在為測試性驗證故障樣本的分配工作提供合理、可行的實施方法。
GJB8895將測試性驗證試驗定義為,為確定產(chǎn)品是否達到規(guī)定的測試性要求,由訂購方認可的試驗機構(gòu)按選定的驗證試驗方案,進行故障抽樣并在產(chǎn)品實物或試驗件上開展的故障注入或模擬試驗。故障檢測率作為測試性最重要的定量指標之一,其定義為用規(guī)定的方法正確檢測到的故障數(shù)與故障總數(shù)之比,常用百分數(shù)表示。測試性指標體系中,故障檢測率、故障隔離率、虛警率的抽樣檢驗都是計數(shù)抽樣檢驗。以故障檢測率的驗證試驗為例,其檢驗?zāi)P涂捎脠D1表示[12]。
圖1 故障檢測率抽樣檢驗?zāi)P?/p>
取c為合格判定數(shù),當(dāng)檢測失敗數(shù)r≤c認為檢測試驗通過,反之,不通過。
其中,總體為N次故障檢測結(jié)果構(gòu)成的集合,q為總體的故障檢測率。需要注意的是,測試性試驗中總體容量N不明確,由于認知的有限以及故障發(fā)生的不確定性,測試性試驗實施時無法獲得故障檢測所構(gòu)成的全集。相比于樣本容量n,認為N?n。q可以假定為一個定值,只能通過樣本對其進行統(tǒng)計推斷。做出正確推斷的前提是:所選取的n個樣本具有代表性和獨立性,即要求每一個個體都有同等機會被選入樣本,而且被選樣本的取值不受其他樣本取值的影響。
故障樣本分配策略即研究如何把抽樣檢驗方案確定的故障樣本總量n分配到裝備具體的各個單元、不同的故障模式中,以便于開展故障注入實驗。研究樣本的分配問題,必須先考慮已知信息,根據(jù)不同驗證方案的應(yīng)用條件限制,進行驗證方案選取。此外,雖然各個驗證方案都會對樣本量n提出要求,但是如果只從確定的故障樣本量n出發(fā)進行樣本量分配問題的研究,而不明確指出所采取的驗證方案,就會造成忽略n值需要調(diào)整時(如小數(shù)取整)判別數(shù)c主導(dǎo)變化的情況,如文獻[13]樣本量需要調(diào)整時簡單的將試驗樣本量43(故障判據(jù)為5)調(diào)整為44,但是若前期是由表1中方案3確定的樣本量n,這個調(diào)整至少變?yōu)?8(不改變雙方風(fēng)險相等的前提),對應(yīng)故障判據(jù)為6。即沒有考慮到所做調(diào)整與驗證方案的匹配性。因此,必須先進行驗證方案研究。
通過對文獻[14-15]中提及的各種測試性驗證方案(故障檢測率均適用)進行梳理,按照其樣本量的確定方法進行分類,總結(jié)所涉及標準中采取的樣本分配策略[16-18]。
根據(jù)前面的抽樣檢驗?zāi)P?,試驗中注入一個故障能被檢測(隔離)的概率為q,前提是這個故障樣本是從N中隨機抽取的,對于一個已知故障的注入來說,注入后被檢測與否是確定的,即概率為1或0(不考慮有時成功有時失敗情況)。因此,標準中提及的樣本分配方案都強調(diào)了基于故障率(主體)的分配方案:故障率高的組成單元(故障模式頻數(shù)比大的故障模式)在故障檢測總體樣本N構(gòu)成中會提供較大比重的樣本,其被隨機抽取的概率相應(yīng)也高,即這些影響分配比例計算的因子(如故障率、工作時間、單元數(shù)量)的考慮,更全面的反映了真實的構(gòu)成數(shù)量比,有助于對實際FDR做出推斷。但是,當(dāng)驗證樣本量要求能對所有故障模式多次覆蓋時,按故障率為主體來制定故障樣本分配方案,考慮因素比較單一,導(dǎo)致故障率高的故障樣本重復(fù)選取而故障率相對發(fā)生頻率小但嚴酷度或其他指標高的故障模式,分配樣本數(shù)量常為0,驗證效果不夠令人滿意。標準雖然提及了分配集中不合理問題,但并未提出統(tǒng)一的分配調(diào)整方案。因此,在不破壞抽樣檢驗樣本抽取代表性獨立性要求前提下,合理考慮其他影響因素,并明確調(diào)整方案,是十分有必要的。
根據(jù)相關(guān)輸入信息,確定驗證方案,按照方案確定樣本量n,明確其可調(diào)整的值。首先進行故障樣本量分配,完成單元樣本數(shù)量分配后再進行故障模式選取。根據(jù)分配結(jié)果,按照故障模式種類數(shù)、驗證樣本量及預(yù)選樣本量(取驗證樣本量的3或4倍)數(shù)量關(guān)系區(qū)分不同情況。嚴酷度反映的是故障發(fā)生后對裝備的影響程度,訂購方重點關(guān)注嚴酷度高的故障模式的檢測[19]。因此需要綜合考慮嚴酷度進行樣本選取。記故障模式種類數(shù)為DZ、預(yù)選樣本數(shù)量為DY、驗證樣本數(shù)量為DP,根據(jù)三者的數(shù)量關(guān)系區(qū)分如下情況:
1)DZ≥DY時(d1=DZ-DY),按不同故障模式的故障率占比抽取預(yù)選樣本數(shù)量,計算結(jié)果大于等于1的按四舍五入進行取整。根據(jù)故障模式嚴酷度,對于分配數(shù)量不足1的,按嚴酷度排序從低到高依次刪除d1個。其余計算結(jié)果不足1的調(diào)整為1,需要時在剩余故障模式嚴酷度排序低于自身的且數(shù)量大于等于2的故障模式集合中取嚴酷度最低的做減法(數(shù)量不少于1),至分配完畢。隨后,不同個數(shù)不同故障模式種類組成預(yù)選樣本集,進行簡單隨機抽樣,遵循每個樣本個體只使用一次原則抽樣,當(dāng)隨機抽樣選中一個故障模式時,該種故障模式被選中概率相應(yīng)變化(變小),其它樣本被選中概率不變。
2)DY>DZ≥DP時(d2=DY-DZ),首先每個故障模式種類擁有1個樣本,d2按照故障率進行分配,分配后按嚴酷度進行排序,兩端向中間進行兩兩配對,嚴酷度高的向上取整,低的向下取整,奇數(shù)情況中間的按照補差策略取整。隨后,由不同個數(shù)不同故障模式種類組成預(yù)選樣本集,抽樣方法同上。
3)DP>DZ時,(d3=DP-DZ),首先每個故障模式種類擁有1個樣本,d3處理方法同d2。即可完成故障模式選取。
針對FDR檢驗問題,考慮多因子加權(quán)后按比例計算的方法無法合理的處理故障模式種類數(shù)接近所需故障樣本數(shù)量的情況。例如某單元有10種故障模式需要選8個故障樣本,若按多因子加權(quán)比例進行計算,這種情況下由于取整問題實際上人為破壞了所考慮因素規(guī)定的比例。因此,考慮故障樣本的充分性要求[20],當(dāng)出現(xiàn)情況2或3時,預(yù)選樣本集中每種故障模式都有至少1個。當(dāng)出現(xiàn)情況1時,刪除故障率低且嚴酷度低的故障模式后使得預(yù)選樣本集中每種故障模式都有1個。考慮嚴酷度排序采用不同的取整策略,通過借助預(yù)選樣本集,由數(shù)量占比體現(xiàn)被選中概率,保證了故障率差距,而采用被選中概率動態(tài)變化(單獨減小),合理增加了嚴酷度高樣本被選中概率。進行隨機抽樣達到驗證個數(shù)即可,不用再考慮取整問題。
上述方法中的核心是解決故障模式的嚴酷度排序問題。文獻[19]關(guān)于嚴酷度評分將其分為4個等級的方法,不能很好地解決更多故障模式時的一致排序問題。本文采用模糊證據(jù)推理和灰色關(guān)聯(lián)理論對嚴酷度進行評價并排序。依據(jù)標準提及的打分規(guī)則,將嚴酷度評價的考慮因素細分為人員危害性(P)、財產(chǎn)危害性(W)、任務(wù)危害性(T)三項。采用5個模糊評價等級進行評價,即H={H11,H22,H33,H44,H55}={很低,低,中等,高,很高}。假設(shè)單個模糊評價等級之間相互獨立且只有相鄰等級之間存在相交情況。利用梯形分布來表示5個評價等級[21],如圖2所示。
圖2 梯形分布表示評價等級的模糊隸屬函數(shù)
利用模糊證據(jù)推理置信結(jié)構(gòu)可以使評估人員靈活的對故障模式的三項危害性進行評價:如任務(wù)危害性評價為“高”,對應(yīng){(H44,1.0)}這種置信結(jié)構(gòu)。{(H15,1.0)}表示對應(yīng)評估等級可以是從“很低”到“很高”中的任意一個。設(shè)L名專家(E1,…,El)對K個故障模式(FM1,…,FMk)關(guān)于N個評價因素項(D1,…,Dn)進行評價。
i=1,…,5;j=1,…,5;k=1,…,K;n=1,…,N
(1)
K個故障模式關(guān)于N個評價因素的綜合置信結(jié)構(gòu)構(gòu)成模糊置信決策矩陣:
利用文獻[21-22]提及方法進行去模糊化,即:
i=1,…,5;j=1,…,5
(2)
hij為Hij去模糊化后的明確值。此時計算故障模式FMk關(guān)于評價因素Dn的總置信度為具體值。
k=1,…,K;n=1,…,N
(3)
模糊置信決策矩陣也可完成去模糊化。
采用灰色關(guān)聯(lián)理論確定故障模式ESR排序。建立由所有評價因素的理想水平確定的參考矩陣X0,計算灰色關(guān)聯(lián)度R。有下式:
(4)
式中λn為綜合評價因素權(quán)重。r{X0(n),Xk(n)}為所有故障模式的各個評價因素的灰色關(guān)聯(lián)系數(shù)。
r{X0(n),Xk(n)}=
(5)
α稱為分辨系數(shù),常取0.5。不改變ESR排序。若理想水平取0值,則X0為K*N階的零矩陣。計算灰色關(guān)聯(lián)度值越小,ESR排序值越小。
表1 專家評價信息模糊置信結(jié)構(gòu)
表2 專家評價信息綜合置信結(jié)構(gòu)
表3 模糊評價等級去模糊化
X1(3)}=0.540 3,由公式(4)求得R1=0.570 6。同理計算得到各故障模式的關(guān)聯(lián)度見表4第五列。最終得到關(guān)聯(lián)度排序為FM4 表4 ESR總評價信息 故障率數(shù)據(jù)及嚴酷度排序見表5。按照本文做法,需要先選取24個樣本個體構(gòu)成預(yù)選隨機抽樣樣本集,屬于情況2。d2=13,將其按照故障率分配給每種故障模式,根據(jù)其嚴酷度等級排序進行取整調(diào)整,得到24個預(yù)選樣本集,進行隨機抽樣抽選8個驗證樣本即可,一次抽樣結(jié)果見表5。與按照故障率分配,及文獻[4]考慮故障率、危害度、擴散度、環(huán)境因子、被檢測難度的綜合加權(quán)分配方案結(jié)果進行對比。 表5 幾種分配方案結(jié)果對比 由分配結(jié)果可知,按故障率的分配方案不會兼顧嚴酷度等級高的故障模式,綜合加權(quán)比例分配方法結(jié)果受因素權(quán)重影響較大,其取整策略一定程度上破壞了實際考慮的加權(quán)比例。本文方法避免了按比例分配樣本時的取整問題,合理改善了僅僅考慮故障率的分配方案造成的樣本分配集中不合理的狀況。 本文結(jié)合測試性驗證試驗故障樣本分配的工程實踐特點,提出了一種用于FDR驗證中故障樣本分配工作的故障模式選取方法,提出了基于模糊證據(jù)推理的故障模式嚴酷度排序解決策略。該方法借助預(yù)選樣本集隨機抽樣、考慮嚴酷度的取整策略,以及動態(tài)概率調(diào)整,有效改善了故障率主體分配方案進行故障模式選取時樣本分配集中不合理的狀況。針對現(xiàn)行一些故障樣本分配方法的不足之處,進行了合理改進。經(jīng)過實例驗證,所提方法對于科學(xué)、有效地開展裝備測試性驗證試驗工作具有一定參考價值。4 結(jié)束語