International Journal of Electrical and Computer Engineering (IJECE) Vol. No. October 2015, pp. ISSN: 2088-8708 Optimal Design of Switched Reluctance Motor Using PSO Based FEM-EMC Modeling Mouellef Sihem. Bentounsi Amar. Benalla Hocine LGEC. Laboratory of Electrotechnics Dept. Mentouri University of Constantine 1. Algeria Article Info ABSTRACT Article history: This paper aims to optimize the design of a prototype of a 6/4 Switched Reluctance Motor (SRM) using the Particle Swarm Optimization (PSO) The geometrical parameters to optimize are the widths of the stator and rotor teeth due to their significant effects on the prototype design and the performances in terms of increased average torque and reduced torque ripple. The studied 3kW SRM is modeled using a numerical-analytical approach based on a coupled Finite Element Method with Equivalent Magnetic Circuit (FEM-EMC). The simulations are performed under MATLAB environment with user-friendly software. The optimal results found are discussed, compared against those obtained by the Genetic Algorithms (GA) and showed a significant improvement in average torque. Received Mar 25, 2015 Revised Jun 3, 2015 Accepted Jun 25, 2015 Keyword: Equivalent Magnetic Circuit Finite Element Method Genetic Algorithms Particle Swarm Optimization Switched Reluctance Motor Copyright A 2015 Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Mouellef Sihem. LGEC. Laboratory of Electrotechnics Dept. Mentouri University of Constantine 1, 47 Rue E-A-K. Khroub. Constantine. Algeria. Email: mouellef_sihem@yahoo. INTRODUCTION To meet challenging requirements, new design and more efficient structures of electrical machines are investigated by manufacturers and researchers. In this context. Permanent Magnet Synchronous Motor (PMSM). Brushless dc motor (BLDC). Linear Switched Reluctance Motor (LSRM) and rotary Switched Reluctance Motors (SRM. have been explored in the literature as they are an attractive alternative to induction and synchronous machines . Due to their robustness, reliability, high performance and reduced cost, the SRM found numerous applications at high speed drive or low speed generator . lectrical vehicles, air-conditioners, extractors, centrifuges, flywheel energy storage, shipbuilding, aeronautics, and gearless wind generato. Metaheuristic methods are general optimizing algorithms applicable to a wide variety of problems. They appeared in the 1980s, with a common ambition: to solve efficiently the difficult optimization problems, for which there is no known most effective classical method . , . New techniques inspired by artificial intelligence have emerged and developed to offer as potential alternative techniques to improve the quality of the solution, namely Genetic Algorithms (GA). Particle Swarm Optimization (PSO), and so on. The PSO is a still relatively unknown and relatively young technique in the field of design . , . It is analogous to GA in the sense that the system is initialized with a random population of solutions. is compared to all its neighbors by maintaining each time the best result . Unlike the GA and other metaheuristic algorithms. PSO has the flexibility to control the balance between global and local exploration of the search space . , . The PSO has achieved rapid development following advantages . simple concept, easy implementation, robustness and computational efficiency. In . , the torque production is improved using PSO algorithm to optimize the stator and rotor angles of a 8/6 SRM. In . , the PSO is applied to the rotor pole arc of a 4/2 SRM to minimize the torque ripple. Journal homepage: http://iaesjournal. com/online/index. php/IJECE A ISSN: 2088-8708 One aspect of the contribution of this work lies with the application of the Particle Swarm Optimization method (PSO) for optimizing the average torque of a 6/4 SRM through various geometric The other aspect is the comparative study of the performance of PSO and GA algorithms applied to the machine GA . The results show that the PSO-based approach gives the best performance in terms of solution quality, accuracy and convergence time. The main contribution of this work is related to the numerical-analytical approach used to model the studied SRM using a user-friendly program carried out under MATLAB. The paper is organized as follows: Section 2 describes the FEM-EMC approach modeling of the studied SRM. In Section 3 the PSO algorithm is presented with the formulation of the problem. The results obtained are discussed in Section 4. The paper concludes in Section 5. FEM-EMC MODELING OF THE SRM Modeling of electric machines can be classified into three categories: analytical models, finite element analysis (FEA) and equivalent magnetic circuits (EMC), which can be considered as a semianalytical method . Modeling using EMC has been chosen for further investigation because it seemed a good technique with great speed and acceptable accuracy. The model produced will be used later in an optimization process that aims to find the best system parameters. The machine topology studied is a double saliency three-phase 6/4 SRM with Ns=6 stator teeth and Nr=4 rotor teeth as represented "Figure 1". Its operating principle, similar to the stepper motor, has long been known: by exciting successively the three stator phases, the rotor teeth are positioned to maximize the inductance of the power phase, under the rule of 'maximum flux' . ligned positio. by turning off the power, the motor will continue its movement until it reaches a position corresponding to the minimum value of inductance or flux . naligned positio. On the linked flux (A)-current . characteristics, the area between the previous two extreme positions represents the electrical energy converted into mechanical energy per cycle. W=Wa-Wu, as shown in "Figure 2". As described in . , to determine analytically the relations flux-At from only seven characteristics equal-flux lines traced by the finite element method (FEM) and corresponding to seven magnetic equivalent circuits (EMC), we implemented a program in MATLAB package software for the iterative calculation of the saturated aligned and unaligned inductances, respectively La and Lu, and the corresponding energies. Wa and Wu, from which one can deduce the average torque as depicted in "Figure 3": T av A qN r AW a A W u A 2A W a A E A 1 A A 2 A . A A n E * A i E Wu A AuI p Ai A . The integration step Ai is the ratio of the peak value of current Ip on the number of intervals n. Figure 1. Cross-section of the studied 6/4 SRM IJECE Vol. No. October 2015 : 887 Ae 895 IJECE ISSN: 2088-8708 Aligned Position Magnetic Flux . Converted Energy W = (Wa-W. Unaligned Position 1000 2000 3000 4000 5000 6000 7000 8000 9000 Stator Excitation (A. Figure 2. Extremes magnetic characteristics flux vs. excitation mmf Figure 3. Flowchart of main program of simulation under MATLAB Optimal Design of Switched Reluctance Motor Using PSO Based FEM-EMC Modeling (Mouellef Sihe. A ISSN: 2088-8708 PARTICLE SWARM OPTIMIZATION ALGORITHM Basic Concepts of PSO The Particle Swarm Optimization method (PSO) is a relatively recent heuristic proposed by Eberhart and Kennedy for the first time in the early 90s . , based on a stochastic population candidates solutions to develop an optimal solution to the problem presented. This method is particularly suitable for non-linear systems. it does not require the calculation of the first and second derivative, unlike the gradient type methods. Its basic idea is inspired from the actions of animal groups . in their search for the best subsistence areas. Thus, each individual in the population has the memory of its previous experience and the information provided by the group on the most promising regions. This contribution to the overall experience, in addition to personal experience is one of the features of PSO which ensure it success in global A swarm of particles, which are potential solutions to the optimization problem, "flies" the search space, in the search for the global optimum. The movement of a particle is influenced by three components . A component of inertia: the particle tends to follow its current travel direction. A cognitive component: the particle tends to rely on its own experience, and thus to move towards the best site in which it has already passed. A social component: the particle tends to rely on the experience of its congeners, and thus to move towards the best sites already reached collectively by the swarm. Basic Principle of PSO In a search space of dimension D , the algorithm starts with a random initialization of the particle A Particle i of the swarm is modeled by its position vector xi A A xi1 , xi 2 ,. , xiD A and the velocity A vector vi A Avi1 , vi 2 ,. , viD A . The quality of its position is determined by the value of the objective function at that point. This remembers the best position in which it has already passed, which is noted A P best i A A pbest i1 , pbest i 2 ,. , pbest iD A . The best position achieved by its neighboring particles is A noted Gbest A A gbest1 , gbest 2 ,. , gbest D A . Indeed, at iteration t 1, the velocity vector and the position vector are calculated from the equation . and equation . , respectively. Au Ay Au vit,Aj1 A wvit, j A c1r1t, j pbestit, j A xit, j A c2 r2t, j gbesttj A xit, j xit,Aj1 A xit, j A vit,Aj1 , j Ea A1,2,. where : v it,Aj1 , v it, j : Are the speed of the particle to t and t 1 iterations. Pbest : Is the best position of the particle. Gbest : Is the best position of a neighbor at iteration t. x it, j : Is the position of the particle at iteration t. w : Is generally called a constant factor of inertia, it keeps a balance between exploration and c1 and c2 : are two constants called acceleration coefficients, they keep the balance between individual and social behavior of the particle when they are equal . r1 and r2 : are two randomly generated numbers with a range of . , for each iteration and for each dimension j. Problem Formulation The objective function f x, u used to formulate the SRM problem represents a maximizing average torque. In the case of an optimization problem where the objective is to be maximized, the function is considered with the opposite sign A f x, u . the equality constraints expressed by the function g x, u are represented by the equations of the maximum and minimum inductances La and Lu as well as Wa and Wu energies of the two extreme positions of which the average torque and inequality constraints will be deduced which reflect the lower and upper dental openings which are given by equations limits . IJECE Vol. No. October 2015 : 887 Ae 895 IJECE ISSN: 2088-8708 A sm A EE 2A qN EE A 30 A C A s C 45 A A EE A N EE A rm A A sm A 30 A C A r C 60 A A AA r A A sm A . AA s A A r A A a r A EE 2A N EE A 90 A . r E r E Therefore, the proposed solutions must take the constraints of construction into account. These constraints are taken into account by penalizing proportionally the objective function for constraint In the context of taking into account the constraints, it is to degrade the performance of infeasible individuals in function of their proximity to the feasible area of the search space. For each element of the search space, its proximity to the feasible region can be measured through the level of violation of each Using this measure of infeasibility of the individual x from each constraint, the penalty function in the general form can be introduced: Min Fobj A A f A x A A F penalty xi1 C xi C xiu , i A 1,. , n . SIMULATION RESULTS AND DISCUSSION A comparative study with Genetic Algorithms (GA) has been made to verify the performance of the proposed algorithm. The PSO and GA parameters used for simulation are summarized in Table 1. For the implementation of PSO, several parameters must be specified, such as acceleration factors ( c1 and c2 ), the inertia factor ( w ), the size of the swarms and the stop criterion. The PSO algorithm has been applied to the objective function according to the flowchart in Figure 4. Figure 4. The flowchart of adaptive PSO for SRM Optimal Design of Switched Reluctance Motor Using PSO Based FEM-EMC Modeling (Mouellef Sihe. A ISSN: 2088-8708 Table 1. Simulation parameters Population size:20 Generations:100 PSO Number of Particles:20 Iterations:100 Crossover rate:0. c1 = c2 =2 w = -0. Mutation rate:0. To confirm the performance of this method, a comparison of its results with the results of genetic algorithms has been made. The comparison is shown in Table 2. The convergence characteristics of the two methods: PSO and GA for a variable reluctance motor are shown in Figures 5, 6, and 7. According to the results, it can be noticed that the PSO explores a solution superior to the genetic algorithm for the same number of population and generation. Table 2. Comparison results of PSO and GA for Fa= 1691At Prototype Stator pole arc . Rotor pole arc . average torque (N. Optimization PSO It can be seen from figure 5, firstly, that the PSO algorithm converges toward the global optimum from the thirty sixth iteration . , while the convergence of the GA algorithm is reached at iteration . with an optimal value lower compared to the PSO algorithm. This proves that the power of convergence to the global optimum in the PSO method exceeds that of the method of genetic algorithms (GA), this will have a direct impact on the time required for convergence of the two methods. Furthermore, the robustness of the PSO algorithm is more remarkable. The difference in average torque between the two optimization methods . 6178Nm compared with 15. is virtually insignificant or negligible . slight difference of about 48%). This will confirm our findings in terms of robustness of the PSO convergence. The results presented in Figures 6 and 7 show variations of AAs and Ar around their optimum values. The safety constraints are also checked for these two angles. These are qualified in their ranges. Fitness Value PSO Generations Figure 5. Objective function IJECE Vol. No. October 2015 : 887 Ae 895 IJECE ISSN: 2088-8708 PSO Pole arc angle of stator Generations Figure 6. Change of As around the optimal value PSO Pole arc angle of rotor Generations Figure 7. Change of Ar around the optimal value CONCLUSION The paper proposes a permeance network modeling and presents a direct coupling between the finite element method (FEM) and the equivalent magnetic circuit method (EMC) to model the switched reluctance The modeling tool presented is designed to be integrated into an optimization process that modifies the geometry of the engine. The optimization method chosen is the particle swarm optimization (PSO), the stochastic nature, metaheuristics, allows the application to difficult and non-linear problems. The principle of the method is explained as well as the different coefficients of the algorithm and the influence they have on the evolution of the algorithm. The optimization procedure of the design of SRM using PSO is presented with the aim of maximizing the average torque for an efficient solution all by acting on the tooth geometry which has a great influence on motor performance. The difference in the average torque estimated by the two algorithms. PSO and GA, is negligible suggesting the simultaneous convergence to the same quasi-optimal solution. From the simulation results, it can be found that the PSO can lead to optimal feasible solution, and that is the relative ease of implementation and ability to provide reasonably good solutions for many combinatorial problems. REFERENCES