Junction Temperature Cycle Recording of an IGBT Power Module in a Hybrid Car

The accuracy of the lifetime calculation approach of IGBT power modules used in hybrid-electric powertrains suffers greatly from the inaccurate knowledge of application typical load-profiles. To verify the theoretical load-profiles with data from the field this paper presents a concept to record all junction temperature cycles of an IGBTpowermodule during its operation in a test vehicle. For this purpose the IGBT junction temperature ismeasuredwith amodified gate driver that determines the temperature sensitive IGBT internal gate resistor by superimposing the negative gate voltage with a high-frequency identification signal. An integrated control unit manages the TJ measurement during the regular switching operation, the exchange of data with the system controller, and the automatic calibration of the sensor system. To calculate and store temperature cycles on a microcontroller an online Rainflow counting algorithmwas developed.The special feature of this algorithm is a very accurate extraction of lifetime relevant information with a significantly reduced calculation and storage effort. Until now the recording concept could be realized and tested within a laboratory voltage source inverter. Currently the IGBT driver with integrated junction temperature measurement and the online cycle recording algorithm is integrated in the voltage source inverter of first test vehicles. Such research will provide representative load-profiles to verify and optimize the theoretical load-profiles used in today’s lifetime calculation.


Introduction
The combination of an internal combustion engine and an electric machine enables the improvement of the efficiency and the performance of the drivetrain of personal cars, busses, and utility vehicles [1].In view of the reliability and the lifetime of voltage source inverters used in hybrid-electric powertrains the IGBT power module can be considered as the most lifetime critical component.This is especially true if power modules with conventional linking and packaging technology are used.Those modules are characterized by a bond-wire connection, a direct copper bonded Al 2 O 3substrate where the chip is soldered on, and a copper baseplate.This results in a complex structure whose materials have different coefficients of thermal expansion CTE.In case of temperature cycles this CTE mismatch causes thermomechanical stresses in the modules interconnections and leads to the lift-off or the heel-cracking of bond wires or the degradation of the die-attach or the substrate solder joint [2].
To estimate the lifetime of an IGBT power module in a hybrid car a simple lifetime calculation approach has become dominant in recent years [3].This calculation approach is derived from the lifetime estimation of mechanical parts and demands the linkage of an application typical load-profile with a lifetime model of the IGBT power module using a cycle counting algorithm and a linear damage accumulation rule.For mechanical parts like shafts or gearwheels in transmissions this lifetime calculation approach could be verified over the years and today it is possible to design their lifetime with a high accuracy.On the contrary the lifetime calculation of IGBT power modules in hybrid cars is in a very early stage and currently it is not possible to quantify the accuracy of the lifetime calculation approach.What is known, however, is that the calculation approach suffers from different factors of uncertainty like the interaction of different failure mechanisms [4] and the information loss due to the cycle counting.However, the most critical point in lifetime calculation is the limited representation accuracy of

State of the Art
The lifetime calculation of an IGBT power module in a hybrid car requires basically the linkage of an application typical load-profile with an empirical lifetime model of the power module using a cycle counting algorithm.Today this load-profile originates from an application typical velocity profile like the New European Driving Cycle (NEDC) and the simulation of the entire hybrid-electric powertrain.A subsequent electrothermal model of the power module delivers the transient IGBT junction temperature over a certain timespan.This transient temperature curve is called the loadprofile of the power module.It is analyzed with a cycle counting algorithm and valued with an empirical lifetime model.In recent years the following lifetime models and counting algorithms have been presented.

Empirical Lifetime Models.
Empirical lifetime models originate from the accelerated ageing of power modules.They specify the number of temperature cycles a power module can bear until a failure criterion is reached.In recent years various lifetime models were publicized that differ primarily in the number of parameters used to describe a temperature cycle.The elementary lifetime model is a simple Coffin-Manson law [8] that states that the number of temperature cycles to failure   depends solely on the size of the amplitude Δ  of a temperature cycle: Today there are extended lifetime models [9][10][11] which consider additional parameters to describe a temperature cycle.
In 1997 the LESIT [9] project investigated the temperature cycle stability of power modules with conventional packaging technology from European and Japanese suppliers.It was found that the medium cycle temperature  ,med has a notable influence on the sustainable number of cycles.For this reason the Coffin-Manson law was extended by an Arrhenius term.Equation (2) shows the LESIT model, where the number of cycles to failure   is a function of the cycle amplitude Δ  and the medium cycle temperature  ,med .The parameters  = 640 and  = −5 were derived from accelerated ageing and  = 8.314 J/mol⋅K is the gas constant and  = 7.8 ⋅ 10 4 J⋅mol −1 is the activation energy: ) . ( Since the technologies of conventional IGBT power modules have been improved, in 2008 the number of sustainable temperature cycles to failure was reinvestigated by power cycling of several Infineon IGBT modules.It became apparent that many additional parameters have an impact on the module lifetime.The developed CIPS08 [10] lifetime model describes the number of cycles to failure   as a function of the amplitude Δ  , the minimum temperature  ,min =  ,med − (1/2) ⋅ Δ  , and the heating time  on of a temperature cycle: Moreover the current per bond wire   , the nominal voltage   , and the bond wire diameter  were taken into account.
The parameters  and  1 to  6 and their validity ranges are given in [12].For instance the heating time  on of a temperature cycle must be set to  on = 15 s for  on > 15 s.Equation (3) shows the CIPS08 lifetime model.For lifetime calculation the current per bond wire can be set to   = 10 A. The diameter of the bond wire  and the voltage class of the power module   are constants, so that the lifetime of the power module depends solely on the temperature cycles the power module is exposed to during its operation.This comparison of different lifetime models shows that the accuracy of the empirical models used for lifetime calculation could be improved due to the more accurate specification of a temperature cycle.The present state of the art is the parameterization of a temperature cycle with its amplitude Δ  , its minimum temperature  ,min , and its heating time  on .

Cycle Counting Algorithm.
Counting algorithms enable the evaluation of an application typical load-profile that consists of several different temperature cycles with an empirical lifetime model.For this purpose they extract and parameterize all temperature cycles within the load-profile and store them in a data vector.Widely accepted counting methods are the half-cycle counting, the maximum-edge counting, and the Rainflow counting [13,14].Figure 1 shows the application of these counting algorithms on an exemplary temperature profile.
In the half-cycle counting all rising and falling edges were counted as half temperature cycles.Their amplitudes are calculated as differences of two consecutive extreme values.The minimum temperature of each cycle equals the smallest cycle temperature.This applies to all mentioned algorithms.Half-cycles with rising edges consist of a heating time, which is the time span between their occurrences.The cooling time of half-cycles with falling edges is not considered in the reviewed lifetime models, so that these cycles do   The unique feature of the Rainflow method is the counting of closed temperature cycles, which appear as closed hysteresis curves in the stress-strain diagram [15].This physical background distinguishes the Rainflow method from the previous described algorithms.For this reason the Rainflow algorithm became the state of the art in the lifetime calculation of power modules.The principle of the Rainflow method is standardized by ASTM E-1049 [16].
To calculate closed temperature cycles with a Rainflow algorithm the extreme values of the entire load-profile must be known.Today there are different versions of Rainflow algorithms, like fast Rainflow, range-counting, or graphical algorithms [17].In the following the working principle of a simple range-counting Rainflow should be outlined.Figure 2 shows an exemplary temperature profile, whose extreme values are stored in a one-dimensional vector that is scanned by a moving window.A closed temperature cycle exists when the absolute difference of the extreme values 2 and 3 is greater or equals the absolute difference of the extreme values 1 and 2.Otherwise, the window has to be shifted.At the starting position the extreme values in the window do not meet the cycle condition so that the window is shifted by one extreme value.The new values in the window fulfill the cycle condition and the temperature cycle Δ 1 = 1 − 2 = 6 ∘ C is calculated.Finally the extreme values 1 and 2 were deleted and the scanning window is shifted back to the beginning of the vector.
Figure 3 shows the complete algorithm to convert a transient temperature profile into temperature cycles that are parameterized with their amplitudes Δ  , their minimum temperatures  ,min , and their heating times  on .Firstly a simple three-point algorithm determines all extreme values of the load-profile and stores them with a timestamp in a two-column vector.Once all extreme values of the loadprofile were determined, the vector is scanned by the rangecounting Rainflow algorithm in search of full temperature cycles.Whenever a closed temperature cycle is found it is parameterized and stored in a three-column data vector.
On the basis of the stored temperature cycles the lifetime  of the power module can be estimated according to (4), where   is the number of sustainable temperature cycles with a certain amplitude Δ  , minimum temperature  ,min , and heating time  on and   is the number of cycles with similar parameters that were found in the load-profile.In the case of a vector storage each temperature cycle has to be considered with   = 1.The lifetime  of the power module results from the linear accumulation [17] of the damage of all temperature cycles   and the scaling with the duration  LP of the application typical load-profile.Consider  = (

Temperature Cycle Recorder
The objective of this work is the verification of the state of the art to calculate the lifetime of an IGBT power module used in hybrid-electric powertrains.Therefore the temperature cycles of the power module should be determined during the operation in the field and stored on the system controller of the voltage source inverter.Basically this requires the measurement of the IGBT junction temperature during the regular inverter operation and the online temperature cycle calculation and storage on the system controller or the IGBT driver.Figure 4 shows the block diagram of the recording concept.The junction temperature is measured with a modified IGBT driver and transmitted to the system controller with an optical fiber.The online cycle recording algorithm calculates closed temperature cycles and stores them in a compact data matrix.In the following the IGBT driver with integrated junction temperature measurement and the online cycle recording algorithm are presented in detail.
3.1.IGBT Driver with   Measurement.Due to experiencing curve effects voltage source inverters used within hybrid cars contain standard IGBT power modules with conventional linking and packaging technology.The measurement of the IGBT junction temperature of those conventional power modules during their operation in a way that is suitable for series production is a very challenging task.Today's sensor concepts are based on either the installation of temperature sensors on the chip surfaces [18] or the utilization of temperature sensitive electrical parameters of the IGBT.In the laboratory the most popular temperature sensitive parameter is the saturation voltage   = (  ,   ) that is measured at a constant collector current of   = 100 mA.In this case the saturation voltage has a temperature sensitivity of approximately 2.3 mV/K [19].To realize the   -method during inverter operation a switching sequence with attached measurement phase could be used [20,21], but this also results in an interruption of the motor current and deteriorates the system properties.In [22,23] the temperature sensitive quasi-threshold voltage is used for junction temperature measurement.Thereby the induced voltage across the parasitic inductance between the power emitter and an auxiliary emitter is used to trigger the measurement of the gate emitter voltage.However, the method requires an increased measurement and calibration effort and can be considered to be very noise sensitive.The feasibility of other temperature sensitive parameters for online junction temperature measurement, such as switching times [24], the current slope [25], or the width of the miller plateau [26], is rather low.Reasons for this are their small temperature sensitivity in the range of about 1 ns/K, their poor selectivity, the limited resolution of affordable sensors, and the need to add subsequent compensating procedures.
A further temperature sensitive parameter that can be found in conventional power modules, where each IGBT is built up of several paralleled IGBT single chips, is the internal gate resistor  Gi .Often this resistor is located directly under the gate bond in the center of each chip and therefore closes to the junction of the semiconductors.To use this resistor as temperature sensor in [27] a modified IGBT single chip with double-sided bond connection of the internal gate resistor is presented.This allows the impression of a small DC sense current and the determination of the temperature sensitive internal gate resistance by measuring the voltage drop across the  Gi .Although the internal gate resistor was found to be well suited for junction temperature measurement the sensor system is inappropriate for paralleled chips and nonmodified, conventional IGBT power modules.
In summary, today there is no sensor concept suitable for series production to measure the junction temperature of a conventional power module during the real inverter operation.Existing solutions either are based on the modification  of the power module or require a large intervention in the power electronic system.In the following an IGBT driver concept with integrated junction temperature measurement [5] is presented that is easy to apply, does not affect the real inverter operation, and is applicable to conventional IGBT power modules.

Measurement Method.
To realize an easy to apply sensor system it is advantageous to integrate the junction temperature measurement in the gate circuit or in the IGBT driver, respectively.It can be seen in Figure 5 that the temperature sensitive internal gate resistor  Gi is in series with the parasitic capacitance   of the IGBT, the inductance   of the gate connection, and the external gate resistor   .To determine the resistance of  Gi utilizing the available terminal elements of a conventional power module the negative gate voltage is superimposed with the sinusoidal identification signal  GAC .
During the feeding in period the voltage drop  RG across the external gate resistor   is measured.To ensure a high measurement accuracy an external gate resistor with a low temperature dependency has to be used.Figure 6 sketches a short section of the gate voltage, where the negative voltage  Goff is superimposed with the sinusoidal, high-frequency identification signal  GAC ().To realize the measurement during the switching operation of the IGBT with a switching frequency of  sw = 5 kHz the sensor response time has to be smaller than  off = 100 s.The frequency of the identification signal is in the range of The measurement of the amplitude ûRG and the phase angle  of the voltage drop  RG enables the calculation of the internal gate resistor according to (5).If the identification frequency   equals the resonance frequency   of the gate circuit the phase angle becomes  = 0 ∘ .In this case the cos() = 1 is constant and solely the amplitude ûRG of the voltage drop across the external gate resistor   contains the temperature information and has to be measured.The measurement of the voltage amplitude ûRG can be done with a very simple sensor system.Consider Within a feasibility study the measurability of the internal gate resistor  Gi in the series connection was investigated.
Thereby it was assumed that the internal gate resistor can be measured if it has a sufficient large share in the total impedance  = || of the series connection.Figure 7 illustrates the external and the internal gate resistors   and  Gi as well as the reactances   and   of the parasitic components of the gate circuit in a vector diagram.
If the frequency   =   /(2) of the identification signal equals the resonance frequency   of the gate circuit both reactances   =   distinguish themselves and the ratio of  Gi / reaches its maximum value, where  =   +  Gi .The power module [28,29] used in further investigations has an  Gi / ratio of 62% at a resonance frequency of about   = 25.1 MHz.In the following a feeding-in method is presented that consists of a parallel path whose parasitic inductance   decreases the resonance frequency to   = 4.5 MHz.The driver circuit inductance   and therefore the switching behavior of the IGBT remain uninfluenced.

Feeding-in of the Identification Signal.
The superimposition of the negative gate voltage with a high-frequency identification signal is challenging.It is especially important not to alter the switching behavior of the IGBT, to consider safety aspects, and to ensure a continuous low-resistant connection of the gate to avoid an unwanted turn-on of the IGBT.In this context the serial feeding-in of a high-frequency signal using a transformer, whose secondary coil is integrated in series to the driver circuit, is an unsatisfactory solution [30].On the contrary Figure 8 illustrates the parallel feeding-in of a negative gate voltage  Goff that is superimposed with the sinusoidal identification signal  GAC .For this purpose a parallel auxiliary MOSFET is used.During the feedingin period this auxiliary MOSFET is turned on and the conventional IGBT driver becomes disabled so that its output is forced to a high impedance state.Due to the parallel feeding-in the original driver circuit and thus the IGBT switching behavior remain uninfluenced.To simplify the feeding-in process during the inverter operation the driver concept consists of an integrated control unit that receives the switching pattern and transmits the measured junction temperature via an optical fiber.Figure 9  Initially the IGBT driver is enabled EN = 1 and its highside-switch (HSS) is conductive (IN = 1), so that the positive voltage  Gon is applied at the gate and the IGBT is in on-state.At the time  0 the IGBT is switched off through the negative gate voltage  Goff that is applied by the low-side-switch (LSS) by setting IN = 0. To apply the negative, superimposed gate voltage  GAC at the time  1 the IGBT driver is disabled, setting EN = 0, so that its output is forced to a high impedance state.At the same time the feeding-in of the identification voltage is activated by setting HF = 1.To sustain a lowresistant gate connection the signal voltage source  GAC incorporates an impedance converter with a very low output impedance.During the feeding-in of the identification signal the measurement setup is analyzing the voltage drop  RG to generate an ADC compatible sensor voltage.At the time  2 the driver concept received the instruction to reactivate the IGBT.Therefore the identification voltage is disabled (HF = 0) and the positive gate voltage is applied, setting EN = 1 and IN = 1.In the case of failure the current limited buffer becomes deactivated, so that the gate signal of the IGBT driver is dominant.

Generation of the Sensor Output Voltage.
As it was outlined earlier, in case of resonance the amplitude of the voltage drop  RG across the external gate resistor   depends linearly on the junction temperature of the IGBT.In Figure 10 a measurement setup is presented that converts the amplitude change into an ADC compatible direct voltage.Therefore the voltage drop across   is picked up differentially with the subtractor S1, rectified with the Schottky diode , and smoothed by the capacitor .The level of the resulting direct voltage changes upon a certain offset linearly with the junction temperature.To adjust the voltage change to the input voltage range of an ADC, as a first step, the offset voltage  0 is subtracted in a way that the ADC input voltage equals  ADC = 1 V at an IGBT junction temperature of   = 20 ∘ C. Later, this offset voltage is used for the calibration of the junction temperature measurement.
Secondly, an operational amplifier adjusts the ADC input voltage to be  ADC = 3V at   = 120 ∘ C.This results in U GAC a temperature sensitivity of 20 mV/ ∘ C. When using a 10bit analog to digital converter with a reference voltage of 5 V the junction temperature can theoretically be measured with a resolution of 0.24 ∘ C/bit.The response time of the measurement setup to generate a valid sensor voltage is primarily defined by the capacitor  and is set to   < 70 s.

Experimental Results.
To qualify the presented measuring method a prototype of the IGBT driver with junction temperature measurement was developed and integrated in the double pulse experiment with   = 600 V and   = 250 A. A tempering system enables the homogenous heating of the power module to different junction temperatures   .
Figure 11 shows the gate voltage   at a switching frequency of  sw = 5 kHz.After the IGBT is switched off the negative gate voltage  Goff is superimposed with the identification signal  GAC = 0.6 V with a frequency of   ≈ 4.5 MHz.
The sensor output voltage  ADC reaches, after a response time of   ≈ 70 s, a constant voltage level that is sampled by the analog to digital converter of the control unit during the timespan   ≈ 6 s.It can be seen that the sensor output voltage increases linearly with an increasing junction temperature.
To determine the calibration curve of the driver concept the sensor output voltage was measured at four different junction temperatures.At each temperature level the sensor output voltage was measured several times to determine the standard deviation and the remaining noise of the temperature measurement.The measurement result in Figure 12 shows that the sensor output voltage  ADC is a very linear function of the junction temperature   .The standard deviation of the voltage measurement was found to be  UADC ≤ 24 mV that corresponds to a remaining noise of the temperature measurement of less than ±1.0 K. To realize an easy to apply IGBT driver, this calibration curve is programmed on the control unit of the driver concept, so that the measured sensor output voltage can be translated immediately in the junction temperature of the IGBT and transmitted to the system controller of the voltage source inverter with an optical fiber.
Because of process variations of the internal gate resistor the IGBT driver has to be calibrated at a known junction temperature, for example, at room temperature.Thereby the offset voltage  0 is set in a way that the sensor output voltage at a homogeneous calibration temperature complies with the reference calibration curve on the control unit.The automatic calibration process is organized by the control unit and enables the pairing of the IGBT driver to a power module that is still installed in the voltage source inverter.
To qualify the IGBT driver during the regular inverter operation a test setup consisting of one inverter phase that is connected to a virtual electric machine was developed.All IGBTs were controlled by a dSpace system.A MATLAB Simulink model allows the emulation of different loadprofiles of a hybrid-electric powertrain.During the regular switching of the IGBTs the developed IGBT driver measures the junction temperature of one IGBT and reports it to the dSpace system, where it is displayed in real time.To verify the temperature measurement with an IR camera a coated power module was used.Figure 13 shows the motor phase current  AC and the junction temperature   that is measured with the developed IGBT driver  ,DR and an IR camera  ,IR .The timespan of the load-profile is one minute.It can be seen that the temperature variations measured with the IGBT driver are in a very good agreement with those obtained with the IR camera.For this reason the developed IGBT driver is suitable for the temperature cycle recording of a power module during the inverter operation in a hybrid car.
Figure 14 shows the temperature profile of the IGBT for a timespan of 45 minutes.It can be seen that the loadprofile consists of a majority of operational temperature cycles caused by the transient power flow over the power module.Those active temperature cycles depend on the driver command and the superior hybrid strategy and lead to the heating of the coolant and the cold plate temperature during the first 15 minutes of the driving.The cool-down in the middle of the load-profile represents a time period without electric power flow.Consequently the coolant and the power module temperature decline until the subsequent operation cycles will cause the reheating of the system.After about 42 minutes the mission ends and the entire system cools down.In view of the homogeneous temperature variation of the power module during the start and the cooldown, these cycles should be labeled as passive ones.
Because of the limited storage capacity on the system controller it is not possible to store the whole temperature profile over a representative long period of time.For this purpose, an online recording algorithm is presented that processes the measured junction temperature into a compact set of information that is relevant for the verification of the theoretical load-profiles and the lifetime calculation approach.In view of the state of the art the verification can be done on the basis of temperature cycles that were calculated with a Rainflow algorithm and stored with their amplitude, their minimum temperature, and their heating time.

Online Temperature Cycle Recording Algorithm.
The online calculation and storage of temperature cycles on the system controller require a modification of the state-of-theart algorithms used on desktop computers with regard to their online ability, their computing power, and the required storage capacity.Thereby the overall objective is to reach a compromise between a minimum calculation and storage effort and a maximum accuracy regarding lifetime relevant information.Figure 15 shows the block diagram of the developed online cycle recording algorithm that consists of an extreme value filter algorithm, an online Rainflow algorithm, and a classified storage of the determined temperature cycles in a four-dimensional frequency distribution.In the following the working principles of the algorithms are presented.

Extreme Value Filter
Algorithm.The calculation of temperature cycles using a cycle counting algorithm requires previous identification of extreme values.Therefore a threepoint algorithm [31] with a sliding average filter can be used.The filter averages the temperature profile and enables the filtering of noise.It is not suitable to filter out small temperature variations in the range of Δ  ≈ 5 ∘ C, as they were caused by the hybrid strategy.The three-point algorithm picks out three subsequent temperature values.If the second value is higher than the first and third value a maximum is identified.Because of the limited filter level this algorithm finds a lot of extreme values, which have to be further processed in the online Rainflow.To reduce the computing power a logical filter algorithm was developed that enables the filtering of these small temperature variations but still delivers the exact temperature of valid extreme values.The principle of the filtering is outlined in Figure 16.In the filter routine two variables   and   were increased with the temperature   ().Only after the temperature declines, a predefined filter temperature   , a valid extreme value is detected.It occurred at the time   and has the temperature value   .
In Figure 17 the filter algorithm is applied to a cutout of the presented load-profile of a hybrid car.The filter temperature is set to   = 5 ∘ C. It can be seen that the loadprofile consists of many small temperature variations, which does not contribute notably to the damage of the power module.These temperature variations are smaller than the specified filter temperature so that they are not identified as valid extreme values.In a sequence of small temperature variations the logical filter algorithm determines the lowest extreme value.In this view the developed extreme value filter meets all requirements and combines an adjustable filter level and a high data reduction but still delivers the accurate temperature of valid extreme value.
The filter temperature should be set in a way that it enables the highest possible data reduction without the rejection of lifetime relevant information.To parameterize the filter for its use in a hybrid car Figure 18 shows the number of identified Optimum filter temperature: more than 85 percent.The filtering of the small temperature variations that are caused by the hybrid strategy results in a deviation of the calculated lifetime of  = +0.002percent.In view of the major reduction of extreme values and therefore the required computing power, this minor calculation error is acceptable.Consequently the optimum filter temperature to analyze the presented load-profile of a hybrid car is set to Finally the result of the filtering in the lifetime calculation should be demonstrated.Therefore Figure 19 shows the number   of temperature cycles with similar amplitudes Δ  that were determined with a filter level of   = 0 ∘ C and   = 5 ∘ C. Additionally the lifetime consumption due to the temperature cycles of each class was calculated with the CIPS08 lifetime model and normalized to one operating hour LC/h.It is obvious that cycles with amplitudes of Δ  ≤ 5 ∘ C make up a majority of all cycles but do not contribute notably to the lifetime consumption of the power module, so that the accumulated lifetime consumption is nearly identical.
In summary the filter algorithm prevents the calculation of temperature cycles that are smaller than a filter temperature of   = 5 ∘ C, so that the computing power of the online Rainflow algorithm can be significantly reduced.In Figure 19 all greyed temperature cycles with Δ  ≤   , which make up  a share of about 85 per cent, will not be calculated anymore.Despite the filtering the exact temperature of valid extreme values is determined.

Online Cycle Counting Algorithm.
The ability of counting algorithms to calculate cycles during the inverter operation depends on their working principle.As it was illustrated in Figure 1 the half-cycle or the maximum-edge method calculates a temperature cycle on the basis of two or three subsequent extreme values.For this reason the online implementation of those simple counting methods would be rather simple.Online Rainflow algorithms are described in [32][33][34].The feature of the Rainflow method is the counting of closed temperature cycles.For this reason an online Rainflow always demands the intermediate storage of extreme values in a working storage, so that its online implementation on the system controller is complicated.

Substitutability of the Rainflow Method.
In the following investigations it should be analyzed whether the Rainflow method can be substituted with an easy to implement halfcycle or maximum-edge algorithm [6].For this purpose the capability of the three algorithms to convert the presented load-profile into a distribution of defined temperature cycles is investigated.Figure 20 shows the number of temperature cycles  with similar amplitudes Δ  that were extracted with a Rainflow algorithm.Beside the majority of operational temperature cycles with amplitudes of Δ  < 70 ∘ C two passive cycles with amplitudes of Δ  > 100 ∘ C were calculated.These cycles are caused by the homogenous heating of the entire power module during the start and the cool-down.All cycles were valued with the CIPS08 lifetime model and normalized to their lifetime consumption per operating hour LC/h.It turns out that the lifetime consumption of the two passive cycles is greater than the accumulated lifetime consumption of all operational cycles.In this view it is very important to extract and parameterize passive cycles with a high accuracy.Because of its physical background the results of the Rainflow  counting should be taken as reference for the evaluation of half-cycle and maximum-edge counting.
The half-cycle method extracts twice as many operational cycles.Their valuing as half temperature cycles leads to a reduced operational lifetime consumption compared to the Rainflow algorithm [6].This applies equally to the passive cycles.Additionally the passive cycle amplitudes were found significantly smaller than in the Rainflow counting, so that their lifetime consumption is drastically undervalued.The maximum-edge counting compensates the undervaluing of the passive cycles partially due to the interpretation of the half temperature cycles as full ones.Nevertheless the cycle amplitudes differ from the Rainflow counting and cause a smaller lifetime consumption.
In summary, simple counting methods calculate a significant lower lifetime consumption than the Rainflow method so that they predict a longer lifetime of the power module.In this view it is not possible to substitute the Rainflow method with a simple counting algorithm and there is a need to implement an online Rainflow.
Online Rainflow Algorithm.For the calculation of temperature cycles with a state-of-the-art Rainflow algorithm all extreme values of the load-profile have to be stored in a vector.In the case of longer test drives this also results in a large quantity of data that has to be stored temporarily in a working storage.Contrarily to the presented rangecounting Rainflow that scans the working storage with a window in search of closed temperature cycles Figure 21 shows an online version of this principle.Thereby the window is fixed at the beginning of a working storage with a variable size.Once a new extreme value appears the working storage is shifted and the new extreme value is written to the 3 position.The Rainflow cycle condition is checked and if it is true the temperature cycle is calculated and exported to the storage.The accompanied extreme values were deleted and the working storage is shifted for two positions.Since one new extreme value can close more temperature cycles the Rainflow cycle condition is checked again until it becomes false.To calculate the cycle heating times a corresponding  working storage for the time stamps of the extreme values is available.To enable online cycle recording for an unlimited driving time with a storage capacity of a few bytes only the timespans between the extreme values were stored.This enables a reconstruction of heating times of up to 255 seconds that is enough compared with the maximum specified heating time of the CIPS08 model.

Working Storage Overflow. A critical parameter of an online
Rainflow algorithm that should be implemented on a microcontroller is the size of the required working storage.In the following the number of extreme values in the working storage during the analysis of the load-profile should be examined.Therefore Table 1 outlines four exemplary temperature profiles   (), the number of extreme values   in the working storage, and the extracted temperature cycles with the amplitudes Δ  .It becomes clear that solely the temperature profile 3 causes an increase of the number of extreme values in the working storage above the overflow level of  ,max = 10.After the start temperature is reached all cycles were closed simultaneously and the working storage is emptied.All other temperature profiles enable the instantaneous calculation of cycles and do not increase the number of extreme values in the storage.
In a hybrid car temperature profile 3 typically arises during the voltage-and frequency-controlled starting procedure of the electric machine [35].Figure 22 sketches an exemplary junction temperature profile during the start-up of the electric machine in a hybrid car.It consists of many small temperature cycles, whose amplitude and frequency depend on the size and the fundamental frequency of the inverter output current.
The oscillations of the junction temperature during the start-up of the electric machine cause an increase of the number of extreme values   within the working storage.This leads to an overflow of the available working storage with a maximum size of 17 extreme values and a loss of information.To avoid an overflow the developed online Rainflow algorithm was expanded by an overflow routine, which is scanning the working storage in search of the starting point of profile 3 once the overflow level is reached.Subsequently all temperature cycles of the detected sequence were calculated and deleted from the working storage.For this reason the overflow routine allows the use of a small working storage and improves the implementation of the Rainflow counting on a microcontroller.

Temperature Cycle Storage.
In view of the state of the art the lifetime calculation on the basis of the presented load-profile demands the storage of each cycle with its amplitude Δ  , its minimum temperature  ,min , and its heating time  on [6].To store the calculated temperature cycles on the system controller with a minimum storage capacity a four-dimensional frequency distribution should be used [7].Figure 23 outlines the storage of   cycles with similar amplitudes Δ  and minimum temperatures  ,min .
To consider cycle specific heating times, each temperature class is linked with a histogram containing the number   of cycles with equal time stamps.In Figure 23 the class width  CW to group cycles with similar amplitudes is set to  CW = 5 ∘ C. To reach a compromise between the required storage capacity and the achievable accuracy in lifetime calculation, the resolutions of the temperature and the time axis were varied.
As an example of this variation Figure 24 diagrams the storage capacity  and the lifetime calculation error  in dependency of an increasing class width  CW to group cycles with similar amplitudes Δ  .It turns out that a class width of  CW = 5 ∘ C enables a reduction of the storage capacity to   = 51.6 kilobytes.A further increase of the class width would lead to an increased error, which is caused by the temperature cycles that were rounded up into the next higher temperature class and their nonlinear overvaluing in the CIPS08 lifetime model.In the case of a class width of  CW = 5 ∘ C the calculated lifetime is affected by less than  = −0.5 percent, so that it is a suitable compromise between a minimum storage capacity and a good storage accuracy.

Validation of the Online Cycle Recording Algorithm.
Finally the properties of the developted online cycle recording algorithm should be evaluated.The major advantage of a conventional algorithm, which was presented before, is that each temperature cycle becomes valued with its exact individual cycle parameters.For this reason a conventional algorithm enables the most accurate lifetime calculation, so that it can be used as reference.In Figure 25 the required computing power to determine temperature cycles in the load-profile of a hybrid car and the required storage capacity to store them for lifetime calculation are diagrammed for the conventional algorithm and the developed online cycle recording algorithm.
It turns out that the required computing power to calculate the temperature cycles within the presented load-profile could be reduced by 85 percent due to an extreme value filter algorithm.Moreover the storage of the determined cycles in a four-dimensional frequency distribution and the variation of the class widths enable a reduction of the required storage capacity by 99 percent compared to minimum class widths of one degree and one second.Despite these major reductions the calculated lifetime on the basis of the temperature cycles determined by the online algorithm deviates less than 0.8 percent from the lifetime calculated with the conventional extracted cycles.Because of its properties the developed algorithm can be implemented on the system controller of the voltage source inverter to record the temperature cycles of the IGBT power module during its operation [24].

Conclusions
The analysis of the state of the art to design the lifetime of IGBT power modules used in hybrid cars showed that there is a need to verify the theoretical load-profiles with data from the field.For this purpose a temperature cycle recorder was presented that records the temperature cycles of an IGBT power module during its operation in a voltage source inverter.For the measuring of the junction temperature during inverter operation a modified IGBT gate driver was presented.The driver concept determines the temperature sensitive IGBT internal gate resistor by superimposing the negative gate voltage with a high-frequency identification signal.For this purpose the conventional IGBT driver was extended by an auxiliary MOSFET and a control unit that manages the feeding-in process during the regular switching operation, the transmission of the measured junction temperature, and the automatic calibration of the sensor system.Finally, it is shown that the driver enables realtime junction temperature measurement during the regular inverter operation.
In order to reduce the quantity of data that has to be stored on the system controller an online cycle recording algorithm was developed that processes the measured junction temperature into a frequency distribution of temperature cycles.A review of empirical lifetime models showed that each temperature cycle has to be parameterized with its amplitude, its minimum temperature, and its heating time.The algorithm consists of a filter that identifies all extreme values that are relevant for the lifetime calculation.An online rangecounting Rainflow algorithm processes them into closed temperature cycles and stores them in a four-dimensional frequency distribution.In relation to the cycle counting algorithm used on desktop computers the computing power could be decreased by 85 percent and the required storage capacity could be decreased to 8 kilobytes.Despite these reductions the compliance of the online algorithm and the conventional desktop computer algorithm amounts to 99%.
Currently the IGBT driver with junction temperature measurement and the online temperature cycle recording algorithm are integrated in the voltage source inverter of first test vehicles.On the basis of the recorded load-histories it will be possible to verify and adjust the theoretical load-profiles used in today's lifetime calculation approach.For this reason the presented temperature cycle recorder makes an important contribution to improve the accuracy of the lifetime calculation of IGBT power modules and the reliability of future voltage source inverters.

Figure 2 :
Figure 2: Working principle of a conventional range-counting Rainflow algorithm that scans a one-dimensional extreme value vector.

Figure 3 :Figure 4 :
Figure 3: Conventional algorithm used on desktop computers to determine and store the temperature cycles that are incorporated into the load-profile.

Figure 5 :
Figure 5: AC equivalent gate circuit of an IGBT power module consisting of the internal gate resistor  Gi and the parasitic parameters   and   .

Figure 6 :
Figure 6: The negative gate voltage  Goff is superimposed with the high-frequency identification signal  GAC () during the IGBT is in off-state.

Figure 7 :
Figure 7: Vector diagram of the resistors   and  Gi and the reactances   and   of the parasitic components   and   at the frequency   .

Figure 8 :
Figure 8: IGBT driver concept with an auxiliary MOSFET to feed in the identification voltage  GAC .The junction temperature of the IGBT is calculated on the basis of the measured voltage drop  RG .

Figure 9 :
Figure 9: Signals to control the IGBT driver and to superimpose the negative gate voltage  Goff with the sinusoidal identification signal  GAC ().

Figure 10 :
Figure 10: Circuit to measure the temperature dependent amplitude of the voltage drop  RG and to convert it into an ADC compatible direct voltage.

Figure 11 :
Figure 11: Measurement result of the driver concept where the identification voltage  GAC leads to a temperature dependent sensor voltage  ADC that is sampled by an analog to digital converter during the time span   .

Figure 12 :
Figure 12: Calibration curve of the IGBT driver with junction temperature measurement.

Figure 13 :Figure 14 :
Figure 13: Junction temperature during inverter operation measured with the IGBT driver and an IR camera.

Figure 15 :Figure 16 :Figure 17 :
Figure 15: Online cycle counting algorithm to calculate and store temperature cycles with less effort during the operation of the IGBT power module.

Figure 18 :
Figure 18: The number   of extreme values that were found in the presented load-profile of a hybrid car decreases with an increasing filter temperature   .

Figure 19 :
Figure 19: Number   of temperature cycles Δ  determined with the filter temperatures   = 0 ∘ C and   = 5 ∘ C and a subsequent Rainflow algorithm.

Figure 20 :
Figure 20: Number of temperature cycles  with similar amplitudes Δ  that were determined with a Rainflow counting algorithm.The cycles were valued with the CIPS08 lifetime model; LC/h = (Δ  ,  ,min ,  on ).

Figure 21 :
Figure 21: Working principle of the online Rainflow algorithm that consists of a variable working storage for extreme values and their time information and a routine to check the Rainflow cycle condition (|2 − 3| ≤ |2 − 1|).

Figure 22 :
Figure 22: The junction temperature profile during the start-up of the electric machine complies with the exemplary temperature profile number 3.

Figure 23 :
Figure 23: Storage of temperature cycles with similar (a) amplitudes Δ  and minimum cycle temperatures  ,min in a matrix.(b) Each matrix element contains a histogram that gives information about the cycle heating times  on .

Figure 24 :
Figure 24: Required storage capacity  and the error  in the calculated module lifetime in dependency of an increasing class width  CW .

Figure 25 :
Figure 25: Comparison of the conventional algorithm used on desktop computers with the developed online cycle recording algorithm that should be implemented on the system controller of the voltage source inverter.

Table 1 :
Number of extreme values   in the working storage during the analysis of four exemplary temperature profiles with an online Rainflow.