Next Article in Journal
Baseline-Free Structural Damage Identification for Beam-Like Structures Using Curvature Waveforms of Propagating Flexural Waves
Previous Article in Journal
Privacy-Preserving IoT Data Aggregation Based on Blockchain and Homomorphic Encryption
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on Real Time IGBT Junction Temperature Estimation Using the NTC and Calculation of Power Losses in the Automotive Inverter System

Electric Motor Control Laboratory, Graduate School of Automotive Engineering, Kookmin University, Seoul 02707, Korea
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(7), 2454; https://doi.org/10.3390/s21072454
Submission received: 11 February 2021 / Revised: 11 March 2021 / Accepted: 24 March 2021 / Published: 2 April 2021
(This article belongs to the Section Vehicular Sensing)

Abstract

:
This paper proposes a junction temperature estimation algorithm for the insulated gate bipolar transistor (IGBT) based on a power loss calculation and a thermal impedance model for inverter systems. The Simulink model was designed to calculate the power losses of power semiconductor devices and to estimate the junction temperature with a simplified thermal impedance model. This model can estimate the junction temperature up to the transient state, including the steady state. The parameters used to calculate the power losses, the thermal resistance, and the thermal capacitance were optimized for a given inverter to be tested for improving the accuracy. The simulation results and experimental measurement data were compared to verify the proposed junction temperature estimation algorithm. Finally, the algorithm was installed on the inverter controller, and the performance was verified by comparing the real time estimation result with the measured temperature.

1. Introduction

The temperature management strategy of an insulated gate bipolar transistor (IGBT), which is the power module of an inverter, has become an important issue with the increasing development of inverters for driving traction motors in electric vehicle application. Generally, current IGBTs are limited by the junction temperature, and the possibility of failure due to overheating tends to greatly increase unless the junction temperature is properly managed [1,2,3]. Recently, there have been numerous studies on methods for measuring the temperature of the IGBT indirectly by configuring the hardware to manage the junction temperature of the IGBT. One method measures temperature by attaching a negative temperature coefficient thermistor (NTC) near the junction. Most of the packaging modules provide a built-in NTC, but there is thermal resistance and thermal capacitance due to the physical separation distance, resulting in a difference from the junction temperature. Accordingly, correlation between the NTC temperature and the junction temperature is required [1,4]. Another method indirectly measures the junction temperature by monitoring the collector-emitter voltage during the on-state of the IGBT (i.e., the saturation voltage of IGBT). Another example is Mitsubishi’s J1 series. Companies such as semiconductor manufacturers provide a function to check the junction temperature by attaching a temperature sensor at a very close position to the junction [5]. In addition to the methods indirectly measuring the HW as described above, methods of estimating the junction temperature through thermal models have also been studied [6].
Previous studies are largely divided into analytical models [7,8,9], numerical models [10,11,12], and thermal network models [13,14,15,16,17,18] according to the method of designing the thermal model [19]. The thermal network models covered in this paper are divided into the Cauer and Foster thermal impedance models [20]. The Cauer thermal impedance model is based on the correlation between the IGBT module packaging structure and the physical property, and the Foster thermal impedance model is an equivalent model fitted from measurement results or simulation results. In the case of the thermal network model, valid results can be expected only within a limited range. It is difficult to consider temperature-dependent variables under conditions such as wide-range ambient temperature variations or transient conditions during operation at overloads or short-circuits [3]. Accordingly, many studies have been conducted to design an effective thermal impedance model. For example, there are studies about correcting the thermal impedance model considering temperature-dependent variables [3] or finite element method (FEM) thermal impedance models designed based on the physical structure of the IGBT module [3,21]. These studies improve the reliability of the estimated junction temperature. However, because most studies are simulation-based, the more complex the thermal impedance model, the more difficult it is to implement an algorithm that can estimate the junction temperature in the inverter in real time, although the complex thermal impedance model has higher accuracy.
For this reason, in the field of mass-producing inverters for driving automobiles, to manage the junction temperature in real time when driving the inverter, the relationship between NTC and junction temperature is functionalized or tabled for various driving conditions to find the junction temperature. Due to the limitation of the memory of the inverter controller and the limitation of the time resource required to obtain experimental data, the area in which the junction temperature can be checked is limited to within a few hundred rpm compared to the operating range of the inverter for driving a vehicle from 12,000 to 20,000 rpm. In addition, even if the relationship between NTC and junction temperature is functionalized or tabled, it is difficult to respond to rapid changes in ambient temperature or changes in cooling conditions. This naturally leads to requiring a design margin according to the junction temperature when selecting the specifications of the IGBT module, which is directly connected to the increase in the inverter cost.
Therefore, this study focused on overcoming these limitations and designing a junction temperature estimation algorithm that can take into account variations in ambient and coolant temperatures, and can be used in inverter controllers. It also aimed at ensuring the junction temperature estimation performance over the entire operating range of the inverter. The algorithm designed for this should be as simple as possible while ensuring accuracy for various load patterns. First, the method of calculating the instantaneous power losses for each pulse width modulation (PWM) period is described in the current paper using the operating parameters of the IGBT provided in the datasheet and the variables sampled in real time by the inverter controller. Then, a simplified thermal impedance model is presented based on the NTC temperature of the IGBT module sampled by the inverter controller. The simulation results of the proposed algorithm are compared with the measured data to verify the reliability. Finally, the algorithms were implemented in the C language and applied to the inverter controller. The reliability of the estimation results was then verified by comparison with the measured temperature by performing under the load experiment on the dynamometer.

2. Calculation Power Losses of IGBT Modules

The IGBT module used in this study is a FF1000R17IE4 shown in Figure 1, which provides a package for one phase. Specifically, a substrate with two IGBTs and two freewheeling diodes is arranged in six parallel structures. The power losses of the IGBT module can be explained by the sum of the IGBT power losses and the diode power losses:
P t o t a l = P I G B T + P D i o d e
The power losses of IGBT in Equation (1) can be divided into the switching losses occurring in the transient state (i.e., turn-on and turn-off state) and the conduction losses caused by the saturation voltage during the on-state of IGBT:
P t o t a l [ I G B T ] = P s w [ I G B T ] + P c o n d I G B T
P c o n d I G B T = V c e s a t · I c
P s w I G B T = E o n + E o f f · f s w
Like the IGBT, the power losses of the diode can be divided into the conduction losses and the switching losses:
P t o t a l [ D i o d e ] = P s w [ D i o d e ] + P c o n d [ D i o d e ]
P c o n d D i o d e = V f · I f
P s w D i o d e = E r e c · f s w
In this study, the parameters related to the power losses on the IGBT datasheet were optimized for application. The method of optimizing the power loss parameters provided in the datasheet was applied to improve the accuracy of the power losses calculated in real time [22]. Optimized power loss parameters are provided as look-up tables, allowing parameters to be calculated in real time according to dependent variables (e.g., current, voltage, and junction temperature). The proposed method can calculate the power losses with the sampling variables (e.g., the terminal voltage, the phase current, the duty of PWM, and the switching frequency) according to the driving conditions and the optimized parameters of power losses. This proposed method was implemented as a Simulink model. The simulation results were compared with the experimental measurement data, and its validity is proved. A block diagram of the proposed method is shown in Figure 2.

2.1. Optimizing the Operating Parameters of IGBT

The collector-emitter saturation voltage of the IGBT depends on the collector-emitter voltages during the turn-off state of the IGBT, the collector current flowing through the IGBT, the junction temperature, and a given gate-emitter voltage. Semiconductor manufacturers already provide the saturation voltage ( V c e s a t ) characteristic curve in this regard. Therefore, based on this characteristic curve, the saturation voltage ( V c e s a t ) optimized for a given inverter to be tested was recalculated and applied as a look-up table. The Simulink model was constructed to calculate the saturation voltage ( V c e s a t ) according to the collector current flowing through the IGBT and junction temperature in real time. The saturation voltage ( V c e s a t ) optimized by recalculation using the saturation voltage ( V c e s a t ) characteristic curve can be expressed as follows:
V c e s a t . T j I = V c e s a t . T j · V c e s a t . V g e _ I n v V c e s a t . V g e _ n o m · V d c . I n v V d c . n o m
P c o n d I G B T = V c e s a t . T j I · I c
Switching losses of the IGBT depend on the collector current flowing through the IGBT, the junction temperature, the given gate resistor, the gate-emitter voltages, and the collector-emitter voltages. The optimized switching losses ( E o n , E o f f ) were recalculated using the characteristic curve of switching losses ( E o n , E o f f ) provided in the datasheet. The Simulink model was then implemented using the look-up table to calculate the switching losses ( E o n , E o f f ) according to the collector current flowing through the IGBT and the junction temperature in real time. The recalculated switching losses ( E o n , E o f f ) are given by the following Equations (10) and (11). The switching losses of the IGBT are finally calculated as Equation (12) and the results are shown in Figure 3.
E o n . T j I = E o n . T j I · E o n . R g _ I n v E o n . R g _ n o m · V d c . I n v V d c . n o m
E o f f . T j I = E o f f . T j I · E o f f . R g _ I n v E o f f . R g _ n o m · V d c . I n v V d c . n o m
P s w I G B T = E o n + E o f f · f s w

2.2. Optimizing the Operating Parameters of Diode

The conduction losses of the diode are caused by the voltage drop ( V f ) during the on-state of the diode, which depends on the current flowing through the diode and the temperature of the junction. Therefore, the characteristic curve of the voltage drop ( V f ) based on the current and the junction temperature is applied via the look-up table to calculate the voltage drop ( V f ) according to the driving conditions. The voltage drop ( V f ) is shown as Equation (13). The conduction losses of the diode are given by Equation (14).
V f . T j I = V f . T j I
P c o n d D i o d e = V f . T j I · I f
Switching losses of the diode are caused by reverse current in the reverse recovery area of the diode during turn-off. This is dependent on the value of the gate resistor for turn-on, the collector-emitter voltages during turn-off state, and the current flowing through the diode at the forward bias state. In addition, the switching losses energy ( E r e c ) optimized for the given inverter was recalculated by utilizing the characteristic curve of the switching losses energy ( E r e c ) provided in the datasheet and applied as the look-up table, so it can be calculated according to the current and the junction temperature in real time. The optimized switching losses ( E r e c ) is given by Equations (15) and (16) and the results are shown in Figure 4.
E r e c . T j I = E r e c . T j I · E r e c . R g _ I n v E r e c . R g _ n o m · V d c . I n v V d c . n o m
P s w D i o d e = E r e c · f s w

2.3. The Simulation Result of IGBT Power Loss Calculation

To verify the algorithm that calculates the power losses in real time by applying the optimized power loss parameters as look-up tables, the inverter composed with the black paint module of IGBT was tested in the dynamometer environment. The data of power losses according to operating conditions were measured using a power analyzer (i.e., WT3000), and the simulation was performed by applying the same operating conditions. Considering the junction temperature over time to enhance the accuracy of simulation, the measured data with an IR camera were applied in the simulation. The power loss data obtained on the power analyzer were calculated using the root mean square (RMS) value as the current unit. The operating condition was thus selected as a stall condition to reduce errors due to the limitations of the measuring equipment. Other conditions were given randomly to verify the reliability of the power losses calculated within the wide operating range. The error of each experiment was calculated as the difference between the result of the simulation and the measured data. As a result, it was verified that the maximum error was 1.22%, which indicates high accuracy. Thus, it is possible to obtain the power losses for estimating the junction temperature using the optimized power loss parameters look-up tables and variables (e.g., the phase current, the PWM duty, and the DC-link voltage) of the inverter controller. The simulation results for each experiment are shown in Table 1.

3. Design of the Simplified Thermal Impedance Model Using NTC

The described research results are based on the experimental method to design the Foster model. The thermal impedance model proposed in this paper is the thermal impedance model that has an NTC, which is packaged in an IGBT module to provide the temperature of the baseplate as the ground node, rather than a model that has the coolant temperature (or ambient temperature) as the ground node. In the case of the existing thermal circuit model, it is common to select a node based on a vertical structure ranging from coolant to thermal grease, baseplate, solder, substrate, solder, and junction. This can, however, differ largely from the actual thermal impedance because the factors to be considered in the design process are complex and it is not easy to obtain accurate data. In addition, it is difficult to consider the effect of the change in temperature, and flow rate of the cooling water and the thermal grease. In the case of the NTC, it is affected by the coolant temperature (or ambient temperature) to the same degree as the junction of the IGBT module. This means that this model has an advantage of being able to react indirectly to a changing coolant temperature (or ambient temperature). When the NTC is equalized to the junction with one thermal impedance, the effect of cooling water and thermal grease is affected by both the NTC and the junction, such that the influence can be indirectly reflected by the temperature change of the NTC. In this case, even if there is a difference in the influence of the junction and the NTC on the cooling water and the thermal grease, it can be assumed that the error in the thermal impedance caused by the difference is very small compared to the thermal impedance due to the separation distance between the NTC and the junction. In addition, because both the NTC and the junction are located on the substrate, the temperature data can be easily measured. To implement the C language-based algorithm, the thermal impedance model was designed with one thermal resistance and one thermal capacitance by simplifying it as much as possible. The proposed thermal impedance model is shown in Figure 5. Figure 5a is the proposed circuit model shown with a schematic cross section of an IGBT module. A heatsink with a water jacket is a commonly used cooling structure in the traction inverter. The baseplate is attached on the heatsink with thermal compound. The substrate is soldered onto the baseplate. The IGBT and the diode is soldered onto the substrate and the NTC is also soldered. The packaged circuit is connected with bond wire. The thermal impedance between the NTC and the junction are expressed as just one thermal impedance value. Figure 5b shows the simplified thermal circuit diagram. The ground node sets with temperature of the NTC.
The experimental data were obtained by measuring the NTC temperature and the junction temperature under actual load conditions according to the reference method for the Foster model [17] to calculate the thermal resistance and the thermal capacitance of the proposed thermal impedance model. The temperature data of the two measuring points are applied to Equation (17), which was used to obtain the thermal impedance over time. Finally, the thermal impedance model of the junction and the NTC can be obtained by extracting the thermal resistance and the thermal capacitance through curve fitting with the thermal impedance obtained over time. The form of curve fitting is equal to Equation (18).
Z t h _ j N t = Δ T j N t P l o s s = T j u n c t i o n t T N T C t P l o s s
Z t h _ j N t = R t h _ j N ( 1 e t R t h _ j N · C t h _ j N )

3.1. The Thermal Impedance Model between Junction and NTC

To design the thermal impedance model from the junction to the NTC, the temperature of the NTC and the junction temperature were acquired in real time under the load test conditions in a dynamometer environment. The temperature of the NTC was measured with a temperature measuring instrument with a sampling time of 10 ms. The temperature of the junction was measured with an IR camera using the black paint module shown in Figure 6. To obtain an accurate thermal impedance, it is first important to measure each temperature accurately. The purpose of this study is to design the simplest thermal impedance model to secure junction temperature estimation performance, so the process of designing thermal impedance is very important. To ensure the junction temperature estimation performance in all operating areas, various experimental cases were selected. The load test conditions were selected within the operating range of the inverter required by the system. Each experimental condition was selected according to the DC-link voltage (V), the switching frequency (kHz), the coolant temperature (°C), load (%), and fundamental frequency (Hz). The FF1000R17IE4 used in this experiment has a 6-parallel structure. As shown in Figure 7, six high-side IGBTs and six low-side diodes are arranged on the upper side, and six low-side IGBTs and six high-side diodes are arranged on the lower side. Therefore, the thermal impedance model for stall condition and driving condition were independently extracted because the thermal impedance has a difference according to whether the fundamental frequency in the operating condition exists. The remainder of the experiment conditions were randomly selected and conducted within the operating range of the inverter. A total of 40 experiments were conducted according to each condition to obtain experimental data. The thermal impedance graph was obtained using the temperature data measured for each experimental condition through Equation (17). The thermal resistance and thermal capacitance were then extracted by fitting the thermal impedance graph. It was confirmed that the standard deviation of the calculated thermal resistance and thermal capacitance for all experimental cases was 1.16 for the constraint condition and 0.17 for the driving condition. Therefore, the final thermal resistance and thermal capacitance were selected as average values. The results in Table 2 were selected as average values. In the stall condition, the average thermal resistance of the IGBT was 41.13 (°C/kW) and the average thermal resistance of the diode was 102.1 (°C/kW). The average thermal capacitance of the IGBT was 11.21 (W/°C) and the average thermal capacitance of the diode was 3.36 (W/°C). This is because the area of the IGBT is bigger than the area of the diode. In the driving condition, the average thermal resistance of the IGBT was 54.76 (°C/kW) and the average thermal resistance of the diode was 155.59 (°C/kW). The average thermal capacitance of the IGBT was 13.56 (W/°C) and the average thermal capacitance of the diode was 5.16 (W/°C). The reason why the result of the driving condition is larger than that of the stall condition is due to the switching pattern along the current path. In the stall condition, the inverter’s switching state is fixed for DC current and heat is generated in six IGBTs and six diodes. However, in the driving condition, the inverter’s switching state changes for AC current and heat is generated in 12 IGBTs and 12 diodes.

3.2. The Simulation Result of Junction Temperature Estimation

Simulation was performed by implementing the thermal impedance model. The thermal resistance and the thermal capacitance obtained through the experiment were applied in Simulink. The temperature of the NTC was applied as raw data experimentally acquired for each operating condition, and the power losses were calculated from the power loss calculation model designed in Section 2 and applied as the input. There is a disadvantage in that the error present in the power loss calculation model is reflected as the error of the estimated junction temperature, but simulation verification was performed in consideration of this to confirm the reliability of the junction temperature estimation algorithm during the final implementation stage. The thermal impedance model implemented in Simulink is shown in Figure 7, and the simulation results for the operating conditions are shown in Table 3. The error of the estimated temperature for each case was calculated based on the maximum temperature over time, and the maximum error was 1.5% for the IGBT and 1.3% for the diode. The high accuracy of these results indicates that the designed thermal impedance value was close to the actual impedance value and also proved that even a simple thermal impedance model with the NTC as a ground node can obtain sufficient estimation performance.

4. Experimental Verification of the Junction Temperature Estimation Algorithm

The power loss calculation algorithm verified through simulation and the thermal impedance model based on the NTC must be implemented as C language-based software to mount the inverter controller. When designing software, an operation period must be selected in consideration of the execution time, and an error occurs when there is a decrease in resolution according to the operation period. In addition, an error that is not confirmed in the simulation may additionally occur due to the potential error of the sampling variables used in the calculation. Therefore, it is necessary to verify whether the software is properly designed by comparing the result of the junction temperature estimation performed by the micro controller unit (MCU) and the actual temperature data through an experiment. In the experimental case, to check the temperature estimation performance for the overall driving range of the inverter, the entire driving speed range including the stall condition was considered, and the load pattern was also selected in various ways.

4.1. Experimental Setup

A black paint module was set up on the inverter shown in Figure 8. The junction temperature was measured by the IR camera. The data of the inverter (the estimated junction temperature, the calculated power losses, and the temperature of the NTC) were obtained by the controller area network (CAN) Analyzer with a sampling time of 100 ms.
The designed software consists of two functions: a power loss calculation algorithm and a thermal impedance model-based junction temperature estimation algorithm. The power loss parameter of the power loss calculation algorithm was applied as a lookup table and calculated through interpolation according to the driving conditions. At this time, the execution time could not be converged within the PWM period due to interpolation, so the operation period was set to 2 ms. Because the operation time was longer than the PWM period, the resolution was lowered, which in turn caused an error in the power loss calculation result. Finally, it was reflected as an error of the estimated temperature. To reduce the error due to the software limitation, the integral value of the power losses per 2 ms were applied in the stall condition because the integral value of the power losses due to the instantaneous current was the same as the integral value of the power losses due to the average current of 2 ms. By comparison, in the driving condition, as the frequency increases, an error increases between the integral value of the power losses due to the instantaneous current and the integral value of the power losses due to the average current of 2 ms. To solve this problem, the software was designed to predict the instantaneous power loss value by multiplying the sine wave of the fundamental frequency by the maximum power loss value.
Experimental conditions were selected with the sequence of torque command, the switching frequency, the fundamental frequency, the temperature of coolant, and the DC-link voltage. The block diagram of the proposed junction estimation algorithm is shown in Figure 9. The algorithm has two main functions. One is the thermal impedance model consisting of the thermal resistance and the thermal capacitance from the junction to the NTC. This function needs the results of the power loss calculation function, with each IGBT and diode as the inputs. The other is the power loss calculation function. This function needs the inputs that have the information for power loss parameters (e.g., the phase current, the duty of the PWM, the temperature of the NTC, the switching frequency, and the feedback, as a result of the estimated junction temperature). Finally, the results of this function will be the power losses of the IGBT and diode. These results will be fed back to the function of the thermal impedance model.

4.2. Experimental Results

The experimental results were classified according to the presence or absence of the fundamental frequency and analyzed by dividing them into the stall condition and the driving condition.

4.2.1. The Results of the Experiment at the Stall Condition

Table 4 shows the experiment results. For each experimental case, the error between the estimated value and the measured value of each case was calculated based on the maximum temperature over all data. As a result, the error was up to 2%. A value of 2% means that the maximum error is very similar to that of the simulation result, which was 1.3% for the IGBT and 1.5% for the diode. To mount the designed algorithm on the inverter controller, the software must be designed in consideration of the problem that the operation cycle cannot be shortened due to the execution time problem.
In addition, when considering the potential errors of the sampled variables, an error of 2% proves that the software is designed properly. Figure 10 shows graphs with the result of the experiment (i.e., the estimated temperature, the measured temperature, and the temperature feedback from the NTC) for Case 1.
With this result, it can be known that the slope at which the temperature of the NTC rises in the transient period during which the current rises is lower than the junction temperature. The junction temperature has a very steep rising slope. This means that it takes time for heat to be conducted due to the physical separation distance from the IGBT and diode chip, which are the heating points, to the NTC. In addition, the temperature of the junction can be two times the temperature of the NTC. This makes it difficult to manage the junction temperature only by observing the temperature of the NTC.
By comparison, the junction temperature estimated by the proposed algorithm converges with very little error in the entire section compared to the actual temperature. Furthermore, the estimated junction temperature converges to the actual temperature even in the state where the current does not flow and the temperature drops to the steady state. After cooling sufficiently, the junction temperature converges close to the temperature of the NTC. The reason why the temperature of the NTC is slower than the cooling rate seems to be the difference due to the shape of the channel, where the coolant flows directly under the baseplate where the IGBT and diode chip are located. The fact that the junction temperature converges to the temperature of the NTC when no current flows means that it is reasonable to assume that the junction temperature and the temperature of the NTC are affected equally by the ambient and the coolant temperature because the junction and NTC are on the same layer. Therefore, by designing the thermal impedance model based on the NTC, the external factors affecting the junction temperature can be considered indirectly.

4.2.2. The Results of the Experiment at the Driving Condition

The specifications of this traction motor-inverter system are listed in Table 5. Table 6 shows the experiment results in the driving condition. As in the stall condition, the error for each experiment was calculated based on the maximum temperature. The maximum error was 5.46%. For Cases 1, 2, and 3, the junction temperature was estimated while driving along the speed-torque (S-T) curve of the interior permanent magnet synchronous motor (IPMSM). The S-T curve ranged from 100 to 12,000 rpm based on the maximum output torque required by the system. Other test cases had a fixed speed. Through these experiments, it was possible to verify the performance of the proposed algorithm while operating with the maximum load applied up to the maximum fundamental frequency required for the inverter in the automotive traction system. The estimation results for Case 1 are shown in Figure 11. The results of estimating the junction temperature while operating along the S-T curve show a larger error than the result when the rpm is fixed at a low speed. Looking at Figure 11, it can be seen that a certain error continues to exist in the region where the fundamental frequency is high. The errors that occurred mostly fell within 3%. In addition to the transient region in which the junction temperature rises rapidly, the result of tracking the behavior of the actual junction temperature in the entire region can be confirmed. When the proposed algorithm is applied, it is advantageous in that the thermal behavior of the junction temperature can be estimated with high accuracy in the inverter controller, even when the required operating range is wide, as in an automotive traction system.

5. Conclusions

The power loss calculation algorithm of this paper calculates the instantaneous power losses for each switching period, which were used as an input to the thermal impedance model. Therefore, the junction temperature can be estimated in real time and in the steady state or the transient state using the proposed algorithm. In addition, the thermal impedance model of the NTC to the junction has meaningful results in estimating the junction temperature. These results were verified experimentally by comparing the estimated and measured temperatures. The error of the constraint condition was within 2%, and the error of the driving condition was confirmed to be within 6%. This is proof that the method of designing the power loss calculation algorithm and the temperature estimation algorithm using the NTC-based thermal circuit model is reasonable. In particular, it showed that it is possible to achieve not only a similar performance to the generally thermal impedance model, but also practicality. Above all, the simplified thermal circuit model based on the NTC calculated each thermal impedance for various load patterns during the design process, and the thermal impedance was selected as an average value taking into account the deviation of each case. As a result, it was shown that the temperature estimation performance can be secured by showing a narrow error distribution, with only the initial design values of the thermal circuit model without a separate correction process. Therefore, it can be applied to mass production sufficiently if additional verification is performed in consideration of the error distribution that may occur when applying the temperature estimation algorithm to mass-produced products. In addition, compared to conventional methods, the design margin due to junction temperature can be reduced, which is expected to reduce inverter costs. On this point, the proposed algorithm of this paper would be useful in estimating the junction temperature of a water-cooled inverter for automotive systems in real time, and preventing the failure of power modules by overly high temperatures.

Author Contributions

Conceptualization, H.L. and G.L.; data curation, H.L., J.H. and S.K.; formal analysis, H.L., S.K. and H.B.; investigation, H.L. and J.H.; methodology, H.L.; project administration, G.L.; software, H.L., S.K. and H.B.; validation, H.B. and J.U.; visualization, H.L. and J.H.; writing—original draft, H.L.; writing—review and editing, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (Innovative Incubation Center for Autonomous xEV Technology) (No. 5199990814084). This work was also supported by the Korea Evaluation Institute of Industrial Technology (KEIT) (No. 20011866).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interests.

List of Symbols

C t h the thermal capacitance
C t h _ j N the thermal capacitance between junction and NTC
E o n the energy of turn-on losses per switching cycle
E o f f the energy of turn-off losses per switching cycle
E r e c the energy of reverse-recovery losses per switching cycle
f e the frequency of electric angle
f s w the switching frequency
i d e * the command of d-axis current
I c the collector current
I f the forward current
i q e * the command of q-axis current
P c o n d [ D i o d e ] the conduction losses of the diode
P c o n d I G B T the conduction losses of the IGBT
P D i o d e the power losses of the diode
P I G B T the power losses of the IGBT
P s w [ D i o d e ] the switching losses of the diode
P s w [ I G B T ] the switching losses of the IGBT
P t o t a l the total power losses of the IGBT module
P t o t a l [ D i o d e ] the total power losses of the diode
P t o t a l [ I G B T ] the total power losses of the IGBT
R t h the thermal resistance
R t h _ j N the thermal resistance between junction and NTC
T e s t the temperature of estimation
T j u n c t i o n the temperature of junction
T m e a s u r e the temperature of measurement by IR camera
T N T C the temperature of NTC
T r e f the torque reference
V c e s a t the collector-emitter saturation voltage
V d c the DC-link voltage
V d e * the voltage command of d-axis current controller
V f the forward voltage
V q e * the voltage command of q-axis current controller
Z t h _ j N the thermal impedance between junction and NTC
Indices
ccollector
condconduction
ded-axis of synchronously rotating reference frame
eelectric angle
estestimation
fforward bias
fdbfeedback
Invinverter
jjunction
jNjunction to NTC
nomnominal
qeq-axis of synchronously rotating reference frame
recreverse-recovery
refreference
Rggate resistor
satsaturation
swswitching
ththermal
Vgegate-emitter voltage

References

  1. Infineon Technologies. Using the NTC inside a Power Electronic Module; Application Note AN2009-10; Infineon Technologies: Neubiberg, Germany, 2009. [Google Scholar]
  2. Hu, Z.; Ge, X.; Xie, D.; Zhang, Y.; Yao, B.; Dai, J.; Yang, F. An Aging-Degree Evaluation Method for IGBT Bond Wire with Online Multivariate Monitoring. Energies 2019, 12, 3962. [Google Scholar] [CrossRef] [Green Version]
  3. Wu, R.; Wang, H.; Ma, K.; Chimire, P.; Iannuzzo, F.; Blaabjerg, F. A temperature-dependent thermal model of IGBT modules suitable for circuit-level simulations. In Proceedings of the IEEE Energy Conversion Congress and Exposition, Pittsburgh, PA, USA, 14–18 September 2014. [Google Scholar]
  4. Martin, S.; Ma, X. Correlating NTC-Reading and Chip-Temperature in Power Electronic Modules. In Proceedings of the PCIM Europe 2015, International Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management, Nuremberg, Germany, 19–21 May 2015. [Google Scholar]
  5. Seiichiro, I.; Mikio, I. J-Series IPM and T-PM for EV and HEV Applications; TECHNICAL REPORTS; Mitsubishi Electric ADVANCE: Tokyo, Japan, March 2015. [Google Scholar]
  6. Kim, Y.S.; Sul, S.K. On-Line Estimation of IGBT Junction Temperature Using On-State Voltage Drop. In Proceedings of the Conference Record of 1998 IEEE Industry Applications Conference. Thirty-Third IAS Annual Meeting, St. Louis, MO, USA, 12–15 October 1998; pp. 441–444. [Google Scholar]
  7. Mantooth, H.A.; Hefner, A.R., Jr. Electro-thermal simulation of an IGBT PWM inverter. IEEE Trans. Power Electron. 1997, 12, 474–484. [Google Scholar] [CrossRef]
  8. Ishiko, M.; Kondo, T.; Usui, M.; Tadano, H. A compact calculation method for dynamic electro-thermal behavior of IGBTs in PWM inverters. In Proceedings of the 2007 Power Conversion Conference-Nagoya, Nagoya, Japan, 2–5 April 2007; pp. 1043–1048. [Google Scholar]
  9. Du, B.; Hudgins, J.L.; Santi, E.; Bryant, A.T.; Palmer, P.R.; Mantooth, H.A. Transient electrothermal simulation of power semiconductor devices. IEEE Trans. Power Electron. 2010, 25, 237–248. [Google Scholar]
  10. Lakhsasi, A.; Hamri, Y.; Skorek, A. Partially coupled electro-thermal analysis for accurate prediction of switching devices. Proc. Can. Conf. Electr. Comput. Eng. 2001, 1, 375–380. [Google Scholar]
  11. Riccio, M.; Irace, A.; Breglio, G.; Spirito, P.; Napoli, E.; Mizuno, Y. Electro-thermal instability in multi-cellular trench-IGBTs in avalanche condition: Experiments and simulations. In Proceedings of the 2011 IEEE 23rd International Symposium on Power Semiconductor Devices and ICs (ISPSD), San Diego, CA, USA, 23–26 May 2011; pp. 124–127. [Google Scholar]
  12. Greco, G.; Vinci, G.; Bazzano, G.; Raciti, A.; Cristaldi, D. Layered electro-thermal model of high-end integrated power electronics modules with IGBTs. In Proceedings of the IECON 2014-40th Annual Conference of the IEEE Industrial Electronics Society (IECON), Dallas, TX, USA, 29 October–1 November 2014; pp. 1575–1580. [Google Scholar]
  13. Hefner, R.; Blackburn, D.L. Simulating the dynamic electrothermal behavior of power electronic circuits and systems. IEEE Trans. Power Electron. 1993, 8, 376–385. [Google Scholar] [CrossRef]
  14. Yun, C.S.; Malberti, P.; Ciappa, M.; Fichtner, W. Thermal component model for electrothermal analysis of IGBT module systems. IEEE Trans. Adv. Packag. 2001, 24, 401–406. [Google Scholar]
  15. Kojima, T.; Yamada, Y.; Nishibe, Y.; Torii, K. Novel RC compact thermal model of HV inverter module for electro-thermal coupling simulation. In Proceedings of the 2007 Power Conversion Conference-Nagoya (PCC), Nagoya, Japan, 2–5 April 2007; pp. 1025–1029. [Google Scholar]
  16. Castellazzi, A.; Ciappa, M.; Fichtner, W.; Batista, E.; Dienot, J.; Mermet-Guyennet, M. Electro-thermal model of a high-voltage IGBT module for realistic simulation of power converters. In Proceedings of the ESSDERC 2007-37th European Solid State Device Research Conference, Munich, Germany, 11–13 September 2007; pp. 155–158. [Google Scholar]
  17. Gragger, J.V.; Fenz, C.J.; Kernstock, H.; Kral, C. A fast inverter model for electro-thermal simulation. In Proceedings of the 2012 Twenty-Seventh Annual IEEE Applied Power Electronics Conference and Exposition (APEC), Orlando, FL, USA, 5–9 February 2012; pp. 2548–2555. [Google Scholar]
  18. Batard, C.; Ginot, N.; Antonios, J. Lumped dynamic electrothermal model of IGBT module of inverters. IEEE Trans. Compon. Packag. Manuf. Technol. 2015, 5, 355–364. [Google Scholar] [CrossRef]
  19. Qian, C.; Fan, J.; Tang, H.; Sun, B.; Ye, H.; Zhang, G. Thermal management on IGBT power electronic devices and modules. IEEE Access 2018, 6, 12868–12884. [Google Scholar] [CrossRef]
  20. Infineon Technologies. Transient Thermal Measurements and Thermal Equivalent Circuit Models; Application Note AN2015-10; Infineon Technologies: Neubiberg, Germany, 1995. [Google Scholar]
  21. Arash, N.; Ashkan, N.; Osama, A.M. A Physics-Based, Dynamic Electro-Thermal Model of Silicon Carbide Power IGBT Devices. In Proceedings of the 2013 Twenty-Eighth Annual IEEE Applied Power Electronics Conference and Exposition (APEC), Long Beach, CA, USA, 17–21 March 2014; pp. 513–518. [Google Scholar]
  22. Infineon Technologies. Calculation of Major IGBT Operating Parameters; Application Note ANIP9931E; Infineon Technologies: Neubiberg, Germany, 1999. [Google Scholar]
Figure 1. Structures of FF1000R17IE4. (a) Circuit diagram of FF1000R17IE4. (b) Six-parallel package structure of FF1000R17IE4.
Figure 1. Structures of FF1000R17IE4. (a) Circuit diagram of FF1000R17IE4. (b) Six-parallel package structure of FF1000R17IE4.
Sensors 21 02454 g001
Figure 2. Block diagram of the power loss calculation algorithm.
Figure 2. Block diagram of the power loss calculation algorithm.
Sensors 21 02454 g002
Figure 3. The result of optimization with power loss parameters of an insulated gate bipolar transistor (IGBT). (a) The optimized collector-emitter saturation voltage as a function of collector current at each different junction temperature. (b) The optimized turn-on switching losses as a function of collector current at each different junction temperature. (c) The optimized turn-off switching losses as a function of collector current at each different junction temperature.
Figure 3. The result of optimization with power loss parameters of an insulated gate bipolar transistor (IGBT). (a) The optimized collector-emitter saturation voltage as a function of collector current at each different junction temperature. (b) The optimized turn-on switching losses as a function of collector current at each different junction temperature. (c) The optimized turn-off switching losses as a function of collector current at each different junction temperature.
Sensors 21 02454 g003
Figure 4. The result of optimization with power loss parameters of the diode. (a) The optimized forward voltage as a function of the forward current at each different junction temperature. (b) The optimized switching losses as a function of the forward current at each different junction temperature.
Figure 4. The result of optimization with power loss parameters of the diode. (a) The optimized forward voltage as a function of the forward current at each different junction temperature. (b) The optimized switching losses as a function of the forward current at each different junction temperature.
Sensors 21 02454 g004
Figure 5. The proposed thermal impedance model from the junction to the NTC. (a) The proposed thermal circuit model shown in a schematic cross section of an IGBT module including the heatsink and water jacket. (b) The proposed thermal equivalent circuit diagram.
Figure 5. The proposed thermal impedance model from the junction to the NTC. (a) The proposed thermal circuit model shown in a schematic cross section of an IGBT module including the heatsink and water jacket. (b) The proposed thermal equivalent circuit diagram.
Sensors 21 02454 g005
Figure 6. Black paint module of the FF1000R17IE4.
Figure 6. Black paint module of the FF1000R17IE4.
Sensors 21 02454 g006
Figure 7. Simulink model of the proposed the thermal impedance model.
Figure 7. Simulink model of the proposed the thermal impedance model.
Sensors 21 02454 g007
Figure 8. Experimental setup to verify the proposed algorithm. (a) Inverter with black paint module assembled. (b) Inverter installed in a dynamo environment, IR camera, and temperature recorder.
Figure 8. Experimental setup to verify the proposed algorithm. (a) Inverter with black paint module assembled. (b) Inverter installed in a dynamo environment, IR camera, and temperature recorder.
Sensors 21 02454 g008
Figure 9. Block diagram of the proposed junction temperature estimation algorithm.
Figure 9. Block diagram of the proposed junction temperature estimation algorithm.
Sensors 21 02454 g009
Figure 10. Comparison of estimation results in the stall condition (Case 1). (a) The estimated junction temperature of the IGBT (red), the measured junction temperature of the IGBT (blue), and the measured temperature of the NTC (green). (b) The estimated junction temperature of the diode (red), the measured junction temperature of the diode (blue), and the measured temperature of the NTC (green).
Figure 10. Comparison of estimation results in the stall condition (Case 1). (a) The estimated junction temperature of the IGBT (red), the measured junction temperature of the IGBT (blue), and the measured temperature of the NTC (green). (b) The estimated junction temperature of the diode (red), the measured junction temperature of the diode (blue), and the measured temperature of the NTC (green).
Sensors 21 02454 g010
Figure 11. Comparison of estimation results in the driving condition (Case 1). (a) The estimated junction temperature of the IGBT (red), the measured junction temperature of the IGBT (blue), and the measured temperature of the NTC (green). (b) The estimated junction temperature of the diode (red), the measured junction temperature of the diode (blue), and the measured temperature of the NTC (green).
Figure 11. Comparison of estimation results in the driving condition (Case 1). (a) The estimated junction temperature of the IGBT (red), the measured junction temperature of the IGBT (blue), and the measured temperature of the NTC (green). (b) The estimated junction temperature of the diode (red), the measured junction temperature of the diode (blue), and the measured temperature of the NTC (green).
Sensors 21 02454 g011
Table 1. The simulation results of the power loss calculation model.
Table 1. The simulation results of the power loss calculation model.
Condition of ExperimentResults of
Power Loss Calculation
CaseCoolant
(°C)
fsw
(kHz)
fe
(Hz)
V d c
V
T r e f
%
MeasurementSimulationError
Total
K w
Total
K w
IGBT
K w
Diode
K w
%
16540550661.621.6261.140.4860.38
26580550301.641.6281.180.4480.73
36540680581.531.5391.0980.4410.59
46540810501.7291.7471.2610.4861.04
525405501002.432.4391.7180.7210.37
Table 2. R t h and C t h from the junction to the NTC extracted by experiments in the stall condition and the driving condition.
Table 2. R t h and C t h from the junction to the NTC extracted by experiments in the stall condition and the driving condition.
ResultIGBTDiode
R t h j N _ I G B T
(°C/kW)
C t h j N _ I G B T
(W/°C)
R t h j N _ D i o d e
(°C/kW)
C t h j N _ D i o d e
(W/°C)
Stall condition41.1311.21102.13.36
Driving condition54.7613.56155.595.16
Table 3. The simulation results of the junction temperature estimation with the proposed thermal impedance model in various conditions of the inverter.
Table 3. The simulation results of the junction temperature estimation with the proposed thermal impedance model in various conditions of the inverter.
Condition of Inverter DrivingMeasurement vs. Simulation Result
CaseCoolant
(°C)
f s w
(kHz)
f e
(Hz)
V d c
(V)
T r e f
(%)
IGBTDiode
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
16540550100153.8158.91.1155.5158.70.25
26540550100153157.21.5155.3158.50.8
36540550100155.4160.21.3157.4160.70.3
46540550100157161.81.2158.9162.20.25
56540550100155159.11.0156.6159.90.5
665833.3550100138.2135.30.94139.5136.71.03
765833.3680100146.1139.40.54146.91411.14
865833.3810100146.1136.40.75147.2138.21.3
Table 4. The experiment results of each case for verification of the proposed junction estimation algorithm in the stall condition.
Table 4. The experiment results of each case for verification of the proposed junction estimation algorithm in the stall condition.
Condition of Experiment CaseMeasurement vs. Estimation
CaseCoolant
(°C)
f s w
(kHz)
f e
(Hz)
V d c
(V)
T r e f
(%)
IGBTDiode
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
16540550100153.7152.30.91157156.70.19
2654068080144.9146.30.97145.5145.60.07
3654071080149.5150.10.4149.2148.40.54
4654081060134.21372.0132.9132.30.45
56540550100154.4153.60.52158.1157.90.13
66540550100154.5152.41.36157.9156.80.7
Table 5. The specifications of the traction motor-inverter system.
Table 5. The specifications of the traction motor-inverter system.
SpecificationsValue
DC-link voltageMinimum550 (V)
Nominal710 (V)
Maximum810 (V)
Motor speedRange0 to 12,000 (rpm)
TorqueMaximum500 (Nm)
PowerMaximum120 (kW)
Phase currentMaximum636 (A)
Switching frequency0 to 500 (rpm)4 (kHz)
500 to 12,000 (rpm)8 (kHz)
Table 6. The experimental results of each case for verification about the proposed junction estimation algorithm in the driving condition.
Table 6. The experimental results of each case for verification about the proposed junction estimation algorithm in the driving condition.
Condition of Experiment CaseExperiment vs. Estimation
CaseCoolant
(°C)
f s w
(kHz)
f e
(Hz)
V d c
(V)
T r e f
(%)
IGBTDiode
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
T m e a s u r e
(°C)
T e s t
(°C)
Error
(%)
1658TN550100135.61435.46124.3126.11.45
2658TN68080134138.53.36124.6125.60.8
3658TN81080151.5155.62.71144.3146. 61.59
465833.355080109.9113.23.0109.2109.50.27
565833.368080122.9124.91.63114.61182.97
665820.0710809594.90.7295.897.92.19
765833.371080126.7127.50.63117.5119.41.62
865833.371080144.2147.40.84130.9133.43.16
965833.381080139.5138.52.22126.5126.51.91
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lim, H.; Hwang, J.; Kwon, S.; Baek, H.; Uhm, J.; Lee, G. A Study on Real Time IGBT Junction Temperature Estimation Using the NTC and Calculation of Power Losses in the Automotive Inverter System. Sensors 2021, 21, 2454. https://doi.org/10.3390/s21072454

AMA Style

Lim H, Hwang J, Kwon S, Baek H, Uhm J, Lee G. A Study on Real Time IGBT Junction Temperature Estimation Using the NTC and Calculation of Power Losses in the Automotive Inverter System. Sensors. 2021; 21(7):2454. https://doi.org/10.3390/s21072454

Chicago/Turabian Style

Lim, Heesun, Jaeyeob Hwang, Soonho Kwon, Hyunjun Baek, Juneik Uhm, and GeunHo Lee. 2021. "A Study on Real Time IGBT Junction Temperature Estimation Using the NTC and Calculation of Power Losses in the Automotive Inverter System" Sensors 21, no. 7: 2454. https://doi.org/10.3390/s21072454

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop