劉 棟，宋 濤，劉子凌，黃 凱，黃一峰

（1.江蘇大學(xué) 能源與動(dòng)力工程學(xué)院，江蘇鎮(zhèn)江 212013；2.江蘇永一泵業(yè)科技集團(tuán)有限公司，江蘇揚(yáng)州 225831）

0 引言

高溫泵廣泛應(yīng)用于核電、火電、冶金等領(lǐng)域[1-2]，在運(yùn)行過程中，高溫介質(zhì)對(duì)泵的結(jié)構(gòu)影響較大，容易出現(xiàn)葉輪變形、口環(huán)變形、磨損等問題[3-7]，嚴(yán)重影響泵運(yùn)行的可靠性。在溫度載荷和壓力載荷的綜合作用下，口環(huán)極易發(fā)生變形，造成動(dòng)靜部件的磨損，嚴(yán)重時(shí)甚至導(dǎo)致轉(zhuǎn)子部件的咬死。適當(dāng)增大口環(huán)間隙可以避免上述問題，但會(huì)導(dǎo)致泵的容積效率下降，降低口環(huán)間隙內(nèi)流體的支承剛度，影響高溫泵轉(zhuǎn)子運(yùn)行的穩(wěn)定性。

目前關(guān)于口環(huán)間隙的研究，主要集中在離心泵性能和結(jié)構(gòu)強(qiáng)度等方面。DAQIQSHIRAZI等[8]采用數(shù)值模擬與試驗(yàn)相結(jié)合的方法研究口環(huán)間隙對(duì)離心泵性能的影響，發(fā)現(xiàn)間隙增大時(shí)，泵的效率明顯減小。張景等[9]研究了不同口環(huán)間隙對(duì)離心泵性能及壓力脈動(dòng)的影響，發(fā)現(xiàn)隨著口環(huán)間隙增大離心泵的效率及壓力脈動(dòng)幅值均有所下降。ZHAO等[10]采用數(shù)值計(jì)算方法，研究了3種口環(huán)間隙方案對(duì)離心泵性能的影響，認(rèn)為前口環(huán)間隙對(duì)性能的影響比后口環(huán)間隙的影響大。崔哲等[11]通過實(shí)驗(yàn)測(cè)試研究了口環(huán)間隙對(duì)離心泵性能和振動(dòng)的影響，發(fā)現(xiàn)隨著間隙的增加，泵的揚(yáng)程下降，流致振動(dòng)能量逐漸往高頻移動(dòng)。牟介剛等[12]研究了不同口環(huán)結(jié)構(gòu)對(duì)離心泵性能的影響，通過對(duì)比分析，發(fā)現(xiàn)鋸齒形口環(huán)結(jié)構(gòu)可以提升離心泵的水力效率。盡管口環(huán)間隙的研究較多，但關(guān)于其對(duì)轉(zhuǎn)子運(yùn)行穩(wěn)定性的研究卻鮮有報(bào)道，而在其他領(lǐng)域有不少相關(guān)研究[13-14]，如汽輪機(jī)的汽封和口環(huán)作用一致。趙三星等[15]研究了汽封力對(duì)軸承靜態(tài)平衡點(diǎn)及動(dòng)特性系數(shù)的影響，發(fā)現(xiàn)汽封力對(duì)系統(tǒng)穩(wěn)定性的影響較為復(fù)雜，與密封進(jìn)出口壓力、進(jìn)口預(yù)旋速度、密封參數(shù)等因素有關(guān)。吳立明等[16]用等效剛度和阻尼對(duì)間隙汽封力進(jìn)行描述，分析汽封力對(duì)汽輪機(jī)組軸系穩(wěn)定性的影響規(guī)律，發(fā)現(xiàn)考慮汽封力作用，可以明顯提高轉(zhuǎn)子的穩(wěn)定性。

本文將通過數(shù)值模擬研究口環(huán)間隙對(duì)高溫泵水力性能的影響規(guī)律，以此為基礎(chǔ)，基于流熱固耦合，分析口環(huán)間隙對(duì)葉輪口環(huán)處形變量的影響；計(jì)算口環(huán)和軸承動(dòng)力特性系數(shù)，研究口環(huán)間隙對(duì)轉(zhuǎn)子臨界轉(zhuǎn)速和位移響應(yīng)的影響，分析口環(huán)間隙對(duì)轉(zhuǎn)子穩(wěn)定性的影響規(guī)律，為高溫泵口環(huán)的設(shè)計(jì)計(jì)算提供參考。

1 數(shù)值計(jì)算方法與實(shí)驗(yàn)驗(yàn)證

1.1 計(jì)算模型

高溫泵的設(shè)計(jì)參數(shù)：流量Q＝23.6 m3/h，轉(zhuǎn)速 n＝1 450 r/min，工作溫度 t＝230 ℃，葉輪進(jìn)口直徑 D1＝80 mm，出口直徑 D2＝209 mm，出口寬度b1＝12 mm，葉片數(shù)Z為5。230 ℃水的密度為826.8 kg/m3，黏度為 0.000 108 6 Pa·s，比熱容為4 642 J/（kg·K）；25 ℃水的密度為 997 kg/m3，黏度為 0.000 902 75 Pa·s，比熱容為 4 179 J/（kg·K）。圖1示出流體計(jì)算域，包括進(jìn)水管、葉輪、蝸殼和出水管水體域。

圖1 流體計(jì)算域Fig.1 Computational domains of fluid

固體計(jì)算域?yàn)檗D(zhuǎn)子部件，由螺母、葉輪、軸組成，其中軸承1為導(dǎo)軸承、軸承2為滾動(dòng)軸承，且轉(zhuǎn)子部件一端與聯(lián)軸器相連，轉(zhuǎn)子部件的幾何模型如圖2所示。

圖2 固體計(jì)算域Fig.2 Computational domains of solid

為研究口環(huán)間隙對(duì)高溫泵運(yùn)行穩(wěn)定性的影響，基于API 610標(biāo)準(zhǔn)確定口環(huán)最小間隙為0.38 mm，在此基礎(chǔ)上，通過適當(dāng)增加間隙來進(jìn)行數(shù)值模擬，間隙分別為 b＝0.4，0.6，0.8，1.0 mm。

1.2 控制方程與邊界條件

采用CFX軟件對(duì)高溫泵內(nèi)部流場(chǎng)進(jìn)行數(shù)值計(jì)算，控制方程為三維定常不可壓縮雷諾時(shí)均N-S方程，湍流模型采用RNG k-ε湍流模型，流體計(jì)算區(qū)域熱傳遞模型為Thermal Energy模型，根據(jù)泵的流量確定進(jìn)口速度，出口設(shè)置為壓力出口，近壁面區(qū)域流動(dòng)采用標(biāo)準(zhǔn)壁面函數(shù)法處理，計(jì)算的收斂精度為10-4。

轉(zhuǎn)子結(jié)構(gòu)主要包括葉輪和泵軸兩部分，采用ANSYS Workbench軟件對(duì)流場(chǎng)、溫度場(chǎng)和固體結(jié)構(gòu)進(jìn)行耦合分析，轉(zhuǎn)子部件的材料為35CrMo，材料特性參數(shù)：密度7 870 kg·m3，楊氏模量213 GPa，泊松比0.286。建立流場(chǎng)分析模塊CFX與熱應(yīng)力分析模塊Thermal-Stress之間的數(shù)據(jù)傳遞，將流體域的壓力載荷與溫度載荷施加給固體域，實(shí)現(xiàn)流熱固耦合。對(duì)轉(zhuǎn)子系統(tǒng)施加重力載荷，離心力載荷通過在轉(zhuǎn)子部件施加旋轉(zhuǎn)角速度實(shí)現(xiàn)，整個(gè)轉(zhuǎn)子部件受到軸承的約束作用，軸承接觸施加剛度、阻尼，在與聯(lián)軸器相連的軸端設(shè)置為固定約束。圖3示出轉(zhuǎn)子部件載荷和約束的設(shè)置。

圖3 轉(zhuǎn)子部件載荷和約束的設(shè)置Fig.3 Setting of rotor component load and constraints

1.3 試驗(yàn)驗(yàn)證

為了驗(yàn)證數(shù)值計(jì)算方法的可靠性，本文完成了高溫泵輸送常溫清水的外特性驗(yàn)證，并將模擬結(jié)果與試驗(yàn)值進(jìn)行對(duì)比，結(jié)果如圖4所示。

圖4 高溫泵外特性曲線Fig.4 External characteristic curve of high temperature pump

模型泵口環(huán)間隙為0.4 mm，數(shù)值計(jì)算的模型與試驗(yàn)?zāi)Ｐ鸵恢?。由圖4可以看出，數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果相比存在一些差異，這是由于數(shù)值模擬忽略了機(jī)械損失，導(dǎo)致與試驗(yàn)結(jié)果有所差異，但模擬結(jié)果與試驗(yàn)結(jié)果的趨勢(shì)吻合較好，因此采用CFX數(shù)值模擬的方法具有可靠性。以此為基礎(chǔ)，研究模型泵輸送高溫介質(zhì)時(shí)的性能，結(jié)果表明：高溫泵的揚(yáng)程和效率略有增加，不同工況下的變化規(guī)律與常溫一致。

2 軸承及口環(huán)的動(dòng)力特性系數(shù)

2.1 導(dǎo)軸承的剛度阻尼

導(dǎo)軸承結(jié)構(gòu)參數(shù)：軸頸直徑24 mm，寬徑比0.8，寬度19.2 mm，半徑間隙0.014 mm。通過求解 Reynolds方程[17]得到：

式中 R，φ，z——柱坐標(biāo)；

h——油膜厚度；

η——潤(rùn)滑油黏度；

p——油膜壓力；

Ω——軸頸轉(zhuǎn)速。

考慮溫度的影響，計(jì)算導(dǎo)軸承高溫與常溫工況下的動(dòng)力特性系數(shù)，如表1所示。

表1 導(dǎo)軸承的剛度、阻尼Tab.1 Stiffness and damping of guide bearing

本文中的滾動(dòng)軸承為7307B角接觸球軸承，查閱機(jī)械設(shè)計(jì)手冊(cè)第四卷[18]，計(jì)算滾動(dòng)軸承的剛度K＝267 000 N/mm，由于滾動(dòng)軸承的阻尼較小，本文忽略阻尼對(duì)軸承動(dòng)力特性的影響。

2.2 口環(huán)動(dòng)力特性系數(shù)

在流體域口環(huán)間隙處設(shè)置速度、壓力監(jiān)測(cè)點(diǎn)，運(yùn)用上述公式計(jì)算口環(huán)動(dòng)力特性系數(shù)。如表2所示，隨著口環(huán)間隙增大，口環(huán)剛度、阻尼均減小，這是由于口環(huán)泄漏量增大，導(dǎo)致口環(huán)間隙壓差和流速減小。

表2 口環(huán)動(dòng)力特性系數(shù)Tab.2 Dynamic characteristic coefficient of wear-ring

3 結(jié)果分析

3.1 口環(huán)間隙對(duì)高溫泵外特性的影響

本文采用 CFX 分別在 0.6Q，0.8Q，1.0Q，1.2Q，1.4Q 5個(gè)工況下進(jìn)行數(shù)值模擬，圖5為不同口環(huán)間隙變化方案下，泵的揚(yáng)程和效率隨工況的變化曲線。從圖中可以看出，泵的揚(yáng)程和效率均隨著間隙值的增大有所下降，當(dāng)口環(huán)間隙從0.4 mm變化至0.6 mm時(shí)，泵的揚(yáng)程和效率變化較??；進(jìn)一步增大到1.0 mm時(shí)，泵的揚(yáng)程和效率明顯降低，在額定工況下分別下降5.95%和6.26%。造成上述結(jié)果的原因是口環(huán)間隙變大，泵腔內(nèi)高壓流體泄漏量增加，導(dǎo)致高溫泵的容積損失增大。

