# Model Predictive Controller-Based Optimal Slip Ratio Control System for Distributed Driver Electric Vehicle

Authors: Qingxian Li; Liangjiang Liu; Xiaofang Yuan
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The slip ratio control is an important research topic in in-wheel-motored electric vehicles (EVs). Traditional control methods are usually designed for some specified modes. Therefore, the optimal slip ratio control cannot be achieved while vehicles work under various modes. In order to achieve the optimal slip ratio control, a novel model predictive controller-based optimal slip ratio control system (MPC-OSRCS) is proposed. The MPC-OSRCS includes three parts, a road surface adhesion coefficient identifier, an operation mode recognizer, and an MPC based-optimal slip ratio control. The current working road surface is identified by the road surface adhesion coefficient identifier, and a modified recursive Bayes theorem is used to compute the matching degree between current road surfaces and reference road surfaces. The current operation state is recognized by the operation mode recognizer, and a fuzzy logic method is applied to compute the matching degree between actual operation state and reference operation modes. Then, a parallel chaos optimization algorithm (PCOA)-based MPC is used to achieve the optimal control under various operation modes and different road surfaces. The MPC-OSRCS for EV is verified on simulation platform and simulation results under various conditions to show the significant performance.


---

## Body

## 1. Introduction

With the development of social, the environment is getting worse [1]. Compared with traditional vehicle, EVs have great advantage in decreasing environment pollution. Therefore, more and more scholars are involved in the research of EVs [2, 3].In the recent years, the slip ratio control is researched by scholars [4]. A model predictive control-based slip control was proposed [5];the wheel slip ratio was controlled on a stable zone rather than an optimal value. Aiming at improving EV safety, a sliding mode framework control system was extended [6]. In order to achieve the multiobjective optimization control, an MPC-based slip control system for EVs was proposed [7]. A wheel slip control algorithm combined on wheel slip ratio and wheel acceleration regulation was proposed [8]. In order to achieve traction control, a decoupling state feedback controller based on the uncertain frictional coefficient was derived [9]. A robust and fast wheel slip control based on the moving sliding surface technique was proposed [10]. A model predictive controller-based multimodel system by optimal slip ratio control was proposed [11]. The torque is assumed as a constant and the effect of the road surface adhesion coefficient is not considered in this paper. However, the EVs usually work under various operation modes and different road surfaces, while the present research studies of EVs focus on some typical modes and cannot achieve good performance under various operation modes and different road surfaces.In order to solve this problem, a novel MPC-OSRCS is proposed in this paper. The MPC-OSRCS includes three parts, a road surface adhesion coefficient identifier, an operation mode recognizer, and an MPC based-optimal slip ratio control. The current working road surface is identified by the road surface adhesion coefficient identifier. In order to accurately describe the state of road surface, five road surface reference models are established. The matching degree between the actual road surfaces and five reference road surfaces is computed. A modified recursive Bayes theorem is used to calculate the matching degree. The operation mode recognizer is used to recognize the current operation state. In order to accurately describe the operation state, three operation mode reference models are established. The matching degree between the actual operation state and three operation mode reference models is computed by the fuzzy method. The PCOA [12] is used to realize the optimal design of MPC. Finally, the control output of MPC-OSRCS is computed by the weighted output of each model to achieve optimal slip ratio control under various operation states and different road surfaces.The main contributions of this paper cover the following points. (1) The state of EV is divided into fifteen kinds of typical modes, and the reference model is established for these fifteen modes, which separately represent the link of five road surfaces and three operation modes. (2) The identifier is designed for recognizing road surface adhesion coefficient and operation mode, respectively. The output of the identifier presents the matching degree between the state of actual EV and each typical model. (3) Aiming at obtaining the control output under each state of EV, each type of matching coefficient is substituted into the controller and the weighted output of each state of EV makes up the output of MPC-OSRCS. (4) The optimal design of MPC is achieved by PCOA because of the global optimization with fast and accurate performance can be achieved by it.


## 2. System Model and Problem Statement

In order to achieve the optimal slip ratio control, a model of EV is established in this section, and the model of EV mainly consists of vehicle, tyre, and motor torque. Tyre model and vehicle are established according to the EV dynamics. The symbols in this paper and their physical meanings are displayed in Table1.Table 1
Symbols of the vehicle model.
Definition
Symbol
Unit
Vehicle mass
m
kg
Vehicle yaw moment of inertia
I
z
kgm
2
Axle tread
d
m
Distance from c.g. to front and rear axle
l
f
,
l
r
m
Front and rear track width
D
f
,
D
r
m
Height of c.g.
h
m
Front and rear cornering stiffness
C
f, Cr
N/rad
Resistance and inductance of motor
R, L
Ω
Flux of motor
ρ
f
Wb
Tyre radius
r
m
Sideslip angle in vehicle body
β
rad
Sideslip angle of tyre
α
rad
Yaw rate in vehicle body
γ
rad/s
Angular velocity of tyre
ω
rad/s
Pitch rate in vehicle body
ρ
rad/s
Front wheel steering angle from controller and driver
δ
f, δ
rad
Tyre longitudinal slip ratio
κ
—
Later force of front or rear tyres
F
y
f
,
F
y
r
N
Moment of inertia of each wheel
J
kgm
2
Longitudinal or lateral acceleration
a
x
,
a
y
m
/
s
2
Friction coefficient of road
μ
—
Longitudinal torque of each wheel
T
e
Nm
Total drive torque from driver pedal command
T
t
Nm
Additional yaw moment of each wheel
M
^
z
Nm
Fitting coefficients
K
a
,
K
b
—


### 2.1. Vehicle Model

The rotational movement of EV and longitudinal motion was included in the vehicle model in this paper, and the vehicle model was established based on two degree-of-freedom (2DOF) plane and shown in Figure1. Many terms in this section are explained in detail in Rajamanis book [13]. The vehicle model can be described as follows.Figure 1
2DOF vehicle model.Longitudinal model:(1)m⋅v˙=Fxfl+Fxfr+Fxrl+Fxrr,where Fxfl, Fxfr, Fxrl and Fxrr are longitudinal forces of four wheels, respectively.Rotation movement:(2)J⋅ω˙=Te−Tb−r·Fx,where Tb is motor braking torque.


### 2.2. Tyre Model with Various Actuators

With the Dugoff tyre model [14], the longitudinal force of tyre is shown as follows:(3)Fx=Cx⋅κ1+κ⋅fS,fS=2−SS,S≤1,1,S>1,S=μ⋅Fz⋅1+κ2Cx⋅κ2+Cy⋅tanα2,where Fx is the longitudinal force, Cx and Cy are the tire stiffness of longitudinal and lateral, respectively, μ is the road surface adhesion coefficient, κ is the longitudinal slip ratio, α is the tire slideslips angle, and Fz is the vertical load. The sideslip of tyre is not considered; then, the sideslip angle of tyre is considered as zero. Therefore, equation (3) can be rewritten as follows:(4)Fx=μ⋅Fz−μ2⋅Fz2⋅1+κ4⋅Cx⋅κ,μ⋅Fz⋅1+κ2⋅Cx⋅κ≤1,Cx⋅κ1+κ,μ⋅Fz⋅1+κ2⋅Cx⋅κ>1.The value of the wheel normal force could be calculated as follows:(5)Fzfl=mLrLf+Lrg−hcgLf+Lrax12−hcgd⋅gay,Fzfr=mLrLf+Lrg−hcgLf+Lrax12+hcgd⋅gay,Fzrl=mLrLf+Lrg+hcgLf+Lrax12−hcgd⋅gay,Fzrr=mLrLf+Lrg+hcgLf+Lrax12+hcgd⋅gay.Theκ represents the wheel slip ratio, which relates to the wheel speed and the vehicle speed. The relationship can be calculated as follows:(6)κi=vx−ωi⋅rmaxvx,ωi⋅r,ψ,i=fl,fr,rl,rr.In order to avoidvx⟶0⇒κi⟶∞, the ψ is set to 0.1 m/s. The longitudinal motion is only considered in this paper. In order to realize torque control, a torque balance equation is established based on vehicle characteristic for each wheel. In this case, the weight of EV is assigned to four tyres and every part equals to one quarter of the weight of the vehicle. Therefore, the single tyre model can be described as follows:(7)14mv˙=Fx.Driving motion: the EV works under uniform speed or acceleration operation mode whiler⋅ω>v. According to equations (1)–(5), (11), and (12), the κ can be described as follows, where vx = v:(8)κ˙=−r2Jv−κ+11/4mv⋅Fx+rJvTe−Tb.Braking motion: the EV works under braking operation mode whilev>r·ω. According to equations (1)–(5), (11), and (12), the κ can be described as follows:(9)κ˙=r2Jv+κ+11/4mv⋅Fx−rJvTe−Tb.According to equations (6) and (12), dμ/dκ is shown as follows [15]:(10)dμdκ=1/4⋅m⋅d2v/dt2+Fz⋅μ2/4⋅Cx⋅κ21−2⋅Fz⋅μ/4⋅Cx+Fz⋅μ/4⋅Cx⋅κ.


