透平末級(jí)葉片材料抗水蝕特性的數(shù)值研究

邸娟,王順森,蔣希航,顏曉江,周偉

(1.太原科技大學(xué)機(jī)械工程學(xué)院,030024,太原;2.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;3.紐約州立大學(xué)石溪分校機(jī)械工程系,11794,美國(guó)紐約;4.中南大學(xué)軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,410075,長(zhǎng)沙)

火電和核電凝汽式汽輪機(jī)的末級(jí)葉片長(zhǎng)期工作在濕蒸氣區(qū),在運(yùn)行過(guò)程中常常遭受嚴(yán)重的水滴沖蝕。沖蝕不僅會(huì)導(dǎo)致葉片的氣動(dòng)性能惡化,降低級(jí)效率,而且威脅汽輪機(jī)的安全運(yùn)行,是目前亟需解決的關(guān)鍵問(wèn)題。隨著汽輪機(jī)末級(jí)葉片發(fā)展得越來(lái)越長(zhǎng),如不采取有效措施,水蝕現(xiàn)象將變得更嚴(yán)重。

目前常用的防水蝕措施主要為2種:主動(dòng)式和被動(dòng)式防水蝕。主動(dòng)式防水蝕是指從汽輪機(jī)結(jié)構(gòu)設(shè)計(jì)方面減少有害水滴的產(chǎn)生。被動(dòng)式防水蝕是指采用更抗水蝕的葉片材料或針對(duì)葉片易水蝕的區(qū)域進(jìn)行強(qiáng)化。被動(dòng)式防護(hù)的核心問(wèn)題為高速固液撞擊引起的材料損傷。研究高速液固撞擊機(jī)理可為提高葉片水蝕防護(hù)能力提供關(guān)鍵性的意見(jiàn)或建議。關(guān)于汽輪機(jī)葉片水蝕及高速液固撞擊問(wèn)題,國(guó)內(nèi)外許多學(xué)者進(jìn)行了大量研究。

在數(shù)值模擬方面,國(guó)內(nèi)外許多研究者嘗試?yán)没谟邢拊臄?shù)值方法來(lái)討論材料的水蝕問(wèn)題。數(shù)值方法不僅能夠靈活便捷地針對(duì)不同因素(如材料類型、顆粒速度、尺寸、形狀和沖擊角等)的水蝕行為進(jìn)行模擬,而且可以大大節(jié)約試驗(yàn)周期和成本。Mabrouki等使用任意歐拉-拉格朗日耦合方法(ALE)對(duì)彈塑性材料鋁合金(A2024T3)進(jìn)行了數(shù)值模擬分析,模擬結(jié)果很好地展示了撞擊過(guò)程中射流和靶材的運(yùn)動(dòng)[11]。Haller等著重模擬了速度為500 m/s,直徑為200 μm的液滴撞擊剛體平面的過(guò)程,闡述了撞擊過(guò)程中液體的壓縮過(guò)程、激波的傳遞過(guò)程及側(cè)向射流的產(chǎn)生[12]。Ma將有限元法和光滑粒子流體動(dòng)力學(xué)(SPH)方法相結(jié)合,研究了高速水射流沖擊靶材的過(guò)程并描述了水蝕機(jī)理;通過(guò)與試驗(yàn)結(jié)果的比較,證明了該耦合方法的有效性,并比較了有限元耦合光滑粒子流體動(dòng)力學(xué)方法和ALE方法,指出在相同計(jì)算精度條件下,有限元耦合光滑粒子流體動(dòng)力學(xué)方法比ALE方法更節(jié)省計(jì)算資源[13]。Hsu等使用耦合拉格朗日-歐拉法(CEL),模擬了速度為570 m/s的水射流撞擊PMMA靶材的過(guò)程,發(fā)現(xiàn)數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果吻合,且數(shù)值模擬對(duì)試驗(yàn)難以得到的細(xì)節(jié)做了很好的補(bǔ)充[14]。汪勇等通過(guò)SPH耦合有限元算法(FEM)模擬計(jì)算了高速液滴撞擊PMMA的過(guò)程,指出撞擊后產(chǎn)生的高速側(cè)向射流是固體發(fā)生破壞的主要原因,同時(shí)對(duì)撞擊過(guò)程的激波動(dòng)態(tài)結(jié)構(gòu)進(jìn)行了分析,得到了撞擊過(guò)程中壓力分布的合理解釋[15]。

