Abstract
To improve the control precision of the drive system of hydraulic tracked vehicles, we established a mathematical model of the drive system based on the analysis of structural characteristics of the high-clearance hydraulic tracked vehicles and the dual-pump dual-motor drive system and developed a control strategy based on the quantitative feedback theory. First, the mutual independence of the two motor channels was achieved through channel decoupling. Then, the loop-shaping controller and the pre-filter were designed for the two channels. The result of a simulation experiment indicates that the proposed control method is very effective in suppressing external uncertainties and smoothening the speed-switching process of the hydraulic motor. Finally, an hydraulic tracked vehicle steering experimental test was carried out. The results show that under two different steering modes, the maximum standard deviation of the output speeds of the inner and outer motors of the hydraulic tracked vehicle is only 0.42, which meets the performance requirement on the hydraulic motor speed. The average steering track radii of the geometric centers of the inner and outer tracks are 1.828 and 0.033 m, respectively, and the relative errors are 1.56% and 3.19%, respectively. This demonstrates that the proposed control method achieves satisfactory results in the robust control of the hydraulic tracked vehicle drive system. It provides some references for the future control research of the hydraulic servo drive system of the high-clearance hydraulic tracked vehicles.
Keywords
Introduction
Because of the strong-grade ability, large ground-contact area, mild soil-compaction effect, and high stability, the use of high-clearance hydraulic tracked vehicles (HTVs) has been widely increasing in small patches of farmland and hilly areas for fertilization, weeding, and pesticide spraying purposes. High labor intensity and poor working environment characterize vehicle-based farming activities related to high stalk crops. In particular, the health of the operator is put at risk during the pesticide spraying operation due to the easy-to-spread nature of pesticide sprays. 1 Therefore, developing highly automated tracked vehicles has become a trend in the research field of agriculture machinery. Simple and compact structure and high reliability are the basic requirements of vehicle automation. To ensure the controllability and operability of the vehicles, a coupling control strategy of force and displacement,2–4 a H-extension controller based on active front steering system in extreme conditions, 5 and optimally adapt automatic control of intelligent electric vehicles to driving styles 6 had been proposed successively by some experts. With the development of technology and the ever-higher expectation of the users on the control performance of tracked vehicles, some techniques like these will be increasingly used in agricultural tracked vehicles.
The electro-hydraulic servo system, combining the advantages of electronic and hydraulic control, is characterized by high control precision, fast response, high torque output, and flexible signal processing, and the ease of implementing feedbacks of various parameters, and it is a good choice to realize the automation of agricultural tracked vehicle drive system. However, because of several factors such as the oil compressibility, nonlinear friction, time variation, possible load disturbances, and unpredictable parameter perturbation, the hydraulic servo system may not fulfill its purpose of high-precision control.7,8 To explore how to enable robust control with the hydraulic servo system under the negative influence of model uncertainty, some scholars have developed adaptive fuzzy proportional-integral-derivative (PID) algorithm for online real-time adjustment of PID parameters, 9 the genetic algorithm for adaptive adjustment of PID control parameters, 10 the anti-interference adaptive control scheme based on full-state feedback, 11 the sliding mode controller which is immune from the influence of parameter variation and external disturbances of electro-hydraulic servo system, 12 the model reference adaptive controller (MRAC) based on mathematical model, 13 and the adaptive neural network controller, 14 effectively improving the steady-state control precision and dynamic response speed. There also have been developed some hybrid controllers such as an adaptive integral back-stepping controller, 15 a novel sliding mode controller, 16 and predictive controllers. 17 However, hydraulic servo systems are time-varying systems. As the conventional control methods are usually suitable for slow time-varying systems, most of them are too complicated to be directly applied to hydraulic servo systems18,19 due to poor transient response performance under complicated working conditions.20,21 Therefore, it is especially necessary to develop control strategies suitable for tracked vehicle hydraulic servo systems.
Quantitative feedback theory (QFT) is a robust controller design method especially suitable for systems with uncertainties.22,23 In the practice of control system design, applying QFT helps maintain the balance between various performance indicators of the system, which provides an effective method for adjusting the controller. Initially, the QFT was used in the design of the flight control system; 24 at present, its application has expanded to many fields such as process control, 25 vehicle control,26,27 robot control, wind power system control,28,29 and spacecraft control. 30
This study applies the QFT to develop a drive system control strategy based on the analysis of structural characteristics of the high-clearance HTVs and the dual-pump dual-motor drive system. The proposed control strategy aims at solving some steering-related control problems caused by parameter uncertainties, time-varying, and nonlinear nature of the hydraulic drive system. First, a mathematical model of the drive system was established, and the mutual independence of the two motor channels was achieved through channel decoupling. Then, the loop-shaping controller and the pre-filter were designed for the two channels. The simulation results prove that the proposed control method is very effective in suppressing external uncertainties and outperforms the conventional PID control method in smoothening the speed-switching process of a hydraulic motor. The results of an HTV steering experiment show that the proposed control method enables high-precision steering in two steering modes, fulfilling its purpose of providing a robust control for the drive system of tracked vehicle.
System modeling and analysis
Description of the drive system of the tracked vehicle
Figure 1 shows the drive system of the tracked vehicle. The drive system uses a dual-pump dual-motor throttle governing circuit; it is a typical electro-hydraulic servo valve-controlled motor system. The drive system mainly consists of engine, double pump, hydraulic motor, servo valve, balance valve, and three-way compensator. According to the figure, two hydraulic motors (1) are connected with the two corresponding tracks (7). When the tracked vehicle is moving, the movement direction can be controlled by adjusting the speeds of the two motors, and the speed of each motor can be adjusted by changing the opening of the corresponding servo valve. The drive system can drive the tracked vehicle to perform various actions ranging from moving in a straight line to in-situ steering.

