城軌輔助變流器熱設(shè)計(jì)方法研究

楊　寧

(中國(guó)鐵道科學(xué)研究院機(jī)車(chē)車(chē)輛研究所，北京　100081)

?

城軌輔助變流器熱設(shè)計(jì)方法研究

楊寧

(中國(guó)鐵道科學(xué)研究院機(jī)車(chē)車(chē)輛研究所，北京100081)

針對(duì)輔助變流器熱設(shè)計(jì)的需求，建立有限體積法的迭代方程，并結(jié)合城軌輔助變流器的熱設(shè)計(jì)，提出熱仿真的設(shè)計(jì)步驟與方法。依照此方法建立北京城軌13號(hào)線輔助變流器箱體與內(nèi)部主要電氣部件的數(shù)學(xué)模型，采用仿真與公式計(jì)算結(jié)合的IGBT熱損耗計(jì)算方法，使用六面體網(wǎng)格離散時(shí)間與空間域的模型，借助流體力學(xué)軟件迭代有限體積法方程進(jìn)行熱仿真的計(jì)算，計(jì)算結(jié)果滿(mǎn)足溫升要求。依據(jù)此熱設(shè)計(jì)結(jié)果試制成樣機(jī)，完成穩(wěn)態(tài)溫升試驗(yàn)，并將主要電氣部件試驗(yàn)結(jié)果與仿真結(jié)果進(jìn)行對(duì)比。結(jié)果滿(mǎn)足設(shè)計(jì)初衷，驗(yàn)證了有限體積法的熱仿真分析方法的可行性。

城軌輔助變流器；熱設(shè)計(jì)；建模；有限體積法

對(duì)于城軌列車(chē)輔助變流器的設(shè)計(jì)而言，其箱體密封性好，運(yùn)行環(huán)境惡劣，安裝在車(chē)下等特定的需求，對(duì)發(fā)熱設(shè)計(jì)帶來(lái)更高的要求。傳統(tǒng)的試驗(yàn)結(jié)合經(jīng)驗(yàn)[1-2]的方式難免會(huì)帶來(lái)設(shè)計(jì)周期的加長(zhǎng)和成本的浪費(fèi)。目前國(guó)外結(jié)合計(jì)算機(jī)輔助軟件，將有限元數(shù)值分析引入熱設(shè)計(jì)的仿真中[3]，雖然其計(jì)算精度高，但求解時(shí)間較長(zhǎng)。也有使用有限差分的方法替代導(dǎo)數(shù)求解，但僅限于對(duì)規(guī)則的模型求解。本文綜合2種方法的優(yōu)點(diǎn)使用有限體積法對(duì)城軌輔助變流器進(jìn)行熱仿真設(shè)計(jì)方法的研究。

1　有限體積法方程的建立

建立的熱仿真方程必須遵守流體流動(dòng)的物理守恒定律，包括：質(zhì)量守恒定律、動(dòng)量守恒定律和能量守恒定律，控制方程就由三大定律構(gòu)成的基本方程，稱(chēng)作納維-斯托克斯方程(Navier-Stokes，N-S)，簡(jiǎn)稱(chēng)N-S方程。

N-S方程被離散化成一系列非重疊、連續(xù)的有限體積，溫度方程的微分形式為[4]

(1)

根據(jù)一維方程，擴(kuò)展到三維模型。圖1是一個(gè)一維x方向上均勻的網(wǎng)格。3個(gè)相鄰網(wǎng)格中，處于網(wǎng)格中心節(jié)點(diǎn)P與相鄰節(jié)點(diǎn)lx和hx控制體積的界面為lxf和hxf。

圖1　一維水平相鄰的網(wǎng)格

單元體積為

(2)

水平一維x方向的面積為

(3)

則式(1)可寫(xiě)為瞬態(tài)項(xiàng)、對(duì)流項(xiàng)、傳導(dǎo)項(xiàng)與源項(xiàng)的方程

(4)

式(4)可簡(jiǎn)化為

(5)

整理后寫(xiě)成線性形式

(6)

式中：

式(6)的含義是流出熱量等于流入熱量加上熱源。由一維方程可以推廣到三維方程，對(duì)于一般變量φ有方程

(7)

由式(7)可以看出：當(dāng)前溫度值φp是由上一步φt的值與其相鄰6個(gè)單元的溫度值迭代而成的結(jié)果。因此有限體積法不僅增加了運(yùn)算的速度，同時(shí)對(duì)運(yùn)算的資源要求也不高。

2　控制方程仿真的實(shí)現(xiàn)

根據(jù)N-S方程建立溫度控制方程后，采用有限體積法中心格式網(wǎng)格的劃分對(duì)控制方程進(jìn)行離散化，可以通過(guò)編寫(xiě)程序或者利用計(jì)算機(jī)流體力學(xué)(簡(jiǎn)稱(chēng)CFD)軟件對(duì)輔助變流器的箱體的模型進(jìn)行熱仿真計(jì)算。

使用CFD軟件熱仿真流程需要建立準(zhǔn)確的模型，模型的建立包括幾何模型的建立、材料參數(shù)的設(shè)定、器件功耗的計(jì)算。模型建立后，可對(duì)模型進(jìn)行網(wǎng)格的劃分與穩(wěn)態(tài)的求解。

本研究結(jié)合北京城軌13號(hào)線輔助變流器的熱設(shè)計(jì)過(guò)程對(duì)熱設(shè)計(jì)方法進(jìn)行了具體分析應(yīng)用。

3　輔助變流器模型的建立

3.1幾何模型的建立

幾何模型的建立是熱仿真結(jié)果準(zhǔn)確與否的關(guān)鍵。相對(duì)于城軌的牽引變流器，輔助變流器的幾何模型較復(fù)雜，內(nèi)部的機(jī)械部件和電氣部件種類(lèi)繁多。在系統(tǒng)級(jí)別的仿真下，為了得到便于計(jì)算又不失真的模型，需要對(duì)幾何模型機(jī)械部件和電氣部件的機(jī)械結(jié)構(gòu)進(jìn)行簡(jiǎn)化處理。簡(jiǎn)化處理的內(nèi)容包括：去除不必要的零件(螺絲、螺母、扎線桿、標(biāo)簽等)，合并小零件，鈑金件的改造(將圓倒角改為垂直角)，薄板的處理(薄板改成面或增加厚度)，封堵多余的縫隙與過(guò)孔，電氣部件的模型修改為便于賦值的幾何結(jié)構(gòu)等。

輔助變流器簡(jiǎn)化后的幾何模型見(jiàn)圖2。

圖2　輔助變流器簡(jiǎn)化后幾何模型

3.2材料參數(shù)的設(shè)定

將輔助變流器的幾何模型簡(jiǎn)化完成導(dǎo)入到CFD熱仿真軟件中后，對(duì)材料參數(shù)進(jìn)行逐個(gè)設(shè)定，包括材料密度、熱傳導(dǎo)率、電阻率等參數(shù)。其中主要材料的類(lèi)型見(jiàn)表1。其中熱管的模型按表2的參數(shù)進(jìn)行建模。

表1　主要部件材料

表2　熱管模型主要參數(shù)

3.3器件功耗計(jì)算

器件功率損耗計(jì)算的正確與否，決定著熱仿真結(jié)果是否對(duì)箱體設(shè)計(jì)具有指導(dǎo)意義。城軌輔助變流器發(fā)熱的主要器件有IGBT、工頻變壓器、輸入電抗器和輸出電抗器。

(1)IGBT的損耗

IGBT的損耗分為導(dǎo)通損耗與開(kāi)關(guān)損耗，以及內(nèi)部反并聯(lián)二極管的損耗。IGBT損耗的計(jì)算方法有很多種，為了得到較準(zhǔn)確的計(jì)算結(jié)果，采用公式求解和仿真2種方法進(jìn)行損耗計(jì)算，評(píng)估出合理的熱計(jì)算損耗。城軌輔助變流器采用雙極性的調(diào)制方式，根據(jù)調(diào)制方式得到IGBT的導(dǎo)通損耗計(jì)算公式為

(8)

式中，M為調(diào)制度；φ為電壓電流相位差；VCEO為門(mén)檻電壓；ICP為電流峰值；rCE為IGBT通態(tài)等效電阻。VCEO和rCE可以通過(guò)IGBT手冊(cè)的參數(shù)計(jì)算得出。

IGBT的開(kāi)關(guān)損耗計(jì)算公式為

(9)

式中，fsw為開(kāi)關(guān)頻率；Vdc為輸入直流電壓；Eon為在額定電流ICN和額定電壓VCEN下，IGBT開(kāi)通一次損失的能量；Eoff為在額定電流ICN和額定電壓VCEN下，IGBT關(guān)斷一次損失的能量。

反并聯(lián)二極管(簡(jiǎn)稱(chēng)VD)的損耗分為VD的導(dǎo)通損耗和關(guān)斷的反向恢復(fù)損耗。在雙極性下的導(dǎo)通損耗為

(10)

式中，VF0為VD門(mén)檻電壓；rF為VD通態(tài)等效電阻。VF0和rF可以通過(guò)IGBT手冊(cè)的參數(shù)計(jì)算得出。

VD關(guān)斷反向恢復(fù)損耗為

(11)

式中，Eoff為在額定電流ICND和額定電壓VCEND下VD關(guān)斷一次損失的能量。

IGBT的損耗還可根據(jù)廠家提供的選型軟件進(jìn)行近似的計(jì)算，并與公式計(jì)算進(jìn)行對(duì)比，評(píng)估出合理的損耗值。

(2)變壓器及電抗器的損耗

變壓器屬于電磁元件，由鐵芯和繞組組成。為了便于賦值，分別計(jì)算變壓器的鐵耗和銅耗。工頻變壓器的鐵耗分為基本鐵耗和雜散鐵耗，它們主要是由磁滯和渦流造成，受負(fù)載電流變化較小，可以設(shè)為定值。工頻變壓器的銅耗分為基本銅耗和雜散銅耗，銅耗中占大部分比例的損耗與負(fù)載電流有關(guān)。它與額定電流下的損耗有如下關(guān)系

(12)

式中，I0為計(jì)算電流；PCu為在I0電流下的銅耗；PC為在額定電流Irms下的損耗。

輸出電抗器的計(jì)算可以類(lèi)比工頻電抗器進(jìn)行計(jì)算。輸入電抗器的損耗分為鐵耗與銅耗。鐵耗與直流電源的品質(zhì)有關(guān)，對(duì)于直流輸入電流波動(dòng)較小，直流電壓比較平穩(wěn)，其鐵耗可以忽略。輸入電抗器的銅耗可根據(jù)計(jì)算工況的直流等效電阻，直流電流計(jì)算出電抗器的銅耗。

(3)其他損耗

輔助變流器的功率損耗除去主要電氣件的損耗，還包括功率電纜、母排、電容、雜散電感、接觸器觸點(diǎn)等的損耗，因其所占比重較小且不用考核，仿真時(shí)可忽略不計(jì)。

城軌13號(hào)線單臺(tái)輔助變流器的額定工況：輸入電壓750V，輸入電流153.6A，輸出電壓3AC383.7V，輸出電流3AC195.3A，功率因數(shù)0.849。得出單個(gè)IGBT的損耗為354W，工頻變壓器損耗為1 900W，輸入電抗器損耗為500W，輸出電抗器損耗為1 080W。效率95%，計(jì)算結(jié)果滿(mǎn)足預(yù)期設(shè)計(jì)要求。

4　網(wǎng)格的劃分與求解

在求解計(jì)算前需要對(duì)模型進(jìn)行網(wǎng)格的劃分。網(wǎng)格的劃分是把原來(lái)在時(shí)間域與空間域上連續(xù)物理量的場(chǎng)，用一系列有限個(gè)離散點(diǎn)上變量值的集合來(lái)代替，通過(guò)一定的原則和方式建立起關(guān)于這些離散點(diǎn)上場(chǎng)變量之間關(guān)系的代數(shù)方程組，然后求解代數(shù)方程組獲得場(chǎng)變量的近似值。網(wǎng)格的劃分就是規(guī)定計(jì)算域的集合，較密的網(wǎng)格更加準(zhǔn)確的描述參數(shù)的變化，但會(huì)加長(zhǎng)求解的時(shí)間，稀疏的網(wǎng)格對(duì)參數(shù)的變化描述的較粗糙，但會(huì)縮短求解的時(shí)間。因此，網(wǎng)格劃分的原則不是越密越好，而是選擇密度合理的網(wǎng)格。依據(jù)有限體積法方程的推導(dǎo)過(guò)程，采用六面體進(jìn)行網(wǎng)格的合理劃分。

采用數(shù)值方法進(jìn)行穩(wěn)態(tài)的求解分析，需要多次的迭代計(jì)算才能獲得近似的溫度分布結(jié)果，因此迭代結(jié)果的收斂有一定的準(zhǔn)則。模型第i個(gè)控制體的控制方程可寫(xiě)成

(13)

計(jì)算的新值為φr，則

(14)

(15)

式中，nc為總共的控制體數(shù)量；R為殘差。

收斂的準(zhǔn)則是：當(dāng)殘差小于一定值時(shí)計(jì)算收斂并且結(jié)束。一般取值為10。

網(wǎng)格劃分后可采用有限體積法進(jìn)行穩(wěn)態(tài)求解。文章采用Flotherm軟件，依據(jù)北京城軌13號(hào)線輔助變流器的額定工況參數(shù)得出仿真結(jié)果，見(jiàn)圖3。

圖3　輔助變流器熱仿真結(jié)果

5　樣機(jī)的驗(yàn)證

根據(jù)仿真結(jié)果得到各個(gè)部件的穩(wěn)態(tài)溫升在合理的范圍內(nèi)，將此模型試制成樣機(jī)，并以相同的額定工況進(jìn)行溫升試驗(yàn)，采用布點(diǎn)的方式對(duì)發(fā)熱器件進(jìn)行溫度的測(cè)量。鑒于樣機(jī)的溫升測(cè)試受限于采樣點(diǎn)的數(shù)量和位置，無(wú)法得到機(jī)箱內(nèi)部發(fā)熱器件準(zhǔn)確的溫度云圖，下面對(duì)測(cè)量點(diǎn)的結(jié)果與仿真結(jié)果進(jìn)行如下對(duì)比。實(shí)測(cè)點(diǎn)溫度見(jiàn)表3。

表3　實(shí)測(cè)點(diǎn)溫度

5.1對(duì)比分析

圖4給出了散熱器內(nèi)側(cè)的仿真溫度云圖，從圖中可以直觀地看出，散熱器內(nèi)側(cè)最高溫度59.5 ℃，與表3散熱器內(nèi)測(cè)點(diǎn)1溫度較接近，最低溫度41.7 ℃，與散熱器內(nèi)測(cè)點(diǎn)2溫度較接近。仿真結(jié)果與機(jī)組實(shí)測(cè)結(jié)果存在誤差，但測(cè)試點(diǎn)溫度在仿真結(jié)果的溫度范圍內(nèi)。

圖5給出了變壓器鐵心對(duì)應(yīng)測(cè)點(diǎn)附近的溫度，從圖中可以直觀地看出，圖中右下側(cè)最高溫度87.3 ℃，與變壓器鐵芯測(cè)點(diǎn)1溫度較接近。圖中左上側(cè)最低溫度80.3 ℃，與變壓器鐵芯測(cè)點(diǎn)2附近最高溫度較接近，仿真結(jié)果與機(jī)組實(shí)測(cè)存在誤差，但測(cè)試點(diǎn)溫度在仿真結(jié)果的溫度范圍內(nèi)。

圖6給出了輸入電抗器線包仿真溫度云圖，從圖中可以直觀地看出，輸入電抗器對(duì)應(yīng)測(cè)點(diǎn)1附近溫度118 ℃左右。輸入電抗器對(duì)應(yīng)測(cè)點(diǎn)2附近溫度109 ℃，仿真結(jié)果與機(jī)組實(shí)測(cè)基本一致，但存在誤差。

圖5　變壓器鐵芯溫度

圖6　輸入電抗器線包溫度

變壓器線包測(cè)點(diǎn)1與2溫度仿真值穩(wěn)定在132～133 ℃。仿真與實(shí)測(cè)較為一致。

輸出電抗器鐵心對(duì)應(yīng)實(shí)測(cè)點(diǎn)1附近溫度121 ℃。輸出電抗器鐵心對(duì)應(yīng)實(shí)測(cè)點(diǎn)2附近溫度117 ℃，仿真與實(shí)測(cè)較一致。

5.2誤差分析

城軌輔助變流器電氣主要部件的實(shí)測(cè)結(jié)果基本滿(mǎn)足仿真設(shè)計(jì)的預(yù)期，從定量的對(duì)比上存在誤差，產(chǎn)生這些誤差的原因是多方面的，主要如下：

(1)由于求解方程迭代造成的誤差；

(2)由于源項(xiàng)(功率損耗)計(jì)算結(jié)果造成的誤差，以及在此結(jié)果對(duì)幾何模型的分配上造成的誤差；

(3)對(duì)于網(wǎng)格劃分的疏密程度帶來(lái)的誤差；

(4)由于實(shí)測(cè)點(diǎn)位置的偏差，給測(cè)量帶來(lái)的誤差。

以上幾種誤差都會(huì)帶來(lái)實(shí)測(cè)點(diǎn)溫度與仿真溫度的誤差，但仿真結(jié)果與實(shí)測(cè)基本一致，可以指導(dǎo)樣機(jī)的研制。

6　結(jié)語(yǔ)

本研究采用有限體積法進(jìn)行熱仿真求解，以六面體網(wǎng)格簡(jiǎn)化迭代過(guò)程，加快了求解過(guò)程，較有限元和有限差分法迭代簡(jiǎn)單，具有較快的求解速度，節(jié)省了計(jì)算資源。根據(jù)北京城軌13號(hào)線輔助變流器的熱設(shè)計(jì)應(yīng)用實(shí)例，建立模型，劃分網(wǎng)格與求解得到仿真結(jié)果，依據(jù)仿真試制成樣機(jī)進(jìn)行溫升試驗(yàn)，試驗(yàn)結(jié)果表明主要電氣部件的溫升滿(mǎn)足仿真設(shè)計(jì)要求，驗(yàn)證了此熱設(shè)計(jì)方法的可行性。

[1]陳恩.電子設(shè)備熱設(shè)計(jì)研究[C]∥第十三屆全國(guó)熱泵與系統(tǒng)節(jié)能技術(shù)大會(huì)論文集.敦煌：中國(guó)制冷空調(diào)工業(yè)協(xié)會(huì)，2009:53-58.

[2]段生記，張敏強(qiáng). 電子機(jī)箱抗惡劣環(huán)境熱設(shè)計(jì)方法[C]∥中國(guó)電子學(xué)會(huì)電子機(jī)械工程分會(huì)2009年機(jī)械電子學(xué)學(xué)術(shù)會(huì)議論文集.太原：中國(guó)電子學(xué)會(huì)電子機(jī)械工程分會(huì)，2009：349-355.

[3]吳根盛.電動(dòng)車(chē)用DC/DC電源模塊有限元熱分析研究[J].汽車(chē)實(shí)用技術(shù)，2013(3):10-16.

[4]吳子牛.計(jì)算流體力學(xué)基本原理[M].北京:科學(xué)出版社，2001.

[5]楊寧，宋術(shù)全，李紅.基高速動(dòng)車(chē)組輔助變流器箱體的熱仿真設(shè)計(jì)方法[J].中國(guó)鐵道科學(xué)，2013(3):87-92.

[6]胡建輝，李錦庚.變頻器中的IGBT模塊損耗計(jì)算及散熱系統(tǒng)設(shè)計(jì)[J].電工技術(shù)學(xué)報(bào)，2009(3):159-163.

[7]楊寧，李巖磊，馬穎濤.城軌列車(chē)牽引變流器箱體熱仿真設(shè)計(jì)方法[J].鐵道機(jī)車(chē)車(chē)輛，2015(3):91-95.

[8]魏建華，李亞利.高溫環(huán)境應(yīng)用內(nèi)燃機(jī)車(chē)輔助變流器的設(shè)計(jì)[J].大功率變流技術(shù)，2016(1):17-20.

[9]左鵬，李紅，程建華.地面水冷卻循環(huán)裝置的研制[J].中國(guó)鐵道科學(xué)，2011，32(4):139-142.

[10]葛新石，等譯.傳熱學(xué)[M].北京：科學(xué)出版社，2002.

[11]呂棟，申雯雯，施科研，等.模塊化多電平變流器種子模塊散熱設(shè)計(jì)[J].電源學(xué)報(bào)，2012(2):42-46.

[12]徐殿國(guó)，李向榮.極限溫度下的電力電子技術(shù)[J].電工技術(shù)學(xué)報(bào)，2006，21(3):15-23.

[13]張舟云，徐國(guó)卿，沈祥林.牽引逆變器散熱系統(tǒng)的分析與設(shè)計(jì)[J].同濟(jì)大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，32(6):775-778.

[14]祝效華，余志樣.ANSYS高級(jí)工程有限元分析范例精選[M].北京:電子工業(yè)出版社，2004.

[15]何文志，丘東元，肖文勛，等.高頻大功率開(kāi)關(guān)電源結(jié)構(gòu)的熱設(shè)計(jì)[J].電工技術(shù)學(xué)報(bào)，2013，28(2):185-191，218.

Research on Thermal Design Method for Auxiliary Converter of Urban Rail Transit

YANG Ning

(Locomotive&CarResearchInstitute，ChinaAcademyofRailwaySciencesBeijing100081，China)

Inresponsedoauxiliaryconverterthermaldesignrequirements，theiterativeequationoffinitevolumemethodisestablishedanddesignproceduresandmethodforthermalsimulationarepresentedincombinationwiththeurbanrailtransitauxiliaryconverterthermaldesign.Basedonthismethod，mathematicalmodelsoftheauxiliaryconvertercaseandtheinternalmainelectricalcomponentsforBeijingSubwayLine13areestablished.WithIGBTheatlosscalculationmethodincombinationofsimulationandformulacalculation，thehexahedralgriddiscretizingtimeandspacedomainmodel，andbyvirtueoftheCFDsoftwareiteratingfinitevolumemethod，thethermalsimulationcalculationisconducted.Theresultssatisfytherequirementsofthetemperaturerise.Accordingtothethermaldesignresults，prototypeismade，thesteady-statetemperaturerisetestisconducted，andactualtestresultsandsimulationresultsarecompared.Theresultssatisfythedesignexpectationandverifythefeasibilityofthethermalsimulationfinitevolumeanalysismethod．

Auxiliaryconverterofurbanrailtransit;Thermaldesign;Modeling;Finitevolumemethod

2016-02-23；

2016-05-06

國(guó)家科技部科研院所基金(2012YJ009)

楊寧(1983—)，男，助理研究員，2013年畢業(yè)于中國(guó)鐵道科學(xué)研究院載運(yùn)工具與運(yùn)用工程專(zhuān)業(yè)，工學(xué)博士，主要從事機(jī)車(chē)車(chē)輛牽引與輔助供電方面的設(shè)計(jì)與研究工作，E-mail：yangning@zemt.cn。

1004-2954(2016)08-0165-05

U264.3+7

A< class="emphasis_italic">DOI

:10.13238/j.issn.1004-2954.2016.08.035