# The Design of LQG Controller for Active Suspension Based on Analytic Hierarchy Process

Authors: Lingjiang Chai; Tao Sun
Journal: Mathematical Problems in Engineering
 (2010)

---

## Abstract

A full vehicle model with seven degrees of freedom based on active suspension control is established, and linear quadratic gaussian (LQG) is designed by applying optimal control theory. Especially, the methodology of Analytic Hierarchy Process (AHP) is used to make the selection of weighted coefficients of performance indexes, which can reduce ineffective job in contrast with experience method. From the simulation results, it is shown that ride quality of the vehicle with active suspension has been effectively improved in comparison with the vehicle of passive suspension by the methodology of AHP applying to the selection of the weights.


---

## Body

## 1. Introduction

Active suspension plays a major role for improving ride comfort of ground vehicle, which depends on its control algorithm design to a great extent, due to its effective isolation from road surface unevenness. Therefore, many experts have done researching for suspension control algorithms, such as optimal control theory, fuzzy control and neural net technology, and robust control [1–3]. Among them, optimal control theory is widely used, in which the linear quadratic gaussian control can ascertain weighted matrix of state variable and control variable by giving the optimal performance indexes. So the designer is offered an approach associated with weights selection of overall consideration of various factors of suspension system, to improve ride quality of ground vehicle.The approach of selecting weighted coefficients is conventionally designed by experience, and results in a time-consuming and inefficient job. In [4, 5], it does not provide a method with the selection of weighting coefficients. In [6], AHP applying to the selection of the weights is only used in quarter vehicle model with 2 degrees of freedom (DOF), which makes the limitations to design and analysis of vehicle performance. This paper makes use of AHP to optimize the weights of performance indexes in the design of active suspension LQG controller based on full vehicle model with 7 DOF. The simulation results illustrate that the body acceleration of active suspension has been effectively reduced, and ride quality of the vehicle  has been improved by AHP applying to the selection of the weights.


## 2. Modeling of Active Suspension with 7 Degrees of Freedom

Building a dynamic model of vehicle suspension is the basis of system design and performance analysis. Quarter vehicle model with 2 DOF and half vehicle model with 4 DOF have limitations to design and analyze vehicle performance. In order to reflect vertical vibration of the vehicle, quarter vehicle model is often used to conceptual design of suspension and research of control theory. Half vehicle model can only reflect vertical and pitch vibration of the vehicle. Therefore, the completed vehicle model with 7 DOF in this paper is used to grasp comprehensive ride quality including body vertical jump, pitch and roll vibration of the vehicle. The model is shown in Figure1.Figure 1
Model of completed vehicle with seven degrees freedom of suspension.The completed vehicle suspension itself is shown to consist of four springsk1,k2,k3,andk4, four dampers c1,c2,c3,andc4 and four active forces actuator u1, u2, u3, and u4. The unsprung mass m1, m2, m3, and m4 represent the equivalent mass due to the axle and tire. The vertical stiffness of the tires is represented by the springs k5, k6, k7, and k8. The variables z1, z2, z3, and z4 represent the vertical displacements from static equilibrium of the unsprung mass. The variables z5, z6, z7, and z8 represent the vertical displacements from static equilibrium of the sprung mass. The variables z01, z02, z03, and z04 represent the vertical displacements from static equilibrium of the road. zb is vertical displacement of body mass. mb is body mass. θ which encloses x-axis rotation is roll angle. ϕ which encloses y-axis rotation is pitch angle. The displacement from center of mass to front axle is a. The displacement from center of mass to rear axle is b. The displacement from center of mass to the center of left tire is ll. The displacement from center of mass to the center of right tire is lr.The mathematical model is set up as follows:(2.1)mbz̈b=c1(ż1-ż5)+k1(z1-z5)+u1+c2(ż2-ż6)+k2(z2-z6)+u2+c3(ż3-ż7)+k3(z3-z7)+u3+c4(ż4-ż8)+k4(z4-z8)+u4,IPϕ̈=[c3(ż3-ż7)+k3(z3-z7)+u3+c4(ż4-ż8)+k4(z4-z8)+u4]b-[c1(ż1-ż5)+k1(z1-z5)+u1+c2(ż2-ż6)+k2(z2-z6)+u2]a,Irθ̈=[c2(ż2-ż6)+k2(z2-z6)+u2+c4(ż4-ż8)+k4(z4-z8)+u4]lr-[c1(ż1-ż5)+k1(z1-z5)+u1+c3(ż3-ż7)+k3(z3-z7)+u3]ll,m1z̈1=k5(z01-z1)+k1(z5-z1)+c1(ż5-ż1)-u1,m2z̈2=k6(z02-z2)+k2(z6-z2)+c2(ż6-ż2)-u2,m3z̈3=k7(z03-z3)+k3(z7-z3)+c3(ż7-ż3)-u3,m4z̈4=k8(z04-z4)+k4(z8-z4)+c4(ż8-ż4)-u4.Whenθ, ϕ is very small, the displacement equations of four end point with sprung mass are(2.2)z5=zb+llθ-aϕ,z6=zb-lrθ-aϕ,z7=zb+llθ+bϕ,z8=zb-lrθ+bϕ.The road input model is the basis of automobile dynamic response and control. In general, the mathematical models of the road disturbances in time domain include random road disturbances model and discrete event road disturbances model. This paper will use them as road excitation input. There are four road inputs corresponding to left front tire, right front tire, left rear tire and right rear tire in this paper. So tread correlated of left and right tires and wheelbase lag of front and rear tires are considered. Based on these, random road input model [7] is built.Firstly, the tread correlated equation of left and right tires is set up as follows:(2.3)[η̇t1η̇t2]=Aη1[ηt1ηt2]+Bη1W1,W2=Cη1[ηt1ηt2]+Dη1W1.Secondly, the wheelbase lag equation of front and rear tires is set up and second-order Pade approximation algorithm is used to express time delay. So the model is(2.4)[η̇t3η̇t4]=[-6t-12t210][ηt3ηt4]+[12t20]W1,W3=[-t0][ηt3ηt4]+[0]W1.Lastly, a filtered white noise model of the four wheels random road inputs is(2.5)ż0=Fwz0+B0ηηt+B0W,η̇t=Aηηt+BηW,z=[z01z02z03z04]T, z01, z02, z03, z04 represent road inputs of left front, right front, left rear, right rear tire, respectively. ηt is the state transition vector, ηt=[ηt1ηt2ηt3ηt4]T. W is the white noise, W=[W1W2W3W4]T. t is time delay.Besides random road profile, the model [8] of discrete event road disturbances is also set up as follows: (2.6)z={-h2cos[2πvtxl]+h2,0≤t≤xlv,0,t>xlv.
The length of arc road is xl. The height from road surface to convex closure is h. v is velocity of vehicle running.


## 3. The Optimal Controller Design of Active Suspension

### 3.1. The State Equation of the System

For the completed vehicle model, the roll, pitch, jump vibration, suspension and tire deflection are very important performance items and should be chosen as output variables. So the state variablesXT, and output variables YT are shown as below:
(3.1)XT=[zbżbθθ̇ϕϕ̇z1ż1z2ż2z3ż3z4ż4],YT=[z̈bθ̈ϕ̈z5-z1z6-z2z7-z3z8-z4z1-z01z2-z02z3-z03z4-z04].The state space model of the completed vehicle system can be written as(3.2)Ẋ=AX+Bu+Fz,Y=CX+Du.In the above description, the state vector is denoted asX. u is the control vector. u=[u1u2u3u4]T. A is the system matrix. B is control matrix. F is disturbance matrix. C is the output matrix. D is transfer matrix.


### 3.2. The LQG Controller Design

The aim of optimal control of active suspension is to improve ride comfort and handling stability. The body acceleration, roll angle acceleration, pitch angle acceleration and tire deflection, therefore, should be depressed as much as possible. To avoid suspension bottoming, suspension working space is also limited. Considering the completed model is complex and has many trade-off performance items, the performance function is defined as follows:(3.3)J=12∫0∞[q1z̈b2+q2θ̈2+q3ϕ̈2+q4(z5-z1)2+q5(z6-z2)2+q6(z7-z3)2+q7(z8-z4)2+q8(z1-z01)2+q9(z2-z02)2+q10(z3-z03)2+q11(z4-z04)2+r1u12+r2u22+r3u32+r4u42]dt=12∫0∞(YTQY+uTRu)dt,q1 is the weighted coefficient of body vertical acceleration. q2 is the weighted coefficient of roll angle acceleration. q3 is the weighted coefficient of pitch angle acceleration. q4, q5, q6, and q7 are the weighted coefficients of suspension deflection. q8, q9, q10, and q11 are the weighted coefficients of tire deflection. r1, r2, r3, and r4 are the weighted coefficients of control vector. Q is the weight matrix of the state. R is the weight matrix of the control.The formula (3.2) is substituted into formula (3.3), and the performance function is then written as(3.4)J=12∫0∞[XTQ̂X+2XTNu+uTR̂u]dt.There areQ̂=CTQC, N=CTQD,  and R̂=R+DTQD in formula (3.4).The solution to the optimal control problem that minimizes this given performance index is a state feedback lawu=-KX where the feedback gain K is determined by solving the following Riccati equation:(3.5)AP+PAT+Q-PBR-1BTP=0,K=R-1(NT+BTP).
So the state equation of closed-loop is(3.6)Ẋ=(A-BK)X+Fz.


### 3.3. The Selection of the Weights Based on AHP

AHP is a multicriteria decision aiding method based on a solid axiomatic foundation. It is a systematic procedure for dealing with complex decision making problems in which many competing alternatives exist. The procedure can be described as the following three steps [9].(1) Making Pairwise Comparison
The evaluation matricesH=(hij)n×n are built up through pairwise comparing each decision factor under the topmost goal
(3.7)H=(hij)n×n=[1h12⋯⋯h1n1h121h23⋯h2n⋯1h23⋯⋯⋯⋮⋮⋮1h1n1h2n⋯⋯1].n is the number of factors.
The fundamental 1 to 9 scale can be used to rank the judgments as shown in Table1.
According toH matrices, the subject weights [6] are calculated as follows(1)
calculatingH multiplying vector of every row
(3.8)M=[M1M2⋯Mn]T,Mi=∏j=1nhij,(i,j=1,2,…,n),(2)
calculatingMn(3.9)E̅=[E̅1E̅2⋯E̅n]T,E̅i=Min,(i=1,2,…,n),(3)
calculating positive vector ofE̅(3.10)E=E̅∑i=1nE̅i(i=1,2,…,n),E is the subjective weighted coefficient of each evaluating index.Table 1
A fundamental scale of 1 to 9.
i/jEquallyModeratelyStronglyVeryExtremelyhij135792, 4, 6, and 8 indicate the medium value of above pairwise comparison.(2) Calculating Maximum Eigenvalue of Matrix H
(3.11)λmax=∑i=1n(HE)inEi(i=1,2,…,n).(3) Checking for Consistency
If every element in matrixH satisfies the equation hij=1/hji and hij=hik·hkj, the matrix H is the consistency matrix. The evaluation matrices are often not perfectly consistent due to people's random judgments. These judgment errors can be detected by a consistency ratio CR(3.12)CR=λmax-nRI(n-1)RI is random index. When n is 11, RI equals 1.51. The value CR calculated is less than 1, the result is right. If CR is more than 1, the evaluation matrix needs to be revised on the basis of [10].The formula fourteen has many weights which are calculated by AHP.Firstly, according to the performance data of a passive antitype suspension under the given working condition, objective weights of evaluating indexes can be determined. Passive suspension to the specific demand for ride comfort under a given operation condition hasσ. The objective weight of body acceleration is supposed to 1. According to [6], the equation of calculating other objective weights is(3.13)σa2×1=σ2×β.Secondly, it needs to calculate subjective weighted proportion coefficient of evaluating indexes. The subjective weighted proportion coefficient of body vertical acceleration is supposed to be 1. The subjective weighted proportion coefficients of suspension and tire deflection are(3.14)E1=Eiri(i=1,2,…,n).Lastly, the general weights of evaluating indexes of related to ride comfort can be obtained based on the objective weights and subjective weighted proportion coefficients(3.15)q=β×γ.The weights of the active suspension are derived as follows:(3.16)q1=1,q2≈1,q3≈1,q4=1567,q5=1538q6=764,q7=305,q8=4279,q9=22821,q10=2867,q11=2423.


## 3.1. The State Equation of the System

For the completed vehicle model, the roll, pitch, jump vibration, suspension and tire deflection are very important performance items and should be chosen as output variables. So the state variablesXT, and output variables YT are shown as below:
(3.1)XT=[zbżbθθ̇ϕϕ̇z1ż1z2ż2z3ż3z4ż4],YT=[z̈bθ̈ϕ̈z5-z1z6-z2z7-z3z8-z4z1-z01z2-z02z3-z03z4-z04].The state space model of the completed vehicle system can be written as(3.2)Ẋ=AX+Bu+Fz,Y=CX+Du.In the above description, the state vector is denoted asX. u is the control vector. u=[u1u2u3u4]T. A is the system matrix. B is control matrix. F is disturbance matrix. C is the output matrix. D is transfer matrix.


## 3.2. The LQG Controller Design

The aim of optimal control of active suspension is to improve ride comfort and handling stability. The body acceleration, roll angle acceleration, pitch angle acceleration and tire deflection, therefore, should be depressed as much as possible. To avoid suspension bottoming, suspension working space is also limited. Considering the completed model is complex and has many trade-off performance items, the performance function is defined as follows:(3.3)J=12∫0∞[q1z̈b2+q2θ̈2+q3ϕ̈2+q4(z5-z1)2+q5(z6-z2)2+q6(z7-z3)2+q7(z8-z4)2+q8(z1-z01)2+q9(z2-z02)2+q10(z3-z03)2+q11(z4-z04)2+r1u12+r2u22+r3u32+r4u42]dt=12∫0∞(YTQY+uTRu)dt,q1 is the weighted coefficient of body vertical acceleration. q2 is the weighted coefficient of roll angle acceleration. q3 is the weighted coefficient of pitch angle acceleration. q4, q5, q6, and q7 are the weighted coefficients of suspension deflection. q8, q9, q10, and q11 are the weighted coefficients of tire deflection. r1, r2, r3, and r4 are the weighted coefficients of control vector. Q is the weight matrix of the state. R is the weight matrix of the control.The formula (3.2) is substituted into formula (3.3), and the performance function is then written as(3.4)J=12∫0∞[XTQ̂X+2XTNu+uTR̂u]dt.There areQ̂=CTQC, N=CTQD,  and R̂=R+DTQD in formula (3.4).The solution to the optimal control problem that minimizes this given performance index is a state feedback lawu=-KX where the feedback gain K is determined by solving the following Riccati equation:(3.5)AP+PAT+Q-PBR-1BTP=0,K=R-1(NT+BTP).
So the state equation of closed-loop is(3.6)Ẋ=(A-BK)X+Fz.


## 3.3. The Selection of the Weights Based on AHP

AHP is a multicriteria decision aiding method based on a solid axiomatic foundation. It is a systematic procedure for dealing with complex decision making problems in which many competing alternatives exist. The procedure can be described as the following three steps [9].(1) Making Pairwise Comparison
The evaluation matricesH=(hij)n×n are built up through pairwise comparing each decision factor under the topmost goal
(3.7)H=(hij)n×n=[1h12⋯⋯h1n1h121h23⋯h2n⋯1h23⋯⋯⋯⋮⋮⋮1h1n1h2n⋯⋯1].n is the number of factors.
The fundamental 1 to 9 scale can be used to rank the judgments as shown in Table1.
According toH matrices, the subject weights [6] are calculated as follows(1)
calculatingH multiplying vector of every row
(3.8)M=[M1M2⋯Mn]T,Mi=∏j=1nhij,(i,j=1,2,…,n),(2)
calculatingMn(3.9)E̅=[E̅1E̅2⋯E̅n]T,E̅i=Min,(i=1,2,…,n),(3)
calculating positive vector ofE̅(3.10)E=E̅∑i=1nE̅i(i=1,2,…,n),E is the subjective weighted coefficient of each evaluating index.Table 1
A fundamental scale of 1 to 9.
i/jEquallyModeratelyStronglyVeryExtremelyhij135792, 4, 6, and 8 indicate the medium value of above pairwise comparison.(2) Calculating Maximum Eigenvalue of Matrix H
(3.11)λmax=∑i=1n(HE)inEi(i=1,2,…,n).(3) Checking for Consistency
If every element in matrixH satisfies the equation hij=1/hji and hij=hik·hkj, the matrix H is the consistency matrix. The evaluation matrices are often not perfectly consistent due to people's random judgments. These judgment errors can be detected by a consistency ratio CR(3.12)CR=λmax-nRI(n-1)RI is random index. When n is 11, RI equals 1.51. The value CR calculated is less than 1, the result is right. If CR is more than 1, the evaluation matrix needs to be revised on the basis of [10].The formula fourteen has many weights which are calculated by AHP.Firstly, according to the performance data of a passive antitype suspension under the given working condition, objective weights of evaluating indexes can be determined. Passive suspension to the specific demand for ride comfort under a given operation condition hasσ. The objective weight of body acceleration is supposed to 1. According to [6], the equation of calculating other objective weights is(3.13)σa2×1=σ2×β.Secondly, it needs to calculate subjective weighted proportion coefficient of evaluating indexes. The subjective weighted proportion coefficient of body vertical acceleration is supposed to be 1. The subjective weighted proportion coefficients of suspension and tire deflection are(3.14)E1=Eiri(i=1,2,…,n).Lastly, the general weights of evaluating indexes of related to ride comfort can be obtained based on the objective weights and subjective weighted proportion coefficients(3.15)q=β×γ.The weights of the active suspension are derived as follows:(3.16)q1=1,q2≈1,q3≈1,q4=1567,q5=1538q6=764,q7=305,q8=4279,q9=22821,q10=2867,q11=2423.


## 4. Simulation Results and Analysis

The parameters [11] of a car are used in this paper as shown in Table 2. The velocity of the car is 20 m/s. The simulation time is 20 seconds in all. The  road profile is a mixed random road excitation (smooth and rough) with the road roughness G0=5×10-6 at the first 10 seconds, and 16×10-6 at the rest. The system simulation model of active suspension and passive suspension is built and simulated as in Table 2.Table 2
Parameter of a car.
ParameterValueParameterValueParameterValuea/m1.3k1, k2, k3, k4/N/m20000Ip/kg·m22440b/m1.5k5, k6, k7, k8/N/m200000Ir/kg·m2380ll, lr/m0.75c1, c2,c3,c4/N·s/m1000mb/kg1380m1, m2,m3, m4/kg40


### 4.1. Root Mean Square (RMS) of the Active Suspension and Passive Suspension

From Table3, it is shown that the active suspension with LQG controller reduces body acceleration and roll angle acceleration remarkably on two grade roads in contrast with the passive one. The RMS value of body acceleration of the active suspension on both grade roads reduces 13% and 8%, respectively. The RMS value of roll angle acceleration in two grade roads reduces 14% and 15%, respectively. But the pitch angle acceleration is not noticeable to depress. The simulation results also show that the RMS value of right front axle of the active suspension working space reduces 23% and 12%, respectively, on the two grade roads. But the RMS value of left rear axle working space reduces 10% and 2%, respectively. It is obviously that control effect for left rear axle working space is even worse than right front axle working space. It can be seen that left rear wheel dynamic deflection of the passive suspension is better than that of the active suspension in Table 3, which shows that handling stability and ride comfort are incompatible to some extent.Table 3
Comparison of RMS for the active and passive suspensions.
IndexUnitRMSG0=5×10-6G0=16×10-6activepassiveactivepassivez̈bm/s21.1721.3421.69731.8356θ̈rad/s21.30021.50832.07462.4298ϕ̈rad/s20.862720.880791.34691.3609z6-z2mm12.9716.89820.2823.031z2-z02mm2.35752.75893.70354.1664z7-z3mm20.3522.54927.13627.752z3-z03mm6.85686.786110.91410.714


### 4.2. Time Response of the Two Different Systems in Mixed Road Profile

From Figures2 and 3, it can be seen that the active suspension with LQG controller remarkably reduces body acceleration and roll angle acceleration on both grade roads in contrast with the passive one. For Figures 4 and 5, it can be seen that the axles working space are controlled at the scope of the design requirements (±100 mm). According to Gaussian characteristics of road profile distribution, the RMS of the largest suspension deflection is 33 mm. It is less than theoretical value and the design is reasonable. Compared to right front axle working space, the control effect of left rear axle working space is relatively worse. In addition, wheel dynamic deflection related to handling stability impacts the effect of tire and road adhesion. For Figures 6 and 7, the control has no effect on the left rear wheel dynamic deflection on two grade roads with comparison to the passive suspension, which results from the interaction and coupling of the indexes.Figure 2
Body acceleration.Figure 3
Roll angle acceleration.Figure 4
Right front axle working space.Figure 5
Left rear axle working space.Figure 6
Right front wheel dynamic deflection.Figure 7
Left rear wheel dynamic deflection.


### 4.3. Time Response of the Two Different Systems in Random Road Excitation with Variations of Suspension Stiffness and Damping Parameters

Figures8 and 9 are right front axle working space with suspension stiffness parameters decreasing 10 percent and increasing 10 percent, respectively. Figures 10 and 11 are right front axle working space with suspension damping coefficients decreasing 10 percent and increasing 10 percent, respectively. From Figures 8–11, it can be seen that active suspension control has shown a good effect at mixed road excitation in comparison with the passive suspension. It illustrates that LQG controller suspension system based on AHP has a good stability of performance when suspension stiffness and damping coefficients vary.Figure 8
Right front axle working space (stiffness-18000).Figure 9
Right front axle working space (stiffness-22000).Figure 10
Right front axle working space (damping-900).Figure 11
Right front axle working space (damping-1100).


### 4.4. Active Suspension Control Force at Mixed Random Road Excitation

It can be seen that, from Figures12 and 13, suspension control force is within the control of actuator(±1000 N). Under the effect of suspension control force, active suspension can effectively reduce to vibration transmission from uneven roads to the body of the vehicle. Therefore, active suspension can availably depress values of body acceleration, roll angle acceleration and so on in contrast with passive suspension. So the vehicle of active suspension can effectively improve ride comfort.Figure 12
Right front suspension control force.Figure 13
Left rear suspension control force.


### 4.5. Power Spectral Density of the Two Different Systems in Conditions of the Same Road Roughness Coefficient (16×10-6)

From Figures14, 15, and 16, it can be illustrated that vibration-isolation effect of active suspension is satisfactory at low frequency related to ride comfort of ground vehicle to a large extent, including body bounce vibration, roll and pitch mode and human sensitivity band. Meanwhile, the amplitudes near the tire resonance frequency associated with handling and stability of ground vehicle have no noticeable changes.Figure 14
Power spectral density of BA.Figure 15
Power spectral density of RA.Figure 16
Power spectral density of PA.


### 4.6. Time Response of the Two Different Systems in Discrete Event Road Excitation

It can be seen that, from Figures17, 18, 19, and 20, the body acceleration, pitch angle acceleration, right front axle working space and left rear wheel dynamic deflection of the active suspension have improved noticeably, the amplitudes of performance indexes become small and the oscillation time are very short in contrast with the passive suspension.Figure 17
Body acceleration.Figure 18
Pitch angle acceleration.Figure 19
Right front axle working space.Figure 20
Left rear wheel dynamic deflection.


## 4.1. Root Mean Square (RMS) of the Active Suspension and Passive Suspension

From Table3, it is shown that the active suspension with LQG controller reduces body acceleration and roll angle acceleration remarkably on two grade roads in contrast with the passive one. The RMS value of body acceleration of the active suspension on both grade roads reduces 13% and 8%, respectively. The RMS value of roll angle acceleration in two grade roads reduces 14% and 15%, respectively. But the pitch angle acceleration is not noticeable to depress. The simulation results also show that the RMS value of right front axle of the active suspension working space reduces 23% and 12%, respectively, on the two grade roads. But the RMS value of left rear axle working space reduces 10% and 2%, respectively. It is obviously that control effect for left rear axle working space is even worse than right front axle working space. It can be seen that left rear wheel dynamic deflection of the passive suspension is better than that of the active suspension in Table 3, which shows that handling stability and ride comfort are incompatible to some extent.Table 3
Comparison of RMS for the active and passive suspensions.
IndexUnitRMSG0=5×10-6G0=16×10-6activepassiveactivepassivez̈bm/s21.1721.3421.69731.8356θ̈rad/s21.30021.50832.07462.4298ϕ̈rad/s20.862720.880791.34691.3609z6-z2mm12.9716.89820.2823.031z2-z02mm2.35752.75893.70354.1664z7-z3mm20.3522.54927.13627.752z3-z03mm6.85686.786110.91410.714


## 4.2. Time Response of the Two Different Systems in Mixed Road Profile

From Figures2 and 3, it can be seen that the active suspension with LQG controller remarkably reduces body acceleration and roll angle acceleration on both grade roads in contrast with the passive one. For Figures 4 and 5, it can be seen that the axles working space are controlled at the scope of the design requirements (±100 mm). According to Gaussian characteristics of road profile distribution, the RMS of the largest suspension deflection is 33 mm. It is less than theoretical value and the design is reasonable. Compared to right front axle working space, the control effect of left rear axle working space is relatively worse. In addition, wheel dynamic deflection related to handling stability impacts the effect of tire and road adhesion. For Figures 6 and 7, the control has no effect on the left rear wheel dynamic deflection on two grade roads with comparison to the passive suspension, which results from the interaction and coupling of the indexes.Figure 2
Body acceleration.Figure 3
Roll angle acceleration.Figure 4
Right front axle working space.Figure 5
Left rear axle working space.Figure 6
Right front wheel dynamic deflection.Figure 7
Left rear wheel dynamic deflection.


## 4.3. Time Response of the Two Different Systems in Random Road Excitation with Variations of Suspension Stiffness and Damping Parameters

Figures8 and 9 are right front axle working space with suspension stiffness parameters decreasing 10 percent and increasing 10 percent, respectively. Figures 10 and 11 are right front axle working space with suspension damping coefficients decreasing 10 percent and increasing 10 percent, respectively. From Figures 8–11, it can be seen that active suspension control has shown a good effect at mixed road excitation in comparison with the passive suspension. It illustrates that LQG controller suspension system based on AHP has a good stability of performance when suspension stiffness and damping coefficients vary.Figure 8
Right front axle working space (stiffness-18000).Figure 9
Right front axle working space (stiffness-22000).Figure 10
Right front axle working space (damping-900).Figure 11
Right front axle working space (damping-1100).


## 4.4. Active Suspension Control Force at Mixed Random Road Excitation

It can be seen that, from Figures12 and 13, suspension control force is within the control of actuator(±1000 N). Under the effect of suspension control force, active suspension can effectively reduce to vibration transmission from uneven roads to the body of the vehicle. Therefore, active suspension can availably depress values of body acceleration, roll angle acceleration and so on in contrast with passive suspension. So the vehicle of active suspension can effectively improve ride comfort.Figure 12
Right front suspension control force.Figure 13
Left rear suspension control force.


## 4.5. Power Spectral Density of the Two Different Systems in Conditions of the Same Road Roughness Coefficient (16×10-6)

From Figures14, 15, and 16, it can be illustrated that vibration-isolation effect of active suspension is satisfactory at low frequency related to ride comfort of ground vehicle to a large extent, including body bounce vibration, roll and pitch mode and human sensitivity band. Meanwhile, the amplitudes near the tire resonance frequency associated with handling and stability of ground vehicle have no noticeable changes.Figure 14
Power spectral density of BA.Figure 15
Power spectral density of RA.Figure 16
Power spectral density of PA.


## 4.6. Time Response of the Two Different Systems in Discrete Event Road Excitation

It can be seen that, from Figures17, 18, 19, and 20, the body acceleration, pitch angle acceleration, right front axle working space and left rear wheel dynamic deflection of the active suspension have improved noticeably, the amplitudes of performance indexes become small and the oscillation time are very short in contrast with the passive suspension.Figure 17
Body acceleration.Figure 18
Pitch angle acceleration.Figure 19
Right front axle working space.Figure 20
Left rear wheel dynamic deflection.


## 5. Conclusion

This paper applies the optimal control theory to design a Linear Quadratic Gaussian controller for the active suspension. Especially, the methodology of AHP is used to make the selection of weighted coefficients of performance indexes, which can reduce ineffective job in contrast with experience method. From the simulation results, it can be seen that ride quality of the vehicle with active suspension has been effectively improved in comparison with the vehicle of passive suspension. At the same time, it can also illustrate that active suspension control system based on AHP to select weighted coefficients of the LQG controller has a good stability when suspension damping coefficients and stiffness parameters vary. Furthermore, the future work will apply AHP to investigate stability of the integrated control system that consists of an active control of vehicle suspension and a steering system in cornering on curved by means of the optimal control theory.



---

*Source: 701951-2010-07-11.xml*