Purpose

This paper proposes a novel non-singular fast terminal sliding mode (NFTSM) control approach to address challenges in autonomous underwater vehicles (AUVs) caused by external disturbances and modeling uncertainties. The goal is to achieve finite-time convergence and avoid singular problems in control systems.

Design/methodology/approach

A NFTSM control is designed to ensure finite-time stability. A control law is derived based on the sliding mode control rate obtained from the sliding mode surface. Lyapunov stability analysis is conducted to verify the convergence of tracking errors to a range near zero within a finite time period.

Findings

Simulation results demonstrate the effectiveness of the proposed control scheme. The approach achieves faster convergence and improved performance than two other finite-time tracking controllers, validating its suitability for AUV applications.

Originality/value

The proposed control scheme introduces a NFTSM approach, which avoids singular problems and ensures finite-time convergence. This innovation enhances the robustness and efficiency of control systems for AUVs, offering significant contributions to underwater vehicle control research.

Autonomous underwater vehicles (AUVs) have become vital and cost-effective tools in various marine applications, such as underwater mining, underwater rescue and marine research (Li et al., 2023; Chen et al., 2023; Zhu et al., 2023). Their significant importance has been greatly enhanced by their outstanding capability to navigate perilous oceanic environments and gather dependable, long-term data. To guarantee the successful execution of AUV tasks and missions, it is essential to provide a tracking control system that offers an optimal tracking performance (Karkoub et al., 2017).

To date, various control strategies have been utilized for AUVs to guarantee efficient mission execution despite external disturbances (Yan et al., 2023; Bessa et al., 2010; Guerrero et al., 2023; Wen et al., 2023; Wang et al., 2023a; González-García et al., 2020; Liu et al., 2023). Among these control schemes, siding mode control has been extensively studied for AUV control design because of its robustness, simplicity and insensitivity to parameter variations, as mentioned in references Qin et al. (2023) and Wang et al. (2023b). Linear sliding mode control is commonly used in the development of control schemes for AUV owing to its simplicity. An example of a proposed trajectory tracking control scheme for underwater vehicles based on the linear sliding mode surface, where external disturbances were approximated using an adaptive Takagi–Sugeno–Kang fuzzy model (Li et al., 2018). However, these linear sliding mode control schemes only ensure asymmetrical stability, not finite-time stability. To achieve improved tracking performance for AUV, terminal sliding mode (TSM) control is employed to design the control in this field owing to its ability to accelerate the convergence speed near the equilibrium point. An innovative terminal sliding mode control was proposed to effectively tackle the trajectory tracking challenge of AUV in the horizontal plane (Elmokadem et al., 2017). Yu and colleagues Yu and Zhihong (2002) introduced fast terminal sliding mode (FTSM) control, which demonstrates faster convergence speeds than TSM when the initial position is either distant from or close to the equilibrium point. However, upon analyzing the control design, it was discovered that both TSM and FTSM face the issue of singularity. To solve this challenge, a novel concept called nonsingular terminal sliding mode (NTSM) was implemented to eliminate singularity without requiring any additional procedures (Feng et al., 2002). Furthermore, Elmokadem et al. (2017), Hou et al. (2020) and Yang and Yang (2011) introduced a nonsingular fast terminal sliding mode (NFTSM) control that adapts the terminal sliding mode surface to tackle singularities.

Apart from singularity avoidance, chattering reduction is also a crucial concern in the design of sliding mode control (Li et al., 2024). Several algorithms have been developed to reduce chattering. The boundary layer technique, combined with a saturation function, is commonly employed to suppress the effects of chattering.

In this paper, we propose a novel finite-time NFTSM control, along with its corresponding control law. This control scheme effectively minimizes the impact of chattering and mitigates the influence of uncertain internal parameters and external disturbances on AUV pose convergence. By utilizing the Lyapunov function, we illustrate that the proposed NFTSM guarantees the finite-time convergence of the AUV system while ensuring a nonsingular controller.

The remainder of this paper is organized as follows. In the next section, a brief overview of the dynamic AUV model is provided, followed by a presentation of the problem statement. Section 3 presents a finite-time tracking control law for an AUV, utilizing NFTSM control. The proposed control scheme is validated through the simulation results, presented in Section 4, followed by concluding remarks and potential avenues for future research.

The main features and contributions of this article are summarized as follows:

  1. This paper proposes a nonsingular fast terminal sliding mode control with a nonsingular sliding mode surface. The key is to restructure the sliding mode surface function so that its derivatives no longer contain negative fractional powers. This ensures that the control input is smooth and bounded throughout the state space, especially near the equilibrium point.

  2. A novel sliding mode control method is proposed for an underwater vehicle system. This method takes into account the characteristics of the underwater vehicle dynamic model and designs a controller that ensures finite-time convergence.

  3. An adaptive coefficient is proposed that dynamically adjusts the switching gain in the control input, reducing system jitter and eliminating its interference.

  4. The proposed nonsingular fast terminal sliding mode control is applied to a six-degree-of-freedom underwater robot system. Compared with previously proposed control methods, the proposed method achieves better control performance.

This section first explains the AUV model, followed by an explanation of the relevant lemmas, consistencies and definitions.

A general description of the six-DOF nonlinear equation of the AUV motion is described as follows (Geranmehr and Nekoo, 2015):

(1)

where η ∈ R6×1 = [η1,η2]T = [x,y,z,ϕ,θ,ψ]T represents the position and attitude of AUV in the earth frame; ν ∈ R6×1 = [ν1,ν2]T=[u,v,w,p,q,r]T represents the velocity and angular velocity in the body-fixed frame; J(η)is the transformation matrix between the inertial and body-fixed frames, and it is assumed that the spatial transformation matrix and its inverse exist, which is reasonable for many AUV applications (Fischer et al., 2014), especially for the fully actuated AUV; τR6×1=[Fe,Me]T=[X,Y,Z,K,M,N]T represents the control forces and moments, which includes surge thrust X, sway thrust Y, heave thrust Z, roll thrust K, pitch thrust M, yaw thrust N.M ∈ R6×6 is the inertia matrix, C(ν) ∈ R6×6 represents the centripetal and Coriolis terms for the rigid-body and hydrodynamics, D(ν) is the hydrodynamic damping matrix, G(η) is the gravity and buoyancy forces and moments and d(η, ν, t) represents the external disturbances acting on AUV.

The position and attitude of AUV are usually described in Earth's frame. The transformation matrix between the body-fixed and earth frames is represented by Guerrero et al. (2024).

(2)

Where

Remark 1.

The pitch angle θ is bounded by |θ|<π2 to avoid singularity problems in J(η) owing to rotation (Qiao et al., 2017).

Remark 2.

In practical systems, uncertainties d(η, ν, t) such as hydrodynamic coefficients, external disturbances and modeling errors, result in inaccurate system parameters that is d(η,ν,t)=M̃(ν)ν̇+C̃(ν)ν+D̃(ν)ν+G̃(η)+dis, where M̃(ν),C̃(ν),D̃(ν),G̃(η) represent the hydrodynamic coefficients and modeling errors, and dis represents the external disturbances. In practical systems, the velocities, accelerations and uncertain parameters are all bounded, |dn(η, ν, t)| ≤ Ln, n = 1, 2, 3, 4, 5, 6, where Ln is the upper bound of the n dimension.

Remark 3.

τ represents the thrust and moments of the propellers decomposed into the six degrees of freedom direction, and it satisfies a matrix relation with the actual propeller forces of the AUV model.

The two equations in Eq.(1) can be expressed as follows:

(3)

where M(η) = MJ1, C(η) = −MJ2J−1(η) + C(ν)J−1(η) + D(ν)J−1(η),

To conveniently design a finite-time tracking controller, Eq.(3) can be rewritten in a concise form:

(4)

Where, A=(MJ1)1(MJ2J1(η)C(ν)J1(η)D(ν)J1(η)), B=(MJ1)1,D=(MJ1)1d(η,t), f=Aη̇BG(η).

The objective of this study is to propose an NFTSM and provide a control law based on a nonsingular sliding mode surface for an AUV to guarantee that the tracking error converges to a small region near zero in finite time.

The objective of this paper is to propose an NFTSM and provide a control law based on a nonsingular sliding mode surface for an AUV to guarantee that the tracking error converges to a small region near zero in a finite time. This section introduces the symbol representation, related lemmas and definitions that will be used in the subsequent analysis and proof.

Definition 1.

Hu and Jiang (2018), given a sliding surface s, moving on the sliding surface is called an ideal sliding mode, which can be described as:

Definition 2.

Hu and Jiang (2018), given a sliding surface s, moving within the δ range near the sliding surface is called an actual sliding mode, which can be described as:

Lemma: Shen and Huang (2009) and Plestan et al. (2010) assume that there exists a continuously differentiable and positive definite function V(x), which satisfies:

(5)

where c > 0, b > 0, 0 < α < 1. Then V(x) reaches V(x) = 0 in finite time and the stable time Tx(x0) satisfies

(6)

The FTSM is described by the following differential equation (Mishra et al., 2016):

(7)

where: e=ηηd,ė=η̇η̇d,ηd and η̇d are the desired position vector and desired velocity vector, respectively, m = diag(m1, , m6), n = diag(n1, …, n6), mi, ni > 0, i = 1, , 6, p and q are positive odd integers, which satisfy the following condition: p > q.

When the system state variable reaches the sliding mode surface (s = 0), for properly chosen q and p, given an initial state e(0) ≠ 0, dynamics (7) will reach e = 0 in finite time, as shown in Eq.(8). The physical interpretation is that when e is far away from zero, the approximate dynamics become ė=me whose fast convergence when far away from zero is well understood. When close to e = 0, the approximate dynamics become ė=neqp which is a terminal attractor.

From any initial state e(0) ≠ 0, the system state converges to the equilibrium point in time (Mishra et al., 2016):

(8)

To make dynamical model (4) reach equilibrium quickly along the given surface, the control τ is chosen as:

(9)

where γ1, γ2, m, n and qp represent the parameter diagonal matrices and all diagonal parameters are positive. When ė0 and e = 0, eqppė becomes extremely large for the control input τ that, leading to a singular problem. In ideal sliding mode control, as defined in Definition 1, as long as the system enters the sliding mode (s = 0) and q < p < 2p, then eqppė=meqpne2qpp and no singular problem occurs.

In the control design of FTSMs, a singularity problem often exists. This is because the designed controllers have nonlinear terms eα (α < 0), which results in unboundness of the control input when x diminishes to zero (Xu et al., 2015). To overcome this difficulty, the NFTSM is proposed as follows:

Definition NFTSM is described as

(10)

where: α = diag(α1, , α6), β = diag(β1, , β6) and αi, βi > 0, i = 1, , 6.0 < qi < pi < 2qi for represent internal parameters, which are odd integers,

Based on the requirement of finite-time stability, the control law is designed as:

(11)

where: α, β, p and q are defined in Equation (10). σ=[σ1,,σ6]T represents internal parameters with σi > 0, |Di| < Li, i = 1, , 6, where Li represents the upper bound of D in the i-th dimension.

This subsection demonstrates the finite-time stability of both the sliding mode surface and tracking error e through Lyapunov stability analysis.

Theorem: For system (4), when the sliding surface is selected as (10) and its control law τ is expressed in (11), the system trajectory rapidly converge to zero within a finite time. In addition, the occurrence of singularity is guaranteed to be avoided throughout the entire process.

Proof: Consider the Lyapunov function candidate as

(12)

whose time derivative is

Let k=pq(1β(ė+αe))pqq then

(13)

Because 0 < p < 2q and p is an odd integer, p − q and p + q are even integers,and (1β(ė+αe))pqq>0,which implies k > 0 and sq+pp. Therefore, when V̇0, the condition for Lyapunov stability is satisfied.

Eq.(13) can be rewritten as

(14)

The Eq.(14) is further reformulated as

where α=2α,k=2k and α>0,kσ>0,q+p2p(0,1).This inequality has the same form as Lemma 1, and the coefficients satisfy the corresponding coefficient requirements of Lemma 1. This implies that the NFTSM converges in finite time. Once the system enters the sliding mode manifold, we have e+(1β(ė+αe))pq=0,and the tracking error of the AUV converges to the equilibrium point in the same time as in Equation (6). This indicates that the tracking error converges in a finite time. Therefore, Theorem 1 is proved.

Remark 4.

The NFTSM can only remain nonsingular when 0 < q < p < 2q, p and q are odd integers.

Remark 5.

Owing to the model uncertainties and time-varying external disturbances, the sign( ) function in the control law is used to counteract the model uncertainties and external disturbances. As the tracking error approaches the equilibrium point, parameter k, which is related to the tracking error, tends to zero. This weakens the impact of (L|s| − sD(x, t)), effectively suppressing the chattering. If we assume that the system can achieve the ideal sliding mode manifold as per Definition 1, k becomes zero when the system state variable reaches the equilibrium point. Consequently, the presence of model uncertainties and external disturbances will have no impact on the system, resulting in a chattering-free sliding mode control.

To analyze the convergence of the tracking error, consider the following Lyapunov function:

(15)

from Eq.(8),when the system reaches the equilibrium point, i.e s = 0, it lead to

(16)

Eq.(16) is reorganized as follows:

(17)

The time derivative of Eq.(15) is

It is evident from the above that V2 ≥ 0 and V̇20, therefore, it is guaranteed that e asymptotically converges to zero, that is, η = ηd. Thus, it can be concluded that the convergence of the AUV’s state to the desired state ensures the asymptotic convergence of the position tracking errors e to zero.

In order to verify the effectiveness of the proposed control scheme, this paper established a kinematic and dynamic model in Matlab, added control input and disturbance and used the Euler method to update the model iteratively. In order to more convincingly illustrate the performances of the proposed finite-time tracking control scheme, two additional finite-time tracking control scheme are employed for comparison: classical NFTSM control (Hou et al., 2020) and nonsingular terminal sliding mode control (NTSMC; Mu and He, 2018). The simulation results, were recorded as [18] and [29].

In these simulations, the initial states of the system are set as η0=[x0,y0,z0,ϕ0,θ0,ψ0]T = [10,5,5,1,1,0]T, The desired trajectory in the earth-fixed frame is described by

where xd = 4(1 − cost); yd = 2sint; zd = −0.2 t; ϕd = 0; θd = 0 and ψd = 90. The relevant parameters of the control method proposed in this paper are as follows: α = diag(10, 10, 10, 1, 1, 1), β = diag(10, 10, 10, 1, 1, 1) and q = [3,3,3,3,3,3]T, p = [5,5,5,5,5,5]T. The relevant parameters of the control method proposed in Hou et al. (2020) are as follows: α=diag(0.1,0.1,0.1,0.1,0.1,0.1),β=diag(1,1,1,1,1,1),K1=diag(2,2,2,2,2,2)T,K2=diag(1,1,1,1,1,1),γ2=[53,53,53,53,53,53]T,γ1=[2,2,2,2,2,2]T. The relevant parameters of the control method proposed in [29] are as follows: β = diag(10, 10, 10, 1, 1, 1), q = [3,3,3,3,3,3]T, p = [5,5,5,5,5,5]T and D accounts for external disturbances and internal parameter uncertainties.

D = [2sin(200t), 2sin(t), − 0.02 t, 0.2sin(200t), 0.2sin(200t), 0.2sin(200t)].

In terms of parameter selection, the singular fast terminal sliding mode control method, expressed as Equation (7), and the NFTSM control proposed in this paper, expressed as Equation (10), have the same convergence characteristics when the sliding surface s reaches 0, as shown in Equations (17) and (18). To compare the convergence characteristics of the two methods before reaching the sliding surface through numerical simulation, the same parameters were used for the simulations.

In Figure 1, we present the convergence trajectories of the initial state of the AUV model towards the desired state for three distinct sliding mode control approaches.It is clear that NFTSM in this paper and NTSM demonstrate remarkable precision in achieving the desired poses. In contrast, the classical NFTSM showed a slightly lower overall convergence accuracy. When the NFTSM in this study and the NTSM are configured with identical parameters, the sliding mode control method proposed in this study shows faster convergence, with a convergence error of nearly zero.

Figure 1
Six plots comparing position and angle errors over time for three methods against a zero reference line.The six line plots are arranged in a 3 by 2 grid. Each subplot compares four curves labeled in the legend as “in this paper” (blue solid curve), “in open square bracket 18close square bracket ” (orange solid curve), “in open square bracket 29close square bracket ” (yellow solid curve), and “reference” (red dashed curve). All subplots share a Time (seconds) on the horizontal axis from 0 to 10 seconds, and each plot tracks a different error variable over time. In the top-left subplot, the vertical axis is labeled Position x error (meter), and the legend represents x subscript e. All curves begin with positive error values near 8 to 10 meters on the vertical axis, then oscillate downward toward zero, crossing into negative values at different times, reaching a minimum of negative 4. The red reference line remains constant near zero. The top-right subplot has the vertical axis labeled Position y error (meter). The three method curves again start from high positive values near 5 meters and decay with oscillations reaching a minimum of negative 2, settling close to zero. The red reference line is flat near zero. In the middle-left subplot, the vertical axis is Position z error (meter). All method curves start from approximately 6 meters, rapidly decay toward zero without oscillation, and converge at the bottom by around 5 seconds. The reference line is again constant near zero. The middle-right subplot shows the Angle roll error (degrees). The curves begin near 1 degree, decreasing quickly toward zero, with the blue curve showing more oscillatory behavior. A small inset zooms into the time interval from approximately 4.8 to 6.1 seconds, where small high-frequency oscillations are shown for all three methods, with values on the order of 10 to negative 3 power degrees. The bottom-left subplot plots the Angle pitch error (degrees). All curves start near 1 degree, rapidly decrease toward zero, and converge by about 3 seconds. A small inset shows oscillations between 5.4 and 6.1 seconds with values around 10 to negative 3 degrees. The bottom-right subplot displays the Angle yaw error (degrees). All curves start near zero, increasing over time. The blue curve rises quickly and levels near 1 degree, the orange curve reaches a slightly higher value, and the yellow curve increases less gradually. The red reference line remains at a constant high value across the entire time axis. Note: All numerical data values are approximated.

Six degrees of freedom trajectory tracking error

Figure 1
Six plots comparing position and angle errors over time for three methods against a zero reference line.The six line plots are arranged in a 3 by 2 grid. Each subplot compares four curves labeled in the legend as “in this paper” (blue solid curve), “in open square bracket 18close square bracket ” (orange solid curve), “in open square bracket 29close square bracket ” (yellow solid curve), and “reference” (red dashed curve). All subplots share a Time (seconds) on the horizontal axis from 0 to 10 seconds, and each plot tracks a different error variable over time. In the top-left subplot, the vertical axis is labeled Position x error (meter), and the legend represents x subscript e. All curves begin with positive error values near 8 to 10 meters on the vertical axis, then oscillate downward toward zero, crossing into negative values at different times, reaching a minimum of negative 4. The red reference line remains constant near zero. The top-right subplot has the vertical axis labeled Position y error (meter). The three method curves again start from high positive values near 5 meters and decay with oscillations reaching a minimum of negative 2, settling close to zero. The red reference line is flat near zero. In the middle-left subplot, the vertical axis is Position z error (meter). All method curves start from approximately 6 meters, rapidly decay toward zero without oscillation, and converge at the bottom by around 5 seconds. The reference line is again constant near zero. The middle-right subplot shows the Angle roll error (degrees). The curves begin near 1 degree, decreasing quickly toward zero, with the blue curve showing more oscillatory behavior. A small inset zooms into the time interval from approximately 4.8 to 6.1 seconds, where small high-frequency oscillations are shown for all three methods, with values on the order of 10 to negative 3 power degrees. The bottom-left subplot plots the Angle pitch error (degrees). All curves start near 1 degree, rapidly decrease toward zero, and converge by about 3 seconds. A small inset shows oscillations between 5.4 and 6.1 seconds with values around 10 to negative 3 degrees. The bottom-right subplot displays the Angle yaw error (degrees). All curves start near zero, increasing over time. The blue curve rises quickly and levels near 1 degree, the orange curve reaches a slightly higher value, and the yellow curve increases less gradually. The red reference line remains at a constant high value across the entire time axis. Note: All numerical data values are approximated.

Six degrees of freedom trajectory tracking error

Close modal

The proposed sliding mode control method demonstrates clear superiority across all six pose variables, except for the pitch angle. The convergence variances of the pitch angle, calculated using Eq.(18), are 1.8139e−05, 5.9779e−04 and 0.0172, respectively. This analysis confirms the advantages of the approach presented in this study.

(18)

This paper quantitatively validates the performance of our method by calculating the root mean square error of approximation (RMSEA) of the pitch tracking error from the experimental results. The proposed method achieves a minimum root mean square error (RMSE) value of 1.8139e−05, which is smaller than the other two methods (5.9779e−04 and 0.0172). The proposed method maintains extremely low tracking error in the presence of model uncertainty and external disturbances, demonstrating its robustness. Figures 1 and 2 show that our method converges faster. The lower RMSE value reflects the system’s ability to reach a steady state more quickly, confirming its rapid convergence.

Figure 2 illustrate the trajectory tracking curves for the three sliding mode control methods, and Figure 3 shows the necessary control input profiles. A comparative analysis between the NFTSM and classical NFTSM reveals that the proposed sliding mode control exhibits not only a faster convergence rate but also a lower control input requirement, thereby making it more suitable for practical use.

Figure 2
Six plots comparing x, y, z, roll, pitch, and yaw errors over time for three methods and a zero slash error reference.The figure is arranged as a 3 by 2 grid of line plots, each comparing four curves: “in this paper” (blue solid line), “in open square bracket 18close square bracket ” (orange solid line), “in open square bracket 29close square bracket ” (yellow solid line), and “reference” (red dashed line). All plots share a Time (seconds) on the horizontal axis from 0 to 10 seconds. The left column presents position errors (x, z, and pitch), and the right column presents position y error, angle roll error, and angle yaw error. In the top-left subplot, the vertical axis is labeled Position x (meters). The curves begin near positive values around 10 meters, decrease, and oscillate with two major peaks at 3 seconds and 9 seconds, and troughs at 0.5 seconds and 6.5 seconds before converging toward zero around 10 seconds. The top-right subplot shows Position y (meters). Each curve starts near 5 meters, then decays with oscillatory motions, crossing zero at different times. The middle-left subplot shows Position z (meters). The curves start near 5 meters and monotonically decay toward zero with no oscillations, converging by approximately 5 seconds. The reference curve begins at zero meters on the vertical axis and decreases to 10 seconds on the horizontal axis. The middle-right subplot is labeled Angle roll (degree). The curves begin near 1 degree, descend toward zero, and the red curve displays the largest oscillations. A small inset at roughly 5.8 to 7.4 seconds zooms in on the residual oscillations, which appear at amplitudes around plus or minus 0.01 degrees. In the bottom-left subplot, the vertical axis is Angle pitch (degrees). All curves start nearly at 1 degree, then rapidly decay toward zero. A small inset focuses on the narrow region from about 6.6 to 7.4 seconds, showing minor oscillations on the order of plus or minus 0.01 degrees. The bottom-right subplot shows Angle yaw (degrees). All curves start from zero and rise steadily. The blue and orange curves rise fastest toward values between 80 degrees and 100 degrees, the yellow curve ascends more slowly, and the red dashed reference line remains at a constant high value near 90°. Note: All numerical data values are approximated.

Six degrees of freedom state trajectory

Figure 2
Six plots comparing x, y, z, roll, pitch, and yaw errors over time for three methods and a zero slash error reference.The figure is arranged as a 3 by 2 grid of line plots, each comparing four curves: “in this paper” (blue solid line), “in open square bracket 18close square bracket ” (orange solid line), “in open square bracket 29close square bracket ” (yellow solid line), and “reference” (red dashed line). All plots share a Time (seconds) on the horizontal axis from 0 to 10 seconds. The left column presents position errors (x, z, and pitch), and the right column presents position y error, angle roll error, and angle yaw error. In the top-left subplot, the vertical axis is labeled Position x (meters). The curves begin near positive values around 10 meters, decrease, and oscillate with two major peaks at 3 seconds and 9 seconds, and troughs at 0.5 seconds and 6.5 seconds before converging toward zero around 10 seconds. The top-right subplot shows Position y (meters). Each curve starts near 5 meters, then decays with oscillatory motions, crossing zero at different times. The middle-left subplot shows Position z (meters). The curves start near 5 meters and monotonically decay toward zero with no oscillations, converging by approximately 5 seconds. The reference curve begins at zero meters on the vertical axis and decreases to 10 seconds on the horizontal axis. The middle-right subplot is labeled Angle roll (degree). The curves begin near 1 degree, descend toward zero, and the red curve displays the largest oscillations. A small inset at roughly 5.8 to 7.4 seconds zooms in on the residual oscillations, which appear at amplitudes around plus or minus 0.01 degrees. In the bottom-left subplot, the vertical axis is Angle pitch (degrees). All curves start nearly at 1 degree, then rapidly decay toward zero. A small inset focuses on the narrow region from about 6.6 to 7.4 seconds, showing minor oscillations on the order of plus or minus 0.01 degrees. The bottom-right subplot shows Angle yaw (degrees). All curves start from zero and rise steadily. The blue and orange curves rise fastest toward values between 80 degrees and 100 degrees, the yellow curve ascends more slowly, and the red dashed reference line remains at a constant high value near 90°. Note: All numerical data values are approximated.

Six degrees of freedom state trajectory

Close modal
Figure 3
Six plots showing x, y, z, roll, pitch, and yaw control inputs over time for three different control methods.The figure is arranged as a 3 by 2 grid of line plots, each showing control input versus Time (seconds) from 0 to 10 seconds. In every subplot, three curves are compared: “control input in this paper” (blue solid line), “control input in open square bracket 18close square bracket ” (orange solid line), and “control input in open square bracket 29close square bracket ” (yellow solid line). In the top-left subplot, the vertical axis is labeled Position x control input (Newtons), with values ranging between approximately negative 20 and 25 Newtons. All curves start with higher amplitude oscillations around 20 newton that gradually decrease, transitioning into smaller sustained oscillations toward the end of the time interval. The top-right subplot is labeled Position y control input (Newtons). Values range from roughly negative 25 to 50 newtons. The curves show an initial large positive input around 50 newton that quickly decays into oscillatory behavior. The oscillations gradually diminish, with different amplitudes across the three methods. In the middle-left subplot, the vertical axis is Position z control input (Newtons). Values span from negative 125 to about 25 newton. The curves start with a sharp negative input around negative 30 newton for two methods and a sharp positive input for one method. All then converge into low-amplitude oscillations after the first few seconds. The middle-right subplot displays Angle roll control input (newton meter). Values range from negative 5 to 25 newton meters. One curve exhibits a distinct peak near 5 to 25 newton meter at around 1 s, while the others rise more modestly. All curves eventually converge toward zero with slight oscillations. The bottom-left subplot shows Angle pitch control input (newton meter), ranging from negative 5 to 30 newton meters. One curve rises sharply to a peak of around 30 newton meters before descending, while the other curves peak at lower values. All eventually settle toward zero with small-amplitude oscillations. The bottom-right subplot presents Angle yaw control input (newton meters). Values extend from negative 25 to 10 newton meters. One curve features an early negative dip before rising toward zero, while the other two curves rise directly from zero. All three converge into nearly constant, small values between 0 and 5 newton meters. Note: All numerical data values are approximated.

