Abstract
This paper presents a new approach to estimating suspension state information and parameter in realtime. 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.
1. Introduction
In recent years, there has been growing interest in the development of suspension system control strategy [13], such as active and semiactive 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 integratorbased 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 realtime; 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 [68].
The performance of the controllable suspension can be further improved if accurate state information could be obtained in realtime. 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 leastsquares (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 sensorbased 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 [1214] 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 realtime. 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. 1Nonlinear 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. ${m}_{u}$ is the unsprung mass which is supported by the tire modeled as linear spring with stiffness coefficient ${k}_{t}$. The displacement of the sprung mass and unsprung mass are denoted as ${z}_{s}$ and ${z}_{u}$ respectively; and the road profile is denoted as $q$. So the dynamic equations of the nonlinear suspension can be described as:
where, ${F}_{s}$ is the air spring force; which is the function of damper stroke, as shown in Fig. 2; ${F}_{d}$ 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 3D lookup table model is created to describe the nonlinear characteristics of the adjustable damper.
Fig. 2The characteristic of the air spring force
Fig. 3Damping forcevelocity 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 suboptimal 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 $\mathbf{x}$, the mean is $\stackrel{}{\mathbf{x}}$ and the covariance of $\mathbf{x}$ is denoted as ${\mathbf{P}}^{\mathbf{x}}$. The sigma points are chosen so that their mean and covariance are exactly equal to $\stackrel{}{\mathbf{x}}$ and ${\mathbf{P}}^{\mathbf{x}}$ respectively. Define ${\chi}_{k}$ as a set of $2N+1$ sigma points [26]:
where, $N$ is the dimension of the state vector $\mathbf{x}$; $\lambda ={\alpha}^{2}\left(N+\kappa \right)N$ is a scaling parameter, the constant $\alpha $ determines the spread of the sigma around $\stackrel{}{\mathbf{x}}$, and is usually set to a small value (1$\le \alpha \le $10^{4})); $\kappa \ge 0\text{,}$ it can make sure that the covariance matrix is positive definite. ${\left(\sqrt{\left(N+\lambda \right){\mathbf{P}}^{\mathbf{x}}}\right)}_{k}$ is the $i$th column of the matrix square root.
Fig. 4Flowchart of KF
3.1.2. Nonlinear system function
These sigma vectors are propagated through the nonlinear system function:
where, $\mathrm{\Delta}T$ is the sampling time; $f$ is the nonlinear system function; ${u}_{k+1}$ is the system control input; ${\mathbf{w}}_{k+1}$ is the process noise.
3.1.3. Measurements update
where, $h$ is measurement function; ${\mathbf{w}}_{k+1}$ is the measurement noise.
The mean of the state vector ${\mathbf{x}}_{k+1}$ and the measurement value ${\mathbf{y}}_{k+1}$ can be approximated by using a weighted sample mean of the sigma vectors:
Fig. 5Flowchart of UKF
3.1.4. Covariance update
The covariance of the state vector and the measurement vector can be updated by:
${\mathbf{P}}_{{\mathrm{}}^{{\mathrm{}}_{k+1}}}^{\mathbf{y}}=\sum _{i=0}^{2N}\left[{\mathbf{W}}_{i}^{c}\right(\left({Y}_{k+1}{)}_{i}{\mathbf{y}}_{k+1}\right)\left({\left({\mathbf{Y}}_{k+1}{)}_{i}{\mathbf{y}}_{k+1}\right)}^{\mathrm{\text{'}}}\right]+\mathbf{R},$
${\mathbf{P}}_{{\mathrm{}}^{{\mathrm{}}_{k+1}}}^{\mathbf{x}\mathbf{y}}=\sum _{i=0}^{2N}\left[{\mathbf{W}}_{i}^{c}\right(\left({\chi}_{k+1}{)}_{i}{\mathbf{x}}_{k+1}\right)\left({\left({\mathbf{Y}}_{k+1}{)}_{i}{\mathbf{y}}_{k+1}\right)}^{\mathrm{\text{'}}}\right].$
The weight can be updated by:
where, assuming that the system process noise and the measurement noise are Gaussian white noises, and the covariance of them are $\mathbf{Q}$ and $\mathbf{R}$ respectively. $\beta $ considers the high order moment of the prior distribution, for Gaussian distribution, $\beta =\mathrm{}$2 is optimal.
3.1.5. Correction
${\mathbf{P}}_{k+1k+1}^{x}={\mathbf{P}}_{k+1k}^{x}{\mathbf{K}}_{k+1}{\mathbf{P}}_{k+1}^{y}({\mathbf{K}}_{k+1}{)}^{\mathrm{\text{'}}},$
where, ${\mathbf{z}}_{k+1}$ is the measurement signal from the sensors.
3.2. Suspension state estimation
Define the state vector as $\mathbf{x}=[{\dot{z}}_{s},{\dot{z}}_{u},{z}_{s}{z}_{u},{z}_{u}q,{z}_{s},{z}_{u}{]}^{T}$, and choose suspension stroke rate and stroke as observer measurements, i.e. $h\left(\mathbf{x}\right)=[{\dot{z}}_{s}{\dot{z}}_{u},{z}_{s}{z}_{u}{]}^{T}$. The system dynamics can be rewritten in the format of statespace as:
where:
$\mathbf{G}=[\begin{array}{llllll}0& 0& 0& 1& 0& 0\end{array}{]}^{{\rm T}},$
represent the nonlinear system function and road disturbance matrix, respectively; $h\left(\mathbf{x}\right)={\left[\begin{array}{ll}{x}_{1}{x}_{2},& {x}_{3}\end{array}\right]}^{{\rm T}}$ represents the measurement outputs. $\mathbf{w}=\dot{q}$ is the process noise, which is the derivative of road disturbance and it is a white noise with covariance $\mathbf{Q}$; $v$ is the measurement noise assumed to be zero mean Gaussian white noise with covariance $\mathbf{R}$.
Where, ${x}_{1}={\dot{z}}_{s}$ – vertical velocity of sprung mass; ${x}_{2}={\dot{z}}_{u}$ – vertical velocity of unsprung mass; ${x}_{3}={z}_{s}{z}_{u}$ – relative displacement of suspension; ${x}_{4}={z}_{u}q$ – tire deflection; ${x}_{5}={z}_{s}$ – vertical displacement of sprung mass; ${x}_{6}={z}_{u}$ – vertical displacement of unsprung mass.
And the continuoustime nonlinear system Eq. (13) can be written as discretetime nonlinear system function by using the firstorder approximation, as shown in Eq. (14):
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. 6The 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 realtime. 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 $\theta =1/\widehat{m}$ as the system state and the sprung mass acceleration as measurement. The estimation function can be described as:
where, $\mathbf{H}$ is the suspension force from UKF estimator; $\mathbf{y}={\ddot{z}}_{s}$, is the measurement; ${\mathbf{w}}_{1k}$, ${\mathbf{v}}_{1k}$ are the process noise and the measurement noise in the sprung mass identification system, respectively.
Fig. 7The 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 1Suspension parameters of offroad 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:
where, $n$ is the space frequency in m^{1}; ${n}_{0}$ is the reference space frequency, ${n}_{0}=\mathrm{}$0.01 m^{1}; ${G}_{q}\left({n}_{0}\right)$ is the road roughness coefficient, here, ${G}_{q}\left({n}_{0}\right)=\mathrm{}$256×10^{6}; $w$ is frequency index, which reflects the frequency structure of the pavement, usually $w=\mathrm{}$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:
where, $W\left(t\right)$ is the Gauss white noise; ${n}_{00}$ is the low cutoff space frequency, ${n}_{00}=\mathrm{}$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. 9Fig. 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. 8The random road profile
Fig. 9Comparison of the sprung mass velocity estimation results
Fig. 10Comparison of the unsprung mass velocity estimation results
Fig. 11Comparison 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. 1314.
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. 12The speed bump
Fig. 13Comparison of the sprung mass velocity results
Fig. 14Comparison of the unsprung mass velocity results
Fig. 15Errors analysis of estimation results
a) Random road condition
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:
where, ${\Vert \Vert}_{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

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. 4457.

Kong Y., Zhao D., Yang B., et al. Nonfragile multiobjective static output feedback control of vehicle active suspension with timedelay. Vehicle System Dynamics, Vol. 52, Issue 7, 2014, p. 948968.

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. 18941917.

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. 13851388.

Heth EicTims Vehicle Active Suspension System Sensor Reduction. Ph.D. Thesis, the University of Texas at Austin, 2005.

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. 22812296.

Tang G., Li J., Ding C., et al. Sprung Mass Identification of Suspension in a Simplified Model. SAE Technical Paper, 2014010051.

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. 8693.

HernandezAlcantara D., AmezquitaBrooks L., MoralesMenendez R. StateObservers for SemiActive Suspension Control Applications with Low Sensitivity to Unknown Road Surfaces. SAE Technical Paper, 2014010867.

Best M. C., Gordon T. J. Suspension system identification based on impulsemomentum equations. Vehicle System Dynamics, Vol. 29, Issue S1, 1995, p. 598618.

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. 8693.

Jonathan Daniel Ziegenmeyer Estimation of Disturbance Inputs to a Tire Coupled Quartercar Suspension Test Rig. Master Thesis, Virginia Polytechnic Institute and State University, 2007.

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. 19231948.

Zhang X., Guo K., et al. Adaptive realtime estimation on road disturbances properties considering load variation via vehicle vertical dynamics. Mathematical Problems in Engineering, 2013, p. 19.

Yamamoto A., Sugai H., Kanda R., Buma S. Preview Ride Comfort Control for Electric Active Suspension (eActive3). SAE Technical Paper 2014010057.

Gerasimos Rigatos, Pierluigi Siano, Sesto Pessolano Design of active suspension control system with the use of Kalman filterbased disturbances estimator. Cybernetics and Physics, Vol. 1, Issue 4, 2012, p. 279294.

Ren H., Chen S., Shim T., et al. Effective assessment of tyreroad friction coefficient using a hybrid estimator. Vehicle System Dynamics, Vol. 52, Issue 8, 2014, p. 10471065.

Ren H., Chen S., Liu G., et al. Vehicle state information estimation with the unscented Kalman filter. Advances in Mechanical Engineering, 2014, p. 111.

Lu Fan Study on Vehicle Vibration State Observation Algorithm Based on the Nonlinearity of Suspension. Ph.D. Thesis, Beijing Institute of Technology, 2014.

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. 416433.

Luque P., Mántaras D. A., Fidalgo E., et al. Tyreroad 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. 18721893.

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. 814819.

Dixit R. K., Buckner G. D. Sliding mode observation and control for semiactive vehicle suspensions. Vehicle System Dynamics, Vol. 43, Issue 2, 2005, p. 83105.

Eric Wan, Rudolph van der Merwe Kalman Filtering and Neural Networks. Chapter 7: The Unscented Kalman Filter. Wiley Publishing, 2001.

Kandepu R., et al. Applying the unscented Kalman ﬁlter for nonlinear state estimation. Journal Process Control, Vol. 18, 2008, p. 753768.

Haykin Simon Kalman Filtering and Neural Networks. Wiley, New York, 2001.

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. 493502.

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. 795798.