# Harmonic Modeling of a Diode-Clamped Multilevel Voltage Source Converter for Predicting Uncharacteristic Harmonics

Authors: Tian Hao Huang; Kuo Lung Lian
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In response to the growing demand for medium- and high-power trends, multilevel voltage source converters (VSCs) have been attracting growing considerations. One of the widely used VSCs are the diode-clamped multilevel VSC (DCM-VSC). As these converters proliferate, their harmonic impact may become significant. Nevertheless, a harmonic model for the DCM-VSC is currently lacking in the literature. In this paper, the ABCD matrix, mapping the input harmonics to the output harmonics of DCM-VSC, is derived. As the matrix is formulated in the time-domain, the output harmonics are exact and do not suffer from harmonic truncation errors. As the paper will demonstrate, the derived ABCD matrix can be easily applied to a microgrid system and users can easily predict all the uncharacteristic harmonics when a microgrid is subjected to various conditions of imbalance. In addition to all the results being validated with those of PSCAD/EMTDC, the computation time of the proposed method is in contrast much shorter.


---

## Body

## 1. Introduction

Voltage source converters (VSCs) have been widely used in the modern power system [1]. They are the main interconnection devices for distributed generators (DGs) and energy storage systems in a microgrid system [2]. At the transmission and distribution levels, they are also the main building blocks of the flexible ac transmission system devices and custom power controllers, respectively. While many benefits may be realized from proliferation of VSCs in the power grid, these come at the expense of generating harmonic distortion. It is therefore essential to obtain accurate harmonic models for VSCs to predict their generated current and voltage harmonic under various conditions and to understand how these harmonics will interact with the rest of the power system.Electromagnetic transient programs such as PSCAD/EMTDC can provide accurate harmonic analysis results. However, its disadvantage is that an initialization transient must die out before steady-state is reached when Fourier analysis can be performed [3]. The simulation time step must be sufficiently small for the VSCs to capture the switching dynamics, and this leads to excessively long simulation times. Alternatively, one can rely on harmonic power flow to obtain the harmonic solutions. However, a frequency map (FM) between converters’ input and output harmonics usually requires deriving. Hence, the main objective of this paper is to derive the FM of a multilevel VSC.Depending on the number of the voltage levels, VSCs can be classified into two-level and multilevel (more than two-level) converters. The multilevel VSCs have the advantages of reduced voltage stress on power semiconductor switches, lower distortions of the input current, and lower dv/dt, as compared with their two-level counterparts. Consequently, multilevel VSCs have become attractive even for the low- and medium-power range such as in the distributed energy conversion application [4]. Among many various types of multilevel VSCs, diode-clamped multilevel VSCs (DCM-VSCs) are the most widely used for utility and high-power industry applications [5] ever since their first appearance in 1981 [6]. As these converters become more dispersed, their harmonic impact on the system could be significant. There have been numerous publications [7–9] addressing harmonic modeling of a two-level VSC. Nevertheless, a generalized harmonic model of the DCM-VSC has not been developed yet. Thus, the focus of this paper is on the DCM-VSCs.As VSCs may be subjected to various power quality disturbance, the derived harmonic models of DCM-VSCs should be general and able to account for various conditions such as unbalanced grid voltage [10, 11] and unbalanced dc voltages [12, 13]. The imbalance in the dc voltage comes from (i) the DGs with different voltage levels connected to the dc link of the converters [14, 15] and (ii) the inherent dc voltage imbalance of DCM-VSCs [16]. It will be demonstrated in the paper that the two cases will lead to drastically different harmonic profiles. The first dc voltage imbalance may result in low-order even current harmonics at the ac side [14], which in turn may cause undesired dc component and unequal positive and negative peak values [17]. However, in the second case, no even harmonics are generated, but some uncharacteristic odd harmonics are generated. Moreover, as will be demonstrated in the paper, the harmonic profile of converters will also alter significantly in a microgrid environment or in the presence of neighboring converters. Therefore, it is essential to have an accurate harmonic model, capable of predicting uncharacteristic harmonics under various conditions.The current models for DCM-VSCs are inadequate for harmonic modeling. Khajehoddin et al. [4] proposed a current flow model for the DCM-VSCs. They formulated the capacitor currents as well as the input dc current as a function of the switching functions and output currents. Nevertheless, they ignored the high-order harmonics of the switching functions, which may result in inaccurate side band harmonic arising with converter modulation. Generalized averaged models of DCM-VSCs have been proposed in [18]. However, all the high-order harmonics are effectively eliminated in averaging operations, thereby producing inaccurate harmonic results. Jayaram et al. [19] derived a five-level DCM-VSC in abc frame for which the switching functions are obtained using curve fitting. However, such a model becomes very complex when the number of level is high. He et al. [20] have modeled the DCM-VSC as a Norton equivalent circuit. However, they have assumed that the capacitors are large enough so that the dc link voltage can be assumed to be constant. Such models are only accurate in predicting characteristic harmonics and are not able to account for the uncharacteristic harmonics under unbalanced conditions.Based on the work of Larsen et al. [21], the FM of an AC-DC converter is more conveniently modeled by means of an ABCD matrix. A DCM-VSC can be visualized in Figure 1. The harmonic injections caused by external networks can be modeled as ac voltage harmonics Vac (in this paper, capital letters represent harmonic domain variables, whereas lower-case letters represent time-domain variables) and dc current harmonics Idc at the point of common coupling. The converter in turn produces ac current harmonics and dc voltage harmonics Vdc via the ABCD matrix as stated in (1).(1)IacVdc=ABCDVacIdc.Figure 1
Equivalent external network of DCM-VSC.Therefore, the main goal of this paper is to derive an accurate ABCD matrix for a DCM-VSC under various conditions. The proposed method essentially consists of two steps: the first is to derive the state-space equations of the DCM-VSC, and the second is to treat the harmonics of interest as state variables and solve the entire augmented state-space equations piecewise with the periodicity constraint imposed. Then the ABCD matrix is obtained by mapping the resulting submatrix of the state transition matrix to the harmonic domain by a transformation matrix. The advantage of the proposed method is that(i)
calculation times are much shorter than the brute-force time-domain method as it directly bypasses the transient state;(ii)
different from the harmonic domain methods, there is no need to know how many harmonics need to be included to obtain the harmonics of interest;(iii)
the predicted harmonics are exact and do not suffer from harmonic truncation errors [7];(iv)
it can be easily expanded to study a microgrid system containing multiple numbers of DCM-VSCs.To validate the proposed model, all the predicted harmonics are verified with those obtained from PSCAD/EMTDC.The rest of the paper is organized as follows: Section2 describes the state-space model of a DCM-VSC. Section 3 shows how the harmonic ABCD coupling matrix is derived. In Section 4, various scenarios of imbalance for a DCM-VSC and DCM-based microgrid are studied, and the predicted results are verified with those of PSCAD/EMTDC, demonstrating the validity of the proposed method. Finally, a conclusion is given in Section 5.


## 2. Generalized State-Space Representation of a DCM-VSC

Different from other conventional multilevel converters such as flying capacitors or cascaded multilevel converters, DCM-VSCs have forbidden states [22], in which the output voltage is different if the current direction is different under the same switching function state. Figure 2 shows the three-level VSC where the upper switches are assumed to be complementary with the lower switches in each leg.Figure 2
Three-phase, three-level DCM-VSC.To illustrate the forbidden states, Figure3 is referred to, and without loss of generality, only phase A leg of the VSC is shown. When s11 turns on and s12 turns off and if the current flows as shown in Figure 3(a), vao=-vc2. However, when the current flows in the reverse direction, the circuit becomes Figure 3(b) and vao=vc1. Consequently s11=1 and s12=0 constitute a forbidden state because the output voltage is different when the current direction is different. In contrast, when s11=0, s12=1, the circuit remains that shown in Figure 4, regardless of the direction of current flow. Such as a state is a legitimate state. Table 1 lists all the possible legitimate switching states for phase A.Table 1
Legitimate switching states of phase A.
s11 s 12 vao 0 0 vC1  0 1 0  1 1 -vC2Figure 3
DCM-VSC forbidden states: (a) Left: current out of converter; (b) Right: current into the converter.Figure 4
DCM-VSC legitimate states: (a) current out of converter; (b) current into the converter.
(a) (b)To have a better representation of all these states, Figure2 is replaced with Figure 5 where all the bidirectional switches are replaced with single-pole triple-throw switches [23]. sa in Figure 5 can be represented by three switching states, sa1, sa2, and sa3. When sa is connected to ea1 in Figure 5, sa1=1, sa2=0, and sa3=0. Table 2 summarizes the switching status for sa. Similar rules apply to sb and sc.Table 2
Switching status ofsa.
s a 1 s a 2 sa3 s a → e a 1 1 0 0  s a → e a 2 0 1 0  s a → e a 3 0 0 1Figure 5
The equivalent single-pole, triple-throw model for a three-level DCM-VSC.Based on Table2, one can derive a generalized m-level (m is assumed to be odd.), three-phase, four-wire model (see Figure 6) of DCM-VSC with parameters such as inductance and resistance for each phase that may not necessarily be the same. In addition, without loss of generality, the stimuli are assumed to take the form of unej(nωt+θn), where n is the number of harmonics and θn the phase angle at the n-th harmonic.(2)ddtxu=FG0Ωxuwhere(3)x=iaibicvC1cC2⋯vCm-1Tu=vsavsbvsciLTSubmatrices, F, G, and Ω are stated in the Appendix. Moreover, the derivation of a simple three-level DCM-VSC can also be found in [24].Figure 6
The equivalent single-pole, m-throw model for a three-phase, four-wire m-level DCM-VSC.


## 3. Derivation of the ABCD Matrix for Harmonic Evaluation

Equation (2) is a time-domain state-space representation. To step into the harmonic domain, one can treat the harmonics of interest as harmonic states and augment them to (2) by recasting a Fourier integral into a differential equation [7] as stated in (4).(4)dxhdt=jnωxh+1Te-jnωTxt,where xh is the harmonic state, T is the period, and xh(0)=0.If one wants to know the harmonics for the three-phase current and all them-1 capacitor voltages, the overall equation is as stated in (5).(5)ddtxuxh=F˘xuxhwhere(6)F˘=FG00Ω0D0JD=D1D2⋮D2+mD1=ejnωTT0⋯0ejn-1ωTT0⋯0⋮⋮⋮⋮1T0⋯0⋮⋮⋮⋮e-jn-1ωTT0⋯0e-jnωTT0⋯02∗n+1×2+mD2+m=0⋯0ejnωTT0⋯0ejn-1ωTT⋮⋮⋮⋮0⋯01T⋮⋮⋮⋮0⋯0e-jn-1ωTT0⋯0e-jnωTT2∗n+1×2+mJ=diagJ1,J2,⋯,J2+mJ1=J2=⋯=J2+m=diag-jnω,-jn-1ω,⋯,0,⋯,jn-1ω,jnωEvaluating (5) over T, one gets (7).(7)xTuTxhT=Φˇx0u0xh0where(8)Φˇ=eFˇnswT-tnsw⋯eFˇ2t2-t1eFˇ1t1-0,nsw is the number of switching times.As proved in [7], Φˇ in general takes the following form.(9)Φˇ=αβ00γ0ςηξImposing the periodicity constraints(10)xt+T=xt,where t is generally set to 0 and the values of xh(T) are essentially the harmonics of interests [7]. In other words, the output harmonics, such as AC current harmonics Iabc and DC voltage harmonics Vdc, are equal to xh(T) and can be obtained analytically as stated in (11).(11)IabcVdc=Mvsabc0idc0where(12)M=ςI-α-1β+ηThe matrix M in (11) relates the steady-state initial conditions of the ac voltages and dc currents (i.e., vsabc(0) and idc(0)) to the ac current harmonics Iac and dc voltage harmonics Vdc. Because the input is the time-domain variable, one needs to find the relation between input steady-state initial conditions (vsabc(0) and idc(0))) and input harmonics (i.e., Vsabc and Idc). As the stimuli are in the form of unej(nωt+θn),(13)VsabcIabc=vsabc0idc0.Hence,(14)IabcVdc=ABCDabcVsabcIdc,where ABCDabc=[M]. Alternatively, one can also formulate (14) in terms of the sequential components (zero, positive, and negative sequence components are designated as 0, 1, and 2, respectively). Since(15)Iabc=P·I012,(16)Vsabc=P·Vs012where(17)P=1111a2a1aa2a=1∠120∘.Substituting (15) into (16) and converting harmonic phase currents Iabc into harmonic sequential currents I012, one obtains(18)P00II012Vdc=ABCDabcP00IVa012Idc(19)I012Vdc=ABCD012V012Idc,where ABCD012=P00I-1ABCDabcP00I and I is the identity matrix.


## 4. Case Studies

To validate the uncharacteristic harmonics predicted by the proposed method, five case studies are investigated. In the first case, the characteristic harmonics produced by a three-level DCM-VSC under a balanced ac grid are first evaluated. Then, the uncharacteristic harmonics due to the negative sequence ac voltage injection are studied. The second case is to investigate the uncharacteristic harmonics of a three-level DCM-VSC with unequal dc voltage levels, caused by external DGs. The third case is to show that the proposed model can reproduce the inherent unbalanced dc voltage of a five-level DCM-VSC. In the fourth case, the proposed model is shown to be able to predict the uncharacteristic harmonics arising from an imbalance in system parameters. The converter parameters are extracted from [25] and listed in Table 3. Note that for Cases 1 to 3, RN and LN are assumed to be infinity. Thus, a three-phase, three-wire DCM-VSC system is investigated. In contrast, RN and LN for Case 4 are finite, and a three-phase, four-wire DCM-VSC system is investigated. Finally, in case 5, a microgrid consisting of the converters presented in Cases 2, 3, and 4 is modeled using the proposed method.Table 3
The parameter list of the three-level DCM-VSC.
Parameter Value Parameter Value R a,Rb,Rc 0.1Ω firing angle,δ 0 ∘  L a,Lb,Lc 0.004H modulation index,ma 0.833  f 60 Hz modulation factor,mf 15  C 1,C2 1000μF source voltage 380VThe modulation strategy for multilevel converters can be classified into three types, specifically, carrier-based pulse width modulation (PWM), space vector modulation, and selective harmonic elimination. Among them, the carrier-based modulation strategy is widely used in industry [26]. The carrier-based PWM can be further classified into phase-shifted PWM (PS-PWM) and level-shifted PWM (LS-PWM). PS-PWM is rarely used for DCM-VSC as it has no modular structure and the carrier cannot be associated with a particular cell or treated independently [27]. The LS-PWM, on the other hand, can be easily adapted to a DCM. There are three types of LS-PWMs [28]—(i) in-phase disposition (IPD), where all carrier are in phase; (ii) alternative phase opposite disposition (APOD), where all carriers are alternatively in opposite disposition; and (iii) phase opposite disposition (POD), where all the positive carriers are in phase with each other and in opposite phase of the negative carriers. For the case studies, the IPD-LS-PWM is used because of its better power quality [28]. The switching times are obtained by solving the transcendental equations of the modulating and carrier signals.


### 4.1. Case1: Balanced and Unbalanced AC Grid (Negative Sequence Voltage Injection)

The three-level DCM-VSC is first connected to a balanced AC grid. Figure7 shows the maximum difference (Id) of the ac currents obtained by the proposed method and those by PSCAD/EMTDC under various time steps. As the figure shows, the simulation step tstep needs to reduce to 0.1μs in order for PSCAD/EMTDC model to reach the same level of accuracy as the proposed method. Figure 8 shows the resulting balanced three-phase currents for both PSCAD/EMTDC (tstep=0.1μs) and the proposed method. Nevertheless, the calculation time for PSCAD/EMTDC with time step of 0.1μs is about 85 times longer than the proposed method. Figure 9 shows the current harmonic spectra in terms of its positive and negative sequence components. It clearly shows that under balanced condition, the DCM-VSC only produces characteristic harmonics.Figure 7
I d versus solution time step.Figure 8
Balanced steady-state AC current waveforms.Figure 9
(a) The harmonic spectra of ac currents of Figure8 in terms of positive sequence components; (b) the harmonic spectra of ac currents of Figure 8 in terms of negative sequence components.
(a) (b)It is to be noted that input stimuli in (2) can be easily extended to other stimuli such as the cosine function by means of Euler’s formula, which is stated in (20).(20)cos⁡ωt+θ=ejωt+θ+e-jωt+θ2The following, as an example, are the step-by-step procedures for calculating the positive sequence of the fundamental current and dc voltage harmonic under the balanced condition:(1)
The input stimuli for phases A, B, and C are{(vs/2)e-jωt+θ,(vs/2)ej(ωt+θ)}, {(vs/2)e-j(ωt+θ-2π/3),(vs/2)ej(ωt+θ-2π/3)}, and {(vs/2)e-j(ωt+θ+2π/3),(vs/2)ej(ωt+θ+2π/3)}. Note that vs is the magnitude of the source voltage.(2)
Substitute the parameter values listed in Table3 into the form, similar to (5)(3)
Evaluate (7) to obtain Φˇ.(4)
Extract the submatrices,A, B, C, and D, which are as follows.(21)A=1.6490+0.0489j1.6709-0.2624j-0.7822-1.4525j-1.0626+1.5416j-0.8668+1.4036j-0.6083-1.2792j-0.7822-1.4525j-0.6083-1.2792j-0.8668+1.4036j1.6709+0.2624j1.6490+0.0489j-1.626+1.5416j-0.8668+1.4036j-1.0626+1.5416j1.6490+0.0489j-0.6083-1.2792j-0.7822-1.4525j1.6709-0.2624j(22)B=0.1887+0.0148j-0.1072+0.1560j-0.0815-0.1708j(23)C=5.8125-0.3137j5.8125+0.3137j-3.178-4.8769j-3.178+4.8769j-2.6345+5.1907j-2.6345-5.1907j5.8125-0.3137j5.8125+0.3137j-3.178-4.8769j-3.178+4.8769j-2.6345+5.1907j-2.6345-5.1907j(24)D=00(5)
The ac current and dc voltage harmonics can be obtained by evaluating (14). Note that the ac current harmonics need to be multiplied by two to obtain the correct phasor forms(25)Ia1Ib1Ic1Vc10Vc20=200000200000200000100000122.76-0.34j11.09+19.88j11.67-19.54j292.03292.03=-45.53-0.68j22.18+39.77j23.35-39.09j292.03292.03,whereIa1, Ib1, and Ic1 denote the fundamental currents of phases A, B, and C, respectively, and Vc10 and Vc20 are the two dc voltages across C1 and C2(6)
Evaluate (19), and the fundamental positive (I11), negative (I21), and zero sequence (I01) currents and the two dc voltages are stated as (26).(26)I01I11I21Vc10Vc20=045.53-0.68j0292.03292.03Next, a5% negative sequence of the fundamental voltage is added to the AC grid, and the voltage takes the following form.(27)Vsa=Vsejωt+0.05VsejωtVsb=Vsejωt-120∘+0.05Vsejωt+120∘Vsc=Vsejωt+120∘+0.05Vsejωt-120∘Figure10 shows the harmonic spectra of this unbalanced case in terms of positive and negative sequence components. In this instance, the uncharacteristic harmonics such as the negative sequence fundamental and third harmonics are accurately predicted by the proposed method.Figure 10
Unbalanced ac grid case: (a) the harmonic spectra of ac currents of in terms of positive sequence components; (b) the harmonic spectra of ac currents of in terms of negative sequence components.
(a) (b)


### 4.2. Case2: Unequal DC Voltage due to DGs with Different Voltage Levels

In this case, we try to simulate the case where the split dc link capacitors of the DCM-VSC can be controlled independently. A commonly seen situation is to have a different PV panel connected to each capacitor, with each panel being controlled independently to achieve its own point of maximum power [14]. Such a case leads to unequal dc voltage level for different capacitors with respect to the neutral point, N. To simulate such a case, two DC sources with vDC1=350 V connected in parallel with C1 and vDC1=230 V connected in parallel with C2 are connected to the dc link. The resulting current spectra, shown in Figure 11, feature undesired negative sequence even-order current harmonics such as the second and 14-th harmonics. The values of these harmonics are accurately predicted by the proposed method when compared with those of PSCAD/EMTDC.Figure 11
Case 2: (a) the harmonic spectra of ac currents in terms of positive sequence components; (b) the harmonic spectra of ac currents in terms of negative sequence components.
(a) (b)


### 4.3. Case3: Inherent Unbalanced Capacitor Voltage for A 5-Level DCM-VSC

Different from the previous case, this case study investigates a five-level DCM-VSC. The parameter settings are the same as those in Table3. As stated in [16, 29], when the level number is greater than 3, an inherent dc voltage imbalance becomes apparent. In addition, as demonstrated in [22], if the condition(28)macos⁡δ>0.55is satisfied, no balancing strategy exists to balance the capacitor voltages. Hence, for the chosen ma and δ in Table 3, voltage imbalance among capacitors is expected.Figures12(a) and 12(b) show the DC voltage waveforms for each capacitor (in this paper, the capacitor labels follow those indicated in Figure 6), obtained by the proposed method and PSCAD/EMTDC. Their corresponding harmonic spectra are in great agreement as shown in Figure 13. Figure 14 shows the current waveforms of the five-level DCM-VSC. Figure 15 shows the corresponding current spectra. Compared with Figure 9, one can see that the positive sequence 7-th, 13-th, 19-th, and 25-th harmonics are much greater for the 5-level DCM-VSC. Moreover, for the negative sequence, the value of the 5-th harmonic is about 5 times larger than that of 11-th harmonic for the 5-level case, whereas the 11-th harmonic for the 3-level converter is 3 times higher than that of the 5-th harmonic as seen in Figure 9. Nevertheless, different from Case 2, there is no even current harmonic generated at the ac side. The reason is that, with respect to N, there exists a symmetry between the upper and lower capacitor voltages. In other words, vc1=vc4 and vc2=vc3. Consequently, such a symmetry with respect to N will not generally lead to any even current harmonics at the ac side.Figure 12
Case 3: (a) the capacitor voltage waveforms obtained by the proposed method; (b) the capacitor voltage waveforms obtained by the PSCAD/EMTDC.
(a) (b)Figure 13
Harmonic spectra of each capacitor in the five-level DCM-VSC.Figure 14
Steady-state AC current waveforms of five-level DCM-VSC.Figure 15
Case 3: (a) the current harmonic spectra of five-level DCM-VSC in terms of positive sequence components; (b) the current harmonic spectra of five-level DCM-VSC in terms of negative sequence components.
(a) (b)


### 4.4. Case4: Unbalanced Parameters of a Three-Phase, Four-Wire DCM-VSC

In this case, we investigate a three-phase, four-wire DCM-VSC system.RN and LN are set equal to Rb and Lb, respectively. The inductance of phase A, La, is increased by 20% whereas other parameters remain fixed. Figure 16 shows the portion of [ABCD]012, which relates the positive sequence input voltage harmonics Vs1 to the zero sequence output current harmonics I0. Figure 17 shows the ac current harmonic spectra in terms of the positive, negative, and zero sequence components. The figure clearly shows the negative sequence second harmonic due to the unbalanced inductance values. In addition, it also clearly indicates that the derived ABCD matrix can accurately predict the zero sequence current harmonic at various frequency, and their values are in strong agreement with those obtained by PSCAD/EMTDC.Figure 16
The ABCD matrix, relatingVs1 to I0.Figure 17
Case 4: (a) the current harmonic spectra of in terms of positive sequence components; (b) the current harmonic spectra of in terms of negative sequence components; (c) the current harmonic spectra of in terms of zero sequence components.
(a) (b) (c)


### 4.5. Case5: A Multiconverter Unbalanced Microgrid System

Figure18 shows a microgrid system, consisting of multiple VSCs, which is similar to the Aichi microgrid system in Japan [30]. With the aid of the proposed ABCD matrix, the overall microgrid formulation can be derived.Figure 18
A multiconverter microgrid system.From Figure18, one can see the following.(29)Vpcc=Vs-ZsIs(30)Is=Ivsc1+Ivsc2+Ivsc3+⋯+IvscnAlso, each converter can be expressed as(31)IvsckVdck=AkBkCkDkVpccIdck,where k=1,2,3,⋯,n.Combining (29), (30), and (31), one can get (32).(32)IsVdc1Vdc2⋮Vdcn=M-1NVsIdc1Idc2⋮Idcnwhere(33)M=I+A1+A2+⋯+AnZs00⋯0C1ZsI0⋯0C2Zs0I⋯0⋮⋮⋮⋱0CnZs00⋯I,N=A1+A2+⋯+AnB1B2⋯BnC1D10⋯0C20D2⋯0⋮⋮⋮⋱0Cn00⋯DnFigure19 shows a multi-DCM-VSC microgrid system. Assume that the ac grid is unbalanced and takes the form of (27). VSC1 in Figure 19 is a three-level DCM-VSC, having unequal DC voltage due to DGs with different voltage level, as described in Case 2. VSC2 is a 5-level DCM-VSC with inherent unbalanced capacitor voltage as in Case 3. Finally, VSC3 is a three-level VSC with unbalanced parameters as delineated in Case 4. Employing (32), one can obtain the ac currents and capacitor voltages for each converter. Figures 20, 21, and 22 show the ac current spectra for VSC1, VSC2, and VSC3, respectively. To show that the presence of the neighboring converters can have a significant impact on the harmonic spectra of the converter, we resimulate Case 4 with the grid voltage taking the form of (27), and its current spectra are shown in Figure 23. Comparing Figure 22 with Figure 23, one can see that VSC3 exhibits a positive sequence, negative sequence, and zero sequence of the fundamental current of 45.5 A, 7.7 A, and 1.7 A, respectively, when it is the only converter connected to the grid. On the other hand, those components reduce to 37 A, 5.3 A, and 1.22 A, which correspond to 18%, 31%, and 26% reduction, respectively, under the microgrid environment. Moreover, as can be observed from Figures 22(b) and 23(b), one can notice that a negative sequence of second harmonic is injected from other converters under the microgrid environment whereas no even harmonics exist in Figure 23(b). All the results predicted by the proposed method are consistent with those of PSCAD/EMTDC.Figure 19
An unbalanced microgrid.Figure 20
Current spectra of VSC1: (a) positive sequence ofIvsc1; (b) negative sequence of Ivsc1.
(a) (b)Figure 21
Current spectra of VSC2: (a) positive sequence ofIvsc2; (b) negative sequence of Ivsc2.
(a) (b)Figure 22
Current spectra of VSC3: (a) positive sequence ofIvsc3; (b) negative sequence of Ivsc3; (c) zero sequence of Ivsc3.
(a) (b) (c)Figure 23
Only VSC3 is connected to the unbalanced ac grid: (a) positive sequence ofIvsc3; (b) negative sequence of Ivsc3; (c) zero sequence of Ivsc3.
(a) (b) (c)


## 4.1. Case1: Balanced and Unbalanced AC Grid (Negative Sequence Voltage Injection)

The three-level DCM-VSC is first connected to a balanced AC grid. Figure7 shows the maximum difference (Id) of the ac currents obtained by the proposed method and those by PSCAD/EMTDC under various time steps. As the figure shows, the simulation step tstep needs to reduce to 0.1μs in order for PSCAD/EMTDC model to reach the same level of accuracy as the proposed method. Figure 8 shows the resulting balanced three-phase currents for both PSCAD/EMTDC (tstep=0.1μs) and the proposed method. Nevertheless, the calculation time for PSCAD/EMTDC with time step of 0.1μs is about 85 times longer than the proposed method. Figure 9 shows the current harmonic spectra in terms of its positive and negative sequence components. It clearly shows that under balanced condition, the DCM-VSC only produces characteristic harmonics.Figure 7
I d versus solution time step.Figure 8
Balanced steady-state AC current waveforms.Figure 9
(a) The harmonic spectra of ac currents of Figure8 in terms of positive sequence components; (b) the harmonic spectra of ac currents of Figure 8 in terms of negative sequence components.
(a) (b)It is to be noted that input stimuli in (2) can be easily extended to other stimuli such as the cosine function by means of Euler’s formula, which is stated in (20).(20)cos⁡ωt+θ=ejωt+θ+e-jωt+θ2The following, as an example, are the step-by-step procedures for calculating the positive sequence of the fundamental current and dc voltage harmonic under the balanced condition:(1)
The input stimuli for phases A, B, and C are{(vs/2)e-jωt+θ,(vs/2)ej(ωt+θ)}, {(vs/2)e-j(ωt+θ-2π/3),(vs/2)ej(ωt+θ-2π/3)}, and {(vs/2)e-j(ωt+θ+2π/3),(vs/2)ej(ωt+θ+2π/3)}. Note that vs is the magnitude of the source voltage.(2)
Substitute the parameter values listed in Table3 into the form, similar to (5)(3)
Evaluate (7) to obtain Φˇ.(4)
Extract the submatrices,A, B, C, and D, which are as follows.(21)A=1.6490+0.0489j1.6709-0.2624j-0.7822-1.4525j-1.0626+1.5416j-0.8668+1.4036j-0.6083-1.2792j-0.7822-1.4525j-0.6083-1.2792j-0.8668+1.4036j1.6709+0.2624j1.6490+0.0489j-1.626+1.5416j-0.8668+1.4036j-1.0626+1.5416j1.6490+0.0489j-0.6083-1.2792j-0.7822-1.4525j1.6709-0.2624j(22)B=0.1887+0.0148j-0.1072+0.1560j-0.0815-0.1708j(23)C=5.8125-0.3137j5.8125+0.3137j-3.178-4.8769j-3.178+4.8769j-2.6345+5.1907j-2.6345-5.1907j5.8125-0.3137j5.8125+0.3137j-3.178-4.8769j-3.178+4.8769j-2.6345+5.1907j-2.6345-5.1907j(24)D=00(5)
The ac current and dc voltage harmonics can be obtained by evaluating (14). Note that the ac current harmonics need to be multiplied by two to obtain the correct phasor forms(25)Ia1Ib1Ic1Vc10Vc20=200000200000200000100000122.76-0.34j11.09+19.88j11.67-19.54j292.03292.03=-45.53-0.68j22.18+39.77j23.35-39.09j292.03292.03,whereIa1, Ib1, and Ic1 denote the fundamental currents of phases A, B, and C, respectively, and Vc10 and Vc20 are the two dc voltages across C1 and C2(6)
Evaluate (19), and the fundamental positive (I11), negative (I21), and zero sequence (I01) currents and the two dc voltages are stated as (26).(26)I01I11I21Vc10Vc20=045.53-0.68j0292.03292.03Next, a5% negative sequence of the fundamental voltage is added to the AC grid, and the voltage takes the following form.(27)Vsa=Vsejωt+0.05VsejωtVsb=Vsejωt-120∘+0.05Vsejωt+120∘Vsc=Vsejωt+120∘+0.05Vsejωt-120∘Figure10 shows the harmonic spectra of this unbalanced case in terms of positive and negative sequence components. In this instance, the uncharacteristic harmonics such as the negative sequence fundamental and third harmonics are accurately predicted by the proposed method.Figure 10
Unbalanced ac grid case: (a) the harmonic spectra of ac currents of in terms of positive sequence components; (b) the harmonic spectra of ac currents of in terms of negative sequence components.
(a) (b)


## 4.2. Case2: Unequal DC Voltage due to DGs with Different Voltage Levels

In this case, we try to simulate the case where the split dc link capacitors of the DCM-VSC can be controlled independently. A commonly seen situation is to have a different PV panel connected to each capacitor, with each panel being controlled independently to achieve its own point of maximum power [14]. Such a case leads to unequal dc voltage level for different capacitors with respect to the neutral point, N. To simulate such a case, two DC sources with vDC1=350 V connected in parallel with C1 and vDC1=230 V connected in parallel with C2 are connected to the dc link. The resulting current spectra, shown in Figure 11, feature undesired negative sequence even-order current harmonics such as the second and 14-th harmonics. The values of these harmonics are accurately predicted by the proposed method when compared with those of PSCAD/EMTDC.Figure 11
Case 2: (a) the harmonic spectra of ac currents in terms of positive sequence components; (b) the harmonic spectra of ac currents in terms of negative sequence components.
(a) (b)


## 4.3. Case3: Inherent Unbalanced Capacitor Voltage for A 5-Level DCM-VSC

Different from the previous case, this case study investigates a five-level DCM-VSC. The parameter settings are the same as those in Table3. As stated in [16, 29], when the level number is greater than 3, an inherent dc voltage imbalance becomes apparent. In addition, as demonstrated in [22], if the condition(28)macos⁡δ>0.55is satisfied, no balancing strategy exists to balance the capacitor voltages. Hence, for the chosen ma and δ in Table 3, voltage imbalance among capacitors is expected.Figures12(a) and 12(b) show the DC voltage waveforms for each capacitor (in this paper, the capacitor labels follow those indicated in Figure 6), obtained by the proposed method and PSCAD/EMTDC. Their corresponding harmonic spectra are in great agreement as shown in Figure 13. Figure 14 shows the current waveforms of the five-level DCM-VSC. Figure 15 shows the corresponding current spectra. Compared with Figure 9, one can see that the positive sequence 7-th, 13-th, 19-th, and 25-th harmonics are much greater for the 5-level DCM-VSC. Moreover, for the negative sequence, the value of the 5-th harmonic is about 5 times larger than that of 11-th harmonic for the 5-level case, whereas the 11-th harmonic for the 3-level converter is 3 times higher than that of the 5-th harmonic as seen in Figure 9. Nevertheless, different from Case 2, there is no even current harmonic generated at the ac side. The reason is that, with respect to N, there exists a symmetry between the upper and lower capacitor voltages. In other words, vc1=vc4 and vc2=vc3. Consequently, such a symmetry with respect to N will not generally lead to any even current harmonics at the ac side.Figure 12
Case 3: (a) the capacitor voltage waveforms obtained by the proposed method; (b) the capacitor voltage waveforms obtained by the PSCAD/EMTDC.
(a) (b)Figure 13
Harmonic spectra of each capacitor in the five-level DCM-VSC.Figure 14
Steady-state AC current waveforms of five-level DCM-VSC.Figure 15
Case 3: (a) the current harmonic spectra of five-level DCM-VSC in terms of positive sequence components; (b) the current harmonic spectra of five-level DCM-VSC in terms of negative sequence components.
(a) (b)


## 4.4. Case4: Unbalanced Parameters of a Three-Phase, Four-Wire DCM-VSC

In this case, we investigate a three-phase, four-wire DCM-VSC system.RN and LN are set equal to Rb and Lb, respectively. The inductance of phase A, La, is increased by 20% whereas other parameters remain fixed. Figure 16 shows the portion of [ABCD]012, which relates the positive sequence input voltage harmonics Vs1 to the zero sequence output current harmonics I0. Figure 17 shows the ac current harmonic spectra in terms of the positive, negative, and zero sequence components. The figure clearly shows the negative sequence second harmonic due to the unbalanced inductance values. In addition, it also clearly indicates that the derived ABCD matrix can accurately predict the zero sequence current harmonic at various frequency, and their values are in strong agreement with those obtained by PSCAD/EMTDC.Figure 16
The ABCD matrix, relatingVs1 to I0.Figure 17
Case 4: (a) the current harmonic spectra of in terms of positive sequence components; (b) the current harmonic spectra of in terms of negative sequence components; (c) the current harmonic spectra of in terms of zero sequence components.
(a) (b) (c)


## 4.5. Case5: A Multiconverter Unbalanced Microgrid System

Figure18 shows a microgrid system, consisting of multiple VSCs, which is similar to the Aichi microgrid system in Japan [30]. With the aid of the proposed ABCD matrix, the overall microgrid formulation can be derived.Figure 18
A multiconverter microgrid system.From Figure18, one can see the following.(29)Vpcc=Vs-ZsIs(30)Is=Ivsc1+Ivsc2+Ivsc3+⋯+IvscnAlso, each converter can be expressed as(31)IvsckVdck=AkBkCkDkVpccIdck,where k=1,2,3,⋯,n.Combining (29), (30), and (31), one can get (32).(32)IsVdc1Vdc2⋮Vdcn=M-1NVsIdc1Idc2⋮Idcnwhere(33)M=I+A1+A2+⋯+AnZs00⋯0C1ZsI0⋯0C2Zs0I⋯0⋮⋮⋮⋱0CnZs00⋯I,N=A1+A2+⋯+AnB1B2⋯BnC1D10⋯0C20D2⋯0⋮⋮⋮⋱0Cn00⋯DnFigure19 shows a multi-DCM-VSC microgrid system. Assume that the ac grid is unbalanced and takes the form of (27). VSC1 in Figure 19 is a three-level DCM-VSC, having unequal DC voltage due to DGs with different voltage level, as described in Case 2. VSC2 is a 5-level DCM-VSC with inherent unbalanced capacitor voltage as in Case 3. Finally, VSC3 is a three-level VSC with unbalanced parameters as delineated in Case 4. Employing (32), one can obtain the ac currents and capacitor voltages for each converter. Figures 20, 21, and 22 show the ac current spectra for VSC1, VSC2, and VSC3, respectively. To show that the presence of the neighboring converters can have a significant impact on the harmonic spectra of the converter, we resimulate Case 4 with the grid voltage taking the form of (27), and its current spectra are shown in Figure 23. Comparing Figure 22 with Figure 23, one can see that VSC3 exhibits a positive sequence, negative sequence, and zero sequence of the fundamental current of 45.5 A, 7.7 A, and 1.7 A, respectively, when it is the only converter connected to the grid. On the other hand, those components reduce to 37 A, 5.3 A, and 1.22 A, which correspond to 18%, 31%, and 26% reduction, respectively, under the microgrid environment. Moreover, as can be observed from Figures 22(b) and 23(b), one can notice that a negative sequence of second harmonic is injected from other converters under the microgrid environment whereas no even harmonics exist in Figure 23(b). All the results predicted by the proposed method are consistent with those of PSCAD/EMTDC.Figure 19
An unbalanced microgrid.Figure 20
Current spectra of VSC1: (a) positive sequence ofIvsc1; (b) negative sequence of Ivsc1.
(a) (b)Figure 21
Current spectra of VSC2: (a) positive sequence ofIvsc2; (b) negative sequence of Ivsc2.
(a) (b)Figure 22
Current spectra of VSC3: (a) positive sequence ofIvsc3; (b) negative sequence of Ivsc3; (c) zero sequence of Ivsc3.
(a) (b) (c)Figure 23
Only VSC3 is connected to the unbalanced ac grid: (a) positive sequence ofIvsc3; (b) negative sequence of Ivsc3; (c) zero sequence of Ivsc3.
(a) (b) (c)


## 5. Conclusion

An AC-DC converter has been more conveniently characterized by the ABCD matrix for harmonic modeling and analysis. This paper presented a time-domain-based ABCD matrix for a diode-clamped multilevel voltage source converter. The derived model produced accurate predictions of characteristic and all the uncharacteristic harmonics for various unbalanced conditions. Also, the method yielded accurate results and does not suffer from harmonic truncation errors. Moreover, all the predicted results are in great agreement with those obtained by PSCAD/EMTDC with the calculation time being much shorter.



---

*Source: 3250761-2018-12-05.xml*