Sliding Mode Control for Active Suspension System with Data Acquisition Delay

This paper addresses the problem of control of an active suspension system accomplished using a computer. Delay in the states due to the acquisition and transmission of data from sensors to the controller is taken into account. The proposed control strategy uses state predictors along with sliding mode control technique. Two approaches are made: a continuous-time and a discrete-time control. The proposed designs, continuous-time and discrete-time, are applied to the active suspension module simulator from Quanser. Results from computer simulations and experimental tests are analyzed to show the effectiveness of the proposed control strategy.


Introduction
There are different approaches to control active suspension systems:  ∞ control [1,2], proportional integral sliding mode control [3,4], linear quadratic regulator [5], fuzzy controllers [6], artificial neural networks [7], and so forth.In this paper focus is given to problems that can arise when computers do the control.In this context, the sampling time involved in the process and possible delays sampled signals may damage the control performance.These signal delays can arise in Networked Control Systems (NCS) where the communication between the elements of the control loop is done by a nonideal shared digital network [8,9].Wireless Networked Control Systems (WNCS) usually show delays of many sampling periods and cannot be ignored [10].
The main advantage of Variable Structure Control with Sliding Mode (VSC/SM) is its robustness with respect to matched uncertainties, representing disturbances in control input, and parametric uncertainties or even system modeling uncertainties, provided that they belong to the image of the control input matrix [11,12].When the control signal is generated by microprocessors, besides sampling time, the system performance is greatly compromised when delays in both data acquisition and data transmission occur.In VSC/SM, this issue is even more important.By using a high speed switching control law in order to take the states trajectory to a sliding surface dependent on the current states, if the states used are delayed, the control law may not direct the states to this surface, which can also generate performance loss or even lead system to instability [13,14].The damage caused by delays to sliding mode control motivates several studies [15][16][17][18].
This paper proposes a strategy that uses state predictors along with VSC/SM in order to perform control of an active suspension system in presence of input control disturbance and data acquisition delay of many sampling periods.With purpose of comparison, it proposes two designs: one continuous-time control and one discrete-time control.In both cases, the control signal is generated from a digital computer.

Mathematical Problems in Engineering
For the purpose of illustration, the results from proposed controls, both continuous and discrete, are compared with the results of a linear quadratic regulator [19,20] proposed by Quanser, the manufacturer of the active suspension system used in the practice tests.

Control Strategies Proposed
The objective of variable structure control with sliding mode, by means of a high speed switching control law, is to lead the state trajectory of the system to a surface in the state space and ensure that it holds to this surface.The surface is designed so that the system presents a desired dynamic when in sliding mode.Thus, the VSC/SM design is made up of two steps: the design of the switching surface and the design of the control law [12].
This section presents two controllers designs for state delayed systems: one uses continuous-time equations and the other uses discrete-time equations; that is, the sampling period of the computer is taken into account.In both projects, the proposed control strategy is to use a state predictor to estimate the current states, which forms the VSC/SM law.In this work the sampling period and the data acquisition delay are considered as known.Also, all the states of the plant are considered available.

Design of Sliding
Surface.Consider a continuous-time system in the regular form given by (1), where 1) , and  22 ∈ R × are constant matrices: The linear sliding surface proposed is given by (2), where  1 ∈ R ×(−) ,  2 ∈ R × , and () ∈ R  : Under these conditions, the dynamics of the states in sliding mode is described by (3) [12].Note that the dynamics of the system, which previously had dimension , is replaced by dimension  −  when in sliding mode: The dynamics described by (3) has state feedback structure; thus one can use standard techniques to find the gain, as pole placement or quadratic optimal control [19].
Note that the sliding surface design in discrete-time is made in an analogous manner to the continuous-time one, and when the original system is not in the regular form, a nonsingular coordinate transformation can be used to put it in this form [21].

Continuous-Time Control
Design.Consider the following system:

Plant
Data acquisition where () ∈ R  is the state vector; () ∈ R  is the control vector; the matrices  and  are constants and with appropriate dimensions.
In this paper, it is considered that the control signals are generated by microprocessor and the states feedback is done by a data acquisition system, including analog/digital converters.Time delays may occur in data acquisition and data transmission; these delays are considered known and constant.Figure 1 is a schematic in which is illustrated this condition through emulation, that is, using small sampling periods so that the system can be approximated as a continuous-time.
If the control signals are generated from the delayed-time sampled states ( − ℎ), then the controller will generate the control signal ( − ℎ); that is, controlled plant (4) that was free of delay becomes a plant with delay described by (5): Generally delays adversely affect the performance of control systems, whatever control method is used.In particular, under the approach of the VSC/SM, eliminating the data acquisition delay becomes important, because it is very sensitive to any kind of delay.The estimator presented in (6) [22] will be used in order to estimate the current state: where   () is the estimate of the current state vector (( − ℎ + ℎ) = ()).
The dynamics of the predictive state vector   () is free of delay as can be viewed by (8), from the time derivative of   (), as follows: Mathematical Problems in Engineering 3 Substituting (5) into (7), one has that shows that the state predictive dynamics is free of delay.Now, consider an uncertain system undergoes delay during data acquisition; in this case the dynamics of the system is represented by Applying the continuous predictor proposed, one has the following predictive state: Again, its dynamics are calculated as But replacing (9) in (11) ẋ  () =  ℎ  ( − ℎ) +  ℎ  ( − ℎ) +  ℎ Δ () Therefore, in this case, the initial uncertainty of the system Δ() becomes  ℎ Δ() and may influence the system dynamics in sliding even though it was originally a matched uncertainty.Thus, even when the system presents uncertainties in the plant, the dynamics of predictive state remain free of delays.
Therewith, one can develop the conventional VSC/SM design from predictive state vector   ().The sliding surface is given by where  ∈ R × is one constant matrix, which establishes the system dynamics in sliding mode.The purpose of the VSC/SM is to find a control law which ensures the existence of the sliding mode.A control law that achieves this goal consists of one continuous part along with one discontinuous part, as proposed in Decarlo et al. [12]: where  eq () is the continuous part called equivalent control and   () is the discontinuous part.The equivalent control is that which determines the dynamics of the system in sliding mode.It is found using the condition described by Given that in sliding mode   () = 0 and considering the product  nonsingular, we have The discontinuous control part (  ()) is responsible for bringing and maintaining the system state trajectory to the sliding surface.It is found using a generalized Lyapunov function which ensures the convergence of state trajectory to the sliding surface.Consider the following Lyapunov function candidate: The time derivative of () is negative definite if In view of the control given in ( 14) and ( 17) it follows that A control law that satisfies (19) considering (20) and assuming  =  is given in where   () is given by Control law (22) suffers from the chattering problem.A modification to reduce this effect is proposed in Spurgeon and Davies [23] and presented in This paper will use only the expression continuous sliding mode control (CSMC) to refer to the control described by ( 14), ( 17), (21), and (23).

Discrete-Time Control Design.
In the control scheme shown in Figure 1, it is assumed now that the signals are sampled in the sampling period .Thus, the dynamic representation of the controlled plant in discrete-time state space is given by where () is the state vector sampled at time  and V() is the control vector created at the same time.The matrices Φ ∈ R × and Γ ∈ R × are constants and they depend on the sampling period .Now, the time delay in the system states caused by data acquisition and/or data transmission is taken into account.This delay is given by ℎ = .Similar to the continuous case, if the control signals are generated from the delayed sampled states ( − ) then the control law will be V( − ); that is, controlled plant (24), that was free of delay, becomes a plant with control delayed in  samples.So the representation of the system dynamics sampled with delay is given by In order to estimate the current states sampled, an estimator in the form of (26), adapted from Xia et al. [24], is used: where   () is the current vector state estimate (( −  + ) = ()).
The dynamics of the predictive state is free of delay as can be seen from ( 28) by calculating   ( + 1), as follows: Substituting ( 25) into (27), it comes to If the digital system with data acquisition delay considered is uncertain, as described by ( 29) the predictive state is given by Calculating   ( + 1) it follows that Using ( 29) one reaches Thus, the dynamics of predictive state is free from delays even in the presence of uncertainties in the plant.Note that Δ() has been multiplied by Φ  ; that is, assuming that Δ() was initially a matched uncertainty it lets this condition bypass through the predictor.
Therefore, as the predictive state dynamics is free of delay one can develop the VSC/SM design from the predictive state vector   ().A discrete sliding surface is given by where  ∈ R × is a constant matrix, which establishes the system's dynamic in sliding mode.
The discrete control law that establishes a sliding mode is proposed by Garcia et al. [25]: where V eq () is the equivalent control that determines the dynamics of the system in sliding mode and V  () is responsible for taking the system to the sliding mode.The equivalent control is found when it is guaranteed that the state trajectory remains on the surface, that is, Replacing the system (25) into (35), considering Γ nonsingular, and recalling that in sliding mode V  () = 0, one finds where  denotes the identity matrix with appropriate dimensions.
The control V  () is chosen to ensure the convergence of the state trajectory to the sliding surface.This problem is similar to a stability problem, and one can use the discretetime Lyapunov's second method to solve this.
This paper utilizes a state space linear model of a quarter of vehicle corresponding to Quanser's bench active suspension system.Its diagram is shown in Figure 2 [20].
The active suspension system can be modeled as a double mass-spring-damper system [20].Thus, the two system inputs are control signal   () and the derivative of road surface ż  ().
The coordinate   () represents the displacement of the tire, which has mass   , and the coordinate   () represents the displacement of the vehicle body, which has mass   .The movements are related to the movement imposed by the road surface on which the vehicle is traveling   ().
The motion equations of the system are obtained using the free body diagram method and can be described in the state space in the form of (46).The Quanser Innovate Educate [20] presents a detailed deduction of this model: In (46), the  1 is the related matrix to the road surface derivative ( ż  ()) and  2 is the related matrix to the control signal (  ()).
The four states considered accessible are defined in The first state is the suspension travel.The second state is the vehicle body vertical velocity.The third state is the tire deflection.The fourth state is the wheel vertical velocity.The first input, (), is the road surface velocity ( ż  ()).The second input, (), is the control action (  ()).

Application of the Proposed Control to Active Suspension
with Data Acquisition Delay.The numerical values of the state space matrices ,  1 , and  2 , in (46) and (48), are calculated from data provided by the manufacturer's manual [27], shown in Table 1.
The calculations involved in the design of the controllers were performed using Matlab language and based on the nominal values given by the manufacturer.However, in Section 4 uncertainties were included in order to demonstrate robustness, shown in Section 2.  For the continuous-time control design, ( 14), ( 17), (21), and ( 23), the  matrix from sliding surface (13)  The sampling period for this control was chosen to be 1 ms so that the emulation presented good results.
The practical experiments were performed in Quanser's bench active suspension system [27] using Quanser's AMPAQ current amplifier [29] and QPID data acquisition board [30], according to the scheme shown in Figure 4(a).In this scheme, the computer (PC Dell with Intel Core Quad 2.40 GHz), provided with Matlab/Simulink Real-Time and Quanser Quarq v 2.2, implements the controls used: continuous-time and discrete-time LQR, continuous-time predictor CSMC and discrete-time predictor DSMC.Analogto-digital converters are used in communication with the bench system.The equipment used in the experiments is presented in Figure 4(b), realizing that the equipment remains fixed on a bench, and the road surface is emulated by the movement of the lower plate of the equipment [20].
The format of the road surface was chosen based on the work of Sam and Osman [3] and Sam et al. [4], which represent a typical road surface by (51).It used  = 0.03 m for 0.5 ≤  ≤ 0.75 and 5 ≤  ≤ 5.25 and  = 0.01 m for 3 ≤  ≤ 3.25.In Figure 5 is presented the road surface profile used: However, other road surfaces possibilities exist in the literature and could have been used like square wave [26] and random road profile [6,7], among others [1,5].
To check the robustness of the CSMD and DSMD controls, a disturbance in the control signal actuator was introduced in most tests.The disturbance used was a sine wave with frequency 0.2 Hz and amplitude equal to 4 N.

Obtained Results and Discussion
4.1.Predictor Necessity.In Figures 6 and 7 are presented simulation results for the CSMC and DSMC controls with sampling periods of 1 ms and 3 ms, respectively, and with a delay of 12 sample periods (12 ms and 36 ms) without predictors.As shown in these figures, data acquisition delays let CSMC and DSMC controls be ineffective.
In the Sections 4.2 and 4.3 is shown in simulation and practical experiments that the use of state predictors is effective means that enable use of variable structure control when there is data acquisition delay.

Continuous-Time Control Results
. In Figures 8 and 9 is shown the response of the control system using CSMC with continuous predictor (CP-CSMC) compared with the LQR controller results, proposed by the manufacturer.In these graphs the system suffers delay of 12 sampling periods (12 ms) and disturbances in the control input.These graphs contain simulation and experimental results (using Quanser's bench active suspension system) on which one notes that the CP-CSMC rejected the disturbance.
In         shown.Even in this extreme condition, the control model composed of the pair continuous predictor-CSMC remained with good performance, while the LQR control showed very bad performance.The graphs in Figures 8-11 demonstrate the effectiveness of continuous predictor along with CSMC to control systems with data acquisition delay when are used small sampling periods.
When this sampling period increases, the continuous control displays deterioration in their performance.It is no different from the proposed control, as can be seen in Figures 12 and 13 where, even in the absence of disturbance in the control input, the control performance of continuous predictor-CSMC is unsatisfactory when the sampling period is 3 ms.
To deal with larger sampling periods, where the emulation becomes inefficient, one uses the discrete control approach, whose results are presented in Section 4.3.

Discrete-Time Control Results
. This test stage uses sampling periods of 3 ms, during which continuous control showed poor results.In Figures 14,15,16,and 17 are presented simulation and bench results to delays in data acquisition from 12 and 60 sampling periods (36 ms and 180 ms) and disturbance in the control input of the system with discrete predictor-DSMC (DP-DSMC).The results are compared with the discrete LQR results (DLQR) designed with matrices proposed by Quanser, the manufacturer of the active suspension system bench.
Note that the pair discrete predictor-DSMC remained with good results even when the system suffers large delays in data acquisition and the LQR controller, despite its lower sensitivity to such delays, cannot resist the delay of 60 sampling periods of delay, and bench results were not possible in this case.

Conclusions
In this paper, the problem of data acquisition delays in continuous-time and discrete-time sliding mode control when performed by a computer was analyzed.
An approach containing a state predictor has been successfully employed in order to overcome the problem of data acquisition delay, even in the presence of large delays, along with either continuous or discrete sliding mode control.
In the absence of state predictor, it was shown that controls CSMC and DSMC have poor results when there is data acquisition delay in the system.Moreover, it was shown that the use of predictors is able to overcome this difficulty, in simulations and practical applications.So the CP-CSMC and DP-DSMC controllers reject the matched uncertainty, even with delay time presence.

Figure 1 :
Figure 1: Representation of data acquisition delay using emulation.

Figure 3 :
Figure 3: Comparison scheme between control methods.

Figure 6 :
Figure 6: Simulation results for CSMC control with 12 sampling periods' delay without predictor.(a) Car body displacement and (b) control signal.

Figure 7 :
Figure 7: Simulation results for DSMC control with 12 sampling periods' delay without predictor.(a) Car body displacement and (b) control signal.

Figure 8 :
Figure 8: Car body displacement with delay of 12 sampling periods in data acquisition and disturbance in the control input.

Figure 9 :
Figure 9: Control signal with delay of 12 sampling periods in data acquisition and disturbance in the control input.

Figure 10 :
Figure 10: Car body displacement with delay of 60 sampling periods in data acquisition and disturbance in the control input.

Figure 11 :
Figure 11: Control signal with delay of 60 sampling periods in data acquisition and disturbance in the control input.

Figure 12 :
Figure 12: Comparison between car body displacements with delay of 12 sampling periods in data acquisition without disturbance in the control input using different sampling periods.

Figure 13 :Figure 14 :Figure 15 :Figure 16 :
Figure 13: Comparison between control signals with delay of 12 sampling periods in data acquisition without disturbance in the control input using different sampling periods.

Figure 17 :
Figure 17: Control signal with delay of 60 sampling periods in data acquisition and disturbance in the control input.