CN110988641A - 一种igbt模块结温在线估计电路系统及方法 - Google Patents

一种igbt模块结温在线估计电路系统及方法 Download PDF

Info

Publication number
CN110988641A
CN110988641A CN201911240970.8A CN201911240970A CN110988641A CN 110988641 A CN110988641 A CN 110988641A CN 201911240970 A CN201911240970 A CN 201911240970A CN 110988641 A CN110988641 A CN 110988641A
Authority
CN
China
Prior art keywords
igbt
state
junction temperature
voltage
module
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.)
Pending
Application number
CN201911240970.8A
Other languages
English (en)
Inventor
何怡刚
王晨苑
时国龙
许水清
李猎
杜博伦
曾昭瑢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201911240970.8A priority Critical patent/CN110988641A/zh
Publication of CN110988641A publication Critical patent/CN110988641A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/26Testing of individual semiconductor devices
    • G01R31/2607Circuits therefor
    • G01R31/2608Circuits therefor for testing bipolar transistors
    • G01R31/2619Circuits therefor for testing bipolar transistors for measuring thermal properties thereof

Landscapes

  • Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Individual Semiconductor Devices (AREA)

Abstract

本发明公开了一种IGBT模块结温在线估计电路系统及方法,包括以下步骤:步骤1、在线测量IGBT功率模块的IGBT通态压降VCE(ON),并建立IGBT结温Tj查找表;步骤2、计算IGBT功率模块中IGBT和二极管的功率损耗;步骤3、采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型;步骤4、根据IGBT等效热网络模型以结温Tj和IGBT、二极管功率损耗为状态变量建立状态空间模型,并更新状态空间模型;步骤5、根据状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure DDA0002306216510000011
本发明融合了模型估计和在线测量的优点,克服了通态压降VCE(ON)测量不准确的缺点,保持了信号的连续性,降低了噪声干扰,且估计精度高,能实时、非干预式的估计IGBT模块结温。

Description

一种IGBT模块结温在线估计电路系统及方法
技术领域
本发明涉及电力电子关键器件可靠性技术领域,尤其涉及一种IGBT模块结温在线估计电路系统及方法。
背景技术
目前,以IGBT模块为核心的电力电子装置广泛应用于电机控制、电动汽车、新能源、智能电网等领域,在实际运行中,器件长期承受的复杂的功率、温度循环可能导致巨大的结温以及结温波动。由于半导体器件中各层材料热膨胀系数不一致,结温的波动引起各层膨胀程度不同,使得各层之间产生热应力,导致焊料层疲劳老化或键合线脱落,从而造成器件电、热参数等外特性变化,最终导致器件老化失效。因此,对IGBT进行结温估计十分重要。精确的结温估计在IGBT模块结温控制、健康管理以及可靠性中具有重要的意义。
结温的测量已有很多方法,直接测量主要采用的是红外测温仪和光导纤维,但对于封装好的模块并不适用;采用集成传感器在设计和制造过程中需要特别考虑,并且需要额外的外部引脚,可能会增加制造成本并带来新的可靠性问题;RC等效热模型也为结温估算提供了一种简单有效的方法,但该方法的准确度和精度依赖于模型的准确度和模型参数精确度,并且没有考虑热路径的退化即焊料层的疲劳累积效应。而温敏电参数法测结温具有响应速度快、测量原理简单和不需要破坏结构的优点而被广泛应用。但由于温敏电气参数对结温的灵敏度较低,测量具有一定误差,并且受电力电子装备负载条件影响,因此直接利用温敏电气参数估计结温具有一定难度。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供一种IGBT模块结温在线估计电路系统及方法。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种IGBT模块结温在线估计电路系统,用于在线测量IGBT 功率模块的结温数据,该电路系统包括:功率变换器、通态压降VCE(ON)在线测量电路、功率损耗计算模块、结温查找表模块、计及焊料层疲劳累积效应热网络模块、自适应卡尔曼滤波模块;其中:
功率变换器,其核心器件为IGBT功率模块,用于获取负载电流IL
通态压降VCE(ON)在线测量电路,用于在线测量IGBT通态压降VCE(ON)
功率损耗计算模块,用于根据负载电流IL和IGBT通态压降VCE(ON),计算IGBT功率模块中IGBT和二极管的功率损耗;
结温温度查表模块,用于建立结温Tj的查找表,通过测量获取不同负载电流IL下的通态压降VCE(ON)和结温Tj,并选取一定数量的点进行拟合,得到结温Tj与负载电流IL、通态压降VCE(ON)之间的关系;
计及焊料层疲劳累积效应热网络模块,用于建立IGBT等效热网络模型,并根据IGBT等效热网络模型以及结温Tj和IGBT、二极管功率损耗为状态变量,建立状态空间模型;采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型,进而实现对状态空间模型的更新;
自适应卡尔曼滤波模块,用于根据状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure BDA0002306216490000021
进一步地,本发明的通态压降VCE(ON)在线测量电路包括分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4、运算放大器U1、反馈电阻 R5、ESD二极管、下降沿检测器、单稳态多谐振荡器、电容器、运算放大器 U2、模数转换器ADC;其中:
分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;
运算放大器U1是差分放大器,当IGBT器件处于关断状态时,通态压降 VCE(ON)不需要测量,运算放大器U1输出饱和,Vout=VH_U1,其中Vout为运放 U1的输出,VH_U1为运算放大器U1的饱和电压;当IGBT处于导通状态时,运算放大器U1从饱和状态恢复到VCE(ON)值,即Vout=VCE(ON);
反馈电阻R5用于限制运算放大器U1的放大倍数,稳定工作点;
ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;
并联于分压电阻器的电容器用于改善分压器的频率响应;
下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,并触发单稳态多谐振荡器以启用采样电路,延迟时间能通过调节多振子元器件的值来调整,以为VCE提供足够的稳定时间;
运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;
电容器用于隔离所述运算放大U2的高频信号;
模数转换器用于在测量端测量通态压降VCE(ON)
本发明提供一种IGBT模块结温在线估计方法,该方法包括以下步骤:
步骤1、采用IGBT模块结温在线估计电路系统测量IGBT功率模块的 IGBT通态压降VCE(ON),并建立IGBT结温Tj查找表;
步骤2、计算IGBT功率模块中IGBT和二极管的功率损耗;
步骤3、建立IGBT等效热网络模型,并采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型;
步骤4、根据IGBT等效热网络模型以步骤1中结温Tj和步骤2中IGBT、二极管功率损耗为状态变量建立状态空间模型,并按步骤3更新IGBT等效热网络模型,进而更新状态空间模型;
步骤5、根据步骤4中的状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure BDA0002306216490000031
进一步地,本发明的步骤1中在线测量IGBT通态压降VCE(ON)的具体方法如下:
分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;
运算放大器U1是差分放大器,当IGBT器件处于关断状态时,通态压降 VCE(ON)不需要测量,运算放大器U1输出饱和,Vout=VH_U1,其中Vout为运放 U1的输出,VH_U1为运算放大器U1的饱和电压;当IGBT处于导通状态时,运算放大器U1从饱和状态恢复到VCE(ON)值,即Vout=VCE(ON)
反馈电阻R5用于限制运算放大器U1的放大倍数,稳定工作点;
ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;
并联于分压电阻器的电容器用于改善分压器的频率响应;
下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,并触发单稳态多谐振荡器以启用采样电路,延迟时间能通过调节多振子元器件的值来调整,以为VCE提供足够的稳定时间;
运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;
电容器用于隔离所述运算放大U2的高频信号;
模数转换器用于在测量端测量通态压降VCE(ON)
进一步地,本发明的步骤1中建立IGBT结温Tj查找表的具体方法如下:
通过所测通态压降VCE(ON)和负载电流IL建立IGBT结温Tj的查找表,通过测量获取不同负载电流IL下的通态压降VCE(ON)和结温Tj,并选取一定数量的点应用Matlab软件进行拟合,得到Tj=(VCE(ON),IL)。
进一步地,本发明的步骤2中计算IGBT功率模块IGBT和二极管的功率损耗的方法为:
根据通态压降、负载电流、占空比、IGBT数据手册开关损耗曲线计算IGBT 和二极管功率损耗。
进一步地,本发明的步骤3中计及焊料层疲劳累积效应的热网络更新策略为:
建立四阶IGBT热网络模型,采用阶段性适时更新热网络方式,选取热阻增大10%作为热网络更新系数,以热阻增大50%作为器件完全失效判据,将模块老化过程离散为Rth、1.1×Rth、1.2×Rth、1.3×Rth、1.4×Rth五个阶段。
进一步地,本发明的步骤4中状态空间模型的建立方法如下:
步骤41、自热建模:根据IGBT热网络模型,选择选择IGBT和二极管功率损耗以及结温作为状态变量,列写状态空间方程:
状态方程:
Figure BDA0002306216490000041
输出方程:Tj(t)=Cx(t)+Du(t)
Figure BDA0002306216490000051
C=[1 1 1… 1];D=[0 1];u(t)=[PD(t),Ta]
其中,An×n为系统矩阵,Bn×2为输入矩阵,Cl×n为输出矩阵,D1×2为前馈矩阵,状态向量
Figure BDA0002306216490000052
表示RC单元间的温差,u(t)为系统输入向量,其中PD(t)为功率损耗,Ta为环境温度;
步骤42、考虑二极管交叉耦合效应:将自热阻抗的n阶状态空间模型与交叉耦合阻抗的m阶模型相结合,构造出同时表示两种阻抗的扩展模型,得到完整的IGBT结温状态空间模型:
Figure BDA0002306216490000053
Figure BDA0002306216490000054
其中,xs1,…,xsn为自热阻抗状态,ps1,…,psn和ks1,…,ksn分别为自热阻抗的极点和留数,类似地,xc1,…,xcm为交叉耦合阻抗的状态,pc1,…,pcm和kc1,…,kcm分别为自热阻抗的极点和留数,Tj为IGBT结温,Ta为环境温度,PIGBT和PD分别为IGBT和二极管的功率损耗;
步骤43、模型参数辨识:利用仿真软件中拟合工具箱拟合器件的瞬态热阻抗曲线求解参数值得健康状态下的模型参数,再根据步骤3更新步骤42状态空间模型中对应参数。
进一步地,本发明的步骤43中采用的仿真软件为MATLAB或者Saber。
进一步地,本发明的步骤5中建立自适应卡尔曼滤波器系统,利用步骤1 中IGBT结温Tj、步骤3中IGBT和二极管功率损耗为滤波器输入,计算得到结温Tj估计值
Figure BDA0002306216490000061
本发明产生的有益效果是:本发明的IGBT模块结温在线估计电路系统及方法,(1)通态压降VCE(ON)的测量是在大电流条件下,基于模块外部测量,不需要破坏模块完整性,不需要修改栅极驱动电路或干扰IGBT的正常开关周期;(2)计及焊料层疲劳累积效应,提出自适应热模型来量化热阻的变化,并相应地更新模型参数,使得结温Tj估计精度高;(3)在结温Tj的测量中引入卡尔曼滤波器来提高测量精度,消除噪声和干扰,对结温Tj进行最优估计,具有预测速度快、准确性高等优点。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明实施例的IGBT模块结温在线估计方法的流程图;
图2为本发明实施例的IGBT通态VCE(ON)的测量电路;
图3为本发明实施例的Foster热网络等效电路。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明实施例的IGBT模块结温在线估计电路系统,用于在线测量IGBT功率模块的结温数据,该电路系统包括:功率变换器、通态压降 VCE(ON)在线测量电路、功率损耗计算模块、结温查找表模块、计及焊料层疲劳累积效应热网络模块、自适应卡尔曼滤波模块;其中:
功率变换器,其核心器件为IGBT功率模块,用于获取负载电流IL
通态压降VCE(ON)在线测量电路,用于在线测量IGBT通态压降VCE(ON)
功率损耗计算模块,用于根据负载电流IL和IGBT通态压降VCE(ON),计算IGBT功率模块中IGBT和二极管的功率损耗;
结温查找表模块,用于建立结温Tj的查找表,通过测量获取不同负载电流IL下的通态压降VCE(ON)和结温Tj,并选取一定数量的点进行拟合,得到结温Tj与负载电流IL、通态压降VCE(ON)之间的关系;
计及焊料层疲劳累积效应热网络模块,用于建立IGBT等效热网络模型,并根据IGBT等效热网络模型以及结温Tj和IGBT、二极管功率损耗为状态变量,建立状态空间模型;采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型,进而实现对状态空间模型的更新;
自适应卡尔曼滤波模块,用于根据状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure BDA0002306216490000072
如图2所示,通态压降VCE(ON)在线测量电路包括分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4、运算放大器U1、反馈电阻R5、ESD 二极管、下降沿检测器、单稳态多谐振荡器、电容器、运算放大器U2、模数转换器ADC;其中:
分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;
运算放大器U1是差分放大器,当IGBT器件处于关断状态时,通态压降 VCE(ON)不需要测量,运算放大器U1输出饱和,Vout=VH_U1,其中Vout为运放 U1的输出,VH_U1为运算放大器U1的饱和电压;当IGBT处于导通状态时,运算放大器U1从饱和状态恢复到VCE(ON)值,即Vout=VCE(ON)
反馈电阻R5用于限制运算放大器U1的放大倍数,稳定工作点
ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;
并联于分压电阻器的电容器用于改善分压器的频率响应;
下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,并触发单稳态多谐振荡器以启用采样电路,延迟时间能通过调节多振子元器件的值来调整,以为VCE提供足够的稳定时间;
运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;
电容器用于隔离所述运算放大U2的高频信号;
模数转换器用于在测量端测量通态压降VCE(ON)
本发明实施例的IGBT模块结温在线估计方法包括以下步骤:
步骤1、采用IGBT模块结温在线估计电路系统测量IGBT功率模块的 IGBT通态压降VCE(ON),并建立IGBT结温Tj查找表;
步骤2、计算IGBT功率模块中IGBT和二极管的功率损耗;
步骤3、建立IGBT等效热网络模型,并采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型;
步骤4、根据IGBT等效热网络模型以步骤1中结温Tj和步骤2中IGBT、二极管功率损耗为状态变量建立状态空间模型,并按步骤3更新IGBT等效热网络模型,进而更新状态空间模型;
步骤5、根据步骤4中的状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure BDA0002306216490000081
步骤1中,以功率变换器为例,通态压降VCE(ON)的在线测量和结温Tj查找表的建立的具体方法如下:
步骤11、在门极电压驱动下,IGBT工作在饱和区,将测量电路连接在IGBT 集电极(C)与发射极(E)之间测量IGBT通态压降VCE(ON)。如图2所示,通态压降VCE(ON)测量电路包括分压电阻器R1-R2、分压电阻器R3-R4、运算放大器U1、ESD二极管、下降沿检测器、单稳态多谐振荡器、电容器、运算放大器U2、模数转换器(ADC)。
分压电阻器R1-R2和分压电阻器R3-R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;运算放大器U1是具有差分结构的精密放大器,当IGBT器件处于关断状态时,通态压降VCE(ON)不需要测量, U1输出饱和,Vout=VH_U1,其中Vout为运放U1的输出,VH_U1为运放的饱和电压;运算放大器U1具有较短的恢复时间,当IGBT处于导通状态时,U1从饱和状态恢复到VCE(ON)值,Vout=VCE(ON);ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;并联于分压电阻器的电容器用于改善分压器的频率响应;所述下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,并触发单稳态多谐振荡器以启用采样电路,延迟时间可以通过多振子元件选择适当的值来调整,以为VCE提供足够的稳定时间;运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;电容器用于隔离所述运算放大U2的高频信号;模数转换器用于再测量端测量通态压降VCE(ON)。此通态压降VCE(ON)测量电路可实现在大电流条件下无干扰、无破坏性地在线测量通态压降VCE(ON)参数。
步骤12、IGBT结温Tj的查找表(LUT)主要通过所测通态压降VCE(ON)和负载电流IL建立。改变功率变换器负载电流IL,通过通态压降VCE(ON)测量电路,以及红外测温仪测量IGBT结温Tj,建立结温Tj-通态压降VCE(ON)-负载电流IL三维数据库,并选取一定数量的点,使用Matlab工具箱进行拟合,得到结温Tj-通态压降VCE(ON)-负载电流IL函数关系Tj=(VCE(ON),IL)。
步骤2中,计算IGBT功率模块IGBT和二极管的功率损耗的具体方法如下:
PIGBT=Pcond_T+PSW_T PDiodePcon+d_PD
Pcond_T=VCE(ON)×IL×δk Pcond_D=VCE(ON)×IL×(1-δk)
PSW=(Eon+Eoff)×fSW
其中:PIGBT为IGBT功率损耗;Pcond_T为IGBT通态损耗;PSW_T为IGBT 开关损耗;PDiode为二极管功率损耗;Pcond_D为二极管通态损耗;PSW_D为二极管开关损耗;δk为占空比;fSW为开关频率;Eon为开通过程能量损耗;Eoff为关断过程能量损耗。
步骤3中,计及焊料层疲劳累积效应的热网络更新策略的具体方法如下:
热阻抗Zθja(t)由结温Tj(t)与恒定环境温度Ta的差与阶跃输入功率PD之比来描述,即Zθja(t)=(Tj(t)-Ta)/PD。等效RC网络可以用一系列RC参数来描述Zθja(t)函数,形成如图3所示的Foster网络。Foster网络的响应用一系列指数形式表示为
Figure BDA0002306216490000101
对其作拉普拉斯(Laplace)变换,得到热阻抗在频域内传递函数的部分分式展开式为
Figure BDA0002306216490000102
其中ki和pi分别为传递函数的留数和极点,s 为复变量。通过代数运算,可以发现留数和极点与RC分量的关系公式为
Figure BDA0002306216490000103
建立4阶RC热网络行为模型,等效热阻R1与等效热容C1并联连接记为第1阶RC网络,等效热阻R2与等效热容C2并联连接记为第2阶RC网络,等效热阻R3与等效热容C3并联连接记为第3阶RC网络,等效热阻R4与等效热容C4并联连接记为第4阶RC网络。采用阶段性适时更新热网络方式,选取热阻增大10%作为热网络更新系数,以热阻增大50%作为器件完全失效判据,将模块老化过程离散为Rth、1.1×Rth、1.2×Rth、1.3×Rth、1.4×Rth五个阶段,其中 Rth为IGBT在健康状态下的热阻,每个阶段的疲劳累计为ΔD=0.2。
步骤4中,状态空间模型的建立的具体方法如下:
步骤41、自热建模:根据步骤4中热网络行为模型,选择IGBT和二极管功率损耗以及结温作为状态变量,列写状态空间方程:
状态方程:
Figure BDA0002306216490000104
输出方程:Tj(t)=Cx(t)+Du(t)
Figure BDA0002306216490000105
C=[1 1 1…1];D=[0 1];u(t)=[PD(t),Ta]
其中,An×n为系统矩阵,Bn×2为输入矩阵,Cl×n为输出矩阵,D1×2为前馈矩阵。状态向量
Figure BDA0002306216490000111
表示RC单元间的温差,u(t)为系统输入向量,其中PD(t)为功率损耗,Ta为环境温度,Ri为第i阶热网络等效热阻,Ci为第i阶热网络等效热容,本例中选用4阶热网络模型,n取4。
步骤42、考虑二极管交叉耦合效应:与步骤41相同,建立m阶二极管状态空间方程,将自热阻抗的n阶状态空间模型与交叉耦合阻抗的m阶模型相结合,构造出同时表示IGBT和二极管两种阻抗的扩展模型,得到完整的IGBT 结温状态空间模型:
Figure BDA0002306216490000112
Figure BDA0002306216490000113
其中,xs1,…,xsn为自热阻抗状态,ps1,…,psn和ks1,…,ksn分别为自热阻抗的极点和留数。类似地,xc1,…,xcm为交叉耦合阻抗的状态,pc1,…,pcm和kc1,…,kcm分别为自热阻抗的极点和留数。Tj是IGBT结温。Ta为环境温度,PIGBT和PD分别为IGBT和二极管的功率损耗。
步骤43、模型参数辨识:利用MATLAB或者Saber仿真软件中拟合工具箱拟合器件的瞬态热阻抗曲线求解参数值得健康状态下的IGBT和二极管的模型参数Ri、Ci。根据步骤4更新步骤43中状态空间模型的参数,只需对应地将自热阻抗的极点ps1,…,psn和交叉耦合阻抗的极点pc1,…,pcm更新。对应于五个老化阶段Rth、1.1×Rth、1.2×Rth、1.3×Rth、1.4×Rth分别更新为ps1,…,psn
Figure BDA0002306216490000121
Figure BDA0002306216490000122
步骤5中,计算IGBT结温Tj估计值
Figure BDA0002306216490000123
的具体方法如下:
根据步骤(4)中的状态空间模型编写自适应卡尔曼滤波程序,卡尔曼滤波分为预测和校正。在预测步骤中,模型是利用当前所算功率损耗计算预测状态
Figure BDA0002306216490000124
修正状态
Figure BDA0002306216490000125
是利用上一步功率损耗。预测误差协方差Pk -计算同理。当测量值Tk已知时,校正步骤是可行的,剩余量ek是通过测量值Tk和估计值
Figure BDA0002306216490000126
所决定。修正状态
Figure BDA0002306216490000127
的计算通过更新预测状态,剩余量ek的计算利用卡尔曼滤波增益矩阵Kk
其递推公式为:
状态估计时间更新:
Figure BDA0002306216490000128
误差协方差时间更新:
Figure BDA0002306216490000129
卡尔曼增益矩阵:
Figure BDA00023062164900001210
状态测量更新:
Figure BDA00023062164900001211
误差协方差测量更新:
Figure BDA00023062164900001212
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (10)

1.一种IGBT模块结温在线估计电路系统,其特征在于,用于在线测量IGBT功率模块的结温数据,该电路系统包括:功率变换器、通态压降VCE(ON)在线测量电路、功率损耗计算模块、结温查找表模块、计及焊料层疲劳累积效应热网络模块、自适应卡尔曼滤波模块;其中:
功率变换器,其核心器件为IGBT功率模块,用于获取负载电流IL
通态压降VCE(ON)在线测量电路,用于在线测量IGBT通态压降VCE(ON)
功率损耗计算模块,用于根据负载电流IL和IGBT通态压降VCE(ON),计算IGBT功率模块中IGBT和二极管的功率损耗;
结温查找表模块,用于建立结温Tj的查找表,通过测量获取不同负载电流IL下的通态压降VCE(ON)和结温Tj,并选取一定数量的点进行拟合,得到结温Tj与负载电流IL、通态压降VCE(ON)之间的关系;
计及焊料层疲劳累积效应热网络模块,用于建立IGBT等效热网络模型,并根据IGBT等效热网络模型以结温Tj和IGBT、二极管功率损耗为状态变量,建立状态空间模型;采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型,进而实现对状态空间模型的更新;
自适应卡尔曼滤波模块,用于根据状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure FDA0002306216480000011
2.根据权利要求1所述的IGBT模块结温在线估计电路系统,其特征在于,通态压降VCE(ON)在线测量电路包括分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4、运算放大器U1、反馈电阻R5、ESD二极管、下降沿检测器、单稳态多谐振荡器、电容器、运算放大器U2、模数转换器ADC;其中:
分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;
运算放大器U1是差分放大器,当IGBT器件处于关断状态时,通态压降VCE(ON)不需要测量,运算放大器U1输出饱和,Vout=VH_U1,其中Vout为运放U1的输出,VH_U1为运算放大器U1的饱和电压;当IGBT处于导通状态时,运算放大器U1从饱和状态恢复到VCE(ON)值,即Vout=VCE(ON)
反馈电阻R5用于限制运算放大器U1的放大倍数,稳定工作点;
ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;
并联于分压电阻器的电容器用于改善分压器的频率响应;
下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,VCE为IGBT集电极C与发射极E之间的电压,并触发单稳态多谐振荡器以启用采样电路,延迟时间能通过调节多振子元器件的值来调整,以为VCE提供足够的稳定时间;
运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;
电容器用于隔离所述运算放大U2的高频信号;
模数转换器用于在测量端测量通态压降VCE(ON),通态压降VCE(ON)为IGBT处于导通状态且为额定电流时,集电极C与发射极E之间的电压。
3.一种IGBT模块结温在线估计方法,采用权利要求1所述的IGBT模块结温在线估计电路系统,其特征在于,该方法包括以下步骤:
步骤1、采用IGBT模块结温在线估计电路系统测量IGBT功率模块的IGBT通态压降VCE(ON),并建立IGBT结温Tj查找表;
步骤2、计算IGBT功率模块中IGBT和二极管的功率损耗;
步骤3、建立IGBT等效热网络模型,并采用计及焊料层疲劳累积效应热网络更新策略,阶段性适时更新IGBT等效热网络模型;
步骤4、根据IGBT等效热网络模型以步骤1中结温Tj和步骤2中IGBT、二极管功率损耗为状态变量建立状态空间模型,并按步骤3更新IGBT等效热网络模型,进而更新状态空间模型;
步骤5、根据步骤4中的状态空间模型设计自适应卡尔曼滤波算法程序获得IGBT结温Tj估计值
Figure FDA0002306216480000021
4.根据权利要求1所述的IGBT模块结温在线估计方法,其特征在于,步骤1中在线测量IGBT通态压降VCE(ON)的具体方法如下:
分压电阻器R1、分压电阻器R2、分压电阻器R3、分压电阻器R4起分压的作用,将通态压降VCE(ON)降至运算放大器U1所允许的输入电压的范围;
运算放大器U1是差分放大器,当IGBT器件处于关断状态时,通态压降VCE(ON)不需要测量,运算放大器U1输出饱和,Vout=VH_U1,其中Vout为运放U1的输出,VH_U1为运算放大器U1的饱和电压;当IGBT处于导通状态时,运算放大器U1从饱和状态恢复到VCE(ON)值,即Vout=VCE(ON);
ESD二极管用于保护运算放大器U1的输入,抵抗高压尖峰;
并联于分压电阻器的电容器用于改善分压器的频率响应;
下降沿探测器用于比较运算放大器U1的输出和一个预定义的电压,以检测VCE的下降沿,并触发单稳态多谐振荡器以启用采样电路,延迟时间能通过调节多振子元器件的值来调整,以为VCE提供足够的稳定时间;
运算放大器U2为隔离放大器,用于实现电容隔离,并使用高频信号调制;
电容器用于隔离所述运算放大U2的高频信号;
模数转换器用于在测量端测量通态压降VCE(ON)
5.根据权利要求1所述的IGBT模块结温在线估计方法,其特征在于,步骤1中建立IGBT结温Tj查找表的具体方法如下:
通过所测通态压降VCE(ON)和负载电流IL建立IGBT结温Tj的查找表,通过测量获取不同负载电流IL下的通态压降VCE(ON)和结温Tj,并选取一定数量的点应用Matlab软件进行拟合,得到Tj=(VCE(ON),IL)。
6.根据权利要求1所述的IGBT模块结温在线估计方法,其特征在于,步骤2中计算IGBT功率模块IGBT和二极管的功率损耗的方法为:
根据通态压降、负载电流、占空比、IGBT数据手册开关损耗曲线计算IGBT和二极管功率损耗。
7.根据权利要求1所述的IGBT模块结温在线估计方法,其特征在于,步骤3中计及焊料层疲劳累积效应的热网络更新策略为:
建立四阶IGBT热网络模型,采用阶段性适时更新热网络方式,选取热阻增大10%作为热网络更新系数,以热阻增大50%作为器件完全失效判据,将模块老化过程离散为Rth、1.1×Rth、1.2×Rth、1.3×Rth、1.4×Rth五个阶段。
8.根据权利要求7所述的IGBT模块结温在线估计方法,其特征在于,步骤4中状态空间模型的建立方法如下:
步骤41、自热建模:根据IGBT热网络模型,选择选择IGBT和二极管功率损耗以及结温作为状态变量,列写状态空间方程:
状态方程:
Figure FDA0002306216480000044
输出方程:Tj(t)=Cx(t)+Du(t)
Figure FDA0002306216480000041
C=[1 1 1 … 1];D=[0 1];u(t)=[PD(t),Ta]
其中,An×n为系统矩阵,Bn×2为输入矩阵,Cl×n为输出矩阵,D1×2为前馈矩阵,状态向量
Figure FDA0002306216480000045
表示RC单元间的温差,u(t)为系统输入向量,其中PD(t)为功率损耗,Ta为环境温度;
步骤42、考虑二极管交叉耦合效应:将自热阻抗的n阶状态空间模型与交叉耦合阻抗的m阶模型相结合,构造出同时表示两种阻抗的扩展模型,得到完整的IGBT结温状态空间模型:
Figure FDA0002306216480000042
Figure FDA0002306216480000043
其中,xs1,…,xsn为自热阻抗状态,ps1,…,psn和ks1,…,ksn分别为自热阻抗的极点和留数,类似地,xc1,…,xcm为交叉耦合阻抗的状态,pc1,…,pcm和kc1,…,kcm分别为自热阻抗的极点和留数,Tj是IGBT结温,Ta为环境温度,PIGBT和PD分别为IGBT和二极管的功率损耗;
步骤43、模型参数辨识:利用仿真软件中拟合工具箱拟合器件的瞬态热阻抗曲线求解参数值得健康状态下的模型参数,再根据步骤3更新步骤42状态空间模型中对应参数。
9.根据权利要求8所述的IGBT模块结温在线估计方法,其特征在于,步骤43中采用的仿真软件为MATLAB或者Saber。
10.根据权利要求1所述的IGBT模块结温在线估计方法,其特征在于,步骤5中建立自适应卡尔曼滤波器系统,利用步骤1中IGBT结温Tj、步骤3中IGBT和二极管功率损耗为滤波器输入,计算得到结温Tj估计值
Figure FDA0002306216480000051
CN201911240970.8A 2019-12-06 2019-12-06 一种igbt模块结温在线估计电路系统及方法 Pending CN110988641A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911240970.8A CN110988641A (zh) 2019-12-06 2019-12-06 一种igbt模块结温在线估计电路系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911240970.8A CN110988641A (zh) 2019-12-06 2019-12-06 一种igbt模块结温在线估计电路系统及方法

Publications (1)

Publication Number Publication Date
CN110988641A true CN110988641A (zh) 2020-04-10

Family

ID=70090726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911240970.8A Pending CN110988641A (zh) 2019-12-06 2019-12-06 一种igbt模块结温在线估计电路系统及方法

Country Status (1)

Country Link
CN (1) CN110988641A (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111398766A (zh) * 2020-04-13 2020-07-10 武汉大学 一种igbt模块健康状态在线监测方法及系统
CN111505475A (zh) * 2020-04-28 2020-08-07 合肥阳光电动力科技有限公司 一种功率半导体模块电热模型参数的标定方法及装置
CN111783287A (zh) * 2020-06-19 2020-10-16 合肥工业大学 三相igbt功率模块的在线结温计算方法
CN111812479A (zh) * 2020-07-17 2020-10-23 中国科学院电工研究所 一种igbt通态压降在线测量电路及结温检测方法
CN111856233A (zh) * 2020-07-03 2020-10-30 武汉大学 基于键合引线退化的igbt模块可靠性评估方法及装置
CN111931402A (zh) * 2020-07-26 2020-11-13 上海电力大学 一种焊料老化状态下igbt模块结温估计方法
CN112014708A (zh) * 2020-07-27 2020-12-01 西安中车永电电气有限公司 基于FPGA的SiC功率器件在线结温计算方法
CN112364591A (zh) * 2020-09-16 2021-02-12 南京航空航天大学 SiC MOSFET结温实时动态预测方法及电压检测电路
CN113447787A (zh) * 2021-05-11 2021-09-28 上海大学 一种功率半导体器件老化在线诊断方法
CN113514747A (zh) * 2021-04-15 2021-10-19 华电(烟台)功率半导体技术研究院有限公司 一种测量电力电子器件温度分布的电学方法
CN114006353A (zh) * 2021-11-05 2022-02-01 珠海格力节能环保制冷技术研究中心有限公司 Igbt结温识别方法及空调机组
CN115436771A (zh) * 2022-09-20 2022-12-06 瑶芯微电子科技(上海)有限公司 一种结温测量电路及系统
CN116148624A (zh) * 2023-04-04 2023-05-23 无锡钊源电力电子有限公司 一种在线igbt结温测量方法
CN116484777A (zh) * 2023-06-19 2023-07-25 湖南大学 热网络模型建模及芯片温度计算方法、设备及存储介质

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636291A (zh) * 2011-02-15 2012-08-15 三一电气有限责任公司 一种igbt结温检测装置及其方法
CN103809098A (zh) * 2014-01-26 2014-05-21 中国科学院微电子研究所 一种检测igbt功率器件可靠性的系统和方法
CN104155587A (zh) * 2014-07-18 2014-11-19 浙江大学 一种igbt模块工作结温的在线检测系统及检测方法
CN107219016A (zh) * 2017-05-24 2017-09-29 湖南大学 计算igbt模块瞬态结温的方法和系统
CN109738773A (zh) * 2018-06-19 2019-05-10 北京航空航天大学 一种非平稳工况下igbt模块寿命预测方法
CN109871591A (zh) * 2019-01-24 2019-06-11 武汉大学 一种igbt功率模块在线估算结温的方法
CN110082660A (zh) * 2019-04-16 2019-08-02 合肥工业大学 基于卡尔曼滤波器的变流器igbt模块结温估计方法
CN110108999A (zh) * 2019-05-07 2019-08-09 阳光电源股份有限公司 一种igbt模块工作结温在线检测系统及方法
CN110133465A (zh) * 2019-05-16 2019-08-16 上海金脉电子科技有限公司 Igbt模块结温的计算方法及系统

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636291A (zh) * 2011-02-15 2012-08-15 三一电气有限责任公司 一种igbt结温检测装置及其方法
CN103809098A (zh) * 2014-01-26 2014-05-21 中国科学院微电子研究所 一种检测igbt功率器件可靠性的系统和方法
CN104155587A (zh) * 2014-07-18 2014-11-19 浙江大学 一种igbt模块工作结温的在线检测系统及检测方法
CN107219016A (zh) * 2017-05-24 2017-09-29 湖南大学 计算igbt模块瞬态结温的方法和系统
CN109738773A (zh) * 2018-06-19 2019-05-10 北京航空航天大学 一种非平稳工况下igbt模块寿命预测方法
CN109871591A (zh) * 2019-01-24 2019-06-11 武汉大学 一种igbt功率模块在线估算结温的方法
CN110082660A (zh) * 2019-04-16 2019-08-02 合肥工业大学 基于卡尔曼滤波器的变流器igbt模块结温估计方法
CN110108999A (zh) * 2019-05-07 2019-08-09 阳光电源股份有限公司 一种igbt模块工作结温在线检测系统及方法
CN110133465A (zh) * 2019-05-16 2019-08-16 上海金脉电子科技有限公司 Igbt模块结温的计算方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIANBO ZHOU等: "Electro-Thermal-Mechanical Multiphysics Coupling Failure Analysis Based on Improved IGBT Dynamic Model", 《IEEE ACCESS》 *
安宁: "功率器件参数自适应热网络模型的研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111398766A (zh) * 2020-04-13 2020-07-10 武汉大学 一种igbt模块健康状态在线监测方法及系统
CN111398766B (zh) * 2020-04-13 2021-05-28 武汉大学 一种igbt模块健康状态在线监测方法及系统
CN111505475A (zh) * 2020-04-28 2020-08-07 合肥阳光电动力科技有限公司 一种功率半导体模块电热模型参数的标定方法及装置
CN111783287A (zh) * 2020-06-19 2020-10-16 合肥工业大学 三相igbt功率模块的在线结温计算方法
CN111783287B (zh) * 2020-06-19 2022-10-04 合肥工业大学 三相igbt功率模块的在线结温计算方法
CN111856233B (zh) * 2020-07-03 2021-04-27 武汉大学 基于键合引线退化的igbt模块可靠性评估方法及装置
CN111856233A (zh) * 2020-07-03 2020-10-30 武汉大学 基于键合引线退化的igbt模块可靠性评估方法及装置
CN111812479A (zh) * 2020-07-17 2020-10-23 中国科学院电工研究所 一种igbt通态压降在线测量电路及结温检测方法
CN111812479B (zh) * 2020-07-17 2023-03-28 中国科学院电工研究所 一种igbt通态压降在线测量电路及结温检测方法
CN111931402A (zh) * 2020-07-26 2020-11-13 上海电力大学 一种焊料老化状态下igbt模块结温估计方法
CN111931402B (zh) * 2020-07-26 2022-05-17 上海电力大学 一种焊料老化状态下igbt模块结温估计方法
CN112014708A (zh) * 2020-07-27 2020-12-01 西安中车永电电气有限公司 基于FPGA的SiC功率器件在线结温计算方法
CN112014708B (zh) * 2020-07-27 2023-02-07 西安中车永电电气有限公司 基于FPGA的SiC功率器件在线结温计算方法
CN112364591A (zh) * 2020-09-16 2021-02-12 南京航空航天大学 SiC MOSFET结温实时动态预测方法及电压检测电路
CN113514747A (zh) * 2021-04-15 2021-10-19 华电(烟台)功率半导体技术研究院有限公司 一种测量电力电子器件温度分布的电学方法
CN113447787A (zh) * 2021-05-11 2021-09-28 上海大学 一种功率半导体器件老化在线诊断方法
CN114006353A (zh) * 2021-11-05 2022-02-01 珠海格力节能环保制冷技术研究中心有限公司 Igbt结温识别方法及空调机组
CN114006353B (zh) * 2021-11-05 2024-01-30 珠海格力节能环保制冷技术研究中心有限公司 Igbt结温识别方法及空调机组
CN115436771A (zh) * 2022-09-20 2022-12-06 瑶芯微电子科技(上海)有限公司 一种结温测量电路及系统
CN116148624A (zh) * 2023-04-04 2023-05-23 无锡钊源电力电子有限公司 一种在线igbt结温测量方法
CN116484777A (zh) * 2023-06-19 2023-07-25 湖南大学 热网络模型建模及芯片温度计算方法、设备及存储介质
CN116484777B (zh) * 2023-06-19 2023-10-13 湖南大学 热网络模型建模及芯片温度计算方法、设备及存储介质

Similar Documents

Publication Publication Date Title
CN110988641A (zh) 一种igbt模块结温在线估计电路系统及方法
US10705133B2 (en) Method and device for estimating level of damage or lifetime expectation of power semiconductor module
CN109871591A (zh) 一种igbt功率模块在线估算结温的方法
JP6742154B2 (ja) 電解コンデンサ状態をオンラインで監視する方法及びシステム
US9151725B2 (en) Humidity measurement device and method
US7400482B2 (en) Circuit breaker and method for sensing current indirectly from bimetal voltage and determining bimetal temperature and corrected temperature dependent bimetal resistance
CN111505475B (zh) 一种功率半导体模块电热模型参数的标定方法及装置
CN112989569B (zh) 一种带温度约束条件的锂电池传感器故障滤波诊断方法
CN107562978B (zh) 热阻抗拓扑结构以及热功率滤波器
Han et al. Junction temperature estimation of SiC MOSFETs based on extended Kalman filtering
Codecasa Canonical forms of one-port passive distributed thermal networks
JP2003185719A (ja) 二次電池の充電率推定装置
CN213274628U (zh) 一种提高ntc温度采样精度的电路
Lu et al. Electrothermal Modeling Based Digital Twin Method for Degradation Parameters Identification of DC-DC Converter
US11336225B2 (en) Method for determining a corrected current-voltage characteristic curve of an electrical system
CN116522773A (zh) 一种基于数据驱动的wpt系统建模方法
CN113608011B (zh) 一种具有自校准功能的直流高电压测量装置及其操作方法
CN116754915A (zh) 一种半导体开关器件的工作结温监测系统和方法
CN116380266A (zh) 一种温度检测电路和温度控制系统
KR20170119507A (ko) 인버터 스위칭 소자의 온도추정을 위한 파라미터 결정장치
CN106482859B (zh) 用于获取控制器的开关装置的多个开关单元的温度的方法和装置以及控制器的开关装置
CN112685908B (zh) 一种基于复频域的igbt结温估算快速迭代方法
CN110083992B (zh) 一种基于多新息递推最小二乘的Boost变换器故障诊断方法
CN117787055A (zh) 一种电力半导体器件寿命周期安全运行域标定方法及系统
CN116865675A (zh) 温度漂移补偿方法、装置及相关设备

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200410