Multileg Interleaved Buck Converter for EV Charging: Discrete ‐ Time Model and Direct Control Design

: This paper presents the modeling and the implementation of the digital control of a multileg interleaved DC ‐ DC buck converter for electrical vehicle (EV) charging. Firstly, we derive a discrete averaged model of an n ‐ leg interleaved buck converter (IBC). Secondly, we present a direct tuning procedure for one primary discrete PIDF (PID + filter) and multiple secondary PI controller. The objective of the control system is to regulate the current flow in each leg of the converter. This task is accomplished by introducing a novel control paradigm that simultaneously addresses two aims: on the one hand, the control scheme must guarantee an acceptable level of robustness under load variations; while on the other, an even distribution of power on each leg must be ensured at any operational condition. The proposed strategy hinges on a technique that combines simplicity and precision in the fulfillment of design frequency specifications. We use simulations and a digital signal processor (DSP) based experimental implementation of the design technique to validate the proposed methodology.


Introduction
Electric vehicle (EV) technologies are facing huge developments in the last years because of their potential to reverse the carbon emissions trend and therefore to lead to a green shift in the coming years [1]. However, EV manufacturers are still addressing relevant problems such as charging times and users' range anxiety. In the field of power electronics, these targets can be achieved mainly by enhancing the performance of the power converter, and at the same time, by keeping the cost as low as possible. Electric vehicle charging facilities provided by fast charger feeders are expected to be one of the most significant players in building a reliable and efficient charging network to significantly increase range capabilities through a sharp reduction of current EVs' charging times [2]. Several configurations such as the Vienna rectifier [3], multilevel neutral point clamped (NPC) choppers [4], and interleaved converters have been proposed in the literature. Among them, the interleaved converter topology appears to fit the economic and technical constraints mentioned above. Because of its inherent modular structure, it is possible to take advantage of the low cost and highly standardized modules, while at the same time increasing the performance in terms of carried power and electric power quality. Interleaved topology is used for both AC/DC [5][6][7] and DC/DC [8,9] power converters. The DC/DC buck is the last active stage before the battery. It is, therefore, crucial to keep the output current (in terms of current ripple) within battery constraints. For this reason, this paper focuses on the chopper control system. One of the most well-known benefits of interleaved topology is the spontaneous output current ripple aforementioned frequency domain specifications to be exactly satisfied, thus ensuring not only the satisfaction of the steady-state performance, but also the shaping of the transient (avoidance of overshoot/undershoot, velocity of the convergence, etc.) and, at the same time, guarantees an acceptable level of robustness under load variations during battery charge.
The paper is organized as follows. The continuous and discrete-time models of standard multileg interleaved buck converter are presented in Section 2. In Section 3, we propose the control structure and PI/PIDF tuning procedure. Simulation results and analysis in comparison with the method introduced in [13] and results of DSP-based experimental implementation are presented in Section 4. Conclusions and possible future works finally end the paper.

Continuous and Discrete-Time Model of n-Leg Interleaved Buck Converter
Let us consider the n-leg interleaved buck converter shown in Figure 1. In this scheme, V in denotes the input voltage, V o represents the output voltage, while L 1 , L 2 , . . . , L n and i 1 , i 2 , . . . , i n denote the leg inductances and the corresponding currents, respectively. Here, C is the filter capacitance introduced to filter the high-frequency ripple produced by the switching operation, while R L and R C denote the parasitic series resistances of the inductors and capacitor, respectively. The load represents the battery to be charged. The n duty cycle signals d 1 , d 2 , . . . , d n modulate the PWM signals to drive the buck converter switches according to the proposed control design. Using the Kirchhoff's current law (KCL), the total current i t (t) can be written as: The model of the buck converter is obtained using the average modeling technique described in [14,24], while the droop control design method is used to model the battery as a power load R [13,25]. This simplified model, compared to other dynamical models, such as Thévenin or non-linear models [26], is generally used in the design of the controller since the product of its voltage and current vary very slowly with respect to the converter dynamics [27]. In particular, the load droop coefficient R is selected by considering the allowable battery voltage range ∆V o and the constant reference current I ref , which is R = ∆V o /I ref . In this contest, the state space equations of the interleaved converter are where and R Swk represents k-leg power switch conduction losses. It is assumed that the value of R Sw is equal in both switches of the same leg. and RSwk represents k-leg power switch conduction losses. It is assumed that the value of RSw is equal in both switches of the same leg. Neglecting the effects introduced by the capacitor resistance RC, the output voltage can be approximated to VC, while the output current iout(t) can be approximated to the total current it(t). In this way, the total current can be obtained by considering the summation of inductors' currents ik(t) only, and no additional sensor is required.
Let us now define the average duty cycle as .. 2 1 , and the average model transfer function of the n-leg interleaved buck converter is The Bode diagrams of the frequency response G (jω) when n is equal to 3, 6, and 9 are shown in Figure 2. Neglecting the effects introduced by the capacitor resistance R C , the output voltage can be approximated to V C , while the output current i out (t) can be approximated to the total current i t (t). In this way, the total current can be obtained by considering the summation of inductors' currents i k (t) only, and no additional sensor is required.
Let us now define the average duty cycle as Assuming the average model transfer function of the n-leg interleaved buck converter is where The Bode diagrams of the frequency response G (jω) when n is equal to 3, 6, and 9 are shown in Figure 2.  The transfer function of the discrete model of the n-leg interleaved buck converter is given by where a = e −ξω n T s cos ω n T s 1 − The transfer function (7) can be obtained by applying the definition of the Z-transform to the product of the continuous-time transfer function of the converter G(s) and the transfer function of the zero-order hold Notice that G(z) is characterized by two complex conjugate poles.
which gives rise to well-documented detrimental oscillatory effects, see [28].

The Proposed Control
The battery charging profile for general battery chargers was described in [16]. After a pre-charging mode, during which the output current of the converter gradually ramps up with a staircase wave current control, a constant current control leads to a fast charge up to the maximum battery voltage. Then, the charging current gradually decreases under constant voltage control up to the final value of the output current, see Figure 3.  The proposed control strategy for the interleaved buck converter topology shown in Figure 1 can be used both for pre-charging mode control and for constant current control. The main idea is to address the total output control problem and the inductor current sharing control separately. In particular, the average duty cycle dt(z) is computed by a discrete controller to track the total current reference signal, while leg duty cycles are computed by separate digital controllers using the current balancing technique as in [13].

Average Duty Cycle Control
Let us consider the discrete control scheme shown in Figure 4, where G(s) and G(z) represent the buck converter transfer Functions (5) and (7), respectively. Moreover, r(z) is the discrete reference The proposed control strategy for the interleaved buck converter topology shown in Figure 1 can be used both for pre-charging mode control and for constant current control. The main idea is to address the total output control problem and the inductor current sharing control separately. In particular, the average duty cycle d t (z) is computed by a discrete controller to track the total current reference signal, while leg duty cycles are computed by separate digital controllers using the current balancing technique as in [13].

Average Duty Cycle Control
Let us consider the discrete control scheme shown in Figure 4, where G(s) and G(z) represent the buck converter transfer Functions (5) and (7), respectively. Moreover, r(z) is the discrete reference signal, while i t (z) denotes the sampled and hold current to be controlled. The value of the average duty cycle d t (z) is generated by the compensator C(z) to track the error signal e(z) = r(z) − i t (z). In addition, let L(z) denote the loop-gain transfer function, i.e., L(z) = C(z) G(z).
can be used both for pre-charging mode control and for constant current control. The main idea is to address the total output control problem and the inductor current sharing control separately. In particular, the average duty cycle dt(z) is computed by a discrete controller to track the total current reference signal, while leg duty cycles are computed by separate digital controllers using the current balancing technique as in [13].

Average Duty Cycle Control
Let us consider the discrete control scheme shown in Figure 4, where G(s) and G(z) represent the buck converter transfer Functions (5) and (7), respectively. Moreover, r(z) is the discrete reference signal, while it(z) denotes the sampled and hold current to be controlled. The value of the average duty cycle dt(z) is generated by the compensator C(z) to track the error signal e(z) = r(z) − it(z). In addition, let L(z) denote the loop-gain transfer function, i.e., L(z) = C(z) G(z). The proposed discrete-time controller is the following biquadratic filter which can be directly implemented on a microcontroller. With the values the controller (10) can be written in the equivalent form The proposed discrete-time controller is the following biquadratic filter which can be directly implemented on a microcontroller. With the values the controller (10) can be written in the equivalent form with complex conjugate zeros. The proposed design procedure for the tuning of the parameters a 1 , a 2 , b 0 , b 1 , b 2 is based on the conjunction of the classical pole/zero cancellation and the so-called inversion formulae methods introduced in [25] for a standard single-leg buck converter. In this paper, the design procedure is summarized by giving the logical steps of the algorithm for the solution of the following design problem.

Control Problem
Find the values of the parameters a 1 , a 2 , b 0 , b 1 , b 2 to guarantee zero position error, and to assign the phase margin Φ m and the gain crossover frequency ω g of the open-loop frequency response L (e jωTs ), where T s is the sampling period: Step 1: Calculate the values of ω n , ω o , ξ using (6) on the basis of the parameters of the circuit and determine the discrete transfer function G(z) using (7) and (8).
Step 2: Compute the values of ω d , δ d by placing the zeros of (12) in the same location of the complex conjugate poles of G(z) to cancel their detrimental effects using.
Step 3: Evaluate the magnitude and the phase of the frequency response G(e j2ωTs ) of the plant multiplied by the frequency response of the factor of the controller (12) that has already been determined Energies 2020, 13, 466 7 of 18 at the desired gain crossover frequency ω g .
Step 4: Compute the magnitude and the phase that the controller should introduce to exactly satisfy the given specification on the phase margin Φ m .
Step 5: Calculate the values of the remaining degrees of freedom of the controller (12) using the inversion formulae.
The given control problem has a feasible solution if and only if β d > 0 and K i > 0, see Remark 4.1 in [25].
Step 6: Compute the parameters a 1 , a 2 , b 0 , b 1 , b 2 using (11) to write the controller in the form (10) which is directly implementable on a microcontroller board.

Circulating Current Control
For the minimization of the ripple in the output current, the delay angle of the PWM carriers of each buck converter leg is 360 • /n. An example of inductor currents in 3-leg interleaved buck converter is shown in Figure 5.

Circulating Current Control
For the minimization of the ripple in the output current, the delay angle of the PWM carriers of each buck converter leg is 360°/n. An example of inductor currents in 3-leg interleaved buck converter is shown in Figure 5. It is well known that applying the average duty cycle dt(z) to each PWM is not enough to guarantee a proper operation when the converter is exposed to variations of the inductor parameter. It follows that other control loops have to be considered in order to minimize the differences of the inductor currents both at steady-state and during the transient. For this purpose, let us consider the control system shown in Figure 6.
In particular, the signals to be controlled are It is well known that applying the average duty cycle d t (z) to each PWM is not enough to guarantee a proper operation when the converter is exposed to variations of the inductor parameter. It follows that other control loops have to be considered in order to minimize the differences of the inductor currents both at steady-state and during the transient. For this purpose, let us consider the control system shown in Figure 6. In particular, the signals to be controlled are The control signals denote the differences of duty cycles of the buck converter in each considered mesh. The transfer function of the discrete plant to be controlled is It can be proved that

L Ts
When R s can be neglected, we obtain For the control problem considered in this paper, the following form of a discrete proportional-integral (PI) compensator can be taken into account: The values of parameters K P and K I can be selected to exactly satisfy design specifications on the gain crossover frequency ω G and on the phase margin Φ M , while the specification of zero position error is automatically fulfilled by the structure of the controller. The analytical solution of the control problem can be obtained using the following inversion formulae: where The pulse width modulation (PWM) duty cycles d 1 , d 2 , . . . , d n can be computed solving Equations (3) and (19), i.e.,

Numerical and Experimental Results
Let us consider the parameter configurations of the buck converter shown in Table 1. Case (a) refers to the simulated control system used to validate the proposed model and to evaluate the effectiveness of the proposed control procedure in comparison with the method described in [13]. For this reason, the parameters a) were selected to be equal to the parameters considered in [13] (in particular, see Table III. Case (b) refers to the parameters of the hardware device used for the experimental validation. According to the limited power capabilities of the devices available in our laboratory, the input nominal voltage and the output reference current were scaled down by a factor of 7 and 10, respectively. The A/D converter resolution (12 bits) and the sampling/switching frequency (3 × 20 kHz) were maintained the same in order to have a direct correspondence with the control system implemented in simulations (Case a).

Proposed Control Procedure
The discrete models (7) and (8) of the converter with the parameters a) in Table 1 is The average duty cycle control shown in Figure 4 can be designed to yield zero position error, a phase margin Φ m and a gain crossover frequency ω g of the open-loop frequency response equal to 80 • and 3000 rad/s, respectively.
Using the proposed procedure, the complex poles of G(z) can be cancelled by selecting δ d = 0.934 and ω d = 0.87 rad/s in (12). The gain and the phase that the controller (12) has to introduce at frequency ω g to satisfy the design specification on the phase margin are M g = 0.0023 and ϕ g = 339.6 • , respectively. Using (17) and (18) the remaining parameters of the PIDF controller (12) are β d = 1.01, K i = 3.346 × 10 −4 . The discrete-time controller (10) thus obtained is which follows directly from (11). The step responses of the output current (0-125 A) of the open and closed-loop system are shown in Figure 7. Notice that the steady-state error of the uncompensated system was reduced to zero and the overshoot has been eliminated by the control. Also, the given dynamical control specifications were exactly satisfied, as highlighted by the Nyquist plot and the Bode diagrams of the frequency response of the controlled system shown in Figures 8 and 9.
Regarding the design of the discrete controller for the circulating current control, the zero position error was automatically achieved by the structure of the controller (20), while the specifications on the gain crossover frequency ω G = 8000 rad/s and phase margin Φ M = 50 • can be exactly met using the presented design procedure. The gain and the phase that the controller (20) has to introduce at frequency ω G to satisfy the given dynamic specifications are M G = 0.0044 and ϕ G = 323.8 • , respectively. From (21), the controller (20) that solves the problem is given by The Nyquist plot of the open-loop frequency response of the control system is shown in Figure 10. Notice that the given specifications are exactly met.
Energies 2020, 13, x FOR PEER REVIEW 11 of 19 to introduce at frequency ωG to satisfy the given dynamic specifications are MG = 0.0044 and φG = 323.8°, respectively. From (21), the controller (20) that solves the problem is given by

Numerical Comparison
The proposed control system has been simulated using Matlab R2018a and Simulink 9.1 (solver: ode45), and compared with the method proposed in [13] by adopting the same system parameters

Numerical Comparison
The proposed control system has been simulated using Matlab R2018a and Simulink 9.1 (solver: ode45), and compared with the method proposed in [13] by adopting the same system parameters (see Table III in [13]). In particular, all the blocks of the interleaved buck converter have been selected

Numerical Comparison
The proposed control system has been simulated using Matlab R2018a and Simulink 9.1 (solver: ode45), and compared with the method proposed in [13] by adopting the same system parameters (see Table III in [13]). In particular, all the blocks of the interleaved buck converter have been selected from the Simscape 4.4 Electrical toolset. The 3-leg buck converter has been controlled in [13] with a PI regulator to provide a phase margin of 71.3 • and a gain margin of 4.83 dB in the continuous-time domain. The corresponding transfer function in the z domain is, see Equation (24) in [13] C PI (z) = 1.219 · 10 −3 z + 433.4601 · 10 −6 z − 1 Notice that the resulting phase and gain margin of the loop-gain transfer function were reduced to 18 • and 6.53 dB, respectively. The consequent degradation in dynamic performance was due to the use of an indirect design procedure carried out in the continuous-time domain and was avoided with the direct method proposed in this paper. The step responses of the inductor currents and of the total current obtained with the method described in [13] and the method presented in our paper are shown in Figure 11. The significant improvement of the dynamic behavior was the elimination of the oscillation during the transient, as well as the reduction of the settling time. In both cases, the total current flowing to the load presents a significant reduction of the ripple compared to the inductor currents. Moreover, the overlap of the total current with the output signal of the discrete system (22) confirmed the accuracy of the presented discrete model. [ ] were reduced to 18° and 6.53 dB, respectively. The consequent degradation in dynamic performance was due to the use of an indirect design procedure carried out in the continuous-time domain and was avoided with the direct method proposed in this paper.
The step responses of the inductor currents and of the total current obtained with the method described in [13] and the method presented in our paper are shown in Figure 11. The significant improvement of the dynamic behavior was the elimination of the oscillation during the transient, as well as the reduction of the settling time. In both cases, the total current flowing to the load presents a significant reduction of the ripple compared to the inductor currents. Moreover, the overlap of the total current with the output signal of the discrete system (22) confirmed the accuracy of the presented discrete model. The robustness to parameter variations and model uncertainties were considered to take into account discrepancies between the system model and the physical system. For this reason, the behavior of the inductor and load currents under 20% output load resistance variation was simulated, analyzed, and compared with the method presented in [13]. The results are shown in Figure 12. While the load variation leads to quasi-persistent oscillations using the control described in [13], the considered variation was promptly stabilized in less than 1 ms using the method described in Section 3. Notice that the average value of the currents returns to the setpoint value, while the current ripple changes its value at steady-state due to the consequent duty cycle variation. The robustness to parameter variations and model uncertainties were considered to take into account discrepancies between the system model and the physical system. For this reason, the behavior of the inductor and load currents under 20% output load resistance variation was simulated, analyzed, Energies 2020, 13, 466 13 of 18 and compared with the method presented in [13]. The results are shown in Figure 12. While the load variation leads to quasi-persistent oscillations using the control described in [13], the considered variation was promptly stabilized in less than 1 ms using the method described in Section 3. Notice that the average value of the currents returns to the setpoint value, while the current ripple changes its value at steady-state due to the consequent duty cycle variation. Moreover, the designed control was tested under 50% inductance parameter variation and compared with the method presented in [13], see the corresponding step responses in Figure 13. Both controls compensated the induction variation and tracked the expected steady-state current values. However, the proposed procedure led to a significant improvement in the transient.   Moreover, the designed control was tested under 50% inductance parameter variation and compared with the method presented in [13], see the corresponding step responses in Figure 13. Both controls compensated the induction variation and tracked the expected steady-state current values. However, the proposed procedure led to a significant improvement in the transient. Moreover, the designed control was tested under 50% inductance parameter variation and compared with the method presented in [13], see the corresponding step responses in Figure 13. Both controls compensated the induction variation and tracked the expected steady-state current values. However, the proposed procedure led to a significant improvement in the transient.    Finally, Figure 14 shows the total and the inductor currents under 10% V in variation. Even in this case, the proposed control promptly compensates the input voltage variation reaching a new steady-state condition. In particular, the settling time has been reduced from 1.5 ms to 0.5 ms with respect to the method presented in [13].
Energies 2020, 13, x FOR PEER REVIEW 15 of 19 Finally, Figure 14 shows the total and the inductor currents under 10% Vin variation. Even in this case, the proposed control promptly compensates the input voltage variation reaching a new steadystate condition. In particular, the settling time has been reduced from 1.5 ms to 0.5 ms with respect to the method presented in [13]. Regarding the battery charging profile, see Figure 3, other reference current profiles have been considered. Output and inductor currents under reference current variation from 0 A to 120 A within 10 A steps are shown in Figure 15. This behavior simulated the pre-charging mode control of the battery charging profile, showing that the presented control structure is suitable also with this charging profile. Moreover, Figure 16 shows the rising and falling system step response from 0 A to 125 A and from 125 A to 5 A.  Regarding the battery charging profile, see Figure 3, other reference current profiles have been considered. Output and inductor currents under reference current variation from 0 A to 120 A within 10 A steps are shown in Figure 15. This behavior simulated the pre-charging mode control of the battery charging profile, showing that the presented control structure is suitable also with this charging profile. Moreover, Figure 16 shows the rising and falling system step response from 0 A to 125 A and from 125 A to 5 A. Finally, Figure 14 shows the total and the inductor currents under 10% Vin variation. Even in this case, the proposed control promptly compensates the input voltage variation reaching a new steadystate condition. In particular, the settling time has been reduced from 1.5 ms to 0.5 ms with respect to the method presented in [13]. Regarding the battery charging profile, see Figure 3, other reference current profiles have been considered. Output and inductor currents under reference current variation from 0 A to 120 A within 10 A steps are shown in Figure 15. This behavior simulated the pre-charging mode control of the battery charging profile, showing that the presented control structure is suitable also with this charging profile. Moreover, Figure 16 shows the rising and falling system step response from 0 A to 125 A and from 125 A to 5 A.

Experimental Results
The experimental tests were carried out using a small-scale hardware setup composed by a three-phase module employed in an interleaved configuration having parameters shown in Table 1 as "Case (b)". In particular, the Mitsubishi PS22A79 intelligent power IGBT module (1200 V, 50 A) was driven by a Texas Instruments Delfino F28379D DSP control card running the TMS320F28379D microcontroller unit (MCU) via optical interface links ( Figure 17). Inductor currents were sensed employing Hall effect sensors LEM LA 55-P (50 A) and forwarded to the ADC inputs of the DSP board. The control scheme visible in Figure 6 (light blue background) was implemented on Simulink 9.1 and executable C code was generated using Simulink Coder, Embedded Coder, and Texas Instruments Support from Embedded Coder. Finally, employing Code Composer Studio interface, the generated code was deployed to the DSP target. In particular, the DSP board was in charge of running the proposed control method considering sampled currents as input digital signals and PWM waveforms as output signals.
The step response from 0 A to 10 A and zoom at the steady-state of the inductor current sensor outputs and output voltage are shown in Figure 18. In particular, the output current, which can be calculated by the output voltage over resistive load, reached the setpoint value in 5 ms, without oscillations. The zoom of the output signal of the inductor current sensor and the output voltage at steady-state highlights the PI control effect in the balancing of the inductor current.

Experimental Results
The experimental tests were carried out using a small-scale hardware setup composed by a three-phase module employed in an interleaved configuration having parameters shown in Table 1 as "Case (b)". In particular, the Mitsubishi PS22A79 intelligent power IGBT module (1200 V, 50 A) was driven by a Texas Instruments Delfino F28379D DSP control card running the TMS320F28379D microcontroller unit (MCU) via optical interface links ( Figure 17). Inductor currents were sensed employing Hall effect sensors LEM LA 55-P (50 A) and forwarded to the ADC inputs of the DSP board. The control scheme visible in Figure 6 (light blue background) was implemented on Simulink 9.1 and executable C code was generated using Simulink Coder, Embedded Coder, and Texas Instruments Support from Embedded Coder. Finally, employing Code Composer Studio interface, the generated code was deployed to the DSP target. In particular, the DSP board was in charge of running the proposed control method considering sampled currents as input digital signals and PWM waveforms as output signals.
The step response from 0 A to 10 A and zoom at the steady-state of the inductor current sensor outputs and output voltage are shown in Figure 18. In particular, the output current, which can be calculated by the output voltage over resistive load, reached the setpoint value in 5 ms, without oscillations. The zoom of the output signal of the inductor current sensor and the output voltage at steady-state highlights the PI control effect in the balancing of the inductor current.

Experimental Results
The experimental tests were carried out using a small-scale hardware setup composed by a three-phase module employed in an interleaved configuration having parameters shown in Table 1 as "Case (b)". In particular, the Mitsubishi PS22A79 intelligent power IGBT module (1200 V, 50 A) was driven by a Texas Instruments Delfino F28379D DSP control card running the TMS320F28379D microcontroller unit (MCU) via optical interface links ( Figure 17). Inductor currents were sensed employing Hall effect sensors LEM LA 55-P (50 A) and forwarded to the ADC inputs of the DSP board. The control scheme visible in Figure 6 (light blue background) was implemented on Simulink 9.1 and executable C code was generated using Simulink Coder, Embedded Coder, and Texas Instruments Support from Embedded Coder. Finally, employing Code Composer Studio interface, the generated code was deployed to the DSP target. In particular, the DSP board was in charge of running the proposed control method considering sampled currents as input digital signals and PWM waveforms as output signals.
The step response from 0 A to 10 A and zoom at the steady-state of the inductor current sensor outputs and output voltage are shown in Figure 18. In particular, the output current, which can be calculated by the output voltage over resistive load, reached the setpoint value in 5 ms, without oscillations. The zoom of the output signal of the inductor current sensor and the output voltage at steady-state highlights the PI control effect in the balancing of the inductor current. Step response from 0A to 10A (a) and zoom in steady-state condition (b) of inductor currents (current transducer outputs, blue, yellow, magenta) and output voltage (green).
The outputs of the current sensors and the output voltage under 20% output load variation, i.e., from 4.75 Ω to 5.94 Ω and from 5.94 Ω to 4.75 Ω, are shown in Figure 19a,b, respectively. Despite the load perturbation, the total current it(t) reaches the reference value (10 A), while the inductor currents balance each other since they have the same average value in both cases.

Conclusions
In this work, we introduced a new design/tuning procedure for the control of an interleaved buck converter for EV charging. The challenge to guarantee good dynamic performances, even current share among the legs, as well as an acceptable level of robustness under load variations, made rules of thumb and trial-and-error methods inapplicable to the control of these converters. The approach that we proposed overcomes this problem since it directly and exactly satisfied standard design constraints such as gain crossover frequency, phase margin, and zero position error for the PIDF main controller and circulating current PIs. In order to prove the effectiveness of the proposed approach, we compared our strategy with the method introduced in [13]. Moreover, taking advantage of the generic formulation of the discrete-time model and the circulating current control, it was possible to quickly scale the system on interleaved topology employing a higher number of legs staked together. Experimental results carried out on a laboratory prototype proved the control system capabilities. Since the load resistance (representing the battery) varied during charging Step response from 0A to 10A (a) and zoom in steady-state condition (b) of inductor currents (current transducer outputs, blue, yellow, magenta) and output voltage (green).
The outputs of the current sensors and the output voltage under 20% output load variation, i.e., from 4.75 Ω to 5.94 Ω and from 5.94 Ω to 4.75 Ω, are shown in Figure 19a,b, respectively. Despite the load perturbation, the total current i t (t) reaches the reference value (10 A), while the inductor currents balance each other since they have the same average value in both cases. (a) (b) Figure 18.
Step response from 0A to 10A (a) and zoom in steady-state condition (b) of inductor currents (current transducer outputs, blue, yellow, magenta) and output voltage (green).
The outputs of the current sensors and the output voltage under 20% output load variation, i.e., from 4.75 Ω to 5.94 Ω and from 5.94 Ω to 4.75 Ω, are shown in Figure 19a,b, respectively. Despite the load perturbation, the total current it(t) reaches the reference value (10 A), while the inductor currents balance each other since they have the same average value in both cases.

Conclusions
In this work, we introduced a new design/tuning procedure for the control of an interleaved buck converter for EV charging. The challenge to guarantee good dynamic performances, even current share among the legs, as well as an acceptable level of robustness under load variations, made rules of thumb and trial-and-error methods inapplicable to the control of these converters. The approach that we proposed overcomes this problem since it directly and exactly satisfied standard design constraints such as gain crossover frequency, phase margin, and zero position error for the PIDF main controller and circulating current PIs. In order to prove the effectiveness of the proposed approach, we compared our strategy with the method introduced in [13]. Moreover, taking advantage of the generic formulation of the discrete-time model and the circulating current control, it was possible to quickly scale the system on interleaved topology employing a higher number of legs staked together. Experimental results carried out on a laboratory prototype proved the control system capabilities. Since the load resistance (representing the battery) varied during charging

Conclusions
In this work, we introduced a new design/tuning procedure for the control of an interleaved buck converter for EV charging. The challenge to guarantee good dynamic performances, even current share among the legs, as well as an acceptable level of robustness under load variations, made rules of thumb and trial-and-error methods inapplicable to the control of these converters. The approach that we proposed overcomes this problem since it directly and exactly satisfied standard design constraints such as gain crossover frequency, phase margin, and zero position error for the PIDF main controller and circulating current PIs. In order to prove the effectiveness of the proposed approach, we compared our strategy with the method introduced in [13]. Moreover, taking advantage of the generic formulation of the discrete-time model and the circulating current control, it was possible to quickly scale the system on interleaved topology employing a higher number of legs staked together. Experimental results carried out on a laboratory prototype proved the control system capabilities.
Since the load resistance (representing the battery) varied during charging procedures, future work might investigate the possibility of improving the versatility of the proposed procedure by employing a real-time measure/estimation of the load. In this way, the parameters of the PIDF controller could be updated by live DSP computations using the inversion formulae described in Section 3 turning the proposed method into an adaptive control system design procedure. Moreover, findings readily available in this paper might be used to directly tune the control system of any third party interleaved buck converter regardless of power, number of modules, switching frequency, coupling reactors, and filters.