# SIMSCAPE Electrical Modelling of the IGBT with Parameter Optimization Using Genetic Algorithm

Authors: Mohamed Baghdadi; Elmostafa Elwarraki; Naoual Mijlad; Imane Ait Ayad
Journal: Journal of Electrical and Computer Engineering
 (2021)

---

## Abstract

The concept introduced by MathWorks in the Simscape product is the link representation between the SIMSCAPE library components that correspond to physical connections transmitting power. In this paper, a power insulated-gate bipolar transistor (IGBT) model using MATLAB graphical software is reproduced. An electrical IGBT behavior model using the Simscape Electronics library components is developed and analyzed. This model is parameterized using the constructor datasheet to ensure a good representation of the dynamic and static IGBT behaviors. An extraction and optimization studies of the IGBT model parameters using a stochastic algorithm implemented in Matlab are presented. The proposed method is based on the Genetic Algorithm (GA) to perfectly extract and optimize the model parameters using the mathematical model circuit equations and the provided datasheet characteristics. A simulation in the Matlab/Simulink environment and a comparison with the experimental results for an IGBT device example are carried out to demonstrate the proposed model accuracy.


---

## Body

## 1. Introduction

For reducing design costs and increasing electronic power circuits’ reliability, the modeling of power semiconductor devices becomes indispensable. The model must capture the dynamic and static behavior of the power semiconductor devices to ensure the safe global circuit operating area. The IGBT is one of the most important power semiconductor devices widely used in different power electronic applications.In the literature, several IGBT models were developed at different complexity levels and for various uses. At this level, an accurate physical device description has been required in [1, 2]. Hefner’s model and a model list based on the Hefner work are examples of suitable mathematical models for circuit simulator implementation, which have proved to have a good accuracy and high reliability [3–5]. In [6, 7], the authors have proposed a complete and complex mathematical model using a complex equation for each model part, including the wide-base PNP section, the MOSFET channel section, capacities, the JFET effect, and the thermal effect. The difference between the analytical models is in the simplifications made in order to reduce the simulation time. To have more accurate results, numerical methods were employed to solve the 2 D and 3 D ambipolar diffusion equations. However, the numerical models’ implementation requires special simulators. These model types are specially developed to optimize IGBT structure technologies [8].For macromodels, they are used when high accuracy is not needed and the representation of some phenomena is ignored. Thus, the physical IGBT mechanisms are not considered in the models. Indeed, the fitted characteristics of the IGBT component are translated to an equivalent electrical circuit [9, 10].The literature reported several works on the automatic extraction and optimization of IGBT model parameters. In [11], the authors developed the most suitable physical models and reference for the IGBT, and the IMPACT software has been elaborated to extract the IGBT model parameters. The IMPACT software has included five programs, which have extracted the 20 physical and structural parameters from the most recent version of the developed IGBT model. Each program has provided automatic algorithms to perform the extraction, using experimental measurements for each parameter set. The behavioral model of Oh and El Nokali implemented in Saber [12] has been reproduced by Kang et al. who have developed a new extraction algorithm to determine its parameters [13]. The algorithm is based on the availability of experimental switching loss data and the manufacturer-supplied characteristics, while the Matlab optimization toolbox is used to find model parameters automatically.In [14], the IGBT model based on the finite-elements (FEM) physics has been implemented in the circuit simulator IsSPICE. The proposed extraction procedure is the numerical optimization algorithm called “Simulated Annealing” (SA). Its implementation is carried out in the Matlab programming environment, which allows a data exchange with IsSPICE.This paper aims to model the IGBT device using the graphical Matlab software. The proposed electrical IGBT model is developed using Pspice [15, 16] and Saber simulators [17] in the previous contributions. The main contribution of this work is the automatic identification and optimization of the model parameters based on a stochastic algorithm. The nonlinear equations describing the IGBT model for each operating zone are rewritten as an optimization problem for their model parameters, which is solved using the Genetic Algorithm (GA) with the datasheet characteristics. The proposed technique is easily implemented in the Matlab environment, and the simulation results proved its high reliability. The rest of this paper is organized as follows: the electrical IGBT model is developed in Section 2, while Section 3 describes the model parameter optimization using GA. In Section 4, simulation results of the proposed model are discussed. Finally, conclusions are given in Section 5.


