Issue 
Sci. Tech. Energ. Transition
Volume 79, 2024
Power Components For Electric Vehicles



Article Number  13  
Number of page(s)  10  
DOI  https://doi.org/10.2516/stet/2024008  
Published online  08 March 2024 
Regular Article
Study on the impact of uncertain design parameters on the performances of a permanent magnetassisted synchronous reluctance motor
^{1}
IFP Energies nouvelles, Institut Carnot IFPEN Transports Energie, 1 et 4 avenue de BoisPréau, 92852 RueilMalmaison, France
^{2}
Paris Saclay University, CNRS, SATIE, GifsurYvette, France
^{3}
CY Cergy Paris University, CNRS, SATIE, Cergy, France
^{*} Corresponding author: adan.reyesreyes@ifpen.fr
Received:
15
September
2023
Accepted:
29
January
2024
In this paper, deterministic and robust design optimizations of a permanent magnetassisted synchronous reluctance motor were performed to study the impact of different uncertain input parameters on the design. These optimizations were carried out using a surrogate model based on 2D finite element simulations. Different robust optimizations considering geometric and magnetic uncertain parameters were compared to the deterministic optimization. It was noticed that both geometrical and magnetic properties tolerances greatly impact the machines’ mean torque and torque ripple, whereas the magnetic properties tolerances had a more significant impact on the mean torque. In such a case, robust optimization is essential to find optimal and robust electric motor designs.
Key words: Synchronous machines / Robust design optimization / Manufacturing uncertainties / Surrogate models / Finite element analysis
© The Author(s), published by EDP Sciences, 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
With the increasing concerns over climate change, many measures have been adopted to reduce greenhouse gas emissions. For transportation systems, in order to replace internal combustion engine vehicles, Electric and Hybrid Vehicles (EV, HEV) have been intensively developed. In these vehicles, the electrical machine is one of their main components.
Among the different types of electrical machines used in electric vehicles, Permanent Magnet Assisted Synchronous Reluctance Machines (PMaSynRelM) are one of the most used machines nowadays [1–3]. Unlike SurfaceMounted Permanent Magnet Synchronous Machines (SMPMSM), PMaSynRelM exploits two types of torque: hybrid torque generated using permanent magnets and, reluctance torque which makes a profit of the machine’s saliency. Since SMPMSMs only generate hybrid torque, they depend solely on permanent magnets to achieve the required performances. For all these machines, optimization is often used to find the best design respecting all the required specifications.
Many optimization methodologies applied to electrical machines can be found in the literature [4–7]. Most of these methodologies can be described as deterministic since they do not consider any uncertainties in the input parameters. However, in practice, there are many discrepancies between the nominal (obtained via data sheets) and real (measurements) values of these parameters. These differences can be caused by manufacturing and assembly tolerances in the prototype as well as by the lack of precision in the magnetic properties of used materials. These variabilities can lead to degraded performances compared to the nominal performances simulated in the design phase. To reduce such deviations, the parameter uncertainties should be considered in the optimization procedure.
In opposition to deterministic optimization, robust optimization considers two types of input parameters: certain parameters also known as controllable parameters, and uncertain parameters. Controllable parameters are the same ones used in deterministic optimization whereas uncertain parameters are specific to robust optimization techniques. This type of parameter can take varying values due to the associated uncertainties: it is then modeled by a random variable and an associated Probability Distribution Function (PDF). The presence of random input variables for the simulator leads to random output variables and then, random objective and constraint functions. Various formulations of the resulting optimization problem are proposed in the literature [8–10] based on the expectation, probability, or quantiles of these random variables.
ReliabilityBased Design Optimization (RBDO) is a method used to obtain optimal and safe designs in the sense that the outputs of certain functions are inside a security domain, described by constraints. A robust or reliable design has therefore a high probability of respecting these constraints. Examples of this approach can be found in [11].
Worstcase optimization considers the extreme values as objective functions and/or constraints i.e., the maximum or minimum value of the outputs caused by the uncertainty propagation [12].
There is another very common formulation that was also adopted in this work: Robust Design Optimization (RDO). In this methodology, the expectation (average) of the objective function is optimized. To limit extreme values, a second objective based on the variances of the objective function can be added [13].
The computation of robustness metrics such as expectations or variances requires a large sample of the uncertain input variables and thus a large number of simulations. To limit this high computational cost, especially with the use of finite element simulations, metamodeling techniques coupled with the design of experiments are used to replace the costly simulations with predictions using the resulting surrogate models [14].
We will present in this paper a study on the impact of uncertain design parameters on the performances of a PMaSynRelM motor. It is a 3phase 10pole 60slot PMaSynRelM with a Machaon rotor structure (Fig. 1). It has an outer stator diameter of 220 mm and an active length of 200 mm. Each pole has 3 flux barriers and 7 PMs. This machine topology will be used to compare the results of deterministic and robust optimizations in terms of robustness. To do so, techniques like Design Of Experiments (DOE), Finite Element Method (FEM) surrogate modeling, sensitivity analysis, quasiMonte Carlo methods, and optimization algorithms were used.
Fig. 1 Geometry of the PMaSynRelM motor. 
The remainder of this paper is as follows: We will first introduce the surrogate modelsbased methodology used in this work. Next, we will analyze the obtained results via deterministic and robust formulation optimizations. Finally, we will verify some solutions’ performances obtained by the metamodels through Finite Element simulations.
2 Optimization workflow
As mentioned before, we will perform in this paper deterministic and robust optimizations in order to study the sensitivity of a PMaSynRelM machine to design parameter uncertainties. Figure 2 shows the design parameters for the stator and the rotor considering one layer. Table 1 lists all the design parameters to be used in the optimizations as well as their lower and upper bounds. The manufacturing tolerance is also given.
Fig. 2 Design parameters for one layer (# is the number of the layer). 
Optimization variables.
Geometrical parameters are not the only design parameters that can be uncertain. The characteristics of the magnetic materials used to build a machine, like permanent magnets and electrical steel, can also present some uncertainties and their characteristics can deviate from those given in the datasheet.
The permanent magnet model used in our electromagnetic simulations is a linear model defined by two quantities: B_{r} which is the remanent induction and H_{cb} which is the coercitive magnetic field (Fig. 3). We will consider that the values of B_{r} and H_{cb} can vary between their nominal characteristics and, in the worst case, lower characteristics representing a 6.5% degradation. A coefficient β between 0% and 6.5% will be used as a random variable describing the level of degradation in the robust optimizations.
Fig. 3 Nominal and degraded PM characteristic. 
As for the B(H) characteristic, the degradation will be limited to its “knee” part. A random variable α between 0 (representing fully degraded characteristic) and 1 (nominal characteristic) will be used to define the degradation level of the B(H) curve. At maximum degradation (α = 0), the induction level is reduced by 25%. Figure 4 shows a comparison between the initial and degraded B(H) curves.
Fig. 4 Nominal and degraded B(H) curve. 
Figure 5 shows the workflow used to carry out the optimizations. At first, a DOE was built with the upper and lower bounds of the input parameters shown in Table 1 to fit a surrogate model for each of the considered objective functions: the mean torque and torque ripple. These two quantities are computed using the finite element method with the opensource software FEMM/XFEMM [15]. Different simulations at a fixed maximum current (600 A) and for a range of the control angle (between 0° and 90°) were launched. The optimal control angle (giving maximum torque) was found, and the mean torque and torque ripple were calculated. The built surrogate models were used to perform a global sensitivity analysis to detect the most impacting parameters on the objective functions. This analysis will allow us to limit the number of parameters considered as uncertain. At last, and after performing the metamodelbased deterministic and robust optimizations, FEM simulations were carried out to verify the results.
Fig. 5 Optimization workflow. 
2.1 Surrogate models
To reduce computation time, surrogate models have been built for each of the objective functions. To build such models, there are three steps to follow: build a DOE, train the metamodel, and check its predictivity with a test set. The chosen DOE is a maximin Latin Hypercube Sampling (LHS) which is a technique that covers well the search space while preserving good projection properties [16]. This DOE was built with 234 points using the bounds described in Table 1. As for the surrogate model, it is a universal Kriging with a linear trend function. For mean torque, a tensorized Matérn 5/2 covariance function has given the best predictivity. As for Torque ripple, a tensorized absolute value exponential kernel was used since the Torque ripple function is not that smooth. Kriging was chosen as a metamodel since it is very good at learning nonlinear objective functions and has demonstrated good performances in electrical machine optimization [17]. Finally, to evaluate the accuracy of the metamodel, a Normalized Root Mean Square Error (NRMSE) was computed on a validation test set:(1)zero indicates a good model fit. The obtained NRSME of the Kriging model for mean torque is 0.2% and for torque ripple is 8%. These results were obtained with a train and test sets composed of 175 and 59 samples, respectively. Torque ripple depends not only on mean torque but also on torque amplitude which makes this function more difficult to model than mean torque. We consider those metamodels sufficiently accurate for performing the sensitivity analysis and the optimization procedures.
2.2 Sensitivity analysis
For the robust optimizations, we have decided not to consider all the geometric parameters given before as uncertain parameters. Only the most impacting ones will be considered. To do so, a sensitivity analysis will be performed [18]. For this work, the Sobol Indices were chosen as they measure the global impact of the input variables on the output functions. Such indices represent the amount (or percentage) of the total output variance attributable to each subset of input variables. For instance, in the simplest case of one output variable Y and two inputs X_{1} and X_{2}, we find three factors causing the total variance of Y: due to the variation of X_{1} alone, due to the variation of X_{2} alone, and due to the variation of X_{1} and X_{2} simultaneously. Dividing these three quantities by the total variance of Y, we can obtain three percentages that can be directly considered as sensitivity measures. The largest is the value of these indices, the largest is their importance. The commonly used ones are the first order indices (S_{i}) representing the contribution to the variance of Y due to the variation uniquely of X_{ i }, and the total indices (S_{TOTAL,i }) representing the contribution to the variance of Y due to the variation of X_{ i } and all its interactions with the remaining input variables:(2) (3)where X_{ ~i } = X_{1},…,X_{ i−1}, X_{ i+1},…,X_{ Nx } and N_{ x } is the number of optimization parameters. We computed the indices using the Kriging surrogate models. The results of the sensitivity analysis applied to mean torque and torque ripple are presented in Figures 6 and 7. Only the most important geometrical inputs are displayed for better visibility. It was found that the stator slot width opening angle (Slot_angle) and the flux barrier opening angles for barriers 1 and 2 (Beta_L1_P1, Beta_L1_P2, Beta_L2_P1, and Beta_L2_P2) have the biggest impacts on the mean torque as well as on the torque ripple. We can also see in Figure 6 that the total order Sobol indices for torque ripple are higher than the firstorder ones. This means that torque ripple is more sensitive to the interaction between the variables (Beta_L1_P1 and Beta_L1_P2 for example) than to one variable alone. This is not the case for the mean torque, where the predominant indices are firstorder ones.
Fig. 6 Mean torque Sobol’ indices. 
Fig. 7 Torque ripple Sobol’ indices. 
Finally, the most impacting variables 5 parameters: Slot_angle, Beta_L1_P1, Beta_L1_P2, Beta_L2_P1, and Beta_L2_P2 will be considered as uncertain variables in the robust optimizations. The dispersion of these variables will be then integrated into the robust optimizations by considering a perturbation vector U.
2.3 Robustness metrics
To perform a robust optimization, one or several robustness measures can be used as objectives or constraints. In this section, we will recall some of them. Consider Figure 8. This figure shows a possible Probability Density Function (PDF) of the torque ripple considering perturbations on one (or some) input design parameters. The alphaquantile q_{ α } is the output value for which a given sample has a probability equal to α of being less than or equal to q_{ α }. For example, in Figure 8, the 10% quantile (q_{10%}) has a torque ripple value around 5.5%. This means that 10% of the considered samples have a torque ripple less than or equal to 5.5%. The expectation (E) is equal to the average value of the output variable and would coincide with the 50% quantile in a perfectly symmetric PDF. We can quantify how spread out the output values are by means of the Standard Deviation (STD), which measures the mean square dispersion around the expectation.
Fig. 8 Example of a torque ripple PDF. 
3 Deterministic and robust optimizations
We will present in this section the results of different optimizations: a deterministic and a robust optimization. For the deterministic optimization problem, we have:(4)where f_{1} is the opposite of the mean torque (in order to maximize it) and f_{2} the torque ripple; X is the controllable parameters space defined in Table 1. For the robust optimization problems, two formulations have been considered:

Expectations optimization:

Worstcase optimization:
u_{ uj } is the manufacturing tolerance of the parameter number j, and x_{ lj } and x_{ uj } are the lower and upper bounds for the optimization variable x_{ j }, respectively. The tolerance for each geometrical parameter is given in Table 1. Based on the sensitivity study in the previous section, only 5 geometrical parameters will be considered as uncertain. Their uncertainties were considered following uniform distributions, i.e., U_{ j } ~ Unif(−u_{ uj }, u_{ uj }). For parameters with no considered uncertainties, u_{ u } is simply equal to 0.
The goal of the first robust formulation is to optimize the mean torque and torque ripple’s expectations in a Pareto sense. In the worstcase formulation, the objective is to limit the worst possible value of the mean torque and torque ripple caused by uncertainties (equivalent to q_{100%} and q_{0%} respectively). By adopting such a formulation, the designer can be sure that all the manufactured machines will exceed the performances found on the Pareto front. On the other hand, a formulation using expectations does not guarantee that.
To solve these optimization problems, the genetic algorithm NSGA 2 was used [19]. This algorithm has shown good performances in other studies of electrical machine optimization as in [20]. We used a DOE maximin LHS to compute samples of x + U to calculate the objective functions’ expectations with a quasiMonte Carlo method. When it comes to the Worstcase formulation, we have two options: computing samples of x + U and taking the maximum value of these samples as an estimator of max_{uϵΩ}f(x + u) or obtaining the absolute maximum value with an optimization algorithm. In this work, the last solution was adopted using a Particle Swarm Optimization (PSO) [21] algorithm for the embedded monoobjective optimizations. By doing so, we will obtain more accurate worstcase estimates.
To study the impact of the two different types of uncertainties (related to geometrical parameters and material properties) on the performances of the PMaSynRelM machine, we will perform in the following sections two groups of optimizations: the first one will only consider geometrical parameters as uncertain parameters (U_{g}) while the second one will also consider the magnetic material characteristics as uncertain as described before (U_{g}, U_{m}).
3.1 Robust optimizations considering uncertainties on geometrical parameters
Figure 9 shows a comparison between the deterministic (blue) and the expectations optimization (red) Pareto fronts. The expected performances of the deterministic Pareto front have been reevaluated (pink): The design variables were perturbed by adding sampled values of the uncertain variables. These expected values represent the average mean torque and average torque ripple for each machine obtained by the deterministic Pareto optimization considering a posteriori uncertainty on the input parameters. As we can notice, a deterministic optimization does not guarantee a robust design: For the same average mean torque, the average torque ripple of a batch of machines issued from the deterministic optimization (pink) shows higher values than that of a batch of machines issued from the robust optimization (red).
Fig. 9 Pareto fronts obtained by deterministic (front in blue and expectations obtained by posterior perturbations in pink), and expectations optimization (red) where uncertainties come from geometrical perturbations U_{g}. Dark gray zones highlight points with similar mean torque expectation values. 
The Pareto front of the worstcase optimization (green) is presented in Figure 10. As in Figure 9, the deterministic Pareto front was also added (blue). The worstcase performances of the deterministic Pareto front have been evaluated (light green) thanks to a posteriori uncertainty on the input parameters and PSO maximization. Once again, these results show the importance of robust optimization in limiting the performance degradation that a sample of machines can have.
Fig. 10 Pareto fronts obtained by deterministic (front in blue and worstcases obtained by posterior perturbations in light green), and worstcases optimization (green) where uncertainties come from geometrical perturbations U_{g}. Dark gray zones highlight points with similar mean torque worstcase values. 
To go deeper into this analysis, we empirically compared the distribution of different designs. For this purpose, we show in Figures 11 and 12 boxplots of a subset of points selected from the deterministic and robust Pareto fronts shown in Figures 9 and 10, respectively. Each pair of boxplots in Figure 11 represents a comparison of the distribution of torque ripple values between a deterministic machine (pink in Fig. 9) and a robust machine (red in Fig. 9) falling in one of the zones (, , , , ). The “Exp” index means that the design is issued from a robust optimization based on the expectance formulation. The “Ug” is to say that the considered uncertainties are geometrical parameters. The same boxplots are represented in Figure 12. However, this time, the designs were selected from the deterministic Pareto front in Figure 10 (light green) and the robust one (dark green). All the boxplots in Figures 11 and 12 show the values of q_{25%} (q1), and q_{50%}, which is the median of the sample (q2) and q_{75%} (q3) as represented in Figure 11 (boxplot zone ).
Fig. 11 Boxplots showing comparisons of predicted torque ripple values between deterministic (blue) and robust (red) machines from Figure 9 with similar predicted mean torque expectation values. 
Fig. 12 Boxplots showing comparisons of predicted torque ripple values between deterministic (blue) and robust (green) machines from Figure 10 with similar predicted mean torque worstcase values. 
We can observe in Figure 11 that for each pair of machines, the robust one shows better overall performance than the deterministic one. For example, for the machines in zone , we can notice that the median of the deterministic machine is almost the same as the q3 value of the robust machine (4.5%). This means that there is a 50% chance for a deterministic design to have a torque ripple value higher than 4.5%. For the robust design, the probability is only 25%. We can also notice that the expectation and Standard Deviation (STD) of the torque ripple associated with the robust optimization solution (4.3%, and 0.35%, respectively) outperform the expectation and STD obtained by the posterior analysis of the deterministic solution (4.6%, 0.56%, respectively). These results stress the importance of robust optimizations when dealing with uncertainties.
Figure 12 also shows pairs of boxplots comparing predicted torque ripple values between deterministic (light green) and robust (dark green) machines from Figure 10 with the same minimum (worstcase) predicted mean torque value in each zone. For all the selected zones, the worstcase torque ripple value (maximum value for a batch) of a robust design is lower (or equal) than a deterministic one. Besides this, robust solutions also have much lower STD values than deterministic solutions, especially for zones with low worstcase mean torque. In zone for example, the robust machine (first green boxplot in Fig. 12) has an STD of 0.26% and a worstcase torque ripple of 6.7% compared to 0.54% and 7.2%, respectively for its deterministic counterpart (first blue boxplot in Fig. 12). Both optimizations lead practically to the same design for high values of worstcase mean torque like in zone . This can also be seen in Figure 11 with the Pareto fronts of the deterministic and robust optimizations getting very close with increasing torque.
The worstcase optimization has allowed to limit the performance degradation of the least performant machine in a sample, it has led to worse expectance values of torque ripple than those obtained by deterministic optimization. Instead of using the worstcase as an objective function, it could be used as a constraint in a constrained optimization problem while still using the expectations as an objective. Such formulation allows having good machine samples (expectationwise) while limiting the worst performances we can have. Given that the quantile and worstcase are related (as the worstcase is equal to the 100% quantile in a minimization scheme), and we have found that the worstcase value could be set as a constraint; a quantile formulation (like 95% quantile) could be used as a more permissive constraint. This is important because finding samples having a 100% chance of satisfying a constraint could be difficult, especially when the tail of the PDF of the quantity of interest is large.
3.2 Robust optimizations considering uncertainties on geometrical parameters and material properties
In order to evaluate the impact of magnetic material uncertainties on the performances of the PMaSynRelM motor, Figures 13 and 14 show the results of two robust optimizations performed using respectively the expectations and worstcase formulations. This time, the uncertain parameters are both the geometric parameters and the magnetic material properties. We added the same figures the Pareto fronts from Figures 9 and 10 for comparison.
Fig. 13 Pareto fronts obtained by expectations optimization (brown) and expectations obtained by posterior perturbations from deterministic optimization results (orange) where uncertainties come from geometrical and materials’ properties perturbations (U_{g}, U_{m}). Red and pink sets of points come from Figure 9. Dark gray zones highlight points with similar mean torque expectation values. 
Fig. 14 Pareto fronts obtained by worstcases optimization (darkcyan) and worstcase obtained by posterior perturbations from deterministic optimization (cyan) where uncertainties come from geometrical and materials’ properties perturbations (U_{g}, U_{m}). Green and light green sets of points come from Figure 10. Dark gray zones highlight points with similar mean torque worstcase values. 
Compared to results obtained considering only U_{g} as random variables, we can notice that the introduction of magnetic material uncertainties in the robust optimizations only affects the mean torque while the torque ripple stays almost the same. We can observe that the main effect of degradation of the material’s magnetic properties is a horizontal translation towards lower values of the mean torque compared to the solutions obtained when we only considered U_{g} as random variables. In order to look deeper into this result, we present in Figure 15 the evolution of the estimated mean torque (metamodel) for machine in Figure 9 in respect to α and β. As it can be seen, the mean torque metamodel has a linear behavior with respect to these two properties, and its maximum value is found when no degradation is applied (α = 1 and β = 0). At α = 0 and β = 1 (full degradation), the torque is around 411 N∙m with a 20 N∙m decrease.
To confirm this observation all over the search space (X_{g}), a subset of 5000 points was analyzed (this subset was created using LHS maximin over X_{g} including as a reference to the machine related to zone and coming from the expectations Pareto front, for which the relationship is linear (as seen in Fig. 15). For all these points, the second partial derivative with respect to U_{m} was evaluated. The derivative results at some U_{m} values are given in Figure 16. Each boxplot in this figure represents the derivative values for all 5000 points at a fixed α and β (values given on the horizontal axis). We can notice little variation for each boxplot, which means that the second derivatives are very similar all over the X_{g} search space. Moreover, the derivative values are all close to zero. This confirms that the mean torque can be considered as a linear function with respect to U_{m}:(8)
Fig. 16 Boxplots showing the variability (with respect to a batch of X_{g} samples) of the mean torque metamodel second partial derivative (with respect to U_{m}) evaluated at a given value of U_{m}. The blue horizontal lines stand for the machine related to zone and coming from the expectations Pareto front. 
Having shown that the metamodel constructed for the mean torque can be approximated to a linear function with respect to U_{m}, let us calculate the variability of the operator A(X_{g}) with respect to X_{g}. According to our linear approximation, the operator A(X_{g}) is equivalent to the partial derivative with respect to U_{m} of the mean torque’s metamodel. Figure 17 shows the results in the same manner as Figure 16, this time for the first derivative with respect to U_{m}. Once again, no significative variations are observed for all the boxplots meaning that the first partial derivative with respect to U_{m} is almost independent of X_{g}. The matrix A(X_{g}) can be considered as independent of X_{g}:(9)
Fig. 17 Boxplots showing the variability (with respect to a batch of X_{g} samples) of the mean torque metamodel first partial derivative (with respect to U_{m}) evaluated at a given value of U_{m}. The dark cyan horizontal lines stand for the machine related to zone and coming from the expectations Pareto front. 
In the same way, we observe that this first derivative is also almost independent of U_{m} (the derivative values are very close for all the boxplots). This is consistent with what was mentioned about the second derivative in Figure 16.
So far, we have shown the linear behavior of the mean torque (f_{1}) with respect to U_{m}. However, this does not explain the results obtained in Figures 13 and 14. According to our approximation, the expectation of f_{1} with respect to U_{m} and U_{g} as well as its expectation with respect to U_{g} alone are written as follows:(10) (11)with μ_{Um} the mean of U_{m} uncertain variables. And, for the worstcase:(12) (13)
We note that the approximate expression for the expectation of f_{1} with respect to U_{g} (11) where U_{m} is equal to the nominal material properties U_{nom} and the same expectation, but with respect to U_{g} and U_{m} (10) are the same except for terms that do not depend on X_{g}. The same is true for the worstcase expressions of f_{1} (12) and (13). This means that their maximum values will be for the same X_{g} but with shifted torque values (for the case of U_{g} and U_{m} with respect to the case of U_{g} only).
3.3 Results verification by FEM simulation
All the results presented in the previous sections were based on surrogate models. We will therefore verify in this section some results using FEM simulations. Figures 18 and 19 show the same boxplots based on the metamodels as in Figures 11 and 12. However, this time, FEM simulations were used to produce the same boxplot for comparison. Only two zones are represented in these figures : and . The STD values were also added for comparison.
Fig. 18 Boxplots showing comparisons of torque ripple between FEM simulations and surrogate model predictions: deterministic (blue) and robust (red) machines with similar predicted mean torque/mean torque expectation values. 
Fig. 19 Boxplots showing comparisons of torque ripple between FEM simulations and surrogate model predictions: deterministic (blue) and robust (green) machines with similar predicted mean torque worstcase values. 
The results from Figure 18 confirm what we have already noticed using the metamodels. For both zones and , the robust optimization presents similar or better solutions than the deterministic optimization in terms of robustness. For zone , both machines have practically the same FEMcomputed mean torque ripple with a smaller STD for the robust design. As for zone , the robust design presents a lower expected mean torque. The FEMcomputed q3 value of the robust design is also equal to the median of the deterministic one (5.8%), meaning that 75% of the produced machines using the robust design would have a torque ripple lower than 5.8% vs. only 50% for the deterministic design. Although the results in Figure 18 were obtained from the robust optimization based on the expectance formulation, the design in zone also present a particularly lower torque ripple in a worstcase scenario, with 8.4% vs. 9.9% for the deterministic design.
As seen in Section 3.1, using the optimization formulation with the worstcase scenario tends to reduce the STD of the robust designs compared to the deterministic ones, especially at low torque. This is also confirmed by FEM simulations in Figure 19. For the machines in zone , the robust design has a notably lower FEMcomputed STD compared to the deterministic one (0.4% vs. 0.78%). As for the worstcase torque ripple value, both designs possess similar FEMcomputed performances, around 7.3%. As for zone , similar results were obtained with deterministic and worstcase optimizations with a torque ripple of 11.8% and 11.6% respectively.
Regarding the precision of the metamodels compared to FEM simulations, some differences can be noticed. This lack of precision, especially for torque ripple, is somehow expected with a strategy using fixed metamodels in optimization, especially with the difficulty of fitting a torque ripple metamodel as seen in Section 2.1. A vast number of simulations is needed in this case to obtain an acceptable level of accuracy. An alternative approach could be to use an adaptative strategy to update the surrogate model with additional simulations during the optimization. The additional simulations will guide the algorithm towards the optimal zones while improving the metamodels’ accuracy. By using this kind of approach, additional FEM simulations are performed only in promising points in the search space, limiting computational time and increasing the precision for optimal designs. This will be the subject of a future study.
Finally, we note that surrogatebased optimizations require much less computational time compared to FEMbased optimizations. For example, 17 h were needed to finalize the DOE simulations (Intel(R) Xeon(R) W2195 CPU @ 2.30 GHz and 18 cores). The deterministic optimization using the metamodel took only 13 s for 300 iterations and 150 particles. Doing the same optimization using FEM simulations would have taken around 3 months to complete.
4 Conclusions
We have presented in this paper a comparison between different optimizations performed on a Permanent Magnet Assisted Synchronous Reluctance Machine. The first optimization used a deterministic formulation considering the design parameters as certain. In order to study the sensitivity of the machine to parameters uncertainties, robust optimizations were performed considering some geometric and magnetic parameter as uncertain. Two different formulations were adopted for the robust optimizations. In the first one, we used the expectations of the mean torque and torque ripple as objective functions. In the second formulation, the objective functions were defined as the worstcase values. In order to reduce computation time, surrogate models have been built for each of the objective functions. These surrogate models have been also used to perform a Sobol indicesbased sensitivity analysis to detect the most impacting input parameters. Objective functions’ expectations were computed with a quasiMonte Carlo scheme while worst cases were calculated with an embedded Particle Swarm Optimization algorithm. Both geometrical and magnetic property tolerances were shown to substantially influence the machines’ performances, with the magnetic property tolerances having a greater impact on mean torque. More precisely, we found that the magnetic material properties tolerances affect the mean torque linearly. We have shown that the considered magnetic material degradation has almost the same impact on the performances over all the search space: The Pareto front of the robust optimization considering both types of uncertainties (geometrical U_{g} and magnetic material properties U_{m}) is simply obtained by shifting horizontally (toward lower torque values) the Pareto front of the robust optimization with only the geometrical parameters (U_{g}) as uncertain. This result cannot be generalized since the search space in this paper has been restricted to a small region around a reference design.
It should also be noted that while the predicted values of mean torque were consistent with Finite Element Method simulations, some differences were observed for torque ripple. This problem will be addressed in future work by using an adaptive strategy, like Bayesian Optimization, to update the surrogate models with additional simulations during the optimization. Nevertheless, the Finite Element Method simulations have globally confirmed the trend of the predicted results using the metamodels. The comparison of Pareto fronts has shown that robust solutions outperform deterministic solutions in terms of different robustness criteria such as quantiles, expectation, worstcase, and standard deviation values. However, in a constrained optimization problem, it is recommended to use the worstcase as a constraint and not as an objective. In a real problem like a production factory, the worstcase, which is very conservative, is replaced by a quantile. Such formulation allows increasing the average performance of a batch of prototypes for example, while ensuring that only a small percentage of them would not respect the defined requirements.
The aforementioned aspects show the importance of developing and applying new optimization techniques to electrical machines when dealing with uncertainties.
Acknowledgments
We would like to thank the Mexican National Science and Technology Council CONAHCYT (Consejo Nacional de Humanidades Ciencias y Tecnologías) for its support through the Ph.D. scholarship 834161.
References
 Levent A.H., Lordoglu A., Aydeniz M.G. (2020) Design and optimization of permanent magnet synchronous motor for electric vehicle applications, in 2nd Global Power Energy and Communication Conference (GPECOM), Izmir, Turkey, pp. 148–151. https://doi.org/10.1109/GPECOM49333.2020.9247885. [Google Scholar]
 Choi G., Bramerdorfer G. (2022) Comprehensive design and analysis of an interior permanent magnet synchronous machine for lightduty passenger EVs, IEEE Access 10, 819–831. https://doi.org/10.1109/ACCESS.2021.3137897. [CrossRef] [Google Scholar]
 Balasubramanian L., Bhuiyan N.A., Javied A., Fahmy A.A., Belblidia F., Sienz J. (2023) Design and optimization of interior permanent magnet (IPM) motor for electric vehicle applications, CES Trans. Electr. Mach. Syst. 7, 2, 202–209. https://doi.org/10.30941/CESTEMS.2023.00021. [CrossRef] [Google Scholar]
 Zhu X., Xiang Z., Zhang C., Quan L., Du Y., Gu W. (2017) Coreduction of torque ripple for outer rotor fluxswitching PM motor using systematic multilevel design and control schemes, IEEE Trans. Ind. Electron. 64, 2, 1102–1112. [CrossRef] [Google Scholar]
 Sun X., Wan B., Lei G., Tian X., Guo Y., Zhu J. (2021) Multiobjective and multiphysics design optimization of a switched reluctance motor for electric vehicle applications, IEEE Trans. Energy Convers. 36, 4, 3294–3304. [CrossRef] [Google Scholar]
 Faltakh N., Hlioui S., Gabsi M., Benlamine R., Vangraefschèpe F. (2017) Design of a low permanent magnet mass hybrid excited flux switching machine for a PHEV application, in 15th International Conference on Electrical Machines, Drives and Power Systems (ELMA), Sofia, Bulgaria, pp. 271–276. https://doi.org/10.1109/ELMA.2017.7955447. [Google Scholar]
 Cisse K.M., Hlioui S., Belhadi M., Mermaz Rollet G., Gabsi M., Cheng Y. (2021) Design optimization of multilayer permanent magnet synchronous machines for electric vehicle applications, Energies 14, 7116. https://doi.org/10.3390/en14217116. [CrossRef] [Google Scholar]
 Ren Z., Zhang D., Koh C.S. (2013) A new sensitivitybased reliability calculation algorithm in the optimal design of electromagnetic devices, J. Electr. Eng. Technol. 8, 2, 331–338. [CrossRef] [Google Scholar]
 Yang Y., Zhang C., Bramerdorfer G., Bianchi N., Qu J., Zhao J., Zhang S. (2022) A computationally efficient surrogate model based robust optimization for permanent magnet synchronous machines, IEEE Trans. Energy Convers. 37, 3, 1520–1532. [Google Scholar]
 Ren Z., Zhang D., Koh C.S. (2013) New reliabilitybased robust design optimization algorithms for electromagnetic devices utilizing worst case scenario approximation, IEEE Trans. Magn. 49, 5, 2137–2140. [CrossRef] [Google Scholar]
 Liang J., Mourelatos Z.P., Tu J. (2007) A singleloop method for reliabilitybased design optimization, J. Mech. Des. 129, 12, 1215–1224. [CrossRef] [Google Scholar]
 Lei G., Bramerdorfer G., Liu C., Guo Y., Zhu J. (2021) Robust design optimization of electrical machines: a comparative study and space reduction strategy, IEEE Trans. Energy Convers. 36, 1, 300–313. [CrossRef] [Google Scholar]
 Wu J., Zhu X., Fan D., Xiang Z., Xu L., Quan L. (2022) Robust optimization design for permanent magnet machine considering magnet material uncertainties, IEEE Trans. Magn. 58, 2, 1–7. [Google Scholar]
 Jang G.U., Kim C.W., Bae D., Cho Y., Lee J.J., Cho S. (2020) Reliabilitybased robust design optimization for torque ripple reduction considering manufacturing uncertainty of interior permanent magnet synchronous motor, J. Mech. Sci. Technol. 34, 3, 1249–1256. [CrossRef] [Google Scholar]
 Crozier R., Mueller M. (2016) A new MATLAB and octave interface to a popular magnetics finite element code, in Proc. – 2016 22nd International Conference on Electrical Machines ICEM, pp. 1251–1256. https://doi.org/10.1109/ICELMACH.2016.7732685. [Google Scholar]
 Pronzato L., Müller W.G. (2012) Design of computer experiments: space filling and beyond, Stat. Comput. 22, 3, 681–701. [CrossRef] [MathSciNet] [Google Scholar]
 Park H.J., Yeo H.K., Jung S.Y., Chung T.K., Ro J.S., Jung H.K. (2018) A robust multimodal optimization algorithm based on a subdivision surrogate model and an improved sampling method, IEEE Trans. Magn. 54, 3, 1–4. [Google Scholar]
 Sun X., Wan B., Lei G., Tian X., Guo Y., Zhu J. (2021) Multiobjective and multiphysics design optimization of a switched reluctance motor for electric vehicle applications, IEEE Trans. Energy Convers. 36, 4, 3294–3304. [CrossRef] [Google Scholar]
 Deb K., Pratap A., Agarwal S., Meyarivan T. (2002) A fast and elitist multiobjective genetic algorithm: a fast and elitist multiobjective genetic algorithm: NSGAII, IEEE Trans. Evol. Comput. 6, 2, 182–197. [CrossRef] [Google Scholar]
 Lei G., Bramerdorfer G., Ma B., Guo Y., Zhu J. (2021) Robust design optimization of electrical machines: multiobjective approach, IEEE Trans. Energy Convers. 36, 1, 375–401. [Google Scholar]
 Kennedy J., Eberhart R. (1995) Particle swarm optimization, in Proceedings of ICNN’95 – International Conference on Neural Networks, Perth, WA, Australia, 27 Nov.–1 Dec. 1995, pp. 1942–1948. [Google Scholar]
All Tables
All Figures
Fig. 1 Geometry of the PMaSynRelM motor. 

In the text 
Fig. 2 Design parameters for one layer (# is the number of the layer). 

In the text 
Fig. 3 Nominal and degraded PM characteristic. 

In the text 
Fig. 4 Nominal and degraded B(H) curve. 

In the text 
Fig. 5 Optimization workflow. 

In the text 
Fig. 6 Mean torque Sobol’ indices. 

In the text 
Fig. 7 Torque ripple Sobol’ indices. 

In the text 
Fig. 8 Example of a torque ripple PDF. 

In the text 
Fig. 9 Pareto fronts obtained by deterministic (front in blue and expectations obtained by posterior perturbations in pink), and expectations optimization (red) where uncertainties come from geometrical perturbations U_{g}. Dark gray zones highlight points with similar mean torque expectation values. 

In the text 
Fig. 10 Pareto fronts obtained by deterministic (front in blue and worstcases obtained by posterior perturbations in light green), and worstcases optimization (green) where uncertainties come from geometrical perturbations U_{g}. Dark gray zones highlight points with similar mean torque worstcase values. 

In the text 
Fig. 11 Boxplots showing comparisons of predicted torque ripple values between deterministic (blue) and robust (red) machines from Figure 9 with similar predicted mean torque expectation values. 

In the text 
Fig. 12 Boxplots showing comparisons of predicted torque ripple values between deterministic (blue) and robust (green) machines from Figure 10 with similar predicted mean torque worstcase values. 

In the text 
Fig. 13 Pareto fronts obtained by expectations optimization (brown) and expectations obtained by posterior perturbations from deterministic optimization results (orange) where uncertainties come from geometrical and materials’ properties perturbations (U_{g}, U_{m}). Red and pink sets of points come from Figure 9. Dark gray zones highlight points with similar mean torque expectation values. 

In the text 
Fig. 14 Pareto fronts obtained by worstcases optimization (darkcyan) and worstcase obtained by posterior perturbations from deterministic optimization (cyan) where uncertainties come from geometrical and materials’ properties perturbations (U_{g}, U_{m}). Green and light green sets of points come from Figure 10. Dark gray zones highlight points with similar mean torque worstcase values. 

In the text 
Fig. 15 Mean torque as a function of α and β for the design in zone in Figure 9. 

In the text 
Fig. 16 Boxplots showing the variability (with respect to a batch of X_{g} samples) of the mean torque metamodel second partial derivative (with respect to U_{m}) evaluated at a given value of U_{m}. The blue horizontal lines stand for the machine related to zone and coming from the expectations Pareto front. 

In the text 
Fig. 17 Boxplots showing the variability (with respect to a batch of X_{g} samples) of the mean torque metamodel first partial derivative (with respect to U_{m}) evaluated at a given value of U_{m}. The dark cyan horizontal lines stand for the machine related to zone and coming from the expectations Pareto front. 

In the text 
Fig. 18 Boxplots showing comparisons of torque ripple between FEM simulations and surrogate model predictions: deterministic (blue) and robust (red) machines with similar predicted mean torque/mean torque expectation values. 

In the text 
Fig. 19 Boxplots showing comparisons of torque ripple between FEM simulations and surrogate model predictions: deterministic (blue) and robust (green) machines with similar predicted mean torque worstcase values. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.