綜上所述,目前已有很多學(xué)者針對(duì)液固撞擊過(guò)程進(jìn)行了大量理論、試驗(yàn)和數(shù)值研究?？傮w而言,由于水蝕過(guò)程的復(fù)雜性,加之其在不同靶材和不同水蝕環(huán)境下的差異性,因此對(duì)材料水蝕規(guī)律和機(jī)理還缺乏深入的理解,如靶材表面形貌在水蝕過(guò)程中的演化特征,水蝕失效模式及不同表面強(qiáng)化工藝影響其抗水蝕性能的內(nèi)在深層原因,材料表面粗糙度、射流速度及角度對(duì)水蝕性能的影響等。本文基于FEM-SPH耦合算法,分析射流角度、射流速度及材料表面粗糙度對(duì)汽輪機(jī)典型末級(jí)葉片基材17-4PH抗水蝕性能的影響,有助于進(jìn)一步理解葉片材料水蝕規(guī)律,指導(dǎo)工藝參數(shù)優(yōu)化。

1 液固撞擊理論和數(shù)值方法

1.1 液固撞擊理論

Cook[16]假定被撞擊固體是剛體并且撞擊是一維的,在此情況下,根據(jù)能量守恒原理,液體在撞擊前的動(dòng)能轉(zhuǎn)化為液體的壓縮能,由此可推導(dǎo)出“水錘”的壓力公式為

P=ρ0C0V0

(1)

式中:P為液體撞擊固體所產(chǎn)生的撞擊壓力,即“水錘”壓力,Pa;ρ0為未受擾動(dòng)時(shí)液體的密度,kg/m3;C0是未受擾動(dòng)時(shí)液體的聲速,m/s;V0為射流速度,m/s。

由式(1)可知,撞擊產(chǎn)生的壓力在發(fā)生撞擊之后會(huì)急速下降,這是因?yàn)樵谶@之后從撞擊區(qū)域向液體和靶材中釋放了壓力波。如果撞擊持續(xù)進(jìn)行,由伯努利方程可推得,撞擊壓力將會(huì)達(dá)到一個(gè)穩(wěn)定值。平穩(wěn)階段的撞擊壓力為

(2)

(3)

(4)

其中k為考慮液體性質(zhì)的常數(shù),在液體是水的情況下取2。

考慮到在數(shù)值計(jì)算中通常采用三次方關(guān)系的激波關(guān)系式[18-19],它比Heymann線性擬合的關(guān)系式(3)準(zhǔn)確度高,因此之后在討論水錘壓力理論值時(shí)都用激波三次方關(guān)系式計(jì)算得出的數(shù)值,即

(5)

1.2 計(jì)算模型

圖1是水柱射流沖擊靶材的計(jì)算模型,取射流直徑為0.15 mm,射流高度為1 mm。射流與靶材之間取間隔距離為0.1 mm,以防止在計(jì)算前模型有初始穿透。設(shè)置靶材長(zhǎng)、寬均為1 mm,高度為0.3 mm。將固體靶材底面完全約束,并在其底面和4個(gè)側(cè)面上施加無(wú)反射邊界條件以模擬無(wú)窮大平板的射流撞擊過(guò)程。

為與高速水柱射流沖擊試驗(yàn)條件匹配,射流參數(shù)選取和試驗(yàn)條件一致,如表1所示,具體試驗(yàn)參數(shù)及試驗(yàn)步驟見(jiàn)文獻(xiàn)[20-21]。采用17-4PH作為靶材材料,射流速度取為650 m/s,靶材在試驗(yàn)過(guò)程中繞軸旋轉(zhuǎn),輪盤轉(zhuǎn)速n0為1 500 r/min,根據(jù)式(6)計(jì)算其在碰撞點(diǎn)的旋轉(zhuǎn)線速度為

Vr=wr=2πn0r

(6)

式中:r是與碰撞點(diǎn)射流相對(duì)應(yīng)的圓盤半徑。本試驗(yàn)中r=144 mm,求得旋轉(zhuǎn)線速度Vr≈22.62 m/s。

表1 模擬采用射流參數(shù)

在數(shù)值模擬中,將靶材的旋轉(zhuǎn)速度等效為22.62 m/s的平動(dòng)約束,計(jì)算時(shí)間為2.3 μs。通過(guò)在靶材表面開(kāi)設(shè)6個(gè)凹槽結(jié)構(gòu)的方式模擬靶面初始粗糙度,凹槽深度分別為2 μm和10 μm,見(jiàn)圖2。本文共模擬計(jì)算了4個(gè)角度射流撞擊靶材的過(guò)程。

1.3 物性參數(shù)及邊界條件

考慮到汽輪機(jī)末級(jí)的實(shí)際工作環(huán)境,將試驗(yàn)箱體內(nèi)氣壓設(shè)置為0.1～0.15個(gè)大氣壓,故可認(rèn)為射流從噴嘴中噴出后空氣對(duì)其影響很小,在模擬過(guò)程中忽略空氣的存在,僅考慮水柱射流和固體靶材的相互作用。

(a)2 μm凹槽90°模型圖

本文選用馬氏體基材17-4PH作為靶材,對(duì)其進(jìn)行高速水柱射流水蝕模擬。靶材的本構(gòu)模型選用隨動(dòng)塑性材料模型,具體材料參數(shù)如表2所示。

表2 靶材的材料參數(shù)

由于水滴在撞擊時(shí)壓力較高,受到較大壓力,因此需考慮其壓縮性。采用Mie-Gruneisen狀態(tài)方程[19]描述水在撞擊過(guò)程中的狀態(tài)。試驗(yàn)中射流由蒸餾水箱中經(jīng)過(guò)高速射流發(fā)生器產(chǎn)生,故考慮其溫度為室溫。水射流采用常溫常壓下的物性參數(shù)。水的初始物性參數(shù)按標(biāo)準(zhǔn)大氣壓下20 ℃時(shí)選取,C為Us-Up曲線(沖擊波速度Us和質(zhì)點(diǎn)速度Up曲線)的截距,即當(dāng)?shù)芈曀?γ0為Gruneisen常數(shù);μ為水的動(dòng)力黏性系數(shù),取1.002×10-3Pa·s;a是γ0的一階修正系數(shù);S1、S2、S3為Us-Up曲線的斜率系數(shù),也稱為水狀態(tài)方程系數(shù)。水的物性參數(shù)和Mie-Gruneisen狀態(tài)方程參數(shù)見(jiàn)表3[19]。

表3 水的物性參數(shù)和Mie-Gruneisen狀態(tài)方程參數(shù)[19]

1.4 模型驗(yàn)證

水柱射流沖擊剛性靶材的初始沖擊壓力P為1.954 GPa,如圖3所示。與理論水錘壓力相比,計(jì)算誤差為5.4%,在誤差允許范圍內(nèi)。

圖3 水射流垂直沖擊剛性靶材的初始沖擊壓力云圖

將靶材材料設(shè)置為剛體,在保證SPH粒子均勻分布的情況下,逐漸增加粒子密度,直到粒子撞擊靶材過(guò)程中撞擊壓力網(wǎng)格無(wú)關(guān)。最終選擇的SPH粒子密度為每層209個(gè),共107層。確認(rèn)SPH粒子密度后,將靶材網(wǎng)格密度從長(zhǎng)×寬×高為150×150×15增加到320×320×32,計(jì)算所得初始撞擊壓力隨網(wǎng)格密度的變化如圖4所示。從圖中可以看出,網(wǎng)格劃分大于240×240×24以后初始撞擊壓力趨于穩(wěn)定,300×300×30時(shí)初始撞擊壓力為1.954 GPa,320×320×32時(shí)初始撞擊壓力為1.96 GPa,相差0.31%,可認(rèn)為靶材網(wǎng)格無(wú)關(guān),取網(wǎng)格尺寸為300×300×30。有限元模型的網(wǎng)格拓?fù)潢P(guān)系見(jiàn)圖5,采用三維六面體8節(jié)點(diǎn)實(shí)體顯式單元solid164對(duì)靶材進(jìn)行劃分網(wǎng)格。

圖4 初始撞擊壓力隨網(wǎng)格密度變化曲線

圖5 有限元模型的網(wǎng)格拓?fù)潢P(guān)系圖(放大倍數(shù)約為5)

2 數(shù)值計(jì)算結(jié)果與分析

圖6以30°水柱射流撞擊10 μm凹槽靶材為例,選取6個(gè)典型時(shí)刻展示水柱射流過(guò)程。

(a)0 μs

從圖6可以看出,水柱與靶材發(fā)生初始接觸時(shí)并不會(huì)產(chǎn)生碰撞,需要射流繼續(xù)向前運(yùn)動(dòng)與靶材相互作用才會(huì)開(kāi)始碰撞過(guò)程。水柱與靶材碰撞瞬間,沒(méi)有產(chǎn)生側(cè)向射流,水柱前端為壓縮狀態(tài)。隨后靶材表面應(yīng)力呈環(huán)狀由撞擊中心區(qū)域向外部區(qū)域傳播,水柱射流在發(fā)生初始撞擊后出現(xiàn)側(cè)向射流,且撞擊后由撞擊表面產(chǎn)生壓力波向射流內(nèi)部傳播,應(yīng)力最大值出現(xiàn)在射流水柱外圈與靶材接觸的環(huán)狀區(qū)域。隨后不再有新的高應(yīng)力區(qū)出現(xiàn),之前產(chǎn)生的應(yīng)力波繼續(xù)向靶材內(nèi)部傳播。隨著射流繼續(xù)撞擊靶材,靶材表面逐步被破壞,形成沖坑。

圖7為不同粗糙度(光滑靶材、2 μm及10 μm凹槽)靶材內(nèi)部最大有效應(yīng)力隨沖擊過(guò)程的變化,可以看到,2 μm和10 μm凹槽靶材內(nèi)部最大有效應(yīng)力最終均穩(wěn)定在730 MPa附近,而光滑靶材內(nèi)部最大有效應(yīng)力最終穩(wěn)定在280 MPa左右,比帶凹槽結(jié)構(gòu)靶材的最終穩(wěn)定應(yīng)力小得多。

圖7 不同粗糙度光滑靶材內(nèi)部最大有效應(yīng)力隨沖擊過(guò)程的變化

(a)光滑表面30°沖擊

2.1 表面粗糙度和射流角度對(duì)水蝕的影響

(a)2 μm凹槽30°沖擊

(a)10 μm凹槽30°沖擊

圖8～圖10分別為光滑表面、2 μm及10 μm凹槽表面在不同射流角度(30°、60°、80°、90°)水柱射流沖擊下,靶材的微觀形貌和有效應(yīng)力(即von Mises應(yīng)力)云圖?？梢钥闯?對(duì)于每一種表面結(jié)構(gòu),90°沖擊時(shí)形成的沖坑最嚴(yán)重,質(zhì)量損失也最大,30°射流沖擊下試樣沖坑最淺,尤其是光滑表面,基材僅沖出非常淺的沖坑;在其他條件相同的情況下,靶材表面粗糙度越大,水蝕越嚴(yán)重;10 μm凹槽試樣比2 μm凹槽試樣更容易形成沖坑,沖蝕損傷也更嚴(yán)重。即材料表面初始粗糙度越大,水蝕破壞越嚴(yán)重,抗水蝕性能越差。

圖11給出了不同粗糙度(光滑表面、2 μm凹槽、10 μm凹槽)馬氏體基材的質(zhì)量損失隨水柱射流沖擊角的變化?？梢钥闯?隨著基材表面粗糙度的增加,水蝕累積質(zhì)量損失增大,水蝕破壞加劇,且最大質(zhì)量損失發(fā)生在沖擊角為90°;10 μm凹槽時(shí)基材的質(zhì)量損失平均為2 μm凹槽試樣的1.6倍,2 μm凹槽試樣的平均質(zhì)量損失約為光滑表面試樣的1.4倍。90°水柱射流沖擊試驗(yàn)中,表面未拋光的線切割基材質(zhì)量損失平均為其拋光樣的1.86倍。結(jié)合數(shù)值模擬結(jié)果,充分說(shuō)明表面粗糙度對(duì)材料水蝕性能影響顯著,隨著靶面粗糙度的增加,水蝕損傷加劇。Hancox也報(bào)道了類似的研究[9]:以18/8不銹鋼為研究對(duì)象,射流直徑為1.3 mm,試樣的表面粗糙度從37 μm變?yōu)? μm,測(cè)量了表面凹陷形成和材料去除所需的沖擊次數(shù),發(fā)現(xiàn)當(dāng)表面粗糙度從12 μm變?yōu)? μm時(shí),水蝕破壞減少了40%。

為進(jìn)一步分析材料表面粗糙度對(duì)水蝕性能影響的深層機(jī)理,對(duì)不同粗糙度下靶材最大沖擊壓力進(jìn)行了統(tǒng)計(jì)。圖12展示了不同粗糙度的馬氏體基材17-4PH的最大沖擊壓力隨水柱射流沖擊角的變化?？梢钥闯?材料的最大沖擊壓力隨水柱射流沖擊角的增加而增大,當(dāng)射流角度為90°時(shí),撞擊壓力最大;隨著表面粗糙度的增加,各射流角度下,靶面最大沖擊壓力升高,靶材質(zhì)量流失加劇,因此抗水蝕性能較差。基材最大撞擊壓力與靶面粗糙度的關(guān)系,充分說(shuō)明提高材料表面光潔度能促進(jìn)水滴壓力波的自由擴(kuò)展,降低材料最大撞擊壓力可以提高材料的抗水蝕性能。

圖12 不同粗糙度馬氏體基材的最大沖擊壓力隨水柱射流角度的變化

圖13給出了30°射流角度的水柱沖擊光滑表面,馬氏體基材最大壓力隨沖擊時(shí)間的分布?？梢钥闯?靶材最大沖擊壓力在水柱與靶面的接觸瞬間迅速增加,隨著射流沖擊過(guò)程的繼續(xù)進(jìn)行,最大壓力在250～700 MPa范圍內(nèi)波動(dòng),在1.9 μs后,壓力最終穩(wěn)定在200 MPa附近。

圖13 光滑表面馬氏體基材在30°水柱射流沖擊下最大壓力隨沖擊時(shí)間的變化

2.2 射流速度對(duì)水蝕的影響

以2 μm凹槽粗糙度為例,選擇4種水柱射流沖擊速度325、418、563、650 m/s分析水柱射流速度對(duì)水蝕的影響。

水蝕質(zhì)量損失Δm和射流速度V0具有以下指數(shù)關(guān)系

(7)

式中:n為速度指數(shù)。使用對(duì)數(shù)坐標(biāo)圖擬合得到馬氏體基材17-4PH速度指數(shù)n≈3.83,見(jiàn)圖14。從圖14中可以看出,射流速度越高,水蝕累積質(zhì)量損失越大,材料水蝕損害越嚴(yán)重。Ahmad等研究結(jié)果顯示,對(duì)于高屈服強(qiáng)度鋼合金,n=3.3～4.5;對(duì)于Ti6Al4V,n≈5[4]。本文研究的馬氏體基材17-4PH屬于高屈服強(qiáng)度鋼合金,擬合所得速度指數(shù)n≈3.83,與文獻(xiàn)[4]報(bào)道的結(jié)果吻合。

圖14 馬氏體基材累積質(zhì)量損失隨水柱射流速度的變化

3 結(jié) 論

本文采用SPH-FEM耦合算法建立了水柱射流撞擊模型,模擬了高速射流與靶材的撞擊過(guò)程。通過(guò)數(shù)值方法分析得到了固體靶材在不同角度、不同速度射流撞擊下的應(yīng)力變化及表面粗糙度對(duì)射流撞擊過(guò)程的影響,主要結(jié)論如下。

(1)相同靶材表面結(jié)構(gòu)下,射流角度越大,靶材水蝕越嚴(yán)重。

(2)靶材表面粗糙度對(duì)水蝕影響較大。10 μm凹槽時(shí)材料的質(zhì)量損失平均為2 μm凹槽試樣的1.6倍,2 μm凹槽試樣的平均質(zhì)量損失約為光滑表面試樣的1.4倍。

(3)提高靶材表面光潔度,能促進(jìn)水滴壓力波的自由擴(kuò)展,顯著降低撞擊壓力,從而提高材料的抗水蝕性能。

(4)沖擊速度越大,材料水蝕失重越多。靶材的水蝕累積質(zhì)量損失與其射流速度成指數(shù)關(guān)系,擬合得到速度指數(shù)n≈3.83。