Design and optimization of Halbach magnetic gear

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Design and Optimization of Halbach Magnetic Gear'''</big></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Peng Wang<span id="fnc-1"></span><span style="text-align: center; font-size: 75%;">[[#fn-1|<sup>1</sup>]]</span>, Junyue Yang2, Yanjun Ge3</div>
-->

'''Abstract''':This article proposes a new structure to enhance the torque and torque density of permanent magnet gears. The outer magnet is designed with a dual-layer permanent magnet configuration, where the first layer adopts a Halbach array, and the other layer utilizes a conventional radial magnetization structure. With the optimization objectives defined as maximum torque and torque density, and guided by parameter sensitivity, optimization analysis was conducted using response surface method and multi-objective genetic algorithm. The optimized parameters were determined, and finite element simulations were performed to validate the optimization results. The result indicate a 22.22% increase in torque and a 27.42% increase in torque density compared to the previous the traditional magnet gear. The performance of the proposed new structure of magnetic gear has shown significant improvement.

'''Keywords''': Mangentic gear; Halbach; RSM; Electromagnetic torque; Torque density

=1. Introduction=

As the demand for energy increases in society, wind energy has gained significant attention as a clean and renewable energy source [1]. The wind turbine gearbox is a crucial component in the doubly-fed wind power generation system, playing a vital role in power transmission. However, the mechanical gearbox is constrained by factors such as the continuous impact of wind direction changes, variable loads, and noise vibrations, making it the component with the highest failure rate in wind turbine generators. While direct-drive systems in wind power generation, which eliminate the need for mechanical gearboxes by directly utilizing the wind turbine blades to drive the generator, simplify the structure, the increasing single-machine power leads to a growing weight and volume of these direct-drive systems. Numerous studies indicate that the drive chain of direct-drive systems is limited by factors such as weight, cost, efficiency, resulting in various shortcomings [2]. Considered as a potential alternative to mechanical gears, magnet gears present several advantages, including low noise, minimal vibration, maintenance-free operation, and built-in overload protection. However, early limitations, such as the parallel-axis topology and lower magnetic material performance, curtailed their ability to achieve higher torque densities, impeding widespread development and adoption. Since the inception of the field modulated permanent magnet gear by D. Howe in 2001[3], there has been a substantial improvement in the performance and application of magnetic materials in permanent magnet gears. This advancement has resulted in a significant upswing in both output torque and torque density. As a result, scholars have introduced a range of topological structures, use comprehensive depiction of the evolving  in a renewed research focus on permanent magnet gear [4-5].

Presently, lots of scholars have contributed noteworthy research findings in the analysis of magnetic gears. In [6], the incorporation of Halbach arrays into magnetic gears is introduced, resulting in a significant augmentation of their torque capabilities. In [7], the modulation effect of the demagnetization block is equated to the air gap magnetic conductivity. This model calculates torque and air gap magnetic density using scalar magnetic potential theory, rendering it more amenable to parameterized optimization computations. Additionally, in [8] puts forth magnetic gears with variable transmission ratios, enhancing their suitability for variable-speed electric vehicles and wind turbines, thereby increasing the practicality of magnetic gears. Correspondingly, in [9] conducts a meticulous examination of the material properties of the modulators in magnet gears, aiming to minimize losses and elevate transmission efficiency. In [10], an optimization design is executed on axial permanent magnet gears utilizing a two-dimensional magnetic equivalent circuit model and a genetic algorithm. This design not only diminishes ineffective harmonics in the air gap, augmenting air gap flux density for heightened torque density but also undergoes optimization through a genetic algorithm.

To enhance the output torque and torque density of magnet gears, this article proposes a Halbach Array enhanced Field Modulated Permanent Magnetic Gear (HFMPMG) based on the Field Modulated Permanent Magnetic Gear (FMPMG). The HFMPMG is parametrically modeled using ANSYS software. The Response Surface Method (RSM) is employed to analyze the optimization of structural parameters for the HFMPMG, resulting in a three-dimensional response surface. Subsequently, the Multi-Objective Genetic Algorithm (MOGA) is applied to obtain the optimal solution. Finally, the proposed HFMPMG is validated through confirmation tests.

=2. Topology and Operation Principle=

Fig. 1 depicts the topological structure of the FMPMG, which is made up of an outer rotor, an inner rotor and modulators which can modulate the inner and outer air-gap.


{|
|-
| [[Image:Draft_Wang_788937672-image1.png|270px]]
| [[Image:Draft_Wang_788937672-image2.png|center|270px]]
|-
|Fig. 1. FMPMG topology 
|Fig. 2.  HFMPMGtopology
|}


In order to efficiently couple the magnetic fields generated by the internal and exterior permanent magnets and achieve torque transfer, the modulators primarily modifies the inner and outer air-gap magnetic fields. The harmonic logarithm of magnetic gear air-gap fields is expressed as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>p_{m,j}=\vert mp+jn\vert \mbox{ }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}


Where <math>\mbox{ }m=1,3,5...j=0,\pm 1,\pm 2,\pm 3,...</math> , ''p'' is the number of pole pairs for the inner or outer PMs , n is the number of modulators.

If the magnetic pole pair number of the inner rotor permanent magnet in FMPMG is  <math>p_1</math> and the corresponding rotational speed is  <math display="inline">{\omega }_1</math> ; the number of modulator is  <math>n_2</math> , and the corresponding rotational speed is  <math>{\omega }_2</math> ; the magnetic pole pair number of the outer rotor permanent magnet is  <math>p_3</math> , and the corresponding rotational speed is  <math>{\omega }_3</math> , the harmonic magnetic field pole pair numbers in the air gap and the corresponding angular velocities of the harmonics are given by:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\omega }_{m,j}=\frac{mp_1}{mp_1+jn_2}{\omega }_1+</math><math>\frac{jn_2}{mp_1+jn_2}{\omega }_2\mbox{   }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}


Where  <math>{\omega }_{m,j}</math> is the angular velocity of the spatial harmonic component.

In coupling the magnetic fields of the inner and outer rotors in the Field Modulated Permanent Magnetic Gear (FMPMG), it becomes imperative to align the radially harmonic component of the modulated air gap flux generated by the outer rotor's permanent magnet with the radial harmonic component produced by the inner rotor's permanent magnet.

According to the theory of magnetic field wave angular velocity, two rotating magnetic fields can only achieve stable speed and torque transmission when the magnetic pole numbers and rotational speeds are equal. Therefore, to couple the magnetic fields of the inner and outer rotors in FMPMG, it is necessary to match the radially harmonic component of the modulated air gap flux produced by the outer rotor permanent magnet with the radial harmonic component produced by the inner rotor permanent magnet. Among all the harmonic magnetic fields, when  <math>m=1,\, j=-1</math> , the synchronous harmonic generated by the inner and outer rotor permanent magnets at the air gap is the largest. Therefore, this harmonic can be considered as the main working harmonic of FMPMG. Other harmonics have relatively smaller amplitudes, and their impact on the average torque of FMPMG is smaller. To achieve the stable operation of the three rotors in FMPMG, the pole pair number of the outer rotor permanent magnet should be represented as  <math>p_3=n_2-p_1</math> .

When  <math>{\omega }_2\mbox{=}0</math> , the gear ratio  <math>G_r</math> is given by:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>G_r=-\frac{p_3}{p_1}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}


This article proposes a improved Halbach Array Field Modulated Permanent Magnetic Gear (HFMPMG), depicted in Fig. 2.

The inner rotor is the same as the conventional FMPMG topology which has one layer of permanent magnets with the arrangement of N-pole-S-pole-N-pole. However, a modification is introduced in the structure of the outer rotor's permanent magnets, employing a double-layer configuration. The first layer adheres to a Halbach array arrangement, while the second layer follows the conventional FMPMG arrangement. This dual-layer design serves dual purposes: it mitigates magnetic flux leakage and adjusts the peak value of radial magnetic density, thereby enhancing the sinusoidal nature of the waveform. Essential structural parameters for the HFMPMG are outlined in Table 1.

<span id='_Hlk153395243'></span><div id="_Hlk153102970" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 1  Design parameters of HFMPMG'''</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Parameter
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Description
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Values
|-
|  style="border-top: 1pt solid black;text-align: center;"|''R''<span style="text-align: center; font-size: 75%;">''i''</span>
|  style="border-top: 1pt solid black;text-align: center;"|Inner radius of inner rotor core /mm
|  style="border-top: 1pt solid black;text-align: center;"|40
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''1''</span>
|  style="text-align: center;"|Thickness of inner rotor core /mm
|  style="text-align: center;"|20
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''2''</span>
|  style="text-align: center;"|Thickness of inner PMs /mm
|  style="text-align: center;"|10
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''3''</span>
|  style="text-align: center;"|Thickness of modulator /mm
|  style="text-align: center;"|10
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''4''</span>
|  style="text-align: center;"|Thickness of outer PMs /mm
|  style="text-align: center;"|10
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''5''</span>
|  style="text-align: center;"|Thickness of outer rotor core /mm
|  style="text-align: center;"|10
|-
|  style="text-align: center;"|''g''<span style="text-align: center; font-size: 75%;">''1''</span>
|  style="text-align: center;"|Thickness of the inner air gap/mm
|  style="text-align: center;"|1
|-
|  style="text-align: center;"|''g''<span style="text-align: center; font-size: 75%;">''2''</span>
|  style="text-align: center;"|Thickness of the outer air gap /mm
|  style="text-align: center;"|1
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|''Al''
|  style="border-bottom: 2pt solid black;text-align: center;"|Axial length /mm
|  style="border-bottom: 2pt solid black;text-align: center;"|60
|}

