Design of Static Output Feedback Suspension Controllers for Ride Comfort Improvement and Motion Sickness Reduction

: This paper presents a method to design a static output feedback active suspension controller for ride comfort improvement and motion sickness reduction in a real vehicle system. Full-state feedback controller has shown good performance for active suspension control. However, it requires a lot of states to be measured, which is very di ﬃ cult in real vehicles. To avoid this problem, a static output feedback (SOF) controller is adopted in this paper. This controller requires only three sensor outputs, vertical velocity, roll and pitch rates, which are relatively easy to measure in real vehicles. Three types of SOF controller are proposed and optimized with linear quadratic optimal con-trol and the simulation optimization method. Two of these controllers have only three gains to be tuned, which are much smaller than those of full-state feedback. To validate the performance of the proposed SOF controllers, a simulation is carried out on a vehicle simulation package. From the results, the proposed SOF controllers are quite good at improving ride comfort and reducing motion sickness.


Introduction
To date, ride comfort and road handling have been regarded as two objectives to be satisfied in suspension design and control [1,2].The vertical acceleration, az, of a sprung mass (SPM) has been used as a measure for ride comfort.On the other hand, tire deflection has been used as a measure for road handling.Figure 1 shows ISO2631-1, which represents the sensitivity of human bodies to frequency ranges in the vertical and horizontal directions [3,4].As shown in Figure 1, the frequency range of az related to ride comfort, Wk: vertical, is 4.0~10.0Hz.For this reason, az within this range should be reduced by suspension control for ride comfort [4].
Since the early 2010s, autonomous driving has been intensively studied [5].Autonomous driving or advanced driving assistance systems (ADAS) can make passengers engage in non-driving-related visual tasks such as having a coffee, checking emails, and reading a book or display devices [6][7][8].This is a key advantage of autonomous driving.However, it can make motion sickness severe compared to conventional manual driving.For this reason, different from the studies on conventional suspension control for ride comfort, motion sickness should be taken into account when developing autonomous driving function.
For more than a decade, it has been known that motion sickness is caused by vertical vibration within the range of 0.1~0.2Hz, shown as Wf: Motion Sickness in Figure 1.This was obtained from the vibration on a ship.In the early 2000s, the studies carried out by Griffin showed that this range covers the range of 0.03~0.2Hz, shown as Wf: Griffin in Figure 1 [7,[9][10][11][12][13].For the last decade, studies investigating motion sickness on driving have shown that motion sickness is easily caused by combined az, pitch and roll rates (  and   ) in the 0.8~8.0Hz range, shown as Wf: DiZio in Figure 1  [14,15].This is overlapped with Wk: vertical, 4.0~10.0Hz.A comprehensive literature review on motion sickness can be found in reference [16].For this reason, az,   and   of an SPM over the range of 0.8~8.0Hz should be reduced in order to improve ride comfort and to reduce motion sickness.For the purpose of the previous study, those variables of an SPM were controlled with active suspension and the SOF control method [17].To date, several types of actuators have been used to control the suspension in vehicles such as active suspension and magneto-rheological (MR) damper.Among those actuators, active suspension control (ASC) has been intensively studied because it is quite an effective actuator for improving ride comfort and road handling.A literature survey on ASC was given in [1,2,[18][19][20].Recent developments over the last decade in the area of ASC are reviewed in [21].In this paper, active suspension is adopted as an actuator, which is needed to generate a vertical force in a suspension.
To improve ride comfort and to reduce motion sickness, the vertical, roll and pitch (VRP) motions ought to be controlled by a controller.For this purpose, a vehicle model describing those motions is needed.To date, quarter-, half-and full-car models have been selected for controller design [19].Among those models, a full-car model can describe the VRP motions of an SPM.For this reason, a full-car model is used as vehicle one in this paper.From this model, a linear state-space equation (SSE) is derived.In case there are nonlinear springs and dampers in a suspension, a Simulink model is built to describe their nonlinear behavior for simulation optimization (SLOM) [17].
To date, for ASC, various controller design methods such as linear quadratic optimal control (LQOC), adaptive and nonlinear control methods have been applied [22][23][24][25].Among those methods, the linear quadratic regulator (LQR) is the most commonly used method for ASC because it is systematic and easy to tune [20].LQR is a full-state feedback controller, which needs all system states to be measured or estimated for feedback.For instance, quarter-, half-and full-car models have 4, 8 and 14 state variables and one, two, and four control inputs, respectively.As a result, the gain matrices of LQR for quarter-, half-and full-car models have the dimensions of 1  4, 2  8, and 4  14, respectively.For the full-car model, LQR is too large to be implemented in a real vehicle.Moreover, it is very difficult to measure all system states in a real vehicle, and, consequently, LQR is very hard to implement on a real vehicle.To resolve this problem, two solutions have been adopted.The first is to use a state observer or Kalman filter to estimate state variables from sensor signals [22][23][24].However, the parameters of a model should be known exactly a priori for a state observer.Moreover, a state observer requires an extra design procedure besides LQR.Instead of a state observer, the second is to use a static output feedback (SOF) control with available sensor signals in real vehicles [17,[26][27][28][29]. Especially important, previous studies have only used two gains in the linear quadratic SOF (LQSOF) controller for Magnitude the full-car model [28,29].For this reason, SOF control is adopted as a controller structure in this paper.
In this paper, the vertical velocity (vz), roll and pitch rates (  and   ) of the SPM are selected as an available output for the SOF controller.In the full-car model, there are four control inputs.As a result, 4  3 gains are needed for the SOF controller.Besides the SOF controller with 12 gains in this paper, from the symmetry between roll and pitch motions of the SPM, two types of SOF controllers are proposed, which have only three gains [17].These SOF controllers are much easier to implement in a real vehicle because they have a much smaller number of gains.Moreover, it is also easier to optimize these controllers.
To design or optimize the SOF controller, linear quadratic optimal control (LQOC) and the simulation optimization method (SLOM) are adopted in this paper.LQOC uses the SSE and LQ cost function for controller design.In LQOC, an LQ cost function is derived from the state variables and minimized by a heuristic optimization method [28,29].SLOM uses the Simulink model with nonlinear springs and dampers and a cost function defined for improving ride comfort and reducing motion sickness.The Simulink model is built from the full-car model, and the cost function is evaluated from the results of the Simulink model and optimized by a heuristic optimization method (HOM) [17,30].As a HOM for optimization, the covariance matrix adaptation evolution strategy (CMA-ES) is selected in this paper [31].To check the control performance of the controllers designed by LQOC and SLOM, a simulation is performed on the vehicle simulation package, Car-Sim.By analyzing the simulation results, it is identified which controller is the best for ride comfort improvement and motion sickness reduction.
The objective of this paper is to design SOF active suspension controllers with LQOC and SLOM for the purpose of improving ride comfort and reducing motion sickness.The contributions of this paper can be summed up as follows: 1. Three types of structures of SOF controllers are presented.With vz,   and   of an SPM as an output, three SOF control structures are presented.This paper comprises four sections.In Section 2, a full-car model is presented and its SSE is derived.From the geometry of the SPM, three types of SOF control structures are proposed and designed with LQOC and SLOM.In Section 3, frequency response analysis and a simulation are conducted on CarSim.The conclusions are drawn in Section 4.

Design of Static Output Feedback Controllers
In this section, the SSE for the full-car models is derived, as presented in previous studies [17,28,29,32,33].With the SSE, LQR and LQSOF controllers are designed.For a vehicle with nonlinear spring and damper, a vehicle model is built with MATLAB/Simulink (version: 9.6.0.1472908,R2019a, The MathWorks Inc., Natick, MA, USA), and SLOM is applied to design SOF controllers.

Full-Car Model and State-Space Equation
Figure 2 shows a free-body diagram of the full-car model, which has four suspensions, S1, S2, S3 and S4.In the model, there are seven motions: three motions for the SPM and four motions of the unsprung mass (USPM).For the SPM, the VRP motions are described by three corresponding variables, i.e., the vertical displacement, zc, the roll angle, ϕ, and the pitch angle, θ, respectively.In the SPM, the vertical displacements at each corner are described by zs, which is determined by the geometry of the SPM and three variables, zc, ϕ and θ.For the UPSM, the vertical motions are described by the vertical displacement, zu.
As shown in Figure 2, each suspension Si has a spring of stiffness ksi, a damper of damping coefficient bsi and control input ui generated by an actuator.An active actuator is installed alongside a spring and a damper in a suspension.Those four actuators can generate the control inputs, u1, u2, u3 and u4, at each suspension.A controller calculates those control inputs from the variables describing the motions of the SPM and USPM.
In the full-car model, three motions of the SPM and four motions of the USPM are excited by four disturbances, i.e., the road profiles, zr1, zr2, zr3 and zr4, which are applied to the USPM.Generally, the disturbances exciting on the front USPMs do act on the rear ones in some delay.The delay between the front and rear suspensions depends on vehicle speed.In the suspension Si, the suspension force fi is calculated as in Equation ( 1).In Equation (1), ui is the control input generated by an actuator in Si.To derive a linear equation, the cubic terms in (1) are neglected.With those, three and four equations of motions for the SPM and USPM are derived as Equations ( 2) and (3), respectively.Assuming that there are no longitudinal and lateral motions for the SPM, i.e., vx = vy = 0, and that the cross moments of inertia, Ixy, Ixz and Iyz, can be neglected, Equation ( 2) is linearized as (4).Equation ( 4) is obtained as Equation ( 5).Hereafter, the matrix H plays an important role in this paper.
In Equation ( 1), the vertical displacements and velocities of the SPM at each corner, si z and si z  , are not a state variable.So, it is to be represented with the state variables, i.e., zc, ϕ and θ. Figure 3 shows how to calculate the vertical displacements of the SPM at four corners, i.e., , ,  and  in Figure 2. As shown in Figure 3, zs1, zs2, zs3 and zs4 are derived as in Equation ( 6) from the geometry of the SPM.In Equation ( 6), the nonlinear terms are approximated as sinθ  θ and sinϕ  ϕ assuming that ϕ and θ are small.With this approximation, Equation ( 6) is rearranged into Equation ( 7) [28,29,32,33].In Equation ( 7), the matrix H stands for the geometric relationship among four variables, zs1, zs2, zs3 and zs4, and the three state variables, zc, ϕ and θ.

Design of LQR
When designing LQR, an LQ cost function needs to be defined.The LQ cost function, J, with the state variables is given as in Equation (19).The weights ζi are determined by Bryson's rule, i.e., ζi = 1/ξi 2 , where ξi is the maximum allowable value on each term [34].For ride comfort improvement, ξ1, ξ2, and ξ3 ought to be set low.On the other hand, for road handling, ξ8 and ξ9 ought to be set low.Another objective of suspension control is to reduce motion sickness in this paper.According to recent studies, motion sickness is caused by combined vibrations along vertical and pitch directions in the 0.8~8.0Hz frequency range [14,15].For motion sickness reduction, the az and   of an SPM ought to be reduced.For the purpose, ξ1 and ξ7 ought to be set lower.
The LQ cost function J is rewritten into Equation (20).The matrices Q, R and N are derived from the matrices A, B2 and the weights in J. LQR is a full-state feedback controller, Equation (21), which minimizes J.It is easy to calculate the gain matrix of LQR, KLQR, from the solution of the Riccati equation, P, for A, B2, Q, R and N. The dimension of KLQR is 414, which is too large to implement on a real vehicle.Moreover, to implement KLQR on a real vehicle, 14 state variables should be measured or estimated by a state estimator.For this reason, it is very hard to implement KLQR on a real vehicle.To overcome this problem, a static output feedback (SOF) control is adopted in this paper.

Sensor Signal Processing
A SOF controller uses available signals measured with sensors in a real vehicle.In this paper, the vertical velocity, the pitch, and the roll rates of the SPM are assumed to be available for SOF control.Those signals are obtained from four accelerometers located at the four corners of the SPM, i.e., , ,  and , as shown in Figure 2 [35][36][37].The signals measured from those accelerometers are filtered sequentially by high-pass and low-pass filters in order to reject noise and DC blocking, respectively.Those filters are given in Equations ( 22) and ( 23), respectively [17,35].As pointed in the reference [35], Equation ( 22) plays a differentiator role below 0.1 Hz and an integrator role above 0.1 Hz, which excludes a possible DC offset.As a result, the vertical velocities at four corners of the SPM are calculated.With those velocities, the vertical velocity, roll and pitch rates of the SPM are obtained by Equation (11).This procedure of sensor signal processing is shown in Figure 4.  High-pass filter , 0.707, 0.1Hz 2

Design of SOF Controller with LQOC
Different from Figure 4, the vector of available outputs, y, ought to be derived from the state vector x and the matrices A, B1 and B2 for LQOC.The vector y is defined from the state vector x as in Equation (24).In this paper, three types of SOF controller are proposed.The first type of SOF controller is given in Equation (25).Since there are three sensor outputs and four control inputs, this type is natural when designing the SOF controller.As shown in Equation (25), there are twelve gains in KSOF.Let this controller be denoted as the SOF one.

 
The second type of SOF controller is given in Equation ( 26).This type is called the structured SOF or SSOF controller.This type is based on the fact that there are symmetries between the front/rear and left/right suspensions, as shown in Figure 3 [28].For this reason, the number of gains in KSSOF is three, which is much smaller than KSOF.The first column of KSSOF corresponds to the vertical control force.The second and third columns of KSSOF correspond to the roll and pitch moments, respectively.For this reason, the structure of KSSOF is identical to that of H T in Equation ( 7).
There are three motions in the SPM: VRP directions.If there are a vertical control force, Fzc, and roll and pitch moments, Mϕ and Mθ, used to control the VRP motions independently, these can be transformed into four control inputs in u with the matrix H, as shown in Equations ( 5) and (27).This is called Lotus modal control [38][39][40][41][42].In ( 27), H + is the pseudo-inverse of H. Fzc, Mϕ and Mθ can be determined by several methods such as LQR, LQSOF and sliding mode control (SMC) [42].In the previous study, two state variables are needed to calculate the control force and moments for each direction [42].For this reason, six gains are needed to generate the control inputs.In this paper, Fzc, Mϕ and Mθ are simply calculated by multiplying the constant gain matrix, KLSOF, by the vector of the vertical velocity, the roll and pitch rates, y [38].This is a simple derivative control which combines the Lotus modal decomposition with Karnopp's skyhook damper [38].Then, these three force/moments are converted into four control inputs at each corner with H + , as given in Equation (28).The number of gains in KLSOF is three, which is identical to KSSOF.This is the third type of SOF controller, called LSOF controller.
The gain matrices KSOF, KSSOF and KLSOF of the SOF controllers are optimized to minimize J.Those SOF controllers are called LQSOF, LQSSOF and LQLSOF controllers, respectively.To date, it has been shown that there are no analytical methods to optimize KSOF or KSSOF or KLSOF in terms of J [43].Moreover, it has also been shown that there are no analytic methods to find an initial gain which stabilizes a given system.If one of KSOF or KSSOF or KLSOF is selected, let this be Ks.With Ks, the closed-loop system Ac is calculated as in Equation (29).To optimize Ks, the optimization problem is formulated as in Equation (30) [17,28,29].In this paper, CMA-ES is applied as a HOM [31].

Design of SOF Controller with SLOM
The SOF controllers designed by LOQC, i.e., LQSOF, LQSSOF and LQLSOF controllers, are designed with the SSE of ( 18) and the LQ cost function of (19).Generally, the forces of springs and dampers are nonlinear with respect to the suspension stroke and suspension velocity in a real vehicle.If springs and dampers in the full-car model are nonlinear, the SSE cannot be derived.As a consequence, LQOC, (30), cannot be applied to design the LQSOF, LQSSOF and LQLSOF controllers.
To design three SOF controllers for nonlinear systems, a simulation optimization (SLOM) is applied in this paper [17,30].Simulation optimization is a method to find an optimum solution through simulation on a target system [44,45].For SLOM, the Simulink model is built from Equations ( 3) and ( 4), and the SOF controllers, Equations ( 25), ( 26) and ( 28).In the Simulink model, the springs and dampers are nonlinear.SLOM optimizes twelve elements in KSOF or three elements in KSSOF or KLSOF in terms of a certain cost function.
As mentioned earlier, the az,   and   of an SPM ought to be reduced for ride com- fort improvement and motion sickness reduction.For this purpose, the cost function of the SLOM, JSLOM, is defined as in Equation ( 31) by combining these values.For a given particular gain matrix, each term in JSLOM is obtained from a simulation carried out over the horizon T from t0 to t1.In Equation (31),  is the conversion constant from radian to degree, i.e., 57.2957795, and α and β are tuning parameters for   and   , respectively.Hereafter, α and β are set to 0.1.For optimization, CMA-ES is applied [29,31].When applying CMA-Es, no limits are imposed on each element in the gain matrices, KSOF, KSSOF and KLSOF.The optimization procedure of the SLOM with CMA-ES is given in Figure 5.Let the SOF controllers, SOF, SSOF and LSOF, designed by SLOM be called SOSOF, SOSSOF and SOLSOF, respectively.For SLOM, a road profile ought to be selected.In this paper, a twisted sine wave road (TSWR) is selected as a road profile, as given in Figure 6.TSWR has the period of 12.2 m and the amplitude of 0.05 m.In TWSR, the road profile on the right wheels is moved by the quarter period from that on the left ones.As a result, TSWR can simultaneously excite the VRP motions of the SPM.When using TSWR for SLOM, the vehicle speed is set to 15.0 m/s for the purpose of making severe conditions.

Simulation and Discussion
In this section, a simulation is carried out to evaluate the performance of six SOF controllers, LQSOF, LQSSOF, LQLSOF, SOSOF, SOSSOF and SOLSOF.With the simulation results, those controllers are compared to one another in terms of ride comfort and motion sickness.

Simulation Environment
Table 1 shows the parameters of the full-car model, which are references to an E-Class sedan given in CarSim [46].The maximum allowable values in J are given in Table 2.As shown in Table 2, for ride comfort improvement and motion sickness reduction, ξ1, ξ5 and ξ7, were set low.This set tries to reduce the az,   and   of an SPM.In this paper, the bandwidth of the actuator is assumed to be 20.0Hz.For SLOM, the force characteristic of nonlinear damper in the Simulink model are given in Figure 7.In CarSim, a linear spring was used, and its stiffness was set to 154,000.0N/m.Those data were references to those of the E-Class sedan model given in CarSim.As given in Section 2.5, SOSOF, SOSSOF and SOLSOF were designed with SLOM on MATLAB/Simulink.Six SOF controllers, LQSOF, LQSSOF, LQLSOF, SOSOF, SOSSOF and SOLSOF, were implemented and simulated on a co-simulation environment built with MATLAB/Simulink and CarSim.

Simulation with the SOF Controllers Designed with LQOC
With the parameters and the weights given in Tables 1 and 2, the gain matrices of LQSOF, LQSSOF and LQLSOF are given in Table 3.As shown in Table 3, those gain matrices are nearly identical to one another.This is a natural consequence because the structures of those controllers are nearly identical to one another and there are symmetries between the front/rear and left/right suspensions in view of the SPM.As a result, those controllers are expected to show equivalent performances to one another.The performance of those SOF controllers can be improved by setting the maximum allowable values given in Table 2 smaller except ξ10 on the control input.Generally, ξ10 on the control input has a significant effect on the control performance.However, this has little effect on the control performance because it is very small.Instead, ξ1, ξ5 and ξ7 should be set smaller.In the previous work, these values were not set small enough to give a better performance [17].In this paper, these values were set much smaller than those in the previous work.

KSOF
Let the LQLSOF controllers with the three sets of the maximum allowable values, (ξ1, ξ5, ξ7) of (0.2, 0.2, 0.2), (0.1, 0.1, 0.1) and (0.05, 0.05, 0.05) be denoted as Weight#1, Weight#2 and Weight#3, respectively.The other maximum allowable values are set as given in Table 2. Figure 8 shows the frequency responses of those three weight sets from the front left road profile, zr1, to az,   and   of an SPM.These responses are obtained with the SSE, Equation (18).As shown in Figure 8, the control performance is significantly improved as ξ1, ξ5 and ξ7 are set smaller.From the results, it is recommended to set ξ1, ξ5 and ξ7 as small as possible for ride comfort improvement and motion sickness reduction.Three SOF controllers, LQSOF, LQSSOF and LQLSOF, are simulated on a co-simulation environment with MATLAB/Simulink and CarSim.The controller gains given in Table 3 are used.The road profile used in CarSim is identical to that in Figure 6.The vehicle speed was set to 15 m/s and this is maintained as constant as possible by a speed controller.
Figure 9 shows the simulation results of those controllers designed with the weights given in Table 2. Table 4 summarizes the simulation results.As shown in Figure 9, those controllers show good performance in terms of ride comfort and motion sickness because az and   are significantly decreased by those controllers.These results can be expected from the controller gains given in Table 3 and the frequency responses given in Figure 8.The spring and damper forces obtained from the simulation results of four cases coincide with those in Figure 7, as given in Figure 9e.As shown in Table 4, there are rare differences among those controllers in terms of the az,   and   of the SPM.This is also expected from the controller gains given in Table 3.For this reason, only the control input of LQLSOF is given in Figure 9d.With the simulation results of the LQSOF, LQSSOF and LQLSOF controllers, the frequency responses are obtained as given in Figure 10.As shown in Figure 10, the az,   and   of an SPM were significantly reduced compared to the uncontrolled case below the frequency of 10 Hz.Especially important, the roll and pitch rates were reduced over the entire frequency range.These results mean that those three SOF controllers significantly improve ride comfort and motion sickness.Figure 11 shows the frequency responses of those controllers designed by SLOM on the TSWR.These plots were drawn with the SSE, as in Equation (18).As shown in Figure 11, SOSOF shows the best performance in terms of the roll and pitch rates.This is a natural consequence because SOSOF has twelve elements, which can cover a larger range of solution space than SSOF and LSOF. Figure 11b,c shows that there are a small differences between SOSOF and SOSSOF/SOLSOF.This is caused by the fact that there are no symmetries in SOSOF compared to SOSSOF and SOLSOF.In other words, the symmetries between the front/rear and left/right suspensions are not held for SOF.On the contrary, SSOF and LSOF have symmetries in their structures.As shown in Figure 11, the differences between SOSOF and SOSSOF/SOLSOF in terms of the pitch and roll rates are small, enough to be neglected because SOSSOF/SOLSOF shows quite good performance in controlling the pitch and roll motions of the SPM.For this reason, SOF is not recommended as a controller in terms of the number of gain elements.Figure 12 shows the simulation results of SOSOF, SOSSOF and SOLSOF on CarSim.Table 5 summarizes the simulation results of those controllers.As expected from Figure 11, those controllers show good performance in terms of az,   and   of an SPM.How- ever, as shown in Figure 12a, there is chattering in the vertical acceleration of SOLSOF, which cannot be applied to real vehicles.This is caused by large gains, which means large damping or derivative control.To remove this chattering, α and β in JSO should be tuned finely.As shown in Table 5, there are differences in the three measures for each controller.On the contrary, the controllers, LQSOF, LQSSOF and LQLSOF, show nearly identical and consistent results, as shown in Table 4.This is caused by the cost function JSO of SLOM.If α and β in JSO are tuned finely, then the SOSSOF and SOLSOF controllers can give better results than the LQSOF, LQSSOF and LQLSOF ones.With the simulation results of the SOSOF, SOSSOF and SOLSOF controllers, the frequency responses are obtained as given in Figure 13.As shown in Figure 13, SOLSOF shows the best performance, and SOSOF shows the worst performance.These results mean that SOLSOF is recommend as a suspension controller for ride comfort improvement and motion sickness reduction under the condition that it is tuned finely.In view of computation time, the SOF controllers designed by LQOC are superior to those designed by SLOM because SLOM requires a larger amount of computation time than LQOC.Moreover, by comparing the results given in Tables 4 and 5, SOF controllers designed by LQOC are recommended as a suspension controller for ride comfort improvement and motion sickness reduction.SOLSOF gains significantly improved the control performance in terms of ride comfort and motion sickness by reducing the vertical acceleration, and the roll and pitch rates of the sprung mass.From the experimental results, it was verified that the sensor signals for SOF controllers can be obtained from the accelerometers located at corners the SPM, and that the accelerometers on the wheel centers are not needed.In view of this controller design methodology, with the SOF controller structure, LQOC is preferred to SLOM because the former is more systematic than the latter and the latter is more case-dependent than the former.Moreover, the former can generate more consistent results than the latter regardless of vehicle type, parameters, and road profiles.For those reasons, the SOF controllers calculated with LQOC can be used as a starting point in experiments on real vehicles.This was confirmed by the experimental results.In further research, an experimental investigation will be carried out with the proposed controller on real vehicles, and the experimental results will be given within 6 months.

Figure 3 .
Figure 3. Calculation of vertical displacements at four corners of the SPM.(a) Front view; (b) side view.

Figure 4 .
Figure 4. Block diagram for velocity and angular rate calculation.

Figure 7 .
Figure 7. Damper force with respect to suspension stroke rate in the Simulink model.

Figure 8 .
Figure 8. Frequency responses of the SOF controllers designed with LQOC and three sets of weights.(a) Vertical accelerations; (b) roll rates; (c) pitch rates.

3. 3 .
Simulation with SOF Controllers Designed by SLOM Three SOF controllers, SOSOF, SOSSOF and SOLSOF, are designed by SLOM on MATLAB/Simulink.Different from the SOF controllers designed by LQOC, there are large Magnitude (dB) Magnitude (dB) Magnitude (dB) differences among the gains of those controllers.Those controllers were simulated on a co-simulation environment with MATLAB/Simulink and CarSim.

Figure 11 .Figure 6 .
Figure 11.Frequency responses of SOSOF, SOSSOF and SOLSOF designed with SLOM on TSWR.(a) Vertical accelerations; (b) roll rates; (c) pitch rates.Three SOF controllers designed by SLOM, i.e., SOSOF, SOSSOF and SOLSOF controllers, were simulated on CarSim.The road profile used in CarSim is identical to that of

T
, t0, t1 simulation horizon from the start time t0 to the end time t1 ui forces generated by an actuator at i-th suspension (N) vz =  c z vertical velocity of SPM (m/s) zc heave displacement at center of gravity of a sprung mass (m) zri road elevation acting on i-th tire (m) zsi vertical displacement of i-th corners of a sprung mass (m) zui vertical displacement of i-th wheel center (m) ξi maximum allowable value of weight in LQ cost function ζi weight in LQ cost function ϕ,   roll angle and rate of a sprung mass (rad, rad/s) θ,   pitch angle and rate of a sprung mass (rad, rad/s) ,

Table 2 .
Maximum allowable values in LQ cost function.

Table 4 .
Summary of the simulation results of LQSOF, LQSSOF and LQLSOF on CarSim.

Table 5 .
Summary of the simulation results of SOSOF, SOSSOF and SOLSOF on CarSim.