计及流-热耦合热网络模型的IGBT结温计算*
Calculation of IGBT Junction Temperature with Thermal Network Model Considering Flow-thermal Coupling
通讯作者: 葛兴来,男,1979年生,博士,教授。主要研究方向为电力牵引传动系统故障建模、诊断及稳定性分析。E-mail:xlgee@163.com
收稿日期: 2022-02-28 修回日期: 2022-05-15
| 基金资助: |
|
Received: 2022-02-28 Revised: 2022-05-15
作者简介 About authors
许智亮,男,1999年生,硕士研究生。主要研究方向为功率变流器中IGBT器件可靠性。E-mail:
绝缘栅双极型晶体管(Insulated gate bipolar transistor,IGBT)的结温计算对于评估其剩余寿命、保障电力电子变换器可靠运行具有重要意义。传统的热网络法忽略了模块与散热器之间真实的传热过程,难以准确有效计算结温变化,为此,提出一种计及流-热耦合的热网络模型的IGBT结温计算方法。该方法通过考虑流-热耦合的有限元仿真提取得到IGBT模块与散热器的热网络模型,根据易于获取的散热器水温作为温度参考点计算IGBT结温。为验证模型有效性,搭建了功率循环加速老化试验平台提取结温变化曲线。结果表明,所提的热网络模型能够更加真实准确地计算IGBT结温,可以为实际应用中IGBT模块的可靠性评估、剩余寿命预测提供更加真实准确的结温变化数据。
关键词:
The junction temperature calculation of insulated gate bipolar transistor(IGBT) is of great significance for evaluating its remaining life and ensuring the reliable operation of power electronic converters. The traditional thermal network method ignores the real heat transfer process between the module and the radiator, and it is difficult to calculate the junction temperature change accurately and effectively. To address this, a calculation method of IGBT junction temperature with thermal network model considering flow-thermal coupling is proposed. In this method, the thermal network model of IGBT module and the radiator is extracted by the finite element simulation considering the flow-thermal coupling, and the junction temperature of the IGBT is calculated based on the easily obtained radiator water temperature as the temperature reference point. In order to verify the validity of the model, a power cycle accelerated aging test platform is built to extract the junction temperature change curve. The results show that the proposed thermal network model can calculate the IGBT junction temperature more realistically and accurately. And it can provide more real and accurate junction temperature change data for the reliability evaluation and remaining life prediction of IGBT modules in practical applications.
Keywords:
本文引用格式
许智亮, 葛兴来, 李金, 柯倩霞, 朱丹.
XU Zhiliang, GE Xinglai, LI Jin, KE Qianxia, ZHU Dan.
1 引言
目前国内外关于利用热网络法进行IGBT结温计算已有大量研究[10⇓⇓-13]。其中,文献[10]利用IGBT模块的数据手册拟合得到器件热网络模型,用以计算牵引变流器中的IGBT结温,但是数据手册中的热网络数据与实际器件具有一定误差。进一步,文献[11-12]对IGBT模块进行有限元分析(Finite element analysis,FEA),同时考虑模块内各芯片自热和耦合热阻抗,提取更精准的热网络参数。文献[13]更加深入地考虑实际工况中的IGBT全桥模块的热耦合作用,利用FEA评估了芯片之间、芯片与桥臂以及相与相之间的热耦合作用,但是缺乏对散热条件的考虑。为此,文献[14]单独考虑散热器与环境的对流热阻,利用小波包变换和马尔科夫链对对流热阻实现随机模拟,并将对流热阻与实际器件损耗匹配。然而利用统计学方法模拟得到对流热阻依赖于算法准确性,未曾从实际传热过程中的热耦合进行考虑。尽管热网络法结温监测已考虑热耦合效应和散热器条件,但是与真实模拟模块内热耦合以及模块与散热器热耦合仍存在一定差距。
综上,本文在已有研究的基础上,利用考虑IGBT模块与散热器之间传热影响的有限元仿真,提出一种计及流-热耦合的热网络模型的IGBT结温计算方法。第2节对1 200 V/50 A模块及相应水冷散热器进行有限元仿真,考虑流体与温度场耦合对IGBT结温的影响,分别提取IGBT模块的热网络模型和计及流-热耦合的IGBT模块和散热器热网络模型。第3节搭建功率循环试验平台,提取壳温及散热器水温分别代入热网络模型计算结温,并与试验结温对比验证。结果表明,所提方法计算所得IGBT结温更贴近真实结温,且更易于实际应用。
2 有限元仿真和热网络模型搭建
2.1 IGBT模块及散热器
图1
图2
表1 IGBT模块三维尺寸参数
| 结构层 | 材料 | 长/mm | 宽/mm |
|---|---|---|---|
| IGBT芯片 | Si | 7.24 | 6.90 |
| 二极管芯片 | Si | 5.55 | 5.55 |
| IGBT焊层 | Sn-Ag | 7.24 | 6.90 |
| 二极管焊层 | Sn-Ag | 5.55 | 5.55 |
| 上铜层 | Cu | 28.5 | 25.8 |
| 绝缘衬板 | Al2O3 | 30.65 | 28 |
| 下铜层 | Cu | 28.5 | 25.8 |
| 基板 | Cu | 91.4 | 31.4 |
表2 IGBT模块及散热器的材料特性
| 材料 | 热导率/ [W/(m·K)] | 密度/ (kg/m3) | 常压比热容/ [J/(kg·K)] |
|---|---|---|---|
| Si | 118 | 2 330 | 700 |
| Sn-Ag | 35 | 7 300 | 226 |
| Cu | 398 | 8 960 | 225 |
| Al2O3 | 15 | 3 780 | 830 |
| Al | 237 | 2 700 | 897 |
图3
2.2 有限元仿真
根据IGBT模块与散热器的结构及三维尺寸,利用有限元仿真软件ANSYS Workbench搭建IGBT及散热器的三维有限元模型。其中IGBT模块及散热器的网格划分选用Mechanical选项,折中考虑仿真精度与仿真时间成本,本文在IGBT芯片等对仿真结果影响较大的部分采用精细化密集网格划分,在模块其他部分采用较粗略的扫掠式网格划分,流道及流体采用CFD物理选项进行网格划分,网格划分的结果如图4所示。
图4
本文分别对IGBT模块进行双向流-热耦合仿真和瞬态热仿真,其中双向流-热耦合仿真为散热器施加1.5 m/s、15 ℃的冷却水,根据数据手册中IGBT的功率损耗及芯片三维尺寸,直接分别为IGBT芯片加载30 A和60 A负载电流下的生成热载荷4.5×109 W/m3、9×109 W/m3,同时借助Workbench中的System Coupling模块,实现流体仿真与瞬态热仿真之间传热界面数据的双向传输更新。瞬态热仿真中,需要为IGBT模块基板底部设置与实际散热器相当的散热条件。
对流换热过程中的定性温度tm为
式中,tw为水冷散热基板的平均壁温,参考环境温度设置为25 ℃;tf为冷却水的温度,设置为15 ℃。
根据定性温度tm查阅水物性参数表可知,水的运动黏度υm为1.008 7×10-6 m2/s,水的导热系数为λm=0.599 W/(m·K),普朗特数Prm=7,为了确定管道内的冷却水流体状态,需要计算散热器中冷却水纵掠管道内壁的雷诺数。
式中,v为冷却水流速,根据已选型的水冷机计算得出为1.5 m/s;L为纵掠平壁情况下平壁板长,本文中该参数为散热基板中水道的长度,即100 mm。
在根据工程经验临界雷诺数取Rem=5×105,由于Rem<5×105且0.6<Prm<60,故整个纵掠平壁可均视为层流强迫,纵掠平壁时整个平壁表面的努谢尔数Num计算公式按POHLHAUSEN[16]给出的平均对流换热系数关联式计算
将式(2)和式(3)代入对流换热系数求解公式可得,纵掠平壁时的对流换热系数为
故水冷基板中的冷却水与IGBT基板的对流换热系数h为2 937 W/(m2·K)。因此,在瞬态热仿真中为IGBT基板底部设置3 000 W/m2的对流换热 系数。
IGBT芯片在导通过程中表现为焦耳热效应,因此可将IGBT芯片等效为内热源,根据热力学传热微分方程,可得在有限元仿真软件中模块温度分布求解公式为
式中,λ为材料热导率;T为模块内部温度;q为热流量;ρ为材料密度;cp为材料热容。
图5
图6
根据不同电流等级下的温度分布云图可知,流-热耦合仿真和瞬态热仿真最高结温存在较大差异,根据IGBT模块的加热过程提取模块的升温曲线如图7所示。
图7
对比流-热耦合仿真与瞬态热仿真中的升温曲线可知,由于流-热耦合仿真中存在IGBT模块与底部散热器的实时热交换,而瞬态热仿真中将散热条件利用集中参数法等效为一个对流换热系数,缺乏对真实环境的模拟,两种仿真中的IGBT模块升温曲线变化趋势存在较大差异,这将影响提取得到的热网络参数,进而影响结温计算的准确性。
2.3 热网络提取
为了根据壳温或环境温度等易于获取的温度参考点计算结温,需要对IGBT模块或IGBT模块及散热器整体建立热网络模型。
常用的热网络模型有Foster热网络和Cauer热网络模型,相较于需要根据IGBT模块实际结构及各层材料物理参数搭建的Cauer热网络模型,本文选用更易于在实际中应用的四阶Foster热网络模型,Foster热网络模型不具备实际物理意义,根据IGBT模块的温度变化曲线与热网络定义式即可拟合得出,图8为IGBT模块与散热器的热网络模型示意图。
图8
通常对IGBT模块与散热器分别建立热网络,但是实际应用中壳温变化数据不易提取,因此本文拟利用散热器中的冷却水温作为温度参考点,在考虑流-热耦合的基础上,建立包含IGBT模块和散热器的四阶热网络模型。四阶热网络模型中热阻抗计算公式如式(6)所示
式中,Zth为热阻抗;Ri为第i阶热阻;Ci为第i阶热容;t为时间变量。
结合第2.2节中提取的IGBT结温升温曲线,同时提取流-热耦合仿真中的冷却水温变化曲线,瞬态热仿真中的壳温变化数据,利用Matlab软件拟合得到流-热耦合仿真和瞬态热仿真中的四阶热网络模型,具体参数如表3所示。
表3 流-热耦合仿真和瞬态热仿真中的热网络模型参数
| 热网络参数 | 流-热耦合仿真 | 瞬态热仿真 |
|---|---|---|
| R1/(K/W) | 0.101 7 | 0.136 |
| R2/(K/W) | 0.169 9 | 0.086 12 |
| R3/(K/W) | 0.146 3 | 0.204 8 |
| R4/(K/W) | 0.080 17 | 0.050 59 |
| C1/(J/K) | 0.377 | 0.225 9 |
| C2/(J/K) | 0.024 61 | 0.208 6 |
| C3/(J/K) | 0.562 9 | 0.423 5 |
| C4/(J/K) | 0.466 7 | 0.440 6 |
3 试验验证及结果分析
3.1 加速老化试验平台搭建
本文根据IEC标准(IEC60749-34)搭建直流功率循环加速老化试验平台,试验原理如图9所示。
图9
待测器件(Device under test,DUT)栅极的直流电源产生足够的驱动电压与驱动电流使其开通,加热电流利用陪跑IGBT模块S1控制,当陪跑IGBT模块开通时,较大的加热电流通过待测DUT的IGBT芯片,而后陪跑IGBT模块S1转换为截止状态,进而使待测IGBT模块不再流经负载电流。同时对待测IGBT芯片通入较小的测量微电流(100 mA),通过测量IGBT芯片的热敏电参数获取IGBT芯片结温。当加热电流断开时,对IGBT芯片进行散热处理,当IGBT芯片下降至指定最低结温时,使陪跑IGBT模块S1再次导通,开始新一轮的结温波动,以此模拟IGBT芯片在实际工作环境中受到的热冲击。
图10
图11
当两个二极管D1和D2具有相同的特性时,选择合适的放大器增益时,可利用放大器输出电压Vop表示IGBT饱和导通压降Vce。
加速老化试验平台装置如图12所示,包括可编程直流源Chorma,待测模块FF50R12RT4和陪跑模块FF75R12RT4,工业水冷机和红外测温仪。
图12
3.2 结温波动曲线提取
基于Vce(T)法对IGBT模块进行结温监测时,需要建立IGBT芯片结温Tj与IGBT饱和导通压降Vce之间的线性关系式,故在对IGBT模块进行结温在线监测前,需要对IGBT结温Tj与IGBT饱和通态压降Vce的线性关系进行校正,并拟合得到两者之间的结温系数。对DUT进行离线结温系数校正的试验如图13所示。
图13
本文利用恒温加热平台将IGBT模块加热至指定温度,并利用使IGBT芯片处结温保持稳定,此时IGBT驱动模块输出栅极控制信号,并利用IGBT结温监测电路板为IGBT模块通入100 mA微电流,确保IGBT模块处于开通状态且不会产生较大热损耗对结温测量造成影响,并利用红外测温仪得到IGBT结温,最后利用Matlab软件将IGBT结温与饱和通态压降进行数据拟合,以一阶线性函数为作为拟合函数,拟合得到结温校正曲线,拟合结果表明IGBT结温与饱和导通压降Vce具有较好的线性函数关系。
3.3 热网络法结温计算验证
利用可编程直流源分别输出30 A电流和60 A电流,并分别采集升温阶段5 s内待测IGBT模块的结壳温变化曲线,利用热电偶提取试验平台中IGBT模块的壳温变化曲线,同时与FEA中的壳温变化曲线对比。
由图14可知,试验中的壳温变化曲线与FEA中的相接近,因此将试验中IGBT模块的壳温变化数据输入第2.3节中根据瞬态热仿真所建立的热网络中,同时将散热器水温输入根据流-热耦合建立的热网络中,计算得到IGBT模块的结温变化曲线。
图14
由图15可知,计及流-热耦合的热网络模型计算的结温更接近实物试验结温,基于瞬态热仿真建立的热网络模型的结温计算高于实际结温,这是由于集总参数法等效散热条件时,对真实参数存在一定程度的简化,导致搭建的热网络模型与真实情况存在一定差异。此外,基于瞬态热仿真的热网络模型需要以壳温作为温度参考点进行计算,而计及流-热耦合的热网络模型可以选用散热器内冷却水温作为温度参考点,在电气化交通领域无疑具有更大的应用前景。
图15
4 结论
为了方便准确地获取IGBT模块结温,本文在传统热网络法的基础上,考虑散热器流体散热效果,提出一种计及流-热耦合的热网络结温计算方法。基于有限元仿真软件实现流-热耦合仿真并提取热网络模型,搭建了功率循环加速老化试验平台,对所提方法与传统方法进行对比验证。对仿真及试验结果分析得到以下结论。
(1) 所提计及流-热耦合的热网络模型在不同负载电流下对IGBT模块的结温计算均具有较高的 精度。
(2) 所提计及流-热耦合的热网络法基于散热器冷却水温作为参考温度计算IGBT结温,在实际工程应用中更易于选取温度参考点,具有较好的工程应用前景。
参考文献
轨道交通用高压IGBT技术特点及其发展趋势
[J].
Technical features and development trend of high-voltage IGBT for rail transit traction application
[J].
Condition monitoring of IGBT module for reliability improvement of power converters
[C]//
Study and handling methods of power IGBT module failures in power electronic converter systems
[J].DOI:10.1109/TPEL.2014.2373390 URL [本文引用: 1]
功率IGBT模块的寿命预测
[J].
Lifetime prediction of power IGBT module
[J].
地铁车辆牵引逆变器IGBT模块的结温与疲劳寿命计算
[J].
Junction temperature and fatigue life calculation of IGBT module for metro vehicle traction inverter
[J].
Temperature measurements of semiconductor devices-A review
[C]//
An adaptive thermal equivalent circuit model for estimating the junction temperature of IGBTs
[J].DOI:10.1109/JESTPE.2018.2796624 URL [本文引用: 1]
基于热敏感电参数法的大容量IGBT模块动态结温在线检测研究
[J].
Dynamical junction temperature online extraction with thermal sensitive electrical parameters for high power IGBT modules
[J].
Thermal coupling analysis in a multichip paralleled IGBT module for a DFIG wind turbine power converter
[J].DOI:10.1109/TEC.2016.2614526 URL [本文引用: 1]
基于SVPWM调制方式的列车牵引逆变器功率器件寿命预测
[J].
Life prediction of power devices of train traction inverter based on SVPWM modulation
[J].
Thermal networks generation and application in IGBT module packaging
[C]//
Improved thermal couple impedance model and thermal analysis of multi-chip paralleled IGBT module
[C]//
电动汽车用IGBT全桥模块(6 in 1)的热耦合作用机制
[J].
Thermal coupling mechanism of IGBT full-bridge module(6 in 1)for electric vehicles
[J].
计及环境对流随机性的功率器件结-环境热网络模型
[J].
Junction-to-ambient thermal network model of power devices considering randomness of thermal convective environment
[J].
Der Wärmeaustausch Zwischen festen Körpern und Flüssigkeiten mit kleiner Reibung und Wärmeleitung. Z angew Math. Mech 1:115-121
[J].
The heat exchange between solid bodies and liquids with low friction and heat conduction
[J].
一种基于反向串联稳压二极管钳位的IGBT导通压降在线监测电路
[J/OL].DOI:10.13334/j.0258-8013.pcsee.211259
An online monitoring circuit for IGBT conduction voltage drop based on reverse series zener diodes clamping
[J/OL].DOI:10.13334/j.0258-8013.pcsee.211259
/
| 〈 |
|
〉 |
