1. Introduction
Electrical machines are of great interest thanks to their key role in systems electrification and their ever-increasing performance requirements for industrial and e-mobility applications. Among many types of electric machines, Synchronous Reluctance (SynRel) machines have been recently recognized to be a potential candidate for a number of applications, as they present cost-effective energy conversion without the use of permanent magnets, a wide speed range of operation and robust structure. Furthermore, Switched Reluctance Motors (SRM) present similar overall characteristics, but they do require a custom converter depending on the stator and rotor salient teeth combination and polarity. Considering their drive requirements, SynRel machines benefit from relatively cheaper converter topology and a more common control technique when compared to SRM, which is another popular motor type without permanent magnets that has been widely investigated.
Throughout the years, different methodologies have been adopted to design SynRel machines, focusing mainly on the rotor geometry. Early studies used analytical approaches to form flux barriers for maximum torque production and minimum torque fluctuation [
1,
2,
3]. These works provided a strong basis for further research where the finite element method (FEM) is adapted and used with computational algorithms to effectively optimize the flux barrier shape.
Traditional optimization of SynRel machines aims to find optimum values of parameters describing a specific rotor shape. Different flux barrier topologies, most commonly circular [
4], fluid-shaped [
5] and straight-segmented barrier structures [
6,
7], are used in the literature to design SynRel machines. In [
8], these flux barrier topologies are compared in terms of performance and computation time, showing that each of them presents advantages and disadvantages.
Parametric optimization, as described above, is proven to be a very effective approach for reaching the objectives in a reasonable amount of computation time. However, parametric optimizations naturally limit the exploration of the research space of possible geometries. TO aims to find the optimal structure with a set of materials without imposed shapes (template-based geometries), as is the case with parametric optimization where degrees of freedom are naturally limited.
Starting with the work of Dyck [
9] in 1996, TO found its place in electromagnetic design. Since then, many methodologies have been developed and combined together to design electrical machines. Among the different methods, On/Off, density and level-set methods are the most popular. The On/Off method, thanks to its simple theory and implementation, has been widely applied to the design of electric machines. In [
10,
11], a surface permanent magnet machine and a SynRel machine are optimized using a heuristic and stochastic approach, respectively. A modified version of this method, incorporating normalized Gaussian networks, has been used to design SynRel machines and internal permanent magnet (IPM) machines in [
12]. The same method is employed in [
13] for the optimization of rotor bars in a line-start IPM machine.
Another TO technique, the density method (DM), has also been investigated, as its implementation with sensitivity-based methods and gradient optimization algorithms lead to faster convergence. This method has been applied to wound field synchronous machines in [
14] and compared with the design achieved by parametric optimization. In a similar study presented in [
15], rotors of SynRel machines are optimized using a combined magneto-structural analysis. The same method is applied to two case studies in [
16], which include the minimization of the rotor mass for a surface permanent magnet machine and the minimization of torque ripple for a V-type IPM rotor.
Although in the literature there are a number of works using DM and exploring its adoption for the design of electric machines, there are just a few studies focusing on SynRel machines. Considering the complex anisotropy requirements to improve torque production of such machine type, an optimal design is even more challenging using TO compared to parametric optimization. Therefore, to better understand the potential advantages and drawbacks of using TO, for this type of electrical machine, a more detailed investigation is required. This paper aims to investigate the use of DM through the analysis of its optimization settings and their effects on the optimization results, which represents a novel aspect with respect to the state-of-the-art. The findings of this work are discussed in terms of the performance of the optimal solutions and computation time required for their achievement.
First, brief background information about DM is given in
Section 2. In
Section 3, design studies are shown with different optimization settings. Optimizations are performed on a benchmark machine to evaluate the effect of the mesh structure, initial density, penalization coefficient and optimization objectives. Results are then discussed in terms of performance, feasibility and computation time. Finally, TO is performed on an existing design to discuss its usage as a design refinement tool.
3. Investigation of Optimization Settings
A commercial finite element software, JMAG Designer, was used in this study. Mesh elements are used as cells that make up the design space, and they are given material properties so that optimization variables are created. Since the focus is SynRel machines, only iron and air were considered. Interpolation between these two materials was performed by using Equation (
1).
3.1. Details of the Benchmark Machine
A benchmark SynRel machine available and designed in former projects for light traction applications was considered. The specifications of this motor are summarized in
Table 1. The maximum current density was limited to 10 A/mm
, and a distributed single layer winding was used. A water jacket cooling system with a flow rate of 3 L/min was implemented. The lamination material was M290-50A for both the stator and rotor.
The benchmark machine is the result of a parametric multi-objective optimization design using a genetic algorithm with the objectives of maximizing the average torque and minimizing the torque ripple. It was selected as one of the optimum machines having a good compromise between average torque and torque ripple. The details of this optimization study are given in [
27], while the benchmark machine is shown in
Figure 2. The average torque of the optimal solution is 100.8 Nm with 7.3% ripple.
3.2. Optimization Starting from an Empty Design Space
In this part, TO was performed with different settings to optimize the material distribution in the rotor. An empty design space was used respecting the outer dimensions of the benchmark motor’s rotor. During the analyses, symmetrical poles were considered. In order to prevent and reduce the number of unfeasible results, an outer layer surrounding the rotor was used with a 0.5 mm thickness. The stator, winding configuration and peak current of the original machine was unchanged. The effect of the mesh resolution, starting material density, penalization factor and objective functions were investigated in terms of objective values and computation time. Later, an existing design space was used as the initial design space to perform a TO.
3.2.1. Effect of Discretization
The discretization of the design space directly affects TO in terms of computation time. A higher resolution can be obtained by increasing the total number of cells in the design region; however, the computation time is adversely affected. The sides of the mesh elements are limited to maximum lengths of 0.25, 0.5, 1 and 2 mm. The minimum and maximum mesh densities are shown in
Figure 3a,b, respectively.
TO was performed with the objective of maximizing the average torque with a penalization coefficient of 1 and initial density of 0.9. The number of iterations were limited to 300 for each model, even though the computation time was quite different due to the number of mesh elements. The optimization results are shown in
Figure 4. The models with a maximum mesh length of 1 and 2 mm reach their knee points (around 30th iteration) faster than the models with 0.25 and 0.5 mm mesh lengths (around 40th iteration). This might be insignificant in terms of iteration number but the model with 0.25 mm mesh takes 4.4 h to reach the knee point, whereas the one with 2 mm takes 9 min.
The maximum objective values were obtained as 97.5, 97.2, 94.7 and 93.2 Nm for decreasing mesh intensity. This suggests that with the increased mesh resolution, it is possible to obtain slightly better average torque by compromising the computation time. Quantitatively speaking, between the two extreme cases, there is 31 h of difference in computation time and 4.5% difference in average torque.
The material and flux density distributions of the final designs are given in
Figure 5 to show the convergence to a reluctance-like structure. Even though a bounding outer frame is used, models with 0.5 and 2 mm mesh do not have a feasible mechanical connection to the frame. This points out the need to add further manufacturing feasibility constraints to the optimization process or the need for a final refinement stage in order to make the optimal solution feasible. Nevertheless, analyzing all the optimal solutions shown in
Figure 5a–g, it is evident that there are negligible amounts of intermediate materials that do not represent any major challenges for shape definition.
3.2.2. Effect of Initial Density
The effect of initial density has been studied in [
21] for a magnetic shield, and only slight changes were observed in terms of objective value. Although the target problem is relatively complex in that case, similar results were obtained hereafter through optimizations conducted with initial densities of 0.2, 0.5 and 0.9 using the 1 mm mesh model. Slightly different average torque values of 94.3, 95.2 and 94.7 Nm were obtained for final designs with increasing starting density, as shown in
Figure 6. Analyzing the results, it is clear that starting from a design space with more air, the time required for convergence is greatly increased. This can be justified by the low overall sensitivities obtained in the early iterations.
The distribution of material densities is depicted in
Figure 7. Investigating the material distributions, it is seen that two models with starting densities of 0.2 and 0.5 have more intermediate materials compared to the initial density of 0.9.
Based on the material distributions and objective values, it is clear that the initial starting point has only a little effect on the optimization results. Considering that the optimal material distribution is affected by the complexity and nature of the problem, this observation should not be considered general.
3.2.3. Effect of Penalization Coefficient
Different types of interpolation schemes are used in the literature, and a summary can be found in [
26]. The exponential term in Equation (
1) is usually referred to as a penalization coefficient because a steep interpolation drives the material property towards higher permeability. In other words, a sharper change in the material properties occur.
To assess the dependency of the optimization results on the penalization parameter, a number of optimizations has been performed on the 1 mm mesh model. The average torque versus the iteration number is given in
Figure 8. The results indicate that the higher the penalization coefficient, the lower the convergence speed. In addition, there is a substantial drop in average torque when the penalization coefficient is higher than one.
When analyzing the material distributions, as shown in
Figure 9, it is noticeable that the number of intermediate materials are increasing with the penalization coefficient. Although the intermediate materials shown between red and blue have magnetic properties closer to air, it is hard to differentiate between material borders. It might be possible to further decrease the number of intermediate materials by increasing the number of iterations; however, it is evident that increasing the penalization coefficient worsens the convergence rate.
3.3. Optimization of Torque Ripple
The optimization of torque ripple is generally more challenging than optimizing average torque alone, as it actually requires a uniform distribution of flux in the air-gap. In DM, this requires the calculation of each cell’s sensitivity with respect to rotor position. Taking into account this complexity, the number of iterations for this study is increased to 1000. Using models with maximum mesh lengths of 1 and 2 mm, two optimizations are run with the objectives of increasing average torque and decreasing torque ripple, starting with an initial density of 0.9 and penalization coefficient of 1.
The best of both optimizations, considered as the one with the highest average torque, is shown in
Figure 10. The model with the 2 mm mesh has 92.3 Nm average torque with a 4.6% ripple, and the other model has 92.1 Nm average torque with a 6.3% ripple. One of the most important distinctions of torque ripple objectives is seen as the formation of an additional air segment, resembling a three barrier structure not seen in the previous sections. As a remark, despite the higher computational time, a more feasible magnetic structure is obtained.
When material distributions of two models, as shown in
Figure 10a,b, are compared, the model with the 1 mm mesh exhibits more intermediate materials than the model with the 2 mm mesh. This can be explained by the fact that the model had a higher number of variables and increased complexity of the objective functions, leading to slow convergence.
3.4. Optimization Starting from a Non-Empty Design Space
In the previous subsections, various optimizations have been performed on the benchmark machine starting from an empty design space, and the effect of some critical optimization settings have been investigated. In this part, the TO is applied to a non-empty design space, i.e., starting from an already optimized rotor structure. For this purpose, the rotor of the benchmark machine is used with its three flux barrier geometry. A single objective optimization of maximizing average torque is used, limiting the maximum number of iterations to 300.
The obtained material distribution and corresponding flux density distribution are shown in
Figure 11b,c. Electromagnetic torque versus output position is shown in
Figure 11a for both the original and optimized machine. The average torque of the optimized machine is slightly increased from 100.8 to 102.3 Nm (+1.5%), while the torque ripple reduces from 7.3% to 6.5%. From
Figure 11b, it can be clearly seen that the outermost and middle barriers have minor shape modifications while the innermost barrier is enlarged towards the shaft. This increase in the total ratio of air to iron thickness in the rotor is the reason for the torque increments due to better electromagnetic exploitation of the soft magnetic material. As a result, a minor performance increment is obtained by applying TO as a refinement tool. The underlying significance of this change in the performance is the fact that when using the parametric expression of the geometry, it is not possible to achieve such a shape, unless certain modifications are made to increase the geometrical degrees of freedom.
4. Conclusions
In this work, TO is applied to a SynRel machine using DM with particular focus on the rotor design. Effect of design space resolution, starting point, penalization and objective functions are investigated by performing TO starting from an empty rotor design space. These critical parameters and their effects on the material distribution, performance and computation time are investigated in depth.
Based on the results, it is apparent that using a high penalization coefficient to reduce the number of intermediate materials causes convergence problems for material distribution, which makes it hard to distinguish boundaries of material areas. Furthermore, any of the designs with intermediate materials were not advantageous in terms of performance.
In addition, it is found that the inclusion of torque ripple within the objective function, despite increasing the computation burden, helps to obtain a more uniform material distribution. The results show the importance of structural consideration, as in its absence, mechanical integrity is not guaranteed. The effect of discretization and initial densities has more of an effect on computation time than the objective value when compared to the other optimization settings. Nevertheless, starting TO with a medium number of cells is considered to be good practice when time is a constraint.
In the case where TO was applied starting from an already optimized rotor structure, a minor improvement is reached. The final design shows improved saliency through better electromagnetic exploitation of the rotor material, which is not possible to obtain using parametric-based optimizations.
The results of the reported design exercise suggest that the use of topology optimization can be used as a refinement design tool once the machine geometry has been largely defined and optimized with standard parametric-based approaches.