CN105335602A - 一种功率igbt模块的寿命预测方法 - Google Patents
一种功率igbt模块的寿命预测方法 Download PDFInfo
- Publication number
- CN105335602A CN105335602A CN201410393391.8A CN201410393391A CN105335602A CN 105335602 A CN105335602 A CN 105335602A CN 201410393391 A CN201410393391 A CN 201410393391A CN 105335602 A CN105335602 A CN 105335602A
- Authority
- CN
- China
- Prior art keywords
- life
- omega
- igbt module
- temperature
- formula
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012544 monitoring process Methods 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 3
- 238000012935 Averaging Methods 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 description 7
- 230000032683 aging Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Landscapes
- Power Conversion In General (AREA)
Abstract
一种功率IGBT模块的寿命预测方法,该方法基于实时温度监测,首先,由实时监测的温度获得其离散有限序列,将此序列进行线性拟合,得到连续函数,再经傅里叶变换为正弦量,按电路解析方法,通过寿命损耗调节函数C,求解在寿命消耗函数R(T)和Rc上寿命损耗,进而预测出功率IGBT模块的剩余寿命。
Description
技术领域
本发明涉及新能源发电系统可靠性技术领域,特别是涉及到一种核心变流装置功率IGBT模块的寿命预测。
背景技术
近年来,以功率IGBT模块为核心的变流装置在新能源发电、电动汽车领域广泛应用。然而在实际应用中,负载的不稳定性和剧烈波动性,加快了作为功率变流装置核心器件的IGBT的老化,系统可靠性降低。因此对功率IGBT模块的寿命预测显得十分重要。在影响导致功率IGBT模块寿命的因素中,结温及结温变化是导致功率IGBT模块损坏最直接、最重要的因素。但现有的基于结温监控预测IGBT寿命的方法,均需基于设定工况作出,只能预测出设定工况下IGBT的寿命,而对于电动汽车、新能源发电等负载不稳定、且变化不规律的使用场合,采用现有预测方法的预测寿命与实际寿命存在较大误差。无法对其进行准确的实时剩余寿命预测。
中国专利CN101087125B公开了一种具有寿命预测功能的电动汽车电机驱动系统,该文件中指出可以通过对IGBT结温及结温变化实时预测其剩余寿命,公式及方法如下:
LIGBT=N(Tj,ΔTj)·tcycle(5)
式中:N(Tj,ΔTj)是在IGBT结温Tj和结温变化ΔTj下循环功率曲线,tcycle是功率循环时间,单位:s,这两者是由IGBT制造商给出。
该方法中计算剩余寿命需依靠IGBT制造商给出的功率循环曲线和功率循环时间,但这两个参数的同样是基于某种设定工况作出的,当实际工况与设定工况不符时,其预测出的剩余寿命也并不准确,因此,提供一种能够适应各种工况,基于结温及结温变化准确实时预测IGBT剩余寿命的方法,成为现有技术中亟待解决的问题。
发明内容
为了克服上述现有技术的不足,提供一种能够适应各种工况,基于结温实时检测准确实时预测IGBT剩余寿命的方法。
本发明提供了一种功率IGBT模块的寿命预测方法,其特征在于该方法基于实时温度监测,
首先,由实时监测的温度获得其离散有限序列,将此序列进行线性拟合得到连续函数,再经傅里叶变换为正弦量,按电路解析方法,通过寿命损耗调节函数C,求解在寿命消耗函数R(T)和Rc上寿命损耗,进而预测出功率IGBT模块的剩余寿命。
所述的功率IGBT模块的寿命预测方法,其特征在于IGBT模块结温以一定频率和幅值波动时,功率IGBT模块的瞬时寿命损耗P依下(1.1)式计算
所述式(1)中T(t)为将实时监测得到结温T的离散有限序列进行线性拟合得到表征T随时间t变化的连续函数;ω为角频率,ω=2πfT,fT为IGBT模块结温变化频率,
所述功率IGBT模块累积消耗的寿命θs根据式(2.1)计算
当IGBT以恒定结温运行时,P依下(1.2)式计算,
累积消耗寿命根据式(2.2)计算
ts为IGBT模块实际运行时间,tz是实际工况下的最终可运行时间
模块剩余寿命Lr如式(3)tz是实际工况下的最终可运行时间。
Lr=L-θs(3),L为功率IGBT总寿命,
所述式(1)的参数中,R(T)、Rc和C的获取过程如下:
需要选取相同型号的功率IGBT模块3组,每组m个,m≥1,
将3组模块分别在按以下步骤操作。
1)模块组1运行在恒定结温T1下,得到n1个采样时间点,对应的温度序列和剩余寿命值,每组得到m个总寿命值,将总寿命值求平均值得L1。
2)让其他两组模块围绕结温T1,分别在正弦波动下运行,所述正弦波幅值固定,频率分别为f1、f2,分别得到对应的温度序列,m个寿命值,各自将总寿命值求平均值得L2,L3。
L1带入(1.2),L2,L3带入式(1.2)联立可得到R(T)、Rc和C的离散值,改变T1,得到不同温度状态下的R(T)、Rc和C参数,将这3个参数矩阵进行拟合即可得到R(T)、Rc及C。
其中R(T)、Rc的单位为T2/t,C的单位为t/T2,T为温度单位,t为时间单位,
与现有技术相比,本发明的有益效果是:
本发明提供的功率IGBT模块的寿命预测方法,基于实时结温监测数据,有效推测功率IGBT模块的剩余寿命,与现有的结温检测方法只能在结温或结温变化幅值过高时向控制系统发出警告,并不能判断模块是否真正损坏,更不能由此推测模块的剩余寿命相比,能够基于基于实时结温监测的功率IGBT模块的剩余寿命。本发明充分考虑到了结温这一影响功率IGBT模块寿命的最重要因素,将模型构建在功率IGBT模块结温的实时监测之上,综合考虑到了恒定高结温以及结温频率与幅值波动这两种影响因素。且与CN101087125B中的寿命预测方法相比,本发明提供的寿命预测方法其关键参数R(T)、Rc、C建立在对现有所检测功率IGBT模块通过简单方法实测并推导得来,参数硬件实现简单,经济成本低,因此在复杂工况下对功率IGBT剩余寿命的进行实施检测时,其预测结果准确度更高。
本发明中各个参数的获取计算简单,过程明确,参数一旦获取,即可直接由实时监测结温数据预测功率IGBT模块的寿命,具有相当的实用性。
附图说明
图1为实施例1所示功率IGBT模块寿命预测方法步骤1)操作等效电路模型图
图2为实施例1所示功率IGBT模块寿命预测方法步骤2)操作等效电路模型图。
具体实施方式
实施例1功率IGBT模块预测寿命方法推导
所述功率IGBT模块的瞬时寿命损耗P依下(1.1)式计算
所述式(1)中T(t)为将实时监测得到结温T的离散有限序列进行线性拟合得到表征T随时间t变化的连续函数;ω为角频率,ω=2πfT,fT为IGBT模块结温变化频率,
所述功率IGBT模块累积消耗的寿命θs根据式(2.1)计算
当IGBT以恒定结温运行时,P依下(1.2)式计算,
累积消耗寿命根据式(2.2)计算
ts为IGBT模块实际运行时间,tz是实际工况下的最终可运行时间
模块剩余寿命Lr如式(3)tz是实际工况下的最终可运行时间。
Lr=L-θs(3),L为功率IGBT总寿命,
所述式(1)的参数中,R(T)、Rc和C的获取过程如下:
需要选取相同型号的功率IGBT模块3组,每组m个,m≥1,
将3组模块分别在按以下步骤操作。
1)模块组1运行在恒定结温T1下,得到n1个采样时间点,对应的温度序列和剩余寿命值,每组得到m个总寿命值,将总寿命值求平均值得L1。
2)让其他两组模块围绕结温T1,分别在正弦波动下运行,所述正弦波幅值固定,频率分别为f1、f2,分别得到对应的温度序列,m个寿命值,各自将总寿命值求平均值得L2,L3。
L1带入(1.2),L2,L3带入式(1.2)联立可得到R(T)、Rc和C的离散值,改变T1,得到不同温度状态下的R(T)、Rc和C参数,将这3个参数矩阵进行拟合即可得到R(T)、Rc及C。
其中R(T)、Rc的单位为T2/t,C的单位为t/T2,T为温度单位,t为时间单位,
在本实施例中各参数单位及其描述如下
参数名称 | 符号 | 单位 | 参数描述 |
总寿命 | L | 年 | 功率IGBT模块的总寿命 |
结温 | T | 摄氏度 | 功率IGBT模块最高温度 |
耗寿命函数1 | R(T) | 摄氏度2/年 | 模拟结温对功率IGBT模块寿命的影响 |
耗寿命函数2 | Rc | 摄氏度2/年 | 模拟结温波动对功率IGBT模块寿命的影响 |
寿命损耗调节函数 | C | 年/摄氏度2 | 调节结温波动的幅值和频率 |
进行步骤1)操作时等效电路模型如图1所示,耗寿命函数1等效为阻性负载R(T)、R(T)作为寿命消耗模块与等效为待测IGBT模块的温度源T串联组成等效电路。
进行步骤2)操作时等效电路模型如图2所示,耗寿命函数1、2分别等效为阻性负载R(T)、Rc,寿命损耗调节函数等效为可调节电容C,可调节电容C与阻性负载Rc件串联后,再与R(T)并联组成寿命消耗模块,寿命消耗模块与等效为待测IGBT模块的温度源T串联组成等效电路。
在运行时
由温度源T向整个装置提供实际结温T,由R(T)和Rc作为负载不断消耗寿命,并C来调节总负载变化来从实现工况变化,并可检测由结温T波动所引起的功率IGBT模块老化速度的变化。
实施例2方法验证
以下通过设计工况变化,验证实施例1提供的功率IGBT模块预测寿命方法。
工况1采用某型号功率IGBT模块,其设计指标为运行在恒温75℃时,使用寿命20年,即T=75℃,L75=20年,
恒温运行时Rc和C为常量,因此在推导参数时相当于仅进行步骤1)得到L1=L75=20年,且θs=L1,且带入式(1)得到R(T)即可
通过计算得到:参数R(T)=5625。
该型号功率IGBT模块在恒温105℃时,可使用寿命为:
工况2采用工况1同型号功率IGBT模块按如下表格设计工况持续工作,通过步骤1)的模型计算可得每阶段的剩余寿命,剩余寿命值如表格所示:
上述剩余寿命均为是以原设计使用工况下的寿命预测值。经过加速试验验证,上述预测值与实际检测值一致。
工况3采用某型号功率IGBT模块,其设计指标为运行在恒温120℃时,使用寿命1年,按照该设计指标,采用实施例1的方法推导,R(T)、Rc及C,m=1
步骤1)按照已知设计指标带入即可,即T1=120℃,L1=θs=1年,带入式(1.2);
步骤2)另取同型号IGBT模块2个,分别以幅值不变,频率为f1和f2的正弦波动下运行,其结温变化曲线函数分别为
分别带入式(1.1)中,通过计算得到:参数R(T)=14400,Rc=-3947,C2=-9.35*10-3。
当该型IGBT模块工作时实际结温变化曲线为时,通过计算得其总寿命为:L=1.9年。
当改变T1时,L1可以按照工况1中的方法计算得到,通过得到不同温度状态下的R(T),Rc,C参数,将这3个参数矩阵进行拟合即可得到R(T)、Rc及C,并据此可以对各种复杂工况下的实时剩余寿命进行预测。
Claims (2)
1.一种功率IGBT模块的寿命预测方法,其特征在于该方法基于实时温度监测,首先,由实时监测的温度获得其离散有限序列,将此序列进行线性拟合,得到连续函数,再经傅里叶变换为正弦量,按电路解析方法,通过寿命损耗调节函数C,求解在寿命消耗函数R(T)和Rc上寿命损耗,进而预测出功率IGBT模块的剩余寿命。
2.如权利要求1所述的功率IGBT模块的寿命预测方法,其特征在于IGBT模块结温以一定频率和幅值波动时,功率IGBT模块的瞬时寿命损耗P依下(1.1)式计算
所述式(1)中T(t)为将实时监测得到结温T的离散有限序列进行线性拟合得到表征T随时间t变化的连续函数;ω为角频率,ω=2πfT,fT为IGBT模块结温变化频率,
所述功率IGBT模块累积消耗的寿命θs根据式(2.1)计算
当IGBT以恒定结温运行时,P依下(1.2)式计算,
累积消耗寿命根据式(2.2)计算
ts为IGBT模块实际运行时间,tz是实际工况下的最终可运行时间
模块剩余寿命Lr如式(3)tz是实际工况下的最终可运行时间,
Lr=L-θs(3),L为功率IGBT总寿命;
所述式(1)的参数中,R(T)、Rc和C的获取过程如下:
需要选取相同型号的功率IGBT模块3组,每组m个,m≥1,
将3组模块分别在按以下步骤操作
1)模块组1运行在恒定结温T1下,得到n1个采样时间点,对应的温度序列和剩余寿命值,每组得到m个总寿命值,将总寿命值求平均值得L1;
2)让其他两组模块围绕结温T1,分别在正弦波动下运行,所述正弦波幅值固定,频率分别为f1、f2,分别得到对应的温度序列,m个寿命值,各自将总寿命值求平均值得L2,L3;
L1带入(1.2),L2,L3带入式(1.2)联立可得到R(T)、Rc和C的离散值,改变T1,得到不同温度状态下的R(T)、Rc和C参数,将这3个参数矩阵进行拟合,即得到R(T)、Rc及C。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410393391.8A CN105335602B (zh) | 2014-08-12 | 2014-08-12 | 一种功率igbt模块的寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410393391.8A CN105335602B (zh) | 2014-08-12 | 2014-08-12 | 一种功率igbt模块的寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105335602A true CN105335602A (zh) | 2016-02-17 |
CN105335602B CN105335602B (zh) | 2017-11-17 |
Family
ID=55286123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410393391.8A Expired - Fee Related CN105335602B (zh) | 2014-08-12 | 2014-08-12 | 一种功率igbt模块的寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105335602B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169650A (zh) * | 2016-12-06 | 2018-06-15 | 深圳市蓝海华腾技术股份有限公司 | 一种检测igbt使用寿命是否达标的方法及装置 |
CN109284543A (zh) * | 2018-09-04 | 2019-01-29 | 河北工业大学 | 基于最优尺度高斯过程模型的igbt剩余寿命预测方法 |
CN109387779A (zh) * | 2018-09-28 | 2019-02-26 | 河北工业大学 | 一种基于统计数据驱动的万能式断路器操作附件剩余寿命预测方法 |
US10288696B2 (en) | 2016-11-16 | 2019-05-14 | Industrial Technology Research Institute | Intelligent diagnosis system for power module and method thereof |
CN109831102A (zh) * | 2017-11-21 | 2019-05-31 | 河北工业大学 | 一种损耗计算的方法、装置及计算机可读存储介质 |
CN110703066A (zh) * | 2019-11-12 | 2020-01-17 | 云南电网有限责任公司电力科学研究院 | 预测晶体管故障的方法及其系统 |
CN111208404A (zh) * | 2020-01-17 | 2020-05-29 | 西安工程大学 | 一种SiC MOSFET模块的寿命预测方法及系统 |
CN112098797A (zh) * | 2020-09-17 | 2020-12-18 | 南方电网科学研究院有限责任公司 | 一种SiC功率模块热阻测量方法 |
US12081208B2 (en) | 2022-09-02 | 2024-09-03 | Industrial Technology Research Institute | High power multiple frequency coupling generator and driving method thereof |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101087125A (zh) * | 2007-06-25 | 2007-12-12 | 中国科学院电工研究所 | 一种具有寿命预测功能的电动汽车电机驱动系统 |
US20120078559A1 (en) * | 2010-09-27 | 2012-03-29 | Hitachi, Ltd. | Fast steady state field analysis method, fast steady state field analysis apparatus, fast steady state field analysis program, and computer-readable recording medium of storing its program |
-
2014
- 2014-08-12 CN CN201410393391.8A patent/CN105335602B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101087125A (zh) * | 2007-06-25 | 2007-12-12 | 中国科学院电工研究所 | 一种具有寿命预测功能的电动汽车电机驱动系统 |
US20120078559A1 (en) * | 2010-09-27 | 2012-03-29 | Hitachi, Ltd. | Fast steady state field analysis method, fast steady state field analysis apparatus, fast steady state field analysis program, and computer-readable recording medium of storing its program |
Non-Patent Citations (3)
Title |
---|
方鑫等: "IGBT模块寿命预测模型综述", 《电源学报》 * |
王彦刚等: "功率IGBT模块的寿命预测", 《机车电传动》 * |
王瑞: "IGBT功率模块寿命预测的研究与设计", 《电子测试》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10288696B2 (en) | 2016-11-16 | 2019-05-14 | Industrial Technology Research Institute | Intelligent diagnosis system for power module and method thereof |
CN108169650B (zh) * | 2016-12-06 | 2020-04-14 | 深圳市蓝海华腾技术股份有限公司 | 一种检测igbt使用寿命是否达标的方法及装置 |
CN108169650A (zh) * | 2016-12-06 | 2018-06-15 | 深圳市蓝海华腾技术股份有限公司 | 一种检测igbt使用寿命是否达标的方法及装置 |
CN109831102A (zh) * | 2017-11-21 | 2019-05-31 | 河北工业大学 | 一种损耗计算的方法、装置及计算机可读存储介质 |
CN109284543A (zh) * | 2018-09-04 | 2019-01-29 | 河北工业大学 | 基于最优尺度高斯过程模型的igbt剩余寿命预测方法 |
CN109284543B (zh) * | 2018-09-04 | 2023-05-23 | 河北工业大学 | 基于最优尺度高斯过程模型的igbt剩余寿命预测方法 |
CN109387779A (zh) * | 2018-09-28 | 2019-02-26 | 河北工业大学 | 一种基于统计数据驱动的万能式断路器操作附件剩余寿命预测方法 |
CN109387779B (zh) * | 2018-09-28 | 2020-08-04 | 河北工业大学 | 一种基于统计数据驱动的万能式断路器操作附件剩余寿命预测方法 |
CN110703066A (zh) * | 2019-11-12 | 2020-01-17 | 云南电网有限责任公司电力科学研究院 | 预测晶体管故障的方法及其系统 |
CN110703066B (zh) * | 2019-11-12 | 2021-11-12 | 云南电网有限责任公司电力科学研究院 | 预测晶体管故障的方法及其系统 |
CN111208404A (zh) * | 2020-01-17 | 2020-05-29 | 西安工程大学 | 一种SiC MOSFET模块的寿命预测方法及系统 |
CN112098797A (zh) * | 2020-09-17 | 2020-12-18 | 南方电网科学研究院有限责任公司 | 一种SiC功率模块热阻测量方法 |
CN112098797B (zh) * | 2020-09-17 | 2021-12-14 | 南方电网科学研究院有限责任公司 | 一种SiC功率模块热阻测量方法 |
US12081208B2 (en) | 2022-09-02 | 2024-09-03 | Industrial Technology Research Institute | High power multiple frequency coupling generator and driving method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN105335602B (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105335602A (zh) | 一种功率igbt模块的寿命预测方法 | |
CN111758036B (zh) | 用于监测运行中的电力设备的运行状态的系统和方法 | |
EP2551689B1 (en) | Device for calculating impedances of battery cell and battery impedance measuring system | |
US9448287B2 (en) | Battery monitoring device | |
CN103399241B (zh) | 基于温升与负荷关系的配电变压器故障诊断系统及方法 | |
CN103278771B (zh) | 工业现场三相异步电机异常状态的辨识方法 | |
CN103104509A (zh) | 变频水泵全工况运行状态获取方法 | |
WO2014122705A1 (ja) | 太陽電池パネル監視プログラム、太陽電池パネル監視装置及び太陽電池パネル監視方法 | |
CN103324863B (zh) | 使用实测湍流强度修正风力发电机的测量功率曲线的方法 | |
CN111211719B (zh) | 一种永磁同步电机转子磁钢温度估算方法及系统 | |
US20150198144A1 (en) | Method of automatically calculating power curve limit for power curve monitoring of wind turbine | |
CN103592528B (zh) | 一种基于动态轨迹灵敏度的光伏逆变器模型参数辨识方法 | |
CN103235239B (zh) | 多机系统中基于无功功率积分的振荡中心定位分析方法 | |
CN104198138A (zh) | 一种风力发电机组异常振动的预警方法及系统 | |
CN104597132A (zh) | 一种基于共振声学的瓷支柱绝缘子带电检测方法 | |
CN104833935A (zh) | 一种针对有功功率方向法定位谐波源局限性的检验方法 | |
CN101566663A (zh) | 一种配电系统电压跌落源定位方法 | |
CN104074687A (zh) | 一种用于兆瓦级风电机组的载荷及性能测试方法及装置 | |
CN103353926A (zh) | 一种电机温度分布实时监测方法 | |
CN106682407A (zh) | 基于戴维南等值和支路传输功率极限的电压稳定评估方法 | |
CN105760669A (zh) | 一种风电变流器功率模块故障率评估方法及系统 | |
CN105572455A (zh) | 一种基于谐波功率监测的谐波电压责任的测定方法 | |
CN114033617B (zh) | 一种控制参量自适应调整的可控风力发电方法及系统 | |
CN104716902A (zh) | 光伏组件的最大功率测量装置、方法及发电量测量方法 | |
CN103389457B (zh) | 通过确定熔体激活能检测熔断器寿命的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171117 |