CN111400888A - 一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 - Google Patents
一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 Download PDFInfo
- Publication number
- CN111400888A CN111400888A CN202010166126.1A CN202010166126A CN111400888A CN 111400888 A CN111400888 A CN 111400888A CN 202010166126 A CN202010166126 A CN 202010166126A CN 111400888 A CN111400888 A CN 111400888A
- Authority
- CN
- China
- Prior art keywords
- thermal
- heat dissipation
- igbt module
- heat
- time constant
- 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
Links
- 230000017525 heat dissipation Effects 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000001816 cooling Methods 0.000 claims abstract description 6
- 239000004519 grease Substances 0.000 claims description 5
- 229920001296 polysiloxane Polymers 0.000 claims description 4
- 239000000758 substrate Substances 0.000 claims description 3
- 230000032683 aging Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，包括如下步骤：(1)根据IGBT模块和散热系统的物理结构，建立三阶考尔型热网络模型；(2)通过改变三次散热工况，建立四组散热条件下热网络时间常数和热网络参数之间的约束方程组；(3)测量四组散热条件下，IGBT模块在降温过程中的壳温以及环境温度；(4)通过拟合壳温相对于环境温度的曲线，获得热网络时间常数；(5)最后将拟合得到的热网络时间常数代入建立的约束方程组，辨识出IGBT模块和散热系统的热阻和热容值。本发明无需测量IGBT模块的结温和功率损耗，就可以同时辨识出IGBT模块和散热系统的热参数，降低了硬件的成本以及操作的复杂性。
Description
技术领域
本发明涉及热参数辨识技术领域，尤其是一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法。
背景技术
IGBT模块及其散热系统是功率变流器的核心组件，然而工业界的调查结果表明，IGBT模块及其散热系统是功率变流器中最容易发生失效的两个组件。由于功率变流器处理功率的波动性和间歇性，功率变流器内部容易形成温度循环，加之功率模块和散热系统内部物理材料热膨胀系数的不匹配，在运行过程中容易产生热疲劳老化，降低其可靠性。相关的研究表明，随着IGBT模块及其散热系统老化程度的加剧，其对应的热参数会发生改变，因此辨识出热参数能够反映IGBT模块及其散热系统的健康状态。
目前来说，热参数的辨识方法主要有以下三类：解析计算法、有限元仿真法以及实验测量法。解析计算法是根据IGBT模块及其散热系统的机械尺寸以及材料物理参数，计算出热阻和热容值。有限元仿真法则是通过在软件中搭建IGBT模块及其散热系统的物理模型，利用仿真得到的温度响应来计算出热参数值。尽管这两类方法不需要硬件电路，成本较低，但是由于无法计及老化进程中材料的热疲劳损伤，因而辨识出的热参数只能反映IGBT模块以及散热系统的初始健康状态，不适用于状态监测应用。
实验测量法是指利用实验来辨识出IGBT模块以及散热系统的热阻和热容值。常见的实验测量法主要有IEC-60747-9-2007测量标准和JESD51-14测量标准。前者通过测量热平衡状态下的结温，壳温，环境温度以及功率损耗，计算出IGBT模块和散热系统的热阻。后者通过测量加热功率损耗以及降温阶段的结温降温曲线来辨识出IGBT模块以及散热系统的热参数。尽管上述测量标准可以较为准确地辨识出热参数，但是由于需要测量结温以及功率损耗，并且要满足热平衡条件，所以仅适用于离线的情形。少量文献根据易于提取的其它物理特征量来间接辨识出热参数值，这些特征量包括变流器输出电流五次谐波含量、壳温的二维分布、IGBT模块壳到环境的温差等等。虽然这些物理特征量易于在线测量，且能够给出IGBT模块总热阻的数值，但是这些特征量容易受散热系统老化的影响，并不能满足实际应用需求。
因此，需要发明一种无需测量结温和功率损耗就可以同时辨识出IGBT模块和散热系统热参数的方法。
发明内容
本发明所要解决的技术问题在于，提供一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，
为解决上述技术问题，本发明提供一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，包括如下步骤：
(1)根据IGBT模块和散热系统的物理结构，建立三阶考尔型热网络模型；
(2)通过改变三次散热工况，建立四组散热条件下热网络时间常数和热网络参数之间的约束方程组；
(3)测量四组散热条件下，IGBT模块在降温过程中的壳温以及环境温度；
(4)通过拟合壳温相对于环境温度的曲线，获得热网络时间常数；
(5)最后将拟合得到的热网络时间常数代入建立的约束方程组，辨识出IGBT模块和散热系统的热阻和热容值。
优选的，步骤(1)中，将IGBT模块和散热系统等效为三阶考尔型热网络模型，并将IGBT模块设置为热网络模型的第一阶RC参数，导热硅脂设置为热网络模型的第二阶RC参数，散热器设置为热网络模型的第三阶RC参数。
优选的，步骤(2)中，四组散热条件下热网络时间常数和热网络参数之间的约束方程组为：
其中R1，C1为IGBT模块的热阻和热容值，R2，C2为导热硅脂的热阻和热容值，R3，C3为散热条件FS1下散热器的热阻和热容值，R′3，C′3为散热条件FS2下散热器的热阻和热容值，R″3，C″3为散热条件FS3下散热器的热阻和热容值，R″′3，C″′3为散热条件FS4下散热器的热阻和热容值。τ1，τ2，τ3、τ′1，τ′2，τ′3、τ″1，τ″2，τ″3和τ″′1，τ″′2，τ″′3为四组散热条件下热网络时间常数。
优选的，步骤(3)中，利用热电偶测量壳温Tc(t)以及环境温度Ta，壳温测量位置在芯片中心正下方IGBT模块基板上，环境温度则为散热器进风口的空气温度。
优选的，步骤(4)中，利用非线性最小二乘法，通过拟合壳温相对于环境温度Tca(t)＝Tc(t)-Ta的曲线来获得热网络时间常数，拟合形式为
其中A1A2A3为拟合系数，t为测量壳温过程中的时间。
优选的，步骤(5)中，将四组散热条件下拟合得到的热网络时间常数带入约束方程组中，求解出IGBT模块和散热系统的热阻和热容值。
本发明的有益效果为：本发明无需测量IGBT模块的结温和功率损耗，就可以同时辨识出IGBT模块和散热系统的热参数，降低了硬件的成本以及操作的复杂性。
附图说明
图1为本发明的方法流程示意图。
图2(a)为本发明中IGBT模块和散热系统的三阶考尔型热网络散热工况FS1时的模型示意图。
图2(b)为本发明中IGBT模块和散热系统的三阶考尔型热网络散热工况FS2时的模型示意图。
图2(c)为本发明中IGBT模块和散热系统的三阶考尔型热网络散热工况FS3时的模型示意图。
图2(d)为本发明中IGBT模块和散热系统的三阶考尔型热网络散热工况FS4时的模型示意图。
具体实施方式
如图1所示，一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，包括如下步骤：
(1)根据IGBT模块和散热系统的物理结构，建立三阶考尔型热网络模型；
本实施例中，步骤(1)具体包括：根据IGBT模块和散热系统的物理结构，将IGBT模块和散热系统等效为三阶考尔型热网络模型，如图2(a)、图2(b)、图2(c)和图2(d)所示，其中第一阶参数R1,C1和IGBT模块相关，第二阶参数R2,C2和导热硅脂相关，第三阶参数和散热器相关，R3,C3为散热条件FS1下散热器的热参数，R′3,C′3为散热条件FS2下散热器的热参数，R″3,C″3为散热条件FS3下散热器的热参数，R″′3,C″′3为散热条件FS4下散热器的热参数。
(2)利用电热比拟理论以及特征方程根与系数的关系，建立热网络时间常数和热网络参数之间的约束方程组；
在本实施例中，步骤(2)具体包括：建立的热网络自然频率和热网络参数之间的约束方程组，如下所示：
其中τ1，τ2，τ3、τ′1，τ′2，τ′3、τ″1，τ″2，τ″3和τ″′1，τ″′2，τ″′3为四组散热条件下热网络时间常数。
(3)测量四组散热条件下，IGBT模块在降温过程中的壳温以及环境温度；
在本实施例中，步骤(3)具体包括：利用热电偶测量壳温Tc(t)以及环境温度Ta，壳温测量位置在芯片中心正下方IGBT模块基板上，环境温度则为散热器进风口的空气温度。
(4)借助拟合壳温相对于环境温度的曲线，获得热网络的自然频率；
在本实施例中，步骤(4)具体包括：利用非线性最小二乘法，通过拟合壳温相对于环境温度Tca(t)＝Tc(t)-Ta的曲线来获得热网络时间常数，拟合形式为
其中A1A2A3为拟合系数，t为测量壳温过程中的时间。
(5)将四组散热条件下拟合得到的热网络时间常数代入约束方程组中，辨识出IGBT模块和散热系统的热阻和热容值。
在本实施例中，步骤(5)具体包括：将步骤(4)中拟合得到的τ1，τ2，τ3、τ′1，τ′2，τ′3、τ″1，τ″2，τ″3和τ″′1，τ″′2，τ″′3带入步骤(2)中建立的约束方程组中，辨识出IGBT模块和散热系统的热阻和热容值。
在本实施例中，IGBT模块及其散热系统热参数辨识方法主要包括三阶等效考尔型热网络模型的建立、热网络时间常数和热网络参数之间的约束方程组、利用壳温以及环境温度实现热网络时间常数的提取，本发明可以对IGBT模块及其散热系统热参数同时实现辨识，而且无需测量功率损耗以及结温，降低了硬件的成本以及操作的复杂性。
Claims (6)
1.一种基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，包括如下步骤：
(1)根据IGBT模块和散热系统的物理结构，建立三阶考尔型热网络模型；
(2)通过改变三次散热工况，建立四组散热条件下热网络时间常数和热网络参数之间的约束方程组；
(3)测量四组散热条件下，IGBT模块在降温过程中的壳温以及环境温度；
(4)通过拟合壳温相对于环境温度的曲线，获得热网络时间常数；
(5)最后将拟合得到的热网络时间常数代入建立的约束方程组，辨识出IGBT模块和散热系统的热阻和热容值。
2.如权利要求1所述的基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，步骤(1)中，将IGBT模块和散热系统等效为三阶考尔型热网络模型，并将IGBT模块设置为热网络模型的第一阶RC参数，导热硅脂设置为热网络模型的第二阶RC参数，散热器设置为热网络模型的第三阶RC参数。
3.如权利要求1所述的基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，步骤(2)中，四组散热条件下热网络时间常数和热网络参数之间的约束方程组为：
其中R1,C1为IGBT模块的热阻和热容值，R2,C2为导热硅脂的热阻和热容值，R3,C3为散热条件FS1下散热器的热阻和热容值，R′3,C′3为散热条件FS2下散热器的热阻和热容值，R″3,C″3为散热条件FS3下散热器的热阻和热容值，R″′3,C″′3为散热条件FS4下散热器的热阻和热容值。τ1,τ2,τ3、τ′1,τ′2,τ′3、τ″1,τ″2,τ″3和τ″′1,τ″′2,τ″′3为四组散热条件下热网络时间常数。
4.如权利要求1所述的基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，步骤(3)中，利用热电偶测量壳温Tc(t)以及环境温度Ta，壳温测量位置在芯片中心正下方IGBT模块基板上，环境温度则为散热器进风口的空气温度。
5.如权利要求1所述的基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，步骤(4)中，利用非线性最小二乘法，通过拟合壳温相对于环境温度Tca(t)＝Tc(t)-Ta的曲线来获得热网络时间常数，拟合形式为
其中A1A2A3为拟合系数，t为测量壳温过程中的时间。
6.如权利要求1所述的基于热网络时间常数的IGBT模块和散热系统的热参数辨识方法，其特征在于，步骤(5)中，将四组散热条件下拟合得到的热网络时间常数带入约束方程组中，求解出IGBT模块和散热系统的热阻和热容值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010166126.1A CN111400888A (zh) | 2020-03-11 | 2020-03-11 | 一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010166126.1A CN111400888A (zh) | 2020-03-11 | 2020-03-11 | 一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111400888A true CN111400888A (zh) | 2020-07-10 |
Family
ID=71436111
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010166126.1A Pending CN111400888A (zh) | 2020-03-11 | 2020-03-11 | 一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111400888A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114599A (zh) * | 2020-09-25 | 2020-12-22 | 中车株洲电力机车研究所有限公司 | 一种功率器件的温度控制方法及其装置 |
| CN112327124A (zh) * | 2020-09-30 | 2021-02-05 | 河海大学 | 一种监测igbt模块热疲劳老化的方法及不均匀老化的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105224819A (zh) * | 2015-11-10 | 2016-01-06 | 重庆大学 | 基于数值迭代的功率器件散热器动态响应性能优化方法 |
| CN105718694A (zh) * | 2016-01-28 | 2016-06-29 | 重庆大学 | 基于igbt结温信息的热网络参数辨识方法 |
| CN106443400A (zh) * | 2016-09-14 | 2017-02-22 | 河北工业大学 | 一种igbt模块的电‑热‑老化结温计算模型建立方法 |
-
2020
- 2020-03-11 CN CN202010166126.1A patent/CN111400888A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105224819A (zh) * | 2015-11-10 | 2016-01-06 | 重庆大学 | 基于数值迭代的功率器件散热器动态响应性能优化方法 |
| CN105718694A (zh) * | 2016-01-28 | 2016-06-29 | 重庆大学 | 基于igbt结温信息的热网络参数辨识方法 |
| CN106443400A (zh) * | 2016-09-14 | 2017-02-22 | 河北工业大学 | 一种igbt模块的电‑热‑老化结温计算模型建立方法 |
Non-Patent Citations (2)
| Title |
|---|
| JUN ZHANG 等: "Thermal Parameter Monitoring of IGBT Module Using Case Temperature", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 * |
| SHUAI ZHENG 等: "Measurement of Thermal Parameters of SiC MOSFET Module by Case Temperature", 《IEEE JOURNAL OF EMERGING AND SELECTED TOPICS IN POWER ELECTRONICS》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112114599A (zh) * | 2020-09-25 | 2020-12-22 | 中车株洲电力机车研究所有限公司 | 一种功率器件的温度控制方法及其装置 |
| CN112114599B (zh) * | 2020-09-25 | 2021-07-13 | 中车株洲电力机车研究所有限公司 | 一种功率器件的温度控制方法及其装置 |
| CN112327124A (zh) * | 2020-09-30 | 2021-02-05 | 河海大学 | 一种监测igbt模块热疲劳老化的方法及不均匀老化的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lineykin et al. | Analysis of thermoelectric coolers by a spice-compatible equivalent-circuit model | |
| CN108614005B (zh) | 一种基于多层复合材料的热物性参数的测试方法及系统 | |
| CN111400888A (zh) | 一种基于热网络时间常数的igbt模块和散热系统的热参数辨识方法 | |
| Evans et al. | Automated fast extraction of compact thermal models for power electronic modules | |
| CN106960085A (zh) | 一种在线预估igbt元件结温的方法 | |
| CN106124955B (zh) | 液冷冷板热阻的瞬态电学测试方法 | |
| CN109815596B (zh) | 基于温控散热器的半导体器件环境温度模拟系统及方法 | |
| CN110991123A (zh) | 基于l-m算法优化的变压器绕组热点瞬态温度计算方法 | |
| CN114544699B (zh) | 一种材料热阻和导热系数的测试方法 | |
| US8078438B2 (en) | Method for simulating thermal resistance value of thermal test die | |
| Han et al. | A thermal resistance network model based on three-dimensional structure | |
| Free et al. | Computational heat transfer platform for frequency domain reliability design | |
| CN104569065A (zh) | 一种大功率led器件固晶层散热性能的快速评估方法 | |
| CN217901039U (zh) | 一种测温装置及测温系统 | |
| Gao et al. | Beat the deviations in estimating maximum power of thermoelectric modules | |
| CN108846171A (zh) | 仿真mosfet温度电学特性的子电路模型的建立方法 | |
| CN111103319B (zh) | 一种导热硅脂和散热风扇的健康状态监测方法 | |
| CN114509185A (zh) | 低温环境声表面波温度传感器测量精度分区系统及方法 | |
| CN114398809A (zh) | 一种考虑散热器的igbt模块结温计算方法 | |
| Bazzo et al. | Digital control system using a thermoelectric cell for temperature electronic devices testing | |
| CN113420407A (zh) | 一种igct水冷散热器建模及其结温计算方法 | |
| CN111400884A (zh) | 一种基于热网络自然频率的碳化硅模块热参数测量方法 | |
| Sprunck et al. | Thermal Conductivity Measurement Setup for Pad and Paste Thermal Interface Materials | |
| Möller et al. | Building reliable FE simulation models for a better behavior prediction of power electronic systems | |
| Sridhar et al. | An accurate experimental apparatus for measuring losses in magnetic components |
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: 20200710 |