Next Article in Journal / Special Issue
Multi-Electrode Architecture Modeling and Optimization for Homogeneous Electroporation of Large Volumes of Tissue
Previous Article in Journal
Carbon Nanotube Films for Energy Applications
Previous Article in Special Issue
EMC Component Modeling and System-Level Simulations of Power Converters: AC Motor Drives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of a Gate Distribution Layout to Compensate the Current Imbalance Generated by the 3D Geometry of a Railway Inverter

by
Andressa Nakahata-Medrado
1,2,
Jean-Luc Schanen
1,*,
Jean-Michel Guichon
1,
Pierre-Olivier Jeannin
1,
Emmanuel Batista
2 and
Guillaume Desportes
2
1
G2ELab, University Grenoble Alpes, 21 Av des Martyrs, 38000 Grenoble, France
2
Alstom Transport, 21 rue Dr Guinier, 65600 Semeac, France
*
Author to whom correspondence should be addressed.
Energies 2021, 14(7), 1891; https://doi.org/10.3390/en14071891
Submission received: 18 February 2021 / Revised: 23 March 2021 / Accepted: 25 March 2021 / Published: 29 March 2021
(This article belongs to the Special Issue Electromagnetic Modeling in Power Electronics)

Abstract

:
The impact of the stray inductances originated from interconnects in power electronics becomes crucial with the next generation of SiC devices. This paper shows that the existing layout of a railway inverter, operating with Si IGBTs already exhibits a dynamic current imbalance between paralleled modules. This will not allow using this geometry with SiC MOSFETs. A complete investigation of the electromagnetic origin of this issue has been performed. A generic circuit model has been proposed to establish a cabling rule to design a Gate Distribution Printed Circuit Board (PCB) in such a way that it compensates the power dissymmetry. An optimization strategy has been used to obtain a new geometry of this PCB, which has been validated with a time domain simulation.

1. Introduction

The European project “Shift2rail” [1] is investigating the advantages of the SiC wide bandgap devices in railway applications. The potential benefits of using SiC MOSFETs in future converters are explored, but also the mandatory impacts on the 3D layout design. With the high speed commutation of these new switches, there is a strong need of properly designing a geometrical arrangement of the various parts of a converter, in order to prevent all consequences of the stray elements originated from the cabling. Voltage overshoots, ringing, and power drive interactions are major concerns of power electronics designers in the “wide bandgap era” [2,3,4].
Among the issues arising from high switching speed in high current applications, the current division among paralleled power modules is a key topic, which has to be addressed carefully to avoid any derating of the devices [5].
Several parameters may generate a bad current division among power modules [6]:
  • difference between device’s parameters due to industrial processes or different operating temperatures,
  • the power layout which provides various current paths among the paralleled devices and therefore different loop inductances, leading to a different current commutation speed (dI/dt),
  • the gate circuit layout, including the common emitter and the mutual coupling between gate and power circuits: both effects modify the gate-to-emitter voltage, and consequently the dI/dt.
All these phenomena have been investigated for many years for Si IGBTs, and more recently for SiC MOSFETs. The impact of power layout dissymmetry has been reported in [7,8,9,10]. It was shown that an electrical symmetry (i.e., equivalent loop inductance seen from each devices) can be obtained even with non-symmetrical geometries, taking into account the effect of mutual inductances. However, this necessitates many degrees of freedom for the 3D placement of the power modules, which are not always available in a railway inverter. Indeed, compact solutions are requested, and the cold plate of the cooling system often imposes the modules arrangement.
Another very important origin of the bad dynamic current division is due to the gate driver. Many papers have addressed the issue of the common inductance between power and gate circuits on the emitter (or source) pin [11,12,13,14]. Of course, this phenomena has mainly been studied for the internal layout of the power module, and has more impact on the current division among the dies than the difference between several modules. However, these papers also described the impact of the mutual coupling between gate and power circuits, which is another origin of power-drive interaction. In other words, reducing the common emitter impedance is not sufficient to avoid the feedback effect from power to gate during transients, and the impact of the magnetic couplings has to be considered. This has been done in [15], where a metric has been defined in order for all the dies inside the power module to be submitted to the same feedback voltage from the power circuit. Therefore, using mutual couplings, and despite different values of common emitter inductances, all gate to emitter voltages during transients were identical, leading to the same dI/dt. Defining a gate circuit layout according to this metric was obtained, using an optimization process, for a fixed power layout.
Nevertheless, in a general situation, the current imbalance among power devices during transients combines both effects of power layout and power-drive interaction. To elucidate these joint effects, it is necessary to add the device’s behavior. This was done in [16] to account for electro-thermal effects, the device’s parameters being modified with temperature.
The present paper will propose a full investigation of the impact of all inductive stray effects generated by both power and gate layout. It will be shown that a dynamic current imbalance generated by the power layout can be compensated by properly designing a gate circuit layout. The electrical equations developed in this paper are more or less similar to the recent work reported in [17] and will be used in an optimization process to obtain an improved current division, despite the current imbalance caused by the power layout.
The approach to get the circuit equations must start with a clear link between the 3D geometry and the electrical representation. Partial Element Equivalent Circuit (PEEC) [18] is today the classical modeling approach to reach this goal. It allows for the representation of any conductor with a partial inductance, and for mutual couplings with the rest of the circuit. However, if the transfer from a simple interconnection linking two points to its electrical model is straightforward, multiple port connections, as busbars, must be clarified. This paper will present a generic electrical representation for this kind of multiport conductors, which is implemented in many PEEC software but not reported in the literature so far.
To support the proposed method, a 1700 V 1200 A railway inverter, ONIX 671, will be used as an example. Even if this inverter is not based on SiC devices, it exhibits some current imbalance, whose origin will be investigated thanks to simulation (Section 2). In Section 3, the generic circuit model derived from the 3D layout will be presented. Equations got from the model will be used to obtain a cabling rule to guarantee balanced currents during transients, despite a non-symmetric power layout. This cabling rule will be the objective function of the optimization process carried out in Section 4.

2. Inverter Description

The 1700 V 1200 A studied inverter is a 3 phases + 1 brake chopper topology. The electrical representation of the inverter is given in Figure 1, and the 3D layout in Figure 2. Several parts can be identified:
  • An Electrolytic Capacitor Bank made of 54 devices (18 groups of three series capacitors), connected to the two HT+ and the two HT− terminals (DC bus);
  • two identical subparts, each of them realizing two legs. One sub part is composed of six power modules (three paralleled devices for top and bottom switches), eight decoupling capacitors, and three busbars for linking the aforementioned elements: one power busbar (DC bus) and two phase busbars for connection to the output phases.
  • Each subpart contains drivers for controlling the IGBTs. A specific Printed Cicuit Board (PCB) carries the driver signal to the three paralleled devices for the top and bottom switch of each phase. Two Gate Distribution PCBs are thus used for each subpart.
  • A cold plate is supporting all 12 power modules (dry natural cooling system).
Figure 2 shows the 3D representation of the inverter. The Electrolytic Capacitor Bank associates three groups of capacitors realized with simple copper bars and is linked to the two inverter subparts through a dedicated busbar. The HT+ and HT− connections are therefore duplicated. The IGBT module is also displayed in Figure 2. It offers separate Emitter and Collector pins for top and bottom switches (E1, C1 and E2, C2), the series connection being made by the phase busbar. A Kelvin connection is of course provided for each emitter of each IGBT (Ec1 and Ec2). Gate pins G1 and G2 as well as separate collectors Cc1 and Cc2 for control purposes are also indicated in the figure. On the bottom of Figure 2, one can see the inverter power layout, which is split into two parts, each of them associating two inverter legs. For each subpart, six power modules are fixed on the cold plate and interconnected with a power busbar in the center and two phase busbars on both sides, which also provide the connection between all top emitters and bottom collectors of the three paralleled IGBTs. The power busbar associates the power module and eight decoupling capacitors, and provides links to the DC capacitor tank. The gate drivers are on the top, and a specific PCB is used for interfacing them with the three paralleled IGBTs. These four Gate Distribution PCBs are localized directly over the phase busbars.
The cold plate is not much larger than the 12 power modules. Therefore, there are almost no degrees of freedom for their geometrical implementation. The DC power busbar, linking the decoupling capacitors Cdec to the power modules, has to provide both small loop inductance and as much symmetry as possible for ensuring an equal current division between power modules. As mentioned before, the power modules are not providing the series association of the two IGBTs of the leg, and the electrical link between the Emitter of switch 1 and the Collector of switch 2 is performed through a part of the phase busbar. The switching cell is thus composed of Decoupling Capacitors, Power Busbar, and a part of the Phase Busbar.
It is worth noting that all phase busbars exhibit a slit, visible on Figure 3. This one has been inserted for managing the static current division (DC resistance) between the three paralleled IGBTs. Indeed, the current path from the power modules to the phase output will have almost the same length, whereas without the slit the module close to the power connections would see a much shorter path than the most distant one.
In the experiment, the operations of this converter were not fully satisfying regarding the current division among the three paralleled modules. A huge experimental investigation was carried out to check if this current imbalance was due to the devices’ properties, but even when carefully sorting IGBTs with the same characteristics, the issue remained.
The 3D layout is thus at the origin of the bad current distribution among the three switches. Moving to SiC MOSFET with this kind of layout is not possible because the influence of the cabling will be much more important with a higher commutation speed. Understanding the root cause of the issue is therefore mandatory. For this purpose, an investigation using a simulation has been carried out, using Ansoft Q3D [19] and the associated time domain simulator Simplorer.
For this purpose, it was decided to focus on the bottom switch of phase W, composed of 3 paralleled IGBTs. The power layout is composed of the DC busbars and the phase W busbar. IGBT power modules have been simplified and represented by simple lumped inductances, to avoid overloading the 3D simulation model. This simplification is not a real issue, since all modules are identical: this is not the origin of the current dissymmetry. The eight decoupling capacitors have been identified through impedance measurement (simple esl–esr model [20]) and their equivalent circuit is connected to the appropriate ports of the 3D PEEC model. The DC capacitor bank has also be simplified: from Figure 2, it can be seen that the parallel association of 18 groups of three series capacitors is realized with three blocks of 6 * 3 series components. An equivalent model of each block has been identified (again, equivalent series inductance and resistance), and only the DC capacitor busbar (Figure 2) has been modeled in Ansoft Q3D. The power busbar and phase W busbar are also included in the 3D description. It is worth noting that to save memory space, the DC capacitor busbar has been modeled in a separate file, since its mutual coupling with the rest of the power structure is negligible.
The gate driver is replaced by a perfect voltage source, common to each of the three paralleled IGBTs. The link between this voltage source and each of the three devices is made by the Gate Distribution PCB (Figure 4 represents the part of this PCB corresponding to the bottom switches). It is a pure passive PCB, connecting the drive signal to each IGBT, and providing external gate resistances (some additional feedback and protection circuits are not studied here). It was included in the full 3D model, together with power and phase W busbars.
An illustration of the electrical simulation with all aforementioned blocks is provided in Figure 5.
Figure 6 shows the simulation results for the three IGBTs. It is obvious in this figure that IGBT 3 (far from the phase output and close to the driver input) is slower than the two other ones, resulting in higher switching losses at turn off. IGBT 2 (in the middle) is the quickest one, and consequently takes the major part of the diode recovery current.
To understand the root cause of the imbalance and determine if was due to the power layout or to the power-drive interaction, a second simulation was done, removing the Gate Distribution PCB from the 3D representation and directly controlling the IGBTs by the ideal driver voltage source. Therefore, only the power layout was impacting the current division. Simulation results are given in Figure 7. Again, current division is not perfect, but this time the slowest device is not IGBT3, as in Figure 6, but IGBT 1. Therefore, the current imbalance is necessarily affected by both effects, power and gate circuits.
A last simulation was also carried out to confirm the impact of mutual couplings in the power-drive interaction. The Gate Distribution PCB was described in a separate file in Q3D, therefore keeping the effects of its impedance, but not considering the mutual couplings with the power layout. The results of Figure 8 show that IGBT1 is still the slowest, as with the perfect gate circuit, which clearly shows that the most impacting effect is due to this coupling between Power and Gate Distribution PCB. This is not surprising considering the position of this PCB, directly on the top of the phase busbar, which is part of the switching cell, so carrying high dI/dt.
From these simulations, it can be concluded that:
  • the power layout generates a current imbalance, IGBT 3 being slower than the two others
  • the couplings between this PCB and the power part actually modify the current division, since taking them into account leads to a real change in the slowest device (IGBT 1 instead of IGBT 3). They must therefore be considered.
The power layout is hard to modify, due to a lack of space and the mandatory position of the power modules on the cold plate. It has been shown that it is generating a current imbalance, but that this one is reversed when considering the Gate Driver PCB. Therefore, there is a possibility to compensate the imbalance caused by the power layout, using the coupling with the gate circuit. However, to guarantee this compensation, it is first necessary to consider both effects in a circuit model, combining an electrical representation of the power layout, the gate layout, their couplings, and the IGBT model. This is the aim of the following section.

3. Electrical Model of Paralleled IGBTs

Obtaining the circuit model from a 3D layout is quite straightforward using the PEEC approach for a simple conductor linking two points: the conductor is replaced by a resistance, an inductance, and mutual coupling with the rest of the circuit. However, in the kind of 3D geometry used in Figure 2, some conductors—as busbars, for instance—are not linking only two points but several points. Consequently, the electrical representation becomes tricky. In some cases, the geometrical layout is organized in such a way that the user can “guess” a circuit representation almost corresponding to the current path and the layout (Figure 9 top). However, in the case of plates, for instance, the current path is unknown, and the equivalent circuit topology not straightforward. A possibility would be to connect all points with R-L coupled circuits, as in Figure 9 bottom, but this is clearly not the representation using the minimal number of elements. The first subsection will propose a generic approach to get an electrical circuit from the 3D layout, based on a “terminal behavior” vision. This generic representation is used in several existing PEEC softwares, as Ansoft Q3D [19] or Altair Flux-PEEC [21], but has not been explained in such detail so far.
This circuit representation of both power and gate circuits will be combined with an IGBT conventional model in order to investigate the current division in the simplified case of two paralleled IGBTs. From this analysis, a rule on the impedance matrix will be derived to reach balanced currents in each IGBT.

3.1. Generic Circuit Representation

The results of the previous sections show the importance of considering all stray inductive elements involved in the dynamic behavior of the switching, including power, layout, and mutual couplings with the gate circuits. It is thus mandatory to take all of them into account. For simplification and illustration purposes, only two paralleled power modules will be considered in this section, even if, obviously, the generalization to three and more is quite straightforward with the proposed equations and matrix writing. We will focus on the IGBT 1 and the IGBT 3 of the bottom switch of phase W to illustrate the geometry.
The chosen electrical circuit representation of a conductor linking N_port > 2 will be composed of N_port-1 inductors and resistors. One point among the I/O ports will be considered as a reference and linked with the others with R-L circuits (and couplings). The reference point can be chosen arbitrarily. Figure 10 shows a simple illustration for the basic case of two plates associating three capacitors in parallel.
It is worth noting that this method is the one used in most PEEC-based softwares to generate Pspice-compatible netlists. One further remark is that in this representation, the full PEEC matrix exhibits both a real and an imaginary part. In other words, not only the inductors but also the resistors are coupled. Physically, this corresponds to a common resistive path within the same conductor and eddy currents induced in the rest of the geometry. These “resistive couplings” can be represented in circuit equations with current-controlled voltage sources (Figure 11), or by writing differential equations in a programming language (VHDL AMS or MAST for Saber software).

3.2. Application to Two Paralleled IGBTs

Let’s consider two paralleled IGBTs connected between phase W and the HT− negative potential (IGBT 1 and IGBT 3 of the full geometry). The PEEC method was applied to model the power busbar as well as the phase busbar W. The power module was also described using a PEEC model.
For the power busbar, the phase W connecting terminal was chosen as reference point for linking it to both collectors C1 and C2, and the HT− potential for linking to both power emitters, E1 and E2. For the power module, collectors and gate pins are simple connections linking the die to the external pin. Only Emitters have to be linked to the Die, Control Emitter, and Power Emitter pins. In this case, the reference point was the emitter die. Therefore, the circuit of Figure 12 left can be obtained, which can be reduced to the circuit of Figure 12 right by associating all impedances in series.
An important note is that this arbitrary circuit representation does not correspond to the geometrical layout: for instance, a part of the circuit can be common to the power and gate parts in the module (the contribution of a portion of DBC or bondings). This is not seen as a common inductance in the equivalent circuit of Figure 12, but it is taken into account in the value of mutual impedance between stray elements. This is illustrated further in Figure 13 for more clarity. In this figure, a part of the emitter path is common to the power and gate circuits. A geometrical representation would be as in Figure 13 left, but the generic representation is always the same whatever the layout geometry.

3.3. Circuit Equations

The method starts with considering the power circuit only, including a semiconductor model composed of a current source and the output capacitor, as illustrated in Figure 14.
I m o s = g m V G E V t h
To simplify writing, only inductive elements are displayed and considered in the equations, but all resistive elements and resistive couplings could be easily taken into account.
The expressions of the two power currents, IP1 and IP2, can be written as follows.
I P 1 = gm 1 V G E 1 V t h 1 + C C E 1 a d 2 I P 1 d t 2 b d 2 I P 2 d t 2
I P 2 = gm 2 V G E 2 V t h 2 + C C E 2 c d 2 I P 1 d t 2 d d 2 I P 2 d t 2
with
a = L C 1 + L E p 1 + 2 M C 1 E p 1
b = M C 1 C 2 + M C 1 E p 2 + M E p 1 C 2 + M E p 1 E p 2
c = M C 2 C 1 + M C 2 E p 1 + M E p 2 C 1 + M E p 2 E p 1
d = L C 2 + L E p 2 + 2 M C 2 E p 2
according to the notations of Figure 10 right, keeping the power elements only.
Considering identical devices, CCE1 = CCE2, gm1 = gm2, and Vth1 = Vth2. Furthermore, instead of solving the equation in the time domain, the frequency domain is preferred. Therefore, it leads to
I P 1 = g m V G E 1 V t h C C E b ω 2 I P 2 1 + C C E a ω 2
I P 2 = g m V G E 2 V t h C C E c ω 2 I P 1 1 + C C E d ω 2
Finally, combining (8) and (9), we obtain
I P 1 = g m C C E ω 2 ( a d b c d V G E 1 V t h b V G E 2 V t h
I P 2 = g m V G E 2 V t h C C E ω 2 d c d g m C C E ω 2 ( a d b c d V G E 1 V t h b V G E 2 V t h
Considering that the dynamic current share phenomena involves only high frequencies, CDS.a.ω2 >> 1 and CDS.d.ω2 >> 1. If we want to have IP1 = IP2, it leads to
V G E 1 = a + b c + d V G E 2 + V t h 1 a + b c + d
This means that to obtain an equal dynamic share despite dissymmetric layouts, it is possible to act on the gate voltages in the following way:
V G E 1 = α V G E 2 + V t h 1 α
with α = a + b c + d
It is worth noting that if the power layout is symmetrical, a = d and b = c; therefore, α = 1 and thus VGS1 = VGS2: it is not necessary to compensate the power currents from the gate side.
Once the power layout has been defined, the current dissymmetry is defined (i.e., a, b, c, d, and thus α), and it is possible to modify the gate circuit, in order to obtain different values of VGE1 and VGE2.
These gate voltages can be expressed as a function of gate and power currents and the various stray impedances of the circuit from Figure 12 right:
V G E 1 V G E 2 = V d r i v e r V d r i v e r j ω Z G | Z P G I G 1 I G 2 I P 1 I P 2
ZG being the stray impedance matrix on the gate distribution PCB, Vdriver being the driver output voltage (common for all devices), and ZPG being the coupling matrix between gate distribution PCB and power circuit. This coupling matrix contains all mutual impedances the between power and the gate circuit, as already explained to introduce Figure 12 right.
Neglecting gate currents in comparison with power currents leads to focus on ZPG only. Therefore, the following expressions can be obtained for VGE1 and VGE2:
V G E 1 = V d r i v e r j ω Z P G   l i n e 1 I P 1
V G E 2 = V d r i v e r j ω Z P G   l i n e 2 I P 2
If we want to obtain IP1 = IP2, combining (11) and (12) with (9) leads to:
j ω I P Z P G   l i n e 1 α Z P G   l i n e 2 V t h V d r i v e r 1 α   = 0
the summation term is just here to show that more than two paralleled devices can be considered.
Therefore, to reach equal currents, the ZPG matrix has to be designed in such a way that (13) is fulfilled, in order to compensate the power dissymmetry, expressed by a non-unity value of α. In other words, all couplings between gate PCB and power layout should be designed in order to compensate the non-ideal behavior of the power part. This is a major improvement in comparison to previous works where only power dissymmetry [7,8,9,10] or power-drive interaction [11,12,13,14,15] were considered, but not the possibility of compensating power dissymmetry through the gate circuit layout, using the device’s model, as in [16,17]. This target will be reached through optimization in the following section.

4. Optimization of the Gate Distribution PCB

Finding the layout of the Gate Distribution PCB which, coupling with the power circuit, fits (13) is not possible with a human brain, even with years of trial-and-error work! Therefore, an optimization algorithm has been set up. The principle is illustrated in Figure 15: a home-made optimization environment (RSM, developed by the Alstom group, Semeac France), controls the Ansoft Q3D PEEC software, in which a parametrized geometry of the Gate Distribution PCB and the power layout is described. By varying the geometrical parameters of the PCB, the objective function defined by (13) is minimized. The optimization used in the RSM environment is based on genetic algorithms NSGA2 [22] and will not be described in this paper. The parametrized geometry of the Gate Distribution PCB is given in Figure 16. For simplicity, a multilayer description has been chosen, in order to draw the two gate tracks and the two emitter tracks independently. Track crossing in each layer is also allowed, to simplify the description. After optimization, the number of layers can be easily reduced, and track crossing can be avoided using vias. The geometrical description of the gate connections involves four parameters per gate, whereas the source connections are simpler and defined using two parameters per IGBT. Parameter description and an example of corresponding geometrical constraints are given in Figure 16 and Table 1, respectively. Choosing these 12 parameters allows the layout building some “loops”, which can provide positive or negative couplings with the power part, and thus allow sufficient degrees of freedom for the algorithm to reach the goal defined by the objective function. All connection points to the power modules are obviously fixed, since the power layout and the module location remain unchanged.
The optimization function is frequency-dependent: first of all, because equations were written in the frequency domain, and also because impedance evaluation in the PEEC software is carried out at a specific frequency. Since the goal is to balance currents during the switching phase, an equivalent frequency corresponding to the switching time, 0.35/tswitching, has been used [23]. It has been evaluated to 50 MHz. The impact of this choice will be discussed later.
The main interest of the proposed method is that it does not imply any time simulation in the optimization loop, which is really time-consuming and often generates convergence issues: some configurations generated by the optimizer may indeed not be able to converge in the time domain, and the optimization process may thus fail.
After several hours of optimization, which is quite a reasonable amount of time, since the process is automatic, a result has been obtained. It is displayed in Figure 17. The geometry displayed on the top of the figure is compatible with the imposed constraints (fixed connection points and maximum size of PCB).
On the bottom of Figure 17, the evolution of the objective function is plotted as a function of the number of iterations of the algorithm. This representation helps in evaluating the convergence of the NSGA2 algorithm. It can be seen, for instance, that after iteration 400, the minimal value of the objective function has been reached and is not really modified between 400 and 1050, where it has been stopped. This minimal value converges up to 225,000, which is lower than the initial value, around 260,000. This objective function is the left part of (13). It has not reached zero due to all the constraints, but it has been reduced compared to the initial situation.
Since some assumptions have been made to build the objective function, and since it has been evaluated at a specific frequency (50 MHz), there is a strong need to validate the obtained result in the time domain. For this purpose, the PEEC model of the 3D layout (power and gate circuit in the same file for including couplings) has been combined with precise models of IGBTs in the Simplorer time domain simulator (Figure 18).
The initial geometry of the gate distribution PCB and the power layout has been simulated as well (light/dark blue curves in Figure 19). Results are compared to the optimized gate distribution PCB (magenta/red in Figure 19). The improvement is really obvious, which validates both the methodology and the optimized geometry.

5. Conclusions

This paper has presented a 3D layout of a power inverter used in a railway application. Due to the high-speed switching of candidate SiC MOSFETs, specific care has to be taken with regard to stray inductances and mutual couplings generated by the interconnections. In particular, it has been shown that the present 3D layout of the existing Si inverter was exhibiting uneven current division among the paralleled IGBTs. It was therefore mandatory to improve this layout before going to the next step with SiC devices. Using the PEEC approach, it has been demonstrated that the power layout was at the origin of a dynamic current imbalance, and that the power-gate interaction caused another imbalance. Modifying the power layout was not possible due to several other constraints. However, it was possible to design a Gate Distribution PCB in such a way that the couplings between this PCB and the power layout compensate the non-ideal behavior of the power part. For this purpose, a detailed analysis of the circuit equations has been carried out, using a generic representation independent from the geometry. This has led to an objective function to be minimized, which has been implemented in an optimization environment controlling a PEEC based software. The optimized layout has been validated with a time domain simulation.

Author Contributions

Conceptualization, J.-L.S., J.-M.G. and P.-O.J.; methodology, J.-M.G. and E.B.; software, J.-M.G. and E.B.; validation, P.-O.J. and A.N.-M.; formal analysis, J.-M.G. and A.N.-M.; investigation, A.N.-M.; resources, A.N.-M. and E.B.; data curation, A.N.-M. and E.B.; writing—original draft preparation, A.N.-M. and J.-L.S.; writing—review and editing, J.-L.S.; visualization, A.N.-M.; supervision, J.-L.S., J.-M.G. and P.-O.J.; project administration, G.D. and J.-L.S.; funding acquisition, G.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ALSTOM, within the European Project Shift2rail.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shift2Rail Project. European Funding for Research & Innovation. Available online: https://shift2rail.org/ (accessed on 7 March 2021).
  2. Rabkowski, J.; Peftitsis, D.; Nee, H. Silicon Carbide Power Transistors: A New Era in Power Electronics Is Initiated. IEEE Ind. Electron. Mag. 2012, 6, 17–26. [Google Scholar] [CrossRef]
  3. Seal, S.; Mantooth, H.A. High Performance Silicon Carbide Power Packaging—Past Trends, Present Practices, and Future Directions. Energies 2017, 10, 341. [Google Scholar] [CrossRef] [Green Version]
  4. Schanen, J.; Jeannin, P. Integration solutions for clean and safe switching of high speed devices. In Proceedings of the CIPS 2018 10th International Conference on Integrated Power Electronics Systems, Stuttgart, Germany, 20–22 March 2018; pp. 1–11. [Google Scholar]
  5. Hatanaka, A.; Kageyama, H.; Masuda, T. A 160-kW high-efficiency photovoltaic inverter with paralleled SiC-MOSFET modules for large-scale solar power. In Proceedings of the 2015 IEEE International Telecommunications Energy Conference (INTELEC), Osaka, Japan, 18–22 October 2015; pp. 1–5. [Google Scholar] [CrossRef]
  6. Schrader, R.; Cordt, D.; Münster, P.; Eckel, H. Influences of Gate-Circuit and Parasitic Inductances on Turn-Off Current Imbalances of Paralleled IGBTs Due to Differences in Their Switching Behaviour. In Proceedings of the 2018 20th European Conference on Power Electronics and Applications (EPE’18 ECCE Europe), Riga, Latvia, 17–21 September 2018; pp. P.1–P.10. [Google Scholar]
  7. Jeannin, P.; Schanen, J.; Clavel, E. Original cabling conditions to insure balanced current during switching transitions between paralleled semiconductors. IEEE Trans. Ind. Appl. 2002, 38, 181–188. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Zang, Z.; Hu, B.; Alsaif, F.; Alkhalid, K.; Li, X.; Wang, J. Current Sharing Analysis of SiC Power Modules in Parallel Operation. In Proceedings of the2020 IEEE Energy Conversion Congress and Exposition (ECCE), Detroit, MI, USA, 15–20 October 2020; pp. 711–716. [Google Scholar] [CrossRef]
  9. Wissen, M.; Domes, D.; Brekel, W.; Yilmaz, K. How asymmetric busbar design causes symmetric switching behavior of paralleled IGBT modules. In Proceedings of the CIPS 2018, 10th International Conference on Integrated Power Electronics Systems, Stuttgart, Germany, 20–22 March 2018; pp. 1–3. [Google Scholar]
  10. Mainali, K.; Wang, R.; Sabate, J.; Klopman, S. Current Sharing and Overvoltage Issues of Paralleled SiC MOSFET Modules. In Proceedings of the 2019 IEEE Energy Conversion Congress and Exposition (ECCE), Baltimore, MD, USA, 29 September–3 October 2019; pp. 2413–2418. [Google Scholar] [CrossRef]
  11. Li, H.; Zhou, W.; Wang, X.; Nielsen-Munk, S.; Li, D.; Wang, Y.; Dai, X. Influence of Paralleling Dies and Paralleling Half-Bridges on Transient Current Distribution in Multichip Power Modules. IEEE Trans. Power Electron. 2018, 33, 6483–6487. [Google Scholar] [CrossRef]
  12. Bęczkowski, S.; Jørgensen, A.B.; Li, H.; Uhrenfeldt, C.; Dai, X.; Munk-Nielsen, S. Switching current imbalance mitigation in power modules with parallel connected SiC MOSFETs. In Proceedings of the 2017 19th European Conference on Power Electronics and Applications (EPE’17 ECCE Europe, Warsaw, Poland, 11–14 September 2017; pp. P.1–P.8. [Google Scholar] [CrossRef] [Green Version]
  13. Li, H.; Munk-Nielsen, S.; Bȩczkowski, S.; Wang, X. A Novel DBC Layout for Current Imbalance Mitigation in SiC MOSFET Multichip Power Modules. IEEE Trans. Power Electron. 2016, 31, 8042–8045. [Google Scholar] [CrossRef]
  14. Li, H.; Munk-Nielsen, S.; Bęczkowski, S.; Wang, X.; Eni, E. Effects of auxiliary source connections in multichip power module. In Proceedings of the 2016 IEEE Applied Power Electronics Conference and Exposition (APEC), Long Beach, CA, USA, 21–24 March 2016; pp. 3101–3106. [Google Scholar] [CrossRef]
  15. Martin, C.; Guichon, J.; Schanen, J.; Pasterczyk, R. Gate Circuit Layout Optimization of Power Module Regarding Transient Current Imbalance. IEEE Trans. Power Electron. 2006, 21, 1176–1184. [Google Scholar] [CrossRef]
  16. Tang, Y.; Ma, H. Dynamic Electrothermal Model of Paralleled IGBT Modules with Unbalanced Stray Parameters. IEEE Trans. Power Electron. 2017, 32, 1385–1399. [Google Scholar] [CrossRef]
  17. Zeng, Z.; Zhang, X.; Li, X. Layout-Dominated Dynamic Current Imbalance in Multichip Power Module: Mechanism Modeling and Comparative Evaluation. IEEE Trans. Power Electron. 2019, 34, 11199–11214. [Google Scholar] [CrossRef]
  18. Ruehli, A.E. Equivalent Circuit Models for Three-Dimensional Multiconductor Systems. IEEE Trans. Microw. Theory Tech. 1974, 22, 216–221. [Google Scholar] [CrossRef]
  19. Ansys Q3D. Available online: https://www.ansys.com/products/electronics/ansys-q3d-extractor (accessed on 7 March 2021).
  20. Nunnally, W.C.; Kristiansen, M.; Hagler, M.O. Differential Measurement of Fast Energy Discharge Capacitor, Inductance, and Resistance. IEEE Trans. Instrum. Meas. 1975, 24, 112–114. [Google Scholar] [CrossRef]
  21. Altair Flux PEEC. Available online: https://www.altair.com/resource/conductor-impedance-and-near-field-simulation-using-altair-flux (accessed on 7 March 2021).
  22. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  23. Brooks, D. Signal Integrity Issues and Printed Circuit Board Design; Prentice Hall: Upper Saddle River, NJ, USA, 2004; ISBN 139780133359473. [Google Scholar]
Figure 1. Railway inverter composition. The global function is a three-phase inverter with one additional non-reversible leg for braking. Each leg uses three paralleled IGBT modules. The full converter is built from two identical subparts composed of two legs and decoupling capacitors (bottom of the figure). A common DC capacitor tank is used to provide the energy storage at DC-bus level.
Figure 1. Railway inverter composition. The global function is a three-phase inverter with one additional non-reversible leg for braking. Each leg uses three paralleled IGBT modules. The full converter is built from two identical subparts composed of two legs and decoupling capacitors (bottom of the figure). A common DC capacitor tank is used to provide the energy storage at DC-bus level.
Energies 14 01891 g001
Figure 2. (Top: left), DC capacitor bank with duplicated DC bus connections, one for each subpart. (Center), overview of the global 3D layout. (Right), IGBT power Module. (Bottom): 3D layout of the power part.
Figure 2. (Top: left), DC capacitor bank with duplicated DC bus connections, one for each subpart. (Center), overview of the global 3D layout. (Right), IGBT power Module. (Bottom): 3D layout of the power part.
Energies 14 01891 g002
Figure 3. Zoom of the slit in the phase busbar to provide almost equal length for the current path between the phase output and the power modules.
Figure 3. Zoom of the slit in the phase busbar to provide almost equal length for the current path between the phase output and the power modules.
Energies 14 01891 g003
Figure 4. Gate Distribution PCB, collecting the driver signal and linking it to each of the three paralleled IGBTs.
Figure 4. Gate Distribution PCB, collecting the driver signal and linking it to each of the three paralleled IGBTs.
Energies 14 01891 g004
Figure 5. Time domain simulation scheme. Two Partial Element Equivalent Circuit (PEEC) models used from separate files: power and phase busbar plus Gate Distribution PCB and DC capacitor tank busbar. Equivalent esr–esl models of the decoupling capacitors and the three blocks of the DC Capacitor tank. “Home made” IGBT and diode models are used. Idealized drive signal is a simple voltage source.
Figure 5. Time domain simulation scheme. Two Partial Element Equivalent Circuit (PEEC) models used from separate files: power and phase busbar plus Gate Distribution PCB and DC capacitor tank busbar. Equivalent esr–esl models of the decoupling capacitors and the three blocks of the DC Capacitor tank. “Home made” IGBT and diode models are used. Idealized drive signal is a simple voltage source.
Energies 14 01891 g005
Figure 6. Simulated results of current division among the three bottom switches of phase W. Full 3D model.
Figure 6. Simulated results of current division among the three bottom switches of phase W. Full 3D model.
Energies 14 01891 g006
Figure 7. Simulated results of current division among the three bottom switches of phase W with power layout only, identical idealized gate signals.
Figure 7. Simulated results of current division among the three bottom switches of phase W with power layout only, identical idealized gate signals.
Energies 14 01891 g007
Figure 8. Simulated results of current divisions among the three bottom switches of phase W with power layout and Gate Distribution PCB layout—no coupling considered.
Figure 8. Simulated results of current divisions among the three bottom switches of phase W with power layout and Gate Distribution PCB layout—no coupling considered.
Energies 14 01891 g008
Figure 9. Electrical representation of conductors linking more than two points. (Top), the current path is quite obvious, and an electrical representation can be “guessed”. (Bottom), the case of a plate with no straightforward representation, except for the “heavy” full connection of each point with R-L coupled elements (mutual coupling not displayed for clarity).
Figure 9. Electrical representation of conductors linking more than two points. (Top), the current path is quite obvious, and an electrical representation can be “guessed”. (Bottom), the case of a plate with no straightforward representation, except for the “heavy” full connection of each point with R-L coupled elements (mutual coupling not displayed for clarity).
Energies 14 01891 g009
Figure 10. Generic representation, based in an input/output port choosing an arbitrary reference port. This representation is NOT linked with any geometrical current path.
Figure 10. Generic representation, based in an input/output port choosing an arbitrary reference port. This representation is NOT linked with any geometrical current path.
Energies 14 01891 g010
Figure 11. A possible way to reproduce “resistive couplings” in Pspice-compatible representations.
Figure 11. A possible way to reproduce “resistive couplings” in Pspice-compatible representations.
Energies 14 01891 g011
Figure 12. (Left): individual models of each part of the system. (Right): reduction of the model by associating series elements.
Figure 12. (Left): individual models of each part of the system. (Right): reduction of the model by associating series elements.
Energies 14 01891 g012
Figure 13. Example of a PEEC representation in the case of classical common emitter coupling. The (left) representation is based on the geometry and the proposed generic representation; (right) is based on terminals only and is independent from the geometry. Both account for exactly the same phenomena.
Figure 13. Example of a PEEC representation in the case of classical common emitter coupling. The (left) representation is based on the geometry and the proposed generic representation; (right) is based on terminals only and is independent from the geometry. Both account for exactly the same phenomena.
Energies 14 01891 g013
Figure 14. Representation of the power part only, including a simple IGBT model.
Figure 14. Representation of the power part only, including a simple IGBT model.
Energies 14 01891 g014
Figure 15. Optimization methodology.
Figure 15. Optimization methodology.
Energies 14 01891 g015
Figure 16. Gate Distribution PCB: parametrized geometry.
Figure 16. Gate Distribution PCB: parametrized geometry.
Energies 14 01891 g016
Figure 17. (Top): Optimized Gate Distribution PCB corresponding to the optimal point. (Bottom): evolution of the objective function vs iterations. The green point corresponds to the minimum of the objective function.
Figure 17. (Top): Optimized Gate Distribution PCB corresponding to the optimal point. (Bottom): evolution of the objective function vs iterations. The green point corresponds to the minimum of the objective function.
Energies 14 01891 g017
Figure 18. Post-processing of the optimization results: after geometric optimization, the full geometry is exported to Simplorer for a time domain simulation.
Figure 18. Post-processing of the optimization results: after geometric optimization, the full geometry is exported to Simplorer for a time domain simulation.
Energies 14 01891 g018
Figure 19. Results of the time domain simulation: in blue, the initial geometry of the gate distribution PCB; in red/pink, the optimized PCB.
Figure 19. Results of the time domain simulation: in blue, the initial geometry of the gate distribution PCB; in red/pink, the optimized PCB.
Energies 14 01891 g019
Table 1. Example of geometrical constraints associated to parametrized geometry: example for A, B, C, D.
Table 1. Example of geometrical constraints associated to parametrized geometry: example for A, B, C, D.
Parameter NameMin ValueMax Value
A0Amax = PCB max size X
B0Bmax = PCB max size Y
C0Cmax = PCB max size X
D0Dmax = PCB max size X
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nakahata-Medrado, A.; Schanen, J.-L.; Guichon, J.-M.; Jeannin, P.-O.; Batista, E.; Desportes, G. Optimization of a Gate Distribution Layout to Compensate the Current Imbalance Generated by the 3D Geometry of a Railway Inverter. Energies 2021, 14, 1891. https://doi.org/10.3390/en14071891

AMA Style

Nakahata-Medrado A, Schanen J-L, Guichon J-M, Jeannin P-O, Batista E, Desportes G. Optimization of a Gate Distribution Layout to Compensate the Current Imbalance Generated by the 3D Geometry of a Railway Inverter. Energies. 2021; 14(7):1891. https://doi.org/10.3390/en14071891

Chicago/Turabian Style

Nakahata-Medrado, Andressa, Jean-Luc Schanen, Jean-Michel Guichon, Pierre-Olivier Jeannin, Emmanuel Batista, and Guillaume Desportes. 2021. "Optimization of a Gate Distribution Layout to Compensate the Current Imbalance Generated by the 3D Geometry of a Railway Inverter" Energies 14, no. 7: 1891. https://doi.org/10.3390/en14071891

APA Style

Nakahata-Medrado, A., Schanen, J. -L., Guichon, J. -M., Jeannin, P. -O., Batista, E., & Desportes, G. (2021). Optimization of a Gate Distribution Layout to Compensate the Current Imbalance Generated by the 3D Geometry of a Railway Inverter. Energies, 14(7), 1891. https://doi.org/10.3390/en14071891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop