Next Article in Journal
Economic Dispatch of BESS and Renewable Generators in DC Microgrids Using Voltage-Dependent Load Models
Previous Article in Journal
Factors Affecting Fuelwood Consumption and CO2 Emissions: An Example from a Community-Managed Forest of Nepal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Ripple Current in the Capacitors of Active Power Filters

Institute of Control and Industrial Electronic, Warsaw University of Technology, 75 Koszykowa St., 00-662 Warsaw, Poland
Energies 2019, 12(23), 4493; https://doi.org/10.3390/en12234493
Submission received: 18 October 2019 / Revised: 17 November 2019 / Accepted: 18 November 2019 / Published: 25 November 2019

Abstract

:
This article provides formulae to determine the root mean square (rms) value of a capacitor current in an inductive-capacitive-inductive (LCL) filter used in a parallel active power filter (PAPF) circuit. The article presents an analysis of three components of the capacitor current: a component forced by the usually distorted voltage of the grid; a component forced by the nonlinear load current harmonics and harmonics in the output current of the PAPF that compensates them (a novel aspect presented in this document); and a component forced by the inverters of the PAPF containing carrier and sideband harmonics. The article also presents formulae for determining the rms value of current harmonics in dc-link capacitors forced by the ripple of the ac output current without load of the filter inverters (also novel to this document). The results of the analysis have been confirmed by simulation and experimental research of a commercial active filter consisting of two parallel interleaved voltage inverters. Elements of the LCL filter of the PAPF have been selected according to dependencies available in scientific and technical literature. In addition, the formulae presented in the article are used to verify the correctness of selection of capacitors from the point of view of their catalogue acceptable rms value of capacitor current.

Graphical Abstract

1. Introduction

The parallel active power filter (PAPF) contains three groups of passive power elements: dc-link capacitors, inductances, and capacitors in the LCL-ripple filter. Specified in the data sheet, the permissible capacitor rms current is one of its basic parameters. Liserre et al. [1] propose a step-by-step procedure for designing the LCL filter of a front-end three-phase active rectifier. Jeong et al. in [2] and Asiminoei et al. in [3] give formulae allowing for calculation of rms current ripple in the LCL inductors. In [4], Jalili and Bernet applied Bessel functions to select the LCL filter parameters for defined maximum grid current harmonics.
This article presents an analysis leading to the determination of the rms current of LCL-ripple filter capacitors used at the output of the PAPF. There are three components of this current:
(1)
Component forced by the usually distorted mains voltage;
(2)
Component forced by the nonlinear load current harmonics and harmonics in the output current of the PAPF that compensates them;
(3)
Component forced by the inverters of the PAPF and containing carrier and sideband harmonics (for equal and different inductances on the converter side).
In systems with a single inverter, the dominant harmonic frequency of the output voltage is equal to the frequency of the auxiliary carrier harmonic. Sahoo et al. [5] give the relation determining the rms current of the LCL filter capacitor in a one-grid connected inverter for the first component (for an undistorted grid) and the third component of this current. In PAPF systems with two parallel interleaved inverters, “carrier harmonics” disappear; thus, the share of the second component in the rms current of the capacitor increases.
Studies describing the rms current of the LCL-capacitor ripple of output filter installed at the output of parallel interleaved inverters are not common in the scientific literature. Tang et al. in [6] and Vodyakho et al. in [7] analyze the impact of capacitance values of LCL filters on the stability of the PAPF.
This article presents the equations describing the rms value of the second component, independent of the number of parallel interleaved inverters, and the rms value of the third component for two parallel interleaved inverters. Formulae are given for the PAPF containing a transformer with a Dy5 connection group.
The ripple currents in the LCL capacitors, LCL inductors, and dc-link capacitors depend on each other [3]. This article provides formulae to determine the rms value of the ripple current in dc-link capacitors, which comes from the ripple current in the LCL inductors. McGrath et al. in [8], Kolar et al. in [9], Zhang et al. in [10,11], provide an algorithm for calculating the value of the harmonic current in the dc-link capacitor for the voltage source when the load current is assumed to be a fundamental single-frequency sine. The analysis presented in these works, applied to the PAPF system with the deactivated reactive power compensation for the basic harmonic, will not reveal the harmonic values determined in this article, even if it takes into account the 5-, 7-, 11-, and 13-order ac harmonic in the output current. Zhang et al. in [11], Ye et al. in [12], and Quan at al. [13] proposed an interleaving scheme to reduce dc-link current harmonics, but the optimization of the interleaving angle from the point of view of harmonic currents in the dc-circuit also influences the ripple current in the ac capacitors.
Many articles publish theoretical analysis and results of simulation and experimental tests of current components in LCL filter capacitors and dc-link circuit. However, they do not provide dependencies which can be easily used to easily calculate the rms values of currents in these capacitors. This article gives simplified equations with an error of about 10% on the basis of which one can determine both the rms value of the capacitor current in the dc-link circuit and the three current components of the LCL filter capacitor.
In the present paper, the rms values of currents in LCL filter capacitors and dc-link capacitors was determined analytically for the control system of the filter output current with asymmetrical regular-sampled PWM. In addition, simulation tests for control with symmetrical regular-sampled PWM have been conducted.

2. Description of the Input Circuit of PAPF

Figure 1a,b shows the schematic diagrams of the analyzed medium-voltage power supply system. A controlled twelve-pulse converter is powered via a transformer Tr2 with a vector group of Dd0y5 connections. A PAPF containing two voltage inverters and a medium voltage transformer Tr1 with a vector group of Dy5 compensates harmonic current drawn by the non-linear load in the form of a 12Th converter (Th1–Th12). Two capacitors with capacitance C/2 at the output of each phase function as a capacitor with capacitance C of the LCL filter reducing common mode voltage (CMV) [14]. This connection of capacitors ensures the free flow of harmonics of their zero sequence currents to the dc-circuit. The chokes L_1, L_2 are installed at the output of the inverters and, together with the capacitors with the resultant capacitance C and the leakage inductance of the transformer Tr1, create an LCL type filter [6,15,16,17,18,19] in each phase. Rd is LCL-filter resistor damping.
The PAPF has two inverters (Q1_1–Q6_2) [3] with FF600R12IS4F hybrid modules (combination of Insulated Gate Bipolar Transistor (IGBT) with silicon carbide Schottky diodes).
Voltages and currents ug, u, ig, iL, is, ip, ic, io_1, io_2, and i shown in Figure 1a,b are vectors: ug = [ug1,ug2,ug3]T; u = [u1,u2,u3]T; uc = [uc1,uc2,uc3]T; ig = [ig1,ig2,ig3]T; iL = [iL1,iL2,iL3]T; is = [is1,is2,is3]T; ip = [ip1,ip2,ip3]T; ic = [ic1,ic2,ic3]T; io_1 = [io1_1,io2_1,io3_1]T; io_2 = [io1_2,io2_2,io3_2]T; and i = [i1,i2,i3]T.
Table 1 shows the relation of the primary and secondary vectors of the Dy5 transformer voltage (where υz is the transformer turns ratio). Vectors ug, u, ig, iL, and is are values brought to transformer Tr1 primary side terminals 1, 2, 3. Reference [19] gives the relations between the primary and the secondary vectors of the Dy11 transformer voltage.
The algorithm described in the literature [20,21] is implemented in the control system with a voltage regulator in the dc-circuit. The control block diagram, criteria for selection of the current regulator gain factor, and criteria for determining the inductance L at the outputs of the PAPF (L is the nominal value of L_1 and L_2) are presented in [22]. The analysis carried out in this article does not influence the impact of PAPF control methods on the quality of compensation of higher harmonics of non-linear load current.
Commonly applied in power electronics, a TMS320F28335 Digital Signal Processor was used to implement the chosen filter control.

3. Inverter Output Voltage with Asymmetrical Regular-Sampled PWM

Equations describing the output voltage of the transistor branch uoi_j produced in the system with asymmetrical regular-sampled PWM are given in [23]:
u o i _ j = n = 1 F 0 n cos ( n ω g t + n θ g i ) + m = 1 n = F m n cos [ ( m ω C / ω g + n ) ω g t + m θ C i _ j + n θ g i ]
where indexes i ϵ {1, 2, 3}, j ϵ {1, 2} denote the inverter’s branch and number, respectively; θ g i , θCi_j are the phase angles of the fundamental harmonic of the mains voltage ug brought to transformer Tr1 primary side terminals 1, 2, 3, and a carrier wave, respectively; ωC is the pulsation of the triangular carrier; ωg is the angular frequency of the fundamental harmonic of line voltage ug; and m is the carrier index variable. Harmonic coefficients F m n are given in [23] by:
F m n = 2 U d c π J n ( [ m + n ω g ω C ] π M 2 ) m + n ω g ω C sin [ ( m + n ) π 2 ]
where Jn(x) denotes a Bessel function of the first type, with order n and argument x; U d c is dc-link voltage; M is modulation depth ( M 2 U ^ g , 1 / U d c , 0 ≤ M ≤ 1); and U ^ g , 1 is the amplitude of the fundamental harmonic voltages u g 1 , 1 , u g 2 , 1 and u g 3 , 1 . The form of Equation (2) indicates that Fmn coefficient takes a zero value if the conditions of m + n = 0, ±2, ±4… are fulfilled.
From Equations (1) and (2) it can be concluded that the carrier and sideband harmonics of the same order set for θCi_j = 0 or θCi_j = π have the same amplitudes and their phase values are equal or differ by π rad. If the active filter includes two inverters connected as shown in Figure 1a, from the point of view of the higher harmonics in the output current of the filter it is advantageous to offset the carriers from each other in modulator systems that control the transistor in the respective phase legs of the inverters by an angle π rad [3,24]. After insertion of θCi_1 into Equation (1) for j = 1 and then substituting θCi_2 = θCi_1 + π in Equation (1) for j = 2, we can obtain the following:
u o i _ 1 = n = 1 U o _ 1 , 0 n cos ( n ω g t + n θ g i ) + m = 1 n = U o _ 1 , m n cos { [ m ( ω C / ω g ) + n ] ω g t + m θ C i _ 1 + n θ g i }
u o i _ 2 = n = 1 U o _ 2 , 0 n cos ( n ω g t + n θ g i ) + m = 1 n = U o _ 2 , m n cos { [ m ( ω C / ω g ) + n ] ω g t + m θ C i _ 1 + n θ g i }
Harmonic coefficients U o _ 1 , m n and U o _ 2 , m n are given by
U o _ 1 , m n = F m n
U o _ 2 , m n = F m n ( 1 ) m
If we assume that ωC/ωg ≥ 40, then the coefficients of subharmonics U o _ 1 , m n or U o _ 2 , m n present for non-integer values of quotient ωC/ωg values will be very small and their impact on the effective value of the capacitor current will be negligible. Therefore, further considerations can be carried out assuming that the waveforms uoi_1 and uoi_2 are the periodic functions, which means that the size of ξ, defined by formula ξ = ωC/ωg is an integer. Harmonic coefficients U o _ 1 , m n denoted by U o _ 1 , k , U o _ 2 , m n will also be denoted by U o _ 2 , k where k = + n. The choice of ωC/ωg ≥ 40 is justified by the limit values of inequalities (17) given in Section 4.3. Figure 2 shows the harmonics spectrum of the output voltages expressed in the p.u. system of each inverter with asymmetrical regular-sampled PWM. Amplitudes of voltage harmonic coefficients have been determined in relation to half the voltage of the dc-link circuit ( U ^ o , k ( p u ) = 2 U ^ o , k / U d c ) .
The spectrum of uoi_j shows that significant harmonics of the output currents of the branches of inverters occur for m = 1 if n = ±2, and for m = 2 if n = ±1, ±3 (∣n∣ ≤ 3).

4. LCL System Capacitor Current

The analysis in this article will be carried out assuming full symmetry of both the load and the three-phase supply voltage.
For clarity of the following considerations, the vector magnitudes were replaced by their components without indices “1”, “2”, and “3”.
In the power supply system with the PAPF containing a transformer, determining the rms values of the LCL-ripple filter capacitor current requires calculation of amplitudes and initial phases of mains voltage harmonics ug,k and of load current iL,k brought to transformer Tr1 primary side terminals 1, 2, and 3.
The following relations can describe the k-th order harmonic of voltage ug,k, and current iL,k:
u g , k = U ^ g , k sin ( k ω g t + ψ g , k )
i L , k = I ^ L , k sin ( k ω g t + φ L , k )
where ψ g , k , φ L , k , U ^ g , k , I ^ L , k are the initial phases and amplitudes of the k-th harmonic voltage u g , k and current i L , k , respectively.
The formulae presented in Table 1 allow determination of the amplitudes ( U ^ g , k , I ^ L , k ) and the initial phases ( ψ g , k , φ L , k ) of the harmonic voltage u g , k and current i L , k , respectively, for the system with a transformer of a group of Dy5 connections. These relationships are given in Table 2.
We can select three components of the LCL filter capacitor current based on the sources that force them:
ic = icG + icL + icC
where the icG component is forced by the usually distorted voltage of mains supply; the icL component is forced by the current harmonics of nonlinear load and the harmonics in the output current of the PAPF that compensates them and the component icC is also forced by the PAPF and contains the carrier (for L_1L_2) and sideband harmonics.
In real conditions, the number of measured (by typical Total Harmonic Distortion measurement) harmonics does not exceed 40. For the determination of the first two components of the rms current of the capacitor, it was assumed that the voltage inverters with inductances L_1 and L_2 at the output can be replaced with current sources for harmonics of an order below 40.

