CFM56-7B發(fā)動(dòng)機(jī)高壓渦輪葉片疲勞壽命預(yù)測(cè)

曹惠玲 張昊

摘要：基于CFM56-7B發(fā)動(dòng)機(jī)運(yùn)行產(chǎn)生的QAR數(shù)據(jù)，通過(guò)建立模型、確立邊界條件、有限元軟件仿真，在最大起飛功率狀態(tài)下完成了高壓渦輪葉片溫度場(chǎng)以及應(yīng)力、應(yīng)變場(chǎng)的計(jì)算?？紤]到實(shí)際葉片在多軸非比例循環(huán)載荷下工作，選用SWT模型作為疲勞壽命的預(yù)測(cè)模型，得到葉片的疲勞壽命。結(jié)果表明，仿真計(jì)算得到的葉片疲勞壽命為14778個(gè)循環(huán)，與實(shí)際葉片的平均壽命的誤差率僅為5.09%，此方法可用于監(jiān)控渦輪葉片的剩余壽命、維修計(jì)劃的制定。

關(guān)鍵詞：QAR數(shù)據(jù)；高壓渦輪葉片；有限元；SWT模型；疲勞壽命

Keywords：QAR data；high pressure turbine blade；finite element；SWT model；fatigue life

0 引言

高壓渦輪葉片是航空發(fā)動(dòng)機(jī)里最重要的部件之一，其工作環(huán)境極其惡劣，一方面承受著高溫、高壓燃?xì)獾淖饔茫硪环矫媸艿秸駝?dòng)等多種載荷的作用。根據(jù)航空工業(yè)集團(tuán)失效分析中心發(fā)布的資料[1]顯示，2001年到2012年期間共處理97件有關(guān)渦輪葉片的失效故障，占整個(gè)失效故障分析工作的1/10左右，且據(jù)實(shí)際統(tǒng)計(jì)，在由發(fā)動(dòng)機(jī)導(dǎo)致的飛行事故中渦輪葉片失效占比達(dá)70%以上[2]，因此進(jìn)行高壓渦輪葉片的結(jié)構(gòu)強(qiáng)度分析及壽命預(yù)測(cè)具有重要意義。

快速存取記錄器（QAR）記錄了發(fā)動(dòng)機(jī)從啟動(dòng)到停車(chē)過(guò)程的大量氣動(dòng)站位參數(shù)，這些參數(shù)真實(shí)地反映了發(fā)動(dòng)機(jī)實(shí)際的運(yùn)行工況，因此將QAR數(shù)據(jù)應(yīng)用到高壓渦輪葉片的結(jié)構(gòu)強(qiáng)度分析和壽命預(yù)測(cè)中能夠得到更加精確的結(jié)果。以國(guó)內(nèi)某航空公司某航班CFM56-7B發(fā)動(dòng)機(jī)運(yùn)行產(chǎn)生的QAR數(shù)據(jù)為例，通過(guò)飛行高度、高壓渦輪轉(zhuǎn)速、EGT（發(fā)動(dòng)機(jī)排氣溫度）以及燃油流量篩選出起飛最大功率狀態(tài)，該狀態(tài)下高壓渦輪轉(zhuǎn)速高達(dá) 97.6%，EGT達(dá)到804℃，燃油流量達(dá)到3607kg/h，對(duì)比其他工況，該狀態(tài)下的葉片承受著更高的溫度以及離心力作用，因此以此狀態(tài)作為高壓渦輪葉片壽命計(jì)算的工況。在起飛最大功率狀態(tài)下，利用QAR記錄的高壓壓氣機(jī)進(jìn)出口總溫、EGT、低壓壓氣機(jī)進(jìn)口總壓等參數(shù)，結(jié)合熱力計(jì)算可得到用于葉片壽命預(yù)測(cè)的邊界條件。

本文利用有限元仿真計(jì)算的方法，通過(guò)建立模型和計(jì)算邊界條件，使用Fluent軟件求解燃?xì)饬鲌?chǎng)以及Ansys軟件進(jìn)行熱分析和結(jié)構(gòu)強(qiáng)度分析，得到葉片的溫度、應(yīng)力應(yīng)變分布，綜合考慮對(duì)葉片影響較大的疲勞損傷，最終計(jì)算得到葉片的疲勞壽命。

1 有限元仿真計(jì)算

1.1 模型的建立

三維建模是有限元數(shù)值仿真的基礎(chǔ)，精確的幾何模型能夠極大地提高有限元仿真的精度，本文針對(duì)CFM56-7B高壓渦輪葉片結(jié)構(gòu)，參照航空發(fā)動(dòng)機(jī)設(shè)計(jì)手冊(cè)—渦輪分冊(cè)[3]，通過(guò)實(shí)體結(jié)構(gòu)的測(cè)量，利用UG完成了三維模型的建立，具體模型如圖1、圖2所示。根據(jù)建立的高壓渦輪葉片模型，選取中截面葉型、葉片高度、葉柵柵距等參數(shù)建立圖3所示簡(jiǎn)化的光滑葉片葉柵通道模型。利用ICEM CFD完成模型的非結(jié)構(gòu)網(wǎng)格劃分，葉片模型的網(wǎng)格數(shù)為35萬(wàn)，葉柵通道的網(wǎng)格數(shù)為196萬(wàn)，網(wǎng)格質(zhì)量均在0.2以上且進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證。

1.2 邊界條件的計(jì)算

QAR記錄的參數(shù)包括高壓壓氣機(jī)進(jìn)出口總溫和、EGT、低壓壓氣機(jī)進(jìn)口總壓、燃油流量以及油氣比f(wàn)，這些參數(shù)在起飛最大功率狀態(tài)下的數(shù)值如表1所示。在進(jìn)行有限元仿真計(jì)算的過(guò)程中需要首先確定高壓壓氣機(jī)出口的總溫、高壓渦輪進(jìn)口的總溫、總壓以及高壓渦輪出口的總溫、靜壓。高壓壓氣機(jī)出口的總溫可由QAR數(shù)據(jù)直接獲取，高壓渦輪進(jìn)出口的溫度、壓力可通過(guò)表1所給出的QAR數(shù)據(jù)結(jié)合熱力計(jì)算的方式得到，具體計(jì)算方法參見(jiàn)文獻(xiàn)[4]。根據(jù)計(jì)算結(jié)果可得到有限元仿真計(jì)算所需要的邊界參數(shù)，具體數(shù)值如表2所示。

CFM56-7B高壓渦輪葉片由高壓壓氣機(jī)排氣進(jìn)行冷卻，高壓壓氣機(jī)排氣從葉根進(jìn)入葉片內(nèi)部的冷卻通道，首先對(duì)轉(zhuǎn)子葉片的根部進(jìn)行冷卻，然后氣流在葉片的各冷卻腔室內(nèi)迂回流動(dòng)，對(duì)葉片進(jìn)行充分的冷卻，最后冷卻氣體從葉片表面的氣膜孔以及尾緣劈縫流出，在葉片表面形成氣膜，隔絕高溫燃?xì)馀c高壓渦輪葉片。為表征冷卻氣體與葉片的對(duì)流換熱情況，引入對(duì)流換熱系數(shù)，冷卻氣體與葉片的平均對(duì)流換熱系數(shù)取29124W/K·m2。

1.3 仿真計(jì)算結(jié)果

根據(jù)得到的有限元模型以及邊界條件，將其導(dǎo)入Fluent中完成流場(chǎng)的計(jì)算，采用壓力進(jìn)口、壓力出口邊界條件，選取k-SST湍流計(jì)算模型。為了使計(jì)算更加符合實(shí)際情況，考慮渦輪進(jìn)口溫度的不均勻性[5]，取進(jìn)口溫度沿徑向按拋物線分布，并編制UDF導(dǎo)入Fluent進(jìn)行求解計(jì)算，圖4為溫度不均勻時(shí)流場(chǎng)進(jìn)口溫度的分布。

將Fluent計(jì)算得到的數(shù)據(jù)導(dǎo)入Ansys中進(jìn)行熱固耦合的計(jì)算，葉片材料為定向凝固高溫鎳基合金DZ125，在熱分析模塊中完成溫度場(chǎng)的計(jì)算；考慮離心載荷的影響，在結(jié)構(gòu)分析模塊中進(jìn)行應(yīng)力、應(yīng)變的計(jì)算。仿真流程如圖5所示。

通過(guò)仿真計(jì)算得到高壓渦輪葉片在起飛最大功率狀態(tài)下的溫度場(chǎng)以及應(yīng)力、應(yīng)變場(chǎng)，具體如圖6至圖8所示。

根據(jù)葉片溫度、應(yīng)力應(yīng)變的計(jì)算結(jié)果可以看出，在冷卻通道穿過(guò)的區(qū)域以及氣膜孔、尾緣劈縫處的溫度明顯更低，且葉片前緣溫度明顯高于尾緣，由于冷卻通道是從葉根區(qū)域穿入，葉片溫度的最大值出現(xiàn)在葉頂區(qū)域。葉片應(yīng)力、應(yīng)變值沿葉高逐漸升高，這是因?yàn)槿~片在高速轉(zhuǎn)動(dòng)過(guò)程中，葉身部分產(chǎn)生的離心力完全作用在葉根區(qū)域，導(dǎo)致葉根區(qū)域應(yīng)力、應(yīng)變值較大，而葉根區(qū)域氣膜孔處由于應(yīng)力集中導(dǎo)致葉片的最大應(yīng)力出現(xiàn)在該處。

2 疲勞壽命預(yù)測(cè)

根據(jù)得到的葉片溫度、應(yīng)力應(yīng)變的計(jì)算結(jié)果，在Ansys結(jié)構(gòu)分析模塊中插入疲勞壽命預(yù)測(cè)模塊，并選擇SWT模型，得到圖9所示的葉片疲勞壽命分布云圖。

由圖9可以看出，葉片疲勞壽命的最小值出現(xiàn)在靠近葉根的氣膜孔處，與葉片等效應(yīng)力的最大點(diǎn)為同一點(diǎn)，其值為14778循環(huán)。根據(jù)實(shí)際的外場(chǎng)統(tǒng)計(jì)數(shù)據(jù)可知，該型發(fā)動(dòng)機(jī)高壓渦輪葉片平均壽命為14062循環(huán)[7]，預(yù)測(cè)的壽命值與實(shí)際的葉片壽命比較接近，這表明真實(shí)QAR數(shù)據(jù)與有限元仿真相結(jié)合計(jì)算得到的預(yù)測(cè)壽命是吻合的。

3 結(jié)論

本文基于發(fā)動(dòng)機(jī)運(yùn)行產(chǎn)生的QAR數(shù)據(jù)，選取起飛最大功率狀態(tài)，利用有限元仿真的方法完成了高壓渦輪葉片流場(chǎng)、溫度場(chǎng)以及應(yīng)力應(yīng)變場(chǎng)的計(jì)算，根據(jù)得到的計(jì)算結(jié)果完成了疲勞壽命的預(yù)測(cè)。仿真計(jì)算得到的葉片疲勞壽命為14778循環(huán)，而葉片的實(shí)際平均壽命為14062循環(huán)，誤差率僅為5.09%，通過(guò)誤差分析驗(yàn)證了該方法的可行性，該方法可作為高壓渦輪葉片疲勞壽命預(yù)測(cè)的一種工程計(jì)算方法。

參考文獻(xiàn)

[1] 馬楠楠，陶春虎，何玉懷，等.航空發(fā)動(dòng)機(jī)葉片多軸疲勞試驗(yàn)研究進(jìn)展[J].航空材料學(xué)報(bào)，2012，32（6）：44-49.

[2] 蘇清友. 航空渦噴、渦扇發(fā)動(dòng)機(jī)主要零部件定壽指南[M].北京：航空工業(yè)版社，2004.

[3] 航空發(fā)動(dòng)機(jī)設(shè)計(jì)手冊(cè)總編委會(huì).航空發(fā)動(dòng)機(jī)設(shè)計(jì)手冊(cè)：第10冊(cè) 渦輪[M].北京：航空工業(yè)出版社，2001.

[4] 汪晗. HPT葉尖間隙的熱固耦合分析及對(duì)EGTM的影響研究[D].天津：中國(guó)民航大學(xué)，2019.

[5] 殷宇陽(yáng)，王宏光. 進(jìn)口溫度不均勻?qū)u輪葉片熱應(yīng)力的影響[J]. 能源研究與信息，2014，30（2）：113-117.

[6]《中國(guó)航空材料手冊(cè)》編輯委員會(huì). 中國(guó)航空材料手冊(cè)[M]. 北京：中國(guó)標(biāo)準(zhǔn)出版社，1992.

[7] 袁鍇，左洪福，孫見(jiàn)忠.基于統(tǒng)計(jì)方法的高壓渦輪轉(zhuǎn)子葉片壽命分析[J].飛機(jī)設(shè)計(jì)，2015，35（2）：31-36.