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.
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).
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.
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.
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.,
and
, (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.
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.
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).
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).
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.
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).
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).
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).
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.
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.
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.
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.