王志平,馬旻昱,張春曉,蘇景新

(中國民航大學(xué) 天津市民用航空器適航與維修重點實驗室，天津 300300)

基于馬爾可夫鏈的飛機結(jié)構(gòu)腐蝕狀態(tài)預(yù)測

王志平,馬旻昱,張春曉,蘇景新

(中國民航大學(xué) 天津市民用航空器適航與維修重點實驗室，天津 300300)

腐蝕損傷是最常見的飛機結(jié)構(gòu)損傷形式之一，會降低結(jié)構(gòu)的剩余強度和耐久性。針對在役民用典型機型飛機易腐蝕結(jié)構(gòu)，首先將腐蝕狀態(tài)劃分為5個等級；其次擬合遭受腐蝕的飛機結(jié)構(gòu)剩余厚度遞減函數(shù)，給出腐蝕狀態(tài)轉(zhuǎn)移概率矩陣，建立馬爾可夫狀態(tài)預(yù)測模型。最后收集某型飛機后貨艙隔框腐蝕數(shù)據(jù)作為算例，預(yù)測了隔框結(jié)構(gòu)在潛伏期之后的腐蝕狀態(tài)。結(jié)果發(fā)現(xiàn)使用26個月之后，沿海地區(qū)飛機隔框結(jié)構(gòu)腐蝕已達(dá)到三級狀態(tài)。預(yù)測結(jié)果與工程實踐相符，說明馬爾可夫模型適用于飛機結(jié)構(gòu)腐蝕狀態(tài)預(yù)測。

飛機結(jié)構(gòu)；腐蝕；馬爾可夫模型；狀態(tài)預(yù)測

腐蝕損傷是最常見的飛機結(jié)構(gòu)損傷形式之一。嚴(yán)重的結(jié)構(gòu)腐蝕不但直接影響飛機的飛行安全，降低飛機的服役期限[1],還會加重機務(wù)維修工作負(fù)擔(dān)，帶來高額的維修費用。因此，針對在役民用典型機型飛機易腐蝕結(jié)構(gòu)，探索腐蝕檢修數(shù)據(jù)的規(guī)律，進行飛機結(jié)構(gòu)腐蝕退化性能預(yù)測，對減少腐蝕帶來的經(jīng)濟損失、確保飛機結(jié)構(gòu)的安全運行有著重要的工程實踐意義。

目前，國內(nèi)外關(guān)于飛機結(jié)構(gòu)腐蝕的研究主要集中在飛機主體材料的腐蝕損傷原理與檢測、腐蝕防護與修復(fù)方面[2-7]，而基于民航在役飛機的腐蝕數(shù)據(jù)，對飛機結(jié)構(gòu)性能腐蝕退化預(yù)測的研究相對較少。如張騰等人[8]提出了通過預(yù)腐蝕疲勞試驗確定飛機結(jié)構(gòu)壽命包線的腐蝕影響系數(shù)法，建立了典型服役環(huán)境中飛機金屬結(jié)構(gòu)腐蝕疲勞關(guān)鍵件的剩余壽命預(yù)測方法。

從作者掌握的文獻來看，現(xiàn)有理論模型未充分考慮腐蝕維修工程實踐的時變性和樣本采集的局限性，使得模型參數(shù)帶有很大的不確定性，難以滿足結(jié)構(gòu)腐蝕預(yù)測和可靠性動態(tài)管理的要求。馬爾可夫理論能夠很好地處理變量的隨機變化過程并進行動態(tài)預(yù)測，在道路、橋梁、管道劣化模型方面都有成功應(yīng)用的案例[9-11]。Provan和Rodriguez[12]最先使用非均勻馬爾可夫過程模擬腐蝕坑深度的發(fā)展，Hong[13]等人利用馬爾可夫過程模擬了點蝕的增長階段。因此，本工作采用馬爾可夫理論進行飛機結(jié)構(gòu)腐蝕狀態(tài)的預(yù)測研究。

1　飛機結(jié)構(gòu)腐蝕狀態(tài)

馬爾可夫過程[14-16]是具有無后效性的隨機過程。后效性是指系統(tǒng)未來所處的狀態(tài)只與現(xiàn)在系統(tǒng)狀態(tài)有關(guān),而與系統(tǒng)過去的狀態(tài)無關(guān)。通過對某航空公司飛機結(jié)構(gòu)腐蝕檢測數(shù)據(jù)分析可知：飛機結(jié)構(gòu)腐蝕是動態(tài)系統(tǒng)的狀態(tài)轉(zhuǎn)移過程，該過程具有隨機性；飛機結(jié)構(gòu)未來的腐蝕狀態(tài)只與系統(tǒng)所處的腐蝕狀態(tài)有關(guān)，而與過去的腐蝕狀態(tài)無關(guān)。可見飛機結(jié)構(gòu)腐蝕具有馬爾科夫性。并且飛機結(jié)構(gòu)腐蝕狀態(tài)間的轉(zhuǎn)移概率僅與腐蝕狀態(tài)有關(guān)，而與時刻n無關(guān)，因此將飛機結(jié)構(gòu)腐蝕看作齊次馬爾可夫鏈。

結(jié)構(gòu)腐蝕狀態(tài)可以用多種特征量來表示，例如腐蝕深度、面積、體積等。由于本工作研究的飛機結(jié)構(gòu)腐蝕形態(tài)以點蝕為主，所以選取腐蝕深度作為度量腐蝕損傷的代表值[17-18]。航空器結(jié)構(gòu)持續(xù)完整性大綱[19-20]通過不影響結(jié)構(gòu)件極限設(shè)計強度可以除去的腐蝕最大量(通常指材料厚度)將飛機結(jié)構(gòu)腐蝕劃分為三個等級，即1級腐蝕、2級腐蝕、3級腐蝕。根據(jù)民航飛機結(jié)構(gòu)維修的施工慣例，對不同程度的腐蝕損傷采取的維修措施又有區(qū)別，因此可進一步將三個等級的結(jié)構(gòu)腐蝕劃分為5個狀態(tài)，即腐蝕狀態(tài)空間S={1，2，3，4，5}，如表1所示。

表1　腐蝕深度的狀態(tài)劃分Tab. 1　Corrosion states divided by maintenance measurements

2　馬爾可夫預(yù)測模型

飛機不同結(jié)構(gòu)的極限強度不同，而且腐蝕環(huán)境也存在差異，因此需針對不同結(jié)構(gòu)，不同站位的結(jié)構(gòu)腐蝕數(shù)據(jù)分別進行腐蝕深度預(yù)測。本工作在對飛機結(jié)構(gòu)腐蝕狀況進行預(yù)測時，主要是基于收集到的剩余厚度數(shù)據(jù)(剩余厚度等于原厚度與腐蝕深度的差值)建立飛機結(jié)構(gòu)腐蝕狀況預(yù)測模型。

2.1飛機結(jié)構(gòu)剩余厚度遞減模型

依據(jù)飛機維修時得到的腐蝕數(shù)據(jù)確定腐蝕發(fā)展趨勢及狀態(tài)分布建立馬爾可夫預(yù)測模型，但是在實際腐蝕檢測過程中，缺乏充分的檢測數(shù)據(jù)，需要通過檢測得到的數(shù)據(jù)近似擬合出飛機結(jié)構(gòu)剩余厚度的降低趨勢線。飛機結(jié)構(gòu)腐蝕屬于電化學(xué)腐蝕[19]，初始階段腐蝕速率較快，腐蝕速率隨著時間逐漸減慢，因而可假設(shè)剩余厚度按以下指數(shù)形式發(fā)展[20]：

(1)

式中:G,h均為系數(shù);t為前一次維修到本次維修的間隔時間(month);δ為剩余厚度(mm)。根據(jù)飛機維修得到的腐蝕數(shù)據(jù)，利用最小二乘法，擬合成式(1)指數(shù)函數(shù)形式，分析剩余厚度的發(fā)展趨勢，確定結(jié)構(gòu)腐蝕狀態(tài)。