## 2. IGBT Design

### 2.1. Electrical IGBT Model

The insulated-gate bipolar transistor (IGBT) is a hybrid transistor, which includes the Metal Oxide Semiconductor Field Effect Transistor (MOSFET) at the input and the Bipolar Junction Transistor (TBJ) at the output. Consequently, the IGBT has low power control as an advantage of the MOSFET and presents the benefits of the TBJ which are low conduction losses and high-voltage handling [18]. The established model can predict the IGBT behavior for both dynamic and steady state. Besides, it can be easily parameterized using manufacturer data, and it represents an optimum compromise between accuracy and simulation complexity [17].The static MOSFET model is constituted by(i)
DG: a diode which allows the passage from the blocked state to the on state.(ii)
VTH: DC voltage source. It represents the threshold voltage.(iii)
RGE: leakage gate-emitter resistance.(iv)
Jmos: a controlled current source which represents the drain current (Imos).(v)
RDN: drain resistance which characterizes the resistive aspect of the MOSFET channel.(vi)
DDC: an antiparallel body diode which allows the passage from the linear zone to the saturated zone. It is represented by a binary resistance RONDC/ROFFDC.The MOSFET current is controlled byVGE voltage, according to three cases:(i)
VGE < VTH: the transistor is off, so the channel is not constructed(ii)
VGE + VCE0 − VTH > VCE:(1)Imos=KpVGE−VTH∗VCE−VCE0−VCE−VCE02/21+θVGE−VTH(iii)
VGE + VCE0 − VTH < VCE:(2)Imos=Kp2VGE−VTH21+θVGE−VTH,where Kp is the transconductance and θ is the correction factor, which takes into consideration the transverse field in the MOSFET channelFor the bipolar part, the static TBJ model is composed by(i)
DE and DP: the base-emitter and base-collector junctions(ii)
VCE0: conduction threshold voltage of the base-emitter junction(iii)
JPNP: a controlled current source which represents the current gain of the bipolar transistor(iv)
β: the bipolar gain


### 2.2. Dynamic IGBT Model

To model the dynamic IGBT behavior, other elements are added to the static model. These elements include nonlinear terminal capacitances and the voltage-controlled current source that models the load responsible of the tail currentIC at the turn off IGBT state.Constructors provide the measurement capacitors between two terminals by eliminating one or two capacitors. The input capacitance (CIE), reverse transfer capacitance (CRE), and output capacitance (COE) allow us to determine the interelectrode capacitance values CGE, CCE, and CGC using the following equations [19]:(3)CIE=CGE+CGCCOE=CCE+CGCCRE=CGC⟹CGE=CIE−CRECCE=COE−CRECGC=CRE.Figure1 shows the evolution of CGE, CCE, and CGC capacitances versus VCE. We note that the CGE capacity remains almost constant while the CGC and CCE capacities vary strongly nonlinear following the device operating area. These two capacities change values around the same threshold voltage VCET0.Figure 1
Curves of capacitances between terminals for the commercial IGBT (IRGBC20U) [18]. CGE capacity: its value is given by the system of equations (3).CGC capacity: this capacity can be constituted of two capacities in series: COXD oxide capacity that is fixed by the constructor and CGCJ depletion capacity. In Figure 1, the CGC capacity changes value from a certain VCE threshold voltage VCET. We can use the diode model [16] to represent the nonlinearity of the CGC capacity between two values CGCmin and CGCmax (Figure 2).Figure 2
Model of the nonlinearCGC capacity [15].TheCGCmax capacity is represented by a voltage-controlled current source, with a gain given by the following expression:(4)GGC=RonGCCGCmaxLGC.CCE capacity: this capacity, as it seen in Figure 1, varies also as nonlinear. It is represented in the same manner as that for the CGC capacity. The coupling coefficient GCE is given by the following expression:(5)GCE=RonGCCCEmaxLGC.The bipolar part presents the responsible charge of the tail currentIC, and these charges can be represented only by the PNP base-emitter junction storing charges. Thereafter, these charges are modelled by a controlled source (Figure 3), whose expression is as follows:(6)Gt=RonDECBELBE,where CBE represents the bipolar part storing charges, with CBERonDE = toff, where toff is deduced from the device datasheets.Figure 3
IGBT electrical model.The studied IGBT model is shown in Figure3.