### 2.3. Control Problem Formulation

In order to realize the optimal slip ratio control, some problems should be considered in this paper:(1)
Various operation modes and different road surfaces problem: EVs usually work under different operation modes and various road surfaces. A mount of research studies show that the slip ratio is related to the adhesion coefficient and operation mode [16–18]. Therefore, different operation modes and various adhesion coefficients should be considered.(2)
Constraints problem: In order to keep the EV in a stability state, all the control parameters should be kept at constraint range. Theκ should be limited within the stable slip ratio range for each road surface. κp is the optimal for each road surface:(11)−κmax≤κ≤κmax,κmax=κp.In addition, the motor torque commandT should be limit within the maximum motor output torque Tmax [7], and it is shown as follows:(12)−Tmax≤T≤Tmax.


## 2.1. Vehicle Model

The rotational movement of EV and longitudinal motion was included in the vehicle model in this paper, and the vehicle model was established based on two degree-of-freedom (2DOF) plane and shown in Figure1. Many terms in this section are explained in detail in Rajamanis book [13]. The vehicle model can be described as follows.Figure 1
2DOF vehicle model.Longitudinal model:(1)m⋅v˙=Fxfl+Fxfr+Fxrl+Fxrr,where Fxfl, Fxfr, Fxrl and Fxrr are longitudinal forces of four wheels, respectively.Rotation movement:(2)J⋅ω˙=Te−Tb−r·Fx,where Tb is motor braking torque.


## 2.2. Tyre Model with Various Actuators

With the Dugoff tyre model [14], the longitudinal force of tyre is shown as follows:(3)Fx=Cx⋅κ1+κ⋅fS,fS=2−SS,S≤1,1,S>1,S=μ⋅Fz⋅1+κ2Cx⋅κ2+Cy⋅tanα2,where Fx is the longitudinal force, Cx and Cy are the tire stiffness of longitudinal and lateral, respectively, μ is the road surface adhesion coefficient, κ is the longitudinal slip ratio, α is the tire slideslips angle, and Fz is the vertical load. The sideslip of tyre is not considered; then, the sideslip angle of tyre is considered as zero. Therefore, equation (3) can be rewritten as follows:(4)Fx=μ⋅Fz−μ2⋅Fz2⋅1+κ4⋅Cx⋅κ,μ⋅Fz⋅1+κ2⋅Cx⋅κ≤1,Cx⋅κ1+κ,μ⋅Fz⋅1+κ2⋅Cx⋅κ>1.The value of the wheel normal force could be calculated as follows:(5)Fzfl=mLrLf+Lrg−hcgLf+Lrax12−hcgd⋅gay,Fzfr=mLrLf+Lrg−hcgLf+Lrax12+hcgd⋅gay,Fzrl=mLrLf+Lrg+hcgLf+Lrax12−hcgd⋅gay,Fzrr=mLrLf+Lrg+hcgLf+Lrax12+hcgd⋅gay.Theκ represents the wheel slip ratio, which relates to the wheel speed and the vehicle speed. The relationship can be calculated as follows:(6)κi=vx−ωi⋅rmaxvx,ωi⋅r,ψ,i=fl,fr,rl,rr.In order to avoidvx⟶0⇒κi⟶∞, the ψ is set to 0.1 m/s. The longitudinal motion is only considered in this paper. In order to realize torque control, a torque balance equation is established based on vehicle characteristic for each wheel. In this case, the weight of EV is assigned to four tyres and every part equals to one quarter of the weight of the vehicle. Therefore, the single tyre model can be described as follows:(7)14mv˙=Fx.Driving motion: the EV works under uniform speed or acceleration operation mode whiler⋅ω>v. According to equations (1)–(5), (11), and (12), the κ can be described as follows, where vx = v:(8)κ˙=−r2Jv−κ+11/4mv⋅Fx+rJvTe−Tb.Braking motion: the EV works under braking operation mode whilev>r·ω. According to equations (1)–(5), (11), and (12), the κ can be described as follows:(9)κ˙=r2Jv+κ+11/4mv⋅Fx−rJvTe−Tb.According to equations (6) and (12), dμ/dκ is shown as follows [15]:(10)dμdκ=1/4⋅m⋅d2v/dt2+Fz⋅μ2/4⋅Cx⋅κ21−2⋅Fz⋅μ/4⋅Cx+Fz⋅μ/4⋅Cx⋅κ.


## 2.3. Control Problem Formulation

In order to realize the optimal slip ratio control, some problems should be considered in this paper:(1)
Various operation modes and different road surfaces problem: EVs usually work under different operation modes and various road surfaces. A mount of research studies show that the slip ratio is related to the adhesion coefficient and operation mode [16–18]. Therefore, different operation modes and various adhesion coefficients should be considered.(2)
Constraints problem: In order to keep the EV in a stability state, all the control parameters should be kept at constraint range. Theκ should be limited within the stable slip ratio range for each road surface. κp is the optimal for each road surface:(11)−κmax≤κ≤κmax,κmax=κp.In addition, the motor torque commandT should be limit within the maximum motor output torque Tmax [7], and it is shown as follows:(12)−Tmax≤T≤Tmax.


## 3. Design of Control System

In order to realize the optimal slip ratio control when EV works on various operation modes and different road surfaces, a novel MPC-OSRCS is proposed in this paper. Figure2 shows the control structure of MPC-OSRCS. In the MPC-OSRCS, two recognizers are established according to the characteristics of EV, and the optimal slip ratio control is realized by a MPC under various operation modes and different road surfaces. The MPC-OSRCSs can be explained as follows. (1) The road surface adhesion coefficient identifier is used to identify the current working road surface. It includes five reference models and represents five typical road surfaces, respectively. Each reference model can accurately represent the state of road surfaces. A modified recursive Bayes theorem is used to calculate the matching degree between the actual road surfaces and five reference road surfaces. (2) The operation mode recognizer is used to recognize the current operation state. In order to accurately describe the state of operation state, three operation mode reference models are established. A fuzzy logic method is used to calculate the matching degree between the actual operation state and three reference operation modes. (3) The optimal slip ratio control under various operation modes and different road surfaces is realized by a MPC, and the optimal design of MPC is achieved by PCOA. (4) In order to obtain the control output under each state of EV, each type of matching coefficient is substituted into the controller and the weighted output of each state of EV makes up the output of MPC-OSRCS.Figure 2
The structure of MPC-OSRCS for EV.


### 3.1. The identifier of Road Surfaces Adhesion Coefficient

The slip ratioκ is related to the road surface adhesion coefficient μ, different road surfaces have various characteristic between adhesion coefficient μ, and slip ratio κ. μ−κ function curve equation is proposed by Zhang et al. [19]; μ−κ function is easy and accurate to describe the mathematical relationship between the slip ratio κ of tyre based on various roads and the adhesion coefficient μ as follows:(13)μκ=c1⋅1−exp−c2⋅κ−c3⋅κ,(14)dμdκ=c1⋅c2⋅exp−c2⋅κ−c3,where c1, c2, and c3 are the fitting coefficients for various road surfaces according to via experiment research. The parameters of five typical road surface are shown in Table 2.Table 2
Parameters of road surface adhesion coefficient.
Road surface
c
1
c
2
c
3
Dry asphalt
1.28
23.99
0.52
Wet asphalt
0.86
33.82
0.35
Wet pebbles
0.40
33.71
0.12
Snow
0.19
94.13
0.06
Ice
0.05
306.4
0.001In the design of road surface adhesion coefficient identifier, it is difficult to compute the matching degree between the actual road surfaces and five reference road surfaces. The switching idea between two operation modes has been introduced in many literature studies [20–22]. Based on the same method, this theory can also be used to identify road conditions.A modified recursive Bayes theorem is used to compute the weight coefficientξj, where ξj represents the matching degree between the actual road surfaces and each reference road surface. The former relative error is applied to the calculation of ξj and then the smooth switching can be achieved. The posterior probability is evaluated by a modified recursive Bayes theorem for the jth model (or value) at the kth time instant and can be described as follows:(15)Pj,k=H⋅exp−1/2εj,kT⋅G⋅εj,k⋅Pj,k−1∑j=15H⋅exp−1/2εj,kT⋅G⋅εj,k⋅Pj,k−1,j=1,2,3,4,5,where εj,k=ymk−yjk/ymk represents the relative error between the actual output state yjk and the reference value ymk. Pjk represents the posterior probability in the jth linearized model at k moment. H=h10000h20000h30000h4 represents the weight matrix between state variable kj, hee=1,2,3,4 represents the influence factor and relates to control performance, and G represents a time invariant weight matrix and usually is selected to a diagonal. According to the normal distribution, G represents the inverse matrix of the residual covariance. The higher parameter of G means that the residual variance is small and the confidence in the residual of each model is greater. The higher the value of the element of G is, the ability of rejection of model with large residual is more stronger. The weight coefficient of road surfaces ξj is described as follows:(16)ξj=Pj,k∑j=13Pj,k,j=1,2,3,4,5.Here,ξj is a value between 0 and 1. The summation of ξj is 1:(17)∑j=15ξj=1.


### 3.2. The Recognizer of Operation Modes

The T-S fuzzy controller is applied to recognize various operation modes in this paper, and it is proposed by the scholars Takagi and Sugeno in 1985. The actual of pedal opening (PO) and the change rate of pedal opening (RO) are the inputs of the T-S fuzzy controller. If the EV is under driving modes, the pedal opening is defined as positive. Otherwise, the PO is defined as negative. The PO and RO are different under various operation modes. The output of the T-S fuzzy controller is membership valueηi,i=1,2,3, which represents three typical operation modes: acceleration, uniform speed, and deceleration.The values of PO and RO are limited within−100,100. The value of PO is defined as seven kinds of states, there are zero (Z), very small (VS), small (S), medium (M), big (B), very big (VB), and great (G), respectively. At the same time, the value of RO is defined as five kinds of states, and there are very small (VS), small (S), medium (M), big (B), and very big (VB), respectively. The value of η is defined as four kinds of states, and there are very small (VS), small (S), medium (M), and big (B). The membership functions for PO and RO are set, as shown in Figure 3.Figure 3
The membership functions. (a) PO. (b) RO. (c)ηi.
(a)
(b)
(c)Table3 shows the rule set in the fuzzy inference system (FIS). The rule set includes 7∗5=35 rules.Table 3
Fuzzy logic rule set.
PO
RO
η
1
η
2
η
3
Z
VS
B
S
S
Z
S
B
S
S
Z
M
B
S
S
Z
B
B
S
S
Z
VB
B
S
S
VS
VS
B
S
S
VS
S
B
S
S
VS
M
B
S
S
VS
B
B
S
S
VS
VB
B
S
S
S
VS
B
S
S
S
S
B
M
S
S
M
B
M
S
S
B
B
M
S
S
VB
B
S
S
M
VS
M
B
S
M
S
M
B
S
M
M
M
B
M
M
B
S
M
B
M
VB
S
M
B
B
VS
B
M
S
B
S
M
B
S
B
M
M
B
M
B
B
S
M
B
B
VB
S
S
B
VB
VS
B
S
S
VB
S
S
B
S
VB
M
M
B
M
VB
B
S
M
B
VB
VB
S
S
B
G
VS
B
S
S
G
S
M
B
S
G
M
S
B
S
G
B
S
M
B
G
VB
S
S
BRule1: if PO is “Z,” RO is “VS,” thenη1 is “B,” η2 is “S,” and η3 is “S.”…Rule35: if PO is “G,” RO is “VB,” thenη1 is “S,” η2 is “S,” and η3 is “M.”In this way, the whole rules can be obtained.Here,∑i=13λi=1 and λi∈0,1, where i=1,2,3 Normalization processing is carried out to calculate the weight coefficient λi as follows:(18)λi=ηi∑j=13ηi,where ηi is membership value of three typical operation modes.


### 3.3. Design of the MPC

The optimal and constraint problem are efficiently solved by MPC. MPC is selected to optimumκ for each operation mode and each road surface in this paper. MPC has a great advantage to achieve satisfy diverse and even conflicting requirements of vehicles [23] as its characteristic. In addition, the effect of model mismatch and unmeasured disturbance can be attenuated by MPC [24, 25]. Moreover, the effect of switching between various road surfaces and different operation modes can be attenuated by the weighted output of each MPC.


#### 3.3.1. Model Predictive Control Law

x represents state vector, u represents control vector, and y represents output vector, and they can be expressed as follows:(19)x=κflκfrκrlκrr,u=TflTfrTrlTrr,y=κflκfrκrlκrr.In order to achieve the optimal MPC design, time state-space models in equations (12)-(13) are discrete by the Euler method. k represents the sample time and is defined as k=intt/Ts, where running time is defined as t and fixed step size is defined as Ts. At sample time k, the predictive state is calculated as follows:(20)xk+1=fkx+1,uk⋅Ts+xk,yk=C⋅xk.Then, the discretization of the state space model of equation (12) is described as equation(27), and the discretization of state space of equation (13) is described as equation (28):(21)xek+1=−r2Jv−xek+11/4mv⋅Fx⋅Ts⋅xek+rJv·uek⋅Ts⋅xek+xek,e=1,2,3,4,xek+1=r2Jv+xek+11/4mv⋅Fx−rJv⋅uek⋅Ts+1⋅xek.According to the principle of MPC,p represents a predictive horizon and m represents a control horizon. m=p=6 in this paper. Uk represents the optimization vectors and Yk represents the predictive control output. They can be expressed as follows:(22)Uk=ukkuk+1k⋅uk+m−1km×1,Yk=ykkyk+1k⋅yk+mkp×1,where Uk is defined as an array of control input u and Yk represents the output vector from sampling time k to sampling time k+i. Moreover, Rk represents reference sequence. uk is defined as the control input change and can be calculated by Δuk=uk−uk−1. They can be defined as follows:(23)Rk=rkrk⋅rkp×1,ΔUk=△ukk△uk+pk⋅△uk+m−1km×1.According to the principle of MPC, the relationship ofuk and Uk can be described as follows:(24)uk=1,0,…,0⋅Uk,where uk represents optimal vectors, which is used to compute the optimal control action of the EV system.


#### 3.3.2. Design of the Objective Function

In order to realize the optimal control, the multiple objective function of the MPC is established and it is achieved by minimization of this multiple objective function. The optimization function includes two parts. One is the motor torqueT; it includes drive torque Te and brake torques Tb; the other is slip ratio κ of the four wheels. The cost function is described as follows:(1)
The main importance is to make the actual slip ratioκ track the optimal slip ratio κp, whatever the conditions of EV are. Therefore, the first cost function can be described as follows:(25)J1=Yk−RkQ2,whereQ represents a positive weight factors and can be used to adjust the performance of tracking.(2)
In order to save energy, the sum of torque command squares should be controlled as small as possible. Therefore, the second cost function can be calculated as follows:(26)J2=UkR2,whereR represents a positive weight factors and can be used to adjust the performance of T.(3)
Aiming at achieving the optimal slip ratio control, the actual slip ratioκ needs to track the optimal slip ratio κp as accurate as possible. Hence, the third cost function can be calculated as follows:(27)J3=κk−κpkF2.(4)
The tyre longitudinal performance is related toT. Therefore, in order to ensure the longitudinal stability of EV, T is limited within Tmax, and the forth cost function can be calculated as follows:(28)J4=Uk−TmaxG2.Therefore, the total cost function can be defined as follows:(29)J=J1+J2+J3+J4.


#### 3.3.3. The Process of Optimal Design

The optimal design is realized by a novel PCOA, and we have carried out some works about PCOA. In addition, the PCOA has great advantage to realize global optimum. Objective function in the PCOA is described in equation (30) and it can be rewritten as follows:(30)minfX=fx1,x2,x3,x4,xi∈Li,Pi,where X=fx1,x2,x3,x4∈W4 is a vector in the 4-dimensional decision variables space, xi represents different variables parameters and xi∈Li,Pi, Li represents the lower bound of the ith variable, and Pi represents the upper bound of the ith variable. Then, PCOA evolves a stochastic population of N candidate individuals with 4-dimensional parameter vectors; the candidate individuals will experience twice carrier wave mechanism. The different chaotic traces are obtained by the first carrier wave; the search precision is enhanced by the second carrier wave.


### 3.4. The Control Output of MPC-OSRCS

As described above, the weight coefficient of different reference models and the output of each MPC are calculated. At each sampling point, the control output of MPC-OSRCS is constituted by the output of each MPCui·jk and its weight coefficient. Therefore, the control output of MPC-OSRCS can be described as follows:(31)uk=∑i=13∑j=15λi⋅ξj⋅ui·jk,where uk represents the control output of MPC-OSRCS, the ui·jk represents the i⋅jth reference model, and λi and ξj are the weight coefficients.


## 3.1. The identifier of Road Surfaces Adhesion Coefficient

The slip ratioκ is related to the road surface adhesion coefficient μ, different road surfaces have various characteristic between adhesion coefficient μ, and slip ratio κ. μ−κ function curve equation is proposed by Zhang et al. [19]; μ−κ function is easy and accurate to describe the mathematical relationship between the slip ratio κ of tyre based on various roads and the adhesion coefficient μ as follows:(13)μκ=c1⋅1−exp−c2⋅κ−c3⋅κ,(14)dμdκ=c1⋅c2⋅exp−c2⋅κ−c3,where c1, c2, and c3 are the fitting coefficients for various road surfaces according to via experiment research. The parameters of five typical road surface are shown in Table 2.Table 2
Parameters of road surface adhesion coefficient.
Road surface
c
1
c
2
c
3
Dry asphalt
1.28
23.99
0.52
Wet asphalt
0.86
33.82
0.35
Wet pebbles
0.40
33.71
0.12
Snow
0.19
94.13
0.06
Ice
0.05
306.4
0.001In the design of road surface adhesion coefficient identifier, it is difficult to compute the matching degree between the actual road surfaces and five reference road surfaces. The switching idea between two operation modes has been introduced in many literature studies [20–22]. Based on the same method, this theory can also be used to identify road conditions.A modified recursive Bayes theorem is used to compute the weight coefficientξj, where ξj represents the matching degree between the actual road surfaces and each reference road surface. The former relative error is applied to the calculation of ξj and then the smooth switching can be achieved. The posterior probability is evaluated by a modified recursive Bayes theorem for the jth model (or value) at the kth time instant and can be described as follows:(15)Pj,k=H⋅exp−1/2εj,kT⋅G⋅εj,k⋅Pj,k−1∑j=15H⋅exp−1/2εj,kT⋅G⋅εj,k⋅Pj,k−1,j=1,2,3,4,5,where εj,k=ymk−yjk/ymk represents the relative error between the actual output state yjk and the reference value ymk. Pjk represents the posterior probability in the jth linearized model at k moment. H=h10000h20000h30000h4 represents the weight matrix between state variable kj, hee=1,2,3,4 represents the influence factor and relates to control performance, and G represents a time invariant weight matrix and usually is selected to a diagonal. According to the normal distribution, G represents the inverse matrix of the residual covariance. The higher parameter of G means that the residual variance is small and the confidence in the residual of each model is greater. The higher the value of the element of G is, the ability of rejection of model with large residual is more stronger. The weight coefficient of road surfaces ξj is described as follows:(16)ξj=Pj,k∑j=13Pj,k,j=1,2,3,4,5.Here,ξj is a value between 0 and 1. The summation of ξj is 1:(17)∑j=15ξj=1.


## 3.2. The Recognizer of Operation Modes

The T-S fuzzy controller is applied to recognize various operation modes in this paper, and it is proposed by the scholars Takagi and Sugeno in 1985. The actual of pedal opening (PO) and the change rate of pedal opening (RO) are the inputs of the T-S fuzzy controller. If the EV is under driving modes, the pedal opening is defined as positive. Otherwise, the PO is defined as negative. The PO and RO are different under various operation modes. The output of the T-S fuzzy controller is membership valueηi,i=1,2,3, which represents three typical operation modes: acceleration, uniform speed, and deceleration.The values of PO and RO are limited within−100,100. The value of PO is defined as seven kinds of states, there are zero (Z), very small (VS), small (S), medium (M), big (B), very big (VB), and great (G), respectively. At the same time, the value of RO is defined as five kinds of states, and there are very small (VS), small (S), medium (M), big (B), and very big (VB), respectively. The value of η is defined as four kinds of states, and there are very small (VS), small (S), medium (M), and big (B). The membership functions for PO and RO are set, as shown in Figure 3.Figure 3
The membership functions. (a) PO. (b) RO. (c)ηi.
(a)
(b)
(c)Table3 shows the rule set in the fuzzy inference system (FIS). The rule set includes 7∗5=35 rules.Table 3
Fuzzy logic rule set.
PO
RO
η
1
η
2
η
3
Z
VS
B
S
S
Z
S
B
S
S
Z
M
B
S
S
Z
B
B
S
S
Z
VB
B
S
S
VS
VS
B
S
S
VS
S
B
S
S
VS
M
B
S
S
VS
B
B
S
S
VS
VB
B
S
S
S
VS
B
S
S
S
S
B
M
S
S
M
B
M
S
S
B
B
M
S
S
VB
B
S
S
M
VS
M
B
S
M
S
M
B
S
M
M
M
B
M
M
B
S
M
B
M
VB
S
M
B
B
VS
B
M
S
B
S
M
B
S
B
M
M
B
M
B
B
S
M
B
B
VB
S
S
B
VB
VS
B
S
S
VB
S
S
B
S
VB
M
M
B
M
VB
B
S
M
B
VB
VB
S
S
B
G
VS
B
S
S
G
S
M
B
S
G
M
S
B
S
G
B
S
M
B
G
VB
S
S
BRule1: if PO is “Z,” RO is “VS,” thenη1 is “B,” η2 is “S,” and η3 is “S.”…Rule35: if PO is “G,” RO is “VB,” thenη1 is “S,” η2 is “S,” and η3 is “M.”In this way, the whole rules can be obtained.Here,∑i=13λi=1 and λi∈0,1, where i=1,2,3 Normalization processing is carried out to calculate the weight coefficient λi as follows:(18)λi=ηi∑j=13ηi,where ηi is membership value of three typical operation modes.


## 3.3. Design of the MPC

The optimal and constraint problem are efficiently solved by MPC. MPC is selected to optimumκ for each operation mode and each road surface in this paper. MPC has a great advantage to achieve satisfy diverse and even conflicting requirements of vehicles [23] as its characteristic. In addition, the effect of model mismatch and unmeasured disturbance can be attenuated by MPC [24, 25]. Moreover, the effect of switching between various road surfaces and different operation modes can be attenuated by the weighted output of each MPC.


### 3.3.1. Model Predictive Control Law

