DC link parameters.

## Abstract

A harmonic analysis for a wind energy conversion system (WECS) based in a doubly fed induction generator (DFIG) is presented. A back-to-back frequency converter as voltage source for excite to rotor winding is used whereas than the electrical network excites to stator winding. The analysis is based on the induction machine model and the steady-state frequency converter model. Additionally, dynamic-state and steady-state models are developed and compared validating the results of the simulations obtaining a harmonic model, in steady-state, clear, and precise of the wind energy conversion system.

### Keywords

- harmonic analysis
- doubly fed induction generator
- back-to-back frequency converter
- electrical network
- harmonic and inter-harmonics

## 1. Introduction

Energy demand is increasing day by day around the world and renewable energy has proven to be an excellent choice of competition against conventional energy sources, which produce environmental pollution and contribute to climate change and global warming. An effective and proven type of renewable energy is wind energy, which achieved an increase of 63 GW in Mexico, reaching 433 GW in total capacity worldwide. This energy has been used by the human civilization for millennia to navigate and pump water for agricultural activities. Nowadays, the wind energy is used to produce electricity that leads the world of technology. This energy is used to impulse, mechanically, the generators that produce electrical energy. The device that allows this conversion is called wind turbine. There are different concepts regarding the coupling of the generator with the wind turbine being the most well-known fixed velocity and variable velocity wind turbines. The squirrel-cage induction generator was the first fixed velocity concept developed during the early years of wind turbines. This configuration is shown in Figure 1. The generator is directly connected to the electrical network but has the disadvantage of operating at low efficiency producing poor power quality. On the other hand, the double-fed induction generator (DFIG) is the technology that has the largest presence in the wind industry today as it has the advantage of have a higher energy efficiency compared to the squirrel-cage induction generator, a reduction in mechanical loads, a simple system of pitch control, and greater control of active and reactive power (see Figure 2). Due to the unpredictable nature of wind, wind energy still presents several challenges such as fluctuations in voltage generated and complications to achieve efficient couplings with the electrical network.

These problems have led to develop strategies to solve these mishaps. However, the devices used by these strategies produce certain alterations in the energy signal required by the generator to function properly. This implies variations in the energy delivered by the generator, which leads to a further challenge in the area of wind energy: ** the power quality**. This concept describes how close the energy is to the appropriate standards to operate and allow the proper use of equipment that employ this energy. Variations in voltage can result in malfunction of sensitive components of system, a reduction in generator efficiency and in general, degradation of power quality, while harmonics, which are current and/or voltages signals occurring in electrical systems with multiple frequencies of the fundamental, can cause severe problems such as generator induction effects and torque oscillations reducing the life of generator. Therefore, it is very important to obtain a clear harmonic analysis of the wind energy conversion system (WECS) when they are excited by non-sinusoidal voltage sources since the wind industry currently lacks models for harmonic propagation studies; this is done using the Fourier transform tool described later in the chapter.

## 2. Harmonic analysis

In the beginning, Fourier series emerged as a method for the explanation of physics phenomena, at a time when there was no knowledge of electronics or telecommunications, but they are also very useful in our times applied in computer science and communication. This is possible since electrical phenomena are also physical phenomena and their behaviors explained through periodic functions or periodic wave movements, where a Fourier transform allows obtaining discontinuous waves. In addition, the harmonic analysis is the process of calculating the magnitudes and phases of the fundamental frequency and all the higher order of a periodic waveform. The resulting series is known as the Fourier series [1], which for the periodic function

This constitutes a representation in the time-domain of the periodic function. In this expression,

*is the period of the function*T

*harmonic. These signals can be visualized in a 3-D system in which their magnitude, frequency location, and over time are represented (see Figure 3).*n-th

The constant term of the Fourier series is given by:

which represents the area under the curve of

*, that is the average value of the function. While*T

The above equations are simplified when the waveform has symmetry. If the waveform has odd symmetry, this is

*, and*n

The Fourier series for an odd function will have only sine terms. If the waveform has even symmetry, this is

*, and*n

### 2.1. Transform Fourier