### 2.3. The IGBT Model Parameter Extraction

This section aims to develop and show the extraction method of the IGBT model parameters. To identify the model parameters in the static case, the main idea is to take out the circuit equations according to the diode states (Figure4). Indeed, the diodes states define the different IGBT operation zones.Figure 4
Decomposition of the static characteristic.For the Ohmic zone, all diodes are turned on. The resolution of the equations system described the static model gives the expression of theVCE voltage versus the IC current:(7)VCE=VCE0−RONDP∗RONDCRONDP+RONDC+RDN∗Imos+RONDE+RONDP∗RONDC+RDN∗1−GP∗RONDERONDP+RONDC+RDN∗IC.For the quasisaturation zone, all the diodes are turned on except theDP diode. The expression of the VCE voltage versus the IC current becomes(8)VCE=VCE0−RONDC∗ROFFDPRONDC+ROFFDP+RDN∗Imos+RONDE+ROFFDP∗RONDC+RDN∗1−GP∗RONDEROFFDP+RONDC+RDN∗IC.For the saturation zone, it is characterized by the blocking diodeDDC of the MOSFET part and the blocking diode DP of the bipolar part. The expression of the VCE voltage becomes(9)VCE=VCE0−ROFFDC∗ROFFDPROFFDC+ROFFDP+RDN∗Imos+RONDE+ROFFDP∗ROFFDC+RDN∗1−GP∗RONDEROFFDP+ROFFDC+RDN∗IC.The parameters of the studied model are deduced from the technical datasheet and the equations established in [15]. The IGBT parameters are recorded in Table 1 for the commercial IGBT (IRGBC20U).Table 1
IGBT electrical parameters (IRGBC20U).
ParametersValuesStatic parametersVCE01.9 VRONDE8.66 10−2 ΩRDN2.3073 10−3 ΩRONDC2.479 10−1 ΩROFFDC108 ΩRONDP3.42 10−3 ΩROFFDP2.373 107 ΩGP1.154Ω−1VTH5.3 VKP1.766Θ8.47 10−2RGE2 108Β10−1Dynamic parametersCGE3 10–10 FCGCmin, RGC, LGC,RONGCVCET6 10−12 F, 10−1 Ω, 10−10 H, 103 Ω, 20VGGC2.44 10−3FCCE, GCE5.9 10−13 F, 3.06 103 Ω−1Gt910Ω−1


## 2.1. Electrical IGBT Model

The insulated-gate bipolar transistor (IGBT) is a hybrid transistor, which includes the Metal Oxide Semiconductor Field Effect Transistor (MOSFET) at the input and the Bipolar Junction Transistor (TBJ) at the output. Consequently, the IGBT has low power control as an advantage of the MOSFET and presents the benefits of the TBJ which are low conduction losses and high-voltage handling [18]. The established model can predict the IGBT behavior for both dynamic and steady state. Besides, it can be easily parameterized using manufacturer data, and it represents an optimum compromise between accuracy and simulation complexity [17].The static MOSFET model is constituted by(i)
DG: a diode which allows the passage from the blocked state to the on state.(ii)
VTH: DC voltage source. It represents the threshold voltage.(iii)
RGE: leakage gate-emitter resistance.(iv)
Jmos: a controlled current source which represents the drain current (Imos).(v)
RDN: drain resistance which characterizes the resistive aspect of the MOSFET channel.(vi)
DDC: an antiparallel body diode which allows the passage from the linear zone to the saturated zone. It is represented by a binary resistance RONDC/ROFFDC.The MOSFET current is controlled byVGE voltage, according to three cases:(i)
VGE < VTH: the transistor is off, so the channel is not constructed(ii)
VGE + VCE0 − VTH > VCE:(1)Imos=KpVGE−VTH∗VCE−VCE0−VCE−VCE02/21+θVGE−VTH(iii)
VGE + VCE0 − VTH < VCE:(2)Imos=Kp2VGE−VTH21+θVGE−VTH,where Kp is the transconductance and θ is the correction factor, which takes into consideration the transverse field in the MOSFET channelFor the bipolar part, the static TBJ model is composed by(i)
DE and DP: the base-emitter and base-collector junctions(ii)
VCE0: conduction threshold voltage of the base-emitter junction(iii)
JPNP: a controlled current source which represents the current gain of the bipolar transistor(iv)
β: the bipolar gain


