Next Article in Journal
Economic Analysis of the Application of the Technological System for Removing Suspended Solids from Mine Drainage Waters
Next Article in Special Issue
Innovative Energy-Saving Propulsion System for Low-Speed Biomimetic Underwater Vehicles
Previous Article in Journal
Energy Management Model for HVAC Control Supported by Reinforcement Learning
Previous Article in Special Issue
Automatic Clutch Engagement Control for Parallel Hybrid Electric Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Static Output Feedback and Structured Controllers for Active Suspension with Quarter-Car Model

1
Department of Electronic Engineering, Korea National University of Transportation, Chungju-si 27469, Korea
2
Research Center for Electrical and Information Technology, Seoul National University of Science and Technology, Seoul 01811, Korea
*
Author to whom correspondence should be addressed.
Energies 2021, 14(24), 8231; https://doi.org/10.3390/en14248231
Submission received: 16 November 2021 / Revised: 6 December 2021 / Accepted: 6 December 2021 / Published: 7 December 2021
(This article belongs to the Special Issue Vehicle Dynamics and Control)

Abstract

:
This paper presents a method to design active suspension controllers for a 7-Degree-of-Freedom (DOF) full-car (FC) model from controllers designed with a 2-DOF quarter-car (QC) one. A linear quadratic regulator (LQR) with 7-DOF FC model has been widely used for active suspension control. However, it is too hard to implement the LQR in real vehicles because it requires so many state variables to be precisely measured and has so many elements to be implemented in the gain matrix of the LQR. To cope with the problem, a 2-DOF QC model describing vertical motions of sprung and unsprung masses is adopted for controller design. LQR designed with the QC model has a simpler structure and much smaller number of gain elements than that designed with the FC one. In this paper, several controllers for the FC model are derived from LQR designed with the QC model. These controllers can give equivalent or better performance than that designed with the FC model in terms of ride comfort. In order to use available sensor signals instead of using full-state feedback for active suspension control, LQ static output feedback (SOF) and linear quadratic Gaussian (LQG) controllers are designed with the QC model. From these controllers, observer-based controllers for the FC model are also derived. To verify the performance of the controllers for the FC model derived from LQR and LQ SOF ones designed with the QC model, frequency domain analysis is undertaken. From the analysis, it is confirmed that the controllers for the FC model derived from LQ and LQ SOF ones designed with the QC model can give equivalent performance to those designed with the FC one in terms of ride comfort.

1. Introduction

The main objective in designing suspension is to improve ride comfort and road adhesion. Ride comfort can be evaluated with the vertical acceleration of a sprung mass of a vehicle. Figure 1 shows the frequency weightings representing the sensitivity of the human body to vibration, specified in ISO2631-1 [1]. According to Figure 1, human bodies are considered to be most sensitive to frequencies in the 0.5–2 Hz range in horizontal direction and in the 4–10 Hz range in vertical direction. So, it is necessary to reduce the vertical acceleration of the sprung mass near 4–10 Hz range for the ride comfort [2]. Road adhesion can be evaluated with suspension stroke and tire deflection. Generally, it has been known that there is trade-off between the ride comfort and road adhesion [3].
There have been several devices for controlling the suspension in vehicles, such as active suspension, continuous damping control and air spring, etc. Among them, active suspension has been known to be the most effective for enhancing ride comfort and road adhesion. As a result, there have been a lot of studies on active suspension control. Comprehensive surveys on active suspension control can also be found in [3,4,5,6,7]. Recent advances for the last decade in active suspension control are summarized in [8].
Several suspension models have been used in designing a controller for the active suspension. Typical models are 2-Degree-of-Freedom (DOF) quarter-car (QC), 4-DOF half-car (HC) and 7-DOF full-car (FC) ones [6]. Nonlinear elements or actuators have been included in suspension models [9,10]. A typical example including nonlinear elements and actuators is an in-wheel motor system, which has been developed for the last decade [11,12,13]. To capture nonlinearities in suspension models, a Takagi–Sugeno (T-S) or interval fuzzy model has been adopted [14,15,16]. Among these models, 7-DOF linear FC model have been preferred because it can represent the vertical, roll and pitch motions of a sprung mass or vehicle body [3]. However, the FC model is so complex that it is not easy to implement a controller for active suspension.
Several methodologies such as linear optimal control, H control, nonlinear control, and adaptive control theories have been applied in designing an active suspension controller [3,4,5,6,7,8]. Among these methodologies, linear quadratic regulator (LQR) has been widely adopted with 2-DOF QC and 7-DOF FC models because it can provide a systematic way to design a controller for active suspension. LQR is a full-state feedback controller, which requires all state variables to be measured. However, some of state variables of 7-DOF FC model cannot be measured in real vehicles. For example, the vertical displacement, the roll and pitch angles of a sprung mass are too hard to be measured with commercial sensors. For this reason, a state observer has been adopted for the purpose of estimating state variables of the FC model. Generally, the Kalman filter has been widely adopted as a state observer for LQR. With the Kalman filter and LQR, the linear quadratic Gaussian (LQG) for active suspension control has been extensively studied [17,18,19].
Another drawback in using LQR or LQG with the FC model for active suspension control is that there are so many state variables and control inputs. For example, the FC model has 14 state variables and four control inputs. If LQR is adopted for active suspension control, the dimension of gain matrix is 4 × 14. As a result, it is not easy to implement LQR or LQG on real vehicles. Therefore, it is necessary to design a simple controller, which requires a smaller number of state variables and is easy to implement in real vehicles.
In this paper, a 2-DOF QC model describing the vertical motions of sprung and unsprung masses is used to design a controller for active suspension. LQR is adopted to design the controller with the QC model. Let this controller be denoted as LQRq. LQRq has four elements in the gain matrix. Let the LQR designed with the FC model be denoted as LQRf. LQRf has 56 elements in the gain matrix. Based on the fact that the FC model consists of four QC ones, a full-state feedback controller for the FC model can derived from LQRq with some matrix manipulations. Let this controller be denoted as LQRfq. This is the key contribution of this paper. In LQRfq, the number of distinct gain elements is four, which is borrowed from LQRq. So, LQRfq is much simpler than LQRf and requires a much smaller number of gain elements in the gain matrix. LQRfq can give an equivalent or better performance than LQRf in terms of the ride comfort. This idea was presented in the previous work [20]. In the work, a proportional-integral-derivative (PID) control was adopted with the suspension stroke signal. With this method, it is not easy to evaluate and validate performances of the designed controller. In view of the controller for the FC model, LQRfq is a static output feedback (SOF) or structured controller. Therefore, this controller can be directly designed using LQ SOF control methodology.
Although the number of state variables needed for LQRq is smaller than those for LQRf, it is still not easy to measure some state variables such as vertical displacement of the sprung mass. To cope with the problem, a LQ SOF controller is designed with the QC model and available sensor signals [21,22]. Let this controller be denoted as LQSOFq. Following the identical procedure from LQRq to LQRfq, a full-state feedback controller for the FC model is derived from LQSOFq, designed with the QC model. Let this controller be denoted as LQSOFfq. The number of gain elements needed for LQSOFfq is just 2, which is borrowed from LQSOFq. Moreover, LQSOFfq can give equivalent or better performance than LQRf in terms of the ride comfort.
Another method to cope with the problem is to use an observer-based output feedback controller. LQG is the typical example. In this paper, LQG is designed with the QC and FC models. Let the LQG designed with the QC and FC models be denoted as LQGq and LQGf, respectively. When designing an observer-based output feedback controller for the FC model, the designed controllers from the QC models, i.e., LQRfq and LQSOFfq, are adopted as a full-state feedback controller.
To verify the performance of the controllers derived from ones designed with the QC model, frequency-domain analysis with singular value plot is conducted for the FC model. The full-state feedback, static-output feedback and observer-based output feedback controllers are compared in terms of the ride comfort.
This paper is organized as follows: In Section 2, 2-DOF QC and 7-DOF FC models are presented and modeled as a state-space equation. Based on the QC model, the controller for active suspension is design with LQR. SOF and LQG controllers are also designed with the QC and FC models. In Section 3, frequency-domain analysis with singular value plot is done. The conclusions are given in Section 4.

2. Controller Design for Active Suspension

In this section, the state-space equations for 2-DOF quarter-car and 7-DOF full-car modes following the procedure presented in the previous research [23,24]. In these studies, vector-matrix forms of equations of motions were derived and merged into a single state-space equation. Hereafter, the subscripts q and f mean controllers designed with the QC and FC models, respectively. The subscript fq means a controller for the FC model derived from that designed with the QC model.

2.1. Controller Design with Quarter-Car Model

Figure 2 shows a 2-DOF QC model, which describes the vertical motions of sprung and unsprung masses. As shown in Figure 2, the control input u generated by an active suspension is located between the sprung and unsprung masses, ms and mu. The disturbance is the road profile, zr.
The force acting on suspension is derived as (1). The equations of vertical motions of the sprung and unsprung masses are given in (2). By combining (1) and (2), and changing it into vector-matrix form, (3) is obtained.
f = k s ( z s z u ) b s ( z ˙ s z ˙ u ) + u
{ m s z ¨ s = f m u z ¨ u = f k t ( z u z r )
[ m s 0 0 m u ] [ z ¨ s z ¨ u ] = [ k s k s k s k s k t ] [ z s z u ] + [ b s b s b s b s ] [ z ˙ s z ˙ u ] + [ 1 1 ] u + [ 0 k t ] z r
New vector and matrices are defined as (4) and (5), respectively [23,24]. With those definitions, (3) is rewritten as (6). The state vector of the QC model is defined as (7). With the definitions of new matrices as given in (8), the state-space equation for QC model is obtained as (9).
z = [ z s z u ] T
M q [ m s 0 0 m u ] ,   K q [ k s k s k s k s k t ] ,   B q [ b s b s b s b s ] ,   U q [ 1 1 ] ,   L q [ 0 k t ]
M q z ¨ = K q z + B q z ˙ + U q u + L q z r
x q = [ z z ˙ ] = [ z s z u z ˙ s z ˙ u ] T
A q [ 0 2 × 2 I 2 × 2 M q 1 K q M q 1 B q ] ,   B 1 q [ 0 2 × 1 M q 1 L q ] ,   B 2 q [ 0 2 × 1 M q 1 U q ]
x ˙ q = A q x q + B 1 q z r + B 2 q u
The LQ objective function for active suspension control is given as (10). The weights ρi can be set to the inverse of the square of the maximum allowable value by Bryson’s rule, i.e., ρi = 1/ηi2 [25]. For ride comfort, the weight ρ1 on the vertical acceleration of the sprung mass should be set to higher values while maintaining the other weights constant. On the other hand, for road adhesion or cornering, the weight ρ3 on the tire deflection should be set to higher. Generally, these two objectives, ride comfort and road adhesion, conflict with each other. The regulated output zq is defined as (11) with the state vector xq and the control input u. With (11), the weighting matrices in Jq are obtained as (12). With (12), the LQ objective function (10) can be rewritten into the vector-matrix form, (13). LQR is a controller with the form of full-state feedback, u = −Kqxq, which minimizes Jq. It is easy to obtain Kq from Riccati equation for a given set of weights, i.e., the matrices Qq, Nq and Rq, as given in (12). Let denote this controller Kq as LQRq.
J q = 0 { ρ 1 z ¨ s 2 + ρ 2 ( z s z u ) 2 + ρ 3 z u 2 + ρ 4 u 2 } d t
z q = C q x q + D q u
Q q = C q T C q ,   N q = C q T D q , R q = D q T D q
J q = 0 z q T z q d t = 0 { [ x q u ] T [ Q q N q N q T R q ] [ x q u ] } d t

2.2. Controller Design with Full-Car Model

Figure 3 shows a 7-DOF FC model, which describes the vertical, roll and pitch motions of a sprung mass, and the vertical motions of four unsprung masses. As shown in Figure 3, the indices of a suspension are in order of front left, front right, rear left and rear right ones, represented as ①, ②, ③ and ④, respectively. In the FC model, there are four external disturbances, i.e., the road profiles, zr1, zr2, zr3 and zr4, on the unsprung mass. Actuators located between the sprung and unsprung masses can generate the control inputs, u1, u2, u3 and u4, at each suspension.
The suspension forces in the FC model are derived as (14). In (14), ui is the control input or force acting on the i-th suspension, which is identical to (1). The equations of motion for the sprung and unsprung masses are given in (15) and (16), respectively.
f i = k s i ( z s i z u i ) b s i ( z ˙ s i z ˙ u i ) + u i , i = 1 4
{ m s z ¨ c = f 1 + f 2 + f 3 + f 4 I x ϕ ¨ = t f f 1 t f f 2 + t r f 3 t r f 4 I y θ ¨ = l f ( f 1 + f 2 ) + l r ( f 3 + f 4 )
{ m u 1 z ¨ u 1 = f 1 + k t 1 ( z u 1 z r 1 ) m u 2 z ¨ u 2 = f 2 + k t 2 ( z u 2 z r 2 ) m u 3 z ¨ u 3 = f 3 + k t 3 ( z u 3 z r 3 ) m u 4 z ¨ u 4 = f 4 + k t 4 ( z u 4 z r 4 )
In (14), the vertical displacement of each corner, zsi, is not a state variable. Therefore, it should be derived from state variables. Figure 4 shows the roll and pitch motions of the sprung mass. As shown in Figure 4, the vertical displacements of the sprung mass, i.e., zs1, zs2, zs3 and zs4, at each corner are derived as (17) from the geometrical relationship under the assumption that the sprung mass is a rigid body. As shown in (17), there is a sine function, which is nonlinear. Generally, the roll and pitch angles are less than 5 deg for most cases of vehicle driving conditions. Under the condition that the roll and pitch angles are small or less than 15 deg, i.e., sin θ θ and sin ϕ ϕ , (17) can be rewritten as the vector-matrix form of (18) [23,24]. In (18), the matrix G represents the geometric relationship between the vertical motions of four corners and the vertical, roll and pitch ones of the sprung mass.
{ z s 1 = z c + t f sin ϕ l f sin θ z s 2 = z c t f sin ϕ l f sin θ z s 3 = z c + t f sin ϕ + l r sin θ z s 4 = z c t f sin ϕ + l r sin θ
[ z s 1 z s 2 z s 3 z s 4 ] = [ 1 t f l f 1 t f l f 1 t r l r 1 t r l r ] G T [ z c ϕ θ ] = G T [ z c ϕ θ ]
For further derivation, the vectors of dynamic variables and the matrices of parameters are defined as (19) and (20), respectively [23,24]. In (20), diag(a,b,c,…) represents the diagonal matrix consisting of the elements, a,b,c,… With those definitions, the Equations (15) and (16) are represented as the vector-matrix forms, (21). Moreover, the Equations (18) and (14) are represented as the vector-matrix forms, (22) and (23), respectively.
{ z s [ z s 1 z s 2 z s 3 z s 4 ] T ,   z u [ z u 1 z u 2 z u 3 z u 4 ] T ,   z r [ z r 1 z r 2 z r 3 z r 4 ] T   p [ z c ϕ θ ] T ,   u f [ u 1 u 2 u 3 u 4 ] T ,   f [ f 1 f 2 f 3 f 4 ] T
{ M s diag ( m s , I x , I y ) ,   M u diag ( m u 1 , m u 2 , m u 3 , m u 4 ) K s diag ( k s 1 , k s 2 , k s 3 , k s 4 ) ,   K t diag ( k t 1 , k t 2 , k t 3 , k t 4 ) B s diag ( b s 1 , b s 2 , b s 3 , b s 4 )
{ M s p ¨ = G f M u z ¨ u = f + K t ( z u z r )
z s = G T p
f = K s ( z s z u ) B s ( z ˙ s z ˙ u ) + u f   = K s ( G T p z u ) B s ( G T p ˙ z ˙ u ) + u f
By replacing f of (21) with (23), (21) is converted into (24). (24) is rearranged into (25) as the vector-matrix form. The vectors and matrices in (25) are defined as (26) and (27), respectively. With those definitions, (25) is rewritten into (28).
{ M s p ¨ = G K s ( G T p z u ) G B s ( G T p ˙ z ˙ u ) + G u f M u z ¨ u = K s ( G T p z u ) + B s ( G T p ˙ z ˙ u ) + K t ( z u z r ) u f
[ M s 0 0 M u ] [ p ¨ z ¨ u ] = [ G K s G T G K s K s G T K s + K t ] [ p z u ] + [ G B s G T G B s B s G T B s ] [ p ˙ z ˙ u ] + [ G I ] u f + [ 0 K t ] z r
q [ p z u ] ,   w f z r
{ M e [ M s 0 0 M u ] ,   K e [ G K s G T G K s K s G T K s + K t ] B e [ G B s G T G B s B s G T B s ] ,   U e [ G I ] ,   L e [ 0 K t ]
M e q ¨ = K e q + B e q ˙ + U e u f + L e w f
The state vector of the FC model is defined as (29). With the definitions of matrices as given in (30), the state-space equation of the 7-DOF FC model is obtained as (31).
x f [ q q ˙ ]
A f = [ 0 7 × 7 I 7 × 7 M e 1 K e M e 1 B e ] ,   B 1 f = [ 0 7 × 4 M e 1 L e ] ,   B 2 f = [ 0 7 × 4 M e 1 U e ]
x ˙ f = A f x f + B 1 f w f + B 2 f u f
The LQ objective function for active suspension control with the FC model is given as (32). The weights ζi can be set to the inverse of the square of the maximum allowable value by Bryson’s rule, i.e., ζi = 1/ξi2 [25]. For ride comfort, the weights ξ1, ξ2, and ξ3 on the vertical, roll and pitch accelerations should be set to higher values. On the other hand, for road adhesion or cornering, the weights ξ8 and ξ9 on the suspension stroke and tire deflection should be set to higher. Generally, these two objectives, ride comfort and road adhesion, conflict with each other.
The regulated output zf is defined as (33) with the state and input vectors, xf and uf. With (33), the weighting matrices in Jf are obtained as (34). With (34), the LQ objective function (32) can be rewritten into the vector-matrix form of (35). LQR is a controller with the form of full-state feedback, uf = −Kfxf, which minimizes Jf. It is easy to obtain Kf from Riccati equation for a given set of weights, i.e., the weighting matrices Qf, Nf and Rf, as given in (34). Let denote this controller Kf as LQRf1.
J f = 0 { ξ 1 z ¨ c 2 + ξ 2 ϕ ¨ 2 + ξ 3 θ ¨ 2 + ξ 4 ϕ 2 + ξ 5 ϕ ˙ 2 + ξ 6 θ 2 + ξ 7 θ ˙ 2 + ξ 8 i = 1 4 ( z s i z u i ) 2 + ξ 9 i = 1 4 z u i 2 + ξ 10 i = 1 4 u i 2 } d t
z f = C f x f + D f u f
Q f = C f T C f ,   N f = C f T D f , R f = D f T D f
J f = 0 z f T z f d t = 0 { [ x f u f ] T [ Q f N f N f T R f ] [ x f u f ] } d t
The LQ objective function for the FC model, (32), can be converted into that for the QC one, (36), just like (10). This function, (36), has the identical form to (10). The weights σi can be set to the inverse of the square of the maximum allowable value by Bryson’s rule, i.e., σi = 1/χi2 [25]. With the new definitions of matrices as given in (37), the output and the weighting matrices are defined as (38) and (39), respectively. With these matrices, the LQ objective function (36) is represented as the vector-matrix form of (40). As mentioned earlier, it is easy to obtain Kf from the Riccati equation for a given set of weights, i.e., the matrices Qfq, Nfq and Rfq, as given in (39). Let this controller be denoted as LQRf2. LQRf1 and LQRf2 are the full-state feedback controllers, which has different LQ objective functions, (32) and (36).
J f q = 0 { σ 1 i = 1 4 z ¨ s i 2 + σ 2 i = 1 4 ( z s i z u i ) 2 + σ 3 i = 1 4 z u i 2 + σ 4 i = 1 4 u i 2 } d t
{ C z s a G T [ A f , 8 A f , 9 A f , 10 ] T C z s b G T [ B 2 f , 8 B 2 f , 9 B 2 f , 10 ] T C s s q [ G T I 4 × 4 0 4 × 7 ] C z u q [ 0 3 × 4 I 4 × 4 0 4 × 7 ]
z f q = C f q x f + D f q u f   = [ C z s a C s s q C z u q 0 4 × 14 ] x f + [ C z s b 0 4 × 4 0 4 × 7 I 4 × 4 ] u f
Q f q = C f q T C f q ,   N f q = C f q T D f q , R f q = D f q T D f q
J f q = 0 z f q T z f q d t = 0 { [ x f u f ] T [ Q f q N f q N f q T R f q ] [ x f u f ] } d t
LQR provides the systematic way to design the active suspension controller. Generally, LQR has the structure of full-state feedback. As shown in (29) and (19), i.e., xf and uf, there are 14 state variables and 4 control inputs in the 7-DOF FC model, respectively. So, the dimension of the gain matrix Kf of LQR is 4 × 14. To implement the LQR in real vehicles, it is necessary to precisely measure or estimate 14 state variables. However, it is too hard to measure or estimate those variables in real vehicles. For this reason, it is necessary to design a controller which has the smaller number of state variables and can be used for the FC model.

2.3. How to Use LQR for Quarter-Car Model as a Controller for Full-Car One

As shown in Figure 2 and Figure 3, it can be regarded that the FC model consists of four QC models. The vertical displacements of the sprung and unsprung masses, zsi and zui, in the FC model are identical to those of the QC one. Therefore, LQRq can be used to design a full-state feedback controller having the identical structure to LQRf. This subsection explains how to use LQRq in designing the controller.
Some of the state variables in the FC model can be grouped into those of the QC one, xqi, as given in (41). Kq in (42) represents the gain matrix of LQRq. By multiplying (42) and (41) together, the control input of active suspension at each corner is obtained as (43). With Kq in (42), the control input for FC model, uf, is calculated as (44), which has the identical form to (43). (44) is converted into (45) with the definition of (22). The four vectors in (45) can be represented by the state vector of the FC model, xf, as given in (46). With (46), the control input for the FC model is represented as (47) with the gain elements of Kq. This means that the full-state feedback controller for the FC model can be derived from the gain elements of Kq. As shown in (47), the number of gain elements needed to derive the full-state feedback controller for the FC model is just four, which is much smaller than that of Kf. This is the key contribution of this paper.
As given in (48), the gain matrix H has a special structure. Therefore, it is different from the LQR gain matrix, Kf. Let this controller H be denoted as LQRfq. It is expected that LQRfq can give an equivalent performance to LQRf2 because these adopt the identical LQ objective function, (10) and (36).
x q i = [ z s i z u i z ˙ s i z ˙ u i ] T , i = 1 , 2 , 3 , 4
K q = [ k 1 k 2 k 3 k 4 ]
u i = K q x q i = ( k 1 z s i + k 2 z u i + k 3 z ˙ s i + k 4 z ˙ u i ) ,   i = 1 , 2 , 3 , 4
u f = k 1 z s + k 2 z u + k 3 z ˙ s + k 4 z ˙ u
u f = k 1 G T p + k 2 z u + k 3 G T p ˙ + k 4 z ˙ u
{ p = C 1 x f , C 1 [ I 3 × 3 0 3 × 11 ] p ˙ = C 2 x f , C 2 [ 0 3 × 7 I 3 × 3 0 3 × 4 ] z u = C 3 x f , C 3 [ 0 4 × 3 I 4 × 4 0 4 × 7 ] z ˙ u = C 4 x f , C 4 [ 0 4 × 10 I 4 × 4 ]
u f = k 1 G T C 1 x f + k 2 C 3 x f + k 3 G T C 2 x f + k 4 C 4 x f   = ( k 1 G T C 1 + k 2 C 3 + k 3 G T C 2 + k 4 C 4 ) x f   = H x f
H = [ r 1 r 1 r 1 r 1 r 2 r 2 r 2 r 2 r 3 r 3 r 4 r 4 r 5 I 4 × 4 r 6 r 6 r 6 r 6 r 7 r 7 r 7 r 7 r 8 r 8 r 9 r 9 r 10 I 4 × 4 ]
Indeed, the controller with the form of (47), i.e., LQRfq, is LQ static output feedback (SOF) or structured controller (SC). LQRfq is LQ SOF controller in that some of state variables, i.e., available sensor signals, are used for the feedback. It is also structured controller because H has a special structure, as given in (48). In (48), ri are distinct gain elements. All the controllers derived from (47) have the identical structure of (48). Although there are 10 gain elements in H, only four gain ones are needed to calculate it, as given in (47).
In (42) and (47), the controller gains, k1, k2, k3 and k4, are already determined from LQRq. It should be noticed that LQRq has the LQ objective function of (10). Instead of using LQRq to derive H, the four controller gains in H can be directly determined by the optimization method as presented in the previous work [26].
Let the controller gain matrix be defined as Sq as (49). With Sq, the full-state feedback gain matrix Sf for the FC model is obtained as (50). As given in (50), the gain matrix Sf has the form identical to (47). In other words, the full-state feedback controller for the FC model can be designed with four gain elements in Sq. The problem of determining the controller gains of Sq can be formulated as an optimization one, as given in (51). In (51), Re[] is the function that returns real parts of complex numbers. It should be noticed that the LQ objective function of (36) or the weighting matrices of (39) is adopted for this optimization problem, (51). This problem is non-convex and nonlinear. So, it is hard to apply a convex optimization method. To solve this problem, the heuristic optimization method, i.e., the derandomized evolution strategy with covariance matrix adaptation (CMA-ES), is adopted [27]. With the optimum gains obtained by solving (51), the full-state feedback gain matrix Sf is obtained as (50). Let denote this controller Sf as LQSCfq. It is expected that LQSCfq is equivalent to LQRf2 and LQRfq because these have the identical LQ objective function, (36).
S q = [ q 1 q 2 q 3 q 4 ]
u f = ( q 1 G T C 1 + q 2 C 3 + q 3 G T C 2 + q 4 C 4 ) x f = S f x f
min q 1 , q 2 , q 3 , q 4 J q s = 1 2 trace ( P f ) , P f = P f T > 0 s . t . { max ( Re [ A f + B 2 f S f ] ) < 0 ( A f + B 2 f S f ) T P f + P f ( A f + B 2 f S f ) + Q f q + S f T N f q T + N f q S f + S f T R f q S f = 0

2.4. Static Output Feedback Control with Quarter-Car Model

When using LQR for active suspension control, it is too hard to measure the state variables. For example, the state variables in the QC model, i.e., the vertical displacements and velocities of the sprung and unsprung masses are hard to measure. To cope with this problem, SOF is adopted. SOF uses available sensor signals for feedback control. For the QC model, the typical sensor signals are the suspension stroke and its rate.
The SOF controller has the form of (52). With the definition of the state vector, (7), the vector of sensor outputs, yq, in the QC model are defined as (53). As shown in (53), the sensor outputs are the suspension stroke and its rate. Therefore, there are two elements in KSOF while Kq has four elements. By replacing yq in (52) with the output definition of (53), the control input u is obtained as (54). From (54), the gain matrix of the full-state feedback from LQ SOF controller for the QC model is given in (55). Let this controller be denoted as LQSOFq.
u = K S O F y q
y q = [ z s z u z ˙ s z ˙ u ] = [ 1 1 0 0 0 0 1 1 ] x q   = C s x q
u = K S O F y q = K S O F C s x q
V q K S O F C s
LQ SOF control is to find KSOF which gives a minimum of LQ objective function. So, this is formulated as an optimization problem, as given in (56). To find the optimum gain KSOF, the heuristic optimization method, CMA-ES, is applied [27]. After finding the gain matrix Vq or KSOF, the active suspension controller for the FC model can be obtained as (57). In (57), v1, v2, v3 and v4 are the elements of Vq. Let denote this controller as LQSOFfq.
min K S O F J f s = 1 2 trace ( P s ) , P s = P s T > 0 s . t . { max ( Re [ A q + B 2 q V q ] ) < 0 ( A q + B 2 q V q ) T P s + P s ( A q + B 2 q V q ) + Q q + V q T N q T + N q V q + V q T R q V q = 0
u f = ( v 1 G T C 1 + v 2 C 3 + v 3 G T C 2 + v 4 C 4 ) x f = W x f
Equation (57) has the identical form to (47) and (50). In other words, LQRfq, LQSCfq, and LQSOFfq have the identical form. The difference between these controllers is how to obtain gain elements. The four gain elements of LQRfq are identical to those of LQRq. The four gain elements of LQSCfq and LQSOFfq are obtained from the optimization problem (51) and (56), respectively. The difference between LQSCfq and LQSOFfq is the model used in the optimization. While LQSCfq uses the FC model and the weighting matrices of Qfq, Nfq and Rfq, LQSOFfq uses the QC one with the weighting matrices Qq, Nq and Rq.

2.5. Linear Quadratic Gaussian Control with Quarter-Car Model

As mentioned earlier, it is not easy to measure some state variables in LQR. Instead of using SOF, an observer-based controller, especially LQG, has been adopted. In this section, LQG is designed with QC and FC models. In LQG, the gain of LQR for QC model is used to design the full-state feedback controller, as given in (47).
In this paper, it is assumed that the available sensor outputs are the suspension stroke and its rate in the QC model. So, the vector of available sensor outputs is given in (53), which is identical to that of the LQ SOF controller. The state observer or Kalman filter for the QC model has the form of (58). The covariance matrices of system uncertainty and measurement noise are given as Mq and Nq, respectively. The gain matrix of the observer, Lq, is calculated by solving the filter algebraic Riccati equation (FARE) with the matrices Aq, Cs, Mq and Nq. With the estimated state and the controller gain Kq, the control input for the QC model is obtained as (59). The state-space equation of the controlled system and the state observer for the QC model is given as (60). This is called LQG or compensator [28]. Let this controller (60) be denoted as LQGq.
x ^ ˙ q = A q x ^ q + B 2 q u + B 1 q z r + L q ( y q C s x ^ q )
u = K q x ^ q
[ x ˙ q x ^ ˙ q ] = [ A q B 2 q K q L q C s A q B 2 q K q L q C s ] [ x q x ^ q ]
For the FC model, the vector of available sensor outputs, i.e., the suspension stroke and its rate, given in (53) is represented as (61). In (61), Af,i represents the i-th row of the matrix Af. The state observer or Kalman filter for the FC model has the form of (62). The covariance matrices of system uncertainty and measurement noise are given as Mf and Nf, respectively. The gain matrix of the observer, Lf, is calculated by solving FARE with the matrices Af, Ce, Mf and Nf. With the estimated state and the controller gain Kf, the control input for the FC model is obtained as (63). The state-space equation with the controlled system and the observer for the FC model is given as (64). This is LQG for the FC model.
Let the LQG with the LQ objective functions of (35) and (40), i.e., LQRf1 and LQRf2, be denoted as LQGf1 and LQGf2, respectively. If the controller gain Kf is replaced with H in (64), the controller gain Kq can be used for LQG with FC model, as given in (65). Let the LQG of (65) be denoted as LQGfq. If the controller gain H is replaced with Sf in (65), the controller gain Sq can be used for LQG with FC model. Let this LQG be denoted as LQGSCfq. It is expected that LQGf2 is equivalent to LQGfq and LQGSCfq because these have the identical LQ objective function and the identical sensor signals.
y f = [ z s z u z ˙ s z ˙ u ] = [ G T p z u G T p ˙ z ˙ u ]   = [ G T I 4 × 4 0 4 × 7 0 4 × 7 G T I 4 × 4 ] x f   = C e x f
x ^ ˙ f = A f x ^ f + B 2 f u f + B 1 f w f + L f ( y f C e x ^ f )
u f = K f x ^ f
[ x ˙ f x ^ ˙ f ] = [ A f B 2 f K f L f C e A f B 2 f K f L f C e ] [ x f x ^ f ]
[ x ˙ f x ^ ˙ f ] = [ A f B 2 f H L f C e A f + B 2 f H L f C e ] [ x f x ^ f ]

3. Frequency-Domain Analysis on Designed Controllers

In this section, frequency-domain analysis is done to check the performance of the designed controllers. Through frequency-domain analysis, full-state feedback, static output feedback and observer-based output feedback controllers, i.e., LQR, LQSOF, and LQG, are compared to one another.
Table 1 shows the parameters of 2-DOF QC and 7-DOF FC models, which were referred from E-Class sedan given in CarSim [29]. Weights in LQ objective functions, (10), (32), and (36), are given in Table 2. Those weights emphasize the ride comfort, i.e., the reduction of the vertical acceleration of the sprung mass. As a consequence, the road adhesion will deteriorate. For this reason, frequency responses within the range from 4 to 10 Hz should be focused, as given in Figure 1. In case of the roll and pitch angles, frequency responses within the range from 0.5 to 2 Hz should be focused. In this paper, it is assumed that the actuator used to generate the control input has infinite bandwidth. In other words, delay and saturation of an actuator are not taken into account. The covariance matrices of system uncertainty and measurement noise used in designing LQG are given in Table 3.

3.1. Full-State Feedback Case

The first type of frequency-domain analysis is done for the full-state feedback controllers. For frequency-domain analysis, the FC model is used. Figure 5 shows the frequency responses of each output from the road profile zr1 with the FC model for each controller. The gain matrices of Kq and Sq of LQRq and LQSCfq are given in Table 4. LQRq was obtained by solving Riccati equation with the state-space Equation (9), the LQ objective function (10) and Table 2. LQSCfq was obtained by solving (51) with the state-space Equation (31), LQ objective function (36) and Table 2. As expected in the Section 2.3, LQRfq and LQSCfq has nearly identical elements as given in Table 4.
As shown in Figure 5a, LQRf2 is better than LQRf1 in terms of the ride comfort. This means that higher emphasis on the vertical acceleration at each corner in (36) is better than higher emphasis on the vertical, roll and pitch accelerations in (32). On the other hand, LQRf1 is better than LQRf1 in terms of the roll and pitch angles. This is caused by the fact that the LQ objective function of LQRf2, (36), has no terms on the roll and pitch angles.
As shown in Figure 5, LQRf2, LQRfq and LQSCfq show nearly identical performances in terms of the ride comfort and the road adhesion. The measures of the ride comfort, i.e., the vertical acceleration between 4 and 10 Hz, the roll and pitch angle between 0.5 and 2 Hz are improved by the controllers, as shown in Figure 5a–c, respectively. This is natural because these were designed with the identical LQ objective functions. From these results, it can be concluded that the controllers derived from LQRq, i.e., LQRfq and LQSCfq, can give equivalent performance to that designed with the FC one, LQRf2. This is quite important that a much smaller number of gains are needed in designing the active suspension controller for the FC model. In other words, only four gain elements of LQRq can give nearly identical performance to LQRf2 which needs 56 gain elements in the active suspension controller for the FC model.
Another feature of the controllers is that the suspension strokes of the controlled suspensions are nearly identical to one another within the frequency range from 0.1 to 6 Hz, and larger than the passive suspension as shown in Figure 5d. Moreover, the tire deflections of these controllers were smaller than that of the passive system, as shown in Figure 5e. This means that the active suspension controllers enhance the ride comfort by generating larger motion of the suspension stroke, and that the road adhesion was deteriorated by those controllers as a result. This is a natural consequence of the active suspension control. Hereafter, the frequency responses of the suspension stroke and the tire deflection are not presented because it will show identical tendency regardless of controllers.
Table 5 shows the decibel scale magnitudes of frequency responses for each controller. The values in Table 5 are calculated from the frequency responses in Figure 5. The magnitudes of the vertical acceleration, the suspension stroke and the tire deflection are sampled at 6 Hz, which is the peak frequency of Wk in Figure 1. The magnitudes of the roll and pitch angles are sampled at 1 Hz, which is the peak frequency of Wd in Figure 1. As shown Table 5, the LQRf2, LQRfq and LQSCfq show nearly similar performance to one another. Moreover, the ride comfort enhancement achieved by these controllers are better than LQRf1. For example, LQRfq reduced the vertical acceleration from 38.9 dB and 29.9 dB the passive system and LQRf1 to 23.7 dB, respectively. These results mean that the vertical acceleration of the passive system and LQRf1 were reduced by 39% and 20% by virtue of LQRfq, respectively. Among these controllers, LQSCfq shows the best performance in terms of the vertical acceleration. On the other hand, the roll and pitch angles of LQRf1 are better than the three controllers because the LQ objective function of LQRf1 has the terms emphasizing the roll and pitch angles. As shown in Figure 5 and Table 5, the suspension stroke increases, and the tire deflection decreases as the vertical acceleration is reduced by the controllers. The smaller the tire deflection, the smaller the vertical and longitudinal/lateral tire forces. As a result, the road adhesion is deteriorated.

3.2. Static Output Feedback Case

The second type of frequency-domain analysis is done for the LQ SOF controllers. Figure 6 shows the frequency responses of each output from the road profile zr1 with the FC model for each controller. The gain matrices of KSOF and Vq of LQSOFq are given in Table 6. The gain KSOF of LQSOF is obtained by solving (56) with the LQ objective function (10) and Table 2.
As shown in Figure 6, all the measures of LQSOFfq designed with the QC model are nearly identical to LQRf2 and LQRfq within the frequency range from 2 to 10 Hz. This means that only two gain elements of KSOF can give the equivalent performance to LQRf2 and LQRfq, the active suspension controller designed with the FC model. The first reason why LQSOFfq gives good performance equivalent to LQRf2 and LQRfq is that the controller structure of LQSOFfq is identical to the suspension force (1). In other words, the first and second gain elements of KSOF correspond to ks and bs of (1), respectively. If those of KSOF coincide to ks and bs of (1), the forces generated from tire cannot be transmitted to the sprung mass. The second reason is that these controllers have identical LQ objective functions of (10) and (36).
As shown in Figure 6, the drawback of LQSOFfq is that its performance is deteriorated in the frequency range below 1 Hz. As shown in Figure 1, the motion sickness near 0.2 Hz and the horizontal motion near 1 Hz are important in terms of the ride comfort. The vertical acceleration of the sprung mass of LQSOFfq is worse than the passive system below 0.5 Hz, as shown in Figure 6a. Moreover, the roll and pitch angles of the sprung mass of LQSOFfq are worse than the passive system below 0.7 Hz, as shown in Figure 6b,c. This is critical because it deteriorates the motion sickness and the ride comfort along the horizontal direction. For this reason, LQSOFfq is not recommended for the ride comfort in spite of its good performance along the vertical direction and the simplicity in controller design procedure.

3.3. LQG Case

The third type of frequency-domain analysis is done for the observer-based controllers, i.e., LQG. Figure 7 shows the frequency responses of each output from the road profile zr with the QC model for LQRq, LQSOFq and LQGq, as given in (60).
As shown in Figure 7a, the performance of LQGq is slightly worse than that of LQRq and LQSOFq. Moreover, it becomes worse than LQRq and LQSOFq below 4 Hz and above 40 Hz in terms of the ride comfort. This tendency of LQGq is nearly same as that of LQSOFq, as given in Figure 6. This is caused by the fact that LQSOFq and LQGq use identical sensor signals. Nevertheless, the performance of LQGq is satisfactory in terms of the ride comfort because it shows nearly identical performance to LQRq within the frequency range from 4 to 10 Hz. Moreover, the performance of LQGq is better than LQSOFq for the frequency range below 0.6 Hz, as shown in Figure 7a. The frequency responses of the suspension stroke and the tire deflection in Figure 7b,c show nearly identical tendency to those of Figure 5d,e, respectively.
These results are varied according to the selection of the covariance matrices, i.e., Mq and Nq under the condition that an identical set of weights in LQ objective function is used. It also depends on the selection of the available sensor outputs, as given in (53). If LQSOFq and LQGq use identical sensor signals for output feedback, these can give nearly identical performances, as shown in Figure 7. From these results, it can be concluded that LQG can be adopted as an output feedback controller for the purpose of enhancing the ride comfort.
Figure 8 shows the frequency responses of each output from the road profile zr1 with the FC model for each controller. As shown in Figure 8, the relationship between LQRq and LQGq given in Figure 7 is nearly identical to that of LQRfq and LQGfq. In other words, the performance of LQGfq becomes worse than LQRfq below 4 Hz and above 40 Hz in terms of the ride comfort. This can be checked in Figure 7a and Figure 8a. Moreover, the ride comfort of LQGq becomes identical to the passive system below 0.5 Hz, as shown in Figure 7a. This is also confirmed by the frequency responses of the passive system and LQGfq in Figure 8a. Nevertheless, the LQG controllers, i.e., LQGf2, LQGfq and LQGSCfq, show good performance within the frequency range from 2 to 10 Hz, in terms of the ride comfort.
As shown in Figure 8, the frequency responses of LQGf2, LQGfq and LQGSCfq show identical tendency as LQRf2, LQRfq and LQSCfq, as given in Figure 5. This is confirmed by the fact that the frequency responses of these controllers are nearly identical to each other. This is quite natural because these have the identical LQ objective function and the identical Kalman filter. Moreover, LQGfq and LQGSCfq have the identical controller structure and nearly identical controller gains, as given in Table 4.
Figure 7 did not show the frequency responses of the roll and pitch angles because the QC model has no roll and pitch motions. As shown in Figure 8b,c, the LQG controllers show identical performance to the passive system below 0.6 Hz. This means that the LQG controllers, i.e., LQGf2, LQGfq and LQGSCfq, are not effective in controlling the roll and pitch motions for the ride comfort along the horizontal direction and the motion sickness.

4. Conclusions

In this paper, the controllers for active suspension were designed with the 2-DOF QC and 7-DOF FC models. Two LQRs, i.e., LQRq and LQRf2, were designed with the QC and FC models based on the identical LQ objective functions. To cope with the problem in measuring the state variables of the FC model, LQRq designed with the QC model was used to derive the full-state feedback controller for the FC model, i.e., LQRfq. In deriving LQRfq, only four gain elements of LQRq are needed. To cope with the identical problem, SOF and LQG controllers, i.e., LQSOFq and LQGq, are designed with the QC models and used to derive the controllers for the FC model, i.e., LQSOFfq and LQGfq. In deriving LQSOFfq, only two gain elements of LQSOFq are needed. To verify the fact that the full-state feedback controllers, i.e., LQRfq and LQGfq, derived from LQRq can give equivalent or better performance than LQRf2 and LQGf2 designed with the FC model in terms of the ride comfort, frequency-domain analysis was conducted for full-state feedback, static output feedback and observer-based output feedback controllers. From the analysis, it can be concluded that the controllers with much simpler structure or much smaller gain elements designed with the QC model, i.e., LQRfq, LQSOFfq and LQGfq, can give equivalent performance to LQRf2. The vertical acceleration of the passive system and LQRf1 were reduced by 39% and 20% by virtue of LQRfq, respectively. Further research can include comparison between robust control methodologies in designing active suspension controllers with the QC model for the FC one. Moreover, it can also include test results on a vehicle simulation package and experimental validation on real vehicles.

Author Contributions

Conceptualization, M.P. and S.Y.; methodology, S.Y.; software, S.Y.; validation, M.P. and S.Y.; formal analysis, S.Y.; investigation, M.P.; resources, M.P.; data curation, M.P.; writing—original draft preparation, S.Y.; writing—review and editing, M.P.; visualization, S.Y.; supervision, M.P.; project administration, M.P.; funding acquisition, S.Y. 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

All the experimental data are generated during the study at a local computer. Therefore, it cannot be made publicly available for other researchers.

Acknowledgments

This study was supported by the Research Program funded by the Seoul National University of Science and Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ISO 2631-1. Mechanical Vibration and Shock–Evaluation of Human Exposure to Whole-Body Vibration–Part 1: General Requirements; International Organization for Standardization: Geneva, Switzerland, 1997. [Google Scholar]
  2. Rimell, A.N.; Mansfield, N.J. Design of digital filters for frequency weightings required for risk assessments of workers exposed to vibration. Ind. Health 2007, 45, 512–519. [Google Scholar] [CrossRef] [Green Version]
  3. Tseng, H.E.; Hrovat, D. State of the art survey: Active and semi-active suspension control. Veh. Syst. Dyn. 2015, 53, 1034–1062. [Google Scholar] [CrossRef]
  4. Sharp, R.S.; Crolla, D.A. Road vehicle suspension system design—A review. Veh. Syst. Dyn. 1987, 16, 167–192. [Google Scholar] [CrossRef]
  5. Hrovat, D. Survey of advanced suspension developments and related optimal control applications. Automatica 1997, 33, 1781–1817. [Google Scholar] [CrossRef]
  6. Cao, D.; Song, X.; Ahmadian, M. Editors’ perspectives: Road vehicle suspension design, dynamics, and control. Veh. Syst. Dyn. 2011, 49, 3–28. [Google Scholar] [CrossRef] [Green Version]
  7. Theunissen, J.; Tota, A.; Gruber, P.; Dhaens, M.; Sorniotti, A. Preview-based techniques for vehicle suspension control: A state-of-the-art review. Annu. Rev. Control 2021, 51, 206–235. [Google Scholar] [CrossRef]
  8. Al-Ashmori, M.; Wang, X.A. Systematic literature review of various control techniques for active seat suspension systems. Appl. Sci. 2020, 10, 1148. [Google Scholar] [CrossRef] [Green Version]
  9. Sun, W.; Pan, H.; Zhang, Y.; Gao, H. Multi-objective control for uncertain nonlinear active suspension systems. Mechatronics 2014, 24, 318–327. [Google Scholar] [CrossRef]
  10. Chen, S.-A.; Wang, J.-C.; Yao, M.; Kim, Y.-B. Improved optimal sliding mode control for a non-linear vehicle active suspension system. J. Sound Vib. 2017, 395, 1–25. [Google Scholar] [CrossRef]
  11. Wang, R.; Jing, H.; Yan, F.; Karimi, H.R.; Chen, N. Optimization and finite-frequency H control of active suspensions in in-wheel motor driven electric ground vehicles. J. Frankl. Inst. 2015, 352, 468–484. [Google Scholar] [CrossRef]
  12. Liu, M.; Gu, F.; Zhang, Y. Ride comfort optimization of in-wheel-motor electric vehicles with in-wheel vibration absorbers. Energies 2017, 10, 1647. [Google Scholar] [CrossRef] [Green Version]
  13. 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]
  14. Li, H.; Yu, J.; Hilton, C.; Liu, H. Adaptive sliding-mode control for nonlinear active suspension vehicle systems using T-S fuzzy approach. IEEE Trans. Ind. Electron. 2013, 60, 3328–3338. [Google Scholar] [CrossRef]
  15. Yatak, M.O.; Sahin, F. Ride comfort-road holding trade-off improvement of full vehicle active suspension system by interval type-2 fuzzy control. Eng. Sci. Technol. Int. J. 2021, 24, 259–270. [Google Scholar] [CrossRef]
  16. Mrazgua, J.; Chaibi, R.; Tissir, E.H.; Ouahi, M. Static output feedback stabilization of T-S fuzzy active suspension systems. J. Terramech. 2021, 97, 19–27. [Google Scholar] [CrossRef]
  17. Wilson, D.A.; Sharp, R.S.; Hassan, S.A. The application of linear optimal control theory to the design of active automobile suspensions. Veh. Syst. Dyn. 1987, 15, 105–118. [Google Scholar] [CrossRef]
  18. Hac, A. Optimal linear preview control of active vehicle suspension. Veh. Syst. Dyn. 1992, 21, 167–195. [Google Scholar] [CrossRef]
  19. Abdel-Hady, M.B.A.; Crolla, D.A. Active suspension control algorithms for a four wheel vehicle model. Int. J. Veh. Des. 1992, 13, 144–158. [Google Scholar]
  20. Mudduluru, S.R.; Chizari, M. Quarter and full car models optimisation of passive and active suspension system using genetic algorithm. arXiv 2021, arXiv:2101.12629. Available online: https://arxiv.org/abs/2101.12629 (accessed on 30 August 2021).
  21. Camino, J.F.; Zampieri, D.E.; Peres, P.L.D. Design of a vehicular suspension controller by static output feedback. In Proceedings of the American Control Conference, San Diego, CA, USA, 2–4 June 1999; pp. 3168–3171. [Google Scholar]
  22. Elmadany, M.M.; Al-Majed, M.I. Quadratic synthesis of active controls for a quarter-car model. J. Vib. Control 2001, 7, 1237–1252. [Google Scholar] [CrossRef]
  23. Lu, J.; DePoyster, M. Multiobjective optimal suspension control to achieve integrated ride and handling performance. IEEE Trans. Control Syst. Technol. 2002, 10, 807–821. [Google Scholar]
  24. Jung, Y.H.; Choi, J.W.; Seo, Y.B. Overlapping decentralized EA control design for an active suspension system of a full car model. In Proceedings of the 39th SICE Annual Conference, Iizuka, Japan, 28 July 2000; pp. 85–90. [Google Scholar]
  25. Bryson, A.E., Jr.; Ho, Y. Applied Optimal Control; Hemisphere: New York, NY, USA, 1975. [Google Scholar]
  26. Yim, S. Design of a robust controller for rollover prevention with active suspension and differential braking. J. Mech. Sci. Technol. 2012, 26, 213–222. [Google Scholar] [CrossRef]
  27. Hansen, N.; Muller, S.D.; Koumoutsakos, P. Reducing the time complexity of the derandomized evolution strategy with covariance matrix adaptation (CMA-ES). Evol. Comput. 2003, 11, 1–18. [Google Scholar] [CrossRef] [PubMed]
  28. Stengel, R.F. Stochastic Optimal Control: Theory and Application; A Wiley-Interscience Publication: New York, NY, USA, 1986. [Google Scholar]
  29. Mechanical Simulation Corporation. CarSim Data Manual; Version 8; Mechanical Simulation Corporation: Ann Arbor, MA, USA, 2009. [Google Scholar]
Figure 1. Frequency weightings specified in ISO 2631-1 standard [2].
Figure 1. Frequency weightings specified in ISO 2631-1 standard [2].
Energies 14 08231 g001
Figure 2. 2-DOF quarter-car model.
Figure 2. 2-DOF quarter-car model.
Energies 14 08231 g002
Figure 3. 7-DOF full-car model.
Figure 3. 7-DOF full-car model.
Energies 14 08231 g003
Figure 4. Roll and pitch motions of the sprung mass. (a) Roll motion (front view); (b) Pitch motion (side view).
Figure 4. Roll and pitch motions of the sprung mass. (a) Roll motion (front view); (b) Pitch motion (side view).
Energies 14 08231 g004
Figure 5. Frequency responses from the road profile zr1 to each output for full-state feedback case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass; (d) Suspension stroke; (e) Tire deflection.
Figure 5. Frequency responses from the road profile zr1 to each output for full-state feedback case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass; (d) Suspension stroke; (e) Tire deflection.
Energies 14 08231 g005aEnergies 14 08231 g005bEnergies 14 08231 g005c
Figure 6. Frequency responses from the road profile zr1 to each output for static output feedback case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass.
Figure 6. Frequency responses from the road profile zr1 to each output for static output feedback case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass.
Energies 14 08231 g006aEnergies 14 08231 g006b
Figure 7. Frequency responses from the road profile zr to each output for the quarter-car model. (a) Vertical acceleration of the sprung mass; (b) Suspension stroke; (c) Tire deflection.
Figure 7. Frequency responses from the road profile zr to each output for the quarter-car model. (a) Vertical acceleration of the sprung mass; (b) Suspension stroke; (c) Tire deflection.
Energies 14 08231 g007
Figure 8. Frequency responses from the road profile zr1 to each output for LQG case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass.
Figure 8. Frequency responses from the road profile zr1 to each output for LQG case. (a) Vertical acceleration of the sprung mass; (b) Roll angle of the sprung mass; (c) Pitch angle of the sprung mass.
Energies 14 08231 g008aEnergies 14 08231 g008b
Table 1. Parameters of the 7-DOF full-car model [29].
Table 1. Parameters of the 7-DOF full-car model [29].
ms1653.0 kgmu45.0 kg
Ix614.0 kg·m2Iy2765.0 kg·m2
lf1.402 mlr1.646 m
tf, tr0.8 mkt230,000 N/m
ks34,000 N/mbs3500 N·s/m
Table 2. Maximum allowable values in LQ cost function.
Table 2. Maximum allowable values in LQ cost function.
ξ10.2 m/s2ξ210 deg/s2ξ310 deg/s2
ξ43 degξ510 deg/sξ63 deg
ξ710 deg/sξ80.2 mξ90.2 m
ξ103000 Nη1, χ1ξ1η2, χ2ξ8
η3, χ3ξ9η4, χ4ξ10
Table 3. Covariance matrices in LQG for quarter-car and full-car models.
Table 3. Covariance matrices in LQG for quarter-car and full-car models.
Mq104·I4×4Nq10−4
Mf104·I14×14Nf10−4·I8×8
Table 4. Gain matrices of LQR and LQSq.
Table 4. Gain matrices of LQR and LQSq.
Kq [ 28,929.0 31,583.0 1538.4 3017.5 ]
Sq [ 30,152.0 32,070.0 1968.9 3102.7 ]
Table 5. Magnitudes of frequency responses for each controller (dB scale).
Table 5. Magnitudes of frequency responses for each controller (dB scale).
z ¨ c ϕθzs1zu1zu1zr1
Passive38.9−7.8−12.40.6−3.1
LQRf129.9−34.9−35.11.9−9.8
LQRf224.4−23.1−28.51.9−9.8
LQRfq23.7−20.6−27.32.8−8.0
LQSCfq22.0−21.0−28.22.8−8.2
Table 6. Gain matrices of LQSOF and LQSOFq.
Table 6. Gain matrices of LQSOF and LQSOFq.
KSOF [ 31,285.0 2972.2 ]
Vq [ 31,285.0 31,285.0 2972.2 2972.2 ]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Park, M.; Yim, S. Design of Static Output Feedback and Structured Controllers for Active Suspension with Quarter-Car Model. Energies 2021, 14, 8231. https://doi.org/10.3390/en14248231

AMA Style

Park M, Yim S. Design of Static Output Feedback and Structured Controllers for Active Suspension with Quarter-Car Model. Energies. 2021; 14(24):8231. https://doi.org/10.3390/en14248231

Chicago/Turabian Style

Park, Manbok, and Seongjin Yim. 2021. "Design of Static Output Feedback and Structured Controllers for Active Suspension with Quarter-Car Model" Energies 14, no. 24: 8231. https://doi.org/10.3390/en14248231

APA Style

Park, M., & Yim, S. (2021). Design of Static Output Feedback and Structured Controllers for Active Suspension with Quarter-Car Model. Energies, 14(24), 8231. https://doi.org/10.3390/en14248231

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