某型航空發(fā)動(dòng)機(jī)整體葉盤強(qiáng)度分析*

劉 濤，鄧 強(qiáng)，劉 源，龍 超

(西南交通大學(xué)機(jī)械工程學(xué)院，四川成都 610031)

0 引言

葉盤是發(fā)動(dòng)機(jī)的核心部件，它的性能水平對發(fā)動(dòng)機(jī)的可靠性和安全性具有決定性的影響。隨著航空發(fā)動(dòng)機(jī)推重比的提高，導(dǎo)致葉盤的負(fù)荷越來越大，對葉盤工作溫度、轉(zhuǎn)速、性能的要求越來越高。隨著第四代戰(zhàn)斗機(jī)的出現(xiàn)，開始采用整體葉盤結(jié)構(gòu)，其優(yōu)點(diǎn)是省去了榫頭減小了輪緣載荷、有效地減輕了葉盤的重量，可以增加發(fā)動(dòng)機(jī)的推重比和改善冷卻效果，提高發(fā)動(dòng)機(jī)的效率［1－3］。

國內(nèi)外對發(fā)動(dòng)機(jī)葉盤強(qiáng)度問題展開了大量的研究工作，特別是美國在1980年后隨著HOST計(jì)劃的推進(jìn)，在提高渦輪葉盤等熱端部件耐久性和可靠性方面取得了豐碩成果［4－5］;我國在 80年代后對發(fā)動(dòng)機(jī)熱端部件展開的研究工作，也取得了一定的成果［6－8］?，F(xiàn)如今隨著對發(fā)動(dòng)機(jī)推重比要求的提高，渦輪葉盤研究的一大趨勢是既要保證葉盤的高強(qiáng)度又要降低葉盤的質(zhì)量［1］?；诖耍P者以質(zhì)量更輕的整體葉盤為研究對象，對兩種不同斷面結(jié)構(gòu)的整體葉盤的強(qiáng)度進(jìn)行綜合分析，該計(jì)算對整體葉盤的結(jié)構(gòu)設(shè)計(jì)具有借鑒意義。

1 計(jì)算模型的建立

1.1 葉盤的結(jié)構(gòu)特點(diǎn)

基于預(yù)先研究的某型航空發(fā)動(dòng)機(jī)上的整體葉盤，該整體葉盤上共有30片葉片，省去了榫頭連接形式，葉片通過先進(jìn)焊接技術(shù)與渦輪盤連接在一起，渦輪盤中心開有一個(gè)孔，作為冷卻氣流流通用。忽略葉—盤焊接處的局部影響，在UG中建立實(shí)心葉盤模型，如圖1(a)所示。將UG中建立的葉盤模型導(dǎo)入ABAQUS中，采用四節(jié)點(diǎn)線性四面體單元對模型進(jìn)行離散，得到葉盤的有限元模型，如圖1(b)所示。由于葉盤模型形狀不規(guī)則需要手工劃分網(wǎng)格，來防止所分網(wǎng)格包含奇異角單元，在計(jì)算過程中造成剛度矩陣奇異，導(dǎo)致計(jì)算失敗。

按質(zhì)量更輕準(zhǔn)則［1］，對實(shí)心葉盤的斷面進(jìn)行優(yōu)化，得到空心葉盤，優(yōu)化后的空心葉盤和實(shí)心葉盤的斷面結(jié)構(gòu)對比，如圖2所示。

圖1 整體葉盤計(jì)算模型

圖2 實(shí)心、空心葉盤斷面結(jié)構(gòu)

1.2 葉盤的材料參數(shù)

葉盤的材料是某高溫合金，它是以鎳－鉻為機(jī)體并添加其它元素形成彌散強(qiáng)化的合金，在高溫下具有足夠的強(qiáng)度［9］，在各溫度下材料的彈性模量、泊松比、線膨脹系數(shù)、強(qiáng)度極限和密度等性能參數(shù)，如表1。

表1 某型高溫葉盤材料力學(xué)性能參數(shù)

1.3 邊界條件與計(jì)算載荷

(1)邊界條件 采用整體葉盤進(jìn)行分析計(jì)算時(shí)，整個(gè)葉盤受到剛體約束，為防止葉盤軸向發(fā)生剛體位移，還需要對葉盤中心孔處的節(jié)點(diǎn)進(jìn)行軸向約束。

(2)載荷施加 發(fā)動(dòng)機(jī)工作時(shí)，葉盤承受的載荷主要有自身的離心力載荷、溫度梯度載荷，氣動(dòng)載荷、葉盤振動(dòng)時(shí)產(chǎn)生的振動(dòng)載荷［3］。而其中氣動(dòng)載荷和振動(dòng)載荷數(shù)據(jù)比較有限，對葉盤強(qiáng)度影響不大，所以在葉盤強(qiáng)度計(jì)算時(shí)，主要考慮離心載荷和溫度載荷的影響。

本型整體葉盤的設(shè)計(jì)轉(zhuǎn)速為10 000 r/min(即1 046.67 rad/s)，葉盤自身離心載荷以角速度的形式施加在整個(gè)計(jì)算模型上，加載后的模型離心載荷分布，如圖3所示。

通過經(jīng)驗(yàn)公式與實(shí)驗(yàn)數(shù)據(jù)相擬合，得到與半徑相關(guān)的溫度梯度載荷表達(dá)式為:

式中:r為葉盤半徑，其溫度載荷分布如圖4所示。

圖3 離心載荷分布

圖4 溫度載荷分布

2 葉盤應(yīng)力計(jì)算與分析

2.1 單溫度載荷作用下的應(yīng)力計(jì)算

單溫度場作用下實(shí)心整體葉盤和空心整體葉盤的等效應(yīng)力分布如圖5所示。

計(jì)算結(jié)果表明，在單溫度場作用下，當(dāng)葉緣溫度高于盤心溫度時(shí)，由溫度梯度引起的熱應(yīng)力，實(shí)心整體葉盤和空心整體葉盤均沿半徑方向由外向內(nèi)逐漸擴(kuò)大，并在中心孔某位置處達(dá)到最大值，實(shí)心葉盤和空心葉盤等效熱應(yīng)力最大值分別為959.2 MPa和1 000 MPa。

圖5 單溫度場作用下的等效應(yīng)力分布云圖

2.2 單離心載荷作用下的應(yīng)力計(jì)算

單離心力作用下實(shí)心整體葉盤和空心整體葉盤的等效應(yīng)力分布如圖6所示。

圖6 單離心載荷作用下的等效應(yīng)力分布云圖

計(jì)算結(jié)果表明，在單離心載荷作用下實(shí)心整體葉盤和空心整體葉盤，沿半徑方向的等效應(yīng)力變化趨勢相同，均出現(xiàn)峰值波動(dòng)，在葉緣和中心孔處應(yīng)力值最低，應(yīng)力沿半徑方向由外向內(nèi)，在葉片部位應(yīng)力逐漸增大，到盤體部位在逐漸減小后又增大，到中心孔附近再逐漸減小。不同的是實(shí)心葉盤應(yīng)力最大值出現(xiàn)在葉－盤結(jié)合處，大小為565.3 MPa，而空心葉盤的應(yīng)力最大值出現(xiàn)在空腔位置處，大小為791.2 MPa。

2.3 離心力和溫度場耦合作用下的應(yīng)力計(jì)算

在離心力和溫度場耦合作用下，實(shí)心葉盤和空心葉盤的等效應(yīng)力分布如圖7所示。應(yīng)力－半徑變化曲線如圖8所示。

圖7 等效應(yīng)力分布云圖

圖8 應(yīng)力－半徑變化曲線圖

計(jì)算結(jié)果表明，無論是實(shí)心葉盤還是空心葉盤，葉緣處的應(yīng)力水平都是最低的，而在葉－盤結(jié)合處，應(yīng)力具有較高的水平，在制造時(shí)應(yīng)注意葉－盤結(jié)合部的強(qiáng)度要求，該處對焊接技術(shù)要求極高，因?yàn)榧纫獫M足接合部位的強(qiáng)度要求，又要保證焊接過程中不發(fā)生重熔現(xiàn)象而影響壽命，這正是整體葉盤相較于傳統(tǒng)渦輪制造工藝復(fù)雜的原因。實(shí)心整體葉盤的最大應(yīng)力發(fā)生在中心孔處，其最大應(yīng)力值1 056 MPa?？招恼w葉盤的最大應(yīng)力發(fā)生在中心孔和空腔部位，其最大應(yīng)力值為1 081 MPa，所以在設(shè)計(jì)過程中要特別重視中心孔處的強(qiáng)度分析，對于空心葉盤還須注意空腔部位的強(qiáng)度需求。

從應(yīng)力－半徑變化曲線圖可以看出，實(shí)心整體葉盤和空心整體葉盤，均滿足工作溫度下合金材料屈服極限和強(qiáng)度極限的要求。結(jié)構(gòu)優(yōu)化后的空心葉盤與實(shí)心葉盤相比，除了質(zhì)量更輕之外，在中心孔附近區(qū)域的應(yīng)力值也有所下降，提高了葉盤中心孔這一應(yīng)力最大區(qū)域的可靠性，進(jìn)一步表明了優(yōu)化后空心葉盤的結(jié)構(gòu)更為合理。

3 結(jié)論

(1)整體葉盤應(yīng)力水平低，能夠滿足葉盤強(qiáng)度要求;與實(shí)心整體葉盤相比，質(zhì)量更輕的空心葉盤性能更優(yōu)，對提高發(fā)動(dòng)機(jī)的功率效果更加明顯。

(2)由于葉盤溫度分布不均勻，具有較大的溫度梯度，產(chǎn)生了較大的熱應(yīng)力，在計(jì)算過程中不可忽略溫度的影響。

(3)實(shí)心葉盤最大應(yīng)力集中發(fā)生在渦輪盤中心孔處;空心葉盤最大應(yīng)力集中發(fā)生在空腔部位及中心孔處，因此在設(shè)計(jì)過程中要重視對最大應(yīng)力集中部位的強(qiáng)度計(jì)算。

(4)與實(shí)心葉盤相比，結(jié)構(gòu)優(yōu)化后，質(zhì)量更輕的空心葉盤，在同一工作載荷下能夠滿足工作強(qiáng)度的要求，且在中心孔附近區(qū)域的應(yīng)力值也有所下降，提高了葉盤中心孔這一應(yīng)力最大區(qū)域的可靠性。

［1］ 王 營，余朝蓬.航空發(fā)動(dòng)機(jī)渦輪盤結(jié)構(gòu)優(yōu)化設(shè)計(jì)［J］.機(jī)械設(shè)計(jì)與制造，2010(5):4－6.

［2］ 李水姣.某型燃?xì)廨啓C(jī)低壓渦輪盤強(qiáng)度分析［J］.機(jī)械工程與自動(dòng)化，2012，5(174):37－39.

［3］ 趙 霞.某渦漿發(fā)動(dòng)機(jī)渦輪轉(zhuǎn)子強(qiáng)度與壽命計(jì)算［D］.南京:南京航空航天大學(xué)，2005.

［4］ Behrooz Farshi，Hamid Jahed，Abdochossein Mehrabian.Optimimum Design of Inhomogeneous Non－uniform Rotating Discs［J］.Computers and structures，2004(82):773－779.

［5］ Koop W.The Integrated High Performance Turbine Engine Technology(IHPTET)Program［R］.ISABE 97－7175.

［6］ 王占學(xué)，劉增文，蔡元虎，等.推重比15一級發(fā)動(dòng)機(jī)關(guān)鍵技術(shù)及分析［J］.航空發(fā)動(dòng)機(jī)，2010，36(1):58－62.

［7］ 陳 光.航空發(fā)動(dòng)機(jī)結(jié)構(gòu)設(shè)計(jì)分析［M］.北京:北京航空航天大學(xué)出版社，2006.

［8］ 梁春華.通用的經(jīng)濟(jì)可承受的先進(jìn)渦輪發(fā)動(dòng)機(jī)研究計(jì)劃的主要特點(diǎn)［J］.航空發(fā)動(dòng)機(jī)，2011，37(5):58－62.

［9］ 中國航空材料手冊編輯委員會(huì).中國航空材料手冊［M］.北京:中國標(biāo)準(zhǔn)出版社，2002.