=2. Multiobjective Optimization of HFMPMG=

Given the both structural disparities, the HFMPMG necessitates the consideration of more parameters in the optimization process, because the parameters significantly influence the output torque and torque density. The optimization goal is also increased from single objective to multiobjective. To obtain accurate design parameters for achieving higher torque and torque density, this article advocates a comprehensive approach that integrates parameter sensitivity analysis with response surface method. Based on the sensitivity analysis results, the MOGA algorithm is employed to optimize the modeling parameters of HFMPMG. Finally, finite element analysis is conducted to compare the performance differences between HFMPMG and FMPMG. The whole process of optimization design is shown in Fig. 3.

<div id="_Hlk153104148" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Wang_788937672-image18.png|216px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 3. Optimization design process of HFMPMG</div>

==2.1 Optimization Objectives and Optimization Parameters==

FMPMG is typically used to transmit mechanical torque, so the magnitude of the output torque per unit volume is a crucial criterion for assessing its performance. Additionally, depending on the application scenario, parameters such as torque peak values and fluctuation are standards for evaluating FMPMG performance. In the optimization of HFMPMG in this article, the primary objectives for optimization will be its maximum output torque and torque density.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Wang_788937672-image19.png|270px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 4. Optimization parameters model</div>


The main optimization parameters of HFMPMG are illustrated in Fig. 4 ''R''<span style="text-align: center; font-size: 75%;">''i''</span>is the inner radius of inner yoke iron, ''L''<span style="text-align: center; font-size: 75%;">''1''</span> is the thickness of inner rotor core, ''L''<span style="text-align: center; font-size: 75%;">''2''</span> is the thickness of inner PMs, ''L''<span style="text-align: center; font-size: 75%;">''3''</span> is the thickness of modulator,''L''<span style="text-align: center; font-size: 75%;">''4''</span>is the thickness of outer PMs, ''L''<span style="text-align: center; font-size: 75%;">''5 ''</span>is the thickness of outer rotor core. <math display="inline">\alpha </math> is the central angle corresponding to the modulator, and  <math display="inline">\beta </math> is the central angle corresponding to the nonmagnetic epoxy resin. The modulator fill factor  <math display="inline">\lambda </math> represents the proportion of the angle occupied by the magnetic silicon steel block and nonmagnetic epoxy resin in the modulator. It can be expressed as:'' ''

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">\lambda \mbox{=}\frac{\alpha }{\alpha \mbox{+}\beta }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}