2.2狀態(tài)轉(zhuǎn)移矩陣

這里將飛機結(jié)構(gòu)腐蝕狀態(tài)劃分為5個狀態(tài)，則轉(zhuǎn)移概率矩陣P的一般形式為：

(2)

式中，轉(zhuǎn)移概率值Pmn反映腐蝕程度從狀態(tài)1一直到狀態(tài)5的變化過程。假設(shè)飛機結(jié)構(gòu)腐蝕不采取防腐措施，隨著服役時間的增長，最大腐蝕深度將越來越大，其腐蝕狀態(tài)等級保持不變或變差。當(dāng)腐蝕狀態(tài)達(dá)到5時，便不能轉(zhuǎn)移到其他狀態(tài)，停留在此狀態(tài)。此時轉(zhuǎn)移概率矩陣P的形式為：

(3)

通過公式(1)，求解各時間間隔t(1,2,…,ti…,tn)下的剩余厚度，并按照表1確定腐蝕狀態(tài)。當(dāng)時間間隔為tn時，腐蝕狀態(tài)從4變?yōu)?，此時腐蝕狀態(tài)將一直停留在狀態(tài)5，表明飛機結(jié)構(gòu)腐蝕狀態(tài)已達(dá)到最嚴(yán)重狀態(tài)，終止時間t的取值。從而得到一系列隨時間變化的腐蝕狀態(tài)數(shù)據(jù)。利用統(tǒng)計分析法計算飛機結(jié)構(gòu)腐蝕系統(tǒng)在ti時由狀態(tài)i到狀態(tài)j的轉(zhuǎn)移概率值Pij，計算過程如下：①統(tǒng)計飛機結(jié)構(gòu)腐蝕狀態(tài)為i的次數(shù)gi；②計算出處經(jīng)過一步轉(zhuǎn)移從狀態(tài)i到狀態(tài)j的次數(shù)gij，則Pij的計算公式如下：

(3)

進而求出轉(zhuǎn)移概率矩陣P，反映出從開始腐蝕到五級腐蝕的變化規(guī)律。

2.3腐蝕狀態(tài)預(yù)測模型

應(yīng)用馬爾可夫鏈對飛機某結(jié)構(gòu)腐蝕剩余厚度進行預(yù)測時，需要確定剩余厚度的初始狀態(tài)和轉(zhuǎn)移概率矩陣。將最后一次檢測得到的剩余厚度作為初始狀態(tài)，用向量K=[K1,K2,K3,K4,K5]表示，其中K1,…,K5為結(jié)構(gòu)腐蝕狀況在各個狀態(tài)的概率值，范圍為[0,1]。當(dāng)初始狀態(tài)K和轉(zhuǎn)移概率矩陣P確定后，根據(jù)齊次馬爾可夫鏈性質(zhì)，在最后一次檢測之后第m次(m≥1)檢測時的腐蝕狀態(tài)分布預(yù)測值為：

(5)

3　計算實例

3.1確定飛機結(jié)構(gòu)腐蝕狀態(tài)

利用以上模型對沿海地區(qū)某型飛機后貨艙767站位隔框腐蝕深度進行預(yù)測。一般需要使用同一架飛機結(jié)構(gòu)件在不同時間測得的結(jié)構(gòu)腐蝕數(shù)據(jù)。但是在對飛機結(jié)構(gòu)的例行維護和檢查中，一旦發(fā)現(xiàn)腐蝕損傷，會立即進行修理(即去除腐蝕層，并恢復(fù)結(jié)構(gòu)強度和防護層)，因此很難在同架飛機上得到同一位置的腐蝕深度的連續(xù)發(fā)展數(shù)據(jù)。本工作使用的飛機結(jié)構(gòu)件腐蝕數(shù)據(jù)取自同一個服役機場的同一機型飛機，且腐蝕構(gòu)件的材料、加工工藝、防護體系相同，將這些測量到的腐蝕損傷數(shù)據(jù)作為一個母體來考慮，亦可以得到飛機在不同服役時間狀態(tài)下的腐蝕深度的連續(xù)發(fā)展數(shù)據(jù)。算例中采用的數(shù)據(jù)如表2所示，數(shù)據(jù)來自某機場同一型號飛機后貨艙區(qū)域的767站位隔框，腐蝕發(fā)展時間為距離上次維護檢查的時間(即維護間隔)。

表2　某型飛機767站位隔框在不同腐蝕發(fā)展時間 的腐蝕數(shù)據(jù)Tab. 2　Corrosion data of 767 aft cargo frame at different corrosion time

根據(jù)表2中時間和剩余厚度數(shù)據(jù)，通過最小二乘法，運用SPSS進行曲線擬合，擬合成公式(1)的指數(shù)函數(shù)形式。擬合結(jié)果如圖1所示。

圖1　指數(shù)函數(shù)擬合效果圖Fig. 1　Fitting effect chart of exponential function

由表3可知擬合函數(shù)的擬合優(yōu)度R2接近于1，說明擬合效果很好，而且該遞減模型的顯著性檢驗(F檢驗)P值小于0.05，進一步說明擬合的模型顯著，該模型可用作飛機結(jié)構(gòu)剩余厚度遞減模型，因此剩余厚度隨時間發(fā)展的函數(shù)關(guān)系式為：

(6)

表3　擬合函數(shù)的檢驗結(jié)果及參數(shù)Tab. 3　Test results and parameters of fitting function

已知隔框結(jié)構(gòu)的原始厚度為1.6 mm，根據(jù)飛機結(jié)構(gòu)中隔框的極限強度[21]進行腐蝕狀態(tài)劃分，如表4所示。

表4　隔框腐蝕狀態(tài)劃分Tab. 4　Corrosion states of frame divided into 5 levels by corrosion depth

表4中劃分了隔框腐蝕狀態(tài)，由表2可知最后一次維修的間隔時間為19個月，最大腐蝕深度為0.183 mm，腐蝕狀態(tài)為3，則初始狀態(tài)為K=[0,0,1,0,0]。

從檢測到腐蝕深度開始對時間進行取值，帶入公式(6)，計算剩余厚度，得到腐蝕狀態(tài)變化情況，如表5所示。

表5　不同腐蝕發(fā)展時間的腐蝕數(shù)據(jù)及對應(yīng)狀態(tài)Tab. 5　Corrosion data and corresponding states at different corrosion times

3.2確定腐蝕狀態(tài)轉(zhuǎn)移矩陣

依據(jù)表5中的飛機結(jié)構(gòu)腐蝕數(shù)據(jù)及其對應(yīng)狀態(tài)，統(tǒng)計處于各腐蝕狀態(tài)的次數(shù)gi和腐蝕狀態(tài)之間的轉(zhuǎn)移次數(shù)gij，根據(jù)公式(4)求解轉(zhuǎn)移概率矩陣中各狀態(tài)轉(zhuǎn)移概率的值Pij，得到轉(zhuǎn)移概率矩陣P：

(7)

3.3腐蝕狀態(tài)預(yù)測及驗證

根據(jù)公式(5)進行馬爾可夫鏈進行預(yù)測，由第19個月的腐蝕狀態(tài)開始，運用MATLAB計算得到腐蝕狀態(tài)分布的預(yù)測值，如表6所示。

通過對表6中數(shù)據(jù)進行分析可知，飛機后貨艙隔框腐蝕狀態(tài)在26個月的狀態(tài)向量K26=[0，0，0.133 5，0.381 2，0.485 4]，其中狀態(tài)5的概率值最大，為0.485 4。可見飛機后貨艙隔框腐蝕狀態(tài)最可能進入狀態(tài)5，達(dá)到三級腐蝕，因此需要報告并進行隔框更換。

表6　腐蝕狀態(tài)的概率值在不同腐蝕發(fā)展時間的 預(yù)測結(jié)果Tab. 6　Predictive results of probability values of corrosion states at different corrosion times

