Research on the Modeling, Control, and Calibration Technology of a Tracked Vehicle Load Simulation Test Bench

Featured Application: The device mentioned in this paper can be used for the road simulation test of a power system of a heavy haul vehicle. The method mentioned in this paper can be used to model, control, and calibrate the vehicle load simulator. Abstract: The load simulation test bench plays an important role in tracked vehicle development. The stability and accuracy of the system have a vital impact on the experimental results. To accurately reproduce the power performance of a tracked vehicle on the test platform, this paper aims to investigate the model, control, and calibration method of the test bench. Firstly, the dynamic model of a tracked vehicle under complex driving conditions is analyzed and established, which takes driving torque as the input and driving wheel speed as the output. Then, considering the uncertainties and disturbances in the system model, a 2-degree-of-freedom (2-DOF) control method combined with a disturbance observer is proposed to solve the stability problem of the system. Furthermore, in order to investigate the accuracy of the simulation on the test bed, a method of calibrating the system by a ﬂywheel set with standard inertia is proposed. In the calibration process, the inﬂuence of the system resistance torque and the original mechanical inertia on the results is considered, and the response time of the inertia simulation is analyzed in both a steady and dynamic state. Finally, the load simulation test is carried out with the corrected system. The test results show that the system has a high load simulation accuracy under various load simulation tests.


Introduction
Tracked vehicles are widely used in agriculture, weapons, and engineering construction because of their superior performance under complex conditions. The load simulation technology of a tracked vehicle is rapidly evolving because it can provide a convenient, effective, and lifelike test environment [1][2][3]. By modifying different parameters on the test bench, the power output, emission performance, and braking performance of the vehicle under different temperature, humidity, and ground conditions can be tested. Additionally, only when the body inertia and ground load simulated by the loading system are consistent with the real situation can the load simulation research test be meaningful [4][5][6][7]. At present, in order to overcome the problem that the traditional mechanical simulation method has a narrow simulation range and cannot achieve stepless simulation, most 1 -electric motor, 2 -reducer, 3 -torque meter and encoder, 4 -tested vehicle, 5 -torque meter and encoder, 6reducer, 7 -electric motor.

Dynamic Model of a Tracked Vehicle
In order to achieve accurate load simulation, firstly, the dynamics model of the tracked vehicle needs to be analyzed to obtain the relationship between the torque and rotation speed of the driving wheel.
Tracked vehicles are subjected to the interaction of driving force, inertia force, and driving resistance in the course of driving. The force exerted on the body reaches the moment balance at the active wheel [Error! Reference source not found.]. By identifying the equivalent of the force to the driving wheel, the torque formula of the driving wheel can be obtained as follows: = 21.45 (4) where, is the angular acceleration of the vehicle output axle, is the vehicle inclination, is the wind resistance coefficient, is the vehicle forward projection area, is the vehicle driving moment, is the vehicle inertia moment, is the slope resistance moment, is the wind resistance moment, is the ground resistance moment, T is the road resistance moment, and R is the effective radius of the driving wheel.
Because tracked vehicles often travel on complex roads, it is necessary to study their driving conditions when switching between different slopes in body dynamics modeling. In this paper, the driving conditions of tracked vehicles are divided into six processes, namely, the horizontal slope, increase of slope, downhill slope, decrease of slope, increase of slope, and decrease of slope, respectively, corresponding to stages 1 to 6 in

Dynamic Model of a Tracked Vehicle
In order to achieve accurate load simulation, firstly, the dynamics model of the tracked vehicle needs to be analyzed to obtain the relationship between the torque and rotation speed of the driving wheel.
Tracked vehicles are subjected to the interaction of driving force, inertia force, and driving resistance in the course of driving. The force exerted on the body reaches the moment balance at the active wheel [31]. By identifying the equivalent of the force to the driving wheel, the torque formula of the driving wheel can be obtained as follows: T = T f + T s (5) where, α d is the angular acceleration of the vehicle output axle, γ is the vehicle inclination, C D is the wind resistance coefficient, A is the vehicle forward projection area, T e is the vehicle driving moment, T a is the vehicle inertia moment, T s is the slope resistance moment, T w is the wind resistance moment, T f is the ground resistance moment, T is the road resistance moment, and R is the effective radius of the driving wheel. Because tracked vehicles often travel on complex roads, it is necessary to study their driving conditions when switching between different slopes in body dynamics modeling. In this paper, the driving conditions of tracked vehicles are divided into six processes, namely, the horizontal slope, increase of slope, downhill slope, decrease of slope, increase of slope, and decrease of slope, respectively, corresponding to stages 1 to 6 in Figure 2.

Dynamic Model of a Tracked Vehicle
In order to achieve accurate load simulation, firstly, the dynamics model of the tracked vehicle needs to be analyzed to obtain the relationship between the torque and rotation speed of the driving wheel.
Tracked vehicles are subjected to the interaction of driving force, inertia force, and driving resistance in the course of driving. The force exerted on the body reaches the moment balance at the active wheel [Error! Reference source not found.]. By identifying the equivalent of the force to the driving wheel, the torque formula of the driving wheel can be obtained as follows: 21.45 (4) where, is the angular acceleration of the vehicle output axle, is the vehicle inclination, is the wind resistance coefficient, is the vehicle forward projection area, is the vehicle driving moment, is the vehicle inertia moment, is the slope resistance moment, is the wind resistance moment, is the ground resistance moment, T is the road resistance moment, and R is the effective radius of the driving wheel.
Because tracked vehicles often travel on complex roads, it is necessary to study their driving conditions when switching between different slopes in body dynamics modeling. In this paper, the driving conditions of tracked vehicles are divided into six processes, namely, the horizontal slope, increase of slope, downhill slope, decrease of slope, increase of slope, and decrease of slope, respectively, corresponding to stages 1 to 6 in  Detailed slope decomposition for each stage is shown in Figure 3. In this figure, s means the up (down) slope distance, θ i (i=1,2,3,4) means the slope angle, h c means the centroid height, and l d means the relative horizontal distance between the centroid and geometric center.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 5 of 26 Detailed slope decomposition for each stage is shown in Figure 3. In this figure, s means the up (down) slope distance, (i=1,2,3,4) means the slope angle, ℎ means the centroid height, and means the relative horizontal distance between the centroid and geometric center.
(1) horizontal slope (2) increase of slope(climbing) According to the above figure, the inclination angle and road resistance of the tracked vehicle at a certain instantaneous stage can be obtained. The judgment process is shown in Figure 4. According to the above figure, the inclination angle and road resistance of the tracked vehicle at a certain instantaneous stage can be obtained. The judgment process is shown in Figure 4.
By Laplace transformation of Formula (1) to Formula (5), the vehicle dynamic equation is obtained as follows: where ω em (s) means the speed of the driving wheel and J e means the vehicle equivalent inertia.  By Laplace transformation of Formula (1) to Formula (5), the vehicle dynamic equation is obtained as follows: where ( ) means the speed of the driving wheel and means the vehicle equivalent inertia.

