戴順安，王　燁，蔡　景

(1.中國(guó)民航上海航空器適航審定中心，上海　200335； 2.南京航空航天大學(xué)，南京　211106)

?

民用飛機(jī)隱蔽故障風(fēng)險(xiǎn)的定量評(píng)估方法研究

戴順安1，王燁1，蔡景2

(1.中國(guó)民航上海航空器適航審定中心，上海200335； 2.南京航空航天大學(xué)，南京211106)

摘要：對(duì)隱蔽性故障開(kāi)展風(fēng)險(xiǎn)評(píng)估是保持民機(jī)持續(xù)適航的重要方面，針對(duì)隱蔽性故障的特點(diǎn)，結(jié)合ARP 5150的風(fēng)險(xiǎn)評(píng)估理念，開(kāi)展了民用飛機(jī)隱蔽故障風(fēng)險(xiǎn)的定量評(píng)估方法研究，采用帶截尾數(shù)據(jù)的威布爾分布極大似然估計(jì)法，利用實(shí)際使用維修數(shù)據(jù)確定了隱蔽性故障的故障規(guī)律；分析了影響隱蔽性故障風(fēng)險(xiǎn)的因素，建立了新的隱蔽故障風(fēng)險(xiǎn)的定量評(píng)估方法；最后通過(guò)實(shí)例，對(duì)隱蔽性故障的風(fēng)險(xiǎn)進(jìn)行定量評(píng)估及比較分析，結(jié)果表明了本方法的有效性和適用性。

關(guān)鍵詞：民用飛機(jī)；風(fēng)險(xiǎn)；隱蔽性故障；定量評(píng)估

本文引用格式：戴順安，王燁，蔡景.民用飛機(jī)隱蔽故障風(fēng)險(xiǎn)的定量評(píng)估方法研究[J].兵器裝備工程學(xué)報(bào)，2016(6):162-165.

Citationformat:DAIShun-an,WANGYe,CAIJing.ResearchonQuantitativeRiskAssessmentMethodforHiddenFailureofCivilAircraft[J].JournalofOrdnanceEquipmentEngineering,2016(6):162-165.

由于飛機(jī)空勤組在履行正常職責(zé)時(shí)，無(wú)法發(fā)現(xiàn)飛機(jī)的隱蔽性故障，因此隱蔽性故障對(duì)飛機(jī)而言，存在著重大的安全隱患，成為民機(jī)關(guān)注的焦點(diǎn)[1-3]。為了避免隱蔽性故障的發(fā)生，在以ATAMSG-3為指導(dǎo)的維修大綱制定中規(guī)定：“必須選擇有效和適用的維修工作，避免隱蔽故障的發(fā)生，否則必須重新設(shè)計(jì)”。對(duì)隱蔽故障的重視不僅體現(xiàn)在維修大綱的制定過(guò)程中，而且在飛機(jī)投入運(yùn)營(yíng)后，需持續(xù)評(píng)估隱蔽故障的風(fēng)險(xiǎn)，在ARP5150的附錄D1.6中對(duì)隱蔽性故障進(jìn)行了解釋，并給出了簡(jiǎn)單的定量評(píng)估隱蔽性故障風(fēng)險(xiǎn)的例子。例子中假定故障服從指數(shù)分布，而且假定每天都進(jìn)行針對(duì)隱蔽故障的“故障發(fā)現(xiàn)檢查”(FailureFindingInspection，F(xiàn)FI)。雖然該例子簡(jiǎn)單易懂，但滿足以上假定條件的隱蔽故障很少，因此，僅憑ARP5150附錄D1.6中提供的方法，難以在實(shí)際持續(xù)適航管理中開(kāi)展針對(duì)隱蔽性故障的風(fēng)險(xiǎn)定量評(píng)估。

為此，本文以ARP5150附錄D1.6中提供的方法為基礎(chǔ)，結(jié)合多年從事持續(xù)適航風(fēng)險(xiǎn)評(píng)估的經(jīng)驗(yàn)，首先利用實(shí)際使用的維修數(shù)據(jù)，確定隱蔽性故障的故障規(guī)律；其次，根據(jù)隱蔽故障的“故障發(fā)現(xiàn)”維修策略，分析影響隱蔽性故障風(fēng)險(xiǎn)的主要因素；最后，建立隱蔽性故障風(fēng)險(xiǎn)的定量評(píng)估方法，并進(jìn)行實(shí)例驗(yàn)證。

1　隱蔽性故障的分布確定

ARP5150附錄D1.6的例子中，假定隱蔽故障服從指數(shù)分布。然而，根據(jù)美國(guó)宇航局統(tǒng)計(jì)數(shù)據(jù)表明，航空設(shè)備故障率大致分為6種類型，其中指數(shù)分布只占14%。目前在我國(guó)的民用飛機(jī)可靠性管理中，威布爾分布因能替代指數(shù)分布等，具有更強(qiáng)的通用性而被廣泛采用[4-5]。根據(jù)比例風(fēng)險(xiǎn)模型(proportionalhazardsmodel,PHM)或健康管理的理念可知，相同部件在不同使用環(huán)境下表現(xiàn)出來(lái)的故障規(guī)律各不相同，甚至差異很大，因此，必須根據(jù)部件實(shí)際的使用維修數(shù)據(jù)，擬合得到部件的故障分布參數(shù)[6-10]。

在實(shí)際的使用維修中，航空公司為了減少隱蔽性故障的風(fēng)險(xiǎn)，通常在實(shí)施多次“故障發(fā)現(xiàn)檢查”工作后，如未發(fā)現(xiàn)故障，就進(jìn)行預(yù)防性更換，這其實(shí)就是體現(xiàn)了ATAMSG-3中的“綜合維修”工作的思想。隱蔽性故障的故障發(fā)現(xiàn)檢查和預(yù)防性更換的時(shí)序圖如圖1所示。

從圖1可以看出，部件有的是因失效而更換，這種失效是由檢查發(fā)現(xiàn)的，因此部件失效更換的時(shí)間不是部件真實(shí)的確切失效時(shí)間，但兩者非常接近，因此，可以近似為完全數(shù)據(jù)；而有的部件則是因預(yù)防性維修而更換的，屬于右截尾壽命數(shù)據(jù)[11]。假定收集了n組部件的更換數(shù)據(jù)，其中r組為完全數(shù)據(jù)，分別為t(1)≤t(2)≤…≤t(r)，n-r組為右截尾數(shù)據(jù)，預(yù)防性更換時(shí)間為tp。

已知兩參數(shù)威布爾分布的密度函數(shù)[12-13]為

因此，其對(duì)應(yīng)的似然函數(shù)為

對(duì)上式取對(duì)數(shù)并求導(dǎo)，再整理簡(jiǎn)化后可得：

(1)

這是兩個(gè)超越方程，需用數(shù)值方法迭代求解，以得到參數(shù)m和η。

2　隱蔽性故障風(fēng)險(xiǎn)的定量評(píng)估

2.1隱蔽性故障的風(fēng)險(xiǎn)因素分析

根據(jù)圖1分析可知，影響隱蔽性故障的發(fā)生概率來(lái)自兩個(gè)方面，已使用時(shí)間T，即隨著使用時(shí)間的推移，部件的故障率水平就越來(lái)越高，因此故障發(fā)生概率越來(lái)越高；如果檢查間隔τ越大，那么在一個(gè)檢查周期內(nèi)發(fā)生故障的概率就越大。但是對(duì)于隱蔽性故障而言，故障的發(fā)生概率不等于安全風(fēng)險(xiǎn)，安全風(fēng)險(xiǎn)是隱蔽性故障概率與相關(guān)故障概率的結(jié)合。關(guān)于這一點(diǎn)，在ATAMSG-3中進(jìn)行了詳細(xì)描述：隱蔽性故障是否對(duì)安全造成影響的判據(jù)是“一個(gè)隱蔽的功能故障和另一個(gè)相關(guān)的系統(tǒng)或者備用功能的故障之綜合對(duì)使用安全性有直接有害的影響嗎？”，此判據(jù)告訴我們，隱蔽性故障本身不會(huì)對(duì)安全造成影響，只有和其他故障結(jié)合了才會(huì)造成安全問(wèn)題。譬如：燃油泄漏是一個(gè)隱蔽性故障，但該故障只有和點(diǎn)火源的失效同時(shí)發(fā)生才能引起火災(zāi)，燃油泄漏的時(shí)間越長(zhǎng)引起火災(zāi)的概率自然就越高。這也就是ARP5150附錄D1.6中所描述的，隱蔽性故障的暴露時(shí)間(ExposurePeriod)越長(zhǎng)，對(duì)應(yīng)的風(fēng)險(xiǎn)就越大。

因此，對(duì)于給定的部件(故障服從參數(shù)為m和η的威布爾分布)，已使用時(shí)間T、檢查間隔τ、故障暴露時(shí)間和相關(guān)故障(即和隱蔽性故障結(jié)合對(duì)安全有害的故障)的失效率λ，共同決定了安全風(fēng)險(xiǎn)。為此，隱蔽性故障的安全風(fēng)險(xiǎn)可以定義為關(guān)于m,η,T,τ,λ的函數(shù)H(m,η,T,τ,λ)。

2.2隱蔽性故障的風(fēng)險(xiǎn)評(píng)估

根據(jù)2.1的分析，結(jié)合圖1的描述，可知H(m,η,T,τ,λ)由兩大部分構(gòu)成：隱蔽性故障的發(fā)生概率P1，在隱蔽性故障暴露時(shí)間內(nèi)相關(guān)故障的發(fā)生概率P2。

1) 蔽性故障的發(fā)生概率。P1是一個(gè)條件概率，即在時(shí)間點(diǎn)T之前部件不發(fā)生故障的條件下，在[T,T+τ]內(nèi)發(fā)生故障的概率。因此

2) 在隱蔽性故障暴露期間內(nèi)相關(guān)故障的發(fā)生概率P2。在[T,T+τ]內(nèi)的t時(shí)刻發(fā)生故障，那么隱蔽性故障的暴露時(shí)間為T+τ-t，對(duì)應(yīng)的隱蔽性故障暴露時(shí)間的期望為

再結(jié)合相關(guān)故障的失效率λ，因此

綜合以上的分析可得，隱蔽性故障的安全風(fēng)險(xiǎn)H(m,η,T,τ,λ)為

(2)

3　實(shí)例分析

某機(jī)型燃料管泄漏屬于隱蔽功能系統(tǒng)安全影響類，表1給出了燃料管的漏油失效更換和預(yù)防性更換的壽命記錄，單位為飛行小時(shí)，其中帶“*”的為預(yù)防性更換數(shù)據(jù)。

表1　燃料管的壽命數(shù)據(jù)

采用式(2)，對(duì)表1中的壽命數(shù)據(jù)進(jìn)行威布爾參數(shù)估計(jì)，得到m=1.89，η=15 114，如將預(yù)防性更換數(shù)據(jù)簡(jiǎn)單的視為失效數(shù)據(jù)，不采用基于截尾數(shù)據(jù)的威布爾參數(shù)極大似然法，得到參數(shù)為m=3.09，η=14 739，是否采用基于截尾的威布爾極大似然法的參數(shù)結(jié)果比較如圖2所示。

圖2　是否采用基于截尾的威布爾極大似然法的結(jié)果比較