The optimization range structural parameters for HFMPMG is shown in Table 2.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 2   Parameter optimization range of HFMPMG'''</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Parameter
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Description
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Values
|-
|  style="border-top: 1pt solid black;text-align: center;"|''R''<span style="text-align: center; font-size: 75%;">''i''</span>
|  style="border-top: 1pt solid black;text-align: center;"|Inner radius of inner rotor core/mm
|  style="border-top: 1pt solid black;text-align: center;"|[39-41]
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''1''</span>
|  style="text-align: center;"|Thickness of inner rotor core/mm
|  style="text-align: center;"|[19-21]
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''2''</span>
|  style="text-align: center;"|Thickness of inner PMs/mm
|  style="text-align: center;vertical-align: top;"|[9-11]
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''3''</span>
|  style="text-align: center;"|Thickness of modulator/mm
|  style="text-align: center;vertical-align: top;"|[8-11]
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''4''</span>
|  style="text-align: center;"|Thickness of outer PMs/mm
|  style="text-align: center;vertical-align: top;"|[9-11]
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''5''</span>
|  style="text-align: center;"|Thickness of outer rotor core/mm
|  style="text-align: center;vertical-align: top;"|[9-11]
|-
|  style="border-bottom: 2pt solid black;text-align: center;width: 25%;"|<math display="inline">\lambda </math> 
|  style="border-bottom: 2pt solid black;text-align: center;"|Modulator fill factor
|  style="border-bottom: 2pt solid black;text-align: center;"|[0.4-0.6]
|}

=2. 2 Parameter sensitivity analysis=

The optimization process for HFMPMG encompasses a total of 7 parameters, leading to a geometric expansion in the number of samples. Thus, it becomes imperative to initially scrutinize the influence of optimization parameters on the objectives. Table 2 outlines the impact of each parameter on the optimization objectives. This analysis facilitates the identification of effective ranges for each optimization variable and allows for the prioritization of parameters significantly influencing the optimization objectives. This approach significantly shortens the optimization design time and improves optimization efficiency.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 3  Sensitivity Analysis of HFMPMG'''</div>

{| style="width: 86%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;"|Parameter
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''1''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''2''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''3''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''4''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''5''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|''R''<span style="text-align: center; font-size: 75%;">''i''</span>
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;width: 13%;"|<math display="inline">\lambda </math> 
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;"|Sensitivity
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.02098
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.2351
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.2416
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.1612
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.1485
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.0285
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.1629
|}

=2. 3 Rsm optimization=

Utilizing response surface method for the optimization analysis of magnetic gears does not require explicit expressions of the relationship between various structural parameters and optimization objectives. It only involves fitting a surface equation that closely approximates the actual model and establishing the corresponding response surface model. Based on the response surface method, the response model for HFMPMG is given as:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(x)={\beta }_0+\sum_{i=1}^3{\beta }_ix_i+\sum_{i=1}^3{\beta }_{ii}x_i^2+</math><math>\sum_{i=1}^2\sum_{j\geqslant 1}^3{\beta }_{ij}x_iy_i+</math><math>\epsilon </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}


<span id='_Hlk153353620'></span>Where <math>f(x)</math> is the predicted value of the outer torque or torque density,  <math>{\beta }_0</math> is a constant term;  <math>{\beta }_i</math> is the linear coefficient of the design variable;  <math>{\beta }_{ii}</math> is the quadratic coefficient of structural parameters;  <math>{\beta }_{ij}</math> is the product coefficient of different structural parameters. In addition, and  <math>x_i</math> , <math>x_j</math> are two different structural parameters,  <math>\epsilon </math> is the statistical error.


