# Design of an Integrated Vehicle Chassis Control System with Driver Behavior Identification

Authors: Bing Zhu; Yizhou Chen; Jian Zhao; Yunfu Su
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

An integrated vehicle chassis control strategy with driver behavior identification is introduced in this paper. In order to identify the different types of driver behavior characteristics, a driver behavior signals acquisition system was established using the dSPACE real-time simulation platform, and the driver inputs of 30 test drivers were collected under the double lane change test condition. Then, driver behavior characteristics were analyzed and identified based on the preview optimal curvature model through genetic algorithm and neural network method. Using it as a base, an integrated chassis control strategy with active front steering (AFS) and direct yaw moment control (DYC) considering driver characteristics was established by model predictive control (MPC) method. Finally, simulations were carried out to verify the control strategy by CarSim and MATLAB/Simulink. The results show that the proposed method enables the control system to adjust its parameters according to the driver behavior identification results and the vehicle handling and stability performance are significantly improved.


---

## Body

## 1. Introduction

In recent years, improving the active safety performance of vehicles is what researchers have been working on. With the development of electronic control technology, many active safety systems such as antilock braking system (ABS), traction control system (TCS), electronic stability control (ESC), active front steering (AFS), and four-wheel steering (4WS) have been widely equipped on vehicles to ensure safer and more stable driving experience. However, the potential conflicts among the active control systems increase when they are combined without coordination [1–3]. Thus, there have been plenty of attempts to integrate the chassis control subsystems, for instance, the Integrated Chassis Control (ICC), Unified Chassis Control (UCC), and Vehicle Dynamics Management (VDM) [4–6]. With the implementation of integrated control systems, the interference and coupling among dynamic subsystems are effectively eliminated and the stability performance of vehicle is significantly improved. Nevertheless, integrated control systems nowadays are generally designed in the uniform mode, lacking the consideration of the influence that drivers exert on control systems. As a matter of fact, the driver and the active control system have strong coupling on each other while controlling the vehicle. According to a survey sponsored by National Highway Traffic Safety Administration (NHTSA) of USA, driver’s mistake accounts for 90% or more of all the crashes recorded [7]. Therefore, human driver’s characteristics should be involved in the process where the integrated control system is developed.As the primary control element within the traditional driver-vehicle system, the role of human driver has been studied by a plenty of scientists and researchers. The driver behavior includes driver’s sensing, judging, reasoning, deciding, and finally operating the vehicle to turn, accelerate, and brake with strong randomicity, adaptivity, discreteness, and variability. Currently, research is mainly carried out in the aspects of driver behavior modeling, the identification of fatigue, intelligent traffic control, and advanced driver assistance systems (ADAS) [8].Koh et al. present a tire slip-angle based speed control race driver model through analyzing the vehicle-driver interaction at limit handling [9]. Miyajima et al. modeled driver’s behavior such as turning the steering wheel or hitting the pedal with Gaussian mixture model (GMM) to identify different kinds of drivers. Compared with the method using raw pedal operation signals, spectral analysis method proposed shows a far better accuracy [10]. Bolia et al. proposed a two-level preview driver steering control model. The outer loop focuses on path following while the inner one tries to capture driver’s physical behavior [11]. Lin et al. used several sophisticated artificial neural network architectures to develop driver models in a Driver-Vehicle-Environment (DVE) system [12]. Sathyanarayana et al. proposed a “context and driver aware” (CDA) active vehicle system combining GMM, universal background model, and likelihood maximization based on information fusion to identify driver status and predict driver’s distracted behavior [13]. Profound and enlightening works above endeavored to explain what kind of a role driver is playing in a driver-vehicle closed-loop system. With further research, early efforts have been made to design the active control system on the basis of driver’s characteristics. Macadam studied the physiological limits and physical characteristics of drivers in detail and then established the lateral and longitudinal manipulating model of drivers [14]. Hoult and Cole built a neuromuscular linear model which considers the coactivation of neuromuscular system, muscle, body, and vehicle [15]. Chai et al. described a method to adjust the parameters of steer-by-wire (SBW) system according to the driver’s steer characteristics, which could be estimated from experimental data based on the general driver’s model [16]. Raksincharoensak et al., Japanese researchers, designed the direct yaw moment control system based on the identification of driver’s intention and the performance of the vehicle is significantly improved [17]. In the literature [18], a driver behavior signal capturing system was introduced, based on which the riding comfort and safety of the vehicle are enhanced. Fu et al. integrated a driver model with a run-off-road recovery controller considering driver’s target planning, pursuit behavior, compensate behavior, and physical limitations [19]. Keen and Cole proposed a steering controller based on linear model predictive control. A formal system identification procedure is applied to avoid bias from the closed-loop operation of the driver-vehicle system [20]. Although the driver behavior identification methods have made remarkable progress, in practical applications, the driver behavior observation equipment mentioned above is seldom used due to its disadvantages in prize and portability.The active chassis control systems and the driver have mutual influence for vehicle control with strong couplings. Due to the complexity, polytrope, and uncertainty of driver behavior, it is extremely essential to fully understand the characteristics of drivers and control systems, build the collaborative optimization mechanism, and optimize the performance of the Driver-Vehicle-Environment system.To design the personalized vehicle control system, one of the most important prerequisites is identifying the driver’s individualities. On the basis of the literature [14], Professor Guo proposed a “Preview-Following” system theory, which defines the decision procedure of drivers as preview and compensation. According to this theory, Preview Optimal Curvature Model is established. As it demonstrates clearly a relationship between the vehicle stability and driver’s characteristics, this theory has been widely used in the research of vehicle stability control and adaptive cruise control due to its high precision and independence of special equipment [21]. Simulation and test results show that this theory precisely depicts driver’s steering procedure. The key to applying this theory is the identification of driver’s steering characteristics.In this paper, an integrated chassis control strategy based on the identification of driver characteristics is proposed. First, a real-time driver operation signal acquisition system based on dSPACE is built. Using this system, 10 skilled, 10 normal, and 10 novice drivers’ manipulation signals on the same path are collected. Then, genetic algorithm is applied to analyzing the steering characteristics of different kinds of drivers based on the Preview Optimal Curvature Model. An online driver behavior identification method is established using BP neural network. An integrated chassis control (ICC) system with active front steering (AFS) and direct yaw moment control (DYC) is designed using model predictive control (MPC) and considering the driver behavior characteristics. Finally, simulations are carried out to validate the proposed method through the cosimulation of MATLAB/Simulink and CarSim. Test results show that the proposed method enables the control system to adjust its parameters according to the driver behavior identification results and the vehicle stability is effectively enhanced.


## 2. Driver Behavior Data Acquisition

### 2.1. Driver Behavior Data Acquisition System

In order to precisely analyze the driver behavior, the driver behavior data acquisition system is designed and established based on dSPACE real-time simulation platform, as shown in Figure1. The dSPACE real-time simulation system could realize the seamless connection to Matlab/Simulink through the automatic code generation and downloading Real-time Interface (RTI). And with the test software ControlDesk, the dSPACE can accomplish the visual management and automatic control of the tests.Figure 1
Structure of the driver behavior data acquisition system.In this system, dSPACE DS1006 is used as real-time simulation processor, which collects the opening of accelerator pedal and brake pedal along with the steering wheel angel in real-time and meanwhile feeds back the signals to the CarSim vehicle model which runs on an industrial personal computer (IPC). CarSim generates a virtual scene projected to a large screen, according to which the driver perceives the driving environment and operates the vehicle.In order to enhance the immersion, the torque steering wheel made by SENSODRIVE is adopted in the system. The steering wheel, possessing adjustable variables of damp and friction, can simulate the torque feedback of the real road and greatly enhances the driving experience. The steering wheel is connected to dSPACE DS1006 through CAN bus. The linear potentiometer displacement sensors are used for the acquisition of the accelerator pedal and the brake pedal signals through the AD ports.


### 2.2. Driver Behavior Data Acquisition

In order to identify the characteristics of different drivers, the usage modes of different kinds of drivers are analyzed. In normal working conditions, the usage modes of different drivers tend to be similar and difficult to distinguish. Thus, driving data acquisition experiments are carried out under the double lane change (DLC) working conditions. The target path is shown in Figure2.Figure 2
Working condition used in driver behavior data acquisition.In the experiments, to guarantee the validation of the test data, 20 male drivers and 10 female drivers are chosen as the test samples. Among the drivers, 10 are skilled drivers, 10 are normal ones, and the last 10 are novice in driving. Before the experiment started, all drivers are trained to be acquainted with the experiment process and the data collection system. When the data collection begins, the driver accelerates from the same start point to the speed of 60 Km/h, 80 Km/h, 100 Km/h, 110 Km/h, and 120 Km/h, respectively, at the end of lane line 1. Then, the lane change is carried out according to the target path. When the car passes lane line 3, the driver decelerates the car to stop.Typical collection signals of three different kinds of drivers are shown in Figure3. It is seen that, facing the panic lane change, the skilled driver shows shorter response time and could implement lane change smoothly. In contrast, the novice driver tends to operate the vehicle sharply and the vehicle may lose its stability.Figure 3
Typical signals of different kinds of drivers.
(a)
Accelerator pedal input
(b)
Brake pedal input
(c)
Steering wheel input
(d)
Vehicle speed


## 2.1. Driver Behavior Data Acquisition System

In order to precisely analyze the driver behavior, the driver behavior data acquisition system is designed and established based on dSPACE real-time simulation platform, as shown in Figure1. The dSPACE real-time simulation system could realize the seamless connection to Matlab/Simulink through the automatic code generation and downloading Real-time Interface (RTI). And with the test software ControlDesk, the dSPACE can accomplish the visual management and automatic control of the tests.Figure 1
Structure of the driver behavior data acquisition system.In this system, dSPACE DS1006 is used as real-time simulation processor, which collects the opening of accelerator pedal and brake pedal along with the steering wheel angel in real-time and meanwhile feeds back the signals to the CarSim vehicle model which runs on an industrial personal computer (IPC). CarSim generates a virtual scene projected to a large screen, according to which the driver perceives the driving environment and operates the vehicle.In order to enhance the immersion, the torque steering wheel made by SENSODRIVE is adopted in the system. The steering wheel, possessing adjustable variables of damp and friction, can simulate the torque feedback of the real road and greatly enhances the driving experience. The steering wheel is connected to dSPACE DS1006 through CAN bus. The linear potentiometer displacement sensors are used for the acquisition of the accelerator pedal and the brake pedal signals through the AD ports.


## 2.2. Driver Behavior Data Acquisition

In order to identify the characteristics of different drivers, the usage modes of different kinds of drivers are analyzed. In normal working conditions, the usage modes of different drivers tend to be similar and difficult to distinguish. Thus, driving data acquisition experiments are carried out under the double lane change (DLC) working conditions. The target path is shown in Figure2.Figure 2
Working condition used in driver behavior data acquisition.In the experiments, to guarantee the validation of the test data, 20 male drivers and 10 female drivers are chosen as the test samples. Among the drivers, 10 are skilled drivers, 10 are normal ones, and the last 10 are novice in driving. Before the experiment started, all drivers are trained to be acquainted with the experiment process and the data collection system. When the data collection begins, the driver accelerates from the same start point to the speed of 60 Km/h, 80 Km/h, 100 Km/h, 110 Km/h, and 120 Km/h, respectively, at the end of lane line 1. Then, the lane change is carried out according to the target path. When the car passes lane line 3, the driver decelerates the car to stop.Typical collection signals of three different kinds of drivers are shown in Figure3. It is seen that, facing the panic lane change, the skilled driver shows shorter response time and could implement lane change smoothly. In contrast, the novice driver tends to operate the vehicle sharply and the vehicle may lose its stability.Figure 3
Typical signals of different kinds of drivers.
(a)
Accelerator pedal input
(b)
Brake pedal input
(c)
Steering wheel input
(d)
Vehicle speed


## 3. Analysis of Driver Behavior Characteristics

### 3.1. Preview Optimal Curvature Model

As is shown in Figure3, among the drivers’ operating signals, the steering wheel angles show the most significant difference among three kinds of drivers. Thus, in this research, steering wheel angle is chosen to identify behavior characteristics of different kinds of drivers. Preview optimal curvature model, which is widely used in driver-vehicle closed-loop system research as mentioned above, is introduced here to analyze the driver steering behavior [22].It is widely believed that the driver plans for the desired path when he/she is driving a vehicle. Assume that the desired path isY=f(X) and the preview time of driver is a constant, as is shown in Figure 4. At time t, the state of the vehicle is(1)y=yt,y˙=y˙t.Figure 4
Frame of driver’s steering decision.On the condition that the preview distance isd, namely, the driver always gazes at preview point which is d away in front of the vehicle, the preview time could be defined as(2)T=dVx,where Vx is the longitudinal velocity of the vehicle.The lateral coordinate of the preview point isf(t+T). If the driver operates a steering angel δsw, in accordance with which the curvature of the vehicle is 1/R, and the lateral acceleration is y¨(t); in time T, the lateral displacement of the vehicle is(3)ft+T=yt+T·y˙t+T22y¨t.Considering(4)y¨=Vx2R=δswiLVx2,among which L is the wheel base and i is the steering ratio.We get the optimal steering wheel angle:(5)δsw=2iLd2ft+T-yt-T·y˙t.Drivers always intend to operate a steering wheel angleδsw under which the displacement of the vehicle could fit the desired path in time T. Considering the lag of driver’s perception and action together with the dynamic response of the vehicle under a sharp steering, the steering wheel angle δsw that the driver is supposed to operate to track the target path f is(6)δswfs=2C01+Tcse-tdss2Tp21+Thss2+2C01+Tcs1+Tpse-tdsay/δsw,where Tp is the driver’s preview time, td is the driver’s neural delay time, Th is the driver’s muscle delay time, Tc is the correction parameter of steering, and C0 is the proportional gain correction parameter. The relationships among the parameters are(7)Tc=td+Th+Ta-aTp3,C0=1Gay,where Ta is a parameter of the vehicle and a is the following order. Gay is the steady-state gain of lateral acceleration of ay to δsw whose dynamic response is defined as the vehicle model V(s), which will be described in the next section.Therefore, the frame of driver preview optimal curvature model could be demonstrated as in Figure5, from which it is clearly seen that the parameters Tp, td, and Th could characterize different kinds of drivers’ behavior, and, namely, the task of identification of different kinds of drivers’ steering behavior is confirming the three parameters.Figure 5
Driver preview optimal curvature model.


### 3.2. Identification of Driver Behavior Characteristics Parameters

#### 3.2.1. Identification of the Parameters of Vehicle Model

As shown in Figure5, before identifying the driver’s parameters, the characteristics of the vehicle should be specified.According to the preview optimal curvature model,V(s) could be described as the dynamic response of ay to δsw:(8)Vs=ayδsws=Gay1+Ty1s+Ty2s21+T1s+T2s2.In the ideal situation,V(s) could be directly described as Gay. Considering that the vehicle’s dynamic response could not be ignored when it is taking a sharp turning, Ty1, Ty2, T1, and T2 are constants to be identified and Ta in (7) is defined as (T1-Ty1).The aim of parameters identification is to describe the mathematic relation between the input and output of the system. Since the mechanism of the vehicle dynamic response is not necessarily required in this research, ARX model, which does not demand the explicit physical relationship of the system and works well when dealing with the higher-order system, is adopted here to identify the parameters of the vehicle model. The expression of the ARX model is shown below:(9)AzAyz=Bzuz+ez,where(10)Az=1+a1z-1+a2z-2+⋯+anz-nBz=1+b1z-1+b2z-2+⋯+bmz-m.Here, n,m=2, e(z) is white noise, u(z) is δsw, and Ay(z) is ay.Using the MATLAB toolbox to compute, the real lateral acceleration and the output of the identification model are shown in Figure6 which demonstrates that the vehicle model built here has a high accuracy and the identification rate is 97.2%.Figure 6
Output of vehicle model and real lateral acceleration.


#### 3.2.2. Offline Driver Model Parameters Identification for Different Kinds of Drivers

On the basis of the identification of the vehicle modelV(s), driver behavior characteristic parameters, Tp, td, and Th, are to be identified. Genetic algorithm, which simulates the mechanism of natural selection and biological evolutionism, has been widely used in pattern recognition, computing science, automatic control, and so forth. Here, genetic algorithm is adopted to identify the driver model parameters offline.Firstly, the three model parameters of drivers should be constrained to appropriate scope to guarantee the accuracy of identification results as well as accelerate the identification procedure. As is demonstrated in the literature [22], the parameters are supposed to meet the requirements:(11)Tj∈Ujmin,Ujmax,where Tj is Tp, td, and Th. Ujmin and Ujmax are the minimum and maximum of Tp, td, and Th, respectively. j=1,2,3.The parameters are constrained to a certain scope shown in Table1.Table 1
Constraints of driver behavior characteristic parameters.
T
p
t
d
T
h
[
0.58
,
2.072
]
[
0.17
,
0.53
]
[
0.08
,
0.36
]To ensure the precision of identification of driver behavior characteristic and that the dynamic responses of the driver model accord well with the experimental data, the target function is defined as the weighted value of the error square sum of steering wheel angel, lateral acceleration, and lateral displacement:(12)J=w1∗J1+w2∗J2+w3∗J3,where J1, J2, and J3 are the deviation indices of the steering wheel angel, the lateral acceleration, and the lateral displacement, respectively. Their expressions are shown as(13)J1=1t0∫0t0δsw-δ^sw2dt,J2=1t0∫0t0ay-a^y2dt,J3=1t0∫0t0y-y^2dt,where δ^sw, a^y, and y^ are the theoretical data based on driver model and vehicle model. w1, w2, and w3 are weight coefficients and t0 is the sampling time.Since there are three variables, the fitness function is determined according to the simple adaptive function method. Besides, the binary encoding method, proportion selection strategy, single point switching strategy, and simple mutation strategy are adopted [23].As mentioned in Section2, the drivers can be divided into three categories, the skilled ones, normal ones, and novice ones. With genetic algorithm, the driver’s characteristics could be identified quantitatively and the results of the identification of parameters are listed in Table 2.Table 2
Identification results of driver behavior characteristic parameters.
Driver
Average
Standard deviation
Target variance
Novice drivers
Tp
0.7100
0.0964
4.7016
e
-
003
td
0.1489
0.0834
Th
0.3740
0.0637
Skilled drivers
Tp
1.6954
0.0961
5.5144
e
-
003
td
0.3763
0.0720
Th
0.3632
0.0560
Normal drivers
Tp
1.1534
0.1015
7.5145
e
-
003
td
0.2388
0.1132
Th
0.2577
0.0873When the offline identification is finished, the results are adopted as the training output of the online BP neural network model used in online characteristics identification.


#### 3.2.3. Online Driver Behavior Characteristic Parameters Identification

Due to the genetic algorithm’s relative time-consuming, low efficiency, and absence of meeting the real-time requirement of the control system, in this research, BP neural network algorithm is designed for online identification of driver behavior characteristic parameters. BP neural network is established as in Figure7.Figure 7
Structure of the BP neural network model.The BP neural network model has 3 layers; the number of nodes of the input layer, hidden layer, and output layern, l, and m is 5, 8, and 3, respectively. The five nodes x1~x5 in input layer are δsw, f, y, ay, and Vx, respectively, δsw is driver’s steering wheel angle, f is the lateral displacement of the target path, y is the real lateral displacement of the vehicle, ay is the lateral acceleration of the vehicle, and Vx is the longitudinal velocity of the vehicle. The three nodes z1~z3 in output layer are Tp, td, and Th, which are the drive’s behavior characteristics to be identified.The output of the hidden layer is(14)Hj=g∑i=1nωijxi-hj,j=1,2,…,l.The output of the output layer is(15)Ok=∑j=1lωjkHj-ok,k=1,2,…,m,where ωij and ωjk are weighting coefficients of input layer to hidden layer and hidden layer to output layer, respectively. The thresholds of the hidden layer and output layer are h and o, respectively. The g(x) is the excitation function:(16)gx=11+e-x.The deviation calculation function is used to calculate the deviation of the desired output and the calculated output of the model:(17)ek=12Yk-Ok2,k=1,2,…,m,where Yk is the desired output.In order to minimize the variance of the output, back-propagation modification weight matrix is used as the learning rule of the model:(18)ωijnew=ωij+ηHj1-Hjxi∑k=1mωjkeki=1,2,…,n;j=1,2,…,l,ωjknew=ωjk+ηHjek,j=1,2,…,l;k=1,2,…,mhjnew=hj+ηHj1-Hj∑k=1mωjkek,j=1,2,…,loknew=ok+ek,k=1,2,…,m,where ωijnew and ωjknew are weighting coefficients after update and hjnew and oknew are thresholds after update. η is the learning rate.The offline identification results obtained from genetic algorithm are used to train the BP neural network model. And when the training performance is as shown in Figure8, the training is stopped. The datasheet trained from neural network model could be directly used as driver behavior characteristic parameters identification model. Therefore, the driver model whose characteristics fit the current driver’s behavior is achieved.Figure 8
Performance of neural network model.Thus, when a driver is driving a car, his/her behavior characteristic parameters could be identified rapidly, on the basis of which the control system could be designed to fit the driver’s characteristics.


## 3.1. Preview Optimal Curvature Model

As is shown in Figure3, among the drivers’ operating signals, the steering wheel angles show the most significant difference among three kinds of drivers. Thus, in this research, steering wheel angle is chosen to identify behavior characteristics of different kinds of drivers. Preview optimal curvature model, which is widely used in driver-vehicle closed-loop system research as mentioned above, is introduced here to analyze the driver steering behavior [22].It is widely believed that the driver plans for the desired path when he/she is driving a vehicle. Assume that the desired path isY=f(X) and the preview time of driver is a constant, as is shown in Figure 4. At time t, the state of the vehicle is(1)y=yt,y˙=y˙t.Figure 4
Frame of driver’s steering decision.On the condition that the preview distance isd, namely, the driver always gazes at preview point which is d away in front of the vehicle, the preview time could be defined as(2)T=dVx,where Vx is the longitudinal velocity of the vehicle.The lateral coordinate of the preview point isf(t+T). If the driver operates a steering angel δsw, in accordance with which the curvature of the vehicle is 1/R, and the lateral acceleration is y¨(t); in time T, the lateral displacement of the vehicle is(3)ft+T=yt+T·y˙t+T22y¨t.Considering(4)y¨=Vx2R=δswiLVx2,among which L is the wheel base and i is the steering ratio.We get the optimal steering wheel angle:(5)δsw=2iLd2ft+T-yt-T·y˙t.Drivers always intend to operate a steering wheel angleδsw under which the displacement of the vehicle could fit the desired path in time T. Considering the lag of driver’s perception and action together with the dynamic response of the vehicle under a sharp steering, the steering wheel angle δsw that the driver is supposed to operate to track the target path f is(6)δswfs=2C01+Tcse-tdss2Tp21+Thss2+2C01+Tcs1+Tpse-tdsay/δsw,where Tp is the driver’s preview time, td is the driver’s neural delay time, Th is the driver’s muscle delay time, Tc is the correction parameter of steering, and C0 is the proportional gain correction parameter. The relationships among the parameters are(7)Tc=td+Th+Ta-aTp3,C0=1Gay,where Ta is a parameter of the vehicle and a is the following order. Gay is the steady-state gain of lateral acceleration of ay to δsw whose dynamic response is defined as the vehicle model V(s), which will be described in the next section.Therefore, the frame of driver preview optimal curvature model could be demonstrated as in Figure5, from which it is clearly seen that the parameters Tp, td, and Th could characterize different kinds of drivers’ behavior, and, namely, the task of identification of different kinds of drivers’ steering behavior is confirming the three parameters.Figure 5
Driver preview optimal curvature model.


## 3.2. Identification of Driver Behavior Characteristics Parameters

### 3.2.1. Identification of the Parameters of Vehicle Model

As shown in Figure5, before identifying the driver’s parameters, the characteristics of the vehicle should be specified.According to the preview optimal curvature model,V(s) could be described as the dynamic response of ay to δsw:(8)Vs=ayδsws=Gay1+Ty1s+Ty2s21+T1s+T2s2.In the ideal situation,V(s) could be directly described as Gay. Considering that the vehicle’s dynamic response could not be ignored when it is taking a sharp turning, Ty1, Ty2, T1, and T2 are constants to be identified and Ta in (7) is defined as (T1-Ty1).The aim of parameters identification is to describe the mathematic relation between the input and output of the system. Since the mechanism of the vehicle dynamic response is not necessarily required in this research, ARX model, which does not demand the explicit physical relationship of the system and works well when dealing with the higher-order system, is adopted here to identify the parameters of the vehicle model. The expression of the ARX model is shown below:(9)AzAyz=Bzuz+ez,where(10)Az=1+a1z-1+a2z-2+⋯+anz-nBz=1+b1z-1+b2z-2+⋯+bmz-m.Here, n,m=2, e(z) is white noise, u(z) is δsw, and Ay(z) is ay.Using the MATLAB toolbox to compute, the real lateral acceleration and the output of the identification model are shown in Figure6 which demonstrates that the vehicle model built here has a high accuracy and the identification rate is 97.2%.Figure 6
Output of vehicle model and real lateral acceleration.


### 3.2.2. Offline Driver Model Parameters Identification for Different Kinds of Drivers

On the basis of the identification of the vehicle modelV(s), driver behavior characteristic parameters, Tp, td, and Th, are to be identified. Genetic algorithm, which simulates the mechanism of natural selection and biological evolutionism, has been widely used in pattern recognition, computing science, automatic control, and so forth. Here, genetic algorithm is adopted to identify the driver model parameters offline.Firstly, the three model parameters of drivers should be constrained to appropriate scope to guarantee the accuracy of identification results as well as accelerate the identification procedure. As is demonstrated in the literature [22], the parameters are supposed to meet the requirements:(11)Tj∈Ujmin,Ujmax,where Tj is Tp, td, and Th. Ujmin and Ujmax are the minimum and maximum of Tp, td, and Th, respectively. j=1,2,3.The parameters are constrained to a certain scope shown in Table1.Table 1
Constraints of driver behavior characteristic parameters.
T
p
t
d
T
h
[
0.58
,
2.072
]
[
0.17
,
0.53
]
[
0.08
,
0.36
]To ensure the precision of identification of driver behavior characteristic and that the dynamic responses of the driver model accord well with the experimental data, the target function is defined as the weighted value of the error square sum of steering wheel angel, lateral acceleration, and lateral displacement:(12)J=w1∗J1+w2∗J2+w3∗J3,where J1, J2, and J3 are the deviation indices of the steering wheel angel, the lateral acceleration, and the lateral displacement, respectively. Their expressions are shown as(13)J1=1t0∫0t0δsw-δ^sw2dt,J2=1t0∫0t0ay-a^y2dt,J3=1t0∫0t0y-y^2dt,where δ^sw, a^y, and y^ are the theoretical data based on driver model and vehicle model. w1, w2, and w3 are weight coefficients and t0 is the sampling time.Since there are three variables, the fitness function is determined according to the simple adaptive function method. Besides, the binary encoding method, proportion selection strategy, single point switching strategy, and simple mutation strategy are adopted [23].As mentioned in Section2, the drivers can be divided into three categories, the skilled ones, normal ones, and novice ones. With genetic algorithm, the driver’s characteristics could be identified quantitatively and the results of the identification of parameters are listed in Table 2.Table 2
Identification results of driver behavior characteristic parameters.
Driver
Average
Standard deviation
Target variance
Novice drivers
Tp
0.7100
0.0964
4.7016
e
-
003
td
0.1489
0.0834
Th
0.3740
0.0637
Skilled drivers
Tp
1.6954
0.0961
5.5144
e
-
003
td
0.3763
0.0720
Th
0.3632
0.0560
Normal drivers
Tp
1.1534
0.1015
7.5145
e
-
003
td
0.2388
0.1132
Th
0.2577
0.0873When the offline identification is finished, the results are adopted as the training output of the online BP neural network model used in online characteristics identification.


### 3.2.3. Online Driver Behavior Characteristic Parameters Identification

Due to the genetic algorithm’s relative time-consuming, low efficiency, and absence of meeting the real-time requirement of the control system, in this research, BP neural network algorithm is designed for online identification of driver behavior characteristic parameters. BP neural network is established as in Figure7.Figure 7
Structure of the BP neural network model.The BP neural network model has 3 layers; the number of nodes of the input layer, hidden layer, and output layern, l, and m is 5, 8, and 3, respectively. The five nodes x1~x5 in input layer are δsw, f, y, ay, and Vx, respectively, δsw is driver’s steering wheel angle, f is the lateral displacement of the target path, y is the real lateral displacement of the vehicle, ay is the lateral acceleration of the vehicle, and Vx is the longitudinal velocity of the vehicle. The three nodes z1~z3 in output layer are Tp, td, and Th, which are the drive’s behavior characteristics to be identified.The output of the hidden layer is(14)Hj=g∑i=1nωijxi-hj,j=1,2,…,l.The output of the output layer is(15)Ok=∑j=1lωjkHj-ok,k=1,2,…,m,where ωij and ωjk are weighting coefficients of input layer to hidden layer and hidden layer to output layer, respectively. The thresholds of the hidden layer and output layer are h and o, respectively. The g(x) is the excitation function:(16)gx=11+e-x.The deviation calculation function is used to calculate the deviation of the desired output and the calculated output of the model:(17)ek=12Yk-Ok2,k=1,2,…,m,where Yk is the desired output.In order to minimize the variance of the output, back-propagation modification weight matrix is used as the learning rule of the model:(18)ωijnew=ωij+ηHj1-Hjxi∑k=1mωjkeki=1,2,…,n;j=1,2,…,l,ωjknew=ωjk+ηHjek,j=1,2,…,l;k=1,2,…,mhjnew=hj+ηHj1-Hj∑k=1mωjkek,j=1,2,…,loknew=ok+ek,k=1,2,…,m,where ωijnew and ωjknew are weighting coefficients after update and hjnew and oknew are thresholds after update. η is the learning rate.The offline identification results obtained from genetic algorithm are used to train the BP neural network model. And when the training performance is as shown in Figure8, the training is stopped. The datasheet trained from neural network model could be directly used as driver behavior characteristic parameters identification model. Therefore, the driver model whose characteristics fit the current driver’s behavior is achieved.Figure 8
Performance of neural network model.Thus, when a driver is driving a car, his/her behavior characteristic parameters could be identified rapidly, on the basis of which the control system could be designed to fit the driver’s characteristics.


## 3.2.1. Identification of the Parameters of Vehicle Model

As shown in Figure5, before identifying the driver’s parameters, the characteristics of the vehicle should be specified.According to the preview optimal curvature model,V(s) could be described as the dynamic response of ay to δsw:(8)Vs=ayδsws=Gay1+Ty1s+Ty2s21+T1s+T2s2.In the ideal situation,V(s) could be directly described as Gay. Considering that the vehicle’s dynamic response could not be ignored when it is taking a sharp turning, Ty1, Ty2, T1, and T2 are constants to be identified and Ta in (7) is defined as (T1-Ty1).The aim of parameters identification is to describe the mathematic relation between the input and output of the system. Since the mechanism of the vehicle dynamic response is not necessarily required in this research, ARX model, which does not demand the explicit physical relationship of the system and works well when dealing with the higher-order system, is adopted here to identify the parameters of the vehicle model. The expression of the ARX model is shown below:(9)AzAyz=Bzuz+ez,where(10)Az=1+a1z-1+a2z-2+⋯+anz-nBz=1+b1z-1+b2z-2+⋯+bmz-m.Here, n,m=2, e(z) is white noise, u(z) is δsw, and Ay(z) is ay.Using the MATLAB toolbox to compute, the real lateral acceleration and the output of the identification model are shown in Figure6 which demonstrates that the vehicle model built here has a high accuracy and the identification rate is 97.2%.Figure 6
Output of vehicle model and real lateral acceleration.


## 3.2.2. Offline Driver Model Parameters Identification for Different Kinds of Drivers

On the basis of the identification of the vehicle modelV(s), driver behavior characteristic parameters, Tp, td, and Th, are to be identified. Genetic algorithm, which simulates the mechanism of natural selection and biological evolutionism, has been widely used in pattern recognition, computing science, automatic control, and so forth. Here, genetic algorithm is adopted to identify the driver model parameters offline.Firstly, the three model parameters of drivers should be constrained to appropriate scope to guarantee the accuracy of identification results as well as accelerate the identification procedure. As is demonstrated in the literature [22], the parameters are supposed to meet the requirements:(11)Tj∈Ujmin,Ujmax,where Tj is Tp, td, and Th. Ujmin and Ujmax are the minimum and maximum of Tp, td, and Th, respectively. j=1,2,3.The parameters are constrained to a certain scope shown in Table1.Table 1
Constraints of driver behavior characteristic parameters.
T
p
t
d
T
h
[
0.58
,
2.072
]
[
0.17
,
0.53
]
[
0.08
,
0.36
]To ensure the precision of identification of driver behavior characteristic and that the dynamic responses of the driver model accord well with the experimental data, the target function is defined as the weighted value of the error square sum of steering wheel angel, lateral acceleration, and lateral displacement:(12)J=w1∗J1+w2∗J2+w3∗J3,where J1, J2, and J3 are the deviation indices of the steering wheel angel, the lateral acceleration, and the lateral displacement, respectively. Their expressions are shown as(13)J1=1t0∫0t0δsw-δ^sw2dt,J2=1t0∫0t0ay-a^y2dt,J3=1t0∫0t0y-y^2dt,where δ^sw, a^y, and y^ are the theoretical data based on driver model and vehicle model. w1, w2, and w3 are weight coefficients and t0 is the sampling time.Since there are three variables, the fitness function is determined according to the simple adaptive function method. Besides, the binary encoding method, proportion selection strategy, single point switching strategy, and simple mutation strategy are adopted [23].As mentioned in Section2, the drivers can be divided into three categories, the skilled ones, normal ones, and novice ones. With genetic algorithm, the driver’s characteristics could be identified quantitatively and the results of the identification of parameters are listed in Table 2.Table 2
Identification results of driver behavior characteristic parameters.
Driver
Average
Standard deviation
Target variance
Novice drivers
Tp
0.7100
0.0964
4.7016
e
-
003
td
0.1489
0.0834
Th
0.3740
0.0637
Skilled drivers
Tp
1.6954
0.0961
5.5144
e
-
003
td
0.3763
0.0720
Th
0.3632
0.0560
Normal drivers
Tp
1.1534
0.1015
7.5145
e
-
003
td
0.2388
0.1132
Th
0.2577
0.0873When the offline identification is finished, the results are adopted as the training output of the online BP neural network model used in online characteristics identification.


## 3.2.3. Online Driver Behavior Characteristic Parameters Identification

Due to the genetic algorithm’s relative time-consuming, low efficiency, and absence of meeting the real-time requirement of the control system, in this research, BP neural network algorithm is designed for online identification of driver behavior characteristic parameters. BP neural network is established as in Figure7.Figure 7
Structure of the BP neural network model.The BP neural network model has 3 layers; the number of nodes of the input layer, hidden layer, and output layern, l, and m is 5, 8, and 3, respectively. The five nodes x1~x5 in input layer are δsw, f, y, ay, and Vx, respectively, δsw is driver’s steering wheel angle, f is the lateral displacement of the target path, y is the real lateral displacement of the vehicle, ay is the lateral acceleration of the vehicle, and Vx is the longitudinal velocity of the vehicle. The three nodes z1~z3 in output layer are Tp, td, and Th, which are the drive’s behavior characteristics to be identified.The output of the hidden layer is(14)Hj=g∑i=1nωijxi-hj,j=1,2,…,l.The output of the output layer is(15)Ok=∑j=1lωjkHj-ok,k=1,2,…,m,where ωij and ωjk are weighting coefficients of input layer to hidden layer and hidden layer to output layer, respectively. The thresholds of the hidden layer and output layer are h and o, respectively. The g(x) is the excitation function:(16)gx=11+e-x.The deviation calculation function is used to calculate the deviation of the desired output and the calculated output of the model:(17)ek=12Yk-Ok2,k=1,2,…,m,where Yk is the desired output.In order to minimize the variance of the output, back-propagation modification weight matrix is used as the learning rule of the model:(18)ωijnew=ωij+ηHj1-Hjxi∑k=1mωjkeki=1,2,…,n;j=1,2,…,l,ωjknew=ωjk+ηHjek,j=1,2,…,l;k=1,2,…,mhjnew=hj+ηHj1-Hj∑k=1mωjkek,j=1,2,…,loknew=ok+ek,k=1,2,…,m,where ωijnew and ωjknew are weighting coefficients after update and hjnew and oknew are thresholds after update. η is the learning rate.The offline identification results obtained from genetic algorithm are used to train the BP neural network model. And when the training performance is as shown in Figure8, the training is stopped. The datasheet trained from neural network model could be directly used as driver behavior characteristic parameters identification model. Therefore, the driver model whose characteristics fit the current driver’s behavior is achieved.Figure 8
Performance of neural network model.Thus, when a driver is driving a car, his/her behavior characteristic parameters could be identified rapidly, on the basis of which the control system could be designed to fit the driver’s characteristics.


## 4. Design of Integrated Chassis Control System

The architecture of the proposed ICC system for AFS and DYC integration with driver behavior identification is shown in Figure9. The control system is mainly consisted of driver identification module, preview optimal curvature driver model, 2-Degree-of-Freedom (DOF) vehicle reference model, and model prediction control (MPC) based controller module. Different kinds of drivers’ characteristics are marked by the three parameters Tp, td, and Th as is shown in Section 3; the parameters Tp, td, and Th are identified online according to the driver’s operation and the state of the vehicle. When the controller is on, the three parameters are used to build the driver model. δfc is defined as predicted front wheel angle based on driver model. δfc of the certain driver could be predicted at next sample time according to the identification result along with vehicle motion state. Based on the predicted δfc, 2-DOF reference model could output the desired yaw rate and side slip angle which are the target of the MPC controller. The MPC controller optimizes online the active control front wheel steering angle Δδf and yaw moment ΔMz.Figure 9
ICC system configurations with driver behavior identification.


### 4.1. Linear 2-DOF Reference Model

The 2-DOF reference vehicle model which considers both accuracy and simplicity is used for target outputs calculation, as shown in Figure10. The side-slip angle β and yaw rate ωr are described in the model.Figure 10
2-DOF reference model.The vehicle state space equation is(19)x˙=A0x+B0u,y=C0x,where(20)
x
=
β
ω
r
T
,
u
=
δ
f
M
Z
T
,
A
0
=
a
11
a
12
a
21
a
22
=
-
c
f
+
c
r
m
·
V
x
-
1
+
c
r
l
r
-
c
f
l
f
m
·
V
x
2
c
r
l
r
-
c
f
l
f
I
z
-
c
r
l
r
2
+
c
f
l
f
2
I
z
·
V
x
,
B
0
=
b
11
b
12
b
21
b
22
=
c
f
m
·
V
x
0
c
f
l
f
I
z
1
I
z
,
C
0
=
c
11
c
12
c
21
c
22
=
1
0
0
1
,where cf, cr are front and rear axle cornering stiffness; m is mass; lf, lr describe the distances from the vehicle CG to the front and rear axle, respectively; Iz is yaw inertia of the vehicle; and δf and Mz are the front wheel steering angle and yaw moment of the vehicle.To guarantee the lateral stability of the vehicle, the yaw rate should be restricted within a stable field. The desired yaw rate could be obtained from the steady-state gain of the yaw rate of the reference model:(21)ωrd=Vx·δflf+lr·1+Vx/vch2,where(22)vch2=cf·cr·lf+lr2mcrlr-cflf.Also, the desired yaw rate should be restricted according to the road friction coefficientμ:(23)ωrd≤μ·gVx.To maintain lateral stability, it is important to sustain driver’s control authority, which can be achieved when the vehicle sideslip angle is small. According to some literatures [24], the desired sideslip angle can be chosen as(24)βd=0.


### 4.2. Design of Integrated Chassis Controller Based on MPC

As a novel methodology with prediction of future states to minimize the deviation between the ideal and actual output of the system, the model predictive control (MPC) shows good stability and robust performance [20, 25–27]. So it draws attentions by the design of integrated chassis control system with consideration of driver behaviors.In this paper, the actual and expected side slip angle and yaw rate are selected as the inputs of the MPC controller, and the outputs are expected active yaw moment and front steering angle. In the controller, the predictive model based on the 2-DOF reference model is built to predict the controlled output firstly. Then, the quadratic programming is developed for the optimal solution.First, the state space function (21) is discretized as(25)xk+1=Axk+Buk,yk=Cxk.And the increment equation could be written as(26)Δxk+1=AΔxk+BΔuk,yk=CΔxk+yk-1,where(27)Δxk=xk-xk-1Δuk=uk-uk-1.A=I+TsA0; B=TsB0; and C=C0, where Ts is the discrete sampling time.Then the predictive controlled output is(28)Yk+1∣k=SΔxk+Eyk+FΔUk,where(29)Yk+1∣k=yk+1∣kyk+2∣k⋮yk+p∣k,ΔUk=ΔukΔuk+1⋮Δuk+c-1,S=CACA2+CA⋮∑i=1pCAiF=CB00⋯0∑i=12CAi-1BCB0⋯0⋮⋮⋮⋯⋮∑i=1CCAi-1B∑i=1C-1CAi-1B⋯CB⋮⋮⋮⋯⋮∑i=1PCAi-1B∑i=1P-1CAi-1B⋯∑i=1P-C+1CAi-1B,E=II⋮I,where p and c are the predictive and control horizons, respectively, and c≤p.The cost function could be described as(30)Jyk,ΔUk=λyYk+1∣k-Ydk+12+λuΔUk2,where λy is the weight coefficient of the output deviation, λu is the weight coefficient of control input increment, and Yd(k+1) is the desired output in the prediction horizon; namely,(31)Ydk+1=βdωrdT.At each time step the following optimization problem is solved:(32)minΔUk⁡Jyk,ΔUk.With the constraints of the inputs, the variation of inputs and the outputs, respectively, is(33)umink+j≤uk+j≤umaxk+j,j=0,1,…,c-1,(34)-Δumink+j≤Δuk+j≤Δumaxk+j,j=0,1,…,c-1,(35)ymink+j≤yk+j≤ymaxk+j,j=0,1,…,p-1,where inequalities (34) limit the control inputs and (35) constrain the changes of the control input, while (36) are constraints on system output variables.Obviously, it is a typical constraint optimization problem and can be transformed into a quadratic programming (QP) problem:(36)minΔUk12ΔUkTWHΔUk+WGTΔUks.t.AconstΔUk≥bconst,where(37)WH=2STλyTλyS+λuTλuWG=-2STλyTλyEpk+1∣kEpk+1∣k=Ydk+1-SΔxk-Eyk,where Aconst and bconst are the constrained matrices.Then, the first sample of the results is used to compute the optimal steering angle and direct yaw moment by the following feedback control law:(38)Δuk=100⋯0010⋯02×2CΔUk,uk=uk-1+Δuk.


## 4.1. Linear 2-DOF Reference Model

The 2-DOF reference vehicle model which considers both accuracy and simplicity is used for target outputs calculation, as shown in Figure10. The side-slip angle β and yaw rate ωr are described in the model.Figure 10
2-DOF reference model.The vehicle state space equation is(19)x˙=A0x+B0u,y=C0x,where(20)
x
=
β
ω
r
T
,
u
=
δ
f
M
Z
T
,
A
0
=
a
11
a
12
a
21
a
22
=
-
c
f
+
c
r
m
·
V
x
-
1
+
c
r
l
r
-
c
f
l
f
m
·
V
x
2
c
r
l
r
-
c
f
l
f
I
z
-
c
r
l
r
2
+
c
f
l
f
2
I
z
·
V
x
,
B
0
=
b
11
b
12
b
21
b
22
=
c
f
m
·
V
x
0
c
f
l
f
I
z
1
I
z
,
C
0
=
c
11
c
12
c
21
c
22
=
1
0
0
1
,where cf, cr are front and rear axle cornering stiffness; m is mass; lf, lr describe the distances from the vehicle CG to the front and rear axle, respectively; Iz is yaw inertia of the vehicle; and δf and Mz are the front wheel steering angle and yaw moment of the vehicle.To guarantee the lateral stability of the vehicle, the yaw rate should be restricted within a stable field. The desired yaw rate could be obtained from the steady-state gain of the yaw rate of the reference model:(21)ωrd=Vx·δflf+lr·1+Vx/vch2,where(22)vch2=cf·cr·lf+lr2mcrlr-cflf.Also, the desired yaw rate should be restricted according to the road friction coefficientμ:(23)ωrd≤μ·gVx.To maintain lateral stability, it is important to sustain driver’s control authority, which can be achieved when the vehicle sideslip angle is small. According to some literatures [24], the desired sideslip angle can be chosen as(24)βd=0.


## 4.2. Design of Integrated Chassis Controller Based on MPC

As a novel methodology with prediction of future states to minimize the deviation between the ideal and actual output of the system, the model predictive control (MPC) shows good stability and robust performance [20, 25–27]. So it draws attentions by the design of integrated chassis control system with consideration of driver behaviors.In this paper, the actual and expected side slip angle and yaw rate are selected as the inputs of the MPC controller, and the outputs are expected active yaw moment and front steering angle. In the controller, the predictive model based on the 2-DOF reference model is built to predict the controlled output firstly. Then, the quadratic programming is developed for the optimal solution.First, the state space function (21) is discretized as(25)xk+1=Axk+Buk,yk=Cxk.And the increment equation could be written as(26)Δxk+1=AΔxk+BΔuk,yk=CΔxk+yk-1,where(27)Δxk=xk-xk-1Δuk=uk-uk-1.A=I+TsA0; B=TsB0; and C=C0, where Ts is the discrete sampling time.Then the predictive controlled output is(28)Yk+1∣k=SΔxk+Eyk+FΔUk,where(29)Yk+1∣k=yk+1∣kyk+2∣k⋮yk+p∣k,ΔUk=ΔukΔuk+1⋮Δuk+c-1,S=CACA2+CA⋮∑i=1pCAiF=CB00⋯0∑i=12CAi-1BCB0⋯0⋮⋮⋮⋯⋮∑i=1CCAi-1B∑i=1C-1CAi-1B⋯CB⋮⋮⋮⋯⋮∑i=1PCAi-1B∑i=1P-1CAi-1B⋯∑i=1P-C+1CAi-1B,E=II⋮I,where p and c are the predictive and control horizons, respectively, and c≤p.The cost function could be described as(30)Jyk,ΔUk=λyYk+1∣k-Ydk+12+λuΔUk2,where λy is the weight coefficient of the output deviation, λu is the weight coefficient of control input increment, and Yd(k+1) is the desired output in the prediction horizon; namely,(31)Ydk+1=βdωrdT.At each time step the following optimization problem is solved:(32)minΔUk⁡Jyk,ΔUk.With the constraints of the inputs, the variation of inputs and the outputs, respectively, is(33)umink+j≤uk+j≤umaxk+j,j=0,1,…,c-1,(34)-Δumink+j≤Δuk+j≤Δumaxk+j,j=0,1,…,c-1,(35)ymink+j≤yk+j≤ymaxk+j,j=0,1,…,p-1,where inequalities (34) limit the control inputs and (35) constrain the changes of the control input, while (36) are constraints on system output variables.Obviously, it is a typical constraint optimization problem and can be transformed into a quadratic programming (QP) problem:(36)minΔUk12ΔUkTWHΔUk+WGTΔUks.t.AconstΔUk≥bconst,where(37)WH=2STλyTλyS+λuTλuWG=-2STλyTλyEpk+1∣kEpk+1∣k=Ydk+1-SΔxk-Eyk,where Aconst and bconst are the constrained matrices.Then, the first sample of the results is used to compute the optimal steering angle and direct yaw moment by the following feedback control law:(38)Δuk=100⋯0010⋯02×2CΔUk,uk=uk-1+Δuk.


## 5. Analysis and Simulation

In order to verify the integrated vehicle chassis control strategy based on driver behavior identification, the cosimulation of MATLAB/Simulink and CarSim is adopted, and the primary parameters of vehicle are shown in Table3.Table 3
Vehicle parameters for simulation.
Parameters
Symbol
Values
Mass (kg)
m
2210
Front axle cornering stiffness (N/rad)
c
f
62800
Rear axle cornering stiffness (N/rad)
c
r
68000
Distance from vehicle CG to the front axle (m)
l
f
1.07
Distance from vehicle CG to the rear axle (m)
l
r
2.23
Vehicle moment of inertia about thez-axis (kgm2)
I
z
4331.6The simulation is carried out under a typical double lane change working condition. The driver model in CarSim is used and its characteristic is identified by the algorithm online. To verify the robust performance of the control strategy when the speed varies, the simulation is started at the longitudinal speed of 60 Km/h and 80 Km/h, respectively.According to the online identification results, the three steering parameters of driver model in CarSim are shown in Table4.Table 4
Identification results of driver behavior characteristic parameters.
T
p
t
d
T
h
1.3886
0.4176
0.1589The three parameters are adopted in the preview optimal curvature driver model built in MATLAB/Simulink. The test results of the control strategy are shown in Figures11-12, in which (a) is the trajectories, (b) and (c) are yaw rate and lateral acceleration respectively, and (d) and (e) are the controlled steering angle and direct yaw moment, respectively. It is seen that both vehicles can track the target path acceptably; moreover, the path-following performance of the vehicle is clearly further improved by the usage of driver behavior identification (DBI) and MPC compared to the vehicle with only MPC. In the meanwhile, both the yaw rate and the lateral acceleration of the vehicle with DBI and MPC track the desired value more closely than vehicle with only MPC, which indicates that the handling and stability are further improved by the introduction of DBI into MPC. Furthermore, both actuators execute reduced output magnitude. It is implied that their energy consumption are optimized compared to the vehicle equipped with only MPC controller; besides, as normally implemented by active braking control, a reduced direct yaw moment control will also lessen the losses of vehicle speed and engine power.Figure 11
Test results at speed of 60 Km/h.
(a)
(b)
(c)
(d)
(e)Figure 12
Test results at speed of 80 Km/h.
(a)
(b)
(c)
(d)
(e)Therefore, on double lane change, the integrated chassis control system based on MPC can effectively ensure the vehicle driving stability. When combined with DBI, the integrated controller further regulates the lateral displacement, yaw rate, and lateral acceleration to a desired region with reduced output of actuators.


## 6. Conclusion

This paper presents an integrated control strategy with driver behavior identification.Firstly, the driver behavior data acquisition system is designed and established based on dSPACE real-time simulation platform, and the driver inputs of different kinds of drivers were collected under the double lane change test condition. As is shown in collected data, among the drivers’ operating signals, the steering wheel angles show the most significant difference among different kinds of drivers.Based on this, preview optimal curvature model is introduced to analyze the driver steering behavior. The genetic algorithm for offline usage and a neural network algorithm for online usage are designed to recognize different patterns of drivers, specifically, to identify the three steering parameters of driver model.Then, an integrated control strategy for active steering angle and direct yaw moment is proposed. The identified driver model is used to predict the steering angle next sample time according to the target path and vehicle state. A linear 2-DOF reference model is adopted to calculate the desired state of the vehicle and model predictive control (MPC) is used to calculate the active control quantity to be exerted on vehicle.Finally, simulations are carried out with the cosimulation of MATLAB/Simulink and CarSim. The results indicate that the proposed method could identify the driver behavior characteristic parameters and enable the control system to adjust its parameters according to the driver behavior identification results. Meanwhile, the integrated control considering the driver’s characteristics further enhanced the vehicle’s handling and stability performance compared to individual MPC control.



---

*Source: 954514-2015-09-03.xml*