Next Article in Journal
Pedestrian Gender Recognition by Style Transfer of Visible-Light Image to Infrared-Light Image Based on an Attention-Guided Generative Adversarial Network
Next Article in Special Issue
An Algebraic Approach for Identification of Rotordynamic Parameters in Bearings with Linearized Force Coefficients
Previous Article in Journal
High Persistence and Nonlinear Behavior in Financial Variables: A More Powerful Unit Root Testing in the ESTAR Framework
Previous Article in Special Issue
Adaptive Robust Motion Control of Quadrotor Systems Using Artificial Neural Networks and Particle Swarm Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Suspension Control Using an MPC-LQR-LPV Controller with Attraction Sets and Quadratic Stability Conditions

by
Daniel Rodriguez-Guevara
1,
Antonio Favela-Contreras
1,*,
Francisco Beltran-Carbajal
2,
David Sotelo
1 and
Carlos Sotelo
1
1
Tecnologico de Monterrey, School of Engineering and Sciences, Ave. Eugenio Garza Sada 2501, Monterrey 64849, Mexico
2
Departamento de Energía, Universidad Autónoma Metropolitana, Unidad Azcapotzalco, Av. San Pablo No. 180, Col. Reynosa Tamaulipas, Mexico City 02200, Mexico
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(20), 2533; https://doi.org/10.3390/math9202533
Submission received: 26 July 2021 / Revised: 25 August 2021 / Accepted: 25 August 2021 / Published: 9 October 2021

Abstract

:
The control of an automotive suspension system by means of a hydraulic actuator is a complex nonlinear control problem. In this work, a linear parameter varying (LPV) model is proposed to reduce the complexity of the system while preserving the nonlinear behavior. In terms of control, a dual controller consisting of a model predictive control (MPC) and a Linear Quadratic Regulator (LQR) is implemented. To ensure stability, quadratic stability conditions are imposed in terms of Linear Matrix Inequalities (LMI). Simulation results for quarter-car model over several disturbances are tested in both frequency and time domain to show the effectiveness of the proposed algorithm.

1. Introduction

A vehicle can experiment different road disturbances while maneuvering in normal conditions such as bumps or bends. The suspension system of a car is designed to attenuate those disturbances to preserve comfort for the passengers while maintaining safe driving conditions to control the car’s direction. However, when road conditions are harsh, passive suspension systems may fail to preserve both comfort and road holding.
The active suspension system has been used to improve road-holding conditions while improving the comfort of passengers by means of a hydraulic actuator. To provide active suspension control, several control strategies have been proposed in the literature such as PID controller [1,2,3,4], H2 and H∞ control [5,6,7,8,9], fuzzy logic control [10,11,12,13], and sliding mode control [14,15,16,17]. All these controllers have exhibited a trade-off between comfort and road holding, with specific tuning conditions to manage each one of the design specifications according to the desired performance.
Another control strategy widely used in active suspensions is model predictive control (MPC). MPC approaches encompass several MIMO control strategies involving the prediction of the future behavior of the system along a prediction horizon N p and finding an optimal control solution subject to constraints in inputs, outputs and states. The general structure of the MPC strategy is to solve at each step an optimization problem where a cost function is minimized subject to constraints to find the optimal input sequence to be introduced in the system.
Some MPC approaches for active suspension systems are the following. In [18] an autoregressive with exogenous variable (ARX) model-based predictive control is presented to improve passenger comfort and road holding in a vehicle using a semi-active suspension with a Bouc–Wen representation. The results showed improvement when compared with passive suspension; however, the results were limited by the type of suspension used in this application. In [19] a full-car suspension model is controlled by a linear MPC with 6 degrees of freedom (6-DOF). In this approach, the control goal is to achieve a desired tilt angle to preserve comfort and road holding. The actuator is considered to be ideal, which results in a linear behavior of the system. This allows real-time implementation due to the short optimization time; however, the performance of the control algorithm may not be as effective as the one shown in simulation due to the linear design of the controller. Other MPC approaches using linear models are presented in [20,21,22].
Another MPC approach for the active suspension system is presented in [23]. In previous work, the active suspension quarter-car system was modeled as nonlinear by considering the nonlinear effect of the actuator. To comply with the MPC approach, the system is modeled as a Takagi–Sugeno model (T–S) by a fuzzy representation consisting of two sub-models. The model predictive control is also designed as a fuzzy MPC where there exist two interconnected linear models and the switching between one model and the another is performed by fuzzy logic. Additionally, terminal equality constraints are included in terms of Linear Matrix Inequalities (LMI) to ensure stability. In [24] a robust model predictive controller (RMPC) for an active suspension full-car system is presented. This approach considers both the nonlinearities of the hydraulic actuator and the nonlinearities presented by the relationship of the movement of each one of the four corners of the car. Therefore, to design the MPC, the model is simplified into a linear fuzzy logic system. To add robustness to the control approach, an adaptive control law is proposed based on the MPC and a fractional PID controller.
As shown in the previous works, the nonlinear model of an active suspension system is represented by linear representations to comply with the MPC strategy. In this research work, a linear parameter varying (LPV) representation is proposed. This kind of representation is common in semi-active suspension control approaches [25,26,27]. In this approach, the LPV representation is done using one scheduling parameter embedding the nonlinearities of the hydraulic actuator.
Therefore, the proposed control strategy consists of a model predictive controller for an active suspension system with an electro-hydraulic actuator with a servo spool valve. The model of the system is constructed as a linear parameter varying model using one scheduling variable ρ 1 . quadratic stability conditions are included in the MPC algorithm as LMI, as presented in [28]. To improve performance, a terminal cost using attraction sets is included, as shown in [29]. Finally, the inclusion of a terminal set and a Linear Quadratic Regulator (LQR) controller in the terminal set is included.
The rest of the paper is organized in the following structure. Section 2 presents the active suspension with electro-hydraulic actuator model. Section 3 shows a state–space LPV representation of the active suspension. Section 4 describes the MPC-LPV control algorithm. Section 5 shows a Recursive Least Squares (RLS) algorithm for the prediction of the scheduling parameter along the prediction horizon. Section 6 presents the quadratic stability conditions for the MPC-LPV approach. Section 7 describes the attraction sets and terminal set for control switching. Results and simulations are presented in Section 8 and conclusions are discussed in Section 9.

2. Quarter-Car Active Suspension Model

active suspension systems add an actuator to the passive system mainly consisting of the wheel mass and the chassis mass. Figure 1 presents a schematic model of an active suspension system as found in [7]. In this model, the actuator produces a force f s which reduces the vertical movement of both masses m s , which is the sprung mass representing the chassis body, and m u s , which is the unsprung mass representing the suspension unit and wheel of the quarter-car.
The dynamic force equations of the system are the following:
m s z ¨ s + k s ( z s z u s ) + b s ( z ˙ s z ˙ u s ) f s = 0
m u s z ¨ u s k s ( z s z u s ) b s ( z ˙ s z ˙ u s ) + k t z u s r ( t ) + f s = 0
with k s being the constant of the spring between the two masses and b s being the damping coefficient. k t represents the tire elastic constant and r ( t ) represents the road disturbances. The force f s is generated by an electro-hydraulic actuator with a servo spool valve. A schematic of the electro-hydraulic actuator is shown in Figure 2. By means of this actuator, the force f s is generated by the pressure supplied to the system P l and the area of the moving piston A. Therefore, the force can be expressed using the following equation.
f s = A P l
The pressure P l dynamics are modeled using the following differential Equation (4).
V t 4 β e P l ˙ = Q C t p P l A ( z s ˙ z u s ˙ )
with Q = sgn P s sgn ( z v ) P l C d w x v 1 ρ | P s sgn ( z v ) P l | , where V t is the total actuator volume, Q the load flow, B e the effective bulk modulus, C t p the piston leakage coefficient, C d the discharge coefficient, w the spool valve area gradient, ρ the hydraulic fluid density and P s the pressure supply. The displacement of the spool valve z v is proportional to the control action u ( k ) which is a voltage signal. The valve displacement is expressed by the following equation.
z v ˙ = z v τ + k v u ( k )
where k v represents the valve gain and τ is a scaling factor.

3. LPV-SS Representation of the Quarter-Car Active Suspension Model

To comply with the MPC strategy, Equations (1) through (5) will be used to build a linear parameter varying state–space (LPV-SS) model with a scheduling parameter ρ 1 of the form:
x ˙ ( t ) = A ρ 1 ( t ) x ( t ) + B u ( t ) + B r r ( t )
with the discrete LPV-SS obtained by a Zero-Order Hold (ZOH) represented as:
x ( k + 1 ) = A ρ 1 ( k ) x ( k ) + B u ( k ) + B r r ( k )
where x ˙ ( t ) = [ x 1 x 2 x 3 x 4 x 5 x 6 ] T = [ z s z s ˙ z u s z u s ˙ P l z v ] T . With A being the state matrix and B the input matrix, B r represents the input disturbance matrix while u ( k ) is the control input. Therefore, the space-state matrices can be defined as the following:
A ρ 1 ( k ) = 0 1 0 0 0 0 k s m s b s m s k s m s b s m s A 0 0 0 0 1 0 0 k s m u s b s m u s k s + k t m u s b s m u s A m u s 0 0 α A 0 α A β γ ρ 1 0 0 0 0 0 1 τ
B = 0 0 0 0 0 k v τ
B r = 0 0 0 k t m u s 0 0
With α = 4 β e V t , β = α C t p , γ = α C d 1 ρ and ρ 1 = sgn P s sgn ( x 6 ) x 5 | P s sgn ( x 6 ) x 5 | . The inclusion of the scheduling variable ρ 1 allows the system to be expressed as an LPV-SS representation which allows the MPC law to be computed in a compact matrix form.

4. LPV-MPC Controller

To develop a MPC scheme for LPV models, the future states need to be formulated so a trajectory can be formed along the prediction horizon. The i-steps-ahead prediction can be structured as the following:
x ( k + i | k ) = j = 0 i 1 A ρ 1 ( k + j ) x ( k ) + s = 1 i 1 l = s i 1 A ρ 1 ( k + l ) B u ( k + s 1 ) + B u ( k + i 1 )
Prediction of the future states needs to be performed for the future N p time steps, thus using (11) the following matrix equation can be deduced.
X = Φ x ( k ) + Ψ U
where:
X = x ( k + 1 | k ) x ( k + 2 | k ) : . x ( k + N p | k )
Φ = A ( ρ 1 ( k ) ) j = 0 1 A ρ 1 ( k + j ) j = 0 N p 1 A ρ 1 ( k + j )
Ψ = B 0 n x × n u 0 n x × n u A ρ 1 ( k + 1 ) B B 0 n x × n u A ρ 1 ( k + 2 ) A ρ 1 ( k + 1 ) B A ρ 1 ( k + 2 ) B 0 n x × n u i = 1 N p 1 A ρ 1 ( k + i ) B i = 1 N p 1 A ρ 1 ( k + i + 1 ) B B
U = u ( k ) u ( k + 1 ) : . u ( k + N p 1 )
with X R N p · n x , Φ R N p · n x × n x , Ψ R N p · n x × N p · n u and U R N p · n u where n x is the number of states and n u the number of inputs. With the state prediction equation, we can construct a cost function to minimize the deviation from the equilibrium states and the energy used by the inputs, so that the cost function is defined as:
J = X T Q c X + U T R c U
where Q c and R c are weight matrix of appropriate dimensions. To find the optimal control trajectory U , (17) needs to be minimized subject to the constraints in the inputs (18) and the states (19).
u m i n U u m a x
x m i n X x m a x
with both u m i n & u m a x R N p × n u and both x m i n & x m a x R N p × N p . However, to properly solve the MPC problem, the future values of ρ 1 , which are unknown, must be estimated. To obtain the values we must estimate them using an RLS approach to obtain an approximate value of the scheduling parameter based on its previous behavior and the system response.

5. Scheduling Parameter Prediction Using RLS

To obtain an estimation of the future scheduling parameter, an RLS approach is used as presented by Sename, Morato & Normey-Rico in [26]. The scheduling parameter is assumed to be measurable at instant k and all previous values can be stored; however, the future parameters will be estimated based on the previous measurements of the scheduling parameters as well as the previous inputs and outputs to consider the behavior of the system.
The behavior of the scheduling parameter of a LPV system can be approximated by a linear ARX model, which is a function of the previous scheduling parameter values, the previous inputs and the previous outputs. This ARX model can be represented as:
ρ 1 ( k + N p ) = a 0 ρ 1 ( k ) + + a N p ρ 1 ( k N p ) + b 0 u ( k 1 ) +
+ b N p u ( k N p 1 ) + c 0 y ( k ) + + c N p y ( k N p )
Afterwards, (20) can be expressed in a compact form and be dependent only on known values to be suitable for MPC design. To find a solution to the RLS, parameters a 0 to c N p need to be calculated. These parameters will be grouped into the following vector:
Θ ( k ) = [ a 0 c N p ] T
resulting in:
ρ 1 ( k ) = γ ( k ) T Θ ( k )
with:
γ ( k ) T = ρ 1 ( k N p ) , . . . , ρ 1 ( k 2 N p ) , u ( k N p 1 ) , u ( k 2 N p 1 ) , y ( k N p ) , . . . , y ( k 2 N p )
with (22) and (23) a direct solution can be built and used to find ρ 1 in an online RLS algorithm as presented in [30]:
Θ ( k ) = Θ ( k 1 ) + σ ( k ) ρ 1 ( k ) γ ( k 1 ) T Θ ( k 1 )
Q ^ ( k ) = I σ ( k ) γ ( k ) T Q ^ ( k 1 ) μ
with μ [ 0 , 1 ] being a forgetting factor that gives exponentially less weight to older error samples of the RLS algorithm and σ ( k ) being a vector defined as:
σ ( k ) = 1 μ c ( k ) Q ^ ( k 1 ) γ ( k )
and c ( k ) is a scalar defined by:
c ( k ) = 1 + γ ( k ) T Q ^ ( k 1 ) μ γ ( k )
Therefore, the RLS algorithm for estimating the future scheduling parameters is shown as Algorithm 1:
After solving the RLS algorithm for the N p future scheduling parameters, they will be considered to be known and exact to build a vector P ^ ( k ) = ρ 1 ( k ) , , ρ 1 ( k + N p ) T which contains all of them; therefore, (12) is no longer an equation with unknown variables and can be solved through LMI optimization.
Algorithm 1
Offline
Step 1—Initialize Θ ( 0 ) and Q ^ ( 0 )
Online
Step 2—Obtain ρ 1 ( k ) , y ( k ) and u ( k )
Step 3—Construct γ T ( k ) vector
Step 4—Calculate scalar c
Step 5—Obtain vector σ ( k )
Step 6—Obtain Θ ( k )
Step 7—Obtain Q ^ ( k )
Step 8—Obtain ρ 1 ( k )
Step 9—Set k = k + 1 , If k < N p go to Step 10, else, go back to step 3
Step 10—Construct P ^ ( k ) = ρ 1 ( k ) , ρ 1 ( k + 1 ) , , ρ 1 ( k + N p )

6. Quadratic Stability in MPC-LPV Approach

To ensure quadratic stability in the MPC-LPV approach, system (7) can be considered to be a parametric uncertain system. In parametric uncertain systems, the scheduling variable is limited to vary in a range Δ ρ 1 m i n Δ ρ 1 k Δ ρ 1 m a x . To ensure stability in parametric uncertain systems, the following condition needs to be met as presented in [31].
A ( ρ 1 ) + B K T P A ( ρ 1 ) + B K P < 0
which is the Riccati Equation for parametric uncertain systems where P > 0 is a positive definite matrix of appropriate dimensions and K a static feedback gain matrix. Then, (28) can be pre- and post-multiplied by a matrix Q = P 1 and K Q = R to obtain:
Q A T ( ρ 1 ) + R T B T Q 1 A ( ρ 1 ) Q + B R Q < 0
To cope with the MPC paradigm, the Schur complement is applied to (29) to obtain the following LMI:
Q Q A T ( ρ 1 ) + R B A ( ρ 1 ) Q + B R Q > 0
for every possible value of ρ 1 at time instant k which leads to an infinite number of LMI. However, as system (7) is considered to be a parametric uncertain system, (30) can be evaluated on the vertex of matrix A to consider the worst-case scenarios. Therefore, (30) can be written as:
Q Q A i , j T ( ρ 1 ) + R B A i , j ( ρ 1 ) Q + B R Q > 0
The previous condition must be met j [ k , k + N p ] and i [ 1 , 2 l ] , where l is the number of scheduling variables ρ 1 , Q > 0 is a positive definite stability matrix to be determined, and K Q = R , where K is the static feedback gain matrix. With these adjustments, the number of LMI to be solved is now finite and equal to 2 l N p . Since there is the consideration of a static feedback gain, the control law is determined as u ( k ) = K x ( k ) , but to comply with the MPC paradigm, the previous expression can be considered to be an inequality as u ( k ) < K x ( k ) . This leads to a conservative MPC performance due to the limitations of the input variable. However, this problem will be addressed in Section 7 with the inclusion of terminal sets. Therefore, using (17)–(19) and (31) the optimization problem needs to find the optimal control sequence at each time step k is the following:
min U J s . t . ( 18 ) , ( 19 ) & ( 31 )

7. MPC-LQR for LPV Models

7.1. Attraction Sets and Terminal Set

The inclusion of LMIs to ensure robust stability to the MPC paradigm often leads to a conservative performance of the control of the system. Therefore, to steer the system into a desired equilibrium state in the presence of disturbance or uncertainty, a series of terminal sets can be defined. In [32] a set of shrinking ellipsoids is determined using a decay rate, which can vary the speed of the system and the stability determined by similar stability conditions to the ones shown in Section 6, to steer the states to the equilibrium point. However, the determination of the decay rate and the constructions of the ellipsoids make this algorithm too slow for real-time applications and is rather implemented as an offline algorithm. In [33] a set of ellipsoidal sets are defined to predict the behavior of the system in the presence of bounded disturbances and uncertain bounded parameter changes. In [34] a path of ellipsoids is defined to predict the possible behavior of the scheduling parameter along the prediction horizon. In all three approaches, the goal of the ellipsoidal sets is that the states reach a terminal set or a terminal point, where a state-dependent stationary gain is applied to the system instead of the MPC law.
In this work, the future scheduling parameter is not known but predicted using the RLS algorithm presented in Section 3; therefore, the ellipsoids to build do not consider a variation on the scheduling parameter but rather the prediction error generated by the RLS algorithm. To generate the optimal desired trajectory to the setpoint, a path must be defined from every possible initial state to the terminal ellipsoidal set.
To steer the system into the desired terminal set, a term J T S is added to the cost function J presented in (17). J T S is defined as the following:
J T S = x ( k + N p ) ( x d s + x d i s t ) T L x ( k + N p ) ( x d s + x d i s t ) E ( ρ 1 )
where E ( ρ 1 ) = i = k k + N p ρ 1 ( i ) γ ( i 1 ) T Θ ( i 1 ) 2 represents the sum of the squared errors of the prediction of the future parameter values. x ( k + N p ) are the predicted states at the end of the prediction horizon, x d s represents the desired state after N p steps, x d i s t is the predicted effect of the disturbance on the states N p steps ahead and it was obtained by performing an open loop simulation of every possible disturbance from every initial set of states. Both x d s and x d i s t were computed offline and stored in a lookup table. L is a weighing matrix of appropriate dimensions. Therefore, (17) is redefined as:
J = X T Q c X + U T R c U + J T S
However, the computation of every desired trajectory for every state needs to be computed offline and stored in a lookup table before the implementation of the MPC algorithm to increase execution speed.

7.2. MPC-LQR Dual Controller

To reduce the computational load of the algorithm, when the current states reach a terminal invariant set around the equilibrium point, the MPC algorithm does not need to be computed. Instead, an LQR gain can be computed based on the value of the actual prediction parameter to cope with the small error that may be present inside the terminal invariant set. The control law is then presented as:
u ( k ) = U mpc x ( k ) T K LQR ( ρ 1 ) x ( k ) x ( k ) T
Where K LQR ( ρ 1 ) is the LQR gain dependent on the scheduling parameter ρ 1 and T is the terminal invariant set defined around the equilibrium point of the system.
Figure 3 presents the block diagram for the proposed LPV-MPC-LQR control strategy. Additionally, the LPV-MPC-LQR algorithm is shown in the flowchart presented in Figure 4.

8. Results and Discussion

The following simulations are made to observe the advantages and performance of implementing the proposed LQR-MPC-LPV algorithm described in Section 7. The algorithm was tested in the active suspension system described in Section 2. Table 1 shows the specifications of the active suspension system obtained from [7].
A discretization is made to comply with the MPC paradigm using a sampling time of T s = 10 ms. A prediction horizon of N p = 3 was defined after several tests using different prediction horizons were made. Using a larger prediction horizon resulted in longer optimization time and more inexact variable scheduling predictions while it does not exhibit a significant improvement in control performance. The control objective is to steer all the states to the origin while complying with the following constraints.
12 V u ( k ) 12 V
1 cm z v 1 cm
The results will be divided into frequency-domain results and time-domain results.

8.1. Frequency-Domain Results

To obtain a frequency analysis in the nonlinear active suspension system, an algorithm similar to the process of defining the system response as a describing function is used. To produce these results, Algorithm 1 presented in [35] is issued. Figure 5 presents the frequency response of the active suspension deflection gain using the proposed LPV-MPC-LQR algorithm. Figure 6 presents the frequency response of the acceleration of the chassis mass using the proposed LPV-MPC-LQR algorithm; also, the frequency response plots are compared with the ones presented in [36] which use an LPV gain scheduling approach.
The results show how the suspension deflection is attenuated at every frequency, which results in better road holding and driving conditions. Additionally, the chassis acceleration stays in values which guarantee passenger comfort. Compared to the frequency responses of the work of Fialho et al. [36] the MPC-LQR-LPV approach presents an improvement especially in terms of road holding, shown in Figure 5, without affecting the passenger comfort.

8.2. Time-Domain Results

To obtain time-domain results using the proposed LPV-MPC-LQR control algorithm, the system was simulated using two different disturbances. Figure 7, Figure 8 and Figure 9 present the suspension behavior when a bump disturbance of 5 cm is introduced. Figure 10, Figure 11 and Figure 12 show the suspension behavior when driving through a sinusoidal road. The system was simulated using Matlab®; also, the software package YALMIP [37] using QP-solver SDPT3 was used for the MPC optimization. The results presented by [7] are included to make a comparison. Additionally, the results using the MPC with a frozen scheduling parameter approach without using the RLS to show the effect of the scheduling variable prediction in control performance are included.
The results of both displacement and deflection show a better performance, which results in better road holding while maintaining passenger comfort. Additionally, the comfort exhibits improvement in terms of chassis acceleration as shown in Figure 9. Additionally, to express the results numerically, both the RMS value and the maximum value of the displacement of the chassis, the suspension deflection and the acceleration of the chassis are presented in Table 2 and Table 3 respectively.
Similar to the bump disturbance case, the proposed LPV-MPC-LQR control strategy exhibits better performance in both displacement and deflection, which results in better road holding. In terms of comfort, the acceleration of the chassis presented in Figure 12 shows a major improvement. Table 4 presents the peak values for the displacement of the chassis, the suspension deflection, and the acceleration of the chassis.
As shown in the previous figures, the proposed LPV-MPC-LQR control algorithm presents a better performance when compared with the H2 control strategy in both disturbance cases (bump disturbance and sinusoidal road disturbance). The RLS prediction of the future scheduling parameters have improved the control performance as well. Additionally, the proposed algorithm shows an appropriate optimization time with a worst optimization time of 930 ms and an average optimization time of 93 ms.

9. Conclusions

In this paper, a novel LPV-MPC-LQR control algorithm ensuring quadratic stability and with the inclusion of attraction sets was presented. This method runs an RLS algorithm to obtain the prediction of the future scheduling parameter values, which simplifies the prediction of the future states while ensuring quadratic stability. This application can cope with nonlinear systems that can be embedded into LPV representation and therefore reduce the complexity of the algorithm and allow fast execution times. This control strategy was designed and tested on a nonlinear active suspension system. The results show improvements to the performance of the active suspension in terms of road holding and passenger comfort. Future research works should deal with recursive feasibility analysis based on stability conditions, and robustness analysis. Optimization of the LPV-MPC-LQR algorithm to achieve faster execution times using techniques of the embedded systems will also be considered in future works as well.

Author Contributions

All Authors D.R.-G., A.F.-C., F.B.-C., D.S. and C.S. have contributed as follows: Conceptualization, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; Methodology, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; Software, D.R.-G., D.S. and C.S.; Validation, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; Formal analysis, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; Investigation, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; Writing—original draft preparation, D.R.-G. and A.F.-C.; Writing—review and editing, D.R.-G., A.F.-C., F.B.-C., D.S. and C.S.; supervision, A.F.-C. and F.B.-C.; project administration A.F.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank Consejo Nacional de Ciencia y Tecnología (CONACyT) and Tecnológico de Monterrey for the financial support to conduct the present research. Additionally, thanks go to the Sensors and Devices Research Group and the Robotics Research Group from the School of Engineering and Sciences of Tecnológico de Monterrey for the support given to develop this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mouleeswaran, S. Design and development of PID controller-based active suspension system for automobiles. In PID Controller Design Approaches-Theory, Tuning and Application to Frontier Areas; Vagia, M., Ed.; Intech: Rijeka, Croatia, 2012; pp. 71–98. [Google Scholar]
  2. Ahmed, A.E.N.S.; Ali, A.S.; Ghazaly, N.M.; Abd el-Jaber, G.T. PID controller of active suspension system for a quarter car model. Int. J. Adv. Eng. Technol. 2015, 8, 899. [Google Scholar]
  3. Bello, M.M.; Shafie, A.A.; Khan, M.R. Electro-hydraulic pid force control for nonlinear vehicle suspension system. Int. J. Eng. Res. Technol. 2015, 4, 517–524. [Google Scholar]
  4. Phu, D.X.; An, J.-H.; Choi, S.-B. A Novel Adaptive PID Controller with Application to Vibration Control of a Semi-Active Vehicle Seat Suspension. Appl. Sci. 2017, 7, 1055. [Google Scholar] [CrossRef]
  5. Yu, S.; Wang, F.; Wang, J.; Chen, H. Full-car active suspension based on H2/generalised H2 output feedback control. Int. J. Veh. Des. 2015, 68, 37–54. [Google Scholar] [CrossRef]
  6. Afshar, K.K.; Javadi, A.; Jahed-Motlagh, M.R. Robust H∞ control of an active suspension system with actuator time delay by predictor feedback. IET Control. Theory Appl. 2018, 12, 1012–1023. [Google Scholar] [CrossRef]
  7. Ghazaly, N.M.; Ahmed, A.E.N.S.; Ali, A.S.; Abd, E.-J.G.T. H.∞ Control of active suspension System for a Quarter Car Model. Int. J. Veh. Struct. Syst. (IJVSS) 2016, 8, 35–40. [Google Scholar] [CrossRef]
  8. Kaleemullah, M.; Faris, W.F.; Hasbullah, F. Control of active suspension system using robust H∞ control with genetic algorithm. Int. J. Adv. Sci. Technol. 2019, 28, 763–782. [Google Scholar]
  9. Jin, X.; Wang, J.; Sun, S.; Li, S.; Yang, J.; Yan, Z. Design of Constrained Robust Controller for active suspension of In-Wheel-Drive Electric Vehicles. Mathematics 2021, 9, 249. [Google Scholar] [CrossRef]
  10. Wen, S.; Chen, M.Z.; Zeng, Z.; Yu, X.; Huang, T. Fuzzy control for uncertain vehicle active suspension systems via dynamic sliding-mode approach. IEEE Trans. Syst. Man Cybern. Syst. 2016, 47, 24–32. [Google Scholar] [CrossRef]
  11. Mustafa, G.I.; Wang, H.; Tian, Y. Model-free adaptive fuzzy logic control for a half-car active suspension system. Stud. Inform. Control 2019, 28, 13–24. [Google Scholar] [CrossRef]
  12. Palanisamy, S.; Karuppan, S. Fuzzy control of active suspension system. J. Vibroeng. 2016, 18, 3197–3204. [Google Scholar] [CrossRef] [Green Version]
  13. Alfadhli, A.; Darling, J.; Hillis, A.J. The Control of an Active Seat Suspension Using an Optimised Fuzzy Logic Controller, Based on Preview Information from a Full Vehicle Model. Vibration 2018, 1, 20–40. [Google Scholar] [CrossRef] [Green Version]
  14. Rath, J.J.; Defoort, M.; Karimi, H.R.; Veluvolu, K.C. Output feedback active suspension control with higher order terminal sliding mode. IEEE Trans. Ind. Electron. 2016, 64, 1392–1403. [Google Scholar] [CrossRef]
  15. Chokor, A.; Talj, R.; Charara, A.; Shraim, H.; Francis, C. active suspension control to improve passengers comfort and vehicle’s stability. In Proceedings of the IEEE 19th International Conference on Intelligent Transportation Systems (ITSC), Rio de Janeiro, Brazil, 1–4 November 2016; pp. 296–301. [Google Scholar]
  16. Taghavifar, H.; Mardani, A.; Hu, C.; Qin, Y. Adaptive robust nonlinear active suspension control using an observer-based modified sliding mode interval type-2 fuzzy neural network. IEEE Trans. Intell. Veh. 2019, 5, 53–62. [Google Scholar] [CrossRef]
  17. Aljarbouh, A.; Fayaz, M. Hybrid Modelling and Sliding Mode Control of Semi-active suspension Systems for Both Ride Comfort and Road-Holding. Symmetry 2020, 12, 1286. [Google Scholar] [CrossRef]
  18. Piñón, A.; Favela-Contreras, A.; Félix-Herrán, L.C.; Beltran-Carbajal, F.; Lozoya, C. An ARX Model-Based Predictive Control of a Semi-Active Vehicle Suspension to Improve Passenger Comfort and Road-Holding. Actuators 2021, 10, 47. [Google Scholar] [CrossRef]
  19. Yao, J.; Wang, M.; Li, Z.; Jia, Y. Research on model predictive control for automobile active tilt based on active suspension. Energies 2021, 14, 671. [Google Scholar] [CrossRef]
  20. Theunissen, J.; Sorniotti, A.; Gruber, P.; Fallah, S.; Ricco, M.; Kvasnica, M.; Dhaens, M. Regionless explicit model predictive control of active suspension systems with preview. IEEE Trans. Ind. Electron. 2019, 67, 4877–4888. [Google Scholar] [CrossRef]
  21. Enders, E.; Burkhard, G.; Munzinger, N. Analysis of the Influence of Suspension Actuator Limitations on Ride Comfort in Passenger Cars Using model predictive control. Actuators 2020, 9, 77. [Google Scholar] [CrossRef]
  22. Narayan, J.; Gorji, S.A.; Ektesabi, M.M. Power reduction for an active suspension system in a quarter car model using MPC. In Proceedings of the IEEE International Conference on Energy Internet (ICEI), Sydney, Australia, 24–28 August 2020; pp. 140–146. [Google Scholar]
  23. Bououden, S.; Chadli, M.; Karimi, H.R. A robust predictive control design for nonlinear active suspension systems. Asian J. Control 2016, 18, 122–132. [Google Scholar] [CrossRef]
  24. Wang, D.; Zhao, D.; Gong, M.; Yang, B. Research on robust model predictive control for electro-hydraulic servo active suspension systems. IEEE Access 2017, 6, 3231–3240. [Google Scholar] [CrossRef]
  25. Morato, M.M.; Sename, O.; Dugard, L. LPV-MPC fault tolerant control of automotive suspension dampers. IFAC-PapersOnLine 2018, 51, 31–36. [Google Scholar] [CrossRef]
  26. Morato, M.M.; Normey-Rico, J.E.; Sename, O. Novel qLPV MPC design with least-squares scheduling prediction. IFAC-PapersOnLine 2019, 52, 158–163. [Google Scholar] [CrossRef]
  27. Morato, M.M.; Normey-Rico, J.E.; Sename, O. Sub-optimal recursively feasible Linear Parameter-Varying predictive algorithm for semi-active suspension control. IET Control Theory Appl. 2020, 14, 2764–2775. [Google Scholar] [CrossRef]
  28. Boyd, S.; Balakrishnan, V.; Feron, E.; ElGhaoui, L. Control system analysis and synthesis via linear matrix inequalities. In Proceedings of the American Control Conference, San Francisco, CA, USA, 2–4 June 1993; pp. 2147–2154. [Google Scholar]
  29. Ferramosca, A.; Limón, D.; González, A.H.; Odloak, D.; Camacho, E.F. MPC for tracking target sets. In Proceedings of the 48h IEEE Conference on Decision and Control (CDC) Held Jointly with 2009 28th Chinese Control Conference, Shanghai, China, 15–18 December 2009; pp. 8020–8025. [Google Scholar]
  30. Franklin, G.F.; Powell, J.D.; Workman, M.L. Digital Control of Dynamic Systems; Addison-Wesley: Reading, MA, USA, 1988; Volume 3. [Google Scholar]
  31. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1994; Chapter 5; pp. 61–76. [Google Scholar]
  32. Longge, Z.; Yan, Y. Robust shrinking ellipsoid model predictive control for linear parameter varying system. PLoS ONE 2017, 12, e0178625. [Google Scholar] [CrossRef] [Green Version]
  33. Ping, X.B.; Wang, P.; Zhang, J.F. A Multi-step Output Feedback Robust MPC Approach for LPV Systems with Bounded Parameter Changes and Disturbance. Int. J. Control Autom. Syst. 2018, 16, 2157–2168. [Google Scholar] [CrossRef]
  34. Suzukia, H.; Sugie, T. MPC for LPV systems with bounded parameter variation using ellipsoidal set prediction. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; p. 6. [Google Scholar]
  35. Poussot-Vassal, C.; Savaresi, S.M.; Spelta, C.; Sename, O.; Dugard, L. A methodology for optimal semi-active suspension systems performance evaluation. In Proceedings of the 49th IEEE Conference on Decision and Control (CDC), Atlanta, GA, USA, 15–17 December 2017; pp. 2892–2897. [Google Scholar]
  36. Fialho, I.; Balas, G.J. Road adaptive active suspension design using linear parameter-varying gain-scheduling. IEEE Trans. Control. Syst. Technol. 2002, 10, 43–54. [Google Scholar] [CrossRef] [Green Version]
  37. Lofberg, J. Automatic robust convex programming. Optim. Methods Softw. 2012, 27, 115–129. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Active suspension System.
Figure 1. Active suspension System.
Mathematics 09 02533 g001
Figure 2. Electro-Hydraulic Actuator.
Figure 2. Electro-Hydraulic Actuator.
Mathematics 09 02533 g002
Figure 3. Block diagram of the proposed LPV-MPC-LQR control strategy for the active suspension system.
Figure 3. Block diagram of the proposed LPV-MPC-LQR control strategy for the active suspension system.
Mathematics 09 02533 g003
Figure 4. Flow diagram of the LPV-MPC-LQR control strategy.
Figure 4. Flow diagram of the LPV-MPC-LQR control strategy.
Mathematics 09 02533 g004
Figure 5. Frequency response of the active suspension deflection gain.
Figure 5. Frequency response of the active suspension deflection gain.
Mathematics 09 02533 g005
Figure 6. Frequency response of the chassis acceleration.
Figure 6. Frequency response of the chassis acceleration.
Mathematics 09 02533 g006
Figure 7. Chassis displacement—Bump Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Figure 7. Chassis displacement—Bump Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Mathematics 09 02533 g007
Figure 8. Suspension Deflection—Bump Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Figure 8. Suspension Deflection—Bump Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Mathematics 09 02533 g008
Figure 9. Chassis Acceleration—Bump disturbance (Blue—MPC-LQR, Red—Passive, Yellow—H2, Purple—MPC-Frozen).
Figure 9. Chassis Acceleration—Bump disturbance (Blue—MPC-LQR, Red—Passive, Yellow—H2, Purple—MPC-Frozen).
Mathematics 09 02533 g009
Figure 10. Chassis displacement—Sinusoidal Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Figure 10. Chassis displacement—Sinusoidal Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Mathematics 09 02533 g010
Figure 11. Suspension Deflection—Sinusoidal Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Figure 11. Suspension Deflection—Sinusoidal Disturbance (Blue—disturbance, Red—MPC-LQR, Yellow—Passive, Purple—H2, Green—MPC-Frozen).
Mathematics 09 02533 g011
Figure 12. Chassis Acceleration—Sinusoidal disturbance (Blue—MPC-LQR, Red—Passive, Yellow—H2, Purple—MPC-Frozen).
Figure 12. Chassis Acceleration—Sinusoidal disturbance (Blue—MPC-LQR, Red—Passive, Yellow—H2, Purple—MPC-Frozen).
Mathematics 09 02533 g012
Table 1. Constant Values of the active suspension system.
Table 1. Constant Values of the active suspension system.
VariableValueUnits
m s 250kg
m u s 50kg
m u s 50kg
k t 190,000N/m
k s 16,812N/m
b s 1000N/(m/s)
P s 10,300,000Pa
τ 1/30s
A 3.35 × 10 4 m 2
β 1s 1
α 4.515 × 10 13 N/m 5
k v 1 × 10 4 m/V
Table 2. RMS Values performance.
Table 2. RMS Values performance.
VariableMPC-LQR-LPVH2 (Ghazaly, 2016)PassiveMPC-Frozen
Chassis Displacement (m)0.00790.00910.01420.0107
Suspension Deflection (m)0.00890.01490.02400.0122
Chassis Acceleration (m/s 2 )0.07130.11040.10410.0838
Table 3. Max Values performance.
Table 3. Max Values performance.
VariableMPC-LQR-LPVH2 (Ghazaly, 2016)PassiveMPC-Frozen
Chassis Displacement (m)0.02930.02840.03800.0367
Suspension Deflection (m)0.03550.04990.04640.0439
Chassis Acceleration (m/s 2 )0.26440.39780.29250.2899
Table 4. Peak Values performance.
Table 4. Peak Values performance.
VariableMPC-LQR-LPVH2 (Ghazaly, 2016)PassiveMPC-Frozen
Chassis Displacement (m)0.00270.00440.01510.0055
Suspension Deflection (m)0.00340.00400.00970.0051
Chassis Acceleration (m/s 2 )0.01640.320.1620.0312
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodriguez-Guevara, D.; Favela-Contreras, A.; Beltran-Carbajal, F.; Sotelo, D.; Sotelo, C. Active Suspension Control Using an MPC-LQR-LPV Controller with Attraction Sets and Quadratic Stability Conditions. Mathematics 2021, 9, 2533. https://doi.org/10.3390/math9202533

AMA Style

Rodriguez-Guevara D, Favela-Contreras A, Beltran-Carbajal F, Sotelo D, Sotelo C. Active Suspension Control Using an MPC-LQR-LPV Controller with Attraction Sets and Quadratic Stability Conditions. Mathematics. 2021; 9(20):2533. https://doi.org/10.3390/math9202533

Chicago/Turabian Style

Rodriguez-Guevara, Daniel, Antonio Favela-Contreras, Francisco Beltran-Carbajal, David Sotelo, and Carlos Sotelo. 2021. "Active Suspension Control Using an MPC-LQR-LPV Controller with Attraction Sets and Quadratic Stability Conditions" Mathematics 9, no. 20: 2533. https://doi.org/10.3390/math9202533

APA Style

Rodriguez-Guevara, D., Favela-Contreras, A., Beltran-Carbajal, F., Sotelo, D., & Sotelo, C. (2021). Active Suspension Control Using an MPC-LQR-LPV Controller with Attraction Sets and Quadratic Stability Conditions. Mathematics, 9(20), 2533. https://doi.org/10.3390/math9202533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop