Modeling Single-Phase Inverter and Its Decentralized Coordinated Control by Using Feedback Linearization

It is a very crucial problem to make a microgrid operated reasonably and stably. Considering the nonlinear mathematics model of inverter established in this paper, the input-output feedback linearization method is used to transform the nonlinear mathematics model of inverters to a linear tracking synchronization and consensus regulation control problem. Based on the linear mathematics model and multiagent consensus algorithm, a decentralized coordinated controller is proposed to make amplitudes and angles of voltages from inverters be consensus and active and reactive power shared in the desired ratio. The proposed control is totally distributed because each inverter only requires local and one neighbor’s information with sparse communication structure based on multiagent system. The hybrid consensus algorithm is used to keep the amplitude of the output voltages following the leader and the angles of output voltage as consensus. Then the microgrid can be operated more efficiently and the circulating current between DGs can be effectively suppressed. The effectiveness of the proposed method is proved through simulation results of a typical microgrid system.


Introduction
More and more AC and DC distributed generations (DGs) are connected to the systems which requires more flexible and modular control systems.Hence, microgrids can play this role when mass grid independency is unfitted.In a microgrid, DGs are controlled to deliver desired power to the system; meanwhile, to minimize circulating current is another important control objective in order to keep microgrid be operated stably.In general, many control methods for such a microgrid system are proposed which can be categorized as centralized control [1,2] and decentralized control [3][4][5][6].Both of the two control methods are based on the conventional droop control [7][8][9][10].
The centralized control requires a complex communication network with two-way communication that can adversely affect the stability of the system.The droop control achieves the active and reactive power being shared proportionally by controlling system frequency and fundamental voltage amplitudes conventionally [11].Since no communication is needed between inverters by this method, it is attractive.However, the conventional droop controller has several drawbacks that limit its application area, such as loaddependent frequency and deviation of voltage amplitude, high dependency on the inverter output-impedance, and large circulating current between inverters [9,[12][13][14].In [15], a strategy to improve the control technique is proposed by changing the droop coefficient and estimating the effect of the line impedance value.However, this strategy is quite complicated and sensitive [16].Meanwhile, the conventional droop control is a linear method, but the inverter is a nonlinear model [10]; thus, the microgrid cannot be effectively controlled only by droop control.Thus, the nonlinear factors should be considered in the control method.The input-output feedback linearization is used to transform the nonlinear heterogeneous dynamics of DGs to linear system [17].Meanwhile, in order to make power shared without deviation in either the frequency or the amplitude of the output voltage, some communication should also be added in the system.
Considering the above limitations, the decentralized coordinated control with minimal communication is an Mathematical Problems in Engineering effective way to control DGs.Meanwhile, the distance between the inverters may make an intercommunication between all inverters impractical in a real system.
Recently, multiagent systems are extensively studied due to their flexibility and computational efficiency.The authors of [18] design the consensus algorithm considering an active leader agent followed by other agents.Considering that every DG unit is an agent in the multiagent system, the decentralized coordinated control with multiagent consensus protocol is used in this paper.In addition, because slight differences in phases and amplitudes of the inverters output voltages can cause large circulating currents [19], circulating currents are not well suppressed only by controlling the output voltage amplitude.A new droop control based on the angle of the output voltage, referred to as angle droop, is proposed and has more advantages than frequency droop control [20,21].Thus, in order to regulate the angle and amplitudes of voltage at the same time to suppress the circulating currents, the angleactive droop control is used to combine with the reactivevoltage control in this paper.Then the amplitudes of output voltages are controlled by the leader consensus algorithm to track the leader's information which comes from one DG in the microgrid.Meanwhile, the angles of output voltages are controlled by the leaderless consensus algorithm to achieve the average consensus.
In this paper, hybrid multiagent consensus algorithm is used to address the tracking synchronization problem about amplitudes of output voltages and the consensus regulation problem about angles by using the input-output feedback linearization.This method is superior to the scalar control by which only amplitudes of the voltage are controlled.A decentralized coordinated controller based on angle-active power droop control and voltage-reactive power droop control is designed to make active and reactive power shared proportionally and suppress circulating currents between DGs.The multiagent consensus algorithm is used to calculate the deviation of angle and amplitudes between inverters by using the leaderless consensus algorithm and the leader consensus algorithm.Then the calculation results are used as the feedback values.The main contributions of this paper are given as follows.
(1) By using feedback linearization method combined with hybrid multiagent consensus algorithm, the power can be shared proportionally and the angles and amplitudes of output voltage can be kept consensus.(2) By using hybrid multiagent consensus algorithm, the amplitude of output voltages and angles can be controlled more reasonably and economically.(3) The control method requires only a sparse communication structure which can be more reliable.(4) The circulating current between inverters can be effectively suppressed; thus, energy conservation and the stability operation of the microgrid can be achieved.
This paper is organized as follows.In Section 2, the dynamical model of single-inverter-based DGs is established.Section 3 discusses the process of the combination of the input-output feedback linearization and hybrid multiagent consensus algorithm.Section 4 discusses the analysis of circulating current between parallel inverters.In Section 5, the effectiveness of the proposed method is verified by comparing simulation results of conventional control method and proposed control method.Section 6 concludes the paper.

Nonlinear Model of Single-Phase Inverter
The proposed control method is designed based on the large signal dynamical model of single-phase inverters in DGs.There are power calculation, voltage and current controllers, and second order generalized integrator in the control loops to adjust the amplitude and angle of output voltage.The control block diagram of a DG based on a single-phase inverter is presented in Figure 1.
Because the single-phase inverters are studied in this paper, the output voltage and current from inverters should be divided into the - components first by using a second order generalized as shown in Figure 2 [22].
The phase of output voltage from the th DG   can be expressed as where   is the angle of voltage from th DG and   is the frequency of th DG.So, instead of using the common reference frame, regulating angles of output voltages can make the - reference frames rotated synchronously.
There are second order generalized integrator, park transformation, power calculation, low-pass filter, and droop controller in the power controller shown in Figure 3.The power controller can provide the amplitude reference of voltage V *  for the first stage bridge and the angle reference of output voltage  *  for the second stage bridge.The  park which is used to make the park transformation is calculated by making the reference V *  ≡ 0. The instantaneous active and reactive power can be expressed as where V  and V  are - axis output voltage components for th DG,   and   are - axis output current components for th DG, and p and q are instantaneous active and reactive power for th DG.The differential equations of the active and reactive power can be expressed as where   and   are the average active and reactive power for th DG and   is cutoff frequency of power filter.The angle-active power and amplitude-reactive power droop controllers are presented as where   is the angle-active droop parameter for th DG,   is the voltage-reactive droop parameter for th DG,  set is the set value of angle-active droop control for th DG, V set is the set value of voltage-reactive droop control for th DG,  *  is the reference of voltage angle for th DG, and V *  and V *  are the - axis reference of voltage amplitude in the voltage controller for th DG.Equations ( 3) are inserted into (2) and it yields The block diagram of the voltage controller is shown in Figure 4.The input differential equations of the voltage controller are presented as where   and   are the - axis auxiliary state variables in the voltage controller for th DG.
Then the assist terms in the voltage controller are as follows: where  is the feedforward parameter in the voltage controller,   is the filter capacitor for th DG, and ASS  and ASS  are the - axis auxiliary variables used to be the decoupling terms and the feedforward terms.
The outputs of the voltage controller are presented as where  *  and  *  are - axis current reference in the current controller for th DG,   and   are the - axis proportion parameters in the voltage controller for th DG, and   and   are the - axis integral parameters in the voltage controller for th DG.
The block diagram of the current controller is shown in Figure 5.The input differential equations of the current controller are presented as    where   and   are the - axis auxiliary state variables in the current controller for th DG and   and   are the - axis filter current components for th DG.The outputs of the current controller are presented as where V *  and V *  are - axis output of current controller for th DG,   is the filter inductance for th DG,   and   are the - axis proportion parameters in the current controller for th DG, and   and   are the - axis integral parameters in the current controller for th DG.
Assumption 1.The parameters in connection lines should meet the following conditions in (13) as shown in Figure 1.Both resistive and inductor in connected lines are considered [14].Consider where   and   are the line parameters in the microgrid.
The differential equations for the output LCL filter and output connector considering the - line parameters are as follows: From ( 1), ( 7)-( 8), (11), and ( 14), the model of th DG can be rewritten into a matrix formation as follows: where the state vector is The detailed expression of (), , and  can be found in the Appendix.

Designing Decentralized Coordinated Controller by Feedback Linearization
In this subsection, the hybrid multiagent consensus algorithm is used to control the amplitudes and angles of output voltages.To make amplitudes of output voltages follow the leader is considered the tracking synchronization problem and to make the angles of output voltages consensus is considered consensus regulation problem.Every DG can be seen as an agent with two control inputs and two outputs.The control signals are calculated by using the error between local and one of the neighbors information.The information about leader can come from one DG in the system.In this section, the graph theory and feedback linearization designation based on the multiagent protocol are illustrated, respectively.

Graph Theory.
The decentralized coordinated controller should use a communication network called directed graph which can be expressed as   = (  ,   ,   ) [23].The set of nodes in the network can be expressed as   = {V 1 , V 2 , . . ., V  }, the set of edges can be expressed as   ⊆   ×   , and the weighted adjacency matrix can be expressed as   = [  ] × with nonnegative adjacency element   .An edge rooted at node  and ended at node  is denoted by (V  , V  ), which means that information can be acquired from node  to node .For a graph with 0-1 adjacency elements, the in-degree and outdegree of node V  are defined as The degree matrix of diagraph   is a diagonal matrix Δ = [Δ  ], where Δ  = 0 for all  ̸ =  and Δ  = deg out (V  ).The graph Laplacian associated with the digraph   is defined as 3.2.Designation of Feedback Linearization.The  set and V set are selected as the control inputs to keep the angle   and the amplitude V  of the th DG consensus.Because V  is zero, the amplitude of output voltage is The inputs of the system are  = [ set V set ]  and the outputs of the system are  = [  V  ]  .From the above discussion, the model established is a MIMO nonlinear system.Then the feedback linearization method is used to design the decentralized coordinated controller.
The  set is used as the input to control angles of all output voltages to be consensus between DGs and the V set is used as the input to control the amplitude of output voltage to be consensus with the leader's voltage  leader,2 ≡ V ref .
Because it is difficult to find the relationship between input and output directly in such a nonlinear system, the input-output feedback linearization is used to transform the nonlinear system into linear system and find the relationship.
The process of feedback linearization in this MIMO nonlinear system described in ( 15) is as follows [24].
The control variables are shown up after one derivative of output; thus, the relative degree of this nonlinear system is [1,1].
Second, according to the relative degree, a decoupling matrix () which should be nonsingular should be calculated as given in (20) and another vector () is calculated at the same time as in (21).The   ℎ  () in ( 20) and ( 21) is the Lie derivative of ℎ  () with the respect to  and is defined as Consider ) , ( 20) Because the relative degree of this nonlinear system is [1, 1], the decoupling matrix is calculated as It is obvious that matrix   () is nonsingular.Additionally, matrix   () is given as where Third, the control law of a MIMO nonlinear system is defined as where 22) and ( 23), the control law   () can be defined as in (25).The appropriate V  () can be calculated as in (26) to provide the synchronization: where .
The nonlinear system can be transformed into linear system as where ,  ,1 =   , and  ,2 = V  , for all .
Using the feedback linearization, the nonlinear system ( 15) is transformed to the linear system (27).Equation ( 27) presents the th DG.It is obvious that the V ,1 and V ,2 can be designed, respectively.
From Section 3.1, DGs can communicate with their neighbors based on the digraph   and in this paper only one can communicate with the leader.The hybrid consensus algorithm consists of two kinds of consensus algorithms.One of the algorithms is called leaderless consensus algorithm which is used to control the angles of output voltages and the other one of the algorithm is called leader consensus algorithm which is used to control the amplitudes of output voltages.
The control of the angles from DGs is a consensus regulation problem which means the leaderless consensus algorithm is used.Thus, the errors about the angles between the DGs can be calculated as where  ,1 represents the error of angles of output voltage about the th DG.
To make all of the amplitude of output voltages equal to the leader's information, the tracking error about the amplitude of the output voltage and the reference information  leader,2 is given as where  ,2 represents the error of amplitudes of output voltage about the th DG.If node  observes the leader's amplitude information, edge (V leader,2 , V  ) is said to exist with weighting gain  ,2 = 1.Denote the pinning matrix by  2 = diag{ ,2 } ∈  × .Because a microgrid includes  DGs, the errors can be written into matrix as where (32 The system (28) can be rewritten as where The relationship between e 1 , e 2 and v 1 , v 2 should be designed in order to keep the system stable with the control input.The relationship can be established by the coupling gain and the feedback control gain.The design of the coupling gain and feedback control gain can be chosen according to the theory about optimal and nonlinear control in [24].For the sake of space, the designation of the feedback control gain and coupling gain is not included in this paper.
In this paper, only the stability of the system (33) is proved as follows.
Proof.Because the angle controller is the same as the amplitude controller, if one of the controllers is proved to be stable, the other one should be stable.In this section, the Ẏ 2 =   v 2 is proved to be stable under some conditions.
Inserting (31) into (33), the system can be as Defining  2 =  2 − leader,2 , choose the Lyapunov function as where   =  and  is the positive and definite matrix.Then the derivative of ( 35) is expressed as Then, Then, according to [25], the matrix −( +  2 ) can be Hurwitz.Therefore, given any positive number , the positive definite matrix  can be chosen as Then, inserting (38) into (37), Then, (39) shows that V ≤ 0. Thus, the internal dynamics are asymptotically stable.The proof is completed.

Circulating Current Analysis
In order to illustrate the effectiveness of proposed method focusing on the problem of circulating current between inverters, the relationship between a circulating current and amplitude and angle of output voltage is analyzed in this section.The simplified circuit schematic is shown in Figure 6.From the circuit scheme, the following circuit equations can be expressed for two inverters: where The load current is the phasor summation of the current of all the parallel inverters as given in (41).The load voltage is related to the load current and the impedance of load.From (41)-(44), the bus voltage of load V  can be presented as The circulating current between two inverters is equal to half of the difference current given in Inserting ( 45) into (46), the current can be presented by using the angle and amplitude of output voltage from inverters as Assuming that the equivalent load impedance of bus is resistive, the impedance   is replaced by   .Meanwhile, even though it is a microgrid in a small area, it is reasonable to illustrate that the impedance of the connection between inverters is inductive through designing the inductive virtual impedance [15]; in this paper, it is not discussed in detail.Thus, it means that   ≫   ; then the circulating current is given as follows: In order to analyze the active and reactive circulating current, respectively, (48) is decomposed into two parts as 2Δ İ =   −   . (48) For convenience, the following variables are defined: Then active and reactive current between inverters can be presented, respectively, as follows: In order to analyze the circulating current, it is assumed that  1 =  2 =   and then (51) and (52) can be rewritten as From ( 52) and (53), it is found that if the amplitude and angle of output voltage can be kept as consensus, the active and reactive current can be suppressed effectively.Thus, the proposed method can suppress the circulating current.

Simulation Results
Figure 7 shows the diagram of the typical microgrid test system.The configuration of DGs which are connected in series is considered in this microgrid.Because it is the typical microgrid system, the simulation results can be more convincing.The system data and simulation parameters are listed in Table 1.It is assumed that the total power demand in the microgrid can be supplied by the DGs and no load shedding is required.

Case 1: Conventional Controller.
In this case, the desired active and reactive ratio of DG-1, DG-2, and DG-3 is 1 : 1.33 : 2. The DG-1 is chosen as the system leader.The simulation process is as follows which is divided into 4 steps: first, at the beginning, load 1 , load 3 , load 4 , and load 5 are connected in the microgrid; second, at 0.3 s, load 2 is connected in the microgrid; third, at 0.6 s, load 1 and load 5 are disconnected from the system; fourth, at 0.8 s, load 5 is connected in the microgrid again.The load power is shared by DG-1, DG-2, and DG-3 all the time.The system can be operated stably.However, due to the weak controller, from 0.3 s to 1 s, the load cannot be shared in the desired ratio of 1 : 1.33 : 2 as shown in Figures 8 and 9.Meanwhile, the amplitudes and angles of output voltage from different DGs cannot be kept as consensus as shown in Figures 10 and 11.Because the power cannot be shared in the desired ratio shown in Figures 8 and 9, the amplitudes of output voltages are not in the desired ratio and the deviation angles between four outputs voltages are obvious, as shown in Figures 10 and 11.Thus, the circulating currents between DGs are very large in the system.In order to illustrate the problem of the conventional controller, Figure 11 is partially enlarged from 0.2 s to 1 s, which is overlaid on the low left of the original figures.

Case 2: Proposed Controller.
The performance of the proposed controller can be verified though the simulation in Case 2. In order to compare with the simulation results in Case 1, the microgrid is operated at similar situation as described in Case 1.The proposed controller is started at 0.2 s and the communication topology is shown in Figure 12.Compared with results in Case 1, the load can be shared in the desired ratio by DGs as shown in Figures 13 and 14.The amplitudes can be controlled to track with the leader's information and the angles of output from DGs can be kept as consensus as shown in Figures 15 and 16, respectively.Meanwhile, the deviation angles between four output voltages are very small as shown in Figure 15.From the discussion in Section 4, the circulating current between DGs can be Active power (W)  reduced largely.In addition, from Figures 15 and 16, it can also be concluded that the amplitude of output voltage can track the leader quickly.In addition, when the amplitudes of voltage are consensus, the power flow can be shared more reasonably.In order to illustrate the effectiveness of the proposed controller, Figure 16 is partially enlarged from

Figure 1 :Figure 2 :
Figure 1: Block diagram of a DG based on a single-phase inverter.

Figure 4 :
Figure 4: The block diagram of the voltage controller.

Figure 5 :
Figure 5: The block diagram of the current controller.

Figure 6 :
Figure 6: Simplified circuit diagram for two inverters.

Figure 12 :
Figure 12: The communication topology for proposed controller.