Optimization Design of a New Type of Interior Permanent Magnet Generator for Electric Vehicle Range Extender

Aiming at the disadvantages of large leakage flux and low magnetic flux density of radial magnetic circuit and tangential magnetic circuit, a new type of permanent magnet (PM) rotor with parallel tangential and radial magnetic circuits is proposed. Based on Ohm’s law and Kirchhoff’s law of magnetic circuits, equivalent magnetic circuits for rotor poles are developed. *e structure parameters of the generator are preliminarily determined. At the same time, by means of the Taguchi method and employing finite element analysis, the rotor poles of generator are optimized to improve air gap magnetic density, the cogging torque, and the distortion of back-EMF waveform. Finally, the validity of proposed design methods is validated by the analytical and experimental results.


Introduction
Electric vehicles have been widely used in people's daily life due to their advantages of low pollution, low noise, and zero emission.But, the limited range has become a bottleneck restricting its development.Before electric vehicle battery storage technology made a major breakthrough, adding an extended-range generator for electric vehicles was one of the important ways to increase the endurance mileage [1,2].So, the generator device is the key component of the power supply system of electric vehicles.With the perfection of the electric automobile, the power consumption is increasing and the silicon-rectifying generator has been unable to meet the demand of electrical equipment for power consumption.However, the permanent magnet generator is excited by PM, without electric excitation winding, and has the advantages of simple structure, high power density, high reliability, and so on, so it has a broad market prospect in the power supply system of electric vehicles [3,4].
e literature [5] has proposed a novel surface-mounted and interior double radial PM synchronous generator, through the equivalent magnetic circuit (EMC) method.e magnetic permeance and magnetic leakage permeance of the generator are analyzed, and the optimal structural parameters of the generator are determined.One study [6] showed that using the level set method to optimize the PM synchronous generator, the cogging torque of the generator is reduced and the performance of the motor is improved.e literature [7] has proposed an axial magnetic circuit generator in which rotor pole is composed of core and PM and by optimizing the shape of the magnetic pole, harmonic content is reduced.e literature [8] presented a new type of rotor structure, in which main magnetic flux and leakage flux of the PM are calculated by the EMC method, the iteration does not converge when iron core is severely saturated.e literature [9] has established the open-circuit EMC model of the surface-mounted and interior synchronous PM motor and the EMC model under the armature reaction.e air gap magnetic field and the armature back EMF of the motor are solved; the accuracy of the analytical calculation is verified by the finite element results.e literature [10] makes use of genetic algorithm and particle swarm optimization to improve the comprehensive performance of the motor, but establishment and solution analysis of objective function in these algorithms are quite complex, which not only makes it difficult to realize the fast calculation of the optimal parameters but also has certain localization in the multiobjective optimization design.However, the Taguchi method is a local optimization algorithm which can realize the multiobjective optimization design of a generator.It not only can realize the rapid design of the generator but also has high design precision.It has been widely used in the field of generator development and design in recent years [11,12].
In order to accurately design and analyze the relationship between the output characteristics of the tangential and radial parallel magnetic circuit permanent magnet generator (ITQPMG), the parameters of magnetic pole are determined.By establishing the model of EMC, the parameters of magnetic pole are preliminarily determined, and then, they are optimized by the Taguchi method.In this study, the parameters of magnetic pole have been selected as the horizontal variables.
e multiobjective optimization of output characteristics such as the peak value of air gap magnetic density, waveform distortion rate of no-load back EMF, and the peak value of cogging torque is realized, so the optimal parameters combination can be obtained to improve the performance of ITQPMG.

Determination of Main Parameters
e magnetic circuit of the conventional PM generator is divided into radial directions and tangential directions.Although the generator of radial magnetic circuit has less leakage and back-EMF waveform has good sinusoidal property, the peak value of air gap flux density is lower.e interior radial permanent magnet generator (IRPMG) is shown in Figure 1(a); when the leakage of tangential magnetic circuit generator is large, tangential magnetic poles can produce a certain degree of magnet congregate effect.So the peak value of air gap density is high and the distortion rate of back-EMF waveform is low [13].e interior tangential permanent magnet generator (ITPMG) is shown in Figure 1(b).In this study, the advantages of two kinds of PM generator of magnetic circuit are discussed and a new type of PM generator is proposed, which reduces the magnetism leakage between magnetic poles and tends back-EMF waveforms to be sinusoidal.e radial and tangential flux together provides the flux and synthesizes it in air gap, which has remarkable magnetic concentration effect, makes up the depression of back-EMF waveforms, and improves the efficiency of the generator.e structure of ITQPMG is shown in Figure 1(c).Initial design parameters are listed in Table 1.

Establishment and Analysis of EMC Model
According to the PM layout form and the characteristics of rotor topological structure, the EMC model is established under the no-load condition of the generator.e model consists of two independent flux paths: the first magnetic flux path is a closed magnetic path independently formed by tangential PM and the second flux paths are closed magnetic circuits formed by two series radial rectangle PM connections.e main flux path and leakage flux path of the new type structure of PM generator are shown in Figure 2, and the equivalent magnetic circuit of the ITQPMG is shown in Figure 3.
When the generator is working in no-load state, the direct axis component of armature reaction is F d � 0. According to Ohm's law and Kirchhoff's law of magnetic circuits, we can establish the EMC model as follows: where F mT is the equivalent magnetomotive force of tangential PM steel, H c is the coercivity of PM, b mT is the thickness of tangential PM steel in the magnetization direction, F mQ is the equivalent magnetomotive force of radial PM steel, F d is d axis component of the armature reaction, G mT is the permeance of tangential PM steel, G mQ is the permeance of radial PM steel, G Tl is the leakage permeance of tangential PM steel, G Ql is the leakage permeance of radial PM steel, G r1 is the permeance between rotor core and air gap, G r2 is the permeance of rotor core between two radial PM steels, G g is the permeance of air gap between rotor and stator, G t is the permeance of stator tooth, G ts is the permeance of stator boot, G y is the permeance of stator yoke, φ mT is the flux provided by tangential PM steel, φ mQ is the flux provided by radial PM steel, φ Ql is the leakage flux of radial PM steel, φ Tl is the leakage flux of radial PM steel, and φ U is the effective flux.
e calculation formula of the phase's back EMF of generator: where f is the frequency of the generator; K w is the winding coefficient of the armature winding; ), where q is the number of slots per pole of each phase; q � (Z s /2pm), where Z s is the number of stator slots and m is the number of phases, m � 3; K p is the short range factor; K q � sin(180β/2), β � (y/mq), where y is the pitch of winding; K φ is the waveform coefficient of air gap flux, K φ � 1.11; and N is the armature winding turns of each phase.e air gap circumferential flux produced by PM steel in rotor can be expressed in the following calculation formula: Journal of Electrical and Computer Engineering where D is the rotor diameter, α j is the mechanical pole arc coefficient of the generator, α m is the magnetic pole arc coefficient, B r is the remanent flux density, and L is the generator's axial length.e wye connection method is used in this design, so the no-load back EMF is equal to times the rated voltage.erefore, the width of the tangential PM steel is 11 mm, the thickness of the magnetization direction is 4 mm, the width of the radial PM steel is 6 mm, the thickness of the magnetization direction is 2 mm, and the implanting depth of the radial PM is 14 mm.

Optimal Design of Rotor Pole
e size and position of the rotor pole have a significant influence on the performance of the PM generator.However, the parameters of magnetic pole calculated by the EMC method may not be the best performance parameter, so it is necessary to optimize the design of the rotor pole.e Taguchi method is a local optimization algorithm, which can optimize multiple objectives at the same time.By establishing the orthogonal test table, the optimal combination of multiobjective optimization design parameters can be calculated with the least experiment times.In this study, the Taguchi method is used to optimize the rotor poles of ITQPMG.

Test Schemes.
For PM generator with parallel magnetic circuit, the parameters of PM and implanting radial depth of PM have great influence on the performance [14].erefore, there are three optimization objectives in this study, such as the peak value of the cogging torque (T), the peak value of the air gap flux density (G), and the distortion rate of the noload back EMF (K r ).e width of tangential PM (h mT ), the thickness of the magnetization direction of the tangential PM (b mT ), the width of the radial PM (h mQ ), the thickness of the magnetization direction of the radial PM (b mQ ), and the implanting depth of the radial PM (b) are selected as variables.e experimental matrix and finite element results are shown in Table 2.   e calculation formula of back-EMF waveform distortion rate for PM generator [15][16][17]: where U N is the amplitude of Nth harmonics.
According to the factors and levels in Table 2, the experimental matrix of L16(4 5 ) is established.If the traditional single-variable and single-objective optimization method uses 4 5 � 1024 experiments, the optimization design of the multiparameter and multitarget of permanent magnet generator can be completed by the Taguchi method only with 16 experiments.e orthogonal experimental matrix of influencing factors is established, and the experimental results are solved with finite element analysis software.
e results of the experimental orthogonal table and the simulation solution are shown in Table 3.
According to the formula, the average value of the solution in Table 3 is calculated: where n is the experiment times and S i is the average value of the target performance of the ith times.e calculation results are shown in Table 4. en, the average value of each factor level of the parameter to the specific target performance is calculated: where M xi is the average value of the performance index under the ith factor of parameter x, M x is the performance index of parameter x under a certain experiment, and J, k, l, and n are the serial numbers of the experiment.
In order to facilitate the comprehensive analysis, the change of the indicators with the level of factors is represented by broken lines, as shown in Figures 4-6.
From Figures 4-6, it is observed that these three sets of horizontal combinations are designed for optimization of single performance indices.For example, considering the effect of three performance indices on ITQPMG, it is necessary to analyze the variance to further analyze the influence of the changes of each parameter on the different performance indices and to determine the parameters.Changing the proportion of the influence on its performance indices and obtaining the optimization results, we obtain where s is the influence of factors, such as b 1 , h 1 , b 2 , h 2 , and b; S(s) is the variance of performance index under the parameter s; M(s) is the total average value of a performance index.e calculation results of variance reflect the proportion of different level factors to each test index.e result of S(s) is shown in Table 5.

Determination of the Final Optimization Scheme.
As can be seen from Table 5, the size of variance can directly reflect the proportion of the impact of optimization parameters on performance indicators [18,19].e parameter b mT has the greatest influence on the air gap magnetic density peak and the back-EMF waveform of the PM generator.Compared with the proportion of S Kr and S G , the width of tangential PM steel h mT has great influence on the cogging torque of the PM generator.e magnetization direction thickness of the radial PM steel b mQ and the implanting depth of PM steel b have the greatest influence on the peak value of air gap magnetic density.

Journal of Electrical and Computer Engineering
According to the above analysis, the selection of all factors should be based on the optimization standard of the unloaded back EMF-waveform distortion rate, peak value of air gap magnetic density maximum, and peak value of cogging torque minimum.Finally, the parameter combination [h mT (1) b mT (4) h mQ (1) b mQ (3) b( 2)] is determined to be the best combination.e magnetization direction of tangential PM steel is 5 mm, the width of radial rectangular PM steel is 4 mm, the thickness of the magnetization direction of radial rectangular PM steel is 4 mm, and the radial rectangular PM steel is implanted deep.
e implanting depth of PM steel is 12 mm.

Result of Finite Element Analysis
In order to verify the superiority of the optimized PM generator, the initial design structure calculated by the EMC method and the optimized structure parameters obtained by the Taguchi method are compared and analyzed.In the case of the same stator structure, winding mode, rotor diameter, silicon steel sheet, and PM material of the two PM generators mentioned above, both are compared and analyzed by the finite element method.
Figure 7 is the no-load back-EMF waveform of the two cycles before and after optimization.It can be seen that the peak value of the no-load back-EMF wave before optimization has obvious depression and the no-load back-EMF waveform of the optimized ITQPMG is closer to the sinusoidal wave.It is shown that the sinusoidal EMF waveform can be obtained by the optimized design.
Figure 8 shows the amplitude of each harmonic content in the no-load back EMF before and after the optimization.It can be seen that the amplitude of the fundamental wave before optimization is 36.77V, the waveform distortion rate is 45%, the amplitude of the fundamental wave after optimization is 38.63 V, and the waveform distortion rate is 33%.
Figure 9 shows the cogging torque before and after the optimization of the generator.It can be seen that the optimized cogging torque decreases from 1.12 N•m to 0.55 N•m after optimization, which is reduced by about 51%.Journal of Electrical and Computer Engineering Figure 10 shows the variation pattern of the peak value of air gap magnetic density along the axial and radial directions of the rotor before and after optimization.It can be seen that the peak value of air gap magnetic density before optimization is 0.83 T and it appeared at the position of 12 mm in the axial direction and 94 mm in the radial direction, while the optimized peak value of air gap magnetic density is 1.12 T, it appears at the position of 5 mm axial distance and 27 mm circumferential distance.Compared to the optimization, the peak value of air gap magnetic density has been improved by about 35%.
In order to verify the advantages of the ITQPMG, it is compared with the efficiency of two other traditional PM generators, as shown in Figure 11.With the increase of speed, the efficiency of ITQPMG is higher than that of the other two traditional PM generators.
e maximum efficiency of ITPMG is 86.5%, and the maximum efficiency of IRPMG is 88.6%, while the maximum efficiency of ITQPMG can reach 91.2%. is is because ITQPMG not only has less magnetic leakage but also has significant magnet congregate effect.

Performance Test
e photographs are shown in Figure 12, and the generator test system and its principle block are shown in Figure 13.e working process of the generator test system is mainly controlled by IPC to control the speed of frequency motor, so as to change the speed of ITQPMG and electronic load for testing.Data acquisition chip and signal processing unit are used to collect and analyze the signals of voltage, current, and speed.en, the test data are fed back to the electrical control unit to complete the automatic adjustment of the electronic load, thus completing the performance test of the generator.
e experimental platform is driven by a 11 kW variable frequency motor as the drive motor to operate the generator, and the no-load and load experiments are carried out, respectively.
When the ITQPMG prototype is running at the rated speed of 4000 r/min, the no-load back-EMF waveform before and after optimization is shown in Figure 14.It can be seen that the optimized no-load back-EMF waveform is closer to the sine wave, the peak value of the no-load back-EMF is increased, and the distortion rate of the waveform is reduced.e experimental test is basically consistent with the result of the finite element analysis method.
Finally, in the case of load power of 480 W, 500 W, and 520 W, the load performance test of the optimized ITQPMG is carried out and the experimental results are shown in Table 6.
It is found from Table 6 that the output voltage maintains around 26.2 V∼28.6 V for the generator speed of 2000 r/ min∼4800 r/min and load power of 480 W∼520 W, which meets the design requirements.

Conclusion
A new type of rotor topology of PM generator with tangential magnetic circuit and radial magnetic circuit which provide a common air gap magnetic field is developed.is structure has less magnetic flux, lower cogging torque, and significant magnetic congregate effect.e structural parameters of the PM generator are preliminarily determined by the analytical calculation of the EMC method, and the Taguchi method is introduced to optimize the multiobjective design of the tangential and radial parallel magnetic circuits of the PM generator.e optimum combination of ITQPMG is obtained.Combined with the finite element method, the optimization results are compared and analyzed.e cogging torque is reduced by 51%, the air gap magnetic density is increased  by 35%, the peak value of no-load back-EMF waveform is reduced by 12%, and the waveform is closer to the sine.e performance test is performed with various rotary speeds and load powers.When the generator speed changes from 2000 r/min∼4800 r/min and the load power varies from 480 W∼520 W, the output voltage maintains around 26.2 V∼28.6 V, which has an excellent regulator performance.

Figure 2 :
Figure 2: Diagram of magnetic flux of ITQPMG.(1) e main magnetic flux of tangential PM steel.(2) e main magnetic flux of radial PM steel.(3) e leakage flux of tangential PM steel.(4) e leakage flux of tangential PM steel.

Figure 4 :Figure 5 :
Figure 4: Impact of various factors on K r .

Figure 6 :
Figure 6: Impact of various factors on G.

Figure 14 :
Figure 14: Comparison experimental results of no-load back-EMF waveforms.(a) No-load back-EMF waveforms before optimization.(b) No-load back-EMF waveforms after optimization.

Table 1 :
Initial design parameters of the generator.

Table 2 :
Test table of optimization parameters and factor level.

Table 3 :
Experimental matrix and finite element analysis.

Table 4 :
Average value of performance index.

Table 5 :
Variance and proportion of the performance indices of each parameter under 4 level factors.

Table 6 :
Test results of the generator output voltage.