程琦　路義萍　孫海鋒　艾麗昆　韓家德　謝增林

摘要：屏蔽電機重要部件的峰值溫度是影響屏蔽電機安全穩(wěn)定運行的關(guān)鍵因素。對內(nèi)置螺旋管的屏蔽電機而言，螺旋管的布置方式影響著屏蔽電機的溫度場。本文對某內(nèi)置螺旋管冷卻的屏蔽電機，分別將螺旋管布置成叉排及順排，并采用基于計算流體力學原理（CFD）的有限體積法對兩種不同布置方式下屏蔽電機的溫度場進行數(shù)值模擬和對比分析。計算結(jié)果表明，電機機殼和一、二次水的溫度分布呈現(xiàn)上部溫度低，下部溫度高的特點;雙層螺旋管叉排布置方式時，冷卻水的換熱量大于其順排布置方式下冷卻水的換熱量，且屏蔽電機的整體平均溫度高;為屏蔽電機的內(nèi)置換熱器布置方式的選擇提供參考。

關(guān)鍵詞：屏蔽電機;螺旋管布置方式;溫度場;有限體積法;數(shù)值模擬

DOI：10.15938/j.jhust.2022.04.002

中圖分類號： TK124

文獻標志碼： A

文章編號： 1007-2683（2022）04-0010-06

Effect of Spiral Tube Arrangement on Temperature Field

of Canned Motor for Nuclear Pump

CHENG Qi LU Yi-ping SUN Hai-feng AI Li-kun HAN Jia-de XIE Zeng-lin

（1.School of Electric and Engineering， Harbin University of Science and Technology， Harbin 150080，China;

2.Harbin Electric Power Equipment Co.， Ltd.， Harbin 150040， China）

Abstract：The maximum temperature of main components of a canned motor is the key factor affecting the safe and stable operation of the canned motor.?Especially for a canned motor with built-in spiral tube， the arrangement of spiral tube will have a great influence on the temperature field of the canned motor.?In this paper， the physical model of cooling canned motor with built-in spiral tube was established and its internal spiral tube was arranged into two ways of staggered and parallel.?Then， the temperature fields of the canned motor under two different arrangements were numerically simulated and contrastively analyzed by the finite volume method based on the CFD principle.?The temperature distribution of the motor casing and the primary and secondary water shows the characteristics of low temperature in the upper part and high temperature in the lower part.?The heat transfer of cooling water in the arrangement of staggered of double-layer spiral tubes is greather than that in the arrangement of parallel， and the overall average temperature of the canned motor is higher.?It provides a reference for the selection of the layout of the built-in heat exchanger arrangement of the canned motor.

Keywords：canned motor; spiral tube arrangement; thermal fields; finite volume method; numerical simulation

0引言

隨著“十四五”規(guī)劃和碳中和理念的不斷實施，核能發(fā)電作為一種清潔高效發(fā)電技術(shù)目前處于重要地位[1-3]。核主泵屏蔽電機作為核電站的重要驅(qū)動設(shè)備，直接影響電廠運行的安全性及可靠性，為避免核電站產(chǎn)生事故，在設(shè)計階段，進行核主泵屏蔽電機結(jié)構(gòu)變化的熱設(shè)計及其數(shù)值模擬研究，具有十分重要的意義[4-7]。

目前，國內(nèi)外學者對各種屏蔽電機溫度場的影響因素和螺旋管的流動傳熱分別進行了各種研究。丁樹業(yè)等[8]對額定工況下5500kW核主泵屏蔽電機的熱流場進行數(shù)值模擬，揭示了屏蔽電機內(nèi)的溫度分布規(guī)律。仲維濱等[9]建立了屏蔽電機模型，通過改變定子端腔表面發(fā)射率，發(fā)現(xiàn)定子端腔表面發(fā)射率的提高使得定子繞組峰值溫度降低。胡家祺[10]通過對屏蔽電機的熱流場進行數(shù)值模擬，研究了核主泵冷卻劑溫度及連接部件對屏蔽電機溫度場的影響。曹力等[11]對電機定子屏蔽套的渦流損耗和溫升進行分析，同時對屏蔽套損耗的敏感性因素進行研究。LU Y P等[12]研究了不同的水摩擦損耗計算方法對屏蔽電機溫度場模擬結(jié)果的影響。畢曉舜等[13]分析了冷卻水流量、水道寬度、水道數(shù)及冷卻水溫對電機升溫的影響，為電機冷卻系統(tǒng)的選擇提供參考。Bornschlegell A S等[14]對大功率屏蔽電機進行了熱優(yōu)化問題的研究。Hu X等[15]研究了定轉(zhuǎn)子之間的氣隙厚度對電機升溫和電磁性能的影響，確定了氣隙厚度的最佳取值。黃云云等[16]對三種截面的螺旋管的流動及傳熱進行模擬，確定在傳熱性能方面，圓管最佳。Wang Y等[17]對帶有反向回路和線圈插入件的螺旋管的熱性能進行研究。岳清雯等[18]采用熱流耦合的數(shù)值模擬方法分析了殼側(cè)流體的流量和溫度變化對水平螺旋管換熱器熱流場的影響。董珊珊[19]對某屏蔽電機水冷夾套內(nèi)螺旋管的強化換熱進行研究。甘劉意等[20]通過對兩種螺旋板式換熱器進行數(shù)值模擬，分析對比得出雷諾數(shù)對換熱器傳熱影響因子、阻力影響因子以及綜合性能評價系統(tǒng)的變化規(guī)律。Nabeel J等[21]對單層螺旋管換熱器進行數(shù)值模擬，并對螺旋管換熱器沿長度和徑向方向的溫度進行測定，對于螺旋管換熱器的設(shè)計及溫度場的研究具有指導意義。

綜上所述，大多數(shù)學者對屏蔽電機和螺旋管換熱器內(nèi)部溫度場及其影響因素的研究主要是關(guān)于自身因素，關(guān)于將螺旋管換熱器和屏蔽電機進行整體研究未見文獻報道。本文中，所涉及的核主泵位于屏蔽電機下方，電機機殼外側(cè)布置雙層螺旋管作為一次水的換熱器，用于小型研究堆，是近期核電的重要發(fā)展方向之一。

1模型建立

1.1物理模型

本屏蔽電機轉(zhuǎn)速由工廠提供，為2922r/min，通過一、二次冷卻水循環(huán)流動來傳遞熱量進而冷卻降溫。一次冷卻水循環(huán)包括一條主循環(huán)回路流動和兩條支路循環(huán)流動：主循環(huán)冷卻水從定子端蓋一次冷卻水入口進入沿軸心孔向下流動，流經(jīng)輔葉輪、定轉(zhuǎn)子屏蔽套間隙、轉(zhuǎn)子下部空腔和下導軸承定子封頭徑向孔匯集到匯流腔，最終通過內(nèi)置螺旋管進行冷卻，從一次水出口流出進入一次水入口完成主循環(huán)。兩條支路分別為繞經(jīng)上導軸承和繞經(jīng)下導軸承及下飛輪后匯合到主循環(huán)回路。二次冷卻水循環(huán)上部端腔二次水入口進入，流經(jīng)循環(huán)槽道到達下端腔，最終從下端腔外罩板處的二次水出口排出，見圖1。

本屏蔽電機沿周向?qū)ΨQ分布，考慮物理模型、網(wǎng)格劃分和數(shù)值計算的高效性，故選取整機周向四分之一建立物理模型，Y軸與電機轉(zhuǎn)軸的中心線重合，X軸和Z軸沿半徑方向，見圖2。

1.2網(wǎng)格劃分

本文運用Gambit軟件對屏蔽電機物理模型進行網(wǎng)格劃分，見圖3。為了保證網(wǎng)格質(zhì)量，采用結(jié)構(gòu)性網(wǎng)格和非結(jié)構(gòu)性網(wǎng)格混合使用的方法。Gambit軟件自帶網(wǎng)格評價標準，通過0～1的數(shù)值大小進行判斷，數(shù)值計算的網(wǎng)格質(zhì)量要求Equisize skew小于0.97，最終對全部網(wǎng)格質(zhì)量進行檢查，Equisize skew最差為0.94，多數(shù)為0.1，同時考慮到一次冷卻水路存在間隙流，并為了滿足數(shù)值計算的二階迎風格式，對間隙流處網(wǎng)格進行加密，網(wǎng)格數(shù)目為4。最后，為分析螺旋管布置方式對屏蔽電機溫度場的影響，將螺旋管由叉排布置改為順排布置，為清晰可見，圖3中給出了整機模型和重要位置定轉(zhuǎn)子間隙、定子繞組、內(nèi)外套板中間螺旋管采用叉排布置及順排布置時的局部放大網(wǎng)格圖。

2數(shù)學模型及求解條件

2.1基本假設(shè)

1）屏蔽電機內(nèi)流體的馬赫數(shù)（Ma）小于1，故可認為不可壓縮流體。

2）流體的物性參數(shù)如密度、比熱容等均為常數(shù)。

3）設(shè)每個部件熱源按照體積均勻分布。

4）相鄰部件之間無接觸熱阻存在。

5）一次冷卻水和二次冷卻水的雷諾數(shù)較大（Re>2300），流體的流動處于湍流狀態(tài)。

2.2控制方程

根據(jù)計算流體力學理論，對屏蔽電機流體場、溫度場進行數(shù)值模擬時，模型中的介質(zhì)滿足質(zhì)量、動量和能量守恒方程，控制方程的通用形式為

2.3計算條件

影響仿真計算結(jié)果準確性的重要因素是計算條件。根據(jù)實際運行情況，屏蔽電機熱態(tài)額定工況下的計算條件;

1）一、二次冷卻水入口速度分別為1.67m/s、1.471m/s，入口溫度分別為48℃、40℃，出口表壓力為0Pa。

2）機殼表面為同時考慮自然對流與輻射的混合邊界條件，對流換熱系數(shù)為10W/（m²·k），主泵下方冷卻劑溫度為285.1℃。

3）電機內(nèi)部具有多種損耗，額定工況下，部件損耗考慮了定轉(zhuǎn)子屏蔽套、定轉(zhuǎn)子鐵心、定子繞組、轉(zhuǎn)子銅條、轉(zhuǎn)子端環(huán)、錐形環(huán)和定子齒壓板，數(shù)值較大的熱源強度見表1。

3結(jié)果與分析

3.1電機溫度分布特征

為保證屏蔽電機的安全穩(wěn)定運行，對屏蔽電機的溫度場進行數(shù)值模擬。為了便于對屏蔽電機整體溫度進行分析，同時考慮到在屏蔽電機定子端腔處繞組溫度較高、鐵心中部處繞組溫度較低[2]，故在屏蔽電機Y=0.8m、Y=1.3m和整機45°處做剖面，并給出溫度分布云圖如圖4所示。

由圖4可見，該電機沿軸向溫度總體分布呈現(xiàn)從上到下溫度逐漸升高的趨勢，定子繞組端腔中溫度較高，且下端腔的峰值溫度高于上端部。一方面原因是冷卻水沿流動方向逐漸升溫，冷卻能力逐漸降低，另一方面，屏蔽電機下部存在主泵溫度為285.1℃的高水溫，高溫水向四周進行熱量傳遞。由于定子繞組端部被氮氣腔包裹，氮氣與繞組端部、氮氣與機殼發(fā)生自然對流換熱，而自然對流換熱系數(shù)遠遠小于定子繞組中段與冷卻水發(fā)生的強制對流換熱系數(shù)，故定子繞組端部溫度高，中間段溫度低;屏蔽電機沿徑向的溫度分布呈現(xiàn)出以定子為中心向半徑增大及半徑減小的方向逐漸降低的趨勢，這是因為定轉(zhuǎn)子屏蔽套中的冷卻水和機殼外部的冷卻水溫度較低，熱量總是由高溫傳向低溫。

3.2叉排與順排布置溫度對比

3.2.1定轉(zhuǎn)子部分溫度分布對比

雙層螺旋管換熱器的布置方式對屏蔽電機內(nèi)部部件的溫度分布具有影響，叉排布置方式下定子繞組的峰值溫度為169.5℃，比順排布置方式下定子繞組的峰值溫度高0.3℃。不同布置方式下定子鐵心齒部及軛部溫度也有差異，均為峰值溫度叉排時高，但是高出的數(shù)值有限，都不超過1℃。這是由于叉排布置方式下一次冷卻水的換熱充分，導致叉排布置時二次冷卻水的溫度高，進而導致傳熱溫差小，而轉(zhuǎn)子鐵心和轉(zhuǎn)子銅條與機殼換熱器相距較遠，并且與一次水相隔，故不同布置方式對轉(zhuǎn)子部分的峰值溫度影響不明顯。

3.2.2電機內(nèi)部冷卻水溫度分布對比

屏蔽電機內(nèi)部具有復雜的一次冷卻水回路，一次冷卻水在電機內(nèi)部循環(huán)流動，并且位于水路的不同位置具有不同的溫度特征，為方便監(jiān)測屏蔽電機不同位置的冷卻水溫度，在劃分網(wǎng)格時，將一次冷卻水根據(jù)位置進行分組，以便于后期監(jiān)測。不同布置方式下，電機上部一次冷卻水溫度差異很小，表2給出了具有差異下部位置的一、二次冷卻水溫度數(shù)值。從表中可以看出螺旋管叉排布置時，二次水出口、下導軸承、匯流腔和下飛輪上部的一次冷卻水高于順排布置方式下相同位置的溫度，說明叉排布置換熱效果好;而下飛輪周圍及下飛輪底部溫度反而低于順排布置方式下相同部位的溫度;順排布置時的匯流腔水溫低，則更易于實現(xiàn)二次水斷供故障5min之內(nèi)，匯流腔水溫不超過95℃這一指標。

3.2.3機殼溫度分布對比

由于機殼與上下端部存在焊縫，是熱應力計算重點監(jiān)測位置。為比較雙層螺旋管不同布置方式對屏蔽電機機殼溫度分布的影響，在屏蔽電機機殼徑向1/2處做一條軸向溫度采樣線，圖5為該采樣線溫度分布圖。對比可知，雙層螺旋管不同布置方式對屏蔽電機的總體溫度分布趨勢無影響，都是上部溫度低，下部溫度高，中部鐵心段相同位置處，基本呈線性變化，叉排時機殼溫度高于順排時，叉排與順排布置方式下機殼的最高溫度分別為60.9℃、58.8℃，兩者相差2.1℃。機殼溫度數(shù)值可為后續(xù)熱應力計算提供基礎(chǔ)數(shù)據(jù)。

3.2.4換熱器內(nèi)冷卻水溫度分布及換熱量對比

屏蔽電機的冷卻主要通過內(nèi)外冷卻回路，一次冷卻水主要冷卻屏蔽電機內(nèi)部部件并吸收熱量，二次冷卻水主要冷卻一次冷卻水和機殼，為了說明雙層螺旋管換熱器的布置方式對螺旋管內(nèi)一次冷卻水的冷卻和螺旋槽內(nèi)二次冷卻水的升溫影響，同時，考慮到雙層螺旋管換熱器的布置方式的區(qū)別主要在于外圈螺旋管，故分別選取雙層螺旋管/槽的內(nèi)圈螺旋管/槽的45°截面中心處及螺旋管的入出口直管截面中心、上端腔及下端腔中的螺旋管/槽45°截面中心處做采樣監(jiān)測點，比較雙層螺旋管順/叉排布置方式下螺旋管/槽內(nèi)冷卻的水的溫度分布，見圖6。

由圖6可知，對于螺旋管內(nèi)的一次冷卻水來說，雙層螺旋管的布置方式對螺旋管內(nèi)一次冷卻水的溫度分布趨勢無影響，都是沿著向上流動方向溫度逐漸降低，最高溫度位于換熱器下部端腔螺旋管內(nèi)水的入口處，但是叉/順排布置方式下螺旋管內(nèi)冷卻水的最高溫度數(shù)值不同，分別為69.7和65.7℃;對于二次冷卻水來說，都是沿著向下流動方向溫度逐漸升高，二次冷卻水的峰值溫度叉排比順排布置方式下高出3.7℃，說明換熱效果更好。

為了進一步對比雙層螺旋管不同布置方式對換熱的影響，對不同布置方式下一次冷卻水的放熱量Q₁和二次冷卻水的吸熱量Q₂進行計算并分析，見式3。

