Motion control of valve-controlled hydraulic actuators with input saturation and modelling uncertainties

Abstract

This article mainly concerns the high-performance motion control of valve-controlled hydraulic actuators with input saturation and modelling uncertainties. The nonlinear mathematic model including a continuously differentiable static friction model is constructed, and then adaptive robust design framework is adopted to cope with the modelling uncertainties, which always impede the progress of high-performance motion controller. Input saturation, which frequently exists in most physical systems, has been found to be prone to performance decay. To address this specific issue, an embedded anti-windup block containing two adjusting mechanisms is properly designed to improve the motion controller to ensure the stability and performance preservation in circumstance of input saturation, which is proved via rigorous Lyapunov analysis. Typical simulation is implemented to illustrate the availability of the proposed control method.

Keywords

Hydraulics motion control dynamic modeling nonlinear system robust design

Introduction

Hydraulic actuators have been widely utilized in modern industry and defence because of the advantages of high power density, low calorific value and the ability to apply large force/torques compared to electro-motor counterparts,¹ such as aircraft actuators,² robots and manipulators,³ rolling mill,⁴ punching machine⁵ and hardware-in-the-loop equipment.⁶ However, hydraulic actuators always have a number of modelling uncertainties and inherent nonlinear behaviours, which inevitably complicate the design and development of the high-performance motion controller. Generally, the modelling uncertainties include structured and unstructured uncertainties. Examples of structured uncertainties include the change in load and the parametric variations combined with known basis functions. Other common uncertainties, such as the unpredictable external disturbances, unmodelled friction force, dead zone, hysteresis and external leakage, whose exact describing models are not known, are called unstructured uncertainties. When the control strategy is synthesized based on the nominal mathematic model, the modelling uncertainties will inevitably cause performance degradation, or even make the system unstable.

The approaches to cope with modelling uncertainties have attracted great attentions and various practical control methods are developed. Typical examples are adaptive control,⁷ sliding mode control⁸ and intelligent control.^9,10 Besides, in the last 20 years, adaptive robust control (ARC),¹¹ which is characterized with rigorous theoretical basis, has achieved great development^3,12,13 and also provides a practical design framework for derived high-performance controller.

Although all the control strategies mentioned above have been extensively utilized and realized excellent motion tracking performance, a notable fact is that input saturation, which frequently exists in most physical systems, is rarely considered. Input saturation has been found to be prone to cause performance degradation. Generally, input saturation may be inescapable when system operates under some critical conditions, such as unmatched initial condition, serious external disturbances and physical restriction. At this time, the input value calculated by the controller may overstep the certain physical limit; control effort will hard to act on the actuators and eventually performance deterioration will be unavoidable. To address this problem, many available methods are developed, such as bounded function methods,¹⁴ model predictive control¹⁵ and anti-windup compensation.^16–18 The first method can achieve asymptotic tracking performance in the input saturation circumstance. But specific assumption about unstructured uncertainties is needed, which makes the controller conservative in some cases. The second method is capable of coping with hard control constraint. But its effectiveness may face challenges in the presence of modelling uncertainties. The last method employs an anti-windup compensator, which comes into operation once saturation occurs, to adjust the performance degradation and force the input amplitude to an acceptable level. However, it does not consider the mismatched uncertainties. Although the last method has been tried on hydraulic actuators by Guo et al.,¹⁸ the inevitable existing parametric uncertainties are neglected. Other approaches such as in Zhu et al.,¹⁹ the saturation nonlinearity is regarded as disturbance, and estimation methods are adopted to compensate the disturbance, which is passive and may cause high-gain feedback risk when the saturation is severe.

Motivated by the above analysis, an interesting idea is to develop a novel controller which can simultaneously cope with input saturation and modelling uncertainties trying to achieve excellent motion control performance. This is because in many cases, the two issues may be encountered at the same time. For this consideration, the inspiration naturally rises for this article. The main contribution of this study is to propose an integrated controller based on the ARC framework. Two adjusting filter dynamics, which make up the anti-windup block, are purposefully adopted for input saturation circumstance and can be integrated with the designed adaptive robust controller, while the stability is proved via rigorous Lyapunov analysis.

The remainder of the article is organized as follows. In section ‘Nonlinear model and preliminaries’, the nonlinear model of hydraulic linear actuators and preliminaries are given. Then, the saturated adaptive robust control (SARC) with anti-windup compensator is designed in detail and the performance results are given in section ‘Controller design and main results’. The availability of the proposed control methods is demonstrated via simulation results in section ‘Simulation experiment and performance analysis’. Section ‘Conclusion’ gives some conclusions and future work plan.

Nonlinear model and preliminaries

The considered hydraulic linear actuator is the same as in Yao et al.,⁶ and the principle diagram is illustrated in Figure 1.

Figure 1.

The principle diagram of the considered hydraulic linear actuator.

The dynamics of the inertia mass can be described as follows

m \overset{\cdot\cdot}{y} = P_{1} A_{1} - P_{2} A_{2} - F_{fr} (\overset{\cdot}{y}) - f (t, y, \overset{\cdot}{y})

(1)

where m is the inertial mass; y and $\overset{\cdot}{y}$ represent the displacement and velocity of the inertial mass, respectively; P₁ and P₂ represent the inner pressures of the two cylinder chambers; A₁ and A₂ represent the working ram areas of two sides; F_fr ( $\overset{\cdot}{y}$ ) is the friction force; and f (t, y, $\overset{\cdot}{y}$ ) represents the total unstructured uncertainties, which contains hard-to-model friction, external disturbances and so on.

For the friction dynamic in model equation (1), a continuously differentiable model is adopted as²⁰

F_{fr} (\overset{\cdot}{y}) = a_{1} [\tanh (c_{1} \overset{\cdot}{y}) - \tanh (c_{2} \overset{\cdot}{y})] + a_{2} \tanh (c_{3} \overset{\cdot}{y}) + a_{3} \overset{\cdot}{y}

(2)

where a₁, a₂ and a₃ represent the different friction levels; and c₁, c₂ and c₃ denote the various shape coefficients to approximate various friction effects. The nominal values of those friction parameters and the detailed identification experiment can be found in previous works.²¹ The static mapping relation between velocity and friction force is shown in Figure 2.

Figure 2.

Static fiction force and its curve fitting.

The pressures dynamic in the actuators can be described as²²

\begin{matrix} {\overset{\cdot}{P}}_{1} = \frac{β_{e}}{V_{1}} (- A_{1} \overset{\cdot}{y} - C_{t} P_{L} + Q_{1}) \\ {\overset{\cdot}{P}}_{2} = \frac{β_{e}}{V_{2}} (A_{2} \overset{\cdot}{y} + C_{t} P_{L} - Q_{2}) \end{matrix}

(3)

where β_e denotes the effective bulk modulus; C_t denotes the internal leakage coefficient of the cylinder; V₁ and V₂ represent the total control volume of the two chambers respectively; P_L = P₁−P₂ is the load pressure of the actuator; and Q₁ and Q₂ are the supplied flow rate and the return flow rate, respectively, which can be modelled as²²

\begin{matrix} Q_{1} = k_{q} x_{v} [s (x_{v}) \sqrt{P_{s} - P_{1}} + s (- x_{v}) \sqrt{P_{1} - P_{r}}] \\ Q_{2} = k_{q} x_{v} [s (x_{v}) \sqrt{P_{2} - P_{r}} + s (- x_{v}) \sqrt{P_{s} - P_{2}}] \end{matrix}

(4)

where x_v is the spool displacement; P_s and P_r are the supply pressure and the return pressure, respectively; $k_{q} = C_{d} w \sqrt{2 / ρ}$ ; C_d and w are the discharge coefficient and the area gradient of spool valve, respectively; ρ is the density of hydraulic oil; and s(*) is a discontinuous function defined as

s (*) = {\begin{matrix} 1, if * ⩾ 0 \\ 0, if * < 0 \end{matrix}

(5)

The valve dynamics can be neglected as in Yao et al.²³ and the function relationship of the spool displacement x_v and the control input u is a known static proportional mapping,²⁴ that is, x_v = k_i u, where k_i is a positive servo valve gain. The main reason is that the dynamic frequency of servo valve is higher than the overall system frequency. Therefore, if the dynamic of servo valve is included, additional sensors will be required to measure the spool displacement, while performance improvement may be less for motion tracking

Therefore, equation (4) can be rewritten as

\begin{matrix} Q_{1} = k_{gt} u [s (u) \sqrt{P_{s} - P_{1}} + s (- u) \sqrt{P_{1} - P_{r}}] \\ Q_{2} = k_{gt} u [s (u) \sqrt{P_{2} - P_{r}} + s (- u) \sqrt{P_{s} - P_{2}}] \end{matrix}

(6)

where k_gt = k_qk_i is the total flow gain with respect to the control input.

Define the state vector x = [x₁, x₂, x₃]^T = [y, $\overset{\cdot}{y}$ , (P₁A₁−P₂A₂)/m]^T. To facilitate the controller design, the state-space form of the system is expressed as

\begin{matrix} {\overset{\cdot}{x}}_{1} = x_{2} \\ {\overset{\cdot}{x}}_{2} = x_{3} - θ_{1} S_{f} - θ_{2} P_{f} - θ_{3} x_{2} - d \\ {\overset{\cdot}{x}}_{3} = θ_{4} gu - θ_{5} f_{1} - θ_{6} f_{2} \end{matrix}

(7)

where θ₁ = a₁/m, θ₂ = a₂/m, θ₃ = a₃/m, θ₄ = (k_gtβ_e)/m, θ₅ = β_e/m, θ₆ = (C_tβ_e)/m, S_f = tanh( $c_{1} \overset{\cdot}{y}$ ) − tanh( $c_{2} \overset{\cdot}{y}$ ), P_f = tanh( $c_{3} \overset{\cdot}{y}$ ), d = f(t, y, $\overset{\cdot}{y}$ )/m, and

\begin{matrix} R_{1} = s (u) \sqrt{P_{s} - P_{1}} + s (- u) \sqrt{P_{1} - P_{r}} > 0 \\ R_{2} = s (u) \sqrt{P_{2} - P_{r}} + s (- u) \sqrt{P_{s} - P_{2}} > 0 \\ g = \frac{A_{1}}{V_{1}} R_{1} + \frac{A_{2}}{V_{2}} R_{2} \\ f_{1} = (\frac{A_{1}^{2}}{V_{1}} + \frac{A_{2}^{2}}{V_{2}}) x_{2}, f_{2} = (\frac{A_{1}}{V_{1}} + \frac{A_{2}}{V_{2}}) P_{L} \end{matrix}

Assumption 1

The unknown system parameters and the unstructured term d meet the following features

$\hat{θ} \in Ω_{θ} at {θ : {\hat{θ}}_{min} ⩽ \hat{θ} ⩽ {\hat{θ}}_{max}}, d ⩽ δ$
(8)

where θ _min = [θ_1min, θ_2min, θ_3min, θ_4min, θ_5min, θ_6min]^T and θ _max = [θ_1max, θ_2max, θ_3max, θ_4max, θ_5max, θ_6max]^T are known constant bound vectors of θ , and δ is a known upper bound constant for d.

Assumption 2

The states x₁, x₂, P₁ and P₂ are available for measurement.

In this article, the hydraulic actuator considered is likely to be subjected to input saturation and the saturation amplitude bound is u_max, that is, $| u | ⩽ u_{max}$ . Two saturation functions used for this study are defined as

$\begin{matrix} sa t_{u_{max}} (u) = sign (u) {u_{max}, | u |} \\ sa t_{a, b} (s) = {\begin{matrix} a if s < a \\ s if a ⩽ s ⩽ b \\ b if s > b \end{matrix} \end{matrix}$
(9)

Given a desired position trajectory x_1d, the design objective is to synthesize an effective control input u such that the system output y = x₁ tracks x_1d as accurately as possible while input saturation and modelling uncertainties simultaneously exist.

Controller design and main results

Controller design

Define z₁ = x₁−x_1d as the output position tracking error, and the intermediate error signals as follows¹⁷

$\begin{matrix} z_{2} = x_{2} - α_{1} + ξ_{1} \\ z_{3} = x_{3} - α_{2} + ξ_{2} \\ {\overset{\cdot}{ξ}}_{1} = - k_{ξ 1} ξ_{1} + k_{Δ 1} Δ u \\ {\overset{\cdot}{ξ}}_{2} = - k_{ξ 2} ξ_{2} + k_{Δ 2} Δ u \end{matrix}$
(10)

where Δu = u−sat_u_max(u); ξ₁ and ξ₂ are designed to adjust the error when saturation occurs and their values are undated by filter dynamic equation (10); and k_ξ₁, k_ξ₂, k_Δ1 and k_Δ2 are positive parameters to be designed.

For the positive semi-definite function $V_{1} = z_{1}^{2} / 2$ , we can obtain

${\overset{\cdot}{V}}_{1} = z_{1} {\overset{\cdot}{z}}_{1} = z_{1} (z_{2} + α_{1} - ξ_{1} - {\overset{\cdot}{x}}_{1 d})$
(11)

Let α₁ = ${\overset{\cdot}{x}}_{1 d}$ + ξ₁–k₁z₁, then we have

${\overset{\cdot}{V}}_{1} = - k_{1} z_{1}^{2} + z_{1} z_{2}$
(12)

For the positive semi-definite function $V_{2} = V_{1} + (z_{2}^{2} / 2) + (ξ_{1}^{2} / 2)$ , we can obtain

$\begin{matrix} {\overset{\cdot}{V}}_{2} = {\overset{\cdot}{V}}_{1} + z_{2} (z_{3} + α_{2} - θ_{1} S_{f} - θ_{2} \\ P_{f} - θ_{3} x_{2} - ξ_{2} - d - {\overset{\cdot}{α}}_{1} - k_{ξ 1} ξ_{1} + k_{Δ 1} Δ u) + ξ_{1} {\overset{\cdot}{ξ}}_{1} \end{matrix}$
(13)

From equation (13), we can see that the ARC design techniques in Yao and Tomizuka¹¹ can be used to construct the visual control input α₂ as follows

$\begin{matrix} α_{2} = α_{2 a} + α_{2 s 1} + α_{2 s 2} \\ α_{2 a} = {\hat{θ}}_{1} S_{f} + {\hat{θ}}_{2} P_{f} + {\hat{θ}}_{3} x_{2} + {\overset{\cdot}{α}}_{1} + ξ_{2} - z_{1} + k_{ξ 1} ξ_{1} \\ α_{2 s 1} = - k_{s 2} z_{2} \end{matrix}$
(14)

where k_s₂ is the positive design parameter; and ${\hat{θ}}_{1}, {\hat{θ}}_{2} and {\hat{θ}}_{3}$ denote the adaptive estimations of uncertain parameters θ₁, θ₂ and θ₃.

Substituting equation (14) into equation (13), we have

$\begin{matrix} {\overset{\cdot}{V}}_{2} = - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} + z_{2} z_{3} \\ + (α_{2 s 2} + {\tilde{θ}}^{T} φ_{1} - d + k_{Δ 1} Δ u) z_{2} + ξ_{1} {\overset{\cdot}{ξ}}_{1} \end{matrix}$
(15)

where θ ^T = [θ₁, θ₂, θ₃, θ₄, θ₅, θ₆]; $\tilde{θ} = \hat{θ} - θ$ is the parametric estimate error vector; and $φ_{1}^{T}$ = [S_f, P_f, x₂, 0, 0, 0]. The nonlinear robust term α_2s2 should guarantee the following properties^6,11

$\begin{matrix} (α_{2 s 2} + {\tilde{θ}}^{T} φ_{1} - d + k_{Δ 1} Δ u) z_{2} ⩽ ε_{1} \\ α_{2 s 2} z_{2} ⩽ 0 \end{matrix}$
(16)

where ε₁ is the positive design parameter and α_2s2 is given as follows^6,11

$α_{2 s 2} = - \frac{h_{1}^{2} z_{2}}{4 ε_{1}}, h_{1} ⩾ ‖ θ_{M} ‖ ‖ φ_{1} ‖ + δ + k_{Δ 1} ‖ Δ u ‖$
(17)

where θ _M = θ _max− θ _min denotes the maximum perturbation of the uncertain parameters.

From equations (15) and (16), we have

${\overset{\cdot}{V}}_{2} = - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} - k_{ξ 1} ξ_{1}^{2} + z_{2} z_{3} + ε_{1} + k_{Δ 1} ξ_{1} Δ u$
(18)

Then, given the final positive semi-definite function $V_{3} = V_{2} + (z_{3}^{2} / 2) + (ξ_{2}^{2} / 2)$ , its derivative is obtained as

$\begin{matrix} {\overset{\cdot}{V}}_{3} = - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} - k_{ξ 1} ξ_{1}^{2} - k_{ξ 1} ξ_{2}^{2} + z_{2} z_{3} \\ + ε_{1} + z_{3} (θ_{4} gu - θ_{5} f_{1} - θ_{6} f_{2} \\ - {\overset{\cdot}{α}}_{2 c} - {\overset{\cdot}{α}}_{2 u} - k_{ξ 2} ξ_{2} + k_{Δ 2} Δ u) \\ + k_{Δ 1} ξ_{1} Δ u + k_{Δ 2} ξ_{2} Δ u \end{matrix}$
(19)

where ${\overset{\cdot}{α}}_{2 c}$ and ${\overset{\cdot}{α}}_{2 u}$ can be calculated as follows

$\begin{matrix} {\overset{\cdot}{α}}_{2} = {\overset{\cdot}{α}}_{2 c} + {\overset{\cdot}{α}}_{2 u} \\ {\overset{\cdot}{α}}_{2 c} = \frac{\partial α_{2}}{\partial t} + \frac{\partial α_{2}}{\partial x_{1}} x_{2} + \frac{\partial α_{2}}{\partial x_{2}} \hat{\overset{\cdot}{x_{2}}} + \frac{\partial α_{2}}{\partial \hat{θ}} \overset{\cdot}{\hat{θ}} + \frac{\partial α_{2}}{\partial ξ_{1}} {\overset{\cdot}{ξ}}_{1} + \frac{\partial α_{2}}{\partial ξ_{2}} {\overset{\cdot}{ξ}}_{2}, \\ {\overset{\cdot}{α}}_{2 u} = \frac{\partial α_{2}}{\partial x_{2}} \tilde{\overset{\cdot}{x_{2}}} \\ \hat{\overset{\cdot}{x_{2}}} = x_{3} + \overset{\cdot}{ξ} - {\hat{θ}}_{1} S_{f} - {\hat{θ}}_{2} P_{f} - {\hat{θ}}_{3} x_{2}, \\ \tilde{\overset{\cdot}{x_{2}}} = {\tilde{θ}}_{1} S_{f} + {\tilde{θ}}_{2} P_{f} + {\tilde{θ}}_{3} x_{2} - d \end{matrix}$
(20)

