高壓開關(guān)柜溫度場數(shù)值分析

陳荊洲, 杜志葉, 王 棟， 任君鵬， 阮江軍

(武漢大學(xué)電氣工程學(xué)院，湖北 武漢 430072)

1 引言

高壓開關(guān)柜作為電力系統(tǒng)中重要的終端執(zhí)行元件，其安全運行決定著電力系統(tǒng)供電的安全性和可靠性[1-3]。在實際運行中，發(fā)熱問題是影響開關(guān)柜安全運行的關(guān)鍵因素，近年來由于過熱問題引起開關(guān)柜故障的案例越來越多，已經(jīng)引起了生產(chǎn)運營單位和研究機(jī)構(gòu)的廣泛關(guān)注[4-8]。高壓開關(guān)柜主要采用封閉式結(jié)構(gòu)，散熱性能較差，當(dāng)工作在高電壓、大電流條件下就會產(chǎn)生溫升，嚴(yán)重時將會嚴(yán)重影響開關(guān)柜的運行狀態(tài)[9]。目前，國內(nèi)外的工作主要集中在開關(guān)柜在線測溫裝置的研發(fā)和改進(jìn)上，應(yīng)用較廣泛的測溫方法包括傳統(tǒng)的溫度傳感器測溫法、紅外測溫法和光纖光柵測溫法[10-15]。

為進(jìn)一步探究開關(guān)柜溫度分布特點，合理選擇溫度監(jiān)測傳感器布置位置，提高過熱監(jiān)測的準(zhǔn)確性，需要獲得開關(guān)柜內(nèi)溫度的特征。目前在開關(guān)柜內(nèi)溫度分布計算方面，主要有平均溫升簡易算法、熱路方法和溫度場數(shù)值模擬算法三種。其中平均溫升簡易算法[16]是一種較粗略的估算方法，該方法形式簡單且易于使用，但是誤差略大于其余兩種方法，在實際應(yīng)用中有一定的指導(dǎo)價值；熱路方法[17，18]利用熱電類比法建立開關(guān)柜某部件的熱路模型然后求解溫度場，但是該方法求解的精度不是很高，需要進(jìn)一步修正完善；隨著計算機(jī)技術(shù)的飛速發(fā)展，通過多物理場耦合計算的數(shù)值模擬算法逐漸成為主流[19-21]。數(shù)值模擬算法雖然一定程度上提高了開關(guān)柜溫度場求解的準(zhǔn)確度，但目前的研究分析大多對開關(guān)柜模型作了極大的簡化或者單獨考慮某個部件的溫度場，沒有對開關(guān)柜的整體模型作全面的分析計算，同時也未考慮開關(guān)柜的風(fēng)機(jī)作用，會影響開關(guān)柜溫度場的求解精確度。

本文以10kV KYN-28A高壓開關(guān)柜為研究對象，通過Solidworks建立開關(guān)柜的實體模型，利用ANSYS Icepak對開關(guān)柜的溫度場、流體場進(jìn)行數(shù)值計算和分析，得到了開關(guān)柜整體的溫度分布特征，并且將仿真結(jié)果與溫升試驗數(shù)據(jù)進(jìn)行對比，驗證了計算方法的準(zhǔn)確性。

2 理論基礎(chǔ)及控制方程

開關(guān)柜內(nèi)的熱量主要來自于回路載流導(dǎo)體的焦耳產(chǎn)熱，導(dǎo)體產(chǎn)生的焦耳熱主要以熱傳導(dǎo)的方式在固體間傳遞，同時通過熱對流和輻射換熱傳遞到柜體內(nèi)空氣中，因此數(shù)值計算需要同時求解溫度場方程和流體場方程。

2.1 溫度場控制方程

開關(guān)柜內(nèi)的熱量傳遞包含三種基本方式：熱傳導(dǎo)、熱對流和熱輻射。三種傳熱方式的方程分別如下所示：

(1)

Q2=hcA(tw-tf)

(2)

Q3=δ0Aεxt(T14-T24)

(3)

式中，Q1、Q2、Q3分別為熱傳導(dǎo)熱量、熱對流換熱量和熱輻射換熱量；A為換熱面積；λ為導(dǎo)熱系數(shù)；hc為對流換熱系數(shù)；tw為固體表面溫度；tf為周圍空氣溫度；δ0為斯蒂芬—玻爾茲曼常數(shù)，其值為5.669×10-8W/(m2·k4)；εxt為系統(tǒng)發(fā)射率；T1為物體表面熱力學(xué)溫度；T2為環(huán)境溫度的熱力學(xué)溫度。

2.2 流體動力學(xué)控制方程

計算流體動力學(xué)一般包含三組控制方程，分別是質(zhì)量守恒方程、動量守恒方程和能量守恒方程。

質(zhì)量守恒方程：

(4)

動量守恒方程：

(5)

(6)

(7)

能量守恒方程：

(8)

式中，V為流體流速；u、v、w為流速在x，y，z方向上的分量；ρ為流體密度；μ為流體動力粘度；T為流體溫度；p為流體壓力；Cp為定熱容；Su、Sv、Sw為動量守恒方程廣義源項；ST為粘性耗散項。

3 仿真計算與分析

3.1 計算模型

本文利用Solidworks軟件針對10kV KYN-28A開關(guān)柜建立三維模型，實際開關(guān)柜的模型極其復(fù)雜，直接利用實物結(jié)構(gòu)模型進(jìn)行仿真計算非常困難，需要根據(jù)研究問題的側(cè)重點對實際開關(guān)柜模型進(jìn)行簡化，刪除一些對熱分析影響不大的部件，同時對一些復(fù)雜的結(jié)構(gòu)建立簡化的等效模型，建立的模型如圖1所示。

圖1 開關(guān)柜三維模型Fig.1 Three dimensional model of switchgear

開關(guān)柜模型包含四個獨立的隔離室，分別為母線室、架空進(jìn)線室、斷路器室和儀表室。斷路器和電流互感器部分均采用簡化模型，同時刪除了外殼及隔板上的安裝孔。在Workbeach平臺下，Icepak的標(biāo)準(zhǔn)CAD接口為DesignMolder(DM)，需要把模型導(dǎo)入到DM進(jìn)行修復(fù)處理，將所有的部件轉(zhuǎn)換為Icepak認(rèn)可的幾何體后，由DM導(dǎo)入到Icepak。開關(guān)柜實際模型中有6個風(fēng)機(jī)，其中5個位于頂部外殼，1個位于斷路器底部，Icepak具有獨立的風(fēng)機(jī)模塊，模型導(dǎo)入后需要在Icepak里添加風(fēng)機(jī)模型。

3.2 熱源分析

開關(guān)柜發(fā)熱的熱源主要包括載流導(dǎo)體和電接觸的焦耳損耗。文獻(xiàn)[22]表明電接觸部分的產(chǎn)熱對于開關(guān)柜溫升有重要影響，本文電接觸部分的發(fā)熱量將單獨計算。整個載流回路發(fā)熱功率為：

P=I2R

(9)

式中，I為流過載流導(dǎo)體的電流；R為載流導(dǎo)體電阻。載流導(dǎo)體的電阻由定義式計算，電接觸部分的接觸電阻受接觸表面狀況、接觸壓力等諸多因素的影響，文獻(xiàn)[23]給出了接觸電阻計算的經(jīng)驗公式：

Rc=Kc/(Fk/9.8)m

(10)

式中，Rc為接觸電阻；Kc為接觸材料系數(shù)；Fk為接觸壓力；m為與接觸形式有關(guān)的系數(shù)。

開關(guān)柜內(nèi)電接觸部分包括母排搭接處、母排與靜觸頭連接處和動觸頭梅花觸指。其中母排搭接處及母排與靜觸頭連接處通過高強(qiáng)螺栓進(jìn)行連接，均為面接觸，m值取1，同時接觸壓力可通過緊固力矩計算公式獲得；梅花觸指為點接觸，m值取0.5，接觸壓力通過查找開關(guān)說明技術(shù)參數(shù)獲得。經(jīng)過計算，當(dāng)開關(guān)柜流過額定電流4000A時各部分的熱量計算結(jié)果見表1。

表1開關(guān)柜各部分發(fā)熱量
Tab.1 Heating power of each part of switchgear

發(fā)熱部位發(fā)熱功率/W母線3309架空進(jìn)線7911斷路器7223電流互感器117

表1中所示發(fā)熱量均為三相總發(fā)熱量，也包括電接觸的發(fā)熱量。在仿真計算中，母排搭接處的發(fā)熱量換算成體密度加載在母排上，母排與靜觸頭連接處及動觸頭梅花觸指處的電接觸發(fā)熱量加載在接觸面上。

3.3 仿真計算

Solidworks建好后的實體模型導(dǎo)入 ANSYS Icepak后進(jìn)行模型預(yù)處理，開始網(wǎng)格剖分，通過精細(xì)控制散熱路徑和母排表面的網(wǎng)格尺寸和數(shù)量，來準(zhǔn)確模擬熱流的傳熱特性和流動特性，剖分后網(wǎng)格數(shù)量約為250萬單元。

母線室進(jìn)線母排、架空進(jìn)線母排、斷路器靜觸頭、動觸頭、觸臂的材質(zhì)為銅，絕緣套管、靜觸頭盒材質(zhì)為環(huán)氧樹脂，柜體外殼為鍍鋅鋼板。銅和鍍鋅鋼板的物理參數(shù)參考文獻(xiàn)[24]得到，環(huán)氧樹脂的物理參數(shù)參考文獻(xiàn)[25-27]整理得到。各種材質(zhì)的物理參數(shù)見表2。

表2材料物理參數(shù)
Tab.2 Physical parameters of materials

材料物理參數(shù)熱導(dǎo)率/(W/(M·K))密度/(g/cm3)比熱率/(J/(G·K))電阻率/(10-6Ω·m)銅3928903900179環(huán)氧樹脂027609814-鍍鋅鋼板4678050137

對各個部件定義物理參數(shù)后，按表1計算出的發(fā)熱功率進(jìn)行熱源加載。環(huán)境溫度設(shè)為25℃，根據(jù)軟件自動計算出的雷諾數(shù)和瑞利數(shù)選擇湍流模型，外殼設(shè)置對流換熱系數(shù)，采用經(jīng)驗值10W/(m2·K)。輸入實際的風(fēng)機(jī)P-Q曲線，計算模型較復(fù)雜，網(wǎng)格數(shù)量很多，為了加快求解速度采用雙核并行計算。

3.4 計算結(jié)果分析

3.4.1 溫度場分析

額定工況下開關(guān)柜的溫度分布如圖2所示，可以看出，斷路器本體是開關(guān)柜中溫度較高的部位，最高溫度出現(xiàn)在中間上動觸頭處，可達(dá)61℃，是開關(guān)柜溫度監(jiān)控的關(guān)鍵部位；母排的溫度在45℃左右，A相母排靠近風(fēng)扇，散熱面積大，溫度較B相和C相母排較低；母排與靜觸頭接頭由于與斷路器直接接觸，同時由于接觸電阻的影響，出現(xiàn)了較高溫升，其余載流回路的溫度在45℃-50℃之間。具體溫度分布見表3。

圖2 開關(guān)柜三維溫度分布圖Fig.2 Three dimensional temperature distribution of switchgear

部位名稱平均溫度/℃部位名稱平均溫度/℃上動觸頭565上靜觸頭560下動觸頭559下靜觸頭556母線462架空進(jìn)線425電流互感器485出線連接處496

3.4.2 流體場分析

開關(guān)柜三維流場中三個不同截面的流速矢量圖如圖3所示。其中，圖3(a)為z軸中心的xy平面圖；圖3(b)和3(c)分別為x=0.55和x=1.45時的yz平面圖。從流速矢量圖中可以得出流體流動的方向和速度大小。

圖3 開關(guān)柜截面速度矢量圖Fig.3 Velocity vector diagram of cross section

在柜頂風(fēng)機(jī)和斷路器室底部風(fēng)機(jī)的共同作用下，空氣從左部外殼進(jìn)風(fēng)口和前后外殼進(jìn)風(fēng)口進(jìn)入柜體內(nèi)。 從左板進(jìn)風(fēng)口進(jìn)入的空氣初始平均速度約為5.7m/s；而從前后板進(jìn)風(fēng)口進(jìn)入的空氣初始平均速度約為4.6m/s。從前后板進(jìn)風(fēng)口進(jìn)入的空氣大部分流入母線室和架空進(jìn)線室，然后從母線室頂部風(fēng)機(jī)流出，母線室和架空進(jìn)線室風(fēng)路較通暢，熱量散失較快，因而溫升較低。從左板進(jìn)風(fēng)口進(jìn)入的空氣在斷路器底部風(fēng)機(jī)的作用下大部分進(jìn)入斷路器室， 然后一部分會從斷路器外側(cè)繞經(jīng)儀表室，最后從斷路器室頂部流出，然而斷路器室結(jié)構(gòu)較為封閉，風(fēng)路狹窄，而且一部分氣流在儀表室下側(cè)形成環(huán)流，導(dǎo)致熱流路徑變長，風(fēng)速降低；另一小部分則通過上靜觸頭盒流至母線室，最后從母線室頂部流出，然而由于隔板的阻流影響，風(fēng)速也很低?？梢园l(fā)現(xiàn)，斷路器室封閉復(fù)雜的結(jié)構(gòu)延長了熱流路徑，降低了風(fēng)速，極大地降低了該處的散熱能力，導(dǎo)致斷路器部分整體溫升較高。

4 試驗對比

為了驗證仿真計算結(jié)果的準(zhǔn)確性，針對本開關(guān)柜開展了溫升試驗，溫升試驗中，電流為額定電流，即4000A，環(huán)境溫度為24.8℃。由于斷路器及觸頭部分是開關(guān)柜溫度較高的部位，因此測點主要布置在斷路器處，溫度通過溫度傳感器測得。共選取了9個溫度測點，分別為后出線連接處、電流互感器處、進(jìn)線與靜觸頭連接處、下靜觸頭、下動觸頭、上動觸頭、上靜觸頭、母線與靜觸頭連接處、上分支母線。測點布置如圖4所示。取三相平均溫度作為測點溫度，數(shù)字1-9依次代表上述9個測點，仿真計算結(jié)果與試驗結(jié)果對比見表4，對比如圖5所示。

圖4 開關(guān)柜溫度測點布置圖Fig.4 Arrangement of temperature measuring points of switchgear

測點編號平均溫度值/℃實測值計算值相對誤差(%)1474496462462485493529554474532556455555559076579565-2475495602853856244945246222

圖5 溫升計算值與實測值對比圖Fig.5 Experiment result compared with simulation result

從表4可以得出，仿真計算與試驗實測的溫度分布趨勢一致，最高溫度均在上動觸頭處。各個測點溫度值與計算溫度值也能較好地吻合，誤差均在5% 以內(nèi)，工程上可以接受，驗證了本文計算結(jié)果的正確性。

5 結(jié)論

(1)建立了開關(guān)柜有限元分析三維模型，計算了開關(guān)柜各個部位的發(fā)熱功率，并對熱流路徑和流固交界面的網(wǎng)格進(jìn)行了精細(xì)控制，在此基礎(chǔ)上利用Icepak對開關(guān)柜的溫度場和流體場進(jìn)行了數(shù)值分析。計算結(jié)果表明上動觸頭是開關(guān)柜溫度最高的部位，最高溫升達(dá)到36℃，是開關(guān)柜溫升控制和溫度監(jiān)控的關(guān)鍵。

(2)將仿真計算結(jié)果與溫升試驗結(jié)果進(jìn)行對比，計算結(jié)果與試驗結(jié)果吻合的很好，誤差在5%以內(nèi)，工程上可以接受，驗證了仿真計算的準(zhǔn)確性。

(3)對流體場進(jìn)行了分析，發(fā)現(xiàn)由于斷路器室結(jié)構(gòu)較為封閉，導(dǎo)致熱流路徑延長，風(fēng)速降低，極大的降低了斷路器室的散熱能力。在接下來的研究中，應(yīng)根據(jù)開關(guān)柜實際情況，優(yōu)化斷路器室的風(fēng)路結(jié)構(gòu)，減小熱流路徑，來控制該處溫升。

[1] 丁健(Ding Jian). 金屬封閉高壓開關(guān)柜觸頭發(fā)熱機(jī)理分析及預(yù)防過熱故障措施探討(Heating mechanism analysis on metal-enclosed high voltage switchgear contact with measures to prevent overheating)[J]. 高壓電器(High Voltage Apparatus)，2012，48(8)：114-122.

[2] 朱根良(Zhu Genliang). 淺議中壓開關(guān)柜事故調(diào)查中的故障分析(Simple explanation of faults analysis in incident investigation of medlum voltage switchgear)[J]. 高壓電器(High Voltage Apparatus)，2002，38(5)：62-63.

[3] 王有元，李寅偉，陸國俊，等(Wang Youyuan，Li Yinwei，Lu Guojun，et al.). 開關(guān)柜局部放電暫態(tài)對地電壓傳播特性的仿真分析(Simulation of transient earth voltage aroused by partial discharge in switchgears)[J]. 高電壓技術(shù)(High Voltage Engineering)，2011，37(7)：1683-1688.

