Handling Enhancement of Autonomous Emergency Steering for Reduced Road Friction Using Steering and Differential Braking

: Steering has more potential than braking to prevent rear-end collisions at higher relative velocities. A path tracking controller based on multi-input multi-output (MIMO) model predictive control (MPC) is proposed to enhance the handling performance of autonomous emergency steering in this paper. A six-state MIMO bicycle model including actuator dynamics of steering and differential braking is used for model prediction. Two control inputs are front wheel steering angle and direct yaw moment. Two model outputs are lateral displacement and heading angle. According to the work load ratios at four wheels, control allocation is used to determine the optimal braking force distribution to prevent tire force saturation. The performance of a single-input single-output (SISO) MPC that uses only steering angle control to track the lateral displacement of the desired path is employed to benchmark the performance of the proposed algorithm. Simulation results show that both SISO MPC and MIMO MPC can track the path on nominal road surface with high road friction coefﬁcient of 0.9. For a road surface with medium road friction coefﬁcient of 0.7, the SISO MPC is unable to track the path and loses directional stability. However, the MIMO MPC can still track the path and demonstrate robust path tracking and handling enhancement against model uncertainty due to reduced road friction.


Introduction
According to Traffic Safety Facts [1], rear-end collisions account for 18% of the fatal collisions between motor vehicles. Singh et al. [2] pointed out that driver error is the primary cause for most accidents. The main factor is the recognition error caused by driver inattention, internal and external distractions, and inadequate surveillance. Advanced driver assistance systems (ADAS) can assist the driver to prevent possible collisions with the obstacle in front of the vehicle via braking or steering. Bosch [3] reported that the autonomous emergency braking (AEB) system would reduce 30% of rear-end collisions. Alfred et al. [4] showed that the necessary minimum distance to avoid a collision by steering is smaller than that by braking for maneuvers with large relative velocities. Especially for a suddenly appearing obstacle, AEB is not able to avoid a collision due to insufficient braking distance. The European New Car Assessment Programme introduced safety ratings for AEB in 2014, and has announced that [5] emergency steering assist (ESA) and autonomous emergency steering system (AES) tests will start in 2020 and 2022, respectively. The goal of the tests is to deliver a further significant reduction in crashes and casualties.

Related Work
Both ESA and AES can make a sudden lane change or evasive maneuver via steering actions to avoid a collision. The main difference between ESA and AES is the trigger mechanism. When the ESA [6][7][8][9] detects the risk of collision and the steering torque from the driver, the system provides additional overlay torques to adjust the steering angle to track the evasive path without collision or directional instability. However, Shimizu et al. [10] showed that 37.9% of drivers did not take any avoidance maneuvers when encountering the approaching danger. Therefore, collision avoidance cannot be guaranteed for inattentive maneuvers with late steering actions from the driver. Unlike ESA, the AES [11][12][13][14][15] can steer automatically when the detected risk of collision is over the designed threshold.
Three types of planning approaches, search-based, optimization-based, and parameterbased, are commonly used for generating a desired path to avoid the collision. Search-based methods [16][17][18] are designed to generate an action set of multiple drivable trajectories, allowing an adjacent behavior planner to pick the most appropriate action for the global state in the scene. Optimization-based methods [19][20][21] plan the path from a mapping run and offline compute an optimal trajectory around the track using the trajectory optimization module to maximize the progress rate, and minimize slip angle cost and input rate. However, the first and second types of planning approaches are not easy to implement on production control units due to high computation loads. Parameter-based methods [7] including circular, sine/cosine, s-shape, polynomial, and trapezoidal-acceleration-profile (TAP) trajectories are designed by specifying the desired lateral displacement for a lane change maneuver that can also be employed for collision avoidance. Among the parameterbased methods, only TAP can consider transition time and ride comfort at the same time [22].
Two types of path tracking controllers, model-free and model-based, are commonly used for tracking an evasive path. If the dynamic model of the plant to be controlled is not needed for the controller design, it is classified as the model-free approach. The conventional model-free controller calculates the steering angle command based on the lateral displacement error and the heading error at the previewed point [23][24][25]. Machine learning can also be classified as the model-free approach. A neural network model can be trained to generate the steering action using the dataset of human drivers [26] or deep reinforcement learning [27]. However, the major challenge is the difficulty in proving the safety issue of these methods [28]. As for model-based approaches, the single-input single-output (SISO) model predictive control (MPC) is often used to obtain the vehicle steering angle based on a four-state bicycle model, which is used to predict the future lateral displacement of a vehicle [29][30][31].
Differential braking via electric hydraulic braking (EHB) can be used to generate the desired direct yaw moment (DYM) to improve path tracking performance. Falcone et al. [32] proposed a multi-input multi-output (MIMO) MPC controller that uses four-wheel longitudinal slip ratios and front wheel steering angle as the control input to minimize errors in heading angle, yaw rate and lateral displacement. Their simulation results showed that path error can be reduced by 50%. However, intensive computations are required for online linearization and optimization due to the complexities of the prediction model with nonlinear tire dynamics. Hajiloo et al. [33] proposed a MIMO MPC controller to determine the optimal front lateral force and the yaw moment to improve the yaw rate responsiveness for collision avoidance without decreasing the stability. Lin et al. [34] proposed a MIMO MPC controller for AES using front wheel steering angle and direct yaw-moment control (DYC) to enhance the performance of path tracking. The pseudo inverse control allocation (CA) is used to allocate the four-wheel brake forces for generating the desired DYM. Their simulation results showed that the path tracking performance of the MIMO MPC is better than that of the conventional SISO MPC. However, the actuator dynamics of steering and differential braking are not discussed in their proposed approach. The control performance might be diminished or control lost if the actuator dynamics are not fast enough to be negligible.
Electric power steering (EPS) can be used to generate the desired front wheel steering angle. However, the actuator dynamics might cause the actual steering angle to lag slightly behind the steering angle command. Chen et al. [35] proposed a MPC controller by considering EPS actuator dynamics as a first-order model to enhance the tracking performance and ride comfort. Kim et al. [36] proposed a MPC controller by considering the EPS as a second-order model. Simulation results showed that tracking performance is improved by incorporating steering actuator dynamics in the prediction model. Thornton et al. [37] proposed a MPC controller with steering input to evaluate the average and maximum yaw rates of different actuator dynamics. Their experiment results showed that the lumped first-order dynamics can meet design objectives with low computational complexity.

Main Contributions
The proposed path-tracking controller for AES is an extension of the previous work [34] in which the actuator dynamics were neglected in the vehicle model.

•
The dynamic responses of the vehicle model in the previous work might deviate from those of the actual vehicle model with actuator dynamics for severe maneuvers. The EPS and EHB dynamics are integrated with the four-state bicycle model to enhance the accuracy of the prediction model. • A diagonal matrix with components of tire workload ratios is used to replace the identity weighting matrix of CA to prevent saturating tire forces while generating the desired DYM. • Different road friction coefficients are used for sensitivity analysis to evaluate the robustness and handling enhancement of the proposed algorithm at different vehicle velocities via simulation studies.

Organization
The remainder of this paper is organized as follows: system modeling for controller design is briefly introduced in Section 2, followed by the proposed algorithm in Section 3. Simulation results are presented in Section 4. Finally, conclusions are made in Section 5.

Actuator Dynamics
First-order actuator dynamics of EPS and EHB as shown below are included in the plant model for MPC. .
where δ f ,r and M z,r are the control commands of the front wheel steering angle and DYM, respectively; δ f and M z are the responses of front wheel steering angle and DYM, respectively; τ EPS and τ EHB are the time constants of EPS and EHB, respectively. τ EPS is set to be 0.125 s [37] and τ EHB is set to be 0.1 s [38] in this paper.

Bicycle Model
The six-state bicycle model is obtained by integrating the actuator dynamics with the four-state bicycle model that is linearized using small angle approximation and constant longitudinal velocity [39]. The state-space representation can be expressed as follows.
where x = y v y ψ r δ f M z T is the state vector, u = δ f ,r M z,r T is the input vector, and y = y ψ T is the output vector; v y and r are the lateral velocity and yaw rate of the vehicle, respectively; y and ψ are the lateral displacement and heading angle of the vehicle, respectively; A m , B m and C m are system matrix, input matrix and output matrix as shown below.
where m is the vehicle mass; v x is the vehicle longitudinal velocity; C α f and C αr are the cornering stiffness of the front and rear axles, respectively; l f and l r are the distances from the C.G. to the front and rear axles, respectively; I z is the yaw-plane rotational inertia of vehicle. The ranks of controllability matrix C and observability matrix O of the six-state bicycle model can be calculated as follows.
Since both C and O are full rank of 6, the proposed model is fully controllable and observable.

