Simulation Tool for IGBT modules

R. Schnell ABB Switzerland Ltd, Semiconductors, Fabrikstrasse 3, CH 5600 Lenzburg, Switzerland e-mail: [email protected] Abstract:
ABB has recently released a new simulation tool based on Microsoft Excel for dimensioning of two level voltage source IGBT converters. The tool allows the calculation of switching and conduction losses under the assumption of sinusoidal output currents. Additionally static and transient temperature rises in the module and an optionally specified heatsink can be calculated. The losses in the IGBT and companion freewheeling diode are iteratively adjusted to the temperature rises. Therefore the calculation results allow a quick selection of ABB IGBT modules for the planned application.

losses are given by the integration of the forward losses The calculation of the conduction and switching losses (VCE0, rCE including RCC'+REE') up to T0/2: T /2 of the IGBT and the freewheeling diode are done based 1 0 Pcond IGBT = VCE 0 + rCE ( sin(t )) 2 (t )dt on the application parameters. T0 0 As application parameters the DC-link voltage, output frequency, switching frequency, modulation index and with (t) as the function of pulse pattern with IGBT the power factor cos can be specified. turned-on=1 and IGBT turned-off=0. (t) can be If the maximum permissible output current has to be substituted by a function of modulation (m) and phase calculated the junction and case or heatsink temperature angle (). needs to be specified. In this case the permissible output With an infinite switching frequency we get the duty current can be read out of a diagram or directly cycle variation over time (PWM pattern). calculated by using the solver. If the output current is given the simulation tool 1 calculates with 10 iterations the losses and temperature (t ) = (1 + m sin(t + )) 2 rise of the module. Additionally it is possible to calculate the temperature rise in an optionally specified cooler. In addition special Inserting (t) into the formula and solving the integral overload conditions can be specified and the simulation we get the conduction losses: tool calculates the temperature excursion during overload in the heatsink and the module. Again the 2 1 Pcond = 1 (VCE0 + rCE ) + m cos (VCE0 + rCE 2 ) losses are iteratively adjusted to the junction 2 4 8 3 temperature.

Introduction:

Loss Calculation:
The simulation tool offers a relatively exact and fast method for loss calculation. The model data for the IGBT modules are based on their respective data-sheet values. The calculations are done with a linear approximation of the devices forward characteristics and with a polynomial function of the IGBT and diode switching energy. Additionally the forward and switching-energy characteristics are temperature dependent. The power dissipation calculation of the IGBT and diode are done by an average calculation for the conduction and switching losses over one period T0 of the output frequency. This method does not require adding up the switching energies [1]. IGBT loss calculation: Since the IGBT of one switch conducts only over one half period the conduction

The simulation tool restricts the modulation index to m12345675368395

36
3
3395
3 The switching losses are the sum of all turn-on and turnoff energies at the switching events. The measured turn-on and turn-off energies given in the data-sheet can be described as a polynomial function (Esw=f(I)):

E sw = Eon + Eoff = (a + b I + c I 2 )

Since the DC-link voltage can vary in different applications the dependence of the switching energy on the DC-voltage needs to be considered. Within certain limits this dependence can be assumed linear:

PDiode = Pcond + Psw

Thermal Calculation:
The loss calculation yields in average losses over one period of the output frequency. In fact the losses per switch only occurs during one half-period whereas during the other half period the antagonist switch is producing losses. Figure 1 shows in the upper graph a PWM chopped current as it appears in case of a two level VSI. In addition the resulting junction temperature is shown (calculated with THIMP V3.08). It is obvious that the junction temperature oscillates with the frequency of the output current. In the lower graph the corresponding losses Psmoothed (THIMP) are shown. As a comparison the calculated average losses (Pav) from simulation tool are given (dashed line). If the junction temperature is calculated with Tjav = Pav * Rth evidently the peak value of junction temperature exceeds the result of Tj calculated with the average losses Pav. Therefore the calculation of Tj with average losses yields in a to optimistic value. In order to minimize this inaccuracy the simulation tool calculates with Pmod which is two times the average losses Pav dissipated during one half period of the phase current. The resulting junction temperature Tjmod is shown in the lower graph of figure 1 and matches in terms of the peak and bottom values quite well with the real value of Tj.
800 125 120 600 115 400 Ic [A] Tj [C] 200 110 Tj 105 100 0.100 125 120 115 Tj 110 105 0.065 0.070 0.075 0.080 t [s] 0.085 0.090 0.095 100 0.100 Ic

E sw = (a + b + c 2 )

VDC Vnom

To calculate the switching losses the switching energies have to be added up.

Psw =

1 E sw ( ) T0 n

whereas n depends on the switching frequency. Therefore we can calculate the switching losses as a function of phase-current and switching frequency.

Psw =

V 1 f sw (a + b + c 2 ) DC Vnom

The total IGBT losses are the sum of the conductionand switching losses:

PIGBT = Pcond + Psw

Diode loss calculation: The diode losses can be nearly calculated the same way as the IGBT. Since the freewheeling diode conducts when the IGBT is turned-off the function of the pulse pattern has to be negated:

2 1 Pcond = 1 (VF0 + rT ) m cos (VT 0 + rT 2 ) 2 4 8 3

In case of the diode the turn-on energy can be neglected. Therefore only the recovery energy counts. The recovery energy given in the data-sheet diagram can be described as a polynomial function:

0 0.060

0.065

0.070

0.075

0.080 t [s]

0.085

0.090

0.095

1600 1400 1200 1000 800 600 400 200 P

Psmoothed [W] Pav [W] Pmod [W] Tjmod [C] Tjav [C]

E rec = (a + b I + c I 2 )
The recovery losses as a function of phase-current and switching frequency and VDC can be written as:

0 0.060

Prec =

V 1 f sw (a + b + c 2 ) DC Vnom

Figure 1: junction temperature as a function of fO

The maximum junction temperature Tvj max in function of the phase output current frequency fO can be calculated if the transient thermal resistance is known:

The total Diode losses are the sum of the conductionand switching losses:

Tvj max = 2 PAV Ri

i =1

1 e

2 2 f O i 1 f O i

1 e

+ Tref

With Tref as the reference heatsink temperature Th for shows the thermal equivalent block diagram for base-less modules or case temperature Tc for based modules with base-plate: modules.

Heatsink and Transient Calculations:

In addition to the module calculation the simulation tool allows the user to calculate the temperature rises and module losses of the IGBT-module together with a specified cooler. The simulation tool uses a simplified approach to calculate temperature rises and cross-talk effects. This has the advantage that the calculations can be done with the given values in the module- and heatsink datasheets. More accurate methods based on finite element calculation require specific knowledge about module and cooler construction and need lots of time for parameter extraction. On the other hand the simplified method of the simulation tool allows quick and accurate simulation without detailed knowledge on cooler and module construction.
Switch 1

Based Modules
Switch n

PIGBT ZthJC

PDiode ZthJC

PIGBT ZthJC

PDiode ZthJC

Reference: Tc RthCH

ZthHA

Figure 3: thermal equivalent block diagram for devices with baseplate

As in case of base-less modules the static calculation of the temperature rise in the heatsink can be calculate with the thermal resistance of the cooler and the number of dissipating heat sources mounted on it:

Thermal equivalent block diagram: TH = T A + ( PIGBT + PDiode ) nS RthHA The simplification done in the simulation tool is the assumption of common reference temperature point where the temperature is assumed to be homogenous Additionally the temperature rise in the interface TCH over the full area. In case of base-less devices this needs to be calculated. Since it is possible more than one module is mounted on the heatsink and that the heat reference is the heatsink temperature (figure 2): sources are distributed in several modules, as for example in a three phase inverter built with three halveBase-less Modules bridge modules mounted on a single cooler, it becomes necessary to scale the interface resistance accordingly:
Switch 1 Switch n PIGBT ZthJH PDiode ZthJH PIGBT ZthJH PDiode ZthJH

TCH = ( PIGBT + PDiode ) n S

Reference: Th ZthHA

RthCH nM

