李彥軍 陳 超 裴 吉 王文杰 吳天澄
(江蘇大學(xué)流體機械工程技術(shù)研究中心, 鎮(zhèn)江 212013)
壓力脈動是影響泵站安全穩(wěn)定的重要因素之一[1]。目前國內(nèi)外對軸流泵的壓力脈動[2-4]已進行了較多的研究,但對泵裝置[5]空化工況下的壓力脈動研究較少。張德勝等[6]在軸流泵葉輪段和導(dǎo)葉段外壁面布置了 6 個壓力脈動監(jiān)測點,對多個流量工況的壓力脈動進行了動態(tài)測量,揭示了軸流泵內(nèi)部不同位置處壓力脈動規(guī)律,發(fā)現(xiàn)由于受到葉輪內(nèi)壓力梯度的交替變化影響,葉輪進口、葉輪中間和葉輪出口監(jiān)測點的時域脈動周期與葉片旋轉(zhuǎn)周期一致,在小流量工況下葉輪內(nèi)部渦流誘導(dǎo)了明顯的二次諧波,導(dǎo)葉進口、導(dǎo)葉中間和導(dǎo)葉出口監(jiān)測點的壓力脈動頻域中出現(xiàn)了撞擊和回流誘導(dǎo)的低頻信號,同時也存在葉輪的主頻及其諧頻。陸榮等[7]采用高精度數(shù)字壓力采集系統(tǒng),對軸流泵模型的葉輪進口、導(dǎo)葉流道內(nèi)和導(dǎo)葉出口進行壓力測試,發(fā)現(xiàn)受葉輪與導(dǎo)葉的動靜干涉影響,導(dǎo)葉內(nèi)部及出口均出現(xiàn)了小峰值的二次諧波。在非穩(wěn)定工況下各點的時域信號均出現(xiàn)較大峰值的二次諧波,由于回流和葉頂泄漏渦等因素的影響,各監(jiān)測點的頻率成分復(fù)雜,脈動幅值大于其余工況。楊帆等[8]開展了立式軸流泵裝置非定常數(shù)值模擬,通過設(shè)置監(jiān)測點獲得了不同位置的壓力脈動結(jié)果,并進行了頻譜分析,發(fā)現(xiàn)大流量工況時葉輪進口的脈動相對幅值較大,且渦帶誘發(fā)的脈動主頻與軸頻相同,各工況下箱涵式進水流道底部水力脈動均為低頻渦帶脈動,大流量工況時,箱涵式進水流道底部的脈動壓力主頻與軸頻相同,次主頻與葉頻相同,扭矩變化主頻與葉輪的軸頻相同,而渦帶工況時葉輪進口的脈動主頻已不再受軸頻控制。
由于軸流泵試驗系統(tǒng)龐大,試驗成本較高[6],目前國內(nèi)外主要采用數(shù)值模擬研究軸流泵壓力脈動[9-13]。但是數(shù)值模擬往往誤差較大,尤其是在非設(shè)計點工況和發(fā)生空化的情況下[14-15],模擬結(jié)果與試驗有較大的差別。本文針對派河口泵站軸流泵站裝置模型在不同空化工況下進行試驗和分析,揭示非空化和發(fā)生空化工況下軸流泵裝置內(nèi)部的壓力脈動特性。
壓力脈動試驗在江蘇大學(xué)流體機械工程技術(shù)研究中心水力機械四象限多功能試驗臺進行。試驗裝置示意圖如圖1所示,試驗設(shè)備主要由高壓罐、低壓罐、電磁流量計、蝶閥、循環(huán)泵等組成;試驗儀器主要包括:智能電磁流量計1支、EJA型智能差壓變送器1支、EJA型智能絕對壓力變送器1支、JCL2/500N·m智能型轉(zhuǎn)矩轉(zhuǎn)速傳感器 1支、CY200 型數(shù)字壓力傳感器3支、485-20型集線器1臺。數(shù)據(jù)采集系統(tǒng)采用成都泰斯特公司研制的系統(tǒng),其融合了高精密度、信號采集處理、通訊、總線等技術(shù)。
試驗?zāi)P捅眯阅軈?shù)為:設(shè)計揚程Hd=4.5 m,額定流量Qd=303.7 L/s,額定轉(zhuǎn)速n=1 430 r/min,比轉(zhuǎn)數(shù)ns=887;葉片數(shù)z=3,導(dǎo)葉葉片數(shù)zd=6,葉輪直徑D2=300 mm,葉頂間隙C=0.15 mm。模型泵裝置結(jié)構(gòu)圖如圖2所示,性能曲線如圖3所示。
圖1 泵試驗臺布置圖Fig.1 Layout of pump test system1.穩(wěn)流罐 2.高壓罐 3.蝶閥 4.電磁流量計 5.循環(huán)泵 6.低壓罐 7.模型泵裝置
圖2 模型泵裝置結(jié)構(gòu)示意圖Fig.2 Schematic diagram of structure of model pump device1.肘形進水流道 2.葉輪 3.導(dǎo)葉 4.泵軸 5.駝峰形出水流道
圖3 模型泵裝置效率與揚程實測特性曲線Fig.3 Measurement characteristic curves of efficiency and head of model pump device
試驗依據(jù)GB/T 17189—2007《水力機械振動和脈動現(xiàn)場測試規(guī)程》進行。壓力脈動監(jiān)測點選擇了葉輪進口C1、葉輪出口C2、導(dǎo)葉出口C3共3個壓力脈動監(jiān)測點,如圖4所示。
圖4 監(jiān)測點位置示意圖Fig.4 Diagram of monitoring point positions
試驗對0°葉片安放角采集了在2.5、3.5、4.5、5.4 m揚程下分別在未發(fā)生空化(效率無下降)、臨界空化(效率下降1%)、深度空化(效率下降3%)3種空化工況附近的壓力脈動信號,共12組數(shù)據(jù),采樣頻率為2 000 Hz,采樣時間為15 s。
圖5 不同工況下壓力脈動時域圖Fig.5 Time domain diagrams of pressure pulsation under different working conditions
各監(jiān)測點在不同工況下的壓力脈動時域圖如圖5所示。其中,壓力脈動幅值以壓力脈動系數(shù)Cp表征,定義為
(1)
式中p——瞬時壓力,Pa
ρ——流體密度,kg/m3
u2——葉頂圓周速度,m/s
定義葉輪旋轉(zhuǎn)周期數(shù)
N=t/T
(2)
式中t——信號時間長度,s
T——葉輪旋轉(zhuǎn)一周的時間,s
本文對葉輪旋轉(zhuǎn)6個周期的壓力脈動值作時域分析。
葉輪進口C1處的壓力脈動曲線為平滑穩(wěn)定的曲線,一個周期內(nèi)有2~3個主波峰和波谷;除低揚程工況下壓力脈動曲線近似正弦曲線外,其余工況下均有二次波峰的出現(xiàn),但是二次波峰的出現(xiàn)沒有明顯的規(guī)律性,有的工況甚至出現(xiàn)了一個周期內(nèi)有的主波峰沒有諧波而其余主波峰卻有二次波峰的情況,比如在4.5 m揚程未發(fā)生空化的工況下,這說明流場受葉輪旋轉(zhuǎn)的影響較??;在未發(fā)生空化和深度空化工況下,隨著揚程的升高,諧波出現(xiàn)的趨勢逐漸明顯,具體表現(xiàn)為:一個周期內(nèi)帶有二次波峰的主波出現(xiàn)的個數(shù)增加,二次波峰與主波的差別如幅值大小的差異、波形的陡峭程度等逐漸明顯,最終諧波明顯生成,如5.4 m揚程未發(fā)生空化的工況,此時每個周期內(nèi)有3個主波,每個主波帶有一個二次諧波,呈現(xiàn)明顯的規(guī)律性。原因可能是隨著流量的減小,葉輪進口處的脫流回流造成漩渦的現(xiàn)象逐漸嚴重,而漩渦破滅形成的壓力波是引起諧波的主要原因;隨著流量的減小最終在葉輪進口區(qū)域附近形成穩(wěn)定的漩渦區(qū)從而造成二次諧波周期性有規(guī)律的出現(xiàn)。
葉輪出口C2處的壓力脈動幅值與葉輪進口相比在低揚程工況時明顯增大且一個周期內(nèi)有明顯的3個主波峰和波谷,表明受葉輪旋轉(zhuǎn)作用明顯;但是,諧波出現(xiàn)的情況與葉輪進口類似:每個周期內(nèi)諧波的個數(shù)沒有明顯的規(guī)律性,有的工況壓力脈動曲線為平滑的未出現(xiàn)諧波的曲線,如2.5 m揚程發(fā)生臨界空化的工況,有的工況每個周期內(nèi)出現(xiàn)了穩(wěn)定周期性的主波和二次波峰,如5.4 m揚程未發(fā)生空化的工況,這說明葉輪出口處壓力脈動諧波的產(chǎn)生和葉輪與導(dǎo)葉之間動靜干涉的影響關(guān)系不大,空化產(chǎn)生的渦以及空泡的潰滅對諧波的產(chǎn)生影響較大,泄漏渦的影響在文獻[7]也有提及,但文獻[7]并沒有說明空化渦對諧波的影響;在同一種揚程下,隨著空化情況的加深,每個周期內(nèi)波峰峰頂和波谷谷底的形狀由尖銳向圓滑趨勢發(fā)展,且諧波的趨勢越來越不明顯,即二次波峰與主波的幅值、形狀等的差異越來越小。
導(dǎo)葉出口C3處壓力脈動曲線均為平滑的曲線,一個周期內(nèi)主波峰和波谷的數(shù)量為2~4個,每個周期內(nèi)諧波的出現(xiàn)也無明顯的規(guī)律性,有的工況完全沒有諧波出現(xiàn),如4.5 m揚程臨界空化工況,有的工況一個周期內(nèi)的3個主波峰卻只有一個出現(xiàn)諧波,如2.5 m揚程發(fā)生深度空化工況;在同一個揚程工況下,隨著空化余量的降低,導(dǎo)葉出口壓力脈動曲線的幅值逐漸增加,具體表現(xiàn)為:在未發(fā)生空化時,導(dǎo)葉出口處的壓力脈動幅值比葉輪進口處略低,在4.5 m揚程和5.4 m揚程工況時表現(xiàn)得更加明顯;在空化余量下降到臨界空化時,導(dǎo)葉出口處的壓力脈動曲線幅值與葉輪進口處大小相似;在空化余量下降到深度空化時,導(dǎo)葉出口處的壓力脈動幅值比葉輪進口處大,但仍比葉輪出口處小。
為了對各個測點壓力脈動程度有一個總體估量,統(tǒng)計了各工況下C1~C3監(jiān)測點 97%置信度壓力脈動系數(shù)峰峰值,如圖6所示。在葉輪進口C1處,深度空化工況下,各個揚程下的峰峰值均最小且隨著揚程的升高緩慢增長;未發(fā)生空化時各個揚程下的壓力脈動系數(shù)峰峰值均比深度空化時大,而發(fā)生臨界空化時除了最低揚程2.5 m工況外,其余工況下的峰峰值均為最大;葉輪出口C2處不同工況下的峰峰值除了2.5 m深度空化的工況下,各個工況下的峰峰值均比導(dǎo)葉出口C3處大,原因可能是旋轉(zhuǎn)的葉輪與靜止的導(dǎo)葉形成的動靜干涉作用對流體的壓力脈動影響所致;導(dǎo)葉出口C3處的峰峰值變化情況在深度空化和未發(fā)生空化時與葉輪進口C1處相反:深度空化時各個工況下的峰峰值均比未發(fā)生空化時大。
圖6 各監(jiān)測點97%置信度壓力脈動系數(shù)峰峰值對比Fig.6 Comparisons of peak to peak values with 97% confidence of each monitoring point
為了更好地分析軸流泵裝置內(nèi)部的壓力脈動規(guī)律,對時域信號進行快速傅里葉變換(FFT)并進行頻域分析,如圖7所示。為了便于分析,定義葉片轉(zhuǎn)頻倍數(shù)f/fbpf表示采樣周期
(3)
式中f——脈動數(shù)據(jù)經(jīng)過FFT變換后對應(yīng)的頻率,Hz
fbpf——葉片通過頻率,Hz
從圖7中可以看出,監(jiān)測點C1、C2和C3在不同的揚程和空化工況下壓力脈動的主頻均在葉片通過頻率整數(shù)倍頻附近,有的工況為1倍葉頻,有的工況為2倍葉頻。在同一種空化情況下,隨著揚程的增加,葉輪進口C1處的壓力脈動諧頻逐漸明顯,具體表現(xiàn)為:頻帶變寬、高階諧波個數(shù)增加及其幅值增大,原因可能是隨著流量的減小、揚程的增加,葉輪進口處流體不穩(wěn)定造成的脫流回流等不良流態(tài)越來越嚴重,流體對壁面不同程度的撞擊形成不同頻率的壓力波從而使諧波的頻帶變寬;在同一個揚程工況下,隨著空化余量的降低,葉輪進口C1處的壓力脈動諧波趨勢逐漸減弱,具體表現(xiàn)為:高階諧頻的減少及其幅值的降低并且主頻向低頻移動,4.5 m揚程下未發(fā)生空化和臨界空化兩種工況形成明顯的對比。
圖7 不同工況下壓力脈動頻域圖Fig.7 Frequency domain diagrams of pressure pulsation under different working conditions
葉輪出口C2處的壓力脈動的諧頻現(xiàn)象與其他兩個監(jiān)測點相比最明顯:在未發(fā)生空化時,各個揚程工況下葉輪出口壓力脈動諧波主頻均出現(xiàn)在2倍葉頻附近,當空化余量下降至達到臨界空化和深度空化工況時,除5.4 m揚程工況外其他工況的主頻均出現(xiàn)在1倍葉頻附近,而高揚程5.4 m下的不同空化工況壓力脈動主頻卻都出現(xiàn)在2倍葉頻附近,表明葉頻是影響壓力脈動頻率的主要因素;在同一個揚程工況下時,隨著空化余量的降低,主頻及諧頻也向低頻區(qū)域移動,原因可能是空化情況的加劇造成了更多空泡的破滅,而空泡的破滅產(chǎn)生了較多的低頻信號,這與文獻[6]的結(jié)果有相似之處。但文獻[6]并沒有分析2倍葉頻與3倍葉頻出現(xiàn)的原因,在本文中:在不同的工況下,在發(fā)生臨界空化和深度空化時,葉輪進口C1處幾乎沒有出現(xiàn)2倍葉頻,在葉輪出口C2和導(dǎo)葉出口C3處出現(xiàn)了,原因可能是葉輪與導(dǎo)葉的動靜干涉,動靜干涉的影響不僅出現(xiàn)了2倍葉頻,3倍葉頻也時有出現(xiàn);從圖中可以看出C3處的2倍葉頻與3倍葉頻均比C2處弱,原因可能是動靜干涉產(chǎn)生的壓力波在經(jīng)過導(dǎo)葉的整流作用后并沒有完全消失,而是通過流道壁等因素傳遞到了C3處,故幅值有衰減。
導(dǎo)葉出口C3處的壓力脈動情況與葉輪進口C1處相似:在2.5 m和3.5 m揚程各空化工況下,主頻均出現(xiàn)在1倍葉頻附近;4.5 m揚程未發(fā)生空化以及5.4 m揚程未發(fā)生空化和臨界空化的工況下,出現(xiàn)了高階諧波頻率;導(dǎo)葉出口高階諧頻與葉輪進口有一定的類似關(guān)系,即:葉輪進口出現(xiàn)諧頻的頻率附近,導(dǎo)葉出口處也出現(xiàn)了諧頻,原因可能是葉輪進口和導(dǎo)葉出口諧頻的出現(xiàn)與泵段外部因素有關(guān),比如輔助泵的影響、管道共振等。
(1)低揚程工況下一個旋轉(zhuǎn)周期內(nèi)主波峰和波谷數(shù)為2~4個,而高揚程工況下時域圖中一個旋轉(zhuǎn)周期內(nèi)有3個周期性的主波峰和波谷,與葉片數(shù)相一致,表明葉輪的旋轉(zhuǎn)作用對葉輪進口流場有較大的影響,并且二次諧波有規(guī)律地周期性出現(xiàn)。
(2)除個別工況外,葉輪出口處壓力脈動時域圖中一個周期內(nèi)并未出現(xiàn)明顯的有規(guī)律性的諧波,表明葉輪與導(dǎo)葉的強動靜干涉作用并不明顯。
(3)同一種揚程工況下,葉輪進口處幅值最小,隨著流體向出口流動,導(dǎo)葉出口處壓力脈動幅值隨著空化余量的減小而逐漸增加,在深度空化時幅值大于葉輪進口,葉輪出口處脈動始終最強烈,幅值最大。