Published: 31 March 2014

Controlling and optimization of vehicle based on genetic algorithm method while encountering an emergency collision avoidance

Wei Wang1
You-qun Zhao2
Hui Yang3
Yue-qiao Chen4
Wen-tao Yang5
1, 2, 3, 4, 5College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, China
Corresponding Author:
You-qun Zhao
Views 37
Reads 11
Downloads 1187

Abstract

In consideration of the nonlinear tire lateral deviation, a nonlinear two degrees of freedom model is established in Matlab/Simulink, which can be used for emergency avoidance steering. The correctness of the model is verified by real vehicle test. Based on vehicle handling stability evaluation and single neuron self-adaptive Proportional-Integral-Differential control theory, a control algorithm to make vehicle move along the given single route with different speed is proposed. And the result of experiment is good. Next, taking advantage of multi-island genetic algorithm, an optimal control system is achieved to optimize the controlling parameters. The same experiment is carried out again with the optimized system. The results show that the optimized control algorithm is more accurate than the previous one. And it also has good robustness and self-adaption.

1. Introduction

With the continuous development of the automobile industry, the number of vehicle accidents grows rapidly, especially in traffic accidents involving pedestrians and cyclists. In some cases, vehicle accidents can be seen as a collision between vehicles and obstacles. To avoid such kind of accidents, the emergency avoidance problem is proposed. Today, people pay more and more attention to the problem of high-speed emergency avoidance [1-2].

A. B. Kurzhanski emphasized the importance of emergency avoidance problem (Kurzhanski, 2004), and some details were discussed at a groundbreaking speech [3].

The obstacle avoidance problem has been studied in the literatures. A brief review is presented in this section.

Richards and How (2002) introduced a method searching for the optimal trajectory of multiple aircraft collision avoidance [4]. It described the emergency avoidance problem as the linear programming meeting mixed integer constraints. Ryu et al. (2002) presented an emergency obstacle avoidance control strategy used in automated vehicles [5].

Hagenaars et al. (2004) pointed out that the objective function of continuous time and fixed terminal state were not suitable for the optimal control problem of emergency avoidance under the constraint of non convex state, and then put forward a new objective function including intermediate target state [6].

Mukai et al. (2007) transformed the problem of generating an optimal path without a collision between an automobile and obstacles as a mixed integer programming problem [7].

Darguzis A. et al. (2011) researched dynamic processes of a vehicle moving over step-shaped obstacles [8]. Ghaffari et al. (2011) selected lane change maneuver as the object behavior and proposed novel adaptive neuro-fuzzy inference models for the behavior [9]. Jin et al. (2011) redefined the equation of time-to-collision (TTC) using visual angle information by taking lateral separation into account [10]. So far, there are only a few researches which take into account path tracking optimization of vehicle while encountering an emergency collision avoidance.

The objective of this paper is to improve tracking accuracy. Aiming at the difficulty of identification of vehicle steering parameters which is caused by nonlinear dynamic characteristics, the single neuron self-adaptive Proportional-Integral-Differential (PID) controller is used and the vehicle emergency avoidance control algorithm is established. The optimization of control parameters is realized by using multi-island genetic algorithm. Results show that vehicle can track the given path well and have advantages in robustness and adaptability at different speed.

2. Nonlinear dynamics model of vehicle

2.1. Tire model

Considering the pure tire lateral characteristics and ignoring the impact of the tire roll characteristics, the lateral force semiempirical model which was proposed by Guo Konghui academician is used in this paper [11-13]:

1
F-y=1-exp-φ-E1φ2-E12+112φ3,
2
Fzn=FzFz_rated,
3
φ=Kytan(α)+Kyrsin(γ)μyFz,F-y=FyμyFz,Fz_rated=7110N,

where, Fz_rated is nominal load, E1 is curvature factor, Ky is lateral stiffness, Kyr is roll stiffness, μy is side friction coefficient, α is tire sideslip angle, and r is tire camber angle. E1, Ky, and Kyr, which will change along with the vertical load and the inclination angle of the tires, are defined as the structure parameters of the unified model:

4
E1=12+s12exp-Fzns22,
5
Ky=Fzs3+s4Fzn+s5Fzn2,
6
μy=s6+s7Fzn+s8Fzn2,

where s1~s7 are all identification parameters [14], and the values are s1=777.88, s2=0.39502, s3=0.06075, s4=–0.0368, s5=0.03755, s6=1.0227, s7=0.22898, s8=–0.20741.

Bringing the above parameters to the unified model and using the Matlab program, the relationship graphics between cornering force and sideslip angle of car’s front and rear wheels are shown in Fig. 1 and Fig. 2 respectively. Fig. 1 shows the relationship between cornering force and sideslip angle of front wheel. Fig. 2 shows the relationship between cornering force and sideslip angle of rear wheel.

From Fig. 1 and Fig. 2, when the sideslip angle is small, the relationship between cornering force and sideslip angle is linear. When the sideslip angle is large, the cornering force grows slowly. When the sideslip angle is 8 degree, the cornering force reaches saturation and almost stops growing.

Fig. 1The relationship between cornering force and sideslip angle of front wheel

The relationship between cornering force  and sideslip angle of front wheel

Fig. 2The relationship between cornering force and sideslip angle of front wheel

The relationship between cornering force  and sideslip angle of front wheel

2.2. Full vehicle model

In the model, the vehicle travels at a constant speed u, the inertia and damping of the steering system, the effect of suspension, the cornering properties change of left and right tires due to different vertical loads, and the influence of tire aligning torque are all ignored. As shown in Fig. 3, two generalized coordinates, yawing angle ψ and centroid sideslip angle β, are utilized to represent the vehicle’s motion state.

According to D’Alembert principle, the force equilibrium equations of the vehicle’s motion are as follows:

7
Muωr+β˙=-Py1+Py2,
8
IZω˙r=-aPy1+bPy2,

where, IZ is rotational inertia of the vehicle around Z-axis and Py1, Py2 are cornering force of front and rear wheels respectively.

Considering nonlinear tire cornering properties and combining with vehicle dynamics equations, a two degrees of freedom (DOF) nonlinear vehicle steering motion model is built in Matlab/Simulink, as shown in Fig. 4.

Fig. 32DOF vehicle steering model

2DOF vehicle steering model

Fig. 42DOF nonlinear vehicle Simulink model

2DOF nonlinear vehicle Simulink model

2.3. Experimental verification

A single lane-change real vehicle test of a car is done in this paper, and the test data is collected to compare with the results of simulation test by using the above nonlinear vehicle steering motion model. Fig. 5 shows the comparison results.

As shown in Fig. 5, the tendency of the two curves is consistent, and the amplitude of the two curves is slightly different, but the total difference is in the allowable range and there is a little difference of time between the two curves. It shows that the nonlinear vehicle steering motion model conforms to the reality, so the results from the model are correct.

Fig. 5Comparison figure between a single lane-change test results and simulation results

Comparison figure between a single lane-change test results and simulation results

3. Nonlinear vehicle steering control implementation

3.1. Basic theory

Single neuron self-adaptive PID control, which is based on quadratic form, is by means of quadratic performance index of optimal control in the single neuron learning algorithm index. Quadratic performance index is introduced into the adjustment of weighting coefficient. The weighting sum of squares for output error and controlled quantity is minimized to adjust weighting coefficient, in this way constrained control of output error and controlled increment is realized indirectly [15].

The performance index is supposed as follows:

9
J=12Prk-yk2+QΔ2uk,

where, P, Q are weighting coefficient of output error and controlled increment respectively, rk, y(k) are reference input and output at the moment k respectively, and Δu(k) is output increment of the controller.

The neurons output is:

10
uk=uk-1+KwP'kxPk+wI'kxIk+wD'kxPk,

where, b0 is the first value of the output response, and wi(k) is the weighted value corresponding to the neurons input, and K is the proportion of neurons coefficient, and K>0.

11
wP'k=wPkwPk+wIk+wDk,wI'k=wIkwPk+wIk+wDk,wD'k=wDkwPk+wIk+wDk,
12
xPk=ek-ek-1,xIk=ek,xDk=Δ2ek=ek-2ek-1+ek-2,zk=ek.

3.2. Vehicle steering control algorithm

At a constant speed, the driver’s desired trajectory of vehicle, in fact, is determined by the decision and track of the ideal lateral acceleration, and there is a one-to-one correspondence between the steering wheel angle changes and lateral acceleration changes. In this paper, two quadratic performance indexes, the trajectory following error and the busyness degree of steering wheel, are used to reflect vehicle handling stability, and their expressions are:

13
J1=12[y¨*(k)-y¨(k)]2,
14
J2=12[δsw(k)-δsw(k-1)]2.

Therefore, the performance index function for vehicle steering single neuron controller is:

15
J=12P[y¨*(k)-y(k)]2+Q[δsw(k)-δsw(k-1)],

where, y¨*(k), y¨(k) are the ideal lateral acceleration input and the actual lateral acceleration output at the moment k respectively, δsw(k), δsw(k-1) are the steering wheel angle at the moment k and k-1 respectively, and P, Q are the weighting coefficient of the trajectory following error and the busyness degree of steering wheel.

Gradient descent method is a kind of iterative optimizing algorithms which is seeking the minimum of a particular performance index. Its basic idea is to make the correction of neurons’ connective weights wi(k) decrease along with direction of reduction of J, that is:

16
Δwik-1=wik-wik-1=-ηiJwik-1.

According to the established performance index function J, gradient descent method is adopted to make optimizing design for the performance index function, which objective is to minimize J, so as to realize the correction of neurons’ connective weights. The control rule and the learning algorithm of single neuron self-adaptive controller are as follows:

17
ek=y¨*k-y¨k,xPk=ek-ek-1,xIk=ek,xDk=Δ2ek=ek-2ek-1+ek-2,
18
Twx(k-1)=wI(k-1)xI(k-1)+wP(k-1)xP(k-1)+wD(k-1)xD(k-1),wP(k)=wP(k-1)+ηPK[Pb0e(k)xP(k-1)-QKTwx(k-1)xP(k-1)],wI(k)=wI(k-1)+ηIK[Pb0e(k)xI(k-1)-QKTwx(k-1)xI(k-1)],wD(k)=wD(k-1)+ηDK[Pb0e(k)xD(k-1)-QKTwx(k-1)xD(k-1)],
19
w'P(k)k=wPkwIk+wPk+wDk,wI'k=wIkwIk+wPk+wDk,wD'k=wDkwIk+wPk+wDk,
20
δswk=δswk-1+KwP'kxPk+wP'kxIk+wD'kxDk,

where, ηP, ηI, ηD are proportional, integral and differential learning rate respectively, K is neuron proportion coefficient, and b0 is the output response value of the control system at the initial state.

3.3. Vehicle steering control model

Transfer function cannot describe single neuron PID controller directly, and the simulation cannot be carried on by applying Matlab/Simulink simply, so the function S is introduced to define the system dynamic characteristics, and pid_s module is built. A nonlinear 2DOF vehicle module and a single neuron self-adaptive PID module are connected to constitute the complete control system diagram, and computer simulation model is established by using Matlab/Simulink, as shown in Fig. 6. The control of the nonlinear 2DOF vehicle system by applying the single neuron self-adaptive control algorithm which is based on quadratic form is achieved, and precise tracking for single lane-change path comes true [16].

Fig. 6Computer simulation model

Computer simulation model

3.4. Simulation analysis

In order to simulate the vehicle emergency avoidance at high speed, driver’s preview-follower method which is aimed to determine the ideal lateral acceleration is directly used, and vehicle single lane-change condition (simulating emergency avoidance process) simulation calculation are conducted with 2DOF vehicle model which contains tire nonlinear cornering characteristics at different speeds from low to high. Meanwhile, in order to reflect the robustness of the algorithm while speed changes, proportional, integral and differential learning rate ηP, ηI, ηD at different speeds. The neuron proportion coefficient K, the weighting coefficient of the trajectory following error P and the busyness degree of steering wheel Q are all constant. Vehicle parameters are presented in Table 1.

The situations of tracking the specified path at the speed u=60 km/h,u=80 km/h,u=100 km/h are presented in Fig. 7-9.

Table 1Parameters of vehicle

Parameter
Value
Vehicle mass (M)
1250 kg
Suspension mass (Ms)
900 kg
Front axle to CG (a)
1.170 m
Rear axle to CG (b)
1.195 m
Moment of inertia around Z-axis (Iz)
1800 kg·m2
Effective cornering stiffness of front wheel (Kf)
–60042 N/rad
Effective cornering stiffness of rear wheel (Kr)
–60053 N/rad
Transmission ratio of steering system (i)
20

Fig. 7Tracking the specified path and lateral displacement error of 60 km/h before optimization

Tracking the specified path and lateral displacement error of 60 km/h before optimization

a)

Tracking the specified path and lateral displacement error of 60 km/h before optimization

b)

Fig. 8Tracking the specified path and lateral displacement error of 80 km/h before optimization

Tracking the specified path and lateral displacement error of 80 km/h before optimization

a)

Tracking the specified path and lateral displacement error of 80 km/h before optimization

b)

Fig. 9Tracking the specified path and lateral displacement error of 100 km/h before optimization

Tracking the specified path and lateral displacement error of 100 km/h before optimization

a)

Tracking the specified path and lateral displacement error of 100 km/h before optimization

b)

From Fig. 7-9, the maximum lateral displacement error at the speed of 60 km/h, 80 km/h and 100 km/h is about 0.1216 m, 0.1615 m and 0.1909 m, respectively. On the whole, the vehicle tracks the given trajectory well, which proves that the vehicle has good steering handling trajectory following performance, and single neuron PID control has stronger robustness and adaptability for speed changes. But with the increase of speed, trajectory following error also increases. The main reasons are the following two points.

(1) While emphasizing tracking accuracy, the trajectory following error and the busyness degree of steering wheel are considered comprehensively.

(2) With the increase of speed, the nonlinear characteristics of vehicle become more obvious, but the precision of the single neuron self-adaptive control algorithm declines slightly.

4. Nonlinear vehicle steering control optimization

4.1. Multi-island genetic algorithm

In this paper, control parameter optimizing design through multi-island genetic algorithm is used to realize the optimal control of the control system [14]. Multi-island genetic algorithm is a kind of pseudo parallel genetic algorithm. Compared with traditional genetic algorithms, it divides the whole evolution population into several subgroups which are called “islands”. Then “migration” operation to subgroups in every island is done, which transfers those subgroupsto other islands. The above measure can keep the diversity of population, so the precocious phenomenon of evolution is suppressed. The optimizing procedures are as follows:

(1) Determine the optimal design variables, and then establish the optimizing model.

(2) Determine the coding method, genetic control parameters, population’s size, terminal conditions, and fitness function.

(3) Initialize the population, and conduct fitness evaluation.

(4) Judge whether the terminal conditions meet the requirements, if meet, output optimal results, otherwise, turn to step (5) and conduct genetic operations to create new species.

(5) Genetic operations, according to the previous genetic control parameters, are that selecting individuals randomly to do crossover and mutation operations in sub populations. Then optimal save operation for the previous generation’s optimal individual will be done, and finally produce the next generation.

4.2. Simulation analysis

Based on the “inverse problem” solving, the optimizing model is:

21
minJ,
22
s.t. b_ibib-i, (i=1,2,,6)

where, bi is the control parameter.

In the process of optimization, some high speed evaluation indexes need to be punished, so a punishing item 1+(Tc/T^c)2+(Tu/D^)4 is added. The vehicle handling performance evaluation index takes the following form:

23
J=121+Tc/T^c2+Tu/D^4Py¨*k-yk2+Qδswk-δswk-12.

In the optimizing process, for every set of optimized parameter values, a given track is followed. After optimization by using genetic algorithm, the situation of tracking the specified path at the speed of u=60 km/h,u=80 km/h,u=100 km/h, and path tracking error comparison before and after optimization are shown in Fig. 10-12.

From Fig. 10 to Fig. 12, maximum lateral displacement error at the speed u=60 km/h, u=80 km/h, u=100 km/h are 0.0454 m, 0.0610 m, 0.0878 m respectively. Comparing the above values with the pre-optimized values 0.1216 m, 0.1615 m, 0.1909 m, they reduce significantly. To ensure the maneuverability of the vehicle, the comparison between the optimized Fig. 10-12 and the pre-optimized Fig. 10-12 show that the trajectory error decreases, so the tracking accuracy improves.

Fig. 10Tracking the specified path and lateral displacement error of 60 km/h after optimization

Tracking the specified path and lateral displacement error of 60 km/h after optimization

a)

Tracking the specified path and lateral displacement error of 60 km/h after optimization

b)

Fig. 11Tracking the specified path and lateral displacement error of 80 km/h after optimization

Tracking the specified path and lateral displacement error of 80 km/h after optimization

a)

Tracking the specified path and lateral displacement error of 80 km/h after optimization

b)

Fig. 12Tracking the specified path and lateral displacement error of 100 km/h after optimization

Tracking the specified path and lateral displacement error of 100 km/h after optimization

a)

Tracking the specified path and lateral displacement error of 100 km/h after optimization

b)

5. Conclusions

(1) Aiming at the present situation of less research on nonlinear vehicle high-speed emergency avoidance, tire nonlinear cornering model is established and the validity of the model is verified by the real vehicle test in this paper. On the basis of the above mentioned research, avoidance conditions for a single lane change is set and single neuron self-adaptive PID algorithm is put forward which is based on quadratic form. Through the above measures, the trajectory tracking is realized, and the relatively accurate control performance is reflected, as well as robustness and adaptability to the change of speed.

(2) In accordance with related theories, this paper applies multi-island genetic algorithm to optimize single neuron PID control algorithm, and the results show that the method can improve the accuracy of path tracking effectively at different speeds. The research not only proves the feasibility of genetic algorithm in the study of nonlinear vehicle high-speed emergency avoidance, but also shows that the measure in this paper has certain reference value to the problem of auxiliary driving.

References

  • Mas A., Merienne F., Kemeny A. Lateral control assistance and driver behavior in emergency situations. Advances in Transportation Studies, Special, 2011, p. 149-158.
  • Hattori Y., Ono E., Hosoe S. Optimum vehicle trajectory control for obstacle avoidance problem. IEEE/ASME Transactions on Mechatronics, Vol. 11, Issue 5, 2006, p. 507-512.
  • Kurzhanski A. B. Dynamic optimization for nonlinear target control synthesis. The 6th IFAC Symposium on Nonlinear Control Systems, Stuttgart, Germany, 2004, p. 20-34.
  • Richards A., How J. Aircraft trajectory planning with collision avoidance using mixed integer linear programming. American Control Conference, Anchorage, AK, USA, 2002, p. 1936-1941.
  • Ryu J., Kim H. S., Kim J. H. An emergency obstacle avoidance control strategy for automated highway vehicles. Vehicle System Dynamics, Vol. 38, Issue 5, 2002, p. 319-339.
  • Hagenaars H., Imura J., Nijmeijer H. Approximate continuous-time optimal control in obstacle avoidance by time/space discretization of non-convex state constrains. IEEE Conference on Control Applications, Taipei, Taiwan, China, 2004, p. 878-883.
  • Mukai M., Kawabe T., Nishira H., Takagi Y., Deguchi Y. On vehicle path generation method for collision avoidance using mixed integer programming. The 16th IEEE International Conference on Control Applications Part of IEEE Multi-conference on Systems and Control, Singapore, 2007, p. 1371-1374.
  • Darguzis A., Sapragonas J., Pilkauskas K. Dynamic processes of a vehicle moving over step-shaped obstacles. Journal of Vibroengineering, Vol. 13, Issue 3, 2011, p. 536-343.
  • Ghaffari A., Khodayari A., Arvin S., Alimardani F., Lagnemma K. An ANFIS design for prediction of future state of a vehicle in lane change behavior. IEEE International Conference on Control System, Computing and Engineering, Penang, Malaysia, 2011, p. 156-161.
  • Jin S., Huang Z., Tao P., Wang D. Car-following theory of steady-state traffic flow using time-to-collision. Journal of Zhejiang University-Science A, Applied Physics & Engineering, Vol. 12, Issue 8, 2011, p. 645-654.
  • Konghui G., Ye Z., Dang L. A study on speed-dependent tyre-road friction and its effect on the force and the moment. Vehicle System Dynamics, Vol. 43, 2005, p. 329-340.
  • Konghui G., Dang L., ShiKen C. The unitire model: a nonlinear and non-steady-state tyre model for vehicle dynamics simulation. Vehicle System Dynamics, Vol. 43, 2005, p. 341-358.
  • Konghui G., Nan X., Dang L. A model for combined tire cornering and braking forces with anisotropic tread and carcass stiffness. SAE International Journal of Commercial Vehicles, Vol. 4, Issue 1, 2011, p. 84-95.
  • Umsrithong A., Sandu C. Stochastic transient tyre model for vehicle dynamic simulations on rough rigid ground. International Journal of Vehicle Systems Modeling and Testing, Vol. 7, Issue 4, 2012, p. 351-371.
  • Wang W., Zhou J., Fan, Y. Parameters optimization of laser shot peening based on multi-island genetic algorithm. Applied Mechanics and Materials, Vol. 43, 2011, p. 387-390.
  • Hanxu S., Hou K., Qingxuan J. Development, analysis and control of a spherical aerial vehicle. Journal of Vibroengineering, Vol. 15, Issue 2, 2013, p. 1069-1080.

About this article

Received
25 September 2013
Accepted
21 December 2013
Published
31 March 2014
Keywords
genetic algorithm
nonlinear
emergency avoidance
steering control
parameter optimization
Acknowledgements

This work was supported in part by the National Science Foundation of China (Grant No. 11072106).