Performance Description of Brushless Doubly-Fed Induction Machine in Its Asynchronous and Variable Speed Synchronous Modes

Brushless Doubly-Fed Machine has attracted considerable attention in recent years due to its advantages. It has the robustness of the squirrel cage induction machine, and the speed and power factor controllability of the synchronous machine as well as the absence of brushes and slip rings, and using a fractionally rated frequency converter. Hence, there are considerable benefits over the conventional machines, when the machine is applied to applications such as a wind turbine generator or high power adjustable speed drive. However, these benefits are obtained in slightly more complex structure, higher cost and larger dimensions in comparison to the conventional induction machine. This paper presents fundamental aspects of the three modes of operation of brushless doubly fed machine, i.e. simple induction mode, cascade induction mode, and synchronous mode. The investigation is performed by analyzing the spatial harmonic contents of the rotor magnetic flux density. The direct cross couplings between stator and rotor fields as well as, indirect cross coupling between stator fields by the special rotor of this machine is described. Furthermore, loss analysis of the machine in various modes is presented and the torque-speed curves for asynchronous modes are obtained. A 2-D magnetodynamic finite element model based on the D-180 4/8 pole prototype machine is extracted and simulated to verify the results.


Introduction
The conceptual idea of the machines with two stator windings can be traced to some publications in the 1900s [1][2][3].From the stator winding point of view, dual-stator winding machines have been categorized as split-wound and self-cascaded [4].The split-wound dual-winding machine was proposed to improve the power capability of large synchronous generators and to achieve better drive reliability resulting from its inherent redundancy [5].The second type, the self-cascaded machine or the brushless doubly-fed machine (BDFM), with two sets of stator windings, which do not couple directly but via the intermediate action of a special rotor, was introduced by Hunt in 1907 [2].The special rotor structure increases the manufacturing cost and the complexity of the machine and produces undesirable spatial harmonics which decrease machine efficiency.But the brushless doubly fed machine shows commercial advantages because of the fractionally rated converter, and absence of the brush gear.The major application of BDFM is in wind generation, especially in harsh conditions, where it requires lower maintenance.Furthermore, it can be used as an adjustable speed drive [6,7].
The BDFM has two balanced three-phase windings on its stator.One of them is the power winding with P p pole pairs, which is directly connected to the grid, and the control winding with P c pole pairs, which is connected to the grid via a bidirectional fractional rated frequency converter (Figure 1).There are some rules on the selection of the stator windings pole number to avoid electromagnetic coupling between them [8].
The rotor should couple the power and control windings magnetic fields, in the appropriate manner.There are some suggestions about the rotor configuration [9].and Variable Speed Synchronous Modes The most popular structure is nested loop rotor, which was proposed by Broadway and Burbridge [10] (Figure 2).In this structure, the rotor consists of a number of nests and some loops in each nest which are short circuited through an end ring.The number of nests should be equal to the summation of the number of power winding and control winding pole pairs to provide indirect cross coupling between the fields of the stator windings (P r = P p + P c ).
The modeling and performance analysis of the BDFM has been examined in many publications.However, there is no comprehensive discussion on the electromagnetic phenomena that takes place and results in the three modes of operation.An alternative formulation of BDFM operation in synchronous mode has been published recently [11].The paper presents a precise characterization of all magnetic field components, which can be used in the machine design process and investigate the appropriate rotor configuration.
This paper presents the characteristics of the BDFM in simple induction, cascade induction and synchronous modes of operation and also, describes the direct and indirect cross coupling by simple electromagnetic relations.To simplify the analysis procedure and extraction of the relations, it is supposed that the nested loop rotor has one loop in each nest.As in the case of squirrel cage induction machine, increasing the number of rotor bars sdecreases the spatial harmonic contents of the rotor flux and its leakage inductance, and therefore, results in better performance characteristics.A 2-D magnetodynamic finite element model of the D-180 frame prototype BDFM is used for verification.

Magnetic Fields
The magnetic field density of the power winding (B p ) in its reference and in the rotor reference frames are stated in Equations ( 1) and ( 2), respectively.

  
where θ p and θ r are the mechanical angle around the air gap in power winding and rotor reference frames, respectively.ω p and ω r represent the angular frequencies of the power winding supply and rotor, respectively.The rotating magnetic field induces an emf in the rotor bars according to e   v Bl , which has the following distribution in the rotor reference frame, where l and r are the machine axial length and the air gap radius, respectively.As stated, the number of rotor nests equals to P p + P c .The arrangement of the rotor with one loop per nest and the nest span of are depicted in Figure 3 where N rb represents the number of rotor slots.

P N 
The emf induced in one loop can be calculated by subtracting two induced bar emfs as Equation (4).
With some simplifications, , 2  where This voltage produces a loop current through the loop impedance, where R l and L l are the resistance and leakage inductance of the loop, respectively.The maximum value of the loop current depends only on the rotor angular velocity, that is: ,max where The rotor currents produce magnetic motive force in the air gap, which can be calculated by the production of the current and turn function of rotor.The distribution of the k th nest MMF as a function of the rotor angular position in the rotor reference frame is shown in Figure 4.
Fourier series representation of this spatial MMF function has the following coefficients,

  
Fourier series of the rotor MMF can be obtained by the summation of MMF Fourier series of all nests as (13).
The production of the two terms containing θ k in Equation (13) can be written in summation form as, Substituting n = P c + hP r , h = 0, 1, 2, •••, in the first term of the right hand side of (14), the coefficient of k will be −2π(h + 1), which eliminates the dependency on the k factor, and therefore, the result of the summation will be non zero.Similarly, if n = P p + hP r , h = 0, 1, 2, •••, in the second term of the right hand, the coefficient of k will be 2πh, which also leads to non zero value for the summation.Hence, space harmonic orders n = P p + hP r and n = P c + hP r appear in the rotor produced MMF in the air gap due to power winding rotating field.All harmonics have the same time frequency.However, the rotation direction of n = P p + hP r and n = P c + hP r orders are in the same and the opposite direction, respectively, in comparison with the power winding flux.and Variable Speed Synchronous Modes Similarly, it can be shown that if the power winding is opened and in turn, the control winding is supplied, the rotor will produce an MMF with spatial harmonic orders n = P p + hP r and n = P c + hP r , which rotate in the opposite and same direction with respect to the control winding flux, respectively.

Direct and Indirect Cross Coupling
Among the rotor field space harmonics due to the field of the power (/control) winding, the P p (/P c ) pole pair harmonic field couples with the power (/control) winding field.This mechanism is called "direct cross coupling".Similarly, the rotor field harmonic which has P c (/P p ) pole pairs due to the field of the power (/control) winding, can couple with the control (/power) winding, if they have also same time frequency.However, this mechanism is called "indirect cross coupling".
The rotor MMF harmonics with P p and P c pole pairs corresponding to the power winding magnetic field are presented in Equations ( 15) and (16), respectively. , ) The amplitudes of these components are related just to the shaft rotational speed.Considering Equation (17), where g is the air gap width, the produced magnetic field densities can be written as Equations ( 18) and (19).
The mechanical angular frequency of B rp,Pc (t, θ r ) equals rotor.Hence, the field mechanical angular frequency from the control winding's point of view is, Considering N s as the number of stator slots, the distance between two consecutive slots is 2π/N s , and the number of slots in each pole and phase is N s /6P c .Considering one layer full pitch winding, the induced voltage in each phase of the control winding can be represented as Equation ( 23).
It is evident from Equation (23) that the amplitude and fre on occurs when one of duction mode, only the quency of the induced voltage are zero at the natural speed (ω r = ω n = ω p /p r ).If the voltage can produce a current, the air gap flux produced by the currents for ω r < ω n and ω r > ω n will rotate in the opposite and the same direction with respect to the rotor, respectively.

Simple Induction Mode
Simple induction mode of operati the stator windings is open-circuited and the other is supplied through a three phase source.This situation can take place in converter failures, for instance.In the following analysis, it is assumed that the control winding is opened, and it is proven that the machine can act approximately like an induction machine, however with poorer performance.

Torque
In the simple in   , , p r rp P B t com-ction torq ponent of the rotor field produces indu ue and leads the machine to operate as an induction machine.The other harmonic orders appear as rotor leakage flux with no average torque.As a result, performance of the BDFM is poorer than the squirrel cage induction machine of the same dimension.The electromagnetic torque can be calculated as follows, field behavior and the finite element method (FEM) to discrete the domain.The formulation uses the magnetic vector potential as an unknown parameter, the Galerkin method to obtain the set of equations that can be solved numerically, the Euler recurrence method to discrete the temporal derivatives, and the successive approximation and Newton-Raphson method to consider the nonlinear characteristic of the magnetic material.Movement is taken into account by means of the moving band technique, the Maxwell stress tensor and the mechanical oscillation equation.The simulation details are mentioned in [12].
The power and control windings have 2 and 4 pole pairs, respectively.The rotor has 6 nests with 3 loops per nest.Other specifications of the simulated machine can be found in [8].The rotor bars have been considered as a winding, they are connected with each other in an external electric circuit to form the nested loop rotor configuration, and it is then linked to the FE model to perform the simulations (Figure 5).
The control winding is open circuit and the power winding is supplied by a 50 Hz, 220 V rms voltage source.The applied load torque and the machine initial speed are zero.The machine starts to rotate and reaches a speed close to the synchronous speed of the power winding.
T is a function of the rotor speed as in the case of indu esults te element model is extracted 2-D distribution of the flux lines, magnetic field density in the air gap and its space harmonic components during starting are shown in Figure 6, respectively.It can be seen that there are all predicted pole pairs in the air gap field, including P c . Figure 7 shows the space harmonic spectrum of the air gap field for the steady state operation near the synchronous speed of the power winding (1500 rpm).The air gap flux density has clearly a 4 pole distribution.In fact, the speed of the power winding produced flux with respect to the rotor is very small, and therefore the emf induced in this speed is small, too.Hence, the amplitude of other space harmonic fluxes at the air gap is smaller than the start up period.ction machine.Finite element simulation results show that the starting torque is very low in this mode and the starting pr re takes a long time compared to the induction machine.Finally, the machine reaches a speed close to its synchronous speed, i.e. 1500 rpm.In the steady state, the machine torque ripple is higher than an induction machine because of the low number of rotor bar per pole which results in higher harmonic content.The steady state produced electromagnetic torque in this mode is transformer.Figure 10 shows the torque-speed characteristic in this mode.

Simulation
3) Induced EMF in the Control Winding The induced emf aft speed increases to the synchronous wer winding from standstill, is shown in Figure 11.From Equation (23), it was predictable that rotor speed reaches the natural speed (500 rpm), frequency and consequently, amplitude of the emf become zero.Then both these parameters start to increase again by increasing the rotor speed from the natural speed; however, the amplitude decreases under steady-state operation, because the emf induced in the rotor becomes very low, as mentioned

4) Losses
High harmonic content of the rotor flux density leads to high core loss and low efficiency.Figure 13 shows the core, rotor and stator copper losses.In the natural speed, the core loss is decreased, because of the zero speed of the P c pole pair flux produced by the rotor with respect to the stator.Similarly, very low speed of the power winding flux with respect to the rotor leads to reduction of copper loss under steady state conditions.

Cascade Induction Mode
This mode of operation occurs when one of the stator windings are supplied, while another is short circuited.In fact, BDFM in the cascade mode can be considered as two induction machines with rotors are connected together both electrically and mechanically, and behaves as an induction machine with P p + P c pole pairs.Furthermore, if both stator windings are excited, but the rotor speed is not appropriate to create indirect cross coupling between their fields and put the machine in its desirable mode of operation (synchronous mode), BDFM will operate as two separate cascade machines.This condition happens when the speed controller fails to stabilize the machine in synchronous speed.If the rotational speed of the rotor space harmonic component in response to one of the stator winding field (e.g.po r winding) with the same pole pairs as the other wi ng (e.g.control Induced emf in the control winding in simple induction mode when the sh we ndi and Variable Speed Synchronous Modes Supposing that the power winding is supplied and the control winding is short circuited, the The voltages generate current in shorted windings, and thus a rotating magnetic field in the air gap.The rotating field distribution, in the control winding and the rotor rence frames are presented in Equations ( 29) and ( 30 The phase and amplitude of the magnetic field are functions of the rotor speed.The amplitude is zero at ω p /P r and ω p /P p in the control winding and the rotor reference frames, respectively.For ω < ω n the magnetic field is clockwise (negative sequence), and for ω > ω n rotates anti-clockwise (positive sequence).Initially, the machine is supposed to operate below the natural speed.
The rotor bar induced voltage due to e magnetic field th of Equation ( 30) is in the form of, The corresponding current in the k th loop is presented in Equat The Fourier series of the MMF arises from the roto cu r nest Production of the two terms containing θ k in Equ ) can be expanded as,  ation (34 Similar to the previous section, according to Equation (3 n = P p + hP r (anti-clockwise) space harmonic components in its air gap magnetic field.It should be noted that the directions in the parenthesis are for speeds below natural speed.Otherwise, the rotation direction of the co re For ω < ω n , the P p pole pair harmonic of the rotor flux due to control winding field couples with the power winding field.Since they have identical direction and speed, it produces acceleration torque.However, for ω > ω n , the field rotates in the opposite direction of the power wi tio harmonic fluxes of the rotor and control winding rotate in opposite and same direction with respect to the P P pole pair harmonic fluxes, respectively, from the stationary oint of view.This is the reason why BDFM operates e ), the rotor can produce n = P c + hP r (clockwise) and ntrol winding magnetic field and these terms would be versed.
nding field, which produces braking torque.As menned before, for ω < ω n and ω > ω n the P c pole pair p n ar natural speed in this mode, and it operates similar to an induction machine with P p + P c pole pairs.The torque can be calculated similar to the previous section as Equation (36), As expected, the torque for ω r < ω p /P r is positive, for s negative, and for ω r = ω p /P r is zero.Furthermore, the torque will be zero if ω r = ω p /P p , because of zero emf induced in the rotor.
As stated, when both power and control windin fed while rotor does not rotate at synchronous speed, with synchronous speeds machine operating point would be between ω c /P r and ω p /P r and, one of the machines will act such as an applied load to the other.
Figures 14 and 15 illustrate the torque respectively, when the control winding is short-circuited.As observed, the rotor rotates near natural speed.
The magnetic flux distribution is shown in Figure 16.e machine has 6 poles.Furthermore, it is obvious that the air gap magnetic field contains 2 and 4 pole pair space harmonics, while other harmonics have noticeable amplitude, because the relatively high speed of the power winding with respect to the leads to high flux produced by the rotor.
The machine losses are divided into core as well as It is interesting to see the control winding current waveform while it is short circuited as shown in Figure 17.The control winding current amplitude and frequency are decreased when the speed increases to the natural speed, as expected.
Figure 18 shows the rotor current.The rotor current frequency (ω p − P p ω r ) has been decreased with increasing rotor speed.The rotor frequency in this condition is (ω p − P p ω r )/2π = 33.33,which is depicted in Figure 19.rotor and stator copper losses, which are shown in Fig- .
In this mode, the co erefore the currents are higher than simple induction mode.However, due to zero load torque, the steady state current and its corresponding losses are small.
The rotation speed is related to the load torque as in the case of conventional induction machine.Figures 21  and 22 show the electromagnetic torque and rotor speed in −10 and +20 Nm load torques, respectively.
The torque-speed characteristic in cascade mode of operation (ω n = 500 rpm) is depicted in Figure 23 using Equation (36).The two operating points of Figures 21 and 22 which are indicated on this figure confirm the accuracy of the analytical procedure.Furthermore, the machine behavior is like an inductio ine at the eed, and motoring and generating regions are clearly observable.The maximum torque in this mode is more than 50 times greater than the simple induction mode, which is the result of the high currents that was described before.
To show that the BDFM operates in asynchronous mode when both stator windings are excited, in the next simulation, the power winding is fed by a 100 V rms , 50 Hz and the control winding by a 220 V rms , 30 Hz voltage sources.As mentioned before, when the machine does not operate in synchronous mode, it works as two cascade machines, and the machine becomes stable between the natural speeds of two aforementioned machines.The

Synchronous Mode
The In this mode, the stator windings are supplied, and the rotor rotates at a certain speed, which enable it to couple the two stator windings magnetic fields.It should be noted that the two stator windings do not have any direct coupling with each other because of their different pole numbers.Initially, it is assumed that the machine operates as a motor.The two stator windings produce two magnetic fields in the air gap in the form of Equations (37) and (38), ro Each stator ma the rotor loops the air gap.The rotor flux in response to each stator field consists of two main terms with P p and P c pole pairs as well as other harmonics with relatively small amplitude.The part of the rotor field in response to one of the stator windings field couples with the field of the other stator winding with same pole pairs and produce electromagnetic torque.In fact, the fields of the two stator windings couple with each other indirectly through the rotor, which rotates at a specific speed.The second main part of the rotor flux due to one of the stator winding has equal pole pair numbers to it and couples with the field of that winding.So, there are synchronous and induction torques in this mode of operation.However, if the time frequency of the part of the rotor fluxes corresponding to a stator winding does not match the other stator filed frequency, no average torque will be produced, and the BDFM will work as two cascaded machines in one frame, which was described earlier.
The rotor magnetic field corresponding to the power winding field, which has P c pole pairs, and control winding field in the rotor reference frame are stated in Equations (39) and ( 40 Synchronous speed can be obtained by equalization of mechanical speed of the power winding field and the part of rotor field produced in response to the control winding, which has P p pole pair.
In synchronous mode, the torque is sum of the induction torque, which is produced due to the existence of direct couplings, and synchronous torque, resulting from the indirect cross coupling process.where θ c0 is the phase difference between the power and control windings in stator stationary reference frame. (44) The produced torque due to the coupling between the power winding field and P p pole pair part of the rotor field, which is produced in response to the power winding field (power winding direct coupling), is: The rotor shaft applied force due to coupling of the power winding field and P p pole pair part of the rotor field corresponding to the control winding field (control winding to power winding indirect cross coupling) can be calculated as, The resulting torque from this electromagnetic force can be evaluated from Equation (48).
Also, the torque related to coupling between the control winding field and P c pole pair part of the rotor field corresponding to control winding field (control winding irect coupling) is, d ding field and P c pole pair component of the rotor field corresponding to the power winding field (power winding to con indirect cross coupling) can be calculated as f Finally, the force applied to rotor shaft and its torque due to coupling of the control win trol winding ollows, rotate at the synchronous speed i.e. 800 rpm.
Figure 25 shows the air gap magnetic field, which has 6 poles distribution.It can be seen that the air gap magnetic flux density contains P p and P c components, i.e. harmonic orders of 2 and 4. High harmonic content of the air gap is acceptable, because both the stator windings are supplied and the rotor produces the fluxes in response to them.It should be noted that the stator fluxes electrical speeds with respect to the rotor are same and not small in general.
Figures 26 and 27 show the electromagnetic torque and speed, respectively.
Figure 28 shows the BDFM core, rotor, and stator copper losses.The copper losses are smaller than the cascaded induction mode.This is due to non-zero control winding voltage which results in fewer currents in the stator and rotor.
If the wind power leads the shaft to rotate at the speed the BDFM is used as a wind turbine genera- Clearly, the total torque is equal to summation of T 1 to T 4 .
To simulate synchronous mode, the stator power winding is fed by a 220 V rms , 50 Hz, and the control winding by a 100 V rms , 30 Hz voltage source.By regulating the load torque, the machine will be stabilized and tor, the frequency and amplitude of the induced voltage in power winding, which is directly connected to the grid, must remain fixed in spite of the wind speed variation, by changing the control winding frequency according to ω c = P r ω r − ω p .When BDFM is used as a wind turbine generator, the frequency of control winding should be regulated equal to ω c = P r ω r −ω p to achieve the desired rotor speed which results in maximum power tracking from the wind.
The control winding rotating field induces voltage and currents in the rotor loops, and acts similar to the field winding of the conventional synchronous machine.The rotor induced currents produce the rotor air gap flux with P p + hP r and P c + hP r space harmonic orders.The rotor field harmonic with P p pole pairs induces a voltage in power winding according to the following equations, Hence, the amplitude and phase of the induced emf in power winding is directly dictated by the control winding voltage via converter.Then, the active and reactive power flow between power winding and grid can be controlled by its phase and amplitude, respectively.
For generating operation simulation, the machine's shaft rotates at fixed speed of 800 rpm and the control winding is excited via a 30 Hz source.The induced emf in the power winding has a 50 Hz oscillation as shown in Figure 29.

Conclusions
Brushless doubly-fed machine has an interesting but slightly complex structure and performance with desirable characteristic in variable speed applications such as ve and mplete analytical peration of BDFM is presented based on simple electromagnetic relations.The mechanism of direct and indirect cross couplings for each mode is explained in details and simulated by means of a 2D magnetodynamic finite element model.The components of machine losses are obtained by FEM for each mode and their behavior with speed variation are explained.Also, the variations of losses portions are compared in three modes of operation.The torque-speed characteristics of asynchronous modes are depicted and the behavior compared with each other.Furthermore, the spatial harmonics in the air gap in three modes are illustrated and described.The performance of BDFM in synchronous mode relies on the pole number changing action by its nested loop rotor while the control winding field can affect the active and reactive powers of the PW.causing undesirable results in the machine behavior, for instance, increase iron loss, torque ripples, magnetizing currents, etc.Therefore, optimizing the structure of the machine is necessary to make it commercially applicable and bring it in

Figure 3 .
Figure 3. Nested loop rotor with one loop per nest.

Figure 4 .
Figure 4. Spatial distribution of the k th nest ampere turn.
first the induced emf in one coil with N c seri itioned at θ c angle should be calculated.

Figure 5 .Figure 6 .Figure 7 .
Figure 5.The configuration of rotor bars in one nest used in FE simulations.

Figure 8 .Figure 9 .
Figure 8. Electromagnetic torque in the simple induction mode of operation.

Figure 12 .
Figure 12.Phase sequence changing of the induced emf in control winding after the natural speed.

Figure 13 .
Figure 13.Core, rotor and stator copper losses in the simple induction mode.

Figure 16 .Figure 19 .Figure 20 .Figure 21 Figure 22 .Figure 23 .
Figure 16.Steady state operation in cascade induction mode (a) flux lines distribution; (b) magnetic field vs. θ r ; (c) space harmonic components of air gap field.er of (c) these two fields occurs, when the mechanical rotating speeds are equal.The mechanical speed of gnetic field induces voltage and current .These currents generate rotor flux in in clockwise direction.To provide cross coupling criteria between them, the following equality should be satisfied,

Figure 29
Figure 29.Ind as generator.uced emf in power winding at steady state operation