Research Article  Open Access
Alireza Khadem Abbasi, Mohd Wazir Mustafa, "Mathematical Model and Stability Analysis of InverterBased Distributed Generator", Mathematical Problems in Engineering, vol. 2013, Article ID 195038, 7 pages, 2013. https://doi.org/10.1155/2013/195038
Mathematical Model and Stability Analysis of InverterBased Distributed Generator
Abstract
This paper presents a mathematical (smallsignal) model of an electronically interfaced distributed generator (DG) by considering the effect of voltage and frequency variations of the prime source. Dynamic equations are found by linearization about an operating point. In this study, the dynamic of DC part of the interface is included in the model. The stability analysis shows with proper selection of system parameters; the system is stable during steadystate and dynamic situations, and oscillatory modes are well damped. The proposed model is useful to study stability analysis of a standalone DG or a Microgrid.
1. Introduction
Distributed generation (DG) systems have been expected to be an important electric power supply system for next generation. DGs are able to be installed near the loads, so they can increase the power quality and reliability of electricity delivered to sensitive loads. Some of the DG technologies require a power electronics interface in order to convert the energy into the grid compatible AC power. These interface devices make the sources more flexible in their operation and control compared to the conventional electrical machines. However, due to their negligible physical inertia, they also make the system potentially susceptible to oscillation resulting from network disturbances [1].
The coordinated operation and control of DGs together with loads and storage devices are central to the concept of microgrid [2]. The analysis of the dynamic stability of conventional power systems is well established, but for microgrid there is a need to investigate how circuit and control features give rise to particular oscillatory modes, and which of these have poor damping. Finding an exact mathematical model by considering DGs and their control is needed to investigate dynamic stability of the microgrid under transient events such as islanding from main grid and smallsignal deviation like slow changing in load.
Reference [3] presented a smallsignal model for inverter in standalone AC supply system by using only the droop controller's variables as state variables of the DG model. An averaged current source model has been suggested for the converter in [4]. In this investigation, the high frequency converter current dynamics have been neglected in order to focus exclusively on the dynamics and control of the islanded microgrid. Regardless of the type of the DG, an equivalent RLC circuit including output filter and transformers impedances has been modeled as a power circuit of DG in most of the literatures [1, 5–8]. In these cases, output currents and voltages are considered as state variables. References [1, 6, 7] by adding controller equations to the DG system made an accurate model for small signal stability analysis of a microgrid.
In most articles, DC part of the DG (DC filter’s voltage and current) is not considered in dynamic studies, or it is assumed to be a constant value [6], so its effect is neglected in stability analysis. In the reference [8], the input DC voltage variation of the inverter has been represented as the external perturbation in the openloop model of the DG, but it has been set as a constant during control loop design and frequencydomain analysis. The DC voltage has been related to the input and output equations of electronically interface of the DG by using switching function of rectifier and inverter in [5]. Moreover, two types of models have been presented in this paper for the prime source of the DG. While the proposed models have been used in steadystate and load flow analysis, but their linear forms can be used also in dynamic studies.
The objective of this paper is to find a comprehensive dynamic model of DG, including the prime source, power electronically interfaces, output filter, and controller. The proposed model represents all components of the DG in a dq0 reference frame, thus it ensures any application such as steadystate and dynamic analysis that meets requirements and constraints of both AC and DC parts of the system.
The dynamic stability of the DG is investigated by the small signal and step response analysis. In this paper, the DG is connected only to the static load, but it can be extended into the network or typical microgrid for more applications.
2. System Studied
2.1. Power Circuit
A simplified singleline diagram of the studied DG is illustrated in Figure 1. DG is included of a threephase source, a threephase rectifier, an inverter, and a threephase static load on the output bus. The DG source can be representatives of wind turbine or microturbine generator for medium and low voltage levels, respectively. Regardless of the power measurement and control loops, the system shown in Figure 1 can be studied as a power circuit of DG in microgrid stability analysis.
2.2. Control System
As shown in Figure 1, the control of DG included power, voltage, and current controllers. The power control has been achieved by applying droop control for both active and reactive powers. The droop control for an electronically interfaced DG has the role of the governor for a synchronously generator. Where, by an increase in load, the reference frequency is decreased. In the similar way, the reactive power is controlled by a droop characteristic in the voltage magnitude.
In this study, the sintriangle modulation strategy is used to prepare applied voltage in the converter. The next parts of the control system are the outer voltage and inner current control loops, which are designed to reject highfrequency disturbances and provide sufficient damping for the output filter by using the proportional integral (PI) compensators.
3. System Equations
3.1. Power Circuit
The dynamic model of AC side of DG in threephase abc frame is obtained from where , , and are vectors of the instantaneous amounts of the output currents, bus voltages, and the converter output voltages, respectively. Then (1) is transformed to the  reference frame of DG that rotates at frequency related to the angular velocity of the voltage space vector of the bus. The transferred equations are where is the operator .
3.2. DC Interface Equations
The average voltage and current of capacitor and inductance of dc filter will change with time during transient. For the rectifier under operation without phase delay and with commutating inductance, we have [9] We also have From (3) and (4), Equation (5) expresses the relationship between dynamic parts of the dc filter and DG source voltage. The is the frequency of the source, and is the network side voltage’s component. The voltage is wrote in the form of Apart from the harmonics in the voltage waveform and considering sintriangle modulation for inverter, voltages become where is duty cycle.
We will see that the component of converter voltage is controlled by controller and is considered to be zero. The instantaneous power balance between two sides of the inverter is established by According to (7), we have Rewriting (6) with regard to (9), The dynamic equations of DC part are represented by (5) and (10).
3.3. Control System Equations
Figure 2 shows inverter power controller, where and are measured values of inverter active and reactive output powers. Moreover, and are filtered values of power after passing through the lowpass filter. As shown in Figure 3, the inverter output voltage and inverter output frequency are adjusted by the droop controller characteristics described by (11) and (12), respectively, as follows: Here, and represent the nominal frequency and axis output voltage set points, respectively. is the duty cycle ().
The power measurement outputs are given by (13) and (14), where is the cutoff frequency of the lowpass filter and is the Laplace operator
3.4. Common Reference Frame Theory
All equations of DG have been written in its reference frame. In power system analysis, it is necessary to relate all components’ reference frames to a common frame. This common reference frame would be represented as constant voltages in the synchronous reference frame. In microgrid studies, DG with the biggest nominal power is considered as a common reference frame.
Figure 3 shows the common reference frame of a microgrid by and axis rotating at the synchronous angular velocity of . The axis is in orientation of the voltage vector of the reference DG, . By assuming on the axis of the DG2 reference frame, , is the angle between the DG2 reference frame and common reference frame. In order to transfer variables in the DG2 reference frame to common reference frame, the transformation matrix is defined as To connect an inverter to the whole system the output variables that have the connection with the network need to be converted to the common reference frame. In this case the output variables are the output currents and the network side voltages . Because the dynamic of one DG is studied in this paper, the is zero, and transformation matrix is the unit.
3.5. SmallSignal Model
By liberalization of nonlinear differential equations around the operating steadystate point, the smallsignal model of the system has been created [9].
For inverter equations, first output variables must be transferred to the common reference frame and then be linearized. By using (15) inverter output currents and voltages in the common reference frame can be achieved as Then with linearization technique, we have where The subscript 0 denotes steadystate quantities. The could be rewritten as a function of output voltages in common reference frame as follows With linearization of (21), where For more information, see [3, 7]. Furthermore, the linear equation of the inverter output voltage could be written as By replacing (7) and (18) into the inverter equations and linearization with considering (19), (22), and (24) the linear equations of the DG have been arranged in the basic style of the linear differential equations as where and are control and state variables respectively. The state matrix of the system is expressed by that is used in dynamic analysis of the DG.
The output vector of the system is defined as a linear combination of control and state variables in the form of In the power system small signal stability analysis, the current injection into the network from the device is considered as the output signal [10]. For the DG, the output current is one of the state variables, so we have All matrices are available in Appendix A in detail.
Finally, the transfer function is defined by
3.6. Eigenvalue Analysis
In power systems smallsignal stability analysis is aimed to determine the properties of operation parameter variations that are independent from disturbance intensity. Eigenvalue analysis is used to show the information of different stability modes for power system smallsignal stability problems. A system is stable when all of its modes are stable. Furthermore, it is required that all oscillations are well and quickly damped. This technique is declared in [11] in details.
In this paper, the participation factors and eigenvalues of the state matrix are computed and analyzed. Eigenvalues with negative real parts show that the system is stable. the state variable with the highest normalized participation factor is the best choice for feedback signal if it may be a measurable physical variable.
4. Case Study and Results
To study smallsignal stability of the model, a 10 kVA DG connected to a 9 kW static load is considered as the case study. All DG information is taken from [1]. The output impedance of DG is . The cutoff frequency is considered 31.41 Rad/s.
The rectifier’s input inductance is used to amend the input current. The inductor with low value may not restrict the peak current, and its high amount may reduce the input voltage [10]. Hence, a 1.06 (mH) inductor is chosen as an optimized value.
The DC filter is involved in a small series inductance, , (here 0.03 mH) and a shunt capacitance, , to smooth the diode rectifier output. With regard to [10], a 165 μF/kW capacitance is selected in this study. For DG, we have Hence The steadystate values of the system are calculated by load flow program.
In order to realize the sensitivity and dynamic behaviour, eigenvalues trajectory plots are drawn as a function of some system parameters. Figure 4 shows the locations of eigenvalues as a function of variation of from 0.00000004 to 0.000094 with the rate of 0.000001. In this case, is equal to 0.0013.
Figure 5 shows the loci of eigenvalues corresponding to variation of from 0.00003 to 0.1 with the rate of 0.001. is considered to be 0.000015.
Trajectory plots show instability in modes of 1, 2, 3, and 4 if the and increase to big values.
With regard to and as constant values of 0.000015 and 0.0013, Figure 6 shows the root locus plot for the model as a function of the filter cutoff frequency from 37 to 25 with the rate of −0.5. It can be seen that eigenvalues 5 and 6 move to right side of the plot.
Results are used to adjust the gains of the power controller. With optimal values of controller gains and appropriate selection of the DG parameters, all oscillatory modes are well damped, and the proposed model is stable.
By using QRdecomposition method in the MATLAB software environment, the eigenvectors and eigenvalues of the state matrix have been computed in operating point. The DG system along with its controller is described by six eigenvalues: , , and .
The normalized participation factors of the dynamic model are represented in Table 1. It is a useful indication of the role of a state on the mode. As seen in oscillatory modes of 1, 2, 3, and 4, the state variables of DC part have the main role.

Figures 7 and 8 show step responses of the output current according to the source voltage and frequency. The small steadystate error in is caused by the effects of them on DC voltage and DC current. The variation in is due to coupling between and components of current. By control of rectifier and making a constant DC voltage this effects will be removed. Error specifications also could be adjusted by adding current and voltage control loops.
5. Conclusion
The mathematical model of the inverterbased DG has been proposed in this paper. The proposed model has included power and control circuits of the DG. The dynamics of the DC filter and the effects of the prime source have been also considered in the DG model. By using eigenvalue analysis, the stability of proposed model has been studied in detail. This study presented that all oscillatory modes are well damped with suitable selection of system parameters. Since the proposed model is stable, it can be used in the dynamic and steadystate study of a microgrid.
Appendices
Linear equations of the DG.
A. Power Circuit
A.1. AC Part
By substituting (18) into the (2) and linearization where By (23) and
A.2. DC Part
By rearranging (5) and (10) in the frame linear differential equation By replacing (19) and (22) into the (A.4),
A.3. Differential Equation of Power Circuit
B. Control Circuit
B.1. Active Power Droop Controller
Linear form of (11) and (13), By combining (B.1) and (B.2),
B.2. Reactive Power Droop Controller
In the same way of the active power control, from (12) and (14) we have
B.3. Differential Equation of Control Circuit
With writing the active and reactive power in terms of voltage and current
C. Linear Differential Equations of DG (25)
References
 N. Pogaku, M. Prodanović, and T. C. Green, “Modeling, analysis and testing of autonomous operation of an inverterbased microgrid,” IEEE Transactions on Power Electronics, vol. 22, no. 2, pp. 613–625, 2007. View at: Publisher Site  Google Scholar
 P. R. H. Lasseter and P. Paigi, “Microgrid: a conceptual solution,” in Proceedings of the IEEE 35th Annual Power Electronics Specialists Conference (PESC '04), pp. 4285–4290, June 2004. View at: Google Scholar
 E. Antônio, A. Coelho, P. C. Cortizo, P. Francisco, and D. Garcia, “Smallsignal stability for parallelconnected inverters in standalone AC supply systems,” IEEE Transactions on Industry Applications, vol. 38, no. 2, pp. 533–542, 2002. View at: Google Scholar
 C. K. Sao and P. W. Lehn, “Intentional islanded operation of converter fed microgrids,” in Proceedings of the IEEE Power Engineering Society General Meeting, Quebec, Canada, 2006. View at: Google Scholar
 H. Nikkhajoei and R. Iravani, “Steadystate model and power flow analysis of electronicallycoupled fistributed resource units,” IEEE Transactions on Power Delivery, vol. 22, no. 1, pp. 721–728, 2007. View at: Publisher Site  Google Scholar
 R. Majumder, A. Ghosh, G. Ledwich, and F. Zare, “Stability analysis and control of multiple converter based autonomous microgrid,” in Proceedings of the IEEE International Conference on Control and Automation (ICCA '09), pp. 1663–1668, Christchurch, New Zealand, December 2009. View at: Publisher Site  Google Scholar
 F. Katiraei, M. R. Iravani, and P. W. Lehn, “Smallsignal dynamic model of a microgrid including conventional and electronically interfaced distributed resources,” IET Generation, Transmission and Distribution, vol. 1, no. 3, pp. 369–378, 2007. View at: Publisher Site  Google Scholar
 Y. Zhang, Z. Jiang, and X. Yu, “Smallsignal modeling and analysis of parallelconnected voltage source inverters,” in Proceedings of the IEEE 6th International Power Electronics and Motion Control Conference (IPEMC '09), pp. 377–383, Wuhan, China, 2009. View at: Google Scholar
 P. C. Krause, O. Wasynczuk, and S. D. Sudhoff, Analysis of Electric Machinery and Drive Systems, John Wiley & Sons, New York, NY, USA, 2002.
 S. C. Vegunta, J. V. Milanović, and S. Z. Djokić, “Modelling of V Hz and vector controlled ASDs in PSCAD/EMTDC for voltage sag studies,” Electric Power Systems Research, vol. 80, no. 1, pp. 1–8, 2010. View at: Publisher Site  Google Scholar
 P. Kundur, Power System Stability and Control, McGrawHill, New York, NY, USA, 1994.
Copyright
Copyright © 2013 Alireza Khadem Abbasi and Mohd Wazir Mustafa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.