Next Article in Journal
Research on Underground Track Tracking of a New Type Shaft Boring Machine
Next Article in Special Issue
Numerical Modeling and Structure Optimization for Magnetic Levitation Planar Machine Using PCB Coils
Previous Article in Journal
A Multi-Source Domain Adaptation Method for Bearing Fault Diagnosis with Dynamically Similarity Guidance on Incomplete Data
Previous Article in Special Issue
Linear Model Predictive Control and Back-Propagation Controller for Single-Point Magnetic Levitation with Different Gap Levitation and Back-Propagation Offline Iteration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Optimization of Support and Drive System for Magnetic Levitation Air Compressor for Fuel Cells

1
School of Electromechanical and Automotive Engineering, Yantai University, Yantai 264005, China
2
Tianji Kinetic Energy (Beijing) Maglev Technology Development Co., Ltd., Beijing 102600, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Actuators 2025, 14(1), 26; https://doi.org/10.3390/act14010026
Submission received: 11 November 2024 / Revised: 2 January 2025 / Accepted: 8 January 2025 / Published: 13 January 2025
(This article belongs to the Special Issue Actuators in Magnetic Levitation Technology and Vibration Control)

Abstract

:
The 5-degree-of-freedom active magnetic bearings (5-DOF AMB) and high-speed permanent magnet synchronous motor (HPMSM) were combined and applied to energy-recovery-type air compressors for fuel cells, which gives full play to the advantages of both and meets the design requirements for air compressors in fuel cells. Based on the energy recovery air compressor for fuel cells with a power of 30 kW and a rated speed of 100,000 rpm, this paper combined 5-DOF AMB with HPMSM and used it as its support and drive system. Multi-physics field and multi-objective optimization were carried out by integrating the multi-physics field with the Multi-objective Grey Wolf Algorithm (MOGWO), and the feasibility of the design of the system and its reliability were verified using finite element software.

Graphical Abstract

1. Introduction

At present, proton exchange membrane fuel cells (PEMFCs) are most widely used in fuel cell vehicles due to their high efficiency, zero emission, low noise, and high energy density [1]. Among them, the fuel cell compressor is one of the core components of the fuel cell engine air supply system, which is called the “lung” of the fuel cell vehicle. It is used to supply gas to the cathode of PEMFCs, and its performance directly affects the whole system’s efficiency [2,3]. However, it is found that the air compressor consumes a large amount of power, which accounts for about 80% of the power consumption of the hydrogen fuel cell auxiliary system, 20–30% of the total output of the fuel cell stack, and 20% of the cost of the hydrogen fuel cell system [4,5,6]. This not only reduces the efficiency of the overall fuel cell system but also increases the cost of the hydrogen fuel cell system, which hinders the wider application of PEMFCs.
Relevant studies show that the exhaust gas produced by the electrochemical reaction of PEMFCs has the characteristics of a low temperature, and medium and high pressure. The temperature can generally reach 80–120 °C, and the pressure is about 80~90% of the PEMFC inlet pressure (i.e., the outlet pressure of the air compressor), which still has a high energy recovery value [7,8,9]. Therefore, experts at home and abroad have carried out in-depth research on how to recover energy from PEMFC exhaust characteristics. At present, energy recovery air compressors using expanders are considered to be an ideal solution to reduce the power consumption of fuel cell air compressors [10,11].
In addition, fuel cell compressors have the merits of high efficiency and power density, good dynamic response, and the ability to achieve ultra-high rotor speeds aside from not producing any oil contamination [12] Therefore, traditional mechanical bearings and low-speed motors are increasingly unable to meet the requirements of high speed, high efficiency, and low power consumption required by future fuel cell air compressors. In this regard, various experts at home and abroad have carried out in-depth research on the support and drive system of fuel cell air compressors. Among them, HPMSM has outstanding advantages such as oil-free operation, high power density, high efficiency, and large high-efficiency speed range, making it the best choice among the drive motors for fuel cell air compressors [13]. As a component supporting the stable and high-speed operation of the rotor system of fuel cell compressors, the bearing plays an essential role in ensuring the stability and safety of fuel cell compressors. The magnetic bearings and gas bearings, which are frictionless and do not require an oil lubrication system, are just suitable for fuel cell compressors.
The gas bearing uses the high-pressure gas film between the bearing and the shaft neck to support the rotor. According to the principle, it can be categorized into dynamic pressure gas bearings and static pressure gas bearings. The former will cause serious wear and tear during low-speed stages such as the start-stop stage, which will affect the service life of bearings. In addition, dynamic pressure gas bearings also have problems such as instability, limited bearing capacity, and high requirements for bearing materials and machining accuracy during high-speed stages [14]. For the latter, the bearing capacity of the hydrostatic gas bearing is greater than that of the dynamic gas bearing, but it requires an air supply system to ensure stability and reliability of the gas bearing, and its structure and processing are more complicated [15,16].
The magnetic bearing uses the magnetic field force to levitate the rotor in the air, which has the characteristics of no friction, low energy consumption, low noise, long service life, no need for lubrication, no oil pollution, and adjustable dynamic performance [17,18]. Its bearing capacity is much higher than that of gas bearings, which can meet the requirements of complex application scenarios. Moreover, it is very convenient to maintain, repair, and replace. Therefore, magnetic bearings are considered the ideal solution for high-speed rotor support, showing unparalleled advantages in a variety of applications [19,20]. According to the different degrees of freedom and bearing capacity, magnetic bearings can be divided into five degrees of freedom active magnetic bearings (5-DOF AMB), passive magnetic bearings, and hybrid magnetic bearings [21]. Among them, the 5-DOF AMB is characterized by active controllability and has the potential to become a key unit in future intelligent machinery with a broad application prospect [19]. Therefore, this paper takes 5-DOF AMB as the object of study. All the magnetic bearings involved below refer to 5-DOF AMBs unless otherwise specified.
The combination of 5-DOF AMB and HPMSM in air compressors, which gives full play to the advantages of both, meets the design requirements of the air compressor in the fuel cell [22]. Although the performance of the 5-DOF AMB is closely related to the design of the control system, the design of its mechanical structure also has a crucial impact on it [23]. From the perspective of increasing the bearing capacity against vibration and shock, the larger the axial length and volume of the 5-DOF AMB, the better. However, increasing the axial length of the 5-DOF AMB will increase the axial length of the motor rotor, reduce the critical rotational speed of the rotor, and increase the windage loss [24]. Therefore, to improve the performance index, the structural dimensions of the 5-DOF AMB should be optimally designed. With the development of computer technology and the proposal of various intelligent algorithms, a variety of intelligent optimization algorithms have been successfully applied to the optimization design process of 5-DOF AMB.
Zhang SS et al. [25] and Smirnov, A et al. [26] optimized the overall structural parameters of 5-DOF AMB based on the Multi-objective Genetic Algorithm (MOGA), and the optimized simulation results showed that the bearing capacity of 5-DOF AMB was greatly improved compared with that before optimization. Santosh N et al. took the minimum coil weight and copper loss as the objective functions and used MOGA to optimize the structure of the RMB. After optimization, the coil weight of the RMB decreased by 37.19%, while the power loss increased by 43.72% [27]. Zhou J et al. [28] and Zhou Y et al. [29] used the multi-island genetic algorithm (MIGA) to optimize the structural design of a 5-DOF AMB with the objectives of maximizing the electromagnetic force and minimizing the rotor amplitude, respectively. As a result, the bearing capacity of the TMB was improved by 11.88% and the imbalance amplitude was reduced by nearly 50% compared with the pre-optimization period. Wang et al. used the fruit fly optimization algorithm to perform multi-variate and multi-objective optimization for the bearing capacity, volume, and axial length of the radial magnetic bearing (RMB) in the 5-DOF AMB; through the optimization, the bearing capacity of the RMB was increased by 50%, and the axial length and volume were reduced by 30.6% and 19.3%, respectively [30]. Moreover, optimization based on the NGA-II algorithm was implemented by Wang et al. [31], taking maximum first-order intrinsic frequency and minimum volume as the objective function; by Feng C [32] from Shandong University of Technology, taking maximum electromagnetism and minimum volume as the objective function; and by Cao Z [19] from Southeast University, taking minimum copper loss and volume as the objective function. Through these optimizations, the first-order intrinsic frequency of the 5-DOF AMB was increased from 443.14 Hz to 623.91 Hz; the minimum space volume of the thrust magnetic bearing (TMB) and radial electromagnetic bearing (RMB) was reduced by 9.37% and 16.58%, respectively; and the copper loss of the TMB and RMB was reduced by 49% and 29.1%, respectively. Betancor J et al. minimized the volume of RMBs by comparing the Genetic Algorithm (GA) and Pattern Search (PS) methods. After optimization, the volume was reduced by about 35%. The PS method has a larger diameter but a shorter bearing length as compared to GA. Nonetheless, GA produces thicker AMBs and smaller outer diameters [33]. Yadav et al. optimized the 5-DOF AMB based on the heat transfer search (HTS) algorithm with the minimum volume as the optimization goal. The study found that the volume of the HTS-optimized 5-DOF AMB was reduced by about 23% compared with the volume of the 5-DOF AMB obtained by more popular optimization techniques, such as Genetic Algorithm (GA) and Pattern Search (PS) [34].
The 5-DOF AMB is a typical mechatronic product, which involves a variety of technical fields such as mechanical engineering, electromagnetism, material science, rotor dynamics, control engineering, and computer science, and its design is relatively complex [19]. Therefore, compared with the wide application of gas bearings in fuel cell compressors, the application of 5-DOF AMB in fuel cell compressors has received little attention at home and abroad, and there exists a research gap in this aspect. This paper introduces the electromagnetic design of the magnetic levitation air compressor support and drive system for fuel cells and its algorithmic optimization process and points out some comprehensive considerations in its system design. The contributions of this paper are summarized below:
  • Combining 5-DOF AMB with HPMSM and realizing integrated design and application in fuel cell air compressors fills the research gap in China.
  • The multi-physics optimization design of HPMSM was carried out. After optimization, the end length of HPMSM was reduced by 37%. The critical speed of the rotor system has been greatly improved.
  • Based on the MOGWO algorithm, the structural parameters of 5-DOF AMB are optimized with the maximization of bearing capacity and the minimization of volume and drag loss as the objective functions. After optimization, the volume of RMBs is reduced by 36.4%, the drag loss is reduced by 14.3%, and the bearing capacity is increased by 10%. The volume of the TMB has been reduced by 17.5%, the drag loss has been reduced by 12.5%, and the load-bearing capacity has been increased by 4%.
  • Through the finite element software, the feasibility and reliability of the system design are ensured from the aspects of electromagnetism and rotor dynamics.
