Open AccessArticle

A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators

Zengcheng Zhou

Guoyuan Tang

^1,*,

Ruikun Xu

¹,

Lijun Han

¹ and

Maolin Cheng

^2,3,4

School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

CCCC Second Harbor Engineering Co., Ltd., Wuhan 430074, China

Key Laboratory of Large–Span Bridge Construction Technology, Wuhan 430074, China

⁴

Research and Development Center of Transport Industry of Intelligent Manufacturing Technologies of Transport Infrastructure, Wuhan 430040, China

Author to whom correspondence should be addressed.

J. Mar. Sci. Eng. 2021, 9(3), 269; https://doi.org/10.3390/jmse9030269

Submission received: 8 February 2021 / Revised: 22 February 2021 / Accepted: 23 February 2021 / Published: 3 March 2021

(This article belongs to the Section Ocean Engineering)

Download

Browse Figures

Versions Notes

Abstract

In this paper, the tracking control problem of underwater robot manipulators is investigated under the influence of the lumped disturbances, including unknown ocean current disturbances and parameter uncertainties. The proposed novel continuous nonsingular finite–time (CNFT) control method is twofold. Firstly, the modified adaptive super–twisting algorithm (ASTA) is proposed with a nonsingular fast terminal sliding mode (NFTSM) manifold to guarantee the finite–time convergence both in the sliding mode phase and the reaching phase. Secondly, a higher–order super–twisting disturbance observer (HOSTDO) is exploited to attenuate the effects of the lumped disturbances. Considering the time–varying gain matrix of the closed–loop control system, the bounded stability is strictly proved via the Lyapunov theory. Hence, the superiority of the proposed controller is singularity–free, fast convergence, chattering–free, high steady–state tracking performance, and good robustness by resorting to the methods of CNFT control and ASTA in combination with a disturbance observer. Finally, numerical simulations are conducted on a two degree–of–freedom (DOF) underwater robot manipulator to demonstrate the effectiveness and high tracking performance of the designed controller.

Keywords:

underwater robot manipulators; continuous finite–time control; adaptive super–twisting algorithm; disturbance observer; time–varying gain matrix

1. Introduction

The last few decades have seen giant undersea technological improvements on underwater robot manipulators in many scientific fields, such as oceanographic observation, sub–sea detection, and some military applications [1]. Considering the complex underwater circumstances, it is quite difficult to complete tasks with high control precision and a fast dynamic response [2]. Hence, the impact of the hydrodynamic forces, such as drag torques and buoyancy torques, should be considered [3]. Minglu Zhang et al. performed the modeling analysis and the simulation of a viscous hydrodynamic model of single–DOF underwater manipulators [4]. In addition, the hydrodynamic regulation was analyzed for underwater manipulators [5]. However the tracking problem was not included. Considering the roughly modeled hydrodynamics, some robust control strategies have been proposed for underwater manipulators in the last few decades to achieve satisfactory tracking performance. A proportional–integral–derivative (PID) control scheme was designed for the underwater cable–driven hyper–redundant manipulator, which experimentally achieved underwater grasping and load carrying [6]. A model–reference adaptive control was proposed for a hydrodynamic underwater manipulator to deal with the nonlinear dynamics and hydrodynamics [7]. Besides, sliding mode control was designed for the underwater manipulator but the hydrodynamic effect was only regarded as the external disturbances [8]. However, the robustness and the precision can not be ensured at the same time. In 2020, an adaptive sliding mode control was proposed for the underwater manipulator to achieve finite time tracking with the asymmetric saturation [9]. Chao Yang et al. combined the adaptive control, sliding mode and PID control to address the joint tracking control problem [10]. However, the proposed sliding mode controllers still suffer from the chattering problem, which decreases tracking performance. Hence, it is still a great challenge to obtain high precision and strong robustness tracking results.

Sliding mode (SM) control has been wildly considered in many nonlinear systems with significant uncertainties [11]. In recent years, different kinds of sliding mode control, including robust sliding mode control [12], multiple sliding mode control [13], and dynamic sliding mode control [14], have been widely investigated. Specially, the nonsingular fast terminal sliding mode (NFTSM) control was proposed to achieve fast convergence and singularity–free at the same time [15,16,17]. However, the discontinuous nature of the control law always leads to an undesirable chattering problem [18]. Thus the higher–order SM (HOSM) method was proposed to reduce the chattering effect [19,20]. The core idea of HOSM control is to design a continuous sliding variable which is differentiable up to the higher order, so as to ensure the continuity of control signals [21]. In the SM control community, the super–twisting algorithm (STA) is one of the most effective HOSM techniques to achieve a high control accuracy [22,23]. However, these control methods require the prior information of the lumped disturbances, which is not easily accessible for the underwater robot manipulator system in a submarine environment. In order to tackle this problem, the adaptive method for the switching function, which is known as adaptive STA (ASTA), can be used [24,25].

Although TSM control with STA can be demonstrated stable through strict Lyapunov theory [23], it is difficult to prove the stability when combining NFTSM control with STA because of the time–varying gain matrix. Some researches directly utilized STA to eliminate the chattering without strict proof [26], while others adopted an integral sliding surface to avoid this issue [27]. Therefore, it is meaningful to take the time–varying gain matrix into consideration when proving the stability of the control system through strict Lyapunov theory.

It is worth noting that the STA can also be used for system observation [28]. J. A. Moreno demonstrated the stability and robustness of a second–order super–twisting observer for a nonlinear system by using strict Lyapunov functions [29]. Besides, a super–twisting observer was proposed for a Stewart platform to solve the forward kinematics problem [30]. However, most of these super–twisting observers are used to estimate the state variables, and a minority of them are considered to approximate the lumped disturbances.

Motivated by the above mentioned studies, our target is to propose a novel CNFT control for the underwater robot manipulators. The contributions of our paper are concluded as follows: (1) A novel HOSTDO–based NFTSM control with the modified ASTA scheme is firstly proposed for the underwater robot manipulators considering the unknown lumped disturbances. It provides a new theoretical way to resolve the chattering problem, which is the main obstacle for the application of SM control. (2) The modified ASTA is formulated to reserve the merits of traditional STA while prior information of the lumped disturbances is not needed. Additionally, it has a faster convergence rate than the standard STA. (3) The disturbance attenuation is achieved through the HOSTDO, which can strongly improve robustness of the control system. (4) Considering the time–varying gain matrix due to the combination of NFTSM and modified ASTA, the corresponding closed–loop stability of the control system is strictly verified through the Lyapunov theory.

The rest of our paper is organized as follows. The dynamic model of the underwater robot manipulators and some fundamental facts are presented in Section 2. The detailed structural procedure of the control strategy is reported in Section 3. In Section 4, the proposed control is applied to a two–DOF underwater robot manipulator. The trajectory tracking performances are also compared with those of the other three controllers. Finally, concluding remarks are given in Section 5.

2. Problem Formulation

2.1. Underwater Robot Manipulator Model and Properties

The n–joint underwater robot manipulators can be described by the following Euler–Lagrange equation [31,32]:

M (q) \ddot{q} + C (q, \dot{q}) \dot{q} + G (q) + τ_{w} (q, \dot{q}) = τ_{e d i s} + τ

(1)

where

q, \dot{q}, \ddot{q} \in R^{n}

denote the position, velocity, and acceleration vectors of joints, respectively;

M (q) \in R^{n \times n}

C (q, \dot{q}) \in R^{n}

G (q) \in R^{n}

are the inertia matrix, Coriolis/centrifugal matrix, and gravitational vector, respectively;

τ_{w} (q, \dot{q}) \in R^{n}

denotes the hydrodynamic force vector, which contains drag torque

τ_{D}

and buoyancy torque

τ_{B}

and satisfies

τ_{w} = τ_{D} + τ_{B}

;

τ_{e d i s} \in R^{n}

is the unknown external disturbance vector; and

τ \in R^{n}

denotes the control input torque vector.

Assumption 1.

The model parameters can be partitioned into different parts which can be given as follows:

M (q) = M_{0} (q) + M_{Δ} (q), C (q, \dot{q}) = C_{0} (q, \dot{q}) + C_{Δ} (q, \dot{q}), G (q) = G_{0} (q) + G_{Δ} (q), τ_{w} (q, \dot{q}) = τ_{w 0} (q, \dot{q}) + τ_{w Δ} (q, \dot{q})

, where

M_{0} (q)

C_{0} (q, \dot{q})

G_{0} (q)

, and

τ_{w 0} (q, \dot{q})

are the nominal parts and

M_{Δ} (q)

C_{Δ} (q, \dot{q})

G_{Δ} (q)

and

τ_{w Δ} (q, \dot{q})

denote the uncertain bounded parts.

Assumption 2.

The target trajectory and its derivatives

q_{d}, {\dot{q}}_{d}, {\ddot{q}}_{d} \in R^{n}

are known and bounded.

To simplify the design and analysis of the proposed controller, we define

x_{1} = q

x_{2} = \dot{q}

, and considering Assumption 1, the dynamic model (1) can be rewritten as:

\begin{matrix} {\dot{x}}_{1} & = x_{2} \\ {\dot{x}}_{2} & = {M_{0}}^{- 1} (x_{1}) u (t) + f (x_{1}, x_{2}) + Δ (x_{1}, x_{2}, t) \\ y & = x_{1} \end{matrix}

(2)

where

f (x_{1}, x_{2}) = {M_{0}}^{- 1} [- C_{0} x_{2} - G_{0} - τ_{w}]

u (t) = τ

is the time–varying control input, and

Δ \in R^{n}

represents the lumped disturbance vector which is defined as

Δ (x_{1}, x_{2}, t) = M_{0}^{- 1} [- M_{Δ} x_{2} - C_{Δ} x_{2} - G_{Δ} - τ_{w Δ} + τ_{e d i s}]

The control objective is to design a proper controller so that high tracking precision and fast convergence can be obtained without chattering and singularity problems.

Remark 1.

The fluid in the underwater circumstance should be irrational and unbounded, thus the rotation of the fluid could be very small compared with that of the rigid body.

Remark 2.

The added mass torque and the vortex influence are not considered in this paper. Considering the fact that the acceleration of the underwater manipulator body is always slow, the added mass torque and the vortex influence are small enough to be neglected in practical underwater robotic applications.

2.2. Fundamental Facts

Our control method will be accomplished based on some fundamental facts to achieve the finite–time convergence and system stability.

Definition 1.

For the sake of brevity and clarity, we define that

d i a g \{\cdot\} \in R^{n \times n}

stands for the positive definite diagonal matrix, and the vectors

s g n (ς) \in R^{n}

and

ς^{[r]} \in R^{n}

are as follows:

s g n (ς) = {[s i g n (ς_{1}), \dots, s i g n (ς_{n})]}^{T}

ς^{[r]} = [{|ς_{1}|}^{r} s i g n (ς_{1}), \dots,

{{|ς_{n}|}^{r} s i g n (ς_{n})]}^{T}

, where

ς = {[ς_{1}, \dots, ς_{n}]}^{T} \in R^{n}

is the vector variable,

r > 0

stands for a known constant, and when

r \geq 1

, it can be verified that

\frac{d}{d t} ς^{[r]} = r {|ς|}^{r - 1} \dot{ς}

. Besides,

s i g n (\cdot)

is a standard symbolic function.

Definition 2.

The powers of error vectors are defined as:

\begin{matrix} e^{[φ]} & = {[e_{1}^{[φ]}, \dots, e_{n}^{[φ]}]}^{T} \in R^{n} \\ {\dot{e}}^{[l / p]} & = {[{\dot{e}}_{1}^{[l / p]}, \dots, {\dot{e}}_{n}^{[l / p]}]}^{T} \in R^{n} \\ {\dot{e}}^{[2 - l / p]} & = {[{\dot{e}}_{1}^{[2 - l / p]}, \dots, {\dot{e}}_{n}^{[2 - l / p]}]}^{T} \in R^{n} \\ {|e|}^{φ - 1} & = d i a g \{{|e_{1}|}^{φ - 1}, \dots, {|e_{n}|}^{φ - 1}\} \in R^{n \times n} \\ {|\dot{e}|}^{(l / p) - 1} & = d i a g \{{|{\dot{e}}_{1}|}^{(l / p) - 1}, \dots, {|{\dot{e}}_{n}|}^{(l / p) - 1}\} \in R^{n \times n} . \end{matrix}

(3)

Lemma 1

([33]). Considering a system described in (2), and any real numbers

a > 0

and

0 < b < 1

so that the Lyapunov function

V (x)

is positive definite and

\dot{V} (x) + a V^{b} (x)

is negative semi–definite, the control system will be stable and all the signals will be bounded in finite time, and the settling time can be estimated by

T_{0} \leq \frac{1}{a (1 - b)} V^{1 - b} (x_{0})

3. The Design of CNFT Control Based on HOSTDO and ASTA

In this section, a novel CNFT control method based on ASTA is proposed for underwater manipulators. First, a NFTSM manifold is designed to be free from the singularity problem. Second, the modified STA is derived from the adaptive method to alleviate the chattering phenomenon and obtain good dynamical performance. Besides, the parameter adaptation method ensures that there is no need for prior information of the upper bound of the lumped disturbances. Finally, HOSTDO is utilized to estimate and compensate the lumped disturbances.

3.1. The Design of HOSTDO

STA is an important method for the control and the observation. The HOSTDO dynamics is given in the following form:

\begin{matrix} {\dot{\hat{x}}}_{1} & = {\hat{x}}_{2} + z_{1} \\ {\dot{\hat{x}}}_{2} & = {\hat{x}}_{3} + M_{0}^{- 1} u + f + z_{2} \\ {\dot{\hat{x}}}_{3} & = {\hat{x}}_{4} + z_{3} \\ {\dot{\hat{x}}}_{4} & = z_{4} \end{matrix}

(4)

where the state vectors are given as

x_{3} = Δ (x_{1}, x_{2}, t)

x_{4} = \dot{Δ} (x_{1}, x_{2}, t)

;

{\hat{x}}_{1}

{\hat{x}}_{2}

{\hat{x}}_{3}

{\hat{x}}_{4}

are the estimations of state vectors

x_{1}

x_{2}

x_{3}

x_{4}

, respectively; and

z_{1}

z_{2}

z_{3}

, and

z_{4}

are the correction terms. And for simplicity, the initial values of

z_{1}

z_{2}

z_{3}

z_{4}

are assumed to be zero. Then we define the state estimate error variables

ε_{1} = x_{1} - {\hat{x}}_{1}

ε_{2} = x_{2} - {\hat{x}}_{2}

ε_{3} = x_{3} - {\hat{x}}_{3}

, and

ε_{4} = x_{4} - {\hat{x}}_{4}

. The correction terms are defined as

z_{1} = ξ_{1} {ε_{1}}^{[3 / 4]}

z_{2} = ξ_{2} {ε_{1}}^{[1 / 2]}

z_{3} = ξ_{3} {ε_{1}}^{[1 / 4]}

z_{4} = ξ_{4} s g n (ε_{1})

, in which

ξ_{1}

ξ_{2}

ξ_{3}

, and

ξ_{4}

are positive constants.

Using (2), we obtain the observer dynamic error model as:

\begin{matrix} {\dot{ε}}_{1} & = - ξ_{1} {ε_{1}}^{[3 / 4]} + ε_{2} \\ {\dot{ε}}_{2} & = - ξ_{2} {ε_{1}}^{[1 / 2]} + ε_{3} \\ {\dot{ε}}_{3} & = - ξ_{3} {ε_{1}}^{[1 / 4]} + ε_{4} \\ {\dot{ε}}_{4} & = - ξ_{4} s g n (ε_{1}) + \ddot{Δ} . \end{matrix}

(5)

By selecting the gains

ξ_{1}

ξ_{2}

ξ_{3}

, and

ξ_{4}

according to [34], we can obtain that the proposed HOSTDO reaches finite–time stability and

ε_{i} = 0, i = 1, 2, 3, 4

, which implies

{\dot{ε}}_{i} = 0, i = 1, 2, 3, 4

. Then the estimates of disturbances are given as follows

\hat{Δ} = - ξ_{2} {ε_{1}}^{[1 / 2]} + v_{1}, {\dot{v}}_{1} = - ξ_{3} {ε_{1}}^{[1 / 4]} + v_{2}, {\dot{v}}_{2} = - ξ_{4} s g n (ε_{1})

, where the initial value of

\hat{Δ}

v_{1}

v_{2}

are supposed to be 0 for simplicity.

Remark 3.

The proposed disturbance observer generates a more smooth output to estimate and compensate the unknown lumped disturbances which can obviously alleviate the chattering issue. Furthermore, Equation (4) can be used to approximate the system state vectors in other research [35,36]. The convergence and robustness of the HOSTDO can be verified by Lyapunov theory, and the detailed proof can be referred to [34].

3.2. HOSTDO Based CNFT Control Design with ASTA Method and Time–Varying Gain Matrix

To ensure the chattering–free ability and eliminate the need for the prior information of the upper bound, the modified STA scheme is utilized with the parameter adaptation method to obtain good comprehensive performances. A novel HOSTDO–based CNFT control scheme with the ASTA method is proposed in this section.

Define the tracking error as

e = q - q_{d}

with

q_{d}

being the desired trajectory. To ensure the finite–time convergence ability and singularity–free performance, a NFTSM manifold is designed as [26]:

s = e + σ_{1} e^{[φ]} + σ_{2} {\dot{e}}^{[l / p]}

(6)

where

s = {[s_{1}, \dots, s_{n}]}^{T} \in R^{n}

is the sliding vector;

σ_{1}

σ_{2}

are the designed positive constants; l and p are positive odd numbers satisfying

1 < l / p < 2

; and

φ

is a designed constant satisfying

φ > l / p

The derivative of (6) is:

\begin{matrix} \dot{s} & = \dot{e} + σ_{1} φ {|e|}^{φ - 1} \dot{e} + σ_{2} \frac{l}{p} {|\dot{e}|}^{l / p - 1} (\ddot{q} - {\ddot{q}}_{d}) \\ = \dot{e} + σ_{1} φ {|e|}^{φ - 1} \dot{e} + σ_{2} \frac{l}{p} {|\dot{e}|}^{l / p - 1} ({M_{0}}^{- 1} (x_{1}) u + f (x_{1}, x_{2}) + Δ (x_{1}, x_{2}, t) - {\ddot{q}}_{d}) . \end{matrix}

(7)

It can be deduced from (7) that

φ - 1 > 0

and

l / p - 1 > 0

are valid considering

1 < l / p < 2

and

φ > l / p

, so the singularity phenomenon can be avoided successfully. Meanwhile, the system’s dynamic convergence speed is enhanced by combining the merits of the NFTSM manifold.

According to the NFTSM manifold and adaptive modified STA scheme, the CNFT control law based on the HOSTDO is designed as:

u = u_{e q} + u_{r e}

(8)

where

u_{e q} = M_{0} (x_{1}) [\begin{matrix} - \frac{1}{σ_{2}} \frac{p}{l} ({\dot{e}}^{[2 - l / p]} + σ_{1} φ {|e|}^{φ - 1} \cdot {\dot{e}}^{[2 - l / p]}) - f (x_{1}, x_{2}) + {\ddot{q}}_{d} - \hat{Δ} (x_{1}, x_{2}, t) \end{matrix}]

is designed to control the nominal parts, and

u_{r e} = - M_{0} (x_{1}) [α \frac{s}{{∥s∥}^{1 / 2}} + k_{1} s + β \int_{0}^{t} \frac{s}{∥s∥} d t + k_{2} \int_{0}^{t} s d t]

is used to compensate the lumped disturbances, where

k_{1}

and

k_{2}

are the designed positive constants,

α

and

β

are the adaptation parameters, and the adaptive gains are designed as:

\begin{matrix} \dot{α} & = \{\begin{matrix} ω_{1} \sqrt{\frac{Φ_{1}}{2}} s i g n (|s_{i}| - ρ_{0}) & i f α > α_{m} \\ r_{0} & i f α \leq α_{m} \end{matrix} \\ \dot{β} & = 2 κ \dot{α} \end{matrix}

(9)

where

κ = \frac{ω_{2}}{2 ω_{1}} \sqrt{\frac{Φ_{2}}{Φ_{1}}} > 0

;

α_{m}

is a small positive constant;

ω_{1}

ω_{2}

Φ_{1}

Φ_{2}

ρ_{0}

, and

r_{0}

are arbitrary positive constants. For the sake of system stability in the control process, the gain parameters

α

β

k_{1}

, and

k_{2}

should satisfy the following inequities:

\begin{matrix} α > 0, β > 4 α^{2} \frac{λ_{M}^{4}}{λ_{m}^{3}}, k_{1} > \frac{3}{5} \frac{λ_{M}}{λ_{m}} δ, k_{2} > 16 k_{1}^{2} \frac{λ_{M}^{4}}{λ_{m}^{3}} \end{matrix} .

(10)

At the same time, they should also satisfy:

\begin{matrix} α > 0, β > 0, k_{1} > max \{2 \frac{λ_{M}}{λ_{m}} δ, δ \sqrt{\frac{9}{14} \frac{λ_{M}^{3}}{λ_{m}^{3}}}\}, \\ k_{2} > \frac{3 k_{1}^{2} λ_{M}^{4} δ + \frac{1}{2} k_{1} λ_{M}^{4} δ^{2} {+ k}_{1}^{3} λ_{M} (2 λ_{M}^{3} - λ_{m}^{3})}{λ_{m}^{2} (k_{1} λ_{m} - 2 λ_{M} δ)}, \end{matrix}

(11)

thus the robustness and stability of the control system can be ensured.

The block diagram of our newly proposed control method is given in Figure 1.

It can be observed that the proposed CNFT control scheme (8) contains no discontinuous terms compared with traditional SMC that uses the discontinuous function

s g n (s)

thanks to the utilization of ASTA. Hence, the chattering issue can be successfully avoided. However, prior information of the lumped disturbances is not always accessible in real systems. To overcome this control limitation, the adaptation method is applied to approximate the designed parameters of the STA reaching control law, so that the prior information of the upper bound is not needed. On the other hand, NFTSM surface is designed in the proposed controller, thus the high control accuracy and fast convergence speed can be achieved and the singularity problem can be eliminated at the same time. The stability of the proposed control scheme is presented in Theorem 1.

Theorem 1.

For the underwater robot manipulator system (2), if an HOSTDO (4) and (5), NFTSM manifold (6), and control input (8) are proposed, the designed parameters for the ASTA are selected to satisfy (10) and (11), then the system position signal q will reach the given desired position signal

q_{d}

in finite time, and it can be ensured that there is no occurrence of singularity and chattering problem during the whole tracking process.

Proof.

Several similar control methods have been briefly proven in some other studies [37].

Substituting the CNFT control law (8) into (7), we have:

\begin{matrix} \dot{s} = σ_{2} \frac{l}{p} {|\dot{e}|}^{l / p - 1} (- α \frac{s}{{∥s∥}^{1 / 2}} - k_{1} s - β \int_{0}^{t} \frac{s}{∥s∥} d t - k_{2} \int_{0}^{t} s d t + \tilde{Δ} (x_{1}, x_{2}, t)) \end{matrix} .

(12)

We define the time–varying gain matrix P as follows:

P = σ_{2} \frac{l}{p} {|\dot{e}|}^{l / p - 1} .

(13)

Assumption 3.

The time–varying gain matrix is symmetric and positive definite. The largest and smallest eigenvalues of the time–varying gain matrix P are assumed to be bounded by upper (

λ_{M}

) and lower (

λ_{m}

) bound values respectively.

Assumption 4.

The estimation error of the lumped disturbances also satisfies the following inequality such that

∥D∥ \leq δ ∥s∥

, where

D = \tilde{Δ} (x_{1}, x_{2}, t)

δ > 0

is a positive known constant, and

\tilde{Δ} (x_{1}, x_{2}, t) = Δ (x_{1}, x_{2}, t) - \hat{Δ} (x_{1}, x_{2}, t)

is the bounded observer estimate error vector.

Considering (13) and Assumption 4, then (12) can be rewritten as:

\begin{matrix} \dot{s} & = - α P \frac{s}{{∥s∥}^{1 / 2}} - k_{1} P s + P v + P D \\ \dot{v} & = - β \frac{s}{∥s∥} - k_{2} s \end{matrix}

(14)

where

v = - β \int_{0}^{t} \frac{s}{∥s∥} d t - k_{2} \int_{0}^{t} s d t

Remark 4.

It should be noted that in the NFTSM control method, the time–varying gain matrix P cannot be eliminated during the design process compared with the nominal STA–based sliding mode control [38], which makes the demonstration more difficult.

Then, select the Lyapunov function candidate as:

V = V_{1} (s, v) + \frac{1}{2 μ_{1}} {(α - α^{*})}^{2} + \frac{1}{2 μ_{2}} {(β - β^{*})}^{2}

(15)

where

α^{*}, β^{*}

are the upper bounds of

α, β

, respectively and

μ_{1}, μ_{2} > 0

are positive constants. The function

V_{1} (s, v)

is given as:

\begin{matrix} V_{1} (s, v) = (2 β + \frac{α^{2}}{2} λ_{M}) ∥s∥ + (k_{2} + \frac{k_{1}^{2}}{2} λ_{M}) s^{T} s \\ + v^{T} P v + α k_{1} λ_{M} \frac{s^{T} s}{{∥s∥}^{1 / 2}} - k_{1} s^{T} P v - α \frac{v^{T} P s}{{∥s∥}^{1 / 2}} . \end{matrix}

(16)

It can be demonstrated that

V_{1} (s, v)

is positive definite since:

V_{1} \geq 2 β ∥s∥ + k_{2} s^{T} s + \frac{1}{2} v^{T} P v + \frac{1}{2} χ^{T} P χ

(17)

where

χ = α \frac{s}{{∥s∥}^{1 / 2}} + k_{1} s - v

The time derivative of Lyapunov function (16) can be expressed as:

\begin{matrix} {\dot{V}}_{1} & = (2 β + \frac{α^{2}}{2} λ_{M}) \frac{s^{T} \dot{s}}{∥s∥} + (2 k_{2} + k_{1}^{2} λ_{M}) s^{T} \dot{s} + 2 v^{T} P \dot{v} + \frac{3}{2} α k_{1} λ_{M} \frac{s^{T} \dot{s}}{{∥s∥}^{1 / 2}} \\ - k_{1} (v^{T} P \dot{s} + s^{T} P \dot{v}) - α \frac{v^{T} P \dot{s}}{{∥s∥}^{1 / 2}} + \frac{α}{2} \frac{(v^{T} P s) (s^{T} \dot{s})}{{∥s∥}^{5 / 2}} - α \frac{s^{T} P \dot{v}}{{∥s∥}^{1 / 2}} . \end{matrix}

Lemma 2 ([37]).

As we have the diagonal matrix

P = P^{T} \geq 0

, the following inequality holds:

\frac{(v^{T} P s) (v^{T} P s)}{{∥s∥}^{5 / 2}} \leq \frac{v^{T} P^{2} v}{{∥s∥}^{1 / 2}} .

(18)

Substituting (14) and considering the inequality in Lemma 2, the equation becomes:

\begin{matrix} {\dot{V}}_{1} & \leq - \frac{α}{{∥s∥}^{1 / 2}} λ_{m} (β + λ_{M} \frac{α^{2}}{2}) ∥s∥ - λ_{m} (k_{1} β + 2 λ_{M} α^{2} k_{1}) ∥s∥ + \frac{2 α^{2}}{{∥s∥}^{1 / 2}} λ_{M}^{2} ∥v∥ {∥s∥}^{1 / 2} \\ - \frac{λ_{m}}{{∥s∥}^{1 / 2}} (\frac{5}{2} α k_{1}^{2} λ_{M} + α k_{2}) {∥s∥}^{2} - λ_{m} (k_{1} k_{2} + k_{1}^{3} λ_{M}) {∥s∥}^{2} + 2 k_{1}^{2} λ_{M}^{2} ∥s∥ ∥v∥ \\ + \frac{4 α k_{1} λ_{M}^{2}}{{∥s∥}^{1 / 2}} ∥s∥ ∥v∥ - k_{1} λ_{m}^{2} {∥v∥}^{2} - \frac{α}{2} λ_{m}^{2} \frac{{∥v∥}^{2}}{{∥s∥}^{1 / 2}} + λ_{M} (2 β + \frac{α^{2}}{2} λ_{M}) δ ∥s∥ \\ + λ_{M} (k_{1}^{2} λ_{M} + 2 k_{2}) δ {∥s∥}^{2} + \frac{3}{2} α k_{1} λ_{M}^{2} \frac{δ {∥s∥}^{2}}{{∥s∥}^{1 / 2}} + k_{1} λ_{M}^{2} δ ∥v∥ ∥s∥ + \frac{3}{2} α λ_{M}^{2} δ ∥v∥ {∥s∥}^{1 / 2} . \end{matrix}

Define a new vector

ζ = {[\begin{matrix} {∥s∥}^{1 / 2} & ∥s∥ & ∥v∥ \end{matrix}]}^{T}

, thus it can be obtained that when

ζ

converges to zero, s and v will also converge to zero at the same time.

Re–express the last inequality as:

{\dot{V}}_{1} \leq - \frac{1}{{∥s∥}^{1 / 2}} ζ^{T} Ω ζ - ζ^{T} Ψ ζ

(19)

where

Ω = [\begin{matrix} Ω_{11} & 0 & Ω_{13} \\ 0 & Ω_{22} & Ω_{23} \\ Ω_{13} & Ω_{23} & Ω_{33} \end{matrix}]

Ψ = [\begin{matrix} ψ_{11} & 0 & ψ_{13} \\ 0 & ψ_{22} & ψ_{23} \\ ψ_{13} & ψ_{23} & ψ_{33} \end{matrix}]

, with elements:

Ω_{11} = λ_{m} (α β + \frac{α^{3}}{2} λ_{M})

Ω_{13} = - α^{2} λ_{M}^{2}

Ω_{22} = λ_{m} (α k_{2} + \frac{5}{2} α k_{1}^{2} λ_{M}) - \frac{3}{2} α k_{1} λ_{M}^{2} δ

Ω_{23} = - 2 α k_{1} λ_{M}^{2}

Ω_{33} = \frac{α}{2} λ_{m}^{2}

and

ψ_{11} = λ_{m} (k_{1} β + 2 α^{2} k_{1} λ_{M}) - λ_{M} (2 β + \frac{α^{2}}{2} λ_{M}) δ

ψ_{13} = - \frac{3}{4} α λ_{M}^{2} δ

ψ_{22} = λ_{m} (k_{1} k_{2} + k_{1}^{2} λ_{M}) - λ_{M} (k_{1}^{2} λ_{M} + 2 k_{2}) δ

ψ_{23} = - k_{1}^{2} λ_{M}^{2} - \frac{k_{1}}{2} λ_{M}^{2} δ

ψ_{33} = k_{1} λ_{m}^{2}

The function

{\dot{V}}_{1}

will be negative definite if the matrices

Ω = Ω^{T} > 0

Ψ = Ψ^{T} > 0

. By selecting the parameters, it can be ensured that the diagonal elements and the determinants are positive. Therefore, the parameters should satisfy the following inequalities to ensure

Ω = Ω^{T} > 0

if:

\begin{matrix} α > 0, β > 4 α^{2} \frac{λ_{M}^{4}}{λ_{m}^{3}}, k_{1} > \frac{3}{5} \frac{λ_{M}}{λ_{m}} δ, k_{2} > 16 k_{1}^{2} \frac{λ_{M}^{4}}{λ_{m}^{3}}, \end{matrix}

(20)

and

Ψ = Ψ^{T} > 0

if:

\begin{matrix} α > 0, β > 0, k_{1} > max \{2 \frac{λ_{M}}{λ_{m}} δ, δ \sqrt{\frac{9}{14} \frac{λ_{M}^{3}}{λ_{m}^{3}}}\}, \\ k_{2} > \frac{3 k_{1}^{2} λ_{M}^{4} δ + \frac{1}{2} k_{1} λ_{M}^{4} δ^{2} {+ k}_{1}^{3} λ_{M} (2 λ_{M}^{3} - λ_{m}^{3})}{λ_{m}^{2} (k_{1} λ_{m} - 2 λ_{M} δ)} . \end{matrix}

(21)

Hence, the obtained matrices

Ω

Ψ

are symmetric positive definite when the given conditions (20) and (21) hold. Then we have:

{\dot{V}}_{1} \leq - \frac{λ_{min} (Ω)}{{∥s∥}^{1 / 2}} {∥ζ∥}^{2} - λ_{min} (Ψ) {∥ζ∥}^{2}

(22)

where

λ_{min} (Ω)

and

λ_{min} (Ψ)

are the smallest eigenvalues of the matrices

Ω

and

Ψ

, respectively. Thus, the candidate Lyapunov function in (16) can be rewritten as

V_{1} (s, v) = {\bar{ζ}}^{T} Q \bar{ζ}

, where

\bar{ζ} = {[\begin{matrix} \frac{s^{T}}{{∥s∥}^{1 / 2}} & s^{T} & v^{T} \end{matrix}]}^{T}

, and

Q = Q^{T} > 0

, given as:

Q = [\begin{matrix} 2 β + \frac{α^{2}}{2} λ_{M} & \frac{α}{2} k_{1} λ_{M} & - \frac{α P}{2} \\ \frac{α}{2} k_{1} λ_{M} & k_{2} + \frac{k_{1}^{2}}{2} λ_{M} & - \frac{k_{1} P}{2} \\ - \frac{α P}{2} & - \frac{k_{1} P}{2} & P \end{matrix}]

then the following inequality holds:

λ_{min} (Q) {∥\bar{ζ}∥}^{2} \leq V_{1} \leq λ_{max} (Q) {∥\bar{ζ}∥}^{2}

(23)

where

∥\bar{ζ}∥ = ∥ζ∥

, and

λ_{min} (Q), λ_{max} (Q) > 0

are the lowest and highest eigenvalues of the matrix Q, respectively.

Besides, consider the following fact:

{∥s∥}^{1 / 2} \leq ∥ζ∥ \leq \frac{{V_{1}}^{1 / 2}}{\sqrt{λ_{min} (Q)}} .

(24)

Then we have:

{\dot{V}}_{1} \leq - ρ_{1} {V_{1}}^{1 / 2} - ρ_{2} V_{1} \leq - ρ_{1} {V_{1}}^{1 / 2}

(25)

where

ρ_{1} = \frac{\sqrt{λ_{min} (Q)} λ_{min} (Ω)}{λ_{max} (Q)}

ρ_{2} = \frac{λ_{min} (Ψ)}{λ_{max} (Q)}

Since

α^{*}, β^{*}

are the upper bounds of

α, β

, it can be obtained that

\tilde{α} = α - α^{*} < 0

\tilde{β} = β - β^{*} < 0

. Then the time derivative of the Lyapunov function candidate (15) can be given as:

\begin{matrix} \dot{V} & = {\dot{V}}_{1} + \frac{1}{μ_{1}} \tilde{α} \dot{α} + \frac{1}{μ_{2}} \tilde{β} \dot{β} \\ \leq - ρ_{1} {V_{1}}^{1 / 2} - \frac{ω_{1}}{\sqrt{2 μ_{1}}} |\tilde{α}| - \frac{ω_{2}}{\sqrt{2 μ_{2}}} |\tilde{β}| + \frac{1}{μ_{1}} \tilde{α} \dot{α} + \frac{1}{μ_{2}} \tilde{β} \dot{β} + \frac{ω_{1}}{\sqrt{2 μ_{1}}} |\tilde{α}| + \frac{ω_{2}}{\sqrt{2 μ_{2}}} |\tilde{β}| . \end{matrix}

(26)

It should be noted that the following inequality holds:

{(x^{2} + y^{2} + z^{2})}^{1 / 2} \leq |x| + |y| + |z| .

(27)

Thus this leads to:

\dot{V} \leq - ρ_{3} V^{1 / 2} + \frac{1}{μ_{1}} \tilde{α} \dot{α} + \frac{1}{μ_{2}} \tilde{β} \dot{β} + \frac{ω_{1}}{\sqrt{2 μ_{1}}} |\tilde{α}| + \frac{ω_{2}}{\sqrt{2 μ_{2}}} |\tilde{β}|,

(28)

where

ρ_{3} = min \{ρ_{1}, ω_{1}, ω_{2}\}

. Then it can be written as:

\dot{V} \leq - ρ_{3} V^{1 / 2} + Θ,

(29)

where

Θ = - |\tilde{α}| (\frac{\dot{α}}{μ_{1}} - \frac{ω_{1}}{\sqrt{2 μ_{1}}}) - |\tilde{β}| (\frac{\dot{β}}{μ_{2}} - \frac{ω_{2}}{\sqrt{2 μ_{2}}})

Supposing that

α > α_{m}

and

|s_{i}| \geq ρ_{0}

, and we choose

Φ_{1} = μ_{1}

Φ_{2} = μ_{2}

, then we have:

\dot{α} = ω_{1} \sqrt{\frac{Φ_{1}}{2}}, \dot{β} = 2 k \dot{α} = ω_{2} \sqrt{\frac{Φ_{2}}{2}},

Θ = - |\tilde{α}| (\frac{\dot{α}}{μ_{1}} - \frac{ω_{1}}{\sqrt{2 μ_{1}}}) - |\tilde{β}| (\frac{\dot{β}}{μ_{2}} - \frac{ω_{2}}{\sqrt{2 μ_{2}}}) .

If we define

κ = \frac{ω_{2}}{2 ω_{1}} \sqrt{\frac{Φ_{2}}{Φ_{1}}} > 0

, we can obtain

Θ = 0

\dot{V} \leq - ρ_{3} V^{1 / 2} .

(30)

It should also be noted that the finite–time convergence can be ensured when the adaptive gains satisfy the inequalities (20) and (21). Thus the positive definiteness of the matrices

Ω

Ψ

and the convergence of whole processing system can be guaranteed.

Supposing

α > α_{m}

and

|s_{i}| < ρ_{0}

, the adaptive gains can be given as:

\dot{α} = - ω_{1} \sqrt{\frac{Φ_{1}}{2}}, \dot{β} = - ω_{2} \sqrt{\frac{Φ_{2}}{2}},

and we have:

Θ = 2 |α - α^{*}| \frac{ω_{1}}{\sqrt{2 μ_{1}}} + 2 |β - β^{*}| \frac{ω_{2}}{\sqrt{2 μ_{2}}} .

Since

α, β

are bounded by

α^{*}, β^{*}

, and

ω_{1}, ω_{2}, μ_{1}, μ_{2}

are chosen to be positive constants, then we have

Θ > 0

. According to [39], V is bounded which means all the signals have bounds. Therefore, the tracking errors can converge to a small neighborhood of the origin.

When

α \leq α_{m}

, the adaptive gains are:

\dot{α} = r_{0}, \dot{β} = 2 κ r_{0}, α = α_{0} + r_{0} t, β = β_{0} + 2 κ r_{0} t,

where

α_{0}, β_{0}

are the initial value of

α

and

β

, which are assumed to be 0 for simplicity.

Then we have the parameter

Θ

given as:

\begin{matrix} Θ & = - |\tilde{α}| (\frac{r_{0}}{μ_{1}} - \frac{ω_{1}}{\sqrt{2 μ_{1}}}) - |\tilde{β}| (\frac{2 κ r_{0}}{μ_{2}} - \frac{ω_{2}}{\sqrt{2 μ_{2}}}) \\ = - |r_{0} t - α^{*}| (\frac{r_{0}}{μ_{1}} - \frac{ω_{1}}{\sqrt{2 μ_{1}}}) - |2 κ r_{0} t - β^{*}| (\frac{2 κ r_{0}}{μ_{2}} - \frac{ω_{2}}{\sqrt{2 μ_{2}}}) . \end{matrix}

It can be observed that when

α \leq α_{m}

Θ

in (29) is negative. Hence,

\dot{V}

would be sign indefinite [40]. However,

α

will increase with time in this case. As soon as

α

increases over the small constant

α_{m}

, then (29) is valid and V starts decreasing. In order to accelerate the deceleration process, one can choose a larger

r_{0}

In conclusion, when

|s_{i}| > ρ_{0}

which means the sliding parameter

s_{i}

is far away from the sliding mode surface, it can reach the domain

|s_{i}| \leq ρ_{0}

within limited time. During the adaptation process,

s_{i}

may leave the domain in finite time due to the increase of

α

. However, it can always remain in a larger domain

|s_{i}| \leq {\bar{ρ}}_{0}

{\bar{ρ}}_{0} > ρ_{0}

. Thus, the whole control system can remain stable and bounded in finite time according to Lemma 1. Therefore, the proof of Theorem 1 is completed. □

4. Simulation Results

To demonstrate the effectiveness and advantage of our newly proposed control method, the simulations on a 2–DOF underwater robot manipulator, which is shown in Figure 2, are performed in this section. The rigid body dynamics of the manipulator system are taken from [31], and the underwater force

τ_{w} (q, \dot{q})

is taken from [32], thus the dynamics of the 2–DOF underwater robot manipulator in the form of (1) is given as:

M (q) \ddot{q} + C (q, \dot{q}) \dot{q} + G (q) + τ_{w} (q, \dot{q}) = τ_{e d i s} + τ

(31)

with matrices:

\begin{matrix} M (q) & = [\begin{matrix} (m_{1} + m_{2} {) L}_{1}^{2} + m_{2} L_{2}^{2} + 2 m_{2} L_{1} L_{2} c_{2} & m_{2} L_{2}^{2} + m_{2} L_{1} L_{2} c_{2} \\ m_{2} L_{2}^{2} + m_{2} L_{1} L_{2} c_{2} & m_{2} L_{2}^{2} \end{matrix}], \\ C (q, \dot{q}) & = [\begin{matrix} - m_{2} L_{1} L_{2} s_{2} {\dot{q}}_{2} & - m_{2} L_{1} L_{2} s_{2} ({\dot{q}}_{1} + {\dot{q}}_{2}) \\ m_{2} L_{1} L_{2} s_{2} {\dot{q}}_{1} & 0 \end{matrix}], \\ G (q) & = [\begin{matrix} (m_{1} + m_{2}) L_{1} g c_{2} + m_{2} L_{2} g {c_{1}}_{2} \\ m_{2} L_{2} g {c_{1}}_{2} \end{matrix}], \\ τ_{B} & = [\begin{matrix} τ_{B 1} \\ τ_{B 2} \end{matrix}], τ_{D} = [\begin{matrix} τ_{D 1} \\ τ_{D 2} \end{matrix}] \end{matrix}

in which

s_{2} = sin (q_{2})

c_{2} = cos (q_{2})

c_{12} = cos (q_{1} + q_{2})

;

m_{1}

m_{2}

are the masses;

L_{1}

L_{2}

are the lengths; and g is the gravity acceleration. The hydrodynamic forces

τ_{B}

and

τ_{D}

are given in Appendix A. In order to easily compare the performances of the proposed controller, traditional SM, NSM, and NFTSM control, which are referred to as controller 1, controller 2, and controller 3 respectively, are given in Appendix B. The parameters for the manipulator system and the model simulation are shown in Table 1.

The referred trajectory for both links is given as

q_{d} = {[\begin{matrix} q_{d 1} & q_{d 2} \end{matrix}]}^{T}

q_{d 1} = q_{d 2} = sin (0.5 t) + cos (0.15 t) - 1

. The time–varying lumped disturbances for the 2–DOF underwater robot manipulator are designed as

τ_{e d i s} = {[τ_{e d i s 1} τ_{e d i s 2}]}^{T}

with

τ_{e d i s 1} = τ_{e d i s 2} = 0.5 + 1.5 sin (t)

. Simulations are performed to compare the differences of the controllers in terms of the position precision, response speed, and the chattering phenomenon. The obtained simulation results are given in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7.

Figure 3 shows the trajectory tracking performances of four controllers for the 2–DOF underwater robot manipulator under the lumped disturbances. It can be seen that good tracking performance can be ensured by all four controllers, while the proposed control method approaches the desired trajectories with the highest reaching speed. Thanks to the utilization of ASTA, no noticeable chattering problem occurs in our proposed controller as shown in Figure 4. The performances in three comparing controllers encounter with the instantaneous jump in the control torques. However, the joint servo–motors can not reverse their rotation direction immediately which may cause failure or severe damage. In comparison, the proposed controller generates a smooth control input, thus the strongest robustness of our proposed controller is demonstrated.

Furthermore, in Figure 5, the fastest convergence and highest precision are obtained by our proposed controller. It can be seen that controller 1 shows the worst tracking precision when reaching the steady state. The steady tracking errors of controller 2 and 3 oscillate near the origin with a relatively higher amplitude than the proposed control. The estimations of the lumped disturbances are given in Figure 6, which nearly overlap the given disturbance signal. The compensation capability of the HOSTDO should be valued in the real applications since the actual currents and waves are impossible to calculate and proper estimations can significantly improve the disturbance–attenuation ability. The adaptation gains of the ASTA are shown in Figure 7, which release the limitation of the need for upper bound prior information of the lumped disturbances.

To further analyze the position tracking errors and energy consumption of the four controllers, the following integrated absolute error (IAE) [41], root mean square error (RMSE) [27], and energy of control input (ECI) [42] are taken into consideration, which are defined as:

\begin{matrix} {I A E}_{i} = \sum_{k = 1}^{N} |e_{i} (k)|; {R M S E}_{i} = \sqrt{\frac{1}{N} \sum_{k = 1}^{N} e_{i}^{2} (k)}; {E C I}_{i} = \sum_{k = 1}^{N} |τ_{i}^{2} (k)| \end{matrix}

(32)

where i is the link number; N is the number of samples; and

e_{i} (k)

τ_{i} (k)

are the trajectory tracking error and control input of link i at the

k^{t h}

sampling instant, respectively.

IAE and RMSE are used to evaluate the tracking performance while ECI is adopted to evaluate the control input energy consumption. Thus, the control target is to have a lower IAE, RMSE, and ECI so that better performances can be obtained with less energy consumption. The comparison results of the four controllers are shown in Table 2.

It should be noted that the ECI varies significantly due to the existence of the lumped disturbances. The ECI values of our proposed control method are only 50.2% and 52.7% for Controller 1, 45.9% and 47.6% for Controller 2, and 48.3% and 51.4% for Controller 3, respectively, which indicate that our proposed control consumes less input energy in the presence of the lumped disturbances. Furthermore, our proposed control still guarantees the best tracking precision in this case when considering the values of IAE and RMSE. Thus, considering the unknown lumped disturbances, our proposed control method can obtain the best tracking precision with the least control input energy consumption among the four controllers.

5. Conclusions

In this paper, a novel HOSTDO–based CNFT control with the ASTA method was proposed for underwater robot manipulators under time–varying lumped disturbances. The proposed controller utilized the NFTSM manifold and the ASTA method in the CNFT control scheme to ensure high control precision, singularity–free, chattering–free, and fast convergence. Meanwhile, the HOSTDO was used in this paper to estimate and compensate the lumped disturbances, which strengthen the robustness of our newly proposed controller. It should also be noted that the adaptation method was presented in the ASTA scheme so that prior information of the upper bound of the lumped disturbances was not needed. The stability of closed–loop system was demonstrated by the Lyapunov theory considering the ASTA dynamics and NFTSM surface. Finally, the effectiveness and robustness of our newly proposed control method were verified by comparing with the other three controllers on a 2–DOF underwater manipulator system. Although the proposed control is based on the underwater manipulators, it can also be utilized in some other robot system such as the manipulators on land, autonomous underwater vehicles, and so on.

In the proposed control method, the uncertain system modeling and the unknown external disturbances in the underwater manipulator system have been considered. However, it is still not clear how the robot interacts within underwater circumstances. In the future, we will intend to conduct underwater experiments to investigate the specific values of hydrodynamic coefficients, which will help us to build a more accurate mathematical dynamic model. To improve the application ability, the effects of time delays and actuator faults will be studied.

Author Contributions

Conceptualization, methodology, simulation, writing—original draft Z.Z.; conceptualization, funding acquisition, supervision, resources, writing—review & editing, G.T.; visualization, investigation, validation, R.X.; validation, writing—review & editing, L.H.; formal analysis, data processing, M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (Grant number 51979116), Wuhan Science and Technology Plan (Grant number 2020010602012052), the HUST Interdisciplinary Innovation Team Project and the Fundamental Research Funds for the Central Universities (Grant number HUST: 2020JYCXJJ063).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the anonymous reviewers for their valuable comments and suggestions that helped improve the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

It can be noted that the direction of the buoyancy torque

τ_{B}

stays opposite with that of the gravity

G (q)

. The buoyancy torque can be expressed as

τ_{B} = - (ρ_{f} / ρ_{m}) G (q)

, where

ρ_{f}

is the fluid density;

ρ_{m}

is the manipulator arm density.

By neglecting the tangential direction drag torque for slender cylinder rod which is assumed to be very small, the normal direction drag torque

τ_{D}

is considered in this paper. And the ith element of

τ_{D}

is given as [32]:

\begin{matrix} τ_{D i} & =^{i} {N_{i}}^{T}^{i} z_{i} \\ ^{i} N_{i} & =_{i + 1}^{i} R^{i + 1} N_{i + 1} +^{i} p_{i + 1} \times_{i + 1}^{i} R^{i + 1} f_{i + 1} \\ + \frac{1}{2} ρ_{f} C_{D} D_{i} {\int_{0}^{L_{i}} [\begin{matrix} x & 0 & 0 \end{matrix}]}^{T} \times^{i} v^{n} {(x)}_{i} ∥^{i} v^{n} {(x)}_{i}∥ d x \end{matrix}

(A1)

where

i = 1, 2

;

^{i} z_{i}

is the normal direction unit vector on ith link;

^{i} N_{i}

^{i + 1} N_{i + 1}

are the total drag torques on ith and

i + 1

th links, respectively;

_{i + 1}^{i} R

is the transformation matrix toward ith link from

i + 1

th link;

^{i} p_{i + 1}

is the directional vector toward ith link from

i + 1

th link;

^{i + 1} f_{i + 1}

is the drag force on

i + 1

th link;

D_{1}

D_{2}

describe the diameters of link 1 and link 2, respectively;

C_{D}

denotes the drag coefficient; and

^{i} v^{n} {(x)}_{i}

is the normal directional velocity vector of ith link.

Appendix B

Three controllers are simulated in this section for comparison. The first one is our newly proposed one (Proposed) and the second one is the conventional SM controller (Controller 1) which is given as [43]:

\begin{matrix} u_{1} = M_{0} (x_{1}) [- f (x_{1}, x_{2}) - σ_{0} \dot{e} + {\ddot{q}}_{d}] - M_{0} (x_{1}) [{k^{'}}_{0} s + ({γ^{'}}_{0} + {η^{'}}_{0}) sgn (s)] \end{matrix}

(A2)

with the SM surface as

s = \dot{e} + σ_{0} e

, where

σ_{0} > 0

k_{0}^{'}, γ_{0}^{'}, η_{0}^{'} > 0

are the designed parameters. And the third one is the NSM controller (Controller 2) as [44]:

\begin{matrix} u_{2} = M_{0} (x_{1}) [- f (x_{1}, x_{2}) - \frac{p^{'}}{σ_{0}^{'}} l^{'} {\dot{e}}^{[2 - l^{'} / p^{'}]} + {\ddot{q}}_{d}] - M_{0} (x_{1}) [k_{0}^{″} s + (γ_{0}^{″} + η_{0}^{″}) sgn (s)] \end{matrix}

(A3)

with the NSM manifold

s = e + σ_{0}^{'} {\dot{e}}^{[l^{'} / p^{'}]}

, where

σ_{0}^{'} > 0

k_{0}^{″}, γ_{0}^{″}, η_{0}^{″} > 0

and

1 < l^{'} / p^{'} < 2

are the known constants. It should be noted that the FTSM type reaching law was designed in [44], which was replaced by the normal exponential reaching law in this paper for a fair comparison. The last one (Controller 3) uses the NFTSM surface (15) for comparison, which is given as follows:

u_{0} = u_{e q 0} + u_{r e 0}

(A4)

where

u_{e q 0} = M_{0} (x_{1}) [\begin{matrix} - \frac{1}{σ_{2}} \frac{p}{l} ({\dot{e}}^{[2 - l / p]} + σ_{1} φ {|e|}^{φ - 1} \cdot {\dot{e}}^{[2 - l / p]}) \\ - f (x_{1}, x_{2}) + {\ddot{q}}_{d} \end{matrix}]

(A5)

is used to control the nominal parts, and

u_{r e 0} = - M_{0} (x_{1}) [k_{0} s + (γ_{0} + η_{0}) s g n (s)]

(A6)

is used to compensate for the lumped disturbances, where

γ_{0}

is the upper bound of the lumped disturbances,

k_{0} > 0

, and

η_{0}

is a small positive constant.

References

Qiao, L.; Zhang, W. Double–loop integral terminal sliding mode tracking control for UUVs with adaptive dynamic compensation of uncertainties and disturbances. IEEE J. Ocean. Eng. 2019, 44, 29–53. [Google Scholar] [CrossRef]
Leabourne, K.N.; Rock, S.M. Model development of an underwater manipulator for coordinated arm–vehicle control. In Proceedings of the OCEANS’98, Nice, France, 28 September–1 October 1998; Volume 2, pp. 941–946. [Google Scholar]
Ioi, K.; Itoh, K. Modelling and simulation of an underwater manipulator. Adv. Robot. 1989, 4, 303–317. [Google Scholar] [CrossRef]
Zhang, M.; Liu, X.; Tian, Y. Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single–DOF Manipulator. J. Mar. Sci. Eng. 2019, 7, 261. [Google Scholar] [CrossRef] [Green Version]
Wang, T.; You, Z.; Song, W.; Zhu, S. Dynamic Analysis of an Underwater Cable–Driven Manipulator with a Fluid–Power Buoyancy Regulation System. Micromachines 2020, 11, 1042. [Google Scholar] [CrossRef] [PubMed]
Tang, J.; Zhang, Y.; Huang, F. Design and Kinematic Control of the Cable–Driven Hyper–Redundant Manipulator for Potential Underwater Applications. Appl. Sci. 2019, 9, 1142. [Google Scholar] [CrossRef] [Green Version]
Yao, J.; Wang, C. Model reference adaptive control for a hydraulic underwater manipulator. J. Vib. Control 2012, 18, 893–902. [Google Scholar] [CrossRef]
Barbalata, C.; Dunnigan, M.W.; Petillot, Y. Coupled and Decoupled Force/Motion Controllers for an Underwater Vehicle–Manipulator System. J. Mar. Sci. Eng. 2018, 6, 96. [Google Scholar] [CrossRef] [Green Version]
Zhou, Z.; Tang, G.; Huang, H.; Han, L.; Xu, R. Adaptive nonsingular fast terminal sliding mode control for underwater manipulator robotics with asymmetric saturation actuators. Control Theory Technol. 2020, 18, 81–91. [Google Scholar] [CrossRef]
Yang, C.; Yao, F.; Zhang, M.; Zhang, Z.; Wu, Z.; Dan, P. Adaptive Sliding Mode PID Control for Underwater Manipulator Based on Legendre Polynomial Function Approximation and Its Experimental Evaluation. Appl. Sci. 2020, 10, 1728. [Google Scholar] [CrossRef] [Green Version]
Xu, R.; Tang, G.; Han, L.; Huang, H.; Xie, D. Robust finite–time attitude tracking control of a CMG–based AUV with unknown disturbances and input saturation. IEEE Access 2019, 7, 56409–56422. [Google Scholar] [CrossRef]
Le, T.H.; Thanh, H.L.N.N.; Huynh, T.T.; Van, M.; Hoang, Q.D.; Do, T.D. Robust Position Control of an Over–actuated Underwater Vehicle under Model Uncertainties and Ocean Current Effects Using Dynamic Sliding Mode Surface and Optimal Allocation Control. Sensors 2021, 21, 747. [Google Scholar]
Thanh, H.L.N.N.; Vu, M.T.; Mung, N.X.; Nguyen, N.P.; Phuong, N.T. Perturbation Observer–Based Robust Control Using a Multiple Sliding Surfaces for Nonlinear Systems with Influences of Matched and Unmatched Uncertainties. Mathematics 2020, 8, 1371. [Google Scholar] [CrossRef]
Vu, M.T.; Thang, H.T.T.L.; Thang, Q.; Duc, T.; Hoang, Q.D.; Le, T.H. Station–Keeping Control of a Hovering Over–Actuated Autonomous Underwater Vehicle Under Ocean Current Effects and Model Uncertainties in Horizontal Plane. IEEE Access 2021, 9, 6855–6867. [Google Scholar] [CrossRef]
Yang, L.; Yang, J. Nonsingular fast terminal sliding–mode control for nonlinear dynamical systems. Int. J. Robust Nonlinear Control 2011, 21, 1865–1879. [Google Scholar] [CrossRef]
Zhao, L.; Jia, Y. Decentralized adaptive attitude synchronization control for spacecraft formation using nonsingular fast terminal sliding mode. Nonlinear Dyn. 2014, 78, 2779–2794. [Google Scholar] [CrossRef]
Xu, S.S.; Chen, C.; Wu, Z. Study of nonsingular fast terminal sliding–mode fault–tolerant control. IEEE Trans. Ind. Electron. 2015, 62, 3906–3913. [Google Scholar] [CrossRef]
Zhang, L.; Liu, L.; Wang, Z.; Xia, Y. Continuous finite–time control for uncertain robot manipulators with integral sliding mode. IET Control Theory Appl. 2018, 12, 1621–1627. [Google Scholar] [CrossRef]
Van, M.; Kang, H.J.; Suh, Y.S. Second order sliding mode–based output feedback tracking control for uncertain robot manipulators. Int. J. Adv. Robot. Syst. 2013, 10, 16. [Google Scholar] [CrossRef] [Green Version]
Levant, A.; Alelishvili, L. Integral high–order sliding modes. IEEE Trans. Autom. Control 2007, 52, 1278–1282. [Google Scholar] [CrossRef] [Green Version]
Li, X.; Yu, X.; Han, Q. Stability analysis of second–order sliding mode control systems with input–delay using poincaré map. IEEE Trans. Autom. Control 2013, 58, 2410–2415. [Google Scholar] [CrossRef]
Li, J.; Zhang, Q.; Yan, X.; Spurgeon, S.K. Integral sliding mode control for Markovian jump T–S fuzzy descriptor systems based on the super–twisting algorithm. IET Control Theory Appl. 2017, 11, 1134–1143. [Google Scholar] [CrossRef] [Green Version]
Zhao, Z.; Gu, H.; Zhang, J.; Ding, G. Terminal sliding mode control based on super–twisting algorithm. J. Syst. Eng. Electron. 2017, 28, 145–150. [Google Scholar] [CrossRef]
Rehman, F.U.; Mufti, M.R.; Din, S.U.; Afzal, H.; Qureshi, M.I.; Khan, D.M. Adaptive Smooth Super–Twisting Sliding Mode Control of Nonlinear Systems With Unmatched Uncertainty. IEEE Access 2020, 8, 177932–177940. [Google Scholar] [CrossRef]
Shtessel, Y.; Taleb, M.; Plestan, F. A novel adaptive–gain supertwisting sliding mode controller: Methodology and application. Automatica 2012, 48, 759–769. [Google Scholar] [CrossRef]
Van, M.; Ge, S.S.; Ren, H. Finite time fault tolerant control for robot manipulators using time delay estimation and continuous nonsingular fast terminal sliding mode control. IEEE T. Cybern. 2017, 47, 1681–1693. [Google Scholar] [CrossRef]
Wang, Y.; Yan, F.; Ju, F.; Chen, B.; Wu, H. Optimal nonsingular terminal sliding mode control of cable–driven manipulators using super–twisting algorithm and time–delay estimation. IEEE Access 2018, 6, 61039–61049. [Google Scholar] [CrossRef]
Lin, C.; Sun, S.; Walker, P.; Zhang, N. Accelerated adaptive second order super–twisting sliding mode observer. IEEE Access 2019, 7, 25232–25238. [Google Scholar] [CrossRef]
Moreno, J.A.; Osorio, M. A Lyapunov approach to second–order sliding mode controllers and observers. In Proceedings of the Decision and Control, Cancun, Mexico, 9–11 December 2008; pp. 2856–2861. [Google Scholar]
Kumar, P.R.; Behera, A.K. Bandyopadhyay, B. Robust finite–time tracking of Stewart platform: A super–twisting like observer–based forward kinematics solution. IEEE Trans. Ind. Electron. 2017, 64, 3776–3785. [Google Scholar] [CrossRef]
Cuong, P.V.; Nan, W.Y. Adaptive trajectory tracking neural network control with robust compensator for robot manipulators. Neural Comput. Appl. 2016, 27, 525–536. [Google Scholar] [CrossRef]
Jia, P. Reach on Grasp Planning and Optimization of Grasping Forces for Robotic Dexterous Hands. Ph.D. Thesis, Mechanical and Electronic Engineering, Harbin Engineering University, Harbin, China, 2011. (In Chinese). [Google Scholar]
Hong, Y.; Huang, J.; Xu, Y. On an output feedback finite–time stabilization problem. IEEE Trans. Autom. Control 2001, 46, 305–309. [Google Scholar] [CrossRef] [Green Version]
Angulo, M.T.; Moreno, J.A.; Fridman, L. Robust exact uniformly convergent arbitrary order differentiator. Automatica 2013, 49, 2489–2495. [Google Scholar] [CrossRef]
Han, Y.; Liu, X. Continuous higher–order sliding mode control with time–varying gain for a class of uncertain nonlinear systems. ISA Trans. 2016, 62, 193–201. [Google Scholar] [CrossRef] [PubMed]
Capisani, L.M.; Ferrara, A.; Ferreira, A.; Fridman, L. Higher order sliding mode observers for actuator faults diagnosis in robot manipulators. In Proceedings of the IEEE ISIE, Bari, Italy, 4–7 July 2010; pp. 2103–2108. [Google Scholar]
Vidal, P.V.N.M.; Nunes, E.V.L.; Hsu, L. Output–feedback multivariable global variable gain super–twisting algorithm. IEEE Trans. Autom. Control 2017, 62, 2999–3005. [Google Scholar] [CrossRef]
Wang, Y.; Yan, F.; Chen, J.; Chen, B. Continuous nonsingular fast terminal sliding mode control of cable–driven manipulators with super–twisting algorithm. IEEE Access 2018, 6, 49626–49636. [Google Scholar] [CrossRef]
Song, H.; Zhang, T. Fast Robust Integrated Guidance and Control Design of Interceptors. IEEE Trans. Control. Syst. Technol. 2016, 24, 349–356. [Google Scholar] [CrossRef]
Plestan, F.; Shtessel, Y.; Brégeault, V.; Poznyak, A. New methodologies for adaptive sliding mode control. Int. J. Control 2010, 83, 1907–1919. [Google Scholar] [CrossRef] [Green Version]
Mondal, S.; Mahanta, C. Adaptive second order terminal sliding mode controller for robotic manipulators. J. Franklin Inst. Eng. Appl. Math. 2014, 351, 2356–2377. [Google Scholar] [CrossRef]
Skogestad, S. Simple analytic rules for model reduction and PID controller tuning. J. Process Control 2003, 13, 291–309. [Google Scholar] [CrossRef] [Green Version]
Islam, S.; Liu, X.P. Robust sliding mode control for robot manipulators. IEEE Trans. Ind. Electron. 2011, 58, 2444–2453. [Google Scholar] [CrossRef]
Yu, S.; Yu, X.; Shirinzadeh, B.; Man, Z. Continuous finite–time control for robotic manipulators with terminal sliding mode. Automatica 2005, 41, 1957–1964. [Google Scholar] [CrossRef]

Figure 1. Block diagram of the proposed control method.

Figure 2. Two–link underwater robot manipulator control system.

Figure 3. Position tracking performance with disturbances: (a) Link 1. (b) Link 2.

Figure 4. Control inputs with disturbances: (a) Link 1. (b) Link 2.

Figure 5. Tracking errors with disturbances: (a) Link 1. (b) Link 2.

Figure 6. Disturbance estimations: (a) Link 1. (b) Link 2.

Figure 7. Adaptation gains with disturbances: (a)

α

. (b)

β

Figure 7. Adaptation gains with disturbances: (a)

α

. (b)

β

Table 1. Parameters of the manipulator system and the model simulation. Sliding mode (SM); adaptive super twisting algorithm (ASTA); higher–order super–twisting disturbance observer (HOSTDO).

Parameter	Value	Parameter	Value	Parameter	Value
Parameters of the manipulator system
$m_{1}$	3.39 kg	$D_{1}$	0.04 m	$L_{1}$	1 m
$m_{2}$	3.39 kg	$D_{2}$	0.04 m	$L_{2}$	1 m
$ρ_{f}$	1000 kg/m $^{3}$	$C_{D}$	0.6	g	9.8 m/s $^{2}$
$ρ_{m}$	2700 kg/m $^{3}$
Parameters of the model simulation comparison
Parameters of SM surfaces		$σ_{0}$	1	$σ_{0}^{'}$	1
		$σ_{1}$	1	$σ_{2}$	1
		l	7	$l^{'}$	7
		p	9	$p^{'}$	9
		$φ$	1.3
Gains of SM reaching laws		$k_{0}$	20	$k_{1}^{'}$	20
		$k_{1}^{″}$	20	$η_{0}$	0.1
		$η_{0}^{'}$	0.1	$η_{0}^{″}$	0.1
		$γ_{0}$	2	$γ_{0}^{'}$	2
		$γ_{0}^{″}$	2	$η_{0}$	0.1
Parameters of the ASTA		$k_{1}$	15	$k_{2}$	0.1
		$ω_{1}$	5	$Φ_{1}$	2
		$α_{m}$	0.5	$r_{0}$	0.5
		$ρ_{0}$	0.05	$κ$	0.5
Parameters of the HOSTDO		$ξ_{1}$	8	$ξ_{2}$	22
Parameters of the HOSTDO		$ξ_{3}$	8	$ξ_{4}$	2

Table 2. Tracking performance comparisons.

Type of Controller	Link	IAE	RMSE	ECI
Controller 1	1	18.120	0.022718	4.1162 $\times 10^{6}$
Controller 1	2	15.813	0.021560	7.4115 $\times 10^{5}$
Controller 2	1	7.7463	0.016877	4.5057 $\times 10^{6}$
Controller 2	2	7.7831	0.016592	8.1977 $\times 10^{5}$
Controller 3	1	5.7875	0.014082	4.2817 $\times 10^{6}$
Controller 3	2	5.5934	0.013763	7.5878 $\times 10^{5}$
Proposed	1	3.4521	0.011636	2.0673 $\times 10^{6}$
Proposed	2	5.2355	0.012449	3.9025 $\times 10^{5}$

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhou, Z.; Tang, G.; Xu, R.; Han, L.; Cheng, M. A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators. J. Mar. Sci. Eng. 2021, 9, 269. https://doi.org/10.3390/jmse9030269

AMA Style

Zhou Z, Tang G, Xu R, Han L, Cheng M. A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators. Journal of Marine Science and Engineering. 2021; 9(3):269. https://doi.org/10.3390/jmse9030269

Chicago/Turabian Style

Zhou, Zengcheng, Guoyuan Tang, Ruikun Xu, Lijun Han, and Maolin Cheng. 2021. "A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators" Journal of Marine Science and Engineering 9, no. 3: 269. https://doi.org/10.3390/jmse9030269

APA Style

Zhou, Z., Tang, G., Xu, R., Han, L., & Cheng, M. (2021). A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators. Journal of Marine Science and Engineering, 9(3), 269. https://doi.org/10.3390/jmse9030269

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

A Novel Continuous Nonsingular Finite–Time Control for Underwater Robot Manipulators

Abstract

1. Introduction

2. Problem Formulation

2.1. Underwater Robot Manipulator Model and Properties

2.2. Fundamental Facts

3. The Design of CNFT Control Based on HOSTDO and ASTA

3.1. The Design of HOSTDO

3.2. HOSTDO Based CNFT Control Design with ASTA Method and Time–Varying Gain Matrix

4. Simulation Results

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

Appendix A

Appendix B

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI