Design Optimization of an Automotive Permanent-Magnet Synchronous Motor by Combining DOE and NMGWO

: This study proposes an optimization methodology for automotive permanent-magnet synchronous motors (PMSMs) to achieve maximum efﬁciency, maximum average torque, and minimum torque ripple. Many geometrical parameters can be used to deﬁne the PMSM of an automobile. To identify the most signiﬁcant parameters for optimization, the fractional factorial design of the design of experiment (DOE) was employed for screening, considering the interaction effects. The central composite design was used to construct the proxy model between the optimization target and optimization variable, and the effectiveness of the model was judged. Aiming at the multi-objective optimization problem of a motor, a new mechanism for grey wolf optimizer (NMGWO) algorithm combining an elite reverse learning strategy, a local search strategy, and a nonlinear control parameter strategy is innovatively proposed. This algorithm was applied to solve the multi-objective optimization model. The numerical calculation results show that this is an effective optimization design method that can improve the performance of automotive PMSMs. The effectiveness of the NMGWO algorithm on the optimization results of permanent-magnet synchronous motors is veriﬁed by the experimental results.


Introduction
A PMSM has the advantages of a simple structure, small size, low noise, and high efficiency, which have attracted widespread attention in the automotive field [1,2].For electric instrument vehicles, the complex road conditions in oil fields require sufficient output torque during operation, high reliability, strong environmental adaptability, and high operating efficiency.Thus, the driving motor for electric instrument vehicles must have high efficiency, large output torque, and small torque ripple.It is necessary to continuously strive to increase investment in material research, torque improvement, increased efficiency, and optimized control [3].Thus, the optimization of motor performance is very important.
With the rapid development of permanent-magnet motor technology, Gillon and others have applied the DOE to optimize a permanent-magnet motor to improve the motor's electromagnetic torque and no-load back EMF [4,5].Kurihara and other scholars used methods based on the time-stepped finite element method and the response surface method of the DOE in the single-phase, capacitor circuit, permanent-magnet synchronous motor to maximize efficiency [6].The permanent-magnet synchronous motor design process involves several parameters, each having a different degree of influence on the optimization objective.Furthermore, the DOE can quantify the degree of influence of each parameter on the objective function and screen out the significant factors [7][8][9].Response surface methodology is often used in the optimal design of electromagnetic devices to fit a more accurate model [10][11][12][13].This is followed by optimization of the resulting target response surface.The above studies have achieved certain results in optimizing permanentmagnet synchronous motors by using the ability of the DOE to solve complex parameter designs.However, the interaction effect between factors and the effect of mixing situations were not considered in the analysis process, especially in screening significant factors.
The definition of a permanent-magnet synchronous machine is a multi-parameter problem.Meanwhile, optimizing the machine's parameters is very complex and requires considerable computing and experimental time.To avoid this issue, a fractional factorial design was used to reduce the number of experiments, conduct significant variable screening, and determine which variables significantly impact the performance indicators [14,15].The design of experiment (DOE) was adopted.The DOE test design method is a scientific method to deal with multi-factor testing [16,17].It can obtain better process optimization parameters [18] through reduced test times and is widely used in process and structure optimization [19].
In 2014, Mirjalili et al. proposed the grey wolf optimization (GWO) algorithm by simulating the leadership and hunting mechanism of grey wolf populations in nature.The simulation results showed that the GWO algorithm is highly competitive in exploration, exploitation, local optimum avoidance, and convergence [20].However, the GWO algorithm also has the same defects common to many population intelligence optimization algorithms, such as a low convergence accuracy, slow late convergence, and easy premature convergence.The simplex method has been used [21] to improve the GWO algorithm, based on the reflection, expansion, outer contraction, and inner contraction operations of the simplex method, to improve the current population of poorer individuals to avoid falling into the local optimum.The literature [22] demonstrates an improved grey wolf position update formula by introducing the particle swarm algorithm and combining the position information of the individual optimal value and the population optimal value to avoid premature convergence.A previous study [23] introduced the Cauchy variation operator to improve the probability of premature convergence of the GWO algorithm and used it for the optimal allocation of independent microgrid capacity.Regarding the multi-objective optimization problem of electric motors, the traditional grey wolf optimization (GWO) algorithm has the common shortcomings of numerous swarm intelligence optimization algorithms, including low convergence accuracy, slow convergence speed in the later stage, and easy premature convergence [20,24].The article proposes a multi-objective optimization design method for the PMSM based on the improved grey wolf optimization algorithm.
The optimization objectives are to maximize motor efficiency, average torque, and minimum torque ripple [25,26], primarily utilizing elite reverse learning strategies to initialize the population and increase its diversity.The nonlinear control parameter strategy has been adopted to balance the local and global search ability and accelerate the rate of convergence speed of the algorithm; the local search strategy of the downhill simplex algorithm was introduced to avoid the premature convergence stagnation of the population [27].
In this study, an optimization procedure based on the DOE and NMGWO methods was employed to obtain an optimum design of the automotive PMSM to achieve maximum efficiency, maximum average torque, and minimum torque ripple simultaneously.Then, a fractional factorial design was used to screen activity and identify the main factors that significantly influence the responses (efficiency, average torque, and torque ripple).To ensure a good convergence for the performance of the automotive PMSM, the improved grey wolf algorithm was innovatively proposed.All simulation experiments were performed using FEM.Finally, the validity of the optimization solution was experimentally verified.