4.1. Rms of Capacitor C Current Component I c G , r m s

The frequency spectrum of the capacitor current component i c G contains fundamental and baseband harmonics k = n = 1,2… (m = 0).
The first of these components can be determined from the equivalent circuit shown in Figure 3a. Based on this diagram, we can determine the amplitude I ^ c G , k :
I ^ c G , k = U ^ g , k [ k ω g ( L l q + L g ) 1 / ( k ω g C ) ] 2 + R d 2
where L l q is the dispersion inductance Tr1, and L g is the grid-side inductance related to the low voltage side of Tr1. The damping resistor Rd was determined on the basis of the relation given in [25]:
R d = 2 ζ p ( L l q + L g ) L [ 2 ( L l q + L g ) + L ] C
where ζp is the resonant-pole damping factor, and L is the nominal value of L_1 and L_2 (L_1L_2).
The damping factor has been considered critical in some theoretical works [26].

4.2. Rms of Capacitor C Current Component I c L , r m s

The frequency spectrum of the capacitor current component icL contains fundamental and baseband harmonics of the order k = n = 1, 2… (m = 0).
The amplitudes of particular harmonics of the second component can be determined on the basis of the scheme shown in Figure 3b:
I ^ c L , k = | k 2 ω g 2 C [ I ^ k ( L g + L l q ) I ^ L , k L g ] [ 1 k 2 ω g 2 C ( L l q + L g ) ] 2 + ( k ω g C R d ) 2 |
where I ^ L , k and I ^ k are the amplitudes of the k-order harmonic of the load current and the component of the filter output current compensating this harmonic, respectively. If the harmonic attenuation factor λ of the k-th harmonic of load current ( λ k = I ^ k / I ^ L , k ) is known, we obtain:
I ^ c L , k = | k 2 ω g 2 C [ λ k ( L g + L l q ) L g ] [ 1 k 2 ω g 2 C ( L l q + L g ) ] 2 + ( k ω g C R d ) 2 | I ^ L , k
λk is usually a set point parameter in the PAPF.
Harmonic components of icG and icL may be of the same order. The resulting harmonic of this order depends on the amplitude and initial phase of harmonic components. In practical systems, the measurement of the amplitude and the initial phase of mains voltage harmonics should be made with the load and PAPF switched off.
Current icG,k forced by the k-th harmonic voltage ug,k is capacitive if k satisfies the following inequality:
k < 1 / ( ω g ( L l q + L g ) C )
The rms value of IcGL,rms, considering the icG and icL components, is described by the formula
I c G L , r m s = k = 1 [ 0.5 ( I ^ c G , k 2 + I ^ c L , k 2 ) ± s i g n [ λ k ( L g + L l q ) L g ] I ^ c G , k I ^ c L , k sin ( Δ φ k ) ]
where Δ φ k   =   ψ g , k φ L , k     =   ψ g , k φ L , k . In Equation (13) “+” symbolizes the inductive and “−” the capacitive nature of the current icG,k (respectively). If the initial phases ψ g , k and φ L , k are not known the equation used is
I c G L , r m s 0.5 k = 1 ( I ^ c G , k + I ^ c L , k ) 2

4.3. Rms of Capacitor C Current Ripple I c C , r m s

The frequency spectrum of the capacitor ripple current contains carrier harmonics of the k = ξ + n order (for L_1L_2) and sideband harmonics of the k = + n (m = 2, 4, 6…) order.
Figure 4a,b show carrier waveforms for two PWM techniques: (Figure 4a) double interleaved PWM carriers [3] and (Figure 4b) double three interleaved PWM carriers [13,22].
The duration time of the microprocessor interruption service routine is denoted by τw (PWM computation delays [22,27]). The choice of the PWM technique is justified by the results of the simulation research collected in Section 6.3. These studies show that double three interleaved PWM carriers ensure minimum values of common mode voltage (CMV), minimum values of the number of operations in the microprocessor interruption service routine (N), and minimum values of rms current dc-circuit capacitors in the near-full range of the depth modulation index M.
After insertion of θC1_1 in Equation (1) and then substituting phases angle θC2_1 = θC1_1 −2π/3, θC3_1 = θC1_1 + 2π/3, θCi_2 = θCi_1 + π (for i = 1, 2, 3), θ g 1 = 0, θ g 2 = 2π/3, θ g 3 = 2π/3, respectively, we obtain six equations, from which the following conclusion can be drawn: If the fundamental harmonics u g 1 , 1 , u g 2 , 1 ,   u g 3 , 1 of voltages and carrier signals of the first, second and third branch of each inverter create the same, three-phase positive sequences, then the individual harmonics of the output voltages of each inverter (uo1_1, uo2_1, uo3_1 and uo1_2, uo2_2, uo3_2) make a positive, negative or zero sequence if the conditions in Table 3 are met. The same rules also apply if the fundamental harmonics u g 1 , 1 , u g 2 , 1 ,   u g 3 , 1 of voltages and carrier signals simultaneously form the same, three-phase negative sequences.
If we assume θC2_1 = θC1_1 + 2π/3, θC3_1 = θC1_1 − 2π/3, θCi_2 = θCi_1 + π (for i = 1, 2, 3), θ g 1 = 0, θ g 2 = −2π/3, θ g 3 = 2π/3 then the fundamental harmonics u g 1 , 1 , u g 2 , 1 ,   u g 3 , 1 of voltages form a three-phase positive sequence and carrier signals form a three-phase negative sequence. In the case of the opposite three-phase sequences u g 1 , 1 , u g 2 , 1 ,   u g 3 , 1 of voltages to the sequence of carrier signals, the other rules are presented in Table 3. In this article, the relations for which the sequence of fundamental harmonics u g 1 , 1 , u g 2 , 1 ,   u g 3 , 1 of voltages opposite to the sequence of carrier signals will be marked by superscript (op).
The same rules apply to the harmonic of capacitor currents ic1, ic2, and ic3 and output currents of branches of each inverter io1_1, io2_1, io3_1, io1_2, io2_2, and io3_2, because the harmonic of the output voltages of each inverter uo1_1, uo2_1, uo3_1 and uo1_2, uo2_2, and uo3_2 directly enforce the harmonics of these currents.
Figure 5a shows the equivalent circuit of the output circuit of the respective branches of the two inverters, based on which the ÎcC,k amplitude of the k-th harmonic of the C capacitor current can be determined if the harmonics of the currents icC1,k, icC2,k and icC3,k form three-phase positive or negative sequences. As harmonics with three-phase zero sequence of io1,k and io2,k currents cannot flow in windings without the neutral wire of the transformer Tr1, the equivalent circuit shown in Figure 5a can be simplified for harmonics with zero sequence as shown in Figure 5b.
The aim of the discussion conducted below is to design a relation determining the rms value of the icC component of the current of the C capacitor. Figure 6 shows the harmonic spectrum of the current component icC expressed in the p.u. system in the LCL circuit in a system with two inverters.
The amplitude of harmonic ÎcC,k(pu) is described by the relation: ÎcC,k(pu) = ÎcC,k/Ib. Base current Ib is defined by relation I b = U d c / ( 8 L f C ) [17]. Ib is equal to the maximum output current ripple of each branch of the inverter for uci = 0, where i ϵ {1, 2, 3}.
The capacitance C and inductors L_1 and L_2 of L nominal inductance are selected so that the resonant frequency fres of the LCL system is given by
f r e s = 1 2 π L + 2 ( L l q + L g ) L ( L l q + L g ) C
and approximately satisfies the conditions given in [28,29]
10 f g < f r e s < 0.5 f C
where fg = ωg/(2π) and fC is the carrier frequency. When the value of the inductances of the chokes L_1 and L_2 are within the tolerance of 5%, then the value of fres determined by Equation (15) yields an error lower than ±2.5%.
The lower limit of the condition in Equation (16) was adopted for an LCL filter in a three-phase grid-connected inverter [28]. For LCL filters used in an PAPF, the lower limit given in [6] depends on the maximum value of the compensated harmonic order (kmax) and is 4 kmax fg, so is usually greater than 10 fg.
4 k max f g < f r e s < 0.5 f C
Further analysis will repeatedly be based on simplifications resulting from a sufficiently high value of ξ = fC/fg. The required value ξ = 8 kmax can be determined from the limit values of fres. If we assume that the three-wire three-phase PAPF compensates only the fifth and lower harmonics, then the least favorable value of fC/fg = 40 is obtained (from the point of view of analysis accuracy). This value is in accordance with the assumption in Section 3 (fC/fg ≥ 40).
In systems consisting of parallel interleaved inverters the upper limit of fres given by Equations (16) and (17) can have higher values [30]. This is due to the fact that for L_1 = L_2 the minimum frequency of the higher harmonics of the LCL filter capacitor current is approximately equal to the double frequency of the carrier harmonic fC (Figure 6). On the other hand, with the limited accuracy of the L_1 and L_2 inductances, a harmonic of frequency fC may appear in the current of the capacitor C. As can be seen from calculations supported by simulation research, when the values of L_1 and L_2 are held to within 5%, the value of amplitude ÎcC,ξ for modulation depth M = 0.9 is close to 45% of the value of amplitude ÎcC,2ξ − 1. The influence of the harmonic of the carrier frequency on the risk of creating a current resonance in the LCL filter is high. Due to the above, further analysis of the effective current in the LCL filter capacitor assumes that Equation (17) is fulfilled. Meeting the upper limit of the inequality of Equation (17) is equivalent to meeting the fC/fres > 2 condition.
The coefficients of current components with a positive and negative harmonic sequence I c C , k ( p , n ) and with zero harmonic sequence I c C , k ( 0 ) can be determined on the basis of the equivalent circuit of the diagrams shown in Figure 5a,b respectively. Resulting from Equations (3a), (3b) (4a) and (4b) the harmonics of voltages uo_1,k and uo_2,k for m = 2, 4, 6… are equal. The time waveform current harmonic i c C , i , k is described by:
i c C , i , k = I c C , k cos [ k ω g t + m θ C i _ 1 + n θ g i + φ c , k ] for m = 2 , 4 , 6
where φ c , k is the phase shift for the k-th harmonic of the current icC,k forced by the circuit’s impedance.
In order to determine the relationships describing the harmonic coefficients IcC,k, it is convenient to present the current icC,k, voltages uo_1,k and uo_2,k and impedances in circuits in complex forms I _ c C , k = I c C , k e j ( π / 2 + m θ C 1 _ 1 + φ c , k ) , U _ o _ 1 , k = U o _ 1 , k e j ( π / 2 + m θ C 1 _ 1 + n θ g ) ,   U _ o _ 2 , k = U o _ 2 , k e j ( π / 2 + m θ C 1 _ 1 + n θ g ) , X _ 1 = j k ω g L _ 1 , X _ 2 = j k ω g L _ 2 , X _ 3 = j k ω g ( L q l + L g ), and X _ c = 1 / ( j k ω g C ) , and then apply the Thevenin method. After a few transformations we get
I c C , k ( p , n ) = U o _ 1 , k L _ 2 + U o _ 2 , k L _ 1 L _ 2 + L _ 1 L l q + L g L e q + L l q + L g 1 Z c , k for m ± n = ± 1 , ± 2 , ± 4 , ± 5 , ± 7
I c C , k ( 0 ) = U o _ 1 , k L _ 2 + U o _ 2 , k L _ 1 L _ 2 + L _ 1 1 Z c , k ( 0 ) for m ± n = 0 , ± 3 , ± 6
where Leq = L_1L_2/(L_1 + L_2) while Zc,k and Z c , k ( 0 ) are the equivalent impedances of the circuit for positive and negative sequence harmonics (Zc,k) and zero sequence harmonics ( Z c , k ( 0 ) ), respectively, given by
Z c , k = [ k ω g ( L l q + L g ) L e q L e q + L l q + L g 1 k ω g C ] 2 + R d 2
Z c , k ( 0 ) = [ k ω g L e q 1 / ( k ω g C ) ] 2 + R d 2
Equations (19a) and (19b) is true for systems in which voltages u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 and the voltages u T 1 _ 1 , 1 , u T 2 _ 1 , 1 , and u T 3 _ 1 , 1 form the same or the opposite sequence. In the first case, the conditions in the formulae are related to the sum of m + n, in the second case to the difference of mn.
For m ≥ 2, ξ ≥ 40, and ∣n∣ ≤ 3, the approximation + n has an error of no more than 3.75%. After taking the approximation k = and taking into account in Equations (20a) and (20b) the relations given in Equations (9) and (15), as well as ξ = fC/fg, we obtain
Z c , m ξ = m ω C ( L l q + L g ) L e q L e q + L l q + L g [ 1 f r e s 2 m 2 f C 2 ] 2 + ( 2 ξ p f r e s m f C ) 2
Z c , m ξ ( 0 ) = m ω C L e q [ 1 f r e s 2 m 2 f C 2 ( 1 L e q L Σ ) ] 2 + [ 2 ξ p f r e s m f C ( 1 L e q L Σ ) ] 2
where L = Leq + Llq + L g .
Reference [6] recommends for a PAPF equal inductances on the grid and converter-side to produce the lowest resonance frequency, while in [28] it is recommended for the three-phase active rectifier that the ratio between converter-side and grid-side inductances is in the range of 3–7. If we assume the limit values L/Leq = 2, m = 2, fC/fres = 2, ζp = 0, then the value of the square root in Equation (21b) maximally deviates from unity and reaches the value of 0.97. The value of the square root in Equation (21a) maximally deviates from the unity for m = 2, fC/fres = 2, ζp = 0 and is equal to 0.94. Therefore, for m ≥ 2, it is possible to determine Z c , m ξ and Z c , m ξ ( 0 ) with errors not exceeding 6% and 3%, respectively, using the following dependencies:
Z c , m ξ = m ω C ( L l q + L g ) L e q L e q + L l q + L g
Z c , m ξ ( 0 ) = m ω C L e q
The harmonic spectrum presented in Figure 6 for an PAPF with two inverters shows that the value ÎcC depends on harmonics whose order is close to the value of 2ξ. To calculate the rms value of current I c C , 2 ξ ( r m s ) ( p , n ) , we need to find the square root of the sum of its squares for n = −3, −1, 3. Using Equation (19a) and assuming U o _ 2 , 2 ξ + n = U o _ 1 , 2 ξ + n we obtain the relation, without introducing a significant error, for the value of ξ ≥ 40:
I c C , 2 ξ ( r m s ) ( p , n ) = 2 2 Z c , 2 ξ L l q + L g L e q + L l q + L g n = 3 , 1 , 3 U o _ 1 , 2 ξ + n 2
If L_2L_1 then IcC,ξ is different than zero. By substituting U o 2 , ξ = U o _ 1 , ξ into Equation (19a) we get
I c C , ξ ( r m s ) ( p , n ) = 2 2 Z c , ξ | L _ 2 L 1 L _ 2 + L 1 | L l q + L g L e q + L l q + L g U o _ 1 , ξ
By substituting U o _ 2 , 2 ξ + 1 = U o _ 1 , 2 ξ + 1 into Equation (19b), the rms value of current I c C ( 0 ) is given by
I c C , 2 ξ + 1 ( r m s ) ( 0 ) = U o _ 1 , 2 ξ + 1 / ( 2 Z c , 2 ξ ( 0 ) )
where Zc,2ξ and Z c , 2 ξ ( 0 ) are the impedances for the harmonic of order 2ξ. Equations (23a) and (23c) are based on the assumption that Z c , 2 ξ   Z c , 2 ξ 1   Z c , 2 ξ ± 3 and Z c , 2 ξ ( 0 ) Z c , 2 ξ + 1 ( 0 ) . After substituting m = 2 and ωC/ωg = ξ into Equations (2) and (4a), we obtain for each pair of corresponding branches of the inverters:
U o _ 1 , 2 ξ + n = U d c π J n ( [ 1 + n / ( 2 ξ ) ] π M ) 1 + n / ( 2 ξ ) sin ( n π 2 )
The harmonic spectrum justifies limiting the scope of the order of harmonics included in Equation (23a) to the range (2ξ − 3, 2ξ + 3). Since for n = −3, −1, 1, 3, the (2ξ >> n) inequality applies, Equation (24) can be simplified to
U o _ 1 , 2 ξ + n = ( U d c / π ) J n ( π M ) sin ( n π / 2 )
After introducing Equations (22a) and (25) into Equation (23a) we get
I c C , 2 ξ ( r m s ) ( p , n ) = 2 U d c 4 π ω C L e q J 1 2 ( π M ) + 2 J 3 2 ( π M )
Similarly, we get
I c C , 2 ξ + 1 ( r m s ) ( 0 ) = 2 U d c J 1 ( π M ) / ( 4 π ω C L e q )
Because Equation (22b) can be used for m ≥ 2, the impedance Z c , ξ in Equation (23b) specified for m = 1 has been replaced by another, simplified relation (Equation (21a)) to give
I c C , ξ ( r m s ) ( p , n ) = 2 U d c π ω C L e q | L _ 2 L _ 1 L _ 2 + L _ 1 | 1 1 ( f r e s / f C ) 2 J 0 ( π M 2 )
To calculate the rms value IcC,rms for the inductance with a non-zero value tolerance, we need to find the square root:
I c C , r m s = U d c 2 π ω C L e q J 1 2 ( π M ) + J 3 2 ( π M ) + 8 [ ( % Δ L ) 100 J 0 ( π M / 2 ) 1 ( f r e s / f C ) 2 ] 2
where (%∆L) is the tolerance of the inductance of chokes L_1 and L_2. If L_1 = L_2 = L, then rms value IcC,rms can be calculated from the equation:
I c C , r m s U d c π ω C L J 1 2 ( π M ) + J 3 2 ( π M )
An equation that describes the rms of capacitor C current takes into account all the above components and takes the following form:
I c , r m s = I c G L , r m s 2 + I c C , r m s 2
The values of the functions used in Equations (29a) and (29b) are presented in Table 4.

5. The Ripple Current in Dc-Link Capacitors

Because the simulation research showed a small influence of L_1 and L_2 choke inductance asymmetry on the rms value of the total current in dc-link capacitors forced by the ripple of output currents of the inverter branch, the analysis concerning the rms value of this current will be carried out with the assumption of L_1 = L_2 = L.
The harmonics of current iv forced by the harmonics of current flowing in the ac-circuit for orders typical of PAPF (5, 7, 11, 13…) can be determined on the basis of the analysis given in [8,31].
I v , 3 n ( r m s ) = 3 2 U ^ g , 1 U d c ( λ 3 n 1 I L , 3 n 1 ( r m s ) sin φ L , 3 n 1 + λ 3 n + 1 I L , 3 n + 1 ( r m s ) sin φ L , 3 n + 1 ) 2 + ( λ 3 n 1 I L , 3 n 1 ( r m s ) cos φ L , 3 n 1 + λ 3 n + 1 I L , 3 n + 1 ( r m s ) cos φ L , 3 n + 1 ) 2
where I L , k ( r m s ) is the rms value of the compensated k-th harmonic load current, brought to the primary side of the transformer.
If the PAPF does not compensate the reactive component of the fundamental harmonic, the share of harmonics of current iv forced by the ripple currents in the inductors, L_1 and L_2 increase significantly. The current denoted in Figure 1a as i c p is half of the sum of currents of LCL filter capacitors flowing in dc-circuits.
i c p = ( i c 1 + i c 2 + i c 3 ) / 2
The equation for the current in the dc-circuit is as follows:
i d c = i v i c p
where iv is the current consumed by two inverters (iv = iv1 + iv2) and idc is the capacitor Cdc current. The following relation can be used to determine iv current:
i v = s 1 _ 1 i o 1 _ 1 + s 2 _ 1 i o 2 _ 1 + s 3 _ 1 i o 3 _ 1 + s 1 _ 2 i o 1 _ 2 + s 2 _ 2 i o 2 _ 2 + s 3 _ 2 i o 3 _ 2
where si_j denotes the switching function for the i-th branch of the j-th inverter.
Below, the harmonic coefficients of the dc-link current component ivC, forced only by the carrier and sideband harmonics of the inverter branch output currents will be determined, with the absence of fundamental and baseband harmonics in this current. Thus, the following analysis applies to harmonic currents in the inductors L_1 and L_2, for which the carrier index variable meets the assumption of m ≥ 1.

5.1. Current icp

The harmonics of voltages uoi_j of 2ξ + 1, 2ξ − 5, 4ξ − 1, 4ξ − 7, 4ξ + 5 order in the three-phase system make zero sequence components. If the relation of Equation (17) is fulfilled, m ≥ 1, ξ ≥ 40 and ∣n∣ ≤ 3 then impedance Z c , m ξ + n ( 0 ) for the mentioned harmonics is of inductive nature and increases significantly as its order increases. Thus, the order of the dominating harmonic is 2ξ + 1. The amplitude of the harmonics sum of the currents icC1,k + icC2,k + icC3,k of this order, assuming that Z c , 2 ξ + 1 ( 0 ) Z c , 2 ξ ( 0 ) , φ c , 2 ξ + 1 ( 0 ) φ c , 2 ξ ( 0 ) , is triple the amplitude I ^ c , 2 ξ + 1 . The harmonic coefficient I c , 2 ξ + 1 can be determined from Equations (18), (19b) and (25) for m = 2 and n = 1 taking into account the equality m θ C i _ 1 + n θ i = m θ C 1 _ 1 is true for i = 1, 2, 3 if m + n = 3. The time waveform of the current harmonic i c p , 2 ξ + 1 is described by
i c p , 2 ξ + 1 = I c p , 2 ξ + 1 cos [ ( 2 ξ + 1 ) ω g t + m θ C 1 _ 1 + φ c , 2 ξ ( 0 ) ]
where coefficient I c p , 2 ξ + 1 and phase angle φ c , 2 ξ ( 0 ) are described by the following relations:
I c p , 2 ξ + 1 = 3 U d c J 1 ( π M ) / ( 2 π Z c , 2 ξ ( 0 ) )
φ c , 2 ξ ( 0 ) = a r c t g { [ ω C L 1 / ( 2 ω C C ) ] / R d }
For real values of parameters ξ, L, C, and Rd in a power system with an PAPF, phase angle φ c , 2 ξ ( 0 ) π / 2 rad.

5.2. Dc-Link Capacitor Current Component idcC

