摘要：為了提高汽輪機葉片的可靠性預(yù)防電廠中汽輪機葉片疲勞斷裂的發(fā)生首先提出了葉片應(yīng)力和動態(tài)特性的數(shù)值分析方法然后基于動態(tài)應(yīng)力分析結(jié)果考慮了影響汽輪機葉片失效的各種因素引入Neuber 準(zhǔn)則和雨流計數(shù)法建立了汽輪機葉片的疲勞壽命預(yù)測模型并分析了一個5片成組的851mm 葉片疲勞壽命狀況．結(jié)果表明：提出的汽輪機葉片動態(tài)應(yīng)力三維有限元分析模型具有良好的工程精度葉片疲勞壽命預(yù)測可以定量地分析影響葉片安全運行的各種因素和葉片的疲勞壽命為現(xiàn)實中葉片的設(shè)計制造及汽輪機運行提供了參考．

關(guān)鍵詞：汽輪機；葉片振動；動態(tài)應(yīng)力；高周疲勞；疲勞壽命；數(shù)值分析

汽輪機葉片的運行環(huán)境十分復(fù)雜葉片事故時有發(fā)生葉片的檢修費用占汽輪機檢修費用的505％［1］其中大部分是由于葉片的疲勞斷裂所致［2］．目前汽輪機葉片設(shè)計和安全性考核主要采用的是建立在材料力學(xué)基礎(chǔ)上的許用應(yīng)力為基準(zhǔn)的靜強度準(zhǔn)則和以安全倍率為標(biāo)準(zhǔn)的動強度準(zhǔn)則．這些準(zhǔn)則對于汽輪機葉片的設(shè)計發(fā)揮了重要的作用但許用應(yīng)力方法無法精細(xì)考慮葉片局部的應(yīng)力和由于交變載荷產(chǎn)生的疲勞失效而安全倍率方法是基于葉片事故統(tǒng)計數(shù)據(jù)的對于影響葉片安全性的許多因素如激振力、結(jié)構(gòu)特點、阻尼等無法定量估算其方法是比較粗糙的．隨著社會發(fā)展對于汽輪發(fā)電機組制造和運行的進(jìn)一步要求以上準(zhǔn)則已經(jīng)不能完全滿足需要．本文以葉片的三維動態(tài)應(yīng)力分析為基礎(chǔ)考慮了多種對葉片疲勞壽命有影響的因素建立了汽輪機葉片的疲勞壽命分析模型最后分析了5片成組的851mm 葉片的疲勞壽命．

1 葉片三維有限元動態(tài)應(yīng)力分析

動態(tài)應(yīng)力分析是葉片疲勞壽命分析的基礎(chǔ)有許多的學(xué)者如朱寶田［3］、鄭潤生［4］ 等人發(fā)展了各種葉片的動態(tài)應(yīng)力分析方法并且取得了很大的成果．但是他們或者是對葉片的一個部分進(jìn)行分析或者人為地忽略了葉片結(jié)構(gòu)的某些因素不能反映葉片的整體應(yīng)力狀況分析模型距離真實葉片尚存在一些差距．因此有必要發(fā)展出一種適合于葉片動態(tài)應(yīng)力分析的真實三維有限元分析模型．

1．1 葉片計算的三維有限單元

線性八節(jié)點協(xié)調(diào)單元在承受彎曲載荷時由于寄生剪切效應(yīng)計算變形比實際結(jié)構(gòu)應(yīng)該產(chǎn)生的變形要小有時可能會使計算結(jié)果發(fā)生嚴(yán)重失真因此對于汽輪機葉片尤其是長葉片這種單元模型是不太合適的．三維八節(jié)點非協(xié)調(diào)單元［5］ 是在線性八節(jié)點協(xié)調(diào)單元的基礎(chǔ)上在單元內(nèi)部增加了9個自由度大大改善了單元的抗彎性能保證了使用此單元對葉片進(jìn)行應(yīng)力分析結(jié)果的正確性．本文采用非協(xié)調(diào)單元進(jìn)行分析獲得了良好的計算精度．

1．2 葉片動態(tài)應(yīng)力場的求解

轉(zhuǎn)動葉片的運動微分方程［6］為

（1）

式中：［ M］是質(zhì)量矩陣；［ C］是阻尼矩陣；［ MG ］是陀螺矩陣；［ Kc ］是離心力剛度矩陣；［ K0］是線性剛度矩陣；VeNL是葉片單元考慮二次項后的非線性應(yīng)變能；｛Qc｝是葉片的離心力載荷向量；｛Qp｝是葉片表面上的氣流力向量；｛F（ t）｝是作用在葉片上的非穩(wěn)態(tài)力；｛δ¨｝、｛δ· ｝、｛δ｝是加速度、速度和位移向量．忽略葉片幾何大變形的影響上式中葉片的非線性應(yīng)變能可以寫為

（2）

