# Analysis of the Rotor Magnetomotive Force of Built-In Radial Permanent Magnet Generator for Vehicle

Authors: Huihui Geng; Xueyi Zhang; Tao Si; Lanian Tong; Qingzhi Ma; Mingjun Xu
Journal: International Journal of Rotating Machinery
 (2021)

---

## Abstract

Permanent magnet generator (PMG) for vehicles has attracted more and more attention because of its high efficiency, high power density, and high reliability. However, the weak main air-gap magnetic field can affect the output performance and the normal use of electrical equipment. Aiming at the problem, this paper took the rotor magnetomotive force (MMF), the direct influencing parameter of the main air-gap magnetic field, as the research object, deduced the analytical expression of rotor MMF of the built-in radial PMG in detail, and analyzed its main influencing factors in analytical expression, including the permanent magnet steel (PMS) material, the thickness of PMS in magnetizing direction, the vertical length of the inner side of PMS, and the effective calculation length of PMS. Based on this, the rotor parameters were optimized to obtain the best values. After that, the finite element simulation and prototype test of the optimized generator were carried out. The comparative analysis results showed that the optimized rotor parameters could effectively improve the rotor MMF and optimize the output performance of the generator.


---

## Body

## 1. Introduction

The performance of the vehicle generators directly affects the performance and service life of the electrical equipment. Permanent magnet generator (PMG) has a broad application prospect because of its advantages of high efficiency, high power density, and high reliability [1–3]. Compared with the built-in tangential PMG, the permanent magnet steel (PMS) of the surface-mounted PMG and the built-in radial PMG is closer to the main air gap, and the distance between the PMSs is closer too, so there is less magnetic leakage between them. Moreover, compared with the surface-mounted PMG, the built-in radial PMG has high structural strength and high reliability and therefore has become the future development trend [4–7]. However, the rotor magnetic field of this generator is provided by two adjacent radial PMSs in series and therefore needs to pass through two PMSs. Due to the large magnetic resistance of the PMS, the rotor magnetic field is weakened, which can lead to the weakening of the air-gap magnetic field and affect the output performance of the generator [8]. Therefore, it is of great significance to accurately analyze the influencing factors of rotor magnetic field and improve the main air-gap magnetic field and the output characteristics of the generator.At present, there are two main methods to study the main air-gap magnetic field. One is the finite element method, which mainly uses the finite element software to analyze the variation law of the main structural parameters and the main air-gap magnetic field parameters, for example, the study [9] used the finite element method to analyze the air-gap magnetic flux density and its Fourier variation parameters of a permanent magnet motor with two rotors, studied the variation law of design parameters and air-gap magnetic flux density, and obtained the optimal stator parameters. In [10], the influence of the coupling effect of permanent magnet linear motor on the air-gap magnetic field was analyzed by finite element method, and the key parameters were analyzed and optimized. The other research method for the main air-gap magnetic field is the equivalent magnetic circuit method, which method analyzes the influence of different structures on the main air-gap magnetic field by establishing a simple equivalent model. For example, in [11], the equivalent magnetic circuit model of permanent magnet linear motor considering magnetic saturation effect was derived. The expression of main magnetic flux was calculated, and the model was verified and calibrated by experiments. Ref. [12] established an equivalent magnetic circuit model by analyzing the change of air-gap reluctance of a switched reluctance motor at different rotor positions and analyzed the output performance of the motor. This method can analyze the magnetic circuit distribution law of different topologies, but most usually use simple reluctance and magnetic conductivity as parameters, and the accuracy is limited. Even if the test method is used for correction as in [10], it cannot adapt to the motors of other sizes.Therefore, some scholars began to study the detailed rotor magnetomotive force (MMF) model and substitute it into the equivalent magnetic circuit model to improve the accuracy. Rotor MMF is the MMF that the rotor magnetic field can provide for the main air gap. It is the main target parameter of rotor magnetic circuit design and is mainly used for the analysis of rotor magnetic field. Rotor MMF was used to calculate the main air-gap magnetic field in the early stage. Some scholars calculated the main air-gap magnetic field as the product of the rotor MMF function and the main air-gap length function [13]. Then, the rotor MMF was used to assist in calculating the electromagnetic parameters and related parameters of the generator, such as main air-gap magnetic flux, induced electromotive force and its harmonics, and mutual inductance between rotor and stator magnetic circuit. For example, Ref. [14] proposed an equivalent magnetic circuit model of a single rotor circuit of a doubly fed induction motor, took the spatial harmonic of rotor MMF as the influence parameter of rotor magnetic field quality into the model, and calculated the spatial processes of self-induction, mutual inductance, and back electromotive force of generator. Ref. [15] studied the changes of rotor MMF and its harmonics with different stacking addenda of PMS in rotor core of the built-in multilayer PMG. Then, it analyzed the stator iron loss law based on the optimization of rotor MMF. In [16], the rotor MMF of hybrid excitation salient pole motor was derived based on the principle of minimum reluctance. Based on this, the magnetic permeability function of the main magnetic field considering nonuniform air gap was calculated by using the equivalent magnetic circuit method. At present, the research on rotor MMF still focuses on the simple equivalent methods of magnetic conductivity and reluctance and less on the analytical calculation of magnetic circuit distribution [17–22]. However, the rotor structure parameters directly affect the magnitude of the rotor MMF, so it is necessary to find out the relationship between them. On this basis, this paper deduced the analytical expression of the rotor MMF of the built-in radial PMG in detail, analyzed the influencing factors in the analytical expression, and took a 3-phase 8-pole 36-slot built-in radial PMG as an example to analyze the influence rules of the rotor parameters on the rotor magnetic field and main air-gap magnetic field.


## 2. Equivalent Magnetic Circuit Analysis and Analytical Expression of Rotor MMF

Taking a 3-phase 8-pole 36-slot radial PMG as an example, its technical indicators are shown in Table1. In this generator, the PMSs are evenly embedded in the rotor core in the radial direction, and to reduce the magnetic leakage at both ends of the PMS, magnetic separation slots are set at both ends. In addition, the rotor is laminated and welded with 0.5 mm silicon steel sheets in order to reduce eddy current loss. The generator structure and the rotor punching are shown in Figure 1, and the magnetic circuit of the main flux and leakage magnetic flux is shown in Figure 2.Table 1
The technical indicators of the generator.
Main parametersParameter valueRated power1000 WRated voltage14 VRated speed4000 r/minCooling modeSelf-cooling fan IC01Working temperature-40°C~75°CFigure 1
The structure of the built-in radial PMG and the rotor punching.
(a)
The structure of PMG(b)
The structure of rotor punchingFigure 2
The magnetic circuit analysis diagram of generator.As shown in Figure2, the main flux of the generator passes through the rotor core of the two poles and the PMS. Then, it passes through the main air gap and the stator core to form a circuit, whereas the leakage magnetic flux mainly passes through the two ends of the PMS and forms a loop through the rotor core. According to the magnetic circuit analysis, the equivalent magnetic circuit diagram is drawn and shown in Figure 3.Figure 3
Equivalent magnetic circuit model of generator.In Figure3, the magnetic field of PMS is equivalent to the series connection of permanent magnetic MMF and internal permeance. And in Figure 3, Fpm1 and Fpm2 are the MMFs generated by the adjacent two pole PMS, respectively. Gpm is the internal permeance of PMS. Gl1, Gl2, and Gl3 are, respectively, the permeance of the magnetic circuit inside the PMS of two adjacent magnetic poles, the permeance of the outer magnetic circuit of the PMS, and the leakage permeance of the outer magnetic circuit of the PMS. Gδ is the permeance of main air gap. Gst is the permeance of the magnetic circuit of the stator. Φpm is total magnetic flux produced by two PMS. Φr is the leakage magnetic flux at the end of PMS. Φu is the effective magnetic flux of the generator, and Fad is the direct axis reaction potential of the generator. Figure 3 shows that the main air-gap magnetic field is mainly related to the MMF of PMS, the internal permeance of PMS, the permeance of the magnetic circuit inside the rotor core, the permeance of the main air gap, the permeance of the magnetic circuit inside the stator, and the direct axis reaction potential of the generator. Among them, due to the application requirements of the generator, the main technical parameters are basically unchanged. That is, the permeance of the main air gap, the permeance of the magnetic circuit inside the stator, and the direct axis reaction potential of the generator are basically unchanged, and the main air-gap magnetic field is mainly related to the magnetic circuit inside the rotor core, the MMF, and the internal permeance of PMS. The magnetic circuit part of all these parameters can be simplified and equivalent to the rotor MMF, and the simplified equivalent magnetic circuit model is shown in Figure 4. Therefore, analyzing and improving the rotor MMF can effectively improve the main air-gap magnetic field and improve the output performance of the generator.Figure 4
Simplified equivalent magnetic circuit model of generator.When the built-in radial PMG expands along the outer diameter of the rotor in the circumferential direction and ignores the influence of the stator slot, the distribution of the rotor MMF is shown in Figure5.Figure 5
The circumferential distribution of rotor MMF.In Figure5, the coordinate axis x is the circumferential length coordinate along the outer diameter of the rotor, the axis y is the radial length coordinate of the generator, F is the rotor MMF in the air gap, Fa is the amplitude of the rotor MMF, tp is the effective calculation length of the PMS, and ts is the rotor pole moment. It can be seen from Figure 5 that the rotor MMF in the main air gap is approximately rectangular due to the leakage magnetic flux at the ends of the PMS and the magnetic flux loss of the rotor core. To simplify the analysis, the rotor MMF is equivalent to the standard rectangle, as shown in the lower coordinate system of Figure 5. Moreover, assuming that the axial direction is infinite and ignoring the axial length of the generator, the rotor MMF with circumferential distribution Fx can be calculated as follows:
(1)Fx0,kts−ts−tp2≤x<kts+ts−tp2,−1kFa,kts+ts−tp2≤x<kts+ts+tp2.The Fourier decomposition ofFx can be obtained as follows:
(2)Fx=∑n=1∞−1n+3/24Fanπsinnπtp2tssinnπtsxn=1,3,5,⋯,where n is the number of harmonics.The amplitude of the rotor MMFFa is related to the magnetic circuit of the generator. For analytical calculation, the internal magnetic circuit of the rotor core is equivalent to a circle, as shown in Figure 2. And it can be seen from Figure 2 and the equivalent magnetic circuit model that the rotor MMF of each pole of the rotor is jointly provided by the N-pole and S-pole PMS and the MMF of two parts through the rotor core to the outer circular of the rotor. Therefore, the rotor MMF is the MMF generated by two PMSs minus the MMF of two parts of the rotor core. And Fa can be calculated as follows:
(3)Fa=2Fpm−2Fr,where Fpm is the MMF on the surface of PMS and Fr is the MMF when the MMF provided by PMS passes through the rotor core.Considering the leakage magnetic flux at the ends of the PMS, the intersection point between the inner radial edge of the PMS and the connecting line between the top vertex of the end of the PMS and the rotor center is taken as the dividing point of the leakage magnetic flux and the main magnetic flux, and it can be seen in Figure2. The part between the dividing points on both sides of PMS is recorded as the effective calculation length of PMS tp. The arc length of the outer diameter of the rotor corresponding to the two dividing points is the effective calculation length of the rotor pole and is recorded as te. The relationship between them is as follows:
(4)te=rarctantp2rc,where r is the outer diameter of the rotor and rc is the distance from the center of the inner edge of PMS to the rotor center and is the vertical length of the inner side of PMS.According to the geometric relationship, the radial effective calculation length of PMStp can be calculated as follows:
(5)tp=rcaprc+bp,where ap is the radial length of PMS and bp is the thickness of PMS in magnetizing direction.Taking one-fourth of the effective radial calculation length of the PMS as the equivalent calculation point of the rotor magnetic potential and according to the geometric relationship, the magnetic circuit lengthlc1 of the rotor core outside the PMS can be calculated as follows:
(6)lc1=rsinarccostp4r−rc−bp.In Figure2, the intersection O of the tangential extension line in the inner side of PMS and the pole division line is taken as the center of the magnetic circuit of the rotor core in the inner side of PMS, and the distance from the center of the inner side of PMS to the center of the magnetic circuit O is recorded as ro; it can be calculated as
(7)ro=rctanπp.The length of magnetic circuitlc2 between the inner calculation points of N-pole PMS and S-pole PMS can be calculated as
(8)lc2=π4pp−24ro−tp.Therefore, the MMF of rotor coreFr can be calculated as
(9)Fr=Brlc1+lc2μr,where Br is the magnetic induction intensity of the rotor core, μr is the permeability of the rotor core, and μr=μr0μ0, where μr0 is the relative permeability of the rotor core, because the rotor core is made of silicon steel laminated, and μr0=8000Hm−1, and μ0 is the vacuum permeability, and μ0=4π×10−7Hm−1.The MMF of PMSFp can be calculated as follows:
(10)Fp=Bpbpμp,where Bp is the magnetic induction intensity of the PMS, μp is the permeability of the PMS, and μp=μp0μ0, where μp0 is the relative permeability of the rotor core; usually, μp0 is 1.05~1.1.Taking the above formula into Equation (3), the amplitude of the rotor MMF Fa can be calculated as follows:
(11)Fa=2Bpbpμp−2Brμrrsinarccosrcap8rrc+bp−rc−bp+π4pp−22rctanπp−rcap2rc+bp.Taking Equation (11) into Equation (2), the Fourier decomposition of Fx is as follows:
(12)Fx=∑n=1∞−1n+3/24nπsinnπtp2tssinnπtsx2Bpbpμp−2Brμrrsinarccosrcap8rrc+bp−rc−bp+π4pp−22rctanπp−rcap2rc+bpn=1,3,5,⋯.According to Equation (12), the main influencing parameters of the rotor MMF of the built-in radial PMG are the magnetic induction intensity of PMS Bp, the thickness of PMS in magnetizing direction bp, the permeability of PMS μp, the magnetic induction intensity of the rotor core Br, the permeability of the rotor core μr, the vertical length of the inner side of PMS rc, the effective calculation length of PMS tp, and the rotor pole moment ts. Among them, the magnetic induction intensity and the permeability are the material parameters, and the rotor core of the generator is usually made of 0.5 mm silicon steel sheet; therefore, only the PMS material is considered. The rotor polar moment is mainly determined by the overall structure of the generator and is not considered temporarily. In conclusion, without considering the overall structural parameters of the generator and the rotor core material, the rotor MMF is mainly affected by the rotor parameters of the PMS materials, hp, rc, and tp.


## 3. Influencing Rules of Rotor Parameters on Rotor MMF

### 3.1. The PMS Materials

There are many kinds of PMS materials used in vehicle generators, and the properties are quite different. Presently, three types of PMS materials are widely used: the ferrite PMS material, the rare earth cobalt PMS material, and the neodymium iron boron (Nd-Fe-B) PMS material [23]. Their main magnetic properties are shown in Table 2.Table 2
The main magnetic properties of PMS materials.
NameModel numberMagnetic induction strength (T)Maximum coercivity (kA/m)Maximum magnetic energy product (kJ/m3)Ferrite PMS materialY300.421030Rare earth cobalt PMS materialXG196/960.96690183Nd-Fe-B PMS materialNdFe351.17868287Taking a 3-phase 8-pole 36-slot built-in radial PMG as an example, the outer diameter of the stator is 140 mm, the inner diameter of the stator is 105 mm, the axial length of the stator is 30 mm, the main air-gap length is 0.4 mm, and the inner diameter of the rotor is 20 mm. Using the three kinds of PMS materials, the generator simulation model is established and simulated. The rotor MMF is calculated as the product of the rotor magnetic flux of each pole and the reluctance of the main magnetic circuit of the rotor. The rotor magnetic flux of each pole can be calculated by the field processor of the finite element software. The calculated results of rotor MMF with different PMS materials are shown in Figure6, and the no-load induction electromotive force (NLIEF) waveforms are shown in Figure 7.Figure 6
The rotor MMF histogram with different PMS materials.Figure 7
The NLIEF waveforms of different PMS materials.Figure6 shows that the rotor MMFs with different PMS materials are quite different. Among them, the maximum value is 367.39 A, which uses the Nd-Fe-B PMS material. A little lower is rare earth cobalt PMS material. And the minimum value is 81.36 A, which uses the ferrite PMS material. It can also be seen from Figure 7 that the maximum value of NLIEF produced by rare earth cobalt PMS material can reach 50 V, which can meet the requirement of rated voltage 14 V. Meanwhile, rare earth cobalt PMS material is cheaper than Nd-Fe-B PMS material and therefore becomes the best PMS material in this paper.


### 3.2. The Thickness of PMS in Magnetizing Directionhp

Using the finite element model, change the thickness of PMS in magnetizing directionhp from 3 mm to 6 mm in steps of 0.5 mm. The variation diagrams of the rotor MMF and the NLIEF waveform are simulated as shown in Figures 8 and 9, respectively.Figure 8
The rotor MMF curve with different thicknesses of PMS in magnetizing direction.Figure 9
The NLIEF waveforms with different thicknesses of PMS in magnetizing direction.Figure8 shows that with the increase of the thickness of PMS in magnetizing direction hp, the rotor MMF increases continuously, but the increase decreases when hp increases to 4.5 mm; it can also be seen in Figure 9, when hp is from 4.5 mm to 6 mm, the NLIEF waveforms change very little. Therefore, considering the rotor MMF and material utilization ratio, it is more appropriate to select the thickness of PMS in magnetizing direction from 4 mm to 4.5 mm.


### 3.3. The Vertical Length of the Inner Side of PMSrc

When the vertical length of the inner side of PMSrc is from 43.2 mm to 45.2 mm and the step is 0.5 mm, the simulation results of the rotor MMF curve and the NLIEF waveform are simulated and shown in Figures 10 and 11, respectively.Figure 10
The rotor MMF curve with different vertical lengths of the inner side of PMS.Figure 11
The NLIEF waveforms with different vertical lengths of the inner side of PMS.As shown in Figures10 and 11, with the increase of vertical length of the inner side of PMS, the rotor MMF and the peak of NLIEF waveform increase. However, if rc is too large, the PMS will be too close to the outer circle of the rotor, and the structural strength of the rotor will be reduced, whereas if rc is too small, the leakage magnetic flux at the ends of the PMS will increase, and the rotor MMF will decrease. When rc is from 44 mm to 44.4 mm, the rotor MMF is from 279.4 A to 309.99 A, and the peak value of NLIEF is about 55 V, which can meet the output requirements. Therefore, it is more appropriate to take rc from 44 mm to 44.4 mm.


