Modeling and Parameter Identification of MR Damper considering Excitation Characteristics and Current

Smart structures such as damping adjustable dampers made of magnetorheological (MR) fluid can be used to attenuate vibration transmission in vehicle seat suspension. -e main research content of this paper is the nonlinearity and hysteresis characteristics of the MR damper. A hysteretic model considering both excitation characteristics and input current is proposed to fit the damper force-velocity curve for the MR damper under different conditions. Multifactor sensitivity analysis based on the neural network method is used to obtain importance parameters of the hyperbolic tangent model. In order to demonstrate the fitting precision of the different models, the shuffled frog-leaping algorithm (SFLA) is employed to identify the parameters of MR damper models. -e research results indicate that the modified model can not only describe the nonlinear hysteretic behavior of the MR damper more accurately in fixed conditions, compared with the original model, but also meet the fitting precision under a wide range of magnitudes of control current and excitation conditions (frequency and amplitude). -e method of parameter sensitivity analysis and identification can also be used to modify other nonlinear dynamic models.


Introduction
Vehicle ride comfort plays an important role in vehicle dynamics research, in which the transmission of vibration through seat exerts a crucial influence. Traditional passive seat suspensions are very limited in terms of ride comfort and cannot meet the increasingly stringent comfort requirements. In recent years, the semiactive seat suspension system is becoming more and more popular due to simple structure, low power consumption, good vibration damping behavior, and high reliability [1][2][3][4][5][6].
ere are many types of semiactive suspensions such as throttle variable damper suspension, air suspension, electrorheological damper suspension, and magnetorheological damper suspension [7,8]. erefore, an MR damper is a fluid damping device, widely used in vibration control of civil engineering structures and mechanical systems because of its wide damping adjustable range, short adjustment time, high reliability, and low cost [9][10][11][12]. e structure of the magnetorheological damper is similar to the structure of an ordinary vehicle cylinder damper. e medium inside the MR damper is a magnetorheological fluid. e MR fluid can reversibly change from free flowing, linear viscous liquids to non-Newtonian fluid with a controllable yield strength when exposed to a magnetic field [13].
However, the magnitude of the generated damper force is related to many factors that involve many nonlinear elements and increase model complexity of the MR damper. e causes of nonlinearity are listed as follows: (1) different magnetoresistivity of electromagnetic circuits; (2) hydrodynamics of MR fluid flowing in piston gap; (3) the chemical properties of MR fluid; and (4) the relationship between the yield stress generated by MR fluid and applied magnetic field strength.
Various papers have established MR damper mechanical models by various proposed methods. According to relevant literature, MR damper force models are divided into three main groups [10,[14][15][16]: pseudostatic models, parametric models, and nonparametric models. In pseudostatic models, the damper force is calculated by the pressure difference between the two ends of the damper piston, structural parameters, MR fluid parameters, and magnetic field strength [17]. is type of model can generally reflect the forcedisplacement characteristics properly, but the fitting of the force-velocity characteristics is not accurate enough. Pseudostatic models are mainly used for the design of MR dampers, but not suitable for the research of semiactive suspension control. Parametric models and nonparametric models are mainly distinguished according to whether the parameters make sense in physical. Nonparametric models generally consider that the damper force consists of a polynomial with input current and damper velocity, or neural network model [18]. And, in parametric models, the damper force is calculated by physically meaningful parameters.
ere are many types of parametric models, among which the Bingham model is the most basic one. Although this model can describe the relationship between the damper force and displacement well, it cannot reflect the hysteresis relationship of force-velocity. At present, the Bouc-Wen model is widely used, which can exactly describe the hysteretic characteristics of the MR damper force versus piston velocity by designing the hysteresis variable z. However, the nonlinear differential equations contained in the Bouc-Wen model inevitably increases difficulty to parameter identification. In 1997, Spencer and Dyke [19] proposed a modified Bouc-Wen model, which further improved the fitting accuracy of the force-velocity hysteresis curve. ere is a wide range of parametric models: nonlinear hysteretic biviscous model, generalized sigmoid hysteresis model, hyperbolic tangent function-based model, and so on. However, the previous models have shortcomings in terms of model accuracy, complexity, or limitation [14,16,[20][21][22]. e purpose of this paper is to build a detailed model of the MR damper used for semiactive suspension control. e paper is structured as follows: in Section 2, the mechanical performance test of the MR damper was performed to study the characteristics of damper force in different excitation inputs. In Section 3, several classic MR damper models were introduced in detail, and the more suitable model was selected. In Section 4, the sensitivity analysis method of parameters of the hyperbolic tangent model was applied to illustrate the importance of every model parameter, as references for model extensions in subsequent chapters. e shuffled frog-leaping algorithm used for parameter identification is explained in Section 5. Subsequently, the identification results of damper hysteretic model parameters are discussed and an extended hysteresis model is achieved and verified in Section 6. In the end, conclusions are presented in Section 7.

Experimental Setup.
Two types of MR dampers, shown in Figure 1, were employed for dynamic tests in this study. e strokes of the pistons of these two dampers are 74 mm (long stroke, RD-8041-1) and 55 mm (short stroke, RD-8040-1), respectively. Except the length of stroke, the two dampers have the same structure parameters and MR fluid. e damper consists of two liquid chambers that are separated by a piston with small orifices. e rated value of the input current is 1 A in the electromagnet coil around the piston head. Figure 2 shows the experimental setup including actuator, computer closed-loop control electro-hydraulic servo system, force sensor, displacement sensor, fixture, and constant current power supply. e damper was driven by a Schenck actuator, and the position sensor and force transducer were used to collect the measured displacement and damper force of the two dampers.
Sine vibrations, with the frequency range from 0.5 Hz to 9 Hz and vibration amplitude of 5 mm, 10 mm, 15 mm, and 20 mm, were provided by the actuator. And, the current inputs ranged from 0 to 1 A with 0.2 A intervals. e specification of excitations is shown in Table 1.

Experimental Results.
In the following analysis, only the results of the RD-8040-1 damper were discussed, in which the mechanical behaviors of the two dampers were similar, as was shown in the experiment. e responses of the MR damper applied a 20 mm 2 Hz sinusoid excitation which are shown in Figure 3, at six constant current values, 0 A, 0.2 A, 0.4 A, 0.6 A, 0.8 A, and 1 A. Figures 3(a) and 3(b) illustrate the force-displacement (f-d) and the force-velocity (f-v) characteristics, respectively. It is obvious that the f-d curve is approximately rounded rectangle. ere are evident bumps near the zero-displacement area, which is caused by the elasticity of the damper. e effect of changing current value is readily observed, and the damper force curve rises monotonically as the current value increases. e increase of current value is found to be nonlinear in nature. Under 20 mm 2 Hz sinusoid excitation, the maximum damper force increases by 450.71N when current value increases from 0.2 A to 0.4 A. Similarly, for increase from 0.4 A to 0.6 A and from 0.6 A to 0.8 A, the maximum damper force increases by 297.02 N and 324.06 N severally. When the current value increases from 0.8 A to 1 A, the increment of maximum damping force is 64.52N. It can be also observed that the magnetic field intensity generated by the electromagnetic coil is close to the magnetic saturation intensity once the   Shock and Vibration current value increases above 0.8 A, reflected in which the increment of the damper force of the MR damper becomes smaller. e damper force-velocity curve is symmetrical and nonlinear, with apparent hysteresis characteristics. When the velocity of the piston increases, the damper force increases along the lower branch curve. e hysteresis curve of damper force moves counterclockwise with damper velocity. When the input current value of the coil increases, the hysteresis loop of the damper force-velocity of the MR damper gradually increases. Figure 4(a) shows the force-displacement characteristics with 0.8 A current applied under 5 mm sine vibration and various frequency (1 Hz, 2 Hz, 3 Hz, 4 Hz, 5 Hz, 6 Hz). It can be observed that the maximum damper force increases when increasing frequency, keeping the current and amplitude invariable. As the frequency grows higher, the area enclosed by the damper force-displacement loops increases, and the energy dissipation in each vibration cycling increases. As is shown in Figure 4(b), the damper force at the same velocity is approximately equal and shows the Bingham plastic behavior of the MR damper when the velocity and acceleration are in the same direction. However, the damper force in the same velocity decreases when frequency grows higher, in case the velocity and acceleration are reversed. is phenomenon of the MR damper is because the MR fluid presents viscoelasticity at low velocity and viscoplasticity at a high velocity. Furthermore, the size of the hysteresis loop mainly depends on the inertia effect of damper. e area of hysteresis loops increases due to the large inertia resulting from the high acceleration as the velocity direction changes.
In Figures 5(a) and 5(b), the area enclosed by the damper force-displacement curve increases as the excitation amplitude increases under the same input current and frequency, and the area of force-velocity hysteresis loop becomes larger as well.
Based on the characterizations of the MR damper above, it can be observed that the MR damper behaves strong nonlinearity. All of the excitation amplitudes, frequencies, and input currents have impacts on the mechanical characteristics of the MR damper.
It is difficult to accurately describe the MR damper model by conventional methods. e establishment of a comprehensive dynamic model of the MR damper has an important influence on the control of the semiactive seat suspension system. In this study, a detailed hysteretic model was established to describe the mechanical characteristics of the MR damper for the purpose of vehicle suspension control.

The Mechanical Models of the MR Damper
As a controllable fluid, MR fluid is a Newtonian fluid that can freely flow in the initial state and turns into non-Newtonian fluid when exposed to a magnetic field. e change of MR fluid viscous is reversible and reflected by a significant variation in the damper output force.
A lot of models have been proposed to describe the nonlinear behavior of the damper force versus velocity of MR dampers.
Bouc-Wen model is widely used in the modeling of hysteretic systems ( Figure 6). It can exhibit various hysteretic behaviors by changing its model parameters and has very good adaptability. e expression is shown as follows: where, α, β, c, and n are the model parameters and z is a hysteresis variable. e bias force generated by the accumulator inside the MR damper is converted to the initial offset x 0 of the spring k 0 . e Bouc-Wen model has good fitting accuracy for the damper force-velocity characteristics of MR dampers. However, it makes great difficulty to the parameter identification and may affect the robustness of the control system, due to the existence of the differential term _ z. Research shows that the Bouc-Wen model can better predict the force-displacement and force-velocity characteristics of MR dampers, compared with the Bingham model. But, in the low-speed area, the model is not accurate enough to fit the velocity hysteresis characteristics.
Spencer proposed a modified Bouc-Wen model in order to more exactly predict the damper response in low velocity. e damper force is given as  where c 0 and c 1 represent the viscous damping at low velocity and high velocity, respectively; k 1 represents stiffness at high speed; and x in is the internal variate of displacement, as shown in Figure 7. Other parameters are consistent with the Bouc-Wen model. Although the accuracy of the model is improved, the complexity of the model also increases significantly. e newly brought differential intermediate variable in this model brings greater challenges in parameter identification and the robustness of control system. Kwok and Ha [30] proposed a damper model in 2006. ey used a hyperbolic tangent function to represent the hysteresis element z (see Figure 8). e viscosity and stiffness elements are expressed by linear functions c 0 _ x and k 0 x. Due to lack of differential equations, the calculation efficiency of this model is promoted, and parameter identification becomes simple to operate. e model is expressed as follows: rough calculation and analysis, it can be found that the hyperbolic tangent model can accurately fit the velocitydamping force hysteresis curve of the MR dampers. Figure 9 shows the comparison of the calculation results of the two models. e fitting accuracy of the hyperbolic tangent model is higher than the Bouc-Wen model especially in the velocity transition zone. e output damping force error RMS of the two models compared to the test damping force are 30.88N and 49.49N, respectively.
By comparing different hysteretic models, it can be found that the hyperbolic tangent model is simple in structure and does not contain differential terms. Meanwhile, it can fit the velocity-damping force hysteretic

Parameter Sensitivity Analysis
Sensitivity analysis of the mathematical model parameters is not only critical to model validation but also can provide guidance for research of parameters analysis. e analysis of hyperbolic tangent hysteresis model shows that it can predict response for fixed current input and vibration excitation. However, the model parameters under different currents and excitation conditions are different. In order to establish a general model suitable for different working conditions, a hysteresis model including current, amplitude, and frequency characteristics is needed. e method of sensitivity analysis can distinguish the important parameters from unimportant parameters in the mathematical model [31]. Due to the large number of parameters in the hyperbolic tangent model, the sensitivity analysis about the contribution of the parameters to the MR damper model output is carried out to filter out the parameters with poor sensitivity. e relationships between the sensitive parameters and the excitation characteristics or input current are established, so as to obtain the extensional MR damper model that is suitable for different excitation conditions.
In sensitivity analysis, the sensitivity coefficient can be expressed as follows: where F i and F i ′ represent the model output force values at the ith sampling spot of parameters before and after changing; N � total number of sampling spot; and S � sensitivity of parameters. e sensitivity of the parameter increases as the value of S becomes larger. e main parameters that have relatively significant influence on the sensitivity coefficient of damper force were obtained via adjusting each parameter within its allowable range, so as to study the degree of influence of these parameters on the accuracy of the model. Furthermore, in subsequent model extensions, model design was based on the ranking of parameter sensitivity, focusing on parameters with larger sensitivity coefficients, and establishing the relationships between model parameters and excitation characteristics or input current. In this way, the number of model parameters is minimized as much as possible on the premise of ensuring model accuracy. erefore, the complexity of the model is reduced, and the computational efficiency of the model is improved.
ere are plenty of sensitivity techniques in previous research, and the results of sensitivity ranking may exist slightly disparity for different methods [32]. In this paper, multiparameter sensitivity analysis based on BP neural network was put in use. e neural network (NN) models are well known for strong self-learning adaptability, high robustness, and excellent fault tolerance. It can reflect the implicit feature rules by the connection weights between neurons and solve many nonlinear problems including complex information, unclear rules, and noise pollution. e neural network models can evaluate the sensitivity values of model parameters in the situation that multiple parameters change simultaneously in complex models [33].
In this paper, the BP neural network was employed to quantify the sensitivity of parameters of the damper force model. e weights and thresholds of the input layer to the hidden layer and of the hidden layer to the output layer were extracted through training the BP neural network with sample database. e BP neural network model and weights represented the mapping relationship between each parameter and output. e method of Garson (equation (5)) was applied to compute the sensitivity values of model parameters.
where W ij � the connection weight from input unit to hidden unit and W jk � the connection weight of hidden unit to output unit. e sine excitation at amplitude 10 mm and frequency 4 Hz and input current 0.8 A were loaded on the MR damper to analyse parameters sensitivity. e main steps of sensitivity analysis method based on BP neural network are as follows: Step 1. Structure design of neural network: a three-layer neural network was created, including input layer, hidden layer, and output layer, to calculate the sensitivity values of each model parameter. ere are six units in input layer: α, β, δ, F 0 , c 0 , and k 0 , respectively.
e sensitivity values of parameters were calculated by equation (4). e number of units in the hidden layer was calculated by equation (6) where, n, n 1 , and m represent separately the number of input unit, hidden unit, and output unit; a � any  constant between 1 and 10. After calculation, the number of units in the hidden layer was determined to be 10. e structure of neural network is illustrated in Figure 10.
Step 3. Analysis of the neural network model: by analysing the calculation results of the neural network, conclusions can be made that there are more than 90% samples whose errors are within 5%. Obviously, the mapping established through the neural network can get accurate relationship between model parameters and damper force output. e weights (W ij and W jk ) are listed in Table 2. Table 3 lists the sensitivity coefficients of each model parameter calculated by the methods of Garson (equation (14)). It is noteworthy that the parameter a has the greatest impact on the damper force output, with the sensitivity value of 39.43%. And, β, c 0 , and d rank behind, with values of 30.68%, 14.76%, and 14.71%, respectively. e remaining two parameters have only slight effect on the damper force output. e calculation result of sensitivity analysis show that the effects of parameter α, β, δ, and c on the model output are significant ones; meanwhile, the parameters F 0 and k 0 can be set to constants due to the low sensitivity.

The Shuffled Frog-Leaping Algorithm (SFLA)
To describe the force-velocity hysteretic characteristics of the MR damper, the parameters in the MHT model need to be identified.
Parameter identification is a key issue in this research due to the complexity and nonlinearity of the model. In order to accurately determine the parameters, researchers have successively proposed a series of identification methods aimed at nonlinear models, such as genetic algorithm (GA), particle filter algorithm (PA), and least square method (LSC) [34].
In this study, the processes of parameter identification utilizing shuffled frog-leaping algorithm and the experiment data fitting were carried out [35,36]. e shuffled frogleaping algorithm was firstly proposed in 2003 by Lansey and Eusuff [37] to solve the combinatorial optimization problem.
is algorithm combines the advantages of particle swarm algorithm and memetic algorithm and is a new metaheuristic group evolution algorithm. It has been successfully applied in many fields especially in solving multiobjective complex function optimization problems. Figure 11 shows the parameter identification process of the model.
Although there are some resemblances between SFLA and genetic algorithm (both are updated iteratively through fitness sorting), the SLFA has a faster spread speed. Memes in SFLA can flexibly use different mechanisms to spread information from one member of the memeplex to another, while genes in GA can only pass from parents to offspring. e evolution in genetic algorithm is limited by the number of offspring that a single parent can produce.
Steps of the shuffled frog-leaping algorithm are as follows: Step 1. Initialize parameters and population: set t � 1; determine the population size F and randomly generate the initial population X (t); determine the number of subpopulations m and the number of individuals in each subpopulation n; F � m * n.
Step 2. Rank frogs: calculate the fitness values of each individual frog in the current population X (t).
Step 3. Partial update iteration: t � t + 1. Sort the entire frog group according to their fitness and perform a local search in each subpopulation to update the worst frog.
where SW and NL are the regeneration step size and the updated position of the worst frog, respectively; SW max � the maximum step size of frog (3) If the two substeps above cannot meet the iteration requirements, the frog that is in the worst position in the subpopulation will jump to a random position (rp) in the solution space.
Step 4. e iteration ends after all the subpopulations have completed the partial search. Mix all the frogs together and regroup according to the updated fitness value.
Step 5. Determine whether the current result meets the convergence demand. If not, go to step 2 and resort and divide the memeplex to perform the local search process. If it does, the process of algorithm ends.

Parameter Identification.
e dynamic performance of the MR damper is not only related to input current but also influenced by the characteristics of the excitation (amplitude and frequency) according to the research in previous sections. e identified parameters extracted from a fixed working condition cannot accurately fit for other excitation conditions. A general model that is suitable for different excitation conditions should be built by establishing the relationship between model parameters and loading conditions including input current, vibration frequency, and amplitude.
Although several studies have made effort to introduce the applied load to the damper model, most of these papers only study the influence of input current on the model output, without considering excitation characteristics [38]. Cheng and Chen [39] proposed a model containing amplitude and frequency.
Noteworthy, the amplitude and frequency of vibration excitations cannot be measured in real time, especially under random excitation. erefore, it is necessary to explore the relationship between excitation characteristics and measurable parameters during modeling. Wang and Ma [21] proposed a relational expression between movement states (displacement, velocity, and acceleration) of an MR damper and excitation characteristics. e maximum velocity (v m ) of the damper has a linear relationship with the sine excitation amplitude and frequency. v m can thus be expressed by where c 0 is a constant coefficient and a and f represent the amplitude and frequency, respectively. To verify the equation, the maximum velocity values under a series of working conditions were sorted out and then observed the relationship between the v m and amplitude or frequency (see Figure 12). It is evident that the v m amplitude curve and v m frequency curve both show linear relationships. As for v m variation under fixed amplitude (5 mm, 10 mm, 15 mm, and 20 mm), it increases with the increase of frequency straightly. Similarly, the linear increase relationship of v m versus amplitude is observed under each fixed frequency (1)(2)(3)(4)(5). Figure 13 demonstrates that the value of v m does not depend on the input of current approximately under each constant amplitude and frequency.
Because the variable v m cannot be directly measured in real time either, and the damper instantaneous displacement and acceleration are connected with v m as follows:

Shock and Vibration
where _ x is the velocity of the damper and € x and x represent instantaneous displacement and acceleration, respectively. erefore, the excitation frequency and amplitude are expressed by measurable dynamics variables of the damper (when the vehicle is operating in real time) by taking advantage of equations (10) and (11).
Here, we firstly identified the parameters of the damper model from experiment data under a wide range of exciting conditions. e shuffled frog-leaping algorithm was utilized to obtain the accurate model parameters. And then, the correlation rules between each parameter and the excitation characteristics (current and v m ) were excavated in purpose of building the general damper model. e relationships between parameter α and the intensity of the input current or variable v m are shown in Figure 14. It is clearly observed that the parameter a can be linearly represented with the input current (see Figure 14(a)). e quadratic polynomial has good fitting effect on the a-v m curve in Figure 14(b), described as e relationships between parameter ß and excitation characteristics or control current are shown in Figure 15.
e results show that ß decreases nonlinearly as the excitation velocity amplitude increases. It can be found that the accuracy of the power function fitting is adequate. In the meanwhile, ß almost remains unchanged with the variation of control current. erefore, ß can be expressed as follows: Noticeably, in equations (12) and (13), the coefficients α 1 , α 2 , α 3 , α 4 , α 5 , β 1 , and β 2 are constants. Figures 16(a) and 16(b) reveals that parameter δ depends on only the current value, rather than vibration characteristic v m . For fixed vibration input (v m ), the parameter δ increases as the control current increases. A first-order polynomial can approximate the relationship between δ and current precisely.
Observations from Figures 16(c) and 16(d) are that the parameter c 0 linearly rises as current increases, while the effect of v m on parameter c 0 can be neglected. e parameter c 0 is obtained as follows: e variations of other parameters, k and F 0 , have nothing to do with the input current and excitation characteristics (Figure 17), and the sensitivity values of these parameters to the model output are quite low without exception. erefore, it can be regarded as constants.
Plugged the calculated results in the hyperbolic tangent model, equations are shown as follows: e number of parameters to be identified increases from 6 (α, β, δ, F 0 , c 0 , and k 0 ) to 13 (a 1 , a 2 , a 3 , a 4 , a 5 , ß 1 , ß 2 , d 1 , d 2 , F 0 , c 01 , c 02 , and k 0 ). e parameters were reidentified under every input current and vibration condition, using the shuffled frog-leaping algorithm. e recalculated parameters are shown in Table 4.  e hyperbolic tangent model only considering current variable was established through the same process as above: (17) Figure 18 shows velocity-damping force characteristics of different models under 5 mm, 2 Hz sine excitation. e black curve represents the damping force measured by the MR damper mechanical test, the red curve represents the hyperbolic tangent model only considering the input currents, and the blue curve represents the hyperbolic tangent model considering both excitations and currents.
We can notice that the damper force calculated with the hyperbolic tangent model which only considers current variable has a large margin of error whether in low or high velocity area. e reason leading to great error is that the hysteretic model parameters do not include the elements related to excitation characteristics, while the parameters α, β, and c 0 are significantly influenced by the excitation frequency and amplitude as well. Shock and Vibration erefore, it can be found that the improvement effect of the model is better as the input current increases.

Model Validation.
A two-degree-of-freedom seat suspension model is established to verify the damping control of the proposed model applied to the semiactive suspension. e model structure is shown in Figure 19. e motion differential equation of the seat suspension model shown in Figure 19 is shown as follows: (18) e parameter meanings in equation (18) Figure 20 shows the simulation results under random excitation. e included frequency range is 1-10 Hz. e acceleration of the seat upper surface is reduced when the semiactive control is applied. When the modified MR damper model is applied to the control system, the acceleration is further reduced, as the black curve in the figure shows.
In order to see the vibration reduction effect more clearly, the vibration transmissibility under different frequencies is shown in Figure 21. Semiactive suspension of the MR damper can effectively reduce the transmission of vibration, especially at frequencies sensitive to the human body (4)(5)(6)(7)(8). At the resonance peak frequency, the original tanh model reduces the transmissibility by 0.79, and the modified model further reduces the transmissibility by 0.18. Figure 22 shows the comparison of the control current of the two models. In order to clearly see the difference, the     sampling time of the current is set as 0.1 s. It can be seen that the output current of the two models is different. According to , it can be determined that the control current calculated by the modified model is more accurate. e accurate model can further improve the vibration reduction effect of the whole semiactive suspension control system.  Figure 19: Semiactive seat suspension system model.

Conclusions
is paper proposed a detailed, accurate damper model, which was based on the hyperbolic tangent model, to solve the problems of nonlinearity and hysteresis characteristics of an MR damper.
By employing the mechanical test of the MR damper, the fitting effect of the different hysteresis model was verified. e hyperbolic tangent model was selected as the research object. e calculation result of sensitivity method was analysed comprehensively to confirm the parameters that have important impacts on the output damper force of the hyperbolic tangent hysteretic model, in order to provide a guidance for the design of the extended model. It is proved that the neural network model is appropriate to analyse parameter sensitivity.
After analysing and comparing different hysteretic model, we identified the model parameters using the shuffled frog-leaping algorithm. By establishing the relationship between sensitive parameters and loading characteristics (input current, excitation amplitude, and frequency), the modified model and its parameters were determined. e results show that the modified model is an efficient hysteretic model to characterize the damping property of the MR damper under a wide range of control current and excitation conditions. e precise model can greatly improve the damping effect of the whole semiactive suspension control system.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.