王成斌，金 雷，張 帥，謝 磊，王 勇,3

(1.江蘇大學(xué)國(guó)家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013；2.合肥恒大江海泵業(yè)股份有限公司，合肥 231131；3.江蘇大學(xué)鎮(zhèn)江流體工程裝備技術(shù)研究院，江蘇 鎮(zhèn)江 212009)

0 引 言

徑向?qū)~式潛水泵屬于多級(jí)泵的范疇，具有揚(yáng)程高，體積小等優(yōu)點(diǎn)，被廣泛用于農(nóng)村灌溉、礦井排水等領(lǐng)域。對(duì)于離心泵，尤其是低、中比轉(zhuǎn)速離心泵而言，其軸功率隨著流量的增加而連續(xù)上升[1]，在大流量工況下容易出現(xiàn)過(guò)載，這種現(xiàn)象在徑向?qū)~式潛水泵中更為常見(jiàn)。另一方面，為了響應(yīng)國(guó)家“節(jié)能減排”的號(hào)召，人們對(duì)于泵的效率與穩(wěn)定性等方面越來(lái)越重視，針對(duì)徑向?qū)~式潛水泵的研究具有較高的實(shí)用價(jià)值。

隨著計(jì)算流體力學(xué)的發(fā)展，越來(lái)越多的學(xué)者和研究人員采用數(shù)值計(jì)算方法對(duì)潛水泵進(jìn)行研究[2-5]。數(shù)值計(jì)算可以較準(zhǔn)確地預(yù)測(cè)水力性能，對(duì)完善泵設(shè)計(jì)理論以及提高設(shè)計(jì)水平也有著重要意義。張憶寧[6]基于數(shù)值計(jì)算，對(duì)兩級(jí)潛水泵進(jìn)行了優(yōu)化設(shè)計(jì)，最終選取了葉輪葉片數(shù)Z=3，葉片出口角β2=12°的方案。崔寶玲等[7]基于CFD數(shù)值模擬，對(duì)多級(jí)潛水泵的效率進(jìn)行了優(yōu)化，優(yōu)化后的模型效率提高了4.12%，單級(jí)揚(yáng)程提高1.449 m。程效銳等[8]以5級(jí)井用潛水泵為研究對(duì)象，通過(guò)改變導(dǎo)葉出口邊位置提高了泵的效率和揚(yáng)程。曹衛(wèi)東等[9]結(jié)合正交試驗(yàn)設(shè)計(jì)法與數(shù)值技術(shù)對(duì)導(dǎo)葉不同的兩級(jí)離心泵進(jìn)行了研究，得到了較優(yōu)方案，其揚(yáng)程和效率分別提高1.24m和4%。施衛(wèi)東等[10]采用數(shù)值模擬與性能試驗(yàn)相結(jié)合的方法，驗(yàn)證了等比例縮放法在潛水泵的應(yīng)用。李文全等[11]基于FLUENT對(duì)井筒式潛水軸流泵出水管道進(jìn)行了數(shù)值研究，獲得了進(jìn)水口流量對(duì)出水管道水力特性的影響規(guī)律。Boncinelli等[12]通過(guò)數(shù)值計(jì)算得到潛水泵在不同導(dǎo)葉結(jié)構(gòu)下的性能，從而確定導(dǎo)葉性能較優(yōu)時(shí)的相關(guān)參數(shù)。葛明亞等[13]設(shè)計(jì)了9組不同的潛水模型泵，并與原型泵進(jìn)行了數(shù)值模擬對(duì)比分析，得到了泵性能最優(yōu)時(shí)的葉輪和導(dǎo)葉葉片數(shù)匹配方案。

目前對(duì)徑向?qū)~式潛水泵的內(nèi)流及非定常特性的研究相對(duì)較少。本文以一臺(tái)單級(jí)徑向?qū)~式潛水泵為研究對(duì)象，對(duì)其內(nèi)部流場(chǎng)、葉輪徑向力、軸向力以及壓力脈動(dòng)等特性進(jìn)行分析，以期為潛水泵的優(yōu)化設(shè)計(jì)提供借鑒。

1 計(jì)算模型與計(jì)算設(shè)置

1.1 計(jì)算建模

本文所選用的研究對(duì)象為單級(jí)徑向?qū)~式潛水泵，介質(zhì)為常溫清水，潛水泵結(jié)構(gòu)如圖1所示；主要設(shè)計(jì)參數(shù)為流量Qd=100 m3/h，揚(yáng)程Hd=21 m，轉(zhuǎn)速nd=1 450 r/min；比轉(zhuǎn)速ns=76.6。葉輪的主要參數(shù)如下：進(jìn)口直徑D1=125 mm，出口直徑D2=300 mm，出口寬度b2=13，葉片數(shù)Z=4；另外正、反導(dǎo)葉的葉片數(shù)都為7。

采用Creo2.0軟件對(duì)計(jì)算域水體進(jìn)行三維造型，為了克服邊界條件對(duì)內(nèi)部流場(chǎng)的影響，分別在進(jìn)口處加水槽，在出口處進(jìn)行了延長(zhǎng)。其中水槽為柱體，外徑為葉輪進(jìn)口直徑的8倍，出口段長(zhǎng)度為出口直徑的3倍。圖2為計(jì)算域水體圖，從左到右依次為出口延長(zhǎng)段、反導(dǎo)葉、正導(dǎo)葉、葉輪、進(jìn)水段、水槽。

圖1 泵結(jié)構(gòu)示意圖Fig.1 Pump structure diagram

圖2 計(jì)算域水體Fig.2 Water body of calculation domain

1.2 計(jì)算設(shè)置

1.2.1 網(wǎng)格劃分及網(wǎng)格無(wú)相關(guān)性分析

采用ICEM CFD軟件對(duì)計(jì)算域水體進(jìn)行網(wǎng)格劃分。為了獲得更準(zhǔn)確的計(jì)算結(jié)果，縮短計(jì)算時(shí)間，同時(shí)考慮葉輪的后續(xù)優(yōu)化設(shè)計(jì)，現(xiàn)葉輪水體采用四面體非結(jié)構(gòu)化網(wǎng)格，其余計(jì)算域水體采用六面體結(jié)構(gòu)化網(wǎng)格。同時(shí)為了選取合適的網(wǎng)格數(shù)，需要對(duì)網(wǎng)格進(jìn)行無(wú)關(guān)性檢驗(yàn)。本文通過(guò)調(diào)節(jié)重要計(jì)算面拓?fù)渚€上的網(wǎng)格節(jié)點(diǎn)數(shù)，在CFX前處理設(shè)置完全相同的前提下，得到了5套不同網(wǎng)格數(shù)的網(wǎng)格，并以模擬揚(yáng)程為評(píng)價(jià)標(biāo)準(zhǔn)評(píng)判，如表1所示。

表1 網(wǎng)格無(wú)關(guān)性檢查Tab.1 Inspection of grid independence

由表1中可以看出，隨著網(wǎng)格數(shù)的逐漸增加，揚(yáng)程的計(jì)算值逐漸趨于穩(wěn)定。其中網(wǎng)格數(shù)方案4和方案5中的揚(yáng)程計(jì)算結(jié)果較為接近，這說(shuō)明再增加網(wǎng)格數(shù)對(duì)數(shù)值計(jì)算的結(jié)果影響很小。同時(shí)考慮到計(jì)算機(jī)的運(yùn)算時(shí)間以及其運(yùn)算能力，確定方案4為本文的計(jì)算網(wǎng)格方案。圖3是所確定的各計(jì)算流體域的網(wǎng)格。

圖3 網(wǎng)格劃分Fig.3 Mesh generation

1.2.2 控制方程

在離心泵數(shù)值計(jì)算中所需的控制方程包括質(zhì)量方程、動(dòng)量方程以及能量方程。本文中流體介質(zhì)為清水，是不可壓縮介質(zhì)，只需滿足質(zhì)量方程和動(dòng)量方程[14]。

(1)質(zhì)量守恒方程：

(1)

(2)動(dòng)量守恒方程：

(2)

式中：ρ為流體的密度，kg/m3；u為流體速度，m/s；x為空間坐標(biāo)，i，j為不同方向分量；μ為流體動(dòng)力黏度，Pa.s；Si為源項(xiàng)。

1.2.3 CFX數(shù)值計(jì)算設(shè)置

采用ANSYS CFX軟件進(jìn)行潛水泵的數(shù)值計(jì)算，將葉輪水體設(shè)置為旋轉(zhuǎn)域，其余水體均設(shè)置為靜止域。定常計(jì)算時(shí)，葉輪與其他部件的動(dòng)靜交界面采用凍結(jié)轉(zhuǎn)子模型；非定常計(jì)算時(shí)，動(dòng)靜交界面采用瞬態(tài)動(dòng)靜轉(zhuǎn)子模型。

本文采用對(duì)高應(yīng)變率及流線彎曲程度較大的流動(dòng)有較好自適應(yīng)性的RNGk-ε湍流模型。邊界條件采用壓力進(jìn)口和質(zhì)量流量出口邊界條件組合，進(jìn)口設(shè)置為總壓1 atm，出口處的質(zhì)量流量的數(shù)值根據(jù)流量計(jì)算得到。固壁處設(shè)置為無(wú)滑移壁面，近壁處采用Scalable壁面函數(shù)。

定常計(jì)算時(shí)，控制方程的對(duì)流離散項(xiàng)采用二階高精度，將平均殘差值RMS作為計(jì)算收斂依據(jù)，設(shè)置收斂精度為10-4。非定常數(shù)值計(jì)算以定常計(jì)算結(jié)果作為初始值，時(shí)間步長(zhǎng)設(shè)置為1.149 4×10-4s，即葉輪每旋轉(zhuǎn)1°計(jì)算一次，葉輪共計(jì)算6圈，最后一圈的數(shù)據(jù)用于本文的分析。

2 外特性計(jì)算結(jié)果及試驗(yàn)驗(yàn)證

為了驗(yàn)證數(shù)值計(jì)算的準(zhǔn)確性，在江蘇大學(xué)流體中心實(shí)驗(yàn)室的開(kāi)式試驗(yàn)臺(tái)上進(jìn)行試驗(yàn)，試驗(yàn)臺(tái)如圖4所示。得到了不同工況下外特性的試驗(yàn)值，并與模擬值進(jìn)行對(duì)比，如圖5所示。

由圖5可以看出，數(shù)值計(jì)算得到的性能曲線與試驗(yàn)值變化趨勢(shì)相同。在設(shè)計(jì)流量下，揚(yáng)程的計(jì)算值為20.15 m，與試驗(yàn)值的偏差為2.6%；效率的計(jì)算值為75.42%，與試驗(yàn)值的偏差為1%。除卻大流量工況(流量大于1.2Qd)，效率計(jì)算值在1.2Qd工況下的偏差最大，為4.7%，揚(yáng)程計(jì)算值在0.6Qd工況下的偏差最大，為6.4%。可以看出在流量不超過(guò)1.2Qd時(shí)，對(duì)于潛水泵的計(jì)算是較為準(zhǔn)確的；對(duì)于大流量工況下計(jì)算不準(zhǔn)確在下文中會(huì)進(jìn)行分析。同時(shí)可以看出，該潛水泵的功率曲線在大流量工況下趨于平緩，有較好的無(wú)過(guò)載性能?；谳^為準(zhǔn)確的計(jì)算結(jié)果，選取三個(gè)典型工況(0.8Qd、1.0Qd、1.2Qd)進(jìn)行潛水泵內(nèi)流分布、徑向力、軸向力和壓力分布分析。

圖4 試驗(yàn)裝置圖Fig.4 Diagram of test bench

圖5 計(jì)算與試驗(yàn)結(jié)果對(duì)比Fig.5 Comparison between calculation and experiment

3 計(jì)算結(jié)果分析

3.1 內(nèi)流場(chǎng)分析

選取潛水泵的葉輪出口所在的中間截面為基準(zhǔn)面，對(duì)潛水泵的主要部件葉輪和導(dǎo)葉的靜壓分布、速度分布進(jìn)行分析。選取的中間截面如圖6所示。

圖6 中間截面示意圖Fig.6 Midsection diagram

(1)靜壓分布云圖。圖7分別是3個(gè)工況下，葉輪和導(dǎo)葉在基準(zhǔn)面處的靜壓云圖。從圖7可以看出，在不同的流量下壓力都呈現(xiàn)出從葉輪進(jìn)口到正導(dǎo)葉出口(過(guò)渡段)逐漸增大的規(guī)律。其中，在葉輪進(jìn)口處壓力最低，隨著流量的增加進(jìn)口低壓區(qū)面積逐漸增大。在1.2Qd工況下葉輪流道在稍微偏離進(jìn)口處存在極低壓區(qū)，結(jié)合圖5在流量大于1.2Qd時(shí)揚(yáng)程、效率的計(jì)算值與試驗(yàn)值偏差較大，可以考慮在大流量工況下發(fā)生了汽蝕；正是由于空泡的出現(xiàn)，泵內(nèi)的流動(dòng)變得復(fù)雜，損失加劇，導(dǎo)致?lián)P程下降幅度較大。在0.8Qd工況下，在葉輪出口處及正導(dǎo)葉內(nèi)的流道壓差比較大，液流不穩(wěn)定，隨著流量的增加，該處壓差逐漸減小，液流逐漸穩(wěn)定。在0.8Qd與1.2Qd工況下，正導(dǎo)葉流道與導(dǎo)葉過(guò)渡段流道內(nèi)存在著兩種較為穩(wěn)定的壓力量級(jí)，其中0.8Qd工況下的壓力分別高于1.2Qd工況下的壓力，這與小流量下的揚(yáng)程高于大流量下的揚(yáng)程相對(duì)應(yīng)。整體而言，設(shè)計(jì)工況下葉輪與導(dǎo)葉流道的壓力梯度較小，這也說(shuō)明了額定工況下的潛水泵的內(nèi)流比較穩(wěn)定，其損失較小。

圖7 基準(zhǔn)面靜壓分布云圖Fig.7 Static pressure distribution in midsection

(2)速度分布與流線圖。圖8給出了在基準(zhǔn)面處的葉輪和導(dǎo)葉的速度分布與流線圖，其中葉輪中的速度為相對(duì)速度。從圖中可以看出，不同流量工況下，葉輪和導(dǎo)葉流道內(nèi)均出現(xiàn)不同大小的低速區(qū)。對(duì)于葉輪，低速區(qū)均出現(xiàn)在各個(gè)流道的靠近葉片的工作面處，從進(jìn)口覆蓋到約半個(gè)葉片。葉輪流道的低速區(qū)面積在小流量工況下很明顯，隨著流量的增加，低速區(qū)面積越來(lái)越小，至1.2Qd工況時(shí)已經(jīng)幾乎消失；葉輪流道葉片背面處速度相對(duì)較大，且隨著流量的增加該處的速度逐漸增大。對(duì)于導(dǎo)葉，流道中出現(xiàn)大面積的低速區(qū)，結(jié)合圖7可以看出正導(dǎo)葉起到了動(dòng)能轉(zhuǎn)化為壓力能的作用。在0.8Qd工況下低速區(qū)面積最大，以至于幾乎充滿整個(gè)流道，隨流量增加，低速區(qū)面積逐漸減小，并且低速區(qū)偏向于正導(dǎo)葉流道的右側(cè)葉片處。從流線圖看出，在較低流量下，在葉輪出口靠近正導(dǎo)葉葉片處、正導(dǎo)葉流道與過(guò)渡段內(nèi)均存在漩渦及回流，流動(dòng)紊亂；在設(shè)計(jì)流量下流動(dòng)最穩(wěn)定，除了個(gè)別正導(dǎo)葉有漩渦外，只有在過(guò)渡段流道內(nèi)存在尺寸較小的渦；至較大流量工況下，導(dǎo)葉內(nèi)的渦增多，過(guò)渡段每個(gè)流道內(nèi)依舊存在大小不一的渦。由此可見(jiàn)，過(guò)渡段內(nèi)的流動(dòng)復(fù)雜，損失很大，相對(duì)而言，設(shè)計(jì)工況下的流動(dòng)較穩(wěn)定，與效率最高相對(duì)應(yīng)。

圖8 基準(zhǔn)面速度分布與流線圖Fig.8 Velocity distribution and streamline in midsection

3.2 徑向力與軸向力分析

葉輪是泵的旋轉(zhuǎn)部件，其性能對(duì)泵的穩(wěn)定性有很大的影響，接下來(lái)對(duì)潛水泵葉輪在典型工況0.8Qd、1.0Qd、1.2Qd下的徑向力與軸向力進(jìn)行分析。

(1)葉輪的徑向力。圖9是潛水泵的葉輪在不同工況下的徑向力時(shí)域圖?？梢钥吹皆诓煌r下，葉輪的徑向力都是圍繞一基準(zhǔn)水平線(平均徑向力)上下有規(guī)律的波動(dòng)。三個(gè)工況下，葉輪在一個(gè)旋轉(zhuǎn)周期內(nèi)所受的徑向力都出現(xiàn)28個(gè)波峰與波谷，考慮到葉輪葉片數(shù)為4，導(dǎo)葉葉片數(shù)為7，說(shuō)明葉輪和導(dǎo)葉的動(dòng)靜干涉對(duì)徑向力脈動(dòng)有較大影響。設(shè)計(jì)工況下的平均徑向力最小，0.8Qd工況下次之，1.2Qd工況最大，其分別為34 N、39 N與67.6 N。隨著流量的增加，葉輪所受徑向力的波動(dòng)幅度逐漸增加。

圖9 葉輪徑向力時(shí)域圖Fig.9 Time domain diagram of impeller radial force

圖10是不同工況下葉輪的徑向力分布矢量圖，其中橫、縱坐標(biāo)分別代表葉輪徑向力在中間截面上相互垂直方向的分力?？梢詮膱D中看出，三個(gè)工況下葉輪徑向力的矢量圖都呈現(xiàn)出四角星狀分布，這與葉輪的葉片數(shù)為4相對(duì)應(yīng)。在0.8Qd工況下，矢量圖近似于一個(gè)中空的正方形，說(shuō)明在該工況下徑向力變化比較穩(wěn)定。設(shè)計(jì)工況下的徑向力極值大小與0.8Qd工況相近，1.2Qd工況下徑向力明顯大于其他兩個(gè)工況。

圖10 葉輪徑向力矢量圖Fig.10 Vector diagram of impeller radial force

(2)葉輪的軸向力。圖11給出了葉輪在旋轉(zhuǎn)一圈過(guò)程中不同工況下軸向力的時(shí)域圖。可以看出，不同工況下葉輪軸向力也出現(xiàn)28個(gè)波峰波谷，與徑向力一樣，葉輪、導(dǎo)葉的動(dòng)靜干涉對(duì)葉輪軸向力也有較大影響。相比較徑向力的波動(dòng)而言，葉輪軸向力波動(dòng)振幅較小，比較穩(wěn)定。隨著流量的增加，葉輪的軸向力逐漸增加，軸向力平均值分別為900 N，915 N，931 N，分別大于徑向力；同時(shí)可以看出相鄰工況下的軸向力差值幾乎相等。

圖11 葉輪軸向力時(shí)域圖Fig.11 Time domain diagram of impeller axial force

3.3 壓力脈動(dòng)分析

泵內(nèi)的壓力脈動(dòng)情況對(duì)于泵的運(yùn)行穩(wěn)定性、可靠性有很大的影響，有必要對(duì)其進(jìn)行分析。本文對(duì)潛水泵的主要部件葉輪和導(dǎo)葉上設(shè)置壓力脈動(dòng)的監(jiān)測(cè)點(diǎn)，通過(guò)數(shù)值計(jì)算得到其非定常特性。由于葉輪和導(dǎo)葉流道具有對(duì)稱性，選取單個(gè)流道內(nèi)的相應(yīng)位置設(shè)置監(jiān)測(cè)點(diǎn)。監(jiān)測(cè)點(diǎn)的設(shè)置如圖12所示，其中在葉輪中間截面上的葉片工作面和背面各設(shè)置了4個(gè)監(jiān)測(cè)點(diǎn)，葉輪流道的中間流線上設(shè)置4個(gè)監(jiān)測(cè)點(diǎn)，彼此均勻分布；在正導(dǎo)葉流道、過(guò)渡段流道位于中間截面上分別設(shè)置2個(gè)監(jiān)測(cè)點(diǎn)與1個(gè)監(jiān)測(cè)點(diǎn)；在反導(dǎo)葉流道的中間截面上設(shè)置3個(gè)監(jiān)測(cè)點(diǎn)。

圖12 壓力脈動(dòng)監(jiān)測(cè)點(diǎn)Fig.12 Monitoring points of pressure fluctuation

模擬計(jì)算的結(jié)果得到的壓力值為各個(gè)監(jiān)測(cè)點(diǎn)在不同時(shí)刻的靜壓值，為了更加形象比較壓力大小，引入壓力系數(shù)Cp：

(3)

圖13是0.8Qd工況下各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域圖?？梢钥吹?，各監(jiān)測(cè)點(diǎn)處的壓力脈動(dòng)呈現(xiàn)明顯的周期性。葉輪葉片工作面和背面處的各點(diǎn)處的壓力脈動(dòng)曲線出現(xiàn)7個(gè)主峰，葉輪流道中間流線處的壓力脈動(dòng)除了葉輪出口外，其余主峰數(shù)也為7，與導(dǎo)葉葉片數(shù)7對(duì)應(yīng)，可見(jiàn)導(dǎo)葉對(duì)于葉輪內(nèi)部流道的壓力脈動(dòng)的干涉較大。葉輪流道中間出口處至正導(dǎo)葉中間流道處的脈動(dòng)曲線出現(xiàn)4個(gè)主峰，從過(guò)渡段處至反導(dǎo)葉流道主峰數(shù)變至28個(gè)，可以看出對(duì)于葉輪出口處及正導(dǎo)葉內(nèi)部流場(chǎng)，葉輪對(duì)其脈動(dòng)較大；對(duì)于正導(dǎo)葉出口至反導(dǎo)葉出口的內(nèi)部流道，脈動(dòng)比較復(fù)雜，主要體現(xiàn)出葉輪及導(dǎo)葉在此產(chǎn)生的干涉作用。

圖13 0.8Qd各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域圖Fig.13 Time domain diagram of pressure fluctuation at 0.8Qd

由圖13也可以看出a4處波動(dòng)異常劇烈，這是由于葉片工作面直接對(duì)工作介質(zhì)做功，從前文的分析可知介質(zhì)進(jìn)入導(dǎo)葉后壓力降低，故在出口處涉及能量轉(zhuǎn)化，情況比較復(fù)雜；同時(shí)參考圖7(a)可以看到該處處于幾個(gè)不同等級(jí)的壓力區(qū)的匯集處，所以壓力波動(dòng)比較大。但是如上文所述葉輪在轉(zhuǎn)過(guò)不同流道的時(shí)間內(nèi)也是呈現(xiàn)出來(lái)一定的周期性。對(duì)于葉輪葉片的工作面、背面以及中間流線，從葉輪進(jìn)口至出口壓力脈動(dòng)幅值都是由大到小變化；對(duì)于葉輪的單個(gè)流道的壓力脈動(dòng)，從葉輪葉片的工作面至背面，其幅值沿周向逐漸減小。正導(dǎo)葉內(nèi)的壓力脈動(dòng)幅值沿徑向逐漸減小，至過(guò)渡段及反導(dǎo)葉段脈動(dòng)穩(wěn)定，并且特征趨于一致。

圖14與圖15分別是設(shè)計(jì)工況與1.2Qd工況下的壓力脈動(dòng)時(shí)域圖?？梢钥闯鰞煞N工況下各點(diǎn)的壓力脈動(dòng)的波動(dòng)趨勢(shì)與設(shè)計(jì)工況保持一致，各點(diǎn)處進(jìn)入高壓區(qū)與進(jìn)入低壓區(qū)的時(shí)刻點(diǎn)相同。下面選取1.0Qd工況，也就是設(shè)計(jì)工況，來(lái)作為典型樣本進(jìn)一步分析。

基于快速傅里葉變換(Fast Fourier Transform)將時(shí)域數(shù)據(jù)變換得到頻域數(shù)據(jù)，從而進(jìn)行頻譜分析，如圖16所示。其中定義軸頻倍數(shù)NF：

NF=60F/n=F/Fn

(4)

式中：F為傅里葉變換后的頻率；Fn為相應(yīng)轉(zhuǎn)速下泵軸的轉(zhuǎn)動(dòng)頻率(APF)，本文Fn=24.17 Hz。

圖14 1.0Qd各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域圖Fig.14 Time domain diagram of pressure fluctuation at 1.0Qd

圖15 1.2Qd各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域圖Fig.15 Time domain diagram of pressure fluctuation at 1.2Qd

圖16 1.0Qd各監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)頻域圖Fig.16 Frequency domain diagram of pressure fluctuation at 1.0Qd

由圖16可見(jiàn)，葉輪中無(wú)論葉片的工作面、背面還是中間流線，壓力脈動(dòng)Cp幅值都是在NF為7的時(shí)候達(dá)到最大值，也就是說(shuō)葉輪的壓力脈動(dòng)的主頻為軸頻的7倍，即為7APF(169.17 Hz)。同時(shí)也可以發(fā)現(xiàn)，葉輪內(nèi)除了中間流線出口處的監(jiān)測(cè)點(diǎn)外，其余監(jiān)測(cè)點(diǎn)的次頻都是7APF的倍數(shù)。導(dǎo)葉內(nèi)的主頻由正導(dǎo)葉內(nèi)的4APF逐漸變化到過(guò)渡段與反導(dǎo)葉內(nèi)的7APF。4APF對(duì)應(yīng)葉輪葉片數(shù)，7APF對(duì)應(yīng)導(dǎo)葉葉片數(shù)，由此可見(jiàn)葉輪和導(dǎo)葉內(nèi)的壓力脈動(dòng)主要是葉輪與導(dǎo)葉的動(dòng)靜干涉作用導(dǎo)致。

4 結(jié) 語(yǔ)

(1)在設(shè)計(jì)流量下，揚(yáng)程與效率的計(jì)算值為20.15 m、75.42%，與試驗(yàn)值的偏差分別為2.6%、1%，流量不大于1.2Qd時(shí)數(shù)值計(jì)算較為準(zhǔn)確；在流量大于1.2Qd工況下，潛水泵發(fā)生汽蝕，導(dǎo)致計(jì)算偏差較大。

(2)不同流量下，導(dǎo)葉內(nèi)的過(guò)渡段總存在大小不一的漩渦，損失較大；其中設(shè)計(jì)工況下的液流漩渦規(guī)模最小，相對(duì)最穩(wěn)定。

(3)葉輪的徑向力與軸向力受葉輪與導(dǎo)葉的動(dòng)靜干涉影響，在旋轉(zhuǎn)一圈的時(shí)間內(nèi)均出現(xiàn)28個(gè)波峰波谷。其中軸向力的波動(dòng)振幅小于徑向力，但是其平均值大于徑向力。

(4)葉輪的壓力脈動(dòng)受葉輪與導(dǎo)葉的干涉，其主頻為軸頻的7倍(7APF)，導(dǎo)葉內(nèi)的壓力脈動(dòng)的主頻由正導(dǎo)葉內(nèi)的4APF逐漸變化到過(guò)渡段與反導(dǎo)葉內(nèi)的7APF。

□