Modeling Thyristor and Diodes; On-State Voltage and Transient Thermal Impedance, Effective Tools in Power Electronic Design

John W. Motto Jr., William H. Karstaedt, Jerry M. Sherbondy Sr., Scott G. Leslie

Abstract: Circuit modeling is an essential tool in the design of power electronic applications. Both the on-state forward voltage drop and transient thermal impedance of high power SCR’s and Diodes are complex functions requiring, in the past, tedious, copious, inefficient, semi-mechanical calculations to evaluate one or two specific conditions for a given power electronic application. Today, with the availability of personal computers, greater efficiency and accuracy can be realized using mathematical models of these basic characteristics. The widespread use of these models, depends upon; the understanding of how to use these models, the trade-off of the accuracy and simplicity of the models and the availability of parametric data for devices. This main objective of this paper is how to use the modeling equations to evaluate a given Power Electronic Application.

Previous papers by the authors1,2,3 and others6,8 have addressed this subject. As has been described, the forward on-state forward voltage drop can be modeled by both the classical ABCD and the new MNOP… parameters. The transient thermal impedance has been shown to be well represented via four or five exponential terms representing the significant transient thermal time constants of the device. The methodology to quickly and accurately calculate the forward voltage drop and transient thermal impedance modeling parameters has been discussed and described in detail1,2,3. This paper will describe the techniques for using the basic device modeling equations to evaluate a given power electronic application. Both the digital (hard code) simulation and "analog" PSPICE simulation techniques will be described. Also, a new tool, called STARSim, will be described. STARSim permits comparison of the older superposition method and new modeling method by performing the calculations using both approaches. A simple SCR model for PSPICE incorporating the on-state forward voltage drop and transient thermal impedance models is also described and evaluated. This information is considered to be important in providing the tools required for power electronic design engineers to quickly evaluate a given device in proposed or existing power electronic application.

I. INTRODUCTION


Circuit modeling is an important tool in the design of power electronic applications. Power distribution and transmission applications require explicit analysis to optimize cost of prototype and production units. Industrial power applications must minimize cost and time to market. The models chosen, however, must be accurate and closely aligned with the manufacturer’s specifications and rating system to be effective. This paper will describe the techniques in using the basic device modeling equations to evaluate a given power electronic application. Both the digital (hard code) simulation and "analog" PSPICE simulation techniques will be described. Also, a new tool, called STARSim, will be described. This program will permit the confidence which has been established with the older more tedious methods such as curve look up and superposition to be transferred to the new modeling methods by performing the calculations of both the old and new approaches, permitting comparison of the two methods. Tests will be made on two devices the Powerex (50mm) T9G0 and the (75mm) TBK7. This information is considered to be an important advancement in providing the tools required for power electronic design engineers to quickly evaluate a given device in proposed or existing power electronic applications.



II. THE CLASSICAL RATING METHOD (Curve Lookups & Superposition) and the NEW
Vtm, TRANSIENT THERMAL IMPEDANCE MODELS

The output of the STARS Visual System for the Powerex TBK7 (600V, 3000A) SCR is illustrated in Figure 1. This is the standard five points per decade table evaluated in Reference [5] where the validity of linear interpolation was carefully investigated. Note the wide range of the independent variables (On-State Current and time) required to permit surge and overload ratings to be evaluated and a reasonable range of applications simulated.

This forward voltage drop and transient thermal impedance are represented in STARS by Data Tables as illustrated in Figure 1. For a given instantaneous anode current the corresponding On-State Forward Voltage Drop is looked up and linearly interpolated. For a given point in time the Transient Thermal Impedance is looked up and also linearly interpolated. The calculation of junction temperature in STARS and in the past has been accomplished with superposition[6]. Here one evaluation point on the power waveform is chosen and the waveform prior to this point is divided into a number of time slices. The power dissipation for each segment is calculated and the temperature rise for that section is calculated by looking up the transient thermal impedance from the evaluation point to the leading edge (heating) and falling edge (cooling). All of the sections are then summed to get the junction temperature rise for the one chosen evaluation point in time.

The output of a superposition calculation module in STARS for an Arbitrary Current Waveshape is illustrated in Figure 2 for the Powerex T9G0 (2400V, 1000A) SCR. Here the first two (or three) columns is a text input file and the junction temperature is calculated and placed in the 4th column.

   
Figure 1.
STARS Visual Window with TBK7 Tabular Data
 

 
Figure 2. Convolution of Arbitrary Waveform 50mmT9G0 SCR Itm(peak)=17,000A Tj(peak)=369°C

 
Figure 3.
Parametric Data for On-State Vtm and Transient Thermal Impedance 77mm TBK7 SCR

A. SCR and Diode On-State Modeling:

The model parameters described in References 1, 2 and 3. are illustrated in Figure 3. Under the Data menu in STARS, there is a choice of the Vf, Vtm Model Parameters i.e., the ABCD or MNOPQ… . The parameters are quickly evaluated and placed in a parameter file with the format as illustrated. The range of anode current can also be chosen (last two numbers in first row).

B. Transient Thermal Impedance Models:

The transient thermal impedance can be modeled from physical parameters of the SCR or Diode as described in [3]. The resulting ladder network however is too complex to be used as a model but can be simplified with good accuracy by regression of the resulting transient thermal impedance curve by four or five series connected parallel RC stages.[2] The R(n), Tau(n) for the Transient Thermal Impedance Model Parameters are also automatically calculated in STARS for the selected device and placed into a file as displayed at the bottom of Figure 3.



III. THE STARSIMulation SYSTEM

The newest generation of the STARS Visual Rating System, called STARSim, provides a transition from curve or table based data to the Vf, Vtm and Transient Thermal Impedance Models. The idea is to make calculations with the old tables and the new modeling equations to transfer the confidence from the old standard calculations to the new equation modeling methods.

The first step is to compare the model generated curve, for transient thermal impedance, to the table of values curve as illustrated in Figure 4. This is actually two curves which provides, by the overlap of the curves, a quick evaluation as to how well the table generated curves match the model equation curves.

 
Figure 4.
Transient Thermal Impedance Compared in STARSim TBK7

Figures 5a. and 5b. illustrate the STARSim Arbitrary Waveform Evaluation Module for the Powerex T9G0 SCR. This module takes an arbitrary wave shape and evaluates the junction temperature rise using the Vf,  Vtm and transient thermal Modeling Equations.

The arbitrary waveform data can be viewed in Figure 5b. by clicking on the curve in Figure 5a, Note this is the same current waveform used in Figure 2. which used the old Vf or Vtm and transient thermal impedance tables and superposition.

Comparison of the peak junction temperatures which occurs at 4.5 msec is 369°C in STARS (Figure 2) compared to 353°C in STARSim (Figure 5b) indicating agreement of the two methods within 4%.

 
Figure 5a.
STARSim Arbitrary Current Waveform Tj Evaluation T9G0

 
Figure 5b.
Arbitrary Current Waveform Table, for T9G0 Tj(peak)=353°C

Next we will describe how the junction temperature is calculated in STARSim . This is done by hard coding in Microsoft Visual BASIC. While some readers may not be interested in tackling the time consuming task of hard coding this section will provide insight into: 1. How the modeling procedures are implemented and 2. How the SPICE modeling works which will be described in the next section of this paper.

IV. HARD CODE SCR DIODE MODELING

Hard coding for circuit simulation consist of iterating very small steps of time, resulting in small steps in the applied source voltage. Next, depending upon the state of the SCR(s), (see the state diagram; Figure 6), calculations are made as to how this voltage step affects the other voltage and currents in the circuit.

Equations will be presented in next section. Of interest here are the instantaneous anode current through the SCR(or Diode) and the resulting on-state forward voltage drop, VTM or VFM, using the ABCD or MNOPQ.. models. The product of the device voltage and current then represents the instantaneous power dissipation in the device.

Finally, this power dissipation is used as a current in the equivalent analog thermal circuit of the device (several stages of parallel RC circuits) and the calculated voltage across the RC stages is analogous to the virtual junction temperature of the device.



A. Modeling the Electrical Circuit:
Modeling the circuit consists of first solving the differential equations which describe the load current change with a step change in the applied voltage to the load. The initial load current must also be included. Considering the load resistance (RL) and reactance (LL) the equation is:

(1)
As we will iterate with very small time increments, each time increment will result in a new instantaneous load increment by the equation:

(2)
The SCR Model subroutine follows the State Diagram given in Figure 6. The state of the thyristor is switched to the on-state based upon the voltage across the device and the gate signal and to the off-state by the current decreasing to zero through the device.

The code to model the SCR is given as Code Fragment 1 in the Appendix. Note that line 1400 is the calculation of instantaneous current as described by Equation (2). The EX variable is the exponential term value for the fixed time step t =K which is calculated prior to the time iteration loop. The computer code calculation of instantaneous current is:

IL=IL*Ex + (VL/RL)*(1! - Ex) ‘SCR (Load) Current (3)

B. Modeling the Thermal Circuit:
The SCR or Diode junction temperature is modeled using the electrical - thermal analog circuit of 4 or 5 series connected parallel RC sections. As described in Reference [2] this provides a very accurate thermal model over a wide range of time. In this model the voltage is analogous to temperature and current is analogous to power dissipation in the device. The resistors and capacitors represent the major thermal time constants of the device. As described [2] The R and C parameters can be found by regression of the transient thermal impedance curve of the device. The equation for a given RC section (n) is then:

As we iterate with very small time increments, each time increment results in a new set of instantaneous temperatures(voltages), and switching from the electrical analog to thermal symbols the equation for the temperature contribution of a given RC section is:

(4)
The resulting BASIC code fragment shown below

(5)
was described briefly in Reference [3] is given in Code Fragment 1 (lines 1636 to 1690) of the Appendix. Note how Equation (5) is implemented in code as Equation 6:
‘Note T(n) = Initial(delta) Temp of Tau(n) Section
'Note S(1) = 0 = OFF S(1) = 1 = ON
Tjc = 0
FOR N = 1 TO 4
IF S(1) = 0 THEN T(N) * EXT(N) ‘Cooling Only
IF S(1) = 1 THEN T(N) = PT(1) * RT(N) * (1 - . _ EXT(N)) + T(N) * EXT(N) ’(6)
Tjc = Tjc + T(N) ‘Sum of the four sections
NEXT N 'Note S(1)=1 Switch on S(1)=0 Switch off
Tj = Tjc + Pave * (RTCS + RTSA) + Tamb

The modeling waveforms in STARSim for an SCR in an AC switch are illustrated in Figure 7. The device is the TBK7 77mm 600 Volt SCR, The conditions are a two cycle 60,000 Ampere surge at a 90 degree conduction angle The source is 600 (peak) Volts and the load a 0.06 Ohm resistor in series with a 6µHenry di/dt inductor .The waveforms include: the applied AC line voltage(VdVr), SCR current(IL), SCR On-State Voltage drop(Vtm), instantaneous power dissipation(Pt) and last, but certainly not least: the instantaneous virtual junction temperature of the SCR (Tj).

Note that the peak junction temperature occurs around 23.6 msec and is 121.6°C. This compares well with 122.9°C in STARSim Arbitrary Waveform Module at 23.9msec. The old superposition table look up in STARS resulted in 117.6°C. The agreement of the two approaches is 3.3%.

 
Figure 7.
STARSimulator "Hard Code" 60kA Waveforms TBK7

Next we switch our attention from Hard Coding in a higher level language to the simulation of the SCR and Diode in Power Electronic Circuits with PSPICE.

V. PSPICE SCR (Diode) MODELING

A. PSPICE Model Electrical Circuit:
The electrical circuit in Figure 9. consists of a sinusoidal source, load resistor and the SCR Vtm model. The SCR model consist of a Diode (D_SCR), a voltage controlled voltage source (E_Vtm), and a zero voltage source current sensor (V_Pwr_Sen). The SCR current is: Itm = I*(V_Pwr_Sen).

The SCR Electrical Model provides the correct instantaneous on-state forward voltage drop for any value of instantaneous anode current. This is achieved by sensing the anode current I*(V_Pwr_Sen) and using the ABCD or MNOPQ… Vtm Models to force the correct instantaneous Vtm by the Value statement in (E_Vtm). This is illustrated by the circuit in Figure 9. and the Net List Code Fragment 2. in the Appendix.



B. PSPICE Thermal "Circuit":
The thermal circuit is also illustrated in Figure 9. It consists of a voltage controlled current source (G_PowerM). The value of this current source is controlled by the calculated instantaneous power dissipation in the SCR i.e., V(A) * I*(V_Pwr_Sen). The current source then drives the analog electrical to thermal equivalent circuit for the device i.e., four or five sections of R-C parallel stages. This is the transient thermal impedance model of the SCR or Diode.

This circuit is easy to follow in Figure 9. and the Net List in Code Fragment 2. where R1,C1 connect from the Tj Node via Tj1 to R2,C2 etc. Note that the same forced charging and natural discharging of the thermal capacitors is occurring as described in the hard code section but now PSPICE is taking care of the differential equations.

The PSPICE waveforms for the Powerex T9G0 withstanding a 17,000 Ampere single cycle surge is illustrated in Figure 8. Note that not only is the current waveshape through the SCR, I(Rload) displayed but also the Forward Voltage Drop, Vtm , Power Dissipation, Pt and the resulting Virtual Junction Temperature Waveforms. The power dissipation is displayed as a current where (1 Amp = 1 Watt). Temperature is shown as the analog voltage Vtj (1 Volt=1°C) at the node Tj . The peak junction temperature is 356°C at 4.25msec versus 353°C for the hard code Arbitrary Waveforms in Figures 5a and 5b resulting in an difference of 0.8%.

Figure 8.
PSPICE SCR Model Waveforms T9G0 Tj(Peak)=353°C

C. PSPICE Model Sub Circuit:
The model for the SCR or Diode is diagrammed in Figure 9, and the net list presented in Code Fragment 2, is implemented as a Subcircuit in PSPICE. Note in the listing an asterisk at the beginning of the line makes that line a comment. I.e., the Vtm ABCD Model is active and the MNOPQ… model is "asterisk out". The thermal circuit has five stages where the values of the analog components are quickly determined from the transient thermal impedance as described earlier. All of the parameter values for A, B, M, N, … Rn, TAUn Cn=TAUn/Rn] are inserted via PSPICE Parameter Statements noted by the curly brackets around the label e.g., {C}. Figure 10. is the TBK7 SCR PSPICE waveforms with the same current waveform as the hard code waveform Figure 7. Note the peak junction temperature occurs at 23.8 msec vs 23.6 msec with a value of 122.1°C vs 121.6°C. This is excellent agreement with an difference of 0.8% and 0.4% respectfully.

Figure 9.
PSPICE SCR Model Circuit; Including Gate
Control and Junction Temperature

Figure 10.
PSPICE SCR Model Waveforms TBK7 2 Cycle 90° Conduction Angle 60 kA Surge Tj(Peak)=122.6°C

VI. CONCLUSIONS

Circuit modeling is an important tool in the design of power electronic applications. This paper has demonstrated that the new mathematical models of forward voltage drop and transient thermal impedance are quite capable of providing quick accurate simulations of SCRs and Diodes in power electronic applications.

The methodology to use these models, in both "digital" Hard Code, and "analog" PSPICE simulations have been described in detail including graphs of the output waveforms and quantitative evaluations of the overall accuracy. The agreement between the new mathematical models method and old, classical, but tedious, curve look up, superposition methods of calculating virtual junction was good. There was also shown to be excellent agreement between the Hard Code and the new Electro Thermal PSPICE SCR Model implementing the on-stat drop and thermal models.

This information hopefully will assist the power electronic engineer achieve more innovative, challenging and reliable designs in the many areas where power the electronics can provide useful functions to the public.



VII. REFERENCES

[1] J. W. Motto Jr., William H. Karstaedt, Jerry M. Sherbondy, Scott G. Leslie, "Thyristor(Diode) On-State Voltage, The ABCD Modeling Parameters Revisited Including Isothermal Overload and Surge Current Modeling" IEEE-IAS Annual Meeting, San Diego, California, October 1996

[2] J. W. Motto Jr., William H. Karstaedt, Jerry M. Sherbondy, Scott G. Leslie, "Thyristor(Diode) Transient Thermal Impedance Modeling Including the Spatial Temperature Distribution During Surge and Overload Conditions", IEEE-IAS Annual Meeting, Orlando Florida, October 1995

[3] J. W. Motto Jr., "Thyristor(Diode) Transient Thermal Impedance Modeling and Verification for Inductive Load Applications" , IEEE-IAS Annual Meeting, Denver, Colorado, October 1994

[4] J. W. Motto Jr., "Thyristor Steady State Current Ratings Past, Present and Future" , IEEE-IAS Annual Meeting, Toronto, Canada, October 1993

[5] J. W. Motto Jr., "Computer Aided Analysis of Thyristor Current Ratings" I&GA, IEEE Group Meeting Pittsburgh Pennsylvania, October 1967

[6] F. W, Gutzwiller and T. P. Sylvan "Power Semiconductors Under Transient and Intermittent Loads", AIEE Winter General Meeting, New York, New York January 31, 1960

[7] W. E. Newell "Transient Thermal Analysis of Solid State Power Devices - Making the Dreaded Process Easy", IEEE Transactions on Industry and General Applications Vol IA-12 July August, 1976

[8] D.E. Piccone, L.D. Eriksson, Dr. D.J. Urbanek, W.H. Tobin and I.L. Somos, "A Thermal Analog of Higher Accuracy and Factory Test Method for Predicting Thyristor Fault Suppression Ratings" IEEE-IAS Annual Meeting October 1988

[9] A. R. Hefner Jr. "A Dynamic Electro.-Thermal Model for the IGBT" IEEE Transactions on Industry Applications Vol 30 No 2 March/April 1994

[10] C. D. Mohler, "Digital Computer Calculation of Rectifier and Silicon Controlled Rectifier Ratings", AIEE Winter General Meeting, , New York, New York, January 30, 1962

VIII. APPENDIX: CODE FRAGMENTS of HARD CODE and PSPICE MODELS

1000 Rem ***** Switch Model { S(N)=0 Switch OFF - S(N)=1 Switch ON } ****
1010 'Switch State & Current Condition - For AC Switch Check Turn-off First
1040 If S(1) = 1 Then If IL < 0 Then S(1) = 0: VS(1) = E '**Switch #1 ON -> OFF
1050 If S(2) = 1 Then If IL > 0 Then S(2) = 0: VS(2) = E '**Switch #2 ON -> OFF
1060 If S(1) = 0 Then If S(2) = 0 Then VS(1) = E: VS(2) = E
1110 'Switch State & Voltage & Gate Condition
1130 If S(1) = 0 Then If VS(1) > 0 Then If EC > 0 Then S(1) = 1: VS(1) = 5: VS(2) = -VS(1) '**Switch #1 OFF -> ON
1140 If S(2) = 0 Then If VS(2) < 0 Then If EC < 0 Then S(2) = 1: VS(2) = -5: VS(1) = -VS(2) '**Switch #2 OFF -> ON
1150 If S(1) = 0 Then If S(2) = 0 Then VS(1) = E: VS(2) = E: VL = 0
1300 Rem *** Calculate Load Voltages and Currents Based on Switch Conditions **
1350 If S(1) = 1 Then VL = E - VS(1)
1360 If S(2) = 1 Then VL = E - VS(2)
1400 IL = IL * Ex + (VL / RL) * (1! - Ex) '<=============<<<< Transient Solution for Load Current
1430 If S(1) = 1 Then If IL > 0 Then VS(1) = ABCD(1) + ABCD(2) * Log(IL) + ABCD(3) * IL + ABCD(4) * Sqr(IL): VS(2) = VS(1)
1450 If S(2) = 1 Then If IL < 0 Then VS(2) = -ABCD(1) - ABCD(2) * Log(Abs(IL)) - ABCD(3) * IL - ABCD(4) * Sqr(Abs(IL)): VS(1) = VS(2)
1510 If S(1) = 1 Then Vtms = VS(1): PT(1) = VS(1) * IL 'Thyristor Power Disipation
1520 If S(1) = 0 Then Vtms = 0: PT(1) = 0
1530 PL = Abs(IL * VL): PLsum = PLsum + PL 'Sum Load Power
1600 PTsum = PTsum + PT(1) 'Sum of thyristor pwr disp for PAVE
1602 If S(1) = 1 Then ITave = ITave + IL: ITrms = ITrms + IL ^ 2
1604 ILrms = ILrms + IL ^ 2
1620 If Tim = ts Then If InStr(O$, "t") Then Print "P"; Int(Pave);: Rem QB LOCATE 54 / LD, 9:
1622 If Tim = ts Then If InStr(O$, "l") Then Print "PLOAD="; Int(PLave);: Rem LOCATE 54 / LD, 9:
1630 If Tim = ts Then If ST$ = "s" Then For n = 1 To Mthz: T(n) = RT(n) * Pave: Next n 'Note Initial Conditions for SS Valid after 1st run TS->TE
1632 If Tim = ts Then If ST$ = "t" Then For n = 1 To Mthz: T(n) = 0: Next n
1636 Tjc = 0
1640 For n4 = 1 To Mthz 'NOTE T(n4)=Initial(delta)Temp of Tau(n4) 1ma+1ma per 100v
1650 If S(1) = 0 Then If S(2) = 0 Then T(n4) = (0.001 + 0.00001 * Abs(VS(1))) * Abs(VS(1)) * RT(n4) * (1 - EXT(n4)) + T(n4) * EXT(n4)
1660 If S(1) = 0 Then If S(2) = 1 Then T(n4) = T(n4) * EXT(n4) 'Note PT(1)=0 cooling only
1670 If S(1) = 1 Then T(n4) = PT(1) * RT(n4) * (1 - EXT(n4)) + T(n4) * EXT(n4)
1680 Tjc = Tjc + T(n4)
1690 Next n4
1700 Rem xxx ST$ = "s" Then tj = Tjc + Pave * Dcy * (RTcs + RTsa) + Tamb(0)
1701 If ST$ = "s" Or ST$ = "e" Then tj = Tjc + Pave * Dcy * (RTcs + NDHS * RTsa * (1 - Exp(-Tim / TauSA))) + Tamb(0) 'if RTsa=0 then RTcs=0
1710 If ST$ = "t" Then tj = Tjc + Tamb(0)
1990 Return
Code Fragment 1. Visual BASIC "Hard Code" for SCR Model

 

******************* POWEREX TBK7 SCR MODEL ****************************
*ANODE ANODE input node
* CATHODE CATHODE input node
* Tj Junction temperature output node
* AMB_T Starting ambient temperature **************************************************************************
.SUBCKT PRX_TBK7 A K G Kg Tj PARAMS: AMB_T=27

***************** Vtm Model *** ****** April 15, 1997 ************************
D_SCR A A1 Dmod ;Provide SCR with Reverse Blocking [Designed for Vf = K = 0.4Volts]
EVtm A1 K1 Value= {{A} + {B} * log(Abs(I(V_Pwr_Sen))) + {C} * I(V_Pwr_Sen) + {D} *PWR(I(V_Pwr_Sen),.5)-0.4}
**EVtm A1 K1 Value= {{M} + {N} *PWR(I(V_Pwr_Sen),.25) + {O} *PWR(I(V_Pwr_Sen),.5)+ {P} *PWR(I(V_Pwr_Sen),.75)+ {Q} * I(V_Pwr_Sen)-0.4}
Sscr K1 k2 G2 KG SMOD ;Voltage Controlled Switch
Rg G KG 1 ;Gate Resistor
V_Pwr_Sen K2 K
************ **Thermal Model ****** March 11, 1997 ****************************
G_PowerM 0 Tj Value= {I(V_Pwr_Sen) * (V(A1,K1)+0.4)} ; <=====I(power) for ThZ Model
R1 Tj Tj1 {R1}
C1 Tj Tj1 {TAU1/R1}
R2 Tj1 Tj2 {R2}
C2 Tj1 Tj2 {TAU2/R2}
R3 Tj2 Tj3 {R3}
C3 Tj2 Tj3 {TAU3/R3}
R4 Tj3 Tj4 {R4}
C4 Tj3 Tj4 {TAU4/R4}
R5 Tj4 0 {R5}
C5 Tj4 0 {TAU5/R5}
* * * * MODELS * * * *
.Model Dmod D(Eg = 0.9 Cjo = 1.0uf) ; ["Ideal" except for: Vf = K = 0.4Volts]
.Model SMOD VSWITCH(RON=1E-7 ROFF=1E+5 VON=1 VOFF=0)
.ENDS ;----- END OF SCR SUBCIRCUIT MODEL**********************************
Code Fragment 2. Net List for PSPICE SCR Modeling using Sub Circuit

 

News Flash
POWEREXtra
Technical Spotlight