劉 濤,鄧 強(qiáng),劉 源,龍 超
(西南交通大學(xué)機(jī)械工程學(xué)院,四川成都 610031)
葉盤是發(fā)動(dòng)機(jī)的核心部件,它的性能水平對(duì)發(fā)動(dòng)機(jī)的可靠性和安全性具有決定性的影響。隨著航空發(fā)動(dòng)機(jī)推重比的提高,導(dǎo)致葉盤的負(fù)荷越來越大,對(duì)葉盤工作溫度、轉(zhuǎn)速、性能的要求越來越高。隨著第四代戰(zhàn)斗機(jī)的出現(xiàn),開始采用整體葉盤結(jié)構(gòu),其優(yōu)點(diǎn)是省去了榫頭減小了輪緣載荷、有效地減輕了葉盤的重量,可以增加發(fā)動(dòng)機(jī)的推重比和改善冷卻效果,提高發(fā)動(dòng)機(jī)的效率[1-3]。
國(guó)內(nèi)外對(duì)發(fā)動(dòng)機(jī)葉盤強(qiáng)度問題展開了大量的研究工作,特別是美國(guó)在1980年后隨著HOST計(jì)劃的推進(jìn),在提高渦輪葉盤等熱端部件耐久性和可靠性方面取得了豐碩成果[4-5];我國(guó)在 80年代后對(duì)發(fā)動(dòng)機(jī)熱端部件展開的研究工作,也取得了一定的成果[6-8]?,F(xiàn)如今隨著對(duì)發(fā)動(dòng)機(jī)推重比要求的提高,渦輪葉盤研究的一大趨勢(shì)是既要保證葉盤的高強(qiáng)度又要降低葉盤的質(zhì)量[1]?;诖?,筆者以質(zhì)量更輕的整體葉盤為研究對(duì)象,對(duì)兩種不同斷面結(jié)構(gòu)的整體葉盤的強(qiáng)度進(jìn)行綜合分析,該計(jì)算對(duì)整體葉盤的結(jié)構(gòu)設(shè)計(jì)具有借鑒意義。
基于預(yù)先研究的某型航空發(fā)動(dòng)機(jī)上的整體葉盤,該整體葉盤上共有30片葉片,省去了榫頭連接形式,葉片通過先進(jìn)焊接技術(shù)與渦輪盤連接在一起,渦輪盤中心開有一個(gè)孔,作為冷卻氣流流通用。忽略葉—盤焊接處的局部影響,在UG中建立實(shí)心葉盤模型,如圖1(a)所示。將UG中建立的葉盤模型導(dǎo)入ABAQUS中,采用四節(jié)點(diǎn)線性四面體單元對(duì)模型進(jìn)行離散,得到葉盤的有限元模型,如圖1(b)所示。由于葉盤模型形狀不規(guī)則需要手工劃分網(wǎng)格,來防止所分網(wǎng)格包含奇異角單元,在計(jì)算過程中造成剛度矩陣奇異,導(dǎo)致計(jì)算失敗。
按質(zhì)量更輕準(zhǔn)則[1],對(duì)實(shí)心葉盤的斷面進(jìn)行優(yōu)化,得到空心葉盤,優(yōu)化后的空心葉盤和實(shí)心葉盤的斷面結(jié)構(gòu)對(duì)比,如圖2所示。
圖1 整體葉盤計(jì)算模型
圖2 實(shí)心、空心葉盤斷面結(jié)構(gòu)
葉盤的材料是某高溫合金,它是以鎳-鉻為機(jī)體并添加其它元素形成彌散強(qiáng)化的合金,在高溫下具有足夠的強(qiáng)度[9],在各溫度下材料的彈性模量、泊松比、線膨脹系數(shù)、強(qiáng)度極限和密度等性能參數(shù),如表1。
表1 某型高溫葉盤材料力學(xué)性能參數(shù)
(1)邊界條件 采用整體葉盤進(jìn)行分析計(jì)算時(shí),整個(gè)葉盤受到剛體約束,為防止葉盤軸向發(fā)生剛體位移,還需要對(duì)葉盤中心孔處的節(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ù)比較有限,對(duì)葉盤強(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 溫度載荷分布
單溫度場(chǎng)作用下實(shí)心整體葉盤和空心整體葉盤的等效應(yīng)力分布如圖5所示。
計(jì)算結(jié)果表明,在單溫度場(chǎng)作用下,當(dāng)葉緣溫度高于盤心溫度時(shí),由溫度梯度引起的熱應(yīng)力,實(shí)心整體葉盤和空心整體葉盤均沿半徑方向由外向內(nèi)逐漸擴(kuò)大,并在中心孔某位置處達(dá)到最大值,實(shí)心葉盤和空心葉盤等效熱應(yīng)力最大值分別為959.2 MPa和1 000 MPa。
圖5 單溫度場(chǎng)作用下的等效應(yīng)力分布云圖
單離心力作用下實(shí)心整體葉盤和空心整體葉盤的等效應(yīng)力分布如圖6所示。
圖6 單離心載荷作用下的等效應(yīng)力分布云圖
計(jì)算結(jié)果表明,在單離心載荷作用下實(shí)心整體葉盤和空心整體葉盤,沿半徑方向的等效應(yīng)力變化趨勢(shì)相同,均出現(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。
在離心力和溫度場(chǎng)耦合作用下,實(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)度要求,該處對(duì)焊接技術(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)度分析,對(duì)于空心葉盤還須注意空腔部位的強(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)更為合理。
(1)整體葉盤應(yīng)力水平低,能夠滿足葉盤強(qiáng)度要求;與實(shí)心整體葉盤相比,質(zhì)量更輕的空心葉盤性能更優(yōu),對(duì)提高發(fā)動(dòng)機(jī)的功率效果更加明顯。
(2)由于葉盤溫度分布不均勻,具有較大的溫度梯度,產(chǎn)生了較大的熱應(yīng)力,在計(jì)算過程中不可忽略溫度的影響。
(3)實(shí)心葉盤最大應(yīng)力集中發(fā)生在渦輪盤中心孔處;空心葉盤最大應(yīng)力集中發(fā)生在空腔部位及中心孔處,因此在設(shè)計(jì)過程中要重視對(duì)最大應(yīng)力集中部位的強(qiáng)度計(jì)算。
(4)與實(shí)心葉盤相比,結(jié)構(gòu)優(yōu)化后,質(zhì)量更輕的空心葉盤,在同一工作載荷下能夠滿足工作強(qiáng)度的要求,且在中心孔附近區(qū)域的應(yīng)力值也有所下降,提高了葉盤中心孔這一應(yīng)力最大區(qū)域的可靠性。
[1] 王 營(yíng),余朝蓬.航空發(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一級(jí)發(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] 中國(guó)航空材料手冊(cè)編輯委員會(huì).中國(guó)航空材料手冊(cè)[M].北京:中國(guó)標(biāo)準(zhǔn)出版社,2002.