# Behavior of a Dual Stator Induction Machine Fed by Neutral Point Clamped Multilevel Inverter

Authors: Mohamed Arbi Khlifi
Journal: Journal of Energy
 (2018)

---

## Abstract

Use of multilevel inverters with induction machines has become popular in most energy conversion and management systems. The present paper discusses the behavior of the dual stator-winding induction machine (DSIM) for the power source, which is a power supply with Neutral Point Clamped (NPC). Multilevel inverter with PWM technique control is analyzed. The DSIM control is obtained by the PWM technique to multicarrier PWM technique and after a comparative study of various characteristics of the DSIM taking into account the different electrical offsets between the two stars (0°, 30°, 60°). The gap between the two stars was considered, and it is impacted on torque and rates harmonic distortion of circulation currents.


---

## Body

## 1. Introduction

Induction machines with n-phase symmetrical have an equal angle between the successive phases of the stator a=2p/n. It is assumed that the windings are distributed sinusoidal, so that all higher spatial harmonics of the magnetomotive force can be neglected. The machines with odd number of phase have the distinction of eliminating the first harmonic of the larger current, harmonic “5” with pentaphase machine, and the harmonic “7” with heptaphase machine [1–4].One of the main objectives of the multiphase machines and the reliability are known; for example, in the case of multistar machinery for failure in a phase we just need to disconnect the star that contains the faulty phase, which reduces torque, deployed 50% with a significant increase in torque ripples and magnetic Song [5–8]. Through the history of multiphase machines, the most popular multiphase machine type has been a dual three-phase machine. Dual three-phase machines are characterized by the multiphase structure with two sets of three-phase stator windings in the same stator frame. As a disadvantage, dual three-phase machines can suffer from problems with undesired stator current harmonics [9–11].When power is increased, problems arise both in the voltage inverter and in the machine itself. The static switches of the voltage inverter must switch high currents and it is often necessary to place several structures in parallel, given power, reducing switching currents through the voltage increase. The inverter voltage PWM technique imposes high voltage gradients, thus causing accelerated aging of insulating [4, 12–14].To avoid this, one among the possible solutions is the large number of phases leading to distribute the power supplied to the machine. Another solution is to conduct multilevel inverters providing a supply of better quality while requiring smaller sizes switches (IGBT). This allows for simple structures and tested for inverters, with reduced thermal problems and electromagnetic interference [3, 15, 16].The idea of a machine with several windings is born in order to have two separate windings. Compared with the conventional machine, having a second winding facilitates control of the speed of the machine and transmits significant power through a static converter.Multiphase machines are machines whose number of phases q is greater than three and can provide more compared to three-phase machines because they cannot obey the requirements requested in specific application areas; they are developed mainly in the field of training for variable speed high power because the multiplication of the number of phases allows on the one hand improving the safety of the operation until at least three phases are active; there may be up to three phases open without the neutral connection which is required on the other hand to reduce the size of components in power modulators of energy.These constraints make the multiphase machine a very interesting concept; it is present in the fields of marine, rail traction, the petrochemical industry, avionics, automotive, etc. [17–20].The best known multiphase machine is probably the dual stator-winding induction machine whose two stars are identical and share the same stator and are phase shifted by an angle of 30°. Its windings have the same number of poles and are supplied at the same frequency. The structure of the rotor is identical to that of a three-phase machine, so it can be either squirrel cage or wound to form a three-phase winding.


## 2. Mathematical Model of Dual Stator Induction Motor

A dual three-phase machine is the most common multiphase machine structure. Dual three-phase machines have two sets of three-phase stator windings in the same stator frame. Displacement between the winding sets can take different values. However, only 0, 30, or 60 electrical degrees are actually encountered in practice. If the winding sets are not spatially shifted, the resulting machine is essentially a conventional three-phase machine with two parallel winding sets. On the other hand, displacing the winding set by 60 electrical degrees results in a symmetrical six-phase machine. Although both of these alternatives can offer some advantages, the most popular solution is that the star which connected three-phase stator windings is spatially shifted by 30 electrical degrees, and the neutral points of the sets are galvanic ally isolated from each other.(i)
The phases of the first star aresa1, sb1, and sc1, phases which take the second star are sa2,sb2, and sc2, and rotor phases are marked by ra,rb, and rc.(ii)
α expresses the angle of displacement between the two stars.(iii)
θ1 expresses the position of the rotor (phase ra) compared to a star (phase sa1).(iv)
θ2 expresses the position of the rotor (phase ra) compared to a star (phase sa2).(v)
The operation of the machine is assumed without saturation of magnetic circuit and neglecting the effect of hysteresis and Foucault currents.(vi)
The construction of the machine is assumed to be homogeneous.(vii)
The magnetomotive force created by each phase of the two frames is sinusoidal distribution.(viii)
The two-phase stator windings are balanced and identical.Taking into account the simplifying assumptions mentioned above and notation of vectors of quantities of voltage, current, and flux, we can write the following vectors:For star 1,(1)Vs1=Vas1Vbs1Vcs1TIs1=Ias1Ibs1Ics1Tϕs1=ϕas1ϕbs1ϕcs1TFor star 2,(2)Vs2=Vas2Vbs2Vcs2TIs2=Ias2Ibs2Ics2Tϕs2=ϕas2ϕbs2ϕcs2TFor rotor,(3)Vr=VarVbrVcrTIr=IarIbrIcrTϕr=ϕarϕbrϕcrTThe electric equations of star 1 and star 2 and of rotor are, respectively, expressed by(4)Vs1=Rs1.Is1+ddtϕs1(5)Vs2=Rs2.Is2+ddtϕs2(6)Vr=Rr.Ir+ddtϕrExpressions flux is as follows:(7)ϕs1=Ls1,s1Is1+Ms1,s2Is2+Ms1,rIr(8)ϕs2=Ms2,s1Is1+Ls2,s2Is2+Ms2,rIr(9)ϕr=Mr,s1Is1+Mr,s2Is2+Lr,rIrwhere Rs1,Rs2,Rr are the stator resistance matrix (star 1 and 2) and rotor andLs1,Ls2,Lr are the leakages inductances of one phase of star 1, star 2, and the rotor.(10)Ms2,s1=Ms1,s2T=Mr,s1=Ms1,rT=Mr,s2=Ms2,rTDuring the application of the d-q transformation and making the necessary manipulations, (4), (5), and (6) in d-q become(11)Vds1=Rs1.ids1+dϕds1dt-waϕqs1(12)Vqs1=Rs1.iqs1+dϕqs1dt+waϕds1(13)Vds2=Rs2.ids2+dϕds2dt-waϕqs2(14)Vqs2=Rs2.iqs2+dϕqs2dt+waϕds2(15)0=Rr.idr+dϕdrdt-wa-wϕqr(16)0=Rr.iqr+dϕqrdt+wa-wϕdrwith w=d(ψ-θm)/dt being speed of rotation of the coordinate (d, q) relative to the rotor and wa=dψ/dt being speed of rotation of the coordinate (d, q) relative to star 1.Applying the Park transformation on the equations of flux, we obtain(17)ϕds1=Ls1ids1+32Lmsids1+32Lmsids2+32Msridrϕqs1=Ls1iqs1+32Lmsiqs1+32Lmsiqs2+32Msriqrϕds2=Ls2ids2+32Lmsids2+32Lmsids1+32Msridrϕqs2=Ls2iqs2+32Lmsiqs2+32Lmsiqs1+32Msriqrϕdr=Lridr+32Lmridr+32Msrids1+32Msrids2ϕqr=Lriqr+32Lmriqr+32Msriqs1+32Msriqs2(18)32Lms=32Lmr=32Msr=Lmwith Lm being the cyclic mutual inductance between star 1, star 2, and the rotor.The previous system of equations becomes(19)ϕds1=Ls1.ids1+Lmids1+ids2+idr(20)ϕqs1=Ls1.iqs1+Lmiqs1+iqs2+iqr(21)ϕds2=Ls2.ids2+Lmids1+ids2+idr(22)ϕqs2=Ls2.iqs2+Lmiqs1+iqs2+iqr(23)ϕdr=Lr.idr+Lmids1+ids2+idr(24)ϕqr=Lr.iqr+Lmiqs1+iqs2+iqrwithLs1+Lm being the cyclic inductance of the star 1; Ls2+Lm being the cyclic inductance of the star 2; Lr+Lm being the cyclic inductance of the rotor.The electromagnetic torque can be expressed by(25)Cem=pϕds1iqs1-ϕqs1ids1+ϕds2iqs2-ϕqs2ids2By introducing the equations of flux (19) and (20) in (25) we obtain(26)Cem=pLmiqs1+iqs2idr-ids1+ids2iqr)Introducing the equations of current idr and iqr in (25) we obtain(27)Cem=pLmLm+Lriqs1+iqs2ϕdr-ids1+ids2ϕqrFor a reference frame related to the stator (wa=0) we can obtain the state equation of the form(28)X˙=AX+BUwith X=ϕds1ϕqs1ϕds2ϕqs2ϕdrϕqr being state vector and U=Vds1Vqs1Vds2Vqs2 being input vector.After a calculation, we obtain the following matrices:(29)A=La-Ls1Ts1Ls10LaTs1Ls20LaTs1Lr00La-Ls1Ts1Ls10LaTs2Ls10LaTs1LrLaTs2Ls10La-Ls2Ts2Ls20LaTs2Lr00LaTs2Ls10La-Ls2Ts2Ls20LaTs2LrLaTrLs10LaTrLs20La-LrTrLr-w0LaTrLs10LaTrLs2wLa-LrTrLr(30)B=100001000010000100000000(31)With:La=11/Ls1+1/Ls2+1/Lr+1/Lm(32)Ts1=Ls1Rs1;Ts2=Ls2Rs2;Tr=LrRrThe dual stator-winding induction machine model, the NPC inverter model, and the controller were implemented in MATLAB. The reference reactive power of the outer stator was set to zero to ensure unity power factor of the grid.


## 3. Two Three-Level NPC VSI-DSIM Cascades

### 3.1. Two Three-Level NPC VSI Structures

The operating mode is to supply a DSIM by an inverter dual stator-winding, that is to say, six arms controlled by a PWM control, as shown in Figure1.Figure 1
Structural diagram of an inverter double star.In this case the first three arms (A, B, and C) supply the first star which are shifted 120° between them, the same for the other three arms (D, E, and F), but the two stars are shifted by an electrical angleα. Figure 1 shows the structural diagram of an inverter double star that supplies a double star load (standard R-L) [6].The matrix below shows the simple voltages for the two stars with shifted angel (0°):(33)
V
A
N
1
t
U
B
N
1
t
U
C
N
1
t
U
D
N
2
t
U
E
N
2
t
U
F
N
2
t
=
1
3
2
-
1
-
1
0
0
0
-
1
2
-
1
0
0
0
-
1
-
1
2
0
0
0
0
0
0
2
-
1
-
1
0
0
0
0
0
0
-
1
-
1
2
-
1
-
1
2
V
A
O
t
V
B
O
t
V
C
O
t
V
D
O
t
V
E
O
t
V
F
O
t
=
E
3
2
-
1
-
1
0
0
0
-
1
2
-
1
0
0
0
-
1
-
1
2
0
0
0
0
0
0
2
-
1
-
1
0
0
0
0
0
0
-
1
-
1
2
-
1
-
1
2
C
1,1
C
2,1
C
3,1
C
4,1
C
5,1
C
6,1Writing the following matrix represents the composed voltages for the two stars:(34)UABtUBCtUCAtUDEtUEFtUFDt=1-1000001-1000-1010000001-100000000-110-11VAOtVBOtVCOtVDOtVEOtVFOt=E1-1000001-1000-1010000001-100000000-110-11C1,1C2,1C3,1C4,1C5,1C6,1Voltages between the two nodesN1 and N2 are(35)VON1t=-13(VAOt+VBOt+VCOt(36)VON2t=-13(VDOt+VEOt+VFOt(37)VN1N2=13VAO+VBO+VCO-VDO-VEO-VFOThese voltages are expressed in terms of control signals and the DC voltage where(38)VN1N2=13E2C1,1+2C2,1+2C2,1-2C4,1-2C5,1-2C6,1


### 3.2. PWM Strategy for the Two Three-Level NPC VSI

We adopted this technique (multicarrier PWM technique) to control the multilevel inverter. This technique is based on conventional sinusoidal modulation. It is the comparison of a sinusoidal signal called modulating wave over a triangular signal called carrier.The PWM technique is to treat separately the case where the set point is positive and those where it is negative, as shown in Figure2.Figure 2
Principle of PWM technique (multicarrier PWM technique) for a seven-level inverter.For an inverter with m levels, the command requires(m-1) triangular signals at the same frequency fp and the same magnitude Ap (peak to peak).The control signals of switches for the NPC inverter are obtained from the three reference sinusoidal signals phase shifted by2π/3, frequency fm, and magnitude Am.PWM technique is characterized by two reports, the rate of modulation and the frequency ratio, defined as follows:ma=Am/(m-1)Ap:the rate of modulation.mc=fp/fm:the frequency ratio.


### 3.3. Simulation of the DSIM Supplied by a Multilevel Inverter

#### 3.3.1. Principle

To identify the simulation results, we will look at the evolution of electromagnetic, mechanical, electrical, and magnetic characteristics for the various levels of the inverter used, taking into account a shiftα = 30°.The simulation of the inverter at m levels has been established by the software “Matlab Simulink”.


#### 3.3.2. Results of Simulation

Figures3–6 illustrate the evolution of the output voltage of an arm of the inverter when increasing the number of levels.Figure 3
Waveform of output voltageVas1 for m=2.Figure 4
Waveform of output voltageVas1 for m=3.Figure 5
Waveform of output voltageVas1 for m=5.Figure 6
Waveform of output voltageVas1 for m=7.According to the results found there, when the level of the inverter voltage is m = 2 to m = 7, the output voltage approaches more and more perfect sinusoidal form.Figures7–10 illustrate the evolution of the electromagnetic torque and rotational speed for various levels of the inverter used.Figure 7
Electromagnetic torque and rotational speed for m = 2.Figure 8
Electromagnetic torque and rotational speed for m = 3.Figure 9
Electromagnetic torque and rotational speed for m = 5.Figure 10
Electromagnetic torque and rotational speed for m = 7.Table1 shows the rate of undulation on the following number of levels corresponding to.Table 1
relative torque ripple rate as function of the inverter levels number.
Levels number Relative torque ripple rate in % m = 2 45.2% m = 3 10.52% m = 5 5.1% m = 7 4.96%The influence of the waveform of the supply voltage appears at the level of electromagnetic torque in steady state and has no influence on the shape of the characteristic speed of rotation.


## 3.1. Two Three-Level NPC VSI Structures

The operating mode is to supply a DSIM by an inverter dual stator-winding, that is to say, six arms controlled by a PWM control, as shown in Figure1.Figure 1
Structural diagram of an inverter double star.In this case the first three arms (A, B, and C) supply the first star which are shifted 120° between them, the same for the other three arms (D, E, and F), but the two stars are shifted by an electrical angleα. Figure 1 shows the structural diagram of an inverter double star that supplies a double star load (standard R-L) [6].The matrix below shows the simple voltages for the two stars with shifted angel (0°):(33)
V
A
N
1
t
U
B
N
1
t
U
C
N
1
t
U
D
N
2
t
U
E
N
2
t
U
F
N
2
t
=
1
3
2
-
1
-
1
0
0
0
-
1
2
-
1
0
0
0
-
1
-
1
2
0
0
0
0
0
0
2
-
1
-
1
0
0
0
0
0
0
-
1
-
1
2
-
1
-
1
2
V
A
O
t
V
B
O
t
V
C
O
t
V
D
O
t
V
E
O
t
V
F
O
t
=
E
3
2
-
1
-
1
0
0
0
-
1
2
-
1
0
0
0
-
1
-
1
2
0
0
0
0
0
0
2
-
1
-
1
0
0
0
0
0
0
-
1
-
1
2
-
1
-
1
2
C
1,1
C
2,1
C
3,1
C
4,1
C
5,1
C
6,1Writing the following matrix represents the composed voltages for the two stars:(34)UABtUBCtUCAtUDEtUEFtUFDt=1-1000001-1000-1010000001-100000000-110-11VAOtVBOtVCOtVDOtVEOtVFOt=E1-1000001-1000-1010000001-100000000-110-11C1,1C2,1C3,1C4,1C5,1C6,1Voltages between the two nodesN1 and N2 are(35)VON1t=-13(VAOt+VBOt+VCOt(36)VON2t=-13(VDOt+VEOt+VFOt(37)VN1N2=13VAO+VBO+VCO-VDO-VEO-VFOThese voltages are expressed in terms of control signals and the DC voltage where(38)VN1N2=13E2C1,1+2C2,1+2C2,1-2C4,1-2C5,1-2C6,1


## 3.2. PWM Strategy for the Two Three-Level NPC VSI

We adopted this technique (multicarrier PWM technique) to control the multilevel inverter. This technique is based on conventional sinusoidal modulation. It is the comparison of a sinusoidal signal called modulating wave over a triangular signal called carrier.The PWM technique is to treat separately the case where the set point is positive and those where it is negative, as shown in Figure2.Figure 2
Principle of PWM technique (multicarrier PWM technique) for a seven-level inverter.For an inverter with m levels, the command requires(m-1) triangular signals at the same frequency fp and the same magnitude Ap (peak to peak).The control signals of switches for the NPC inverter are obtained from the three reference sinusoidal signals phase shifted by2π/3, frequency fm, and magnitude Am.PWM technique is characterized by two reports, the rate of modulation and the frequency ratio, defined as follows:ma=Am/(m-1)Ap:the rate of modulation.mc=fp/fm:the frequency ratio.


## 3.3. Simulation of the DSIM Supplied by a Multilevel Inverter

### 3.3.1. Principle

To identify the simulation results, we will look at the evolution of electromagnetic, mechanical, electrical, and magnetic characteristics for the various levels of the inverter used, taking into account a shiftα = 30°.The simulation of the inverter at m levels has been established by the software “Matlab Simulink”.


### 3.3.2. Results of Simulation

Figures3–6 illustrate the evolution of the output voltage of an arm of the inverter when increasing the number of levels.Figure 3
Waveform of output voltageVas1 for m=2.Figure 4
Waveform of output voltageVas1 for m=3.Figure 5
Waveform of output voltageVas1 for m=5.Figure 6
Waveform of output voltageVas1 for m=7.According to the results found there, when the level of the inverter voltage is m = 2 to m = 7, the output voltage approaches more and more perfect sinusoidal form.Figures7–10 illustrate the evolution of the electromagnetic torque and rotational speed for various levels of the inverter used.Figure 7
Electromagnetic torque and rotational speed for m = 2.Figure 8
Electromagnetic torque and rotational speed for m = 3.Figure 9
Electromagnetic torque and rotational speed for m = 5.Figure 10
Electromagnetic torque and rotational speed for m = 7.Table1 shows the rate of undulation on the following number of levels corresponding to.Table 1
relative torque ripple rate as function of the inverter levels number.
Levels number Relative torque ripple rate in % m = 2 45.2% m = 3 10.52% m = 5 5.1% m = 7 4.96%The influence of the waveform of the supply voltage appears at the level of electromagnetic torque in steady state and has no influence on the shape of the characteristic speed of rotation.


## 3.3.1. Principle

To identify the simulation results, we will look at the evolution of electromagnetic, mechanical, electrical, and magnetic characteristics for the various levels of the inverter used, taking into account a shiftα = 30°.The simulation of the inverter at m levels has been established by the software “Matlab Simulink”.


## 3.3.2. Results of Simulation

Figures3–6 illustrate the evolution of the output voltage of an arm of the inverter when increasing the number of levels.Figure 3
Waveform of output voltageVas1 for m=2.Figure 4
Waveform of output voltageVas1 for m=3.Figure 5
Waveform of output voltageVas1 for m=5.Figure 6
Waveform of output voltageVas1 for m=7.According to the results found there, when the level of the inverter voltage is m = 2 to m = 7, the output voltage approaches more and more perfect sinusoidal form.Figures7–10 illustrate the evolution of the electromagnetic torque and rotational speed for various levels of the inverter used.Figure 7
Electromagnetic torque and rotational speed for m = 2.Figure 8
Electromagnetic torque and rotational speed for m = 3.Figure 9
Electromagnetic torque and rotational speed for m = 5.Figure 10
Electromagnetic torque and rotational speed for m = 7.Table1 shows the rate of undulation on the following number of levels corresponding to.Table 1
relative torque ripple rate as function of the inverter levels number.
Levels number Relative torque ripple rate in % m = 2 45.2% m = 3 10.52% m = 5 5.1% m = 7 4.96%The influence of the waveform of the supply voltage appears at the level of electromagnetic torque in steady state and has no influence on the shape of the characteristic speed of rotation.


## 4. Comparative Study of Various Characteristics of the DSIM Taking into Account the Different Electrical Offsets between the Two Stars (0°, 30°, 60°)

Figures11, 12, and 13, respectively, represent the characteristic electromagnetic torque, rotational speed, and torque-speed characteristic of the DSIM supplied by a multilevel inverter (by 7 levels) for α=0°,30°,60° (see Tables 2 and 3).Table 2
relative torque ripple rate as function different shifts.
α Δ C % α = 0 ° 8.48% α = 3 0 ° 4.96% α = 6 0 ° 6.54%Table 3
Effect of an angular offset on the operating characteristics of the DSIM.
Forα=0° Forα=30° Forα=60° Peak torque 23.75 N.m 22.05 N.m 17.85 N.m Average torque 12.37 N.m 11.36 N.m 9.22 N.m settling time for torque 2.63 s 3.01 s 3.38 s settling time for speed 2.49 s 3.02 s 3.43 sFigure 11
Characteristic electromagnetic torque.Figure 12
Characteristic of rotational speed.Figure 13
Torque-speed characteristic.


## 5. Discussion

The simulation results found when combining inverter-DSIM are similar to results found by simulating the DSIM supplied by a balanced source direct, but with the presence of harmonics at the electromagnetic torque, the torque-speed characteristic, of the stator and rotor currents and flux.These harmonics are due to the presence of voltage inverters; they are decreasing gradually as the number of levels of the inverter increases.In reality, static converters can provide a voltage (and current) cut (e), because the power electronics can be an electronic switching. To reduce the adverse effects of cutting the output voltage and thus tend slightly more to the "ideal converter", we shall increase the number of available output levels of the static converter. This will then reduce the amplitude of the voltage cut fronts, so the amplitude of the harmonic lines is induced by cutting.


## 6. Conclusion

This paper presents the impact of the multilevel NPC inverter sources on the circulating currents and voltage between the two stars of the DSIM. On the other hand, a higher level of output voltage of the multilevel inverter helps to improve a system of operation of the induction machine mainly by reducing the distortion of current consumption and electromagnetic torque ripple.The minimization of voltage and current induces the decreasing of the THD and improving of the energy quality, which means the reduction of the torque ripples and the cleaning of the stator currents. The PWM give the better results. The results obtained by simulation (mechanical, electrical, and magnetic characteristics) confirm the validity of theoretical studies on the machine and established models.



---

*Source: 6968023-2018-10-10.xml*