以第12個月的腐蝕狀態(tài)為初始狀態(tài)，運用MATLAB計算得到13個月到19個月的飛機后貨艙隔框腐蝕狀態(tài)馬爾可夫預(yù)測值，與13個月到19個月的實際數(shù)據(jù)對比以驗證馬爾可夫模型準(zhǔn)確性。如表7所示。

表7　馬爾可夫模型的驗證結(jié)果Tab. 7　Verification results of Markov modeling

通過與實際腐蝕深度狀態(tài)變化進行對比，發(fā)現(xiàn)利用馬爾可夫鏈預(yù)測的結(jié)果與實際狀態(tài)變化基本一致，說明馬爾可夫鏈適用于腐蝕條件下飛機結(jié)構(gòu)性能退化預(yù)測。

4　結(jié)論

基于飛機結(jié)構(gòu)腐蝕的基本原理，通過恰當(dāng)方法選擇腐蝕數(shù)據(jù)，應(yīng)用馬爾可夫鏈對飛機結(jié)構(gòu)腐蝕狀態(tài)進行預(yù)測。通過計算實例，發(fā)現(xiàn)運用馬爾可夫鏈預(yù)測與實際情況基本符合，表明馬爾可夫鏈可應(yīng)用于飛機結(jié)構(gòu)腐蝕狀態(tài)預(yù)測。主要結(jié)論如下：

(1) 某沿海地區(qū)典型機型后貨艙隔框在腐蝕環(huán)境下的剩余厚度遞減函數(shù)關(guān)系式為δ=1.818 exp(-0.013t)。

(2) 該機型飛機后貨艙隔框在沿海地區(qū)在26個月進入腐蝕狀態(tài)5，達(dá)到三級腐蝕。

(3) 以第11個月的飛機結(jié)構(gòu)腐蝕為初始狀態(tài)，對12到19個月的結(jié)構(gòu)腐蝕狀態(tài)進行預(yù)測，發(fā)現(xiàn)與實際情況基本一致。

因此,應(yīng)用馬爾可夫鏈進行飛機結(jié)構(gòu)腐蝕狀態(tài)預(yù)測是合理可行的。

[1]張有宏,呂國志,李仲,等. 鋁合金結(jié)構(gòu)腐蝕疲勞裂紋擴展與剩余強度研究[J]. 航空學(xué)報,2007,2(28):332-335.

[2]JONES K,DAVID W H. Prior corrosion and fatigue of 2024-T3 aluminum alloy[J]. Corrosion science,2006,48(10):3109-3122.

[3]RUSSO S,SHARP P K,MILLS T B,et al. The influence of the environment and corrosion on the structural integrity of aircraft materials[J]. Fatigue and Fracture of Engineering Materials and Structures,2009,32(6):464-472.

[4]賀崇武,蔡新鎖,李素強. 飛機典型連接件腐蝕及腐蝕疲勞試驗研究[J]. 腐蝕與防護,2006,27(3):118-121.

[5]陳群志,孫詐東,陸維忠,等. SEBF/SLF重腐蝕防護涂層應(yīng)用于典型飛機結(jié)構(gòu)中防腐性能綜合評定[J]. 中國腐蝕與防護學(xué)報,2005,25(6):365-368.

[6]張蕾,劉慧叢,朱立群,等. 局部修復(fù)對典型環(huán)境中飛機外部結(jié)構(gòu)件腐蝕的影響及腐蝕當(dāng)量關(guān)系修正[J]. 腐蝕與防護,2014,35(8):833-837.

[7]耿德平,宋慶功. 航空材料腐蝕疲勞研究進展[J]. 腐蝕與防護,2011,32(3):206-209.

[8]張騰,何宇廷,張海威,等. 基于壽命包線的飛機金屬結(jié)構(gòu)壽命預(yù)測方法[J]. 南京航空航天大學(xué)學(xué)報,2014,46(3):413-418.

[9]袁賡,王樹剛,黃一. 基于馬爾可夫鏈的埋地燃?xì)怃摴芄鼙诟g預(yù)測[J]. 哈爾濱工業(yè)大學(xué)學(xué)報,2010,42(8):1328-1331.

[10]吳明亮,郝點,劉錦昆. 基于灰色馬爾可夫理論的油氣管道腐蝕剩余壽命預(yù)測[J]. 管道技術(shù)與設(shè)備,2008(5):43-45.

[11]CALEYO F,VELZQUEZ J C,VALOR A,et al. Markov chain modelling of pitting corrosion in underground pipelines[J]. Corrosion Science,2009,51(9):2197-2207.

[12]PROVAN J W,RODRIGUEZ E S. Development of a Markov description of pitting corrosion[J]. Corrosion,1989,45(3):178-192.

[13]HONG H P. Inspection and maintenance planning of pipeline under external corrosion considering generation of new defects[J]. Structural Safety,1999,21(3):203-222.

[14]胡迪鶴. 隨機環(huán)境中的馬爾可夫過程[M]. 北京:高等教育出版社,2011.

[15]錢衛(wèi)東,劉志強. 基于灰色馬爾可夫的道路交通事故預(yù)測[J]. 中國安全科學(xué)學(xué)報,2008,18(3):33-36.

[16]王永剛,呂學(xué)梅. 民航事故征候的灰色馬爾可夫預(yù)測[J]. 安全與環(huán)境學(xué)報,2008,8(1):163-165.

[17]AC-121-65航空器結(jié)構(gòu)持續(xù)完整性大綱[S].

[18]陳躍良,楊曉華,秦海勤. 飛機結(jié)構(gòu)腐蝕損傷分布規(guī)律研究[J]. 材料科學(xué)與工程,2002,20(3):378-380.

[19]穆志韜,熊玉平. 高強度鋁合金的腐蝕損傷分布規(guī)律研究[J]. 機械工程材料,2002,26(4):14-16.

[20]蘇維國,穆志韜. 基于損傷檢測的腐蝕疲勞壽命預(yù)測概率模型[J]. 裝備環(huán)境工程,2009,6(5):33-38.

[21]張行,田秀云,翟月英,等. 飛機結(jié)構(gòu)修理中的極限強度分析方法[J]. 中國民航學(xué)院學(xué)報,1986(1):1-8.

Application of Markov Chain to Prediction of Corrosion Conditions for Aircraft Structures

WANG Zhi-ping, MA Min-yu, ZHANG Chun-xiao, SU Jing-xin

(Tianjin Key Laboratory for Civil Aircraft Airworthiness and Maintenance, Civil Avation Unversity of China,Tianjin 300300, China)

The corrosion damage is one of the most common injuries in aircraft structures, which reduces the residual strength and durability of aircraft structures. Therefore, this study focused on the prediction of corrosion states of a typical civil aircraft structure based on Markov theory. Firstly, the corrosion states were divided into five levels. Secondly, the decreasing function of remaining thickness of corroded structures were fitted. The corrosion state transition probability matrix was calculated, and the corrosion depth of Markov predictive models was given. The corrosion states of aft cargo frame in the typical civil aircraft were adopted in the example. The corrosion state of the frame at the time after incubation period was predicted, and the result indicated that the frame corroded into the most serious state at coastal site after 26 months. The prediction of Markov model is conformed to the engineering practice, indicating that Markov model is applicable to the prediction of aircraft structural corrosion state.

aircraft structure; corrosion; Markov modeling; state prediction

10.11973/fsyfh-201510009

2014-10-15

國家自然科學(xué)基金(21303261); 中央高校基本科研業(yè)務(wù)費中國民航大學(xué)專項(3122013D004); 天津市民用航空器適航與維修重點實驗室開放課題

王志平(1963-)，教授，博士，從事材料表面改性研究，022-24092074，zpwang@cauc.edu.cn

TG174; V215.5

A

1005-748X(2015)10-0942-05