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

. 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.


Introduction
For reducing design costs and increasing electronic power circuits' reliability, the modeling of power semiconductor devices becomes indispensable. e model must capture the dynamic and static behavior of the power semiconductor devices to ensure the safe global circuit operating area. e 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][4][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. e 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. ese 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.
us, 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].
e 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.
e 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. e 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]. e 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. e 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.
is paper aims to model the IGBT device using the graphical Matlab software.
e proposed electrical IGBT model is developed using Pspice [15,16] and Saber simulators [17] in the previous contributions. e main contribution of this work is the automatic identification and optimization of the model parameters based on a stochastic algorithm. e 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. e proposed technique is easily implemented in the Matlab environment, and the simulation results proved its high reliability. e 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.

Electrical IGBT Model.
e insulated-gate bipolar transistor (IGBT) is a hybrid transistor, which includes the Metal Oxide Semiconductor Field Effect Transistor (MOS-FET) 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]. e 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]. e static MOSFET model is constituted by  e MOSFET current is controlled by V GE voltage, according to three cases: where K p is the transconductance and θ is the correction factor, which takes into consideration the transverse field in the MOSFET channel For the bipolar part, the static TBJ model is composed by (i) D E and D P : the base-emitter and base-collector junctions (ii) V CE0 : conduction threshold voltage of the baseemitter junction (iii) J PNP : a controlled current source which represents the current gain of the bipolar transistor (iv) β: the bipolar gain

Dynamic IGBT Model.
To model the dynamic IGBT behavior, other elements are added to the static model. ese elements include nonlinear terminal capacitances and the voltage-controlled current source that models the load responsible of the tail current I C at the turn off IGBT state.
Constructors provide the measurement capacitors between two terminals by eliminating one or two capacitors. e input capacitance (C IE ), reverse transfer capacitance (C RE ), and output capacitance (C OE ) allow us to determine the interelectrode capacitance values C GE , C CE , and C GC using the following equations [19]: (3) Figure 1 shows the evolution of C GE , C CE , and C GC capacitances versus V CE . We note that the C GE capacity remains almost constant while the C GC and C CE capacities vary strongly nonlinear following the device operating area.

2
Journal of Electrical and Computer Engineering ese two capacities change values around the same threshold voltage V CET0 .
C GC capacity: this capacity can be constituted of two capacities in series: C OXD oxide capacity that is fixed by the constructor and C GCJ depletion capacity. In Figure 1, the C GC capacity changes value from a certain V CE threshold voltage V CET . We can use the diode model [16] to represent the nonlinearity of the C GC capacity between two values C GCmin and C GCmax (Figure 2). e C GCmax capacity is represented by a voltage-controlled current source, with a gain given by the following expression: C CE capacity: this capacity, as it seen in Figure 1, varies also as nonlinear. It is represented in the same manner as that for the C GC capacity. e coupling coefficient G CE is given by the following expression: e bipolar part presents the responsible charge of the tail current I C , and these charges can be represented only by the PNP base-emitter junction storing charges. ereafter, these charges are modelled by a controlled source (Figure 3), whose expression is as follows: where C BE represents the bipolar part storing charges, with C BE R onDE � t off , where t off is deduced from the device datasheets. e studied IGBT model is shown in Figure 3.

e IGBT Model Parameter Extraction.
is 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 ( Figure 4). Indeed, the diodes states define the different IGBT operation zones.
For the Ohmic zone, all diodes are turned on. e resolution of the equations system described the static model gives the expression of the V CE voltage versus the I C current: For the quasisaturation zone, all the diodes are turned on except the D P diode. e expression of the V CE voltage versus the I C current becomes For the saturation zone, it is characterized by the blocking diode D DC of the MOSFET part and the blocking diode D P of the bipolar part. e expression of the V CE voltage becomes e parameters of the studied model are deduced from the technical datasheet and the equations established in [15]. Figure 2: Model of the nonlinear C GC capacity [15]. e IGBT parameters are recorded in Table 1 for the commercial IGBT (IRGBC20U).

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. erefore, 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. e 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. e 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) e 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 rules A 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. e 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. Figure 5 shows the genetic algorithm diagram. e 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. e selection based on the fitness function allows improving the population [20].
is 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. en, the algorithm manipulates the inputs in order to improve the outputs.   Journal of Electrical and Computer Engineering

Model Parameter Optimization Using GA.
e 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. ey are based on natural selection principle and heredity. is concept is applied to an initial population of potential solutions P(t) 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.
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]. e flowchart of the proposed optimization approach is presented in Figure 7. e fitness function measures the error between manufacturer and simulated data. e 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). e fitness function is given as follows: where I Cdata is the value of the desired current given by the manufacturer and I CM 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 K P , V TH , θ, and the bipolar gain β.
e second program will use equations (7)-(9) with the mathematical equation of the output    Tables 2 and 3, respectively. Figure 8 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. e comparison between the obtained results and the simulation results without optimization is illustrated in Figure 9. e improvement brought by the model in the reproduction of the output and transfer characteristics in the saturation regime is shown. e error between the curves obtained by optimization and the curves given by the manufacturer is low.   Figure 10 illustrates the experimental circuit used to perform the dynamic and static characteristics of the studied IGBT. e electrical IGBT environment is modelled by a wiring inductance L P and a wiring resistance R P . A single parasitic inductance is sufficient to represent the total circuit inductances [24]. e IGBTdriver causes the rise and fall times of the V G impulsion. is interaction between the control section and the device is also included in the test circuit simulation.

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 voltage V GE is maintained at 15 V, and the voltage V CE 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.    Journal of Electrical and Computer Engineering To perform the transfer characteristic I C (V GE ), the V CE is fixed at 100 V while the V GE varies from 5 V to 20 V. Figure 11 presents the manufacturer and simulation results of static characteristics. In the two curves, simulation results show good agreements compared to manufacturer data.

Validation of the Dynamic IGBT Model Behavior.
e switching IGBT characteristics are affected by parasitic impedances, especially parasitic inductances. e I C current decreases and V CE 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 I C and the V CE voltage peak caused by these parasitic inductances increase the duration of the switching operation. en, the switching losses at the turnoff phase are increased. e 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.

Conclusions
is 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. e established method used electrical equations of the model and measurement data provided by the manufacturer to perform the model parameter optimization. e 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.

Data Availability
e data used to support the findings of this study have not been made available because they are confidential.

Conflicts of Interest
e authors declare that they have no conflicts of interest.