The content of this paper is structurally organized as follows: Section 2 describes the integrated design of the support and drive scheme for fuel cell air compressors and its integrated optimization design of HPMSM considering multiple physical fields such as rotor dynamics; Section 3 illustrates the detailed design of the initial theoretical parameters of the 5-DOF AMB based on the actual requirements of fuel cell air compressor; Section 4 defines and compares its optimization objective with the structural parameter optimization process based on the MOGWO algorithm; and Section 5 incorporates extensive simulation analysis and calibration with finite element software to ensure the feasibility of its system design and its reliability. Finally, Section 6 summarizes the whole paper. This paper can provide some theoretical value and practical significance for the application of 5-DOF AMB in fuel cell systems.

2. Support and Drive Scheme and Their HPMSM Design

2.1. Support and Drive Scheme Design

The principle of the energy-recovery-type air compressor for fuel cells is as follows:
As shown in Figure 1, the outside air enters the compressor from the inlet at the left and is discharged from the outlet after being compressed. Then, it goes through the intercooler and humidifier to enter the fuel cell stack for electrochemical reaction. The electrochemical reaction at the fuel cell stack produces low-temperature, and medium- and high-pressure exhaust gas [7,8,9]. The exhaust gas goes through the humidifier and enters from the inlet of the expander, where the energy of the exhaust gas is recovered via the expansion of the expander turbine and is finally discharged from the outlet of the expander [35].
The energy-recovery-type magnetic levitation air compressor initially designed in this study has its internal rotor system shown in Figure 2. As shown in the figure, the compressor impeller, high-speed motor, and expander turbine are connected coaxially and supported by the 5-DOF AMB; the HPMSM is used to drive the rotation of the rotor shaft, which is placed in the center; the compressor impeller and the expander turbine are located at both ends of the main shaft to balance the axial forces so that the bearing capacity required by the TMB can be decreased, thereby reducing the volume of the TMB and the diameter of its thrust disc. In addition, eddy current sensors are adopted for the system, and it is necessary to retain a certain safety distance among the sensors to prevent electromagnetic interference from affecting the sensor accuracy.

2.2. HPMSM Electromagnetic Design and Analysis

2.2.1. HPMSM Design Considering Multi-Physics Fields

The energy recovery magnetic levitation air compressor system for fuel cells is a typical mechatronic system [19], so it requires comprehensive consideration and optimized design. The mass of the rotor system and the external disturbances to which it is subjected determine the selection of the maximum bearing capacity of the 5-DOF AMB, and the diameter of the rotor journal of the motor determines the selection of the inner diameter of the RMB rotor and the inner diameter of the TMB thrust disc. Meanwhile, the design of the 5-DOF AMB system is also inseparable from the dynamic characteristics of the rotor of the supported object [36]. Therefore, it is necessary to optimize the design of HPMSM before designing 5-DOF AMB.
Figure 3 shows the schematic diagram of the power and recovery efficiency of the energy recovery air compressor system designed by our research group [35]. It can be seen from the figure that both the power of the compressor and expander and the energy recovery efficiency of the energy recovery efficiency increase with the increase in rotational speed. The energy recovery efficiency of the expander reached its maximum at a speed of 100,000 rpm and the power required for the energy-recovery-type air compressor in this study reached its maximum at 26.58 kW. Considering factors such as overload, the rated power of the motor in this study was set at 30 kW, and a high-power density design was adopted. The main design indexes of HPMSM in this study are shown in Table 1.
HPMSM can have an impact on the 5-DOF AMB system. In particular, the volume of the 5-DOF AMB system is much larger than that of the air bearing, thus affecting the critical speed of the rotor as well as the reliability and safety of the system as a result. Considering this factor, before the algorithmic optimization of 5-DOF AMB, multi-physics field analysis was carried out for the optimized design of HPMSM to ensure the safety and reliability of the whole system, as shown in Figure 4.
The specific analysis is as follows.
(1)
Stator–rotor structure design considering electromagnetic performance: As shown in Figure 5, the 1-pair pole rotor structure suitable for HPMSM has 1-pair solid and ring-shaped permanent magnet structures, whose magnetizing directions are parallel with each other. Both structures have their advantages and disadvantages. This study chose a solid structure. The solid type uses a whole-piece permanent magnet, which has high structural strength and a smaller radial length of HPMSM and is easy to fabricate industrially [37,38]. In terms of the stator structure of HPMSM, the stator groove structure is usually divided into three types: no groove, few grooves, and multiple grooves [39]. For the HPMSM with a pole pair of 1 designed in this paper, the coordination of integer slots is mostly adopted, and to reduce harmonics, the multi-slot structure is generally adopted. This also contributes to the reduction in the peak cogging torque and the magnetic flux density ripple on the rotor sheath surface. However, the increase in the number of slots will reduce the width of the stator teeth and the area in the slots, which is not conducive to the offline operation of the stator winding, affects the structural strength, and increases the manufacturing cost [39,40]. Therefore, combined with the parallel magnetization method of solid permanent magnet in this paper, the slotted stator structure with 2 poles and 24 slots was finally selected.
(2)
Selection of permanent magnet and core material considering heat generation: To ensure the performance of HPMSM, reduce the size of the motor, and avoid demagnetization at high temperatures, the permanent magnet material with high magnetic energy product and high-temperature resistance was preferred. Given the high operating temperature in this paper, a samarium cobalt (SmCo30, which is manufactured by Hangzhou Permanent Magnet Group Co., Ltd., Hangzhou, China) permanent magnet was used. Meanwhile, the high frequency of HPMSM leads to high eddy current loss. Therefore, it is necessary to use core materials with excellent properties such as ultra-thinness, high permeability, high resistivity, and low loss coefficient to effectively suppress the eddy current loss in the core [41]. The 0.1 mm 10JNEX900 silicon steel sheet (10JNEX900, which is manufactured by Kawasaki in Tokyo, Japan) has a lower specific iron loss, but the material is too thin and not easy to process. Considering that the structure of the stator punching piece in this paper is complex and the tooth width is small, the 0.2 mm B20AT1500 silicon steel sheet (B20AT1500, which is manufactured by China Baowu Iron & Steel Group Co., Shanghai, China) is selected as the core material. The main advantages are a low unit loss value, easy machining, and a high tensile strength [42].
(3)
Dimensional design and sheath material selection considering rotor strength: Unlike normal speed motors, the selection of a rotor outer diameter and length for HPMSM was prioritized to ensure that the permanent magnets are not damaged by centrifugal tension. A preliminary estimation of the rotor outer diameter can be made from the following equations and in conjunction with the electromagnetic design [43,44].
F = m v 2 / r = A ρ r 2 ω 2  
σ t = F / A = ρ r 2 ω 2
σ t σ t / S  
v m a x = σ t / S ρ  
D m a x = 2   v m a x / ω  
where m is the mass of the permanent magnet rotor, v is the rotor linear velocity, r is the radius of the rotor, A is the outer surface area of the material, ρ is the density of the material, ω is the angular velocity of the rotation, σ t is the allowable stress of the material, S is the safety margin. v m a x is the maximum linear velocity of the rotor surface, and D m a x is the maximum outer diameter of the rotor. In addition, to ensure reliable operation of the rotor, it is necessary to add a sheath outside the permanent magnet to protect the permanent magnet. The commonly used rotor sheaths are mainly made of carbon fiber composite materials, metal materials, etc. Since the rotor structure in this paper is a solid permanent magnet structure, the alloy material GH4169(GH4169, which is manufactured by China’s Jiangxi Province Baoshunchang Special Alloy Manufacturing Co., Xinyu, China) with higher strength was selected to stably connect the spindle and permanent magnets.
(4)
Winding design considering rotor dynamics: As shown in Figure 6, the stator winding suitable for HPMSM mainly has two structures: double-layer short-pitch distributed winding and toroidal winding [45]. In particular, when the number of pole pairs is 1, the toroidal winding can significantly reduce the axial length of the rotor relative to the double-layer short-pitch distributed winding, thereby increasing the critical speed of the rotor system [37,45]. Table 2 shows the end length values of double-layer short-pitch distributed winding and toroidal winding. From the table, it can be seen that the end length is 27 mm for the double-layer short-pitch distributed winding scheme, while the end length is approximately equal to the slot depth of the motor (17 mm) for the toroidal winding scheme, which is a reduction of 37% compared to the former.
Combined with the multi-physics field analysis in the previous sections, the structure and materials of the motor stator and rotor were selected to finalize the electromagnetic scheme for the HPMSM with a rated speed of 100,000 r/min and a rated power of 30 kW, as required in this paper. Some of its parameters are shown in Table 3.

2.2.2. Finite Element Analysis of HPMSM

The 2D model was built according to the motor design scheme in Section 2.2.1, and the motor rotor shaft was omitted due to the use of solid permanent magnets in this paper, as shown in Figure 7.
To verify the electromagnetic performance of the design scheme, the performance of the motor under no-load and load conditions is analyzed in this section, and the main indicators include back electromotive force under no-load conditions, cogging torque, torque under load conditions, magnetic dense cloud diagram, and magnetic field line distribution diagram, as shown in Figure 8, Figure 9, Figure 10 and Figure 11.
As can be seen from Figure 8 above, the no-load back EMF symmetry of HPMSM is good, the peak value of line back EMF is close to 417 V, close to the bus DC voltage of 470 V, which meets the design requirements of the motor, and the back EMF sinusoidality of the 2-pole motor is generally very good, and its harmonic distortion rate is negligible. Figure 9 below is a load torque waveform diagram, the average value of torque is 2.8755 N·m, the maximum value is 2.8821 N·m, and the minimum value is 2.8665 N·m at steady state. The torque fluctuation is 0.31%, which meets the design requirements. Due to the grooving of HPMSM, there will be more or less cogging torque in HPMSM, as shown in Figure 10 below, the maximum value of cogging torque is 0.00172 N·m, which is much less than the value of rated torque, that is, the influence of cogging torque on the motor is very small.
Figure 11 shows the magnetic dense cloud and magnetic field line distribution map of HPMSM under the rated load. It is not difficult to find that most of the magnetic density under the load condition of HPMSM is about 1.3 T, and only the stator yoke exceeds 1.8 T. At the same time, the magnetic flux leakage at the stator notch is small, and the position of the winding is biased towards the bottom of the groove, which helps to inhibit the skin effect of the wire and reduce the AC loss of the winding. In addition, the large gap at the top of the stator slot contributes to the heat dissipation of HPMSM. In summary, the electromagnetic characteristics of the motor meet the technical requirements.

3. Initial Theoretical Design of 5-DOF AMB

The goal of the electromagnetic design and optimization of the 5-DOF AMB is as follows: The designed 5-DOF AMB is expected to provide sufficient support for the rotor and resist external disturbances. Thus, the volume and loss should be reduced as much as possible. Figure 12 shows the electromagnetic design and optimization flow of the 5-DOF AMB in this paper.

3.1. Theoretical Design of RMB

(1)
Determination of main design indexes of RMB
Among them, the static load of 5-DOF AMB mainly comes from the gravity of the rotor itself, while the dynamic load mainly comes from the unbalanced vibration of the rotor and the disturbance of the external load. The static bearing capacity of the 5-DOF AMB should not be less than F s gravity of the rotor itself, while the dynamic bearing capacity of F d should include the centrifugal force caused by the unbalanced vibration of the rotor and the impact force caused by external bumps and other shocks [9,19,46]. Therefore, the required bearing capacity of the RMB needs to meet the following:
  F = F s + F d m g + m e Ω 2 + m a  
  e = G / Ω = 9549 G / n  
where F is the required bearing capacity of the RMB; F s is static bearing capacity; F d is the dynamic carrying capacity; m is the mass of the rotor system; e is the allowable amount of residual unbalance; a is the impact acceleration, here taken as 5 g; G is the rotor dynamic balance level; and n is the rotor speed.
For example, according to the rotor balancing requirements, expanders, compressors, and generator rotors with speeds above 950 rpm have a dynamic balancing accuracy class of G2.5 [40]. When the rotor dynamic balance class is G2.5, for the rotor with a mass of 1.6 kg and a rated speed of 100,000 rpm, e , F s , and F d can be calculated as 23.87 μm, 15.68 N , and 120.29 N , respectively, and the required bearing capacity in the radial direction should not be less than 136 N. According to the design experience, the maximum bearing capacity of two RMBs needs to meet the following:
F m a x = 1 2 × F
Taking a safety margin of about 1.5, the total maximum load capacity required for two RMBs is approximately 200 N, that is, the maximum load capacity required for a single RMB should be not less than 100 N. In addition, according to the HPMSM rotor diameter, the outer rotor diameter of the RMB was taken as 35 mm. The input design indexes for the electromagnetic design of the RMB in this paper are shown in Table 4:
(2)
Selection of RMB structure, and stator and rotor core materials
Compared with the homopolar, the anisotropic RMB has smaller circumferential dimensions and a simple structure and is easy to process, and its application is more widespread [19,47]. Among them, the distribution of magnetic poles of the heteropolar structure can be categorized into NSSN and NSNS. When the rotor rotates, the direction of the magnetic lines of force distributed by the NSSN changes fewer times than that of the NSNS in one cycle, which results in smaller losses. Therefore, it becomes a better choice. Therefore, in this paper, the magnetic pole NSSN distribution of the heteropolar radial magnetic levitation bearing is selected. In addition, the core material of a magnetic levitation-bearing stator rotor needs to have the characteristics of high saturation magnetic density, high relative permeability, low iron dissipation, easy processing, and high mechanical strength [19]. A B20AT1500 silicon steel sheet as the core material better meets these requirements, so in this paper, radial magnetic levitation bearing stator rotor core is made by stacking B20AT1500 silicon steel sheets.
(3)
Initial theoretical parameter design
In the actual design process, the number of magnetic poles of the RMB and the value of the air gap between the stator and the rotor were selected based on the outer diameter of the rotor [48], while taking into consideration the limitations of machining accuracy, sensor accuracy, and control system performance [19]. In this paper, the number of magnetic poles of RMBs was selected as 8 and the air gap value was selected as 0.4 mm.
The geometric dimensions of RMB are shown in Figure 13, and the following Formulas (9)–(22) were used to give theoretical calculations for the structural parameters in the figure.
D 2 = D + 2 c
  P = π D 2 / 2 N p  
A = μ 0 F m a x / ( B m a x 2 cos α )  
W p = 1 1.5 × P
W r = P
D 1 = D 2 W r  
  L r = A / P  
N 1 2 B m a x c / ( μ 0 I m a x )  
d c 4 I m a x / ( π J )
A c u = N π d c 2 / ( 4 η )  
H = D 2 2 + 8 N p A c u / π D 2 / 2  
  D 3 = D 2 + 2 H  
  D 4 = D 3 + 2   W p  
I 0 = 0.5 I m a x  
where   D 2 is the inner diameter of the stator; D is the inner diameter of the stator, which was taken as 35 mm; P is the magnetic pole width; c is the unilateral air gap value, which was taken as 0.4 mm; N p is the number of magnetic poles, which was taken as 8; A is the magnetic pole area; μ 0 is the air permeability; α is the coefficient of the force caused by the magnetic pole dispersion, and usually, 8 poles are obtained, i.e., α = π / 8 ; and B max   is the saturated magnetic density of the core material. According to the design experience, the saturated magnetic density at the air gap was taken as 1 T; F max   is the maximum bearing capacity required for RMB; W p is the width of the stator yoke; W r is the width of the rotor yoke;   D 1 is the inner diameter of the rotor; L r is the thickness of the stator; N is the number of coil windings in a pair of magnetic poles; d c is the minimum wire diameter of the coil winding wire; I m a x is the maximum allowable current; and J is the maximum allowable current density in the conductor under certain heat dissipation conditions. Under natural air cooling conditions, the value is J generally taken as 4 8   A / m m 2 ; A c u is the area of the stator coil cavity; η is the slot full factor of the coil winding, which was taken as 0.6; H is the stator window depth; D 3 is the medium diameter of the stator; D 4 is the outer diameter of the stator; and I 0 is the bias current.
The structural parameters of the RMB were calculated as shown in Table 5.

3.2. Theoretical Design of TMB

(1)
Determination of main design indexes of TMB
The rotor system of the energy-recovery-type magnetic levitation air compressor for fuel cells should not only prevent the rotor from touching the protection bearing in the radial direction but also prevent the rotor from generating a large displacement in the axial direction during operation. The axial displacement of the rotor is mainly caused by the aerodynamic force generated by the impeller impacting the airflow in high-speed rotation. The energy-recovery-type magnetic levitation air compressor designed in this paper has its compressor and expander to balance out part of the axial force, but the remaining axial force still needs to be balanced by the thrust disc and TMB.
Figure 14 [35] shows a schematic diagram of the impeller under force, where the impeller inlet wheel cover is subjected to the pressure of the gas flow and generates a gas force F 1 directed towards the root of the impeller. At the same time, part of the airflow at the impeller outlet will flow into the back gap of the blade, generating a gas force F 2 directed towards the top of the blade.
After analysis and calculation, the overall axial force distribution of the rotor system of the energy-recovery-type magnetic levitation air compressor is shown in Figure 15.
According to the calculation, when the air compressor operates at the design speed after the axial force of the main shaft is balanced by the compressor and the expander, the total remaining axial force of the shaft system is 75.82 N [35]. The force that points to the compressor from the expander is borne by the thrust disc near the compressor side. The TMB is mounted on both sides of the thrust disk to share the applied axial force and ensure that the shaft system does not experience large axial displacements. In addition, according to the design experience, taking twice the safety margin, the maximum load capacity required for the TMB is 152 N.
At the same time, considering the strength and dynamics problems caused by the large diameter of the thrust disc, the diameter of the thrust disc is limited to ≤80 mm. The input design indexes of the TMB electromagnetic design in this paper are shown in Table 6.
(2)
Selection of TMB structure, and stator and rotor core materials
Under the condition of the same structural parameters, the heteropolar TMB has serious magnetic circuit coupling and more magnetic leakage [19]. Therefore, in this paper, the homopolar TMB was used. Electrotechnical pure iron has the advantages of good electromagnetic characteristics, easy access to materials, good toughness, convenient processing, inexpensiveness, and so on, but the disadvantage is low strength. Considering that the rotor system speed is up to 100,000 rpm, it cannot be used as the material for the thrust disc of the TMB in this paper. Because of the processing difficulty, mechanical strength, and loss of the materials used in the TMB, Armco Pure Iron (99.85% Fe) was chosen for the stator core of the TMB, and 4340 was used for the thrust disc.
(3)
Initial theoretical parameter design
The geometric dimensions of TMB are shown in Figure 16. Similarly, the initial design parameters can first be determined by the following empirical equations [19,42].
  A 1   = A 2   = A 0  
A 0 = μ 0 F max   / B m a x 2  
s 7 c  
D s 1 = D + 2 s  
D s 2 = 4 A 0 / π + D s 1 2  
N 2 2 B m a x c / ( μ 0 I m a x )  
d c 4 I m a x / ( π J )  
A w = N π d c 2 / 4 η  
H = H 1 = H 2 = A w  
D s 3 = D s 2 + 2 H
  D s 4 = 4 A 0 / π + D s 3 2  
h A 0 / ( π D s 2 )  
L = 2 h  
D s 5   D s 4  
I 0 = 0.5 I m a x
where A1 and A2 are the magnetic pole areas of the inner and outer rings of TMB, which are equal; μ 0 is the air permeability; F max   is the maximum bearing capacity required by TMB, which was taken as 152 N; B max   is the saturation magnetic density of the iron core material, and the saturation magnetic density at the air gap was set to be 0.8 T according to design experience; s is the distance between the TMB stator and the rotating shaft; c is the single-side air gap value, which was taken as 0.4 mm; D s 1 is the inner diameter of the stator; D is the outer diameter of the rotor, which was taken as 35 mm; D s 2 is the inner diameter of the stator winding slot; N 2 is the number of coil windings; I m a x is the maximum allowable current; A w is the area of the stator coil cavity; η is the slot fullness rate of the coil winding, which was taken as 0.6; d c is the minimum diameter of the coil winding wire; J is the maximum current density allowed in the conductor under certain heat dissipation conditions, where under natural air-cooling conditions, the value of is generally taken as 4 8   A / m m 2 ; H 1   and   H 2 are, respectively, the width and height of the coil slot, which are equal; D s 3 is the outer diameter of the stator winding slot; h is the thickness of the yoke; L is the thickness of the thrust disc; D s 4 is the outer diameter of the stator; and D s 5 is the diameter of the thrust disc.
The final structural parameters of TMB. were calculated as shown in Table 7.

4. Structural Parameter Optimization of 5-DOF AMB

4.1. Mogwo Principles

MOGWO was improved by Mirjallii et al. [11] based on the Grey Wolf Algorithm. On the premise of maintaining the excellent performance of the Grey Wolf Algorithm, the external population archive was introduced to store the non-dominated solutions and integrate the grid mechanism to improve the quality of the non-dominated solutions. The algorithm process is as follows:
  D α =   C 1 × X α ( t ) X ( t ) | D β =   | C 2 × X β ( t ) X ( t ) | D δ =   C 3 × X δ ( t ) X ( t ) |    
X 1 = X α ( t ) A 1 × D α X 2 = X β ( t ) A 2 × D β X 3 = X δ ( t ) A 3 × D δ    
  X t + 1 = ( X 1 + X 2 + X 3 ) / 3
A = 2 a × r 1 a , C = 2 r 2  
where   D α ,   D β , and D δ are, respectively, the distances between the individual grey wolves and the three leader wolves, i.e., α ,   β ,   a n d   δ ; C 1 , C 2 , a n d   C 3 represent random vectors; X 1 , X 2 , and X 3 are, respectively, the position vectors under the guidance of α ,   β ,   a n d   δ ; X is the current position of the grey wolves; a is the convergence coefficient, which has an initial value of 2 and decreases to 0 with the increase in the number of iterative steps; and r 1 and r 2 are random numbers with a value range of [0, 1].

4.2. Objective Function

The main function of the 5-DOF AMB is to provide a sufficiently large and controllable electromagnetic suction force to keep the rotor stably suspended and withstand certain vibrations and shocks, so the bearing capacity of the 5-DOF AMB should be as large as possible. At the same time, to increase the critical speed of the rotor of the high-speed motor, the volume of the AMB should be shortened as much as possible [15]. In addition, considering the extremely high rotational speed in this study, where the rated speed reached up to 100,000 rpm and the heat generation caused by air friction loss on the surface of the rotor was very high, the calculation of windage loss cannot be neglected.
In summary, the optimization objectives selected in this paper are to increase the bearing capacity of the 5-DOF AMB and reduce the volume and windage loss of the 5-DOF AMB. Among them, the RMB objective function is shown in Equations (42)–(46). The TMB objective function is shown in Equations (47)–(51).
(1)
Objective function of RMB
  m a x F r m i n V R   m i n P a i r c  
where F r is the electromagnetic force of RMB, V R is the overall volume of the RMB stator–rotor core, and P a i r c is the air gap windage loss caused by the tangential flow generated by the rotation of the RMB rotor. The specific expansion is shown in Equations (43)–(46).
  F r = B m a x 2 / μ 0 cos ( π / 8 )  
  V R = π / 4 D 2 D 1 2 L r + π / 4 D 2 + 2 H + 2 w p 2 D 2 2 8 A c u L  
  P a i r c = k C f 1 π ρ ω 3 D 4 L r  
where k is the surface roughness of the rotor, which is generally taken as 1–1.4, and for the smooth rotor, the value is taken as 1 [40,49,50]; C f 1 is the friction coefficient between the rotor surface and the air; and ρ is the air density, which was taken as 1.29   k g / m 3 . Among them, C f 1 can be estimated by the following empirical Equation (46) [51,52]:
  C f 1 = 0.515 ( c / D ) 0.3 / R e c 0.5 ,   500 < R e c < 10 4 C f 1 = 0.0325 ( c / D ) 0.3 / R e c 0.2 ,   10 4 < R e c R e c = ρ ω D c / μ  
where R e δ is the Reynolds number of the rotor surface, which is related to the flow state of the gas, and μ is the dynamic viscosity coefficient of the air, which was taken as 17.9 × 10 6     P a · s .
(2)
Objective function of TMB
  m a x F a m i n V A m i n P a i r d  
where F a is the TMB electromagnetic force of TMB, V A is the overall volume of the TMB stator–rotor core, and P a i r d is the windage loss generated by the high-speed rotation of the airflow in the gap at the end of the TMB thrust disc. The specific expansion is shown in Equations (48)–(52).
  F a = B m a x 2 A 0 / μ 0  
  V A = π / 4   D s 4 2 D s 1 2 H 2 + h π / 4 D s 3 2 D s 2 2 H 2 + π / 4   D s 4 2 D 2 L
P a i r d = 0.5 C f 2 ρ ω 3   D s 5 5   D 5  
where C f 2 is the friction coefficient between the surface of the thrust disc and air and ρ is the air density, which was taken as 1.29   k g / m 3 . C f 2 can be estimated from the following empirical Equation (51) [51,52,53]:
C f 2 = 64 / ( 3 R e d ) ,   R e d < 30 C f 2 = 3.87 / R e d 0.5 ,   30 < R e d < 3 × 10 5 C f 2 = 0.146 / R e d 0.2 ,   3 × 10 5 < R e d R e d = ρ ω D s 5 2 / μ  
where R e d is the Reynolds number at the end of the thrust disc, which is related to the flow state of the gas, and μ is the dynamic viscosity coefficient of the air, which was taken as 17.9 × 10 6   P a · s .

4.3. Optimization Variables

Based on the objective function determined above, five optimization variables with strong correlation were selected. Among them, the RMB optimization variable is shown in Equation (52). The TMB optimization variable is shown in Equation (53).
X 1 =   P   W p     D 1   H   L r  
  X 2 =   D s 1 D s 2 D s 4   H 2   h    

4.4. Constraints

The range of constraint selection is determined based on the bearing capacity, air gap magnetic density, geometric dimension, and design experience. The RMB constraints are shown in Equation (54) and the TMB constraints are shown in Equation (55).
  100   N   F m a x 150   N D 4 110   m m P a i r c 11   W A 83.3   mm 2 6 P 12.3 P   W p 1.2 P 18   D 1 22 20 H 22 12 L r 14  
152   N   F m a x 200   N P a i r c 60   W A 0 191   mm 2 40.6 D s 1 43 44.6 D s 2 54 D s 4 80   mm 14 H 16 2   h 5  

4.5. Optimization Results

Based on the above multi-objective optimization design process, after 300 iterative operations, a Pareto solution set containing 100 sets of optimal solutions was generated. Figure 17 shows the Pareto optimal front obtained by RMB optimization, and Figure 18 shows the Pareto optimal front obtained by TMB optimization.
Comprehensively taking into consideration the values of the three optimization objectives, a satisfactory solution was selected from the set of Pareto solutions according to the actual needs. The comparison between the optimized design and the initial theoretical design is shown in Table 8 and Table 9.
Compared with the pre-optimization, the volume of the optimized RMB is reduced by 36.4%, the windage loss is reduced by 14.3%, and the bearing capacity is increased by 10%. The optimized TMB has a volume reduction of 17.5%, a windage loss reduction of 12.5%, and a 4% increase in bearing capacity.

