^{1}

^{1}

^{2}

^{1}

^{2}

A method is designed to optimize the weight matrix of the LQR controller by using the simulated annealing algorithm. This method utilizes the random searching characteristics of the algorithm to optimize the weight matrices with the target function of suspension performance indexes. This method improves the design efficiency and control performance of the LQR control, and solves the problem of the LQR controller when defining the weight matrices. And a simulation is provided for vehicle active chassis control. The result shows that the active suspension using LQR optimized by the genetic algorithm compared to the chassis controlled by the normal LQR and the passive one, shows better performance. Meanwhile, the problem of defining the weight matrices is greatly solved.

The suspension system is an important component of the vehicle, and its performance is a major factor that affects the vehicle ride comfort, operation, and stability. The traditional passive suspension is generally composed of the elastic component and damping components with the fixed parameters. And it is generally designed to meet a particular road, so it limits the vehicle performance. Therefore, industrial countries have already begun to study the active suspension and the semiactive suspension systems based on the active vibration control since the 1970s. In recent years, with the rapid development of the electronic technology, testing techniques, and system dynamics theories, the vibration isolation for the vehicle suspension system has changed from the traditional passive to the active vibration control. There is a lot of research in the suspension control areas in China. The application of the modern control theory in the vehicle suspension control also has a variety of algorithms, such as the PID, skyhook, optimal, and sliding-mode controllers [

The conventional fuzzy control lacks the control precision, and it is difficult to establish fuzzy rules. And the high complexity of the neural network, as a representative of the intelligent control technology, makes it impossible to analyze the various performance indexes accurately. Therefore, the system versatility cannot be satisfied.

Among the modern control theories, the optimal control is widely used, and its theoretical foundation is also perfect. As a kind of optimal control, the linear quadratic regulator (LQR) can provide the required optimal performance indexes, by determining the weighting matrix of the system state and control variables. And it can also give sufficient design space for designers. When applied to the vehicle active suspension design, the LQR control can put forward different objective functions for the different performance need.

That is to say, the LQR control can improve the vehicle performance by the comprehensive consideration of the various factors in the vehicle suspension system.

As to the LQR control, the choice of the weight matrix for the performance indexes has great influence on the control system performance. Usually, the weight matrix is determined manually according to the system’s physical process, so it needs full understanding and sufficient experience of the system. Therefore, such “optimal” results are completely subjective and artificial. If the weight matrix is improperly selected, the optimal solution of no practical significance. In addition, the choice of the weighted matrix relies too much on the designer’s experience. The designers must adjust the weighted matrix step by step according to the system output until a satisfactory response value of output is acquired. Thus, the obtaining of feedback coefficient for the optimal control cannot guarantee the system to achieve the optimal aims [

At present, some random search methods such as the Genetic Algorithm (GA) and the Simulated Annealing (SA) have already caused the interest of many scholars. And they are already applied in the active suspension control field successfully. The simulated annealing algorithm is such a kind of the random search methods, and it is an approximation algorithm suitable and effective for solutions to the optimization problems in the mass general combinatorial fields. Compared with the previous approximate algorithm, the SA has the advantages of simple description, flexibility in use, universality, high efficiency, little constraint by the initial conditions, and so forth.

Due to the problems of the LQR in determining the weight matrix, a vehicle active suspension is selected as the control object, and the SA is applied in the LQR controller design optimized by the objective function of the active suspension performance index, aiming at improving the design efficiency and property.

For research convenience, a simplified quarter-vehicle model is selected as the research object, as Figure

The quarter-vehicle body model of 2 DOF.

In Figure

The state variables and the output variables are selected, respectively, as follows:

The control input is the active force

Then, the system state equation can be achieved:

The main performance for automobile suspension design includes three indexes: the body vertical vibration acceleration which represents the car ride, the suspension dynamic travel which represents the body posture and the suspension structure, and the dynamic load of the tyre which represents the tire grounding characteristic.

Therefore, the target performance index

In the equation, the symbols

Then, formula (

The symbols

The feedback matrix

According to the feedback variables

The

According to the above formulas, the optimal control design depends on the choice of weighting coefficients totally. Previous designs often rely on the designer’s own experience through repeated trials. Though corresponding optimal results can be got, a large subjectivity is met.

Therefore, the simulated annealing algorithm is used to optimize the weights coefficient matrix in order to reduce the controller design time and avoid the subjectivity at the same time.

The simulated annealing (SA) algorithm is put forward by Metropolis at the earliest. His starting point is based on the similarity of the solid material annealing process to the general combinatorial optimization problem. The SA algorithm is a general optimization algorithm with its annealing process composed of three processes: the heating process, the isothermal process, and the cooling process [

The optimum controller design based on simulated annealing algorithm is a combination of the SA and LQR control method. The SA is used to optimize and solve the weight matrix, and the LQR is used to calculate the system performance index.

The weights coefficient is taken as the optimized variable of the SA:

Because there is difference in the order of magnitude and the units, the optimization aim of the objective function is taken as follows:

In the equations, the symbols are noted as follows:

: the root mean square value of relevant data;

The optimization process of the SA is shown in Figure

Creation of the initial feasible solution:

Then, new solutions are created through the neighborhood function according to the current solution

Make sure that the criteria are met at current temperature. If ok, then the process continues or returns to step

Then, the cooling process is operated and

Make sure the algorithm meets the principle of termination. If it is ok, then the process is terminated or returns to step

The optimization process for the LQR controller by the SA.

Under the environment of the Matlab/Simulink, simulations are performed as to the active suspensions controlled by the LQR controller optimized by the SA and other controlling manners. The B-grade road in the standard file-ISO/TC108/SC2N67 is selected according to the previous model. A random road input is generated at the speed of 20 km/h, and ^{3}/cycle. And

Vertical acceleration of the active and passive suspension.

The LQR weights that the matrix gets from the simulation are

In order to validate the superiority of the SALQR in the suspension control, the performance results of the suspension in Simulink simulation are compared and are shown in Table

Performance contrast of the suspension of different control measures.

Performance indexes | Unit | The mean square root value | ||
---|---|---|---|---|

Passive suspension | Conventional LQR | SALQR | ||

Vertical acceleration of suspension | m/s^{2} |
2.8358 | 1.6342 | 0.3326 |

Dynamic travel of suspension | mm | 40.7465 | 47.7899 | 35.4843 |

Dynamic deflection of the tyre | mm | 11.8446 | 6.5684 | 24.7046 |

It can be seen from Figures

The SALQR controller significantly reduces the vertical acceleration of the suspension. In this regard, its effect is very well.

In regard of the suspension’s dynamic travel, the SALQR also gives good effect compared to the passive suspension and the conventional LQR suspension.

In the regard of the tire’s dynamic deflection, the effect of the SALQR is no good compared to other controlling manners. But the weakness is not so significant to affect the whole performance of the suspension.

Dynamic travel of the active and passive suspension.

Dynamic deflection of the active and passive suspension.

In view of the reliance on the designer’s experience of the conventional LQR controller, a controller combining the SA and LQR is designed, and it has the following advantages.

The dependence on the designer’s experience is reduced and the design convenience and efficiency are improved.

The controller performance is enhanced compared to the conventional LQR controller; therefore, the SA has a wide application prospect in the active control of automobile suspension.

The authors declare that there is no conflict of interests regarding the publication of this paper.

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (no. 51275060) and the Natural Science Foundation of Jiangsu Province (no. BK2011367).