Equivalent Inertia Modeling of Tracked Vehicles
Formula (6) shows that in order to obtain ( ), the inertia of the whole vehicle must be equivalent to the driving wheel. Equivalent analysis of the inertia of the tracked vehicle shows that where, means the vehicle mass (excluding track), means the vehicle speed, means the track mass, means the track speed, means the number of load wheels, means the load wheel rotation inertia around the axle center, means the load wheel angular speed, means the drive wheel rotation inertia around the axle center, means the drive wheel angular speed, means the induced wheel rotation around the axle center Inertia, means the angular velocity of the inducer, means the number of supporting roller, means the inertia of the supporting roller rotating around the axle center, and means the angular velocity of the supporting roller. Equation 7 shows that the translational inertia and rotational inertia of all parts of the vehicle can be expressed by the inertia revolving around the driving wheel. In order to calculate the equivalent inertia accurately, it should be divided into three parts, as shown in Formula 8.
Here, means the vehicle body equivalent inertia (excluding the track and wheel rotational inertia); means the equivalent inertia of the track; and means other equivalent inertia, including that of the driving wheel, induced wheel, and supporting roller. These three kinds of inertia are analyzed separately below. Figure 5 shows the relationship between body speed and driving wheel speed.

Equivalent Inertia Modeling of Tracked Vehicles
Formula (6) shows that in order to obtain ω em (s), the inertia of the whole vehicle must be equivalent to the driving wheel. Equivalent analysis of the inertia of the tracked vehicle shows that where, m v means the vehicle mass (excluding track), v v means the vehicle speed, m t means the track mass, v t means the track speed, n l means the number of load wheels, J l means the load wheel rotation inertia around the axle center, ω l means the load wheel angular speed, J d means the drive wheel rotation inertia around the axle center, ω d means the drive wheel angular speed, J i means the induced wheel rotation around the axle center Inertia, ω i means the angular velocity of the inducer, n b means the number of supporting roller, J b means the inertia of the supporting roller rotating around the axle center, and ω b means the angular velocity of the supporting roller. Equation (7) shows that the translational inertia and rotational inertia of all parts of the vehicle can be expressed by the inertia revolving around the driving wheel. In order to calculate the equivalent inertia accurately, it should be divided into three parts, as shown in Formula (8).
Here, J v means the vehicle body equivalent inertia (excluding the track and wheel rotational inertia); J t means the equivalent inertia of the track; and J o means other equivalent inertia, including that of the driving wheel, induced wheel, and supporting roller. These three kinds of inertia are analyzed separately below. Figure 5 shows the relationship between body speed and driving wheel speed.  Figure 5. The relationship between body speed and driving wheel speed.
According to Figure 5, we can get where δ means the slip rate between the track and ground. So for the body, there is By combining the two forms, we can get When δ = 0, the original format becomes The biggest difference between a tracked vehicle and traditional vehicle is that a tracked vehicle is equipped with high-quality tracks [Error! Reference source not found.]. Because the track is irregular and the absolute speed of each part when driving is different, it should be divided into four parts, including the upper, lower, front, and rear parts, when calculating. Assuming that the tracked vehicle is moving at a uniform speed on the ground with no relative slip, that is, the slip rate δ = 0, the inertia of the grounding part of the track can be ignored. However, in a normal driving process, due to bad ground conditions (such as marsh, snow, etc.), the slip rate is usually not stable, so it is necessary to consider the total inertia of the track. To analyze it, the speed of each part is decomposed into relative speed and frame speed, which is shown in Figure 6. By understanding three velocities as three vectors, we can get According to Figure 5, we can get where δ means the slip rate between the track and ground. So for the body, there is By combining the two forms, we can get When δ = 0, the original format becomes The biggest difference between a tracked vehicle and traditional vehicle is that a tracked vehicle is equipped with high-quality tracks [32]. Because the track is irregular and the absolute speed of each part when driving is different, it should be divided into four parts, including the upper, lower, front, and rear parts, when calculating. Assuming that the tracked vehicle is moving at a uniform speed on the ground with no relative slip, that is, the slip rate δ = 0, the inertia of the grounding part of the track can be ignored. However, in a normal driving process, due to bad ground conditions (such as marsh, snow, etc.), the slip rate is usually not stable, so it is necessary to consider the total inertia of the track. To analyze it, the speed of each part is decomposed into relative speed and frame speed, which is shown in Figure 6. According to Figure 5, we can get where δ means the slip rate between the track and ground. So for the body, there is By combining the two forms, we can get When δ = 0, the original format becomes The biggest difference between a tracked vehicle and traditional vehicle is that a tracked vehicle is equipped with high-quality tracks [Error! Reference source not found.]. Because the track is irregular and the absolute speed of each part when driving is different, it should be divided into four parts, including the upper, lower, front, and rear parts, when calculating. Assuming that the tracked vehicle is moving at a uniform speed on the ground with no relative slip, that is, the slip rate δ = 0, the inertia of the grounding part of the track can be ignored. However, in a normal driving process, due to bad ground conditions (such as marsh, snow, etc.), the slip rate is usually not stable, so it is necessary to consider the total inertia of the track. To analyze it, the speed of each part is decomposed into relative speed and frame speed, which is shown in Figure 6.  By understanding three velocities as three vectors, we can get The upper part of the track is shown in Figure 7. The absolute speed of the upper part v tup consists of the traction speed v tupe and relative speed v tupr .
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 26 R Figure 5. The relationship between body speed and driving wheel speed.
According to Figure 5, we can get where δ means the slip rate between the track and ground. So for the body, there is By combining the two forms, we can get When δ = 0, the original format becomes The biggest difference between a tracked vehicle and traditional vehicle is that a tracked vehicle is equipped with high-quality tracks [Error! Reference source not found.]. Because the track is irregular and the absolute speed of each part when driving is different, it should be divided into four parts, including the upper, lower, front, and rear parts, when calculating. Assuming that the tracked vehicle is moving at a uniform speed on the ground with no relative slip, that is, the slip rate δ = 0, the inertia of the grounding part of the track can be ignored. However, in a normal driving process, due to bad ground conditions (such as marsh, snow, etc.), the slip rate is usually not stable, so it is necessary to consider the total inertia of the track. To analyze it, the speed of each part is decomposed into relative speed and frame speed, which is shown in Figure 6.  By understanding three velocities as three vectors, we can get By understanding three velocities as three vectors, we can get v tup = v tupe + v tupr (13) According to the Vector Superposition Rule, where π is the vector angle between v tupe and v tupr .
In addition, according to the law of the conservation of energy, By simultaneously considering formulae (13) to (17), we can get For the front part of the track, the velocity relationship is shown in Figure 8. In the same way, we can derive the equivalent inertia Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 26 According to the Vector Superposition Rule, where is the vector angle between ⃑ and ⃑.
For , there is Furthermore, for , there is In addition, according to the law of the conservation of energy, By simultaneously considering formulae (13) to (17), we can get For the front part of the track, the velocity relationship is shown in Figure 8. In the same way, we can derive the equivalent inertia.    It is easy to calculate the inertia of other parts of the track by the same method.
Here, v tup , v tb , v t f , and v tr mean the absolute speed of the upper, lower, front, and rear parts of the track; v tupe , v tbe , v t f e , and v tre mean the traction velocity of the upper, lower, front, and rear parts of the track; v tupr , v tbr , v t f r , and v trr mean the relative velocity of the upper, lower, front, and rear parts of the track to driving wheel; m tup , m tb , m t f , and m tr mean the mass of the upper, lower, front, and rear parts of the track; and J tup , J tb , J t f , and J tr mean the equivalent inertia of the upper, lower, front, and rear parts of the track. Therefore, the equivalent inertia of the tracked part on the driving wheel can be presented as J o includes the inertia of the loading wheel, driving wheel, induced wheel, and supporting roller. Among them, the structure of the wheel is relatively special, because it has an I-shaped cross section, where R 1 is the outer diameter and R 2 is the inner diameter. It is shown in Figure 9.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 26 includes the inertia of the loading wheel, driving wheel, induced wheel, and supporting roller. Among them, the structure of the wheel is relatively special, because it has an I-shaped cross section, where R is the outer diameter and R is the inner diameter. It is shown in Figure 9. Using the same deduction method as above, we can get the deduction process as described in Formula (28-31), as follows: In the same way, we can get the equivalent inertia of other rollers: The equivalent inertia of the whole vehicle in the driving wheel can be obtained by synthesizing the formula above, that is, According to the output torque of the vehicle and the equivalent inertia of the tested vehicle, the target speed of the vehicle can be deduced.

Control Strategy
When testing the power system of a vehicle on the bench, it is the primary task of the bench to exert a load equal to the actual working condition. For the test bench in this paper, the task of the control system is to accurately simulate the driving resistance of the tracked vehicle [Error! Reference source not found.].
In order to achieve this function, the driving moment of the tested vehicle is measured by the torque sensor, and the theoretical rotational speed of the driving wheel can be obtained by substituting into Formula (6). When the combined torque of and output torque of the load Using the same deduction method as above, we can get the deduction process as described in Formula (28)-(31), as follows: v In the same way, we can get the equivalent inertia of other rollers: The equivalent inertia of the whole vehicle in the driving wheel can be obtained by synthesizing the formula above, that is, According to the output torque of the vehicle and the equivalent inertia of the tested vehicle, the target speed of the vehicle can be deduced.

Control Strategy
When testing the power system of a vehicle on the bench, it is the primary task of the bench to exert a load equal to the actual working condition. For the test bench in this paper, the task of the control system is to accurately simulate the driving resistance of the tracked vehicle [33].
In order to achieve this function, the driving moment T e of the tested vehicle is measured by the torque sensor, and the theoretical rotational speed ω em of the driving wheel can be obtained by substituting T e into Formula (6). When the combined torque of T e and output torque of the load analog motor makes the system rotate, the current rotational speed ω is obtained by the sensors and a controller C is designed to make ω infinitely approach ω em , so that the simulation of vehicle driving resistance can be realized. The system control block diagram is shown in Figure 10, where G(s) and J(s) are the transfer function of the motor and transmission mechanism, respectively.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 26 driving resistance can be realized. The system control block diagram is shown in Figure 10, where G(s) and J(s) are the transfer function of the motor and transmission mechanism, respectively. + -+ Figure 10. Control block diagram of the load simulation system.
As can be seen from the above figure, the control method is actually used to track the speed response of the system under the current torque. However, in the actual operation process, because the system is affected by environmental factors and friction coefficients at different velocities, a nonlinear uncertainty ∆(s) exists in its transfer function. In addition, the system loop is also affected by the driving torque , which is an external additional torque for the whole load simulation system. [Error! Reference source not found.,Error! Reference source not found.]. However, in practice, the system is driven by the combined torque of and . For the control loop, is not controlled by the system and is controlled by the external controller on the test vehicle or throttle. Therefore, for the load simulation system, is an uncontrollable random variable, and when it changes dramatically, it is very likely that the system will deviate from the target speed. Because of this, when ∆(s) changes or mutations occur, the stability of the system is problematic. In this regard, it is necessary to design a controller C which can not only adjust the loading and transmission system, but also respond quickly and follow under the action of system disturbance torque. There is an intuitive idea that the robust control method can be used to design the controller C(s) by taking the maximum disturbance and the maximum uncertainty of the test-bed system as the disturbance rejection index of the system and combining it with the tracking requirement of [Error! Reference source not found.,Error! Reference source not found.]. However, for a single controller, the responsiveness and strong anti-jamming ability are contradictory. Therefore, according to the structural invariance principle, the system is designed as a 2-degree-of-freedom loop (2-DOF), as shown in Figure 11. In the figure, ξ is the sensor noise, ( ) is the nominal transfer function of the system, ( ) = ( ) • ( ), and Q(s) is the filter to be designed. The disturbance observer (DOB) composed of Q(s) and ( ) in the dashed line frame is used to suppress the disturbance of the system, enhance the robustness of the system, and make the system transfer function close to the nominal transfer function ( ). The DOB and the controlled system constitute the inner loop of the system [Error! Reference source not found.]. The outer loop is controlled by controller C. Since the disturbance observer in the inner loop has suppressed the disturbance and noise of the system, the outer loop controller only needs to control the response speed and steadystate error of the system. According to the design requirements and test experience, controller C needs to control the response time of the system within 100ms and the tracking error within 1%. As can be seen from the above figure, the control method is actually used to track the speed response of the system under the current torque. However, in the actual operation process, because the system is affected by environmental factors and friction coefficients at different velocities, a non-linear uncertainty ∆(s) exists in its transfer function. In addition, the system loop is also affected by the driving torque T e , which is an external additional torque for the whole load simulation system. [34,35]. However, in practice, the system is driven by the combined torque of T e and T e . For the control loop, T e is not controlled by the system and is controlled by the external controller on the test vehicle or throttle. Therefore, for the load simulation system, T e is an uncontrollable random variable, and when it changes dramatically, it is very likely that the system will deviate from the target speed. Because of this, when ∆(s) changes or T e mutations occur, the stability of the system is problematic. In this regard, it is necessary to design a controller C which can not only adjust the loading and transmission system, but also respond quickly and follow ω em under the action of system disturbance torque.
There is an intuitive idea that the robust control method can be used to design the controller C(s) by taking the maximum disturbance and the maximum uncertainty of the test-bed system as the disturbance rejection index of the system and combining it with the tracking requirement of ω em [36,37]. However, for a single controller, the responsiveness and strong anti-jamming ability are contradictory. Therefore, according to the structural invariance principle, the system is designed as a 2-degree-of-freedom loop (2-DOF), as shown in Figure 11. In the figure, ξ is the sensor noise, G n (s) is the nominal transfer function of the system, G n (s) = G(s)·J(s), and Q(s) is the filter to be designed. The disturbance observer (DOB) composed of Q(s) and G n −1 (s) in the dashed line frame is used to suppress the disturbance of the system, enhance the robustness of the system, and make the system transfer function close to the nominal transfer function G n (s). The DOB and the controlled system constitute the inner loop of the system [38]. The outer loop is controlled by controller C. Since the disturbance observer in the inner loop has suppressed the disturbance and noise of the system, the outer loop controller only needs to control the response speed and steady-state error of the system. According to the design requirements and test experience, controller C needs to control the response time of the system within 100ms and the tracking error within 1%.  For the disturbance observer in the inner loop, because ( ) is usually regular, ( ) is irregular and not realizable, so it is necessary to introduce a filter Q(s) to make Q(s) ( ) realizable. Therefore, it is required that the relative order of Q(s) be greater than or equal to that of ( ), that is, ≥ . The closed-loop transfer function of the inner loop can be obtained when the effects of and ∆(s) on the system are collectively referred to as disturbances and defined as d.
It can be seen that the effect of disturbance and noise on the system results is determined by Q(s) and 1-Q(s). Here, Q (s) is the sensitivity function of the system and 1-Q(s) is the complementary sensitivity function of the system. For the systems, since disturbances usually occur in the lowfrequency band and noise usually occurs in the high-frequency band, it is necessary to design a lowpass filter Q(s) to suppress ξ in the high-frequency band and 1-Q (s) to suppress d in the lowfrequency band.
Using the above control method, a set of load simulation equipment consisting of a direct current motor (DC) and mechanical transmission equipment is simulated in this paper. The transfer function of the DC motor is known as The transfer function of the mechanical transmission equipment is For the disturbance observer in the inner loop, because G n (s) is usually regular, G n −1 (s) is irregular and not realizable, so it is necessary to introduce a filter Q(s) to make Q(s)G n −1 (s) realizable. Therefore, it is required that the relative order of Q(s) be greater than or equal to that of G n (s), that is, k Q ≥ k G . The closed-loop transfer function of the inner loop can be obtained when the effects of T e and ∆(s) on the system are collectively referred to as disturbances and defined as d.
It can be seen that the effect of disturbance and noise on the system results is determined by Q(s) and 1-Q(s). Here, Q(s) is the sensitivity function of the system and 1-Q(s) is the complementary sensitivity function of the system. For the systems, since disturbances usually occur in the low-frequency band and noise usually occurs in the high-frequency band, it is necessary to design a low-pass filter Q(s) to suppress ξ in the high-frequency band and 1-Q(s) to suppress d in the low-frequency band.
Using the above control method, a set of load simulation equipment consisting of a direct current motor (DC) and mechanical transmission equipment is simulated in this paper. The transfer function of the DC motor is known as The transfer function of the mechanical transmission equipment is By simultaneously considering the two formulas, we can get Based on the above conditions, the DOB for the system with the transfer function G n (s) is designed by using the H ∞ mixed sensitivity method. The optimal filter transfer function is The outer loop controller C(s) is designed as the PI controller, and the transfer function is According to the above model, a tracked vehicle with an equivalent inertia of 3180 kg·m 2 is simulated. The control model is established in Simulink and the speed response under a random torque signal is simulated, and the results are shown in Figure 12. Because of the large amplitude transformation of random signals, this type of signal also acts as the interference signal affecting the system, and the output speed deviates greatly from the theoretical speed under the control of the traditional 1DOF loop. Furthermore, by applying the 2DOF system with DOB described in this paper, it is found that the speed deviation of the system decreases significantly. Therefore, it can be seen that the control loop has an obvious effect on stabilizing the system and reducing the deviation. According to the above model, a tracked vehicle with an equivalent inertia of 3180 kg • m is simulated. The control model is established in Simulink and the speed response under a random torque signal is simulated, and the results are shown in Figure 12. Because of the large amplitude transformation of random signals, this type of signal also acts as the interference signal affecting the system, and the output speed deviates greatly from the theoretical speed under the control of the traditional 1DOF loop. Furthermore, by applying the 2DOF system with DOB described in this paper, it is found that the speed deviation of the system decreases significantly. Therefore, it can be seen that the control loop has an obvious effect on stabilizing the system and reducing the deviation.

