基于固定界面模態(tài)綜合法的轉(zhuǎn)子實(shí)際機(jī)組低壓轉(zhuǎn)子軸承系統(tǒng)非線性動(dòng)力學(xué)

基于固定界面模態(tài)綜合法的轉(zhuǎn)子實(shí)際機(jī)組低壓轉(zhuǎn)子軸承系統(tǒng)非線性動(dòng)力學(xué)

200mw汽輪發(fā)電是我國(guó)最大的兩種服務(wù)裝置。據(jù)統(tǒng)計(jì)，200mw汽輪發(fā)電的低頻振動(dòng)率很低。在低壓動(dòng)力裝置和電機(jī)的壓力容量中經(jīng)常發(fā)生非線性油膜失穩(wěn)定事故。因此，有必要研究該裝置的非線性振動(dòng)失穩(wěn)定的機(jī)理及其限制方法。由于200mw汽輪發(fā)電壓縮式元壓縮式模型中單元數(shù)較少，因此采用直接數(shù)值解算方法計(jì)算的計(jì)算出的計(jì)算時(shí)間相對(duì)較短。在這項(xiàng)工作中，我們使用了基于短橢圓形軸承框架的固定界面模型，并采用了固定界面綜合法，計(jì)算了實(shí)際200mw汽輪發(fā)電壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂壓裂的剖面。目前,考慮近年來(lái)人們對(duì)單盤及具有較少自由度的轉(zhuǎn)子一軸承系統(tǒng)穩(wěn)定性進(jìn)行了研究120轉(zhuǎn)子—0MW汽輪機(jī)低壓—轉(zhuǎn)子軸承系統(tǒng)振動(dòng)方程200MW汽輪機(jī)低壓轉(zhuǎn)子的結(jié)構(gòu)如圖1所示,對(duì)轉(zhuǎn)子—軸承系統(tǒng)采用有限元法建模,得到的系統(tǒng)振動(dòng)方程為:式中:[M]為質(zhì)量陣;[D]為阻尼陣;[K]為剛度陣;ω2軸承直徑d基于短軸承假設(shè),采用修正的圓瓦油膜壓力的表達(dá)式其中:D為軸承直徑;L為軸承長(zhǎng)度.通過(guò)積分(3)中表達(dá)式,可以得到油膜力表達(dá)式如(4)所示.式中:為分析計(jì)算方便,通常將轉(zhuǎn)子一軸承系統(tǒng)的振動(dòng)方程無(wú)量綱化.3系統(tǒng)方程降維對(duì)固定轉(zhuǎn)速下200MW汽輪發(fā)電機(jī)組低壓段轉(zhuǎn)子—軸承系統(tǒng),采用固定界面界面模態(tài)綜合法對(duì)系統(tǒng)方程降維,由表1低壓轉(zhuǎn)子分為27段,28個(gè)節(jié)點(diǎn),有兩個(gè)軸承支承,分別作用在第4和第24兩個(gè)節(jié)點(diǎn)上.首先計(jì)算時(shí)假定在低壓轉(zhuǎn)子的第11個(gè)節(jié)點(diǎn)處有偏心距為5μ3.1轉(zhuǎn)子系統(tǒng)floquet乘子隨潤(rùn)滑動(dòng)力黏度的變化規(guī)律固定轉(zhuǎn)子轉(zhuǎn)速為3780r/min,改變軸承潤(rùn)滑油的動(dòng)力黏度,圖2顯示了轉(zhuǎn)子振動(dòng)的振幅隨潤(rùn)滑油黏度變化的曲線,潤(rùn)滑油黏度在μ=0.016Pa·s～0.030Pa·s時(shí),振動(dòng)一直為同步運(yùn)動(dòng),隨著潤(rùn)滑油黏度繼續(xù)增大,轉(zhuǎn)子的動(dòng)力行為發(fā)生變化,出現(xiàn)分岔現(xiàn)象,潤(rùn)滑油黏度增大到μ=0.030Pa·s,轉(zhuǎn)子開(kāi)始出現(xiàn)油膜振蕩.主導(dǎo)Floquet乘子隨潤(rùn)滑油動(dòng)力黏度變化規(guī)律,如圖3可知,當(dāng)μ=0.016Pa·s時(shí),系統(tǒng)的主導(dǎo)Floquet乘子的模大于1,系統(tǒng)失穩(wěn);在μ=0.030Pa·s時(shí),系統(tǒng)的主導(dǎo)Floquet乘子為-0.939376-0.382486i、-0.939376+0.382486i,模為1.014260大于1,主導(dǎo)Floquet乘子為一對(duì)共軛復(fù)特征乘數(shù)穿過(guò)單位圓,而其他Floquet乘子的模均小于1,此時(shí)系統(tǒng)此時(shí)產(chǎn)生明顯的Hopf分岔.從以上分析可知,適當(dāng)?shù)臏p小潤(rùn)滑油黏度可以提高轉(zhuǎn)子運(yùn)行的穩(wěn)定性,并且利用Floquet理論分析的同一轉(zhuǎn)速系統(tǒng)失穩(wěn)時(shí)的油黏度要小一些,比較準(zhǔn)確地分析系統(tǒng)的穩(wěn)定性和分岔類型.3.2轉(zhuǎn)子振動(dòng)隨偏心距的變化規(guī)律在固定轉(zhuǎn)速為3360r/min時(shí),改變轉(zhuǎn)子在第11個(gè)節(jié)點(diǎn)處的偏心距.如圖4所示,當(dāng)偏心距e=0.2～1.7mm范圍內(nèi)變化時(shí),轉(zhuǎn)子的振動(dòng)主要受不平衡力的影響.增大偏心距,轉(zhuǎn)子的振動(dòng)發(fā)生了倍周期分岔,出現(xiàn)半速渦動(dòng),圖4顯示了轉(zhuǎn)子振動(dòng)的振幅隨偏心距變化的規(guī)律,當(dāng)偏心距達(dá)到e=1.8mm時(shí),油膜渦動(dòng)轉(zhuǎn)變?yōu)橛湍ふ袷?轉(zhuǎn)子的振幅急劇增加.圖5顯示了主導(dǎo)Floquet乘子隨偏心距變化規(guī)律,當(dāng)e=1.8mm時(shí),系統(tǒng)的主導(dǎo)Floquet乘子為-1.012993+0.000000i,模為1.012993為大于1,系統(tǒng)失穩(wěn);e=3.0mm時(shí)系統(tǒng)的Floquet乘子,主導(dǎo)Floquet乘子從(-1,0)穿過(guò)單位圓,而其他Floquet乘子的模均小于1,此時(shí)系統(tǒng)產(chǎn)生了倍周期分岔.3.3轉(zhuǎn)子失穩(wěn)轉(zhuǎn)速隨軸承長(zhǎng)徑比的變化規(guī)律圖6為在轉(zhuǎn)速3600r/min時(shí)轉(zhuǎn)子的振幅隨軸承的長(zhǎng)徑比變化規(guī)律,當(dāng)軸承的長(zhǎng)徑比小于0.80時(shí),轉(zhuǎn)子的振動(dòng)幅值較小,均在100μm以下;當(dāng)軸承的長(zhǎng)徑比大于0.80時(shí)轉(zhuǎn)子的振動(dòng)幅值會(huì)增大得很快,開(kāi)始出現(xiàn)油膜失穩(wěn)現(xiàn)象,可見(jiàn),對(duì)于短圓柱軸承,適當(dāng)?shù)臏p小長(zhǎng)徑比,可以提高失穩(wěn)轉(zhuǎn)速,以保證轉(zhuǎn)子在正常工作轉(zhuǎn)速下不發(fā)生失穩(wěn).圖7顯示了主導(dǎo)Floquet乘子隨軸承長(zhǎng)徑比變化規(guī)律,當(dāng)L/D=0.80時(shí),系統(tǒng)的主導(dǎo)Floquet乘子為-0.981266-0.244307i、-0.981266+0.244307i,模為1.011221為大于1,系統(tǒng)失穩(wěn).L/D=0.85時(shí)主導(dǎo)Floquet乘子為一對(duì)共軛復(fù)特征乘數(shù)穿過(guò)單位圓,而其他Floquet乘子的模均小于1,此時(shí)系統(tǒng)產(chǎn)生了明顯的Hopf分岔.4轉(zhuǎn)子失穩(wěn)轉(zhuǎn)速隨參數(shù)變化的規(guī)律本文提出利用Floquet理論分析降維計(jì)算得到的轉(zhuǎn)子—軸承系統(tǒng)的主導(dǎo)Floquet乘子,當(dāng)主導(dǎo)Floquet乘子的模大于1時(shí),得到的失穩(wěn)轉(zhuǎn)速作為判斷系統(tǒng)穩(wěn)定性準(zhǔn)則.采用轉(zhuǎn)子參數(shù)及圓柱瓦軸承參數(shù)如表3所示,通過(guò)改變軸承的潤(rùn)滑油動(dòng)力黏度、轉(zhuǎn)子的偏心距、軸承間隙比、軸承長(zhǎng)徑比來(lái)計(jì)算轉(zhuǎn)子發(fā)生油膜振蕩的失穩(wěn)轉(zhuǎn)速,為汽輪發(fā)電機(jī)組的設(shè)計(jì)和運(yùn)行提供了定量和定性的理論依據(jù).4.1失穩(wěn)轉(zhuǎn)速降低說(shuō)改變低壓轉(zhuǎn)子前軸承的潤(rùn)滑油動(dòng)力黏度,如圖8可見(jiàn),隨著轉(zhuǎn)子前軸承潤(rùn)滑油的動(dòng)力黏度的增加,失穩(wěn)轉(zhuǎn)速降低,即振動(dòng)穩(wěn)定性隨軸承潤(rùn)滑油的動(dòng)力黏度的增加而降低.降低油黏度增強(qiáng)系統(tǒng)穩(wěn)定性的一種方法是提高油溫度,當(dāng)平均油溫很高時(shí),繼續(xù)提高進(jìn)油溫度可提高溫定性,但在實(shí)際系統(tǒng)運(yùn)行中,對(duì)油溫度是有一定限制的.由Floquet理論分析判斷,在失穩(wěn)轉(zhuǎn)速時(shí)分岔類型均為Hopf分岔.4.2失穩(wěn)轉(zhuǎn)速隨轉(zhuǎn)子t理論變化的規(guī)律改變轉(zhuǎn)子在第11個(gè)節(jié)點(diǎn)的偏心距,計(jì)算得到轉(zhuǎn)子前軸承軸頸在不同偏心距下的失穩(wěn)轉(zhuǎn)速,如圖9所示.可見(jiàn),系統(tǒng)在e=2.0mm以內(nèi),低壓轉(zhuǎn)子上的不平衡量的大小對(duì)發(fā)電機(jī)前軸承的穩(wěn)定性影響并不大,隨著不平衡量的增大,失穩(wěn)轉(zhuǎn)速降低.當(dāng)e=2.1mm時(shí),失穩(wěn)轉(zhuǎn)速突然降低,即較大的不平衡量可以激起油膜振蕩,不利于系統(tǒng)的穩(wěn)定.由Floquet理論分析判斷,失穩(wěn)轉(zhuǎn)速隨偏心距變化時(shí)分岔類型,當(dāng)e=0.1～2.0mm范圍內(nèi),在失穩(wěn)轉(zhuǎn)速時(shí)系統(tǒng)分岔類型為Hopf分岔,當(dāng)e=2.1～3.0mm范圍內(nèi)時(shí)在失穩(wěn)轉(zhuǎn)速時(shí)系統(tǒng)產(chǎn)生倍周期分岔.4.3軸承間隙比對(duì)失穩(wěn)轉(zhuǎn)速的影響改變發(fā)電機(jī)前軸承的軸承間隙比,計(jì)算得到的低壓轉(zhuǎn)子前軸承軸頸的失穩(wěn)轉(zhuǎn)速如圖10所示長(zhǎng)徑比為0.58時(shí)的變化曲線,隨軸承間隙比增加,失穩(wěn)轉(zhuǎn)速先降低,在間隙比為1.4‰時(shí),轉(zhuǎn)子的失穩(wěn)轉(zhuǎn)速最低,穩(wěn)定性最差;繼續(xù)增加間隙比,轉(zhuǎn)子的失穩(wěn)轉(zhuǎn)速隨之不斷增高,穩(wěn)定性增強(qiáng).由Floquet理論分析判斷,在失穩(wěn)轉(zhuǎn)速時(shí)系統(tǒng)分岔類型均為Hopf分岔.4.4失穩(wěn)轉(zhuǎn)速隨軸承長(zhǎng)徑比變化的規(guī)律改變發(fā)電機(jī)前軸承的軸承長(zhǎng)徑比,計(jì)算得到的失穩(wěn)轉(zhuǎn)速如圖11所示,當(dāng)軸承長(zhǎng)徑比小于0.45時(shí),隨軸承長(zhǎng)徑比增加,失穩(wěn)轉(zhuǎn)速升高,當(dāng)軸承長(zhǎng)徑比大于0.45時(shí),隨軸承長(zhǎng)徑比增加,失穩(wěn)轉(zhuǎn)速降低.由Floquet理論分析判斷,當(dāng)L/D=0.15～0.45時(shí),在失穩(wěn)轉(zhuǎn)速時(shí)系統(tǒng)分岔類型為倍周期分岔,當(dāng)L/D=0.45～1.0時(shí),在失穩(wěn)轉(zhuǎn)速時(shí)系統(tǒng)產(chǎn)生Hopf分岔.5失穩(wěn)轉(zhuǎn)速隨潤(rùn)滑動(dòng)力黏度的變化規(guī)律采用固定界面模態(tài)綜合法結(jié)合Runge-Kutta數(shù)值積分法的降維計(jì)算方法,并利用Floquet理論進(jìn)行對(duì)比分析計(jì)算出在固定轉(zhuǎn)速和給定轉(zhuǎn)子—軸承參數(shù)下,圓柱瓦軸承支撐的轉(zhuǎn)子振動(dòng)穩(wěn)定性隨軸承的潤(rùn)滑油動(dòng)力黏度、轉(zhuǎn)子偏心距、軸承長(zhǎng)徑比、軸承間隙比的變化規(guī)律;并分析了低壓轉(zhuǎn)子前軸承軸頸的失穩(wěn)轉(zhuǎn)速隨軸承的潤(rùn)滑油動(dòng)力黏度、轉(zhuǎn)子偏心距、軸承長(zhǎng)徑比、軸承間隙比的變化規(guī)律.計(jì)