### 3.4. The Effective Calculation Length of PMStp

When the effective calculation length of PMStp is from 20 mm to 30 mm and the step is 2 mm, the simulation results of the rotor MMF curve and the NLIEF waveform are shown in Figures 12 and 13, respectively.Figure 12
The rotor MMF curve with different effective calculation lengths of PMS.Figure 13
The NLIEF waveforms with different effective calculation lengths of PMS.As shown in Figures12 and 13, with the increase of the effective calculation length of PMS tp, the rotor MMF increases, and the peak value of the NLIEF waveform increases too, whereas when tp increases to 26 mm, the increasing trend of rotor MMF and the NLIEF becomes slower. And with the increase of tp, the consumption of PMS also increases. So comparing comprehensively, it is more appropriate to select tp from 24 mm to 26 mm.


### 3.5. Multiparameter Optimization Analysis

Taking the thickness of PMS in magnetizing direction, the vertical length of the inner side of PMS, and the effective calculation length of PMS as optimization parameters, the multiparameter optimization analysis of rotor MMF is carried out by finite element software. Each parameter analysis value is set up to refer to the above optimization results, and all parameters are subdivided into 27 groups, as shown in Table3.Table 3
Grouping table of optimization parameters.
Group numberhp (mm)rc (mm)tp (mm)Group numberhp (mm)rc (mm)tp (mm)144424154.344.226244425164.344.424344426174.344.4254444.224184.344.4265444.225194.644246444.226204.644257444.424214.644268444.425224.644.2249444.426234.644.225104.34424244.644.226114.34425254.644.424124.34426264.644.425134.344.224274.644.426144.344.225Three parameters with different values are input into finite element software for parametric simulation, the rotor MMFs of different groups are obtained, and its histogram is drawn in order from big to small, as shown in Figure14.Figure 14
The rotor MMF change histogram of different combinations of rotor parameters.Figure14 shows that the rotor MMF varies with the different combinations of rotor parameters. And the first six groups of high rotor MMF, including group 20, group 21, group 23, group 24, group 26, and group 27, are selected to analyze the NLIEF waveform as shown in Figure 15.Figure 15
The NLIEF waveforms of different combinations of rotor parameters.Figure15 shows that the NLIEF waveforms of group 21, group 24, and group 27 almost overlap, as well as group 20, group 23, and group 26. The peaks of NLIEF waveforms in the second overlap group are slightly lower than those in the first overlap grouping, but the smoothness of the curve increases, the distortion rate at the top of the peak decreases, and the sinusoidal property of the whole curve increases. Therefore, group 20, group 23, and group 26 have higher application values. And the data analysis of these three groups shows that only the vertical lengths of the inner side of PMS are different, which are 44 mm, 44.2 mm, and 44.4 mm, respectively. Because the larger the vertical length of the inner side of PMS, the closer the PMS is to the outer circle of the rotor and the lower the structural strength of the rotor. Therefore, the design parameters of group 26 are selected as the optimal design parameters.


## 3.1. The PMS Materials

There are many kinds of PMS materials used in vehicle generators, and the properties are quite different. Presently, three types of PMS materials are widely used: the ferrite PMS material, the rare earth cobalt PMS material, and the neodymium iron boron (Nd-Fe-B) PMS material [23]. Their main magnetic properties are shown in Table 2.Table 2
The main magnetic properties of PMS materials.
NameModel numberMagnetic induction strength (T)Maximum coercivity (kA/m)Maximum magnetic energy product (kJ/m3)Ferrite PMS materialY300.421030Rare earth cobalt PMS materialXG196/960.96690183Nd-Fe-B PMS materialNdFe351.17868287Taking a 3-phase 8-pole 36-slot built-in radial PMG as an example, the outer diameter of the stator is 140 mm, the inner diameter of the stator is 105 mm, the axial length of the stator is 30 mm, the main air-gap length is 0.4 mm, and the inner diameter of the rotor is 20 mm. Using the three kinds of PMS materials, the generator simulation model is established and simulated. The rotor MMF is calculated as the product of the rotor magnetic flux of each pole and the reluctance of the main magnetic circuit of the rotor. The rotor magnetic flux of each pole can be calculated by the field processor of the finite element software. The calculated results of rotor MMF with different PMS materials are shown in Figure6, and the no-load induction electromotive force (NLIEF) waveforms are shown in Figure 7.Figure 6
The rotor MMF histogram with different PMS materials.Figure 7
The NLIEF waveforms of different PMS materials.Figure6 shows that the rotor MMFs with different PMS materials are quite different. Among them, the maximum value is 367.39 A, which uses the Nd-Fe-B PMS material. A little lower is rare earth cobalt PMS material. And the minimum value is 81.36 A, which uses the ferrite PMS material. It can also be seen from Figure 7 that the maximum value of NLIEF produced by rare earth cobalt PMS material can reach 50 V, which can meet the requirement of rated voltage 14 V. Meanwhile, rare earth cobalt PMS material is cheaper than Nd-Fe-B PMS material and therefore becomes the best PMS material in this paper.


## 3.2. The Thickness of PMS in Magnetizing Directionhp

Using the finite element model, change the thickness of PMS in magnetizing directionhp from 3 mm to 6 mm in steps of 0.5 mm. The variation diagrams of the rotor MMF and the NLIEF waveform are simulated as shown in Figures 8 and 9, respectively.Figure 8
The rotor MMF curve with different thicknesses of PMS in magnetizing direction.Figure 9
The NLIEF waveforms with different thicknesses of PMS in magnetizing direction.Figure8 shows that with the increase of the thickness of PMS in magnetizing direction hp, the rotor MMF increases continuously, but the increase decreases when hp increases to 4.5 mm; it can also be seen in Figure 9, when hp is from 4.5 mm to 6 mm, the NLIEF waveforms change very little. Therefore, considering the rotor MMF and material utilization ratio, it is more appropriate to select the thickness of PMS in magnetizing direction from 4 mm to 4.5 mm.


## 3.3. The Vertical Length of the Inner Side of PMSrc

When the vertical length of the inner side of PMSrc is from 43.2 mm to 45.2 mm and the step is 0.5 mm, the simulation results of the rotor MMF curve and the NLIEF waveform are simulated and shown in Figures 10 and 11, respectively.Figure 10
The rotor MMF curve with different vertical lengths of the inner side of PMS.Figure 11
The NLIEF waveforms with different vertical lengths of the inner side of PMS.As shown in Figures10 and 11, with the increase of vertical length of the inner side of PMS, the rotor MMF and the peak of NLIEF waveform increase. However, if rc is too large, the PMS will be too close to the outer circle of the rotor, and the structural strength of the rotor will be reduced, whereas if rc is too small, the leakage magnetic flux at the ends of the PMS will increase, and the rotor MMF will decrease. When rc is from 44 mm to 44.4 mm, the rotor MMF is from 279.4 A to 309.99 A, and the peak value of NLIEF is about 55 V, which can meet the output requirements. Therefore, it is more appropriate to take rc from 44 mm to 44.4 mm.


## 3.4. The Effective Calculation Length of PMStp

When the effective calculation length of PMStp is from 20 mm to 30 mm and the step is 2 mm, the simulation results of the rotor MMF curve and the NLIEF waveform are shown in Figures 12 and 13, respectively.Figure 12
The rotor MMF curve with different effective calculation lengths of PMS.Figure 13
The NLIEF waveforms with different effective calculation lengths of PMS.As shown in Figures12 and 13, with the increase of the effective calculation length of PMS tp, the rotor MMF increases, and the peak value of the NLIEF waveform increases too, whereas when tp increases to 26 mm, the increasing trend of rotor MMF and the NLIEF becomes slower. And with the increase of tp, the consumption of PMS also increases. So comparing comprehensively, it is more appropriate to select tp from 24 mm to 26 mm.


## 3.5. Multiparameter Optimization Analysis

Taking the thickness of PMS in magnetizing direction, the vertical length of the inner side of PMS, and the effective calculation length of PMS as optimization parameters, the multiparameter optimization analysis of rotor MMF is carried out by finite element software. Each parameter analysis value is set up to refer to the above optimization results, and all parameters are subdivided into 27 groups, as shown in Table3.Table 3
Grouping table of optimization parameters.
Group numberhp (mm)rc (mm)tp (mm)Group numberhp (mm)rc (mm)tp (mm)144424154.344.226244425164.344.424344426174.344.4254444.224184.344.4265444.225194.644246444.226204.644257444.424214.644268444.425224.644.2249444.426234.644.225104.34424244.644.226114.34425254.644.424124.34426264.644.425134.344.224274.644.426144.344.225Three parameters with different values are input into finite element software for parametric simulation, the rotor MMFs of different groups are obtained, and its histogram is drawn in order from big to small, as shown in Figure14.Figure 14
The rotor MMF change histogram of different combinations of rotor parameters.Figure14 shows that the rotor MMF varies with the different combinations of rotor parameters. And the first six groups of high rotor MMF, including group 20, group 21, group 23, group 24, group 26, and group 27, are selected to analyze the NLIEF waveform as shown in Figure 15.Figure 15
The NLIEF waveforms of different combinations of rotor parameters.Figure15 shows that the NLIEF waveforms of group 21, group 24, and group 27 almost overlap, as well as group 20, group 23, and group 26. The peaks of NLIEF waveforms in the second overlap group are slightly lower than those in the first overlap grouping, but the smoothness of the curve increases, the distortion rate at the top of the peak decreases, and the sinusoidal property of the whole curve increases. Therefore, group 20, group 23, and group 26 have higher application values. And the data analysis of these three groups shows that only the vertical lengths of the inner side of PMS are different, which are 44 mm, 44.2 mm, and 44.4 mm, respectively. Because the larger the vertical length of the inner side of PMS, the closer the PMS is to the outer circle of the rotor and the lower the structural strength of the rotor. Therefore, the design parameters of group 26 are selected as the optimal design parameters.


## 4. The Analysis of Optimization Results

After optimization and analysis, the optimized rotor pole parameters are shown in Table4. The finite element simulation model is established and simulated by using the optimized parameters, and the comparison of the NLIEF curves before and after optimization is shown in Figure 16.Table 4
The optimized rotor pole parameters.
Main parametersParameter valuesPMS materialRare earth cobalt PMS material (XG196/96)The thickness of PMS in magnetizing directionhp4.6 mmThe vertical length of the inner side of PMSrc3.6 mmThe effective calculation length of PMStp25 mmFigure 16
The NLIEF waveform before and after optimization.Figure16 shows that the peak value of the NLIEF waveform of the optimized PMG is obviously improved, the smoothness of the peak top is improved, and the sinusoidal property is enhanced. The rotor MMFs before and after optimization are 291.72 A and 321.84 A, respectively, which increases by 10.3% after optimization. This will effectively increase the magnetic flux density in the main air gap and improve the output performance of the generator.To verify the effectiveness of finite element simulation, the prototype is trial produced with the optimized parameters, and the bench test is carried out with the generator performance test bench. The test bench is shown in Figure17, and the comparison of NLIEF waveforms obtained by the test and the finite element simulation is shown in Figure 18.Figure 17
The test bench.Figure 18
The NLIEF waveform obtained by the test and the finite element simulation.It can be seen from Figure18 that the peak value of the NLIEF waveform in the test is slightly lower than the simulation result, but the NLIEF waveform is basically the same. Therefore, the finite element analysis result is effective.


## 5. Conclusions

In this paper, the analytical expression of the rotor MMF of the built-in radial PMG is derived in detail, and the main influence parameters are analyzed. Moreover, the variation rules of rotor MMF with the PMS material, the thickness of PMS in magnetizing direction, the vertical length of the inner side of PMS, and the effective calculation length of PMS are analyzed to optimize the rotor parameters of a 3-phase 8-pole 36-slot built-in radial PMG. After optimization, the rotor MMF increases by 10.3%, the main air-gap flux density increases significantly, the sinusoidal character of NLIEF becomes better, and the output characteristics of the generator improve obviously.



---

*Source: 5319615-2021-11-08.xml*