圖5 不同口環(huán)間隙下?lián)P程、效率曲線對(duì)比Fig.5 Comparison of head and efficiency with different clearances

3.2 口環(huán)間隙對(duì)轉(zhuǎn)子部件形變量的影響

如圖6所示，葉輪口環(huán)處的形變量最大，這是由于口環(huán)間隙處的激振流體壓力的作用，導(dǎo)致易發(fā)生應(yīng)力集中。

圖6 口環(huán)間隙為0.4 mm時(shí)轉(zhuǎn)子形變量云圖Fig.6 Nephogram of rotor deformation for clearance of 0.4 mm

當(dāng)口環(huán)的形變量大于口環(huán)間隙時(shí)，易導(dǎo)致高溫泵無法運(yùn)行。圖7為不同工況下葉輪口環(huán)處最大形變量隨口環(huán)間隙的變化曲線，由圖可以看出，隨著口環(huán)間隙增大，葉輪口環(huán)處的形變量不斷減小，且間隙值為0.6 mm時(shí)，口環(huán)處的形變量降低最明顯，這是由于口環(huán)泄漏量增大導(dǎo)致泵腔內(nèi)流體壓力減小，從而減小了葉輪口環(huán)處的形變量。當(dāng)口環(huán)間隙b大于0.4 mm時(shí)，葉輪口環(huán)的形變量都小于0.5 b。

圖7 不同口環(huán)間隙下葉輪口環(huán)處的形變量Fig.7 Deformation of impeller wear-ring with different clearances

由上述結(jié)果可知，口環(huán)間隙為1.0 mm時(shí)，泵的揚(yáng)程和效率下降明顯，而葉輪口環(huán)處的形變量較0.8 mm時(shí)變化不大，因此，綜合性能和口環(huán)形變量?jī)煞矫婵紤]，認(rèn)為口環(huán)間隙大于0.8 mm時(shí)的研究意義不大，故下文只對(duì)口環(huán)間隙為0.4，0.6，0.8 mm時(shí)的模型泵進(jìn)行研究。

3.3 轉(zhuǎn)子部件的臨界轉(zhuǎn)速和振型

通過Workbench對(duì)高溫泵轉(zhuǎn)子部件進(jìn)行濕模態(tài)的分析，考慮口環(huán)和導(dǎo)軸承的支撐效應(yīng)，以及轉(zhuǎn)速、重力和陀螺效應(yīng)的影響，計(jì)算了不同口環(huán)間隙下轉(zhuǎn)子部件的臨界轉(zhuǎn)速。

由表3可知，考慮口環(huán)間隙后，轉(zhuǎn)子部件的臨界轉(zhuǎn)速變化主要體現(xiàn)在前兩階，一階臨界轉(zhuǎn)速增加了5.1%，二階臨界轉(zhuǎn)速增加了5.97%。隨著口環(huán)間隙的減小，轉(zhuǎn)子部件的臨界轉(zhuǎn)速有所提高，這是由于口環(huán)動(dòng)力特性系數(shù)增大，增強(qiáng)了轉(zhuǎn)子部件的剛度，提高了轉(zhuǎn)子部件的臨界轉(zhuǎn)速。由此可見，在高溫泵的臨界轉(zhuǎn)速計(jì)算時(shí)，考慮口環(huán)間隙的影響是必要的。

表3 不同口環(huán)間隙下轉(zhuǎn)子部件的臨界轉(zhuǎn)速Tab.3 Critical speed of rotor components with different clearances

不同口環(huán)間隙下高溫泵轉(zhuǎn)子的振型無明顯變化，本文只分析口環(huán)間隙為0.6 mm時(shí)轉(zhuǎn)子部件的前四階振型，從圖8中可以看出，各階轉(zhuǎn)子部件振型存在明顯的差異，第1階和第2階呈現(xiàn)以泵軸為中心的擺動(dòng)變形，且擺動(dòng)方向分別沿x、y軸互相垂直；第3階呈現(xiàn)以泵軸為中心的扭轉(zhuǎn)變形，此時(shí)變形量的分布呈中心對(duì)稱且葉輪邊緣處的變形量達(dá)到最大值；第4階振型呈現(xiàn)出以泵軸為中心的一彎擺動(dòng)變形，且葉輪在彎曲方向的邊緣處變形量達(dá)到最大。

圖8 口環(huán)間隙為0.6 mm時(shí)轉(zhuǎn)子部件前四階振型Fig.8 First to fourth order vibration modes of rotor components for clearance of 0.6 mm

3.4 轉(zhuǎn)子部件的諧響應(yīng)分析

由于機(jī)械加工誤差或高溫變形，會(huì)導(dǎo)致轉(zhuǎn)子部件運(yùn)行過程中產(chǎn)生不平衡力，在這種不平衡載荷作用下產(chǎn)生振動(dòng)響應(yīng)稱為不平衡諧響應(yīng)。在諧響應(yīng)求解設(shè)置過程中，設(shè)置頻率范圍為100 Hz，在葉輪上施加10 N的不平衡載荷，方向沿Y軸方向，葉輪出口輸出位移響應(yīng)數(shù)據(jù)，圖9示出不同口環(huán)間隙下葉輪在y方向的位移－頻率響應(yīng)曲線。

圖9 葉輪在Y方向的位移－頻率響應(yīng)曲線Fig.9 Displacement-frequency response curve of impeller in Y direction

如圖9所示，不考慮口環(huán)間隙時(shí)，葉輪的最大位移遠(yuǎn)大于考慮口環(huán)間隙時(shí)，且隨著間隙的減小葉輪的最大位移明顯減小，當(dāng)口環(huán)間隙減小至0.6 mm時(shí)，位移幅值降低至0.088 3 mm，較0.8 mm時(shí)候降低了52.6%，這說明口環(huán)間隙減小使口環(huán)動(dòng)力特性系數(shù)增大，降低了轉(zhuǎn)子部件的振動(dòng)位移；隨著口環(huán)間隙的減小，葉輪發(fā)生最大位移時(shí)所對(duì)應(yīng)的頻率有所增大，這是由于轉(zhuǎn)子部件的臨界轉(zhuǎn)速增大，而葉輪發(fā)生最大位移時(shí)的頻率為轉(zhuǎn)子部件臨界轉(zhuǎn)速所對(duì)應(yīng)的頻率，故轉(zhuǎn)子部件在臨界轉(zhuǎn)速所對(duì)應(yīng)的頻率附近易發(fā)生共振，為提高轉(zhuǎn)子運(yùn)行穩(wěn)定性，轉(zhuǎn)子部件的外部激勵(lì)應(yīng)避免該頻率區(qū)域。

4 結(jié)論

（1）口環(huán)間隙值從0.4 mm增到1.0 mm，模型泵的外特性有明顯變化，揚(yáng)程下降了5.95%，效率下降了6.26%，且口環(huán)間隙越大性能降低越明顯。

（2）采用流熱固耦合的方法，對(duì)轉(zhuǎn)子部件進(jìn)行靜力學(xué)分析，最大形變量位于葉輪流道中部蓋板外緣位置和口環(huán)處，口環(huán)間隙增大，口環(huán)泄漏量增大導(dǎo)致泵腔內(nèi)流體壓力減小，葉輪口環(huán)處的最大變形量減小。

（3）隨著口環(huán)間隙的減小，口環(huán)動(dòng)力特性系數(shù)增大，增強(qiáng)了轉(zhuǎn)子剛度從而提高了轉(zhuǎn)子部件的臨界轉(zhuǎn)速，但主要體現(xiàn)在前兩階，而不同口環(huán)間隙下轉(zhuǎn)子部件的振型無明顯變化。

（4）口環(huán)間隙減小時(shí)，葉輪在y方向的最大位移明顯減小，而發(fā)生最大位移時(shí)對(duì)應(yīng)的頻率有所增大，這是由于口環(huán)間隙減小提高了轉(zhuǎn)子部件臨界轉(zhuǎn)速，葉輪發(fā)生最大位移時(shí)的頻率為轉(zhuǎn)子部件臨界轉(zhuǎn)速所對(duì)應(yīng)的頻率。