A Model-Based Method for Estimating the Attitude of Underground Articulated Vehicles

This paper presents a novel model-based method for estimating the attitude of underground articulated vehicles (UAV). We selected the Load–Haul–Dump (LHD) vehicle as our application object, as it is a typical UAV. First, we established the involved models of the LHD vehicle, including a kinematic model, the linear and angular constraints of a center articulation model, and a dynamic four degrees-of-freedom (DOF) yaw model. Second, we designed a Kalman filter (KF) to integrate the kinematic and constraint models with the data from an inertial measurement unit (IMU), overcoming gyroscope drift and disturbances in external acceleration. In addition, we designed another KF to estimate the yaw based on the dynamic yaw model. The accuracy of the estimations was further enhanced by data fusion. Then, the proposed method was validated by a simulation and a field test under different dynamic conditions. The errors in the estimation of roll, pitch, and yaw were 3.8%, 2.4%, and 4.2%, respectively, in the field test. The estimated longitudinal acceleration was used to obtain the velocity of the LHD vehicle; the error was found to be 1.2%. A comparison of these results to those of other methods showed that the proposed method has high precision. The proposed model-based method will greatly benefit the location, navigation, and control of UAVs without any artificial infrastructure in a global positioning system (GPS)-free environment.


Introduction
With the increase in depth of underground mines, a hazardous environment and safety problems are becoming increasingly prominent. Autonomous vehicles (AVs) are urgently needed to ensure that underground mines are safe and efficient [1,2]. As a typical and important kind of underground vehicle, underground articulated vehicles (UAVs) are widely used in underground mines due to their advantages of a higher maneuverability and efficiency [3]. Therefore, automated UAVs (AUAVs) have first priority in underground mines. However, the automation of their precise control and navigation requires related information on the target vehicle [4][5][6], which includes the attitude, velocities, and even accelerations relative to different directions [7].
The acquisition of the attitude and states of vehicles has been studied intensively for many years [8][9][10][11]. What needs to be emphasized is that, although a global positioning system (GPS) can provide the related information with high accuracy, obstruction of the signal can cause the system to be invalid in some urban areas and tunnels, which is an especially serious problem in underground mines [9]. Two kinds of methods were designed in order to cope with this issue with UAVs according to the specific characteristics of underground mines: infrastructure-based methods [8] and sensor-based methods [12,13]. In addition, a mixture of both is sometimes used [14]. As efficient vehicles for the special structure of articulated vehicles. The states of the target vehicle were obtained by data fusion. A field test showed that the errors in the pitch, roll, and yaw angles in a static state were 0.027 • , 0.025 • , and 0.426 • , respectively. The pitch and roll motions were found to vary between ±20 • , and the yaw motion was found to vary between 0 • and 60 • . The corresponding errors were all less than 4%. However, the results were obtained on pavement, and the system did not cope well with gyroscope drift. Drift in the gyroscope and interference with the accelerometer remain key problems with IMUs, and some methods have been proposed to improve the accuracy of estimations. These include the threshold-based switching approach [25], the adaptive filter method [26], a model-based algorithm [10], and a comprehensive algorithm [27]. Furthermore, Ahmed et al. [28] investigated a novel algorithm for attitude estimation using low-cost IMU sensors, which was applied to a land vehicle. A kinematic model and the vector norm property of gravity were used to estimate the pitch and roll angles. After that, the obtained results were used to estimate the yaw angle in combination with a magnetometer. The results of a simulation and field tests showed that the proposed algorithm was robust and accurate under different dynamic conditions.
Most of the target vehicles that were studied in previous papers were four-wheeled commercial vehicles [29], which can be considered to be single rigid bodies. However, different from the target vehicles in the abovementioned studies, UAVs contain two parts, which are connected to each other with center articulation [24]. Although some methods have been applied to articulated vehicles, the dynamics, the kinematic model, and the defects in IMU sensors have not been considered. Moreover, the roll and pitch angles of UAVs have been neglected in related research [9,22], as the motion of the target vehicle was always considered to be on a plane. The motions of UAVs involve violent pitch and roll motions that occur during normal operation between different levels of underground mines or on an uneven road. Estimations of roll and pitch are essential to the navigation and dynamical control of UAVs in an underground mine environment.
In addition, as the most important piece of information about a target vehicle, the yaw motion is measured by a magnetometer in an IMU. The magnetometer measures the magnetic force between different yaw motions to obtain the yaw angle. Unfortunately, UAVs are usually put to work in an underground metal mine, where the magnetic field may be disordered by metallic minerals. In addition, in the commonly used method, the accuracy of the yaw angle's estimation cannot be improved by using data on acceleration. Therefore, the yaw angle of a UAV cannot be obtained by an IMU alone. Inspired by the work in [9], we use the information from measurements and a kinematic model of center articulation to improve the accuracy of yaw angle estimation, which is obtained primarily by a dynamic model. Motivated by the above discussion, we make a slight improvement in the accuracy of the estimation of the state of a UAV. The three original contributions that distinguish the present work from previous works can be summarized as follows: (1) A novel model-based method for estimating the attitudes of UAVs is proposed. In the first step, the kinematic model and the constraint of center articulation models are used to overcome the drift in the IMU.
(2) A model-based algorithm that combines the dynamic model with a KF is developed to estimate the yaw motion of an LHD vehicle. (3) A fusion of the data from different models and sensors is carried out to improve the accuracy of attitude estimations.
The remainder of this paper is organized as follows. In Section 2, the involved models are established, these models are embedded into the KFs, and algorithms are described for the pitch and roll motion and yaw motion. In Section 3, a simulation of an LHD vehicle is carried out under different dynamic conditions, and data from the simulation are used to verify the method. Section 4 presents the results of a field test and a discussion. Finally, Section 5 summarizes our conclusions and future work.

Method
The proposed method is briefly described in the diagram shown in Figure 2. Some sensors are necessary, including IMUs, an encoder, a pressure sensor, and a length sensor. In Figure 2, the green line denotes the attitude of the roll and the pitch. KFs are the major algorithms. The yaw and other attitudes are obtained separately. In Figure 2, the black line after the angular model denotes the yaw of the LHD vehicle, and the green line denotes the other attitudes. First, the data from the IMU that is located on the front frame are used to obtain the dynamic acceleration and attitude (roll and pitch) of the front frame. During this process, the kinematic model of the front frame and the norm characteristics of the gravitational acceleration vector are used by the KF. Second, the information on acceleration is combined with linear constraints to calculate the acceleration of the rear frame, which can be used to calculate the attitude of the pitch and roll of the rear frame. Simultaneously, the attitude of the front frame is used to calculate the attitude of the rear frame with angular constraints. Additionally, an improvement is carried out with data fusion to obtain the final attitude of the LHD vehicle. The dynamic model of yaw is used to obtain the yaw of the front frame with the help of the KF and data from the sensors. The corresponding part of the rear frame is obtained by the model of angular constraints.

Kinematic Model of the LHD Vehicle
The directions of the linear and angular velocities of the LHD vehicle are defined as shown in Figure 3. The front and rear frames are considered to be rigid bodies. According to the motion of a rigid body, the kinematic model of the front frame and the rear frame of the LHD vehicle can be expressed as: where the subscripts x, y, and z indicate the different directions in Figure 3. The subscripts 1 and 2 indicate the front frame and the rear frame of the LHD vehicle, respectively. The kinematic constraints were obtained between the two frames, which contain linear and angular constraints. More specifically, the linear constraints in different directions can be expressed as: where ϕ is the steering angle of center articulation, and L f2 and L r2 are the distance from the center of gravity of the front frame and the rear frame to the point of articulation, respectively, as shown in Figure 3. Similarly, the angular constraints in different directions can be expressed as: where θ, ϕ, and Ψ are the pitch, roll, and yaw angle, respectively. According to [10], the external acceleration must be removed to improve the accuracy of pitch and roll estimation. The constraints (2a)-(2c) may be useful for the removal of external acceleration. In addition, Equations (3a)-(3c) can also be used to estimate the attitude of the rear frame. However, neither of them can be the only source of the attitude's estimation because of the overly ideal assumptions and the errors in the measurement.

Dynamic Model of Yaw for the LHD Vehicle
The dynamic model of yaw is used to restrain the drift in the yaw information from the gyroscope in the IMU. Considering the efficiency of the method, a simplified dynamic model of the LHD vehicle was adopted, which is shown in Figure 4. The directions of the motions are the same as the directions that are defined in Figure 3. According to Newton's second law, the simplified dynamic model can be given by: where I z1 and I z2 are the moment of inertia of the front frame and the rear frame, respectively. F xf , F xr , F yf , and F yr are the force of front and rear tire in X and Y directions, respectively. T 0 is the steering moment, and although it can be modeled as an equation of φ in a steady state [9], it normally operates in an unsteady state. Here, we consider the steering moment as the input of the dynamic model, and it can be indirectly measured by the pressure and displacement sensors [30]. It is sufficient to focus on one of the two yaw motions in Equation (4) considering the constraints in Equation (3c). In addition, the cornering force of a tire can be given by: where k i and α i are the side-slip stiffness and the angle of the tire, respectively. They can be described by: Therefore, the state equation of the dynamic model of yaw for the front frame can be given by: .
where the sideslip angle β can be expressed as: In addition to the kinematic model and the dynamic model, a model of the sensors should also be taken into account.

Model of the Sensors
Two kinds of sensor are applied in the proposed method, i.e., an IMU sensor and an angular encoder. Considering the noise in the signals from the involved sensors, which include a tri-axis gyroscope, a tri-axis accelerometer, and an angular encoder, the models of the involved sensors can be described in the corresponding frame as follows: where the subscript i = 1,2 denotes the front frame and the rear frame of the LHD vehicle, respectively, y i,G is a column vector [ω i,Gx ω i,Gy ω i,Gz ] T , which contains three angular velocities of the X, Y, and Z axes in the frame of installation, vector ω i,G is the corresponding true value [ω i,x ω i,y ω i,z ] T , and n i,G is white Gaussian noise with a zero mean about each axis. In Formula (9b), y i,A is also a column vector [g i,Ax g i,Ay g i,Az ] T , which contains the acceleration of gravity g i,A = [g i,x g i,y g i,z ] T , the true acceleration of the dynamics of the target a A = [a i,Ax a i,Ay a i,Az ] T , and the corresponding white Gaussian noise with a zero mean n i,A . Similarly, in Formula (9c), y v is a scalar value of the longitudinal velocity from the angular encoder including the true value v x and noise n v . In Formula (9d), y c is also a scalar value of the angular transducer, which contains the true angle φ c and noise n c .

State Vector of Pitch and Roll
Considering the conditions of an LHD vehicle, the Euler angle method is a good choice for attitude representation. Then, the states of the LHD vehicle can be defined as pitch (θ i ), roll (φ i ), and yaw (ψ i ). The transformation relationship between the vehicle coordinate frame (VX) and the navigation coordinate frame (NX) can be expressed as [27]: where T i is the transformation matrix, and the subscripts i = 1,2 represent the front frame and the rear frame, respectively. The transformation matrix can be expressed as: where c and s are the trigonometric functions cosine and sine, respectively. We note that the last row of the transformation matrix only contains the pitch and roll angles, which can be used to calculate these angles with: where T ij k is the (j, k)-th element of T i .
As mentioned in [10], the state vector of a target vehicle at time t, which is also actual state, can be expressed as: The acceleration of gravity in Equation (9b) can be given by: where g is the gravitational acceleration constant, which is 9.81 m/s 2 .

State Vector of Yaw
According to Equation (7), the yaw motion of the front frame can be estimated by a state equation, and the corresponding angle of the rear frame can be obtained by Equation (3c). Therefore, the state vector of yaw motion can be defined as:

