Observer design based on nonlinear suspension model with unscented Kalman filter

Ying Fan1 , Hongbin Ren2 , Sizhong Chen3 , Yuzhuang Zhao4

1, 2, 3, 4School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China

4Corresponding author

Journal of Vibroengineering, Vol. 17, Issue 7, 2015, p. 3844-3855.
Received 17 June 2015; received in revised form 16 August 2015; accepted 4 September 2015; published 15 November 2015

Copyright © 2015 JVE International Ltd. 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.
Creative Commons License
Abstract.

This paper presents a new approach to estimating suspension state information and parameter in real-time. An observer with unscented Kalman filter is designed based on a nonlinear quarter car model. The proposed observer could estimate the sprung mass, vertical velocity of sprung and unsprung mass for the nonlinear suspension systems with vehicle load variation. The designed observer has low sensitivity and robust to unknown road surfaces. The efficiency of the estimator is validated through the simulations with two different types of road excitation and payload variations. The simulation results clearly indicate that compared with the extended Kalman filter estimator, the unscented Kalman filter is more accurate and robust. The estimated state information and parameters could be used in the design of suspension control systems.

Keywords: nonlinear suspension, unscented Kalman filter (UKF), observer.

1. Introduction

In recent years, there has been growing interest in the development of suspension system control strategy [1-3], such as active and semi-active suspension systems, because of its high performance of isolating the passengers and goods from the road disturbance. The control effect greatly relies on the accuracy of the suspension state information and parameters. Some of the required state information is easy to be measured by the sensors which already exist in modern vehicle, such as the sensor of vehicle body vertical acceleration; but others are difficult to be measured due to the expense, complexity and technological limitations, such as suspension vibration velocity or displacement. In practical application, a common technique to obtain the sprung mass vertical velocity is to integrate signals from accelerometers. However, it is well known that integrator-based estimations present larger errors due to drift [4]. Eliminating the number of sensors is a potential approach to cutting down the cost, improving the reliability and performances of the suspension controller [5]. State observer is a system which provides the internal state of a given system in real-time; it inputs the online measurements to the observer and outputs the estimated values. So it has gained increasing attentions in the field of automobile industry [6-8].

The performance of the controllable suspension can be further improved if accurate state information could be obtained in real-time. Suspension state estimation based on Romberg observer with low sensitivity to unknown road surfaces was designed in [9]; the results indicated that low gain Romberg observers presented the best overall performance when considering a wide series of perturbation scenarios. Best and Gordon [10] proposed an online system identification to identify mass and stiffness by using the recursive least-squares (RLS). Hedrick [11] designed an adaptive observer to estimate state and parameters in the active suspension system. System state and parameter identification is mainly based on two different principles. One is estimation algorism; the other is based on physical sensor, such as 3D camera. But the sensor-based technical is not widely applied, because of it is heavily depending on the high speed processor and cost. Suspension vibration is mainly from road disturbance, some researchers have tried to obtain the state information indirectly by designed observer [12-14] or directly from vision sensors. A system with preview controller employing a vehicle mounted preview sensor was proposed in [15]. But the road information is difficult and expensive to be measured directly in practice.

Kalman filtering (KF) is an optimal estimator which estimates unknown variables from measured noisy signals. The KF has numerous applications in the technology. Its common application is identifying the vehicle state and parameters, particularly in aircraft and spacecraft. Whereas, the linear KF cannot be directly used in nonlinear systems, so some improved techniques are proposed, such as extended Kalman filter (EKF) [16], UKF [17, 18] and feedback linearization observer [19] et al.

The main motivation of this paper is to investigate the state and parameter estimation of the nonlinear suspension systems. An observer algorithm is designed to estimate the sprung and unsprung mass vertical velocity by measuring the suspension stroke and the stroke rate based on the UKF theory; and another observer to identify the vehicle payload by measuring the sprung mass vertical acceleration based on the linear KF theory is proposed. Both observers exchange the estimated results in real-time. These estimated states and parameter are necessary for controllable suspension systems to determine the control signal.

The paper is organized as follows: a nonlinear quarter car model and the dynamic equations are explained in Section 2; while in Section 3, the state observer is designed based on UKF and payload identification is also introduced in detail. The efficiency of the proposed observer is illustrated by simulation results in Section 4. At the end of this paper, the conclusions and future work are given.

2. Modeling of nonlinear quarter car

The linear model can replace the nonlinear one around the operation conditions, out of this level, the linear model is not valid, and a linear representation of the system dynamics is not sufficient [20]. The quarter car model of the nonlinear suspension system is shown in Fig. 1.

Fig. 1. Nonlinear quarter car model

 Nonlinear quarter car model

m s is the sprung mass of the car body, including passengers and internal components; it varies according to the number of passengers and the payload condition of the car. It is supported by the suspension system including air spring and damping adjustable absorber, and both of them have nonlinear characteristics. The shock absorber is a continuous adjustable damper. The damping force can be adjusted by controlling the current of solenoid valve according to a given logic. mu is the unsprung mass which is supported by the tire modeled as linear spring with stiffness coefficient kt. The displacement of the sprung mass and unsprung mass are denoted as zs and zu respectively; and the road profile is denoted as q. So the dynamic equations of the nonlinear suspension can be described as:

(1)
m s z ¨ s = - F s - F d ,
(2)
m u z ¨ u = F s + F d - k t z u - q ,
(3)
F s = F s z s - z u ,
(4)
F d = F d z s - z u , z ˙ s - z ˙ u , I ,

where, Fs is the air spring force; which is the function of damper stroke, as shown in Fig. 2; Fd is the damping force of the continuously adjustable damper; which is the function of damper stroke, stroke rate and the control current, as shown in Fig. 3. In the simulation, a 3-D look-up table model is created to describe the nonlinear characteristics of the adjustable damper.

Fig. 2. The characteristic of the air spring force

The characteristic of the air spring force

Fig. 3. Damping force-velocity characteristics for different control currents

Damping force-velocity characteristics  for different control currents

3. UKF observer design

Some suspension control algorithms assume that all states are measureable. But in practice, some states are either difficult to be measured or cannot be measured at all due to the high cost or limitations of technology. So it is necessary to design the observer to identify the states or parameters of the suspension systems. Many estimation algorithms for the nonlinear systems have been proposed in recent years, such as EKF [21], UKF [22], sliding mode observer [23] et al.

EKF is a nonlinear version of KF to optimally estimate the state information and parameters for the nonlinear systems. It is based on the Taylor series expansion theory and first order approximation to linearize the nonlinear system. This method induces a large error in the actual values of mean and covariance, which may lead to sub-optimal performance and sometimes even cause divergence of the filter [24]. Furthermore, it is difficult to calculate the Jacobi matrix during the system linearization especially for the strong nonlinear systems, such as the suspension with adjustable damper.

UKF is a nonlinear Kalman Filter which avoids calculating the Jacobi matrix and shows superior accuracy compared with the EKF which works with a linearized nonlinear model. It is based on the Unscented Transform (UT) theory and statistical linearization technique. This technique is a new approximation method by propagating means and covariance through nonlinear transformations. The UKF is better than the EKF in terms of robustness and speed of convergence; the computational load in applying the UKF is similar to the EKF [25].

3.1. UKF approach

Fig. 4 is the flowchart of linear KF; and Fig. 5 is the flowchart of UKF theory. The unscented transfer is a mathematical tool for calculating the statistics of a random variable which undergoes a nonlinear transformation. The procedures of the algorithm can be organized as follows [18].

3.1.1. Sigma points calculation

The sigma points are a set of points which have mean and covariance equal to the given mean and covariance. The elements of sigma points are discrete in probability distribution. This distribution can be propagated exactly by applying the nonlinear function to each point.

In the paper, we use the symmetrically sampling method to pick up the sigma points. For the random variable state vector x, the mean is x- and the covariance of x is denoted as Px. The sigma points are chosen so that their mean and covariance are exactly equal to x- and Px respectively. Define χk as a set of 2N+1 sigma points [26]:

(5)
χ k = x - , k = 0 , x - + ( N + λ P x ) k , k = 1,2 , , N , x - - ( N + λ P x ) k - N , k = N + 1 , ,   2 N ,

where, N is the dimension of the state vector x; λ=α2N+κ-N is a scaling parameter, the constant α determines the spread of the sigma around x-, and is usually set to a small value (1α10-4)); κ0, it can make sure that the covariance matrix is positive definite. N+λPxk is the ith column of the matrix square root.

Fig. 4. Flowchart of KF

 Flowchart of KF

3.1.2. Nonlinear system function

These sigma vectors are propagated through the nonlinear system function:

(6)
χ k + 1 = χ k + f χ k , u k + 1 , w k + 1 Δ T ,

where, ΔT is the sampling time; f is the nonlinear system function; uk+1 is the system control input; wk+1 is the process noise.

3.1.3. Measurements update

(7)
Y k + 1 = h x k + 1 , u k + 1 , v k + 1 ,

where, h is measurement function; wk+1 is the measurement noise.

The mean of the state vector xk+1 and the measurement value yk+1 can be approximated by using a weighted sample mean of the sigma vectors:

(8)
x k + 1 | k i = 0 2 N ( W   i m ( χ k + 1 ) i ) ,       y k + 1 i = 0 2 N ( W   i m ( Y k + 1 ) i ) .

Fig. 5. Flowchart of UKF

Flowchart of UKF

3.1.4. Covariance update

The covariance of the state vector and the measurement vector can be updated by:

(9)
P     k + 1 | k x = i = 0 2 N [ W i c ( χ k + 1 ) i - x k + 1 χ k + 1 ) i - x k + 1 ' + Q ,
P     k + 1 y = i = 0 2 N [ W i c ( Y k + 1 ) i - y k + 1 Y k + 1 ) i - y k + 1 ' + R ,
P     k + 1 x y = i = 0 2 N [ W i c ( χ k + 1 ) i - x k + 1 Y k + 1 ) i - y k + 1 ' .

The weight can be updated by:

(10)
W i m = λ N + λ , 1 2 N + λ , , 1 2 N + λ 2 N ,    
(11)
W i c = λ N + λ + 1 - α 2 + β , 1 2 N + λ , , 1 2 N + λ 2 N ,

where, assuming that the system process noise and the measurement noise are Gaussian white noises, and the covariance of them are Q and R respectively. β considers the high order moment of the prior distribution, for Gaussian distribution, β= 2 is optimal.

3.1.5. Correction

(12)
K k + 1 = P k + 1 x y P k + 1 y - 1 ,       x k + 1 | k + 1 = x k + 1 | k + K k + 1 z k + 1 - y k + 1 ,
P k + 1 | k + 1 x = P k + 1 | k x - K k + 1 P k + 1 y ( K k + 1 ) ' ,

where, zk+1 is the measurement signal from the sensors.

3.2. Suspension state estimation

Define the state vector as x=[z˙s,z˙u,zs-zu,zu-q,zs,zu]T, and choose suspension stroke rate and stroke as observer measurements, i.e. h(x)=[z˙s-z˙u,zs-zu]T. The system dynamics can be rewritten in the format of state-space as:

(13)
x ˙ = f x + G w , y = h x + v ,

where:

f x = - F d m s - F s m s , F d m u + F s m u - k t m u x 4 , x 1 - x 2 , x 2 , x 1 , x 2 Τ ,
G = [ 0 0 0 - 1 0 0 ] Τ ,

represent the nonlinear system function and road disturbance matrix, respectively; h(x)=x1-x2,x3Τ represents the measurement outputs. w=q˙ is the process noise, which is the derivative of road disturbance and it is a white noise with covariance Q; v is the measurement noise assumed to be zero mean Gaussian white noise with covariance R.

Where, x1=z˙s – vertical velocity of sprung mass; x2=z˙u – vertical velocity of unsprung mass; x3=zs-zu – relative displacement of suspension; x4=zu-q – tire deflection; x5=zs – vertical displacement of sprung mass; x6=zu – vertical displacement of unsprung mass.

And the continuous-time nonlinear system Eq. (13) can be written as discrete-time nonlinear system function by using the first-order approximation, as shown in Eq. (14):

(14)
x k + 1 = x k + 1 + f x k + G w k Δ T , y k + 1 = h x k + 1 + v k + 1 .

Based on the above system function and the UKF estimation theory, the state observer of nonlinear suspension system is designed, as shown in Fig. 6. Here the road disturbance is considered as system process noise. So the designed observer is of low sensitivity to unknown random road disturbances.

Fig. 6. The flowchart of suspension state observer

The flowchart of suspension state observer

3.3. Sprung mass identification

The sprung mass often changes with the number of the passengers, loaded or unloaded conditions. Thus the resonant frequency of the vibration system will change [27], and the ride comfort will reduce greatly if the corresponding adjustment of the suspension cannot be made. So it is meaningful to identify the sprung mass in real-time. The sprung mass estimator is designed based on the linear KF. And combined with the state observer which designed in Section 3.2; the two estimators exchange the estimated results simultaneously, as shown in Fig. 7. Assuming that the sprung mass is constant when vehicle is driving, we choose θ=1/m^ as the system state and the sprung mass acceleration as measurement. The estimation function can be described as:

(15)
θ k + 1 = θ k + w 1 k , y k + 1 = H θ k + 1 + v 1 k , w 1 c o v a r i a n c e Q 1 , v 1 c o v a r i a n c e R 1 ,

where, H is the suspension force from UKF estimator; y=z¨s, is the measurement; w1k, v1k are the process noise and the measurement noise in the sprung mass identification system, respectively.

Fig. 7. The flowchart of suspension observer

 The flowchart of suspension observer

4. Simulation and analysis

The efficiency of designed observer is verified through numerical simulations on two different typical road profiles, the random road and speed bump road, respectively. The random road represents consistent excitations with a wide range of frequencies; and the speed bump road represents the discrete events of relatively short duration and high intensity. The simulation parameters of the quarter car are listed in Table 1.

Table 1. Suspension parameters of off-road vehicle

Description
Symbol
Value
Sprung mass (unloaded)
m s
2500 kg
Sprung mass (full load)
m s
3000 kg
Unsprung mass
m u
250 kg
Tire stiffness
k t
1501200 N/m

4.1. Case 1 random road excitation

In this section, the estimation efficiency is validated on random road excitation. The ISO has proposed a series of standards of road roughness classification by using the power spectral density (PSD) values (ISO 1982). The PSD of the random road excitation can be expressed as the following form:

(16)
G q n = G q n 0 n n 0 - w ,

where, n is the space frequency in m-1; n0 is the reference space frequency, n0= 0.01 m-1; Gqn0 is the road roughness coefficient, here, Gqn0= 256×10-6; w is frequency index, which reflects the frequency structure of the pavement, usually w= 2.

The model of road random excitation is built through integrated Gaussian white noise [28]. The equation of the random road in time domain can be expressed as:

(17)
q ˙ t = 2 π n 0 G q n 0 u W t - 2 π n 00 u q t ,

where, W(t) is the Gauss white noise; n00 is the low cut-off space frequency, n00= 0.001 m-1; u is the vehicle speed.

The irregular road profile is shown in Fig. 8. The vehicle speed is kept constant as 20 m/s. The estimation results are shown in Fig. 9-Fig. 11.

Fig. 9 and Fig. 10 are the comparison of UKF and EKF estimation results when vehicle is driving on the random road excitation. It is clearly indicated that the estimation results of UKF observer are more precise than those of EKF. Even though both of them have a short time delay, which is acceptable in engineering applications. Fig. 11 plots the results comparison of the sprung mass estimation when the vehicle is in unloaded and full load conditions, respectively. It can be found that both estimation methods work well.

Fig. 8. The random road profile

The random road profile

Fig. 9. Comparison of the sprung mass velocity estimation results

 Comparison of the sprung mass velocity estimation results

Fig. 10. Comparison of the unsprung mass velocity estimation results

 Comparison of the unsprung mass velocity estimation results

Fig. 11. Comparison of the sprung mass estimation results

Comparison of the sprung mass  estimation results

4.2. Case 2 speed bump road excitation

The speed bump road profile is shown in Fig. 12. The vehicle speed is kept constant as 5 m/s. The simulation results are shown in Figs. 13-14.

Fig. 13 and Fig. 14 are the comparison of suspension state estimation results when the vehicle is driving over a speed bump road. From Fig. 13, we can find that estimation result of sprung mass velocity from UKF observer is more precise compared with that from the EKF method.

Fig. 12. The speed bump

The speed bump

Fig. 13. Comparison of the sprung mass velocity results

 Comparison of the sprung mass  velocity results

Fig. 14. Comparison of the unsprung mass velocity results

Comparison of the unsprung mass  velocity results

Fig. 15. Errors analysis of estimation results

Errors analysis of estimation results

a) Random road condition

Errors analysis of estimation results

b) Speed bump road condition

4.3. Estimation errors analysis

In order to evaluate the performance of the proposed estimation algorithm, the index equation of the estimation performance is introduced:

(18)
P = R O b s e r v e r i R R e f e r e n c e i ,

where,  i is i-norm, and i is root mean square, maximum or minimum of the estimation results. It means that when the index value is more close to 1, the estimation result is more precise.

Fig. 15 plots the analysis results of the estimation errors of suspension states and the sprung mass. The estimation errors are analyzed with three indexes: root mean square, maximum and minimum of the estimation results according to Eq. (18). It can be seen that the performance indexes of the UKF’s state estimation are more close to 1 compared with the EKF estimation results; this means that the precision of the UKF is better than that of the EKF. The sprung mass estimation results are similar for the two estimation methods.

5. Conclusions

An observer with low sensitivity to road disturbance is designed based on UKF theory for nonlinear suspension systems. The UKF observer is successfully in estimating the state information of the suspension vibration with good accuracy. The efficiency of the proposed estimator is validated on two typical road excitations and in the condition of sprung mass variation. Compared with the EKF observer, the simulation results indicated that the performance of the designed UKF observer is more precise, and the sprung mass estimator is also proposed combined with the designed suspension state estimator. The performance of the sprung mass estimator is validated on random road excitations and the estimation results are very satisfying. The future work will be concentrating on developing the states observer for full vehicle suspension considering pitch and roll motion and validating the designed observer algorithms in real vehicle tests.

6. Acknowledgements

The present work is supported by the National Nature Science Foundation of China (Grant Nos. 51375046, 51205021).

References

  1. Koch G., Kloiber T. Driving state adaptive control of an active vehicle suspension system. IEEE Transactions on Control Systems Technology, Vol. 22, Issue 1, 2014, p. 44-57. [Search CrossRef]
  2. Kong Y., Zhao D., Yang B., et al. Non-fragile multi-objective static output feedback control of vehicle active suspension with time-delay. Vehicle System Dynamics, Vol. 52, Issue 7, 2014, p. 948-968. [Search CrossRef]
  3. Anubi O. M., Crane III C. D. Roll stabilization of road vehicles using a variable stiffness suspension system. Vehicle System Dynamics, Vol. 51, Issue 12, 2013, p. 1894-1917. [Search CrossRef]
  4. Tan U. X., Veluvolu K., Latt W. T., Shee C. Y., et al. Estimating displacement of periodic motion with inertial sensors. IEEE Sensor Journal, Vol. 8, 8, p. 1385-1388. [Search CrossRef]
  5. Heth EicTims Vehicle Active Suspension System Sensor Reduction. Ph.D. Thesis, the University of Texas at Austin, 2005. [Search CrossRef]
  6. Deshpande V. S., Mohan B., Shendge P. D., et al. Disturbance observer based sliding mode control of active suspension systems. Journal of Sound and Vibration, Vol. 333, Issue 11, 2014, p. 2281-2296. [Search CrossRef]
  7. Tang G., Li J., Ding C., et al. Sprung Mass Identification of Suspension in a Simplified Model. SAE Technical Paper, 2014-01-0051. [Search CrossRef]
  8. Rajamani R., Hedrick J. K. Adaptive observers for active automotive suspensions: theory and experiment. IEEE Transactions on Control Systems Technology, Vol. 3, Issue 1, 1995, p. 86-93. [Search CrossRef]
  9. Hernandez-Alcantara D., Amezquita-Brooks L., Morales-Menendez R. State-Observers for Semi-Active Suspension Control Applications with Low Sensitivity to Unknown Road Surfaces. SAE Technical Paper, 2014-01-0867. [Search CrossRef]
  10. Best M. C., Gordon T. J. Suspension system identification based on impulse-momentum equations. Vehicle System Dynamics, Vol. 29, Issue S1, 1995, p. 598-618. [Search CrossRef]
  11. Rajamani R., Hedrick J. K. Adaptive observers for active automotive suspensions: theory and experiment. IEEE Transactions on Control Systems Technology, Vol. 3, Issue 1, 1995, p. 86-93. [Search CrossRef]
  12. Jonathan Daniel Ziegenmeyer Estimation of Disturbance Inputs to a Tire Coupled Quarter-car Suspension Test Rig. Master Thesis, Virginia Polytechnic Institute and State University, 2007. [Search CrossRef]
  13. Rahman M., Rideout G. Using the lead vehicle as preview sensor in convoy vehicle active suspension control. Vehicle System Dynamics, Vol. 50, Issue 12, 2012, p. 1923-1948. [Search CrossRef]
  14. Zhang X., Guo K., et al. Adaptive real-time estimation on road disturbances properties considering load variation via vehicle vertical dynamics. Mathematical Problems in Engineering, 2013, p. 1-9. [Search CrossRef]
  15. Yamamoto A., Sugai H., Kanda R., Buma S. Preview Ride Comfort Control for Electric Active Suspension (eActive3). SAE Technical Paper 2014-01-0057. [Search CrossRef]
  16. Gerasimos Rigatos, Pierluigi Siano, Sesto Pessolano Design of active suspension control system with the use of Kalman filter-based disturbances estimator. Cybernetics and Physics, Vol. 1, Issue 4, 2012, p. 279-294. [Search CrossRef]
  17. Ren H., Chen S., Shim T., et al. Effective assessment of tyre-road friction coefficient using a hybrid estimator. Vehicle System Dynamics, Vol. 52, Issue 8, 2014, p. 1047-1065. [Search CrossRef]
  18. Ren H., Chen S., Liu G., et al. Vehicle state information estimation with the unscented Kalman filter. Advances in Mechanical Engineering, 2014, p. 1-11. [Search CrossRef]
  19. Lu Fan Study on Vehicle Vibration State Observation Algorithm Based on the Nonlinearity of Suspension. Ph.D. Thesis, Beijing Institute of Technology, 2014. [Search CrossRef]
  20. Chamseddine Abbas, Hassan Noura Control and sensor fault tolerance of vehicle active suspension. IEEE Transactions on Control Systems Technology, Vol. 16, Issue 3, 2008, p. 416-433. [Search CrossRef]
  21. Luque P., Mántaras D. A., Fidalgo E., et al. Tyre-road grip coefficient assessment, Part II: online estimation using instrumented vehicle, extended Kalman filter, and neural network. Vehicle System Dynamics, Vol. 51, Issue 12, 2013, p. 1872-1893. [Search CrossRef]
  22. Hamann H., Hedrick J. K., Rhode S., et al. Tire force estimation for a passenger vehicle with the unscented Kalman filter. IEEE Intelligent Vehicles Symposium Proceedings, 2014, p. 814-819. [Search CrossRef]
  23. Dixit R. K., Buckner G. D. Sliding mode observation and control for semiactive vehicle suspensions. Vehicle System Dynamics, Vol. 43, Issue 2, 2005, p. 83-105. [Search CrossRef]
  24. Eric Wan, Rudolph van der Merwe Kalman Filtering and Neural Networks. Chapter 7: The Unscented Kalman Filter. Wiley Publishing, 2001. [Search CrossRef]
  25. Kandepu R., et al. Applying the unscented Kalman filter for nonlinear state estimation. Journal Process Control, Vol. 18, 2008, p. 753-768. [Search CrossRef]
  26. Haykin Simon Kalman Filtering and Neural Networks. Wiley, New York, 2001. [Search CrossRef]
  27. Song X., Ahmadian M., Southward S., et al. An adaptive semiactive control algorithm for magnetorheological suspension systems. Journal of Vibration and Acoustics, Vol. 127, Issue 5, 2005, p. 493-502. [Search CrossRef]
  28. Wu Z. C., Chen S. Z., Yang L., Zhang B. Model of road roughness in time domain based on rational function. Transaction of Beijing Institute of Technology, Vol. 29, Issue 9, 2009, p. 795-798. [Search CrossRef]