根據(jù)ARP 5150附錄D1.6中給的實(shí)例，假定燃料管泄漏的相關(guān)故障點(diǎn)火源的失效率λ=0.000 05，采用式(2)，可得燃料管泄漏在相同檢查間隔τ下的不同使用時(shí)刻T的風(fēng)險(xiǎn)的變化，如圖3所示。以及相同使用時(shí)刻T下的不同檢查間隔τ的風(fēng)險(xiǎn)變化，如圖4所示。

圖3　相同檢查間隔τ下不同使用時(shí)刻T的風(fēng)險(xiǎn)

圖4　相同使用時(shí)刻T下不同檢查間隔τ的風(fēng)險(xiǎn)

圖3中共有5條曲線，其中曲線“——”代表在不同使用時(shí)間，采用500h的檢查間隔的風(fēng)險(xiǎn)。除5條曲線外，圖3中給出了一條豎直線，該豎直線與5條曲線分別相交，從下而上的5個(gè)交點(diǎn)，分別表示在T=7 000h時(shí)，如檢查間隔分別采用100h、200h、300h、400h、500h的對(duì)應(yīng)風(fēng)險(xiǎn)，風(fēng)險(xiǎn)依次為：7.91×10-8、6.36×10-7、2.16×10-6、5.14×10-6、1.01×10-5。圖4中也有5條曲線，其中曲線“——”代表當(dāng)使用時(shí)間為9 000h時(shí)，采用不同檢查間隔的風(fēng)險(xiǎn)。圖4中的豎直線與5條曲線的交點(diǎn)，從下而上分別表示，當(dāng)使用時(shí)間為5 000h、6 000h、7 000h、8 000h、9 000h，如采用400h的檢查間隔所對(duì)應(yīng)風(fēng)險(xiǎn)，風(fēng)險(xiǎn)依次為：3.22×10-6、4.18×10-6、5.14×10-6、6.05×10-6、6.87×10-6。根據(jù)圖3、圖4分析可以看出，使用時(shí)間越長(zhǎng)風(fēng)險(xiǎn)越高，檢查間隔越長(zhǎng)風(fēng)險(xiǎn)也越高，這一結(jié)果符合故障發(fā)生的普遍規(guī)律。如果按照ARP5150 附錄D1.6中的方法進(jìn)行計(jì)算的話，那么得到的安全風(fēng)險(xiǎn)是一個(gè)定值，不僅不符合燃料管隨著使用時(shí)間推移，故障發(fā)生概率不斷增加的情況(如圖2所示，不可靠度隨著時(shí)間的推移在減低，即故障發(fā)生概率在增加)；而且沒(méi)有考慮檢查間隔的影響，與實(shí)際脫鉤。同時(shí)，本文提供的方法可以方便評(píng)估在任何時(shí)刻采取某種檢查間隔所帶來(lái)的定量風(fēng)險(xiǎn)，為快速維修決策，以及確保持續(xù)適航提供了有力的手段。

4　結(jié)語(yǔ)

民用飛機(jī)的隱蔽性故障對(duì)飛機(jī)構(gòu)成了重大潛在風(fēng)險(xiǎn)，對(duì)隱蔽性故障開(kāi)展定量風(fēng)險(xiǎn)評(píng)估是保持飛機(jī)持續(xù)適航的重要手段。本文參考APR5150理念，結(jié)合實(shí)際需求，建立了一種更實(shí)用的隱蔽故障風(fēng)險(xiǎn)的定量評(píng)估方法，對(duì)APR5150中的方法進(jìn)行了擴(kuò)充與完善，具有重要的理論意義和實(shí)際應(yīng)用價(jià)值。

參考文獻(xiàn)：

[1]蘇義.機(jī)載設(shè)備隱蔽故障及維修策略[J].航空制造技術(shù)，2010(12):86-88.

[2]賈寶惠，劉濤，楊杭，等.民機(jī)隱蔽故障維修間隔優(yōu)化方法研究[J].航空制造技術(shù)，2015(S1):20-23.