5. Multi-Physics Field Analysis of Optimization Results

5.1. Electromagnetic Characterization

To verify whether the 5-DOF AMB model designed and optimized in this paper meets the technical requirements, a finite element simulation analysis of the electromagnetic characteristics was carried out by Maxwell 2020.

5.1.1. Electromagnetic Characterization of RMB

The modeling of the RMB was completed based on the structural parameters calculated in the previous section. The main structures of the model include a stator, a rotor, a coil winding, and a main shaft. The static magnetic field was selected as the simulation environment of the model, and material properties were assigned to each part of the model in Maxwell, in which the stator and rotor were made of silicon steel (model B20AT1500); the spindle was made of 17-4PH and the coil winding was made of copper. After setting the material properties, the model was meshed according to the material of each part. A bias current of 2 A and a maximum current of 4 A were applied to the coil winding. After completing the relevant settings, the calculation was solved and the relevant results were viewed to analyze the magnetic field distribution and bearing capacity characteristics of the RMB, as shown in Figure 19, Figure 20 and Figure 21.
Figure 19a and Figure 19b, respectively, show the magnetic field line distribution and magnetic induction intensity distribution when the RMB coil winding was loaded with the maximum current of 4 A. From Figure 19a, it can be seen that there are four magnetic flux circuits. It is obvious that the magnetic flux density of the RMB is mostly around 1.4 T at the maximum current and exceeds 1.7 T only at a few places. Overall, there is no magnetic saturation, and the design structure is reasonable.
Figure 20 shows the variation curve of magnetic induction intensity at the RMB air gap, where the horizontal coordinate is the length of the air gap line in the circumferential direction. From Figure 20, it can be seen that at the maximum current, the maximum value of magnetic induction strength at the air gap is about 0.93 T, which approximates the design value of 1 T and meets the design requirements.
Figure 21 below shows the three-dimensional electromagnetic characteristic curve of RMB electromagnetic force versus displacement and control current obtained from the Maxwell simulation. As shown in Figure 21, the electromagnetic force that RMB can generate, i.e., the load-carrying capacity, increases with the increase in displacement and control current, which is proportional. In this case, the maximum load-bearing force provided by RMB is 108 N in stable suspension, i.e., when the displacement is 0, which is consistent with the theoretically calculated value of 110 N.

5.1.2. Electromagnetic Characterization of TMB

The finite element analysis method for TMB is the same as that for RMB. The modeling of the TMB is completed based on the structural parameters calculated in the previous section. The main structures of the model include a stator, a thrust disc, and a coil winding. The static magnetic field was selected as the simulation environment for the model and material properties were assigned to each part of the model in Maxwell2020, where the stator was Armco Pure Iron (99.85%Fe), the thrust disc was made of 4340 , and the coil winding was made of copper. After completing the setting of material properties, the model was meshed according to the material of each part, and the coil winding was loaded with a bias current of 2 A and a maximum current of 4 A. After completing the relevant settings, the calculations were solved and the relevant results were viewed to analyze the magnetic field distribution and bearing capacity characteristics of the TMB, as shown in Figure 22, Figure 23 and Figure 24.
Figure 22a and Figure 22b are, respectively, the magnetic field line distribution diagram and the magnetic induction intensity distribution diagram when the TMB coil winding was loaded with the maximum current. As can be seen from Figure 22a, the magnetic field strength of the stator, rotor, and air gap are all less than 1.2 T. Only some of the corners exceed the maximum allowable magnetic induction intensity, so the design structure is reasonable. From Figure 22b, it can be seen that the magnetic field lines are concentrated around the coil winding, and there is some magnetic leakage at the top of the thrust disc and the gap between the TMB and the spindle.
The variation curve of magnetic induction strength at the TMB air gap is shown in Figure 23. It can be seen that the magnetic induction strength at the TMB air gap is 0.75 T at the maximum current of 4 A, which approximates the design value of 0.8 T and meets the design requirements.
Figure 24 below shows the three-dimensional electromagnetic characteristic curves of the TMB electromagnetic force versus displacement and control current obtained from the Maxwell simulations. Similarly, as shown in Figure 24, the electromagnetic force generated by the TMB increases proportionally with the increase in displacement and control current. In particular, at the maximum displacement and control current, the TMB produces a maximum load-bearing force of 246 N. Even in extreme cases, it can still provide sufficient axial equilibrium force without large axial deflections.

5.2. Rotor Dynamics Analysis

When the rotor system of an air compressor operates at high speed, the rotor system is prone to malignant resonance due to the influence of the unbalanced mass of the rotor and external disturbances, which destroys the stable operation of the system and has a serious impact on the safety of the system [54]. To ensure the stable operation of the ultra-high-speed rotor, it is necessary to carry out a detailed rotor dynamics analysis for the air compressor rotor system.

5.2.1. Analysis of Critical Rotor Speed

For a rotor dynamics analysis, the first step is to construct a finite element model of the rotor. In this paper, the dynamics module of Workbench 2021 was used to perform the relevant calculations. After importing the finite element model of the rotor system, the corresponding material properties were assigned to the model, as shown in Table 10.
The model was meshed. The impeller model was relatively complex, for which unstructured tetrahedral mesh was used, and the mesh element was 2 mm. The rest of the structure of the rotor system was relatively regular, so the structured hexahedral mesh was enough, and the mesh element was 4 mm. A total of 141,622 nodes and 75,746 units were obtained. For the contact settings, the actual assembly of the rotor system was simulated, and the contact between the main shaft and the thrust disc was set to be a 0.05 mm interference fit with a friction coefficient of 0.2; the contact between the outer surface of the permanent magnet and the rotor sheath was set to be a 0.03 mm interference fit with a friction coefficient of 0.2; and the bonded type was used for the rest of the parts. The profiled model is shown in Figure 25.
In terms of constraints, the 5-DOF AMB for the rotor support belongs to the elastic-type support, which can be similar to an elastic damped support with stiffness coefficients and damping coefficients for modal analysis in Workbench 2021, and the constraints are generally simulated using spring units, where the damping coefficient of the equivalent spring is small and can usually be ignored [55]. It was easy to obtain from Equations (56) and (57) [56,57] that the support stiffness was 300 N/mm for RMB and 500 N/mm for TMB.
k 1 x = ( 4 μ 0 A N 1 2 I 0 2 c o s α / c 3 )
  k 2 x = ( μ 0 A 0 N 2 2 I 0 2 / c 3 )
where   k 1 x is the displacement stiffness coefficient of RMB, and   k 2 x is the displacement stiffness coefficient of TMB. N 1 is the number of turns of coil winding for RMB and N 2 is the number of turns of coil winding for TMB.
In addition, the rotor will generate a vortex in actual operation. To obtain more accurate calculation results, the gyroscopic effect was turned on in the setting of Workbench for analysis while taking full consideration of the influence of the vortex on the critical speed. Meanwhile, the Campbell diagram was drawn, as shown in Figure 26 and Figure 27.
The rotor resonance mainly results from the unbalanced force of the rotor, which excites the precession mode with the same eddy direction and speed direction, so the intersection between the synchronous excitation line and the forward precession mode frequency curve represents the critical speed of the rotor. As can be seen in Figure 26 and Figure 27, the first-order bending critical speed of the rotor is 145,660 r/min. The safety margin between it and the rated working speed of the rotor is more than 25%, which meets the requirements of general engineering applications [54].

5.2.2. Analysis of Rotor Unbalance Response

High-speed magnetic levitation air compressors for fuel cells have higher requirements on rotor vibration characteristics during operation. Due to the small gap between the 5-DOF AMB, protective bearings and rotor, impeller, and volute, larger rotor vibration displacements will lead to stator–rotor contact wear, which will lead to the instability of the 5-DOF AMB in serious cases, resulting in serious safety accidents [40]. The unbalanced response of the rotor system refers to the rotating mechanical vibration caused by the existence of an unbalanced mass of the rotor. The study of the unbalanced response of the rotor can help to further improve the rotor structure and avoid the vibration impact caused by its structural defects. In addition, the unbalanced response analysis can also provide a theoretical basis for subsequent dynamic balancing experiments and vibration suppression [54,58]. According to the ISO standard, the formula for the unbalanced force due to rotor unbalance is as follows [59]:
F = m e ω 2 = U ω 2  
where F is the unbalanced force; U is the rotor unbalanced; m is the rotor mass; e is the unbalanced degree; and ω is the angular velocity of the rotor.
The allowable imbalance under the G2.5 standard is 0.00038 kg·mm. The frequency range was set to 0~2500 Hz in Workbench to cover the full rotational speed and first-order bending mode conditions, and the number of analysis sub-steps was set to 1000. The number of analysis sub-steps was set to 1000, and an imbalance of 0.00038 kg-mm was applied to the rotor. The vibration displacement distribution of the rotor of the magnetic levitation compressor under different calculation frequency bands obtained through simulation is shown in Figure 28.
As shown in Figure 28, the rotor vibration shows significant peaks at the rigid body and first-order bending mode, and the support section of the protective bearing does not exceed 0.00001 mm.
According to the standard area division in the vibration evaluation standard of magnetic bearings, the vibration displacement is in the A/B area, i.e., the area that is usually considered to be used for a long time without restriction, and the maximum vibration displacement at the minimum clearance is required to be [40]
δ m a x < 0.3 δ min    
where δ min   is the minimum radial clearance between the stator and rotor. The minimum assembly clearance of the magnetic levitation compressor designed in this paper is at the protective bearing, and the radial clearance = 0.2 mm, where the maximum allowable vibration displacement is 0.06 mm according to analysis, which is much larger than 0.00001 mm. Therefore, the vibration displacements of the rotor of the magnetic levitation air compressor designed in this paper can be located in the A/B region in the full speed range (0~2500 Hz), which can ensure the long-term stable operation of the machine.

6. Conclusions

To solve the problems of high energy consumption, short service life of air bearings, and insufficient bearing capacity of fuel cell air compressors in our country, an energy recovery magnetic levitation air compressor with coaxial connection of compressor, HPMSM, and expander and 5-DOF AMB support was proposed. This paper points out some comprehensive considerations in the design of the energy recovery magnetic levitation air compressor support drive system for fuel cells. It also details the electromagnetic design of its system and its algorithmic optimization process, specifically including the following:
(1)
The integrated and optimized design of HPMSM was carried out considering multi-physics fields such as rotor dynamics.
(2)
The initial theoretical design of the 5-DOF AMB was optimized in terms of structural parameters based on MOGWO, with the objective function of maximizing the bearing capacity of the 5-DOF AMB and minimizing its volume and windage loss.
(3)
Through the finite element software, the feasibility and reliability of the system design were verified by analyzing the electromagnetism, rotor dynamics, and unbalanced responses.
Compared with the wide application and research of gas bearings in fuel cell air compressors, the application and research of 5-DOF AMB in fuel cell air compressors have received little attention at home and abroad, and there is a large research gap. Based on the system requirements of the energy recovery magnetic levitation air compressor for fuel cells, this paper integrates multi-physics analysis and optimizes the design based on the MOGWO algorithm, and the optimized results show that the axial length of the HPMSM winding and its related performance of 5-DOF AMB have been greatly improved. Finally, the finite element software analyzes and verifies the feasibility and reliability of its performance and comprehensive optimization design. Through this study, it can provide certain theoretical value and practical significance for the application of 5-DOF AMB in fuel cell systems.

Author Contributions

E.X. contributed to the conception and design of the study, the drafting of the manuscript, and the acquisition of funding. Q.G. contributed to the implementation of the system, the analysis and interpretation of data, and the review and editing of manuscripts. Y.D. contributed to the data collection. W.B. contributed to the software simulation of the data. All authors have read and agreed to the published version of the manuscript.

Funding

Key R&D Program of Shandong Province, Award Number: 2019GGx105015, Recipient: Enhui Xing.

Data Availability Statement

The data underlying this article will be shared on reasonable request to the corresponding author.

Conflicts of Interest

Author Wenxin Bai was employed by the company Tianji Kinetic Energy (Beijing) Maglev Technology Development Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Ma, Y.Y. Application of Proton Exchange Membrane Fuel Cells in Electric Vehicle. Mod. Veh. Power 2024, 2, 6–9+23. [Google Scholar]
  2. Fang, S.Y. Research on Control Method of Air Supply System for Vehicle Fuel Cells. Master’s Thesis, Dalian University of Technology, Dalian, China, 2019. [Google Scholar]
  3. Wu, Y.; Bao, H.; Fu, J.; Wang, X.; Liu, J. Review of recent developments in fuel cell centrifugal air compressor: Comprehensive performance and testing techniques. Int. J. Hydrogen Energy 2023, 48, 32039–32055. [Google Scholar] [CrossRef]
  4. Guo, A.; Li, Q.; Chen, W.R.; Liu, Z.X. Characteristics of Cathode System in Fuel Cells for Electric Vehicles. J. Southwest Jiao Tong Univ. 2013, 48, 1052–1058. [Google Scholar] [CrossRef]
  5. Pei, F.L.; Hou, M.T.; He, J.L.; Wu, B.; Chen, F.X. Modeling of air compressor for proton exchange membrane fuel cell. Energy Storage Sci. Technol. 2019, 8, 1247–1252. [Google Scholar] [CrossRef]
  6. Gao, Z.W.; Li, X.G.; Liu, Z.X.; Rao, J. Development trend and research status of the hydrogen fuel cell vehicle. Mater. Rep. 2022, 36, 74–81. [Google Scholar] [CrossRef]
  7. García, A.; Monsalve-Serrano, J.; Martinez-Boggio, S.; Gaillard, P. Emissions reduction by using e-components in 48 V mild hybrid trucks under dual-mode dual-fuel combustion. Appl. Energy 2021, 299, 117305. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Bao, P.; Wan, Y.; Xu, S. Modeling and analysis of air supply system of polymer electrolyte membrane fuel cell system. Energy Procedia 2017, 142, 1053–1058. [Google Scholar] [CrossRef]
  9. Li, H.Q. Design and Research of High-Speed Centrifugal Compressor for Hydrogen Fuel Cell Vehicle. Master’s Thesis, Zhejiang University, Zhejiang, China.
  10. Yu, W.; Sichuan, X.; Ni, H. Air Compressors for Fuel Cell Vehicles: A Systematic Review. SAE Int. J. Altern. Powertrains 2015, 4, 115–122. [Google Scholar] [CrossRef]
  11. Zhang, Y. Analysis of the development status and trend of air compressors for the fuel cell vehicle. Intern. Combust. Engine Parts 2019, 2, 201–202. [Google Scholar]
  12. Bao, P.L.; Zhang, D.B.; Xu, S.C.; Fang, Y. Development status and trend of air compressor for fuel cell vehicles. Chin. J. Power Sources 2016, 40, 1731–1734. [Google Scholar] [CrossRef]
  13. Meng, H.K. High-speed permanent Magnet Synchronous Motor Control System for Hydrogen Fuel Cell Centrifugal Fans. Master’s Thesis, Jiangsu University of Science and Technology, Zhenjiang, China, 2023. [Google Scholar]
  14. Xiong, W.L.; Wang, J.; Chen, Z.Y.; Sun, W.B.; Xue, H.N.; Tang, X.Q. Review of Research Status and Development of Foil Air Bearings. J. Mech. Eng. 2022, 58, 92–113. [Google Scholar] [CrossRef]
  15. Sun, S. Collaborative Control Study of the Air-Magnetic Suspension Drive System for Fuel Cell Compressors. Master’s Thesis, Yangzhou University, Yangzhou, China, 2023. [Google Scholar]
  16. Peng, C.; He, J.; Zhu, M.; Deng, Z.; Zhen, Z.; Liu, Q. Optimal synchronous vibration control for magnetically suspended centrifugal compressor. Mech. Syst. Signal Process. 2023, 2, 776–789. [Google Scholar] [CrossRef]
  17. Zhou, J.Y.; Zhu, H.Q. Development and Research Summary of Magnetic Bearings. Micromotors 2022, 55, 93–98. [Google Scholar] [CrossRef]
  18. Zhang, J. Structure Design and Performance Analysis of Gas Magnetic Bearing. Master’s Thesis, Yangzhou University, Yangzhou, China, 2020. [Google Scholar]
  19. Cao, Z. Research On Key Problems Of active Magnetic Bearing Aring System for High-Speed Motors. Master’s Thesis, Southeast University, Nanjing, China, 2022. [Google Scholar]
  20. Ravikovich, Y.A.; Ermilov, Y.I.; Kholobtsev, D.P.; Sukhomlinov, I.Y.; Golovin, M.V. Experience of improving rotor system having petal-shaped gas-dynamic bearings for centrifugal refrigeration compressor. Chem. Pet. Eng. 2016, 51, 838–843. [Google Scholar] [CrossRef]
  21. Su, H.N.; Song, J.S.; Zhang, T.X.; Tian, Y.; Liu, Z.H. Research progress of magnetic bearings. J. Shanghai Electr. Technol. 2023, 16, 60–63. [Google Scholar] [CrossRef]
  22. Fang, Q.H. Analysis of Loss and Temperature Field of High-Speed Maglev Permanent Magnet Motor. Master’s Thesis, Wuhan University of Technology, Wuhan, China, 2021. [Google Scholar]
  23. Xiao, L.J.; Zhang, X.S.; Chang, L.; Guo, H. Design and Optimization of Magnetic Bearings Based on ANSYS and ISIGHT. Bearing 2012, 5, 5–8. [Google Scholar] [CrossRef]
  24. Huang, Z.; Han, B. Effective approach for calculating critical speeds of high-speed permanent magnet motor rotor-shaft assemblies. IET Electr. Power Appl. 2015, 9, 628–633. [Google Scholar] [CrossRef]
  25. Zhang, S.S.; Zhou, J.; Zhang, F.P. Structural Optimization Design for Magnetic Bearings Based on Multi-objective Genetic Algorithm. Mach. Electron. 2013, 08, 3–6. [Google Scholar] [CrossRef]
  26. Smirnov, A.; Uzhegov, N.; Sillanpää, T.; Pyrhonen, J.; Pyrhonen, O. High-speed electrical machine with active magnetic bearing system optimization. IEEE Trans. Ind. Electron. 2017, 64, 9876–9885. [Google Scholar] [CrossRef]
  27. Shelke, S.N.; Chalam, R.V. Optimum power loss in eight poles radial magnetic bearing: Multi-objective genetic algorithm. In Proceedings of the International Conference on Computational Intelligence and Information Technology, Surat, India, 7–8 November 2011; Springer: Berlin/Heidelberg, Germany; pp. 72–77. [Google Scholar] [CrossRef]
  28. Zhou, J.; Gao, T.Y.; Dong, J.Y.; Chen, Y.; Gao, S.M. Optimal Design for Structure of Radial Magnetic Bearings Based on ISIGHT. Bearing 2018, 7, 6–11. [Google Scholar]
  29. Zhou, Y.; Zhou, J.; Zhang, Y. Optimum Structural Design of Active Magnetic Bearing Based on RBF Approximation Model. J. Southwest Jiao Tong Univ. 2022, 57, 682–692. [Google Scholar] [CrossRef]
  30. Wang, X.Y.; Zhang, D.Z.; Gao, P. Multi-variable and Multi-objective Optimization Design of Active Radial Magnetic Bearing Based on Fruit Fly Optimization Algorithm. Mot. Technol. 2017, S1, 174–181. [Google Scholar]
  31. Wang, X.; Xu, Y. Multi-Physics Multi-Objective Optimal Design of 5-DOF Magnetic Bearing System with High-speed Permanent Magnet Synchronous Motor. In Proceedings of the 2021 24th International Conference on Electrical Machines and Systems (ICEMS), Gyeongju, Republic of Korea, 31 October–3 November 2021; pp. 457–461. [Google Scholar] [CrossRef]
  32. Feng, C. Optimized Design and Performance Analysis of Active Electromagnetic Bearings for Flywheel Batteries. Master’s Thesis, Shandong University of Technology, Zibo, China, 2023. [Google Scholar]
  33. Betancor, J.; Sahinkaya, M.N.; Zweiri, Y.H. Radial active magnetic bearing design optimization. In Proceedings of the 10th International Conference on Rotor Dynamics–IFToMM: Vol.1 10; Springer International Publishing: Berlin/Heidelberg, Germany, 2019; pp. 321–334. [Google Scholar] [CrossRef]
  34. Yadav, V.K.; Kumar, P.; Bhushan, G. A novel approach for design optimization of radial active magnetic bearing. Aust. J. Mech. Eng. 2023, 21, 1461–1473. [Google Scholar] [CrossRef]
  35. Liu, F.Y. Design and Performance Analysis of Key Components for Fuel Cell Centrifugal Air Compressor with Energy Recovery. Master’s Thesis, Yantai University, Yantai, China, 2024. [Google Scholar]
  36. Fu, Z.S. Research and Design of Active Magnetic Bearing for High-Speed Motors. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2017. [Google Scholar]
  37. Yu, Z.J. Design and Optimization of High-Speed Permanent Magnet Synchronous Motor for Air Compressor. Master’s Thesis, Harbin Institute of Technology, Harbin, China, 2021. [Google Scholar]
  38. Qiu, M.S.; Hua, Q.S. Electromagnetic optimization design of high-speed permanent magnet synchronous motor based on Taguchi iterative method. J. Beijing Norm. Univ. (Nat. Sci.) 2022, 58, 713–719. [Google Scholar] [CrossRef]
  39. Chen, M.X. Design and Magnetocaloric Coupling Analysis of High-Speed Permanent Magnet DRIVE Motor. Master’s Thesis, Shandong University of Technology, Zibo, China, 2023. [Google Scholar]
  40. Gao, T.Y. Research on the Design and Comprehensive Optimization Technology of Maglev High-Speed Permanent Magnet Generator. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2019. [Google Scholar]
  41. Zhu, L.F. Design and Optimization of a Gramme-Ring-Winding High-Speed Slot Less Permanent-Magnet Motor. Master’s Thesis, Nanjing University of Science & Technology, Nanjing, China, 2021. [Google Scholar]
  42. Zhang, M.X. Design and Rotor Loss Optimization of the High-Speed Permanent Magnet Machine for Compressor. Master’s Thesis, Harbin University of Science and Technology, Harbin, China, 2023. [Google Scholar]
  43. Zhong, Y. Design and Research of High-Speed Permanent Magnet Synchronous Motor for Aeration Fan. Master’s Thesis, Harbin University of Science and Technology, Harbin, China, 2022. [Google Scholar]
  44. Shao, Y.Y. Design of Ultra-High-Speed Permanent Magnet Synchronous Motor. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2020. [Google Scholar]
  45. Du, K.L. DESIGN and Experimental Study of Toroidal Winding Flywheel Energy STORAGE motor. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2023. [Google Scholar]
  46. Voigt, A.J. Towards Identification of Rotor Dynamic Properties for Seals in Multiphase Flow Using Active Magnetic Bearings. Design and Commissioning of a Novel Test Facility. Ph. D. Thesis, Technical University of Denmark, Kongens Lyngby, Denmark, 2016. Available online: https://findit.dtu.dk/en/catalog/581cbeaac1ff02154f00008b (accessed on 1 October 2024).
  47. Du, Y.H. Design and Research of Radial Magnetic Bearing Structure for High-Speed Motor. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2021. [Google Scholar]
  48. Wang, Y. Design and Control of Magnetic Bearing for Flywheel Battery. Master’s Thesis, Shandong University of Technology, Zibo, China, 2022. [Google Scholar]
  49. Su, H. Research on Vibration and Control Caused by Thermal Bending of Magnetic Levitation Rotor. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2022. [Google Scholar]
  50. Yan, X. Study on the Mechanism of the Morton Effect of Magnetic Bearings. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2020. [Google Scholar]
  51. Liu, B.; Yu, Z.J.; Fu, J. Cooling Design and Temperature Field Analysis of Radial-axial Magnetic Bearings. Proc. CSEE 2024, 44, 4945–4956. [Google Scholar] [CrossRef]
  52. Xu, Z.Y. Research on Key Technologies of Rotor of High-Speed Permanent Magnet Motor. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2019. [Google Scholar]
  53. Han, B.; Liu, X.; Huang, Z.; Zhang, X.; Zhou, Y. Loss calculation, thermal analysis, and measurement of magnetically suspended PM machine. IEEE Trans. Ind. Electron. 2017, 65, 4514–4523. [Google Scholar] [CrossRef]
  54. Li, H.; Xu, X.B.; Chen, S.; Bi, Z.W. Research on dynamic analysis of three-stage rotor of high-speed magnetic levitation motor. Small Spec. Electr. Mach. 2024, 52, 6–10. [Google Scholar] [CrossRef]
  55. Wang, M.Y. Design and Characteristic Analysis of Ultra-High-Speed Flywheel Energy Storage Rotor and Its Support System. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2022. [Google Scholar]
  56. Hu, Z.D. Research on the Control System of Active Magnetic Bearing. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2020. [Google Scholar]
  57. Yang, J.C. Research on Key Components of Active Magnetic Bearing CONTROL system. Master’s Thesis, Huazhong University of Science and Technology, Wuhan, China, 2022. [Google Scholar]
  58. Wang, C.T.; Guo, H.; Sun, F. Simulation analysis of the automatic balance of unbalanced magnetic levitation rotor. Mech. Eng. Autom. 2021, 5, 60–62. [Google Scholar] [CrossRef]
  59. Ge, Y.Z.; Lei, X.M.; Zhang, Y.C.; Fang, Z.X.; Lu, Z. Analysis and test of the allowable residual unbalance of the rotor. J. Vib. Meas. Diagn. 2023, 43, 371–377+414–415. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of exhaust gas recovery of energy recovery air compressor.
Figure 1. Schematic diagram of exhaust gas recovery of energy recovery air compressor.
Actuators 14 00026 g001
Figure 2. Schematic diagram of rotor system of magnetic levitation air compressor. (a) Overall appearance of the rotor system assembly. (b) Sectional view of the inner part of the rotor system.
Figure 2. Schematic diagram of rotor system of magnetic levitation air compressor. (a) Overall appearance of the rotor system assembly. (b) Sectional view of the inner part of the rotor system.
Actuators 14 00026 g002
Figure 3. Schematic diagram of the power consumption curve of the air compressor. (a) The power of the air compressor system. (b) The recovery efficiency curve of the expander.
Figure 3. Schematic diagram of the power consumption curve of the air compressor. (a) The power of the air compressor system. (b) The recovery efficiency curve of the expander.
Actuators 14 00026 g003
Figure 4. Influence of HPMSM structural design and material selection on multi-physics field.
Figure 4. Influence of HPMSM structural design and material selection on multi-physics field.
Actuators 14 00026 g004
Figure 5. Schematic diagram of rotor structure.
Figure 5. Schematic diagram of rotor structure.
Actuators 14 00026 g005
Figure 6. Ref. [40] Structural diagram of different types of winding. (a) Double-layer short-pitch distributed winding. (b) Toroidal winding.
Figure 6. Ref. [40] Structural diagram of different types of winding. (a) Double-layer short-pitch distributed winding. (b) Toroidal winding.
Actuators 14 00026 g006
Figure 7. HPMSM model.
Figure 7. HPMSM model.
Actuators 14 00026 g007
Figure 8. Back EMF under no-load conditions.
Figure 8. Back EMF under no-load conditions.
Actuators 14 00026 g008
Figure 9. Load torque waveform diagram.
Figure 9. Load torque waveform diagram.
Actuators 14 00026 g009
Figure 10. Cogging torque waveform diagram.
Figure 10. Cogging torque waveform diagram.
Actuators 14 00026 g010
Figure 11. Cloud view of magnetic density and magnetic field line distribution under load conditions. (a) Magnetic flux density. (b) Magnetic field line distribution.
Figure 11. Cloud view of magnetic density and magnetic field line distribution under load conditions. (a) Magnetic flux density. (b) Magnetic field line distribution.
Actuators 14 00026 g011
Figure 12. Electromagnetic design and optimization process of 5-DOF AMB.
Figure 12. Electromagnetic design and optimization process of 5-DOF AMB.
Actuators 14 00026 g012
Figure 13. Schematic diagram of RMB structure.
Figure 13. Schematic diagram of RMB structure.
Actuators 14 00026 g013
Figure 14. Ref. [35]. Schematic diagram of the force on the impeller of the air compressor.
Figure 14. Ref. [35]. Schematic diagram of the force on the impeller of the air compressor.
Actuators 14 00026 g014
Figure 15. Schematic diagram of axial forces of rotor system of magnetic levitation air compressor.
Figure 15. Schematic diagram of axial forces of rotor system of magnetic levitation air compressor.
Actuators 14 00026 g015
Figure 16. Schematic diagram of TMB geometric dimensions.
Figure 16. Schematic diagram of TMB geometric dimensions.
Actuators 14 00026 g016
Figure 17. Pareto optimal solution set of RMB.
Figure 17. Pareto optimal solution set of RMB.
Actuators 14 00026 g017
Figure 18. Pareto optimal solution set of TMB.
Figure 18. Pareto optimal solution set of TMB.
Actuators 14 00026 g018
Figure 19. Cloud view of RMB magnetic flux density and magnetic field line distribution. (a) Magnetic field line distribution of RMB. (b) Magnetic flux density of RMB.
Figure 19. Cloud view of RMB magnetic flux density and magnetic field line distribution. (a) Magnetic field line distribution of RMB. (b) Magnetic flux density of RMB.
Actuators 14 00026 g019
Figure 20. Variation curve of magnetic induction at RMB air gap.
Figure 20. Variation curve of magnetic induction at RMB air gap.
Actuators 14 00026 g020
Figure 21. RMB 3-dimensional electromagnetic characterization curve.
Figure 21. RMB 3-dimensional electromagnetic characterization curve.
Actuators 14 00026 g021
Figure 22. Cloud view of TMB magnetic flux density and magnetic field line distribution. (a) TMB magnetic field line distribution. (b) TMB magnetic flux density.
Figure 22. Cloud view of TMB magnetic flux density and magnetic field line distribution. (a) TMB magnetic field line distribution. (b) TMB magnetic flux density.
Actuators 14 00026 g022
Figure 23. Variation curve of magnetic induction strength at TMB air gap.
Figure 23. Variation curve of magnetic induction strength at TMB air gap.
Actuators 14 00026 g023
Figure 24. TMB 3-dimensional electromagnetic characteristic curve.
Figure 24. TMB 3-dimensional electromagnetic characteristic curve.
Actuators 14 00026 g024
Figure 25. The meshing of air compressor rotor system.
Figure 25. The meshing of air compressor rotor system.
Actuators 14 00026 g025
Figure 26. Air compressor rotor system first-order bending formation (Mode-6-FW).
Figure 26. Air compressor rotor system first-order bending formation (Mode-6-FW).
Actuators 14 00026 g026
Figure 27. Campbell diagram of critical speed of air compressor rotor system.
Figure 27. Campbell diagram of critical speed of air compressor rotor system.
Actuators 14 00026 g027
Figure 28. Rotor system vibration displacement at full speed.
Figure 28. Rotor system vibration displacement at full speed.
Actuators 14 00026 g028
Table 1. Main design indexes of HPMSM.
Table 1. Main design indexes of HPMSM.
ParameterValue
Rated power30 kw
DC bus voltage470 V
Rated rotational speed100,000 rpm
Efficiency>96%
Power factor>0.9
Table 2. End lengths of two winding schemes.
Table 2. End lengths of two winding schemes.
Winding StructureEnd Length (mm)
Double-layer short-pitch distributed winding27
Toroidal winding17
Table 3. The design scheme of HPMSM.
Table 3. The design scheme of HPMSM.
TypeParameterTypeParameter
Permanent   magnet   outer   diameter   ( m m ) 32Stator outer diameter (mm)96
Air   gap   ( m m ) 3.5Axial length (mm)61
Pole pairs1Number of stator slots24
Permanent magnet materialSmco30Stator core materialB20AT1500
Sheath materialGH4169Winding structureToroidal winding
Table 4. Main design indexes of RMBs.
Table 4. Main design indexes of RMBs.
ParameterValue
Maximum   bearing   capacity   F m a x 100   N
Rated rotational speed n 100 , 000   r p m
Winding heat dissipation conditionsNatural air cooling
Outer diameter of stator D4≤110 mm
Outer diameter of rotor D35 mm
Stator thickness Lr≤16 mm
Maximum   allowable   current   I m a x 4 A
Table 5. Initial design parameters of RMB.
Table 5. Initial design parameters of RMB.
ParameterValue
Number   of   poles   N p 8
Inner   diameter   of   rotor   D 1 21 mm
Outer diameter of rotor D 35 mm
Single-side   air   gap   c 0.4 mm
Inner   diameter   of   stator   D 2 35.8 mm
Intermediate   diameter   of   stator   D 3 84.8 mm
Outer   diameter   of   stator   D 4 100.3 mm
Stator   pole   width   P 7.03 mm
Stator   yoke   width   W p 7.73 mm
Stator   thickness   L r 14 mm
Individual   pole   area   A 94.5 mm2
Area   of   sator   coil   cavity   A c u 290 mm2
Number   of   coils   of   a   pair   of   magnetic   poles   N 1 190
Wire   diameter     d c 1.12 mm
Maximum   current   I m a x 4 A
Bias   current   I 0 2 A
Maximum   air   gas   permeability   B max   1.2 T
Maximum   load   capacity   F max   100 N
Table 6. Main design indexes of TMB.
Table 6. Main design indexes of TMB.
ParameterValue
Maximum   bearing   capacity   F m a x ≥152 N
Rated rotational speed n 100 , 000   r p m
Winding heat dissipation conditions Natural air-cooling
Diameter   of   thrust   disc   D s 5 ≤80 mm
Outer diameter of rotor D 35 mm
Maximum   allowable   current   I m a x 4 A
Table 7. Initial design parameters of TMB.
Table 7. Initial design parameters of TMB.
ParameterValue
Outer diameter of rotor D 35 mm
Single-side air gap c 0.4 mm
Inner   diameter   of   stator   D s 1 41 mm
Inner   diameter   of   stator   winding   slot   D s 2 45 mm
Outer   diameter   of   stator   winding   slot   D s 3 74.8 mm
Outer   diameter   of   stator   D s 4 77 mm
Thrust   disc   diameter D s 5 77.5 mm
Winding   slot   depth   H 14.9 mm
Magnetic yoke thickness h2.5 mm
Thrust disc thickness L5 mm
Area   of   stator   magnetic   pole   A 0 264.4 mm2
Area   of   stator   coil   cavity   A w 221.7 mm2
Number   of   the   coils   of   a   pair   of   magnetic   poles   N 2 135
Wire   diameter   d c 1.12 mm
Maximum   current   I m a x 4 A
Bias   current   I 0 2 A
Maximum   air   gap   magnetic   density   B max   0.85 T
Maximum   bearing   capacity   F max   152 N
Table 8. Comparison of initial design parameters with optimized design parameters of RMB.
Table 8. Comparison of initial design parameters with optimized design parameters of RMB.
ParameterInitial ValueOptimized Value
P 7.03 mm8.6 mm
    W p   7.73 mm8.6 mm
      D 1   21 mm22 mm
    H       24.5 mm20 mm
    L r 14 mm12 mm
F r 100   N 110   N
V R 7.26 × 10 4   mm 3 4.62 × 10 4   mm 3
  P a i r c 11.61   W 9.95   W
Table 9. Comparison of initial design parameters with optimized design parameters of TMB.
Table 9. Comparison of initial design parameters with optimized design parameters of TMB.
ParameterInitial ValueOptimized Value
  D s 1 41 mm40.6 mm
D s 2   45 mm44.7 mm
    D s 4   77 mm75 mm
    H 1     14.9 mm14 mm
  h 2.5 mm2 mm
F a 152   N 158   N
V A 3.48 × 10 4   mm 2 2.87 × 10 4   mm 2
  P a i r d 61.5   W 53.8   W
Table 10. Material properties of the rotor system model.
Table 10. Material properties of the rotor system model.
Name of Material Density   ρ /
( k g / m 3 )
Elastic   Modulus   E /
( G P a )
Poisson s   Ratio   μ Application Components
17 - 4   P H 77501960.28Spindle
B 20 A T 1500 77001970.26Radial magnetic bearing rotor core
G H 4169 82401990.3Motor rotor sheath
4340 78402100.278Thrust disc
T C 4 45101100.34Pressure sleeve, expander impeller, lock nut
7075 2810710.33Compressor impeller, inspection ring, baffle
S m C o 84001200.24Permanent magnet
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xing, E.; Gao, Q.; Dong, Y.; Bai, W. Design and Optimization of Support and Drive System for Magnetic Levitation Air Compressor for Fuel Cells. Actuators 2025, 14, 26. https://doi.org/10.3390/act14010026

AMA Style

Xing E, Gao Q, Dong Y, Bai W. Design and Optimization of Support and Drive System for Magnetic Levitation Air Compressor for Fuel Cells. Actuators. 2025; 14(1):26. https://doi.org/10.3390/act14010026

Chicago/Turabian Style

Xing, Enhui, Qi Gao, Yuanqi Dong, and Wenxin Bai. 2025. "Design and Optimization of Support and Drive System for Magnetic Levitation Air Compressor for Fuel Cells" Actuators 14, no. 1: 26. https://doi.org/10.3390/act14010026

APA Style

Xing, E., Gao, Q., Dong, Y., & Bai, W. (2025). Design and Optimization of Support and Drive System for Magnetic Levitation Air Compressor for Fuel Cells. Actuators, 14(1), 26. https://doi.org/10.3390/act14010026

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop