230 American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) ISSN (Print) 2313-4410, ISSN (Online) 2313-4402 © Global Society of Scientific Research and Researchers http://asrjetsjournal.org/ Wind Energy Conversion System for Harmonics and Inter-harmonic Analysis Emmanuel Hernándeza*, Miguel Ángel Hernándezb, Hugo Jorge Cortina Marreroc, Reynaldo Irachetad a,bUniversidad del Istmo, Calle Sta. Cruz s/n Av. Universidad, Oaxaca. C.P. 70760, México aEmail: emanuel.mayoral7@gmail.com; bEmail: mahluni2000@yahoo.com.mx cEmail: hcmlhd@yahoo.com; dEmail: riracheta@conacyt.mx Abstract In this research, a harmonic and inter-harmonic analysis is presented for two types of wind energy conversion systems (WECS): a squirrel cage induction machine and doubly-fed induction machine because of the incorporation of power electronics in these systems, such as soft-starter and power converter that distort the voltage and current waveform. The analysis is based on the well-known steady-state induction machine model, additionally; a dynamic-state model is developed to compare the steady-state model validating the results of the simulations obtaining a harmonic and inter-harmonic model, in steady-state, clear and precise of the wind energy conversion system. Keywords: Doubly-fed induction machine; harmonic and inter-harmonic analysis; wind energy conversion system. 1. Introduction The increase in power generated in wind farms makes it necessary to study the power quality injected into the electrical network, transmission systems and electrical machinery. This power generated can present problems mainly due to three factors: 1. Introduction of voltage fluctuations and the presence of harmonics and inter-harmonic in the electrical network due to natural oscillations in wind velocity, and those caused by the presence of the mast. ------------------------------------------------------------------------ * Corresponding author. http://asrjetsjournal.org/ American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 231 2. Stability problems of the wind turbines due to contingencies in the electrical network (short-circuits, lightning, maneuvers, etc.), as well as those due to the great variability of the wind. 3. Problems of prediction and planning in power generated. Currently, one of the problems that has the wind turbine connected to the electrical networks is the intermittence of the injected power to network by the wind dependence and, on the other hand, the injection of harmonic currents to the electrical network due to the use of power electronics for their connection. The type of technology of wind turbines conversion system used in wind farms determines the power quality. For example, the fixed-velocity wind turbines based in squirrel cage induction machine, the element that causes distortion of the current is the soft-starter, formed by a thyristors 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 electronic power converters to be able to vary the velocity and to maintain the power factor within the desired limits. This is done with a switching control that applies a pulse width modulation (PWM) for velocity and power factor control, as mentioned. The magnitude of the harmonic distortion 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 kHz to 3 kHz. The new turbines that now reach up to 5 MW of nominal power, the converters used are medium voltage, with a frequency of only 150 Hz to 250 Hz. As consequence, the frequency decrease being greater the distortion of the wave, or as it is the same, the emission of harmonic currents to the system is greater. These technologies cause harmonics (integer multiples of the fundamental frequency) and inter-harmonics (voltages or currents of frequencies that are not integer multiples of the system fundamental frequency which may appear as discrete frequencies or as a broadband spectrum) because they are loads that the dependence between the voltage and the current is not linear, therefore, it affects quality of the waveform, and can negatively affect other consumers connected to the electrical network. Nowadays, a large variety of harmonic, inter-harmonic and subharmonic methods are employed to be synchronized with the fundamental component of the electrical network, thus avoiding the "leakage" effect with techniques based on the Discrete Fourier Transform (DFT), providing an optimum solution for the measurement of harmonics, inter-harmonics and sub-harmonics present in an electrical signal. 2. Wind Energy Conversion System – Squirrel cage induction machine Depending on the electric machine type being used as a generator, there are several configurations for the connection to the electrical network, these configurations have had an evolution over the years. Considering the different types of WECS configuration, the first one is based on a squirrel cage induction machine connected directly to the electrical network via a transformer as shown to Figure 1. The simplest configuration used in this type of wind turbine is the one that used soft-starter [1], [2] for the synchronization and connection to the network, being directly connected to it in a permanent operation. The squirrel cage induction machine is the cheapest machine, simple, robust and with less maintenance, which makes it very interesting for its application in wind generation. However, it has the disadvantage of the reactive power generation capacity, moreover, the American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 232 need to consume reactive power for its magnetization, in no-load and a full-load, which makes it necessary to use capacitor’s bank Figure 1: Asynchronous wind turbine connected to the electrical network This configuration allows the power generation at constant frequency even when the velocity varies between 1 and 2% of the nominal velocity so this configuration is called “fixed-velocity” [3]. The machines with two windings are sometimes used, each with different number of pole pairs to work at two velocities. One of the windings with the lowest synchronization velocity is used for low power and the other with higher synchronization velocity for higher power [4]. Fixed-velocity wind turbines convert wind velocity fluctuations into mechanical fluctuations and these fluctuations in power regardless of the power control system used. In the case of weak networks, these fluctuations can cause voltage fluctuations at the connection point. For this reason, this configuration supposes a low insertion of the wind generation in electrical network, since it does not help to maintain the frequency stability in the network and leaves that work to the conventional synchronous generation. 2.1. Steady-state model of squirrel cage induction machine at fundamental frequency The general configuration of induction machine in the stationary reference-frame is shown below which is obtained of a well-accepted model, with all the parameters referred to the stator, described by the circuit of Figure 2 [5] which shows the model of a three-phase symmetrical induction machine [6]. Squirrel-cage induction machine is used in many industrial applications. Figure 2: Steady-state equivalent circuit of the induction machine The equation that describes the equivalent circuit is given by, 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅 𝐶𝐶𝐶𝐶𝐶𝐶 𝐶𝐶𝑅𝑅𝐶𝐶𝐶𝐶𝑅𝑅𝑅𝑅 𝑆𝑆𝐶𝐶𝑆𝑆𝑅𝑅 𝑇𝑇𝑇𝑇𝑅𝑅𝐶𝐶𝐶𝐶𝑆𝑆𝐶𝐶𝑇𝑇𝐶𝐶𝑅𝑅 𝐸𝐸𝐸𝐸𝑅𝑅𝑅𝑅𝑅𝑅𝑇𝑇𝑅𝑅𝑅𝑅𝑅𝑅𝐸𝐸 𝑁𝑁𝑅𝑅𝑅𝑅𝑁𝑁𝐶𝐶𝑇𝑇𝑁𝑁 𝐺𝐺𝑅𝑅𝑅𝑅𝑇𝑇 𝑊𝑊𝑅𝑅𝐶𝐶𝑊𝑊 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑅𝑅𝐶𝐶𝑅𝑅 𝑆𝑆𝑆𝑆𝑇𝑇𝑅𝑅𝑇𝑇𝑇𝑇𝑅𝑅𝐸𝐸 𝐶𝐶𝑅𝑅𝐶𝐶𝑅𝑅 𝑉𝑉𝑠𝑠 𝑇𝑇𝑠𝑠 𝐼𝐼𝑠𝑠 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿′𝑙𝑙𝑙𝑙 𝐼𝐼′𝑙𝑙 𝑇𝑇′𝑙𝑙 𝐶𝐶� + - American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 233 �𝑉𝑉𝑠𝑠0 � = � 𝑇𝑇𝑠𝑠 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑠𝑠 + 𝐿𝐿𝑀𝑀) 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑇𝑇𝑙𝑙′/𝐶𝐶 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑙𝑙′ + 𝐿𝐿𝑀𝑀)� � 𝐼𝐼𝑠𝑠 𝐼𝐼𝑙𝑙′ � (1) Where 𝑉𝑉𝑠𝑠 and 𝑉𝑉𝑙𝑙′ are the stator and rotor phasor voltages, 𝐼𝐼𝑠𝑠 and 𝐼𝐼𝑙𝑙′ are the stator and rotor phasor currents, respectively, and 𝐿𝐿𝑀𝑀 is the magnetizing inductance given by 𝐿𝐿𝑀𝑀 = 3 2⁄ 𝐿𝐿𝑚𝑚𝑠𝑠 and 𝑗𝑗𝑒𝑒𝑠𝑠 = 2𝜋𝜋𝑆𝑆𝑒𝑒𝑠𝑠 is the electric angular velocity and 𝑆𝑆𝑒𝑒𝑠𝑠 is the frequency of the excitation source in the stator. The expression for the internal electromagnetic torque is, 𝑇𝑇𝑒𝑒 = 𝑉𝑉𝑠𝑠2(𝑇𝑇′𝑙𝑙 𝐶𝐶⁄ ) (𝑇𝑇𝑠𝑠 + 𝑇𝑇′𝑙𝑙 𝐶𝐶⁄ ) + 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑠𝑠 + 𝐿𝐿𝑙𝑙𝑙𝑙) (2) The electromagnetic torque and rotor velocity are related as follows 𝑇𝑇𝑒𝑒 = 𝐽𝐽 2 𝑃𝑃 𝑊𝑊𝑗𝑗𝑒𝑒𝑙𝑙 𝑊𝑊𝑅𝑅 + 𝐵𝐵𝑚𝑚 2 𝑃𝑃 𝑗𝑗𝑒𝑒𝑙𝑙 + 𝑇𝑇𝐿𝐿 (3) where J is the rotor inertia and of the connected load, 𝑗𝑗𝑒𝑒𝑙𝑙 is the electric angular velocity, 𝐵𝐵𝑚𝑚 is the damping coefficient associated with system load and the mechanical load. This parameter is typically small and often negligible, and 𝑇𝑇𝐿𝐿 is the torque load. Additionally, s is the slip defined as 𝐶𝐶 = 𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 𝑗𝑗𝑒𝑒𝑠𝑠 (4) The slip of the induction machine is defined as the relative velocity between the magnetic field produced by the currents injected into the stator and the mechanical velocity of the rotor per unit of field velocity. In general, the field velocity is called the synchronous velocity of the machine and the slip indicates how close the machine is to this velocity. Figure 3 shows the relationship between the rotor velocity and the induced torque generates by the machine. Figure 3: Maximum slip in relation to rotor velocity and induced torque. The slip values describe the operation of the induction machine as follows: the slip is positive if 𝑗𝑗𝑒𝑒𝑠𝑠 > 𝑗𝑗𝑙𝑙 which means that the machine is operating as a motor, the slip is negative if 𝑗𝑗𝑒𝑒𝑠𝑠 < 𝑗𝑗𝑙𝑙 when the machine is operating as a generator, the slip is 1 if 𝑗𝑗𝑙𝑙 = 0 when the machine has the rotor blocked and works as a transformer, and finally the slip is zero if 𝑗𝑗𝑒𝑒𝑠𝑠 = 𝑗𝑗𝑙𝑙 when the machine is rotating but the torque is zero (un- loaded). 0 500 1000 0 5 10 Rotor velocity [rpm] In du ce d To rq ue [N .m ] Pull out Torque Slip (Max) Starting Torque American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 234 If it is a squirrel cage, then the rotor is short-circuited, at that moment a current is induced in the rotor with an angular velocity 𝑗𝑗𝑒𝑒𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 or an electric frequency of 𝑆𝑆𝑒𝑒𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑙𝑙 2𝜋𝜋⁄ . This frequency can be expressed as a function of the slip and the stator frequency: 𝑆𝑆𝑒𝑒𝑙𝑙 = 𝐶𝐶𝑆𝑆𝑒𝑒𝑠𝑠 (4) However, if the rotor is fed with a source of frequency 𝑆𝑆𝑒𝑒𝑙𝑙 it will induce a current in the stator of frequency, 𝑆𝑆𝑒𝑒𝑠𝑠 = 𝑆𝑆𝑒𝑒𝑠𝑠 𝐶𝐶⁄ (5) The phasors current are obtained solving (1), i.e., 𝐼𝐼𝑠𝑠 = |𝐼𝐼𝑠𝑠|∠𝜑𝜑𝑠𝑠 and 𝐼𝐼′𝑙𝑙 = |𝐼𝐼′𝑙𝑙|∠𝜑𝜑𝑙𝑙. Their representations in the time domain in their respective windings are, 𝑅𝑅𝑠𝑠 = |𝐼𝐼𝑠𝑠| cos(𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠) (6) 𝑅𝑅𝑙𝑙 = |𝐼𝐼𝑙𝑙′ |cos (𝐶𝐶𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙′ − 𝜃𝜃𝑒𝑒𝑒𝑒) where 𝜃𝜃𝑒𝑒𝑒𝑒 is the phase angle between the rotor and stator in steady-state. The relation between the angular velocity and the angular displacement is the angle obtained from 𝜃𝜃𝑙𝑙 = ∫ 𝑗𝑗𝑙𝑙(𝜉𝜉)𝑊𝑊𝜉𝜉 + 𝜃𝜃𝑙𝑙(0)𝑡𝑡 0 where 𝜉𝜉 is a dummy variable of integration. In addition, in steady-state results in 𝜃𝜃𝑙𝑙 = 𝑗𝑗𝑙𝑙𝑅𝑅 + 𝜃𝜃𝑒𝑒𝑒𝑒 and then: 𝜃𝜃𝑒𝑒𝑒𝑒 = 𝜃𝜃𝑙𝑙 − 𝑗𝑗𝑙𝑙𝑅𝑅 (7) If the analysis is developed seen from rotor, then the machine inductances are in function of rotor velocity, which describe the machine behavior to be varying in time. For this is necessary a change of variables to the desired reference-frame to avoid the complexity of the differential equations. The general transformation refers the machine variables to a reference-frame, which rotates at any arbitrary angular velocity. Assuming that the machine operates under balanced conditions, the stator and rotor quantities are zero knowing that for steady- state conditions the variables d and q are sinusoidal in all reference-frames except the rotating synchronous reference frame where they are constant. The general configuration of induction machine seen from rotor in the stationary reference-frame is shown in Figure 4 with all the parameters referred to the rotor, which shows the model of a three-phase symmetrical induction machine. Figure 4: Steady-state equivalent circuit induction machine seen from rotor American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 235 The equation that describes the system is given by (as seen from the rotor) �0 𝑉𝑉𝑙𝑙 � = �𝑇𝑇𝑠𝑠 ′/𝐶𝐶𝑙𝑙 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑙𝑙(𝐿𝐿𝑙𝑙𝑠𝑠′ + 𝐿𝐿𝑀𝑀) 𝑗𝑗𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑀𝑀 𝑗𝑗𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑀𝑀 𝑇𝑇𝑙𝑙 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑙𝑙(𝐿𝐿𝑙𝑙𝑙𝑙 + 𝐿𝐿𝑀𝑀)� � 𝐼𝐼𝑠𝑠′ 𝐼𝐼𝑙𝑙 � (8) Since the rotor is moving, the slip seen from the rotor is given by: 𝐶𝐶𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑙𝑙 + 𝑗𝑗𝑙𝑙 𝑗𝑗𝑒𝑒𝑙𝑙 (9) Now, the phasors current are obtained of (8), i.e., 𝐼𝐼𝑙𝑙 = |𝐼𝐼𝑙𝑙|∠𝜑𝜑𝑙𝑙 and 𝐼𝐼′𝑠𝑠 = |𝐼𝐼′𝑠𝑠|∠𝜑𝜑𝑠𝑠. Their representations in the time domain, in their respective windings, are 𝑅𝑅𝑙𝑙 = |𝐼𝐼𝑙𝑙| cos(𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙) (10) 𝑅𝑅𝑠𝑠 = |𝐼𝐼𝑠𝑠′|cos (𝐶𝐶𝑙𝑙𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠′ + 𝜃𝜃𝑒𝑒𝑒𝑒) 2.2. Steady-state model of squirrel cage induction machine at harmonics frequency When the induction machine is fed by non-sinusoidal voltages and currents, the magnetic field in the air-gap and the rotor currents contain components at harmonic frequencies. This harmonics can be positive, negative and zero sequence. The positive sequence harmonics (1, 4, 7, 10, 13, etc.) produce magnetic fields and currents that rotate in the same direction as the fundamental frequency. For the current system (consisting of the fundamental and 5th harmonic components where the 5th harmonic system rotates in counterclockwise direction) the fifth time harmonic has the angular velocity𝑊𝑊𝜃𝜃5 𝑊𝑊𝑅𝑅⁄ = −5𝑗𝑗𝑒𝑒. Phasor representation for forward rotating fundamental and forward rotating 5th harmonic current system is given in Figure 5. Figure 5: Forward (+) rotating fundamental mmf of an induction machine superposed with forward (+) rotating 5th time harmonic rotating in the same direction. 𝑊𝑊𝜃𝜃5 𝑊𝑊𝑅𝑅 = −5𝑗𝑗𝑒𝑒 𝐹𝐹𝑒𝑒 𝑊𝑊𝜃𝜃1 𝑊𝑊𝑅𝑅 = 𝑗𝑗𝑒𝑒 𝐹𝐹5 = 5𝑡𝑡ℎ 𝑇𝑇𝑅𝑅𝑅𝑅𝑁𝑁𝑁𝑁𝑅𝑅𝑇𝑇𝑊𝑊 𝑇𝑇𝐶𝐶𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝐶𝐶𝐶𝐶 Time harmonic 𝐹𝐹5 American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 236 The negative sequence harmonics (2, 5, 8, 11, 14, etc.) develop magnetic fields and currents that rotate in the opposite direction of the fundamental. Zero sequence harmonics (3, 9, 15, 21, etc.) do not develop a useful torque, but cause additional losses in the induction machine. Only wye-connections with insulated neutral and delta are used to connect the stator winding. From this, it follows that in these induction machines with zero sequence there are no homopolar currents in the stator winding. Similar analysis can be performed to determine the rotating directions of each individual time harmonic and to define the positive-, negative-, and zero sequence harmonic orders, as listed in Table 1. Table 1: Positive, negative and zero sequence time harmonics Time-harmonic sequence + − 0 Time-harmonic order 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 ⋮ ⋮ ⋮ Considering that the machine is fed from the stator with a harmonic voltage source while the rotor is short- circuited. This machine has a well-adopted steady-state model [7], with all parameters referred to the stator, described by the circuit of Figure 6. Figure 6: Steady-state harmonic equivalent circuit of the induction machine The equation that then represents the circuit is: �𝑉𝑉𝑠𝑠ℎ 0 � = � 𝑇𝑇𝑠𝑠 + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑠𝑠 + 𝐿𝐿𝑀𝑀) 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑇𝑇𝑙𝑙′/𝐶𝐶ℎ + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑙𝑙′ + 𝐿𝐿𝑀𝑀)� � 𝐼𝐼𝑠𝑠ℎ 𝐼𝐼𝑙𝑙ℎ′ � (11) The equation that determinates the slip of the induction machine at harmonic frequency is, 𝐶𝐶ℎ = ±ℎ𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 ±ℎ𝑗𝑗𝑒𝑒𝑠𝑠 (12) The sign (±) is used for positive or negative sequence, respectively. In this case the harmonics with positive sequence behavior are ℎ = 3𝑁𝑁 + 1 and negative ℎ = 3𝑁𝑁 − 1 for 𝑁𝑁 = 1,2,3, … where the most common harmonics are the 5,7,11,13,15,17,… known as the characteristic harmonics. The harmonic phasors current are obtained of (11), i.e. 𝐼𝐼𝑠𝑠ℎ = |𝐼𝐼𝑠𝑠ℎ|∠𝜑𝜑𝑠𝑠ℎ and 𝐼𝐼𝑙𝑙ℎ′ = |𝐼𝐼𝑙𝑙ℎ′ |∠𝜑𝜑𝑙𝑙ℎ′ . Their representations in the time domain, in their respective windings are, 𝑉𝑉𝑠𝑠ℎ 𝑇𝑇𝑠𝑠 𝐼𝐼𝑠𝑠ℎ 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿′𝑙𝑙𝑙𝑙 𝑇𝑇′𝑙𝑙 𝐶𝐶ℎ� 𝐼𝐼′𝑙𝑙ℎ 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 + American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 237 𝑅𝑅𝑠𝑠ℎ = |𝐼𝐼𝑠𝑠ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ) (13) 𝑅𝑅𝑙𝑙ℎ = |𝐼𝐼𝑙𝑙ℎ′ |cos (𝐶𝐶ℎℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ′ ∓ 𝜃𝜃𝑒𝑒𝑒𝑒) The mathematical analysis developed is implemented in an algorithm in MATLAB/Simulink. To validate the proposed model for harmonic analysis, in steady state, of the squirrel cage induction machine, it is necessary to compare it with the induction machine model in dynamic-state. This same analysis and modeling are developed in the MATLAB/Simulink computer platform. It is necessary to mention that both results are compared once the model in dynamic-state reaches its stable value. Considering that, the induction machine is fed from the stator with a sinusoidal voltage source 𝑉𝑉𝑠𝑠 and the rotor is short-circuited, then the phasor currents in the stator and rotor winding are obtained. The synchronous velocity of the machine and the electromagnetic torque are also determined. The characteristics of the machine are 50 HP/37.3 kW, 4-pole, at 60 Hz and 46.8 A, at a nominal velocity of 1705 rpm. Figure 7: Stator current The parameters of the induction machine are shown in Table 2. The voltage at which the machine is fed is 120 V at 60 Hz. The stator and rotor currents modeled in both states are shown in Figure 7 and 10, respectively. In this case, a mechanical torque of 198 N × m is considered. Table 2: Induction machine parameters Parameters 50 HP/37.3 kW Number of poles 4 Inertia 1.662 kg·m2. Nominal line current 46.8 Amps Nominal line-to-line voltage 460 Vrms Nominal torque 198 N.m Nominal frequency 60 Hz Stator resistance, 𝑇𝑇𝑠𝑠 0.087Ω Stator inductance, 𝐿𝐿𝑙𝑙𝑠𝑠 8 H Rotor resistance, 𝑇𝑇𝑙𝑙 0.034 Ω Rotor inductance, 𝐿𝐿𝑙𝑙𝑙𝑙 8 H Magnetizing inductance, 𝐿𝐿𝑚𝑚𝑙𝑙 = 𝐿𝐿𝑚𝑚𝑠𝑠 155 H Rotor speed 1705 rpm 0 0.2 0.4 0.6 0.8 1 1.2 1.4 -200 0 200 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 238 Figure 8: Zoom of the Figure 11 Figure 8 shows the stator current, in dynamic-state, reaching the stable moment (red line) and comparing with the stator current modelled in steady-state (blue line). In 0.8 seconds, both waveforms of the stator current are compared in magnitude and phase angle and it is observed that they behave in the same way by validating the proposed model in steady-state of the induction machine. The system frequency is 60 Hz, which are observed in Figure 9. Figure 9: Stator frequencies In this figure the frequency of the induction machine, in the dynamic-state (red bar), and steady-state (blue bar) is observed where the similarity in the magnitude of the frequency in both states are analyzed. Additionally, for this study case, the values of the simulation are: slip, 𝐶𝐶 = 0.0464 and the stator power is, 𝑃𝑃𝑠𝑠 = 3,780 𝑊𝑊. This same procedure is performed for the rotor winding where in Figure 10 the rotor current of phase A is shown in dynamic-state, this being compared to the rotor current in steady-state. Both behaviors coincide with magnitude and phase when the current in dynamic-state reaches its stable value. The Figure 11 shows the system frequency in the rotor. Figure 10: Rotor current 0.66 0.68 0.7 0.72 0.74 0.76 0.78 0.8 0.82 -100 0 100 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 0 50 100 150 200 250 300 350 20 40 60 80 100 Frequency [Hz] I st at or [A ] Steady-state Dynamic-state 0.2 0.4 0.6 0.8 1 1.2 1.4 -400 -200 0 200 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 239 Figure 11: Rotor frequencies The electromagnetic torque and the velocity of the induction machine are shown in Figure 12 and Figure 13, respectively and it is observed that the torque that develops in the machine tends to zero in a very habitual behavior when the machine is only fed by the stator. It is also observed that the machine quickly reaches its rated velocity (0.7 sec) because it is operating under normal operating conditions. Figure 12: Torque of the induction machine With the increase in the use of non-linear loads (from power electronics), problems have started to occur, due to the generation of harmonic currents and voltages in the electrical power system. These include: overheating of cables, transformers and induction machines, excessive currents in the neutral, phenomena of resonance between circuit elements, etc. Increasing harmonic voltage distortion can cause malfunctioning of many equipment (especially the less robust ones) that have been designed to operate under normal conditions. Figure 13: Velocity of the induction machine In power systems, induction machines are a very representative component of the load and are widely used in 0 10 20 30 40 50 60 70 80 90 0 50 100 Frequency [Hz] I ro to r [A ] Steady-state Dynamic-state 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -500 0 500 1000 1500 Time [s] To rq ue [N xm ] 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 500 1000 1500 2000 Time [s] V el oc ity [r pm ] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 240 industrial and commercial installations. Induction machines are sensitive to harmonics and are subjected to all kinds of variations of the power source, which affects their characteristics of operation. If it is considered that, the induction machine is fed from the stator with a harmonic voltage source 𝑉𝑉𝑠𝑠ℎ with frequency 𝑆𝑆𝑠𝑠ℎ being the rotor short-circuited. The voltage source is 460 V at 60 Hz in the stator, which includes the third, fifth, and seventh harmonics. The magnitudes and angles of the components of the harmonic voltage are shown in Table 3. Table 3: Harmonic components of the stator voltage Fundamental Third Fifth Seventh Magnitude 460 153.3 92 65.714 Angle 30° 113° 42.85° 137.15° The three-phase harmonic voltage with which it is excited to the stator winding of the induction machine is shown in Figure 14. This voltage waveform is obtained of (11) and the harmonic equivalent circuit of the induction machine shown in Figure 6. Figure 14: Stator voltage Figure 15: Stator current The stator and rotor currents modeled in both states are shown in Figure 15 and Figure 17, respectively. In this case, a mechanical torque of 198 N × m is considered. Figure 16 shows an enlargement of Figure 15 in which the moment is observed where both modeling of the machine coincide in magnitude and angle verifying the 0 0.2 0.4 0.6 0.8 1 -400 -200 0 200 400 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 -500 0 500 1000 Time [s] V ol ta ge [V ] Phase C Phase A Phase B Phase C American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 241 efficiency of the proposed model, in steady-state, for studies of propagation of harmonics. Figure 16: Zoom of the Figure 15 Figure 17 shows the rotor current in dynamic-state. Once this reaches its stable value is compared to the rotor current of the proposed model in steady-state observing that after 1 second both models can be compared in magnitude and phase angle observing the same waveform. Figure 18 shows an enlargement of the Figure 17. Figure 17: Stator current Figure 18: Zoom of the Figure 17 For this study case, the harmonic slip for the positive, negative and zero sequences components, will be equal to: Positive sequence: 𝐶𝐶ℎ(+) = ℎ𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 ℎ𝑗𝑗𝑒𝑒𝑠𝑠 (14) 1.965 1.97 1.975 1.98 1.985 1.99 1.995 2 -100 0 100 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 -200 0 200 400 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 -100 -50 0 50 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 242 Negative sequence: 𝐶𝐶ℎ(−) = ℎ𝑗𝑗𝑒𝑒𝑠𝑠 + 𝑗𝑗𝑙𝑙 ℎ𝑗𝑗𝑒𝑒𝑠𝑠 (15) In according to eqs (14) and (15), the harmonic slip, in positive and negative sequence, will be: 𝐶𝐶ℎ(−) = ℎ𝑗𝑗𝑒𝑒𝑠𝑠 + 𝑗𝑗𝑙𝑙 ℎ𝑗𝑗𝑒𝑒𝑠𝑠 = 5 × 1800 + 1200 5 × 1800 = 𝟏𝟏.𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐶𝐶ℎ(+) = ℎ𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 ℎ𝑗𝑗𝑒𝑒𝑠𝑠 = 7 × 1800 − 1200 7 × 1800 = 𝟎𝟎.𝟗𝟗𝟎𝟎𝟗𝟗𝟗𝟗 The frequencies induced in the winding of the rotor are obtained with: 𝑆𝑆𝑙𝑙 = 𝐶𝐶ℎ × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 (16) Then, the fundamental frequency that the stator winding induces the rotor winding, according to (16), 𝐶𝐶ℎ × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 0.332 × 1 × 377 2𝜋𝜋 = 𝟐𝟐𝟎𝟎 𝑯𝑯𝑯𝑯 The fifth harmonic inducing the stator winding to the rotor winding will be 𝐶𝐶5 × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 1.333 × 5 × 377 2𝜋𝜋 = 𝟏𝟏𝟑𝟑𝟎𝟎 𝑯𝑯𝑯𝑯 Finally, the seventh harmonic, inducing the stator winding to the rotor winding, will be 𝐶𝐶7 × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 0.9047 × 7 × 377 2𝜋𝜋 = 𝟏𝟏𝟑𝟑𝟎𝟎 𝑯𝑯𝑯𝑯 These frequencies can be observed in Figure 19 and Figure 20 in which are shown the induced frequencies of the stator to the rotor, respectively, in according with the previous analysis. Figure 19: Stator frequencies 0 50 100 150 200 250 300 350 400 450 500 50 100 150 Frequency [Hz] I st at or [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 243 Figure 20: Rotor frequencies These frequencies induced in the rotor cannot be called harmonic frequencies since they are not integral multiples of the fundamental frequency but it should be noted that they are the 17th and 19th harmonic of the fundamental frequency of the rotor. The electromagnetic torque and the velocity of the induction machine under harmonics operating conditions are shown in Figure 21 and Figure 22, respectively and it is observed that the torque that develops in the machine present harmonic content when the machine is fed by the harmonic voltage. It is also observed that the machine it takes time to reach its rated velocity in this conditions operating. Figure 21: Electromagnetic torque of the induction machine in harmonic conditions of operation. Figure 22: Velocity of the induction machine 3. Wind Energy Conversion System – Doubly-fed induction machine In this configuration, known as a double-fed wind turbine, the rotor is connected to the electrical network via a 0 50 100 150 200 250 300 350 400 450 0 20 40 60 80 Frequency [Hz] I ro to r [A ] Steady-state Dynamic-state 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -500 0 500 1000 1500 Time [s] To rq ue [N xm ] 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 500 1000 1500 Time [s] V el oc ity [r pm ] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 244 frequency converter while the stator is connected directly to the electrical network via a transformer as shown in Figure 23. With this configuration, it is possible to control the active and reactive power [8], with low levels of harmonics both on the generator and on the electrical network side. The performance is limited by the size of the equipment. For a power equipment of 25% of the rated power of the machine, the speed variation is around ± 33% [9]. In addition, it has the advantage of delivering over-rated power to the electrical network at super- synchronous operating velocities [8]. The size required for the converter in this configuration makes it very economically attractive. Its main disadvantage is the use of brushes and the need for additional protection in case of faults in the electrical network. A dip of voltage in the network causes an increase in stator current, and due to the magnetic coupling between stator and rotor windings, this current will also flow through the rotor and the power equipment can destroy it [10]. This requires additional equipment such as the so-called "active crowbar”. In the technology, a double-fed induction generator with wound rotor (DFIG) is used. Although its axis rotates at variable velocity, the DFIG can generate a constant frequency voltage when the rotor windings are fed with a frequency converter. This converter consists of two reversible AC/DC electronic converters connected to each other via a DC bus. The converter connected to the rotor regulates the amplitude, frequency and phase sequence of the voltage applied to the rotor, allowing a vector control of the machine to regulate the electromagnetic torque and the power factor of the generator over a wide velocity range of rotation [11], [12]. Figure 23: Asynchronous wind turbine connected to the electrical network 3.1. Steady-state model of doubly-fed induction machine at fundamental frequency In DFIG, the rotor is connected to the electrical network via a frequency converter and the stator is connected directly to the network. This topology can generate a stator voltage 𝑉𝑉𝑠𝑠 of angular frequency 𝑗𝑗𝑒𝑒𝑠𝑠 is constant, whereas by the rotor the voltage 𝑉𝑉𝑙𝑙 of angular frequency 𝑗𝑗𝑙𝑙 variable. At velocity 𝑗𝑗𝑙𝑙 is known as relative 𝐺𝐺𝑅𝑅𝐶𝐶𝑅𝑅𝑇𝑇𝑅𝑅𝑅𝑅𝐶𝐶𝑇𝑇 𝑇𝑇𝑇𝑇𝑅𝑅𝐶𝐶𝐶𝐶𝑆𝑆𝐶𝐶𝑇𝑇𝐶𝐶𝑅𝑅𝑇𝑇 𝐸𝐸𝐸𝐸𝑅𝑅𝑅𝑅𝑅𝑅𝑇𝑇𝑅𝑅𝑅𝑅𝑅𝑅𝐸𝐸 𝑁𝑁𝑅𝑅𝑅𝑅𝑁𝑁𝐶𝐶𝑇𝑇𝑁𝑁 𝐺𝐺𝑅𝑅𝑅𝑅𝑇𝑇 − 𝑇𝑇𝐶𝐶𝑏𝑏 𝑊𝑊𝑅𝑅𝐶𝐶𝑊𝑊 𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑅𝑅𝐶𝐶𝑅𝑅 𝐼𝐼𝐶𝐶𝑊𝑊𝑇𝑇𝑅𝑅𝑅𝑅𝑅𝑅𝐶𝐶𝐶𝐶 Machine 𝐹𝐹𝑇𝑇𝑅𝑅𝑆𝑆𝑇𝑇𝑅𝑅𝐶𝐶𝑅𝑅𝐹𝐹 𝑅𝑅𝐶𝐶𝐶𝐶𝑅𝑅𝑅𝑅𝑇𝑇𝑅𝑅𝑅𝑅𝑇𝑇 American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 245 velocity or slip. According to the fundamental principle of rotating electric machines with respect to the angular frequencies of magnetic fields between stator and rotor windings, to develop a pair with a value other than zero [13], and considering that the number of pole pairs P is same in both windings, it is, 𝑗𝑗 = 𝑗𝑗𝑒𝑒𝑠𝑠 ± 𝑗𝑗𝑙𝑙 (17) Since 𝑗𝑗 is the electric angular velocity of the rotor, which is defined as, 𝑗𝑗 = 𝑃𝑃 ± 𝑗𝑗𝑚𝑚 (18) where 𝑗𝑗𝑚𝑚 is the mechanical angular velocity or rotation of the rotor. The velocity 𝑗𝑗 may be greater or less than𝑗𝑗𝑒𝑒𝑠𝑠. In the case where 𝑗𝑗 < 𝑗𝑗𝑒𝑒𝑠𝑠, the velocity 𝑗𝑗𝑙𝑙 is positive in (17), and the DFIG is capable to operate at sub-synchronous velocity. On the other hand, if in 𝑗𝑗 > 𝑗𝑗𝑙𝑙 the velocity 𝑗𝑗𝑙𝑙 is negative then the DFIG will operate at super-synchronous velocity. A general configuration of a DFIG uses a wound rotor with slip rings to transmit current between the stator and rotor voltage sources. This machine has a well-adopted steady-state model, with all the parameters referred to the stator described by Figure 24. The voltage equation that describes the equivalent circuit is given by: � 𝑉𝑉𝑠𝑠 1 𝐶𝐶 𝑉𝑉′𝑙𝑙 � = � 𝑇𝑇𝑠𝑠 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑠𝑠 + 𝐿𝐿𝑀𝑀) 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 𝑇𝑇𝑙𝑙′/𝐶𝐶 + 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠(𝐿𝐿𝑙𝑙𝑙𝑙′ + 𝐿𝐿𝑀𝑀)� � 𝐼𝐼𝑠𝑠 𝐼𝐼𝑙𝑙′ � (19) Figure 24: Equivalent circuit in steady-state of the induction machine The parameters and variables of this equivalent circuit were described in the previous section. Now, if the rotor is short-circuited, a current is induced in the rotor at an angular frequency 𝑗𝑗𝑒𝑒𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑠𝑠 − 𝑗𝑗𝑙𝑙 or at an electric frequency of 𝑆𝑆𝑒𝑒𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑙𝑙/2𝜋𝜋. This is important when the induction machine is double fed, for example if a machine has a slip of 0.2 and is fed from the stator at a frequency of 60 Hz, it induces a frequency of 12 Hz in the rotor and if the rotor is also fed with a frequency of 7 Hz, it will also induce a frequency of 35 Hz in the stator. Therefore, the stator will have a current with two components, one at 60 Hz and other at 35 Hz and the rotor will have one of 6 Hz and other of 12 Hz. It is common for a doubly fed induction machine with a constant frequency in the rotor, which will induce different frequencies in the stator depending on the slip values. These frequencies are called inter-harmonic, since they are not integer multiples of the fundamental frequency 𝑆𝑆𝑒𝑒𝑠𝑠. Considering that, the doubly fed induction machine is fed from the stator with a voltage source 𝑉𝑉𝑠𝑠 of frequency 𝑆𝑆𝑒𝑒𝑠𝑠 with the short-circuited rotor. The phasors current, in their representations in the time domain, for their respective windings are [14], 𝑉𝑉𝑠𝑠 𝑇𝑇𝑠𝑠 𝐼𝐼𝑠𝑠 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿′𝑙𝑙𝑙𝑙 𝑇𝑇′𝑙𝑙 𝐶𝐶� 𝑉𝑉′𝑙𝑙 𝐶𝐶 𝐼𝐼′𝑙𝑙 𝑗𝑗𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 + - + - American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 246 𝑅𝑅𝑠𝑠 = |𝐼𝐼𝑠𝑠| 𝑅𝑅𝐶𝐶𝐶𝐶(𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠) (20) 𝑅𝑅𝑙𝑙 = |𝐼𝐼𝑙𝑙′|𝑅𝑅𝐶𝐶𝐶𝐶 (𝐶𝐶𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙′ − 𝜃𝜃𝑒𝑒𝑒𝑒) Considering the opposite case, the DFIG is fed from the rotor with a voltage source 𝑉𝑉𝑙𝑙 of frequency 𝑆𝑆𝑒𝑒𝑙𝑙 and stator in short-circuit. If the rotor is in motion, the slip seen from the rotor is given by, 𝐶𝐶𝑙𝑙ℎ = ±ℎ𝑗𝑗𝑒𝑒𝑙𝑙 + 𝑗𝑗𝑙𝑙 ±ℎ𝑗𝑗𝑒𝑒𝑙𝑙 (21) The sign (+) is used for positive and �– � for negative sequence, respectively. A simple analysis to verify this, for example, if the rotor is fed with a frequency 𝑗𝑗𝑒𝑒𝑙𝑙 = 𝐶𝐶𝑗𝑗𝑒𝑒𝑠𝑠 from where 𝑗𝑗𝑙𝑙 = 𝑗𝑗𝑒𝑒𝑠𝑠(1 − 𝐶𝐶), then from this results in 𝐶𝐶𝑙𝑙 = 1/𝐶𝐶, that is, the induced frequency in the stator is 𝑆𝑆𝑒𝑒𝑠𝑠 = 𝐶𝐶𝑙𝑙𝑆𝑆𝑒𝑒𝑙𝑙 = 𝑆𝑆𝑒𝑒𝑙𝑙/𝐶𝐶 which is in agreement with (5). The phasors current in their representations in the time domain, in their respective windings, are 𝑅𝑅𝑙𝑙 = |𝐼𝐼𝑙𝑙| cos(𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙) (22) 𝑅𝑅𝑠𝑠 = |𝐼𝐼𝑠𝑠′|cos (𝐶𝐶𝑙𝑙𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠′ + 𝜃𝜃𝑒𝑒𝑒𝑒) The complete solution for the DFIG is given by the sum of the two effects, i.e., 𝑅𝑅𝑠𝑠 = |𝐼𝐼𝑠𝑠| cos(𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠) + |𝐼𝐼𝑠𝑠′| cos�𝐶𝐶𝑙𝑙𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠′ + 𝜃𝜃𝑒𝑒𝑒𝑒� (23) 𝑅𝑅𝑙𝑙 = |𝐼𝐼𝑙𝑙| cos(𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙) + |𝐼𝐼𝑙𝑙′ | cos�𝐶𝐶𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙′ − 𝜃𝜃𝑒𝑒𝑒𝑒� It should be noted that both currents are sinusoidal only if the voltage source in the rotor has a frequency 𝑆𝑆𝑒𝑒𝑙𝑙 = 𝐶𝐶𝑆𝑆𝑒𝑒𝑠𝑠. 3.2. Steady-state model of doubly-fed induction machine at harmonic frequency In an doubly-fed induction machine, the harmonics are generated on both windings: harmonics in the stator voltage source with frequencies 𝒇𝒇𝒔𝒔𝒔𝒔 = 𝒔𝒔𝒇𝒇𝒆𝒆𝒔𝒔 and harmonics in the rotor voltage source with frequencies 𝒇𝒇𝒓𝒓𝒔𝒔 = 𝒔𝒔𝒇𝒇𝒆𝒆𝒓𝒓, where h is an integer. The induced frequencies in the rotor, however, due to harmonics in the stator, are inter- harmonics of the rotor fundamental frequency and vice versa. Figure 25: Harmonic equivalent circuit in steady-state of DFIG seen from the rotor 𝑉𝑉𝑠𝑠ℎ 𝑇𝑇𝑠𝑠 𝐼𝐼𝑠𝑠ℎ 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿′𝑙𝑙𝑙𝑙 𝑇𝑇′𝑙𝑙 𝐶𝐶ℎ� 𝐼𝐼′𝑙𝑙ℎ 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑀𝑀 + - American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 247 Considering that the machine is fed from the stator with a source of harmonic voltage of frequency ℎ𝑆𝑆𝑒𝑒𝑠𝑠, being the rotor in short-circuit. This machine has a well-adopted steady-state model, with all parameters referred to the stator, described by the circuit of Figure 25. The equation that represents the circuit is, � 0 𝑉𝑉𝑙𝑙ℎ � = � 𝑇𝑇𝑠𝑠′ 𝐶𝐶𝑙𝑙ℎ + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙(𝐿𝐿𝑙𝑙𝑠𝑠′ + 𝐿𝐿𝑀𝑀) 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑀𝑀 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑀𝑀 𝑇𝑇𝑙𝑙 + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙(𝐿𝐿𝑙𝑙𝑙𝑙 + 𝐿𝐿𝑀𝑀) � �𝐼𝐼𝑠𝑠ℎ ′ 𝐼𝐼𝑙𝑙ℎ � (24) The phasors current in their representations in the time domain, in their respective windings, are 𝑅𝑅𝑙𝑙ℎ = |𝐼𝐼𝑙𝑙ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ) (25) 𝑅𝑅𝑠𝑠ℎ = |𝐼𝐼𝑠𝑠ℎ′ |cos (𝐶𝐶𝑙𝑙ℎℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ′ ± 𝜃𝜃𝑒𝑒𝑒𝑒) The complete solution for the DFIG is given by the sum of the two effects, i.e., 𝑅𝑅𝑠𝑠ℎ = |𝐼𝐼𝑠𝑠ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ) + |𝐼𝐼𝑠𝑠ℎ′ | cos�𝐶𝐶𝑙𝑙ℎℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ′ ± 𝜃𝜃𝑒𝑒𝑒𝑒� (26) 𝑅𝑅𝑙𝑙ℎ = |𝐼𝐼𝑙𝑙ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ) + |𝐼𝐼𝑙𝑙ℎ′ | cos�𝐶𝐶ℎℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ′ ∓ 𝜃𝜃𝑒𝑒𝑒𝑒� The slip’s equation, (12) and (21), are used are only valid under balanced conditions, for positive and negative sequences. To an analysis of zero sequence, the induction machine works as two-decoupled windings, resulting in the circuit of Figure 26, Figure 26: Induction machine triplex harmonic model seen from stator and rotor The equation that represents the circuit is 𝑉𝑉𝑠𝑠ℎ = (𝑇𝑇𝑠𝑠 + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠)𝐼𝐼𝑠𝑠ℎ (27) 𝑉𝑉𝑙𝑙ℎ = (𝑇𝑇𝑙𝑙 + 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑙𝑙𝑙𝑙)𝐼𝐼𝑙𝑙ℎ The phasors current are obtained of (27) in their representations in the time domain, in their respective windings, are 𝑅𝑅𝑠𝑠ℎ = |𝐼𝐼𝑠𝑠ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ) (28) 𝑅𝑅𝑙𝑙ℎ = |𝐼𝐼𝑙𝑙ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ) 𝐼𝐼𝑙𝑙ℎ 𝑉𝑉𝑠𝑠ℎ 𝑇𝑇𝑠𝑠 𝐼𝐼𝑠𝑠ℎ 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝐿𝐿𝑙𝑙𝑠𝑠 𝑗𝑗ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝐿𝐿𝑙𝑙𝑙𝑙 𝑇𝑇𝑙𝑙 𝑉𝑉𝑙𝑙ℎ + + - - American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 248 Now the general solution for balanced conditions including voltage sources at fundamental and harmonic frequencies in the stator and rotor is given by 𝑅𝑅𝑠𝑠 = ∑ |𝐼𝐼𝑠𝑠ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ) + ∑ |𝐼𝐼𝑠𝑠ℎ′ | cos�𝐶𝐶𝑙𝑙ℎℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑠𝑠ℎ′ + 𝜃𝜃𝑒𝑒𝑒𝑒�𝐻𝐻 ℎ=1,3𝑘𝑘+1 𝐻𝐻 ℎ=1 + ∑ |𝐼𝐼𝑠𝑠ℎ′ | cos�𝐶𝐶𝑙𝑙ℎℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 +𝐻𝐻 ℎ=3𝑘𝑘−1 𝜑𝜑𝑠𝑠ℎ′ − 𝜃𝜃𝑒𝑒𝑒𝑒� (29) 𝑅𝑅𝑙𝑙 = ∑ |𝐼𝐼𝑙𝑙ℎ| cos(ℎ𝑗𝑗𝑒𝑒𝑙𝑙𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ) + ∑ |𝐼𝐼𝑙𝑙ℎ′ | cos�𝐶𝐶ℎℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 + 𝜑𝜑𝑙𝑙ℎ′ − 𝜃𝜃𝑒𝑒𝑒𝑒�𝐻𝐻 ℎ=1,3𝑘𝑘+1 𝐻𝐻 ℎ=1 + ∑ |𝐼𝐼𝑙𝑙ℎ′ | cos�𝐶𝐶ℎℎ𝑗𝑗𝑒𝑒𝑠𝑠𝑅𝑅 +𝐻𝐻 ℎ=3𝑘𝑘−1 𝜑𝜑𝑙𝑙ℎ′ + 𝜃𝜃𝑒𝑒𝑒𝑒� (30) These equations are given for the proposed doubly-fed induction machine model. Below is a brief explanation: all the current harmonics included positive-, negative-, and zero-sequence harmonics (first summation of 29). All the current harmonics because of the induction effect of the positive-sequence voltage source harmonics in the rotor (second summation of 29) and all the current harmonics because of induction effect of the negative- sequence voltage source harmonics in the rotor (third summation of 29). The same interpretation is for (30). Considering that, the induction machine is fed from the stator with a sinusoidal voltage source 𝑉𝑉𝑠𝑠 and the rotor is fed from the rotor with a sinusoidal voltage source 𝑉𝑉𝑙𝑙 too. The characteristics of the machine are: 50 HP/37.3 kW, 4-pole, at 60 Hz and 46.8 A, at a nominal velocity of 1705 rpm. The parameters of the induction machine are shown in Table 2. The voltage at which the machine is fed is 460 V at 60 Hz and the rotor voltage is 50 V. The stator and rotor currents modeled in both states are shown in Figure 27 and Figure 29, respectively, where they show the behavior of the stator current, in dynamic-state and steady-state. At the end, both waveforms of the stator current are compared in magnitude and phase angle and it is observed that they behave in the same way by validating the proposed model in steady-state of the induction machine. A mechanical torque of 198 N × m is also used for this case study. Figure 27: Stator current Figure 28: Zoom of the Figure 27 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 -100 0 100 200 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 0.68 0.7 0.72 0.74 0.76 0.78 0.8 -50 0 50 100 Time [s] St at or c ur re nt [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 249 Figure 29: Rotor current Figure 30: Zoom of the Figure 33 Figure 31 and Figure 32 show the induced frequencies of the stator to the rotor and vice versa. According to the analysis previously demonstrated, the machine is excited, on the stator side, to a voltage 𝑉𝑉𝑠𝑠 at 60 Hz. Figure 31: Stator frequency Considering that the machine has a slip, 𝐶𝐶 = 0.1, then a current component is generated at 6 Hz, as 𝑆𝑆𝑙𝑙 = 𝐶𝐶𝑆𝑆𝑒𝑒𝑠𝑠 = 0.1 × 60 = 6 𝐻𝐻𝐻𝐻 verifying that both, stator and rotor, generate current with two components each: the own frequency and the induced frequency of the other winding. The electromagnetic torque and the velocity of the induction machine under normal operating conditions are shown in Figure 33 and Figure 34, respectively. Figure 32: Rotor frequency 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 -100 0 100 200 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state 0.5 0.6 0.7 0.8 0.9 1 1.1 -50 0 50 100 150 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state 0 50 100 150 200 250 300 20 40 60 80 100 120 0 Frequency [Hz] I st at or [A ] Steady-state Dynamic-state 0 5 10 15 20 25 30 35 40 45 50 20 40 60 80 100 Frequency [Hz] I ro to r [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 250 Figure 33: Electromagnetic torque of the induction machine Figure 34: Velocity of the induction machine Considering that the induction machine is fed from the stator with a harmonic voltage source 𝑉𝑉𝑠𝑠ℎ with frequency 𝑆𝑆𝑠𝑠ℎ and in the rotor a harmonic voltage source 𝑉𝑉𝑙𝑙ℎ with frequency 𝑆𝑆𝑙𝑙ℎ. The voltage source is 460 V at 60 Hz in the stator, which includes the third, fifth, and seventh harmonics. The voltage source in rotor is 50 V at 45 Hz, which includes the third, fifth, and seventh harmonics. The magnitudes and angles of the components of the harmonic voltages are shown in Table 4. Table 4: Harmonic components of the stator and rotor voltage STATOR Fundamental Third Fifth Seventh Magnitude 460 153.3 92 65.714 Angle 30° 113° 42.85° 137.15° ROTOR Fundamental Third Fifth Seventh Magnitude 50 16.7 10 7.142 Angle 13.6° 8.63° 114.5° 47.32° The stator and rotor currents modeled in both states are shown in Figure 35 and Figure 37, respectively. A mechanical torque of 198 N × m is also used for this case study. 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -500 0 500 1000 1500 Time [s] To rq ue [N xm ] 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 500 1000 1500 Time [s] V el oc ity [r pm ] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 251 Figure 35: Stator current Figure 36: Zoom of the Figure 35 Figure 37: Rotor current Figure 38: Zoom of the Figure 37 Considering a slip value of 0.96, the fundamental frequency that the stator winding induces the rotor winding, according to (16), 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 -100 0 100 200 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 1.95 1.955 1.96 1.965 1.97 1.975 1.98 1.985 1.99 1.995 2 -50 0 50 100 150 Time [s] S ta to r c ur re nt [A ] Steady-state Dynamic-state 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -400 -200 0 200 400 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 -50 0 50 100 Time [s] R ot or c ur re nt [A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 252 𝐶𝐶ℎ × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 0.96 × 1 × 377 2𝜋𝜋 = 𝟓𝟓𝟗𝟗.𝟑𝟑𝟎𝟎 𝑯𝑯𝑯𝑯 The fifth harmonic inducing the stator winding to the rotor winding will be 𝐶𝐶5 × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 1.15 × 5 × 377 2𝜋𝜋 = 𝟏𝟏𝟗𝟗𝟓𝟓.𝟏𝟏𝟑𝟑 𝑯𝑯𝑯𝑯 Finally, the seventh harmonic, inducing the stator winding to the rotor winding, will be 𝐶𝐶7 × ℎ × 𝑗𝑗𝑒𝑒𝑠𝑠 2𝜋𝜋 = 0.86 × 7 × 377 2𝜋𝜋 = 𝟏𝟏𝟑𝟑𝟏𝟏.𝟗𝟗 𝑯𝑯𝑯𝑯 These frequencies can be observed in Figure 39 and Figure 40 in which are shown the induced frequencies of the stator to the rotor, respectively, in according with the previous analysis. Figure 39: Stator frequencies Figure 40: Rotor frequencies Table V summarizes the harmonic currents in the induction machine for the case studies. Note that the waveform current presented in the paper has been attained from current shown in this table, which have been obtained from the solution of the equations mentioned in the previous section. 0 50 100 150 200 250 300 350 400 450 500 20 40 60 80 100 Frequency [Hz] I st at or [A ] Steady-state Dynamic-state 0 50 100 150 200 250 300 350 400 450 20 40 60 80 Frequency [Hz] I ro to r [ A ] Steady-state Dynamic-state American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 253 Table 5: Summary of the harmonic currents for this case study STATOR ROTOR Sequence Magnitude Angle Frequency Sequence Magnitude Angle Frequency + 460 180.6° 60 Hz + 45 204.9° 45 Hz 0 153.3 183.2° 180 Hz 0 15 -78.3° 193.5 Hz − 260 181.2° 300 Hz − 9 -85.1° 607.8 Hz + 71.42 179.2° 420 Hz + 6.42 -28.4° 1122 Hz − 41.81 18.2° 610.2 Hz − 2.65 120.2° 457.65 Hz + 35.38 -5.8° 267 Hz + 2.36 117° 200.25 Hz 4. Conclusion This research proposed a steady-state model for the squirrel cage induction machine and doubly-fed induction machine for harmonic and inter-harmonic studies. The currents response of the different wind energy conversion system are obtained of the mathematical model when it is fed with sinusoidal and nonsinusoidal voltage sources. Depending on the slip and the fundamental frequency of both voltage sources the current harmonics and non- harmonics can exist on both sides of the machine. The simulation results obtained with the steady-state proposed model were compared with those obtained from dynamic simulation, where the comparison between them is in complete agreement. Finally, this research gives a clear analysis of the frequencies generated by the wind energy conversion system, in steady-state, resulting in a clear and precise model for harmonic and inter- harmonic analysis, which can be used for “harmonic” analysis in a power system. References [1] T. Ackermann. “Wind power in power systems”, Wiley, 2005. [2] J. Martínez García, M. García-García, M. P. Comech, D. García-García, “Modelling and simulation of an asynchronous wind turbine of squirrel cage”, ICREPQ’04 Proceeding, 31 March – 02 April 2004, Barcelona. [3] J. G. Slootweg, W. L. Kling, “Is the answer blowing in the wind?” IEEE Power & Energy Magazine, pp 26-33, Nov/Dec 2003. [4] R. Datta, V. T. Ranganathan, “Variable-speed wind power generation using doubly-fed wind rotor induction machine a comparison with alternative schmes.”, IEEE Transactions on Energy Conversion, vol. 17, no. 3, p.p. 414-421, Sept. 2002. [5] P. C. Krause, Analysis of Electric Machinery. New York, NY, USA: McGraw-Hill, 1987. [6] E. H. Camm, M. R. Behnke, O. Bolado, M. Bollen, M. Bradt, C. Brooks, W. Dilling, M. Edds, W. J. Hejdak, D. Houseman, S. Klien, F. Li, J. Li, P. Maibach, T. Nicolai, J. Patino, S. V. Pasupulati, N. Samaan, S. Saylors, T. Seibert, T. Smith,M. Starke, and R.Walling,“Wind power plant grounding, overvoltage protection and insulation coordination: IEEE PES Wind Plant Collector System Design American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2017) Volume 33, No 1, pp 230-254 254 Working Group,” in Proc. IEEE Power and Energy Soc. General Meet., 2009, pp. 1–8. [7] L. Fan, S. Yuvurajan and R. Kavasseri, “Harmonic Analysis of a DFIG for a Wind Energy Conversion System” IEEE Trans. Energy Convers., vol. 25, no. 1, pp. 181-190, March. 2010. [8] A. Tapia, G. Tapia, J.S. Ostolaza, J.R. Saenz, .“Modeling and control of a wind turbine driven doubly fed induction generator.”, IEEE Transactions on Energy Conversion, vol. 18, no. 2, pp. 194-204, June 2003. [9] J. Morren, S.W.H. Haan, .“Ride through of wind turbines with doubly fed induction generator during a voltage dip.”, IEEE Transactions on Energy Conversion, vol. 20, no. 2, pp. 435-441, June 2005. [10] I. Erlich, W. Winter, A. Dittrich,“Advanced Grid Requirements for the Integration of Wind Turbines into the German transmission system.”, 2006 IEEE Power Engineering Society General Meeting, 18-22 June 2006, pp. 7. [11] Peña R., Clare J.C. y Ascher G.M., “Doubly Fed Induction Generator using back to back PWM converters and its application to variable-speed wind energy generation”, IEE Proceeding Electric Power Applications, Vol. 43, No. 3, pp. 231-241, Mar 1996. [12] Müller S., Deicke M. y De Doncker R., “Doubly Fed Induction Generator Systems for Wind Turbines”, IEEE Industry Applications Magazine, pp. 26-33, May-June 2002. [13] Kundur Prabha, Power System Stability and Control. New York: McGraw-Hill, 1994. [14] Krause Paul, Wasynczuk Oleg and Sudhoff Scott, Analysis of Electric Machinery. New York: IEEE Press, 1995