{|
|-
| [[Image:Draft_Wang_788937672-image34.png|270px]]
| [[Image:Draft_Wang_788937672-image35.png|center|270px]]
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
''              ''</div>

<span id='_Hlk153353591'></span><div id="_Hlk153353597" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(a)'' L''<span style="text-align: center; font-size: 75%;">''2''</span> And'' L''<span style="text-align: center; font-size: 75%;">''4''</span>                                                                    (b)'' '' <math display="inline">\lambda </math> and ''L''<span style="text-align: center; font-size: 75%;">''3''</span>  </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 5. Response of structural parameters to optimization objectives</div>

In accordance with formula (5), it can obtain the results of RSM, which offering a comprehensive depiction of the evolving relationship between the optimization objective and the chosen structural parameters. The outcomes are depicted in Fig. 5. From Fig. 5 (a), it can be seen that as the ''L''<span style="text-align: center; font-size: 75%;">''2''</span> and ''L''<span style="text-align: center; font-size: 75%;">''4''</span> decreases, the torque density also gradually decreases. The decrease in the thickness of the PMs can reduce the volume of HFMPMG, but due to the reduction of the magnetic field coupling strength at the air gap, it will greatly reduce the torque density of HFMPMG. From Fig. 5 (b), it can be seen that when the thickness of the modulator is constant and the modulator fill factor  <math display="inline">\lambda </math> is around 0.5, the torque density also reaches its maximum. When  <math display="inline">\lambda </math> increases or decreases, the torque density decreases. Because when  <math display="inline">\lambda </math> decreases, the magnetic circuit saturates, the magnetic density of the inner and outer air gaps decreases, and the modulator is not significant. However, as  <math display="inline">\lambda </math> increases, the air gap magnetic density gradually increases. When  <math display="inline">\lambda </math> is too large, a large amount of magnetic leakage will lead to a decrease in the performance of HFMPMG. Similarly, when  <math display="inline">\lambda </math> is fixed, the torque density increases as ''L''<span style="text-align: center; font-size: 75%;">''3''</span> decreases. From this, it can be seen that the modulation effect of the modulator in the entire HFMPMG is very obvious. After obtaining the appropriate response surface, the MOGA algorithm was selected for intelligent optimization to determine the final optimization value with high sensitivity, as shown in Table 4.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 4  Optimization results of structural parameters'''</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Parameters
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Description
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Value
|-
|  style="border-top: 1pt solid black;text-align: center;"|''R''<span style="text-align: center; font-size: 75%;">''i''</span>
|  style="border-top: 1pt solid black;text-align: center;"|Inner radius of inner rotor core/mm
|  style="border-top: 1pt solid black;text-align: center;"|39.447
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''1''</span>
|  style="text-align: center;"|Thickness of inner rotor core/mm
|  style="text-align: center;"|19.798
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''2''</span>
|  style="text-align: center;"|Thickness of inner PMs/mm
|  style="text-align: center;vertical-align: top;"|10.441
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''3''</span>
|  style="text-align: center;"|Thickness of modulator/mm
|  style="text-align: center;vertical-align: top;"|9
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''4''</span>
|  style="text-align: center;"|Thickness of outer PMs/mm
|  style="text-align: center;vertical-align: top;"|10.174
|-
|  style="text-align: center;"|''L''<span style="text-align: center; font-size: 75%;">''5''</span>
|  style="text-align: center;"|Thickness of outer rotor core/mm
|  style="text-align: center;vertical-align: top;"|9
|-
|  style="border-bottom: 2pt solid black;text-align: center;width: 25%;"|<math display="inline">\lambda </math> 
|  style="border-bottom: 2pt solid black;text-align: center;"|Modulator fill factor
|  style="border-bottom: 2pt solid black;text-align: center;"|0.555
|}

=3. Performance Analysis and Comparison=

==3.1 Magnetic field characteristics==

After conducting finite element analysis, the gap magnetic field profiles of the FMPMG and HFMPMG are revealed, as depicted in Fig. 6. A careful examination of the magnetic field line comparison diagram elucidates that after modulation both configurations exhibit 8 peaks extending towards the outer rotor and 17 peaks reaching the inner rotor. Notably, these peaks align precisely with the pole pairs of the inner and outer rotors. In comparison to the FMPMG , the HFMPMG used the refined Halbach magnetization method applied to the outer rotor, achieves a heightened ability to guide the radial magnetic permanent magnet's orientation, thereby augmenting the magnetic focusing effect.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[Image:Draft_Wang_788937672-image36.png|246px]]
| [[Image:Draft_Wang_788937672-image37.png|center|246px]]
|}
</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(a) FMPMG                                                                       (b) HFMPMG</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 6. Flux line distributions</div>

Fig. 7 depicts the radial and tangential distribution of magnetic density in the outer air gap for both the FMPMG and HFMPMG, while Fig. 8 presents the Fourier harmonic spectrum of magnetic density within the outer air gap. A clear insight gleaned from both Fig. 7 and Fig. 8 reveals a noteworthy trend: the presence of the double-layer outer permanent magnet in the HFMPMG substantially amplifies the peak values of radial effective harmonics, specifically harmonics 4 and 17. This augmentation holds the promise of fostering an increase in torque, consequently enhancing torque density.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[File:Draft_Wang_788937672_1281_fig7a.jpg|282px]]
| [[File:Draft_Wang_788937672_8064_fig7b.jpg|center|282px]]
|}
</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(a) Radial component                                        (b) Tangential component</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 7. Flux density in outer air gap</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[File:Draft_Wang_788937672_4790_fig8a.jpg|282px]]
| [[File:Draft_Wang_788937672_5037_fig8b.jpg|center|282px]]
|}
(a) Radial component                                        (b) Tangential component</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 8. Harmonic spectra of flux density in outer air gap</div>

==3.2 Torque Characteristics==

In accordance with the formula (3) provided above, when the modulator remains stationary, the inner rotor and outer rotor exhibit motion in opposite directions in accordance with the ratio of their pole numbers. In the case of the HFMPMG, specific speeds are assigned to the inner and outer rotors, set at 170rpm and 40rpm, respectively. Employing finite element software simulation facilitates the derivation of stable torque output waveforms for both the inner and outer rotors.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
| [[File:Draft_Wang_788937672_2288_fig9a.jpg|282px]]
| [[File:Draft_Wang_788937672_9015_fig9b.jpg|center|282px]]
|}
</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(a) Dynamic torque curve                                      (b) Static torque curve</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 9. Torque characteristics</div>

Examining Fig. 9 (a), it is evident that both the FMPMG and HFMPMG demonstrate the capability to generate nearly linear torque waveforms during motion. However, following structural optimization and refinement, the output torque of the HFMPMG reaches 297N·m , surpassing the output torque of the FMPMG, which stands at 243N·m. This represents a notable increase of 22.22%. Similarly, as illustrated in Fig. 9(b), the static torque curves of both the FMPMG and HFMPMG exhibit sine characteristics consistent with a transmission ratio of 4.25. In summation, based on the insights derived from both figures, it becomes apparent that the performance of the HFMPMG outpaces that of the FMPMG.

==3.3 Torque density==

To elucidate the electromagnetic characteristics of the HFMPMG, we conducted calculations comparing the volume, PM mass, maximum output torque and torque density of both HFMPMG and FMPMG in Table 5. It can be observed that HFMPMG utilizes a smaller volume and less permanent magnet material, yet it achieves an increase in output torque, resulting in a higher torque density.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 5  Comparison of the characteristics of HFMPMG and FMPMG'''</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Parameter
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|HFMPMG
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|FMPMG
|-
|  style="border-top: 1pt solid black;text-align: center;"|Magnetic gear volume ( <math>\mbox{L}^3</math> )
|  style="border-top: 1pt solid black;text-align: center;"|1.879
|  style="border-top: 1pt solid black;text-align: center;"|1.961
|-
|  style="text-align: center;"|Magnets Weight ( <math>\mbox{kg}</math> )
|  style="text-align: center;"|4.01
|  style="text-align: center;"|4.24
|-
|  style="text-align: center;"|Outer torque ( <math>\mbox{N}\cdot \mbox{m}</math> )
|  style="text-align: center;"|297
|  style="text-align: center;"|243
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|Torque density ( <math>\mbox{N}\cdot m{\mbox{/L}}^3</math> )
|  style="border-bottom: 2pt solid black;text-align: center;"|158
|  style="border-bottom: 2pt solid black;text-align: center;"|124
|}

=4. Conclusion=

(1) This article introduces an enhanced Halbach Array Field Modulated Permanent Magnetic Gear (HFMPMG). Its outer magnet is designed with a dual-layer permanent magnet configuration, where the first layer adopts a Halbach array, and the other layer utilizes a conventional radial magnetization structure.

(2) An exhaustive optimization analysis of critical structural parameters for the HFMPMG is undertaken employing a methodology encompassing parameter sensitivity analysis, response surface analysis, and a multi-objective genetic algorithm. This systematic approach ensures the swift, efficient, and precise attainment of optimized parameters.

(3) The results of magnetic and torque characteristics between the HFMPMG and FMPMG unequivocally affirm that HFMPMG demonstrates superior output torque and higher torque density. The outer torque increased by 22.22% and torque density increased by 27.42%. respectively. thereby validating the rationale and correctness of the optimized parameter acquisition process.

==REFERENCES==

[1]. Kjaer A B, Korsgaard S, Nielsen S S, et al, Design fabrication, test, and benchmark of a magnetically geared permanent magnet generator for wind power generation. IEEE Transactions on Energy Conversion, 35(2019), No.1, 24-32.

[2]. Padinharu D K K, Li G J, Zhu Z Q, et al, Permanent magnet vernier machines for direct-drive offshore wind power: Benefits and Challenges[J], IEEE Access, 10(2022), 20652-20668.

