Neutral networkPID control algorithm for semiactive suspensions with magnetorheological damper
Cheng Lin^{1} , Wei Liu^{2} , Hongbin Ren^{3}
^{1, 2, 3}School of Mechanical Engineering, Beijing Institute of Technology, Beijing, 100081, China
^{1, 2}Collaborative Innovation Center of Electric Vehicles in Beijing, Beijing, 100081, China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 8, 2015, p. 44324444.
Received 1 July 2015; received in revised form 18 August 2015; accepted 28 August 2015; published 30 December 2015
JVE Conferences
In this paper, a semiactive suspension control system based on MagnetoRheological (MR) damper is designed for a commercial vehicle to improve the ride comfort and driving stability. A mathematical model of MR damper based on the BoucWen hysteresis model is built. The mathematical model could precisely describe the characteristics of MR damper compared with the bench test results. The neural networkPID controller is designed for the semiactive suspension systems. According to the numerical results, the proposed controller can constrain vehicle vibrations and roll angle significantly. A detailed multibody dynamic model of the light vehicle with four semiactive suspensions are established, and an actual vehicle handling and stability tests are carried out to verify the control performances of the proposed controller. It can be concluded that MR semiactive suspension systems can play a key role in coordination between the ride comfort and handling stability for the commercial vehicle.
Keywords: magnetorheological (MR), semiactive suspension, neutral networkPID.
1. Introduction
A good vehicle suspension system is to isolate the passengers from external road disturbances and internal vibration from engine to improve the ride comfort. The classical suspension is made up of spring, damper elements and a set of mechanical elements which link the suspended vehicle body to wheels. Most current suspension systems are conventional passive types [1], however, vehicle ride comfort and handing stability are conflicting for the traditional passive suspension systems. Good vehicle ride comfort performance can be achieved with soft damping characteristics especially in the range of 48 Hz, which is known to be a sensitive frequency range of human body. But the soft damping cannot keep the tire in contact with the road surface and will deteriorate the vehicle safety and stability [2]. Good tune and design of a passive suspension characteristics can to some extend optimize and coordinate the ride comfort quality and stability, but cannot eliminate these compromises thoroughly [3]. With the development of automobile techniques, the controllable suspension attracts more and more attentions in recent years [46], such as active and semiactive suspension systems. Advanced active suspensions can offer opportunities for substantial improvements in ride and handling stability as well as overall vehicle posture by exerting independent forces to vehicle body through a computercontrolled force generator, where additional power sources are needed, such as pumps and compressors. This also implies high energy consumption, complexity and additional expense. These drawbacks have limited its development and application in modern automobile industry. Semiactive suspension is a good compromise between active and passive suspension performance, complexity, and expense. The semiactive damper generates force in a passive manner, but the damping forces can be actively controlled. The damping force is modulated in accordance with the operational conditions, which is controlled by certain logic from sensors connected to CPU. It can be adjusted by changing the orifice size or oil viscosity separately or continuously, such as continuously variable dampers, magnetorheological suspension. In addition, semiactive suspension is more stable and failsafe; because it still can work under pure passive mode in case of control system failed. Due to its less energy consumption and better control priorities under extreme driving conditions, it has been considered as a good alternative [7].
A wide range of control techniques have been developed for active and semiactive suspension systems in last two decades. The skyhook control strategy is proposed by Karnopp, et al. in 1974 [8] and since then many control strategies are proposed based on skyhook logic. An adaptive semiactive control algorithm combined with the online sprung mass estimation was proposed by Song [9]. Moradi, et al. [10] designed an adaptive PIDslidingmode faulttolerant control for full car suspension systems. Some intelligent approaches are also applied in the control because of the nonlinear and uncertainty characteristics of the vehicle, such as genetic algorithms [11]. Although these intelligent approaches have been applied, the mathematical proof for stability is still not to be demonstrated now; and the system stability is especially important for the active suspension systems [12].
In recent years, the controller design for suspension systems with complex nonlinear characteristics and uncertain disturbances becomes a challenging work, and the applications of neural network in the automotive industry have been achieved many remarkable achievements [1315]. The approximate neural network model can substitute the actual model in many applications especially in control system design. With appropriate offline training, the neural network controller can achieve a good coordination between ride comfort and handing stability for the semiactive suspensions [16, 17].
The contribution of this paper is to propose a neural network controller for semiactive suspension system with MR dampers. The fullcar model of the light vehicle including suspension and steering system, is built in SIMPACK. A mathematical model of MR damper is established and verified combined with the bench test. Compared with the experimental data, the proposed model is precise. Combined with the neural networkPID control algorithm and mathematical model of MR damper, a closed loop control system for semiactive suspension is proposed.
The rest of paper is organized as following: a nonlinear vehicle model and mathematical model of MR damper are presented in Section 2; while in Section 3 the skyhook reference model and neutral networkPID control algorithm are introduced in details; the efficacy of the proposed controller is validated by simulation presents in Section 4. Vehicle tests are carried out in Section 5; at the end of the paper, the conclusions are given.
2. Nonlinear dynamic modeling
2.1. Full vehicle model
The full scale vehicle suspension dynamic model is illustrated in Fig. 1. The sprung mass consists of the chassis, which are including passenger and internal components. It connects the suspension systems to four wheels (unsprung masses). The vehicle body is free to heave, pitch, and roll. The four wheels are free to bounce vertically relatively to the vehicle body. The damping forces can be adjusted by controlling the current of MR damper according to designed control logic. ${m}_{u}$ is the unsprung mass, which is supported by tire which is modeled as linear spring with stiffness coefficient ${K}_{t}$. The displacement of sprung mass and unsprung mass are denoted as ${z}_{s}$ and ${z}_{u}$ respectively. A road disturbance is integrated by Gaussian white noise, denoted as $q$.
The differential equations of the vertical, pitch and roll motions are expressed as:
where, ${F}_{sus}$ is the suspension force, which consists of the spring force and damping force, can be described as:
where, ${u}_{ij}$ is damping force control current; $fl$, $fr$, $rl$, $rr$ denote frontleft, frontright, rearleft, rearright of the vehicle body.
Fig. 1. Fullvehicle semiactive suspension model
The dynamic equations of wheels’ vertical motion are as follows:
Assuming that vehicle body is rigid, and the pitch angle $\theta $ and roll angle $\phi $ are small, the suspension deflections of four corners are:
$\mathrm{\Delta}{{z}_{s}}_{rl}={z}_{s}+b\mathrm{s}\mathrm{i}\mathrm{n}\theta +\frac{B}{2\mathrm{s}\mathrm{i}\mathrm{n}\phi}{z}_{url},\mathrm{\Delta}{{z}_{s}}_{rr}={z}_{s}+b\mathrm{s}\mathrm{i}\mathrm{n}\theta \frac{B}{2\mathrm{s}\mathrm{i}\mathrm{n}\phi}{z}_{urr}.$
2.2. Modeling of MR dampers
A MR damper is a damper filled with magnetorheological fluid, which can be controlled by a magnetic field. MR fluid is a material that responds to an applied magnetic field with a significant change in its rheological behavior. When the magnetic field is applied, the properties of such a fluid can change from a freeflowing, low viscosity fluid to a near solid, and this change in properties takes place in a few milliseconds and is fully reversible. The damping characteristics of the shock absorber can be continuously controlled by varying the control current of the electromagnet. The MR shock absorber is notably applied in semiactive vehicle suspensions which may adapt to road conditions. The shock absorber is in the form of double tube structure, as shown in Fig. 2. The magnetic intensity of the valve can be adjusted by controlling the current of the coil, and the fluid rheological effect will occur, the damping force could be adjusted continuously. To study the damping characteristics of MR shock absorber, the bench test is carried on under sine excitation, the amplitude is 30 mm, the frequency is 0.3 Hz, and the max velocity is 0.057 m/s. The field current, 0 A, 0.2 A, 0.4 A, 0.6 A, 0.8 A, 1.0 A, 1.5 A, 2 A are applied, as shown in Fig. 3.
Fig. 2. MR damper schematic diagram and physical contrast diagram
Fig. 3. Bench test of MR damper
The nonlinear model of MR damper has been described by many mathematical models, such as Bingham model and nonlinear hysteretic biviscous model, et al. [18]. One of the most popular is the BoucWen hysteresis model, originally deduced by Bouc in 1967 and later modified by Wen in 1976 [19]. Later on, Spencer proposed a phenomenological model for MR dampers, a more accurate one based on the BoucWen model [20], shown in Fig. 4. This model predicts the forcedisplacement behavior of the damper well, and it describes forcevelocity behavior that more closely resembles the experimental data. According to the nonlinear forcevelocity and forcedisplacement responses of MR damper, the mathematical model of MR damper can be expressed as:
where, ${F}_{d}$ is the damping force of MR damper, and the values of unknown factors ($A$, $\beta $, $\gamma $, $n$, ${c}_{0}$, ${c}_{1}$, $\alpha $, ${k}_{0}$, ${k}_{1}$) should be determined by test data. The results of parameter identification are obtained by using the optimization algorithm shown in the Table 1.
Fig. 4. BoucWen model of MR damper
Table 1. Parameters of Phenomenological model
Symbols

Values

$A$

31

$\beta $

2.6 mm^{2}

$\gamma $

2.9 mm^{2}

$n$

2

${c}_{0}$

29 N·s·mm^{1}

${c}_{1}$

1400 N·s·mm^{1}

$\alpha $

430 mm^{2}

${k}_{0}$

–300 Nm

${k}_{1}$

–298 Nm

Fig. 5 shows the displacementforce characteristics of the MR shock absorber with the increasing intensity of control current, which implies that the dissipation energy in a vibration cycle would be increasing. Fig. 6 describes the relationship between the damping force and the velocity of the pistonrod of the absorber, when there is no magnetic field applied, the MR fluid would be the Newtonian flow, and the damping force proportionate to the pistonrod velocity; under the same velocity of the pistonrod condition, the damping force would be nonlinear increasing with the current increasing.
Fig. 5. Displacementforce characteristics of MR damper
Fig. 6. Velocityforce characteristics of MR damper
Fig. 7 shows the comparison of theoretical and experimental results. It can also be found that the theoretical predictions are in close agreement with the experimental results. It means that the mathematical model can be used to precisely describe the characteristics of the prototype MR damper.
Fig. 7. Comparison of theoretical model and experimental data for MR damper
Fig. 8. The configuration of skyhook reference model
3. Design of semiactive controller
3.1. Skyhook reference model
The ideal skyhook is a comfortoriented control policy. The skyhook damper control is assumed to be approximately realized by a state feedback control scheme. The skyhook reference control model is as shown in Fig. 8.
The onoff skyhook control logic can be expressed as:
where, ${C}_{sky}$ is the skyhook damping coefficient.
3.2. Neutral networkPID Control
The aim of semiactive controller is to control the motion of the vehicle body and the vertical vibration of unsprung masses to adapt different operating conditions. Generally speaking, the ride comfort and handing stability are conflicting for the traditional passive suspension. The main purpose of semiactive controller is to adjust the damping force according to the designed control algorithm determined by the transient suspension states, so as to reduce the vertical vibration of suspension and maintaining proper vehicle body posture. PID control is a more classic control algorithm, composed by the proportional, integral and differential links. When the vehicle roll rate is too large to cause the instability, the incremental PID controller can adjust the damping of suspension to produce a corrected roll moment. However, there are couplings among heave, roll and pitch motion of vehicle, and it is very difficult to improve the ride comfort and handing stability of automobile simultaneously. Therefore, the neural network methodology is adopted here to identify the different driving conditions and adjust the proportional, integral and differential coefficients (${K}_{p}$, ${K}_{i}$ and ${K}_{d}$) of PID controllers. After being trained by using a large amount of sample data, neural network algorithm can predict the state of the vehicle, and make timely and effective control strategy adjust to cope with the road excitation and the state of vehicle changes. Finally, a hybrid damping force control scheme based on PID and neural network control algorithm is proposed, shown in Fig. 9.
According to Eqs. (1)(7), the system function can be rewritten as:
where,$\mathbf{u}=[\begin{array}{llll}{u}_{fl}& {u}_{fr}& {u}_{rl}& {u}_{rr}\end{array}{]}^{{\rm T}}$ are the damping control forces from control algorithm; $\mathbf{q}$ is the road disturbance input.
Fig. 9. PID neural network control algorithm block diagram
The input of this control algorithm, shown in Fig. 9, is the tracking errors of ideal value from reference model and actual value, and the output of this control algorithm is the added value of the controlled object. The incremental PID control methodology is adopted in this paper, which provides the damping forces for semiactive suspensions. The PID control algorithm can be expressed as:
where, $e\left(k\right)$ is the deviation of the actual and reference value; $u\left(k\right)$ is the output of PID controller.
Artificial neural network is developed with the biological neural network, and the basic unit of biological neural network is biological neuron, while the artificial neuron is the basic unit of artificial neural networks. Neural network controllers have been expressed in many forms, a frequent representation is a multilayer feedforward network. In neural network representation, it can be easy to visualize and analyze the signal flow though the incremental PID controllers. A simple neural network control system prototype is shown in Fig. 10. There are two inputs (roll angle and vertical acceleration of vehicle), two input membership functions, and six outputs (proportional, integral and differential coefficients of PID controllers).
Fig. 10. Neural network control system prototype
Every node in the input layer, which is a membership function, is an adaptive node. There are six input nodes in this layer. ${A}_{1}=e\left(\phi \right)$ is the deviation of the ideal and the actual heeling angle; ${A}_{2}=de\left(\phi \right)/dt$ reflects the rate of change of the deviation; ${A}_{3}=$1 means the biasing logic. ${B}_{1}=e\left(a\right)$ is the deviation of the ideal and the actual acceleration of unsprung mass; ${A}_{2}=de\left(a\right)/dt$ reflects the rate of change of the deviation; ${A}_{3}=$1 means the biasing logic. The output of this node is a matching degree of an input to the corresponding membership functions in the PID set:
where, ${o}_{1,j}$ is the membership grade of a PID set.
On the hidden layer, each node represents an “AND” operator, and is a fixed node whose output is the product of the entire incoming signals. The inputs and outputs of hidden layer can be expressed as:
where, ${\theta}_{i}$ is the valve value of the hidden layer. ${w}_{ij}$ is the connected weights between the input layer and the hidden layer. $g\left(\right)$ is Sigmoid type active function, that can be expressed as:
where, ${\theta}_{1}$ is the parameter represents the biasing value. ${\theta}_{0}$ is used for regulating the shape of the Sigmoid function.
The input of the output layer is:
where, ${w}_{li}$ is the connection weights of the hidden layer neuron $i$ and the output layer neuron $l$; ${\theta}_{l}$ is the valve value of the output layer neuron.
The outputs of the output layer correspond to the three parameters of PID controller, ${K}_{p}$, ${K}_{i}$ and ${K}_{d}$. The basic principle of learning algorithm is the gradient steepest descent method, that is, adjust the weights of the network to make the error of network to a minimum. Select the performance index function as follows:
According to gradient method, the correction formula of the connection weights of output layer neurons is:
where, $\eta $ is the learning rate, $\eta >$0. The connection weights of output layer can be expressed as:
According to gradient method, the correction formula of the connection weights of hidden layer neurons is:
The connection weights of hidden layer can be expressed as:
The formula of the connection weights of output and hidden layer neurons $l$ and $i$ under the training of sample $p$ can be expressed as:
4. Simulation and discussion
The virtual prototyping models of front and rear suspension are established in SIMPACK software. MR damper is an important part of this model, and the damping force can be controlled by semiactive suspension. The Magic Formula tire force model is used in SIMPACK. The outs of the tire are obtained from the lookuptable provided by tire tests. The light vehicle model including driver model, suspension and steering systems is shown in Fig. 11.
To simulate the ride comfort and handing stability of the light bus, a neural network controller and virtual prototyping model of automobile are established respectively in MATLAB and SIMPACK software. By defining the data exchange interface between the integrated controller and the virtual prototyping model, the control of semiactive suspensions is achieved in simulation environment. The SIMPACK software, contains various kinds of random roads, has been used in simulation tests. And the interface standard provided by SIMPACK/SIMAT makes the cosimulation with MATLAB software. The four MR dampers with adjustable damping coefficient are controlled by the neural network PID control algorithm. The communication interface between the software makes the connection between physical models of the dampers and controller possible. In the cooperating applications, there are three signal transmission paths. In the first path, the physical model of fullcar passes the values required by neural network controller through the SIMAT module of SIMPACK. Then, the neural network controller analyzes the input signals and makes an output signal to control the intensity of currents. The magnetic intensity of MR damper model is controlled by the input currents, in the end, the forces of dampers reactions outputted by MR damper model is entered into the semiactive suspension of the fullvehicle model, as shown in Fig. 12.
Fig. 11. Dynamic multibody model of vehicle
Fig. 12. Block diagram of MATLAB and SIMPACK cosimulation
In order to verify the control performances of neural networkPID control strategy, the passive and semiactive suspension vehicle are compared. In ride comfort simulation conditions, the speed is kept constant as 80 km/h driving on grade B road. In this situation, the vertical, pitch and roll acceleration of vehicle body are tested to verify the control performance of semiactive suspensions.
The power spectral density (PSD) comparison of the vertical, pitch and roll acceleration curves are plotted in Fig. 13. We can find that the PSD amplitude of semiactive suspensions is much smaller than passive suspension. It clearly indicates that the semiactive suspensions controlled by neural network controller could be better coordinate the ride comfort and handling stability compared with the passive suspension.
In order to evaluate the handling performance of the control algorithm, a sine wave test is carried out. The vehicle speed is kept constant as 50 km/h, and the roll angle and yaw angular velocity of vehicle body are tested to verify the effect of semiactive suspensions (see Fig. 14, red line – time domain curve of passive suspension, black line  time domain curve of semiactive suspension). As seen from the simulation results, with the application of semiactive suspensions, the magnitude of heeling angle and yaw angular velocity of vehicle body has been effectively constrained.
5. Test results analysis
The prototype vehicle is equipped with four MR dampers. The control systems include four acceleration sensors, A/D and D/A boards, dSPACE AutoBox control system with the control software, four MR dampers and four electronic current drivers which provide control current and a data acquisition instrument. The acceleration sensors are used to measure the vertical accelerations of the front and rear sprung masses. The location of acceleration sensors placements are shown in Fig. 15.
Fig. 13. Time history and PSD analysis results of vehicle body vibration (dotted line –passive suspension PSD results, continuous line – semiactive suspension PSD results)
a) Vehicle body vertical acceleration
b) Vehicle body pitch acceleration
c) Vehicle body roll acceleration
Fig. 14. Time history of snaking motion simulation (dotted line – passive; continuous line – semiactive)
a) Heeling angle
b) Yaw angular velocity
Fig. 15. Location of acceleration sensors
Fig. 16. Vertical acceleration test results ((dotted line – passive; continuous line – semiactive)
a) PSD comparsion of driver seat vertial acceleration
b) PSD comparsion of vehicle body vertial acceleration
Fig. 17. Handing stability test
The vehicle speed is kept constant as 80 km/h driving on a random highway road, and the vertical acceleration comparison of driver seat and vehicle body center are shown in Fig. 16.
The handing stability testing of the light vehicle with semiactive suspensions is shown in Fig. 17. In addition to the testing instruments used in the ride comfort experiment, an equipment called Gyro instrument was used to test the roll and yaw angular velocity of vehicle body. In the handing stability test, six stakes were layout on the test site (30 m spacing), and the test vehicle passed through all stake at the constant speed of /h.
The heeling angle and yaw angular velocity of vehicle body are shown in Fig. 18 (dotted line – passive suspension simulation results, continuous line – semiactive suspension simulation results)). As seen from the Fig. 18, the semiactive controller applied on MR dampers effectively reduces the roll and yaw angle of vehicle body, and the effect of semiactive suspension has been reflected in the handing stability test.
Fig. 18. Experimental results of handing stability ((dotted line – passive; continuous line – semiactive)
a) Heeling angle
b) Yaw angular velocity
6. Conclusion
A fullvehicle semiactive suspension control system equipped with four MR dampers has been developed in the paper. The Neural networkPID control algorithm is proposed for the semiactive suspension systems. The numerical simulation and experiments test results demonstrated that semiactive suspension can effectively attenuate the vertical vibration and the pitch, roll motion of vehicle body. It can be concluded from the obtained results that the semiactive suspension with MR damper can improve the ride comfort and handling performance of automobile simultaneously.
Acknowledgements
The project is supported by the National SciTech Support Plan (Grant No. 2014BAG02B02).
References
 Tseng H. Eric, Davor Hrovat State of the art survey: active and semiactive suspension control, Vehicle System Dynamics, Vol. 53, Issue 7, 2015, p. 10341062. [Search CrossRef]
 Hong K. S., Sohn H. C., Hedrick J. K. Modified skyhook control of semiactive suspensions: a new model, gain scheduling, and hardwareintheloop tuning. Journal of Dynamic Systems, Measurement, and Control, Vol. 124, Issue 1, 2012, p. 158167. [Search CrossRef]
 Ahmadian M. E. H. D. I., David E. Simon An analytical and experimental evaluation of magneto rheological suspensions for heavy trucks. Vehicle System Dynamics, Vol. 37, 2003, p. 3849. [Search CrossRef]
 Lajqi Shpetim, Stanislav Pehan Designs and optimizations of active and semiactive nonlinear suspension systems for a terrain vehicle. Strojniški vestnik – Journal of Mechanical Engineering, Vol. 58, Issue 12, 2012, p. 732743. [Search CrossRef]
 Tchamna Rodrigue, Edward Youn, Iljoong Youn Combined control effects of brake and active suspension control on the global safety of a fullcar nonlinear model. Vehicle System Dynamics, Vol. 52, Issue 1, 2014, p. 6991. [Search CrossRef]
 Assadsangabi B., Eghtesad M., Daneshmand F., et al. Hybrid sliding mode control of semiactive suspension systems. Smart Materials and Structures, Vol. 18, Issue 12, 2009, p. 110. [Search CrossRef]
 Kashem Saad Bin Abul, Mehran Ektesabi, Romesh Nagarajah Comparison between different sets of suspension parameters and introduction of new modified skyhook control strategy incorporating varying road condition. Vehicle System Dynamics, Vol. 50, Issue 7, 2012, p. 11731190. [Search CrossRef]
 Karnopp D., Crosby M. J., Harwood R. A. Vibration control using semiactive force generators. Journal of Manufacturing Science and Engineering, Vol. 96, Issue 2, 1974, p. 619626. [Search CrossRef]
 Song X., Ahmadian M., Southward S., et al. An adaptive semiactive control algorithm for magnetorheological suspension systems. Journal of Vibration and Acoustics, Vol. 127, Issue 5, 2005, p. 493502. [Search CrossRef]
 Moradi Morteza, Afef Fekih Adaptive PIDsliding mode fault tolerant control approach for vehicle suspension systems subject to actuator faults. IEEE Transactions on Vehicular Technology, Vol. 63, Issue 3, 2014, p. 10411054. [Search CrossRef]
 Sun L., Cai X., Yang J. Genetic algorithmbased optimum vehicle suspension design using minimum dynamic pavement load as a design criterion. Journal of Sound and Vibration, Vol. 301, Issue 1, 2007, p. 1827. [Search CrossRef]
 Chiang H., Lee L. Optimized virtual model reference control for ride and handling performanceoriented semiactive suspension systems. IEEE Transactions on Vehicular Technology, Vol. 64, Issue 5, 2015, p. 16791690. [Search CrossRef]
 Cong S., Liang Y. PIDlike neural network nonlinear adaptive control for uncertain multivariable motion control systems. IEEE Transactions on Industrial Electronics, Vol. 56, Issue 10, 2009, p. 38723879. [Search CrossRef]
 Hahm D., Koh H.M., Ok S.Y., Park W., Chung C., Park K.S. Costeffectiveness evaluation of MR damper system for cablestayed bridges under earthquake excitation. Proceedings of the 3rd International Conference on Bridge Maintenance, Safety and Management – Bridge Maintenance, Safety, Management, LifeCycle Performance and Cost, 2006, p. 301302. [Search CrossRef]
 Guo A. X., Cui L. L., Li H. Structural control of seismically induced pounding of elevated bridges by using magnetorheological dampers. Proceedings of the 3rd International Conference on Bridge Maintenance, Safety and Management – Bridge Maintenance, Safety, Management, LifeCycle Performance and Cost, 2006, p. 685686. [Search CrossRef]
 Ramli R., Pownall M., Levesley M., Crolla D. A. Dynamic analysis of semiactive suspension systems using a cosimulation approach. MultiBody Dynamics: Monitoring and Simulation TechniquesIII, Vol. 32, Issue 4, 2004, p. 391399. [Search CrossRef]
 Lee S. H., Hwang Y. S. A study on a scenario using the PID method. Progress in Nuclear Energy, Vol. 51, Issue 2, 2009, p. 253257. [Search CrossRef]
 Hong S. R., Wereley N. M., Choi Y. T., Choi S. B. Analytical and experimental validation of a nondimensional Bingham model for mixedmode magnetorheological dampers. Journal of Sound and Vibration, Vol. 312, Issue 3, 2008, p. 399417. [Search CrossRef]
 Wen Y. K. Method for random vibration of hysteretic systems. Journal Engineering Mechanics, Vol. 102, Issue 2, 1976, p. 249263. [Search CrossRef]
 Spencer B. F., Jr. Dyke S. J., Sain M. K., Carlson J. D. Phenomenological model of magnetorheological dampers. Journal of Engineering Mechanics, Vol. 123, Issue 2, 1996, p. 230238. [Search CrossRef]