Based on the ARC design techniques¹¹ and equation (19), we can get the final controller as follows

$\begin{matrix} u = u_{a} + u_{s 1} + u_{s 2} \\ u_{a} = \frac{1}{{\hat{θ}}_{4} g} ({\hat{θ}}_{5} f_{1} + {\hat{θ}}_{6} f_{2} + {\overset{\cdot}{α}}_{2 c} + k_{ξ 2} ξ_{2} - z_{2}), \\ u_{s 1} = - (k_{s 3} - μ k_{μ}) z_{3} / g \\ \overset{\cdot}{\hat{θ}} = Proj (Γ τ) = {\begin{matrix} 0, if {\hat{θ}}_{i} > θ_{i max} and Γ τ > 0 \\ 0, if {\hat{θ}}_{i} < θ_{i min} and Γ τ < 0 \\ Γ τ \end{matrix} \\ τ = - φ_{1} z_{2} - φ_{2} z_{3} \end{matrix}$
(21)

where Proj() denotes the discontinuous mapping which can ensure that the parameter estimation process is controllable;¹¹ $\overset{\cdot}{\hat{θ}} θ$ denotes the parameter estimation vector; Г is the diagonal adaptation rate matrix; $φ_{2}^{T}$ = [S_f∂α₂/∂x₂, P_f∂α₂/∂x₂, x₂∂α₂/∂x₂, −gu_a, f₁, f₂] represents the parameter regress vector; k_s₃ and k_μ* are the positive design parameters; and μ is the adjusting parameter which can be updated from the following dynamics¹⁶

$\begin{matrix} \overset{\cdot}{λ} = sa t_{- ρ^{-}, ρ^{+}} (k_{λ} (sa {t_{0,}}_{ω} (k_{λ} | Δ u |) - λ)) \\ μ = sa t_{0, 1} (λ) \end{matrix}$
(22)

where ρ⁻, ρ⁺, ω and k_λ are positive chosen parameters.

Substituting equation (21) into equation (19), we have

$\begin{matrix} {\overset{\cdot}{V}}_{3} = - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} - k_{ξ 1} ξ_{1}^{2} - k_{ξ 1} ξ_{2}^{2} - θ_{4} (k_{s 3} - μ k_{μ}) z_{3}^{2} + ε_{1} \\ + (θ_{4} g u_{s 2} + {\tilde{θ}}^{T} φ_{2} - \frac{\partial α_{2}}{\partial x_{2}} d + k_{Δ 2} Δ u) z_{3} \\ + k_{Δ 1} ξ_{1} Δ u + k_{Δ 2} ξ_{2} Δ u \end{matrix}$
(23)

The nonlinear robust term u_s₂ should guarantee the following properties^6,11

$\begin{matrix} (θ_{4} g u_{s 2} + {\tilde{θ}}^{T} φ_{2} - \frac{\partial α_{2}}{\partial x_{2}} d + k_{Δ 2} Δ u) z_{3} ⩽ ε_{2} \\ θ_{4} g u_{s 2} z_{3} ⩽ 0 \end{matrix}$
(24)

where ε₂ is the positive design parameter and u_s₂ is chosen as^6,11

$\begin{matrix} u_{s 2} = - \frac{h_{2}^{2} z_{3}}{4 ε_{2}}, \\ h_{2} ⩾ \frac{1}{‖ θ_{4} ‖ ‖ g ‖} (‖ θ_{M} ‖ ‖ φ_{2} ‖ + σ ‖ \frac{\partial α_{2}}{\partial x_{2}} ‖ + k_{Δ 2} ‖ Δ u ‖) \end{matrix}$
(25)

Through the above design, we can further obtain

$\begin{matrix} {\overset{\cdot}{V}}_{3} ⩽ - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} - k_{ξ 1} ξ_{1}^{2} - k_{ξ 1} ξ_{2}^{2} - θ_{4} (k_{s 3} - μ k_{μ}) z_{3}^{2} \\ + ε_{1} + ε_{2} + k_{Δ 1} ξ_{1} Δ u + k_{Δ 2} ξ_{2} Δ u \end{matrix}$
(26)

It is obvious that

$\begin{matrix} k_{Δ 1} ξ_{1} Δ u ⩽ \frac{1}{4} k_{Δ 1} ξ_{1}^{2} + k_{Δ 1} ‖ Δ u ‖^{2}, \\ k_{Δ 2} ξ_{2} Δ u ⩽ \frac{1}{4} k_{Δ 2} ξ_{2}^{2} + k_{Δ 2} ‖ Δ u ‖^{2} \end{matrix}$
(27)

Furthermore, from equations (26) and (27), we have

$\begin{matrix} {\overset{\cdot}{V}}_{3} ⩽ - k_{1} z_{1}^{2} - k_{s 2} z_{2}^{2} - θ_{4} (k_{s 3} - μ k_{μ}) z_{3}^{2} \\ - (k_{ξ 1} - \frac{1}{4} k_{Δ 1}) ξ_{1}^{2} - (k_{ξ 2} - \frac{1}{4} k_{Δ 2}) ξ_{2}^{2} + ψ \end{matrix}$
(28)

where ψ = ε₁ + ε₂ + k_Δ1||Δu||² + k_Δ2||Δu||². By choosing λ_m = min{k₁, k_s₂, θ₄(k_s₃–k_μ), (k_ξ₁–k_Δ1/4), (k_ξ₂–k_Δ2 /4)}, we can infer that

${\overset{\cdot}{V}}_{3} ⩽ - 2 λ_{m} V_{3} + ψ$
(29)

Main results

If we define the compact set $Ω = {(z_{1}, z_{2}, z_{3}, ξ_{1}, ξ_{2}) : V_{3} ⩽ q}$ , it is notable that on the compact set Ω, $‖ Δ u ‖$ and ψ have the maximum values. Then by integrating both sides of equation (29), we can obtain that V₃ is upper bounded as

$V_{3} (t) ⩽ V_{3} (0) \exp (- 2 λ_{m} t) + \frac{ψ}{2 λ_{m}} [1 - \exp (- 2 λ_{m} t)]$
(30)

Thus, we know that V₃(t) is bounded and its upper bound depends on its initial value and the control gains. Furthermore, from the definition of V₃(t), we know that the tracking error z₁ and auxiliary signal z₂, ξ₁ and ξ₂ are all bounded. Because the position trajectory x_1d and its derivatives are bounded, we can infer that the control input defined in equation (21) is bounded. Finally, the proposed controller can guarantee that all the signals are bounded in the presence of input saturation and modelling uncertainties, that is, the closed-loop stability is achieved.

The selected dynamics in equation (22) can achieve the designated properties:¹⁶ when saturation occurs, the adjustable parameter μ rises to 1 rapidly to propel the amplitude of u into acceptable region while preserving the position tracking performance; when the control input is small enough, μ goes back to zero and the normal closed-loop tracking performance recovers.

Remark 1

The key point of the proposed control method is the two adjusting mechanisms designed in equations (10) and (22). They are all triggered by Δu and have different functions while cooperating with each other. The first can urge the adjusting parameters ξ₁ and ξ₂ to change towards the opposite direction of the errors and thus weaken the error fluctuation caused by input saturation. The second in equation (22) can generate the adjusting parameter μ needed by anti-windup design in the final controller.

Simulation experiment and performance analysis

The availability of the proposed high-performance motion controller for hydraulic actuator was verified via a simulation model shown in Figure 3. Specifications of hydraulic system components are listed in Table 1.

Figure 3.
The simulation model utilized for performance verification.

Table 1.
Specifications of the hydraulic actuators.

m (kg) 30 A₁ = A₂ (m²) 9 × 10⁻⁴

β_e (MPa) 700 k_gt (m⁴/(s·V·N^−1/2)) 1.18 × 10⁻⁸

V₀₁ = V₀₂ (m³) 4 × 10⁻⁵ C_t (m⁵/(N·s)) 9 × 10⁻¹²

p_s (MPa) 20 p_r (MPa) 0

a₁ (N·s/m) 177 a₂ (N·s/m) 63

a₃ (N·s/m) 4000

The following typical controllers are used for comparison:

SARC: This is the proposed saturated adaptive robust control with anti-windup compensation and the detailed form is given in section ‘Nonlinear model and preliminaries’. Parameters of the controller are chosen as in Table 2.

ARC: This is the known adaptive robust control proposed in Yao and Tomizuka,¹¹ which is commonly used in industries. Based the constructed system model in equation (7), the controller has the following form

$\begin{matrix} z_{2} = x_{2} - υ_{1}, z_{3} = x_{3} - υ_{2}, υ_{1} = {\overset{\cdot}{x}}_{1 d} - k_{1} z_{1}, υ_{2} = {\hat{θ}}_{1} S_{f} + {\hat{θ}}_{2} P_{f} + {\hat{θ}}_{3} x_{2} + {\overset{\cdot}{υ}}_{1} - z_{1} - k_{s 2} z_{2} \\ u = \frac{1}{{\hat{θ}}_{4} g} ({\hat{θ}}_{5} f_{1} + {\hat{θ}}_{6} f_{2} + {\overset{\cdot}{υ}}_{2}) - k_{s 3} z_{3} / g, \overset{\cdot}{\hat{θ}} = Pro j_{\hat{θ}} (Γ τ), τ = - ψ_{2} z_{2} - ψ_{3} z_{3} \\ ψ_{2} = {[S_{f}, P_{f}, x_{2}, 0, 0, 0]}^{T}, ψ_{3} = {[S_{f} \frac{\partial υ_{2}}{\partial x_{2}}, P_{f} \frac{\partial υ_{2}}{\partial x_{2}}, x_{2} \frac{\partial υ_{2}}{\partial x_{2}}, - gu, f_{1}, f_{2}]}^{T} \end{matrix}$
(31)

Table 2.
Parameters of the proposed controller.

ρ ⁻ 0.2 ρ ⁺ 30

ω 1 k_λ 10

k_ξ ₁ 3 k_ξ ₂ 12

k _Δ1 1500 k _Δ2 1000

k ₁ 50 k_s ₂ 100

k_s ₃ 120 k_μ 119

Meaning of the ARC parameters is the same as the corresponding parameters in SARC. The desired tracking trajectory is chosen as x_1d = 100 sin(πt). The saturation amplitude u_max is set as 2 V. The simulation results are shown in Figures 4 –6.

Figure 4.
Tracking errors.

Figure 5.
Control inputs.

Figure 6.
Curve of Δu and parameter μ.

The position tracking errors of the two compared controllers are exhibited in Figure 4. Obviously, for both the controllers, the tracking error shows large fluctuation once the control input touches the saturation bound. What is more notable is that the proposed SARC controller, especially the filter designs in equation (10), can adjust the tracking error timely and effectively, which can be further revealed by comparing the maximum tracking error calculated via M_e = max{|z₁|}. The maximum tracking errors of SARC and ARC are 1.427 and 1.172 mm, respectively. This means a 17.8% reduction, which verified the desirability and positive significance of the proposed controller for which the maximum tracking error is a key performance indicator in many hydraulic applications. It is easy to see that all of them gradually return to normal error level after saturation disappears.

The control inputs are exhibited in Figure 5. Obviously, the saturation occurs once the amplitude exceeds 2 V, which corresponds to the error fluctuation in Figure 4. The difference is that for SARC, the control input will recover to acceptable region faster than ARC, which benefits from the ingenious design of u_s₁. The change curve of Δu and the updated dynamics given in equation (22) are shown in Figure 6. Evidently, the parameter μ rises to 1 rapidly when saturation occurs and goes back to 0 after saturation disappears, thus ensuring a timely adjusting manner for control input.

What needs illustration is that the designed adjusting anti-windup compensator will inevitably influence the tracking error dynamics, especially slow down the rate of convergence after saturation disappears (e.g. 1–2 and 4.5–5.5 time intervals) and cause the unexpected peaks (e.g. 3.5–4 time interval). How to improve them is certainly worth further concerning and automatically becomes the research direction for our future efforts.

Conclusion

In this article, a novel saturated motion tracking controller for hydraulic actuators is proposed to address the issues of input saturation and modelling uncertainties. The modelling uncertainties are handled well via ARC based on the nonlinear mathematic model including a continuously differentiable static friction model. The integrated anti-windup compensator containing two adjusting filter dynamics, which respectively adjust the control input and tracking error, is utilized to cope with the saturation circumstance. The stability is proved via Lyapunov analysis. In comparison with traditional ARC, the maximum tracking error value decreases significantly when saturation occurs, which verified the desirability and positive significance of the proposed controller. The future works may focus on how to improve the rate of convergence after saturation disappears and the unexpected peaks, as well as applying the proposed control strategy to other practical systems, such as vehicles, robot manipulators and hardware-in-the-loop systems.

Footnotes

Handling Editor: Zengtao Chen

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: This work was supported by the National Natural Science Foundation of China (grant no. 51605145); the Henan Provincial Key Scientific Research Project of Colleges and Universities of China (grant no. 18A460020); and the Henan Provincial Science and Technology Project of China (grant no. 172102410071).

ORCID iD

Zhenle Dong

References

Tao

Feng

Research on dynamics and stability in the stairs-climbing of a tracked mobile robot. Int J Adv Robot Syst 2012; 9: 52850.

Ouyang

Zhao

et al . Integral variable PI control on flow load simulator of aircraft hydraulic system. J Zhejiang Univ 2017; 51: 1111–1118.

Mohanty

Yao

Indirect adaptive robust control of hydraulic manipulators with accurate parameter estimates. IEEE T Contr Syst T 2011; 19: 567–575.

Zhao

Fang

Zhong

Design of adaptive robust output feedback dynamic controller for thickness control in a cold strip rolling mill. Control Theory Appl 2008; 25: 787–790.

Song

Liu

Pan

et al . Investigation on the design of hydraulic high-speed punching press with a large press force and a long press stroke. Adv Mech Eng 2016; 8: 83601.

Yao

Jiao

Yao

et al . Nonlinear adaptive robust force control of hydraulic load simulator. Chin J Aeronaut 2012; 25: 766–775.

Garagic

Srinivasan

Application of nonlinear adaptive control techniques to an electrohydraulic velocity servomechanism. IEEE T Contr Syst T 2002; 12: 303–314.

Guan

Pan

Adaptive sliding mode control of electro-hydraulic system with nonlinear unknown parameters. Control Eng Pract 2008; 16: 1275–1284.

Truong

Ahn

KK.

Force control for hydraulic load simulator using self-tuning grey predictor–fuzzy PID. Mechatronics 2009; 19: 233–246.

10.

Hanh

Ahn

Kha

et al . Trajectory control of electro-hydraulic excavator using fuzzy self tuning algorithm with neural network. J Mech Sci Technol 2009; 23: 149–160.

11.

Yao

Tomizuka

Adaptive robust control of SISO nonlinear systems in a semi-strict feedback form. Automatica 1997; 33: 893–900.

12.

Yao

Reedy

et al . Adaptive robust motion control of single-rod hydraulic actuators: theory and experiments. IEEE/ASME T Mech 2000; 5: 79–91.

13.

Qiu

Performance-oriented asymptotic tracking control of hydraulic systems with radial basis function network disturbance observer. Adv Mech Eng 2016; 8: 2016.

14.

Dong

Yao

Robust asymptotic tracking control of a chain of integrator nonlinear systems with input constraint and hysteresis nonlinearity. T I Meas Control 2016; 38: 1500–1508.

15.

Cuzzola

Geromel

Morari

An improved approach for constrained robust model predictive control. Automatica 2002; 38: 1183–1189.

16.

Teel

Zaccarian

Marcinkowski

JJ.

An anti-windup strategy for active vibration isolation systems. Automatica 1994; 30: 1869–1883.

17.

Sun

Zhao

Gao

Saturated adaptive robust control for active suspension systems. IEEE T Ind Electron 2013; 60: 3889–3896.

18.

Guo

Wang

Liu

Anti-windup control of electro-hydraulic system with load disturbance and modelling uncertainty. IEEE T Ind Inform 2018; 14: 3097–3108.

19.

Zhu

Wang

Zhou

et al . Adaptive neural network saturated control for MDF continuous hot pressing hydraulic system with uncertainties. IEEE Access 2018; 6: 2266–2273.

20.

Makkar

Dixon

Sawyer

et al . A new continuously differentiable friction model for control systems design. In: Proceedings of IEEE/ASME international conference on advanced intelligent mechatronics, Monterey, CA, 24–28 July 2005, pp.600–605. New York: IEEE.

21.

Yue

Yao

Integral robust based asymptotic tracking control of electro-hydraulic load simulator. Acta Aeronaut Astronaut Sin 2017; 38: 289–298.

22.

Merritt

HE.

Hydraulic control system. New York: John Wiley & Sons, 1968.

23.

Yao

Jiao

et al . High-accuracy tracking control of hydraulic rotary actuators with modeling uncertainties. IEEE/ASME T Mech 2014; 19: 633–641.

24.

Chen

Meguid

SA.

The transient response of a piezoelectric strip with a vertical crack under electromechanical impact load. Int J Solid Struct 2000; 37: 6051–6062.

m (kg)	30	A₁ = A₂ (m²)	9 × 10⁻⁴
β_e (MPa)	700	k_gt (m⁴/(s·V·N^−1/2))	1.18 × 10⁻⁸
V₀₁ = V₀₂ (m³)	4 × 10⁻⁵	C_t (m⁵/(N·s))	9 × 10⁻¹²
p_s (MPa)	20	p_r (MPa)	0
a₁ (N·s/m)	177	a₂ (N·s/m)	63
a₃ (N·s/m)	4000

ρ ⁻	0.2	ρ ⁺	30
ω	1	k_λ	10
k_ξ ₁	3	k_ξ ₂	12
k _Δ1	1500	k _Δ2	1000
k ₁	50	k_s ₂	100
k_s ₃	120	k_μ	119