## 2.2. Dynamic IGBT Model

To model the dynamic IGBT behavior, other elements are added to the static model. These elements include nonlinear terminal capacitances and the voltage-controlled current source that models the load responsible of the tail currentIC at the turn off IGBT state.Constructors provide the measurement capacitors between two terminals by eliminating one or two capacitors. The input capacitance (CIE), reverse transfer capacitance (CRE), and output capacitance (COE) allow us to determine the interelectrode capacitance values CGE, CCE, and CGC using the following equations [19]:(3)CIE=CGE+CGCCOE=CCE+CGCCRE=CGC⟹CGE=CIE−CRECCE=COE−CRECGC=CRE.Figure1 shows the evolution of CGE, CCE, and CGC capacitances versus VCE. We note that the CGE capacity remains almost constant while the CGC and CCE capacities vary strongly nonlinear following the device operating area. These two capacities change values around the same threshold voltage VCET0.Figure 1
Curves of capacitances between terminals for the commercial IGBT (IRGBC20U) [18]. CGE capacity: its value is given by the system of equations (3).CGC capacity: this capacity can be constituted of two capacities in series: COXD oxide capacity that is fixed by the constructor and CGCJ depletion capacity. In Figure 1, the CGC capacity changes value from a certain VCE threshold voltage VCET. We can use the diode model [16] to represent the nonlinearity of the CGC capacity between two values CGCmin and CGCmax (Figure 2).Figure 2
Model of the nonlinearCGC capacity [15].TheCGCmax capacity is represented by a voltage-controlled current source, with a gain given by the following expression:(4)GGC=RonGCCGCmaxLGC.CCE capacity: this capacity, as it seen in Figure 1, varies also as nonlinear. It is represented in the same manner as that for the CGC capacity. The coupling coefficient GCE is given by the following expression:(5)GCE=RonGCCCEmaxLGC.The bipolar part presents the responsible charge of the tail currentIC, and these charges can be represented only by the PNP base-emitter junction storing charges. Thereafter, these charges are modelled by a controlled source (Figure 3), whose expression is as follows:(6)Gt=RonDECBELBE,where CBE represents the bipolar part storing charges, with CBERonDE = toff, where toff is deduced from the device datasheets.Figure 3
IGBT electrical model.The studied IGBT model is shown in Figure3.


## 2.3. The IGBT Model Parameter Extraction

This section aims to develop and show the extraction method of the IGBT model parameters. To identify the model parameters in the static case, the main idea is to take out the circuit equations according to the diode states (Figure4). Indeed, the diodes states define the different IGBT operation zones.Figure 4
Decomposition of the static characteristic.For the Ohmic zone, all diodes are turned on. The resolution of the equations system described the static model gives the expression of theVCE voltage versus the IC current:(7)VCE=VCE0−RONDP∗RONDCRONDP+RONDC+RDN∗Imos+RONDE+RONDP∗RONDC+RDN∗1−GP∗RONDERONDP+RONDC+RDN∗IC.For the quasisaturation zone, all the diodes are turned on except theDP diode. The expression of the VCE voltage versus the IC current becomes(8)VCE=VCE0−RONDC∗ROFFDPRONDC+ROFFDP+RDN∗Imos+RONDE+ROFFDP∗RONDC+RDN∗1−GP∗RONDEROFFDP+RONDC+RDN∗IC.For the saturation zone, it is characterized by the blocking diodeDDC of the MOSFET part and the blocking diode DP of the bipolar part. The expression of the VCE voltage becomes(9)VCE=VCE0−ROFFDC∗ROFFDPROFFDC+ROFFDP+RDN∗Imos+RONDE+ROFFDP∗ROFFDC+RDN∗1−GP∗RONDEROFFDP+ROFFDC+RDN∗IC.The parameters of the studied model are deduced from the technical datasheet and the equations established in [15]. The IGBT parameters are recorded in Table 1 for the commercial IGBT (IRGBC20U).Table 1
IGBT electrical parameters (IRGBC20U).
ParametersValuesStatic parametersVCE01.9 VRONDE8.66 10−2 ΩRDN2.3073 10−3 ΩRONDC2.479 10−1 ΩROFFDC108 ΩRONDP3.42 10−3 ΩROFFDP2.373 107 ΩGP1.154Ω−1VTH5.3 VKP1.766Θ8.47 10−2RGE2 108Β10−1Dynamic parametersCGE3 10–10 FCGCmin, RGC, LGC,RONGCVCET6 10−12 F, 10−1 Ω, 10−10 H, 103 Ω, 20VGGC2.44 10−3FCCE, GCE5.9 10−13 F, 3.06 103 Ω−1Gt910Ω−1


## 3. Optimization of the IGBT Model Parameters Using Genetic Algorithm

### 3.1. Genetic Algorithm

In classical information treatment, the problems are solved in fixed ways before the introduction of heuristic methods using random processes. In fact, during the design phase of a system, it receives all the necessary data for its operating conditions as known at the time of its conception, which prevents its adaptation to unknown changing environmental conditions. Therefore, computer researchers are studying methods to allow systems to adapt spontaneously to new conditions. Hence, the heuristic and stochastic methods emerge such as evolutionary programming based on the Darwin theory. The theory of the species evolution is exposed under the external constraints influence, and the fixed systems’ existence is rejected.Among the evolutionary algorithms, the Genetic Algorithm (GA) is developed for optimization purposes. The genetic process is based on the natural selection theory and the evolution genetics. John Holland introduced the first formal model of genetic algorithms in 1976, which added an intelligence to a computer program with the crossing notion (genetic material exchange) and the mutation notion (genetic diversity source) [20]. Based on John Holland’s work, several research studies have been developed.Unlike the different optimization techniques, the GA is characterized by(i)
The use of an encoding describing the parameters(ii)
GA works on a population of points, instead of a single point(iii)
It uses only the function values under study(iv)
It uses probabilistic, not deterministic, transition rulesA genetic algorithm offers acceptable solutions for difficult problems in a reasonable calculation time. It tries to optimize the defined solution for a number of iterations and chooses the most optimized solution available for a problem at the program end. The initial creation of the population, the evaluation of the optimization condition (adaptation function or fitness function), selection, crossing, and mutation are the five fundamental functions of the algorithm.Figure5 shows the genetic algorithm diagram. The genetic algorithm process simulates the evolution of a population, which is subjected to a selection at each generation, and the recombination and mutation operators are applied. The selection based on the fitness function allows improving the population [20]. This algorithm does not require knowledge of the problem, which can be represented by a black box with inputs representing the variables and outputs representing the objective functions. Then, the algorithm manipulates the inputs in order to improve the outputs.Figure 5
Flowchart of the standard genetic algorithm (GA).


### 3.2. Model Parameter Optimization Using GA