Calibration of Electrical Inertia Simulation Accuracy
In order to study the accuracy of the electrical inertia simulation, the driving shaft and gearbox were used to connect the two ends of the tested motor to remove the intermediate vehicle, as shown in Figure 13. At this time, the right analog end simulates the inertia load by using the speed tracking method, and the left driving motor imposes a preset load. According to Formula 44, when the system torque is known, instantaneous acceleration of the system can be obtained by deriving the system speed, and the inertia of the instantaneous motor simulation can be obtained. Comparing the measured inertia with the standard inertia, the electric inertia can be calibrated.

Calibration of Electrical Inertia Simulation Accuracy
In order to study the accuracy of the electrical inertia simulation, the driving shaft and gearbox were used to connect the two ends of the tested motor to remove the intermediate vehicle, as shown in Figure 13. At this time, the right analog end simulates the inertia load by using the speed tracking method, and the left driving motor imposes a preset load. According to Formula (44), when the system torque is known, instantaneous acceleration of the system can be obtained by deriving the system speed, and the inertia of the instantaneous motor simulation can be obtained. Comparing the measured inertia with the standard inertia, the electric inertia can be calibrated.

Calibration of Standard Inertia
In the electrical inertia simulation test, the main factors affecting the accuracy of the test results are the accuracy of the inertia simulation and the system response time [Error! Reference source not found.]. In this paper, a set of calibrated high-precision standard inertia boxes are used to calibrate

Calibration of Standard Inertia
In the electrical inertia simulation test, the main factors affecting the accuracy of the test results are the accuracy of the inertia simulation and the system response time [39]. In this paper, a set of calibrated high-precision standard inertia boxes are used to calibrate the electric inertia.
Suppose that the inertia of one stage in the flywheel box and driving components to be calibrated is J s and J ba , respectively. Then, the inertia to be simulated by the motor is J a . To achieve J ba ± J a = J s , the inertia box and the tested motor are driven with the same torque to obtain the speed-up curve. By comparing the experimental results of the two groups, the calibration of the electric inertia simulation is completed.
The mechanical structure of the standard inertia (flywheel box) test is shown in Figure 14. According to Formula (44), when the input torque is constant, the speed of the system should rise uniformly in a constant angular acceleration α. Therefore, the standard inertia of the system can be calculated by the input torque and the measured velocity curve. In practical experiments, the formula is modified to T t − T d = Jα because of the resistance of the system, which is mainly related to the rotation speed [25,28]. When the rotation speed of the system is at a proper level, the resistance of the system is basically stable. Therefore, the angular acceleration can be measured by changing different driving torques, and the standard inertia of the system can be obtained by two sets of equations. 1 -electric motor, 2-reducer, 3-reducer, 4 -torque sensor, 5 -reducer, and 6 -tested motor Figure 13. Structure of the electric inertia calibration test bench.