KF for Pitch and Roll
The model of a linear KF, which includes the model of process and measurement, can be expressed as: where X i,t is the state vector at time t as defined in Equation (13), Φ t−1 is the state transition matrix, H i,t is the observation matrix, and w and v are the white Gaussian noise in the process and measurement, respectively. The transformation matrix can be calculated with an integration from t − 1 to t [13]: where I 3 is a 3 × 3 identity matrix, ∆t is the interval of the sampling, and the symbol-denotes a cross-product of the target vector [27]. Not all of the terms in T i,t are involved, as only the last row of T i,t is useful for the estimation of attitude. We replace the transformation matrix T i,t in Equation (17) with the state vector X i,t , which is defined in Equation (13). The new integration can be given by where ω can only be measured by the gyroscope in the IMU, which is modeled in Equation (9a). We substitute Equation (9a) into Equation (18) to obtain Comparing Equation (19) with Equation (16a), the state transition matrix and the noise in the process can be expressed as: When combined with Equation (15), the covariance of the noise in the process can be expressed as: where E is the expectation operator, and the covariance of the noise from the gyroscope [27] is equal to σ i,G 2 I 3 , which is based on the covariance being the same in different directions.
In order to improve the accuracy of estimations, we identify the external acceleration with the kinematic model of the target vehicle [10]. Equations (1a)-(1c) can be used to determine the external acceleration of the front frame while neglecting some minor terms. Finally, the kinematic model can be given by: where the longitudinal velocity v x1 is vital to the model, and the angular velocity ω x1 and ω z1 can be replaced by data from the IMU. In addition, . v y1 and . v z1 can be described with a first-order low-pass filter [10]: where c 1y,a and c 1z,a are constant between 0 and 1 and determine the cut-off frequency [10], which is considered to be equal to c a , and ε t is the error in this model, which is time-varying. We substitute Equations (22a) and (22b) into Equations (23a) and (23b), respectively, to obtain: where a t,1 = [a 1y,t a 1z,t ] T , y Gt,1 = [G 1z,t G 1x,t ] T , n G,1 = [n Gz n Gx ] T , n v,1 = [n v n v ] T , and ε t = [ε yt ε zt ] T .
Then, we substitute Equation (5b) into Equation (24) to obtain: where y At,1 = [A y,t A z,t ] T , H m is a matrix for coping with our disregard of the acceleration along the X-axis, and n A,1 = [n Ay n Az ] T . In order to correspond with the measurement model in Equation (16b), Equation (25) can be written as: y At,1 − c a a t−1 − y Gt,1 v t = H m gX t + n G,1 n v,1 + ε t + n A,1 .
Then, we can obtain Then, we can obtain the improved estimate of X 1,t (2) and X 1,t (3) in the state vector X 1,t . All the terms of the vector satisfy: This provides an approach to calculating X 1,t (1) in combination with the measurement model, which can be given by [27] y Ax,t = 1 − X 1,t (2) 2 + X 1,t (3) 2 + γ t .
Then, we can obtain z t,2 = y Ax,t (30a) After the external acceleration of the front frame has been obtained, on the one hand, the data regarding acceleration from the IMU can be used to estimate the attitude of the front frame; on the other hand, these data can also be used to estimate the corresponding acceleration of the rear frame by the derived linear kinematic constraints. The attitude of the front frame can also be used to estimate the attitude of the rear frame with the angular constraints.
According to Equations (3a)-(3c), the acceleration of the rear frame can be expressed as where y At,2 = [A x,t A y,t A z,t ] T and o t (a t,1 ) = [A x,t A y,t A z,t ] T . Then, we can obtain z t,3 = y At,2 − o t (a t,1 ) (32a) Similarly, the pitch and roll of the rear frame can be obtained with Equations (12a) and (12b). According to Equations (3a)-(c), the attitude vector of the rear frame can also be expressed as a function of the front frame: In a general way, we can obtain the attitude of the rear frame by substituting the estimation of the attitude of the front frame into Equation (33).

KF for Yaw
According to Equations (7), (12a) and (12b), the KF for yaw can be expressed as , and u z,k = T 0,k β k T , H z = 1. Then, the yaw of the front frame can be obtained by: The corresponding part of the rear frame can be obtained with Equation (3c).

Fusion of States
The attitude of the rear frame, which can be obtained by the linear and angular constraints, should be fused to improve the accuracy of estimations. The attitudes from different sources can be expressed as: where X t and X t are the estimated value from the angular and linear constraints, respectively.
The fusion of the data can be given bŷ where U −1 1,t and U −1 2,t are the corresponding estimation errors, respectively.

Simulation Verification
To verify the validity and accuracy of the proposed method, a simulation was performed with a verified prototype model of the LHD vehicle in ADAMS, and the relevant parameters were shown in Appendix A. The data were processed with the proposed method to estimate the attitude of the LHD vehicle.

Simulation Setup
According to [15,25], the sampling frequency of the simulation was set to 100 Hz. The speed of the LHD vehicle was set to 14 km/h under normal running conditions and 9.5 km/h during the constant radius steering process. The noise variances of the accelerometer and the gyroscope were set to 0.4 mg/
As we can see from Figure 5, during C1, the LHD vehicle was static, the accelerations in the lateral and longitudinal directions were zero, and the accelerations of two of the frames in the vertical direction were the same as gravity. During C2, the LHD vehicle accelerated in a longitudinal direction according to the acceleration data in Figure 5a, and the corresponding angular velocity was very small because the LHD vehicle has no suspension. During C3, the left and right front wheels of the LHD vehicle simultaneously encountered a bump; as a result, the vertical accelerations in vertical motion and pitch motion changed dramatically. During C4, the LHD vehicle ran over a sine bump, and the pitch motion changed with the rise and fall of the road. After that, the LHD vehicle ran at a constant speed, C5, which was similar to the speed during C1.
According to Figure 6, the velocity of the LHD vehicle oscillated during C4 and C5. This was caused by the rough control strategy for velocity during the simulation. However, the range remained within the LHD vehicle's normal running velocity range.
Then, the left wheel of the LHD vehicle encountered a sine bump, C6, which brought about changes in the roll motion, and the lateral acceleration also changed with the rise and fall of the road. During C7, the yaw motion of the front frame and the rear frame changed obviously. During C8, the acceleration in the longitudinal direction, in contrast to C2, was negative, and the speed of the LHD vehicle was also reduced to 9.5 km/h. Then, the LHD vehicle steered in a constant radius. At the start of C9, the angular velocities of the front frame and the rear frame were different, which implies that the angles of central articulation had changed.

Estimation Results
The results of the estimation of the pitch, roll, and yaw angles of the LHD vehicle are shown in Figures 7 and 8, respectively, and the time history of the error covariance matrix during the estimation was shown in Figure 9. The upper part in the figures is a comparison between the estimated values and the true values, and the lower part shows the errors in the estimation. We found that the roll and pitch of the front frame and the rear frame were almost the same, so only the attitude of the rear frame is shown in Figure 7; however, we included the yaw of the two frames. As we can see from Figures 7-9, the estimated values and the true values are in good agreement. The RMS of errors in the roll, pitch, and yaw were 0.2 • (2.4%), 0.19 • (2.1%), and (1.1 • ) 0.23%, respectively. The percentage was obtained with the maximum RMS and the largest values of corresponding angles. The maximum error in the estimations appeared in C9 (the constant radius steering condition), which was caused by the slip of the LHD vehicle during steering. Although small differences appeared during the dynamic process according to the error, the estimations were found to be satisfactory.   Then, the estimated longitudinal accelerations of the rear frame were integrated to obtain the velocity of the LHD vehicle. Figure 10a compares the estimated velocity with the reference velocity. Figure 10b shows the error in the estimation of the velocity. During 15-18 s condition, the drastic changes of the pitch angle and characteristic of the tire would increase the estimation error of acceleration and the dynamic model, and the error of acceleration caused the deviation of velocity. Therefore, a deviation occurred. The error in the estimated speed was less than 0.27 km/h (3%). This means that the proposed method can be used in the control and navigation of an LHD vehicle without any infrastructure.

Experimental Setup
The field test of the LHD vehicle was carried out on an uneven road that was similar to the road in underground mines. The test system is shown in Figure 11. The LMS SCADAS ( 1 , SCM05) was used to record the data during the field test, and the sample frequency was set to 100 Hz according to [15]. Two IMUs ( 2 , 3 , Beijingsanchi, SC-AHRS-200A) were installed in the front frame and the rear frame, respectively. Two kinds of encoders were used to measure the articulated angle of a frame ( 4 , MIRAN, WOA-C) and the speed of transmission ( 5 , SCHMEASAL IFL 5-18M-10P), respectively. A gyroscope ( 6 ) with high precision (the bias stability is less than 0.05 • /h, and the random walk coefficient is less than 0.005 • / √ Hz) was installed to provide references. Different conditions were adopted to simulate the working conditions of the LHD vehicle. Usually, the characteristics of the route of an LHD vehicle are special in underground mines, and always include a square turn and a change in operation direction. Thus, the route for the field test was planned as shown in Figure 12, where the direction of the arrow represents the direction of the LHD vehicle and the solid and hollow lines correspond to the forward and backward directions, respectively. The key dimensions are L1 = 50 m and L2 = 20 m. The operations between Point A and Point B were used to simulate the operations between the blast point and the dump point of the underground mines [2]. Although the trajectory in the field test was not so long as the other autonomous ground vehicle, the operation and dynamic characteristics in this route can cover all of characteristics of the LHD during the normal operations.

Experimental Results
The results of the estimation of the front and rear frames were almost the same according to the data from the field test; thus, only the results of the estimation of the rear frame's attitude angle are shown and were analyzed. The results of the estimation corresponding to the roll, pitch, and yaw are shown in Figures 13-15, respectively. The upper part (a) in each figure is a comparison between the reference values and the estimation, and the lower part (b) shows the error in the estimation, and the dramatic or typical changes of errors were detailed. The reference values were sourced from the data of a high-precision reference gyroscope.

Discussion of Results
As we can see from Figures 13 and 14, the estimations of roll and pitch are consistent with the reference values, while some errors fluctuate sharply. This is due to the special structure of the LHD vehicle's chassis, which is also a common characteristic of UAVs. The frames of an LHD vehicle are rigidly connected to the front and rear drive axle, which means that these vehicles do not have suspension that can filter vibrations from the road. Although rubber pads are used to eliminate sharp fluctuations, the improvement is limited. Hence, vibrations bring about sharp fluctuations in the errors in the estimation.
As we can see from Figure 15, the estimation of yaw is highly consistent with the reference value, and the error is less than 5 degrees and does not increase with time. The error only obviously changed after the square turn (100-150 s, 260-310 s). During the square turn, the tire of the LHD vehicle seriously slipped, so error would be brought about by the change in Equations (8a) and (8b). In addition, according to Figure 16, ripples also occurred between 164 s and 166 s. This error could be traced back to the input of the estimation model, which is the steering moment. The steering moment is decided by the pressure in the cylinder. The pressure is shown in Figure 16, where "right" and "left" refer to the rod-less cavity pressure in the right and left steering cylinder, respectively. As we can see from the Figure 16, the two pressures rippled alternately-between 12.5 MPa and 0.2 MPa-after a steering maneuver. This means that the oscillatory yaw motion in the frames occurred, which is always caused by instability in the hydraulic steering system. Although different errors occurred during the estimation, the accuracy was found to meet the requirements of practical applications. During the estimation, the estimated value of the longitudinal dynamic acceleration is converted to the velocity of the LHD vehicle by integration. Figure 17 compares the estimated velocity obtained from the shaft encoder with the reference value. As we can see from the figure, not only are the trends consistent, but the estimated and reference values are also very close according to the details in Figure 17. Despite all this, the accuracy of the estimation is high enough to indicate the attitude and speed of articulated vehicles, which can be verified by the root mean square (RMS) of the error in the estimated values, which are listed in Table 1. The RMS of errors for the roll, pitch, and yaw were 0.19 • (3.8%), 0.1 • (2.4%), and 2.08 • (4.2%), respectively. The error in the estimated velocity was 0.21 km/h (1.2%). Compared with the estimation error of 5% in [9], the proposed model-based method improves the accuracy of estimations and involves more information on attitude, which is useful for the navigation and control of autonomous vehicles [28].

Conclusions
A model-based method was presented for estimating the attitude of articulated vehicles in a GPS-free environment. A kinematic model, a two-constraint model of articulation, and a dynamic four degrees-of-freedom (DOF) yaw model of UAVs were established and integrated with rough data from IMUs by KFs to obtain information on attitude. Then, the method was verified by a simulation under different dynamic conditions. A field test of the LHD vehicle was carried out to prove the method's practical efficacy.
The results of the field test showed that the method could estimate the yaw and velocity accurately. The errors in the estimated roll, pitch, and yaw were 3.8%, 2.4%, and 4.2%, respectively. The error in the estimated velocity was 1.2%. The characteristics of these errors were analyzed and interpreted. The errors in the estimated pitch and roll were brought about by the direct vibrations from the ground. The errors in the estimated yaw included trend errors and ripples. The trend errors occurred at the square turn, when the constraints between the tires and the road changed. Local ripples were induced by instability in the hydraulic steering system.
A trajectory of the typical work cycle for LHD was designed during the validation field test. The field test can prove the validation of this method during the work cycle, which is most of working conditions of the LHD vehicle. For the rarely long-traveled transition, this method still needs further verification. Therefore, motivated by these limitations and sources of error, future research will focus on improvements in accuracy. The method would be verified during the long-traveled transition, the dynamic yaw model will be improved by covering tire slips, and instability in the hydraulic steering system will be studied and controlled to eliminate ripples. In addition, the algorithm will also be improved to decrease the effect that vibrations have on the frame.

Acknowledgments:
The authors would like to acknowledge Shandong Gold Group Co., Ltd. for the experimental articulated vehicle and test sites.

Conflicts of Interest:
The authors declare no conflicts of interest.