Model Discretization
A D-Class Sedan in CarSim is selected as the host vehicle to obtain the parameter values of m, a, b, and I z . Nominal values of C α f and C αr are identified on a road surface with the nominal friction coefficient of 0.9 by minimizing the path errors between the trajectory responses of CarSim and the bicycle model. The vehicle parameters used in CarSim are shown in Table 1.
For the discrete controller design, the discretized plant model can be expressed as follows.
where k denotes for the kth sample. The sample time T s is selected to be 0.04 s in this paper; u k is the discretized input vector; Φ m , Γ m and H m can be expressed as follows.

Controller Design
The block diagram of the proposed algorithm is shown in Figure 1. The path planning module generates the reference path R s using the trapezoidal acceleration profile [7]. The six-state bicycle model is used to predict the future path. The PGC calculation module calculates the PGC index from R s . The optimizer module is used to compute the best set of future control commands. The first two elements of the set are used as the desired control commands of front wheel steering angle and DYM. The CA module distributes brake pressure at the four wheels to generate the desired DYM. After the AES is triggered, the driver is disengaged from the vehicle until the driver's intent to override the system is detected.

Model Predictive Controller
The predicted state vectors in the prediction horizon N p with the future control commands in the control horizon N c can be expressed as follows.
where the initial condition of the state vector x k,i can be expressed as follows.
where the initial lateral displacement and heading angle are zero in the body coordinates of the vehicle; r and δ f can be measured using yaw rate and steering angle sensors, respectively;v y is the estimated lateral velocity that can be estimated using localization-based [40] or inertial-measurement-unit-based [41] approaches;M z can be estimated using estimated logitudinal tire forces as follows.M whereF x is the estimated longitudinal tire force that can be estimated using a closed-loop disturbance estimator [42]; B CA is the effective matrix as shown below.
where d is the track width of the vehicle. The predicted output vectors in the prediction horizon N p can be expressed as follows.
From Equations (11) and (15), the augmented prediction model can be formulated as follows.
where U p is the input vector; Y p is the output vector; F and G are the system and input matrices, respectively. They are defined as follows.
The optimizer calculates the optimal U p to minimize the cost function J p as shown below.
where R s is the reference input in the prediction horizon; Q and R are weighting matrices. They are defined as follows.
where y r is the desired lateral displacement; ψ r is the desired heading angle; Q 1 is the weighting of lateral displacement error; Q 2 is the weighting of heading angle error; R 1 is the weighting of the steering angle; R 2 is the weighting of the DYM. These weightings are designed using Bryson's rule [43] as shown below.
where y e,max is the maximum allowable lateral displacement error; ψ e,max is the maximum allowable heading angle error; δ f ,max is the maximum allowable steering angle; M z,max is the maximum allowable DYM. The optimal U p can be obtained as follows.
The current optimal control input u k,i is equal to the vector that consists of the first and second components of U p as follows.

Evasive Path Generation
The desired evasive path is shown in Figure 2. A Trapezoidal Acceleration Profile (TAP) is used to design the desired global lateral displacement of the evasive path [7]. The lateral acceleration can be obtained by specifying the maximum allowable a y,max , the maximum allowable jerk J y,max , and the evasive lateral displacement y eva . t 1 and t 2 of the TAP can be expressed as follows [7].
where a y,max =μg; g is the gravity;μ is the estimated road friction coefficient that is assumed to be obtained from other approaches [44][45][46]. The desired global lateral position Y d of the evasive path can then be obtained as follows by integrating the global lateral acceleration .. Y d twice [7].
where 5]. By assuming the vehicle longitudinal velocity v x is constant along the evasive path, the desired global longitudinal velocity . X d of the evasive path can then be obtained as follows.
where v x,0 is the initial longitudinal velocity of vehicle; . Y d is the global lateral velocity that can be obtained by integrating .. (31). We can then obtain the desired global X d by integrating . X d . The station number of the evasive path S d that is the distance along the designed evasive path can be defined as follows.

Y d in Equation
where j denotes the jth sample point of the evasive path. The initial values of X d , Y d , .
Y d , and S d are set to be zero. The desired evasive path can then be expressed as functions of the station number as follows.
The station number in the prediction horizon can be defined as follows with a constant v x .
S p,k+i+1 = S p,k+i + v x iT s , i ∈ 0, 1, · · · , N p − 1 where the initial station number S p,k is the current station number S c . The desired evasive path in the preview horizon can then be expressed as follows.
The desired evasive path in the body coordinates of the vehicle can be calculated using coordinate transformation as shown below.
where Ψ c is the current global heading angle of vehicle; X c and Y c are the current global position of the vehicle, respectively; x r and y r are the local coordinates of the evasive path. X c , Y c , Ψ c and S c are assumed to be available from the localization system. The desired lateral displacement of the evasive path can be expressed as y r = f r (x r ) in the body coordinates of the vehicle. The heading angle reference, which is the first derivative of the desired path, can then be approximated as follows.

