CN113435090B - 一种基于工况的igbt模块电-热-流体多场耦合仿真方法 - Google Patents
一种基于工况的igbt模块电-热-流体多场耦合仿真方法 Download PDFInfo
- Publication number
- CN113435090B CN113435090B CN202110728570.2A CN202110728570A CN113435090B CN 113435090 B CN113435090 B CN 113435090B CN 202110728570 A CN202110728570 A CN 202110728570A CN 113435090 B CN113435090 B CN 113435090B
- Authority
- CN
- China
- Prior art keywords
- loss
- chip
- simulation
- igbt module
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于工况的IGBT模块电‑热‑流体多场耦合仿真方法,首先建立IGBT模块内部芯片级的损耗模型;其次建立描述IGBT模块内部热传导过程的热模型;然后建立描述散热系统中流体的热对流过程的共轭传热模型;最后基于COMSOL Multiphysic with MATLAB编写电‑热‑流体之间耦合变量的交互程序,本发明能够准确仿真出模块在不同工况下的瞬态结温变化,本发明方法充分考虑了不同物理场应力变化的时间尺度,能够在保证计算精度的同时有效地提高仿真效率,能够为电力电子系统的热设计提供有力的支持。
Description
技术领域
本发明涉及功率半导体模块领域,具体涉及一种基于工况的IGBT模块电-热-流体多场耦合仿真方法。
背景技术
绝缘栅双极型晶体管(IGBT)作为电能转换的“CPU”,在大功率变换器应用中得到越来越广泛的应用,例如风力涡轮机,铁路牵引逆变器和高压直流输电系统中。随着系统功率等级的提升,对电力电子器件的安全性提出了更高的要求。器件内部的温度分布是影响功率半导体器件安全运行的关键参数,过高的结温会使得器件损耗上升,发生热崩溃的风险加剧,严重的会导致器件烧毁。在工程中,设计者往往采用更大容量的器件以提供足够的安全裕量,造成了不必要的成本浪费。因此,在设计中对器件结温进行准确的计算,对保证系统的安全运行和成本控制有着重要的意义。
然而,基于实际工况准确计算结温的难度很大:一方面,器件的损耗与结温具有很强的耦合关系;另一方面,器件的温度分布受散热系统中流体的对流传热影响。因此,单一物理场很难反映器件温度变化的机理,需要考虑模块内部的多物理场耦合效应。
目前的技术存在以下缺陷:
1)现有技术多只考虑单一物理场或某几个物理场对IGBT模块温度的影响。例如应用热—流体耦合仿真对散热器的结构进行优化设计;或应用热场仿真分析模块各层结构的热阻分布。
2)现在技术多针对理想的变换器工况:如在热仿真采用理想的等效换热系数,忽略流体传热的不均匀性;在损耗仿真中采用理想的电路参数,忽略寄生参数、控制策略等实际参数的影响。
3)现在技术多针对整体模块,未考虑模块内部多芯片间的损耗及传热的不均匀性。现有平均化的损耗模型使得结温的仿真结果偏低,不利于器件的安全应用。
4)现在技术多只应用单一时间尺度进行仿真。而不同物理场的时间尺度差异很大:如器件的损耗随着器件的开关周期变化;由于封装材料的热容效应,器件内部的温度不会随着损耗一同剧烈变化;相比较之下,流体流动交换热量的时间尺度更长。采用统一的时间尺度仿真会带来一系列问题:时间尺度过小会造成计算量剧增,严重提高计算成本;时间尺度过大则会忽略细节,无法准确仿真出器件的温度变化。
发明内容
本发明的目的在于提供一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,以解决现有技术存在的问题,本发明能够准确仿真出模块在不同工况下的瞬态结温变化,本发明方法充分考虑了不同物理场应力变化的时间尺度,能够在保证计算精度的同时有效地提高仿真效率,为电力电子系统的热设计提供有力的支持。
为达到上述目的,本发明采用如下技术方案:
一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,包括以下步骤;
步骤一:建立IGBT模块内部芯片级的损耗模型;
步骤二:建立描述IGBT模块内部热传导过程的热模型;
步骤三:建立描述散热系统中流体的热对流过程的共轭传热模型;
步骤四:基于COMSOL Multiphysic和MATLAB编写电-热-流体之间耦合变量的交互程序,实现基于实际工况的功率半导体模块电-热-流体多物理场耦合仿真。
进一步地,步骤一具体为:一方面,利用ANSYS Q3D中提取IGBT模块的杂散电感矩阵,将提取结果导入到LTspice中,结合芯片的Spice模型搭建双脉冲仿真模型,进而得到各芯片开关损耗随负载电流及结温变化的规律;另一方面,在Simulink中搭建电路仿真模型,得到特定控制方式和调制策略下的负载电流和IGBT模块的开关时序;将上述两方面得到的器件信息导入到MATLAB工作区中,即得到IGBT模块内部芯片级的损耗模型。
进一步地,利用ANSYS Q3D中提取IGBT模块的杂散电感矩阵,将提取结果导入到LTspice中,结合芯片的Spice模型搭建双脉冲仿真模型,进而得到各芯片开关损耗随负载电流及结温变化的规律,具体为:
在ANSYS Q3D中建立IGBT模块的杂散电感提取模型,提取频率设置为20MHz,将提取后的杂散电感导入到Ltspice中,结合芯片的Spice模型搭建双脉冲仿真模型,仿真得到不同结温和负载电流下并联各芯片损耗占总损耗的比例,并通过下式拟合:
ki=a+bIC+cTjav+dIC 2+eICTjav+fTjav 2 (1)其中:a~f为拟合系数;IC为集电极电流;Tjav为并联芯片的平均结温;
最终得到的芯片损耗模型,芯片的开通损耗表示为:
Eon_i=ki·Eon_ref+kTj·△Tj_i,i=1,2,3 (2)
其中:Eon_ref为并联芯片总开通损耗;ki不同结温和负载电流下并联各芯片损耗占总损耗的比例,由公式(1)决定;其中kTj为开通损耗的温度系数,表示为:
芯片的关断损耗与杂散电感分布相关性较低,平均化表示:
通过IGBT模块的数据手册,提取单个芯片的正向特性及并联芯片的总开关损耗如下:
Eoff_ref(Tjav,IC)=KT1·IC+KT2 (5)
Vce(Tj_i,IC_i)=KT3·3 IC_i+KT4 (6)
Erec_ref(Tdjav,If)=KT5·If 2+KT6·If+KT7 (7)
Eon_ref(Tjav,IC)=KT8·IC 2+KT9·IC+KT10 (8)
Vf(Tdj_i,If_i)=KT11·(3 If_i)2+KT12·(3 If_i)+KT13 (9)
其中,Eoff_ref为并联芯片总关断损耗;KT1-KT13为与温度相关的拟合系数;IC为并联芯片的集电极电流;Vce为并联芯片的集射极电压;IC_i为单个芯片的集电极电流;Erec_ref为反向恢复的总损耗;If为二极管的总反向恢复电流;Eon_ref为并联芯片的总开通损耗;Vf为二极管的正向电压;If_i为单个二极管的反向恢复电流。
进一步地,在Simulink中搭建电路仿真模型,得到特定控制方式和调制策略下的负载电流和IGBT模块的开关时序,具体为:
在Simulink中仿真得到指定拓扑及控制策略下器件的桥臂电流、母线电压及开关时序,通过公式(10)-(16)计算开关周期内的瞬时损耗Psw;
其中,Econ_i为单个IGBT芯片的传导损耗;Edcon_i为单个反并联二极管芯片的传导损耗;th_on为高压侧芯片的开通时间;th_off高压侧芯片的关断时间;Tji为IGBT芯片温度为Ti时的结温;Tdji为续流二极管温度为Ti时的结温;Tp为开关周期;Vce为并联芯片的集射极电压;Ici为单个芯片的集电极电流;Vf为二极管的正向电压;IC为并联芯片的集电极电流;Esw_i为单个开关周期内单个IGBT芯片的总损耗;Eon_i为单个芯片的开通损耗;Eoff_i为单个芯片的关断损耗;U为集射极电压;Uref为数据手册给出的阻断电压参考值;Kv和Kdv分别为IGBT芯片和续流二极管的电压系数;Edsw_i单个开关周期内单个续流二极管的总损耗;Erec_i为单个芯片的反向恢复损耗;
通过公式(17)计算1ms内的平均损耗,即等效损耗Peq;
其中,Pswi和Pdswi为公式(15)、(16)确定的瞬时损耗,Pdeq_i为二极管的等效损耗;
通过公式(18)计算工频周期内的平均损耗,即平均损耗Pav;
其中,Peq_i和Pdeq_i为公式(17)确定的等效损耗,即得到瞬时损耗、等效损耗、平均损耗的关系,即IGBT模块内部芯片级的损耗模型。
进一步地,步骤二具体为:在Solidworks中建立IGBT模块的第一CAD模型,忽略不在导热路径上的几何体;将第一CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的功率参数P1…Pn及二极管芯片的功率参数Pd1…Pdn,其中n为芯片的编号,同时设置仿真初始时间t0及仿真步长dt;在组件中设置提取各芯片最高结温的域探针Tj1…Tjn及Tdj1…Tdjn;设置插值函数Tc(x,y)及T(x,y,z);物理场选择固体传热,分别并将各芯片设置为热源,大小为在全局定义中设置的功率参数P1…Pn及Pd1…Pdn,将底面设置为固定温度Tc(x,y),其余面设置为绝热,并将温度初值设置为T(x,y,z);添加瞬态研究,将时间步设置为range(t0,dt,t0+dt)。
进一步地,步骤三具体为:在Solidworks中建立IGBT模块及其散热器的第二CAD模型,忽略不在导热路径上的几何体;将第二CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的平均功率参数Pav1…Pavn及二极管芯片的功率参数Pavd1…Pavdn,其中n为芯片的编号;物理场选择固体传热和湍流,多物理场选择非等温流动;将各芯片设置为热源,大小为在全局定义中设置的功率参数Pav1…Pavn及Pavd1…Pavdn,在散热器入口处分别设置边界条件流入和入口,出口处设置边界条件流出和出口,其余面设置为绝热和壁;添加稳态研究。
进一步地,所述的不在导热路径上的几何体包括键合线和端子。
进一步地,首先设定仿真初始温度,在初始温度下通过损耗模型计算各芯片在工频周期下的平均损耗Pav;其次运行共轭传热模型,仿真得到稳态下IGBT模块底面的温度分布,将各点的平面坐标(x,y)及对应的温度Tc,用于构建插值函数Tc(x,y);仿真得到稳态下IGBT模块温度的三维分布,将各点的平面坐标(x,y,z)及对应的温度T,用于构建插值函数T(x,y,z);然后运行瞬态热模型,设定初始仿真时间t0=0,仿真步长dt=1ms,根据稳态仿真出的各芯片结温计算1ms的等效损耗P,刷新插值函数Tc(x,y)及T(x,y,z);第一次瞬态仿真结束后,更新仿真时间为t0=t0+dt,根据仿真出的各芯片结温计算下次瞬态仿真的等效损耗P,同时将仿真得到温度分布作为下次瞬态仿真的初值;当仿真时间达到指定时间后,瞬态仿真结束,计算瞬态仿真过程中各芯片的平均损耗Pav_out,比较Pav_out和Pav值的大小,当二者的差值大于1%Pav_out时,用Pav_out代替Pav重新计算Tc(x,y),迭代计算直至平均损耗收敛。
与现有技术相比,本发明具有以下有益的技术效果:
本发明将所有能影响IGBT结温的物理场及其间的耦合关系都考虑在内,使得计算结果更接近实际;本发明基于实际的变换器工况,充分考虑散热性能的不均匀性及电路参数的非理想特性,相比于现有技术更能满足实际工程设计的需要;本发明针对并联芯片的损耗分布,建立了不均匀损耗的仿真模型,进而能够预测出模块内部温度最高的芯片以进行安全性评估,对大功率的多芯片IGBT模块在实际工况下的安全运行有着重要的意义;本发明针对不同物理场的特点,采用变步长的混合仿真策略,解决了单一步长下仿真精度和仿真效率不能兼得的矛盾。
具体如下:
(1)仿真精度高:本发明将所有影响结温的因素全都考虑在内,包括芯片损耗受温度的影响;流体对壳温不均匀分布的影响等,因此仿真的精度与实际情况吻合度更高。
(2)仿真效率高:本发明根据不同物理场应力的变化尺度,采用多时间尺度的混合仿真策略,具体为损耗计算根据电路仿真的输出采用ns级别的步长;热传导采用1ms的仿真步长;由于壳温对损耗中的交流成分不敏感,因此对热-流体共轭传热过程采用稳态仿真。混合仿真策略有利于在保证仿真精度的同时极大地提升仿真效率。
(3)适用性广:本发明的应用不受限于具体电路拓扑、控制策略、器件型号及散热方式的限制,能够适用于不同实际工况的仿真计算。
附图说明
图1为本发明实施步骤示意图;
图2(a)为IGBT模块内部示意图;
图2(b)为Q3D仿真示意图;
图2(c)为双脉冲仿真示意图;
图3为芯片级开通损耗模型示意图;
图4为瞬时损耗、等效损耗、平均损耗示意图;
图5为迭代计算过程中的变量的传递情况;
图6为瞬态结温仿真结果,其中(a)为IGBT芯片,(b)为二极管芯片;
图7为瞬态三维结温分布仿真结果;
图8为仿真与实验的对照图。
具体实施方式
下面结合具体实例对本发明做进一步详细说明,所述内容均为对本发明的解释,但这些说明不能被理解为限制了本发明的范围,本发明的保护范围由随附的权利要求书限定,任何在本发明权利要求基础上的改动就是本发明的保护范围。
一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其中“电”指的是模块内部芯片级的损耗模型;“热”指描述模块封装中的热传导过程的有限元模型;“流体”指描述散热系统中流体的热对流过程的有限元模型。
本发明的实施方式如附图1所示,具体步骤为:
步骤一:在ANSYS Q3D中建立IGBT模块的杂散电感提取模型,提取频率设置为20MHz,如附图2(b)。将提取后的杂散电感导入到Ltspice中,结合芯片的Spice模型搭建双脉冲仿真模型,如附图2(c)所示。仿真得到不同结温和负载电流下并联各芯片损耗占总损耗的比例,并通过下式拟合
ki=a+bIC+cTjav+dIC 2+eICTjav+fTjav 2 (1)其中:a~f为拟合系数;IC为集电极电流;Tjav为并联芯片的平均结温。
最终得到的芯片损耗模型如附图3所示。芯片的开通损耗可以表示为:
Eon_i=ki·Eon_ref+kTj·△Tj_i,i=1,2,3 (2)
其中:Eon_ref为并联芯片总开通损耗;ki不同结温和负载电流下并联各芯片损耗占总损耗的比例,由公式(1)决定;kTj为开通损耗的温度系数,表示为:
芯片的关断损耗与杂散电感分布相关性较低,可以平均化表示:
通过IGBT模块的数据手册,可以提取单个芯片的正向特性及并联芯片的总开关损耗如下:
Eoff_ref(Tjav,IC)=KT1·IC+KT2 (5)
Vce(Tj_i,IC_i)=KT3·3 IC_i+KT4 (6)
Erec_ref(Tdjav,If)=KT5·If 2+KT6·If+KT7 (7)
Eon_ref(Tjav,IC)=KT8·IC 2+KT9·IC+KT10 (8)
Vf(Tdj_i,If_i)=KT11·(3 If_i)2+KT12·(3 If_i)+KT13 (9)
其中Eoff_ref为并联芯片总关断损耗;KT1-KT13为与温度相关的拟合系数;IC为并联芯片的集电极电流;Vce为并联芯片的集射极电压;IC_i为单个芯片的集电极电流;Erec_ref为反向恢复的总损耗;If为二极管的总反向恢复电流;Eon_ref为并联芯片的总开通损耗;Vf为二极管的正向电压;If_i为单个二极管的反向恢复电流;
步骤二:在Simulink中仿真得到指定拓扑及控制策略下器件的桥臂电流、母线电压及驱动信号,并导入到MATLAB工作区。通过公式(10)-(16)计算开关周期内的瞬时损耗Psw,Tp为开关周期。
其中,Econ_i为单个IGBT芯片的传导损耗;Edconi为单个反并联二极管芯片的传导损耗;th_on为高压侧芯片的开通时间;th_off高压侧芯片的关断时间;Tji为IGBT芯片温度为Ti时的结温;Tdji为续流二极管温度为Ti时的结温;Tp为开关周期;Vce为并联芯片的集射极电压;Ici为单个芯片的集电极电流;Vf为二极管的正向电压;IC为并联芯片的集电极电流;Esw_i为单个开关周期内单个IGBT芯片的总损耗;Eon_i为单个芯片的开通损耗;Eoff_i为单个芯片的关断损耗;U为集射极电压;Uref为数据手册给出的阻断电压参考值;Kv和Kdv分别为IGBT芯片和续流二极管的电压系数;Edsw_i单个开关周期内单个续流二极管的总损耗;Erec_i为单个芯片的反向恢复损耗。
通过公式(17)计算1ms内的平均损耗,即等效损耗Peq。
其中,Psw_i和Pdsw_i为公式(15)、(16)确定的瞬时损耗。
通过公式(18)计算工频周期内的平均损耗,即平均损耗Pav。
其中,Peq_i和Pdeq_i为公式(17)确定的等效损耗。
瞬时损耗、等效损耗、平均损耗的关系如附图4所示。
步骤三:建立描述IGBT模块内部热传导过程的热模型:在Solidworks中建立IGBT模块的CAD模型,忽略键合线和端子等不在导热路径上的几何体;将CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的功率参数P1…Pn及二极管芯片的功率参数Pd1…Pdn,其中n为芯片的编号,同时设置仿真初始时间t0及仿真步长dt;在组件中设置提取各芯片最高结温的域探针Tj1…Tjn及Tdj1…Tdjn;设置插值函数Tc(x,y),数据源为文件,文件名设置为“C:\临时文件路径\TC_data.txt”,数据格式选择电子表格,变元数设置为2;设置插值函数T(x,y,z),数据源为文件,文件名设置为“C:\临时文件路径\T_data.txt”,数据格式选择电子表格,变元数设置为3;物理场选择固体传热,并将各芯片设置为热源,大小为在全局定义中设置的功率参数P1…Pn及Pd1…Pdn,将底面设置为固定温度Tc(x,y),其余面设置为绝热,并将温度初值设置为T(x,y,z);添加瞬态研究,将时间步设置为range(t0,dt,t0+dt)。
步骤四:建立描述散热系统中流体的热对流过程的共轭传热模型:在Solidworks中建立IGBT模块及其散热器的CAD模型,忽略键合线和端子等不在导热路径上的几何体;将CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的平均功率参数Pav1…Pavn及二极管芯片的功率参数Pavd1…Pavdn,其中n为芯片的编号;物理场选择固体传热和湍流,多物理场选择非等温流动;将各芯片设置为热源,大小为在全局定义中设置的功率参数Pav1…Pavn及Pavd1…Pavdn,在散热器入口处分别设置边界条件“流入”和“入口”,出口处设置边界条件“流出”和“出口”,其余面设置为“绝热”和“壁”;添加稳态研究。
步骤五:基于COMSOL Multiphysic with MATLAB编写“电-热-流体”之间耦合变量的交互程序。首先设定仿真初始温度,在初始温度下通过损耗模型计算各芯片在工频周期下的平均损耗Pav;其次运行共轭传热模型,仿真得到稳态下IGBT模块底面的温度分布,将各点的平面坐标(x,y)及对应的温度Tc,导入到路径为“C:\临时文件路径\TC_data.txt”的记事本文件中;仿真得到稳态下IGBT模块温度的三维分布,将各点的平面坐标(x,y,z)及对应的温度T,导入到路径为“C:\临时文件路径\T_data.txt”的文件中;然后运行瞬态热模型,设定初始仿真时间t0=0,仿真步长dt=1ms,根据稳态仿真出的各芯片结温计算1ms的等效损耗P,刷新插值函数Tc(x,y)及T(x,y,z);第一次瞬态仿真结束后,更新仿真时间为t0=t0+dt,根据仿真出的各芯片结温计算下次瞬态仿真的等效损耗P,同时将仿真得到温度分布作为下次瞬态仿真的初值;当仿真时间达到指定时间后,瞬态仿真结束,计算瞬态仿真过程中各芯片的平均损耗Pav_out,比较Pav_out和Pav值的大小,当二者相差过大时,用Pav_out代替Pav重新计算Tc(x,y),迭代计算直至损耗收敛。迭代计算过程中的变量的传递情况如附图5所示。最终的仿真结果如附图6和图7所示,可以看出该仿真方法可以计算出瞬时损耗下各芯片的三维温度分布,进而得到IGBT模块在实际工况中的温度薄弱点,能够帮助IGBT模块作进一步散热及可靠性优化。仿真与实验的对照如附图8所示,实验结果能够充分验证该仿真方法的准确性。
以上结合具体实例对本发明的具体实施方式对本发明作了进一步详细说明,所述内容均为对本发明的解释但这些说明不能被理解为限制了本发明的范围,本发明的保护范围由随附的权利要求书限定,任何在本发明权利要求基础上的改动就是本发明的保护范围。
Claims (5)
1.一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其特征在于,包括以下步骤;
步骤一:建立IGBT模块内部芯片级的损耗模型;
具体为:一方面,利用ANSYS Q3D中提取IGBT模块的杂散电感矩阵,将提取结果导入到LTspice中,结合芯片的Spice模型搭建双脉冲仿真模型,进而得到各芯片开关损耗随负载电流及结温变化的规律;另一方面,在Simulink中搭建电路仿真模型,得到特定控制方式和调制策略下的负载电流和IGBT模块的开关时序;将上述两方面得到的器件信息导入到MATLAB工作区中,即得到IGBT模块内部芯片级的损耗模型;
其中,利用ANSYS Q3D中提取IGBT模块的杂散电感矩阵,将提取结果导入到LTspice中,结合芯片的Spice模型搭建双脉冲仿真模型,进而得到各芯片开关损耗随负载电流及结温变化的规律,具体为:
在ANSYS Q3D中建立IGBT模块的杂散电感提取模型,提取频率设置为20MHz,将提取后的杂散电感导入到Ltspice中,结合芯片的Spice模型搭建双脉冲仿真模型,仿真得到不同结温和负载电流下并联各芯片损耗占总损耗的比例,并通过下式拟合:
ki=a+bIC+cTjav+dIC 2+eICTjav+fTjav 2 (1)
其中:a~f为拟合系数;IC为集电极电流;Tjav为并联芯片的平均结温;
最终得到的芯片损耗模型,芯片的开通损耗表示为:
Eon_i=ki·Eon_ref+kTj·△Tj_i,i=1,2,3 (2)
其中:Eon_ref为并联芯片总开通损耗;ki不同结温和负载电流下并联各芯片损耗占总损耗的比例,由公式(1)决定;其中kTj为开通损耗的温度系数,表示为:
芯片的关断损耗与杂散电感分布相关性较低,平均化表示:
通过IGBT模块的数据手册,提取单个芯片的正向特性及并联芯片的总开关损耗如下:
Eoff_ref(Tjav,IC)=KT1·IC+KT2 (5)
Vce(Tj_i,IC_i)=KT3·3IC_i+KT4 (6)
Erec_ref(Tdjav,If)=KT5·If 2+KT6·If+KT7 (7)
Eon_ref(Tjav,IC)=KT8·IC 2+KT9·IC+KT10 (8)
Vf(Tdj_i,If_i)=KT11·(3 If_i)2+KT12·(3 If_i)+KT13 (9)
其中,Eoff_ref为并联芯片总关断损耗;KT1-KT13为与温度相关的拟合系数;IC为并联芯片的集电极电流;Vce为并联芯片的集射极电压;IC_i为单个芯片的集电极电流;Erec_ref为反向恢复的总损耗;If为二极管的总反向恢复电流;Eon_ref为并联芯片的总开通损耗;Vf为二极管的正向电压;If_i为单个二极管的反向恢复电流;
在Simulink中搭建电路仿真模型,得到特定控制方式和调制策略下的负载电流和IGBT模块的开关时序,具体为:
在Simulink中仿真得到指定拓扑及控制策略下器件的桥臂电流、母线电压及开关时序,通过公式(10)-(16)计算开关周期内的瞬时损耗Psw;
其中,Econ_i为单个IGBT芯片的传导损耗;Edcon_i为单个反并联二极管芯片的传导损耗;th_on为高压侧芯片的开通时间;th_off高压侧芯片的关断时间;Tji为IGBT芯片温度为Ti时的结温;Tdji为续流二极管温度为Ti时的结温;Tp为开关周期;Vce为并联芯片的集射极电压;Ici为单个芯片的集电极电流;Vf为二极管的正向电压;IC为并联芯片的集电极电流;Esw_i为单个开关周期内单个IGBT芯片的总损耗;Eon_i为单个芯片的开通损耗;Eoff_i为单个芯片的关断损耗;U为集射极电压;Uref为数据手册给出的阻断电压参考值;Kv和Kdv分别为IGBT芯片和续流二极管的电压系数;Edsw_i单个开关周期内单个续流二极管的总损耗;Erec_i为单个芯片的反向恢复损耗;
通过公式(17)计算1ms内的平均损耗,即等效损耗Peq;
其中,Psw_i和Pdsw_i为公式(15)、(16)确定的瞬时损耗,Pdeq_i为二极管的等效损耗;
通过公式(18)计算工频周期内的平均损耗,即平均损耗Pav;
其中,Peq_i和Pdeq_i为公式(17)确定的等效损耗,即得到瞬时损耗、等效损耗、平均损耗的关系,即IGBT模块内部芯片级的损耗模型;
步骤二:建立描述IGBT模块内部热传导过程的热模型;
步骤三:建立描述散热系统中流体的热对流过程的共轭传热模型;
步骤四:基于COMSOL Multiphysic和MATLAB编写电-热-流体之间耦合变量的交互程序,实现基于实际工况的功率半导体模块电-热-流体多物理场耦合仿真。
2.根据权利要求1所述的一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其特征在于,步骤二具体为:在Solidworks中建立IGBT模块的第一CAD模型,忽略不在导热路径上的几何体;将第一CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的功率参数P1…Pn及二极管芯片的功率参数Pd1…Pdn,其中n为芯片的编号,同时设置仿真初始时间t0及仿真步长dt;在组件中设置提取各芯片最高结温的域探针Tj1…Tjn及Tdj1…Tdjn;设置插值函数Tc(x,y)及T(x,y,z);物理场选择固体传热,分别并将各芯片设置为热源,大小为在全局定义中设置的功率参数P1…Pn及Pd1…Pdn,将底面设置为固定温度Tc(x,y),其余面设置为绝热,并将温度初值设置为T(x,y,z);添加瞬态研究,将时间步设置为range(t0,dt,t0+dt)。
3.根据权利要求2所述的一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其特征在于,步骤三具体为:在Solidworks中建立IGBT模块及其散热器的第二CAD模型,忽略不在导热路径上的几何体;将第二CAD模型导入到COMSOL有限元仿真模型中进行网格划分;在全局定义中设置各IGBT芯片的平均功率参数Pav1…Pavn及二极管芯片的功率参数Pavd1…Pavdn,其中n为芯片的编号;物理场选择固体传热和湍流,多物理场选择非等温流动;将各芯片设置为热源,大小为在全局定义中设置的功率参数Pav1…Pavn及Pavd1…Pavdn,在散热器入口处分别设置边界条件流入和入口,出口处设置边界条件流出和出口,其余面设置为绝热和壁;添加稳态研究。
4.根据权利要求3所述的一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其特征在于,所述的不在导热路径上的几何体包括键合线和端子。
5.根据权利要求2所述的一种基于工况的IGBT模块电-热-流体多场耦合仿真方法,其特征在于,首先设定仿真初始温度,在初始温度下通过损耗模型计算各芯片在工频周期下的平均损耗Pav;其次运行共轭传热模型,仿真得到稳态下IGBT模块底面的温度分布,将各点的平面坐标(x,y)及对应的温度Tc,用于构建插值函数Tc(x,y);仿真得到稳态下IGBT模块温度的三维分布,将各点的平面坐标(x,y,z)及对应的温度T,用于构建插值函数T(x,y,z);然后运行瞬态热模型,设定初始仿真时间t0=0,仿真步长dt=1ms,根据稳态仿真出的各芯片结温计算1ms的等效损耗P,刷新插值函数Tc(x,y)及T(x,y,z);第一次瞬态仿真结束后,更新仿真时间为t0=t0+dt,根据仿真出的各芯片结温计算下次瞬态仿真的等效损耗P,同时将仿真得到温度分布作为下次瞬态仿真的初值;当仿真时间达到指定时间后,瞬态仿真结束,计算瞬态仿真过程中各芯片的平均损耗Pav_out,比较Pav_out和Pav值的大小,当二者的差值大于1%Pav_out时,用Pav_out代替Pav重新计算Tc(x,y),迭代计算直至平均损耗收敛。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110728570.2A CN113435090B (zh) | 2021-06-29 | 2021-06-29 | 一种基于工况的igbt模块电-热-流体多场耦合仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110728570.2A CN113435090B (zh) | 2021-06-29 | 2021-06-29 | 一种基于工况的igbt模块电-热-流体多场耦合仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113435090A CN113435090A (zh) | 2021-09-24 |
CN113435090B true CN113435090B (zh) | 2022-10-25 |
Family
ID=77757677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110728570.2A Active CN113435090B (zh) | 2021-06-29 | 2021-06-29 | 一种基于工况的igbt模块电-热-流体多场耦合仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113435090B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113985122A (zh) * | 2021-11-01 | 2022-01-28 | 苏州亿马半导体科技有限公司 | 一种基于SolidWorks Flow Simulation的SiC功率分析方法 |
CN114564850B (zh) * | 2022-04-24 | 2022-08-02 | 成都复锦功率半导体技术发展有限公司 | 一种电力电子系统电热联合仿真方法、系统及终端 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015078367A1 (zh) * | 2013-11-27 | 2015-06-04 | 国家电网公司 | 一种基于模块化多电平电压源型换流器的损耗确定方法 |
CN105825019A (zh) * | 2016-03-22 | 2016-08-03 | 三峡大学 | 一种绝缘栅双极晶体管igbt模块温度求解算法 |
CN106326532A (zh) * | 2016-08-12 | 2017-01-11 | 中国人民解放军海军工程大学 | Igbt传热多时间尺度模型建模方法 |
CN107453615A (zh) * | 2016-05-31 | 2017-12-08 | 西门子公司 | 模块化多电平变换器及电力电子变压器 |
CN108664675A (zh) * | 2017-03-31 | 2018-10-16 | 上海北京大学微电子研究院 | 一种电热耦合瞬态仿真模型的建立方法 |
CN109918700A (zh) * | 2018-12-04 | 2019-06-21 | 中国人民解放军海军工程大学 | Igbt多时间尺度结温预测模型建模方法 |
CN110245414A (zh) * | 2019-06-11 | 2019-09-17 | 南方电网科学研究院有限责任公司 | 一种压接式igbt模块多物理场耦合仿真方法 |
CN110765601A (zh) * | 2019-10-12 | 2020-02-07 | 北京北方华德尼奥普兰客车股份有限公司 | 一种基于igbt热电耦合模型的igbt结温估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109871591A (zh) * | 2019-01-24 | 2019-06-11 | 武汉大学 | 一种igbt功率模块在线估算结温的方法 |
-
2021
- 2021-06-29 CN CN202110728570.2A patent/CN113435090B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015078367A1 (zh) * | 2013-11-27 | 2015-06-04 | 国家电网公司 | 一种基于模块化多电平电压源型换流器的损耗确定方法 |
CN105825019A (zh) * | 2016-03-22 | 2016-08-03 | 三峡大学 | 一种绝缘栅双极晶体管igbt模块温度求解算法 |
CN107453615A (zh) * | 2016-05-31 | 2017-12-08 | 西门子公司 | 模块化多电平变换器及电力电子变压器 |
CN106326532A (zh) * | 2016-08-12 | 2017-01-11 | 中国人民解放军海军工程大学 | Igbt传热多时间尺度模型建模方法 |
CN108664675A (zh) * | 2017-03-31 | 2018-10-16 | 上海北京大学微电子研究院 | 一种电热耦合瞬态仿真模型的建立方法 |
CN109918700A (zh) * | 2018-12-04 | 2019-06-21 | 中国人民解放军海军工程大学 | Igbt多时间尺度结温预测模型建模方法 |
CN110245414A (zh) * | 2019-06-11 | 2019-09-17 | 南方电网科学研究院有限责任公司 | 一种压接式igbt模块多物理场耦合仿真方法 |
CN110765601A (zh) * | 2019-10-12 | 2020-02-07 | 北京北方华德尼奥普兰客车股份有限公司 | 一种基于igbt热电耦合模型的igbt结温估计方法 |
Non-Patent Citations (4)
Title |
---|
Electrothermal Modeling of the Insulated Gate Bipolar Transistor;Wasma Hanini,Moez Ayadi;《Journal of Control, Automation and Electrical Systems》;20210103;第507–521页 * |
IGBT模块开关损耗计算方法综述;李志刚等;《电子技术应用》;20160112(第01期);第10-14页 * |
一种IGBT模块的新型建模与瞬态仿真方法研究;沈天浩等;《机电工程》;20180720(第07期);第111-116页 * |
压接型IGBT 在MMC 系统中的电热耦合仿真;侯婷,苟浪中,李岩等;《南方电网技术》;20200531;第14卷(第5期);第28-36页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113435090A (zh) | 2021-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113435090B (zh) | 一种基于工况的igbt模块电-热-流体多场耦合仿真方法 | |
Xu et al. | Power loss and junction temperature analysis of power semiconductor devices | |
CN103324843A (zh) | 一种适用于不同子模块类型的mmc阀损耗计算方法 | |
Wang et al. | A transient 3-D thermal modeling method for IGBT modules considering uneven power losses and cooling conditions | |
CN115629553A (zh) | 一种高时空分辨率igbt模块电热应力计算方法 | |
Li et al. | A 1200 V/200 a half-bridge power module based on Si IGBT/SiC MOSFET hybrid switch | |
Li et al. | Using multi time-scale electro-thermal simulation approach to evaluate SiC-MOSFET power converter in virtual prototyping design tool | |
CN112560318A (zh) | 一种不同工况下igbt器件各物理层温度监测仿真方法 | |
Reichl et al. | Design optimization of hybrid-switch soft-switching inverters using multiscale electrothermal simulation | |
Zhong et al. | Design & analysis of a novel IGBT package with double-sided cooling | |
Davoudi et al. | Averaged-switch modeling of fourth-order PWM DC–DC converters considering conduction losses in discontinuous mode | |
CN115688509A (zh) | 一种多芯片SiC MOSFET并联功率半导体模块的电热联合仿真方法 | |
Wang et al. | Fast and accurate temperature estimation of three-level IGBT converter based on 3-D coupled thermal model | |
Peng et al. | Analytical model for predicting the junction temperature of chips considering the internal electrothermal coupling inside SiC metal–oxide–semiconductor field‐effect transistor modules | |
CN113536627B (zh) | 一种多芯片igbt模块热安全运行域刻画方法 | |
CN203504514U (zh) | Mos芯片并联均流集成开关及其封装模块 | |
CN107863781B (zh) | 一种模块化多电平换流阀损耗确定方法及简化方法 | |
Sartori et al. | Power density comparative analysis concerning to three transistor technologies applied to a CCM PFC BOOST converter using optimization techniques | |
Cheng et al. | Dynamic modeling framework for evaluating electromagnetic-electro-thermal behavior of power conversion system during load operation | |
Bulut et al. | Thermal modeling for power MOSFETs in DC/DC applications | |
Wang et al. | Heat dissipation system modeling of three-level converter considering thermal coupling effects | |
Xin et al. | Electro-thermal based junction temperature estimation model and thermal performance analysis for IGBT module | |
Nejadpak et al. | Effects of different switching algorithms on the thermal behavior of IGBT modules under pulse-load conditions | |
Ishiko et al. | A compact calculation method for dynamic electro-thermal behavior of IGBTs in PWM inverters | |
Fang | Design and implementation of a thermal hardware-in-the-loop (THIL) platform for evaluating physical cooling systems in real-time simulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |