Parameter Identification and Control Algorithm of Electrohydraulic Servo System for Robotic Excavator Based on Improved Hammerstein Model

In view of the nonlinearity and time-varying characteristics of the electrohydraulic servo system of the robotic excavator, a nonlinear adaptive identification and control algorithm based on improved Hammerstein model is proposed. +e Hammerstein algorithm model can approximate the nonlinear system with enough precision, but for the time-varying systems is not satisfactory. In order to compensate for the influence of time-varying factors, the fuzzy control module is designed to adaptively update the forgetting factor. +e experimental results show that the improved Hammerstein model error is about 40.11% less than the classical Hammerstein model error. +is proves that the improved Hammerstein model is feasible and effective to describe the electrohydraulic servo system of the robotic excavator.


Introduction
In order to continuously improve the operation efficiency and operation quality of excavator, and to expand the application areas of the excavator, autonomous operation has become an important research direction [1,2]. A key problem of the excavator's autonomous operation is how to establish the scientific effective model of the electrohydraulic servo system and improve the trajectory tracking accuracy of the robotic excavator [3,4]. Due to the existence of dead zone, saturation, nonlinear friction, and the nonsymmetry of hydraulic cylinder, the electrohydraulic servo system has strong nonlinear and time-varying characteristics [5][6][7][8]. Traditionally, the mechanism modeling method [9] and the linear identification algorithm [10] are used to establish the mathematical model of the nonlinear system. In the modeling process, the mechanism modeling method makes some assumptions and simplifies the system, resulting in a certain error between the established model and the actual system. e linear model identification algorithm establishes the model of the approximation system according to the input and output data, but the model cannot describe the nonlinear characteristics such as friction and dead zone. Many researchers have done a lot of works in this field. Tri et al. [11] proposed and designed an adaptive control algorithm based on a modified backstepping algorithm with the iterative learning scheme for the trajectory control of an electrohydraulic actuator; experimental results showed that the proposed controller can provide excellent tracking response. Yung et al. [12] designed a controller based on a proportional feedback, a disturbance compensator, and a relay controller for a hydraulic cylinder in mobile hydraulics applications; the controller shows its robustness to different lowering velocities by comparing with a PI controller. Guo et al. [13] carried out a study on the nonlinear cascade trajectory tracking control for an electrohydraulic system using the sliding mode control and the backstepping technique, and the stability of the controller was proved based on Lyapunov theory. Both the simulation and experimental results proved the excellent tracking performance of proposed approaches. Boaventura et al. [14] conducted an investigation on impedance control for a hydraulic robot actuator, many relevant aspects regarding the control algorithm were presented, and the experiment showed that tracking capabilities of the system were improved. Yao et al. [15] carried out a study on high-accuracy tracking control for hydraulic actuators, an improved LuGre friction model was derived, and then, an adaptive backstepping controller was proposed. Experimental results showed the controller can not only achieve the excellent asymptotic tracking performance but also ensure the robustness. Besides, the hybrid synchronization problem of chaotic systems was also investigated [16][17][18][19]. From all the research studies above, we can find that many valuable results have been obtained in the research field of hydraulic servo system. However, existing studies focused on the nonlinear characteristics, and little research on the timevarying characteristics has been done.
According to the strong nonlinear and time-varying characteristics of the electrohydraulic servo system, an adaptive identification algorithm based on the Hammerstein model is proposed. Hammerstein model [20,21] consists of static nonlinear module and dynamic linear module. Its linear module is described by discrete autoregressive model, and nonlinear module is described by piecewise polynomial basis function.
e Hammerstein model can effectively describe the nonlinear system, and the output error is significantly smaller than the linear identification algorithm. However, when the algorithm is used for nonlinear systems with time-varying parameters, the identification results are not satisfactory. erefore, on the basis of the Hammerstein nonlinear identification algorithm, the limited memory method is introduced, fuzzy control module is designed, and the forgetting factor is adaptively adjusted by the fuzzy algorithm. e identified model can be updated in real time to better approximate the actual system. e rest of the paper is organized as follows. In Section 2, the mathematical model of the electrohydraulic servo system is established by using the mechanism modeling method, and the nonlinearity and time-varying characteristics of the system are analyzed. In Section 3, the Hammerstein model structure is presented, and the identification algorithm of the electrohydraulic servo system based on the Hammerstein algorithm is developed. In Section 4, the forgetting factor is introduced, and the fuzzy control model is established to tune the forgetting factor in real time. In Section 5, the experimental setup and experimental procedure are explained, and the experimental results are discussed. In Section 6, conclusions are presented in the final section.

Model Analysis of Electrohydraulic
Servo System e electrohydraulic servo system of the robotic excavator is mainly composed of pilot electrohydraulic proportional valve, valve-controlled asymmetric hydraulic cylinder, sensor, and control system, etc. e system model can be represented by the dynamic equation of electrohydraulic proportional valve, flow equation of spool valve, the flow continuity equation of hydraulic cylinder, the hydraulic cylinder, and the load force balance equation. e electrohydraulic proportional valve is a conversion unit between the electric signal and the hydraulic output, and the corresponding pressure output is related to the input signal, which is used to move the spool of the multiway directional valve. is can be expressed by the first-order linear differential equations: where x v is the main spool displacement, τ v is the time constant, K 1 is the gain of the electrohydraulic proportional valve, and u is the system input. e movement of the robotic excavator's working device can be simplified as the valve control asymmetric hydraulic cylinder system. e pressure fluid and load can be seen as a mass-spring-damped oscillation system. e system schematic diagram is shown in Figure 1. A 1 and A 2 are the effective areas of the cylinder chambers; p 1 and p 2 are the pressures in cylinder chambers; q 1 and q 2 are inlet flow and return flow; p s and p 0 are system pressure and return pressure; y is the cylinder piston displacement; m is the equivalent mass of the piston and workload; B p is the viscous damping coefficient; K is the load elasticity coefficient; F L is the disturbance force.
Take the forward movement of the hydraulic cylinder as an example to establish the mathematic model of the valvecontrolled cylinder system. e load pressure equation of the hydraulic cylinder in steady state is as follows: where p L is the load pressure; n � A 2 /A 1 is the effective area ratio. e flow of the directional valve is related to the pressure of inlet and outlet and the speed of the piston. e flow calculation formula of the directional valve is as follows: where C d is the flow coefficient, ω is the area gradient of the valve, and ρ is the liquid density. According to equations (2) and (3), the load flow equation of the spool valve can be obtained: Its linearized expression is as follows: where K q is the flow gain and K c is the pressure-flow gain.
Ignoring the influence of the pipeline on the hydraulic system, and assuming that the pressure in the working chamber of the hydraulic cylinder is equal everywhere, the flow continuity equation of the hydraulic cylinder is as follows [22]: 01 and V 02 are the initial volumes of the hydraulic cylinder; C ip and C ep are internal leakage and external leakage coefficients; β e is liquid bulk elastic modulus. At this point, the calculation formula of load pressure is as follows: From equation (3) and equation (6), the load pressure equation can be obtained: where V t is the equivalent volume, V t � 4V 1 /1 + n 3 , and e load force of the hydraulic cylinder includes inertial force, viscous damping force, elastic force, and random load force. According to Newton's second law, the force equilibrium equation can be obtained: e piston displacement of the hydraulic cylinder under the action of the system input u and the external load force F L can be obtained based on the above equations: It can be seen from the transfer function of the electrohydraulic servo system that some parameters are difficult to be calculated. Due to the nonlinear factors such as dead zone, saturation, and friction, the system has strong nonlinear and time-varying characteristics; at the same time, the linear terms are cross-coupled with the nonlinear terms. If the parameters of the electrohydraulic servo system are identified directly, several state variables need to be measured, in addition to increasing the test cost; the measurement noise is also introduced to produce larger identification error. erefore, it is necessary to design an identification algorithm which can accurately describe the nonlinear and time-varying characteristics of the system.

Nonlinear Identification Algorithm
Based on Hammerstein e Hammerstein algorithm model describes the nonlinear system by the combination of static nonlinear module and dynamic linear module, and its system structure diagram is shown in Figure 2, where u(k), x(k), v(k), and y(k) are identification input, intermediate state, noise sequence, and system output, respectively. N(·) and L(z) are nonlinear module and linear module. e linear and nonlinear characteristics are distinguished by the intermediate state x(k), which can be used to decouple the linear and nonlinear models. e Hammerstein algorithm only requires the input and output data, which can effectively reduce the measurement noise and improve the identification precision [23].
Due to the asymmetric characteristics of friction in the electrohydraulic servo system, the nonlinear part of the Hammerstein algorithm is represented by piecewise polynomial function [24], and the input and output relationship of the nonlinear part is as follows: Figure 1: Schematic diagram of servo hydraulic cylinder system.
where f k and g k are nonlinear coefficients and r is the order of the polynomial. e linear part of the Hammerstein model is represented by a discrete autoregressive model, namely: where ; n k is system delay; a i (i � 1, 2, · · · , n a ) and b j (j � 0, 1, · · · , n b ) are linear coefficients. Take the case of u(k) > 0. Substituting equation (11) into equation (12), the least square format based on the Hammerstein model can be simplified as follows: e autoregressive variable ϕ(k) and the estimated parameter vector θ are as follows: where C ij � f i b j . e estimated value of θ can be obtained by the least square method: where Using the normalization method, and let f 1 � 1. e following equation can be obtained: After matrix transformation, the expression of f k is as follows: j � 0, 1, . . . , n b ; k � 0, 1, . . . , r . (17) According to equation (16), the intermediate variable x(k) can be calculated in the iterative process.

Adaptive Dynamic Identification Based on Fuzzy Algorithm
Due to the time-varying characteristics of the electrohydraulic servo system, the problem of model mismatch may occur over the long run, which requires the system model has to be frequently updated. At the same time, it is necessary to consider the phenomenon of "data saturation" in the update process, that is, the impact of new data on parameter estimation becomes weaker and weaker over time.
When the system is highly variable, the Hammerstein model will not be able to quickly track the changes in the system, causing the identification to fail. erefore, the "forgetting factor" is introduced to solve the problem of data saturation. According to the limited memory recursive least square algorithm [25], the following equation can be obtained: where P(k) � (Φ T Φ) − 1 . e forgetting factor λ is the weight of the old data in the identification process. e bigger the forgetting factor λ, the smaller the influence of the new data on the identification results. e smaller the forgetting factor λ, the greater the impact of the new data on the identification results. erefore, in order to solve the problem of model mismatch, when the system changes slowly, it can adopt the larger forgetting factor λ to reduce the influence of the sensor's interference error on the model; when the system changes rapidly, it should adopt the smaller forgetting factor λ to make the model match the system more quickly.
In order to adapt to this change and obtain the optimal forgetting factor in real time, a fuzzy controller [26,27] is designed to update the forgetting factor. e input of the controller is the displacement error e and the error variation e linear function can adjust the factors to reduce the error rapidly. e curve function is relatively smooth and is conducive to the stability of control. e membership function of fuzzy controller is designed based on the advantages of above two functions. When the error is large, the speed of adjustment is mainly considered and the trigonometric membership function is adopted; when the error is small, the stationarity is considered mainly and the curve membership function is adopted. e curve membership function is shown in Figure 3, and the corresponding fuzzy control rules are shown in Table 1.
e methods of defuzzification include weighted average method, maximum membership function method, and gravity method. According to the advantages and disadvantages of these methods, the gravity method is chosen for defuzzification, its calculation formula is as follows: where λ * is the updated value of forgetting factor λ. e overall control structure and data flow of the electrohydraulic servo system of the robotic excavator are shown in Figure 4. In the figure, y r is the desired displacement and y is the actual displacement.
In summary, the steps of adaptive identification and control for robotic excavator based on Hammerstein are as follows: Step 1: according to the mechanism characteristics of the electrohydraulic servo system of the robotic excavator, set the values of n a , n b , and n k Step 2: set the values of θ(0), P(0), and forgetting factor λ Step 3: offline identification; calculate the nonlinear module coefficients of the Hammerstein model Step 4: collect the experimental data of the electrohydraulic servo system of the robotic excavator Step 5: update the forgetting factor by fuzzy control algorithm Step 6: update the linear module coefficients of the Hammerstein model Step 7: calculate the next output u(k + 1) of the controller Step 8: return to Step 4 and continue the loop

Experimental Analysis
e experimental prototype of the robotic excavator is developed from the existing backhoe mountain excavator, which is composed of backhoe mountain excavator, sensor control system, and the 3D laser radar. e working device of the robotic excavator is composed of boom, arm, telescopic arm, and bucket in series, as shown in Figure 5. Each joint is driven by servo hydraulic cylinder, and the displacement of the piston rod is measured by using the cable type of absolute encoder, as shown in Figure 6. e structure of electrohydraulic proportional position control system of the robotic excavator's working device is the same, which is mainly composed of the electrohydraulic pilot proportional pressure reducing valve, LUDV multiway valve, hydraulic cylinder, and the cable type absolute encoder used to measure the displacement of the hydraulic cylinder piston rod, as shown in Figure 6. e DSP controller is connected with main control computer and encoder through CAN bus. e main control computer gives the target position of the piston rod of the hydraulic cylinder, and the controller collects the displacement of the piston rod as the feedback, takes the difference between the target position and the feedback as the input of the controller, and then sends the control amount to the driver. e driver generates the PWM signal PWM and changes the voltage of the electrohydraulic pilot proportional reducing valve by adjusting the duty ratio of the PWM signal value. erefore, the input pressure of LUDV multiway valve can be controlled to control the flow of the hydraulic cylinder, and then the movement of the working device can be driven by the piston rod of the hydraulic cylinder.
In order to study the nonlinear dynamic characteristics of the electrohydraulic servo system of the robotic excavator, take the boom control system for example, and the combination of low amplitude sinusoidal signal is used as excitation signal (as shown in Figure 7(a)). e stable lowspeed output signal of the boom is shown in Figure 7(b).
It can be seen from Figure 7(b) that when the excitation signal is near 0, the output remains almost constant due to dead zone characteristics; when the input is sinusoidal symmetric signal, the output shows a downward trend, which is that the response characteristic of the positive and negative direction is asymmetrical.
In order to verify the effectiveness of the adaptive Hammerstein nonlinear identification algorithm, the system model of the electrohydraulic servo system of the robotic excavator is identified by the classical Hammerstein nonlinear identification algorithm (CHNLIA) and adaptive Hammerstein nonlinear identification algorithm (AHN-LIA), respectively, which is combined with the control algorithm to form the controller, and the control performance of two algorithms is compared by experiments.
Based on the hydraulic control theory, the linear module orders of the Hammerstein model structure are set to n a � 3, n b � 1, and n k � 1; the polynomial order of the nonlinear module is set to r � 3. e basic universe of error e is [−4, 4], and the basic universe of error variation Δe is [−2, 2]. Let f 1 � 1.0, based on the adaptive Hammerstein nonlinear identification algorithm, and the model parameters of the boom control system of the robotic excavator are obtained, as shown in Table 2.
For the same test data, the model parameters of the boom control system of the robotic excavator are identified by the classical Hammerstein nonlinear identification algorithm again. e identified system model is compared with the actual system, and servo tracking experiment is carried out. e experimental results of the two algorithms are shown in Figure 8.
It can be seen from Figure 8 that the adaptive Hammerstein nonlinear identification algorithm (AHNLIA) has Mathematical Problems in Engineering better approximation accuracy than the classical Hammerstein nonlinear identification algorithm (CHNLIA). e calculated standard deviation of the CHNLIA is σ 1 � 5.8921, and the error of the identified model is obvious when the hydraulic cylinder changes its motion direction. e calculated standard deviation of the AHNLIA is σ 2 � 3.5287, the      Mathematical Problems in Engineering error of the algorithm fluctuates within a small range, and the hydraulic cylinder always stays near the designated position. It can be seen that the identification error of the AHNLIA is about 40.11% lower than that of the CHNLIA, and it can describe the dynamic characteristics of the electrohydraulic servo system of the robotic excavator more accurately.

Conclusion
On the basis of establishing the nonlinear model of the electrohydraulic servo system of the robotic excavator, the complex dynamic characteristics of the system are analyzed. Due to the nonlinear and time-varying characteristics of the electrohydraulic servo system, it is difficult to establish a precise mathematical model. e Hammerstein model with dynamic linear module and static nonlinear module is applied to describe the system; a forgetting factor is introduced to adjust the model, and the fuzzy control module is designed to adaptively update the forgetting factor. e identification results provide a reference model for the adaptive control of the robotic excavator, and the model predictive control can be used to track the predetermined trajectory accurately.
Experimental results show that the tracking accuracy of the adaptive Hammerstein nonlinear identification algorithm is improved significantly compared with the classical Hammerstein identification algorithm, which can describe the dynamic characteristics of the electrohydraulic servo system of the robotic excavator more accurately. e method provides a reference for the identification and control of other nonlinear time-varying systems in the industrial field.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.