Multi-Objective Optimization of Yokeless and Segmented Armature Machine for In-Wheel Traction Applications Based on the Taguchi Method

: For electrical machines with complex structures, the design space of parameters can be large with high dimensions during optimization. Considering the calculation cost and time consumption, it is hard to optimize all the design parameters at the same time. Therefore, in that situation, sensitivity analysis of these design parameters is usually used to sort out crucial parameters. In this paper, the sensitivity analysis-based Taguchi method is applied to optimize the axial-flux permanent magnet (AFPM) machine with yokeless and segmented armature (YASA) topology for an in-wheel traction system. According to the key parameters and their sensitivity analysis, the optimal machine scheme to meet the performance requirements can be formed. In this case study, the machine performance is improved significantly after optimization. Lastly, the experimental results verify the accuracy of the model used in this study.


Introduction
With the development of electrical vehicles, in-wheel traction has recently become an important research hotspot [1,2].This application requires a small un-sprung weight, and the traction machine should meet the performance requirements in the limited installation volume [3,4].Compared with radial flux permanent magnet (RFPM) machines, axial-flux permanent magnet (AFPM) machines with yokeless and segmented armature (YASA) topology have the advantages of less mass and more efficiency, which are suitable for in-wheel traction application, especially in limited axial dimensions [5][6][7].
In order to obtain the optimal performance of the YASA machine, multi-objective optimization is essential and indispensable.The calculation of the AFPM machine requires a 3D solver [8,9].Traditional global optimization methods (such as genetic algorithm [10], particle swarm optimization [11], etc.) may not be suitable for a 3D model due to their calculation cost [12][13][14].Differently, the local optimization method like the Taguchi method reduces the calculation amount and improves the optimization efficiency through appropriate parameter selection [15,16].However, in the traditional Taguchi method, optimization parameters still varied over a wide range, and it could not be conducive when finding the optimal solution [17].
Some improved Taguchi methods have been proposed to effectively optimize electrical machines.Reference [17] proposed an improved fuzzy-based sequential Taguchi optimization method applied to an inter permanent magnet synchronous machine (PMSM) to in- crease the optimization accuracy by optimizing the structural parameters at multiple levels.To save computing time, reference [18] utilizes the fuzzy-based sequential Taguchi robust optimization method to improve the comprehensive performance of a five-phase PMSM.In [19], to address the inherent limitation of the Taguchi method in solving multi-response optimization problems, an improved regression rate methodology and a weighted factor multi-objective technique are incorporated to form a Taguchi method-based multi-objective design framework of line-start PMSM.Reference [16] proposed a new system-level sequential Taguchi method to achieve the optimal solution with high robustness for switched reluctance machines with the consideration of manufacturing tolerances.
However, the Taguchi method for YASA machines has not been studied yet.For the existing optimization algorithms used for YASA, the selected optimization parameters are usually not well selected and, thus, may not be appropriate.In this paper, the improved Taguchi method optimization model is proposed for the YASA machine.First, analysis of the sensitivity of the structural parameters to the optimization objectives should be conducted to find the main optimization variables, which have a high influence proportion on the optimization objective.Based on that, the optimal solution of the YASA machine can be obtained.
This paper explores the possibility of using the Taguchi method in YASA machine optimization, which is the main contribution of this paper.In addition, the proposed optimization model provides a feasible reference for other electric machines.
In this paper, the topology, design requirements of the in-wheel traction system, and initial design of the YASA machine are introduced in Section 2. In Section 3, the optimization process of the YASA machine using the improved Taguchi method is described.Section 4 analyzes the performances of the no-load and load of the optimized YASA machine.Based on the optimization results, a prototype is manufactured, and the performances are verified by experimental tests.Finally, the main findings are summarized in the conclusion in Section 6.

Initial Design
The topology of the YASA machine is shown in Figure 1, where the YASA machine consists of a yokeless and segmented stator and double external rotors.The yokeless and segmented stator combined with concentrated armature winding can effectively reduce the mass and improve its efficiency [1].And the main magnetic flux starts from the permanent magnet-N and passes through the stator core to the permanent magnet-S on the other side.After passing through the rotor core on the second side, the flux starts from the N-pole on the second side and passes through the stator core to the S-pole on the first side, and the main magnetic flux forms a closed loop [20].
Some improved Taguchi methods have been proposed to effectively optimize electrical machines.Reference [17] proposed an improved fuzzy-based sequential Taguchi optimization method applied to an inter permanent magnet synchronous machine (PMSM) to increase the optimization accuracy by optimizing the structural parameters at multiple levels.To save computing time, reference [18] utilizes the fuzzy-based sequential Taguchi robust optimization method to improve the comprehensive performance of a five-phase PMSM.In [19], to address the inherent limitation of the Taguchi method in solving multiresponse optimization problems, an improved regression rate methodology and a weighted factor multi-objective technique are incorporated to form a Taguchi methodbased multi-objective design framework of line-start PMSM.Reference [16] proposed a new system-level sequential Taguchi method to achieve the optimal solution with high robustness for switched reluctance machines with the consideration of manufacturing tolerances.
However, the Taguchi method for YASA machines has not been studied yet.For the existing optimization algorithms used for YASA, the selected optimization parameters are usually not well selected and, thus, may not be appropriate.In this paper, the improved Taguchi method optimization model is proposed for the YASA machine.First, analysis of the sensitivity of the structural parameters to the optimization objectives should be conducted to find the main optimization variables, which have a high influence proportion on the optimization objective.Based on that, the optimal solution of the YASA machine can be obtained.
This paper explores the possibility of using the Taguchi method in YASA machine optimization, which is the main contribution of this paper.In addition, the proposed optimization model provides a feasible reference for other electric machines.
In this paper, the topology, design requirements of the in-wheel traction system, and initial design of the YASA machine are introduced in Section 2. In Section 3, the optimization process of the YASA machine using the improved Taguchi method is described.Section 4 analyzes the performances of the no-load and load of the optimized YASA machine.Based on the optimization results, a prototype is manufactured, and the performances are verified by experimental tests.Finally, the main findings are summarized in the conclusion in Section 6.

Initial Design
The topology of the YASA machine is shown in Figure 1, where the YASA machine consists of a yokeless and segmented stator and double external rotors.The yokeless and segmented stator combined with concentrated armature winding can effectively reduce the mass and improve its efficiency [1].And the main magnetic flux starts from the permanent magnet-N and passes through the stator core to the permanent magnet-S on the other side.After passing through the rotor core on the second side, the flux starts from the N-pole on the second side and passes through the stator core to the S-pole on the first side, and the main magnetic flux forms a closed loop [20].The main design objectives of the 5 kW YASA machine applied on the in-wheel traction system are shown in Table 1.The design volume of the YASA machine is limited, especially the axial dimension.According to the performance requirements and size limitation, a 36 stator slots-32 rotor poles (36s/32p) YASA machine is initially designed by the traditional design method.The corresponding structural parameters are listed in Table 2, and the parameters are labeled in a partial model of the YASA, as seen in Figure 2. Some of the parameters originate empirically from previous similar projects and can be changed for other applications.The main design objectives of the 5 kW YASA machine applied on the in-wheel traction system are shown in Table 1.The design volume of the YASA machine is limited, especially the axial dimension.According to the performance requirements and size limitation, a 36 stator slots-32 rotor poles (36s/32p) YASA machine is initially designed by the traditional design method.The corresponding structural parameters are listed in Table 2, and the parameters are labeled in a partial model of the YASA, as seen in Figure 2. Some of the parameters originate empirically from previous similar projects and can be changed for other applications.However, using the traditional design method, it is hard to meet the multiple demanding performance requirements.Therefore, the initial design with the structural parameters of the YASA machine should be optimized using the Taguchi method in this pa- However, using the traditional design method, it is hard to meet the multiple demanding performance requirements.Therefore, the initial design with the structural parameters of the YASA machine should be optimized using the Taguchi method in this paper; see the following sections.

Multiple Objectives Optimization
In this section, the proposed optimization procedure using the Taguchi method is given.At the beginning, the type of in-wheel traction machine should be determined.Then, the optimization objectives/variables should be confirmed before giving the experimental matrix and FEA results.Based on the FEA and data analysis, the optimization levels of the variables will be determined.Furthermore, the optimal solution can be confirmed according to the sensitivity analysis.Figure 3 presents the optimization flowchart for the proposed Taguchi method, and the following subsections will discuss the steps in detail.
Then, the optimization objectives/variables should be confirmed before giving the experimental matrix and FEA results.Based on the FEA and data analysis, the optimization levels of the variables will be determined.Furthermore, the optimal solution can be confirmed according to the sensitivity analysis.Figure 3 presents the optimization flowchart for the proposed Taguchi method, and the following subsections will discuss the steps in detail.

Confirmation of Optimization Objectives and Variables
Considering the performance requirements of the in-wheel traction system, the output torque should be as large as possible within the volume and temperature limitation.The torque ripple of output torque should also be minimized to ensure the smoothness and stability of the YASA machine during operation [21].At the same time, to meet the requirements of heavy load operating conditions such as vehicle acceleration and climbing, etc., the overload torque drop rate of the YASA machine also needs to be reduced.Therefore, the average rated torque Tave, torque ripple coefficient krip, and overload torque drop rate ka are selected as optimization objectives.The torque ripple coefficient krip is the difference between the maximum and minimum output torque divided by the average value.The formula of krip can be expressed as: where Tmax and Tmin are the maximum and minimum values of rated torque, respectively; Tave is the rated torque of the YASA machine.The torque drop rate ka is related to the rated torque and overload torque, which is used to measure the overload capacity.ka is defined as follows where T5ave is the output torque under five times the rated current.

Confirmation of Optimization Objectives and Variables
Considering the performance requirements of the in-wheel traction system, the output torque should be as large as possible within the volume and temperature limitation.The torque ripple of output torque should also be minimized to ensure the smoothness and stability of the YASA machine during operation [21].At the same time, to meet the requirements of heavy load operating conditions such as vehicle acceleration and climbing, etc., the overload torque drop rate of the YASA machine also needs to be reduced.Therefore, the average rated torque T ave , torque ripple coefficient k rip , and overload torque drop rate k a are selected as optimization objectives.The torque ripple coefficient k rip is the difference between the maximum and minimum output torque divided by the average value.The formula of k rip can be expressed as: where T max and T min are the maximum and minimum values of rated torque, respectively; T ave is the rated torque of the YASA machine.The torque drop rate k a is related to the rated torque and overload torque, which is used to measure the overload capacity.k a is defined as follows where T 5ave is the output torque under five times the rated current.
A suitable selection of optimization variables is critical to the efficient machine optimization.When the inner and outer diameters are determined, the structural parameters of the YASA machine mainly include the rotor core thickness h r , permanent magnet thickness h m , pole arc coefficient α, air gap h g , pole shoot height h p , pole shoot width w p , slot height h s , and slot width w s .According to the design requirements, the axial length of the YASA machine should be less than 45 mm, which can be expressed as: Since the variation range of the air gap is small, it can be regarded as a fixed value, which is set as 0.9 mm here.The above formula can be written as: Then, there are three independent height variables left to be optimized, i.e., rotor core thickness h r , permanent magnet thickness h m , and pole shoot height h p .The slot height h s is a dependent variable, derived by Equation ( 4).In addition, to ensure these size parameters do not overlap with each other, the pole shoot width coefficient is defined as k p with a range of 0~0.5.Then, the w p can be expressed as: In summary, six optimization variables are listed in Table 3, where each parameter is divided into three levels.The selection of levels determined the number of experiments and the amount of calculation.For the above six optimization variables, the traditional design method needs to calculate 3 6 = 729 schemes, which is too much.In contrast, the Taguchi method can establish the experimental matrix where only 27 schemes are needed, which reduces the calculation amount by 96.3%.At the same time, the step size of each optimization variable can be reasonably selected, which makes the optimization result more accurate.

Experimental Matrix and FEA Results
The established experimental matrix is shown in Table 4.The advantage of the orthogonal experiment is that the most representative experiments can be selected from a large number of experiments to obtain reliable results, and the analysis and calculation are convenient [15].It can be seen from Table 4 that there are combinations of different optimization variable levels corresponding to the order of the experiments.And the results of the optimization objective are obtained by the finite element analysis (FEA), which will be used as the basis for parameter sensitivity analysis, as shown in Table 4.

Sensitivity Analysis of Optimization Variables
This section analyzes the influence of each optimization parameter on the FEA results.Firstly, the average value of the optimization objectives is analyzed by Equation ( 6), and the calculation results are shown in Table 5.
where n is the number of experiments, and S i is the FEA result of the i-th experiment.Each optimization variable corresponds to the average value of optimization objectives under different levels.For example, the average value of rated torque T ave of variable h m under level 2 is shown in Equation (7).Similarly, the average values of optimization objectives under different levels of other optimization variables can be calculated.The results are shown in Table 6.
where T ave (i) is the result of rated torque T ave in the i-th experiment, and m 2 hm (T ave ) is the average value of rated torque T ave about the optimization variable h m under level 2. Next, according to the results of the average values, the proportion of the optimization variables on the optimization objectives can be obtained using variance analysis SS.The calculation formula is as follows: where x represents the optimization variables.S stands for the optimization objective.m x i (S) is the average value of the optimization objectives S under level i. m(S) is the average value of the optimization objectives.The calculation results are shown in Table 7.According to the results of variance analysis, it is known that the rotor core thickness, pole shoot height, and slot width have a great influence on the rated torque T ave , among which the rotor core thickness accounts for the largest proportion, 68.03%.The pole arc coefficient and the pole shoot width coefficient have the greatest influence on the torque ripple coefficient k rip , accounting for more than 75%.The rotor core thickness, permanent magnet thickness, and pole shoot width coefficient have a great influence on the torque drop rate k a .
In the final stage, the parameter variables can be determined intuitively using the above results.To meet the requirements of high torque density, low torque ripple, and high overload capacity of in-wheel traction application, the final optimization variables are confirmed as h r (2), h m (2), α(2), h p (1), w s (2), k p (1).