Required control input for six degrees of freedom

Figure 3
Six plots showing x, y, z, roll, pitch, and yaw control inputs over time for three different control methods.The figure is arranged as a 3 by 2 grid of line plots, each showing control input versus Time (seconds) from 0 to 10 seconds. In every subplot, three curves are compared: “control input in this paper” (blue solid line), “control input in open square bracket 18close square bracket ” (orange solid line), and “control input in open square bracket 29close square bracket ” (yellow solid line). In the top-left subplot, the vertical axis is labeled Position x control input (Newtons), with values ranging between approximately negative 20 and 25 Newtons. All curves start with higher amplitude oscillations around 20 newton that gradually decrease, transitioning into smaller sustained oscillations toward the end of the time interval. The top-right subplot is labeled Position y control input (Newtons). Values range from roughly negative 25 to 50 newtons. The curves show an initial large positive input around 50 newton that quickly decays into oscillatory behavior. The oscillations gradually diminish, with different amplitudes across the three methods. In the middle-left subplot, the vertical axis is Position z control input (Newtons). Values span from negative 125 to about 25 newton. The curves start with a sharp negative input around negative 30 newton for two methods and a sharp positive input for one method. All then converge into low-amplitude oscillations after the first few seconds. The middle-right subplot displays Angle roll control input (newton meter). Values range from negative 5 to 25 newton meters. One curve exhibits a distinct peak near 5 to 25 newton meter at around 1 s, while the others rise more modestly. All curves eventually converge toward zero with slight oscillations. The bottom-left subplot shows Angle pitch control input (newton meter), ranging from negative 5 to 30 newton meters. One curve rises sharply to a peak of around 30 newton meters before descending, while the other curves peak at lower values. All eventually settle toward zero with small-amplitude oscillations. The bottom-right subplot presents Angle yaw control input (newton meters). Values extend from negative 25 to 10 newton meters. One curve features an early negative dip before rising toward zero, while the other two curves rise directly from zero. All three converge into nearly constant, small values between 0 and 5 newton meters. Note: All numerical data values are approximated.

Required control input for six degrees of freedom

Close modal

In Figure 4, when the initial state is not on the target trajectory, the control method proposed in this paper can converge to the target trajectory in a shorter time, while the other two methods converge more slowly, which has the advantage of convergence speed. In the convergence process, the control method proposed in this study has a higher convergence accuracy. The two control methods, classical NFTSM and NTSM, have obvious tracking errors in the figure, which intuitively shows that the method proposed in this study has the advantage of convergence accuracy.

Figure 4
A 3 D plot comparing four trajectories from an initial high-z point to convergent final positions.The 3 D plot showing four trajectories in an x, y, and z space. The horizontal axes, on the right, represent x (meters) and ranges from 0 to 10, and on the left represent y (meters) and ranges from negative 2 to 4. The vertical axis represents z (meters) and ranges from negative 2 to 6. A legend identifies the plotted curves: a red curve labeled “trajectory in this paper”, a green curve labeled “trajectory in open square bracket 18 close square bracket ”, a blue curve labeled “trajectory in open square bracket 29 close square bracket ”, and a purple dashed curve labeled “reference”. All trajectories begin from an annotated “Initial position” located at a high positive z-value and near the origin in x and y. The curves descend in spiraling or looping shapes as they move outward along both horizontal axes. Near the lower z-region (approximately negative 2 meter), several “Final position” labels appear in red, green, blue, and purple, marking where each method’s trajectory converges. Note: All numerical data values are approximated.

Underwater robot trajectory

Figure 4
A 3 D plot comparing four trajectories from an initial high-z point to convergent final positions.The 3 D plot showing four trajectories in an x, y, and z space. The horizontal axes, on the right, represent x (meters) and ranges from 0 to 10, and on the left represent y (meters) and ranges from negative 2 to 4. The vertical axis represents z (meters) and ranges from negative 2 to 6. A legend identifies the plotted curves: a red curve labeled “trajectory in this paper”, a green curve labeled “trajectory in open square bracket 18 close square bracket ”, a blue curve labeled “trajectory in open square bracket 29 close square bracket ”, and a purple dashed curve labeled “reference”. All trajectories begin from an annotated “Initial position” located at a high positive z-value and near the origin in x and y. The curves descend in spiraling or looping shapes as they move outward along both horizontal axes. Near the lower z-region (approximately negative 2 meter), several “Final position” labels appear in red, green, blue, and purple, marking where each method’s trajectory converges. Note: All numerical data values are approximated.

Underwater robot trajectory

Close modal

This paper introduces an innovative approach to NFTSM control, which stands out for its robust performance in dynamic environments. The proposed control scheme is characterized by its nonsingular nature, ensuring that the control inputs remain well-defined even in the presence of uncertainties. A comprehensive proof of the finite-time convergence property was established through the application of a Lyapunov function, demonstrating the effectiveness of this method in achieving stability. One of the significant advantages of this approach is its ability to effectively mitigate the adverse effects of imprecise internal parameters and external disturbances, which often challenge control systems. By addressing these issues, this method enhances the reliability and accuracy of system performance.

The practical implementation of this control strategy is demonstrated through experiments conducted on an underwater robot system, specifically designed to achieve precise pose control. The experimental results revealed that the robot can rapidly converge all six degrees of freedom to nearly zero, indicating high levels of accuracy and responsiveness.

Furthermore, a comparative analysis with other representative NFTSM control methods highlights the superior convergence accuracy achieved by this approach. The findings suggest that this method not only improves performance in terms of speed but also ensures a higher degree of precision in control, making it a significant advancement in the field of control systems for underwater robotics. Due to its inherent robustness, nonsingularity, fast convergence and weak chattering characteristics, it will be applied to actual underwater robot systems in the future.

Bessa
,
W.M.
,
Dutra
,
M.S.
and
Kreuzer
,
E.
(
2010
), “
An adaptive fuzzy sliding mode controller for remotely operated underwater vehicles
”,
Robotics and Autonomous Systems
, Vol. 
58
No. 
1
, pp. 
16
-
26
, ISSN:
[PubMed]
, doi: .
Chen
,
Q.
,
Wang
,
M.
,
Nan
,
Y.
and
Wei
,
C.
(
2023
), “
Unknown system dynamic estimator-based two-phase power reaching law control for DC–DC buck converter
”,
IEEE Transactions on Circuits and Systems I: Regular Papers
, Vol. 
70
No. 
11
, pp. 
4641
-
4649
, doi: .
Elmokadem
,
T.
,
Mohamed
,
Z.
and
Youcef-Toumi
,
K.
(
2017
), “
Terminal sliding mode control for the trajectory tracking of underactuated Autonomous Underwater Vehicles
”,
Ocean Engineering
, Vol. 
129
, pp. 
613
-
625
, ISSN:
[PubMed]
, doi: .
Feng
,
Y.
,
Yu
,
X.
and
Man
,
Z.
(
2002
), “
Non-singular terminal sliding mode control of rigid manipulators
”,
Automatica
, Vol. 
38
No. 
12
, pp. 
2159
-
2167
, ISSN:
[PubMed]
, doi: .
Fischer
,
N.
,
Hughes
,
D.
,
Walters
,
P.
,
Schwartz
,
E.M.
and
Dixon
,
W.E.
(
2014
), “
Nonlinear RISE-based control of an autonomous underwater vehicle
”,
IEEE Transactions on Robotics
, Vol. 
30
No. 
4
, pp. 
845
-
852
, doi: .
Geranmehr
,
B.
and
Nekoo
,
S.R.
(
2015
), “
Nonlinear suboptimal control of fully coupled non-affine six-DOF autonomous underwater vehicle using the state-dependent Riccati equation
”,
Ocean Engineering
, Vol. 
96
, pp. 
248
-
257
, ISSN:
[PubMed]
, doi: .
González-García
,
J.
,
Gómez-Espinosa
,
A.
,
Cuan-Urquizo
,
E.
,
García-Valdovinos
,
L.G.
,
Salgado-Jiménez
,
T.
and
Cabello
,
J.A.E.
(
2020
), “
Autonomous underwater vehicles: localization, navigation, and communication for collaborative missions
”,
Applied Sciences
, Vol. 
10
No. 
4
, p. 
1256
, ISSN:
[PubMed]
, doi: .
Guerrero
,
J.
,
Chemori
,
A.
,
Torres
,
J.
and
Creuze
,
V.
(
2023
), “
Time-delay high-order sliding mode control for trajectory tracking of autonomous underwater vehicles under disturbances
”,
Ocean Engineering
, Vol. 
268
, 113375, ISSN:
[PubMed]
, doi: .
Guerrero
,
J.
,
Chemori
,
A.
,
Creuze
,
V.
and
Torres
,
J.
(
2024
), “
Improved adaptive high-order sliding mode-based control for trajectory tracking of autonomous underwater vehicles
”,
IEEE Journal of Oceanic Engineering
, Vol. 
49
No. 
4
, pp. 
1337
-
1349
, doi: .
Hou
,
H.
,
Yu
,
X.
,
Xu
,
L.
,
Rsetam
,
K.
and
Cao
,
Z.
(
2020
), “
Finite-time continuous terminal sliding mode control of servo motor systems
”,
IEEE Transactions on Industrial Electronics
, Vol. 
67
No. 
7
, pp. 
5647
-
5656
, doi: .
Hu
,
Q.
and
Jiang
,
B.
(
2018
), “
Continuous finite-time attitude control for rigid spacecraft based on angular velocity observer
”,
IEEE Transactions on Aerospace and Electronic Systems
, Vol. 
54
No. 
3
, pp. 
1082
-
1092
, doi: .
Karkoub
,
M.
,
Wu
,
H.-M.
and
Hwang
,
C.-L.
(
2017
), “
Nonlinear trajectory-tracking control of an autonomous underwater vehicle
”,
Ocean Engineering
, Vol. 
145
, pp. 
188
-
198
, ISSN:
[PubMed]
, doi: .
Li
,
H.
,
Wang
,
J.
,
Du
,
H.
and
Karimi
,
H.R.
(
2018
), “
Adaptive sliding mode control for Takagi–Sugeno fuzzy systems and its applications
”,
IEEE Transactions on Fuzzy Systems
, Vol. 
26
No. 
2
, pp. 
531
-
542
, doi: .
Li
,
J.
,
Xiang
,
X.
,
Dong
,
D.
and
Yang
,
S.
(
2023
), “
Prescribed time observer based trajectory tracking control of autonomous underwater vehicle with tracking error constraints
”,
Ocean Engineering
, Vol. 
274
, 114018, ISSN:
[PubMed]
, doi: .
Li
,
Z.
,
Tian
,
B.
,
Zhang
,
X.
and
Zong
,
Q.
(
2024
), “
Design and implementation of continuous bounded finite-time attitude control of reusable launch vehicle
”,
IEEE Transactions on Industrial Electronics
, Vol. 
71
No. 
2
, pp. 
1917
-
1926
, doi: .
Liu
,
X.
,
Zhang
,
M.
,
Chu
,
Z.
and
Rogers
,
E.
(
2023
), “
A sphere region tracking control scheme for underwater vehicles
”,
IEEE Transactions on Vehicular Technology
, Vol. 
72
No. 
8
, pp. 
9835
-
9844
, doi: .
Mishra
,
J.P.
,
Yu
,
X.
,
Jalili
,
M.
and
Feng
,
Y.
(
2016
), “
On fast terminal sliding-mode control design for higher order systems
”,
IECON 2016-42nd Annual Conference of the IEEE Industrial Electronics Society
, pp. 
252
-
257
, doi: .
Mu
,
C.
and
He
,
H.
(
2018
), “
Dynamic behavior of terminal sliding mode control
”,
IEEE Transactions on Industrial Electronics
, Vol. 
65
No. 
4
, pp. 
3480
-
3490
, doi: .
Plestan
,
F.
,
Shtessel
,
Y.
,
Brégeault
,
V.
and
Poznyak
,
A.
(
2010
), “
New methodologies for adaptive sliding mode control
”,
International Journal of Control
, Vol. 
83
No. 
9
, pp. 
1907
-
1919
, doi: .
Qiao
,
L.
,
Yi
,
B.
,
Wu
,
D.
and
Zhang
,
W.
(
2017
), “
Design of three exponentially convergent robust controllers for the trajectory tracking of autonomous underwater vehicles
”,
Ocean Engineering
, Vol. 
134
, pp. 
157
-
172
, ISSN:
[PubMed]
, doi: .
Qin
,
H.
,
Si
,
J.
,
Wang
,
N.
and
Gao
,
L.
(
2023
), “
Fast fixed-time nonsingular terminal sliding-mode formation control for autonomous underwater vehicles based on a disturbance observer
”,
Ocean Engineering
, Vol. 
270
, 113423, ISSN:
[PubMed]
, doi: .
Shen
,
Y.
and
Huang
,
Y.
(
2009
), “
Uniformly observable and globally lipschitzian nonlinear systems admit global finite-time observers
”,
IEEE Transactions on Automatic Control
, Vol. 
54
No. 
11
, pp. 
2621
-
2625
, doi: .
Wang
,
H.-Du
,
Zhai
,
Y.X.
,
Shah
,
U.H.
,
Karkoub
,
M.
and
Li
,
M.
(
2023a
), “
Adaptive fuzzy control of underwater vehicle manipulator system with dead-zone band input nonlinearities via fuzzy performance and disturbance observers
”,
Ocean Engineering
, Vol. 
277
, 114194, ISSN:
[PubMed]
, doi: .
Wang
,
X.
,
Xu
,
J.
and
Liu
,
P.
(
2023b
), “
Adaptive non-singular integral terminal sliding mode-based fault tolerant control for autonomous underwater vehicles
”,
Ocean Engineering
, Vol. 
267
, 113299, doi: .
Wen
,
L.
,
Yu
,
S.
,
Zhao
,
Y.
and
Yan
,
Y.
(
2023
), “
Adaptive dynamic event-triggered consensus control of multiple autonomous underwater vehicles
”,
International Journal of Control
, Vol. 
96
No. 
3
, pp. 
746
-
756
, doi: .
Xu
,
S.S.-D.
,
Chen
,
C.-C.
and
Wu
,
Z.-L.
(
2015
), “
Study of nonsingular fast terminal sliding-mode fault-tolerant control
”,
IEEE Transactions on Industrial Electronics
, Vol. 
62
No. 
6
, pp. 
3906
-
3913
, doi: .
Yan
,
Z.
,
Yan
,
J.
,
Cai
,
S.
,
Yu
,
Y.
and
Wu
,
Y.
(
2023
), “
Robust MPC-based trajectory tracking of autonomous underwater vehicles with model uncertainty
”,
Ocean Engineering
, Vol. 
286
, 115617, ISSN:
[PubMed]
, doi: .
Yang
,
L.
and
Yang
,
J.
(
2011
), “
Nonsingular fast terminal sliding-mode control for nonlinear dynamical systems
”,
International Journal of Robust and Nonlinear Control
, Vol. 
21
No. 
16
, pp. 
1865
-
1879
, doi: .
Yu
,
X.
and
Zhihong
,
M.
(
2002
), “
Fast terminal sliding-mode control design for nonlinear dynamical systems
”,
IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications
, Vol. 
49
No. 
2
, pp. 
261
-
264
, doi: .
Zhu
,
Z.
,
Duan
,
Z.
,
Qin
,
H.
and
Xue
,
Y.
(
2023
), “
Adaptive neural network fixed-time sliding mode control for trajectory tracking of underwater vehicle
”,
Ocean Engineering
, Vol. 
287
, 115864, ISSN:
[PubMed]
, doi: .
Published in Journal of Intelligent Manufacturing and Special Equipment. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at Link to the terms of the CC BY 4.0 licence.

or Create an Account

Close Modal
Close Modal