于喜偉 孟大偉 李瓊

摘 要：以一臺6.5 MW正壓型防爆電機為例，首次將多組分流體理論應用于正壓型防爆電機，采用了有限元法對電機的安全性能進行了仿真分析，分別完成了“模擬爆炸環(huán)境瞬態(tài)仿真計算”和“吹掃瞬態(tài)仿真計算”，得到了吹掃過程中爆炸性氣體及清潔空氣各組成成分的濃度變化云圖，通過各部位的濃度云圖能夠得知電機內(nèi)部是否存在吹掃死角，增強了安全性能。仿真結(jié)果經(jīng)過防爆試驗驗證，計算數(shù)據(jù)滿足精度要求。進一步優(yōu)化了進氣管道的布置方案，通過合理的調(diào)整進氣管道開孔的位置、數(shù)量及大小提高了吹掃效率，節(jié)約了試驗成本，同時也可以使電機快速啟動，更好的服務工況。

關(guān)鍵詞：正壓型防爆;多組分;流體場;瞬態(tài)分析

中圖分類號：TM 315

文獻標志碼：A

文章編號：1007-449X（2019）04-0049-07

0 引 言

防爆電機廣泛應用于石油、化工等行業(yè)，按照防爆類型主要分為隔爆型、增安型、無火型及正壓型，隔爆型電機是通過加厚機殼厚度及減小防爆間隙來實現(xiàn)防爆，這使得材料利用率很低[1-3]，且由于需要通過較高的水壓試制，因此難以在大型電機里得到應用[4-6];增安型和無火花型電機是通過減小定轉(zhuǎn)子電暈來實現(xiàn)防爆，電機為了滿足定子繞組脈沖試驗，繞組設計極為復雜且成本很高，在行業(yè)里應用較少[7];正壓型電機起動前向電機內(nèi)部吹入干燥的空氣或惰性氣體（氦氣、二氧化碳、氮氣等），用來置換電機內(nèi)部可能存在的爆炸性氣體，使電機內(nèi)部爆炸性氣體的濃度降到爆炸極限以下，這一過程稱為“吹掃”過程。在電機運行過程中也需要向電機內(nèi)部不斷通入少量保護氣體，用以保持電動機內(nèi)部壓力大于外部壓力，防止外部爆炸性氣體再次進入電機內(nèi)部，從而達到防爆的目的。正壓型電機克服了其它防爆電機的主要缺點，近年來，隨著正壓保護裝置逐漸完善及國際標準IEC60079-2第6版正式發(fā)布，各大電機公司陸續(xù)將正壓型防爆電機推向了市場。

對正壓型防爆電機，目前主要采用試驗的方式來保證其安全性能[8-10]，如此不僅成本高，且覆蓋性不好，必須每臺電機都要進行實驗同時缺乏理論基礎，難以為今后的設計提供參考依據(jù)。由于正壓型防爆電機的應用歷史較短，且分析的區(qū)域是由空氣和爆炸性氣體構(gòu)成的多組分流體場，因而還沒有開展其安全性能的模擬仿真研究。目前對多組分流體場的分析主要集中在石油化工領域，不但考慮了多組分的密度和粘度等自身特性影響[11-14]，有些還分析了多組分之間的相變、溶解、化學反應等因素[15-16]，在環(huán)境污染領域也偶有建樹[17]。而在電機領域的應用多為單一物質(zhì)的流體場分析，文獻[18-19]對電機的定子流體場進行了分析，文獻[20]計算了電機轉(zhuǎn)子流體場，文獻[21]對無刷勵磁機通風冷卻進行了模擬，文獻[22]對電機的冷卻器的流體場進行了優(yōu)化，文獻[23]對電機的全域流體進行了預測，而對電機進行多組分流體場分析的文獻尚未發(fā)現(xiàn)。

本文以一臺6.5 MW正壓型防爆電機為例，結(jié)合實際運行情況給出合理的基本假設與邊界條件，并根據(jù)多組分流體力學原理，采用有限元法對電機的安全性能進行了仿真分析，得到了吹掃過程中不同時間下的各組分氣體的濃度分布情況、防爆試驗氣體消耗量及總計吹掃時間;同時與防爆試驗進行了對比，為今后正壓型防爆電機安全性能的仿真分析提供了理論基礎;最后對進氣管道的布置進行了改進，為后續(xù)類似產(chǎn)品設計提供參考。

1 物理模型

電動機的基本參數(shù)如表1所示。

樣機實物如圖1所示，圖2為根據(jù)電機實物確立的求解模型，因分析吹掃過程只與電機內(nèi)部氣體流體有關(guān)，所以建模時略去了強迫風機、軸承、接線盒等部件，而對電機內(nèi)部的定子、轉(zhuǎn)子和冷卻管等部件都進行了詳細繪制。

在圖2中定義的橫向截面和縱向截面及點A、B、C、D、E是用于仿真分析說明，在軸向截面上與點A、B、C、D、E相對應點F、G、H、I、J，由于空間有限在此圖中未示意。

2 有限元數(shù)值仿真分析

2.1 網(wǎng)格劃分

采用結(jié)構(gòu)網(wǎng)格對模型進行劃分，由于電機整體尺寸很大，內(nèi)部又具有很多狹小空間，為了得到較好的網(wǎng)格質(zhì)量需要將單元尺寸設置很小，共計單元數(shù)1993萬，節(jié)點數(shù)352萬。圖3為網(wǎng)格剖分圖。

2.2 計算模型建立

2.2.1 數(shù)學描述

計算區(qū)域內(nèi)由多種成分組成，流體的物性參數(shù)由每個組分流體的物性及其質(zhì)量分率決定，各組分質(zhì)量分率之和等于 1 ，即：

采用Species Transport模型，求解時利用質(zhì)量守恒方程對流體流動進行控制，如式（4）所示;且由于流體粘性的作用，流體流動滿足動量守恒方程，如式（5）所示;吹掃過程的流動狀態(tài)由雷諾數(shù)進行判斷，如式（6）所示，經(jīng)過計算雷諾數(shù)為90 855，遠大于4 000，流動狀態(tài)屬于湍流，采用了標準的湍流k-ε模型，如式（7）所示。

2.2.2 基本假設

1）不考慮多組分之間化學反應;

2）氣體在整個運動過程中溫度不變;

3）電機內(nèi)氣體處于常壓下，忽略流體域內(nèi)流體浮力和重力的影響;

4）流體流速遠小于聲速，即把電機內(nèi)流體作為不可壓縮流體處理。

2.2.3 邊界條件

2.2.3.1 模擬爆炸環(huán)境瞬態(tài)仿真計算

1）電機安裝場所為Zone 2（IIC），主要爆炸氣體為H2和C2H2，模擬爆炸氣體采用He（防爆試驗時常使用He模擬爆炸氣體，因He密度與氫氣相近，且為惰性氣體）。

2）電機內(nèi)部初始為空氣（由21%的O2和79%的N2組成）。

3）從電機入口通入濃度100%的He，采用速度入口，速度為30 m/s。

4）采用壓力出口。

5）計算的結(jié)束條件為：He在各個取樣點的濃度不低于70%（按照防爆標準要求）。

2.2.3.2 吹掃瞬態(tài)仿真計算

1）結(jié)合真實試驗過程，“模擬爆炸環(huán)境”后將對電機進行“吹掃”試驗，因此以“模擬爆炸環(huán)境”的結(jié)束狀態(tài)作為吹掃瞬態(tài)計算的初始狀態(tài)。

2）向電機內(nèi)部吹入空氣（由21%的O2和79%的N2組成），采用速度入口，速度為30 m/s。

3）采用壓力出口。

4）計算的結(jié)束條件為：He在各個取樣點的濃度不高于1%（按照防爆標準要求）。

2.3 計算結(jié)果

本文的所有分析均為瞬態(tài)過程，為了獲取不同時刻的各組分的濃度情況，每隔240 s進行一次數(shù)據(jù)存儲，圖4為模擬爆炸環(huán)境瞬態(tài)計算的t=1 680 s時He濃度分布云圖，當t=1 680 s時計算域內(nèi)所有位置He濃度都達到70%以上，完成了模擬爆炸環(huán)境，所消耗的He總量為Q=71 m3。根據(jù)保存下來的不同時刻濃度數(shù)據(jù)，可以獲得電機不同取樣點的各組分濃度隨著時間變化曲線，如圖5所示。

圖6為吹掃瞬態(tài)計算t=5 040 s時的He濃度分布云圖，圖中所有位置He濃度都降低到1%以下，且無吹掃死角，吹掃時間為3 360 s，最小換氣量142 m3。圖7為吹掃瞬態(tài)計算時電機不同取樣點的各組分濃度隨著時間變化曲線。

3 防爆試驗

樣機裝配后進行了防爆試驗，圖8（a）為試驗時在電機外表面設置濃度測試點，機座上有6處，冷卻器上有4處，圖8（b）是在通過儀器進行各點濃度測量。進行爆炸環(huán)境試驗：當所有測試點He濃度都達到70%以上時用新鮮空氣進行吹掃試驗：當所有測試點He濃度都達到1%以下時，吹掃時間為58 min，最小換氣量150 m3。仿真計算結(jié)果與試驗數(shù)據(jù)相比，吹掃時間與最小換氣量的相對誤差分別為5.33%和3.44%，滿足工程精度要求。

4 改進

由仿真結(jié)果可以看出，在吹掃過程中圖7的C點He濃度降低最慢，圖6縱向截面圖的中部濃度降低較慢（即轉(zhuǎn)子內(nèi)部的濃度降低較慢）。濃度降低速度主要與進氣管道布置有關(guān)，因此對進氣管道進行了改進。改進前后的進氣管道如圖9所示，改進后A端比B端的開孔略大，可以提高C點側(cè)的濃度降低速度，改為兩排進氣使中部的進氣量增多，提高中部濃度降低速度。

圖10為改進進氣管道后模擬爆炸環(huán)境瞬態(tài)計算t=1 260 s時的He濃度分布云圖，所有位置He濃度都達到70%以上，所消耗的He總量Q=53 m3。與改進前相比，He用量降低了25%。圖11為改進進氣管道后模擬爆炸環(huán)境瞬態(tài)計算時各組分濃度隨著時間變化曲線。

圖12為改進進氣管道后吹掃瞬態(tài)計算t=3 780 s時的He濃度分布云圖，圖中所有位置He濃度都降低到1%以下，吹掃時間為2 520 s，較改進前減少了25%，極大提高了吹掃效率。同時，轉(zhuǎn)子內(nèi)部濃度降低速度已經(jīng)有所提高。圖13為改進進氣管道后吹掃瞬態(tài)計算時各組分濃度隨著時間變化曲線，從曲線中可以看出C點側(cè)的濃度降低速度已經(jīng)明顯提高。

5 結(jié) 論

本文對通過對6.5 MW正壓型防爆電機安全性能研究，可以得到如下結(jié)論：

1）仿真結(jié)果與防爆試驗數(shù)據(jù)吻合，說明采用多組分流體場有限元法對電機安全性能進行仿真是可行的。

2）通過該仿真方式能夠得到任一時間任一位置的氣體濃度分布情況，可以得知電機內(nèi)部是否存在吹掃死角，增強了電機安全性。

3）合理布置進氣管道，能夠節(jié)約試驗成本，提高吹掃效率，也加速了電機啟動。

參 考 文 獻：

[1]杜愛娟.礦用隔爆型變頻一體機水冷結(jié)構(gòu)設計及散熱研究[D].濟南：山東大學，2015.

[2]何惠明.隔爆變頻一體式電機關(guān)鍵技術(shù)研究[D].沈陽：沈陽工業(yè)大學，2012.

[3]郭儉.隔爆電機三維渦流場分析與計算[D].沈陽：沈陽工業(yè)大學，2009.

[4]DING Jie，ZHANG Ping.Coupled 3D fluid field & thermal field calculation of mine-used explosion-proof integrative variable-frequency motor[J].Electric Machines and Control，2015，19（7）：27.

[5]LIU Chenyang，ZHENG Xiaoquan，YIN Mo.The design of H level thermal-conductivity composite insulation structure for explosion-proof motor with high efficiency and low voltage[C]//2015 IEEE 11th International Conference on the Properties and Applications of Dielectric Materials，July 19-22，2015，Sydney，Australia. 2015：624-627.

[6]WANG Yu，HE Huiming，WANG Jiayin，et al.Research of novel water cooling jacket for explosion-proof motor[C]//2013 International Conference on Electrical Machines and Systems，Oct. 26-29，2013，Busan，South Korea. 2013：691-694.

[7]朱孟華.高壓增安型電機繞組防電暈試驗研究[J].防爆電機，2009，44（150）：44.

ZHU Menghua.Study on the anti corona test of high voltage increased safety motor windings[J].Explosion Proof Electric Machine，2009，44（150）：44.

[8]應一，魏桂梅.正壓型防爆電機的相關(guān)試驗[J].防爆電機，2011，46（159）：31.

YING Yi，WEI Guimei.Brief discussion on the related test of pressurized explosion-proof motor[J].Explosion Proof Electric Machine，2011，46（159）：31.

[9]許春家.正壓型防爆電機的防爆原理與設計[J].防爆電機，2008，43（143）：8.

XU Chunjia.Explosion proof principle and design of pressurized explosion proof motor[J].Explosion Proof Electric Machine，2008，43（143）：8.

[10]張文.15 MW正壓型電機的設計與研制[D].上海：上海交通大學，2014.

[11]邸元，吳玉樹，鞠斌山.裂縫性油藏二氧化碳驅(qū)的多組分數(shù)值計算模型[J].西南石油大學學報，2015，37（2）：93.

DI Yuan，WU Yushu，JU Binshan.Model formulation for compositional numerical simulation of CO2 flooding in fractured reservoirs[J].Journal of Southwest Petroleum University，2015，37（2）：93.

[12]李軼.多相流測量技術(shù)在海洋油氣開采中的應用與前景[J].清華大學學報，2014，54（1）：88.

LI Yi.Application and perspective of multiphase flow metering technologies for ocean oil and gas exploitation [J].Journal of Tsinghua University，2014，54（1）：88.

[13]何淼，柳貢慧，李軍，等.多相流全瞬態(tài)溫度壓力場耦合模型求解及分析[J].石油鉆探技術(shù)，2015，43（2）：25.

HE Miao，LI Gonghui，LI Jun，et al.Solution and analysis of fully transient temperature and pressure coupling model for multiphase flow[J].Petroleum Drilling Techniques，2015，43（2）：25.

[14]鄧慧，白焰.用擴散層模型與多相流方法對空冷凝汽器逆流管內(nèi)凝結(jié)與流動的數(shù)值模擬[J].中國電機工程學報，2017，37（3）：826.

DENG Hui，BAI Yan.Numerical simulation for flow and condensation in air cooled condenser reflux tube with diffusion layer model and multiphase-flow approach[J].Proceedings of the CSEE，2017，37（3）：826.

[15]付宇航，趙述芳，王文坦，等.多相/多組分LBM模型及其在微流體領域的應用[J].化工學報，2014，65（7）：2535.

FU Yuhang，ZHAO Shufang，WANG Wentan，et al.Application of lattice Boltzmann method for simulation of multiphase/multicomponent flow in microfluidics [J].CIESC Journal，2014，65（7）：2535.

[16]陳益峰，周創(chuàng)兵，童富果，等.多相流傳輸THM全耦合數(shù)值模型及程序驗證[J].巖石力學與工程學報，2009，28（4）：649.

CHEN Yifeng，ZHOU Chuangbing，TONG Fuguo，et al.A numerical model for fully coupled THM processes with multiphase flow and code validation[J].Chinese Journal of Rock Mechanics and Engineering，2009，28（4）：649.

[17]黃璐，趙成剛.基于復合混合物理論的熱-水力-力學污染物輸運模型[J].北京交通大學學報，2009，33（3）：93.

HUANG Lu，ZHAO Chenggang.A thermo hydro mechanical contaminant transport model based on hybrid mixture theory [J].Journal of Beijing Jiaotong University，2009，33（3）：93.

[18]韓家德，杜鵬，路義萍.凸極電機定子風路變化對熱流場影響[J].電機與控制學報，2016，20（12）：59.

HAN Jiade，DU Peng，LU Yiping.Effect of the stator ventilation ducts changes on thermal and flow field of salient synchronous motor [J].Electric Machines and Control，2016，20（12）：59.

[19]溫嘉斌，侯健，于喜偉.定子通風槽鋼對中型高壓電機內(nèi)溫度場的影響[J].電機與控制學報，2016，20（8）：40.

WEN Jiabin，HOU Jian，YU Xiwei.Influence of stator ventilation channel on the temperature field in the middle-size high voltage motor [J].Electric Machines and Control，2016，20（8）：40.

[20]李偉力，李丹，李金陽.水輪發(fā)電機轉(zhuǎn)子通風道堵塞時的耦合場分析[J].電機與控制學報，2017，21（1）：30.

LI Weili，LI Dan，LI Jinyang.Analysis of hydro generator rotor coupling field on rotor ventilation duct blocked[J].Electric Machines and Control，2017，21（1）：30.

[21]路義萍，張東學，孫博，等.某無刷勵磁機通風冷卻數(shù)值模擬研究[J].電機與控制學報，2016，20（6）：26.

LU Yiping，ZHANG Dongxue，SUN Bo，et al.Numerical simulation of ventilation cooling system for brushless exciter[J].Electric Machines and Control，2016，20（6）：26.

[22]溫嘉斌，鄭軍，于喜偉.YKK中型高壓異步電動機換熱分析與冷卻器優(yōu)化[J].電機與控制學報，2015，19（9）：34.

WEN Jiabin，ZHENG Jun，YU Xiwei.Heat transfer analysis and cooler optimization of YKK medium size high voltage asynchronous motor[J].Electric Machines and Control，2015，19（9）：34.

[23]孟大偉，何金澤，夏云彥，等.箱式Y(jié)JKK緊湊型中型高壓電動機全域流體預測[J].電機與控制學報，2012，16（12）：66.

MENG Dawei，HE Jinze，XIA Yunyan.All-around fluid forecast research of the compact medium size high voltage motor[J].Electric Machines and Control，2012，16（12）：66.

（編輯：劉琳琳）