基于流動模擬的多級離心泵葉輪應力變形分析*

倪羽鳳 覃 磊 杜俠明 王 君

(中國石油大學(華東)化學工程學院)

多級離心泵應用領域非常廣泛,由于近年來工業的飛速發展,離心泵向著高轉速、大容量發展，對葉輪的強度要求越來越高。葉輪的輕微破壞(磨損)將導致葉輪的入口角、出口角變化，影響泵的水利性能；嚴重破壞(斷裂)將直接影響泵的穩定運行，因此對多級離心泵葉輪進行強度分析，對保證泵的穩定運行和指導泵結構設計具有重要意義。文獻[1,2]對葉輪內液體流動情況進行數值模擬；文獻[3]對葉輪受到離心力作用的情況進行了分析，文獻[4,5]對葉輪在受到介質壓力和離心力情況下進行了受力變形分析；文獻[6，7]對葉輪在受到介質壓力，離心力和剪應力的情況下進行了分析；但是以上文獻沒有對葉輪的全部受力進行耦合分析，載荷分析過于簡化，與實際工作狀態相差較大。

筆者考慮多級離心泵末級葉輪承受的介質力最大，所以對十級離心泵末級葉輪在空轉和實際運行兩種工況進行分析。分別考慮葉輪只受離心力和在介質力、離心力、前后蓋板摩擦力耦合作用下的受力變形,提出了一種基于FLUENT和ABAQUS的多級離心泵葉輪有限元分析方法，能夠準確得到葉輪的應力變形規律。

1 葉輪介質力載荷和有限元模型

1.1葉輪介質力載荷

在FLUENT中，通過對多級離心泵工作過程進行流動模擬，得到葉輪內的介質力分布如圖1所示。葉輪前后蓋板的壓力分布相似，壓力由葉輪入口到葉輪出口逐漸變大；葉片的壓力分布由入口到出口逐漸增大，且工作面的壓力高于非工作面的壓力。

圖1 葉輪介質力分布

1.2葉輪有限元模型的建立

葉輪主要設計參數為：

材料 ZG1Cr13NiMo

揚程 122m

設計流量 158m3/h

轉速 2 980r/min

葉片進口直徑 101mm

葉片出口直徑 308mm

葉片厚度 5mm

葉片進口角 10°

葉片出口角 26°

葉片包角 130°

在ABAQUS中建立葉輪模型，對葉輪進行網格劃分，網格單元為C3D4，設置單元尺寸為2，共將模型離散為單元147 962個，節點238 942個，如圖2所示。

圖2 葉輪模型

2 離心泵葉輪受力變形分析

2.1葉輪空轉情況下的應力變形分布

葉輪空轉時重力影響很小，只考慮離心力的作用。在葉輪輪轂內側表面施加全約束，根據葉輪的轉速2 980r/min，繞軸向給整個葉輪施加轉速為ω=312.1r/s的載荷。

2.1.1應力分析

空轉情況下的應力分布如圖3所示，圖片顯示為放大400倍的效果。最大的應力為周向應力，位于靠近葉輪入口，葉片與前蓋板交界處，大小為15.6MPa；徑向應力在葉片與前、后蓋板交界處較其它地方明顯；周向應力分布比較均勻，在葉片與前蓋板交界處明顯；軸向應力在葉片與后蓋板交界處靠近葉輪入口位置較大。

圖3 葉輪空轉情況下的應力分布

2.1.2變形分析

葉輪空轉工況下的變形分布如圖4所示，圖片顯示為放大400倍的效果。葉輪的變形由葉輪入口到出口逐漸變大，到達葉輪外緣時達到最大，最大變形為4.3×10-3mm。葉輪前后蓋板的變形明顯不同，前蓋板變形較后蓋板變形明顯。

圖4 空轉情況下的葉輪變形分布

2.2葉輪實際運轉工況下應力變形分布

葉輪運轉時，受到介質力、離心力和摩擦力的共同作用。對葉輪施加旋轉離心力、葉輪內部介質力載荷、摩擦力，對其進行耦合受力變形分析。

摩擦力的計算式為[8]：

式中λ——摩擦系數

γ——介質重度。

得到摩擦力分布圖如圖5所示，摩擦力的方向與介質流動方向相同，即沿著流道的切線方向由葉輪入口指向葉輪出口。其中，p1=9.83kPa，p2=3.17kPa，p3=1.04kPa。

圖5 摩擦力分布

2.2.1應力分析

實際運轉工況的應力分布如圖6所示，圖片顯示為放大300倍的效果。軸向應力最大，最大應力點位于葉片與前蓋板交界處，大小為24MPa；周向應力最??；徑向應力在葉片與前后蓋板交界處較明顯，其他部位應力分布均勻；軸向應力在葉片靠近前蓋板的部位比靠近后蓋板的部位明顯；前后蓋板的應力分布規律不同。

圖6 運轉工況下的葉輪應力分布

2.2.2變形分析

葉輪運轉工況下的變形如圖7所示。圖片顯示為放大400倍的效果。葉輪的總變形由葉輪入口到出口逐漸變大，在葉輪外緣達到最大，最大變形分別位于前蓋板靠近葉片非工作面處和后蓋板與葉片交界處，大小為16.80μm；軸向最大變形量為16.77μm，與總變形量近似，對總變形的影響最大；徑向和周向的最大變形分別為1.20μm和1.70μm，對總變形的影響較??；葉輪前后蓋板的變形明顯不同，前蓋板變形比后蓋板變形明顯。

圖7 運轉工況下的葉輪變形分布

3 結論

3.1通過對多級離心泵內部流動情況的數值模擬，得到了多級離心泵內的壓力分布，以流場為載荷邊界條件，分析了葉輪在運轉工況下的應力和變形分布，得到了一種基于流場模擬的多級離心泵葉輪受力變形分析方法。

3.2葉輪運行工況下，葉片與前后蓋板連接處的應力值較大，其他部位應力分布均勻，不存在應力集中現象；葉輪軸向應力較大，變形主要由軸向變形引起；離心力對葉輪的應力影響較大。

3.3葉輪空轉時，只受離心力的作用，應力在葉片與前后蓋板交界處靠近葉輪入口位置明顯，存在應力集中現象，此處葉輪結構有待進一步改進；葉輪的變形從葉輪入口到葉輪出口逐漸變大，在葉輪邊緣達到最大；前后蓋板變形分布明顯不同，前蓋板的變形較后蓋板明顯。

[1] 王志堅, 佟亮, 李璐璐, 等.基于CFD的離心泵內部三維流動數值模擬和性能預測[J].流體機械，2012，40(6)：14～18.

[2] Cui B L, Lin Y G, Jin Y Z .Numerical Simulation of Flow in Centrifugal Pump with Complex Impeller[J].Journal of Thermal Science,2011,20(1):47～52.

[3] 賈寧寧，楊昌明，張圣,等．基于 ANSYS Workbench 的離心泵葉輪強度分析[J]．西華大學學報(自然科學版)，2011，30(6)：52～55.

[4] 周玉豐.基于ANSYS的離心泵葉輪結構有限元分析[J].自動化與儀器儀表， 2011，(4)：136 ～137.

[5] 于保敏, 劉宏杰.離心葉輪的有限元分析[J].機械工程師，2007，(7)：60～61.

[6] 董曉嵐.基于FLUENT和ABAQUS離心泵葉輪強度和剛度的研究[J].機械設計與制造，2011，(6)：219～220.

[7] 高丙朋，陳飛.廢油處理系統中的離心泵應用研究[J].化工自動化及儀表,2012,39(10): 1296～1299.

[8] 關醒凡.現代泵設地手冊[M].北京：中國宇航出版社，1995.