Abstract

The estimation of road excitation profile is important for evaluation of vehicle stability and vehicle suspension performance for autonomous vehicle control systems. In this work, the nonlinear dynamics of the active automotive system that is excited by the unknown road excitation profile are considered for modeling. To address the issue of estimation of road profile, we develop an adaptive supertwisting observer for state and unknown road profile estimation. Under Lipschitz conditions for the nonlinear functions, the convergence of the estimation error is proven. Simulation results with Ford Fiesta MK2 demonstrate the effectiveness of the proposed observer for state and unknown input estimation for nonlinear active suspension system.

1. Introduction

The automotive vehicle suspension dynamics contribute significantly in evaluating the effective performance with regard to passenger comfort, road handling, and stability of the vehicle [1]. The design of suspension systems evolved from passive suspensions to the active suspension system [2] adding more control capabilities. Suspension dynamics predominantly provide information regarding vertical stability [2, 3] of the vehicle. Effective analysis of the suspension performance provides information regarding the vertical load acting on the vehicle, a critical component in determining the effective tractive force [3]. In such scenarios, the suspension system dynamics that replicate the behavior of mass spring damper system [3] are complex to analyze when nonlinear behavior of the spring and the damper systems are considered. Road profile that replicates the randomness of road surface in form of cleats or troughs affects the suspension system performance [4]. The different levels of road excitation necessitate the continuous regulation of damping force generated by the suspension for maintaining the stability of the vehicle. As the vehicle operating range varies, analysis of the nonlinear dynamics of the suspension system excited by the road profile is one of the major domains of research for suspension systems.

Road profile is a critical parameter that results in undesirable vertical vibrations for the vehicle if it is not compensated by an adequate control effort in the suspension. These unfavorable vibrations result in dynamic variations in the vertical load of the vehicle affecting its stability. The design of suitable controllers to compensate for these variations is dependent on the effective measurement of the road profile. For experimental purposes, expensive instruments called profilographs are used for measuring the road profile [4, 5]. In [6], a real time conditioning algorithm was designed to measure the road profile based on measurement of the vertical acceleration. As an alternative to these expensive instrumentation and sensor technology that are affected by noise, estimation of the road profile by use of observers has been an important issue.

In [7], a Kalman filter-based estimation of the road profile measured using a road profilometer was performed. The designed estimation worked on the vertical dynamics of an active suspension system and was experimentally validated. In [8], a neural network-based approach was adopted to estimate the random road profile. This work employed accelerometers to estimate the road profile modeled as a function of road roughness coefficient. In [9], a minimum order observer was designed for a linearized model of passive suspension dynamics to estimate the road profile.

Sliding mode theory [10, 11] has evolved over time as an effective tool for estimation of unknown inputs in control domain. In [12], a first order sliding mode observer was designed to estimate the states of the suspension system under the influence of different road profiles. To remove the chattering effect that is inherent in first order sliding mode [1317], higher order sliding mode theory was developed. Recently, HOSM observers have been popular for state and unknown input estimation in uncertain nonlinear systems [1820]. For systems with relative degree one, the super-twisting algorithm (STA) served as an ideal tool for the estimation of unknown inputs [21]. In [2224], the STA based observer was used to estimate the road profile acting as unknown input to the system. By measurement of the vertical velocities, the road profile and the tire forces were estimated. All these works [79] did not consider the nonlinear dynamics of the suspension system. Other works [2224] did not consider the randomness of the road profile.

To address these issues, we consider the nonlinear dynamics of the active suspension system for a quarter vehicle excited by a random road profile. The suspension dynamics considered in this paper effectively replicate the nonlinear behavior of the spring and damper of the suspension. The road excitation profile is considered as an unknown input in this work and is estimated with an adaptive STA observer [25]. For analysis, the road roughness values are based on power spectral density (PSD) values as proposed by International Organization for Standardization (ISO) [4, 5]. For the design of the observer, under the rank conditions for the output matrix, the system is then partitioned into two subsystems where the unknown input appears in one subsystem. For the subsystem affected by the unknown input, an adaptive STA based observer is then designed to ensure the stability of the error dynamics of the subsystem in finite time. For the subsystem without unknown inputs, a nonlinear observer is designed under Lipschitz conditions to ensure the stability of the system in sliding mode. The application of the proposed method to the modelled vehicle dynamics is validated through simulations.

Throughout this paper, denotes the maximum eigenvalue of matrix A, denotes the 2-norm of a matrix A, and represents the minimum singular value of matrix A. For any vector and any scalar , we denote

2. Modeling Active Suspension Dynamics

The active suspension system in vehicles incorporates an active controlled force actuator instead of the shock absorber generally found in a passive suspension system. The nonlinear dynamics that govern the active suspension system are given as [26] where is the sprung mass or the vehicle mass, is the unsprung mass or the wheel mass, is the sprung mass displacement, is the unsprung mass displacement, is the controlled actuator force, is the tire stiffness, and is the road excitation profile. The nonlinear damping force, , and spring force, , for the suspension dynamics can be described as [26] where is the linear spring stiffness constant, is the linear damper constant, is the nonlinear spring stiffness, and is the nonlinear damping constant. The motion of the vehicle over a bump that restricts the wheel travel within a given range and prevents contact between the tyre and the vehicle body is effectively modeled by the nonlinear spring force, . Similarly for the damper, the damping force generated while the wheel traverses in vertical direction owing to road profile is a nonlinear effect. This nonlinear effect is well approximated by the nonlinear dynamics as depicted in (5). In this work, the damping force provided by the tyre which is very complex to model has been neglected. The active suspension dynamics for a quarter wheel vehicle model together with the modeling parameters are shown in Figure 1. The active suspension dynamics are affected by the road profile requiring control of the effective damping force needed to be provided by the actuator for good handling of the vehicle [1], ride performance, and road stability. Road profile is often modeled as a sinusoidal disturbance or a trapezoidal disturbance to identify with crests/trough or cleats that appear on practical roads. This type of modeling however does not represent the typical roughness profile of roads and the resulting effects it has on the suspension. The standards of road roughness according to ISO [4, 5] can be classified into different road classes as shown in Table 1.

The integrated dynamics (2)–(5) for the active suspension system can be represented in state space as With where In the modeled system dynamics, the active actuator control force is the control input for the system. In the system dynamics (6), the displacement of the sprung mass, , and velocity of the unsprung mass, , are considered as outputs, and the output matrix can be defined by The unknown input for the system is the road profile denoted as . The measurement of road profile is an extremely complex task that requires the use of complex measuring instruments such as profilographs [4, 5] that are expensive and impractical. Hence our focus is on the development of an approach to estimate the random road profile, , for the active suspension systems.

3. Observer Design

In this section, we discuss the design of the observer for the active suspension system. A combination of nonlinear Lipschitz observer and adaptive super-twisting observer is employed. To facilitate the design of the observer, the following assumptions are required.

Assumption 1. All invariant zeros of the triple must lie in the left half plane and rank = rank .

Assumption 2. The nonlinear functions in satisfiy the Lipschitz conditions.

Assumption 3. The function and its first derivative are bounded.

Assumption 4. The control input is bounded and the system is assumed to be bounded input bounded state stable (BIBS).

For the dynamics defined in (6), it can be easily verified that the . Further, the triple does not contain any invariant zeros that lie on the right hand plane. The nonlinear function in the system dynamics (6) can be divided into and to analyze the Lipschitz continuity as follows: With , that is, and , that is, representing the physical states of practically realizable automotive suspension system, the function can be determined to be locally Lipschitz with a Lipschitz constant obtained as . Similarly for the nonlinear function , with , that is, and , that is, bounded, local Lipschitz continuity can be established. The Lipschitz constant for is obtained as . For the active automotive suspension system modeled in (6), by use of the actuator control force, , it can be asserted that such that the function maintains its Lipschitz continuity with a Lipschitz constant, .

In the modeled system dynamics (6), the unknown input is the road profile that is considered as function of road roughness coefficient and other physical parameters relating to the road conditions which are bounded. It can be thus deduced from the dynamics of the road profile model that the road excitation profile and its derivative are both bounded. For the modeled active suspension dynamics (6), the control input to the system is the actuator force, , that is bounded.

In order to design the combined observer, the original system dynamics will be divided into two subsystems, such that one of the subsystems will be free from unknown inputs. With the Assumptions 14 for the active suspension system, (6) being satisfied, we can directly partition system (6) into two subsystems and as follows: where

With the system (6) partitioned as above, the objective is to design an adaptive STA based observer to estimate the states and unknown input for the -subsystem (11) and a nonlinear Lipschitz observer (NLO) to estimate the states for the -subsystem (12). The overview of the design is shown in Figure 2.

The estimation error can be defined as where is the observed state and and are the errors for the subsystems and .

3.1. Adaptive Super-Twisting Observer Design for S1-Subsystem

For -subsystem (11) satisfying the above assumptions, the following observer based on the adaptive STA can be designed to estimate the states and the unknown input: where is the robust sliding term based on the adaptive STA [25] and defined as The adaptive gains , in (16) are designed as where , , , and are positive constants. The parameter is an arbitrary small positive constant.

To establish the convergence of the observer dynamics, the error dynamics (14) can be obtained as , which serves as the sliding surface for the designed adaptive STA based observer. The objective of the designed observer is to ensure that the error converges to zero and to reconstruct the unknown road excitation profile, , from the robust term (16).

Theorem 5. For system (11) satisfying the Assumptions 14, the observer system (15) with the robust term (16) will ensure that the error dynamics () will converge to zero in finite time.

Proof. The error dynamics of the system (11) can be obtained from (11) and (15) as where includes the perturbation terms. The matrices and are known, and hence they are bounded. As system satisfies Assumptions 24, the boundedness of the nonlinear function and the unknown input, , can be easily established. For the obtained error dynamics (18), it can be proved that is locally bounded by a constant, as is twice differentiable on a compact set. This is not restrictive as the active suspension system dynamics are bounded at least locally. It will be shown later that the subsystem is asymptotically stable. Further, as the system is free from unknown inputs, under the Assumptions 24, the subsystem boundedness can be established. Based on the above arguments, the boundedness of the perturbation is obtained as where is a constant (not necessarily known). With the perturbation terms , satisfying the condition (19) required for the adaptive STA (15), the convergence of the error dynamics (18) can now be proved with the following Lyapunov function: where and , , , and are positive constants and is a positive definite matrix. Similar to the results in [25] with and satisfying (17), can be shown to be a negative definite and the error converges to zero in finite time. The sliding surface is thus reached in finite time and maintained thereafter.

3.2. Nonlinear Lipschitz Observer for S2-Subsystem

For subsystem (12), a NLO is designed as follows to estimate the states of the system: where the feedback is to be discussed in Theorem 6 later. The error dynamics (14) of the subsystem can be obtained as The following theorem establishes the stability of the -subsystem.

Theorem 6. For system (12) satisfying the Assumptions 14, the observer (21) ensures that the state estimation error is asymptotically stable provided that the gain L satisfies where is the Lipschitz constant for (Assumption 2) and Q is a positive definite matrix.

Proof. With the convergence of the subsystem error () to zero in the sliding mode, the error dynamics (14) can be written as With the system satisfying Assumption 2, the Lipschitz constant for is evaluated as . With the choice of the Lyapunov function as and differentiating with respect to time, one has From the above results, its can be deduced that In the sliding mode as , we have It can thus be written as follows: Further, one can obtain where is satisfied. If the design of the feedback gain, , is such that (23) is satisfied, then . The error dynamics will thus be asymptotically stable.

Remark 7. Equation (23) can be written as an algebraic Riccati equation in the following form: for some . The following condition [27, 28] ensures the asymptotic stability of the system (24): where represents the minimum singular value of a matrix. If the above condition (23) is satisfied and if there exists a stable matrix, then there exists a symmetric positive definite (SPD) solution for the Riccati equation (30).

3.3. Estimation of Unknown Input: Road Excitation Profile

In the sliding mode, with the error converging to zero in finite time , the equivalent control [29] can be obtained from (18) as follows: As in finite time, and the nonlinearities satisfy Lipschitz assumptions, we have from (18) the following: The unknown road excitation profile when can be thus obtained as

Remark 8. The design of the adaptive STA observer and the NLO considered in this work can be easily extended for estimation of multiple unknown inputs. The general class of nonlinear systems that is similar to (6) is represented by where , , , , with are the unknown inputs/uncertainties under similar assumptions, general class of nonlinear systems is defined in (35); a linear transformation [10] can always be employed to obtain the required structure for the design of the observers.

4. Results

For the performance evaluation of the proposed observer, we select the following active suspension system parameters of a Ford Fiesta MK2 [30] vehicle:  kg,  kg,  N/m,  Ns/m, and  N/m. The nonlinear spring stiffness, , and damping constant, , values are taken as 10 of the original linear values and ,  respectively. To design the adaptive STA, we choose the gains as , , , , and . The initial conditions for the plant and the observer were chosen as () = and = . The feedback gain for the subsystem (12) and the positive definite matrix, , satisfying (23) were computed as The Lipschitz constant for subsystem (12) was evaluated as . The road profile model [31] can be obtained as where is the vehicle longitudinal velocity, is the road roughness coefficient, is the reference space frequency, and is the Gaussian white noise. With changes in road roughness coefficient keeping the longitudinal velocity of motion fixed, different excitation profiles can be obtained. The increase in roughness coefficient leads to poorer road conditions. For simulation purposes, the class C road profile was considered as an unknown input. The vehicle was considered to be travelling at a speed of 30 Km/hr with road roughness coefficient being considered as shown in Table 1. The generated road profile is shown in Figure 3. The simulation results obtained for state estimation are shown in Figure 4. The system dynamics are nonlinear and affected by the road profile, , which is a function of white Gaussian noise. It can be deduced that the estimation of the states even under the effect of is good. The norm of the estimation error of the states is shown in Figure 5. In Figure 6, the unknown road excitation profile estimated with (34) is shown. A smooth estimation of the unknown road profile is obtained without any low-pass filtering.

5. Conclusions

In this work, an adaptive super-twisting observer was proposed for state and unknown input estimation for the active suspension system. The paper considered the nonlinear model of the active suspension system excited by the random road profile as an unknown input. Under the Lipschitz condition for the nonlinear functions, the convergence of the system errors is proven. The proposed adaptive super-twisting observer accurately estimates the road excitation profile for an average road without the use of low-pass filter.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. 2011-0023999).