Impact of Chaos on MOSFET Thermal Stress and Lifetime

: The reliability of power electronic switching components is of great concern for many researchers. For their usage in many mission profiles, it is crucial for them to perform for the duration of their intended lifetime; however, they can fail because of thermal stress. Thus, it is essential to analyze their thermal performance. Non-linear switching action, bifurcation and chaotic events may occur in DC-DC power converters. Consequently, they show different behaviors when their parameters change. However, this is an opportunity to study these bifurcation phenomena and the existence of chaos, e.g., in boost converters, on their performance as the effects of load variations (mission profiles) on the system’s behavior. These variations generate many non-linear phenomena such as periodic behavior, repeated period-doubling bifurcations and chaos in the MOSFET drain-source current. Thus, we propose, for the first time, an analysis of the influence of chaos on the junction temperature. First of all, this paper provides a step-by-step procedure to establish an electrothermal model of a C2M0080120D MOSFET with integrated power loss. Then, the junction temperature is estimated by computing the power losses and a thermal impedance model of the switch. Additionally, this model is used to investigate the bifurcation and chaotic behavior of the MOSFET junction temperature. The paper contributes by providing a mathematical model to calculate several coefficients based on experimental data and thermal oscillations. Estimation of the number of cycles to failure is given by the Coffin–Manson equation, while temperature cycles are counted using the rainflow counting algorithm. Further, the accumulated damage results are calculated using the Miner’s model. Finally, a comparison is made between the damage accumulated during different mission profiles: significant degradation of the MOSFET’s lifetime is pointed out for chaotic currents compared to periodic ones.


Introduction
Today, power switching semiconductor devices such as IGBTs and MOSFETs are widely used in aerospace, transport and industrial applications.The trend towards high power density is accompanied by increasingly severe operational constraints in terms of current, voltage and temperature.Indeed, it is important to maintain the functional performance of the components in their new environment, but above all to ensure their operational safety to avoid any power device failure that could impact the entire system [1].
Among the components, power semiconductor devices are the most fragile in converters and inverters and are strongly affected by the mission profiles [2].Thus, one approach is to adopt a realistic mission profile in order to obtain a reliable assessment of the MOSFET's lifetime.There are daily solar irradiance profiles for photovoltaic systems [3,4], reference speed [5,6] and reference torque [5] profiles for drive motor applications and reference voltage [7] and load mission profiles [8] for electric bus applications.
Due to the strong non-linear modeling of switching converters [9], these converters can exhibit complex dynamical behaviors such as bifurcation and chaos with variation of a bifurcation parameter [10].This bifurcation parameter can be the load [11], the reference voltage [12], the reference current [13,14], the inductance [15] of the converter, or a control parameter [16,17].There is a bifurcation route towards chaos when increasing or decreasing these parameters, forcing the MOSFET to operate under harsh conditions.In this paper, the mission profiles are the different values of the load bifurcation parameter: load variation frequently leads to a loss of stability if the duty cycle is outside a certain range.
Generally, a MOSFET is considered an ideal switch in fast simulation in order to analyze the efficiency and reliability of power electronic systems, as in [18][19][20][21][22]. On the other hand, it is difficult to calculate the MOSFET's power losses (which include the conduction, switching losses, diode conduction and reverse recovery losses [23]) and to take into account the effect of the junction temperature on the device parameters.In the literature, various electrothermal modeling strategies have been reported.The lookup table method is applied [24] to approximate switching as well as on-state losses of the IGBT.In [25], a more complete electrothermal model is proposed, which takes into account the junction temperature and also the electrical load as input variables.The model of [26] considers the interactions between the junction temperature and the voltages and currents of the transistor.In [27], the power losses are averaged every several switching cycles using a high-speed electrothermal model.The authors of [28,29] determine a realistic model of a DC-DC converter wherein the parameters are identified based on the component datasheets, the printed circuit board (PCB) layout and the length and diameter of the cables.Some parasitic inductances and capacitances result from the geometric structure of the circuit (e.g., the printed circuit and the cables used).
To analyze the MOSFET's thermal behavior [30], power losses are used to deduce the junction temperature.The heat transfer model can use the RC Foster [31,32] or Cauer [33] networks, the finite element method or the finite difference method.In [34], a third-order Foster circuit is used to fit the impedance curve of a MOSFET given by the datasheet.In [35], 450 different electrothermal (Foster) variants are proposed to simulate an electric vehicle inverter.In [36], the junction temperature is estimated using the finite element method to reduce the duration of accelerated aging tests.With on-state voltage measurement, ref. [37] uses a Kalman filter for the junction temperature estimation of a full-bridge inverter.In [38], the resistances of the linear Cauer network are replaced by controlled voltage sources.In this way, a more realistic analysis is realized that considers the thermal coupling between power devices containing non-linear convection terms.
A real-time monitoring strategy of the junction temperature is proposed by [39].The junction temperature distribution is not uniform, and control solutions depend on the load profile: the controller is adjusted to ensure the desired performance [40].Active thermal control is used by [41] to adjust the thermal stress of the power device.The authors of [42] improve the thermal management in an IGBT's structure, thus reducing the leakage current.
According to [43], changes in the power device junction temperature affect their reliability and lifetime estimations.From the thermal and aging points of view, the R ds−on drain-source MOSFET resistance is influenced by the junction temperature, as in [44].The two most common failures for power electronics devices are die-attach solder fatigue [45] and bond-wire damage caused by temperature cycling.Bond-wire damage is serious degradation of a MOSFET, while solder fatigue does not destroy the device but affects its durability.Therefore, mean junction temperature and junction temperature variation are both required for lifetime estimation.
In order to evaluate the lifetime of a MOSFET, the Coffin-Manson law is the most used model to evaluate device failure cycles.For example, in [46][47][48], the lifetime of an inverter [49] and of a multi-level converter is estimated with the rainflow counting algorithm with the mean junction temperature and the junction temperature variation as inputs.
To illustrate our method, we chose the C2M0080120D MOSFET [50] to study the influence of the junction temperature on its lifetime when the boost converter loses its stability.Therefore, a bifurcation diagram is defined to understand and analyze the stability of boost converters under different operating conditions.The load (representing the mission profile) generates many non-linear phenomena as periodic behavior, repeated period-doubling bifurcations and chaos of the MOSFET drain-source current, impacting the junction temperature.Then, a step-by-step procedure for establishing the electrothermal model of the MOSFET with integrated power loss is proposed in this paper.We highlight not only a reduction in the stable range of the junction temperature variation but also an extension of the chaos range.Additionally, we propose a mathematical model to calculate several coefficients based on experimental data and thermal swings.The estimation of the number of cycles to failure is given by the Coffin-Manson equation, while temperature cycles are counted using the rainflow counting algorithm.Accumulated damage results from the Miner's model.This paper presents for the first time the significant impact of the overall magnitude of the current ripple for the chaos behavior, which leads to significant thermal stress.Finally, a comparison is made between the damage accumulation by different mission profiles: an important degradation of the MOSFET's lifetime is pointed out for a chaotic current compared to a periodic one.
The main aspects of this paper are as follows: Section 2 describes the boost converter system.In Section 3, bifurcation analysis is carried out for different mission profiles.Then, in Section 4, the electrical model is detailed, followed by the power loss calculations for the MOSFET and diodes in Section 5. Thermal model simulation results are presented whereby a cascade of bifurcations doubling the stable junction temperature period is observed in Section 6. Section 7 details the lifetime model that estimates the damage to a C2M0080120D MOSFET.We end with a conclusion in Section 8.

Boost Converter
The schematic diagram of the boost converter circuit that uses a current loop is shown in Figure 1.The target power device is a silicon carbide (SiC) C2M0080120D MOSFET transistor.The circuit has two states determined by the state of the transistor.The MOSFET is controlled by an R − S latch.A clock pulse of period T sets the latch at the beginning of every cycle, turning on the MOSFET: the inductance current i increases.I re f denotes the reference peak inductor current i L .When i exceeds the value of I re f , the latch is reset and the switch is turned off.The switch again turns on at the beginning of the next clock pulse.The capacitor voltage v and the inductor current i are the state variables of the circuit.The boost converter operating in continuous mode is described by two sets of linear differential equations as:

•
MOSFET switch on: ẋ(t) = A on x(t) + B on V in (t) Equations ( 1) and (2) solve for the the currents and voltages across the components.An electrothermal model is developed based on the power module configuration.For the electrical model of the boost converter, Figure 1 depicts the conduction and switching power loss blocs of the MOSFET and diode generated for a selected mission profile.Then, losses are injected into the thermal model for the calculation of the device junction temperature T j .The estimated T j provides feedback for power loss computations, which it affects at each next time step.The simulations are made with the SimScape toolbox of MATLAB R2020b Simulink and include the boost converter with its control strategy (peak current control); the power losses and thermal model are MATLAB Simulink subsystems.

Non-Linear Behavior and Bifurcation Analysis of the Boost Converter
In this section, we present the dynamical behavior and bifurcation diagram of the boost converter with the current mode control.Figures 2-4 show the simulation results of inductance current obtained with the load resistance R as the bifurcation parameter.R is varied over the range 5 Ω-20 Ω without a thermal effect.
A cascade of period-doubling bifurcations of stable cycles is observed in the system.Figure 2a presents a stable limit cycle (period-one T response).With a further increase in R parameter values (from R = 6 Ω to 8 Ω), a stable limit cycle of a double period 2T is observed in Figure 2b.A stable cycle of period four emerges at R = 13.1 Ω, and at R = 13.8Ω, there is a stable cycle of period eight (Figure 3).For higher values of the bifurcation parameter R, the system reveals a sequence of stable cycles of a subharmonic cascade of bifurcations.For R = 15.2Ω, the system is attracted to the period-six T orbit, as in Figure 4a.If the value of R is further increased, the system is stroked by a border collision before it goes into chaos (Figure 4b).Obviously, in a chaotic state, the system operates with large values of inductance current i L .
A bifurcation diagram represents all the plots of the steady-state orbit as a function of the bifurcation parameters.We vary the load resistance R parameter and plot hundreds of consecutive values of the inductance current i L at each clock period T. The bifurcation diagram is presented in Figure 5: the period-doubling bifurcation occurs at R = 7.32 Ω and changes the steady-state behavior of the system.Along with the increase in R, the converter undergoes several period-doubling bifurcations and eventually evolves into a chaotic state.For different values of R, one can easily identify the boost converter behavior: the system is periodic for small values of R while it is chaotic for large values of R. It is essential to study the occurrence of chaos due to the variation of system parameters to predict non-linear phenomena and to size components for correct operation.These six values of R (6 Ω, 8 Ω, 13.1 Ω, 13.8 Ω, 15.2 Ω and 20 Ω) are used to generate the load profiles.

Electrical MOSFET Model
Power semiconductors are non-ideal transistors because of their parasitics, such as source inductance L s , drain inductance L d , and gate inductance L g and parasitic capacitances C gs , C gd and C ds , as shown in Figure 6.Their real electric model takes into consideration the source resistance R s , drain resistance R d , gate resistance R g and drain-source resistance R ds .
where U th is the threshold voltage, k is the transistor gain, and λ is the channel length modulation.Based on the static current-voltage characteristics shown in Figure 7a, the following polynomial fitting functions are used to obtain fitted simulation curves of I ds -U ds at two different junction temperatures T j : The coefficients of the I ds -U ds curves are fitted using a MATLAB curve-fit tool.Figure 7b displays the static current-voltage characteristics I ds -U gs described by the following polynomial fitting functions: In order to verify the accuracy of the MOSFET electrical model, these static characteristics obtained by simulation are compared with the ones from the datasheet; the characteristics are reported together in Figure 7a,b, respectively, for 25 • C and 150 • C. Indeed, these curves are very close to each other over the full operating range.

Power Loss Computation
The MOSFET produces power losses, which can be distinguished as conduction loss and switching loss.Figure 1 presents also the power losses for the diode.Addition of all losses computes the total power loss.

MOSFET Power Loss Computation
The conduction power loss at each time step can be calculated as the product of the drain-source current I ds through the MOSFET during conduction and the drain-source voltage U ds across it, as shown in: Using the static current-voltage characteristics I ds -U gs (Figure 7a), the fittype and fit MATLAB functions determine the voltage U ds at two different temperatures, 25 • C and 150 • C, as in the following equations: For many manufacturer datasheets, the dependence of switching energy losses on temperature, current or voltage is passed over.For example, the C2M0080120D datasheet gives the energy switching loss for two values of drain-source voltage (600 V and 800 V) at 20 • C and for a large variation of drain-source currents.The manufacturer datasheet also provides the energy switching loss for a large variation of the junction temperature at a drain-source current of 20 A and drain-source voltage of 800 V.The lookup table interpolates these know values in order to generate the switching energy losses at different data points as a function of the mission profiles.
The lookup table method is applied to approximate conduction as well as switching losses of the MOSFET and diode.A 2-D lookup table having vectors of temperature T j and current I ds as breakpoint inputs interpolates the voltage U ds at any I ds and T j values.Based on Equations ( 9) and ( 10), the 2-D lookup table of U ds is represented as Figure 8a.Then, the instantaneous conduction power loss P Mcond of the MOSFET is calculated as the product between I ds and U ds .The C2M0080120D MOSFET datasheet presents two switching loss energies: E sw−on (µJ) and E sw−o f f (µJ)-the device turn-on and turn-off energies, respectively, for junction temperatures of both T j = 25 • C and U ds = 800 V, as shown in Figure 9a.These two energies are calculated using the rise and fall times of the drain-source current and voltage.
The fitting functions of the datasheet points represented in this figure are The datasheet also presents the variation of the switching loss energies E sw−on and E sw−o f f as a function of the junction temperature T j (Figure 9b).For I ds = 20 A and U ds = 800 V, the polynomial fitting functions give a linear variation as: A comparison of the simulation results (based on Equations ( 11) and ( 12)) and the datasheet points is shown in Figure 9a for varying drain-source current I ds .Figure 9b presents the same switching loss energies with T j variation.A good match is found for both turn-on and turn-off energy losses for the case of datasheet points and simulation results (based on Equation ( 13)) with T j variation.E sw−on and E sw−o f f depend on three variable functions of I ds , U ds and T j .Based on Equations ( 11) and ( 13), a 3-D lookup table having the temperature T j and the current I ds as inputs (MATLAB breakpoints of 3-D lookup table) and U ds as the output interpolates the switching loss energies E sw−on at any I ds , U ds and T j values for the state MOSFET turn-on (Figure 8b).
Another 3-D lookup table (Figure 8c) is used to predict the switching loss energies E sw−o f f with T j , I ds and U ds as inputs based on Equations ( 12) and ( 13).E sw−on and E sw−o f f at 25 • C and 150 • C, presented in Figure 10, are used to estimate switching losses.The instantaneous switching power loss P Msw of the MOSFET is then calculated as the product of the loss energies and the switching frequency.

Diode Power Loss Computation
Calculating the conduction loss of a body diode is similar to its calculation for a MOSFET.The following equations present the diode conduction loss as: where U f represents the forward saturation voltage, and I f represents the forward saturation current.These two variables can be extracted from the datasheet and are indicated as datasheet points in Figure 11 for two different temperatures (T j = 25 • C and T j = 150 • C). Figure 11a shows a comparison of datasheet points and fitted curves obtained using the following polynomial fitting functions: A 2-D lookup table (represented in Figure 11) has the temperature T j and the current I f as breakpoint inputs and U f as the output: it interpolates the voltage U f for any I f and T j values.Then, the instantaneous MOSFET conduction power loss P Dcond is calculated as the product of I f and U f .
The diode contributes to switching energy E rr as a result of reverse recovery charge Q rr during turn-off of the MOSFET.In [23], the analytical model of the diode switching energy depends on the peak reserve recovery current I rr , the junction T j temperature and the diode forward saturation voltage U f .Furthermore, the slope of the diode turn-off current depends on the switching rates.The peak reverse recovery current I rr is a function of the turn-off switching rate (di f /dt), the reverse recovery charge Q rr and the snappiness factor S as follows: The snappiness factor S is calculated with specific information from the datasheet: Q rr (25 • C) = 152 nC and di f /dt = 1950 A/µs, giving S = 4.928.S is considered constant for different operating conditions.In order to calculate Q rr at T j = 125 • C, the following relation is used: where di f /dt = 1090 A/µs.The reverse recovery charge is 271.93 nC.Equation ( 19) is used to simulate the effect of the temperature T j on the reverse recovery charge Q rr : Then, the energy switching loss E rr is defined as the product of Q rr and U f .Figure 12a represents E rr at 25 • C and 150 • C. A 3-D lookup table with breakpoint inputs T j , I f and U f calculates by interpolation the switching loss energies E rr , as shown in Figure 12b.The instantaneous diode switching power loss P rr is then calculated as the product of the loss energies and the switching frequency.Finally, the total power loss P out is the sum of all power losses.

Thermal Model
The thermal behavior of a MOSFET can be predicted with a resistor and capacitor thermal network.There are two conventional thermal networks that can equivalently represent the thermal performance of a MOSFET: Foster and Cauer.Due to its simplicity, almost all semiconductors use the Foster network.The thermal model of a SiC C2M0080120D MOSFET represented as an RC network is derived from the thermal impedance.The total power losses P out are passed through the Foster network and result in the device junction temperature.The junction temperature T j is estimated by: where Z th−jc (t) is the thermal impedance, and T a (t) is the ambient temperature.According to Equation (20), computation of the junction temperature T j requires knowledge of the thermal impedance Z th−jc (t).Z th−jc (t) is graphically provided in the datasheet of the C2M0080120D MOSFET, as shown in Figure 13.The Foster thermal network is achieved by implementing a curve-fitting identification.The parameters of the Foster model are summarized in Table 1: it contains four RC branches to maintain consistency with the original data from the manufacturers.Different load values represent a typical load mission profile of the electrical boost converter.After evaluation of the corresponding power losses, the thermal model of the MOSFET produces the thermal profile shown in Figures 14-16.Figure 14a presents a stable limit cycle of junction temperature T j for a load profile R = 6 Ω.The T j time waveform is a 1T limit cycle for inductance current i L (Figure 2a) when the thermal conditions are passed over.With a further increase in parameter R (from 6 Ω to R = 8 Ω), a stable limit cycle of a double period 2T is observed in Figure 12b, as for i L in Figure 2b.On the other hand, if the load increases to R = 13.1 Ω, the junction temperature T j does not follow the same behavior as i L (stable period-4T operation as in Figure 3a), instead keeping a stable limit cycle with a 2T period.This demonstrates that the temperature influences the behavior of the converter.At R = 13.8Ω, the stable cycle of period 8T (Figure 3b) is divided by two.Indeed, T j is attracted to the period 4-T orbit, as can be seen in Figure 15b.If the value of R is further increased, the system is stroked by a border collision before it goes into chaos (Figure 16a,b      Therefore, the temperature has an influence on the behavior of the converter: the period-doubling bifurcation is restrictive, causing chaos to appear earlier, i.e., for lower values of R. Additionally, for the six mission profiles, the junction temperature variations increase from an extremely small ∆T j (less than 0.5 • C) to a large ∆T j (greater than 18 • C).It is the same for the mean junction temperature: T jm increases from 88.17 • C to 127.39 • C.

Lifetime Analytical Model
An accelerated power cycling test has been designed and executed to collect failure data.For this purpose, the jumps in junction temperature ∆T j and the mean junction temperature T jm are the main stress factors of the power device.It can be noted that device performance degradation differs from one device to another due to small discrepancies in device manufacturing parameters.In order to determine the qualitative behavior of the C2M0080120D, three MOSFETs are exposed to multitude power cycles under the same conditions.The experimental test represents a long cycle time: the MOSFET is turned on for 4 s and turned off for 12 s.The constant current value I ds is intended to maintain the temperature at 125 • C. The datasheet characteristics I ds -U ds and I ds -U gs (Figure 7) present variations between different devices of the same type.Therefore, variability exists between the mean and swing junction temperatures of the tested MOSFETs.The temperatures obtained during three experimental tests to determine the heating behaviors were: T jm1 = 127 • C and ∆T j1 = 16 • C, T jm2 = 126.5 • C and ∆T j2 = 14.5 • C and T jm3 = 114.2• C and ∆T j3 = 12.5 • C. According to the test conditions, device failure was recorded after N f 1 = 8640 cycles, N f 2 = 12,270 cycles and N f 3 = 25,400 cycles.
These results are used to derive the experimental coefficients A and δ and also the activation energy E a .The Coffin-Manson equation is used to estimate the device failure cycles.The lifetime model of the power device for the first test can be expressed as where N f 1 is the number of cycles to failure for a certain thermal stress, and k is the Boltzmann constant (8.617 × 10 −5 eV/K).For the second and third tests, the Coffin-Manson models are The coefficients A, δ and E a are calculated from Equations ( 21) and (22).Applying the natural logarithm, these exponential functions ( 21) and ( 22) can be expressed in linear forms.Eliminating the unknown constant A, we can write and Then, log Consequently, the development of δ can be represented by Hence, the thermal activation energy can be calculated from The last coefficient can be written as Using the three experimental test datasets and according to relations ( 26)-(28), the following numerical values are obtained: δ = −4.4887,E a = 0.0667 and A = 2.8823 × 10 +8 .Hence, the calculated values are: for the first profile ∆T j1 = 0.52 • C, T jm1 = 88.17In order to predict the lifetime of a power MOSFET, a thermal cycle counting algorithm is required, especially when the MOSFET operates with a variable load.Due its low error compared to other counting methods, rainflow counting is the most widely adopted method.A rainflow counting algorithm translates complicated and irregular stress profiles (time series of the junction temperature) into a set of organized cycles to facilitate the lifetime prediction for MOSFETs, especially when the converter operates with varying loads.
The basic sequence of the rainflow algorithm includes the four following steps.Firstly, the peaks and valleys of the junction temperature are extracted in order to construct linear connections.For the second step, the junction temperature undergoes a rotation of the time axis by 90 • such that the plot has the shape of a pagoda, as in Figure 17.This figure plots the virtual water traces (the red and blue lines) with the origins on each peak and valley.Then, the end of the trace occurs by merging it with a virtual water trace from a higher peak or lower valley or by reaching a peak for which its opposite peak is higher and the valley is lower.The final step is to couple half cycles to full cycles.The number of cycles to failure for each mission profile is extracted by the rainflow algorithm and is presented in Figures 18-20.
For simple periodic mission profiles as shown in Figure 18a, the counting algorithm is not very useful since the cycles are easily determined.However, in cases of aperiodic mission profiles, the use of counting algorithms such as rainflow is very efficient for extracting the cycles.The complete thermal profiles contain a large number of thermal cycles from the first mission profile (Figure 18), on the order of 8.5 × 10 +5 cycles, whereas the last mission profile (Figure 20) has a relatively small number of thermal cycles (1956 cycles).Another important factor in lifetime estimation is the accumulated stress Q of the MOSFET, which evaluates the accumulated stress after a number of thermal cycles.Q is appraised using the Miner's rule.It is evaluated as: where N c represents the number of cycles calculated using the rainflow counting algorithm, and N f is the number of cycles to failure at a particular stress condition using Equation (21).
With the help of the Miner's rule, degradation of a semiconductor and its remaining lifetime can be extracted, as in Table 2.The accumulated damage is insignificant for lower values of load (6 Ω and 8 Ω), indicating a small impact due to the mission profile on the MOSFET's lifetime.However, only a small damage increase is observed when the load is raised to 13.1 Ω.The mission profile begins to have a strong impact on the reliability assessment of the MOSFET for a load bigger than 13.1 Ω: for example, at 13.8 Ω, the junction temperature behavior overbalances to a stable period-4T operation.The repercussions of the chaotic mission profile are significant for the thermal stress and the reliability performance of a MOSFET: the accumulated fatigue is 10% to 65% more than the results from the reduced profile.

Conclusions
In this paper, the impacts of mission profiles on the thermal stress of a MOSFET and on its lifetime have been investigated.Several mission profiles from stable 1T-period to chaotic behavior have been considered.The results indicate that chaotic mission profiles have a strong impact on the swing current amplitude of the MOSFET and on the junction temperature (especially for large values of load).Taking into account the impact of large load values into the bifurcation diagram and the junction temperature, a reduction in the stable range of the load variation was highlighted together with an extension of the chaos range.We also proposed a mathematical model to calculate several coefficients based on experimental data and thermal swings and then used it to calculate the failure cycles.This paper presents the significant impact of the overall magnitude of the current ripple for the chaos behavior, leading to large thermal stress.Finally, a comparison is made between the amounts of damage accumulation for different mission profiles: we point out a non-negligible degradation of the MOSFET's lifetime for chaotic (reducing it by half) compared to periodic current behavior.

Figure 1 .
Figure 1.Block diagram of electrical and thermal model of boost converter: L = 12 mH, C = 47 µF, I re f = 20 A, T = 100 µs, and V in = 48 V.

Figure 5 .
Figure 5. Bifurcation diagram of the boost converter for inductance current i L versus R without thermal effect.

Figure 6 .
Figure 6.MOSFET electrical model.Some MOSFET parameters have been extracted from electrical characteristics (from static current-voltage characteristics I ds -U ds and I ds -U gs ) and the reverse diode from the C2M0080120D MOSFET datasheet (the point curves of Figure 7).All these parameters (resistances, inductances and capacitances) are required by the SimScape MOSFET model.Herein, the MOSFET model using MATLAB Simulink is based on the Shichman-Hodges model.The Shichman-Hodges model takes channel length modulation into account and uses a controlled current source to describe the drain-source current I ds .In order to reproduce these two characteristics, the Shichman-Hodges model describes the three

Figure 7 .
Figure 7. Fitting curves of the static current-voltage characteristics: (a) I ds -U ds at 25 • C and 150 • C and U gs = 20 V; (b) I ds -U gs at 25 • C and 150 • C and U ds = 20 V.

Figure 8 .
Figure 8.(a) 2-D lookup table of the output U ds (I ds , T j ); (b) 3-D lookup table of E sw−on (I ds , U ds , T j ); (c) 3-D lookup table of E sw−o f f (I ds , U ds , T j ).

Figure 9 .
Figure 9. MOSFET (C2M0080120D) energy power losses: (a) varying the drain-source current I ds , T j = 25 • C and U ds = 800 V; (b) varying the junction temperature T j , I ds = 20 A and U ds = 800 V.

Figure 10 .
Figure 10.(a) E sw−on varying the drain-source current I ds and the drain-source voltage U ds at 25 • C and 150 • C; (b) E sw−o f f varying the drain-source current I ds and the drain-source voltage U ds at 25 • C and 150 • C.

Figure 11 .
Figure 11.(a) Body diode characteristics at T j = 25 • C and T j = 150 • C; (b) 2-D lookup table of the output U f .

Figure 12 .
Figure 12.(a) E rr when varying the forward saturation current I f and the forward saturation voltage U f at 25 • C and 150 • C. (b) 3-D lookup table of E rr (U f , I f , T j ).

Figure 13 .
Figure 13.Comparison between the transient thermal impedance Z th−jc curves obtained via simulation and from the datasheet.

Figure 15 .
Figure 15.Junction temperature T j : (a) Stable period-2T operation for a load profile R = 13.1 Ω; (b) Stable period-4T operation for a load profile R = 13.8Ω.

Figure 16 .
Figure 16.Junction temperature T j : (a) Chaotic behavior for a load profile R = 15.2Ω; (b) Chaotic behavior for a load profile R = 20 Ω.

Figure 17 .Figure 18 .Figure 19 .Figure 20 .
Figure 17.Application of the rainflow counting algorithm for the cycle extraction of a complex junction temperature.

Table 2 .
Accumulated fatigue results for the six load mission profiles.