# Using Genetic Algorithms for Optimal Electromagnetic Parameters of SPM Synchronous Motors

<figure><img src="image_1.png" /></figure>

Check for
updates

Trinh Cong Truong<fnref n="1" /> ID, Thanh Nguyen Vu<fnref n="1" />,<fnref n="2" />* ID, Hung Bui Duc<fnref n="1" />,<fnref n="2" /> ID, Vuong Dang Quoc<fnref n="1" />,<fnref n="2" /> ID

Corresponding Author Email: thanh.nguyenvu@hust.edu.vn

Copyright: ©2023 IIETA. This article is published by IIETA and is licensed under the CC BY 4.0 license
(http://creativecommons.org/licenses/by/4.0/).

https://doi.org/10.18280/jesa.560601

**Received:** 18 September 2023

**Revised:** 20 November 2023

**Accepted:** 4 December 2023

**Available online:** 28 December 2023

**Keywords:**

surface-mounted permanent magnet synchronous motor (SPMSM), cogging torque, torque ripple, genetic algorithm (GA), finite element method

## ABSTRACT

The permanent magnet synchronous motors (PMSMs) have been widely used in industrial applications due to the high efficiency, reliable performance and different shapes and sizes. Based on the arrangement of permanent magnets (PMs), the PMSM can be split into two primary types, i.e., surface-mounted permanent magnet (SPM) motors and interior permanent magnet (IPM). For the SPM motor, PMS are mounted on the rotor surface, while the IPM has the magnets embedded into the rotor. The use of PMs for the PMSMs has eliminated the necessity for excitation currents, thanks to the high flux density and significant coercive force. The resulting absence of excitation losses plays a key role in enhancing overall efficiency. This research, introduces a multi-objective optimal design strategy for a surface-mounted PMSM, with the primary goal of achieving maximum efficiency while minimizing material costs. The optimization is carried out through the application of a genetic algorithm. In addition, a finite element method is proposed to validate a comprehensive assessment and comparison of the variances between the initial design and optimal design. The proposed methods are applied to the practical problem of 5.5 kW SPMSM. The FEM and calculation results showed that the motor's efficiency increased 0.5% and material cost decreased 15.2$ after the optimization process, both are the expected results.

## 1. INTRODUCTION

As we have known the permanent magnet synchronous motor (PMSM) has gained widespread adoption in various industrial applications owing to the reliability, compact size, high efficiency, and its adaptability to a range of shapes and sizes [1]. Especially, at low speeds, the PMSMs find extensive applications in fields such as ship propulsion, lifting, mining, and oil field exploitation [2, 3]. The PMSMs can be categorized into two main types based on the arrangement of the permanent magnet (PM): interior permanent magnet (IPM) motors, where magnets are embedded into the rotor, and surface-mounted PM (SPM) motors, where the PMs are mounted on the rotor surface. These motors are excited by PMs without using the excitation currents, leveraging their high flux density and significant coercive force. This unique characteristic eliminates excitation losses, leading to a substantial increase in efficiency [4].

In order to obtain the high torque, the fractional-slot concentrated winding is a commonly utilized configuration. This winding arrangement, termed a fractional-slot winding, involves a non-integer number of slots per pole per phase, less than one. When each coil of the winding is wound around each stator tooth, it is referred to as a concentrated winding. Motors featuring this type of winding offer significant advantages in

low-speed direct drive transmission systems, including high efficiency, reduced end turn length, high slot fill factor, and low copper loss. However, the fractional-slot configurations come with drawbacks, such as a slightly lower winding factor and higher harmonic content in the magnetomotive force (MMF) distribution compared to integral slot machines [5, 6]. Consequently, the motor design must carefully consider the appropriate combination of the number of pole pairs and slots based on its intended applications.

The genetic algorithm (GA) is a search heuristic commonly utilized to produce efficient solutions for optimization and search problems. Drawing inspiration from natural evolution, the GA incorporates techniques such as inheritance, mutation, selection, and crossover to generate solutions for optimization problems. Its versatility renders it a robust tool for global optimization and the analysis of extensive datasets [7]. Many researchers have extensively applied GAs to investigate optimal designs for PMSMs [8-13]. However, the discrete variation of variables and their incomplete correlation with each other and other parameters can result in alterations to output power following the optimization process. The magnetic strength of Neodymium-iron-boron magnets (NdFeB) are stronger than Samarium Cobalt (SmCo) magnets and Ferrite magnets. Moreover, they are relatively inexpensive.

<footnote n="1">¹ School of Electrical and Electronic Engineering, Hanoi University of Science and Technology, Hanoi 100000, Vietnam</footnote>

<footnote n="2">² Laboratory of High Performance Electric Machines (HiPems), Hanoi University of Science and Technology, Hanoi 100000, Vietnam</footnote>
In this paper, a surface-mounted permanent magnet synchronous motor (SPMSM) using the magnet of NdFeB N35 is designed. Subsequently, the GA is applied to determine the optimal dimensions that achieve maximum efficiency while minimizing material costs (such as electrical steel, copper, and magnets) with power conservation.

## 2. ANALYTICAL DESIGN

The electromagnetic torque ($T_e$) is computed via the expression as

$$T_e = \frac{\pi}{4} D_{is}^2 L \cdot TRV \quad (1)$$

where, $D_{is}$ is the inner diameter of stator, $L$ is the length of machine, TRV is the torque density following to the volume (kNm/m³), its value is from 35 kNm/m³ to 85 kNm/m³ [14-16]. However, to determine the rotor volume ($V_o$), it relies on a coefficient concept that depends on the cooling systems, as presented in [17], that is,

$$V_o = \frac{D_{is}^2 L}{T_e} \quad (2)$$

It should be noted that the value of $V_o$ will be taken depending on the output power and the method of cooling.

The magnet thickness can be calculated as:

$$h_{mag} = \frac{\mu_r g B_{air-g} \pi}{B_r 4 \sin\left(\frac{\rho_{mag}}{2}\right) - B_{air-g} \pi} \quad (3)$$

where, $\rho_{mag}$ is the electrical angle, $B_r$ is the magnetic flux density in rotor, $\mu_r$ is the magnet permeability, $g$ is the length of air gap, and $B_{air-g}$ is the magnetic flux density in the air gap between stator and rotor. The value on $B_{air-g}$ can be defined as:

$$B_{air-g} = \frac{4}{\pi} \sin\left(\frac{\rho_{mag}}{2}\right) B_{mag} \quad (4)$$

where, $B_{mag}$ is the magnetic flux density of PM.

The iron loss ($p_{iron}$) of the PMSM due to the hysteresis eddy current losses and additional loss can be defined as [18]:

$$p_{iron} = k_h f B_{mag}^{\alpha} + k_e f^2 B_{mag}^2 + k_a f^{1.5} B_{mag}^{1.5} \quad (5)$$

where, factors $k_h$, $k_e$, $k_a$ and $\alpha$ are given in the study [19].

In this paper, the iron loss and copper loss accounted for the main proportion of the total loss of the motor. Thus, we assumed that the total loss is the sum of the iron loss and the copper loss, which is given as below:

$$P_{total loss} = P_{Fe} + P_{Cu} \quad (6)$$

where, $P_{Fe}$ is the loss appearing in the stator, rotor and teeth of the motor, and $P_{Cu}$ is the copper loss which can be defined as follows:

$$P_{Cu} = N_s I^2 R \quad (7)$$

where, $N_s$ is the number of slots, $I$ is the phase current and $R$ is the resistance of the winding in each tooth.

## 3. OPTIMIZATION PROCESS

The optimization process of the PMSM is performed via the GA diagram as shown in Figure 1. It involves a multi-objective problem with numerous variables and constraints. The careful selection of variables is paramount, as it exerts a direct influence on the outcome. In this investigation, a meticulous evaluation resulted in the incorporation of seven variables within the optimization procedure as given in Table 1, where variables $x_1, x_2, x_3, x_4, x_5, x_6, x_7$ are respectively the inner stator diameter, air gap length, height of the stator yoke, height of the rotor yoke, width of the tooth, current density and turn number per tooth.

<figure><img src="image_2.png" /><figcaption>Figure 1. GA for a multi-objective genetic algorithm</figcaption></figure>

<table><caption>Table 1. Upper and lower boundaries of optimal variables</caption><thead><tr><th>Variable</th><th>Lower Boundary</th><th>Upper Boundary</th></tr></thead><tbody><tr><td>x<sub>1</sub> (mm)</td><td>120</td><td>140</td></tr><tr><td>x<sub>2</sub> (mm)</td><td>0.7</td><td>1</td></tr><tr><td>x<sub>3</sub> (mm)</td><td>12</td><td>17</td></tr><tr><td>x<sub>4</sub> (mm)</td><td>12</td><td>17</td></tr><tr><td>x<sub>5</sub> (mm)</td><td>12</td><td>17</td></tr><tr><td>x<sub>6</sub> (A/mm<sup>2</sup>)</td><td>4.5</td><td>5.5</td></tr><tr><td>x<sub>7</sub> (turns)</td><td>55</td><td>65</td></tr></tbody></table>

The variables $x_1, x_3, x_4$ and $x_5$ are related to steel volume and iron loss, whereas $x_2, x_6$, and $x_7$ are associated to the magnet dimensions and copper loss. The precise formulation of proposed constraints is essential to guarantee that the optimization outcomes align with the operational requirements of the PMSM. This study centers on three key parameters: stator volume, rotor volume, and magnet volume. The stator volume is proportional to the power of the motor. Thus, this value was kept constant in order to ensure the power conservation of the motor. From the Eq. (1), it can be seen that the output torque of the motor is proportional to the rotor volume, so this value also remains to keep the output torque as a constant. The PM is the part which produced magnetic flux in the air gap. Therefore, only a minor change in the magnet volume will cause a change in the magnetic flux in the air gap, thereby changing the power and the torque of the motor. Thus, this value was also kept as a constant. Each volume of the three parts (i.e., stator, rotor, PM) accompanied by an equality constraint: $1T \le B_{xy} \le 1.5T$; $1.3T \le B_{ry} \le 1.6T$; $1.5T \le B_{tooth} \le$
2T; 0.45 ≤ $k_{fill}$ ≤ 0.55, where $B_{sy}$, $B_{sy}$ and $B_{tooth}$ are the magnetic flux densities in the stator, rotor and tooth, respectively. The $k_{fill}$ is the slot fill factor. It should be noted that the volumes of the stator ($V_{stator}$), rotor ($V_{rotor}$) and the magnet ($V_{magnet}$) are considered as a constant. In addition, the values of flux density of the stator yoke, rotor yoke, and tooth of the PMSM are the four inequality constraints already given in the study [20].

Additionally, it is essential to define the objective functions, with a particular focus on minimizing total loss and material cost. The material cost includes expenses related to electrical steel, copper, and magnets, outlined as follows:

$$C = c_{Fe}(M_{stator} + M_{rotor} + M_{teeth}) + c_{Cu}M_{Cu} + c_{mag}M_{mag} \quad (8)$$

where, $c_{Fe}$, $c_{Cu}$, $c_{mag}$ are the cost per kilogram of electrical steel, copper and PM, respectively. The terms $M_{stator}$, $M_{rotor}$, $M_{teeth}$, $M_{Cu}$ and the $M_{mag}$ are respectively the masses of stator, rotor, teeth, copper and PM. These values are defined as follows:

$$M_{stator} = Y_1 \frac{4Y_2x_1x_3 - 4x_3^2}{x_1^2} \quad (9)$$

$$M_{rotor} = Y_1 \frac{(x_1 - 2K_3x_2)^2 - (x_1 - 2K_4x_2 - 2x_4)^2}{x_1^2} \quad (10)$$

$$M_{teeth} = \frac{Y_1}{x_1^2} \left( \frac{\pi}{4} (Y_2x_1 - 2x_4)^2 - Q_s \left( \frac{Y_2-1}{2} x_1 - x_3 - Y_5 \right) \left( \frac{Y_2+1}{2} x_1 - \frac{Y_7}{Y_5} x_3 - x_5 + Y_7 \right) + Q_s Y_8 + Y_8 \frac{h_w}{2} \left( \frac{Y_7}{Y_5} x_1 - \frac{Y_9}{x_6} 2x_6 Y_5 + \frac{Y_{10}}{x_6} 2x_6 \frac{V_{rotor}}{x_1^2} + \frac{Y_9}{x_6} 2x_6 \frac{Y_8}{\sqrt{x_7}} \left( \frac{x_7}{\left( \frac{Y_2-1}{2}x_1-x_3-Y_5 \right) \sqrt{1+a_1^2} \sqrt{x_6}} \right) \right) \right) \quad (11)$$

$$M_{Cu} = \frac{Y_9}{x_6} 2x_6 Y_5 + \frac{Y_{10}}{x_6} 2x_6 \frac{V_{rotor}}{x_1^2} + \frac{Y_9}{x_6} 2x_6 \frac{Y_8}{\sqrt{x_7}} \left( \frac{x_7}{\left( \frac{Y_2-1}{2}x_1-x_3-Y_5 \right) \sqrt{1+a_1^2} \sqrt{x_6}} \right) \quad (12)$$

where, the factors $Y_1$, $Y_2$, $Y_3$, $Y_4$, $Y_5$, $Y_6$, $Y_7$, $Y_8$ and $Y_9$ are respectively given as:

$$Y_1 = \frac{\pi}{4} 10^{-9} \gamma_{Fe} V_{rotor}, Y_2 = \sqrt{\frac{V_{stator}}{V_{rotor}}}, Y_3 = \frac{\mu_r}{B_{rmag}} - 1 + 1,$$

$$Y_4 = \frac{1}{4}(Y_2^2 - 1) - \frac{K_2+1}{2} Y_5, Y_5 = h_w + h_{so}, Y_6 = \frac{\pi K_5}{Q_s},$$

$$Y_7 = h_{so} b_{so} + \frac{1}{2} h_w b_{so} + \frac{\pi h_w Y_5}{Q_s},$$

$$Y_8 = 1,095 \sqrt{\frac{4I}{\pi}}, Y_{10} = \frac{\gamma_{Cu} l}{10^{-6}}.$$

The other parameters $V_{stator}$, $V_{rotor}$, $B_{rmag}$, $B_{mag}$, $h_w$, $h_{so}$, $b_{so}$, $\gamma_{Fe}$ and $\gamma_{Cu}$ represent the volume of stator, rotor, the operating remanence of the magnets at working temperature, the flux density above the magnets, the wedge height, the tooth tip depth, the opening width of the semi-closed slot; the mass density of electrical steel and copper, respectively. In addition, the mass of the PM can be determined via the PM volume ($V_{mag}$).

$$M_m = 10\gamma_m V_m 10^{-9} \quad (13)$$

The cost is considered as a function of the variables, i.e.,

$$C = f_1(x) = f_1(x_1, x_2, x_3, x_4, x_5, x_6, x_7) \quad (14)$$

The losses in the stator, rotor and tooth are computed as:

$$P_{stator} = \frac{Y_{11}}{(x_2 x_3 x_1^{-2})^{1,6}} + \frac{Y_{12}}{(x_2 x_3 x_1^{-2})^2} \quad (15)$$

$$P_{rotor} = \frac{Y_{11}}{(x_2 x_4 x_1^{-2})^{1,6}} + \frac{Y_{12}}{(x_2 x_4 x_1^{-2})^2} \quad (16)$$

$$P_{tooth} = \frac{Y_{13}}{(x_2 x_5 x_1^{-2})^{1,6}} + \frac{Y_{14}}{(x_2 x_5 x_1^{-2})^2} \quad (17)$$

for

$$Y_{11} = k_h f \left( \frac{B_m V_m (Y_3-1)}{2 V_{rotor} k_j} 10^{-6} \right)^{1,6},$$

$$Y_{12} = k_e f^2 \left( \frac{B_m V_m (Y_3-1)}{2 V_{rotor} k_j} 10^{-6} \right)^2,$$

$$Y_{13} = k_h f \left( \frac{\pi B_m V_m (Y_3-1)}{Q_s V_{rotor} k_j} 10^{-6} \right)^{1,6},$$

$$Y_{14} = k_e f^2 \left( \frac{\pi B_m V_m (Y_3-1)}{Q_s V_{rotor} k_j} 10^{-6} \right)^2,$$

where, $f$ is the frequency stacking factor and. and $k_j$ is the iron lamination. The winding resistance in each tooth can be defined as:

$$R = 2Y_{15}x_6x_7 \left( x_5 + \frac{V_{rotor}}{x_1^2} \right) + Y_{15}x_6x_7 \frac{Y_9}{\sqrt{x_6}} \mathrm{ceil} \left( \frac{x_7}{\mathrm{floor} \left( \frac{\frac{Y_2-1}{2}x_1-x_3-k_5}{Y_9} + \frac{1+a_1^2}{\sqrt{x_6}} \right)} \right) \quad (18)$$

where, $Y_{15} = \frac{\rho_{Cu}}{110^{-3}}$ and $\rho_{Cu}$ is the electrical resistivity of copper. The minimum of total loss means the maximum efficiency. By combining Eqs. (5), (6), and (7) together, the total loss is considered as a function of variables, that is,

$$P_{\text{total loss}} = f_2(x) = f_2(x_1, x_2, x_3, x_4, x_5, x_6, x_7) \quad (19)$$

## 4. OPTIMAL RESULTS

The test problem is a practical SPMSM of 5.5kW as presented in Figure 2. The main parameters are given in Table 2.

The distribution of magnetic flux densities in the SPMSM before and after optimization is pointed out in Figure 3 and Figure 4. It can be seen that the value of B focuses on the tooth with a maximum value of 2.084 (T), and for 2.091 (T) after optimization. This implies a marginal rise in flux density on the tooth, accompanied by a reduction in total loss. Consequently, this alteration will not adversely impact the operational parameters of the proposed motor. In addition, the output power is also increase from 5530W upto 5548W. The distribution of the back electromotive force (EMF) before and after optimization is shown in Figure 5. It shows that before optimization, the RMS of back EMF is 311.3V, and 311.9V for after optimization, showing an increase of 0.6V. However, both of these values are still below the supply voltage of 353.6V. It should be noted that the efficiency of the machine
increases by 0.33%, rising from 93.97% before optimization to 94.3% after optimization.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** 2D Model of SPMSM of 5.5kW, 12slots, 2p=8 (top) and winding diagram (bottom)</figcaption></figure>

<table><caption>**Table 2.** Main parameters of SPMSM</caption><thead><tr><td>Parameter</td><td>Value</td><td>Unit</td></tr></thead><tbody><tr><td>Rated Power</td><td>5500</td><td>W</td></tr><tr><td>DC voltage (DC bus)</td><td>500</td><td>V</td></tr><tr><td>Number of slot</td><td>12</td><td>slot</td></tr><tr><td>Number of pole pairs</td><td>4</td><td>pole</td></tr><tr><td>Rate torque</td><td>70</td><td>N.m</td></tr><tr><td>Frequency</td><td>50</td><td>Hz</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Magnetic flux density distribution in the SPMSM before optimization.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Magnetic flux density distribution in the SPMSM after optimization</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Back EMF distribution for before and after optimizations</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Output torque before and after optimizations</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Pareto front results</figcaption></figure>

The waveform of the output torque before and after optimizations is presented in Figure 6. The torque at the motor shaft is also increased after optimization, going from 70.4 Nm to 70.65 Nm. Therefore, after the optimization process, both power, efficiency, and torque at the motor shaft are elevated. The torque ripple remains nearly unchanged, at 25.9% before optimization and 25.8% after optimization. However, whether before or after optimization, the torque ripple values are relatively high, indicating a need for adjustments to reduce the torque ripple.

Figure 7 presents the Pareto front results of the objective functions. It can be observed that the two objective functions, namely material cost and total loss, exhibit an inverse relationship: as material cost decreases, total loss increases, and vice versa. In this paper, the author selects the point of the
minimum value of the sum of these two optimized functions to carry out the process of calculating parameters for the motor after optimization. The optimal results on main parameters of the SPMSM is given in Table 3. It should be noted that the efficiency after optimization is smaller that that before optimization.

<table><caption>**Table 3.** Optimal results on main parameters of the SPMSM</caption><thead><tr><th>Parameters</th><th>Before Optimization</th><th>After Optimization</th></tr></thead><tbody><tr><td>Stator outer diameter (mm)</td><td>204.68</td><td>203.29</td></tr><tr><td>Stator inner diameter (mm)</td><td>136.58</td><td>135.64</td></tr><tr><td>Tooth width (mm)</td><td>14.43</td><td>14.3</td></tr><tr><td>Slot height (mm)</td><td>17.91</td><td>17.93</td></tr><tr><td>Magnet thickness (mm)</td><td>3</td><td>3</td></tr><tr><td>Magnet electrical angle</td><td>140</td><td>139</td></tr><tr><td>Active length of iron core (mm)</td><td>136.58</td><td>138.46</td></tr><tr><td>Number of turns</td><td>56</td><td>56</td></tr><tr><td>Conductor diameter (mm)</td><td>1.72</td><td>1.7</td></tr><tr><td>Slot fill factor</td><td>0.54</td><td>0.55</td></tr><tr><td>Power (W)</td><td>5530</td><td>5548</td></tr><tr><td>Phase back EMF (V)</td><td>179.7</td><td>180.1</td></tr><tr><td>Efficiency (%)</td><td>93.97</td><td>94.3</td></tr><tr><td>Torque ripple (%)</td><td>25.9</td><td>25.8</td></tr><tr><td>Power factor</td><td>0.9</td><td>0.9</td></tr></tbody></table>

In order to improve the efficiency without changing the output torque and power, a skewing technique for the PM is proposed in this part.

To apply the cross-magnetization technique, magnets need to be divided into smaller sections. In this study, each magnet is divided into five parts. Subsequently, it is necessary to choose the inclination angles for each part of the magnet. Moreover, these angles should be selected to be symmetrical across the 0-degree angle, as indicated in the table. In this research, these angles are randomly chosen to demonstrate the effectiveness of the cross-magnetization technique in smoothing the output torque waveform of the motor. However, these angles can also be optimized in other processes to achieve the best results for reducing torque ripple to the minimum value. The skewing PM with different angles is given in Table 4.

<table><caption>**Table 4.** Skewing PM with different angles</caption><thead><tr><th>Slice</th><th>Proportional Length</th><th>Angle</th></tr></thead><tbody><tr><td>1</td><td>1</td><td>-6</td></tr><tr><td>2</td><td>1</td><td>-3</td></tr><tr><td>3</td><td>1</td><td>0</td></tr><tr><td>4</td><td>1</td><td>3</td></tr><tr><td>5</td><td>1</td><td>6</td></tr></tbody></table>

The model of skewing PMs and cogging torques of skewing PM are presented in Figure 8 and Figure 9. Obtained results after using skewing PM technique are given in Table 5. We observe that torque ripple has been significantly reduced after implementing the cross-magnetization technique. This result can be explained as follows: It is known that when the magnet is divided into segments, the total cogging torque is equal to the cogging torque caused by each magnet as it passes through each tooth.

<figure><img src="image_9.png" /><figcaption>**Figure 8.** Model of skewing PMs</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 9.** Model of cogging torque using skewing PMs
Notes: $T_{total} = T_{cog\_slice1} + T_{cog\_slice2} + T_{cog\_slice3} + T_{cog\_slice4} + T_{cog\_slice5}$</figcaption></figure>

<table><caption>**Table 5.** Obtained results after using skewing PM technique</caption><thead><tr><th>Parameters</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>Output power</td><td>5313</td><td>W</td></tr><tr><td>Shaft torque</td><td>67.65</td><td>Nm</td></tr><tr><td>Efficiency</td><td>94.3</td><td>%</td></tr><tr><td>Torque ripple</td><td>6</td><td>%</td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>**Figure 10.** Waveform of cogging torques in each slice after using the skew PM</figcaption></figure>

The waveform of cogging torques in each slice after using the skew PM is pointed out in Figure 10. It can be seen that each cogging torque generated by a magnet in one cycle will have a different waveform depending on the magnet skew angle we choose. When these cogging torques are synthesized, they will cancel each other out, significantly reducing the total cogging torque, thereby smoothing the output torque curve of the motor. It can be observed that if we can choose the optimal magnet skew angle, the cogging torque can be nearly eliminated. However, after implementing magnet skewing, the
average torque is reduced [3]. Furthermore, as the magnet skew angle increases, the torque ripple is reduced, but the average torque also decreases [4].

To maintain constant power, or in other words, preserve the required output torque, two methods can be employed: increasing the thickness of the magnets or reducing the air gap. While increasing the magnet thickness is effective, it raises material costs. Therefore, the option of reducing the air gap is chosen and implemented. As we know, the magnetic flux generated by the magnets is constant. However, depending on the size of the air gap, the leakage flux and stray flux will vary. A larger air gap results in larger leakage and stray flux, leading to greater power loss. Conversely, a smaller air gap reduces these flux components, resulting in a smaller power loss. Additionally, reducing the air gap does not impact material costs, meaning it does not affect the optimization objective mentioned earlier. The value of the air gap cannot be too small, as it would pose challenges in manufacturing. After the calculation process, a new air gap value of 0.65 mm is applied for simulation. The simulation results are presented in Table 6.

<table><caption>**Table 6.** Simulated results after using skewing PM technique with air gap of 0.65mm</caption><thead><tr><td>Parameters</td><td>Value</td><td>Unit</td></tr></thead><tbody><tr><td>Power</td><td>5524</td><td>W</td></tr><tr><td>Shalf torque</td><td>70.3</td><td>Nm</td></tr><tr><td>Efficiency</td><td>94.4</td><td>%</td></tr><tr><td>Torque ripple</td><td>6.8</td><td>%</td></tr></tbody></table>

The output torque before and after airgap length adjustment is presented in Figure 11. So, after adjustment, the motor power has been increased to the required value, with efficiency rising by an additional 0.1%, from 94.3% to 94.4%. The output torque also reaches 70.3 N.m, slightly exceeding the required value of 70 N.m. The torque ripple, after reducing the air gap length, has increased slightly from 6% to 6.8%. These values are all within an acceptable range and do not significantly impact the machine operation.

The back EMF before and after airgap length adjustment is presented in Figure 12. It shows that before adjustment, the back EMF is 297.7 and it increases to 310.6V for after optimization. Both of these values are still lower than the supply voltage of 353.6V.

<figure><img src="image_12.png" /><figcaption>**Figure 11.** Output torque before and after airgap length adjustment</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 12.** Back EMF before and after airgap length adjustment</figcaption></figure>

Figure 13 shows the harmonic components of back EMF. We can see that after using the skew PM technique, the harmonic components of 5th, 7th and 11th order of the EMF have significantly reduced. Therefore, the waveform of the EMF after using the skew PM technique has become much more sinusoidal. The magnetic flux density distribution at the air gap for different positions with the skew PM is presented in Figure 14. The maximum value is obtained 1.2 (T). The fux linkage density at the air gap for no and full loads with the skew PM is shown in Figure 15. The waveform is quite sinusoidal and smooth for both cases.

<figure><img src="image_14.png" /><figcaption>**Figure 13.** Harmonic components of back EMF</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 14.** Magnetic flux density distribution at the air gap for different positions with the skew PM</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>Figure 15. Flux linkage density at the air gap for no and full loads with the skew PM</figcaption></figure>

## 5. CONCLUSIONS

This paper has successfully proposed the GA to optimize the SPMSM of 5.5kW. The primary goal of the optimization process was to minimize two critical factors: total loss and material cost. To assess the effectiveness of the chosen variables and constraints in the optimization method, the motor efficiency was evaluated using the FEM. The results indicate that the optimized design not only significantly reduces material costs compared to the original design while maintaining the desired efficiency level but also induces favorable changes in other parameters like cogging torque and torque ripple. In addition, the skew PM technique has been also developed to improve the efficiency, output torque, cogging torque and the back EMF.

In the future work, it could be explored the development and optimization of additional objective functions, including torque output, cogging torque, and torque ripple. Additionally, other optimization methods such as particle swarm optimization, cultural algorithm, bee algorithm, and more could be employed to address optimization challenges and achieve even more optimal outcomes.

## ACKNOWLEDGMENT

This research is funded by Hanoi University of Science and Technology under project number T2022-PC-008.

## REFERENCES

[1] Zhang, H., Dong, Z., Zhou, J. (2017). Optimization design and analysis of permanent magnet synchronous motor based on VC. In 2017 20th International Conference on Electrical Machines and Systems (ICEMS) Sydney, Australia, pp. 1-4.
https://doi.org/10.1109/ICEMS.2017.8055957

[2] Crider, J.M., Sudhoff, S.D. (2015). An inner rotor flux-modulated permanent magnet synchronous machine for low-speed high-torque applications. IEEE Transactions on Energy Conversion, 30(3): 1247-1254.
https://doi.org/10.1109/TEC.2015.2412547

[3] Feng, L., Yu, S., Zhang, F., Jin, S., Sun, Y. (2021). Study on performance of low-speed high-torque permanent magnet synchronous motor with dynamic eccentricity rotor. Energy Reports, 8: 1421-1428.
https://doi.org/10.1016/j.egyr.2022.03.018

[4] Qinghua, L. (2005). Analysis, design and control of permanent magnet synchronous motors for wide-speed operation. Doctoral Thesis, National University of Singapore, Singapore.

[5] Hwang, C.C., Lyu, L.Y., Liu, C.T., Li, P.L. (2008). Optimal design of an SPM motor using genetic algorithms and taguchi method. IEEE Transactions on Magnetics, 44(11): 4325-4328.
https://doi.org/10.1109/TMAG.2008.2001526

[6] Zhao, W., Shen, H., Chai, W., Wang X., Kwon, B.I. (2018). Optimal design and experimental test of a SPM motor with cost-effective magnet utilization to suppress torque pulsations. IEEE Transactions on Magnetics, 54(11): 1-5.
https://doi.org/10.1109/TMAG.2008.2001526

[7] Mendaci, S., Allag, H., Mekideche, M.R. (2015). Multi-objective optimal design of surface-mounted permanent magnet motor using NSGA-II. The Applied Computational Electromagnetics Society Journal, 30(5): 519-526.

[8] Braiwish, N. (2016). Design optimisation of brushless permanent magnet synchronous motor for electric vehicles. Doctoral Thesis, Wolfson Centre for Magnetics, School of Engineering, Cardiff University, UK.

[9] Gao, J., Dai, L., Zhang, W. (2018). Improved genetic optimization algorithm with subdomain model for multi-objective optimal design of SPMSM. CEs Transactions on Electrical Machines and Systems, 2(1): 160-165.
https://doi.org/10.23919/TEMS.2018.8326463

[10] Tun, E.K.A.N., Tin, T. (2014). Design of conventional permanent magnet synchronous motor used in electric vehicle. International Journal of Scientific Engineering and Technology Research, 3(16): 3289-3293.

[11] Gu, W., Zhu, X., Quan L., Du, Y. (2015). Design and optimization of permanent magnet brushless machines for electric vehicle applications. Energies, 8(12): 13996-14008.
https://doi.org/10.3390/en81212410

[12] Liu, H. (2015). Design of high-efficiency rare-earth permanent magnet synchronous motor and drive system. Doctoral Thesis, University of Central Florida, USA.

[13] Liu, G., Liu, M., Zhang, Y., Wang, H., Gerada, C. (2020). High-speed permanent magnet synchronous motor iron loss calculation method considering multi-physics factors. IEEE Transactions on Industrial Electronics, 67(7): 5360-5368.

[14] Sundaram, M., Anand, M., Chelladurai, J., Varunraj, P., Smith, S.J., Sharma, S., Assad, M.E.H., Alayi, R. (2022). Design and FEM analysis of high-torque power density permanent magnet synchronous motor (PMSM) for two-wheeler e-vehicle applications. International Transactions on Electrical Energy System.
https://doi.org/10.1155/2022/1217250

[15] Seo, U.J., Chun, Y.D., Choi, J.H., Chung, S.U., Han, P.W., Koo, D.H. (2013). General characteristic of fractional slot double layer concentrated winding synchronous machine. Journal of Electrical Engineering and Technology, 8(2): 282-287.

[16] Jussila, H., Salminen, P., Niemelä, M., Pyrhönen, J. (2007). Guidelines for designing concentrated winding fractional slot permanent magnet machines. In 2007 International Conference on Power Engineering, Energy
and Electrical Drives, pp. 191-194.
https://doi.org/10.1109/POWERENG.2007.4380186

[17] Mishra, S., Sahoo S., Das, M. (2017). Genetic algorithm: An efficient tool for global optimization. Advances in Computational Sciences and Technology, 10(8): 2201-221.

[18] Hong, G., Wei, T., Ding, X. (2018). Multi-objective optimal design of permanent magnet synchronous motor for high efficiency and high dynamic performance. IEEE

Access, 6: 23568-23581.
https://doi.org/10.1109/ACCESS.2018.2828802

[19] Pluta, W.A. (2011). Core loss models in electrical steel sheets with different orientation. Przegląd Elektrotechniczny, 87(9b): 37-4.

[20] Ponomarev, P. (2013). Tooth-coil permanent magnet synchronous machine design for special applications. Doctoral Thesis, Lappeenranta University of Technology, Finland.