Integrated Path Following and Lateral Stability Control of Distributed Drive Autonomous Unmanned Vehicle

: Intelligentization is the development trend of the future automobile industry. Intelligentiza-tion requires that the dynamic control of the vehicle can complete the trajectory tracking according to the trajectory output of the decision planning the driving state of the vehicle and ensure the driving safety and stability of the vehicle. However, trajectory limit planning and harsh road conditions caused by emergencies will increase the difficulty of trajectory tracking and stability control of unmanned vehicles. In view of the above problems, this paper studies the trajectory tracking and stability control of distributed drive unmanned vehicles. This paper applies a hierarchical control framework. Firstly, in the upper controller, an adaptive prediction time linear quadratic regulator (APT LQR) path following algorithm is proposed to acquire the desired front-wheel-steering angle considering the dynamic stability performance of the tires. The lateral stability of the DDAUV is determined based on the phase plane, and the sliding surface, in the improved sliding mode control (SMC), is further dynamically adjusted to obtain the desired additional yaw moment for coordinating the path following and lateral stability. Then, in the lower controller, considering the slip and the working load of four tires, a comprehensive cost function is established to reasonably distribute the driving torque of four in-wheel motors (IWMs) for producing the desired additional yaw moment. Finally, the proposed control algorithm is verified by the hardware-in-the-loop (HIL) experiment platform. The results show the path following and lateral stability can be coordinated effectively under different driving conditions.


Introduction
The autonomous unmanned vehicle (AUV) plays an essential role in an intelligent traffic system (ITS).It can liberate the driver's control of the steering wheel and improve the safety and fluidity of the traffic system [1][2][3].The distributed drive autonomous unmanned vehicle (DDAUV) is provided with four independent in-wheel motors (IWMs) to drive wheels, which makes the DDAUV have a drive redundancy [4,5].
Path following is an important research field of the DDAUV, which has been studied by many scholars.Considering the dynamic characteristics of tires, an autonomous front-wheel-steering system with the coupling control of displacement and force achieves good path following performance and road feeling through a robust control of the yaw rate [6].In reference [7], sliding mode control (SMC) is utilized to deal with nonlinearities and uncertainties due to the complicated vehicle model, and an adaptive fuzzy SMC is employed to reduce the jitter to enhance the path following capability.Falcone et al. propose a control method using nonlinear model predictive control (NMPC) for path following [8], in which the current optimal front-wheel-steering angle can be obtained by status variables and control variables of the DDAUV combined with the prediction model [9].But the complex nonlinear model and the process of rolling optimization will cause a heavy computational burden for the vehicle controller.The nonlinear vehicle model is linearized to improve computational efficiency [10].However, in [11], many experiments show that this linear dynamics MPC method is suitable for path following at a low velocity.To solve this singularity at high velocity, a MPC method with switching following error is proposed [12], in which the controller can keep a small following error at a relatively high velocity.However, this control method fails to guarantee the lateral stability of DDAUV at this moment of switching.In addition, a linear quadratic regulator (LQR) with an adaptive Q matrix is used in [13,14] to overcome some of the shortcomings of the linear time-varying MPC (LTV-MPC) [15], which improves the computational efficiency of the algorithm, while the adaptive Q matrix is applied to effectively solve problems of disturbance generation.However, the proposed LQR controller does not have a predictive function to reduce the drastic motion of the steering wheels for a future moment, which may lead to a poor lateral stability [15,16].Fitri Yakub demonstrates that the MPC controller has a prediction function for a future moment.The control of lateral stability and path following accuracy can be comparable to the LTV-MPC when a prediction module is introduced into the LQR controller [17].In the above studies, scholars have paid more attention on steering maneuverability and path following accuracy.However, lateral stability is a critical factor in path following control, which should be investigated to guarantee driving safety.
For stability control, some scholars have studied vehicle stability by dealing with driving redundancy.References [18][19][20] apply a hierarchical architecture for stability control, where the upper controller solves for the additional yaw moment through the MPC, SMC or LQR method, and the lower controller is utilized to calculate the torque for IWMs under multiple constraints.Chen et al. design the status observers for estimating tire forces and vehicle velocity under different driving conditions [21].The literature, in recent years, has started to investigate multiple actuators coupling and path following integrated control strategies to improve the lateral dynamic control.Reference [22] proposes a novel sliding mode control (SMC) method with Takagi-Sugeno (T-S) fuzzy rules.Xie et al. applies the LTV-MPC to obtain the steering angle.The weighting coefficient of following deviation on the sliding mode surface (SMS) is employed to obtain the additional yaw moment that can be produced through torque vector control.In reference [23], the SMC is used for coordinated control between path following and lateral stability, but the weight coefficient is not dynamically adjusted along with the change in the vehicle status.Reference [24] designs the objective function without considering the minimum slip of tires.Actuator failure is considered in references [10,25], so the path following model and torque allocation model are reconstructed.A monolithic control algorithm is utilized in reference [26] to achieve the purpose of integrated control based on the LTV-MPC method with the steering angle of steering wheels and torque of four IWMs as control variables.
Lateral control of the DDAUV is important for safe driving.The coordination control of path following and lateral stability should be emphasized due to the characteristics of vehicle nonlinearity and time-varying tire parameters.
In this paper, the coordinated control algorithm of path following and lateral stability for the DDAUV is proposed.The main contributions of this paper are given as follows: (1) The traditional LQR algorithm can solve the steady-state error problem of steering trajectory tracking by feedforward control, but it does not encompass roll optimization, and the control sequence within the control step can not be solved according to the predicted state deviation sequence, which makes the vehicle steering become radical.
To solve this problem, this paper adds an adaptive prediction module based on LQR (Linear Quadratic Regulator) trajectory tracking control and judges the steering steady state according to the dynamic characteristics of the tire, so as to adjust the prediction time to make the steering more stable; (2) A stability compensation controller is designed by the additional yaw moment calculated by an improved SMC controller.The phase plane diagram of the sideslip angle and the sideslip angular velocity are used to judge the current stability state of the vehicle, so as to adaptively adjust the weight ratio of the sliding mode surface in the sliding mode control to coordinate the path following and lateral stability of the DDAUV; (3) Considering the slip and the working load of four tires, a comprehensive cost function is established and solved by the Karush-Kuhn-Tuckert (KKT) condition to allocate torque for IWMs.Moreover, the proposed algorithm is verified by the hardware-inthe-loop (HIL) experiments.
The rest of the paper is organized as follows: Section 2 introduces the framework of the proposed coordinated control algorithm and vehicle system model.Section 3 discusses the design of autonomous steering controller based on the APT LQR.Section 4 presents the design of stability compensation controller.In Section 5, HIL experiments are applied to verify the proposed coordinated control algorithm.Finally, a summary conclusion is given in Section 6.

Algorithm Framework
Figure 1 shows the hierarchical control architecture of the proposed coordinated control algorithm.The upper controller is composed of the APT LQR path following controller and lateral stability compensation controller.The lower controller is composed of a torque distribution controller and a calculation module of total driving torque.In APT LQR path following controller, the adaptive prediction time (APT) module based on tires dynamics adjusts the prediction time dynamically, and the APT LQR controller with prediction and feedforward functions is proposed to calculate the steering angle of autonomous front steering.In the stability compensation controller, the phase plane of side-slip angle is utilized to dynamically determine the stability of the DDAUV coefficient and adjust the weight coefficient of SMS.The improved SMC is employed to coordinate the additional yaw moment.According to the calculated total driving force and additional yaw moment, the torque distribution controller distributes the torque reasonably for four independent IWMs considering the minimum slip and working load of four tires, so that four wheels have enough stability margin to prevent the risk of extreme operation, which can develop a good closed-loop system.Furthermore, the 2-DOF vehicle model is used to calculate the reference yaw rate and reference side-slip angle.The calculation module of driving torque is applied to obtain the desired total driving torque to reach the target velocity for the DDAUV.

Vehicle Dynamic Model
A 2-degree-of-freedom (2-DOF) four-wheel vehicle model is utilized to describe the dynamic characteristics of the DDAUV in Figure 2. The basic parameters and variable symbols of the DDAUV are listed in Table 1.

Vehicle Dynamic Model
A 2-degree-of-freedom (2-DOF) four-wheel vehicle model is utilized to describe the dynamic characteristics of the DDAUV in Figure 2. The basic parameters and variable symbols of the DDAUV are listed in Table 1.

Vehicle Dynamic Model
A 2-degree-of-freedom (2-DOF) four-wheel vehicle model is utilized to describe the dynamic characteristics of the DDAUV in Figure 2. The basic parameters and variable symbols of the DDAUV are listed in Table 1.The wheels are numbered 1, 2, 3 and 4, which represent the front left wheel, front right wheel, rear left wheel and rear right wheel, respectively.The lateral dynamic differential equation of the DDAUV is expressed [27]: where m is equipment quality, v x is the longitudinal velocity, v y is the lateral velocity of the DDAUV, .φ denotes yaw rate, and .. φ is the yaw angular acceleration.F xi (i = 1, 2, 3, 4) is the longitudinal force of the i-th wheel, F yi (i = 1, 2, 3, 4) is the lateral force of the i-th wheel, δ f denotes the steering angle of the front wheels, I z is rotational inertia, a is the distance mass center to front wheels, b is the distance mass center to rear wheels, l f and l r are front wheel width and rear wheel width, respectively.
where ∆M z is additional yaw moment.In the 2-DOF vehicle model, the lateral force could be regarded as a linear relationship with tire side-slip angle (Deng et al., 2021) [18].
where α i is tire side-slip angle of the i-th wheel, and C αi is lateral tire stiffness of the i-th wheel.α i could be described as: Considering that the front-wheel angle δ f is relatively small, let sinδ f = 0, cosδ f = 1.Combining equations from Equation (1) to Equation ( 4), the lateral dynamic differential equation of the DDAUV model can be described as:

Path Following Model
The lateral deviation ed and heading deviation e φ can be employed to evaluate the path following performance when completing the unmanned task.The lateral deviation is defined as the deviation of the distance from the real position to the reference position.The heading deviation is the angle deviation from vehicle direction to tangential direction at the desire position [28,29].As is shown in Figure 3.The reference position P means the closest point on the desired path from the center of vehicle.
World Electr.Veh.J. 2024, 15, x FOR PEER REVIEW 6 of 19 Therefore, the vehicle path following model can be described as:

Design of LQR Path following Controller
The linear quadratic regulator uses a linear feedback control law to form a closedloop optimal control system, and its control object is a state space equation expressed in a linear form.In the trajectory tracking control, the state space equation required for LQR trajectory tracking control is constructed by using the dynamic differential equations of vehicle lateral and yaw degrees of freedom and the trajectory tracking model.In this paper, the position prediction module is added on the basis of LQR trajectory tracking, and the steering intensity is judged according to the dynamic characteristics of the tire, and the prediction time is adaptively adjusted to make the steering more stable.
In an autonomous steering controller, the additional yaw moment is unnecessary to be considered in the active steering control.∆Mz is removed from Equation ( 1), and it can When the vehicle is following the desired path, heading deviation can be expressed as [14]: where φ is real heading angle of the DDAUV, and φ ref is reference heading angle.
From the geometric relationship shown in Figure 3, it can be found that the change rate of lateral deviation meets the following equation: Therefore, the vehicle path following model can be described as: The linear quadratic regulator uses a linear feedback control law to form a closed-loop optimal control system, and its control object is a state space equation expressed in a linear form.In the trajectory tracking control, the state space equation required for LQR trajectory tracking control is constructed by using the dynamic differential equations of vehicle lateral and yaw degrees of freedom and the trajectory tracking model.In this paper, the position prediction module is added on the basis of LQR trajectory tracking, and the steering intensity is judged according to the dynamic characteristics of the tire, and the prediction time is adaptively adjusted to make the steering more stable.
In an autonomous steering controller, the additional yaw moment is unnecessary to be considered in the active steering control.∆M z is removed from Equation (1), and it can be transformed into: Combining Equations ( 6)-( 9), the path following model can be further expressed as: In general, the curvature of the road changes relatively smooth.In Equation (10), .φ re f can be ignored.Equation (10) can be rewritten as the form of state space: To reduce the deviation in the state variables e d , .
e φ , the LQR controller is used to make the state variables as small as possible by solving the current control variable U. Before solving that, as inspired by the literature [14], an objective function of the system is designed as: where Q lq represents the weight matrices for the state error, and R lq represents the weight matrices for the control variable error in LQR status feedback control module.A large Q lq means that the system places more weight on the rapid attenuation of state deviation X, and a large R lq means the system attaches more importance to the reduction in control variable U.
The optimal control variable can be expressed as: where is the coefficient matrix of the LQR status feedback control module.P is a constant matrix, which can be calculated by solving the Riccati equation [2].
To eliminate the steady-state error, the state feedforward control module is introduced.Therefore, the control variable can be calculated by: where δ f and δ q are calculated by feedback control module and feedforward control module.
Equation ( 14) is substituted into Equation ( 11), so .X = 0 X = 0 can be made when the system is stable.δ q can be expressed as:

LQR Path Following Controller with Adaptive Prediction Time
Traditional LQR without a prediction module fails to reduce the over-drastic motion of the steering wheels according the road condition of the future period [21,30,31].In Figure 4a, the lateral deviation and heading deviation are large, so the front steering wheel angle requires to be turned quickly to follow the path.However, the DDAUV is moving towards the correct path.Good path tracking and vehicle stability will be reached even not to adjust the steering angle.In Figure 4b, the lateral deviation and heading deviation of the DDAUV from the projection point are both small, but the vehicle needs to change the front wheel angle to adapt to the future path.
In this paper, APT LQR with adaptive prediction function is designed for path following, in which the prediction module is used to emulate the driving habits of humans by replacing the current vehicle position with the prediction position.That can be explained as follows.
When the prediction time is short, the deviation between the prediction position and the reference position is large, and the track tracking control pays more attention to the track tracking accuracy.However, when the prediction time is longer and the deviation between the predicted position and the reference position is small, the control is more focused on steering stationarity.Therefore, the prediction time should not be constant, but should be automatically adjusted according to the current stability state of the vehicle.When the stability state is good, the prediction time should be shorter to pay attention to the trajectory tracking accuracy, but when the stability state is poor, the prediction time should be longer to pay attention to the steering smoothness.
where x p and y p are the information of the prediction position, φ p is the predicted heading angle, and T s means prediction time.Moreover, the velocity and yaw rate of the DDAUV are both considered constants when the prediction time is short enough.ure 4a, the lateral deviation and heading deviation are large, so the front steering wheel angle requires to be turned quickly to follow the path.However, the DDAUV is moving towards the correct path.Good path tracking and vehicle stability will be reached even not to adjust the steering angle.In Figure 4b, the lateral deviation and heading deviation of the DDAUV from the projection point are both small, but the vehicle needs to change the front wheel angle to adapt to the future path.In this paper, APT LQR with adaptive prediction function is designed for path following, in which the prediction module is used to emulate the driving habits of humans by replacing the current vehicle position with the prediction position.That can be explained as follows.
When the prediction time is short, the deviation between the prediction position and the reference position is large, and the track tracking control pays more attention to the track tracking accuracy.However, when the prediction time is longer and the deviation between the predicted position and the reference position is small, the control is more focused on steering stationarity.Therefore, the prediction time should not be constant, but should be automatically adjusted according to the current stability state of the vehicle.When the stability state is good, the prediction time should be shorter to pay attention to the trajectory tracking accuracy, but when the stability state is poor, the prediction time should be longer to pay attention to the steering smoothness.T s should be dynamically adjusted along with the state of the DDAUV.When the stability state is good, the short T s is design to enhance the path following ability.But when the stability state is poor, Ts should be increased to ensure the lateral ability.References [24,26,32] proved that the stability state of the DDAUV and tire dynamic characteristics are closely related.The tire side-slip characteristic is shown in Figure 5.The blue area represents the stable region, yellow areas represent the critical region, and orange areas denote the instable region.


where t1 is the based prediction time from tests and looking-up table, cs is a coefficient to limit the prediction time, θi denotes the instability coefficient of the i-th tire, θ can be seen as a function of θi, αU is the upper boundary of critical instability, and αL is the lower boundary of critical instability.

Description of Stability Compensation Controller
Sliding mode control is a kind of variable structure control system.The control structure of the system can change according to the set rules according to the current state of the system.It has the advantages of a fast response speed and strong robustness.Therefore, SMC is widely used in the control field of complex systems such as nonlinearity and uncertainty.When the controlled system performs sliding mode control in any initial state, the control process includes two stages.One is the reaching mode, that is, under the action of the control law, where the system state error moves from the initial point to the designed sliding mode surface within a certain time.The second is the sliding mode, that is, the system state error slides along the sliding mode surface to the flat under the control law.However, due to the influence of system inertia, uncertain disturbance, and control lag, the chattering phenomenon will inevitably occur in the sliding process.
As the DDAUV is an overdrive redundant system, the four IWMs can generate the indicated torque independently [33,34].A stability compensation controller is designed to Therefore, the prediction time T s can be dynamically adjusted according to following equations.
where t 1 is the based prediction time from tests and looking-up table, c s is a coefficient to limit the prediction time, θ i denotes the instability coefficient of the i-th tire, θ can be seen as a function of θ i , α U is the upper boundary of critical instability, and α L is the lower boundary of critical instability.

Description of Stability Compensation Controller
Sliding mode control is a kind of variable structure control system.The control structure of the system can change according to the set rules according to the current state of the system.It has the advantages of a fast response speed and strong robustness.Therefore, SMC is widely used in the control field of complex systems such as nonlinearity and uncertainty.When the controlled system performs sliding mode control in any initial state, the control process includes two stages.One is the reaching mode, that is, under the action of the control law, where the system state error moves from the initial point to the designed sliding mode surface within a certain time.The second is the sliding mode, that is, the system state error slides along the sliding mode surface to the flat under the control law.However, due to the influence of system inertia, uncertain disturbance, and control lag, the chattering phenomenon will inevitably occur in the sliding process.
As the DDAUV is an overdrive redundant system, the four IWMs can generate the indicated torque independently [33,34].A stability compensation controller is designed to improve lateral stability, then the minimum slip and working load of four tires are used as the control objectives to ensure that the vehicle can maintain a balance between maneuverability and lateral stability [35].
The yaw rate and side-slip angle of the DDAUV are applied to describe the lateral stability [26,36].Therefore, the comprehensive error e can be designed as: φ + τe β (20) where e .φ is error of the yaw rate, e β is error of side-slip angle, and τ is stability weight coefficient.
As the linear SMS fails to guarantee that the initial state is located on the SMS [37], the integral term is introduced, and the SMS is shown as follows: where c is a positive constant.The large c means that it takes short time to reach the stability state.But the large c will cause severe jitter of the SMC.In this paper, c is designed as 0.6 through repeated tests.
The first order derivative of Equation ( 21) can be shown as: .
Choosing a suitable convergence law will improve the dynamic quality of the improved SMC, which can suppress jitter of the system.In this paper, the saturation function is used as the convergence law to suppress jitter further.
.. φ in Equation ( 5) is substituted into Equation ( 22), and letting .s = −K v sat(s) The control variable ∆M z can be obtained as: sat(s) is the saturation function, expressed as: , and where s 0 is the positive parameter of saturation boundary.

The Adaptive Adjustment of Stability Weight Coefficient
The stability weight coefficient τ needs to be adjusted dynamically according to the stability state of the DDAUV, and the phase plane is applied to analyze the lateral stability state [24,26].When the stability of the DDAUV motion system is favorable, the small τ makes the torque distribution mainly focus on the path following.On the contrary, the large τ makes it mainly focus on lateral stability when the stability is weak.In the scene, with a velocity v x 60 km/h, a front-wheel angle δ f 3 deg, and a road adhesion coefficient µ 0.5, the phase plane is plotted in Figure 6.In Figure 6, the balance point is the convergence point of the vehicle motion in the stability state.The stability boundary is the boundary line between the stable region and the unstable region.The stability boundary can be described by: where β is the side-slip angle, • β is the side-slip angle rate, B 1 and B 2 are the parameters of stability boundary corresponding to the road adhesion coefficient µ, e w denotes system disturbance, which may be ignored to simplify the problem due to its small value, and the relationship between B 1 , B 2 and µ can refer to reference [18].

Design of Torque Distribution Controller
The purpose of the torque distribution controller is to distribute total driving torque to produce the desired additional yaw moment.The DDAUV in the stable region can converge to the balance point by its own regulation, but it should be stabilized by applying the additional yaw moment in the instable region.In Figure 6, the indicator d, which means the deviation from the state point of the DDAUV to the lower boundary of stability, is introduced to describe the stability state.A small d leads to poor stability, for which a large stability weight factor is required.A large d implies a strong self-regulation ability.Therefore, τ can be designed as a small value to enhance the ability of path following.However, τ needs to be designed no less than 0.3 to ensure the emergency capacity of the DDAUV.τ can be calculated by: where W denotes the distance between two boundaries of stability, d > 0 indicates the state point above the lower boundary, and d < 0 indicates the state below the lower boundary.

Design of Torque Distribution Controller
The purpose of the torque distribution controller is to distribute total driving torque to produce the desired additional yaw moment.
where R is the tire working radius, T i (i = 1, 2, 3, 4) is the driving torque, and T need is the total driving torque needed to reach the target velocity, which can be obtained through the vehicle driving equation.
Torque distribution coefficients h and l are introduced to further distribute the torque of the front and rear wheels as shown in Equation (27).
where h is the distribution coefficient of left IWMs, l is the distribution coefficient of right IWMs.The comprehensive cost function includes the following two objectives.
The tire utilization adhesion coefficient means the proportion of tire force to maximum tire force provided by the road, which is used to describe the tire working load [32,38].The tire utilization adhesion coefficient and the first objective function can be expressed as: The second objective function is designed to control driving torques and make the longitudinal slip and lateral slip have the minimum value.The second objective function can be described as: where s xi is slip ratio of the longitudinal direction.
Before integrating the two control objectives, it is necessary to normalize the slip power of loss to the range of [0, 1].At the same time, the DDAUV are inevitably constrained by ground adhesion, motion characteristics, and motor working capacity.The problem of torque vector control can be described as the optimization involving h and l, and the comprehensive cost function with many constraints can be shown following: (31) where Γ is weight coefficient.f is the cost function involving h and l.After testing by many simulations and HIL experiments, Γ = 0.4 is advisable, F ximax is the longitudinal adhesion, and T imax is the peak torque the IWM can provide.It is obvious that the problem is a convex nonlinear quadratic programming (QP) problem.Therefore, KKT method can be used to solve the equation [39].

HIL Experiment Results and Analysis
The HIL experiment is carried out to test the effectiveness in a real controller as shown in Figure 7.The HIL experiment platform mainly consists of the main I/O module, the power supply module, the rapid prototyping controller, and the real-time simulator.The software mainly includes the lower software and the upper control strategy.The lower software includes communication protocols and diagnosis modules, etc.During the experiment, the DDAUV model in CarSim is compiled into dynamic chain file and downloaded to the real-time simulator.The proposed control strategy is compiled and downloaded to the real controller.The communication between the controller and real-time simulator is defined in NI Veristand software 2015.

DLC Manoeuvre on Low Adhesion Road
The experiment is set up to perform the double lane change (DLC) maneuv the road adhesion coefficient µ = 0.5 and the velocity 60 km/h.Three different contr egies are designed for comparative analysis."Compensation control" denotes compensation controller with adaptive stability weight coefficient is introduced APT LQR autonomous steering controller."Autonomous steering" describes the A autonomous steering controller without the compensation controller."Stability means the compensation controller with a fixed stability weight coefficient τ = 1.indicates the experiment results under DLC maneuver on a low adhesion road.
Figure 8a shows the path following of the DDAUV with three different contr egies.The path following accuracy in autonomous steering is more favorable tha stability compensation control.And the path following accuracy with compensat trol is more favorable than that in stability control.The path following accuracy is

DLC Manoeuvre on Low Adhesion Road
The experiment is set up to perform the double lane change (DLC) maneuver with the road adhesion coefficient µ = 0.5 and the velocity 60 km/h.Three different control strategies are designed for comparative analysis."Compensation control" denotes that the compensation controller with adaptive stability weight coefficient is introduced into the APT LQR autonomous steering controller."Autonomous steering" describes the APT LQR autonomous steering controller without the compensation controller."Stability control" means the compensation controller with a fixed stability weight coefficient τ = 1. Figure 8 indicates the experiment results under DLC maneuver on a low adhesion road.
Figure 8a shows the path following of the DDAUV with three different control strategies.The path following accuracy in autonomous steering is more favorable than that in stability compensation control.And the path following accuracy with compensation control is more favorable than that in stability control.The path following accuracy is slightly decreased under compensation control than autonomous steering control, but the lateral stability can be significantly improved, which explains the possibility of stability compensation control while performing path following.Figure 8b shows the torque distribution coefficients in stability compensation control.The driving torque of four IWMs with the stability compensation is indicated in Figure 8c.The torque between left IWMs and right IWMs may be different to ensure stability when the DDAUV plans to change lane at high velocity, and even some IWMs generate reverse torque to reach the target additional yaw moment and improve wheels stability margins.In Figure 8d,e, the comparison results about the yaw rate and side-slip angle show that the stability of the DDAUV is improved widely under stability compensation control.The phase trajectory portrait for the three control strategies is shown in Figure 8f, where the phase trajectory deviates severely from the balance point under the autonomous steering control.However, the phase trajectory under the compensation control is around the balance point, which also explains that the stability compensation makes a favorable effect on the lateral stability and path following control.
sation control while performing path following.Figure 8b shows the torque distributi coefficients in stability compensation control.The driving torque of four IWMs with t stability compensation is indicated in Figure 8c.The torque between left IWMs and rig IWMs may be different to ensure stability when the DDAUV plans to change lane at hi velocity, and even some IWMs generate reverse torque to reach the target additional ya moment and improve wheels stability margins.In Figure 8d,e, the comparison resu about the yaw rate and side-slip angle show that the stability of the DDAUV is improv widely under stability compensation control.The phase trajectory portrait for the thr control strategies is shown in Figure 8f, where the phase trajectory deviates severely fro the balance point under the autonomous steering control.However, the phase trajecto under the compensation control is around the balance point, which also explains that t stability compensation makes a favorable effect on the lateral stability and path followi control.

DLC Maneuver on Joint Road
A joint road condition is utilized to further verify the maneuverability in the compensation control, wherein the joint road can be shown in Figure 9.When the DDAEV enters a low adhesion road from a high adhesion road, the ground cannot provide enough steering force for the steering wheels, which will lead to understeer.The proposed coordinated algorithm can help steering through differential tire forces as the following experiment results on joint road at the velocity 60 km/h under DLC maneuver, in Figure 10.
slip angle increase rapidly when entering a low adhesion road.Stability control is unfavorable for the lateral stability instead in this scene, but the proposed compensation control method can effectively coordinate stability and maneuverability, which can also be confirmed through phase trajectory portrait Figure 10f."Stability control" has a good performance in trajectory tracking, but compared with "Autonomous steering", the side deflection angle of the center of mass is increased by about 18.1%.Therefore, the "Compensation control" method is more suitable for the coordinated control in this paper.

Conclusions
The integrated control of path following and lateral stability control is proposed to coordinate lateral stability and path following for the DDAUV.In the upper controller, the path following controller based on APT LQR method is designed, in which the position prediction module is designed and the prediction time is dynamically adjusted according to the tire dynamic characteristics.Based on the reference path, vehicle position, In Figure 10a, since stability control and compensation control both have auxiliary action, the steering behavior is significantly compensated.Figure 10b shows the torque distribution coefficients under the stability compensation control.Figure 10c indicates the torque fluctuation of four wheels under proposed coordinated algorithm.When the DDAUV passes through the road with an adhesion coefficient from 0.6 to 0.1 at about 1.2 s, the right IWMs are significantly strengthened to generate auxiliary steering yaw moment that can help avoid obstacles.Figure 10d,e show the yaw rate and side-slip angle.The yaw rate stops increasing when entering the joint part.That is because the low adhesion road fails to provide enough steering force.Although the stability control is more effective in a single adhesion road, the excessive additional yaw moment makes the side-slip angle increase rapidly when entering a low adhesion road.Stability control is unfavorable for the lateral stability instead in this scene, but the proposed compensation control method can effectively coordinate stability and maneuverability, which can also be confirmed through phase trajectory portrait Figure 10f."Stability control" has a good performance in trajectory tracking, but compared with "Autonomous steering", the side deflection angle of the center of mass is increased by about 18.1%.Therefore, the "Compensation control" method is more suitable for the coordinated control in this paper.

Conclusions
The integrated control of path and lateral stability control is proposed to coordinate lateral stability and path following for the DDAUV.In the upper controller, the path following controller based on APT LQR method is designed, in which the position prediction module is designed and the prediction time is dynamically adjusted according to the tire dynamic characteristics.Based on the reference path, vehicle position, and status information, the APT LQR path following controller can calculate the appropriate frontwheel angle.The improved SMC is applied to seek the desired additional yaw moment.The SMS can be dynamically adjusted according to the lateral stability determined by the phase plane method to maintain a balance between lateral stability and path following.In the lower controller, considering the minimum slip and working load of four tires, a torque distribution algorithm designed to allocate driving torque for producing the desired additional yaw moment.The results of HIL experiments show that the stability compensation controller plays a key role in coordination between the lateral stability and path following.The results show that the trajectory tracking and lateral stability of the hierarchical control architecture with the adaptive prediction module are well coordinated in different working conditions.The limitation of this paper is that it does not consider the influence of driving torque on the analysis of tire slip loss energy and the influence of vehicle speed change on vehicle trajectory tracking and lateral stability.Future research will be carried out from the following aspects.
(1) Considering the great influence of driving torque on the loss energy from tire slip and the working efficiency of IWMs, future research will focus on the integrated control of stability and energy saving through torque vectoring control; (2) In this paper, the fixed speed condition is used in the HIL test, but the speed of the unmanned vehicle changes in the trajectory planned by the decision planning layer, and the speed has a certain influence on the vehicle trajectory tracking and lateral stability.Therefore, the influence of time-varying vehicle speed on its control performance should be considered when designing control strategies in the future.

Figure 1 .
Figure 1.Algorithm framework of integrated path following and lateral stability control.

Figure 1 .
Figure 1.Algorithm framework of integrated path following and lateral stability control.

Figure 1 .
Figure 1.Algorithm framework of integrated path following and lateral stability control.

) 3 .
Design of Autonomous Steering Controller Based on the APT LQR 3.1.Design of LQR Path Following Controller
Path following of the vehicle (b) Torque distribution coefficients (c) Torque distribution (d) Yaw rate (e) Side-slip angle (f) Phase trajectory

Figure 8 .
Figure 8. Experiment results under DLC maneuver on a low adhesion road.Figure 8. Experiment results under DLC maneuver on a low adhesion road.

Figure 8 .
Figure 8. Experiment results under DLC maneuver on a low adhesion road.Figure 8. Experiment results under DLC maneuver on a low adhesion road.

Figure 9 .
Figure 9. Adhesion coefficient of a designed joint road.

Figure 10 .
Figure 10.Experiment results under DLC maneuver on a joint road.

Figure 10 .
Figure 10.Experiment results under DLC maneuver on a joint road.

Table 1 .
Parameters of the DDAUV.

Table 1 .
Parameters of the DDAUV.