A Systematic Control Design Method with Active Damping Control in Voltage Source Converters

: This paper proposes a systematic control design method for active damping control of grid-connected voltage source converters (VSCs). The proposed control method considers the conventional cascaded control loops and improves them by including additional states feedback-based active damping. In such a way, all control gains are lumped into one control gain matrix based on the proposed formulation. The lumping of all control gains into one matrix leads to a linear optimization problem, so different techniques can be used to calculate control gains. This work calculates them by using a simple but effective optimal control theorem as a noteworthy feature. The proposed control method can overcome the challenges of designing multiple control loops, evaluating wide time scale dynamics, and tuning required control parameters. Moreover, direct relationships between the proposed tuning parameters and system well-known stability and performance indicators such as maximum damping factor, minimum damping ratio, and the control efforts are identiﬁed, providing good physical insight. Finally, the proposed control structure and optimal gain calculations ensure power converter robustness against uncertainties in the grid’s short-circuit ratio (SCR) and different operating-point conditions. When the grid’s SCR changes from 10 (strong grid condition) to 1 (ultra-weak grid condition), the system under the proposed control method maintains good stability margins and simultaneously provides a fast dynamic response by facilitating the implementation of a high-bandwidth phase-locked loop (PLL). The performance of the proposed control strategy was investigated analytically and practically by conducting eigenvalue analysis, simulations, and experiments.


Introduction
Grid-connected voltage source converters (VSCs) have become one of the main components of the modern power system due to their fast dynamic response, full controllability, and flexibility.They can be found in many power-system applications such as renewableenergy-based generation (wind turbines (WT) and photovoltaic (PV) power plants), electric power transmission systems, and load and power quality equipment.Multiple control loops are employed in the control of VSCs to improve the performance, including active power control loop or dc-link voltage controller (DVC), reactive power control loop or ac voltage magnitude controller (AVC), inverter current control loop (CC) and protection, voltage feedforward controller (VFF) to damp LCL filter resonances, and a phase-locked loop (PLL) for grid synchronization.However, it has been revealed that these control loops are strongly coupled and lead to complex, asymmetric, and wide timescale dynamics.Therefore, an inaccurate control loop design and the possible interaction among different control loops may cause performance degradation or even system instability [1][2][3][4].These stability concerns become more severe under weak grids, where the power grid at the connection point has a lower short circuit ratio (SCR) due to the high penetration level of power converters or long distance between distributed generators and the main grid.
Many efforts have been made to present more efficient design strategies and enhance the control system's performance by including additional active damping terms to address the mentioned problems.Conventional design strategies such as the generalized Nyquist theorem or passivity theorem are based on the impedance model to investigate the system's stability and control-parameter selection [5,6].For instance, in Reference [7], a recursive adjustment procedure for the current controller's gains is suggested based on the frequency response.This work considers only the control delay effect.A comprehensive stability analysis and the current control design of a grid-connected LCL-filtered VSC are discussed in Reference [8].To damp the LCL filter resonances and the negative impact of the control delay, a grid-side current control-based active damping and virtual flux-based active damping are proposed in References [9,10], respectively.Moreover, to limit the negative impact of digital filters (commonly used in the CC and VFF) on the dissipativity of the output converter admittance, some design guidelines are proposed by Reference [11].However, calculating the impedance model and performing a stability analysis are challenging and require high expertise and physical insight.Moreover, designing a control system and selecting the optimal values of the control parameters are not trivial tasks, especially for a high-order multi-input, multi-output (MIMO) system such as grid-connected VSCs.Therefore, the control gain selection is generally an iterative process.
A full state-space feedback control method to improve the resonance damping of the LCL filter and increase the system robustness against the grid SCR variations is presented in References [12,13].The proportional-integral (PI) current control and active damping gains in Reference [12] are calculated by using the pole placement technique.Despite the simplicity in concept, selecting the desired closed-loop poles is a challenging task; it becomes more difficult or even impossible for a high-order MIMO system.A systematic design procedure based on the optimal control theorem is suggested in Reference [13] to avoid direct pole placement.However, weighting factors' selection is required [13,14].
In References [15][16][17][18][19], a robust proportional-resonant (PR) current control design based on the H-infinity theorem and linear matrix inequalities (LMIs) is proposed to consider the impact of the grid SCR variations in the design procedure.However, the optimal robust H-infinity control leads to high control gains, making some difficulties for practical implementation [15,16,18].
Although the previous works provide valuable results, they did not consider the impact of PLL and outer control loops.Therefore, they suffer from some conservatism in their design.Moreover, active damping based on the feedforward PCC voltage or current feedback, which is widely suggested, is only efficient to damp oscillations in the LCL filter resonance range and not low-frequency oscillations due to the outer control loops interaction.
In Reference [20], a state feedback-based active damping control to dampen lowfrequency oscillations is proposed.However, a low-order and simplified model is considered for the fast dynamics parts (LCL filter, delay, and current controller) to simplify the control system design.Moreover, it suffers from shortcomings of the pole-placement techniques and recursive process.Finally, the possibility of implementing a high-performance PLL in weak grids was not discussed.
Research work in Reference [21] introduced an improved current control design to reduce the negative impact of the PLL on the current controller in weak grids, which need extra capacitor current sensors.
In Reference [22], an optimal control method is proposed to enhance the inner current control loop performance and reduce the negative impact of PLL in weak grids.However, many weighting factors are added to the design procedure, transparent relationships between tuning parameters and stability and performance indicators (e.g., damping factor and ratio) are not identified, and the outer control loops are not considered.
The authors of Reference [23] presented a loop-at-a-time stability analysis (LAAT), and based on LAAT, a design algorithm for outer control loops (DVC and AVC) is also suggested.The proposed method is closely correlated to the characteristic loci in the generalized Nyquist criterion for MIMO systems and tries to overcome challenges regarding the MIMO control system design in the frequency domain.However, to reduce the complexity of the method, the inner current control loop and PLL are embedded inside the plant.
In Reference [24], a fixed-structure H-infinity control design is presented to calculate the outer loops' control gains.However, the design procedure needs dynamic weighting function calculations for each input and output relation.In this regard, many parameters are required to be tuned, significantly affecting the system response and feasibility of the H-infinity optimization problem.
Overall, the prior state-of-the-art methods mentioned above partially address grid uncertainties, low and high-frequency oscillations, control loops interaction, high-order MIMO structure, and recursive control system design.However, none of them succeeds in presenting a comprehensive and efficient solution, as summarized in Table 1.For example, all previous works consider only inner or outer control loops and ignore the other ones in the design step.When all control loops are considered simultaneously, the couplings of some inner and outer control parameters have appeared, leading to a nonlinear system and optimization problem.In this respect, this paper tries to propose a new systematic design approach in order to simplify the control gains calculation by providing a linear optimization problem.Based on the new formulation, all control gains are lumped into one control gain matrix, and then different techniques can be used to calculate the system control gain matrix.Among others, this work employs the linear-quadratic regulator (LQR) design framework due to its efficiency and applicability.
The main contributions of this paper are listed as follows: • Proposing a direct design approach with the minimum recursive process to calculate the control parameters,

•
Providing a clear insight between tuning parameters and stability and performance indicators,

•
Increasing system robustness against uncertainties in the grid impedance and operating point conditions,

•
Improving system transient response and fault ride-through (FRT) capability by implementing a PLL with higher bandwidths.
The remainder of this paper is organized as follows.Section 2 describes the understudied system, a grid-connected LCL-filtered VSC with all control loops.The proposed nonlinear and linearized dynamics of a grid-connected VSC are also presented in this section.The proposed robust and optimal active damping control is introduced in Section 3. Furthermore, some guidelines for selecting the state and input weighting matrixes are provided.This section is completed with a comprehensive robustness analysis based on the eigenvalue analysis to consider the impact of different grid SCRs and operating point conditions.A comparison with the conventional one is also given.The analytical results were verified by several simulations and experimental tests under different operational conditions, as recounted in Section 4. Finally, Section 5 concludes this paper.

System Description and Modeling
Figure 1 shows the principal scheme of the studied system, which includes primary energy sources (i.e., the solar power plant or the wind turbine systems), the power grid, back-to-back three-phase voltage source converter (VSC), an LCL output filter, and a control system.The power grid is represented by a Thevenin model seen from the point of common coupling (PCC).

•
Proposing a direct design approach with the minimum recursive process to calculate the control parameters,

•
Providing a clear insight between tuning parameters and stability and performance indicators,

•
Increasing system robustness against uncertainties in the grid impedance and oper ating point conditions,

•
Improving system transient response and fault ride-through (FRT) capability by im plementing a PLL with higher bandwidths.
The remainder of this paper is organized as follows.Section 2 describes the under studied system, a grid-connected LCL-filtered VSC with all control loops.The proposed nonlinear and linearized dynamics of a grid-connected VSC are also presented in this sec tion.The proposed robust and optimal active damping control is introduced in Section 3 Furthermore, some guidelines for selecting the state and input weighting matrixes are provided.This section is completed with a comprehensive robustness analysis based on the eigenvalue analysis to consider the impact of different grid SCRs and operating poin conditions.A comparison with the conventional one is also given.The analytical results were verified by several simulations and experimental tests under different operationa conditions, as recounted in Section 4. Finally, Section 5 concludes this paper.

System Description and Modeling
Figure 1 shows the principal scheme of the studied system, which includes primary energy sources (i.e., the solar power plant or the wind turbine systems), the power grid back-to-back three-phase voltage source converter (VSC), an LCL output filter, and a con trol system.The power grid is represented by a Thevenin model seen from the point o common coupling (PCC).The grid inductance (Lg2) can also vary widely to mimic a weaker or stronger connec tion point.The remaining parameters in Figure 1 are as follows: if, ig, vf, v, vg, vdc, and idc are the converter and grid currents, voltage across the filter capacitor, inverter and grid volt ages, dc-link voltage, and dc-side current, respectively.Moreover, Lf, Lg1, Cf, Lg2, and Cd represent the converter and grid-side filter inductances, filter capacitance, equivalent grid inductance, and dc-link capacitance, respectively; δ is the phase angle obtained by the PLL; and the subscript ref gives the reference values for different control loops.

Derivation of Proposed Control Method
The main aims of the grid-side VSC are to inject the produced power by renewable energy systems into the power grid with the sinusoidal current.However, the voltage compensation of the PCC is also necessary when connected to a weak grid.The grid inductance (L g2 ) can also vary widely to mimic a weaker or stronger connection point.The remaining parameters in Figure 1 are as follows: i f , i g , v f , v, v g , v dc , and i dc are the converter and grid currents, voltage across the filter capacitor, inverter and grid voltages, dc-link voltage, and dc-side current, respectively.Moreover, L f , L g1 , C f , L g2 , and C dc represent the converter and grid-side filter inductances, filter capacitance, equivalent grid inductance, and dc-link capacitance, respectively; δ is the phase angle obtained by the PLL; and the subscript ref gives the reference values for different control loops.

Derivation of Proposed Control Method
The main aims of the grid-side VSC are to inject the produced power by renewable energy systems into the power grid with the sinusoidal current.However, the voltage compensation of the PCC is also necessary when connected to a weak grid.
The conventional control systems of grid-side VSC are shown in Figure 2. The conventional one consists of five main parts, including the dc-link voltage controller (DVC), AC voltage magnitude controller (AVC), current controller (CC), voltage feedforward controller (VFF), and a synchronous reference frame phase-locked loop (SRF-PLL).It is worth remarking that they all (except for VFF) employ a proportional-integral (PI) controller to ensure good dynamic response and zero steady-state tracking error in the synchronous reference frame.Moreover, grid-side current feedback and voltage feedforward control through a high-pass filter are provided to damp the LCL filter resonances and better disturbance rejection.In this work, a first-order high-pass filter is considered for VFF, which is expressed as follows [25]: AC voltage magnitude controller (AVC), current controller (CC), voltage feedforward controller (VFF), and a synchronous reference frame phase-locked loop (SRF-PLL).It is worth remarking that they all (except for VFF) employ a proportional-integral (PI) controller to ensure good dynamic response and zero steady-state tracking error in the synchronous reference frame.Moreover, grid-side current feedback and voltage feedforward control through a high-pass filter are provided to damp the LCL filter resonances and better disturbance rejection.In this work, a first-order high-pass filter is considered for VFF, which is expressed as follows [25]: As it can be seen, the conventional structure employs many control loops with at least eleven control gains which should be designed carefully to provide a proper response and stability margins and also avoid control loops' interactions.Moreover, there are many states in the system that can be appropriately utilized.It seems proper that employing them may remarkably increase the system's robustness and performance, as shown in Figure 3. Here, Kf and x are the active damping gain matrix and system state's vector, respectively.However, it may add extra feedback gains and increase complexity to calculate the optimal control gains.
The control structure in Figure 3 can be simplified as follows.Since kpd and Kf1(1,15) are multiplied by Δvdc and provide the same damping, one can be ignored to reduce the number of tuning parameters.The same can also be concluded for CC, AVC, and integral gains.The resultant simplified control structure is depicted in Figure 4.
By comparing Figures 3 and 4, it can be seen that duplicate coefficients are ignored, and all control gains of different control loops and state feedback-based active damping are lumped in the control gain matrix (Kf).
Therefore, the proposed control method can utilize the maximum capability of the circuit and control systems and employs all available states (x) to improve the system's response.As it can be seen, the conventional structure employs many control loops with at least eleven control gains which should be designed carefully to provide a proper response and stability margins and also avoid control loops' interactions.Moreover, there are many states in the system that can be appropriately utilized.It seems proper that employing them may remarkably increase the system's robustness and performance, as shown in Figure 3. Here, K f and x are the active damping gain matrix and system state's vector, respectively.However, it may add extra feedback gains and increase complexity to calculate the optimal control gains.
The control structure in Figure 3 can be simplified as follows.Since k pd and K f1 (1,15) are multiplied by ∆v dc and provide the same damping, one can be ignored to reduce the number of tuning parameters.The same can also be concluded for CC, AVC, and integral gains.The resultant simplified control structure is depicted in Figure 4.
By comparing Figures 3 and 4, it can be seen that duplicate coefficients are ignored, and all control gains of different control loops and state feedback-based active damping are lumped in the control gain matrix (K f ).
Therefore, the proposed control method can utilize the maximum capability of the circuit and control systems and employs all available states (x) to improve the system's response.
Moreover, a high-pass filter-based VFF, which is necessary to damp LCL filter responses under weak grid conditions, is also replaced with the simpler constant feedback controller.It is shown here that, although a dynamic controller is not considered for the VFF, the proposed control structure can still provide better robustness under weak grid conditions, as another noteworthy feature.
The remaining question is how the control gain matrix can be obtained and based on which criteria.In the following, at first, a proper representation of the system equations is provided.Then, based on the proposed system model, a design procedure is presented to obtain the control gain matrix optimally and simply.
VFF, the proposed control structure can still provide better robustness under weak conditions, as another noteworthy feature.
The remaining question is how the control gain matrix can be obtained and base which criteria.In the following, at first, a proper representation of the system equatio provided.Then, based on the proposed system model, a design procedure is presente obtain the control gain matrix optimally and simply.VFF, the proposed control structure can still provide better robustness under weak conditions, as another noteworthy feature.The remaining question is how the control gain matrix can be obtained and base which criteria.In the following, at first, a proper representation of the system equatio provided.Then, based on the proposed system model, a design procedure is presente obtain the control gain matrix optimally and simply.

Proposed Nonlinear and Linearized State-Space Model of a Grid-Connected VSC for Optimal Control System Design
Considering the active power-balance equation, Park transformation, and control delay in the digital implementation in Figures 1 and 4, the proposed state-space nonlinear model of the overall system is determined as follows: .
where L g = L g1 + L g2 , and γ q , γ dc , γ g , γ id , and γ iq are integral states for the PLL, DVC, AVC, and current controller (CC), respectively.It is worth noting that the converter is assumed to be lossless, and the delay is approximated by the following transfer function, using first-order Pade approximation [26]: Furthermore, two dq-frames are considered to include the impact of PLL in the system model and analysis.The global dq-frame is aligned with the infinite bus voltage, and the converter or local dq-frame is aligned with the PCC voltage.Signals in the global and the converter dq frames are represented with and without the apostrophe superscript ('), respectively [27].
The nonlinear model can be rewritten in the following compact form: .
Here the system state (x) and the control and disturbance inputs (u and d) are as follows: A linearized state-space model around an equilibrium point (x e , u e , and d e ) is required, which facilitates designing a state feedback control for a grid-connected VSC.The following equations can be used to calculate the linearized model:

Design of the Control Gain Matrix
In the existing literature [12][13][14][15][16][17][18][19], it is evident that using more states in the control structure is better for the dynamic performance and stability margin under a weak grid.However, how to design the control gain matrix is an unsolved question.This section presents a linear-quadratic regulator (LQR) design framework to achieve an optimal, robust, and multivariable control of grid-connected VSCs.
Given the system dynamics in ( 6), the static full-state-feedback law is as follows: The optimal gain matrix, K f , is calculated to minimize the following infinite-horizon quadratic cost function: where Q and R are constant weighting matrixes that set relative weights of state deviation and input usage; the optimal gain matrix is given by the following: where the positive definite matrix, P, is a solution of the following algebraic Riccati equation: The solution leads to a stable closed-loop transfer function A cl = A − BK f and minimizes the defined performance quadratic cost function (J).Furthermore, it is wellknown that the equation has a single positive-definite solution if, and only if, the pairs of (A, B) and √ Q, A are stabilizable and detectable, respectively.Moreover, there are suitable numerical methods to find the solution, even for MIMO systems or high -dimensional problems.

Weighting Matrixes' Selection
As the previous subsection can conclude, the optimal state-feedback control design of grid-connected VSC is reduced to the proper choice of weighting matrixes (Q and R).Bryson's rule can be used to select weighting matrixes as a primary choice.Moreover, References [28,29] have provided some guidelines to choose weighting matrixes.For instance, a procedure based on minimizing the infinity norm of selected transfer functions is suggested in Reference [29].However, there are two main issues for the understudied system that should be appropriately answered: (1) Is it possible to reduce the number of tuning parameters?(Q and R have 130 elements to tune.) (2) Can physical meaning be provided for them, and how can they be related to well-known system stability and performance indicators such as the maximum damping factor (σ max ) and minimum damping ratio (ξ min )?
Many recommendations for properly tuning weighting matrixes in general and, after that, for the specific understudied system are presented in this section.For simplicity in tuning, diagonal matrixes are usually selected.It is worth remarking that the diagonal elements of weighting matrixes (q i or r i ) are the maximum expected or acceptable value of the related state x i 2 −1 or input u i 2 −1 .A higher Q results in a smaller steady-state error, higher bandwidth, faster dynamic response, and better disturbance rejection.However, it causes higher sensitivity to noise and unmodeled dynamics and increases the norm of control gain matrix, leading to chattering and control input saturation.Therefore, the selection of Q and R compromises a fast dynamic response, adequate stability margin, and acceptable control matrix norm.The compromised values can be obtained systematically with a simple computer-aided program.For example, for different values of weighting matrixes, the damping factor and damping ratio of the critical mode and the norm of the control gain matrix (norm(K f ) = K f ) are calculated.Then the stability and performance criteria are checked.If the desired performance is fulfilled, the process is stopped; otherwise, proper actions should be taken.
Although the previous process can also find proper weighting factors, nineteen tuning parameters (fifteen parameters in Q and four in the R) should be selected, which increases the design process time.Therefore, further reduction of the number of individual tuning parameters is still desirable.Moreover, physical inspiration about their impact is not presented yet.Thus, a more detailed discussion on the weighting factors tuning for the understudied system is presented as follows.
Since there are four control inputs (u 1d , u 1q , u 2d , and u 2q ) and three important stability and performance indicators (σ max , ξ min , and norm (K f )), three or four individual parameters in the weighting matrixes may provide the requirements mentioned above.
Therefore, the following suggestions are provided.First, the input cost matrix is selected as R = rI 4×4 ; therefore, only one r is needed to be chosen.Moreover, the different characteristics of the states can be divided into three different groups: integrals of the tracking errors γ id , γ iq , γ dc , γ g , γ q , dc-link voltage tracking error (∆v dc ), and the remaining states x dd , x dq , δ, i f d , i f q , v f d , v f q , i gd , i gq .
For the first group, these states are directly related to the control errors, which represent the steady-state error.For the second group, it can affect the transient response of dc-link voltage tracking error and, consequently, inverter current.For the last one, the states are just intermediate states which are directly related to the control performance.Based on the classification of the states, the states' cost matrix is selected as Q = diag(q 1 , q 1 , q 2 , q 2 , q 1 , q 1 , q 1 , q 2 , q 2 , q 2 , q 2 , q 2 , q 2 , q 3 ).
First, stability and performance definitions are provided.Then the impact of different values of tuning parameters (q 1 , q 2 , q 3 , and r) on the system stability and performance and the control gain matrix are studied.The results are summarized in Table 2.
Table 2. Impact of tuning parameters (q 1 , q 2 , q 3 , and r) on the control gain matrix and the system stability and performance.The following equations can represent the ith eigenvalue of the system (λ i ) and its damping factor (σ i ), frequency of oscillations (ω di ), and damping ratio (ξ i ): It is worth remarking that the system is stable for σ < 0 and ξ > 0, oscillating for σ = ξ = 0, and unstable for σ > 0 and ξ < 0. The unstable or weakest damped eigenvalue among all eigenvalues, called a critical mode, causes significant oscillations or even system instability.To guarantee stability margin and limit the system overshoot, the damping factor and ratio of the critical mode should fulfill some requirements (σ ≤ σ & min max ).The de- sired damping ratio depends on the modal frequency and accepted oscillation settling time by different utilities and standards [30,31].It is usually selected to be something between 5 and 20 percent.In this work, a desired 10% damping ratio is accepted, guaranteeing an acceptable overshoot and decay rate of oscillations.Moreover, placing the maximum damping factor far from the imaginary axis can ensure system stability and improve response speed in terms of the rise time.The system's rise time is assigned by the outer voltage control loops that have lower bandwidths among different control loops.For example, a rise time (t r ) lower than 100 ms is usually required for AVC, resulting in a maximum damping factor less than 22 rad/s: (σ max ≈ −2.2/t r = −2.2/100ms = −22 rad/s) [23].

Selection of q 1 :
As it can be seen in Table 2, this parameter has a significant impact on the maximum damping factor.For example, for q 1 = 10, the system is close to the imaginary axis and hazardous area (σ max = −1.5 rad/s).By increasing this parameter to q 1 = 1 × 10 4 , the critical mode moves far from the imaginary axis to obtain a good stability margin (σ max = −40 rad/s).However, by increasing the q 1 more, the critical mode is not changed anymore, while the control gain norm increases.Thus, at this point, if the obtained damping factor is acceptable, q 1 can be used for the control gain calculation in (10).Otherwise, the new initial values for all tuning parameters should be selected, and the process should be repeated.
A good decision is to set a lower value for r to allow more control input usage.This action may provide a better stability margin in terms of the damping factor.Finally, as it can be seen, the damping ratio is not affected, and it remained constant for different values of q 1 .Selection of q 2 : Unlike the q 1 , q 2 does not affect the damping factor, while it affects the damping ratio considerably.For a lower value, the ξ min is close to zero, which is undesirable and causes higher overshoots.Increasing this parameter increases the system damping ratio remarkably in the cost of higher control gains.Therefore, the selection of this parameter is a compromise between damping ratio and control gains norm.
Selection of q 3 : Overall, q 3 has a significant impact on the tracking error of dc-link voltage.A higher q 3 causes the faster tracking of dc-link voltage at the price of higher overshoot in the inverter current.The inverter current (especially d-axis current) acts as a control input for DVC.Therefore, it can be expected that the faster convergence rate of the dc-link voltage may lead to more current variations and overshoots.Inverter overcurrent protection limits the current overshoot.Thus, there is also a compromise here that q 3 can obtain.
Figure 5 shows system response for two different values of q 3 .For a higher value, dc-link voltage tracks the reference one quickly and without oscillations at the expense of the higher current overshoot and vice versa.

Selection of r:
This parameter affects the norm of the control gain matrix.Lower values mean permission for higher control gains and vice versa.As expected, and as seen in Table 2, for the lower value of r, the stability and performance indicators are superior, but the norm of the control gain matrix is high.High control gains cause many difficulties for practical implementation and increase the sensitivity to noise and time delay.Figure 6 shows the system output for two different norm values of the control gain matrix.For the higher one, chattering and oscillations in outputs are seen.By increasing r and reducing the norm of the control gain matrix, these oscillations disappear.A summary of the weighting factors design procedure: As shown so far, in the proposed design algorithm, an almost direct relation between tunning parameters and desired stability and performance requirements is identified that Figure 5. Impacts of tuning parameter q 3 on the system outputs (q 1 = 1 × 10 4 , q 2 = 1, and r = 1.5), (q 3 = 1, norm (K) = 9.7 × 10 3 , σ max = −39, and ξ min = 0.28) in blue, and (q 3 = 10, norm (K) = 9.9 × 10 3 , σ max = −31, and ξ min = 0.28) in orange.

Selection of r:
This parameter affects the norm of the control gain matrix.Lower values mean permission for higher control gains and vice versa.As expected, and as seen in Table 2, for the lower value of r, the stability and performance indicators are superior, but the norm of the control gain matrix is high.High control gains cause many difficulties for practical implementation and increase the sensitivity to noise and time delay.Figure 6 shows the system output for two different norm values of the control gain matrix.For the higher one, chattering and oscillations in outputs are seen.By increasing r and reducing the norm of the control gain matrix, these oscillations disappear.

Selection of r:
This parameter affects the norm of the control gain matrix.Lower values mean permission for higher control gains and vice versa.As expected, and as seen in Table 2, for the lower value of r, the stability and performance indicators are superior, but the norm of the control gain matrix is high.High control gains cause many difficulties for practical implementation and increase the sensitivity to noise and time delay.Figure 6 shows the system output for two different norm values of the control gain matrix.For the higher one, chattering and oscillations in outputs are seen.By increasing r and reducing the norm of the control gain matrix, these oscillations disappear.A summary of the weighting factors design procedure: As shown so far, in the proposed design algorithm, an almost direct relation between tunning parameters and desired stability and performance requirements is identified that considerably simplifies the control system's design.Figure 6.Impacts of tuning parameter r and norm of control gain matrix (K) on the system outputs (q 1 = 1 × 10 4 , q 2 = 1, and q 3 = 5): (a) (r = 0.5, norm (K) = 1.7 × 10 4 , σ max = −41, and ξ min = 0.4) and (b) (r = 1.5, norm (K) = 9.1 × 10 3 , σ max = −41, and ξ min = 0.28).

A summary of the weighting factors design procedure:
As shown so far, in the proposed design algorithm, an almost direct relation between tunning parameters and desired stability and performance requirements is identified that considerably simplifies the control system's design.
In such a way, by adjusting q 1 , q 2 , r, and q 3 , the desired maximum damping factor, minimum damping ratio, control gain matrix norm, and transient response are obtained.
Hence, the design process can be summarized as follows: Step 1: Small values for q 1 , q 2 , q 3 , and r are selected.
Step 2: Increasing q 1 adjusts the maximum damping factor.
Step 3: Increasing q 2 adjusts the minimum damping ratio.
Step 4: Increasing r reduces control gain matrix norm.
Step 5: Increasing q 3 changes dc-link voltage and inverter current transient response.
As it can be seen, the design procedure is very straightforward, and repetitive actions are reduced considerably.

Robustness Analysis
The previous design process can provide sufficient stability margins.However, these excellent stability margins are guaranteed for the nominal values of system parameters.However, some of these parameters may not be constant during the inverter operation.For example, their operating-point conditions continuously change during the day, and the grid SCR can vary over a wide range due to variations of the connected load and generators.
In this respect, a robustness assessment based on the eigenvalue calculation was carried out to examine the system's robustness under these uncertainties.In such a way, the variations in the critical mode are studied for different values of grid SCR and operating point conditions.In addition, a comparison between the proposed and the conventional control methods is also given.
The stability margins of the proposed control method are calculated based on the linearized model in (6) and feedback control gain in (10) A cl = A − BK f .It is worth remarking that the control gain matrix in (10) is calculated for the nominal power (P n ), the average value of grid SCR (SCR = 5.5), and f PLL = 10 Hz.Then the control gain matrix is kept constant during the robustness assessment.Moreover, the critical mode of a conventional control system, as shown in Figure 2, is determined based on the linearized model represented in Reference [25].Therefore, it is not repeated, and the interested reader can refer to it for more details.Figure 7 investigates the impact of grid SCR uncertainties (1 < SCR < 10), different PLL bandwidths (8 Hz < f PLL < 80 Hz), and inverter power levels (0.01 P n < P < P n ) on the critical mode under both proposed and conventional control methods.The following conclusions can be drawn from this figure: Step 3: Increasing q2 adjusts the minimum damping ratio.
Step 4: Increasing r reduces control gain matrix norm.
Step 5: Increasing q3 changes dc-link voltage and inverter current transient response.
As it can be seen, the design procedure is very straightforward, and repetitive actions are reduced considerably.

Robustness Analysis
The previous design process can provide sufficient stability margins.However, these excellent stability margins are guaranteed for the nominal values of system parameters.However, some of these parameters may not be constant during the inverter operation.For example, their operating-point conditions continuously change during the day, and the grid SCR can vary over a wide range due to variations of the connected load and generators.
In this respect, a robustness assessment based on the eigenvalue calculation was carried out to examine the system's robustness under these uncertainties.In such a way, the variations in the critical mode are studied for different values of grid SCR and operating point conditions.In addition, a comparison between the proposed and the conventional control methods is also given.
The stability margins of the proposed control method are calculated based on the linearized model in (6) and feedback control gain in (10) ( ) . It is worth remarking that the control gain matrix in (10) is calculated for the nominal power (Pn), the average value of grid SCR (SCR = 5.5), and fPLL = 10 Hz.Then the control gain matrix is kept constant during the robustness assessment.Moreover, the critical mode of a conventional control system, as shown in Figure 2, is determined based on the linearized model represented in Reference [25].Therefore, it is not repeated, and the interested reader can refer to it for more details.Figure 7 investigates the impact of grid SCR uncertainties (1 < SCR < 10), different PLL bandwidths (8 Hz < fPLL < 80 Hz), and inverter power levels (0.01 Pn < P < Pn) on the critical mode under both proposed and conventional control methods.The following conclusions can be drawn from this figure: (1) In Figure 7a, the grid SCR varies from 10 (strong grid condition) to 1 (ultra-weak grid condition); thus, the assessment results can cover all possible situations in practice.As it can be seen, under a lower grid SCR, the critical mode moves to a hazardous area under the conventional control method.Moreover, it covers a larger area, meaning higher sensitivity to grid SCR variation.(1) In Figure 7a, the grid SCR varies from 10 (strong grid condition) to 1 (ultra-weak grid condition); thus, the assessment results can cover all possible situations in practice.As it can be seen, under a lower grid SCR, the critical mode moves to a hazardous area under the conventional control method.Moreover, it covers a larger area, meaning higher sensitivity to grid SCR variation.
For example, as shown in Figure 7a, the system equipped with the conventional control system becomes unstable under a SCR lower than 1.37, even for a low PLL bandwidth (f PLL = 10 Hz).In contrast, the proposed control method employs a higher PLL bandwidth (f PLL = 50 Hz) and provides excellent stability margins under a wide range of SCR variations, from 10 to 1.The maximum damping factor is lower than −30 rad/s under different grid conditions.
It is worth mentioning that employing a high PLL bandwidth provided by the proposed control method is essential for good FRT capability and disturbance rejection.
(2) Figure 7b investigates the impact of different PLL bandwidths under the SCR = 1.5.The system under the conventional control method becomes unstable for PLL bandwidths higher than 17 Hz, while the proposed control method remains stable for all considered ones.Moreover, the critical mode covers a large area under the conventional ones, meaning higher sensitivity to the PLL bandwidths and severe control loops interaction.In contrast, the impact of the synchronization block on the system's critical mode is considerably reduced under the proposed control method.
It is worth noting that, unlike the grid SCR, the PLL bandwidth is a known and, usually, fixed parameter; thus, the design procedure can be performed for the desired value of PLL bandwidth to guarantee the required stability margin and the optimal control gain calculation, as shown in Figure 7a,c.
However, additional functionalities are recently suggested for VSCs to support power grids, e.g., transient stability enhancement, virtual inertia implementation, and fault ridethrough capability.To obtain these objectives, the VSC may be equipped with a PLL designed with adaptive or gain scheduling techniques.In addition, it would be interesting to know how the PLL bandwidth can affect the system response under the proposed control method and the possibility of implementing a fast PLL.They can be considered reasons to provide Figure 7b.
(3) The operational point variations on the system stability should also be investigated since the linearized model around the specific operating point is used to calculate the control gains.As shown in Figure 7c, the inverter power level under a weak grid condition (SCR = 1.5) affects the damping factor and frequency of oscillations; however, the system remains stable under different power levels and both control methods.Under the conventional control method, increasing the inverter power level reduces the stability margin due to higher coupling and interaction between different control loops, VSC, and the power grid.Moreover, the frequency of oscillation is higher at lower powers and vice versa.It can be supposed that the control system's high-frequency parts contribute to the critical mode at lower loading.Moreover, low-frequency control parts affect the critical mode at higher loading.Finally, note that the stability margin is linearly improved under the proposed control method by increasing the inverter power level from lower values to nominal.The reason is that the optimal control gain is calculated for the nominal power, and calculated control gain can provide the optimal performance under this condition.Therefore, we can expect stability to deteriorate as the power level moves away from the nominal value.

Control Gain Matrix Simplification
The state-feedback-based active damping can considerably improve system performance and robustness, as presented in previous subsections.However, grid-side current sensors are required, and computations are increased due to employing all states.
Therefore, it would be worth investigating how the number of sensors and calculations can be reduced.A rational solution is employing a full-or reduced-order Luenberger observer and Kalman filter to eliminate the grid-side sensors.However, it increases the computations and adds additional observer gains that should be tuned, leading to the higher complexity of the control system design and implementation.In this subsection, a more straightforward solution is proposed.First, it is assumed that all states are available, and the optimal control gain matrix is calculated.
Then it is assumed that the mth state is not available, so the related feedback gains are set to zero x m : off → K f (:, m) = 0 .Finally, the impact of eliminating this state on the system stability and performance is investigated by eigenvalues' analysis.If the critical mode and stability requirements have remained unchanged, this state can be simply excluded.Otherwise, this state and the related gains must be used.For example, if an eliminating grid-side current i gd , i gq is required, the thirteenth and fourteenth columns of K f must be set to zeros x 13 &x 14 : off → K f (:, 13 : 14) = 0 .Moreover, based on the new control gain matrix, the analysis of the eigenvalues can be carried out again.
Figure 8 shows the critical mode location when one of the states and its related gains are kept at zeros.The number next to each mode expresses the state number that is omitted.Moreover, 00 is equal to the case that all states are used.For example, as shown, eliminating states 8, 9, or 15 causes system instability, whereas when states 11, 13, or 14 are eliminated, the system remains stable.To examine the analytical results, time-domain simulations are presented in Figure 9.
Appl.Sci.2022, 12, x FOR PEER REVIEW 14 of 21 excluded.Otherwise, this state and the related gains must be used.For example, if an eliminating grid-side current (  ′ ,   ′ ) is required, the thirteenth and fourteenth columns of Kf must be set to zeros ( 13 & 14 :off →   (: ,13: 14) = 0).Moreover, based on the new control gain matrix, the analysis of the eigenvalues can be carried out again.
Figure 8 shows the critical mode location when one of the states and its related gains are kept at zeros.The number next to each mode expresses the state number that is omitted.Moreover, 00 is equal to the case that all states are used.For example, as shown, eliminating states 8, 9, or 15 causes system instability, whereas when states 11, 13, or 14 are eliminated, the system remains stable.To examine the analytical results, time-domain simulations are presented in Figure 9.
This figure illustrates the system outputs under three different control gain matrixes, where all states are available, x13 is not available ( 13 :off →   (: ,13) = 0), and x9 is not available ( 9 :off →   (: ,9) = 0).As shown, eliminating x13 does not lead to system instability.On the other hand, eliminating x9 does.
As shown in Figure 8, there are seven states that eliminating them does not lead to system instability.It would be worth investigating the changes of critical mode locations where all of them are off simultaneously.This action reduces the number of employed states from 15 to 8.Moreover, 46% of the control gain matrix elements become zeros, reducing computations:  3 ,  4 ,  7 ,  11 ,  12 ,  13 & 14 :off →   (: ,3: 4) =   (: ,7) =   (: ,11: 14) = 0.
The impact of grid SCR variation on the critical mode under the conventional control method and the proposed full and reduced state-feedback control methods are examined in Figure 10.The system robustness is a little worse than full state feedback by reduced state feedback; as expected, however, it is still considerably better than the conventional one.In summary, the system remains stable under the reduced state feedback control, with a good stability margin over a wide range of grid SCR variations.At the same time, the possibility of implementing a high-performance PLL is achieved.excluded.Otherwise, this state and the related gains must be used.For example, if an eliminating grid-side current (  ′ ,   ′ ) is required, the thirteenth and fourteenth columns of Kf must be set to zeros ( 13 & 14 :off →   (: ,13: 14) = 0).Moreover, based on the new control gain matrix, the analysis of the eigenvalues can be carried out again.
Figure 8 shows the critical mode location when one of the states and its related gains are kept at zeros.The number next to each mode expresses the state number that is omitted.Moreover, 00 is equal to the case that all states are used.For example, as shown, eliminating states 8, 9, or 15 causes system instability, whereas when states 11, 13, or 14 are eliminated, the system remains stable.To examine the analytical results, time-domain simulations are presented in Figure 9.
This figure illustrates the system outputs under three different control gain matrixes, where all states are available, x13 is not available ( 13 :off →   (: ,13) = 0), and x9 is not available ( 9 :off →   (: ,9) = 0).As shown, eliminating x13 does not lead to system instability.On the other hand, eliminating x9 does.
As shown in Figure 8, there are seven states that eliminating them does not lead to system instability.It would be worth investigating the changes of critical mode locations where all of them are off simultaneously.This action reduces the number of employed states from 15 to 8.Moreover, 46% of the control gain matrix elements become zeros, reducing computations:  3 ,  4 ,  7 ,  11 ,  12 ,  13 & 14 :off →   (: ,3: 4) =   (: ,7) =   (: ,11: 14) = 0.
The impact of grid SCR variation on the critical mode under the conventional control method and the proposed full and reduced state-feedback control methods are examined in Figure 10.The system robustness is a little worse than full state feedback by reduced state feedback; as expected, however, it is still considerably better than the conventional one.In summary, the system remains stable under the reduced state feedback control, with a good stability margin over a wide range of grid SCR variations.At the same time, the possibility of implementing a high-performance PLL is achieved.This figure illustrates the system outputs under three different control gain matrixes, where all states are available, x 13 is not available x 13 : off → K f (:, 13) = 0 , and x 9 is not available x 9 : off → K f (:, 9) = 0 .As shown, eliminating x 13 does not lead to system instability.On the other hand, eliminating x 9 does.As shown in Figure 8, there are seven states that eliminating them does not lead to system instability.It would be worth investigating the changes of critical mode locations where all of them are off simultaneously.This action reduces the number of employed states from 15 to 8.
Moreover, 46% of the control gain matrix elements become zeros, reducing computations: x 3 , x 4 , x 7 , x 11 , x 12 , x 13 &x 14 : off → K f (:, The impact of grid SCR variation on the critical mode under the conventional control method and the proposed full and reduced state-feedback control methods are examined in Figure 10.The system robustness is a little worse than full state feedback by reduced state feedback; as expected, however, it is still considerably better than the conventional one.In summary, the system remains stable under the reduced state feedback control, with a good stability margin over a wide range of grid SCR variations.At the same time, the possibility of implementing a high-performance PLL is achieved.

Simulation and Experimental Results
In this section, we recount how simulation and experimental tests were carried out to validate the performance of the proposed control method and compare the results with the conventional one.The simulation model was prepared in MATLAB/Simulink with parameters as close as possible to the real power system.The control and system parameters are given in Table 3.Based on that, various analyses and tests were carried out to reveal system performance under different conditions, which might not be possible by experiments due to some practical limitations (e.g., large inductors values, power, and current limitations).Finally, analytical and simulation results were validated by conducting experiments on a laboratory setup.It is worth noting that a downscaled power system is designed to emulate weak grid conditions and reduce the required inductors in practice.Table 4 shows the parameters that differ from the main circuit parameters.The common parameters by the main circuit are not repeated in this table.
The experimental setup includes two 5 kW Danfoss voltage source converters operating back-to-back, LCL filters, a grid simulator (Chroma 61845), and a DS1007 dSPACE system, as shown in Figure 11.
The first VSC acts as a constant power source, emulating a primary source such as

Simulation and Experimental Results
In this section, we recount how simulation and experimental tests were carried out to validate the performance of the proposed control method and compare the results with the conventional one.The simulation model was prepared in MATLAB/Simulink with parameters as close as possible to the real power system.The control and system parameters are given in Table 3.The series resistance of C f (r c ) 0.5 (mΩ) Proposed control method The series resistance of L f (r f ) 1 (mΩ) 1.5, 1 × 10 4 r, q 1 The series resistance of L g (r g ) 1 (mΩ) 1, 5 q 2 , q 3 DC-link capacitor (C dc ) 1. Based on that, various analyses and tests were carried out to reveal system performance under different conditions, which might not be possible by experiments due to some practical limitations (e.g., large inductors values, power, and current limitations).Finally, analytical and simulation results were validated by conducting experiments on a laboratory setup.It is worth noting that a downscaled power system is designed to emulate weak grid conditions and reduce the required inductors in practice.Table 4 shows the parameters that differ from the main circuit parameters.The common parameters by the main circuit are not repeated in this table.
The experimental setup includes two 5 kW Danfoss voltage source converters operating back-to-back, LCL filters, a grid simulator (Chroma 61845), and a DS1007 dSPACE system, as shown in Figure 11.The second VSC is connected to a grid simulator through an LCL filter and transfers the received power to the power grid.It contains CC, VFF, PLL, AVC, and DVC.A DS1007 dSPACE system platform is used to realize both control methods and drive VSCs. Figure 12 shows the steady-state performance of the proposed control method under strong and ultra-weak grid conditions (SCR = 10 and 1).This figure confirms a robust and stable operation over a wide range of grid SCR variations, especially critical scenario (SCR = 1), while a high-performance PLL is implemented.
The system responses of the conventional and the proposed control methods under different grid SCRs are shown in Figure 13.As it can be seen, the performance of the conventional control method is degraded under a weaker grid, and it becomes unstable under the weakest one (SCR = 1.25).
In comparison, the proposed control method can provide a stable response under weak grids.The simulation results also confirm the robustness analysis in the previous section.
The system performance under different PLL bandwidths is investigated in Figure 14.As shown, under the conventional control method, it is impossible to implement fast PLL in weak grids.In contrast, the proposed control method can facilitate the implementation of high-performance and fast PLL.Moreover, the system is stable over a wide range of PLL bandwidth selections.Therefore, the PLL bandwidth can be selected based on the other performance requirements, such as fault ride-through (FRT) capability, large-signal stability, and ancillary services such as virtual inertia, and not only restricted limitations imposed by small-signal stability.The first VSC acts as a constant power source, emulating a primary source such as PV or wind turbine and supplies the second three-phase grid-connected VSC.
The second VSC is connected to a grid simulator through an LCL filter and transfers the received power to the power grid.It contains CC, VFF, PLL, AVC, and DVC.A DS1007 dSPACE system platform is used to realize both control methods and drive VSCs.
Figure 12 shows the steady-state performance of the proposed control method under strong and ultra-weak grid conditions (SCR = 10 and 1).This figure confirms a robust and stable operation over a wide range of grid SCR variations, especially critical scenario (SCR = 1), while a high-performance PLL is implemented.
The system responses of the conventional and the proposed control methods under different grid SCRs are shown in Figure 13.As it can be seen, the performance of the conventional control method is degraded under a weaker grid, and it becomes unstable under the weakest one (SCR = 1.25).
In comparison, the proposed control method can provide a stable response under weak grids.The simulation results also confirm the robustness analysis in the previous section.The dynamic response of the VSC under 20% voltage sag at the grid voltage is depicted in Figure 15.The dynamic response of the VSC under 20% voltage sag at the grid voltage is depicted in Figure 15.The system performance under different PLL bandwidths is investigated in Figure 14.As shown, under the conventional control method, it is impossible to implement fast PLL in weak grids.In contrast, the proposed control method can facilitate the implementation of high-performance and fast PLL.Moreover, the system is stable over a wide range of PLL bandwidth selections.Therefore, the PLL bandwidth can be selected based on the other performance requirements, such as fault ride-through (FRT) capability, large-signal stability, and ancillary services such as virtual inertia, and not only restricted limitations imposed by small-signal stability.The dynamic response of the VSC under 20% voltage sag at the grid voltage is depicted in Figure 15.The dynamic response of the VSC under 20% voltage sag at the grid voltage is depicted in Figure 15.The conventional control method has a larger rise time.However, it can be expected when outer voltage control loops are considered, and the system operates under a very weak grid connection, as discussed in References [23,32], where a rise time equal to or smaller than 100 ms was suggested for AVC as a compromise between the system stability and performance.In contrast, the proposed control method shows an excellent transient performance and disturbance-rejection capability that are even higher than the required levels by standards.
Furthermore, the transient performance of both control methods under step changes in the dc voltage reference is shown in Figure 16.Both control methods can provide a good transient response; however, the proposed one can respond a little faster with lower oscillations.
(a) The conventional control method has a larger rise time.However, it can be expected when outer voltage control loops are considered, and the system operates under a very weak grid connection, as discussed in References [23,32], where a rise time equal to or smaller than 100 ms was suggested for AVC as a compromise between the system stability and performance.In contrast, the proposed control method shows an excellent transient performance and disturbance-rejection capability that are even higher than the required levels by standards.
Furthermore, the transient performance of both control methods under step changes in the dc voltage reference is shown in Figure 16.Both control methods can provide a good transient response; however, the proposed one can respond a little faster with lower oscillations.
Finally, the feasibility of the proposed control method is examined in practice.The performance of both control methods under weak grid conditions (SCR = 1.21) is shown in Figures 17 and 18.Experiments show their capability under different PLL bandwidths and 20% voltage sag.As shown, the proposed control method can employ a fast and high-performance PLL under weak grids, while the conventional one cannot.Moreover, the conventional one becomes unstable under a 20% voltage sag.In contrast, the proposed control method has an excellent disturbance rejection and fault ride-through capability.and performance.In contrast, the proposed control method shows an excellent transient performance and disturbance-rejection capability that are even higher than the required levels by standards.
Furthermore, the transient performance of both control methods under step changes in the dc voltage reference is shown in Figure 16.Both control methods can provide a good transient response; however, the proposed one can respond a little faster with lower oscillations.

Conclusions
This paper proposes a new design approach for robust and optimal active damping control of voltage source converters (VSCs) for grid applications.Based on the proposed formulation, all control gains are considered in one control gain matrix, leading to a linear optimization problem.Then different state-of-the-art techniques can be employed to calculate control gains in such a way the control system meets the desired stability and performance criteria.Among others, this paper used the linear-quadratic regulator (LQR) to calculate all control gains optimally and straightforwardly.In this respect, almost direct relationships between the tuning parameters and stability and performance indices such

Conclusions
This paper proposes a new design approach for robust and optimal active damping control of voltage source converters (VSCs) for grid applications.Based on the proposed formulation, all control gains are considered in one control gain matrix, leading to a linear optimization problem.Then different state-of-the-art techniques can be employed to calculate control gains in such a way the control system meets the desired stability and performance criteria.Among others, this paper used the linear-quadratic regulator (LQR) to calculate all control gains optimally and straightforwardly.In this respect, almost direct relationships between the tuning parameters and stability and performance indices such as maximum damping factor, minimum damping ratio and control gain matrix are identified, which considerably reduces the iterative process in the conventional control system design.
Thanks to the proposed control structure and optimal control gains calculation, the proposed method provides a fast dynamic response and significantly increases the system robustness and disturbance rejection capability under ultra-weak grid conditions.For example, the possibility of implementing a PLL bandwidth of more than 50 Hz under a grid SCR lower than two is obtained, which effectively suppresses disturbances on the control system and conducts fault ride-through of grid-connected VSCs.
The performance of the proposed control method under different operating conditions is comprehensively investigated by simulations and experiments and compared with the conventional one.The obtained results demonstrate the performance of the proposed control method under various operating conditions.
For future works, an extension of the proposed methodology to other VSC-based power system applications, such as active power filters and uninterruptable power supplies, is recommended.For these applications, extra resonators at desired harmonics are required to compensate for nonlinear loads, making the control system design difficult.

Figure 1 .
Figure 1.Single line diagram of an LCL-filtered grid-connected three-phase voltage source con verter (VSC).

Figure 1 .
Figure 1.Single line diagram of an LCL-filtered grid-connected three-phase voltage source converter (VSC).

Figure 3 .
Figure 3.The proposed solution is to include state-feedback-based active damping in the inner outer control loops.

Figure 4 .
Figure 4.The proposed control method.

Figure 3 .
Figure 3.The proposed solution is to include state-feedback-based active damping in the inner and outer control loops.

Figure 3 .
Figure 3.The proposed solution is to include state-feedback-based active damping in the inne outer control loops.

Figure 4 .
Figure 4.The proposed control method.

Figure 4 .
Figure 4.The proposed control method.

Figure 8 .
Figure 8. Critical mode locations of a grid-connected VSC in the complex s-plane (a) SCR = 1.5, and (b) SCR = 5.

Figure 9 .
Figure 9. Time-domain simulations showing the system response under different control gain matrixes.

Figure 8 .
Figure 8. Critical mode locations of a grid-connected VSC in the complex s-plane (a) SCR = 1.5, and (b) SCR = 5.

Figure 8 .
Figure 8. Critical mode locations of a grid-connected VSC in the complex s-plane (a) SCR = 1.5, and (b) SCR = 5.

Figure 9 .
Figure 9. Time-domain simulations showing the system response under different control gain matrixes.

Figure 9 .
Figure 9. Time-domain simulations showing the system response under different control gain matrixes.

Figure 10 .
Figure 10.Impact of grid SCR variations on critical mode locations of a grid-connected VSC in the complex s-plane (P = Pn and fPLL = 10 Hz (conventional controller) and 50 Hz (proposed controllers)).

Figure 10 .
Figure10.Impact of grid SCR variations on critical mode locations of a grid-connected VSC in the complex s-plane (P = P n and f PLL = 10 Hz (conventional controller) and 50 Hz (proposed controllers)).

Figure 13 .
Figure 13.Time-domain simulations showing the system response under different grid SCR (SCR = 1.5, 1.35, and 1.25; P = 10 kW): (a) the conventional control method and (b) the proposed control method.

Figure 13 .Figure 14 .
Figure 13.Time-domain simulations showing the system response under different grid SCR (SCR = 1.5, 1.35, and 1.25; P = 10 kW): (a) the conventional control method and (b) the proposed control method.

Figure 14 .
Figure 14.Time-domain simulations showing the system response under different PLL bandwidths (P = 10 kW, SCR = 1.5):(a) the conventional control method and (b) the proposed control method.

Figure 15 .
Figure 15.Time-domain simulations showing the system response under 20% voltage sag at the grid voltage (P = 10 kW, SCR = 2): (a) the conventional control method and (b) the proposed control method.

Figure 15 .
Figure 15.Time-domain simulations showing the system response under 20% voltage sag at the grid voltage (P = 10 kW, SCR = 2): (a) the conventional control method and (b) the proposed control method.

Figure 16 .Figure 16 .
Figure 16.Time-domain simulations showing the system response under several step changes in the dc voltage reference (P = 10 kW, SCR = 2): (a) the conventional control method and (b) the proposed control method.Finally, the feasibility of the proposed control method is examined in practice.The performance of both control methods under weak grid conditions (SCR = 1.21) is shown in Figures17 and 18.Experiments show their capability under different PLL bandwidths and 20% voltage sag.As shown, the proposed control method can employ a fast and highperformance PLL under weak grids, while the conventional one cannot.Moreover, the

Figure 17 .
Figure 17.Experimental results showing the system response under different PLL bandwidths (SCR = 1.21,Lg = 15.5 mH and P = 5 kW): (a,b) the conventional control method and (c,d) the proposed control method.
and Current (A)V g = 1 pu V g = 0.8 pu

Figure 17 .
Figure 17.Experimental results showing the system response under different PLL bandwidths (SCR = 1.21,L g = 15.5 mH and P = 5 kW): (a,b) the conventional control method and (c,d) the proposed control method.

Figure 17 .Figure 18 .
Figure 17.Experimental results showing the system response under different PLL bandwidths (SCR = 1.21,Lg = 15.5 mH and P = 5 kW): (a,b) the conventional control method and (c,d) the proposed control method.

Figure 18 .
Figure 18.Experimental results showing the system response under 20% voltage sag at the grid voltage (SCR = 1.21,L g = 15.5 mH, and P = 5 kW): (a) the conventional control method (f PLL = 10 Hz) and (b) the proposed control method (f PLL = 50 Hz).

Table 1 .
State-of-the-art control-system design techniques.

Table 3 .
Main system and control parameters.

Table 3 .
Main system and control parameters.

Table 4 .
Downscaled system and control parameters.

Table 4 .
Downscaled system and control parameters.