Path following control of an underactuated unmanned marine vehicle with model asymmetry in the presence of ocean current disturbances

Abstract

Path following control system of an underactuated unmanned marine vehicle is an important guarantee for its successful operation, the path following control problem of an underactuated and asymmetrical unmanned marine vehicle sailing in the presence of ocean current disturbances is addressed. Asymmetry of the unmanned marine vehicle model and higher order velocity coupling terms of damping coefficients of the unmanned marine vehicle are discussed, which would improve the accuracy of the path following control system. A kind of global differential homeomorphism transformation design is proposed to solve the difficulty of nonzero term appearing in the non-diagonal elements of the system inertia coefficient matrix and damping coefficient matrix, which is caused by the model asymmetry of unmanned marine vehicle. An improved line-of-sight guidance algorithm is presented by introducing longitudinal position error and tracking error weight factor into traditional line-of-sight algorithm, which could speed up the path following process, meanwhile the method could be extended to the application of curve path following. Virtual velocity in the tangent direction of the path to be followed is designed for the control system, by increasing a virtual control input, the underactuated control system is simplified, and the higher order velocity coupling terms of damping coefficients are integrated considered in the virtual control law. Stability of the path following control algorithm proposed for unmanned marine vehicle is proved by Lyapunov theory, and some simulation experiments are carried out to verify the effectiveness of the path following control system designed.

Keywords

Unmanned marine vehicle path following line of sight underactuated control system ocean current disturbances

Introduction

Unmanned marine vehicle (UMV), usually contained unmanned/autonomous underwater vehicle (UUV/AUV) and unmanned/autonomous surface vehicle (USV/ASV), can be used to perform a multitude of different tasks, such as mineral resources sampling, offshore oil and gas operations, ocean engineering maintenance, and military reconnaissance, and it is attracting more and more interest from the scientific, commercial, and naval sectors.^1

–5 Although much advancements have been realized in this area, the demand for more advantages in navigation, guidance, and control system for UMVs continues to grow, as more and more vehicle autonomy is required.^6

–10 In practical implementation, many UMVs are designed of underactuated configurations due to practical considerations, such as reducing weight and/or cost and improving structural strength.^11
–13 Path following control system is one of the most fundamental and critical control systems during the whole operation process of UMV; after the scheduled task is determined, the UMV needs to accurately sail on the specified path according to the requirements of mission planning, and inaccurate path following may lead to task failure.^14

–18 In recent years, more and more studies have been conducted on ocean exploration by UMV, and there are many research studies on the path following control of an underactuated UMV for ocean exploration, as path following control system is one the most important control systems for an UMV.^19

–22

Accurate path following control plays an important role for UMV in the oceanic surveys and exploration, an adaptive nonlinear second-order sliding mode controller is proposed to eliminate the chattering motion through a sliding surface during the path following control by Zhang et al.,²³ and David et al. proposed a novel improved proportional–integral–derivative control approach that enables variable ballast systems on an UMV to complement the underwater vehicle’s hydroplanes when altitude keeping over a variable seabed,²⁴ while subject to complex unknown nonlinearities including un-modeled dynamics, uncertainties, and unknown disturbances, a novel fuzzy uncertainty observer-based path following control scheme for an underactuated marine vehicle is proposed in the study by Wang et al.²⁵ The cooperative path following problem of multiple underactuated UMVs involves two tasks,²⁶ which may be used to cooperative topographic survey, and the first one is to force each UMV to converge to the desired parameterized path, while the second one is to satisfy the requirement of a cooperative behavior along the paths. A path following control scheme for deep-sea work class remotely operated vehicles subjected to disturbing forces and moments generated by parameter variations, ocean currents, umbilical cables, and other unknown disturbances is designed,²⁷ while based on a multi-body system concept, a model is designed based on an adaptive controller and a disturbance model is analyzed, which compensates the dynamics of a finite element tether and multi-body underwater vehicle manipulator.²⁸ Peng and Wang present a design method for output feedback path following control of underactuated UMV moving in a vertical plane without using surge, heave, and pitch velocities, and an extended state observer is developed to recover the unmeasured velocities as well as to estimate total uncertainty induced by internal model uncertainty and external disturbance²⁹, while Jiang et al. discuss the task of both optimizing and following routes that present a ladder shape; a planning method is proposed to determine a nearly optimal coverage route is designed, and then the track spacing is optimized considering the seabed type and the sonar range for the purpose of increasing detection probability.³⁰ The problem of three-dimensional (3D) path following control for an underactuated UMV in the presence of parameter uncertainties and external disturbances is addressed,³¹ while a teach and repeat path following method is proposed by King et al.³² for an UMV navigating long distances in environments where external navigation aides are denied, and the method utilizes sonar images to construct a series of reference views along a path, stored as a topological map. Aiming to fuse bathymetric data from both systems for a full coverage survey with a confident map of these data, the selected model for the fusion process is derived from belief functions.³³ Extension to fuzzy subsets allows us to integrate the precision and the reliability of the source data, and the targeted estimator of the bathymetry error associated with every sounding involved in the fusion process should give an objective quality of the fused bathymetry.

By analyzing the existing research results described above, we can find that, at present, most of the research results on path following of an UMV for ocean exploration are based on fully symmetric UMV model, meanwhile the damping coefficient higher order velocity coupling terms of UMV model are not considered.^34

–38 While considering the higher accuracy requirements of the UMV’s path following control system for the ocean exploration, a novel path following control system, which considers the asymmetry and the damping coefficient higher order velocity coupling terms of the UMV model, is designed based on an improved line-of-sight (LOS) guidance algorithm for and underactuated UMV used for ocean exploration in the presence of ocean current disturbances. Compared with the existing research results on path following control algorithm of an underactuated UMV for ocean exploration, the main contributions of this article are as follows: (i) Asymmetry of the UMV model is considered, and a kind of global differential homeomorphism transformation design is proposed to solve the difficulty of nonzero term appearing in the non-diagonal elements of the system inertia coefficient matrix and damping coefficient matrix, which is caused by the model asymmetry of UMV. (ii) Higher order velocity coupling terms of damping coefficients of the UMV are discussed, which include both linear damping and nonlinear damping, and an ingenious intermediate process design is presented for dealing with the complex nonlinear terms. (iii) An improved LOS guidance algorithm for path following control of UMV is proposed, which combines with the longitudinal tracking error and lateral tracking error in the designed guidance algorithm. (iv) The virtual velocity in the tangent direction of the path to be followed is designed for the control system, which increases a virtual control input and leads to simplify the underactuated control system. (v) Path following control algorithm designed in this article could achieve both straight line path and curve path following, while some algorithm could only achieve straight line path or curve path following. (vi) Ocean current disturbance is considered in this article, and the ocean current relative coordinate system and ocean current relative velocity are designed in this article to integrate the ocean current disturbance into the kinematics equations and dynamics equations of UMV.

The remainder of the article is organized as follows. UMV modeling and the path following problem of the underactuated UMV is formulated in the “Problem formulation” section, meanwhile a kind of global differential homeomorphism transformation design is proposed. Path following control laws are presented in the “Path following controller design for an underactuated UMV” section, an improved LOS guidance algorithm is proposed and the virtual velocity in the tangent direction of the path to be followed is designed. Stability of the control laws designed for the underactuated UMV is proved by Lyapunov stability theory in the “Stability analysis of path following control algorithm” section. In the “Simulation experiments” section, some simulation experiments of path following control during ocean exploration are carried out to verify the effectiveness of the control algorithm proposed and the path control system design for UMV. In the “Discussion” section, some discussions about our research results on path following control of UMV are made. At last, in the “Conclusions” section, some main conclusions of path following control of an underactuated UMV for ocean exploration in the presence of ocean current disturbances are summarized.

Problem formulation

In this section, UMV modeling and the path following problem of the UMV are formulated, and a kind of global differential homeomorphism transformation design is proposed.

UMV modeling

As shown in Figure 1, Cartesian coordinate system ${I} - x o y$ is a North-East-Down inertial coordinate system, and the path $Θ (s)$ is a preplanned path which UMV needs to follow.

Figure 1.

Reference coordinate systems of UMV path following, including Serret–Frenet coordinate system, inertial coordinate system and UMV body-fixed coordinate system. UMV: unmanned marine vehicle.

The point $o_{s f}$ is a path following reference point on path $Θ (s)$ , establishing Serret–Frenet coordinate system ${S F} - x_{s f} o_{s f} y_{s f}$ based on point $o_{s f}$ , where the $x_{s f}$ -axis direction is the tangent direction of point $o_{s f}$ on the path $Θ (s)$ , while the $y_{s f}$ -axis direction is the normal direction. The angle between the $x_{s f}$ -axis and the x-axis is $φ_{s f}$ .

Establishing UMV body-fixed coordinate system ${b} - x_{b} o_{b} y_{b}$ based on the UMV’s gravity center o_b, u is the longitudinal velocity of UMV, while v is the lateral velocity, and $U = \sqrt{u^{2} + v^{2}}$ , β is the drift angle. The position of o_b in the coordinate system ${S F} - x_{s f} o_{s f} y_{s f}$ is $(x_{e}, y_{e})$ , which is also the path following error.

According to Fossen,³⁹ UMV’s motion control system, usually could be divided into the horizontal motion control and vertical motion control, and those two subsystems are similar in the way. So horizontal plane control is usually considered first, then the vertical motion control could be the similar way, and then the horizontal plane path following control is discussed in this article, kinematics and dynamics equations of the underactuated UMV in the presence of constant ocean currents can be described as follows

{\begin{matrix} \dot{η} = J (η) υ \\ M {\dot{υ}}_{r} + C (υ_{r}) υ_{r} + D υ_{r} = τ \end{matrix}

with

η = {[x, y, φ]}^{Τ}, υ = {[u, v, r]}^{Τ}, υ_{r} = {[u_{r}, v_{r}, r]}^{Τ} = [u - u_{c}, v - v_{c}, r], τ = {[τ_{u}, 0, τ_{r}]}^{Τ},

J (η) = [\begin{matrix} cos φ & - sin φ & 0 \\ sin φ & cos φ & 0 \\ 0 & 0 & 1 \end{matrix}], C (υ_{r}) = [\begin{matrix} 0 & 0 & c_{13} \\ 0 & 0 & c_{23} \\ c_{31} & c_{32} & 0 \end{matrix}] = [\begin{matrix} 0 & 0 & - m_{22} v_{r} - m_{23} r \\ 0 & 0 & m_{11} u_{r} \\ m_{22} v_{r} + m_{23} r & - m_{11} u_{r} & 0 \end{matrix}]

M = [\begin{matrix} m_{11} & 0 & 0 \\ 0 & m_{22} & m_{23} \\ 0 & m_{32} & m_{33} \end{matrix}], D (υ_{r}) = D_{linear} (υ_{r}) + D_{nonlinear} (υ_{r}) = [\begin{matrix} d_{11}^{l} & 0 & 0 \\ 0 & d_{22}^{l} & d_{23}^{l} \\ 0 & d_{32}^{l} & d_{33}^{l} \end{matrix}] + [\begin{matrix} d_{11}^{n l} & 0 & 0 \\ 0 & d_{22}^{n l} & d_{23}^{n l} \\ 0 & d_{32}^{n l} & d_{33}^{n l} \end{matrix}]

where u_c and v_c are the ocean current velocities, $τ_{u}$ and $τ_{r}$ are the surge force and yaw moment of UMV, m₁₁, m₂₂, m₃₃, m₂₃, and m₃₂ are UMV’s inertia coefficients including added mass effects, and $d_{11}^{l}$ , $d_{22}^{l}$ , $d_{33}^{l}$ , $d_{23}^{l}$ , $d_{32}^{l}$ are all hydrodynamic linear damping coefficients, while $d_{11}^{n l}$ , $d_{22}^{n l}$ , $d_{33}^{n l}$ , $d_{23}^{n l}$ , $d_{32}^{n l}$ are all hydrodynamic nonlinear damping coefficients.

Expanding system (1) above leads to the UMV’s kinematics equations (2) and dynamics equations (3)

{\begin{matrix} \dot{x} = u cos φ - v sin φ \\ \dot{y} = u sin φ + v cos φ \\ \dot{φ} = r \end{matrix}

{\begin{cases} \dot{u} = \frac{m_{22}}{m_{11}} v_{r} r - \frac{d_{11}^{l} + d_{11}^{n l}}{m_{11}} u_{r} + \frac{m_{23}}{m_{11}} r^{2} + \frac{1}{m_{11}} τ_{u} \\ \dot{v} = - \frac{m_{11} m_{33} - m_{23}^{2}}{m_{22} m_{33} - m_{23}^{2}} u_{r} r - \frac{(m_{11} - m_{22}) m_{23}}{m_{22} m_{33} - m_{23}^{2}} u_{r} v_{r} - \frac{m_{33} (d_{22}^{l} + d_{22}^{n l}) - m_{23} (d_{32}^{l} + d_{32}^{n l})}{m_{22} m_{33} - m_{23}^{2}} v_{r} \\ - \frac{m_{33} (d_{23}^{l} + d_{23}^{n l}) - m_{23} (d_{33}^{l} + d_{33}^{n l})}{m_{22} m_{33} - m_{23}^{2}} r - \frac{m_{23}}{m_{22} m_{33} - m_{23}^{2}} τ_{r} \\ \dot{r} = \frac{(m_{11} - m_{22}) m_{22}}{m_{22} m_{33} - m_{23}^{2}} u_{r} v_{r} + \frac{(m_{11} - m_{22}) m_{23}}{m_{22} m_{33} - m_{23}^{2}} u_{r} r + \frac{m_{23} (d_{22}^{l} + d_{22}^{n l}) - m_{22} (d_{32}^{l} + d_{32}^{n l})}{m_{22} m_{33} - m_{23}^{2}} v_{r} \\ + \frac{m_{23} (d_{23}^{l} + d_{23}^{n l}) - m_{22} (d_{33}^{l} + d_{33}^{n l})}{m_{22} m_{33} - m_{23}^{2}} r + \frac{m_{22}}{m_{22} m_{33} - m_{23}^{2}} τ_{r} \end{cases}

Global differential homeomorphism transformation

To solve the difficulties that, nonzero term appears in the non-diagonal elements of the system inertia coefficient matrix and damping coefficient matrix, which is caused by the model asymmetry of UMV, a kind of global differential homeomorphism transformation is presented here.

Designing the following intermediate state variables

{\begin{cases} z_{1} = x + m_{23} [cos (φ) - 1] / m_{22} \\ z_{2} = y + m_{23} sin (φ) / m_{22} \\ z_{3} = φ \\ z_{4} = u \\ z_{5} = v + m_{23} r / m_{22} \\ z_{6} = r \\ z_{4 r} = u_{r} \\ z_{5 r} = v_{r} + m_{23} r / m_{22} \\ z_{4 c} = u_{c} \\ z_{5 c} = v_{c} \end{cases}

Then after state transformation design, system (2) and system (3) could be equivalently transformed into the following system (5)

{\begin{cases} {\dot{z}}_{1} = z_{4 r} cos (z_{3}) - z_{5 r} sin (z_{3}) + z_{4 c} cos (z_{3}) - z_{5 c} sin (z_{3}) \\ {\dot{z}}_{2} = z_{4 r} sin (z_{3}) + z_{5 r} cos (z_{3}) + z_{4 c} sin (z_{3}) + z_{5 c} cos (z_{3}) \\ {\dot{z}}_{3} = z_{6} \\ {\dot{z}}_{4 r} = m_{22} z_{5 r} z_{6} / m_{11} - d_{11} z_{4 r} / m_{11} + τ_{u} / m_{11} \\ {\dot{z}}_{5 r} = - m_{11} z_{4 r} z_{6} / m_{22} - d_{22} z_{5 r} / m_{22} + (d_{22} m_{23} - d_{23} m_{22}) z_{6} / m_{22}^{2} \\ {\dot{z}}_{6} = (m_{11} - m_{22}) m_{22} z_{4 r} z_{5 r} / (m_{22} m_{33} - m_{23}^{2}) + (d_{22} m_{23} - d_{32} m_{22}) z_{5 r} / (m_{22} m_{33} \\ - m_{23}^{2}) - (d_{22} m_{23}^{2} - d_{23} m_{22} m_{23} \\ - d_{32} m_{22} m_{23} + d_{33} m_{22}^{2}) z_{6} / [m_{22} (m_{22} m_{33} \\ - m_{23}^{2})] + m_{22} τ_{r} / (m_{22} m_{33} - m_{23}^{2}) \end{cases}

UMV path following error system

From Figure 1, UMV’s path following errors could be described as follows

[\begin{array}{l} x_{e} \\ y_{e} \\ φ_{e} \end{array}] = [\begin{matrix} cos (φ_{s f}) & sin (φ_{s f}) & 0 \\ - sin (φ_{s f}) & cos (φ_{s f}) & 0 \\ 0 & 0 & 1 \end{matrix}] [\begin{array}{l} x - x_{s f} \\ y - y_{s f} \\ φ_{f f} - φ_{s f} \end{array}]

where $φ_{f f} = φ + β$ and $β = arctan (z_{5} / z_{4})$ .

Designing the following intermediate state variables

{\begin{cases} z_{1 s f} = x_{s f} + m_{23} [cos (φ_{s f}) - 1] / m_{22} \\ z_{2 s f} = y_{s f} + m_{23} sin (φ_{s f}) / m_{22} \\ z_{3 s f} = φ_{s f} \end{cases}

And then after mathematical calculation and analysis, the path following errors described in equation (6) could be equivalently transformed into the following expression (8)

[\begin{array}{l} x_{e} \\ y_{e} \\ φ_{e} \end{array}] = [\begin{matrix} cos (z_{3 s f}) & sin (z_{3 s f}) & 0 \\ - sin (z_{3 s f}) & cos (z_{3 s f}) & 0 \\ 0 & 0 & 1 \end{matrix}] [\begin{array}{l} z_{1} - z_{1 s f} \\ z_{2} - z_{2 s f} \\ z_{3 f f} - z_{3 s f} \end{array}]

where $z_{3 f f} = φ_{f f}$ .

Differentiating both sides of expression (8) leads to

{\begin{cases} {\dot{x}}_{e} = K (s) \dot{s} y_{e} - \dot{s} + z_{4 r} cos (z_{3} - z_{3 s f}) - z_{5 r} sin (z_{3} - z_{3 s f}) + z_{4 c} cos (z_{3} - z_{3 s f}) - z_{5 c} sin (z_{3} - z_{3 s f}) \\ {\dot{y}}_{e} = - K (s) \dot{s} x_{e} + z_{4 r} sin (z_{3} - z_{3 s f}) + z_{5 r} cos (z_{3} - z_{3 s f}) + z_{4 c} sin (z_{3} - z_{3 s f}) + z_{5 c} cos (z_{3} - z_{3 s f}) \\ {\dot{φ}}_{e} = z_{6} + \dot{β} - K (s) \dot{s} \end{cases}

where s is the length of the path and $K (s)$ means curvature at the reference point.

At last, the problem of path following control of an underactuated UMV in a preplanned route is transformed into the problem of designing control input of $τ_{u}$ and $τ_{r}$ , which could have the path following errors $p_{e} = {[x_{e}, y_{e}, φ_{e}]}^{Τ}$ converge to 0, and make the longitudinal velocity u converge to desired velocity u_d or have the longitudinal velocity error $u_{e} = u - u_{d}$ converge to 0.

Path following controller design for an underactuated UMV

Path following control algorithms for an underactuated UMV in a preplanned route with the disturbance of ocean current is proposed in this section, and the main process of the control algorithm design and its implementation is shown in Figure 2, where $V_{c}$ is the velocity of constant ocean current disturbances and $β_{c}$ is the current direction angle, and $u_{c} = V_{c} cos β_{c}$ , $v_{c} = V_{c} sin β_{c}$ .

Figure 2.

Flow diagram of the underactuated UMV path following control algorithm design process and its implementation. UMV: unmanned marine vehicle.

Flow chart of our laboratory semi-physical simulation platform could be described as shown in Figure 3. It shows the hardware architecture of semi-physical simulation platform of UMV. Inside the dashed box is virtual simulation, while outside is physical simulation. The virtual simulation includes an Hewlett-Packard Company (HP) graphics workstation, an information center. It is used for scene simulation, hydrodynamic calculations, and the sensor simulation of sound and light vision. The information center is equipped with an 8-channel Peripheral Component Interconnect (PCI) serial port card, module Analog Data Acquisition Module (ADAM), and it is used to simulate the motion information sensors, such as Global Positioning System (GPS), Inertial Measurement Unit (IMU), inertial navigator, and data interface of propellers. It is of Windows operation system and it communicates with HP graphics workstation through Transmission Control Protocol/Internet Protocol (TCP/IP) protocol.

Figure 3.

Flow diagram of the semi-physical simulation platform.

The monitor terminal is used to receive condition information through network and serial port protocol to monitor the working condition and record it for further data analysis. Besides the monitor, outside dashed box, they are all physical simulation sections. The physical simulation is composed of vision processing unit, planning unit, and path following control unit. The planning unit and path following control unit adopt multi-plate embedded system based on PC104, which comprises a core module integrated with a Pentium M processor, a power source that can provide DC 5 V and some interface cards for various task requirement. The physical simulation system is used for overall path programming, obstacle avoidance and motion control, and target identification.

Path following control laws design of $τ_{u}$ and $τ_{r}$

A novel look-ahead distance is designed in this article for UMV path following control system, which is different from traditional LOS guidance algorithm:

Δ = \frac{\sqrt{λ_{1} + {(λ_{2} x_{e})}^{2}}}{λ_{3}}

where Δ is the look-ahead distance, λ₁, λ₂, and λ₃ are all positive designed control parameters, then the improved LOS guidance algorithm for path following could be described as follows

φ_{los} = - arctan \frac{y_{e}}{Δ} = - arctan \frac{λ_{3} y_{e}}{\sqrt{λ_{1} + {(λ_{2} x_{e})}^{2}}}

where $φ_{los}$ is the improved LOS guidance angle. To speed up the heading angle error $φ_{e}$ converges to 0, the UMV heading angle error’s virtual control input for path following control system could be designed as follows

α_{φ_{e}} = φ_{los}

To further simplify the path following error system (9), the following new intermediate state variables could be designed

{\begin{cases} ϖ_{1} = z_{4 r} - z_{4 r d} \\ ϖ_{2} = φ_{e} - α_{φ_{e}} = z_{3} - z_{3 s f} + β - α_{φ_{e}} \\ ϖ_{3} = {\dot{ϖ}}_{2} = {\dot{φ}}_{e} - {\dot{α}}_{φ_{e}} = {\dot{z}}_{3} - {\dot{z}}_{3 s f} + \dot{β} - {\dot{α}}_{φ_{e}} \\ = z_{6} - {\dot{z}}_{3 s f} + \dot{β} - {\dot{α}}_{φ_{e}} \end{cases}

where $z_{4 r d} = u_{r d} = u_{d} - u_{c}$ . Then, differentiate both sides of new intermediate state variables of equations (13), which results in the following new error system (14)

{\begin{cases} {\dot{ϖ}}_{1} = m_{22} z_{5 r} z_{6} / m_{11} - (d_{11}^{l} + d_{11}^{n l}) z_{4 r} / m_{11} + τ_{u} / m_{11} - {\dot{z}}_{4 r d} \\ {\dot{ϖ}}_{2} = ϖ_{3} \\ {\dot{ϖ}}_{3} = (m_{11} - m_{22}) m_{22} z_{4 r} z_{5 r} / (m_{22} m_{33} - m_{23}^{2}) + [(d_{22}^{l} + d_{22}^{n l}) m_{23} - (d_{32}^{l} \\ + d_{32}^{n l}) m_{22}] z_{5 r} / (m_{22} m_{33} - m_{23}^{2}) - [(d_{22}^{l} + d_{22}^{n l}) m_{23}^{2} - (d_{23}^{l} + d_{23}^{n l}) m_{22} m_{23} \\ - (d_{32}^{l} + d_{32}^{n l}) m_{22} m_{23} + (d_{33}^{l} + d_{33}^{n l}) m_{22}^{2}] z_{6} / [m_{22} (m_{22} m_{33} - m_{23}^{2})] \\ - m_{22} τ_{r} / (m_{22} m_{33} - m_{23}^{2}) - {\dot{z}}_{3 s f} + \dot{β} - {\dot{α}}_{φ_{e}} \end{cases}

Analyzing the above error system (14) by the state feedback linearization method and nonlinear backstepping technique and combing with the Lyapunov theory, the following control law (15) for $τ_{u}$ could be proposed

τ_{u} = - m_{22} z_{5 r} z_{6} + (d_{11}^{l} + d_{11}^{n l}) z_{4 r} + m_{11} {\dot{z}}_{4 r d} - k_{1} m_{11} ϖ_{1}

where $k_{1}$ is a positive designed control parameter. The control law $τ_{u}$ is the longitudinal thrust provided by propellers, and so on. It is a control quantity calculated by strict mathematical calculation under the condition of ensuring the local stability of the system, first equation of system (14).

According to the second and third equations of error system (14), and also combing with the feedback linearization method, nonlinear backstepping technique, and Lyapunov theory, the path following control law for $τ_{r}$ could be put forward as follows

\begin{matrix} τ_{r} = - (m_{11} - m_{22}) z_{4 r} z_{5 r} - [(d_{22}^{l} + d_{22}^{n l}) m_{23} - (d_{32}^{l} + d_{32}^{n l}) m_{22}] z_{5 r} / m_{22} + [(d_{22}^{l} + d_{22}^{n l}) m_{23}^{2} - (d_{23}^{l} + d_{23}^{n l}) m_{22} m_{23} - (d_{32}^{l} + d_{32}^{n l}) m_{22} m_{23} + (d_{33}^{l} + d_{33}^{n l}) m_{22}^{2}] z_{6} / m_{22}^{2} - k_{2} (m_{22} m_{33} - m_{23}^{2}) ϖ_{2} / m_{22} - k_{3} (m_{22} m_{33} - m_{23}^{2}) ϖ_{3} / m_{22} \end{matrix}

where $k_{2}$ and $k_{3}$ are both positive designed control parameters. The control law $τ_{r}$ is the lateral control moment provided by rudder, and so on. It is also a control quantity calculated by strict mathematical calculation under the condition of ensuring the local stability of the system, third equation of system (14).

Path following control law design of $\dot{s}$

Intermediate error variable ϖ₂ could converge to 0 under the path following control laws of $τ_{u}$ and $τ_{r}$ designed in equations (15) and (16), which would be proved in the ‘Stability analysis of path following control algorithm’ section. Combing with the LaSalle invariance principle in Isidori,⁴⁰ the first equation of system (9) could be described as follows

\begin{matrix} {\dot{x}}_{e} = K (s) \dot{s} y_{e} - \dot{s} + z_{4 r} cos (z_{3} - z_{3 s f}) - z_{5 r} sin (z_{3} - z_{3 s f}) + z_{4 c} cos (z_{3} - z_{3 s f}) - z_{5 c} sin (z_{3} - z_{3 s f}) \\ = K (s) \dot{s} y_{e} - \dot{s} + z_{4 r} cos (α_{φ_{e}} - β) - z_{5 r} sin (α_{φ_{e}} - β) + z_{4 c} cos (α_{φ_{e}} - β) - z_{5 c} sin (α_{φ_{e}} - β) \\ = K (s) \dot{s} y_{e} - \dot{s} + z_{4 r} (cos α_{φ_{e}} cos β + sin α_{φ_{e}} sin β) - z_{5 r} (sin α_{φ_{e}} cos β - cos α_{φ_{e}} sin β) \\ + z_{4 c} (cos α_{φ_{e}} cos β + sin α_{φ_{e}} sin β) - z_{5 c} (sin α_{φ_{e}} cos β - cos α_{φ_{e}} sin β) \end{matrix}

To simplify the system description, define the following variables

{\begin{cases} Ω = \sqrt{λ_{1} + {(λ_{2} x_{e})}^{2}} \\ Λ = \sqrt{{(z_{4 r} + z_{4 c})}^{2} + {(z_{5 r} + z_{5 c})}^{2}} \\ Θ = \sqrt{λ_{1} + {(λ_{2} x_{e})}^{2} + {(λ_{3} y_{e})}^{2}} \end{cases}

As it could be seen that

{\begin{cases} sin α_{φ_{e}} = - λ_{3} y_{e} / \sqrt{λ_{1} + {(λ_{2} x_{e})}^{2} + {(λ_{3} y_{e})}^{2}} = - λ_{3} y_{e} / Θ \\ cos α_{φ_{e}} = \sqrt{λ_{1} + {(λ_{2} x_{e})}^{2}} / \sqrt{λ_{1} + {(λ_{2} x_{e})}^{2} + {(λ_{3} y_{e})}^{2}} = Ω / Θ \\ sin β = (z_{5 r} + z_{5 c}) / \sqrt{{(z_{4 r} + z_{4 c})}^{2} + {(z_{5 r} + z_{5 c})}^{2}} = (z_{5 r} + z_{5 c}) / Λ \\ cos β = (z_{4 r} + z_{4 c}) / \sqrt{{(z_{4 r} + z_{4 c})}^{2} + {(z_{5 r} + z_{5 c})}^{2}} = (z_{4 r} + z_{4 c}) / Λ \end{cases}

Then combing with expressions (18) and (19), equation (17) could be rewritten as follows

\begin{array}{l} {\dot{x}}_{e} = K (s) \dot{s} y_{e} - \dot{s} + z_{4 r} [Ω (z_{4 r} + z_{4 c}) / (Λ Θ) - λ_{3} y_{e} (z_{5 r} + z_{5 c}) / (Λ Θ)] - z_{5 r} [- λ_{3} y_{e} (z_{4 r} + z_{4 c}) \\ / (Λ Θ) - Ω (z_{5 r} + z_{5 c}) / (Λ Θ)] + z_{4 c} [Ω (z_{4 r} + z_{4 c}) / (Λ Θ) - λ_{3} y_{e} (z_{5 r} + z_{5 c}) / (Λ Θ)] \\ - z_{5 c} [- λ_{3} y_{e} (z_{4 r} + z_{4 c}) / (Λ Θ) - Ω (z_{5 r} + z_{5 c}) / (Λ Θ)] \\ = K (s) \dot{s} y_{e} - \dot{s} + Ω (z_{4 r}^{2} + z_{4 r} z_{4 c} + z_{5 r}^{2} + z_{5 r} z_{5 c} + z_{4 c} z_{4 r} + z_{4 c}^{2} + z_{5 c} z_{5 r} + z_{5 c}^{2}) / (Λ Θ) \\ + λ_{3} y_{e} (- z_{4 r} z_{5 r} - z_{4 r} z_{5 c} + z_{5 r} z_{4 r} + z_{5 r} z_{4 c} - z_{4 c} z_{5 r} - z_{4 c} z_{5 c} + z_{5 c} z_{4 r} + z_{5 c} z_{4 c}) / (Λ Θ) \\ = K (s) \dot{s} y_{e} - \dot{s} + Ω [{(z_{4 r} + z_{4 c})}^{2} + {(z_{5 r} + z_{5 c})}^{2}] / (Λ Θ) \\ = K (s) \dot{s} y_{e} - \dot{s} + Λ Ω / Θ \end{array}

Analyzing the above error system (20) by the state feedback linearization method and nonlinear backstepping technique, combing with the Lyapunov theory, the virtual path following control law for $\dot{s}$ could be designed as follows

\dot{s} = Ω Λ / Θ + λ_{4} x_{e} Λ / Θ = (Ω + λ_{4} x_{e}) Λ / Θ

where λ₄ is a positive designed control parameter.

Stability analysis of path following control algorithm

In this section, stability of the path following control algorithm would be proved, Lyapunov stability theory and LaSalle invariance principle are adopted during the whole analysis process.

Stability analysis of intermediate state variables ϖ₁, ϖ₂, ϖ₃

Firstly, we discuss the stability and convergence of the intermediate state variables ϖ₁, ϖ₂, ϖ₃ designed in equation (13).

Substituting the path following control laws designed in equations (15) and (16) into the new error system (14) leads to

{\begin{matrix} {\dot{ϖ}}_{1} = - k_{1} ϖ_{1} \\ {\dot{ϖ}}_{2} = ϖ_{3} \\ {\dot{ϖ}}_{3} = - k_{2} ϖ_{2} - k_{3} ϖ_{3} \end{matrix}

Considering the following Lyapunov function candidate (23) for error system (22)

V_{1} = \frac{1}{2} ϖ_{1}^{2} + \frac{1}{2} k_{2} ϖ_{2}^{2} + \frac{1}{2} ϖ_{3}^{2}

Differentiating both sides of the above Lyapunov function candidate (23), combing with error system (22), it is easy to obtain the following differential expression (24)

\begin{array}{l} {\dot{V}}_{1} = ϖ_{1} {\dot{ϖ}}_{1} + k_{2} ϖ_{2} {\dot{ϖ}}_{2} + ϖ_{3} {\dot{ϖ}}_{3} \\ = - k_{1} ϖ_{1}^{2} + k_{2} ϖ_{2} ϖ_{3} - k_{2} ϖ_{2} ϖ_{3} - k_{3} ϖ_{3}^{2} \\ = - k_{1} ϖ_{1}^{2} - k_{3} ϖ_{3}^{2} \\ \leq 0 \end{array}

According to the above Lyapunov function candidate (23) and its differential expression (24), combing with the Lyapunov stability theory, it is easy to know that all the intermediate state variables ϖ₁, ϖ₂, ϖ₃ would globally asymptotically converge to 0 as $t \to \infty$ .

Stability analysis of path following errors x_e, y_e

In the “Stability analysis of intermediate state variables ϖ₁, ϖ₂, ϖ₃” section, ϖ₂ converges to 0 is proved; according to the LaSalle invariance principle in Isidori,⁴⁰ the second equation of system (9) could be described as follows

\begin{array}{l} {\dot{y}}_{e} = - K (s) \dot{s} x_{e} + z_{4 r} sin (z_{3} - z_{3 s f}) + z_{5 r} cos (z_{3} - z_{3 s f}) + z_{4 c} sin (z_{3} - z_{3 s f}) + z_{5 c} cos (z_{3} - z_{3 s f}) \\ = - K (s) \dot{s} x_{e} + z_{4 r} sin (α_{φ_{e}} - β) + z_{5 r} cos (α_{φ_{e}} - β) + z_{4 c} sin (α_{φ_{e}} - β) + z_{5 c} cos (α_{φ_{e}} - β) \\ = - K (s) \dot{s} x_{e} + z_{4 r} (sin α_{φ_{e}} cos β - cos α_{φ_{e}} sin β) + z_{5 r} (cos α_{φ_{e}} cos β + sin α_{φ_{e}} sin β) \\ + z_{4 c} (sin α_{φ_{e}} cos β - cos α_{φ_{e}} sin β) + z_{5 c} (cos α_{φ_{e}} cos β + sin α_{φ_{e}} sin β) \end{array}

Then combing with expressions (18) and (19), equation (25) could be rewritten as follows

\begin{array}{l} {\dot{y}}_{e} = - K (s) \dot{s} x_{e} + z_{4 r} [- λ_{3} y_{e} (z_{4 r} + z_{4 c}) / (Θ Λ) - Ω (z_{5 r} + z_{5 c}) / (Θ Λ)] + z_{5 r} [Ω (z_{4 r} + z_{4 c}) \\ / (Θ Λ) - λ_{3} y_{e} (z_{5 r} + z_{5 c}) / (Θ Λ)] + z_{4 c} [- λ_{3} y_{e} (z_{4 r} + z_{4 c}) / (Θ Λ) - Ω (z_{5 r} + z_{5 c}) \\ / (Θ Λ)] + z_{5 c} [Ω (z_{4 r} + z_{4 c}) / (Θ Λ) - λ_{3} y_{e} (z_{5 r} + z_{5 c}) / (Θ Λ)] \\ = - K (s) \dot{s} x_{e} + λ_{3} y_{e} (- z_{4 r}^{2} - z_{4 r} z_{4 c} - z_{5 r}^{2} - z_{5 r} z_{5 c} - z_{4 c} z_{4 r} - z_{4 c}^{2} - z_{5 c} z_{5 r} - z_{5 c}^{2}) / (Θ Λ) \\ + Ω (- z_{4 r} z_{5 r} - z_{4 r} z_{5 c} + z_{5 r} z_{4 r} + z_{5 r} z_{4 c} - z_{4 c} z_{5 r} - z_{4 c} z_{5 c} + z_{5 c} z_{4 r} + z_{5 c} z_{4 c}) / (Θ Λ) \\ = - K (s) \dot{s} x_{e} + λ_{3} y_{e} [- {(z_{4 r} + z_{4 c})}^{2} - {(z_{5 r} + z_{5 c})}^{2}] / (Θ Λ) \\ = - K (s) \dot{s} x_{e} - λ_{3} y_{e} Λ / Θ \end{array}

Substituting the virtual path following control laws designed in equation (21) into the UMV path following error equations (20) and (26) results in:

{\begin{cases} {\dot{x}}_{e} = K (s) y_{e} [(Ω + λ_{4} x_{e}) Λ / Θ] - λ_{4} x_{e} Λ / Θ \\ {\dot{y}}_{e} = - K (s) x_{e} [(Ω + λ_{4} x_{e}) Λ / Θ] - λ_{3} y_{e} Λ / Θ \end{cases}

Considering the following Lyapunov function candidate (28) for UMV path following error system (27):

V_{2} = \frac{1}{2} x_{e}^{2} + \frac{1}{2} y_{e}^{2}

Differentiating both sides of above Lyapunov function candidate (28), combing with error system (27), it is easy to obtain the following differential expression (29):

\begin{array}{l} {\dot{V}}_{2} = x_{e} {\dot{x}}_{e} + y_{e} {\dot{y}}_{e} \\ = K (s) x_{e} y_{e} [(Ω + λ_{4} x_{e}) Λ / Θ] - λ_{4} Λ x_{e}^{2} / Θ \\ - K (s) x_{e} y_{e} [(Ω + λ_{4} x_{e}) Λ / Θ] - λ_{3} Λ y_{e}^{2} / Θ \\ = - λ_{4} Λ x_{e}^{2} / Θ - λ_{3} Λ y_{e}^{2} / Θ \\ \leq 0 \end{array}

According to above Lyapunov function candidate (28) and its differential expression (29), combing with the Lyapunov stability theory, it is easy to know that both of the UMV path following errors x_e and y_e would globally and asymptotically converge to 0 as $t \to \infty$ .

Stability analysis of the whole UMV path following control system

As the UMV path following errors x_e and y_e both converge to 0, as proved in the “Stability analysis of path following errors x_e, y_e” section, according to expression (11), LOS guidance angle $φ_{l o s}$ would globally and asymptotically converge to 0 as $t \to \infty$ . Combing with expression (12), the UMV heading angle error’s virtual control input $α_{φ_{e}}$ would also globally and asymptotically converge to 0 as $t \to \infty$ .

In addition, as ϖ₁ converges to 0, as proved above in the “Stability analysis of intermediate state variables ϖ₁, ϖ₂, ϖ₃” section, then the first equation of system (13) could be rewritten as

ϖ_{1} = z_{4 r} - z_{4 r d} = u_{r} - u_{r d} = (u - u_{c}) - (u_{d} - u_{c}) = u - u_{d} = 0

So the longitudinal velocity u converges to the desired velocity u_d is proved here.

Furthermore, as ϖ₂ converges to 0, as proved in the “Stability analysis of intermediate state variables ϖ₁, ϖ₂, ϖ₃” section, and the virtual control input $α_{φ_{e}}$ is also proved to converge to 0, then combing with the second equation of system (13), it is easy to obtain that the UMV path following heading error $φ_{e}$ would also globally and asymptotically converge to 0 as $t \to \infty$ .

In summary, the path following control laws $τ_{u}$ , $τ_{r}$ designed in equations (15) and (16), proposed for an underactuated UMV in a preplanned route with the disturbance of ocean current, could have the path following errors $p_{e} = {[x_{e}, y_{e}, φ_{e}]}^{Τ}$ converge to 0, and make the longitudinal velocity u converges to desired velocity u_d. The UMV path following control system designed based on the control algorithm proposed in this article will have good stability, and it will help improve the accuracy and stability of ocean exploration and survey.

Simulation experiments

To verify and illustrate the effectiveness and reliability of the path following control system designed based on the control algorithm proposed in this article, several computer simulation experiments are carried out, for an underactuated UMV in a preplanned route with the disturbance of ocean current, on the semi-physical simulation platform in our laboratory. The UMV model is adopted with the following hydrodynamic parameters as presented in Table 1.

Table 1.

UMV model hydrodynamic parameters.

Parameter	Value	Parameter	Value	Parameter	Value
m ₁₁	25.8	m ₂₂	33.8	m ₃₃	2.76
m ₂₃	6.2	m ₃₂	6.2	$d_{11}^{l}$	27
$d_{11}^{n l}$	2.5 $\| u \|$	$d_{22}^{l}$	17	$d_{22}^{n l}$	4.5 $\| v \|$
$d_{33}^{l}$	0.5	$d_{33}^{n l}$	0.1 $\| r \|$	$d_{23}^{l}$	0.2
$d_{23}^{n l}$	0	$d_{32}^{l}$	0.5	$d_{32}^{n l}$	0

UMV: unmanned marine vehicle.

The control parameters are designed as presented in Table 2.

Table 2.

UMV path following control parameters.

Parameter	Value	Parameter	Value	Parameter	Value
$k_{1}$	1.5	$k_{2}$	2	$k_{3}$	5
λ ₁	100	λ ₂	2	λ ₃	1
λ ₄	1.5

UMV: unmanned marine vehicle.

A series of simulation experiment cases are carried out as follows:

Case 1

The preplanned route is assuming as a straight line path: $x_{s f} = \sqrt{2} s / 2$ , $y_{s f} = \sqrt{2} s / 2$ and the desired velocity of UMV is assuming as $u_{d} = 5 m/s$ , the velocity of constant ocean current disturbances is assuming as $V_{c} = 1 m/s$ with the current direction angle $β_{c} = 45^{\circ}$ , while the initial state of the underactuated UMV is chosen as $x (0) = 20 m$ , $y = 40 m$ , $φ = 0$ , $u (0) = 1 m/s$ , $v (0) = 0$ , and $r (0) = 0$ . To show the advantages of the control algorithm proposed in this article, a comparative algorithm from Liao and Zhang⁴¹ is chosen as the comparative experiment, which is described as algorithm A in the following results, and the control algorithm proposed in this article is described as algorithm B, all the same in other simulation experiment cases. Simulation results of the path following control of an underactuated UMV in a preplanned straight line path with the disturbance of ocean current are shown in Figure 4, algorithms proposed in this article and in Liao and Zhang⁴¹ are respectively adopted here.

Figure 4.

Path following control results of an underactuated UMV in a preplanned straight line path with the disturbance of ocean current. UMV: unmanned marine vehicle.

And the detail convergence response curve of the path following errors and the actual UMV control input curve in the preplanned straight line path with the disturbance of ocean current are shown in Figure 5.

Figure 5.

Path following errors convergence results and UMV’s actual control input of Case 1. (a) Convergence results of UMV position errors x_e and y_e. (b) Convergence results of UMV heading angle errors and velocity errors $φ_{e}$ and u_e. (c) UMV control input curve. UMV: unmanned marine vehicle.

Case 2

To fully demonstrate the effectiveness and reliability of the control algorithm proposed in this article, assuming that in a 2000 × 2000 m² area which need to be surveyed, the path is preplanned, the algorithm proposed in this article is adopted to achieve the UMV path following control of the preplanned survey path, the initial state of the underactuated UMV is chosen as $x (0) = 0$ , $y = 0$ , $φ = 90^{\circ}$ , the ocean current disturbance and other initial state variables’ values are same with Case 1, and the detail path following control results are shown in Figure 6.

Figure 6.

Path following control results of an underactuated UMV survey in a preplanned 2000 × 2000 m² area with the disturbance of ocean current. UMV: unmanned marine vehicle.

Case 3

Considering some preplanned path maybe curve path during the UMV survey process, a curve path is chosen for the following as $x_{s f} = 200 cos (0.005 s)$ and $y_{s f} = 200 sin (0.005 s)$ , the initial state of the underactuated UMV is chosen as $x (0) = 0$ , $y = 0$ , $φ = 0$ , while the ocean current disturbance and the UMV’s other initial state variables’ values are same with Case 1. Algorithms proposed in this article and in Liao and Zhang⁴¹ are also, respectively, adopted here, and simulation results of the path following control of an underactuated UMV survey in a preplanned curve path with the disturbance of ocean current are shown in Figure 7.

Figure 7.

Path following control results of an underactuated UMV survey in a preplanned curve path with the disturbance of ocean current. UMV: unmanned marine vehicle.

The detailed convergence response curve of the path following errors and the actual UMV control input curve during the survey in the preplanned curve path with the disturbance of ocean current are shown in Figure 8.

Figure 8.

Path following error convergence results and UMV’s actual control input of Case 3. (a) Convergence results of UMV position errors x_e and y_e. (b) Convergence results of UMV heading angle errors and velocity errors $φ_{e}$ and u_e. (c) UMV control input curve. UMV: unmanned marine vehicle.

Case 4

To fully demonstrate the effectiveness and reliability of the control algorithm proposed in this article, assuming that in a 600 × 600 m² area which need to be surveyed, the survey path is preplanned as three circle path as shown in Figure 8, the algorithm proposed in this article is adopted to achieve the UMV path following control of the preplanned path, the ocean current disturbance and the UMV’s initial state variables value are the same with Case 3, the detailed path following control results are shown in Figure 9.

Figure 9.

Path following control results of an underactuated UMV survey in a preplanned 600 × 600 m² area with the disturbance of ocean current. UMV: unmanned marine vehicle.

Discussion

Path following control problem of an underactuated UMV sailing in the presence of ocean current disturbances in a preplanned route is addressed in this article. By analyzing the existing research results described in the Introduction, we can find that, at present, most of the research results on path following of an UMV are based on fully symmetric UMV model, meanwhile the damping coefficient higher order velocity coupling terms of UMV model are not considered. While considering the higher accuracy requirements of the UMV’s path following control system, a novel path following control system, which considers the asymmetry and the damping coefficient higher order velocity coupling terms of the UMV model, is designed in this article based on an improved LOS guidance algorithm for and underactuated UMV in the presence of ocean current disturbances.

Simulation results are shown in Figures 3 to 8, which all verified the effectiveness and reliability of the path following control system designed for an underactuated UMV survey in a preplanned route with the disturbance of constant ocean current. And results in Figures 3 and 6 show that, the control algorithm proposed in this article is better than that of reference,⁴¹ as UMV model’s asymmetry and the damping coefficient higher order velocity coupling terms of UMV model are both not considered in Liao and Zhang.⁴¹ It demonstrates that the UMV model’s asymmetry and the damping coefficient higher order velocity coupling terms of UMV model could affect the control precision of the UMV control system, so they should be considered during the accuracy UMV ocean exploration and survey.

An improved LOS guidance law is proposed, by combing with the longitudinal path following error and lateral path following error, the tracking errors converge faster than traditional LOS method as the control parameters designed could speed up the convergence of the path following control system. In addition, as both longitudinal path following error and lateral path following error are integrated in the LOS guidance law, it means the guidance angle could change with the path following errors quickly, and this leads to that, the improved LOS guidance algorithm could achieve both straight line path following control and curve path following control, while traditional LOS algorithm could only achieve straight line path following control.

Virtual control input design for $\dot{s}$ and other intermediate state variables design during the control algorithm process, simplified the UMV control system, especially when asymmetry of the UMV model is considered, the difficulties that nonzero term appears in the non-diagonal elements of the system inertia coefficient matrix and damping coefficient matrix, meanwhile higher order velocity coupling terms of damping coefficients of the UMV are discussed, which include both linear damping and nonlinear damping. Several ingenious intermediate process designs and a kind of global differential homeomorphism transformation design are presented for dealing with the complex nonlinear terms and difficulties caused by the model asymmetry of UMV. A series of computer simulation experiments carried out for an underactuated UMV survey in a preplanned route with the disturbance of ocean current, on the semi-physical simulation platform in our laboratory, demonstrate the effectiveness and superiority of control algorithm designed.

Future research will be focused on path following control of the underactuated UMV in the presence of random environment disturbances and UMV model uncertainty.

Conclusions

Path following control problem of an underactuated UMV sailing in a preplanned route with the disturbance of ocean current is discussed, an improved LOS guidance algorithm is proposed by integrating both longitudinal path following error and lateral path following error into the guidance angle. Asymmetry of the UMV model is considered and a kind of global differential homeomorphism transformation design is proposed to solve the difficulty of nonzero term appearing in the non-diagonal elements of the system inertia coefficient matrix and damping coefficient matrix, which is caused by the model asymmetry of UMV. Higher order velocity coupling terms of damping coefficients of the UMV are discussed, which include both linear damping and nonlinear damping, an ingenious intermediate process design is presented for dealing with the complex nonlinear terms. Virtual velocity in the tangent direction of the path to be followed is designed for the control system, which increases a virtual control input and leads to simplify the underactuated control system. Path following control algorithm designed in this article could achieve both straight line path and curve path following, while some algorithm could only achieve straight line path or curve path following. Ocean current disturbance is considered in this article, and the ocean current relative coordinate system and ocean current relative velocity are designed in this article to integrate the ocean current disturbance into the kinematics equations and dynamics equations of UMV. Stability of the proposed path following control algorithm is proved by Lyapunov theory and LaSalle invariant principle, and a series of simulation experiments are carried out on the semi-physical simulation platform in our laboratory, which demonstrates the effectiveness and superiority of control algorithm designed. While path following control of the underactuated UMV in the presence of random environment disturbances and UMV model uncertainty will be studied in the near future.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Pre-research of Equipment Project (41412030201), the National Natural Science Foundation of China (51779057, 51509054, 51709061, 51709214), the Innovation Special Zone of National Defense Science and Technology.

ORCID iD

Sun Yushan

References

Wang

Sun

Yin

, et al. Fuzzy unknown observer-based robust adaptive path following control of underactuated surface vehicles subject to multiple unknowns. Ocean Eng 2019; 176: 57–64.

Fzzani

Zerr

Mansour

, et al. Fusion of swath bathymetric data: application to AUV rapid environment assessment. IEEE J Ocean Eng 2019; 44: 111–120.

Dong

Yang

Zheng

, et al. Parameter identification of unmanned marine vehicle manoeuvring model based on extended kalman filter and support vector machine. Int. J Adv Robot Syst 2019; 16: 1–10.

Taheri

Ferdowsi

Danesh

, et al. Closed-loop randomized kinodynamic path planning for an autonomous underwater vehicle. Appl Ocean Res 2019; 83: 48–64.

Chen

Zhang

, et al. Underwater terrain position method using maximum a posteriori estimation and PCNN model. J Navig 2019; 72(5): 1233–1253. DOI: 10.1017/S0373463319000067.

Xiang

Zhang

, et al. Adaptive fuzzy trajectory tracking control of an underactuated autonomous underwater vehicle subject to actuator saturation. Int J Fuzzy Syst 2018; 20: 269–279.

Mohammad

Hamidreza

Alireza

, et al. A low-cost dead reckoning navigation system for an AUV using a robust AHRS: design and experimental analysis. IEEE J Ocean Eng 2018; 43: 927–939.

Dong

Bao

Zheng

, et al. Heading control of unmanned marine vehicles based on an improved robust adaptive fuzzy neural network control algorithm. IEEE Access 2019; 7: 9704–9713.

Freir

Oliveira

Vale

RTS

, et al. Development of an AUV control architecture based on systems engineering concepts. Ocean Eng 2018; 151: 157–169.

10.

Wang

Jiang

Chen

, et al. Trajectory tracking control of underwater vehicle-manipulator system using discrete time delay estimation. IEEE Access 2017; 5: 7435–7443.

11.

Liang

Hou

, et al. Three-dimensional trajectory tracking control of an underactuated autonomous underwater vehicle based on ocean current observer. Int J Adv Robot Syst 2018; 15(5): 1–9.

12.

Belleter

Maghenem

Paliotta

, et al. Observer based path following for underactuated marine vehicles in the presence of ocean currents: a global approach. Automatica 2019; 100: 123–134.

13.

Dong

Wang

Song

, et al. Robust exponential point stabilization control of the high-speed underactuated unmanned marine vehicle with model asymmetry. Int J Adv Robot Syst 2019; 16: 1–12.

14.

Wang

Tong

Wang

. High-gain observer-based line-of-sight guidance for adaptive neural path following control of underactuated marine surface vessel. IEEE Access 2019; 7: 26088–26101.

15.

Lamraoui

Zhu

. Path following control of fully actuated autonomous underwater vehicle based on LADRC. Pol Marit Res 2018; 25: 39–48.

16.

Liu

Dong

, et al. Path following control of the underactuated USV based on the improved line-of-sight guidance algorithm. Pol Marit Res 2017; 24: 3–11.

17.

King

Vardy

Forrest

. Teach-and-repeat path following for an autonomous underwater vehicle. J Field Robot 2018; 35(5): 748–763.

18.

Wang

Zhang

, et al. Finite-time observer based accurate tracking control of a marine vehicle with complex unknowns. Ocean Eng 2017; 145: 406–415.

19.

Zheng

Sun

Xie

. Error-constrained LOS path following of a surface vessel with actuator saturation and faults. IEEE Trans Syst Man Cybern Syst 2018; 48: 1794–1805.

20.

Kelasid

Kohl

Pettersen

, et al. Experimental investigation of locomotion efficiency and path-following for underwater snake robots with and without a caudal fin. Annu Rev Control 2018; 46: 281–294.

21.

Xiang

Lapierre

, et al. Nonlinear guidance and fuzzy control for three-dimensional path following of an underactuated autonomous underwater vehicle. Ocean Eng 2017; 146: 457–467.

22.

Campos

Monroy

Abundis

, et al. A nonlinear controller based on saturation functions with variable parameters to stabilize an AUV. Int J Nav Archit Ocean Eng 2019; 11: 211–224.

23.

Zhang

Huang

Qin

, et al. A novel adaptive second order sliding mode path following control for a portable AUV. Ocean Eng 2018; 151: 82–92.

24.

David

Bauer

Seto

. Coupled hydroplane and variable ballast control system for autonomous underwater vehicle altitude-keeping to variable seabed. IEEE J Ocean Eng 2018; 43: 873–887.

25.

Wang

Sun

, et al. Fuzzy uncertainty observer-based path-following control of underactuated marine vehicles with unmodeled dynamics and disturbances. Int J Fuzzy Syst 2018; 20: 2593–2604.

26.

Wang

Liu

. Globally stable adaptive dynamic surface control for cooperative path following of multiple underactuated autonomous underwater vehicles. Asian J Control 2018; 20: 1204–1220.

27.

Huo

Wang

. Horizontal path-following control for deep-sea work-class ROVs based on a fuzzy logic system. Ships Offshore Struct 2018; 13: 637–648.

28.

Huang

Zhang

, et al. Model based adaptive control and disturbance compensation for underwater vehicles. Chin J Mech Eng 2018; 31: 1–13.

29.

Peng

Wang

. Output-feedback path-following control of autonomous underwater vehicles based on an extended state observer and projection neural networks. IEEE Trans Syst Man Cybern—Syst 2018; 48: 535–544.

30.

Jiang

, et al. Route optimizing and following for autonomous underwater vehicle ladder surveys. Int J Adv Robot Syst 2018; 15: 1–16.

31.

Liang

Wan

, et al. Three-dimensional path following of an underactuated AUV based on fuzzy backstepping sliding mode control. Int J Fuzzy Syst 2018; 20: 640–649.

32.

King

Vardy

Forrest

. Teach-and-repeat path following for an autonomous underwater vehicle. J Field Robot 2018; 35: 748–763.

33.

Fezzani

Zerr

Mansour

, et al. Fusion of swath bathymetric data: application to AUV rapid environment assessment. IEEE J Ocean Eng 2019; 44: 1–10.

34.

Fernandes

Rodrigues

Neto

, et al. Modeling positional uncertainty for hydrographic surveys with AUV. J Surv Eng-ASCE 2018; 145: 04018014.

35.

Wang

Pan

, et al. Yaw-guided trajectory tracking control of an asymmetric underactuated Surface Vehicle. IEEE Trans Ind Inform 2019; 15: 3502–3513.

36.

Wang

Yan

Chen

, et al. A new adaptive time-delay control scheme for cable-driven manipulators. IEEE Trans Ind Inform 2019; 15: 3469–3481.

37.

Wang

Han

, et al. Backpropagating constraints based trajectory tracking control of a quad rotor with constrained actuator dynamics and complex unknowns. IEEE Trans Syst Man Cybern—Syst 2019; 49: 1322–1337.

38.

Wang

Yan

, et al. Practical adaptive fractional-order nonsingular terminal sliding mode control for a cable-driven manipulator. Int J Robust Nonlinear Control 2019; 29: 1396–1417.

39.

Fossen

TI.

Guidance and control of ocean vehicle. New York, NY: Wiley, 1994, pp.53–157.

40.

Isidori

. Nonlinear control systems. Berlin, Germany: Springer Science & Business Media, 1995, pp. 21–79.

41.

Liao

Zhang

. Serret-Frenet frame based on path following control for underactuated unmanned surface vehicles with dynamic uncertainties. J Cent South Univ 2015; 22: 214–223.

Path following control of an underactuated unmanned marine vehicle with model asymmetry in the presence of ocean current disturbances

Abstract

Keywords

Introduction

Problem formulation

UMV modeling

Global differential homeomorphism transformation

UMV path following error system

Path following controller design for an underactuated UMV

Path following control laws design of τ u and τ r

Path following control law design of s ˙

Stability analysis of path following control algorithm

Stability analysis of intermediate state variables ϖ1, ϖ2, ϖ3

Stability analysis of path following errors xe, ye

Stability analysis of the whole UMV path following control system

Simulation experiments

Case 1

Case 2

Case 3

Case 4

Discussion

Conclusions

Footnotes

Declaration of conflicting interests

Funding

ORCID iD

References

Path following control laws design of $τ_{u}$ and $τ_{r}$

Path following control law design of $\dot{s}$

Stability analysis of intermediate state variables ϖ₁, ϖ₂, ϖ₃

Stability analysis of path following errors x_e, y_e