# Vehicle Stability Control Strategy Based on Recognition of Driver Turning Intention for Dual-Motor Drive Electric Vehicle

Authors: Shu Wang; Xuan Zhao; Qiang Yu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Vehicle stability control should accurately interpret the driving intention and ensure that the actual state of the vehicle is as consistent as possible with the desired state. This paper proposes a vehicle stability control strategy, which is based on recognition of the driver’s turning intention, for a dual-motor drive electric vehicle. A hybrid model consisting of Gaussian mixture hidden Markov (GHMM) and Generalized Growing and Pruning RBF (GGAP-RBF) neural network is constructed to recognize the driver turning intention in real time. The turning urgency coefficient, which is computed on the basis of the recognition results, is used to establish a modified reference model for vehicle stability control. Then, the upper controller of the vehicle stability control system is constructed using the linear model predictive control theory. The minimum of the quadratic sum of the working load rate of the vehicle tire is taken as the optimization objective. The tire-road adhesion condition, performance of the motor and braking system, and state of the motor are taken as constraints. In addition, a lower controller is established for the vehicle stability control system, with the task of optimizing the allocation of additional yaw moment. Finally, vehicle tests were carried out by conducting double-lane change and single-lane change experiments on a platform for dual-motor drive electric vehicles by using the virtual controller of the A&D5435 hardware. The results show that the stability control system functions appropriately using this control strategy and effectively improves the stability of the vehicle.


---

## Body

## 1. Introduction

Vehicle stability control, which relies on the antilock braking system (ABS) and traction control system (TCS) of a vehicle, plays a significant role in preventing single-vehicle accidents caused by vehicle instability [1]. The mechanical structure, dynamic characteristics, response characteristics, and actuator complexity of a dual-motor drive electric vehicle differ significantly from those of a single-drive vehicle. Therefore, dual-motor drive electric vehicles also require an efficient and stable stability control system.Depending on the system structure, two types of vehicle stability system controllers are used: centralized and hierarchical. Compared with the centralized controller, the hierarchical controller system has superior extendibility and fault tolerance and allows more convenient system maintenance and debugging [2, 3]. Furthermore, the much higher level of integration of the vehicle electrification chassis makes it possible to realize coordination control with other systems such as the X-wire control system and active safety systems, e.g., electric control braking (ECB), electric brake force distribution (EBD), and active front steering (AFS).Based on the current state of the vehicle and the driver’s intention, the upper layer of the hierarchical vehicle stability controller makes decision regarding the additional yaw moment to restore the stable state of the vehicle. Zhang and Wang proposed a vehicle stability control system using generalized proportion integration (PI) control [4]. Wang et al. developed a stability control strategy based on an integral separation PID controller to eliminate integral accumulation of the control system [5]. Zhao et al. designed a vehicle stability control system based on a T-S fuzzy model [6, 7]. Some scholars also adopted robust control to improve the influence of uncertain factors such as system loss, tire load fluctuation, and side wind [8–12]. However, the PID algorithm guarantees neither optimal control of the system nor its control stability. The rules of fuzzy control are established on the basis of a large number of experiments and expert experience and need to be adjusted according to the driving environment, which is both high time and economic consuming. The most notable disadvantage of slide mode control is chattering, which occurs when the system approaches the slide mode surface and which can only be reduced but not eliminated. Model predictive control (MPC) can directly consider actuator constraints and system state constraints during controller design. In addition, MPC can effectively control a multiobjective constrained system in complex engineering systems and has strong portability [13]. Jalali et al. proposed a vehicle stability controller using MPC [14–17]. And, the desired yaw rate and sideslip angle were obtained through the driver steering wheel angle.The lower layer of a vehicle stability controller selects the appropriate actuator and reallocates the additional yaw moment based on certain rules. At present, the most commonly used method relies on two allocation strategies: rule-based and optimization theory-based [18]. The increasing actuators and growing complexity of the electronic control systems of vehicles have resulted in the allocation rules becoming more complex; hence, the allocation accuracy cannot be guaranteed to the same extent as before. Therefore, the optimization theory-based allocation strategy is more suitable.The allocation strategy based on optimization theory usually takes the slip rate of the tire, the minimum of the sum of the working load of the tire, and the minimum of the sum of the longitudinal force of the tire as the optimization objectives. At the same time, the optimization strategy considering the actuator failure and energy economy are gradually put forward. Yin proposed the additional yaw moment and traction force distribution strategy for the purpose of minimizing the tire load [19]. Zhai et al. proposed an average torque distribution strategy, tire-dynamic-load-based torque distribution strategy, and minimum-objective-function-based optimal torque distribution strategy to control the motor driving torque or regenerative braking torque for vehicle stability enhancement [20]. Park et al. proposed a torque distribution strategy that considers the driving economics, driver’s acceleration demand, and tire slip and used fuzzy logic to select a suitable distribution strategy [21]. Kim and Kim chose the regenerative braking and an electrohydraulic brake as the actuators in the vehicle stability controller. Aiming to minimize EHB energy consumption, the braking force distribution ratio of the front and rear axles were optimized [22].According to the above studies, the hierarchical vehicle stability controller obtains the desired driving state by using the driver steering wheel angle. However, in-depth research showed that the use of only the steering wheel angle is inadequate to express the driver intention and obtain the desired driving state vehicle stability controller.The driver turning intention is initiated based on experience depending on the condition of the road, environment, and vehicle state. Then, the driver manipulates the vehicle actuator such that the vehicle responds consistently with their intention. Therefore, the driver’s turning intention reflects their subjective demand for the vehicle state. When drivers take sharp turning operation, such as go through continuous S-turn, hairpin bend, and the other roads with a smaller radius curvature, as well as emergency obstacle avoidance or sudden overtaking, the driver would be considered to intend rapidly maneuvering the vehicle under current road conditions, which also means the driver has certain expectations in terms of the vehicle yaw rate and sideslip. However, for some drivers, because of constraints such as their driving experience, the driving environment, theoretical knowledge, and other factors, this expectation ignores correct judgment of the lateral motion stability of the vehicle. During the process of vehicle stability control, if it is a significant difference between the actual and expected state of the vehicle, the driver would proceed with the turning operation, and the vehicle stability would deteriorate and even cause the driver to distrust the system.Therefore, accurate interpretation and prediction of the driving intention during the process of vehicle stability control and ensuring the actual vehicle state is maintained as closely as possible to the expected vehicle state would greatly improve the stability and driving safety of the vehicle. Driver intention can be obtained using fuzzy reasoning, support vector machine (SVM), artificial neural network (ANN), and hidden Markov model (HMM) [23–27]. Fuzzy reasoning, SVM, and ANN are used to recognize the intentions of driver at a certain moment. However, neural network and fuzzy logic are difficult to deal with temporal-ordered information. Thus, they are mainly used for static recognition problems. HMM, as a kind of dynamic information processing method based on time-series cumulative probability, considers only the state sequence with the maximum log-likelihood and ignores the possibility of small probability events. Therefore, it is difficult to recognize easily confused intentions using HMM.All the above analysis motivated us to develop a method to detect the driver’s turning intention, establish a steering urgency coefficient for sharp turning, and propose a vehicle stability control reference model considering the driver’s turning intention. Based on the abovementioned results, we proposed a stability control strategy for a dual-motor drive electric vehicle. As a result of this research, the desired control target of stability control is no longer only based on the steering wheel angle, and the vehicle state, in accordance with the driving intention, is under control to a certain extent. This paper is organized as follows. The method to identify the driver’s turning intention is presented in Section1. A modified reference model for vehicle stability considering the driver’s turning intention is presented in Section 2. Construction of the vehicle stability control strategy based on the modified reference model is discussed in Section 3. Section 4 compares the vehicle stability control strategy based on the modified reference model constructed in this work with the traditional strategy. Finally, Section 5 presents the conclusions of this study.


## 2. Driver’s Turning Intention Recognition Hybrid Model of GHMM and GGAP-RBF Neural Network

The turning operation is a complex event that continues for a certain period of time. The observation sequence of the turning process is a set of temporal data. The Gaussian hidden Markov model (GHMM) displays a strong modeling ability for dynamic time sequences. However, the model does not take into account overlaps between different classes, and this is a severe limitation. In contrast, generalized growing and pruning RBF (GGAP-RBF), an engineering model, can simulate the thinking mechanism of the human brain, has strong classification and decision-making abilities and can describe uncertain information. Thus, it compensates for the inadequacies of HMM. In addition, the model allows insignificant neurons to be removed in each iterative training cycle to effectively control the growth of the neural network and simplify the structure of a network with large data capacity. However, the ability of GGAP-RBF to describe dynamic sequential processes is not especially strong [28–30]. This led us to construct a hybrid model in the form of a GHMM/GGAP-RBF neural network. Given the advantages of sequential model building and its nonlinear mapping ability, the hybrid model can obtain newly identified information, thereby considerably increasing the accuracy of the classification of classes with slight differences [27]. At the same time, to improve the real-time performance, this work uses the initial stage of the turning operation to identify the driver’s turning intention.The structure of the driver turning intention recognition system based on the GHMM/GGAP-RBF hybrid model is shown in Figure1. The hybrid model includes a lower layer (the GHMM model) and an upper layer (the GGAP-RBF model). The lower layer of the model includes the sharp turning GHMM and the normal turning GHMM. The vehicle speed, the steering wheel angle, the steering wheel angle velocity, and the steering wheel torque are the inputs of the lower layer. The log-likelihood of the sharp turning GHMM and the normal turning GHMM are the outputs of the lower layer. Based on the data of the turning initial stage, the normal turning and sharp turning GHMMs are designed and trained using the Baum–Welch algorithm to calculate the most likely sequence of states. In addition, a forward algorithm is used to calculate the log-likelihood of the GHMMs.Figure 1
Structure of the GHMM/GGAP-RBF hybrid model.In the upper layer, the log-likelihood of the GHMM, the vehicle yaw rate, and the lateral acceleration form a vectorδT11,δT12,…,δT1N. The nonlinear combination of this vector is regarded as the input of the GGAP-RBF neural network, and the nonlinear mapping ability of neural network-based methods is used to recognize the driver actual turning intention.


### 2.1. Establishment of the GHMM

The feature parameters of the GHMM/GGAP-RBF hybrid model have a considerable effect on the accuracy of the recognition of the driver’s turning intention. The reliefF algorithm is used to collect the appropriate parameters. This study employs the steering wheel angle, steering wheel angle velocity, and steering wheel torque as feature parameters to recognize the driver’s steering intention. Based on the initial stage of the turning operation, sharp turning and normal turning GHMMs are established separately. The observation sequence of the GHMM can be described as a multidimensional vector [31]:(1)Ot=at,bt,ct,where at is the steering wheel angle, bt is the steering wheel torque, and ct is the steering wheel angular velocity.The Baum–Welch algorithm is used to optimize the three parameters of the GHMM, which are described asλ=π,A,B, where π is the initial state distribution, A is the state transition probability matrix, and B is a probability density function.The probability density function of the model is(2)biO=∑j=1MωijNO,μij,σij,where NO,μij,σij is the j-dimensional Gauss probability density of state i, O is the observation sequence, μij is the mean of the Gauss function, and σij is the covariance of the Gauss function.Assuming thatεti,j is the probability of the jth Gaussian mixture function in the state observation sequence i at time t, the probability that the Markov chain is in state j at time t+1 is as follows:(3)εti,j=pst=i,st+1=j,OλpOλ=αtiβti∑αtiβti×ωijPO,μij,σ∑1KωijPO,μij,σ,where μij is the mean matrix of the Gaussian mixture function, σ is the mixed covariance matrix, and ωij is the weight of output probabilities of different Gaussian mixture functions. Based on the Gaussian mixture model, the parameter re-estimation is as follows:(4)ωij′=∑t=1Tεti,j∑t=1T∑t=1Kεti,j,μij′=∑t=1Tεti,jOt∑t=1Tεti,j,σij′=∑t=1Tεti,jOt−μijOt−μij′∑t=1Tεti,j.After optimization of the parameters of the GHMM, the matching between the collected data and GHMM is calculated using the forward-backward algorithm.


### 2.2. Establishment of the GGAP-RBF Neural Network

In order to ensure that small probability events can also happen, function (5) is no longer used to recognize the turning intention. Instead, turning intention is described as a function of the log-likelihood of the two GHMMs given test data (loglik1 and loglik2), the vehicle yaw rate ωr, and the lateral acceleration ay expressed as follows:(5)inention=maxloglik1,loglik2,(6)inention=Floglik1,loglik2,ωr,ay.The input parameter of the input layer isloglik1∗, loglik2∗, ωr∗, and ay∗, as shown in equation (7). The output of the layer is the driver turning intention, as shown in equation (8):(7)xi=loglik1∗,loglik2∗,ωr∗,ay∗,(8)fxi=loglikxi=ω0+∑k=1Kωk⋅e−wk−xi2/σk2,where loglik1∗, loglik2∗, ωr∗, and ay∗ represent the deviation standardization of loglik1, loglik2, ωr, and ay; ωk is the center of the RBF of the kth neuron; and σk is the standard deviation of the Gaussian function, which indicates the width of the Gaussian function.After theith training iteration, RAN is used to optimize the growth of neurons. The parameters of a neuron, upon addition of a new neuron, are given as follows [32–35]:(9)ωn=ei=yi−fi−1x,wn=xi,σn=κxi−wir,where ωn is the weight connecting the new neuron to the output neuron; wn is the center of the new neuron; σn is the width of the new neuron; ei=yi−fi−1x is the error of the a priori estimate; κ is an overlap factor that determines the overlap of the responses of the hidden neurons in the input space; and wir is the value of the center of the hidden neuron nearest to xi.The worth of the new neuron is estimated using the conditions in inequality (10). If the conditions are met, the new training data are valuable to the network, and the network performance is improved effectively by the new neuron. As the result, the new neuron is added, and the training data are accepted. If not, the training data and the new neuron are rejected.(10)xi−wir>εi,ei=yi−fxi>emin,where εi and emin are the threshold of the distance and the error, respectively; wir is the value of the center of the hidden neuron nearest to xi; and ei is the error of the priori estimate.Then, insignificant neuron should be judged and removed using a pruning algorithm. The mean squared error of prediction output after thekth neuron is removed from the network in the ith training iteration as follows:(11)Eqk=εk,1,εk,1,…,εk,nq=1n∑i=1nεqk,i1/q=ωkq1n∑i=1nRkqxi1/q,where Rkxi is the Gaussian RBF, n is the training time, q is the norm, and k is the kth hidden neuron.The inputs of the RBF,loglik1∗, loglik2∗, ωr∗, and ay∗, follow the normal sampling distribution Nμ,σ2, respectively. The sample range of the ith training data xi,yi is X, which is divided into J equal small parts Δj. As Δj tends to infinity, the sum over the sample range becomes approximately equal to the integral value.(12)Esigk=limJ⟶∞Eqk=ωkq∏l=1L∫alble−wk−xi,l2/σk2plxdx1/q,where Esig(k) is the significance of the kth neuron to the network, which means the contribution of that neuron to the entire network and P1x is the probability distribution function. If Esig(k) is less than the learning accuracy emin, the neuron is considered to be insignificant and is removed. Otherwise, the neuron is significant and should be retained. Because insignificant neurons can be removed from the GGAP-RBF neural network, the size of the network can be limited to a reasonable range.


### 2.3. Establishment of GHMM/GGAP-RBF Model

The establishment of GHMM/GGAP-RBF model needs offline training by using test data for which is obtained in real vehicle experiment.Driving experience, gender, and personality can affect driver’s decision-making. In order to eliminate the influence of drivers on the test results, three drivers with different driving experiences are selected. According to the analysis of Specification for Design of Municipal Roads (CJJ37-2012), Specification for Design of Interactions on Urban Road (CJJ152-2010), and Vehicle Handling Stability Test Method (GB/T6323-2 2014), the radius of turning the test road is set as 10 m, 25 m, 40 m, and 60 m. The radius of turning the test road is set as 10 m, 25 m, 40 m, and 60 m. And, the tests speed is set as 20 km/h, 30 km/h, and 40 km/h. For distinguishing between the straight driving and turning, straight driving tests are also conducted at 20 km/h, 30 km/h, and 40 km/h. The distribution of test data is shown in Table1.Table 1
The distribution of test data.
20 km/h
30 km/h
40 km/h
Sharp turning
80
80
60
Normal turning
88
88
68
Straight driving
30
21
21The vehicle parameters are shown in Table2. Because of the noise in the sensor data collected by the data acquisition instrument, the data must be preprocessed. T-testing is used to remove abnormal data. The mixed Gaussian clustering method is then used to extract data pertaining to the initial stage of the turning operation as part of the entire turning process. These preprocessed data can be divided into two parts. 75% of the test data is used for model offline training, and the other is used for model online verification. After the offline training, the initial state matrix, the state transition matrix, the weight of each Gaussian function in the GHMM, the mean and covariance of each Gaussian function, and the parameters of the GGAP-RBF could be gotten.Table 2
Vehicle parameters.
Parameter
Value
Vehicle mass
443.9 kg
Axle load distribution
45 : 55
Wheelbase
1200 mm
Axle base
1550 mm
Height of mass center
281.9 mm
Distance from the center of mass to the front axis
854.2 mm
Distance from the center of mass to the rear axis
695.8 mm
Rated power (kW)
32
Peak power (kW)
80
Rated torque (Nm)
80
Peak torque (Nm)
160


## 2.1. Establishment of the GHMM

The feature parameters of the GHMM/GGAP-RBF hybrid model have a considerable effect on the accuracy of the recognition of the driver’s turning intention. The reliefF algorithm is used to collect the appropriate parameters. This study employs the steering wheel angle, steering wheel angle velocity, and steering wheel torque as feature parameters to recognize the driver’s steering intention. Based on the initial stage of the turning operation, sharp turning and normal turning GHMMs are established separately. The observation sequence of the GHMM can be described as a multidimensional vector [31]:(1)Ot=at,bt,ct,where at is the steering wheel angle, bt is the steering wheel torque, and ct is the steering wheel angular velocity.The Baum–Welch algorithm is used to optimize the three parameters of the GHMM, which are described asλ=π,A,B, where π is the initial state distribution, A is the state transition probability matrix, and B is a probability density function.The probability density function of the model is(2)biO=∑j=1MωijNO,μij,σij,where NO,μij,σij is the j-dimensional Gauss probability density of state i, O is the observation sequence, μij is the mean of the Gauss function, and σij is the covariance of the Gauss function.Assuming thatεti,j is the probability of the jth Gaussian mixture function in the state observation sequence i at time t, the probability that the Markov chain is in state j at time t+1 is as follows:(3)εti,j=pst=i,st+1=j,OλpOλ=αtiβti∑αtiβti×ωijPO,μij,σ∑1KωijPO,μij,σ,where μij is the mean matrix of the Gaussian mixture function, σ is the mixed covariance matrix, and ωij is the weight of output probabilities of different Gaussian mixture functions. Based on the Gaussian mixture model, the parameter re-estimation is as follows:(4)ωij′=∑t=1Tεti,j∑t=1T∑t=1Kεti,j,μij′=∑t=1Tεti,jOt∑t=1Tεti,j,σij′=∑t=1Tεti,jOt−μijOt−μij′∑t=1Tεti,j.After optimization of the parameters of the GHMM, the matching between the collected data and GHMM is calculated using the forward-backward algorithm.


## 2.2. Establishment of the GGAP-RBF Neural Network

In order to ensure that small probability events can also happen, function (5) is no longer used to recognize the turning intention. Instead, turning intention is described as a function of the log-likelihood of the two GHMMs given test data (loglik1 and loglik2), the vehicle yaw rate ωr, and the lateral acceleration ay expressed as follows:(5)inention=maxloglik1,loglik2,(6)inention=Floglik1,loglik2,ωr,ay.The input parameter of the input layer isloglik1∗, loglik2∗, ωr∗, and ay∗, as shown in equation (7). The output of the layer is the driver turning intention, as shown in equation (8):(7)xi=loglik1∗,loglik2∗,ωr∗,ay∗,(8)fxi=loglikxi=ω0+∑k=1Kωk⋅e−wk−xi2/σk2,where loglik1∗, loglik2∗, ωr∗, and ay∗ represent the deviation standardization of loglik1, loglik2, ωr, and ay; ωk is the center of the RBF of the kth neuron; and σk is the standard deviation of the Gaussian function, which indicates the width of the Gaussian function.After theith training iteration, RAN is used to optimize the growth of neurons. The parameters of a neuron, upon addition of a new neuron, are given as follows [32–35]:(9)ωn=ei=yi−fi−1x,wn=xi,σn=κxi−wir,where ωn is the weight connecting the new neuron to the output neuron; wn is the center of the new neuron; σn is the width of the new neuron; ei=yi−fi−1x is the error of the a priori estimate; κ is an overlap factor that determines the overlap of the responses of the hidden neurons in the input space; and wir is the value of the center of the hidden neuron nearest to xi.The worth of the new neuron is estimated using the conditions in inequality (10). If the conditions are met, the new training data are valuable to the network, and the network performance is improved effectively by the new neuron. As the result, the new neuron is added, and the training data are accepted. If not, the training data and the new neuron are rejected.(10)xi−wir>εi,ei=yi−fxi>emin,where εi and emin are the threshold of the distance and the error, respectively; wir is the value of the center of the hidden neuron nearest to xi; and ei is the error of the priori estimate.Then, insignificant neuron should be judged and removed using a pruning algorithm. The mean squared error of prediction output after thekth neuron is removed from the network in the ith training iteration as follows:(11)Eqk=εk,1,εk,1,…,εk,nq=1n∑i=1nεqk,i1/q=ωkq1n∑i=1nRkqxi1/q,where Rkxi is the Gaussian RBF, n is the training time, q is the norm, and k is the kth hidden neuron.The inputs of the RBF,loglik1∗, loglik2∗, ωr∗, and ay∗, follow the normal sampling distribution Nμ,σ2, respectively. The sample range of the ith training data xi,yi is X, which is divided into J equal small parts Δj. As Δj tends to infinity, the sum over the sample range becomes approximately equal to the integral value.(12)Esigk=limJ⟶∞Eqk=ωkq∏l=1L∫alble−wk−xi,l2/σk2plxdx1/q,where Esig(k) is the significance of the kth neuron to the network, which means the contribution of that neuron to the entire network and P1x is the probability distribution function. If Esig(k) is less than the learning accuracy emin, the neuron is considered to be insignificant and is removed. Otherwise, the neuron is significant and should be retained. Because insignificant neurons can be removed from the GGAP-RBF neural network, the size of the network can be limited to a reasonable range.


## 2.3. Establishment of GHMM/GGAP-RBF Model

The establishment of GHMM/GGAP-RBF model needs offline training by using test data for which is obtained in real vehicle experiment.Driving experience, gender, and personality can affect driver’s decision-making. In order to eliminate the influence of drivers on the test results, three drivers with different driving experiences are selected. According to the analysis of Specification for Design of Municipal Roads (CJJ37-2012), Specification for Design of Interactions on Urban Road (CJJ152-2010), and Vehicle Handling Stability Test Method (GB/T6323-2 2014), the radius of turning the test road is set as 10 m, 25 m, 40 m, and 60 m. The radius of turning the test road is set as 10 m, 25 m, 40 m, and 60 m. And, the tests speed is set as 20 km/h, 30 km/h, and 40 km/h. For distinguishing between the straight driving and turning, straight driving tests are also conducted at 20 km/h, 30 km/h, and 40 km/h. The distribution of test data is shown in Table1.Table 1
The distribution of test data.
20 km/h
30 km/h
40 km/h
Sharp turning
80
80
60
Normal turning
88
88
68
Straight driving
30
21
21The vehicle parameters are shown in Table2. Because of the noise in the sensor data collected by the data acquisition instrument, the data must be preprocessed. T-testing is used to remove abnormal data. The mixed Gaussian clustering method is then used to extract data pertaining to the initial stage of the turning operation as part of the entire turning process. These preprocessed data can be divided into two parts. 75% of the test data is used for model offline training, and the other is used for model online verification. After the offline training, the initial state matrix, the state transition matrix, the weight of each Gaussian function in the GHMM, the mean and covariance of each Gaussian function, and the parameters of the GGAP-RBF could be gotten.Table 2
Vehicle parameters.
Parameter
Value
Vehicle mass
443.9 kg
Axle load distribution
45 : 55
Wheelbase
1200 mm
Axle base
1550 mm
Height of mass center
281.9 mm
Distance from the center of mass to the front axis
854.2 mm
Distance from the center of mass to the rear axis
695.8 mm
Rated power (kW)
32
Peak power (kW)
80
Rated torque (Nm)
80
Peak torque (Nm)
160


## 3. Reference Model of Vehicle Stability Control System considering the Driver’s Intention

This paper proposes a reference model for vehicle stability control. The model, which takes the driver’s turning intention into consideration, is shown in Figure2. First, the GHMM/GGAP-RBF hybrid model is used to recognize the driver’s turning intention on the basis of data relating to the steering wheel operation: angle, angular velocity, and steering wheel torque. When, at the initial stage of the turning operation, the driver’s intention is identified as sharp turning, the vehicle is made to respond quickly to the desired yaw rate for the current steering wheel angle (i.e., the yaw rate can follow the driver’s steering wheel operation successfully) by correcting the reference yaw rate with the aid of the established steering urgency coefficient. When the steering operation enters the turning keeping stage and the turning reversal stage, vehicle stability is ensured by no longer modifying the reference yaw rate. The reference model is not modified when normal turning is intended.Figure 2
Modified model of reference yaw rate.


### 3.1. Reference Model Based on 2-DOF Linear Vehicle Dynamic Model

The most commonly used vehicle stability control reference model is the 2-DOF linear vehicle dynamic model, which only considers the lateral motion and yaw motion of the vehicle [36]. The state equations are as follows:(13)mv˙−uωr=k1+k2β+1uak1−bk2ωr−k1δf,(14)IZω˙r=ak1−bk2β+1ua2k1+b2k2ωr−ak1δf.Thus, the ideal yaw rate and sideslip angle are as follows:(15)ωrdes=uR=uδfa+b+mu2bk2−ak1/2k1k2a+b,βdes=b−a/2a+bk2mu2a+b+mu2bk2−ak1/2k1k2a+bδf.When the vehicle is driving on a road with a low adhesion coefficient, e.g., when the road surface is wet or covered by snow or sand, the adhesion force allowed by the adhesion conditions between the road surface and tires decreases and cannot produce the high yaw rate required by the vehicle. Therefore, when the 2-DOF linear vehicle dynamic model is adopted as the ideal model, it must be limited by the conditions under which the tires adhere to the road.The upper boundary of the ideal yaw rate is(16)ωrupper_bound=0.85μgu.Therefore, the reference value of the vehicle yaw rate for steady-state steering is(17)ωrref=ωrdes,ωrdes≤ωrupper_bound,ωrupper_boundsgnωrdes,ωrdes>ωrupper_bound.Similarly, the upper bound of the ideal sideslip angle must be specified such that it is not too large. The upper boundary of the sideslip angle is(18)βupper_bound=tan−10.02μg.Therefore, the reference value of the sideslip angle for steady-state steering is(19)βref=βdes,βdes≤βupper_bound,βupper_boundsgnβdes,βdes>βupper_bound.This enables the basic control target of the stability control system, namely, the reference yaw rateωrref and sideslip angle βref, to be obtained.


### 3.2. Reference Model Considering Driver’s Intention

As mentioned in Section1, for the same angle, angular velocity and torque of steering wheel can be used as characteristic parameters to recognize the driver’s turning intention. Therefore, the steering wheel angular velocity and torque are taken as parameters to establish the steering urgency coefficient for sharp turning intention.When the vehicle is undergoing steady-state turning, the relationship between the steering torque and the gradient of lateral acceleration satisfies [37](20)dMLdv2/ρ=mnvlHiLVLl=1g1iLVLFzvnv.Therefore, according to the actual lateral acceleration, the ideal steering wheel torque for steady-state turning can be calculated as(21)MLdes=∫1g1iLVLFzvnvdv2ρ.According to equation (21) and the model of the steering system, the actual yaw rate and speed can be used to calculate the ideal steering wheel angular velocity by differentiating(22)δwdes=L1+Ku2uωr,K=mL2ak2−bk1.δswdes=fδwdes.The deviation between the actual torque and angular velocity and the ideal torque and angular velocity of the steering wheel is used to reflect the urgency of the driver’s sharp turning intention. Finally, the steering urgency coefficient for sharp turning intention is calculated by the following equation:(23)τ=12δ¯sw−δ¯swdes2+M¯L−M¯Ldes2,where δ¯sw, δ¯swdes, M¯L, and M¯Ldes is the normalized actual steering wheel torque, actual angular velocity, ideal steering wheel torque, and ideal angular velocity using the following equation to avoid the impact caused by the different dimensions of each parameter.(24)x¯=x−minmax−min.Whenωrdes≤ωr≤ωrupper_bound, the reference yaw rate is determined by the dynamic characteristics of the vehicle, and the limitation generated by the adhesion condition of the road is not reached. If the driver intention is sharp turning at the same time, the reference yaw rate would be modified as follows:(25)ωrref=ωrdes+Δωr,where Δωr=τωr−ωrdes is the correction of the reference yaw rate considering the driver turning intention and is related to the urgency of the driver turning intention.Whenωrdes≤ωrupper_bound≤ωr, the reference yaw rate is determined by the dynamic characteristics of the vehicle, and the limitation generated by the adhesion condition of the road is also reached. If equation (25) is still adopted for modification, the reference yaw rate may exceed the limitation. Thus, under this condition, the reference yaw rate would be modified as follows:(26)Δωr=τωrupper_bound−ωrdes.Whenωrupper_bound≤ωrdes, the reference yaw rate is determined by the adhesion condition of the road. In order to ensure the driving safety, the reference yaw rate is no longer modified.Thus, for the sharp turning intention, the reference yaw rate is(27)ωrref∗=ωrdes+τωr−ωrdes,whenωrdes≤ωr≤ωrupper_bound,ωrdes+τωrupper_bound−ωrdes,whenωrdes≤ωrupper_bound≤ωr,ωrref,whenωrupper_bound≤ωrdes.When the turning operation reaches the second stage—keeping stage and the third stage—returning stage, the reference yaw rate is no longer modified for the sake of the driving safety. The conversion from the modified model to the reference model causes the reference yaw rate to change suddenly, and this is most likely to influence the stability of the vehicle. In order to make the reference yaw rate smoother, the S-shaped acceleration and deceleration curve is chosen as the transition function, which is shown in Figure3 [38].Figure 3
Switching transition function of reference yaw rate.When the inequality(28)ωrref∗−ωrref>ε,is satisfied, the modified reference yaw rate is(29)ωrref∗t=ωrref∗t0+12Jt2,t0≤t≤t1,ωrref∗t0−Jt12+2Jt1t−12Jt2,t1<t≤t2.The flow of the calculation to obtain the reference yaw rate is shown in Figure4.Figure 4
Flowchart of reference yaw rate.


## 3.1. Reference Model Based on 2-DOF Linear Vehicle Dynamic Model

The most commonly used vehicle stability control reference model is the 2-DOF linear vehicle dynamic model, which only considers the lateral motion and yaw motion of the vehicle [36]. The state equations are as follows:(13)mv˙−uωr=k1+k2β+1uak1−bk2ωr−k1δf,(14)IZω˙r=ak1−bk2β+1ua2k1+b2k2ωr−ak1δf.Thus, the ideal yaw rate and sideslip angle are as follows:(15)ωrdes=uR=uδfa+b+mu2bk2−ak1/2k1k2a+b,βdes=b−a/2a+bk2mu2a+b+mu2bk2−ak1/2k1k2a+bδf.When the vehicle is driving on a road with a low adhesion coefficient, e.g., when the road surface is wet or covered by snow or sand, the adhesion force allowed by the adhesion conditions between the road surface and tires decreases and cannot produce the high yaw rate required by the vehicle. Therefore, when the 2-DOF linear vehicle dynamic model is adopted as the ideal model, it must be limited by the conditions under which the tires adhere to the road.The upper boundary of the ideal yaw rate is(16)ωrupper_bound=0.85μgu.Therefore, the reference value of the vehicle yaw rate for steady-state steering is(17)ωrref=ωrdes,ωrdes≤ωrupper_bound,ωrupper_boundsgnωrdes,ωrdes>ωrupper_bound.Similarly, the upper bound of the ideal sideslip angle must be specified such that it is not too large. The upper boundary of the sideslip angle is(18)βupper_bound=tan−10.02μg.Therefore, the reference value of the sideslip angle for steady-state steering is(19)βref=βdes,βdes≤βupper_bound,βupper_boundsgnβdes,βdes>βupper_bound.This enables the basic control target of the stability control system, namely, the reference yaw rateωrref and sideslip angle βref, to be obtained.


## 3.2. Reference Model Considering Driver’s Intention

As mentioned in Section1, for the same angle, angular velocity and torque of steering wheel can be used as characteristic parameters to recognize the driver’s turning intention. Therefore, the steering wheel angular velocity and torque are taken as parameters to establish the steering urgency coefficient for sharp turning intention.When the vehicle is undergoing steady-state turning, the relationship between the steering torque and the gradient of lateral acceleration satisfies [37](20)dMLdv2/ρ=mnvlHiLVLl=1g1iLVLFzvnv.Therefore, according to the actual lateral acceleration, the ideal steering wheel torque for steady-state turning can be calculated as(21)MLdes=∫1g1iLVLFzvnvdv2ρ.According to equation (21) and the model of the steering system, the actual yaw rate and speed can be used to calculate the ideal steering wheel angular velocity by differentiating(22)δwdes=L1+Ku2uωr,K=mL2ak2−bk1.δswdes=fδwdes.The deviation between the actual torque and angular velocity and the ideal torque and angular velocity of the steering wheel is used to reflect the urgency of the driver’s sharp turning intention. Finally, the steering urgency coefficient for sharp turning intention is calculated by the following equation:(23)τ=12δ¯sw−δ¯swdes2+M¯L−M¯Ldes2,where δ¯sw, δ¯swdes, M¯L, and M¯Ldes is the normalized actual steering wheel torque, actual angular velocity, ideal steering wheel torque, and ideal angular velocity using the following equation to avoid the impact caused by the different dimensions of each parameter.(24)x¯=x−minmax−min.Whenωrdes≤ωr≤ωrupper_bound, the reference yaw rate is determined by the dynamic characteristics of the vehicle, and the limitation generated by the adhesion condition of the road is not reached. If the driver intention is sharp turning at the same time, the reference yaw rate would be modified as follows:(25)ωrref=ωrdes+Δωr,where Δωr=τωr−ωrdes is the correction of the reference yaw rate considering the driver turning intention and is related to the urgency of the driver turning intention.Whenωrdes≤ωrupper_bound≤ωr, the reference yaw rate is determined by the dynamic characteristics of the vehicle, and the limitation generated by the adhesion condition of the road is also reached. If equation (25) is still adopted for modification, the reference yaw rate may exceed the limitation. Thus, under this condition, the reference yaw rate would be modified as follows:(26)Δωr=τωrupper_bound−ωrdes.Whenωrupper_bound≤ωrdes, the reference yaw rate is determined by the adhesion condition of the road. In order to ensure the driving safety, the reference yaw rate is no longer modified.Thus, for the sharp turning intention, the reference yaw rate is(27)ωrref∗=ωrdes+τωr−ωrdes,whenωrdes≤ωr≤ωrupper_bound,ωrdes+τωrupper_bound−ωrdes,whenωrdes≤ωrupper_bound≤ωr,ωrref,whenωrupper_bound≤ωrdes.When the turning operation reaches the second stage—keeping stage and the third stage—returning stage, the reference yaw rate is no longer modified for the sake of the driving safety. The conversion from the modified model to the reference model causes the reference yaw rate to change suddenly, and this is most likely to influence the stability of the vehicle. In order to make the reference yaw rate smoother, the S-shaped acceleration and deceleration curve is chosen as the transition function, which is shown in Figure3 [38].Figure 3
Switching transition function of reference yaw rate.When the inequality(28)ωrref∗−ωrref>ε,is satisfied, the modified reference yaw rate is(29)ωrref∗t=ωrref∗t0+12Jt2,t0≤t≤t1,ωrref∗t0−Jt12+2Jt1t−12Jt2,t1<t≤t2.The flow of the calculation to obtain the reference yaw rate is shown in Figure4.Figure 4
Flowchart of reference yaw rate.


## 4. Vehicle Stability Control Strategy

The proposed stability controller for dual-motor drive electric vehicle adopts a hierarchical structure, which includes the upper layer controller—with a decision layer for the additional yaw moment—and the lower layer controller—with a distribution layer for the additional yaw moment. The structure of this system is shown in Figure5.Figure 5
Structure of stability control system for dual-motor drive electric vehicle.The upper layer controller includes the vehicle state estimator model, the driver turning intention recognition model, the modified stability control reference model and the additional yaw moment decision model. This controller chooses the angle, angular velocity, and torque of the steering wheel as inputs to recognize the driver’s turning intention. The reference yaw rate and reference sideslip angle under steady-state steering are determined by using the 2-DOF linear vehicle dynamics model and by simultaneously considering the adhesion conditions. According to the turning intention, the reference model is modified as the final stability control target. Then, based on the difference between the actual and the reference yaw rate and sideslip angle, using the additional yaw moment decision model, the additional yaw moment that needs to be applied to restore the vehicle to the stable state is determined and used as input for the lower controller.The lower layer controller—the additional yaw moment distribution layer—includes the longitudinal force distribution model and the actuator model. The generation of additional yaw moment requires the longitudinal force of the tire to be controlled, and this should take into consideration the vehicle drive form, performance of the motor and hydraulic brake system, and the road conditions. The optimization algorithm allocates additional yaw moment to get additional torque of the motor and hydraulic braking system which will eventually improve the stability of the vehicle.


### 4.1. Design of Upper Controller for Vehicle Stability Control

MPC is used to make decision of additional yaw moments. And, linear 3-DoF vehicle dynamics model including longitudinal motion, lateral motion, and yaw motion is selected as the predictive model, which is shown in Figure6. Taking the yaw rate and the sideslip angle as the state variables, the vehicle state can be obtained as shown in equation (30)–(32) [39, 40]. For both accuracy and efficiency, the tire cornering stiffness estimator based on recursive least square method with forgetting factor (FFRLS) is used to estimate the tire stiffness which will be used for calculating the tire lateral force in real time. Then, the estimated tire cornering stiffness is applied to the linear dynamic model. This improves the effectiveness of the control system in the nonlinear region of the tire [41].(30)V˙=cosβmcosδflFxrl−sinδflFyfl+cosδfrFxfr−sinδfrFyfr+Fxrl+Fxrr+sinβmsinδflFxrl+cosδflFyfl+sinδfrFxfr+cosδfrFyfr+Fyrr+Fyrl,(31)ω¨r=1IzcosδfrFxrl−cosδflFxflBf2+Fxrr−FxrlBr2+asinδfrFxrl+bsinδflFxfl+Bf2sinδfl+acosδflFyfl+−Bf2sinδfr+acosδfrFyfr+aFyrl+bFyrr,(32)β˙=cosβmVsinδflFxrl+cosδflFyfl+sinδfrFxfr+cosδfrFyfr+Fyrl+Fyrr−sinβmVcosδflFxrl−sinδflFyfl+cosδfrFxfr−sinδfrFyfr+Fxrr+FxrlFyfl−ωr.Figure 6
3-DOF linear vehicle dynamics model.The state equation can be expressed in standard form as follows:(33)X˙=AcX+BucU+BdcD,Yc=CcX,where X=ωr;β;V, Yc=ωr;β, Cc=1,0,0;0,1,0, U is the additional yaw moment, and D=δf.Equation (33) is discretized and converted into incremental form to obtain(34)ΔXk+1=AΔXk+BuΔUk+BdΔDk,Yk=CΔXk+Yk−1,where ΔXk=Xk−Xk−1, ΔUk=Uk−Uk−1, ΔDk=Dk−Dk−1, A=eAcT, Bu=∫0TeAcτBucdτ, C=1,0;0,1, and T is the control period.The objective of vehicle stability control is to enable the actual vehicle yaw rate and sideslip angle to follow the reference value through the action of additional yaw moment. Therefore, at timek, the MPC optimization problem based on the linear model can be described as(35)minJYk,Uk=Y¯k−YrefkQ2+ΔUkR2=∑iNpY¯k+ik−Yrefk+ikQ2+∑iNc−1ΔUk+ikR2.The constraints of control variables, the increment of variables, and output of the model are as follows:(36)Umin≤Uk≤Umax,ΔUmin≤ΔUk≤ΔUmax,Ymin≤Yk≤Ymax.In this model, the predictive horizonNp=5, the control horizon Nc=3, and the reference of the predictive model is Yref=ωrref∗,βref∗.The desired control performance can be adjusted by using the weight matrixQ and R, where Q reflects the accuracy requirements of the system and R reflects the size of the action required by the controller.In this research,(37)Q=100015,R=5000005000008000.As shown in equation (38), ΔUk is the increment of the sequence of control inputs, which is obtained by using the optimization objective and constraints during the sampling time of k, and Y¯k is the control output sequence which is obtained from the prediction model.(38)Y¯k=Yk+1kYk+2kYk+3kYk+4kYk+5k,ΔUk=ΔMkkΔMk+1kΔMk+2k,(39)Yk+1k=CAxk+CBΔUk,Yk+2k=CA2xk+CABΔUk+CBΔUk+1,⋮Yk+5k=∑i=15CAixk+∑i=15CAi−1BΔUk+⋯+∑i=13CAi−1BΔUk+4+∑i=13CAi−1BΔdk+Yk.Finally, only the first element of the optimized solution is applied to the system. Thus, the additional yaw moment is expressed by the following equation:(40)Mk=Mk−1+ΔMk−1k.


### 4.2. Design of the Lower Controller for Vehicle Stability Control

The lower layer controller distributes the longitudinal forces of each wheel according to the output of additional yaw moment decision layer. The distribution is restricted by the adhesion condition between the tire and the road. Excessive additional longitudinal force, which causes longitudinal slip and further deteriorates the vehicle stability, should be avoided. This force is restricted by the performance of the motor and braking system. Excessive additional torque causes the motor and mechanical braking system to overload. The torque is also restricted by the working state of the motor; therefore, a state of system failure, in which the motor cannot provide braking force, should be avoided. Under the above constraints, the distributed longitudinal force needs to meet the demand of additional yaw moment.


#### 4.2.1. Optimization Objective

The longitudinal and lateral forces, the tire can provide, are limited by the vertical load. As the longitudinal force of each tire should be distributed according to the vertical load, wheels with higher adhesion would be expected to play a greater role. Therefore, the optimization objective is to minimize the sum of the square of operating working load rate of each tire. Because of the limitation imposed by practical conditions, the lateral force of wheels cannot be directly controlled. In this study, the longitudinal force of tires is controlled to generate additional yaw moment. Thus, the optimization objective is expressed by the following equation:(41)minJ=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2.


#### 4.2.2. Optimization Constraints

(1) Equality Constraints. The strategy for distributing additional yaw moment should not only minimize the sum of the operating working load rate but also ensure that the longitudinal force meets the requirements of braking and acceleration operations of the driver and that the additional yaw moment meets the requirements of the upper controller.(42)maxq=Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr,ΔM=ΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2.(2) Inequality Constraint. (1)
Adhesion condition constraint: when additional yaw moment is distributed, it is necessary to ensure that the longitudinal and lateral forces of the tire are within the tire adhesion ellipse to avoid longitudinal and lateral vehicle slip.(43)Fxij2+Fyij2≤μijFzij.(2)
Motor and brake system performance constraints: the additional longitudinal force is also limited by the performance of the actuators. The additional longitudinal force of the two rear driving wheels is limited by the performance of the driving motor. The two nondriving front wheels can only provide braking force; thus, the additional longitudinal force of the four wheels is limited by the performance of the mechanical braking system.The braking torque generated by motor is limited by the characteristics of the motor, which cannot exceed the maximum torque limit determined by the power generated at the current speed. The large dynamic fluctuation of the braking torque of the motor and the inverse proportion between the maximum braking torque and speed prevent the motor from generating an adequate amount of reverse electromotive force, such that the braking force generated by the motor is ineffective. Due to the large dynamic fluctuation of motor braking torque and the inverse proportion between the maximum braking torque and speed, the reverse electromotive force generated by the motor is too little when the motor is at low speed. The motor cannot generate effective braking force. Therefore, 500 r/min is set as the speed threshold. When the speed is lower than threshold, the system no longer uses motor braking. Therefore, the motor braking torque is required meet the limitation provided in the following equation: [42].(44)Tdmax=0,0<n≤500,Te,500<n≤nN,9550Pe/ηn,nN<n≤nmax.At the same time, control signal failure, power converter failure, insulation failure, and other types of failures of the motor may occur due to the design defects, the service environment and service life of the motor. Thus, the driving motor is no longer suitable to provide longitudinal force, so the motor failure coefficientτij is introduced.(45)τij=0,motor is failed,1,motor is working properly.Therefore, considering the above factors comprehensively, additional longitudinal force of wheels would have to meet the constraints.(46)Tbflmaxr−Fxfl≤ΔFxfl≤−Fxfl,Tbfrmaxr−Fxfr≤ΔFxfr≤−Fxfr,Tbrlmaxr+Tdmaxiτrlr−Fxrl≤ΔFxrl≤Tdmaxiτrlr−Fxrl,Tbrrmaxr+Tdmaxiτrrr−Fxrr≤ΔFxrr≤Tdmaxiτrrr−Fxrr.Therefore, the optimal distribution of additional yawing moment can be expressed as follows:(47)minJΔFxfl,ΔFxfr,ΔFxrl,ΔFxrr=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2,s.t.Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr=maxqΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2=ΔMFxij2+Fyij2≤μijFzijij=fl,fr,rl,rrTbflmaxr−Fxfl≤ΔFxfl≤−FxflTbfrmaxr−Fxfr≤ΔFxfr≤−FxfrTbrlmaxr−Fxrl≤ΔFxrl≤Tdmaxir−FxrlTbrrmaxr−Fxrr≤ΔFxrr≤Tdmaxir−Fxrr.


## 4.1. Design of Upper Controller for Vehicle Stability Control

MPC is used to make decision of additional yaw moments. And, linear 3-DoF vehicle dynamics model including longitudinal motion, lateral motion, and yaw motion is selected as the predictive model, which is shown in Figure6. Taking the yaw rate and the sideslip angle as the state variables, the vehicle state can be obtained as shown in equation (30)–(32) [39, 40]. For both accuracy and efficiency, the tire cornering stiffness estimator based on recursive least square method with forgetting factor (FFRLS) is used to estimate the tire stiffness which will be used for calculating the tire lateral force in real time. Then, the estimated tire cornering stiffness is applied to the linear dynamic model. This improves the effectiveness of the control system in the nonlinear region of the tire [41].(30)V˙=cosβmcosδflFxrl−sinδflFyfl+cosδfrFxfr−sinδfrFyfr+Fxrl+Fxrr+sinβmsinδflFxrl+cosδflFyfl+sinδfrFxfr+cosδfrFyfr+Fyrr+Fyrl,(31)ω¨r=1IzcosδfrFxrl−cosδflFxflBf2+Fxrr−FxrlBr2+asinδfrFxrl+bsinδflFxfl+Bf2sinδfl+acosδflFyfl+−Bf2sinδfr+acosδfrFyfr+aFyrl+bFyrr,(32)β˙=cosβmVsinδflFxrl+cosδflFyfl+sinδfrFxfr+cosδfrFyfr+Fyrl+Fyrr−sinβmVcosδflFxrl−sinδflFyfl+cosδfrFxfr−sinδfrFyfr+Fxrr+FxrlFyfl−ωr.Figure 6
3-DOF linear vehicle dynamics model.The state equation can be expressed in standard form as follows:(33)X˙=AcX+BucU+BdcD,Yc=CcX,where X=ωr;β;V, Yc=ωr;β, Cc=1,0,0;0,1,0, U is the additional yaw moment, and D=δf.Equation (33) is discretized and converted into incremental form to obtain(34)ΔXk+1=AΔXk+BuΔUk+BdΔDk,Yk=CΔXk+Yk−1,where ΔXk=Xk−Xk−1, ΔUk=Uk−Uk−1, ΔDk=Dk−Dk−1, A=eAcT, Bu=∫0TeAcτBucdτ, C=1,0;0,1, and T is the control period.The objective of vehicle stability control is to enable the actual vehicle yaw rate and sideslip angle to follow the reference value through the action of additional yaw moment. Therefore, at timek, the MPC optimization problem based on the linear model can be described as(35)minJYk,Uk=Y¯k−YrefkQ2+ΔUkR2=∑iNpY¯k+ik−Yrefk+ikQ2+∑iNc−1ΔUk+ikR2.The constraints of control variables, the increment of variables, and output of the model are as follows:(36)Umin≤Uk≤Umax,ΔUmin≤ΔUk≤ΔUmax,Ymin≤Yk≤Ymax.In this model, the predictive horizonNp=5, the control horizon Nc=3, and the reference of the predictive model is Yref=ωrref∗,βref∗.The desired control performance can be adjusted by using the weight matrixQ and R, where Q reflects the accuracy requirements of the system and R reflects the size of the action required by the controller.In this research,(37)Q=100015,R=5000005000008000.As shown in equation (38), ΔUk is the increment of the sequence of control inputs, which is obtained by using the optimization objective and constraints during the sampling time of k, and Y¯k is the control output sequence which is obtained from the prediction model.(38)Y¯k=Yk+1kYk+2kYk+3kYk+4kYk+5k,ΔUk=ΔMkkΔMk+1kΔMk+2k,(39)Yk+1k=CAxk+CBΔUk,Yk+2k=CA2xk+CABΔUk+CBΔUk+1,⋮Yk+5k=∑i=15CAixk+∑i=15CAi−1BΔUk+⋯+∑i=13CAi−1BΔUk+4+∑i=13CAi−1BΔdk+Yk.Finally, only the first element of the optimized solution is applied to the system. Thus, the additional yaw moment is expressed by the following equation:(40)Mk=Mk−1+ΔMk−1k.


## 4.2. Design of the Lower Controller for Vehicle Stability Control

The lower layer controller distributes the longitudinal forces of each wheel according to the output of additional yaw moment decision layer. The distribution is restricted by the adhesion condition between the tire and the road. Excessive additional longitudinal force, which causes longitudinal slip and further deteriorates the vehicle stability, should be avoided. This force is restricted by the performance of the motor and braking system. Excessive additional torque causes the motor and mechanical braking system to overload. The torque is also restricted by the working state of the motor; therefore, a state of system failure, in which the motor cannot provide braking force, should be avoided. Under the above constraints, the distributed longitudinal force needs to meet the demand of additional yaw moment.


### 4.2.1. Optimization Objective

The longitudinal and lateral forces, the tire can provide, are limited by the vertical load. As the longitudinal force of each tire should be distributed according to the vertical load, wheels with higher adhesion would be expected to play a greater role. Therefore, the optimization objective is to minimize the sum of the square of operating working load rate of each tire. Because of the limitation imposed by practical conditions, the lateral force of wheels cannot be directly controlled. In this study, the longitudinal force of tires is controlled to generate additional yaw moment. Thus, the optimization objective is expressed by the following equation:(41)minJ=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2.


### 4.2.2. Optimization Constraints

(1) Equality Constraints. The strategy for distributing additional yaw moment should not only minimize the sum of the operating working load rate but also ensure that the longitudinal force meets the requirements of braking and acceleration operations of the driver and that the additional yaw moment meets the requirements of the upper controller.(42)maxq=Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr,ΔM=ΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2.(2) Inequality Constraint. (1)
Adhesion condition constraint: when additional yaw moment is distributed, it is necessary to ensure that the longitudinal and lateral forces of the tire are within the tire adhesion ellipse to avoid longitudinal and lateral vehicle slip.(43)Fxij2+Fyij2≤μijFzij.(2)
Motor and brake system performance constraints: the additional longitudinal force is also limited by the performance of the actuators. The additional longitudinal force of the two rear driving wheels is limited by the performance of the driving motor. The two nondriving front wheels can only provide braking force; thus, the additional longitudinal force of the four wheels is limited by the performance of the mechanical braking system.The braking torque generated by motor is limited by the characteristics of the motor, which cannot exceed the maximum torque limit determined by the power generated at the current speed. The large dynamic fluctuation of the braking torque of the motor and the inverse proportion between the maximum braking torque and speed prevent the motor from generating an adequate amount of reverse electromotive force, such that the braking force generated by the motor is ineffective. Due to the large dynamic fluctuation of motor braking torque and the inverse proportion between the maximum braking torque and speed, the reverse electromotive force generated by the motor is too little when the motor is at low speed. The motor cannot generate effective braking force. Therefore, 500 r/min is set as the speed threshold. When the speed is lower than threshold, the system no longer uses motor braking. Therefore, the motor braking torque is required meet the limitation provided in the following equation: [42].(44)Tdmax=0,0<n≤500,Te,500<n≤nN,9550Pe/ηn,nN<n≤nmax.At the same time, control signal failure, power converter failure, insulation failure, and other types of failures of the motor may occur due to the design defects, the service environment and service life of the motor. Thus, the driving motor is no longer suitable to provide longitudinal force, so the motor failure coefficientτij is introduced.(45)τij=0,motor is failed,1,motor is working properly.Therefore, considering the above factors comprehensively, additional longitudinal force of wheels would have to meet the constraints.(46)Tbflmaxr−Fxfl≤ΔFxfl≤−Fxfl,Tbfrmaxr−Fxfr≤ΔFxfr≤−Fxfr,Tbrlmaxr+Tdmaxiτrlr−Fxrl≤ΔFxrl≤Tdmaxiτrlr−Fxrl,Tbrrmaxr+Tdmaxiτrrr−Fxrr≤ΔFxrr≤Tdmaxiτrrr−Fxrr.Therefore, the optimal distribution of additional yawing moment can be expressed as follows:(47)minJΔFxfl,ΔFxfr,ΔFxrl,ΔFxrr=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2,s.t.Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr=maxqΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2=ΔMFxij2+Fyij2≤μijFzijij=fl,fr,rl,rrTbflmaxr−Fxfl≤ΔFxfl≤−FxflTbfrmaxr−Fxfr≤ΔFxfr≤−FxfrTbrlmaxr−Fxrl≤ΔFxrl≤Tdmaxir−FxrlTbrrmaxr−Fxrr≤ΔFxrr≤Tdmaxir−Fxrr.


## 4.2.1. Optimization Objective

The longitudinal and lateral forces, the tire can provide, are limited by the vertical load. As the longitudinal force of each tire should be distributed according to the vertical load, wheels with higher adhesion would be expected to play a greater role. Therefore, the optimization objective is to minimize the sum of the square of operating working load rate of each tire. Because of the limitation imposed by practical conditions, the lateral force of wheels cannot be directly controlled. In this study, the longitudinal force of tires is controlled to generate additional yaw moment. Thus, the optimization objective is expressed by the following equation:(41)minJ=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2.


## 4.2.2. Optimization Constraints

(1) Equality Constraints. The strategy for distributing additional yaw moment should not only minimize the sum of the operating working load rate but also ensure that the longitudinal force meets the requirements of braking and acceleration operations of the driver and that the additional yaw moment meets the requirements of the upper controller.(42)maxq=Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr,ΔM=ΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2.(2) Inequality Constraint. (1)
Adhesion condition constraint: when additional yaw moment is distributed, it is necessary to ensure that the longitudinal and lateral forces of the tire are within the tire adhesion ellipse to avoid longitudinal and lateral vehicle slip.(43)Fxij2+Fyij2≤μijFzij.(2)
Motor and brake system performance constraints: the additional longitudinal force is also limited by the performance of the actuators. The additional longitudinal force of the two rear driving wheels is limited by the performance of the driving motor. The two nondriving front wheels can only provide braking force; thus, the additional longitudinal force of the four wheels is limited by the performance of the mechanical braking system.The braking torque generated by motor is limited by the characteristics of the motor, which cannot exceed the maximum torque limit determined by the power generated at the current speed. The large dynamic fluctuation of the braking torque of the motor and the inverse proportion between the maximum braking torque and speed prevent the motor from generating an adequate amount of reverse electromotive force, such that the braking force generated by the motor is ineffective. Due to the large dynamic fluctuation of motor braking torque and the inverse proportion between the maximum braking torque and speed, the reverse electromotive force generated by the motor is too little when the motor is at low speed. The motor cannot generate effective braking force. Therefore, 500 r/min is set as the speed threshold. When the speed is lower than threshold, the system no longer uses motor braking. Therefore, the motor braking torque is required meet the limitation provided in the following equation: [42].(44)Tdmax=0,0<n≤500,Te,500<n≤nN,9550Pe/ηn,nN<n≤nmax.At the same time, control signal failure, power converter failure, insulation failure, and other types of failures of the motor may occur due to the design defects, the service environment and service life of the motor. Thus, the driving motor is no longer suitable to provide longitudinal force, so the motor failure coefficientτij is introduced.(45)τij=0,motor is failed,1,motor is working properly.Therefore, considering the above factors comprehensively, additional longitudinal force of wheels would have to meet the constraints.(46)Tbflmaxr−Fxfl≤ΔFxfl≤−Fxfl,Tbfrmaxr−Fxfr≤ΔFxfr≤−Fxfr,Tbrlmaxr+Tdmaxiτrlr−Fxrl≤ΔFxrl≤Tdmaxiτrlr−Fxrl,Tbrrmaxr+Tdmaxiτrrr−Fxrr≤ΔFxrr≤Tdmaxiτrrr−Fxrr.Therefore, the optimal distribution of additional yawing moment can be expressed as follows:(47)minJΔFxfl,ΔFxfr,ΔFxrl,ΔFxrr=CflFxfl+ΔFxfl2μflFzfl2+CfrFxfr+ΔFxfr2μfrFzfr2+CrlFxrl+ΔFxrl2μrlFzrl2+CrrFxrr+ΔFxrr2μrrFzrr2,s.t.Fxfl+ΔFxflcosδwfl+Fxfr+ΔFxfrcosδwfr+Fxrl+ΔFxrl+Fxrr+ΔFxrr=maxqΔFxflcosδwfl2−ΔFxfrcosδwfr2+ΔFxrl2−ΔFxrr2=ΔMFxij2+Fyij2≤μijFzijij=fl,fr,rl,rrTbflmaxr−Fxfl≤ΔFxfl≤−FxflTbfrmaxr−Fxfr≤ΔFxfr≤−FxfrTbrlmaxr−Fxrl≤ΔFxrl≤Tdmaxir−FxrlTbrrmaxr−Fxrr≤ΔFxrr≤Tdmaxir−Fxrr.


## 5. Test Verification

Considering the practical difficulties associated with controller development, i.e., a long development cycle and high cost, we built a test platform for a dual-motor drive electric vehicle based on the A&D5435 semiphysical simulation system and rapid prototyping technology. The tests were carried out under both double-lane and single-lane change conditions, respectively, to verify the feasibility and accuracy of the proposed stability control strategy.


### 5.1. Dual-Motor Drive EV Test Platform Based on the A&D5435 Semiphysical Simulation System

In this study, a dual-motor drive electric vehicle test platform was built based on the A&D5435 hardware in the loop simulation system and a dual-motor drive electric test vehicle. The stability control system proposed in this paper was tested and verified by using this platform in which the A&D5435 replaces the vehicle control unit. The input signals of the controller include: angle, angular velocity and torque of the steering wheel, accelerator pedal opening, brake pedal opening, vehicle speed, wheel speed, motor torque, and motor power. The output signals include the following: the motor drive torque and the brake torque of both the motor and the hydraulic braking system. The input signal mentioned above can be obtained by using the following sensors. The angle, angular velocity, and torque of the steering wheel can be collected using the steering wheel torque and angle sensor manufactured by SensorWay, which is mounted on the steering column. A hall noncontact speed sensor is used to record the speed of the four wheels. A Passat B5 accelerator pedal sensor, which has two channels, is used to measure the accelerator pedal opening. Likewise, the brake pedal opening is measured by the brake pedal sensor. The longitudinal acceleration, lateral acceleration, and yaw rate of the vehicle are acquired by a three-axis gyroscope. The motor speed, torque, and power are obtained from the CAN signal of the motor controller. The vehicle parameters are shown in Table2. The dual-motor drive electric vehicle test platform is shown in Figure 7.Figure 7
The dual-motor drive electric vehicle test platform.


### 5.2. Test Verification: Road Tests

The stability control strategy of dual-motor drive EV proposed in this paper was verified by carrying out vehicle road tests involving both double-lane and single-lane change, respectively. At the same time, in order to demonstrate the effect of the proposed control strategy, the slide mode control (SMC) with the vehicle reference model considering the driver’s intention is chosen as the comparative controller.The vehicle road test was conducted on the stability performance test square of the Chang’an University vehicle comprehensive performance proving ground, as shown in Figure8.Figure 8
Chang’an University vehicle comprehensive performance proving ground.


#### 5.2.1. Double-Lane Change Condition

The test was conducted under ISO 3888-1 which specifies the standard double-lane change condition [43]. The target driving path is shown in Figure 9. The adhesion coefficient of the test road was 0.6.Figure 9
Schematic diagram for double-lane change test (d1=1.1L+0.25; d2=1.3L+0.25; L is the vehicle width).First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 70 km/h, during which time the steering wheel was not to be operated. The driver could then drive freely by negotiating the cone track.The identification result of driver turning intention is shown in Figure10. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 11–14. The comparison with the control experiment is presented in Table 3.Figure 10
Driver turning intention recognition result.Figure 11
Yaw rate with control and without control.Figure 12
Sideslip angle with control and without control.Figure 13
Lateral acceleration with control and without control.Figure 14
Vehicle track with control.Table 3
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−35.517°/s
9.947°
5.733 m/s2
With control
−18.744°/s
4.391°
4.638 m/s2
With SMC
−18.953°/s
4.568°
4.895 m/s2
Reference
−18.660°/s
4.036°
—
Decrease
47.22%
55.85%
19.10%Figure11 shows that the maximum yaw rate without control is −35.517°/s, whereas the reference yaw rate is −18.66°/s, indicative of poor vehicle stability. At 7.3 s, the current initial stage of sharp turning intention is identified, whereupon the reference yaw rate is modified and the stability control system activated. After 8.45 s, the model identified that the current intention is normal turning, thus the reference yaw rate was no longer modified. As a result of the control provided by the vehicle stability system, the maximum yaw rate decreased by 47.22%, to 18.744°/s, the maximum sideslip angle decreased by 55.26%, to −4.391°, and the maximum lateral acceleration decreased by 19.1%, to 4.638 m/s2. Although dangerous conditions such as sideslip and spin without control did not arise, the actual sideslip angle and yaw rate are much higher than the reference. The stability system clearly reduces the sideslip angle and yaw rate to effectively improve the vehicle stability. Thus, the vehicle stability control system enables the vehicle to work much more smoothly. The actual driving trajectory is shown in Figure 14. Meanwhile, it can be seen from Table 3 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


#### 5.2.2. Single-Lane Change Condition

The test was conducted by performing a single-lane change. The target driving path is shown in Figure15. The adhesion coefficient of the test road was 0.4. First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 40 km/h, during which time the steering wheel was not to be operated. The driver was then allowed to drive freely by following the cone track.Figure 15
Schematic diagram for single-lane change test (h=3.5m; b=1.2L+0.25m; L is the vehicle width; d1=50m; l2=100m; l3=15u).The identification result of driver turning intention is shown in Figure16. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 17–20. The comparison with the control effect is presented in Table 4.Figure 16
Driver turning intention recognition result.Figure 17
Yaw rate with control and without control.Figure 18
Sideslip angle with control and without control.Figure 19
Lateral acceleration with control and without control.Figure 20
Vehicle track with control.Table 4
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−62.779°/s
70.740°
5.761 m/s2
With control
21.913°/s
5.157°
2.705 m/s2
With SMC
21.926°/s
4.566°
3.124 m/s2
Reference
21.76°/s
5.208°
—
Decrease
65.1%
92.7%
53.04%Figures17 and 18 show that the vehicle slips and spins after 7.5 s without stability control. The maximum yaw velocity reached −62.779°/s, and the maximum sideslip angle reached 70.740°, which are much more than the reference model. Thus, the vehicle stability control system was activated. Furthermore, the model identified that the current intention was normal turning; hence, the reference yaw rate was no longer modified. The vehicle stability system reduced the maximum yaw rate by 65.1%, and the maximum sideslip angle by 92.7%. The stability system therefore improved the vehicle stability effectively. The vehicle stability control system enables the vehicle to run much more smoothly in comparison, preventing dangerous conditions such as sideslip and spin. The actual driving trajectory is shown in Figure 20. Meanwhile, it can be seen from Table 4 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


## 5.1. Dual-Motor Drive EV Test Platform Based on the A&D5435 Semiphysical Simulation System

In this study, a dual-motor drive electric vehicle test platform was built based on the A&D5435 hardware in the loop simulation system and a dual-motor drive electric test vehicle. The stability control system proposed in this paper was tested and verified by using this platform in which the A&D5435 replaces the vehicle control unit. The input signals of the controller include: angle, angular velocity and torque of the steering wheel, accelerator pedal opening, brake pedal opening, vehicle speed, wheel speed, motor torque, and motor power. The output signals include the following: the motor drive torque and the brake torque of both the motor and the hydraulic braking system. The input signal mentioned above can be obtained by using the following sensors. The angle, angular velocity, and torque of the steering wheel can be collected using the steering wheel torque and angle sensor manufactured by SensorWay, which is mounted on the steering column. A hall noncontact speed sensor is used to record the speed of the four wheels. A Passat B5 accelerator pedal sensor, which has two channels, is used to measure the accelerator pedal opening. Likewise, the brake pedal opening is measured by the brake pedal sensor. The longitudinal acceleration, lateral acceleration, and yaw rate of the vehicle are acquired by a three-axis gyroscope. The motor speed, torque, and power are obtained from the CAN signal of the motor controller. The vehicle parameters are shown in Table2. The dual-motor drive electric vehicle test platform is shown in Figure 7.Figure 7
The dual-motor drive electric vehicle test platform.


## 5.2. Test Verification: Road Tests

The stability control strategy of dual-motor drive EV proposed in this paper was verified by carrying out vehicle road tests involving both double-lane and single-lane change, respectively. At the same time, in order to demonstrate the effect of the proposed control strategy, the slide mode control (SMC) with the vehicle reference model considering the driver’s intention is chosen as the comparative controller.The vehicle road test was conducted on the stability performance test square of the Chang’an University vehicle comprehensive performance proving ground, as shown in Figure8.Figure 8
Chang’an University vehicle comprehensive performance proving ground.


### 5.2.1. Double-Lane Change Condition

The test was conducted under ISO 3888-1 which specifies the standard double-lane change condition [43]. The target driving path is shown in Figure 9. The adhesion coefficient of the test road was 0.6.Figure 9
Schematic diagram for double-lane change test (d1=1.1L+0.25; d2=1.3L+0.25; L is the vehicle width).First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 70 km/h, during which time the steering wheel was not to be operated. The driver could then drive freely by negotiating the cone track.The identification result of driver turning intention is shown in Figure10. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 11–14. The comparison with the control experiment is presented in Table 3.Figure 10
Driver turning intention recognition result.Figure 11
Yaw rate with control and without control.Figure 12
Sideslip angle with control and without control.Figure 13
Lateral acceleration with control and without control.Figure 14
Vehicle track with control.Table 3
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−35.517°/s
9.947°
5.733 m/s2
With control
−18.744°/s
4.391°
4.638 m/s2
With SMC
−18.953°/s
4.568°
4.895 m/s2
Reference
−18.660°/s
4.036°
—
Decrease
47.22%
55.85%
19.10%Figure11 shows that the maximum yaw rate without control is −35.517°/s, whereas the reference yaw rate is −18.66°/s, indicative of poor vehicle stability. At 7.3 s, the current initial stage of sharp turning intention is identified, whereupon the reference yaw rate is modified and the stability control system activated. After 8.45 s, the model identified that the current intention is normal turning, thus the reference yaw rate was no longer modified. As a result of the control provided by the vehicle stability system, the maximum yaw rate decreased by 47.22%, to 18.744°/s, the maximum sideslip angle decreased by 55.26%, to −4.391°, and the maximum lateral acceleration decreased by 19.1%, to 4.638 m/s2. Although dangerous conditions such as sideslip and spin without control did not arise, the actual sideslip angle and yaw rate are much higher than the reference. The stability system clearly reduces the sideslip angle and yaw rate to effectively improve the vehicle stability. Thus, the vehicle stability control system enables the vehicle to work much more smoothly. The actual driving trajectory is shown in Figure 14. Meanwhile, it can be seen from Table 3 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


### 5.2.2. Single-Lane Change Condition

The test was conducted by performing a single-lane change. The target driving path is shown in Figure15. The adhesion coefficient of the test road was 0.4. First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 40 km/h, during which time the steering wheel was not to be operated. The driver was then allowed to drive freely by following the cone track.Figure 15
Schematic diagram for single-lane change test (h=3.5m; b=1.2L+0.25m; L is the vehicle width; d1=50m; l2=100m; l3=15u).The identification result of driver turning intention is shown in Figure16. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 17–20. The comparison with the control effect is presented in Table 4.Figure 16
Driver turning intention recognition result.Figure 17
Yaw rate with control and without control.Figure 18
Sideslip angle with control and without control.Figure 19
Lateral acceleration with control and without control.Figure 20
Vehicle track with control.Table 4
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−62.779°/s
70.740°
5.761 m/s2
With control
21.913°/s
5.157°
2.705 m/s2
With SMC
21.926°/s
4.566°
3.124 m/s2
Reference
21.76°/s
5.208°
—
Decrease
65.1%
92.7%
53.04%Figures17 and 18 show that the vehicle slips and spins after 7.5 s without stability control. The maximum yaw velocity reached −62.779°/s, and the maximum sideslip angle reached 70.740°, which are much more than the reference model. Thus, the vehicle stability control system was activated. Furthermore, the model identified that the current intention was normal turning; hence, the reference yaw rate was no longer modified. The vehicle stability system reduced the maximum yaw rate by 65.1%, and the maximum sideslip angle by 92.7%. The stability system therefore improved the vehicle stability effectively. The vehicle stability control system enables the vehicle to run much more smoothly in comparison, preventing dangerous conditions such as sideslip and spin. The actual driving trajectory is shown in Figure 20. Meanwhile, it can be seen from Table 4 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


## 5.2.1. Double-Lane Change Condition

The test was conducted under ISO 3888-1 which specifies the standard double-lane change condition [43]. The target driving path is shown in Figure 9. The adhesion coefficient of the test road was 0.6.Figure 9
Schematic diagram for double-lane change test (d1=1.1L+0.25; d2=1.3L+0.25; L is the vehicle width).First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 70 km/h, during which time the steering wheel was not to be operated. The driver could then drive freely by negotiating the cone track.The identification result of driver turning intention is shown in Figure10. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 11–14. The comparison with the control experiment is presented in Table 3.Figure 10
Driver turning intention recognition result.Figure 11
Yaw rate with control and without control.Figure 12
Sideslip angle with control and without control.Figure 13
Lateral acceleration with control and without control.Figure 14
Vehicle track with control.Table 3
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−35.517°/s
9.947°
5.733 m/s2
With control
−18.744°/s
4.391°
4.638 m/s2
With SMC
−18.953°/s
4.568°
4.895 m/s2
Reference
−18.660°/s
4.036°
—
Decrease
47.22%
55.85%
19.10%Figure11 shows that the maximum yaw rate without control is −35.517°/s, whereas the reference yaw rate is −18.66°/s, indicative of poor vehicle stability. At 7.3 s, the current initial stage of sharp turning intention is identified, whereupon the reference yaw rate is modified and the stability control system activated. After 8.45 s, the model identified that the current intention is normal turning, thus the reference yaw rate was no longer modified. As a result of the control provided by the vehicle stability system, the maximum yaw rate decreased by 47.22%, to 18.744°/s, the maximum sideslip angle decreased by 55.26%, to −4.391°, and the maximum lateral acceleration decreased by 19.1%, to 4.638 m/s2. Although dangerous conditions such as sideslip and spin without control did not arise, the actual sideslip angle and yaw rate are much higher than the reference. The stability system clearly reduces the sideslip angle and yaw rate to effectively improve the vehicle stability. Thus, the vehicle stability control system enables the vehicle to work much more smoothly. The actual driving trajectory is shown in Figure 14. Meanwhile, it can be seen from Table 3 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


## 5.2.2. Single-Lane Change Condition

The test was conducted by performing a single-lane change. The target driving path is shown in Figure15. The adhesion coefficient of the test road was 0.4. First, the driver was required to keep the vehicle parallel with the road and keep the steering wheel facing forward. Subsequently, the driver controlled the vehicle to rapidly proceed ahead, such that the vehicle speed before entering the target path reached 40 km/h, during which time the steering wheel was not to be operated. The driver was then allowed to drive freely by following the cone track.Figure 15
Schematic diagram for single-lane change test (h=3.5m; b=1.2L+0.25m; L is the vehicle width; d1=50m; l2=100m; l3=15u).The identification result of driver turning intention is shown in Figure16. The yaw rate, sideslip angle, lateral acceleration, and vehicle trajectory are shown in Figures 17–20. The comparison with the control effect is presented in Table 4.Figure 16
Driver turning intention recognition result.Figure 17
Yaw rate with control and without control.Figure 18
Sideslip angle with control and without control.Figure 19
Lateral acceleration with control and without control.Figure 20
Vehicle track with control.Table 4
Control effect comparison.
Yaw rate
Sideslip angle
Lateral acceleration
Without control
−62.779°/s
70.740°
5.761 m/s2
With control
21.913°/s
5.157°
2.705 m/s2
With SMC
21.926°/s
4.566°
3.124 m/s2
Reference
21.76°/s
5.208°
—
Decrease
65.1%
92.7%
53.04%Figures17 and 18 show that the vehicle slips and spins after 7.5 s without stability control. The maximum yaw velocity reached −62.779°/s, and the maximum sideslip angle reached 70.740°, which are much more than the reference model. Thus, the vehicle stability control system was activated. Furthermore, the model identified that the current intention was normal turning; hence, the reference yaw rate was no longer modified. The vehicle stability system reduced the maximum yaw rate by 65.1%, and the maximum sideslip angle by 92.7%. The stability system therefore improved the vehicle stability effectively. The vehicle stability control system enables the vehicle to run much more smoothly in comparison, preventing dangerous conditions such as sideslip and spin. The actual driving trajectory is shown in Figure 20. Meanwhile, it can be seen from Table 4 that both the controller proposed in this paper and SMC with the vehicle reference model considering the driver’s intention can track the control target well, and the control effect of some areas is better than that of the contrast controller.


## 6. Conclusion

This paper proposes a vehicle stability control strategy considering the driver’s turning intention for dual-motor drive electric vehicle.(1)
The upper controller of the hierarchical vehicle stability control system was constructed with the modified reference model as the control target. This model is modified by using turning urgency coefficient which is calculated on the basis of the recognition results of the GHMM/GGAP-RBF hybrid turning intention model. The lower layer controller takes the minimized sum of the square of operating working load rate of each tire as optimization objective and takes motor and road adhesion conditions as constraints to optimize the allocation additional yaw torque. The results show that the proposed vehicle stability control strategy can work satisfactorily and effectively improve the vehicle stability.(2)
Further studies could, apart from the turning intention, also consider the driver’s acceleration intention, braking intention, and complex intention, to improve the stability control strategy,X-by-wire system, and advanced assisted driving system. This could be expected to improve the safety and comfort of vehicle operation.



---

*Source: 3143620-2020-01-13.xml*