[3]郭中偉.淺談民用適航管理與航空安全[J].江蘇航空，2014(2):21-22.

[4]LIX,CAIJ,ZUOH.ResearchonReliabilitySequentialComplianceMethodforRollingBearings[J].JournalofInformationandComputationalScience，2015,12(14):5309-5318.

[5]WARDM,MCDONALDN,MORRISONR,etal.APerformanceImprovementCaseStudyinAircraftMaintenanceandItsImplicationsforHazardIdentification[J].Ergonomics，2010,53(2):247-267.

[6]王少萍.大型飛機(jī)機(jī)載系統(tǒng)預(yù)測(cè)與健康管理關(guān)鍵技術(shù)[J].航空學(xué)報(bào)，2014(6):1459-1472.

[7]于宏軍，韓建軍，張華，等.航空發(fā)動(dòng)機(jī)健康管理系統(tǒng)標(biāo)準(zhǔn)探討[J].航空標(biāo)準(zhǔn)化與質(zhì)量，2012(4):9-11.

[8]楊洲，景博，張劼，等.機(jī)載系統(tǒng)故障預(yù)測(cè)與健康管理驗(yàn)證與評(píng)估方法[J].測(cè)控技術(shù)，2012(3):101-104.

[9]BAUSSARONJ,MIHAELAB,LéOG,etal.ReliabilityAssessmentBasedonDegradationMeasurements:HowtoCompareSomeModels?[J].ReliabilityEngineering&SystemSafety，2014,131:236-241.

[10]RYANPC,STEWARTMG,SpencerN,etal.ReliabilityAssessmentofPowerPoleInfrastructureIncorporatingDeteriorationandNetworkMaintenance[J].ReliabilityEngineering&SystemSafety，2014,132:261-273.

[11]蔣仁言，胡浩然.對(duì)由預(yù)防維修活動(dòng)所產(chǎn)生的截尾數(shù)據(jù)的初步分析[C].中國(guó)運(yùn)籌學(xué)會(huì)可靠性分會(huì)第九屆可靠性學(xué)術(shù)會(huì)議論文集.長(zhǎng)沙:2013136-141.

[12]鄢偉安，宋保維，段桂林，等.威布爾部件的經(jīng)驗(yàn)貝葉斯評(píng)估[J].系統(tǒng)工程理論與實(shí)踐，2013(11):2980-2985.

[13]王繼霞，申培萍.定時(shí)截尾下Weibull分布參數(shù)估計(jì)的EM算法[J].河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2009(2):9-11.

(責(zé)任編輯楊繼森)

doi:10.11809/scbgxb2016.06.038

收稿日期：2015-12-17；修回日期：2016-01-22

基金項(xiàng)目:民航局重點(diǎn)科技項(xiàng)目(MHDR20150103)

作者簡(jiǎn)介：戴順安(1962—)，男，高級(jí)工程師，主要從事民用飛機(jī)適航審定研究。

中圖分類號(hào)：V328

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：2096-2304(2016)06-0162-04

ResearchonQuantitativeRiskAssessmentMethodforHiddenFailureofCivilAircraft

DAIShun-an1,WANGYe1,CAIJing2

(1.ShanghaiAircraftCertificationCenterofChinaCivilAviation,Shanghai200335,China; 2.NanjingUniversityofAeronauticsandAstronautics,Nanjing211106,China)

Abstract:Risk assessment of hidden failure should be carried out to maintain the continued airworthiness. With the concept of ARP 5150, a quantitative risk assessment method was researched based on the characteristics of hidden failure. A maximum likelihood estimation method with censored data for Weibull distribution was adopted, and the characteristic of hidden failure was determined by maintenance data; The factors affecting the risk of hidden failure were analyzed and quantitative risk assessment method was established. Finally, the method was applied by an example, and the results show the effectiveness and applicability of the proposed method.

Key words:civil aircraft; risk; hidden failure; quantitative assessment