Vibrational control of air suspension system using PID controller

Prabu, K.; Jancirani, J.; John, Dennie; Arun, B.

Journal of Vibroengineering

Browse Journal

Submit article

Published: 31 March 2013

Check for updates

Vibrational control of air suspension system using PID controller

K. Prabu¹

J. Jancirani²

Dennie John³

B. Arun⁴

^{1, 2, 3, 4}Department of Automobile Engineering, Anna University, Chennai-44, India

Cite the article Download PDF

Downloads 1935

Abstract

This paper deals with modeling and evaluation of suspension system with a pneumatic actuator controlled by Proportional Integral Derivative (PID) controller. A non-linear mathematical model of the dynamic suspension system with two degrees of freedom is developed. The controller is designed by setting proper gain values obtained by comparing three tuning methods - Ziegler Nicolas, Refined Ziegler Nicolas and Optimal control. The time response of the air suspension system is contrasted with the passive suspension system due to the road disturbance modeled as a single bump input. The proposed model limits suspension travel, minimizes passenger acceleration and keeps body displacement within bound.

1. Introduction

The suspension provides vibration isolation of the passengers from road irregularities. It also influences the traction and road holding abilities of the wheel. The vehicle passive suspension system is coil or leaf springs in combination with hydraulic or pneumatic shock absorbers. A softer spring provides good comfort and stiff spring offers better handling characteristics whereas good suspension design requires compromise between passenger ride comforts and handling. These conflicting demands have been subject of many previous studies and led to development of active suspension [1-5]. Two degree-of-freedom quarter car model is preferred to study vibrations in the vertical direction [6]. The control algorithms such as PID, Linear Quadratic (LQ) and Fuzzy logic are based on feedback control strategies [7]. The PID offers advantages of automatic tuning, gain scheduling, continuous adaptation and design simplicity. The three modes of the PID controller, proportional gain, Integral time constant and Derivative time constant, are used to get desired performance and robustness characteristics of the control system [8-10].

In this study, design and control of pneumatic actuator for nonlinear vehicle suspension system is considered. The model is approximated as a first order system. Controller design is based on Ziegler Nicolas tuning method and Optimal control. The response of the PID controller is compared for air suspension system and passive suspension system.

2. Air suspension modeling

The air suspension provides improved stability and superior comfort under different loads and road conditions. Suspension design is initiated by developing an appropriate mathematical model of the system. Certain assumptions are to be made for modeling: both vehicle body and axle assembly are rigid bodies, the tyre is a linear spring with equivalent stiffness, damping force of the tyre is neglected, passive suspension system is considered as a spring and vibration damper in parallel connection, the stiffness of the spring and damping of the vibration damper is considered and friction in the damper is neglected. The dynamically controlled pneumatic actuator is attached in parallel with the passive suspension system.

Modeling of pneumatic actuator. The pneumatic actuator model is represented in Fig. 1. Pneumatic actuator is a rubber chamber with elastic barrier which contains air under pressure.

Fig. 1Quarter car air suspension model

The vertical force exerted by actuator is the gauge pressure $(P_{g})$ multiplied with the effective area $(A_{e})$ and it can be expressed as:

1

F_{p a} = P_{g} A_{e} .

The force developed by the pneumatic actuator ( $F_{p a}$ ) can be evaluated by differentiating the above equation with respect to vertical displacement:

2

{\dot{F}}_{p a} = \frac{d {(P}_{p a} {* A}_{e})}{d y} = P_{g} \frac{{d A}_{e}}{d y} + \frac{{d p}_{p a}}{d y} A_{e} .

Consider adiabatic process for air, dynamic pressure can be given as:

3

\frac{d ({P V}^{γ})}{d y} = P γ V^{γ - 1} \frac{d (V)}{d y} + \frac{d (P)}{d y} V^{γ} = 0 .

As a result pressure increases with displacement. Due to compression, positive pressure arises, which results in decrease of volume. In the dynamic stiffness equation (4), the reservoir pressure $(P_{r})$ , gauge pressure $(P_{g})$ and atmospheric pressure $(P_{a})$ are taken into account:

4

{\dot{F}}_{p a} = γ ({P_{r} + P}_{a} + P_{g}) \frac{{A_{e}}^{2}}{V} + P_{g} \frac{d A_{e}}{d y} .

Air suspension dynamics. The pneumatic actuator suspension model is presented in Fig. 1. The model consists of sprung mass $m_{s}$ which represents a portion of the total body mass and unsprung mass $m_{u s}$ representing one of the wheel and axle. The assumption is that the mass $m_{s}$ carried by each wheel corresponds to the static weight supported by suspension system. The road input $y_{r}$ is assumed to be irregular road surface and actuator force ‘ $F_{p a}$ ’ acts on the vehicle body mass upwards. $y_{s}$ and $y_{u s}$ are the vertical displacement of the vehicle body and axle respectively. If $k_{t}$ represents the tyre radial stiffness, $k_{s}$ suspension spring stiffness. The $c_{s}$ represents the damping force available in the damper.

By using Newton’s law, the dynamic equation of motion for a suspension system is derived by the following differential equation:

5

m_{s} {\ddot{y}}_{s} + c_{s} ({\dot{y}}_{s} - {\dot{y}}_{u s}) + k_{s} (y_{s} - y_{u s}) - F_{p a} = 0,

6

m_{u s} {\ddot{y}}_{u s} + k_{t} (y_{u s} - y_{r}) - c_{s} ({\dot{y}}_{s} - {\dot{y}}_{u s}) - k_{s} (y_{s} - y_{u s}) + F_{p a} = 0 .

When the following states are introduced, the input and output vectors are:

7

x^{T} = [y_{s} {\dot{y}}_{s} y_{u s} {\dot{y}}_{u s}],

8

y^{T} = [{\ddot{y}}_{s} y_{u s} - y_{r} y_{s} - y_{u s}],

9

L^{T} = [F_{p a}, y_{r}] .

The suspension model is represented as:

10

[\begin{matrix} \begin{matrix} {\dot{y}}_{s} \\ {\ddot{y}}_{s} \\ {\dot{y}}_{u s} \end{matrix} \\ {\ddot{y}}_{u s} \end{matrix}] = [\begin{matrix} 0 & 1 & 0 & 0 \\ - \frac{k_{s}}{m_{s}} & - \frac{c_{s}}{m_{s}} & \frac{k_{s}}{m_{s}} & \frac{c_{s}}{m_{s}} \\ 0 & 0 & 0 & 1 \\ \frac{k_{s}}{m_{u s}} & \frac{c_{s}}{m_{u s}} & \frac{- (k_{s} + k_{t})}{m_{u s}} & - \frac{c_{s}}{m_{u s}} \end{matrix}] [\begin{matrix} \begin{matrix} y_{s} \\ {\dot{y}}_{s} \\ y_{u s} \end{matrix} \\ {\dot{y}}_{u s} \end{matrix}] + [\begin{matrix} o & 0 \\ \frac{1}{m_{s}} & 0 \\ 0 & 0 \\ - \frac{1}{m_{u s}} & \frac{k_{t}}{m_{u s}} \end{matrix}] [\begin{matrix} F_{p a} \\ y_{r} \end{matrix}],

11

[\begin{matrix} {\ddot{y}}_{s} \\ y_{u s} - y_{r} \\ y_{s} - y_{u s} \end{matrix}] = [\begin{matrix} - \frac{k_{s}}{m_{s}} & - \frac{c_{s}}{m_{s}} & \frac{k_{s}}{m_{s}} & \frac{c_{s}}{m_{s}} \\ 0 & 0 & 1 & 0 \\ 1 & 0 & - 1 & 0 \end{matrix}] [\begin{matrix} \begin{matrix} y_{s} \\ {\dot{y}}_{s} \\ y_{u s} \end{matrix} \\ {\dot{y}}_{u s} \end{matrix}] + [\begin{matrix} \frac{1}{m_{s}} & 0 \\ 0 & - 1 \\ 0 & 0 \end{matrix}] [\begin{matrix} F_{p a} \\ y_{r} \end{matrix}] .

The suspension travel is the relative displacement between the sprung mass and the unsprung mass. Tyre displacement is the difference between the unsprung mass displacement and the road surface at any instant.

3. PID control system

The performance of PID controller for air suspension is enhanced by proper selection of the parameters proportional gain (K_P), integral gain (K_i), derivative gain (K_d) values. This paper attempts to design PID controller using Ziegler Nicolas, refined Ziegler Nicolas, and optimal control tuning methods. A PID control system shown in Fig. 2.

Fig. 2PID controlled suspension system

Error signal $e (s)$ is estimated from the differences in reference input signal $r_{i} (s)$ to the system output $y (s)$ . The control signal $u_{i} (s)$ from the controller to the plant model is equal to the proportional gain ( $K_{p}$ ) times the magnitude of the error plus the integral gain ( $K_{i}$ ) times the integral of the error plus the derivative gain ( $K_{d}$ ) times the derivative of the error. The PID controller is mathematically represented as:

12

u_{i} (s) = K_{p} [e (t) + K_{i} \int_{0}^{t} e (τ) d τ + K_{d} \frac{d e (s)}{d t}] .

Best control performance is achieved when choosing appropriate values of proportional gain, integral time, derivative gain. The selection of the PID parameters is called as tuning and the best tuned method is identified by comparing the performances of three different tuned systems, Ziegler-Nichols, refined Ziegler-Nichols and optimal tuning methods. The Error signal is denoted by the following equation:

13

e (s) = y (s) - r_{i} (s) .

The suspension plant model is approximated to a first order plus dead time (FOPDT) model using [9] standard transfer function approach. From the FOPDT model transfer function the transportation lag $L$ , time constant $T$ and the static gain $k$ are obtained. The above parameters are used to determine the controller parameters. Evaluating the equation at $s =$ 0, yields the model parameters obtained as: $T =$ 7.0134, L = 2.9867, $k =$ 0.547 and $a = \frac{k L}{T}$ .

The Ziegler-Nichols PID tuning is a commonly used feedback loop in control systems. This approach is utilized for comparison purpose only.

Refined Ziegler-Nichols tuning. The refined Ziegler-Nichols settings are introduced by adding a new parameter ‘ $β$ ’ in the proportional action of the controller [11]. The new parameter ‘ $β$ ’ can reduce the overshoot to acceptable levels, and thus gives good set point response. When ‘ $β$ ’ value is less than unity the control equation becomes:

14

u_{i} (s) = K_{p} (β e (t) + \frac{1}{T_{i}} \int_{0}^{t} e (τ) d τ) - K_{p} ((1 - β) y + T_{d} \frac{d y}{d t}) .

The PID gain values are defined and on further simplification the equation is refined as:

15

K_{p} = (\frac{6 + 0.5 c}{9 + 8.4 c}), T_{i} = (\frac{4 c + 15}{75}) .

Optimal control method tuning. The general form of the optimal PID controller design should have good rejection performance on the given disturbance:

16

J_{n} (θ) = \int_{0}^{\infty} {[t^{n} e (θ, t)]}^{2} d t .

Tuning methods are proposed [12] for the disturbance signal $d (t)$ and the error signal is $e (θ, t)$ . The gain values $K_{P}$ , $T_{i}$ , $T_{d}$ of the PID controller are evaluated from the following relations:

17

K_{p} = \frac{a_{1}}{T} {(\frac{L}{T})}^{b_{1}}, T_{i} = \frac{T}{a_{2}} {(\frac{L}{T})}^{b_{2}}, T_{d} = a_{1} T {(\frac{L}{T})}^{b_{2}} .

4. Simulation

The vehicle parameters [13] (Table 1) of quarter car model are used in this analysis. The simulation experiment was carried for with the single bump disturbance input.

Table 1Vehicle parameters

Vehicle body mass in kilograms	285.3	Tyre stiffness in N/m	190000
Axle and tyre mass in kilograms	59.5	Suspension damping in Ns/m	1000
Suspension spring stiffness in N/m	16812

5. Single bump road input

A road surface is an unsystematic with sudden speed bumps. Here air suspension system response is analyzed by a single speed bump (Fig. 3). To limit the effort of the random input, simpler road model is used [14]. It consists of a crest of height 5 cm each. A road input, y_r of height 0.05 m, is described by setting ‘ $a$ ’ as 0.025 in the below equation and it is used to study the system response.

Fig. 3Road bump

The road bump is mathematically described [15] as:

18

y_{r} = \{\begin{matrix} a (1 - c o s 8 π t), 0.5 \leq t \leq 0.75, \\ 0, o t h e r w i s e . \end{matrix}

6. Results and discussions

From the results, it is found that addition of pneumatic actuator plays an extensive work in the reduction of vehicle body displacement, suspension travel and body acceleration when compared to that of passive suspension system. The road bump is shown in Fig. 4.

Fig. 4. Single bump road input

Fig. 5Vehicle body displacement for 50 mm single bump road input

For given road input the peak vehicle body displacement (Fig. 5) in optimal (0.0197 m) is lower than passive (0.0370 m) suspension system. The Suspension travel (Fig. 6) of optimal PID system is better than passive when considering the peak overshoot, system response and settling time. The vehicle tyre displacement (Fig. 7) for optimal PID system shows marginal improvement when compared to passive system.

Based on the observation from vehicle body acceleration (Fig. 8), the Ziegler-Nichols method offers minimum overshoot but the settling time is too long when compared to other control methods. From the values it is observed that when both the body acceleration and the settling time are considered, the optimum result is obtained by optimal PID system.

From the Fig. 9, the body acceleration is reduced by 41.4 % and an improvement in settling time of 11.4 % for optimal control method when compared with a passive system. The vehicle body displacement, suspension travel and tyre displacement (Fig. 10) are compared, and it is found that the settling time and the response for the optimum PID controller gives favorable results.

Fig. 6Vehicle suspension travel for 50 mm single bump road input

Fig. 7Vehicle tyre displacement for 50 mm single bump road input

Fig. 8Vehicle body acceleration for 50 mm single bump road input

Fig. 9Maximum body acceleration for different PID controls

Fig. 10Maximum overshoot for various parameters of the system

7. Conclusions

The objective of the study is to optimize the performance of suspension system with a pneumatic actuator controlled by PID controller. Considering the dynamics of the actuator, the different control methods are evaluated by a single speed bump rather than complicated road model as the scope of study is to analyze the controller response to vibration. It is concluded that:

• The body acceleration is reduced and an improvement in settling time for optimal control method.

• The vehicle body displacement, suspension travel and tyre displacement are compared and it is found that the settling time and the response for the optimum PID controller gives better results than passive system.

• Based on the overall performance, it is concluded that the optimal PID is the most suitable control method for the air suspension system.

References

J. Y. Wong Theory of Ground Vehicles. John Wiley & Sons, 2001.

Search CrossRef
Uisoy A. G., D. Hrovat Stability robustness of LQG active suspensions. Proceedings of the American Control Conference, San Diego, California, 1990, p. 1347-1356.

Search CrossRef
Alleyne A., Hedrick J. K. Nonlinear adaptive control of active suspensions. IEEE Trans. on Control Systems Technology, Vol. 3, No. 1, 1995.

Search CrossRef
Hrovat D. Applications of optimal control to advanced automotive suspension design. Journal of Dynamic Systems, Measurement and Control, Vol. 115, No. 2 (B), 1993, p. 328-342.

Search CrossRef
Rakheja S., Afework Y., Sankar S. An analytical and experimental investigation of the driver seat suspension system. Vehicle System Dynamics, Vol. 23, 1994, p. 501-524.

Search CrossRef
K. Prabu, J. Jancirani Development of a new suspension system for passengers to improve occupants comfort. Journal of Mechanical Engineering, Vol. 1, No. 4, 2011, p. 7-13.

Search CrossRef
V. Gavriloski, D. Danev, K. Angushev Mechatronic approach in vehicle suspension system design. 12^th IFToMM World Congress, Besançon (France), 2007.

Search CrossRef
Danlel E. Rlvera, Manfred Morarl, Slgurd Skogestad Internal model control. 4. PID controller design. Ind. Eng. Chem. Process Des. Dev., Vol. 25, 1986, p. 252-265.

Search CrossRef
Li Mei, Li Zhongxing, Shen Xufeng, Guo Jiwei Study on PID control for semi-active air suspension for riding comfort. 2nd WRI Global Congress on Intelligent Systems, 2010, p. 47-50.

Search CrossRef
C. C. Hang, K. J. Wstrom, W. K. Ho Refinement of the Ziegler-Nichols tuning formula. IEEE Proceedings D, Vol. 138, No. 2, 1991, p. 111-118.

Search CrossRef
Dingyu Xue, Yang Quan Chen, Derek P. Atherton Linear feedback control: analysis and design. The Society for Industrial and Applied Mathematics, 2007, p. 187-209.

Search CrossRef
O. A. Dahunsi, J. O. Pedro, O. T. Nyandoro System identification and neural network based PID control of servo-hydraulic vehicle suspension system. South African Institute of Electrical Engineers, Vol. 101(3), 2010, p. 93-105.

Search CrossRef
R. Rajamani, J. K. Hedrick Performance of active automotive suspensions with hydraulic actuators: theory and experiment. Proceedings of American Control Conference, Maryland, 1994.

Search CrossRef
K. Prabu, J. Jancirani Investigation on the performance of the quarter car model suspension system using electromagnetic actuator. IJMAE, Vol. 20(1), 2012, p. 1-8.

Search CrossRef
Jung-Shan Lin, Ioannis Kanellakopoulos Nonlinear design of active suspensions. IEEE Control Systems Magazine, Vol. 17(3), 1997, p. 49-59.

Search CrossRef

About this article

Received

11 September 2012

Accepted

28 February 2013

Published

31 March 2013

Keywords

dynamic suspension design

PID controller

air suspension system

road model

air spring model

Acknowledgements

Research work has been supported by Anna University Chennai and the authors are grateful to the University.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.