式中：｛εl｝是非線性應(yīng)變；［ Kσ］是初應(yīng)力剛度矩陣．對于汽輪機葉片的靜力計算由式（1）可以得到以下方程

（3）

 又由靜力計算方程可以得到轉(zhuǎn)動葉片的動頻計算方程為

（4）

式中：［Φ］是特征向量；［Λ］是特征值即葉片的動頻．

在分析得到葉片的動頻和其特征向量以后就可以使用振型迭加法求得葉片的穩(wěn)態(tài)響應(yīng)．通過傅里葉分解處理作用于汽輪機葉片上的激振力則汽輪機葉片的穩(wěn)態(tài)響應(yīng)為

（5）

其中

（6）

式中：i 代表某一階模態(tài)；K 代表某階激振力；ζi 為模態(tài)阻尼比．

由葉片的響應(yīng)值很容易得到葉片的動態(tài)應(yīng)力．各個位移響應(yīng)就是某個時刻葉片所處的位移場根據(jù)位移場可以由下式

（7）

求得各時刻葉片有限元模型各個節(jié)點處的應(yīng)力而葉片各個節(jié)點應(yīng)力隨時間的變化就是結(jié)構(gòu)的動應(yīng)力．

2 葉片疲勞壽命分析方法

從本質(zhì)上來說疲勞斷裂是由于承受動態(tài)載荷所產(chǎn)生的即使在動應(yīng)力數(shù)值遠(yuǎn)遠(yuǎn)小于材料的靜強度的情況下仍然可能造成疲勞損壞［7］．一般將零件或部件的疲勞壽命劃分為疲勞裂紋初始壽命（簡稱為始裂壽命）和疲勞裂紋擴(kuò)展壽命．汽輪機葉片的大部分壽命都是消耗在疲勞裂紋起始上面所以研究中一般都把裂紋的起始壽命定義為葉片的疲勞壽命［6］在疲勞壽命分析中使用的“局部應(yīng)力應(yīng)變法”所研究的就是葉片的致裂壽命．在一般的情況下雖然名義應(yīng)力處于彈性范圍但從局部微觀角度看塑性變形仍然存在而使用有限元方法要分析到很細(xì)致的局部微觀應(yīng)力其人力和物力的消耗是極其巨大的所以本文采用具有良好的工程精度、使用起來十分方便的Neuber 準(zhǔn)則［8］ 來分析葉片的局部應(yīng)力應(yīng)變分析中用到的名義應(yīng)力就是由本文前面三維有限元模型計算得出的．在分析完局部應(yīng)力應(yīng)變后本文采用雨流計數(shù)法進(jìn)行計數(shù)．這種方法從計數(shù)原則上看與材料的應(yīng)力應(yīng)變遲滯回線相一致能夠反映隨機載荷的全過程．

局部應(yīng)力應(yīng)變法是在低周疲勞壽命估算方法的基礎(chǔ)上發(fā)展起來的在對汽輪機葉片進(jìn)行高周疲勞壽命分析時需要進(jìn)行一些必要的改進(jìn)．本文中的模型考慮了尺寸效應(yīng)（使用尺寸系數(shù)εd 表示）、葉片表面加工和表面處理（使用表面狀態(tài)系數(shù)β來表示β＝βsβcβpβs 是表面加工系數(shù)、βc 是腐蝕系數(shù)、βp 是表面強化系數(shù)）、微動磨損等方面的影響（使用微動疲勞系數(shù)αf 表示）．以上系數(shù)的一些具體數(shù)值可以參考文獻(xiàn)［9］選取．考慮到以上疲勞強度和平均應(yīng)力的影響對Manson－Coffin 曲線的彈性線進(jìn)行修正［10］,取

（8）

將彈性線中的b 以b′代替引入Morrow 對Manson－Coffin 曲線的修正可以得到修正的應(yīng)變－壽命曲線用下式表示

（9）

本文在進(jìn)行疲勞累積損傷分析時使用以下修正線性累積損傷公式進(jìn)行計算

（10）

考慮到葉片的實際情況根據(jù)文獻(xiàn)［9］的建議取α為068．

3 5片成組的851mm 葉片疲勞壽命分析

作為實例本文對某國產(chǎn)300MW 汽輪機組的末級5片成組的851mm 葉片進(jìn)行了分析．5片葉片通過葉頂Z 型圍帶和相對葉高為055處的Z 型拉筋焊接成組．葉片組的有限元計算網(wǎng)格見圖1．中間葉片的靜態(tài)應(yīng)力場分布如圖2．本文分析得出葉片的前3階動頻見表1同時給出了實驗值進(jìn)行對比．因為動頻的計算與葉片的靜態(tài)初應(yīng)力場及單元模型的選取是密切相關(guān)的從表1可看出動頻良好的計算精度這也從另一方面說明了本文所提供的分析葉片靜態(tài)應(yīng)力場模型的精確性．葉片的第2和第3階振型如圖3所示．根據(jù)表1動頻分析結(jié)果本文對葉片組進(jìn)行了三維動態(tài)應(yīng)力分析獲得了葉片型線、拉金和圍帶上危險點的等效應(yīng)力隨時間的變化．圖4是型線部分危險點的應(yīng)力狀態(tài)隨時間變化的曲線．獲得動應(yīng)力狀態(tài)后就可以進(jìn)行葉片的疲勞壽命分析了．

圖1 葉片組網(wǎng)格圖 

圖2 中間葉片靜態(tài)應(yīng)力分布

表1 5片成組的851mm 葉片的動頻值

（a）第2階       （b）第3階

圖3 851mm 葉片組的第2階和第3階振型

圖4 葉型上危險點的等效應(yīng)力時域圖

  根據(jù)851mm 葉片的材料特性數(shù)據(jù)和所處的工作環(huán)境對Manson－Coffin 曲線進(jìn)行修正在其中取尺寸系數(shù)εd ＝05表面加工系數(shù)βs ＝10表面強化系數(shù)βp ＝10因葉片位于濕蒸汽區(qū)取腐蝕疲勞系數(shù)βc＝06得到葉片的應(yīng)變－壽命曲線如下

（11）  

在汽輪機的運行過程中如果由于水質(zhì)不良將使葉片的腐蝕疲勞系數(shù)進(jìn)一步下降因此又選取腐蝕疲勞系數(shù)βc ＝0 5進(jìn)行了分析分析結(jié)果見表2．

表2 851mm 成組葉片壽命分析結(jié)果

從表2的分析結(jié)果可以看出對于851mm 的成組葉片來說型線部分壽命最短正常運行狀況時為691 1a可見正常使用的851mm 成組葉片是一個使用壽命比較長的葉片．

以上的分析涉及到了影響葉片壽命的各個因素并且通過具體的壽命計算分析了它們對葉片壽命的影響從中可以看出葉片材料的性能、葉片表面加工及熱處理工藝、葉片的運行工況及環(huán)境是影響葉片壽命的重要因素．

4 結(jié) 論

（1）本文所提出的葉片動態(tài)應(yīng)力三維有限元分析模型具有良好的工程精度能夠得到令人滿意的結(jié)果．

（2）本文建立的汽輪機葉片疲勞壽命預(yù)測模型考慮了汽輪機葉片制造和運行中的多種因素（如腐蝕環(huán)境、表面加工狀況等）對于葉片壽命的影響葉片壽命的分析結(jié)果不但可以得到一些定性方面的概念還具有一定程度的量化指導(dǎo)意義．

（3）本文基于葉片三維有限元動態(tài)應(yīng)力結(jié)果和金屬材料疲勞理論所提出的汽輪機葉片疲勞壽命預(yù)測方法不僅為汽輪機葉片設(shè)計提供了一個從壽命消耗方面評價葉片安全性的新評估方法并且為在汽輪機運行方式日益復(fù)雜的情況下進(jìn)行葉片壽命管理提供了新思路及實用工具．

參考文獻(xiàn)：

［1］　丁有宇周宏利徐　鑄等．汽輪機強度計算［M］．北京：水利電力出版社1985．220～226．

［2］　吳厚鈺．透平零件結(jié)構(gòu)和強度計算［M］．北京：機械工業(yè)出版社1982．56～60．

［3］　朱寶田．汽輪機葉片振動、激振力和動應(yīng)力特性及其優(yōu)化研究［D］．西安：西安交通大學(xué)能源與動力工程學(xué)院1993．

［4］　鄭潤生吳厚鈺．采用扭桿單元分析汽輪機長葉片振動的有限元法［J］．西安交通大學(xué)學(xué)報199125（3）：99～110．

［5］　謝永彗安　寧王樂天等．汽輪機葉片三維有限元應(yīng)力分析系統(tǒng)［J］．西安交通大學(xué)學(xué)報199730（8）：56～61．

［6］　謝永慧．汽輪機葉片疲勞失效、壽命預(yù)測及設(shè)計分析系統(tǒng)的研究［D］．西安：西安交通大學(xué)能源與動力工程學(xué)院1997

［7］ 　Naeem MSingh RProbert D．Implication of enginedeterioration for a high-pressure turbine blade′s low-cyclefatigue （LCF） life-consumption ［J］．Interna-tional Journalof Fatigue199921：831～847．

［8］　Fuchs H OStephens R I．Metal fatigue in engineering［M］．New York：John Wiley ＆ Sons Inc1980．141～144．

［9］　趙少汴．抗疲勞設(shè)計［ M ］．北京：機械工業(yè)出版社1994．114～120．

［10］ Newport AGlinka G．Effect of notch-strain calculationmethod on fatigue-crack-initiation life predictions ［J ］．Experimental Mechanics199030（2）：208～216．