Jalili 2009
Jalili 2009
5, MAY 2009
Abstract—This paper introduces a new iterative design proce- M , M0 , Mn Modulation depth, no load modulation
dure of L and LCL filters for low-voltage active-front-end PWM depth, and nominal load modulation depth
two-level voltage source converters. The analytical expression of (M = reference signal amplitude/carrier
the converter harmonic voltages by Bessel functions is applied to √
design the filter parameters for defined maximum grid current signal amplitude, M0 = 2 2Vg,1 /Vdc ).
harmonics. Different filter designs are derived for various reso- Lf ,RLf Inductance and resistance of the L filter
nance frequencies and inductance split factors of LCL filter. The inductor.
minimum amount of stored energy of passive components is used Lf,max Upper limit of L filter inductance.
to select a final filter design. A voltage-oriented control scheme, Lf,g , RLf,g Inductance and resistance of grid-side in-
including active damping, is applied in the (400 V and 50 kVA) ex-
perimental setup. Both simulation and experimental investigations ductor of LCL filter.
are presented to verify the accuracy of the filter design procedure. Lf,conv , RLf,conv Inductance and resistance of converter-
Finally, the steady state and transient performance of the active side inductor of LCL filter.
front end with different LCL filters are depicted. Cf LCL filter capacitance.
Index Terms—Active damping, active front end, LCL filter. fres LCL filter resonance frequency.
r LCL filter split factor, r = Lf,g /Lf,conv .
Cdc DC-link voltage capacitance.
N OMENCLATURE
Vdc,reserve,% DC-link voltage reserve.
x Instantaneous value of quantity x. Iˆg,(mf −2),desired,% Desired amplitude of the (mf − 2)th grid
X, X̂ RMS value and amplitude of quantity x. current harmonic in percent of rated fun-
X1 , X̂1 RMS value and amplitude of fundamental damental component (e.g., according to
component of quantity x. IEEE-519).
Xh ,X̂h RMS value and amplitude of harmonic Iˆg,(mf −2),0,% Amplitude of the (mf − 2)th grid current
components of quantity x. harmonic at no load condition in percent
x Complex vector in α–β coordinate. of rated fundamental component.
x1 Complex vector of fundamental compo- Iˆg,(mf −2),n,% Amplitude of the (mf − 2)th grid current
nents in α–β coordinate. harmonic in percent of fundamental com-
g, conv Subscripts for grid and converter ponent at rated load condition.
quantities. Iˆripple,max Maximum amplitude of current ripple
ll Subscripts for line-to-line quantities. (Iˆripple,max = max |i(t) − i1 (t)|).
fg ,Lg Grid fundamental frequency and grid
stray inductance.
S-PWM Asymmetrical regular sampled sinus tri- I. I NTRODUCTION
angle pulsewidth modulation [29].
ZSS-PWM Asymmetrical regular sampled sinus tri-
angle pulsewidth modulation with added
third harmonic [29].
E MERGING applications like regenerative energy sources,
increasing energy prices, and strict international grid stan-
dards are the reason, therefore, that pulsewidth-modulation
fC , fs Carrier frequency and sampling (PWM) voltage-source converters (VSCs) are increasingly ap-
frequency. plied as active-front-end converters in low-voltage drives. Si-
mf Modulation index (mf = carrier nusoidal input currents, an adjustable power factor at the point
frequency/reference frequency). of common coupling (PCC), ride-through capability, and an
Vdc , Vdc,min DC-link voltage and minimum required adjustment of the dc-link voltage are further advantages of such
dc-link voltage. a converter topology (see, e.g., [1]–[9]).
An inductance is the most simple filter configuration between
PWM voltage source converter and grid. However, the limited
Manuscript received January 29, 2008; revised November 11, 2008. First
published January 6, 2009; current version published April 29, 2009. maximum dc voltage and dynamic performance of a converter
K. Jalili is with the Converteam GmbH, 12277 Berlin, Germany (e-mail: with common switching frequencies, as well as the substantial
kamran.jalili@converteam.com). costs and size of the inductance, avoid the use of a pure induc-
S. Bernet is with Dresden University of Technology, 01069 Dresden,
Germany (e-mail: steffen.bernet@tu-dresden.de). tive filter for medium- and high-power converters (S > x kVA)
Digital Object Identifier 10.1109/TIE.2008.2011251 [10], [11].
0278-0046/$25.00 © 2009 IEEE
JALILI AND BERNET: DESIGN OF LCL FILTERS OF ACTIVE-FRONT-END TWO-LEVEL VSCs 1675
where
vg,1 = Vg,1 ∠0 vconv,1 = Vconv,1 ∠θvconv,1
ig,1 = Ig,1 ∠0 ωg = 2πfg .
Vg,1
Iˆripple,max,pu ≈ √ (11)
2 6fC Lf,t
Vdc − Vdc,min
Vdc,reserve,% = 100 (12)
Vdc
Vdc,min = Vg,1 m 1 + L2f,t,pu (13)
ig (s) 1
= , Lf,g,t = Lg +Lf,g , s = jωg .
Fig. 5. Single-phase representation of the active rectifier with LCL filter. iconv (s) 1+s2 Lf,g,t Cf
(17)
B. LCL Filter
As shown in the previous section, the L filter is only a Substituting the filter inductance values and the resonance
useful solution at very high switching frequencies. An LCL frequency for the filter capacitance, the attenuation factor (17)
filter enables a distinctly cheaper and smaller filter solution at becomes
low switching frequencies (e.g., fC = 4 kHz–8 kHz) which are
Iˆg,(mf −2) kf2
usually applied in industrial low-voltage drives [10], [20].
= 2 (18)
An iterative design procedure of an LCL filter is proposed ˆ
Iconv,(mf −2) kf − r − 1
in this section. The main goal is to limit the (mf − 2)th grid
current harmonic at nominal load to the values defined in for the (mf − 2)th current harmonic, with
IEEE-519. A single-phase representation of the active rectifier
fres
with LCL filter is shown in Fig. 5. The filter split factor r and kf = . (19)
the filter resonance frequency fres as the LCL filter parameters (mf − 2)fg
are given by
Applying (18), the converter-side filter inductance is derived
Lf,g,t for given filter parameters r and fres , on the basis of (7), to
r= , Lf,g,t = Lg + Lf,g (15)
Lf,conv
kf2
Lf,conv,pu = a1 2 (20)
1 Lf,conv + Lf,g,t kf − r − 1
fres = . (16)
2π Lf,conv Lf,g,t Cf
where
The filter resonance frequency (16) is considered to be in a √
V
range between ten times of the grid frequency and one-half of √ J2 π 2q1,−2 Vg,1
dc
a1 = 100 2 V
the carrier frequency in order to avoid a filter excitation by the πq1,−2 (mf − 2)Iˆg,(mf −2),0,% Vg,1
dc
converter voltage harmonics [10].
In the proposed LCL filter design procedure, the converter- Lf,conv Vg,1
Lf,conv,pu = , Lb = .
side filter inductance Lf,conv is determined at no load con- Lb ωg Ig,1
dition M = M0 in a first step. The required inputs of the
design procedure are the filter parameters r and fres , the Iˆg,(mf −2),0,% is the amplitude of the grid harmonic current of
IEEE-519-defined (mf − 2)th grid current harmonic amplitude (mf − 2)th order at no load conditions in percent of the rated
Iˆg,(mf −2),desired,% , and the active rectifier basic parameters grid current amplitude.
(see, e.g., Table I). The other passive filter components Lf,g Applying the determined converter-side filter inductance (20)
and Cf are calculated then based on Lf,conv and the filter and the filter split factor (15), the grid-side filter inductance is
parameters defined in (15) and (16). The Thevenin equivalent
kf2
circuit of the converter-side LC-grid is used in the next step
Lf,g,t,pu = a1 r 2 , Lf,g,pu = Lf,g,t,pu − Lg,pu .
to model the LCL filter as an L filter. The (mf − 2)th grid kf − r − 1
current harmonic amplitude Iˆg,(mf −2),n,% for rated grid current (21)
is determined by applying this model and the calculated passive
filter components in order to verify the filter effectiveness at By using (16), the filter capacitance is
nominal load. Finally, the required passive filter components
r2 +r(2−kf2 )+(1−kf2 ) 1
for nominal load conditions are calculated in an iterative man- Cf,pu = a2 , a2 = (22)
ner by applying the calculated deviation ΔIˆg,(mf −2),0,% = rkf 4 (mf −2)2 a1
Iˆg,(mf −2),n,% − Iˆg,(mf −2),desired,% .
for a chosen filter resonance frequency and split factor with
The converter-side filter inductance Lf,conv is determined
on the basis of (7) to limit the (mf − 2)th converter current Cf Ig,1
harmonic at no load condition like described before for the Cf,pu = , Cf,b = .
Cf,b ωg Vg,1
design of the L filter. The filter capacitors are assumed to have
zero impedance for the current harmonics at the switching fre- The modulation depth at nominal load condition Mn can
quency. The amplitude of the (mf − 2)th converter-side current be specified using the calculated passive component values
1680 IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS, VOL. 56, NO. 5, MAY 2009
Vdc,reserve,% = 100
Vg,1
× 1− 2
m XLC,f +(Lf,g,t,pu XLC,fg +Lf,conv,pu )2
Vdc g
(26)
3 2
Wt = I (Lf,conv + Lf,g,t ) + Cf Vg2 . (28)
2 g,1
Fig. 8. Filter components for various values of r and fres to fulfill the IEEE-519 condition; Vg,ll = 400 V, Vdc = 700 V, mf = 81. (a) The converter-side
inductance, (b) the total grid-side inductance, (c) the filter total inductance, (d) the filter capacitance, (e) the amplitude of the grid harmonic current, and (f) the
converter current ripple.
Fig. 9. Stored energy of the reactive components of the LCL filters shown in Fig. 8; Vg,ll = 400 V, Ig,1 = 70 A, Vdc = 700 V, mf = 81. (a) Energy of the
total inductance, (b) energy of the filter capacitance, and (c) the filter total energy.
JALILI AND BERNET: DESIGN OF LCL FILTERS OF ACTIVE-FRONT-END TWO-LEVEL VSCs 1683
TABLE IV
LCL FILTER PARAMETERS AND REACTIVE COMPONENT VALUES CORRESPONDING TO THE MINIMUM FILTER ENERGY FOR DIFFERENT RATED GRID
CURRENTS, DC-LINK VOLTAGE, AND LINE VOLTAGE (MODULATION: S-PWM, Ig,(mf −2),desired,% = 0.2)
TABLE V
LCL FILTER PARAMETERS AND REACTIVE COMPONENT VALUES CORRESPONDING TO THE MINIMUM FILTER ENERGY FOR DIFFERENT RATED GRID
CURRENTS, DC-LINK VOLTAGE, AND LINE VOLTAGE (MODULATION: ZSS-PWM, Ig,(mf −2),desired,% = 0.2)
procedure of Fig. 7 which induces different control reserves. active rectifier with LCL filter in Fig. 5, the current of the
The total filter inductance increases, as shown in Fig. 8, for converter-side inductor can be described by
filters with increasing resonance frequencies, and consequently,
d
the control reserve decreases, leading to a decreasing dynamic Lf,conv iconv + RLf,conviconv = vCf − vconv . (29)
performance. The dc-link voltage can be changed considering dt
the control reserve (26) in order to achieve a defined dc-link Equation (29) is transformed into the rotating frame. The
voltage reserve (Fig. 10). As mentioned in the L filter design, it corresponding d–q components are given in (30), shown at the
should be noticed that the maximum value of the dc-link voltage bottom of the page. The following decoupling terms should be
is limited by the blocking characteristics of the semiconductors added to the output of the current controllers to get a first-order
and the converter design. current control loop:
Δvd = vCf ,d + Lf,conv ωg iconv,q
IV. C ONTROL S YSTEM (31)
Δvq = vCf ,d − Lf,conv ωg iconv,d .
In order to evaluate the proposed design procedure of LCL
filter with simulation and experimental results, a closed-loop The current control loop with corresponding processing and
control system with active damping is applied. The voltage- PWM delays [17] is shown in Fig. 12. The parameters of the
oriented control (VOC) applying a PI current control was cho- current control loop are
sen as a control scheme of the active rectifier [20]–[22], [25].
1 Lf,conv
The converter currents are transformed into the rotating frame K= , τ=
synchronous to the voltage vector of the filter capacitors. The RLf,conv RLf,conv
converter-side inductor current is controlled using the capacitor Lf,conv
Kp = , Ti = τ, a≥2
voltages and the converter currents as feedback. The main 1.5aTs
advantage of such a control is the robustness of the current 1
controller against the grid inductance variation and the effect Ts = : sampling time. (32)
fs
of the active damping based on the capacitor voltages [20].
A block diagram of the applied control system is shown Fig. 12 shows that it is not possible to achieve a completely
in Fig. 11. Considering the single-phase representation of the decoupled system due to the current control loop delay. The
d
Lf,conv dt iconv,d + RLf,conv iconv,d = −vconv,d + vCf ,d + Lf,conv ωg iconv,q
(30)
d
Lf,conv dt iconv,q + RLf,conv iconv,q = −vconv,q + vCf ,q − Lf,conv ωg iconv,d
1684 IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS, VOL. 56, NO. 5, MAY 2009
Fig. 12. Converter current control loop including the processing and PWM
delays.
Fig. 14. Grid and converter currents and the corresponding harmonic spectra. (a) Lf,conv = 0.711 mH, Lf,g = 0.703 mH, Cf = 69.418 μF, fres = 1016 Hz,
(b) Lf,conv = 1.33 mH, Lf,g = 1.33 mH, Cf = 20.92 μF, fres = 1350 Hz. (a-1) and (b-1) The grid and the converter currents. (a-2) and (b-2) The
corresponding harmonic spectra.
TABLE VII
SIMULATION RESULTS FOR THE AMPLITUDE OF THE CONVERTER
CURRENT HARMONICS AND THE CONTROL RESERVE FOR
TWO SETS OF THE LCL FILTER
Fig. 16. Experimental results for the grid and the converter currents and the corresponding harmonic spectra for S-PWM. (a) Lf,conv = 0.6 mH, Lf,g =
0.6 mH, Cf = 88 μF, fres = 980 Hz: (a-1) The grid and the converter currents and (a-2) the harmonic spectra of the grid and the converter currents.
(b) Lf,conv = 1.2 mH, Lf,g = 1.2 mH, Cf = 22 μF, fres = 1385 Hz: (b-1) The grid and the converter currents and (b-2) the corresponding harmonic spectra.
TABLE VIII
SIMULATION AND EXPERIMENTAL RESULTS FOR THE AMPLITUDE OF THE GRID CURRENT HARMONICS FOR TWO SETS OF LCL FILTERS
VII. C ONCLUSION
This paper has proposed an iterative procedure to design L
and LCL filters for active rectifiers. The analytical expression
The experimental results for the step response of the current of the converter voltage harmonics based on Bessel functions
control loop are shown in Fig. 17. The performance of the is applied to determine the filter parameters which enable suffi-
current controller for a step of 100 A on the q-axis for the filters cient damping of the grid current harmonics (e.g., according to
with 980 and 1385 Hz as resonance frequencies is shown in IEEE-519-1992).
Fig. 17(a) and (b), respectively. The capacitor voltage feedback The LCL filter is modeled as L filter. The analytical expres-
is considered to be without compensator for waveform 1, with sion of the amplitude of the significant grid current harmonics
a lead–lag compensator and α = 4/3π for waveform 2, and at nominal load and unity power factor is presented. It is
with a lead–lag compensator and α = 2/3π for waveform 3. shown that LCL filters with different parameters (fres , r) and
1688 IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS, VOL. 56, NO. 5, MAY 2009
[28] M. Malinowski and S. Bernet, “A simple voltage sensorless active damp- Steffen Bernet (M’97) received the Diploma de-
ing scheme for three-phase PWM converters with an LCL-filter,” IEEE gree in electrical engineering from Dresden Uni-
Trans. Ind. Electron., vol. 5, no. 4, pp. 1876–1880, Apr. 2008. versity of Technology, Dresden, Germany, in 1990,
[29] D. G. Holmes and T. A. Lipo, Pulse Width Modulation for Power and the Ph.D. degree in electrical engineering
Converters. New York: Wiley, 2003. from Ilmenau University of Technology, Ilmenau,
[30] U. Nicolai, Application Manual Power Modules. Ilmenau, Germany: Germany, in 1995.
SEMIKRON Int., 2000. During 1995 and 1996, he was a Postdoctoral
[31] V. Blasko, “Adaptive filtering for selective elimination of higher harmon- Researcher with the Electrical and Computer En-
ics from line currents of a voltage source converter,” in Conf. Rec. IEEE gineering Department, University of Wisconsin,
IAS Annu. Meeting, 1998, vol. 2, pp. 1222–1228. Madison. In 1996, he was with ABB Corporate
[32] L. Moran, P. D. Ziogas, and G. Joos, “Design aspects of synchro- Research, Heidelberg, Germany, where he led the
nous PWM rectifier-inverter systems under unbalanced input voltage Electrical Drive Systems Group. From 1999 to 2000, he was responsible for
conditions,” IEEE Trans. Ind. Appl., vol. 28, no. 6, pp. 1268–1293, ABB research in the areas of power electronics systems, drives, and electric
Nov./Dec. 1992. machines. From 2001 to 2007, he was a Professor of power electronics at Berlin
[33] P. Rioual, H. Pouliquen, and J. Louis, “Regulation of a PWM rectifier in University of Technology, Berlin, Germany. Since June 2007, he has been a
the unbalanced network state using a generalized model,” IEEE Trans. Professor at Dresden University of Technology.
Ind. Electron., vol. 11, no. 3, pp. 495–502, May 1996.