Force Compensation Control for Electro-Hydraulic Servo System with Pump– Valve Compound Drive via QFT–DTOC

Each joint of a hydraulic-driven legged robot adopts a highly integrated hydraulic drive unit (HDU), which features a high power–weight ratio. However, most HDUs are throttling-valve-controlled cylinder systems, which exhibit high energy losses. By contrast, pump control systems offer a high efficiency. Nevertheless, their response ability is unsatisfactory. To fully utilize the advantages of pump and valve control systems, in this study, a new type of pump– valve compound drive system (PCDS) is designed, which can not only effectively reduce the energy loss, but can also ensure the response speed and response accuracy of the HDUs in robot joints to satisfy the performance requirements of robots. Herein, considering the force control requirements of energy conservation, high precision, and fast response of the robot joint HDU, a nonlinear mathematical model of the PCDS force control system is first introduced. In addition, pressure–flow nonlinearity, friction nonlinearity, load complexity and variability, and other factors affecting the system are considered, and a novel force control method based on quantitative feedback theory (QFT) and a disturbance torque observer (DTO) is designed, which is denoted as QFT–DTOC herein. This method improves the control accuracy and robustness of the force control system, reduces the effect of the disturbance torque on the control performance of the servo motor, and improves the overall force control performance of the system. Finally, experimental verification is performed using the PCDS performance test platform. The experimental results and quantitative data show that the QFT–DTOC proposed herein can significantly improve the force control performance of the PCDS. The relevant force control method can be used as a bottom-control method for the hydraulic servo system to provide a foundation for implementing the top-level trajectory planning of the robot.


Introduction
In recent years, mobile robots have been increasingly used in various industries and investigated extensively.Based on the motion mode, robots can be classified into wheeled, crawler, legged, serpentine, and spherical robots [1][2][3][4][5].In contrast with other types of robots, legged bionic robots offer noncontinuous support and good adaptability to unknown and complex environments.In recent years, bionic legged robots combined with hydraulic drives have been implemented to replace human beings to perform detection, transportation, rescue, military assistance, and other tasks in a complex environment.These robots feature a high power-toweight ratio, a large loading capacity, and fast response; they are widely used in military and civilian applications.
Currently, a highly integrated electro-hydraulic servo system is primarily adopted in the leg joint drivers of advanced hydraulic-driven legged robots.The typically used electro-hydraulic servo system is available two forms: the valve-controlled cylinder system [6][7][8] and pump-controlled cylinder system [9,10].The valve-controlled cylinder system is controlled using an electro-hydraulic servo valve.Meanwhile, for the pumpcontrolled cylinder system, the output flow and pressure of the pump are matched with the load requirements by changing the position or speed of the pump.The valvecontrolled cylinder system is a throttle system that features high energy loss.Moreover, the different pressure or flow requirements of each joint of the robot will cause a significant energy loss and reduce the endurance of the legged robots.The pump-controlled cylinder system is a direct-drive system.Although it features a higher energy efficiency than that of the valve-controlled cylinder system, its response speed and control accuracy are lower.Therefore, to overcome the shortcomings of the two systems when they are used alone and to combine their respective advantages, an electro-hydraulic servo system with fast response, high precision, and low energy loss must be designed.
Both force and position control are important bottomcontrol methods in robot control systems.In recent years, hydraulic and motor force control systems have been investigated extensively, and various control methods have been adopted, such as fuzzy logic control [11][12][13][14], sliding mode control [15,16], robust control [17,18], and intelligent control [19][20][21][22].The optimization of the control performance of force control systems has been demonstrated previously.However, in the practical application of robot systems, to ensure the stability of the body, algorithms of the core controller based on PID are still necessitated [23][24][25].In particular, an electrohydraulic servo system with a pump-valve compound drive, known as a pump-valve compound drive system (PCDS), features a fast response, high precision, and high energy utilization.However, when the system is disrupted by a complex and changeable external position, the high-performance force control requirements of the robot system will be difficult to satisfy.In addition, the hydraulic actuator of the system exhibits an asymmetric cylinder structure, and a few problems arise, such as the inconsistency of forward and negative dynamic characteristics; therefore, its control accuracy and robustness must be improved.
Based on a previous study as a foundation [26], the authors performed a mathematical modeling and parameter sensitivity analysis of the valve-controlled cylinder force control system and discovered the effects of the parameters in the system on the force control performance, thus providing a foundation for the optimization of the mathematical model.In addition, a force-control-performance optimization algorithm for force-based impedance control was investigated, and a compliance-eliminated controller was designed to improve the disturbance rejection performance of a valve-controlled cylinder force control system [27].Consequently, a preliminary basis was provided for mathematical modeling, system design, and control algorithm research.
This study is performed as follows: First, a nonlinear mathematical model of the PCDS force closed-loop control system is constructed.Second, QFT-DTOC, a forcecontrol method that combines quantitative feedback theory (QFT) and a disturbance torque observer (DTO), is designed to improve the control performance of the conventional PID controller.Finally, the effectiveness of the system performance and force control method is verified on the experimental platform of the PCDS, and the control effect is quantitatively analyzed.

System Composition and Structure
A hydraulic schematic diagram of the PCDS force control performance test platform is shown in Figure 1; it is primarily composed of a position closed-loop control system and a tested PCDS force control system.
To satisfy the control performance requirements of the PCDS force control while providing energy conservation, fast response, and high precision, the system should comprise three closed-loop control systems: the force control system of the pump-controlled loop, the pressure control system of the asymmetrical cylinder rodless cavity, and the pressure control system of the asymmetric cylinder rod cavity.A schematic illustration of PCDS force control is shown in Figure 2.

System Mathematical Modeling
To investigate the force compensation control method of the PCDS, an overall mathematical model of PCDS force control was established based on Figure 2.
In particular, a mathematical model of the key components of the PCDS, including the servo motor, gear pump, filling loop, and asymmetrical cylinder, was established.
A mathematical model of the servo motor was constructed, and the voltage balance equation of the servo motor is as follows: (1) where E is the back electromotive force, i the armature current, L the armature inductance, and R the armature winding resistance.
The reverse electromotive force of the servo motor is expressed as follows: where K c is the electromagnetic torque system, and w is the servo motor speed.
where K t is the torque coefficient of the servo motor, and T e the electromagnetic torque of the servo motor.
The servo motor and gear pump are directly connected, and the torque balance equation is as follows: where J is the conversion of the moment of inertia of the servo motor shaft, B p the viscous friction coefficient of the servo motor, and T L the external load torque of the servo motor.
To control the servo motor accurately, three types of closed-loop control are typically used: current closed loop, speed closed loop, and position closed loop from inside to outside.The torque control of the servo motor adopts a current closed loop, and the dynamic response is rapid.This study primarily focuses on the force control of the PCDS and considers the response of the system; thus, the torque control mode is adopted in the servo motor of the system.The current loop controller K PI can set the control parameters of the current loop to accurately control the servo motor torque.The transmission block diagram of the servo motor torque control system is shown in Figure 3.
In Figure 3, K V represents a voltage and current conversion system, where the effect of the external interference torque of the servo motor is disregarded.The open-loop transfer function between the input voltage and the output speed of the servo motor is expressed as follows: (2) .

Mathematical Model Construction of Gear Pump
When the gear pump is propagating forward, considering the internal and external leakage of the gear pump, the flow-pressure relationship is as shown in Figure 4.
The flow equations of the gear pumps are as follows: where D p is the displacement of the gear pump, w the speed of the gear pump, Q ip the internal leakage flow of the gear pump, Q ep1 the external leakage flow of the gear pump at point a, and Q ep2 the external leakage flow of the gear pump at point b.
Considering the nonlinear problems of leakage flow and pressure, the following expressions are derived: where K ip is the internal leakage coefficient of the gear pump, K im the external leakage coefficient of the gear pump, p 1 the pressure of the gear pump at point a, p 2 the pressure of the gear pump at point b, p 3 the pressure of the gear pump at point c (Pa), and 0 bar the pressure directly related to the external tank.

Mathematical Model Construction of Filling Loop
Mathematical models of a pressurized fuel tank and oneway valves were constructed in this study.
The PCDS is a closed system that requires oil to be replenished to the system through an oil-compensating (6) loop composed of a fuel tank, check valve, and servo valve.Its main functions are as follows: 1) Maintain a constant pressure in the low-pressure cavity of the asymmetric cylinder, 2) supply the hydraulic oil leakage of the force control system, and 3) prevent cavitation problems caused by high-frequency vibrations of the system.
In the modeling of a conventional pump control system, the oil replenishment link is typically disregarded, which results in a certain deviation in the analysis results.Therefore, the relationship between the flow rate and pressure in this study was analyzed, as shown in Figure 5.
The flow equation of each link is as follows: where Q cq1 is the check valve flow connected to the asymmetrical cylinder rodless cavity, Q cq2 the check valve flow connected to the asymmetrical cylinder rod cavity, Q sq1 the servo valve flow connected to the asymmetrical cylinder rodless cavity, Q sq1 the servo valve flow connected to the asymmetrical cylinder rod cavity, and Q tq the flow out from the pressurized fuel tank.
Considering the check valve as the throttle port model, the following expression can be obtained: (11)  where p tp is the output pressure of the pressurized fuel tank, and K cp is the flow pressure coefficient of the check valve.
Owing to the complexity of modeling pressurized fuel tanks and the fact that pressurized fuel tanks and energy accumulators exhibit similar functions, the mathematical model of an energy accumulator was used in this study instead of that of a pressurized fuel tank.The input flow and output pressure of the supercharged fuel tank are correlated as follows: where p gp is the initial pressure of the pressurized fuel tank; V gυ is the initial volume of gas in the pressurized fuel tank; k is the polytropic index of nitrogen, which is generally 1.0-1.4.

Derive Flow Equations of Servo Valve
The servo valves installed on the two sides of the asymmetric cylinder replenish or drain oil from the two chambers of the asymmetric cylinder.The transfer function of the servo valve can be simplified to that of a second-order oscillation system.The transfer function of the servo valve in terms of the input voltage and displacement of the valve core is as follows: where K a is the gain of the servo valve power amplifier, K xv the gain of the servo valve, ζ the damping ratio of the servo valve, and ω the natural frequency of the servo valve.
To simplify calculation, the slide valve level of the servo valve was set to that of an ideal zero-opening four-way slide valve, and a flow equation for the four-way slide valve was derived.The flow of the asymmetrical cylindrical rodless cavity is expressed as follows: The inlet flow rate of the servo valve connected to the rod cavity is expressed as where K d is the equivalent flow coefficient, x v1 the dis- placement of the servo valve connected to the rodless (15) cavity of the asymmetric cylinder, x v2 the displacement of the servo valve connected to the rod cavity of the asymmetric cylinder, p s the supply oil pressure of the servo valve, p 1 the pressure of the asymmetrical cylinder rodless cavity, p 2 the pressure of the asymmetrical cylinder rod cavity, and p 0 the return oil pressure of the system (0 bar).
The equivalent flow coefficient is expressed as follows: where C d is the flow coefficient of the servo valve throt- tle orifice, W the area gradient, and ρ the hydraulic oil density.

Mathematical Model Construction of Asymmetrical Cylinder
(1) Derivation of Flow Continuity Equation The inlet flow of the asymmetrical cylinder rodless cavity and the oil-in-cavity volume are expressed as follows: The inlet flow of the asymmetrical cylinder rod cavity and the oil-out cavity volume are expressed as follows: where A 1 is the effective area of the asymmetrical cylinder rodless cavity, A 2 the effective area of the asymmetrical cylinder rod cavity, x p the piston displacement of the asymmetrical cylinder, C im the internal leakage coefficient of the asymmetrical cylinder, β ε the elastic modu- lus of the effective volume, V 01 the initial volume of the asymmetrical cylinder rodless cavity, and V 02 the initial volume of the asymmetrical cylinder rod cavity.
The initial volumes of the rodless cavity and rod cavity of the asymmetric cylinder are expressed as follows, considering the volume variation of the operating cavity, connecting pipe, and asymmetric cylinder of the gear pump: where V p1 is the volume of the operating cavity of the gear pump connected to an asymmetrical cylinder rodless cavity, V p2 the volume of the operating cavity of the gear pump connected to the asymmetrical cylinder rod (20) , cavity, V g1 the pipe volume connected to the gear pump and the asymmetric cylinder rodless cavity, V g2 the pipe volume connected to the gear pump and the asymmetric cylinder rod cavity, L the total piston displacement of the asymmetrical cylinder, and L 0 the initial position of the asymmetrical cylinder.
(2) Derivation of Force Balance Equation Because the asymmetric cylinder is affected by the inertial force, viscous damping force, elastic force, and any external load force, the balance equations of the output and load forces of the asymmetric cylinder are derived as follows: where m t is the total mass converted to the asymmetric cylinder piston, including the load mass block, piston, connection pipe, and oil in the asymmetric cylinder; K is the load stiffness of the asymmetric cylinder; B p is the damping coefficient of the load and asymmetric cylinder; F f is the Coulomb friction of the load and asymmetric cylinder; X L is the arbitrary load location on the asym- metric cylinder piston.

Mathematical Model Construction of Pressure Sensor
The sampling frequency of the selected pressure sensor is more than five times that of the control system.The transfer function of the sensor is equivalent to a proportional link, and the transfer function between the feedback voltage and oil pressure is expressed as follows: (24) where K ps is the gain of the pressure sensor, and P the pressure signal detected by the pressure sensor.By combining Eqs. ( 1)-( 25), the entire mathematical model of the PCDS force control can be established, in which a PID controller is adopted in the pump control and valve control loops, as shown in Figure 6.

Simplification of Transfer Function of Pump-controlled Asymmetric Cylinder
To facilitate the design of the pump control loop controller, the transfer function of the pump control loop was simplified without considering the effects of valve control and oil replenishment on the system.Additionally, the low-pressure cavity of the asymmetric cylinder was assumed to be zero.The PCDS force-control transfer function is relatively complex and difficult to simplify.In addition, the difference between the forward and backward dynamic characteristics of the asymmetric cylinder was considered.Therefore, to facilitate the modeling, the asymmetric cylinder was modeled separately based on the positive and negative motions of the piston rod, and different proportional gains were used to reduce or eliminate the difference.In the positive motion of the piston rod of the asymmetric cylinder, the gear pump absorbs hydraulic oil from the rod cavity to output high pressure oil to the rodless cavity; its positive and negative motions are illustrated in Figure 7.
When the pump-controlled asymmetric cylinder is undergoing positive motion, its transfer function is derived as follows.The flow of the gear-pump high-pressure cavity is expressed as (25)  The flow continuity equation for the rodless cavity of the asymmetric cylinder is Assuming that the low-pressure cavity of the asymmetric cylinder is zero, the external load and friction forces are disregarded, and the force balance equation is simplified as follows: Combining Eqs. ( 26)-( 28), the open-loop transfer function of the pump-controlled asymmetric cylinder system with positive motion is expressed as follows: Similarly, the open-loop transfer function of the pump-controlled asymmetric cylinder system in negative motion is expressed as follows: (26) .
By comparing Eqs. ( 29) and ( 30), the unified expression for the open-loop transfer function of the pumpcontrolled asymmetric cylinder can be obtained as follows: where, After simplifying and integrating Eq. ( 31), the simplified transfer function of the pump-controlled asymmetric cylinder is expressed as follows: where ω m is the natural frequency of the load; ξ m is the damping ratio of the load; ω r is the ratio of the series cou- pling stiffness to the damping coefficient of the hydraulic spring and load spring; ω 0 is the natural frequency of the system composed of a hydraulic spring, load spring, and mass; ξ 0 is the damping ratio.
The parameters in Eq. ( 32) are expressed as follows: where K h is the hydraulic spring stiffness.
( A control block diagram of the entire force control system after processing is shown in Figure 8.

Design of PCDS Force Controller
In this section, a force controller (QFT-DLOC) based on QFT and a DTO is proposed for the force control system of the pump control loop to improve the force control performance of the system.

QFT Controller
A QFT controller [28][29][30] is a robust controller based on the frequency domain.When the structure of a control object is uncertain and disturbance is present, a highly robust controller is to be designed based on the frequency domain theory of the control system.The design structure of the QFT controller is shown in Figure 9.
In Figure 9, r(t) is the system input, and y(t) is the system output, both of which can be measured independently.G(s) is the controlled object exhibiting uncertainty, d 1 (s) and d 2 (s) are the unknown external dis- turbances of the system, and P(s) is the controller, whose role is to ensure the robust performance of the system.F (s) is a prefilter that guarantees the tracking perfor- mance of the system.By designing the appropriate transfer functions P(s) and F (s) , the control system can satisfy the design performance index, and the bandwidth of the QFT controller can be minimized.The open-loop transfer function L(s) and closed-loop transfer function T R (s) are subsequently defined for the control system structure shown in Figure 9.
Based on the structure of the control system shown in Figure 9, the controller P(s) and prefilter F (s) are designed to stabilize the closed-loop system and satisfy the frequency domain index.The typical frequencydomain indices are as follows: (1) Index of robust stability

Determination of Parameter Uncertainty Set and System Template
The QFT controller was designed based on the transfer block diagram of the pump-controlled asymmetric cylinder force-control system shown in Figure 5, i.e., without considering the disturbance torque of the servo motor.
The output force of the pump control loop F p and the open-loop transfer function of the desired input F r are expressed as follows: where k a is the equivalent proportionality coefficient.By calculating the PCDS parameters, the formula of the open-loop transfer function (Eq.( 46)) and the initial values (as shown in Table 1) are obtained.The equivalent proportionality coefficient is a typical control parameter.To improve the robustness of the system, the equivalent proportional coefficient must be adjusted dynamically based on the different operating conditions.In the sampling system, because of the nonlinear and high-order characteristics of the servo valve, after the servo valve becomes equivalent to the second-order oscillation link, it cannot fully reflect its actual dynamic characteristics.Consequently, the natural frequency and damping ratio of the servo valve are not constant because of the input signal and operating parameters.Therefore, ±7% of each parameter was specified as the uncertainty of the openloop transfer function for force control, and the variation range of each parameter is shown in Table 1.
In the QFT controller, the object template [31,32] is adopted to describe the uncertainty of objects at different frequency points.To obtain an appropriate object template while considering the frequency response of the system, a set of frequency points, i.e., ω = {0.01,0.05, 0.25, 1, 2, 5, 10, 20, 40}, was selected.At each frequency point, the amplitude range and the phase (46) angle of the controlled object is displayed on the Nichols diagram, which forms the object template.The object template of the corresponding frequency points of the system was established using the quantitative feedback toolbox in MATLAB, as shown in Figure 10.Each point in Figure 10 represents the amplitude-frequency and phase-frequency characteristics of the system at different frequencies.Uncertain variations exist in the parameters, although the range is small.

System Index Requirement and Solution
In designing the QFT controller, the closed-loop robust stability and system tracking performance were regarded as the performance indices. (

1) Robust System Stability
When designing the controller in the frequency domain, the steady-state characteristics of the system can be determined based on the amplitude and phase  To ensure the stability of the controlled system and that the resonance peak of the controlled system does not exceed 1.2, the amplitude and phase margins were as G M = 5.27 dB and φ M = 49.25,respectively, based on cal- culations using Eqs.(47) and (48), respectively.
Based on the closed-loop robust stability constraint specified in Eq. ( 45), the robust stability boundary of the nominally controlled object of the pump control loop can be obtained using the QFT robust control toolbox.
As shown in Figure 11, when designing the QFT controller, provided that the Nichols diagram of the open-loop transfer function of the force control system is outside the robust stability boundary, the system can yield a robust stability index. (

2) Index of System Tracking Performance
To facilitate the design of the upper and lower bound functions of the system tracking performance, a typical second-order link was selected as the basic form of the system tracking performance, and the closed-loop transfer function of the second-order system is (47 To stabilize the bandwidth of the PCDS to 3-10 Hz, the boundary functions were obtained when the natural frequencies were 3.18 and 9.55 Hz, and damping ratios were 0.9 and 0.6.The boundary functions were designed as follows: The boundary functions shown in Eqs. ( 49) and (50) were adopted to obtain the trajectory tracking boundary of the nominally controlled object in the pump control loop through the QFT robust control toolbox, as shown in Figure 12.
As shown in Figure 12, when designing the QFT controller, provided that the point of the open-loop transfer function of the force control system on the Nichols diagram is outside (for the closed tracking performance boundary) or above (for the open tracking performance constraint boundary) the tracking performance boundary of each operating frequency, the system can satisfy the tracking performance index.
The two performance index requirements above were converted into a series of boundaries that constrain the open-loop frequency response curve in the Nichols diagram.The different constraint boundaries at various frequency points were integrated to form comprehensive performance boundaries, as shown in Figure 13. (49) (50) T u (s) = 400 s 2 + 36s + 400 , (51) T 1 (s) = 3600 s 2 + 72s + 3600 .As shown in Figure 13, the open-loop Nichols curve of the system comprises a frequency range inside (for the closed constraint boundary) or below (for the open constraint boundary) the comprehensive performance constraint boundary area.
To satisfy the two performance indexes, the open-loop Nichols curve of the system was adjusted to outside or above the comprehensive performance constraint boundary area by adding zeros and poles as well as adjusting the gain, which is analogous to designing the controller's P(s).
The design process is performed such that each frequency point is the closest possible to the comprehensive performance boundary, and the high-frequency gain decreases rapidly to suppress the effects of noise and the high-frequency resonance of the system at high frequencies.After drawing the composite index boundary on the Nichols diagram, the open-loop frequency response curve of the nominal object was plotted on the Nichols diagram using loop shaping.By continuously adjusting the zeros, poles, and controller gain, the response curve can be adjusted to above the performance index boundary at a low frequency, where its closer proximity to the boundary is desired.At this time, the bandwidth of the controller is minimum and does not intersect the stable boundary at high frequencies.When the response curve satisfies the boundary requirements, the zeros, poles, and gains used for adjustment constitute the designed controller.All characteristics are reflected in the Nichols diagram, and the process of adjusting the curve is equivalent to the process of obtaining the transfer function.The appropriate P(s) is selected through simulation analysis, and the P(s) of the QFT controller can be obtained as follows: Figure 14 shows the open-loop frequency response of the system after calibration.As shown, the frequency range of the system open-loop Nichols curve is inside or below the comprehensive performance constraint boundary region.

Design of QFT-DTO Controller
In a PCDS, the disturbance torque generated by the load pressure of the gear pump affects the torque control performance of the servo motor, and the load pressure of the gear pump tends to fluctuate.In particular, when the system is disrupted by external locations, the pressure fluctuation is more intense, which adversely affects the force control accuracy and response speed of the system.When designing the QFT controller described in Section 2.1, the disturbance torque of the servo motor is disregarded to reduce the design difficulty.However, the load torque of the servo motor in the actual system is relatively high, particularly when controlling the HDU of a high-precision robot joint; as such, it should not be disregarded.Therefore, the disturbance torque was introduced into the force control system, and through appropriate control compensation, the disturbance rejection ability of the servo motor and the control precision of the force control system were improved.
The disturbance torque of the servo motor is a function of the gear pump's connection shaft angle and the load pressure of the gear pump.Additionally, it is affected Figure 13 Open-loop frequency response of system before calibration.

Figure 14
Open-loop frequency response of system after calibration by factors such as friction resistance, which cannot be obtained using simple calculation theory.Furthermore, no experimental conditions are available for detecting the disturbance torque.Therefore, the disturbance torque of the servo motor must be measured indirectly via the DTO controller to design the appropriate compensation control experiments.

Design of Servo Motor DTO Controller
To obtain the disturbance torque of the servo motor more accurately, the disturbance torque of the servo motor was estimated using the DTO controller, and compensation was performed to eliminate the effect of the external torque disturbance on the control performance of the servo motor.The transfer function in Figure 8 was simplified to the product form of the basic transfer function link.The transfer block diagram of the servo motor torque control system is shown in Figure 15.
To simplify the transfer function of the servo motor, the transfer block diagram of the torque control system of the servo motor was deformed, and the transfer block diagram of the torque control system of the servo motor after deformation was obtained, as shown in Figure 16.
A simplification of the transfer block diagram shown in Figure 16 yields the block diagram shown in Figure 17.
In Figure 17, G 1 (s) is the transfer function of the servo motor from the control voltage input to the torque output, and G 2 (s) is the transfer function of the servo motor from the torque input to the speed output.G 1 (s) and G 2 (s) are expressed as follows: When designing the QFT controller, the disturbance torque of the servo motor was not considered.By disregarding the effect of the external disturbance torque, the open-loop transfer function of the servo motor can be obtained as follows: The open-loop transfer function of the servo motor was simplified to obtain the following expression: where K sm is the open-loop gain of the servo motor (zero dimension), ω sm the natural frequency of the servo motor, and ζ sm the damping ratio of the servo motor (zero dimension).
Each parameter is expressed as follows: Based on the simplified torque control transfer block diagram of the servo motor, the DTO control block diagram of the servo motor can be obtained, as shown in Figure 18, where T L is the disturbance torque observer.
In Figure 18, G 1 (s) is the transfer function from the control voltage input to the torque output of the servo motor, and G 2 (s) is the transfer function from the torque input to the speed output of the servo motor.Their expressions are shown in Eqs. ( 53) and (54) above.
Using the principle of disturbance invariance, the following were deduced: Based on Figure 18, the estimated value of the equivalent load torque can be obtained as follows: Therefore, the load torque of the servo motor can be accurately estimated using the method above.
Owing to the noise of the measured signal in the actual hydraulic control and the inability to establish an accurate mathematical model of the servo motor, the accuracy of the load torque observation decreased.Therefore, (56) a low-pass filter, G L (s) was introduced after the observed value.And the DTOC compensation control block diagram of service motor load can be obtauned, as shown in Figure 19.
In Figure 19, G L (s) is critical to the design of the DTO, and a positive G L (s)/G 2 (s) must be ensured, i.e., the relative order of G L (s) should not be less than that of G 2 (s).Additionally, the bandwidth design of G L (s) should comprehensively account for the robust stability and disturbance rejection suppression ability of the load observer.
As G 2 (s)=1/Js, a low-pass filter with a thirdorder denominator and a first-order numerator was adopted, namely, N = 3, M = 1, and k = 0.1.Based on can be obtained.The low-pass filter is expressed as follows: where τ = 0.001.

Design of QFT-DTO Controller
Combining the QFT controller described in Section 2.1 and the DTO controller described in this section, the QFT-DTO control block diagram of the pump control loop was obtained, as shown in Figure 20.
The model parameters and initial values used in the simulation of the pump-valve composite driving force control system are listed in Table 2.

Operating Principle of PCDS Experimental Platform
To investigate the performance of the PCDS force control, a PCDS performance test platform was designed, as shown schematically in Figure 21.
The test platform for evaluating the PCDS forcecontrol performance is shown in Figure 21.A counter (61) cylinder was adopted in the experimental platform.
The section shown on the right is the PCDS, which was used to perform an experiment to evaluate the force control system of the HDU of the robot joint.The section shown on the left is the valve-controlled cylinder position control system, which was used to simulate the disturbance load during its movement.However, in the present experiment using the PCDS force control performance test platform, the servo valve pressure port P was connected to the volume of the two cavities of the asymmetric cylinder in the force control system, and the servo   Operating chamber volume of gear pump V p1 valve could only be used for oil discharge.The system was operated in the pump-drain valve compound drive mode.

Experimental Plan
The effectiveness of the QFT-DLO control method was evaluated.The valve control loop was controlled using the same PID parameters, and the specific experimental conditions are listed in Table 3.   3, and 5, when the force control system is provided with a sinusoidal force input and the disturbance position is based on the sinusoidal input, the QFT and DTO can improve the forcefollowing accuracy.In particular, in the case involving load position disturbance, the DTO compensation effect was better.When the sinusoidal force-following signal of the system was 1000 N with a frequency of 0.5 Hz and the disturbance position featured an 8 mm amplitude with 1 Hz frequency, the maximum force error using PID control was 444.11N, and the maximum force error of the force control system using QFT was 330.28 N. The maximum reduction rate of force error was approximately 25.63%.The maximum force error using the QFT-DTOC was 243.19 N, and the maximum force error elimination rate was approximately 45.24%, which was 19.61% less than the maximum force error rate using the QFT system only.

Experimental Result
When only QFT was adopted in the system, the response curve of the force control system featured a certain phase angle lag, which is attributed to the following reason: Because the control method is designed based on the amplitude-frequency characteristic curve and the corresponding requirements for the phase angle is not satisfied when the QFT is adopted, a certain degree of phase angle lag is inevitable.However, the phase angle lag in the force control curve improved after using the QFT-DTOC.This is because by compensating for the disturbance load of the servo motor, the disturbance rejection ability of the servo motor and its response ability are improved.Hence, the response of the pump control loop is improved, and the issue regarding the phase angle lag is mitigated.
In general, under different experimental conditions, the QFT control method can improve the control accuracy by more than 20% compared with PID control, and

Conclusions
To satisfy the requirements of energy conservation, high precision, and fast response of the joint HDU of a legged robot combined with a pump-valve compound drive system, a mathematical model was established for each hydraulic component.Subsequently, the force control method was analyzed, and a force control method combining QFT and a DTO was devised.Before incorporating the QFT-DTO controller into the system, i.e., when only a simple PID controller was used, the force input signal was followed successfully, although the control performance was unsatisfactory.When the amplitude of the sine signal was increased, the system exhibited a significant following error.When only QFT was adopted in the system, the control accuracy was improved by more than 20% compared with that afforded by PID control.However, a phase angle lag was observed in the response curve of the force control system because QFT is a control method that is based on the amplitude-frequency characteristic curve, and it does not satisfy the corresponding requirements for the phase angle.Thus, a certain degree of phase-angle lag is inevitable.After incorporating the QFT-DTOC compensation controller designed in this study, the robustness and control accuracy of the force control system improved, and the following accuracy increased by more than 41%.
The PCDS test and experimental platform used in this study comprised an internal gear pump, which featured low-speed crawling, weak flow pressure pulsation, and other disadvantages, which resulted in significant fluctuations in the force control curve.Therefore, the lowspeed characteristics of the internal gear pump should be investigated in the future such that relevant control compensation methods can be devised to reduce its effect on force control performance.

Figure 4 Figure 5
Figure 4 Diagram showing flow-pressure relationship of gear pump

Figure 6
Figure 6 Overall mathematical model of PCDS force control

Figure 7
Figure 7 Diagrams illustrating positive and negative motions of pump-controlled asymmetric cylinder

Figure 8 Figure 9
Figure 8 Transfer block diagram of pump-controlled asymmetric cylinder force control system

Figure 10
Figure 10 Object templates at different frequency points

Figure 11 Figure 12
Figure 11 Robust stability boundary of pump control loop force control system

Figure 15 Figure 16
Figure 15 Transfer block diagram of the servo motor torque control system

Figure 18
Figure 18 DTO control block diagram of servo motor Figure 19 DTOC compensation control block diagram of service motor load

Figure 20
Figure 20 Control block diagram of QFT-DLO , 24, 25, where PID in the figure indicates the output force of the pump control loop with a simple PID controller, QFT indicates the output force of the pump control loop adopting the force control method based on QFT, and QFT-DTOC indicates the output force of the pump control loop adopting the force control method based on QFT and the DTO.To facilitate the analysis, Tables 4 and 5 are constructed based on the evaluation indices.As shown in Figures 24, 25, Tables

Figure 21
Figure 21 Test platform for evaluating PCDS performance

Figure 22 Figure 23 Figure 24
Figure 22 Force input exhibits 1000 N amplitude with 0.5 Hz frequency and no disturbance position

Figure 25
Figure 25 Force input exhibits 2000 N amplitude with 0.5 Hz frequency, and disturbance position exhibits 8 mm amplitude with 1 Hz frequency

Table 1
Parameter variation range of force control simulation model

Table 2
Simulation model parameters and initial values of pump-valve composite driving force control system

Table 4
lists the o23put force of the pump control loop inFigures 22,23

Table 3
Experimental conditions of QFT-DTOC

Table 4
One of QFT-DTOC performance indicators (%, indicated in two decimal digits)