Optimal switching angles family with genetic algorithms
Vector control has been widely used for the high-performance drive of the induction motor. As in DC motor, torque control of the induction motor is achieved by controlling torque and flux components independently. Vector control techniques can be separated into two categories: direct and indirect flux vector orientation control schemes. For direct control methods, the flux vector is obtained by using stator terminal quantities, while indirect methods use the machine slip frequency to achieve field orientation.
The overall performance of field-oriented-controlled induction motor drive systems is directly related to the performance of current control. Therefore, decoupling the control scheme is required by compensation of the coupling effect between q-axis and d-axis current dynamics (Jung et al., 1999; Lin et al., 2000; Suwankawin et al., 2002).
The PWM is the interface between the control block of the electrical drive and its associated electrical motor (fig.1). This function controls the voltage or the current inverter (VSI or CSI) of the drive. The performance of the system is influenced by the PWM that becomes therefore an essential element of the system. A few problems of our days concerning the variable speed system are related to the conventional PWM: inverter switching losses, acoustical noise, and voltages harmonics (fig.2).
Harmonic elimination and control in inverter applications have been researched since the early 1960’s (Bouchhida et al., 2007, 2008; Czarkowski et al., 2002; García et al., 2003; Meghriche et al., 2004, 2005; Villarreal-Ortiz et al., 2005; Wells et al., 2004). The majority of these papers consider the harmonic elimination problem in the context of either a balanced connected load or a single phase inverter application. Typically, many papers have focused on finding solutions and have given little attention to which solution is optimal in an application context.
A pre-calculated PWM approach has been developed to minimize the harmonic ratio within the inverter output voltage (Bouchhida et al., 2007, 2008; Bouchhida, 2008, 2011). Several other techniques were proposed in order to reduce harmonic currents and voltages. Some benefit of harmonic reduction is a decrease of eddy currents and hysterisis losses. That increase of the life span of the machine winding insulation. The proposed approach is integrated within different control strategies of induction machine.
A novel harmonic elimination pulse width modulated (PWM) strategy for three-phase inverter is presented in this chapter. The torque ripple of the induction motor can be significantly reduced by the new PWM technique. The three-phase inverter is associated with a passive LC filter. The commutation angles are predetermined off-line and stored in the microcontroller memory in order to speed up the online control of the induction motors. Pre-calculated switching is modelled to cancel the greater part of low-order harmonics and to keep a single-pole DC voltage across the polarized capacitors. A passive LC filter is designed to cancel the high-order harmonics. This approach allows substantial reduction of the harmonic ratio in the AC main output voltage without increasing the number of switches per period. Consequently, the duties of the semiconductor power switches are alleviated. The effectiveness of the new harmonic elimination PWM technique for reducing torque-ripple in inverter-fed induction motors is confirmed by simulation results. To show the validity of our approach, DSP-based experimental results are presented.
2. New three-phase inverter model
Figure 3 shows the new structure of the three-phase inverter, with E being the dc input voltage and U12-out=VC1-VC2, U23-out=VC2-VC3 and U31-out=VC3-VC1 are the ac output voltage obtained via a three LC filter. R represents the internal inductors resistance. Qi and Q’i (i=1,2,3) are the semiconductor switches. It is worth to mention that transistors Qi and Q’i undergo complementary switching states. VC1, VC2 and VC3 are the inverter filtered output voltages taken across capacitors C1, C2 and C3 respectively.
2.1. Harmonic analysis
In ideal case, the non filtered three inverter output voltage Vo1, Vo2 and Vo3 is desired to be:
With: α=ωt, and ω is the angular frequency.
The relative Fourier harmonic coefficients of (1), with respect to E, are given by (2.1) or more explicitly by (2.2).
With index i is the phase number.
However, in practice, the non filtered inverter output voltage Vo1 (Vo2, Vo3) is a series of positive impulses (see Fig. 4): 0 when Q’1 (Q’2, Q’3) is on and E when Q’1(Q’2, Q’3) is off, so that voltage of capacitor C1 (C2, C3) is always a null or a positive value. In this case, the relative Fourier harmonic coefficients, with respect to E, are given by (3).
k is the harmonic order
αi are the switching angles
Nα is the number of αi per half period
The other inverter outputs Vo2 and Vo3 are obtained by phase shifting Vo1 with 2/3 π, 4/3 π, respectively as illustrated in fig. 4 for Nα=5.
The objective is to determine the switching angles αi so as to obtain the best possible match between the inverter output Vo1 and Vo1-ideal.
For this purpose, we have to compare their respective harmonics. A perfect matching is achieved only when an infinite number of harmonics is considered as given by (4)., for
In practice, the number of harmonics N that can be identical is limited. Thus, a nonlinear system of N+1 equations having Nα unknowns is obtained as:, for
|Family||symbol||Angles (radians)||Angles (degrees)|
2.2. Dynamic LC filter behavior
Considering the inverter direct output fundamental, the LC filter transfer function is given by:
From (6), one can notice that for ω=0, =1, meaning that the mean value (dc part) of the input voltage is not altered by the filter. Consequently, the inverter dc output part is entirely transferred to capacitor C1. The same conclusion can be drawn for capacitor C2 and C3.
Letting and, the filter transfer function can rewritten as:
For a given harmonic component of order k, the LC filter transfer function Tk is obtained by replacing ω with kω as:
Assuming that the filter L and C components are not saturated, and using the superposition principle, we obtain the inverter filtered output voltages VC1, VC2 and VC3, taken across capacitors C1, C2 and C3 as given by (9.1), (9.2) and (9.3) respectively:
Tk and φk are the kth order magnitude and phase components of the LC filter transfer function respectively.
Each k harmonic term of the inverter output voltage has a frequency of kω and an amplitude equal to akTk, Where ak is the amplitude of the kth order harmonic of Vo1.
The transfer function magnitude and phase are given, respectively, by:
The maximum value Tmax of Tk, obtained when, can be expressed by:
In which case, (12) corresponds to a maximum angular frequency, this last is given by:
The maximum angular frequency ωmax exists if, and only if,. The filter transfer function will exhibit a peak value then decreases towards zero. As consequence, the fundamental, as well as the harmonics, are amplified, this leads to undesirable situation, as illustrated in fig.5.
If (14) is satisfied, the filter transfer function will exhibit a damped behaviour.
This condition matches both practical convenience and system objectives.
2.3. Harmonic rate calculation
Using (1) to (3), the non filtered inverter output voltages Vo1, Vo2 and Vo3 can be expressed as
Taking into consideration the filter transfer function, we get the expressions (16), (17) and (18) for VC1, VC2 and VC3 respectively.
Using (15), (16), (17) and (18), we get the three inverter filtered output voltages expressions as:
These equations are put in following matric form:
with : ;;
3. Indirect field- oriented induction motor drive
The dynamic electrical equations of the induction machine can be expressed in the d-q synchronous reference frame as:
|Vds , Vqs||: d-axis and q-axis stator voltages;|
|ids , iqs||: d-axis and q-axis stator currents;|
|ψdr ,ψqr||: d-axis and q-axis rotor flux linkages;|
|Rs , Rr||: stator and rotor resistances;|
|Ls , Lr||: stator and rotor inductances;|
|Msr||: mutual inductance|
|ωs ,ωr||: electrical stator and rotor angular speed|
|ωg||: slip speed ωg=(ωs-ωr)|
|Ωr||: mechanical rotor angular speed|
|Cr , Cem||: external load torque and motor torque|
|j , f||: inertia constant and motor damping ratio|
|p||: number of pole pairs|
|σ||: leakage coefficient, (|
The bloc diagram of the proposed indirect field-oriented induction motor drive is shown in fig.6. Speed information, obtained by encoder feedback, enables computation of the torque reference using a PI controller. The reference flux is set constant in nominal speed. For higher speeds, rotor flux must be weakened.
4. Simulation results
To demonstrate the performance of a new tree-phase inverter, we simulated three filtered inverter output voltages. Two frequency values are imposed on the inverter, starting with a frequency of 50 Hz, then at time t=0.06(s) the frequency is changed to 60Hz. The three filtered inverter output voltages are illustrated in (figure 7.a) which clearly shows that the three voltages are perfectly sinusoidal and follow the ideal values with a transient time of 0.012(s). The harmonic spectrum of the filtered output voltage is shown in (figure 7.b). In order to compare the performance of the new three-phase inverter with the conventional PWM inverter, the output voltage of the latter, where a modulation index of 35 was used, is shown in (figure 8.a) and its harmonic spectrum is presented in (figure 8.b).
We carried out two simulations of the field-oriented control for induction motor drives with speed regulation using the new structure of the three phase inverter in the first simulation (figure 9.a) and the conventional PWM (figure 9.b) inverter in the second simulation. The instruction speed is set to 100 (rad/sec) for both simulations. During the period between 1.3(s) and 2.3(s), a resistive torque equal to 10 (N.m) (i.e the nominal torque) is applied.
In order to illustrate the effectiveness of the proposed inverter, the torque response obtained by using the proposed and the conventional PWM inverters are shown in (figure 10.a) and (figure 10.b), respectively. The obtained results clearly show that the conventional PWM inverter generates more oscillations in the torque than the proposed structure (figure 11). Moreover, the switching frequency of the proposed inverter is dramatically reduced (see (figure 12.a)) when compared to its counterpart in the conventional PWM inverter (see figure 12.b). Therefore, the proposed inverter gives a better dynamic response than the conventional PWM inverter.
5. Experimental setup
The experimental setup was realized based on the DS1103 TMS320F240 dSPACE kit (dSPACE, 2006a, 2006b, 2006c, 2006d, 2006e). Figure 13 gives the global scheme of the experimental setup. This kit allows real time implementation of inverter and induction motor IM speed drive, it includes several functions such as Analog/Digital converters and digital signal filtering. In order to run the application the control algorithm must be written in C language. Then, we use the RTW and RTI packages to compile and load the algorithm on processor. To visualize and adjust the control parameters in real time we use the software control-desk which allows conducting the process by the computer.
The novel single phase inverter structure for pre-calculated switching is based on the use of IGBT (1000V/25A) with 10 kHz as switching frequency. The switching angles are predetermined off-line using Genetic Algorithms and stored in the card memory in order to speed up the programme running. The non-filtered inverter output voltages are first designed in Simulink/Matlab, then, the Real-Time Workshop is used to automatically generate optimized C code for real time application. Afterward, the interface between Simulink/Matlab and the Digital Signal Processor (DSP) (DS1103 of dSPACE) allows the control algorithm to be run on the hardware.
The master bit I/O is used to generate the required 2 gate signals, and a several Analog-to-Digital converters (ADCs) are used for the sensed line-currents, capacitors voltage, and output voltage. An optical interface board is also designed in order to isolate the entire DSP master bit I/O and ADCs. The block diagram of the experimental plant is given in figure 14
6. Experimental evaluation
Figure 15 shows the experimental filtered inverter output voltage (VC1-VC2) for frequency value equal 50 Hz. The filtered inverter output voltage is perfectly sinusoidal. The experimental result in Figure 16 shows the torque response obtained by using the proposed PCPWM inverter: during the period between 0.65 sec and 1.95 sec, a load torque equal to 13 (N.m) is applied. The torque ripple of the induction motor is dramatically reduced.
Moreover, as shown in Figure 17, the experimental switching frequency of the proposed PCPWM inverter is very less compared to the conventional PWM inverter one Figure 18. As a consequence, the proposed inverter provides higher dynamic response than the conventional PWM inverter in vector controlled induction motor applications.
A three-phase inverter model was developed by combining pre-calculated switching angles and a passive filter to eliminate inverter output harmonics. The inverter model needs a nonlinear system of equations for the switching angles computation. The proposed inverter model succeeds to substantially reduce the harmonics while using polarized capacitors. The reduced number of switching angles provides more reliability and increases system components life time. Moreover, the proposed inverter design and control simplicity could be used as a cost effective solution to harmonics reduction problem. The torque ripple of the induction motor is dramatically reduced by the PCPWM inverter. The global scheme of the experimental setup has been implemented. The obtained experimental results exhibit good matching with the theoretical values. It is shown that the proposed PCPWM has better tracking performance as compared with the conventional PWM.
The squirrel cage induction motor data are:
C1, C2, C3
Rated line voltage
Rated load torque
No. of pole pair
Stator leakage inductance
Rotor leakae inductance
Moment of inertia
Internal inductor resistance