Control Allocation
Weighted pseudo inverse with constraints [47] is used to obtain the distributed brake forces at four wheels. The DYM command M z,r can be expressed as follows.
where F x is the tire force vector that consists of four tire forces.
where the subscripts fl, rl, fr, and rr denote quantities at front-left, rear-left, front-right, and rear-right tires, respectively. The pseudo inverse solution is a two-norm solution to the CA problem that can be formulated as follows.
subject to: where C CA is the offset vector used to represent off-nominal condition with one or more tire forces; F x,min is the vector used to represent the minimum allowable tire forces. The maximum allowable tire forces are set to be 0 such that only braking forces are used. The optimal distribution can then be obtained as follows.
where B # CA is the pseudo inverse of B CA . The C CA is zero in the normal condition, so (44) can be simplified as follows.
where the diagonal weighting matrix W CA is designed with elements equal to the estimated tire workload ratios as follows [42] to prevent saturating tire forces. The larger the ratio is, the closer to saturation the corresponding tire force is. Thus, tire forces that are needed for generating the desired DYM are distributed to tires with lower tire work load ratios.
whereρ is the estimated tire workload ratio of each tire that can be obtained as follows [42].
whereF x ,F y andF z are the estimated tire forces at longitudinal, lateral and normal directions, respectively. If the ith element of F x exceeds the boundary values, the corresponding ith element in C CA is set to be the opposite number of the boundary value and the corresponding ith element in B CA is set to be zero. The adjusted vectors B CA,a and C CA,a can be obtained using the Algorithm 1 as shown below.

Algorithm 1 Saturation Prevention
After obtaining the adjusted vectors, we can replace B CA and C CA with B CA,a and C CA,a , respectively, and obtain the adjusted tire force vector using Equations (44) and (45).
The four-wheel brake pressure command P b,r can be expressed as follows.
where G b, f and G b,r are the brake gains of front and rear wheels, respectively; r w is the wheel radius; diag is a function that creates a square diagonal matrix with the input elements on the main diagonal. The EHB module is used to generate four-wheel brake pressure P b .

Simulation Results
The proposed algorithm was evaluated using CarSim in MATLAB/Simulink. Figure 3 shows the test scenario used in this paper. There was one moving vehicle and one stationary vehicle in front of the host vehicle. After the front moving vehicle changed lanes at the last moment [48], the stationary vehicle was detected as a suddenly appearing obstacle by the CarSim Sensor module of the host vehicle at 5 m ahead the last point to steer (LPS) [4], which was defined as the last location to avoid a collision by steering. The longitudinal vehicle velocity was set at 80 km/h which is the highest longitudinal velocity for testing ESA. Since the stationary vehicle suddenly appeared in front of the host vehicle, the last point to brake (LPB), which was defined as the last location to avoid the collision by braking, was already behind the host vehicle. A single lane-change path generated using TAP was used as the evasive path to avoid the collision. When the host vehicle arrived at the LPS, the AES was trigged automatically. The maximum angular rate of the front wheel steering angle was limited to 42 • /s due to the EPS according to the requirement of the fishhook test [49]. The front wheel steering angle was limited between −35 • and 35 • due to the pinion-rack configuration. Simulation results of the proposed MIMO MPC were compared with those of the SISO MPC [35]. For the first test scenario, the nominal road friction coefficient was 0.9 which was used to obtain the nominal cornering stiffness of the prediction model in MPC. The second test scenario was used to test the robustness of the proposed algorithm with a medium road friction coefficient of 0.7 which is 22% less than the nominal value; i.e., the actual cornering stiffness should be smaller than the nominal values used in the prediction model.

Comparison of Prediction Models
An open-loop test was used to compare the responses of three bicycle models. The structures of these models are listed in Table 2. The four-state bicycle model [34] has two inputs and two outputs that are the same as those of the proposed six-state bicycle model in (3). However, the actuator dynamics of EPS and EHB were not considered for the four-state bicycle. The five-state bicycle model that includes the EPS actuator dynamics has one steering angle control input and the lateral displacement output. The DYM input was not considered in the five-state bicycle model. Dynamic responses of these models are shown in Figure 4. The step steering angle input was applied at the beginning of simulation, and the DYM step input was applied at 1 s. CarSim denotes the dynamic responses of CarSim with two inputs of steering angle and DYM. Before 1 s, the main difference between the four-state model and other models was the transient response due the lack of actuator dynamics. The steady-state responses of three bicycle models were almost the same. The responses of the five-state and the six-state bicycle models were similar to those of CarSim. After the DYM was applied at 1 s, only the responses of the six-state bicycle model were close to those of CarSim. Since the six-state model with actuator dynamics can achieve both transient and steady-state responses closest to those of CarSim, it was used as the prediction model in the proposed algorithm.  The simulation results with accurately estimated road-friction coefficient of 0.9, which is used to obtain the maximum allowable lateral acceleration for designing the desired evasive path, are shown in Figures 5-7. The AES is assumed to be triggered at 3 s. SISO MPC and MIMO MPC can track the desired evasive path and complete the lane change maneuver.
Path tracking responses are shown in Figure 5. SISO MPC and MIMO MPC have similar responses of global lateral displacement. However, their global heading angle responses are less similar to each other. As can be seen from the vehicle dynamic responses in Figure 6, the values of y e of MIMO MPC are more negative than those of SISO MPC between 3.1 and 3.7 s. Thus, the lane change response of MIMO MPC is faster than that of SISO MPC at the beginning of lane change due to the additional DYM input. Meanwhile, the heading angle errors of MIMO MPC are smaller than those of SISO MPC during the lane change maneuver. Figure 7 shows control commands and actuator responses. At the beginning of AES, MIMO MPC and SISO MPC have the same steering angle responses due to the angular rate limitation and the actuator dynamics of EPS. However, MIMO MPC has an additional DYM input, which can accelerate the response of r and reduce the maximum steering angle to be smaller than that of SISO MPC. Since the heading angle is the integration of yaw rate, faster response of r can lead to smaller ψ e . The comparison of the β − r phase plane plots is shown in Figure 8. MIMO MPC has better convergences due to a smaller circling area around the origin. This means that MIMO MPC can reduce β in order to enhance handling performance and directional stability. Although DYM is generated using braking forces, the vehicle longitudinal velocity of MIMO MPC is only reduced slightly compared to that of SISO MPC.

Nominal Road Friction Coefficient of 0.7
The simulation results with a medium road friction coefficient of 0.7, which is 22% less than the nominal value, are shown in Figures 9-13. Assuming the AES is not aware of the overestimation of the road friction coefficient, the nominal value of 0.9 is still used to design the evasive path and the nominal cornering stiffnesses are still used in the prediction model.
Path tracking responses are shown in Figure 9. The AES is assumed to be triggered at 3 s. SISO MPC cannot complete the lane change maneuver and loses the directional stability. Conversely, MIMO MPC can track the desired evasive path and complete the lane change maneuver successfully. Similar to the analysis of Figure 8 in Section 4.1, we found that MIMO MPC has a lane change response faster than that of SISO MPC at the beginning of lane change, as shown in Figure 10. Meanwhile, the peak heading angle error of MIMO MPC is smaller than that of SISO MPC at about 4 s due to minimizing y e and ψ e simutaneously via both steering and DYM. From the yaw rate response, we can observe that the vehicle with SISO MPC starts to spin continuously after 4 s.
As can be seen from Figure 11, SISO MPC requires steering angles much larger than those of MIMO MPC after 4 s. Although SISO MPC tries to correct the steering angle as fast as possible, the steering angle response is not fast enough due to the angular rate limitation and actuator dynamics. Figure 12 shows that the side-slip angle of SISO MPC increases rapidly with large yaw rate after passing β = 0. It also confirms the directional instability of SISO MPC. On the other hand, the β − r phase plane plot of MIMO MPC converges to the origin after completing the lane change maneuver. Figures 13 and 14 show the responses of tire slip angle and lateral tire forces at four wheels. The dashed lines indicate the thresholds of the tire slip angles to saturate lateral tire forces. The front lateral tire forces of SISO MPC are found to be saturated at 4.1 s. Even though the steering angle of SISO MPC in Figure 11 is still increasing after 4 s, the corresponding lateral tire forces are not increased due to saturation as shown in Figure 14.
Since the lateral tire force is no longer increased with the increment of the tire slip angle, the prediction model used in MPC is not accurate anymore. Thus, SISO MPC cannot effectively control the vehicle to track the desired evasive path and prevent directional instability. Although the rear tire slip angles of MIMO MPC are saturated around 4.5 s, the front tire slip angles are not saturated. Thus, the front tires can still provide more lateral tire forces via increasing the steering angle input if necessary. Meanwhile, MIMO MPC has additional DYM which can improve the responses of r and prevent saturating the reference command of the steering angle. Figure 15 shows the responses of tire workload ratios. The CA distributes the longitudinal tire forces as shown in Figure 16 to the wheels with lower tire workload ratios in order to generate the required DYM without saturating the tire forces; thus, demonstrating better, more robust control performance and handling enhancement against model uncertainty due to reduced road friction coefficient.

Sensitivity Analysis
The estimated road friction coefficient be used not only for evasive path generation, but also for obtaining the equivalent front and rear cornering stiffnesses of the bicycle model used in MPC via a look-up table. However, the estimation might not be accurate due to the lack of enough excitation. Different road friction coefficients deviated from the estimated nominal value of 0.9 are used for sensitivity analysis to evaluate the robustness of the proposed algorithm. The minimum distance D L [50] between the host vehicle and the stationary vehicle to be avoided as shown in Figure 17 can be defined as follows.
where P h (x 1 , y 1 ) is the front center point of host vehicle; P s (x 2 , y 2 ) is the front center point of stationary vehicle; W h is the width of host vehicle; W s is the width of stationary vehicle. D L at different road friction coefficients are shown in Table 3. The collision can be avoided for all simulated scenarios. The vehicle sideslip angle can be used to evaluate the yaw stability of the vehicle. The root-mean-square (rms) and maximum (max) values of sideslip angle responses at different road friction coefficients are shown in Table 4. Those values are obtained using the responses from the beginning of control at 3-10 s. The more the road friction coefficient deviates from the nominal value of 0.9, the larger the rms and max values of sideslip angles are. MIMO can stabilize the vehicle in the new target lane for all scenarios. However, SISO exhibits large lateral responses as shown in Figure 10    In order to evaluate the robustness of the proposed algorithm at different velocities, the longitudinal velocity of host vehicle is increased to 120 km/h in Figure 18. The suddenly appearing obstacle is assumed to be detected at 7.5 m ahead of the LPS. D L at different road friction coefficients is shown in Table 5. Similar to the results in Table 3, the collision can be avoided for all simulated scenarios. The rms and max values of sideslip angle responses at different road friction coefficients are shown in Table 6. Similar to the results in Table 4, the more the road friction coefficient deviates from the nominal value of 0.9, the larger the rms and max values of sideslip angles are. MIMO can stabilize the vehicle in the new target lane for all scenarios. However, the rms and max values of sideslip angle suddenly increase for µ ≤ 0.65 which are indications of yaw instabilities. Thus, MIMO demonstrates better robustness of yaw stability and handling performance for reduced road friction coefficients even at a much higher longitudinal vehicle velocity of 120 km/h.

Conclusions
MIMO MPC was employed to design the path tracking controller for autonomous emergency steering using steering and differential braking. The prediction model was based on a six-state road following bicycle model including the actuator dynamics of EPS and EHB. For nominal road friction, the path tracking responses of SISO MPC (using steering only) and MIMO MPC (using both steering and differential braking) are very similar, but the directivity of MIMO MPC can be improved by DYM. The usage of steering angle can be reduced, and the handling performance of MIMO MPC is better than SISO MPC. For medium road friction, SISO MPC is unable to track the path and loses directional stability. However, MIMO MPC can still reduce the maximum required steering angle and accelerate the response of r via DYM. In order to prevent tire force saturation, the CA can distribute the longitudinal tire forces to wheels with lower tire work load ratios for generating the desired DYM. Thus, MIMO MPC can achieve robust path tracking and handling enhancement against model uncertainty due to road friction reduced by 22% and 33% at 80 and 120 km/h, respectively. In this paper, the handling performance was mainly discussed for avoiding the collision with a stationary vehicle via lateral dynamics control. For a slow moving vehicle in front of the host vehicle, collision avoidance via AES should be further investigated to demonstrate the effectiveness the proposed algorithm. Since the lateral tire forces at front steering wheels are affected by load transfer that can be created via acceleration or braking, coordinated longitudinal and lateral dynamics control might further enhance the handling performance for collision avoidance.