The dc-link component of the capacitor current and the iv1, iv2, iv, and icp current components forced by the ripple of the output currents of inverter branches will be denoted idcC, iv1C, iv2C, ivC and icpC, respectively. Figure 7 shows the harmonics spectrum of currents iv1C, iv2C, ivC, icpC and idcC determined using the numerical model shown in Figure 1a for the modulation factor M = 0.9
Assuming sinusoidal voltages uc1, uc2, and uc3, the switching function for the i-th branch of the j-th inverter can be expressed in Fourier form by
s i _ j = 0.5 + n s = 1 S 0 n s cos ( n s ω g t + n s θ g i ) + m s = 1 n s = S m s n s cos [ ( m s ξ + n s ) ω g t + m s θ C i _ j + n s θ g i ]
where S m s n s = F m s n s / U d c .
The order of harmonic switching function is ks = msξ + ns. The output current of the i-th branch of the j-th inverter also includes the harmonics of the same order. The individual harmonics of the switching function and the currents in the three-phase system make positive, negative, or zero sequence depending on the respective algebraic sums ms ± ns or m ± n according to the equations given in Section 4.3. Figure 8 shows the harmonics spectrum of switching function si_j.
According to the considerations made in [31], certain pairs of current harmonic and switching function harmonic can force the harmonic of current iv to be of corresponding order, if certain conditions are met.
For large values of ξ (ξ ≥ 40 was assumed), harmonics coefficients of S k s can be described by the equation
S k s = { 1 / 2 for m s = 0 and n s = 0 2 ξ n s π J n s ( n s π M 2 ξ ) sin ( n s π 2 ) for m s = 0 and n s 0 2 m s π J n s ( m s π M 2 ) sin ( π [ m s + n s ] 2 ) for m s 0
For ms ≠ 0 and ξ ≥ 40 the approximation of ms + ns/ξms was adopted. Similarly for ξ ≥ 40 the approximation ( 2 ξ / π ) J 1 ( π M / 2 ξ ) M / 2 is true, which results in
S 01 M / 2
Labels S k s and S m s n s are the same ( S k s = S m s n s ) .
The source of the inductor L current ripple is the carrier and sideband harmonics forced by the inverter of the PAPF. The first component (for harmonics of any sequence) of the iok flows between the outputs of the branches of the same phase of the PAPF inverters. This component is sourced from harmonics of the difference of the output voltages of the corresponding pair of inverter branches (m = 1, 3…, n = ±2, ±4…). The second (for zero sequence harmonics) and third (for positive and negative sequence harmonics) components of the inductor ripple current contains sideband harmonics (m = 2, 4…, n = ±1, ±3…).
The time waveform of the output current harmonic of the i-th branch of the j-th inverter is given by
i o , i _ j , k = I o , k cos [ ( m ξ + n ) ω g t + m θ C i _ j + n θ g i + φ o , m ξ + n ]
For m ≠ 0 (analysis pertains to the impact of ripple currents in the inductors of equal L_1 and L_2 values on the rms value of capacitor Cdc current), harmonic coefficients of three components of currents io,i_j are given by
I o , k = 2 U d c m 2 π ω C L J n ( m π M 2 ) sin π ( m + n ) 2 for m = 1 , 3 , 5
I o , k = U d c m π Z c , m ξ ( 0 ) J n ( m π M 2 ) sin π ( m + n ) 2 for m = 2 , 4 , 6 m ± n = ± 3 , ± 6 , ± 9
I o , k = 2 U d c m π Z o , m ξ J n ( m π M 2 ) sin π ( m + n ) 2 for m = 2 , 4 , 6 m ± n = ± 1 , ± 2 , ± 4 , ± 5 , ± 7
respectively, where Z c , m ξ ( 0 ) is given by Equation (22b) for Leq = L/2, Z o , m ξ = Z o , k for k = . Equations (42a) and (42b) result from the diagram shown in Figure 5b, and the fulfilment of the conditions for their applicability means I c C , k = 0 and I o , k = I c C , k ( 0 ) /2, respectively.
Impedance Z o , m ξ (for m = 2, 4, 6…), seen from the output terminal of the VSI branch ( Z o , m ξ = | U o _ 1 , k / I o _ 1 , k | = | U o _ 2 , k / I o _ 2 , k | ) , has been determined considering the equivalent circuit shown in Figure 5a assuming L_1 = L_2 = L. Due to m = 2, 4, 6 and the equations U o _ 1 , k = U o _ 2 , k and Io_1,k = Io_2,k being met, the impedance Zo, is equal to double the value of the impedance of the circuit containing the elements between the shorted “out_1” and “out_2” outputs and the point with potential „0”.
Z o , k = 2 k ω g [ L Σ 0.5 k 2 ω g 2 C L ( L l q + L g ) ] 2 + ( k ω g C R d L Σ ) 2 [ 1 k 2 ω g 2 C ( L l q + L g ) ] 2 + ( k ω g C R d ) 2
where L = L/2 + Llq + L g . Having included the relations k = , ξ = fC/fg, and Equations (9) and (15) in Equation (43) we get
Z o , m ξ = m ω c L ( 2 L Σ L ) 2 ( 1 m 2 f C 2 / f r e s 2 ) 2 + ( 2 ξ p m f C / f r e s ) 2 ( 1 2 L Σ L m 2 f C 2 / f r e s 2 ) 2 + ( 2 ξ p m f C / f r e s ) 2
If we take the limit values 2L/L = 2, m = 2, fC/fres = 2, ζp = 0, then the value of the square root in Equation (44) deviates maximally from unity and reaches the value of 0.97. This means that on the basis of the relation in Equation (45), we can determine Z o , m ξ with an error of no more than 3%.
Z o , m ξ = m ω C L
The form of Equation (45) shows the inductive impedance character, and, therefore, φ o , m ξ ≈ −π/2 rad. Similarly the inductive character of impedances Z c , m ξ ( 0 ) and φ c , m ξ ( 0 ) ≈ −π/2 rad results from Equation (22b).
As can be seen from Equation (34), the dc-link current is the result of superposition of all six switch branches currents. Each dc-link current i v l , harmonic of the l order can be a product of many harmonic pairs of functions switching any order of ks and output currents of branches of any order of k but for which k s ± k is constant [10,31].
i v , l = k s , k i v , k s ± k ( c o m p ) for | k s ± k | = c o n s t
where i v , k s ± k ( c o m p ) is the component of the i v , l harmonic. In order to determine one harmonic current component i v , k s ± k ( c o m p ) , we can substitute into Equation (34) a harmonic of output current io,i_j,k of k order which is described by Equation (41) and a harmonic of switching function s i _ j , k s of ks order described below:
s i _ j , k s = { 1 / 2 for m s = 0 and n s = 0 S m s n s cos [ ( m s ξ + n s ) ω g t + m s θ C i _ j + n s θ g i ] for m s 0
Then we obtain the equations describing i v , k ( c o m p ) , i v , k s + k ( c o m p ) and i v , k s k ( c o m p ) in the following forms:
i v , k ( c o m p ) = 1 2 I o , k j = 1 , 2 i = 1 , 2 , 3 cos [ ( m ξ + n ) ω g t + m θ C i _ j + n θ g i + φ o , m ξ + n ]
i v , k s ± k ( c o m p ) = S k s I o , k j = 1 , 2 i = 1 , 2 , 3 { cos [ ( m s ξ + n s ) ω g t + m s θ C i _ j + n s θ g i ] cos [ ( m ξ + n ) ω g t + m θ C i _ j + n θ g i + φ o , m ξ + n ] }
where the phase angles fulfill the following relations: θC2_1 = θC1_1 ± 2π/3, θC3_1 = θC1_1 2π/3, θCi_2 = θCi_1 + π (for i = 1, 2, 3), θ g 1 = 0, θ g 2 = −2π/3, θ g 3 = 2π/3. After considering the trigonometric identity cosαcosβ = 0.5[cos(α + β) + cos(α–β)] and the equality Ci_2 = Ci_1 for even x we obtain the sum of three components in the form of
i v , k ( c o m p ) = 1 + ( 1 ) m 2 I o , k i = 1 , 2 , 3 cos [ ( m ξ + n ) ω g t + m θ C 1 _ 1 + ( m ± n ) θ g i + φ o , m ξ + n ]
i v , k s + k ( c o m p ) = 1 + ( 1 ) m s + m 2 S k s I o , k i = 1 , 2 , 3 cos { [ ( m s + m ) ξ + ( n s + n ) ] ω g t + ( m s + m ) θ C 1 _ 1 + [ ( m s ± n s ) + ( m ± n ) ] θ g i + φ o , m ξ + n }
i v , k s k ( c o m p ) = 1 + ( 1 ) m s m 2 S k s I o , k i = 1 , 2 , 3 cos { [ ( m s m ) ξ + ( n s n ) ] ω g t + ( m s m ) θ C 1 _ 1 + [ ( m s ± n s ) ( m ± n ) ] θ g i φ o , m ξ + n }
Thus, harmonic current components i v , k s ± k ( c o m p ) in combination with two inverters and with double three interleaved PWM carriers only occur when m is even (in the first component) or the sum of m + m s is even (in the second component) or the difference m s m is even (in the third component).
The interpretation of symbols “±” used in the expressions (m ± n) and ( m s ± n s ) should be that “+” is valid if the voltages u g 1 , 1 , u g 2 , 1 , u g 3 , 1 create the same sequence as the u T 1 _ 1 , 1 , u T 2 _ 1 , 1 , u T 3 _ 1 , 1 signals, while “−” is valid for the opposite sequence.
For m ≥ 1, ξ ≥ 40, and n = ±1, ±2, ±3 we assume the approximations + n and φ o , m ξ + n φ o , m ξ   φ c , 2 ξ ( 0 ) . For further considerations, the harmonic current components i v , k ( c o m p ) and i v , k s ± k ( c o m p ) will be denoted by i v C , k ( c o m p ) and i v C , k s ± k   ( c o m p ) , respectively:
i v C , k ( c o m p ) = 3 1 + ( 1 ) m 2 I o , k cos [ ( m ξ + n ) ω g t + ϕ o , m ξ ] for m s = 0 and n s = 0 m ± n = 0 , ± 3 , ± 6
where ϕ o , m ξ = m θ C 1 _ 1 + φ c , m ξ ( 0 ) .
i v C , k s + k ( c o m p ) = 3 1 + ( 1 ) m + m s 2 S k s I o , k cos { [ ( m s + m ) ξ + ( n s + n ) ] ω g t + ϕ o , m ξ } for m s 0 and n s 0 , ( m s ± n s ) + ( m ± n ) = 0 , ± 3 , ± 6
where ϕ o , m ξ = ( m s + m ) θ C 1 _ 1 + φ o , m ξ .
i v C , k s k ( c o m p ) = 3 1 + ( 1 ) m s m 2 S k s I o , k cos { [ ( m s m ) ξ + ( n s n ) ] ω g t + ϕ o , m ξ } for m s 0 and n s 0 , ( m s ± n s ) ( m ± n ) = 0 , ± 3 , ± 6
where ϕ o , m ξ = ( m s m ) θ C 1 _ 1 φ o , m ξ . The forms of the first components of Equations (48a), (49a) and (50a) formed on the basis of a constant value of the switching function (s00 = 0.5) indicate that the harmonic coefficient IvC,l is proportional to the harmonic coefficient Io,k:
I v C , l = I v C , k = 3 1 + ( 1 ) m 2 I o , k
As a product of the switching function, the harmonic coefficient of current ivC,l for ks > 0 can be determined as a sum of harmonic coefficients of individual components of this harmonic
I v C , l = 3 k s , k 1 + ( 1 ) m s + m 2 S k s I o , k for k s + k = c o n s t
or
I v C , l = 3 k s , k 1 + ( 1 ) m s m 2 S k s I o , k for | k s k | = c o n s t
The harmonic coefficient Io,2ξ + 1 resulting from Equation (42b) has the form
I o , 2 ξ + 1 = U d c 2 π J 1 ( π M ) Z c , 2 ξ ( 0 )
Based on Equations (50a) and (52), applicable for the constant component of the switching function S 00 = 0.5 , we determine the time waveform ivC,2ξ + 1:
i v C , 2 ξ + 1 = 3 U d c 2 π J 1 ( π M ) Z c , 2 ξ ( 0 ) cos ( [ 2 ξ + 1 ] ω g t + m θ C 1 _ 1 + φ c , 2 ξ ( 0 ) )
We obtain the relation describing ivC,2ξ + 1 identical to Equation (35). This means that this component does not load the capacitor Cdc. Similarly, it can be demonstrated that the harmonics of 4ξ + 5, 4ξ − 7… order also do not load the capacitor Cdc.
According to research simulation, the order of the dominant harmonic current iv forced by the ripple current in the inductor L is 2ξ – 2. In order to determine the amplitude of the ivC,2ξ − 2 harmonic, several dominant harmonic pairs of switching functions si_j and currents io,i_j were selected on the basis of Figure 8 and Table 5.
The phase angles of the selected harmonic current components i v C , k ± k s ( c o m p ) are equal or close to 2θC1_1 ± π/2 rad.
Table 6 shows the coefficients of the selected harmonic currents and switching functions, and the components of the harmonic coefficient IvC,2ξ-2 determined on the basis of Equations (39), (40), (42a), (42b), (42c), (51b) and (51c). Harmonic coefficient   I v C , 2 ξ 2 ( I v C , 2 ξ + 2 ( o s ) = I v C , 2 ξ 2 ) is described by the equation
I v C , 2 ξ 2 = 3 2 M U d c π Z o , 2 ξ [ J 3 ( π M ) J 1 ( π M ) ] + 3 2 U D C π 2 J 1 ( π M ) J 3 ( 2 π M ) ( 1 Z o , 4 ξ 1 Z o , 2 ξ ) + 8 U d c 3 π 2 ω C L { 9 J 0 ( π M 2 ) J 2 ( π M 2 ) + J 2 ( π M 2 ) J 4 ( 3 π M 2 ) + J 0 ( π M 2 ) J 2 ( 3 π M 2 ) + J 0 ( 3 π M 2 ) J 2 ( π M 2 ) }
If Equation (45) is used and both Z o , 2 ξ and Z o , 4 ξ are shown in simplified forms 2ωCL and 4ωCL, we obtain
I v C , 2 ξ 2 = U d c ω C L H 2 ξ 2
where
H 2 ξ 2 = 8 3 π 2 { J 2 ( π M 2 ) [ 9 J 0 ( π M 2 ) + J 4 ( 3 π M 2 ) + J 0 ( 3 π M 2 ) ] + J 0 ( π M 2 ) J 2 ( 3 π M 2 ) + 9 32 π M [ J 1 ( π M ) J 3 ( π M ) ] + 9 64 J 1 ( π M ) J 3 ( 2 π M ) }
The other harmonic coefficients may be determined likewise.
I v C , l = U d c ω C L H l
On the basis of simulation tests for M = 0.05, 0.225, 0.45, 0.675, and 0.9, the IvC,k harmonic spectra for k ϵ (1–20ξ) were determined and the harmonic coefficients of these orders were selected, for which the maximum rms values were not lower than 15% of the maximum rms value of the idcC current.
Table A1 (see Appendix A) contains the values of the H l ( M ) coefficients (depending on the depth of modulation M) occurring in the above equations.
The rms value of current component idcC can be written as
I d c C , r m s 1.1 2 2 ( I ^ v C , 2 ξ 2 ) 2 + ( I ^ v C , 2 ξ + 4 ) 2 + ( I ^ v C , 4 ξ 4 ) 2 + ( I ^ v C , 4 ξ + 2 ) 2 + 2 ( I ^ v C , 6 ξ 6 ) 2 + ( I ^ v C , 6 ξ ) 2 + + ( I ^ v C , 8 ξ 2 ) 2 + ( I ^ v C , 8 ξ + 4 ) 2 + ( I ^ v C , 8 ξ 10 ) 2 + ( I ^ v C , 10 ξ + 2 ) 2 + ( I ^ v C , 12 ξ ) 2 + ( I ^ v C , 18 ξ ) 2 + ( I ^ v C , 14 ξ + 4 ) 2
or in simplified form,
I d c C , r m s 0.78 U d c ω C L G ( M )
where
G ( M ) = ( H v C , 2 ξ 2 ) 2 + ( H v C , 2 ξ + 4 ) 2 + ( H v C , 4 ξ 4 ) 2 + ( H v C , 4 ξ + 2 ) 2 + 2 ( H v C , 6 ξ 6 ) 2 + ( H v C , 6 ξ ) 2 + + ( H v C , 8 ξ 2 ) 2 + ( H v C , 8 ξ + 4 ) 2 + ( H v C , 8 ξ 10 ) 2 + ( H v C , 10 ξ + 2 ) 2 + ( H v C , 12 ξ ) 2 + ( H v C , 18 ξ ) 2 + ( H v C , 14 ξ + 4 ) 2
The factor with a value of 1.1 in Equations (58) and (59), determined from simulation tests, takes into account a limited number of harmonics in Equations (58) and (60).

6. Simulation Tests of the Power Supply System

Simulation data was performed in the Turbo Pascal package. Simulation tests of the power supply system were performed with the parameters given in Table 7. The value of the L inductance (LL_1L_2) should be small enough to provide the required slew rate of the PAPF output current. On the other hand, the minimum value of L is a function of parameters such as fC, U d c , and τw, and the gain of the output current regulator of PAPF inverters [22]. The capacitance value C (twice the capacitance value of the capacitors connected to dc-link) must ensure that the condition of Equation (17) is met. In [1], for the three-phase active rectifier system, and in [2], for a system with a three-parallel power converter, the value of the damping resistor Rd = 1/(6πfresC) was adopted, from which the resonant pole damping factor ζp = 1/6 and Rd = 180 mΩ can be concluded.
Figure 9a,b shows voltage and current waveforms obtained from simulation tests for the compensation system shown in Figure 1a with the control system described in [22]. The load current, the line current, the mains voltage source and the rectified voltage ua are brought to transformer Tr1 primary side terminals. The simulation (Figure 9a,b) was performed for the same sequence u g 1 , 1 , u g 2 , 1 and u g 3 , 1 voltages and wave carriers with programmed compensation 11, 13, 23 and 25 harmonics. The amplitudes and initial phases of individual harmonic voltages and currents, which are shown by the waveforms in Figure 9a,b, are given in Table 8.
The results shown in Table 8 allow for the verification of Equations (8), (11), (13) and (29b) that determine amplitudes ÎcG,k and ÎcL,k, and rms values IcGL,rms and IcC,rms, respectively. Tests will be used to calculate the value of Ic,rms input data for simulation (L_1 = L_2 = L, C, Llq, Lg, Rd, fC, fg, υz). Amplitude values Ûg1,k and ÎL1,k and initial phases ψg1,k and ϕL1,k should be measured under real conditions. For Ûg1,k = 326.6 V and Udc = 720 V we obtain M ≈ 0.9. After placing data into Equation (29b) we obtain IcC,rms = 14.8 A (calculated values based on the harmonic amplitudes of orders 597, 599, 601 and 603, are shown in Table 8 as 13.3 A).

6.1. Rms of Capacitor C Current

Table 9 contains further results of simulation research φ g 1 , 1 is the initial phase of the fundamental harmonic current ig1). Table 10 combines the results of simulations for compensated and uncompensated (k = 1) harmonics of the load current. For compensated harmonics, the value ÎcL,k was determined based on Equation (11) and λk > 0, but for uncompensated harmonics the same equation was used with λk = 0. On the basis of Equation (13) and the data given in Table 8, IcGL1,k(rms) was calculated for k = 1, 11, 13, 23, and 25 (Table 10). Based on these harmonics, IcGL,rms = 5.9 A was obtained.
On the basis of Equations (13), (29b), and (30), a current of Ic,rms = 15.9 A was determined (the value given in Table 9 calculated in the numerical model for the course of the instantaneous capacitor current is 18.4 A). The value calculated is 15% lower than the value determined numerically for the instantaneous course of the capacitor current. To achieve good accuracy of Equation (13), which takes into account that the number of harmonics in the voltage and load current is limited to just 40, a good damping of the LCL-filter for a PAPF is required, eliminating the phenomenon of resonance.

6.2. Rms Value of Ripple Current in Dc-Link Capacitors

Based on Equation (57) and Table A1, harmonic coefficients IvC,2ξ-2, IvC,2ξ + 4, IvC,4ξ − 4, IvC,4ξ + 2, IvC,6ξ − 6, IvC,6ξ, and IvC,6ξ + 6 can be determined for different coefficients of modulation depth M. The values of the amplitudes calculated in this way and the amplitudes determined in the numerical model are presented in Table 11. The main source of error in the calculations given in Table 11 is that the H l ( M ) harmonic factors describing these harmonics were determined on the basis of a limited number of harmonic pairs of the switching function and the output currents of the inverter branches. The second source of error is the lack of consideration of harmonic components forced by the active current taken from the grid and related to thermal losses in the inverters.
Table 12 contains the results of simulation tests carried out for the control system of the PAPF output current with asymmetrical and symmetrical regular-sampled PWM.
The IdcC,rms value determined from Equation (59) for fC = 15 kHz and M = 0.9 is 35.2 A (the value given in Table 12 calculated in the numerical model for the course of the instantaneous capacitor Cdc current is 36.2 A).

6.3. The Currents of LCL-Ripple Filter and Dc-Link Capacitors in a PAPF with One and Two Voltage Source Inverters for Positive and Negative Sequence Carriers

Simulation tests were carried out with inductors L_1 and L_2 with identical inductance values equal to L in a system with two inverters and inductors with inductance values equal to half of L in a system with one inverter. This choice of inductance values means the same resonance frequency of the LCL filter in systems with two and one inverters.
Simulation research shows that the harmonic spectrum of currents idcC changes with the change of carrier sequence, but the rms current value of the Cdc capacitor does not change (Table 12). This is due to the equality of the relevant harmonic coefficients: I v C , 2 ξ + 2 ( o s ) = I v C , 2 ξ 2 , I v C , 4 ξ + 4 ( o s ) = I v C , 4 ξ 4 , I v C , 2 ξ 4 ( o s ) = I v C , 2 ξ + 4 , I v C , 4 ξ 2 ( o s ) = I v C , 4 ξ + 2 , I v C , 6 ξ + 6 ( o s ) = I v C , 6 ξ 6 , I v C , 6 ξ 6 ( o s ) = I v C , 6 ξ + 6 , I v C , 6 ξ ( o s ) = I v C , 6 ξ . I v C , k ( o s ) is the harmonic coefficient of current idc forced in a system with the opposite sequence of the carrier.
Table 13 contains results of PAPF simulation tests without load for fC = 10 kHz with inductance L_1 and L_2 of equal or different value and symmetrical or asymmetrical supply voltages. The research shows that if the values of L_1 and L_2 and the values of the phase voltages of the grid are held to within ±5%, the rms value of current idcC does not exceed 106% of the value calculated for the system with inductances L_1 and L_2 of equal value and symmetrical three-phase supply voltage.
Simulation tests, the results of which are presented in Table 11, Table 12 and Table 13, were performed for sinusoidal (not containing higher harmonic) voltages u g 1 , u g 2 and u g 3 .

7. Experimental Results

In order to verify the correctness of the analysis of components of the currents of the LCL-ripple filter capacitor and inductors, studies were conducted with the FAS-400k-400 and FA3-100k-400 active filters produced by the MEDCOM company based in Warsaw Jutrzenki 78A str., Poland. Tests with FAS-400k-400 were conducted in industrial conditions, in a coal mine where the PAPF works as a compensator of high harmonic currents consumed by a lifting machine in the mine with preset compensation of the 11th, 13th and 23rd harmonics. Filter parameters U d c , Cdc, C, L, fC, and tdead correspond with those assumed for the simulation research. Parameters of Tr1 are: transformer turns ratio υz is 26; rated low voltage is 400 V, rated high voltage is 6 kV, rated power is 1.6 MVA, and short-circuit voltage is 5.68%.
Figure 10a–c shows the research results of the PAPF for fC = 15 kHz, while Figure 11a–c shows the results for fC = 10 kHz. The rms value of the LCL capacitor current shown in Figure 10c is 17.9 A, while this current shown in Figure 11a is 28 A. Figure 10a shows the oscilloscope waveform of output current io1_1 (CH2). Figure 9b shows the waveform of current io1_1 obtained in the numerical model. The similarity of these waveforms (Δio1_1,max ≈ 100A) proves the high precision of the numerical model and thus also the parameters listed in Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13. The rms value of the LCL capacitor current (16.2 A) determined from the equations derived in the article, taking into account the 5% tolerance of inductance of converter-side inductors, is 10% less than in the real system (17.9 A).
Figure 11b shows the current idc harmonic spectrum for a negative sequence carrier and fC = 10 kHz in an PAPF-type FAS-400k-400 with inductances L_1 = L_2 = L = 80 µH, while Figure 11c shows current idc harmonic spectrum for a positive sequence carrier and fC = 10kHz in an PAPF-type FA3-100k-400 with inductances L_1 = L_2 = L = 200 µH.
On the basis of the oscillogram shown in Figure 11b, the current was obtained as I v C , 2 ξ + 2 ( o s ) = 35.2 A. The value calculated by Equation (57) for the parameters of this filter ( U d c = 720 V, fC = 10 kHz and L = 80 μH) and the coefficients given in Table A1 for M = 0.9 is 36.2 A. Based on the oscillogram shown in Figure 11c, the following results were obtained for FA3-100k-400: I v C , 2 ξ 2 ( r m s ) = 14.7 A; I v C , 2 ξ + 4 ( r m s ) = 7.98 A. The calculated values (for U d c = 720 V, fC = 10 kHz and L = 200 µH) are 14.5 A and 7.82 A, respectively. The oscillograms shown in Figure 11b,c were made in the PAPF system without load. Figure 12 shows a photo of the current model of the FA3-400-400 filter.
The PAPF reduces the THDi factor of the consumed current of the 6 kV grid from 14.7% to 4.8%.
It is assumed that all elements of the examined energy system are normal or functional, which is not in line with the reality [33]. PAPF is an element of the system, the turn off of which from work results only in the deterioration of the quality of the voltage supplying the energy system. Switching the APF off in the system is safe because power of the device is 10 times lower than the power of the compensated load.
PAPF has systems for measuring the capacitance of dc-link circuit capacitors and LCL-filter capacitors when switching the device on and off, and continuous control of the rms value of LCL filter capacitor currents. The symmetry of current distribution of the respective branches of both inverters is also controlled. Exceeding any PAPF limit parameter results in switching off the device from the power supply system.

8. Conclusions

Theoretical analysis and simulation research confirm a significant participation in PAPF systems with two parallel interleaved inverters of the LCL filter capacitor current component forced by non-linear load current and harmonics in the output current of a PAPF that compensates harmonics in the load. Simulation research supported by examples of calculations and experimental research confirm the validity of formulae describing different components of the capacitor current of the LCL-ripple filter, allowing for election of the capacitor from catalogues from the point of view of its effective current limit value. The analysis and relations presented allow calculation of the rms value of the LCL filter capacitor current.
If the PAPF does not compensate the reactive component of the fundamental harmonic, the share of harmonics current in the dc-link capacitor forced by the ripple currents in the inductors L_1 and L_2 increases significantly. In PAPF systems with two symmetrical inverters with output inductors L and double three interleaved carriers, the dc-link capacitor current does not contain carrier harmonics of odd orders or their sideband harmonics.
Simulation studies show that if the values of inductance L_1 and L_2, and the phase voltages of the mains, are held within ±5%, the rms value of the current in the dc-link forced by the ripple current in the branches of the inverters does not exceed 106% of the value calculated for the symmetrical system.
The rms values of ripple currents in the LCL capacitors, LCL inductors and dc-link capacitors, which are sourced from the ripple current in the L_1 and L_2 inductors, significantly depend on transistors’ switching frequencies and on the value of the L_1 and L_2 inductances; however, they are not significantly dependent on the level of compensated power.
Supported by simulations, the analysis shows that if the sequence of three triangular carrier waves in the individual branches of the inverter is different from the harmonic sequence of the basic output voltages of the individual branches of the inverter, then the distribution of the spectrum of the harmonic currents in the dc-link capacitor forced by ripples of the output current of the inverter branches is changed; but this does not change the total effective value of this current.
The formulae describing the rms current of the LCL filter capacitor and the dc-link capacitor current given in the article may be applied to other topologies, such as a three-phase grid-connected converter (excluding the formula for the second mentioned component of the current of the capacitor of the LCL-filter), assuming that the fC/fg ≥ 40 calculation error will not change.
On the basis of the simulation tests (only some results are given in the article), it can be concluded that for fC/fg ≥ 40:
  • the rms values of currents in LCL filter capacitors in PAPF systems with double interleaved PWM carriers and double three interleaved PWM are similar, with accuracy of 1.5%; and
  • the rms values of total currents in LCL filter capacitors and dc-link capacitors in PAPF systems with a current control system with asymmetrical and symmetrical regular-sampled PWM are similar, with accuracy of 6%.
It seems advisable to carry out a similar analysis for a PAPF containing one voltage inverter.

Funding

The author gratefully acknowledges the financial support provided by Institute of Control and Industrial Electronics (WUT).

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

F m n , F k I m n , I k , U o 1 , k , U o 2 , k Fourier harmonic coefficients
F ^ m n , F ^ k , I ^ m n , I ^ k , U ^ o 1 , k , U ^ o 2 , k ,Amplitudes of harmonic coefficients
U ^ g , k , I ^ L , k Amplitudes of voltage and current harmonics brought to terminals 1, 2, 3 of Tr1 transformer primary side
I m n , r m s ,   I m n ( r m s ) , I k , r m s ,   I k ( r m s ) rms value of current harmonics
Z c , k , Z c , k ( 0 ) The circuit impedances for harmonics of positive and negative sequence and of zero sequence, respectively
i v , k s ± k ( c o m p ) The harmonic current component of the ks ± k order
φ o , 2 ξ ,   φ c , 2 ξ ( 0 ) The phase angles for harmonics of positive and negative sequence and of zero sequence, respectively
NThe number of operations in the microprocessor interruption service routine
λHarmonic attenuation factor
υzTransformer turns ratio

Appendix A

Table A1. The values of H l ( M ) and G(M) coefficients.
Table A1. The values of H l ( M ) and G(M) coefficients.
MH2ξ − 2H2ξ + 4H4ξ − 4H4ξ + 2H6ξ±6H6ξH8ξ − 8H8ξ − 2H8ξ + 4H8ξ + 10H10ξ + 2H12ξH14ξ + 4H18ξG(M)
0.000.0000.0000.0000.0000.0000.4860.0000.0000.0000.0000.000−0.2460.0000.1580.567
0.050.0050.0000.0000.0090.0000.4580.0000.018−0.0000.000−0.022−0.1950.0010.0920.507
0.100.0190.0000.0000.0360.0000.3810.0000.063−0.0020.000−0.073−0.0740.010−0.0230.403
0.150.0420.0000.0020.0760.0000.2680.0000.118−0.0110.000−0.1220.0400.036−0.0540.338
0.200.0730.0010.0050.1260.0010.1400.0000.158−0.0280.000−0.1370.0850.0680.0020.312
0.250.1100.0020.0120.1780.0020.0180.0000.167−0.0560.000−0.1050.0550.0770.0380.311
0.300.1510.0050.0230.2270.005−0.0790.0010.139−0.0880.000−0.039−0.0080.0470.0090.335
0.350.1950.0080.0390.2670.011−0.1400.0020.082−0.1150.0000.030−0.051−0.006−0.0270.393
0.400.2380.0140.0600.2940.020−0.1590.0050.014−0.1280.0010.072−0.047−0.042−0.0140.446
0.450.2800.0220.0840.3040.034−0.1430.005−0.044−0.1210.0030.071−0.011−0.0360.0170.474
0.500.3170.0330.1120.2970.051−0.1020.017−0.076−0.0920.0060.0350.025−0.0020.0160.485
0.550.3490.0460.1400.2750.071−0.0500.027−0.078−0.0490.011−0.0090.0350.026−0.0090.493
0.600.3740.0630.1660.2400.0910.0000.038−0.055−0.0010.019−0.0380.0190.025−0.0150.504
0.650.3920.0810.1890.1970.1100.0370.050−0.0210.0380.030−0.042−0.0060.0030.0010.517
0.700.4010.1020.2060.1500.1230.0570.0500.0110.0600.043−0.025−0.019−0.0160.0120.529
0.750.4010.1250.2150.1050.1270.0600.0640.0280.0610.056−0.001−0.015−0.0150.0040.531
0.800.3940.1480.2170.0650.1210.0500.0630.0300.0460.0660.017−0.003−0.002−0.0060.521
0.850.3790.1710.2090.0340.1050.0340.0570.0200.0230.0720.0200.0050.007−0.0060.501
0.900.3580.1930.1940.0130.0830.0190.0460.006−0.0010.0720.0130.0050.0060.0010.475
0.950.3320.2140.1710.0040.0590.0120.046−0.002−0.0160.0650.0040.0020.0010.0040.445
1.000.3010.2320.1430.0070.0380.0140.018−0.001−0.0190.0530.0010.0000.0010.0020.414

References

  1. Liserre, M.; Blaabjerg, F.; Hansen, S. Design and Control of an LCL-Filter-Based Three-Phase Active Rectifier. IEEE Trans. Ind. Appl. 2005, 41, 1281–1291. [Google Scholar] [CrossRef]
  2. Jeong, H.G.; Yoon, D.-K.; Lee, K.B. Design of an LCL-Filter for Three-Parallel Operation of Power Converters in Wind Turbines. J. Power Electron. 2013, 13, 437–445. [Google Scholar] [CrossRef]
  3. Asiminoaei, L.; Aeloiza, E.; Enjeti, P.; Blaabjerg, F. Shunt Active-Power_Filter Topology Based on Parallel Interleaved Inverters. IEEE Trans. Ind. Electron. 2008, 55, 1175–1189. [Google Scholar] [CrossRef]
  4. Jalili, K.; Bernet, S. Design of LCL Filters of Active-Front-End Two-Level Voltage-Source Converters. IEEE Trans. Ind. Electron. 2009, 56, 1674–1689. [Google Scholar] [CrossRef]
  5. Sahoo, A.K.; Shahani, A.; Basu, K.; Mohan, N. LCL filter design for grid-connected inverters by analytical estimation of PWM ripple voltage. In Proceedings of the 2014 Twenty-Ninth Annual Applied Power Electronics Conference and Exposition (APEC), Fort Worth, TX, USA, 16–20 March 2014; pp. 1281–1286. [Google Scholar]
  6. Tang, Y.; Loh, P.C.; Wang, P.; Choo, F.H.; Gao, F.; Blaabjerg, F. Generalized design of high performance shunt active power filter with output LCL filter. IEEE Trans. Ind. Electron. 2012, 59, 1443–1452. [Google Scholar] [CrossRef]
  7. Vodyakho, O.; Mi, C.C. Three-level inverter-based shunt active power filter in three-phase three-wire and four-wire system. IEEE Trans. Power Electron. 2009, 24, 1350–1363. [Google Scholar] [CrossRef]
  8. McGrath, B.P.; Holmes, D.G. A general analytical method for calculating inverter DC-link current harmonics. IEEE Trans. Ind. Appl. 2009, 45, 1851–1859. [Google Scholar] [CrossRef]
  9. Kolar, J.W.; Wolbank, T.M.; Schrödl, M. Analytical calculation of the RMS current stress on the DC link capacitor of voltage DC link PWM converter systems. IEE Proc. Electr. Power Appl. 2006, 153, 535–543. [Google Scholar] [CrossRef]
  10. Zhang, H.; Wheeler, N.; Grant, D. Switching harmonics in the DC link current in a PWM AC-DC-AC converter. In Proceedings of the IAS ′95. Conference Record of the 1995 IEEE Industry Applications Conference Thirtieth IAS Annual Meeting, Orlando, FL, USA, 8–12 October 1995; pp. 2649–2655. [Google Scholar]
  11. Zhang, D.; Wang, F.; Burgos, R.; Lai, R.; Boroyevich, D. DC-link ripple current reduction for paralleled three-phase voltage-source converters with interleaving. IEEE Trans. Power Electron. 2011, 26, 1741–1753. [Google Scholar] [CrossRef]
  12. Ye, H.; Emadi, A. An interleaving scheme to reduce DC-link current harmonics of dual traction inverters in hybrid electric vehicles. In Proceedings of the 2014 IEEE Applied Power Electronics Conference and Exposition-APEC 2014, Fort Worth, TX, USA, 16–20 March 2014; pp. 3205–3211. [Google Scholar]
  13. Quan, Z.; Li, Y. Impact of Carrier Phase Shift PWM on the DC Link Current of Single and Interleaved Three-Phase Voltage Source Converters. In Proceedings of the 2017 IEEE Energy Conversion Congress and Exposition (ECCE), Cincinnati, OH, USA, 1–5 October 2017; pp. 3851–3855. [Google Scholar]
  14. Smolenski, R.; Jasinski, M.; Jarnut, M.; Bojarski, J.; Cecati, C. Compensation of CM voltage in systems consisting of interleaved AC-DC converters. In Proceedings of the 2014 IEEE 23rd International Symposium on Industrial Electronics (ISIE), Istanbul, Turkey, 1–4 June 2014; pp. 1996–2001. [Google Scholar]
  15. Malinowski, M.; Bernet, S. A Simple Voltage Sensorless Active Damping Scheme for Three Phase PWM Converters with an LCL Filter. IEEE Trans. Ind. Electron. 2008, 55, 1876–1880. [Google Scholar] [CrossRef]
  16. Wu, D.; Chen, Y.; Hong, S.; Zhao, X.; Luo, J.; Gu, Z. Mathematical model analysis and LCL Filter Design of VSC. In Proceedings of the 7th International Power Electronics and Motion Control Conference, Harbin, China, 2–5 June 2012; pp. 2799–2804. [Google Scholar]
  17. Lang, Y.; Xu, D.; Hadianamrei, S.R.; Ma, H. A novel design method of LCL type utility interface for three-phase voltage source rectifier. In Proceedings of the 2005 IEEE 36th Power Electronics Specialists Conference, Recife, Brazil, 16 June 2005; pp. 313–317. [Google Scholar]
  18. Shen, G.; Zhu, X.; Zhang, J.; Xu, D. A New Feedback Method for PR Current Control of LCL-Filter-Based Grid-Connected Inverter. IEEE Trans. Ind. Electron. 2010, 57, 2033–2041. [Google Scholar] [CrossRef]
  19. Pei, X.; Zhou, W.; Kang, Y. Analysis and calculation of DC-link current and voltage ripples for three-phase inverter with unbalanced load. IEEE Trans. Power Electron. 2015, 30, 5401–5412. [Google Scholar] [CrossRef]
  20. Lascu, C.; Asiminoaei, L.; Boldea, L.; Blaabjerg, F. Frequency response analysis of current controllers for selective harmonic compensation in active power filters. IEEE Trans. Ind. Electron. 2009, 56, 337–347. [Google Scholar] [CrossRef]
  21. Ben-Sheng, C.; Yuan-Yih, H. A minimal harmonic controller for a STATCOM. IEEE Trans. Ind. Electron. 2008, 55, 655–664. [Google Scholar]
  22. Platek, T.; Osypinski, T. Current Control with Asymmetrical Regular Sampled PWM Modulator Applied for Parallel Active Filter. Bull. Pol. Acad. Sci. Tech. 2016, 64, 287–300. [Google Scholar]
  23. Holmes, D.G.; Lipo, T.A. Pulse Width Modulation For Power Converters; IEEE PRESS: Hoboken NJ 07030 USA, Willey-Interscience; 2003. [Google Scholar]
  24. Zhang, D.; Wang, F.; Burgos, R.; Lai, R.; Boroyevich, D. Impact of interleaving on ac passive components of paralleled three-phase voltage-source converters. IEEE Trans. Ind. Appl. 2010, 46, 1042–1054. [Google Scholar] [CrossRef]
  25. Rockhill, A.A.; Liserre, M.; Teodorescu, R.; Rodriguez, P. Grid-Filter Design for a Multimegawatt Medium-Voltage Voltage-Source Inverter. IEEE Trans. Ind. Electron. 2011, 58, 1205–1217. [Google Scholar] [CrossRef]
  26. Shang, Y. Synchronization in networks of coupled harmonic oscillators with stochastic perturbation and time delays. Ann. Acad. Rom. Sci. Ser. Math. Appl. 2012, 4, 44–58. [Google Scholar]
  27. Holmes, D.G.; Lipo, T.A.; McGrath, B.P.; Kong, W.Y. Optimized design of stationary frame three phase AC current regulators. IEEE Trans. Power Electron. 2009, 24, 2417–2426. [Google Scholar] [CrossRef]
  28. Liu, F.; Zha, X.; Zhou, Y.; Duan, S. Design and research on parameter of LCL filter in three-phase grid-connected inverter. In Proceedings of the 2009 IEEE 6th International Power Electronics and Motion Control Conference, Wuhan, China, 17–20 May 2009; pp. 2174–2177. [Google Scholar]
  29. Liu, Y.; Lai, C.M. LCL Filter Design with EMI Noise Consideration for Grid-Connected Inverter. Energies 2018, 11, 1646. [Google Scholar] [CrossRef]
  30. Abusara, M.A.; Jamil, M.; Sharkh, S.M. Repetitive current control of an interleaved grid connected inverter. In Proceedings of the 2012 3rd IEEE International Symposium on Power Electronics for Distributed Generation Systems (PEDG), Aalborg, Denmark, 25-28 June 2012; pp. 558–563. [Google Scholar]
  31. Jiang, Y.; Ekström, A. General Analysis of Harmonic Transfer Through Converters. IEEE Trans. Power Electron. 1997, 12, 287–293. [Google Scholar] [CrossRef]
  32. Kimball, J.W.; Zawodniok, M. Reducing Common-Mode Voltage in Three-Phase Sine-Triangle PWM with Interleaved Carriers. In Proceedings of the 2010 Twenty-Fifth Annual IEEE Applied Power Electronics Conference and Exposition (APEC), Palm Springs, CA, USA, 21–25 February 2010; Volume 26, pp. 2229–2236. [Google Scholar]
  33. Shang, Y. Resilient Multiscale Coordination Control against Adversarial Nodes. Energies 2018, 11, 1844. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Diagram of circuit with PAPF and (b) its simplified version.
Figure 1. (a) Diagram of circuit with PAPF and (b) its simplified version.
Energies 12 04493 g001
Figure 2. Harmonics spectrum of output voltages of each branch of inverters (M = 0.9).
Figure 2. Harmonics spectrum of output voltages of each branch of inverters (M = 0.9).
Energies 12 04493 g002
Figure 3. Equivalent diagram of LCL circuit in a filter system with two inverters: (a) for the current component forced by the ac mains and (b) for the current component forced by load.
Figure 3. Equivalent diagram of LCL circuit in a filter system with two inverters: (a) for the current component forced by the ac mains and (b) for the current component forced by load.
Energies 12 04493 g003
Figure 4. Carrier waveforms for double PWM techniques. (a) Double Interleaved PWM Carriers [3] and (b) Double Three Interleaved PWM Carriers [13,22].
Figure 4. Carrier waveforms for double PWM techniques. (a) Double Interleaved PWM Carriers [3] and (b) Double Three Interleaved PWM Carriers [13,22].
Energies 12 04493 g004
Figure 5. (a) Equivalent diagram of LCL circuit in an PAPF with two inverters for the current component icC,k for the described control for positive and negative sequence harmonics and (b) for zero sequence harmonics.
Figure 5. (a) Equivalent diagram of LCL circuit in an PAPF with two inverters for the current component icC,k for the described control for positive and negative sequence harmonics and (b) for zero sequence harmonics.
Energies 12 04493 g005
Figure 6. Harmonics spectrum of current component icC (M = 0.9, L_1 = 76 μH, L_2 = 84 μH), where (n), (p), (0) denote negative, positive and zero sequence harmonics, respectively (for the same sequence of u g 1 , 1 , u g 2 , 1 , u g 3 , 1 voltages and wave carriers).
Figure 6. Harmonics spectrum of current component icC (M = 0.9, L_1 = 76 μH, L_2 = 84 μH), where (n), (p), (0) denote negative, positive and zero sequence harmonics, respectively (for the same sequence of u g 1 , 1 , u g 2 , 1 , u g 3 , 1 voltages and wave carriers).
Energies 12 04493 g006
Figure 7. Harmonics spectrum of currents iv1C, iv2C, ivC, icpC and idcC in the system without load for the same sequence u g 1 , 1 , u g 2 , 1 and u g 3 , 1 voltages and wave carriers (M = 0.9).
Figure 7. Harmonics spectrum of currents iv1C, iv2C, ivC, icpC and idcC in the system without load for the same sequence u g 1 , 1 , u g 2 , 1 and u g 3 , 1 voltages and wave carriers (M = 0.9).
Energies 12 04493 g007
Figure 8. Harmonics spectrum of switching function si_j for the same sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers (M = 0.9).
Figure 8. Harmonics spectrum of switching function si_j for the same sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers (M = 0.9).
Energies 12 04493 g008
Figure 9. (a) Simulation of waveforms of currents: iL1, i1, ig1, and ic1 and voltages ua and u12; and (b) simulation of the inverter output currents io1_1 and io1_2 and their sum i1, output current of the PAPF ip1, load current iL1, and filter capacitor current ic1.
Figure 9. (a) Simulation of waveforms of currents: iL1, i1, ig1, and ic1 and voltages ua and u12; and (b) simulation of the inverter output currents io1_1 and io1_2 and their sum i1, output current of the PAPF ip1, load current iL1, and filter capacitor current ic1.
Energies 12 04493 g009
Figure 10. (a) Waveforms of load current iL1 (CH1:45 A/div), output current io1_1 (CH2:200 A/div), setpoint signal of current i1 (CH3:1 V/div) and setpoint signal of current component i1 to regulate the dc-link voltage udc (CH4:2 V/div); (b) waveforms of load current iL1 (CH1:90 A/div) and input mains current ig1 (CH2:120 A/div); (c) waveforms of output current io1_1 (CH1:500 A/div), capacitor current ic1 (CH2:500 A/div) and PAPF output current ip1 (CH3:500 A/div).
Figure 10. (a) Waveforms of load current iL1 (CH1:45 A/div), output current io1_1 (CH2:200 A/div), setpoint signal of current i1 (CH3:1 V/div) and setpoint signal of current component i1 to regulate the dc-link voltage udc (CH4:2 V/div); (b) waveforms of load current iL1 (CH1:90 A/div) and input mains current ig1 (CH2:120 A/div); (c) waveforms of output current io1_1 (CH1:500 A/div), capacitor current ic1 (CH2:500 A/div) and PAPF output current ip1 (CH3:500 A/div).
Energies 12 04493 g010
Figure 11. (a) Waveforms of output current io1_1 (CH1:500A/div), capacitor current ic1 (CH2:500A/div) and PAPF output current ip1 (CH3:500A/div); (b) harmonics spectrum of current idc (8Arms/div) for the opposite sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers; (c) harmonics spectrum of current idc (4Arms/div) for the same sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers (in an PAPF-type FA3-100k-400).
Figure 11. (a) Waveforms of output current io1_1 (CH1:500A/div), capacitor current ic1 (CH2:500A/div) and PAPF output current ip1 (CH3:500A/div); (b) harmonics spectrum of current idc (8Arms/div) for the opposite sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers; (c) harmonics spectrum of current idc (4Arms/div) for the same sequence u g 1 , 1 , u g 2 , 1 , and u g 3 , 1 voltages and wave carriers (in an PAPF-type FA3-100k-400).
Energies 12 04493 g011
Figure 12. FA3-400k-400 installed in a coal mine.
Figure 12. FA3-400k-400 installed in a coal mine.
Energies 12 04493 g012
Table 1. Relation of primary and secondary vectors of Dy5 transformer voltage.
Table 1. Relation of primary and secondary vectors of Dy5 transformer voltage.
Relation   for   u g Relation   for   u Relation   for   i g Relation   for   i L Relation   for   i s Matrix A
u g = 1 υ z A u g u = 1 υ z A u i g = υ z 3 A i g i L = υ z 3 A i L i s = υ z 3 A i s A = [ 1 1 0 0 1 1 1 0 1 ]
Table 2. The relation of the primary and the secondary values of Dy5 transformer voltage.
Table 2. The relation of the primary and the secondary values of Dy5 transformer voltage.
Relation   for   U ^ g , k Relation   for   ψ g , k Relation   for   I ^ L , k Relation   for   φ L , k
U ^ g , k = 3 U ^ g , k υ z ψ g , k = ψ g , k + 5 π 6 I ^ L , k = υ z I ^ L , k 3 φ L , k = φ L , k + 5 π 6
Table 3. Conditions for the type of three-phase harmonic sequence of inverter output voltages.
Table 3. Conditions for the type of three-phase harmonic sequence of inverter output voltages.
( u g 1 , 1 ,   u g 2 , 1 ,   u g 3 , 1   )   and   ( u T 1 _ ,   2 _ ,   3 _   )   The Sequence of the Individual Harmonic Voltages of the Three-Phase Outputs uo1_1, uo2_1, uo3_1 and uo1_2, uo2_2, uo3_2
-positive sequence (p)negative sequence (n) zero sequence (0)
the same sequencem + n = −5, −2,1,4,7…m + n = −4, −1, 2, 5…m + n = −6, −3, 0, 3, 6…
opposite sequencemn = –4, −1,2,5…mn = −5, −2, 1, 4, 7…mn = −6, −3, 0, 3, 6…
Table 4. The values of functions used in formulae.
Table 4. The values of functions used in formulae.
FunctionsValues
M0.000.10.20.30.40.50.60.70.80.91
0.050.150.250.350.450.550.650.750.850.95
J0(πM/2)1.0000.9940.9750.9450.9040.8520.7900.7200.6430.5590.472
0.9980.9860.9620.9260.8790.8220.7560.6820.6020.516
J 1 2 ( π M ) + J 3 2 ( π M ) 0.0000.1550.2990.4210.5140.5710.5920.5790.5400.4870.438
0.0780.2290.3630.4710.5470.5860.5890.5620.5150.461
Table 5. Conditions for harmonic components in dc-link of two parallel interleaved inverters.
Table 5. Conditions for harmonic components in dc-link of two parallel interleaved inverters.
Input, Output ParametersExisting Conditions for Current Harmonics iv,kNo Harmonics
Sequence of ioi_j,k(p)(n)(0)(p)(n)(0)(0)(p), (n)
Sequence of si_j,ks(p)(n)(0)(n)(p)(0)(p), (n)(0)
Conditionsms ± m = −6, −4, −2, 0, 2, 4, 6…Do not exist
(ms ± ns) − (m ± n) = 0, ±3, ±6…(ms ± m) + (ns ± n) = 0, ±3, ±6…
Order of harm. of i v C , k ± k s ( c o m p ) abs(kks)k + ks
Table 6. Coefficients of the i v C , 2 ξ 2 ( c o m p ) and i v C , 2 ξ + 2 ( c o m p ) components for the same (ss) and opposite (os) sequences of u g i , 1 voltages and wave carriers, respectively.
Table 6. Coefficients of the i v C , 2 ξ 2 ( c o m p ) and i v C , 2 ξ + 2 ( c o m p ) components for the same (ss) and opposite (os) sequences of u g i , 1 voltages and wave carriers, respectively.
S k s I o , k Coefficients   of   i v C , 2 ξ 2 ( c o m p )   and   i v C , 2 ξ + 2 ( c o m p )   Components ϕ v C , 2 ξ 2 = ϕ v C , 2 ξ + 2
ms = 1, ns = 0 (p)m = 1, n = −2 (n)ms + m = 2, ns + n = −2 (ss)
ms = 1, ns = 0 (n)m = 1, n = 2 (p)ms + m = 2, ns + n = 2 (os)
( 2 / π ) J 0 ( π M / 2 ) 2 U d c J 2 ( π M / 2 ) / ( ω C L π ) 12 U d c J 2 ( π M / 2 ) J 0 ( π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 1, ns = −2(n)m = 1, n = 0 (p)ms + m = 2, ns + n = −2 (ss)
ms = 1, ns = 2 (p)m = 1, n = 0 (n)ms + m = 2, ns + n = −2 (os)
( 2 / π ) J 2 ( π M / 2 ) 2 U d c J 0 ( π M / 2 ) / ( ω C L π ) 12 U d c J 2 ( π M / 2 ) J 0 ( π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 0, ns = 1 (p)m = 2, n = −1 (p)mms = 2, nns = −2 (ss)
ms = 0, ns = 1 (p)m = 2, n = 1 (n)m + ms = 2, n + ns = 2 (os)
M / 2 U d c J 1 ( π M ) / ( π Z o , 2 ξ ) 3 M U d c J 1 ( π M ) / ( 2 π Z o , 2 ξ ) 2 θ C 1 _ 1 π / 2
ms = 0, ns = 1 (p)m = 2, n = −3 (n)m + ms = 2, n + ns = −2 (ss)
ms = 0, ns = 1 (p)m = 2, n = 3 (p)mms = 2, nns = 2 (os)
M / 2 U d c J 3 ( π M ) / ( π Z o , 2 ξ ) 3 M U d c J 3 ( π M ) / ( 2 π Z o , 2 ξ ) 2 θ C 1 _ 1 π / 2
ms = 1, ns = 0 (p)m = 3, n = −2 (p)mms = 2, nns = −2 (ss)
ms = 1, ns = 0 (n)m = 3, n = 2 (n)mms = 2, nns = 2 (os)
( 2 / π ) J 0 ( π M / 2 ) 2 U d c J 2 ( 3 π M / 2 ) / ( 9 π ω C L ) 4 U d c J 0 ( π M / 2 ) J 2 ( 3 π M / 2 ) / ( 3 ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 3, ns = −4(n)m = 1, n = −2 (n)msm = 2, nsn = −2 (ss)
ms = 3, ns = 4 (p)m = 1, n = 2 (p)msm = 2, nsn = 2 (os)
2 J 4 ( 3 π M / 2 ) / ( 3 π ) 2 U d c J 2 ( π M / 2 ) / ( π ω C L ) 4 U d c J 2 ( π M / 2 ) J 4 ( 3 π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 + π / 2
4 U d c J 2 ( π M / 2 ) J 4 ( 3 π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 1, ns = −(n)m = 3, n = −4 (n)mms = 2, nns = −2 (ss)
ms = 1, ns = 2 (p)m = 3, n = 4 (p)mms = 2, nns = 2 (os)
( 2 / π ) J 2 ( π M / 2 ) 2 U d c J 4 ( 3 π M / 2 ) / ( 9 π ω C L ) 4 U d c J 2 ( π M / 2 ) J 4 ( 3 π M / 2 ) / ( 3 ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 3, ns = −2(p)m = 1, n = 0 (p)msm = 2, nsn = −2 (ss)
ms = 3, ns = 2 (n)m = 1, n = 0 (n)msm = 2, nsn = 2 (os)
2 J 2 ( 3 π M / 2 ) / ( 3 π ) 2 U d c J 0 ( π M / 2 ) / ( π ω C L ) 4 U d c J 0 ( π M / 2 ) J 2 ( 3 π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 + π / 2
ms = 3, ns = 0 (0)m = 1, n = 2 (0)msm = 2, nsn = −2; (ss)
ms = 3, ns = 0 (0)m = 1, n = −2 (0)msm = 2, nsn = 2; (os)
2 J 0 ( 3 π M / 2 ) / ( 3 π ) 2 U d c J 2 ( π M / 2 ) / ( π ω C L ) 4 U d c J 0 ( 3 π M / 2 ) J 2 ( π M / 2 ) / ( ω C L π 2 ) 2 θ C 1 _ 1 + π / 2
ms = 1, ns = 2 (0)m = 3, n = 0 (0)mms = 2, nns = −2 (ss)
ms = 1, ns = −2 (0)m = 3, n = 0 (0)mms = 2, nns = 2 (os)
( 2 / π ) J 2 ( π M / 2 ) 2 U d c J 0 ( 3 π M / 2 ) / ( 9 π ω C L ) 4 U d c J 0 ( 3 π M / 2 ) J 2 ( π M / 2 ) / ( 3 ω C L π 2 ) 2 θ C 1 _ 1 π / 2
ms = 2, ns = −1(p)m = 4, n = −3 (p)mms = 2, nns = −2 (ss)
ms = 2, ns = 1 (n)m = 4, n = 3 (n)mms = 2, nns = −2 (os)
J 1 ( π M ) / π U d c J 3 ( 2 π M ) / ( 2 π Z o , 4 ξ ) 3 U d c J 1 ( π M ) J 3 ( 2 π M ) / ( 2 π 2 Z o , 4 ξ ) 2 θ C 1 _ 1 π / 2
ms = 4, ns = −3 (p)m = 2, n = −1 (p)msm = 2, nsn = −2 (ss)
ms = 4, ns = 3 (n)m = 2, n = 1 (n)msm = 2, nsn = 2 (os)
J 3 ( 2 π M ) / ( 2 π ) U d c J 1 ( π M ) / ( π Z o , 2 ξ ) 3 U d c J 1 ( π M ) J 3 ( 2 π M ) / ( 2 π 2 Z o , 2 ξ ) 2 θ C 1 _ 1 + π / 2
Table 7. Basic grid and converter parameters.
Table 7. Basic grid and converter parameters.
Parameter Value
fg , L g . (grid dispersion inductance)50 Hz, 3 μH
Tr1, υzDy5 (400 V/400 V), √3
Tr2Dd0y5 (400 V/400 V/400 V)
Llq2 (Th1Th12 rectifier side dispersion inductance of Tr2)80 μH
L a , R a ,   I a ,   U a 300 μH, 3 mΩ, 1035 A, 644 V
Firing angle of SCR converter α50°
U d c ,Cdc720 V, 2.7 mF
L_1, L_2, L76 μH, 84 μH, 80 μH
C, Rd, Llq (APF side dispersion inductance of Tr1)0 μF, 180 mΩ, 20 μH
Control system:
fC (carrier frequency)10 kHz, 15 kHz
fsamp (carrier frequency and sampling frequency)20 kHz, 30 kHz
Table 8. Amplitudes and initial phases of individual harmonics in the system with full load.
Table 8. Amplitudes and initial phases of individual harmonics in the system with full load.
k,l U ^ g 1 , k ψ g 1 , k I ^ g 1 , k I ^ L 1 , k φ L 1 , k I ^ 1 , k φ 1 , k I ^ 01 _ 1 , k I ^ c 1 , k I ^ c p , k I ^ υ , l I ^ d c , l
VradAAradAradAAAAA
1326.60.0002162.32156.6−0.9265.91.3982.85.10.00.20.2
20.00.0000.10.21.7100.31.1730.20.00.00.20.2
30.00.0001.01.12.1541.22.7690.60.20.31.31.2
40.00.0000.30.11.5240.52.1930.20.00.00.40.4
50.00.00017.918.72.7162.7−2.901.40.00.00.30.4
60.00.0000.40.10.0000.52.7270.20.00.04.44.4
70.00.0005.810.5−0.9004.6−0.3912.40.00.00.30.4
80.00.0000.70.10.0000.62.5100.40.00.10.30.4
90.00.0002.01.5−0.3982.41.7611.11.01.52.82.5
100.00.0000.80.2−1.5070.71.9810.50.00.10.20.2
1112.0−2.96753.5206.7−0.719169.1−1.46284.64.00.10.50.5
120.00.0000.60.10.0000.1−0.2750.20.10.2163.5163.6
138.0−0.3495.8129.60.440124.30.90962.43.40.20.60.7
140.00.0000.40.1−2.3980.71.8580.30.10.10.60.7
150.00.0001.31.4−2.7622.2−0.7571.01.82.83.82.5
160.00.0000.30.12.9150.21.7780.20.10.10.70.8
170.00.00014.510.4−1.8394.00.3762.00.10.00.90.9
180.00.0001.30.00.0001.2−2.7400.60.10.14.54.6
190.00.00012.18.5−0.4233.73.0702.00.30.31.10.9
200.00.0000.90.10.0000.7−2.5590.30.40.51.10.6
210.00.0001.81.31.0391.9−1.1550.91.21.92.51.6
220.00.0000.50.2−0.0921.4−1.1770.50.30.51.00.9
230.00.00010.070.50.66356.60.15928.23.10.22.02.2
240.00.0001.30.10.0001.2−1.4420.70.30.533.133.0
250.00.00014.757.91.95740.62.56820.33.30.71.92.4
2980.00.0000.00.00.0000.11.17113.10.10.00.10.1
2990.00.0000.40.3−2.0410.10.0001.20.00.00.00.0
3000.00.0000.00.00.0000.00.00032.90.00.00.20.2
3010.00.0000.50.4−0.8160.10.0001.20.00.10.20.0
3020.00.0000.00.00.0000.10.00011.50.10.10.20.1
5500.00.0000.00.00.0000.20.1700.20.20.06.36.2
5620.00.0000.00.00.0000.6−0.3300.30.60.29.49.6
5680.00.0000.00.00.0000.5−1.2470.30.50.37.47.1
5740.00.0000.00.00.0000.33.1230.20.30.36.96.9
5800.00.0000.00.00.0000.2−0.9670.10.20.26.46.5
5890.00.0000.00.00.0002.60.9771.32.64.24.70.5
5920.00.0000.00.00.0000.1−0.0170.10.20.417.016.8
5950.00.0000.10.10.0002.2−0.4101.12.23.52.80.6
5970.00.0000.40.00.0007.5−2.4833.87.80.41.21.5
5980.00.0000.00.00.0001.21.2690.61.30.530.330.8
5990.00.0000.50.00.0009.81.8624.910.10.62.72.1
6000.00.0000.10.00.0000.20.4070.10.20.92.21.3
6010.00.0000.10.10.00011.42.8075.711.417.918.51.1
6020.00.0000.10.00.0001.40.0580.71.41.20.90.4
6030.00.0000.30.00.0007.60.7673.87.90.50.60.9
6040.00.0000.00.00.0000.5−1.3950.20.50.212.112.3
6090.00.0000.10.00.0002.31.2661.22.40.11.71.7
6100.00.0000.00.00.0000.3−2.5000.20.30.025.125.1
6130.00.0000.10.10.0003.0−0.7201.53.04.53.90.7
6160.00.0000.00.00.0000.2−2.0640.10.20.322.322.3
6220.00.0000.00.00.0000.21.6510.10.30.221.521.4
6280.00.0000.00.00.0000.2−0.4420.10.20.05.35.3
6340.00.0000.00.00.0000.32.5510.20.30.18.68.6
6400.00.0000.00.00.0000.32.5460.20.30.16.87.0
6460.00.0000.00.00.0000.10.0000.10.10.17.27.3
6520.00.0000.00.00.0000.10.5190.10.10.15.85.8
Table 9. Further results of simulation research for the system with full load.
Table 9. Further results of simulation research for the system with full load.
I g 1 ( r m s ) I 1 ( r m s ) φ g 1 , 1 I c 1 ( r m s ) I d c ( r m s ) I L 1 ( r m s ) U 12 ( r m s ) U g ( r m s )
1531 A155.1 A−0.93 rad18.4 A138 A1536 A397 V231 V
THDiL1 = 12.3%THDig1 = 4.36%THDu12 = 4.62%THDug1 = 4.42%
Table 10. Results of simulation research and calculations for the system with full load.
Table 10. Results of simulation research and calculations for the system with full load.
Results of Simulation ResearchResults of Calculations
kÛ’g1,kΨ’g1,kÎ’L1,kφ’L1,kIc1,k(rms)λkÎcG1,kÎcL1,kIcGL1,k(rms)
- VradAradA - AAA
132702156.6−0.933.605.130.0323.64
1112−2.97206.7−0.722.80.822.11.982.72
138−0.35129.60.442.40.961.662.12.46
230 70.50.662.20.803.02.12
250 57.92.002.330.702.521.78
Table 11. Results of simulation research and calculations in the system without load.
Table 11. Results of simulation research and calculations in the system without load.
Results of Simulation Research/Results of Calculations
fC [kHz]1510151015101510151015101510
I ^ v C , k I ^ v C , 2 ξ 2 [ A ] I ^ v C , 2 ξ + 4 [ A ] I ^ v C , 4 ξ 4 [ A ] I ^ v C , 4 ξ + 2 [ A ] I ^ v C , 6 ξ 6 [ A ] I ^ v C , 6 ξ [ A ] I ^ v C , 6 ξ + 6 [ A ]
M=
0.9
34.9/
34.2
50.1/
51.3
19.0/
18.4
29.1/
27.7
17.7/
18.5
26.5/
27.8
2.8/
1.2
2.8/
1.9
10.0/
7.9
14.3/
11.9
2.6/
1.8
3.7/
2.7
9.8/
7.9
13.8/
11.9
M=
0.45
25.7/
26.7
39.8/
40.1
2.1/
2.1
3.6/
3.2
7.7/
8.0
12.5/
12.0
30.0/
29.0
45.4/
43.5
3.5/
3.3
5.4/
4.9
14.7/
13.7
22.1/
20.5
3.5/
3.3
5.1/
4.9
M=
0.05
1.06/
0.5
1.4/
0.7
0.07/
0
0.1/
0
0.07/
0
0.0/
0
1.35/
0.9
1.7/
1.3
0.2/
0
0.5/
0
45.5/
43.7
67.7/
65.6
0.2/
0
0.5/
0
Table 12. The currents of LCL-filter and dc-link capacitors and common mode voltage in an PAPF with one and two inverters without load for a positive sequence of fundamental voltages u g 1 , 1 , u g 2 , 1 and u g 3 , 1 .
Table 12. The currents of LCL-filter and dc-link capacitors and common mode voltage in an PAPF with one and two inverters without load for a positive sequence of fundamental voltages u g 1 , 1 , u g 2 , 1 and u g 3 , 1 .
fC = 15 kHzResults for Asymmetrical Regular Sampled PWM/Results for Symmetrical Regular Sampled PWM
PAPF with One Three-Phase Voltage Source InverterPAPF with Two Paralleled Three-Phase Voltage Source Inverters
Converter-Side InductancesL/2L
Carrier Based PWM StrategiesSingle PWM CarrierThree Interleaved PWM Carriers [32]Single PWM CarrierDouble Interleaved PWM Carriers [3]Double Three Interleaved PWM Carriers [13]
Positive Sequence CarrierNegative Sequence Carrier
M0.90.90.6750.450.050.90.6750.450.050.9
Ic,rms [A]61.9/
61.8
67.37/
67.24
62.0/
61.7
17.0/
16.3
19.2/
19.2
18.0/
18.0
2.88/
3.4
17.1/
16.3
19.0/
19.0
17.8/
17.8
2.88/
3.4
17.0/
16.3
Idc,rms [A]53.1/
52.4
56.12/
55.06
53.3/
52.4
31.3/
30.2
41.2/
41.3
64.9/
64.9
123/
121
36.2/
35.2
40.43/
39.38
36.8/
36.9
37.7/
37.4
36.4/
36.2
CMVrms [V]14.4/
15.2
6.32/
7.99
14.4/
15.2
4.4/
8.9
1.17/
8.26
0.60/
8.22
0.04/
8.25
4.33/
8.79
2.57/
8.6
2.5/
8.56
0.39/
8.25
4.3/
8.89
N31662
Table 13. The rms value of currents of LCL-filter and dc-link capacitors for asymmetrical LCL-filter or the grid in an PAPF without load.
Table 13. The rms value of currents of LCL-filter and dc-link capacitors for asymmetrical LCL-filter or the grid in an PAPF without load.
M U g 1 ( r m s ) = U g 2 ( r m s ) = U g 3 ( r m s ) U g 1 ( r m s ) = U g 2 ( r m s ) ,   U g 3 ( r m s ) = 0.9 U g 1 ( r m s )
L_1 = 80µH, L_2 = 80µHL_1 = 76µH, L_2 = 84µHL_1 = 76µH L_2 = 84 µH
IcC1,rmsIdcC,rmsIcC1,rmsIdcC,rmsIcC1,rmsIcC2,rmsIcC3,rmsIdcC,rms
Results of Calculations/Results of Simulations StudiesResults of Simulations Studies
0.922.2/25.1553.0/53.522.88/26.0253.7724.1826.4226.5156.77
0.4524.94/27.5853.0/55.4826.41/29.1355.6329.7728.9928.9956.24
0.053.55/4.4256.5/56.5810.47/11.3056.7511.4111.2711.2656.64

Share and Cite

MDPI and ACS Style

Platek, T. Analysis of Ripple Current in the Capacitors of Active Power Filters. Energies 2019, 12, 4493. https://doi.org/10.3390/en12234493

AMA Style

Platek T. Analysis of Ripple Current in the Capacitors of Active Power Filters. Energies. 2019; 12(23):4493. https://doi.org/10.3390/en12234493

Chicago/Turabian Style

Platek, Tadeusz. 2019. "Analysis of Ripple Current in the Capacitors of Active Power Filters" Energies 12, no. 23: 4493. https://doi.org/10.3390/en12234493

APA Style

Platek, T. (2019). Analysis of Ripple Current in the Capacitors of Active Power Filters. Energies, 12(23), 4493. https://doi.org/10.3390/en12234493

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