Optimization Analysis of Automotive Asymmetric Magnetic Pole Permanent Magnet Motor by Taguchi Method
In order to improve the air-gap flux density of the permanent magnet synchronous motor and reduce the cogging torque, a novel structure with asymmetric magnetic poles for automobile was proposed. Based on the characteristics of the parallel magnetic circuit, the magnetic flux path diagram is established. And the equivalent magnetic circuit model is established by the equivalent magnetic circuit method. The Taguchi method is used to be a multiobjective optimization algorithm. The total harmonic distortion of the air-gap flux density is the first optimization goal. The second and third optimization goals are the cogging torque and the average of output torque, respectively. And the torque ripple is a constraint condition. The optimized parameter combination is obtained by the Taguchi method. Finite element simulation analysis and prototype test are carried out for the optimized motor structure. The results show that the total harmonic distortion of air-gap flux density is reduced by 36.7% comparing with the initial structure. The cogging torque is reduced by 26.0%. And the average output torque is increased by 4.8%.
Electric automobile has been widely concerned around the world because of the superiority of energy saving and environmental protection. As one of the core technologies of electric automobiles, the driving motor system has become the key component. High power density, high efficiency, and wide speed range have become current research directions. At present, driving motors of electric automotive include induction motor, direct current motor, and permanent magnet synchronous motor (PMSM). PMSM has become the most promising drive motor due to its advantages of high power density, high efficiency, and simple structure [1, 2].
According to the relative position of the permanent magnet and the rotor core, PMSM can be divided into surface mount type and interior type. The interior type can be divided into three types: radial type, tangential type, and hybrid type. Their magnetization directions of the permanent magnets are different. At present, the interior PMSM is used widely in industrial engineering. Many automotive companies such as Daimler, Toyota, and BYD have widely applied interior PMSM to their electric automobile [3, 4].
However, the air-gap flux density waveform of PMSM is similar to a square wave. It has many higher harmonics. And high cogging torque is not conducive to the output performance .
In view of above disadvantages, domestic and foreign scholars carried out studies on the influence of output performance and different magnetic pole shapes. It is determined that the air-gap flux density and the cogging torque are greatly affected by the shape of the magnetic poles. Appropriate magnetic pole shape can improve the air-gap flux density and reduce cogging torque [6, 7]. Toyota adopted an interior PMSM structure of V-type and tangential type on Prius hybrid automobile. Oak Ridge National Laboratory proposed a U-type motor structure. The air-gap flux density is improved, and the cogging torque is reduced. But the amount of permanent magnets is increased . Biasing the magnetic poles and changing the pole arc coefficient can also improve the motor performance [9–11]. In , the outer edge of the magnet of the surface mount PMSM was cut to improve air-gap magnetic field. But it makes the manufactural process more difficult. And it can cause sudden changes in magnetic field. For optimization methods, domestic scholars proposed the Taguchi method for multiobjective optimization [13, 14]. It can reduce the number of experiments and improve efficiency. But the traditional Taguchi method has the same number of level values. It is difficult to obtain the best optimized solution. And the optimal motor structure cannot be obtained.
In order to improve the air-gap flux density and the cogging torque, this paper proposes a PMSM structure with asymmetric magnetic poles. The equivalent magnetic circuit model is established. And the effective magnetic flux is analyzed by the equivalent magnetic circuit method. The rotor structure is improved. Optimization goals are selected. Some key parameters are selected to be optimization variables. Using the Taguchi method, a multiobjective optimization scheme is proposed. And a two-dimensional finite element method and prototype experiment are used to verify.
2. Motor Structure and Magnetic Circuit Analysis
2.1. Motor Structure
This paper proposes PMSM structure asymmetric magnetic poles. Its three-dimensional structure is shown in Figure 1. It is an interior-type PMSM with 8 poles and 12 slots. Its winding is concentrated winding. Each of its magnetic poles is a single-layer arrangement. The shape of magnetic poles is similar with “.” Each magnetic pole is combined with two pieces of rare earth permanent magnet materials. The thickness and length of them are different. The relevant structure parameters are shown in Table 1.
2.2. Motor Effective Flux Analysis
This paper takes a traditional V-type structure as comparing structure. Each pole effective magnetic flux of two motor structures is analyzed. Figure 2 shows the comparison diagrams of the flux paths. It can be seen that the effective flux of the traditional V-type structure is a single path. It is provided by two adjacent permanent magnets. The effective flux of asymmetric magnetic pole structure consists of two paths in parallel. A path is provided by a single long permanent magnet. And another path is provided by two adjacent permanent magnets.
Schematic diagrams of flux paths have been shown above. The equivalent magnetic circuit models of two motor structures are established, respectively, by the equivalent magnetic circuit method. And they are shown in Figure 3.
In Figure 3, and are the magnetomotive force. They are, respectively, provided by the long permanent magnet and the short permanent magnet. is the direct axis component of the armature reaction. and are the internal equivalent permeance of the long permanent magnet and the short permanent magnet, respectively. is the permeance between the long permanent magnet and the rotor core. is the permeance between the short permanent magnet and the rotor core. is the leakage permeance of long permanent magnet. And is the leakage permeance of short permanent magnet. is the air-gap permeance. is the permeance of stator teeth. is the permeance of the stator yoke. is the permeance of the rotor core between long permanent magnet and the air gap. And is the performance of the rotor core between short permanent magnet and the air gap. is the permeance of the rotor core between the long and short permanent magnets. is the total magnetic flux in a magnetic circuit. It is provided by the long and short permanent magnet. is the effective magnetic flux. and are the leakage flux of long permanent magnet and short permanent magnet, respectively.
In Figure 3(b), is an extra magnetic flux in the equivalent magnetic circuit model. And it only appears in an asymmetric magnetic pole structure. and are connected in parallel. is used to represent the total permeance of , , , , and .
According to the superposition principle, meets the following conditions:
is used to represent the total permeance of and . It meets the following conditions:Above all, the structure with asymmetric magnetic poles can reduce the total permeance in the magnetic circuit. Its total effective magnetic flux is more than a traditional V-type structure if the magnetomotive force is the same.
3. Motor Structure Optimization
In order to obtain a better geometric motor shape, this paper makes some improvements on the rotor structure. The shape of rotor edge is changed. The auxiliary slots are used inside the rotor. The length of the air gap becomes nonuniform. The parameters are optimized by an appropriate optimization method. And the performance of the rated operating point is selected for research and analysis.
3.1. Optimization Ideas and Optimization Methods
The method of changing the air-gap length is used in this paper. The air-gap length becomes nonuniform. And the effective magnetic flux direction is adjusted. The specific measures are as follows. The auxiliary slots are used inside the rotor. It can enhance the magnetization ability and reduce the rotational inertia. The distance between auxiliary slots and a permanent magnet is . The outer edge of the rotor is changed to be eccentric. The eccentric distance is . Straight line is used in the junction of the eccentric circle. Its length is . It can eliminate sharp corners on the outer edge, as shown in Figure 4.
The above improvement method can improve the air-gap flux density waveform. The schematic diagram of the expected waveform can be represented in Figure 5. It is more sinusoidal than the initial structure.
This paper uses the Taguchi method to optimize the parameters. The Taguchi method was first proposed by Dr. Genichi Taguchi who is a Japanese quality control expert. It is a local optimization method based on orthogonal experiments. It can quickly determine the combination of multiobjective optimization with the least number of experiments. The basic design process is as follows:(1)Determine optimization goals according to requirements(2)Determine optimization variables and corresponding level values(3)Establish an orthogonal experimental matrix(4)Use the finite element method to solve the orthogonal experiment matrix(5)Analyze the influence of each optimization variable and determine the optimal parameters(6)Verify the optimization results by finite element analysis
3.2. Determination of Optimization Goals and Optimization Variables
In this paper, the Taguchi method is used to optimize the air-gap flux density waveform, the cogging torque, and the output torque. And it must ensure that the torque ripple does not increase. Therefore, the total harmonic distortion of the air-gap flux density (THD) is the first optimization goal. The cogging torque () is the second optimization goal. The average of output torque () is the third optimization goal. And the torque ripple is a constraint condition. It requires that THD and reduces and increases as much as possible. And the torque ripple does not increase.
In addition, is chosen as the optimization variable . is chosen as the optimization variable . And is chosen as the optimization variable . The initial values of the three optimization variables are (2, 6.5, 4). On this basis, appropriate values are selected as the level values. And optimization parameters and level values are shown in Table 2. Especially, “none” means no auxiliary slot.
3.3. Orthogonal Experimental Designs and Solution
The optimization goals, optimization variables, and corresponding level values have been determined above. The number of level values is not exactly the same. The traditional standard orthogonal experiment matrix cannot be established. According to principles of the Taguchi method experiment designing, the nonstandard orthogonal experimental matrix is established as shown in Table 3.
The rated operating point performance of each experiment can be calculated by using the two-dimensional finite element analysis method. The orthogonal experimental results are shown in Table 4. The calculation formula of THD is shown in formula (3). And the calculation of is shown in formula (5). can be obtained by the simulation analysis.where is the effective value of the fundamental wave, is the effective value of higher harmonics, and and are the maximum and minimum output torque, respectively. And they can be obtained by the simulation analysis.
3.4. Results Processing and Analysis
In the previous section, the results of each experiment were obtained by the two-dimensional finite element method. The overall average value of the experimental results can be calculated by formula (5). And the results are shown in Table 5.
The average value of each optimization variable can be calculated, respectively, at a certain level value. To give an example, the average value can be obtained by the following formula when the optimization variable is selected as the level 1.Similarly, the average values of THD, , and can be obtained with reference to the calculation method of the above formula. And the results of THD, , and are shown in Table 6, Table 7, and Table 8, respectively.
If a single goal is optimized, the optimization combination can be selected in Table 6, Table 7, and Table 8. To optimize THD, a combination of A4, B1, and C2 can be selected. To optimize , a combination of A4, B2, and C3 can be selected. And a combination of A4, B3, and C3 can be selected to optimize .
It is necessary to analyze the influence of each optimization variable when multiobjective optimization is carried out. Based on the previous calculation, the proportions can be calculated by formula (7).where is the value of each optimization goal under a certain level; is the average value of each optimization goal; is the number of levels, and SS is the variance of the optimization variable under a certain performance. And the calculation results are shown in Table 9.
In the above, the experimental results of the Taguchi method optimization are analyzed. And the proportions of the influence of each variable on the optimization goal are obtained. It can be seen that optimization variable has a great influence on the three optimization goals. THD and are the smallest when level 4 is selected. And is the largest. Optimization variable has a great influence on . Therefore, optimization variable chooses level 2. And optimization variable has a little influence on each optimization goal and chooses level 3. In other words, the eccentric distance is 8 mm. The length of the straight line is 6.5 mm. And the distance between auxiliary slots and a permanent magnet is 4.5 mm.
4. Simulation Analysis and Experimental Verification
4.1. Finite Element Simulation Analyses
According to the above optimization scheme, the motor structure after optimization is obtained and designed. The structure is shown in Figure 6. And the static magnetic field analysis is shown in Figure 7.
The optimized motor structure is analyzed by finite element simulation. The air-gap flux density waveform and harmonic analysis are compared and analyzed in Figure 8. The waveform after optimization meets an excepted waveform.
Cogging torque comparison is shown in Figure 9. And output torque comparison is shown in Figure 10. With the two-dimensional finite element method simulation analysis, it can be calculated that THD is 28.9%. The is 284.3 mN∙m. And the is 13.7 N∙m.
It can be seen from Figure 10 that the is improved. And the range of the output torque ripple is reduced. The torque ripple rate () represents the torque ripple of the output torque. Its calculation is as follows:It can be calculated that before optimization is 41.1% by the above formula. And after optimization is 24.7%. It decreases by 16.4%.
The performance comparison before and after optimization is shown in Figure 11. Compared with before optimization, THD is reduced by 36.7%. is reduced by 26.0%. And is increased by 4.8%.
4.2. Prototype Manufacture and Experiment
In order to verify the effectiveness of the simulation analysis, the rotor structure after optimization is manufactured. It is shown in Figure 12. The prototype is assembled and placed on a cogging torque test platform for experiment. The cogging torque test platform is shown in Figure 13.
The cogging torque experimental curve is shown in Figure 14. The maximum value is 280.8 mN∙m. It is similar to the simulation value. Maximum friction torque is 346.8 mN∙m.
In order to obtain the maximum torque and efficiency of asymmetrical magnetic pole PMSM, it is necessary to test the characteristic curve of the motor. The prototype is placed on the dynamometer test platform for a full load test. It is shown in Figure 15. And the characteristic curve of the motor is obtained as shown in Figure 16.
In Figure 16, the maximum speed of the motor is 5000 revolutions per minute. The maximum output torque is 62 N∙m. And maximum efficiency exceeds 90%.
This paper proposes a novel asymmetric magnetic pole PMSM for automobiles. The equivalent magnetic circuit model is established. Many improvements have been made to the rotor structure. The Taguchi method is used to obtain the optimal solution of each optimization variable. Combined with finite element analysis and prototype experiment, the effectiveness of the structure and optimization method is verified. In addition, some conclusions are drawn as follows:(1)The novel asymmetric magnetic pole permanent magnet synchronous motor has parallel effective magnetic flux. Compared with the traditional V-type structure, the total permeance is reduced. And the effective magnetic flux content is increased(2)The Taguchi method is used to optimize the motor structure. A nonstandard orthogonal experimental matrix is established. The original 54 sets of experiments are replaced by 18 combinations of experiments. The number of experiments is reduced obviously. And the finite element method verifies the effectiveness of the method(3)The eccentric distance of the rotor has the greatest influence on the total harmonic distortion of the air-gap flux density and the average output torque. The eccentric distance and the length of the straight line both have influence on the cogging torque. The proportions are approximately 50%(4)Compared with the performance before optimization, the total harmonic distortion of the air-gap flux density reduced by 36.7%. The cogging torque is reduced by 26.0%. The average output torque is increased by 4.8%. And the output torque ripple is reduced by 16.4%
All data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
This research work is partially supported by the National Natural Science Foundation of China (Grant Nos. 51975340 and 51875327).
W. Cao, B. C. Mecrow, G. J. Atkinson, J. W. Bennett, and D. J. Atkinson, “Overview of electric motor technologies used for more electric aircraft (MEA),” IEEE Transactions on Industrial Electronics, vol. 59, no. 9, pp. 3523–3531, 2011.View at: Publisher Site | Google Scholar
Y. Guo, L. Shi, W. Zhang, C. Wenping, and K. Byung II, “Research on 3/6-phase winding switching integrated starter generator for range extender electric vehicles,” International Conference on Power Electronics and ECCE Asia, pp. 2492–2497, 2019.View at: Google Scholar
R. Xu and C. Qin, “Development status and prospect of electric vehicle and its driving motor,” Southern Power System Technology, vol. 10, no. 71, pp. 82–86, 2016.View at: Google Scholar
S. Ye, “Design and performance analysis of an iterative flux sliding-mode observer for the sensorless control of PMSM drives,” ISA Transactions, vol. 94, pp. 255–264, 2019.View at: Publisher Site | Google Scholar
K. Liu and J. Xu, “Life cycle evaluation of pure electric vehicle driving moto,” Acta Scientiae Circumstantiae., vol. 36, no. 9, pp. 3456–3463, 2016.View at: Google Scholar
X. Liu, H. Chen, J. Zhao, and A. Belahcen, “Research on the performances and parameters of interior PMSM used for electric vehicles,” IEEE Transactions on Industrial Electronics, vol. 63, no. 6, pp. 3533–3545, 2016.View at: Publisher Site | Google Scholar
W. Hu, X. Zhang, Y. Lei, Q. du, L. Shi, and G. Liu, “Analytical model of air-gap field in hybrid excitation and interior permanent magnet machine for electric logistics vehicles,” IEEE Access, vol. 8, pp. 148237–148249, 2020.View at: Publisher Site | Google Scholar
Q. Li, C. Sun, T. Fan, X. Wen, and Y. Li, “Rotor pole optimal design of double-layer interior permanent magnet machines for electric vehicle,” Electric Machines and Control, vol. 24, no. 9, pp. 56–64, 2020.View at: Google Scholar
W. Ren, “A technique of torque ripple compensation for interior permanent magnet machines by using an asymmetrical V-type rotor configuration,” in Annual Conference of the Industrial Electronics Society, pp. 8744–8749, Beijing, China, 2017.View at: Publisher Site | Google Scholar
Q. Chen, G. Xu, G. Liu, W. Zhao, L. Liu, and Z. Lin, “Torque ripple reduction in five-phase IPM motors by lowering interactional MMF,” IEEE Transactions on Industrial Electronics, vol. 65, no. 11, pp. 8520–8531, 2018.View at: Publisher Site | Google Scholar
W. Ren, Q. Xu, Q. Li, and L. Zhou, “Reduction of cogging torque and torque ripple in interior PM machines with asymmetrical V-type rotor design,” IEEE Transactions on Magnetics, vol. 52, no. 7, pp. 1–5, 2016.View at: Publisher Site | Google Scholar
F. Chai, P. Liang, Y. Pei, and S. Cheng, “Magnet shape optimization of surface-mounted permanent-magnet motors to reduce harmonic iron losses,” IEEE Transactions on Magnetics, vol. 52, no. 7, pp. 1–4, 2016.View at: Publisher Site | Google Scholar
J. He, G. LI, R. Zhou, and Q. Wang, “Optimization of permanent-magnet spherical motor based on Taguchi method,” IEEE Transactions on Magnetics, vol. 52, no. 2, pp. 1–7, 2020.View at: Google Scholar
C. Xia, L. Guo, Z. Zhang, T. Shi, and H. Wang, “Optimal designing of permanent magnet cavity to reduce iron loss of interior permanent magnet machine,” IEEE Transactions on Magnetics, vol. 51, no. 12, pp. 1–9, 2015.View at: Google Scholar