Research on Magnetic Characteristics and Fuzzy PID Control of Electromagnetic Suspension

: This paper proposes an electromagnetic suspension with an electromagnetic actuator, which can improve the riding comfort and stability of the vehicle without changing the safety of traditional MacPherson suspension. First, the electromagnetic suspension structure is introduced, and the principle of the proposed actuator is described in detail. Second, a magnetic ﬂux density model of a single PM ring (permanent magnetic ring) and a magnet assembly are built, and a theoretical analysis of the magnetic ﬂux density is carried out for comparison. Then, the magnetic ﬂux distribution of the magnetic ﬁeld is simulated and analyzed using the ﬁnite element method (FEM), and is compared with theoretical and other experimental data. Finally, a vehicle dynamics model with 7 DOF is built, and vehicle simulations based on the fuzzy PID algorithm are carried out on a C-grade road surface and a deceleration strip. The theoretical results and simulation analyses of the FEM indicate that compared with the MacPherson suspension, the root mean square values of the acceleration of centroid acceleration for the electromagnetic suspension are increased by 59.08% and 33.34%, respectively, on a C-grade road surface and a deceleration strip, and other physical quantities have also been improved. The structure and characteristics of the proposed electromagnetic suspension that improve the riding comfort of the suspension and enhance the stability of the MacPherson suspension are feasible.


Introduction
The rapid growth in the number of vehicles, driven by the development of transportation and the improvement of people's quality of life, has attracted increasing attention on body vibration and tire wear.This concerns not only the safety and integrity of the goods transported by automobiles but also the riding comfort and handling stability of the vehicle [1].As a key component of the vehicle's body, the chassis not only supports the vehicle but also accepts the thrust generated by the engine, ensuring the normal operation of the vehicle [2,3].Suspension systems, as key components of automobile chassis, have garnered significant research interest.
Suspension systems play a key role in supporting the vehicle body and ensuring the comfort, stability, and safety of vehicles.The research and development of comprehensive suspension system performance is a key factor affecting vehicle performance, body direction and attitude control, handling stability, and the driving and riding experience.
McPherson-independent suspension, widely used in the front suspension of modern cars, has several shortcomings, including small brake resistance, weak suspension stiffness, poor stability, and a noticeable turning roll.Therefore, designing a high-performance automobile suspension system can fundamentally improve the performance of an automobile [4,5].Suspension systems are developing towards the direction of lighter weight, smaller volume, faster responses, more stable running state, and greater safety and reliability.With the development of cities and increasing road pressure, vehicles, particularly their suspension systems, are more susceptible to problems [6].The performance of a suspension system mainly depends on the response ability of its actuator or the motor-to-road conditions.Therefore, it is of great significance to research and discuss automobile suspension systems and their related components.
Compared with passive suspension systems, semi-active suspension systems can change the damping of shock absorbers by controlling the input of a single variable so as to adapt to different road conditions [7][8][9].At present, there are two types of adjustable shock absorbers in semi-active suspension systems, which differ in terms of the amount of control they exert.One adjusts the damping size of the damper solenoid valve by adjusting the effective area of the throttle, but this method cannot avoid suspension system noise and stability and reliability problems.The other adjusts the spring rigidity of the whole shock absorber by changing the spring stiffness of the variable spring, but this may produce bumps during the driving process.
Because there are no components with controllable forces installed in the system, a semi-active suspension system cannot control the stiffness of the spring.Therefore, research on semi-active suspension is greatly limited, and the improvement of vehicle performance when using this suspension type is also greatly limited [10].
In recent years, researchers have designed various active suspension systems, including those with rack and pinion rotating motors [11][12][13][14], ball screw mechanisms [15][16][17][18], hydraulic actuating quarter-car active suspension systems [19], and electromagnetic suspension [20][21][22], among others.With the development of electromagnetic control technology and PM material, electromagnetic suspension as a new type of suspension structure has become of interest for researchers, such as traditional hydraulic dampers with a linear motor [23], mechanisms with a ball screw mechanism and two one-way clutches [24], coreless linear tube actuators with a special topology [25,26], electromagnetic and hydraulic suspensions [27], and rotating semi-active electromagnetic variable stiffness devices [28], and so on.They use these mechanisms instead of the dampers or springs of traditional suspensions, or even both can be used.However, safety was not considered in the event of suspension failure or power failure.
Suspension using a motor, magnetorheological fluid damper [29,30], and suspension using an electromagnetic actuator are examples of active electromagnetic suspension.Compared with the first two types of active suspension, electromagnetic actuator suspension is simpler and more efficient in terms of its structure, performance, magnetic force, processing, and so on.An electromagnetic actuator can greatly improve the reaction speed of a suspension system, and thus, it has become the research focus of more and more scholars [31,32].In addition, electromagnetic actuators also have the advantages of simple system configuration, fast response speed, and energy recycling, and as a result, they are more and more widely used.
Integrating control methods into active suspension systems can significantly enhance their performance.In [33], the semi-active suspension of automotive magnetorheological fluids was controlled by a fuzzy control algorithm, and the overall performance of the suspension was improved.In [34], a semi-active 1/4 suspension with three degrees of freedom utilizing a magneto-rheological shock absorber was controlled by a hybrid fuzzy and fuzzy PID controller.In [35], a fuzzy PID controller was applied in a half vehicle active suspension system to enhance the vibration levels of the vehicle chassis and passenger seat.Therefore, the fuzzy PID control algorithm combines the advantages of fuzzy control and PID control, which makes it suitable for research on controlling electromagnetic suspension.This paper aims to make up for the ride comfort and stability of the vehicle without changing the safety of the McPherson suspension.In this paper, we propose an electromagnetic suspension structure that can improve the stability and ride comfort performance of vehicles and make up for the defects of the McPherson suspension.Firstly, the structure of the proposed suspension is introduced, and its principle is described in detail.Secondly, the magnetic flux density model of a single PM ring and a magnet component is established, and the theoretical analysis of the magnetic flux density is carried out.Then, the magnetic flux distribution of the magnetic field is analyzed via the two-dimensional finite element method, and the magnetic flux density is simulated via finite element analysis.Finally, a vehicle dynamics model with seven degrees of freedom is established, and a vehicle simulation on C-grade roads and a deceleration strip is carried out based on the fuzzy PID algorithm.After improving the MacPherson suspension system, active control, four-wheel ride comfort, and stability issues were studied with C-grade roads and a deceleration strip.