With nM as the number of modules mounted on a single heatsink. Figure 2: thermal equivalent block diagram for base-less devices In order to calculate the temperature dependent semiconductor losses the simulation tool adjusts the The static calculation of the temperature rise in the junction temperature and the corresponding losses in heatsink can be calculate with the thermal resistance of several iterations. Depending on the module type and the cooler and the number of dissipating heat sources the calculation (with / without heatsink) the calculation mounted on it: starts with the initial conditions Tvj = TA, Tvj = TH or Tvj = TC .

TH = T A + ( PIGBT + PDiode ) nS RthHA

With nS as the number of switches mounted on the cooler. For based devices the reference is the case temperature. The thermal characteristics of devices with a base-plate (copper or AlSiC) is normally specified with a ZthJC (thermal impedance junction case) for the IGBT and Diode part. Additionally the interface resistance caseheatsink RthCH is given for the whole module. Figure 3

Transient Overload Calculation: The simulation tool offers the possibility to additionally calculate the transient thermal behavior of the IGBT module and the heatsink. For the transient calculation overload conditions can be specified. In order to have full flexibility the start conditions of the temperature excursion calculation need to be specified. This includes the overload duration, the start junction temperature and the start case or heatsink temperature. Therefore it is possible to study the transient thermal behavior from any temperature level and even the

thermal recovery behavior from an overload can be The approach is to calculate the temperature ripple at the end of the overload pulse and add it to the calculated computed. The transient temperature rise can be calculated with the average temperature: thermal impedance.

Z th (t ) = Ri (1 e t / i )
i =1

Tvjmax(t) =Tvjav(t) + (2 P(T) Ri

i=1

1 e 1 e

0.5 fOi 1 fOi

+Tref ) P(T) Rth

The temperature rise T(t) is a function of Zth(t) and the temperature dependent power dissipation in the This is valid as long as the Zth(t) of the module is close to its Rth value. For most modules this is the case semiconductor P(T). between 0.5 and 1s. Below this time span the values for Tvj max(t) are slightly optimistic. Therefore the T simulation tool does not allow transient calculations T (t ) = ( start ) Z th (t ) + P(T ) Z th (t ) Rth with duration below 1s. In order to include the starting conditions for T the first term is introduced. Tstart/Rth describes the initial constant power that resulted in Tstart. This deposited power influences the thermal behavior until Zth(t)=Rth(t). The second term describes the heating with the temperature dependant power P(T). Since T is as well a function of the dissipated power and the thermal impedance iterations are necessary to get an accurate result. For the transient heatsink temperature we can write:

Usage:

The simulation tool consists of two excel sheets. In the output sheet (figure 4) the application specific parameters such as output current, DC-voltage, output frequency, switching frequency, modulation, load power factor and the temperatures can be specified. Additionally a device has to be selected in the dropdown menu. The selected device is confirmed in the blue field and the electrical and thermal models of the device are loaded. The calculation results are shown below and include conduction and switching losses as well as the maximum junction temperature and the output T T TH (t ) = TA + ( Hstart A ) Z thHA (t ) + P(T ) Z thHA (t )frequency dependent ripple of the junction temperature. RthHA With the button solve for Iomax the solver calculates the maximum permittable current if the case/heatsink and whereas P(T) depends on the number of switches junction temperatures are specified. This feature requires the properly installed solver plug-in. mounted on the heatsink
Simulation Tool for ABB IGBTs
V 4.00, Sept. 2002
1232456557895855

P(T ) = ( PIGBT (T ) + PDiode (T )) n s

For the transient average junction temperature we can write:

Select Device: Enter application parameters: Output current Io rms [A]: DC link voltage Vdc [V]: Output frequency f0 [Hz]: Switching frequency fsw [Hz]: Modulation index m: Load power factor cos: Th [C] Junction temperature Tvj [C]: (only for diagrams and solver) input 200 600 50 4000 1.00 0.85 80 125

LoPak5, SPT, 300A, 1200V

limits Io < 2 x Icn/sqrt(2) =424A 360V < Vdc < 840V 1Hz < f0 < 1000Hz 5 x f0 < fsw < 30000Hz 0<m<1 -1 < cos <1 5C < Th [C] < 125 5C < Tvj < 150C Gate resistor RG [Ohm]:

used values 200 600 50 4000 1.00 0.85 80 125 3.3

Solve for Io max Calculated results Io max [A]: Losses per IGBT & Diode [W]: IGBT total losses [W]: Tvj max IGBT [C]: Diode total losses [W]: Tvj max Diode [C]: 250.2 391.5 318.7 125 72.9 96

Tvj av (t) = TH (t) + (

Tvj start TH start RthJH

) ZthJH(t) + P(T) ZthJH(t)

Calculated Results: Total switch losses(IGBT&Diode) [W]: IGBT: IGBT total losses [W]: IGBT conduction losses [W]: IGBT switching losses [W]: Tvj max [C]: Ripple Tvj [C]:

281.3

Tvj max [C]: Diode: Diode total losses [W]: Diode conduction losses [W]: Diode switching losses [W]: Tvj max [C]: Ripple Tvj [C]:

112

In case of modules with a base plate TH has to be replaced with TC and RthJH and ZthJH have to be replaced with RthJC and ZthJC respectively. In addition the temperature drop across the interface TCH has to be calculated. Since the interface has no heatcapacity this can be done similar to the static calculations (page 3). Therefore we are able to calculate the transient average junction temperature. As already mentioned the average junction temperature yields in too optimistic value if the temperature ripple in function of the output current needs to be considered. To avoid a too high complexity the simulation tool uses a simplified approach that is valid in most relevant application cases.

228.0 156.3 71.7 112 5.1

53.2 26.2 27.0 91 2.6

Figure 3: Output sheet

Alternatively the maximum output current can be read out of the first diagram. The second diagram shows the temperature gradient and the ripple of the junction temperatures inside the module (figure 4).

Simulation Tool for ABB IGBTs

V 4.02, Mar. 2003

Average losses for sinusoidal output current at Tvj = 125C 5SNS 0300U120100
900
Total switch losses

Heatsink and Transient Calculations

Select Heatsink: Selected Device: Enter application parameters: Output current Io rms [A]: DC link voltage Vdc [V]: Output frequency f0 [Hz]: Switching frequency fsw [Hz]: Modulation index m: Load power factor cos: input 150 650 50 6000 1.00 0.85 limits Io < 2 x Icn/sqrt(2) =424A 360V < Vdc < 840V 1Hz < f0 < 1000Hz 5 x f0 < fsw < 30000Hz 0<m<1 -1 < cos <1 used values 150 650 50 6000 1.00 0.85 Heatsink Characteristics Rth [K/W] 8.61E-04 -1.73E-03 1.17E-01 -8.34E-02 RthHA [K/W] Enter ambient temperature and configuration of devices TA [C] 40 Number of switches on heatsink 6 Number of devices on heatsink 1 Calculated Results: Total switch losses(IGBT&Diode) [W]: 0C < TA < 100C 40 6 1 th [s] 1.05E+04 8.98E+02 6.74E+01 6.30E+01 0.033
112345678997
78

forced air cooling, 200x215mm, 33K/kW LoPak5, SPT, 300A, 1200V

800 700
Average Losses [W]

IGBT total losses Diode total losses max IGBT losses max Diode losses

600 500 400 300 200 100 0 0

248.1

Tvj max [C]: Diode: Diode total losses [W]: Diode conduction losses [W]: Diode switching losses [W]: Tvj max [C]: Ripple Tvj [C]:

116

50

100

150

200

250

300

350

400

450

Output Current Io rms [A]

Temperature gradient and ripple 2 x (Tvj max - Tvj av)

120 100 80
T [C] DTvj-h IGBT DTvj-h Diode ripple IGBT DTvj ripple Diode DTvj

IGBT: IGBT total losses [W]: IGBT conduction losses [W]: IGBT switching losses [W]: Tvj max [C]: Ripple Tvj [C]: Th [C]

190.6 102.6 88.0 116 4.3 89.4

57.5 17.7 39.8 102 2.8

Transient Overload Calculations

Enter overload conditions: Output current Io rms [A]: DC link voltage Vdc [V]: Output frequency f0 [Hz]: Switching frequency fsw [Hz]: Modulation index m: Load power factor cos: Tj start IGBT [C] * Tj start Diode [C] * Th [C] start* Overload duration [s] 200 650 5 6000 0.50 0.50 40 40 40 120 360V Io < 2 x Icn/sqrt(2) =424A < Vdc < 840V 1Hz < f0 < 1000Hz 5 x f0 < fsw < 30000Hz 0<m<1 -1 < cos <1 40 < Tj start IGBT [C] * < 150 40 < Tj start Diode [C] *< 150 40 < Th [C] < 40 1 < t [s] < 1000 TA [C] Tvj max [C]: Diode: Diode total losses [W]: Diode conduction losses [W]: Diode switching losses [W]: Tvj max [C]: 200 650 5 6000 0.50 0.50 40 40 40 120 40 137

60 40 20 0 0 50 100 150 200 250 300 350 400 450

Output Current Io rms [A]

* as default values select results from heatsink calculation results

Calculated Results after overload t =120 s Total switch losses(IGBT&Diode) [W]: 351.0 IGBT: IGBT total losses [W]: IGBT conduction losses [W]: IGBT switching losses [W]: Tvj max [C]: Th [C]

233.0 111.8 121.2 137 95.6

118.0 66.5 51.6 129

Figure 4: Diagrams with calculation results

Temperature Excursion During Overload 5SNS 0300U120100

The sheet heatsink and transient calculates the steady state and transient temperature excursions with the module and the heatsink (figure 5). The application specific parameters can be specified. The device to simulate is taken from the output sheet. In addition the heatsink, the ambient temperature and the number of switches and modules mounted on the cooler have to specified. In order to calculate the transient heatsink behavior it is possible to specify up to n=4 exponential terms for of Zth of the cooler. If the Zth of the cooler is not known the RthHA can be specified in the first term (Rth) and all other terms have to be set to zero (except the that will be automatically set to 1 if not specified). The heatsink parameters have to be entered into the sheet heatsink_spec. If the necessary values are specified the simulation tool calculates the junction temperature and the corresponding losses of the IGBT and the freewheeling diode. For the transient overload calculation a separate electrical load can be specified such as a possible overload or motor start conditions. To have full flexibility the starting temperatures for the Tvj of the IGBT and the Diode as well as the heatsink or case temperature can be entered. Therefore it is possible to simulate the transient behavior from any start condition and even thermal relaxation from overload can be computed. The simulation tool reports the IGBT and Diode losses as well as the temperatures that have been computed after the specified overload duration has been elapsed. Additionally the transient temperature excursion over time will be shown in a diagram (figure 6).

160 140 120 100

T [C]

80 60 40 20 0 0 20 40 60
t overload [s]

Tj IGBT [C] Tj Diode [C] Th [C]

80

100

120

140

Figure 6: heatsink and transient calculations

Conclusion:
The Simulation Tool presented offers a fast and simple but still accurate method to simulate the IGBT modules available from ABB under various application conditions. This even includes the possibility to calculate the transient thermal behavior together with a user-specified heatsink. Right now the thermal and electrical models of the LoPak line up are available, this includes as well models for parallel connected LoPaks (as dual) that already have a derating included. Therefore the Simulation Tool permits the user to select the right device for his application and furthermore gives already valuable application specific information on the behavior of the module in the planned application. Nevertheless the described simplifications done by the simulation tool have to be considered in the simulation result judgement.

References
1. D.Srajber, W.Lukasch, "The calculation of the power dissipation for the IGBT and the inverse diode in circuits with sinusoidal output voltage", electronica'92, Proc, pp. 51-58 M.Rahimo, A.Kopta, R.Schnell, G.Debled, J.de Oate, S.Linder, "New Industry Standard LoPak Modules with 1700V SPT-IGBT and Diode chip set", Proc. PCIM'02, Nrnberg 2002.

2.

3.

M.Rahimo, W.Lukasch, C. von Arx, A.Kopta, R.Schnell, S.Dewar, S.Linder Novel Soft-Punch-Through (SPT) 1700V IGBT Sets Benchmark on Technology Curve, Proc. PCIM01, Nrnberg, 2001.