流固耦合作用對(duì)螺旋離心泵流場(chǎng)影響的數(shù)值分析

　　為研究流固耦合作用對(duì)螺旋離心泵內(nèi)部流場(chǎng)的影響，以ZJ200-25型螺旋離心泵為研究對(duì)象，應(yīng)用計(jì)算流體力學(xué)軟件CFX12.1和有限元軟件ANSYS Workbench對(duì)螺旋離心泵進(jìn)行了考慮內(nèi)部流場(chǎng)和結(jié)構(gòu)相互作用的兩場(chǎng)交替聯(lián)合求解，并以相同的設(shè)置對(duì)不考慮流固耦合作用的內(nèi)部流場(chǎng)進(jìn)行了計(jì)算。對(duì)比分析兩種計(jì)算方法得到的流場(chǎng)發(fā)現(xiàn)，考慮流固耦合作用后:泵進(jìn)口處的壓力波動(dòng)加劇，出口壓力波動(dòng)強(qiáng)度降低，進(jìn)出口波動(dòng)頻率不變，但相位有改變；速度分布滿足共同的規(guī)律，但在小流量工況下差別最大，同時(shí)在葉片與隔舌處于相干位置時(shí)，隔舌附近流動(dòng)也受到明顯影響；徑向力大小隨時(shí)間變化規(guī)律一致，但波動(dòng)范圍減��；考慮流固耦合作用后葉輪流道內(nèi)流動(dòng)更加不對(duì)稱。

　　螺旋離心泵是一種傳統(tǒng)的無(wú)堵塞泵，其核心部件為半開(kāi)式螺旋形葉輪。這種泵葉片少，葉片薄而相對(duì)扭曲，包角大，結(jié)構(gòu)特殊。螺旋離心泵運(yùn)行過(guò)程中，在流場(chǎng)內(nèi)復(fù)雜的水力激勵(lì)作用下，葉輪會(huì)發(fā)生變形及振動(dòng)，變形后的葉輪又會(huì)引起流場(chǎng)的改變。以往關(guān)于螺旋離心泵內(nèi)部流動(dòng)的數(shù)值研究中，都沒(méi)有考慮到結(jié)構(gòu)變形和振動(dòng)對(duì)流場(chǎng)的反作用，而雙向流固耦合方法已經(jīng)被用于普通離心泵內(nèi)流場(chǎng)的研究中，但還沒(méi)有在泵領(lǐng)域得到推廣，成果有限。

　　目前應(yīng)用較多的是將流場(chǎng)載荷加載到結(jié)構(gòu)上的單向流固耦合方法，這種方法忽略了阻尼作用，準(zhǔn)確度低，常用于求解結(jié)構(gòu)響應(yīng)的靜力學(xué)分析。本文使用CFD軟件CFX12.1及有限元分析軟件ANSYS Workbench對(duì)螺旋離心泵進(jìn)行雙向流固耦合數(shù)值模擬，對(duì)比研究耦合與不耦合各工況下流場(chǎng)壓力、速度及葉輪受力的差別，分析流固耦合作用對(duì)流場(chǎng)產(chǎn)生的影響。

1、計(jì)算模型與模擬方法

　　1.1、模型參數(shù)

　　螺旋離心泵模型設(shè)計(jì)參數(shù)為:流量Qd=200m3/h，揚(yáng)程H=25m，轉(zhuǎn)速n=1450r/min，比轉(zhuǎn)數(shù)ns=116，葉片數(shù)Z=2。葉輪進(jìn)口直徑D1=160mm，出口直徑D2=328mm，出口寬度b2=70mm，蝸殼基圓直徑D3=360mm，葉輪外緣與吸水殼體內(nèi)壁間隙為0.9mm。葉輪材料為不銹鋼316L，密度8030kg/m3，彈性模量175GPa，泊松比0.4。

　　1.2、網(wǎng)格生成

　　將螺旋離心泵的三維模型導(dǎo)入網(wǎng)格劃分軟件ICEM生成計(jì)算區(qū)域的網(wǎng)格�？紤]到螺旋離心泵結(jié)構(gòu)特殊，葉片邊緣存在一尖角區(qū)域，且計(jì)算區(qū)域中存在比較復(fù)雜的流動(dòng)，如二次流、交叉流以及小流量工況下的分離旋渦等流動(dòng)現(xiàn)象，采用六面體網(wǎng)格單元不能保證網(wǎng)格質(zhì)量和更好地反映這些流動(dòng)。故采用適應(yīng)范圍更廣的四面體網(wǎng)格單元對(duì)全流道模型進(jìn)行網(wǎng)格劃分，并對(duì)間隙內(nèi)流體進(jìn)行局部網(wǎng)格加密處理，流體網(wǎng)格總數(shù)為1399178(圖1)。轉(zhuǎn)子固體有限元網(wǎng)格在WorkbenchTransientStructure模塊中生成，單元數(shù)為125662(圖2)。

圖1 流體區(qū)域網(wǎng)格

(a)葉輪計(jì)算區(qū)域網(wǎng)格(b)蝸殼計(jì)算區(qū)域網(wǎng)格

圖2 固體區(qū)域網(wǎng)格

　　1.3、數(shù)學(xué)模型及模擬方法

　　數(shù)值模擬在ANSYSCFX12.1全隱式耦合多網(wǎng)格線性求解器中進(jìn)行。采用SSTk-ω湍流模型進(jìn)行非定常全流場(chǎng)計(jì)算，采用適宜于瞬態(tài)求解的SIMPLEC算法進(jìn)行求解。進(jìn)出口邊界條件分別定義為速度進(jìn)口和自由出流，設(shè)置參考?jí)毫�為大氣壓，給定不同工況下的進(jìn)口速度進(jìn)行計(jì)算。固體壁面定義無(wú)滑移邊界條件，粗糙度0.025mm。將進(jìn)口與葉輪、葉輪與蝸殼間的交界面設(shè)置為T(mén)ransientRotorStator。以葉輪旋轉(zhuǎn)9°作為一個(gè)時(shí)間步，每個(gè)時(shí)間步長(zhǎng)ΔT=0.001035s，每個(gè)旋轉(zhuǎn)周期包含40個(gè)時(shí)間步。選擇5個(gè)旋轉(zhuǎn)周期作為采樣對(duì)象，200步的總時(shí)間T=0.207s。

　　在雙向流固耦合計(jì)算的基礎(chǔ)上，以相同的時(shí)間步長(zhǎng)及設(shè)置，進(jìn)行了不考慮流固耦合作用的解耦計(jì)算，以此結(jié)果作為參照，能夠更直觀地對(duì)比得出流固耦合作用對(duì)流場(chǎng)所產(chǎn)生的影響。

結(jié)論

　　(1)流固耦合作用使得進(jìn)口壓力波動(dòng)強(qiáng)度增大，出口壓力波動(dòng)強(qiáng)度減弱，隔舌處波動(dòng)紊亂，但流固耦合前后的流場(chǎng)也滿足共同的規(guī)律，耦合作用并未改變整體流動(dòng)趨勢(shì)。

　　(2)不同工況下，流固耦合作用在非設(shè)計(jì)流量下影響最明顯，小流量下對(duì)壓力影響最明顯，大流量下對(duì)徑向力影響最明顯；一個(gè)旋轉(zhuǎn)周期內(nèi)，當(dāng)葉片與隔舌處于相干位置時(shí)，耦合前后的流動(dòng)差異最大。

　　(3)由于葉輪變形以及質(zhì)心偏移，破壞了轉(zhuǎn)子結(jié)構(gòu)的對(duì)稱性，對(duì)泵進(jìn)口處液流形成擾動(dòng)，并對(duì)下游流道內(nèi)流動(dòng)形成不同程度的影響，這種影響隨時(shí)間及空間積累，會(huì)在流道內(nèi)局部產(chǎn)生明顯的差別。