航空發(fā)動(dòng)機(jī)限壽件疲勞可靠度計(jì)算新方法

游令非，張建國(guó)，周霜，杜小松

1. 北京航空航天大學(xué) 可靠性與系統(tǒng)工程學(xué)院, 北京 100083 2. 北京航空航天大學(xué) 可靠性與環(huán)境工程技術(shù)國(guó)防重點(diǎn)實(shí)驗(yàn)室, 北京 100083

在航空發(fā)動(dòng)機(jī)適航要求中，將原發(fā)失效能夠引起發(fā)動(dòng)機(jī)危害性影響的部件定義為發(fā)動(dòng)機(jī)限壽件(Engine Life Limited Part,ELLP)，如旋轉(zhuǎn)輪盤和大型旋轉(zhuǎn)封嚴(yán)裝置等，其可靠水平直接決定著發(fā)動(dòng)機(jī)整體可靠性和工作性能，設(shè)計(jì)中主要通過降低ELLP的失效概率來提高整機(jī)的可靠性和安全性[1-2]。而對(duì)航空發(fā)動(dòng)機(jī)造成危害性后果的事件有很大部分是由其機(jī)械結(jié)構(gòu)的疲勞所導(dǎo)致。對(duì)壓氣機(jī)盤/渦輪盤進(jìn)行基于疲勞可靠度計(jì)算的概率風(fēng)險(xiǎn)評(píng)估成為航空發(fā)動(dòng)機(jī)型號(hào)適航取證過程以及保證其結(jié)構(gòu)完整性的關(guān)鍵技術(shù)和重要實(shí)施步驟之一。

目前國(guó)內(nèi)外針對(duì)以輪盤為代表的ELLP疲勞可靠性做了大量研究，美國(guó)軍方[3]首先將概率統(tǒng)計(jì)理論引入到傳統(tǒng)輪盤疲勞壽命預(yù)測(cè)方法中，形成了輪盤壽命可靠性設(shè)計(jì)方法。Melis等[4]對(duì)航空燃?xì)廨啓C(jī)壓氣機(jī)盤的壽命進(jìn)行了預(yù)測(cè)和可靠性分析。Gao等[5]應(yīng)用分布式協(xié)同響應(yīng)面法對(duì)航空發(fā)動(dòng)機(jī)渦輪進(jìn)行了疲勞損傷的可靠性分析。Zhu 等[6-7]分別建立了渦輪葉片輪盤疲勞可靠性評(píng)估的計(jì)算試驗(yàn)框架和概率分析框架。將超速試驗(yàn)與隨機(jī)有限元分析相結(jié)合，量化了試驗(yàn)數(shù)據(jù)、材料性能和載荷的不確定性以及材料變化和荷載變化對(duì)可靠性結(jié)果的影響。Hu 等[8]針對(duì)合金渦輪盤的鎳基高溫合金，建立了基于外加機(jī)械功密度的概率模型，利用線性異方差函數(shù)對(duì)蠕變疲勞壽命中的非常數(shù)偏差進(jìn)行評(píng)價(jià)。Wang等[9]探索了渦輪盤低周疲勞壽命預(yù)測(cè)的區(qū)域可靠性方法以及裂紋的閉合行為。李巖等[10]提出了一種基于Kriging和蒙特卡羅半徑外自適應(yīng)重要抽樣(Monte Carlo Radius-Outside Adaptive Importance Sampling,MCROAIS)方法混合的結(jié)構(gòu)概率風(fēng)險(xiǎn)評(píng)估方法，構(gòu)建了壓氣機(jī)盤風(fēng)險(xiǎn)概率模型。陳志英等[11]結(jié)合應(yīng)力應(yīng)變場(chǎng)強(qiáng)法和響應(yīng)面法建立了輪盤疲勞可靠性模型并進(jìn)行了可靠度計(jì)算。除此以外，隨機(jī)過程、神經(jīng)網(wǎng)絡(luò)和支持向量機(jī)等方法也廣泛地應(yīng)用到相關(guān)問題中[12-14]。

目前開展ELLP結(jié)構(gòu)概率風(fēng)險(xiǎn)評(píng)估的難點(diǎn)在于非線性、小失效概率事件的定量要求驗(yàn)證問題。但是對(duì)于高度非線性極限狀態(tài)函數(shù)以及小失效概率事件，以上方法均有一定的局限:① 階矩法對(duì)于非線性較高的問題計(jì)算精度不高；② 響應(yīng)面法與數(shù)值模擬相結(jié)合使計(jì)算步驟增加，特別是針對(duì)小概率失效事件計(jì)算效率低下。因此，一些基于方差抽樣縮減技術(shù)的抽樣方法[15-17]因具有較高的精度和魯棒性成為解決該問題的選擇。

作為在可靠度計(jì)算中廣泛應(yīng)用的數(shù)值模擬方法，重要抽樣方法(Important Sampling, IS)[18]將抽樣密度函數(shù)的抽樣中心移到設(shè)計(jì)驗(yàn)算點(diǎn)(以下簡(jiǎn)稱驗(yàn)算點(diǎn))，可以使更多的樣本點(diǎn)集中在失效域，提高了抽樣的效率。但在未知驗(yàn)算點(diǎn)(或可靠度指標(biāo))的情況下，只能通過蒙特卡羅仿真(Monte Carlo Simulation,MCS)法或解析近似法求解。這極大地限制了它們的適用范圍。為了克服重要抽樣法對(duì)驗(yàn)算點(diǎn)和可靠度指標(biāo)已知的依賴性，產(chǎn)生了蒙特卡羅自適應(yīng)重要抽樣(Monte Carlo Adaptive Important Sampling,MCAIS)法的思想。通過自動(dòng)迭代不斷改進(jìn)MCS的重要抽樣密度函數(shù)，使得抽樣點(diǎn)越來越趨近于失效域，使抽樣中心逐漸趨近于驗(yàn)算點(diǎn)。Au和Beck[19]采用基于Metropolis算法的馬爾科夫鏈模擬方法生成失效域樣本，而后運(yùn)用自適應(yīng)密度估計(jì)方法構(gòu)造重要抽樣密度函數(shù)。Sankaran和Prakasb[20]將自適應(yīng)重要抽樣法應(yīng)用在大結(jié)構(gòu)系統(tǒng)的可靠性分析中。黃毅等[21]采用蒙特卡羅方法對(duì)中心柱設(shè)計(jì)過程中的主要影響因素——低周疲勞進(jìn)行可靠性評(píng)價(jià)。陳向前等[22]提出了一種基于樣本概率密度加權(quán)的采樣中心確定方法，該方法在自適應(yīng)重要抽樣的基礎(chǔ)上增加了有效抽樣中對(duì)失效概率貢獻(xiàn)大的樣本出現(xiàn)的概率。馬紀(jì)明等[23]提出一種改進(jìn)的自適應(yīng)重要抽樣方法，該方法通過失效樣本估計(jì)IS密度函數(shù)的初始參數(shù)，提高了算法的收斂時(shí)間。戴鴻哲等[24]提出了一種自適應(yīng)Metropolis算法和快速高斯變換技術(shù)的結(jié)構(gòu)可靠性分析高效自適應(yīng)重要抽樣方法。但在一些情況下，自適應(yīng)重要抽樣法并不能得到理想的驗(yàn)算點(diǎn)近似結(jié)果，比如變量的聯(lián)合概率密度函數(shù)的標(biāo)準(zhǔn)差比較大時(shí)，或者概率密度函數(shù)的梯度和極限狀態(tài)函數(shù)的梯度在某一區(qū)域大致相同時(shí)等。將會(huì)使近似驗(yàn)算點(diǎn)和真實(shí)驗(yàn)算點(diǎn)間具有較大的誤差。

為了克服以上問題，本文提出一種具有自更新機(jī)制的半徑外自適應(yīng)重要抽樣(Auto Updating Monte Carlo Radius-Outside Adaptive Importance Sampling,AUMCROAIS)疲勞可靠性分析方法，本方法可以兼顧搜索效率和計(jì)算精度，在高效地尋找近似驗(yàn)算點(diǎn)的同時(shí)，能較為精確地計(jì)算可靠度。本文提出的方法總體思路如下：首先利用MCAIS快速逼近驗(yàn)算點(diǎn)附近，在MCAIS算法停止后，在以此時(shí)的近似驗(yàn)算點(diǎn)為中心進(jìn)行極坐標(biāo)抽樣，并依次構(gòu)造主動(dòng)學(xué)習(xí)函數(shù)進(jìn)行評(píng)判，即對(duì)極坐標(biāo)抽樣點(diǎn)先進(jìn)行極限狀態(tài)函數(shù)的主動(dòng)學(xué)習(xí)函數(shù)的構(gòu)造，根據(jù)判據(jù)自動(dòng)選擇那些離極限狀態(tài)面近的抽樣點(diǎn)，而后針對(duì)這些點(diǎn)進(jìn)行最優(yōu)半徑的主動(dòng)學(xué)習(xí)函數(shù)的構(gòu)造并找到最小值點(diǎn)，此時(shí)，近似驗(yàn)算點(diǎn)得到了更新，而后以更新近似驗(yàn)算點(diǎn)為中心，利用半徑外重要抽樣計(jì)算可靠度，繼續(xù)適當(dāng)縮減方差并重復(fù)上述步驟，直到滿足收斂條件，得到最終的可靠度。本方法通過不斷的更新確定出最優(yōu)抽樣半徑，加速失效概率計(jì)算的收斂，不受限于自適應(yīng)抽樣的收斂條件，保證了效率同時(shí)提高了精度。AUMCROAIS疲勞可靠性分析方法解決了疲勞可靠性問題中小概率失效和極限狀態(tài)非線性較強(qiáng)的情況，支撐了適航規(guī)章要求(FAR33.75條款)[25]的實(shí)施，及對(duì)發(fā)動(dòng)機(jī)失效概率事件分析的要求。最后在工程案例上利用本方法分別和蒙特卡羅方法、自適應(yīng)抽樣法和一階可靠性方法進(jìn)行對(duì)比，說明了本方法的適用性、精確性和高效性。

1 MCROAIS方法原理

圖1 MCROAIS法流程圖
Fig.1 Flow chart of MCROAIS method

(1)

MCROAIS雖然相比MCS提高了抽樣效率，但是MCROAIS依然存在一些缺陷，表現(xiàn)在：① 小概率情形下且變量的標(biāo)準(zhǔn)差較大時(shí)，將會(huì)導(dǎo)致β在近驗(yàn)算點(diǎn)附近變化緩慢，將提前結(jié)束迭代過程；② 概率密度函數(shù)的梯度和極限狀態(tài)函數(shù)的梯度在某一區(qū)域大致相同，導(dǎo)致β變化緩慢，將提前結(jié)束迭代過程；③ 若非線性較高時(shí)同時(shí)收斂偏離驗(yàn)算點(diǎn)，每次迭代時(shí)標(biāo)準(zhǔn)差的減小將不會(huì)使近似驗(yàn)算點(diǎn)的收斂路線得到修正，反而將加速迭代結(jié)束，得到不精確的結(jié)果；④ 若收斂準(zhǔn)則嚴(yán)格(即提高精度要求)，則計(jì)算效率將大大增加。為此，在第2節(jié)中將討論如何基于MCROAIS法加以改進(jìn)并構(gòu)造新的抽樣方法的以同時(shí)兼顧抽樣效率和計(jì)算精度。

2 基于MCROAIS法的極坐標(biāo)抽樣和主動(dòng)學(xué)習(xí)函數(shù)

通過第1節(jié)分析可知，MCROAIS法并不能在任何極限狀態(tài)函數(shù)和變量分布，特別是小概率事件下得到理想的結(jié)果，本節(jié)將通過極坐標(biāo)抽樣和主動(dòng)學(xué)習(xí)函數(shù)的構(gòu)造，使其在滿足自適應(yīng)抽樣收斂條件下，繼續(xù)構(gòu)造和識(shí)別離驗(yàn)算點(diǎn)近的樣本點(diǎn)，從而達(dá)到自更新的目的。極坐標(biāo)抽樣即以每次自適應(yīng)迭代的結(jié)果得出的近似驗(yàn)算點(diǎn)為中心繼續(xù)構(gòu)造抽樣，進(jìn)而用主動(dòng)學(xué)習(xí)函數(shù)依次對(duì)貼近極限狀態(tài)和最優(yōu)半徑進(jìn)行篩選，更新出最優(yōu)點(diǎn)，用以構(gòu)造最優(yōu)半徑或者下一次迭代的起始點(diǎn)。極坐標(biāo)抽樣和主動(dòng)學(xué)習(xí)函數(shù)相結(jié)合的方法提高了抽樣效率，同時(shí)也考慮了計(jì)算精度。

2.1 極坐標(biāo)抽樣

以單次自適應(yīng)抽樣的近似驗(yàn)算點(diǎn)為中心進(jìn)行極坐標(biāo)抽樣，極坐標(biāo)抽樣基于布朗運(yùn)動(dòng)的軌跡概率分布，其特點(diǎn)是分別基于極角和極徑分布進(jìn)行抽樣，且樣本點(diǎn)數(shù)控制在很小的范圍內(nèi)，大大地提高了抽樣效率同時(shí)也使樣本點(diǎn)更多地接近近似驗(yàn)算點(diǎn)一側(cè)。下面以二維布朗運(yùn)動(dòng)加以簡(jiǎn)要說明，高維同理。設(shè)x方向和y方向的位移X和Y都是獨(dú)立的隨機(jī)變量，符合相同的正態(tài)分布N(0,σ2)。X和Y的聯(lián)合分布是fX和fY的乘積：

(2)

式中：σ2為方差。

位移R2=X2+Y2是隨機(jī)變量的函數(shù)，其分布函數(shù)為

(3)

式中：r為極徑；θ為極角；D是x2+y2≤r2的區(qū)域，通過變成極坐標(biāo)求得結(jié)果。對(duì)式(3)求導(dǎo)即得到概率密度分布:

(4)

(5)

(6)

則極徑r和極角θ的聯(lián)合分布為

pRΘ(r,θ)=

(7)

由于無規(guī)行走的極徑和極角是相互獨(dú)立的隨機(jī)量，故2個(gè)邊緣分布可以直接求積分，則在近似驗(yàn)算點(diǎn)x*對(duì)應(yīng)的近似最優(yōu)半徑β*下極徑和極角的分布為

(8)

(9)

圖2 極坐標(biāo)抽樣
Fig.2 Polar coordinate sampling

2.2 主動(dòng)學(xué)習(xí)函數(shù)的構(gòu)造

(10)

3 AUMCROAIS可靠度計(jì)算方法

AUMCROAIS法通過對(duì)最優(yōu)樣本的不斷更新，達(dá)到對(duì)小概率事件和強(qiáng)非線性極限狀態(tài)函數(shù)的真實(shí)驗(yàn)算點(diǎn)的近似，從而得到最優(yōu)半徑，繼而用MCROAIS法進(jìn)行可靠度求解。其包括內(nèi)循環(huán)和外循環(huán)2部分。

3.1 自更新機(jī)制

外循環(huán)包括以下幾個(gè)部分：以每次改進(jìn)的自適應(yīng)抽樣的近似驗(yàn)算點(diǎn)為中心進(jìn)行的極坐標(biāo)抽樣(pR(β*),pΘ(θ))，而后依次針對(duì)“近極限狀態(tài)面”和“最優(yōu)半徑”進(jìn)行主動(dòng)學(xué)習(xí)函數(shù)的構(gòu)造和判別，首先建立近極限狀態(tài)面主動(dòng)學(xué)習(xí)函數(shù)LG(xi)并判別找出那些近極限狀態(tài)面的點(diǎn)xi(i=1,2,…,N)，再建立最優(yōu)半徑主動(dòng)學(xué)習(xí)函數(shù)Lβ*(xi)并根據(jù)判別，找出近極限狀態(tài)面的點(diǎn)中半徑最小的點(diǎn)，即最優(yōu)驗(yàn)算點(diǎn)x*，并基于更新的驗(yàn)算點(diǎn)進(jìn)行半徑外重要抽樣可靠度計(jì)算R。最終進(jìn)行可靠度收斂判別，若不滿足要求，則進(jìn)行下一次外循環(huán)直至滿足外循環(huán)的可靠度收斂判別。基于以上步驟可同時(shí)兼顧效率和精度，同時(shí)構(gòu)建雙主動(dòng)學(xué)習(xí)函數(shù)，給極限狀態(tài)建立了統(tǒng)一的評(píng)判標(biāo)準(zhǔn)，有助于快速高效地識(shí)別出最優(yōu)驗(yàn)算點(diǎn)。

內(nèi)循環(huán)即根據(jù)本算法的特點(diǎn)和需要改進(jìn)的自適應(yīng)抽樣，旨在快速使近似驗(yàn)算點(diǎn)x*收斂于真實(shí)驗(yàn)算點(diǎn)x*附近。在迭代過程中，為使近似驗(yàn)算點(diǎn)快速收斂，內(nèi)循環(huán)中自適應(yīng)抽樣的收斂準(zhǔn)則可適當(dāng)放大，同時(shí)，由于第1次外循環(huán)已經(jīng)到達(dá)真實(shí)驗(yàn)算點(diǎn)附近，所以第2次以后的循環(huán)中(如果需要)，近似驗(yàn)算點(diǎn)的收斂步長(zhǎng)將適當(dāng)減小，則：在第1次循環(huán)中σhV(v)i=σhV(v)i-1×1，hV(v)為重要抽樣密度函數(shù)，第2次及以后的循環(huán)中σhV(v)i=σhV(v)i-1×0.9。同時(shí)，為加速內(nèi)循環(huán)的收斂，將γ設(shè)為1.02，以此來提高內(nèi)循環(huán)的計(jì)算精度。

3.2 算法流程

步驟4若沒有樣本點(diǎn)落在失效域內(nèi)，令βi=β**并把hV(v)移到Y(jié)**，否則比較β*>γβi-1，若不等式成立，令βi=β**，σhV(v)i=σhV(v)i-1×0.9(k≠1時(shí))。i=i+1，當(dāng)|βi-βi-1|<0.01時(shí)，轉(zhuǎn)到步驟2，直到滿足要求，記錄此時(shí)的Y=x*，σhV(v)k=σhV(v)i×0.9。

步驟5以x*為中心，進(jìn)行極坐標(biāo)抽樣(pR(r),pΘ(θ))，樣本點(diǎn)記為xi(i=1,2,…,N)。

步驟10根據(jù)下式確定第k次外循環(huán)下ELLP的失效概率:

4 應(yīng)用實(shí)例

將本文提出的方法應(yīng)用于某型發(fā)動(dòng)機(jī)低壓壓氣機(jī)輪盤，通過與MCS、MCROAIS法和FORM法進(jìn)行對(duì)比，驗(yàn)證本文方法的高效率、魯棒性和仿真精度。

針對(duì)影響發(fā)動(dòng)機(jī)安全性的關(guān)鍵件壓氣機(jī)輪盤的主要失效模式—低循環(huán)疲勞斷裂，本文選取某型發(fā)動(dòng)機(jī)低壓三級(jí)輪盤為典型關(guān)鍵件，開展發(fā)動(dòng)機(jī)輪盤疲勞可靠性計(jì)算。

4.1 壓氣機(jī)輪盤疲勞試驗(yàn)及仿真驗(yàn)證

采用立式旋轉(zhuǎn)試驗(yàn)系統(tǒng)開展疲勞壽命試驗(yàn)，試驗(yàn)在3 mmHg真空度下進(jìn)行[10]，如圖4所示。立式旋轉(zhuǎn)系統(tǒng)控制輪盤的轉(zhuǎn)速由低速500 r/min至高速9 550 r/min之間呈三角波式周期性變化，周期為30 s，如圖5所示。

通過該試驗(yàn)可確定其疲勞危險(xiǎn)位置，試驗(yàn)轉(zhuǎn)速控制在500-9 550-500 r/min，共完成55 300次試驗(yàn)器循環(huán)，約461 h，完成114 265次試驗(yàn)器循環(huán)后，對(duì)該盤進(jìn)行無損探傷檢查，在銷釘孔6點(diǎn)鐘方向處出現(xiàn)了穿透性裂紋并延伸至盤體，分析可知低壓三級(jí)輪盤銷釘孔裂紋屬于疲勞斷裂，起源于銷釘孔的內(nèi)端面。

同時(shí)，利用有限元方法對(duì)轉(zhuǎn)速為9 550 r/min的試驗(yàn)結(jié)果進(jìn)行仿真驗(yàn)證。渦輪盤的可靠性主要取決于結(jié)構(gòu)危險(xiǎn)點(diǎn)的應(yīng)力應(yīng)變分布。由于結(jié)構(gòu)上的形狀和負(fù)載完全對(duì)稱，因此本例中考慮輪盤的1/37，整體模型見圖6，1/37模型見圖7。為還原真實(shí)約束，仿真時(shí)添加空心銷軸，軸材料為3Cr13，輪盤材料為TC11。結(jié)合結(jié)構(gòu)循環(huán)對(duì)稱的幾何特性，對(duì)有限元模型施加以下位移約束：

1) 以發(fā)動(dòng)機(jī)輪盤的軸線為z軸建立圓柱坐標(biāo)系，在該圓柱坐標(biāo)系內(nèi)，如圖8所示，約束A、B兩面上全部節(jié)點(diǎn)沿旋轉(zhuǎn)方向θ的位移。

2) 在總體笛卡爾坐標(biāo)系下，約束圖8中銷釘兩端面——C、D面，以及輪盤E、F面上全部節(jié)點(diǎn)沿x方向——即沿銷釘軸向的位移。

葉片載荷施加在銷軸上，載荷合力為29 189.6 N(這一載荷為航空發(fā)動(dòng)機(jī)研究所相關(guān)部門給出的葉片產(chǎn)生的等效載荷)，方向?yàn)榈芽栕鴺?biāo)系y方向。參數(shù)詳情見表1。

圖3 AUMCROAIS方法流程
Fig.3 Flowchart of AUMCROAIS method

圖4 低壓壓氣機(jī)輪盤試驗(yàn)組裝
Fig.4 Test assembly of low pressure compressor disk

圖5 500～9 550 r/min轉(zhuǎn)速的變化示意圖
Fig.5 Schematic diagram of change of rotation speed between 500 and 9 550 r/min

圖6 低壓壓氣機(jī)輪盤模型
Fig.6 Low pressure compressor disk model

圖7 1/37輪盤和空心銷模型
Fig.7 Model of 1/37 disk and hollow pin

圖8 各面約束
Fig.8 All-sided constraints

選取表1中的數(shù)值作為仿真參數(shù)，利用ANSYS18.1進(jìn)行仿真。模擬結(jié)果如圖9所示。由結(jié)果可知應(yīng)力應(yīng)變水平最高處位于輪盤與銷軸連接處(位置在銷釘孔下部?jī)?nèi)端面，節(jié)點(diǎn)號(hào)1535)，此位置即結(jié)構(gòu)的破壞點(diǎn)，最大應(yīng)力為791.64 MPa。從而驗(yàn)證了試驗(yàn)的正確性和準(zhǔn)確性，下面將用試驗(yàn)結(jié)果進(jìn)行發(fā)動(dòng)機(jī)輪盤第一級(jí)疲勞應(yīng)力作用下(9 550 r/min)的疲勞可靠性計(jì)算。

表1 渦輪盤參數(shù)Table 1 Turbine disk parameters

圖9 9 550 r/min轉(zhuǎn)速下輪盤的應(yīng)力云圖
Fig.9 Stress cloud chart of disk at 9 550 r/min rotate speed

4.2 壓氣機(jī)輪盤疲勞可靠性分析

lgNlife=lgK-mlg(Smax-703.84)

(11)

式中：Nlife為疲勞載荷循環(huán)次數(shù)；σa為循環(huán)疲勞應(yīng)力幅值；Smax為應(yīng)力場(chǎng)強(qiáng)；m和K為材料參數(shù)；壓氣機(jī)輪盤材料為TC11，由《飛機(jī)結(jié)構(gòu)金屬材料力學(xué)性能手冊(cè)》[26]得到其材料參數(shù)K=6.278 4×1015，m=4.736，利用應(yīng)力-強(qiáng)度干涉理論構(gòu)建模型為

K(Smax-703.84)-m-nmax

(12)

K(Smaxi-703.84)-m-nmaxi

(13)

同時(shí)，由TC11進(jìn)行了材料級(jí)試驗(yàn)[8]，統(tǒng)計(jì)得到TC11材料參數(shù)m和K的均值、變異系數(shù)分別為μm=4.628、COVm=0.01，μK=6.135 26×1015，COVK=0.015，服從正態(tài)分布。

第1級(jí)應(yīng)力場(chǎng)強(qiáng)的均值和變異系數(shù)為μSmax1=791.64，COVSmax1=0.1，服從極值分布。利用本文中提出的AUMCROAIS疲勞可靠性分析方法，對(duì)第一級(jí)循環(huán)應(yīng)力作用下輪盤疲勞失效概率與應(yīng)力循環(huán)次數(shù)的關(guān)系進(jìn)行求解，同時(shí)用MCS, FORM法與MCROAIS法進(jìn)行比較。其中由于MCS抽樣規(guī)模較大，認(rèn)為其結(jié)果是最接近真實(shí)結(jié)果的并以MCS的結(jié)果作為精確度對(duì)比。

在發(fā)動(dòng)機(jī)輪盤第1級(jí)疲勞應(yīng)力作用下得到疲勞失效概率Pr隨載荷循環(huán)次數(shù)的關(guān)系如圖10所示，載荷循環(huán)從10 000～55 000，4種方法的可靠度指標(biāo)計(jì)算結(jié)果見表2，可以看出隨著疲勞載荷循環(huán)次數(shù)的增加，失效概率逐漸升高。

通過將本文方法的仿真效率、魯棒性以及仿真精度與MCS方法進(jìn)行比較，從而驗(yàn)證本文提出方法的適用性。

1) 仿真效率：通過MCROAIS收斂速度對(duì)比，可以驗(yàn)證本方法的仿真效率。如圖11所示，抽取12 000次、24 000次2種情況作為對(duì)比其可靠度指標(biāo)β的收斂過程，可以看出本方法在抽樣效率上的優(yōu)越性，由于結(jié)合了極坐標(biāo)抽樣和主動(dòng)學(xué)習(xí)函數(shù)，使近似驗(yàn)算點(diǎn)能迅速接近真實(shí)驗(yàn)算點(diǎn)x*，不斷更新最優(yōu)抽樣半徑；另外，在抽樣次數(shù)上與MCS和MCROAIS進(jìn)行對(duì)比，在每一級(jí)循環(huán)下，MCS需108次抽樣；MCROAIS分別需32次和25次迭代，每次1 000個(gè)樣本點(diǎn)，則MCROAIS共需32 000和25 000個(gè)樣本點(diǎn)；而本方法則分別需要2次外循環(huán)，11次內(nèi)循環(huán)； 2次外循環(huán)，8次內(nèi)循環(huán)；每次內(nèi)循環(huán)需抽樣1 000個(gè)樣本點(diǎn)，外循環(huán)36個(gè)樣本點(diǎn)，則本方法分別需要11 072和8 072個(gè)樣本點(diǎn)。同時(shí), 在與FORM法收斂速度對(duì)比上可以看出，AUMCROAIS的收斂速度與FORM基本相同，且在一定程度上要先于FORM到達(dá)真實(shí)MPP點(diǎn)附近。綜上所述，本方法在仿真效率上得到了很大的提升。

圖10 第一級(jí)疲勞應(yīng)力作用下輪盤失效概率與 疲勞應(yīng)力循環(huán)次數(shù)的關(guān)系
Fig.10 Relationship between failure probability of disk and number of fatigue stress cycles under the first-order fatigue stress

表2 可靠度指標(biāo)結(jié)果對(duì)比Table 2 Comparisons of reliability index results

圖11 可靠度指標(biāo)收斂速度對(duì)比
Fig.11 Comparison of convergence rate of reliability index

2) 魯棒性：針對(duì)高非線性的風(fēng)險(xiǎn)概率模型的失效概率仿真，驗(yàn)證本方法的魯棒性。如圖10所示，圍繞精確解(MCS結(jié)果)，AUMCROAIS和MCROAIS均有所波動(dòng)，但總體來說本文方法波動(dòng)較小，更具參考性。同時(shí)，F(xiàn)ORM法的結(jié)果隨著應(yīng)力循環(huán)次數(shù)的增加越來越偏離真實(shí)值，而本方法依然在小范圍內(nèi)波動(dòng)，說明本方法將在因高非線性而導(dǎo)致FORM偏離真實(shí)值或者不適用的情況下依然能給出近似解。

3) 仿真精度：通過仿真結(jié)果誤差分析，來驗(yàn)證本文方法的仿真精度。仿真精度見表2，可以看出可靠度指標(biāo)總體誤差不大于0.052 84，總體結(jié)果和FORM法基本相近，相比MCROAIS精度得到了很大提升。

從計(jì)算結(jié)果對(duì)比中可以看出，在不降低仿真精度和仿真效率得情況下，AUMCROAIS得到了較滿意的結(jié)果，解決了在小概率失效和非線性情況下MCROAIS并不能得到穩(wěn)定的精確解的問題，同時(shí)由于極坐標(biāo)抽樣和主動(dòng)學(xué)習(xí)函數(shù)的高效的抽樣和判別能力，使本方法的仿真效率大大提高，因此通過上述分析驗(yàn)證了本文所提方法的較高的計(jì)算效率、魯棒性和計(jì)算精度。

5 結(jié) 論

本文針對(duì)小概率失效事件和強(qiáng)非線性可靠性分析問題，綜合考慮了抽樣效率和抽樣精度，針對(duì)航空發(fā)動(dòng)機(jī)結(jié)構(gòu)疲勞可靠性分析問題，提出了AUMCROAIS方法，經(jīng)具體案例驗(yàn)證表明：

1) 本文提出的AUMCROAIS方法實(shí)施方便，計(jì)算誤差較小，對(duì)類似的小概率失效事件和強(qiáng)非線性可靠性分析問題具有一定的指導(dǎo)意義，例如本文中航空發(fā)動(dòng)機(jī)結(jié)構(gòu)疲勞可靠性分析案例中本方法的計(jì)算誤差最大不超過0.052 84，和MCS方法計(jì)算結(jié)果貼合度較高。

2) 本文提出的AUMCROAIS方法可以結(jié)合仿真，在完善隨機(jī)參數(shù)信息的情況下，可得到較精確的疲勞可靠度計(jì)算結(jié)果，同時(shí)計(jì)算效率大大提高，總算次數(shù)約為MCS方法的1/3左右。