The GA is a robust stochastic process that provides accurate solutions to the optimization problem. Different domains of power electronic field have been benefited from the exploitation of GA. However, for optimization of semiconductors models, the GA method has been rarely used [21, 22].Genetic algorithms belong to the evolutionary algorithms. They are based on natural selection principle and heredity. This concept is applied to an initial population of potential solutionsPt that gives several generations. For each generation, all individuals are evaluated according to an objective function [20]. Figure 6 shows the basic GA algorithm program.Figure 6
A simple genetic algorithm.In order to optimize the static parameters of the IGBT model, an approach is used to combine the manual extraction method that allows obtaining a first estimation of the set of parameters and then using the numerical optimization by the GA to extract the optimal set of parameters using the Toolbox “GA” in MATLAB [23].The flowchart of the proposed optimization approach is presented in Figure7.Figure 7
Organizational chart of the optimization approach applied to the IGBT model.The fitness function measures the error between manufacturer and simulated data. The used fitness function, for parameter evaluation, is defined by the sum of the deviation between the equation representing the manufacturer’s curve (desired value) and the equation of the curve containing the model parameters (calculated values). The fitness function is given as follows:(10)Ffitness=∑ICdata−ICM2,where ICdata is the value of the desired current given by the manufacturer and ICM is the current value calculated by the simulated model.Following that mentioned above, two programs are developed on MATLAB; the first one will use the mathematical equation of the transfer characteristic with equations (1) and (2) to optimize the coefficients KP, VTH, θ, and the bipolar gain β. The second program will use equations (7)–(9) with the mathematical equation of the output characteristic to calculate the other parameters VCE0, RONDE, RDN, RONDC, ROFFDC, RONDP, ROFFDP, and GP.The parameter values obtained by optimization and GA numerical information (data points and steps) are summarized in Tables2 and 3, respectively.Table 2
IGBT electrical parameters (IRGBC20U).
ParametersValuesVCE01.6119 VRONDE8.04 10−2 ΩRDN2.0 10−3 ΩRONDC3.893 10−1 ΩROFFDC4.690106 ΩRONDP8.0 10−3 ΩROFFDP9.361 105 ΩGP1.3.830 Ω−1VTH5.3 VKP1.6166Θ8.10 10−2RGE2 108Β1.34910−1Table 3
GA numerical information.
Data pointsSteps (s)Transfer characteristic14512.283353Output characteristic8510.748621 to 0.860017Figure8 shows the optimization results obtained for the two studied cases. A good convergence of the solution giving by the genetic algorithm in the treated cases can be seen.Figure 8
Convergence between manufacturer characteristics and GA curves giving optimal parameters. (a) Transfer characteristic. (b) Output characteristic.
(a)(b)The comparison between the obtained results and the simulation results without optimization is illustrated in Figure9. The improvement brought by the model in the reproduction of the output and transfer characteristics in the saturation regime is shown. The error between the curves obtained by optimization and the curves given by the manufacturer is low.Figure 9
Simulation results after GA optimization. Transfer characteristic and output characteristic.
(a)(b)


## 3.1. Genetic Algorithm

In classical information treatment, the problems are solved in fixed ways before the introduction of heuristic methods using random processes. In fact, during the design phase of a system, it receives all the necessary data for its operating conditions as known at the time of its conception, which prevents its adaptation to unknown changing environmental conditions. Therefore, computer researchers are studying methods to allow systems to adapt spontaneously to new conditions. Hence, the heuristic and stochastic methods emerge such as evolutionary programming based on the Darwin theory. The theory of the species evolution is exposed under the external constraints influence, and the fixed systems’ existence is rejected.Among the evolutionary algorithms, the Genetic Algorithm (GA) is developed for optimization purposes. The genetic process is based on the natural selection theory and the evolution genetics. John Holland introduced the first formal model of genetic algorithms in 1976, which added an intelligence to a computer program with the crossing notion (genetic material exchange) and the mutation notion (genetic diversity source) [20]. Based on John Holland’s work, several research studies have been developed.Unlike the different optimization techniques, the GA is characterized by(i)
The use of an encoding describing the parameters(ii)
GA works on a population of points, instead of a single point(iii)
It uses only the function values under study(iv)
It uses probabilistic, not deterministic, transition rulesA genetic algorithm offers acceptable solutions for difficult problems in a reasonable calculation time. It tries to optimize the defined solution for a number of iterations and chooses the most optimized solution available for a problem at the program end. The initial creation of the population, the evaluation of the optimization condition (adaptation function or fitness function), selection, crossing, and mutation are the five fundamental functions of the algorithm.Figure5 shows the genetic algorithm diagram. The genetic algorithm process simulates the evolution of a population, which is subjected to a selection at each generation, and the recombination and mutation operators are applied. The selection based on the fitness function allows improving the population [20]. This algorithm does not require knowledge of the problem, which can be represented by a black box with inputs representing the variables and outputs representing the objective functions. Then, the algorithm manipulates the inputs in order to improve the outputs.Figure 5
Flowchart of the standard genetic algorithm (GA).


