Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 4526 (2024 ) Cite this article High Voltage PMSM
Advanced slot and winding designs are imperative to create future high performance electrical machines (EM). As a result, the development of methods to design and improve slot filling factor (SFF) has attracted considerable research. Recent developments in manufacturing processes, such as additive manufacturing and alternative materials, has also highlighted a need for novel high-fidelity design techniques to develop high performance complex geometries and topologies. This study therefore introduces a novel physics-informed machine learning (PIML) design optimization process for improving SFF in traction electrical machines used in electric vehicles. A maximum entropy sampling algorithm (MESA) is used to seed a physics-informed Bayesian optimization (PIBO) algorithm, where the target function and its approximations are produced by Gaussian processes (GP)s. The proposed PIBO-MESA is coupled with a 2D finite element model (FEM) to perform a GP-based surrogate and provide the first demonstration of the optimal combination of complex design variables for an electrical machine. Significant computational gains were achieved using the new PIBO-MESA approach, which is 45% faster than existing stochastic methods, such as the non-dominated sorting genetic algorithm II (NSGA-II). The FEM results confirm that the new design optimization process and keystone shaped wires lead to a higher SFF (i.e. by 20%) and electromagnetic improvements (e.g. maximum torque by 12%) with similar resistivity. The newly developed PIBO-MESA design optimization process therefore presents significant benefits in the design of high-performance electric machines, with reduced development time and costs.
Among all types of synchronous motors without electric current in their rotor, permanent magnet synchronous machines (PMSM)s gained significant consideration in the electric vehicle (EV) market recently1,2. This is due to their high performance, such as high torque at low operational speeds and efficient and reliable inverters. However, they are still far from ideal because of performance and manufacturing issues. From a performance perspective, there are risks associated with: (i) large braking torque when the traction inverter fails; (ii) maximum speed limitations caused by permanent magnets (PM) which can induce high voltages in stator windings; and (iii) no control over the current (i.e. power) in the rotor. From a manufacturing perspective, (i) the usage of rare-earth (RE)-based PMs is critical due to their environmental and sustainability concerns. Thereby enormous research interests are focused on PM-free motor topologies; (ii) expensive; and (iii) mechanically weak (e.g. breakable and damageable due to temperature rise and demagnetization). Researchers are nowadays working on many alternative alternating current (AC) motors which do not use PMs (e.g. reluctance motor (RM))3 or hybrid options with a lesser volume of PMs (e.g. PM-assisted synchronous reluctance motor4). So far, these are not competitive with PMSMs in electromagnetic performance5,6,7. A promising magnet-free and brushless design topology is AC electrically excited synchronous motors (EESM)s with an inductive power transfer system as presented in Fig. 1. The power transfer system consists of rotary and static coils to inject controllable direct current (DC) in the rotor windings. The EESMs were not being widely investigated yet. After the electrically switching excitation circuit invention8, the EESM is revived as a promising candidate for further development. Limited works highlight their capabilities9,10,11,12 for traction applications. The prime benefits of this topology are: (i) cost-effective, the total elimination of magnets which makes EESM much cheaper than other topologies using magnets (e.g. PMSMs), (ii) a more sustainable and environmentally friendly choice due to extremely high demands of PM materials and market limitations; (iii) excellent power factor which allows maximum torque improvement; (iv) increasing magnetic flux density by injecting more current in the rotor winding. The control over the rotor’s current and magnetic field helps with maximum torque density maximization; (v) increasing efficiency by minimizing electromagnetic loss, mainly during low torque operation; and (vii) safer operation, as the magnetic field in the rotor is controllable, allowing fast demagnetization to overcome hazardous braking torques. The brushed separately excited synchronous motor (SSM) topology was discovered more than a century ago, however they were not chosen for automotive applications because of their lower performance, e.g., torque production challenges. To date, many researchers are working on SSMs without brushes and slip rings for injecting electric current into the field (or rotor) windings. This problem was resolved by the researchers, more details are given in8. They invented a switching exciter circuit in the side of the rotor winding to transfer current from the power supply to the rotor of the EESMs. A few researchers investigated the analytical modelling of EESMs9, in which they studied a formalized method to calculate the field distribution and dq-axis model parameters. In another work10, the researchers proposed a low-order lumped parameter thermal network to predict the temperature rise of the EESMs, in which thermal resistances and capacities were identified individually. Stancu et al.11, studied a contactless rotary transformer and field converter for the excitation of the rotor winding. The researchers found that the EESM is capable of considerably more peak torque and power compared to high-performing interior permanent magnet synchronous motors (IPMSM)s. They were thermally limited by the rotor coils to approximately 30 s of performance at that level, but continuous power and torque also exceeded that of IPMSM. Maximum motor efficiency higher than 95%, including the losses of the motor, the rotating transformer, and the field converter was achieved. Overall, excellent efficiencies at high speeds were recorded, outperforming an equivalent optimized IPMSM with a similar size. At lower speeds, the IPMSM was more efficient. Jawad et al.12, introduced a new approach to excite the rotor winding of EESM. They utilized zero-sequence third-harmonic currents to generate an additional third-harmonic magneto-motive force (MMF) pulsating spatially in the stator winding. The proposed topology can be employed either as an inverter or a conventional AC power source to generate the harmonic component of the air-gap flux for rotor excitation and the fundamental component to produce the torque. In13, the authors presented a comprehensive review on low-cost high-performance ferrite magnets which are more sustainable and have a more stable supply chain for EV applications.
Traction electric drive system and additional components for EESMs.
In summary, Table 1 compares the most attractive traction motors qualitatively. Among these, the EESM is the most favourable and researchable topology for future EVs. Because of its inexpensive total material cost and unique power factor (close to one), indicating that most of the inverter power can be transferred to mechanical power at the shaft. The PMSM’s unique feature is performance, high torque at low speeds and premium efficiency at low and high speeds, due to its high power density capability5,6,7. However, they are significantly more expensive than other topologies because of PMs, in addition to their environmental and sustainability concerns. The low safety of PMSMs is mainly due to critical braking torque production if the inverter does not function properly (e.g. short-circuited). The induction motor (IM), without the usage of PMs and any mechanical contact with the moving parts, can be considered a safe and cheap topology. They suffer to provide a constant torque below the nominal speeds, in other words, the IMs fail to offer a constant maximum power over a wide speed range. Hence, they require a variable-speed transmission gearbox for EV applications14. The RMs are the cheapest topology, and yet unpopular because of their poor performance e.g., their maximum and continuous power capability is the lowest and torque density is also very low.
The use of design optimization in electrical machines accomplished significant performance improvements and developments. Many researchers used deterministic methods2,5,15,16,17,18, stochastic19 (e.g. genetic20,21,22 and particle swarm23,24,25 algorithms) methods, and machine learning (ML) methods1,26,27,28,29,30,31,32,33, such as Bayesian optimization (BO) algorithms31,34, to improve different aspects of motor performance. Most probabilistic ML methods offer a promising framework for understanding the uncertainty of the design problem and its performance predictions. Based on35,36, well-regulated BO is essential to gain the best possible outcomes. To find the best Surrogate model for the BO, a comparative study37 is conducted using the most popular surrogate models for calibration reasons, such as random forests (RF)38, deep ensembles (DE)39, Bayesian neural network (BNN)40, and Gaussian processes (GP)41. It is found that GPs can work well with BO-based design optimizations. Applied ML methods are recently utilized to improve different design aspects of electric machines for demagnetization faults42, bearing faults43, short circuit faults44, and control45. For example, Song et al.42 investigated a new meta-learning method under varying operating conditions.
The research hypothesis is: can a magnet-free electric motor provide similar performance to existing PM traction motors? This is under investigation by many researchers today. In this study, a recently invented EESM is electromagnetically investigated as a benchmark motor. A new physics-informed BO (PIBO) design optimization algorithm is developed to improve the slot fill factor (SFF) for the EESMs. The developed PIBO method dealt with the copper filling of stator trapezoidal slots using innovative keystone shaped Litz wires type 8 with as many as possible single insulated strands twisted and compressed into custom conductive profiles. The developed PIBO is driven by ML models, such as GP. Exclusively, a maximum entropy sampling algorithm (MESA) is also introduced for seeding the PIBO-MESA algorithm. The proposed PIBO-MESA could outperform some popular design optimization in terms of computation time by about 45%. The outcome of the proposed PIBO-MESA algorithm and use of keystone wires offer a 20% improvement in the SFF which also resulted in better electromagnetic performance. The developed PIBO-MESA can be used for different electric motors with trapezoidal slots used in their stators.
The EESM topology creates more design complications for the motor compared to other popular traction motors, such as PMSMs and IMs, because of the current injection in the rotor winding by a transformer, which consists of static and rotating parts. In the following sections, the main mathematical equations and design specifications of the studied EESM are given.
The impact of various rotor and stator coil turns is investigated to increase the SFF. The motivation is finding out the optimum design choices for desired inductance and torque improvement. Considering the electromagnetic torque as:
where maximizing the torque density simply depends on magnetic flux \(\Psi\) and power factor \({\text{cos}}\left(\varphi \right)\) for a supplied current \(I\) . Therefore, the torque can be increased if: (1) the resistance of stator and rotor windings are reduced with consideration of the same current baseline, and/or (2) the power factor is improved. The supplied currents based on magnetic flux, power factor, and torque can be calculated as a function of speed:
where dq-axis inductances (\({L}_{d}\) and \({L}_{q}\) ) are a function of dq-axis currents (\({I}_{d}\) and \({I}_{q}\) ) as well as rotor current \({I}_{r}\) . The rotor inductance \({L}_{r}{\prime}\) multiplied by \({I}_{r}\) provides stator flux.
The total (stator and rotor) conduction losses are given as:
The iron loss, based on15, with consideration of hysteresis and eddy current losses, is calculated using:
where \({k}_{h}\) is the hysteresis coefficient at \({f}_{0}\) and \({B}_{0}\) in W/kg, which depends on the core material type and thickness, \({k}_{e}\) is the eddy current coefficient at \({f}_{0}\) and \({B}_{0}\) in W/kg, mainly restricted by the lamination thickness, \(\alpha\) , \(\beta\) , and \(\gamma\) are coefficients (in range of 1–2) that depend on the properties of the iron core material, \(f\) and \({f}_{0}\) represent electrical and reference frequencies, \(B\) and \({B}_{0}\) are magnetic flux density and reference magnetic flux density, computed using a 2D FEA solver, \(\rho\) is volumetric mass density, \(V\) is an element volume, and \({k}_{u}\) is field factor that has a multiplicative effect.
Taking the motor pole number (p = 6) and seven windings in each slot being connected in parallel, this indicates that two windings are in each slot if they are connected in series. Considering the benefits of trapezoidal slot shapes to reduce the winding resistance and leakage inductance, the optimum size of each Litz wires bar is computed for the stator slot:
where h is the bar height, Nc is the number of conductors, B is the magnetic flux density, W is the bat’s width, σc is the conductivity of copper. As B is almost in a tangential direction, the use of fixed dimension bars (consisting of several conductors) is unsuitable. Additionally, the use of single bars has manufacturing challenges because of parallel connection requirements. Hence, the four identical, type eight Litz wires topology was selected as the benchmark model presented in Fig. 2. This Litz wire type contains single insulated magnetic wire strands twisted and flattened into a rectangular shape with optional outer insulation of textile yarn. This topology allows customizing the bar sizes to fill in the slots thoroughly, which offers excellent SFF15,16,17.
Design space feasibility study for the selected variables for lower and upper boundary selection, considering winding temperature, maximum torque, efficiency, and output power. Note that grey and black (variables’ combination) solutions are unfeasible and feasible design variable spaces.
To achieve superior electromagnetic performance (e.g. efficiency, peak torque, and continuous power), the maximization of SFF in the stator and rotor is vital. Particularly, the stator Joule loss (first term in Eq. 3) plays a considerable role in power losses at a wide range of speeds, and the stator resistance has an inverse proportional relationship to SFF. In other words, increased SFF leads to lower stator Joule loss, suppressing the motor’s temperature rise. For this reason, a new design optimization process is introduced, in section B, to further increase SFF and the EESM’s electromagnetic means.
The reference motor design specification is inspired by8; the stator and rotor are made of M270-35A steel laminations with 0.95 stacking factor, dimensions are given in Table 2. The EESM uses distributed winding topology with two layers, in which the windings are limited to two bars inserted in each slot, each bar consists of 65 copper wires with a diameter of 0.8 mm for the stator core. In Fig. 2, normalized available design spaces are presented in grey, only black dots are considered feasible design solutions. The feasible design (variable) spaces satisfied the following constraints: (i) maximum winding temperature < 150 °C, (ii) maximum torque > 250 N m, (iii) efficiency at 7000 rpm > 90%, and (iv) output power > 80 kW for the same motor size. Design space individuals in grey have not satisfied at least one of the above constraints alongside the tooth width and yoke height for avoiding saturation and overheating in the motor. The design variables, as given in Table 3, and their design spaces are studied to find feasible solutions (i.e., variable combinations) for lower/upper boundary selections. The studied motor topology can offer excellent power over a wide range of speeds compared to other attractive traction motors as shown in Fig. 3.
Design optimization parameters, (a) the EESM 2-D model, and (b) the stator slot and four identical, type eight concentric Litz wires topology bars.
To maximize the SFF of the EESM, several related changeable variables are chosen, as given in Table 2. A fast and high-fidelity (PIBO-MESA) optimization algorithm is employed to tailor the conductive bars with four identical, type eight Litz wires. In this design process technique, a randomized maximum entropy sampling algorithm (MESA) is employed to fill the design space of each design variable46,47 illustrated in Table 3 and Fig. 3. Design variable ranges are found using a feasibility analysis, in which the min/max boundaries of each variable are restricted for feasible solutions. The MESA techniques are also used in other applications18,19. For different engineering problems, they can be tied with other local and global searchers to find global optimum solutions20,21,22,23. In this work, the supervised Gaussian-based model is uniquely tied to the MESA using a machine learning algorithm to find the best design selection for the copper in the stator slots of the EESM. In this design problem, the stator trapezoidal slots are assumed to host two bars, thereby the challenge is to find the optimum number of conductors to fit both bars within the slot dimensions considering the current density and temperature rise distributions, as presented in Fig. 3. With the winding topology chosen for the best electromagnetic performance, the windings are manually changed between the two available options—two-layer or four-layer—for the simulations.
In Algorithm 1, the utilization of MESA introduces a new random approach to generate iterative samples for the PIBO in the search space. The objective of MESA is to select a most informative subset of s random design variables from a set of n random design variables, subject to side and/or logical constraints. The randomized variables are mostly assumed to be Gaussian, or that they can be properly transformed in many states. Assuming an optimal solution of \(\widehat{x}\in \left(0,\infty \right)\) , sampling of a s-sized subset \(S\subseteq [n]\) is generated with a probability of:
A GP-based PIBO-MESA optimization algorithm is a probabilistic approach which utilizes a procedure to iteratively fit the probabilistic surrogate model to measured values of an objective function. The use of BO is effective when the objective function is expensive35,36,37. In this work, the MESA sampling mechanism is performed using Algorithm 1 to improve the random input data distributions for a high-fidelity design optimization process, particularly in initialization practice. The implementation and assumptions of the design optimization process in Algorithm 2 are explained next.
In the initialization phase, the samples generated by MESA using Algorithm 1 in the search domains are simulated. Next, based on35, most GPs are initialized by uniform random methods using an instantaneous regret term \({{\text{r}}}_{t}\) to be minimized, see Eq. (7), a \(\frac{{\Lambda }_{0}}{2}\) term is used for multi-level fidelity design optimization processes. In this study, two-level fidelity is considered; in which the first-level process used \(\frac{{\Lambda }_{0}}{2}\) and the second-level process utilized \(\frac{{\Lambda }_{0}}{2}\) .
where \({m}_{t\ge 0}\) is requested at every iteration, \(N\) is the random quantity within all fidelities up to \(\Lambda\) . In this optimization, only when \(f={f}^{M}\) is queried; an instantaneous reward term \({y}_{t}\) is set \(-\infty\) if \({m}_{t}\ne M\) and \({f}^{M}({x}_{t})\) if \({m}_{t}=M\) . Homogenously, \({r}_{t}={f}_{op}-{y}_{t}\) indicates \({r}_{t}=+\infty\) , when \({m}_{t}\ne M\) and \({f}_{op}-{f}^{M}({x}_{t})\) if \({m}_{t}=M\) . The minimum regret \(S\left(\Lambda \right)\) could be achieved using the MESA.
In the experimentation phase, based one35,36,37, the squared exponential (SE) type of kernel is employed. The SE kernel is initialized via maximizing the GPs marginal likelihood on the early sample filling in the SE kernel at every 20 iterations applying marginal likelihood.
The objective was to maximize the SFF subjected to several constraints, such as frequency, skin effect, current density, temperature, and saturation. More details are given in Algorithm 2. For each stator slot, the SFF is improved when the cost function is maximized using the proposed design optimization method, where the cost is:
where j and n are the number of bars and the maximum number of Litz wire bars in each stator slot, i is the number of conductors per bar, and nn is the maximum number of conductors. In this study, n is the total number of bars up to maximum number of n bars (nmax = 200). Of course, the maximization of (6) is restricted to (i) design variable search spaces, and (ii) equal and non-equal performance-related constraints given in Algorithm 2. In this algorithm, a high-fidelity working mechanism of PIBO working with GP is presented. A squared exponential (SE) kernel is employed as they are well-known kernels to work with GP algorithms. Because of their universal nature and compatibility to be integrated with most types of functions. Each function in its prior has many derivatives and two main parameters to compute: (1) the length scale ℓ that defines the length of the fluctuations, in the function, extrapolated within the threshold ℓ of the model. (2) The output variance σ2 (or scale factor) which decides the average distance of the function away from its mean value. The seeding samples in the four- dimensional (4D) space are generated by Algorithm 1, to be searched iteratively and globally. The global search aimed at maximizing the cost-function (line 8) subjected to the constraint functions g given in lines 9–14. \({g}_{1}^{\left({m}_{t},i,j\right)}\left({x}_{t,i,j}\right)\) is required cross-sectional area for the maximum current density \({J}_{s}\) , where a is the number of parallel pathways in the winding, ph is the number of phases in the stator winding, \(V\) is terminal voltage, \(\eta\) indicates efficiency, and \(cos\theta\) is the power factor. \({g}_{2}^{\left({m}_{t},i,j\right)}\left({x}_{t,i,j}\right)\) considers the saturation by computing the fluxleakage corresponding with the field weakening in the stator tooth, where \({\widehat{H}}_{t,i,j}\) is the field intensity in the iron tooth, \({A}_{s}\left({x}_{t,i,j}\right)\) is the slot area without the iron part, \({A}_{s,fe}\left({x}_{t,i,j}\right)\) is the iron part of the stator tooth region. \({g}_{3}^{\left({m}_{t},i,j\right)}\left({x}_{t,i,j}\right)\) computes the saturation factor, where \({\widehat{U}}_{s}\left({x}_{t,i,j}\right)\) , \({\widehat{U}}_{r}\left({x}_{t,i,j}\right)\) , and \({\widehat{U}}_{g}\left({x}_{t,i,j}\right)\) are the magnetic potentials in the stator, rotor, and airgap regions in the EESM, respectively.
Algorithm 1 MESA for Given Design Variables.
Algorithm 2 High-fidelity PIBO-GP algorithm for slot filling.
\({g}_{4}^{\left({m}_{t},i,j\right)}\left({x}_{t,i,j}\right)\) takes into consideration the induced voltage in the conductors. \({g}_{5}^{\left({m}_{t},i,j\right)}\left({x}_{t,i,j}\right)\) calculates the parasitic induced current induced in the conductor based on the main Eddy current loop, in which the current passes across the conductor. The parasitic induced current flows by delay of 90° compared to the current passing \(160 A\) at 12,000 rpm. All queries with x and y values updated the matrix \({D}_{t}^{({m}_{t})}\) , among which the best x combinations are stored in \({{\text{X}}}_{T}^{*}\) .
The working mechanism of the PIBO with the MESA is presented in Fig. 4, in which the start signal introduces the design variables available to influence the cost-function given in Algorithm 2. In the first part of the flowchart, the filling design points are distributed within the search bounds to fill the 4D cubes as cost-effectively as possible, as shown in Fig. 5. As the aim is to reduce the computational time, with assistance of Algorithm 2 (i.e. F1), the sampling distribution shrinks iteratively toward the global optimum region. Faster design optimization is achieved because both Algorithms 1 and 2 are physics-driven and well-calibrated. At the heart of the design optimization process, a 2D FEM solver is employed to simultaneously compute the EESM model. The FEM solver takes about 2 s for every new model simulation. When, the EESM design variables are manipulated to successfully maximize the SFF considering the constraints and search bounds, the design optimization process is ended.
The working mechanism of the proposed PIBO-MSEA design optimization process.
4D improved design optimization sample generation using PIBO-MESA (Algorithms 1 and 2) considering sample’s mass density in the permitted design spaces for (a) x1,x2,and x3, (b) x1,x2,and x4, (c) x1,x2,and x5, (d) x1,x2,and x6, (e) x1,x2,and x7, (f) x1,x2,and x8, (g) x1,x2,and x9, and (h) x1,x2,and x10.
The Bayesian-based design optimization process drives iteratively to fit a probabilistic surrogate model to stated values of a cost function. Therefore, a physics-informed feature can administer the future queries to explore the best areas within the search space to: (i) overcome the uncertainty of the surrogate model, and (ii) minimize the cost function much faster than other population-based methods. The developed PIBO-MESA uses probabilistic ML methods to offer a promising framework for understanding the uncertainty of the design problem and its performance predictions. The BO algorithm’s uncertainty estimation is calibrated based on the calibrator in41, in which the calibrated approach found the global minimum prior to the uncalibrated method, as shown in Fig. 6, for the cost function. Table 4 presents the internal parameters and hyperparameters chosen for the studied algorithms in this paper. For the hyperparameter tuning of neural networks, a grid of the hyperparameters epochs of ([1, 10], 1), dropout rate of ([0, 0.9], 0.08), hidden layer of ([1, 300], 30), and batch size of ([8, 256], 32), and learning rate of [10–5, 0.9] are selected. In Fig. 7a, the regret minimization using PIBO-MESA for these surrogate models is demonstrated. The GP method has shown the best performance on average. The physics-informed design process administers the superior search regions, as shown in Fig. 7b. It shows how the initial 4D filled search regions are shrunk iteratively from green to red, and purple at the end of the optimization. For the design optimization problem defined in section II, the performance of the most popular surrogate models is given in Table 5. In this table, the calibration outcomes of the PIBO-MESA using different ML surrogate models is reported. Considering the regret function, the best performance is reported using GPs and NSGA-II, whereas the GP model is 45.85% faster than the NSGA-II method.
Calibrating the PIBO-MESA using GP.
PIBO-MESA process evaluation, (a) the regret function minimization using different methods, and (b) shrinking iterative 4D search space.
The reported values are averaged by 20 different random initialization processes. The calibration error is ignored while regression is calculated for the recalibrated models, recalibration only is performed during BO. The regret on the objective function is given in the last iteration of the PIBO. Figure 8 indicates the copper SFF for the presented trapezoidal stator slot. AC winding loss can be produced in the copper when the injected current is alternating because of both skin and proximity effects. These loss mechanisms are especially high at high-speed operation, when the AC winding loss can become much higher than the DC winding loss. The typical two layers, four identical type 8 Litz wires benchmark model with 60% SFF is shown in Fig. 8a. The other two slots are optimized SFF with different settings using keystone shaped Litz 8 wires for the best packing density. In Fig. 8b, the maximum SFF of 80% is achieved using keystone shaped wires with two layers bunched Litz wire 8. Due to the benefits of more winding layers to reduce the AC and DC winding losses, a tradeoff should be done to select the best number of layers and SFF. The second optimized solution, as shown in Fig. 8c, is set to four layers; this yields a SFF of 75% but with more flexibility for different winding configurations and lower AC winding loss. The number of parallel paths is set to two, where the impedance for each parallel path is considered to be the same to prevent current unbalance and additional copper loss. The winding copper losses are examined with sinusoidal current excitations.
Slot filling using different methods, (a) benchmark design with four identical type 8 Litz wires, (b) two-layer type 8 Litz, keystone shaped wires winding using the PIBO-MESA with 80% SFF, (c) four-layer type 8 Litz keystone shaped wires winding using the PIBO-MESA with 75% SFF.
The electromagnetic results are presented to verify the outcome of the design optimization using PIBO-MESA, in which GP is chosen as the surrogate model. The average optimal variables after 10 runs are \({x}_{opt}=\{\mathrm{6.52,4.14,22.59,21.55,85.41,122.35,93.75,93}.\mathrm{47,4.52,9.21}\}\) . Figure 9 demonstrates the no-load magnetic flux density distribution of the models. The simulations are done under 200 A and 6000 rpm. For all three models, the maximum magnetic field intensity is observed at the bottom of the stator teeth near the slot opening. In Fig. 9a, the model 1 (M1) with two-layer (L = 2), four identical, type 8 Litz wires and SFF of 60% shows a maximum magnetic intensity of 1.75 T. Using similar number of layers, the SFF is improved by 20% in the second model (M2), as presented in Fig. 9b, where 1.88 T is seen in the same spots. To reduce the AC copper loss and efficiency capabilities at high speeds, a greater number of winding layers (L = 4) is used, however that reduced the SFF by 5%. The third model’s maximum magnetic intensity is 1.81 T at the same spots.
No load magnetic flux density distribution, (a) M1: benchmark motor, (b) M2: two-layer keystone shaped wound motor using the PIBO-MESA with 80% SFF, (c) M3: four-layer keystone shaped wound motor using the PIBO-MESA with 75% SFF.
Table 6 presents the AC and DC losses for the studied EESM motors. The results show that both optimized solutions (M2 and M3) provided a lower copper loss and higher efficiency during both low and high speed operations. Note that the only difference between M2 (with 80% SFF) and M3 (with 75% SFF) is the number of winding layers, in which M2 and M3 have two and four layers, respectively. The M2 motor can achieve a better performance during the low speeds, whereas the M3 motor provided lower AC copper loss at high speed operations. After a trade-off between AC and DC winding losses, the optimized M2 motor is shown the best performance.
Figure 10 illustrates the magnetic iron loss density, i.e. hysteresis and Eddy-current losses, influenced by the new designs with higher SFF. All models show the highest iron losses at the bottom of the stator teeth and rotor outer surface close to the stator core. Figure 10a1–c1 present the hysteresis loss density, in which the maximum hysteresis loss densities of M1, M2, and M3 motors are 435 kW/m3, 571 kW/m3, and 513 kW/m3, respectively. The Eddy-current loss densities, as presented in Fig. 10a2–c2, in which the maximum Eddy-current loss densities are 711 kW/m3, 1195 kW/m3, and 1058 kW/m3, respectively.
No load iron loss density distribution, (a) M1: benchmark motor, (b) M2: two-layer wound motor using the PIBO-MESA with 80% SFF, (c) M3: four-layer wound motor using the PIBO-MESA with 75% SFF.
Figure 11 demonstrates the peak torque and power for the studied models. The results are obtained under a DC-link voltage of 400 V. The maximum electromagnetic torque of M1 with 280 Nm is improved due to higher SFF. As presented in Fig. 11a, the optimized M2 and M3 motors enhanced the maximum torque by 11.78% and 9.28%, respectively. In Fig. 11b, the maximum power of M1 motor with 175 kW is also slightly increased, with the M2 and M3 motors offering 1.14% and 2.28% power increase, respectively.
Comparative peak torque and power profiles, (a) maximum torque, and (b) the maximum power in a wide speed range.
Figure 12 evinces the efficiency maps of all the studied motors under similar conditions for a fair comparison. The efficiency maps are produced at a 400 V DC-link voltage without consideration of inverters power loss. Figure 11a shows the efficiency map of the benchmark motor, in which the highest efficiency is 93.6% achieved at low torque and medium speed (acceleration phase). The M2 motor, as shown in Fig. 11b, achieves a higher efficiency (94.9%) and a larger premium efficiency region. The efficiency is enhanced by approximately 1.39% in the premium efficiency region. Figure 11c indicates the efficiency map of the M3 motor with a peak efficiency of 94.2%. The M3 motor increased the efficiency by about 0.64% compared to M1. The efficiency is improved the most in M2 across a wide range speed. Whereas the M3 sacrificed some SFF in order to have lower AC losses, but as presented in Table 5 this is cancelled by higher DC losses.
Efficiency maps for the studied EESMs, (a) M1: benchmark motor, (b) M2: two-layer, keystone shaped, wound motor using the PIBO-MESA with 80% SFF, (c) M3: four-layer, keystone shaped, wound motor using the PIBO-MESA with 75% SFF.
Increasing SFF eases the amount of air inside the slot which is being replaced by copper. Air is a weak thermal conductor, hence, a decrease in the amount of air in the slot might substantially improve the thermal conduction from the conductors to the cooling liquid ducts. The improved EESMs (M2 and M3) are within the allowed temperature rise, as presented in Fig. 13. The slot area (including winding) temperature computed using the same method as in44, the results show the case where a peak phase current of 240 A is injected in the stator windings. The steady-state temperature of M2 and M3 motors is decreased by 13.4% and 9.5%, respectively.
Transient winding temperature for the studied EESMs.
The electromagnetic contribution of the stator copper winding is vital in most traction electrical machines, such as permanent magnet synchronous and induction machines, used in electric vehicles. One of the promising traction electrical machines for low-cost and high-performance applications is magnet-free brushless electrically-excited synchronous machines (EESM). In this work, a new physics-informed Bayesian design optimization (PIBO) method for improving the slot filling factor (SFF) in such machines is presented. In this design optimization, a maximum entropy sampling algorithm (MESA) is used to seed an iterative PIBO algorithm, where the target function and its approximations are produced by Gaussian processes (GPs). The proposed PIBO-GP-MESA worked exclusively with a FEM solver to perform the GP surrogate to achieve an optimal combination of design variables. Significant computational gains were achieved using the new PIBO-MESA approach, which was 45% faster than existing stochastic methods, such as the non-dominated sorting genetic algorithm II (NSGA-II). The proposed PIBO-MESA was coupled with a 2D finite element model (FEM) to perform a GP-based surrogate and provide the first demonstration of the optimal combination of design variables. The FEM results confirm that the new design optimization process tends to a new approach to stator slot winding filling using keystone shaped wires, leading to a higher SFF (i.e. by 20%) and electromagnetic improvements (e.g. maximum torque by 12%) with similar resistivity. The proposed motor achieved an SFF of 80%, and its maximum power is also increased by 2.28%. In absolute terms, the PIBO-MESA design optimization process presented here represents a significant pathway in the faster design and production of future high-performance electrical machines.
The generated data is partially available upon request from the authors. Data requests can be made via this email: pedram.asef@ucl.ac.uk.
Asef, P., Denai, M., Paulides, J. J. H., Marques, B. R. & Lapthorn, A. A novel multi-criteria local Latin hypercube refinement system for commutation angle improvement in IPMSMs. IEEE Trans. Ind. App. 59(2), 1588–1602 (2022).
Asef, P., Bargallo, R., Moazami, S. & Lapthorn, A. Rotor shape multi-level design optimization for double-stator permanent magnet synchronous motors. IEEE Trans. Energy Conv. 34(3), 1223–1231 (2019).
Mishra, S., Varshney, A., Singh, B. & Parveen, H. Driving-cycle-based modeling and control of solar-battery-fed reluctance synchronous motor drive for light electric vehicle with energy regeneration. IEEE Trans. Ind. App. 58(5), 6666–6675 (2022).
Pescetto, P. & Pellegrino, G. Determination of PM flux linkage based on minimum saliency tracking for PM-SyR machines without rotor movement. IEEE Trans. Ind. App. 56(5), 4924–4933 (2020).
Schweidtmann, A. M. & Mitsos, A. Deterministic global optimization with artificial neural networks embedded. J. Optim. Theory Appl. 180, 925–948 (2019).
Chen, H. & Demerdash, N. A. O. Investigation of a 3D-Magnetic Flux PMSM with high torque density for electric vehicles. IEEE Trans. Energy Conv. 37(2), 1442–1454 (2022).
Asef, P. et al. Assessment of the energy consumption and drivability performance of an IPMSM-driven electric vehicle using different buried magnet arrangements. Energies 14, 1–22 (2021).
Illiano, E. Synchronous machine with switching element in the excitation circuit. US Patent, no. US20130193903A1 (2015).
Laldin, O., Sudhoff, S. D. & Pekarek, S. An analytical design model for wound rotor synchronous machines. IEEE Trans. Energy Conv. 30(4), 1299–1309 (2015).
Wang, E., Grabherr, P., Wieske, P., & Doppelbauer, M. A low-order lumped parameter thermal network of electrically excited synchronous motor for critical temperature estimation. In 2022 International Conference on Electrical Machines (ICEM), Valencia, Spain.
Stancu, C., Ward, T., Rahman, K. M., Dawsey, R. & Savagian, P. Separately excited synchronous motor with rotary transformer for hybrid vehicle application. IEEE Trans. Ind. App. 54(1), 223–232 (2018).
Jawad, G., Ali, Q., Lipo, T. A. & Kwon, B. Novel brushless wound rotor synchronous machine with zero-sequence third-harmonic field excitation. IEEE Trans. Mgn. 52(7), 1–4 (2016).
Luk, P.C.-K., Abdulrahem, H. A. & Xia, B. “Low-cost high-performance ferrite permanent magnet machines in EV applications: A comprehensive review. Transportation 6, 1–13 (2020).
Mun, J. et al. Reliability-based design optimization of a permanent magnet motor under manufacturing tolerance and temperature fluctuation. IEEE Trans. Mgn. 57(6), 1–4 (2021).
Asef, P., Bargallo Perpina, R., Barzegaran, M. R., Lapthorn, A. & Mewes, D. Multiobjective design optimization using dual-level response surface methodology and booth’s algorithm for permanent magnet synchronous generators. IEEE Trans. Energy Conv. 33(2), 652–659 (2018).
Lei, G., Bramerdorfer, G., Ma, B., Guo, Y. & Zhu, J. Robust design optimization of electrical machines: multi-objective approach. IEEE Trans. Energy Conv. 36(1), 390–401 (2021).
Wu, J. et al. Robust optimization of a rare-earth-reduced high-torque-density Pm motor for electric vehicles based on parameter sensitivity region. IEEE Trans. Veh. Tech. 71(10), 10269–10279 (2022).
Zhao, Y., Chen, W. & Yang, X. Adaptive sampling stochastic multigradient algorithm for stochastic multiobjective optimization. J. Optim. Theory Appl. https://doi.org/10.1007/s10957-023-02334-w (2023).
Vatani, M., Ghaffarpour, A., Jalali Kondelaji, M. A. & Mirsalim, M. Study of a modular toothed linear hybrid reluctance motor with permanent magnets in translator slots. IEEE Trans. Transport. Elec. 8(3), 3554–3567 (2022).
Anvari, B., Toliyat, H. A. & Fahimi, B. Simultaneous optimization of geometry and firing angles for in-wheel switched reluctance motor drive. IEEE Trans. Transport. Elec. 4(1), 322–329 (2018).
Sun, X., Xu, N. & Yao, M. Sequential subspace optimization design of a dual three-phase permanent magnet synchronous hub motor based on NSGA III. IEEE Trans. Transport. Elec. 9(1), 622–630 (2023).
Ma, C. & Qu, L. Multiobjective optimization of switched reluctance motors based on design of experiments and particle swarm optimization. IEEE Trans. Trans. Energy Conv. 30(3), 1144–1153 (2015).
Yang, Z., Lu, C., Sun, X., Ji, J. & Ding, Q. Study on active disturbance rejection control of a bearingless induction motor based on an improved particle swarm optimization-genetic algorithm. IEEE Trans. Transport. Elec. 7(2), 694–705 (2021).
Dong, F., Zhao, J., Song, J., Feng, Y. & He, Z. Optimal design of permanent magnet linear synchronous motors at multispeed based on particle swarm optimization combined with SN ratio method. IEEE Trans. Trans. Energy Conv. 33(4), 1943–1954 (2018).
Dong, F., Zhao, J., Song, J., Feng, Y. & He, Z. Zhao-Hua Liu, Xiao-Hua Li, Liang-Hong Wu, Shao-Wu Zhou, and Kan Liu. IEEE Trans . Ind. Inf. 11(5), 1220–1230 (2015).
Song, J. et al. An efficient multiobjective design optimization method for a PMSLM based on an extreme learning machine. IEEE Trans. Ind. Elec. 66(2), 1001–1011 (2019).
Shimizu, Y., Morimoto, S., Sanada, M. & Inoue, Y. Automatic design system with generative adversarial network and convolutional neural network for optimization design of interior permanent magnet synchronous motor. IEEE Trans. Trans. Energy Conv. 38(1), 724–734 (2023).
Asanuma, J., Doi, S. & Igarashi, H. Transfer learning through deep learning: Application to topology optimization of electric motor. IEEE Trans. Mgn. 56(3), 1–4 (2020).
Omar, M. et al. Review of machine learning applications to the modeling and design optimization of switched reluctance motors. IEEE Access 10, 130444–130468 (2022).
Kirchgässner, W., Wallscheid, O., & Böcker, J. Data-driven permanent magnet temperature estimation in synchronous motors with supervised machine learning: A benchmark. IEEE Trans. Energy Conv. 36(3), 2059–2067 (2021).
Lai, C., Feng, G., Mukherjee, K., Tjong, J. & Kar, N. C. Maximum torque per ampere control for IPMSM using gradient descent algorithm based on measured speed harmonics. IEEE Trans. Ind. Inf. 14(4), 1424–1435 (2018).
Dhulipati, H. et al. Advanced design optimization technique for torque profile improvement in six-phase PMSM using supervised machine learning for direct-drive EV. IEEE Trans. Energy Conv. 34(4), 2041–2051 (2019).
Song, Q., Wang, M., Lai, W. & Zhao, S. On Bayesian optimization-based residual CNN for estimation of inter-turn short circuit fault in PMSM. IEEE Trans. Pow. Elec. 38(2), 2456–2468 (2023).
Kandasamy, K., Dasarathy, G., Oliva, J., Schneidar, J. & Poczos, B. Multi-fidelity Gaussian process bandit optimisation. J. Art. Intel. Res. 66, 151–196 (2019).
Huang, X., Song, M., & Zhang, Z. Multi-objective multi-generation Gaussian process optimizer for design optimization. SLAC-PUB, arXiv:1907.00250v2, pp. 1–12 (2020).
Foldager, J., Jordahn, M., Hansen, L. K., & Anderson, M. R. On the role of model uncertainties in Baysian optimization. arXiv:2301.05983v1, pp. 1–14 (2023).
Chai, Z. & Zhao, C. Multiclass oblique random forests with dual-incremental learning capacity. IEEE Trans Neu. Net. Learn. Syst. 31(12), 5192–5203 (2020).
Hasegawa, T. & Kondo, K. Easy ensemble: simple deep ensemble learning for sensor-based human activity recognition. IEEE Inter. Things J. 10(6), 5506–5518 (2023).
Foong, A. Y. K., Burt, D. R., Li, Y., & Turner, R. E. On the expressiveness of approximate inference in Bayesian neural networks. Adv. Neu. Info. Proc. Sys. 33, pp. 1–12 (2020).
Deshpande, S., & Kuleshov, V. Calibrated uncertainty estimation improves Bayesian optimization. arXiv:2112.04620v3, pp. 1–25 (2023).
Song, J. et al. A new demagnetization fault recognition and classification method for DPMSLM. IEEE Trans. Ind. Inf. 16(3), 1559–1570 (2020).
Chen, J. et al. A meta-learning method for electric machine bearing fault diagnosis under varying working conditions with limited data. IEEE Trans. Ind. Inf. 19(3), 2552–2564 (2023).
Ehya, H., Skreien, T. N. & Nysveen, A. Intelligent data-driven diagnosis of incipient interturn short circuit fault in field winding of salient pole synchronous generators. IEEE Trans. Ind. Informatics 18(5), 3286–3294 (2022).
Zhang, M. et al. Vibration reduction controller for a switched reluctance machine based on HW/SW partitioning. IEEE Trans. On Ind. Informatics 17(6), 3879–3889 (2022).
Wu, Z., Sun, L., Zhan, W., Yang, C. & Tomizuka, M. Efficient sampling-based maximum entropy inverse reinforcement learning with application to autonomous driving. IEEE Robat. Auto. Lett. 5(4), 5355–5362 (2020).
Chen, Z., Fampa, M., & Lee, J. Generalized scaling for the constrained maximum-entropy sampling problem. arXiv:2302.04934v1, pp. 1–11, 2023.
Grunditz, E. A., Thiringer, T., Lindstrom, J., Lundmark, S. T. & Alatalo, M. Thermal capability of electric vehicle PMSM with different slot areas via thermal network analysis. eTransportation 8, 1–21 (2021).
Advanced Propulsion Laboratory (APL), Department of Mechanical Engineering, Faculty of Engineering Sciences, University College London, London, E20 3BS, UK
Department of Mechanical Engineering, Institute for Advanced Automotive Propulsion Systems (IAAPS), Faculty of Engineering and Design, University of Bath, Bath, BA2 7AY, UK
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Conceptualization: P.A., C.V.; methodology: P.A.; formal analysis and investigation: P.A.; writing—original draft preparation: P.A., C.V.; writing—review and editing: P.A., C.V.; funding acquisition: P.A., C.V.; resources: P.A., C.V.; supervision: P.A., C.V.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Asef, P., Vagg, C. A physics-informed Bayesian optimization method for rapid development of electrical machines. Sci Rep 14, 4526 (2024). https://doi.org/10.1038/s41598-024-54965-2
DOI: https://doi.org/10.1038/s41598-024-54965-2
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Mining Permanent Magnet Motor Sign up for the Nature Briefing: AI and Robotics newsletter — what matters in AI and robotics research, free to your inbox weekly.