王秀瑾

(上海發(fā)電設(shè)備成套設(shè)計研究院，上海 200240)

1 問題的提出

汽輪發(fā)電機組是核電常規(guī)島的核心設(shè)備，其轉(zhuǎn)子支撐系統(tǒng)的動特性對整個機組的安全、穩(wěn)定運行有著重要影響。而大型汽輪發(fā)電機組的支撐軸承是轉(zhuǎn)子支撐系統(tǒng)的重要組成部分，其軸承動剛度(特別是較高轉(zhuǎn)速時的動剛度)對轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速影響較大。百萬千瓦級核電汽輪發(fā)電機轉(zhuǎn)子多在超過一階臨界轉(zhuǎn)速的區(qū)域內(nèi)使用，因此，汽輪發(fā)電機軸承座剛度的好壞直接影響汽輪發(fā)電機轉(zhuǎn)子－軸承－基礎(chǔ)這一系統(tǒng)能否安全運行［1－3］。隨著設(shè)備安裝和運行情況的變化，軸承座剛性也會有所變化，從而導致汽輪發(fā)電機軸系的動力特性變化。為使轉(zhuǎn)子的工作轉(zhuǎn)速避開臨界轉(zhuǎn)速，往往會采用改變支承剛性的方法。核電汽輪發(fā)電機所采用的滑動軸承支撐主要分為軸承箱支撐和基座端蓋支撐，也就是軸承箱式落地軸承座和端蓋式軸承座。百萬千瓦級核電半速機組汽輪發(fā)電機采用了端蓋式軸承座，由于端蓋式軸承座承擔著汽輪發(fā)電機轉(zhuǎn)子的靜、動載荷以及由于傳遞扭矩而帶來的反作用力，可獲得其端蓋式軸承座的剛度特性，研究其轉(zhuǎn)子動力特性對軸承剛度變化的敏感性，對采用端蓋式軸承座的同類核電汽輪發(fā)電機設(shè)計尤其是動特性設(shè)計是很重要的。

目前，在國內(nèi)百萬千瓦級核電半速汽輪發(fā)電機組軸系動特性設(shè)計計算時，軸承支撐剛度尤其是汽輪發(fā)電機轉(zhuǎn)子軸承支撐剛度通常采用已有的經(jīng)驗值。在工程實際中，可通過計算或用模型試驗方法研究不同支撐剛性對于轉(zhuǎn)子振動特性的影響，但軸承座模型設(shè)計、制造和試驗需要一定的周期，造成設(shè)計周期的延長，而模型試驗時各種原因造成的誤差也不能避免。百萬千瓦級核電機組動力特性研究系統(tǒng)通常由轉(zhuǎn)子、軸承、軸承座及框架基礎(chǔ)等多個子結(jié)構(gòu)組成，其中汽輪發(fā)電機軸承采用端蓋式軸承座。汽輪發(fā)電機轉(zhuǎn)子系統(tǒng)支承邊界模型如圖1所示。

圖1 轉(zhuǎn)子支承系統(tǒng)力學模型

圖1中:kij，cij分別為軸承油膜剛度和阻尼;kp為軸承座支撐剛度;kb為汽輪發(fā)電機基礎(chǔ)軸承座位置的支撐剛度。決定軸系動態(tài)特性的重要因素之一是轉(zhuǎn)子系統(tǒng)的支承特性，而支承特性又與軸承的油膜剛度阻尼特性、軸承座及基礎(chǔ)的支撐剛度有關(guān)［4］。軸承油膜動剛度阻尼可由專門的分析軟件計算獲得，軸承支座及基礎(chǔ)的剛度則通常由經(jīng)驗值獲得，但經(jīng)驗值的準確性有待于進一步研究確認。

本文采用國際通用的大型有限元軟件ANSYS分析了百萬千瓦級核電半速機組汽輪發(fā)電機端蓋式軸承座及其汽輪發(fā)電機基礎(chǔ)的剛度［5－7］，獲得了珍貴的軸承座剛度數(shù)值分析結(jié)果。綜合考慮軸承座計算值和汽輪發(fā)電機基礎(chǔ)支撐剛度計算值，參照汽輪發(fā)電機軸承支撐剛度經(jīng)驗數(shù)據(jù)，本文通過計算來研究軸承座剛度在一定范圍內(nèi)變化時對大型汽輪發(fā)電機轉(zhuǎn)子臨界轉(zhuǎn)速的影響。

2 汽輪發(fā)電機端蓋式軸承座剛度計算

百萬千瓦級核電半速機組汽輪發(fā)電機轉(zhuǎn)子軸承采用了端蓋式軸承座，本文以該軸承座為研究對象，運用大型有限元分析軟件ANSYS分析了該軸承座的靜剛度和動剛度。對汽輪發(fā)電機端蓋式軸承座進行合理簡化，建立軸承座有限元模型。百萬千瓦級核電汽輪發(fā)電機軸承座有限元模型如圖2所示。

圖2 百萬千瓦級核電汽輪發(fā)電機軸承座三維有限元模型

由于汽輪發(fā)電機軸承座為端蓋式，軸承座底部固定，限制了軸承座底板全部自由度?？紤]軸承座上蓋的重量，根據(jù)轉(zhuǎn)子載荷作用下的軸承座位移來計算軸承座支撐動剛度。圖3、圖4分別是汽輪發(fā)電機端蓋式軸承座有限元模型在上述邊界條件下的靜剛度和動剛度計算用位移分布云圖，軸承座剛度計算結(jié)果見表1。

表1 軸承座剛度計算結(jié)果 GN/m

該軸承座在25Hz時的動剛度和靜剛度差異很小，這是因為作用在軸承座的動載荷頻率為25Hz，遠小于軸承座固有頻率;由于受載荷作用的軸承座支撐位置位移不同，計算的剛度值是一個范圍。

3 汽輪發(fā)電機基礎(chǔ)支撐剛度計算

由于彈簧隔振基礎(chǔ)可以降低汽輪發(fā)電機組振動對環(huán)境的影響并隔離地震對汽輪發(fā)電機組的影響，以及第一階固有頻率遠低于核電半速機的工作頻率，彈簧基礎(chǔ)越來越多地應用于核電半速機組中。彈簧隔振器使汽輪發(fā)電機平臺和下面的立柱脫離了動力耦合，立柱剛度遠大于汽輪發(fā)電機平臺板橫梁剛度，可直接建立彈簧基礎(chǔ)汽輪發(fā)電機平臺板模型?？紤]彈簧的支撐作用，百萬千瓦級核電半速機組汽輪發(fā)電機轉(zhuǎn)子端蓋式軸承座所在彈簧基礎(chǔ)支撐剛度計算用位移云圖如圖5所示，負荷作用下汽輪發(fā)電機基礎(chǔ)上各點位移響應不同，計算可得到汽輪發(fā)電機基礎(chǔ)剛度值為2．69～4．26GN/m。

為便于比較計算，軸承座剛度取平均值9．62 GN/m，將汽輪發(fā)電機端蓋式軸承座的剛度值和所配彈簧基礎(chǔ)橫梁的剛度耦合，得到的軸承總支撐剛

圖5 端蓋式軸承座所在彈簧基礎(chǔ)剛度計算用位移云圖

度值(此處主要考慮了垂直方向)見表2，經(jīng)驗數(shù)據(jù)也一并列于表2。

表2 汽輪發(fā)電機轉(zhuǎn)子軸承座剛度分析結(jié)果 GN/m

4 軸承座剛度對軸系臨界轉(zhuǎn)速的影響

決定軸系臨界轉(zhuǎn)速的重要因素之一是轉(zhuǎn)子系統(tǒng)的支撐特性，而支撐特性除了與軸承的油膜剛度、阻尼特性有關(guān)外，還與軸承座及基礎(chǔ)的支撐動剛度有關(guān)［8］。以軸承支撐總剛度經(jīng)驗值計算的汽輪發(fā)電機臨界轉(zhuǎn)速為參考基準，基于計算得出的汽輪發(fā)電機端蓋式軸承座剛度值及在經(jīng)驗值一定范圍內(nèi)變化的汽輪發(fā)電機軸承座剛度值，分別計算分析了汽輪發(fā)電機軸系臨界轉(zhuǎn)速的相對值，計算結(jié)果見表3。

表3 軸承支撐剛度變化時汽輪發(fā)電機臨界轉(zhuǎn)速的相對變化量

從表3可以看出，有限元分析得到的軸承座剛度接近經(jīng)驗數(shù)據(jù)，其對軸系臨界轉(zhuǎn)速有一定影響;同時，汽輪發(fā)電機臨界轉(zhuǎn)速隨著端蓋軸承支撐剛度增減而相應增減。汽輪發(fā)電機轉(zhuǎn)子端蓋式軸承座剛度變化時，汽輪發(fā)電機轉(zhuǎn)子的臨界轉(zhuǎn)速有所變化，汽輪發(fā)電機轉(zhuǎn)子臨界轉(zhuǎn)速在軸承座剛度降低時的相對變化量較剛度升高時大，這為汽輪發(fā)電機轉(zhuǎn)子動特性設(shè)計提供了重要參考。從安全角度考慮，在滿足設(shè)計規(guī)范的前提下，應按較高的軸承支撐剛度來設(shè)計。

5 結(jié)論

(1)軸承支撐剛度對汽輪發(fā)電機動特性(特別是臨界轉(zhuǎn)速)有著重要影響，在汽輪發(fā)電機剛度設(shè)計之初，精確計算剛度有助于機組的安全運行。

(2)本文建立的汽輪發(fā)電機轉(zhuǎn)子端蓋式軸承座剛度分析有限元模型，是獲得端蓋式軸承座剛度的有益探索。計算獲得的端蓋式軸承座綜合剛度與經(jīng)驗數(shù)值較為接近，軸承支撐剛度精確計算中需要考慮汽輪發(fā)電機基礎(chǔ)彈性。

(3)國產(chǎn)百萬千瓦級核電半速機組汽輪發(fā)電機轉(zhuǎn)子端蓋式軸承座變剛度計算分析表明，端蓋式軸承座剛度變化對汽輪發(fā)電機轉(zhuǎn)子臨界轉(zhuǎn)速有一定影響，獲得準確的軸承支撐剛度十分必要。汽輪發(fā)電機轉(zhuǎn)子臨界轉(zhuǎn)速在軸承座剛度降低時的相對變化量較剛度升高時大，應按較大的軸承支撐剛度來設(shè)計。

(4)本文所采用的有限元方法在工程上有很大的實用價值，在滿足實際工程需要的前提下，對工程問題進行合理簡化，以便于分析，符合工程實踐所要求的合理性、準確性及方便性原則。百萬千瓦級核電機組汽輪發(fā)電機轉(zhuǎn)子端蓋式軸承座的剛度計算研究數(shù)據(jù)難得，但要完全反映同級核電半速機組汽輪發(fā)電機轉(zhuǎn)子用軸承座支撐邊界規(guī)律，還需同類型重載軸承座有限元模型的完善和實測研究數(shù)據(jù)的驗證。

［1］鐘一諤．轉(zhuǎn)子動力學［M］．北京:清華大學出版社，1987．

［2］角田英治，武新華．蒸汽透平低壓缸的動態(tài)分析——提高透平安全性［J］．汽輪機技術(shù)，1985(2):33 －38．

［3］Dikmen Emre，Hoogt van der Peter，Boer de André，et al．AFlexible Rotor on Flexible Supports:Modeling＆Experiments［C］//ASME 2009 InternationalMechanical Engineering Congress＆ Exposition IMECE 2009．Lake Buena Vista:ASME 2009 International Mechanical Engineering Congress＆ Exposition IMECE，2009．

［4］楊建剛，黃葆華，高亹．汽輪發(fā)電機組軸承座動力特性識別方法研究［J］．中國電機工程學報，2001，21(3):24－27．

［5］祝效華，余志祥．ANSYS高級工程有限元分析范例精選［M］．北京:電子工業(yè)出版社，2004．

［6］孫慶，王秀瑾，吳茜．超臨界600MW汽輪發(fā)電機組軸系及框架基礎(chǔ)系統(tǒng)動力特性分析研究［C］//中國動力工程學會透平專業(yè)委員會2010年學術(shù)研討會論文集．上海:中國動力工程學會平專業(yè)委員會，2010．

［7］杜連秀，張偉忠．汽輪機低壓缸剛性有限元分析［J］．汽輪機技術(shù)，2005(6):428－429．

［8］SUNQing，WEIQi，HUANGGongwen，etal．DynamicAnalysisofRotorBearingFoundationSystemfor1000MWUltra SupercriticalTurbineGeneratorSet［J］．ChallengesofPower EngineeringandEnvironment，2007(5):297 －302．