Structure and Principle
The proposed structure of the electromagnetic actuator uses a 1/4 MacPherson-type independent suspension system, as shown in Figure 1a.The electromagnetic actuator is inserted into the gap between the hydraulic damper and the spring, the coil assembly of the electromagnetic actuator is connected with the upper cover, and the magnet assembly of the electromagnetic actuator is connected with the cylinder tube of the hydraulic damper, as shown in Figure 1b.The spring, electromagnetic actuator, and hydraulic damper are connected in parallel.The magnet assembly adopts a double-layer PM magnet ring form, and the coil assembly moves in the air gap of the magnet assembly.Compared to the magnet assembly, the coil weighs less; the coils are connected to the upper cover of the shock absorber, which facilitates connection to the vehicle's internal circuit.Consequently, the electromagnetic suspension can be controlled by inputting a current to resolve the disadvantages of the MacPherson-type independent suspension.In this application, it is necessary to maintain a safe distance between the electromagnetic actuator and the spring.When the spring compresses to its lowest point, the distance between the upper cover and lower cover is longer than the axial distance of the electromagnetic actuator, which means that it can avoid the strike of the upper cover and the lower cover at the compression limit.A characteristic of this electromagnetic suspension is that it does not change the structure of the MacPherson suspension.It not only ensures the safety of the MacPherson suspension, but it also improves the stability of the MacPherson suspension.
structure of the proposed suspension is introduced, and its principle is described in Secondly, the magnetic flux density model of a single PM ring and a magnet comp is established, and the theoretical analysis of the magnetic flux density is carried out the magnetic flux distribution of the magnetic field is analyzed via the two-dimen finite element method, and the magnetic flux density is simulated via finite elemen ysis.Finally, a vehicle dynamics model with seven degrees of freedom is establishe a vehicle simulation on C-grade roads and a deceleration strip is carried out based fuzzy PID algorithm.After improving the MacPherson suspension system, active c four-wheel ride comfort, and stability issues were studied with C-grade roads and celeration strip.

Structure and Principle
The proposed structure of the electromagnetic actuator uses a 1/4 MacPherso independent suspension system, as shown in Figure 1a.The electromagnetic actu inserted into the gap between the hydraulic damper and the spring, the coil assem the electromagnetic actuator is connected with the upper cover, and the magnet ass of the electromagnetic actuator is connected with the cylinder tube of the hyd damper, as shown in Figure 1b.The spring, electromagnetic actuator, and hyd damper are connected in parallel.The magnet assembly adopts a double-layer PM m ring form, and the coil assembly moves in the air gap of the magnet assembly.Com to the magnet assembly, the coil weighs less; the coils are connected to the upper co the shock absorber, which facilitates connection to the vehicle's internal circuit.C quently, the electromagnetic suspension can be controlled by inputting a current solve the disadvantages of the MacPherson-type independent suspension.In this ap tion, it is necessary to maintain a safe distance between the electromagnetic actuat the spring.When the spring compresses to its lowest point, the distance between t per cover and lower cover is longer than the axial distance of the electromagnetic ac which means that it can avoid the strike of the upper cover and the lower cover compression limit.A characteristic of this electromagnetic suspension is that it do change the structure of the MacPherson suspension.It not only ensures the safety MacPherson suspension, but it also improves the stability of the MacPherson suspe The structure of the proposed electromagnetic actuator is shown in Figure 2. Th assembly involves a group of coils moving along the axis of the actuator.The magn sembly consists of the outer wall, six small PM rings, and six large PM rings for magnetization, heat-dissipated rings, and a fixed plug.Each PM ring is composed of identical pieces, and their material is Nd-Fe-B, and the PM rings are magnetized as the S pole and the outside as the N pole.The outer wall and the fixed plug are u fix the PM rings.The heat-dissipated ring is used for dissipating heat between tw rings to avoid the phenomenon of heat accumulation.The fixed plug plays a role in the PM rings and heat-dissipated rings.Large PM rings and small PM rings are p side by side; the magnet assembly can generate a quasi-parallel magnetic field.

Design of the Electromagnetic Actuator
For electromagnetic actuators, it is necessary to achieve a compact size and a magnetic force.The influence of the height and radial length of large/small PM ring the height of coils on the magnetic force of the proposed actuator was analyzed orthogonal experimental methods.As shown in Table 1, height and radial leng large/small PM rings and the height of coils were determined as three factors for t The structure of the proposed electromagnetic actuator is shown in Figure 2. The coil assembly involves a group of coils moving along the axis of the actuator.The magnet assembly consists of the outer wall, six small PM rings, and six large PM rings for radial magnetization, heat-dissipated rings, and a fixed plug.Each PM ring is composed of eight identical pieces, and their material is Nd-Fe-B, and the PM rings are magnetized inside as the S pole and the outside as the N pole.The outer wall and the fixed plug are used to fix the PM rings.The heat-dissipated ring is used for dissipating heat between two PM rings to avoid the phenomenon of heat accumulation.The fixed plug plays a role in fixing the PM rings and heat-dissipated rings.Large PM rings and small PM rings are placed side by side; the magnet assembly can generate a quasi-parallel magnetic field.The structure of the proposed electromagnetic actuator is shown in Figure assembly involves a group of coils moving along the axis of the actuator.The m sembly consists of the outer wall, six small PM rings, and six large PM rings magnetization, heat-dissipated rings, and a fixed plug.Each PM ring is compose identical pieces, and their material is Nd-Fe-B, and the PM rings are magneti as the S pole and the outside as the N pole.The outer wall and the fixed plug a fix the PM rings.The heat-dissipated ring is used for dissipating heat betwee rings to avoid the phenomenon of heat accumulation.The fixed plug plays a rol the PM rings and heat-dissipated rings.Large PM rings and small PM rings a side by side; the magnet assembly can generate a quasi-parallel magnetic field.

Design of the Electromagnetic Actuator
For electromagnetic actuators, it is necessary to achieve a compact size a magnetic force.The influence of the height and radial length of large/small PM the height of coils on the magnetic force of the proposed actuator was analy orthogonal experimental methods.As shown in Table 1, height and radial large/small PM rings and the height of coils were determined as three factors thogonal experimental method, with each factor having three levels.The simu rameters are that the coercivity of the PM ring is 868 KA/m and the residual mag

Design of the Electromagnetic Actuator
For electromagnetic actuators, it is necessary to achieve a compact size and a large magnetic force.The influence of the height and radial length of large/small PM rings and the height of coils on the magnetic force of the proposed actuator was analyzed using orthogonal experimental methods.As shown in Table 1, height and radial length of large/small PM rings and the height of coils were determined as three factors for the orthogonal experimental method, with each factor having three levels.The simulation parameters are that the coercivity of the PM ring is 868 KA/m and the residual magnetic flux density of the PM ring is 1.38-1.44T. The L 9 (3 2 × 3 2 ) size optimization simulations of the actuator by the orthogonal experiment method are shown in Table 2.The analysis by the orthogonal experiment method (L 9 (3 2 × 3 2 ) size) is shown in Table 3.As can be observed, the combination (A3B3C2) is the optimal parameters, which are the height of large/small PM rings (12.5 mm), the radial length of large/small PM rings (12 mm), and the height of coils (9.5 mm).Note: K jk (k = 1, 2, 3) is the sum of the simulation results with the same level k in the j th column.K jpk is the mean value of the simulation results with the same level k in the j th column.R j is the range of K jpk .

Sizes and Principles of the Electromagnetic Actuator
According to Section 2.1, the dimension mark and magnetic flux in the air gap are as shown in Figure 3. Figure 3a,b, respectively, show sectional views in the radial direction and in the axial direction.The magnetic flux is emitted from the N pole of the small PM rings through the air gap, the coils, and the coil bobbin.Then, it is gathered at the S level of the large PM ring.This structure for dual PM rings can generate magnetic flux in the same direction as the air gap, which is equivalent to increasing the magnetic focusing effect.Six groups of the dual PM rings generate magnetic flux in the same direction within the air gap.When the coils move in the axial direction, they cause as much magnetic flux as possible.The size parameters of the electromagnetic actuator are shown in Table 4. same direction as the air gap, which is equivalent to increasing the magnetic focusing effect.Six groups of the dual PM rings generate magnetic flux in the same direction within the air gap.When the coils move in the axial direction, they cause as much magnetic flux as possible.The size parameters of the electromagnetic actuator are shown in Table 4.

Model and Analysis of the Magnetic Field
In the electromagnetic actuator, the magnetic field in the air gap of the stator is the main factor for improving performance.The homogeneity of the magnetic field in the air gap along the actuator's axial direction decides the stability and linearity of the actuator's thrust force.The magnetic flux density in the air gap decides the thrust force of an actuator and the ride comfort of a vehicle.First, the model of the magnetic flux density in an air gap is built, and the characteristic of the magnetic field is theoretically analyzed.

Magnetic Flux Density Model of a Single PM Ring
To analyze the quasi-parallel magnetic field in the air gap between the PM ring arrays and the inner wall of the actuator, the mathematical model of the magnetic flux density of the magnet field is built in the y direction according to the Ampere molecular circulation hypothesis, as shown in Figure 4.When the positive direction of the magnet field is +y direction, the spatial magnetic field is equivalent to the magnetic field generated by Surface Currents 1 and 2 [36,37].

Model and Analysis of the Magnetic Field
In the electromagnetic actuator, the magnetic field in the air gap of the stator is the main factor for improving performance.The homogeneity of the magnetic field in the air gap along the actuator's axial direction decides the stability and linearity of the actuator's thrust force.The magnetic flux density in the air gap decides the thrust force of an actuator and the ride comfort of a vehicle.First, the model of the magnetic flux density in an air gap is built, and the characteristic of the magnetic field is theoretically analyzed.

Magnetic Flux Density Model of a Single PM Ring
To analyze the quasi-parallel magnetic field in the air gap between the PM ring arrays and the inner wall of the actuator, the mathematical model of the magnetic flux density of the magnet field is built in the y direction according to the Ampere molecular circulation hypothesis, as shown in Figure 4.When the positive direction of the magnet field is +y direction, the spatial magnetic field is equivalent to the magnetic field generated by Surface Currents 1 and 2 [36,37].A schematic diagram of the magnetic field of a single PM ring is shown in Figure 4.The description of the variables in the formulas is shown in Table 5.The magnetic flux density formula at point P(x, y) in the y direction can be expressed as follows [38][39][40]: The magnetic flux density By11 and By12 generated by the surface currents 1 and 2 can be expressed as (2) and ( 3) [41][42][43]: A schematic diagram of the magnetic field of a single PM ring is shown in Figure 4.The description of the variables in the formulas is shown in Table 5.The magnetic flux density formula at point P(x, y) in the y direction can be expressed as follows [38][39][40]: where The magnetic flux density B y11 and B y12 generated by the surface currents 1 and 2 can be expressed as (2) and (3) [41][42][43]: (2) where Magnetic flux density of small PM rings B y Magnetic flux density of the PM ring array

Magnetic Flux Density Model of the Magnet Assembly
To analyze the magnetic fields of two groups of PM ring arrays, the magnetic flux density of the PM ring array B y in the air gap can be expressed as (4) [43]:

Theoretical Analysis of Magnetic Flux Density
According to the magnetic flux density model of the actuator's magnet assembly, the magnetic flux density was calculated by (4) in the middle of the air gap and along the whole axial length of 90 mm.The calculated results of the magnetic flux density of each PM ring and the superimposed magnetic flux density of the six PM ring arrays are shown in Figure 5.The results of each PM ring are shown as a normal distribution curve, the superimposed magnetic flux density B y is shown as regular fluctuations, and its value varies between 51.73 and 788.49mT.

FEM Analysis Model
To compare with the magnetic field theory model, a 2-D FEM analysis on the magnetic field and magnetic flux density was conducted using a magnetic field analysis software program.A half 2-D FEM model was built with the same dimension parameters as in Table 4, since the structure of the actuator is axisymmetric.In the FEM model, the remanence Br of the permanent magnet was 1.38-1.44T, and the coercivity was −876 KA/m.

FEM Analysis on the Magnetic Field
Based on the cylinder coordinate system, the electromagnetic actuator is simulated to obtain the magnetic density and magnetic flux trend by FEM.The simulation adopts a cylindrical coordinate system.The axle of the actuator is z-axis, shown in Figure 6, and the same graphical scales of 0-100 mm are used.The result of the FEM analysis clarifies the characteristics of the quasi-parallel magnetic field for the actuator using PM ring arrays.Since the S pole of the large PM rings is parallel to the N pole of the small PM rings in the air gap of the electromagnetic actuator, there is less magnetic flux leakage.The magnetic flux is nearly parallel in the air gap and its distribution is relatively even in the air gap except for the two ends.This proves that when it moves in the air gap, the coil moves through, cutting the parallel magnetic flux.
Assumption: Each coil has a diameter of 1 mm and a current of 0.8 A is applied.After applying direct current to each coil, the eddy current loss in the simulation is approximately 3 W. Figure 6 shows the simulation structure of the electromagnetic actuator at 0.04 s.
Figure 6a shows the FEM simulation results for the magnetic flux density.From the legend on the left, the magnetic flux density ranges from 4.9979 × 10 −3 T to 2.500 T. From the magnetic flux density results in the radial direction, it can be seen that the magnetic flux density in the air gap between the large PM rings and the small PM rings is about 550 mT, and the magnetic flux density at both ends is about 670 mT.The magnetic flux density of the large and small heat-dissipated rings is about 1.9 T. The magnetic flux density of the heat-dissipated rings in the middle of the same-size PM rings in the axial direction is relatively large.
Figure 6b shows the FEM simulation results in the direction of magnetic flux.From the legend on the left, the magnetic flux density for the direction of magnetic flux ranges from 4.9979 × 10 −3 T to 2.500 T, which is the same as in Figure 6a.The results show that the magnetic flux between the large PM rings and the small PM rings is perpendicular to the axial direction of the electromagnetic actuator, and the directions of the multiple magnetic

FEM Analysis Model
To compare with the magnetic field theory model, a 2-D FEM analysis on the magnetic field and magnetic flux density was conducted using a magnetic field analysis software program.A half 2-D FEM model was built with the same dimension parameters as in Table 4, since the structure of the actuator is axisymmetric.In the FEM model, the remanence Br of the permanent magnet was 1.38-1.44T, and the coercivity was −876 KA/m.

