丁杰，唐玉兔，忻力，胡昌發(fā)

（南車株洲電力機(jī)車研究所有限公司南車電氣技術(shù)與材料工程研究院，湖南株洲 412001）

輸出頻率對(duì)IGBT元件結(jié)溫波動(dòng)的影響

丁杰，唐玉兔，忻力，胡昌發(fā)

（南車株洲電力機(jī)車研究所有限公司南車電氣技術(shù)與材料工程研究院，湖南株洲 412001）

為分析輸出頻率對(duì)結(jié)溫波動(dòng)的影響，開發(fā)了基于模型降階方法的程序，運(yùn)用該程序?qū)GBT元件的熱模型進(jìn)行了計(jì)算，將得到的溫度計(jì)算結(jié)果與ANSYS穩(wěn)態(tài)計(jì)算結(jié)果進(jìn)行比較，結(jié)果的一致性證明該方法是有效的。在此基礎(chǔ)上，利用該程序計(jì)算得到了不同輸出頻率下IGBT元件的結(jié)溫波動(dòng)曲線。通過(guò)瞬態(tài)溫度與基于平均損耗計(jì)算的穩(wěn)定溫度對(duì)比可知：1 Hz輸出頻率時(shí)半正弦波損耗產(chǎn)生的溫度波動(dòng)范圍為-7.2%～7.6%，方波損耗產(chǎn)生的溫度波動(dòng)范圍為-3.9%～6.5%，輸出頻率對(duì)溫度波動(dòng)的影響顯著；50 Hz輸出頻率時(shí)對(duì)溫度波動(dòng)的影響較小，可以忽略；半正弦波損耗產(chǎn)生的溫度波動(dòng)大于方波損耗。分析結(jié)果可為IGBT元件瞬態(tài)溫度計(jì)算提供參考，基于模型降階的程序可以進(jìn)一步運(yùn)用于其他電力電子器件瞬態(tài)溫度的計(jì)算。

絕緣柵雙極晶體管；輸出頻率；結(jié)溫波動(dòng)；模型降階

1 引言

在大量的工程應(yīng)用中，為了簡(jiǎn)化計(jì)算，通常都是基于一個(gè)調(diào)制周期內(nèi)的平均值來(lái)對(duì)IGBT元件的損耗和結(jié)溫進(jìn)行分析，因此，輸出頻率的改變不會(huì)影響平均損耗和結(jié)溫的計(jì)算值。實(shí)際上，IGBT元件的損耗是隨負(fù)載電流的大小而變化，輸出頻率較高時(shí)產(chǎn)生的結(jié)溫波動(dòng)較小，輸出頻率較低時(shí)的結(jié)溫波動(dòng)會(huì)變大。結(jié)溫的波動(dòng)會(huì)產(chǎn)生交變的熱應(yīng)力，從而對(duì)IGBT元件的熱循環(huán)次數(shù)、可靠性及使用壽命產(chǎn)生一定的影響。

由于IGBT元件結(jié)構(gòu)復(fù)雜、散熱條件多樣，采用有限元法或有限體積法需要解算很大規(guī)模的矩陣方程，導(dǎo)致穩(wěn)態(tài)問(wèn)題的計(jì)算尚需要大量的時(shí)間，瞬態(tài)問(wèn)題則需要更為漫長(zhǎng)的計(jì)算時(shí)間和龐大的計(jì)算機(jī)資源。因此，實(shí)際應(yīng)用與文獻(xiàn)中較少涉及IGBT元件瞬態(tài)溫度的計(jì)算。

為提高計(jì)算速度并保證計(jì)算精度，開發(fā)了基于模型降階方法［1-7］的程序，對(duì)不同輸出頻率下的IGBT元件結(jié)溫波動(dòng)進(jìn)行了快速計(jì)算，分析了芯片最高溫度的變化特點(diǎn)和波動(dòng)情況。分析結(jié)果與開發(fā)程序可為電力電子器件瞬態(tài)溫度計(jì)算提供參考。

2 模型降階方法

圖1示出了基于Matlab軟件、ANSYS軟件和Krylov子空間投影法縮聚自由度進(jìn)行模型降階并快速計(jì)算的程序流程。

圖1 模型降階程序流程Fig.1 Program flow chart of model order reduction

首先通過(guò)CFD模型得到對(duì)流換熱系數(shù)分布，為ANSYS熱傳導(dǎo)有限元模型提供準(zhǔn)確的邊界條件［7］。模型中有m個(gè)熱源的損耗需要改變時(shí)，分別建立損耗均為0、第1個(gè)熱源為單位損耗、第2個(gè)熱源為單位損耗、…、第m個(gè)熱源為單位損耗的（m+1）個(gè)ANSYS有限元模型，再?gòu)倪@些文件中分別提取熱傳導(dǎo)矩陣和熱容矩陣，并保存為文件（Cond_File0.dat，Cap_File0.dat，Cond_File1.dat，Cap_File1.dat，Cond_File2.dat，Cap_File2.dat，…，Cond_Filem.dat，Cap_Filem.dat）。

在Matlab軟件中利用讀取文件命令，可以從Cond_File0.dat文件整理得到新的熱傳導(dǎo)矩陣Conductance0（n×n）、載荷列陣 LoadVector0（n×1），從 Cap_File0.dat文件得到熱容矩陣 Capci-tance0（n×n），n為ANSYS熱傳導(dǎo)有限元模型的節(jié)點(diǎn)數(shù)目。對(duì)其他矩陣文件進(jìn)行相同操作，得到相應(yīng)的熱傳導(dǎo)矩陣、載荷列陣和熱容矩陣。熱傳導(dǎo)矩陣需要進(jìn)行LU分解。載荷列陣用于設(shè)置熱源的損耗數(shù)值，通過(guò)LoadVector=LoadVector0+A*（LoadVectori-LoadVector0）+B*（LoadVectorj-Load-Vector0）+C*（LoadVectork-LoadVector0）+…，可以考慮第i個(gè)熱源的損耗為A單位、第j個(gè)熱源的損耗為B單位、第k個(gè)熱源的損耗為C單位等多個(gè)熱源同時(shí)作用的情況。

采用Krylov子空間投影法縮聚自由度后，可將n×n的矩陣降階為n×r的矩陣v，其中降階階數(shù)r遠(yuǎn)小于原有階數(shù)，因此，計(jì)算效率可比ANSYS有限元計(jì)算提高幾百至幾千倍，使得所需的計(jì)算時(shí)間和資源大大減少。降階階數(shù)對(duì)計(jì)算速度和計(jì)算精度有很大影響。降階階數(shù)過(guò)小時(shí)，計(jì)算速度很快，但計(jì)算精度變差；降階階數(shù)過(guò)大時(shí)，計(jì)算精度會(huì)更高，但會(huì)嚴(yán)重影響計(jì)算效率和存儲(chǔ)資源。降階階數(shù)的確定主要依賴于時(shí)間范圍、時(shí)間步長(zhǎng)、損耗是否變化［2，5］。

根據(jù)大量的工程應(yīng)用經(jīng)驗(yàn)可知，如果后續(xù)步驟中需要對(duì)改變損耗的情況進(jìn)行瞬態(tài)分析，且時(shí)間范圍很?。ㄈ鐣r(shí)間周期為0.001～1 s），相應(yīng)的時(shí)間步長(zhǎng)也很小，降階階數(shù)宜選取20～30。如果后續(xù)步驟中損耗不改變，且時(shí)間范圍很大（如0.001～100 s以上），相應(yīng)的時(shí)間步長(zhǎng)可以較大，降階階數(shù)可取5～10。其余情況，建議取10～20。

時(shí)間步長(zhǎng)可以采用定步長(zhǎng)和對(duì)數(shù)步長(zhǎng)的方法，對(duì)數(shù)步長(zhǎng)方法可以在起始時(shí)間布置較多的時(shí)間點(diǎn)，有利于提高結(jié)果的準(zhǔn)確性，因此，推薦采用對(duì)數(shù)步長(zhǎng)方法。當(dāng)使用某一固定的溫度T0作為初始值時(shí)，初始溫度場(chǎng)通過(guò)命令語(yǔ)句y0=v’*（T0*ones（n，1））來(lái)實(shí)現(xiàn)。降階后的方程可以使用常微分方程求解器（ode15s）進(jìn)行求解，時(shí)間范圍（t1～t2）的計(jì)算取step次，可以得到（step+1）×r的矩陣Y，再通過(guò)命令語(yǔ)句y_reduced=（v*Y’）’將結(jié)果投影到原模型上，得到的y_reduced是（step+1）×n的矩陣。t2時(shí)刻的計(jì)算結(jié)果作為下一時(shí)間范圍（t2～t3）的初始溫度場(chǎng)時(shí)，可以通過(guò)命令語(yǔ)句y0=v’*y_reduced（step+1，：）’來(lái)設(shè)置。重復(fù)上述步驟，即可以實(shí)現(xiàn)熱源不斷改變條件下的循環(huán)計(jì)算。

經(jīng)過(guò)結(jié)果投影到原模型，可以得到整個(gè)模型的節(jié)點(diǎn)溫度。可以利用Matlab的patch命令顯示溫度場(chǎng)的分布云圖，也可以利用max函數(shù)查找出最高溫度及對(duì)應(yīng)的節(jié)點(diǎn)編號(hào)，并以曲線的形式進(jìn)行結(jié)果顯示。

3 仿真條件

IGBT元件工作的過(guò)程中必須采取有效的散熱措施，如強(qiáng)迫風(fēng)冷、水冷、熱管散熱、微通道冷卻等形式，不同散熱方式的能力各不相同，散熱器上的溫度分布亦有很大區(qū)別。由于散熱條件和均溫性的組合甚多，為簡(jiǎn)化分析，可假設(shè)IGBT元件底面為某當(dāng)量換熱系數(shù)的散熱條件［8-9］。假設(shè)某IGBT元件底面當(dāng)量換熱系數(shù)為500 W/m2·℃，IGBT芯片的平均損耗為650 W，二極管芯片的平均損耗為350 W。依據(jù)IGBT元件的幾何結(jié)構(gòu)建立如圖2所示的ANSYS熱傳導(dǎo)有限元模型，網(wǎng)格單元數(shù)目80萬(wàn)，節(jié)點(diǎn)數(shù)目為91萬(wàn)。

圖2 ANSYS熱傳導(dǎo)有限元模型Fig.2 Heat conduction finite element model by ANSYS software

IGBT元件中的IGBT芯片和二極管芯片依次交替產(chǎn)生熱量，需要采用前面所述方法設(shè)置2個(gè)熱源的損耗數(shù)值，然后提取熱傳導(dǎo)矩陣和熱容矩陣，進(jìn)行相應(yīng)的矩陣操作、模型降階、計(jì)算和結(jié)果顯示。

4 基于平均損耗的計(jì)算結(jié)果

圖3是對(duì)ANSYS熱傳導(dǎo)有限元模型進(jìn)行穩(wěn)態(tài)仿真得到的溫度場(chǎng)分布結(jié)果，可以看出二極管芯片的溫度最高，為87.243℃。圖4是利用模型降階程序進(jìn)行時(shí)間范圍0.01～10 000 s的計(jì)算后，得到的最高溫度曲線。由溫度變化曲線可知，溫度先是迅速上升，然后變化趨勢(shì)逐漸變緩，于100 s之后穩(wěn)定在87.243℃。穩(wěn)定后的溫度值與ANSYS穩(wěn)態(tài)計(jì)算的溫度值非常吻合，說(shuō)明了模型降階方法可以得到非常準(zhǔn)確的計(jì)算結(jié)果。從計(jì)算時(shí)間來(lái)看，模型降階程序瞬態(tài)計(jì)算的時(shí)間約為ANSYS穩(wěn)態(tài)計(jì)算的3倍，模型降階方法在計(jì)算效率方面具有非常大的優(yōu)勢(shì)。

圖3 IGBT元件溫度場(chǎng)分布Fig.3 Temperature distribution of IGBT module

圖4 最高溫度變化曲線Fig.4 The highest temperature rise curve

5 輸出頻率對(duì)結(jié)溫的影響

針對(duì)三相橋臂的一個(gè)IGBT元件，輸出電流的正半周通過(guò)IGBT芯片時(shí)，IGBT芯片產(chǎn)生的損耗近似為半正弦波，輸出電流的負(fù)半周通過(guò)二極管芯片時(shí)，二極管芯片產(chǎn)生的損耗亦近似為半正弦波，IGBT芯片和二極管芯片分別在另外半個(gè)周期內(nèi)的損耗則為0。為了簡(jiǎn)化計(jì)算，通常假設(shè)損耗的脈沖波形為矩形或正弦半波［10-12］。

圖5是假設(shè)IGBT芯片和二極管芯片損耗為方波，輸出頻率1 Hz（周期1 s）條件下的最高溫度曲線。圖5a示出的芯片最高溫度變化曲線的規(guī)律與圖4基本一致，由于損耗在矩形波和0之間波動(dòng)，溫度亦產(chǎn)生相應(yīng)的波動(dòng)。圖5b是198～200 s之間IGBT芯片的最高溫度變化曲線。在正半周期的損耗為定值，溫度先是迅速升高，然后變化趨勢(shì)減緩；損耗由定值突變?yōu)?的時(shí)刻，溫度迅速下降；在負(fù)半周期的損耗為0，溫度逐漸下降。圖5c是198～200 s之間二極管芯片的最高溫度變化曲線，由于二極管芯片和IGBT芯片交替產(chǎn)生損耗，在正半周期的損耗為0，溫度逐漸下降，損耗由0突變?yōu)槎ㄖ档臅r(shí)刻，溫度開始上升；在負(fù)半周期的損耗為定值，溫度先是迅速升高，然后變化趨勢(shì)減緩。圖5d表示了198～200 s之間2個(gè)周期的波形，相當(dāng)于圖5b和圖5c疊加后的結(jié)果，最高溫度在83.8～92.9℃范圍內(nèi)變化，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-3.9%～6.5%。

圖5 輸出頻率1 Hz的方波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.5 The highest temperature change curve for the losses of square-wave in 1 Hz output frequency

圖6和圖7分別是輸出頻率5 Hz（周期0.2 s）和50 Hz（周期0.02 s）方波損耗條件下的最高溫度曲線，其曲線變化趨勢(shì)與輸出頻率1 Hz的基本一致，溫度范圍則有較大變化。輸出頻率5 Hz時(shí)，最高溫度在84.5～90.2℃范圍內(nèi)變化，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-3.1%～3.4%；輸出頻率50 Hz時(shí)，最高溫度變化范圍是86.9～88.3℃，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-0.4%～1.2%?？梢哉f(shuō)明溫度波動(dòng)的范圍隨著輸出頻率的增加而逐漸減小。

圖6 輸出頻率5 Hz的方波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.6 The highest temperature change curve for the losses of square-wave in 5 Hz output frequency

圖7 輸出頻率50 Hz的方波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.7 The highest temperature change curve for the losses of square-wave in 50 Hz output frequency

圖8是假設(shè)IGBT芯片和二極管芯片損耗為半正弦波，輸出頻率1 Hz（周期1 s）條件下的最高溫度曲線。圖8a示出的芯片最高溫度變化曲線的規(guī)律與圖5a基本一致，但由于損耗波動(dòng)幅度大，溫度產(chǎn)生相應(yīng)的波動(dòng)亦大于方波的情況。圖8b是198～200 s之間IGBT芯片的最高溫度變化曲線。在正半周期的損耗為半正弦波，溫度先是迅速升高，然后變化趨勢(shì)減緩，溫度在損耗為半正弦波峰值的時(shí)刻仍然繼續(xù)上升，經(jīng)歷一小段時(shí)間后開始下降；在負(fù)半周期的損耗為0，溫度逐漸下降。圖8c為198～200 s之間二極管芯片的最高溫度變化曲線。在正半周期的損耗為0，溫度逐漸下降；在負(fù)半周期的損耗為半正弦波，溫度先上升后下降，溫度的峰值滯后于損耗半正弦波峰值的時(shí)刻。產(chǎn)生芯片結(jié)溫的峰值相對(duì)于半正弦波損耗峰值有一定滯后的原因在于IGBT元件各層材料熱容的影響。圖8d相當(dāng)于圖8b和圖8c疊加后的結(jié)果，在198～200 s之間有2個(gè)周期的波形，最高溫度在81～93.9℃范圍內(nèi)變化，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-7.2%～7.6%。

圖8 輸出頻率1 Hz的半正弦波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.8 The highest temperature change curve for the losses of half-sine wave in 1 Hz output frequency

圖9和圖10分別是輸出頻率5 Hz（周期0.2 s）和50 Hz（周期0.02 s）半正弦波損耗條件下的最高溫度曲線。輸出頻率5 Hz時(shí)，最高溫度在84～90.6℃范圍內(nèi)變化，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-3.7%～3.9%；輸出頻率50 Hz時(shí)，最高溫度則在86.4～88.4℃范圍內(nèi)變化，與基于平均損耗計(jì)算的穩(wěn)定溫度相比，溫度波動(dòng)范圍為-1.0%～1.3%。同樣可以說(shuō)明隨著輸出頻率的提高，芯片最高溫度的波動(dòng)逐漸減小。

圖9 輸出頻率5 Hz的半正弦波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.9 The highest temperature change curve for the losses of half-sine wave in 5 Hz output frequency

圖10 輸出頻率50 Hz的半正弦波損耗對(duì)應(yīng)的最高溫度變化曲線Fig.10 The highest temperature change curve for the losses of half-sine wave in 50 Hz output frequency

由于半正弦波損耗的峰值是方波損耗值的π/2倍，通過(guò)對(duì)比輸出頻率分別為1 Hz，5 Hz和50 Hz時(shí)，方波損耗與半正波損耗導(dǎo)致的芯片最高溫度波動(dòng)范圍可知：半正波損耗產(chǎn)生的溫度波動(dòng)范圍大于方波損耗產(chǎn)生的溫度波動(dòng)；半正弦波損耗產(chǎn)生的最高溫度要低于方波損耗，當(dāng)輸出頻率增加至50 Hz時(shí)，兩者最高溫度波動(dòng)的上限值相等，但半正弦波損耗產(chǎn)生的最高溫度波動(dòng)的下限值仍小于方波損耗。

5 結(jié)論

從方波和半正弦波損耗條件下的芯片最高溫度變化趨勢(shì)和波動(dòng)范圍對(duì)比可知，盡管平均損耗相同，損耗波形不同時(shí)，IGBT元件的瞬態(tài)熱響應(yīng)有差異，這種差異隨著輸出頻率的提高而減小，隨著輸出頻率的降低而增大，對(duì)IGBT元件結(jié)溫波動(dòng)的影響不可忽略。半正弦波損耗條件產(chǎn)生的溫度波動(dòng)要大于方波損耗條件，在簡(jiǎn)化計(jì)算時(shí)，需要予以注意。此外，通過(guò)長(zhǎng)計(jì)算時(shí)間、多計(jì)算周期、損耗不斷變化的復(fù)雜瞬態(tài)問(wèn)題的快速有效計(jì)算，說(shuō)明基于模型降階的程序可以高效準(zhǔn)確地解決瞬態(tài)問(wèn)題，值得進(jìn)一步推廣應(yīng)用。

［1］ 關(guān)樂(lè)，褚金奎，齊東周，等.基于二階Krylov子空間投影法建立MEMS宏模型［J］.功能材料與器件學(xué)報(bào)，2008，14（1）：210-214.

［2］ Wilhelmus H A Schilders，Henk A van der Vorst，Joost Rom-mes.Model Order Reduction：Theory，Research Aspects and Applications［M］.Berlin：Springer，2008.

［3］ 關(guān)樂(lè)，褚金奎，齊東周，等.面向MEMS系統(tǒng)級(jí)仿真的宏模型研究［J］.系統(tǒng)仿真學(xué)報(bào)，2009，21（15）：4561-4567.

［4］ 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.

［5］ Yang Y J，Yu C C，Shen K Y.MEMS Heat Transfer Arnoldi-based Macromodels and the Study of Minimum Required Or-ders［J］.Tamkang Journal of Science and Engineering，2005，8（3）：185-190.

［6］ 齊東周，褚金奎，關(guān)樂(lè).MEMS降階宏模型仿真平臺(tái)研究［J］.微計(jì)算機(jī)信息，2009，25（3）：185-187.

［7］ 蔣耀林.模型降階方法［M］.北京：科學(xué)出版社，2010.

［8］ 丁杰，唐玉兔，忻力，等.網(wǎng)格尺寸對(duì)IGBT水冷散熱器仿真結(jié)果的影響［J］.大功率變流技術(shù)，2012，34（6）：26-30.

［9］ 丁杰，唐玉兔，忻力，等.IGBT模塊封裝的熱性能分析［J］.機(jī)車電傳動(dòng)，2013，52（2）：9-12.

［10］侯少敏.高速動(dòng)車組變流系統(tǒng)損耗分析與散熱研究［D］.北京：北京交通大學(xué)，2011.

［11］唐勇，陳明，汪波.電力電子器件短時(shí)脈沖工作的結(jié)溫特性研究［J］.電力電子技術(shù)，2010，44（3）：88-90.

［12］馮永杰，李斌，黃云.半導(dǎo)體器件功率老化的結(jié)溫控制方法研究［J］.實(shí)驗(yàn)技術(shù)與管理，2010，27（5）：46-49.

修改稿日期：2013-11-01

Effects of Output Frequency on Junction Temperature Ripple for IGBT Module

DING Jie，TANG Yu-tu，XIN Li，HU Chang-fa
（CSR Research of Electrical Technology&Material Engineering，CSR Zhuzhou Institute Co.，Ltd.，Zhuzhou412001，Hunan，China）

In order to analyze the effects of output frequency on the fluctuation of junction temperature，the program based on the model order reduction method was developed，the temperature results was obtained and compared with steady state calculation results by ANSYS software，the consistency of results shows that the method is effective.The temperature change curves of IGBT module in different output frequencies were obtained by the program.According to the comparison between the stable temperature based on the average losses and the transient temperature，the results show that the range of temperature fluctuations produced by half sine wave losses in 1 Hz output frequency is-7.2%～7.6%，and the range of temperature fluctuations produced by square wave losses is-3.9%～6.5%，the effects on the temperature of IGBT module in 1 Hz output frequency is significant and less impact for the temperature of IGBT module while the output frequency is 50 Hz.The temperature fluctuation caused by half sine wave losses is greater than square wave losses.The analysis results can provide reference for IGBT transient temperature calculation，the program based on model order reduction can be further used in other power electronics devices transient temperature calculation.

insulated gate bipolar transistor（IGBT）；output frequency；junction temperature ripple；model order reduction

TM46

A

湖南省自然科學(xué)省市聯(lián)合基金重點(diǎn)項(xiàng)目資助（12JJ8020）

丁杰（1979-），男，碩士，工程師，Email：dj8083@126.com

2013-08-06