MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/6062798 6062798 Research Article Stability Analysis Method of Parallel Inverter Li Jun 1 http://orcid.org/0000-0001-9259-8270 Chen Jie 1 Xue Yaru 1 Qiu Ruichang 2 Liu Zhigang 1 Gajic Zoran 1 School of Electrical Engineering Beijing Jiaotong University Beijing China njtu.edu.cn 2 Beijing Rail Transit Electrical Engineering Technology Research Center Beijing China 2017 2112017 2017 12 05 2017 13 08 2017 2112017 2017 Copyright © 2017 Jun Li et al. 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.

In order to further provide theoretical support for the stability of an auxiliary inverter parallel system, a new model which covers most of control parameters needs to be established. However, the ability of the small-signal model established by the traditional method is extremely limited, so this paper proposes a new small-signal modeling method for the parallel system. The new small-signal model not only can analyze the influence of the droop parameters on the system performance, but also can analyze the influence of the output impedance of the inverter, the unbalanced and nonlinear loads, and the power calculation method and cut-off frequency of the low-pass filter on the system performance and stability. Based on this method, this paper carries out a comprehensive analysis on the performance of a parallel inverter system. And the correctness of the modeling method and analysis process of the system performance and stability are verified by the consistency of the simulation and experimental results.

Basic Scientific Research Project of Central University 2016RC038 National Key R&D Plan 2016YFB1200502
1. Introduction

The stability of a parallel system is one of the important indexes to evaluate the performance of a parallel system. For a high power auxiliary parallel system, a change of the various control parameters has an extremely large impact on the stability of the system. Therefore, it is very important and necessary to analyze the stability of the parallel system.

In general, the stability analysis of a parallel system should be based on a specific parallel control strategy. At present, the wireless parallel system usually adopts the droop method to achieve current sharing automatically , the basic approach of which is to use active power and reactive power to droop the amplitude and frequency of the voltage. The paper also adopts the droop method but has been modified to improve effectively  and the whole cycle adaptive droop method (TWCADM) is proposed. It is pointed out that a new power calculation method without a large delay low-pass filter, combined with the improved droop method, can greatly enhance the dynamic response of the parallel system. However, the author does not demonstrate the reason of the improvement of dynamic response in theory. Therefore, it is necessary to establish a unified model of the parallel system, using the mathematical tools to prove the improvement of the dynamic response in theory.

For the mathematical model of a parallel inverter system, the relevant literature [11, 12] is not abundant. Antonio et al. establish a small-signal model of grid-connected inverter for the first time , but he only considers the outer power loop and regards the inner voltage loop as the ideal state. This method is relatively simple, the established model is of 3rd-order, and, to a certain extent, it can reflect the change trend of a single inverter. However, a big drawback of this method is that the voltage of the parallel point is regarded as an independent variable, so that the parallel inverter automatically shields the influence of another inverter, and most papers adopt this method . Marwali et al. propose a small-signal modeling method which can consider the whole parallel system , and the voltage of a parallel point can be expressed as the physical variable of each inverter. All inverters in a parallel system can be coupled together by the voltage at a parallel point, and the inverters in the model can affect each other, so the established model is more accurate. But in Marwali et al.’s paper, the model established by the discrete modeling method in the Z domain is extremely complex. Actually, it is not greatly significant for the parallel system with large delay link to adopt discrete analysis; besides, this method fails to take the impact of the power calculation method into account.

Aiming at the defects of the methods above, this paper proposes a new method to establish the small-signal model of a parallel system, which can take the factors into account such as the coupling effect between the inverters in the parallel system, the influence of the line parameters, load characteristics, and the power calculation method and the cut-off frequency of the low-pass filter. The model established by this method can effectively analyze the influence of each correlation variable of the parallel inverter on the system performance, has a very wide universality, and is simpler than the model established by Marwali et al.

2. The Control Method of Parallel Inverters

Figure 1 is the topology of the parallel system, which indicates that the inverter adopts the three-phase four-wire system with split capacitor topology, and each inverter is controlled in two stages: parallel control  and inverter control , and the output voltage and output current of each phase need to be sampled.

The topology of auxiliary parallel inverter system.

Figure 2 is the schematic diagram of inverter control and we can obtain the output voltage of the inverter from the diagram:(1)uo=Gsuref-Zsio,where(2)Gs=GvsGis+1e-stpGvsGise-stp+sCGvse-stp+LCs2+rCs+1,Zs=Gise-stp+sL+rGvsGise-stp+sCGise-stp+LCs2+rCs+1.

The schematic diagram of inverter control.

In (1), G(s) is control branch and Z(s) is inverter internal impedance. Equation (1) indicates that both G(s) and Z(s) can be controlled by the voltage controller Gv(s) and the current controller Gi(s). If the appropriate voltage controller and current controller are selected, the internal impedance of the inverter can be designed.

In this paper, the improved resonant controller is used as the voltage controller and a proportional controller is used as the current controller.

Figure 3 is the bode diagram of the internal impedance of the inverter and we can discover an interesting result that the method shown in Figure 2 leads the internal impedance to be resistive. Therefore, the resistive droop method is suitable for the parallel system:(3)Eref=E0-mp+mdsP,φref=φ0-np+nis+ndsQ,where mp is the proportion a droop coefficient of active power, md is the differential coefficient of active power, np is the proportion of a droop coefficient of reactive power, ni is the integral coefficient of reactive power, and nd is the differential coefficient of reactive power. P and Q represent the active power and reactive power, respectively, and both of them are the values after the large delay low-pass filter.

The bode diagram of the internal impedance of the inverter.

Traditionally, the power calculation method is described by the following equation:(4)p~=U˙I˙o,q~=-U˙I˙o-90,where(5)U˙=Umaxcosωt+θu,I˙o=Imaxcosωt+θi.

In (5), Umax is the maximum value of voltage, θu is the initial voltage phase angle, Imax is the maximum value of current, and θi is the initial current phase angle.

In the power calculation method above, the obtained power contains an AC component at two-times the fundamental frequency. Therefore, a large delay low-pass filter is adopted to filter the AC component, which deteriorates the dynamic performance of the parallel system.

In order to improve the dynamic performance of the parallel system, the power calculation method can be improved as the following equation:(6)p~=0.5U˙I˙o+0.5U˙-90I˙o-90,q~=0.5U˙-90I˙o-0.5U˙I˙o-90.

Equation (6) indicates that both the voltage and current need to be delayed by 90 degrees which leads to more storage cost payment. However, the obtained power does not contain an AC component any more, which enables the cut-off frequency to be set higher, and the dynamic performance of parallel system can be effectively enhanced, as shown in(7)P=p~ωcs+ωc,Q=q~ωcs+ωc.

3. The Traditional Small-Signal Modeling Method

According to the above derivation, we can simplify the parallel inverter model, as shown in Figure 4, where Zθ is the output impedance of the inverter, consisting of two parts: internal impedance and line impedance.

The model of a parallel inverter.

From Figure 4, according to the concept of complex power, the output power of the inverter can be obtained:(8)P=1R2+X2UERcosφ+UEXsinφ-U2R,Q=1R2+X2-UERsinφ+UEXcosφ-U2X.

Considering the small disturbance existing near the steady-state value E,Ee,Ue, (8) can be linearised as(9)ΔP=ΔEα-ΔφEβ,(10)ΔQ=-ΔEβ-ΔφEα,where(11)α=1R2+X2URcosφ+UXsinφ,β=1R2+X2URsinφ-UXcosφ.

Based on the droop method mentioned in (3), the differential calculation result can be obtained as follows:(12)ΔE=-mp+mdsΔP,(13)Δφ=np+nis+ndsΔQ.

Combining (9), (10), (12), (13), and (7), Δφ and ΔE can be obtained:(14)Δφ=ωcs+ωcnp+ni1s+nds-ΔEβ-ΔφEα,ΔE=ωcs+ωc-mp-mdsΔEα-ΔφEβ.

Simplifying (14) the small-signal model can be obtained as follows:(15)s3Δφ+As2Δφ+BsΔφ+CΔφ=0,where(16)A=αmp+αωcmd+Eωcα2+β2ndmp+npmdαmd+Eωcα2+β2ndmd,B=αωcmp+Eωcα2+β2nimd+npmpαmd+Eωcα2+β2ndmd,C=Eωcα2+β2nimpαmd+Eωcα2+β2ndmd.

Equation (15) is the traditional small-signal model, and it can be found that the traditional model is of 3rd-order, so the traditional small-signal modeling method is relatively simple and apt to analyze. However, this method is not accurate, mainly for the following reasons:

( 1 ) The voltage at the parallel point is regarded as an independent variable, and the coupling effect between two inverters is eliminated. However, for an actual parallel system, there must be a coupling action between the two inverters, and the change of output characteristics of one inverter will inevitably lead to that of the other inverter. Therefore, the traditional small-signal modeling method is not accurate.

( 2 ) It is well known that the load characteristics have a great impact on the stability of the parallel system, but the traditional small-signal model cannot analyze the influence of the load characteristics.

( 3 ) The traditional small-signal model cannot analyze the influence of different power calculation methods.

4. The New Small-Signal Modeling Method

Aiming at the drawbacks of the traditional small-signal modeling method, this paper proposes a new small-signal modeling method. The model established by this modeling method has a wide enough amount of information and it can not only analyze most of the control parameters outside the voltage loop of the control system but also analyze the influence of the circuit parameters. What is more, the model is simpler than the model established by Marwali et al. .

For convenience, the parallel system shown in Figure 4 can be simplified as shown in Figure 5. The simplification is just to simplify the modeling and description process; readers can restore the output impedance of inverter and the load into any characteristics; it is even possible to consider the influence of two inverters on the parallel system with inconsistent impedance characteristics.

Simplified parallel model.

Based on Figure 5, the output current of the inverter can be obtained as(17)I˙o1=E˙1-U˙R,I˙o2=E˙2-U˙R.

The sum of the output current of two inverters is the load current and the output voltage can be obtained by the load current as(18)U˙=RoI˙o1+RoI˙o2.

Combining (17) and (18), we can obtain(19)U˙=RoE1cosφ1+E2cosφ2+jE1sinφ1+E2sinφ2R+2Ro.

From Figure 5, we can discover that the voltage phase angle at the parallel point is zero. In fact, the parallel voltage is sinusoidal and the phase angel is varied, and, for the convenience of modeling, the phase angel is set to the reference value so the phase angle φ1 and φ2 can be understood as the phase angle difference, rather than a time varying variable. The meaning of the equivalence is that the time varying variable can be transformed into an invariant, which is suitable for analysis at the small-signal model.

According to (19), we can get two equations as follows.(20)U=RoE1cosφ1+E2cosφ2R+2Ro,(21)0=E1sinφ1+E2sinφ2.

Equation (20) specifies the relationship between the parameters of the two inverters and the parallel point voltage. Equation (21) specifies the relationship between the steady-state voltages of the inverters.

Combining (19) and (17), we can obtain(22)I˙o1=E1cosφ1R+Ro-RoE2cosφ2RR+2Ro+jE1sinφ1R+Ro-RoE2sinφ2RR+2Ro,I˙o2=E2cosφ2R+Ro-RoE1cosφ1RR+2Ro+jE2sinφ2R+Ro-RoE1sinφ1RR+2Ro.

According to the traditional power calculation method [29, 30], we can find that before the power calculation we need to obtain I˙o-90°. In practical operation, the current waveform can be stored in electronic memory in real time. However, in the process of stability analysis, it is known that I˙o-90° is actually the result of the exchange between the real part and imaginary part of I˙o, so we can obtain(23)I˙o1-90=E1sinφ1R+Ro-RoE2sinφ2RR+2Ro-jE1cosφ1R+Ro-RoE2cosφ2RR+2Ro,(24)I˙o2-90=E2sinφ2R+Ro-RoE1sinφ1RR+2Ro-jE2cosφ2R+Ro-RoE1cosφ1RR+2Ro.

Combining (20), (23), and (24), we can get the traditional power calculation method:(25)pi=ai+jbi,qi=ci+jdi,where, ai, bi, ci, and di can be expressed by the inverter parameters which are given in Appendix C. Similarly, adopting the same method, a new method of power calculation can be obtained, which is the same as (25). And only the values of ai, bi, ci, and di are different, and they are given in Appendix D.

Using the droop method mentioned in (3):(26)Ei=Ei0-mp+mdsωcs+ωcpi,φi=φi0+nds2+nps+niωcss+ωcqi.

Differentiating (26) near the steady-state value E1e,E2e,φ1e,φ2e, we can obtain (27).(27)ΔE1Δφ1ΔE2Δφ2=-ωcmp+mdss+ωcΔp1ωcnds2+nps+niss+ωcΔq1-ωcmp+mdss+ωcΔp2ωcnds2+nps+niss+ωcΔq2,(28)Δp1Δq1Δp2Δq2=DΔφ1Δφ2ΔE1ΔE2,where D is the coefficient matrix:(29)D=d11d12d13d14d21d22d23d24d31d32d33d34d41d42d43d44.

Combining (28) and (29), we can obtain the small-signal model of the parallel system, as shown in (30). The coefficients in the equation are given in Appendix.

Similarly, we can construct the small-signal model based on the new power calculation method; the model is still as shown in (30), but the coefficients are slightly different; please refer to Appendix. (30)-d11ωcmp+mdss+ωc-d12ωcmp+mdss+ωc-d13ωcmp+mdss+ωc-1-d14ωcmp+mdss+ωcd21ss+ωcωcnds2+nps+ni-1d22ss+ωcωcnds2+nps+nid23ss+ωcωcnds2+nps+nid24ss+ωcωcnds2+nps+ni-d31ωcmp+mdss+ωc-d32ωcmp+mdss+ωc-d33ωcmp+mdss+ωc-d34ωcmp+mdss+ωc-1d41ss+ωcωcnds2+nps+nid42ss+ωcωcnds2+nps+ni-1d43ss+ωcωcnds2+nps+nid44ss+ωcωcnds2+nps+ni=0.

The model established by (30) is a six-order system, and its physical meaning is aimed at the whole parallel system while the traditional small-signal model established by (9) is aimed at the single inverter only. Although the model established by (30) is simpler than the model by Marwali et al., it still cannot be analyzed directly and is required to be solved by mathematical tools such as MATLAB.

5. Analysis of the New Small-Signal Model 5.1. Small-Signal Model of Traditional Power Calculation Method

First, the model established by the traditional power calculation method is analyzed, the circuit parameters are given in Table 1, and it is should be noted that the steady-state values should be satisfied (21). Figure 6 is a root locus of the small-signal model using the traditional power calculation method. (a) is the local root locus when ni is equal to 1 and mp increases from 0 to 1, and (b) is the local root locus when mp is equal to 1 and ni decreases from 1e-2 to 1e-5. Observing Figure 6, we discover that the parallel system, using traditional power calculation method, is not stable in the root locus, for which the traditional power calculation method will introduce an AC variable of double frequency. In order to suppress the oscillation, mp should be increased or ni should be decreased.

Parameters table.

 R/Ω 0.1 R o /Ω 48.4 E 1 e /V 220 E 2 e /V 220 φ 1 e - 1 e - 3 φ 2 e - 1 e - 3 ω c /rad/s 62.8

The root locus of the small-signal model based on the traditional power calculation method.

Figure 7 is a root locus of the small-signal model based on the traditional power calculation method after introducing differential control. The droop coefficients are as follows: ni=1e-3, mp=1e-1, np=1e-4, and md=1e-4, and nd increases from 0 to 8e-3. We can discover from the Table 1 that, after introducing differential control, a pair of conjugate complex roots are generated. This is the basic reason why the dynamic response of the model can be improved.

The root locus of the small-signal model based on the traditional power calculation method.

5.2. Small-Signal Model of the New Power Calculation Method

From the small-signal model, established by the traditional power calculation method, it can be found that the robustness of the parallel inverter is not strong. In order to improve it, the author uses the new power calculation method without a large delay low-pass filter, and this paper regards the model based on the new power calculation method as the main analysis object.

Figure 8 is the root locus of the small-signal model based on the new power calculation method only using proportional droop control. (a) is the root locus when ni is equal to 1e-4 and mp increases from 0 to 1e-2, and (b) is the root locus when mp is equal to 1e-4 and ni increases from 0 to 5e-3. Three important points can be obtained from Figure 8: firstly, the new power calculation method is better than the traditional power calculation method, and the poles always fall on the left side of the imaginary axis; secondly, when only the proportional droop is adopted, the poles fall on the negative real axis, which indicates that the dynamic performance is good; thirdly, with the increase of the proportional droop coefficient, the poles of the system shift to the left and the stability is improved.

The root locus of a small-signal model based on the new power calculation method.

Figure 9 is the root locus of the small-signal model based on the new power calculation method with differential control. In Figure 9, np=md=0, and the four curves correspond to four local root locus, respectively, when ni=mp=1e-1, 1e-2, 1e-3, and 1e-4 and nd increases from 0 to 1e-5. From Figure 9, we can discover that the introduction of differential control generates conjugate complex roots. And the significance of the conjugate complex roots is to improve the dynamic performance of the system, but the excessive increase of the differential coefficient will lead to the instability of the system. In addition, it can be found in Figure 9 that the increase of the proportional droop coefficient will lead the poles to move to the left and increase the stability of the system.

The root locus after a differential controller is used.

The control parameters of Figure 10(a) are as follows: mp=1e-2, ni=1e-2, nd=5e-7, and md=0, and np increases from 0 to 1e-4. The control parameters of Figure 10(b) are as follows: mp=1e-2, ni=1e-2, nd=5e-7, and np=5e-5, and md increases from 0 to 1e-4. Figure 10(a) indicates that np mainly affects the conjugate complex roots and can be modified to fine-tune the dynamic response of the system. md mainly affects the poles on the negative real axis, but this pole is not the main pole. So its influence on the system performance is not very large.

The root locus of the small-signal model based on the traditional power calculation method.

Figure 11 is the root locus when the impedance varies, and all the droop coefficients are given as follows: mp=1e-2, ni=1e-2, nd=5e-7, np=5e-5, and md=1e-8. (a) is the root locus when Ro increases from 1e-2 to 1e2 and (b) is the root locus when Ro increases from 1e-3 to 5e-1. (a) indicates that with the increase of resistance, that is, the load decrease, the poles on the negative real roots move to two sides of the axis, respectively, but the main poles are getting closer to the imaginary axis. So the stability of the system becomes worse when the load becomes smaller. (b) shows that, with a change of the output impedance of the inverter, the pole moves to the right but the main pole has not been changed. So a change of the output impedance has a limited effect on the parallel system.

The root locus when impedance varies.

Figure 12 is the root locus when the cut-off frequency of the low-pass filter increases from 1e-2 to 20, and the droop coefficients used are the same as those in Figure 11. Figure 12 shows that, to a certain extent, the system gets more stable if the cut-off frequency is increased, but the cut-off frequency cannot increase without limit. When the cut-off frequency increases to a certain value, the main pole of the system turns into a pair of conjugate complex roots that move right with the increase of the cut-off frequency. Then the system stability will get worse.

The root locus when ωc varies.

6. The Research on Application of the Model

Based on the small-signal model introduced in the this paper and the droop control that we researched before, the problem of unbalanced current and harmonic current caused by the unbalanced and nonlinear loads can be resolved and the output current can be shared equally. Thus, this part establishes the simulated and experimented platform to verify the theory. The simulated and experimental parameters are all within the range of previous stability analysis so that the correctness of the stability analysis can be well verified.

The simulated model is built by using Simulink/s-function. Power is expressed by per unit value and 1 means 11375 kVA. The simulated parameters are shown in Table 2.

Parameters of simulation and experiment.

 Sampling time Ts/s 1 e - 4 Switching frequency/Hz 5k Filter inductance L/mH 1.5 Parasitic inductance r/mΩ 270 Filter capacitor C/μF 20 Midline inductance Ln/mH 2 Rated voltage/V 220

First, the simulated results and analysis of the control strategy in the parallel inverter system with an unbalanced load are introduced as follows. Figure 13 is the unbalanced load simulated results waveform, which indicates that the unbalanced currents can be shared effectively.

The simulated results with an unbalanced load.

Second, the simulated results and analysis of the control strategy in the parallel inverter system with a nonlinear load are illustrated, and Figure 14 is the nonlinear load simulated waveforms, which indicates that the proposed droop method can share harmonic currents effectively.

The simulated results with a nonlinear load.

The experimental results and analysis of the parallel inverter system with an unbalanced load are described as follows. Figure 15 is the experimental unbalanced load waveforms, and we can discover that the unbalanced current can be accurately shared. The experimental results are consistent with the simulated results.

The experimental results with an unbalanced load.

Similarly, the experimental results and analysis of the nonlinear load are introduced. Figure 16 is the experimental waveforms with a nonlinear load, which indicates that the proposed droop method can share harmonic currents accurately and has good robustness.

The experimental results with a nonlinear load.

Figure 17 is the experimental waveforms with an unbalanced load and nonlinear load at the same time. It can be seen from the figure that the unbalanced and harmonic currents can both be effectively shared.

The experimental results with an unbalanced and nonlinear load.

Figure 18 is a photograph of the experimental platform. All the experiments in this paper were completed using this experimental platform.

The photograph of experiment platform.

7. Conclusions

This paper proposes a new small-signal modeling method for an auxiliary parallel inverter system and analyzes the influence of the parameters on the stability characteristics of the system in detail. The traditional analysis method is based on a small-signal model, but the model established by the traditional small-signal method is not accurate without considering the coupling relationship between parallel inverters, so this model cannot analyze the influence of the load and the output impedance of inverter. The small-signal modeling method proposed in this paper regards the whole parallel system as the research object, including the influences of the control system information outside the voltage loop and circuit parameters of inverter, which provide the most comprehensive analysis method by far.

Appendix A. The Coefficient Matrix <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M131"><mml:mrow><mml:mi>D</mml:mi></mml:mrow></mml:math></inline-formula> of Traditional Method

The coefficient matrix D of the traditional power calculation method is as follows:(A.1)d11=dp1dφ1=-2RoR+RoE12sin2φ1-RoRE1E2sinφ1+φ2RR+2Ro2+j2RoR+RoE12cos2φ1+RoRE1E2cosφ1+φ2RR+2Ro2,d21=dq1dφ1=2RoR+RoE12cos2φ1+RoRE1E2cosφ1+φ2RR+2Ro2+j2RoR+RoE12sin2φ1+RoRE1E2sinφ1+φ2RR+2Ro2,d31=dp2dφ1=-RoRE1E2sinφ1+φ2+2Ro2E12sin2φ1RR+2Ro2+jRoRE1E2cosφ1+φ2-2Ro2E12cos2φ1RR+2Ro2,d41=dq2dφ1=RoRE1E2cosφ1+φ2-2Ro2E12cos2φ1RR+2Ro2+jRoRE1E2sinφ1+φ2-2Ro2E12sin2φ1RR+2Ro2,d12=dp1dφ2=2Ro2E22sin2φ2-RoRE1E2sinφ1+φ2RR+2Ro2+j-2Ro2E22cos2φ2+RoRE1E2cosφ1+φ2RR+2Ro2,d22=dq1dφ2=-2Ro2E22cos2φ2+RoRE1E2cosφ1+φ2RR+2Ro2-j2Ro2E22sin2φ2-RoRE1E2sinφ1+φ2RR+2Ro2,d32=dp2dφ2=-RoRE1E2sinφ1+φ2-2RoR+RoE22sin2φ2RR+2Ro2+jRoRE1E2cosφ1+φ2+2RoR+RoE22cos2φ2RR+2Ro2,d42=dq2dφ2=RoRE1E2cosφ1+φ2+2RoR+RoE22cos2φ2RR+2Ro2+jRoRE1E2sinφ1+φ2+2RoR+RoE22sin2φ2RR+2Ro2,d13=dp1dE1=2RoR+RoE1cos2φ1+RoRE2cosφ1+φ2RR+2Ro2+j2RoR+RoE1sin2φ1+RoRE2sinφ1+φ2RR+2Ro2,d23=dq1dE1=2RoR+RoE1sin2φ1+RoRE2sinφ1+φ2RR+2Ro2-j2RoR+RoE1cos2φ1+RoRE2cosφ1+φ2RR+2Ro2,d33=dp2dE1=RoRE2cosφ1+φ2-2Ro2E1cos2φ1RR+2Ro2+jRoRE2sinφ1+φ2-2Ro2E1sin2φ1RR+2Ro2,d43=dq2dE1=RoRE2sinφ1+φ2-2Ro2E1sin2φ1RR+2Ro2-jRoRE2cosφ1+φ2-2Ro2E1cos2φ1RR+2Ro2,d14=dp1dE2=-2Ro2E2cos2φ2+RoRE1cosφ1+φ2RR+2Ro2+j-2Ro2E2sin2φ2+RoRE1sinφ1+φ2RR+2Ro2,d24=dq1dE2=-2Ro2E2sin2φ2+RoRE1sinφ1+φ2RR+2Ro2-j-2Ro2E2cos2φ2+RoRE1cosφ1+φ2RR+2Ro2,d34=dp2dE2=RoRE1cosφ1+φ2+2RoR+RoE2cos2φ2RR+2Ro2+jRoRE1sinφ1+φ2+2RoR+RoE2sin2φ2RR+2Ro2,d44=dq2dE2=RoRE1sinφ1+φ2+2RoR+RoE2sin2φ2RR+2Ro2-jRoRE1cosφ1+φ2+2RoR+RoE2cos2φ2RR+2Ro2.

B. The Coefficient Matrix <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M134"><mml:mrow><mml:mi>D</mml:mi></mml:mrow></mml:math></inline-formula> of New Method

The coefficient matrix D of the new power calculation method is as follows:(B.1)d11=dp1dφ1=-RoRE1E2sinφ1-φ2RR+2Ro2,d21=dq1dφ1=-RoR+2Ro2E1E2cosφ1-φ2RR+2Ro2,d31=dp2dφ1=-RoRE1E2sinφ1-φ2RR+2Ro2,d41=dq2dφ1=RoR+2Ro2E1E2cosφ1-φ2RR+2Ro2,d12=dp1dφ2=RoRE1E2sinφ1-φ2RR+2Ro2,d22=dq1dφ2=RoR+2Ro2E1E2cosφ1-φ2RR+2Ro2,d32=dp2dφ2=RoRE1E2sinφ1-φ2RR+2Ro2,d42=dq2dφ2=-RoR+2Ro2E1E2cosφ1-φ2RR+2Ro2,d13=dp1dE1=2RoR+RoE1+RoRE2cosφ1-φ2RR+2Ro2,d23=dq1dE1=-RoR+2Ro2E2sinφ1-φ2RR+2Ro2,d33=dp2dE1=RoRE2cosφ1-φ2-2Ro2E1RR+2Ro2,d43=dq2dE1=RoR+2Ro2E2sinφ1-φ2RR+2Ro2,d14=dp1dE2=-2Ro2E2+RoRE1cosφ1-φ2RR+2Ro2,d24=dq1dE2=-RoR+2Ro2E1sinφ1-φ2RR+2Ro2,d34=dp2dE2=RoRE1cosφ1-φ2+2RoR+RoE2RR+2Ro2,d44=dq2dE2=RoR+2Ro2E1sinφ1-φ2RR+2Ro2.

C. The Coefficients of Traditional Method

The coefficients of traditional power calculation method are given as follows: (C.1)a1=RoR+RoE12cos2φ1-Ro2E22cos2φ2+RoRE1E2cosφ1+φ2,b1=RoR+RoE12sin2φ1-Ro2E22sin2φ2+RoRE1E2sinφ1+φ2,c1=RoR+RoE12sin2φ1-Ro2E22sin2φ2+RoRE1E2sinφ1+φ2,d1=-RoR+RoE12cos2φ1+Ro2E22cos2φ2-RoRE1E2cosφ1+φ2,a2=RoRE1E2cosφ1+φ2+RoR+RoE22cos2φ2-Ro2E12cos2φ1,b2=RoRE1E2sinφ1+φ2+RoR+RoE22sin2φ2-Ro2E12sin2φ1,c2=RoRE1E2sinφ1+φ2+RoR+RoE22sin2φ2-Ro2E12sin2φ1,d2=-RoRE1E2cosφ1+φ2-RoR+RoE22cos2φ2+Ro2E12cos2φ1.

D. The Coefficients of New Method

The coefficients of new power calculation method are given as follows: (D.1)p~1=RoR+RoE12-Ro2E22+RoRE1E2cosφ1-φ2RR+2Ro2,q~1=-RoR+2Ro2E1E2sinφ1-φ2RR+2Ro2,p~2=RoRE1E2cosφ1-φ2+RoR+RoE22-Ro2E12RR+2Ro2,q~2=RoR+2Ro2E1E2sinφ1-φ2RR+2Ro2.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Basic Scientific Research Project of Central University (Grant no. 2016RC038) and National Key R&D Plan (Grant no. 2016YFB1200502).

Shafiee Q. Guerrero J. M. Vasquez J. C. Distributed secondary control for islanded microgrids-a novel approach IEEE Transactions on Power Electronics 2014 29 2 1018 1031 2-s2.0-84883235871 10.1109/TPEL.2013.2259506 Jiefeng H. Jianguo Z. Dorrell D. G. Guerrero J. M. Virtual flux droop method - a new control strategy of inverters in microgrids IEEE Transactions on Power Electronics 2014 29 9 4704 4711 10.1109/TPEL.2013.2286159 2-s2.0-84900413413 Yao W. Chen M. Matas J. Guerrero J. M. Qian Z.-M. Design and analysis of the droop control method for parallel inverters considering the impact of the complex impedance on the power sharing IEEE Transactions on Industrial Electronics 2011 58 2 576 588 2-s2.0-78651351020 10.1109/TIE.2010.2046001 Shi Y. Wu W. Wang H. Du Y. Su J. The Parallel Multi-Inverter System Based on the Voltage-Type Droop Control Method IEEE Journal of Emerging and Selected Topics in Power Electronics 2016 4 4 1332 1341 2-s2.0-84995475779 10.1109/JESTPE.2016.2581018 Guerrero J. M. Vasquez J. C. Matas J. de Vicuña L. G. Castilla M. Hierarchical control of droop-controlled AC and DC microgrids: a general approach toward standardization IEEE Transactions on Industrial Electronics 2011 58 1 158 172 10.1109/TIE.2010.2066534 2-s2.0-78650181509 Wei Y. Min C. Matas J. Guerrero J. M. Qian Z.-M. Design and analysis of the droop control method for parallel inverters considering the impact of the complex impedance on the power sharing IEEE Transactions on Industrial Electronics 2011 58 2 576 588 10.1109/TIE.2010.2046001 2-s2.0-78651351020 Chen J. Wang L. Diao L. Du H. Liu Z. Distributed auxiliary inverter of urban rail train-load sharing control strategy under complicated operation condition IEEE Transactions on Power Electronics 2016 31 3 2518 2529 2-s2.0-84949554598 10.1109/TPEL.2015.2427381 Ming H. Haibing H. Xing Y. Guerrero J. M. Multilayer control for inverters in parallel operation without intercommunications IEEE Transactions on Power Electronics 2012 27 8 3651 3663 10.1109/TPEL.2012.2186985 2-s2.0-84860274250 Wei Y. Chen M. Matas J. Guerrero J. M. Qian Z.-M. Design and analysis of the droop control method for parallel inverters considering the impact of the complex impedance on the power sharing IEEE Transactions on Industrial Electronics 2011 58 2 576 588 10.1109/TIE.2010.2046001 2-s2.0-78651351020 Guerrero J. M. Vasquez J. C. Matas J. Castilla M. de Vicuña L. G. Control strategy for flexible microgrid based on parallel line-interactive UPS systems IEEE Transactions on Industrial Electronics 2009 56 3 726 736 10.1109/TIE.2008.2009274 2-s2.0-63149115320 Wu T. Liu Z. Liu J. Wang S. You Z. A Unified Virtual Power Decoupling Method for Droop-Controlled Parallel Inverters in Microgrids IEEE Transactions on Power Electronics 2016 31 8 5587 5603 2-s2.0-84963741872 10.1109/TPEL.2015.2497972 Vasquez J. C. Guerrero J. M. Savaghebi M. Eloy-Garcia J. Teodorescu R. Modeling, analysis, and design of stationary-reference-frame droop-controlled parallel three-phase voltage source inverters IEEE Transactions on Industrial Electronics 2013 60 4 1271 1280 10.1109/tie.2012.2194951 2-s2.0-84870490332 Antonio E. Coelho A. Cortizo P. C. Francisco P. Garcia D. Small signal stability for single phase inverter connected to stiff AC system Proceedings of the 1999 IEEE Industry Applications Conference - 34th IAS Annual Meeting October 1999 2180 2187 2-s2.0-0033331536 Cheng Y. Zha X. Chen H. The small signal model for wireless parallel operation of inverters Proceedings of the 2009 IEEE International Conference on Mechatronics and Automation, ICMA 2009 August 2009 chn 433 437 2-s2.0-77449147689 10.1109/ICMA.2009.5246748 Gao F. Li Z. Li Y. Wang P. Zhu H. Small-signal stability analysis of parallel-connected inverters based on time-varying phasor Proceedings of the 2011 6th IEEE Conference on Industrial Electronics and Applications, ICIEA 2011 June 2011 chn 1239 1244 2-s2.0-80052238298 10.1109/ICIEA.2011.5975775 Coelho E. A. A. Cortizo P. C. Garcia P. F. D. Small-signal stability for parallel-connected inverters in stand-alone ac supply systems IEEE Transactions on Industry Applications 2002 38 2 533 542 2-s2.0-0036505184 10.1109/28.993176 Marwali M. N. Jung J.-W. Keyhani A. Stability analysis of load sharing control for distributed generation systems IEEE Transactions on Energy Conversion 2007 22 3 737 745 2-s2.0-34548260882 10.1109/TEC.2006.881397 Hua M. Hu H. Xing Y. Guerrero J. M. Multilayer control for inverters in parallel operation without intercommunications IEEE Transactions on Power Electronics 2012 27 8 3651 3663 2-s2.0-84860274250 10.1109/TPEL.2012.2186985 Chiang S. J. Chang J. M. Parallel control of the UPS inverters with frequency-dependent droop scheme Proceedings of the 2001 IEEE 32nd Annual Power Electronics Specialists Conference June 2001 can 957 961 2-s2.0-0034785063 Ning B. Tang T. Dong H. Wen D. Liu D. Gao S. Wang J. An introduction to parallel control and management for high-speed railway systems IEEE Transactions on Intelligent Transportation Systems 2011 12 4 1473 1483 2-s2.0-82455199264 10.1109/TITS.2011.2159789 Xu S. Wang J. Xu J. A current decoupling parallel control strategy of single-phase inverter with voltage and current dual closed-loop feedback IEEE Transactions on Industrial Electronics 2013 60 4 1306 1313 2-s2.0-84870473457 10.1109/TIE.2011.2161660 Wei X. Dai K. Fang X. Geng P. Luo F. Kang Y. Parallel Control of Three-Phase Three-Wire Shunt Active Power Filters Proceedings of the 2006 CES/IEEE 5th International Power Electronics and Motion Control Conference August 2006 Shanghai, China 1 5 10.1109/IPEMC.2006.4778099 Kong Q.-J. Li L. Yan B. Lin S. Zhu F. Xiong G. Developing parallel control and management for urban traffic systems IEEE Intelligent Systems 2013 28 3 66 69 2-s2.0-84884133404 10.1109/MIS.2013.59 Guerrero J. M. de Vicuña L. G. Matas J. Castilla M. Miret J. A wireless controller to enhance dynamic performance of parallel inverters in distributed generation systems IEEE Transactions on Power Electronics 2004 19 5 1205 1213 10.1109/TPEL.2004.833451 2-s2.0-4544344304 Wei X. Zhu G. Lu J. Xu X. Instantaneous current-sharing control scheme of multi-inverter modules in parallel based on virtual circulating impedance IET Power Electronics 2016 9 5 960 968 2-s2.0-84966333797 10.1049/iet-pel.2014.0932 Lei Q. Peng F. Z. Yang S. Multiloop control method for high-performance microgrid inverter through load voltage and current decoupling with only output voltage feedback IEEE Transactions on Power Electronics 2011 26 3 953 960 2-s2.0-79957439639 10.1109/TPEL.2010.2096992 Ming H. Hu H. Xing Y. Guerrero J. M. Multilayer control for inverters in parallel operation without intercommunications IEEE Transactions on Power Electronics 2012 27 8 3651 3663 10.1109/TPEL.2012.2186985 2-s2.0-84860274250 Golestan S. Monfared M. Guerrero J. M. Joorabian M. A D-Q synchronous frame controller for single-phase inverters Proceedings of the 2011 2nd Power Electronics, Drive Systems and Technologies Conference, PEDSTC 2011 February 2011 irn 317 323 2-s2.0-79955678868 10.1109/PEDSTC.2011.5742439 Ren Z. Gao M. Mo Q. Liu K. Yao W. Chen M. Qian Z. Power calculation method used in wireless parallel inverters under nonlinear load conditions Proceedings of the 25th Annual IEEE Applied Power Electronics Conference and Exposition, APEC 2010 February 2010 usa 1674 1677 2-s2.0-77952230487 10.1109/APEC.2010.5433456 Sutanto D. Outhred H. R. Lee Y. B. Probabilistic power system production cost and reliability calculation by the z-transform method IEEE Transactions on Energy Conversion 1989 4 4 559 566 2-s2.0-0024943690 10.1109/60.41714