Drive system of the tracked vehicle.
Mathematical model of the tracked vehicle drive system
To design a robust controller for the tracked vehicle drive system, it is necessary to establish a mathematical model of the drive system at first. The model consists of a hydraulic servo system mathematical model and a tracked vehicle dynamics model.
Mathematical model of the hydraulic servo system
The tracked vehicle drive system consists of two symmetrical four-way valve-controlled motor oil circuits. The linearized flow equations for the servo valve ports are
The hydraulic motor flow continuous equations are
The moment balance equations of hydraulic motors are
where QL is the load flow of servo valve, Xv is the spool displacement of servo valve, Kq is the flow gain of servo valve, Kc is the flow pressure coefficient of servo valve, PL is the load pressure drop, Dm is the hydraulic motor’s displacement, θm is the rotation angle of hydraulic motor, Ct is the leakage coefficient of hydraulic motor, Vt is the total volume of the two chambers of the hydraulic motor and the connected pipes, J1 and J2 are the inertia values of hydraulic motors, respectively, βe is the oil volume elastic modulus, and T is the external load moment applied to the hydraulic motor shaft.
Dynamics model of the tracked vehicle
Moving in a straight line and in-situ steering are the critical points of the steering motion of the tracked vehicle, so it is only necessary to establish the steering dynamics model of the tracked vehicle. The force condition of the tracked vehicle is relatively complicated during steering movement.31,32,33 For the convenience of analysis, this study makes the following assumptions
The center of mass of the tracked vehicle coincides with the geometric center.
The pressure of the track is evenly distributed across the ground-contact area, and the influence of the track’s width is ignored.
As the speed of the tracked vehicle is low, the air resistance and the centrifugal force are ignored.
The shearing resistance and bulldozing resistance of the tracked vehicle are ignored. Thus, a simplified steering movement model of the tracked vehicle is established (Table 1).
Parameters of the tracked vehicle.
According to Figure 2, the steering motion of the tracked vehicle can be decomposed into a component parallel to the motion of the vehicle plane center C and a component circling the center C, and steering is realized by creating a speed difference between the two tracks at the two sides. When the inner and outer tracks are rolling in the same direction at the same speed, the tracked vehicle is moving in a straight line. When the inner and outer tracks are rolling in opposite directions with the same speed, the tracked vehicle is doing in-situ steering. Given that the tracked vehicle experiences different force conditions in different driving states, the steering motions can be divided into two types: R ≥ B/2 and 0 ≤ R < B/2.