## 3.2. Model Parameter Optimization Using GA

The GA is a robust stochastic process that provides accurate solutions to the optimization problem. Different domains of power electronic field have been benefited from the exploitation of GA. However, for optimization of semiconductors models, the GA method has been rarely used [21, 22].Genetic algorithms belong to the evolutionary algorithms. They are based on natural selection principle and heredity. This concept is applied to an initial population of potential solutionsPt that gives several generations. For each generation, all individuals are evaluated according to an objective function [20]. Figure 6 shows the basic GA algorithm program.Figure 6
A simple genetic algorithm.In order to optimize the static parameters of the IGBT model, an approach is used to combine the manual extraction method that allows obtaining a first estimation of the set of parameters and then using the numerical optimization by the GA to extract the optimal set of parameters using the Toolbox “GA” in MATLAB [23].The flowchart of the proposed optimization approach is presented in Figure7.Figure 7
Organizational chart of the optimization approach applied to the IGBT model.The fitness function measures the error between manufacturer and simulated data. The used fitness function, for parameter evaluation, is defined by the sum of the deviation between the equation representing the manufacturer’s curve (desired value) and the equation of the curve containing the model parameters (calculated values). The fitness function is given as follows:(10)Ffitness=∑ICdata−ICM2,where ICdata is the value of the desired current given by the manufacturer and ICM is the current value calculated by the simulated model.Following that mentioned above, two programs are developed on MATLAB; the first one will use the mathematical equation of the transfer characteristic with equations (1) and (2) to optimize the coefficients KP, VTH, θ, and the bipolar gain β. The second program will use equations (7)–(9) with the mathematical equation of the output characteristic to calculate the other parameters VCE0, RONDE, RDN, RONDC, ROFFDC, RONDP, ROFFDP, and GP.The parameter values obtained by optimization and GA numerical information (data points and steps) are summarized in Tables2 and 3, respectively.Table 2
IGBT electrical parameters (IRGBC20U).
ParametersValuesVCE01.6119 VRONDE8.04 10−2 ΩRDN2.0 10−3 ΩRONDC3.893 10−1 ΩROFFDC4.690106 ΩRONDP8.0 10−3 ΩROFFDP9.361 105 ΩGP1.3.830 Ω−1VTH5.3 VKP1.6166Θ8.10 10−2RGE2 108Β1.34910−1Table 3
GA numerical information.
Data pointsSteps (s)Transfer characteristic14512.283353Output characteristic8510.748621 to 0.860017Figure8 shows the optimization results obtained for the two studied cases. A good convergence of the solution giving by the genetic algorithm in the treated cases can be seen.Figure 8
Convergence between manufacturer characteristics and GA curves giving optimal parameters. (a) Transfer characteristic. (b) Output characteristic.
(a)(b)The comparison between the obtained results and the simulation results without optimization is illustrated in Figure9. The improvement brought by the model in the reproduction of the output and transfer characteristics in the saturation regime is shown. The error between the curves obtained by optimization and the curves given by the manufacturer is low.Figure 9
Simulation results after GA optimization. Transfer characteristic and output characteristic.
(a)(b)


## 4. Results and Discussion

Figure10 illustrates the experimental circuit used to perform the dynamic and static characteristics of the studied IGBT. The electrical IGBT environment is modelled by a wiring inductance LP and a wiring resistance RP. A single parasitic inductance is sufficient to represent the total circuit inductances [24]. The IGBT driver causes the rise and fall times of the VG impulsion. This interaction between the control section and the device is also included in the test circuit simulation.Figure 10
(a) Simulated circuit. (b) Photography of the IGBT test bench.
(a)(b)


### 4.1. Validation of the Static Behavior of the IGBT Model

To obtain the output characteristic, the simulation is performed as given in the manufacturer specifications in the following terms: the voltageVGE is maintained at 15 V, and the voltage VCE changes linearly from 1 V to 10 V during 25 μs. Figure 11 shows experimental (manufacturer data) and simulated static characteristics of the commercial IGBT IRGBC20U, obtained for test conditions of the manufacturer: (a) transfer characteristic; (b) output characteristic.Figure 11
Experimental (manufacturer data) and simulated static characteristics of the commercial IGBT IRGBC20U, obtained for test conditions of the manufacturer. (a) Transfer characteristic. (b) Output characteristic.
(a)(b)To perform the transfer characteristicIC (VGE), the VCE is fixed at 100 V while the VGE varies from 5 V to 20 V.Figure11 presents the manufacturer and simulation results of static characteristics. In the two curves, simulation results show good agreements compared to manufacturer data.


### 4.2. Validation of the Dynamic IGBT Model Behavior

The switching IGBT characteristics are affected by parasitic impedances, especially parasitic inductances. TheIC current decreases and VCE voltage increases at the turnoff switching phase owing to these parasitic inductances. According to Figures 12–14, it can be shown that the tail current IC and the VCE voltage peak caused by these parasitic inductances increase the duration of the switching operation. Then, the switching losses at the turnoff phase are increased. The model predicts with a good accuracy all dynamic IGBT characteristics, namely, slowly decaying current at turn off (tailing phenomenon characteristic of IGBTs due to bipolar transistor part) and the Miller effect in the gate voltage waveforms.Figure 12
VGE waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behaviour.Figure 13
VCE waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behavior.Figure 14
IC waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behaviour.


## 4.1. Validation of the Static Behavior of the IGBT Model

To obtain the output characteristic, the simulation is performed as given in the manufacturer specifications in the following terms: the voltageVGE is maintained at 15 V, and the voltage VCE changes linearly from 1 V to 10 V during 25 μs. Figure 11 shows experimental (manufacturer data) and simulated static characteristics of the commercial IGBT IRGBC20U, obtained for test conditions of the manufacturer: (a) transfer characteristic; (b) output characteristic.Figure 11
Experimental (manufacturer data) and simulated static characteristics of the commercial IGBT IRGBC20U, obtained for test conditions of the manufacturer. (a) Transfer characteristic. (b) Output characteristic.
(a)(b)To perform the transfer characteristicIC (VGE), the VCE is fixed at 100 V while the VGE varies from 5 V to 20 V.Figure11 presents the manufacturer and simulation results of static characteristics. In the two curves, simulation results show good agreements compared to manufacturer data.


## 4.2. Validation of the Dynamic IGBT Model Behavior

The switching IGBT characteristics are affected by parasitic impedances, especially parasitic inductances. TheIC current decreases and VCE voltage increases at the turnoff switching phase owing to these parasitic inductances. According to Figures 12–14, it can be shown that the tail current IC and the VCE voltage peak caused by these parasitic inductances increase the duration of the switching operation. Then, the switching losses at the turnoff phase are increased. The model predicts with a good accuracy all dynamic IGBT characteristics, namely, slowly decaying current at turn off (tailing phenomenon characteristic of IGBTs due to bipolar transistor part) and the Miller effect in the gate voltage waveforms.Figure 12
VGE waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behaviour.Figure 13
VCE waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behavior.Figure 14
IC waveforms of experimental and simulated dynamic IGBT (IRGBC20U) behaviour.


## 5. Conclusions

This study presents the electrical modeling of the IGBT using the graphical software of Matlab. Parameter optimization of the model based on a metaheuristic algorithm using GA was proposed. The established method used electrical equations of the model and measurement data provided by the manufacturer to perform the model parameter optimization. The method implementation was carried out using the Matlab/Simulink environment. Simulation of the electrical model showed a good precision and high reliability for the static and dynamic behavior of the used IGBT.



---

*Source: 6665384-2021-05-05.xml*