組合轉(zhuǎn)子解析建模與拉桿蠕變導(dǎo)致的退化研究

劉 強，王艾倫，羅 真

（1.中南大學(xué) 高性能復(fù)雜制造國家重點實驗室，長沙，410083；2.中南大學(xué) 機電工程學(xué)院，長沙 410083）

組合轉(zhuǎn)子解析建模與拉桿蠕變導(dǎo)致的退化研究

劉 強1,2，王艾倫1,2，羅 真2

（1.中南大學(xué) 高性能復(fù)雜制造國家重點實驗室，長沙，410083；2.中南大學(xué) 機電工程學(xué)院，長沙 410083）

以鐵木辛柯梁理論為基礎(chǔ)建立集中參數(shù)解析模型，并以結(jié)構(gòu)局部損傷將引起系統(tǒng)退化為切入點，研究拉桿蠕變松弛引起的組合轉(zhuǎn)子退化特性?？紤]拉桿螺栓的應(yīng)力松弛與接觸界面的影響，并對接觸界面等效處理，對組合轉(zhuǎn)子進行動力學(xué)分析，得到組合轉(zhuǎn)子固有頻率隨時間的變化規(guī)律；通過定義組合轉(zhuǎn)子頻率相對降低量作為固有頻率特性退化量參數(shù)，實現(xiàn)了組合轉(zhuǎn)子退化的定量評估，從理論上得出由拉桿蠕變松弛引起的組合轉(zhuǎn)子固有頻率特性退化量隨時間的函數(shù)曲線。

振動與波；拉桿松弛；組合轉(zhuǎn)子；退化；接觸界面

在設(shè)備運行過程中，由于高溫蠕變、螺紋松弛、裂紋等諸多方面的影響，會導(dǎo)致設(shè)備性能逐步退化，這不但使得設(shè)備的整機效率降低，更會使故障發(fā)生的概率升高[1]，從某種意義上可以認為退化是導(dǎo)致設(shè)備失效（或故障）的重要原因之一。燃氣輪機由于長期服役于高轉(zhuǎn)速、高溫等惡劣環(huán)境中，研究其退化特性、掌握退化的趨勢及程度，對保障燃氣輪機長期安全運行具有重要的意義[2]。為此，國內(nèi)外均有相關(guān)研究。Diakunchak指出，在長期高溫高壓的惡劣環(huán)境中運行，燃氣輪機的整機效率會逐漸降低[3]。其性能降低速率主要受運行時間、工況環(huán)境和啟停次數(shù)的影響。Rainer Kurz詳描述了燃氣輪機的主要退化機理，分析了多種缺陷對整機效率的影響規(guī)律[4]。Brook指出，燃氣輪機能量損失的逐年增大主要是由壓氣機部分的性能退化引起的[5]。Aker研究了燃氣輪機整機效能對壓氣機端性能局部退化的敏感性[6]。Brothertom提出燃氣輪機性能退化規(guī)律與常規(guī)復(fù)雜機械系統(tǒng)相似，符合典型故障率曲線，曲線的形狀呈兩頭高、中間低，具有明顯的階段性[7]。Li.Y.G利用通流量作為系統(tǒng)退化指標(biāo)，通過實驗數(shù)據(jù)得出了燃氣輪機系統(tǒng)退化曲線[8]。

綜合以上研究成果，針對蠕變松弛引起的組合轉(zhuǎn)子退化進行的研究中，蠕變松弛是高溫緊固件的主要失效機理。由于組合轉(zhuǎn)子長期在高溫下運行，拉桿在服役期經(jīng)受著蠕變損傷機制的威脅，損傷累積至一定程度，拉桿預(yù)緊力降低，轉(zhuǎn)子結(jié)構(gòu)完整性降低，難以維持正常運轉(zhuǎn)，所以拉桿應(yīng)力松弛是燃氣輪機安全運行必須考慮的問題。為研究拉桿蠕變引起的組合轉(zhuǎn)子退化，首先建立組合轉(zhuǎn)子動力學(xué)集中參數(shù)解析模型，然后定義轉(zhuǎn)子固有頻率作為退化指標(biāo)，得出組合轉(zhuǎn)子的固有頻率特性退化曲線，從而揭示拉桿松弛引起的組合轉(zhuǎn)子退化機理。

1 組合轉(zhuǎn)子解析模型的建立