FEM Analysis on the Magnetic Field
Based on the cylinder coordinate system, the electromagnetic actuator is simulated to obtain the magnetic density and magnetic flux trend by FEM.The simulation adopts a cylindrical coordinate system.The axle of the actuator is z-axis, shown in Figure 6, and the same graphical scales of 0-100 mm are used.The result of the FEM analysis clarifies the characteristics of the quasi-parallel magnetic field for the actuator using PM ring arrays.Since the S pole of the large PM rings is parallel to the N pole of the small PM rings in the air gap of the electromagnetic actuator, there is less magnetic flux leakage.The magnetic flux is nearly parallel in the air gap and its distribution is relatively even in the air gap except for the two ends.This proves that when it moves in the air gap, the coil moves through, cutting the parallel magnetic flux.
Assumption: Each coil has a diameter of 1 mm and a current of 0.8 A is applied.After applying direct current to each coil, the eddy current loss in the simulation is approximately 3 W. Figure 6 shows the simulation structure of the electromagnetic actuator at 0.04 s.
Figure 6a shows the FEM simulation results for the magnetic flux density.From the legend on the left, the magnetic flux density ranges from 4.9979 × 10 −3 T to 2.500 T. From the magnetic flux density results in the radial direction, it can be seen that the magnetic flux density in the air gap between the large PM rings and the small PM rings is about 550 mT, and the magnetic flux density at both ends is about 670 mT.The magnetic flux density of the large and small heat-dissipated rings is about 1.9 T. The magnetic flux density of the heat-dissipated rings in the middle of the same-size PM rings in the axial direction is relatively large.
Figure 6b shows the FEM simulation results in the direction of magnetic flux.From the legend on the left, the magnetic flux density for the direction of magnetic flux ranges from 4.9979 × 10 −3 T to 2.500 T, which is the same as in Figure 6a.The results show that the magnetic flux between the large PM rings and the small PM rings is perpendicular to the axial direction of the electromagnetic actuator, and the directions of the multiple magnetic fluxes are parallel to each other.There will be divergence at the edge of each PM ring.Therefore, a magnetic field will be generated at the heat-dissipated rings.In theoretical calculations, the edge effects of the PM rings are assumed to be 0 mT, and this is reflected in the FEM simulation.
Figure 6c shows the FEM simulation results of the magnetic flux trend.From the legend on the left, the magnetic flux distribution for the magnetic flux trend ranges from −2.214 × 10 −4 T to 2.024 × 10 −4 T. The results show that the parallel magnetic field between the large PM rings and the small PM rings is visible, and the magnetic flux in the large heat-dissipated rings and the small heat-dissipated rings is repulsive; in the simulation boundary, due to the inevitable phenomenon of magnetic leakage, the magnetic flux must exist outside the electromagnetic actuator.fluxes are parallel to each other.There will be divergence at the edge of each PM ring.Therefore, a magnetic field will be generated at the heat-dissipated rings.In theoretical calculations, the edge effects of the PM rings are assumed to be 0 mT, and this is reflected in the FEM simulation.
Figure 6c shows the FEM simulation results of the magnetic flux trend.From the legend on the left, the magnetic flux distribution for the magnetic flux trend ranges from −2.214 × 10 −4 T to 2.024 × 10 −4 T. The results show that the parallel magnetic field between the large PM rings and the small PM rings is visible, and the magnetic flux in the large heat-dissipated rings and the small heat-dissipated rings is repulsive; in the simulation boundary, due to the inevitable phenomenon of magnetic leakage, the magnetic flux must exist outside the electromagnetic actuator.

FEM Analysis of the Magnetic Flux Density
To compare with the theoretical results of the magnetic field model, the superimposed magnetic flux density of the magnet assembly was also analyzed using the 2-D FEM model.The FEM analysis results in Figure 7 clarify the characteristics of the quasi-parallel magnetic field for the magnet assembly using PM ring arrays.The positions of the axial air gap of the PM ring arrays were analyzed, which is the same as the theoretical calculation positions.
Curve By is the magnetic flux density of the PM ring arrays in the axial air gap.Curve By1 to curve By6 are the magnetic flux densities of six magnetic rings from the lower to upper cover.In an air gap magnetic field, the magnetic flux density is roughly in the range of 100 mT-650 mT.Compared to the theoretical result curves, these curves roughly exhibit an irregular sinusoidal shape.At the axial center of the permanent magnet, the magnetic flux density is the highest.At the axial center of the heat-dissipation ring, the magnetic flux density is the smallest.The theoretical calculation curve trend is consistent with the FEM simulation curve trend, with a slightly lower peak magnetic flux density at the middle magnetic rings.The irregular shape of the curve may be due to the consideration of eddy current losses, edge effects of magnets, and other factors in the FEM simulation.

FEM Analysis of the Magnetic Flux Density
To compare with the theoretical results of the magnetic field model, the superimposed magnetic flux density of the magnet assembly was also analyzed using the 2-D FEM model.The FEM analysis results in Figure 7 clarify the characteristics of the quasi-parallel magnetic field for the magnet assembly using PM ring arrays.The positions of the axial air gap of the PM ring arrays were analyzed, which is the same as the theoretical calculation positions.
Curve By is the magnetic flux density of the PM ring arrays in the axial air gap.Curve By1 to curve By6 are the magnetic flux densities of six magnetic rings from the lower to upper cover.In an air gap magnetic field, the magnetic flux density is roughly in the range of 100 mT-650 mT.Compared to the theoretical result curves, these curves roughly exhibit an irregular sinusoidal shape.At the axial center of the permanent magnet, the magnetic flux density is the highest.At the axial center of the heat-dissipation ring, the magnetic flux density is the smallest.The theoretical calculation curve trend is consistent with the FEM simulation curve trend, with a slightly lower peak magnetic flux density at the middle magnetic rings.The irregular shape of the curve may be due to the consideration of eddy current losses, edge effects of magnets, and other factors in the FEM simulation.

FEM Analysis of the Magnetic Flux Density
To compare with the theoretical results of the magnetic field model, the superimposed magnetic flux density of the magnet assembly was also analyzed using the 2-D FEM model.The FEM analysis results in Figure 7 clarify the characteristics of the quasi-parallel magnetic field for the magnet assembly using PM ring arrays.The positions of the axial air gap of the PM ring arrays were analyzed, which is the same as the theoretical calculation positions.
Curve By is the magnetic flux density of the PM ring arrays in the axial air gap.Curve By1 to curve By6 are the magnetic flux densities of six magnetic rings from the lower to upper cover.In an air gap magnetic field, the magnetic flux density is roughly in the range of 100 mT-650 mT.Compared to the theoretical result curves, these curves roughly exhibit an irregular sinusoidal shape.At the axial center of the permanent magnet, the magnetic flux density is the highest.At the axial center of the heat-dissipation ring, the magnetic flux density is the smallest.The theoretical calculation curve trend is consistent with the FEM simulation curve trend, with a slightly lower peak magnetic flux density at the middle magnetic rings.The irregular shape of the curve may be due to the consideration of eddy current losses, edge effects of magnets, and other factors in the FEM simulation.To verify the curve trend of theoretical results and finite element results, the PM arrangement in [44] is similar to that used in this paper.Reference [44] studied an electromagnetic actuator for an active suspension (as shown in Figure 3 of [44]), which uses a single-layer magnetic ring array.The arrangement of magnetic poles is that the outer side is N pole, the inner side is S pole, and each magnetic ring is composed of eight same pieces.The magnetic flux (as shown in Figure 3 of [44]) is emitted from the N pole of the PM rings, passes through the air gap, the outer and inner walls, and then returns to the S pole of the PM rings.This paper adopts a double-layer magnetic ring array, in which the N poles of the small magnetic rings are in face of the S poles of the large magnetic rings in the air gap.Each magnetic ring is also composed of eight same pieces.
In the air gap, compared with the magnet ring array (as shown in Figure 6 of [44]), adding a layer of magnetic ring array will increase the magnetic flux density of the air gap.The theoretical and finite element results of the magnetic ring array are compared and analyzed with the experimental results, as shown in Figure 8.The magnetic field trends of the three curves are consistent, and each curve owns six peaks at the same position.This can prove the correctness of the magnetic flux density trend of the theoretical model derived in this paper.Due to the same position of six peaks, the magnetic flux density will increase as the volume of the magnet rings increases.The experimental results of magnetic field density in [44] are approximately in the range of 50-200 mT, and the theoretical and simulation results of magnetic field intensity in this paper are approximately in the range of 100-650 mT.Therefore, the arrangement of the magnet ring array in this paper with a large magnetic flux density in the air gap is superior to that in [44].To verify the curve trend of theoretical results and finite element results, the PM arrangement in [44] is similar to that used in this paper.Reference [44] studied an electromagnetic actuator for an active suspension (as shown in Figure 3 of [44]), which uses a single-layer magnetic ring array.The arrangement of magnetic poles is that the outer side is N pole, the inner side is S pole, and each magnetic ring is composed of eight same pieces.The magnetic flux (as shown in Figure 3 of [44]) is emitted from the N pole of the PM rings, passes through the air gap, the outer and inner walls, and then returns to the S pole of the PM rings.This paper adopts a double-layer magnetic ring array, in which the N poles of the small magnetic rings are in face of the S poles of the large magnetic rings in the air gap.Each magnetic ring is also composed of eight same pieces.
In the air gap, compared with the magnet ring array (as shown in Figure 6 of [44]), adding a layer of magnetic ring array will increase the magnetic flux density of the air gap.The theoretical and finite element results of the magnetic ring array are compared and analyzed with the experimental results, as shown in Figure 8.The magnetic field trends of the three curves are consistent, and each curve owns six peaks at the same position.This can prove the correctness of the magnetic flux density trend of the theoretical model derived in this paper.Due to the same position of six peaks, the magnetic flux density will increase as the volume of the magnet rings increases.The experimental results of magnetic field density in [44] are approximately in the range of 50-200 mT, and the theoretical and simulation results of magnetic field intensity in this paper are approximately in the range of 100-650 mT.Therefore, the arrangement of the magnet ring array in this paper with a large magnetic flux density in the air gap is superior to that in [44].

