Abstract

This study examines the roll instability mechanism and stability index of articulated steering vehicles (ASVs) by taking wheel loaders as the research object. A seven-degree-of-freedom nonlinear dynamics model of the ASVs is built on the basis of multibody dynamics. A physical prototype model of an ASV is designed and manufactured to validate the dynamic model. Test results reasonably agree with the simulation results, which indicates that the established dynamic model can reasonably describe ASV movements. Detailed analysis of the rollover stability of the wheel loader is performed with the use of the established dynamic model. Analysis results show that rollover will occur when the roll angular velocity exceeds a critical threshold, which is affected by lateral acceleration and slope angle. On this basis, a dynamic stability index applicable to the ASVs is presented.

1. Introduction

Both maneuverability and traction of off-road vehicles are important for enhanced performance in applications on unpredictable and changing terrains. To combine these two important features, articulated frame steering has become popular, and an increasing number of manufacturers choose this simple but efficient steering system. Many wheel loaders, scrapers, and forestry skidders are articulated steering vehicles (ASVs).

Modern ASVs, such as wheel loaders, have two separate parts that are connected by a steering hinge. The wheel axles are unsuspended, such that the rear axle is mounted on a pivoting joint to enable wheels contact on uneven ground. A schematic diagram of a wheel loader is shown in Figure 1. This structure is limited by decreased stability when turning because the center of gravity (CG) is laterally displaced. In addition, ASVs are off-road vehicles that work on steep slopes, unexpected banks, or uneven terrain, which result in a high probability of rollover accident.

Many scholars have focused on the stability and dynamics of such vehicles as cars, heavy trucks, and tractors. However, limited work has been conducted on the stability and handling performance of the ASVs, especially in terms of their motions on slopes and uneven terrain. Azad et al. [1] built a linear ASV model that can be used to analyze snaking instability. Yavin [2] built a simple dynamical model of the ASV on a horizontal plane. The vertical movement and the relative movement between the rear body and the rear axle were not considered in the model, such that it cannot be used to analyze instability when the ASVs pass over obstacles. Corke and Ridley [3] and Dragt et al. [4] studied the steering kinematics of the ASVs on the basis of plane models that do not consider vertical movement.

To evaluate vehicle stability quantitatively, many stability indexes have previously been developed for given conditions or certain prototype vehicles. Some existing stability indexes [57] are only based on acceleration or roll angular velocity threshold values. Given that the ASVs are off-road vehicles that travel on complex roads and work under poor conditions, vehicle roll stability cannot be accurately reflected if only a single lateral acceleration or roll angular velocity is used as a stability index. For example, instant lateral acceleration of ASV is large when passing over obstacles, but the vehicle does not lose stability in this case. As another example, the small roll angular speed of ASVs may cause instability when a vehicle is traveling on sloping road.

Efforts have been exerted to predict vehicle rollover on the basis of vehicle roll energy and tip-over potential energy. Nalecz et al. [8] suggested an energy-based function called the rollover prevention energy reserve (RPER). PRER is defined as the difference between the energy needed to bring the vehicle to its tip-over position and the rotational kinetic energy, which can be transformed into gravitational potential energy to lift the vehicle. Given the relative rotation between the front and rear bodies of the ASVs, the moment of inertia and CG position of the vehicle are changing in real time, such that the PREP of the ASVs is difficult to calculate.

Liu and Ayers [9, 10] proposed a dynamic stability index composed of lateral angle and lateral angular velocity. This index does not consider the effect of lateral acceleration on the ASVs instability. In fact, the ASVs are likely to lose stability because of large lateral acceleration when turning.

Preston-Thomas and Woodrooffe [11] used the lateral load transfer ratio (LTR) to indicate rollover threat. LTR can vary from 0 when the loads carried by the tires are equal to 1 when the tires lift off the road. In the calculation of LTR, the contact pressure of each tire must be known. However, tire contact pressure is difficult to measure. To avoid the direct measurement of tire contact pressure, some scholars approximately calculate LTR on the basis of the suspension system [1214]. However, currently available ASVs have no suspension system.

By summarizing existing research, it can be found that the dynamic models of the ASVs were built without considering the impact of nonstructural terrain. In fact, sudden change in terrain is one of the main causes of instability. The existing models do not consider the impact of the swing bridge on the roll stability of the ASVs. In addition, the existing stability indexes are not entirely suitable for the ASVs.

Therefore, the physical model and stability index for the ASVs still need improvement. This research aims to build a nonlinear dynamic model that considers articulated steering movement, swing bridge movement, and ground environmental parameters to obtain the instability mechanism and dynamic rollover stability index for the ASVs.

2. Modelling of the ASVs

The typical ASVs comprise three parts, namely, the front body, the rear body, and the rear axle, as shown in Figure 1. Each part has three translational DOFs and three revolute DOFs, but all 18 DOFs are not completely independent. The longitudinal movements, the lateral movements, the vertical movements, and the pitch movements of the front body, the rear body, and the rear axle are interrelated. The roll motions of the front body and the rear body are also interrelated. The yaw movements of the rear body and the rear axle are related too. Meanwhile, the yaw movement of the front body is also related to that of the rear body. When the vehicle starts turning, the yaw rate difference between the front body and the rear body is equal to the speed value of the steering angle. When the steering input is finished, the yaw rate of the front body is equal to that of the rear body.

According to the structure characteristics of the swing bridge, the roll motion of the ASVs includes two phases.

Instability of Level 1. The swing bridge does not reach the limit position. The rear body and the rear axle roll over around the oscillating pin. There are relative movements between the two.

Instability of Level 2. The swing bridge reaches the limit position and comes into contact with the rear axle. There are no further relative movements between the rear axle and the rear body.

In addition, the ASVs, such as the wheel loaders, are the heavy mining and construction equipment, which function in loading, transporting, and discharging materials as displayed in Figure 2. Load weight and arm position are changing when the ASVs are loading and discharging materials, which have an impact on the stability of the ASVs. However, the forward speed is very slow or even stops when the ASVs are loading and discharging materials, as shown in Figures 2(a) and 2(c), and the risk of rollover is very low in these conditions. In fact, the ASVs are more likely to lose stability in the transporting state, as shown in Figure 2(b), and, in general, load weight and arm position remain unchanged in the transporting state. In order to simplify the analysis, this paper just considers the stability of the ASVs in the transporting state; that is, the load weight and arm position are assumed to be changeless.

Through the analysis, it is found that the following seven degrees of freedom can describe the movements of the ASVs: (1) the longitudinal velocity of the vehicle ; (2) the lateral velocity of the vehicle ; (3) the vertical movement of the vehicle ; (4) the pitch movement of the vehicle ; (5) the roll movement of the front body and the rear body ; (6) the roll movement of the rear axle ; (7) the yaw movements of the front body, the rear body, and the rear axle .

Figure 3 shows the movements, the structure dimensions of the ASVs, and the forces on tires. The reference point of the coordinate system of the vehicle coincides with the intersection point of the articulated steering pin axis and the oscillating pin axis, which is connected with the rear body. The -axis is parallel to the chassis plane of the rear body. It points to the forward direction of the perpendicular bisector of the rear axle. The positive direction of the -axis is perpendicular to the horizontal plane. The positive direction of the -axis is determined by the right-hand rule.

According to D’Alembert principle, the dynamics equations of the ASVs are written as follows:(1)Vehicle motion along the -axis direction:(2)Vehicle motion along the -axis direction:(3)Vehicle motion along the -axis direction:where , .(4)The roll motion of the front body and the rear body:(5)The roll motion of the rear axle:(6)Vehicle pitch motion around the -axis: (7)Yaw motion of the rear body and the rear axle:(8)Yaw motion of the front body:where

is the contact force which is generated when the swing bridge reaches the limit position. In order to calculate , the contact process is assumed to be connected by a spring and damping system, as shown in Figure 4; can be written as follows:where is contact stiffness and is contact damping. and when the swing bridge does not reach the limit position. is the distance between the oscillating pin and the limit position, .

3. The Tire Model

The Fiala tire model describes the relationship among the lateral force, the vertical force, the longitudinal force, the slip angle, and the slip ratio in detail, and, furthermore, it has been verified by experiments [15]. Therefore, the Fiala tire model is adopted to calculate the longitudinal force, the vertical force, and the lateral force of the tires in this paper.

3.1. Side-Slip Angle

Lateral wheel slip is defined as (11). The side-slip angles of the tires in our model can be calculated using (12). Hence,

3.2. Normal Force

Consider

where is each displacement of the four tires in the -direction of the global coordinate system, which can be calculated as follows:

where is the vertical displacement of the th tire caused by obstacles.

The ASVs may pass over various shapes of obstacles in actual conditions, but we cannot consider every obstacle in the theoretical analysis. In order to facilitate the establishment of a mathematical model of obstacles, regularly shaped obstacles have been chosen, such as triangular, circular, and square-shaped obstacles.

Figure 5 shows that the wheel loader passes over a triangular obstacle. As for the vehicle symmetry, this paper only analyzes roll stability of the ASVs when the left wheels pass over obstacles. can be written as follows.The right front wheel .The left front wheelThe right rear wheel .The left rear wheelwhere , , and .

In the solution process, indicates that the th wheel lifts off the ground; set , , and .

3.3. Longitudinal Force

Considerwhere is the critical longitudinal slip which is used to determine whether tires are sliding or not:

The friction coefficient, depending on static friction (), sliding friction (), longitudinal slip ratio (), and side-slip angle (), is defined as

and can be calculated as follows:

, the longitudinal slip ratio, is calculated as follows:

The dynamics of each individual wheel are given bywhere is driving/braking torque and is rotational inertia of wheel.

This paper focuses on the effects of active braking on roll instability of the ASVs. The time derivative of the th wheel longitudinal slip in braking process is given by

Substituting (22) into (23) gives

3.4. Lateral Force

Considerwhere is critical lateral slip, which is used to judge whether tires are sliding or not, , and .

4. The Model Solution

The established dynamic model encompasses the driving/braking torque and the steering torque . When the vehicle is in the process of driving and steering, the driving torque and the steering torque constantly change with the outside input. If the exact value of the real-time change is desired, it is necessary to establish the model of power transmission (engine, hydraulic torque converter, gearbox, and transfer case) and the model of hydraulic steering system, which will make the whole analysis very complex. In theoretical analysis, the forward speed and the steering angle of the vehicle can be controlled by a classical Proportional-Integral-Derivative control (PID control).

Figure 6 is the block diagram of the speed control in the process of model solution. The difference between the predetermined speed and the output of the model is rendered by a PID to be the control on the driving/braking torque . By adjusting in real time, the forward speed in the dynamic equation (1) is ensured to amount to the set value .

The steering angle control is the same as the speed control. The difference between the predetermined steering angle and the output of the model is rendered by a PID to be the control on the steering torque . By adjusting in real time, the steering angle in the dynamic equation (8) is guaranteed to track the set value .

With the calculation methods of the driving torque, the steering torque, and the tire forces, the equation can be written as a MATLAB program. Figure 7 presents the block diagram of the dynamic model solution. It shows the whole calculation process of the dynamic model. Before calculation, the model parameters and the solution parameters need to be set first. The input of the model mainly includes the steering angle, the forward speed, and the road parameters.

5. Experimental Verification of the Dynamic Model

This paper studies the roll stability of the ASVs. Conducting experiments on real vehicles is risky. Moreover, real ASVs are difficult to modify. For these reasons, we designed an experiment prototype (scaled ASV) that has the main features of the ASVs to verify the mathematical model.

5.1. Introduction to the Experiment Prototype

The experiment prototype has the following features:(1)The prototype is mainly composed of three parts: the front body, the rear body, and the rear axle. The front body is jointed with the rear body by a steering pin. The swing bridge is linked to the rear axle by an oscillating pin, as shown in Figure 8, and the upper part of the swing bridge has a fixed connection with the rear body.(2)The experimental prototype adopts a DC speed regulating motor to drive the two rear wheels. An electric cylinder installed between the front body and the rear body is used to steer. The articulated steering angle is measured by a displacement sensor (WYL-100). The correspondence relationship between the displacement sensor and the steering angle is shown in Figure 9.(3)The experimental prototype contains the necessary sensors to measure vehicular motions. A MEMS-based vertical gyroscope VG440CA-200 (Crossbow Technology, Inc., Milpitas, CA, USA) is installed at the intersection point of the articulated steering pin axis and the oscillating pin axis to provide real-time vehicle kinematics data (velocities, accelerations, and angles). The angle sensor (XW-TS1130) installed on the rear axle is used to measure the roll angle of this axle. The NI data acquisition card (USB-6211) is taken as the acquisition system. A Panasonic Toughbook Ul CF-VEBU11AU (Panasonic Co., Ltd., Osaka, Japan) is installed at the rear frame to collect data and to control forward speed and articulated angle. Table 1 shows the dimensions and specifications of the scaled ASV. Note that the parameters, such as moment of inertia and center of gravity position, are obtained from the three-dimensional model of the ASVs built in software CATIA (Dassault Systèmes SA, France).

5.2. Experiments and Results Analysis

To confirm the relevant simulation conclusions, we conducted experiments under two experimental operating conditions, namely, (1) turning on flat roads and (2) passing over obstacle during turning. The inputs of the steering angle are shown in Figure 10. The base and the height of the triangle obstacle are 0.3 m and 0.1 m, respectively. The shape and size of the obstacle and the process of passing over obstacle are shown in Figure 11. In addition, the parameters of static friction and sliding friction are set as 0.6 and 0.4, respectively.

5.2.1. Turning on Flat Road

Figure 12 illustrates the motions of the scaled ASV when turning on a flat road. The results indicate good agreement between yaw rate and lateral acceleration of the model and the test data. Yaw rate and lateral acceleration are simultaneously increased with increasing forward speed from 1 m/s to 2 m/s, which is entirely consistent with the actual situation. While some deviations are also evident, as it is observed in Figure 12(c), the lateral acceleration of the vehicle model at time range 2.5 s–6 s is above the measured data, which may partially be attributed to the constant speed in simulation. In the dynamic model, the vehicle advances at a constant speed, and the velocity is maintained at  m/s (2 m/s) when turning. In practice, however, maintaining vehicle velocity at  m/s (2 m/s) is difficult because speed often decreases when turning. Therefore, the yaw rate and the lateral acceleration obtained by the simulation are slightly higher than those obtained by the test during cornering.

Given that the vehicle turned on a flat road at low speed, the pitch and vertical motion were not obvious and were therefore not described in this work.

5.2.2. Passing over Obstacle during Turning

Figure 13 summarizes the results of both the field experiments and the simulations of scaled ASV motion when turning on uneven ground. The results show good agreement between the model and the test data. The vehicle has been in the turning state before coming into contact with the obstacle, such that the initial values of the lateral acceleration and yaw rate are not equal to zero. The ASV tilts to the right, and the roll angular velocity increases rapidly when the left front wheel comes into contact with the obstacle. Meanwhile, the roll angle positively increases, while the pitch angle negatively increases. The left rear wheel begins to come into contact with the obstacle at 3 s, such that the roll angle of the rear axle and the pitch angle positively increase. The lateral acceleration of the vehicle will drastically change when the wheels touch and leave the obstacle. Vehicle forward speed will change when passing over the obstacle, and the yaw rate will vary with speed. The scaled ASV will return to the steady state after passing over the obstacle.

6. Study on Roll Instability Mechanism of the ASVs

The dynamic model is applied to analyze the lateral stability of the ZL50 wheel loader (produced by XCMG Construction Machinery Co., Ltd., Xuzhou, China). Table 1 lists the parameters of the ZL50 wheel loader. To determine the roll instability mechanism of the ASVs, the following four conditions have been studied: (1) quick turn on flat road; (2) passing over obstacles when driving in a straight line; (3) passing over obstacles when turning; and (4) passing over obstacles on a sloping road.

Here, we choose the lateral transfer ratio (LTR) as stability index to indicate the vehicle rollover propensity, which is defined as

The LTR varies from 0 to 1, where 0 refers to equal vertical forces on both sides of the vehicle (zero roll) and 1 refers to either the left or the right vehicle tires losing contact with the ground [11].

6.1. Stability Analysis of Condition 1: Quick Turn on Flat Road

To determine the instability mechanism of the ASVs in condition 1, we analyzed the vehicle motions corresponding to 12 kinds of parameters, which are listed in Table 2. According to the analysis results, vehicle instability will occur (the value of LTR reaches 1, as shown in Figure 14(a)) when the steering angle is 20°, 25°, and 30° and the forward speed exceeds 8 m/s, 7 m/s, and 6.5 m/s, respectively.

Figures 14(b), 14(c), and 14(d) show the vehicle motions corresponding to critical rollover instability. The lateral acceleration reaches 5 m/s2, the roll angular velocity is close to 0.6 rad/s, and the yaw rate varies from 0.6 rad/s to 0.8 rad/s when the critical instability occurs in condition 1. Roll instability is mainly caused by the excessive lateral acceleration or the excessive yaw rate when the wheel loader turns on a flat road. The value of the yaw rate is not constant when the wheel loader loses stability under different conditions. Thus, the yaw rate cannot be used as the stability index in condition 1. We therefore take lateral acceleration as the stability index when the ASVs turn on a flat road. In this analysis, the wheel loader will lose stability when the lateral acceleration reaches 5 m/s2.

6.2. Stability Analysis of Condition 2: Passing over Obstacles When Driving in Straight Line

The shapes of obstacles are triangular and circular, as shown in Figure 15. Each shape has two different sizes, such that obstacles are of four types: (1)  m,  m; (2)  m,  m; (3)  m,  m; and (4)  m,  m.

By analyzing the movement of the wheel loader as it passes over the four obstacles at different speeds, we found that the wheel loader will lose stability (the value of LTR reaches 1) when the forward speed reaches 7 m/s, 6 m/s, 6.4 m/s, and 5 m/s. By comparing the wheel loader motion parameters, such as lateral acceleration, roll angular velocity, and vertical acceleration, the roll angular velocity is found to be close to 3 rad/s under the four critical instability conditions, as shown in Figure 16. Therefore, the roll stability of ASVs can be expressed by the roll angular velocity when the ASVs pass over obstacles in a straight line. In this analysis, the wheel loader is in a stable state when the roll angular velocity is less than 3 rad/s; otherwise, the wheel loader will lose stability.

6.3. Stability Analysis of Condition 3: Passing over Obstacles When Turning

The stability of the wheel loader has been analyzed under two conditions: (1) turning on flat road and (2) passing over obstacles when driving in a straight line. The results showed that the stability under condition 1 mainly depends on the lateral acceleration, whereas that under condition 2 mainly depends on the roll angular velocity. We then analyze the condition when the wheel loader passes over obstacles when turning and study the relationship between the lateral acceleration and the roll angular velocity when the ASVs lose stability. The lateral acceleration values are set as 1 m/s2, 1.5 m/s2, 2 m/s2, 2.5 m/s2, 3 m/s2, 3.5 m/s2, and 4 m/s2 by adjusting the steering angle and the forward speed. The wheel loader loses stability by adjusting the size of the obstacles when the vehicle is running with one of the above lateral acceleration values. We obtain the roll angular velocity corresponding to critical rollover instability in each of the above lateral acceleration values.

Taking the lateral acceleration as the -axis and taking the critical roll angular velocity as the -axis, we can draw the relationship between the lateral acceleration and the roll angular velocity, as shown in Figure 17. When the wheel loader is traveling in a straight line, and the lateral acceleration is equal to 0, instability will occur once the roll angular velocity exceeds 3 rad/s, which is consistent with the result in Section 6.2. The roll angular velocity that caused instability will decrease as lateral acceleration gradually increases. The wheel loader will directly lose stability when the lateral acceleration increases to 5 m/s2, which is consistent with the result in Section 6.1. In addition, the roll angular velocity rapidly decreases when the lateral acceleration exceeds 4 m/s2, such that instability of the vehicle is mainly determined by the lateral acceleration at this time.

6.4. Stability Analysis of Condition 4: Passing over Obstacles on Slopes

When the wheel loader is traveling on a sloping road, as shown in Figure 18, the gravitational acceleration generates subgravitational acceleration in the -direction. This condition indicates that the lateral acceleration of the vehicle is not equal to 0 and the roll angular velocity that caused instability will change when the vehicle travels on a sloping road.

In the above analysis, we need to consider the factor of the slope angle and improve the dynamic model established in Section 2. Here, we analyze the condition in which the wheel loader is traveling in a straight line on a sloping road. The improved dynamic model is shown as follows:(1)Vehicle motion along the -axis direction:(2)Vehicle motion along the -axis direction:(3)Vehicle motion along the -axis direction:where , .(4)The roll motion of the front body and the rear body:(5)The roll motion of the rear axle:(6)Vehicle pitch motion around the -axis:

Slope angles are set as 2.5°, 5°, 7.5°, 10°, 12.5°, 15°, and 17.5°. The shape of the obstacle is triangular, the base  m, and the height  m. The wheel loader loses stability by adjusting the forward speed when the vehicle is traveling on a road with one of the above slope angles. We can obtain the roll angular velocity corresponding to the critical rollover instability in each of the above slope angles.

Taking the slope angle as the -axis and the critical roll angular velocity as the -axis, we can draw the relationship between the slope angle and the roll angular velocity, as shown in Figure 19. The slope angle and the roll angular velocity exhibit an exponential relationship. When the slope angle is equal to 0, instability will occur once the roll angular velocity exceeds 3 rad/s, which is consistent with the result in Section 6.2. A larger slope angle results in a smaller roll angular velocity corresponding to critical instability.

7. Stability Index of the ASVs

In summarizing the previous analysis, the following results can be obtained:(1)Stability mainly depends on the lateral acceleration when the wheel loader is turning on a flat road. The wheel loader is in a stable state when the lateral acceleration is less than 5 m/s2; otherwise, the vehicle will lose stability.(2)Stability mainly depends on the roll angular velocity when the wheel loader is passing over obstacles in a straight line. The wheel loader is in a stable state when the roll angular velocity is less than 3 rad/s; otherwise, the vehicle will lose stability.(3)Stability depends on both the lateral acceleration and the roll angular velocity when the wheel loader is passing over obstacles while turning. The roll angular velocity that leads to the ASVs rollover will decrease with the increase of lateral acceleration.(4)Critical roll angular velocity is affected by the slope angle when the wheel loader is traveling on a sloping road. A larger slope angle results in smaller roll angular velocity corresponding to instability.

According to the above results, we can conclude that the roll stability of ASVs is determined by the roll angular speed, lateral acceleration, and slope angle. The relationship between the roll angular velocity and the lateral acceleration can be expressed as in Figure 17. The relationship between the roll angular velocity and the slope angle can be expressed as in Figure 19. After obtaining the relationship among the roll angular speed, lateral acceleration, and slope angle, the of the ASVs can be uniformly expressed with the roll angular velocity, as defined in (33). indicates that the ASVs are in a stable state. A greater value indicates better ASVs stability. indicates that the ASVs are in a critical state of instability. indicates that the ASVs lose stability. A smaller indicates poorer stability. Hence,where is the initial rollover threshold and  rad/s in this study. is the influence coefficient of the lateral acceleration on the rollover threshold, . is the influence coefficient of the slope angle on the rollover threshold, .

The influence coefficients of and can be calculated from Figures 17 and 19, respectively. We transform Figure 17 into the form shown in Figure 20 and linearly fit the curve in Figure 20. The expression of influence coefficient can thus be obtained:

In the same way, we transform Figure 19 into the form shown in Figure 21 and fit the curve in Figure 21. The expression of impact factor can thus be obtained:

8. Conclusions

In this study, a typical ASV, that is, a wheel loader, is selected as the research object. A dynamic ASV model is established and verified with experiments. On this basis, a detailed analysis of the dynamic rollover stability of a wheel loader is conducted. The main conclusions are as follows:(1)According to the structure and the movement characteristics of the ASVs, a 7-DOF nonlinear dynamic model is established. This model considers articulated steering movement, swing bridge movement, and ground environmental parameters. To verify the reliability of the established model, a physical prototype of ASV is designed and fabricated. Experiments are conducted under two types of working conditions: turning on flat roads and passing over obstacle during turning. The experimental results mostly conform to the simulation results of the dynamic model, which indicates that the established dynamic model can truly reflect the movement characteristics of the ASVs.(2)Detailed analysis of the rollover stability of a typical ASV, that is, a ZL50 wheel loader, is conducted under different conditions. The results show that (a) roll instability is mainly caused by excessive lateral acceleration when the ASV is turning on a flat road; (b) roll instability is mainly caused by the large angular velocity when the ASV passes over obstacles in a straight line; (c) roll stability is determined by both the angular velocity and the lateral acceleration when the ASV passes over obstacles during turning; and (d) the slope angle can decrease the ASV stability.(3)In summarizing the various instability conditions, we can conclude that the roll instability of the ASVs is mainly caused by excessive lateral acceleration, large angular velocity, and oversized slope angle. Considering all the factors that cause instability, the of the ASVs is presented. The indicates the ASVs stability by detecting whether the roll angular velocity exceeds the critical threshold. The considers information on roll angular velocity, lateral acceleration, and slope angle. This index can thus be used as early warning index for the ASVs.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grants no. 51505177 and no. 51505176). The authors are also indebted to all members of the laboratory for helping to do experiments in the entire process.