A Magnetoelectric Hybrid Hydraulic Damper for the Mining Robot Suspension System

In order to improve the damping and controllability of the mining robot suspension system, a new magnetoelectric hybrid suspension hydraulic damper, which is a semiactive suspension damper, is proposed based on the traditional hydraulic damper by introducing the magnetic-electric hybrid suspension structure.-e structure and working principle of the damper are introduced, respectively, and the mathematical models of the equivalent stiffness and equivalent damping of the system are calculated by the magnetic circuit method and the oil circuit method, while AMESim/Simulink cosimulation is carried out. In order to test the damping performance, a prototype of the magnetoelectric hybrid suspension hydraulic damper was fabricated. -e results show that the vibration displacement amplitude can be reduced by 20% and the vibration acceleration amplitude can be reduced by 10% by adjusting the stiffness and damping of the system due to the magnetoelectric hybrid suspension structure. Moreover, the experimental results are consistent with the simulation results, which verify the effectiveness and superiority of this type of damper.


Introduction
As a new type of intelligent mining equipment, mining robot has become a research focus in recent years due to its characteristics of automation and high efficiency [1,2]. e working environment in underground coal mines is usually very harsh, and the floor is uneven, which has a great impact on the normal operation of the robot. In particular, it is easy to damage the precision instruments and equipment during the visual data collection, resulting in unnecessary losses [3]. In order to reduce the impact and vibration of mining robot and ensure the stable and efficient work efficiency, it is urgent to research and develop suitable suspension systems for mining robots. e main function of the suspension is to cushion road impact and suppress body vibrations caused by road surface unevenness. Depending on the level of controllability of damping, suspension systems are classified as passive, active, and semiactive [4]. Passive suspension control solutions are low cost and consist of springs and dampers with fixed properties, but these specifications are not adjustable. Active suspension is a system that adds an active controller, which actively generates the corresponding damping force according to different external disturbances as a way to adjust the movement of the body. e practical application of active suspension in vehicles is limited due to its higher power consumption and lower reliability [5]. Semiactive suspension is widely used, which does not require high power semiactive suspension, while it adjusts the damping characteristics by means of low power signals. On the other hand, the semiactive suspension can work as a passive suspension when the controller fails [6].
At present, the most studied mode of vibration reduction is semiactive suspension. Semiactive suspension mainly uses air suspension, pneumatic suspension, and other methods to reduce vibration by improving damping adjustment structure [7,8]. Gokul and Malar [9] studied the static vertical stiffness characteristics of the air spring in the air suspension system. Sang et al. [10] studied the characteristics of doublechamber hydrocarbon suspension. Nieto et al. [11] designed a pneumatic suspension capable of adapting to incoming vibrations, performing a switching control strategy between two different configurations. Maciejewski et al. [12] proposed the use of pneumatic muscles to achieve active vibration control of horizontal seat suspensions. Pintado et al. [13] proposed a fusion index measured by a thermodynamic model of a pneumatic vibration isolator to obtain a simple model of the airflow of the limiter. Zhu et al. [14] studied a new type of semiactive dual-chamber hydropneumatic inerter-based suspension. Zhao et al. [15] proposed an integrated sliding mode control-two-point wheelbase preview strategy for the semiactive air suspension system with a gasfilled adjustable shock absorber. Although air suspension and pneumatic suspension have been applied to some extent, the leakage problem, Jeyasenthil and Choi [16] and Huang et al. [17], has been restricting their development. e electromagnetic suspension system stands out as a new concept of the semiactive or active suspension system, which focuses on stiffness regulation and studies how to reduce the vibration amplitude of the system, which is very consistent with the operating conditions of mining robot [18]. Ding et al. [19] analyzed the three-dimensional numerical calculation and optimization of electromagnetic suspension stiffness in maglev trains. Emre and Pinar [20] proposed a systematic research approach based on improving the lateral sway stiffness of suspensions. Babak et al. [21] studied a hybrid electromagnetic damper to solve a series of problems. Maciejewski et al. [22] studied the modeling process and control strategy of a semiactive seat suspension system with magnetorheological dampers (MR). Mohammed and Alktranee [23] used three strategies to study the behavior of suspension systems. Jamadar et al. [24] proposed a mathematical model for equivalent damping based on magnetorheological dampers (MR). Ding et al. [25] designed a new hybrid electromagnetic actuator and its control system. Gysen et al. [26] deduced several specifications for the design of electromagnetic suspension systems. Although electromagnetic suspension has been deeply studied by scholars for a long time, the application of magnetoelectric hybrid suspension structure is rare. is structure interacts with permanent magnet through electromagnet, which has the advantages of large force and quick response [27,28].
is study designs a magnetoelectric hybrid suspension structure based on the traditional hydraulic damper and proposes a new type of the hydraulic damper for magnetoelectric hybrid suspension.

Principle of Magnetoelectric Hybrid Suspension.
e principle of magnetoelectric hybrid suspension is that the electromagnet and the permanent magnet are formed into a group of corresponding components. Because of the magnetic pole, there is magnetic force between the two poles.
ere are many ways to realize electromagnetic force, among which the electromagnetic-permanent magnetic suspension structure is the simplest and most superior. Figure 1 shows the distribution mode. e main advantages of electromagnetic-permanent magnet suspension structure are as follows: (1) No contact: this is the main advantage of maglev, which can bring frictionless resistance, no mechanical wear, low energy loss, low noise, low maintenance cost, and a series of other advantages with a great application value (2) Easy active control: electromagnetic force can be controlled only by changing the current, which can be applied in various systems that need vibration reduction and support hardness can be changed (3) Energy saving: the coercive force, residual magnetic induction, and maximum magnetic energy product of permanent magnets are relatively large, which can provide a constant magnetic field with a relatively large intensity in a small volume;

Structure and Principle.
In view of the above advantages and for the purpose of changing the system stiffness, this study proposes a magnetoelectric hybrid suspension hydraulic damper. Figure 2 shows the structure model. e main components of the magnetoelectric hybrid suspension hydraulic damper are piston rod, working cylinder, electromagnet, permanent magnet, spring, and so on. e upper and lower supports are, respectively, connected with the car body and the tire. When the car body vibrates, the piston rod drives the electromagnet to move up and down, and the spring then compresses, stretches, and reciprocates. After being squeezed, the hydraulic oil continuously generates damping force through the valve opening of the damping valve to reduce the harmful vibration caused by the spring's own recovery. At the same time, because electromagnet and permanent magnet produce electromagnetic force relative to each other, the spring vibration can be enhanced or weakened. According to Lorenz's law, the electromagnetic force can be adjusted by changing the size and direction of the current in the electromagnet. In order to achieve the purpose of stable vibration reduction, when the spring vibration is larger, the electromagnetic force is reduced to reduce the spring vibration trend; when the spring vibrates less, the electromagnetic force is increased to increase the damping effect.

Electromagnetic.
e stiffness control of the magnetoelectric hybrid suspension hydraulic damper is realized by adjusting the electromagnetic force. e calculation process of the electromagnetic force is relatively complicated, which is usually calculated by the magnetic circuit method. In order to facilitate analysis and calculation, the following assumptions are made based on Figure 1: (1) e leakage flux of the coil, air gap, and permanent magnet is zero 2 Shock and Vibration (2) e magnetic resistance of the core is ignored, and the magnetic potential is considered to be evenly distributed on air gap and permanent magnets. Figure 3 shows the simplified magnetic circuit diagram. According to the magnetic flux continuity theorem, where S is the cross-sectional area of the top of the core, S m is the magnetic pole area of the permanent magnet, B δ is the air gap magnetic flux density, and B m is the magnetic flux density of the permanent magnet.
According to the Ampere loop theorem, High performance of the sintered NdFeB permanent magnetic material N48 is selected, with remanence of B r and residual coerced force of H c . e relationship between the flux density of the permanent magnet and the magnetic field strength and coercivity is as follows: where μ 0 is the vacuum permeability and μ r is the relative permeability of the permanent magnet.
Magnetic flux density at the air gap is where the second term of the molecule on the right side of the equation is another magnetomotive force H c d m generated after the magnetic circuit is connected to a permanent magnet.  At any instant, the magnitude of the electromagnetic force between the permanent magnet and the electromagnet is where x e is the equivalent air gap of the system, U m is the voltage through the electromagnetic coil, and R m is the resistance of the electromagnetic coil. When d m � 0, equation (5) is converted to Equations (5) and (6), respectively, represent the expressions of electromagnetic force in magnetoelectric hybrid suspension (permanent magnet) and pure electromagnetic suspension (ordinary conductor) models. By comparing the two equations, it can be seen that, unlike the pure electromagnetic suspension model, the electromagnetic force in the magnetoelectric hybrid suspension model is no longer proportional to the square of the current and inversely proportional to the square of the air gap due to the thickness of the permanent magnet, but shows a nonlinear relationship. When the structure parameters of the permanent magnet are determined, the electromagnetic force expression of the controlled parameter is only the value of the voltage through electromagnetic coil U m . Due to the residual coercivity of permanent magnets H c is more giant, under the same input voltage value, the magnetoelectric hybrid suspension system of electromagnetic force F than pure electromagnetic suspension of the electromagnetic force F′, the magnetoelectric mixed suspension structure is more energy efficient.

Damping Force.
e analysis of damping force needs to be calculated through the oil circuit diagram. Figure 4 shows the simplified oil circuit diagram of the magnetoelectric hybrid suspension hydraulic damper.
In Figure 4, D 1 is the diameter of the piston rod, D 2 is the diameter of the inner cylinder, d is the diameter of the throttle hole in the damping valve, n is the number of throttle holes, P 1 is the pressure in the inner cylinder, P 2 is the pressure in the outer cylinder, and v is the movement speed of the piston cylinder.
Based on the above simplified oil circuit diagram, the flow through the damping valve is where C is the flow coefficient of the throttle hole, whose value is related to the geometric shape of the damping hole, generally, 0.6-0.8 is taken, C z is the flow correction coefficient between the damping holes, A 1 is the area of the damping hole, A 1 � πd 2 /4, ΔP is the cylinder pressure difference inside and outside, and ΔP � P 1 − P 2 , ρ is the oil density. e volume of piston rod flowing into the external cylinder shall be equal to the flow through the damping valve, i.e., where A 2 is the cross-sectional area of the piston rod, A 2 � πD 2 1 /4. It can be obtained from equations (7) and (8): e damping force expression can be obtained from equation (9): where K is the spring stiffness, x 0 is the initial compression amount of the spring, x is the compression amount of the spring after the movement of the piston rod, A 3 is the crosssectional area of the inner cylinder, and A 3 � πD 2 2 /4.

Road Excitation.
e input of body vibration mainly comes from the unevenness of the road surface, i.e., the effect of road excitation on body vibration. According to whether the pavement excitation is continuous or not, it can be divided into discrete pavement excitation and random pavement excitation. Arched pavement, waveform pavement, and trapezoidal pavement belong to discrete pavement excitation. e rolled pavement belongs to the random type pavement excitation. e research object of this study is a mining robot, and the working environment is often accompanied by gangue, water, pits, and other pavements, and the environment has strong randomness, so the random pavement excitation is mainly studied.
Since the pavement excitation is random, it can be characterized by a spectral density function for better analysis and calculation. In ISO 8608 : 1995(E), the pavement power spectral density G q (n) is recommended to be expressed as follows [29]: where G q is the power spectral density of the pavement in the vertical direction in m 3 , n is the spatial frequency, n 1 � 0.01m −1 , n 0 is the corresponding reference spatial frequency, n 0 � 0.1m −1 , G q (n 0 ) is the pavement unevenness coefficient, and W is the frequency index, dimensionless number, and usually takes the value of 2. e function G q (n) is a statistical characteristic of the pavement, which is a definite value. In the analysis of the impact of road surface vibration of vehicles, only considering the pavement power spectrum density G q (n) is not comprehensive enough but also the role of vehicle speed should be considered. Let the time frequency be f, the vehicle speed be v, and the spatial frequency be n. e relationship is as follows: f � vn. (12) In this study, the time domain research method is chosen, so the frequency domain model of the pavement must be converted into a time domain model. From equation (12), the conversion relationship between time and space spectrum is obtained as When W � 2, equations (11) and (12) are substituted into equation (13), and it can be deduced that e road surface unevenness can be considered as a white noise w(t) generated by a filter, from which the random excitation of the road surface can be obtained as [30] _ Equation (15) is the time domain model of pavement excitation used in the vehicle suspension analysis, but the expression is not accurate enough when the pavement spectrum is in the low frequency band. In order to improve the accuracy and reliability of the model, a constant called the lower cutoff frequency f 0 is introduced here, and then, equation (15) becomes Internationally, the pavement unevenness is divided into 8 levels according to the degree of pavement excitation, as given in Table 1 [31].

System Vibration Mechanical Model.
In order to study the semiactive suspension system in a simple way, a magnetoelectric hybrid suspension system based on a 2-DOF system with a 1/4 travel mechanism was developed with integrity in mind. Figure 5 shows the vertical dynamics modeling of the semiactive suspension system, and Table 2 provides the meaning of the parameters.
According to Newton's second equation, the vibration differential equation of the 1/4 walking mechanism semiactive suspension with 2-DOF is established as shown in the following equation [32].
Assuming that the system state vector is output vector is and input vector is then the state equation is Combine equations (17) and (21) to get the following equation: where A represents the state matrix, B represents the input matrix, C represents the output matrix, and D represents the transfer matrix. Equation (22) is the state equation required by the control system.

Equivalent Stiffness and Equivalent Damping.
e various components of the semiactive suspension system are affected by each other. is effect is mostly nonlinear. Especially after the introduction of magnetoelectric hybrid suspension, the dynamic characteristics of the system will change. In order to better express the stiffness adjustment method of the new type of the magnetoelectric hybrid suspension hydraulic damper, it is necessary to analyze the overall vibration model of the semiactive suspension system and obtain the equivalent stiffness and equivalent damping expressions, respectively, according to the structure and working principle of the damper. e equivalent stiffness of the system is mainly related to the spring stiffness and the electromagnetic force. e equivalent stiffness of the system is mainly related to the spring stiffness and electromagnetic force. Based on equations (5) and (17), the expression of equivalent stiffness k e can be obtained: where spring stiffness k s and vacuum permeability μ 0 are usually taken as 1.0, and core cross-sectional area S, coil turns N, residual coercive force H c , coil resistance R m , and permanent magnet thickness d m are all fixed parameters. Tire displacement x u , fuselage displacement x s , and suspension air gap x m are constantly changing with system vibration and cannot be actively adjusted. e only controllable parameter is the input voltage value U m of the solenoid coil. It can be seen from equation (24) that the equivalent stiffness of the system has a nonlinear relationship with the voltage value, and the increase or decrease of the voltage value will increase or decrease the system stiffness. e equivalent damping of the system is mainly related to the damping force. e expressions of equivalent damping c e can be obtained by combining equations (10) and (17): where the area of the damping hole A 1 , the cross-sectional area of the piston rod A 2 , the cross-sectional area of the inner cylinder A 3 , the initial spring compression x 0 , the spring compression after the piston rod moves x 1 , the oil density ρ, the number of throttle holes n, and the movement speed of the piston rod v are not adjustable.

Establishment of the AMESim Model.
AMESim software can greatly simplify the simulation model on the basis of maintaining the original function, and it is widely used in the automobile suspension system [33,34]. Based on the vibration differential equations and operating principle of the magnetoelectric hybrid suspension and semiactive suspension, Figure 6 shows the AMESim model of the 2-degree-offreedom 1/4 vehicle semiactive suspension system. In this AMESim model, the mechanics module (green part) constitutes the main body of the semiactive suspension system, including springs, damping, tire mass, car body mass, displacement module, force module, and displacement converter module. e electromagnetic module (the purple part) constitutes a magnetic-electric hybrid suspension system, including an electromagnet, a permanent magnet, a power supply, and an electromagnetic force output module. e signal module (red part) constitutes a signal and control system, including road excitation, signal integration and control, and feedback control.
In this model, a finite bandwidth white noise signal is used as the road excitation to generate vibration through a semiactive suspension system consisting of spring-tire-damping, spring, and electromagnetic force-body. e electromagnetic force is provided by electromagnets and permanent magnets. e joint MATLAB/Simulink simulation module collects and analyzes the spring-load displacement (Bdisplacement) and unspring displacement (Wdisplacement) signals and then outputs current signal to the electromagnets to control the electromagnetic force and the body displacement changes accordingly.
is whole process forms a closed-loop feedback system. e main parameter information of the AMESim model of the magnetoelectric hybrid semiactive suspension is given in Table 3.   Road excitation displacement m e established system model includes signal input, acquisition, processing, output, and feedback, which can more accurately simulate the working conditions of the magnetoelectric hybrid suspension, semiactive suspension system, which is beneficial to data analysis and research.

Establishment of Simulink Control
Module. Similar to AMESim, the Simulink visual editor is widely used in the control field and can be used for joint simulation with AMESim with high efficiency and reliability. In the field of semiactive suspension control, the commonly used Simulink control strategies are fuzzy control, PID control, and optimal control.
In the loop control process, due to disturbance, the controlled parameter often deviates from the predetermined value, which will cause a certain system deviation between the measured value and the predetermined value. en, the automatic adjustment unit performs proportional (P), integral (I), and differential (D) operations on the deviation signal that can reduce the error between the measured value and the predetermined value; finally, it will output a standard signal to realize automatic control of the parameter.
Based on the vibration differential equation of the magnetoelectric hybrid semiactive suspension and the control strategy of PID, which are combined with the AMESim joint simulation model, the design of the Simulink control module is shown in Figure 7. In the Simulink controller model, the displacement signals (Bdisplacement and Wdisplacement) are input, and the current signal obtained by adjusting the parameters through the PID controller is used as the output to achieve the purpose of adjusting the electromagnetic force. e main design parameters of the PID controller are given in Table 4.

Simulation Results.
In the vibration simulation process of AMESim, the system vibration is controlled and adjusted jointly with the Simulink control module, while setting the simulation parameters in Table 3, so as to achieve the purpose of simulating the real vibration of the hydraulic damper of the magnetoelectric hybrid suspension.
Here, the three evaluation indexes with the most influential factors are selected: vehicle vertical acceleration, suspension disturbance, and dynamic tire load. Given that this study mainly experiments on the mechanism of the device, it is impossible to consider the comparison of the three performance indicators, so the vertical acceleration of the vehicle body and the dynamic deflection of the suspension are mainly considered, and the vibration displacement curve and acceleration curve of the damper are obtained by simulation. e magnetoelectric hybrid semiactive suspension and ordinary passive suspension are placed in the same simulation environment. e simulation time sets as 10 s, the sampling time sets as 0.01 s, the white noise amplitude sets as  Figure 8 shows the PSD information of the excitation signal. Figure 9 shows the simulation curve of dynamic characteristics of magnetoelectric hybrid suspension. As can be seen from Figure 9, the amplitude of the semiactive suspension vibration curve with the magnetoelectric hybrid suspension hydraulic damper is smaller than that of ordinary passive suspension. On the one hand, for the displacement curve, the amplitude at the first trough is reduced by 33.3%, and the amplitude at the first peak is reduced by 13.5%. e maximum valley difference of the semiactive suspension curve is 6.8 mm, which is 8% lower than the passive suspension curve. On the other hand, for the acceleration curve, the maximum acceleration with passive suspension is reduced from nearly 2.5 m/s 2 to 0.3 m/s 2 with semiactive suspension. In terms of the curve fluctuation, the semiactive suspension curve is flatter, that is, the suspension disturbance and vehicle body vertical acceleration are lower, the shock absorption effect and the ride comfort of the vehicle are better, and the damping effect is more obvious.
e results show that the simulation curve shows some advantages. e next step is to use the experimental device to verify the simulation results.

Establishment of Experimental Platform
According to Lorentz law, the electromagnetic force in the electromagnetic-permanent magnetic structure can be controlled by changing the voltage in the electromagnet. However, the most direct and easily measured parameters in the magnetoelectric hybrid suspension vibration system are not electromagnetic force and voltage, but the vibration displacement of the hydraulic damper. erefore, in this study, the vibration displacement of the hydraulic damper is taken as the collection object, and the relevant experimental platform is built. e experimental design includes hardware and software systems. It is expected to complete the verification of the theoretical model by analyzing the vibration displacement curve of the system.

System
Hardware. Data collection and processing consists of hardware and software, and Figure 10 shows the data acquisition hardware system. e main source of vibration in the vibration system is the shaker. Voltage signals are input to the shaker through the electric control box, and the shaker converts the voltage  Shock and Vibration 9 signal into the displacement signal to make the hydraulic damper vibrate. e laser displacement sensor is installed above the hydraulic damper, which is responsible for collecting the vibration displacement signal of the hydraulic damper and converting the displacement signal into the analog voltage signal, as shown in figure. en, the analog voltage signal is transmitted to the upper computer through the data acquisition card. e upper computer in this study is a personal computer. e data are analyzed and processed by the software program, and the processed voltage signal is transmitted to the driver, which controls the voltage applied to the electromagnet to change the electromagnetic force. After the electromagnetic force is changed, the spring is compressed or stretched by different forces, and the overall vibration absorption performance of the system is changed, resulting in the change of the displacement signal collected by the laser displacement sensor, forming a feedback system. Table 5 provides the main parameters of the system hardware selection.
After the experimental device is designed, the experimental platform can be built. According to the structure model and working principle of the designed magnetoelectric hybrid suspension hydraulic damper, the hydraulic damper can be placed vertically, and from bottom to top are the vibrator-base (acting as a connection between the exciter and the hydraulic damper)-hydraulic damper-vibrating table (i.e., bracket, used to fix and limit the amplitude of the system), the microlaser displacement sensor is fixed on the vibrating table and connected to the data acquisition card. Figure 11 shows the principle of the hardware system, Figure 12 shows the test system, and Figure 13 shows the hydraulic damper experimental device.

Experiment and Simulation Comparison.
After the experimental platform and control system are built, the overall vibration damping performance of the system is verified experimentally. In order to compare with the simulation results, the white noise signal with the same waveform is selected for the excitation signal of the shaker, and Figure 14 shows the comparison curve.
As can be seen from the graph, the experimental curve fluctuations compared with the simulation curve are bigger due to the inevitable error of the experimental device and the delay of data collection and processing. Comparing the upper and lower peak values of the displacement curve obtained by the experimental device and AMESim simulation, it can be concluded that the arithmetic mean of the error is 0.33 mm, and the error is within 10%. For the  acceleration curve shown in figure, the curve obtained by the experiment has a larger error compared with the simulation, which is mainly due to the error caused by the external uncertain disturbance. When the vibration curve is between 0 and 2 s, the vibration amplitude reaches a larger value, and the curve fluctuation tends to slow down as time goes on. e experimental graph shows that the experimental results of the semiactive suspension of the new magneticelectric hybrid suspension hydraulic damper verify the theoretical model. e above two effects are consistent, and both can achieve the vibration reduction effect, and the error is within the allowable range.

Analysis of Experimental Results.
In order to reflect the superiority of the semiactive suspension damping effect of the new hydraulic damper, the damping performance of the semiactive suspension and the passive suspension are compared. Due to the limitations of the experimental device, the vibration experiments of the hydraulic damper under the limited bandwidth white noise excitation signal when the electromagnet current is turned on and when the electromagnet current is not turned on are, respectively, carried out to simulate the actual vibration reduction environment of semiactive suspension and passive suspension. Figure 15 shows the comparison curve, Table 6 provides the comparison of peak and root mean square values of passive suspension and semiactive suspension about displacement curve, and Table 7 provides the comparison of peak and root mean square values of acceleration curves.
As can be seen from Table 6, the semiactive suspension system is obviously superior to the passive suspension in terms of the vibration reduction effect, with a difference of around 20%. In other words, compared with the passive suspension, the semiactive suspension with the new hydraulic damper has a certain advantage in the vibration reduction rate of around 20%.
On the one hand, the acceleration curve can directly reflect the stress state of the system, thus reflecting the changing characteristics of the system stiffness. On the other  Shock and Vibration hand, in the case that the acceleration curve fluctuates smoothly, the ride comfort of the car is better. As can be seen from Table 7, compared with the passive suspension, the semiactive suspension system has less fluctuation in the vertical acceleration of the vehicle body subjected to vibration, and the difference is about 10%. e experiment   shows that compared with passive suspension, the semiactive suspension of the new hydraulic damper has a certain advantage in the comfort of around 10%.
In order to get a better display of all the above curve results, Figure 16 shows the power spectral density diagram of semiactive suspension and passive suspension simulation and experiment. e PSDs show that the control force of the semiactive suspension has a better control effect in both lowand high-frequency regions.

Study on Variable Stiffness of the Damper
e purpose of introducing magnetoelectric hybrid suspension into the hydraulic damper is to make the stiffness of the system controllable. erefore, it is necessary to conduct an experimental comparison of the system with variable stiffness under different voltages, which can prove that changing the voltage of the input electromagnet can change the stiffness of the system and find the most suitable voltage value. Figure 17 shows the control program block diagram.
As can be seen from the Figure 17, the program block diagram for the closed-loop feedback control, control module has two main functions: data processing and preservation; the displacement signal and voltage signal conversion, filtering, and processing: intelligent discriminant of the vibrator vibration signal of different to choose the appropriate voltage amplitude of the signal input in the drive, which can adjust the vibration.
is experiment device of electromagnet rated voltage for 24 V, in order to choose the appropriate voltage signal, a smart choice in the control program module, small vibration amplitude, vibration wave as smooth as the premise, to select the optimal damping effect than that of the voltage signal value −18 V, as same as the above limited bandwidth white noise signal as excitation source vibration experiments. On the premise of stiffness change, two other comparable input    voltage signals, respectively, 12 V and 24 V, are selected and vibration reduction experiments are carried out under the same conditions. Figure 18 shows the results.
As can be seen from the displacement comparison curve in Figure 18, when the input voltage of the electromagnet is 18 V, the displacement and acceleration vibration curve of the damper are located in the relatively middle position, the amplitude of the damper fluctuates more stable, and the response speed is faster. When the input voltage of the electromagnet is reduced to 12 V, it can be seen from the experimental results and the quantitative relationship obtained from the previous analysis that the system stiffness is reduced compared with the former, and the amplitude of the vibration curve fluctuates less, but the response speed is slower. When the voltage of the input electromagnet is increased to 24 V, the stiffness of the system is relatively improved, and the amplitude of the vibration curve is increased. Meanwhile, the power consumption of the solenoid at the rated voltage of 24 V is larger. e experimental results show that when the input voltage signal is the optimal value of the system vibration damping effect, the system stiffness is higher, and the system vibration amplitude becomes larger when the voltage signal increases. When the voltage signal decreases, the system stiffness and amplitude decrease, but the system response speed is slower. e experimental results are consistent with the theoretical results above, that is, the stiffness of the system can be adjusted by adjusting the input voltage value.

Conclusion
A new type of magnetoelectric hybrid suspension hydraulic damper is proposed by introducing magnetoelectric hybrid suspension structure on the basis of the spring-damping suspension system. In the damper, the electromagneticpermanent magnet suspension structure is used to adjust the vibration of the system. According to the established structure model combined with its working principle, the method of joint analysis of magnetic circuit and oil circuit is adopted to obtain the expression of electromagnetic force and damping force. By analyzing the mathematical model of the 2-DOF 1/4 chassis vibration damping system, the mathematical model of the equivalent stiffness and equivalent damping of the system and the space state equation of the system are obtained.
According to the structural model, a two-degree-offreedom quarter vehicle AMESim simulation model is established, and the Simulink control model is established from the state equation in the vibration system. White noise with limited bandwidth is used as the road irregularity signal, PID is used as the control strategy, and AMESim/ Simulink cosimulation is carried out meanwhile. e vibration test bench is set up to verify the rationality of the theoretical model and the simulation model. e experimental results show that the semiactive suspension system equipped with the new magnetoelectric hybrid suspension hydraulic damper has a better damping effect than the ordinary passive suspension system, with the displacement amplitude reduced by 20% and the acceleration amplitude reduced by 10%. By adjusting the input voltage, the vibration curve of the system can be changed, and the optimal voltage value of the system is 18 V, which proves that the stiffness of the designed damper has certain controllability.
Data Availability e dynamic model data used to support the findings of this study are included within the article and references.

Disclosure
Kun Hu and Fei Li are the co-first authors.

Conflicts of Interest
e authors declare that there are no conflicts of interest.

Authors' Contributions
Kun Hu and Fei Li contributed equally to this work.