電動汽車變流器用IGBT水冷散熱器熱仿真分析

丁杰，張平

?

電動汽車變流器用IGBT水冷散熱器熱仿真分析

丁杰1, 2，張平1

(1. 湘潭大學(xué)土木工程與力學(xué)學(xué)院，湖南湘潭，411105；2. 南車株洲電力機車研究所有限公司南車電氣技術(shù)與材料工程研究院，湖南株洲，412001)

以某電動汽車變流器用IGBT水冷散熱器為研究對象，考慮到IGBT水冷散熱器內(nèi)部各槽道平均雷諾數(shù)()處于過渡區(qū)，利用FLUENT軟件分別采用層流、標準?湍流模型和6種低數(shù)湍流模型計算IGBT水冷散熱器的穩(wěn)態(tài)結(jié)果。運用自主開發(fā)的模型降階計算程序?qū)GBT水冷散熱器進行瞬態(tài)問題快速計算，在驗證模型降階方法準確性的基礎(chǔ)上，對基于中國典型城市道路工況的瞬態(tài)熱仿真進行快速計算，得到IGBT元件各芯片在不同時刻的溫度變化曲線。研究結(jié)果表明：該水冷散熱器的散熱性能滿足中國典型城市道路工況的需求。

電動汽車；變流器；IGBT；水冷散熱器；模型降階；瞬態(tài)計算

為緩解能源和環(huán)境對人類生活與社會發(fā)展的壓力，世界各國相繼開發(fā)電動汽車[1]。變流器是電動汽車中的核心部件之一[2?3]，而IGBT水冷散熱器是保證變流器乃至電動汽車安全可靠運行的關(guān)鍵[4?5]。由于電動汽車的運行工況復(fù)雜多變，變流器內(nèi)部的IGBT元件產(chǎn)生的損耗不斷變化，初估的某一額定工況或極限工況具有較大隨意性，導(dǎo)致熱設(shè)計工作有較多不確定性。加之汽車行業(yè)具有嚴格的成本控制與輕量化要求，變流器及其內(nèi)部的IGBT水冷散熱器的結(jié)構(gòu)往往設(shè)計得非常緊湊，這給IGBT水冷散熱器進行精細化熱設(shè)計帶來了很多困難。由于IGBT水冷散熱器結(jié)構(gòu)較復(fù)雜，很難通過理論或經(jīng)驗公式準確計算出隨時間變化的瞬態(tài)溫度變化結(jié)果。若采用有限體積法進行溫度場的瞬態(tài)計算，則需選取較大網(wǎng)格尺寸來限制網(wǎng)格數(shù)量，否則會因龐大的計算機資源需求而導(dǎo)致瞬態(tài)計算相當困難；當選取較大網(wǎng)格尺寸生成較少數(shù)量的網(wǎng)格時，受流體邊界層和熱邊界層的影響，計算結(jié)果將有較大誤差，難以提高仿真結(jié)果的準確性[6]。為此，本文作者以某電動汽車變流器的IGBT水冷散熱器為研究對象，利用FLUENT軟件進行IGBT水冷散熱器的穩(wěn)態(tài)計算，并通過實驗結(jié)果驗證穩(wěn)態(tài)仿真結(jié)果的準確性及確定合適的流動狀態(tài)模型；在此基礎(chǔ)上，運用自主開發(fā)的模型降階計算程序?qū)χ袊湫统鞘械缆饭r的瞬態(tài)熱問題進行仿真計算。

1 變流器結(jié)構(gòu)

某電動汽車變流器的供電來自整車高壓配電和低壓蓄電池，司機指令經(jīng)整車控制器和CAN總線發(fā)送給變流器，從而對電機進行驅(qū)動。電機的速度和溫度信息則由傳感器采集后反饋給變流器。圖1所示為變流器結(jié)構(gòu)示意圖。變流器主要由機殼、外蓋、IGBT元件、水冷散熱器、支撐電容、控制板組件、驅(qū)動板組件、復(fù)合母排、銅排和接線座等構(gòu)成。3個型號為2MBI1400VXB-120P-50的IGBT元件安裝在水冷散熱器上，通過電路連接和控制策略構(gòu)成逆變功能。IGBT元件在工作過程中產(chǎn)生的損耗占據(jù)變流器中所有電氣部件損耗的絕大部分，且IGBT元件的熱量基本上通過水冷散熱器進行散熱，因此，IGBT元件和水冷散熱器是變流器熱設(shè)計工作的重點。

圖1 變流器結(jié)構(gòu)示意圖

2 穩(wěn)態(tài)計算

2.1 仿真模型

在利用FLUENT軟件進行穩(wěn)態(tài)計算之前，需要建立IGBT元件和水冷散熱器的網(wǎng)格模型?？紤]到IGBT元件各層材料的厚度均很小，且水冷散熱器內(nèi)部槽道寬度僅為3 mm，正方體網(wǎng)格邊長選取0.5 mm，以便能夠較準確地模擬水冷散熱器內(nèi)部的流動行為及邊界層的影響，從而獲得較準確的溫度場分布。利用HyperMesh軟件可以劃分出以六面體為主、極少量棱柱體的高質(zhì)量網(wǎng)格模型，劃分后的網(wǎng)格數(shù)量為1 185萬。

通過電氣參數(shù)可以簡單估算單個IGBT元件在額定工況下的損耗為1.127 8 kW，在極限工況下?lián)p耗為2.710 0 kW。冷卻介質(zhì)采用配比為50%(體積分數(shù)，下同)水和50%乙二醇的混合物，入口溫度為65 ℃，入口流量為20 L/min。冷卻介質(zhì)的熱物理參數(shù)如下：密度為1.043 t/m3，導(dǎo)熱系數(shù)為0.378 W/(m×K)，比熱容為3.345 J/(g×K)，動力黏度為1.525 g/(m×s)。通過手工計算可知該水冷散熱器入口的數(shù)為18 142。槽道的截面長×寬為3 mm×15 mm。假定冷卻介質(zhì)均勻分配在并行的10個槽道中，可計算出槽道內(nèi)的數(shù)為 2 533。根據(jù)流體力學(xué)知識可知，對于內(nèi)部流動而言，當數(shù)大于10 000時，可認為是湍流區(qū)；當數(shù)大于2 300而小于10 000時，應(yīng)屬于過渡區(qū)；當數(shù)小于2 300時，則為層流。丁杰等[7]分析當直槽道水冷散熱器內(nèi)部槽道結(jié)構(gòu)相同時，不同部位的流動阻力會有所區(qū)別，使得各槽道的流速不同，相應(yīng)的數(shù)也會不同，導(dǎo)致水冷散熱器內(nèi)部的流動狀態(tài)難以確定。在工程應(yīng)用中，通常采用的標準?模型基于完全湍流的假設(shè)，故又稱為高數(shù)?模型，需要對其修正以適應(yīng)低數(shù)的情況。自第1種低數(shù)?模型被提出以來，先后有10余種模型被提出，不同模型之間的區(qū)別主要在于：1) 考慮壁面對湍流脈動的衰減影響；2)和方程中附加項的構(gòu)造；3) 壁面上和方程邊界條件的處理。由于湍流問題的復(fù)雜性高，這些模型并未給出明確的適應(yīng)范圍[8]。為探討流動狀態(tài)模型對仿真結(jié)果的影響，本文選取常用的層流模型和標準?模型，還有FLUENT軟件提供的6種低數(shù)?模型，即Abid模型、Lan?Bremhorst模型、Launder Shorma模型、Yang?Shih模型、Abe?Kondoh?Nagano模型和Chang?Hsich?Chen模型。

2.2 流場仿真結(jié)果及分析

采用Lan?Bremhorst低數(shù)湍流模型進行計算得到的中間截面流速分布見圖2。從圖2可以看出冷卻介質(zhì)在槽道內(nèi)和轉(zhuǎn)折區(qū)域的流動情況，各槽道的冷卻介質(zhì)流速存在差異，最高流速可達2.32 m/s。冷卻介質(zhì)在槽道內(nèi)流動時產(chǎn)生的壓降為8.46 kPa。

2.3 溫度仿真結(jié)果及分析

圖3所示為采用Lan?Bremhorst低數(shù)湍流模型對額定損耗進行計算得到的帶IGBT元件的水冷散熱器溫度場分布。從圖3(a)可看出：靠近出口側(cè)的IGBT元件芯片溫度最高，為101.5 ℃。圖3(b)所示為水冷散熱器元件安裝面的溫度分布，IGBT元件安裝部位的溫度分布并不相同，IGBT元件端部位置出現(xiàn)熱點，其最高溫度為84.74 ℃。采用Lan?Bremhorst低數(shù)湍流模型對極限損耗進行計算可知：IGBT元件芯片最高溫度為151.97 ℃，該值已超出IGBT元件所能承受的結(jié)溫極限150.00 ℃，而水冷散熱器元件安裝面上的最高溫度可達109.69 ℃。

圖2 中間截面流速分布

(a) 水冷散熱器整體；(b) 水冷散熱器元件安裝面

圖4所示為冷卻介質(zhì)與水冷散熱器固體區(qū)域接觸面(流固耦合面)的對流換熱系數(shù)分布，用于說明冷卻介質(zhì)在不同部位的熱交換能力。從圖4可看出：內(nèi)部槽道的對流換熱系數(shù)分布并不均勻，最高值1.16×104W/(m2×K)出現(xiàn)在冷卻介質(zhì)流速值高的轉(zhuǎn)折區(qū)域。

圖4 對流換熱系數(shù)分布

3 實驗研究

3.1 實驗測試

通過不同的流動狀態(tài)模型計算相同工況(額定損耗或極限損耗)下的仿真結(jié)果。IGBT元件芯片最高溫度相差5.9 ℃，因流動狀態(tài)模型導(dǎo)致的芯片溫升差異可達7%，壓降相差0.9 kPa，這不利于水冷散熱器性能評價，對變流器的熱設(shè)計帶來較大影響。為此，對該水冷散熱器進行實驗研究。

以IGBT元件為熱源進行散熱器熱性能測試時存在以下問題：1) IGBT元件需要驅(qū)動、控制等復(fù)雜的系統(tǒng)組成才能正常工作；2) IGBT元件損耗特性對溫度敏感，溫度波動會帶來一定的試驗誤差；3) IGBT元件屬于精密器件，長時間極限工況運行很容易造成器件損傷；4) 芯片封裝于絕緣材料中，不能直接測量芯片結(jié)溫，在去除絕緣材料過程中很容易使IGBT元件遭受機械損傷；5) IGBT元件的價格高。這說明IGBT元件不宜直接作為水冷散熱器精細實驗的熱源，需要采用模擬熱源來代替。

在實驗準備過程中，在水冷散熱器安裝面上鉆18個直徑為0.8 mm的小孔，熱電偶放置在小孔中(見圖5)并與水冷散熱器緊密接觸；在水冷散熱器的模擬熱源安裝區(qū)域均勻涂敷一層厚約100mm的導(dǎo)熱硅脂，用于減小接觸熱阻；放置模擬熱源，將熱電偶線從模擬熱源底部的線槽中引出并逐一做好標記以示區(qū)別；最后用測力扳手將模擬熱源按照規(guī)定的方式與力矩擰緊。

圖5 熱電偶位置

模擬熱源上的電熱管引線聯(lián)至調(diào)壓電路，可手動調(diào)節(jié)電熱管的功率。水冷散熱器的入口和出口接上軟水管，與水泵、水箱、電磁流量計、壓力表等相連，電磁流量計可以測量實驗裝置的流量，壓力表可以測量管路系統(tǒng)和水冷散熱器的壓力1。熱電偶接入數(shù)據(jù)采集器的端口，并用數(shù)據(jù)線將數(shù)據(jù)采集器的通訊端口與計算機相聯(lián)，實現(xiàn)溫度的自動采集。

受實驗條件限制，很難實現(xiàn)水冷散熱器在設(shè)計條件下所規(guī)定的入口溫度與冷卻液配比，因此，實驗測試時大多采用室溫下的純水進行實驗。測溫點1，6，7，12，13和18在實驗工況1(調(diào)節(jié)單個模擬熱源的發(fā)熱量為1.128 kW，冷卻介質(zhì)采用入口流量20 L/min、溫度25 ℃的純水)得到的溫度曲線見圖6。由圖6可知：在0~100 s時，各測溫點的溫度迅速上升；100~300 s時，各測溫點的溫度上升速度明顯降低；300 s之后，各測溫點溫度上升幅度很小，這說明溫度已基本達到平衡；溫度平衡時，靠近冷卻水出口的測溫點13溫度最高，為34.1 ℃，靠近冷卻水入口側(cè)的測溫點6溫度最低，為30.7 ℃，其余測溫點的溫度處于測溫點6和測溫點13的溫度之間。

測溫點編號(見圖5)：1—1號；2—6號；3—7號；4—12號；5—13號；6—18號；7—入口；8—出口。

在進行水冷散熱器的性能實驗時，還進行了實驗工況2(單個模擬熱源發(fā)熱量為2.050 kW、冷卻介質(zhì)采用入口流量為20 L/min、溫度為25 ℃的純水)、實驗工況3(單個模擬熱源發(fā)熱量為1.143 kW、冷卻介質(zhì)采用入口流量為10 L/min、溫度為25 ℃的純水)、實驗工況4(單個模擬熱源發(fā)熱量為2.066 kW、冷卻介質(zhì)采用入口流量為10 L/min、溫度為25 ℃的純水)、實驗工況5(單個模擬熱源發(fā)熱量為1.153 kW、冷卻介質(zhì)采用入口流量為20 L/min、溫度為25 ℃的50%(體積分數(shù)，下同)乙二醇、50%水)和實驗工況6(實驗工況5進行至590 s，迅速調(diào)節(jié)單個模擬熱源發(fā)熱量至2.055 kW，其余條件不變)的溫升測量工作[9]。由于調(diào)壓電路的變化量需要手動調(diào)節(jié)，很難在每次實驗時保證完全一致，使得上述實驗工況的模擬熱源發(fā)熱量存在一些差異。在上述溫度測量過程中，還記錄了整體水路系統(tǒng)的壓力1。

為了測量水冷散熱器的壓降，待溫升測量工作完成后，拆除水冷散熱器，換成10 cm長的短接水管，此時，可以測出管路系統(tǒng)的壓力為2。因此，水冷散熱器的壓降Δ應(yīng)等于整體水路系統(tǒng)的壓力1減去管路系統(tǒng)的壓力2，即Δ=1?2。通過實驗測量可知冷卻介質(zhì)采用溫度為25 ℃的純水，入口流量采用 10 L/min和20 L/min時的壓降分別為8 kPa和21 kPa。

3.2 仿真結(jié)果與實驗結(jié)果對比

利用上述實驗工況1的條件作為仿真參數(shù)，參照文獻[9]中方法建立模擬熱源和水冷散熱器的網(wǎng)格模型，利用不同的流動狀態(tài)模型進行模擬，得到相應(yīng)仿真結(jié)果。采用Lan?Bremhorst低數(shù)湍流模型得到的帶模擬熱源的水冷散熱器溫度場分布見圖7。從圖7(a)可看出：模擬熱源上的電熱管溫度最高，為47.3 ℃。圖7(b)所示為水冷散熱器上模擬熱源安裝面的溫度場分布，3個模擬熱源安裝部位的溫度分布并不相同，局部熱點的最高溫度為36.44 ℃。

(a) 整體水冷散熱器；(b) 水冷散熱器模擬熱源安裝面

按照測溫點位置提取不同流動狀態(tài)模型的仿真結(jié)果，繪制成如圖8所示的溫度對比曲線。從圖8可知：層流模型和標準?模型的各測溫點溫度是各種流動狀態(tài)模型中的上、下限，低數(shù)湍流模型的溫度稍微高于標準?模型的溫度，Lan?Bremhorst模型的溫度與實驗溫度吻合度最高；層流以及標準?，Abid，Lan?Bremhorst，Launder Shorma，Yang?Shih，Abe?Kondoh?Nagano和Chang?Hsich?Chen模型得到的水冷散熱器壓降分別為7.14，8.11，7.54，7.53，7.74，7.32，7.53和7.56 kPa，這些壓降均小于實驗結(jié)果的主要原因在于：水冷散熱器的入口和出口均安裝了水嘴接頭，由于水嘴接頭會產(chǎn)生較大壓降，仿真時難以建立水嘴接頭的仿真模型而未考慮其影響。從水嘴接頭的產(chǎn)品數(shù)據(jù)手冊可知單個水嘴接頭在20 L/min純水流過時的壓降約為6.6 kPa，由此可知Lan?Bremhorst模型的壓降仿真結(jié)果與實驗結(jié)果接近。對基于實驗工況2~6的仿真結(jié)果分別與實驗結(jié)果進行對比，對比結(jié)果均可說明Lan?Bremhorst模型在該水冷散熱器的穩(wěn)態(tài)計算中具有好的適用性。

1—層流模型；2—Launder Shorma模型；3—標準k?ε模型；4—Lan?Bremhorst模型；5—Abid模型；6—Yang?Shih模型；

4 瞬態(tài)計算

4.1 基于模型降階方法的快速計算程序

由于利用FLUENT軟件針對網(wǎng)格數(shù)量多的模型進行流場與溫度場計算時，需要不斷地迭代計算而消耗大量的計算時間，如在DELL T7600臺式工作站上進行8處理器并行計算穩(wěn)態(tài)問題時，迭代3 000次基本可以達到10?6殘差收斂準則，消耗的時間約為15 h，未采用并行計算時的消耗時間達62 h。為實現(xiàn)快速計算瞬態(tài)問題，人們對模型降階方法進行了研究，并取得了好的應(yīng)用效果[10?15]?；贛atlab軟件、ANSYS軟件和Krylov子空間投影法進行模型降階并快速計算的自主開發(fā)程序流程見圖9[16?17]。首先通過精細網(wǎng)格的FLUENT模型得到對流換熱系數(shù)分布，為粗糙網(wǎng)格的ANSYS熱傳導(dǎo)有限元模型提供準確的邊界條件。當模型中有多個熱源的損耗需要改變時，依次設(shè)置每個熱源的發(fā)熱量為單位值，其余熱源的發(fā)熱量為0，再分別從這些有限元模型中提取熱傳導(dǎo)矩陣和熱容矩陣文件。在Matlab軟件中讀取矩陣文件并整理出新的熱傳導(dǎo)矩陣(×階)、載荷列陣(×1階)和熱容矩陣(×階)(其中，為ANSYS熱傳導(dǎo)有限元模型的節(jié)點數(shù)目)。采用Krylov子空間投影法將×階矩陣降階為×階矩陣(＜)，使用常微分方程求解器進行求解，然后將結(jié)果投影到原模型上，從而得到有限元模型每個節(jié)點在不同時刻的溫度結(jié)果。重復(fù)熱源設(shè)置、求解和投影等步驟，即可實現(xiàn)多熱源變工況的瞬態(tài)計算。降階階數(shù)主要取決于時間步長，通常取5~50，這遠小于原有階數(shù)，從而使計算效率比ANSYS有限元計算效率高幾百至幾千倍，所需的計算時間和資源大大減少。

圖9 快速計算程序流程

4.2 基于定值損耗的溫升計算

針對IGBT元件與水冷散熱器建立ANSYS熱傳導(dǎo)有限元模型時，正方體網(wǎng)格邊長為3 mm，得到的模型網(wǎng)格數(shù)為23萬，節(jié)點數(shù)為26萬。基于前面穩(wěn)態(tài)計算所述的結(jié)果，將圖4所示的對流換熱系數(shù)插值到ANSYS有限元模型中，進行額定損耗的穩(wěn)態(tài)計算可得到溫度場分布，IGBT元件上的芯片最高溫度為101.117 ℃，與圖3中的IGBT元件芯片最高溫度101.5 ℃相比非常接近，說明對流換熱系數(shù)由精細網(wǎng)格向粗糙網(wǎng)格插值處理的有效性。

以額定損耗為輸入，利用自主開發(fā)程序計算0.01~1 000 s時間段的瞬態(tài)問題，經(jīng)過數(shù)十秒計算即可得到有限元模型每個節(jié)點在不同時刻的溫度。為便于說明瞬態(tài)計算結(jié)果，提取IGBT元件上的最高溫度可繪制成如圖10所示的溫升曲線。由圖10可知：0.01~100 s內(nèi)溫度迅速上升，完全穩(wěn)定的溫度為101.117 ℃，所需時間為310 s。穩(wěn)定后的溫度與ANSYS穩(wěn)態(tài)計算溫度完全一致，說明自主開發(fā)程序計算的準確性。

圖10 以額定損耗為輸入的瞬態(tài)最高溫度計算結(jié)果

在以額定損耗為輸入、溫度穩(wěn)定的情況下，以極限損耗為輸入，利用自主開發(fā)程序計算0.01~1 000 s時間段的瞬態(tài)問題，提取IGBT元件上的最高溫度繪制成如圖11所示的溫升曲線。從圖11可見：其溫升變化規(guī)律與圖10所示的相同，完全穩(wěn)定的溫度為151.78 ℃；當溫度由101.117 ℃到達125 ℃需5.5 s，由101.117 ℃到達150 ℃需62.6 s。這些結(jié)果有利于確定IGBT元件由額定損耗工況穩(wěn)定運行情況下切換至極限損耗工況安全可靠運行的時間。

圖11 以極限損耗為輸入的瞬態(tài)最高溫度計算結(jié)果

對模擬熱源與水冷散熱器建立ANSYS熱傳導(dǎo)有限元模型時，正方體網(wǎng)格邊長為3 mm，得到的模型網(wǎng)格數(shù)為24萬，節(jié)點數(shù)為28萬。ANSYS有限元模型的穩(wěn)態(tài)計算結(jié)果與FLUENT穩(wěn)態(tài)計算結(jié)果相近。以實驗工況1的損耗為輸入，快速計算得到的溫升曲線與圖6的實驗溫升曲線規(guī)律基本一致。

4.3 基于中國典型城市道路工況的溫升計算

公交運營具有頻繁起步和停車、發(fā)動機長時間處于怠速工況、低速重載和持續(xù)運行等特點。以中國典型城市道路工況[18]進行主電路仿真，可以得到如圖12所示的不同時刻電氣參數(shù)(如電機相電壓、電機電流、功率因數(shù)和輸出頻率)，由此可進一步計算出不同時刻的IGBT芯片損耗和二極管芯片損耗[19](見圖13)。將IGBT元件的芯片損耗作為自主開發(fā)程序的輸入，經(jīng)過10余min計算可以得到3個IGBT元件的芯片溫度變化曲線(見圖14)。由于不同時刻的電氣參數(shù)發(fā)生變化，產(chǎn)生的損耗相應(yīng)發(fā)生變化，與之對應(yīng)的芯片溫度亦發(fā)生改變。在基于中國典型城市道路工況運行過程中，芯片溫度未超過115 ℃，說明該水冷散熱器的散熱性能可以滿足中國典型城市道路工況的需求。此外，所得的瞬態(tài)溫升結(jié)果可為后續(xù)IGBT模塊疲勞壽命預(yù)估提供依據(jù)。由于在計算過程中得到有限元模型每個節(jié)點在不同時刻的溫度，可為全面深入分析模型不同部位的溫升變化提供依據(jù)。

(a) 電機電流；(b) 輸出頻率；(c) 電機相電壓；(d) 功率因數(shù)

圖13 單個IGBT元件的損耗變化曲線

圖14 3個IGBT元件的芯片溫度變化曲線

5 結(jié)論

1) 通過層流模型、標準?模型和6種低數(shù)湍流模型的穩(wěn)態(tài)計算結(jié)果與實驗結(jié)果對比分析，結(jié)果表明Lan?Bremhorst低數(shù)湍流模型在該水冷散熱器中具有較強適用性。

2) 采用對流換熱系數(shù)插值方法可以保證網(wǎng)格數(shù)目較少的ANSYS熱傳導(dǎo)有限元模型能獲得與網(wǎng)格數(shù)目多的FLUENT模型基本一致的溫度。

3) 基于模型降階方法的快速計算自主開發(fā)程序不僅具有很高的計算效率，而且可準確得到與實驗測試的溫升具有相同變化規(guī)律的計算結(jié)果。

4) 基于中國典型城市道路工況的瞬態(tài)熱仿真計算可以獲得不同時刻的芯片溫度，水冷散熱器的散熱性能滿足中國典型城市道路工況的需求，計算溫度較基于初估的某一額定損耗和極限損耗得到的穩(wěn)態(tài)溫度更具現(xiàn)實意義。

[1] 楊衛(wèi)華. 電動汽車永磁同步電動機及其控制器研究[D]. 武漢: 湖北工業(yè)大學(xué)電氣與電子工程學(xué)院, 2009: 1?12. YANG Weihua. Research on permanent magnet synchronous motor and its controller in electric vehicle[D]. Wuhan: Hubei University of Technology. School of Electrical and Electronic Engineering, 2009: 1?12.

[2] 楊利輝. 電動汽車驅(qū)動變流器的優(yōu)化設(shè)計[D]. 重慶: 重慶大學(xué)電氣工程學(xué)院, 2004: 25?34. YANG Lihui. Optimized design of electric vehicle drive-motor controller[D]. Chongqing: Chongqing University. School of Electrical Engineering, 2004: 25?34.

[3] 邵桂欣, 張承寧. 電動汽車用續(xù)流增磁變流器溫升的改善方法[J]. 微電機, 2007, 40(11): 29?31. SHAO Guixin, ZHANG Chengning. Improvement elevated temperature of enhanced magnetism motor controller applied to electrics vehicles[J]. Micromotors Servo Technique, 2007, 40(11): 29?31.

[4] 張舟云, 徐國卿, 沈祥林. 用于電動汽車的電機和驅(qū)動器一體化冷卻系統(tǒng)[J]. 同濟大學(xué)學(xué)報(自然科學(xué)版), 2005, 33(10): 1367?1371. ZHANG Zhouyun, XU Guoqing, SHEN Xianglin. Uniform cooling system of motor and its controller used in electric vehicles[J]. Journal of Tongji University (Natural Science), 2005, 33(10): 1367?1371.

[5] 樂智. 純電動汽車電機驅(qū)動系統(tǒng)的冷卻系統(tǒng)設(shè)計與研究[D]. 天津: 河北工業(yè)大學(xué)機械工程學(xué)院, 2010: 23?32. YUE Zhi. Design and research on the cooling system of pure electric vehicle motor drive system[D]. Tianjin: Hebei University of Technology. School of Mechanical Engineering, 2010: 23?32.

[6] 丁杰, 唐玉兔, 忻力, 等. 網(wǎng)格尺寸對IGBT水冷散熱器仿真結(jié)果的影響[J]. 大功率變流技術(shù), 2012(6): 26?30. DING Jie, TANG Yutu, XIN Li, et al. Effects of mesh size on simulation results of water-cooling radiator for IGBT module[J]. High Power Converter Technology, 2012(6): 26?30.

[7] 丁杰, 李江紅, 陳燕平, 等. 流動狀態(tài)與熱源簡化方式對IGBT水冷板仿真結(jié)果的影響[J]. 機車電傳動, 2011(5): 21?25. DING Jie, LI Jianghong, CHEN Yanping, et al. Effects of flow regime and heat source simplifications on simulation results for water-cooling plate of IGBT[J]. Electric Drive for Locomotives, 2011(5): 21?25.

[8] 陶文銓. 數(shù)值傳熱學(xué)[M]. 2版. 西安: 西南交通大學(xué)出版社, 2001: 332?409. TAO Wenquan. Numerical heat transfer[M]. 2nd ed. Xi’an: Xi’an Jiaotong University Press, 2001: 332?409.

[9] 丁杰, 張平. 電機控制器用IGBT水冷散熱器溫升實驗與熱仿真[J]. 大功率變流技術(shù), 2015(3): 23?28. DING Jie, ZHANG Ping. Temperature-rise test and thermal simulation for IGBT water-cooled radiator of motor controller[J]. High Power Converter Technology, 2015(3): 23?28.

[10] WILHELMUS H A S, HENK A V, JOOST R. Model order reduction: theory, research aspects and applications[M]. Berlin:Springer, 2008: 1?72.

[11] 關(guān)樂, 褚金奎, 齊東周, 等. 面向MEMS系統(tǒng)級仿真的宏模型研究[J]. 系統(tǒng)仿真學(xué)報, 2009, 21(15): 4561?4567. GUAN Le, CHU Jinkui, QI Dongzhou, et al. Review of research on macromodel for MEMS system-level simulation[J]. Journal of System Simulation, 2009, 21(15): 4561?4567.

[12] RUDNYI E B, KORVINK J G. Model order reduction for large scale engineering models developed in ANSYS[J]. Lecture Notes in Computer Science, 2006, 3732(1): 349?356.

[13] YANG Y J, YU C C, SHEN K Y. MEMS heat transfer arnoldi-based macromodels and the study of minimum required orders[J]. Tamkang Journal of Science and Engineering, 2005, 8(3): 185?190.

[14] 齊東周, 褚金奎, 關(guān)樂. MEMS降階宏模型仿真平臺研究[J]. 微計算機信息, 2009, 25(3): 185?187. QI Dongzhou, CHU Jinkui, GUAN Le. On the tool of MEMS CAD system-level simulation[J]. Microcomputer Information, 2009, 25(3): 185?187.

[15] 蔣耀林. 模型降階方法[M]. 北京: 科學(xué)出版社, 2010: 77?122. JIANG Yaolin. Model order reduction method[M]. Beijing: Science Press, 2010: 77?122.

[16] 丁杰, 唐玉兔. 模型降階方法在瞬態(tài)熱仿真中的應(yīng)用[J]. 機車電傳動, 2014(5): 51?55. DING Jie, TANG Yutu. Application research of model order reduction method in transient thermal analysis[J]. Electric Drive for Locomotives, 2014(5): 51?55.

[17] 丁杰, 唐玉兔, 忻力, 等. 輸出頻率對IGBT元件結(jié)溫波動的影響[J]. 電氣傳動, 2014, 44(6): 76?80. DING Jie, TANG Yutu, XIN Li, et al. Effects of output frequency on junction temperature ripple for IGBT module[J]. Electric Drive, 2014, 44(6): 76?80.

[18] 汪偉, 劉凌. 電動公交客車動力總成技術(shù)的發(fā)展及展望[J]. 大功率變流技術(shù), 2012(2): 1?4. WANG Wei, LIU Ling. Development and prospect of powertrain technology for electric bus[J]. High Power Converter Technology, 2012(2): 1?4.

[19] 丁杰, 忻力, 何海興, 等. 二電平IGBT變頻器載荷工況下的損耗和溫升計算程序開發(fā)[J]. 大功率變流技術(shù), 2014(1): 1?4. DING Jie, XIN Li, HE Haixing, et al. Program development of losses and temperature rise for two-level IGBT converter under loading conditions[J]. High Power Converter Technology, 2014(1): 1?4.

(編輯 陳燦華)

Thermal analysis of IGBT water-cooling radiator for electric vehicle converter

DING Jie1, 2, ZHANG Ping1

(1. College of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105, China;2. CSR Research of Electrical Technology & Material Engineering,CSR Zhuzhou Institute Co. Ltd., Zhuzhou 412001, China)

An IGBT water cooling radiator used in an electric vehicle converter was studied. Considering that the average Reynolds number of each channel was in the transition region, steady state temperature of the radiator was calculated by using the laminar model, the standard-turbulence model and six kinds of different low Reynolds number turbulence models with FLUENT software, respectively. Transient problem of the radiator was calculated by using self-developed model order reduction calculation program. Through verifying the accuracy of model order reduction method, transient thermal simulation of typical Chinese urban road conditions was carried out to obtain the temperature variation curves of each chip in IGBT device at different time. The results show that the heat dissipating performance of the water cooling radiator can meet the demand of typical Chinese urban road conditions.

electric vehicle; converter; IGBT; water-cooling radiator; model order reduction; transient calculation

10.11817/j.issn.1672?7207.2017.02.033

U464.138

A

1672?7207(2017)02?0525?08

2016?02?13；

2016?04?20

湖南省自然科學(xué)省市聯(lián)合基金重點資助項目(12JJ8020)(Project(12JJ8020) supported by the Hunan Province and Municipal Joint Key Natural Science Foundation)

丁杰，博士研究生，高級工程師，從事一般力學(xué)與力學(xué)基礎(chǔ)、變流器結(jié)構(gòu)仿真與熱仿真研究；E-mail：dj8083@126.com