由于本臺屏蔽電機符合熱平衡要求，故一、二次冷卻水的放、吸熱量等于各自的傳熱熱量。表3為兩種布置方式時，換熱器中一、二次冷卻水的放、吸熱量。由表3可知，雙層螺旋管叉排布置方式下一、二次冷卻水的放、吸熱量均大于順排布置方式下的數(shù)值，二次冷卻水的吸熱量大于一次冷卻水的放熱量，原因是二次冷卻水溫度最低，不僅吸收一次冷卻水放出的熱量，同時，還吸收屏蔽電機經(jīng)由機殼傳遞出的熱量及下方285.1℃的高溫冷卻劑向上方傳遞的熱量。

4結(jié)論

1）本文研究的主泵位于下方的核主泵屏蔽電機，溫度分布呈現(xiàn)由上到下溫度逐漸升高的特點。

2）雙層螺旋管叉排布置時，屏蔽電機內(nèi)部重要監(jiān)測部件及冷卻水的峰值溫度高于雙層螺旋管順排布置時內(nèi)部重要部件及冷卻水的峰值溫度。

3）雙層螺旋管叉排布置方式時一、二次冷卻水的放、吸熱量大于其順排布置方式下一、二次冷卻水的放、吸熱量。

參 考 文 獻：

[1]鄧清華， 胡樂豪， 李軍， 等. 大型發(fā)電技術(shù)發(fā)展現(xiàn)狀及趨勢[J].熱力透平，2019，48（3）：175.DENG Qinghua， HU Lehao， LI Jun， et al. State-of-Art and Tendency on Technologies of Large Electric Power Generation[J]. Thermal Tbrbine，2019，48（3）：175.

[2]索文旭， 尉盛強， 付偉達， 等. 一核主泵屏蔽電機熱流場特征分析[J].電站系統(tǒng)工程，2021，37（5）：71.SUO Wenxu， WEI Shengqiang， FU Weida， et al. Analysis of Thermal and Flow Field Characteristics for Shielded Motor of Primary Pump[J]. Power System Engineering，2021，37（5）：71.

[3]HE Y L， TANG S Z， TAO W Q， et al. A General and Rapid Method for Performance Evaluation of Enhanced Heat Transfer Techniques[J]. International Journal of Heat and Mass Transfer，2019，145：118780.

[4]李藏雪， 任智達， 路義萍，等. 一種屏蔽電動機的三維溫度場及影響因素研究[J].電機與控制學報，2015，19（8）：15.LI Cangxue， REN Zhida， LU Yiping， et al. Influence Factors Research on Three Dimensional Temperature Field of the Canned Motor[J]. Electric Machines and Control，2015，19（8）：15.

[5]龍云， 袁壽其， 朱榮生，等. 核主泵內(nèi)部流動研究現(xiàn)狀與技術(shù)發(fā)展綜述[J].排灌機械工程學報，2020，38（11）：1081.LONG Yun， YUAN Shouqi， ZHU Rongsheng， et al. Review on Research Status of Internal Flow and Technological Development of Reactor Coolant Pump[J]. Jonrnal of Drainage and Irrigation Machinery Engineering，2020，38（11）：1081.

[6]LI Z， LI M， ZHANG J， et al. Seismic Analysis of Nuclear Power Plant Canned Motor Pump Unit Based on Integral Calculation Method[J]. Stave Bníobzor-Civil Engineering Journal，2016，25（4）：1.

[7]索文旭，艾麗昆，胡家祺，等.屏蔽電機機殼內(nèi)螺旋槽道截面形狀對流動及傳熱影響[J].大電機技術(shù)，2021（5）：59.SUO Wenxu， AI Likun， HU Jiaqi， et al. Effect of the Cross Section Shape of Spiral Groove in the Canned Motor Case on Water Flow and Heat Transfer[J]. Large Electric Machine and Hydraulic Turbine， 2021（5）：59.

[8]丁樹業(yè)， 孟繁東， 葛云中. 核主泵屏蔽電機溫度場研究[J].中國電機工程學報，2012，32（36）：149.DING Shuye， MENG Fandong， GE Yunzhong. Temperature Field Investigation of Canned Primary Pump Motors in Nuclear Power Stations[J]. Proceedings of the CSEE，2012，32（36）：149.

[9]仲維濱， 劉延浩， 周超，等. 屏蔽電機氮氣腔壁表面吸收率對溫度場的影響[J].大電機技術(shù)，2018（6）：62.ZHONG Weibin， LIU Yanhao， ZHOU Chao， et al. Influence of Wall Absorptivity in Nitrogen Cavity on the Temperature Field of Canned Motor[J]. Large Electric Machine and Hydraulic Tnrbine，2018（6）：62.

[10]胡家祺.連接部件對某屏蔽電機溫度場影響的數(shù)值模擬[D].哈爾濱：哈爾濱理工大學，2021.

[11]曹力， 胡巖， 卓亮. 高速永磁屏蔽電機損耗分析與溫升研究[J].微電機，2021，54（4）：11.CAO Li， HU Yan， ZHUO Liang. Study on Loss Analysis and Temperature Rise of High-speed Permanent Magnet Shielded Motor [J]. Micromotors，2021，54（4）：11.

[12]LU Y P， MUSTAFA A， REHAN M A， et al. The Effects of Water Friction Loss Calculation on the Thermal Field of the Canned Motor[J]. Processes，2019，7（5）：256.

[13]畢曉舜， 安躍軍， 楊維國， 等. 真空干泵用屏蔽電機溫度場仿真分析與試驗[J].電機與控制應用，2020，47（5）：64.BI Xiaoshun， AN Yuejun， YANG Weiguo， et al. Simulation and Experiment of Temperature Field of Canned Motor for Vacuum Dry Pump[J]. Electric Mchines and Control Application， 2020，47（5）：64.

[14]BORNSCHLEGELL A S， PELLE J， HARMAND S， el al. Thermal Optimization of a High-Power Salient-Pole Electrical Machine[J]. IEEE Transactions on Industrial Electronics，2013，60（5）：1734.

[15]HU X， LI Y， LUO L. The Influence of Air Gap Thickness Between the Stator and Rotor on Nuclear Main Pump[J]. Energy Procedia，2017，142：259.

[16]黃云云，黃永曙，王為旺.不同截面形狀螺旋管流動與傳熱的數(shù)值研究[J].計算機與應用化學，2016，33（12）：1261.HUANG Yunyun， HUANG Yongshu， WANG Weiwang. Numerical Research on Flow and Heat Transfer in Helically Coiled Tubes with Different Cross-sections[J]. Computers and Applied Chemistry，2016，33（12）：1261.

[17]WANG Y， ALVARADO J L， TERRELL W. Thermal Performance of Helical Coils with Reversed Loops and Wire Coil Inserts[J]. International Journal of Heat and Mass Transfer，2020，146：118723.

[18]岳清雯， 賴喜德， 陳小明， 等. 水平螺旋管式換熱器的流熱耦合傳熱特性研究[J].熱能動力工程，2021，36（4）：118.YUE Qingwen， LAI Xide， CHENG Xiaoming， et al. Study on Thermal-fluid Coupling Heat Transfer Characteristics of Horizontal Spirally Coiled Tubes Heat Exchanger[J]. Jpurnal of Engineering for Thermal Energy and Power，2021，36（4）：118.

[19]董姍姍.屏蔽電機水冷夾套內(nèi)螺旋管強化換熱研究[D].哈爾濱：哈爾濱理工大學，2020.

[20]甘劉意， 陸怡， 查涵清， 等. 新型螺旋板式換熱器及其傳熱特性研究[J].流體機械，2021，49（6）：36.GAN Liuyi， LU Yi， ZHA Hanqing， et al. Reaearch on New Type Spiral Plate Heat Exchanger and Its Heat Transfer Characteristics[J]. Fluid Machinery，2021，49（6）：36.

[21]JAFFAR N， INAYAT M，WARDAG A. Design a Helical Coil Heat Exchanger Via CFD Simulations[C]// 2017 14th International Bhurban Conference on Applied Sciences and Technology （IBCAST） IEEE，2017：594.

（編輯：溫澤宇）