陳奕宏，李 維

(中國航發(fā)湖南動力機械研究所，湖南 株洲 412002)

0 引言

渦輪是將高溫高壓燃氣中的能量轉(zhuǎn)變?yōu)闄C械能的旋轉(zhuǎn)機械，是燃氣輪機的核心部件，其技術(shù)的先進程度直接與燃氣輪機性能相關(guān)。葉尖泄漏流是葉輪機械轉(zhuǎn)子中最普遍和有影響的流動現(xiàn)象之一，渦輪中大約三分之一的損失與葉尖泄漏流動有關(guān)[1]。葉尖小翼[2]作為一種能有效控制葉尖泄漏的方法，國內(nèi)外學(xué)者對葉尖小翼也開展了諸多研究，以尋找利用葉尖小翼控制葉尖泄漏的最佳途徑。但目前的研究主要集中在葉尖小翼結(jié)構(gòu)的試驗研究[3]和流場控制機理[4]研究，對葉尖小翼葉片的強度影響很少。

本文在設(shè)計點狀態(tài)，對某型發(fā)動機燃氣渦輪一級轉(zhuǎn)子葉片添加葉尖壓力面葉尖小翼、吸力面葉尖小翼、雙側(cè)葉尖小翼結(jié)構(gòu)，并與原始葉片進行定常數(shù)值研究，對比分析了未造型的原始葉片和造型葉片的應(yīng)力分布和位移變化，最后評估了葉尖小翼結(jié)構(gòu)對葉片強度的影響。

1 研究對象與數(shù)值計算方法

1.1 研究對象

本文以某型發(fā)動機燃氣渦輪一級轉(zhuǎn)子作為研究對象，以葉片徑向95%的位置作為造型起始點，對葉片葉頂部分進行周向造型，分別生成了壓力面小翼、吸力面小翼、雙側(cè)小翼三種帶葉尖小翼結(jié)構(gòu)的葉片模型，小翼的寬度統(tǒng)一設(shè)置為2 mm，其中雙側(cè)小翼采用與壓力面?zhèn)群臀γ鎮(zhèn)刃∫硐嗤慕Y(jié)構(gòu)拼接而成，各結(jié)構(gòu)葉片如圖1所示。

表1為不同葉尖結(jié)構(gòu)葉片的重量變化，其中壓力側(cè)小翼由于長度不及壓力面小翼，重量變化較小，雙側(cè)小翼相對于原始葉型重量增加了0.68%。

表1 不同葉尖結(jié)構(gòu)葉片質(zhì)量對比

圖1 不同葉尖結(jié)構(gòu)葉型示意圖

表2為不同葉尖結(jié)構(gòu)葉片的質(zhì)心變化，由于葉尖小翼結(jié)構(gòu)增加了葉片尖部的質(zhì)量，三種葉尖小翼結(jié)構(gòu)葉片質(zhì)心均有Z軸(徑向)的升高，且由于小翼結(jié)構(gòu)在葉片周向的發(fā)展，葉片質(zhì)心在X軸(發(fā)動機軸向)和Y軸(葉片軸向)都有相應(yīng)變化。

表2 不同葉尖結(jié)構(gòu)質(zhì)心相對變化

1.2 數(shù)值計算方法

1.2.1 壓力和溫度載荷計算

對葉片進行強度計算前需要對葉片的溫度載荷和壓力載荷進行計算，該部分計算通過商用軟件CFX求解雷諾平均Navier-Stokes(N-S)方程組實現(xiàn)。其中，湍流模型采用帶有自動壁面處理的剪切應(yīng)力輸運模型(SST模型)。湍流項和對流項的離散均采用高階格式。

計算域分為流體域和固體域兩個部分，計算網(wǎng)格采用icem非結(jié)構(gòu)化網(wǎng)格，固體域網(wǎng)格數(shù)量80萬左右，流體域在葉尖間隙和葉片前尾緣處進行了加密，邊界層第一層網(wǎng)格高度0.0025，層數(shù)設(shè)置為12層，網(wǎng)格y+控制在2附近，網(wǎng)格數(shù)量120萬左右，兩種計算域網(wǎng)格如圖2所示。

圖2 計算域網(wǎng)格

1.2.2 應(yīng)力計算

在取得了葉片的壓力載荷和溫度載荷后，即可以對葉片進行應(yīng)力計算分析。該部分計算采用Workbench軟件三維有限元方法對葉片局部應(yīng)力、變形和位移進行計算。

燃氣渦輪一級工作轉(zhuǎn)子的有限元網(wǎng)格見圖3，模型采用十節(jié)點四面體和二十節(jié)點六面體單元，共劃分有763947個單元，1157992個節(jié)點。

圖3 葉片有限元網(wǎng)格示意圖

2 葉尖小翼結(jié)構(gòu)對葉片強度的影響

2.1 葉尖小翼結(jié)構(gòu)對葉片最大應(yīng)力的影響

圖4和圖5分別為四種不同葉尖結(jié)構(gòu)葉型壓力面和吸力面的當量應(yīng)力分布圖，不同葉尖結(jié)構(gòu)吸力面壓力面的最大當量應(yīng)力值和兩側(cè)最大應(yīng)力差值如表3—表5所示。

圖4 不同葉尖結(jié)構(gòu)葉型壓力面應(yīng)力分布圖

表3 不同葉尖結(jié)構(gòu)葉型吸力面最大應(yīng)力

圖5 不同葉尖結(jié)構(gòu)葉型吸力面應(yīng)力分布圖

表4 不同葉尖結(jié)構(gòu)葉型壓力面最大應(yīng)力

表5 不同葉尖結(jié)構(gòu)葉型兩側(cè)最大應(yīng)力差值

葉片當量應(yīng)力分布圖可以看出，由于葉尖小翼結(jié)構(gòu)在葉片尖部增加了重量，三種葉尖小翼結(jié)構(gòu)葉片尖部的應(yīng)力分布與原始葉型相比均有所變化，原本僅在葉身出現(xiàn)的較高應(yīng)力區(qū)域沿葉尖小翼結(jié)構(gòu)延伸到了葉頂平面，而在葉根處的最大應(yīng)力值也均有所上升，但葉尖小翼結(jié)構(gòu)對葉身的應(yīng)力分布規(guī)律并沒有明顯影響。

葉尖小翼結(jié)構(gòu)同樣對葉片最大應(yīng)力的分布有較大影響，由于葉尖小翼改變了葉片的質(zhì)心，四種葉片結(jié)構(gòu)僅原始葉型最大當量應(yīng)力出現(xiàn)在葉片吸力面。與原始葉型相比，吸力面小翼的質(zhì)心偏移量較小，其對最大應(yīng)力的增量和兩側(cè)的應(yīng)力差值影響都相對較??；而壓力面小翼和雙側(cè)小翼則均較大影響，壓力面小翼由于質(zhì)心偏移量大，葉片兩側(cè)的最大應(yīng)力分布有很顯著的變化。

2.2 葉尖小翼結(jié)構(gòu)對葉片位移的影響

圖6和表6為不同葉尖結(jié)構(gòu)的綜合位移示意圖和最大綜合位移值。葉片的最大綜合位移均出現(xiàn)在葉片尖部尾緣位置，其中吸力面和壓力面小翼的最大位移相比原始葉型有不同程度的增加，而雙側(cè)小翼則與原始葉型變化不大。

圖6 不同葉尖結(jié)構(gòu)葉型總位移分布圖

表6 不同葉尖結(jié)構(gòu)葉型綜合位移

3 結(jié)論

本文對原始葉型、壓力面小翼葉型、吸力面小翼和雙側(cè)小翼葉型進行了應(yīng)力計算，分析了四種不同葉尖結(jié)構(gòu)葉片的最大當量應(yīng)力和最大綜合位移。通過對不同算例的計算結(jié)果進行分析，可以得出以下結(jié)論：

1)渦輪葉尖小翼對葉身的應(yīng)力分布規(guī)律無明顯影響，但是會使葉尖處的高應(yīng)力區(qū)沿葉尖小翼結(jié)構(gòu)延伸至葉頂平面；

2)由于葉尖小翼結(jié)構(gòu)改變了葉片的質(zhì)心，對葉片兩側(cè)的最大應(yīng)力差值有很大影響，其中壓力面小翼由于偏離較大，葉片兩側(cè)最大應(yīng)力相差30 MPa；

3)雙側(cè)葉尖小翼結(jié)構(gòu)對葉片的綜合位移影響較小，壓力面和吸力面小翼結(jié)構(gòu)葉片的綜合位移相比原始葉型分別增加0.029 mm、0.032 mm。

盡管葉尖小翼結(jié)構(gòu)作為一種有效抑制葉尖泄漏的結(jié)構(gòu)被廣泛研究，但其對葉片強度有顯著的影響。本文所論述的內(nèi)容關(guān)系到葉片的可靠性和使用壽命，對帶葉尖小翼的渦輪葉片設(shè)計有積極的意義。