A Trajectory Tracking Control Strategy of 4 WIS / 4 WID Electric Vehicle with Adaptation of Driving Conditions

The Four Wheel Independent Steering/Driving (4WIS/4WID) electric vehicle has the advantage that the rotation angle and driving torque of each wheel can be independently and accurately controlled. In this paper, a trajectory tracking strategy based on the hierarchical control method is designed. In the path tracking layer, the nonlinear state feedback controller is used, and the neural network Proportion Integration Differentiation (NNPID) controller is designed to track the desired path and to obtain the desired yaw rate. By tracking the desired yaw rate and vehicle speed, the terminal sliding mode controller in vehicle dynamics control layer calculates the desired resultant tire force. In the tire force distribution layer, the multiple optimization objectives, including vehicle stability performance objective, energy-saving performance objective, and tire wear energy consumption objectives are determined and the weight coefficient is adaptive to different working conditions based on fuzzy logic theory. Finally, the wheel steering angle and driving torque of each wheel are calculated by the nonlinear three-degree-of-freedom vehicle model. Simulation results show that it realizes the adaptive control of tire force while tracking the desired trajectory, improves the stability and energy saving of the vehicle, and effectively reduces tire wear.


Introduction
Autonomous vehicles can sense the road environment through the on-board sensor system, automatically plan the traffic route and then reach the target position [1].As the major component of Intelligent Transportation System (ITS), autonomous vehicles have become the focus of research these years.One of the most fundamental issues related to autonomous cars is the path following problem, which serves as the basis for the design of control laws that enable a vehicle to follow the desired path with minimum error and behave with good stability against external disturbances [2,3].
The past decades have witnessed many researches in trajectory tracking [4][5][6].Many control structures have been proposed, which are mainly based on EPS (electronic power steering) [7], AFS (active steering system) [8], or SBW (steer by system) [9].A model predictive controller (MPC) is presented for solving the path-tracking problem of terrestrial autonomous vehicles that can better match the target criteria.A new path tracking technique called "vector pursuit" is also proposed.It was shown that the new method could be more robust, resulting in more accurate path tracking.Shin et al. pointed out the limitations of existing schemes that restrict their consideration to kinematic models and showed that it is possible to obtain an increase in performance through the use of approximate dynamical models that capture first-order effects [10][11][12].Nevertheless, due to the structural limits on the front wheel steering vehicle, the control approaches in the aforementioned literature have been proposed to just stabilize the path-following errors, and little research addressed the vehicle stability.
The 4WS (Four wheel steering) vehicle has the structure with independents of both steering positions of the front and rear wheels, which is conducive to the control of vehicle stability [13].In recent years, there have been many studies on the stability control and path tracking control of 4WS vehicles.A model-based control for fast autonomous four-wheel mobile robots on soft soils is developed.This control strategy takes into account slip and skid effects to extend the mobility over planar granular [14].In Reference [15], a nonlinear four-wheel steering (4WS) autonomous Vehicles (AVs) path tracking systems using robust adaptive sliding mode controller is researched.An automatic path-tracking controller of a four-wheel steering (4WS) vehicle based on the sliding mode control theory is proposed in Reference [16].The controller has an advantage in that the front and rear wheel steering can be decoupled at the front and rear control points.In Reference [17], the path-following control problem for four-wheel independently actuated autonomous ground vehicles is investigated through integrated control of active front-wheel Steering and direct yaw-moment control.From the above literature, it can be seen that the 4WS vehicle can solve the stability control problem well.However, as for 4WS vehicle, due to the lateral and longitudinal coupling characteristics and nonlinear characteristics of vehicle dynamics, it is difficult to guarantee the accuracy of path tracking under the complex conditions of high speed or variable speed.
The 4WIS/4WID electric vehicle is driven by in wheel motors and is designed with independent steering actuators.It enables higher flexibility in vehicle motion control by independently controlling the driving force and the rotation angle of each wheel, and has the advantages of environmentally friendliness, saving energy, and so on [18][19][20].In recent years, the research on the trajectory tracking of 4WIS/4WID electric vehicles has attracted more and more attention.In Reference [21], the lateral displacement at the virtual points near front and rear axles as new state variables is used in the path tracking controller design.Tan et al. presents two methods to enable a 4WIS 4WID vehicle to accurately follow a predefined path as well as its reference trajectories including velocity and acceleration profiles [22].In Reference [23], the projected error e p combine the lateral position error e and the heading angle error ∆φ is used to design the lateral controller.An adaptive hierarchical trajectory following the control system is presented in Reference [24] to enhance the tracking performance and improve the lateral stability of autonomous 4WID electric vehicles.These schemes enable the effective following of a desired path for an autonomous 4WIS/4WID vehicle.However, through the integrated chassis control system, the vehicle can have better performance while realizing the path tracking function.Since the four-wheel longitudinal and lateral forces of 4WIS/4WID electric vehicle are independently controllable, the performance of the vehicle can be improved by distributing the tire forces of each wheel.Therefore, the 4WIS/4WID vehicle is a good carrier to realize the function of automatic driving.Furthermore, the application of 4WIS/4WID vehicle to solve the coupling characteristics of lateral and longitudinal motion of vehicle dynamics needs further study.
The main contributions of this paper are threefold.(1) To improve the control precision, the 4WIS/4WID vehicle control structure considering the lateral and longitudinal coupling characteristics of vehicle dynamics is proposed.(2) An adaptive control method is presented for different driving conditions.Three different evaluation indices are used to distribute the tire force of each wheel in order to further improve the vehicle performance while achieving trajectory tracking.(3) The proposed control structure can regulate the lateral velocity, in a reasonable range, for vehicle stability.Furthermore, the simulation experiments are carried out to verify the control method proposed in this paper.
In this paper, a nonlinear NNPID controller is designed.The desired path of the 4WIS/4WID vehicle is tracked under the vehicle speed change condition to obtain the desired yaw rate.The Terminal sliding mode controller is designed to track the speed and the desired yaw rate, to converge the tracking error to zero in the limited time, and to ensure the robustness and global stability of the closed loop system.The desired resultant forces and desired resultant yaw moments are obtained.Considering the motor energy-saving control and tire wear performance index control, the multi-optimization objectives tire force distribution control strategy is designed to distribute the longitudinal and lateral forces of each wheel.An adaptive variable weight coefficient strategy based on vehicle speed and path curvature is designed to achieve optimal distribution of tire force.Finally, the wheel angles and driving torques are calculated by the four-wheel three-degree-of-freedom vehicle model and the tire model.
The rest of the paper is organized as follows.The nonlinear four-wheel vehicle dynamics model and vehicle trajectory mode for the optimization controller is introduced in Section 2. The hierarchical control structure and the controller design method are proposed in Section 3. In Section 4, the tire force distribution method and actuator control method are proposed.Simulation results and analysis are illustrated in Section 5, followed by the conclusion in Section 6.

Vehicle System Modeling
The picture of 4WIS/4WID electric vehicle is shown in Figure 1a.All the wheels of the vehicle can be independently steering and driving.In Figure 1b, the 3-degree of freedom (DOF) nonlinear vehicle dynamics model is applied to present vehicle dynamics features.In this model, the longitudinal degree of freedom is added to the normal two-degree-of-freedom model.Since this research mainly takes into account lateral, longitudinal, and yaw dynamics, other vehicle motions here, such as pitch, roll, and vertical motions are neglected in the presented vehicle model [25].
Appl.Sci.2019, 9, x FOR PEER REVIEW 3 of 27 stability of the closed loop system.The desired resultant forces and desired resultant yaw moments are obtained.Considering the motor energy-saving control and tire wear performance index control, the multi-optimization objectives tire force distribution control strategy is designed to distribute the longitudinal and lateral forces of each wheel.An adaptive variable weight coefficient strategy based on vehicle speed and path curvature is designed to achieve optimal distribution of tire force.Finally, the wheel angles and driving torques are calculated by the four-wheel three-degree-of-freedom vehicle model and the tire model.The rest of the paper is organized as follows.The nonlinear four-wheel vehicle dynamics model and vehicle trajectory mode for the optimization controller is introduced in Section 2. The hierarchical control structure and the controller design method are proposed in Section 3. In Section 4, the tire force distribution method and actuator control method are proposed.Simulation results and analysis are illustrated in Section 5, followed by the conclusion in Section 6.

Vehicle System Modeling
The picture of 4WIS/4WID electric vehicle is shown in Figure 1a.All the wheels of the vehicle can be independently steering and driving.In Figure 1b, the 3-degree of freedom (DOF) nonlinear vehicle dynamics model is applied to present vehicle dynamics features.In this model, the longitudinal degree of freedom is added to the normal two-degree-of-freedom model.Since this research mainly takes into account lateral, longitudinal, and yaw dynamics, other vehicle motions here, such as pitch, roll, and vertical motions are neglected in the presented vehicle model [25].In Figure 1b, CG is the vehicle's center of gravity; 1 B s wheelbase; β is the vehicle sideslip angle; and a and b are the distance from CG to front axle and rear axle, respectively.The vehicle lateral dynamics including the longitudinal, lateral, and yaw motions can be given as follows: where  denotes the mass of the vehicle;  and  are the longitudinal and lateral acceleration, respectively;  denotes the aerodynamic resistance coefficient;  denotes windward area;  denotes the rotational inertia;  is the tire-road friction coefficient;  is the gravitational acceleration; and  and  represent the longitudinal and lateral velocity, respectively. denotes In Figure 1b, CG is the vehicle's center of gravity; B 1 s wheelbase; β is the vehicle sideslip angle; and a and b are the distance from CG to front axle and rear axle, respectively.The vehicle lateral dynamics including the longitudinal, lateral, and yaw motions can be given as follows: where m denotes the mass of the vehicle; . v x and .
v y are the longitudinal and lateral acceleration, respectively; C D denotes the aerodynamic resistance coefficient; A denotes windward area; I z denotes the rotational inertia; µ is the tire-road friction coefficient; g is the gravitational acceleration; and v x and v y represent the longitudinal and lateral velocity, respectively.
. ϕ denotes the yaw rate and .. ϕ denotes the angular acceleration.ΣF x , ΣF y , and ΣM z are the longitudinal tire resultant force, lateral tire resultant force, and yawing moment.τ x , τ y , and τ z represent the external disturbance input.ΣF x , ΣF y , and ΣM z can be written as: where F xi and F yi (i = 1, 2, 3, 4) are the tire forces of each wheel in the x direction and y direction (the subscript i denotes the front left, front right, rear left, and rear right wheels).
As shown in Figure 2a,b, the tire forces of each wheel in the x direction and y direction are related with the wheel steering angle δ i (i = 1, 2, 3, 4) and longitudinal and lateral tire force F mi , F ni (i = 1, 2, 3, 4), respectively, which can be expressed by matrix form: where matrices M, N, and P are given as follows: Appl.Sci.2019, 9, x FOR PEER REVIEW 4 of 27 the yaw rate and  denotes the angular acceleration.Σ , Σ , and Σ are the longitudinal tire resultant force, lateral tire resultant force, and yawing moment. ,  , and  represent the external disturbance input.Σ , Σ , and Σ can be written as: where Fxi and Fyi (i = 1, 2, 3, 4) are the tire forces of each wheel in the x direction and y direction (the subscript i denotes the front left, front right, rear left, and rear right wheels).
As shown in Figure 2a,b, the tire forces of each wheel in the x direction and y direction are related with the wheel steering angle δi (i = 1, 2, 3, 4) and longitudinal and lateral tire force Fmi, Fni (i = 1, 2, 3, 4), respectively, which can be expressed by matrix form: where matrices M, N, and P are given as follows:

Tire Model
As the vehicle lateral dynamics possesses inherent nonlinearities, the linear tire model cannot accurately represent nonlinear characteristics of tire dynamics.Several tire models have been used in the literature to describe nonlinear characteristics of the tire force.Pacejka tire models are one of the most commonly used tire models.Pacejka tire model, also called Magic Formula, contains the combination formula of trigonometric functions to describe the tire characteristics under different working conditions [26].Should the tire operates under the assumption of pure cornering and pure longitudinal slip condition, the Pacejka tire model can be expressed as: Since the form and definition of each group of parameters are similar, the parameters B x , C x , D x , E x B y , C y , D y , E y can be abbreviated as B, C, D, E.Where parameters B, C, D, E are stiffness factors, curve shape factors, peak factors, and curvature factors, respectively.These parameters variables are the functions of tire vertical load, camber angle and tire characteristics; s is the longitudinal slip ratio; α is the tire slip angle; and F x0 and F y0 are the longitudinal and lateral tire force.
The tire model on the combined condition of lateral and longitudinal force can be expressed as follows: The tire slip angle, defined as the angle between the orientation of wheel longitudinal velocity vector and the wheel resultant velocity vector, can be written as: where v mi and v ni are the longitudinal and lateral velocities at the center of i-th wheel, respectively.v mi and v ni can be calculated by the v xi and v yi , which are wheel velocities of each wheel in the x direction and y direction.As shown in Figure 2, the matrix form is obtained as: The wheel velocities of each wheel in the x direction and y direction can be calculated by vehicle longitudinal and lateral velocities: Since the effects of load transfer are caused by lateral and longitudinal accelerations, the tire normal forces are defined as: (20) Appl.Sci.2019, 9, 168 6 of 26 where h CG is the height of mass center.The longitudinal slip ratio s describes as the proportion of sliding motion in wheel motion is defined as: where R w is wheel rolling radius.

Vehicle Trajectory Model
Only considering the plane motion of the vehicle and ignoring the lateral velocity of the vehicle, the vehicle trajectory model can be expressed as follows [27]: where .

X and .
Y are the longitudinal and lateral velocity in the global coordinate system; X and Y are the position coordinates in the global coordinate system; and ϕ is the vehicle heading angle.Combining the first two Equations in Reference (22), the constraint equation can be obtained as follows: .
The physical meaning of the constraint is that the car can only make instantaneous motion in the direction of the vehicle heading angle.

Hierarchical Control Structure
The 4WIS/4WID electric vehicle has four independent steering systems and four in-wheel motors.The chassis integration control of 4WIS/4WID electric vehicle is adopted in this paper, so as to give consideration to other performance requirements while achieving target trajectory tracking.
The study of chassis integrated control falls into centralized control structure and hierarchical control structure.The centralized control structure is based on the vehicle dynamics characteristics of the whole vehicle, and the subsystems are directly controlled by the centralized controller at the global level.The controller design is generally based on linear or nonlinear models and the design method of multi-input multi-output system (MIMO) is adopted to solve the coupling problem of vehicle longitudinal dynamics and lateral dynamics.However, the design of the centralized controller is difficult and the control algorithm is complex [28].
To solve the problem of complexity of centralized control algorithm, a hierarchical control structure is adopted to reduce the complexity of the algorithm, as shown in Figure 3.The desired path and the desired vehicle speed are planned by the path planning layer according to the surrounding environment and other information.The path tracking layer tracks the desired path considering the change of vehicle speed to obtain the desired yaw rate.The vehicle dynamics control layer calculates the longitudinal resultant force demand, lateral resultant force demand, and yaw moment demand to achieve the desired response according to the desired yaw rate and desired vehicle speed.The longitudinal tire force and the lateral tire force of each wheel are obtained from the tire force distribution layer via a certain optimization method.The actuator control layer controls the actuators of each subsystem to achieve the vehicle movement control.
considering the change of vehicle speed to obtain the desired yaw rate.The vehicle dynamics control layer calculates the longitudinal resultant force demand, lateral resultant force demand, and yaw moment demand to achieve the desired response according to the desired yaw rate and desired vehicle speed.The longitudinal tire force and the lateral tire force of each wheel are obtained from the tire force distribution layer via a certain optimization method.The actuator control layer controls the actuators of each subsystem to achieve the vehicle movement control.

Path Tracking Control Layer
Equation ( 22) can be rewritten as: where  is the input of the control model [29].The Path tracking problem can be expressed as: For any given geometric path   ,  = 0, find  =  , to make the system move along the desired path.That is to say, for any given  > 0, exist  > 0 while  > , (, ) < 0.
Define the function  = (, ), the vehicle speed is  .The first and second derivatives of function  can be written as follows: , , ( cos sin 2 sin cos ) ...
where = ,  = ,  , = ,  , = ,  , = setting  = , where  is the new control input.The controller is designed to make  →  = 0,  →  = 0.While  = 0, Equation ( 25) can be written as: where is the tangential direction of the curve  = (, ) at point (, ).The yaw rate  is controlled, while  → 0 and  → 0, the direction of the vehicle is always the tangent of the desired path.The control algorithm can be written as follows:

Path Tracking Control Layer
Equation ( 22) can be rewritten as: where u is the input of the control model [29].The Path tracking problem can be expressed as: For any given geometric path f X de f , Y de f = 0, find u = .
ϕ, to make the system move along the desired path.That is to say, for any given ε > 0, exist T > 0 while t > T, f (X, Y) < 0.
Define the function z = f (X, Y), the vehicle speed is v x .The first and second derivatives of function z can be written as follows: . where where λ is the new control input.The controller is designed to make lim 25) can be written as: where ϕ is controlled, while z → 0 and .z → 0 , the direction of the vehicle is always the tangent of the desired path.The control algorithm can be written as follows: .
In order to facilitate the design of path planning layer, the cubic function is selected as the desired path.
Equation ( 32) can be written as In Equations ( 29)- (31), Path tracking of autonomous vehicles is a non-linear control problem, and it will be disturbed by the change of vehicle speed and external environment.To improve the anti-jamming performance, robust performance and tracking performance of the system, a NNPID controller is designed to solve the above control problems.The neural network has the ability of self-adaptation and self-learning, hence the NNPID controller can effectively improve the unstable control effect caused by the change of the system structure and parameters [30].It also has a wide range of applications and fast computing speed.
The structure of the PID control system based on the BP (Back Propagation) network is shown in Figure 4.The controller consists of two parts:

•
A classic PID controller.The controlled object is directly controlled by closed-loop control, and the three parameters k p , k i , k d are adjusted online.

•
A BP neural network.According to the operating state of the system, the parameters of the PID controller are adjusted so as to achieve the optimization of a certain performance index.
The output value of the output layer neuron corresponds to three adjustable parameters k p , k i , k d of the PID controller.Through self-learning of neural network and weighted coefficient adjustment, the neural network outputs PID controller parameters corresponding to some optimal control law.The BP neural network is designed to have three layers with nodes of three, five, and one.
In order to facilitate the design of path planning layer, the cubic function is selected as the desired path.
Equation ( 32) can be written as In Equations ( 29)- (31), Path tracking of autonomous vehicles is a non-linear control problem, and it will be disturbed by the change of vehicle speed and external environment.To improve the anti-jamming performance, robust performance and tracking performance of the system, a NNPID controller is designed to solve the above control problems.The neural network has the ability of self-adaptation and self-learning, hence the NNPID controller can effectively improve the unstable control effect caused by the change of the system structure and parameters [30].It also has a wide range of applications and fast computing speed.
The structure of the PID control system based on the BP (Back Propagation) network is shown in Figure 4.The controller consists of two parts:

•
A classic PID controller.The controlled object is directly controlled by closed-loop control, and the three parameters  ,  ,  are adjusted online.

•
A BP neural network.According to the operating state of the system, the parameters of the PID controller are adjusted so as to achieve the optimization of a certain performance index.The output value of the output layer neuron corresponds to three adjustable parameters  ,  ,  of the PID controller.Through self-learning of neural network and weighted coefficient adjustment, the neural network outputs PID controller parameters corresponding to some optimal control law.The BP neural network is designed to have three layers with nodes of three, five, and one.The PID control algorithm can be expressed as: The PID control algorithm can be expressed as: where k p is the proportion coefficient, k i is the Integral coefficient, and k d is the differential coefficient.The structure of three-layer BP neural network is proposed in this paper.The neural network consists of one input variable and three output variables.There are five neurons in the hidden layer.The input of the neural network input layer is: where e is the control error.The input and output of the hidden layer are: where w i1 (2) is the weighting coefficient of the hidden layer.
The activation function of neurons in the hidden layer selects the sigmoid function of positive and negative symmetry: The input and output of the network output layer is: The output nodes in the output layer correspond to three adjustable parameters k p , k i , k d .The performance indicator function is determined as: The weight coefficient of the network is modified according to the gradient descent method, that is, the negative gradient direction of the weighted coefficient is searched and adjusted according to E(k), and an inertial term is added to make the search rapidly converge to the global minimum.
The learning algorithm of weight coefficient of network output layer can be expressed as: The learning algorithm of the weight coefficient of a network hidden layer can be expressed as: where The NNPID control algorithm is designed, as shown in Algorithm 1.

8.
w li (1) (k) = w li (1) (k − 1) The block diagram of the path tracking control layer is shown in Figure 5.The inputs of the system are the desired vehicle speed and desired trajectory function.The outputs are the vehicle yaw rate.
Appl.Sci.2019, 9, x FOR PEER REVIEW 10 of 27 The NNPID control algorithm is designed, as shown in Algorithm 1.

The NN PID Control Algorithm
The block diagram of the path tracking control layer is shown in Figure 5.The inputs of the system are the desired vehicle speed and desired trajectory function.The outputs are the vehicle yaw rate.

Vehicle Motion Control Layer
Through the path tracking layer, the required yaw rate when the vehicle tracks the desired path at the desired vehicle speed can be obtained.The yaw motion and longitudinal motion of the vehicle are generated by the tire force.In this section, the vehicle motion control layer is designed to obtain the requirements of longitudinal resultant force, lateral resultant force, and yaw resultant moment.
Since the parameter uncertainty and strong coupling characteristics of MIMO nonlinear system, the terminal sliding mode control strategy is adopted to control the vehicle motion in this paper.The selection of Terminal sliding-mode function ensures that the output tracking error converges to zero in finite time.In addition, the initial state of the system is on the sliding surface, which eliminates the reaching stage of the sliding mode control and ensures the global robustness and stability of the closed-loop system [31].
Equation ( 1) can be written as the form of state equation:

Vehicle Motion Control Layer
Through the path tracking layer, the required yaw rate when the vehicle tracks the desired path at the desired vehicle speed can be obtained.The yaw motion and longitudinal motion of the vehicle are generated by the tire force.In this section, the vehicle motion control layer is designed to obtain the requirements of longitudinal resultant force, lateral resultant force, and yaw resultant moment.
Since the parameter uncertainty and strong coupling characteristics of MIMO nonlinear system, the terminal sliding mode control strategy is adopted to control the vehicle motion in this paper.The selection of Terminal sliding-mode function ensures that the output tracking error converges to zero in finite time.In addition, the initial state of the system is on the sliding surface, which eliminates the reaching stage of the sliding mode control and ensures the global robustness and stability of the closed-loop system [31].
Equation (1) can be written as the form of state equation: where and τ x , τ y , τ z T are uncertain variables.We define the state variables ξ = [x 1 , x 2 , x 3 ] and the input variables u = [u 1 , u 2 , u 3 ].
The First-order MIMO nonlinear system is expressed as: where The uncertainty part ∆ f (X, t) satisfies the following conditions: where F(X, t) is a non-negative function.
Defining the error vector: The sliding surface equation is designed as: where Defining Φ(t) satisfies the following assumptions: Assumption 1: The function Φ i (t) can be expressed as follows: Furthermore .
The Lyapunov function can be designed as follows: where The controller is designed as: Substituting Equation (61) into Equation (60): .
According to Assumption 1 and Terminal sliding surface Equation: The initial state of the system is already on the sliding surface.The arrival stage of synovial control is eliminated, and the global robustness and stability of the closed loop system are ensured.
In order to reduce the chattering, the continuous function vector S δ is adopted to replace C T Γ C T Γ : where δ 0 , δ 1 are two positive constants.In summary, the control strategy of Terminal sliding mode controller is:

Tire Force Distribution and Actuator Control
As shown in Figure 6, the output of the vehicle motion control layer is ∑ F x , ∑ F y , ∑ F z .Since the characteristics of four-wheel independent steering and independent driving/braking, the 4WIS/4WID has at least eight controllable degrees of freedom.The function of the tire force distribution layer is to allocate the tire force reasonably to the lateral and longitudinal direction of each wheel according to the optimal performance function, the actuator constraints and the pavement attachment conditions.In this section, the optimal allocation strategy is designed, including three evaluation indexes related to vehicle performance.The importance of the three evaluation indicators is different in different working conditions.This paper proposes an adaptive weight coefficient adjustment strategy based on fuzzy logic theory in Section 4.2, which adjusts the weight coefficient online.The signal input in Figure 6 represents the vehicle state collected by the sensor as the necessary information input for the tire force distribution layer.
Appl.Sci.2019, 9, x FOR PEER REVIEW 14 of 27 4WIS/4WID has at least eight controllable degrees of freedom.The function of the tire force distribution layer is to allocate the tire force reasonably to the lateral and longitudinal direction of each wheel according to the optimal performance function, the actuator constraints and the pavement attachment conditions.In this section, the optimal allocation strategy is designed, including three evaluation indexes related to vehicle performance.The importance of the three evaluation indicators is different in different working conditions.This paper proposes an adaptive weight coefficient adjustment strategy based on fuzzy logic theory in Section 4.2, which adjusts the weight coefficient online.The signal input in Figure 6 represents the vehicle state collected by the sensor as the necessary information input for the tire force distribution layer.

Tire Force distribution curvature calculation
Actuator control Figure 6.Block diagram of tire force distribution and actuator control.

Minimize Stability Margin
Eiichi proposed the concept of the tire adhesion margin [32].It represents the proportion of the residual adhesion between the tire and the ground to the total adhesion under the current load and the current road adhesion coefficient.In this paper, the lowest tire adhesion margin of each wheel is taken as the cost function of the control target.As the cost function approaches one, it approaches

Minimize Stability Margin
Eiichi proposed the concept of the tire adhesion margin [32].It represents the proportion of the residual adhesion between the tire and the ground to the total adhesion under the current load and the current road adhesion coefficient.In this paper, the lowest tire adhesion margin of each wheel is taken as the cost function of the control target.As the cost function approaches one, it approaches the adhesion limit.

Energy-Saving Distribution Method of Longitudinal Tire Force
The drive efficiency diagram of in-wheel motor is shown in Figure 7.There is a relatively high efficiency area in the operation of the motor.At low rotate speed and small torque, the efficiency of motor is low.At the same rotate speed, the efficiency of motor increases first and then decreases with the increase of motor torque in the low and medium rotate speed range.Therefore, the distribution of longitudinal tire force should make the motor work in the high efficiency area as far as possible, which can effectively achieve the goal of energy saving. .

Minimize Stability Margin
Eiichi proposed the concept of the tire adhesion margin [32].It represents the proportion of the esidual adhesion between the tire and the ground to the total adhesion under the current load and he current road adhesion coefficient.In this paper, the lowest tire adhesion margin of each wheel is aken as the cost function of the control target.As the cost function approaches one, it approaches he adhesion limit.The drive efficiency diagram of in-wheel motor is shown in Figure 7.There is a relatively high fficiency area in the operation of the motor.At low rotate speed and small torque, the efficiency o otor is low.At the same rotate speed, the efficiency of motor increases first and then decreases with he increase of motor torque in the low and medium rotate speed range.Therefore, the distribution f longitudinal tire force should make the motor work in the high efficiency area as far as possible hich can effectively achieve the goal of energy saving.Since the motor efficiency is a function the motor torque and rotate speed, that is: Defining the cost function J 2 : Constraint 1 means that the sum of the longitudinal tire force of four wheels is the resultant longitudinal force required; the constraint 2 means that the motor speed n i is not higher than the maximum speed n max that can be provided; the constraint 3 means that the longitudinal force of the wheels is less than the maximum longitudinal force provided by the in-wheel motor, and the adhesion condition is also satisfied; the constraint 4 means that the longitudinal force of the wheels does not generate additional yaw moment to the vehicle; in constraint 5 and 6, where ∂M w ∂δ is deformation steering coefficient, M w is the torque acting on the suspension, δ amax is maximum allowable steering angle due to deformation steering.In theory, as long as constraint 4 is satisfied, the vehicle can keep running in a straight line, but the torque of the left and right wheels on the same axis is not distributed.The torque difference between the left and right wheels could generate additional yawing moment to the suspension, which could result in the deformation steering effect.Therefore, the torque difference should not be too large.

Minimize Tire Dissipation Energy
In steering, the contact surface between the tire and the ground will appear sliding movement and sliding force, so the tire will dissipate energy due to slippage.The energy dissipation rate is defined as the product of the sliding force and the sliding velocity in the slip region.Suzuki proposes a tire force control method to minimize the energy dissipation rate [33].The cost function is defined as: where v sxi is the longitudinal tire sliding velocity, v syi is the lateral tire sliding velocity, χ sxi and χ syi are, respectively, the ratio of the slippage force to the tire force in the longitudinal and lateral directions.
The following parameters can be expressed as follows: where λ and ξ s can be expressed as: where K s means the longitudinal tire force corresponding to the unit longitudinal slip ratio when α = 0, s → 0 ; K α means the lateral tire force corresponding to the unit tire slip angle when α → 0 , s = 0.
The constraint 1 means the sum of the lateral forces of each wheel is the lateral resultant force of the vehicle; the constraint 2 means the resultant yaw moment of the vehicle is produced by the longitudinal and lateral forces of each wheel.

Weight Coefficient Adaptive Adjustment Strategy
The total cost function is defined as: where g 1 , g 2 and g 3 represent the weight coefficients of the three cost functions and are used to represent the importance of different evaluation indicators.The flow diagram of the tire force optimization algorithm is shown in Figure 8.In the cost function J 1 , the vertical load of each wheel needs to be estimated in real time.The vehicle longitudinal acceleration and lateral acceleration obtained from the sensor are input into the tire load distribution module, and the vertical load of each wheel is calculated.The rotation speed of each wheel needs to be obtained in the cost function J 2 .In this paper, the wheel rotation speeds are estimated by vehicle longitudinal speed.In the cost function J 3 , the longitudinal velocity at the wheel center, the tire slip ratio, and the tire slip angle should be obtained.The above state is estimated by the observer in real time.The sequential quadratic programming method is used to solve the problem.Finally, the longitudinal force and lateral force of each wheel are F xi and F yi , respectively.The tire force optimal allocation strategy proposed in this paper considers three optimization objectives, which can improve vehicle performance comprehensively.However, there is a complex coupling relationship between the three optimization objectives and their importance is also different under different conditions.Therefore, an adaptive weight coefficient adjustment strategy based on fuzzy logic theory is proposed in this paper to adjust the weight coefficient online.In this paper, we used the type of Fuzzy based Mamdani model.An adaptive fuzzy controller is synthesized from a collection of fuzzy IF-THEN rules.The parameters of the membership functions The tire force optimal allocation strategy proposed in this paper considers three optimization objectives, which can improve vehicle performance comprehensively.However, there is a complex coupling relationship between the three optimization objectives and their importance is also different under different conditions.Therefore, an adaptive weight coefficient adjustment strategy based on fuzzy logic theory is proposed in this paper to adjust the weight coefficient online.In this paper, we used the type of Fuzzy based Mamdani model.An adaptive fuzzy controller is synthesized from a collection of fuzzy IF-THEN rules.The parameters of the membership functions characterizing the linguistic terms in the fuzzy IF-THEN rules change according to some adaptive law for the purpose of getting the needed control value [34].
Since in the course of vehicle movement, the vehicle longitudinal speed and the curvature of the driving path will directly represent the motion state of the vehicle, and will affect the function values of the three cost functions introduced in Section 4.1.Therefore, the input of the fuzzy control strategy is selected as the vehicle longitudinal speed v x and the curvature of the desired path τ.The calculation formula for τ can be expressed as: The larger the value of τ, the greater the curvature of the path.
Analyzing the degree of influence of three evaluation indexes.When the vehicle moves along a straight line or the path of small curvature, the vehicle mainly does the longitudinal motion, so the requirements for the energy saving of the motor are high, and the requirements for tire wear control and stability control are low; when the vehicle is driving on the high-curvature track, the control requirements of tire wear energy consumption and vehicle stability are high, so the control requirements of motor energy efficiency can be relaxed.Additionally, when the vehicle is driving at high speed along the bend, the control of the stability margin is the most important.With the increase of vehicle speed and path curvature, the importance of this factor increases gradually.The fuzzy subset of v x and τ is {S, M, B}, in which the S, M, B represent small, medium, and big, respectively.Additionally, the fuzzy subset of g 1 , g 2 , g 3 is {NI, I, VI}, in which the NI, I, VI represent not important, important, and very important.The fuzzy rules are designed through expert rules and written, as shown in Table 1.
Table 1.The fuzzy rules.

The Fuzzy Rules
The research domain selection for the input of the fuzzy controller is as follows: The fuzzy domain of the vehicle speed v x is [0, 30]; the fuzzy domain of the Path curvature τ is [0, 0.2]; and the fuzzy domain for the output of the fuzzy controller g 1 , g 2 , g 3 , which are the weight coefficients of cost functions, is [0, 1].
Since the three subsets of S, M and B are a relatively fuzzy definition of vehicle speed, it is considered that the decay of membership function with the change of vehicle speed is not linear.Therefore, we adopt the Gaussian-type membership function for v x , as shown in Figure 9a.For τ and g 1 , g 2 , g 3 , experience shows that when the membership function changes linearly with the parameters, the control strategy has better control effect.Hence, by utilizing the expertise-based method, the Triangular-type membership function is decided.It is adopted for τ and g 1 , g 2 , g 3 .The membership function curves of τ and g 1 , g 2 , g 3 are shown in Figure 9b,c.
The research domain selection for the input of the fuzzy controller is as follows: The fuzzy domain of the vehicle speed  is [0, 30]; the fuzzy domain of the Path curvature  is [0, 0.2]; and the fuzzy domain for the output of the fuzzy controller  ,  ,  , which are the weight coefficients of cost functions, is [0, 1].
Since the three subsets of S, M and B are a relatively fuzzy definition of vehicle speed, it is considered that the decay of membership function with the change of vehicle speed is not linear.Therefore, we adopt the Gaussian-type membership function for  , as shown in Figure 9a.For  and  ,  ,  , experience shows that when the membership function changes linearly with the parameters, the control strategy has better control effect.Hence, by utilizing the expertise-based method, the Triangular-type membership function is decided.It is adopted for  and  ,  ,  .The membership function curves of  and  ,  ,  are shown in Figure 9b,c.Surface of the fuzzy logical controller output is demonstrated in Figure 10a-c.
The research domain selection for the input of the fuzzy controller is as follows: The fuzzy domain of the vehicle speed   is [0, 30]; the fuzzy domain of the Path curvature  is [0, 0.2]; and the fuzzy domain for the output of the fuzzy controller  1 ,  2 ,  3 , which are the weight coefficients of cost functions, is [0, 1].
Since the three subsets of S, M and B are a relatively fuzzy definition of vehicle speed, it is considered that the decay of membership function with the change of vehicle speed is not linear.Therefore, we adopt the Gaussian-type membership function for   , as shown in Figure 9a.For  and  1 ,  2 ,  3 , experience shows that when the membership function changes linearly with the parameters, the control strategy has better control effect.Hence, by utilizing the expertise-based method, the Triangular-type membership function is decided.It is adopted for  and  1 ,  2 ,  3 .The membership function curves of  and  1 ,  2 ,  3 are shown in Figure 9b,c

Actuator Control Layer
As for real vehicles, the longitudinal and lateral forces of each wheel cannot be directly controlled, and the motion of vehicles is controlled by actuators.The actuators of the 4WIS/4WID electric vehicle include independent steering system and in-wheel motor.By controlling the actuators, the longitudinal and lateral tire forces required for each wheel can be obtained.
Through the tire model under combined working conditions introduced in Section 2.2, the required tire slip angle and tire slip rate under the expected tire lateral force and longitudinal force are calculated.The Newton Raphson algorithm is adopted in this paper.Equations ( 10)-( 16) can be written as: Defining the follow functions: Setting the initial value as: According to the dynamic characteristics of tires, the tire may work in a linear or saturated area in the same lateral tire force.The tire slip angles corresponding to these two conditions are different.Similarly, different longitudinal slip rates may also be obtained in the same tire longitudinal force.Therefore, in order to make the tire work in the linear region, the initial value is set in this paper to ensure the direction of iteration convergence.The initial value is defined as S 0 = 0.001, α 0 = 0.1.
The recurrence formula is written as: The computational procedure of Newton Raphson algorithm is shown in Algorithm 2: Algorithm 2. Newton Raphson algorithm.
1. define initial value x (0) , accuracy value ε The desired longitudinal slip rate s ide f and the desired tire slip angle α i under the expected tire longitudinal force F xi and the lateral force F yi of each wheel can be obtained through the Iterative calculation.The desired wheel steering angle of each wheel is δ i , which can be calculated by Equations ( 17)- (19).
The desired longitudinal slip rate is tracked by controlling drive torque of in-wheel motor.The in-wheel motor adopts DC motor.The desired wheel angular velocity can be expressed as: Designing the PID controller.The controlled quantity is the torque of in-wheel motor M i .The error is defined as e = ω a − ω ide f .
The block diagram of actuator control layer is shown in Figure 11.
The block diagram of actuator control layer is shown in Figure 11.

Experimental Results
In order to appraise the performances of the proposed control architecture, a series of numerical simulation experiments are carried out.In this section, a control system with the decoupled longitudinal and lateral is applied to contrast with the presented control approach.Linear quadratic regulator (LQR) control is one of the most representative control methods since it is a relatively mature control method.Hence, in this paper, LQR control is applied as a comparison to further verify the proposed control strategy.The optimal solution of LQR control can be written as a unified analytic expression, and the closed-loop optimal control system can be reconstructed simply by using the state linear feedback control law, which takes into account multiple performance indicators.Therefore, the lateral vehicle motion controller is selected as a traditional LQR controller.In terms of longitudinal control, the four-wheel drive is distributed evenly.An 11 degrees-of-freedom high fidelity vehicle model is established, which is used for simulation.The vehicle parameters are presented in Table 2.

Experimental Results
In order to appraise the performances of the proposed control architecture, a series of numerical simulation experiments are carried out.In this section, a control system with the decoupled longitudinal and lateral is applied to contrast with the presented control approach.Linear quadratic regulator (LQR) control is one of the most representative control methods since it is a relatively mature control method.Hence, in this paper, LQR control is applied as a comparison to further verify the proposed control strategy.The optimal solution of LQR control can be written as a unified analytic expression, and the closed-loop optimal control system can be reconstructed simply by using the state linear feedback control law, which takes into account multiple performance indicators.Therefore, the lateral vehicle motion controller is selected as a traditional LQR controller.In terms of longitudinal control, the four-wheel drive is distributed evenly.An 11 degrees-of-freedom high fidelity vehicle model is established, which is used for simulation.The vehicle parameters are presented in Table 2.The reference path profile used in simulation is a polynomial curve with small curvature radius.Figure 12 shows the desired acceleration and velocity of the test vehicle.The vehicle starts the run at initial velocity 20 m/s.As shown in Figure 12a,b, the vehicle begins to accelerate from zero, then it runs at a constant speed of 18.2 m/s.The above condition demonstrates that the vehicle runs at a case of high vehicle and small path curvature, which can verify the control effect.The reference path profile used in simulation is a polynomial curve with small curvature radius.Figure 12 shows the desired acceleration and velocity of the test vehicle.The vehicle starts the run at initial velocity 20 m/s.As shown in Figure 12a,b, the vehicle begins to accelerate from zero, then it runs at a constant speed of 18.2 m/s.The above condition demonstrates that the vehicle runs at a case of high vehicle and small path curvature, which can verify the control effect.Desired tire lateral and longitudinal forces from tire force distribution layer for the four wheels during the automatic driving are shown in Figure 13.It is clear that when the autonomous electric vehicles runs at an accelerated velocity in the road with small curvature, the desired longitudinal Desired tire lateral and longitudinal forces from tire force distribution layer for the four wheels during the automatic driving are shown in Figure 13.It is clear that when the autonomous electric vehicles runs at an accelerated velocity in the road with small curvature, the desired longitudinal force of each wheel is about 500 kN, which is devoted to ensure the vehicle motion in longitudinal direction.Simulation results show that the presented control distribution algorithm in the multi-objective hierarchical control architectures can be used to obtain preferable distribution results.Desired tire lateral and longitudinal forces from tire force distribution layer for the four wheels during the automatic driving are shown in Figure 13.It is clear that when the autonomous electric vehicles runs at an accelerated velocity in the road with small curvature, the desired longitudinal force of each wheel is about 500 kN, which is devoted to ensure the vehicle motion in longitudinal direction.Simulation results show that the presented control distribution algorithm in the multi-objective hierarchical control architectures can be used to obtain preferable distribution results.The dynamic behaviors of the proposed control architectures and the traditional LQR control system are compared in Figures 14-16.As shown in Figure 14, the real trajectories are plotted together with the simulated trajectory, and the vehicle follows the path reference with an acceptable performance.It is clear that, compared to the LQR controller, the proposed architectures has less tracking error.The dynamic behaviors of the proposed control architectures and the traditional LQR control system are compared in Figures 14-16.As shown in Figure 14, the real trajectories are plotted together with the simulated trajectory, and the vehicle follows the path reference with an acceptable performance.It is clear that, compared to the LQR controller, the proposed architectures has less tracking error.The dynamic behaviors of the proposed control architectures and the traditional LQR control system are compared in Figures 14-16.As shown in Figure 14, the real trajectories are plotted together with the simulated trajectory, and the vehicle follows the path reference with an acceptable performance.It is clear that, compared to the LQR controller, the proposed architectures has less tracking error.Figure 16a shows the simulation results of changeable weight coefficients.When the vehicle accelerates on a road with a small curvature radius, the coefficient  is larger than  ,  , which means the energy saving control of in-wheel motors is more important.With the increasing of vehicle speed and curvature, the weight coefficient  and  gradually increase to ensure the stability of the vehicle and reduce tire wear.After 6 s, due to the small curvature radius, the weight coefficients  ,  ,  tends to be stable, where  is larger than  and  .Figure 16b-d indicate that, compared with LQR controller, the proposed control architecture can effectively reduce the stability margin, power consumption, and tire energy dissipation.Among them, the effect of motor energy saving control is the most obvious, and the other two have also been optimized, which coincides with the curve of weight coefficients.

Medium Speed and Large Path Curvature
The reference path profile used in simulation is a polynomial curve with small curvature radius.Figure 17 shows the desired velocity of the test vehicle.Figure 15 shows that the vehicle slip angle are controlled in a smaller range by the proposed control architecture than traditional LQR control system.The results show that the proposed control architecture can restrict the vehicle slip angle within 0.8 deg, whereas the traditional LQR control system make the vehicle slip angle be bounded by 2.5 deg.It means that the vehicle has better maneuverability.
Figure 16a shows the simulation results of changeable weight coefficients.When the vehicle accelerates on a road with a small curvature radius, the coefficient g 1 is larger than g 2 , g 3 , which means the energy saving control of in-wheel motors is more important.With the increasing of vehicle speed and curvature, the weight coefficient g 2 and g 3 gradually increase to ensure the stability of the vehicle and reduce tire wear.After 6 s, due to the small curvature radius, the weight coefficients g 1 , g 2 , g 3 tends to be stable, where g 1 is larger than g 2 and g 3 .
Figure 16b-d indicate that, compared with LQR controller, the proposed control architecture can effectively reduce the stability margin, power consumption, and tire energy dissipation.Among them, the effect of motor energy saving control is the most obvious, and the other two have also been optimized, which coincides with the curve of weight coefficients.

Medium Speed and Large Path Curvature
The reference path profile used in simulation is a polynomial curve with small curvature radius.Figure 17 shows the desired velocity of the test vehicle.Figure 16b-d indicate that, compared with LQR controller, the proposed control architecture can effectively reduce the stability margin, power consumption, and tire energy dissipation.Among them, the effect of motor energy saving control is the most obvious, and the other two have also been optimized, which coincides with the curve of weight coefficients.

Medium Speed and Large Path Curvature
The reference path profile used in simulation is a polynomial curve with small curvature radius.Figure 17 shows the desired velocity of the test vehicle.Figure 18 shows the desired tire lateral and longitudinal forces calculated by proposed control architecture.The desired lateral force of each tire is about 1400 kN, which is devoted to ensure the vehicles automatic steering in the curve.Figure 19 shows that the reference path can be tracked effectively.It can be seen from Figure 20 that the vehicle slip angle is reduced by the proposed control method, which means better stability performance.Figure 18 shows the desired tire lateral and longitudinal forces calculated by proposed control architecture.The desired lateral force of each tire is about 1400 kN, which is devoted to ensure the vehicles automatic steering in the curve.Figure 19 shows that the reference path can be tracked effectively.It can be seen from Figure 20 that the vehicle slip angle is reduced by the proposed control method, which means better stability performance.Figure 21a shows the simulation results of changeable weight coefficients in the case of medium vehicle speed and large path curvature.In the path tracking process, there are time periods when the path curvature radius is large, so the weight coefficient g 2 , g 3 is large and g 1 is relatively small.However, the coefficient g 1 is relatively large in the time period of low speed and small path curvature radius.Figure 21b-d show the optimization effect of the proposed control architecture, which further proves the above conclusion.Figure 21a shows the simulation results of changeable weight coefficients in the case of medium vehicle speed and large path curvature.In the path tracking process, there are time periods when the path curvature radius is large, so the weight coefficient  ,  is large and  is relatively small.However, the coefficient  is relatively large in the time period of low speed and small path curvature radius.Figure 21b-d show the optimization effect of the proposed control architecture, which further proves the above conclusion.

Conclusions
This paper proposed a path tracking control strategy of 4WIS/4WID electric vehicle.A hierarchical control method is designed, which is aimed to improve the properties of automatic driving for trajectory tracking.The control strategies proposed in this paper includes a path-tracking

Conclusions
This paper proposed a path tracking control strategy of 4WIS/4WID electric vehicle.A hierarchical control method is designed, which is aimed to improve the properties of automatic driving for trajectory tracking.The control strategies proposed in this paper includes a path-tracking layer, vehicle dynamics control layer, tire force distribution layer, and actuator control layer.Path tracking layer calculates the desired yaw rate under variable vehicle speed; vehicle dynamics control layer calculates the desired resultant tire force to track the desired yaw rate and vehicle speed; tire force distribution layer considers multiple optimization objectives, designs the weight coefficient adaptive adjustment strategy, and calculates the longitudinal and lateral forces of each wheel; the wheel steering angle and driving torque of each wheel are calculated by actuator control layer.Some numerical simulations are carried out.The simulation results show that the control method can not only realize the accurate path tracking, but also realize the adaptive optimization of tire force under different working conditions, thus improving the stability and economy of the vehicle.

Figure 2 .
Figure 2. Schematic of tire force.(a) Schematic of front left wheel tire force; and (b) Schematic of rear left wheel tire force.

Figure 2 .
Figure 2. Schematic of tire force.(a) Schematic of front left wheel tire force; and (b) Schematic of rear left wheel tire force.

Figure 6 .
Figure 6.Block diagram of tire force distribution and actuator control.

2 .
Energy-Saving Distribution Method of Longitudinal Tire Force

Figure 7 .
Figure 7. Drive efficiency diagram of in wheel motor.

Figure 7 .
Figure 7. Drive efficiency diagram of in wheel motor.

Figure 8 .
Figure 8. Block diagram of tire force optimal distribution.

Figure 9 .Figure 10 .
Figure 9. Membership function curves (a) Membership function curves of  ; (b) membership function curves of ; and (c) membership function curves of  ,  ,  .Surface of the fuzzy logical controller output is demonstrated in Figure10a-c.

Figure 9 .
Figure 9. Membership function curves (a) Membership function curves of v x ; (b) membership function curves of τ; and (c) membership function curves of g 1 , g 1 , g 3 .

Figure 9 .Figure 10 .
Figure 9. Membership function curves (a) Membership function curves of   ; (b) membership function curves of ; and (c) membership function curves of  1 ,  1 ,  3 .Surface of the fuzzy logical controller output is demonstrated in Figure 10a-c.

Figure 10 .
Figure 10.Surface of the fuzzy logical controller output: (a) Surface of fuzzy logical controller output g 1 ; (b) surface of fuzzy logical controller output g 2 ; and (c) surface of fuzzy logical controller output g 3 .
controller.The controlled quantity is the torque of in-wheel motor  .The error is defined as  =  −  .

Figure 11 .
Figure 11.Block diagram of actuator control layer.

Figure 11 .
Figure 11.Block diagram of actuator control layer.

1 .
High Speed and Small Path Curvature

Figure 12 .
Figure 12.Desired acceleration and velocity in the high-speed case.(a) Desired acceleration.(b) Desired velocity.

Figure 12 .
Figure 12.Desired acceleration and velocity in the high-speed case.(a) Desired acceleration.(b) Desired velocity.

Figure 12 .
Figure 12.Desired acceleration and velocity in the high-speed case.(a) Desired acceleration.(b) Desired velocity.

Figure 13 .
Figure 13.Desired tire forces provided by the tire Force Distribution layer: (a) Front left tire force; (b) front right tire force; (c) rear left tire force; and (d) rear Right tire force.

Figure 14 .
Figure 14.Path tracking in the x-y plane.

Figure 15
Figure15shows that the vehicle slip angle are controlled in a smaller range by the proposed control architecture than traditional LQR control system.The results show that the proposed control architecture can restrict the vehicle slip angle within 0.8 deg, whereas the traditional LQR control system make the vehicle slip angle be bounded by 2.5 deg.It means that the vehicle has better maneuverability.

Figure 13 .
Figure 13.Desired tire forces provided by the tire Force Distribution layer: (a) Front left tire force; (b) front right tire force; (c) rear left tire force; and (d) rear Right tire force.

Figure 13 .
Figure 13.Desired tire forces provided by the tire Force Distribution layer: (a) Front left tire force; (b) front right tire force; (c) rear left tire force; and (d) rear Right tire force.

Figure 14 .
Figure 14.Path tracking in the x-y plane.

Figure 15
Figure15shows that the vehicle slip angle are controlled in a smaller range by the proposed control architecture than traditional LQR control system.The results show that the proposed control architecture can restrict the vehicle slip angle within 0.8 deg, whereas the traditional LQR control system make the vehicle slip angle be bounded by 2.5 deg.It means that the vehicle has better maneuverability.

Figure 14 .
Figure 14.Path tracking in the x-y plane.

Figure 15
Figure15shows that the vehicle slip angle are controlled in a smaller range by the proposed control architecture than traditional LQR control system.The results show that the proposed control architecture can restrict the vehicle slip angle within 0.8 deg, whereas the traditional LQR control system make the vehicle slip angle be bounded by 2.5 deg.It means that the vehicle has better maneuverability.

Figure 16 .
Figure 16.Experience results in high speed and small curvature case: (a) Weight coefficients; (b) stability margin; (c) electrical efficiency; and (d) tire energy dissipation.

Figure 16 .
Figure 16.Experience results in high speed and small curvature case: (a) Weight coefficients; (b) stability margin; (c) electrical efficiency; and (d) tire energy dissipation.

Figure 16 .
Figure 16.Experience results in high speed and small curvature case: (a) Weight coefficients; (b) stability margin; (c) electrical efficiency; and (d) tire energy dissipation.

Figure 17 .
Figure 17.Desired velocity in medium speed case.Figure 17.Desired velocity in medium speed case.

Figure 17 .
Figure 17.Desired velocity in medium speed case.Figure 17.Desired velocity in medium speed case.

Figure 18 .
Figure 18.Desired tire forces from the tire Force Distribution layer: (a) Front left tire force; (b) front right tire force; (c) rear left tire force; and (d) rear Right tire force.

Figure 18 .Figure 18 .
Figure 18.Desired tire forces from the tire Force Distribution layer: (a) Front left tire force; (b) front right tire force; (c) rear left tire force; and (d) rear Right tire force.

Figure 19 .
Figure 19.Path tracking in the x-y plane.

Figure 19 .
Figure 19.Path tracking in the x-y plane.

Figure 21 .
Figure 21.Experience results in medium speed and large curvature case: (a) Weight coefficients; (b) stability margin; (c) electrical efficiency; and (d) tire energy dissipation.

Figure 21 .
Figure 21.Experience results in medium speed and large curvature case: (a) Weight coefficients; (b) stability margin; (c) electrical efficiency; and (d) tire energy dissipation.

Table 2 .
The parameters of vehicle.

Table 2 .
The parameters of vehicle.