FEM Analysis on Magnetic Force
Based on the magnetic field characteristics of the electromagnetic actuator, it is necessary to analyze the magnetic force characteristics of the electromagnetic actuator, as shown in Figure 9.The magnetic force characteristics are divided into the current-magnetic force characteristics and the displacement-magnetic force characteristics.
The current-magnetic force curve is shown in Figure 9a.The current range is −1.8A to 1.8 A, and the magnetic force range of the electromagnetic actuator is −219.8N to 220.9 N. The current-magnetic force characteristic curve passes through the origin [0,0].When the current is 1 A, the magnetic force is 122.2N. When the current is −1 A, the magnetic force is −122.0N. The slope of the current force curve is approximately 122.2.

FEM Analysis on Magnetic Force
Based on the magnetic field characteristics of the electromagnetic actuator, it is necessary to analyze the magnetic force characteristics of the electromagnetic actuator, as shown in Figure 9.The magnetic force characteristics are divided into the current-magnetic force characteristics and the displacement-magnetic force characteristics.
The current-magnetic force curve is shown in Figure 9a.The current range is −1.8A to 1.8 A, and the magnetic force range of the electromagnetic actuator is −219.8N to 220.9 N. The current-magnetic force characteristic curve passes through the origin [0,0].When the current is 1 A, the magnetic force is 122.2N. When the current is −1 A, the magnetic force is −122.0N. The slope of the current force curve is approximately 122.2.
The displacement-magnetic force curve is shown in Figure 9b, which is the displacementmagnetic force curve when the current is 0.5 A, 1.0 A, and 1.5 A, respectively.With increasing displacement, the magnetic attenuation is 0.187 N/mm.When the current is 1.0 A and 1.5 A, the magnetic attenuation is 0.376 N/mm and 0.618 N/mm, respectively.
The displacement-magnetic force curve is shown in Figure 9b, which is the displacement-magnetic force curve when the current is 0.5 A, 1.0 A, and 1.5 A, respectively.With increasing displacement, the magnetic attenuation is 0.187 N/mm.When the current is 1.0 A and 1.5 A, the magnetic attenuation is 0.376 N/mm and 0.618 N/mm, respectively.

Vehicle Simulation Based on the Fuzzy PID Algorithm
To obtain the dynamic characteristics of the electromagnetic suspension, a seven-DOF dynamics model of the whole vehicle was established, and then the fuzzy PID simulation analysis was carried out.

Vehicle Dynamics Model
The pitch and roll characteristics of the whole vehicle are analyzed through the seven DOF dynamics model.The seven-DOF dynamics model based on fuzzy PID control is shown in Figure 10, and its parameters and descriptions are shown in Table 6.Based on the dynamics model of the passive suspension of the whole vehicle, electromagnetic suspension actuators are introduced at the four shock absorbers, and the magnetic forces FiA, FiB, FiC, and FiD are generated.

Vehicle Simulation Based on the Fuzzy PID Algorithm
To obtain the dynamic characteristics of the electromagnetic suspension, a seven-DOF dynamics model of the whole vehicle was established, and then the fuzzy PID simulation analysis was carried out.

Vehicle Dynamics Model
The pitch and roll characteristics of the whole vehicle are analyzed through the seven DOF dynamics model.The seven-DOF dynamics model based on fuzzy PID control is shown in Figure 10, and its parameters and descriptions are shown in Table 6.Based on the dynamics model of the passive suspension of the whole vehicle, electromagnetic suspension actuators are introduced at the four shock absorbers, and the magnetic forces F iA , F iB , F iC , and F iD are generated.According to the dynamics model of the whole vehicle with seven DOF, the dynamics model of the whole suspension system is established using Newton's law.The seven DOF are centroid acceleration, pitch angle acceleration, roll angle acceleration, and vertical displacements of four wheels.Then, the vertical displacements of the four wheels can be expressed by Equations ( 7)- (10), and the centroid acceleration, pitch angle acceleration, and roll angle acceleration can be expressed by Equations ( 11)-( 13), respectively.The seven-DOF dynamics model is as follows, and the physical quantities in Equations ( 7)-( 13) can be seen in Figure 10.x wB = k tB x gB − x wB + k sB (x bB − x wB ) + c sB .
x wC = k tC x gC − x wC + k sC (x bC − x wC ) + c sC .
x wD = k tD x gD − x wD + k sD (x bD − x wD ) + c sD .
x bA + k sA (x wA − x bA ) + c sB .
x bB + k sB (x wB − x bB ) + c sC .
x bC +k sC x wC − x bC )+c sD . where Vehicle parameters are shown in Table 7.

Fuzzy PID Control Algorithm
The fuzzy PID control block diagram is shown in Figure 11.The external disturbance of the vehicle suspension system is road excitation.The centroid acceleration is set to 0 mm/s 2 , and the input of the fuzzy controller is the deviation e and the variation of the deviation ec.Based on fuzzy control rules, PID parameters are adjusted online to meet different requirements for PID parameters with different deviations e and incremental deviations ec.The simulation employs a variable-step solver with a maximum sampling step of 0.01 s, utilizing the Ode45 solver algorithm.The simulation is set to run for a duration of 5 s.The actives of the simulation are centroid acceleration, acceleration of the roll angle, acceleration of the pitch angle, dynamic deflections for four wheels, and dynamic loads for four wheels.

Fuzzy PID Control Algorithm
The fuzzy PID control block diagram is shown in Figure 11.The external disturbance of the vehicle suspension system is road excitation.The centroid acceleration is set to 0 mm/s 2 , and the input of the fuzzy controller is the deviation e and the variation of the deviation ec.Based on fuzzy control rules, PID parameters are adjusted online to meet different requirements for PID parameters with different deviations e and incremental deviations ec.The simulation employs a variable-step solver with a maximum sampling step of 0.01 s, utilizing the Ode45 solver algorithm.The simulation is set to run for a duration of 5 s.The actives of the simulation are centroid acceleration, acceleration of the roll angle, acceleration of the pitch angle, dynamic deflections for four wheels, and dynamic loads for four wheels.

Vehicle Simulation Results Based on Fuzzy PID under C-Grade Surface
Based on Matlab/Simulink simulation software, the fuzzy PID control simulation analysis of electromagnetic suspension is carried out and compared with the passive suspension.Three evaluation indicators of vehicle ride comfort are analyzed, namely body acceleration, roll angle, and pitch angle.Take the C-grade road surface (as shown in Figure 12) as the input road surface and set the vehicle speed to 20 km/h.The fuzzy PID simulation parameters are adjusted and determined by the parameter setting method.Therefore, the fuzzy PID control law is shown in Table 8.The PID control parameters are kp = 1, ki = 20, and kd = 0.01.

Vehicle Simulation Results Based on Fuzzy PID under C-Grade Surface
Based on Matlab/Simulink simulation software, the fuzzy PID control simulation analysis of electromagnetic suspension is carried out and compared with the passive suspension.Three evaluation indicators of vehicle ride comfort are analyzed, namely body acceleration, roll angle, and pitch angle.Take the C-grade road surface (as shown in Figure 12) as the input road surface and set the vehicle speed to 20 km/h.The fuzzy PID simulation parameters are adjusted and determined by the parameter setting method.Therefore, the fuzzy PID control law is shown in Table 8.The PID control parameters are k p = 1, k i = 20, and k d = 0.01.

Fuzzy PID Control Algorithm
The fuzzy PID control block diagram is shown in Figure 11.The external disturbance of the vehicle suspension system is road excitation.The centroid acceleration is set to 0 mm/s 2 , and the input of the fuzzy controller is the deviation e and the variation of the deviation ec.Based on fuzzy control rules, PID parameters are adjusted online to meet different requirements for PID parameters with different deviations e and incremental deviations ec.The simulation employs a variable-step solver with a maximum sampling step of 0.01 s, utilizing the Ode45 solver algorithm.The simulation is set to run for a duration of 5 s.The actives of the simulation are centroid acceleration, acceleration of the roll angle, acceleration of the pitch angle, dynamic deflections for four wheels, and dynamic loads for four wheels.

Vehicle Simulation Results Based on Fuzzy PID under C-Grade Surface
Based on Matlab/Simulink simulation software, the fuzzy PID control simulation analysis of electromagnetic suspension is carried out and compared with the passive suspension.Three evaluation indicators of vehicle ride comfort are analyzed, namely body acceleration, roll angle, and pitch angle.Take the C-grade road surface (as shown in Figure 12) as the input road surface and set the vehicle speed to 20 km/h.The fuzzy PID simulation parameters are adjusted and determined by the parameter setting method.Therefore, the fuzzy PID control law is shown in Table 8.The PID control parameters are kp = 1, ki = 20, and kd = 0.01.The simulation results of centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle based on the fuzzy PID controller are shown in Figures 13-15, which show the simulation comparison results of three main evaluation indicators.The trend of active suspension and passive suspension of the three curves is the same.The root mean square (RMS) values of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle are 59.08%, 3.06%, and 3.54% higher than those of passive suspension, respectively.The tracking response of these three curves is fast.
In the simulation results of vehicle body acceleration in Figure 13, there is a slight delay, with a delay time of 0.05 s.The acceleration of roll angle curve is shown in Figure 14.At 3 s, the peak value of the electromagnetic suspension is 0.012 rad/s 2 less than that of the passive suspension.The acceleration of the pitch angle curve is shown in Figure 15.At 1.96 s, the maximum value of the electromagnetic suspension is 0.153 rad/s 2 less than that of the passive suspension.The trend of active suspension and passive suspension of the three curves is the same.The root mean square (RMS) values of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle are 59.08%, 3.06%, and 3.54% higher than those of passive suspension, respectively.The tracking response of these three curves is fast.
In the simulation results of vehicle body acceleration in Figure 13, there is a slight delay, with a delay time of 0.05 s.The acceleration of roll angle curve is shown in Figure 14.At 3 s, the peak value of the electromagnetic suspension is 0.012 rad/s 2 less than that of the passive suspension.The acceleration of the pitch angle curve is shown in Figure 15.At 1.96 s, the maximum value of the electromagnetic suspension is 0.153 rad/s 2 less than that of the passive suspension.The simulation results of the dynamic deflections for the four wheels are shown in Figure 16.The simulation results of the suspension's dynamic deflection for wheels A, B, C, and D are analyzed.The simulation results show that the maximum dynamic deflections of the electromagnetic active suspension are slightly lower than those of the passive suspension, and the controller has good tracking performance.The suspension dynamic deflections of wheels A, B, C, and D are 13.20%, 9.84%, 11.98%, and 16.42%, respectively, which are lower than the corresponding root mean square (RMS) value of the passive suspension.The simulation results of the dynamic deflections for the four wheels are shown in Figure 16.The simulation results of the suspension's dynamic deflection for wheels A, B, C, and D are analyzed.The simulation results show that the maximum dynamic deflections of the electromagnetic active suspension are slightly lower than those of the passive suspension, and the controller has good tracking performance.The suspension dynamic deflections of wheels A, B, C, and D are 13.20%, 9.84%, 11.98%, and 16.42%, respectively, which are lower than the corresponding root mean square (RMS) value of the passive suspension.The simulation results of the dynamic load for four wheels are shown in Figure 17.The dynamic load is closely related to the tire operation life and the impact of road damage.The simulation results show that the dynamic load of the four wheels of the electromagnetic suspension is less than the dynamic load of the passive suspension, which can reduce the impact of the tire on the road and the wear of the tire itself.The dynamic loads of the tires for the electromagnetic suspension show partial attenuation near each peak.The RMS values of the dynamic loads of wheels A, B, C, and D are 2.87%, 1.42%, 2.09%, The simulation results of the dynamic load for four wheels are shown in Figure 17.The dynamic load is closely related to the tire operation life and the impact of road damage.The simulation results show that the dynamic load of the four wheels of the electromagnetic suspension is less than the dynamic load of the passive suspension, which can reduce the impact of the tire on the road and the wear of the tire itself.The dynamic loads of the tires for the electromagnetic suspension show partial attenuation near each peak.The RMS values of the dynamic loads of wheels A, B, C, and D are 2.87%, 1.42%, 2.09%, and 5.04% lower than those of the passive suspension, respectively.

Vehicle Simulation Results Based on Fuzzy PID under a Deceleration Strip Surface
Based on the simulation on a C-level road surface in Section 5.3, a road excitation with a deceleration strip is introduced as shown in Figure 18.The deceleration strip adopts a trapezoidal shape with an upper edge of 200 mm, a lower edge of 400 mm, and a height of 30 mm.A fuzzy PID control simulation analysis was conducted on the electromagnetic suspension and compared with the passive suspension with a vehicle speed of 20 km/h.The fuzzy PID simulation parameters are adjusted and determined by setting the parameters.The fuzzy PID control law is shown in Table 9.The PID control parameters are kp = 0.8, ki = 50, and kd = 0.001.This simulation analysis can provide a detailed analysis of the pitch and roll characteristics to showcase the advantages of the electromagnetic suspension.

Vehicle Simulation Results Based on Fuzzy PID under a Deceleration Strip Surface
Based on the simulation on a C-level road surface in Section 5.3, a road excitation with a deceleration strip is introduced as shown in Figure 18.The deceleration strip adopts a trapezoidal shape with an upper edge of 200 mm, a lower edge of 400 mm, and a height of 30 mm.A fuzzy PID control simulation analysis was conducted on the electromagnetic suspension and compared with the passive suspension with a vehicle speed of 20 km/h.The fuzzy PID simulation parameters are adjusted and determined by setting the parameters.The fuzzy PID control law is shown in Table 9.The PID control parameters are k p = 0.8, k i = 50, and k d = 0.001.This simulation analysis can provide a detailed analysis of the pitch and roll characteristics to showcase the advantages of the electromagnetic suspension.The centroid acceleration result shows that the tracking performance of acceleration is good without any delay phenomena.At 2.13 s, the peak value of centroid acceleration using electromagnetic suspension is 0.557 m/s 2 smaller than when using passive suspension.In the simulation results of the acceleration of the roll angle, there is a slight delay, with a delay time of 0.04 s.The peak values of the accelerations of the roll angle and the pitch angle decreased by 0.319 rad/s 2 and 0.015 rad/s 2 , respectively.rs 2023, 12, x FOR PEER REVIEW The centroid acceleration result shows that the tracking pe is good without any delay phenomena.At 2.13 s, the peak valu using electromagnetic suspension is 0.557 m/s 2 smaller than w sion.In the simulation results of the acceleration of the roll an with a delay time of 0.04 s.The peak values of the acceleration pitch angle decreased by 0.319 rad/s 2 and 0.015 rad/s 2 , respectiv The dynamic deflections of the four wheels of the vehicle a RMS values of the dynamic deflection of wheels A, B, C, and 21.40%, 25.01%, and 24.24%, respectively.When the rear wheels the shock absorbers generate a downward dynamic deflection the dynamic deflections of the four dampers decrease by 0.00 and 0.0026 m, respectively.The fluctuation of the rear wheel is to the fluctuation of the front wheel.The dynamic The dynamic deflections of the four wheels of the vehicle are shown in Figure 22.The RMS values of the dynamic deflection of wheels A, B, C, and D decreased by 21.26%, 21.40%, 25.01%, and 24.24%, respectively.When the rear wheels pass through the reducer, the shock absorbers generate a downward dynamic deflection.After fuzzy PID control, the dynamic deflections of the four dampers decrease by 0.0057 m, 0.0047 m, 0.0042 m, and 0.0026 m, respectively.The fluctuation of the rear wheel is delayed by 0.02 s relative to the fluctuation of the front wheel.The dynamic deflections of the four wheels of the passive suspension are −0.0148m, −0.0155 m, −0.0127 m, and −0.0113 m, and it can be seen that the difference between the maximum and minimum is 0.0042 m.The dynamic deflections of wheels A, B, C, and D of the electromagnetic suspension are −0.0091m, −0.0108 m, −0.0085 m, and −0.0087 m, and it can be seen that the difference between the maximum and minimum is 0.0023 m.These results indicate that electromagnetic suspension can be effectively controlled on a deceleration strip or when braking, which makes up for the poor stability of the MacPherson suspension.deflections of the four wheels of the passive suspension are −0.0148m, −0.0155 m, −0.0127 m, and −0.0113 m, and it can be seen that the difference between the maximum and minimum is 0.0042 m.The dynamic deflections of wheels A, B, C, and D of the electromagnetic suspension are −0.0091m, −0.0108 m, −0.0085 m, and −0.0087 m, and it can be seen that the difference between the maximum and minimum is 0.0023 m.These results indicate that electromagnetic suspension can be effectively controlled on a deceleration strip or when braking, which makes up for the poor stability of the MacPherson suspension.
The simulation results of the dynamic load for the four wheels are shown in Figure 23.The RMS values of the dynamic loads of wheels A, B, C, and D are 11.43%, 11.91%, 14.91%, and 15.03% lower than those of the passive suspension, respectively.When passing through the deceleration strip, the dynamic loads of the front wheel tires are larger than those of the rear wheel tires.The results show that when the vehicle passes through the deceleration strip, the dynamic load of the tires can be reduced and the excessive friction of the tires can be effectively prevented via the fuzzy PID control algorithm.The simulation results of the dynamic load for the four wheels are shown in Figure 23.The RMS values of the dynamic loads of wheels A, B, C, and D are 11.43%, 11.91%, 14.91%, and 15.03% lower than those of the passive suspension, respectively.When passing through the deceleration strip, the dynamic loads of the front wheel tires are larger than those of the rear wheel tires.The results show that when the vehicle passes through the deceleration strip, the dynamic load of the tires can be reduced and the excessive friction of the tires can be effectively prevented via the fuzzy PID control algorithm.

Conclusions
This paper proposes an electromagnetic suspension that can ensure the safety of a vehicle and improve the ride comfort and stability of the MacPherson suspension.This suspension does not change the safety of the MacPherson suspension with embedding an electromagnetic actuator.Once the electromagnetic actuator fails, the safety of the electromagnetic suspension can be guaranteed.
This paper mainly studies its magnetic characteristics and fuzzy PID control.First, an electromagnetic suspension is proposed, and its theoretical magnetic field model is established and analyzed.Second, on the basis of the electromagnetic suspension system, a FEM model is established, and the magnetic field, magnetic flux density, and magnetic force are simulated and analyzed.Finally, on the basis of the fuzzy PID control algorithm, a vehicle simulation analysis of the electromagnetic suspension is carried out.The results show that the magnetic field trends of the theoretical results and the finite element results are consistent via a comparison.The direction and the trend of the magnetic flux density in the FEM simulation show that the magnetic fluxes are parallel to each other in the air gap of the magnet assembly.Comparing the theoretical results, finite element results, and experimental results in [44], the electromagnetic suspension proposed in this paper owns a larger magnetic flux density in the air gap.At the same time, it can also prove the correctness of the trend of the theoretical results and the feasibility of the proposed electromagnetic actuator.Based on the fuzzy PID control algorithm, the vehicle dynamics of the electromagnetic suspension are simulated.The simulation results show that compared to a passive suspension system, RMS values of the acceleration, the acceleration of roll angle, and the acceleration of pitch angle for the vehicle body are reduced by 59.08%, 3.06%, and 3.54%, respectively, and the dynamic deflection of the four suspension systems and the dynamic load of the four wheels are reduced.Therefore, this control algorithm is suitable for an electromagnetic suspension system.

Conclusions
This paper proposes an electromagnetic suspension that can ensure the safety of a vehicle and improve the ride comfort and stability of the MacPherson suspension.This suspension does not change the safety of the MacPherson suspension with embedding an electromagnetic actuator.Once the electromagnetic actuator fails, the safety of the electromagnetic suspension can be guaranteed.
This paper mainly studies its magnetic characteristics and fuzzy PID control.First, an electromagnetic suspension is proposed, and its theoretical magnetic field model is established and analyzed.Second, on the basis of the electromagnetic suspension system, a FEM model is established, and the magnetic field, magnetic flux density, and magnetic force are simulated and analyzed.Finally, on the basis of the fuzzy PID control algorithm, a vehicle simulation analysis of the electromagnetic suspension is carried out.The results show that the magnetic field trends of the theoretical results and the finite element results are consistent via a comparison.The direction and the trend of the magnetic flux density in the FEM simulation show that the magnetic fluxes are parallel to each other in the air gap of the magnet assembly.Comparing the theoretical results, finite element results, and experimental results in [44], the electromagnetic suspension proposed in this paper owns a larger magnetic flux density in the air gap.At the same time, it can also prove the correctness of the trend of the theoretical results and the feasibility of the proposed electromagnetic actuator.Based on the fuzzy PID control algorithm, the vehicle dynamics of the electromagnetic suspension are simulated.The simulation results show that compared to a passive suspension system, RMS values of the acceleration, the acceleration of roll angle, and the acceleration of pitch angle for the vehicle body are reduced by 59.08%, 3.06%, and 3.54%, respectively, and the dynamic deflection of the four suspension systems and the dynamic load of the four wheels are reduced.Therefore, this control algorithm is suitable for an electromagnetic suspension system.

Figure 2 .
Figure 2. Structure of the magnet assembly.

Figure 2 .
Figure 2. Structure of the magnet assembly.

Figure 2 .
Figure 2. Structure of the magnet assembly.

Figure 3 .
Figure 3. Dimensional drawing.(a) Sectional view in the radial direction.(b) Sectional view in the axial direction.

Figure 3 .
Figure 3. Dimensional drawing.(a) Sectional view in the radial direction.(b) Sectional view in the axial direction.

Figure 4 .
Figure 4. Schematic diagram of the magnetic field of a single PM ring.

Figure 4 .
Figure 4. Schematic diagram of the magnetic field of a single PM ring.

Figure 5 .
Figure 5. Theoretical results of magnetic flux density.

Figure 5 .
Figure 5. Theoretical results of magnetic flux density.

Figure 6 .
Figure 6.FEM analysis of magnetic field for half of the actuator.(a) The magnetic flux density.(b) The direction of magnetic flux.(c) The trend of magnetic flux.

Figure 6 .
Figure 6.FEM analysis of magnetic field for half of the actuator.(a) The magnetic flux density.(b) The direction of magnetic flux.(c) The trend of magnetic flux.

Figure 6 .
Figure 6.FEM analysis of magnetic field for half of the actuator.(a) The magnetic flux density.(b) The direction of magnetic flux.(c) The trend of magnetic flux.

Figure 7 .
Figure 7. Simulation results of magnetic flux density.

Figure 7 .
Figure 7. Simulation results of magnetic flux density.

Figure 8 .
Figure 8. Three results of magnetic flux density.

Figure 8 .
Figure 8. Three results of magnetic flux density.

Figure 9 .
Figure 9. Simulation results of magnetic force characteristics.(a) The current-magnetic force curve.(b) The displacement-magnetic force curve.

Figure 9 .
Figure 9. Simulation results of magnetic force characteristics.(a) The current-magnetic force curve.(b) The displacement-magnetic force curve.

Figure 10 .
Figure 10.Seven-DOF dynamics model based on fuzzy PID control.
= k tA x gA − x wA + k sA (x bA − x wA ) + c sA .

Figure 12 .
Figure 12.Road excitation on a C-grade road surface for four wheels.The simulation results of centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle based on the fuzzy PID controller are shown in Figures 13-15, which show the simulation comparison results of three main evaluation indicators.

Figure 12 .
Figure 12.Road excitation on a C-grade road surface for four wheels.The simulation results of centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle based on the fuzzy PID controller are shown in Figures 13-15, which show the simulation comparison results of three main evaluation indicators.

Figure 12 .
Figure 12.Road excitation on a C-grade road surface for four wheels.

Figure 14 .
Figure 14.Acceleration of the roll angle on a C-grade road surface.

Figure 14 .
Figure 14.Acceleration of the roll angle on a C-grade road surface.

Figure 15 .
Figure 15.Acceleration of the pitch angle on a C-grade road surface.

Figure 15 .
Figure 15.Acceleration of the pitch angle on a C-grade road surface.

Figure 16 .
Figure 16.Dynamic deflections on a C-grade road surface for four wheels.

Figure 16 .
Figure 16.Dynamic deflections on a C-grade road surface for four wheels.

Actuators 2023 , 26 Figure 17 .
Figure 17.Dynamic loads on a C-grade road surface for four wheels.

Figure 17 .
Figure 17.Dynamic loads on a C-grade road surface for four wheels.

Figure 18 .
Figure 18.Road excitation with the C deceleration strip for four wheels.The simulation results of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle based on a fuzzy PID controller are shown in Figures 19-21, which show the simulation comparison results of the three main evaluation indicators.The RMS values of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle of active suspension are 33.34%,4.82%, and 4.10% higher than those of passive suspension, respectively.

Figure 18 .Figure 18 .
Figure 18.Road excitation with the C deceleration strip for four wheels.The simulation results of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle based on a fuzzy PID controller are shown in Figures 19-21, which show the simulation comparison results of the three main evaluation indicators.The RMS values of the centroid acceleration, acceleration of the roll angle, and acceleration of the pitch angle of active suspension are 33.34%,4.82%, and 4.10% higher than those of passive suspension, respectively.

Figure 20 .Figure 20 .
Figure 20.Acceleration of the roll angle on a deceleration strip.

Figure 21 .
Figure 21.Acceleration of the pitch angle on a deceleration strip.

Figure 21 .
Figure 21.Acceleration of the pitch angle on a deceleration strip.

Figure 22 .
Figure 22.Dynamic deflections on a deceleration strip for four wheels.

Figure 23 .
Figure 23.Dynamic loads on a deceleration strip for four wheels.

Figure 23 .
Figure 23.Dynamic loads on a deceleration strip for four wheels.

Table 1 .
Factors and levels in the orthogonal experiment method.

Table 3 .
Analysis by the orthogonal experiment method.

Table 4 .
Size parameters of the electromagnetic actuator.

Table 4 .
Size parameters of the electromagnetic actuator.
Q(a, b)Horizontal and vertical coordinates of any current source P(x, y)Horizontal and vertical coordinates of any point in a plane jRing number from left to right, j = 1-6 θ 1 , θ 2 Angle between the upper and lower points of surface current 1 θ Angle between the line of point Q and point P and the horizontal line r Distance between point Q and point P B yj1Magnetic flux density of surface current 1 B yj2Magnetic flux density of surface current 2 B y1Magnetic flux density of large PM rings B y2

Table 6 .
Variable table for seven-degree-of-freedom vehicle suspension.

Table 6 .
Variable table for seven-degree-of-freedom vehicle suspension.

Table 8 .
Control law on C-grade road surface for fuzzy PID algorithm.

Table 8 .
Control law on C-grade road surface for fuzzy PID algorithm.

Table 9 .
Control law with a deceleration strip for fuzzy PID algorithm.

Table 9 .
Control law with a deceleration strip for fuzzy PID algorithm.