Next Article in Journal
The Effect of Supercritical CO2 on Shaly Caprocks
Next Article in Special Issue
Nonlinear Characteristics Compensation of Inverter for Low-Voltage Delta-Connected Induction Motor
Previous Article in Journal
Dual Deep Learning Networks Based Load Forecasting with Partial Real-Time Information and Its Application to System Marginal Price Prediction
Previous Article in Special Issue
Optimization of a Cascaded H-Bridge Inverter for Electric Vehicle Applications Including Cost Consideration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Optimization of Three-Phase Dual-Active-Bridge Converters for Electric Vehicle Charging Stations

1
Aichi Institute of Technology, Toyota 470-0356, Japan
2
Hanoi University of Science and Technology, Hanoi 100803, Vietnam
3
Shibaura Institute of Technology, Tokyo 135-8548, Japan
*
Author to whom correspondence should be addressed.
Energies 2020, 13(1), 150; https://doi.org/10.3390/en13010150
Submission received: 27 November 2019 / Revised: 20 December 2019 / Accepted: 25 December 2019 / Published: 27 December 2019
(This article belongs to the Special Issue Electric Vehicle Power Conversion Technologies)

Abstract

:
In this paper, design and optimization method of a three-phase dual-active-bridge DC/DC converter is discussed. Three single phase transformers connected in star-star configuration were designed with large leakage inductance aiming to eliminate the need for external inductors. Switching frequency, peak flux density, number of turns, number of layers, etc., were optimized using non-linear programming technique for minimizing the overall converter loss. Experimental results on a 10 kW prototype show that the optimized converter can operate efficiently an efficacy of up to 98.65% and a low-temperature rise of less than 70 degrees Celsius on both transformers and semiconductor devices.

1. Introduction

Electric vehicles (EVs) play an importance role in reducing CO 2 emission, as they are powered by chemical batteries or super capacitors, rather than by fuel combustion. Recently, some new ideas have been proposed to make use of EVs more effectively. Among those, the idea of using EVs as power system stabilizers via vehicle-to-grid (V2G) or vehicle-to-home (V2H) technologies is very interesting and promising thanks to their significant storage capabilities [1,2]. Accordingly, the EV would be connected to the utility via a so-called EV supply equipment (EVSE) and a charging outlet (such as Chademo, CCS/Combo, etc.) when it is not in used in a quite long time. Depending on the actual needs of the grid, the EV could be charged or discharged to help stabilize the grid. In order for the energy stabilizer mode to function, the EVSE must be a bidirectional converter that allows the power to transfer bidirectionally. An example of such the system is depicted in Figure 1.
The EVSE in Figure 1 consists of two parts: a unidirectional rectifier and a bidirectional DC/DC charger. The output of the rectifier is tied to a 380 V DC bus of a DC micro-grid, from which the bidirectional DC/DC converter charges the electric vehicle via a CHAdeMO outlet. When the vehicle is left unused, the outlet remains connected for charging or for reversing the power back to the DC bus. Thanks to that, the DC micro-grid can be stabilized. This paper focuses on designing and optimization of the bidirectional DC/DC converter of such an EVSE (from now on, the DC charger).
From another aspect, according to [3], high frequency current ripple causes negative effect on the battery lifetime. Therefore, it is necessary to design a DC charger with low output current ripple. It should also have high power density, high efficiency, be galvanic isolated, and be bidirectional. Among many DC/DC converter topologies, dual-active-bridge (DAB) has been found to be the most appropriate choice for this application [4,5,6].
There are two popular types of DAB converter: single-phase and three-phase, both having inherent soft-switching, bidirectional, and isolation characteristics. Unlike the single-phase counter part, the three-phase structure can triple the effective frequency and reduce the amplitude of the ripple. Hence, input and output filters can be downsized; the transformer RMS current is also lower and power density is higher than that of the single-phase version. Up till now, numerous studies on three-phase DAB converter have been published [6,7,8,9,10] regarding various topics, including modulation and control to system design and optimization, etc. While many discussions on designing and optimizing the single-phase DAB can be found in the literature, only a few are available for three-phase topology.
Of course, the design and optimization procedures that are applied for single-phase DAB converter can be applicable for three-phase DAB. For example, the typical approach is starting from loss modeling, formulating the optimization problem, and then solving the problem by using some techniques, such as that reported in [11]. However, in [11], an external inductor was used, and it is preferable to integrate the inductor into the transformer for better space saving and more efficient performance. Moreover, the switching frequency was fixed in that research; only the number of turns and the inductance were optimized. The conduction loss model used in that research did not address the dependency of AC resistance of winding on the proximity effect, which varies correspondingly to the number of turns and layers [12].
In [13], the inductor-integrated transformer in a single-phase DAB converter was designed and optimized by using the particle swarm optimization method. The skin and proximity effects were regarded in loss modeling. However, conduction loss and switching loss of semiconductor devices were not considered. Therefore, the overall converter system was not optimized; only transformers were. The research reported in [14] modeled both electronics loss and transformer loss. The transformer used shell-type winding for better integration of leakage inductance. However, only the control variables (phase shift, switching frequency) were optimized. The transformer structure was fixed and not optimized.
Design and optimization of three-phase DAB converter for DC charging applications are resolved in this paper. Both transformer structure and operation parameters are optimized. Inductors are integrated into the transformers for better space and cost saving. The value of leakage inductance, number of turns, number of layers, etc., are also optimized. For better accuracy, skin and proximity effects are also addressed in modeling the copper loss of transformers. A nonlinear programming technique was employed to find the optimal parameter set. Due to certain facts, including that the window area of the transformer is limited, the semiconductor devices have a certain switching capability, the operating temperature needs to be handled as well, etc., some inequality constraints were considered for the optimization problems. Afterward, the optimized transformer was verified by a finite element method magnetic (FEMM) analysis and by an experimental study.
This paper is actually an improved version of that reported in [15]. The converter was modeled with consideration that number of turns and number of layers are integers. That makes the optimizing function discontinuous and non-differentiable. Therefore, a nonlinear programming method based on gradient, such as the Newton–Raphson, becomes inapplicable. Besides, in order to find the global optima, the Global Optimization Toolbox of MATLAB & Simulink from Mathworks incorporation was employed. Furthermore, operation in the whole power range was examined in this improved version. The paper is structured as follows: analysis in the steady state is presented in Section 2; transformer design method is described in Section 3; the loss model can be found in Section 4; the optimization algorithm is narrated in Section 5; the optimization, simulation, and experimental results are provided in Section 6; and finally, some conclusions are given in Section 7.

2. Steady State Analysis

Figure 2 describes the circuit diagram of a three-phase dual-active-bridge (from now on, DAB3) converter. Three single phase high-frequency transformers are employed for easier fabrication and implementation of the system. It assumes that all transformers have the same leakage inductance of L k and the same winding ratio of n : 1 . According to [8], among some transformer configurations, the Y–Y topology requires smaller leakage inductance to send a specific amount of power. That also means a smaller number of turns of transformer winding is required, resulting in smaller winding resistance. Therefore, in this study, the Y–Y connection is used aiming for performance optimization. Notes here that, delta-type transformer structures can be transformed into star connection by using the following transformation:
L Y a = L D b L D c L D a + L D b + L D c .
Two three-phase inverters are located at two sides of the transformers. The conventional phase shift modulation and 180-degree-conduction-mode are used to control the converter. The switching frequency is f s which needs to be optimized. Voltages at the two DC terminals are V 1 and V 2 , respectively. The voltages across the transformer winding have the four-level form and are ψ degrees shifted from each other as shown in Figure 3, where ψ is the phase shift control angle.
Assuming that the phase shift ψ is smaller than π 3 for better reactive power reduction, according to [8], the transition current I 0 to I 5 at the steady state can be calculated by (1):
I 0 = I M 2 ( 1 M ) + 3 M ψ π I 1 = I M 2 ( 1 M ) + 3 ψ π I 2 = I M ( 1 M ) + 3 M ψ π I 3 = I M ( 1 M ) + 6 ψ π I 4 = I M ( 1 M ) + 6 M ψ π I 5 = I M ( 1 M ) + 3 ψ π
where I M = V 1 18 f s L k ; M is the voltage conversion ratio, M = n V 2 V 1 ; and ψ is the phase shift, ψ < π 3 .
The root-mean-squared (rms) transformer current, rms switch currents, peak flux density (obtained by calculating the volt-second product of each conducting interval and assuming a balanced volt-second distribution of transformers), and output power are calculated by equations from (2) to (6) as follows:
I p r i , r m s = I M 3 3 5 ( 1 M ) 2 + 27 2 ψ π ψ π 2 M
I Q x , r m s = I p r i , r m s 2
I S x , r m s = I Q x , r m s n
B p k = V 1 18 f s N 1 A c 1 + M 1 3 ψ 2 π
P o u t = 3 M V 1 I M 2 ψ 4 3 ψ π ,
where N 1 is the number of turns of the primary winding; A c is the cross section area of the center limb of the magnetic core, assuming that E core is used.

3. Transformer Design

At the first step, some preliminary parameters must be determined:
  • Maximum output power;
  • Nominal input and output voltages;
  • Current density (J = 6 − 10 A/mm 2 );
  • Window utilization factor ( K u = 0.2 − 0.4 for transformer);
  • Initial switching frequency and flux density.
After that, from (2), (5) and (6), the core size and winding wire size can be selected according to Ap or Kg methods [16]. Here, in order to reduce the skin and proximity effects, Litz wire is preferable. Three transformers with the traditional core-typed winding structure, as depicted in Figure 4, were designed. There are only two windings in each transformer, primary and secondary; no interleaving. The centralized wounding technique was used. The insulation between two windings was made thick for leakage inductance integration purpose. The insulation thickness is, thus, a variable for tuning the leakage inductance.
According to [16], the leakage inductance can be calculated by:
L k = μ 0 ( M L T ) N 1 2 a c + b 1 + b 2 3 ,
where μ 0 is the permeability of vacuum, μ 0 = 4 π × 10 7 ; M L T is the mean-length-turn of the core; a is the winding height; b 1 and b 2 are the windings’ thicknesses; c is the insulation thickness; and N 1 is the number of turns of primary winding, from (5). N 1 can be derived as (8):
N 1 = ceil V 1 18 f s B p k A c 1 + M 1 3 ψ 2 π
where ceil ( ) is the ceiling (round-up) function.
Now, let m 1 and m 2 be the number of layers of primary and secondary windings; the winding dimension will be:
a = O D 1 × N 1 m 1
b 1 + b 2 = O D 1 m 1 + O D 2 m 2 ,
where O D 1 and O D 2 are the outer dimensions of winding wires. Let H be the window height of the bobbin. The winding height a must be smaller than H. From (9), we get:
N 1 H O D 1 m 1 .
Due to the round shape of the wire and the necessary space for creepage and clearance distance, the whole core height cannot be used for winding. Therefore, by denoting k c 1 as the utilization factor of the core height, k c 1 0.9 , m 1 and m 2 can be determined by:
m 1 = ceil N 1 O D 1 k c 1 H
m 2 = ceil N 2 O D 2 k c 1 H .
Since b 1 and b 2 depend on the chosen current density J, and have less effect on leakage inductance value, as formulated in (7), the tuning parameter for L k is the insulation thickness c. Note that the summation ( b 1 + b 2 + c ) must be smaller than the total window width W of the bobbin. The leakage inductance is highest at c = c m a x :
c c m a x = W ( b 1 + b 2 ) .
From (7)–(14), the maximum achievable leakage inductance for a given core is determined by:
L k L k , m a x = μ 0 ( M L T ) H O D 1 m 1 2 W O D 1 2 3 m 1 + O D 2 O D 1 m 2 .
The leakage inductance is minimum when there is no insulation between the two windings ( c = 0 ):
L k L k , m i n = μ 0 ( M L T ) N 1 m 1 3 m 1 + O D 2 O D 1 m 2 .
From (6), for a given output power of P m , the maximum phase shift is calculated by:
ψ m = 2 π 3 1 1 f s f s , m a x ,
where f s , m a x = M V 1 2 9 L k P m . Obviously, in order for (17) to be valid, the switching frequency must be less than the maximum value of f s , m a x :
f s f s , m a x = M V 1 2 9 L k P m .
From (8), (17) and (18), we have:
f s V 1 18 N 1 A c B p k 1 + M 1 f s f s , m a x .
Solving (19) for f s , the lower boundary of switching frequency is:
f s f s , m i n = λ f s , m a x 1 λ 2 M 2 + M 1 λ + λ 2 4 M 2 ,
where λ = P 2 M V 1 × L k N 1 A c B p k .

4. Power Loss Modeling

4.1. Power Electronics Loss

Assuming that zero-voltage transition is achieved, turn-on loss can thus be ignored. The overall power dissipation on the Silicon-Carbide Metal–Oxide–Semiconductor Field-Effect Transistor (SiC-MOSFETs) can be estimated by:
Δ P p e = Δ P c o n d + Δ P s w ,
where Δ P c o n d and Δ P s w are conduction loss and turn-off loss respectively:
Δ P c o n d = 6 R d s , O N ( I Q x , r m s 2 + n 2 I S x , r m s 2 ) Δ P s w = 3 V 1 ( I 0 + M I 1 ) t f f s
with R d s , O N and t f being the ON resistance and the falling time of the MOSFET; both are given in the device data-sheet; I 0 and I 1 are the transition currents given in (1); I Q x , r m s and I S x , r m s are the rms currents flowing through primary and secondary switches, which are given in (3) and (4), respectively.

4.2. Transformer Loss

Transformer loss includes copper loss dissipating on windings and ferrite loss on the magnetic core as follows:
Δ P t r = Δ P C u + Δ P F e ,
where Δ P C u and Δ P F e are the winding loss (or copper loss) and core loss (ferrite loss). Using the Steimetz equation, the core loss can be estimated by:
Δ P F e = 3 k c f s α B p k β g ,
where B p k is the peak flux density, and g is the weight of two core halves. k c , a and b are the Steimetz coefficients: k c = 0.00004855 , α = 1.62 , and β = 2.63 [16].
The copper loss is determined by (24):
Δ P C u = 3 R a c I p r i , r m s 2 ,
where R a c is the primary-referred AC resistance. Note that the AC resistance depends strongly on the switching frequency and the winding geometry by the skin and proximity effects. According to [12], the AC resistance can be estimated for Litz wire as follows:
R a c = R d c A s t r sinh 2 A s t r + sin 2 A s t r cosh 2 A s t r cos 2 A s t r + 2 3 m 2 n s t r 1 sinh A s t r sin A s t r cosh A s t r + cos A s t r ,
where A s t r = 1 δ π 4 0.75 d b 3 d w ; δ is the skin depth, δ = 6.62 × 10 2 f s ; d b and d w are the bare and outer diameters of a single strand of the Litz wire; m is the number of layers; n s t r is the number of strands; R d c is the DC resistance of the winding:
R d c = ( M L T ) N r s t r n s t r ,
where N is the number of turns (i.e., N 1 for primary winding and N 2 for secondary winding); r s t r is the resistance density of a strand; and M L T is the mean-length-turn of the magnetic core.

4.3. Capacitor Loss

Since there is ripple in the input and output currents of the converter, the ripple will be filtered out by the input and output capacitors and dissipate some power on the series resistance of the capacitors. The capacitor loss is determined by:
Δ P c a p = r c i n I c i n , r m s 2 + r c o u t I c o u t , r m s 2 ,
where r c i n and r c o u t are the series resistances of capacitors; I c i n , r m s and I c o u t , r m s are the input and output rms capacitor currents which can be easily derived from (1).
Finally, the total power dissipation of the converter is derived as:
Δ P t o t = Δ P p e + Δ P t r + Δ P c a p .

5. Optimization

In order to maximize system efficiency, the overall power dissipation should be minimized. The analysis above shows that the four most basic parameters that affect all others are the switching frequency f s , the peak flux density B p k , the voltage conversion ratio M, and the leakage inductance L k . When f s , B p k , M, and L k are known, all other parameters can be easily derived.
Here, some additional constraints are considered. First, after designing the transformer, according to [16], one can estimate the temperature rise of the core by (28).
Δ T t r = 450 Δ P t r 3 A t 0.826 ,
where A t is the total surface of the transformer; A t can be calculated without difficulty when the core dimensions are known.
The temperature rise of power electronics devices can also be estimated if the thermal resistance of the heatsink R h s is given. Assuming that the heat generated from semiconductor devices is evenly distributed on the heatsink, then the temperature rise is:
Δ T p e = Δ P p e R h s .
The maximum temperature rise should not be greater than a certain allowable value of Δ T m a x .
max Δ T t r , Δ T p e Δ T m a x .
It is also worth noticing that all analyses above are conducted under the assumption that soft-switching is achieved. According to [8] and considering the required deadtime T d in modulation, the lowest soft-switching boundary when ψ < π / 3 is:
ψ ψ m i n = 2 π 3 1 M + 2 π f s T d .
At that operational point, the transferred power is:
P m i n = M V 1 2 12 f s L k ψ m i n 4 3 π ψ m i n .
The rated power P should be greater than P m i n at least k p times for a soft-switching range of k p :1 as (32):
P k p P m i n .
Eventually, the optimization problem is defined as follows:
*) Minimization problem:
Find a set of f s , B p k , L k , M to minimize the total converter loss:
f f s , B p k , L k , M = Δ P t o t m i n
subject to:
g 1 f s , B p k , L k , M = F s , m i n f s 0 g 2 f s , B p k , L k , M = f s F s , m a x 0 g 3 f s , B p k , L k , M = L k L k , m a x 0 g 4 f s , B p k , L k , M = L k , m i n L k 0 g 5 f s , B p k , L k , M = max Δ T t r , Δ T p e Δ T m a x 0 g 6 f s , B p k , L k , M = k p P m i n P 0 ,
where k p : 1 is the power range of the converter, k p should be as great as possible. In this paper, k p is chosen as 5 for a power range of 5:1.
*) Solution:
In order to solve the above optimization problem, there are numerous of available methods. However, note here that the optimizing function is non-continuous and non-differentiable, as it includes a round-up (ceiling) function. Therefore, explicit methods based on gradient are not applicable; instead, heuristic techniques such as the genetic algorithm (GA), differential evolution (DE), particle swarm optimization (PSO), etc., are preferred [17]. In this paper, GA is selected due to its simplicity and feasibility.
In terms of implementation, among many available tools, the Global Optimization Toolbox (GOT) [18] of MATLAB and Simulink appears to be a very good choice, as it can generate the global optima rather than a local one. Besides, its usage is also very simple and does not require deep knowledge on the fitness function. Therefore, the GA function of GOT was employed to solve the above optimization problem. The fitness function was given by (33). The nonlinear inequalities constraints were (34). There were no linear and/or equality constraints. Flow chart diagram of the whole optimization and design process is given in Figure 5.

6. Results and Discussion

6.1. Optimization Results

The proposed design and optimization strategy was applied to design a 10 kW, three-phase dual-active-bridge (DAB) converter for battery charging application. From (2), the primary current can be calculated with a preliminary phase shift angle. Afterward, the current density J is chosen to be, for example, 9 A/mm 2 . Then, the magnetic core of ETD54/28/19 from TDK corporation is selected. The specification and preliminary parameters of the required converter are given in Table 1.
Lower and upper bounds of optimizing variables are given in Table 2. Running the above optimization algorithm, optimization progress is shown in Figure 6. The optimal loss is 210.02 W, achieved at the switching frequency of 72.71 kHz, the peak flux density of 137.6 mT, the leakage inductance of 5.05 mT, and the voltage conversion ratio of 1:1. The best efficiency was obtained as 97.9% at the rated power of 10 kW. The temperature rise of transformer was calculated as 70 degrees—just right—equal to the maximum allowable value given in Table 1. Other parameters are listed in Table 3.
Since 72.71 kHz is not suitable for programming, it was round-up to 75 kHz. This causes a slightly deviation in the total loss and efficiency of the converter compared to the optimal value. However, that deviation is insignificant and can be ignored. The modified optimal parameters and all other design parameters are described in Table 3.
Also listed in Table 3 is the competing parameter set. For the same transformer design, we operated the system at a different switching frequency (and thus, different phase shift), and then compared the performance to confirm whether or not the optimized parameters are better. It was not possible to lower switching frequency, or the control range constraint would have been violated. Therefore, the competing switching frequency was 100 kHz.
Figure 7 illustrates the calculated efficiency and temperature curves obtained by the optimized and competing parameters using the model in the previous sections when varying the phase shift from the minimum value to the nominal one. As shown in Figure 7a, the minimum power was calculated at the minimum phase shift, which is equivalent to the dead-time period. Maximum efficiency of near 98.5% was achieved at around 4 kW. The expected efficiency at the rated power was 97.87%, as mentioned above. The plotted temperature curve indicates the trend of the greater value between temperature of transformers and power electronics devices with the assumption that the ambient temperature is 15 degree Celsius. At the rated power, the hottest spot was about 85 degree Celsius. Note that this calculation assumes there is no forced air cooling for both transformers and semiconductors, only passive cooling. The characteristics are same when considering the competing parameter set. However, the efficiency is expected to be slightly lower while temperature rise is a little bit less than that when applying the optimized parameter set (due to smaller transformer loss), as can be seen in Figure 7b.

6.2. Simulation and Experimental Results

Before conducting experiments, a finite-element analysis (FEA) using the Finite-Element-Method- magnetic (FEMM) version 4.2 [19] was carried out to validate the designed transformers. The flux density obtained by the FEA is illustrated in Figure 8. The flux has higher density at the two outer limbs then at the center one. The peak flux density is 150 mT, which is 12.5% greater than the estimation value. Except for that, other parameters including leakage inductance, AC resistance, and copper and ferrite losses are almost identical to the designed values. Therefore, it is reasonable to fabricate three transformers, namely, Y 1 , Y 2 , and Y 3 , according to the designed parameters. The complete prototype of the converter is shown in Figure 9.
A detailed comparison between calculations, Finite-element analysis (FEA) results, and measurements is given in Table 4. It can be seen that measurements matched calculation and FEA results very well. The simulated leakage inductance was 4.87 μ H, while the expected value was 5.05 μ H, only 3.5% different. The prediction of transformer loss was also very close to the FEA simulation results. However, the simulated peak flux density was about 13% higher than expected. That is because the FEMM simulation is 2D, so two outer limbs of the core are treated as rectangular (while in fact, two outer limbs of the ETD core do not have rectangular shapes), making the effective cross-section area ( A c ) slightly smaller than the real parameter. On the other hand, the center limb which has cylinder shape in deed is also treated as rectangular, and thus, the flux density in the center limb in FEA simulation is smaller than expected. Anyhow, the difference of peak flux density contributes only to the difference on core loss, which is not significant in this design. Since the transformer designed can be confirmed by FEA simulation, three single phase transformers were built accordingly. However, there are small differences between the actual transformers. This distinction may lead to unevenly loss and heat distributions among phases. It is also worth noting that the measured AC resistances were slightly greater than expected due to the non-ideal winding arrangement. In reality, the actual transformer losses would be slightly higher than estimated. Nevertheless, it can be confirmed here that the converter models and optimization strategy presented in the previous sections are reasonable.
As the transformers were validated, two inverters made of SiC MOSFET C2M0040120D from CREE incorporation were placed at primary and secondary sides. A TMS320F28335 control card is employed to generate pulse-width-modulation (PWM) signals. A programmable power supply connects to the input side of Inverter 1, whereas a DC electronic load connects to the DC side of Inverter 2. The experiment system is depicted in Figure 10.
Dead-time was fixed at 100 ns (equivalent to about 2.7 electrical degrees). The input voltage was set at 380 V. The DC electronics load was configured to operate in the constant resistance mode whose resistance was changed in an open-loop manner to keep the input and output voltage equal. The phase shift was then increased gradually from ψ m i n until reaching the maximum power. At each step, efficiency and temperature were recorded. Drain-source voltage of a primary switch and phase current waveform at 9.3 kW are demonstrated in Figure 11. As shown, soft-switching was achieved; both waveforms are smooth as there are no ringings or spikes.
Efficiency in the whole power range at unity voltage conversion ratio when using the optimized parameter set is shown in Figure 12. Due to the current limitation of the power supply, the maximum power that can be investigated is only 9.3 kW. The maximum efficiency of 98.65% is obtained at around 3.6 kW, and the average efficiency over the whole power range is 97.91%, as shown in Figure 12a. The measured results match well with the calculations in the low and mid-power ranges. When the power increases, the converter performance is, however, not as good as expected. For example, the estimated efficiency at 9.3 kW is 97.9% but the measured result shows only 97.15%; in other words, there is a 0.75% error in the determination of the efficiency. The reason may be due to the increment of drain-source resistance of MOSFETs when temperature changes, which was not regarded in the loss model above. Another reason is the mismatch of actual parameters of the real, hand-made transformers when compared to the ideal values. Voltage drops on cables connecting input and output of the converter to the power supply and load also contribute to the error between modeling and experiment. Note that at 10 kW, a 1 V drop on one cable may lead to a 0.25% reduction in the overall efficiency.
The measured temperature also matches well with calculation in the mid and high-power range, as can be seen in Figure 12c. In the low-power range, the temperature was overestimated in calculations, as the measurement values were about 5 degrees cooler than expected. Notes here that the plotted temperatures are the hottest spots on the converter regardless of transformers or inverters being measured. And note also that there is no forced air cooling on either heatsinks or transformers, only passive cooling. For example, as illustrated in Figure 13, the temperature measurement at 8.2 kW using a FLIR camera shows that the hot spots on the transformer and inverter were 71 and 63 degrees Celsius, respectively, whereas the expected value was 70.7 degrees. These results indicate that the converter loss and heat model, and the optimized strategy, are reliable. Furthermore, if the converter has forced-air cooling, it can perform even better.
It is also worth paying attention to Figure 13b. The loss appears to not be distributed evenly among three legs. The temperature of leg c of the secondary inverter seems to be hotter than the others. That is comprehensible, since the three transformers are not identical, as described in Table 4. A slight difference in transformer parameters leads to uneven loss and heat distribution, as shown in Figure 13b. In order to avoid this phenomenon, better transformer fabrication is preferable. An alternative method is to modify the modulation strategy for evenly current distribution between phases. However, this topic is beyond the scope of this paper.
As expected, when using the competing parameter set, the efficiency was not as high as that when switching the converter at the optimized frequency (see Figure 12b). The best performance was only 98.3%, around 4 kW. When increasing the power, the efficiency drops rapidly. At 6.5 kW, the converter can transfer only 97.2% of the input power, while that value is 97.9% when using the optimized switching frequency. However, the system is cooler, as predicted above, when the measured temperature at 6.5 kW is about 58.9 degrees Celsius compared to 60 degrees of the previous case (Figure 12d). Hence, in terms of efficacy, it can be concluded that the converter with the optimized parameters performs better.

7. Conclusions

This paper presented the detail and accurate loss modeling for three-phase, dual-active-bridge DC/DC converters with inductor-integrated transformers. After that, an optimization strategy was proposed to minimize the converter loss. Several nonlinear inequality constraints were regarded. There are four variables that were optimized, including switching frequency, peak flux density, leakage inductance, and the voltage conversion ratio. Due to the highly nonlinear optimizing function, the genetic algorithm was employed to optimize the converter. Finite-element-analysis showed that the converter model is reliable. Experimental results confirmed that the optimized design of the converter had the maximum efficiency of 98.65% and average efficiency in the power range of 97.91%.

Author Contributions

D.-D.N. proposed the main idea, built and performed experiments; N.-T.B. designed the Genetic Algorithm; K.Y. provides supporting ideas, comments and funding; all members contributed significantly in writing and revising the paper. All authors have read and agreed to the published version of the manuscript.

Funding

Knowledge base Aichi Research and development base that supports high value-added manufacturing-Project ID: PI6.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kempton, W.; Tomić, J. Vehicle-to-grid power implementation: From stabilizing the grid to supporting large-scale renewable energy. J. Power Sources 2005, 144, 280–294. [Google Scholar] [CrossRef]
  2. Ota, Y.; Taniguchi, H.; Nakajima, T.; Liyanage, K.M.; Baba, J.; Yokoyama, A. Autonomous distributed V2G (vehicle-to-grid) satisfying scheduled charging. IEEE Trans. Smart Grid 2011, 3, 559–564. [Google Scholar] [CrossRef]
  3. Uddin, K.; Moore, A.D.; Barai, A.; Marco, J. The effects of high frequency current ripple on electric vehicle battery performance. Appl. Energy 2016, 178, 142–154. [Google Scholar] [CrossRef] [Green Version]
  4. He, P.; Khaligh, A. Comprehensive analyses and comparison of 1 kW isolated DC–DC converters for bidirectional EV charging systems. IEEE Trans. Transp. Electrif. 2016, 3, 147–156. [Google Scholar] [CrossRef]
  5. Zhu, L.; Taylor, A.R.; Liu, G.; Bai, K. A Multiple-Phase-Shift Control for a SiC-Based EV Charger to Optimize the Light-Load Efficiency, Current Stress, and Power Quality. IEEE J. Emerg. Sel. Top. Power Electron. 2018, 6, 2262–2272. [Google Scholar] [CrossRef]
  6. Haghbin, S.; Alatalo, M.; Yazdani, F.; Thiringer, T.; Karlsson, R. The design and construction of transformers for a 50 kW three-phase dual active bridge dc/dc converter. In Proceedings of the 2017 IEEE Vehicle Power and Propulsion Conference (VPPC), Belfort, France, 11–14 December 2017; pp. 1–5. [Google Scholar]
  7. van Hoek, H.; Neubert, M.; De Doncker, R.W. Enhanced modulation strategy for a three-phase dual active bridge—Boosting efficiency of an electric vehicle converter. IEEE Trans. Power Electron. 2013, 28, 5499–5507. [Google Scholar] [CrossRef]
  8. Baars, N.H.; Everts, J.; Wijnands, C.G.; Lomonova, E.A. Performance evaluation of a three-phase dual active bridge dc–dc converter with different transformer winding configurations. IEEE Trans. Power Electron. 2016, 31, 6814–6823. [Google Scholar] [CrossRef]
  9. Segaran, D.; Holmes, D.; McGrath, B. Comparative analysis of single-and three-phase dual active bridge bidirectional dc-dc converters. Aust. J. Electr. Electron. Eng. 2009, 6, 329–337. [Google Scholar] [CrossRef]
  10. Mainali, K.; Tripathi, A.; Madhusoodhanan, S.; Kadavelugu, A.; Patel, D.; Hazra, S.; Hatua, K.; Bhattacharya, S. A Transformerless Intelligent Power Substation: A three-phase SST enabled by a 15-kV SiC IGBT. IEEE Power Electron. Mag. 2015, 2, 31–43. [Google Scholar] [CrossRef]
  11. Iyer, V.M.; Gulur, S.; Bhattacharya, S. Optimal design methodology for dual active bridge converter under wide voltage variation. In Proceedings of the 2017 IEEE Transportation Electrification Conference and Expo (ITEC), Chicago, IL, USA, 22–24 June 2017; pp. 413–420. [Google Scholar]
  12. Wojda, R.; Kazimierczuk, M. Winding resistance of litz-wire and multi-strand inductors. IET Power Electron. 2012, 5, 257–268. [Google Scholar] [CrossRef]
  13. Qin, H.; Kimball, J.W.; Venayagamoorthy, G.K. Particle swarm optimization of high-frequency transformer. In Proceedings of the IECON 2010-36th Annual Conference on IEEE Industrial Electronics Society, Glendale, AZ, USA, 7–10 November 2010; pp. 2914–2919. [Google Scholar]
  14. Jauch, F.; Biela, J. Generalized modeling and optimization of a bidirectional dual active bridge DC-DC converter including frequency variation. IEEJ J. Ind. Appl. 2015, 4, 593–601. [Google Scholar] [CrossRef] [Green Version]
  15. Nguyen, D.D.; Kazuto, Y.; Katou, A.; Yoshida, S. Design optimization of a three-phase Dual-Active-Bridge converter for charging stations. In IEEE Vehicle Power and Propulsion Conference—VPPC 2019; IEEE: Hanoi, Vietnam, 2019; in process. [Google Scholar]
  16. McLyman, C.W.T. Transformer and Inductor Design Handbook; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  17. Kachitvichyanukul, V. Comparison of three evolutionary algorithms: GA, PSO, and DE. Ind. Eng. Manag. Syst. 2012, 11, 215–223. [Google Scholar] [CrossRef] [Green Version]
  18. Global Optimization Toolbox. Available online: https://www.mathworks.com/products/global-optimization.html (accessed on 20 November 2019).
  19. Finite Element Method Magnetics 4.2. Available online: http://www.femm.info/wiki/HomePage (accessed on 20 November 2019).
Figure 1. Three-phase dual active bridge converter topology.
Figure 1. Three-phase dual active bridge converter topology.
Energies 13 00150 g001
Figure 2. Three-phase dual active bridge converter topology.
Figure 2. Three-phase dual active bridge converter topology.
Energies 13 00150 g002
Figure 3. Theoretical voltage and current waveforms.
Figure 3. Theoretical voltage and current waveforms.
Energies 13 00150 g003
Figure 4. Transformer with thick insulation configuration.
Figure 4. Transformer with thick insulation configuration.
Energies 13 00150 g004
Figure 5. Design and optimization flow chart.
Figure 5. Design and optimization flow chart.
Energies 13 00150 g005
Figure 6. Optimization process using genetic algorithm.
Figure 6. Optimization process using genetic algorithm.
Energies 13 00150 g006
Figure 7. Calculated efficiency and temperature obtained by: (a) optimized parameters; (b) competing parameters.
Figure 7. Calculated efficiency and temperature obtained by: (a) optimized parameters; (b) competing parameters.
Energies 13 00150 g007
Figure 8. Finite-element-analysis using FEMM 4.2.
Figure 8. Finite-element-analysis using FEMM 4.2.
Energies 13 00150 g008
Figure 9. A 10 kW, three-phase DAB converter prototype.
Figure 9. A 10 kW, three-phase DAB converter prototype.
Energies 13 00150 g009
Figure 10. A 10 kW, three-phase DAB converter prototype.
Figure 10. A 10 kW, three-phase DAB converter prototype.
Energies 13 00150 g010
Figure 11. Voltage and current waveforms at 9.3 kW.
Figure 11. Voltage and current waveforms at 9.3 kW.
Energies 13 00150 g011
Figure 12. Performance in the soft-switching range. (a) Efficiency when using the optimal parameter set; (b) efficiency when using the competing parameter set; (c) temperature when using the optimal parameter set; (d) temperature when using the competing parameter set.
Figure 12. Performance in the soft-switching range. (a) Efficiency when using the optimal parameter set; (b) efficiency when using the competing parameter set; (c) temperature when using the optimal parameter set; (d) temperature when using the competing parameter set.
Energies 13 00150 g012
Figure 13. Temperature at 9.3 kW: (a) transformer and (b) secondary inverter.
Figure 13. Temperature at 9.3 kW: (a) transformer and (b) secondary inverter.
Energies 13 00150 g013
Table 1. Specification and preliminary parameters.
Table 1. Specification and preliminary parameters.
ParameterSymbolValue
Input voltage V 1 380 V
Output voltage V 2 240−480 V
Nominal output voltage V 2 , n o m 380 V
Rated powerP10 kW
Maximum temperature rise Δ T 70 C
Preliminary frequency f s 50 kHz
Preliminary phase shift ψ 30 degrees
Preliminary peak flux density B p k 250 mT
Calculated inductance L k 15.56  μ H
Desired current densityJ9 A/mm 2
Selected wire Litz 1050AWG44
Selected ferrite core ETD54/28/19 (N87)
SiC MOSFET CREE C2M0040120D
Table 2. Lower and upper bounds for genetic algorithm.
Table 2. Lower and upper bounds for genetic algorithm.
ParameterSymbolLower BoundUpper Bound
Switching frequency f s 50 kHz300 kHz
Peak flux density B p k 50 mT250 mT
Leakage inductance L k μ H10  μ H
Voltage conversion ratioM0.52
Table 3. Optimal and modified parameters.
Table 3. Optimal and modified parameters.
ParameterSymbolOptimized ValueModified ValueCompeting ValueUnit
Switching frequency f s 72.7175100kHz
Peak flux density B p k 137.6133.1397.68mT
Leakage inductance L k 5.055.055.05 μ H
Conversion ratioM111
Number of turns N 1 : N 2 15:1515:1515:15
Phase shift angle ψ 13.113.5418.46degrees
AC resistance R A C 31.8631.9933.61m Ω
Transformer temperature rise Δ T t r 70.069.7569.15 C
Heatsink temperature rise Δ T p e 40.741.4449.73 C
Power electronics loss Δ P p e 162.80165.76198.93W
Transformer loss Δ P t r 42.9742.7842.34W
Total loss Δ P t o t 210.02212.94247.24W
Efficiency@ P r a t e d η 97.9%97.87%97.53%
Table 4. Design validation.
Table 4. Design validation.
ParameterSymbolCalculationFEAMeasurementUnit
Y 1 Y 2 Y 3
Leakage inductance L k 5.054.874.924.94.82 μ H
Magnetizing inductance L m 1.21.171.331.131.24mH
AC resistance @ f s , o p t R A C 31.9930.8136.536.535.2m Ω
Peak flux density B p k 133.13149.9---mT
Copper loss Δ P C u 10.805710.19---W
Core loss Δ P F e 3.45553.03---W

Share and Cite

MDPI and ACS Style

Nguyen, D.-D.; Bui, N.-T.; Yukita, K. Design and Optimization of Three-Phase Dual-Active-Bridge Converters for Electric Vehicle Charging Stations. Energies 2020, 13, 150. https://doi.org/10.3390/en13010150

AMA Style

Nguyen D-D, Bui N-T, Yukita K. Design and Optimization of Three-Phase Dual-Active-Bridge Converters for Electric Vehicle Charging Stations. Energies. 2020; 13(1):150. https://doi.org/10.3390/en13010150

Chicago/Turabian Style

Nguyen, Duy-Dinh, Ngoc-Tam Bui, and Kazuto Yukita. 2020. "Design and Optimization of Three-Phase Dual-Active-Bridge Converters for Electric Vehicle Charging Stations" Energies 13, no. 1: 150. https://doi.org/10.3390/en13010150

APA Style

Nguyen, D. -D., Bui, N. -T., & Yukita, K. (2020). Design and Optimization of Three-Phase Dual-Active-Bridge Converters for Electric Vehicle Charging Stations. Energies, 13(1), 150. https://doi.org/10.3390/en13010150

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