A Steady-State Analysis Method for Modular Multilevel Converters Connected to Permanent Magnet Synchronous Generator-Based Wind Energy Conversion Systems

: Modular multilevel converters (MMCs) have shown great potential in the area of multi-megawatt wind energy conversion system (WECS) based on permanent magnet synchronous generators (PMSGs). However, the studies in this area are few, and most of them refer to the MMC used in high-voltage direct current (HVDC) systems, and hence the characteristics of the PMSG are not considered. This paper proposes a steady-state analysis method for MMCs connected to a PMSG-based WECS. In the proposed method, only the wind speed (operating condition) is required as input, and all the electrical quantities in the MMC, including the amplitudes, phase angles and their harmonics, can be calculated step by step. The analysis method is built on the proposed d-q frame mathematical model. Interactions of electrical quantities between the MMC and PMSG are comprehensively considered. Moreover, a new way to calculate the average switching functions are adopted in order to improve the accuracy of the analysis method. Applications of the proposed method are also presented, which includes the characteristic analysis of capacitor voltage ripples and the capacitor sizing. Finally, the accuracy of the method and the correctness of the analysis are veriﬁed by simulations and experiments.


Introduction
Wind energy, as one of the fastest growing renewable energy sources, is gradually becoming mainstream and competitive with conventional energy sources [1][2][3][4]. The installed capacity of wind energy has increased from 6100 MW in 1996 to 432,883 MW by 2015 [2]. The continuous growth in demand for wind energy promotes the development of larger wind turbines, because the large turbines can capture more wind power with lower installation and maintenance costs, compared to the group of small turbines [5]. The size and capacity of commercial wind turbines have exponentially increased over the past 30 years, and has currently reached the level of 10MWs [2,5].
In multi-megawatt wind energy conversion systems (WECS), permanent magnet synchronous generators (PMSGs) with full-scale power converters have been widely adopted, due to their salient features such as high-power density, no need of excitation, and high efficiency. However, as the power capacity of wind turbines has increased up to 10 MW, the traditional two-and three-level converters

•
The analysis method is built on the proposed d-q frame mathematical model of MMC connected to a PMSG. Interactions of electrical quantities between the MMC and PMSG are comprehensively considered. In the d-q frame mathematical model, the time-varying quantities are transformed into constant quantities, which therefore simplifies the derivation. Due to this, the equivalent resistances in the MMC arms can be considered, and the algebraic solution of non-linear equations can be obtained to calculate the unknowns in the average switching functions.

•
Only the wind speed (operating condition) is required as input, and all the electrical quantities in MMC, including the amplitudes, phase angles and their harmonics, can be calculated step by step, when the parameters, such as the capacitance and inductance in the MMC and the flux linkage in the PMSG, are set down.

•
In order to improve the accuracy of the analysis, a new way to calculate the average switching functions in the mathematic model are adopted. In addition, the calculation results are obtainable at one stroke and no need for iteration. Therefore, the accuracy of analysis can be improved without sacrificing a lot calculation speed. • Capacitor voltage ripple is one of the key problems in the MMC connected to variable-speed machines; this problem is further analyzed based on the proposed analysis method. The analysis results show that the MMC connected to PMSG has different characteristics from the MMC used for motor drives. A capacitor sizing method for MMC connected to PMSG is also proposed in this paper.
The rest of the paper is organized as follows: Section 2 builds the d-q frame mathematical model of a MMC connected to a PMSG. The equations of electrical quantities in the MMC are derived in Section 3. Section 4 presents the new method to calculate the average switching functions, and then provides the flow chart of the proposed steady-state analysis method. In Section 5, applications of the proposed method are presented, which includes the characteristic analysis of capacitor voltage ripples and the capacitor sizing. In Sections 6 and 7, the simulations and experiments are carried out to validate the accuracy of the proposed method. Section 8 provides the conclusions of this work. Figure 1 shows the typical structure of a MMC connected to the PMSG-based WECS. At first, the wind energy is converted into mechanical power by blades. Then, the mechanical power is converted into electrical power via PMSG. Finally, the MMC rectifies the ac power into dc power, and transfers the electrical power from the PMSG to the dc bus.

MMC Connected to the PMSG-Based WECS
Energies 2017, 10, x FOR PEER REVIEW 3 of 30 obtainable at one stroke and no need for iteration. Therefore, the accuracy of analysis can be improved without sacrificing a lot calculation speed.

•
Capacitor voltage ripple is one of the key problems in the MMC connected to variable-speed machines; this problem is further analyzed based on the proposed analysis method. The analysis results show that the MMC connected to PMSG has different characteristics from the MMC used for motor drives. A capacitor sizing method for MMC connected to PMSG is also proposed in this paper.
The rest of the paper is organized as follows: Section 2 builds the d-q frame mathematical model of a MMC connected to a PMSG. The equations of electrical quantities in the MMC are derived in Section 3. Section 4 presents the new method to calculate the average switching functions, and then provides the flow chart of the proposed steady-state analysis method. In Section 5, applications of the proposed method are presented, which includes the characteristic analysis of capacitor voltage ripples and the capacitor sizing. In, Sections 6 and 7, the simulations and experiments are carried out to validate the accuracy of the proposed method. Section 8 provides the conclusions of this work. Figure 1 shows the typical structure of a MMC connected to the PMSG-based WECS. At first, the wind energy is converted into mechanical power by blades. Then, the mechanical power is converted into electrical power via PMSG. Finally, the MMC rectifies the ac power into dc power, and transfers the electrical power from the PMSG to the dc bus. During the wind energy conversion process, there is a strong coupling relationship between PMSG and MMC, which is shown in Figure 1   During the wind energy conversion process, there is a strong coupling relationship between PMSG and MMC, which is shown in Figure 1 by the blue dashed lines. The wind energy (P wind ) determines the input mechanical power (P m ), which then affects mechanical torque (T m ) and rotor speed (ω m ). The rotor speed influences the output voltage of PMSG (u a ), which determines the phase current (i a ). The phase current is part of the arm current (i arm ) in MMC, and it affects the capacitor current (i cap ) and voltage (u cap ). The capacitor voltage determines the output voltage of sub-module (u sub ), which afterwards affects the phase currents by influencing the electromotive force (EMF) of MMC (not shown in Figure 1). Finally, the phase current in return determines the electromagnetic torque (T e ), and then the electromagnetic torque in return affects the rotor speed. Because of this strong coupling relationship, the analysis of MMC connected to the PMSG-based WECS becomes complicated. Hence, in this section, the basic configuration of MMC is described first, and then the mathematical model of MMC connected to the PMSG-based WECS is built in the d-q frame. The established d-q frame mathematical model is the basement of the proposed steady-state analysis method, and it can also further help to understand the coupling relationship.

Basic Configuration of MMC
As shown in Figure 1, the three ac terminals of MMC are connected to the positive and negative dc bus bars through three upper and lower arms. There, i arm.j.up and i arm.j.dn (j = a, b, c, similarly hereinafter) denote the arm currents of upper and lower arms, respectively. u cap and i cap are capacitor voltage and current. C denotes the capacitor in submodules. u sub and i sub are submodule current and voltage. L arm is the arm inductor. I dc and U dc denote the dc-side current and voltage.
Different from the MMC connected to the power grid, the frequency in the MMC connected to PSMG varies with the mechanical angular speed of PMSG rotor and is also affected by the number of pole-pairs in the generator, which can be expressed as: where, ω is the angular speed of MMC; ω r and ω m denote the electrical and mechanical angular speed of the rotor in PMSG; p is the number of pole-pairs. Therefore, it is more convenient to take the d-axis of PMSG rotor flux as the datum line. In other words, the MMC model is built in the coordinate system of PMSG, and the phase angle of the ac-side voltage of MMC is no longer zero, as it used to be.
Taking phase A for example, the ac-side voltage of MMC can be defined as: where, U s and φ are the amplitude and the phase angle of the output voltage. According to the previous papers [17,18], the arm current of MMC is consisted of the dc component, phase current component and circulating current component, which can be expressed as: where: i a (t) = I s cos(ωt + ϕ) i cir.a (t) = I cir cos(2ωt + θ cir ) I s and ϕ are the amplitude and the phase angle of phase current. I cir and θ cir denote the amplitude and the phase angle of 2nd harmonic circulating current.
It should be noted that although there are 4th and higher even-order harmonics in the circulating current, they are usually controlled to be small enough to be neglected [20]. Hence, only 2nd harmonics are considered in this paper.
The upper and lower arm voltages of MMC are defined as the sum of sub-module voltages u sub in their respective arm: u arm.a.dn (t) = where, N is the number of submodules in one arm. From Kirchhoff's law and Equations (2)-(8), the following two equations can be obtained: Then, Equations (11) and (12) can be derived based on Equations (9) and (10): where, the EMF e MMC.a and the common-mode voltage u com.a are defined as: The EMF consists of the fundamental component and the odd-order harmonics. However, the odd-order harmonics are usually less than a hundredth of the fundamental component, and hence can be ignored. The EMF of phase A can be defined as: where, E MMC and α are the amplitude and the phase angle of EMF. Equation (11) reflects the external characteristic of MMC, and expresses the relationship between the output voltage and EMF. Equation (12) reflects the internal characteristic of MMC, and describes the relationship between the common-mode voltage and circulating current.

Modelling of MMC Connected to the PMSG-Based WECS
It is assumed that in the steady state, the time delay of phase B and C are T/3 and 2T/3 respectively with phase A as the reference phase, where T is the fundamental period of the electrical quantities. It should be noted that the 2nd harmonic quantities are negative sequence.
Using the standard model, the voltage equations of PMSG can be expressed in the d-q frame as follows [21,22]: where, λ m represents the maximum flux produced by the permanent magnet; L d and L q are the d-axis and q-axis synchronous inductance respectively; R s is the stator resistance; u d and u q are the d-axis and q-axis output voltage of PMSG. In order to obtain the MMC mathematical model in the d-q frame, the Equations (11) and (12) should firstly be rewritten into three-phase form: where, T(θ) (θ = ω, 2ω) is the d-q transformation matrix, and is shown in Equation (23): From Equations (16), (17), (21) and (22), the d-q frame mathematical model of MMC connected to the PMSG-based WECS can be expressed as: Based on the mathematical model, the following equivalent circuit can be obtained. As shown in Figure 2, the upper two circuits are the d-axis equivalent circuits, which are built based on Equations (24) and (26). The lower two circuits are q-axis equivalent circuits, which are built based on Equations (25) and (27). In addition, Equations (24) and (25) describe the relationship between PMSG and MMC. Equations (26) and (27) represent the relationship between 2nd common-mode voltages and circulating currents in MMC.
It should be noted that the four equations are not independent. As is shown by the blue dashed line in Figure 2, the four equivalent circuits (or equations) are coupled with each other by the capacitors installed in MMC. This coupling relationship can be described as follows. The phase currents (i d , i q ) and circulating currents (i cir.d , i cir.q ) make up the ac components of arm currents. Then the arm currents will charge and discharge the capacitors, which causes the capacitor voltage ripple. The voltage ripple mainly contains fundamental frequency, 2nd order and 3rd order components, and will be delivered to the arm voltages. From Equations (13) and (14), the EMF and common-mode voltage can be calculated based on the arm voltages. Consequently, on the one hand, the odd-components of arm voltages will cause changes in the phase angle and amplitude of EMF, and then change in EMF, in return, influence the phase currents. On the other hand, the 2nd components of arm voltages make up the 2nd components of common-mode voltages, and become the voltage sources of 2nd circulating currents. The coupling relationship described in Figure 2 is the electrical part in the coupling relationship described in Figure  1. There are mainly two benefits to build the mathematical model in the d-q frame for MMC connected to the PMSG-based WECS:

•
Modelling MMC in the d-q frame facilitates the establishment of associations between the electrical quantities in PMSG and the electrical quantities in MMC, because the mathematical model of PMSG is built in the d-q frame, and the d-axis and q-axis synchronous inductances are not equal for the salient-pole PMSG.

•
In the d-q frame mathematical model, the time-varying quantities are transformed into constant quantities, which therefore can simplify the derivation. Due to this, the equivalent resistances in the MMC arms can be considered in the analysis, and the algebraic solution of non-linear equations can be obtained to calculate the amplitudes and phase angles of average switching functions.

A More Accurate Method to Calculate the Average Switching Function
IGBT gate signals, which control the on/off state of the IGBTs, directly determine the operation of MMC. In the mathematical model of MMC, the IGBT gate signals are described as the switching function, which can be expressed as follows: 1 ,when is on and is off ( ) 0 ,when is off and is on where, Si(t) denotes the switching function. S1 and S2 is the upper and lower IGBT in each sub-module respectively, which is shown in Figure 1.
The mathematical expression of the switching function in Equation (28) is a discrete function and difficult to be obtained, therefore the average switching function is normally adopted in the mathematical model of MMC [17]. Taking phase A as an example, the average switching functions of upper and lower arms can be expressed as: .
. a d d . From Equations (13) and (14), the EMF and common-mode voltage can be calculated based on the arm voltages. Consequently, on the one hand, the odd-components of arm voltages will cause changes in the phase angle and amplitude of EMF, and then change in EMF, in return, influence the phase currents. On the other hand, the 2nd components of arm voltages make up the 2nd components of common-mode voltages, and become the voltage sources of 2nd circulating currents. The coupling relationship described in Figure 2 is the electrical part in the coupling relationship described in Figure 1. There are mainly two benefits to build the mathematical model in the d-q frame for MMC connected to the PMSG-based WECS:

•
Modelling MMC in the d-q frame facilitates the establishment of associations between the electrical quantities in PMSG and the electrical quantities in MMC, because the mathematical model of PMSG is built in the d-q frame, and the d-axis and q-axis synchronous inductances are not equal for the salient-pole PMSG.

•
In the d-q frame mathematical model, the time-varying quantities are transformed into constant quantities, which therefore can simplify the derivation. Due to this, the equivalent resistances in the MMC arms can be considered in the analysis, and the algebraic solution of non-linear equations can be obtained to calculate the amplitudes and phase angles of average switching functions.

A More Accurate Method to Calculate the Average Switching Function
IGBT gate signals, which control the on/off state of the IGBTs, directly determine the operation of MMC. In the mathematical model of MMC, the IGBT gate signals are described as the switching function, which can be expressed as follows: 1, when S 1 is on and S 2 is off 0, when S 1 is off and S 2 is on (28) where, S i (t) denotes the switching function. S 1 and S 2 is the upper and lower IGBT in each sub-module respectively, which is shown in Figure 1.
The mathematical expression of the switching function in Equation (28) is a discrete function and difficult to be obtained, therefore the average switching function is normally adopted in the mathematical model of MMC [17]. Taking phase A as an example, the average switching functions of upper and lower arms can be expressed as: S dn.a (t) = 1 2 + S re f .a (t) − S add.a (t) (30) where, S up.a (t) and S dn.a (t) denote the average switching functions in the upper and lower arms respectively. S re f .a (t) is a fundamental harmonic sinusoidal function. S add.a (t) is a 2nd harmonic sinusoidal function.
In the traditional methods, the half of the modulation index is used as the amplitude of S re f .a (t) [16,17]. The modulation index is defined as Equation (31). The phase angle of S re f .a (t) is set to be the same as the phase angle of EMF. For the 2nd harmonic components S add.a (t), they are ignored in the traditional methods: For the ideal condition, the capacitor voltage only contains the dc component, and can be expressed as: Therefore, based on Equations (7), (8) and (13), the EMF can be calculated as: The calculation result of Equation (33) equals Equation (15), which indicates that the traditional method is accurate under the ideal conditions. In fact, in addition to the dc component, the capacitor voltage also contains fundamental, 2nd and 3rd order harmonics. These harmonics are also delivered to the arm voltages, and will cause the change in the amplitude and phase angle of EMF. In other words, the calculation result of Equation (33) is not equal to Equation (15) in practice. Therefore, the traditional methods are actually inaccurate.
Consequently, a more accurate method to calculate the value of the average switching function are proposed. Figure 3 shows the overall control scheme of MMC connected to the PMSG-based WECS. The reference signal is outputted by phase current control unit, which realizes the maximum power point tracking (MPPT). The additional signal is outputted by circulating current control unit, which suppresses the circulating current. The arm voltage signal is composed of the dc component, the reference signal and the additional signal. More details of the controller are referred in paper [12], which are out of the scope of this paper.
For the ideal condition, the capacitor voltage only contains the dc component, and can be expressed as: Therefore, based on Equations (7), (8) and (13), the EMF can be calculated as: The calculation result of Equation (33) equals Equation (15), which indicates that the traditional method is accurate under the ideal conditions. In fact, in addition to the dc component, the capacitor voltage also contains fundamental, 2nd and 3rd order harmonics. These harmonics are also delivered to the arm voltages, and will cause the change in the amplitude and phase angle of EMF. In other words, the calculation result of Equation (33) is not equal to Equation (15) in practice. Therefore, the traditional methods are actually inaccurate.
Consequently, a more accurate method to calculate the value of the average switching function are proposed. Figure 3 shows the overall control scheme of MMC connected to the PMSG-based WECS. The reference signal is outputted by phase current control unit, which realizes the maximum power point tracking (MPPT). The additional signal is outputted by circulating current control unit, which suppresses the circulating current. The arm voltage signal is composed of the dc component, the reference signal and the additional signal. More details of the controller are referred in paper [12], which are out of the scope of this paper. It can be seen in Figure 3 that the IGBT gate signal is obtained by inputting the arm voltage signal into the modulation controller. In the modulation controller, since the carrier frequency is high, it It can be seen in Figure 3 that the IGBT gate signal is obtained by inputting the arm voltage signal into the modulation controller. In the modulation controller, since the carrier frequency is high, it makes an ignorable change to the arm voltage signal. In a word, the arm voltage signal has almost the same amplitude and phase angle as the low frequency components (dc, fundamental and 2nd) in IGBT gate signals.
Therefore, in this paper, the arm voltage signal is chosen as the value of the average switching function in the mathematical model. In other word, the reference signal in the arm voltage signal is assumed to be the fundamental component S re f .a (t) in the average switching function; and the additional signal in the arm voltage signal is assumed to be the 2nd harmonic component S add.a (t) in the average switching function. The problem is that the amplitude and phase angle of the arm voltage signal are unknown. Thus, at first, the arm voltage signals are assumed to be sinusoidal functions containing unknown quantities, which are shown in Equations (34) and (35). Then, these unknown quantities will be obtained based on the coupling relationship of electrical quantities in PMSG and MMC.
Then, the average switching function can be rewritten as: where, M re f and β are the amplitude and the phase angle of the reference signal; M add and θ add are the amplitude and the phase angle of the additional signal.
In the average switching functions, although the 2nd harmonic components are smaller than the fundamental components, they are not neglected for the following two reasons: It should be noted that the four quantities (M re f , β, M add and θ add ) in in Equations (34) and (35) are unknown quantities. The solving method will be shown in the followings parts.

Calculation of the AC-Side Phase Current
For simplicity, the harmonics higher than 4th order are ignored in the derivation, because of their lower amplitudes compared with the fundamental components.
The ac-side phase currents in the d-q frame have the following relationship with the electromagnetic torque [5]: Equation (36) means that the electromagnetic torque can be controlled by controlling phase currents. Usually, the d-axis current is controlled to be zero to remove the coupling between the d-axis and q-axis currents in the electromagnetic torque [23,24]. Hence, the phase currents can be calculated as follows: Energies 2018, 11, 461 10 of 31 The mechanical power extracted from wind energy is determined by the wind speed and wind turbine, which is given by [25]: where ρ is the air mass density; A is the area covered by the wind blades; v wind is the wind speed; C p is the performance coefficient, which is determined by the pitch angle β pit and tip-speed ratio λ. Therefore, the electromagnetic torque can be calculated by: where, P 0 denotes the no-load loss of PMSG, which can be obtained from manufacturer's data. When wind energy system works under MPPT condition, the electromagnetic torque should be controlled as follows in order to maximize C p : where: r wind is the radius of wind blades; λ opt is the optimum tip-speed ratio, which is a constant value and can be obtained from the manufacturer's data; T 0 denotes the no-load torque.
The mechanical rotor speed can then be obtained as: Therefore, from Equations (36)-(42), the phase currents in the d-q frame can be derived as:

Calculation of the DC-Side Current
Equations (24)- (27) show the mathematical model of the MMC connected to PMSG, in which R s and R eq represent the losses of PMSG and MMC respectively. The mechanical power of PMSG and dc-side power of MMC have the following relationship: where: P m is the input mechanical power of PMSG; P l.PMSG denotes the stator loss of PMSG; P l.MMC denotes the loss of MMC; P dc is the dc-side power of MMC.
The dc-side current of MMC can then be derived based on Equations (24)- (27), (45) and (46): where, I s.rms and I cir.rms are the root-mean-square (RMS) value of the phase current and circulating current respectively.

Calculation of the Capacitor Voltage and Current
From Equations (3)-(6), (34) and (35), the capacitor currents in upper and lower arms can be calculated as: where: In the steady-state analysis, the converter is considered as working under stable operation condition and the capacitor voltages are balanced. Therefore, the capacitor voltages in the same arm can be simplified to be the same. Then, the capacitor voltage can be derived by the integration of the corresponding capacitor current, which can be expressed as: where: From Equations (54)-(59), in addition to the dc component, the capacitor voltages also contain the fundamental, 2nd and 3rd harmonic components, which are called as the capacitor voltage ripples. The voltage ripples should be limited to be under a certain value, because the excessive ripple voltage may result in the breakdown of IGBTs and reduce the stability of MMC.

Calculation of the Sub-Module Voltage and Arm voltage
The output voltages of sub-modules can be calculated by multiplying the switching function with the capacitor voltages, which can be derived based on Equations (54)-(59), (34) and (35): where: Equations (60)-(65) prove that the capacitor voltage ripples are delivered to the output voltages of sub-modules, which make the sub-module voltages contain not only the dc component but also the fundamental, 2nd, and 3rd harmonic components. Then, these harmonic components propagate to the arm voltages. In the steady state, the arm voltages can be simplified as:

Calculation of the EMS and Common-Mode Voltage
From Equation (14), the common-mode voltage is equal to the sum of upper and lower arm voltages, which can be expressed as: where: The common-mode voltage mainly contains the dc component and 2nd harmonic component, and the 2nd harmonic component will generate the 2nd order circulating current.
The EMF can be calculated by Equation (13): where:

Solution Procedure of the Unknowns in the Average Switching Functions
In Section 3, the formulas are derived for the electrical quantities in MMC. However, in these formulas, the amplitude and phase angle of the average switching functions are unknown quantities. Since there is a coupling relationship among the electrical quantities as described in Figure 2, the unknowns cannot be directly obtained. Therefore, in this part, four non-linear equations are established in the d-q frame in order to calculate the amplitude and phase angle of the average switching functions. Moreover, based on several simplifications, the algebraic solutions of the non-linear equations can be obtained, which can make the calculation results obtainable at one stroke and no need for iteration.
Equation (72) gives the function of the fundamental component in EMF. For simplicity, the items containing the product of the amplitude of reference signals, including M 2 add , M 2 re f and M add M re f , are ignored. These items have subtle influence on the results, but heavily add more complexity to the calculative process. In the simulation section, the results can prove the feasibility of this simplicity. Therefore, Equation (72) can be simplified as: The EMF in the d-q frame can be calculated by Equation (75): Then, Equation (76)  Equation (70) is the function of the 2nd component in common-mode voltage of phase A, and can be transformed into the d-q form by multiplying transform matrix, which is shown in Equation (77): Therefore, the following equation is obtained: Is M re f 8ωC cos(β + ϕ) In the steady state, the phase currents and circulating currents in the d-q frame are constant quantities. Therefore, their time derivatives equal to zero. Then, substituting Equations (76) and (78) into Equations (24)-(27), the following equations can be obtained: It can be seen that based on the aforementioned simplifications, the Equations (79) and (80) no longer contain M add and θ add , which greatly reduce the difficulty of solving equations. Then, the unknown quantities in the switching functions can be obtained by solving the non-linear Equations (79)-(82). The solving results can be expressed as: where: After obtaining the amplitudes and phase angles in the average switching functions, all the electrical quantities in MMC and PMSG can be calculated. It can be seen from Section 3 that, most of the electrical quantities are affected by the switching functions. Hence, compared with the traditional method, the adoption of a more accurate way to calculate the average switching functions will improve the accuracy of the analysis results.

Flow Chart of the Proposed Analysis Method
Based on the previous derivations, a steady-state analysis method is proposed for MMC connected to the PMSG-based WECS in this part. In the proposed method, only the wind speed (operating condition) is required as input, and all the electrical quantities in MMC, including the amplitudes, phase angles and their harmonics, can be calculated step by step, when the parameters, such as the capacitance and inductance in MMC and the flux linkage in PMSG, are set down. The flow chart of the proposed analysis method is summarized as shown in Figure 4.
After obtaining the amplitudes and phase angles in the average switching functions, all the electrical quantities in MMC and PMSG can be calculated. It can be seen from Section 3 that, most of the electrical quantities are affected by the switching functions. Hence, compared with the traditional method, the adoption of a more accurate way to calculate the average switching functions will improve the accuracy of the analysis results.

Flow Chart of the Proposed Analysis Method
Based on the previous derivations, a steady-state analysis method is proposed for MMC connected to the PMSG-based WECS in this part. In the proposed method, only the wind speed (operating condition) is required as input, and all the electrical quantities in MMC, including the amplitudes, phase angles and their harmonics, can be calculated step by step, when the parameters, such as the capacitance and inductance in MMC and the flux linkage in PMSG, are set down. The flow chart of the proposed analysis method is summarized as shown in Figure 4. The total calculation procedure can be divided into five steps:

•
Step 1: the value of the wind speed is inputted, which decides the operating condition of WCES.

•
Step 2: the mechanical power and the mechanical rotor speed of PMSG are calculated. It should be noted that the proposed method is also applicable for the condition, in which the wind turbine is not working under the MPPT control. In that condition, the first step is skipped, and the calculation procedure starts with the second step, which means the inputted operating condition is not the wind speed but the mechanical power and mechanical rotor.

•
Step 3: the angular speed, phase currents and dc-side currents of MMC are calculated. The The total calculation procedure can be divided into five steps:

•
Step 1: the value of the wind speed is inputted, which decides the operating condition of WCES.

•
Step 2: the mechanical power and the mechanical rotor speed of PMSG are calculated. It should be noted that the proposed method is also applicable for the condition, in which the wind turbine is not working under the MPPT control. In that condition, the first step is skipped, and the calculation procedure starts with the second step, which means the inputted operating condition is not the wind speed but the mechanical power and mechanical rotor.

•
Step 3: the angular speed, phase currents and dc-side currents of MMC are calculated. The calculations of these electrical quantities are not relevant to the switching functions.

•
Step 4: the unknown quantities in the average switching functions are calculated based on Equations (83)-(86).

•
Step 5: all the electrical quantities of MMC, including the amplitudes, the phase angles, and their harmonics, can be obtained based on the formulas shown in this paper.
The proposed method can be realized based on the proper calculation software, such as MATLAB. Since there is no iteration process, all the electrical quantities in MMC can be calculated at one stroke according to the proposed method, and hence the calculation speed is fast. In addition, for the wind conversion system, where the PMSG-side converter controls the DC voltage while the grid-side converter regulates active power [26], the proposed analysis method is also applicable.

Analysis of Capacitor Voltage Ripple in the MMC Connected to PMSG
Capacitor voltage ripple is one of the key problems in MMC. For the MMC applied in motor drives, this problem is more severe, because the capacitor voltage ripple becomes larger at low frequencies. This characteristic can also be seen in the capacitor voltage ripple formulas, shown in Equations (57)-(59). The large capacitor voltage ripple may lead to IGBT breakdown and reduce the system reliability. Therefore, for the MMC applied in motor drive, the circulating current injection methods are usually used to suppress the capacitor voltage ripple at start-up and low-frequency operations [25,27].
However, for MMC connected to the PMSG-based WECS, the characteristics are different. When MMC works under the MPPT control mode, according to Equations (39)-(42), the input power and torque can be expressed as: Then, the ac-side currents can be derived from Equations (44): where: When the power losses are ignored, the input power of PMSG approximately equals to the dc-side power of MMC. Then, according to Equations (45), (46) and (93), the dc-side currents can approximately be calculated as follows: where: In Equations (96) and (98), K 1 and K 2 are constants, which will not change with the different operation conditions. Equations (96) and (98) can be brought into Equations (57)-(59), and then the capacitor voltage ripple can be obtained as: Usually, circulating currents are suppressed to be zero or controlled to be a small value. Therefore, according to Equations (99)-(101), it can be seen that there is a positive correlation between the voltage ripple and system frequency. Based on the proposed analysis method, the fluctuation magnitude at different frequencies when the circulating currents are fully suppressed are calculated and shown in Figure 5. The parameters of PMSG and MMC are shown in Tables 1 and 2.
Usually, circulating currents are suppressed to be zero or controlled to be a small value. Therefore, according to Equations (99)-(101), it can be seen that there is a positive correlation between the voltage ripple and system frequency. Based on the proposed analysis method, the fluctuation magnitude at different frequencies when the circulating currents are fully suppressed are calculated and shown in Figure 5. The parameters of PMSG and MMC are shown in Tables 1 and 2.  In summary, for MMC connected to the PMSG-based WECS, there is a positive correlation between the capacitor voltage ripple and the system frequency. Therefore, for MMC connected to the PMSG-based WECS, capacitor voltage ripple suppression is not needed under low-frequency  In summary, for MMC connected to the PMSG-based WECS, there is a positive correlation between the capacitor voltage ripple and the system frequency. Therefore, for MMC connected to the PMSG-based WECS, capacitor voltage ripple suppression is not needed under low-frequency operation, which is different from the MMC applied in motor drives. In addition, the commonly used capacitor voltage ripple suppression method, in which the capacitor voltage ripple is suppressed by injecting high frequency circulating currents and common-mode voltages [25,27], is not applicable for MMC connected to the PMSG-based WECS. This is because that the large capacitor voltage ripple occurs when the wind speed approaches its rated value. The modulation ratio is high at that time, and hence there is only a small margin for injecting common-mode voltage, which leads to the need of very large injecting circulating currents and the over-modulation is highly possible to happen.

Capacitor Sizing For the MMC Connected to PMSG
The capacitors occupy approximately half the volume of MMC, and are the key components determining the performance of MMC.
The capacitor sizing mainly depends on the capacitor voltage ripples [28]. The capacitance should be large enough to limit the maximal voltage of capacitor below a given value; otherwise it may lead to IGBT breakdown and reduce the system reliability. However, choosing large capacitance will lead to high project cost and occupy a large area. Usually the appropriate capacitance should limit the capacitor voltage fluctuation ratio below 10%-15% [29]. Therefore, the calculation of the capacitor voltage is vital in the capacitor sizing.
According to Equations (54) and (55), the analytical expression of capacitor voltage can be expressed as: where, u cap.dc is the DC component and ∆u cap (t) is the fluctuation component. Then, the capacitor voltage fluctuation ratio ε can be defined as follows [29]: In the previous papers, the capacitor sizing method can be divided into two categories. One kind of the methods derives the capacitance requirement through the energy variation in the converter arms [30]. This kind of methods avoids deducing the equation of capacitor voltage ripple. Another kind of the methods sizes the capacitor by directly analyzing the capacitor voltage [29]. Since the mathematical expression of u cap (t) is complicated, the capacitor voltage fluctuation ratio ε is usually calculated by a lot of simplifications, which results in choosing an inappropriate value of capacitance. In addition, the previous capacitor sizing methods are proposed for the MMC used in HVDC transmission. The MMC connected to PMSG has different characteristics.
Since the proposed method can accurately calculate the capacitor voltages, it can be used for capacitor sizing. The paper focuses on the commonly used control strategy, in which the modulation controller uses the sinusoidal modulation method, and the 2nd harmonic circulating current is fully suppressed. There are optimization methods that can reduce the capacitance requirement. For example, in paper [31] an enhanced flat-topped modulation method is used to reduce the capacitance requirement. In paper [12] circulating current injection can also be used to suppress the capacitor voltage ripple. The influences of these optimization methods are out of the scope of this paper and will be studied in our future research.
From the analysis results in Section 5.1, for MMC connected to the PMSG-based WECS, the maximum fluctuation occurs when the wind speed reaches its rated value. Therefore, the capacitor should be sized under the rated wind speed.
The sizing procedure can be divided into two steps:

•
Step 1: Draw the parameter curve based on the proposed mathematical model and analysis method shown in Section 4. The parameter curve describes the relationship between the capacitance and the capacitor voltage fluctuation ratio ε. • Step 2: Select the appropriate value of capacitance according to the obtained parameter curve. Figure 6 shows the flow chart of drawing the parameter curve, where k denotes the loop number. Firstly, the initial value of capacitance C 1 , the step value of capacitance V step and the max value of loop number k max are inputted. Secondly, the loop number k is initialized to be zero. Then, k increases by 1. After that, the k-th value of the capacitance C k is calculated. Next, the capacitor voltage u cap.k (t), when the capacitance is C k , is calculated by using the proposed mathematical model and analysis method; the calculation step is shown in Figure 4. Then, the max value of u cap.k (t) can be obtained, and therefore the capacitor voltage fluctuation ratio ε k can be easily calculated by Equation (103). Finally, after executing the loop k times, the parameter curve between C k and ε k can be dawn.      Tables 1 and 2 (except the  capacitance). It can be seen in the figure that the voltage fluctuation decreases with the increasing of capacitance. Therefore, for the parameters shown in the Tables 1 and 2, if the voltage fluctuations are expected between 10-15%, the capacitance should be 4000-6000 µF; and if the voltage fluctuations are expected less than 10%, the capacitance should be greater than 6000 µF.
In this paper, the capacitance is chosen to be 5000 µF, and hence the capacitor voltage fluctuation will be lower than 12.41%. Figure 7 shows the parameter curve for the parameters shown in Tables 1 and 2 (except the capacitance). It can be seen in the figure that the voltage fluctuation decreases with the increasing of capacitance. Therefore, for the parameters shown in the Tables 1 and 2, if the voltage fluctuations are expected between 10-15%, the capacitance should be 4000-6000 μF; and if the voltage fluctuations are expected less than 10%, the capacitance should be greater than 6000 μF.

Simulation Model
In this section, simulations are carried out in MATLAB/Simulink (8.5.0.197613, MathWorks, Natick, MA, USA) in order to verify the accuracy of the proposed method and the correctness of the analysis. Figure 8 is the overall diagram of the simulation system. The simulation system is composed of three parts, which are the PSMG part, the MMC connected to PSMG part, and the AC grid part.
(a) In the PMSG part, the wind turbine model is built based on the wind power formula. The mechanical shaft model uses a two-mass rotor model, with separate masses for the turbine and generator. The PSMG model is the typical model provided by MATLAB/Simulink. More details of this part can be referred to the paper [32]. (b) In the MMC connected to PSMG part, a five level MMC is built. As is shown in Figure 7, the sub-modules use half-bridge topology. MPPT controller is used to make the PMSG obtain the highest possible power from wind. PI controllers are used to track the reference values in the inner current controller and the circulating current controller. The IGBT gate signals (the switching function in the mathematic model) are obtained by comparing the arm voltage signals with the carrier signals. Details of this part can be referred to the paper [12]. (c) In the AC grid part, the DC/AC converter is the traditional modular multilevel converter used in HVDC [17]. The AC grid is substituted with a three-phase AC source, in which the frequency is 50 Hz.
The simulation parameters of MMC and PMSG are listed in Tables 1 and 2. In this paper, the capacitance is chosen to be 5000 μF, and hence the capacitor voltage fluctuation will be lower than 12.41%.

Simulation Model
In this section, simulations are carried out in MATLAB/Simulink (8.5.0.197613, MathWorks, Natick, MA, USA) in order to verify the accuracy of the proposed method and the correctness of the analysis. Figure 8 is the overall diagram of the simulation system. The simulation system is composed of three parts, which are the PSMG part, the MMC connected to PSMG part, and the AC grid part.
(a) In the PMSG part, the wind turbine model is built based on the wind power formula. The mechanical shaft model uses a two-mass rotor model, with separate masses for the turbine and generator. The PSMG model is the typical model provided by MATLAB/Simulink. More details of this part can be referred to the paper [32]. (b) In the MMC connected to PSMG part, a five level MMC is built. As is shown in Figure 7, the submodules use half-bridge topology. MPPT controller is used to make the PMSG obtain the highest possible power from wind. PI controllers are used to track the reference values in the inner current controller and the circulating current controller. The IGBT gate signals (the switching function in the mathematic model) are obtained by comparing the arm voltage signals with the carrier signals. Details of this part can be referred to the paper [12]. (c) In the AC grid part, the DC/AC converter is the traditional modular multilevel converter used in HVDC [17]. The AC grid is substituted with a three-phase AC source, in which the frequency is 50Hz.
The simulation parameters of MMC and PMSG are listed in Tables 1 and 2.

The Proposed Method to Calculate the Average Switching Functions
As mentioned above, conventionally, the half of the modulation ratio and the phase angle of EMF are taken as the amplitude and phase angle of the fundamental component in the average switching functions. In this paper, a more accurate method is adopted to calculate the average

The Proposed Method to Calculate the Average Switching Functions
As mentioned above, conventionally, the half of the modulation ratio and the phase angle of EMF are taken as the amplitude and phase angle of the fundamental component in the average switching functions. In this paper, a more accurate method is adopted to calculate the average switching functions, in which the "reference signal" (shown in Figure 3) are considered as their fundamental components.
The simulation results for cases with different wind speed are shown in Table 3. In the table, the values of the amplitude and phase angle are obtained by processing the simulation data with the Fast Fourier Transform (FFT). In the first column, the simulation results show the amplitude and phase angle of the fundamental component in the IGBT gate signal, which can be considered as the benchmark. In the second column, the simulation results of the reference signals are shown. In the third column, the simulation results are the half of the modulation ratio and the phase angle of EMF. It can be seen that compare with the results in the third column, the results in the second column are more close to the results in the first column. This means that compared with the traditional method, the reference signal is more close to the fundamental component in the IGBT gate signal. Moreover, when the wind speed approaches its rated value, the errors between the half of the modulation ratio and the amplitude of the IGBT gate signal become larger.  Figure 9 further shows the error between the half of the modulation ratio and the amplitude of the fundamental component in the IGBT gate signal for cases with different wind speed. It can be seen that the error becomes larger as the wind speed increases. However, according to Table 4, the errors of the reference signals are small enough to be ignored under all the wind speeds operating conditions. nergies 2017, 10, x FOR PEER REVIEW 21 o ird column, the simulation results are the half of the modulation ratio and the phase angle of EM can be seen that compare with the results in the third column, the results in the second column a ore close to the results in the first column. This means that compared with the traditional metho e reference signal is more close to the fundamental component in the IGBT gate signal. Moreov hen the wind speed approaches its rated value, the errors between the half of the modulation ra nd the amplitude of the IGBT gate signal become larger.  Figure 9 further shows the error between the half of the modulation ratio and the amplitude e fundamental component in the IGBT gate signal for cases with different wind speed. It can een that the error becomes larger as the wind speed increases. However, according to Table 4, t rrors of the reference signals are small enough to be ignored under all the wind speeds operati onditions.    seen that the simulation and calculation results have similar amplitudes and phase angles, but there still exist small errors. Therefore, more accurate method to calculate the average switching functions can be used to improve the accuracy of steady-state analysis. The calculation results using both of the proposed analysis method and the new method of calculating the average switching functions are shown in Figures 11 and 12 in the next part. that the simulation and calculation results have similar amplitudes and phase angles, but there still exist small errors. Therefore, more accurate method to calculate the average switching functions can be used to improve the accuracy of steady-state analysis. The calculation results using both of the proposed analysis method and the new method of calculating the average switching functions are shown in Figures 11 and 12 in the next part. Figure 10. Comparisons between the simulation and calculation results using the traditional method to calculate the average switching function.

Simulation and Calculation Results for the Proposed Analysis Method
To verify the accuracy of the proposed steady-state analysis method, the simulation and calculation results for the cases with different wind speed are shown in Figure 11. In the subfigures, the red solid lines are the simulation results and the blue dashed lines are the calculation results. The circulating currents are suppressed to be zero. In each figure, from top to down, it presents the comparisons of the capacitor voltage ripples, the reference signals, the additional signals and the arm currents. To verify the accuracy of the proposed steady-state analysis method, the simulation and calculation results for the cases with different wind speed are shown in Figure 11. In the subfigures, the red solid lines are the simulation results and the blue dashed lines are the calculation results. The circulating currents are suppressed to be zero. In each figure, from top to down, it presents the comparisons of the capacitor voltage ripples, the reference signals, the additional signals and the arm currents. It can be seen that all the calculation results in Figure 11 fit well with the simulation results. Moreover, the simulation results in Figure 11c and Figure 10 are all carried out at the wind speed of 12 m/s. However, the calculation results in them, shown by the blue dashed lines, are obtained by using and without using the proposed method to calculate the average switching functions respectively. It can be seen that the calculation results obtained by using the proposed method, in Figure 11c, significantly have smaller errors. This means that by using more accurate way to calculate the average switching functions, the accuracy of the analysis method can be effectively improved.
To further verify the accuracy of the proposed steady-state analysis, the calculation results and their errors of more electrical quantities are shown in Table 4. The simulation and calculation are carried out at the rated wind speed.
It can be seen from Table 4 that the errors in the calculation results of both amplitudes and phase angles are small, and most of the errors, especially for the fundamental components, are less than 1%. This further proves that the proposed method has high accuracy.
In addition, in some papers such as [12], the circulating current is controlled to be a certain value (not zero) to suppress the capacitor voltage ripple. In order to prove that the proposed method also has high accuracy under this condition, Figure 12 presents the simulation and calculation results at the rated wind speed when cir I = 246.87 A and θ cir = 5.5292 rad [12]. It can be seen that all the calculation results in Figure 11 fit well with the simulation results. Moreover, the simulation results in Figure 11c and Figure 10 are all carried out at the wind speed of 12 m/s. However, the calculation results in them, shown by the blue dashed lines, are obtained by using and without using the proposed method to calculate the average switching functions respectively. It can be seen that the calculation results obtained by using the proposed method, in Figure 11c, significantly have smaller errors. This means that by using more accurate way to calculate the average switching functions, the accuracy of the analysis method can be effectively improved.
To further verify the accuracy of the proposed steady-state analysis, the calculation results and their errors of more electrical quantities are shown in Table 4. The simulation and calculation are carried out at the rated wind speed.
It can be seen from Table 4 that the errors in the calculation results of both amplitudes and phase angles are small, and most of the errors, especially for the fundamental components, are less than 1%. This further proves that the proposed method has high accuracy.
In addition, in some papers such as [12], the circulating current is controlled to be a certain value (not zero) to suppress the capacitor voltage ripple. In order to prove that the proposed method also has high accuracy under this condition, Figure 12 presents the simulation and calculation results at the rated wind speed when cir I = 246.87 A and θ cir = 5.5292 rad [12]. In Figure 12, since the circulating current is not zero, the arm current is not sinusoidal. It can be seen that the calculation results still fit well with the simulation results. This is because that the influence of circulating currents to other electrical quantities are considered in the stead-state analysis, and the additional signals, albeit small, are not ignored in the derivation procedure.
Moreover, the comparisons of calculation time between using the traditional method and using

Simulation and Calculation Results for the Proposed Analysis Method
To verify the accuracy of the proposed steady-state analysis method, the simulation and calculation results for the cases with different wind speed are shown in Figure 11. In the subfigures, the red solid lines are the simulation results and the blue dashed lines are the calculation results. The circulating currents are suppressed to be zero. In each figure, from top to down, it presents the comparisons of the capacitor voltage ripples, the reference signals, the additional signals and the arm currents.
It can be seen that all the calculation results in Figure 11 fit well with the simulation results. Moreover, the simulation results in Figures 11c and 10 are all carried out at the wind speed of 12 m/s. However, the calculation results in them, shown by the blue dashed lines, are obtained by using and without using the proposed method to calculate the average switching functions respectively. It can be seen that the calculation results obtained by using the proposed method, in Figure 11c, significantly have smaller errors. This means that by using more accurate way to calculate the average switching functions, the accuracy of the analysis method can be effectively improved.
To further verify the accuracy of the proposed steady-state analysis, the calculation results and their errors of more electrical quantities are shown in Table 4. The simulation and calculation are carried out at the rated wind speed.
It can be seen from Table 4 that the errors in the calculation results of both amplitudes and phase angles are small, and most of the errors, especially for the fundamental components, are less than 1%. This further proves that the proposed method has high accuracy.
In addition, in some papers such as [12], the circulating current is controlled to be a certain value (not zero) to suppress the capacitor voltage ripple. In order to prove that the proposed method also has high accuracy under this condition, Figure 12 presents the simulation and calculation results at the rated wind speed when I cir = 246.87 A and θ cir = 5.5292 rad [12].
In Figure 12, since the circulating current is not zero, the arm current is not sinusoidal. It can be seen that the calculation results still fit well with the simulation results. This is because that the influence of circulating currents to other electrical quantities are considered in the stead-state analysis, and the additional signals, albeit small, are not ignored in the derivation procedure.
Moreover, the comparisons of calculation time between using the traditional method and using the proposed method to calculate the average switching function are shown in Table 5. The calculations are carried out in MATLAB software. The comparison results show that the calculation speed based on the proposed method is also very fast. This mainly benefits from that the unknown quantities in the average switch functions are obtainable at one stroke and no iteration is needed.

Simulation Results of the Capacitor Voltage Ripple
In Section 5, the characteristics of capacitor voltage ripple in the MMC connected to PMSG is analyzed. Different from the MMC applied in motor drives, for the MMC connected to PMSG, there is a positive correlation between the capacitor voltage ripple and the system frequency. In Figure 13, waveforms are shown when the wind speed increases from 6 m/s at 2 s to12 m/s at 6 s. Figure 13a,b present the waveforms of the input power, electromagnetic torque, wind speed and system frequency. In Figure 13c, the circulating current is fully suppressed. The capacitor voltage ripple becomes larger as the wind speed increases, and the maximum fluctuation occurs when the wind speed reaches its rated value at 6 s. The simulation results coincide with the analysis in Section 5, which verifies the correctness of the analysis to capacitor voltage ripples. Moreover, the simulation proves that, for MMC connected to the PMSG-based WECS, capacitor voltage ripple suppression is not needed under low-frequency operation. In Section 5, the characteristics of capacitor voltage ripple in the MMC connected to PMSG is analyzed. Different from the MMC applied in motor drives, for the MMC connected to PMSG, there is a positive correlation between the capacitor voltage ripple and the system frequency. In Figure 13, waveforms are shown when the wind speed increases from 6 m/s at 2 s to12 m/s at 6 s.  Figure 13a,b present the waveforms of the input power, electromagnetic torque, wind speed and system frequency. In Figure 13c, the circulating current is fully suppressed. The capacitor voltage ripple becomes larger as the wind speed increases, and the maximum fluctuation occurs when the wind speed reaches its rated value at 6 s. The simulation results coincide with the analysis in Section 5, which verifies the correctness of the analysis to capacitor voltage ripples. Moreover, the simulation

Simulation Results of the Capacitor Voltage Ripple under Different Capacitance Conditions
In order to verify the correctness of the proposed capacitor sizing method. Comparisons between the simulation and calculation results of the capacitor voltage ripple under the different capacitance are carried out in this section. Figure 14 shows the waveforms of capacitor voltages. The parameters are shown in Tables 1 and 2 except that the capacitances from Figure 14a-d, are 4000 µF, 5000 µF, 6000 µF and 7000 µF respectively. Since the maximum fluctuation of the capacitor voltage occurs at the rated wind speed, all the simulations in Figure 14 are carried out at the rated wind speed. Table 6 presents the comparisons. The voltage fluctuation value ε is obtained by using Equation (103), and the calculation results are obtained by parameter curve shown in Figure 9.
It can be seen from Table 6 that the calculation results meet well with the simulation results, which prove the accuracy of the proposed method. Therefore, for the parameters shown in Tables 1 and 2, if the voltage fluctuations are expected between 10-15%, the capacitance should be 4000-6000 µF; and if the voltage fluctuations are expected less than 10%, the capacitance should be greater than 6000 µF. are carried out in this section. Figure 14 shows the waveforms of capacitor voltages. The parameters are shown in Table 1 and 2 except that the capacitances from Figure 14a-d, are 4000 μF, 5000 μF, 6000 μF and 7000 μF respectively. Since the maximum fluctuation of the capacitor voltage occurs at the rated wind speed, all the simulations in Figure 14 are carried out at the rated wind speed. Table 6 presents the comparisons. The voltage fluctuation value ε is obtained by using Equation (103), and the calculation results are obtained by parameter curve shown in Figure 9.

Experiment Validation
In order to further prove that the proposed method has high accuracy, a down-scaled prototype is implemented. The photograph of the prototype is presented in Figure 15a, and its diagram is shown in Figure 15b.

Experiment Validation
In order to further prove that the proposed method has high accuracy, a down-scaled prototype is implemented. The photograph of the prototype is presented in Figure 15a, and its diagram is shown in Figure 15b.  The experiment part intends to verify that the calculation results in the MMC are accurate, and only focuses on the steady state. Since the composition of torque in PMSG is zero in the steady state, the PMSG can be taken as a voltage source for the MMC. Therefore, the PMSG part is emulated by the programmable power source, in which the wind turbine characteristics are programmed to emulate a generator operating at variable speed and variable voltage [33]. The experiment parameters of the emulated PMSG is shown in Table 2. In the MMC part, each arm contains four sub-modules. The control scheme of MMC is the same as it in the simulation part, which can be referred in paper [12]. For the hardware realization of the control scheme, the modulation parts are realized in Field Programmable Gate Array (FPGA) chip, called EP4CE30F23 (Altera Corporation, San Jose, CA, USA); the inner-current controller and circulating current controller are realized by Digital Signal Processor (DSP) chip, called TMS320f28335 (Texas Instruments, Dallas, TX, USA). The experiment parameters of MMC are shown in Table 1.
In Figures 16-18, the comparisons between the experiment waveforms and the calculation results are shown. The blue lines are the waveforms of benchmark angles, used to observe the phase angle of acquired signals. The red lines are the waveforms of reference and additional signals. The benchmark angles, reference signals and additional signals are first acquired by FPGA, and then they are outputted to oscilloscope via digital-analogue (D/A) converter. The orange and green lines are the waveforms of capacitor voltages and arm currents. The black dashed lines are the calculation results obtained by using the proposed method.
the inner-current controller and circulating current controller are realized by Digital Signal Processor (DSP) chip, called TMS320f28335 (Texas Instruments, Dallas, TX, USA). The experiment parameters of MMC are shown in Table 1.
In Figures 16-18, the comparisons between the experiment waveforms and the calculation results are shown. The blue lines are the waveforms of benchmark angles, used to observe the phase angle of acquired signals. The red lines are the waveforms of reference and additional signals. The benchmark angles, reference signals and additional signals are first acquired by FPGA, and then they are outputted to oscilloscope via digital-analogue (D/A) converter. The orange and green lines are the waveforms of capacitor voltages and arm currents. The black dashed lines are the calculation results obtained by using the proposed method.  the inner-current controller and circulating current controller are realized by Digital Signal Processor (DSP) chip, called TMS320f28335 (Texas Instruments, Dallas, TX, USA). The experiment parameters of MMC are shown in Table 1.
In Figures 16-18, the comparisons between the experiment waveforms and the calculation results are shown. The blue lines are the waveforms of benchmark angles, used to observe the phase angle of acquired signals. The red lines are the waveforms of reference and additional signals. The benchmark angles, reference signals and additional signals are first acquired by FPGA, and then they are outputted to oscilloscope via digital-analogue (D/A) converter. The orange and green lines are the waveforms of capacitor voltages and arm currents. The black dashed lines are the calculation results obtained by using the proposed method.   Figure 16 shows the comparisons when the prototype works under the rated wind speed. The system frequency is 40 Hz. In Figures 17 and 18, the prototype works at the wind speeds of 9 m/s and 6m/s respectively, in which the system frequencies are 30 Hz and 20 Hz respectively. From Figures  16-18, it can be seen that the calculation results fit well with the experiment waveforms, which further validate the accuracy of the proposed analysis method.  Figure 16 shows the comparisons when the prototype works under the rated wind speed. The system frequency is 40 Hz. In Figures 17 and 18, the prototype works at the wind speeds of 9 m/s and 6m/s respectively, in which the system frequencies are 30 Hz and 20 Hz respectively. From Figures 16-18, it can be seen that the calculation results fit well with the experiment waveforms, which further validate the accuracy of the proposed analysis method.

Conclusions
This paper proposes a steady-state analysis method for MMC connected to the PMSG-based WECS. All the electrical quantities in MMC, including the amplitudes, phase angles and their harmonics, can be easily obtained according to the proposed method step by step. The proposed method can help design the parameters and optimize the operation of MMC used in the PMSG-based WECS.
The analysis method is based on the proposed d-q frame mathematical model of the MMC connected to PMSG. By using the d-q frame, the time-varying quantities are transformed into constant quantities, which simplifies the derivation. Due to this, the equivalent resistances in the MMC arms can be considered, and the algebraic solution of non-linear equations can be obtained to calculate the unknowns in the average switching functions.
A more accurate method to calculate the average switching functions are used in this paper. The comparisons between the simulation and calculation results show that the accuracy of analysis can be effectively improved by using the proposed method.
The analysis to capacitor voltage ripple concludes that there is a positive correlation between the capacitor voltage ripple and the system frequency for MMC connected to the PMSG-based WECS. Therefore, the capacitor voltage ripple suppression is not needed under low-frequency operation, which is different from the MMC applied in motor drives. In addition, a capacitor sizing method is also proposed in this paper, and is verified by the simulation results.
Comparisons between the simulation and calculation results have shown that the proposed steady-state analysis method has high accuracy, and most of the calculation errors, especially for the fundamental components, are less than 1%. Experiments carried out by a down-scaled prototype further prove that the proposed method has high accuracy under different operating condition.