Next Article in Journal
Performance Improvement of Polar Codes via UEP Product Coding
Previous Article in Journal
Resolving the Unusual Gate Leakage Currents of Thin-Film Transistors with Single-Walled Carbon-Nanotube-Based Active Layers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Control Scheme to Suppress Circulating Currents in Parallel-Connected Three-Phase Inverters

Grupo de Sistemas Electrónicos Industriales, Departamento de Ingeniería Electrónica, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
*
Author to whom correspondence should be addressed.
Electronics 2022, 11(22), 3720; https://doi.org/10.3390/electronics11223720
Submission received: 6 October 2022 / Revised: 3 November 2022 / Accepted: 10 November 2022 / Published: 13 November 2022

Abstract

:
The parallel operation of inverters has many benefits, such as modularity and redundancy. However, the parallel connection of inverters produces circulating currents that may result in malfunctions of the system. In this work, a control technique for the elimination of the low-frequency components of the circulating currents in grid-connected inverters is presented. The proposed control structure contains n − 1 zero-sequence control loops, with n being the number of inverters connected in parallel. Simulation and experimental results have been carried out on a prototype composed of two 5 kW inverters connected in parallel. The results have been obtained by considering the following mismatches between both inverters: inductance values of the grid filters, unbalance of the delivered power, and the use of different modulation techniques.

1. Introduction

The parallel connection of inverters has been used in various applications, as it has some advantages, such as modularity, redundancy, and easy maintenance [1,2,3]. However, when there are mismatches in the hardware components, such as the inductance values of the filters or the control parameters of the modules connected in parallel, circulating currents happen to lead to undesirable effects in the system [4,5,6]. Circulating currents can produce adverse effects, such as the distortion of the AC currents, an unbalance in the power delivered by each inverter, higher conduction losses in the switching devices, and lower efficiency of the whole system.
One of the simplest methods to eliminate circulating currents uses transformers to electrically isolate the parallel connected modules. This method leads to an expensive and bulky system, with control methods not requiring extra hardware being preferred [7,8,9,10,11].
The circulating current is composed of high-frequency harmonics, close to the switching frequency and its multiples, and low-frequency harmonics placed at the fundamental frequency and its multiples. High-frequency circulating currents are usually eliminated using passive elements, as in [12,13]. This work is focused on low-frequency circulating currents.
To reduce circulating currents, in some works, such as [14,15], the n modules connected in parallel have been considered as a single converter that has as many branches as the branches of each module multiplied by n. This control is very complex for high values of n and does not allow for a modular design. Other works, such as [13], proposed a finite time control method. In [16], a double modulation compensation technique was used; in [8], a harmonic elimination PWM modulation (HEPWM) that needs a high switching frequency, leading to high switching losses, was proposed. Other studies proposed the control of the zero-sequence current of the inverters [17,18]. In [17], a common-mode voltage injection PWM method was described, and [18] proposed a multicarrier PWM for parallel three-phase active front-end converters. In both cases, a control action was added to the three control signals of the PWM modulator to reduce the circulating currents. In those works, the effect that the addition of the proposed control signal had on the stability of the system had not been analyzed. Reference [19] proposed the control of the homopolar component of the phase currents to avoid circulating currents using three-dimensional space vector modulation (3D-SVM).
This paper proposes a control structure based on the fact that the circulating currents among inverters agree with the zero-sequence component of the three-phase currents, so that, by regulating the zero-sequence component to a null setpoint value, circulating currents are eliminated. The complexity of the proposed solution is very low, allowing for a modular design that can be easily extended to a high value of n. Regarding [19], where only differences among inverter inductances were taken into account, the present work considers various causes that produce circulating currents: mismatches in the inductance values of the grid filters among inverters and in each phase of the same inverter, unbalance of the delivered power, and the use of different modulation techniques. In addition, the harmonic contents of the circulating currents have been considered to adjust the resonant controllers in the inverter’s current loops, which improves the control action at the frequencies of interest. The control structure proposed in the current work has been validated in an experimental laboratory set-up. The stability of the whole system has been analyzed and supported by both analytical and experimental results.

2. Modeling of n Inverters Connected in Parallel in the Synchronous Reference Frame

Figure 1 shows a general scheme of n inverters connected in parallel. The grid filter of each inverter is an LCL filter with a damping resistance, Rd, connected in series with the filter capacitance, Cf. The inductances of the filter can be implemented by means of three-phase inductors, being Lf and Lfg, the inverter-side inductance and grid-side self-inductance, respectively, whereas rf and rfg are the corresponding parasitic series resistances. Mf and Mfg stand for the inverter-side and grid-side mutual inductance values, respectively. The grid impedance at the point of common coupling (PCC) is characterized by its self-inductance, Lg, parasitic series resistance, rg, and mutual inductance, Mg.
The goal of this section is to obtain a MIMO model of the system in a dqo synchronous reference frame (SRF), having the grid voltage aligned with the d-axis. Then, the active and reactive powers can be regulated with the d and q components of the AC currents, respectively. It is worth pointing out that the zero-sequence component (channel o) of the current is normally omitted in three-phase inverters because this component is naturally zero in the case of three-wire inverters without a neutral conductor. In the proposed control structure, the zero-sequence component will also be regulated.

2.1. Averaged Model Equations

Figure 2 shows the averaged equivalent circuit of Figure 1 represented in the stationary three-phase reference frame. Equations (2)–(6) are obtained from this scheme. In Equations (2)–(6), the value of i (i = 1, …, n) identifies each module, with n being the total number of modules. The load factor of each module, ci, is defined by Equation (1), with I2i_RMS_OP being the RMS phase current at the operation point provided by inverter #i and I2i_RMS_NOM being its nominal RMS phase current. Fm stands for the gain of the SVM modulation and depends on its implementation inside the digital controller in which the control loops are programmed.
c i = I 2 i _ R M S O P I 2 i _ R M S N O M ,           i = 1 , 2 , n
i ¯ g φ = i ¯ 2 φ 1 + i ¯ 2 φ 2 + + i ¯ 2 φ n = j = 1 n c j c i i ¯ 2 φ i ,           i = 1 , 2 , n ,           φ = a , b , c
L f M f M f M f L f M f M f M f L f · d d t i ¯ 1 a i i ¯ 1 b i i ¯ 1 c i = F m · v ¯ d c d ¯ a i d ¯ b i d ¯ c i r f i ¯ 1 a i i ¯ 1 b i i ¯ 1 c i v ¯ c a i v ¯ c b i v ¯ c c i R d i ¯ 1 a i i ¯ 1 b i i ¯ 1 c i i ¯ 2 a i i ¯ 2 b i i ¯ 2 c i v ¯ N i v ¯ N i ´ v ¯ N i ´ ,           i = 1 , 2 , n
L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g · d d t i ¯ 2 a i i ¯ 2 b i i ¯ 2 c i = v ¯ N i ´ v ¯ N i ´ v ¯ N i ´ ( r f g + j = 1 n c j c i r g ) i ¯ 2 a i i ¯ 2 b i i ¯ 2 c i + v ¯ c a i v ¯ c b i v ¯ c c i + R d i ¯ 1 a i i ¯ 1 b i i ¯ 1 c i i ¯ 2 a i i ¯ 2 b i i ¯ 2 c i v ¯ g a v ¯ g b v ¯ g c v ¯ N i v ¯ N i v ¯ N i ,           i = 1 , 2 , n
d d t v ¯ c a i v ¯ c b i v ¯ c c i = 1 C f i ¯ 1 a i i ¯ 1 b i i ¯ 1 c i i ¯ 2 a i i ¯ 2 b i i ¯ 2 c i ,           i = 1 , 2 , n
d d t v ¯ d c = 1 n C o i ¯ d c F m d ¯ a 1 d ¯ b 1 d ¯ c 1 T i ¯ 1 a 1 i ¯ 1 b 1 i ¯ 1 c 1 F m d ¯ a 2 d ¯ b 2 d ¯ c 2 T i ¯ 1 a 2 i ¯ 1 b 2 i ¯ 1 c 2 F m d ¯ a n d ¯ b n d ¯ c n T i ¯ 1 a n i ¯ 1 b n i ¯ 1 c n
It is worth pointing out that the circulating currents comply with Equation (7), so it is only necessary to regulate the n − 1 zero-sequence currents (io1 to ion) to zero, being null to the zero-sequence duty cycle of the n-th converter. Taking into account that the filter capacitors are not connected to the midpoint of the dc-link, the circulating currents neither flow through them nor to the grid, as it was shown in [12]. In view of this, and following the procedure described in Appendix A, the equations in the SRF have been obtained, Equations (8)–(12).
i ¯ o n = i ¯ o 1 + i ¯ o 2 + i ¯ o n 1
d d t   i ¯ 1 d i i ¯ 1 q i = 0 ω ω 0 i ¯ 1 d i i ¯ 1 q i + F m   v ¯ d c L f M f d ¯ d i d ¯ q i 1 L f M f v ¯ c d i v ¯ c q i v ¯ c o i R d + r f L f M f   i ¯ 1 d i i ¯ 1 q i + R d L f M f i ¯ 2 d i i ¯ 2 q i ,     i = 1 , 2 , n
d d t i ¯ 2 d i i ¯ 2 q i = 0 ω ω 0 i ¯ 2 d i i ¯ 2 q i + 1 L f g M f g + j = 1 n c j c i ( L g M g ) v ¯ c d i v ¯ c q i + R d L f g M f g + j = 1 n c j c i ( L g M g ) i ¯ 1 d i i ¯ 1 q i R d + r f g + j = 1 n c j c i r c L f g M f g + j = 1 n c j c i ( L g M g ) i ¯ 2 d i i ¯ 2 q i 1 L f g M f g + j = 1 n c j c i ( L g M g ) v ¯ g d v ¯ g q ,           i = 1 , 2 , n
d d t i ¯ o 1 i ¯ o 2 i ¯ o n 1 = r f + r f g L f + L f g + 2 M f + 2 M f g 0 0 0 r f + r f g L f + L f g + 2 M f + 2 M f g 0 0 0 r f + r f g L f + L f g + 2 M f + 2 M f g i ¯ o 1 i ¯ o 2 i ¯ o n 1 + ( n 1 ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) ( n 1 ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g ) ( n 1 ) F m v ¯ d c n ( L f + L f g + 2 M f + 2 M f g )
d d t v ¯ c d i v ¯ c q i = 0 ω ω 0 v ¯ c d i v ¯ c q i + 1 C f i ¯ 1 d i i ¯ 1 q i i ¯ 2 d i i ¯ 2 q i ,           i = 1 , 2 , n
d d t v ¯ d c = 1 n C o i ¯ p v F m d ¯ d 1 d ¯ q 1 d ¯ o 1 T i ¯ 1 d 1 i ¯ 1 q 1 i ¯ o 1 F m d ¯ d 2 d ¯ q 2 d ¯ o 2 T i ¯ 1 d 2 i ¯ 1 q 2 i ¯ o 2 F m d ¯ d n 1 d ¯ q n 1 d ¯ o n 1 T i ¯ 1 d ( n 1 ) i ¯ 1 q ( n 1 ) i ¯ o ( n 1 ) F m d ¯ d n d ¯ q n T i ¯ 1 d n i ¯ 1 q n

2.2. Derivation of the Small Signal MIMO Model

The small signal model of the parallel inverters has been obtained in this section. It has been obtained by applying Equation (13) to SRF Equations (8)–(12). In Equation (13), x ¯ stands for each of the variables of the averaged model, expressed as the sum of the value at the operation point, X, and the small signal value, x ^ .
x ¯ = X + x ^
The equations of the small signal model in the SRF frame are expressed in Equations (14)–(21). Because of the model complexity, some auxiliary matrices (A9)–(A25) have been defined to simplify A, B, C, and D. The input vector U defines the control variables (the duty cycle of each branch of the inverter) and the disturbances (the grid voltage). The constant Kvo is expressed by Equation (A19), considering that at the DC side, there is a constant power load (CPL) or a photovoltaic field that can be modelled as a current source.
d d t X = A · X + B · U
Y = C · X + D · U
X = Y = [ i ^ 1 d 1 i ^ 1 q 1 i ^ 2 d 1 i ^ 2 q 1 i ^ o 1 v ^ c d 1 v ^ c q 1 i ^ 1 d n 1 i ^ 1 q n 1 i ^ 2 d n 1 i ^ 2 q n 1 i ^ o n 1 v ^ c d n 1 v ^ c q n 1 i ^ 1 d n i ^ 1 q n i ^ 2 d n i ^ 2 q n v ^ c d n v ^ c q n v ^ d c T
U = d ^ d 1 d ^ q 1 d ^ o 1 d ^ d n 1 d ^ q n 1 d ^ o n 1 d ^ d n d ^ q n v ^ g d v ^ g q v ^ g o T
A = A i 1 ^ 2 d q _ i 12 d q 0 A i ^ 12 d q _ v c d q 0 0 0 0 0 A i ^ 12 d q _ v d c 0 A i ^ o _ i o 0 0 0 0 0 0 A i ^ o _ v d c A v ^ c d q _ i 12 d q 0 A v ^ c d q _ v c d q 0 0 0 0 0 0 0 0 0 A i 1 ^ 2 d q _ i 12 d q 0 A i ^ 12 d q _ v c d q 0 0 A i ^ 12 d q _ v d c 0 0 0 0 A i ^ o _ i o 0 0 0 A i ^ o _ v d c 0 0 0 A v ^ c d q _ i 12 d q 0 A v ^ c d q _ v c d q 0 0 0 0 0 0 0 0 0 A i 1 ^ 2 d q _ i 12 d q A i ^ 12 d q _ v c d q A i ^ 12 d q _ v d c 0 0 0 0 0 0 A v ^ c d q _ i 12 d q A v ^ c d q _ v c d q 0 A v ^ d c _ i 12 d q A v ^ d c _ i o 0 A v ^ d c _ i 12 d q A v ^ d c _ i o 0 A v ^ d c _ i 12 d q 0 A v ^ d c _ v d c
B = B i ^ 12 d q _ d d q 0 0 0 0 B i ^ 12 d q _ v g d q o 0 B i ^ o _ d o i 0 B i ^ o _ d o 0 0 0 0 0 0 0 0 0 0 B i ^ 12 d q _ d d q 0 0 B i ^ 12 d q _ v g d q o 0 B i ^ o _ d o 0 B i ^ o _ d o i 0 0 0 0 0 0 0 0 0 0 0 0 B i ^ 12 d q _ d d q B i ^ 12 d q _ v g d q o 0 0 0 0 0 0 B v ^ d c _ d d q B v ^ d c _ d o B v ^ d c _ d d q B v ^ d c _ d o B v ^ d c _ d d q 0
C = I
D = 0
The expressions of the variables at the operation point are summarized in Table 1. In order to simplify their calculation, the inductor resistances have been neglected, as well as the filter capacitors, so that the currents at the grid side (I1dqi) and at the converter side (I2dqi) agree (IdqiI1dqiI2dqi). In addition, the expressions of the operation point have been obtained by considering the control objectives. Since zero reactive power and zero circulating currents are aimed at, both the reactive and zero-sequence components of the currents are chosen to be null. The value of the load factor for each inverter, ci, can vary between 0 and 1.

3. Control Architecture

Figure 3 shows the scheme of the proposed control scheme. It is composed of n current control loops for d- and q-channels (in both cases, one per inverter) and n − 1 current control loops for regulating the o-channel.
By aligning the d-axis of the reference frame with the direct-sequence fundamental space vector of the grid voltage, the active power may be controlled by regulating the d component of the current, whereas the reactive power can be controlled by acting on the q-axis component. In this case, the reference for the q-axis current is zero, i q * = 0 , in order to obtain a unity power factor. The references in the d- and o- axis are defined by id* and io*, respectively.
In three-phase inverters, the zero-sequence component of the phase currents represents the circulating currents among the inverters. The zero-sequence component is defined as one-third of the sum of the phase currents provided by an inverter, io = (ia + ib + ic)/3. If only one inverter is connected, io1 = 0, as the sum of the three-phase currents is zero by application of Kirchhoff’s first law (ia1 + ib1 + ic1 = 0). If two inverters are connected in parallel, n = 2, io1 = −io2, and the circulating currents can be controlled with only one zero-sequence current control. When there are n inverters connected in parallel, they need n − 1 zero-sequence current control regulators. As has been previously pointed out, there is no need to control the zero-sequence component of a three-wire, three-phase inverter when it is operating as a single inverter. However, that is not true in the case of the parallel connection of inverters due to the existence of circulating currents, which are directly related to the zero-sequence component. Therefore, in the proposed solution, the zero-sequence components are controlled with a reference signal io* = 0 to suppress the circulating currents among inverters.
Although the zero-sequence components of the currents can be controlled by means of a sinusoidal PWM modulation, the use of SVM modulation has been preferred. The reason is that SVM modulations are digital, so they allow for a much more flexible trip sequence of the transistors than analog PWM modulations. Generally, in three-phase, three-wire systems, 2D-SVM, or simply SVM, is used. However, to control the zero-sequence component of the currents, 3D-SVM is needed, since the three components αβo of the phase voltage can be controlled in the whole operation range of the inverter.
When inverters are represented in the synchronous reference frame, coupling terms between the d- and q-axis of the phase currents appear. Equations (24) and (25) express the proposed decoupling terms between both axes, which have been calculated by omitting the effect of the filter capacitors. Under this assumption, the currents at the inverter side (I1dqoi) and at the grid side (I2dqoi) agree (IdqoiI1dqoiI2dqoi). Assuming no disturbances, both at the DC-link voltage and at the grid voltage, v ^ d c s = v ^ g d s = 0 , and Equation (22) is obtained. Since εidi is the current error signal in the d-axis and Gi is the current regulator transfer function, the control action in channel d- can be expressed according to Equation (23). The decoupling term Kqd is obtained from Equations (22) and (23). The expression of Kdq is calculated similarly. The decoupling terms are valid from DC to medium frequencies, where the effect of the filter capacitance in the coupling can be neglected.
Note that the decoupling terms depend on the number of inverters on stream, n. If the grid inductances are negligible, the parameters Lg and Mg could be considered null, so that the decoupling terms would be constantly independent of n. However, for the cases where the value of Lg and Mg could not be considered zero, it is easy to automatically readjust the value of the decoupling terms by stocking the needed value as a function of n.
d ^ d i s · F m · V d c + ω L f M f + L f g M f g + j = 1 n c j c i L g M g · i ^ q i s = L f M f + L f g M f g + j = 1 n c j c i L g M g s · i ^ d i s ,           i = 1 , 2 , n
d ^ d i s = G i s ε ^ i d i s + K q d i · i ^ q i s ,           i = 1 , 2 , n
k q d i = ω L f M f + L f g M f g + j = 1 n c j c i L g M g F m · V d c ,           i = 1 , 2 , n
k d q i = ω L f M f + L f g M f g + j = 1 n c j c i L g M g F m · V d c ,           i = 1 , 2 , n
Equations (26)–(28) express the current control open loop gains in channels d, q, and o, respectively. In those equations, Rs is the gain of the current sensor, Gi(s) is the current regulator transfer function, D(s) represents the transfer function of a delay equal to a switching period, which has been calculated from a second-order Padé approximation, and FPB(s) is the transfer function of a low-pass antialiasing filter implemented in the current sensors. The transfer functions of the current as a function of the duty cycle has been obtained from the state space model of the system using MATLABTM (R2018b). The current sensors have been located at the inverter side of the LCL filter. Placing the sensors at the inverter side, instead of sensing the currents at the grid side, provides higher robustness to the control system [20].
T i d i s = R s · G i s · D s · F P B s · i ^ i 1 d i s d ^ d i s ,           i = 1 , 2 , n
T i q i s = R s · G i s · D s · F P B s · i ^ i 1 q i s d ^ q i s ,           i = 1 , 2 , n
T i o i s = R s · G i s · D s · F P B s · i ^ i 1 o i s d ^ o i s ,           i = 1 , 2 , n

4. Three-Dimension Space Vector Modulator

If the fundamental components of the phase currents are balanced in a three-phase, two-level inverter, the fundamental component of the voltage between the midpoint of the dc-link and the neutral point of the grid is zero. Thus, in that case, it can be considered a virtual connection at the fundamental frequency of the line voltage between the neutral point of the grid and the midpoint of the dc-link [21]. This is an important difference, with respect to three-phase inverters with a four-wire output (three phases + neutral), where the neutral is physically connected. Three-dimensional SVM modulators have been developed for four-wire connected inverters, and, to the best of the authors’ knowledge, the first time 3D-SVM was applied to three-phase inverters without neutral was in this work.
Table 2 shows the phase-to-neutral voltage in abc coordinates using the DC voltage vdc as a reference of the eight possible switching vectors with the considerations described above. It also shows the same voltage in the αβ0 coordinates (Clarke). Figure 4a shows the twelve tetrahedra, which are defined by two consecutive nonzero switching vectors (in Clarke coordinates) and neutral vectors v0 and v7, alternately. The composition of the above tetrahedra defines the whole operation range of the inverter represented in Figure 4b.
Vector modulation consists of constructing the reference voltage vector v* at the inverter output by alternately applying its adjacent switching vectors. First, the tetrahedron containing the reference vector must be obtained. Once the tetrahedron containing the reference vector is defined, the switching vectors adjacent to v* are known. From the projections of the reference vector on the switching vectors that delimit the tetrahedron, the time that each switching vector must be applied was calculated in [21,22,23]. In this work, a symmetric distribution of the switching vectors has been implemented. Table 3 represents the sequence of the switching vectors applied in each prism.

5. Experimental Set-Up

A 10 kW prototype composed of two three-phase 5 kW inverters has been built and experimentally tested. Figure 5 shows the single-line scheme of the converters, whereas Table 4 summarizes the parameters of the system. As can be seen from Figure 5, the grid filter was an LCL grid filter composed of Cf capacitors with Rd damping resistors, three single-phase inductances at the inverter side with a value Lf, and a three-phase inductor Lg, Mg, common to the two inverters connected in parallel. In this case, the inductance Lfg of the general model was not used, so it was considered null.
The experimental results have been obtained by means of the set-up depicted in Figure 6. Both parallel operating inverters have been connected to the same AC and DC buses and have been controlled using the dual-core Texas Instruments TMS320F28379D microcontroller. The AC bus has been supplied by a bidirectional power supply, Cinergia GL and EL-50, emulating an AC grid of 230 V RMS phase–phase and 50 Hz. Anyway, a GSS Regatron programmable power supply has set the DC bus to 500 V. Table 5 summarizes the test equipment of the experimental set-up.

6. Design of the Control Loops

PI controllers are usually used in power converters to implement the regulators of the current control loops in a synchronous reference frame. However, PI controllers have difficulties when the regulated variables have frequency contents close to the crossover frequency of the loops, which is the case of the low-frequency harmonics of the circulating currents. To address this issue, second-order generalized integrators (SOGI) [24] have been added to the o-channel controllers to increase their gain at the frequency of the circulating currents’ low-frequency harmonics, resulting in PI + SOGI o-channel controllers. As it will be seen in the simulation and experimental results section, when the powers provided by the parallel inverters are balanced, the circulating currents present third-order harmonics (150 Hz) and their odd multiples, with the ninth one being the most significant harmonic (450 Hz). Nevertheless, if the power provided by one of the parallel inverters becomes unbalanced, a large fundamental harmonic (50 Hz) appears at the circulating current. Therefore, the SOGIs of the o-channel current controllers are tuned to 50 Hz, 150 Hz, and 450 Hz.
The expressions of the current controllers should be chosen by taking the following issues into account. First, the crossover frequency of the current control loop gains must be higher than the frequency of the harmonics of the variables to be controlled. In addition, the loop gain must be much lower (around 1 decade) than the switching frequency to avoid large signal instabilities and noise problems. In addition, appropriate stability margins must be achieved, typically a phase margin PM > 45° and a gain margin GM > 6 dB. Equations (29) and (34) show the expressions of the chosen current regulators on the d- and q-axis and on the o-axis, respectively. The transfer functions of the d and q current channels were similar, so that their controllers could agree. The current controllers in both channels were proportional integral regulators. As indicated above, the current regulator in the zero-sequence channel was a PI + SOGIs one, being the generalized integrators tuned at 50 Hz, 150 Hz, and 450 Hz. The bandwidth of the 50 Hz SOGI was 10 Hz, and its gain was 4. The gain of the 150 Hz SOGI was 4, whereas that of the 450 Hz SOGI was 0.5. The bandwidths of the 150 Hz and 450 Hz SOGIs were reduced, with regard to the bandwidth of the 50 Hz SOGI. The reason for that was to avoid the risk of instability, due to a phase loss close to the crossover frequency of the o-channel current loops.
G i s = P I i s = 0.1 + 10 s
P I o s = 0.2 + 10 s
R 50 s = 4 · 10 · s s 2 + 10 s + 2 τ ω 2
R 150 s = 4 · 10 3 s s 2 + 10 3 s + 2 τ 3 ω 2
R 450 s = 0.5 · 10 9 s s 2 + 10 9 s + 2 τ 9 ω 2
G i o s = P I o s + R 50 s + R 150 s + R 450 s
Figure 7a–c depict the Bode plot measurements and the theoretical curves (gain in dB, phase in deg) of the current loop gains: Tid () Equation (26), Tiq () Equation (27), and Tio () Equation (28). The theoretical plots have been obtained by Matlab R2018b. The experimental measurements of Figure 7 have been carried out by means of a frequency response analyzer FRA5097 from NF Corporation. Figure 7d depicts the Bode plot measurements of the current loop gain in the o-channel.
As it can be seen from those figures, the current loop frequency response predicted by the developed model was reasonably close to the experimental measurements, so the developed model was validated as a design and analysis tool. Note that stability was guaranteed in all cases. Specifically, the crossover frequency of the current loops ranged from 680 Hz to 800 Hz, having, in the worst case, the phase margin PM = 47° and the gain margin GM = 7.2 dB.
The stability of the system has been also analysed by considering variations of Lf in the range between 5–7 mH. Figure 8a,b represent the pole-zero map in the d- and q-axis, respectively. It can be seen that, for higher values of Lf, the system is slower and more stable. Figure 8c shows that the Lf variations do not affect the o-axis stability. The position of the closed loop pole trajectories has been represented in these figures. It can be seen that, for higher values of Kp, the current control loops become faster, compromising the system stability.

7. Simulation and Experimental Results

In this section, the simulation and experimental results that validate the proposed control scheme for the suppression of circulating currents are shown. The prototype has been tested by considering that low-frequency circulating currents are due to various reasons: an inductance value mismatch in the phases of the same inverter and among different inverters, power unbalance among inverters, and the use of different PWM modulators in each inverter. In practice, the last scenario happens if the inverter that does not need to regulate the homopolar component (#1 in Figure 3) implements 2D-SVM, instead of 3D-SVM.
Simulation results have been obtained by considering two and three inverters with the characteristics of Table 4 connected in parallel. The experimental results with two inverters connected in parallel have been presented. Only one of the previously described reasons for unbalance has been studied in each experiment.

7.1. Simulation Results

7.1.1. Two Inverters Connected in Parallel

The simulation results, when two 5 kW inverters are connected in parallel, have been obtained by means of PSIMTM 10.0 and are depicted in Figure 9. In this figure, the circulating currents are represented before and after the activation of the zero-sequence control at t = 0.25 s. The upper part of each of the graphs corresponds to the phase currents of inverter #1, the lower part depicts the phase currents of inverter #2, and in the central part, the circulating currents are shown.
Figure 9a shows the simulation results when the value of the inverter side inductances of inverter #2, Lf, is modified to 7 mH. Despite the fact that the circulating currents are inherently small, without the need to control the zero-sequence component, it was found that, with the proposed control of the zero-sequence channel, the 150 Hz component decreased significantly. As described in Section 6, when the currents provided by the inverters are balanced, the circulating current is composed of a three-order harmonic and its odd multiples if there are mismatches in the value of Lf among inverters and in other cases, as will be shown by the experimental results shown in the following.
Figure 9b shows the simulation results when the load factor of the two inverters becomes unbalanced, so that inverter #1 handles 25% of its nominal power, whereas inverter #2 handles 50% of its nominal power. In this case, the circulating currents are of a higher value than in the case in which the inductances of both inverters are unbalanced, so it is easier to see that the zero-sequence control allows the circulating currents to be significantly reduced.
Figure 9c shows the results when the inductance of phase A of the second inverter is 7 mH, instead of 5 mH. In this case, the circulating current had a large fundamental harmonic, and the phase currents became unbalanced. After activating the zero-sequence control at t = 0.25 s, the circulating current was significantly reduced, leading to a good balancing of the currents, provided by both inverters.
Finally, in Figure 9d, the cases of using a 2D-SVM modulator in inverter #1 and a 3D-SVM modulator in inverter #2 are depicted. In this case, the amplitude of the 150 Hz fundamental component of the circulating current was approximately 4 A, being drastically reduced after activating the zero-sequence component control. It should also be noted that the currents in the phases of each inverter are no longer strongly distorted and are transformed into 50 Hz sinusoidal signals.

7.1.2. Three Inverters Connected in Parallel

In this section, the proposed technique has been applied to three inverters of 5 kW connected in parallel. Since n − 1 zero-sequence current control loops are needed to control the circulating currents, the zero-sequence component of the phase currents in inverters #2 and #3 have been controlled.
Figure 10a shows the simulation results with a mismatch in the nominal value of the inductances among inverters. Specifically, Lf in inverter #1 had a value of 5 mH, Lf in inverter #2 had a value of 7 mH, and Lf in inverter #3 was 6 mH. Figure 10b represents the phase currents and the circulating currents when inverter #1, inverter #2, and inverter #3 handled 25%, 30%, and 50% of their nominal power, respectively. Figure 10c depicts the case in which there was a mismatch among the inductances of each phase of the inverter. The phase B inductance of inverter #1 was 6 mH, and the phase A inductance of inverter #2 was 7 mH, instead of 5 mH. Finally, in Figure 10d, the simulation results are shown with a 2D-SVM modulator in inverter #1 and 3D-SVM modulators in inverters #2 and #3.
Note that these results are very similar to the previous case, and that the circulating currents have been drastically reduced with the proposed control structure.

7.2. Experimental Results

In the experimental tests, the factors analyzed in Section 7.1 causing circulating currents have been studied. It should be noted, however, that the experimental results were obtained from real power converters, so it was impossible to study each of the factors in a completely independent way. Indeed, the inductance value used in the grid filter was theoretically 5 mH, but inductances from the same manufacturer have tolerances of up to 20%, even higher in some cases. In this work, inductances with nominal values of 5 mH, 7 mH and tolerances of 10% have been used. To determine the real value of each component, the inductance of each inductor has been measured, resulting in the values collected in Table 6.
At the top of the graphs in Figure 11, the currents in phase A (orange), B (green), and C (pink) of inverter #1 are shown. The phase currents of inverter #2 are shown at the bottom. It should be noted that, in the central part of each graph, the sum of the currents in the phases of inverter #1 (ia1 + ib1 + ic1) are shown in blue, and the sum of the currents of the phases of inverter #2 are represented in pink. These currents are proportional to the circulating currents between inverters, whose expression is (ia + ib + ic)/3. A dotted black vertical line marks the activation instant of the proposed control for suppression of the circulating currents.
Figure 11a shows the experimental performance of the proposed control technique, with a mismatch between the inductances of the inverters. In inverter #2, the inductor with a nominal inductance of 5 mH has been replaced by an inductance of 7 mH. The actual values of the inductances used in each phase are summarized in Table 7. The results show that a 50 Hz component appeared in the circulating currents, since the inductance differences among phases did not allow for balancing the phases. A small 150 Hz component also appeared when the grid filter inductances between both inverters were different.
Figure 11b shows the experimental results of the control of circulating currents for a load factor of 25% in inverter #1 and 50% in inverter #2. In this case, the 50 Hz and 150 Hz harmonics also appeared, with the 150 Hz harmonic being greater than in the previous case. The proposed control technique considerably reduced the circulating current.
Figure 11c shows the results in the case in which the inductance connected to phase A of the second inverter was replaced by the inductance Lf_a2 of 7.16 mH, represented in Table 7. It was observed that, before activating the zero-sequence control, and with both inverters working at their rated power, the sum of the phase currents resulted in a 50 Hz component of about 3.5 A, so the circulating current had a value of about 1.2 A. After activation of the zero-sequence current control, the circulating current was reduced considerably.
Finally, in Figure 11d, inverter #1 works with 2D-SVM and inverter #2 works with 3D-SVM. It was observed that, before activating the control of the zero-sequence component, the currents in the phases were strongly distorted, due to the high value of the 150 Hz component, 4.5 A, that appeared in the circulating current. After the activation of the zero-sequence control, the currents did not exhibit a noticeable distortion.

8. Conclusions

This work presents a method for the suppression of low-frequency circulating currents among n three-phase inverters connected in parallel. To this end, a control technique, in which n − 1 regulation loops with three-dimensional space vector modulators (3D-SVM), was used to regulate the zero-sequence component of the currents in each inverter. Thus, the circulating currents were controlled as they coincided with the zero-sequence components (io = (ia + ib + ic)/3). As a control objective, a zero setpoint value for the zero-sequence has been imposed, so that the circulating currents cancel out.
Several factors that produce circulating currents have been analyzed: mismatches in the inductance values of the grid filters among inverters and in each phase of the same inverter, unbalance of the delivered power, and finally, the use of different modulation techniques. In addition, the harmonic contents of the circulating currents have been considered to adjust resonant controllers in the zero-sequence current loops. This improves the control action at the frequency of the harmonics that form part of the circulating currents.
The simulation and experimental results have been obtained to validate the correct operation of the proposed control. The most relevant cases that can produce circulating currents among inverters connected in parallel have been evaluated in the tests. From the results obtained, we can conclude that the value of the circulating currents in the event of a mismatch among the nominal value of the inductances of each inverter is very low. The same has been observed when the load factor of each inverter was different. However, when there were differences among the phase inductances of the same inverter, the circulating current at nominal power had a high 50 Hz harmonic of 1.2A. With the proposed control technique, this harmonic was reduced to 8 mA (99% attenuation rate). Finally, when inverter #1 was controlled with 2D-SVM and inverter #2 with 3D-SVM, the circulating current had a high 150 Hz harmonic of 4.5 A, which reduced to 100 mA with the zero-sequence current control loop (98% attenuation rate).

Author Contributions

M.L., E.F. and G.G. developed the main idea, the model, and the control of the system; M.L., I.P. and R.G.-M. developed the experimental results; M.L., E.F. and G.G. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish “Ministerio de Asuntos Económicos y Transformación Digital” and the European Regional Development Fund (ERDF), under grants RTI2018-100732-B-C21 and PID2021-122835OB-C22.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In order to obtain the equations in the SRF frame, Equations (8)–(12), (A29 were applied to Equations (3)–(6), where T was the SRF transformation matrix expressed by Equation (A3). In these equations, ω stands for the grid angular frequency, and x represents any averaged variable (current, voltage, …).
T = 2 3 cos ω t cos ω t 2 π 3 cos ω t + 2 π 3 sin ω t sin ω t 2 π 3 sin ω t + 2 π 3 1 2 1 2 1 2
x ¯ d x ¯ q x ¯ o = T · x ¯ a x ¯ b x ¯ c x ¯ a x ¯ b x ¯ c = T 1 x ¯ d x ¯ q x ¯ o
T d T 1 d t = 0 ω 0 ω 0 0 0 0 0
T L f M f M f M f L f M f M f M f L f d T 1 d t i ¯ 1 d i i ¯ 1 q i i ¯ 1 o i = 0 ω L f M f 0 ω L f M f 0 0 0 0 0 i ¯ 1 d i i ¯ 1 q i i ¯ 1 o i   ,           i = 1 , 2 , n
T L f M f M f M f L f M f M f M f L f T 1 i ¯ 1 d i i ¯ 1 q i i ¯ 1 o i = L f M f 0 0 0 L f M f 0 0 0 L f + 2 M f i ¯ 1 d i i ¯ 1 q i i ¯ 1 o i ,           i = 1 , 2 , n
T L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g d T 1 d t i ¯ 2 d i i ¯ 2 q i i ¯ 2 o i = 0 ω ( L f g M f g + j = 1 n c j c i ( L g M g ) ) 0 ω ( L f g M f g + j = 1 n c j c i ( L g M g ) ) 0 0 0 0 0 i ¯ 2 d i i ¯ 2 q i i ¯ 2 o i , i = 1 , 2 , n
T L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g M f g + j = 1 n c j c i M g L f g + j = 1 n c j c i L g T 1 i ¯ 2 d i i ¯ 2 q i i ¯ 2 o i = L f g M f g + j = 1 n c j c i ( L g M g ) 0 0 0 L f g M f g + j = 1 n c j c i ( L g M g ) 0 0 0 L f g + 2 M f g + j = 1 n c j c i ( L g + 2 M g ) i ¯ 2 d i i ¯ 2 q i i ¯ 2 o i , i = 1 , 2 , n
T 1 x ¯ d x ¯ q x ¯ o T T 1 = x ¯ d x ¯ q x ¯ o T T 1 T T 1 = x ¯ d x ¯ q x ¯ o T

Appendix B

A i 1 ^ 2 d q _ i 12 d q = ( R d + r f ) L f M f ω R d L f M f 0 ω ( R d + r f ) L f M f 0 R d L f M f R d L f g M f g + j = 1 n c j c i ( L g M g ) 0 ( R d + r f g + j = 1 n c j c i r g ) L f g M f g + j = 1 n c j c i ( L g M g ) ω 0 R d L f g M f g + j = 1 n c j c i ( L g M g ) ω ( R d + r f g + 1 n c i c i ) L f g M f g + j = 1 n c j c i ( L g M g ) , i = 1 , 2 , n
A i ^ o _ i o = r f + r f g L f + L f g + 2 M f + 2 M f g
A i ^ 12 d q _ v c d q = 1 L f M f 0 0 1 L f M f 1 L f g M f g + j = 1 n c j c i L g M g 0 0 1 L f g M f g + j = 1 n c j c i L g M g ,           i = 1 , 2 , n
A i ^ 12 d q _ v d c = F m · D d i L f M f F m · D q i L f M f 0 0 ,           i = 1 , 2 , n
A i ^ o _ v d c = F m · D o i n ( L f + L f g + 2 M f + 2 M f g ) ,           i = 1 , 2 , n
A v ^ c d q _ i 12 d q = 1 C f 0 1 C f 0 0 1 C f 0 1 C f
A v ^ c d q _ v c d q = 0 ω ω 0
A v ^ d c _ i 12 d q = F m · D d i n C o F m · D q i n C o 0 0 ,           i = 1 , 2 , n
A v ^ d c _ i o = F m · D o i n C o ,           i = 1 , 2 , n
A v ^ d c _ v d c = K v o n C o
K v o = I d c V d c
B i ^ 12 d q _ d d q = F m · V d c L f M f 0 0 F m · V d c L f M f 0 0 0 0
B i ^ o _ d o i = n 1 F m   v ¯ d c n ( L f + L f g + 2 M f + 2 M f g )
B i ^ o _ d o = F m   v ¯ d c n ( L f + L f g + 2 M f + 2 M f g )
B i ^ 12 d q _ v g d q o = 0 0 0 0 0 0 1 L f g M f g + j = 1 n c j c i L g M g 0 0 0 1 L f g M f g + j = 1 n c j c i L g M g 0 ,   i = 1 , 2 , n
B v ^ d c _ d d q = F m · I d i n C o F m · I q i n C o ,           i = 1 , 2 , n
B v ^ d c _ d o = F m · I o i n C o ,           i = 1 , 2 , n

References

  1. Xueguang, Z.; Li, W.; Xiao, Y.; Wang, G.; Xu, D. Analysis and Suppression of Circulating Current Caused by Carrier Phase Difference in Parallel Voltage Source Inverters With SVPWM. IEEE Trans. Power Electron. 2018, 33, 11007–11020. [Google Scholar] [CrossRef]
  2. Shen, Z.; Jiang, D.; Chen, J.; Qu, R. Circulating Current Reduction for Paralleled Inverters with Modified Zero-CM PWM Algorithm. IEEE Trans. Ind. Appl. 2018, 54, 3518–3528. [Google Scholar] [CrossRef]
  3. Wang, X.; Qing, H.; Huang, P.; Zhang, C. Modeling and Stability Analysis of Parallel Inverters in Island Microgrid. Electronics 2020, 9, 463. [Google Scholar] [CrossRef] [Green Version]
  4. Xiao, H.; Xie, S.; Chen, Y.; Huang, R. An Optimized Transformerless Photovoltaic Grid-Connected Inverter. IEEE Trans. Ind. Electron. 2011, 58, 1887–1895. [Google Scholar] [CrossRef]
  5. Zhang, M.; Song, B.; Wang, J. Circulating Current Control Strategy Based on Equivalent Feeder for Parallel Inverters in Islanded Microgrid. IEEE Trans. Power Syst. 2019, 34, 595–605. [Google Scholar] [CrossRef]
  6. Zou, Z.-X.; Hahn, F.; Buticchi, G.; Günter, S.; Liserre, M. Interleaved Operation of Two Neutral-Point-Clamped Inverters With Reduced Circulating Current. IEEE Trans. Power Electron. 2018, 33, 10122–10134. [Google Scholar] [CrossRef]
  7. Narimani, M.; Moschopoulos, G. Improved Method for Paralleling Reduced Switch VSI Modules: Harmonic Content and Circulating Current. IEEE Trans. Power Electron. 2014, 29, 3308–3317. [Google Scholar] [CrossRef]
  8. Chen, T.-P. Zero-Sequence Circulating Current Reduction Method for Parallel HEPWM Inverters Between AC Bus and DC Bus. IEEE Trans. Ind. Electron. 2012, 59, 290–300. [Google Scholar] [CrossRef]
  9. Mazumder, S.K. A Novel Discrete Control Strategy for Independent Stabilization of Parallel Three-Phase Boost Converters by Combining Space-Vector Modulation with Variable-Structure Control. IEEE Trans. Power Electron. 2003, 18, 1070–1083. [Google Scholar] [CrossRef] [Green Version]
  10. Pan, C.-T.; Liao, Y.-H. Modeling and Control of Circulating Currents for Parallel Three-Phase Boost Rectifiers with Different Load Sharing. IEEE Trans. Ind. Electron. 2008, 55, 2776–2785. [Google Scholar] [CrossRef]
  11. Ye, Z.; Boroyevich, D.; Lee, F.C. Modeling and Control of Zero-Sequence Current in Parallel Multi-Phase Converters. In Proceedings of the 2000 IEEE 31st Annual Power Electronics Specialists Conference (Cat. No.00CH37018), Galway, Ireland, 23 June 2000; Volume 2, pp. 680–685. [Google Scholar]
  12. Shao, Z.; Zhang, X.; Wang, F.; Cao, R. Modeling and Elimination of Zero-Sequence Circulating Currents in Parallel Three-Level T-Type Grid-Connected Inverters. IEEE Trans. Power Electron. 2015, 30, 1050–1063. [Google Scholar] [CrossRef]
  13. Zhang, P.; Zhang, G.; Du, H. Circulating Current Suppression of Parallel Photovoltaic Grid-Connected Converters. IEEE Trans. Circuits Syst. II Express Briefs 2018, 65, 1214–1218. [Google Scholar] [CrossRef]
  14. Ogasawara, S.; Takagaki, J.; Akagi, H.; Nabae, A. A Novel Control Scheme of a Parallel Current-Controlled PWM Inverter. IEEE Trans. Ind. Appl. 1992, 28, 1023–1030. [Google Scholar] [CrossRef]
  15. Ye, Z.; Boroyevich, D.; Choi, J.-Y.; Lee, F. Control of Circulating Current in Two Parallel Three-Phase Boost Rectifiers. IEEE Trans. Power Electron. 2002, 17, 609–615. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, T. Dual-Modulator Compensation Technique for Parallel Inverters Using Space-Vector Modulation. IEEE Trans. Ind. Electron. 2009, 56, 3004–3012. [Google Scholar] [CrossRef]
  17. Prasad, J.S.S.; Ghosh, R.; Narayanan, G. Common-Mode Injection PWM for Parallel Converters. IEEE Trans. Ind. Electron. 2015, 62, 789–794. [Google Scholar] [CrossRef]
  18. Hou, C. A Multicarrier PWM for Parallel Three-Phase Active Front-End Converters. IEEE Trans. Power Electron. 2013, 28, 2753–2759. [Google Scholar] [CrossRef]
  19. Liberos, M.; González-Medina, R.; Garcerá, G.; Figueres, E. Modelling and Control of Parallel-Connected Transformerless Inverters for Large Photovoltaic Farms. Energies 2017, 10, 1242. [Google Scholar] [CrossRef] [Green Version]
  20. Figueres, E.; Garcera, G.; Sandia, J.; Gonzalez-Espin, F.J.; Calvo, J.; Vales, M. Dynamic Analysis of Three-Phase Photovoltaic Inverters with a High Order Grid Filter. In Proceedings of the 2007 IEEE International Symposium on Industrial Electronics, Vigo, Spain, 4–7 June 2007; pp. 691–696. [Google Scholar]
  21. Mohd, A.; Ortjohann, E.; Hamsic, N.; Sinsukthavorn, W.; Lingemann, M.; Schmelter, A.; Morton, D. Control Strategy and Space Vector Modulation for Three-Leg Four-Wire Voltage Source Inverters under Unbalanced Load Conditions. IET Power Electron. 2010, 3, 323–333. [Google Scholar] [CrossRef]
  22. Albatran, S.; Fu, Y.; Albanna, A.; Schrader, R.; Mazzola, M. Hybrid 2D-3D Space Vector Modulation Voltage Control Algorithm for Three Phase Inverters. IEEE Trans. Sustain. Energy 2013, 4, 734–744. [Google Scholar] [CrossRef]
  23. Saleh, K.; Sumner, M. A SVM-3D Based Encoderless Control of a Fault-Tolerant PMSM Drive. Electronics 2020, 9, 1095. [Google Scholar] [CrossRef]
  24. Abdali Nejad, S.; Matas, J.; Elmariachet, J.; Martín, H.; de la Hoz, J. SOGI-FLL Grid Frequency Monitoring with an Error-Based Algorithm for a Better Response in Face of Voltage Sag and Swell Faults. Electronics 2021, 10, 1414. [Google Scholar] [CrossRef]
Figure 1. Circuit of the parallel inverters.
Figure 1. Circuit of the parallel inverters.
Electronics 11 03720 g001
Figure 2. Averaged equivalent circuit of the parallel inverters.
Figure 2. Averaged equivalent circuit of the parallel inverters.
Electronics 11 03720 g002
Figure 3. Scheme of the proposed control stage.
Figure 3. Scheme of the proposed control stage.
Electronics 11 03720 g003
Figure 4. Spatial representation of the 3D-SVM switching vectors. (a) Two tetrahedra composed of two consecutive nonzero switching vectors and neutral vectors v0 and v7 (b) Spatial representation of the twelve tetrahedra.
Figure 4. Spatial representation of the 3D-SVM switching vectors. (a) Two tetrahedra composed of two consecutive nonzero switching vectors and neutral vectors v0 and v7 (b) Spatial representation of the twelve tetrahedra.
Electronics 11 03720 g004aElectronics 11 03720 g004b
Figure 5. Single-line scheme of two inverters connected in parallel.
Figure 5. Single-line scheme of two inverters connected in parallel.
Electronics 11 03720 g005
Figure 6. Experimental set-up.
Figure 6. Experimental set-up.
Electronics 11 03720 g006
Figure 7. Predicted and measured current control loops gain: (a) d-axis, (b) q-axis, and (c) o-axis. (d) Measured current control loop gain in the o-axis where *f stands for the crossover frequency, and *R and *Ɵ for the gain and the phase of the control loop at this frequency.
Figure 7. Predicted and measured current control loops gain: (a) d-axis, (b) q-axis, and (c) o-axis. (d) Measured current control loop gain in the o-axis where *f stands for the crossover frequency, and *R and *Ɵ for the gain and the phase of the control loop at this frequency.
Electronics 11 03720 g007aElectronics 11 03720 g007b
Figure 8. Pole-zero map considering several Lf values: (a) d-axis, (b) q-axis and (c) o-axis.
Figure 8. Pole-zero map considering several Lf values: (a) d-axis, (b) q-axis and (c) o-axis.
Electronics 11 03720 g008
Figure 9. Simulation results of the control of the circulating currents of two inverters connected in parallel (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor between inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter.
Figure 9. Simulation results of the control of the circulating currents of two inverters connected in parallel (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor between inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter.
Electronics 11 03720 g009aElectronics 11 03720 g009b
Figure 10. Simulation results of the control of the circulating currents among three inverters connected in parallel (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor among inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter.
Figure 10. Simulation results of the control of the circulating currents among three inverters connected in parallel (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor among inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter.
Electronics 11 03720 g010aElectronics 11 03720 g010b
Figure 11. Experimental results of the control of the circulating currents (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor between inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter. * stands for values that the oscilloscope did not calculate. In this case, the RMS value of Channel 8.
Figure 11. Experimental results of the control of the circulating currents (a) with unbalanced inductances in both inverters, (b) with unbalanced load factor between inverters, (c) with unbalanced inductances in the same inverter, (d) with a different modulation in each inverter. * stands for values that the oscilloscope did not calculate. In this case, the RMS value of Channel 8.
Electronics 11 03720 g011
Table 1. Expression of the variables in the operation point.
Table 1. Expression of the variables in the operation point.
VariableExpression
VgdVgRMS (phasephase)
Vgq0
Vgo0
Idi I d c _ n o m n · F m · D d i c i ,   i = 1 , 2 , n
Iqi0
Ioi0
Ddi V g d F m · V d c
Dqi ω I d i L f + L f g M f M f g + j = 1 n c j c i L g M g Fm · Vdc ,   i = 1 , 2 , n
Doi0
Table 2. 3D-SVM switching vectors.
Table 2. 3D-SVM switching vectors.
vφv0-000v1-100v2-110v3-010v4-011v5-001v6-101
v A 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c
v B 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c
v C 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c 1 2 v d c
v α 0 2 3 v d c 1 6 v d c 1 6 v d c 2 3 v d c 1 6 v d c 1 6 v d c
v β 00 1 2 v d c 1 2 v d c 0 1 2 v d c 1 2 v d c
v O 3 2 v d c 1 2 3 v d c 1 2 3 v d c 1 2 3 v d c 1 2 3 v d c 1 2 3 v d c 1 2 3 v d c
Table 3. Sequence of switching vectors in each prism.
Table 3. Sequence of switching vectors in each prism.
PrismSequence of Switching Vectors
Iv7-v2-v1-v0-v1-v2-v7
IIv7-v2-v3-v0-v3-v2-v7
IIIv7-v4-v3-v0-v3-v4-v7
IVv7-v4-v5-v0-v5-v4-v7
Vv7-v6-v5-v0-v5-v6-v7
VIv7-v6-v1-v0-v1-v6-v7
Table 4. Parameters of the inverters.
Table 4. Parameters of the inverters.
ParameterNominal ValueParameterNominal Value
Vg-RMS (phase–phase)230 Vrf50 mΩ
Vdc500 Vrfg50 mΩ
Pn5 kWrg50 mΩ
Co1.2 mFCf9 µF
Lf5 mHRd4.4 Ω
Lfg0fsw10 kHz
Lg320 µHRs1 V/A
Mf0β1 V/V
Mfg0Fm0.5 V/V
Mg−80 µH
Table 5. Instruments and laboratory material.
Table 5. Instruments and laboratory material.
AC voltage sourceCinergia GL and EL-50
DC voltage sourceGSS Regatron
Current probesFluke i400s
Differential voltage probesYokogawa 700924
Frequency response analyzer (FRA)NF Corporation FRA5097
OscilloscopeAgilent MSO6014A (4 channels, 200 MHz)
Yokogawa DLM4038 (8 channels 350 MHz)
MultimeterAgilent 1241A, Fluke 175, Metrix MX0024
Thermal cameraFluke Ti25 IR FUSION TECHNOLOGY THERMAL IMAGER
Table 6. Real values of the inductance Lf of 5 mH.
Table 6. Real values of the inductance Lf of 5 mH.
ParameterReal ValueParameterReal Value
Lf_a15.14 mHLf_a25.1 mH
Lf_b15.14 mHLf_b24.85 mH
Lf_c15.27 mHLf_c25.03 mH
Table 7. Real values of the inductance Lf of 7 mH.
Table 7. Real values of the inductance Lf of 7 mH.
ParameterReal Value
Lf_a27.16 mH
Lf_b26.98 mH
Lf_c27.12 mH
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liberos, M.; González-Medina, R.; Patrao, I.; Garcerá, G.; Figueres, E. A Control Scheme to Suppress Circulating Currents in Parallel-Connected Three-Phase Inverters. Electronics 2022, 11, 3720. https://doi.org/10.3390/electronics11223720

AMA Style

Liberos M, González-Medina R, Patrao I, Garcerá G, Figueres E. A Control Scheme to Suppress Circulating Currents in Parallel-Connected Three-Phase Inverters. Electronics. 2022; 11(22):3720. https://doi.org/10.3390/electronics11223720

Chicago/Turabian Style

Liberos, Marian, Raúl González-Medina, Iván Patrao, Gabriel Garcerá, and Emilio Figueres. 2022. "A Control Scheme to Suppress Circulating Currents in Parallel-Connected Three-Phase Inverters" Electronics 11, no. 22: 3720. https://doi.org/10.3390/electronics11223720

APA Style

Liberos, M., González-Medina, R., Patrao, I., Garcerá, G., & Figueres, E. (2022). A Control Scheme to Suppress Circulating Currents in Parallel-Connected Three-Phase Inverters. Electronics, 11(22), 3720. https://doi.org/10.3390/electronics11223720

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