Structural Design and Design Variables of the PMSM
According to the technical requirements of the PMSM, combined with the design theory of a permanent-magnet motor and the empirical formula, the performance parameters of the motor were determined.The characteristics of various rotor structures were compared and analyzed, and the built-in "V+1" rotor structure was selected, as shown in Figure 1.According to the performance parameters of the motor, the stator and rotor structure, main size, air gap, slot size, permanent-magnet size, and other important structural parameters were preliminarily designed.Meanwhile, to reduce the cogging torque and torque ripple, the stator winding was selected as an integer slot double-layer distributed winding, and a 48-slot/8-pole machine was selected.

Structural Design and Design Variables of the PMSM
According to the technical requirements of the PMSM, combined with the design theory of a permanent-magnet motor and the empirical formula, the performance parameters of the motor were determined.The characteristics of various rotor structures were compared and analyzed, and the built-in "V+1" rotor structure was selected, as shown in Figure 1.According to the performance parameters of the motor, the stator and rotor structure, main size, air gap, slot size, permanent-magnet size, and other important structural parameters were preliminarily designed.Meanwhile, to reduce the cogging torque and torque ripple, the stator winding was selected as an integer slot double-layer distributed winding, and a 48-slot/8-pole machine was selected.The stator slot size will affect the magnetic flux leakage and torque characteristics of the motor, so the stator slot opening, stator slot shoulder height, and stator slot shoulder angle were selected as design variables.As the designed rotor structure adopts the "V + 1" structure, the distance between the two layers of permanent magnets will affect the magnetic field distribution.In addition, the outer bridge thickness of the V-shaped permanent magnet can change the distance between the permanent magnets and affect the no-load magnetic leakage coefficient of the motor, thus it was selected as a design variable.Additionally, given that the air gap length is an important parameter in the motor design, it was selected as a design variable.The length of the motor core greatly influences copper loss, iron loss, and stray loss and was also selected as a design variable.In summary, six parameter variables, including the stator slot opening, stator slot shoulder height, stator slot shoulder angle, V-shaped permanent-magnet outer bridge thickness, air gap length, and iron core length, were selected as design variables.
The initial value of the design variable is not the optimal solution of the optimization objective [28,29].It is, therefore, necessary to optimize the design variables for the optimization objective.Generally, the optimal value range of design variables can be determined according to the following constraints: The stator slot size will affect the magnetic flux leakage and torque characteristics of the motor, so the stator slot opening, stator slot shoulder height, and stator slot shoulder angle were selected as design variables.As the designed rotor structure adopts the "V + 1" structure, the distance between the two layers of permanent magnets will affect the magnetic field distribution.In addition, the outer bridge thickness of the V-shaped permanent magnet can change the distance between the permanent magnets and affect the no-load magnetic leakage coefficient of the motor, thus it was selected as a design variable.Additionally, given that the air gap length is an important parameter in the motor design, it was selected as a design variable.The length of the motor core greatly influences copper loss, iron loss, and stray loss and was also selected as a design variable.In summary, six parameter variables, including the stator slot opening, stator slot shoulder height, stator slot shoulder angle, V-shaped permanent-magnet outer bridge thickness, air gap length, and iron core length, were selected as design variables.
The initial value of the design variable is not the optimal solution of the optimization objective [28,29].It is, therefore, necessary to optimize the design variables for the optimization objective.Generally, the optimal value range of design variables can be determined according to the following constraints: where B t is the mean magnetic density of the stator teeth, B j is the magnetic density of the stator yoke, J is the stator winding current density, and S f is the slot filling rate.In this study, 1.9 T, 1.6 T, and 70% are empirical values.When B t < 1.9 T and B j < 1.6 T, the motor has a good performance and long operating life, and S f < 70% is reasonable for installing windings.
According to the principle of symmetry, a smaller value and a larger value than the initial value were selected; they were set to a low and high level, respectively.The value range of the design variable obtained is shown in Table 1.

Fractional Factor Design and Screening
If we have k variables, and each variable has n levels, the number of simulations or experiments for a full factorial design can be given by: The number of simulations or experiments increases exponentially when the number of parameters increases.To avoid this, the fractional factorial design was used to reduce the number of simulations and determine which factors significantly affect the studied responses.
Confusion is unavoidable when the fractional factorial design is used.Confusion represents a merger among the contributions of the main effects of the parameters, between the contributions of the main effects and the interaction effects, or among the contributions of the different order interaction effects.According to the influence results, the interaction effect between variables can be defined as a no interaction effect, positive interaction effect, or reverse interaction effect.According to the number of variables in the interactive item, it can be divided into second order, third order, fourth order, etc. [16].According to the actual application, most interaction effects are non-existent or negligible, so the specific application is mainly the second-order interaction effect.
Due to the presence of confounding, significant variables are difficult to determine and may be mistakenly excluded without accurate analysis.To avoid this situation, we attempted to mix the main or second-order interaction effect with the higher-order interaction effect when selecting the resolution of the fractional factorial design.Meanwhile, for convenience of analysis, the following assumptions were made [30]: (i) If a group of confounders can be ignored, all effects contained in that group can also be ignored.(ii) If a second-order interaction effect is determined to be a significant term, the two variables that make up the second-order effect can be identified as significant.
The available factor design and resolution are listed in Table 2. "All" indicates the complete factorial design and Roman numerals indicate the resolution of the corresponding fractional factorial design.
According to Table 2, if there are six design variables, fractional factorial design schemes with resolutions of III, IV, and VI are usable.Since the fractional factorial design with a resolution of III is a hybrid of the main and second-order effects, it is not possible to directly determine which variables are significant among those that make up the hybrid.The fractional factorial designs with resolutions of IV and VI can separately estimate the main effect and second-order interaction effect, and the corresponding design scheme for IV has fewer runs.Thus, the 2 IV 6-2 design solution was selected, and the detailed fractional factorial design experimental table is presented in Table 3. A, B, C, D, E, and F represent the six variables.The experimental results can be obtained through finite element simulation.
Compared to actual experiments, it not only achieves the goal of screening significant variables but also saves time and economic costs.For the 2 IV 6-2 design scheme, a resolution of IV indicates confusion between the main and third-order interaction effects and between one second-order interaction effect and another.As higher-order interaction effects above the third order are not considered, only the remaining second-order interaction effects are mixed, as shown in Table 4.In the case of fractional factor analysis, only the "representatives" of the confounding items are usually considered for inclusion in the model.Taking the confusion "AB + CE" as an example, only the term "AB" is chosen during the analysis when the actual result is the joint contribution of "AB + CE." Thus, if the result of the representative AB is remarkable, we must determine which term(s) truly work.If the result of the representative AB is not remarkable, we can determine that all the interaction effect terms comprising the confusion are non-significant.
The Pareto chart can help determine the most significant factors.The abscissa of the Pareto chart is the absolute value of t obtained by the t-test.The critical value of t is determined according to the significance level α s and degree of freedom of error terms.The term with an absolute value of t that exceeds the critical value of t will be selected as the significant factor.The critical value of t value is equal to the margin error, ME, and its value is calculated using the following formula: where t αs (n e ) is the quantile of (1 − α s /2) of the t distribution, which can be obtained from the quantile table of the t-distribution, and n e is an effector.PSE is the pseudo-standard error in the Lawrence method when the degree of freedom of the error term is zero.
To ensure that the results of the significance test are correct, the value of α s is first taken as 0.05.The Pareto effect plot of the efficiency obtained when α s = 0.05 is shown in Figure 2a.For torque ripple, A, D, and E are significant variables in the main effects, and the significant confounding term is AF.The analysis of AF is the same as before; therefore, A, D, E, and F are significant variables.The "DOE mean plots" are used to establish a mutual corroboration with the above analysis to fix significant factors, as shown in Figure 3.The DOE mean plot is ideally suited for determining significant factors.The plot presents mean values for the two or more levels of each plotted factor.The mean values for a single factor are connected by a straight line and can be expressed as: where x represents design variables, i-level represents each level, Rxi-level is the response value of factor x at the i-level, and K is the number of response values of factor x at the ilevel.
The DOE mean plot complements the traditional analysis of the variance of designed experiments.The plot reveals which factors are important and then ranks the list of the important factors.The DOE mean plot is qualitative and does not provide a definitive answer to the question, but it does help categorize factors as "clearly important," "clearly not important," and "borderline important." Figure 3 shows that five variables are identified as significant variables, namely stator slot opening Ws, stator slot shoulder height Ht, external bridge thickness Tb of a V-shaped permanent magnet, air gap length δ, and iron core length Lef.These are also significant variables for efficiency.Moreover, the significant variables for both average torque and torque ripple are stator slot opening Ws, external bridge thickness Tb of a V-shaped permanent magnet, air gap length δ, and iron core length Lef.After screening significant variables, the non-significant variables in the subsequent design process were taken as the initial value; that is, the stator slot shoulder angle At was taken as 20°.For efficiency, A, D, E, and F are significant variables in the main effects, while AC and AF are significant confounding terms.For the confounding term AC, also known as "AC + BE," due to the strong significance of E and the critical significance of B, it can be posited that the confounding term is mainly influenced by BE; therefore, B and E are significant variables.For the confounding term AF, also known as "AF + DE," it is easy to conclude that A, D, E, and F are significant variables.In summary, A, B, D, E, and F are significant variables.
When α s = 0.05, the Pareto effect diagram of the average torque is shown in Figure 2b.For the average torque, A, D, E, and F are significant variables in the main effects, with significant confounding terms including AD, AE, and AF.For the confounding term AD, also known as "AD + EF," A, D, E, and F can be significant variables.For the confounding term AE, also known as "AE + BC," since B and C are not significant variables, it can be considered that the confounding term is mainly influenced by A; therefore, A and E are significant variables.The analysis of AF is the same as before, therefore A, D, E, and F are significant variables.
When α s = 0.05, for torque ripple, there is no significant interaction effect, and only A is a significant item in the main effect; therefore, the only significant variable is A. To obtain sufficient significant variables and make the subsequent response surface fitting process more accurate, the values should be reasonably increased.When α s = 0.10, the Pareto effect diagram of torque ripple obtained is shown in Figure 2c.
For torque ripple, A, D, and E are significant variables in the main effects, and the significant confounding term is AF.The analysis of AF is the same as before; therefore, A, D, E, and F are significant variables.
The "DOE mean plots" are used to establish a mutual corroboration with the above analysis to fix significant factors, as shown in Figure 3.The DOE mean plot is ideally suited for determining significant factors.The plot presents mean values for the two or more levels of each plotted factor.The mean values for a single factor are connected by a straight line and can be expressed as: where x represents design variables, i-level represents each level, R xi-level is the response value of factor x at the i-level, and K is the number of response values of factor x at the i-level.Through fractional factorial design and significance analysis, the screening of significant variables was completed on the basis of considering the interaction effects and confounding effects between variables.This laid the foundation for adopting the improved grey wolf algorithm for multi-objective optimization in the next step.

Fitting Model
The response surface methodology (RSM) is a statistical tool used to build an empirical model by finding the relationship between the design variables and response through a statistical fitting method.The true functional relationship between responses and variables is complicated to determine.Thus, the relationship has to be approximated.The second-order model is used to predict a curvature response more accurately under fewer variables.A second-order response surface model that might describe this relationship is The DOE mean plot complements the traditional analysis of the variance of designed experiments.The plot reveals which factors are important and then ranks the list of the important factors.The DOE mean plot is qualitative and does not provide a definitive answer to the question, but it does help categorize factors as "clearly important," "clearly not important," and "borderline important." Figure 3 shows that five variables are identified as significant variables, namely stator slot opening Ws, stator slot shoulder height H t , external bridge thickness T b of a V-shaped permanent magnet, air gap length δ, and iron core length L ef .These are also significant vari-ables for efficiency.Moreover, the significant variables for both average torque and torque ripple are stator slot opening Ws, external bridge thickness T b of a V-shaped permanent magnet, air gap length δ, and iron core length L ef .After screening significant variables, the non-significant variables in the subsequent design process were taken as the initial value; that is, the stator slot shoulder angle A t was taken as 20 • .
Through fractional factorial design and significance analysis, the screening of significant variables was completed on the basis of considering the interaction effects and confounding effects between variables.This laid the foundation for adopting the improved grey wolf algorithm for multi-objective optimization in the next step.

Fitting Model
The response surface methodology (RSM) is a statistical tool used to build an empirical model by finding the relationship between the design variables and response through a statistical fitting method.The true functional relationship between responses and variables is complicated to determine.Thus, the relationship has to be approximated.The secondorder model is used to predict a curvature response more accurately under fewer variables.A second-order response surface model that might describe this relationship is given by: where β is the regression coefficient, k is the number of design variables, and ε is a statistical error.
In RSM, further analysis can be conveniently performed if the natural design variables are transformed into coded variables.They are usually defined as dimensionless quantities with a mean of zero and the same spread or standard deviation.The regression coefficients do not influence one another and can be compared directly after coding.The coding rule is written as: The fitted model should be examined to ensure it provides an adequate approximation for response.The determination coefficient used for the quality estimate of the fitted model can be expressed as: where SS Model is the sum of squares explained by the fitted model, SS Error is the sum of squares of error, and SS Total is the total sum of squares.
Numerous experimental designs can be used to create a response surface.These designs include the central composite design (CCD) and the Box-Behnken design.The CCD has been widely used to fit a second-order response surface.In this study, the CCD was chosen to estimate the curvature properties of the response surface and to complete the modeling.The CCD comprises three kinds of experimental points, each with a coded coordinate value: The value of m is related to the rotatability and the number of design variables.According to the value of m, the CCD can be classified as follows (Figure 4): (1) When m = 2 k/4 , the design is central composite circumscribed design (CCC), where k is the number of design variables; (2) when the value of m is < 1, the design is central composite inscribed design (CCI); and (3) when m is = 1, the design is central composite face-centered design (CCF).
The value of m is related to the rotatability and the number of design variables.According to the value of m, the CCD can be classified as follows (Figure 4 In this study, the CCF design was used.Thus, the star points are at the center of each face of the cube, such that m = ±1.This type requires only three levels (−1, 0, 1) of each factor.The most significant variables for efficiency and cogging torque differ.Thus, the response surface designs are created correspondingly.In this study, the CCF design was used.Thus, the star points are at the center of each face of the cube, such that m = ±1.This type requires only three levels (−1, 0, 1) of each factor.The most significant variables for efficiency and cogging torque differ.Thus, the response surface designs are created correspondingly.
MINITAB 20 software was used to analyze the response surface design and fitting model.For convenience, x 1 symbolizes Ws, x 2 symbolizes H t , x 3 symbolizes T b , x 4 symbolizes δ, and x 5 symbolizes L ef .The estimated regression coefficients were calculated, and the three fitted second-order polynomial functions for efficiency, average torque, and torque ripple are as follows: The determination coefficients of the three fitting models are listed in Table 5, which indicates that the second-order polynomials have good accuracy.We cannot identify the optimal solution directly from these models; thus, a specific optimization algorithm is needed to search the optimal solutions.

NMGWO Algorithm Process and Validity Verification
The improved grey wolf optimization algorithm in this study combines the elite reverse learning strategy, the local search strategy based on the downhill simplex method, and the nonlinear control parameter strategy.NMGWO has a high convergence accuracy, a fast convergence speed, and is not prone to premature convergence.The specific algorithm process of NMGWO is as follows (Figure 5): To verify the convergence accuracy and speed of the NMGWO algorithm, this study selected nine benchmark functions to evaluate the algorithm's performance via numerical testing, as listed in Table 6.The comparison algorithms include the particle swarm optimization algorithm (PSO), biogeography algorithm (BBO), genetic algorithm (GA), differential evolution (DE), and basic grey wolf optimization (GWO) algorithm.Moreover, the downhill simplex method was replaced with the simplex method to obtain another improved grey wolf optimization (NIGWO) algorithm.This was compared with the numerical test results of NMGWO.Each algorithm was independently run 30 times in Matlab 2021, and the optimization performance of NMGWO was tested by counting the maximum function value, minimum function value, average function value, and standard deviation of the seven algorithms.To verify the convergence accuracy and speed of the NMGWO algorithm, this study selected nine benchmark functions to evaluate the algorithm's performance via numerical testing, as listed in Table 6.The comparison algorithms include the particle swarm optimization algorithm (PSO), biogeography algorithm (BBO), genetic algorithm (GA), differential evolution (DE), and basic grey wolf optimization (GWO) algorithm.Moreover, the downhill simplex method was replaced with the simplex method to obtain another improved grey wolf optimization (NIGWO) algorithm.This was compared with the numerical test results of NMGWO.Each algorithm was independently run 30 times in Matlab 2021, and the optimization performance of NMGWO was tested by counting the maximum function value, minimum function value, average function value, and standard deviation of the seven algorithms.
The population size of all seven algorithms was 30, and the maximum number of iterations was set to 500.For NMGWO, the sensitive parameter α = 3, and the termination tolerance adopted the default value of ε = 1.0 × 10 −4 .For NIGWO, the reflection coefficient ρ = 1, expansion coefficient γ = 2, and internal (external) contraction coefficient β = 0.5 in the simplex method; for PSO, the maximum speed ν max = 6, the maximum and minimum values of inertia weight were ω max = 0.9 and ω min = 0.2, respectively, and the acceleration constant c 1 = c 2 = 2; for BBO, the proportion of fixed habitats was rate = 0.2, and the mutation probability was P ' m = 0.1; for GA, the crossover rate P c = 0.8 and the mutation rate P m = 0.1; for DE, the scaling factor F = 0.5 and the crossover factor C r = 0.3.
For the statistical result to be more intuitive and visual, box plots were used to display the numerical results of these seven algorithms, as shown in Figure 6.Compared with the other six algorithms, NMGWO can find the best objective values for different benchmark functions.The box plots of NMGWO are at the lowest position, which indicates that the average value obtained by NMGWO was the smallest.The box plots of NMGWO show a smaller box thickness and fewer outliers, indicating that the NMGWO algorithm has a strong searching ability and good stability.
When multiple objectives need to be achieved in practical engineering application scenarios, it is difficult to obtain the absolute optimal solution due to the internal conflicts among these objectives.In addition, the optimization of one objective is at the cost of the deterioration of other objectives.The trade-off is to coordinate or compromise between the objectives to find a feasible solution that makes each objective as optimal as possible.Typically, such problems can be described by the following mathematical model: where ] is the objective function vector and x = (x 1 , x 2 , • • • , x n ) are the decision variables.

Expression of Function Dimension
Search Interval [−100, 100] [−5.12, 5.12] The multi-objective optimization problem mentioned above does not have an absolute optimal solution, and the feasible solution obtained by trade-off is called a non-inferior or Pareto solution.Generally, the feasible solution obtained by a traditional optimization method each time is part of the Pareto solution set.More Pareto solutions can be obtained by using the population intelligent optimization algorithm.These solutions form an optimal solution set called the Pareto optimal solution set.In practical problem solving, the Pareto optimal solution set is usually selected according to the demands for the optimization objective.
For the multi-objective optimization problem min f (X), let X 1 , X 2 ∈ X f (the set of feasible solutions), where X 1 weakly dominates X 2 if f (X 1 ) ≤ f (X 2 ) and X 1 dominates X 2 if f (X 1 ) < f(X 2 ).Let X*∈ X f .If X∈ X f that dominates X* does not exist, then X* is a weakly Pareto optimal solution; if X* prevails over X for all ∀X∈ X f , then X* is a full Pareto optimal solution; if X ∈ X f that weakly prevails over X* does not exist, then X* is a strong Pareto optimal solution.All of the solutions that satisfy the Pareto optimal solution condition can constitute the Pareto optimal solution set, which can be generally expressed as [31]: When multiple objectives need to be achieved in practical engineering application scenarios, it is difficult to obtain the absolute optimal solution due to the internal conflicts among these objectives.In addition, the optimization of one objective is at the cost of the deterioration of other objectives.The trade-off is to coordinate or compromise between the objectives to find a feasible solution that makes each objective as optimal as possible.Typically, such problems can be described by the following mathematical model: where is the objective function vector and 1 2 ( , , , ) are the decision variables.The multi-objective optimization problem mentioned above does not have an absolute optimal solution, and the feasible solution obtained by trade-off is called a non-inferior or Pareto solution.Generally, the feasible solution obtained by a traditional optimization method each time is part of the Pareto solution set.More Pareto solutions can be obtained by using the population intelligent optimization algorithm.These solutions form an optimal solution set called the Pareto optimal solution set.In practical problem solving, For the multi-objective motor optimization problem in this study, the efficiency and average torque in the optimization target are maximized.Therefore, the objective function shall be transformed.Multiplying the response model of the two by "−1" can be equivalent to the minimization problem.Meanwhile, the design variable has been coded in Section 4; the value interval of each variable is [−1, 1], and it is a continuous variable.Therefore, the mathematical model of the multi-objective optimization problem is: Referring to the multi-objective gray wolf algorithm, introducing an archive unit component and a leader selection strategy extends the previously proposed NMGWO algorithm into a multi-objective version [32].
(1) Archiving unit components The key aspect of the archive unit component, which saves or retrieves the nondominated Pareto optimal solutions up to the current iteration loop, is the archive controller.Since the archive unit component presets the maximum number of archives, the archive controller controls the presence of feasible solutions when they expect to enter the archive or when the number of archives has been maximized.During the iteration process, the undominated Pareto solution obtained up to the current iteration (the new solution) is compared with the solution in the archive unit.There are generally four possible scenarios: (i) At least one solution in the archive cell that (weakly) prevails over the new solution.
In this case, the new solution cannot be entered into the archive cell.(ii) The new solution (weakly) prevails over one or more solutions in the archive cell.In this case, the new solution can enter the archive unit.(iii) If the new solution and the solutions in the archive cell do not dominate each other, then the new solution should be added to the archive cell.(iv) If the archive cell reaches the maximum number of archives, the grid mechanism should first be run to rearrange the partitioning of the target vector space; then, the most crowded part of the space, i.e., that containing the highest number of solutions, should be found and one of the solutions should be deleted.Finally, new solutions are inserted into the least crowded part to improve the diversity of the Pareto optimal front.
(2) Leader selection strategy In GWO, alpha, beta, and delta wolves can be easily obtained.However, in multiobjective search space summarization, comparisons to derive a leader for the gray wolf population based on the aforementioned concept of Pareto optimality are not readily available.The leader selection strategy tries to take the least crowded portion of the archive unit from which the alpha, beta, and delta wolves are provided.The selection is carried out by roulette, and the probability of each section being selected, obtained by the grid mechanism, is as follows: where c is constant and greater than 1, and N i represents the number of Pareto solutions obtained in part i.This equation shows that the less crowded each segmented part is, the higher the probability of providing a leader.It should be noted that since there are three leaders, there will be cases where leaders are provided from different parts.If there are more than three solutions in the least congested part, three will be randomly assigned as leaders; if there are less than three solutions in the least congested part, leaders will continue to be selected from the second least congested portion; if the number of solutions is insufficient to provide a satisfactory number of leaders, leaders will continue to be provided from the other uncongested parts.
Refer to the previously mentioned parameter settings where the initial population size was set to 150 to increase the diversity of the Pareto solution.The Pareto front obtained after 500 iterations of optimization searching is shown in Figure 7.The point that meets the requirements of high efficiency, high average torque, and smaller torque pulsation should be the origin of the coordinates in the graph.However, there is no optimal solution for the distribution of the origin; therefore, the region close to the origin meets the requirements, and the principle of selection can lead to a significant improvement in each optimization objective.After comparing the optimal solution for each Pareto in the region, the point indicated by the red arrow in the three indicators is improved.Therefore, we make a tradeoff, and the point indicated by the red arrow in the figure is chosen as a compromise solution; the coded and actual values of the design variables are shown in Table 7.After optimization, the stator slot opening Ws is 2.67 mm, the stator slot shoulder height H t is 1.164 mm, the external bridge thickness T b of the V-type permanent magnet is 1.72 mm, the air gap length δ is 0.79 mm, and the core length L ef is 135 mm.For ease of expression and analytical intuition, let x 1 denote the stator slot opening Ws, x 2 denote the stator slot shoulder height H t , x 3 denote the external bridge thickness T b of the V-type permanent magnet, x 4 denote the air gap length δ, and x 5 denote the core length L ef .
compromise solution; the coded and actual values of the design variables are shown in Table 7.After optimization, the stator slot opening Ws is 2.67 mm, the stator slot shoulder height Ht is 1.164 mm, the external bridge thickness Tb of the V-type permanent magnet is 1.72 mm, the air gap length δ is 0.79 mm, and the core length Lef is 135 mm.For ease of expression and analytical intuition, let x1 denote the stator slot opening Ws, x2 denote the stator slot shoulder height Ht, x3 denote the external bridge thickness Tb of the V-type permanent magnet, x4 denote the air gap length δ, and x5 denote the core length Lef.:    After obtaining the optimized values of the design variables, a finite element simulation of the permanent-magnet synchronous motor was carried out, and the optimized performance indexes are shown in Table 8.Compared with the pre-optimization period, the efficiency was improved by 0.31%, the average torque was improved by 7.73%, and the torque pulsation was reduced by 13.54%.The high initial efficiency of the motor results showed a less significant improvement after optimization, but the overall performance indexes were significantly improved.Figure 8 shows the comparison of motor output torque before and after optimization.After obtaining the optimized values of the design variables, a finite element simulation of the permanent-magnet synchronous motor was carried out, and the optimized performance indexes are shown in Table 8.Compared with the pre-optimization period, the efficiency was improved by 0.31%, the average torque was improved by 7.73%, and the torque pulsation was reduced by 13.54%.The high initial efficiency of the motor results showed a less significant improvement after optimization, but the overall performance indexes were significantly improved.Figure 8 shows the comparison of motor output torque before and after optimization.Figure 9 shows the distribution of magnetic force lines and flux density at the moment of t = 0 after optimization.It can be seen that the magnetic force lines were uniformly distributed.The maximum flux density at that moment was 2.2 T. Throughout the simulation, the average flux density of the stator teeth was 1.157 T, the average flux density of the stator tooth tips was 1.226 T, and the average flux density of the stator yoke was 1.012 T, which are all within the constraint conditions.
In the whole optimization process, the role of the agent model is to quickly predict the motor's performance, and the heuristic algorithm's role is to find the optimal solution

PMSM Prototype Manufacturing and Test Verification
To verify the operating performance of the optimized built-in permanent-magnet synchronous motor as well as the validity and accuracy of the simulation results, the optimized motor design scheme was used to carry out the trial production of the physical prototype and build the test bench.In the process of building the test bench, the installation position of the torque-speed sensor is shown in Figure 10.The performance parameters of the motor and controller are shown in Table 9 and Table 10, respectively.

Trial production prototype
Torque and speed sensor Dynamometer

Measuring instrument
Motor controller   In the whole optimization process, the role of the agent model is to quickly predict the motor's performance, and the heuristic algorithm's role is to find the optimal solution based on the prediction information of the agent model and the optimization objective.When considering the parameter changes during the motor operation, these parameters can be used as constraints or optimization variables in the heuristic algorithm.

PMSM Prototype Manufacturing and Test Verification
To verify the operating performance of the optimized built-in permanent-magnet synchronous motor as well as the validity and accuracy of the simulation results, the optimized motor design scheme was used to carry out the trial production of the physical prototype and build the test bench.In the process of building the test bench, the installation position of the torque-speed sensor is shown in Figure 10.The performance parameters of the motor and controller are shown in Tables 9 and 10, respectively.

PMSM Prototype Manufacturing and Test Verification
To verify the operating performance of the optimized built-in permanent-magnet synchronous motor as well as the validity and accuracy of the simulation results, the optimized motor design scheme was used to carry out the trial production of the physical prototype and build the test bench.In the process of building the test bench, the installation position of the torque-speed sensor is shown in Figure 10.The performance parameters of the motor and controller are shown in Table 9 and Table 10, respectively.

Trial production prototype
Torque and speed sensor Dynamometer

Measuring instrument
Motor controller  According to the optimization results, an automotive permanent-magnet synchronous motor prototype was manufactured and the test was carried out, as shown in Figure 11.According to the optimization results, an automotive permanent-magnet synchronous motor prototype was manufactured and the test was carried out, as shown in Figure 11.It can be seen from the figure that the trial production prototype was fixed on the test bench, and the torque-speed sensor was installed between the motor output shaft and the dynamometer.The dynamometer was used to load and control the motor speed for the bench test.
The data for the permanent-magnet synchronous motor were collected through the bench test, and the rated external characteristic curve and efficiency cloud diagram of the test prototype were obtained, as shown in Figure 12 and Figure 13, respectively.The motor efficiency MAP diagram mainly reflects the distribution of motor efficiency under different speeds and torques.The points with the same efficiency were connected into a loop and directly projected onto the plane to form a horizontal curve; the loops with different efficiencies did not coincide.From the efficiency MAP diagram, we can obtain the maximum efficiency value that the motor can achieve, as well as the motor speed and motor torque corresponding to this value.It can be seen from the figure that the trial production prototype was fixed on the test bench, and the torque-speed sensor was installed between the motor output shaft and the dynamometer.The dynamometer was used to load and control the motor speed for the bench test.
The data for the permanent-magnet synchronous motor were collected through the bench test, and the rated external characteristic curve and efficiency cloud diagram of the test prototype were obtained, as shown in Figures 12 and 13, respectively.The motor efficiency MAP diagram mainly reflects the distribution of motor efficiency under different speeds and torques.The points with the same efficiency were connected into a loop and directly projected onto the plane to form a horizontal curve; the loops with different efficiencies did not coincide.From the efficiency MAP diagram, we can obtain the maximum efficiency value that the motor can achieve, as well as the motor speed and motor torque corresponding to this value.
tor efficiency MAP diagram mainly reflects the distribution of motor efficiency under different speeds and torques.The points with the same efficiency were connected into a loop and directly projected onto the plane to form a horizontal curve; the loops with different efficiencies did not coincide.From the efficiency MAP diagram, we can obtain the maximum efficiency value that the motor can achieve, as well as the motor speed and motor torque corresponding to this value.As can be seen from Figure 12, due to the existence of experimental errors and the limitation of limited sampling and other factors, the rated external characteristic curve of the test prototype was not sufficiently smooth.However, the constant torque interval and constant power interval can be vaguely seen.The torque characteristic curve shows that the torque of the motor was approximately 110 N•m-120 N•m in the constant torque interval.The power characteristic curve shows that the motor's output power fluctuated up and down by 55 kW in the constant power interval.These results show that the rated torque of the tested motor was approximately 110 N•m-120 N•m, the rated speed was approximately 4500 rpm, and the rated power was approximately 55 kW, which met the determined design requirements.
Analyzing the measured efficiency MAP diagram of the prototype, it can be seen that in the electric state, when the rotational speed was 4500 rpm and the torque was 118 N•m, the motor's efficiency was within the 93-94% range.When the motor speed was > 1000 rpm, the motor efficiency in most areas reached> 90%.This is consistent with the conclusion of the motor efficiency simulation MAP diagram.Observing the simulated and measured efficiency MAP diagrams, it is apparent that the coherence of the efficiency dividing line and the location of the dividing line corresponding to different efficiencies show a high degree of consistency.Due to the richer amount of data collected in the actual test, the measured efficiency MAP diagram is more coherent in the consistency of the dividing line.Therefore, this diagram can more consistently reflect the changes in efficiency.
The measured output torque fluctuation of the prototype is shown in Figure 14, with As can be seen from Figure 12, due to the existence of experimental errors and the limitation of limited sampling and other factors, the rated external characteristic curve of the test prototype was not sufficiently smooth.However, the constant torque interval and constant power interval can be vaguely seen.The torque characteristic curve shows that the torque of the motor was approximately 110 N•m-120 N•m in the constant torque interval.The power characteristic curve shows that the motor's output power fluctuated up and down by 55 kW in the constant power interval.These results show that the rated torque of the tested motor was approximately 110 N•m-120 N•m, the rated speed was approximately 4500 rpm, and the rated power was approximately 55 kW, which met the determined design requirements.
Analyzing the measured efficiency MAP diagram of the prototype, it can be seen that in the electric state, when the rotational speed was 4500 rpm and the torque was 118 N•m, the motor's efficiency was within the 93-94% range.When the motor speed was > 1000 rpm, the motor efficiency in most areas reached> 90%.This is consistent with the conclusion of the motor efficiency simulation MAP diagram.Observing the simulated and measured efficiency MAP diagrams, it is apparent that the coherence of the efficiency dividing line and the location of the dividing line corresponding to different efficiencies show a high degree of consistency.Due to the richer amount of data collected in the actual test, the measured efficiency MAP diagram is more coherent in the consistency of the dividing line.Therefore, this diagram can more consistently reflect the changes in efficiency.
The measured output torque fluctuation of the prototype is shown in Figure 14, with an average torque of 117.852N•m and a torque pulsation of 12.702%.This is within the error tolerance and is relatively consistent with the simulation results after optimization.

Conclusions
In this study, the maximization of motor efficiency, maximization of average torque, and minimization of torque pulsation were selected as the optimization objectives.Optimization variables were selected, and significance screening was completed.The NMGWO algorithm, which combines the gray wolf algorithm, elite backward learning strategy, downhill simplex method, and nonlinear control parameter strategy, was proposed to improve the accuracy and convergence speed of the traditional gray wolf algorithm.The standard test function was used to test the improvement effect and showed that this method effectively improves the performance of a permanent-magnet synchronous motor.Finally, the prototype test was carried out according to the optimized design scheme, and the measured efficiency MAP cloud chart, rated external characteristic curve, and output torque fluctuation chart of the motor were obtained through the test.The simulation results were compared with the test results to verify the validity of the NMGWO algorithm for the optimization of the permanent-magnet synchronous machine.

Figure 1 .
Figure 1.Initial configuration of the analysis model.

Figure 1 .
Figure 1.Initial configuration of the analysis model.

Figure 2 .
Figure 2. Standardized effect diagram with response as efficiency, average torque, and torque ripple.(a) Pareto chart for efficiency; (b) Pareto chart for average torque; (c) Pareto chart for torque ripple.

Figure 2 .
Figure 2. Standardized effect diagram with response as efficiency, average torque, and torque ripple.(a) Pareto chart for efficiency; (b) Pareto chart for average torque; (c) Pareto chart for torque ripple.

Figure 3 .
Figure 3. DOE mean plots.(a) Mean plot of efficiency; (b) Mean plot of average torque; (c) Mean plot of torque ripple.

Figure 3 .
Figure 3. DOE mean plots.(a) Mean plot of efficiency; (b) Mean plot of average torque; (c) Mean plot of torque ripple.

( 1 )
Cube point (or corner point): The coordinate values of each cube point are 1 or −1.(2) Center point: The three-dimensional coordinate values of the center point are zero.(3) Star point: One of the three-dimensional coordinate values of the star point is m or −m, whereas the others are zero.
MINITAB 20 software was used to analyze the response surface design and fitting model.For convenience, x1 symbolizes Ws, x2 symbolizes Ht, x3 symbolizes Tb, x4 symbolizes δ, and x5 symbolizes Lef.The estimated regression coefficients were calculated, and

Figure 8 .
Figure 8.Comparison of output torque before and after optimization.

Figure 8 .
Figure 8.Comparison of output torque before and after optimization.

Figure 9
Figure9shows the distribution of magnetic force lines and flux density at the moment of t = 0 after optimization.It can be seen that the magnetic force lines were uniformly distributed.The maximum flux density at that moment was 2.2 T. Throughout the simulation, the average flux density of the stator teeth was 1.157 T, the average flux density of the stator

Figure 9 .
Figure 9. Distribution of magnetic field line and magnetic flux density at t = 0.

Figure 10 .
Figure 10.Installation position diagram of torque speed sensor.

Figure 9 .
Figure 9. Distribution of magnetic field line and magnetic flux density at t = 0.

Electronics 2023 ,
12, x FOR PEER REVIEW 17 of 22 Distribution of magnetic field lines Magnetic flux density distribution

Figure 9 .
Figure 9. Distribution of magnetic field line and magnetic flux density at t = 0.

Figure 12 .
Figure 12.Rated external characteristic curve of the test prototype.Figure 12.Rated external characteristic curve of the test prototype.

Figure 12 .
Figure 12.Rated external characteristic curve of the test prototype.Figure 12.Rated external characteristic curve of the test prototype.

22 Figure 13 .
Figure 13.Measured efficiency MAP diagram of the prototype.

Figure 13 .
Figure 13.Measured efficiency MAP diagram of the prototype.

Table 1 .
Design variable value range.

Table 2 .
Available factor design and resolution.

Table 3 .
Fractional factorial design experiment table and simulation results.

Table 4 .
Mixed second-order interaction effects.

Table 5 .
Determination index of fitting model.

Table 7 .
Design variable value corresponding to the compromise point

Table 7 .
Design variable value corresponding to the compromise point.

Table 8 .
Comparison of performance indexes before and after optimization

Table 8 .
Comparison of performance indexes before and after optimization.

Table 9 .
Parameters of permanent-magnet synchronous motor

Table 10 .
Parameters of controller

Table 9 .
Parameters of permanent-magnet synchronous motor.

Table 10 .
Parameters of controller.