Result of Optimal Variable
With the increase in rotor core thickness, the rated torque, five times overload torque, and torque drop rate decrease, and the decrease rate tends to be gentle, as shown in Figure 4a.In order to meet both the torque output and small torque drop rate, the rotor core thickness is taken as 4.8 mm.With the reduction in the permanent magnet thickness, the five times overload torque increases first and then decreases when the thickness is equal to 3.4 mm, reaching the peak.The rated torque and torque drop rate increase with the reduction in the permanent magnet thickness.As shown in Figure 4b, the permanent magnet thickness is taken as 3.4 mm.net thickness is taken as 3.4 mm.
The rated torque, five times overload torque, and torque drop rate show an upward trend with the increase in slot width, and the increasing trend of rated torque and five times overload torque tend to slow down.Although the output torque will increase with the slot width, the stator core is easier to saturate, which increases the iron loss and reduces the overload capacity.In order to obtain a small torque drop rate, the slot width is taken as 9.5 mm.With the increase in pole shoot height, the rated torque and five times overload torque decrease, but the torque decrease rate increases.As shown in Figure 4d, the pole shoot height is selected as 1.5 mm.The selection of the pole arc coefficient and pole shoot width coefficient is more complicated and needs to be determined according to the rated torque, five times overload The rated torque, five times overload torque, and torque drop rate show an upward trend with the increase in slot width, and the increasing trend of rated torque and five times overload torque tend to slow down.Although the output torque will increase with the slot width, the stator core is easier to saturate, which increases the iron loss and reduces the overload capacity.In order to obtain a small torque drop rate, the slot width is taken as 9.5 mm.With the increase in pole shoot height, the rated torque and five times overload torque decrease, but the torque decrease rate increases.As shown in Figure 4d, the pole shoot height is selected as 1.5 mm.
The selection of the pole arc coefficient and pole shoot width coefficient is more complicated and needs to be determined according to the rated torque, five times overload torque, torque ripple coefficient, and torque drop rate with the change in pole arc coefficient and pole shoot width coefficient, as shown in Figure 5.
In order to obtain a larger torque, only the cases where the polar arc coefficient is greater than 0.81 are considered.In order to obtain a small torque ripple, the selection area is limited to the range enclosed by the black line in Figure 5a.Considering the influence of the pole arc coefficient and pole shoot width coefficient on the rated torque, five times overload torque, torque ripple coefficient, and torque drop rate, three candidate points are selected and listed in Table 8.Through comparison, when the pole ARC coefficient is 0.83 and the pole shoe width coefficient is 0.24, the torque ripple is from the minimum and the torque meets the requirement.Table 9 shows the optimized machine parameter.torque, torque ripple coefficient, and torque drop rate with the change in pole arc coefficient and pole shoot width coefficient, as shown in Figure 5.In order to obtain a larger torque, only the cases where the polar arc coefficient is greater than 0.81 are considered.In order to obtain a small torque ripple, the selection area is limited to the range enclosed by the black line in Figure 5a.Considering the influence of the pole arc coefficient and pole shoot width coefficient on the rated torque, five times overload torque, torque ripple coefficient, and torque drop rate, three candidate points are selected and listed in Table 8.Through comparison, when the pole ARC coefficient is 0.83 and the pole shoe width coefficient is 0.24, the torque ripple is from the minimum and the torque meets the requirement.Table 9 shows the optimized machine parameter.

Performance Analysis of the Optimal Machine
After optimization by the Taguchi method, the optimal variables of the YASA machine are obtained.Next, the performance of the optimized YASA machine is analyzed by FEA.

No-Load Characteristics
Figure 6 shows the waveform of no-load back EMF and the Fourier decomposition, respectively.A1 represents the fundamental amplitude.An represents the harmonic amplitude.The THD of no-load back EMF is 9.4%, obtained by Formula (9).
The harmonic content of no-load back EMF is related to the pole arc coefficient, pole shoot width, and winding connection mode.By Fourier decomposition of rotor MMF with different arc coefficients, curves of each harmonic with pole arc coefficients can be obtained, as shown in Figure 7.The fundamental wave increases with the increase in the pole arc coefficients, but the variation regulation of other harmonics is inconsistent.As the third harmonic will not appear in the three-phase winding of the star arrangement, the fifth and seventh harmonics should be suppressed to pursue a good performance.

No-Load Characteristics
Figure 6 shows the waveform of no-load back EMF and the Fourier decomposition, respectively.A1 represents the fundamental amplitude.An represents the harmonic amplitude.The THD of no-load back EMF is 9.4%, obtained by Formula ( 9).The harmonic content of no-load back EMF is related to the pole arc coefficient, pole shoot width, and winding connection mode.By Fourier decomposition of rotor MMF with different arc coefficients, curves of each harmonic with pole arc coefficients can be obtained, as shown in Figure 7.The fundamental wave increases with the increase in the pole arc coefficients, but the variation regulation of other harmonics is inconsistent.As the third harmonic will not appear in the three-phase winding of the star arrangement, the fifth and seventh harmonics should be suppressed to pursue a good performance.The waveform of air gap magnetic density is shown in Figure 8a, where the depression is mainly caused by stator slots.And stator slot introduces a large number of highorder harmonics.Due to the magnetic leakage at both sides of the YASA machine, the radial air gap magnetic density is lower on both sides compared to the middle, as shown  The harmonic content of no-load back EMF is related to the pole arc coefficient, pole shoot width, and winding connection mode.By Fourier decomposition of rotor MMF with different arc coefficients, curves of each harmonic with pole arc coefficients can be obtained, as shown in Figure 7.The fundamental wave increases with the increase in the pole arc coefficients, but the variation regulation of other harmonics is inconsistent.As the third harmonic will not appear in the three-phase winding of the star arrangement, the fifth and seventh harmonics should be suppressed to pursue a good performance.The waveform of air gap magnetic density is shown in Figure 8a, where the depression is mainly caused by stator slots.And stator slot introduces a large number of highorder harmonics.Due to the magnetic leakage at both sides of the YASA machine, the radial air gap magnetic density is lower on both sides compared to the middle, as shown The waveform of air gap magnetic density is shown in Figure 8a, where the depression is mainly caused by stator slots.And stator slot introduces a large number of high-order harmonics.Due to the magnetic leakage at both sides of the YASA machine, the radial air gap magnetic density is lower on both sides compared to the middle, as shown in Figure 8b.The no-load magnetic density of the YASA machine is shown in Figure 9, where the maximum value is 1.8 T. in Figure 8b.The no-load magnetic density of the YASA machine is shown in Figure 9, where the maximum value is 1.8 T.
(a) (b)  In the application of in-wheel traction, the YASA machine needs to run smoothly and In the application of in-wheel traction, the YASA machine needs to run smoothly and have a small torque ripple.Cogging torque is one of the main factors causing torque ripple, which puts forward strict requirements.As shown in Figure 10, the peak-to-peak value of cogging torque, whose cycle is 20 electrical angles, is 1 Nm.
The axial magnetic tension generated by the rotors on both sides of the stator core offsets each other, making that on the stator core nearly zero, as shown in Figure 11.As the rotor position changes, the axial magnetic tension on the stator varies periodically with a maximum value of 1.5 Nm.The axial magnetic tension on the rotor changes little with the current density, as shown in Figure 12.At different current densities, the maximum and minimum values of magnetic tension of the rotor core are 6250 Nm and 6200 Nm, respectively.The reason for this is that the air gap magnetic field is mainly provided by permanent magnets, and the armature magnetic field has little effect on the air gap magnetic field.In the application of in-wheel traction, the YASA machine needs to run smoothly and have a small torque ripple.Cogging torque is one of the main factors causing torque ripple, which puts forward strict requirements.As shown in Figure 10, the peak-to-peak value of cogging torque, whose cycle is 20 electrical angles, is 1 Nm.
The axial magnetic tension generated by the rotors on both sides of the stator core offsets each other, making that on the stator core nearly zero, as shown in Figure 11.As the rotor position changes, the axial magnetic tension on the stator varies periodically with a maximum value of 1.5 Nm.The axial magnetic tension on the rotor changes little with the current density, as shown in Figure 12.At different current densities, the maximum and minimum values of magnetic tension of the rotor core are 6250 Nm and 6200 Nm, respectively.The reason for this is that the air gap magnetic field is mainly provided by permanent magnets, and the armature magnetic field has little effect on the air gap magnetic field.In the application of in-wheel traction, the YASA machine needs to run smoothly and have a small torque ripple.Cogging torque is one of the main factors causing torque ripple, which puts forward strict requirements.As shown in Figure 10, the peak-to-peak value of cogging torque, whose cycle is 20 electrical angles, is 1 Nm.
The axial magnetic tension generated by the rotors on both sides of the stator core offsets each other, making that on the stator core nearly zero, as shown in Figure 11.As the rotor position changes, the axial magnetic tension on the stator varies periodically with a maximum value of 1.5 Nm.The axial magnetic tension on the rotor changes little with the current density, as shown in Figure 12.At different current densities, the maximum and minimum values of magnetic tension of the rotor core are 6250 Nm and 6200 Nm, respectively.The reason for this is that the air gap magnetic field is mainly provided by permanent magnets, and the armature magnetic field has little effect on the air gap magnetic field.

Load Characteristics
When the current density of the YASA machine is 6 A/mm 2 , the corresponding armature current is 17 A (RMS).As shown in Figure 13, the average value of the output torque, whose torque ripple is 2.31%, is 99 Nm.The maximum value of the cogging torque is

Load Characteristics
When the current density of the YASA machine is 6 A/mm 2 , the corresponding armature current is 17 A (RMS).As shown in Figure 13, the average value of the output torque, whose torque ripple is 2.31%, is 99 Nm.The maximum value of the cogging torque is

Load Characteristics
When the current density of the YASA machine is 6 A/mm 2 , the corresponding armature current is 17 A (RMS).As shown in Figure 13, the average value of the output torque, whose torque ripple is 2.31%, is 99 Nm.The maximum value of the cogging torque is around 0.51% of the rated torque.The magnetic density of the YASA machine under rated load is shown in Figure 14.The maximum value of the magnetic density is 1.90 T, and the core saturation is located in the tooth boot of the stator core.

Load Characteristics
When the current density of the YASA machine is 6 A/mm 2 , the corresponding armature current is 17 A (RMS).As shown in Figure 13, the average value of the output torque, whose torque ripple is 2.31%, is 99 Nm.The maximum value of the cogging torque is around 0.51% of the rated torque.The magnetic density of the YASA machine under rated load is shown in Figure 14.The maximum value of the magnetic density is 1.90 T, and the core saturation is located in the tooth boot of the stator core.

Load Characteristics
When the current density of the YASA machine is 6 A/mm 2 , the corresponding armature current is 17 A (RMS).As shown in Figure 13, the average value of the output torque, whose torque ripple is 2.31%, is 99 Nm.The maximum value of the cogging torque is around 0.51% of the rated torque.The magnetic density of the YASA machine under rated load is shown in Figure 14.The maximum value of the magnetic density is 1.90 T, and the core saturation is located in the tooth boot of the stator core.The torque drop rate is an important optimization objective of the YASA machine.As shown in Figure 15, when the current density is less than 30 A/mm 2 , the overload curve of the YASA machine has good linearity.And the torque drop rate is only 5.99%, which is due to reasonable optimization variable selection.The torque drop rate is an important optimization objective of the YASA machine.As shown in Figure 15, when the current density is less than 30 A/mm 2 , the overload curve of the YASA machine has good linearity.And the torque drop rate is only 5.99%, which is due to reasonable optimization variable selection.
In this paper, the YASA machine is a surface-mounted rotor structure, the base speed (480 r/min) is controlled by Id = 0, and all currents are used to generate torque.Above the base speed, the field weakening control is used, and part of the current is used for field weakening.Figure 16 shows the external characteristic curve and corresponding efficiency.Before the base speed, the torque is basically maintained at 99 Nm.After the base speed, the torque decreases gradually with the increase in speed.The efficiency under the rated operating condition is 94.86%, and the maximum efficiency of the YASA machine is 95.12%.In this paper, the YASA machine is a surface-mounted rotor structure, the base speed (480 r/min) is controlled by I d = 0, and all currents are used to generate torque.Above the base speed, the field weakening control is used, and part of the current is used for field weakening.Figure 16 shows the external characteristic curve and corresponding efficiency.Before the base speed, the torque is basically maintained at 99 Nm.After the base speed, the torque decreases gradually with the increase in speed.The efficiency under the rated operating condition is 94.86%, and the maximum efficiency of the YASA machine is 95.12%.Figure 17 shows the variation of loss components with regard to the machine speed.The core loss and eddy current loss show an obvious upward trend, with the speed increasing, but there is a turning point at base speed.The eddy current loss and core loss are not only related to the speed but also magnetic density.Below the base speed, the loss increases with the increase in speed.Above the base speed, magnetic density decreases under field weakening control, so the loss at the point is slightly lower.The stator core loss, whose loss distribution is shown in Figure 18, increases with the speed.For the eddy current loss of the rotor core and permanent magnets, the weakening effect of the weak field on magnetic density counteracts the increasing effect of speed increase by more than half, so the eddy current loss increases smoothly above the base speed.Figure 17 shows the variation of loss components with regard to the machine speed.The core loss and eddy current loss show an obvious upward trend, with the speed increasing, but there is a turning point at base speed.The eddy current loss and core loss are not only related to the speed but also magnetic density.Below the base speed, the loss increases with the increase in speed.Above the base speed, magnetic density decreases under field weakening control, so the loss at the point is slightly lower.The stator core loss, whose loss distribution is shown in Figure 18, increases with the speed.For the eddy current loss of the rotor core and permanent magnets, the weakening effect of the weak field on magnetic density counteracts the increasing effect of speed increase by more than half, so the eddy current loss increases smoothly above the base speed.

Experimental Validation
According to the optimization results by the Taguchi method, the YASA machine is manufactured for experimental validation.The prototype and experimental platform are shown in Figure 19.

Experimental Validation
According to the optimization results by the Taguchi method, the YASA machine is manufactured for experimental validation.The prototype and experimental platform are shown in Figure 19.

Experimental Validation
According to the optimization results by the Taguchi method, the YASA machine is manufactured for experimental validation.The prototype and experimental platform are shown in Figure 19.
The above experimental platform was used to test the back EMF of the prototype, and the test results are shown in Figure 20, showing good consistency with the finite element results.The error between the experimental test and the FEA mainly comes from the difference between the actual permanent magnetic materials and the materials used in the FEA simulation.

Experimental Validation
According to the optimization results by the Taguchi method, the YASA machine is manufactured for experimental validation.The prototype and experimental platform are shown in Figure 19.The above experimental platform was used to test the back EMF of the prototype, and the test results are shown in Figure 20, showing good consistency with the finite element results.The error between the experimental test and the FEA mainly comes from the difference between the actual permanent magnetic materials and the materials used in the FEA simulation.The load test was carried out, and the torque current curve of the prototype was measured as shown in Figure 21.It can be seen from Figure 21, the results of the FEA are in good agreement with the test results.And the error between the experimental test and the FEA, which is about 5%, is relatively constant.The error is within a reasonable range.The measured back EMF and torque verify the accuracy of the model used in this study.

Conclusions
In this paper, a 36s/32p YASA machine applied on an in-wheel traction application is optimized by the improved Taguchi method.According to the design requirements, the average rated torque, torque ripple coefficient, and overload capacity of the YASA machine are set as the optimization objectives.And the optimization variables (the rotor core thickness, pole arc coefficient, etc.) are determined through relationship analysis of the structural parameters.Based on the orthogonal experimental matrix and sensitivity analysis of the FEA results, the high sensitivity levels of the optimization variables are determined.The optimal value of each optimization variable is obtained through further detailed optimization of the optimization variables with high sensitivity levels.The no-load and on-load performances of the YASA machine are calculated by FEA.Finally, the accuracy of the optimization process is experimentally verified by the prototyped machine with optimal variables.For future study, the combination of thermal network models in [20] and the proposed Taguchi method here would be a good direction, and other useful analytical models for machine design can also be included.The load test was carried out, and the torque current curve of the prototype was measured as shown in Figure 21.It can be seen from Figure 21, the results of the FEA are in good agreement with the test results.And the error between the experimental test and the FEA, which is about 5%, is relatively constant.The error is within a reasonable range.The measured back EMF and torque verify the accuracy of the model used in this study.The load test was carried out, and the torque current curve of the prototype was measured as shown in Figure 21.It can be seen from Figure 21, the results of the FEA are in good agreement with the test results.And the error between the experimental test and the FEA, which is about 5%, is relatively constant.The error is within a reasonable range.The measured back EMF and torque verify the accuracy of the model used in this study.

Conclusions
In this paper, a 36s/32p YASA machine applied on an in-wheel traction application is optimized by the improved Taguchi method.According to the design requirements, the average rated torque, torque ripple coefficient, and overload capacity of the YASA machine are set as the optimization objectives.And the optimization variables (the rotor core thickness, pole arc coefficient, etc.) are determined through relationship analysis of the structural parameters.Based on the orthogonal experimental matrix and sensitivity analysis of the FEA results, the high sensitivity levels of the optimization variables are determined.The optimal value of each optimization variable is obtained through further de-

Conclusions
In this paper, a 36s/32p YASA machine applied on an in-wheel traction application is optimized by the improved Taguchi method.According to the design requirements, the average rated torque, torque ripple coefficient, and overload capacity of the YASA machine are set as the optimization objectives.And the optimization variables (the rotor core thickness, pole arc coefficient, etc.) are determined through relationship analysis of the structural parameters.Based on the orthogonal experimental matrix and sensitivity analysis of the FEA results, the high sensitivity levels of the optimization variables are determined.The optimal value of each optimization variable is obtained through further detailed optimization of the optimization variables with high sensitivity levels.The no-load and on-load performances of the YASA machine are calculated by FEA.Finally, the accuracy of the optimization process is experimentally verified by the prototyped machine with optimal variables.For future study, the combination of thermal network models in [20] and the proposed Taguchi method here would be a good direction, and other useful analytical models for machine design can also be included.

Figure 1 .
Figure 1.Topology of the YASA machine.

Figure 1 .
Figure 1.Topology of the YASA machine.

Figure 2 .
Figure 2. Structural parameters of the YASA machine.

Figure 2 .
Figure 2. Structural parameters of the YASA machine.

Figure 3 .
Figure 3.The flowchart of the optimization process.

Figure 3 .
Figure 3.The flowchart of the optimization process.

Figure 4 .
Figure 4. Variation curve of optimization objective with variables.(a) Rotor core thickness hr, (b) permanent magnet thickness hm, (c) slot width ws, (d) pole shoot height hv.

Figure 4 .
Figure 4. Variation curve of optimization objective with variables.(a) Rotor core thickness h r , (b) permanent magnet thickness h m , (c) slot width w s , (d) pole shoot height h v . α

Figure 5 .
Figure 5.The map of the optimization objective with the pole arc coefficient α and pole shoot width coefficient kp.(a) Rated torque Tave, (b) torque ripple krip, (c) five times overload torque T5ave, (d) torque drop rate ka.

Figure 5 .
Figure 5.The map of the optimization objective with the pole arc coefficient α and pole shoot width coefficient k p .(a) Rated torque T ave , (b) torque ripple k rip , (c) five times overload torque T 5ave , (d) torque drop rate k a .

Figure 7 .
Figure 7. Harmonic content of rotor MMF under different pole arc coefficients.

Figure 7 .
Figure 7. Harmonic content of rotor MMF under different pole arc coefficients.

Figure 7 .
Figure 7. Harmonic content of rotor MMF under different pole arc coefficients.

Figure 8 .
Figure 8. No-load air gap flux density.(a) Air gap magnetic density, (b) radial air gap magnetic density.

Figure 8 .
Figure 8. No-load air gap flux density.(a) Air gap magnetic density, (b) radial air gap magnetic density.

Figure 8 .
Figure 8. No-load air gap flux density.(a) Air gap magnetic density, (b) radial air gap magnetic density.

Figure 12 .
Figure 12.Variation curve of rotor core axial magnetic pulling force with current density.

Figure 12 .
Figure 12.Variation curve of rotor core axial magnetic pulling force with current density.

Figure 12 .
Figure 12.Variation curve of rotor core axial magnetic pulling force with current density.

Figure 12 .
Figure 12.Variation curve of rotor core axial magnetic pulling force with current density.

Figure 12 .
Figure 12.Variation curve of rotor core axial magnetic pulling force with current density.

Figure 14 .
Figure 14.Magnetic density distribution under rated load.

Figure 15 .
Figure 15.Variation curve of machine torque with current density.

Figure 15 .
Figure 15.Variation curve of machine torque with current density.

Figure 15 .
Figure 15.Variation curve of machine torque with current density.

Figure 18 .
Figure 18.The loss distribution of stator core.

Figure 18 .
Figure 18.The loss distribution of stator core.

Figure 18 .
Figure 18.The loss distribution of stator core.

Figure 18 .
Figure 18.The loss distribution of stator core.

Figure 20 .
Figure 20.Test results of no-load back EMF.

Figure 21 .
Figure 21.Torque current curve of the prototype.

Figure 20 .
Figure 20.Test results of no-load back EMF.

Figure 20 .
Figure 20.Test results of no-load back EMF.

Figure 21 .
Figure 21.Torque current curve of the prototype.

Figure 21 .
Figure 21.Torque current curve of the prototype.

Table 1 .
Main parameters of the YASA machine.

Table 2 .
Initial structural parameters of the YASA machine.

Table 1 .
Main parameters of the YASA machine.

Table 2 .
Initial structural parameters of the YASA machine.

Table 3 .
Optimization variables and levels.

Table 4 .
Experimental matrix and FEA results.

Table 5 .
Average value of each optimization objective.

Table 6 .
Average value of objective under different levels.

Table 7 .
Proportion of each variable on the optimization objective.

Table 8 .
Different values of pole ARC coefficient and pole shoot width coefficient and optimization objective.

Table 8 .
Different values of pole ARC coefficient and pole shoot width coefficient and optimization objective.

Table 9 .
Structural parameters of YASA machine after optimization.

Table 9 .
Structural parameters of YASA machine after optimization.