[4] 李晶，滿家健，周電波，等(Li Jing，Man Jiajian，Zhou Dianbo，et al.). 金屬封閉高壓開關(guān)柜手車觸頭溫度場仿真分析及實驗研究(Simulation analysis and experimental research on temperature distribution of metal-enclosed high voltage switchgear handcart contactor)[J]. 高壓電器(High Voltage Apparatus)，2015，51(8)：74-79.

[5] 鞏憲鋒，衣紅鋼，王長松，等(Gong Xianfeng，Yi Honggang，Wang Changsong，et al.). 高壓開關(guān)柜隔離觸頭溫度監(jiān)測研究(Research on temperature monitoring of isolators in HV switchgear)[J]. 中國電機(jī)工程學(xué)報(Proceedings of the CSEE)，2006，26(1)：155-158.

[6] 周巖(Zhou Yan). 高壓開關(guān)柜觸頭溫度場的數(shù)值仿真分析(Numerical simulation of high voltage switchgear panel contact temperature)[J]. 傳感器與微系統(tǒng)(Transducer and Microsystem Technologies)，2008，27(12)：102-105.

[7] 宋帆，申春紅，徐建源，等(Song Fan，Shen Chunhong，Xu Jianyuan，et al.). 800kV GIS隔離開關(guān)磁場—溫度場計算與分析(Calculation and analysis on magneto- thermal fields of 800 kV GIS disconnector)[J]. 高電壓技術(shù)(High Voltage Engineering)，2008，34(7)：1383-1388.

[8] 熊蘭，徐敏捷，楊子康，等(Xiong Lan，Xu Mingjie，Yang Zikang，et al.). 高壓開關(guān)柜電纜室溫度場分析及在線監(jiān)測系統(tǒng)構(gòu)建(Analysis of temperature field in cable chamber of high voltage switchgear and construction of on line monitoring system)[J]. 電力自動化設(shè)備(Electric Power Automation Equipment)，2014，34(6)：153-157.

[9] 王新彤(Wang Xintong). 基于光纖光柵的開關(guān)柜溫度在線監(jiān)測系統(tǒng)的開發(fā)與應(yīng)用(The development and application of switchgear temperature monitoring system based on fiber bragg grating)[D]. 北京：華北電力大學(xué)(Beijing：North China Electric Power University)，2010.

[10] 李成榕，馬國明(Li Chengrong，Ma Guoming). 光纖布喇格光柵傳感器應(yīng)用于電氣設(shè)備監(jiān)測的研究進(jìn)展(Research and development of fiber bragg grating monitoring in electrical engineering)[J]. 中國電機(jī)工程學(xué)報(Proceedings of the CSEE)，2013，33(12)：114-122.

[11] 張艷，田競，葉逢春，等(Zhang Yan，Tian Jing，Ye Fengchun，et al.). 基于紅外傳感器的高壓開關(guān)柜溫度實時監(jiān)測網(wǎng)絡(luò)的研制(Development of temperature real-time monitoring network for high voltage switchgear based on IR sensor)[J]. 高壓電器(High Voltage Apparatus)，2005，41(2)：91-94.

[12] 鄔鋼(Wu Gang). 光纖在開關(guān)柜觸頭溫度監(jiān)測中的應(yīng)用(Application of optical fiber in temperature monitoring of the contacts of switchgear)[J]. 高電壓技術(shù)(High Voltage Engineering)，2006，32(2)：122.

[13] 趙炳成(Zhao Bingcheng). 高壓開關(guān)柜在線監(jiān)測系統(tǒng)的開發(fā)與應(yīng)用(Development and application of on-line monitoring of high voltage switchgear)[D]. 杭州：浙江大學(xué)(Hangzhou：Zhejiang University)，2011.

[14] 田毅，黃新波(Tian Yi，Huang Xinbo). 高壓開關(guān)柜溫升在線監(jiān)測系統(tǒng)(On-line monitoring system of temperature inside HV switchgear cabinet)[J]. 高壓電器(High Voltage Apparatus)，2010，46(3)：64-67,71.

[15] 張純璽(Zhang Chunxi). 光纖光柵溫度監(jiān)測系統(tǒng)設(shè)計及其在高低壓開關(guān)柜中的應(yīng)用(Design offiber Bragg grating temperature monitoring system and its application in high-voltage & low-voltage switchgear)[D]. 濟(jì)南：山東大學(xué)(Jinan：Shandong University)，2010.

[16] 黃恒禧(Huang Hengxi). 應(yīng)用高壓開關(guān)柜溫升的簡易算法及降低相應(yīng)發(fā)熱的技術(shù)措施(A study of simple arithmetic for applied high voltage switchboard temperature rise and the technical easures to reduce relevant heat)[J]. 福建建筑(Fu jian Architecture & Construction)，2008，117(3)：91-93.

[17] 趙健康，樊友兵，王曉兵，等(Zhao Jiankang, Fan Youbing, Wang Xiaobing，et al.). 高壓電力電纜金屬護(hù)套下的熱阻特性分析(Thermal resistance properties of the part between metal sheath and conductor in high-voltage power cable)[J]. 高電壓技術(shù)(High Voltage Engineering)，2008，34(11)：2483-2487.

[18] 孫凱祺，劉浩，傅晨釗，等(Sun Kaiqi， Liu Hao，F(xiàn)u Chenzhao，et al.). 220kV 瓷柱式斷路器修正熱路模型研究(Modified thermal circuit model for 220kV porcelain stanchion-type circuit breaker)[J]. 華東電力(East China Electric Power)，2013，41(6)：1215-1219.

[19] 廖才波，阮江軍，劉超，等(Liao Caibo， Ruan Jiangjun，Liu Chao，et al.). 油浸式變壓器三維電磁-流體-溫度場耦合分析方法(3-D coupled electromagnetic-fluid-thermal analysis of oil-immersed triangular wound core transformer)[J]. 電力自動化設(shè)備(Electric Power Automation Equipment)，2015，35(9)：150-155.

[20] 溫志偉，傅德平，顧國彪(Wen Zhiwei，F(xiàn)u Deping，Gu Guobiao). 浸潤式混合強(qiáng)迫內(nèi)冷的蒸發(fā)冷卻汽輪發(fā)電機(jī)定子三維溫度場的仿真分析(Simulation of 3D temperature distribution of immersion evaporative cooling turbo-generator combined with forced inner cooling system)[J]. 電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2006，25(3)：13-17.

[21] 鄒斌(Zou Bin). 基于有限元法的配電柜溫度場研究(Research on temperature field of switchboard based on finite element method)[D]. 蘇州：蘇州大學(xué)(Shuzhou：Shuzhou University)，2009.

[22] 王秉政，江健武，趙靈，等(Wang Bingzheng，Jiang Jianwu，Zhao Ling，et al.). 高壓開關(guān)柜接觸發(fā)熱溫度場數(shù)值計算(Temperature field simulation of contact heating in high voltage switchgear)[J]. 高壓電器(High Voltage Apparatus)，2013，49(12)：42-48.

[23] 郭鳳儀(Guo Fengyi). 電接觸理論及其應(yīng)用技術(shù)(Electrical contact theory and applied technology)[M]. 北京：北京電力出版社(Beijing：China Power Press)，2008.

[24] 徐國政(Xu Guozheng). 高壓斷路器原理與應(yīng)用(Principle and application of high voltage circuit break-er)[M]. 北京：清華大學(xué)出版社(Beijing：Tsinghua University Press)，2010.

[25] 程鵬，楊永明，莊曉蕓，等(Cheng Peng，Yang Yongming，Zhuang Xiaoyun，et al.). 立體卷鐵心干式變壓器三維溫度場計算與分析(Calculation and analysis of three-dimensional temperature field of Dry-Type transformer with tridimensional wound core)[J]. 變壓器(Transformer)，2013，50(9)：12-15.

[26] 蔡定國，姚育成，徐業(yè)彬，等(Cai Dingguo，Yao Yucheng，Xu Yebin，et al.). 基于流固耦合分析的干式變壓器溫度場數(shù)值分析(Numerical simulation of dry type transformer temperature field based on the fluid solid coupling analysis)[J]. 華北電力大學(xué)學(xué)報(Journal of North China Electric Power University)，2012，39(2)：34-38.

[27] 范春橋(Fan Chunqiao). 10kV 干式變壓器溫度場仿真計算與壽命評估(The simulation of temperature field and life assessment of 10kV Dry-Type transformer)[D]. 北京：華北電力大學(xué)(Beijing：North China Electric Power University)，2012.