# Analytical Approach to Circulating Current Mitigation in Hexagram Converter-Based Grid-Connected Photovoltaic Systems Using Multiwinding Coupled Inductors

Authors: Abdullrahman A. Al-Shamma’a; Abdullah M. Noman; Khaled E. Addoweesh; Ayman A. Alabduljabbar; A. I. Alolah
Journal: International Journal of Photoenergy
 (2018)

---

## Abstract

The hexagram multilevel converter (HMC) is composed of six conventional two-level voltage source converters (VSCs), where each VSC module is connected to a string of PV arrays. The VSC modules are connected through inductors, which are essential to minimize the circulating current. Selecting inductors with suitable inductance is no simple process, where the inductance value should be large to minimize the circulating current as well as small to reduce an extra voltage drop. This paper analyzes the utilization of a multiwinding (e.g., two, three, and six windings) coupled inductor to interconnect the six VSC modules instead of six single inductors, to minimize the circulating current inside the HMC. Then, a theoretical relationship between the total impedance to the circulating current, the number of coupled inductor windings, and the magnetizing inductance is derived. Owing to the coupled inductors, the impedance on the circulating current path is a multiple of six times the magnetizing inductance, whereas the terminal voltage is slightly affected by the leakage inductance. The HMC is controlled to work under variable solar radiation, providing active power to the grid. Additional functions such as DSTATCOM, during daytime, are also demonstrated. The controller performance is found to be satisfactory for both active and reactive power supplies.


---

## Body

## 1. Introduction

Recently, photovoltaic (PV) energy systems have gained more attention as distributed generation units, as they offer low cost of generation closer to that of conventional plants, as well as less maintenance and no grid noise [1, 2]. Moreover, PV systems can solve multiple typical problems present in conventional AC power systems. However, PV systems present frequency and voltage fluctuations when islanding operation occurs. Therefore, PV plants should be integrated with the power system in order to maintain the overall frequency and voltage at a stable condition. Several studies suggest interconnection methods of PV systems to the grid through voltage source converters (VSCs), because they provide versatile functions enhancing capabilities of the system [3–5]. The main purpose of the VSC is to connect the PV plant to the grid while guaranteeing power quality (PQ) standards. However, the high-frequency switching of VSCs introduces additional harmonic components to the system, hence creating PQ problems if not implemented accurately [6].Most of the available VSCs for PV systems are traditional two-level three-phase VSCs with low power capacity [7]. In the literature, researchers have suggested numerous multilevel topologies for grid-connected PV plants [8–10]. Multilevel VSCs are considered more attractive than traditional two-level VSCs, as they improve the output voltage quality and reduce electromagnetic interference, voltage stress on IGBTs, and common-mode voltage. Moreover, multilevel VSCs operate at a low switching frequency and hence increase system efficiency [10]. Consequently, multilevel VSCs have been widely used in chemical, oil, and different kinds of plants, as well as in power plants, transmission systems, and PQ compensators [11].Current research is mainly focused on three specific multilevel VSC topologies, namely, the neutral-point clamped (NPC) [12], flying capacitor (FC) [13] and the Cascaded H-bridge (CHB) [14] topologies. The NPC VSC requires a high number of clamping diodes to increase the number of voltage levels, which can cause problems related to switches with different ratings, high-voltage rating in blocking diodes, and capacitor voltage imbalance. However, the FC VSC consists of a large number of capacitors, which leads to complications in regulating capacitor voltages. Finally, the CHB VSC provides isolated DC sources, hence being appropriate for use in PV systems [15]. In addition, other advantages of this topology include its modular structure and the reduced number of components compared to the other multilevel converters (i.e., NPC and FC). Therefore, the CHB VSC can reach the same number of voltage levels with a simpler assembly and maintenance. However, this topology presents three main drawbacks: (i) high overall component count; (ii) high energy storage requirement because the instantaneous power related to each H-bridge varies at twice the fundamental frequency, given its single-phase modular structure; and (iii) difficulty to control the voltage across DC-link capacitors.The concept of interconnecting three traditional three-phase VSCs to produce a multilevel converter was first proposed in [16] and applied to medium-voltage variable-speed drives. In [17], three three-phase two-level VSCs are interconnected using three single-phase transformers with a 1 : 1 turn ratio. These interconnected transformers increase the output voltage and suppress the circulating current inside the converter. The power capacity of the overall converter is three times the capacity of each interconnected converter, whereas the volt-ampere rating of each intermediate transformer is equal to that of each interconnected converter. In [18], another topology known as the hexagram multilevel converter (HMC) is proposed, which combines six two-level converters by using six inductors. This topology shares many advantages of the CHB but uses fewer switches and reduces the size of the DC-link capacitor [19]. The advantages of the proposed topology can be summarized as follows: (1) only six standard three-phase VSCs are necessary to generate multilevel output voltage; (2) each VSC module is balanced in operation, equally loaded, and supplies 1/6 output power; (3) modular construction which facilitates system maintenance and spare management; (4) only six isolated DC links with no voltage unbalance problem; (5) low number of power electronics switches and low DC energy storage requirement; and (6) the output transformer contributes to higher output voltage. Figures 1 and 2 are schematic diagrams illustrating multilevel topologies based on cascaded two-level VSCs according to previous research.Figure 1
Topologies of cascade two-level converters.Figure 2
HMC electrical diagram.However, selecting an inductor with suitable inductance is no simple process, where the inductance value should be adequate to minimize the circulating current as well as to reduce an extra voltage drop on the inductors that affects the terminal voltage. This difficulty can be circumvented by using a multiwinding coupled inductor. Therefore, this paper analyzes the utilization of a multiwinding (e.g., two, three, and six windings) coupled inductor to interconnect the six VSCs instead of six single inductors. Consequently, both goals minimize the circulating current and the minimal effects to the output voltage can be accomplished instantaneously. Then, an analytical model to calculate the total inductance imposed to the circulating current path is derived. The equivalent circuit model of the HMC is detailed in the abc reference frame and then transformed into the orthogonal dq0 reference frame. Moreover, to extract the maximum power from the PV arrays, a control algorithm for maximum power point tracking is also presented.The paper is organized as follows: Section2 describes the HMC. Section 3 presents the mathematical model for a grid-connected PV system. Section 4 demonstrates the benefits of the multiwinding coupled inductor inside the HMC. Section 5 describes the control system and modulation strategy. Section 6 shows the simulation results followed by the corresponding discussion. Finally, the conclusions achieved from the present work are summarized in Section 7.


## 2. System Description

The proposed configuration of the three-phase grid-connected PV plants is shown in Figure3. This configuration is composed of six traditional three-phase two-level VSC modules, which have a hexagonal interconnection to produce higher voltage levels as shown in Figure 2. The circulating current in the obtained loop can be suppressed using six inductors, which present a small impedance at the switching frequency. Each VSC module supplies 1/6 of the converter output power. In addition, each module has one of its AC terminals (i.e., a, b, or c) designated as converter output. The three-phase AC output terminals of the HMC are labeled as follows: A (AC terminal a1 of module 1), B (AC terminal b3 of module 3), and C (AC terminal c5 of module 5). The remaining three-phase AC output terminals are labeled as follows: A′ (AC terminal a4 of module 4), B′ (AC terminal b6 of module 6), and C′ (AC terminal c2 of module 2). The other two AC terminals of each module are, respectively, connected to an adjacent module through an inductor. For instance, AC terminal b1 of module 1 is coupled to AC terminal b2 of module 2 through an inductor Ld; AC terminal c1 of module 1 is coupled to AC terminal c6 of module 6 through another inductor. The remaining modules consist of similar connections, as shown in Figure 2.Figure 3
Complete HMC for a grid-connected PV system.As shown in Figure3, the PV system based on the HMC has six output terminals. Therefore, when linked to the three-phase electrical grid, an open-end winding (OEW) transformer is necessary to provide these six terminals. The secondary windings of the OEW transformer are connected differentially. In the proposed converter, AC terminals A-A′, B-B′, and C-C′ are used to provide phases A, B, and C, respectively. High-voltage windings are arranged in a star configuration and coupled directly to the three-phase grid. Each VSC module is supplied with a separate PV string to produce two-level individual outputs. The proposed converter topology presents several advantages, such as the use of fewer power electronic switches, diodes, and storage capacitors when compared with other topologies, making it suitable for renewable energy and other applications.


## 3. System Modeling

### 3.1. PV System Modeling and MPPT Control

Figure4 shows the equivalent model of a PV array. The PV array is composed of several series-parallel connected solar cells. The basic equation of the PV array is given by [20].
(1)I=Ipv−I0expV+RsIaVt−1−V+RsIRp,where Ipv=NpIpv,cell is the PV current of the array, Ipv,cell is the current generated by incident light (directly proportional to Sun irradiation), Np is the number of cells connected in parallel, I0=NpI0,cell is the saturation current of the array, I0,cell is the reverse saturation or leakage current of the diode, Vt=NskT/q is the thermal voltage of the array, Ns is the number of cells connected in series, q is the electron charge, k is the Boltzmann constant, T is the temperature of the p-n junction, a is the diode ideality constant, Rs is the equivalent series resistance of the array, and Rp is the equivalent parallel resistance. Table 1 shows the nominal parameters of the KC200GT PV array. Figure 5 shows the I-V and P-V curves obtained using (1) of the PV array under variable solar radiation.Figure 4
Equivalent model of a PV array.Table 1
Parameters of the KC200GT PV array at 25°C, AM1.5, and 1000 W/m2.
Imp(A)
Vmp(V)
Pmax(W)
Isc(A)
Voc(V)
Kv(V/K)
KI(A/K)
Nss(No.)
Rp(Ω)
Rs(Ω)
Io,n(A)
a
7.61
26.3
200.143
8.21
32.9
−0.123
0.0032
54
415.405
0.221
9.825 × 10−8
1.3Figure 5
I-V and P-V curves of the KC200GT PV array.


### 3.2. Maximum Power Point Tracking

The proposed system consists of 12 PV modules distributed as six pairs of series-connected modules and coupled directly to a single DC link (see Figure3). Therefore, the power rating of the proposed converter is 2.4 kW. The controller is aimed to regulate the voltage of the PV-module pair at 52.6 V, to guarantee maximum power transfer to the grid.Next, the aim of the MPPT control algorithm used in this paper is to ensure that under any solar radiation and temperature conditions, the maximum power is extracted from the PV modules. This is achieved by matching the PV-array maximum power point to the corresponding operating voltage and current of the converters. The perturb and observe (P&O) algorithm is a commonly used MPPT technique because it is easy to implement [21]. The operating principle of the P&O algorithm is shown in Figure 6. This algorithm measures the PV plant voltage and current, then it varies the operating voltage and compares the power received between the two voltage values. After each perturbation, the algorithm compares the output power from the PV before and after the perturbation. The direction of a new perturbation depends upon the output power: the perturbation will follow the same direction if higher power is measured when the output voltage varies and the opposite direction otherwise. These procedures are repeated continuously, and the reference voltage is generated and fed to the converter controller.Figure 6
Flowchart for P&O MPPT method.The numerical illustration of the P&O algorithm [2] is given as follows:
(2)dPpvkdVpvk=Ppvk−Ppvk−1Vpvk−Vpvk−1.Here,Ppvk and Ppvk−1 stand for current power and previous measured power, while Vpvk and Vpvk−1 stand for current PV voltage and previous one.


### 3.3. Voltage and Current Analysis

Under the symmetrical operation condition, the fundamental components of the six VSC modules are the same. The corresponding open-circuit voltage of the HMC is presented in Figure7.Figure 7
Phasor voltage of the HMC.


#### 3.3.1. Current Relations

The phase currents of each VSC modules fulfill the following expressions:(3)ia1+ib1+ic1ia2+ib2+ic2ia3+ib3+ic3ia4+ib4+ic4ia5+ib5+ic5ia6+ib6+ic6=0.Assuming that the intermediate coupled inductors have large magnetizing inductance, the circulating currents are suppressed to a low value and can be ignored, hence(4)ib1+ia2+ic3+ib4+ia5+ic6=0.Because every pair of the six VSC modules is coupled, the current within the proposed converter has the following expressions:(5)ib1ia2ic3ib4ia5ic6=−ib2ia3ic4ib5ia6ic1.Next, suppose that the HMC is linked to a three-phase grid. Then, the output currents will satisfy the following equation:(6)IAIBIC=ia1ib3ic5=−ia4ib6ic2=I0I−120I+120.Using (2), (3), and (4), the output current of each VSC module can be shown to be
(7)ia1ib1ic1=ia3ib3ic3=ia5ib5ic5=−ia2ib2ic2=−ia4ib4ic4=−ia6ib6ic6=IAIBIC=I0I−120I+120,where IA,IB,ICT and I are the converter output phase currents and their RMS values, respectively. It can be concluded from (5) that each VSC module within the HMC will have an identical current under the symmetrical operation condition.


#### 3.3.2. Voltage Relations

Under the symmetrical operation condition, the fundamental component (RMS) of the phase voltages at the VSC modules can be described as(8)Va1o1Vb1o1Vc1o1=Va3o3Vb3o3Vc3o3=Va5o5Vb5o5Vc5o5=−Va2o2Vb2o2Vc2o2=−Va4o4Vb4o4Vc4o4=−Va6o6Vb6o6Vc6o6=V0V−120V+120,(9)V=122Vdcma,where V, Vdc, and ma are the RMS values of each VSC-module phase voltage, the DC-link voltage of each VSC module, and the amplitude modulation index, respectively; subscripts o1, o2, o3, o4, o5, and o6 represent the virtual neutral points of each VSC.The output voltages of the proposed converter are written as(10)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=Va1o1−Va4o4+Vc4o4−Vc3o3+Va3o3−Va2o2+Vb2o2−Vb1o1Vb3o3−Vb6o6+Va6o6−Va5o5+Vb5o5−Vb4o4+Vc4o4−Vc3o3Vc5o5−Vc2o2+Vb2o2−Vb1o1+Vc1o1−Vc6o6+Va6o6−Va5o5−2Zdia1ib3ic5,where impedance Zd=jωLd.Thus, using (5), (6), (7), and (8), the net output voltages of the HMC under the symmetrical operation condition are given by
(11)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=6V06V−1206V+120−2ZdIAIBIC.Consequently, (9) demonstrates that the three-phase voltage of the HMC is approximately six times the phase voltage of each VSC module. In other words, the voltage stress is reduced by a factor of six.


### 3.4. Equivalent Circuit Model

Based on Figure3, the line voltages of the HMC are expressed by
(12)UABUBCUCA=VAA′−VBB′VBB′−VCC′VCC′−VAA′=63V−3063V−15063V+90−2jωLdIA−IBIB−ICIC−IA=EABEBCECA+jωLfIAIBIC−jωLfIBICIA,where EAB,EBC,ECAT are the RMS values of the line-to-line voltages in the grid, UAB,UBC,UCAT and IA,IB,ICT are the RMS values of line voltages and the RMS values of the output current, respectively, and Lf is the AC-side filtering inductance.When the conversion powers between the VSC modules within the HMC are balanced, the six DC-link voltages are equal toVDCav, hence
(13)VDCav=Vdc1=Vdc2=Vdc3=Vdc4=Vdc5=Vdc6=Vdc1+Vdc2+Vdc3+Vdc4+Vdc5+Vdc66.Moreover, the equivalent DC-link voltage can be expressed as(14)UDC_eq=6VDC_av.If each VSC is driven using the same PWM, the relationship between the DC-link voltages and the phase voltages can be expressed as(15)V=122UDC_eqma.From (12) and (14), the HMC can be modeled by
(16)UAB¯UBC¯UCA¯=EABEBCECA+jωLf+2LdIAIBIC−jωLf+2LdIBICIA,where
(17)UAB¯UBC¯UCA¯=63V−3063V−15063V+90.Using (14) and (17), the HMC can be modeled as a conventional three-phase two-level VSC, as shown in Figure 8.Figure 8
Equivalent circuit model of HMC.Figure8 shows the equivalent circuit model of the HMC. The equivalent circuit is composed of an open-circuit voltage source linked in series to the HMC output impedance. The open-circuit voltage of the HMC is shown in Figure 7. As shown in Figure 8, the equivalent interface inductor of the HMC is 2Ld+Lf. This feature is advantageous, because the filtering inductor Lf could be minimalized or even removed. However, this can be inopportune when the line transformer has sufficient inductance for filtering the output current. In this case, essentially if the HMC supplies reactive power, the output impedance causes a voltage drop, decreasing the power capability. Thus, it might be better to minimize the interconnected total inductances. However, the inductances are essential to suppress the circulating current within the HMC. The aforementioned decision shows that the interconnected inductance value selection is not a simple process, where the inductance should be large enough to minimize the circulating current but small enough to avoid an extra voltage drop. In the following section, this feature is investigated using multiwinding coupled inductors.


## 3.1. PV System Modeling and MPPT Control

Figure4 shows the equivalent model of a PV array. The PV array is composed of several series-parallel connected solar cells. The basic equation of the PV array is given by [20].
(1)I=Ipv−I0expV+RsIaVt−1−V+RsIRp,where Ipv=NpIpv,cell is the PV current of the array, Ipv,cell is the current generated by incident light (directly proportional to Sun irradiation), Np is the number of cells connected in parallel, I0=NpI0,cell is the saturation current of the array, I0,cell is the reverse saturation or leakage current of the diode, Vt=NskT/q is the thermal voltage of the array, Ns is the number of cells connected in series, q is the electron charge, k is the Boltzmann constant, T is the temperature of the p-n junction, a is the diode ideality constant, Rs is the equivalent series resistance of the array, and Rp is the equivalent parallel resistance. Table 1 shows the nominal parameters of the KC200GT PV array. Figure 5 shows the I-V and P-V curves obtained using (1) of the PV array under variable solar radiation.Figure 4
Equivalent model of a PV array.Table 1
Parameters of the KC200GT PV array at 25°C, AM1.5, and 1000 W/m2.
Imp(A)
Vmp(V)
Pmax(W)
Isc(A)
Voc(V)
Kv(V/K)
KI(A/K)
Nss(No.)
Rp(Ω)
Rs(Ω)
Io,n(A)
a
7.61
26.3
200.143
8.21
32.9
−0.123
0.0032
54
415.405
0.221
9.825 × 10−8
1.3Figure 5
I-V and P-V curves of the KC200GT PV array.


## 3.2. Maximum Power Point Tracking

The proposed system consists of 12 PV modules distributed as six pairs of series-connected modules and coupled directly to a single DC link (see Figure3). Therefore, the power rating of the proposed converter is 2.4 kW. The controller is aimed to regulate the voltage of the PV-module pair at 52.6 V, to guarantee maximum power transfer to the grid.Next, the aim of the MPPT control algorithm used in this paper is to ensure that under any solar radiation and temperature conditions, the maximum power is extracted from the PV modules. This is achieved by matching the PV-array maximum power point to the corresponding operating voltage and current of the converters. The perturb and observe (P&O) algorithm is a commonly used MPPT technique because it is easy to implement [21]. The operating principle of the P&O algorithm is shown in Figure 6. This algorithm measures the PV plant voltage and current, then it varies the operating voltage and compares the power received between the two voltage values. After each perturbation, the algorithm compares the output power from the PV before and after the perturbation. The direction of a new perturbation depends upon the output power: the perturbation will follow the same direction if higher power is measured when the output voltage varies and the opposite direction otherwise. These procedures are repeated continuously, and the reference voltage is generated and fed to the converter controller.Figure 6
Flowchart for P&O MPPT method.The numerical illustration of the P&O algorithm [2] is given as follows:
(2)dPpvkdVpvk=Ppvk−Ppvk−1Vpvk−Vpvk−1.Here,Ppvk and Ppvk−1 stand for current power and previous measured power, while Vpvk and Vpvk−1 stand for current PV voltage and previous one.


## 3.3. Voltage and Current Analysis

Under the symmetrical operation condition, the fundamental components of the six VSC modules are the same. The corresponding open-circuit voltage of the HMC is presented in Figure7.Figure 7
Phasor voltage of the HMC.


### 3.3.1. Current Relations

The phase currents of each VSC modules fulfill the following expressions:(3)ia1+ib1+ic1ia2+ib2+ic2ia3+ib3+ic3ia4+ib4+ic4ia5+ib5+ic5ia6+ib6+ic6=0.Assuming that the intermediate coupled inductors have large magnetizing inductance, the circulating currents are suppressed to a low value and can be ignored, hence(4)ib1+ia2+ic3+ib4+ia5+ic6=0.Because every pair of the six VSC modules is coupled, the current within the proposed converter has the following expressions:(5)ib1ia2ic3ib4ia5ic6=−ib2ia3ic4ib5ia6ic1.Next, suppose that the HMC is linked to a three-phase grid. Then, the output currents will satisfy the following equation:(6)IAIBIC=ia1ib3ic5=−ia4ib6ic2=I0I−120I+120.Using (2), (3), and (4), the output current of each VSC module can be shown to be
(7)ia1ib1ic1=ia3ib3ic3=ia5ib5ic5=−ia2ib2ic2=−ia4ib4ic4=−ia6ib6ic6=IAIBIC=I0I−120I+120,where IA,IB,ICT and I are the converter output phase currents and their RMS values, respectively. It can be concluded from (5) that each VSC module within the HMC will have an identical current under the symmetrical operation condition.


### 3.3.2. Voltage Relations

Under the symmetrical operation condition, the fundamental component (RMS) of the phase voltages at the VSC modules can be described as(8)Va1o1Vb1o1Vc1o1=Va3o3Vb3o3Vc3o3=Va5o5Vb5o5Vc5o5=−Va2o2Vb2o2Vc2o2=−Va4o4Vb4o4Vc4o4=−Va6o6Vb6o6Vc6o6=V0V−120V+120,(9)V=122Vdcma,where V, Vdc, and ma are the RMS values of each VSC-module phase voltage, the DC-link voltage of each VSC module, and the amplitude modulation index, respectively; subscripts o1, o2, o3, o4, o5, and o6 represent the virtual neutral points of each VSC.The output voltages of the proposed converter are written as(10)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=Va1o1−Va4o4+Vc4o4−Vc3o3+Va3o3−Va2o2+Vb2o2−Vb1o1Vb3o3−Vb6o6+Va6o6−Va5o5+Vb5o5−Vb4o4+Vc4o4−Vc3o3Vc5o5−Vc2o2+Vb2o2−Vb1o1+Vc1o1−Vc6o6+Va6o6−Va5o5−2Zdia1ib3ic5,where impedance Zd=jωLd.Thus, using (5), (6), (7), and (8), the net output voltages of the HMC under the symmetrical operation condition are given by
(11)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=6V06V−1206V+120−2ZdIAIBIC.Consequently, (9) demonstrates that the three-phase voltage of the HMC is approximately six times the phase voltage of each VSC module. In other words, the voltage stress is reduced by a factor of six.


## 3.3.1. Current Relations

The phase currents of each VSC modules fulfill the following expressions:(3)ia1+ib1+ic1ia2+ib2+ic2ia3+ib3+ic3ia4+ib4+ic4ia5+ib5+ic5ia6+ib6+ic6=0.Assuming that the intermediate coupled inductors have large magnetizing inductance, the circulating currents are suppressed to a low value and can be ignored, hence(4)ib1+ia2+ic3+ib4+ia5+ic6=0.Because every pair of the six VSC modules is coupled, the current within the proposed converter has the following expressions:(5)ib1ia2ic3ib4ia5ic6=−ib2ia3ic4ib5ia6ic1.Next, suppose that the HMC is linked to a three-phase grid. Then, the output currents will satisfy the following equation:(6)IAIBIC=ia1ib3ic5=−ia4ib6ic2=I0I−120I+120.Using (2), (3), and (4), the output current of each VSC module can be shown to be
(7)ia1ib1ic1=ia3ib3ic3=ia5ib5ic5=−ia2ib2ic2=−ia4ib4ic4=−ia6ib6ic6=IAIBIC=I0I−120I+120,where IA,IB,ICT and I are the converter output phase currents and their RMS values, respectively. It can be concluded from (5) that each VSC module within the HMC will have an identical current under the symmetrical operation condition.


## 3.3.2. Voltage Relations

Under the symmetrical operation condition, the fundamental component (RMS) of the phase voltages at the VSC modules can be described as(8)Va1o1Vb1o1Vc1o1=Va3o3Vb3o3Vc3o3=Va5o5Vb5o5Vc5o5=−Va2o2Vb2o2Vc2o2=−Va4o4Vb4o4Vc4o4=−Va6o6Vb6o6Vc6o6=V0V−120V+120,(9)V=122Vdcma,where V, Vdc, and ma are the RMS values of each VSC-module phase voltage, the DC-link voltage of each VSC module, and the amplitude modulation index, respectively; subscripts o1, o2, o3, o4, o5, and o6 represent the virtual neutral points of each VSC.The output voltages of the proposed converter are written as(10)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=Va1o1−Va4o4+Vc4o4−Vc3o3+Va3o3−Va2o2+Vb2o2−Vb1o1Vb3o3−Vb6o6+Va6o6−Va5o5+Vb5o5−Vb4o4+Vc4o4−Vc3o3Vc5o5−Vc2o2+Vb2o2−Vb1o1+Vc1o1−Vc6o6+Va6o6−Va5o5−2Zdia1ib3ic5,where impedance Zd=jωLd.Thus, using (5), (6), (7), and (8), the net output voltages of the HMC under the symmetrical operation condition are given by
(11)VAA′VBB′VCC′=Va1a4Vb3b6Vc2c5=6V06V−1206V+120−2ZdIAIBIC.Consequently, (9) demonstrates that the three-phase voltage of the HMC is approximately six times the phase voltage of each VSC module. In other words, the voltage stress is reduced by a factor of six.


## 3.4. Equivalent Circuit Model

Based on Figure3, the line voltages of the HMC are expressed by
(12)UABUBCUCA=VAA′−VBB′VBB′−VCC′VCC′−VAA′=63V−3063V−15063V+90−2jωLdIA−IBIB−ICIC−IA=EABEBCECA+jωLfIAIBIC−jωLfIBICIA,where EAB,EBC,ECAT are the RMS values of the line-to-line voltages in the grid, UAB,UBC,UCAT and IA,IB,ICT are the RMS values of line voltages and the RMS values of the output current, respectively, and Lf is the AC-side filtering inductance.When the conversion powers between the VSC modules within the HMC are balanced, the six DC-link voltages are equal toVDCav, hence
(13)VDCav=Vdc1=Vdc2=Vdc3=Vdc4=Vdc5=Vdc6=Vdc1+Vdc2+Vdc3+Vdc4+Vdc5+Vdc66.Moreover, the equivalent DC-link voltage can be expressed as(14)UDC_eq=6VDC_av.If each VSC is driven using the same PWM, the relationship between the DC-link voltages and the phase voltages can be expressed as(15)V=122UDC_eqma.From (12) and (14), the HMC can be modeled by
(16)UAB¯UBC¯UCA¯=EABEBCECA+jωLf+2LdIAIBIC−jωLf+2LdIBICIA,where
(17)UAB¯UBC¯UCA¯=63V−3063V−15063V+90.Using (14) and (17), the HMC can be modeled as a conventional three-phase two-level VSC, as shown in Figure 8.Figure 8
Equivalent circuit model of HMC.Figure8 shows the equivalent circuit model of the HMC. The equivalent circuit is composed of an open-circuit voltage source linked in series to the HMC output impedance. The open-circuit voltage of the HMC is shown in Figure 7. As shown in Figure 8, the equivalent interface inductor of the HMC is 2Ld+Lf. This feature is advantageous, because the filtering inductor Lf could be minimalized or even removed. However, this can be inopportune when the line transformer has sufficient inductance for filtering the output current. In this case, essentially if the HMC supplies reactive power, the output impedance causes a voltage drop, decreasing the power capability. Thus, it might be better to minimize the interconnected total inductances. However, the inductances are essential to suppress the circulating current within the HMC. The aforementioned decision shows that the interconnected inductance value selection is not a simple process, where the inductance should be large enough to minimize the circulating current but small enough to avoid an extra voltage drop. In the following section, this feature is investigated using multiwinding coupled inductors.


## 4. The Role of the Coupled Inductors

The function of the multiwinding coupled inductors within the HMC is investigated in this section. The HMC performance analysis can be demonstrated based on the function of the coupled inductors. The six VSC modules are interconnected to each other, creating a hexagon as explained in Section3.4. However, instead of using six inductors, multiwinding coupled inductors are used to interconnect the VSC modules.


### 4.1. HMC with Three Two-Winding Coupled Inductors

Two inductors with an equal number of turns are coupled together; that is, the input to one side will produce an output on both, as shown in Figure9. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of the primary and secondary windings are the same (e.g., Lp=Ls=LB). Applying the voltage relationships of coupled inductors, the following equations are given:
(18)vb1b2=L11dib1dt+L12dib5dt,and
(19)vb5b4=L22dib5dt+L21dib1dt,where L22=L11=LB+Lm, Lm=L21=L12, ib1=ib1+icir, and ib5=−ib1+icir, and where ixx=a,b,c is the line current, LB is the leakage inductance of the inductor windings, which is expected to be identical for all windings, and Lm is the magnetizing inductance of the coupled inductor.Figure 9
Equivalent circuit of the two-winding coupled inductor.It is obvious in Figure10 that the voltage drop on one of the coupled inductors is
(20)v=vb1b2+vb5b4.Figure 10
The HMC with three two-winding coupled inductors.Using (18) and (19), the voltage across the coupled inductor can be expressed as
(21)v=LB+Lmdib1+icirdt+Lmd−ib1+icirdt+LB+Lmd−ib1+icirdt+Lmdib1+icirdt,v=2LB+4Lmdicirdt.Therefore, the total voltage drop across the three coupled inductors can be expressed as(22)vtotal=32LB+4Lmdicirdt=6LB+12Lmdicirdt.


### 4.2. HMC with Two Three-Winding Coupled Inductors

Three inductors with an equal number of turns are coupled together as shown in Figure11. Writing the voltage equations of one coupled inductor with winding inductance Z1, the following equations are given:
(23)va5a6=Ladiadt+Labdibdt+Lacdicdt,vb1b2=Lbdibdt+Lbadiadt+Lbcdicdt,vc3c4=Lcdicdt+Lcadiadt+Lcbdibdt,where
(24)La=Lb=Lc=L+Lm,Lm=Lab=Lba=Lbc=Lcb=Lca=Lac,and
(25)ia=ia5+icir,ib=ib1+icir,ic=ic3+icir,where ixx=a,b,c is the line current, L is the leakage inductance of the inductor windings, and Lm is the magnetizing inductance. It is obvious in Figure 12 that the voltage drop on the coupled inductor is
(26)v=va5a6+vb1b2+vc3c4.Figure 11
Equivalent circuit of the three-winding coupled inductor.Figure 12
The HMC with two three-winding coupled inductors.Using (23), the voltage drop across one of the coupled inductor can be expressed as
(27)v=L+Lmdia5+icirdt+Lmdib1+icirdt+Lmdic3+icirdt+L+Lmdib1+icirdt+Lmdia5+icirdt+Lmdic3+icirdt+L+Lmdic3+icirdt+Lmdia5+icirdt+Lmdib1+icirdt,v=L+Lmdia5+ib1+ic3dt+3L+Lmdicirdt+2Lmdia5+ib1+ic3dt+6Lmdicirdt.The voltage drop across the coupled inductor windings is(28)v=L+3Lmdia5+ib1+ic3dt+3L+9Lmdicirdt.Under balanced conditions,(29)ia5+ib1+ic3=ia+ib+ic=0.Accordingly, by applying (28) and (29), the voltage drop on the coupled inductor windings is
(30)v=3L+9Lmdicirdt.The total voltage drop across the two coupled inductors can be expressed as(31)vtotal=6LB+18Lmdicirdt.


### 4.3. HMC with One Six-Winding Coupled Inductor

Instead of using six inductances, one six-winding coupled inductor can be used to interconnect the VSC modules. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of all the windings are the same as shown in Figure13. The HMC with one six-winding coupled inductor is shown in Figure 14. Applying the voltage equations of the coupled inductor, the following equations are given:
(32)va5a6=Ladia5dt+Lmdib1+dic3+dib5+dia3+dic1dt,vb1b2=Lbdib1dt+Lmdia5+dic3+dib5+dia3+dic1dt,vc3c4=Lcdic3dt+Lmdia5+dib1+dib5+dia3+dic1dt,va3a2=Ladia3dt+Lmdib1+dic3+dib5+dia5+dic1dt,vb5b4=Lbdib5dt+Lmdia5+dic3+dib1+dia3+dic1dt,vc1c6=Lcdic1dt+Lmdia5+dib1+dib5+dia3+dic3dt.Since,(33)ia=ia5=ia3=ia1=−ia2=−ia4=−ia6,ib=ib5=ib3=ib1=−ib2=−ib4=−ib6,ic=ic5=ic3=ic1=−ic2=−ic4=−ic6.Figure 13
Equivalent circuit of the six-winding coupled inductor.Figure 14
The HMC with one six-winding coupled inductor.The voltage drop equations across each coupled inductor can be expressed as(34)va5a6=va3a2=Ladiadt+Lmdia+2dib+2dicdt,vb1b2=vb5b4=Lbdibdt+Lmdib+2dia+2dicdt,vc3c4=vc1c6=Lcdibdt+Lmdic+2dia+2dibdt,where
(35)La=Lb=Lc=L+Lm,and
(36)ia=ia+icir,ib=ib+icir,ic=ic+icir.Using (34), the voltage drop across the coupled inductor can be expressed as
(37)v=va5a6+vb1b2+vc3c4+va3a2+vb5b4+vc1c6,(38)v=2va5a6+2vb1b2+2vc3c4,(39)v=6L+36Lmdicirdt+2L+12Lmdia+ib+icdt.Under balanced conditions,(40)ia+ib+ic=0.Consequently, by applying (39) and (40), the voltage drop on the coupled inductor windings is
(41)v=6L+36Lmdicirdt.It is considered that the magnetizing inductance is much higher than the leakage inductance. Therefore, neglecting the leakage inductance, the two-winding coupled inductor imposes twelve times the magnetizing inductance for the circulating current, while the impedance to the circulating current using three-winding coupled inductors is eighteen as much as the magnetizing inductance. Thanks to the six-winding coupled inductor, the impedance on the circulating current path is thirty six times the magnetizing inductance.In a general way ifk is the number of the coupled inductor windings, the impedance to the circulating current can be expressed as
(42)Zcir=2πf6kLm.


## 4.1. HMC with Three Two-Winding Coupled Inductors

Two inductors with an equal number of turns are coupled together; that is, the input to one side will produce an output on both, as shown in Figure9. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of the primary and secondary windings are the same (e.g., Lp=Ls=LB). Applying the voltage relationships of coupled inductors, the following equations are given:
(18)vb1b2=L11dib1dt+L12dib5dt,and
(19)vb5b4=L22dib5dt+L21dib1dt,where L22=L11=LB+Lm, Lm=L21=L12, ib1=ib1+icir, and ib5=−ib1+icir, and where ixx=a,b,c is the line current, LB is the leakage inductance of the inductor windings, which is expected to be identical for all windings, and Lm is the magnetizing inductance of the coupled inductor.Figure 9
Equivalent circuit of the two-winding coupled inductor.It is obvious in Figure10 that the voltage drop on one of the coupled inductors is
(20)v=vb1b2+vb5b4.Figure 10
The HMC with three two-winding coupled inductors.Using (18) and (19), the voltage across the coupled inductor can be expressed as
(21)v=LB+Lmdib1+icirdt+Lmd−ib1+icirdt+LB+Lmd−ib1+icirdt+Lmdib1+icirdt,v=2LB+4Lmdicirdt.Therefore, the total voltage drop across the three coupled inductors can be expressed as(22)vtotal=32LB+4Lmdicirdt=6LB+12Lmdicirdt.


## 4.2. HMC with Two Three-Winding Coupled Inductors

Three inductors with an equal number of turns are coupled together as shown in Figure11. Writing the voltage equations of one coupled inductor with winding inductance Z1, the following equations are given:
(23)va5a6=Ladiadt+Labdibdt+Lacdicdt,vb1b2=Lbdibdt+Lbadiadt+Lbcdicdt,vc3c4=Lcdicdt+Lcadiadt+Lcbdibdt,where
(24)La=Lb=Lc=L+Lm,Lm=Lab=Lba=Lbc=Lcb=Lca=Lac,and
(25)ia=ia5+icir,ib=ib1+icir,ic=ic3+icir,where ixx=a,b,c is the line current, L is the leakage inductance of the inductor windings, and Lm is the magnetizing inductance. It is obvious in Figure 12 that the voltage drop on the coupled inductor is
(26)v=va5a6+vb1b2+vc3c4.Figure 11
Equivalent circuit of the three-winding coupled inductor.Figure 12
The HMC with two three-winding coupled inductors.Using (23), the voltage drop across one of the coupled inductor can be expressed as
(27)v=L+Lmdia5+icirdt+Lmdib1+icirdt+Lmdic3+icirdt+L+Lmdib1+icirdt+Lmdia5+icirdt+Lmdic3+icirdt+L+Lmdic3+icirdt+Lmdia5+icirdt+Lmdib1+icirdt,v=L+Lmdia5+ib1+ic3dt+3L+Lmdicirdt+2Lmdia5+ib1+ic3dt+6Lmdicirdt.The voltage drop across the coupled inductor windings is(28)v=L+3Lmdia5+ib1+ic3dt+3L+9Lmdicirdt.Under balanced conditions,(29)ia5+ib1+ic3=ia+ib+ic=0.Accordingly, by applying (28) and (29), the voltage drop on the coupled inductor windings is
(30)v=3L+9Lmdicirdt.The total voltage drop across the two coupled inductors can be expressed as(31)vtotal=6LB+18Lmdicirdt.


## 4.3. HMC with One Six-Winding Coupled Inductor

Instead of using six inductances, one six-winding coupled inductor can be used to interconnect the VSC modules. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of all the windings are the same as shown in Figure13. The HMC with one six-winding coupled inductor is shown in Figure 14. Applying the voltage equations of the coupled inductor, the following equations are given:
(32)va5a6=Ladia5dt+Lmdib1+dic3+dib5+dia3+dic1dt,vb1b2=Lbdib1dt+Lmdia5+dic3+dib5+dia3+dic1dt,vc3c4=Lcdic3dt+Lmdia5+dib1+dib5+dia3+dic1dt,va3a2=Ladia3dt+Lmdib1+dic3+dib5+dia5+dic1dt,vb5b4=Lbdib5dt+Lmdia5+dic3+dib1+dia3+dic1dt,vc1c6=Lcdic1dt+Lmdia5+dib1+dib5+dia3+dic3dt.Since,(33)ia=ia5=ia3=ia1=−ia2=−ia4=−ia6,ib=ib5=ib3=ib1=−ib2=−ib4=−ib6,ic=ic5=ic3=ic1=−ic2=−ic4=−ic6.Figure 13
Equivalent circuit of the six-winding coupled inductor.Figure 14
The HMC with one six-winding coupled inductor.The voltage drop equations across each coupled inductor can be expressed as(34)va5a6=va3a2=Ladiadt+Lmdia+2dib+2dicdt,vb1b2=vb5b4=Lbdibdt+Lmdib+2dia+2dicdt,vc3c4=vc1c6=Lcdibdt+Lmdic+2dia+2dibdt,where
(35)La=Lb=Lc=L+Lm,and
(36)ia=ia+icir,ib=ib+icir,ic=ic+icir.Using (34), the voltage drop across the coupled inductor can be expressed as
(37)v=va5a6+vb1b2+vc3c4+va3a2+vb5b4+vc1c6,(38)v=2va5a6+2vb1b2+2vc3c4,(39)v=6L+36Lmdicirdt+2L+12Lmdia+ib+icdt.Under balanced conditions,(40)ia+ib+ic=0.Consequently, by applying (39) and (40), the voltage drop on the coupled inductor windings is
(41)v=6L+36Lmdicirdt.It is considered that the magnetizing inductance is much higher than the leakage inductance. Therefore, neglecting the leakage inductance, the two-winding coupled inductor imposes twelve times the magnetizing inductance for the circulating current, while the impedance to the circulating current using three-winding coupled inductors is eighteen as much as the magnetizing inductance. Thanks to the six-winding coupled inductor, the impedance on the circulating current path is thirty six times the magnetizing inductance.In a general way ifk is the number of the coupled inductor windings, the impedance to the circulating current can be expressed as
(42)Zcir=2πf6kLm.


## 5. Control Scheme and Modulation Strategy

### 5.1. Control Scheme

The main function of the controller is to generate reference currents such that the proposed converter only provides available active power from the DC links to the grid at the point of common coupling (PCC) [22, 23]. Using the equivalent circuit model presented in Figure 7 and applying Kirchhoff’s voltage and current laws at the PCC, the following two equations in the abc frame can be obtained:
(43)EAO′EBO′ECO′=UAO¯UBO¯UCO¯+LfdddtIAIBIC+UO′O,(44)CeqdUDC_avdt=SASBSCIAIBIC,where Lfd=Lf+Ld is the equivalent interface inductor and SA,SB, and SC represent the switching states of the equivalent circuit model under balanced conditions.Assuming that the voltages are balanced and the zero-sequence component is zero, the voltage between the neutral virtual point of equivalent circuit model (O) and the grid neutral point (O′) is given by
(45)UO′O=−UAO¯+UBO¯+UCO¯3,(46)UAO¯UBO¯UCO¯=UDC_avSASBSC.Substituting (46) into (43) and (44), the following relation is obtained:
(47)ddtIAIBIC=1LfdEAO′EBO′ECO′−UDCavLfdSASBSC−13SASBSC111.The dynamic model in the abc frame of the HMC equivalent circuit is represented by (47). The switching state functions, di (i=A, B, C), are defined as
(48)dAdBdC=SASBSC−13SASBSC111.The dynamic model of the equivalent circuit model in the abc frame is achieved by combining (47) and (48) in the following equation:
(49)LfdddtIAIBIC=EAO′EBO′ECO′−UDCavdAdBdC.The DC side differential equation can be written as(50)dUDC_avdt=1CeqIdc=dAdBdCIAIBIC,(51)dUDC_avdt=1Ceq2dA+dBIA+1CeqdA+2dBIB.It can be seen that the model represented by (49) and (51) is time varying. Thus, to facilitate the control algorithm implementation, the model can be expressed in the synchronous reference frame rotating at constant frequency ω. The corresponding conversion matrix is
(52)Cdqabc=23cosθcosθ−2π3cosθ−4π3−sinθ−sinθ−2π3−sinθ−4π3,whereθ=ωt.Applying the coordinate transformation to (49) we obtain
(53)LfdddtIdIq=EdEq+LfωIq−Id−UDCavdddq.Similarly, applying this transformation to (50) we obtain
(54)CeqdUDC_avdt=ddId+dqIq.The obtained model, represented by (53) and (54), is nonlinear owing to the product between the state variables (i.e., Id,Iq,andUDC_av) and the inputs (i.e., ddanddq).Figure11 shows the control principle of the HMC. Because the proposed configuration is a three-wire system, only two phase currents are required to be measured. Currents IAandIB are measured and converted to the dq0 frame to obtain the corresponding currents IdandIq. Accordingly, (53) is rewritten as follows:
(55)Lfduduq=EdEq+LfdωIq−Id−UDCavdddq,where uduq=d/dtIdIq.Demonstrated that currentsId and Iq can be controlled separately by acting on inputs ud and uq, respectively. Hence, the controller is designed using the following expressions:
(56)ud=kpid~+ki∫id~dt,uq=kpiq~+ki∫iq~dt,where id~=Id_REF−Id and iq~=Iq_REF−Iq are the current error signals, with Id_REF and Iq_REF being the reference values for currents Id and Iq, respectively.Using (55), the current control law is given by the following expression:
(57)dddq=LfdωUDC_avIq−Id+1UDC_avEdEq−LfdUDCavuduq.Thed-axis reference current (Id_REF) is produced using the DC-link voltage controller, while the q-axis reference current (Iq_REF) is taken as the load q-axis current (ILq). To aim for a unity power factor, the Iq_REF is set to zero (i.e., iq_REF=0). The active power exchange between the DC links and the grid is proportional to the direct-axis current Id and can be expressed as
(58)Pdc=32EdId+EqIq=32EdId.(58) it is shown that direct-axis current Id is responsible for maintaining the DC-link voltages at a desired value. Thus, using (54), one deduces that
(59)CeqdUDC_avdt=ddId=udc.Consequently, the active current is(60)Id=udcdd=udcVDC_avddVDC_av.Assuming that the current loop is perfect and the HMC works under balanced conditions, the following expressions hold:(61)ddVDC_av=Ed,Id=udcdd=udcVDC_avEd,EdEq=32V^0,where V^ and Ed are the RMS voltage and the direct-axis phase voltage at the PCC, respectively. Thus, the control effort of the DC link voltage loop is given by
(62)Id_REF=udcdd=23VDC_avV^udc.To control the DC-link voltage, a PI controller is used, which is expressed as(63)udc=kpdcvdc~+kidc∫vdc~dt,where vdc~=VDC_REF−VDC_av is the DC-link voltage error, whereas VDC_REF and VDC_av are the DC-link reference and average voltages, respectively.


### 5.2. DC-Link Voltage Controller

With reference to Figure3, the HMC is based on a symmetric configuration, having six converters with identical power capabilities that are supplied by six equal PV strings. The PV strings are directly connected to each converter, Vpv1=Vpv2=Vpv3=Vpv4=Vpv5=Vpv6=Vpv, the MPPT must be achieved by the converters, and the DC-link voltages continuously fluctuate. Because the PV strings are supposed to be identical, being created by a single PV string divided into six identical parts, a single MPPT algorithm can be considered. For this reason, the same DC-link voltage reference for the six converters has been considered. The DC-link voltage reference is compared to the sum of the actual six DC-link voltages, and the error is passed through a PI controller to determine the control parameterudc.


### 5.3. Modulation Strategy

To obtain an output voltage with low total harmonic distortion (THD), a multicarrier phase-shifted PWM (PS-PWM) switching strategy [24, 25] is implemented to drive each IGBT in the HMC. Optimum harmonic cancellation is accomplished by shifting each carrier cell by 2πTs/3T in sequence, where Ts is the switching time and T is the cycle modulation time. A reference signal of 50 Hz is generated using the control algorithm represented in Figure 15, with the switching frequency fixed at 2500 Hz.Figure 15
HMC control block diagram.Figure16 shows the relationship between the modulation waveforms and the three groups of carriers within the HMC. As shown in Figure 16, triangular carriers (i.e., carrier #1, carrier #2, and carrier #3) are phase-shifted 120° to each other and directly compared with the modulation signals to drive the IGBTs within module #1, module #2, and module #3. In order to generate the switching signals used to drive the IGBTs within module #2, module #4, and module #6, the modulation signals are inverted and then compared with the triangular carriers.Figure 16
Proposed PS-PWM diagram.


## 5.1. Control Scheme

The main function of the controller is to generate reference currents such that the proposed converter only provides available active power from the DC links to the grid at the point of common coupling (PCC) [22, 23]. Using the equivalent circuit model presented in Figure 7 and applying Kirchhoff’s voltage and current laws at the PCC, the following two equations in the abc frame can be obtained:
(43)EAO′EBO′ECO′=UAO¯UBO¯UCO¯+LfdddtIAIBIC+UO′O,(44)CeqdUDC_avdt=SASBSCIAIBIC,where Lfd=Lf+Ld is the equivalent interface inductor and SA,SB, and SC represent the switching states of the equivalent circuit model under balanced conditions.Assuming that the voltages are balanced and the zero-sequence component is zero, the voltage between the neutral virtual point of equivalent circuit model (O) and the grid neutral point (O′) is given by
(45)UO′O=−UAO¯+UBO¯+UCO¯3,(46)UAO¯UBO¯UCO¯=UDC_avSASBSC.Substituting (46) into (43) and (44), the following relation is obtained:
(47)ddtIAIBIC=1LfdEAO′EBO′ECO′−UDCavLfdSASBSC−13SASBSC111.The dynamic model in the abc frame of the HMC equivalent circuit is represented by (47). The switching state functions, di (i=A, B, C), are defined as
(48)dAdBdC=SASBSC−13SASBSC111.The dynamic model of the equivalent circuit model in the abc frame is achieved by combining (47) and (48) in the following equation:
(49)LfdddtIAIBIC=EAO′EBO′ECO′−UDCavdAdBdC.The DC side differential equation can be written as(50)dUDC_avdt=1CeqIdc=dAdBdCIAIBIC,(51)dUDC_avdt=1Ceq2dA+dBIA+1CeqdA+2dBIB.It can be seen that the model represented by (49) and (51) is time varying. Thus, to facilitate the control algorithm implementation, the model can be expressed in the synchronous reference frame rotating at constant frequency ω. The corresponding conversion matrix is
(52)Cdqabc=23cosθcosθ−2π3cosθ−4π3−sinθ−sinθ−2π3−sinθ−4π3,whereθ=ωt.Applying the coordinate transformation to (49) we obtain
(53)LfdddtIdIq=EdEq+LfωIq−Id−UDCavdddq.Similarly, applying this transformation to (50) we obtain
(54)CeqdUDC_avdt=ddId+dqIq.The obtained model, represented by (53) and (54), is nonlinear owing to the product between the state variables (i.e., Id,Iq,andUDC_av) and the inputs (i.e., ddanddq).Figure11 shows the control principle of the HMC. Because the proposed configuration is a three-wire system, only two phase currents are required to be measured. Currents IAandIB are measured and converted to the dq0 frame to obtain the corresponding currents IdandIq. Accordingly, (53) is rewritten as follows:
(55)Lfduduq=EdEq+LfdωIq−Id−UDCavdddq,where uduq=d/dtIdIq.Demonstrated that currentsId and Iq can be controlled separately by acting on inputs ud and uq, respectively. Hence, the controller is designed using the following expressions:
(56)ud=kpid~+ki∫id~dt,uq=kpiq~+ki∫iq~dt,where id~=Id_REF−Id and iq~=Iq_REF−Iq are the current error signals, with Id_REF and Iq_REF being the reference values for currents Id and Iq, respectively.Using (55), the current control law is given by the following expression:
(57)dddq=LfdωUDC_avIq−Id+1UDC_avEdEq−LfdUDCavuduq.Thed-axis reference current (Id_REF) is produced using the DC-link voltage controller, while the q-axis reference current (Iq_REF) is taken as the load q-axis current (ILq). To aim for a unity power factor, the Iq_REF is set to zero (i.e., iq_REF=0). The active power exchange between the DC links and the grid is proportional to the direct-axis current Id and can be expressed as
(58)Pdc=32EdId+EqIq=32EdId.(58) it is shown that direct-axis current Id is responsible for maintaining the DC-link voltages at a desired value. Thus, using (54), one deduces that
(59)CeqdUDC_avdt=ddId=udc.Consequently, the active current is(60)Id=udcdd=udcVDC_avddVDC_av.Assuming that the current loop is perfect and the HMC works under balanced conditions, the following expressions hold:(61)ddVDC_av=Ed,Id=udcdd=udcVDC_avEd,EdEq=32V^0,where V^ and Ed are the RMS voltage and the direct-axis phase voltage at the PCC, respectively. Thus, the control effort of the DC link voltage loop is given by
(62)Id_REF=udcdd=23VDC_avV^udc.To control the DC-link voltage, a PI controller is used, which is expressed as(63)udc=kpdcvdc~+kidc∫vdc~dt,where vdc~=VDC_REF−VDC_av is the DC-link voltage error, whereas VDC_REF and VDC_av are the DC-link reference and average voltages, respectively.


## 5.2. DC-Link Voltage Controller

With reference to Figure3, the HMC is based on a symmetric configuration, having six converters with identical power capabilities that are supplied by six equal PV strings. The PV strings are directly connected to each converter, Vpv1=Vpv2=Vpv3=Vpv4=Vpv5=Vpv6=Vpv, the MPPT must be achieved by the converters, and the DC-link voltages continuously fluctuate. Because the PV strings are supposed to be identical, being created by a single PV string divided into six identical parts, a single MPPT algorithm can be considered. For this reason, the same DC-link voltage reference for the six converters has been considered. The DC-link voltage reference is compared to the sum of the actual six DC-link voltages, and the error is passed through a PI controller to determine the control parameterudc.


## 5.3. Modulation Strategy

To obtain an output voltage with low total harmonic distortion (THD), a multicarrier phase-shifted PWM (PS-PWM) switching strategy [24, 25] is implemented to drive each IGBT in the HMC. Optimum harmonic cancellation is accomplished by shifting each carrier cell by 2πTs/3T in sequence, where Ts is the switching time and T is the cycle modulation time. A reference signal of 50 Hz is generated using the control algorithm represented in Figure 15, with the switching frequency fixed at 2500 Hz.Figure 15
HMC control block diagram.Figure16 shows the relationship between the modulation waveforms and the three groups of carriers within the HMC. As shown in Figure 16, triangular carriers (i.e., carrier #1, carrier #2, and carrier #3) are phase-shifted 120° to each other and directly compared with the modulation signals to drive the IGBTs within module #1, module #2, and module #3. In order to generate the switching signals used to drive the IGBTs within module #2, module #4, and module #6, the modulation signals are inverted and then compared with the triangular carriers.Figure 16
Proposed PS-PWM diagram.


## 6. Simulation Results

In order to demonstrate the performance of the HMC and its control algorithm, the complete grid-connected PV system was simulated using the MATLAB/Simulink environment. The technical characteristics of the system parameters and the coupled inductors are provided in Table2. In order to validate the performance of the HMC, a PS-PWM technique has been implemented, as shown in Figure 16. The output active and reactive power supplies in response to the fluctuations in solar radiation value are shown in the following subsections.Table 2
General data of the proposed system.
Parameter
Value
System parameters
Nominal power
P = 2.5 kVA
Phase voltage and frequency
Vph = 110 V (rms), fsys = 50 Hz, fswiching = 2500 Hz
Output current
Iout = 8.8 A (rms)
DC bus voltages
Vdc1 = Vdc2 = Vdc3 = Vdc4 = Vdc5 = Vdc6 = 52.6 V
Open-end winding transformers
Vpri/Vsec = 1, ftran = 50 Hz
Current controller parameters
Kp = 55, Ki = 0.001
Voltage controller parameters
Kp = 3, Ki = 20
Coupled inductor parameters
RMS voltage, current
20 V (rms), 8.8 A (rms)
Magnetizing inductance
3.2 mH
Leakage inductance
150μH


### 6.1. Performance Analysis

The seventeen-level phase voltages of the HMC are generated at the steady state and shown in Figure17. It can be seen that the voltages are balanced, as shown in the voltage phasor diagram. Figure 18 demonstrates the currents of modules 1, 3, and 5, and the currents of modules 2, 4, and 6. Given the configuration of the three-phase HMC, the line currents of every VSC module within the HMC are symmetrical, hence verifying the relationship in (5).Figure 17
Three-phase output voltage of the HMC.Figure 18
Currents inside the HMC.From this figure, it can also be acknowledged that the direction of the currents of modules 2, 4, and 6 is reversed from those of modules 1, 3, and 5. Moreover, the currents inside every VSC module within the HMC are the same as the output currents.The THD values of the output voltage and current are calculated using the following equation.(64)%THDx=100∑h≠1xshxs12,where the subscript x indicates the THD in the signal (voltage or current), xs1 is the fundamental component, xsh is the component at the h harmonic frequency.The harmonic spectra of the output voltage and current are shown in Figures19 and 20, respectively. The THD of the HMC current is 0.68%, which is fewer than 5% and meets the power quality standard. Using the suggested modulation technique, the highest harmonic family of the phase voltage appears at the band of the 100th harmonic order. Consequently, the effective switching frequency of the phase voltage is two times higher than the switching frequency.Figure 19
THD of the HMC output voltage.Figure 20
THD of the HMC output current.As a comparison, the HMC using two-, three-, and six-winding coupled inductors are simulated under the conditions. The magnetizing inductance of each coupled inductor is 3.5 mH. According to the current from the equations provided in Section4, if the voltage of the DC links is unbalanced the circulating current will be introduced on the converter output currents. Therefore, to intentionally produce a circulating current, the DC link voltage of Module 3 is decreased from 52.6 to 26.3 V. The created loop voltage is computed using the following equation:
(65)Vloop,rms=32252.6−26.3ma.Thus, using (42) and (65), the circulating current inside the HMC is expressed as
(66)Icir,peak=324π52.6−26.3mafkLm.Investigations with different magnetizing inductance levels have been carried. Figure21 shows the theoretical values of the peak circulating currents with different magnetizing inductance levels. Figure 22 shows the simulation results with a relatively large magnetizing inductance (3.5 mH). The circulating currents are found by measuring the difference between the output currents of module 1 and module 3. As shown in the waveforms, the circulating currents inside the HMC with two-, three-, and six-winding coupled inductors are 2.3, 1.5, and 0.76 A peaks as calculated in (66). From the waveforms in Figure 22, using a six-winding coupled inductor, the circulating current is efficiently minimized.Figure 21
Peak circulating current versus magnetizing inductance.Figure 22
Circulating current inside the HMC with multiwinding coupled inductors.In order to achieve the same circulating currents (e.g., 0.5 A), the magnetizing inductance of the two-winding, three-winding and six-winding coupled inductors should be increased to 17.5, 11.6, and 6 mH, respectively. As shown in Figure23, the waveform matches the theoretical analysis in Figure 21. The output current of module 1, module 3, and the circulating current inside the HMC with two-winding, three-winding, and six-winding coupled inductors under an unbalanced DC link voltage are shown in Figures 24, 25, and 26, respectively. The waveforms show the simulation result with an extremely small magnetizing inductance of 0.5 mH, which is 1/7 of that in Figure 22. The simulations validated the quantitative relationship that the circulating current will be too big to maintain a normal operation if the coupled inductors with low magnetizing inductance are used. As shown, the current ia3 nearly doubles the circulating current, which evidently proves that the HMC is sensitive to an unbalanced DC voltage. The waveforms match the theoretical analysis in Figure 21.Figure 23
Circulating current inside the HMC.Figure 24
Circulating current with two-winding coupled inductors with a magnetizing inductance of 0.5 mH.Figure 25
Circulating current with three-winding coupled inductors with a magnetizing inductance of 0.5 mH.Figure 26
Circulating current with six-winding coupled inductors with a magnetizing inductance of 0.5 mH.


### 6.2. Active Power Variation

The system is tested for different solar radiation conditions and the results are shown in Table3. Figure 27 shows the solar radiation variation which is considered in this study. The output current and grid voltage variation in response to the changes in solar radiation is shown in Figure 28. The transient behavior of the total DC-link voltage is presented in Figure 29. The parameters are successfully attuned by the controller to keep the DC-link voltage at the desired level of 52.6 V. Fluctuations are observed in the DC-links, due to the step variations in the solar radiation. Nevertheless, the controller brings the voltage to the reference level within 0.02 s.Table 3
Total harmonic distortion under different solar radiation conditions.
Solar radiation
THD (%)
Fundamental value
Voltage
Current
Voltage (V)
Current (A)
1000 (W/m2)
13.64
0.68
184.0
8.80
600 (W/m2)
13.65
1.15
181.8
5.20
400 (W/m2)
12.78
2.80
180.9
3.44
0 (W/m2)
23.94
0.35
148.3
9.98Figure 27
Step changes in solar radiation.Figure 28
Response to changes in solar irradiance in the hexagram converter.Figure 29
Variations in DC-link voltages in response to the changes in solar radiation.The active power supplied by the HMC is directly proportional to the magnitude of direct axis,id, as indicated by (38). The solar radiation at all DC links is reduced by 40% at 0.25 s. The fluctuation in solar radiation is imitated through the reduction in the direct-axis current by approximately 40%, in a step. It decreased to 5.216 A from the initial value of 8.823 A, as shown in Figure 30. This reduction in id was to keep the DC-link voltage at the reference level by decreasing the output power drawn from the PV. Moreover, to guarantee maximum utilization of the PV system, the quadrature axis current was kept at zero. At 0.45 s, the solar radiation is further reduced to 400 W/m2. The direct-axis current is reduced to 3.267 A from the original value of 5.216 A. Later, the direct-axis current is increased to 8.823 A, because of the increment in the solar radiation.Figure 30
Response to changes in solar radiation. (a) Output active and reactive power, and (b) direct and quadrature axis current.
(a)
(b)


### 6.3. Reactive Power Compensation

The described control algorithm permits the HMC to act as DSTATCOM in the absence of solar radiation. The output reactive power can be calculated as(67)Qout=−32IqVd.In this condition, the reactive power is increased by 2700 VAR in a step, in the absence of solar radiation. The influence, for the DSTATCOM mode operation, is presented in Figure31.Figure 31
System performance with the increase in the reactive power supply in DSTATCOM mode.The direct-axis current is found to be zero, indicating the fact that no active power is being transferred to the grid in the absence of solar radiation. However, due to the step change in reactive power, the capacitors consume current from the grid. It is found that the direct-axis current takes 0.02 s to stabilize. The nature of the fluctuation in the DC-link voltage is depicted in Figure29. The DC-link controller effectively keeps the DC-link voltage by regulating the power flow through the capacitor. Hence, it can be mentioned that the HMC effectively operates as DSTATCOM in the absence of solar radiation. Figure 31 shows the source voltage and converter output current in DSTATCOM mode. The waveforms show that the output current increases after the reactive power command comes at 0.8 s. Moreover, Figure 31 approves that the phase difference of the converter output current with grid voltage is 90° in this mode of operation.


## 6.1. Performance Analysis

The seventeen-level phase voltages of the HMC are generated at the steady state and shown in Figure17. It can be seen that the voltages are balanced, as shown in the voltage phasor diagram. Figure 18 demonstrates the currents of modules 1, 3, and 5, and the currents of modules 2, 4, and 6. Given the configuration of the three-phase HMC, the line currents of every VSC module within the HMC are symmetrical, hence verifying the relationship in (5).Figure 17
Three-phase output voltage of the HMC.Figure 18
Currents inside the HMC.From this figure, it can also be acknowledged that the direction of the currents of modules 2, 4, and 6 is reversed from those of modules 1, 3, and 5. Moreover, the currents inside every VSC module within the HMC are the same as the output currents.The THD values of the output voltage and current are calculated using the following equation.(64)%THDx=100∑h≠1xshxs12,where the subscript x indicates the THD in the signal (voltage or current), xs1 is the fundamental component, xsh is the component at the h harmonic frequency.The harmonic spectra of the output voltage and current are shown in Figures19 and 20, respectively. The THD of the HMC current is 0.68%, which is fewer than 5% and meets the power quality standard. Using the suggested modulation technique, the highest harmonic family of the phase voltage appears at the band of the 100th harmonic order. Consequently, the effective switching frequency of the phase voltage is two times higher than the switching frequency.Figure 19
THD of the HMC output voltage.Figure 20
THD of the HMC output current.As a comparison, the HMC using two-, three-, and six-winding coupled inductors are simulated under the conditions. The magnetizing inductance of each coupled inductor is 3.5 mH. According to the current from the equations provided in Section4, if the voltage of the DC links is unbalanced the circulating current will be introduced on the converter output currents. Therefore, to intentionally produce a circulating current, the DC link voltage of Module 3 is decreased from 52.6 to 26.3 V. The created loop voltage is computed using the following equation:
(65)Vloop,rms=32252.6−26.3ma.Thus, using (42) and (65), the circulating current inside the HMC is expressed as
(66)Icir,peak=324π52.6−26.3mafkLm.Investigations with different magnetizing inductance levels have been carried. Figure21 shows the theoretical values of the peak circulating currents with different magnetizing inductance levels. Figure 22 shows the simulation results with a relatively large magnetizing inductance (3.5 mH). The circulating currents are found by measuring the difference between the output currents of module 1 and module 3. As shown in the waveforms, the circulating currents inside the HMC with two-, three-, and six-winding coupled inductors are 2.3, 1.5, and 0.76 A peaks as calculated in (66). From the waveforms in Figure 22, using a six-winding coupled inductor, the circulating current is efficiently minimized.Figure 21
Peak circulating current versus magnetizing inductance.Figure 22
Circulating current inside the HMC with multiwinding coupled inductors.In order to achieve the same circulating currents (e.g., 0.5 A), the magnetizing inductance of the two-winding, three-winding and six-winding coupled inductors should be increased to 17.5, 11.6, and 6 mH, respectively. As shown in Figure23, the waveform matches the theoretical analysis in Figure 21. The output current of module 1, module 3, and the circulating current inside the HMC with two-winding, three-winding, and six-winding coupled inductors under an unbalanced DC link voltage are shown in Figures 24, 25, and 26, respectively. The waveforms show the simulation result with an extremely small magnetizing inductance of 0.5 mH, which is 1/7 of that in Figure 22. The simulations validated the quantitative relationship that the circulating current will be too big to maintain a normal operation if the coupled inductors with low magnetizing inductance are used. As shown, the current ia3 nearly doubles the circulating current, which evidently proves that the HMC is sensitive to an unbalanced DC voltage. The waveforms match the theoretical analysis in Figure 21.Figure 23
Circulating current inside the HMC.Figure 24
Circulating current with two-winding coupled inductors with a magnetizing inductance of 0.5 mH.Figure 25
Circulating current with three-winding coupled inductors with a magnetizing inductance of 0.5 mH.Figure 26
Circulating current with six-winding coupled inductors with a magnetizing inductance of 0.5 mH.


## 6.2. Active Power Variation

The system is tested for different solar radiation conditions and the results are shown in Table3. Figure 27 shows the solar radiation variation which is considered in this study. The output current and grid voltage variation in response to the changes in solar radiation is shown in Figure 28. The transient behavior of the total DC-link voltage is presented in Figure 29. The parameters are successfully attuned by the controller to keep the DC-link voltage at the desired level of 52.6 V. Fluctuations are observed in the DC-links, due to the step variations in the solar radiation. Nevertheless, the controller brings the voltage to the reference level within 0.02 s.Table 3
Total harmonic distortion under different solar radiation conditions.
Solar radiation
THD (%)
Fundamental value
Voltage
Current
Voltage (V)
Current (A)
1000 (W/m2)
13.64
0.68
184.0
8.80
600 (W/m2)
13.65
1.15
181.8
5.20
400 (W/m2)
12.78
2.80
180.9
3.44
0 (W/m2)
23.94
0.35
148.3
9.98Figure 27
Step changes in solar radiation.Figure 28
Response to changes in solar irradiance in the hexagram converter.Figure 29
Variations in DC-link voltages in response to the changes in solar radiation.The active power supplied by the HMC is directly proportional to the magnitude of direct axis,id, as indicated by (38). The solar radiation at all DC links is reduced by 40% at 0.25 s. The fluctuation in solar radiation is imitated through the reduction in the direct-axis current by approximately 40%, in a step. It decreased to 5.216 A from the initial value of 8.823 A, as shown in Figure 30. This reduction in id was to keep the DC-link voltage at the reference level by decreasing the output power drawn from the PV. Moreover, to guarantee maximum utilization of the PV system, the quadrature axis current was kept at zero. At 0.45 s, the solar radiation is further reduced to 400 W/m2. The direct-axis current is reduced to 3.267 A from the original value of 5.216 A. Later, the direct-axis current is increased to 8.823 A, because of the increment in the solar radiation.Figure 30
Response to changes in solar radiation. (a) Output active and reactive power, and (b) direct and quadrature axis current.
(a)
(b)


## 6.3. Reactive Power Compensation

The described control algorithm permits the HMC to act as DSTATCOM in the absence of solar radiation. The output reactive power can be calculated as(67)Qout=−32IqVd.In this condition, the reactive power is increased by 2700 VAR in a step, in the absence of solar radiation. The influence, for the DSTATCOM mode operation, is presented in Figure31.Figure 31
System performance with the increase in the reactive power supply in DSTATCOM mode.The direct-axis current is found to be zero, indicating the fact that no active power is being transferred to the grid in the absence of solar radiation. However, due to the step change in reactive power, the capacitors consume current from the grid. It is found that the direct-axis current takes 0.02 s to stabilize. The nature of the fluctuation in the DC-link voltage is depicted in Figure29. The DC-link controller effectively keeps the DC-link voltage by regulating the power flow through the capacitor. Hence, it can be mentioned that the HMC effectively operates as DSTATCOM in the absence of solar radiation. Figure 31 shows the source voltage and converter output current in DSTATCOM mode. The waveforms show that the output current increases after the reactive power command comes at 0.8 s. Moreover, Figure 31 approves that the phase difference of the converter output current with grid voltage is 90° in this mode of operation.


## 7. Conclusion

The HMC for a grid-connected PV system shares many advantages of the CHB but uses fewer switches and reduces the size of the DC-link capacitor. However, the HMC is sensitive to the loop voltage produced by the probable DC-link voltage unbalance. The unbalanced conditions of the DC-link voltages will initiate a line-frequency circulating current. Therefore, to minimize the circulating current as well as to reduce an extra voltage drop on the inductors that affects the terminal voltage, the inductance value should be adequate. The multiwinding coupled inductors are the key to control the circulating current and ensure the proper operation of the HMC. The equivalent circuit model of the HMC configuration is derived to recognize the control scheme. The two-winding coupled inductor imposes twelve times the magnetizing inductance for the circulating current, while the impedance to the circulating current using three-winding coupled inductors is eighteen times as much as the magnetizing inductance. Owing to the six-winding coupled inductor, the impedance on the circulating current path is thirty-six times the magnetizing inductance. The results show the good performance of the control algorithm in both steady state and transient conditions. Moreover, it is interesting to note that in the absence of solar radiation, the controller acts in DSTATCOM mode to supply reactive power to the grid. The performance of the controller under different solar radiation conditions is found to be satisfactory.



---

*Source: 9164528-2018-05-07.xml*