[3]. Atallah K, Howe D., A novel high-performance magnetic gear[J]. IEEE Transactions on magnetics, 37(2001), No.4, 2844-2846.

[4]. Ge Yanjun, Yuan Zhi, Zhao Peng, et al, Modeling and Analysis on Starting Characteristics of CPMGs [J]. CHINA MECHANICAL ENGINEERING, 29(2018), No.13, 1513.

[5]. Jing Libing, Zhang Ting, Lin Ying, et al, A New Method to Improve the Magnetic Gear Torque [J]. Mechanical Transmission, 3(2018), 144-148.

[6]. Mafi H, Afjei E, Ghaheri A, Design of coaxial magnetic gear utilising a novel permanent magnet Halbach array structure[J]. IET Electric Power Applications, 15(2021), No.3, 299-309.

[7]. Shin H M, Chang J H, Analytical magnetic field calculation of coaxial magnetic gear with flux concentrating rotor[J]. IEEE Transactions on Magnetics, 52(2015), No.7, 1-4.

[8]. Chen M, Chau K T, LEE H T, et al, Design and analysis of a new axial-field magnetic variable gear using pole-changing permanent magnets[J]. Progress In Electromagnetics Research, 153(2015), 23-32.

[9]. Yanghui LU, Shuai LUO,Sujun WU, Influence of Material Properties of Magnetic Adjusting Ring on the Loss of Magnetic Gear [J]. JOURNAL OF MECHANICAL ENGINEERING, 52(2023), No.22, 269-275.

[10]. Ruiz-Ponce G, Arjona M A, Hernandez C, et al, Design Optimization of an Axial Flux Magnetic Gear by Using Reluctance Network Modeling and Genetic Algorithm[J]. Energies, 16(2023), No.4, 1852.

----

1  Lect., Dalian Jiaotong University, Dalian, China, e-mail: [mailto:ppmoon@163.com ppmoon@163.com]

2 Prof., Dalian Jiaotong University, Dalian, China, e-mail: [mailto:xyk40@126.com xyk40@126.com]

3 Prof., Dalian Jiaotong University, Dalian, China, e-mail: [mailto:geyanjun1964@126.com geyanjun1964@126.com]