An Improved Active Damping Method Based on Single-Loop Inverter Current Control for LCL Resonance in Grid-Connected Inverters

This paper investigates active damping of LCL filter resonance in grid-connected inverters with only inverter current feedback control, since it only needs to sample one current to realize both current control and inverter protection. The traditional single-loop inverter current control (SLICC) can damp the LCL filter resonance actively. However, if the control delay is considered in digital control, the system stability will depend on the ratio of the LCL resonance frequency 
 
 
 
 f
 
 
 res
 
 
 
 to the sampling frequency 
 
 
 
 f
 
 
 s
 
 
 
 , and the valid damping region is only up to 
 
 
 
 
 f
 
 
 s
 
 
 
 /
 6
 
 . Considering that the design region of the LCL resonance frequency 
 
 
 
 f
 
 
 res
 
 
 
 is up to 
 
 
 
 
 f
 
 
 s
 
 
 
 /
 2
 
 , the system can easily become unstable due to the LCL resonance frequency shifting. Thus, this paper proposes an improved active damping method based on SLICC, including the asymmetric regular sampling method and delay compensation method. The improved sampling method minimizes the control delay without introducing a switching ripple, and the delay compensation method further compensates for the delay effect. With a proper parameter design, the upper limit of the valid damping region is extended up to 
 
 
 
 
 f
 
 
 s
 
 
 
 /
 2
 
 , which can cover all the possible resonance frequencies, and it has inherent robustness against grid-impedance variation. Finally, a few simulations in MATLAB/SIMULINK and experiments based on a 6 kW prototype are performed to verify the theoretical analysis.


Introduction
Recently, the renewable energy distributed generation technique is attracting more and more attention, to deal with the increasing concern of the exhaustion of fossil fuels and the environmental pollution. Grid-connected inverters are the key devices to connect the distributed systems to the power grid [1]. As the grid-connected inverter is controlled by the pulse width modulation (PWM), the output current contains an amount of switching harmonics. e LCL filter is widely used due to its better attenuation capability of the switching harmonics and cost advantage compared with L filter. However, the LCL filter is a third-order low-damping system, which challenges the stability of the whole system [2].
In order to stabilize the system, a variety of active damping techniques without power losses have been developed. Generally, the traditional active damping techniques can be classified into two modes based on whether an additional sensor is needed, except for the essential sensors of the control current and grid voltage: the single-loop and the multiloop feedback control. e multiloop feedback control methods introduce additional damping terms into the current control loop by sensing additional variables, such as capacitor current [3][4][5][6] or capacitor voltage [7]. Due to their flexible and efficient implementation, these methods have been widely used. However, requiring additional sensors increases the hardware cost and affects the system's reliability. Alternatively, a state observer can be used to simplify the system hardware effectively [8,9], but it is sensitive to system parameter uncertainties and variations [10].
Simultaneously, the single-loop feedback control methods can avoid these problems and have increasingly been studied in [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25], including inverter current control and grid current control. In [11], it is first pointed out that the inverter current contains capacitor current component so that the inverter current control can promise stable operation. Moreover, inverter current control is beneficial to the instantaneous overcurrent protection for the inverter, which is more convenient to be implemented and popular in industrial applications [12]. But good damping can be achieved only when the resonance frequency is low enough. In order to increase the damping, an active damping inner loop control is introduced. It is implemented by extracting the high-frequency components of the inverter current or grid current to simulate the capacitor current [13][14][15], and the first-order High Pass Filter (HPF) is the most popular filter. However, when the control current is the grid current, it may not be able to introduce sufficient damping since the high-frequency components have been filtered by the filter capacitor, which might result in the need to compromise between the resonance frequency offset and the resonance suppression ability [16].
On the other hand, the control delay is introduced in the digital control implementation process, and it has a significant influence on all active damping techniques [17][18][19]. e system stability will depend on the ratio of the LCL resonance frequency f res to the sampling frequency f s . For the single-loop feedback control without active damping inner loop or the multiloop feedback control, the critical frequency is f s /6 [17,18], while it can increase to f s /3 at most when adopting the active damping inner loop [20].
us, the LCL resonance frequency must be away from the critical frequency to enhance the system damping. Unfortunately, considering the potential variation of the grid inductance, the LCL resonance frequency may shift to the critical frequency; thus, the system stability cannot be promised [21][22][23]. Hence, how to decrease the control delay and increase the critical frequency is a research focus in recent years. e control delay is reduced by adjusting the sampling instant [17] or multisampling [12], but the sampling instant may be no longer the switching instant of the power tube. It is easy to introduce the switching frequency harmonics into the control loop and deteriorate the system stability [16]. Alternatively, the filter-based methods are presented in [24][25][26]. By inserting a digital filter in the forward path, the frequency characteristics around the LCL resonance frequency can be changed. But such a solution relies on the accuk6rate filter parameters; the variation of grid inductance may greatly affect the system stability. us, it is very difficult to enhance the system stability only by reducing the control delay or introducing digital filters. erefore, in [22], a robust method including the double-sampling double-loading (DSDL) mode and a lead grid current control part is proposed, and a similar method based on the inverter current control is proposed [16]. e DSDL mode minimizes the control delay and the lead control part further compensates for the delay effect. However, the active damping inner loop can be equivalent to a second-order filter since a first-order HPF is in series with a first-order lead control part in it. us, the order of the system control is increased. e amplitude and phase characteristics of high-order filters vary more sharply and widely, making the parameters design more difficult [15].
In this paper, an improved SLICC is proposed. is method extends the valid damping region up to f s /2, improving the system stability margin greatly without increasing the control algorithm complexity and system cost. And it has inherent robustness against grid-impedance variation. e remainder of this paper is organized as follows. In Section 2, the control scheme and the mathematical model of a three-phase grid-connected inverter with an LCL filter are presented at first. e equivalent impedance model of the SLICC considering the control delay is established, and the detailed system stability analysis is given in Section 3. In Section 4, an improved SLICC for extending the valid damping region up to f s /2 is proposed. Moreover, a step-bystep parameter design procedure is presented. Section 5 shows the simulations and experimental results and Section 6 summarizes the conclusions. Figure 1 shows the circuit diagram and control scheme of a three-phase grid-connected inverter with an LCL filter and a constant DC voltage source U dc , where C is the filter capacitor and L 1 and L 2 are the filter inductance of the inverter side and the grid side, respectively. Considering the worstcase condition, where no passive damping of the resonance exists and the system becomes more unstable, the parasitic series resistances of its filter inductor and filter capacitor are all neglected here. u g is the grid voltage, and it has been assumed as balanced, which then allows per-phase diagrams to be used for analysis. i 1 and i 2 are the inverter and grid current, respectively. In this paper, to reduce the impact on the low-frequency stability of the system, the synchronous reference frame phase-locked loop (SRF-PLL) with low bandwidth has been employed to synchronize the inverter to the point of common coupling (PCC) voltage [20]. e overall control structure for per-phase based on SLICC is depicted in Figure 2. For simplicity, the equivalent gain of the inverter is set to 1. G PR (s) is the current controller, whose notation is given as follows [10]:

System Modeling and Description
where K p and K r are the proportional and resonance gains, respectively, ω i is the bandwidth of the resonance term, and ω 0 is the grid fundamental angular frequency. Generally, the digital signal processor (DSP) is introduced to control the grid-connected inverter, and the most common mode for PWM is the single-update mode. us, the computation and PWM delays should be considered. Specifically, G d1 (s) represents the computation delay, which is one sampling period in the traditional sampling scheme. G d2 (s) is the PWM delay, caused by the zero-order hold (ZOH) effect.
us, the total delay in the system can be expressed as [17] 2 Mathematical Problems in Engineering where T s is the sampling cycle of the digital control system. G LCL (s) is the transfer function of the LCL filter under SLICC strategy, as shown in the following equation:

Impedance-Based Analysis
is section generalizes the physical circuit property of the SLICC by presenting an impedance-based analysis [27]. For demonstrating circuit properties realized by SLICC, the control diagram is redrawn as in Figure 3(a) by sensing the voltage of the inverter-side inductance instead of the inverter current, and the equivalent circuit of LCL filter can be shown as Figure 3(b). is informs that the SLICC is equivalent to a virtual impedance Z v (s) connected in series with the inverter-side inductance. e expression of Z v (s) can be written as us, the equivalent virtual impedance Z v (s) can be shaped by varying G PR (s) and G d (s). Since the resonance term of G PR (s) mainly affects the frequency characteristics of the system at the fundamental frequency [19], the equivalent virtual impedance Z v (s) is simplified as K p e − 1.5sT s for simplicity.
If the control delay G d (s) is neglected, the equivalent virtual impedance is a constant positive resistance with a value of K p , which damps the LCL resonance peak. erefore, the SLICC can provide inherent damping to stabilize the system [18].
When the control delay G d (s) is incorporated, the frequency expression of the equivalent virtual impedance is yielded as It shows that the control delay changes the constant positive resistance to a virtual impedance which is composed of resistance and reactance, whose frequency characteristics are shown in Figure 4, and f s is the sampling frequency. e equivalent virtual resistance R v (ω) could be positive or negative depending on the frequency range, whose critical frequency f R is f s /6. Similarly, the equivalent virtual reactance X v (ω) could be inductive or capacitive and its critical frequency f X is f s /3. e negative R v (ω) will add open-loop right-half-plane (RHP) poles to the system, introducing a nonminimum-phase closed-loop behavior to the system, which can impair the overall system stability and robustness [20]. e existence of X v (ω) will shift the LCL resonance frequency f res [16]. As the damping effect, that is, stability of system, is determined by R v (ω) [22], the critical stable frequency of the traditional digital control system adopting SLICC is f s /6. When the actual LCL resonance frequency f res is up to f s /6, the system will be unstable due to the lack of necessary damping for the LCL resonance peak. erefore, the control delay G d (s) plays a decisive role in the stability of digital control system.
To identify the critical frequency f R of R v (ω) when the control delay varies, G d (s) can be rewritten as e − λsT s (normally 0 < λ ≤ 1.5). When R v (ω) is assumed to zero, it leads to where x is the ratio of the critical frequency f R to the sampling frequency f s . Figure 5 shows the relationship between x and λ according to (6). As shown in Figure 5, the critical frequency gradually increases as the control delay decreases. In particular, when the control delay decreases to e − 0.5sT s , the critical frequency reaches the upper limit of the LCL resonance frequency design value, that is, f s /2. at is, if the control delay can be reduced to e − 0.5sT s , the system stability can be ensured with the SLICC within the design range of the LCL resonance frequency.
Due to the inherent characteristics of the digital control system, reducing the control delay to e − 0.5sT s can only be realized by increasing the sampling frequency. In this case, the sampling frequency should be increased to three times its original value [12]. However, it may lead to problems that introduce the switch ripple into the control loop and multiple intersections between the modulated wave and the carrier [3]. at makes the method impractical.

Proposed Active Damping Method
Based on the above discussions, there is a need to decrease the control delay and compensate for its side effects. In this section, an improved active damping method based on SLICC is proposed, which includes the asymmetric regular Single-loop inverter current control scheme abc αβ Figure 1: ree-phase grid-connected inverter with an LCL filter. Mathematical Problems in Engineering sampling method and the delay compensation method. e proposed active damping method can extend the valid damping region up to f s /2, improving the system stability and robustness. e traditional SLICC has inherent robustness against grid-impedance variation, since the gridimpedance decreases the LCL resonance frequency [18]. e proposed active damping method can inherit this feature due to no additional active damping inner loop in the control loop.

Asymmetric Regular Sampling Method.
In the digital control system, in order to avoid the influence of switching noise, the signal sampling and update are operated at the peak or valley of the triangular carrier. e sampling signals are sent to the DSP; then the modulation waveform is obtained by calculating these signals with the control algorithm. Meanwhile, to avoid the phenomenon of multiple intersections between the modulation waveform and the triangular carrier, the calculated modulation waveform is updated at the next sampling instant, as shown in Figure 6(a). e computation delay is the duration between sampling instant of the signal and update instant of the corresponding calculated modulation waveform, which is equal to e − sT s . e PWM is often modeled as a ZOH, so the PWM delay is considered as half of the sampling period, which is equal to e − 0.5sT s . us, the total control delay in the control loop is e − 1.5sT s .
In order to directly reduce the control delay, the asymmetric regular sampling method is adopted, as shown in Figure 6(b). Instead of sampling and updating the signals one time in a control period, this method samples and updates twice time, which means that the signal sampling and update are operated at the peak and valley of the triangular carrier. Since the sampling frequency and update frequency are both doubled, the control delay is reduced to e − 0.75sT s , without changing the switching frequency and introducing the switching ripple. It is worth mentioning that, since the sampling frequency is doubled, this method requires a higher precision current sensor.

Analysis and Design of Delay Compensation Method.
As shown in Figure 5, when the control delay is e − 0.75sT s , the critical frequency of the equivalent virtual resistance is around f s /3. In this section, the delay compensation method is further introduced to compensate for the delay effect, realizing f R ≥ (f s /2) and ensuring that the system has a sufficient stability margin at f res .
Considering that the controller needs to be discretized, the simplest form of the delay compensator can be expressed as [22] where k � 1, 2, ..., m. Since the sampling frequency has doubled, the minimum delay in the control loop is e − 0.5sT s . en the equivalent virtual impedance can be expressed as where Figure 3: Equivalent control diagram and circuit of control diagram with SLICC strategy.  In order to ensure f R ≥ (f s /2), the equivalent virtual resistance should be positive within the interval of (0, (f s /2)]; thus, the constraints are derived as erefore, k is within the interval of [1,2). Figure 7 shows the frequency-domain characteristics of R v1 (ω) and X v1 (ω) with different k values. As can be seen in Figure 7, when k � 1, R v1 (ω) decreases as the frequency increases and decreases to 0 at f s /2. at means there is no damping at f s /2, and the more the LCL resonance frequency f res closer to f s /2, the smaller the stability margin. However, the LCL resonance frequency f res can be designed close to f s /2, and the change of the filter inductance L and the filter capacitor C may shift the actual LCL resonance frequency toward f s /2, worsening the system stability margin.
In contrast, when k � 1.5, R v1 (ω) stays constant within the interval of (0, f s /2], which means the system has a sufficient stability margin whatever the LCL resonance frequency f res is. However, k � 1.5 cannot be implemented in digital control system, and the only practical case would be k � 1. Nevertheless, the main goal of the above analysis is that there is a value that can satisfy the delay compensation requirements. As illustrated in [4,21,22], the delay compensation method works through compensating the lagging phase. us, we can mimic the delay compensator G C1 by the lead compensator. Figure 8 shows the Bode diagram of the delay compensator G C1 for different k (k � 1, 1.5). e leading phases at f s /2 are 45°and around 66°, respectively. Apparently, the leading phase at f s /2 when k � 1.5 is much larger than it when k � 1 that is why the system can be more stable when k � 1.5. To compare the performance of the lead compensator with the delay compensator, the leading phases at f s /2 are the same; that is, θ 1 � 45°and θ 2 � 66°. And the attenuation of the transmission coefficient of the open loop is fully compensated, as shown in Figure 8. e form of the compensator can be expressed as [12] where α and T are us, the equivalent virtual impedance of the system can be expressed as where     Figure 9(a) shows the frequency-domain characteristics of R v1 (ω) and R v2 (ω); Figure 9(b) shows the frequencydomain characteristics of X v1 (ω) and X v2 (ω). As can be seen from Figure 9(a), through the same leading phase at f s /2, R v2 (ω)| θ�45°i s always greater than R v1 (ω)| k�1 within the interval of (0, f s /2). is is because the lead compensator can provide a bigger leading phase. What is more, R v2 (ω)| θ�45°i s greater than R v1 (ω)| k�1.5 within the interval of (0, 0.4f s ), and the range can be extended to (0, f s /2) by increasing the leading phase. us, the system can achieve more effective damping when the lead compensator is introduced, and the parameters design procedure can be more flexible.
It is worth mentioning that X v2 (ω)| θ�45°i s always greater than X v1 (ω)| k�1 within the interval of (0, f s /2), even greater than X v1 (ω)| k�1.5 , as shown in Figure 9(b). And X v2 (ω) can be greater within high-frequency interval by increasing the leading phase θ. is is because of the amplitude amplification effect of the lead compensator, resulting in shifting the actual resonance frequency f res toward f s /2, decreasing the stability margin of the system. us, the leading phase θ should be designed carefully.

Design Guideline for the Proposed Active Damping
Method. In this section, a step-by-step parameter design procedure considering the phase margin (PM) and the gain margin (GM) of the system is summarized as follows: Step 1. Consider the requirements of the power level and the current limit of the inverter; the switching harmonic suppression ability, the filter inductances L 1 and L 2 , and filter capacitance C are preliminarily determined [1].
Step 2. Set the cutoff frequency f c according to the sampling frequency f s , low-frequency harmonic suppression ability, and overshoot constraint. As reported in [10], to achieve an overshoot within the reasonable range, f 0 < f c < 4%f s is expected with PM 1 at f c over 60°, where f 0 is the fundamental grid frequency.
Step 3. Specifically, K p is related to f c by [5] When PR controller is adopted, the magnitude of the loop gain at f 0 should be bigger than 75 dB, ensuring that the amplitude error is less 1% at the rated power [5]. Once f c is specified, K p and K r can be determined.
Step 4. Draw the Bode diagram of the open-loop transfer function of the overall system with the asymmetric regular sampling method only. Generally, the PM 2 at f res is expected within the interval of [30°, 60°], and GM > 5 dB is expected [5]. en examine the actual PM 2 and GM. If the actual PM 2 or GM is not satisfied, go back to Step 3 and slightly decrease K p and K r .
Step 5. If the actual PM 2 or GM cannot be satisfied only by decreasing K p and K r , introduce the delay compensator and examine the actual PM 2 and GM. Go back to Step 3 and slightly decrease K p and K r if necessary.
Step 6. If the actual PM 2 or GM still cannot be satisfied, introduce the lead compensator and set the leading phase θ at f s /2 to be � π/4. Go back to Step 3 and slightly decrease K p and K r if necessary. If it is not satisfied again, increase the leading phase according to the difference between the actual PM 2 and the preset value, and further slightly decrease K p and K r . If it is still not satisfied, go back to Step 1 and slightly decrease the LCL resonance frequency and repeat the other steps until the system stability requirements are met. e flowchart can be shown in Figure 10. e parameters are designed according to the design procedure, as shown in Table 1. Figure 11 shows Bode diagram of the overall system with the designed parameters, and the stability requirements are satisfied, since PM 1 � 70°, PM 2 � 31°, and GM � 6.5 dB.

Simulation Results.
To verify the above analysis and the proposed active damping method, a simulation model based on Figure 1 is built in the MATLAB/SIMULINK, and the main system parameters are shown in Table 1.
For verifying the analysis of the critical frequency of the SLICC, the LCL resonance frequency f res is selected as 0.24f s by intentionally increasing the filter capacitance C to 9.4 μF. Figure 12 shows the simulated grid current for this case. As shown in Figure 12(a), when f res is bigger than the critical frequency of the traditional SLICC, that is, f s /6, the system cannot be stabilized. After decreasing the control delay by adopting the asymmetric regular sampling method, the system can be stabilized while the grid current still Phase (deg) Figure 8: Bode diagram of G C1 and G C2 with different leading phases. 6 Mathematical Problems in Engineering Frequency (Hz) Figure 9: Frequency-domain characteristics of the equivalent virtual impedance with different leading phases. (a) R v1 (ω) and R v2 (ω); (b) X v1 (ω) and X v2 (ω).
Determine the filter inductance and filter capacitance based on the requirements of the power level, the current limit of the inverter, and the switching harmonic suppression ability Set the cutoff frequency according to sampling frequency, low-frequency harmonic suppression ability, and overshoot constraint

Parameters determined
The compensation effect of delay compensator has already been tested?
Decrease the LCL resonance frequency by changing the filter inductance and the filter capacitance   Mathematical Problems in Engineering contains oscillation, indicating that the system stability margin is insufficient. When further adopting the delay compensator, the current oscillation rapidly decays, indicating that the system stability margin has been increased, as shown in Figure 12(b). It is worth mentioning that, under this condition, there is little difference between the results of the delay compensator and the lead compensator (not shown in Figure 12). at means when the LCL resonance frequency f res is relatively low, the system stability margin under the two methods has little difference, which is consistent with the analysis in Section 4.
When the filter capacitance C is regularly designed according to the design guidance, it equals 4.7 μF, as shown in Table 1. In this case, f res � 0.34 f s , which is close to the critical frequency when the control delay decreases to e − 0.75sT s , that is, around f s /3. Figure 13 shows the simulated grid current for this case. With adopting the asymmetric regular sampling method, the control delay reduces to e − 0.75sT s , and the system still cannot be stabilized, as shown in Figures 13(a) and 13(b). e system can be stabilized when further adopting the delay compensator, as shown in Figure 13(a). However, the current still contains high resonance frequency harmonics, which makes the current cannot meet the requirements of total harmonic distortion (THD). Comparatively, when adopting the lead compensator, there is no obvious oscillation in the steady state, as shown in Figure 13(b). It indicates that, at a high level of the LCL resonance frequency, the system has a higher stability margin when adopting the lead compensator.
In addition, the robustness of the proposed active damping method is verified by shifting the LCL resonance frequency. Figure 14 shows the contrastive simulation results under the proposed method with the delay compensator and the lead compensator. e LCL resonance frequency reaches 0.4 f s by decreasing the filter capacitor C to 75% designed value. e resonance frequency harmonics increase significantly with the delay compensator, as shown in Figure 14(a). Comparatively, with the implementation of the lead compensator, the resonance frequency harmonics increase slightly, as shown in Figure 14(b). Compared to the results shown in Figure 13, it indicates that the higher the resonance frequency, the greater the advantage of the lead compensator, verifying the effectiveness of the analysis and the proposed method.

Experiment
Results. In order to verify the above simulation results, a 6 kW three-phase grid-connected inverter prototype is built in the laboratory, as shown in Figure 15.
In practice, the PR controller is decomposed into a direct integrator and a feedback integrator, and they are discretized with a forward difference and backward difference separately [10]. e delay compensator and the lead compensator are discretized with Tustin transform with prewarping to match the frequency responses at f res [5].
us, the discretized controller can be implemented in a TI TMS320F28335 DSP. e main system parameters are shown in Table 1. Figure 16(a) shows the grid current with the asymmetric regular sampling method adopted only.
ough the grid current contains large oscillation, the system can be stabilized when the actual LCL resonance frequency f res is close to the critical frequency. at is because of the damping effects of the parasitic series resistances of the filter inductor and filter capacitor in the prototype, leading to the difference between Figures 13 and 16(a). Figures 16(b) and 16(c) show the grid current with the delay compensator and the lead compensator, respectively. e system stability margin is improved after introducing the delay compensator, and the oscillation rapidly decays, as shown in  Figure 16(b). But the resonance frequency harmonics still cannot meet the standards of THD. Comparatively, when adopting the lead compensator, there is no obvious oscillation in the steady state, as shown in Figure 16(c). It indicates that when the system has a high level of LCL resonance frequency, a higher stability margin can be achieved by adopting the lead compensator, verifying the effectiveness of the proposed method and the simulation results.

Conclusions
is paper studies active damping method for LCL-filtered three-phase grid-connected inverters based on SLICC.
rough systematic impedance-based analyses, it has been shown that the control delay has changed the equivalent virtual resistance of the SLICC, and the critical frequency of the equivalent virtual resistance locates within the range of the LCL resonance frequency. To ensure the stability of the SLICC, an improved active damping method is proposed. e damping region is widened to (0, f s /2), which can cover all the possible resonance frequencies. A step-by-step design procedure considering the different resonance frequency range is also formulated. Particularly, the delay compensator may be unable to provide sufficient PM within the interval of high resonance frequency, and the lead compensator could further improve the system stability. e simulation and experimental results show the effectiveness of the proposed active damping method. To mimic the delay compensator more precisely when the control delay is equal to e − 0.75sT s is the future work.
Data Availability e data that support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest regarding the publication of this paper.