黃曉宇, 王攀, 李海和, 張政

1.西北工業(yè)大學 力學與土木建筑學院, 陜西 西安710072;2.西北工業(yè)大學 飛行器可靠性研究所, 陜西 西安710072

渦輪盤是航空發(fā)動機的重要組成部分，其主要功能是連接渦輪葉片以及傳輸運送所需功率。航空發(fā)動機渦輪盤長期處于高溫、高壓等惡劣工作環(huán)境下，并且隨著飛行狀態(tài)的不同承受著交叉載荷的作用，其危險點處會發(fā)生局部塑性應變，導致其極易發(fā)生疲勞失效。因此，對渦輪盤進行疲勞可靠性分析是保障發(fā)動機安全的重要手段。傳統(tǒng)的可靠性分析模型認為只有達到累積損傷極限時渦輪盤才會發(fā)生失效，從安全狀態(tài)過度到失效狀態(tài)是一個突變的過程，但是由于多種不確定性因素的存在，渦輪盤在疲勞損傷累積過程中也可能失效，其失效邊界呈現(xiàn)出模糊性。

近些年來，已有許多學者針對模糊可靠性分析進行了大量的研究，Cai等[1]提出利用模糊狀態(tài)假設來代替?zhèn)鹘y(tǒng)的二元狀態(tài)假設，然后結(jié)合概率可靠性的分析方法，形成模糊可靠性理論。Mirakbari等[2]根據(jù)實際的數(shù)據(jù)，確定模糊隸屬度函數(shù)，提出一種基于馬爾可夫過程的方法來表示模糊失效概率。黃洪鐘[3]系統(tǒng)地分析影響結(jié)構可靠性的模糊性因素，提出結(jié)構模糊可靠性計算的統(tǒng)一模型，并定義幾種不同類型的隸屬函數(shù)。針對具有模糊失效狀態(tài)的可靠性分析算法方面，董玉革等[4]利用水平截集將模糊可靠性問題轉(zhuǎn)化為常規(guī)可靠性問題，并且使用數(shù)值的方法來求解模糊失效概率，但是其計算成本過大。李貴杰和宋軍等[5-6]將隨機變量積分區(qū)域離散化，使得每個積分區(qū)域內(nèi)功能函數(shù)對模糊失效區(qū)域隸屬函數(shù)的取值近似相等，從而將模糊失效概率轉(zhuǎn)換成常規(guī)的失效概率，然后提出基于矩方法和靶點逼近的模糊可靠性算法。Feng等[7]將模糊失效概率等價轉(zhuǎn)換成一系列常規(guī)失效概率，并使用子集模擬法來重復計算一組總樣本對應的一系列常規(guī)失效概率，進而求解模糊失效概率。Yun等[8]利用Kriging模型對樣本空間進行分層，對模糊狀態(tài)區(qū)域內(nèi)的點進行精確運算，從而有效地計算模糊失效概率。Zhang等[9]提出一種單循環(huán)Monte Carlo的方法來有效地估計模糊失效狀態(tài)下局部和全局可靠性靈敏度。

本文對渦輪盤進行疲勞可靠性分析，在此基礎上考慮渦輪盤失效邊界的模糊性，建立具有模糊失效邊界的渦輪盤疲勞可靠性分析模型。再利用積分轉(zhuǎn)換和高斯求積公式，將模糊失效概率轉(zhuǎn)化成一系列失效概率的加權和，并基于此建立模糊失效概率的全局靈敏度。最后利用改進的主動學習Kriging和Monte Carlo模擬相結(jié)合的方法(AK-MC法)同時近似多個失效面，提升計算效率。

1 渦輪盤低循環(huán)疲勞壽命可靠性分析

1.1 渦輪盤載荷簡化分析

以某航空發(fā)動機渦輪盤簡化模型為研究對象，依據(jù)實際結(jié)構的尺寸，根據(jù)渦輪盤循環(huán)對稱特性，選取渦輪盤的1/36部分作為研究對象。本文主要考慮渦輪盤運行期間的離心力和溫度載荷。其中離心力主要包括渦輪盤本身的離心力和渦輪葉片對渦輪盤的離心力。 在有限元分析過程中,渦輪葉片對渦輪盤的離心力負載在渦輪盤外側(cè)連接面上。

F=mrω2

(1)

式中:m是渦輪葉片的質(zhì)量；ω表示不同飛行任務下渦輪盤的轉(zhuǎn)速,在最大、巡航、空轉(zhuǎn)下分別為1 200,1 100,600 rad/s。

實際飛行中,不同任務會導致發(fā)動機的工作狀態(tài)不同。本文選定的飛機發(fā)動機載荷譜周期為1 000 h,其環(huán)境混頻后的轉(zhuǎn)速循環(huán)壓縮譜如表1所示。

表1 轉(zhuǎn)速循環(huán)壓縮譜

此外,渦輪盤通常在高溫環(huán)境下工作,盤體上的溫度梯度會產(chǎn)生極大的溫度載荷。在本文分析模型中,只有渦輪盤心和盤緣處的溫度數(shù)據(jù),盤體的溫度分布不明確,一般使用經(jīng)驗公式對渦輪盤體的溫度進行插值。常見的經(jīng)驗公式有線性插值和非線性插值,本文簡化了溫度載荷,使用線性插值來擬合渦輪盤體的溫度分布。

T=T0+KΔT

(2)

式中:T0為盤心溫度400℃;ΔT為盤心盤緣溫度差150℃;K為線性系數(shù)。渦輪盤在最大、巡航、空轉(zhuǎn)下的應力應變分布如圖1所示,最危險位置位于盤心處。

圖1 3種工作狀態(tài)對應的應力應變響應

1.2 渦輪盤低循環(huán)疲勞壽命模型

渦輪盤的失效模式復雜,其中低循環(huán)疲勞失效尤為顯著。本文考慮帶Morrow彈性應力線性修正的Manson-Coffin方程來計算疲勞壽命

(3)

渦輪盤在工作時會受到多級循環(huán)載荷的作用,一般使用Miner線性累積損傷法則來計算累積損傷,渦輪盤的累積損傷可以通過(4)式來計算

(4)

式中：k是循環(huán)載荷級數(shù)；ni是第i級載荷實際循環(huán)數(shù)；Nfi是第i級載荷下對應的疲勞壽命；D是多級載荷下的總損傷。

通常當累積損傷量達到損傷極限時,疲勞破壞發(fā)生。以此建立線性累積損傷模型的極限狀態(tài)方程,如(5)式所示

(5)

式中，a為疲勞破壞參數(shù),一般情況下取值為1。

1.3 渦輪盤低循環(huán)疲勞可靠性

對于在惡劣工況下服役的航空發(fā)動機而言,其材料性能和載荷的分散性會導致疲勞壽命存在分散性。在對渦輪盤進行疲勞壽命分析的基礎上,綜合考慮渦輪盤材料性能和載荷的不確定性,進而可以定量地描述渦輪盤的疲勞可靠性。本文主要考慮材料性能和載荷的不確定性[10],其分布信息如表2所示。

表2 隨機變量

根據(jù)(5)式的極限狀態(tài)函數(shù),渦輪盤的失效概率可以表示為

(6)

對于渦輪盤的可靠性分析問題,采用傳統(tǒng)抽樣方法計算效率低下。本文主要采用自適應代理模型法(active learning Kriging)來進行計算,該方法已經(jīng)過許多文獻驗證[11-13],可以通過較小的計算代價得到精度較高的可靠性計算結(jié)果。因此,當設計壽命Tt=15 000 h時,利用AK-MC的方法只需要調(diào)用166次有限元分析結(jié)果就可以計算出渦輪盤的失效概率,失效概率計算結(jié)果Pf=0.008 56。

2 具有模糊失效狀態(tài)的渦輪盤疲勞可靠性及靈敏度

2.1 具有模糊失效狀態(tài)的渦輪盤疲勞可靠性

在結(jié)構疲勞壽命分析中,當疲勞累積損傷大于疲勞破壞參數(shù)時,疲勞失效發(fā)生。但實際過程中由于載荷加載順序以及小載荷強化損傷等因素的影響[14],疲勞破壞參數(shù)取值可能存在擾動,這將引起渦輪盤的安全域與失效域之間存在一個模糊的過渡邊界,導致載荷循環(huán)未達到預定疲勞壽命時也可能發(fā)生失效。為此本文考慮疲勞破壞參數(shù)具有模糊性,則(6)式中失效域指示函數(shù)可以表示為模糊指示函數(shù)(即,失效隸屬度)

(7)

(8)

(9)

(10)

對于(9)式中的一維積分問題,可以利用高斯求積公式,將模糊失效概率轉(zhuǎn)化成一系列失效概率的加權平均

(11)

式中，ωj和λj分別是高斯積分權重和高斯積分點。

可以看到,對于渦輪盤模糊失效概率的計算需要大量調(diào)用渦輪盤的有限元分析結(jié)果,因此傳統(tǒng)的基于抽樣的可靠性分析方法難以有效計算,需要建立高效高精度的模糊可靠性分析方法。

2.2 基于模糊失效概率的渦輪盤全局靈敏度

基于失效概率的全局靈敏度可以衡量輸入變量在其整個分布域內(nèi)變化時對結(jié)構失效概率的平均影響,能夠識別影響失效概率的關鍵因素。根據(jù)文獻[15]中的全局靈敏度指標定義,本文建立(8)式中的模糊失效概率全局靈敏度指標為

(12)

(13)

根據(jù)(9)式,條件模糊失效概率也可以轉(zhuǎn)換為一系列條件失效概率對于隸屬度的積分

(14)

進而,條件失效概率可以表示為高斯積分的形式

(15)

3 基于AK-MC的模糊失效概率及靈敏度分析方法

3.1 基于AK-MC的模糊失效概率分析方法

(16)

式中,S表示抽取的樣本池。對于所建立的Kriging模型,可以通過選取樣本池中U學習函數(shù)最小的樣本點對模型進行更新。根據(jù)文獻[20]中的收斂準則,通常取Kriging模型更新的收斂條件為:minU(x)≥2。

(17)

(18)

Kriging模型更新完成后,可以利用該模型對渦輪盤的響應進行估計,進而求解渦輪盤的模糊失效概率為

(19)

(20)

利用AK-MC對模糊失效概率進行分析,只需要調(diào)用渦輪盤的真實有限元分析結(jié)果來建立和更新Kriging模型,因此可以高效地得到計算結(jié)果。

3.2 基于AK-MC的模糊失效概率靈敏度分析方法

(21)

根據(jù)(12)式,模糊失效概率的全局靈敏度指標可以估計為

(22)

式中，N1為Xi的樣本個數(shù),即(m=1,2,…N1)。

對于模糊失效概率的靈敏度計算,由于只是利用Kriging模型來估計條件模糊失效概率,沒有任何附加的有限元分析過程,因此其計算量和估計模糊失效概率的計算量相同。

4 分析及討論

4.1 渦輪盤模糊可靠性分析

針對渦輪盤的模糊失效概率,當設計壽命Tt為15 000 h時,采用上述改進的AK-MC可以對其進行分析計算?？紤]到渦輪盤的模糊失效概率會隨著疲勞破壞參數(shù)擾動Δ的變化而變化,因此本文計算得到不同擾動下的模糊失效概率結(jié)果如表3所示。

表3 模糊失效概率結(jié)果對比

表3中的結(jié)果表明:與傳統(tǒng)的失效概率(Pf=0.008 56)相比,渦輪盤的模糊失效概率相對較大,并且隨著疲勞破壞參數(shù)擾動Δ的增大而增大。這表明未考慮渦輪盤失效邊界模糊性的可靠性分析結(jié)果是偏危險的,并且隨著模糊邊界的擴大而增加。因此在開展渦輪盤可靠性分析時,應該考慮失效邊界的擾動情況,這樣得到的可靠性分析結(jié)果更為保守。

此外,由表3的結(jié)果可以看出,AK-MC方法只需要調(diào)用數(shù)百次的有限元分析結(jié)果便可建立Kriging模型,進而估算出渦輪盤的模糊失效概率,并且當樣本池取105時,模糊失效概率的變異系數(shù)都小于5%,這說明數(shù)值計算的結(jié)果已經(jīng)收斂。同時,當擾動參數(shù)增加時,Kriging模型所需的更新樣本點也隨之增加,這也說明Kriging模型需要在更廣的失效邊界上獲得較好的預測精度。

4.2 渦輪盤靈敏度分析

針對(22)式中給出的模糊失效概率全局靈敏度,可以利用改進的AK-MC對其進行分析,并且無需額外的有限元分析過程。表4中給出了疲勞破壞參數(shù)擾動Δ從0.1變化到0.5時渦輪盤基于模糊失效的全局靈敏度指標計算結(jié)果。為了直觀表達,全局靈敏度指標分析結(jié)果對比如圖2所示。

圖2 全局靈敏度指標對比圖

表4 全局靈敏度指標計算結(jié)果

5 結(jié) 論

本文針對某GH4169粉末合金渦輪盤有限元模型，開展其疲勞可靠性分析研究?？紤]渦輪盤失效邊界的模糊性，建立具有模糊失效邊界的渦輪盤疲勞可靠性模型，再利用積分轉(zhuǎn)換和高斯求積公式，將模糊失效概率轉(zhuǎn)化成一系列失效概率的加權和，并基于此建立模糊失效概率的全局靈敏度。

為了高效計算模糊失效概率及其靈敏度，本文利用改進的AK-MC法同時近似多個失效面，進而提升計算效率。與傳統(tǒng)失效概率計算結(jié)果進行對比，模糊失效概率相對較大，這說明傳統(tǒng)可靠性分析結(jié)果偏危險，并且隨著疲勞破壞參數(shù)的擾動 的增大而逐漸增大。此外，全局靈敏度的計算結(jié)果識別出影響渦輪盤結(jié)構安全性的3個關鍵參數(shù)。

盡管本文中考慮的失效邊界的模糊性，但對具體的擾動情況只是粗略給出其范圍并研究其影響規(guī)律，對于具體的疲勞失效機理還需進行更為深入的研究。