*Article* **Ship Motion Planning for MASS Based on a Multi-Objective Optimization HA\* Algorithm in Complex Navigation Conditions**

**Meiyi Wu <sup>1</sup> , Anmin Zhang 1,2, Miao Gao 1,\* and Jiali Zhang 1,\***


**Abstract:** Ship motion planning constitutes the most critical part in the autonomous navigation systems of marine autonomous surface ships (MASS). Weather and ocean conditions can significantly affect their navigation, but there are relatively few studies on the influence of wind and current on motion planning. This study investigates the motion planning problem for USV, wherein the goal is to obtain an optimal path under the interference of the navigation environment (wind and current), and control the USV in order to avoid obstacles and arrive at its destination without collision. In this process, the influences of search efficiency, navigation safety and energy consumption on motion planning are taken into consideration. Firstly, the navigation environment is constructed by integrating information, including the electronic navigational chart, wind and current field. Based on the environmental interference factors, the three-degree-of-freedom kinematic model of USVs is created, and the multi-objective optimization and complex constraints are reasonably expressed to establish the corresponding optimization model. A multi-objective optimization algorithm based on HA\* is proposed after considering the constraints of motion and dynamic and optimization objectives. Simulation verifies the effectiveness of the algorithm, where an efficient, safe and economical path is obtained and is more in line with the needs of practical application.

**Keywords:** motion planning; MASS; multi-objective optimization; complex navigation conditions

## **1. Introduction**

As an unmanned intelligent marine carrier platform, the Unmanned Surface Vehicle (USV) is small in size, flexible to operate and of high security. It can be equipped with different sensors or weapon systems, as required, in order to perform various tasks in military and civilian fields [1,2]. Amid the continuous development of the global marine economy and the intensifying disputes over maritime rights and interests, USV is a technical driver which can not only promote the rapid and sound development of the marine economy, but can also boost the strength of marine equipment and safeguard maritime rights and interests [3,4]. In 2020, the European Maritime Safety Agency (EMSA) released an overview of 13,204 maritime casualties from 2014 to 2019 [5]. The data show that the main cause of the casualties was loss of control of the ships (31.4%), followed by ship collision/contact accidents (30.5%), as shown in Table 1. About 52.3% of all the maritime accidents investigated were caused by personnel misconduct. The research on intelligent collision avoidance decision of USV can effectively reduce the influence of human factors and human errors on the navigation safety of ships [6], and gradually some routine or high-risk manual operations can be replaced by the USV. At present, the study of unmanned surface vehicles has become a significant issue in the field of international maritime affairs, attracting the attention of more and more shipping and shipbuilding countries in the world and emerging as a very important development direction in the future shipbuilding industry.

**Citation:** Wu, M.; Zhang, A.; Gao, M.; Zhang, J. Ship Motion Planning for MASS Based on a Multi-Objective Optimization HA\* Algorithm in Complex Navigation Conditions. *J. Mar. Sci. Eng.* **2021**, *9*, 1126. https:// doi.org/10.3390/jmse9101126

Academic Editors: Haitong Xu, Lúcia Moreira and Carlos Guedes Soares

Received: 22 September 2021 Accepted: 9 October 2021 Published: 14 October 2021

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

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

*J.* 

future shipbuilding industry.

future shipbuilding industry.

future shipbuilding industry.

future shipbuilding industry.


**Table 1.** Distribution of casualty events with a ship. with a

*J. Mar. Sci. Eng.* **2021**, *9*, x FOR PEER REVIEW 2 of 22

*J. Mar. Sci. Eng.* **2021**, *9*, x FOR PEER REVIEW 2 of 22

*J. Mar. Sci. Eng.* **2021**, *9*, x FOR PEER REVIEW 2 of 22

*J. Mar. Sci. Eng.* **2021**, *9*, x FOR PEER REVIEW 2 of 22

Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of the unmanned surface vehicle to a certain extent and is also one of the bottleneck factors that restrict USVs in terms of achieving high autonomy at present [14]. Compared with the Unmanned Ground Vehicle (UGV) and the Unmanned Aerial Vehicle (UAV), when applied, the USV can be interfered with by wind, waves and currents. Complicated environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as such it is easy to make the USV roll over due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion response time are also longer than those of Motion is core achieving high in a highly dynamic and uncertain environment [12,13], represents the of the unmanned surface to a certain and is also one of the factors restrict USVs in of achieving high at Compared with the Unmanned Vehicle (UGV) and Vehicle when by wind, and Complicated environmental has on angle navigation, and such it is make the USV roll due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion longer than those of Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of the unmanned surface vehicle to a certain extent and is also one of the bottleneck factors that restrict USVs in terms of achieving high autonomy at present [14]. Compared with the Unmanned Ground Vehicle (UGV) and the Unmanned Aerial Vehicle (UAV), when applied, the USV can be interfered with by wind, waves and currents. Complicated environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as such it is easy to make the USV roll over due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion response time are also longer than those of Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of the unmanned surface vehicle to a certain extent and is also one of the bottleneck factors that restrict USVs in terms of achieving high autonomy at present [14]. Compared with the Unmanned Ground Vehicle (UGV) and the Unmanned Aerial Vehicle (UAV), when applied, the USV can be interfered with by wind, waves and currents. Complicated environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as such it is easy to make the USV roll over due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion response time are also longer than those of Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of the unmanned surface vehicle to a certain extent and is also one of the bottleneck factors that restrict USVs in terms of achieving high autonomy at present [14]. Compared with the Unmanned Ground Vehicle (UGV) and the Unmanned Aerial Vehicle (UAV), when applied, the USV can be interfered with by wind, waves and currents. Complicated environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as such it is easy to make the USV roll over due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion response time are also longer than those of path planning will be improved on the model generate the constraint conditions. Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of surface vehicle certain extent and one of bottleneck that restrict terms of achieving autonomy present [14]. with Ground Vehicle the Aerial applied, the can be interfered with by waves environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as it is to the over due excessive leeway angle angle rate [15,16]. In addition, mostly under-Motion planning is core to the USV achieving high autonomy in a highly dynamic and uncertain navigation environment [12,13], which represents the intelligence level of the unmanned surface vehicle to a certain extent and is also one of the bottleneck factors that restrict USVs in terms of achieving high autonomy at present [14]. Compared with the Unmanned Ground Vehicle (UGV) and the Unmanned Aerial Vehicle (UAV), when applied, the USV can be interfered with by wind, waves and currents. Complicated environmental disturbance has a great influence on its instantaneous speed and attitude angle during navigation, and as such it is easy to make the USV roll over due to the excessive leeway and drift angle or turning angle rate [15,16]. In addition, USVs use mostly underactuated systems, and their inertia and motion response time are also longer than those of UAVs and UGVs, which brings greater uncertainty to the control and motion planning of USVs in complex navigation environments. Therefore, the autonomy of the

USV is essential. This depends on two complex and changeable environmental parameters: wind and current [17,18]. Ignoring the environmental impact in motion planning would not only lead to a great waste of energy when the USV navigates strong ocean currents, but would also increase the potential risk of hitting obstacles. However, most of the current research results are based on idealized assumptions that do not consider whether there is environmental interference or it is a steady environment, and the expression of constraints for USV motion planning is inaccurate and incomplete, which limits its application in different scenarios.

In the field of USV intelligent planning and control, studies focusing on singleobjective motion planning, such as in terms of length, optimal sailing time, energy consumption, smoothness and safety, have gradually deepened knowledge [19,20]. Planning algorithms can generally be divided into graph search algorithms, random sampling algorithms, curve interpolation algorithms, machine learning and dynamic optimization methods [21]. In the practical application situation, the above methods are usually used in combination in order to complete the motion planning. Sang et al. [22] used the improved A\* algorithm to keep a safe distance and avoid collisions by reducing search points near obstacles. At the same time, the turning cost is added to the heuristic function in order to reduce the turning points of the path, avoid frequent turns of USV and improve the smoothness of the path. Liu et al. [23] introduced a safety parameter into the FMS algorithm, which can adaptively adjust the influence of the obstacle size, so as to ensure obstacle avoidance in a constrained environment and improve navigation safety. Zuo et al. [24] proposed A\*-LSPI hierarchical path planning method, by which the global path based on the A\* algorithm was quickly found and the approximate optimal local planning strategy with LSPI was learned. Although the planning time was long, the generated path length was short. Han et al. [25] formulated the multi-criteria global shortest path planning problem with resource constraints as a single objective linear programming model, and provided a modified label-correcting algorithm to solve this problem within a rationally short time. However, the planning environment was too simplified, and the environmental impact was less considered. Subramani et al. [26] formulated a stochastic optimization method to compute energy-optimal paths from among time-optimal paths of autonomous vehicles navigating in a dynamic flow field. Xu et al. [27] generated the path for an automatic ferry on the basis of AIS historical data. The resulting path was safer and more economic, because the AIS data were recorded from the real-time trajectory of ships. Lei et al. [28] proposed a multi-direction A\* algorithm to iteratively find an optimum neighbor node and APF in scalar mode, which can take into account both computational complexity and efficiency. Xu et al. [29,30] proposed the vector field guidance law for the path-following control problem of the underactuated surface ship, which considered straight-line and curved-path path following scenarios in the presence of ocean currents. R. Zaccone et al. [31] developed and proposed a ship voyage optimization method, aiming to find the voyage which would require minimum fuel consumption within safety and comfort constraints by using 3D Dynamic Programming optimization.

In conclusion, amid a complex navigation environment, motion planning considering multi-constraint control decision and multi-objective optimization has become a difficulty and represents a hotspot in USV motion planning research. Therefore, it is necessary to study the USV motion planning of multi-objective optimization with complex constraints in dynamic navigation environments, that is, to make reliable motion planning under the complex constraints of environments, kinematics, dynamics and optimization objects.

According to the planning process, motion planning involves three steps: the environment model, the motion mathematical model and the search algorithm. The structure of this paper is as follows: Firstly, an environmental model, including wind and current, is constructed in order to provide high quality navigation environmental information for USV navigation. Based on the MMG model, a motion mathematical model suitable for *Dolphin-I* USV is proposed. Combined with the above model, the cost function of multi-objective optimization is analyzed, and a motion planning algorithm MOHA\* is proposed in order

to solve the multi-objective problem in dynamic navigation environment. The algorithm is then used to simulate the motion planning of the USV. Finally, the reliability of the MOHA\* algorithm is proved by analyzing the experimental results. The chapter structure is shown in Figure 1, and the full text structure is shown in Figure 2.

**Figure 1.** Section structure diagram of this paper.

**Figure 2.** Full text structure diagram.

#### **2. Materials and Methods**

#### *2.1. Construction of Navigation Environment Based on ENC*

The establishment of the environment model includes the process of extracting and describing environmental characteristic information. The electronic navigational chart (ENC) classifies and stores relevant elements in the geometric form of points, lines, and planes in data files, which can display and select relevant marine environment information

according to need. This has such advantages as a short storage time, fast display speed and high accuracy [32]. Reading the overall package information in the ENC, and further processing the data in the forms of point, line and plane required for conversion, are the basis of navigation environment modeling. Modeling methods commonly used for USV motion planning mainly include grids [33], topology graphs [34], Voronoi diagrams [35], visibility graphs [36], and others. The environment model in this paper is based on ENC. The number of environmental obstacles with a complicated structure is large. It can be too complicated to describe the obstacle nodes by using topology graphs, Voronoi diagrams and visibility graphs. The grid method has a simpler data structure than other environmental modeling methods, which can reduce the complexity and calculation of the boundary processing of complex-shaped obstacles [37,38]. In this paper, grids are used to divide the ENC information. The consistent expression of grid ENC is the basis for improving the efficiency of the path search algorithm. The size of grid granularity determines the advantages and disadvantages of modeling, to a certain extent. The grid size is set to be 25 m × 25 m, comprehensively taking the minimum turning radius of USV (9.6 m), navigating and positioning error (5 m), safety buffer distance (5 m) and electronic navigation chart error (5 m) into consideration, and making sure it can complete the steering operation in a grid size area.

In this paper, Zhoushan islands are selected as the research area, where the longitude range is E120◦5502600~E123◦2903000, the latitude range is N29◦3301500~N32◦2805900, and the proportional scale is 1:2,000,000. The S-57 ENC in this area was transformed by Mercator projection through ArcMap, and the global static obstacle information, including obstructions such as land, islands and the seabed, was extracted and further transformed into a grid map, as shown in Figure 3. The white grids represent navigable areas, and the black grids represent obstructed areas.

**Figure 3.** S-57 environment modeling map. (**a**) S-57 ENC of the study area; (**b**) Grid navigation environment map.

The second version of the NCEP Climate Forecast System (CFSv2) was made operational at National Centers for Environment Prediction (NCEP) in March 2011 [39]. This version has upgrades to nearly all aspects of the data assimilation and forecast model components of the system. This paper adopts the CFSv2 data set (http://cfs.ncep.noaa.gov, accessed on 10 January 2020) as the input wind field data, which contains 0.2◦ , 0.5◦ , 1.0◦ and 2.5◦ horizontal resolution data, and updates the data at hourly intervals.

The current data are obtained from the global Hybrid Coordinate Ocean Model (HY-COM) and Navy Coupled Ocean Data Assimilation (NCODA) 1/12◦ analysis, with a time update frequency of three hours. This system is configured for the global ocean with HY-COM2.2 as the dynamical model and NCODA for data assimilation [40]. Having gradually become a mainstream global ocean circulation model in recent years [41], HYCOM can select appropriate vertical coordinates according to different navigation environments and thus can better capture the various physical processes of oceans [42]. The NCODA system uses the model forecast as a first guess in a multivariate optimal interpolation scheme and

assimilates available in-situ observations [43]. More details about the assimilation system can be found at http://hycom.org, accessed on 10 January 2020.

Due to the limitation of data resolution, the 0.2◦ × 0.2◦ wind field and the 1/12◦ × 1/12◦ current field which were adopted cannot cover all of the non-obstacle grids in the map. Therefore, it is necessary to perform interpolation based on known environment information (e.g., wind, current) to predict environmental information in non-numerical areas. In this paper, based on the grid map resolution of 25 m, the known data are calculated by bicubic interpolation, and the corresponding wind and flow field data are obtained.

The processed ENC information was fused with the wind field and current field information on 10 January 2020, and the results are shown in Figures 4 and 5. Environmental information was added to each non-obstacle grid of the grid map. Besides terrain information, information about the speed and direction of the wind and current were also stored in the grid. The wind reanalysis products cover the time period of the global HYCOM and NCODA assimilation. Therefore, the two types of data can be updated synchronously at a time interval of three hours, and a dynamic environment map can be obtained.

**Figure 4.** Wind field map of local area.

**Figure 5.** Current field map of local area.

#### *2.2. Mathematical Model of USV under the Influence of Wind and Current*

The mathematical model is the basis of USV motion simulation and control [44,45]. In actual navigation, the speed and direction were controlled by the longitudinal propulsion force generated by two propellers at the tail and the turning moment generated by the differential speed, with no lateral driving force.

Meanwhile, considering that the rolling direction of the hull itself is relatively stable, a three-degree-of-freedom plane motion model was established, including surge (longitudinal motion), sway (sideways motion), and yaw (rotation around the vertical axis). The navigation state is formulated by two coordinate systems: one is the *ooxoy<sup>o</sup>* inertial coordinate system demonstrating the absolute position/speed information of the unmanned surface vehicle, and the other is the *oxy* attached coordinate system which studies change in the status of the unmanned surface vehicle, taking the first-order differential . *x*, . *y*, . *θ* 

as the amount *V* = (*u*, *v*,*r*) *T* of status change of the unmanned surface vehicle, as shown in Figure 6. In this paper, *Dolphin-I* USV of Tianjin University is used as the experimental platform. It adopts a modular design which can install different modules according to the application scenarios. The detailed USV performance parameters are shown in Table 2.

**Figure 6.** Schematic diagram of USV of three degrees of freedom.

**Table 2.** USV performance parameter setting details.


The MMG model mainly works to decompose the hydrodynamic force and torque acting on the ship into the hydrodynamic force and torque acting on the bare hull, openwater propeller and open-water rudder according to the physical meaning, and the mutual interference fluid between them. Based on the hull performance (rudderless, double propellers) of *Dolphin-I* USV, the MMG separation modeling is adopted, and the origin of the appendage coordinate system is taken as the center of gravity of the unmanned surface vehicle, and the motion equation is:

$$\begin{cases} (m+m\_x)\dot{u} - (m+m\_y)vr = \mathcal{X}\_H + \mathcal{X}\_{2P} \\ (m+m\_y)\dot{v} + (m+m\_x)ur = \mathcal{Y}\_H + \mathcal{Y}\_{2P} \\ (I\_{zz} + I\_{zz})\dot{r} = \mathcal{N}\_H + \mathcal{N}\_{2P} \end{cases} \tag{1}$$

where *m* is the mass of USV; *m<sup>x</sup>* and *m<sup>y</sup>* are additional inertial masses in x-axis direction and y-axis direction respectively; *Izz* and *Jzz* are the rotational inertia torque and the additional inertia torque in the z-axis direction respectively; *X*, *Y* and *N* are external forces and torques, and subscripts *H* and *2P* respectively represent bare hull and two propellers. The additional mass, inertia torque, additional inertia torque, forces and torques acting on hull and propellers can be calculated using the calculation methods described in the literature [46].

The *Dolphin-I* USV is differentially driven by two propellers, that is, it's speed and direction are controlled by the speed difference between the double propellers. According to the thrust model of two brushless DC thrusters in USBV power system proposed by Jin et al. [47] and the USV dynamic model of rudderless dual thrusters proposed by Li et al. [48], the resultant thrust vector of differential drive USV is established as follows:

$$F\_{thrust} = \begin{bmatrix} X\_{2P} \\ Y\_{2P} \\ N\_{2P} \end{bmatrix} = \begin{bmatrix} F\_L + F\_R \\ 0 \\ d\_{LR}(F\_L - F\_R) \end{bmatrix} \tag{2}$$

in which *F<sup>L</sup>* and *F<sup>R</sup>* are the thrust produced by the left and right thrusters along the *x*-axis in the attached coordinate system, respectively, and *dLR* is the transverse distance from the centerline of the USV to the centerline of each thruster.

In this paper, the operational performance of the unmanned surface vehicle in the wind is studied, and the mathematical model of the USV under wind disturbance is established by using a wind tunnel test and approximate estimation. When the unmanned surface vehicle is sailing, the superstructure device is affected by the wind, leading to the deviation of course or operational difficulty. When the unmanned surface vehicle is sailing at low speed, it is impacted by the wind quite severely. The interference force of wind on the unmanned surface vehicle can be regarded as the superposition of the average wind pressure and the variable wind pressure. In this paper, only the average wind pressure *Fwind* = [*Xwind*,*Ywind*, *Nwind*] *T* is considered, and the calculation expression is:

$$\begin{cases} X\_{wind} = 0.5 \rho\_a A\_f \mathcal{U}\_R^2 \mathbb{C}\_{w \times}(\boldsymbol{a}\_R) \\ \quad \boldsymbol{Y}\_{wind} = 0.5 \rho\_a A\_s \mathcal{U}\_R^2 \mathbb{C}\_{w \boldsymbol{y}}(\boldsymbol{a}\_R) \\ \quad \boldsymbol{N}\_{wind} = 0.5 \rho\_a A\_s \mathcal{L} \mathcal{U}\_R^2 \mathbb{C}\_{w \boldsymbol{n}}(\boldsymbol{a}\_R) \end{cases} \tag{3}$$

where *ρ<sup>a</sup>* is the air density; *A<sup>f</sup>* and *A<sup>s</sup>* are the orthographic projection area and the side projection area above the waterline of the unmanned surface vehicle respectively; *U<sup>R</sup>* is the relative wind speed; *L* is the total length of the USV; *Cwy*(*αR*), *Cwy*(*αR*) and *Cwn*(*αR*) are, respectively, the wind pressure torque coefficient in the direction of *x*-axis and *y*-axis, and the wind pressure coefficient around the *z*-axis, which are calculated according to Isherwood formula [49].

When the motion model of the unmanned surface vehicle was being established under the current interference force, considering the unevenness of the horizontal upstream, the velocity of the current field in Zhoushan maritime space is divided into fields according to the interval of 0.005 m/s, that is, the uneven current field is divided into uniform current fields in different areas for modeling. The impact of water at any position in the uniform flow is the same, which will cause the unmanned surface vehicle to drift and interfere with its original posture and motion state. In the inertial coordinate system, the relationship between absolute current velocity *Vc*, absolute current direction *ψc*, and the absolute velocity of the unmanned surface vehicle is established, and the expression is:

$$\begin{cases} \ u = u\_r + V\_\mathcal{c} \cos(\psi\_\mathcal{c} - \psi) \\ \ v = v\_r + V\_\mathcal{c} \sin(\psi\_\mathcal{c} - \psi) \end{cases} \tag{4}$$

where *u<sup>r</sup>* and *v<sup>r</sup>* are, respectively, the longitudinal and lateral velocity of the unmanned surface vehicle relative to the current; *V<sup>c</sup>* is the absolute current velocity; and *ψ<sup>c</sup>* is the absolute current direction.

The relative velocity method is used in the hydrodynamic calculation, and considering its additional force, the current disturbance creates the current force *Fcurrent* = [*Xcurrent*,*Ycurrent*, *Ncurrent*] *T* as given below:

$$\begin{cases} X\_{current} = \left(m\_x - m\_y\right) V\_c \sin(\psi\_c - \psi) \\ \quad Y\_{current} = \left(m\_x - m\_y\right) V\_c \cos(\psi\_c - \psi) \\ \quad N\_{current} = 0 \end{cases} \tag{5}$$

Considering the environmental interference factors mentioned above, the overall stress analysis of the unmanned surface vehicle motion system is carried out, and a threedegree-of-freedom MMG kinematic model is established with the hull, two propellers and environmental interference force as a function:

$$\begin{cases} (m+m\_x)\dot{u} - (m+m\_y)vr = \mathbf{X}\_H + \mathbf{X}\_{2P} + \mathbf{X}\_{wind} + \mathbf{X}\_{current} \\ (m+m\_y)\dot{v} + (m+m\_x)ur = \mathbf{Y}\_H + \mathbf{Y}\_{2P} + \mathbf{Y}\_{wind} + \mathbf{Y}\_{current} \\ (I\_{zz} + I\_{zz})\dot{r} = \mathbf{N}\_H + \mathbf{N}\_{2P} + \mathbf{N}\_{wind} + \mathbf{N}\_{current} \end{cases} \tag{6}$$

The Runge-Kutta method is used to solve differential equation (6) to obtain the actual speed *V<sup>e</sup>* = [*u<sup>e</sup>* , *v<sup>e</sup>* ] *T* and corresponding heading angle of the USV under the influence of environmental factors. Under the interference of the wind and current, assuming that the initial surge of the USV is 5 m/s, the given wind speed is 6 m/s, the wind direction is 180◦ , the current velocity is 1 m/s, and the current direction is 0◦ , the gyrating motion experiment and direct speed stability experiment are carried out with MATLAB. The actual measurement of the gyrating motion of the *Dolphin-I* USV in the still water shows that the advance distance is 16.5 m, and the initial diameter of the gyration is 24.5 m, as shown by the yellow dotted lines in Figure 7a,c. Comparing the measured data with the simulation results, it is proved that the kinematic model can accurately describe the navigation motion of the unmanned surface vehicle. Figure 7a,c show that the influence of wind on the cycle trajectory is not significant, while the cycle trajectory under the influence of the ocean current shifts eastward with a larger amplitude. In Figure 7b,d show that the speed changes greatly within 50 s under the influence of environmental disturbance, and then tends to be stable over time.

#### *2.3. Dynamic Model of USV under the Influence of Wind and Current*

The dynamic model of the unmanned surface vehicle represents the change of its posture under the action of force and torque. In this paper, based on the rigid-body dynamics model proposed by Fossen [50], with the effect of environmental interference, a three-degree-of-freedom kinematic model of the *Dolphin-I* USV is established. The specific expression is:

$$F\_{thrust} + F\_{env} = M\dot{V} + \mathcal{C}(V)V + D(V)V$$

$$M = \begin{bmatrix} m + m\_x & 0 & 0 \\ 0 & m + m\_y & 0 \\ 0 & 0 & I\_{zz} + I\_{zz} \\ 0 & 0 & I\_{zz} + I\_{zz} \end{bmatrix} = \begin{bmatrix} M\_{11} & 0 & 0 \\ 0 & M\_{22} & 0 \\ 0 & 0 & M\_{33} \end{bmatrix}$$

$$\mathcal{C}(V) = \begin{bmatrix} 0 & 0 & -M\_{22}v \\ 0 & 0 & M\_{11}u \\ M\_{22}v & -M\_{11}u & 0 \end{bmatrix} \tag{7}$$

$$D(V) = -\begin{bmatrix} X\_u & 0 & 0 \\ 0 & Y\_v & 0 \\ 0 & 0 & N\_r \end{bmatrix}$$

in which *M* represents the inertia matrix, including the added mass parameters. *C* is the Coriolis and centripetal matrix. *D* is hydrodynamic drag matrix. *Xu*, *Yv*, and *N<sup>r</sup>* are collectively referred to as the hydrodynamic derivative, and the specific values of the above three variables are calculated by using the formula of literature [17]. *Fthrust*,*Fenv* are the thruster and environment forces (*Fwind*, *Fcurrent*), respectively, applied on the USV.

**Figure 7.** Motion simulation experiment of *Dolphin-I* USV under environmental interference. (**a**) The gyroscopic trajectory with or without wind; (**b**) The change of direct sailing speed in different wind directions; (**c**) The gyroscopic trajectory with or without current; (**d**) The change of direct navigation speed under different current directions.

#### **3. Algorithm of Ship Motion Planning**

Considering the spatial constraints and the constraints of planning behavior of the objects motion planning combines path planning with motion control [51]. Global motion planning refers to the calculation of the path from the departure point to the target point that meets certain performance requirements according to the established prior environmental map, such as shortest distance or the highest safety.

#### *3.1. Traditional Hybrid A\* Algorithm*

Hybrid A\* (HA\*) is an algorithm for UGV kinematics and it was first proposed by Stanford Laboratory [52] in 2008. It can perform the heuristic search in a continuous coordinate system and guarantee that the generated trajectory meets the vehicle nonholonomic constraint. This algorithm is a variant of the A\* algorithm, which adopts a four-dimensional search space and adds the orientation information of the mobile platform and the fourth dimension representing the forward and backward movement on the basis of the twodimensional plane, considering the final directions of the starting point and the end point. The core of the algorithm is to design a cost function for each node to be searched to determine the accessing sequence of each node in the search:

$$f(\mathbf{x}\_n) = \mathbf{g}(\mathbf{x}\_n) + h(\mathbf{x}\_n) \tag{8}$$

where *f*(*xn*) is the total cost estimate from the starting point *x<sup>s</sup>* through the current node *x<sup>n</sup>* to the goal node *xg*; the actual cost *g*(*xn*) from the starting point *x<sup>s</sup>* to the current node *xn*, and the heuristic cost estimate *h*(*xn*) from the current state *x<sup>n</sup>* to the goal point *xg*.

*g*(*xn*) can be calculated by recursive formula:

$$\mathbf{g}(\mathbf{x}\_n) = \mathbf{g}(\mathbf{x}\_{n-1}) + \varepsilon d(\mathbf{x}\_{n-1}, \mathbf{x}\_n) \times \text{DirectionCost} + \sigma |\mathbf{K}| \tag{9}$$

where *g*(*xn*−1) is the parent node of the current node, and *d*(*xn*−1, *xn*) is the Euclidean distance from the parent node to the current node. *DirectionCost* indicates the change of motion direction where the forward driving value is 1, and the reverse driving value is −1, which is used to ensure the forward driving of USV as much as possible; *K* is the curvature of the motion primitive; and *ε* and *σ* are weights, which are used to unify the order of magnitude of each item.

*h*(*xn*) involves two heuristic functions: (1) The 2D heuristic function with holonomic constraints. When the search node is far away from the goal point, the USV always moves towards the goal point, and its nonholonomic characteristics can be ignored. Based on the obstacles shown in the environment map, the heuristic function with complete constraints is only used to consider the position information of USV (2) and the 3D heuristic function with nonholonomic constraints. The motion range was set as an accessible region and discretized into three-dimensional grids. Only considering the motion constraints of USV, the optimal path from the center point of each grid to the goal point is calculated by using the Reeds-Shepp curve. This heuristic function does not depend on the grid map information at runtime, and can perform off-line calculations in advance. The state of the current point is then matched after simple transformation and rotation, thus improving the calculation efficiency. The current node heuristic value is the larger 3D heuristic value and 2D heuristic value, in order to ensure the reliability of the search algorithm.

#### *3.2. Multi-Objective Optimization Model of USV Motion Planning*

In traditional motion planning, only a single optimization objective is usually considered, which means that it is difficult to generate a high-quality path that satisfies spatial constraints, time constraints and motion constraints. In practical navigation, the global planning of the USV can be regarded as a multi-constraint and multi-objective optimization problem. The shortest search time can ensure the ability of the USV to respond in time. The safest path is the premise that the USV can perform various tasks. The most energyefficient route enables the USV to increase mileage. Therefore, under the condition that

space constraints and USV motion constraints are met, a model is proposed to optimize the search efficiency, safety performance and energy consumption.

#### 3.2.1. Graph Expansion/Search Model Based on Hybrid Motion Primitives

The traditional graph search algorithm uses 4-domain or 8-domain connection to expand nodes [53]. When the search map is expanded, the number of grids will increase rapidly, resulting in a sudden increase in the search time. Meanwhile, the obtained path has a large number of redundant nodes which do not meet nonholonomic constraints. Hence, firstly, the environment map is down-sampled, and the heuristic efficiency can be greatly improved by searching and calculating in the processed grid map. The heuristic function values of nodes in the original map are then restored by up-sampling.

When expanding nodes, continuous motion primitives—that is, trajectory segments that satisfy the motion constraints of USV—are used to ensure that the curvature of the path that is ultimately generated is less than the curvature corresponding to the minimum turning radius. In order to reduce the parameters needed by motion primitives and to reduce the computational complexity, circular arcs and straight lines are used to generate motion primitives. Accordingly, the short primitives are of better flexibility and stronger ability to bypass obstacles, but this will lead to an increase in the number of extended nodes. Long primitives can reach the vicinity of the final point faster, but they present a higher risk of collision with obstacles in complex conditions. Therefore, four kinds of motion primitives (forward and backward) are mixed in this paper. The long primitives and short primitives are divided with the grid size of 25 m as the base length, in which the short primitives are S1, S2, S10 , S20 , the long primitives are M1, M2, M10 , M20 , L1, L2, L10 , L20 , and the linear motion primitives are SL and SL0 . The specific length and curvature are shown in Table 3. This is achieved by taking full respective advantages of the above-mentioned motion primitives and setting different cost coefficients for them. The cost of long primitives is lower and the cost of short primitives is higher. This can make USV approach the goal point quickly in exposed waters, while ensuring that it passes through narrow obstacle areas. The expansion/search model based on hybrid motion primitives can effectively reduce the number of path nodes in planning, thus improving the computational efficiency and achieving drivability.


**Table 3.** Details of hybrid motion primitive parameters.

3.2.2. Risk Degree of Navigation Model Based on Ship Domain

Traditional collision detection takes the USV as a particle and expands the obstacle map. This method has high efficiency, but the expansion scale is difficult to choose and the accuracy cannot be guaranteed, resulting in a big difference between collision detection and real results. In actual navigation, it is necessary to maintain an exclusive domain around the USV which is defined as the navigation safety domain (NSD) [54,55] and is established in order to avoid encroachment by other ships or obstacles. NSD is usually ovalshaped, with its long axis being three to eight times the length of the ship. It is designed in order to delimit enough sea space for the USV to take actions to avoid collisions in advance. In this paper, combined with the quaternion ship domain (QSD) proposed by Wang et al. [56] and the basic navigation safety domain (BNSD) proposed by Zhou et al. [57], the ellipse model with four half axes in different directions is established with the USV as the origin coordinate, where the space domain can be divided into four subdomains Ω = {Ω1, Ω2, Ω3, Ω4}. To calculate and simplify the model's complexity, given the coefficients of overtaking encounters situations *s*(*i*) = 1, *t*(*i*) = 0.2, its mathematical expression is:

$$R = \begin{cases} R\_{force} = L + 1.34\sqrt{AD^2 + \left(DT/2\right)^2} \\ R\_{aff} = L + 0.67\sqrt{AD^2 + \left(DT/2\right)^2} \\ R\_{strab} = B + 1.2DT \\ R\_{port} = B + 0.9DT \end{cases} \tag{10}$$

In the above equations, *Rf ore*, *Ra f t*, *Rstarb* and *Rport* are the radii of the navigation safety domain. *L* and *B* represent the length and breadth of the USV. *AD* is the advance distance, the longitudinal forward distance of the gravity center in the case of the USV turning 90◦ from the start of steering. *DT* is the tactical diameter, the transverse distance of the gravity center in the case of the USV turning 180◦ from the start of steering.

Based on the above model, two safety domains *R*<sup>1</sup> = {30.74, 16.97, 31.60, 24.25} and *R*<sup>2</sup> = {46.74, 32.97, 42.60, 35.25} are constructed with one-time length and three-time length as input parameters, as shown in Figure 8.

**Figure 8.** Ship domain model of *Dolphin-I* USV.

According to the different heading angles, the reasonable radius range of the subsafety zone from USV to static obstacles is obtained with different heading angles and the radiuses of navigation safety zones: when 0 ◦ ≤ *<sup>ϕ</sup><sup>e</sup>* < <sup>90</sup>◦ , *<sup>R</sup>*(*xn*) ∈ <sup>Ω</sup>1; <sup>90</sup>◦ ≤ *<sup>ϕ</sup><sup>e</sup>* < <sup>180</sup>◦ , *<sup>R</sup>*(*xn*) ∈ <sup>Ω</sup>2; <sup>180</sup>◦ ≤ *<sup>ϕ</sup><sup>e</sup>* < <sup>270</sup>◦ , *<sup>R</sup>*(*xn*) ∈ <sup>Ω</sup>3; <sup>270</sup>◦ ≤ *<sup>ϕ</sup><sup>e</sup>* < <sup>360</sup>◦ , *R*(*xn*) ∈ Ω4. Collision risk index (CRI) is an evaluation parameter for the collision avoidance of ships, which is

also a key concept in the field of collision avoidance research and operation. Based on the collision risk index (CRI) calculation method of distance to closest point of approach (DCPA), the shortest safe distance of the USV and the actual distance *d*(*xn*) between the USV and obstacles are used to reflect the navigation risk. The smaller the value is, the more suitable it is for navigation, while the larger the value is, the less suitable it is for navigation. Its mathematical expression is:

$$rd(\mathbf{x}\_{\boldsymbol{n}}) = \begin{cases} 0 & d(\mathbf{x}\_{\boldsymbol{n}}) > \mathcal{R}\_{2}(\mathbf{x}\_{\boldsymbol{n}})\\ 0.5 - 0.5\pi[d(\mathbf{x}\_{\boldsymbol{n}})/\mathcal{R}(\mathbf{x}\_{\boldsymbol{n}}) - 1.5] & \mathcal{R}\_{1}(\mathbf{x}\_{\boldsymbol{n}}) < d(\mathbf{x}\_{\boldsymbol{n}}) \le \mathcal{R}\_{2}(\mathbf{x}\_{\boldsymbol{n}})\\ 1 & d(\mathbf{x}\_{\boldsymbol{n}}) \le \mathcal{R}\_{1}(\mathbf{x}\_{\boldsymbol{n}}) \end{cases} \tag{11}$$

#### 3.2.3. Energy Consumption Model Based on Dynamic Analysis

Energy efficiency is an important characteristic of path planning algorithms for autonomous systems [58]. The consumed power of a USV is divided into two parts: the static power due to static consumption and the dynamic power involving the thrust power. The main USV power consumption is due to thrusters (80–90%), transferred to the mechanical power. To maximize the engine efficiency [15], only dynamic consumption is considered in this paper, and the mechanical power is modeled as a function of USV speed and environmental conditions. Compared with the approximate path cost of kinematics method, the dynamic analysis proposed by Fossen [50] can provide more accurate information on the energy consumption cost. According to the dynamics model of *Dolphin-I* USV created in Section 2.3, the scalar product of the thrust vector obtained by Formula (7) and velocity vector obtained by MMG formula can get the power dissipated by a given force. Assuming that the instantaneous acceleration during the navigation from the current node to the extended node is constant, the instantaneous power during the navigation is integrated. The heuristic function of the energy consumption cost from the current node to the extended node can be obtained as follows:

$$e(\mathbf{x}\_n) = \int F\_{thrust}(\mathbf{x}\_n) \cdot V\_e(\mathbf{x}\_n) dt \tag{12}$$

#### *3.3. Multi-Objective Optimization Algorithm for Ship Motion Planning Based on HA\**

In the process of expanding nodes, the Hybrid A\* algorithm takes a long time to make a search and requires lots of iterations and calculations. Moreover, this algorithm does not consider the energy consumption and multiple objectives optimization under the influence of navigation environment at the same time. In order to resolve the above-mentioned problem of USV motion planning in a steady navigation environment, a multi-objective Hybrid A\* algorithm in a dynamic environment (MOHA\*) is adopted to update the actual speed in the extended grid after being affected by the environment in real time, and to generate a multi-objective optimal path which is more in line with the actual application and meets the requirement to be the most efficient, safest and the most energy-saving. The complexity of motion planning mainly comes from the following two aspects: the influence of complex environment and the kinematics constraints of the USV. In order to ensure that the estimated cost of the optimal path is close to the actual optimal path cost, these two factors should be considered when designing heuristic functions. With the multi-objective optimization model of USV motion planning in Section 3.2, the framework of the MOHA\* algorithm is presented, as shown in the Figure 9. The core idea of designing the optimal strategy is the definition of the cost function, which can be expressed as:

$$f(\mathbf{x}\_{\mathbb{N}}) = \frac{g(\mathbf{x}\_{\mathbb{N}}) - g\_{\min}(\mathbf{x}\_{\mathbb{N}})}{g\_{\max}(\mathbf{x}\_{\mathbb{N}}) - g\_{\min}(\mathbf{x}\_{\mathbb{N}})} + \frac{h(\mathbf{x}\_{\mathbb{N}}) - h\_{\min}(\mathbf{x}\_{\mathbb{N}})}{h\_{\max}(\mathbf{x}\_{\mathbb{N}}) - h\_{\min}(\mathbf{x}\_{\mathbb{N}})} + \mathrm{asd}(\mathbf{x}\_{\mathbb{N}}) + \beta \varepsilon(\mathbf{x}\_{\mathbb{N}}) \tag{13}$$

where *g*max(*xn*) is the maximum of the actual cost, *g*min(*xn*) is the minimum of the actual cost, *h*max(*xn*) is the maximum of the heuristic cost, *h*min(*xn*) is the minimum of the heuristic cost, *sd*(*xn*) is safety cost heuristic value, *e*(*xn*) is energy cost heuristic value, *α*

and *β* are constants greater than 0, which are used to control the weight of safety cost and energy cost in the total cost, respectively, thus controlling their influence on the final path.

**Figure 9.** Framework of MOHA\* algorithm.

#### **4. Simulation**

In this paper, the MMG motion model and motion mathematical model of the USV are built in the environmental model by integrating ENC data and ocean reanalysis data, and expressing the complex constraints under ocean dynamic interference elements. Thus, the efficient solution and optimization of USV motion planning are realized. The proposed approach is simulated using MATLAB R2020b. All simulations are performed on a PC with Microsoft Windows 10 as OS with Intel i5 2.90 GHz quad core CPU and 8 GB RAM.

The grid map used in the planning is 5000 m × 5000 m in size and 25 m × 25 m in resolution. The approach takes *xs*(1000, 6500, 0) as the starting point, and *x*<sup>g</sup> (3500, 2500, −π/2) as the ending point. The initial surge speed is given as 5 m/s, and the number of extended motion primitives is six. The cost function values are calculated with five different optimization objectives, respectively, including only using original Hybrid A\*, efficiency-optimized Hybrid A\*, safety-optimized Hybrid A\*, energy consumption optimized Hybrid A\* and MOHA\*. The paths that are obtained are shown in Figure 10a–f. Figure 10b represents a partial enlarged view of the traditional HA\* algorithm path in Figure 10a. In the simulation result diagram, blue lines and blue dots represent the moving primitives and nodes in the expansion process, respectively, while gray lines and gray dots represent the moving primitives and nodes in the retreat, respectively, and their density represents the number of expanded nodes. The red curve is the final path obtained by planning. Figure 10c shows the efficiency-optimized Hybrid A\* algorithm results. In Table 3, it can be seen that the number of extended nodes decreases and the search time is shortened overall. The comparison between Figure 10a,d shows that the curvature of some moving primitives of USV changes to a certain extent in the area close to obstacles with the

safety-optimized HA\* algorithm, which leads to a relative increase of expanded nodes. At the expense of a certain calculation time, the risk of generating paths is reduced by 41.3%. According to Figure 10a–f, MOHA\* expands fewer nodes and has a larger distance from obstacles. The navigation state of USV is changed by the environmental interference, and the total energy consumption is reduced by 24.18%. The simulation results for the five different scenarios verify the effectiveness of the multi-objective motion planning model, and the number of extended nodes, risk degree, running time and energy consumption are compared. The results are shown in Table 4.

**Figure 10.** Comparison of motion planning paths. (**a**) Path of original HA\* (Test 1); (**b**) A local enlargement of original HA\*; (**c**) Path of efficiency-optimized HA\* (Test 2); (**d**) Path of safetyoptimized HA\* (Test 3); (**e**) Path of energy-optimized HA\* (Test 4); (**f**) Path of MOHA\* (Test 5).


**Table 4.** Comparison of Motion Planning Test.

The MOHA\* algorithm is used to carry out four simulation experiments on motion planning at different starting points. The results show that the extended search model, based on motion primitives and map downsampling, can effectively reduce the number of extended nodes. The specific combination of motion primitives are shown in Table 5. The final paths of the four groups can avoid obstacles and always keep a safe distance from obstacles. Test 5 and Test 6 show that the USV can make good use of environmental interference and consume less energy when heading in the same direction as the environmental interference force. In the opposite navigation situation, the MOHA\* algorithm can also reduce energy consumption while avoiding obstacles, and can select the optimal path, as shown in Figure 11a–d. In summary, the results show that, under different working conditions, using the MOHA\* algorithm can result in the USV simultaneously satisfying multi-constraint and multi-objective optimization in a complex navigation environment, and find the optimal path.

**Table 5.** Combination of motion primitives for Test 6-Test 9 motion planning paths.


**Figure 11.** *Cont*.

**Figure 11.** Comparison of motion planning paths at different starting and goal points. (**a**–**d**) Test 6–Test 9.

#### **5. Conclusions**

Marine environment information is essential for optimal path planning. In order to quantitatively study the influence of wind and current on the navigation of USVs, this paper establishes a high-resolution marine environment model that is updated every three hours based on the reanalysis data of a numerical prediction model, which provides a highquality marine environment field for USVs. On the basis of the MMG ship operation model, the kinematics and dynamics model of the *Dolphin-I* USV is established, the influence of environmental disturbance force on the motion of the USV is considered, and the simulation is carried out, which can quantitatively calculate the speed and direction of the USV. In this paper, a motion planning MOHA\* algorithm is proposed in order to simultaneously optimize three objectives (efficient, safety and energy) in dynamic marine environments and satisfy multiple constraints. As the navigation environment information changes, the MOHA\* algorithm can adjust the input parameters to complete the path update. The simulation experiments of single-objective optimization and multi-objective optimization show that the MOHA\* algorithm can improve planning time, reduce navigation risk and decrease navigation energy consumption simultaneously, which can fully ensure the efficiency of the USV in performing tasks and improve their independent decision-making ability, thus supporting larger-scale coordinated motion planning and control research for USV clusters.

**Author Contributions:** Conceptualization, M.W. and J.Z.; methodology, M.W.; software, M.W.; validation, M.W. and M.G.; formal analysis, M.W.; investigation, M.W.; resources, A.Z.; data curation, M.W.; writing—original draft preparation, M.W.; writing—review and editing, M.G.; visualization, M.W. and M.G.; supervision, A.Z.; project administration, A.Z.; funding acquisition, A.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by National Key Research and Development Program of China, grant number 2018YFC1503606, 2018YFC1407401 and National Natural Science Foundation of China (41704010) and supported by the Special Fund for Basic Scientific Research Business Expenses of Central Public Welfare Scientific Research Institutes (TKS190302, TKS20210103) and supported by School of Marine Science and Technology of Tianjin University.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data that support the findings of this study are available within the article.

**Acknowledgments:** We are especially grateful to the 3I research team of School of Marine Science and Technology, Tianjin University, for providing us with *Dolphin-I* USV as an experimental platform. **Conflicts of Interest:** The authors declare no conflict of interest.

## **References**

