# Research on Discretization PI Control Technology of Single-Phase Grid-Connected Inverter with LCL Filter

Authors: Jianke Li; Jinquan Wang; Ye Xu; Jianting Li; Jingjing Chen; Pengfei Hou; Shuhua Qian
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Compared with L-type filter, LCL-type filter is more suitable for high-power low-switching frequency applications with reducing the inductance, improving dynamic performance. However, the parameter design for the LCL filter is more complex due to the influence of the controller response performance of the converter. If the harmonic current around switching frequency can be fully suppressed, it is possible for inverter to decrease the total inductance as well as the size and the cost. In this paper, the model of the LCL filter is analyzed and numerical algorithms are adopted to analyze the stability of the closed-loop control system and stable regions are deduced with different parameters of LCL filter. Then, the minimum sampling frequencies are deduced with different conditions. Simulation and experimental results are provided to validate the research on the generating mechanism for the unstable region of sampling frequency.


---

## Body

## 1. Introduction

Three types of filter, such as L-type, LC-type, and LCL-type, are widely used in grid-connected inverter. Due to the limit of switching frequency, the filter inductance of L-type grid-connected inverter cannot effectively suppress the harmonic voltage of PWM switching frequency, resulting in grid current with large harmonic current around switching frequency, which should be suppressed by larger filter inductance [1–4]. However, larger filter inductance will cause too large equipment volume and high cost and affect dynamic performance of grid-connected control [5]. Therefore, the inductance-capacitance-inductance (LCL) filter is designed to replace the conventional filter [6, 7]. In order to improve the attenuation rate of harmonic current around switching frequency, LCL filter is adopted to suppress harmonic voltage around switching frequency as well [8]. However, LCL filter brings new problems, such as the scope of application of LCL filter and the impact on control characteristics caused by LCL filter resonance.In order to research the stability of inverter with LCL filter, it is important to use simulation models which are sufficiently detailed to realistically represent their real world physical system behaviors. A detailed power inverter model developed within the DIgSILENT power network simulation package, which can reflect the inverter’s real dynamic response to transient events, is presented in [9]. A new modeling approach for inverter-dominated microgrids using dynamic phasor is presented. The proposed dynamic phasor model is able to predict accurately the stability margins of the system, while the conventional reduced-order small-signal model fails [10]. In order to predict the dynamic behaviors of inverter, new small-signal z-domain models are deduced for digitally controlled grid-connected inverters with converter current control scheme and converter current plus grid current control scheme [11]. The proposed methods allow direct design for controllers in z-domain. The simulation results show that the proposed z-domain models are more effective in predicting instabilities. As to control method, a PI controller with self-tuning parameter based on fuzzy inferring is proposed [12]. This controller is capable of automatically adjusting two parameters (P and I) of PI controller.For an LCL filter based single-phase grid-connected full-bridge inverter system, it is possible to decrease the total inductance as well as the size and the cost, if the harmonic current around switching frequency can be fully suppressed. LCL filters resonance may lead to the instability of the control system. In order to address this issue, passive damping and active damping have been presented to improve system stability [13–15]. In order to cope with the grid inductance variations, a simple tuning procedure for the notch filter is proposed to estimate the resonance frequency by means of Fourier analysis. The Goertzel algorithm, instead of the FFT, is used to reduce the calculation and memory requirements. Thus, the proposed self-commissioning notch filter increases and consumes little computational resources [16]. Improved passive damping which includes double loop control and makes the system less loss and more stable is proposed [17]. An active damping strategy with harmonics compensation which can alleviate the harmonics around the resonance frequency caused by the LCL filters is proposed in [18]. However, whether LCL filter can effectively suppress the harmonic current with different switching frequencies is not deeply analyzed [19].Since stability margin of grid current feedback control is small, double closed-loop control, capacitor current feedback inner loop, and grid current feedback outer loop are adopted to increase the stability margin of the system [20]. The inverter side inductance current feedback, which is inner inductance current feedback, is proposed in [8], and they found that the stability margin using this current feedback is larger than current closed-loop control. Besides, the inner inductance current feedback with advantages of simple control algorithm and less feedback parameters has already been applied in [21–23]. In order to enhance the tracking characteristics of grid current, reference current feed forward control is presented. However, whether inner inductor current feedback control can meet grid current tracking features alone is not discussed [8].This paper will analyze mathematical model of the LCL filter and deeply analyze the stability of inner inductance current control system.


## 2. Mathematical Model of Filter

Figure1 shows the schematic of single-phase full-bridge PWM inverter with LCL filter. T
1, T
2, T
3, and T
4 are power MOSFETs [24]. L and R are the inner inductance and parasitic resistance. L
g and R
g are the grid-side inductance and parasitic resistance. u
d and u
b are the DC bus voltage and the inverter output voltage, respectively. u
s and u
c are the grid voltage and capacity voltage, respectively. i
L and i
g are the inner inductance current and grid current, respectively. The PWM inverter can be equivalent to proportion enlargement link K which is generally normalized to 1 [25]. And the voltage u
b between the two bridges can be substituted by reference wave voltage u
i (Figure 2).Figure 1
Schematic of LCL based grid-connected inverter.Figure 2
The transfer function block diagram of LCL filter.The electrical relationship of the schematic can be described as follows:(1)
[
i
˙
L
u
˙
c
i
˙
g
]
=
A
[
i
L
u
c
i
g
]
+
B
[
u
i
u
s
]
,
where
(2)
A
=
[
-
R
L
-
1
L
0
1
C
0
-
1
C
0
1
L
g
-
R
g
L
g
]
,
B
=
[
1
L
0
0
0
0
-
1
L
g
]
.
The transfer function of LCL filter which can be derived from formula (1) and Figure 1 is shown in Figure 2.The dead zone effect can be equal to dead zone equivalent voltage sourceu
de; the transfer function block diagram and state equation are shown in Figure 3 and formula (3), respectively. Consider
(3)
[
i
˙
L
u
˙
c
i
˙
g
]
=
A
[
i
L
u
c
i
g
]
+
B
[
u
i
u
s
u
de
]
,
where
(4)
A
=
[
-
R
L
-
1
L
0
1
C
0
-
1
C
0
1
L
g
-
R
g
L
g
]
,
B
=
[
1
L
0
1
L
0
0
0
0
-
1
L
g
0
]
.Figure 3
The transfer function block diagram of the LCL filter with dead zone effect.The continuous domain transfer functions of the inner inductor current and the grid current are derived from formula (3) as shown below:
(5)
I
L
(
S
)
=
A
(
S
)
D
(
S
)
(
U
i
(
S
)
+
U
de
(
S
)
)
-
1
D
(
S
)
U
s
(
S
)
,
I
g
(
S
)
=
1
D
(
S
)
(
U
i
(
S
)
+
U
de
(
S
)
)
-
A
(
S
)
D
(
S
)
U
s
(
S
)
,
where
(6)
A
(
S
)
=
L
g
C
S
2
+
R
g
C
S
+
1
,
D
(
S
)
=
L
L
g
C
S
3
+
(
L
g
R
+
L
R
g
)
C
S
2
+
(
R
R
g
C
+
L
+
L
g
)
S
+
R
+
R
g
.The output voltage of inverter keeps a constant valueu
i
(
k
T
) in kth sampling period T, until the next sampling time (
k
+
1
); the output voltage becomes u
i
(
(
k
+
1
)
T
). The discretization transfer function is shown below:
(7)
x
(
k
+
1
)
=
G
x
(
k
)
+
H
u
(
k
)
,
y
(
k
)
=
i
g
(
k
)
,
G
=
e
-
A
T
,
H
=
A
-
1
(
1
-
e
-
A
T
)
B
,
x
(
k
)
=
[
i
L
(
k
)
u
c
(
k
)
i
g
(
k
)
]
,
u
(
k
)
=
[
u
i
(
k
)
u
s
(
k
)
u
de
(
k
)
]
.The block diagram of grid-connected inverter with LCL filter is shown in Figure4.Figure 4
The discretization transfer function block diagram of LCL filter.


## 3. Discretization PI Control Technology

The discrete domain block diagram of PI control technology can be obtained from Figure4, as shown in Figure 5.Figure 5
The discrete domain block diagram.The state space equation of PI control system can be obtained as shown below [26]:
(8)
X
(
k
+
1
)
=
G
I
X
(
k
)
+
H
1
I
i
ref
(
k
)
+
H
2
I
u
2
(
k
)
,
Y
(
k
)
=
C
1
I
X
(
k
)
,
C
1
I
=
[
0
0
1
0
]
,
G
I
=
[
G
-
(
K
p
+
K
I
)
H
1
C
3
H
1
-
K
I
C
3
1
]
,
X
(
k
)
=
[
i
L
(
k
)
u
c
(
k
)
i
g
(
k
)
e
I
(
k
)
]
,
H
1
I
=
[
(
K
p
+
K
I
)
H
1
K
I
]
,
H
2
I
=
[
H
2
0
]
.In order to analyze the stability of closed-loop control system, numerical algorithm is adopted to describe the region where the root is less than 1. Six filters with various ratios between inner inductance and outside inductance are selected to evaluate the range of stability domain, as shown in Table1.Table 1
Six group filters.
Group
L/mH
L
g/mH
C/uF
R/mΩ
R
g/mΩ
1st
1
1
4.7
10
10
2nd
1.5
0.5
6.3
5
15
3rd
0.5
1.5
6.3
5
15
4th
0.5
0.5
9.4
5
15
5th
0.75
0.25
12.6
7.5
2.5
6th
0.25
0.75
12.6
2.5
7.5In these six group filters, the resonant frequencyf
res is 3.3 kHz and closed-loop damping ratio coefficient ξ of closed-loop control system is 0.707. The curves of closed-loop root critical frequency f
n
max
⁡ with different sampling frequencies can be obtained according to formula (9), as shown in Figure 6:
(9)
0
<
f
n
<
1
+
ξ
2
-
ξ
π
T
=
1
+
ξ
2
-
ξ
π
f
s
=
f
n
max
⁡
.Figure 6
Stability domain of filters.As shown in Figure6, the closed-loop root critical frequency f
n
max
⁡ almost equals zero in several groups. These regions are defined as the unstable region of sampling frequency. From Figure 6, the unstable region of sampling frequency (ribbon) is present in all these filters. The curves of closed-loop root critical frequency f
n
max
⁡, respectively, overlap between group 1 and group 4, group 2 and group 5, and group 3 and group 6, which indicates that the ranges of closed-loop root frequency are identical when the resonance frequencies are the same and the ratios between inner inductance and outside inductance are the same.The upper and lower frequency can be obtained from Figure6, as shown in Table 2. The ratios K
ufs between upper frequency and resonant frequency and the ratios K
ufs between lower frequency and resonant frequency are listed in Table 2 as well.Table 2
The unstable region under different parameters.
Unstable region 1/Hz
Unstable region 2/Hz
Unstable region 3/Hz
Unstable region 4/Hz
Upper
Lower
Upper
Lower
Upper
Lower
Upper
Lower
1st
6600
3450
2175
1720
1310
1140
940
855
2nd
6600
3500
2175
1750
1310
1160
940
855
3rd
6600
3350
2175
1680
1310
1120
940
840
4th
6600
3450
2175
1720
1310
1140
940
855
5th
6600
3500
2175
1750
1310
1160
940
855
6th
6600
3350
2175
1680
1310
1120
940
840
K
ufs
2/1
2/2
2/3
2/4
2/5
2/6
2/7
2/8From Table2, the unstable region can be described below:
(10)
2
f
res
2
n
<
f
s
≤
2
f
res
2
n
-
1
,
n
∈
N
.In order to verify whether the stability region satisfies formula (10) with different damping ratio coefficients, a group parameter, L
=
1 mH, L
g
=
1 mH, C
=
4.7 uF, R
=
10 mΩ, and R
g
=
10 mΩ, is selected. The curves of closed-loop root critical frequency f
n
max
⁡ are drawn in Figure 7 when the damping ratio coefficient equals 0, 0.2, 0.4, 0.6, 0.8, 1.0, and 1.2, respectively.Figure 7
Closed-loop root critical frequencyf
n
max
⁡ with different damping ratio coefficients.As shown in Figure7, the frequency ranges of unstable region all satisfy formula (10) with different damping ratio coefficients except when ξ is 0. The ranges of closed-loop root frequency will gradually increase when the sampling frequency is twice as high as the resonant frequency. What is more, the range of closed-loop root frequency will be consistent with discrete PI control of L-type filter when the sampling frequency is three times as high as the resonant frequency.In accordance with quantitative requirements of closed-loop root frequency of grid-connected standard, taking into account the stability margin, according to the method of the minimum sampling frequency, the conclusion can be deduced: the sampling frequency should be greater than 9.68 kHz without the grid voltage feed-forward to satisfy the grid-connected standard indicators, while the sampling frequency should be greater than 3.9 kHz with the grid voltage feed-forward.


## 4. Mechanism Research for the Unstable Region of Sampling Frequency

For the same LCL filter with resonant frequencyf
res, we can know from previous sections that a sampling frequency unstable region exists when the filters are controlled by PI method and this area meets formula (10) and has no relation with the closed-loop damping ratio. It is shown that the unstable strip area is related to the discretization of control object. It will be analyzed as follows according to the discretization process of control object.The discretization process of control objects is mainly manifested in the transfer matrixG of control object, and the transfer matrix H is related to input signal; the relationship between matrix functions G and H is represented below:
(11)
G
=
e
-
A
T
,
H
=
A
-
1
(
1
-
e
-
A
T
)
B
,
where
(12)
A
=
[
-
R
L
-
1
L
0
1
C
0
-
1
C
0
1
L
g
-
R
g
L
g
]
,
B
=
[
1
L
0
1
L
0
0
0
0
-
1
L
g
0
]
.The characteristics root of matrixA can be represented by formula (13); since the parasitic resistance of the internal and external inductance is too small, the influence of parasitic resistance can be ignored; the parameters of a, b, and c in formula (13) are shown in formula (14):
(13)
λ
1
=
a
+
b
i
,
λ
2
=
a
-
b
i
,
λ
3
=
c
,
(14)
a
=
c
=
0
,
b
=
L
+
L
g
L
L
g
C
=
2
π
f
res
.The coefficient relations of minimum polynomial of the transfer matrixG are described in formula (15) and the coefficient expressions of minimum polynomial are shown in formula (16):
(15)
d
0
+
2
π
f
res
d
1
-
4
π
2
f
res
2
d
2
=
e
-
2
π
f
res
T
i
,
d
0
-
2
π
f
res
d
1
i
-
4
π
2
f
res
2
d
2
=
e
2
π
f
res
T
i
,
d
0
=
e
0
,
(16)
d
0
=
1
,
d
1
=
-
sin
⁡
(
2
π
f
res
T
)
2
π
f
res
,
d
2
=
1
-
cos
⁡
⁡
(
2
π
f
res
T
)
4
π
2
f
res
2
.The expression of matrix functionsG and H can be deduced as shown below:
(17)
G
=
I
-
sin
(
2
π
f
res
T
)
2
π
f
res
A
+
1
-
cos
⁡
(
2
π
f
res
T
)
4
π
2
f
res
2
A
2
,
H
=
(
sin
⁡
(
2
π
f
res
T
)
2
π
f
res
I
-
1
-
cos
⁡
⁡
(
2
π
f
res
T
)
4
π
2
f
res
2
A
)
×
[
1
L
0
1
L
0
0
0
0
-
1
L
g
0
]
.From formula (17), the relationship matrix between output voltage u
i of inverter and state quantity of LCL filter parameters can be inferred as shown below:
(18)
H
1
=
[
sin
(
2
π
f
res
T
)
2
π
f
res
L
1
-
cos
⁡
(
2
π
f
res
T
)
4
π
2
f
res
2
L
C
0
]
.From formula (18), the output voltage of inverter is positive and the inner inductance current increases when sin
(
2
π
f
res
T
) is greater than zero. When sin
(
2
π
f
res
T
) is less than zero, the output voltage of inverter is positive, but the inner inductance current decreases. However, when designing the inner inductance current closed-loop control, positive output voltage and the enlargement of inner inductance current are considered. Therefore, the system feedback turns into positive from negative, which results in unstable system when sin
(
2
π
f
res
T
) is less than zero. That is, the discretization closed-loop control system will become unstable when formula (19) is satisfied:
(19)
sin
(
2
π
f
res
f
s
)
<
0
.The expression of the unstable region of sampling frequency can be inferred from formula (19), just as formula (10).Therefore, the research for discretization PI control technology in the above section is demonstrated in theory. The mechanism of the unstable region of sampling frequency can be summarized as follows: due to the discretization of LCL filter, the feedback polarity of closed-loop control will change from negative to positive, which results in unstable control system when the sampling frequency is in unstable region.


## 5. Simulation Verification

In order to verify the theoretical analysis, MATLAB/Simulink software is adopted to the simulation analyses. According to earlier designs, the simulation parameters of LCL filter,L
=
1.2 mH, L
g
=
0.5 mH, C
=
8 uF, R
=
12 mΩ, and R
g
=
5 mΩ, are elected. The resonant frequency f
res is 3 kHz and the network voltage is 220 V at 50 Hz. The inverter bridge can be regarded as a controllable voltage source during the simulation.We can know from formula (10) that when the resonant frequency is 3 kHz, the maximum sampling frequency in the unstable region (ribbon) under discrete PI control is 6 kHz. Therefore, this paper chooses 10 kHz as the typical value in the stable region and 6 kHz as the typical value in the unstable region. The damping coefficient of closed-loop control system is 0.707, and the sampling frequency is 10 kHz. Figure 8 presents the simulation results when the sampling frequency is 6 kHz.The simulation results.
(a)
Without limit of output voltage
(b)
380 V limiter of output voltage
(c)
Output voltage with PWM waveform of bipolar modulation
(d)
Steady-state waveform zoom of (c)
(e)
Output voltage with sampling frequency within unstable regionThe closed-loop root critical frequencyf
n
max
⁡ of stable region is 1647 Hz when the sampling frequency is 10 kHz. In Figures 8(a)–8(c), the closed-loop root frequency closes to 1447 Hz in the time periods (0, 0.1) and (0.2, 0.5) and to 1747 Hz in the time period (0.1, 0.2).Comparing Figure8(a) with Figure 8(b), the wave in Figure 8(c) has a little distortion resulting from the switching frequency harmonic of PWM.The steady-state waveform in Figure8(c) is enlarged as shown in Figure 8(d) from which we can know that there are only harmonics around the switching frequency except for fundamental current in the current wave.The closed-loop root critical frequencyf
n
max
⁡ is 988 Hz when the sampling frequency is 6 kHz. In Figure 8(e), the closed-loop root frequency closes to 329 Hz in the time periods (0, 0.1) and (0.2, 0.5) and to 1088 Hz in the time period (0.1, 0.2). It can be seen from Figure 8(e) that when the sampling frequency is in the unstable region, no matter how much the closed-loop root frequency is, the closed-loop system is unstable, which verifies the theoretical analysis of the unstable region of sampling frequency.


## 6. Experimental Verification

In order to verify the theoretical analysis experimentally, grid-connected photovoltaic inverter 2.5 kW is selected as a control object. The rated voltage and current of grid are 220 V at 50 Hz and 11.4 A, respectively. The inductance is 2 mH and the switching frequency is 10 kHz; the sampling frequency is 20 kHz and the closed-loop damping ratio coefficient is 0.707.All experiment results are presented in Figures9 and 10. Each experiment shows two waveforms; the voltage and current waveforms without any filtering are shown in Figure 9(a), while only the current waveform shown in Figure 9(b) is measured with a 5 kHz low-pass filter. Figure 9 shows the voltage and current waveforms when the closed-loop root frequency is 750 Hz with one-step-delay while Figure 10 shows the results when the closed-loop root frequency is 1400 Hz.The experiment waves with 750 Hz.
(a)
Without filtering
(b)
After 5 kHz low-pass filterThe experiment waves with 1400 Hz.
(a)
Without filtering
(b)
After 5 kHz low-pass filterFigures9 and 10 show that the steady-state characteristic of one-step-delay is good when the closed-loop root frequency is 750 Hz and the output current has begun to oscillate and enter to the critical stability region when the closed-loop root frequency is 1400 Hz. The stable region of PI control can be considered to coincide with the theoretical calculation when considering the error of inherent parameters of inverter.The experimental results show that the theoretical derivation of the stable range of digital control parameter (closed-loop single frequency) is accurate.


## 7. Conclusion

Since the discrete PI control may cause stability problems, this paper analyzes the range of closed-loop root frequency with discrete PI control system to ensure the stability. The following conclusions are obtained: the unstable region of sampling frequency is present in PI control system and the range of closed-loop root frequency will be consistent with discrete PI control of LCL-type filter when the sampling frequency is away from the unstable region. When sampling frequency is in unstable region, the feedback polarity of closed-loop control resulting from the discretization of LCL filter will change from negative to positive, which results in an unstable control system. In practical applications, in order to meet the grid-connected standard indicators, the sampling frequency should be set larger than 9.68 kHz for filters without grid voltage feed forward, but 3.9 kHz for filters with grid voltage feed forward.



---

*Source: 129619-2014-09-29.xml*