Next Article in Journal
Predefined-Time Adaptive Fast Terminal Sliding Mode Control of Aerial Manipulation Based on a Nonlinear Disturbance Observer
Next Article in Special Issue
An Efficient Motion Adjustment Method for a Dual-Arm Transfer Robot Based on a Two-Level Neural Network and a Greedy Algorithm
Previous Article in Journal
Incremental Placement Technology Based on Front-End Design
Previous Article in Special Issue
Prediction of Operational Noise Uncertainty in Automotive Micro-Motors Based on Multi-Branch Channel–Spatial Adaptive Weighting Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wheel Drive Driverless Vehicle Handling and Stability Control Based on Multi-Directional Motion Coupling

School of Automobile and Transportation, Chengdu Technological University, Chengdu 611730, China
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(14), 2744; https://doi.org/10.3390/electronics13142744
Submission received: 31 May 2024 / Revised: 10 July 2024 / Accepted: 10 July 2024 / Published: 12 July 2024
(This article belongs to the Special Issue Applications of Artificial Intelligence in Mechanical Engineering)

Abstract

:
To fully unleash the performance potential of the Wheel Drive Driverless Vehicle (WDDV) and enhance its handling stability across a wide range of extreme operating conditions, this paper proposes a novel approach for designing a multi-directional motion coupling control system. Firstly, an analysis of the unmanned driving modes of the WDDV is conducted, followed by the establishment of a method for defining the control target parameter set for handling stability. Subsequently, a coupled dynamic model that considers the wheel drive counter force is developed. Building this model, a method for estimating the handling stability state is introduced, focusing on improving both handling and stability aspects. Furthermore, by combining the sliding mode control algorithm with the coupled dynamic model, a design methodology for a multi-directional motion coupling control law that adapts to extreme operating conditions is proposed. Finally, through comprehensive simulation experiments and testbed, the effectiveness of the proposed multi-directional motion coupling control system is validated, demonstrating superior handling stability compared to the decoupled control system.

1. Introduction

Improving the transport efficiency of vehicles has been an enduring pursuit, necessitating the rational control and optimal utilization of their dynamic performance advantages. The wheel drive vehicle has multiple independent and controllable power units, which can realize the optimal distribution of longitudinal force, lateral force and vertical force of each wheel, so as to improve the handling and stability of the vehicle. The wheel drive system has introduced a novel implementation approach to enhance the dynamic performance of electric vehicles, positioning it as one of the most promising electric vehicle configurations. The wheel drive system adapts to the development needs of electrification and intelligence of automobiles. It has broad application prospects and has been widely studied [1]. Vehicles equipped with wheel drive systems inherit the inherent dynamic performance advantages, thus assuming a pioneering role in bolstering transport efficiency. Consequently, they inevitably confront the imperative of maneuvering stability control under extreme operating conditions. The control of wheel drive vehicles in extreme conditions presents several challenges: the forceful nonlinearity exhibited by the vehicle’s dynamic response necessitates a broader scope of state parameters for modeling and analysis; the control system of wheel drive vehicles under extreme conditions must strike a delicate balance between complexity and real-time responsiveness. Intelligent technologies imbue vehicles with an array of perceptual sensing units, advanced planning and decision-making systems, and high-performance computing units [2]. Leveraging such cutting-edge hardware platforms enables the utilization of intelligent algorithms [3,4,5] to enhance vehicle performance. Therefore, the fusion of wheel drive vehicles and intelligent technologies in the automotive domain holds considerable potential for resolving the issues pertaining to maneuvering stability control under extreme operating conditions. This integration forms the primary subject of this study: the Wheel Drive Driverless Vehicle (WDDV). The WDDV combines the advantages of wheel drive and driverless vehicle software and hardware, amplifying the complexity of designing its handling and stability control system. To improve the handling and stability of WDDV, addressing the challenges posed by control actuators and controllable degrees of freedom is of paramount importance.
The determination of vehicle handling and stability control objectives usually starts with improving the yaw stability and roll stability of the vehicle. Key state variables for characterizing lateral stability in automobiles include yaw rate and lateral displacement of the center of gravity, which are often directly targeted in the pursuit of lateral stability control [6,7,8,9,10,11,12,13,14]. Additionally, the distribution of tire forces, which is closely related to tire slip angles, plays a vital role in determining the vehicle’s handling stability. Hence, employing the slip angles of the front and rear tires as control targets can also effectively achieve manipulation stability control [15,16]. Restricting the vehicle speed to remain below the maximum permissible speed under given operating conditions serves as an effective means of ensuring driving stability.
Consequently, vehicle roll angle is commonly regarded as a control target for roll stability control [17,18]. In the context of autonomous driving, the maximum achievable speed is also frequently employed as a control target for manipulation stability [19,20]. Furthermore, there is a growing trend to consider multiple parameters as comprehensive control targets [21,22]. Therefore, determining reference values for these control target parameters assumes utmost importance. Common methodologies include the estimation of reference values for control targets under extreme conditions using empirical rules [23,24], as well as techniques such as phase-plane analysis [25,26]. The estimation of reference values for control targets under extreme conditions using empirical parameters or experimental data presents challenges and suffers from issues of accuracy.
Moreover, phase-plane analysis methods [27,28] heavily rely on specific vehicle responses, thereby limiting their applicability across different vehicle models. An ideal approach would involve deriving reference values for control targets under extreme conditions through vehicle dynamics models. However, due to the strong nonlinearity exhibited by vehicle dynamics responses under extreme conditions, modeling becomes challenging, and control systems designed based on complex dynamics models may struggle to meet real-time requirements.
Cooperative control methods for enhancing vehicle handling stability can be categorized into rule-based methods and coupling-based methods. Rule-based methods trigger control strategies based on specific driving conditions and driver input behavior [6,29]. These strategies can be implemented across different electronic control units (ECUs) to reduce computational demands on individual units. However, this approach may introduce oscillations in vehicle states at control activation boundaries, which can undermine control system stability. Moreover, achieving optimal comprehensive vehicle dynamics performance in extreme conditions remains a challenge. Consequently, the increasing computational capabilities of autonomous driving vehicles have prompted a growing emphasis on coupling control methods for multi-actuator systems.
Coupling control strategies, such as four-wheel coupling control [30], four-wheel active steering coupling control [31], coupling control between wheel drive and active steering [32,33], coupling control involving wheel drive, active braking, and active steering [34,35,36,37,38,39,40,41,42], and coupling control incorporating wheel drive, active braking, active steering, and active suspension [28], have been proposed to improve different aspects of vehicle stability and handling. However, addressing the coupling control problem of WDDV space multi-directional motion requires the development of multi-degree-of-freedom coupling dynamic models specific to vehicle types, enabling direct design of multi-directional control laws. In practical design, model simplification is often necessary to achieve complete solutions for control variables, which may hinder the full representation of nonlinear effects in extreme conditions. Additionally, wheel drive vehicles operating at the limit often utilize independent wheel torque control to enhance stability, necessitating the consideration of the coupling effect between longitudinal and roll motions in the dynamic modeling of wheel drive vehicles under extreme conditions.
Summarizing the research status, the deficiencies in the existing research are mainly reflected in the following two aspects:
  • Existing research on manipulation stability control targets exhibits certain limitations: the differences between autonomous vehicles and conventional ones are not adequately addressed; a unified and systematic approach to establishing WDDV manipulation stability control targets is lacking.
  • At present, the WDDV handling and stability control method is mainly based on decoupling control, which is difficult to fully consider the influence of the multi-degree-of-freedom coupling relationship of the vehicle, and the influence of the hub drive reaction force is not considered in the design process of the control system.
In summary, the existing research has not fully addressed the control problem of WDDV for stable maneuvering. Achieving stable maneuvering with coupled control for WDDV remains a challenging task that necessitates further investigation. This paper proposes methodologies for establishing precise control objectives and designing coupled control laws tailored specifically to WDDV. The primary contributions of this study are as follows:
  • Beginning with the classification of automated driving levels, incorporating the unmanned driving mode of WDDV, a comprehensive parameter set for control objectives is established. The approach serves as a valuable reference for quantifying the evaluation indicators of maneuvering stability in unmanned driving scenarios.
  • Based on the established WDDV dynamic model suitable for conventional and extreme conditions, the control objectives of handling and stability is derived. These objectives strike a balance between stability and transportation efficiency, making them adaptable to different operating conditions.
  • A sophisticated coupled control law for WDDV’s multi-directional motion under extreme operating conditions is formulated. This approach is highly effective in meeting the control requirements for stable maneuvering under extreme operating conditions.
The subsequent organization of this paper is as follows: Section 2 provides a concise introduction to the underlying theory. In Section 2.1, we establish a parameter set for control objectives. Section 2.2 focuses on the development of an 8-degree-of-freedom (8-DOF) coupled dynamic model for WDDV. Section 2.3 provides a methodology for estimating the reference values of control objectives. Section 2.4 presents the design of coupled control laws for multi-directional motion in WDDV. Section 3 presents the simulation analysis and experimental verification. Finally, Section 4 offers conclusions.

2. System Framework and Methodology

2.1. Research on Control Target Parameter of WDDV

To establish the target parameter set for handling and stability control, it is essential to define the evaluation indices for handling and stability. Currently, there exists industry-level, national-level, and international-level standards, such as “GB/T 6323-2014 [43], The Handling and Stability Test Method for Vehicle,” “QC/T 480-1999 [44], The Handling and Stability Index Limit and Evaluation Method for Vehicle,” and “ISO 3888-2: 2002 Double-shift International Standard Experiment” [45]. However, given the distinct operational characteristics of driverless vehicles compared to traditional vehicles, the conventional handling and stability evaluation methods are insufficient and require adaptation. Consequently, there is a need to update these methods to address the specific requirements of driverless vehicles.

2.1.1. Establish the Evaluation Paradigm of the Handling Stability

The concept of vehicle handling stability encompasses two fundamental aspects: ma handling and stability. Driverless vehicles exhibit diverse handling patterns, including manned driving, human-machine co-driving, and fully autonomous operation. Moreover, the specific handling and stability requirements vary depending on the intended purpose of different types of driverless vehicles. It is crucial to propose evaluation indices for handling and stability that consider both passengers and vehicles, as the conventional evaluation criteria used in traditional driver handling research are no longer applicable to driverless vehicles. When examining stability, the analysis takes into account both vehicle stability and passenger safety, making the proposed stability evaluation criteria still relevant. However, in the driverless context, the stability index for the vehicle as an electromechanical system becomes less restrictive compared to that of traditional vehicles. For instance, even if the vehicle experiences drifting, fishtailing, or a tire slip on one side, it may still maintain operability.
Nevertheless, these extreme phenomena can significantly impact the passengers’ sense of security. This study focuses on the WDDV system, which represents an L4-level automated passenger vehicle. To determine the target for handling and stability control in WDDV, it is necessary to consider its general application as a passenger vehicle and establish appropriate evaluation criteria that align with the characteristics of its driverless control mode. The absence of a human driver characterizes the driverless maneuvering pattern of the WDDV. Instead, the WDDV integrates a sensory sensing system, decision-making planning system, and intelligent execution system, which collectively emulate the visual perception, sensory capabilities, cognitive processes, and motor functions of a human driver. This integrated system functions as the virtual driver of the WDDV. For virtual drivers, the key performance criterion lies in their ability to seamlessly manipulate and control the vehicle in accordance with the intended driving objectives, known as motion controllability in autonomous driving. In other words, within the established hardware system’s capacity, the WDDV can precisely execute driving actions based on the motion objectives determined by the decision-making planning system, thereby exemplifying its objective motion-tracking ability. Motion tracking in the context of driverless vehicles typically encompasses speed tracking and path tracking, allowing for a clear definition of the handling evaluation indices for the WDDV, specifically in terms of longitudinal speed error and lateral error.
As a passenger vehicle, the evaluation of stability requires consideration of both occupant safety and vehicle endurance. Therefore, the stability evaluation criteria for the WDDV are determined by referring to research findings in traditional vehicle stability assessment. In accordance with the guidelines specified in “QC/T 480-1999, The Handling and Stability Index Limits and Evaluation Methods for Vehicles,” the scoring function for assessing the handling and stability of passenger vehicles is as follows:
N z = 1 / 6 N S + N J + N M + N H + N Q + N W
where N S is the comprehensive score of the serpentine driving test. Through the time domain response analysis of yaw rate and steering wheel angle, the scores formed by the two parameters are calculated, respectively, and then the two scores are weighted to obtain the comprehensive score; N J is the comprehensive score of the steering wheel angle step input test, which is calculated by the time domain response analysis of the yaw rate; N M is the comprehensive score of the steering wheel angle pulse test, which is calculated by the time domain and frequency domain response analysis of the yaw rate; N H is the comprehensive score of the steering returnability test, which is obtained by time domain analysis and calculation of yaw rate; N Q is the comprehensive score of the steering portability test, which is obtained by time domain analysis and calculation of steering wheel steering force; and N W is the comprehensive score of the steady-state rotation test, which is obtained by the time domain analysis of roll angle and yaw rate. The comprehensive scores in Equation (1) are all in the percentage system.
With the advancement of vehicle technology, there is a continuous improvement in vehicle handling and stability, resulting in higher load-carrying efficiency. In 2014, the government introduced the latest test methodology for assessing vehicle handling and stability, documented as “GB/T 6323-2014, Test Method of Vehicle Handling and Stability”. This updated test method includes a handling and stability assessment centered around the steering wheel. The test primarily focuses on high-speed conditions, where the response speed, sensitivity, and accuracy of steering are tested when the vehicle is running at elevated speeds. Under high-speed driving, small disturbances can also easily make the vehicle dynamic response enter a nonlinear state. The higher the speed, the easier it is for the car to enter the extreme driving condition. The test is closely aligned with the investigation of handling and stability control methods under extreme vehicle conditions in this study. Consequently, in accordance with the prescribed test methodology, the scoring Function (1) has been revised as follows:
N z c = 1 / 7 N S + N J + N M + N H + N Q + N W + N Z X
where   N Z X is the comprehensive score of the handling and stability test in the central area of the steering wheel. From the handling and stability test method of the steering wheel center area, the handling and stability of the vehicle is evaluated by the analysis results of the steering wheel torque and yaw rate. N Z X   is also a percentage system.
To facilitate the observation of the relationship between evaluation results of the handling and stability and related evaluation parameters, according to Equation (2) and the influence parameters of each scoring item in the equation, the handling and stability evaluation paradigm of traditional vehicles is as follows:
y t = A t f S γ , δ s + B t f J γ , β + C t f M γ + D t f H γ + E t f Q F Q + F t f W , γ + G t f Z γ , T Q
where y t is the comprehensive evaluation result of passenger vehicle handling and stability; γ is yaw rate; β is sideslip angle; δ s is steering wheel angle; F Q is the steering force of the steering wheel; V is the speed; and T Q is the steering wheel torque; f S γ , δ s is the evaluation function of the serpentine driving test, which describes the calculation method of this result by yaw rate and steering wheel angle in the national standard. f J γ , β is the evaluation function of the steering wheel angle step input test, which describes the calculation method of the result by yaw rate in the national standard; f M γ   is the evaluation function of the steering wheel angle pulse test, which describes the calculation method of the result by yaw rate in the national standard; f H γ is the evaluation function of the steering returnability test, which describes the calculation method of the result by yaw rate in the national standard; and f Q F Q is the evaluation function of the steering portability test, which describes the calculation method of this result by steering wheel steering force in national standards. f W , γ is the evaluation function of the steady-state rotation test, which describes the calculation method of the result through roll angle and yaw rate in the national standard; f Z γ , T Q   is the evaluation function of the test in the central area of the steering wheel. The calculation method of this result by yaw rate and steering wheel torque in the national standard is described. A t , B t , C t , D t , E t , F t , and G t are the weight coefficients of each evaluation result.
According to the national standards, δ s and T Q in   f S γ , δ s and f Z γ , T Q are the parameters used to evaluate the handling of real drivers and f Q F Q is a function used to evaluate the handling of real drivers. There is no real driver in WDDV; the handling evaluation is realized by longitudinal speed error and lateral error. Therefore, combined with the handling evaluation parameters of WDDV, referring to the evaluation paradigm of the handling and stability Equation (3) for traditional vehicles, the evaluation paradigm of WDDV handling and stability was proposed as follows:
y c = A c f S c γ + B c f J γ , β + C c f M γ + D c f C S e s + E c f L J e P + F c f H γ + G c f W , γ + H c f Z c γ
where y c is the comprehensive evaluation result of the handling and stability; e s is the longitudinal speed error; e P is the lateral error; and f C S e s is the evaluation function of vehicle speed tracking test, which describes the method of calculating the result by longitudinal speed error. f L J e P is the evaluation function of path tracking test, which describes the calculation method of this result by lateral error. f J γ , β , f M γ , f H γ , and f W , γ are the same as the evaluation function of the handling and stability of a traditional vehicle; f S c γ is the evaluation function of the serpentine driving test for driverless patterns. The calculation method of the result calculated by the yaw rate is described, which is obtained by referring to the calculation method of f S γ , δ s , without considering the influence of the evaluation parameters for the real driver of the traditional vehicle. f Z c γ is the evaluation function of the test for the steering wheel center area of the driverless pattern, which describes the calculation method of this result by yaw rate. A c , B c , C c , D c , E c , F c , G c , and H c are the weight coefficients of each evaluation result. f C S e s and f L J e P propose corresponding test conditions and evaluation result calculation methods for the driverless vehicle. To adequately evaluate the handling of WDDV and better describe the controllability of the vehicle in the process of vehicle speed tracking and path tracking, f C S e s and f L J e P should be able to make a comprehensive evaluation of the rapidity and accuracy of the control response. In addition, in the process of detailed design, the functions f C S e s and f L J e P should pay attention to the characteristics of virtual drivers and correctly configure the weight coefficients.
In this paper, the target parameter set for the handling and stability control of the WDDV is established by considering the influential factors derived from the evaluation results of handling and stability. Consequently, the investigation does not extensively delve into the evaluation result functions and weight coefficients presented in Equation (4).

2.1.2. Establishment of the Target Parameters Set of Control

The evaluation paradigm for handling and stability, as represented by Equation (4), provides insight into the influence parameters affecting the evaluation results of the WDDV. These parameters encompass lateral error, longitudinal speed error, yaw rate, sideslip angle, and body roll angle. To enhance the handling and stability performance of the WDDV, this research focuses on investigating control methodologies applicable to extreme conditions. Specifically, it explores methods to improve handling and stability by regulating the response of the aforementioned five parameters. Considering the influence parameters derived from the handling and stability evaluation of the WDDV, real-time control processes are affected by factors such as the target path, target speed, target yaw rate, target sideslip angle, target body roll angle, and the corresponding control strategies associated with each target parameter. The path tracking strategy corresponds to the target path, the speed tracking strategy corresponds to the target vehicle speed, the yaw stability control strategy corresponds to the target yaw rate and target sideslip angle, and the roll stability control strategy corresponds to the target body roll angle. Collectively, the target path, target speed, target yaw rate, target sideslip angle, and target body roll angle constitute the control target set for enhancing the handling and stability of the WDDV, as illustrated in Equation (5):
C P = P d , V x d , γ d , β d , d
where the target path P d and the target vehicle’s speed V x d are handling control targets, and the target yaw rate γ d , the target sideslip angle β d , and the target body roll angle d are stability control targets. In the stability control targets, the target yaw rate and the target sideslip angle are the yaw stability control targets, and the target body roll angle is the roll stability control target. In this paper, the principle of cooperative control of handling and stability of WDDV is determined as follows: The primary objective is to enhance handling performance during conventional driving conditions for the WDDV, while also ensuring stability when the WDDV is driving in extreme conditions. The path tracking strategy, vehicle speed tracking strategy, yaw stability control strategy, and roll stability control strategy constitute the core control strategy set of WDDV.
Therefore, the handling and stability control system of WDDV is a multi-objective and multi-strategy complex control system.

2.2. Dynamics System Modeling

To address the research requirements for handling and stability control of the WDDV, a comprehensive 8-DOF vehicle dynamics model was developed. This model accurately captures the longitudinal, lateral, yaw, and body roll motions of the vehicle, as well as the variations in four-wheel rotation. The stability of wheel drive vehicle is improved by wheel torque independent control under extreme conditions. When the torque of the hub motor is transmitted to the body, the vertical component force generated can change the roll dynamic response of the vehicle. Therefore, the dynamic model used in the design of WDDV handling and stability control system needs to focus on the influence of the hub drive reaction force on the vehicle motion.

2.2.1. Vehicle Dynamics Modeling

The 8-DOF vehicle dynamics model is shown in Figure 1. The body dynamics equation of WDDV was established as follows:
The longitudinal motion:
m V ˙ X + V Y γ + m u h u γ ϕ ˙ = F X f r + F X f l cos δ F Y f r +   F Y f l sin δ + F X r l + F X r r m g i + f C D A V x 2 21.15
The lateral motion:
m V ˙ Y + V X γ + a m f b m r γ ˙ m u h u ϕ ¨ = F Y f r + F Y f l cos δ + F X f r +   F X f l sin δ + F Y r l + F Y r r
The yaw motion:
a m f b m r V ˙ Y + V X γ + I z γ ˙ I x z ϕ ¨ = a F Y f r + F Y f l cos δ + d f 2 F Y f l F Y f r sin δ b F Y r l + F Y r r   +   d f 2 F X f r F X f l cos δ + a F X f r + F X f l sin δ + d r 2 F X r r F X r l
The roll motion:
I x ϕ ¨ + m u h u V ˙ Y + V X γ I x z γ ˙ = Δ d f F Y f r + F Y f l cos δ + Δ d r F Y r r + F Y r l D ϕ ˙ C m u g h u ϕ + M z s
  J w ω ˙ i j = T d i j T b i j F X i j R
where m is the vehicle mass; V X is the longitudinal speed; V Y is the lateral speed; m u is the sprung mass of the vehicle; ϕ is roll angle; F X i j is the tire’s longitudinal force; F Y i j is the tire lateral force; δ is the average angle of the front wheel; i is the longitudinal slope of the road; f is the tire rolling resistance coefficient; g is the acceleration of gravity; C D is the air resistance coefficient; A is the windward area of the vehicle; m f is the springing mass of the front suspension; m r is the springing mass of the rear suspension; h u is the distance from the center of mass of the sprung mass to the roll axis; I z is the moment of inertia of the vehicle around the Z axis of the vehicle coordinate system; I x is the rotational inertia of the vehicle around the X-axis of the vehicle coordinate system; I x z is the inertial product of the whole vehicle in the XZ plane of the vehicle coordinate system; a and b are the distances from the front axle and the rear axle to the center of mass, respectively; d f and d r are the front and rear wheel tracks, respectively; M z s is the roll moment generated by wheel drive counter force acting on the body; D is the roll-damping coefficient of the vehicle; C is the roll stiffness of the vehicle; Δ d f is the lateral offset of the front tire caused by unit roll angle; Δ d r is the lateral offset of the rear tire caused by unit roll angle; J w is the rotational inertia of the wheel; T d i j is the wheel driving torque; T b i j is the wheel braking torque; and R is the wheel rolling radius (in X i j or X x i j the footmark i represents f or r , which in turn represents the front wheel or the rear wheel, j represents l or r , which in turn represents the left wheel or the right wheel). The tire force in the model is calculated by the Dugoff nonlinear tire model [46]. The 8-DOF vehicle dynamics model included M z s , which can consider the influence of wheel drive configuration on the roll motion of WDDV. The calculation method of M z s will be studied in detail later.

2.2.2. Analysis of Wheel Drive Counter Force

Previous research by scholars has investigated the influence of wheel drive counterforce on the roll dynamic response of the vehicle body and has established a correlation between the suspension structure and the roll moment M z s [47]. Currently, the WDDV is equipped with a multi-link suspension, but there is limited research available on the wheel drive counterforce specifically related to this suspension type. In most studies, the front and rear suspension systems were projected onto the vehicle’s longitudinal plane and the roll moment M z s was calculated through force analysis. Although this calculation method is straightforward, it overlooks the motion coupling relationship between the front and rear suspension systems, as well as the dynamic changes and spatial relationships within the suspension structure. Consequently, significant estimation errors can occur. To address this, this study utilized simulation data fitting to determine the variation pattern of the roll moment M z s induced by the wheel drive counterforce. By considering the generation mechanism of the roll moment M z s , the established dynamic model incorporates the interaction forces between the wheel motor, steering knuckle, suspension rod system, and vehicle body. Furthermore, this dynamic model facilitates the calculation of the roll moment M z s by obtaining the forces at each connection point between the suspension guide rod system and the body. Therefore, a multi-body dynamics model for the WDDV system was developed, and the variation pattern of the roll moment M z s was obtained through simulation.
By analyzing the generation mechanism of roll moment M z s , the main influence parameters of M z s mainly include the driving torque difference between the left wheel and right wheel in wheel motor, steering wheel angle and the height of each quarter suspension. Therefore, the paradigm of M z s is proposed as follows:
M z s = f T s f , T s r , δ s t e , H s f l , H s f r , H s r l , H s r r
where T s f is the torque difference on the left and right sides of the front axle; T s r is the torque difference between the left and right side of the rear axle; δ s t e is the steering wheel angle; H s f l is the suspension height of the first quarter of the left; H s f r is the suspension height of the right front quarter; H s r l is the suspension height of the left rear quarter; and H s r r is the suspension height of the right rear quarter. The peak torque of the wheel motor is 1250 Nm, so the value range of T s f and T s r is [–2500, 2500] Nm. According to the range of steering wheel angle of the prototype vehicle, the range of δ s t e is determined to be [−540, 540]°. According to the Wheel jump range of the front suspension and the rear suspension of the vehicle, the value range of H s f l ,   H s f r ,   H s r l ,   H s r r is determined as [−70, 50] mm. The negative value indicates suspension compression, and the positive value indicates suspension stretching.
Based on the multi-body dynamics model of WDDV, the experimental design was carried out with T s f , T s r , δ s t e , H s f l , H s f r , H s r l , and H s r r in Equation (11) as the influencing factors. Then, the value of M z s under different influencing factors were obtained by simulation and the response surface model of M z s was finally established. Based on the above ideas, firstly, multi-level simulation tests were carried out on the various influencing factors in the Equation (11), and the following conclusions are drawn: the seven influencing factors showed a weak nonlinear relationship with M z s . Therefore, the relationship between each influencing factor and M z s can be described by using fewer levels for each influencing factor. Based on the above analysis, through the value range of each influencing factor, the level of each influencing factor in the test process are determined as follows: T s f and T s r are taken as 8 levels of equal spacing, δ s t e is taken as 7 levels of equal spacing, H s f l , H s f r , H s r l , and H s r r are taken as −70 mm, 0 mm, 50 mm three levels. Finally, the simulation test is carried out according to the comprehensive test method, and the value of M z s under various influencing factors is recorded, and then the response surface model of M z s is obtained by fitting.
Because the dimension of the response surface model is too high to be fully displayed by a figure, the results of M z s response surface under some assumed conditions are shown: (1) Figure 2 shows the response surface model of the body roll moment M z s generated by the wheel motor torque within the allowable range under the uniform linear condition. There was no steering and suspension deformation under this condition, δ s t e = 0 , H s f l = 0 mm, H s f r = 0 mm, H s r l = 0   mm , and H s r r = 0   mm . The figure shows that the response range of M z s is [−3390, 2835] Nm. (2) Figure 3 shows the response surface model of the body roll moment M z s generated within the allowable range of the wheel motor torque under the right sharp turn condition. In this condition, the steering wheel turned right to the limit, δ s t e = −540°, the left suspension is compressed, and the right suspension is stretched, H s f l = 70 mm, H s f r = −50 mm, H s r l = 70 mm, and H s r r = −50 mm. From the figure, the response range of M z s is [−3530, 3050] Nm. The above data shows that the response range of M z s is large, and it is necessary to consider the influence of M z s on the roll stability of WDDV in stability control.
In the process of real-time control, the simulation data is imported into the control strategy for direct application. However, large volumes of data occupy a lot of storage space. The mathematical model of M z s is established by regression analysis. In theory, polynomials fit any mathematical model. Therefore, the regression model of M z s is constructed by polynomials. Firstly, a 7-variable quadratic polynomial and 7-variable cubic polynomial were constructed as the regression models of M z s , respectively. Then, based on the data obtained from the above simulation test, the parameters of the regression model were calculated by using the solving function of the multivariate nonlinear regression model in Matlab. Furthermore, the results of the solution were analyzed, and the items in the polynomial that had minimal influence on the results were discarded. Finally, the mathematical model of M z s is determined as follows:
M z s = ( 7 H s r l 2 14 H s r r 2 + 3 T s f H s f l + 3 T s f H s f r + 3 T s r H s f l + 3 T s r H s f r + 775 T s f 12824 T s r 222 δ s t e 230 H s f l + 336 H s f r 1753 H s r l + 608 H s r r ) / 10000
The regression model with Equation (12) as M z s is analyzed, and the root mean square error (RMSE) ≤ 94 N/m. Compared with the vehicle, the value is very small and within the acceptable range. Furthermore, in the process of regression analysis, the absolute value of residuals is less than 250 N/m, and the absolute value of residuals less than 100 N/m accounts for 86.5% of the total residuals. In summary, the established mathematical model of M z s has high accuracy and meets the research needs.

2.3. Estimation of Reference Values for Control Targets

The handling and stability state of a vehicle encompasses three distinct states: a stable state with weak nonlinearity, an unstable state with strong nonlinearity, and a critical state characterized by a strong dynamic response between these two states. The stable state corresponds to the vehicle operating under normal conditions, while the critical state pertains to extreme conditions, and the unstable state signifies a safety issue with the vehicle. However, this paper focuses solely on the dynamic aspects of the vehicle in the stable state, excluding the investigation of the vehicle’s behavior in the unstable state. Under normal conditions, virtual drivers anticipate a specific handling and stability state that aligns with their intended driving behavior. This state ensures the virtual driver’s driving intentions are easily accommodated, thereby enabling the driverless vehicle to achieve optimal handling.
Similarly, under extreme conditions, virtual drivers also have an expected handling and stability state. In this state, the driverless vehicle operates at high speeds, pushing the dynamic response of the vehicle to its stability limits. This state represents an ideal condition for enhancing vehicle transportation efficiency. For clarity, we will refer to the aforementioned states as the “expected state of handling” and the “extreme stable state”, respectively. This paper aims to investigate the estimation method for determining the expected handling state and stable limit state of WDDV. The estimation results serve as the reference values for controlling the handling and stability of WDDV.

2.3.1. Estimating the Expected State of the Handling

Under conventional conditions, the determination of target speed and target path in WDDV relies on the decision planning module, considering the driver’s intention and the driving environment. To ensure stability, the target values for yaw rate, body roll angle, and sideslip angle must be established. The target path and target speed define the desired dynamic response envisioned by the virtual driver. Consequently, the reference value for the stability control target is derived to effectively characterize the expected dynamic response of the virtual driver in conventional conditions. Ultimately, this process yields the anticipated handling state, which is characterized by the expected sideslip angle, yaw rate, and roll angle.
The existing research on the anticipated handling state lacks full applicability to WDDV, primarily due to the following three limitations:
(1)
The derived expected state model of vehicle handling is based on a 2-DOF vehicle dynamics model that considers lateral motion and yaw motion. However, it fails to incorporate the expected roll angle, thereby resulting in an incomplete estimation of the anticipated state for WDDV.
(2)
The current research commonly employs the average front wheel angle and real-time vehicle speed as inputs to estimate the expected control state. However, this approach inadequately reflects the anticipated state of the virtual driver, which relies on the target vehicle speed and target path as references.
(3)
The prevailing studies often derive the expected handling state using a linear vehicle model, neglecting the nonlinear characteristics of the tires. Consequently, the estimation of the desired control state becomes inaccurate as the influence of tire nonlinearity is disregarded.
Considering the limitations observed in previous studies, this research aims to address them by simplifying the established 8-DOF dynamic model and deducing the expected handling state of WDDV based on the requirements for conventional conditions.
Under conventional conditions, the driverless vehicle operates smoothly, adhering to the target path and target speed while accounting for real-time control system demands. The dynamic model is based on the following assumptions:
(1)
The vehicle traverses a flat road without any input of vertical road roughness, allowing us to neglect the influence of the tire’s vertical force caused by unevenness.
(2)
The suspension system is taken into consideration, focusing solely on the roll rotation of the vehicle’s sprung mass relative to the unsprung mass. It is assumed that the roll axis of the vehicle remains parallel to the ground, leading to negligible changes in the wheel steering angles ( Δ d f 0 , Δ d r 0 ). This assumption disregards the effects of suspension deformation on wheel steering angles and the vertical displacement of the vehicle’s centroid caused by the relative vertical motion between the partially sprung mass and unsprung mass.
(3)
The vehicle’s speed changes gradually without rapid acceleration or deceleration, maintaining a constant value. Consequently, the dynamic differential equation governing longitudinal motion can be omitted.
(4)
The consideration of aerodynamics is omitted from the model.
(5)
In the context of conventional conditions, emergency sharp turns are not encountered, resulting in a small value for the steering angle ( δ ). Consequently, approximations can be made such that sin δ δ and cos δ 1 .
(6)
Under conventional conditions, WDDV maintains stability, eliminating the need to enhance stability through torque vector control of the wheel motor. Hence, the counter force exerted by the wheel drive has minimal impact on the vehicle body’s roll motion, leading to M z s 0 .
Consequently, WDDV can be simplified into a monorail model, considering both the sprung mass and unsprung mass. This simplified vehicle dynamics model accounts for longitudinal, lateral, yaw, and roll motion, supported by front and rear elastic tires interacting with the ground. Finally, by combining Equations (7)–(9) with the simplified Dugoff tire model [48], which is suitable for conventional conditions, the simplified vehicle dynamics equations for WDDV, describing longitudinal motion, lateral motion, yaw motion, and roll motion can be established as follows:
The lateral motion:
m V ˙ Y + V x d γ + a m f b m r γ ˙ m u h u ϕ ¨ = C Y f f σ f β + a γ V x d δ d + C Y r f σ r β b γ V x d
The yaw motion:
a m f b m r V ˙ Y + V x d γ + I z γ ˙ I x z ϕ ¨ = a C Y f f σ f β + a γ V x d δ d b C Y r f σ r β b γ V x d      
The roll motion:
I x ϕ ¨ + m u h u V ˙ Y + V x d γ I x z γ ˙ = D ϕ ˙ C m u g h u ϕ
where δ d is the average rotation angle of the expected front wheel to achieve the target path tracking calculation. C Y f is the equivalent slip angle stiffness of the front axle tire; C Y r is the equivalent slip angle stiffness of the rear axle tire; and f σ f is the coefficient equation in the Dugoff tire model [48]. The simplified vehicle dynamics model takes the target speed V x d and the δ d calculated according to the target path as input, which is more conducive to describing the expected vehicle dynamics response of the virtual driver and has the characteristics of the driverless vehicle.
In conventional conditions, the curvature of the target speed and target path changes slowly. When the corresponding target speed and target path are reached and stabilized during the driving of the WDDV, V ˙ Y 0 , γ ˙ 0 , ϕ ˙ 0 , ϕ ¨ 0 . Further, combined with the simplified WDDV vehicle dynamics Equations (13)–(15), the expected state of the handling of WDDV is obtained as follows:
γ s μ = V x d L 1 + K μ V x d 2   δ d
β s μ = b + V x d 2 a m + b m r a m f C Y f f σ r L L 1 + K μ V x d 2   δ d
ϕ s μ = V x d 2 m u h u L 1 + K μ V x d 2 m u g h u C   δ d
where γ s μ is the yaw rate expected by the virtual driver; β s μ is the sideslip angle expected by the virtual driver; ϕ s μ is the expected roll angle of the virtual driver; and K μ is the stability factor derived from the estimation of the expected state of the handling of WDDV. Compared with the stability factor K [49], K μ considers the effects of unsprung mass, road adhesion characteristics, and tire nonlinearity on vehicle stability. The expression of K μ is similar to K , as follows:
K μ = 1 L 2 m a C Y r f σ r b C Y f f σ f + b m r a m f 1 C Y r f σ r + 1 C Y f f σ f
where K is expressed as: K = 1 / L 2 m a / C Y r f σ r b / C Y f f σ f . The stability factor of the vehicle is affected by many factors, and the expression of K and K μ is not accurate. From the expression form and derivation process of K and K μ , the difference between K μ and K is that K μ replaces the original tire linear stiffness with nonlinear stiffness. Because K μ is derived from the simplified WDDV dynamics model considering the sprung mass and unsprung mass, respectively, and K is derived from the 2-DOF vehicle dynamics model with the sprung mass and unsprung mass as a whole. Therefore, the expression K μ contains the parameter term composed of unsprung mass. In essence, K and K μ are both stability factors. Their connotations and different values have the same influence on the stability of the vehicle. Compared with K , K μ can consider the influence of the unsprung mass and tire nonlinearity on the stability factor. In theory, the stability factor of the vehicle can be calculated more accurately by the expression of K μ .
The expected state of the handling of WDDV is composed of γ s μ , β s μ and ϕ s μ . Considering the influence of road adhesion characteristics and load changes on tire cornering stiffness through the function f σ i , the influence of tire nonlinearity on vehicle handling under different driving conditions can be considered.

2.3.2. Estimating the Stabilization of Extreme State

Under extreme conditions, the WDDV decision-planning module determines the target path based on driving intention and environmental factors, resulting in the remaining control targets reaching their limits. Estimation of the limit values for target yaw rate, target sideslip angle, target body roll angle, and target speed become crucial. The root cause of vehicle instability lies in the insufficiency of tire-generated forces to sustain stable driving. To ensure the stable operation of WDDV under extreme conditions, it is imperative to describe the limit dynamic response state that adheres to the maximum adhesion force provided by the tire and determine the stability control target.
By considering the maximum tire adhesion as a constraint, the calculation model for determining the limit value of the target yaw rate and target sideslip angle can be derived as follows [28]:
γ m a x = μ g 2 a x + g i + f + C D A V x 2 21.15 m 2 V X
β m a x = γ m a x b V X V X a m + b m r a m f C Y f f σ r L
where γ m a x is the limit value of the target yaw rate; β m a x is the limit value of the target sideslip angle; and τ γ is the safety factor, and the value is greater than 0 and less than 1. Considering many practical factors that are not taken into account in the modeling process, the appropriate safety margin is retained. a x is the longitudinal acceleration; L is the wheelbase.
Under extreme conditions, WDDV frequently enhances lateral stability by employing independent control of torque applied to the left and right wheels. During such instances, the wheel drive counter force becomes significant, necessitating careful consideration of the generated body roll torque. The vehicle’s roll instability during driving is primarily observed in the form of rollover or sideslip. Consequently, the dynamic model equations for WDDV’s rollover and sideslip are established as follows:
m u h c a y + m f + m r R a y = m u g d f + d r 4 h u + m g h c α r o + m f + m r g d f + d r 4 + M z s
m u h c a y m a x + m f + m r R a y m a x = m u g d f + d r 4 h u + m g h c α r o F z f x + F z r x d f + d r 2 + M z s
where h c is the height of the vehicle’s center of gravity; α r o is the lateral slope of the road; a y is the lateral acceleration; a y m a x is the maximum lateral acceleration that the vehicle can achieve in stable driving; F z f x is the vertical force of the front wheel unilateral tire; and F z r x is the vertical force of the rear wheel unilateral tire of the vehicle. When a y m a x points to the left side of the driving direction of the vehicle, F z f x = F z f l , F z r x = F z r l . When a y m a x points to the right side of the driving direction of the vehicle, F z f x = F z f r , F z r x = F z r r . The roll model considers the influence of wheel drive configuration on roll stability through M z s .
The maximum roll angle under the critical state of vehicle rollover is as follows:
m a x 1 = m g d f + d r 4 + h c α r o a y R m f + m r + h c m u + M z s m u g h u
With the maximum adhesion of the tire as the constraint, combined with Equation (23), the maximum roll angle of the vehicle in the critical state of sideslip is as follows:
m a x 2 = m u g d f + d r 4 + m g h c α F z f x + F z r x d f + d r 2 μ g 2 a x + g i + f + C D A V x 21.15 m 2 R m f + m r + h c m u + M z s m u g h u
The maximum target roll angle of WDDV is as follows:
m a x = min m a x 1   ,   m a x 2 τ
where τ is the maximum roll angle safety factor of the vehicle, which is greater than 0 and less than 1. To account for numerous practical factors that are not considered during the modeling process, an appropriate safety margin is maintained.
From the above analysis, the γ m a x , β m a x and m a x are the important part of the extreme state of the vehicle’s stable driving. In the estimation process of γ m a x and β m a x , the longitudinal acceleration a x , lateral acceleration a y and roll angle are obtained by the sensor. The road longitudinal slope i and the road lateral slope α r o are obtained by high-precision map. Based on the camera and radar, the road adhesion coefficient μ is estimated. The longitudinal speed V X is obtained by the sensor (integrated navigation system), or by estimation. The extreme state estimation method of yaw stability plays the intelligent advantage of WDDV.
Owing to the dynamic nature of speed limits under random conditions, their real-time identification and subsequent tracking and control pose significant challenges. Leveraging the intelligent capabilities inherent in WDDV, a predictive approach for estimating the limit speed is proposed. For a detailed exposition of this method, the author refers to their prior research findings [28].

2.4. Design of Multi-Directional Motion Coupling Control Law

2.4.1. Analysis of Coupling Control Principle

The dynamic coupling characteristics of the vehicle system manifest through the interplay among the motion relationship, tire forces, and dynamic load distribution. Equations (6)–(9) elucidate the following aspects:
(1)
The established 8-DOF vehicle dynamics model equation for WDDV encompasses longitudinal speed V X , lateral speed V Y , yaw rate γ , and roll angle ϕ . Alterations in motion parameters in any direction induce corresponding changes in the motion response of the other three directions, highlighting the coupling characteristics inherent in the vehicle’s motion relationship.
(2)
The equations governing motion in each direction incorporate the longitudinal force F X i j and lateral force F Y i j exerted by the tires. The roll motion equation accounts for the body roll torque arising from the wheel drive counter force, which is engendered by the longitudinal force exerted by the wheel motor through the tire and onto the ground. The Dugoff tire model is employed to calculate the longitudinal and lateral forces of the tire, with the model relying on the tire’s vertical force as input. By incorporating the Dugoff tire model, the vehicle dynamics model effectively captures the coupling characteristics of tire forces.
(3)
During driving, braking, steering, and other operational scenarios, the redistribution of the vertical load on the four wheels occurs due to changes in the centroid position and the effects of inertia. The coupling relationship of tire forces leads to alterations in the longitudinal and lateral forces exerted by the tires, subsequently affecting the overall vehicle response. By employing the Dugoff tire model, the vehicle dynamics model takes into account the dynamic variations in wheel vertical load and its impact on the dynamic response of WDDV.
The 8-DOF vehicle dynamics model for WDDV has been devised to comprehensively address the dynamic coupling characteristics inherent in the system. Building upon this dynamic model, an all-direction motion control law is formulated for WDDV, serving as the foundation for achieving multi-directional motion coupling control within the WDDV space. Moreover, the sliding mode control algorithm, renowned for its swift response, robustness against parameter variations and disturbances, independence from online system identification, and straightforward physical implementation, has garnered significant attention from industry scholars for actively regulating vehicle dynamics performance [50]. Therefore, this study adopts the sliding mode algorithm to design the Multi-directional Motion Coupling Control System (MMCCS) for WDDV.
The MMCCS of WDDV is shown in the dotted box in Figure 4. The input parameters of MMCCS are road adhesion coefficient μ, road longitudinal slope i , yaw angle ψ R , and other parameters of the WDDV perception fusion module and decision planning module. Other commonly used input parameters for vehicle dynamic performance control include tire force ( F Y i j , F X i j , F z i j ), longitudinal speed V X , lateral speed V Y , yaw rate γ , sideslip angle β , body roll angle ϕ , etc. These parameters include the control target parameters that have been identified in the previous study, including the target path point set ( X n , Y n ), the target vehicle speed V X d , the expected state of the handling: γ s μ , β s μ , ϕ s μ and yaw stability extreme state: γ m a x , β m a x , and roll stability extreme state: m a x . The output parameters of MMCCS include steering wheel Angle δ s t e , the driving torque of each in-wheel motor T d f l , T d f r , T d r l , T d r r , and the braking torque of each wheel T x f l , T x f r , T x r l , T x r r   and the control force of each quarter suspension F x f l , F x f r , F x r l , F x r r . MMCCS provides parameter input for the control system with the environment perception and road reconstruction components, path planning components and speed planning components in the perception fusion module and decision-planning module and gives full play to the intelligent advantages of WDDV.

2.4.2. Design of Coupled Control Law

The multi-directional motion coupling control law of WDDV encompasses the path tracking law, vehicle speed tracking law, yaw stability control law, and roll stability control law. By incorporating the existing research findings of the author, along with Equations (6)–(9) and the sliding mode control algorithm, the design of the path tracking law, speed tracking law, and yaw stability control law can be achieved [28], which will not be elaborated upon. In this paper, the same approach is adopted to formulate the roll stability control law, considering the wheel drive counter force. The design of the sliding mode for roll stability control is as follows:
s A = e A ˙ + c A e A
where c A is an adjustable parameter, and c A > 0; e A is the roll angle tracking error, which is calculated by the current roll angle and the target body roll angle d μ . According to the cooperative principle of improving the handling of WDDV under conventional conditions and ensuring stability under extreme conditions, the reference value of the roll angle control target is composed of s μ and m a x . The specific expression for d μ is as follows:
e A = d μ
d μ = s μ ,           s μ < m a x m a x ,     s μ m a x
Using the saturation function s a t s A replaces the sign function s g n s A [50], and a new exponential-reaching law is obtained as follows:
s A ˙ = ε A · s a t s A k A s A
where ε A is an adjustable parameter; k A is an adjustable parameter.
Combined with Equations (9), (27), (28) and (30), the sliding mode control law for calculating the roll control torque M X is as follows:
M X = I x   ε A   s a t s A   ¨ d μ + c A ˙ ¨ d μ + k A ˙ ¨ d μ + c A d μ D ˙ + I X Z γ ˙ + C g h u m u + a y h u m u + Δ d f F Y f r + F Y f l cos δ + Δ d r F Y r r + F Y r l + M z s
Using Lyapunov’s second stability criterion, the stability analysis of the control system composed of the designed roll stability control law is carried out. The Lyapunov function is defined as V A = s A 2 / 2 . Then V A ˙ ε A s A k A s A 2 can be obtained from Equation (30). ε A > 0, k A > 0, s A 0 , s A 2 0 , so V A ˙ 0 . Therefore, the roll stability control system is asymptotically stable.
During roll control, the objective is to minimize the discrepancy between the roll angles induced by the roll control force in both the front and rear suspensions, aiming to prevent the emergence of supplementary torsional moments on the vehicle’s front and rear bodies. The computation equation for the supplementary suspension-induced driving force is expressed as follows:
F x f l = M X F z f l + F z f r F z d f
F x f r = F x f l
F x r l = M X F z r l + F z r r F z d r
F x r r = F x r l
where F x f l is the control force of the left front suspension; F x f r is the control force of the right front suspension; F x r l is the control force of the left rear suspension; and F x r r is the control force of the right rear suspension.

3. Experiments and Results Analysis

3.1. Simulation Research

To validate the efficacy of the proposed multi-directional motion coupling control methodology, a simulation analysis is conducted employing a representative WDDV model. The example WDDC model’s principal parameters are enumerated in Table 1, at which point the vehicle dynamics model is established in CarSim, taking into consideration the data provided therein. Through thorough validation against prototype vehicle experimental data and comparative analysis with analogous models, the CarSim model demonstrates a commendable level of simulation accuracy. The MMCCS model is implemented in Matlab 2021a/Simulink, while the performance assessment entails a comparative analysis between the MMCCS and a decoupling control system (DCS). The control law for the DCS, detailed in the author’s prior research endeavors [28], adheres to conventional control methodologies.

3.1.1. Setting Conditions

To fulfill the requirements for evaluating the control effectiveness and advantages of the MMCCS, the simulation scenario design encompasses a comprehensive range of driving conditions for the WDDV, including both conventional and extreme scenarios. The following specifications are defined for the solution configuration, encompassing road conditions and driving regulations within the simulation framework:
(1) The simulation adopts a time step of 0.01 s.
(2) The initial velocity of the WDDV is initialized to 0 km/h.
(3) The selection of the target path aims to encapsulate the complex operational conditions encountered during vehicle maneuvering, thereby facilitating an assessment of the vehicle’s handling and stability control efficacy across normal and extreme scenarios. In adherence to industry-standard procedures for evaluating vehicle handling and stability, the target path is designed as a composite trajectory, combining a double lane change condition with a serpentine condition, as visually depicted in Figure 5.
(4) To conduct a comprehensive evaluation of the control system, a diverse set of test roads was designed. Under each road condition, including extreme and conventional speeds, corresponding test analyses were performed. Due to space limitations, this article focuses on the analysis of dry road and ice-snow road scenarios as exemplars. The adhesion coefficient for the dry road is 0.8, while for the ice and snow road, it is 0.15. The extreme speed planning methodology for each respective road condition is as follows: the principle is to ensure that the vehicle does not exceed the extreme speed and that the vehicle’s speed is reduced to a safe level before entering the corner. The conventional speed is significantly lower than the extreme speed. Under dry road conditions, the corresponding speeds relating to the target path are displayed in Figure 6a. In this figure, V x _ m a x f represents the maximum target speed limit, V x d f denotes the extreme speed planned by the WDDV speed planning component, adhering to the principle of not surpassing V x _ m a x f and reducing the speed to a safe level before turning and V x d l represents the significantly lower conventional speed.
WDDV is driven under the target path illustrated in Figure 5 at both extreme speed V x d f and conventional speed V x d l . The lateral acceleration response is presented in Figure 7a, where a y f represents the lateral acceleration response at extreme speed V x d f , and a y l represents the lateral acceleration response at the normal speed V x d l . As observed in Figure 7a, the maximum value of a y f reaches 7.31 m/s2, which closely approximates the lateral acceleration limit of 7.84 m/s2 that the vehicle can achieve on the dry road with an adhesion coefficient of 0.8. This finding substantiates the reasonableness of the chosen extreme speed V x d f , which enables WDDV to operate under extreme conditions. Furthermore, the maximum value of a y l amounts to 4.56 m/s2, significantly lower than 7.84 m/s2. This demonstrates that the set conventional speed V x d l is rational and enables WDDV to operate under conventional conditions. For the ice-snow road condition, the corresponding speeds associated with the target path are depicted in Figure 6b.
Analogous to the ice-snow road conditions, the lateral acceleration responses corresponding to the respective speeds are presented in Figure 7b. Through the same analytical process applied to each vehicle’s speed and its corresponding lateral acceleration under dry road conditions, it can be concluded that the planned extreme speed and set conventional speed are reasonable and consistent. For brevity, further repetition of the analysis process will be omitted.

3.1.2. Verification of Validity

The control system comprises MMCCS and DCS, with the vehicle dynamics model established in CarSim serving as the controlled object. The simulation test was conducted based on the aforementioned settings. In the subsequent simulation results, the No.1 curve represents the effect results of MMCCS, while the No.2 curve represents the effect results of DCS.
Under conventional conditions, WDDV operates at the conventional speed V x d l on both the dry road and the ice-snow road, following the target path depicted in Figure 5 and Figure 6a,b, respectively. Considering the emphasis on handling under conventional conditions, the lateral error and longitudinal speed error outcomes are presented in Figure 8 and Figure 9. The figures indicate that the absolute values of lateral error and longitudinal speed error for both control systems on the dry road do not exceed 0.03 m and 0.007 m/s, respectively. On the ice-snow road, the absolute value of lateral error remains below 0.013 m, while the absolute value of longitudinal speed error does not exceed 0.006 m/s. To facilitate a quantitative evaluation of the control effectiveness, the simulation results were organized, and a summary of the simulation outcomes under conventional conditions is presented in Table 2.
Figure 8 and Figure 9 demonstrate that the MMCCS and DCS effectively achieve path-tracking control and vehicle speed-tracking control for WDDV in conventional scenarios. Consequently, both MMCCS and DCS exhibit favorable handling control performance on WDDV. Furthermore, MMCCS demonstrates superior performance in minimizing lateral error and longitudinal speed error when compared to DCS, particularly in dry and ice-snow road conditions under conventional circumstances. Table 2 illustrates that the reduction in peak lateral error ranges from 17% to 31%, while the reduction in peak longitudinal speed error ranges from 17% to 43%. It is proved that the control expectation state estimation method proposed in this paper can estimate the expected dynamic response of the virtual driver more accurately than the traditional method. The control strategy is used to control the vehicle to track the estimated control desired state in this paper, which is more conducive to improving the path tracking accuracy and vehicle speed tracking accuracy. Thus, MMCCS proves to be more conducive to enhancing the handling of WDDV than DCS in conventional conditions.
The evaluation of the control effect of WDDV under extreme conditions necessitates careful consideration of stability assessment parameters while ensuring that the handling evaluation parameters remain within an acceptable range. Utilizing the extreme speed V d x f depicted in Figure 6a for both dry and ice-snow road conditions, the response outcomes of lateral error, longitudinal speed error, yaw rate, sideslip angle, and roll angle, resulting from WDDV driving along the target path illustrated in Figure 5, are presented in Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.
Figure 10a,b illustrate that under the extreme conditions corresponding to the dry road, the absolute value of lateral error remains within 0.08 m, while the absolute value of longitudinal speed error remains within 0.03 m/s. Similarly, Figure 12b and Figure 13a demonstrate that under the extreme conditions associated with the ice-snow road, the absolute value of lateral error does not exceed 0.1 m, and the absolute value of longitudinal speed error does not surpass 0.06 m/s. Both control systems effectively ensure that WDDV exhibits the required handling under these extreme conditions.
Figure 11, Figure 12a, Figure 13b, and Figure 14 indicate that MMCCS exhibits a narrower range of control over the yaw rate, sideslip angle, and roll angle of WDDV compared to DCS under both road conditions, thereby enhancing the stability of WDDV. Specifically, under extreme ice-snow road conditions, the stability evaluation parameter range of MMCCS is smaller than that of DCS, whereas this distinction is less pronounced under extreme dry and wet road conditions. These findings suggest that the ability of the handling and stability control system to improve WDDV’s stability becomes limited under extremely low adhesion conditions on icy and snowy roads. To quantitatively evaluate the control effect of MMCCS and DCS under extreme conditions, the simulation results of WDDV during double lane changes and serpentine road sections were analyzed, and a summary of the simulation outcomes is presented in Table 3.
Based on the data in Table 3, the following conclusions can be drawn: (1) Under extreme conditions, MMCCS is more effective than DCS in enhancing WDDV’s handling. The improvement rate for reducing the peak lateral error ranges from 17% to 33%, while the improvement rate for reducing the peak longitudinal speed error ranges from 0% to 34%. (2) Under extreme conditions, compared to DCS, MMCCS achieves an improvement rate for reducing the peak yaw rate ranging from 9% to 11% and an improvement rate for reducing the peak sideslip angle ranging from 8% to 14%. This indicates that MMCCS demonstrates superior yaw stability control. Additionally, the improvement rate for reducing the peak roll angle ranges from 10% to 14%, indicating that MMCCS exhibits better roll stability control. It is proved that the stability limit state estimation method proposed in this paper can more accurately estimate the dynamic response of the vehicle to maintain stable driving under extreme conditions than the traditional method. Through the control strategy, the vehicle is controlled to track the stable extreme state estimated in this paper, which is more conducive to reducing yaw rate, sideslip angle and body roll angle. At the same time, the above conclusions also prove that the coupling control method proposed in this paper takes into account the longitudinal motion, lateral motion, yaw motion and roll motion and coupling effect of WDDV, and avoids the control conflict of each direction motion, so it can improve the yaw and roll stability while reducing the lateral path tracking error and longitudinal speed tracking error. These results confirm that MMCCS is more effective in improving the stability of WDDV under extreme conditions. The above conclusions demonstrate that, in comparison to DCS, MMCCS offers more effective control over the dynamic response of WDDV during extreme driving conditions, thereby enhancing both handling and stability.

3.2. Experimental Verification

The assessment of MMCCS’s actual control effectiveness primarily relies on real vehicle testing, which entails significant costs, lengthy research and development cycles, and substantial security risks. However, the testbed test partially addresses the limitations of pure digital simulation, offering benefits such as a secure testing process, reduced testing costs, shorter debugging and optimization cycles, and better control over the testing environment. To further validate the suitability of MMCCS for extreme conditions in practical applications, experimental verification was conducted as an extension of the initial pure digital simulation verification process.

3.2.1. Building a Testbed

The testbed for MMCCS comprises two distinct components: the software system and the hardware system.
The software system of the testbed encompasses two key elements: the central computing unit and the virtual WDDV. Following industry practices, the central computing unit’s software system consists of the Ubuntu operating system and ROS. The MMCCS C++ code, derived from the automatic conversion of Matlab/Simulink, operates within the ROS framework on the central computing unit of the WDDV. The virtual WDDV terminal model primarily encompasses the controlled object’s WDDV model. To better compare with the pure digital simulation results, the model uses the same controlled WDDV model as the pure digital simulation verification, and the model is established in CarSim. Figure 15 depicts the software platform architecture of the testbed.
During the operation of the autonomous vehicle, the control system’s computational power and the communication interface’s capabilities are subjected to stringent requirements. To meet these demands, an industrial control computer is employed as the central computing platform in the implementation of the driverless vehicle. As a result, the industrial control computer functions as the physical controller in the WDDV’s testing. Internally, the industrial control computer’s functionality is realized through the automatic compilation of program code using Matlab/Simulink. In the experimental setup, a personal laptop serves as the virtual WDDV simulation platform, with the joint simulation of Matlab/Simulink and CarSim facilitating the creation of the virtual WDDV. The data communication between the industrial control computer and the laptop is established through the dual-channel Kvaser Memorator Pro CAN analyzer. Figure 16 illustrates the testbed configuration for MMCCS.

3.2.2. Verifying the Adaptability of Extreme Conditions

The verification of MMCCS’s adaptability to extreme conditions in a real controller environment was accomplished by comparing the outcomes of both pure digital simulation tests and testbed tests. The verification process encompassed the same solution settings, comprehensive road conditions, and driving rule designs as described in the aforementioned pure simulation test. In the subsequent analysis, the simulation test results are denoted as X s i m , while the testbed test results are represented as X h i l .
Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21 illustrate the response results of lateral error, longitudinal speed error, yaw rate, sideslip angle, and roll angle exhibited by the WDDV while driving under extreme conditions on both dry and ice-snow roads. Figure 17, Figure 18 and Figure 19a demonstrate the comparison between the testbed test results and the pure digital simulation test results under extreme conditions encountered on dry roads. The following observations can be made: (1) The response trends of each parameter were consistent between the two test methods. (2) Both lateral error and longitudinal speed error exhibited significant increases. The absolute value of lateral error ranged from 0.15 m, while the absolute value of longitudinal speed error ranged from 0.025 m/s. Despite these increases, the magnitudes of both errors remained within acceptable limits. (3) The yaw rate and sideslip angle displayed slight increases, exhibiting a satisfactory correspondence. (4) Notably, the roll angle demonstrated a significant increase, but its absolute value was effectively controlled within the range of 0.062 rad.
Figure 19b, Figure 20 and Figure 21 present a comparative analysis between the testbed test results and the results obtained from the pure digital simulation test under extreme conditions simulated on ice and snow-covered roads. The following observations can be made: (1) Lateral error response trends were consistent, but the values were quite different. (2) The lateral error and the longitudinal speed error increased significantly. The absolute value of lateral error was in the range of 0.11 m, and the absolute value of longitudinal speed error was in the range of 0.0076 m/s. The values of the two errors are still within the acceptable range. (3) The yaw rate and the sideslip angle increased slightly, and the correspondence is good. (4) Noticeable increase in roll angle: the roll angle exhibited a significant increase, but effective control measures limited its absolute value to within the range of 0.012 rad.
The above phenomena show that in the testbed test process, the control ability of MMCCS decreases, and the path tracking error, vehicle speed tracking error and roll angle increase. At this time, the yaw rate and the sideslip angle should also become larger. However, the turning radius of the vehicle is closely related to the yaw rate and the sideslip angle. Because the path tracking error increases too much, it is difficult to control WDDV to drive strictly according to the target path, the turning radius of the vehicle increases, and the yaw rate and the sideslip angle of the center of mass decrease too much. Finally, the phenomenon of the yaw rate and sideslip angle decreasing is shown. Compared to the pure digital simulation test environment, the handling control ability of MMCCS in the testbed test environment experiences a reduction, yet it still meets the fundamental requirements for handling control. The yaw stability control ability remains unchanged, while the roll stability control ability is diminished, but the preservation of roll stability is still ensured.
The testbed test results indicate the following: (1) In the real controller environment, the handling and stability control ability of MMCCS under extreme conditions experiences a reduction due to code conversion and hardware real-time performance. However, it effectively enables the handling and stability control of WDDV. (2) Within the real controller environment, MMCCS achieves coupled control of longitudinal motion, lateral motion, yaw motion, and roll motion of WDDV under extreme conditions. It ensures stability while satisfying the requirements of handling control, demonstrating commendable real-time control capabilities.

4. Conclusions

This paper analyzes the technical characteristics of vehicle handling stability in the context of driverless systems. It establishes control targets for the handling stability of Wirelessly Driven Driverless Vehicles (WDDV). An 8-DOF coupled dynamic model of WDDV is developed based on its configuration characteristics and control requirements. The study investigates the mechanism of handling and stability, deriving the control expectation state and stable extreme state estimation model for WDDV. Furthermore, the paper designs a Multi-directional Motion Coupling Control System (MMCCS) for WDDV based on the 8-DOF coupled dynamics model, which is verified through digital simulation and testbed verification. In the future work, the contents that need to be improved are as follows: (1) Under extreme conditions, it is difficult to accurately describe the real response through the vehicle model. To verify the application feasibility of the proposed method, real vehicle verification under extreme conditions is needed. (2) The test method of the handling and stability evaluation index and the determination method of index limit for the driverless vehicle need to be further studied. (3) Under the background of electrification and intelligence, the controllable degree of freedom of WDDV is greatly increased. With the complexity of the model, the real-time challenge of the control system will also be brought. How to comprehensively consider the driving conditions and carry out the “degradation” and “evolution” of the model in real-time remains to be further studied.

Author Contributions

Conceptualization, Y.L. and K.W.; formal analysis, L.D.; investigation, Y.L.; methodology, K.W. and Z.W.; resources, H.W.; supervision, Z.W. and K.W.; validation, Y.L. and K.W.; visualization, H.W. and L.D.; writing—original draft, Y.L. and K.W.; writing—review and editing, Y.L. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Talent Project of Chengdu Technological University under Grant 2023RC054, and the Science Laboratory Open Fund Project of Chengdu Technological University under Grant 2023LOF020.

Data Availability Statement

Data are available from authors upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

a Distance from center of gravity to front axle.
a x Longitudinal acceleration.
a y Lateral acceleration.
a y m a x Maximum lateral acceleration that the vehicle can achieve in stable driving.
A Windward area of the vehicle.
A t , B t , C t , D t , E t , F t , G t Weight coefficients of the handling and stability evaluation paradigm of traditional vehicles.
A c , B c , C c , D c , E c , F c , G c , H c Weight coefficients of the handling and stability evaluation paradigm of the wheel drive driverless vehicle.
bDistance from center of gravity to rear axle.
C Roll stiffness of the vehicle.
c A Adjustable parameter.
C D Air resistance coefficient.
C Y f Equivalent slip angle stiffness of the front axle tires.
C Y r Equivalent slip angle stiffness of the rear axle tires.
D Roll-damping coefficient of the vehicle.
DCSDecoupling control system.
d f Front wheel track.
d r Rear wheel track.
e A Roll angle tracking error.
e P Lateral error.
e s Longitudinal speed error.
f σ f Coefficient equation in the Dugoff tire model.
f Tire rolling resistance coefficient.
f C S e s Evaluation function of vehicle speed tracking test.
f H γ Evaluation function of the steering returnability test.
f J γ , β Evaluation function of the steering wheel angle step input test.
f L J e P Evaluation function of path tracking test.
f M γ Evaluation function of the steering wheel angle pulse test.
f Q F Q Evaluation function of the steering portability test.
F Q Steering force of the steering wheel.
f S γ , δ s Evaluation function of the serpentine driving test.
f S c γ Evaluation function of the serpentine driving test for driverless patterns.
f W , γ Evaluation function of the steady-state rotation test.
F X Longitudinal control force.
F x f r Control force of the right front suspension.
F x f l Control force of the left front suspension.
F X i j Tire’s longitudinal force.
F x r l Control force of the left rear suspension.
F x r r Control force of the right rear suspension.
F Y i j Tire lateral force.
f Z γ , T Q Evaluation function of the test in the central area of the steering wheel.
f Z c γ Evaluation function of the test for the steering wheel center area of the driverless pattern.
F z f x Vertical force of the front wheel unilateral tire.
F z r x Vertical force of the rear wheel unilateral tire of the vehicle.
g Acceleration of gravity.
h c Height of vehicle center of gravity.
H s f l Suspension height of the first quarter of the left.
H s f r Suspension height of the right front quarter.
H s r l Suspension height of the left rear quarter.
H s r r Suspension height of the right rear quarter.
h u Distance from the center of mass of the sprung mass to the roll axis.
i Longitudinal slope of the road.
I x Roll moment of inertia of vehicle.
I x z Inertial product of the whole vehicle in the XZ plane of the vehicle coordinate system.
I y Pitch moment of inertia of vehicle.
I z Yaw moment of inertia of vehicle.
J w Rotational inertia of the wheel.
K Stability factor.
K μ Stability factor derived from the estimation of the expected state of the handling of WDDV.
k A Adjustable parameter.
LWheelbase.
m Vehicle mass.
m f Springing mass of the front suspension.
m r Springing mass of the rear suspension.
m u Sprung mass of the vehicle.
MMCCSMulti-directional Motion Coupling Control System.
M X Roll control torque.
M z Yaw control torque.
M z s Roll moment generated by wheel drive counter force acting on the body.
N H Comprehensive score of the steering returnability test.
N J Comprehensive score of the steering wheel angle step input test.
N M Comprehensive score of the steering wheel angle pulse test.
N S Comprehensive score of the serpentine driving test.
N W Comprehensive score of the steady-state rotation test.
N Q Comprehensive score of the steering portability test.
N Z X Comprehensive score of the handling and stability test in the central area of the steering wheel.
P d Target path.
R Wheel rolling radius.
T b i j Wheel braking torque.
T d f l Driving torque of the front left in-wheel motor.
T d f r Driving torque of the front right in-wheel motor.
T d i j Wheel driving torque.
T d r l Driving torque of the rear left in-wheel motor.
T d r r Driving torque of the rear right in-wheel motor.
T Q Steering wheel torque.
T s f Torque difference on the left and right sides of the front axle.
T s r Torque difference between the left and right side of the rear axle.
T x f l Braking torque of the front left wheel.
T x f r Braking torque of the front right wheel.
T x r l Braking torque of the rear left wheel.
T x r r Braking torque of the rear right wheel.
V Vehicle speed.
V X Longitudinal speed.
V x d Target vehicle’s speed.
V Y Lateral speed.
WDDVWheel Drive Driverless Vehicle.
X h i l Testbed test results X denotes lateral error, longitudinal speed error, yaw rate, sideslip angle, and roll angle.
( X n , Y n )Target path point set.
X s i m Simulation test results X denotes lateral error, longitudinal speed error, yaw rate, sideslip angle, and roll angle.
y c Comprehensive evaluation result of the handling and stability.
y t Comprehensive evaluation result of passenger vehicle handling and stability.
α r o Lateral slope of the road.
β m a x Limit value of the target sideslip angle.
β Sideslip angle.
β s μ Sideslip angle expected by the virtual driver.
β d ,Target sideslip angle.
γ m a x Limit value of the target yaw rate.
γ d Target yaw rate.
γ Yaw rate.
γ s μ Yaw rate expected by the virtual driver.
δ Average angle of the front wheel.
δ d Average rotation angle of the expected front wheel to achieve the target path tracking calculation.
Δ d f Lateral offset of the front tire caused by unit roll angle.
Δ d r Lateral offset of the rear tire caused by unit roll angle.
δ s Steering wheel angle.
δ s t e Steering wheel angle.
ε A Adjustable parameter.
μRoad adhesion coefficient.
τ Maximum roll angle safety factor of the vehicle.
τ γ Safety factor.
ϕ s μ Expected roll angle of the virtual driver.
ϕ Roll angle.
d Target body roll angle.
ψ R Yaw angle.

References

  1. Mosin, M.; Popov, N.; Anibroev, V.; Vilberger, M.; Domakhin, E. Engine power distribution system for four-wheel drive autonomous electric vehicle. Energy Rep. 2023, 9, 115–122. [Google Scholar] [CrossRef]
  2. Keqiang, L.I.; Yifan, D.A.I.; Shengbo, L.I.; Mingyuan, B.I.A.N. State-of-the-art and technical trends of intelligent and connected vehicles. Automot. Saf. Energy Saving 2017, 8, 1–14. [Google Scholar]
  3. Huang, H.; Huang, X.; Ding, W.; Yang, M.; Yu, X.; Pang, J. Vehicle vibro-acoustical comfort optimization using a multi-objective interval analysis method. Expert Syst. 2023, 213, 119001. [Google Scholar] [CrossRef]
  4. Huang, H.; Lim, T.C.; Wu, J.; Ding, W.; Pang, J. Multitarget prediction and optimization of pure electric vehicle tire/road airborne noise sound quality based on a knowledge-and data-driven method. Mech. Syst. Signal Process. 2023, 197, 110361. [Google Scholar] [CrossRef]
  5. Huang, H.; Huang, X.; Ding, W.; Yang, M.; Fan, D.; Pang, J. Uncertainty optimization of pure electric vehicle interior tire/road noise comfort based on data-driven. Mech. Syst. Signal Process. 2022, 165, 108300. [Google Scholar] [CrossRef]
  6. Liang, J.; Lu, Y.; Yin, G.; Fang, Z.; Zhuang, W.; Ren, Y.; Li, Y. A Distributed Integrated Control Architecture of AFS and DYC Based on MAS for Distributed Drive Electric Vehicles. IEEE Trans. Veh. Technol. 2021, 70, 5565–5577. [Google Scholar] [CrossRef]
  7. Liu, J.; Song, J.; Li, H.; Huang, H. Direct yaw-moment control of vehicles based on phase plane analysis. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 2459–2474. [Google Scholar] [CrossRef]
  8. Fan, L.; Wang, J.; Deng, M.; Peng, Y.; Bao, X.; Wei, H. Torque coordinated control of the through-the-road (TTR) 4-wheel-drive (4WD) hybrid vehicle under extreme road conditions. Sci. Rep. 2023, 13, 11564. [Google Scholar] [CrossRef] [PubMed]
  9. Zhai, L.; Wang, C.; Hou, Y.; Liu, C. MPC-based integrated control of trajectory tracking and handling stability for intelligent driving vehicle driven by four hub motor. IEEE Trans. Veh. Technol. 2022, 71, 2668–2680. [Google Scholar] [CrossRef]
  10. Zhang, X.; Liang, H.; Wang, X.; Li, Q. Integrated Direct Yaw Control and Antislip Regulation Mixed Control of Distributed Drive Electric Vehicle Using Cosimulation Methodology. Math. Probl. Eng. 2022, 2022, 6749649. [Google Scholar] [CrossRef]
  11. Zhang, F.; Xiao, H.; Zhang, Y.; Gong, G. Distributed Drive Electric Bus Handling Stability Control Based on Lyapunov Theory and Sliding Mode Control. Actuators 2022, 11, 85. [Google Scholar] [CrossRef]
  12. Shen, Y.; Zhao, Y.; Deng, H.; Lin, F.; Shen, H. Coordinated control of stability and economy of distributed drive electric vehicle based on Lyapunov adaptive theory. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2023, 238, 1535–1549. [Google Scholar] [CrossRef]
  13. Guo, N.; Zhang, X.; Zou, Y.; Lenzo, B.; Du, G.; Zhang, T. A supervisory control strategy of distributed drive electric vehicles for coordinating handling, lateral stability, and energy efficiency. IEEE Trans. Transp. Electrif. 2021, 7, 2488–2504. [Google Scholar] [CrossRef]
  14. Huang, C.; Wu, X.; Wu, C.; Wang, J.; Shu, X. Rule-based torque vectoring distribution strategy combined with slip ratio control to improve the handling stability of distributed drive electric vehicles. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2024, 09544070231215926. [Google Scholar] [CrossRef]
  15. Liang, J.; Feng, J.; Lu, Y.; Yin, G.; Zhuang, W.; Mao, X. A Direct Yaw Moment Control Framework Through Robust T-S Fuzzy Approach Considering Vehicle Stability Margin. IEEE/ASME Trans. Mechatron. 2024, 29, 166–178. [Google Scholar] [CrossRef]
  16. Li, J.; Lu, J.; Lu, H.; Wei, H.; Wu, B.; Jiang, P. A novel conservative dynamics envelope of vehicle lateral stability under different tire characteristics. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2023, 09544070231206207. [Google Scholar] [CrossRef]
  17. Jin, Z.L.; Liang, X.H.; Zhao, W.Z. Multi-gains Ratio of Steer-by-wire System and Anti-rollover Control for Vehicle. J. Mech. Eng. 2020, 56, 172–180. [Google Scholar]
  18. Zhang, L.; Ai, W.; Liu, Z.; Zhang, L.; Teng, X. Stability control of energy saving electric vehicle using dynamic nonlinear system state estimation. Fractals 2022, 30, 2240096. [Google Scholar] [CrossRef]
  19. Junru, Y. A Modeling and Factor Analysis of Safe Speed on Curve Segments. J. Transp. Inf. Saf. 2018, 36, 1–8. [Google Scholar]
  20. Liang, Z.; Zhao, J.; Liu, B.; Wang, Y.; Ding, Z. Velocity-based path following control for autonomous vehicles to avoid exceeding road friction limits using sliding mode method. IEEE Trans. Intell. Transp. Syst. 2022, 23, 1947–1958. [Google Scholar] [CrossRef]
  21. Liu, X.M.; Liu, F.; Guo, J.B. Research on stability of high-speed autonomous vehicles based on linear time-varying model predictive control. Int. J. Veh. Des. 2023, 91, 360–383. [Google Scholar] [CrossRef]
  22. Yao, X.; Gu, X.; Jiang, P. Coordination control of active front steering and direct yaw moment control based on stability judgment for AVs stability enhancement. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 59–74. [Google Scholar] [CrossRef]
  23. Xiang, C.; Peng, H.; Wang, W.; Li, L.; An, Q.; Cheng, S. Path tracking coordinated control strategy for autonomous four in-wheel-motor independent-drive vehicles with consideration of lateral stability. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2021, 235, 1023–1036. [Google Scholar] [CrossRef]
  24. Chen, L.; Xie, Y.; Cai, Y.; Sun, X.; Teng, C.; Zou, K. Stable Tracking Control of Autonomous Vehicles at Extreme Conditions. Automob. Eng. 2020, 42, 1016–1026. [Google Scholar]
  25. Zhong, L.F.; Peng, Y.H.; Jiang, M. Stability control of distributed drive electric vehicle based on phase plane. Automob. Eng. 2021, 43, 721–729. [Google Scholar]
  26. Xuan, N.; Liao, X.Y. Handling and Stability Control of Four-wheel Drive Electric Vehicle under Combined Slip Conditions. J. Mech. Eng. 2021, 57, 205–220. [Google Scholar]
  27. Park, G.; Choi, S.B. A model predictive control for path tracking of electronic-four-wheel drive vehicles. IEEE Trans. Veh. Technol. 2021, 70, 11352–11364. [Google Scholar] [CrossRef]
  28. Wang, K.; Yang, M.; Li, Y.; Liu, Z.; Wang, H.; Ding, W. Multidirectional motion coupling based extreme motion control of distributed drive autonomous vehicle. Sci. Rep. 2022, 12, 13203. [Google Scholar] [CrossRef] [PubMed]
  29. Jia, F.J.; Jing, H.H.; Gu, M.Q. Cooperative control of yaw and roll motion for in-wheel motor vehicle with semi-active suspension. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2022, 236, 3–15. [Google Scholar] [CrossRef]
  30. Hajiloo, R.; Khajepour, A.; Zengin, H.; Kasaiezadeh, A.; Chen, S.K. A coupled force predictive control of vehicle stability using front/rear torque allocation with experimental verification. Veh. Syst. Dyn. 2022, 60, 2541–2563. [Google Scholar] [CrossRef]
  31. Jan, B.; Tomáš, H. Decoupling of vehicle lateral dynamics using four-wheel steering system. In Proceedings of the 23rd International Conference on Process Control (PC), Strbske Pleso, Slovakia, 9 June 2021. [Google Scholar]
  32. Liu, H.; Liu, C.; Han, L.; Xiang, C. Handling and stability integrated control of AFS and DYC for distributed drive electric vehicles based on risk assessment and prediction. IEEE Trans. Intell. Transp. Syst. 2022, 23, 23148–23163. [Google Scholar] [CrossRef]
  33. Hang, P.; Xia, X.; Chen, X.B. Handling Stability Advancement With 4WS and DYC Coordinated Control: A Gain-Scheduled Robust Control Approach. IEEE Trans. Veh. Technol. 2021, 70, 3164–3174. [Google Scholar] [CrossRef]
  34. Lei, Y.; Wen, G.; Geng, X.; Fu, Y.; Li, X.; Hu, J. Comprehensive control of steering and braking via a 4WIS-4WID vehicle. J. Chin. Soc. Mech. Eng. 2022, 43, 303–313. [Google Scholar]
  35. Xu, T.; Zhao, Y.; Deng, H.; Guo, S.; Li, D.; Lin, F. Integrated optimal control of distributed in-wheel motor drive electric vehicle in consideration of the stability and economy. Energy 2023, 282, 128990. [Google Scholar] [CrossRef]
  36. Wu, X.; Zhou, B.; Wu, T.; Pan, Q. Research on intervention criterion and stability coordinated control of AFS and DYC. Int. J. Vehicle Des. 2022, 90, 116–141. [Google Scholar] [CrossRef]
  37. Li, W.; Yu, S.; Tan, L.; Li, Y.; Chen, H.; Yu, J. Integrated control of path tracking and handling stability for autonomous ground vehicles with four-wheel steering. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2023, 09544070231204249. [Google Scholar] [CrossRef]
  38. Tang, J.; Li, Q.; Wang, Y.; He, D.; Chen, Y. Tracking Control for Autonomous Four-Wheel Independently Driven Vehicle Based on Deep Reinforcement Learning. In Proceedings of the 6th CAA International Conference on Vehicular Control and Intelligence (CVCI), Nanjing, China, 8 December 2022. [Google Scholar]
  39. Jeong, Y.; Kim, D.H.; Youn, M.H.; Park, S.; Li, X.J.; Lee, T. Model Predictive Control-Based Path Tracking with Four-Wheel Independent Steering, Driving, and Braking Autonomous Vehicles on Low Friction Road. In Proceedings of the 23rd International Conference on Control, Automation, and Systems, Yeosu, Republic of Korea, 20 November 2023. [Google Scholar]
  40. Huang, Y.; Yong, S.Z.; Chen, Y. Stability control of autonomous ground vehicles using control-dependent barrier functions. IEEE Trans. Intell. Veh. 2021, 6, 699–710. [Google Scholar] [CrossRef]
  41. Chen, J.; Shuai, Z.; Zhang, H.; Zhao, W. Path following control of autonomous four-wheel-independent-drive electric vehicles via second-order sliding mode and nonlinear disturbance observer techniques. IEEE Trans. Ind. Electron. 2020, 68, 2460–2469. [Google Scholar] [CrossRef]
  42. Jeong, Y.; Yim, S. Path tracking control with four-wheel independent steering, driving and braking systems for autonomous electric vehicles. IEEE Access 2022, 10, 74733–74746. [Google Scholar] [CrossRef]
  43. GB/T 6323-2014; Vehicle Handling and Stability Test Method. China Standards Press: Beijing, China, 2014.
  44. QC/T 480-1999; Vehicle Handling and Stability Index Limits and Evaluation Methods. China National Machinery Industry Bureau: Beijing, China, 1999.
  45. ISO 3888-2: 2002; Double-Line International Standard Experiment. International Organization for Standardization: Geneva, Switzerland, 2002.
  46. Dugoff, H.P.; Fancher, P.S.; Segel, L. An Analysis of Tire Traction Properties and Their Influence on Vehicle Dynamic Performance. SAE Trans. 1970, 79, 1219–1243. [Google Scholar]
  47. Zhang, L.P.; LI, L.; Qi, B.N. Rollover prevention control for a four in-wheel motors drive electric vehicle on an uneven road. Sci. China (Technol. Sci.) 2018, 61, 934–948. [Google Scholar] [CrossRef]
  48. Wang, K.; Ding, W.; Yang, M.; Zhu, Q. Dynamic-boundary-based lateral motion synergistic control of distributed drive autonomous vehicle. Sci. Rep. 2021, 11, 22644. [Google Scholar] [CrossRef] [PubMed]
  49. Yu, Z.S. Vehicle handling and stability. In Automotive Theory, 6th ed.; Machine Press: Beijing, China, 2009. [Google Scholar]
  50. Liu, J.K. Introduction to sliding mode variable structure control. In Sliding Mode Variable Structure Control MATLAB Simulation: Basic Theory and Design Method, 3rd ed.; Tsinghua University Press: Beijing, China, 2015. [Google Scholar]
Figure 1. Illustration of 8-DOF vehicle dynamics model.
Figure 1. Illustration of 8-DOF vehicle dynamics model.
Electronics 13 02744 g001
Figure 2. Response surface model of roll moment M z s under uniform linear condition.
Figure 2. Response surface model of roll moment M z s under uniform linear condition.
Electronics 13 02744 g002
Figure 3. Response surface model of roll moment M z s under a right sharp turn condition.
Figure 3. Response surface model of roll moment M z s under a right sharp turn condition.
Electronics 13 02744 g003
Figure 4. Architecture of MMCCS for a wheel drive driverless vehicle.
Figure 4. Architecture of MMCCS for a wheel drive driverless vehicle.
Electronics 13 02744 g004
Figure 5. The target path.
Figure 5. The target path.
Electronics 13 02744 g005
Figure 6. (a) The various speeds correspond to the dry road. (b) The relevant speed corresponds to the ice-snow road.
Figure 6. (a) The various speeds correspond to the dry road. (b) The relevant speed corresponds to the ice-snow road.
Electronics 13 02744 g006
Figure 7. (a) The lateral acceleration corresponds to various speeds under dry road conditions. (b) The lateral acceleration corresponds to various speeds under ice-snow road conditions.
Figure 7. (a) The lateral acceleration corresponds to various speeds under dry road conditions. (b) The lateral acceleration corresponds to various speeds under ice-snow road conditions.
Electronics 13 02744 g007
Figure 8. Comparison of the path tracking error. (a) The speed tracking error. (b) Two algorithms on the dry road.
Figure 8. Comparison of the path tracking error. (a) The speed tracking error. (b) Two algorithms on the dry road.
Electronics 13 02744 g008
Figure 9. Comparison of the path tracking error. (a) The speed tracking error. (b) Two algorithms on the ice-snow road.
Figure 9. Comparison of the path tracking error. (a) The speed tracking error. (b) Two algorithms on the ice-snow road.
Electronics 13 02744 g009
Figure 10. The comparison of the path tracking error. (a) The speed tracking error. (b) The dry road.
Figure 10. The comparison of the path tracking error. (a) The speed tracking error. (b) The dry road.
Electronics 13 02744 g010
Figure 11. The comparison of the yaw rate. (a) The sideslip angle (b) The dry road.
Figure 11. The comparison of the yaw rate. (a) The sideslip angle (b) The dry road.
Electronics 13 02744 g011
Figure 12. The comparison of the roll angle. (a) The dry road and the path tracking error. (b) The ice-snow road.
Figure 12. The comparison of the roll angle. (a) The dry road and the path tracking error. (b) The ice-snow road.
Electronics 13 02744 g012
Figure 13. The comparison of the speed tracking error. (a) The yaw rate. (b) The ice-snow road.
Figure 13. The comparison of the speed tracking error. (a) The yaw rate. (b) The ice-snow road.
Electronics 13 02744 g013
Figure 14. The comparison of the sideslip angle. (a) The roll angle. (b) The ice-snow road.
Figure 14. The comparison of the sideslip angle. (a) The roll angle. (b) The ice-snow road.
Electronics 13 02744 g014
Figure 15. Architecture of the testbed software platform.
Figure 15. Architecture of the testbed software platform.
Electronics 13 02744 g015
Figure 16. The testbed of MMCCS.
Figure 16. The testbed of MMCCS.
Electronics 13 02744 g016
Figure 17. The comparison of testbed test results and digital simulation test results. (a) The path tracking error. (b)The speed tracking error on the dry road.
Figure 17. The comparison of testbed test results and digital simulation test results. (a) The path tracking error. (b)The speed tracking error on the dry road.
Electronics 13 02744 g017
Figure 18. The comparison of testbed test results and digital simulation test results. (a) The yaw rate. (b) The sideslip angle on the dry road.
Figure 18. The comparison of testbed test results and digital simulation test results. (a) The yaw rate. (b) The sideslip angle on the dry road.
Electronics 13 02744 g018
Figure 19. The comparison of testbed test results and digital simulation test results. (a) shows the roll angle on the dry road. (b) The path tracking error on the ice-snow road.
Figure 19. The comparison of testbed test results and digital simulation test results. (a) shows the roll angle on the dry road. (b) The path tracking error on the ice-snow road.
Electronics 13 02744 g019
Figure 20. The comparison of testbed test results and digital simulation test results. (a) The speed tracking error. (b) The yaw rate on the ice-snow road.
Figure 20. The comparison of testbed test results and digital simulation test results. (a) The speed tracking error. (b) The yaw rate on the ice-snow road.
Electronics 13 02744 g020
Figure 21. The comparison of testbed test results and digital simulation test results. (a) The sideslip angle. (b) The roll angle on the ice-snow road.
Figure 21. The comparison of testbed test results and digital simulation test results. (a) The sideslip angle. (b) The roll angle on the ice-snow road.
Electronics 13 02744 g021
Table 1. Main WDDV parameters.
Table 1. Main WDDV parameters.
SymbolDescriptionValue (Unit)
mVehicle mass1798 (kg)
m f springing mass of the front suspension150 (kg)
m r springing mass of the rear suspension150 (kg)
a Distance from center of gravity to front axle1.050 (m)
bDistance from center of gravity to rear axle1.620 (m)
d f Front wheel track1.565 (m)
d r Rear wheel track1.565 (m)
h c Height of vehicle center of gravity0.650 (m)
I x Roll moment of inertia of vehicle700.7 (kg·m2)
I y Pitch moment of inertia of vehicle2059.2 (kg·m2)
I z Yaw moment of inertia of vehicle2059.2 (kg·m2)
J w Rotational inertia of the wheel0.85 (kg·m2)
RWheel rolling radius0.347 (m)
Table 2. Summary of simulation results of double-shift line driving and serpentine driving under conventional conditions.
Table 2. Summary of simulation results of double-shift line driving and serpentine driving under conventional conditions.
Control SystemDCSMMCCSImprovement Rate
μ Parameters (Units)MinMaxMinMax
0.8 e P (m)−0.0290.029−0.0240.02417%
e s (m/s)−0.0010.007−0.0010.00443%
0.15 e P (m)−0.0130.013−0.0090.00931%
e s (m/s)0.0040.0060.0030.00517%
The improvement rate in the table is calculated relative to the DCS control results and is calculated according to the larger absolute value of the upper and lower boundaries of the corresponding parameter variation range. For example, e s for the road adhesion coefficient is 0.8, from the data in the table: abs (−0.001) < abs (0.007), and abs (−0.001) < abs (0.004), so the optimized ratio is [abs (0.007) − abs (0.004)]/abs (0.007) 43%. The improvement rate is positive, indicating performance improvement.
Table 3. Summary of simulation results of double-shift line driving and serpentine driving under extreme conditions.
Table 3. Summary of simulation results of double-shift line driving and serpentine driving under extreme conditions.
Control SystemDCSMMCCSImprovement Rate
μ Parameters (Units)MinMaxMinMax
0.8 e P (m)−0.0720.072−0.060.0617%
e s (m/s)−0.0290.014−0.0190.00334%
γ (rad/s)−0.4740.474−0.4330.4329%
β (rad)−0.0510.051−0.0440.04214%
ϕ (rad)−0.0580.058−0.050.0514%
0.15 e P (m)−0.0930.093−0.0620.06133%
e s (m/s)0.0030.0060.0030.0060%
γ (rad/s)−0.1970.197−0.1760.17511%
β (rad)−0.0240.024−0.0220.0228%
ϕ (rad)−0.010.01−0.0090.00910%
The improvement rate in the table is calculated relative to the DCS control results and is calculated according to the larger absolute value of the upper and lower boundaries of the corresponding parameter variation range. For example, e s with a road adhesion coefficient of 0.8, from the data in the table: abs (−0.029) > abs (0.014), and abs (−0.019) > abs (0.003), so the optimization ratio is [abs (−0.029) − abs (−0.019)]/abs (−0.029) 34%. The positive improvement rate indicates improved performance.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, K.; Luo, Y.; Du, L.; Wu, Z.; Wang, H. Wheel Drive Driverless Vehicle Handling and Stability Control Based on Multi-Directional Motion Coupling. Electronics 2024, 13, 2744. https://doi.org/10.3390/electronics13142744

AMA Style

Wang K, Luo Y, Du L, Wu Z, Wang H. Wheel Drive Driverless Vehicle Handling and Stability Control Based on Multi-Directional Motion Coupling. Electronics. 2024; 13(14):2744. https://doi.org/10.3390/electronics13142744

Chicago/Turabian Style

Wang, Kai, Yi Luo, Lifang Du, Zhongping Wu, and Han Wang. 2024. "Wheel Drive Driverless Vehicle Handling and Stability Control Based on Multi-Directional Motion Coupling" Electronics 13, no. 14: 2744. https://doi.org/10.3390/electronics13142744

APA Style

Wang, K., Luo, Y., Du, L., Wu, Z., & Wang, H. (2024). Wheel Drive Driverless Vehicle Handling and Stability Control Based on Multi-Directional Motion Coupling. Electronics, 13(14), 2744. https://doi.org/10.3390/electronics13142744

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop