# Optimal design and analysis of quarter vehicle suspension system by using Matlab.

1. Introduction

Nowadays the focus of vehicle suspension design has switched from pure numerical analysis to the application of different algorithms which are designed based on optimization methods. The development of optimizing methods is connected with the development of an automatic control of systems, which is the case in our study.

Algorithms designed based on optimizing methods enables us to fasten the search process of finding optimal solution expressed by an objective function under the fulfillment of constraint conditions, which evidently has a great influence on the development of the whole group of modern control methods that for the calculation of the control effect utilize a mathematical model of the process determined based on the measured process data (MPC methods and algorithms). At present, a great success has been achieved by agent and colonic control methods as well as methods and genetic algorithms used in the system of off-line techniques applied to the searching for a global extreme of a goal function. The applications of the mentioned methods and their algorithms to nonlinear processes modeling and control, systems, significantly contribute to the improvement of approximation properties of modeling processes as well as to the improvement of a control quality of selected types of a control . Genetic algorithms (GAs) have been used in various applications such as function optimization, system identification and control systems. GAs is general-purpose stochastic optimization methods for solving search problems to seek a global optimum. However, GAs is characterized by a large number of function evaluations .

On the other hand, traditional methods, such as sequential quadratic programming (SQP), are well known to exploit all local information in an efficient way, provided that certain conditions are met and the function to be minimized is 'well-conditioned' in the neighbourhood of a unique optimum. These methods require adequate local information to be known (such as the gradient and Hessian matrix). If the basic requirements are not satisfied, the reliability of the SQP method is greatly jeopardized.

The aim of this study is, for quarter car model to find optimal design parameters, by minimising a nonlinear objective function subjected to a set of constraints with two different optimisation algorithms; SQP and GA in order to improve ride quality. By means of the ride quality analysis in the frequency domain, the vertical vehicle body acceleration (VBA), suspension working space (SWS) and dynamic tire load (DTL) can be obtained .

In this design optimization, the main objective is to minimize the VBA acceleration. In the meantime, the SWS and DTL are constrained. If the SWS is too small, the sprung mass will strike the un-sprung mass and this may lead to damage of the vehicle. If the DTL is greater than the static tire load, the vehicle's tires will bounce off the road  and this will result in unstable modes of vehicle motion. Therefore, we aimed to present the development of a method for obtaining poles, eigenvectors, natural frequencies and mode shapes, in order to perform a vibration analysis for a quarter car model by using three different approaches , ,  and comparison of the natural frequencies which have been obtained for the optimised model with the natural frequencies defined by the laws of the vibration theory. The parametric analysis gave information on how the natural frequency is affected when one of the parameters is perturbed. Parametric studies of dynamical systems are an important task in all major engineering design. This is especially true in the design of control systems. It is important to know when one parameter is changed, how that affects the system; if the system becomes unstable, do the natural frequencies change, etc. In addition, we presented a task of state space modelling, Laplace Transformation, Transfer Function, Stability, Controllability, Observability as one of most important attributes to analyse quarter car vehicle model.

Another aim was to develop a Fuzzy Logic Controller (FLC) for a quarter car model with optimised parameters, with a rule base , , which can be tuned to improve each of three objectives; ride comfort, suspension working travel and handling.

In addition, in this study we have applied in integrative way software programs MATLAB and SIMULINK for optimisation, analysis, design and control of a quarter car model.

2. Dynamical model of quarter car vehicle model

The model of the quarter-car active suspension system used in this paper with two degree of freedom is shown in Fig. 1. The model represents a single wheel of a car in which the wheel is connected to the quarter portion of the car body through a hydro pneumatic suspension. The dynamics of quarter car model in Fig. 1 are described by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

The model parameters and their respective units are summarized in Tab. 1.

Equations (1) can be represented in matrix form as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

where M, C, and K represent the mass, damping, and stiffness matrices described by:

M = diagonal [[m.sub.k], [m.sub.p]] (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

Standard second order matrix form for the passive quarter car suspension system is shown in the following:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

A x z - B x z = 0 (7)

The response of this system takes the form:

z(t) = v x [e.sup.[lambda]-t] (8)

where v is a constant vector.

The model represents a single wheel of a car in which the wheel is connected to the quarter portion of the car body through a hydropneumatic suspension. Substituting (8) in (7) leads to eigenvalue problem;

(A - [lambda] x B) x v = 0 (9)

Vector v is called the eigenvector and X represents eigenvalues.

3. Natural frequencies and mode shapes simulation in Matlab

Parametric studies of dynamical systems are an important task in all major engineering design. This is especially true in the design of control systems, when the focus is for low sensitivity to parameter variations . It is important to know when one parameter is changed, how that affects the system; if the system becomes unstable, do the natural frequencies change, etc.

Since we have obtained matrices M, C, and K, we can calculate natural frequencies and mode shapes in MATLAB; Natural Frequencies are:

w = 6.5612 0 and 0 79.5945

Mode shapes are:

u= -0.0707 -0.0012 -0.0031 0.1741

Through simulation are obtained the values from V1 to V4 of a complex conjugate pairs of vectors using MATLAB with three different programming approaches. The simulation results are shown in Table 2.

First approach:

a=[m1 cs k1]; % poles of the system r = -24.3582 +73.2172i -24.3582 -73.2172i -3.8842 + 5.5424i -3.8842 - 5.5424i

Second approach,

r = -7.1146 -15.0873i -7.1146 +15.0873i 0.0009 + 0.0011i 0.0009 - 0.0011i

p = -24.3582 +73.2172i -24.3582 -73.2172i -3.8842 + 5.5424i -3.8842 - 5.5424i

k = 0.0050 -0.3127

Third approach: Frequency response and poles of the system

v = -0.0011 + 0.0008i -0.0011 - 0.0008i -0.0991 - 0.0395i -0.0991 + 0.0395i -0.0042 - 0.0120i -0.0042 + 0.0120i 0.0002 - 0.0049i 0.0002 + 0.0049i -0.0317 - 0.0984i -0.0317 + 0.0984i 0.6039 - 0.3961i 0.6039 + 0.3961i 0.9834 - 0.0166i 0.9834 + 0.0166i 0.0264 + 0.0200i 0.0264 - 0.0200i

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As we can see for all cases similar results have been obtained.

4. Stability of the model by Lyapunov theorem in Matlab

The Lapunov criterion will be used to check whether or not the given system is stable. Starting point for stability analysis is state-space form of system , . Matrix A must be definite, x is the state vector containing n state variables. The quadratic form of Lyapunov function is:

V(x) = [x.sup.T] (t) x Q x x(t) (10)

where Q must be a positive definite symmetric matrix nxn. General form of Lyapunov matrix equation is:

A x Q + Q x [A.sup.T] = -C (11)

where C is a positive definite symmetric matrix. Usually, matrix Q is taken Q = I.

From MATLAB, for closed loop is:

eig(A-B*K)

ans = -24.3567 +73.2176i -24.3567 -73.2176i -3.8842 + 5.5424i -3.8842 - 5.5424i

5. Full Controllability

If any modes of the open-loop system are unstable, they can be stabilized with a feedback using linear quadratic control. A necessary condition for solving the controller problem is that the system is controllable. This means that all the modes can be excited or controlled by the input u. In terms of the state space matrices, the condition of controllability can be shown by constructing the controllability matrix Qc defined by:

[Q.sub.c] = [B, AB, ..., [A.sup.m-1]B] (12)

Controllability is ensured when the rank of [Q.sub.C] equals m, where m is the order of the system.

6. State Controllability

A system is said to be "(state) controllable" if for any to and any initial state x([t.sub.0])= [x.sub.0] and any final state [x.sub.f], there exists a finite time [t.sub.i]> [t.sub.0] and control u(t), such that x([t.sub.1]) = [x.sub.f].

A LTI system with the state equation x = Ax + Bu is controllable if the controllability matrix

Qc = [B, AB, [A.sup.2]B, ..., [A.sup.m-1]b] has rank m.

Controllability

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Since rank of Qc is 4 equals with size of matrix A, system is controllable.

7. State Observability

A system is said to be "(state) observable" if for any [t.sub.0] and any initial state x([t.sub.0]) = [x.sub.0] there exists a finite time [t.sub.1] > [t.sub.0] such that knowledge of u(t) and y(t) for [t.sub.0] [less than or equal to] t [less than or equal to] [t.sub.1] suffices to determine [x.sub.0] A LTI system with the state space model

[??] = Ax + Bu y = Cx (13)

Only if observability matrix [Q.sub.0] = [[C, CA, C[A.sup.2], ..., C[A.sup.n-1]].sup.T] has rank n. Functions in MATLAB

Observability

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Since rank of Qo is 4 equals with size of matrix A, system is observable.

8. Stochastic Road Modeling and PSD Response

Road irregularity or unevenness represents the main disturbing source for either the rider or vehicle structure itself. The road profile elevation is usually expressed in terms of the power spectral density (PSD). The PSD of the road profile elevation is expressed as:

[G.sub.q] (n) = [G.sub.q] ([n.sub.0])[(n / [n.sub.0]).sup.-w] (14)

For the purposes of design optimization, according to James' principle, the root mean square (RMS) of the sprung mass acceleration [[??].sub.2] can be expressed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

The RMS of the suspension working space [f.sub.d] is:

[[sigma].sub.fd] = [{[pi]RV[([m.sub.1] + [m.sub.2])[([m.sub.1][k.sub.1]).sup.1/2]/2[m.sub.1][c.sub.s]]}.sup.1/2] (16)

The RMS of the relative dynamic tire load can be calculated as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

9. Design Optimisation

In this section, the sprung mass vertical acceleration is minimized, while the design constraints on the suspension working space and dynamic tire load should be satisfied. To implement the design optimization, the two optimization algorithms, i.e, SQP and GA, will be applied, respectively.

Optimization based on SQP algorithm and GA

The SQP algorithm is a non-linear programming technique that is used for the purpose of minimizing a smooth nonlinear function subjected to a set of constraints with upper and lower bounds. The objective function and the constraint functions are assumed to be at least twice continuously differentiable. This algorithm is a gradient-based search method ,. This algorithm is well-suited for constrained design optimizations. The reliability for finding the optimum decreases with the increase of number of design variables when using SQP method. In contrast, whether the number of design variables increase the GA can still reliably find the optimum. This can be explained by the fact that GA works on a population of variables in parallel, not on a unique point. GAs are a global search methods that are based on the Darwin's principle of natural selection and genetic modification and have higher reliability for the global optimum with minimum number of computational operations.

The RMS of the acceleration of a sprung mass [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is frequently used to evaluate the riding quality of a vehicle. A rider's comfort improves as the acceleration decreases. Ride comfort is chosen to be the design criterion. The suspension working space and dynamic tire load [[sigma].sub.fd] are selected as the design constraints. The design variables are [m.sub.1], [m.sub.2], [k.sub.1], [k.sub.2] and [c.sub.s], respectively.

Thus, the design optimization problem can be obtained by minimisation of sprung mass acceleration ,  given in (18) subject to constraints (19).

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

Subject to:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19)

In this sub-section, the optimization results are derived for a vehicle travelling at the speed of 40 m/s on the road with an irregularity coefficient of power spectrum taking the value of 6.5x[10.sup.-6] [m.sup.3].

10. Fuzzy Logic Controller

Fuzzy Logic Control has accelerated in recent years in many areas, including feedback control. By using empirical rules according to the designers knowledge and experience, which are represented linguistically with the conditional statements and resulting assertion.

A fuzzy rule base has a very significant effect on the control strategy in FLC, in other words it defines the strategy of the controller. To the active suspension system there are at least three main objectives, namely ride comfort, suspension travel and handling. The rule base can be tuned to improve each of the above objectives.

The fuzzy logic controller used in the active suspension has three inputs: body acceleration [[??].sub.b], body velocity [[??].sub.b], body deflection velocity [[??].sub.b] - [[??].sub.t] and one output: desired actuator force [f.sub.a], shown in Fig. 2.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20)

The output of the fuzzy controller [[mu].sub.D] is corresponding membership function.

The control system itself consists of three steps: fuzzification, fuzzy inference machine and deffuzification. During the fuzzification process the real numbers (crisp) inputs will be converted into fuzzy values, where after fuzzy interference machine processes the input data and computes in cope with the rule base and database. The obtained outputs (fuzzy values) are converted into real numbers by the defuzzification step. Membership functions are chosen for the inputs and the output variables: NV-negative very big, NB-negative big, NM-negative medium, NS-negative small, N-negative, ZE- zero, P-positive, PS-positive small, PM-positive medium, PB-positive big, PV-positive very big.

The fuzzy rule based system modeled by designer's knowledge and experience is shown in Fig. 3. , , , [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], wliich are characterised by their membership functions: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Simulation results for acceleration of unsprunged mass for passive quarter car model are shown in Fig. 4. Accelerations of the active quarter car model are shown in Fig. 5.

11. Conclusions

A comparative study of two optimization algorithms (genetic algorithms, GAs and sequential quadratic programming, SQP), has been conducted through minimizing the vertical sprung mass acceleration subjected to a suspension working space and dynamic tire load.

By optimizing the design parameters compared with the original design, the sprung mass (body) acceleration decreases. The suspension working space and the dynamic tire load satisfy the specified design constraints. Based on the simulation results the optimum found by GAs at 19 generations, while by using the SQP the optimum is found after 51 iterations. As shown from the numerical simulation results the max amplitude of body displacement using optimized design variables is reduced for 9%, while maximal amplitude of body acceleration is reduced around 22%.

Numerical experiments reveal the fact that to improve vehicle ride quality and satisfy the specified suspension working space and relative dynamic tire load, different vehicle speed and road irregularity have different requirements on the design variables, in particular, the un-sprung mass.

From the parametric analysis it can be concluded that mass is a very sensitive parameter, and must be handled with special attention, while the results of the natural frequencies obtained in both ways shows that the values are almost the same.

Active suspension systems using FL Controllers can reduce vertical accelerations, as shown in Fig. 5. This means, that the main properties; ride comfort and road holding were achieved for the quarter car model. For the design of an FLC, an accurate vehicle model is not needed, but it's a very difficult task to express the knowledge and experience in terms of fuzzy logic.

From the results of this study it became evident to us that the combination of MATLAB and SIMULINK represent a powerful tool on research of parametric analysis, optimisation and optimal control of vehicle suspension system. Since the quarter car model is one the bests and represents the simplest model for the analysis of ride comfort, the same approach will be applied in our future research to analyse half car and full car vehicle suspension system

DOI: 10.2507/27th.daaam.proceedings.012

12. References

 Jadlovska, A.; Katalinic, B.; Hrubina, K., & Wessely, E. (2013). On Stability of Nonlinear Systems and Application to Apm Modeling, Chapter 11 in DAAAM International Scientific Book 2013, pp. 257-276, B. Katalinic & Z. Tekic (Eds.), Published by DAAAM International, ISBN 978-3-901509-94-0, ISSN 1726-9687, Vienna, Austria, DOI: 10.2507/daaam.scibook.2013.11

 Jadlovska, A.; Katalinic, B.;Hrubina; K.; Macurova, A.& Wessely, E. (2012). Solution to the Problem Control of a Distributed Parameter Process, Chapter 15 in DAAAM International Scientific Book 2012, pp. 169-186, B. Katalinic (Ed.), Published by DAAAM International, ISBN 978-3-901509-86-5, ISSN 1726-9687, Vienna, Austria DOI:10.2507/daaam.scibook.2012.15

 Jadlovska, A.; Katalinic, B.;Hrubina, K.; Macurova, A. & Wessely, E. (2011). Optimal Control of Nonlinear Systems with Constraints, Chapter 22 in DAAAM International Scientific Book 2011, pp. 265-282, B. Katalinic (Ed.), Published by DAAAM International, ISBN 978-3-901509-84-1, ISSN 1726-9687, Vienna, Austria DOI:10.2507/daaam.scibook.2011.22

 Likaj R.; Shala A.; Bruqi M.; Qelaj M. (2009). Optimal Design of Quarter Car Vehicle Suspension System, Trends in the Development of Machinery and Associated Technology - ISSN 1840-4944; pp. 417-420; Livorno, Italy

 Likaj R.(2005). Fuzzy Logic Control on Nonlinear Vehicle Suspension System, Phd thesis, Prishtina, Kosovo

 Likaj, R.; Geca, A.; Pajaziti, A.; Osmanaj S. (2005). Parametric and sensitivity analysis of four DOFs vehicle model, 4th DAAAM International Conference on Advanced Technologies for developing Countries, September 21-24, 2005, Slavonski Brod, Croatia

 Likaj R.; Pajaziti A., (2002). Nonlinear and fuzzy Logic Control of Vehicle Suspensions, "Intelligent Manufacturing & Automation: Learning from Nature" 13th International DAAAM Symposium; Annals Proceeding, ISBN 3901509-29

 Cai, B.; Konik, D. (1993). Intelligent Vehicle Active Suspension Control using Fuzzy Logic, IFFAC World Congres, Vol.2, pp. 231-236..

 Zhongzhe Ch., Yuping H., Natarer G.: Design Optimization of suspension with a Quarter- Vehicle Model, CSME/ de la SCGM, Vol.32, No .2, pp 297-312, 2008.

 Yuping H: Design Optimization of a Quarter Vehicle suspension system using a Genetic Algorithm and SQP, Department of Systems Design Engineering, University of Waterloo, 1999

Caption: Fig. 1. Quarter car suspension model of 2 DOF

Caption: Fig. 2. Control scheme of quarter car model

Caption: Fig. 4. Acceleration of unsprunged mass for passive and active quarter car model

Caption: Fig. 5. Acceleration of sprunged mass for active quarter car model
```Table 1. Quarter car vehicle model parameters

Symbol        Description       Value

[m.sub.1]      Body mass        200 kg
[m.sub.2]     Wheel mass        33 kg
[c.sub.s]    Damping ratio    1600 Ns/m
[k.sub.1]   Spring constant    9000 N/m
[k.sub.2]   Spring constant   200000 N/m

Table 2. Poles of the system

Poles         Value

s11     -24.3582 +73.2172i
s22     -24.3582 -73.2172i
s33     -3.8842 + 5.5424i
s44     -3.8842 - 5.5424i

Table 3. Optimal design variables based on the SQP and GA for
minimizing the sprung mass vertical acceleration, vehicle speed
40 m/s

Original values            SQP method

Initial values                            [10 10 10 10 10]
[m.sub.1] [kg]          200                     200.0
[m.sub.2] [kg]          33                      32.0
[k.sub.2] [N/m]       200000                  200000.0
[k.sub.1] [N/m]        9000                    9100.0
cs [Ns/m]              1600                    1580.0
[MATHEMATICAL                              1.2673254442031
EXPRESSION NOT
REPRODUCIBLE IN
ASCII]
Optimum found: 52 iterations

GA

Initial values           [10 10 10 10 10]
[m.sub.1] [kg]           199.918612498565
[m.sub.2] [kg]          32.04226032330942
[k.sub.2] [N/m]         200305.81925474654
[k.sub.1] [N/m]          9099.82776078753
cs [Ns/m]               1591.7473971815848
[MATHEMATICAL             1.26976709377
EXPRESSION NOT
REPRODUCIBLE IN
ASCII]
Optimum found: 51 generations

Fig. 3. Fuzzy rules

[[??].sub.b]   [[??].sub.b]   [[??].sub.b]-[[??].sub.t]   [f.sub.a]

ZE                  P                     P                  NS
ZE                  P                    ZE                  NM
ZE                  P                     N                  NB
ZE                  ZE                    P                  PS
ZE                  ZE                   ZE                  ZE
ZE                  ZE                    N                  NS
ZE                  N                     P                  PB
ZE                  N                    ZE                  PM
ZE                  N                     N                  PS
P or N              P                     P                  NM
P or N              P                    ZE                  NB
P or N              P                     N                  NV
P or N              ZE                    P                  PM
P or N              ZE                   ZE                  ZE
P or N              ZE                    N                  NM
P or N              N                     P                  PV
P or N              N                    ZE                  PB
```