Calibration of Standard Inertia
In the electrical inertia simulation test, the main factors affecting the accuracy of the test results are the accuracy of the inertia simulation and the system response time [Error! Reference source not found.]. In this paper, a set of calibrated high-precision standard inertia boxes are used to calibrate the electric inertia.
Suppose that the inertia of one stage in the flywheel box and driving components to be calibrated is and , respectively. Then, the inertia to be simulated by the motor is . To achieve ± = , the inertia box and the tested motor are driven with the same torque to obtain the speed-up curve. By comparing the experimental results of the two groups, the calibration of the electric inertia simulation is completed.  The mechanical structure of the standard inertia (flywheel box) test is shown in Figure 14. According to Formula (44), when the input torque is constant, the speed of the system should rise uniformly in a constant angular acceleration . Therefore, the standard inertia of the system can be calculated by the input torque and the measured velocity curve. In practical experiments, the formula is modified to − = because of the resistance of the system, which is mainly related to the rotation speed [25,Error! Reference source not found.]. When the rotation speed of the system is at a proper level, the resistance of the system is basically stable. Therefore, the angular acceleration can be measured by changing different driving torques, and the standard inertia of the system can be obtained by two sets of equations.  As shown in Figure 15, the driving torque values T t1 and T t2 are set separately at different stages of the same test, and the acceleration values α 1 and α 2 corresponding to the stage of higher speed curve steady speed rise ( r 1 ∼ r 2 )rpm are selected. Because the acceleration at the same speed is collected at different stages, it is considered that the resistance T di of each stage is equal. The inertia value J s of the standard inertia disc can be calculated by taking the difference between T t1 and T t2 . The accuracy of the standard inertia simulation is thus verified.

Calibration of Standard Inertia
In the electrical inertia simulation test, the main factors affecting the accuracy of the test results are the accuracy of the inertia simulation and the system response time [Error! Reference source not found.]. In this paper, a set of calibrated high-precision standard inertia boxes are used to calibrate the electric inertia.
Suppose that the inertia of one stage in the flywheel box and driving components to be calibrated is and , respectively. Then, the inertia to be simulated by the motor is . To achieve ± = , the inertia box and the tested motor are driven with the same torque to obtain the speed-up curve. By comparing the experimental results of the two groups, the calibration of the electric inertia simulation is completed.  The mechanical structure of the standard inertia (flywheel box) test is shown in Figure 14. According to Formula (44), when the input torque is constant, the speed of the system should rise uniformly in a constant angular acceleration . Therefore, the standard inertia of the system can be calculated by the input torque and the measured velocity curve. In practical experiments, the formula is modified to − = because of the resistance of the system, which is mainly related to the rotation speed [25,Error! Reference source not found.]. When the rotation speed of the system is at a proper level, the resistance of the system is basically stable. Therefore, the angular acceleration can be measured by changing different driving torques, and the standard inertia of the system can be obtained by two sets of equations.  The actual structure of the standard inertia box is a flywheel set with large inertia, and its inertia simulation range is (50-1000) kg·m 2 . In order to improve the upper limit of inertia simulation and provide more reference data for electric inertia simulation calibration, a gearbox equipped with speed ratios of 1:2.36 and 1:6.48 is equipped to improve the equivalent inertia of the equipment being marked. The maximum equivalent inertia of the original 1000 kg·m 2 is 42,006.60 kg·m 2 after considering the gear box deceleration ratio. That is to say, the inertia ranging from 255-42,000 kg·m 2 can be simulated by different combinations of inertia plates in the flywheel box. In this simulation method, although it is a simulation, through data processing, it can completely reflect the data curve within 250~41,000, to meet the requirements of comparison. In this way, the inertia calibration test of 400 kg·m 2 inertia box with a 1:2.36 transmission ratio is carried out. The test results are shown in Figure 16.
The actual structure of the standard inertia box is a flywheel set with large inertia, and its inertia simulation range is (50-1000) kg • m . In order to improve the upper limit of inertia simulation and provide more reference data for electric inertia simulation calibration, a gearbox equipped with speed ratios of 1:2.36 and 1:6.48 is equipped to improve the equivalent inertia of the equipment being marked. The maximum equivalent inertia of the original 1000 kg • m is 42006.60 kg • m after considering the gear box deceleration ratio. That is to say, the inertia ranging from 255-42000 kg • m can be simulated by different combinations of inertia plates in the flywheel box. In this simulation method, although it is a simulation, through data processing, it can completely reflect the data curve within 250~41000, to meet the requirements of comparison. In this way, the inertia calibration test of 400 kg • m inertia box with a 1:2.36 transmission ratio is carried out. The test results are shown in Figure 16.  Observing the two deceleration processes, it can be found that the speed decline curve of the system during the deceleration process is smooth, and there is no obvious fluctuation in the and of the torque curve section. The deceleration of the system is 0.24, 0.23, and 0.23 rad/s2, respectively, when analyzing the deceleration of the system at the intervals of 1, 2, and 3. Therefore, the influence of system resistance moment on the inertia calibration of the test-bed in a high speed range can be neglected.  Figure 16a is the overall data of the test while Figure 16b is the calibration section. In the test, the equivalent inertia of the flywheel is 2227.84 kg·m 2 . The average acceleration of the system is 1.3 rad/s 2 under the condition of average torque T t1 = 3140.09Nm, and 1.54 rad/s 2 under the condition of T t2 = 3674.7 Nm. Therefore, the inertia of the system is 2229.94 kg·m 2 and the measurement deviation is 0.19%.
Observing the two deceleration processes, it can be found that the speed decline curve of the system during the deceleration process is smooth, and there is no obvious fluctuation in the T d1 and T d2 of the torque curve section. The deceleration of the system is 0.24, 0.23, and 0.23 rad/s 2 , respectively, when analyzing the deceleration of the system at the intervals of 1, 2, and 3. Therefore, the influence of system resistance moment on the inertia calibration of the test-bed in a high speed range can be neglected.
The standard inertia of 750 kg·m 2 can be calibrated in the same way at the transmission ratio of 1:2.36, shown in Figure 17. The equivalent inertia is 4177.2 kg·m 2 . Under the condition of average torque T t1 = 4220.38 Nm, the average acceleration of the system is 0.94 rad/s 2 , while under the condition of T t2 = 4919.48 Nm, the average acceleration of the system is 1.11 rad/s 2 . Therefore, the inertia of the system is 4157.79 kg·m 2 and the measurement deviation is 0.46%. The standard inertia of 250 kg·m 2 can be calibrated at the transmission ratio of 1:6.48, shown in Figure 18. The equivalent inertia is 10,497.6 kg·m 2 . Under the condition of average torque T t1 = 7206.53 Nm, the average acceleration of the system is 0.67 rad/s 2 , while under the condition of T t2 = 4919.48 Nm, the average acceleration of the system is 0.73 rad/s 2 . The inertia of the system is 10,522.34 kg·m 2 and the measurement error is 0.23%. The variation of standard inertia deviation with the increase of equivalent inertia is shown in Figure 19. It is evident that the mechanical inertia under large inertia is less than 0.5%, except in the case of small inertia, where the inertia calibration deviation is too large due to fluctuation of the sensor signal.

Deviation(%)
Equivalent inertia(kg·m 2 ) Figure 19. Deviation of standard inertia. The variation of standard inertia deviation with the increase of equivalent inertia is shown in Figure 19. It is evident that the mechanical inertia under large inertia is less than 0.5%, except in the case of small inertia, where the inertia calibration deviation is too large due to fluctuation of the sensor signal. The standard inertia of 750 kg • m 2 can be calibrated in the same way at the transmission ratio of 1:2.36, shown in Figure 17. The equivalent inertia is 4177.2 kg • m 2 . Under the condition of average torque 1 = 4220.38Nm , the average acceleration of the system is 0.94rad/s 2 , while under the condition of 2 = 4919.48Nm, the average acceleration of the system is 1.11 rad/s 2 . Therefore, the inertia of the system is 4157.79 kg • m 2 and the measurement deviation is 0.46%. The standard inertia of 250 kg • m 2 can be calibrated at the transmission ratio of 1:6.48, shown in Figure 18. The equivalent inertia is 10497.6 kg • m 2 . Under the condition of average torque 1 = 7206.53 Nm, the average acceleration of the system is 0.67 rad/s 2 , while under the condition of 2 = 4919.48 Nm, the average acceleration of the system is 0.73 rad/s 2 . The inertia of the system is 10522.34 kg • m 2 and the measurement error is 0.23%.

Calibration of Basic Inertia
The variation of standard inertia deviation with the increase of equivalent inertia is shown in Figure 19. It is evident that the mechanical inertia under large inertia is less than 0.5%, except in the case of small inertia, where the inertia calibration deviation is too large due to fluctuation of the sensor signal.

Calibration of Basic Inertia
Through the standard inertia test, the inertia J = J s of each test system can be determined. Because the inertia simulation equipment itself has inertia J ba (which is called the basic inertia of the system for convenience of expression), it is necessary to define the basic inertia J ba of the equipment before the electrical inertia simulation test, so that the simulation inertia J α of the motor can be defined.
To calibrate the basic inertia, mechanical connection is adopted, as shown in Figure 20. The to-be-calibrated inertia includes the tested motor (6), reducer (5), torque meter (4), and coupling. During the test, the driving end is in torque control mode, and the calibrated end is powered off. Under this structure, the basic inertia J ba of the system can be obtained by the same principle as the standard inertia.
To calibrate the basic inertia, mechanical connection is adopted, as shown in Figure 20. The tobe-calibrated inertia includes the tested motor (6), reducer (5), torque meter (4), and coupling. During the test, the driving end is in torque control mode, and the calibrated end is powered off. Under this structure, the basic inertia of the system can be obtained by the same principle as the standard inertia.
1 electric motor, 2-reducer, 3-reducer, 4 -torque sensor, 5 -reducer, and 6 -tested motor. The test result is shown in Figure 21. According to the design and manufacture drawings, the theoretical basic inertia is 512 kg • m . The average acceleration of the system is 1.79 rad/s 2 during the speed 150~200 r/min under the condition of average torque = 1500Nm, and 2.02 rad/s 2 under the torque of = 1620Nm. The actual inertia is 523.96 kg • m after calibration.

Calibration of Electrical Inertia
In the standard inertia calibration test, the test result of 400 kg • m inertia flywheel set at a 1:2.36 transmission ratio is 2232.05 kg • m . As the basic inertia is known to be 523.96 kg • m , the inertia to be simulated by the electric motor is 1708.09 kg • m . After calibrating in the same way as the standard inertia calibration, the experimental results are shown in Figure 22. The system inertia is 2211.09 kg • m and the percentage deviation of the electrical inertia simulation is 0.95%. The test result is shown in Figure 21. According to the design and manufacture drawings, the theoretical basic inertia is 512 kg·m 2 . The average acceleration of the system is 1.79 rad/s 2 during the speed 150~200 r/min under the condition of average torque T t1 = 1500Nm, and 2.02 rad/s 2 under the torque of T t2 = 1620 Nm. The actual inertia is 523.96 kg·m 2 after calibration.
Time Time Figure 21. Test results of basic inertia calibration.

Calibration of Electrical Inertia
In the standard inertia calibration test, the test result of 400 kg·m 2 inertia flywheel set at a 1:2.36 transmission ratio is 2232.05 kg·m 2 . As the basic inertia is known to be 523.96 kg·m 2 , the inertia to be simulated by the electric motor is 1708.09 kg·m 2 . After calibrating in the same way as the standard inertia calibration, the experimental results are shown in Figure 22. The system inertia is 2211.09 kg·m 2 and the percentage deviation of the electrical inertia simulation is 0.95%.  With the same method, the equivalent inertia of 750 kg·m 2 and 250 kg·m 2 in the standard inertia calibration test is 4157.69 kg·m 2 and 10,522.34 kg·m 2 respectively. Through the calibration shown below, the inertia and percentage deviation of the two groups are 4183.43 kg·m 2 (0.61%) and 10,603.38 kg·m 2 (0.76%), respectively, shown in Figure 23. The variation of electrical inertia simulation deviation with the increase of equivalent inertia is shown in Figure 24. It can be found that in the range of inertia of 1000~30000 kg • m , the deviation is less than 1%. The trend of simulation deviation of electrical inertia is the same as that of standard inertia.

Calibration of Electrical Inertia Response Time
When the motor is used to simulate mechanical inertia, there must be a delay problem, so the response time of the system needs to be measured [Error! Reference source not found.,Error! Reference source not found.]. According to T − T = Jα, it can be seen that when the driving torque T of the system is set as a constant value and the resistance moment T of the system is basically stable, the angular acceleration α is a constant value. The corresponding angular acceleration curve can be obtained by measuring the rotational speed of the system. When the angular acceleration of the system is stable, the electrical inertia simulated by the system is stable. Therefore, the response The variation of electrical inertia simulation deviation with the increase of equivalent inertia is shown in Figure 24. It can be found that in the range of inertia of 1000~30,000 kg·m 2 , the deviation is less than 1%. The trend of simulation deviation of electrical inertia is the same as that of standard inertia. With the same method, the equivalent inertia of 750 kg • m 2 and 250 kg • m 2 in the standard inertia calibration test is 4157.69 kg • m 2 and 10522.34 kg • m 2 , respectively. Through the calibration shown below, the inertia and percentage deviation of the two groups are 4183.43 kg • m 2 (0.61%) and 10603.38 kg • m 2 (0.76%), respectively, shown in Figure 23. The variation of electrical inertia simulation deviation with the increase of equivalent inertia is shown in Figure 24. It can be found that in the range of inertia of 1000~30000 kg • m 2 , the deviation is less than 1%. The trend of simulation deviation of electrical inertia is the same as that of standard inertia.

Calibration of Electrical Inertia Response Time
When the motor is used to simulate mechanical inertia, there must be a delay problem, so the response time of the system needs to be measured [Error! Reference source not found.,Error! Reference source not found.]. According to T t − T d = Jα, it can be seen that when the driving torque T t of the system is set as a constant value and the resistance moment T d of the system is basically stable, the angular acceleration α is a constant value. The corresponding angular acceleration curve can be obtained by measuring the rotational speed of the system. When the angular acceleration of the system is stable, the electrical inertia simulated by the system is stable. Therefore, the response

Calibration of Electrical Inertia Response Time
When the motor is used to simulate mechanical inertia, there must be a delay problem, so the response time of the system needs to be measured [40,41]. According to T t − T d = Jα, it can be seen that when the driving torque T t of the system is set as a constant value and the resistance moment T d of the system is basically stable, the angular acceleration α is a constant value. The corresponding angular acceleration curve can be obtained by measuring the rotational speed of the system. When the angular acceleration of the system is stable, the electrical inertia simulated by the system is stable. Therefore, the response time of the electrical inertia simulation of the measurement system can be obtained by observing the angular acceleration (or torque) curve.
As shown in Figure 25, the electric inertia is calibrated by the steady acceleration test. Under the stable torque condition, the system uses the motor to follow the target speed through the step acceleration test. In the dotted line amplification area of Figure 25a (i.e., Figure 25b), the target velocity jumps from step to 10 r/min at 20.950 s, and then the acceleration of the system starts to respond at 20.952 s, and returns to stability at 21.022 s. By evaluating the acceleration response, the system response time of electrical inertia simulation is 72 ms.
The above calibration method can be used to calibrate the response time of electric inertia under steady state conditions. However, in the process of driving load simulation, the speed and torque of the system are often unsteady, so the above calibration results cannot express the response rate of electric inertia in a real vehicle simulation environment.
According to the speed tracking method, the system deduces the target speed ω em (s) of the system through the actual system model G em (s), and then calculates the required motor torque T l (s) by the speed closed loop. According to a series of standard inertia tests, the change of T d in a continuous speed range can be neglected. When the output of load torque is accurate, the response time of the system can be evaluated by comparing the target acceleration α t with the measured acceleration α r .
At the moment t 0 of a continuous torque fluctuation curve, the target acceleration of the system is α t , assuming that after t seconds α r −α t α r < 1%, the response time of the signal is considered to be t. time of the electrical inertia simulation of the measurement system can be obtained by observing the angular acceleration (or torque) curve. As shown in Figure 25, the electric inertia is calibrated by the steady acceleration test. Under the stable torque condition, the system uses the motor to follow the target speed through the step acceleration test. In the dotted line amplification area of Fig. 25.1 (i.e., Fig. 25.2), the target velocity jumps from step to 10 r/min at 20.950 s, and then the acceleration of the system starts to respond at 20.952 s, and returns to stability at 21.022 seconds. By evaluating the acceleration response, the system response time of electrical inertia simulation is 72 ms. The above calibration method can be used to calibrate the response time of electric inertia under steady state conditions. However, in the process of driving load simulation, the speed and torque of the system are often unsteady, so the above calibration results cannot express the response rate of electric inertia in a real vehicle simulation environment.
According to the speed tracking method, the system deduces the target speed ( ) of the system through the actual system model ( ), and then calculates the required motor torque ( ) by the speed closed loop. According to a series of standard inertia tests, the change of in a continuous speed range can be neglected. When the output of load torque is accurate, the response time of the system can be evaluated by comparing the target acceleration with the measured acceleration α . At the moment 0 of a continuous torque fluctuation curve, the target acceleration of the system is α , assuming that after t seconds | α −α | < 1%, the response time of the signal is considered to be t. In the actual operation process, because the sampling rate of the system is much higher than the response time of the system, the single response time adopts the following logical judgment method, in which ∆t is the sampling frequency of the system. After counting m times, the average response time of the system is 1 ∑ =1 , as shown in Figure 26. In the actual operation process, because the sampling rate of the system is much higher than the response time of the system, the single response time adopts the following logical judgment method, in which ∆t is the sampling frequency of the system. After counting m times, the average response time of the system is 1 m m i=1 t i , as shown in Figure 26. As can be seen in Figure 27, the response time is represented by a gray histogram. According to the logical judgment of a single response time, the average response time of the system is 25.40 ms in a preset spectrum test with continuous load variation. The slowest response of system torque occurs in several periods when the torque fluctuates greatly. At this time, the system speed fluctuates repeatedly. Because of the system lag, it is difficult for the system to respond in time, and the maximum response time is 68 ms.
Using the above methods, the response time of the control system can be calibrated in the simulation process. At present, in order to speed up the response of the system, lag compensation (such as Smith prediction), auto-disturbance rejection (ADRC), and several intelligence algorithms (such as the model reference adaptive algorithm) are generally used in the control system. The real system response speed after using this control method can be calibrated by the method described in this paper. Appl. Sci. 2019, 9,   As can be seen in Figure 27, the response time is represented by a gray histogram. According to the logical judgment of a single response time, the average response time of the system is 25.40ms in a preset spectrum test with continuous load variation. The slowest response of system torque occurs in several periods when the torque fluctuates greatly. At this time, the system speed fluctuates repeatedly. Because of the system lag, it is difficult for the system to respond in time, and the maximum response time is 68 ms.
Using the above methods, the response time of the control system can be calibrated in the simulation process. At present, in order to speed up the response of the system, lag compensation (such as Smith prediction), auto-disturbance rejection (ADRC), and several intelligence algorithms (such as the model reference adaptive algorithm) are generally used in the control system. The real system response speed after using this control method can be calibrated by the method described in this paper.

Correction of the Test Bed
Before the test, the calibration results should be used to correct the system. It is known that there is a tracked vehicle with a mass of 30t. The equivalent inertia J = 2723.6 kg • m is obtained from Formula (34), and the input inertia of the system corrected from Figure 24 should be 2703.05 kg • m . The system is simulated with the inertia before and after correction as the input, and the experimental curve shown in Figure 28 is obtained. The results are shown in Table 1.  As can be seen in Figure 27, the response time is represented by a gray histogram. According to the logical judgment of a single response time, the average response time of the system is 25.40ms in a preset spectrum test with continuous load variation. The slowest response of system torque occurs in several periods when the torque fluctuates greatly. At this time, the system speed fluctuates repeatedly. Because of the system lag, it is difficult for the system to respond in time, and the maximum response time is 68 ms.
Using the above methods, the response time of the control system can be calibrated in the simulation process. At present, in order to speed up the response of the system, lag compensation (such as Smith prediction), auto-disturbance rejection (ADRC), and several intelligence algorithms (such as the model reference adaptive algorithm) are generally used in the control system. The real system response speed after using this control method can be calibrated by the method described in this paper.

Correction of the Test Bed
Before the test, the calibration results should be used to correct the system. It is known that there is a tracked vehicle with a mass of 30t. The equivalent inertia J = 2723.6 kg • m is obtained from Formula (34), and the input inertia of the system corrected from Figure 24 should be 2703.05 kg • m . The system is simulated with the inertia before and after correction as the input, and the experimental curve shown in Figure 28 is obtained. The results are shown in Table 1.

Correction of the Test Bed
Before the test, the calibration results should be used to correct the system. It is known that there is a tracked vehicle with a mass of 30 t. The equivalent inertia J = 2723.6 kg·m 2 is obtained from Formula (34), and the input inertia of the system corrected from Figure 24 should be 2703.05 kg·m 2 . The system is simulated with the inertia before and after correction as the input, and the experimental curve shown in Figure 28 is obtained. The results are shown in Table 1 (1) before (2) after    The results show that under the same simulated object, the corrected input inertia reduces the simulation error from 0.96% to 0.27%, and effectively improves the accuracy of electrical inertia simulation.

Load Test and Shift Performance Test
A load test is used to simulate the torque output of a vehicle under a specific road resistance coefficient. The accuracy of vehicle inertia modeling can be verified by a load test. Load tests were carried out on a tracked vehicle with a mass of 30 t. According to Formula (34), the equivalent inertia is J = 2723.6 and the road resistance coefficient is set to 0.05, so the road resistance moment is −4204 Nm. Figure 29a shows the gearshift process of the tracked vehicle on a flat road. Figure 29b is a local amplification of uniform acceleration in a fixed gear. Its acceleration is 0.65 km/h, i.e., a = 0.63 rad/s 2 . The required torque under the above equivalent inertia and deceleration is as follows: (1) before (2) after Figure 28. Comparison of results before and after correction. show that under the same simulated object, the corrected input inertia reduces the simulation error from 0.96% to 0.27%, and effectively improves the accuracy of electrical inertia simulation.

Load Test and Shift Performance Test
A load test is used to simulate the torque output of a vehicle under a specific road resistance coefficient. The accuracy of vehicle inertia modeling can be verified by a load test. Load tests were carried out on a tracked vehicle with a mass of 30 t. According to Formula (34), the equivalent inertia is J = 2723.6 and the road resistance coefficient is set to 0.05, so the road resistance moment is -4204 Nm. Therefore, the required torque in the above period should be T = 4204 + 1724.03 = 5928.03 Nm While the actual average value is 5957.2 Nm and the relative error is 0.49%. Figure 29c is a local zoom-in of the upshifting. In a normal shifting process, the power is released before shifting, so the system has instantaneous deceleration; after a high gear, due to the increase of the clutch front-end speed, explained by the clutch disc engagement process due to the speed difference between the front and rear discs, instantaneous large torque mutation occurs. Figure 29d is the local enlargement of downshifting. When the clutch is released after shifting, the speed of the front end of the clutch is greater than that of the back end, so the output torque has a sudden change process. In the process of the rapid change of torque, the difference between the actual speed and the target speed becomes larger due to the lag of loading of the control system. However, after the shift is completed, the speed error converges rapidly, and there is no big fluctuation in the shift process of the speed tracking control bench test. The maximum speed error is controlled within 2 km/h. The load simulation test results show that the accuracy of the equipment simulation vehicle engine load is controlled within 1%, and thus has a high accuracy.

Longitudinal Slope Passability Test
A longitudinal slope passability test is used to measure and study the passability of different gradients. Through the test, the load of vehicle transmission equipment on different gradient pavement and whether the vehicle skids due to insufficient adhesion can be studied.
The same tracked vehicle as that from the load test is used. The resistance coefficient of the road surface is set to be f = 0.05. As shown in Figure 30, the designed vehicle passes through the climbing section with the gradient of 4 to 20 degrees and the downhill section with the gradient of 20 to 2 degrees. By enlarging the enlarged area of the 10 to 5 degrees section in the figure, it can be found that the vehicle will enter an adjustment area where the speed rises to slow down and the torque falls to increase when the slope changes. This phenomenon accords with the actual situation of vehicles. When the deviation between the actual speed and the target speed increases in the adjustment area, because the car body is in a dynamic change process at this time, the car body inclination angle is constantly changing. At this time, the system load suddenly falls back and the drop rate increases, which causes the system loading lag to cause the error to increase. The simulated working conditions of vehicles in five simulated area have been evaluated. At this time, the average resistance torque of the road slope is 11,528 Nm. At this time, the average load of vehicles is 11,573 Nm, and the simulation error of equipment is 0.39%, which verifies the accuracy of the dynamic model.
When the deviation between the actual speed and the target speed increases in the adjustment area, because the car body is in a dynamic change process at this time, the car body inclination angle is constantly changing. At this time, the system load suddenly falls back and the drop rate increases, which causes the system loading lag to cause the error to increase. The simulated working conditions of vehicles in five simulated area have been evaluated. At this time, the average resistance torque of the road slope is 11528 Nm. At this time, the average load of vehicles is 11573 Nm, and the simulation error of equipment is 0.39%, which verifies the accuracy of the dynamic model. (2) Test data in enlaged area

Acceleration Test
For general tracked vehicles, the acceleration capability of 0-32 km/h indicates their acceleration characteristics. In this study, in order to study the acceleration capability of tracked vehicles at different speed stages, six experiments at different speeds were carried out respectively. The maximum steady speed of tracked vehicles was evenly distributed from 9.8 km/h to 58.85 km/h. Figure 31 is the speed and torque diagram and Table 2 presents the data.

Acceleration Test
For general tracked vehicles, the acceleration capability of 0-32 km/h indicates their acceleration characteristics. In this study, in order to study the acceleration capability of tracked vehicles at different speed stages, six experiments at different speeds were carried out respectively. The maximum steady speed of tracked vehicles was evenly distributed from 9.8 km/h to 58.85 km/h. Figure 31 is the speed and torque diagram and Table 2  (2) Test data in enlaged area

Acceleration Test
For general tracked vehicles, the acceleration capability of 0-32 km/h indicates their acceleration characteristics. In this study, in order to study the acceleration capability of tracked vehicles at different speed stages, six experiments at different speeds were carried out respectively. The maximum steady speed of tracked vehicles was evenly distributed from 9.8 km/h to 58.85 km/h. Figure 31 is the speed and torque diagram and Table 2 presents the data.    According to the experimental results, the tracked vehicle has the greatest acceleration capacity in the 5-40 km/h stage. When the speed exceeds 40 km/h, the acceleration capacity of the vehicle decreases due to the increase of the gear. The experimental results are in accordance with the actual vehicle acceleration performance.

Conclusions
According to the test requirement of a tracked vehicle, this paper has formulated a method that can be employed to build the equivalent inertia model of a tracked vehicle. Based on the model, a 2DOF control loop combined with a disturbance observer is proposed, which solves the stable problem due to the uncertainty and disturbance in the model. In order to discuss the accuracy of inertia simulation under the control program, this paper innovatively proposes a calibration method of inertia simulation. The test results can be used to draw the following conclusions: (1) The 2DOF control loop combined with a disturbance observer has a good accuracy and fast response speed; (2) The disturbance observer can effectively solve the system stability problem under model uncertainty and system disturbance; (3) The calibration method of electrical inertia simulation using a standard inertia ratio has a good feasibility, quantitative analysis of the accuracy of mechanical inertia point simulation, and a high reference value for evaluating the control system and mechanical system; (4) The inertia simulation accuracy of the calibrated inertia simulation system of tracked vehicles is accurate and reliable, which provides a reliable test basis for the design and development of tracked vehicles.