高速渦輪發(fā)電機(jī)轉(zhuǎn)子振動(dòng)特性試驗(yàn)

馮 巖，王績(jī)德，韓東江

（1.中國(guó)中原對(duì)外工程有限公司，北京 100044；2.中國(guó)科學(xué)院 工程熱物理研究所，北京 100190）

渦輪發(fā)電機(jī)是能源系統(tǒng)中熱-功-電轉(zhuǎn)換的核心動(dòng)力裝置，其性能的好壞直接影響能源動(dòng)力系統(tǒng)的經(jīng)濟(jì)性、安全性與穩(wěn)定性。隨著能源系統(tǒng)的小型化與智能化，渦輪發(fā)電機(jī)正向著高速化與微型化的方向發(fā)展，其典型的應(yīng)用場(chǎng)合為冷熱電分布式能源系統(tǒng)與超臨界二氧化碳發(fā)電系統(tǒng)。在上述系統(tǒng)中，對(duì)渦輪發(fā)電機(jī)提出功率密度大、體積小、重量輕、轉(zhuǎn)速高的要求。氣體軸承具有轉(zhuǎn)速高、功耗低、精度高和壽命長(zhǎng)等優(yōu)點(diǎn)[1]，在高速支承、低摩擦損耗支承、高精度支承和特殊工況下的支承領(lǐng)域占有絕對(duì)的應(yīng)用優(yōu)勢(shì)，從而被廣泛應(yīng)用于高速微型動(dòng)力裝置、精密儀器及航空航天等國(guó)防領(lǐng)域[2]。氣體軸承黏度低使得氣體軸承在高速旋轉(zhuǎn)動(dòng)力裝置中應(yīng)用的關(guān)鍵問題之一在于氣體軸承支承下轉(zhuǎn)子的穩(wěn)定性[3]。

國(guó)內(nèi)外學(xué)者對(duì)氣體軸承-轉(zhuǎn)子系統(tǒng)穩(wěn)定性進(jìn)行了大量的理論與試驗(yàn)研究。侯予等[4]將可傾瓦動(dòng)壓氣體軸承應(yīng)用于高速透平膨脹制冷機(jī)中，并對(duì)氣體軸承承載力與穩(wěn)定性開展了一系列研究。楊金福等[5－6]開展氣體軸承支承的高速永磁電機(jī)、高速渦輪膨脹制冷機(jī)軸系動(dòng)力學(xué)特性實(shí)驗(yàn)研究，探討了阻尼墊、軸承供氣壓力等轉(zhuǎn)子穩(wěn)定性調(diào)整因素，并提出軸承轉(zhuǎn)子流-固耦合調(diào)頻原理及工程穩(wěn)定性判別準(zhǔn)則。張廣輝[7－8]通過數(shù)值分析研究動(dòng)靜壓混合氣體軸承支承的轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性，并在旋轉(zhuǎn)沖壓發(fā)動(dòng)機(jī)試驗(yàn)臺(tái)上進(jìn)行氣體軸承-轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)性能研究。韓東江[9]提出基于圖譜分析的轉(zhuǎn)子非線性振動(dòng)行為分析方法，該方法以氣膜振蕩起始點(diǎn)為界，選用不同的圖譜來分析氣膜振蕩前后的動(dòng)力學(xué)特征，給出了氣膜振蕩發(fā)生前不平衡量引起同頻渦動(dòng)的分析圖譜和氣膜振蕩發(fā)生后氣膜非線性引起低頻渦動(dòng)與振蕩的分析圖譜。Su[10－11]通過數(shù)值分析研究多排研究節(jié)流孔供氣徑向軸承的旋轉(zhuǎn)效應(yīng)，結(jié)果表明，軸承承載能力隨偏心率的增加速率大于其隨轉(zhuǎn)速的增加速率，并使用擾動(dòng)法對(duì)多排節(jié)流孔供氣的動(dòng)靜壓混合軸承動(dòng)態(tài)穩(wěn)定性進(jìn)行分析，通過計(jì)算軸承的剛度和阻尼系數(shù)確定其失穩(wěn)轉(zhuǎn)速。Morosi等[12－13]提出了一種可主動(dòng)控制的氣體潤(rùn)滑軸承，采用壓電制動(dòng)裝置控制通過節(jié)流孔空氣的流量與壓力特性，進(jìn)而實(shí)現(xiàn)氣膜動(dòng)力學(xué)特性的在線控制，并建立該氣體潤(rùn)滑軸承的動(dòng)力學(xué)分析模型，開展軸承及轉(zhuǎn)子動(dòng)力學(xué)特性試驗(yàn)研究，驗(yàn)證該軸承的有效性。

