Abstract
Lateral stability is vital to vehicle handling stability and traffic safety. It is also a crucial factor for the path-tracking ability of the vehicle in the intelligent transportation system (ITS). Most of the research focuses on vehicles with an Ackerman steering system. The ASV (articulated steering vehicle) has a lower steering radius. Thus, it is widely used in some special vehicles, such as mining and construction vehicles. The ASV has weaker lateral stability than the vehicle with an Ackerman steering system. To improve the stability of an ASV, the nonlinear lateral dynamic model is established and validated by field test. With the lateral dynamic model, the steering characteristic of the ASV is analyzed. Based on the stability criteria analysis, the vehicle sideslip angle and angular velocity phase portrait are chosen as the stability indicator. An integrated AASS (active articulated steering system) DYC (direct yaw control) controller based on the adaptive MPC (model predictive control) method is designed according to the track on the phase plane. The double lane-change and 0.7 Hz sine with dwell maneuver are initiated based on the integrated vehicle dynamic model. The results suggest that the provided controller has a better stability performance than the current antiwindup PID control algorithm. It lays a good foundation for the vehicle safety and path tracking of ASV in the ITS.
1. Introduction
ASV (articulated steering vehicle) is a vehicle with an articulated steering system. Normally, the front and rear parts of the ASV form an angle to allow the vehicle to change its direction as the driver’s intention. It is widely used in the engineering area, such as loaders and mining vehicles. Developing the intelligent transportation system (ITS) for these vehicles to ease the burden of the driver or liberate them from crucial work conditions is necessary and meaningful. Without the participation of the pedestrian and with the enclosed working condition, the scenario of these ASV in an intelligent transportation system is simpler than the passenger cars. However, the dynamic characteristics of the ASV are quite different from the vehicle with the Ackerman steering system due to its unique mechanical structure. The advantage of an ASV is that it has a lower radius during the steering process [1]. The weakness of the ASV is that it has relatively weaker lateral stability compared to the vehicle with the Ackerman steering system [2].
Plenty of work has been carried out to strengthen the vehicle’s lateral stability, most of them focusing on the vehicle with an Ackerman steering system. Jin et al. provide an integrated AFS (active front steering) and ASS (active suspension system) controller to improve the vehicle handling stability based on a gain-scheduling state-feedback controller [3]. A full X-by-wire chassis coordinated control scheme utilizing the AFS, DYC (direct yaw control), ARC, and ASR is provided to improve the vehicle’s longitudinal, yaw, and roll stability synthetically [4]. Aouadj et al. set up a controller to coordinate AFS and DYC based on sliding mode. The AFS and DYC system is designed based on the phase plane portrait [5]. The robust Takagi–Sugeno fuzzy control is initiated to combine the AFS and DYC systems [6]. Besides, MPC (model predictive control) [7] and NMPC (nonlinear model predictive control) [8] are also applied in this area. The good lateral stability also lays a good foundation for the autonomous driving vehicle [9]. The vehicle with the distribution driving system equipped with lots of sensors, the influence of poor SNR (signal-to-noise ratio), and the fault of the sensor are also analyzed and handled by the researchers [10, 11].
For the ASV, the current research focuses on the effect of the hydraulic steering system on the vehicle dynamics [12, 13], lateral stability influence factor analysis [14], improving the vehicle stability by the control of ASS [15], and antirollover system [16].
A typical architecture of the controller for vehicle stability control is shown in Figure 1. The reference model is a 2 or 3 DOFs (degree of freedoms) that provides the ideal yaw rate and vehicle sideslip angle of the vehicle [17]. Normally, the multilayer controller is used. The upper controller estimates the vehicle state and generates the active steering angle and torque, and the lower controller distributes the torque to the motors. Gao et al. provide a vehicle-based method to estimate the attitude of underground ASV [18]. For the torque distribution algorithms, generalized inversion [19] and quadratic programming methods [20] are provided. The active steering angle and torque are sent to a vehicle dynamic model with 8, 10, or 14 DOFs to verify the effect of the controller.

To improve the lateral stability of an ASV, the AASS (active articulated steering system) is combined with the DYC system. The reference model with 2 DOFs and an integrated simulation vehicle model is built and verified by the field test result. Based on the nonlinear reference model, the boundary of the ideal yaw rate and vehicle sideslip angle is determined. The integrated controller is designed based on the vehicle stability status judged by the vehicle sideslip angle and sideslip angular velocity phase portrait. The torque of the driving system and active articulated angle is coordinated by the MPC method. The provided controller is verified by comparing the vehicle performance to the reference model, without control, and with a current antiwindup PID controller.
2. Vehicle Lateral Dynamic Model
For the lateral dynamics control of the vehicle, it is not common to use a complex vehicle and dynamic model. Normally, a simplified lateral dynamic model with 2 DOFs is used to design and implement the controller [21]. A typical 2 DOFs model is shown in Figure 2. Since this model mainly focuses on the lateral dynamic features, other motions, such as roll, pitch, and vertical, are neglected. Due to its unique mechanic structure, the lateral and yaw motion of the ASV can be described by equations (1) and (2), respectively.

In equations (1) and (2), means the mass, represents the lateral velocity, is the tire force, means the yaw moment of the inertia, the footnote means the front and rear vehicle body, represents the left and right side of the vehicle, and are the distance between the CG (center of gravity) to the axle and articulation point, is the vehicle yaw rate, is the wheelbase, and means the external yaw moment. The detail of the annotation in this paper is listed in the abbreviations. The relationship between the lateral velocity and longitudinal velocity can be established by the following equation:where is the vehicle sideslip angle. The relationship between the lateral force and tire sideslip angle is nonlinear. Therefore, a nonlinear tire model is necessary to calculate the lateral force. A brush tire model variant of the Fiala nonlinear brush model is used in this model. It assumes one coefficient of friction and parabolic force distribution [22]. In the brush model, the lateral force, , can be expressed as follows:where is the slip rate of the wheel, represents the road adhesion coefficient, and and are the longitudinal and lateral stiffness of the tire. The tire sideslip angle, , can be gained by the following equation:
Here, means the articulation angle. The normal load of the wheels on the front and rear axle, and , can be defined as follows:where is the sideslip angle that corresponds to the full saturation of the lateral tire force.
Based on equations (1)–(10), the state-space representation of the lateral dynamic model can be written as follows:where is the state vector, is the system matrix, and are the control input matrixes, and and are the control input vectors.
3. Field Test and Model Validation
To validate the established lateral vehicle dynamic model and integrated simulation model built by MATLAB/Simulink and Recudyn, a field test is initiated. A speed sensor and IMU (inertial measurement unit) are mounted on the vehicle to detect the velocity and yaw rate during the test. The articulation angle is collected by a distance sensor, and based on the geometry transformation, the steering angle during the test can be acquired. A microcontrol unit, STM 32, is used to control the ASV. The layout of the sensors is demonstrated in Figure 3. The detail of the sensors is demonstrated in Table 1.

The steady steering maneuver is used during the field test, and the yaw rate and vehicle trajectory acquired during the test is compared to the value calculated by the established lateral dynamic model and integrated model (shown in Figures 4 and 5).


According to the demonstration in Figures 4 and 5, the calculated yaw rate and vehicle routine are in good correlation with the test result. In summary, the models could represent the dynamic characteristics accurately. Thus, the established models can be used in further analysis.
4. Integrated Controller Design
The gap between the ideal yaw rate and the actual yaw rate of the vehicle is a vital indicator to access vehicle stability [23]. The ideal yaw rate of the vehicle is determined by the vehicle speed and road adhesion coefficient. Normally, the vehicle’s ideal yaw rate is acquired by the steady yaw rate gain. According to equation (11), the yaw rate gain of ASV can be gained.where is the stability factor, which can be gained by the following equation:
Since the lateral acceleration of the vehicle cannot exceed the maximum friction force between the tire and the road surface. When the vehicle is in a steady steering state, the boundary of lateral acceleration can be expressed as follows:
In summary, the idea yaw rate should meet the following equation:
According to equation (16), the yaw rate of the ASV varies with the articulation angle, and the road adhesion coefficient is illustrated in Figure 6.

As shown in Figure 6, the ideal yaw rate increases with the vehicle speed until the yaw rate reaches the boundary determined by the road adhesion coefficient. For example, when the vehicle runs on the road with a road adhesion coefficient of 0.2, the articulation angle is 2 degrees. The ideal yaw rate increases with the speed when the speed is in the range of 0 to 10 m/s and decreases with it when the speed is higher than 10 m/s.
Besides the yaw rate, the vehicle sideslip angle is another indicator of vehicle stability. According to the work of Eslamian et al. [24], based on equation (11), the steady vehicle sideslip angle during the constant cornering process is obtained as follows:
The maximum sideslip angle is also related to the road adhesion coefficient, it can be defined by the following equation:
The vehicle sideslip angle corresponding to different articulation angles and road adhesion coefficient is shown in Figure 7.

As shown in Figure 7, the vehicle sideslip angle is inversely proportional to the articulation angle. Thus, the ideal vehicle sideslip angle can be expressed as follows:
According to the illustration in Figure 7, the vehicle runs on the road with a lower road adhesion coefficient has a smaller region. To evaluate the vehicle stability, normally, the phase portrait, such as the , , or tire sideslip phase portrait is chosen as the indicator. The phase portrait is widely used in the stability assessment [12, 25–27]. Figure 8 shows the phase portrait of the ASV. In Figure 8, the zone between the red dash lines is the stable area; otherwise, it is the unstable region.

The boundary of the stable region (the red dashed line in Figure 8) can be described by the following equation [28, 29]:
The parameter and are acquired based on the stability boundary of the phase portrait. When the vehicle is approaching the stability region, the vehicle will rush into an unstable state when the AFS system works. Thus, the control algorithm of AASS and DYC can be identified as Table 2.
Normally, the ideal vehicle sideslip angle is set as zero [30, 31]. However, setting the vehicle sideslip angle as zero causes the DYC system initiated frequently. Based on equations (16) and (19), the ideal vehicle yaw rate and sideslip angle can be gained.
The state-space of the prediction model can be expressed as follows:
The state vector of the system is defined as follows:where and are the lateral force of the front and rear axle and is the error between the actual and ideal wheel speed. The input vector and disturbance vector of the system arewhere and are the torque and disturbance force of the motor and is the acceleration in the longitudinal direction.
The matrix to describable the system can be defined as follows:where , , and .
The state-space function of the discrete system can be gained based on the state-space formula of the continuous system.where the footnotes and are the state vectors at the step and . The output of the system, , is defined as follows:
The system matrix can be gained based on equations (24)∼(27).where is a unit matrix and means the time step.
The steering angle and vehicle speed are assumed as constant during the to step. Thus, the optimal solution of the model prediction control can be gained.where means the predicted time step, is the ideal state, is the torque distribution of the driving wheels. , , and are the weight matrix. is the predicted output sequence that meets the minimum objective .
According to the QP (quadratic problem) method and constraint of the motor and tire, can be gained and taken as the output of the MPC controller.
Two independent input is needed to track two system status (, ). There are two inputs in this system, but they are not independent. The importance of the system status is determined by the weight coefficient, and , in the matrix . Normally, the yaw rate focuses on the vehicle handling performance while the vehicle sideslip angle focuses on the stability. When the vehicle is in a stable region, better handling performance is expected. When the vehicle is in an unstable region, stability performance is more important. Thus, to improve the performance of the controller, the weight coefficient, and , is determined by the stability criteria. Therefore, this MPC controller is adaptive to the vehicle stability criteria.
5. Discussion
For the assessment of the vehicle’s lateral stability, there are some classic maneuvers, such as steady steering, snake steering, fishhook, double lane-change, and 0.7 sine with dwell maneuver. Among them, the double lane-change and fishhook are the most common conditions in the relevant research. They are chosen to compare the lateral stability performance with the provided control algorithm and an antiwindup PID controller proposed in the paper [32]. The double lane-change is initiated based on the regulation in ISO-3888-1-2018 (passenger cars-test track for a severe lane-change maneuver: obstacle avoidance), and 0.7 sine with dwell maneuver is regulated in Federal Motor Vehicle Safety Standard (FMVSS) No. 126, Electronic Stability Control (ESC) System.
An HIL (hardware in loop) platform is used to verify the lateral stability control algorithm. As shown in Figure 9, the platform is made up of three parts, data visualization system, dSPACE Micro AutoBox, and SCALEXiO. The data visualization system displays the data from the Micro Auto Box and SCALEXiO and input the driver’s intention. The reference model, state estimator, and stability assessment and control algorithm is integrated in the Micro AutoBox. It acts as the VCU (vehicle control unit) in a practical vehicle. The vehicle dynamic model, motor controller model, road model is set in the SCALEXiO.

5.1. Double Lane-Change
The vehicle should shift from the original to the paralleled one and back to the original one without hitting the barriers during the double lane-change. The barriers are regulated in ISO-3888-1-2018. The vehicle yaw rate and sideslip angle of the reference model, without control method, with the provided adaptive MPC controller and antiwindup PID controller, are compared in Figures 9 and 10. The output of the reference model is taken as the baseline.

In Figures 10 and 11, with the provided adaptive MPC controller, the yaw rate and vehicle sideslip angle have the minimum error. Without the controller, the yaw rate and vehicle sideslip have the maximum gap from the baseline. The performance of the antiwindup PID controller is between the provided controller and without any control.

The vehicle trajectory is illustrated in Figure 12. Without control, the vehicle fails to pass the maneuver for contacting the upper boundary. With the adaptive MPC and antiwindup PID controller, the vehicle can pass the test. The minimum trajectory error occurs in the vehicle with an adaptive MPC controller.

5.2. 0.7 Sine with Dwell Maneuver
The steering input of the 0.7 sine with dwell maneuver is illustrated in Figure 13. It is widely used in this area due to its objectivity, practicability, repeatability, and representativeness.

In this test, the two lateral criteria are defined as the ratio of the vehicle yaw rate at a specified time to the first local peak yaw rate generated by 0.7 Hz sine with dwell steering reversal. They can be represented in the mathematical notations as follows:where means the yaw rate at the time , is the first local peak yaw rate generated by the 0.7 Hz sine with dwell steering input, and represents the time to completion of steering input.
Figure 14 shows the yaw rate of the four models. The absolute value of the reference peak yaw rate is lower than the other three models. Without control, the model has the maximum peak value. The yaw rate in equations (35) and (36) can be extracted from Figure 14 and is shown in Table 3.

According to equations (35) and (36), all of three models meet the lateral criteria. Because of the low value, the and of the adaptive MPC controller is higher than the PID controller. The criteria of the model with PID controller have the minimum criteria, which means it has the fastest yaw rate decrease rate. With the adaptive MPC controller, the model has the minimum peak value but the lowest decrease rate. The variation of vehicle sideslip angle is compared in Figure 15.

As shown in the comparison, the peak value of the sideslip angle of the three models is higher than the reference model. Taking the reference model as the baseline, the model without control has the maximum deviation value, and the model with adaptive MPC has the minimum deviation value. All of the peak values of the three models fluctuate, and the fluctuation range gets smaller.
6. Conclusion
The ASV is widely used in the special and engineering vehicle area due to its unique characteristics. It also has the disadvantage of weak lateral stability. To deal with its shortcoming, the lateral dynamic model of the ASV is established based on the dynamic theory. The built model is validated by the field test result. The state-space formula of the ASV is more complicated than the vehicle with an Ackerman steering system. The corresponding ideal yaw rate and vehicle sideslip angle expression also differ from the regular car, especially for the ideal vehicle sideslip angle. For the ASV, the ideal vehicle sideslip angle and articulated angle have the opposite sign; however, for the regular vehicle, the sideslip angle and steering angle have the same sign.
According to the lateral stability analysis, an integrated AAFS and DYC based on MPC theory is provided based on the vehicle sideslip angle and sideslip angular velocity phase portrait. The vehicle performance with the provided adaptive MPC controller and antiwindup PID controller during two typical lateral stability tests is compared. According to the comparison, in the double lane-change test, the vehicle with MPC and PID controller can help the vehicle to pass the test, the vehicle with MPC has a more ideal routine. In the 0.7 Hz sine with dwell maneuver test, the vehicle with MPC has the minimum peak value and deviation value after the steering process, and the vehicle with PID control has the maximum recovery rate.
With the provided integrated controller, the lateral stability of the ASV is improved. Hence, the vehicle safety and path-tracking ability are strengthened, which will benefit the construction of the ITS with the ASV.
Abbreviations
: | Front vehicle body mass |
: | Rear vehicle body mass |
: | Front vehicle body lateral velocity |
: | Rear vehicle body lateral velocity |
: | Lateral force of front left tire |
: | Lateral force of rear left tire |
: | Lateral force of front right tire |
: | Lateral force of rear right tire |
: | Front vehicle body yaw moment of the inertia |
: | Rear vehicle body yaw moment of the inertia |
: | Distance between the front axle and CG of front vehicle body |
: | Distance between the articulation point and CG of front vehicle body |
: | Distance between the rear axle and CG of rear vehicle body |
: | Distance between the articulation point and CG of rear vehicle body |
: | Vehicle yaw rate |
: | Vehicle sideslip angle |
: | Longitudinal stiffness of the tire |
: | Lateral stiffness of the tire |
: | Slip rate of wheel |
: | Road adhesion coefficient |
: | Side slip angle of tire |
: | Articulation angle |
: | Front axle load |
: | Rear axle load |
: | Sideslip angle that corresponds to the full saturation of the lateral tire force |
: | State vector of state-space |
: | System matrix |
: | Control input matrix for control vector 1 |
: | Control input matrix for control vector 2 |
: | Control vector 1 |
: | Control vector 2 |
: | Ideal yaw rate |
: | Stability factor |
: | Maximum ideal yaw rate |
: | Desired ideal yaw rate |
: | Ideal vehicle sideslip angle |
: | Maximum ideal vehicle sideslip angle |
: | Desired ideal vehicle sideslip angle |
: | Fitted parameter for stability boundary function |
: | Fitted parameter for stability boundary function |
: | State vector of the prediction model state-space |
: | Error between the actual and ideal wheel speed |
: | Input vector of the prediction model state-space |
: | Disturbance vector of the prediction model state-space |
: | System matrix of the prediction model state-space |
: | Control matrix of the prediction model state-space |
: | Motor torque |
: | Motor disturbance force |
: | Predicted output sequence that meets the minimum objective |
: | First sequence of the predicted output sequence. |
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This study was funded by the Innovative Research Team Development Program of Ministry of Education of China (Grant no. IRT_17R83) and the 111 Project (Grant no. B17034) of China.