The two types of steering motion of the tracked vehicle: (a) (R ≥ B/2) and (b) (0 ≤ R < B/2).
When the tracked vehicle is turning, the theoretical steering radius can be calculated using equation (7) if no slipping occurs during steering
The speeds and steering angular speeds of the two tracks during steering are
When R ≥ B/2, both the inner and outer tracks are rolling in the same direction (see Figure 2(a)). The braking force F1 of the inner track and the traction force F2 of the outer track can be expressed as
When 0 ≤ R < B/2, both the inner and outer tracks are rolling in opposite directions, and the rolling direction of the outer track is consistent with that of the vehicle (see Figure 2(b)). At this time, the driving forces applied to the inner and outer tracks are of the same magnitude and opposite directions. The expressions of the driving forces are
When the tracked vehicle is turning at a radius of R ≥ B/2, the dynamic balance equation is
When the tracked vehicle is turning at a radius of 0 ≤ R < B/2, the inner and outer tracks are rolling in opposite directions. As a result, the driving force and the frictional force of the inner track are also reversed. So, the dynamic balance equation of the tracked vehicle turning at a radius of 0 ≤ R < B/2 can be obtained by adding negative signs in front of F1 and FR1
Parameter equation of hydraulic motor
When the tracked vehicle is not experiencing slipping during driving, the torque is
where i = 1, 2. T1 and T2 correspond to the torques of the inner motor and outer motor, respectively (N m), i is the side transmission ratio, η is the efficiency of the section from the motor output shaft to the track, and r is the driving wheel radius (m).
The relationship between the motor speed and the steering angular speed when the tracked vehicle is turning can be obtained using equations (8) and (9)
When the tracked vehicle is turning at a radius of R ≥ B/2, the driving force can be calculated using equations (15) and (16)
The torque balance equations of hydraulic motors can be obtained using equations (20)–(24)
Similarly, when the tracked vehicle is turning at a radius of 0 ≤ R < B/2, the torque equations of the hydraulic motors are
Design of QFT controller
In this study, the controller is designed based on QFT. First, the uncertain parameter disturbances and external disturbances of the drive system during the driving process are identified. Then, the system parameter uncertainties and system performance indicators are presented on the Nichols chart as robust stability boundary, tracking boundary, and interference suppression boundary quantitatively, creating the condition for a comprehensive analysis of the drive system of the tracked vehicle.
Design of the system decoupling controller
Given that a strong coupling effect exists in the drive system, the design of a QFT controller for the drive system is difficult because the coupling effect makes it challenging to determine the boundaries of QFT parameter uncertainties. Besides, converting a multiple-input multiple-output (MIMO) system into a single-input single -output (SISO) system increases the difficulty in the design. Therefore, we decided to decouple the drive system before designing a QFT controller for it. For the R ≥ B/2 steering state of the tracked vehicle, a decoupling control matrix is designed to isolate the interference between the two inputs and the two outputs. The same decoupling method is applied to the 0 ≤ R < B/2 steering state. After the drive system model (equations (25) and (26) undergoes the Laplace transform, the decoupled drive system is obtained (see Figure 3)).

Block diagram of drive system decoupling.
The goal is to design a matrix to satisfy the following equation
Analysis of QFT controller design
The implementation of the above decoupling process allows achieving the mutual independence of the two motor channels. Now, we need to design a QFT controller for each motor channel to meet the performance robustness requirements of the drive system in the presence of parameter perturbation. Thus, a control model that aims at solving the parameter uncertainty and input perturbation uncertainty problems in the tracked vehicle drive system model is designed (Figure 4). A loop-shaping controller and a pre-filter are also designed for the controller.

Structure of robust control model based on QFT.
Design requirements and robust boundaries
Considering the actual operations of the tracked vehicle, three main system performance boundaries are needed, that is, robust stability margin, anti-interference performance, and servo tracking robust performance. The first thing that needs to be considered is that the bandwidth of the entire system should be determined by the power mechanism, and the bandwidth of the servo valve is much wider than that of the motor power mechanism. The natural frequency of the motor of the hydraulic drive system calculated using equation (32) is 18.2658 rad/s. Considering that the actual load includes the track and its transmission mechanism, the bandwidth of the actual power mechanism should be lower.
Robust stability margin
This indicator is used to ensure the stability margin of the system within the entire operating frequency band. Considering the general stability margin of the electro-hydraulic servo control system of the tracked vehicle, the amplitude margin is determined to be around 1.66 dB, and the phase margin is about 50°, which can meet the requirements of most working conditions. According to the empirical formula, the value of Mc is about 1.2
Anti-interference performance
Anti-interference performance indicators include output disturbance robustness and input disturbance robustness. From an analysis of the mathematical model of the output disturbance robustness, the external disturbance is attenuated at the rate of 40 dB per 10 octaves in the high-frequency band. So the transfer function is defined to rise at a rate of 40 dB per 10 octaves. Defining the robust performance boundaries can effectively suppress external disturbances lower than 3 Hz
The input disturbance robustness is used to ensure the robustness of the output under the influence of the system input interference signal. Here
Servo tracking robustness
By considering the actual operation of the tracked vehicle, the step response is used to constrain the servo performance
The step response performance indicators of the system within the frequency bandwidth are
Rise time: Tr =1 s
Overshoot: M = 10%
Adjustment time: Ts = 2 s
Design of controller
Generating a design template
Given the setting of the parameter uncertainty, the QFT design template shown in Figure 5 can be obtained. The frequencies ω = 0.1, 0.5, 1, 5, 10, and 30 rad/s are selected on the Nichols chart.

Template of drive system at selected frequencies on the Nichols chart.
Synthesis of performance boundaries
The overall design boundary of the QFT controller can be obtained by synthesizing the robust stability margin specification, robust output disturbance specification, and robust input disturbance specification. Gains, zero poles, and complex poles are added to the model to ensure (1) the low-frequency open-loop Nichols curve of the system is above and as close as possible to the constraint boundary of the response frequency point and (2) the amplitude–phase curve of the high-frequency section is outside the robust stability boundary. Figure 6 shows the resulting corrected system Nichols chart. The system performance is above the boundary at all frequency points, and the minimum gain is guaranteed. The resulting controller transfer function is

Open-loop frequency response and QFT boundaries.
Performance check
The results are fed into the transfer function model for verification to make the results more reliable. The stability margin robustness, output anti-disturbance robustness, and input anti-disturbance robustness are verified separately. The results are shown in Figures 7–9, respectively. The dotted lines in the graphs represent the performance indicators set in equations (33)–(37); the corrected system can meet the set performance requirements very well.

Stability margin robustness verification diagram.

Output anti-disturbance robustness verification diagram.

Input anti-disturbance robustness verification diagram.
Design of pre-filter
Considering that the system has the tracking capability, it is necessary to design a pre-filter for it. From the Bode diagram of Figure 10, the closed-loop frequency response boundary of the system is within the tracking performance boundary under the parameter disturbance condition. Thus, the transfer function of the pre-filter is

Design curve of pre-filter.
Simulation analysis
A simulation experiment on the motor channels was performed to verify the effectiveness of the proposed control strategy, the results were analyzed, and then the proposed control strategy was compared with the conventional PID control method. The necessary parameters can be calculated using equations (1)–(28). Table 2 lists the key parameters of the drive system.
Key parameters of the tracked vehicle drive system.
As shown in Figures 11 and 12, simulation analysis is performed for R ≥ B/2 and 0 ≤ R < B/2 steering states of tracked vehicles, and the comparison of response state and robust performance is obtained under the action of different controllers. In the simulation, an impulse signal with the size of 50 N m and the duration of 0.1 s was added to one side of motor passage at 5 s. It can be seen from the figure that compared with PID control, the change value of motor speed is small under the action of QFT controller, only about one-fifth of the change value under PID control, and the impulse signal has little impact on the motor speed of another passage; the change of motor torque value is little, which can recover in less time with less overshoot; the instantaneous turning radius change value of the tracked vehicle is only about one-fifth of that under PID control. The simulation results show that under the same parameter perturbation condition, QFT controller can suppress the perturbation in a short time and show a high robust performance.

R ≥ B/2 steering states.

0 ≤ R < B/2 steering states.
Experiment on steering performance of the tracked vehicle
A steering experiment was conducted to test the reliability of the robust control of the tracked vehicle drive system when the vehicle is performing R ≥ B/2 steering and 0 ≤ R < B/2 steering. In order to obtain good test results, the tests were carried out on a patch of brick-paved hard ground to increase the friction between the tracks and the ground, thus avoid slipping. The main instruments and tools used in the test include an encoder, tape measure, marking line, Jisibao G970 high-precision GNSS device, Bluetooth module, embedded board, and portable computer. During the test, the encoders mounted on the rear wheels collected the speeds of the two hydraulic motors, and the embedded board capable of Hampel filtering integrated the data before transferring the information to the computer. A vehicle-mounted portable computer-controlled the servo valve-controlled motor system. The tracked vehicle was positioned in real-time using the Jisibao G970 high-precision GNSS device to obtain its motion track during the test. A Bluetooth module transmitted the position data to the vehicle-mounted computer, and the computer processed the position data, obtaining the motion track of the vehicle. Figure 13 shows the test site.

The actual test site.
Analysis of motor output speed
The R ≥ B/2 steering test was carried out under the following conditions: the inner and outer motors were turning in the same direction at the speeds of 50 and 100 r/min, respectively, and the rotational speeds were 50 and 100 r/min, respectively. The 0 ≤ R < B/2 steering test was conducted under the following conditions: the inner and outer motors were turning in the opposite directions at the same speed of 50 r/min (the theoretical turn radius is 0). According to equations (7)–(10), it would take about 15 s for the tracked vehicle to make a 360° turn in the R ≥ B/2 steering test, and 8 s in the 0 ≤ R < B/2 steering test. When the tracked vehicle was turning, the E6B2-CWZ6C encoder collected the hydraulic motor speed signal that was then output by the embedded board (the signal processing frequency was set to 2 Hz).
The motor output speed is stable under the two steering conditions, and the speeds of the inner and outer motors are close to the theoretical values (Table 3). In the two steering tests, the average output speeds of the inner and outer motors are slightly lower than the set speeds due to the uncertain nature of the road surface resistance, but the discrepancies between the actual and the set values were small (the standard deviation was only 0.42) and the output speed was stable. This result demonstrates that the controller has a good robust performance in the presence of system parameter perturbation and uncertainty.
Speed parameters of the two hydraulic motors during steering test.
Analysis of steering track
In order to ensure that enough information points could be collected during the test, the sampling frequency of the G970 system was set to 5 Hz. Furthermore, the static data collection was performed under the idle condition to verify the positioning accuracy. The test results show that the test data points are all distributed in a circle with a diameter of 2 cm, indicating that the measurement accuracy is high enough to meet the test requirement. The steering position data were transmitted to the vehicle-mounted portable computer via the Bluetooth module and stored in the computer. Later, the latitude and longitude values were extracted from the track information and converted using the Gaussian projection calculation formula and sorted based on which the turning circle was obtained through the fitting with MATLAB. Figure 14 illustrates the installation location of the mobile station, where O and O1 are the geometric center of the ground-contact area of the tracked vehicle and the phase center of the mobile station, respectively. The distances between the phase center and the geometric centers a and b are 0.568 and 0.538 m, respectively, and the phase center is 1.6 m high from the ground. Equations (41) and (42) describe the relationship between the turning track radius of the geometric center of the ground-contact area and the turning track radius of the mobile station phase center.
where a is the longitudinal distance between the phase center of the mobile station and the geometric center of the ground-contact area of the tracked vehicle (a = 0.586 m), b is the transverse distance between the phase center of the mobile station and the geometric center of the ground-contact area of the tracked vehicle (b = 0.538 m), RR≥B/2 is the turning track radius of the geometric center of the ground-contact area when the tracked vehicle is performing R ≥ B/2 steering, RR≤B/2 is the turning track radius of the geometric center of the ground-contact area when the tracked vehicle is performing 0 ≤ R < B/2 steering, Rn1 is the turning track radius of the mobile station phase center when the tracked vehicle is performing R ≥ B/2 steering, and Rn2 is the turning track radius of the mobile station phase center when the tracked vehicle is performing 0 ≤ R < B/2 steering.

Installation of mobile station.
Figure 15 shows the turning tracks of the tracked vehicle; the vehicle slips and slides during in-situ steering process due to the high angular speed.

Turning tracks of the tracked vehicle: (a) R ≥ B/2 and (b) 0 ≤ R < B/2.
Table 4 lists the parameters of the turning track of the vehicle. The average turning track radii of the geometric centers of the inner and outer tracks are 1.828 and 0.033 m, respectively, and the relative errors are 1.56% and 3.19%, respectively. The smallness of the relative error demonstrates that the vehicle turns on the brick-paved hard ground with satisfactory performance, meeting the design requirements of the tracked vehicle.
Turning parameters of the tracked vehicle.
Conclusion
We have developed a QFT-based drive system controller based on an analysis of the structural characteristics of the high-clearance HTVs and the dual-pump dual-motor drive system. The proposed drive system controller is easy to implement, can improve the robust stability and dynamic tracking performance of the tracked vehicle drive system, and has good robust performance in the presence of system parameter perturbation and uncertainty. The simulation results show that the proposed control method barely generates oscillation and has a small overshoot when the control signal changes, outperforming the conventional PID control method. The results of the steering experiment show that the tracked vehicle can turn very well on the brick-paved hard ground, indicating that the controller achieves robust high stability by effectively suppressing the dynamic interference of the road surface. Facts have proved that QFT controller is feasible and stable for the hydraulic servo drive system of tracked vehicles, and it can be used as a reference basis for the research of control algorithms for hydraulic servo drive systems of the high-clearance HTVs.
Further work will be carried out in the following areas: the cross-coupling method is used to synchronize the two motor channels of the tracked vehicle, and the double closed-loop control method is used to further improve the control accuracy of the tracked vehicle during walking.
Footnotes
Handling Editor: Wei-Chiang Hong
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The study was supported by the Special project for the construction of Supported by China Agriculture Research System (grant no. CARS-03), Special project of modern agricultural industry technology system in Henan (grant no. S2019-02-G07), and Key Scientific Research Project of Henan Province Higher Education (grant no. 20A210029). We wish to acknowledge them for their support.
