CN112699579A - IGBT state monitoring method based on finite element simulation - Google Patents
IGBT state monitoring method based on finite element simulation Download PDFInfo
- Publication number
- CN112699579A CN112699579A CN202011462185.XA CN202011462185A CN112699579A CN 112699579 A CN112699579 A CN 112699579A CN 202011462185 A CN202011462185 A CN 202011462185A CN 112699579 A CN112699579 A CN 112699579A
- Authority
- CN
- China
- Prior art keywords
- vce
- bonding wire
- state
- igbt
- equation
- 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
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3308—Design verification, e.g. functional simulation or model checking using simulation
-
- 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)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Testing Of Individual Semiconductor Devices (AREA)
Abstract
The invention discloses an IGBT state monitoring method based on finite element simulation, which comprises the following steps: s1, obtaining historical data of the IGBT module Vce, on; s2, establishing a failure physical model of the bonding wire degradation process; s3, establishing a state space equation of the bonding wire degradation process; s4, determining the unknown quantity of Vce and on correction equation through a parameter learning algorithm; and S5, state estimation. The invention discloses the relation between the increase of the IGBT module Vce, on and the degradation of the bonding wire, establishes an accurate model of the degradation of the bonding wire of the IGBT module, and estimates the current crack length by utilizing an edge resampling moving algorithm based on particles so as to carry out state estimation. Compared with the traditional state monitoring method, the method has the advantages that time cost and economic cost are reduced, monitoring errors are reduced, and the reliability of the IGBT is improved.
Description
Technical Field
The invention relates to the field of semiconductors, in particular to an IGBT state monitoring method based on finite element simulation.
Background
Power electronic converters based on IGBTs (insulated gate bipolar transistors) are widely used in modern industry, including wind generators, electric vehicles, power transmission systems, etc. The reliability problem of the power electronic converter is also increasingly highlighted, and the power electronic components are one of the main failure sources of the system. From the viewpoint of safety and economy, the reliability of the converter, particularly the reliability of the power electronic components, has attracted much attention.
Package failures are one of the most common failures in IGBT modules, due to mismatch in the coefficients of thermal expansion of the different materials inside the assembly. When the IGBT is operated, periodic power losses occur, resulting in temperature fluctuations inside the power module, and significant thermo-mechanical stresses will be introduced at the joints between the different material layers, resulting in higher thermal stresses being experienced by the bond wires and solder joints. In practical applications, the degradation of the bonding wire and the fatigue of the solder layer are two main failure mechanisms, and both of them can reduce the performance of the IGBT module over time, and if no proper measures are taken, the device can be damaged, and the operation and safety of the whole system are affected.
Condition monitoring is an effective method for improving the reliability of the IGBT. The state monitoring is mainly to monitor the parameter index of the IGBT on line to know the current operation condition of the IGBT, and once an unhealthy state or fault is detected, further operation is carried out to avoid catastrophic consequences.
Vce, on (turn-on voltage between the collector and emitter of the IGBT) is the most commonly selected parameter, on which many state monitoring methods have been studied. Although the prior art has achieved certain success, certain challenges exist, such as that as the IGBT ages, the case temperature and the metal restructuring of the solder layer have a large effect on Vce, on measurement, and how to reduce these effects is a problem at present.
Disclosure of Invention
Vce and on are good indexes of state monitoring, a bonding wire is one of the parts of the IGBT module which are easy to damage, and in order to reveal the relation between the increase of the Vce and on of the IGBT module and the degradation state of the bonding wire and estimate the current state of the IGBT, the invention provides an IGBT state monitoring method based on finite element simulation. The method is realized from model establishment to algorithm implementation without carrying out a large number of power cycle tests, and compared with the traditional state monitoring method, the method has the advantages of reducing time cost and economic cost and improving the reliability of the IGBT.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that:
the IGBT state monitoring method based on finite element simulation comprises the following steps:
s1, obtaining historical data of the IGBT module Vce, on;
s2, establishing a failure physical model of the bonding wire degradation process;
s3, establishing a state space equation of the bonding wire degradation process;
s4, determining the unknown quantity of Vce and on correction equation through a parameter learning algorithm;
and S5, state estimation.
Further, the specific method of step S1 is:
and establishing a data acquisition unit, continuously acquiring Vce and on data of the IGBT, and storing the data.
Further, the specific method of step S2 is:
and obtaining the relation between the increase of the bonding wire Vce, on and the crack length of the bonding wire through finite element simulation.
Further, the specific method of step S3 includes the following sub-steps:
s3-1, selecting a Paris equation to describe crack propagation to obtain a state equation of the crack propagation, wherein the state equation is expressed in that the change rate of the crack length on the thermal cycle number and a stress intensity factor form an exponential relation, and the stress intensity factor can be obtained from a stress-strain curve of a bonding line;
s3-2, compensating the influence of temperature change on Vce, on, and correcting temperature quantity delta Vce,on,TExpressed as the product of the temperature coefficient and the current temperature minus the reference temperature;
s3-3, compensating the influence of metal reconstruction on Vce, on, and correcting quantity delta V of metal reconstructionce,on,MExpressed as:
ΔVce,on,M=1.6×10-17kVce,on,rΔT5.165
wherein k is an unknown parameter; vce,on,rIs a Vce,onA reference value of (d);delta T is the difference between the current temperature of the IGBT and the reference temperature;
s3-4, establishing a correction equation of Vce and on, wherein the corrected Vce and on are expressed as the superposition of reference Vce and on, temperature correction and metal reconstruction correction.
Further, the specific method of step S4 includes the following sub-steps:
s4-1, processing the sample by using a particle filter algorithm, and reducing the calculated amount;
s4-2, adjusting the weight of the sample, and approximating the target distribution by the weighted sample;
s4-3, estimating Vce by using a marginalized resampling motion algorithm based on particles, and correcting unknown parameters k of an equation in on mode.
Further, the specific method of step S5 includes the following sub-steps:
s5-1, bringing the determined parameter k back to a correction equation, calculating corrected Vce, on, and estimating the state a of the bond wire crack length, wherein the state a is expressed as the sum of products of each particle in the Vce, on state and a corresponding normalized weight value;
s5-2, calculating the stress intensity F of the bonding wire under the current state, wherein the stress intensity F is expressed as the product of the difference of the thermal expansion coefficients of the aluminum and the silicon and the square root of the state a of the crack length of the bonding wire;
and S5-3, judging whether the stress intensity exceeds the breaking stress intensity of the bonding wire, if so, indicating that the bonding wire reaches the failure standard, and recommending to immediately stop the machine for maintenance, otherwise, returning to the step S1.
The invention has the beneficial effects that: the invention provides an IGBT state monitoring method based on finite element simulation, a system model establishes an accurate model of IGBT module bonding wire degradation from a physical angle rather than simple curve fitting, and estimates the current crack length by using a particle-based edge resampling moving algorithm, thereby carrying out state estimation.
Drawings
FIG. 1 is a schematic flow diagram of the process.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 1, the IGBT state monitoring method based on finite element simulation includes the following steps:
s1, obtaining historical data of the IGBT module Vce, on;
s2, establishing a failure physical model of the bonding wire degradation process;
s3, establishing a state space equation of the bonding wire degradation process;
s4, determining the unknown quantity of Vce and on correction equation through a parameter learning algorithm;
and S5, state estimation.
The specific method of step S1 is:
and establishing a data acquisition unit, continuously acquiring Vce and on data of the IGBT, and storing the data.
The specific method of step S2 is:
and obtaining the relation between the increase of the bonding wire Vce, on and the crack length of the bonding wire through finite element simulation.
The specific method of step S3 includes the following substeps:
s3-1, selecting a Paris equation to describe crack propagation to obtain a state equation of the crack propagation, wherein the state equation is expressed in that the change rate of the crack length on the thermal cycle number and a stress intensity factor form an exponential relation, and the stress intensity factor can be obtained from a stress-strain curve of a bonding line;
s3-2, compensating the influence of temperature change on Vce, on, and correcting temperature quantity delta Vce,on,TExpressed as the product of the temperature coefficient and the current temperature minus the reference temperature;
s3-3, compensating the shadow of metal reconstruction on Vce, onSound and metal reconstruction correction quantity delta Vce,on,MExpressed as:
ΔVce,on,M=1.6×10-17kVce,on,rΔT5.165
wherein k is an unknown parameter; vce,on,rIs a Vce,onA reference value of (d); delta T is the difference between the current temperature of the IGBT and the reference temperature;
s3-4, establishing a correction equation of Vce and on, wherein the corrected Vce and on are expressed as the superposition of reference Vce and on, temperature correction and metal reconstruction correction.
The specific method of step S4 includes the following substeps:
s4-1, processing the sample by using a particle filter algorithm, and reducing the calculated amount;
s4-2, adjusting the weight of the sample, and approximating the target distribution by the weighted sample;
s4-3, estimating Vce by using a marginalized resampling motion algorithm based on particles, and correcting unknown parameters k of an equation in on mode.
The specific method of step S5 includes the following substeps:
s5-1, bringing the determined parameter k back to a correction equation, calculating corrected Vce, on, and estimating the state a of the bond wire crack length, wherein the state a is expressed as the sum of products of each particle in the Vce, on state and a corresponding normalized weight value;
s5-2, calculating the stress intensity F of the bonding wire under the current state, wherein the stress intensity F is expressed as the product of the difference of the thermal expansion coefficients of the aluminum and the silicon and the square root of the state a of the crack length of the bonding wire;
and S5-3, judging whether the stress intensity exceeds the breaking stress intensity of the bonding wire, if so, indicating that the bonding wire reaches the failure standard, and recommending to immediately stop the machine for maintenance, otherwise, returning to the step S1.
In the specific implementation process, the IGBT module in a certain working state in the power electronic system is subjected to fault diagnosis, a bonding wire degradation model is established by acquiring historical data of Vce and on, the crack length and the stress intensity of the bonding wire are calculated, whether the stress intensity exceeds the fracture stress intensity of the IGBT bonding wire of the type is judged, if yes, the crack of the bonding wire is expanded to enter an unstable area, and the bonding wire is indicated to have a fault.
In summary, an IGBT state monitoring method based on finite element simulation is proposed herein. Vce and on are good indexes for state monitoring, the relation between the increase of the bonding wire Vce and on and the degradation state of the bonding wire can be revealed through finite element simulation, and an accurate bonding wire degradation model is established from the theory of a fracture mechanism. Compared with the traditional state monitoring method, the method does not need to carry out a large number of power cycle tests, can greatly reduce time cost and economic cost, has small monitoring error, and improves the reliability of power electronic components and even the whole power electronic system.
Claims (6)
1. The IGBT state monitoring method based on finite element simulation is characterized by comprising the following steps of:
s1, obtaining historical data of the IGBT module Vce, on;
s2, establishing a failure physical model of the bonding wire degradation process;
s3, establishing a state space equation of the bonding wire degradation process;
s4, determining the unknown quantity of Vce and on correction equation through a parameter learning algorithm;
and S5, state estimation.
2. The finite element simulation-based IGBT state monitoring method according to claim 1, wherein the specific method of step S1 is as follows:
and establishing a data acquisition unit, continuously acquiring Vce and on data of the IGBT, and storing the data.
3. The finite element simulation-based IGBT state monitoring method according to claim 1, wherein the specific method of step S2 is as follows:
and obtaining the relation between the increase of the bonding wire Vce, on and the crack length of the bonding wire through finite element simulation.
4. The finite element simulation-based IGBT state monitoring method according to claim 1, wherein the specific method of step S3 comprises the following sub-steps:
s3-1, selecting a Paris equation to describe crack propagation to obtain a state equation of the crack propagation, wherein the state equation is expressed in that the change rate of the crack length on the thermal cycle number and a stress intensity factor form an exponential relation, and the stress intensity factor can be obtained from a stress-strain curve of a bonding line;
s3-2, compensating the influence of temperature change on Vce, on, and correcting temperature quantity delta Vce,on,TExpressed as the product of the temperature coefficient and the current temperature minus the reference temperature;
s3-3, compensating the influence of metal reconstruction on Vce, on, and correcting quantity delta V of metal reconstructionce,on,MExpressed as:
ΔVce,on,M=1.6×10-17k·Vce,on,r·ΔT5.165
wherein k is an unknown parameter; vce,on,rIs a Vce,onA reference value of (d); delta T is the difference between the current temperature of the IGBT and the reference temperature;
s3-4, establishing a correction equation of Vce and on, wherein the corrected Vce and on are expressed as the superposition of reference Vce and on, temperature correction and metal reconstruction correction.
5. The finite element simulation-based IGBT state monitoring method according to claim 1, wherein the specific method of step S4 comprises the following sub-steps:
s4-1, processing the sample by using a particle filter algorithm, and reducing the calculated amount;
s4-2, adjusting the weight of the sample, and approximating the target distribution by the weighted sample;
s4-3, estimating Vce by using a marginalized resampling motion algorithm based on particles, and correcting unknown parameters k of an equation in on mode.
6. The bonding wire-based IGBT life prediction method according to claim 1, characterized in that the specific method of step S5 comprises the following sub-steps:
s5-1, bringing the determined parameter k back to a correction equation, calculating corrected Vce, on, and estimating the state a of the bond wire crack length, wherein the state a is expressed as the sum of products of each particle in the Vce, on state and a corresponding normalized weight value;
s5-2, calculating the stress intensity F of the bonding wire under the current state, wherein the stress intensity F is expressed as the product of the difference of the thermal expansion coefficients of the aluminum and the silicon and the square root of the state a of the crack length of the bonding wire;
and S5-3, judging whether the stress intensity exceeds the breaking stress intensity of the bonding wire, if so, indicating that the bonding wire reaches the failure standard, and recommending to immediately stop the machine for maintenance, otherwise, returning to the step S1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011462185.XA CN112699579A (en) | 2020-12-09 | 2020-12-09 | IGBT state monitoring method based on finite element simulation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011462185.XA CN112699579A (en) | 2020-12-09 | 2020-12-09 | IGBT state monitoring method based on finite element simulation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112699579A true CN112699579A (en) | 2021-04-23 |
Family
ID=75507591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011462185.XA Pending CN112699579A (en) | 2020-12-09 | 2020-12-09 | IGBT state monitoring method based on finite element simulation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112699579A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115081268A (en) * | 2022-05-22 | 2022-09-20 | 北京化工大学 | IGBT module service life prediction method based on power loss mathematical model and finite element method |
CN117195665A (en) * | 2023-11-07 | 2023-12-08 | 湖南大学 | Method for constructing life prediction model of bonding wire of power semiconductor device |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861040A (en) * | 2016-09-22 | 2018-03-30 | 北京航空航天大学 | A kind of IGBT intermittent life test methods based on simulation modeling and short-term test |
CN108287976A (en) * | 2018-03-02 | 2018-07-17 | 电子科技大学 | A kind of screening section verification methods of the IGBT based on faulty physical and finite element simulation |
CN108509741A (en) * | 2018-04-10 | 2018-09-07 | 西南科技大学 | A kind of finite element numerical method for solving of mudstone flow equation |
CN109241609A (en) * | 2018-08-31 | 2019-01-18 | 华东交通大学 | A kind of Bayesian dynamic prediction method based on markov chain Monte-Carlo |
CN109738773A (en) * | 2018-06-19 | 2019-05-10 | 北京航空航天大学 | IGBT module life-span prediction method under a kind of non-stationary operating condition |
CN110298126A (en) * | 2019-07-04 | 2019-10-01 | 北京航空航天大学 | A kind of polynary Copula power device method for evaluating reliability based on the physics of failure |
CN111007379A (en) * | 2019-12-27 | 2020-04-14 | 电子科技大学 | Self-correcting IGBT health monitoring method |
CN111783243A (en) * | 2020-06-18 | 2020-10-16 | 东南大学 | Metal structure fatigue crack propagation life prediction method based on filtering algorithm |
-
2020
- 2020-12-09 CN CN202011462185.XA patent/CN112699579A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861040A (en) * | 2016-09-22 | 2018-03-30 | 北京航空航天大学 | A kind of IGBT intermittent life test methods based on simulation modeling and short-term test |
CN108287976A (en) * | 2018-03-02 | 2018-07-17 | 电子科技大学 | A kind of screening section verification methods of the IGBT based on faulty physical and finite element simulation |
CN108509741A (en) * | 2018-04-10 | 2018-09-07 | 西南科技大学 | A kind of finite element numerical method for solving of mudstone flow equation |
CN109738773A (en) * | 2018-06-19 | 2019-05-10 | 北京航空航天大学 | IGBT module life-span prediction method under a kind of non-stationary operating condition |
CN109241609A (en) * | 2018-08-31 | 2019-01-18 | 华东交通大学 | A kind of Bayesian dynamic prediction method based on markov chain Monte-Carlo |
CN110298126A (en) * | 2019-07-04 | 2019-10-01 | 北京航空航天大学 | A kind of polynary Copula power device method for evaluating reliability based on the physics of failure |
CN111007379A (en) * | 2019-12-27 | 2020-04-14 | 电子科技大学 | Self-correcting IGBT health monitoring method |
CN111783243A (en) * | 2020-06-18 | 2020-10-16 | 东南大学 | Metal structure fatigue crack propagation life prediction method based on filtering algorithm |
Non-Patent Citations (4)
Title |
---|
KETING HU等: "Cost-Effective Prognostics of IGBT Bond Wires with Consideration of Temperature Swing", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 * |
岳亚静等: "基于结温补偿的IGBT模块键合线失效监测方法", 《电力电子技术》 * |
李志星等: "基于PFA的IGBT键合线失效机理及寿命预测", 《半导体制造技术》 * |
陈健等: "基于高斯权值-混合建议分布粒子滤波的疲劳裂纹扩展预测", 《航空学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115081268A (en) * | 2022-05-22 | 2022-09-20 | 北京化工大学 | IGBT module service life prediction method based on power loss mathematical model and finite element method |
CN115081268B (en) * | 2022-05-22 | 2024-08-06 | 北京化工大学 | IGBT module service life prediction method based on power loss mathematical model and finite element method |
CN117195665A (en) * | 2023-11-07 | 2023-12-08 | 湖南大学 | Method for constructing life prediction model of bonding wire of power semiconductor device |
CN117195665B (en) * | 2023-11-07 | 2024-01-26 | 湖南大学 | Method for constructing life prediction model of bonding wire of power semiconductor device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10705133B2 (en) | Method and device for estimating level of damage or lifetime expectation of power semiconductor module | |
Wang et al. | Power electronics reliability: State of the art and outlook | |
CN109738773B (en) | IGBT module service life prediction method under non-steady working condition | |
US11248966B2 (en) | Health monitoring and failure prognosis of power electronics devices | |
CN110161396A (en) | The health evaluating based on thermal model of IGBT | |
CN112699579A (en) | IGBT state monitoring method based on finite element simulation | |
CN112214951B (en) | Sectional type LSTM microgrid inverter IGBT reliability evaluation method and system | |
CN109885884B (en) | IGBT module fatigue analysis processing method and semiconductor device processing method | |
CN107508281B (en) | Dynamic reliability assessment method for power flow controller of in-phase power supply system | |
Wei et al. | Condition monitoring of IGBT modules based on changes of thermal characteristics | |
CN111398766A (en) | IGBT module health state online monitoring method and system | |
CN105550397B (en) | A kind of IGBT module state evaluating method based on damage voltage | |
Ni et al. | Review of SiC MOSFET based three-phase inverter lifetime prediction | |
CN112487651B (en) | Service life detection method of photovoltaic converter power device | |
CN111060797A (en) | IGBT module health state monitoring method based on natural frequency of heat network | |
CN112986784B (en) | Abnormity identification method and device for high-power welding type IGBT module | |
Liu et al. | BP neural network for non-invasive IGBT junction temperature online detection | |
CN117932915A (en) | Method for predicting residual life of power device | |
CN115081268B (en) | IGBT module service life prediction method based on power loss mathematical model and finite element method | |
CN115563754A (en) | IGBT module aging characteristic decoupling parameter monitoring method based on temperature characteristic | |
CN115600423A (en) | Motor controller service life assessment method | |
Hong et al. | Prognosis for insulated gate bipolar transistor based on Gaussian process regression | |
Liu et al. | Reliability evaluation of IGBT power module on electric vehicle using big data | |
CN113919116A (en) | IGBT remaining service life prediction method based on GARCH model | |
Ahmedi et al. | Lifetime estimation of IGBT power modules for reliability study of wind turbine systems |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210423 |