Formerly the process to find the Fourier transform was a bit expensive, due to the calculations that had to be done manually, nowadays with the help of technology and computers, different algorithms have been developed to simplify and make this process easier, besides of the discrete Fourier transform (DFT) and the fast Fourier transform (FFT). The main idea of the Fourier transform is to transform a signal from the time domain or space to the frequency-domain. In addition, through inverse Fourier transform, you can pass functions from the frequency-domain to the time-domain. The Fourier series therefore represents the special case of the Fourier transform applied to a periodic signal. In practice, the data are frequently available in the form of a time-sampled function, represented by a time series of amplitudes separated by fixed time intervals of limited duration. When Fourier analysis is applied to a periodic continuous signal in the time-domain produces a series of discrete frequency components in the frequency-domain [2]. The following equations form the pair of the Fourier transform:

where * t* is time,

*is the frequency,*f

*is the test signal,*x(t)

*is spectrum as a function of frequency*X(f)

*Eq. (7) is called the ‘forward transformation’ and Eq. (8) is the ‘reverse transform’. Generally,*f.

The real part of

In analogous form, the imaginary part of

The amplitude spectrum of the frequency signal is obtained from:

In addition, the phase spectrum is:

### 2.2. Fourier discrete transform

When frequency-domain spectrum is a sampled function, as well as the time-domain function, a Fourier transform pair of discrete components is obtained [3]:

Equation (14) can be written in condensed form as:

In Eq. (15),

*components of the function in the frequency-domain, while*N

*samples of the function in the time-domain. The calculation of the*N

*frequency components of*N

*time samples requires a total of*N

### 2.3. Fourier fast transform

The implementation of Eq. (14) involves a number of complex additions and multiplications that is proportional to

*by*k

*in this term is not normally counted as part of the implementation. The fast Fourier transform is an algorithm that allows calculating the discrete Fourier transform of a harmonic signal with a substantial saving of calculations. For the case, a number of samples equal to*n

*is integer, only*k

*) multiplications [4] are required. The saving in calculations increases as the number of samples increases with the ratio:*N

Thus for

^{3}samples, which reduces the number of operations 205 times. The algorithm that is proposed is based on the method called ‘successive folding’. Knowing that

*is a positive number, so*M

which can be expressed as

If it is defined as

For

For

Also, given that

A careful analysis of the last equations shows some interesting properties of these expressions. Note that an N-point transform can be calculated by dividing the original expression into two parts, as indicated in the last two equations. The calculation of the first half of

*positive integer, it can be shown that the number of complex operations (multiplications and additions) is given by:*n

These expressions indicate the number of multiplications

And

The application of this mathematical tool are varied, for electrical engineering, electronics, telecommunications and other fields of information sciences. All these sciences can to demonstrate that the sinusoidal signals can represent the sending of data by electrical transmission networks. In the following sections, a model of the doubly fed induction generator for harmonic analysis based on the Fourier transform is presented. Initially, the mathematical model of DFIG is presented at the fundamental frequency. Then, the model is performed at harmonic frequencies to obtain the harmonics components generated in the electrical variables of the generator. Finally, the results are compared and summarized in tables.

## 3. Wind turbines technologies

The connection of wind farms to an existing transmission or distribution network can affect the behavior of the same depending on two main variables:

Installed wind power.

‘Strength’ to network to which its connects (expressed as PCC/

).

This analysis must be done both from the point of view of static and dynamic behavior of the system. Unlike other types of generation technologies, it is not usual for wind generation to cause problems in terms of the harmonic, inter-harmonic and sub-harmonic distortion. This is due both to the fact that most of the machines used are of the asynchronous type, that is, the technology type of wind energy conversion system used in wind farms determines the harmonic distortion intensity. For example, the fixed velocity wind turbines based in squirrel-cage induction generator, the element that causes current distortion is the soft-starter, formed by a thyristor bridge, which not in all operations of the turbine is in operation. On the other hand, variable-velocity wind turbines today have in common the use of power electronic converters to be able to vary the velocity and to maintain the power factor within the desired limits. However, these devices produce harmonic signals in the generator and the harmonic distortion magnitude depends mainly on the type of inverter or converter and its control because these use the IGBT thyristors (isolated gates bipolar transistors). The switching frequency in these inverters is limited to the range of 2–3 kHz. The frequency converters affect waveform quality and can negatively affect other consumers connected to the electrical network.

### 3.1. WECS using DFIG with conventional back-to-back frequency converter

Considering the different types of WECS configuration, a wind turbine based on a doubly fed induction generator (DFIG) connected directly to the electrical network via a frequency converter is analysed. Among electronic power converters there are three types widely used in the wind power market: * back-to-back* converter, multilevel converter and the matrix converter [5]. In this section, the revision is made around the network-connected WECS-based DFIG and a

*converter attached at a point common coupling (PCC). A rotor-side converter (RSC), a network-side converter (NSC) and a decoupling capacitor between the two converters constitute this converter. Notice that the harmonic distortion is a phenomenon that causes problems for both the users and the entity responsible for the provision of the electric power service [6], so it is considered necessary an appropriate understanding of the techniques currently used to mitigate this phenomenon. In this sense, the literatures that define the harmonic sources in WECS are classified as follow: Harmonic distortion from the network and harmonic distortion form the WECS:*back-to-back

The harmonic distortion from the electrical network: Since the WECS is connected to the electrical network; it becomes sensitive to any harmonic distortion coming from it. Distorted voltages in the stator generate harmonic currents in the WECS, which in turn can generate ripples voltage, torque pulsations, overheating, and increased losses in stator windings, under power factor, among others [7]. In some cases, these affections can have destructive effects on the WECS [8]. The solutions found in the literature involve the use of custom power devices as well as various control strategies implemented in the frequency converters.

The harmonic distortion from the WECS: The harmonics coming from the WECS is another of the harmonic distortion scenario, which the injection of harmonic components of voltage or current into the electrical network can cause problems such as excessive heating of transformers and conductors [9]. The main causes of harmonics from the WECS are: harmonics produced by nonlinear loads inherent to the generation system and harmonics produced by the WECS itself that is to say, associated to the operation of the components (generator, electronic power converter).

#### 3.1.1. Operating principle of doubly fed induction generator

As in traditional three-phase induction generator, the doubly fed induction generators have a magnetic field formed by the stator phases rotating at synchronous velocity

*is the poles number of generator,*p

*is the synchronous velocity in*n

_{s}

*is defined as the relative velocity of the stator fields with respect to the rotor velocity*s

The concept ‘slip’ is very important because it describes the machine operation. Now, with the above equation, we can establish in

where

As can be seen in the equations just described, the velocity at which the rotor rotates depends exclusively on the network and rotor frequencies, and not on the torque imposed by the turbine, as is the case with other types of generators.

#### 3.1.2. Steady-state model of DFIG at harmonic frequencies

The general scheme of doubly fed induction generator (DFIG) in the stationary reference-frame uses a wound rotor with slip rings to transmit current between the stator and rotor voltage sources is shown in Figure 4. The voltage equation that describes the equivalent circuit is given by

where

*, it is necessary to determine if the harmonic is a positive, negative or zero sequence. All multiples of three in a three-phase system are zero sequence, that is, corresponding to the harmonics*h

The sign

*in a three-phase system to be in sequence positive, the following relation must be fulfilled:*h

Since * h* must be an odd number,

*must belong to natural numbers; this implies that*n

*must be a multiple of 2*k

In order for a harmonic of order * h* in a three-phase system to be in sequence negative, the following relation must be fulfilled

Since * h* must be an odd number,

*must belong to natural numbers; this implies that*n

*+1 must be a multiple of 2,*k

* k* + 1 = 2

*; ∀*m

*∈*m

*→*N

*= 3*n

*→ h = 6*m

*− 1; ∀*m

*∈*m

*,*N

*≠ 0*m

In summary:

Harmonics of positive sequence:

Harmonics of negative sequence:

The stator and rotor phasors currents are obtained from the above equation:

Then,

where

#### 3.1.3. Analysis of back-to-back converter at harmonic frequencies

One of the most used equipment in the industry for the control of doubly fed induction generators (DFIG) is the AC-AC three-phase converter. In such equipment, the voltage generated in the DC link

The controlled rectifiers allow a bidirectional flow of power and through appropriate control and modulation; techniques generate virtually sinusoidal input currents. The power factor of this type of rectifiers can be adjusted at will. It should be mentioned that in these devices, the harmonics present in the current are the same as those present in the voltage. Therefore, the ripple that appears on the voltage, just like in the current, is rippled from high frequency due to the switching function. To determine the voltage that is obtained from the inverter, first obtain the voltage in the capacitor according to the following expression:

Solving, in steady state, with

To obtain the value of the inverter current, it is necessary to know the capacitor current

where the switching function * m* Є {−1, 0, 1} is a discrete variable of the modulation SPWM, which is defined as

whose mathematical representation can be obtained by decomposition in Fourier series by

where ω represents the network frequency;

where * P* represents the power in inverter and

By analyzing only, the harmonic components, previous expression is reduced to

In addition, for * b* and

*phases, the independent products can be calculated by*c

Using the trigonometric identity:

Finally, the sum defined by the above equations is represented by

The same procedure is performed for the network-side converter. The DC link parameter is shown in Table 1. Figures 5 and 6 show the voltage generated and harmonic frequency of the rotor side converter.

Parameters | Value |
---|---|

DC link capacitor | 32.5 F |

DC link voltage | 180 V |

Base voltage | 176.5 V |

Base current | 33 A |

Base power | 22.5 kVA |

Figures 7 and 8 show the voltage generated and the harmonic order of the network-side converter. The voltage generated by the converter and shown in Figure 5 will be the excitation source of the rotor winding. The harmonic slip of rotor-side converter is shown in Table 2.

Harmonic | 5 | 7 | 11 | 13 | 17 | 76 | 80 | 158 | 159 |
---|---|---|---|---|---|---|---|---|---|

Sequence | — | + | — | + | — | + | — | + | — |

Harmonic slip | 1.21 | 0.96 | 1.32 | 0.97 | 1.25 | 0.98 | 1.01 | 0.993 | 1.07 |

## 4. Model validation

To validate the proposed model a three-phase generator of 37.5 kW, 460 V, and 46.8 A is used. The induction generator parameters are shown in Table 3. Noted that the results obtained from the proposed model (steady-state) are compared with those obtained from the simulated complete model (dynamic) ones the steady-state has been reached. For this case study, a

### 4.1. Case study 1: DFIG excited with conventional back-to-back converter

For this case, the generator is excited in the stator windings with a sinusoidal three-phase balanced voltage source of 460 V at 60 Hz shown in Figure 9. Rotor windings are excited with a quasi-sine voltage source of 50 V at 45 Hz. The voltage waveforms of rotor winding are shown in Figure 5. The stator and rotor current waveforms obtained from dynamic-state are shown in Figures 10 and 11, respectively. An expansion of these currents is made once they have reached their steady-state. The electromagnetic torque and the nominal velocity reached by the generator are shown in the Figures 12 and 13, respectively.

The current waveforms of stator and rotor are obtained from steady-state model and shown in Figures 14 and 15, respectively. In addition, an expansion of these currents is made and shown in Figures 16 and 17.

Finally, a comparison is made of the forms obtained in both models, both in magnitude and in phase angle, concluding that both models coincide in these aspects. The graphic result of this is shown in Figures 18 and 19, respectively.

## 5. Incorporation to the electrical network

The integration of wind farms into the network is not an easy problem to solve. Numerous studies have been carried out related to the impact of wind penetration on electric power systems, including those related to stability studies, response of electrical protections and harmonic effects in the electrical network. In some power systems, the penetration of wind energy is rapidly increasing and begins to influence the overall behavior of the system and its stability. In Figure 20, it is possible to observe the equivalent diagram of the system to be simulated, which corresponds to the simplification of an induction generator, a frequency converter and typical distribution system. It is presented a three-phase circuit that feeds an industrial load, consisting of a doubly fed induction generator (DFIG), controlled by a back-to-back converter (which injects the harmonic current components).

The presence of harmonics in an electrical power distribution system can bring serious problems, both for the same equipment or system that produces them as for other equipment connected in the same network. The generated harmonic currents interact with the system impedance and cause the voltage distortion, increase the losses and produce effects of overloads connected to the electrical network. Furthermore, they reduce the power factor, given its reactive nature. The presence of harmonic distortion will be more appreciable in some points than in others. However, we can establish that, in general, the effects of harmonics depend on:

the elements that generate harmonics, their diversity within the system, the spectral content they inject and their operating regime in time-domain

the quantity, configuration and values of the elements in the electrical network, both in the time-domain and frequency-domain;

the sensitivity of the equipment to the presence of harmonic currents or voltages;

the degree of interest shown by electricity companies and commercial and industrial companies for the harmonic problems that may arise;

the establishment of norms to set individual total harmonic distortion limits and their compliance on the part of the sectors involved; and

the effectiveness of standards over time, depending on the technological cases registered in the manufacture of new equipment.

In the study of the effects of harmonics in the distribution systems, it is important to know the propagation paths of the harmonic currents, whose important effect in the system consists in the distortion of the voltage in the different nodes, transferring the harmonic presence to other subsystems. The presence of harmonic voltages results in currents with those harmonic frequencies, which cause additional losses. The induction generator operating with non-sinusoidal voltages cause additional losses due to harmonics. The harmonic losses that occur are losses in stator copper, losses in the core, mechanical losses and losses in the rotor. The largest losses occur on rotor bars (for squirrel cage generator). In addition, harmonic losses are independent of the generator load, but may be much greater than losses due to a pure sine wave, depending on the harmonic content of the applied voltage waveform. These factors raise the temperature of all generator components reducing the efficiency and the useful life. This elevation in temperature may result in excessive heating and may not meet the demands of the load coupled to its axis, or resulting in the destruction of the generator. Harmonic currents cause increased losses in copper for to the skin effect due to their high frequencies, cause insulation stresses, Eddy current losses in the conductive parts and possible resonance at a harmonic frequency between the inductive reactance of the transformer and the capacitance of the line or capacitors for power factor improvement.

The effect of the losses, then, is an increase in the heating of the transformer, which may be significant. Due to behavior of triple harmonic currents in three-phase systems, the three currents of line are in phase, and when they reach the center of a wye connection, they will return through the neutral or ground. This fact, although the three-phase system is perfectly balanced, can cause the heating of cables and connectors when its ampacity is exceeded.

### 5.1. Power-quality problems of WECS connected to the electrical network

The term ‘power quality’ refers to voltage and frequency stability, and the absence of various forms of electrical noise (e.g., flicker or harmonic distortion) in the electrical network. Generally, electric companies (and their customers) prefer an AC with a sinusoidal waveform, free of interruptions or disturbances. The consequences of large-scale supply incidents are well documented. A recent study in the United States has shown that digital industrial and commercial firms are losing $ 45.7 billion per year because of power network disruptions. In all commercial sectors, it is estimated that between 104,000 and 164,000 million dollars are lost due to interruptions and another 15,000–24,000 million dollars because of other power quality problems in the network. Many problems in the supply originate in the electric network that with its thousands of miles of transmission is subject to climatic conditions such as hurricanes, storms with snow, ice and floods, along with equipment failures, traffic accidents and major operations of connection.

The power quality study and ways to control it is a topic of interest for electricity suppliers, large industrial companies, businesses and even residential users. The study has intensified as teams have become increasingly sensitive to even minimal changes in voltage, current and frequency in the supply. The Institute of Electrical and Electronics Engineers (IEEE) has developed a standard that includes definitions of energy disturbances. The standard (Standard 1159–1995 of the IEEE, IEEE Recommended Practice for Monitoring Electric Power Quality) describes many problems of the power quality. Regarding the wind sector, there are regulations that regulate the limits for this type of problems, as indicated by the IEC 61400-21 standard. The IEC–61400-21 standard provides recommendations for preparing the measurements and assessment of power-quality characteristics of wind turbines. The measurement procedures described in the IEC61400-21 standard are valid to test the power-quality characteristic parameters for the full operational range of a wind turbine, connected to a MV or HV network with fixed frequency within ±1 Hz, sufficient active and reactive power regulation capabilities and sufficient load to absorb the wind power production.

#### 5.1.1. Start and stop of the wind turbine

Most of the electronic wind turbine power controllers are programmed so that the turbine runs in no-load at low wind speeds (if it were connected to the network at low wind speeds, in fact it would work as an engine). Once the wind becomes strong enough to rotate the rotor and generator at its rated speed, it is important that the turbine generator be connected to the power network at the appropriate time. If this is not the case, only the mechanical strength of the multiplier and the generator will be present to prevent the rotor from accelerating, and eventually it will be packed (there are various safety devices, including fail-safe brakes, in the case that the of correct start fail).

#### 5.1.2. Soft-started with thyristors

If a high-power wind turbine is connected to the electrical network by means of a normal switch, a partial dimming would be observed (due to the current required to magnetize the generator), followed by a power peak, due to the generator current overloading. Another unpleasant side effect when using ‘hard’ switches would be to apply extra wear to the multiplier, since the generator connection would act as if the mechanical brake of the turbine were suddenly activated. To avoid this situation, modern wind turbines have a soft-start; they are connected and disconnected from the electrical network gradually by means of thyristors, a type of semiconductor continuous switch that can be controlled electronically. Thyristors lose about 1–2% of the energy that passes through them. Thus, modern wind turbines are usually equipped with a switch called ‘drift’, that is, a mechanical switch that is activated after the turbine has made the soft start. In this way, the amount of energy lost is minimized.

#### 5.1.3. Weak electrical network

If a turbine is connected to a weak electrical network (i.e., it is very far in a remote corner of an electrical network with a low-power transport capacity), there may be some problems of partial obscuration and/or power overvoltage. In these cases, it may be necessary to reinforce the electrical network to transport the AC from the wind turbine. The local power company is experienced in dealing with these voltage problems, as they are the exact reflection of what happens when a large user (e.g., a factory with large electric motors) is connected to the electrical network.

#### 5.1.4. Islanding

‘Islanding’ is a situation that can occur if a section of the electrical network is disconnected from the main electrical network, as would happen by accidental or intentional tripping of a large circuit breaker in the network (e.g., due to stoppages in the network electrical supply or short-circuits in the network). If the wind turbines continue to operate in the part of the network that has been isolated, it is very likely that the two separate networks are not in phase after a short interval of time. The restoration of the connection to the main power network can cause huge overcurrent’s in the wind turbine network and generator. This would also cause a great release of energy in the mechanical transmission (i.e., in the axes, the multiplier and the rotor), as would a ‘hard connection’ from the generator of the turbine to the electrical network. Thus, the electronic controller will have to constantly monitor the voltage and frequency of the AC of the network. In the event that the voltage or frequency of the local network are within certain limits for a fraction of a second, the turbine will automatically disconnect from the network, and immediately afterwards it will stop (normally activating the aerodynamic brakes).

The distributed generation (generation of electricity through of solar panels, wind farms, etc.), is increasingly having a greater presence in the world, in the case of the sector of self-consumption and injection of surplus energy into the network. This factor causes that the energy systems supplier has more advanced control means to verify the variations of loads that cause important fluctuations in the parameters of the voltage, injection of harmonics, change in the power factor, greater technical losses, among others. Therefore, it is advisable to make a power-quality study in an industry, commerce or home. This analysis must be done with high-resolution RMS equipment in order to detect and accurately measure the different parameters mentioned.

### 5.2. Case study 2. DFIG excited with three-phase back-to-back converter SPWM to improve the power quality

For this case, the generator is excited in the stator windings with a sinusoidal three-phase balanced voltage source of 460 V at 60 Hz shown in Figure 9 and rotor winding is excited with a quasi-sine voltage source of 50 V at 45 Hz to which the SPWM switching technique was applied. The stator and rotor current waveforms obtained from dynamic-state are shown in Figures 21 and 22, respectively. An expansion of these currents is made once they have reached their steady-state.

The current waveforms of stator and rotor are obtained from steady-state model and shown in Figures 23 and 24, respectively. In addition, an expansion of these currents is made and shown in Figures 25 and 26.

Finally, a comparison is made of the forms obtained in both models, both in magnitude and in phase angle, concluding that both models coincide in these aspects. The graphic result of this is shown in Figures 27 and 28, respectively.

### 5.3. Analysis of simulation results

The doubly fed induction generator (DFIG) interconnected to the electric network is simulated based on a steady-state model. This proposed model consists of the integration of induction generator, frequency converter and the electric network. Each of those sections is then analysed.

#### 5.3.1. Rotor-side non-sinusoidal voltage source

Once the harmonic slip for each harmonic signal is determined, the frequencies generated in the rotor-side non-sinusoidal voltage source are obtained as shown in Table 2. Figure 29 shows the frequency spectrum of rotor side non-sinusoidal voltage source. There are harmonics of the order

#### 5.3.2. Electrical network analysis

The network voltage could only be influenced by non-ideal voltage modulation, if the power of the modulated source is in the same class as the network. However, the harmonic current flow over the network generating always-harmonic voltage, which shows undesired effects in the network, independent from the power of the harmonic current source. This means, that harmonic currents are always of high interest for the resulting power quality of the electrical network. For this case, a harmonic voltage and harmonic current of the network is considered. Figures 30 and 31 present the network current and the network voltage, respectively.

#### 5.3.3. Stator and rotor current analysis

The stator current waveform is shown in Figure 23 where it is observed that it contains considerable harmonic components. The interaction between the rotor flux harmonics and the fundamental flux in the air-gap induces currents in the stator which are not integer multiple of the supply frequency. These are commonly called inter-harmonic with stator current fluctuation. The harmonic in the rotor are induced in the stator due to the dominant rotor harmonics (harmonic numbers are fifth and seventh).

The behavior of the rotor current in steady-state is shown in Figure 25. The ripples on the rotor current and its large rate of variation result in pulsed voltage waveforms. The rotor harmonic current establishes rotating magnetic fields in the air-gap, inducing currents of corresponding frequencies to the stator winding. It is known that the 6 * h* + 1 harmonics of the rotor current set up positive sequence magnetomotive force (MMF) waves, rotating in the same direction as the fundamental rotor MMF, whereas the

*= 0,1,2,3…. Similarly, the rotor side-band frequencies due to the rotor-side converter harmonics also set up MMF waves, rotating at [(1 ± 6*h

*)*h

*± 6*s

*]*m

*, with respect to the rotor. The frequency of the induced harmonic currents to the stator is then: |1 ± 6*f

_{es}

*± 6*hs

*|*

_{m}

*,*f

_{es}

*= 0,1,2,3…. Since the slip*h, m

*is non-integer, the harmonic content of stator current consists mainly of sub- and inter-harmonics, which create undesirable effects on the supply system. For instance, low-frequency sub-harmonics appear as unidirectional components superimposed on the phase currents, while sub- and inter-harmonics neighbouring to the supply frequency may create a beat effect on the stator current magnitude. The rotor harmonic current due to the switching of the transistor appear at frequencies*s

*= (6*f

_{er_h}

*± 1)|*h

*|*s

*,*f

_{es}

*= 0,1,2,3…, that is, at the 5th, 7th, 11th, 13th … of the fundamental rotor frequency |*h

*|*s

f

_{es}

*= 0 for the fundamental.*h

The DC current distortion introduced by the network-side converter (at the frequencies

*= |(6*f

_{es_h}

*± 1)*h

*± 6*s

*|*m

*,*f

_{es}

*= 0,1,2,3….*h,m,

It explains only the harmonic frequencies at 20 Hz (fundamental rotor frequency), 100, 140, 220, 260, 340 and 380 Hz, corresponding to * h* = 1, 5, 7, 11, 13, 17 and 19. For instance, the spectral lines at 280 and 320 Hz result from the interaction of the 20 Hz fundamental rotor frequency with the basic inverter DC harmonic at 300 Hz

### 5.4. Case study 3. High-power DFIG excited with back-to-back converter

Finally, a last case study of the doubly fed induction generator is presented. To validate the proposed model, a three-phase generator of 372.8 kW, 2300 V and 93.8 A at 60 Hz is used. The rotor windings are excited with a quasi-sine voltage source of 500 V at 45 Hz. The induction generator parameters are shown in Table 4.

For this case study, a

The electromagnetic torque and the nominal velocity reached by the generator are shown in the Figures 34 and 35, respectively. The current waveforms of stator and rotor are obtained from steady-state model and shown in Figures 36 and 37, respectively. Also, an expansion of these currents is made and shown in Figures 38 and 39.

Noted that the results obtained from the proposed model (steady-state) are compared with those obtained from the simulated complete model (dynamic) ones the steady-state has been reached (Table 5).

## 6. Conclusion

Based on the results of this research, we propose a wind turbine model based on a double-fed induction generator for harmonics propagation studies. This model consists of three stages: analysis of the induction generator, analysis of the frequency converter and analysis of the electric network. It should be mentioned that the models were developed in steady-state and dynamic-state. Finally, the three stages are integrated for the development of a single model for the harmonic analysis of the double-fed induction generator interconnected to the electrical network. The results indicate harmonic and non-harmonic currents can exist in both generator windings depending on the slip and the fundamental frequency in both voltage sources. The results of the steady-state simulation were compared with the results in dynamic-state model obtaining great similarity in the waveforms of the currents in the stator and rotor, in both magnitude and phase angle resulting in a proper model for harmonic and non-harmonic analysis of the doubly fed induction generator (DFIG) which can be used for ‘harmonic’ analysis in electrical power systems.

## References

- 1.
Arrillaga D, Bradley A, Bodger PS. Power System Harmonics. 2nd ed. California: John Wiley & Sons; 1989. p. 500 - 2.
Randall RB. Application of B and K Equipment to Frequency Analysis. 1st ed. USA: Bruel and Kjaer; 2003. p. 485 - 3.
Kreyzig E. Advanced Enginnering Mathematics. 6th ed. USA: John Wiley & Sons; 2010. p. 806 - 4.
Bergland GD. A Fast Fourier Transform Algorithm for Real Valued Series. Communications of the ACM. 1st ed. England; 1986. p. 710 - 5.
Carrasco JM, Franquelo LG, Bialasiewicz JT, Member S, Galván E, Guisado RCP, Member S, Ángeles M, Prats M, León JI, Moreno-Alfonso N. Power-electronic systems for the grid integration of renewable energy sources: A survey. IEEE Transactions on Industrial Electronics. 2006; 53 (4):1002-1016 - 6.
Arcega FJ, Pardina A. Study of harmonics thermal effect in conductors produced by skin effect. IEEE Latin American Transaction. 2014; 26 (1):1488-1495 - 7.
Hu J, Nian H, Xu H, He Y. Dynamic modeling and improved control of DFIG under distorted grid voltage conditions. IEEE Transaction Energy Conversion. 2011; 26 (1):163-175 - 8.
Amini J. Novel control method of grid connected converter of doubly fed induction generator to achieve disturbances rejection. In: International Conference on Environment and Electrical Engineering. Washington, DC; 2011. pp. 1-5 - 9.
Mora-Barajas M, Bañuelos-Sánchez P. Contaminación armónica producida por cargas no lineales de baja potencia: modelo matemático y casos prácticos [dissertation]. Madrid: 2010. p. 111 - 10.
Comech MP. Análisis y Ensayo de Sistemas Eólicos ante Huecos de Tensión [thesis]. Zaragoza: Universidad de Zaragoza; 2007. p. 100 - 11.
Boldea. Synchronous Generators. 1st ed. USA: Taylor & Francis Group; 2006. p. 256 - 12.
Cuevas M. Control Vectorial de un Motor de Inducción Doble Alimentado [thesis]. Mexico: CENIDET; 2003. p. 113 - 13.
Gamarra. Control de un Generador Doblemente Alimentado para Turbinas Eólicas [thesis]. Universidad Politécnica de Cataluña: Ingeniería Técnica Industrial; 2009. p. 110 - 14.
Sierra E. Análisis del diseño y control de un generador trifásico doblemente alimentado [thesis]. Chile: Ingeniería Civil Electricista; 2012. p. 152 - 15.
Hernandez E, Madrigal M. A step forward in the modelling of the doubly fed induction machine for harmonic analysis. IEEE Transaction Energy Conversion. 2014; 29 (3):149-157