組合轉(zhuǎn)子基本結(jié)構(gòu)為由周向均布的多根拉桿沿軸向緊固在一起的多級輪盤，見圖1，是較為典型的離散轉(zhuǎn)子結(jié)構(gòu)，且整體形狀較粗短，故必須以鐵木辛柯梁理論為基礎(chǔ)建立其彎曲振動集中參數(shù)模型，其物理參數(shù)見表1。

圖1 組合轉(zhuǎn)子結(jié)構(gòu)簡圖

表1 組合轉(zhuǎn)子物理參數(shù)

圖2 鐵木辛柯梁的隔離體劃分

以單個輪盤為基本單元，見圖2，運用隔離體法并增加對應(yīng)阻尼后分別導(dǎo)出各盤的牛頓力學(xué)方程

首盤：

其中k為鐵木辛柯形狀系數(shù)，取k=0.9。

軸段彎曲剛度為

拉桿當(dāng)量彎曲剛度為

mi——第i盤的質(zhì)量；

xi——第i盤的縱向位移；

Ji——第i盤繞質(zhì)心的轉(zhuǎn)動慣量；

Gpi——第i盤與i+1盤之間接觸界面的切變剛度；

Tpi——mi與mi+1之間的軸段彎曲剛度；

T—i—第i盤與第i+1盤之間接觸界面的彎曲剛度；

E——輪盤的彈性模量；

ν——輪盤泊松比；

α——輪盤接觸環(huán)面的直徑比

Ij——第j根拉桿的極慣性矩

Kj——第j根拉桿的拉伸剛度

寫成矩陣形式，有

其中

2 拉桿螺栓的松弛應(yīng)力計算

應(yīng)力松弛是指被聯(lián)接件的相對位置不發(fā)生改變，聯(lián)接件由于長時間的受拉，彈性變形量逐漸降低，導(dǎo)致應(yīng)力隨之減小。聯(lián)接件盡管要受到蠕變、低周疲勞等多種因素作用，但在穩(wěn)定載荷狀態(tài)下，應(yīng)力松弛是其最基本的損傷機制[9]。

在應(yīng)力松弛中有

其中ε0為聯(lián)接件的初始彈性應(yīng)變，E為彈性模量，經(jīng)歷蠕變過程后，彈性應(yīng)變?yōu)棣舉，蠕變應(yīng)變?yōu)棣與。

對式(1)求導(dǎo)得

由Norton公式可得

結(jié)合式(2)、式(3)積分可得

式中σ0——初始應(yīng)力；

C1、C2、C3——材料常數(shù)。

基于單軸蠕變理論，拉桿剩余應(yīng)力σ符合下式

式中Ft為時間為t時刻拉桿預(yù)緊力，Arod為拉桿橫截面積，Δl為拉桿預(yù)緊量，lrod為拉桿總長度。參照工程實際中重型燃氣輪機拉桿的預(yù)緊情況，將拉桿初始預(yù)緊量定為拉桿總長度的千分之一，則σ0=E/1 000。

燃氣輪機組合轉(zhuǎn)子拉桿材料采用GH4169，利用該材料蠕變實驗數(shù)據(jù)，確定Norton蠕變模型的參數(shù)值[10]，其中C1=2.2×10-9，C2=10.2，C3=50 825。由文獻[11]可查得E在873 K溫度下的值，由式（4）可得如圖3所示應(yīng)力松弛曲線。

圖3 在873 K溫度下拉桿應(yīng)力松弛曲線圖

可以看出，在應(yīng)力小于材料屈服極限情況下，拉桿的松弛在初始階段比較明顯，隨著時間的延長，軸向應(yīng)力逐漸降低，松弛速率也隨之降低。

3 組合轉(zhuǎn)子界面剛度的等效

界面接觸剛度主要由接觸面相互作用力與界面表面特性決定[12]。接觸剛度隨接觸面壓力的增大而增大；接觸界面粗糙度、形貌度和波浪度等都是影響界面接觸剛度大小的幾何因素[13]，其中粗糙度為研究各種幾何因素的基礎(chǔ)[14]。文中解析模型僅考慮粗糙度與預(yù)緊力的影響，并用虛擬材料法對接觸界面進行等效處理。

3.1 界面剛度的確定

將接觸界面假設(shè)為按照正態(tài)分布的微元體模型接觸，實現(xiàn)微觀到宏觀的剛度轉(zhuǎn)換，從而計算出各級輪盤間界面剛度[15]：

①計算接觸界面名義接觸面積S1，通過受力分析得到接觸面的壓強P；

② 根據(jù)文獻[16]建立實際面積S2、粗糙度0.4 mm、壓力P作用下的微元體模型，微元體法向接觸剛度為Kn；

③ 聯(lián)立獲得界面總壓力F作用下等效剛度K，見圖4。

圖4 組合轉(zhuǎn)子界面等效剛度曲線

3.2 界面剛度的等效

因接觸剛度本身不具有線性特征，而難以直接對組合轉(zhuǎn)子進行動力學(xué)分析[17]，因此必須對接觸剛度進行線性假設(shè)?？捎汕拔闹械玫降牡刃Ы佑|剛度K根據(jù)應(yīng)變能相等，將接觸剛度替換為彈性模量的實體層，則有

h——等效實體層厚度；

設(shè)實體層厚度為0.5 mm，組合轉(zhuǎn)子的工作溫度為873 K。

表2列出了組合轉(zhuǎn)子在不同時刻的固有頻率。為方便表征組合轉(zhuǎn)子的退化特性，定義頻率相對降低量η如下

表2 組合轉(zhuǎn)子固有頻率隨時間的變化關(guān)系

其中(f0-ft)為組合轉(zhuǎn)子頻率降低量；式中f0為初始時刻組合轉(zhuǎn)子固有頻率，ft為拉桿松弛過程中的t時刻組合轉(zhuǎn)子固有頻率。根據(jù)表2可知，隨著松弛時間的增加，組合轉(zhuǎn)子固有頻率發(fā)生漂移。組合轉(zhuǎn)子固有頻率相對降低量隨時間變化關(guān)系如圖5所示。

圖5 組合轉(zhuǎn)子頻率相對降低量隨時間的變化

4 組合轉(zhuǎn)子的退化軌跡的描述

工程上通常以設(shè)備的主要性能參數(shù)變化率作為系統(tǒng)的退化量，以判斷設(shè)備的退化情況[18]。選取組合轉(zhuǎn)子固有頻率的相對降低量（即η）作為組合轉(zhuǎn)子固有頻率特性退化量參數(shù)。

對不同溫度下組合轉(zhuǎn)子退化數(shù)據(jù)進行分析，得到組合轉(zhuǎn)子不同溫度下的退化軌跡見圖6。

圖6 組合轉(zhuǎn)子在不同溫度下的退化軌跡

組合轉(zhuǎn)子的退化量對溫度比較敏感，退化速率隨溫度升高而升高。所以，為提高組合轉(zhuǎn)子的性能可靠性，降低組合轉(zhuǎn)子的退化量，組合轉(zhuǎn)子運行過程中在保證其他必要條件基礎(chǔ)上應(yīng)盡量降低拉桿的溫度。

5 結(jié)語

(1)受初始預(yù)緊力的影響，拉桿松弛在初始階段比較明顯，退化速率較快；隨著時間的延長，蠕變效應(yīng)累積，拉桿軸向應(yīng)力逐漸減小，松弛速率逐漸降低，退化速率也隨之降低。

(2)隨著松弛時間的增加，組合轉(zhuǎn)子固有頻率發(fā)生漂移。在同一時間，組合轉(zhuǎn)子固有頻率降低量具有不同的分布，1階彎振、2階彎振、3階彎振頻率降低量依次增大；而組合轉(zhuǎn)子固有頻率相對降低量在不同時間具有相近的分布值，且頻率相對降低量隨時間呈遞增趨勢。

(3)以組合轉(zhuǎn)子頻率相對降低量作為組合轉(zhuǎn)子固有頻率特性退化量參數(shù)，實現(xiàn)組合轉(zhuǎn)子固有頻率特性退化量的定量表達，可以從理論上得出轉(zhuǎn)子因拉桿蠕變而退化的函數(shù)曲線。

[1]LEE J.Measurement of machine performance degradation using a neural network model[J].Computers in Industry,1996,30(3):193-209.

[2]蘇永雷，王艾倫，李雪鵬，等.拉桿松弛對燃氣輪機組合轉(zhuǎn)子性能退化的影響[J].工程設(shè)計學(xué)報，2014(3)：245-250.

[3]DIAKUNCHAK IS.Performancedegradation in industrial gas turbines[J].ASME J.of Engineering for Gas Turbines and Power,1992,114:161-168.

[4]BROOKS F J.GE gas turbine performance characteristics [R].2003,GE Power Systems GER-3567H.

[5]AKER G,SARAVANAMUTTOO H.Predicting gas turbine performance degradation due to compressor fouling using computer simulation techniques[J].Journal of Engineering for Gas Turbines and Power,1989,111 (2):343-350.

[6]RAINER KURZ,KLAUS BRUN.Degradation in gas turbine systems[J].Journal of Engineering for Gas Turbines and Power,2001,123(1):70-77.

[7]BROTHERTON T,JAHNS G,JACOBS J,et al.Prognosis of Fault in Gas Turbine Engines[C].IEEE Aerospace Conference Proceedings,2000(6):163-171.

[8]LI Y G,NILKITSARANONT P.Gas turbine performance prognostic for condition-based maintenance[J].Applied Energy,2009(86):573-584.

[9]郭進全，軒福貞，王正東，等.基于蠕變的高溫構(gòu)件應(yīng)力松弛損傷模型[J].核動力工程，2009，30(4)：9-12.

[10]李玲.渦軸發(fā)動機轉(zhuǎn)子連接螺栓的強度與蠕變壽命研究[D].南京：南京航空航天大學(xué)，2008.

[11]于慧臣，吳學(xué)仁.航空發(fā)動機設(shè)計用材料數(shù)據(jù)手冊[M].北京：航空工業(yè)出版社，2010：85-107.

[12]GOERKE D,WILLNER K.Normal contact of fractal surfaces-experimental and numerical investigations[J].Wear,2008,264(7):589-598.

[13]肖會芳，卲毅敏，徐金梧.粗糙界面法向接觸振動響應(yīng)與能量耗散特性研究[J].振動與沖擊，2014，33(4)：149-155.

[14]易均，劉恒，劉意，等.接觸界面對拉桿組合柔性轉(zhuǎn)子軸承系統(tǒng)的非線性動力特性影響[J].振動與沖擊，2012，31 (17)：165-170.

[15]李輝光，劉恒，虞烈.粗糙機械結(jié)合面的接觸剛度研究［J］.西安交通大學(xué)學(xué)報，2011，45(6)：69-74.

[16]蘇永雷.拉桿松弛與組合轉(zhuǎn)子性能退化的關(guān)聯(lián)性研究[D].長沙：中南大學(xué)，2014.

[17]艾延廷，翟學(xué)，王志，等.法向接觸剛度對裝配體振動模態(tài)影響的研究[J].振動與沖擊，2012，31(6)：171-174.

[18]鄧愛民，陳循，張春華，等.基于性能退化數(shù)據(jù)的可靠性評估[J].宇航學(xué)報，2006，27(3)；546-552.

Study onAnalytical Modeling and Performance Degradation Induced by Rod Relaxation of Combined Rotors

LIU Qiang1,WANG Ai-lun1,LUO Zhen2
(1.State Key Laboratory of High Performance and Complex Manufacturing,Central South University, Changsha 410083,China; 2.School of Mechanical and Electrical Engineering,Central South University,Changsha 410083,China)

Based on the theory of Timoshenko beam,an analytical model of parameter lumping is established for analyzing the degradation performance of the combined rotors induced by rod creep relaxation.The influence of the stress relaxation of the rod bolts and the contact interface is considered.By using equivalent processing,the dynamic performance of the combined rotor is analyzed,and the regulation of natural frequency change of the combined rotor with the time is obtained.By defining the reduction rate of the combined rotor natural frequency as the degradation parameter,the quantitative evaluation of the combined rotor degradation performance caused by the rod relaxation is realized.The curve of the natural frequency degradation of the combined rotor caused by rod creep relaxation vs.time is obtained.

vibration and wave;rod relaxation;combined rotor;degradation;contact interface

TH113.1

：A

：10.3969/j.issn.1006-1335.2016.06.004

1006-1355(2016)06-0016-05

2016-02-09

劉強(1992-），男，安徽省亳州市人，碩士生，主要研究方向為系統(tǒng)動力學(xué)、轉(zhuǎn)子動力學(xué)。

王艾倫（1959-），男，博士生導(dǎo)師。E-mail:934958469@qq.com