Open access peer-reviewed chapter

Design of Robust Digital Pole Placer for Car Active Suspension with Input Constraint

By Hisham M. Soliman, Rashid Al-Abri and Mohammed Albadi

Submitted: June 20th 2017Reviewed: August 16th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.70587

Downloaded: 230

Abstract

This chapter deals with the problem of state feedback control for an active quarter-car suspension system with control input constraint. The dynamics of the suspension system is first formed in terms of the control objectives: ride comfort, suspension deflection, and maximum actuator control force. The control task is formulated as robustly placing the closed poles in a desired region against different passenger load. Since digital computers are widely used in the vehicle industry, a new saturated controller design method is presented for regional pole-placement of uncertain discrete time systems. The constraint of control input saturation is considered in the design phase. The desired dynamic performance for uncertain discrete-time systems is represented by the settling time and damping ratio. A sufficient condition is derived to place the poles in a desired region. The design is formulated in terms of linear matrix inequality optimization. The effectiveness of the proposed design is illustrated by applying it to a quarter?car active suspension system. Different road tests for the proposed controller are carried out: step and bump disturbances. The proposed design achieves the desired oscillation damping due to road disturbances in addition to passenger comfort. The results are compared with the passive suspension system.

Keywords

  • constrained control
  • regional pole placement
  • linear matrix inequality (LMI) optimization
  • discrete-time systems
  • vehicle active suspension

1. Introduction

The design of a car suspension system is of paramount importance for improving the ride comfort, maintaining vehicle maneuverability, and retaining the safety of passengers [1]. Vehicle suspension systems include wishbone, spring, and shock absorber (e.g., damper) to transmit and mitigate forces between the car body and the road. This contributes to the passenger comfort and the ride. However, the widely used passive suspension systems which use springs and dampers with fixed parameters (stiffness and damper coefficients) do not achieve satisfactory suspension performance under wider driving maneuvers. This problem is alleviated by the recent developments of semi-active suspension and active suspension [2].

Semi-active suspensions use the springs and dampers whose parameters (e.g., stiffness and damper coefficients) are adjusted corresponding to different vehicle driving scenarios [3]. These varying parameters provide considerable improvements over passive suspension systems. However, special springs and/or dampers (e.g., a magneto rheological fluid, MR, damper) are used in this case, and the construction and tuning of such components are not easy. Moreover, using such elements can lead to potential difficulties, for example, MR dampers exhibit hysteresis dynamics, which are highly nonlinear and difficult to model.

Active suspension has attracted extensive research because in this framework an extra actuator is placed between the car body and the wheel-axle. It is installed in parallel to other suspension mechanisms (e.g., damper and spring [4]). So, this actuator can be controlled to dissipate energy from the road disturbances and thus reduces the impact on the displacement of vehicles. It is well known that active suspension requires high energy demand and cost. Therefore, active suspensions have not been widely used in commercial vehicles. However, it is expected that active suspension techniques will be adopted by the industry due to its potentials to improve the suspension performance.

Extensive research work on active suspension has been carried out in the past few decades [5, 6]. A crucial issue in the active suspension system designs is the control strategies, which should make the actuators pull down or push up the suspension motions and also achieve other suspension requirements. A model-reaching adaptive control is presented in [7] to achieve the ideal isolation of a skyhook target; the vehicle maintains a stable posture. An LPV gain-scheduling controller is proposed in [8] for a quarter-vehicle active suspension system. To manage the trade-off between the conflicting performances, the H∞ control method can manage the trade-off and obtain a compromise performance [9, 10]. An H∞ control is designed for active suspension systems that are subject to actuator time delay which is given in [11]. Sampled-data H∞ control of uncertain active suspension systems via fuzzy control is presented in [12]. Adaptive sliding-mode control for active suspensions using fuzzy approach is shown in [13]. In [14], a linear-quadratic-Gaussian (LQG) control is used to obtain a trade-off between the conflicting suspension requirements.

Many physical systems are inherently nonlinear and subject to variation in the operating point. To overcome such difficulties, the system to be controlled is represented by an uncertain linear time-invariant model. The uncertainty can be cast into either polytopic or norm-bounded form. The powerful, robust control techniques of liner systems can then be applied [15]. The poles of systems without uncertainty can be placed in desired locations so as to achieve good dynamic behavior in terms of settling time and damping ratio [16]. However, for systems with uncertainties, the closed-loop poles can be assigned to a domain (D) or region, rather than specific locations [17]. This is termed D-stability or regional pole placement by using robust controllers against system uncertainty. Regional pole placement for continuous-time systems with polytopic uncertainty using state feedback is presented in [18], while output feedback is presented in [19]. Regional pole placement with guaranteed-cost control of active suspension system of a quarter-car model is given in [20].

In many practical control problems, the actuator has limited output, called saturated (or constrained) control [21]. Combining robust pole placement with saturated control is termed multi-objective control. When the actuator saturation is not considered in the design phase, the performance of the designed control system seriously deteriorates. A state-feedback control achieving guaranteed-cost regional pole placement with the control limits, i.e., the actuator force must not exceed a certain limit, is presented in [22] for continuous-time case.

This chapter is concerned with the digital state-feedback controller design problem for active quarter-car suspension systems. The main contributions of this chapter can be summarized as follows: (1) formulate multi-object control constraints in LMIs via state-feedback control approach for the suspension vehicle system and (2) construct a simple new-type state-feedback controller. The proposed design achieves four constraints: (1) damps effectively the car oscillations due to disturbance of road irregularities, (2) robustness against different car loads, (3) no control limits violation, and (4) optimal passenger comfort against road roughness. It is expected that the proposed research pave the way for implementing the theoretical findings into vehicle suspension systems.

Notations: Capital and small letters denote matrices and vectors, respectively. I and 0 denote the identity matrix and zero matrix, respectively. W′ and W−1 denote the transpose and the inverse of any square matrix W, respectively. W > 0 (W < 0) denotes a symmetric positive (negative)-definite matrix W. The symbol • is as an ellipsis for terms in matrix expressions that are induced by symmetry. For example,

L+W+N+W'+N'NN'M=L+W+N+NME1

Facts: the following facts [15] are used in the sequel:

Fact 1 (congruence transformation):

The definiteness of a matrix W does not change under the congruence transformation H′WH.

Fact 2:

For any real matrices W1, W2 and ∆k with appropriate dimensions and ∆k′∆k ≤ I, ↔‖Δk‖ ≤ 1, it follows that

W1ΔkW2+W2'Δk'W1'ε1W1W1'+εW2'W2,ε>0E2

where ∆k represents system uncertainty at discrete time k. The use of this lemma is to eliminate the uncertainty.

Fact 3 (Schur complement):

This fact is useful in transforming a nonlinear matrix inequality into a linear one.

For constant matrices W1, W2, and W3, where W1′ = W1 and 0 < W2 = W2′, it follows that

W1+W'3W21W3<0W1W3'W3W2<0E3

2. Problem formulation

The quarter-vehicle suspension model is shown in Figure 1 [9].

Figure 1.

Quarter-car model with an active suspension.

The dynamics of the quarter-car model can be described by the state equation as follows [9]:

ẋ=Ax+Bu+Dww,x0=x0E4

where

A=0101ksMscsMs0csMs0001ksmcsmktmcsm,Dw=0010,B=[0usMs0usm]E5

where ks and cs are the parameters of the so-called passive suspension; kt stands for the tire stiffness; and Ms and m represent sprung (body) and unsprung (tire) masses, respectively. Moreover, xs − xu is the suspension stroke, xu − x0 is the tire deflection, and x0 is the vertical ground displacement caused by road unevenness. The state variables of the model (4) are defined as x1 = xs − xu, x2=ẋs, x3 = xu − x0, and x4=ẋu. The disturbance due to the road roughness is ω=ẋ0. The normalized active force, u = uf/us, is the control input, and uf is the active force generated by a hydraulic actuator.

The suspension performance, such as ride comfort, suspension deflection, and road holding, are considered as the control design objectives. Since the ride comfort can be quantified by the body acceleration in the vertical direction, it is reasonable to choose body acceleration as the regulated output, z(t). The less the vertical body acceleration, the more comfort ride results [9]. One of the main objectives of the proposed controller is to minimize the vertical acceleration z(t) to improve vehicle ride comfort. The body acceleration is given by [9]

z=Cx+B1uE6

where

C=1Mskscs0cs,B1=usMsE7

In addition, due to the mechanical structure, the suspension stroke should not exceed the allowable maximum. Or equivalently, the active control force provided for the active suspension system should be confined to a certain range. Due to actuator saturation, it is in addition assumed that the normalized active force is bounded as |u(t)| ≤ 1, Figure 2. Therefore, the bound on the control input is umax = 1, and the active force is bounded by us.

Figure 2.

Feedback-saturated control systems.

Since microprocessors are widely used in car industries, system (4) can be discretized to design a discrete-time controller. Considering the disturbance effect later, the resulting discrete-time uncertain system with saturated control input is

xk+1=A+ΔAkxk+B+ΔBk.satukE8

where xk and uk are the state and control vectors of dimension n and m, respectively. The pair (A,B) is assumed to be controllable. The matrices ΔAk and ΔBk are time varying which represents parametric uncertainty due to changes in passenger load and/or nonlinearities. The uncertainty is assumed to be of the norm-bounded form:

ΔAk,ΔBk==MΔkNNbE9
Δk'ΔkIΔk1E10

where M and N are known constant matrices with appropriate dimensions. Δk is an unknown matrix with Lebesgue measurable elements.

The control signal is constrained due to practical limitation. The saturated controller to develop is assumed to be state feedback and symmetric and normalized as follows:

u=Fx,1<uj<+1,j=1morsatuj=1ifuj1ujif1<uj<11ifuj1E11
uk=FxkE12

The closed-loop system is given by

xk+1=A+ΔAk+B+ΔBkFxkE13

The problem is to develop a saturated discrete-time controller, which robustly stabilizes the closed-loop system (13) and ensures a good dynamic performance, described by maximum settling time and minimum damping ratio despite the system uncertainty. To achieve both constraints of max settling time and min damping ratio, the closed-loop poles must lie in the hatched area, as shown in Figure 3. This is termed D-stability in which the poles must lie inside the region D for all admissible uncertainties. The proposed controller must achieve D-stability in addition to the constraint of maximum actuator control force.

Figure 3.

Desired region of closed-loop poles, hatched.

3. Digital saturated regional pole placer

To design the abovementioned controller, the nonlinear saturated control function is first linearized and approximated by a convex hull. For this, the following two lemmas are used [21]:

Lemma-1

For all u∈ Rm and θ ∈ Rm such that |θj| < 1 , j ∈ [1, m],

satucoDiu+Diθi1η,E14

with co denoting the convex hull.

Eq. (14) has the following equivalent form:

satu=i=1ηγiDiu+Diθ,γi0E15

Here, Di is an m × m diagonal matrix with elements either 1 or 0 and Di= (IDi), which results in η = 2m possible matrices. The matrices Di and Diare introduced to model the saturation function as a linear one. If Di is selected as I, Dibecomes 0, and the resulting controller will be unsaturated. Recall that these controllers (14) work in a linear region and do not allow saturation to occur.

The following sets are defined:

DF=xRn:1Fx1andεPρ=xRn:xTPxρρ>0E16

where P is a symmetric positive-definite matrix. The sets D(F) and ε(P, ρ) represent, respectively, a symmetrical polyhedral and an ellipsoidal one. The following result is recalled:

Lemma-2

For a given positive scalar ρ, if there exist matrices Y ∈ Rmxn and Z ∈ Rmxn and a positive-definite matrix X = X ∈ Rnxn and solutions to the following LMIs:

AX+B(DiY+DiZ]+<0,E17
1/ρZjX>0E18
i=1,,η;j=1,,m,E19

then when ΔA = 0 and ΔB = 0, the closed-loop saturated control system is asymptotically stable at the origin ∀x0 ∈ ε(P, ρ) with

F=YX1E20
P=X1E21

It is worth mentioning that LMI (17) guarantees asymptotic stability, while LMI (17) ensures that the ellipsoidal set ε (P,ρ) is contained inside the polyhedral set D(H), allowing the control to be saturated. As a special case, by selecting Di = I, the control works only in a linear region without reaching saturation. In this case, Di= 0 in LMI (18), and Zi is replaced by Yi in LMI (17) to have ε(P, ρ) ⊂  D(F).

Our control target is to design a saturated controller that maintains asymptotic stability against system uncertainty, as well as to place the closed-loop poles in a desired region D-stability, if possible, so as to achieve a good dynamic response in terms of settling time ts and damping ratio ζ.

Pole placement in the region shown in Figure 3 is difficult to solve. However, the problem can be easily solved by approximating the spiral of the constant damping ratio as a circle, as shown in Figure 4. The problem is thus reduced to placing the poles in between the two circles, one for ts and the other for ςapprox.

Figure 4.

(a) ζ circle approximation and (b) desired region of the poles, hatched.

We will consider two design cases: without and with saturated inputs.

Design case 1: unsaturated control

Consider the uncertain system with unsaturated control

xk+1=A+ΔAxk+B+ΔB.ukE22

Theorem 1

If there is a feasible solution to the following LMIs, then the closed-loop poles lie inside the two circles of center q1 and radius r1 and center q2 and radius r2:

X=X'>0,ε1>0,ε2>0,ν1>0,ν2>0,r12XAX+BY+q1XX+ε1+ν1MM'NX0ε1INbX00ν1I<0E23
r22XAX+BY+q2XX+ε2+ν2MM'NX0ε2INbX00ν2I<0E24

Moreover, the controller is given by

F=YX1E25

Proof

It is well known [17] that the eigenvalues of matrix A lie inside a circle of center −q and radius r if and only if

P=P'>0,r2PA+qIP1<0E26

or equivalently,

P=P'>0,r2PA+BF+qIP1+0MΔtN'0'++0MΔtF'Nb'0'+<0E27

The last matrix inequality is satisfied if

P=P'>0,r2PA+BF+qIP1+ε0M0M'+ε1N'0N'0'+ν0M0M'+ν1F'Nb'0F'Nb'0'<0is satisfied or

r2PA+BF+qIP1+ε+νMM'N0εINbF00νI<0E28

The last matrix inequality can be linearized by pre- and postmultiplying by [P−1, I, I, I], that is, applying fact 1 and substituting P−1 = X, FX = Y.

Note that the above condition is only a sufficient condition for regional pole placement in one circle. Since we have to achieve pole placement in the area between the two circles, one for ts and the other for ζ, we have theorem 1.

Design case 2: saturated control

Consider the uncertain system with saturated control

xk+1=A+ΔAxk+B+ΔB.satukE29

Theorem 2

The regional pole placement with robust saturated state-feedback control (12) for the uncertain system (8) can be achieved if there exist X = X′  >  0 , ε1  >  0 ,  ε2  >  0 , ν1  >  0 ,  ν2  >  0 , Y and a feasible solution to the following LMIs:

r12XAX+BDiY+D¯iZ+q1XX+ε1+ν1MM'NX0ε1INbDiY+D¯iHX00ν1I<0E30
r22XAX+BDiY+D¯iZ+q2XX+ε2+ν2MM'NX0ε2INbDiY+D¯iHX00ν2I<0E31
1/ρZjX>0E32

for i = 1,…, η; j = 1, … , m.

Moreover, the saturated robust pole placer is given by

F=YX1E33

Proof

From (26), the poles of the closed-loop uncertain system (13) lie inside the circle of center (−q, 0) and radius r if and only if there is a feasible solution to the following LMI:

P=P'>0r2PA+ΔA+B+ΔBDiF+DiH+qIP1<0E34

Eq. (34) is satisfied if the following inequality is satisfied:

P=P'>0r2PA+BDiF+DiH+qIP1+ε0M0M'+ε1N'0N0+ν0M0M'+ν1DiF+D¯iH'Nb'0NbDiF+D¯iH0<0E35

or equivalently,

r2PA+BDiF+D¯iH+qIP1+ε+νMM'N0εINbDiF+D¯iH00νI<0E36

By post- and premultiplying (36) by P−1, that is, by applying the congruence transformation and substituting P−1 = X, one gets (30). This completes the proof.

4. Design of vehicle active suspension control

In this section, we present a design example to illustrate the effectiveness of the proposed controller design method. Unlike the previous approaches, the proposed design introduces a digital computer or a microprocessor into the control loop as given in theorem 2. It is used to design an active suspension system for the quarter-vehicle model shown in Figure 1 [9].

Note that active forces generated by hydraulic actuators and considered as control inputs are bounded because of actuator saturation. The system parameters and nominal values are shown in Table 1.

ParameterValue
Ms320 kg
m40 kg
Ks18 kN/m
kt200 kN/m
cs1 kN.s/m
us1.5 KN

Table 1.

Quarter-vehicle active suspension parameters.

The normalized input is bounded as |u(t)| ≤ 1, and the active force is bounded by us = allowable spring stroke (±0.08 m) *spring constant.

Due to different passenger load variations, it is assumed that the mass Ms varies between 250 to 390 kg. The system is discretized with the zero order hold method at a sampling time Ts=0.001 sec. Therefore, the discrete time norm-bounded model is obtained as

A=0.99970.0009860.0024760.00098520.055460.99690.0077830.0031060.00022281.245e50.99750.00098670.44330.024854.9340.9727E37
B=2.089e50.0046221.857e50.03694'E38

With uncertainty matrices

M=0.00040.136700.0016'E39
N=0.07290.00370.01020.0041,E40
Nb=0.0061E41

For passenger comfort, the oscillations due road bumps should be damped out within 1 s with a minimum damping ratio ζ = 0.25 [20]. To achieve the first objective, we select r1 = 0.9960 and q1 = 0 for the first circle, while for the second, r2 = 0.68 and q2 = −0.31.

The LMIs (30) are solved to get the saturated control gain matrix as

F=402.734047.9208546.454116.1682E42

As shown in Figure 5, all the closed-loop poles lie inside the desired domain.

Figure 5.

Closed-loop poles with system uncertainties using saturated control.

In general, road disturbances can be classified as shock or vibration. For shock disturbances, the proposed controller is tested for two cases: step and bump road. The frequency response is studied when the system is subjected to resonant frequency disturbance forces for the vibration case.

Case 1: step road test

For the simulation, we assume that the vehicle is subject to a 500 N unit step input due to a step road change. Two cases are considered: light and heavy passenger loads. With passive suspension, the vehicle will oscillate for an unacceptable long time, about 3.5 s, with large overshoot (Figures 6 and 7). This might damage the suspension system, whereas the active suspension damps the oscillations in about 0.3 s without overshoot (Figures 6 and 7).

Figure 6.

Step response at light load.

Figure 7.

Step response at heavy load.

Case 2: road bump test

In order to study the system response due to road bump, the case of an isolated bump in an otherwise smooth road surface is considered [9]. The corresponding ground displacement in this case is given by

x0=A21cos2πVLt,0tLV0,t>LVE43

where A and L are the height and the length of the bump and V is the vehicle forward velocity. These values are chosen as A = 0.1 m, L = 5 m, and V = 27 km/h. The bump response, namely, the suspension stroke (m) and the normalized control, is shown for two loads: light and heavy (Figures 8 and 9). In comparison with the passive suspension case, the proposed robust saturated control gives a better dynamic performance in terms of less overshoot and faster damping of the suspension stroke. Moreover, the normalized control signal is bounded between ±1 as shown in Figures 8 and 9.

Figure 8.

Road bump response at light load.

Figure 9.

Road bump response at heavy load.

4.1. Passenger comfort

The ride comfort is usually measured by the body acceleration as given in (6). This acceleration is shown in Figure 10.

Figure 10.

Passenger comfort due to road bump: light load (left) and heavy load (right).

It is evident from Figures 610 that the active suspension outperforms the passive one since it satisfies the control objectives (damps out the oscillations in <1 s and damping ratio > 0.25) as well as it does not violate the constraints of suspension stroke (8 cm) and hydraulic actuator force limit of 1.5 KN. Active suspension also provides better passenger comfort than the passive one as it has a lower body acceleration.

Case 3: frequency response

It is well known that ride comfort is frequency sensitive. From the ISO2361, the human body is much sensitive to vertical vibrations in the frequency range of 4–8 Hz. Hence, further research work is needed to evaluate the designed constrained suspension in the frequency domain. In the analysis and design of linear systems, the relationship between the time domain and frequency domain characteristics plays an especially important role. However, in system design the state space in the time domain is more convenient as done in the proposed controller. In system analysis, conditions for the robust stability and robust performance are described as specifications on the frequency response in most cases. Therefore, it is indispensable to transform frequency domain conditions into equivalent time domain conditions in order to establish an effective design theory. Parseval’s theorem gives the relationship between the squared integral of a time function and that of its Fourier transform, namely, the energy in the time domain is equal to the energy in the frequency domain. As a powerful tool bridging the time domain and frequency domain properties of systems, Kalman-Yakubovich-Popov (KYP) lemma plays a fundamental role in robust control.

5. Conclusion

This chapter considers placing the poles in a desired region of uncertain linear digital systems in the presence of input saturation. The new scheme is developed to design robust controllers, taking into consideration the effect of saturation nonlinearity. The proposed controller is based on LMI optimization and requires the uncertainty of model parameters to be cast in the norm-bounded form. The designed controller is applied to an active suspension of a quarter-vehicle model. Analysis and simulation results have confirmed the potential benefit of the proposed constrained active suspension in achieving the best possible ride comfort while keeping suspension strokes and control inputs within bounds.

In the future, it will be interesting to consider the actuator dynamic behavior in the controller design phase and also to develop robust control based on the relation between time domain and frequency domain properties.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Hisham M. Soliman, Rashid Al-Abri and Mohammed Albadi (December 20th 2017). Design of Robust Digital Pole Placer for Car Active Suspension with Input Constraint, Vibration Analysis and Control in Mechanical Structures and Wind Energy Conversion Systems, Francisco Beltran-Carbajal, IntechOpen, DOI: 10.5772/intechopen.70587. Available from:

chapter statistics

230total chapter downloads

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

Bifurcation Makes a Wave Resonant Solid-State Gyro Be Stable

By Svetlana Pavlovna Nikitenkova and Dmitry Anatolyevich Kovriguine

Related Book

First chapter

Adaptive Tuned Vibration Absorbers: Design Principles, Concepts and Physical Implementation

By Philip Bonello

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