張德勝 楊 港 趙旭濤 楊雪琪 高雄發(fā)

(江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心，鎮(zhèn)江 212013)

0 引言

近年來南水北調(diào)、地區(qū)灌溉等長距離調(diào)水和大型灌溉工程緩解了水資源分布不均和農(nóng)業(yè)灌溉用水問題。立式離心泵憑其流量大、揚程高的優(yōu)勢在該領(lǐng)域發(fā)揮著重要作用。流量大、揚程高的特點導(dǎo)致立式離心泵裝機功率大，長時間運行能耗極高，提升其水力性能是減少該類泵能耗的重要手段。然而，為了滿足更多的運行要求和更好的結(jié)構(gòu)要求，該類泵通常帶有徑向?qū)~結(jié)構(gòu)[1-2]；導(dǎo)葉與葉輪和蝸殼的匹配性較低是其產(chǎn)生較大水力損失的重要因素之一，因此提升各過流部件的匹配性是改善其性能的有效手段[3-4]。

近年來，近似模型被廣泛運用于流體機械的優(yōu)化設(shè)計中[5-7]，其對非線性關(guān)系的高精度擬合能力成為提高優(yōu)化設(shè)計精度和效率的關(guān)鍵，人工神經(jīng)網(wǎng)絡(luò)模型在眾多近似模型中表現(xiàn)出更好的對泵性能參數(shù)的預(yù)測能力[8-10]。利用近似模型結(jié)合智能優(yōu)化算法的優(yōu)化設(shè)計方法已成為流體機械領(lǐng)域快速、有效的優(yōu)化設(shè)計方法。文獻[11]對蝸殼的4個設(shè)計變量做了多目標優(yōu)化設(shè)計，采用均勻試驗設(shè)計、流固耦合計算、BP神經(jīng)網(wǎng)絡(luò)和多目標遺傳算法組合的方法提高了雙流道泵的水力性能和結(jié)構(gòu)性能；文獻[12]利用貝塞爾曲線擬合管道泵肘形進水管的形狀，采用人工神經(jīng)網(wǎng)絡(luò)、拉丁方試驗設(shè)計方法和多目標遺傳算法，以3個工況下的效率最優(yōu)為目標優(yōu)化了肘形進水管的形狀；文獻[13]搭建了軸流泵自動優(yōu)化平臺，采用徑向基神經(jīng)網(wǎng)絡(luò)與多島遺傳算法相結(jié)合，將軸流泵的水力效率提升了4.46%；文獻[14]選取了混流泵導(dǎo)葉的4個設(shè)計變量進行優(yōu)化設(shè)計，采用拉丁超立方采樣、徑向基神經(jīng)網(wǎng)絡(luò)和序列二次規(guī)劃算法結(jié)合的方法將混流泵設(shè)計工況下的水力效率提升了9.75%；文獻[15]以設(shè)計工況下效率為目標，采用徑向基神經(jīng)網(wǎng)絡(luò)和遺傳算法相結(jié)合，對一余熱排出泵的葉輪軸面投影圖進行了優(yōu)化設(shè)計，優(yōu)化后的葉輪性能改善明顯；文獻[16]以泵工況和水輪機工況下的效率為目標，采用基于帕累托的遺傳算法對水泵水輪機葉輪進行優(yōu)化，優(yōu)化后的模型在泵和水輪機工況下的效率分別提高了0.27%和16.3%。

上述研究成果可見，針對靜止過流部件及旋轉(zhuǎn)過流部件的設(shè)計，近似模型與優(yōu)化算法相結(jié)合的方法都在流體機械優(yōu)化設(shè)計領(lǐng)域得到了成功的應(yīng)用，但多數(shù)研究都是對單個過流部件進行優(yōu)化設(shè)計，從多個過流部件的匹配性考慮對立式離心泵進行優(yōu)化設(shè)計的研究相對缺乏。因此，本文針對立式離心泵各過流部件間的匹配性，采用BP神經(jīng)網(wǎng)絡(luò)和多島遺傳算法相結(jié)合的方法，重點提高其設(shè)計工況下的水力效率，以降低泵的運行能耗，達到節(jié)能運行的目標。

1 計算模型與數(shù)值模擬方法

1.1 立式離心泵水力模型及參數(shù)

本文的立式離心泵計算模型由原型泵水力模型通過相似換算得到。過流部件由肘形進水管、離心式葉輪、活動導(dǎo)葉、固定導(dǎo)葉以及螺旋形蝸殼組合而成。其基本設(shè)計參數(shù)：設(shè)計流量Qdes為0.214 m3/s，揚程H為18 m，額定轉(zhuǎn)速n為1 150 r/min，比轉(zhuǎn)數(shù)ns為215；其主要結(jié)構(gòu)參數(shù)為：葉輪出口直徑D2為360 mm，葉輪出口寬度b2為70.3 mm，葉輪葉片數(shù)Z為7，固定導(dǎo)葉葉片數(shù)Zs為20，活動導(dǎo)葉葉片數(shù)Zg為20。

采用三維造型軟件構(gòu)建立式離心泵模型機水力模型得到本文的原始計算模型，如圖1所示?；顒訉?dǎo)葉和固定導(dǎo)葉位于葉輪與蝸殼之間的重要位置，是連接葉輪與蝸殼的重要過流部件，其主要的作用是控制流量及調(diào)節(jié)最優(yōu)工況點。

1.2 網(wǎng)格劃分

對進水管、葉輪、活動導(dǎo)葉、固定導(dǎo)葉和蝸殼5個過流部件進行了六面體結(jié)構(gòu)網(wǎng)格劃分，對近壁面處的網(wǎng)格作精細化處理，其y+值(無量綱值，表示距離壁面第1層網(wǎng)格的高度)小于50，滿足湍流模型的要求。劃分后的網(wǎng)格如圖2所示；通過設(shè)置7組不同數(shù)量的流體計算域網(wǎng)格來進行網(wǎng)格無關(guān)性驗證，其結(jié)果如圖3所示。

圖2 立式離心泵計算域網(wǎng)格Fig.2 Grid of vertical centrifugal pump1.蝸殼 2.固定導(dǎo)葉 3.活動導(dǎo)葉 4.葉輪

圖3 網(wǎng)格無關(guān)性驗證Fig.3 Mesh independence

由圖3可知，當網(wǎng)格數(shù)大于700萬時，隨著網(wǎng)格數(shù)的增加，對應(yīng)的揚程和效率浮動趨于穩(wěn)定，即網(wǎng)格數(shù)對最終的性能指標已幾乎不產(chǎn)生影響。考慮到計算時長和計算精度，最終將總網(wǎng)格數(shù)定為700萬左右，其中進水管網(wǎng)格數(shù)約為103萬，葉輪網(wǎng)格數(shù)約為174萬，活動導(dǎo)葉網(wǎng)格數(shù)約為155萬，固定導(dǎo)葉網(wǎng)格數(shù)約為82萬，蝸殼網(wǎng)格數(shù)約為186萬。

1.3 數(shù)值模擬方法及設(shè)置

在商業(yè)軟件ANSYS CFX中求解穩(wěn)定的三維RANS方程來預(yù)測立式離心泵的性能參數(shù)。SST(剪切應(yīng)力傳輸)k-ω湍流模型在標準的k-ω湍流模型基礎(chǔ)上改進了渦黏性的定義來考慮湍流主切應(yīng)力運輸?shù)挠绊懀瑥亩倪M了對逆壓梯度流動的預(yù)測[17-18]。在大型水泵數(shù)值模擬過程中采用此湍流模型來封閉N-S方程，可以較為準確地預(yù)測其內(nèi)部流動特性，因此本文選用此湍流模型來計算立式離心泵中的流動特性。

本文的邊界條件設(shè)置為：入口邊界條件設(shè)置總壓為1個標準大氣壓，出口邊界條件設(shè)置為質(zhì)量流量，根據(jù)工況的不同而改變；設(shè)置葉輪的轉(zhuǎn)速為1 150 r/min，旋轉(zhuǎn)域和靜止域間的交界面設(shè)置為凍結(jié)轉(zhuǎn)子；固體壁面滿足無滑移條件設(shè)置為無滑移壁面；葉輪與進水管和徑向?qū)~的相交面設(shè)置為動-靜交界面，其余各部件交界面均為靜-靜交界面；收斂精度定義為10-4。

1.4 試驗驗證

通過對原始模型泵進行外特性試驗，來驗證本文數(shù)值模擬方法的準確性，原始模型的外特性試驗在水利部天津中水北方勘測設(shè)計研究院的立式水泵閉式試驗臺上進行，圖4為試驗臺的運行系統(tǒng)，圖5為葉輪與立式離心泵閉式試驗臺的實物圖。

圖4 試驗裝置運行系統(tǒng)Fig.4 System diagram of test devices1.高壓罐 2.空氣閥 3.蓄水池 4、5.閥門 6.真空泵 7.低壓罐 8.壓力表 9.試驗泵

試驗所得泵外特性曲線與數(shù)值模擬結(jié)果對比如圖6所示。從圖6可以看出，CFD(計算流體力學(xué))模擬結(jié)果與試驗所得外特性曲線趨勢較為一致，性能參數(shù)的計算值與試驗值最大偏差小于5%，在可接受誤差范圍內(nèi)，說明本文采用的數(shù)值模擬方法較為準確。

圖6 外特性模擬曲線與試驗曲線Fig.6 Performance curves of simulation and experiment

圖6中對揚程和流量進行無量綱化，得到相應(yīng)的揚程系數(shù)和流量系數(shù)，其中本文設(shè)計工況對應(yīng)的流量系數(shù)為0.124，設(shè)計揚程對應(yīng)的揚程系數(shù)為0.751，流量系數(shù)和揚程系數(shù)定義為

(1)

(2)

式中g(shù)——重力加速度，m/s2

Q——泵的流量，m3/s

u2——葉輪出口圓周速度，m/s

2 優(yōu)化設(shè)計方法與應(yīng)用

本文旨在改善過流部件匹配性以提高立式離心泵水力性能，優(yōu)化設(shè)計過程如圖7所示，以設(shè)計工況下的水力效率為優(yōu)化目標，首先對原始模型導(dǎo)葉結(jié)構(gòu)進行了改進，對改進后的模型選取導(dǎo)葉與蝸殼各設(shè)計參數(shù)進行Plackett-Burman篩選試驗，以減少優(yōu)化設(shè)計變量；然后采用最優(yōu)拉丁超立方采樣在設(shè)計范圍內(nèi)采樣106組樣本，并搭建了自動數(shù)值模擬優(yōu)化平臺以快速得到樣本點相應(yīng)的優(yōu)化目標值；通過BP神經(jīng)網(wǎng)絡(luò)模型來構(gòu)建優(yōu)化設(shè)計變量和優(yōu)化目標之間的非線性關(guān)系；最終通過多島遺傳優(yōu)化算法求解非線性函數(shù)得到導(dǎo)葉與蝸殼的最優(yōu)參數(shù)組合。

圖7 立式離心泵優(yōu)化流程圖Fig.7 Optimization flow chart of vertical centrifugal pump

2.1 優(yōu)化目標

為盡可能多地降低能耗，本文的優(yōu)化設(shè)計目標選取Qdes工況下的模型泵水力效率，通過定常數(shù)值模擬計算得到，即

(3)

式中T——葉輪扭矩，N·m

ω——葉輪角速度，rad/s

ρ——密度，kg/m3

2.2 導(dǎo)葉結(jié)構(gòu)改進

立式離心泵中的徑向?qū)~結(jié)構(gòu)又稱為葉片式擴壓器，是將葉輪內(nèi)流出的高速水流過渡到蝸殼中的關(guān)鍵水力部件，本文將原始模型的活動導(dǎo)葉和固定導(dǎo)葉簡化為一層固定導(dǎo)葉，固定導(dǎo)葉依然可以起到引導(dǎo)水流的作用，這不僅簡化了模型結(jié)構(gòu)，還能夠降低水流在兩層導(dǎo)葉中的損失從而提高泵工況運行的效率和穩(wěn)定性。同時考慮到動靜干涉作用，導(dǎo)葉數(shù)和葉輪葉片數(shù)需互為質(zhì)數(shù)，因此將固定導(dǎo)葉的數(shù)目從原來的20片減少到13片來增大導(dǎo)葉間水流的過流面積，上述改進后的模型較原始模型在設(shè)計工況下效率提升1個百分點，達到87.6%；揚程提升0.57 m，達到18.57 m。

改進后的模型較原始模型揚程和效率都有所提高，模型性能得到了明顯改進；改進后模型導(dǎo)葉與蝸殼的基本設(shè)計參數(shù)：導(dǎo)葉進口直徑D3為379.8 mm，導(dǎo)葉出口直徑D4為434 mm，導(dǎo)葉進出口寬度b3保持一致均為75.6 mm，導(dǎo)葉進口安放角β3為27°，固定導(dǎo)葉葉片數(shù)Zs減少為13片，蝸殼基圓直徑D5為440 mm，蝸殼進口寬度b4為118 mm，蝸殼出口直徑D6為330 mm。

此外，改進后模型效率只提升1個百分點左右，尚有較大提升空間，由于改進過程只是將活動導(dǎo)葉結(jié)構(gòu)去除，并未考慮到導(dǎo)葉與葉輪和蝸殼的匹配性，但導(dǎo)葉結(jié)構(gòu)處于特殊位置，其與葉輪和蝸殼的匹配性是影響立式離心泵性能的關(guān)鍵因素之一，因此本文針對過流部件匹配性進行下一階段的優(yōu)化設(shè)計。

2.3 設(shè)計參數(shù)優(yōu)化

導(dǎo)葉的設(shè)計參數(shù)極大地影響其與葉輪和蝸殼的匹配程度，因此本文選取導(dǎo)葉主要的6個設(shè)計參數(shù)作為初始優(yōu)化變量，分別為：D3、D4、導(dǎo)葉包角φ、β3、導(dǎo)葉出口安放角β4、b3；具體設(shè)計參數(shù)如圖8所示。

圖8 導(dǎo)葉主要優(yōu)化參數(shù)定義Fig.8 Definition of main geometry variables of vane diffuser

與此同時，本文選取蝸殼中的4個設(shè)計參數(shù)作為優(yōu)化設(shè)計變量來進一步考慮導(dǎo)葉與蝸殼的匹配性，分別為：S8、D6、蝸殼擴散管高度L、b4；具體設(shè)計參數(shù)如圖9所示。其中第8斷面又稱蝸殼喉部，其面積是蝸殼設(shè)計中的核心參數(shù)，本文在改變蝸殼喉部面積時，其他各斷面面積也隨之改變，具體定義方法為

圖9 蝸殼主要優(yōu)化參數(shù)定義Fig.9 Definition of main geometry variables of volute

(4)

式中Sφ——計算斷面面積，mm2

φ8——計算斷面與第8斷面間夾角，(°)

S8——蝸殼喉部面積，mm2

2.4 Plackett-Burman試驗設(shè)計

Plackett-Burman試驗設(shè)計[19]能夠從多個因素中篩選出對試驗指標有顯著影響的因素，該試驗研究(N-1)個因素僅需要進行N次試驗，其中N為4的倍數(shù)，可以以較少的試驗次數(shù)篩選出對試驗指標具有顯著影響的因素。對立式離心泵導(dǎo)葉和蝸殼的10個設(shè)計參數(shù)選取高低兩組設(shè)計水平進行Plackett-Burman試驗設(shè)計，表1給出了試驗設(shè)計因素及水平。

表1 Plackett-Burman設(shè)計因素及水平Tab.1 Range of different factors investigated with Plackett-Burman

Plackett-Burman試驗設(shè)計共進行12次，采用CFX軟件數(shù)值計算得到各方案對應(yīng)的水力效率作為Plackett-Burman試驗設(shè)計的試驗指標。根據(jù)得出的12組計算結(jié)果，建立了以效率為評價指標的回歸方程。表2顯示了設(shè)計點效率和10個設(shè)計參數(shù)之間的皮爾遜相關(guān)系數(shù)和顯著性檢驗單尾概率P值，皮爾遜相關(guān)系數(shù)的絕對值越大則相關(guān)性越大，而P值在0.05以下為顯著影響，在0.05和0.1之間為次顯著影響。

表2 效率顯著性分析結(jié)果Tab.2 Significance analysis results of efficiency

從表2中可以看出，導(dǎo)葉進口直徑D3的相關(guān)系數(shù)達到0.615，顯著性P值只有0.017，說明D3對效率影響最大，為顯著性影響因素，且為正影響，其次蝸殼喉部面積對效率的相關(guān)系數(shù)為-0.434，且顯著性P值為0.079，是對水力效率影響第二大的因素，為次顯著性影響因素，且為負影響；同理還可以發(fā)現(xiàn)，蝸殼進口寬度和導(dǎo)葉進口安放角對效率影響也較大。因此，本文選取上述4個設(shè)計參數(shù)作為最終的優(yōu)化設(shè)計變量。

2.5 最優(yōu)拉丁超立方采樣(OLHS)

為了構(gòu)建高精度的近似模型，需在樣本設(shè)計空間內(nèi)選取合適的樣本點，最優(yōu)拉丁超立方采樣是一種分層采樣方法[20]，能夠以較少的采樣次數(shù)充分反映設(shè)計空間樣本特性；對4個優(yōu)化設(shè)計變量在設(shè)計空間內(nèi)采樣106組樣本點，樣本點數(shù)量遠大于構(gòu)建近似模型所需的樣本數(shù)量，采樣空間范圍見表3。

表3 采樣空間范圍Tab.3 Range of sampling space

2.6 自動數(shù)值模擬優(yōu)化平臺

在Isight軟件中通過批處理命令將CFturbo、ICEM CFD和CFX軟件進行集成，搭建了立式離心泵從設(shè)計到仿真一體化的自動數(shù)值模擬優(yōu)化平臺，以保證數(shù)值模擬的準確性和一致性，同時縮短了優(yōu)化周期，軟件集成運行框架如圖10所示。

圖10 軟件集成運行框架Fig.10 Software integration framework

通過此自動優(yōu)化平臺得到106組方案對應(yīng)的優(yōu)化目標值。采用旋轉(zhuǎn)機械設(shè)計軟件CFturbo對立式離心泵導(dǎo)葉和蝸殼實現(xiàn)參數(shù)化設(shè)計；將導(dǎo)葉與蝸殼的三維模型分別導(dǎo)入ICEM CFD，啟用腳本程序自動進行網(wǎng)格劃分，最終再一起導(dǎo)入CFX軟件自動完成數(shù)值模擬設(shè)置和計算。

2.7 BP神經(jīng)網(wǎng)絡(luò)模型

構(gòu)建優(yōu)化設(shè)計變量與優(yōu)化目標之間關(guān)系是優(yōu)化設(shè)計過程中的關(guān)鍵步驟，也是關(guān)系到最終優(yōu)化結(jié)果準確性的重要因素之一，人工神經(jīng)網(wǎng)絡(luò)模型(Artificial neural network,ANN)是實現(xiàn)上述目標的有效手段。人工神經(jīng)網(wǎng)絡(luò)按照拓撲結(jié)構(gòu)的不同又分為前饋神經(jīng)網(wǎng)絡(luò)和反饋神經(jīng)網(wǎng)絡(luò)，其中BP神經(jīng)網(wǎng)絡(luò)是前饋神經(jīng)網(wǎng)絡(luò)中的常見類型，近年來在眾多工程領(lǐng)域中被廣泛運用[21-23]。

BP神經(jīng)網(wǎng)絡(luò)具有很強的擬合非線性函數(shù)關(guān)系的能力，其分為輸入層、隱藏層和輸出層，通過逐層映射來構(gòu)建變量與目標值之間的關(guān)系，其網(wǎng)絡(luò)結(jié)構(gòu)如圖11所示。

圖11 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.11 Structure of back propagation neural network model

運用上述得到的106組樣本數(shù)據(jù)，采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建優(yōu)化設(shè)計變量和優(yōu)化目標之間的非線性關(guān)系。因此輸入層中含有4個神經(jīng)元，輸出層中含有1個神經(jīng)元且隱藏層設(shè)置10個神經(jīng)元。將85%的數(shù)據(jù)用于訓(xùn)練神經(jīng)網(wǎng)絡(luò)，15%的數(shù)據(jù)用來檢驗擬合精度，擬合精度由決定系數(shù)R2決定，從圖12中可以看出，本文構(gòu)建的神經(jīng)網(wǎng)絡(luò)模型R2達到0.988 1，擬合精度較高。

圖12 R2誤差分析Fig.12 R2 error analysis

2.8 多島遺傳算法(MIGA)

遺傳算法(Genetic algorithm,GA)是一種經(jīng)典的全局優(yōu)化算法，在算法的運行過程中將優(yōu)化問題的解當作染色體，采用選擇、交叉、變異等方法不斷地迭代出新的種群從而得到最優(yōu)解，這與生物進化過程中的遺傳機制十分類似[24]。

多島遺傳算法則是對傳統(tǒng)遺傳算法進行改進的一種全局優(yōu)化算法[25]，如圖13所示，將一個種群分解為多個小的子種群稱為“島”，子種群依然使用遺傳算法進行求解，但是各個子種群中的個體在不斷地“遷移”使得個體更具多樣性；多島遺傳算法與傳統(tǒng)的遺傳算法相比收斂更快且全局尋優(yōu)能力更強。

圖13 多島遺傳算法Fig.13 Multi-island genetic algorithm

本文運用該算法結(jié)合上文中訓(xùn)練完成的BP神經(jīng)網(wǎng)絡(luò)對立式離心泵的優(yōu)化問題進行求解，算法設(shè)置為：子種群個數(shù)為12，島的個數(shù)為12，交叉率為0.9，島間遷移率和變異概率為0.01，遷移間隔代數(shù)為5代。

3 結(jié)果與分析

3.1 外特性分析

在規(guī)定步數(shù)的迭代計算后，多島遺傳算法給出了優(yōu)化后的變量值，采用優(yōu)化后的變量對導(dǎo)葉與蝸殼重新進行三維造型和數(shù)值計算，得到優(yōu)化后模型設(shè)計工況下的水力性能參數(shù)；優(yōu)化前后的設(shè)計參數(shù)及性能參數(shù)對比如表4所示。從表4可以看出，優(yōu)化后導(dǎo)葉進口直徑增大，即葉輪與導(dǎo)葉間的無葉區(qū)范圍擴大，導(dǎo)葉進口安放角與蝸殼進口寬度也有所增大，而蝸殼喉部面積明顯減小，即其他斷面過流面積也隨之減小；優(yōu)化后模型設(shè)計點效率達到90.21%，較原始模型提高3.61個百分點，較改進后模型提高2.61個百分點。同時，BP神經(jīng)網(wǎng)絡(luò)給出的最優(yōu)效率預(yù)測值為90.53%，與數(shù)值模擬實際值相差0.32個百分點，進一步證明了神經(jīng)網(wǎng)絡(luò)模型的高精度。

表4 優(yōu)化前后參數(shù)對比Tab.4 Comparison of parameters between original and optimized model

為了進一步分析優(yōu)化后模型多工況下性能參數(shù)的變化情況，優(yōu)化后模型、改進后模型和原始模型的外特性曲線如圖14所示。從圖14可以看出，優(yōu)化后模型多工況下的效率和揚程都明顯高于原始模型，在設(shè)計點水力效率得到明顯提升。改進后的模型主要在設(shè)計工況和大流量工況下效率有所提升，但在小流量工況下效率提升不明顯，且最優(yōu)效率點偏向大流量工況；而優(yōu)化后模型相比于改進后模型，在設(shè)計流量和小流量下的效率和揚程有明顯的提升，雖然優(yōu)化后模型在大流量下效率和揚程有所下降，但相比于原始模型依然有顯著提升，并且解決了最優(yōu)效率點向大流量工況偏移的問題。

圖14 優(yōu)化前后模型外特性對比Fig.14 Comparison of pump performance before and after optimization

改進后模型在大流量下的性能提高，可能是由于導(dǎo)葉層數(shù)的減少和導(dǎo)葉葉片數(shù)減少導(dǎo)致導(dǎo)葉內(nèi)過流面積的增大使其更加適應(yīng)大流量工況，而優(yōu)化后模型在設(shè)計流量和小流量工況下的性能提升，可能是由于優(yōu)化后模型導(dǎo)葉與葉輪和蝸殼的匹配性提高；優(yōu)化后的導(dǎo)葉與蝸殼對立式離心泵設(shè)計流量和小流量工況下的性能影響更加顯著。

3.2 內(nèi)部流場分析

為了探究過流部件的匹配性對立式離心泵內(nèi)部流場的影響，對改進后模型與優(yōu)化后模型的內(nèi)部流場進行對比分析；圖15給出了導(dǎo)葉不同高度位置的示意圖，其中Span值0.1為靠近前蓋板位置，Span值0.9為靠近后蓋板位置；改進后模型與優(yōu)化后模型在設(shè)計流量下導(dǎo)葉各高度位置的速度流線分布如圖16所示。

圖15 導(dǎo)葉不同高度位置示意圖Fig.15 Schematic at different spans of diffuser

從圖16可以看出，改進后模型導(dǎo)葉內(nèi)的流動較為不穩(wěn)定，多個導(dǎo)葉流道內(nèi)存在大面積低速區(qū)，且不同導(dǎo)葉高度位置流動分布不均勻，Span值0.5位置部分流道內(nèi)漩渦的產(chǎn)生導(dǎo)致流道堵塞，迫使流體向相鄰流道或?qū)~其他高度位置流動，從而增大了導(dǎo)葉內(nèi)的水力損失。優(yōu)化后模型導(dǎo)葉內(nèi)的速度和流線分布在各流道和不同高度位置分布更加均勻且大面積低速區(qū)和局部漩渦消失，僅在部分導(dǎo)葉后緣存在輕微的流動分離現(xiàn)象；這是因為優(yōu)化后的導(dǎo)葉進口直徑與進口安放角都較改進后模型增大，導(dǎo)葉與葉輪的匹配性提高，從而改善了導(dǎo)葉內(nèi)的流態(tài)。

圖16 導(dǎo)葉不同高度位置速度流線分布Fig.16 Velocity and streamline distributions of vane diffuser at different spans

采用阻塞系數(shù)來描述優(yōu)化前后導(dǎo)葉各流道內(nèi)的阻塞情況，導(dǎo)葉流道示意圖如圖17所示，阻塞系數(shù)小于1時，流道內(nèi)存在阻塞情況，阻塞系數(shù)越小說明流道內(nèi)的阻塞情況越嚴重，定義阻塞系數(shù)小于0.95為嚴重堵塞情況，圖18給出了改進后模型和優(yōu)化后模型在設(shè)計工況下導(dǎo)葉內(nèi)各流道的阻塞系數(shù)對比，其中阻塞系數(shù)定義為

圖17 導(dǎo)葉流道示意圖Fig.17 Definition of vane diffuser channel

圖18 導(dǎo)葉內(nèi)阻塞系數(shù)對比Fig.18 Comparison of blocking coefficients in vane diffuser

(5)

式中BD——導(dǎo)葉內(nèi)阻塞系數(shù)

QD——導(dǎo)葉單個流道質(zhì)量流量，kg/s

Qm——設(shè)計工況質(zhì)量流量，kg/s

從圖18可以看出，優(yōu)化前后導(dǎo)葉內(nèi)的阻塞系數(shù)存在較大變化，改進后模型相鄰導(dǎo)葉流道阻塞系數(shù)差值較大，多個流道存在嚴重堵塞情況，這與圖16a中的速度流線分布相對應(yīng)，流道內(nèi)漩渦的產(chǎn)生使流道堵塞情況嚴重，過流能力嚴重降低，并迫使流體流向相鄰流道，導(dǎo)致相鄰流道阻塞系數(shù)明顯升高；優(yōu)化后的模型導(dǎo)葉流道內(nèi)無明顯嚴重堵塞情況，原本阻塞嚴重的導(dǎo)葉流道阻塞系數(shù)都得到顯著提升，這是由于圖16b中優(yōu)化后導(dǎo)葉內(nèi)的速度流線分布更加均勻；綜合來看，優(yōu)化后導(dǎo)葉內(nèi)的阻塞情況有所改善，導(dǎo)葉與葉輪的匹配性提高，使得導(dǎo)葉內(nèi)的流動更加穩(wěn)定，水力損失減小，從而提高了立式離心泵的水力效率。

為了進一步探究導(dǎo)葉與蝸殼的匹配性對立式離心泵水力性能的影響，給出了改進后模型和優(yōu)化后模型設(shè)計工況下導(dǎo)葉與蝸殼內(nèi)的速度流線分布。如圖19a所示，改進后模型蝸殼內(nèi)速度分布不均且流速較小，蝸殼進口與出口位置速度差不明顯，蝸殼內(nèi)各斷面與出口管處的流線分布也較為不均，蝸殼斷面面積較小側(cè)流線反而更加密集，這是因為部分導(dǎo)葉流道存在低速區(qū)產(chǎn)生漩渦阻塞流道，而相鄰導(dǎo)葉流道為高速流動，導(dǎo)致各導(dǎo)葉流道的出流速度和出流量都存在較大差異，且導(dǎo)葉與蝸殼匹配性較低，導(dǎo)葉內(nèi)的出流方向與蝸殼螺旋形狀不符，從而影響到蝸殼的內(nèi)部流動產(chǎn)生較大水力損失；如圖19b所示，優(yōu)化后模型蝸殼斷面面積減小，蝸殼進口寬度增大，蝸殼內(nèi)速度分布較為均勻且蝸殼進口和出口位置速度差明顯，說明優(yōu)化后的蝸殼能夠更好地將速度能轉(zhuǎn)換為壓力能，同時蝸殼內(nèi)的流線分布與各斷面面積變化規(guī)律較為一致且導(dǎo)葉出口處的流線方向與蝸殼螺旋形狀更加貼合，蝸殼出口處的流線分布也得到了極大改善；上述現(xiàn)象說明優(yōu)化后導(dǎo)葉與蝸殼的匹配性更好，導(dǎo)葉內(nèi)的流動改善，使得蝸殼內(nèi)的流態(tài)也隨之改善，因此導(dǎo)葉與蝸殼的匹配性提升也是立式離心泵性能改善的重要原因之一。

圖19 導(dǎo)葉與蝸殼速度流線分布Fig.19 Velocity and streamline distributions of vane diffuser and volute

4 結(jié)論

(1)進行Plackett-Burman篩選試驗后發(fā)現(xiàn)，導(dǎo)葉進口直徑D3和蝸殼喉部面積S8對立式離心泵水力性能影響較大。

(2)訓(xùn)練的BP神經(jīng)網(wǎng)絡(luò)R2達到0.988 1，且優(yōu)化后的網(wǎng)絡(luò)預(yù)測效率與數(shù)值模擬計算效率僅相差0.32個百分點，說明BP神經(jīng)網(wǎng)絡(luò)能夠較為準確地映射出泵設(shè)計參數(shù)與性能參數(shù)間的內(nèi)在聯(lián)系。

(3)通過外特性分析發(fā)現(xiàn)，優(yōu)化后模型多工況下水力效率和揚程都有所提高，設(shè)計工況下效率達到90.21%，較原始模型提高3.61個百分點，較改進后模型提高2.61個百分點，優(yōu)化后的導(dǎo)葉與蝸殼對設(shè)計流量和小流量工況下的泵性能影響更加明顯。

(4)通過內(nèi)部流場分析發(fā)現(xiàn)，優(yōu)化后模型導(dǎo)葉與葉輪和蝸殼的匹配性顯著提高，優(yōu)化后導(dǎo)葉流道內(nèi)的阻塞情況有所改善，蝸殼進出口位置速度差更大，蝸殼斷面面積較小側(cè)與出口處流線分布更加均勻。