x represents state vector, u represents control vector, and y represents output vector, and they can be expressed as follows:(19)x=κflκfrκrlκrr,u=TflTfrTrlTrr,y=κflκfrκrlκrr.In order to achieve the optimal MPC design, time state-space models in equations (12)-(13) are discrete by the Euler method. k represents the sample time and is defined as k=intt/Ts, where running time is defined as t and fixed step size is defined as Ts. At sample time k, the predictive state is calculated as follows:(20)xk+1=fkx+1,uk⋅Ts+xk,yk=C⋅xk.Then, the discretization of the state space model of equation (12) is described as equation(27), and the discretization of state space of equation (13) is described as equation (28):(21)xek+1=−r2Jv−xek+11/4mv⋅Fx⋅Ts⋅xek+rJv·uek⋅Ts⋅xek+xek,e=1,2,3,4,xek+1=r2Jv+xek+11/4mv⋅Fx−rJv⋅uek⋅Ts+1⋅xek.According to the principle of MPC,p represents a predictive horizon and m represents a control horizon. m=p=6 in this paper. Uk represents the optimization vectors and Yk represents the predictive control output. They can be expressed as follows:(22)Uk=ukkuk+1k⋅uk+m−1km×1,Yk=ykkyk+1k⋅yk+mkp×1,where Uk is defined as an array of control input u and Yk represents the output vector from sampling time k to sampling time k+i. Moreover, Rk represents reference sequence. uk is defined as the control input change and can be calculated by Δuk=uk−uk−1. They can be defined as follows:(23)Rk=rkrk⋅rkp×1,ΔUk=△ukk△uk+pk⋅△uk+m−1km×1.According to the principle of MPC, the relationship ofuk and Uk can be described as follows:(24)uk=1,0,…,0⋅Uk,where uk represents optimal vectors, which is used to compute the optimal control action of the EV system.


### 3.3.2. Design of the Objective Function

In order to realize the optimal control, the multiple objective function of the MPC is established and it is achieved by minimization of this multiple objective function. The optimization function includes two parts. One is the motor torqueT; it includes drive torque Te and brake torques Tb; the other is slip ratio κ of the four wheels. The cost function is described as follows:(1)
The main importance is to make the actual slip ratioκ track the optimal slip ratio κp, whatever the conditions of EV are. Therefore, the first cost function can be described as follows:(25)J1=Yk−RkQ2,whereQ represents a positive weight factors and can be used to adjust the performance of tracking.(2)
In order to save energy, the sum of torque command squares should be controlled as small as possible. Therefore, the second cost function can be calculated as follows:(26)J2=UkR2,whereR represents a positive weight factors and can be used to adjust the performance of T.(3)
Aiming at achieving the optimal slip ratio control, the actual slip ratioκ needs to track the optimal slip ratio κp as accurate as possible. Hence, the third cost function can be calculated as follows:(27)J3=κk−κpkF2.(4)
The tyre longitudinal performance is related toT. Therefore, in order to ensure the longitudinal stability of EV, T is limited within Tmax, and the forth cost function can be calculated as follows:(28)J4=Uk−TmaxG2.Therefore, the total cost function can be defined as follows:(29)J=J1+J2+J3+J4.


### 3.3.3. The Process of Optimal Design

The optimal design is realized by a novel PCOA, and we have carried out some works about PCOA. In addition, the PCOA has great advantage to realize global optimum. Objective function in the PCOA is described in equation (30) and it can be rewritten as follows:(30)minfX=fx1,x2,x3,x4,xi∈Li,Pi,where X=fx1,x2,x3,x4∈W4 is a vector in the 4-dimensional decision variables space, xi represents different variables parameters and xi∈Li,Pi, Li represents the lower bound of the ith variable, and Pi represents the upper bound of the ith variable. Then, PCOA evolves a stochastic population of N candidate individuals with 4-dimensional parameter vectors; the candidate individuals will experience twice carrier wave mechanism. The different chaotic traces are obtained by the first carrier wave; the search precision is enhanced by the second carrier wave.


## 3.3.1. Model Predictive Control Law

x represents state vector, u represents control vector, and y represents output vector, and they can be expressed as follows:(19)x=κflκfrκrlκrr,u=TflTfrTrlTrr,y=κflκfrκrlκrr.In order to achieve the optimal MPC design, time state-space models in equations (12)-(13) are discrete by the Euler method. k represents the sample time and is defined as k=intt/Ts, where running time is defined as t and fixed step size is defined as Ts. At sample time k, the predictive state is calculated as follows:(20)xk+1=fkx+1,uk⋅Ts+xk,yk=C⋅xk.Then, the discretization of the state space model of equation (12) is described as equation(27), and the discretization of state space of equation (13) is described as equation (28):(21)xek+1=−r2Jv−xek+11/4mv⋅Fx⋅Ts⋅xek+rJv·uek⋅Ts⋅xek+xek,e=1,2,3,4,xek+1=r2Jv+xek+11/4mv⋅Fx−rJv⋅uek⋅Ts+1⋅xek.According to the principle of MPC,p represents a predictive horizon and m represents a control horizon. m=p=6 in this paper. Uk represents the optimization vectors and Yk represents the predictive control output. They can be expressed as follows:(22)Uk=ukkuk+1k⋅uk+m−1km×1,Yk=ykkyk+1k⋅yk+mkp×1,where Uk is defined as an array of control input u and Yk represents the output vector from sampling time k to sampling time k+i. Moreover, Rk represents reference sequence. uk is defined as the control input change and can be calculated by Δuk=uk−uk−1. They can be defined as follows:(23)Rk=rkrk⋅rkp×1,ΔUk=△ukk△uk+pk⋅△uk+m−1km×1.According to the principle of MPC, the relationship ofuk and Uk can be described as follows:(24)uk=1,0,…,0⋅Uk,where uk represents optimal vectors, which is used to compute the optimal control action of the EV system.


## 3.3.2. Design of the Objective Function

In order to realize the optimal control, the multiple objective function of the MPC is established and it is achieved by minimization of this multiple objective function. The optimization function includes two parts. One is the motor torqueT; it includes drive torque Te and brake torques Tb; the other is slip ratio κ of the four wheels. The cost function is described as follows:(1)
The main importance is to make the actual slip ratioκ track the optimal slip ratio κp, whatever the conditions of EV are. Therefore, the first cost function can be described as follows:(25)J1=Yk−RkQ2,whereQ represents a positive weight factors and can be used to adjust the performance of tracking.(2)
In order to save energy, the sum of torque command squares should be controlled as small as possible. Therefore, the second cost function can be calculated as follows:(26)J2=UkR2,whereR represents a positive weight factors and can be used to adjust the performance of T.(3)
Aiming at achieving the optimal slip ratio control, the actual slip ratioκ needs to track the optimal slip ratio κp as accurate as possible. Hence, the third cost function can be calculated as follows:(27)J3=κk−κpkF2.(4)
The tyre longitudinal performance is related toT. Therefore, in order to ensure the longitudinal stability of EV, T is limited within Tmax, and the forth cost function can be calculated as follows:(28)J4=Uk−TmaxG2.Therefore, the total cost function can be defined as follows:(29)J=J1+J2+J3+J4.


## 3.3.3. The Process of Optimal Design

The optimal design is realized by a novel PCOA, and we have carried out some works about PCOA. In addition, the PCOA has great advantage to realize global optimum. Objective function in the PCOA is described in equation (30) and it can be rewritten as follows:(30)minfX=fx1,x2,x3,x4,xi∈Li,Pi,where X=fx1,x2,x3,x4∈W4 is a vector in the 4-dimensional decision variables space, xi represents different variables parameters and xi∈Li,Pi, Li represents the lower bound of the ith variable, and Pi represents the upper bound of the ith variable. Then, PCOA evolves a stochastic population of N candidate individuals with 4-dimensional parameter vectors; the candidate individuals will experience twice carrier wave mechanism. The different chaotic traces are obtained by the first carrier wave; the search precision is enhanced by the second carrier wave.


## 3.4. The Control Output of MPC-OSRCS

As described above, the weight coefficient of different reference models and the output of each MPC are calculated. At each sampling point, the control output of MPC-OSRCS is constituted by the output of each MPCui·jk and its weight coefficient. Therefore, the control output of MPC-OSRCS can be described as follows:(31)uk=∑i=13∑j=15λi⋅ξj⋅ui·jk,where uk represents the control output of MPC-OSRCS, the ui·jk represents the i⋅jth reference model, and λi and ξj are the weight coefficients.


## 4. Simulation and Analysis

The computer simulation is used to verify the control performance of the proposed MPC-OSRCS. The simulation is based on an 8DOF model simulation platform, as shown in Figure4.Figure 4
8DOF vehicle model.The simulation is carried under three operation modes and five road surfaces. Aiming at verifying the advantage of the MPC-OSRCS, the performance of the MPC-OSRCS is compared with the tractional MPC [18] and the model predictive controller-based multimodel control system (MPC-MMCS). The simulation time is set 5 s, and the fixed step size Ts is chosen as 5 ms. The simulation parameters are set as follows: Q=diag104,104,104,104, R=diag10,10,10,10, F=diag104,104,104,104, S=diag104,104,104,104, m=1359.8kg, Iz=1992.54kgm2, Lf=1.0628m, Lr=1.4852m, Df=1.0628m, Cy=23540N/rad, Cx=23101N/rad, J=0.3534kgm3, L=0.05Hm, ρf=0.1Wb, r=0.29m, and h=0.512m.


### 4.1. Test 1

The EV works on straight running in test 1, and the start speed of 0 km/h. The change of vehicle velocity and the change of road surfaces are considered in this simulation. The vehicle velocityv increases from 10 km/h to 80 km/h during 0−1.75s, keeps uniform speed during 1.75−3.25s, and decreases from 80 km/h to 10 km/h during 3.25−5s. The road surfaces are wet asphalt during 0−0.75s, the road surface is wet pebbles during 0.75−1.5s, the road surface is dry asphalt during 1.5−2.5s, the road surface is ice during 2.5−4s, the road surface is snow during 4−5s. The traditional method means that only a MPC under various operation modes and different road surfaces [17]. The input variables of the traditional MPC are four motor torques from EV, and the output variables are actual state values of the actual state value of κfl, κfr, κrl, and κrr.


#### 4.1.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure5.Figure 5
The simulation results of conventional MPC on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure5(a), the slip ratio is very big at the beginning, even close to 1. With the increase of velocity, the slip ratio is kept on a steady range. It is clearly shown in Figure 5(a) that the actual slip ratio has not changed following the road surfaces changes; in other words, the slip ratio is only a steady value rather than an optimal value. It is proved that the traditional MPC cannot be suitable to the multimodes. Therefore, the traditional MPC is difficult to achieve good performance and realize the optimal slip ratio control in actual working.We can see from Figure5(b) that the motor torques of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N. The motor torques of four wheels is not following the road surface changes, whatever the road surface is wet asphalt or snowy. The required motor torques of four wheels are different while the road surfaces changes. In the switching process between various operation modes and different road surfaces, the oscillatory behavior may occur.From the results of simulation, we can find that traditional MPC cannot realize the optimal slip ratio control while EV switches on various operation modes and different road surfaces; then, it is difficult to ensure the longitudinal stability of EVs.


#### 4.1.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure6. As Figure 6(a) shows, the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot achieve when the road surfaces changes. It is easy to find in Figure 6(b) that T is also kept within a stable range.Figure 6
The simulation results of MPC-MMCS on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)The wheel velocity is displayed in Figure6(c), the smoothly switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 6(d), the λi can accurately reflect the operation modes of EV.In conclusion, the MPC-MMCS can successfully recognize the operation modes of EV. Yet the optimal slip rate cannot achieve when the EV works under different road surfaces.


#### 4.1.3. The Performance of MPC-OSRCS

The simulation results of MPC-OSRCS are shown in Figure7. Note that the weight coefficient λi and ξj are calculated by the fuzzy method and a modified recursive Bayes theorem. They represent the matching degree between the state of actual EV and three reference operation modes and between the actual road surfaces and reference road surfaces, respectively.Figure 7
The simulation results of MPC-OSRCS on test 1. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torque.
(a)
(b)
(c)
(d)
(e)The biggest challenges in the proposed MPC-OSRCS are the optimal slip ratio control under different operation modes and various road surfaces. As shown in Figure7(a) that the optimal wheel ratio slip can be achieved no matter what the state of actual EV, in addition, the optimal slip rate follows the road surfaces change in electric vehicle operation. While the EV has switched different road surfaces, the optimal slip ratio also can smoothly achieved. We can see from the picture that the optimal slip ratio control can be realizes by proposed MPC-OSRCS when EV switches on various operation modes and different road surfaces.Figure7(b) illustrates the wheel velocity r⋅ω. We can see from Figure 6(b) that EV runs smoothly when it switches on various operation modes and different roads surfaces. We can see from Figures 7(c) and 7(d) that the λi represents the matching degree between the state of actual EV and each operation modes, and the ξj represents the matching degree between the actual roads surfaces and the reference road surfaces. As described above, the matching degree is almost identical to operation conditions.We can see from Figure7(e) that T changes smoothly while the state of EV changing. Whether the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small and the comfortable of EV is well.Compared with the traditional MPC, the proposed MPC-OSRCS can realize the smooth switching and optimal slip ratio control whatever the running state of EV. The simulation clearly proves that the proposed MPC-OSRCS can be suitable for different operation modes and various road surfaces and all can achieved smooth switching and optimal slip ratio control.


### 4.2. Test 2

The EV also is tested on various operation modes and different road surfaces in test 2. The starting speed of EV is 85 km/h and the straight running only is considered. The state of EV is deceleration, its velocity decreases from 80 km/h to 10 km/h during0−1.75s, and then the EV keeps the speed unchanged. At this time, the state of EV is uniform speed during 1.75−3.25s. At the next time, the EV is acceleration and the velocity of EV increase from 10 km/h to 80 km/h during 3.25−5s. At the same time, the road surfaces are also changed over time. At the beginning, the EV is driving on the ice road surface. After 0.75 s, the road surfaces are changed to snow road surface. After 0.75 s, the road surfaces are changed to wet pebbles road surface. After 1 s, the road surfaces are changed to wet asphalt road surface. After 1.5 s, the road surfaces are changed to dry asphalt road surface.


#### 4.2.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure8.Figure 8
The simulation results of conventional MPC on test2. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure8(a), we can see that the slip ratio is mainly affected by operation modes and the road surfaces adhesion coefficient changes can be ignored. However, the slip ratio is affected by operation modes and road surfaces, and the optimal slip ratio control cannot be achieved under different operation modes and various road surfaces. Therefore, it is proved that the conventional MPC cannot be suitable for the multimodes and the conventional MPC cannot achieve optimal slip ratio control performance in actual working.We can see from Figure8(b) that the motor torques remain unchanged. In order to be suitable for different operation modes and various road surfaces, the motor torque should change, following the state of EV. However, the motor torque of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N, whether the road surface is wet asphalt or snowy. Therefore, it is difficult to ensure the longitudinal stability of EV under traditional MPC, while EV works on various operation modes and different road surfaces. The conventional MPC cannot ensure driver comfort and longitudinal stability under this situation.From the simulation of traditional MPC, we can conclude that it is difficult to realize the optimal slip ratio control and smooth switching, while EV works on various operation modes and different road surfaces.


#### 4.2.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure9.Figure 9
The simulation results of MPC-MMCS on test 2. (a) Tire longitude slip ratioκ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)As shown in Figure9(a), the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot be achieved when the road surfaces changes. It is easy to find in Figure 9(b) that T is also kept within a stable range. The wheel velocity is displayed in Figure 9(c), and the smooth switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 9(d), the λi accurately reflects the operation modes of EV.


#### 4.2.3. The Performance of MPC-OSRCS

The simulation results of the proposed MPC-OSRCS are shown in Figure10.Figure 10
The simulation results of MPC-OSRCS on test 2. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torques.
(a)
(b)
(c)
(d)
(e)We can see from Figure10(a) the EV works under different road surfaces and various operation modes, and the slip ratio is controlled at the optimal slip ratio for the corresponding road surface, no matter how the state of EV changes. Compared with the traditional MPC, the slip ratio κ changes quickly when the state of EV changes, and it can be proved that the rapidity of the controller is good. The slip ratio can be effectively achieved by the MPC-OSRCS and can be suitable for various operation modes and different road surfaces.Figure10(b) illustrates the test on wheel velocity r⋅ω. It is observed that EV runs smoothly under various operation modes and different roads surfaces. We can see from Figures 10(c) and 10(d) that the identifier can quickly distinguish when the state of EV changes. As described above, the matching degree change following the state of EV changes and can precisely reflect the state of EV changes under different road surfaces and various operation modes.We can see from Figure10(e) that T changes smoothly while the state of EV change, whether it is the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small, and the comfortable of EV is well.Compared with the conventional MPC, the proposed MPC-OSRCS can identify different operation modes and various road surfaces, and optimal slip ratio control can be realized. The simulation results reveal that the proposed MPC-OSRCS can better ensure the longitudinal stability of EV under different operation modes and various road surfaces.


## 4.1. Test 1

The EV works on straight running in test 1, and the start speed of 0 km/h. The change of vehicle velocity and the change of road surfaces are considered in this simulation. The vehicle velocityv increases from 10 km/h to 80 km/h during 0−1.75s, keeps uniform speed during 1.75−3.25s, and decreases from 80 km/h to 10 km/h during 3.25−5s. The road surfaces are wet asphalt during 0−0.75s, the road surface is wet pebbles during 0.75−1.5s, the road surface is dry asphalt during 1.5−2.5s, the road surface is ice during 2.5−4s, the road surface is snow during 4−5s. The traditional method means that only a MPC under various operation modes and different road surfaces [17]. The input variables of the traditional MPC are four motor torques from EV, and the output variables are actual state values of the actual state value of κfl, κfr, κrl, and κrr.


### 4.1.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure5.Figure 5
The simulation results of conventional MPC on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure5(a), the slip ratio is very big at the beginning, even close to 1. With the increase of velocity, the slip ratio is kept on a steady range. It is clearly shown in Figure 5(a) that the actual slip ratio has not changed following the road surfaces changes; in other words, the slip ratio is only a steady value rather than an optimal value. It is proved that the traditional MPC cannot be suitable to the multimodes. Therefore, the traditional MPC is difficult to achieve good performance and realize the optimal slip ratio control in actual working.We can see from Figure5(b) that the motor torques of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N. The motor torques of four wheels is not following the road surface changes, whatever the road surface is wet asphalt or snowy. The required motor torques of four wheels are different while the road surfaces changes. In the switching process between various operation modes and different road surfaces, the oscillatory behavior may occur.From the results of simulation, we can find that traditional MPC cannot realize the optimal slip ratio control while EV switches on various operation modes and different road surfaces; then, it is difficult to ensure the longitudinal stability of EVs.


### 4.1.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure6. As Figure 6(a) shows, the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot achieve when the road surfaces changes. It is easy to find in Figure 6(b) that T is also kept within a stable range.Figure 6
The simulation results of MPC-MMCS on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)The wheel velocity is displayed in Figure6(c), the smoothly switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 6(d), the λi can accurately reflect the operation modes of EV.In conclusion, the MPC-MMCS can successfully recognize the operation modes of EV. Yet the optimal slip rate cannot achieve when the EV works under different road surfaces.


### 4.1.3. The Performance of MPC-OSRCS

The simulation results of MPC-OSRCS are shown in Figure7. Note that the weight coefficient λi and ξj are calculated by the fuzzy method and a modified recursive Bayes theorem. They represent the matching degree between the state of actual EV and three reference operation modes and between the actual road surfaces and reference road surfaces, respectively.Figure 7
The simulation results of MPC-OSRCS on test 1. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torque.
(a)
(b)
(c)
(d)
(e)The biggest challenges in the proposed MPC-OSRCS are the optimal slip ratio control under different operation modes and various road surfaces. As shown in Figure7(a) that the optimal wheel ratio slip can be achieved no matter what the state of actual EV, in addition, the optimal slip rate follows the road surfaces change in electric vehicle operation. While the EV has switched different road surfaces, the optimal slip ratio also can smoothly achieved. We can see from the picture that the optimal slip ratio control can be realizes by proposed MPC-OSRCS when EV switches on various operation modes and different road surfaces.Figure7(b) illustrates the wheel velocity r⋅ω. We can see from Figure 6(b) that EV runs smoothly when it switches on various operation modes and different roads surfaces. We can see from Figures 7(c) and 7(d) that the λi represents the matching degree between the state of actual EV and each operation modes, and the ξj represents the matching degree between the actual roads surfaces and the reference road surfaces. As described above, the matching degree is almost identical to operation conditions.We can see from Figure7(e) that T changes smoothly while the state of EV changing. Whether the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small and the comfortable of EV is well.Compared with the traditional MPC, the proposed MPC-OSRCS can realize the smooth switching and optimal slip ratio control whatever the running state of EV. The simulation clearly proves that the proposed MPC-OSRCS can be suitable for different operation modes and various road surfaces and all can achieved smooth switching and optimal slip ratio control.


## 4.1.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure5.Figure 5
The simulation results of conventional MPC on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure5(a), the slip ratio is very big at the beginning, even close to 1. With the increase of velocity, the slip ratio is kept on a steady range. It is clearly shown in Figure 5(a) that the actual slip ratio has not changed following the road surfaces changes; in other words, the slip ratio is only a steady value rather than an optimal value. It is proved that the traditional MPC cannot be suitable to the multimodes. Therefore, the traditional MPC is difficult to achieve good performance and realize the optimal slip ratio control in actual working.We can see from Figure5(b) that the motor torques of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N. The motor torques of four wheels is not following the road surface changes, whatever the road surface is wet asphalt or snowy. The required motor torques of four wheels are different while the road surfaces changes. In the switching process between various operation modes and different road surfaces, the oscillatory behavior may occur.From the results of simulation, we can find that traditional MPC cannot realize the optimal slip ratio control while EV switches on various operation modes and different road surfaces; then, it is difficult to ensure the longitudinal stability of EVs.


## 4.1.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure6. As Figure 6(a) shows, the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot achieve when the road surfaces changes. It is easy to find in Figure 6(b) that T is also kept within a stable range.Figure 6
The simulation results of MPC-MMCS on test 1. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)The wheel velocity is displayed in Figure6(c), the smoothly switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 6(d), the λi can accurately reflect the operation modes of EV.In conclusion, the MPC-MMCS can successfully recognize the operation modes of EV. Yet the optimal slip rate cannot achieve when the EV works under different road surfaces.


## 4.1.3. The Performance of MPC-OSRCS

The simulation results of MPC-OSRCS are shown in Figure7. Note that the weight coefficient λi and ξj are calculated by the fuzzy method and a modified recursive Bayes theorem. They represent the matching degree between the state of actual EV and three reference operation modes and between the actual road surfaces and reference road surfaces, respectively.Figure 7
The simulation results of MPC-OSRCS on test 1. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torque.
(a)
(b)
(c)
(d)
(e)The biggest challenges in the proposed MPC-OSRCS are the optimal slip ratio control under different operation modes and various road surfaces. As shown in Figure7(a) that the optimal wheel ratio slip can be achieved no matter what the state of actual EV, in addition, the optimal slip rate follows the road surfaces change in electric vehicle operation. While the EV has switched different road surfaces, the optimal slip ratio also can smoothly achieved. We can see from the picture that the optimal slip ratio control can be realizes by proposed MPC-OSRCS when EV switches on various operation modes and different road surfaces.Figure7(b) illustrates the wheel velocity r⋅ω. We can see from Figure 6(b) that EV runs smoothly when it switches on various operation modes and different roads surfaces. We can see from Figures 7(c) and 7(d) that the λi represents the matching degree between the state of actual EV and each operation modes, and the ξj represents the matching degree between the actual roads surfaces and the reference road surfaces. As described above, the matching degree is almost identical to operation conditions.We can see from Figure7(e) that T changes smoothly while the state of EV changing. Whether the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small and the comfortable of EV is well.Compared with the traditional MPC, the proposed MPC-OSRCS can realize the smooth switching and optimal slip ratio control whatever the running state of EV. The simulation clearly proves that the proposed MPC-OSRCS can be suitable for different operation modes and various road surfaces and all can achieved smooth switching and optimal slip ratio control.


## 4.2. Test 2

The EV also is tested on various operation modes and different road surfaces in test 2. The starting speed of EV is 85 km/h and the straight running only is considered. The state of EV is deceleration, its velocity decreases from 80 km/h to 10 km/h during0−1.75s, and then the EV keeps the speed unchanged. At this time, the state of EV is uniform speed during 1.75−3.25s. At the next time, the EV is acceleration and the velocity of EV increase from 10 km/h to 80 km/h during 3.25−5s. At the same time, the road surfaces are also changed over time. At the beginning, the EV is driving on the ice road surface. After 0.75 s, the road surfaces are changed to snow road surface. After 0.75 s, the road surfaces are changed to wet pebbles road surface. After 1 s, the road surfaces are changed to wet asphalt road surface. After 1.5 s, the road surfaces are changed to dry asphalt road surface.


### 4.2.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure8.Figure 8
The simulation results of conventional MPC on test2. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure8(a), we can see that the slip ratio is mainly affected by operation modes and the road surfaces adhesion coefficient changes can be ignored. However, the slip ratio is affected by operation modes and road surfaces, and the optimal slip ratio control cannot be achieved under different operation modes and various road surfaces. Therefore, it is proved that the conventional MPC cannot be suitable for the multimodes and the conventional MPC cannot achieve optimal slip ratio control performance in actual working.We can see from Figure8(b) that the motor torques remain unchanged. In order to be suitable for different operation modes and various road surfaces, the motor torque should change, following the state of EV. However, the motor torque of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N, whether the road surface is wet asphalt or snowy. Therefore, it is difficult to ensure the longitudinal stability of EV under traditional MPC, while EV works on various operation modes and different road surfaces. The conventional MPC cannot ensure driver comfort and longitudinal stability under this situation.From the simulation of traditional MPC, we can conclude that it is difficult to realize the optimal slip ratio control and smooth switching, while EV works on various operation modes and different road surfaces.


### 4.2.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure9.Figure 9
The simulation results of MPC-MMCS on test 2. (a) Tire longitude slip ratioκ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)As shown in Figure9(a), the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot be achieved when the road surfaces changes. It is easy to find in Figure 9(b) that T is also kept within a stable range. The wheel velocity is displayed in Figure 9(c), and the smooth switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 9(d), the λi accurately reflects the operation modes of EV.


### 4.2.3. The Performance of MPC-OSRCS

The simulation results of the proposed MPC-OSRCS are shown in Figure10.Figure 10
The simulation results of MPC-OSRCS on test 2. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torques.
(a)
(b)
(c)
(d)
(e)We can see from Figure10(a) the EV works under different road surfaces and various operation modes, and the slip ratio is controlled at the optimal slip ratio for the corresponding road surface, no matter how the state of EV changes. Compared with the traditional MPC, the slip ratio κ changes quickly when the state of EV changes, and it can be proved that the rapidity of the controller is good. The slip ratio can be effectively achieved by the MPC-OSRCS and can be suitable for various operation modes and different road surfaces.Figure10(b) illustrates the test on wheel velocity r⋅ω. It is observed that EV runs smoothly under various operation modes and different roads surfaces. We can see from Figures 10(c) and 10(d) that the identifier can quickly distinguish when the state of EV changes. As described above, the matching degree change following the state of EV changes and can precisely reflect the state of EV changes under different road surfaces and various operation modes.We can see from Figure10(e) that T changes smoothly while the state of EV change, whether it is the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small, and the comfortable of EV is well.Compared with the conventional MPC, the proposed MPC-OSRCS can identify different operation modes and various road surfaces, and optimal slip ratio control can be realized. The simulation results reveal that the proposed MPC-OSRCS can better ensure the longitudinal stability of EV under different operation modes and various road surfaces.


## 4.2.1. The Performance of Traditional MPC

The simulation results of the traditional MPC are shown in Figure8.Figure 8
The simulation results of conventional MPC on test2. (a) Tire longitude slip ratio κ. (b) Distribution of motor torques.
(a)
(b)As shown in Figure8(a), we can see that the slip ratio is mainly affected by operation modes and the road surfaces adhesion coefficient changes can be ignored. However, the slip ratio is affected by operation modes and road surfaces, and the optimal slip ratio control cannot be achieved under different operation modes and various road surfaces. Therefore, it is proved that the conventional MPC cannot be suitable for the multimodes and the conventional MPC cannot achieve optimal slip ratio control performance in actual working.We can see from Figure8(b) that the motor torques remain unchanged. In order to be suitable for different operation modes and various road surfaces, the motor torque should change, following the state of EV. However, the motor torque of four wheels Tee=1,2,3,4 are appropriate and equal to 80 N, whether the road surface is wet asphalt or snowy. Therefore, it is difficult to ensure the longitudinal stability of EV under traditional MPC, while EV works on various operation modes and different road surfaces. The conventional MPC cannot ensure driver comfort and longitudinal stability under this situation.From the simulation of traditional MPC, we can conclude that it is difficult to realize the optimal slip ratio control and smooth switching, while EV works on various operation modes and different road surfaces.


## 4.2.2. The Performance of MPC-MMCS

The simulation results of MPC-MMCS are shown in Figure9.Figure 9
The simulation results of MPC-MMCS on test 2. (a) Tire longitude slip ratioκ. (b) Distribution of motor torques. (c) Vehicle velocity. (d) Weight of operation modes.
(a)
(b)
(c)
(d)As shown in Figure9(a), the optimal slip rate follows the road surfaces change in electric vehicle operation. However, the optimal slip rate cannot be achieved when the road surfaces changes. It is easy to find in Figure 9(b) that T is also kept within a stable range. The wheel velocity is displayed in Figure 9(c), and the smooth switching can be achieved when the EV works under different operation modes and various surfaces road surfaces. In Figure 9(d), the λi accurately reflects the operation modes of EV.


## 4.2.3. The Performance of MPC-OSRCS

The simulation results of the proposed MPC-OSRCS are shown in Figure10.Figure 10
The simulation results of MPC-OSRCS on test 2. (a) Tire longitude slip ratioκ. (b) Vehicle velocity. (c) Weight of road surfaces. (d) Weight of operation modes. (e) Distribution of motor torques.
(a)
(b)
(c)
(d)
(e)We can see from Figure10(a) the EV works under different road surfaces and various operation modes, and the slip ratio is controlled at the optimal slip ratio for the corresponding road surface, no matter how the state of EV changes. Compared with the traditional MPC, the slip ratio κ changes quickly when the state of EV changes, and it can be proved that the rapidity of the controller is good. The slip ratio can be effectively achieved by the MPC-OSRCS and can be suitable for various operation modes and different road surfaces.Figure10(b) illustrates the test on wheel velocity r⋅ω. It is observed that EV runs smoothly under various operation modes and different roads surfaces. We can see from Figures 10(c) and 10(d) that the identifier can quickly distinguish when the state of EV changes. As described above, the matching degree change following the state of EV changes and can precisely reflect the state of EV changes under different road surfaces and various operation modes.We can see from Figure10(e) that T changes smoothly while the state of EV change, whether it is the change of operation modes or the road surfaces, the T is also kept within a stable range T=240. The difference between motor torques of four wheels is small, and the comfortable of EV is well.Compared with the conventional MPC, the proposed MPC-OSRCS can identify different operation modes and various road surfaces, and optimal slip ratio control can be realized. The simulation results reveal that the proposed MPC-OSRCS can better ensure the longitudinal stability of EV under different operation modes and various road surfaces.


## 5. Conclusions

Aiming at solving the problem that the traditional MPC cannot realize the optimal slip control while EV switches on various operation modes and different road surfaces, a novel MPC-OSRCS is proposed in this paper. It can not only identify different operation modes but also recognize various road surfaces. The control performance of MPC-OSRCS for EV is verified while EV works on different operation modes and various road surfaces. Simulation results demonstrate the advantage of MPC-OSRCS. Compared with the conventional MPC, the MPC-OSRCS can effectively improve longitudinal stability performance and achieve the optimal slip ratio control under various operation states and different road surfaces.



---

*Source: 8086590-2020-04-29.xml*