操瑞嘉，孔祥序，周以松，孫靜如

(1. 江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013； 2. 中工工程機(jī)械成套有限公司，北京 100070)

為提高艦艇的隱身性能，軍方對(duì)艦船用泵等動(dòng)設(shè)備的振動(dòng)與噪聲指標(biāo)提出了更高的要求.磁力驅(qū)動(dòng)離心泵具有無泄漏的特性，在艦船中有一定應(yīng)用[1].磁力驅(qū)動(dòng)離心泵因其采用水導(dǎo)軸承，機(jī)械作用引起的振動(dòng)與噪聲較常規(guī)的機(jī)械密封離心泵低，而水體作用引起的振動(dòng)和噪聲與常規(guī)的機(jī)械密封離心泵相當(dāng)，是產(chǎn)生噪聲的主要因素.

離心泵工作時(shí)，葉輪與蝸殼隔舌產(chǎn)生動(dòng)靜耦合作用，這將影響泵腔內(nèi)的壓力脈動(dòng).壓力脈動(dòng)幅值隨不同流量工況改變，這是影響泵機(jī)組振動(dòng)與噪聲以及泵機(jī)組運(yùn)行穩(wěn)定性與可靠性的主要因素之一.

國內(nèi)外學(xué)者針對(duì)壓力脈動(dòng)及其引發(fā)的誘導(dǎo)噪聲問題進(jìn)行了大量研究.張帆等[2]針對(duì)雙蝸殼式離心泵研究發(fā)現(xiàn)，內(nèi)圈流道的壓力脈動(dòng)強(qiáng)于外圈流道，且隔舌處出現(xiàn)壓力脈動(dòng)較大值，大流量工況下雙蝸殼隔舌和出口產(chǎn)生一定回流導(dǎo)致蝸殼內(nèi)該處附近的監(jiān)測點(diǎn)壓力脈動(dòng)先減小后增大；江偉等[3]提出離心泵葉輪出口邊傾斜布置會(huì)對(duì)單蝸殼泵的壓力脈動(dòng)幅值產(chǎn)生影響，但隔舌位置區(qū)域的壓力脈動(dòng)的周期數(shù)不會(huì)發(fā)生改變；代翠等[4]提出通過將葉片傾斜布置可有效改善泵作為透平時(shí)流激噪聲；SPENCE等[5]對(duì)離心泵幾個(gè)主要參數(shù)對(duì)壓力脈動(dòng)的影響進(jìn)行了評(píng)估；BERTEN等[6]針對(duì)非設(shè)計(jì)工況下高速泵上葉輪出口邊射流尾跡引起的壓力脈動(dòng)現(xiàn)象進(jìn)行了相應(yīng)研究，表明葉輪與蝸殼間因動(dòng)靜干涉引起的壓力脈動(dòng)現(xiàn)象可歸因于葉輪流道出口邊的流動(dòng)分離現(xiàn)象.

盡管已有文獻(xiàn)[3]研究表明采用傾斜葉片出口邊可有效降低單蝸殼泵隔舌位置處的壓力脈動(dòng)效應(yīng)，并且在文獻(xiàn)[4]中標(biāo)明，傳統(tǒng)出口邊垂直的葉片變?yōu)閮A斜時(shí)，會(huì)因葉片掃掠隔舌位置處“線與線”變?yōu)椤包c(diǎn)與點(diǎn)”的位置關(guān)系，進(jìn)而降低壓力脈動(dòng)及其引發(fā)的誘導(dǎo)噪聲，但對(duì)傾斜葉片出口邊在雙蝸殼泵內(nèi)的具體性能缺少詳盡的研究.

雙蝸殼型離心泵能較好地平衡徑向力，在流量變動(dòng)范圍寬、揚(yáng)程相對(duì)較高的工作環(huán)境下運(yùn)行具有明顯的優(yōu)勢，即使在損管條件下也能應(yīng)急運(yùn)行.因此，雙蝸殼型磁力驅(qū)動(dòng)離心泵在船舶上有較好的推廣應(yīng)用前景[7].

文中以雙蝸殼泵為計(jì)算模型，針對(duì)葉輪出口邊與前后蓋板基本面垂直的葉輪(以下簡稱直葉片葉輪)以及與前后蓋板基本面傾斜的葉輪(以下簡稱斜葉片葉輪)在雙蝸殼內(nèi)的性能進(jìn)行對(duì)比研究，初步揭示傾斜葉片在雙蝸殼泵的壓力脈動(dòng)特性，為船用無泄漏磁力泵應(yīng)用雙蝸殼水力模型實(shí)現(xiàn)減振降噪提供一定理論依據(jù).

1 計(jì)算模型

雙蝸殼泵的主要設(shè)計(jì)參數(shù)分別為額定流量Qd=140 m3/h，額定揚(yáng)程Hd=40 m，泵轉(zhuǎn)速n=3 000 r/min.直葉片根據(jù)原有中南水力模型SM-130-10比轉(zhuǎn)數(shù)為ns=133.1 的泵進(jìn)行換算設(shè)計(jì)得到.斜葉片可參照文獻(xiàn)[8]給出的設(shè)計(jì)方法進(jìn)行設(shè)計(jì).2種葉片結(jié)構(gòu)參數(shù)如表1所示，表中D1為葉輪進(jìn)口直徑，D2為葉輪出口直徑，Φ為葉片包角，β1為葉片進(jìn)口安放角，β2為葉片中間流線出口安放角，Z為葉片數(shù)，Dh為葉輪輪轂直徑，b2為葉片出口寬度.

表1 葉片的主要幾何參數(shù)

使用的同一種蝸殼主要幾何參數(shù)分別為基圓直徑D4=210.0 mm，進(jìn)口直徑D5=125.0 mm，隔舌安放角φ=32.7°，蝸殼斷面為梨形.

采用三維造型軟件Pro/E繪制水體模型如圖1所示.計(jì)算水體區(qū)域包括進(jìn)口段、簡化后的葉輪水體域、間隙域、蝸殼域以及加長的出口段.

圖1 三維水體圖

為清晰表明所研究的直出口邊與傾斜出口邊葉輪的差別，單獨(dú)顯示出2種葉輪的水體圖見圖2.

圖2 葉輪出口水體圖

2 CFD數(shù)值模擬

2.1 前處理?xiàng)l件設(shè)定

應(yīng)用ANSYS 18.0下的CFX模組，選用RNGk-ε湍流模型，對(duì)船用雙蝸殼型磁力驅(qū)動(dòng)離心泵進(jìn)行數(shù)值計(jì)算.非定常計(jì)算時(shí)選用的對(duì)流項(xiàng)為“High Re-solution”，瞬態(tài)項(xiàng)采用“Second Order Backward Euler”.由于泵運(yùn)行時(shí)在第6—8圈時(shí)才達(dá)到穩(wěn)定，因此選擇第8圈計(jì)算結(jié)果進(jìn)行分析.

考慮到泵在實(shí)際運(yùn)行時(shí)因軸端負(fù)載而造成的“掉轉(zhuǎn)”現(xiàn)象，計(jì)算時(shí)設(shè)置的轉(zhuǎn)速為n=2 950 r/min.設(shè)置葉輪每旋轉(zhuǎn)3°為一個(gè)計(jì)算時(shí)間步長，則每一步的計(jì)算時(shí)間為1.695×10-4s，總的計(jì)算步數(shù)為960，選取第840—960步進(jìn)行分析.計(jì)算時(shí)選用總壓進(jìn)口、質(zhì)量流量出口的邊界條件，參考?jí)毫υO(shè)為0.1 MPa.考慮泵過流件實(shí)際鑄造質(zhì)量(過流件表面粗糙度)，設(shè)置壁面粗糙度為0.05 mm.

2.2 網(wǎng)格無關(guān)性

應(yīng)用ANSYS 18.0軟件的ICEM CFD模組進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分.在網(wǎng)格劃分過程中，對(duì)流體的邊界層處進(jìn)行網(wǎng)格局部加密，如圖3所示.

圖3 網(wǎng)格劃分

網(wǎng)格劃分時(shí)，試驗(yàn)了2大類共10種網(wǎng)格劃分方案，用以驗(yàn)證計(jì)算精度.表2為網(wǎng)格劃分的具體方案，分別給出了2種葉輪形式的網(wǎng)格數(shù).

表2 網(wǎng)格數(shù)目

對(duì)網(wǎng)格數(shù)進(jìn)行無關(guān)性驗(yàn)證，計(jì)算不同網(wǎng)格數(shù)N時(shí)泵的揚(yáng)程，如圖4所示.可以看出，當(dāng)網(wǎng)格總數(shù)大于140萬時(shí)，網(wǎng)格數(shù)對(duì)計(jì)算結(jié)果影響可以忽略不計(jì).因此，選擇方案4的網(wǎng)格數(shù)進(jìn)行數(shù)值計(jì)算.

2.3 監(jiān)測點(diǎn)布置

為研究泵內(nèi)的壓力脈動(dòng)變化情況，在雙蝸殼上布置9個(gè)監(jiān)測點(diǎn)，各監(jiān)測點(diǎn)位置如圖5所示.其中，在蝸殼截面的基圓處順時(shí)針布置了監(jiān)測點(diǎn)a—f；在Ob向外射線與蝸殼隔板內(nèi)外交點(diǎn)處，即蝸殼隔板的內(nèi)側(cè)和外側(cè)布置了監(jiān)測點(diǎn)g和h；在蝸殼隔板出口處布置了監(jiān)測點(diǎn)i.

3 數(shù)值模擬結(jié)果與分析

3.1 泵外特性試驗(yàn)及數(shù)值計(jì)算結(jié)果驗(yàn)證

為驗(yàn)證數(shù)值計(jì)算結(jié)果的正確性，進(jìn)行泵的外特性試驗(yàn).試驗(yàn)在江蘇振華泵業(yè)制造有限公司的開式B級(jí)精度試驗(yàn)臺(tái)進(jìn)行.試驗(yàn)過程中，由儀表盤讀出相應(yīng)的流量、揚(yáng)程、效率、功率等數(shù)據(jù)，經(jīng)計(jì)算機(jī)處理后，與數(shù)值計(jì)算的相關(guān)數(shù)據(jù)進(jìn)行比較，如圖6所示.可以看出，泵外特性計(jì)算值與試驗(yàn)值之間誤差在7%內(nèi)，可以認(rèn)為，該數(shù)值計(jì)算方法是可靠的.

圖6 外特性曲線

3.2 泵內(nèi)部壓力脈動(dòng)分析

考慮損管條件等極端工況，艦船用泵的流量區(qū)域都要求比較寬，故本研究的工況點(diǎn)取值為0.4Qd，1.0Qd和1.4Qd，即小流量工況點(diǎn)略大于泵的最小連續(xù)流量點(diǎn)，最大工況點(diǎn)略小于損管條件工作點(diǎn).針對(duì)上述3種工況條件，計(jì)算監(jiān)測點(diǎn)的壓力脈動(dòng)情況.

為便于對(duì)比分析在不同流量工況下的壓力脈動(dòng)情況，將脈動(dòng)幅值轉(zhuǎn)化為同一坐標(biāo)系下進(jìn)行比較，將壓力脈動(dòng)的強(qiáng)度變化定義壓力脈動(dòng)系數(shù)Cp，即

(1)

同時(shí)，為研究各點(diǎn)壓力脈動(dòng)幅值的大小，定義壓力脈動(dòng)幅值CA為

CA=Cpmax-Cpmin，

(2)

式中：Cpmax為葉輪旋轉(zhuǎn)一周內(nèi)脈動(dòng)幅值最大值；Cpmin為葉輪旋轉(zhuǎn)一周內(nèi)脈動(dòng)幅值最小值.

3.2.1 脈動(dòng)幅值分析

圖7為根據(jù)式(2)所繪制的各監(jiān)測點(diǎn)在3種不同流量工況下分別采用直葉片與斜葉片時(shí)泵內(nèi)部壓力脈動(dòng)的幅值大小.

由圖7可以看出：在監(jiān)測點(diǎn)a，d(隔舌處)，泵內(nèi)使用傾斜葉片葉輪時(shí)，壓力脈動(dòng)的幅值較直葉片葉輪明顯降低；對(duì)于蝸殼出口處監(jiān)測點(diǎn)g，h，i，傾斜葉片對(duì)蝸殼內(nèi)的壓力脈動(dòng)影響與直葉片幾乎一致，差別不大，說明葉輪形式對(duì)蝸殼流道的遠(yuǎn)端幾乎無影響；在監(jiān)測點(diǎn)b，e處，斜葉片葉輪較直葉片葉輪在相同流量工況下的脈動(dòng)幅值顯著增加，為后者的3～5倍；泵內(nèi)采用直葉片葉輪時(shí)，壓力脈動(dòng)的最大幅值為0.4Qd工況下的點(diǎn)f處，而采用斜葉片葉輪時(shí)，壓力脈動(dòng)的最大幅值位為1.0Qd工況下的點(diǎn)e處.

由此可見，采用傾斜葉片的葉輪改變了蝸殼腔體內(nèi)葉輪出口附近流體的運(yùn)動(dòng)狀態(tài)，有效降低了隔舌處壓力脈動(dòng)的幅值.雖然采用傾斜葉片的葉輪會(huì)使其他位置的脈動(dòng)幅值最大值增大，以及流量工況點(diǎn)均發(fā)生了改變，但這些變化由于泵腔內(nèi)流體的彈性作用，對(duì)外表現(xiàn)較弱.

3.2.2 壓力脈動(dòng)時(shí)域分析

圖8為泵腔內(nèi)9個(gè)監(jiān)測點(diǎn)在3種流量工況下采用直葉片與斜葉片時(shí)的壓力系數(shù)時(shí)域圖，圖中橫坐標(biāo)φ表示在第8圈內(nèi)，葉輪的轉(zhuǎn)動(dòng)角度.由圖8可以看出：2種葉片形式泵的壓力脈動(dòng)規(guī)律在泵腔內(nèi)各監(jiān)測點(diǎn)處有明顯的周期性，壓力脈動(dòng)的峰值在葉輪每旋轉(zhuǎn)60°出現(xiàn)1次，即葉輪每轉(zhuǎn)1圈其波動(dòng)的周期數(shù)為6，恰好為葉輪的葉片數(shù)；在1個(gè)波動(dòng)周期內(nèi)，存在2～3個(gè)周期數(shù)的小型“駝峰”現(xiàn)象，這種次波動(dòng)的“駝峰”現(xiàn)象在監(jiān)測點(diǎn)a,d(隔舌位置處)表現(xiàn)的最為明顯，當(dāng)采用傾斜葉片葉輪時(shí)，其脈動(dòng)幅度低于直葉片葉輪，其絕對(duì)值小于直葉片葉輪.

圖8 壓力脈動(dòng)時(shí)域分布

綜上所述，采用傾斜葉片出口邊的葉輪無疑優(yōu)于傳統(tǒng)的直葉片布置方式，不僅在關(guān)鍵監(jiān)測點(diǎn)上脈動(dòng)幅度低于直葉片葉輪，同時(shí)抑制了各監(jiān)測點(diǎn)在不同流量工況下脈動(dòng)系數(shù)的小型“駝峰”現(xiàn)象.

3.2.3 壓力脈動(dòng)頻域分析

脈動(dòng)幅值分析時(shí)，在所設(shè)置的監(jiān)測點(diǎn)中，由于蝸殼內(nèi)隔板的作用，壓力脈動(dòng)的波動(dòng)在動(dòng)靜耦合區(qū)域(葉頂間隙區(qū)域)內(nèi)波動(dòng)較為劇烈，遠(yuǎn)離葉輪的區(qū)域相對(duì)平穩(wěn).本節(jié)著重研究在動(dòng)靜耦合區(qū)監(jiān)測點(diǎn)a,b,c,d,e,f的頻域變化情況.

圖9為監(jiān)測點(diǎn)a,b,c,d,e,f在3種工況下分別采用直葉片與斜葉片的壓力脈動(dòng)經(jīng)傅里葉變換得到的頻域圖,圖中fb,fs分別為葉頻、軸頻.葉輪內(nèi)6個(gè)葉片均布，則葉頻為軸頻的6倍.由于監(jiān)測點(diǎn)b,e的壓力脈動(dòng)時(shí)域變化較為特殊，因而在對(duì)b,e進(jìn)行頻域分析時(shí)，采用軸頻作為單位時(shí)間刻度進(jìn)行分析.對(duì)監(jiān)測點(diǎn)a,c,d,f采用葉頻進(jìn)行分析，兩者的縱坐標(biāo)表示脈動(dòng)幅值.

圖9 壓力脈動(dòng)頻域分布

由圖9可以看出：對(duì)于監(jiān)測點(diǎn)b，e，采用直葉片時(shí)，壓力脈動(dòng)能量峰值在各流量點(diǎn)工況下出現(xiàn)在6倍軸頻及其倍頻處(或者說葉頻及其倍頻處)，且在6倍葉頻后，各流量點(diǎn)的脈動(dòng)能量幅值基本趨于0；采用斜葉片時(shí)，壓力脈動(dòng)能量峰值在各流量點(diǎn)工況下相較于使用直葉片時(shí)明顯升高，且峰值出現(xiàn)在2倍軸頻處及其倍頻處，脈動(dòng)能量幅值呈現(xiàn)出周期性降低的態(tài)勢，6倍葉頻后，壓力脈動(dòng)能量幅值的變化規(guī)律與采用直葉片時(shí)的變化規(guī)律逐漸一致，壓力脈動(dòng)幅值能量基本不再發(fā)生改變，能量幅值逐漸趨于0，而監(jiān)測點(diǎn)a,c,d,f的脈動(dòng)能量幅值均出現(xiàn)在葉頻及其倍頻處，由葉頻決定，與使用直葉片相同.

結(jié)合圖8中各監(jiān)測點(diǎn)的時(shí)域脈動(dòng)性能可知，采用斜葉片時(shí)，監(jiān)測點(diǎn)b，e的壓力脈動(dòng)變化規(guī)律與采用直葉片變化規(guī)律產(chǎn)生極大差異，表現(xiàn)為這2點(diǎn)的峰值脈動(dòng)頻率出現(xiàn)在2倍軸頻及其倍頻處.這表明采用斜葉片時(shí)，該點(diǎn)處的射流尾跡較采用直葉片時(shí)產(chǎn)生了相當(dāng)大的改變，考慮到葉輪軸旋轉(zhuǎn)1周，通過上下隔舌各1次，因而這2點(diǎn)的壓力脈動(dòng)變化與隔舌密切相關(guān).

4 結(jié) 論

1) 采用傾斜葉片葉輪較采用直葉片葉輪可以明顯降低泵腔內(nèi)隔舌處的壓力脈動(dòng)幅值，改善壓力脈動(dòng)的次波動(dòng)點(diǎn)的“駝峰”現(xiàn)象，并減小脈動(dòng)峰值的出現(xiàn)頻率，這有益于泵的減振降噪.

2) 同傾斜葉片在單蝸殼泵中所呈現(xiàn)的壓力系數(shù)性能不同，在雙蝸殼泵中采用直葉片葉輪時(shí)，內(nèi)部各點(diǎn)的壓力脈動(dòng)周期數(shù)主要受葉輪葉片數(shù)影響.但當(dāng)采用傾斜葉片時(shí)，會(huì)因?yàn)閺娜~輪甩出的射流尾跡發(fā)生較大的改變，致使在隔舌稍前位置處的監(jiān)測點(diǎn)b,e處的周期數(shù)與脈動(dòng)幅值發(fā)生改變.根據(jù)頻域圖，2點(diǎn)的壓力脈動(dòng)變化主要受隔舌影響.

3)采用傾斜葉片葉輪較采用直葉片葉輪會(huì)使得出現(xiàn)壓力脈動(dòng)的最大幅值的工況流量不同，也會(huì)使得出現(xiàn)壓力脈動(dòng)最大幅值的蝸殼內(nèi)的位置發(fā)生了改變，但這種改變對(duì)泵的振動(dòng)噪聲不構(gòu)成重要影響.