基于泵軸流固耦合的雙向流道軸流泵裝置的數(shù)值分析

陳 佳，裴 吉，袁壽其，李彥軍，孟 凡(江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心，江蘇 鎮(zhèn)江 212013)

0 引 言

泵站作為水利事業(yè)的重要組成部分，在人類改造大自然中發(fā)揮了重要作用，也在建設(shè)高產(chǎn)穩(wěn)產(chǎn)農(nóng)田、跨流域調(diào)水、城鎮(zhèn)供排水、工業(yè)供水、礦山排水、流體輸送等方面發(fā)揮著越來越重要的作用。雙向立式軸流泵站同時具有排澇和灌溉的作用，因此被廣泛應(yīng)用在沿江濱湖地區(qū)，泵站運行中的穩(wěn)定性問題也越來越受到重視。由于雙向流道泵站單向運行時，在進(jìn)出水流道中易形成漩渦，造成較大的壓力波動，而泵站出水流道中連接電機的泵軸直接暴露在流道中，經(jīng)常引起軸的振動和偏向運行，進(jìn)而導(dǎo)致機組運行不穩(wěn)定性增強。因此研究泵裝置出水流道中軸的振動特性具有重要的理論意義和工程價值。

在出水流道中泵軸受流場壓力脈動作用較大，且隨時間不斷變化，因此不考慮流體與結(jié)構(gòu)的相互作用，將二者簡單解耦單獨求解可能導(dǎo)致流場分析結(jié)果與實際流場不符，此時考慮流固耦合作用下泵軸振動特性的數(shù)值計算將是非常必要的。流固耦合主要分為3類，直接耦合、順序耦合和同步耦合。直接耦合對計算機要求很高，目前難以應(yīng)用在流體機械中，順序耦合即按照既定順序?qū)ξ锢韴鲞M(jìn)行分別求解，也就是“弱耦合”，同步耦合即“雙向強耦合”，是指流體域和結(jié)構(gòu)域同時進(jìn)行求解計算。流固耦合問題可追溯到19世紀(jì)初的機翼和葉片的氣動彈性問題[1]。目前，流固耦合方法在水輪機和風(fēng)機等旋轉(zhuǎn)機械中已得到了比較廣泛的應(yīng)用[2-7]，近年來也逐漸在泵領(lǐng)域中得到應(yīng)用。Benra等[8]同時采用了單向耦合和雙向耦合的方法對離心泵轉(zhuǎn)子及流場進(jìn)行了流固耦合分析，并對水力激振位移進(jìn)行了實驗測量。Kato等[9]采用單向耦合的方法對多級離心泵中的振動噪聲進(jìn)行了研究。施衛(wèi)東等學(xué)者[10-14]對軸流泵內(nèi)部流場和葉輪結(jié)構(gòu)響應(yīng)進(jìn)行雙向順序流固耦合聯(lián)合求解，分析了葉輪葉片的應(yīng)力特性。劉厚林等[15]應(yīng)用雙向流固耦合方法對導(dǎo)葉式離心泵的外特性和內(nèi)流場進(jìn)行分析，研究了流固耦合作用對外特性影響的內(nèi)流機理。裴吉等[16,17]應(yīng)用同步求解，對離心泵內(nèi)流場和葉輪結(jié)構(gòu)響應(yīng)進(jìn)行數(shù)值研究。然而，國內(nèi)外研究中流固耦合方法主要集中在葉輪轉(zhuǎn)子上，對基于流固耦合作用的泵裝置中泵軸的振動特性的研究還屬于空白。

本文采用CFX14.5和Workbench14.5雙向同步求解的方法對一臺灌排雙向軸流泵裝置進(jìn)行研究，在考慮泵軸與出水流道流場的流固耦合作用下進(jìn)行瞬態(tài)非定常的數(shù)值模擬計算，并與非流固耦合下的瞬態(tài)流場進(jìn)行對比，分析泵軸與流場的流固耦合作用對內(nèi)流特性的影響。

1 物理模型及計算方法

1.1 幾何模型

以某灌排雙向進(jìn)出水軸流泵站的裝置模型為研究對象，該泵裝置的模型由進(jìn)水流道、葉輪、導(dǎo)葉和出水流道組成，泵裝置中進(jìn)水流道按水流方向可分為單向進(jìn)水流道和雙向進(jìn)水流道，單向進(jìn)水流道按形狀又有肘形彎管型、平面蝸殼(鐘形)型及其他型式，本文選用雙向鐘形進(jìn)水流道，其組成為進(jìn)口段、喇叭管和導(dǎo)水錐，出水流道采用肘形對拼式。泵軸連接導(dǎo)葉和電機并直接裸露在出水流道中，其中泵段模型的基本參數(shù)為：流量QDES=1 425 m3/h，轉(zhuǎn)速n=1 550 r/min，揚程H=2.76 m。葉輪外徑D2=300 mm，葉輪葉片數(shù)為3，導(dǎo)葉葉片數(shù)為5。利用Pro/E進(jìn)行三維建模，得到如圖1所示雙向軸流泵裝置模型。

圖1 雙向進(jìn)出水流道泵裝置三維造型Fig.1 The 3D model of axial-flow pumping system with two-way passage

1.2 網(wǎng)格劃分

水體域包括進(jìn)水流道、葉輪、導(dǎo)葉和出水流道。由于葉輪和導(dǎo)葉為泵的核心部件，其網(wǎng)格質(zhì)量與分布對軸流泵性能預(yù)測有著直接的影響[18-20]，因此采用結(jié)構(gòu)化與非結(jié)構(gòu)化混合網(wǎng)格，利用網(wǎng)格劃分軟件ICEM將葉輪和導(dǎo)葉劃分為六面體結(jié)構(gòu)網(wǎng)格，并將葉輪和導(dǎo)葉葉片表面邊界層及其附近網(wǎng)格進(jìn)行加密處理，其網(wǎng)格劃分情況如圖2(a)所示。結(jié)構(gòu)域只考慮泵軸部分，對泵軸采用ANSYS自動網(wǎng)格劃分功能，如圖2(b)所示。流體域網(wǎng)格單元總數(shù)為4 030 944，固體域網(wǎng)格單元數(shù)為5 073。

圖2 網(wǎng)格劃分Fig.2 The grids of fluid and structure field

1.3 邊界條件設(shè)置

流體域CFD計算中，采用SSTk-ω湍流模型求解雷諾時均方程，時間離散采用二階后向歐拉格式，處理轉(zhuǎn)子、定子問題使用的是包含轉(zhuǎn)動和靜止坐標(biāo)系的多重坐標(biāo)系方法，進(jìn)水流道采用質(zhì)量流量進(jìn)口，出水流道出口采用固定總壓，總壓設(shè)定為101 325 Pa，采用自由出流邊界條件。在進(jìn)水流道前池設(shè)置一個自由液面，設(shè)為對稱邊界條件。以穩(wěn)態(tài)RANS 計算結(jié)果作為非定常計算的初始流場，采用無滑移邊界條件，葉輪與靜止部件的交界面設(shè)置為“Transient Rotor stator”，而靜止部件的交界面設(shè)置為“None”。以葉輪旋轉(zhuǎn)3°為一個時間步長，每個時間步長為0.000 322 581 s。

對于泵軸結(jié)構(gòu)響應(yīng)的計算，首先給定模型泵軸結(jié)構(gòu)的材料為結(jié)構(gòu)鋼，其特性參數(shù)為：彈性模量E=200 GPa、泊松比μ=0.3、密度ρ=7 850 kg/m3。設(shè)置與流體計算相同的時間步長和初始時間。對泵軸與流體接觸的表面設(shè)置動網(wǎng)格，并將相應(yīng)的結(jié)構(gòu)表面與流體表面進(jìn)行對應(yīng)，設(shè)定流體向固體傳遞的數(shù)據(jù)類型為Total Forces，結(jié)構(gòu)影響流體區(qū)域的方式設(shè)定為Total Mesh Displacement。

1.4 雙向流固耦合計算方法

在求解過程中，流場非定常計算以定常計算結(jié)果為其初始條件。在雙向流固耦合問題中，需要同時求解流體域和固體域，并通過耦合作用面進(jìn)行數(shù)據(jù)傳遞。在每個物理時間步長上進(jìn)行若干次偶合迭代，具體計算過程如圖3所示。

圖3 雙向同步耦合計算過程Fig.3 The calculation process of two-way coupling method

2 計算結(jié)果與分析

針對雙向流道泵裝置在4種工況下(Q/QDES=0.8、1.0、1.05、1.1)進(jìn)行數(shù)值模擬計算，分析比較流固耦合作用對內(nèi)流場的影響作用，并選擇有代表性的結(jié)果進(jìn)行分析。

2.1 出水流道內(nèi)流場對比分析

泵軸與出水流道中的流體直接接觸，為了研究出水流道中流場的變化，以設(shè)計流量Q/QDES=1.0下葉輪旋轉(zhuǎn)的最后一個周期為研究對象，研究泵裝置出水流道中間截面上的流場分布情況，中間截面左端為出水流道出水口。

如圖4所示是不同時刻下流固耦合前后出水流道中間截面上的壓力和速度分布云圖對比，其中T是葉輪旋轉(zhuǎn)周期。不同時刻截面壓力耦合前后對比可以看出，當(dāng)時間t一定時，流固耦合前后壓力分布情況基本相似，從導(dǎo)葉出口到出水流道出口靜壓逐漸增大，在導(dǎo)葉出口靠近出水邊拐角處出現(xiàn)低壓區(qū)，在出水流道盲區(qū)靠近泵軸處出現(xiàn)高壓區(qū)，泵軸周圍壓力分布不均勻，這也是導(dǎo)致泵軸振動的原因之一，同時可以看到，流固耦合作用后出水流道中靜壓值相對較大的區(qū)域面積比耦合前要大，說明流固耦合作用使流場中擾動更大。隨著時間t的變化，流場中靜壓分布也不同。

通過對比耦合前后不同時刻截面速度云圖分布可知，對于一定的時間t，流固耦合作用對流速的影響不明顯，只在泵軸周圍發(fā)生了變化，可以看到，考慮流固耦合作用后，泵軸周圍的低速區(qū)域范圍增大。同時在導(dǎo)葉出口靠近出水邊拐角處出現(xiàn)最大流速，而在水流盲區(qū)流速幾乎接近于零，說明此處動能損失較大。不同時刻下，耦合前后流速分布總體趨勢不變。

2.2 壓力脈動對比分析

泵軸一端連接電機，一端連接導(dǎo)葉輪轂，在出水流道中由于泵軸的旋轉(zhuǎn)必然引起導(dǎo)葉出口及出水流道中流場的變化，而考慮到流場與泵軸相互作用前后的變化情況不同，因此在導(dǎo)葉出口截面，如圖2(a)中截面A-A上設(shè)置監(jiān)測點，分別讀取其在最后一個周期內(nèi)的壓力值，如圖5所示為3種工況下(Q/QDES=0.8、1.0、1.1)，流固耦合前后其壓力隨時間的變化情況。

由圖5可以看出，不同工況下，導(dǎo)葉出口處的壓力隨時間的變化不同，同時流固耦合前后壓力隨時間的變化情況也不相同。小流量下，流固耦合前后的壓力均隨時間呈周期性變化，但耦合前后兩曲線存在明顯的相位差，耦合前的曲線在波峰和波谷處相對滯后于耦合后的曲線，而耦合后壓力曲線的波峰位置上移，波谷位置下移，壓力脈動范圍增大，強度加大；設(shè)計流量和大流量下，壓力波動較為紊亂，沒有一定周期性，但同樣存在相位超前或滯后現(xiàn)象，同時耦合后的壓力曲線值高于耦合前的壓力曲線，說明考慮流固耦合作用后，泵軸發(fā)生了變形，其對稱性被破壞，從而導(dǎo)致導(dǎo)葉出口壓力波動增大，導(dǎo)葉出口壓力的分布情況發(fā)生變化。

圖4 出水流道截面中速度和靜壓分布Fig.4 The distribution of velocity and pressure in the outflow-passage section

圖5 不同工況下導(dǎo)葉出口壓力脈動對比Fig.5 Pressure fluctuation in diffuser outlet for different conditions

2.3 湍流動能對比分析

如圖6所示為不同工況下，考慮泵軸與出水流道流場流固耦合前后導(dǎo)葉出口面即截面A-A上的湍流動能TKE對比。湍流動能主要通過雷諾切應(yīng)力做功來給湍流提供能量，湍流動能獲得的能量越大，流體在流道內(nèi)的脈動損失越大[21]。從圖6中可以看出，隨著流量的增大，導(dǎo)葉出口截面的湍動能先減小后增大，小流量下脈動能量最大，達(dá)到了2.278 m2/s-2，且湍動能梯度更大，設(shè)計流量下截面的湍動能最小，說明該流量下泵內(nèi)流動更加平緩。在小流量和設(shè)計流量下，耦合前后截面上TKE的分布呈明顯5片花瓣形，這與導(dǎo)葉數(shù)相對應(yīng)，說明湍動能還受導(dǎo)葉分布的影響，而大流量下這種影響消失。在小流量下，耦合前后的湍動能分布隨相似，但考慮流固耦合后，截面上TKE比耦合前要大，且變化梯度也增大，說明泵軸的變形導(dǎo)致流場更加紊亂，脈動增強。而在設(shè)計流量和大流量下泵軸與流場的耦合作用對湍動能的影響較小。

圖6 不同工況下導(dǎo)葉出口截面的湍動能分布Fig.6 TKE in diffuser outlet for different conditions

2.4 外特性對比

為進(jìn)一步說明模擬的可行性，取4個工況(Q/QDES=0.8、1.0、1.05、1.1)下，流固耦合前后葉輪旋轉(zhuǎn)的最后一個周期的時均揚程，作如圖7中的流量-揚程曲線，并與裝置試驗值進(jìn)行對比，本文試驗值采用10倍模型泵的某實際泵裝置運行時的揚程值，經(jīng)相似定律即模型泵與真機的nD(n為轉(zhuǎn)速，D為葉輪外徑)相等換算而來，具有較高的可信度。

圖7 流量-揚程曲線Fig.7 The flow-head curve

由圖7可知，流固耦合前后揚程的變化趨勢均與試驗一致，說明該數(shù)值模擬計算方法具有一定的可行性，此外，耦合前后的揚程值都大于試驗值，這是因為數(shù)值模擬的泵裝置流場未考慮泄露損失和機械損失等客觀因素；在考慮了泵軸與出水流道流場的流固耦合作用后的揚程略低于非流固耦合下的揚程，即更加接近試驗值，設(shè)計工況下，流固耦合與非流固耦合相對于試驗的誤差分別為5.2%、7.1%，由此說明考慮泵軸與流場的流固耦合作用后的結(jié)果更加接近真實流場。

3 結(jié) 語

采用雙向流固耦合的方法，通過ANASYS軟件首次對一臺灌排雙向軸流泵裝置進(jìn)行中的泵軸與出水流道流場進(jìn)行數(shù)值模擬計算，并分析比較了流固耦合前后其出水流道中壓力和速度以及導(dǎo)葉出口湍流動能的變化情況，并將其外特性與試驗結(jié)果進(jìn)行了對比，得到如下結(jié)論。

(1)考慮流固耦合作用后，出水流道中間截面上的靜壓分布發(fā)生了變化，靜壓相對較大的區(qū)域面積增大，說明泵軸的形變對流場靜壓分布產(chǎn)生了干擾，也反映了泵軸與流場的耦合作用影響確實存在，但耦合作用對流速分布的影響不明顯，只在泵軸周圍產(chǎn)生微小變化。

(2)不同工況下，流固耦合都使得導(dǎo)葉出口處的壓力脈動發(fā)生了相位偏移，且在設(shè)計流量和大流量下壓力脈動的強度也明顯增大；在小流量下，耦合作用增強了導(dǎo)葉出口的紊流，使其湍流動能增大，設(shè)計流量和大流量下這種影響較弱。

(3)流固耦合前后其流場的分布規(guī)律相同，耦合作用并不能改變流場的整體流動趨勢。耦合前后泵裝置的揚程變化趨勢一致，但耦合后的揚程更接近試驗值，設(shè)計工況下，相對于試驗的誤差分別為5.2%、7.1%，從而驗證了數(shù)值模擬的可靠性。

因此，對大型灌排雙向軸流泵裝置的內(nèi)流分析應(yīng)盡量采用流固耦合計算提高分析精度。

□

[1] 邢景棠,周 盛,崔爾杰. 流固耦合力學(xué)概述[J]. 力學(xué)進(jìn)展,1997,27(1):20-39.

[2] 羅永要,王正偉,梁權(quán)偉. 混流式水輪機轉(zhuǎn)輪動載荷作用下的應(yīng)力特性[J]. 清華大學(xué)學(xué)報(自然科學(xué)版),2005,45(2):235-237,257.

[3] 鄭小波,羅興琦,鄔海軍. 軸流式葉片的流固耦合振動特性分析[J]. 西安理工大學(xué)學(xué)報,2005,21(4):342-346.

[4] 毛 軍,楊立國,郗艷紅. 大型軸流風(fēng)機葉片的氣動彈性數(shù)值分析研究[J]. 機械工程學(xué)報,2009,45(11):133-139.

[5] Guerri O, Hamdouni A, Sakout A. Fluid structure interaction of wind turbine airfoils[J]. Wind Engineering, 2008,32(6):539-557.

[6] Pei Ji, Dohmen H J, Yuan S Q, et al. Investigation of unsteady flow-induced impeller oscillations of a single-blade pump under off-design conditions [J]. Journal of Fluids and Structures, 2012,35:89-104.

[7] Xiao Ruofu, Wang Zhengwei, Luo Yongyao. Dynamic stresses in a Francis turbine runner based on fluid-structure interaction analysis[J]. Tsinghua Science & Technology, 2008,13(5):587-592.

[8] Benra F-K, Dohmen H J.Comparison of pump impeller orbit curves obtained by measurement and FSI simulation[ C]∥2007 ASME Pressure Vessels and Piping Division Conference, 2007.

[9] Kato C, Yamade Y, Wang Hong, et al. Prediction of the noise from a multi-stage centrifugal pump[C]∥ASME Fluids Engineering Division Summer Meeting, June 19-23, 2005, TX, USA.

[10] 施衛(wèi)東,王國濤,張德勝,等. 基于流固耦合的軸流泵葉片應(yīng)力特性[J]. 排灌機械工程學(xué)報,2013,31(9):737-740,757.

[11] 潘 旭,李 成,鐵 瑛,等. 軸流泵葉片流固耦合強度分析[J]. 水力發(fā)電學(xué)報,2012,31(4):221-226,237.

[12] 黃浩欽,劉厚林,王 勇,等. 基于流固耦合的船用離心泵轉(zhuǎn)子應(yīng)力應(yīng)變及模態(tài)研究[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(15):98-105.

[13] 張 新,鄭 源,毛秀麗,等. 基于流固耦合的軸流泵葉輪強度分析[J]. 水電能源科學(xué),2014,32(7):137-139,150.

[14] 朱 利,楊昌明,鄭 軍,等. 基于流固耦合的軸流泵葉輪結(jié)構(gòu)分析[J]. 流體機械,2013,41(3):20-23,40.

[15] 劉厚林,徐 歡,吳賢芳,等. 流固耦合作用對離心泵內(nèi)外特性的影響[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(13):82-87.

[16] 裴 吉,袁壽其,袁建平. 流固耦合作用對離心泵內(nèi)部流場影響的數(shù)值計算[J]. 農(nóng)業(yè)機械學(xué)報,2009,40(12):107-112.

[17] Pei Ji, Yuan Shouqi, Benra F-K, et al. Numerical prediction of unsteady pressure field within the whole flow passage of a radial single-blade pump [J]. ASME Journal of Fluids Engineering, 134: 101103, 2012.

[18] 張德勝，施衛(wèi)東，張 華，等.不同湍流模型在軸流泵性能預(yù)測中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(1)：66-70.

[19] 施衛(wèi)東，張 華，陳 斌，等. 不同葉頂間隙下的軸流泵內(nèi)部流場數(shù)值計算[J]. 排灌機械工程學(xué)報，2010，28(5)：374-377.

[20] 張德勝，施衛(wèi)東，張 華，等. 軸流泵葉輪端壁區(qū)流動特性數(shù)值模擬[J]. 農(nóng)業(yè)機械學(xué)報，2012，43(3)：73-77.

[21] 葉道星,王 洋. 離心泵葉輪內(nèi)部湍流動能及耗散率分析[J]. 中國農(nóng)村水利水電,2012,(4):84-88.