Open access peer-reviewed chapter

Control Strategies for Vehicle Suspension System Featuring Magnetorheological (MR) Damper

By Min-Sang Seong, Seung-Bok Choi and Kum-Gil Sung

Submitted: December 8th 2010Reviewed: April 28th 2011Published: September 6th 2011

DOI: 10.5772/24556

Downloaded: 5660

1. Introduction

Vehicle suspension is used to attenuate unwanted vibrations from various road conditions. So far, three types of suspension system have been proposed and successfully implemented; passive, active and semiactive. Though the passive suspension system featuring oil damper provides design simplicity and cost-effectiveness, performance limitations are inevitable due to the lack of damping force controllability. On the other hand, the active suspension system can provides high control performance in wide frequency range. However, this type may require high power sources, many sensors and complex actuators such as servovalves. Consequently, one way to resolve these requirements of the active suspension system is to adopt the semiactive suspension system. The semiactive suspension system offers a desirable performance generally enhanced in the active mode without requiring large power sources and expensive hardware.

One of very attractive and effective semiactive vehicle suspension systems is to utilize magnetorheological (MR) fluid. MR fluids are currently being studied and implemented as actuating fluids for valve systems, shock absorbers, engine mounts, haptic systems, structure damper, and other control systems. The rheological properties of MR fluids are reversibly and instantaneously changed by applying a magnetic field to the fluid domain. Recently, a very attractive and effective semi-active suspension system featuring MR fluids has been researched widely. Carlson et al., 1996 proposed a commercially available MR damper which is applicable to on-and-off-highway vehicle suspension system. They experimentally demonstrated that sufficient levels of damping force and also superior control capability of the damping force by applying control magnetic field. Spencer Jr. et al., 1997 proposed dynamic model for the prediction of damping force of a MR damper. They compared the measured damping forces with the predicted ones in time domain. Kamath et al., 1998 proposed a semi-active MR lag mode damper. They proposed dynamic model and verified its validity by comparing the predicted damping force with the measured one. Yu et al., 2006 evaluated the effective performance of the MR suspension system by road testing. Guo & Hu, 2005 proposed nonlinear stiffness model of a MR damper. They proposed nonlinear stiffness model and verified it using simulation and experiment. Du et al., 2005proposed H-infinity control algorithm for vehicle MR damper and verified its effectiveness using simulation. Shen et al., 2007 proposed load-levelling suspension with a magnetorheological damper. Pranoto et al., 2005 proposed 2DOF-type rotary MR damper and verified its efficiencies. Ok et al., 2007 proposed cable-stayed bridges using MR dampers and verified its effectiveness using semi-active fuzzy control algorithm. Choi et al., 2001 manufactured an MR damper for a passenger vehicle and presented a hysteresis model for predicting the field-dependent damping force. Hong et al., 2008 derived a nondimensional Bingham model for MR damper and verified its effectiveness through experimental investigation. Yu et al., 2009 developed human simulated intelligent control algorithm and successfully applied it to vibration control of vehicle suspension featuring MR dampers. Seong et al., 2009 proposed hysteretic compensator of MR damper. They developed nonlinear Preisach hysteresis model and hysteretic compensator and demonstrated its damping force control performance.

As is evident from the previous research work, MR damper is very effective solution for vibration control of vehicle suspension system. So in this chapter, we formulate various vibration control strategies for vibration control of MR suspension system and evaluate their control performances. In order to achieve this goal, material characteristics of MR fluid are explained. Then the MR damper for vehicle suspension system is designed, modelled and manufactured. The characteristics of manufactured MR damper are experimentally evaluated. For vibration control, the quarter vehicle suspension system featuring MR damper is modelled and constructed. Then, various vibration control strategies such as skyhook control, PID control, LQG control, H control, Sliding mode control, moving sliding mode control and fuzzy moving sliding mode control are formulated. Finally, control performances of the proposed control algorithms are experimentally evaluated and compared.

2. Suspension modelling

2.1. MR fluid

Since Jacob Rabinow discovered MR fluid in the late 1940s, of which yield stress and viscosity varies in the presence of magnetic field, various applications using MR fluid have been developed such as shock absorbers, clutches, engine mounts, haptic devices and structure dampers, etc (Kim et al., 2002). Physical property changes of MR fluid are resultedfrom the chain-like structures between paramagnetic MR particles in the low permeability solvent. At the normal condition, MR fluid shows the isotropic Newtonian behavior because the MR particles move freely as shown in Fig. 1 (a). However, when the magnetic field applied to the MR fluid, MR fluid shows the anisotropic Bingham behavior and resist to flow or external shear force because the MR particles make a chain structure as shown in Fig. 1 (b). From this property, force or torque of application devices can be easily controlled by the intensity of the magnetic field.

2.2. MR damper

The schematic configuration of the cylindrical type MR damper proposed in this work is shown in Fig. 2. The MR damper is composed of the piston, cylinder and gas chamber. The floating piston between the cylinder and the gas chamber is also used in order tocompensate for the volume induced by the motion of the piston. Also the gas chamber

Figure 1.

Phenomenological behavior of MR fluid

which is filled with nitrogen gas acts as an accumulator for absorbing sudden pressure variation of lower chamber of the MR damper induced by the rapid motion of the piston. The MR damper is divided into the upper and lower chambers by piston, and it is fully filled with the MR fluid. By the motion of the piston, the MR fluid flows through the annular duct between inner and outer piston from one chamber to the other. The magnetic poles in the piston head is placed to control the yield stress of the MR fluid by supplying current to the coil. In order effectively to generate the magnetic field in the magnetic pole, the outer cylinder and both ends of inner piston are made of ferromagnetic substance, while the center of the inner piston is a paramagnetic substance. In the absence of a magnetic field, the MR damper produces a damping force caused only by fluid viscous resistance. However, if a certain level of magnetic field is supplied to the MR damper, the MR damper produces an additional damping force owing to the yield stress of the MR fluid. This damping force of the MR damper can be continuously tuned by controlling the intensity of the magnetic field.

In order to simplify the analysis of the MR damper, it is assumed that the MR fluid is incompressible and that pressure in one chamber is uniformly distributed. The pressure drops due to the geometric shape of the annular duct and the fluid inertia are assumed to be negligible. For laminar flow in the annular duct, the fluid resistance is given by (Liu et al., 2006; White, 1994)

Re=8ηLπ[Ro4Ri4(Ro2Ri2)2ln(RoRi)]E1

Figure 2.

Schematic configuration of the proposed MR damper

whereηis the viscosity of the MR fluid and Lis the length of the annular duct. RoandRiare the inner radius of the outer piston and outer radius of the inner piston respectively. By assuming that the gas does not exchange much heat with its surroundings, and hence considering its relation as adiabatic variation, the compliance of the gas chamber is obtained by

Cg=V0P0κE2

where V0and P0are the initial volume and pressure of the gas chamber respectively, and κis the specific heat ratio. On the other hand, the pressure drop due to the increment of the yield stress of the MR fluid is given by

PMR=2cLphτy(B)E3

Where c is a coefficient that depends on flow velocity profile and has a value range from 2.0 to 3.0, Lpis the length of the magnetic pole, his the gap of the annular duct, and τy(B)is the yield stress caused by the magnetic flux densityB. Therefore, the damping force of the proposed MR damper can be written as

FD=kexp+cex˙p+FMRE4

where

ParameterValue
Duct Length (L)82mm
Piston Area (Ap )1661.90mm2
Piston Rod Area (Ar )380.13mm2
Duct Width (b)123.53mm
Maximum Stroke164mm

Table 1.

Design parameters of the MR damper

ke=Ar2Cg,ce=(ApAr)2ReFMR=(ApAr)PMRsgn(x˙p)E5

wherexpand x˙pare the piston displacement and velocity respectively, Apand Arrepresent the piston and piston rod areas respectively, and sgn()is a signum function.

Figure 3.

Photograph of the manufactured MR damper

The photograph of the manufactured MR damper with optimally obtained design parameters are shown in Fig. 3. The principal design parameters of the manufactured MR damper, which can be applied to a mid-sized commercial passenger vehicle, are presented in Table 1. Fig. 4 presents the measured and analysed damping force F D characteristics of the MR damper with respect to the piston velocity at various magnetic fields. This is obtained by calculating the maximum damping force at each velocity. The piston velocity is changed by increasing the excitation frequency from 0.5 to 4.0Hz, while the excitation amplitude is kept constant ±20mm. This type of plot is frequently used in the damper manufacturing industry to evaluate the level of damping performance. It is clearly observed that the damping force is increased as the magnetic field increases, as expected. The damping force response of the MR damper is measured as shown in Fig. 5. It can be found that the time constant is about 23ms, which is obtained by inspecting the required time when the damping force reaches 63.2percent of its final steady state value.

Figure 4.

Field-dependent damping force of the MR damper

Figure 5.

Damping force responses of the MR damper

2.3. Quarter vehicle MR suspension

In order to investigate the effectiveness of vibration control algorithm, a quarter vehicle MR suspension system is constructed as shown in Fig. 6. It shows that the proposed quarter vehicle suspension model with the MR damper has two degrees of freedom. The spring for the suspension is assumed to be linear and the tire is also modelled as linear spring component. From the mechanical model, dynamic equation of system considering time constant can be expressed as follows:

Figure 6.

Mechanical model of a quarter vehicle MR suspension system

msz¨s=kszscsz˙s+kszu+csz˙uFMRmuz¨u=kszs+csz˙s(kt+ks)zucsz˙u+FMR+ktzrF˙MR=1τMRDFMR+1τMRDuE6
ParameterValue
Sprung Mass (ms)373.5kg
Unsprung Mass (mu)40kg
Stiffness Coefficient (ks)27358N/m
Tire Stiffness Coefficient (kt)211625N/m
Damping Coefficient (cs)570Ns/m

Table 2.

System parameters of the quarter vehicle suspension system

where τMRDis time constant of the MR damper, ksis the total stiffness coefficient of the suspension, including the effective stiffness keof the MR damper in equation (4), csis the damping coefficient of the suspension and is assumed to be equal toce, ktis the vertical stiffness of the tire, andzs, zuand zrare the vertical displacements of sprung mass, unsprung mass and road excitation respectively. The state space equation of proposed quarter vehicle suspension can be expressed using dynamic equation (Lee et al., 2011):

x˙=Ax+Bu+Lzr,y=CxE7

where

x=[zsz˙szuz˙uFMR]TA=[01000ksmscsmsksmscsms1ms00010ksmucsmu(kt+ks)mucsmu1mu00001τMRD]B=[00001τMRD]T,L=[000ktmu0]TC=[10000]E8

The system parameters of the quarter vehicle MR suspension system are chosen on the basis of the conventional suspension system for a mid-sized passenger vehicle, and listed in Table 2.

3. Control strategies

In order to evaluate vibration control performance of the quarter vehicle MR suspension system, various control strategies are formulated and experimentally implemented.

3.1. Skyhook controller

Skyhook controller is simple but very effective control algorithm. It is well known that the logic of the skyhook controller is easy to implement in the real field. The principle of skyhook control is to design the active or semiactive suspension control so that the sprungmass is linked to the sky in order to reduce the vertical oscillations of the sprung mass. Fig. 7 shows the conceptual scheme of skyhook controller for vehicle suspension system. The desired damping force is set by

u=Cskyz˙sE9

Figure 7.

Scheme of the skyhook controller for vehicle suspension system

whereCskyis the control gain, which physically indicates the damping coefficient. In this work, the value is chosen as 2470 using trial-and-error method.

3.2. LQG controller

Optimal control of a linear system with respect to the quadratic objective function under incomplete measurements corrupted by white Gaussian noise is generally referred to as the linear quadratic Gaussian (LQG) problem. The optimal control is a linear function of the state estimates obtained from the Kalman-Bucy filter. The LQR (linear quadratic regulation) is a state feedback problem, whereas LQG is an output feedback problem, which is more realistic (Bahram& Michael, 1993). Fig. 8 shows the block diagram of an LQG controller. In this study, the control input is formulated as follows:

u=KLQGx^E10

wherex^is the estimated state. Control gain KLQGis set as [1248.31150.34121.815.50.4]in this work.

3.3. Hcontroller

In reality, the sprung mass of the vehicle is varied by the loading conditions such as the number of riding persons and payload. And it makes the pitch and roll mass moment of inertia to be changed. Therefore, in order to successfully control the vibration, a robust control algorithm is required by considering the parameter variations of the system. From the structured suspension model, a desirable damping force required to effectively suppress the vibration is determined by adopting Hcontroller. Fig. 9 shows the nominal plant, perturbed plant, shaped plant and loop gain graph. As shown in Fig., loop gain and perturbed plant are well matched and it means loop perturbation is successfully constructed.

For reducing the steady state error and suppressing the effectiveness of disturbance, weighting function can be designed as follows (Choi et al., 1999; 2002):

W=180s+88sE11

For the design of controller, design index γis set as 2.99. Fig. 10 shows the sensitivity and complementary sensitivity of the closed loop system. Therefore, control gain can be calculated as follows:

u=KHx,KH=N(s)D(s)E12

Where

N(s)=1.45e2s62.29e4s51.82e6s41.12e8s32.10e9s2+5.20e10s6.91e6,D(s)=s7+92.03s6+8.34e3s5+3.83e5s4+3.47e6s3+3.48e7s2+7.51e2sE13

3.4. Sliding mode controller (SMC)

Fig. 11 shows the conceptual scheme of sliding mode control algorithm. After the initial reaching phase, the system states slides along the sliding surface. The first step to formulate the SMC is to design a stable sliding surface. The stable sliding surface for the control system is defined as follows:

s=c1x1+c2x2+c3x3+c4x4+x5E14

Figure 8.

Block diagram of LQG controller

Figure 9.

Singular value plots of the quarter vehicle MR suspension system

Figure 10.

Frequency domain indicators of the quarter vehicle MR suspension system

Figure 11.

Scheme of the sliding mode control

where ci(i=1,2,3,4)are sliding surface coefficients to be determined so that the sliding surface is stable. Then the sliding mode controller can be formulated which satisfies sliding mode condition ss˙0as follows (Choi et al., 2000):

u=τMRDc5(c1x2+c2ms(ksx1csx2+ksx3+csx4x5)+c3x4+c4mu(ksx1+csx2(kt+ks)x3csx4+x5)+1τMRDx5)ksgn(s),(k0)E15

Where k stands for the discontinuous gain which is a positive number. The discrete gain can be changed to continuous gain for reducing the chattering problem.

sat(s)={s/εfor|s/ε|1sgn(s)for|s/ε|1E16

3.5. Moving sliding mode controller (MSMC)

It is required to reduce reaching phase to improve the performance and robustness of the SMC. Reaching time can be successfully reduced in case that sliding surface is rotated orshifted by considering reaching phase. Fig. 12 illustrates two moving patterns of the sliding surface: rotating and shifting. The sliding surface for moving sliding mode controller can be determined as follows:

s=c1(x1,x2)x1+c2x2+c3x3+c4x4+x5+α(x1,x2)c1(x1,x2)=(Δfrx2)x1,α(x1,x2)=c1ox1+x2ΔfsE17

whereΔfrand Δfsare gain for surface rotating and surface shifting respectively. c1ois initial value ofc1. The moving sliding mode controller can be formulated as follows:

u=τMRDc5(c1(x1,x2)x2+c2ms(ksx1csx2+ksx3+csx4x5)+c3x4+c4mu(ksx1+csx2(kt+ks)x3csx4+x5)+1τMRDx5+α(x1,x2))ksgn(s),(k0)E18

3.6. Fuzzy moving sliding mode controller (FMSMC)

The fuzzy moving sliding mode controller, which can change the coefficients and intercepts of sliding surface by fuzzy tuning which takes into account for location of reaching phase, is developed. Fig. 13 presents the block diagram of the proposed FMSMC. The basic configuration of fuzzy control consists of three components: a fuzzification interface, a decision-making logic and a deffuzification interface. In this study, the coefficient of sprung mass displacement x1is considered with priority while other coefficients of sliding surface are fixed. Rotating algorithm, which can change the sliding surface c1as a function ofx1, is applied to reduce the reaching phase. However, it is hard to expect that the robustness is enhanced in case that initial condition of reaching phase is located in quadrant 1 or quadrant 3. Hence, shifting algorithm, which can change the intercept αas a function ofx1, is adopted and determined according to the relative place of sliding surface, and this will sufficiently reduce the reaching phase. Therefore, the sliding surface can be written as follows (Cho et al., 2007):

s=c1(x1)x1+c2x2+c3x3+c4x4+α(x1)E19

Figure 12.

Scheme of the moving sliding mode control

Figure 13.

Block-diagram of the proposed FMSMC

where c1and αare determined by tuning of fuzzy logic which takes into account for displacement and direction of sprung mass displacement ofx1. In this case, the sliding surface should maintain stability, although c1and αare changed according to time. In order to reduce reaching time in the rotating algorithm, c1should be high in case of the opposite case. In the letter case, the sprung displacement x1can be converted toward 0 more quickly. In the shifting control algorithm, absolute value of αshould be high in case that absolute value of x1is large and vice versa. Furthermore, sliding surface is moved upward for positive x1and sliding surface is moved downward for the opposite case. Consequently, we can formulate the following control input of the proposed FMSMC, which is combined rotating and shifting algorithm:

u=τMRDc5{(c1(x1)+Δc1χ1rsinh(χ1r)cosh2(χ1r)(αmaxσs)cosh2(χ1s))x2+c2ms(ksx1csx2+ksx3+csx4x5)+c3x4+c4mu(ksx1+csx2(kt+ks)x3csx4+x5)1τMRDx5}ksgn(s),(k0)E20

The control input determined from the FMSMC is to be applied to the MR damper depending upon the motion of suspension travel. The detailed control algorithm was described by Cho et al., 2007.

3.7. Semi-active condition

The control input udirectly represents the damping force ofFMR. On the other hand, the damping of the suspension system needs to be controlled depending upon the motion of suspension travel. Therefore, the following semi-active actuating condition is imposed:

u={uforu(z˙sz˙u)00foru(z˙sz˙u)0E21

This semi-active condition can assure the increment of damping characteristic of the suspension system and hence increase the stability of the system.

4.Control performances

Vibration control performances of the quarter vehicle MR suspension system are evaluatedunder two types of excitation (road) conditions. The first excitation, normally used to reveal the transient response characteristic, is a bump described by

zr=Am[1cos(ωt)]E22

where ω=2πf,f=1/TandT=D/V. Am(=0.07m)is the bump height, D(=0.8m)is the width of the bump and Vis the vehicle velocity. In the bump excitation, the vehicle travels the bump with a constant vehicle velocity of 3.08km/h(=0.856m/s). The second excitation is a random road conditions.

Fig. 14 shows the experimental apparatus of the quarter vehicle MR suspension system to evaluate and compare the effectiveness of vibration isolation of the proposed control algorithms. The MR suspension (assembly of the MR damper and spring), sprung mass and tire are installed on the hydraulic excitation system. The sprung mass displacement and excitation displacement are measured by LVDT (linear variable differential transformer) and the suspension travel is measured by a wire sensor. The hydraulic system applies the road profile to the MR suspension system. The current amplifier applies the control current determined from the control algorithm to the MR damper.

Figure 14.

Experimental apparatus of the quarter vehicle MR suspension system

Fig. 15 presents vertical displacement, vertical acceleration of sprung mass and suspension travel versus time responses of the MR suspension system for the bump excitation. As shown in the Fig., it can obviously be found that unwanted vibrations induced from the bump excitation are well suppressed by adopting the control algorithms in the MR suspension system. Fig. 16 shows the performance comparison of vertical acceleration RMS of the quarter vehicle MR suspension system for the bump excitation. The RMS value is calculated from:

RMS=x12+x22++xn2nE23

As shown in Fig., skyhook control algorithm shows quite good vibration control performance. Also, Fig. 17 shows the PSD (power spectrum density) of sprung mass acceleration in frequency domain for the random excitation and Fig. 18 shows the performance comparison of vertical acceleration RMS of the quarter vehicle MR suspension system for the random excitation. As expected, the PSD of the vertical acceleration has been

Figure 15.

Bump responses of the quarter vehicle MR suspension system

Figure 16.

Performance comparison of vertical acceleration RMS of the quarter vehicle MR suspension system (bump)

Figure 17.

Random responses of the quarter vehicle MR suspension system (72km/h)

Figure 18.

Performance comparison of vertical acceleration RMS of the quarter vehicle MR suspension system (random, 72km/h)

considerably reduced in the neighbourhood of body resonance (1~2Hz) by applying control input. The control results presented in Figs. 15~18 indicate that ride comfort of a vehicle system can be substantially improved by employing the MR suspension system associated with the proposed control algorithms.

5. Conclusion

In this chapter, control algorithms for vibration control of the quarter vehicle MR suspension system were proposed and its effectiveness was experimentally verified and compared. In order to achieve this goal, a cylindrical MR damper was designed and manufactured. After evaluating the field-dependent damping characteristics of the MR damper, a quarter vehicle suspension system was then constructed and its governing equations of motion were derived. In order to obtain a favourable control performance of the MR suspension system, skyhook controller, LQG controller, H controller, sliding mode controller, moving sliding mode controller and fuzzy moving sliding mode controller were designed and experimentally realized to the quarter vehicle MR suspension system. It has been experimentally shown that the proposed control algorithms can reduce the unwanted vibration under bump excitation. In addition, it has been demonstrated that using the proposed control methodologies vibration levels such as sprung mass acceleration under random excitation can be significantly reduced at body resonance region. The results presented here are self-explanatory justifying that the control strategies implemented in this work can produce very similar vibration control performance. This directly indicates that ride comfort of the vehicle system can be substantially improved by adopting the MR suspension system associated with the proposed control strategies.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Min-Sang Seong, Seung-Bok Choi and Kum-Gil Sung (September 6th 2011). Control Strategies for Vehicle Suspension System Featuring Magnetorheological (MR) Damper, Vibration Analysis and Control - New Trends and Developments, Francisco Beltran-Carbajal, IntechOpen, DOI: 10.5772/24556. Available from:

chapter statistics

5660total chapter downloads

4Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

A Semiactive Vibration Control Design for Suspension Systems with Mr Dampers

By Hamid Reza Karimi

Related Book

First chapter

Passive, Adaptive, Active Vibration Control, and Integrated Approaches

By Dirk Mayer and Sven Herold

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us