本文基于靜壓氣體軸承支承的高速渦輪發(fā)電機(jī)試驗(yàn)裝置，采用轉(zhuǎn)子非線性振動(dòng)分析方法，通過試驗(yàn)研究高速渦輪發(fā)電機(jī)轉(zhuǎn)子升速過程中非線性振動(dòng)特性，并重點(diǎn)分析軸承供氣壓力對(duì)高速渦輪發(fā)電機(jī)轉(zhuǎn)子振動(dòng)特性的影響。

1 高速渦輪發(fā)電機(jī)試驗(yàn)臺(tái)

高速渦輪發(fā)電機(jī)試驗(yàn)臺(tái)[14]如圖1所示，包括軸承供氣及渦輪驅(qū)動(dòng)氣子系統(tǒng)、振動(dòng)信號(hào)采集與處理子系統(tǒng)、控制子系統(tǒng)和高速渦輪發(fā)電機(jī)。

圖1 高速渦輪發(fā)電機(jī)試驗(yàn)臺(tái)結(jié)構(gòu)圖

軸承供氣及渦輪驅(qū)動(dòng)氣子系統(tǒng)由螺桿式空氣壓縮機(jī)提供壓力為1.20 MPa、流量800 Nm3/h的高壓氣源；高壓氣體通過管路支線1進(jìn)入驅(qū)動(dòng)渦輪，為高速渦輪發(fā)電機(jī)提供動(dòng)力，高壓氣體通過管路支線2進(jìn)入靜壓氣體軸承，為軸承提供潤(rùn)滑與支承；在管路支線上布置壓力傳感器、溫度傳感器、過濾器、穩(wěn)壓閥與調(diào)節(jié)閥等，實(shí)現(xiàn)氣體參數(shù)控制與監(jiān)控。

振動(dòng)信號(hào)采集與處理子系統(tǒng)由電渦流位移傳感器、數(shù)據(jù)采集儀、存儲(chǔ)電腦組成。渦輪端布置3個(gè)電渦流位移傳感器以測(cè)量渦輪端轉(zhuǎn)子振動(dòng)位移及旋轉(zhuǎn)速度，壓氣機(jī)端布置2兩個(gè)電渦流位移傳感器以測(cè)量壓氣機(jī)端轉(zhuǎn)子振動(dòng)位移，所測(cè)量振動(dòng)位移信號(hào)通過前置器轉(zhuǎn)換后進(jìn)入數(shù)據(jù)采集儀，并在DASP-V10軟件中實(shí)現(xiàn)在線監(jiān)測(cè)與離線分析。

控制子系統(tǒng)通過調(diào)節(jié)閥控制軸承供氣及驅(qū)動(dòng)渦輪所需高壓氣體的流量，以實(shí)現(xiàn)控制軸承供氣壓力和轉(zhuǎn)子升速率的目的。

高速渦輪發(fā)電機(jī)是實(shí)現(xiàn)熱-功-電轉(zhuǎn)換的關(guān)鍵動(dòng)力裝置，其轉(zhuǎn)子結(jié)構(gòu)如圖2所示，為壓氣機(jī)、透平與發(fā)電機(jī)同軸的一體化轉(zhuǎn)子結(jié)構(gòu)。

圖2 高速渦輪發(fā)電機(jī)一體化轉(zhuǎn)子結(jié)構(gòu)

發(fā)電機(jī)采用永磁盤式電機(jī)結(jié)構(gòu)，在軸承兩側(cè)分別布置一對(duì)磁盤，每個(gè)磁盤上均勻布置8個(gè)永磁體。在徑向上，磁盤通過鍵與轉(zhuǎn)軸連接，在軸向上通過螺母將葉輪與磁盤固定。每對(duì)磁盤間布置一個(gè)定子線圈，其冷卻方式為自然風(fēng)冷，同時(shí)定子線圈上布置2個(gè)溫度測(cè)點(diǎn)，通過定子線圈溫度變化來觀察永磁盤與定子線圈是否發(fā)生碰磨、定子線圈溫度隨發(fā)電負(fù)載功率的變化情況。轉(zhuǎn)子的結(jié)構(gòu)尺寸如表1所示。

表1 轉(zhuǎn)子的結(jié)構(gòu)參數(shù)

轉(zhuǎn)子結(jié)構(gòu)由一對(duì)靜壓氣體軸承支承，軸承為圓柱型節(jié)流孔式徑向-止推混合軸承，軸承材料為石墨合金，具有耐磨與自潤(rùn)滑的功能；采用橡膠“O”型圈提高轉(zhuǎn)子穩(wěn)定性；徑向潤(rùn)滑由2排節(jié)流小孔提供，沿著周向均布分布，數(shù)量為16個(gè)；止推潤(rùn)滑由止推軸承面上均勻布置的10個(gè)節(jié)流小孔提供，以平衡轉(zhuǎn)子軸向推力，其結(jié)構(gòu)尺寸如表2所示。

表2 軸承的結(jié)構(gòu)參數(shù)

2 試驗(yàn)方案

表3中所示為本文所采用的試驗(yàn)方案。在試驗(yàn)1中，采用升速過程變軸承供氣壓力的方案，渦輪端與壓氣機(jī)端軸承以0.67 MPa的軸承供氣壓力啟動(dòng)，當(dāng)轉(zhuǎn)速到達(dá)16600 r/min時(shí)，渦輪端與壓氣機(jī)端軸承供氣壓力增加為0.80 MPa。

表3 軸承供氣壓力方案

轉(zhuǎn)速點(diǎn)16600 r/min的選取是基于多次試驗(yàn)后歸納總結(jié)得到的。由于轉(zhuǎn)子轉(zhuǎn)速在16600 r/min發(fā)生飛升，為了抑制飛升過程中轉(zhuǎn)子低頻振動(dòng)現(xiàn)象，在飛升點(diǎn)前提高軸承供氣壓力至0.80 MPa。在試驗(yàn)2中，升速過程軸承供氣壓力保持恒定，為0.80 MPa。

3 渦輪發(fā)電機(jī)升速過程中轉(zhuǎn)子振動(dòng)特性與發(fā)電特性分析

3.1 試驗(yàn)1渦輪發(fā)電機(jī)轉(zhuǎn)子振動(dòng)特性

圖3給出了試驗(yàn)1軸承在供氣壓力下渦輪端水平方向轉(zhuǎn)子升速過程時(shí)間-轉(zhuǎn)速-幅值三維譜圖，其中橫坐標(biāo)為頻率，單位為Hz，縱坐標(biāo)為幅值，單位為μm。曲線ABCD代表轉(zhuǎn)子升速過程工頻振動(dòng)曲線。

如圖3中所示，轉(zhuǎn)子在B點(diǎn)，轉(zhuǎn)速為11899 r/min時(shí)出現(xiàn)低頻渦動(dòng)現(xiàn)象，渦動(dòng)頻率為146.63 Hz，隨著轉(zhuǎn)子轉(zhuǎn)速增加，渦動(dòng)頻率隨之增加；轉(zhuǎn)速為14279 r/min時(shí)，低頻渦動(dòng)頻率出現(xiàn)跳躍，由153.85 Hz降為132.21 Hz，該低頻頻率跳躍的直接原因?yàn)闅饽偠冉档?，即氣膜的渦動(dòng)分力對(duì)應(yīng)的剛度降低。轉(zhuǎn)子在C點(diǎn)，轉(zhuǎn)速15505 r/min時(shí)低頻振動(dòng)消失，轉(zhuǎn)子出現(xiàn)碰摩特征，轉(zhuǎn)速為16658 r/min時(shí)，轉(zhuǎn)子進(jìn)入飛升（升速率急劇增加）階段，到達(dá)D點(diǎn)，轉(zhuǎn)速42548 r/min時(shí)穩(wěn)定；在轉(zhuǎn)速45780 r/min時(shí)，轉(zhuǎn)子出現(xiàn)低頻振蕩，低頻頻率為158.65 Hz，低頻振蕩的出現(xiàn)引起轉(zhuǎn)子工頻轉(zhuǎn)速的降低，低頻振蕩消失后轉(zhuǎn)子工頻轉(zhuǎn)速回升并穩(wěn)定。

轉(zhuǎn)子在以上各個(gè)階段頻譜圖及軸心軌跡如圖4所示?？梢钥吹剑贏B階段，如圖4(a)所示，轉(zhuǎn)子振動(dòng)譜圖以工頻振動(dòng)為主，軸心軌跡呈現(xiàn)周期一特征；且在7428 r/min(123.80 Hz)時(shí)出現(xiàn)臨界轉(zhuǎn)速特征，該臨界為轉(zhuǎn)子錐動(dòng)臨界轉(zhuǎn)速特征，如圖4(b)所示。在BC階段，轉(zhuǎn)子出現(xiàn)低頻渦動(dòng)后頻譜圖及軸心軌跡如圖4(c)所示，軸心軌跡呈現(xiàn)多周期特征。圖4(d)給出C點(diǎn)后轉(zhuǎn)子轉(zhuǎn)速為15851 r/min時(shí)頻譜特征，可以看到，頻譜分析圖上，出現(xiàn)類似于點(diǎn)碰磨特征的連續(xù)譜，同時(shí)，在時(shí)域頻譜上間斷性出現(xiàn)幅值的突躍-衰減，而軸心軌跡呈現(xiàn)混沌特征，如圖4(e)所示。

圖3 試驗(yàn)1工況下轉(zhuǎn)子升速過程時(shí)間-轉(zhuǎn)速-幅值三維譜圖

圖4(f)給出飛升過程中轉(zhuǎn)子頻譜分析及軸心軌跡，轉(zhuǎn)子在飛升階段升速率為4309 r/min/s。該階段形成的原因?yàn)樵谵D(zhuǎn)速15500 r/min之前轉(zhuǎn)子發(fā)生低頻振動(dòng)，隨著輸入能量（驅(qū)動(dòng)渦輪流量增加，但轉(zhuǎn)子轉(zhuǎn)速增加緩慢）增加，轉(zhuǎn)子在發(fā)生低頻振動(dòng)區(qū)域積蓄能量，當(dāng)?shù)皖l振動(dòng)消失后，轉(zhuǎn)子出現(xiàn)較大升速率的飛升現(xiàn)象。可以看到，在頻譜分析圖上，工頻的頻率值呈現(xiàn)頻帶的特征，其原因是升速率較快，采樣頻率與升速率不匹配。轉(zhuǎn)子在飛升區(qū)域軸心軌跡呈現(xiàn)周期一特征。在飛升轉(zhuǎn)速后期，轉(zhuǎn)子工頻振動(dòng)幅值出現(xiàn)突變特征，如圖4(g)頻譜圖中所示，其原因?yàn)樵陲w升區(qū)域，無外界能量的輸入（此時(shí)驅(qū)動(dòng)流量不增加），轉(zhuǎn)子在飛升初期，幅值是相對(duì)減小的，隨著飛升進(jìn)入后期，在飛升穩(wěn)定時(shí)，幅值相對(duì)增大，其典型軸心軌跡如圖4(h)所示。

圖4(i)給出轉(zhuǎn)子轉(zhuǎn)速45780 r/min時(shí)出現(xiàn)低頻振蕩（振蕩頻率為158.65 Hz）時(shí)的軸心軌跡。

3.2 試驗(yàn)2渦輪發(fā)電機(jī)轉(zhuǎn)子振動(dòng)特性

從試驗(yàn)1軸承供氣壓力下轉(zhuǎn)子升速實(shí)驗(yàn)結(jié)果可以看到，轉(zhuǎn)子在飛升轉(zhuǎn)速16600 r/min前出現(xiàn)低頻振動(dòng)現(xiàn)象，而在轉(zhuǎn)速為16600 r/min時(shí)提高軸承供氣壓力后，轉(zhuǎn)子低頻振動(dòng)消失，且能夠穩(wěn)定運(yùn)行于45000 r/min。本小節(jié)基于試驗(yàn)1結(jié)果，將升速過程中供氣壓力穩(wěn)定在0.80 MPa，觀察供氣壓力對(duì)轉(zhuǎn)子升速過程低頻振動(dòng)的影響。

圖4 試驗(yàn)1工況下升速過程典型區(qū)域轉(zhuǎn)子頻譜特征及軸心軌跡

圖5 試驗(yàn)2工況下轉(zhuǎn)子升速過程時(shí)間-轉(zhuǎn)速-幅值三維譜圖

圖5給出試驗(yàn)2軸承供氣壓力下渦輪端水平方向轉(zhuǎn)子升速過程時(shí)間-轉(zhuǎn)速-幅值三維譜圖，其中橫坐標(biāo)為頻率，單位為Hz，縱坐標(biāo)為幅值，單位為μm。轉(zhuǎn)子在32690 r/min時(shí)出現(xiàn)低頻振蕩現(xiàn)象，低頻頻率為125 Hz，隨著工頻轉(zhuǎn)速增加，低頻頻率略有增加。試驗(yàn)2工況下轉(zhuǎn)子穩(wěn)定運(yùn)行于40770 r/min。

圖6給出試驗(yàn)2工況下轉(zhuǎn)子升速過程中典型軸心軌跡、頻譜特性與臨界轉(zhuǎn)速區(qū)域伯德圖。如圖6(a)所示，轉(zhuǎn)子錐動(dòng)臨界發(fā)生在轉(zhuǎn)速7544 r/min，比試驗(yàn)1工況下臨界轉(zhuǎn)速值高116 r/min，其原因?yàn)檩S承供氣壓力的增高引起氣膜剛度的增加，進(jìn)而提高臨界轉(zhuǎn)速值。在臨界轉(zhuǎn)速區(qū)域，轉(zhuǎn)子以周期一運(yùn)行，其軸心軌跡如圖6(b)所示。圖6(c)給出轉(zhuǎn)子出現(xiàn)低頻振蕩后軸心軌跡及頻譜圖，轉(zhuǎn)速40743 r/min時(shí)，轉(zhuǎn)子以擬周期狀態(tài)運(yùn)行，且軸心軌跡最大值在60 μm以內(nèi)。

圖6 試驗(yàn)2工況下升速過程典型區(qū)域轉(zhuǎn)子頻譜特征及軸心軌跡

分岔圖描述轉(zhuǎn)子升速過程中從周期一穩(wěn)定運(yùn)轉(zhuǎn)到出現(xiàn)混沌失穩(wěn)的路徑。文中的分岔圖是根據(jù)實(shí)驗(yàn)數(shù)據(jù)繪制出來的：按等采樣點(diǎn)采樣，以采樣開始時(shí)轉(zhuǎn)子上鍵相槽產(chǎn)生的鍵相信號(hào)的上升沿為零基準(zhǔn)，分岔圖上每個(gè)點(diǎn)表示鍵相信號(hào)上升沿相對(duì)于零基準(zhǔn)的位移。圖7給出試驗(yàn)2工況下轉(zhuǎn)子升速過程分岔圖，可以看到，分岔點(diǎn)為低頻振蕩出現(xiàn)起始轉(zhuǎn)速點(diǎn)。低頻振蕩出現(xiàn)后轉(zhuǎn)速大于38000 r/min，轉(zhuǎn)子振動(dòng)幅值邊界呈收斂狀態(tài)，根據(jù)工程穩(wěn)定性判別準(zhǔn)則[6]可知，轉(zhuǎn)子在該轉(zhuǎn)速下能夠穩(wěn)定運(yùn)行。

圖7 試驗(yàn)2工況下轉(zhuǎn)子升速過程分岔圖

將試驗(yàn)2工況下轉(zhuǎn)子振動(dòng)數(shù)據(jù)與試驗(yàn)1工況下轉(zhuǎn)子振動(dòng)數(shù)據(jù)對(duì)比分析，可以得到如下結(jié)論：

(1)軸承供氣壓力的增加使得轉(zhuǎn)子錐動(dòng)臨界轉(zhuǎn)速值相應(yīng)增加；

(2)試驗(yàn)2工況下軸承供氣壓力的變化能夠消除轉(zhuǎn)子在試驗(yàn)1工況下升速過程中（轉(zhuǎn)速區(qū)域?yàn)?1899 r/min～15505 r/min）出現(xiàn)的低頻振動(dòng)現(xiàn)象；

(3)試驗(yàn)2工況下軸承供氣壓力的變化抑制低頻振動(dòng)現(xiàn)象，使得轉(zhuǎn)子不存在低頻蓄能區(qū)域，進(jìn)而消除轉(zhuǎn)子飛升區(qū)域；

(4)試驗(yàn)2工況下轉(zhuǎn)子在32690 r/min時(shí)出現(xiàn)低頻振蕩現(xiàn)象，但低頻振蕩幅值較小，隨著轉(zhuǎn)速增加，轉(zhuǎn)子振動(dòng)幅值在安全裕度以內(nèi)，依據(jù)工程穩(wěn)定性判別準(zhǔn)則可知，轉(zhuǎn)子能夠在40000 r/min穩(wěn)定運(yùn)行。

3.3 高速渦輪發(fā)電機(jī)發(fā)電特性分析

圖8給出試驗(yàn)1工況下升速過程中渦輪發(fā)電機(jī)轉(zhuǎn)速-功率-線圈溫度變化曲線。

圖8 試驗(yàn)1工況下轉(zhuǎn)速-功率-線圈溫度特性曲線

可以看到，隨著轉(zhuǎn)速的增加，功率隨之增加，進(jìn)而線圈溫度隨之增加。發(fā)電功率與線圈溫度在32000 r/min～45000 r/min之間存在一個(gè)滯環(huán)。這是由于轉(zhuǎn)子轉(zhuǎn)速的飛升導(dǎo)致功率迅速增加，而飛升的時(shí)間很短，熱量還沒有顯現(xiàn)出來，因此存在線圈溫度落后于發(fā)電功率的滯環(huán)。伴隨轉(zhuǎn)子在高轉(zhuǎn)速下穩(wěn)定運(yùn)行，隨著時(shí)間的積累，熱量增加引起線圈溫度的遲滯增加。

4 結(jié)語(yǔ)

采用時(shí)間-轉(zhuǎn)速-幅值三維譜圖、頻譜圖、時(shí)域分析、軸心軌跡及分岔圖等非線性振動(dòng)測(cè)試與分析方法，對(duì)高速渦輪發(fā)電機(jī)轉(zhuǎn)子升速過程中振動(dòng)特性開展試驗(yàn)研究，主要結(jié)論如下：

(1)渦輪發(fā)電機(jī)轉(zhuǎn)子臨界轉(zhuǎn)速隨軸承供氣壓力的增加而增加。

(2)軸承供氣壓力的增加能夠抑制轉(zhuǎn)子低頻渦動(dòng)現(xiàn)象，去除轉(zhuǎn)子低頻蓄能區(qū)域，進(jìn)而消除轉(zhuǎn)子飛升，提高轉(zhuǎn)子運(yùn)行穩(wěn)定性與安全性。

(3)轉(zhuǎn)子出現(xiàn)低頻振蕩后，只要轉(zhuǎn)子振動(dòng)幅值邊界隨轉(zhuǎn)速增加而收斂，且在安全裕度以內(nèi)，則轉(zhuǎn)子可在低頻振蕩狀態(tài)下安全穩(wěn)定運(yùn)行。