Next Article in Journal
Fairness-Aware Resource Allocation in Multi-Hop Wireless Powered Communication Networks with User Cooperation
Previous Article in Journal
An Enhanced LoRaWAN Security Protocol for Privacy Preservation in IoT with a Case Study on a Smart Factory-Enabled Parking System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Identification for Vector Propulsion of an Unmanned Surface Vehicle: Three Degrees of Freedom Model and Response Model

School of Marine Electrical Engineering, Dalian Maritime University, Dalian 116026, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(6), 1889; https://doi.org/10.3390/s18061889
Submission received: 18 April 2018 / Revised: 24 May 2018 / Accepted: 6 June 2018 / Published: 8 June 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
This paper presents a complete scheme for research on the three degrees of freedom model and response model of the vector propulsion of an unmanned surface vehicle. The object of this paper is “Lanxin”, an unmanned surface vehicle (7.02 m × 2.6 m), which is equipped with a single vector propulsion device. First, the “Lanxin” unmanned surface vehicle and the related field experiments (turning test and zig-zag test) are introduced and experimental data are collected through various sensors. Then, the thrust of the vector thruster is estimated by the empirical formula method. Third, using the hypothesis and simplification, the three degrees of freedom model and the response model of USV are deduced and established, respectively. Fourth, the parameters of the models (three degrees of freedom model, response model and thruster servo model) are obtained by system identification, and we compare the simulated turning test and zig-zag test with the actual data to verify the accuracy of the identification results. Finally, the biggest advantage of this paper is that it combines theory with practice. Based on identified response model, simulation and practical course keeping experiments are carried out to further verify feasibility and correctness of modeling and identification.

1. Introduction

Ship motion control is an important field that combines navigation science and technology. Its ultimate purpose is to improve the level of ship automation and intelligence, and to ensure the safety, economy and comfort of a ship’s navigation [1,2,3]. Meanwhile, with the continuous development of technology and the extension of the field of ship application, the model of a traditional ship is developing towards miniaturization, speediness and intellectualization, and a new surface carrier has been developed—Unmanned Surface Vehicle (USV) [4,5,6,7]. To allow better maneuverability, the propulsion system of a USV has higher requirements. The vector propulsion system (POD, waterjet and outboard servos, etc.) is a new propulsion device in the field of ship propulsion, which has proven to be more efficient than conventional propeller-rudders [8,9,10].
The mathematical model of a ship is the basis for realizing intelligent control and maneuverability prediction. The history of its development can be traced back to 1946. Davidson and Schiff applied the knowledge of rigid dynamics and fluid mechanics for the first time to propose a maneuvering motion model. After continuous exploration and development, two kinds of mathematical models of ship motion have resulted. One is the Abkowitz model, which is also called the global model. Its essence is to consider the hull, propeller and rudder as a whole and carry out Taylor series expansion at the equilibrium position of the fluid acting on the hull [11,12]; the other is a separate model proposed by a ship manoeuvring mathematical model group (MMG), also known as the MMG model [13,14]. Its essence is to decompose the hydrodynamic forces acting on the bare hull, the thrust of the propeller and the force of the rudder according to the physical meaning, and then to consider the interference between them. In general, the separation modeling theory is more convenient for theoretical analysis than the global modeling theory; therefore, the MMG model is employed to analyze the model structure of vector propulsion USV in this paper. Meanwhile, based on the actual situation, we designed the method for obtaining model parameters. The empirical formula method was employed to obtain the thruster (it is produced by the propeller) model. Based on field experimental data, the model parameters of three degrees of freedom (DOF) model, response model and thruster servo model were obtained by system identification [15,16,17].
In the past few decades, various results have been shown in literature [18,19,20,21]. Fossen et al. put forward a method of modeling and identification for a fully driven model supply ship (mass is 23.8 kg, length is 1.255 m, breadth is 0.29 m), the classic CyberShip II in the Marine Cybernetics Laboratory at Norwegian University of Science and Technology, which is widely used in the study of all kinds of ship motion control [22]. The system identification techniques used in [23] included the extended Kalman filter (EKF) and the constrained least-squares method. In [24], a novel identification scheme for non-linear manoeuvring models based on two steps was proposed. Sonnenburg et al. described planar motion modeling for an USV, including a comparative evaluation of several experimentally identified models over a wide range of speeds and planing conditions [25]. On the premise of obtaining full-scale trial data, a new transformed multi-innovation least squares algorithm was used to identify the model parameters of a four DOF model for a YUKUN ship [26]. In [27], a system-based method was hired to investigate a four DOF ship maneuvering motion in calm water for the ONR tumblehome model. However, most of the literature, including the above, have researched ordinary propeller-rudder propulsion ships. Motivated by the above-mentioned observations, the thruster model, three DOF model and response model, which are obtained using the empirical formula method and system identification method, were established for a vector propulsion USV. First, field experiment data were collected through the related sensors. Then, modeling, identification, and validation work were performed separately. Finally, numerical simulation and real ship experiments for the course keeping of USVs were performed. The main contributions of this paper can be summarized as follows:
(1) According to the force analysis and hypothesis, it is proved that the three degrees of freedom (DOF) model of vector propulsion of USVs is an underactuated system, and its response model still conforms to the classic Norrbin model structure. In addition, this conclusion can be generalized to general vector propulsion ships.
(2) The empirical formula method and system identification method are combined to get the parameters of the three DOF model and response model.
(3) The greatest advantage of this paper is the connection between theory and practice. The actual course keeping field experiment is carried out to further verify the correctness of the modeling and identification results.
The rest of this paper is organized as follows. In Section 2, the Lanxin USV and field experiments are introduced. Section 3 shows the models of thruster thrust and servo. The modeling process of the three DOF model and the response model is shown in Section 4. In Section 5, the identification and validation results are displayed. In Section 6, combining theory with practice, we carry out a real ship course keeping experiment. Section 7 contains the conclusions.

2. Field Experiment of Data Acquisition

In this section, Lanxin USV is first introduced, and then data is collected from field experiments by related sensors.

2.1. Lanxin USV

Lanxin USV is the foundation of the field experiment, which is mainly composed of a power propulsion system, an automatic control system and an information collection system, a communication system, and an image acquisition system. The appearance of the USV is shown in Figure 1, and its specific parameters are displayed in Table 1.

2.2. Vector Propulsion System

The dynamic propulsion system guarantees a continuous output of the USV power and is an important component to ensure the performance of a USV. The dynamic propulsion system is equipped with an electric displacement gasoline engine with a displacement of 5 liters, and its maximum output power can reach to 260 HP. The stern machine is equipped with an 0.46 meter diameter three leaf propeller (For the Lanxin USV, the propeller is its thruster), which can ensure a greater transmission ratio and low speed navigation stability. Vector propulsion system is shown in Figure 2.
Vector propulsion is one of the most promising new technologies in the field of ship propulsion, and has the characteristics of improving the efficiency and maneuverability of the ship. The vector thruster can rotate ±35° around the axis, and can achieve maximum thrust in any direction within the scope of +35° and −35°.

2.3. Platform of Sensor Network

To achieve precise control, the controller first needs to obtain the motion state of the USV and the real-time information of the surrounding environment through a variety of detection sensors. The multi-sensor system of the Lanxin USV includes the following three parts: a GY-86 attitude detection sensor (Simrad, Egersund, Norway), a global position system (GPS) navigation sensor (Simrad, Egersund, Norway) and a yacht equipment network. The attitude detection sensor is communicated with the microcontroller through the inter integrated circuit (IIC) protocol. It mainly provides attitude information about the USV, such as the rolling rate, pitching rate and course angle. The GPS navigation sensor is connected with a microcontroller through a serial port. The communication protocol adopted is the national marine electronics association (NMEA)0183 protocol, which mainly provides latitude and longitude information. The network system of yacht equipment is communicated with the microcontroller through the controller area network (CAN) bus. The communication protocol adopts the NMEA2000 protocol, which mainly provides the state information and information about the surrounding environment of the USV, such as the driving speed, the rotating speed of propeller, the water depth, the power supply voltage and so on. The multi-sensor structure of a USV is shown in Figure 3.
(1) GY-86 Attitude Detection Sensor
The attitude sensor module uses three chips, MPU6050, HMC5883L and MS5611, and four sensors are mounted on board. MPU6050 integrates a microelectro mechanical systems (MEMS) gyroscope and a MEMS accelerometer. The magnetic field intensity sensor used is HMC5883L. The air pressure sensor used is MS5611. They are connected in series and are able to compensate for each other’s defects by using an algorithm.
(2) GPS Navigation Sensor
The GPS navigation sensor adopts the NEO-5Q main chip (U-blox, Zurich, Switzerland). The chip is a multi-functional and independent GPS module, which has the advantages of low cost and small volume. The communication between the GPS navigation module and the motion controller is the NMEA0183 protocol, and the NMEA0183 is now the standard protocol for unified GPS navigation equipment. It is the standard format made by the National Oceanic Association of the United States for marine electronic equipment.
(3) Yacht Equipment Network System
The network system is an important part of the state information collection of unmanned vehicles. The system interconnects the engine, sonar, compass and other craft equipment through the NMEA2000 protocol, and forms an efficient information sharing network. The computer terminal adopts microsoft foundation classes (MFC) visual programming, which can monitor and display various information about the USV in real-time and stores the collected data in the database. The controller adopts the digital modular design based on STM32 as the microcontroller. It mainly completes functions such as data acquisition, information release, and steering control. The controller has a variety of sensor interfaces and integrates the required communication protocols to obtain real-time data from the boat-borne sensor network.

2.4. Field Experiment

The specific contents of the field experiment included the turning test and the zig-zag test. Using a 10° turning test and 10°/10° zig-zag test as examples, their specific processes are as follows:
(1) Turning test: keeping the speed unchanged, rotate the propulsion angle to 10° and wait for the USV to rotate steadily.
(2) Zig-zag test: keeping the speed unchanged, rotate the propulsion angle to 10°, and when the course angle changes to 10°, push the propulsion angle to −10°. After that, until the course angle changes to −10°, rotate the propulsion angle to 10°. Repeat the above steps several times.
Remark 1.
For a traditional propeller-rudder propulsion ship, the direction of navigation can be changed by changing its rudder angle. However, the vector propulsion USV has no rudder, so we call it a propulsion angle.
Remark 2.
The force produced by the propeller when it is rotating clockwise or counterclockwise has different effects on a USV, and the counterclockwise rotation of the propeller leads to reversing, which is not within the scope of this paper.
In this subsection, the contents of the field experiment are described in detail. In order to reduce the influence of external disturbance on the accuracy of the collected data, field experiments are conducted in a relatively calm sea state. The sea state is shown in Figure 4.
The environment and content of the field experiment:
(1)
Experimental sea area: nearby waters (longitude: 121.5548, latitude: 38.8612).
(2)
Sea state: one-level marine conditions. The sea surface was quite calm, and the waves were 0–0.1 m high.
(3)
Weather: the weather was fine and the sea breeze was about a one-level northeasterly wind.
(4)
The driving speed of USV: this was kept at around 10 knots (corresponding to this, the engine speed was about 2800/min).
(5)
The contents of the record were the driving speed, V, course angle, ψ , rotating speed of the propeller, n, surge velocity, u, sway velocity, v, and yaw rate, r.
(6)
Sampling frequency: 0.02 s.
(7)
The specific contents of the field experiment were the turning test (10°, 15°, 25°, 35°) and the zig-zag test (10°/10°).
Remark 3.
The maximum speed of the Lanxin USV is 30 knots. If the driving speed is too fast, the driver’s safety cannot be guaranteed, but slow speeds cannot fully explore the manipulation characteristics of the USV. So, in general, the driving speed of a USV is chosen to be kept around 10 knots.

3. Thruster Thrust Model and Servo Model

3.1. Thruster Thrust Model

The vectorial thrust generated by the thruster (propeller) can be decomposed into two forces: the longitudinal force to keep the USV moving forward and the lateral force moment to change the direction of the USV. In essence, the latter is actually the effect of the rudder. To some extent, this vector propulsion device simplifies the traditional propeller rudder, especially for small ships represented by USVs, which greatly improves the propulsion efficiency. Meanwhile, the vectorial propulsion system has a more concise mechanical structure, which is more suitable for smaller USVs.
According to [28], the thrust mathematical model of the vector thruster can be expressed as
T = ( 1 t P ) ρ n 2 D P 4 K T ( J p ) ,
where t P represents the coefficient of the thrust reduction, ρ indicates the density of the sea water, n is the rotating speed of the propeller, D p is the diameter of the propeller, K T ( J P ) is a propeller thrust coefficient, and J P is the advance coefficient.

3.1.1. Thrust Reduction Factor ( t P )

The increase in hull resistance caused by propellers under working conditions is called an increase of resistance, which is represented by Δ R . If the propulsive force generated by the propeller under working conditions is P, part of the thrust will be used to overcome the resistance R of the forward hull, and another part will be used to overcome the resistance increase, Δ R . Thus, it can be seen that only part of the ( P Δ R ) is used to push the USV forward, so it is called the effective thrust, P e . Customarily, Δ R is called the thrust reduction, which is expressed as Δ P . In other words, P e = R and Δ P = Δ R . Then, the coefficient of thrust reduction, t P , can be obtained by Δ P and P. That is
t P = Δ P P = P P e P = P R P .
The value of t P is determined by a variety of factors, including the shape of the USV, the size of the propeller, the loading load and so on. In this paper, Hollo’s formula for a single boat is used to estimate the t P :
t P = 0.001979 L ( B B C p 1 ) + 1.0585 C 10 0.000524 0.1418 D P 2 B d + 0.0015 C s t e r n ,
where B indicates the maximum transverse width of USV hull, L is the maximum longitudinal length of the USV hull, d represents the average draught depth of the USV. C p 1 is defined as C p 1 = 1.45 C P 0.315 0.225 L c b , where C P is the prismatic coefficient and L c b 0.5 L . C 10 is determined by the ratio of the longitudinal and transverse widths of the USV, and C s t e r n is a coefficient associated with the shape of the tail of the USV. Their specific definitions are shown in (4) and (5).
L / B > 5.2 C 10 = B / L L / B 5.2 C 10 = 0.25 0.003328402 B / L 0.134615385
V section C s t e r n = 10 Conventional C s t e r n = 0 U section C s t e r n = + 10
According to the actual situation of the Lanxin USV and the above formula, its thrust reduction factor is t P 0.05 .

3.1.2. Thrust Coefficient ( K T )

The propeller of the USV is an ordinary three-bladed paddle, and according to [25], the thrust coefficient, K T , of a three-bladed series paddle can be expressed as
K T = i = 0 n 1 j = 0 n 2 A i j P D i J P j ,
where the advance coefficient is J P = 1 w p u n D P .
w p is the propeller flow coefficient. According to the Bamier formula, the expression for w p is
w p = 0.165 C b x 3 D P Δ w ,
where ∇ represents the drainage volume of the USV, and x is an exponent. When x = 1 , it is applicable to the middle line propeller; when x = 2 , it is applicable to the side propeller. Δ w represents the modified value of the flow coefficient, which can be expressed as
Δ w = 0.1 ( F n 0.2 ) , F n > 0.2 0 , F n 0.2 ,
where F n is the Furude coefficient. F n = V / g L , and g is the acceleration of gravity. Based on the conditions of the field experiment in Section 2.2 and the above formula, one can deduce that w p 0.04 .
Then, using Lagrange interpolation, the disk surface ratio is also taken into account as an influencing variable. Based on this, (6) can be rewritten as
K T = K T 0.5 ( θ 0.8 ) ( θ 1.1 ) 0.18 + K T 0.8 ( θ 0.5 ) ( θ 1.1 ) 0.09 + K T 1.1 ( θ 0.5 ) ( θ 0.8 ) 0.18 ,
where K T 0.5 , K T 0.8 , and K T 1.1 are thrust factors with disk surface ratios of 0.5, 0.8, and 1.1, respectively.
In the first quadrant, the propeller’s regression formula is expressed as
K T = A 00 + A 01 J P + A 02 J P 2 + A 10 P D + A 20 P D 2 + A 30 P D 3 + A 11 J P P D + A 12 J P 2 P D + A 21 J P P D 2 + A 22 J P 2 P D 2 .
Meanwhile, the coefficients of regression formulas for K T 0.5 , K T 0.8 and K T 1.1 are shown in Table 2.

3.2. Servo Model

The traditional rudder servo model is considered to be a first order inertial link [28,29,30]. That is,
δ ˙ = 1 T d δ + 1 T d δ d ,
where T d is a time constant, δ d is the target propulsion angle, δ is actual propulsion angle, and δ 35 ° .
The servo system of the USV has a higher sensitivity and faster activity. After analyzing the real ship’s experimental data, we cannot think of the servo model as a first-order inertia link. It should be regarded as a two-order under-damped system.
δ ¨ + 2 ζ ω n δ ˙ + ω n 2 δ = K n ω n 2 δ d
where ω n is the natural frequency, ζ is the damping ratio, and K n is the magnification factor.

4. Modeling

4.1. Three DOF Model

Control is the core of many problems, and the model is the basis of control. The design effect of the USV motion controller depends not only on the selected control algorithm, but also on the accuracy of the mathematical model. The complexity of the USV model is mainly manifested in the viscous hydrodynamic force of the hull, the control of the input force/torque, the high nonlinearity and coupling of the external disturbance. Therefore, the simplification of any structure in the model inevitably ignores some important characteristics of the system. The actual movement of the USV is very complex, with six DOF in general, including the surge velocity, sway velocity, heave velocity, yaw rate, rolling rate, and pitching rate. In previous research on the ship model, to simplify the research difficulty, the heave velocity, rolling rate and pitching rate have often been ignored. In other words, only consider the surge velocity u, the sway velocity v and the yaw rate r have been considered. The relationship between them is shown in Figure 5.
Before that, we first need to define some variables. x g is the distance from the center of the USV to the center of gravity, X u ˙ , Y r ˙ , Y v ˙ , N r ˙ , I z ˙ , X u , Y v , and Y r , N v , N r are the corresponding hydrodynamic coefficients, respectively.
The three degree of freedom dynamic equation established by the Lagrange’s mechanics theory can be expressed as
M υ ˙ + C ( υ ) υ + D ( υ ) υ = τ ,
where M = m X u ˙ 0 0 0 m Y v ˙ m x g Y r ˙ 0 m x g Y r ˙ I z N r ˙ , υ = [ u , v , r ] T , C ( υ ) = 0 0 ( m Y v ˙ ) v ( m x g Y r ˙ ) r 0 0 ( m X u ˙ ) u ( m Y v ˙ ) v + ( m x g Y r ˙ ) r ( m X u ˙ ) u 0 , D ( υ ) = X u 0 0 0 Y v Y r 0 N v N r and τ = [ τ u , τ v , τ r ] T .
M is called the inertia matrix, and C ( υ ) is the Coriolis/centripetal force matrix. D ( υ ) is the hydrodynamic damping matrix. τ u is the longitudinal thrust, τ v is the lateral thrust, and τ r is the yaw moment. For the Lanxin USV, the vector thrust distribution direction is shown in Figure 6.
When propulsion angle is δ , the thrust distribution values in each direction are as shown in (14).
τ u = T cos δ τ v = T sin δ τ r = x l T sin δ ,
where x l is the distance from the center of rotation of the longitudinal arm of the USV to the axial point of the thruster. The effective attack angle of the thruster, α R , is a small value with the unit “rad”. Thus, we have sin α R α R = δ , where the propulsion angle is δ [ 0.5236 r a d , 0.5236 r a d ] . If x l > 1 ; we can think of τ v as being equal to zero. At the same time, m 11 = m X u ˙ , m 22 = m Y v ˙ , m 23 = m x g Y r ˙ , m 32 = m x g N v ˙ , m 33 = I z N r ˙ , d 11 = X u , d 22 = Y v , d 23 = Y r , d 32 = N v , and d 33 = N r are defined. Then, (13) can be changed into
m 11 u ˙ m 22 v ˙ + m 23 r ˙ m 32 v ˙ + m 33 r ˙ + d 11 u m 22 v r m 23 r 2 d 22 v d 23 r + m 11 u r m 22 v u + m 23 r u + d 32 v m 11 u v + d 33 r = τ u 0 τ r
One can get (16) by simplifying (15):
u ˙ = d 11 u m 11 m 22 v r m 11 m 23 r 2 m 11 + τ u m 11 v ˙ = ( m 23 m 22 2 m 11 m 22 m 23 m 11 m 23 m 11 m 22 m 33 m 11 m 23 m 32 ) u v m 22 ( m 22 m 33 m 23 m 32 ) + m 23 2 r v ( m 22 m 33 m 23 m 32 ) + ( m 22 m 23 d 32 + m 23 m 32 d 23 + m 23 m 32 d 22 m 22 m 33 d 22 ) v m 22 ( m 22 m 33 m 23 m 32 ) + ( m 33 d 33 + m 33 d 23 ) r ( m 22 m 33 m 23 m 32 ) m 23 τ r ( m 22 m 33 m 23 m 32 ) r ˙ = ( m 11 + m 11 m 22 m 22 2 ) u v m 22 m 33 m 23 m 32 m 22 m 23 r v m 22 m 33 m 23 m 32 ( m 22 d 32 + m 32 d 23 ) v m 22 m 33 m 23 m 32 ( m 22 d 33 + m 32 d 23 ) r m 22 m 33 m 23 m 32 + m 22 τ r m 22 m 33 m 23 m 32 .
Assumption 1.
It is assumed that a USV is symmetrical and the barycenter of the USV coincides with the center of the body-fixed frame. That is to say, x g = 0 , Y r ˙ = 0 , N v ˙ = 0 , Y r = 0 , and N v = 0 .
Based on Assumption 1, (16) can be reduced to
u ˙ = m 22 m 11 v r d 11 m 11 u + 1 m 11 τ u v ˙ = m 11 m 22 u r d 22 m 22 v r ˙ = m 11 m 22 m 33 u v d 33 m 33 r + 1 m 33 τ r .
It can be seen that for the vector propulsion of a USV, the three DOF model is an underactuated system. In addition, the results of this theoretical study can also be extended to general vector propulsion ships.

4.2. Response Model

In Section 4.1, we concluded that the three DOF model of vector propulsion of a USV is an underactuated system. However, due to the many limitations of underactuated systems, the ship’s underactuation model is used in theoretical studies, in most cases [31,32]. In practical engineering applications, the response model is mainly used for course control and path following.
According to Assumption 1, the mathematical model of planar motion with three DOF can also be expressed as
( m + m x ) u ˙ ( m + m y ) v r = X H + X P ( m + m y ) v ˙ + ( m + m x ) u r = Y H + Y P ( I z z + J z z ) r ˙ = N H + N P ,
where X H , Y H , and N H are the hydrodynamic forces and moments acting on bare hulls. The Taylor series expansion can be expressed as
X H = X ( u ) + X H v v v 2 + X H v r v r + X H r r r 2 Y H = Y H v v + Y H r r + Y N L N H = N H v v + N H r r + N N L ,
where I z z is the moment of inertia.
Remark 4.
The added mass and the added moment of inertia are essentially the same as the acceleration hydrodynamic derivatives. Their mutual correspondence is that m x X u ˙ , m y Y v ˙ , J z z N r ˙ , τ u X p , τ v Y p , and τ r N p .
For convenience of research, X H , Y H and N H needed to be linearized. This so-called linearization means that the USV receives less external disturbance, and its motion is always near the initial equilibrium state. At the same time, all kinds of hydrodynamic terms acting on USV are dominated by linear terms, and the orders of magnitude above the second order are negligible.
Generally, the uniform rectilinear motion of a USV is taken as the initial equilibrium state, and it is assumed that u = u 0 , v = v 0 = 0 , r = r 0 = 0 and δ = δ 0 = 0 . u 0 is the initial longitudinal velocity of a USV. When a USV is subjected to external interference, the variation in its motion state is Δ u , Δ v = v , Δ r = r , and Δ δ = δ , respectively. Then, the motion state of the USV is changed to u = u 0 + Δ u , v = v 0 + Δ v , r = r 0 + Δ r and δ = δ 0 + Δ δ .
To preserve the first order small quantities, Δ u , v, r and δ , and ignoring the two order and high order small quantities on the machine, (19) can be simplified to
X H = X ( u 0 + Δ u ) Y H = Y H v v + Y H r r N H = N H v + N H r r ,
where X ( u 0 + Δ u ) is the direct resistance of USV. It can be further expressed as
X ( u 0 + Δ u ) = 1 2 ρ S C t ( u 0 + Δ u ) 2 ,
where S is the wet area, ρ is the water density, and C t is the total drag coefficient.
X ( u 0 + Δ u ) = 1 2 ρ S [ C t 0 + ( C t Δ u ) u 0 Δ u ] ( u 0 + Δ u ) 2 .
When the speed is u 0 , its total resistance coefficient is C t 0 . The Δ u of (22) is linearized as
X H = 1 2 ρ S C t 0 u 0 2 1 2 ρ S [ 2 C t 0 u 0 + ( C t Δ u ) u 0 u 0 2 ] Δ u .
To define X 0 = 1 2 ρ S C t 0 u 0 2 and X H u = 1 2 ρ S [ 2 C t 0 u 0 + ( C t Δ u ) u 0 u 0 2 ] , X 0 is used to represent the straight line resistance of the unmanned vehicle in the initial state. Then X H = X 0 + X H u Δ u . Thus, (20) can be expressed as
X H = X 0 + X H u Δ u Y H = Y H v v + Y H r r N H = N H v + N H r r .
If δ is small, then sin δ = δ and cos δ = 1 . One can obtain
X P = T Y P = T δ N P = x l T δ .
A very important problem to note is that, in the initial state, the resistance of a USV is balanced with the thrust of the thruster (propeller), which means that X 0 + X P = 0 . Based on the above conditions, (18) can be changed into
( m + m x ) Δ u ˙ = X u Δ u ( m + m y ) v ˙ + ( m + m x ) u 0 r = Y H v v + Y Hr r + X p δ ( I z z + J z z ) r ˙ = N H v v + N H r r + x l X p δ .
It is visible from the upper form that the first equation is decoupled from the second and third equations in linear motion. In other words, the longitudinal motion and the rotation motion of the USV can be considered separately.
Assumption 2.
The influence of longitudinal velocity change is not taken into account, and the external disturbance of the USV received is very weak.
In practical applications, engineers consider the r more than the v. N H v = N v , N H r = N r , and x l X p = N δ are defined, and the third equation of (26) can be reexpressed as
( I z z + J z z ) r ˙ = N v v + N r r + N δ δ .
In order to simplify the problem, it is assumed that the initial state is uniform motion, and all the motion variables have an initial value of zero. Then, v ( 0 ) = 0 , v ˙ ( 0 ) = 0 , r ( 0 ) = 0 , r ˙ ( 0 ) = 0 , δ ( 0 ) = 0 , and δ ˙ ( 0 ) = 0 . After the Laplace transformation, Equation (27) can be changed into
( I z z + J z z ) s r ( s ) = N v v ( s ) + N r r ( s ) + N δ δ ( s ) ,
where v ( s ) = L [ v ( t ) ] , r ( s ) = L [ r ( t ) ] , and δ ( s ) = L [ δ ( t ) ] . The transfer function between the propulsion angle, δ , and the yaw rate, r, can be obtained.
H ( s ) = r ( s ) δ ( s ) = K ( 1 + T 3 s ) ( 1 + T 1 s ) ( 1 + T 2 s ) ,
where T 1 , T 2 , T 3 and K are corresponding parameters. T 1 T 2 = ( L U ) ( m + m y ) ( I z z + J z z ) C , T 1 + T 2 = ( L U ) ( m + m y ) N r ( I z z + J z z ) Y v C , K = ( U L ) N v Y δ N δ Y v C , K T 3 = ( m + m y ) N δ C , C = Y v N r [ Y r ( m + m x ) ] N v , where m is the normalization of m. Similarly, m x , m y and I z z , etc. use superscripts to represent their normalization. Using the Laplace inverse transform, Equation (29) can be transformed into a linear response equation in the time domain:
T 1 T 2 r ¨ + ( T 1 + T 2 ) r ˙ + r = K ( δ + T 3 δ ˙ ) .
It is also important to note that if a USV has the characteristic of course stability, C > 0 ; if the USV does not have the characteristic of course stability, C < 0 ; C 0 is called the critical stability. Therefore, C is the number of stability criteria. The nonlinear change of C was proposed by Nomoto, which is expressed as
C = C 0 + ϖ r 2 ,
where C 0 represents the value of C when r = 0 . ϖ is a newly introduced variable.
Substituting (31) into (30), we obtain
( L V ) 2 ( m + m y ) ( I z z + J z z ) r ¨ + ( L V ) ( m + m y ) N r ( I z z + J z z ) Y v r ˙ + ( C 0 + ϖ r 2 ) r = ( L V ) ( N v Y δ N δ Y v ) δ + ( m + m y ) N δ δ ˙ .
The two sides of (32) are divided by C 0 , and one can obtain
( L V ) 2 ( m + m y ) ( I z z + J z z ) C 0 r ¨ + ( L V ) ( m + m y ) N r ( I z z + J z z ) Y v C 0 r ˙ + C 0 + ϖ r 2 C 0 r = ( L V ) N v Y δ N δ Y v C 0 δ + ( m + m y ) N δ C 0 δ ˙ .
Meanwhile, α = n C 0 is defined, and then (33) can be simplified as
T 1 T 2 r ¨ + ( T 1 + T 2 ) r ˙ + r + α r 3 = K δ + K T 3 δ ˙ .
The nonlinear influence is embodied by α r 3 , and α is the newly introduced constant. In practical applications, (34) is often simplified to a first order form. That is
T r ˙ + r + α r 3 = K δ .
This is Norrbin nonlinear response model used in the field of ship motion control [33,34]. Through the study of this paper, we know that the response model of the class of vector propulsion ships still conforms to the classical Norrbin model. Equation (35) can be reexpressed as T r ˙ + r = K δ plus α r 3 . T r ˙ + r = K δ is the classic Nomoto model [35,36].
When the USV is conducting the turning test, r ˙ = 0 , and δ ˙ = 0 . (35) can be simplified as
r + α r 3 = K δ .
This means that in the case of a known K, the nonlinear term coefficient, α , can be fitted through a series of turning tests.
Remark 5.
During the turning test, ideally, an propulsion angle will correspond to a constant value of r. However, under the influence of various conditions, such as external disturbance, the value of r fluctuates. So when fitting the α value, we need to calculate an average r value.

5. Identification and Verification

5.1. Identification

5.1.1. Three DOF Model

(1)
Data: because the data of the zig-zag model is more able to exert the manoeuvre characteristics of the USV, it is used to identify the three DOF underactuated model. Of course, we mainly used u, v, r, n and δ in the zig-zag test.
(2)
τ u and τ r : based on the modeling of the thruster thrust in the Section 2, the real-time τ u and τ r were calculated based on the rotating speed of the propeller and the propulsion angle.
(3)
Based on the sampling time of 0.02 seconds, u ˙ , v ˙ and r ˙ were calculated.
(4)
Recursive least squares method was used to identify the parameters of underactuated model.
The results of the identification were as follows:
u ˙ = 1.065245 v r 0.3197455 u + 0.000377 τ u v ˙ = 0.938753 u r 3.5 v r ˙ = 0.041187 u v 5.407752 r + 0.000238 τ r .

5.1.2. Response Model

Remark 6.
In order to improve the identification accuracy, we first identified the Nomoto model. Then, the nonlinear parameter, α, was fitted through a series of turning test data.
(1) In the zig-zag test, r and δ were used to identify the Nomoto model. Then, the propulsion angles, δ , of four sets of turning tests and the corresponding average, r, were used to fit α .
(2) The recursive least squares method was used to identify the Nomoto model, and the method of data fitting was used to obtain α .
The identification result from the Nomoto model was
0.332 r ˙ + r = 0.707 δ .
That is to say, K = 0.332 and T = 0.707 . Based on this, Equation (35) was rewritten as
0.332 r ˙ + r + α r 3 = 0.707 δ .
The result of parameter fitting was α = 1.102 . The final identification result of the Norrbin model was
0.332 r ˙ + r + 1.102 r 3 = 0.707 δ .

5.1.3. Servo Model

(1) Data: the target propulsion angle δ d (input) and the actual propulsion angle δ (output) are derived from the zig-zag test data.
(2) The recursive least squares method was used to identify the servo model.
The identification result of servo model was
δ δ d = 1.95 s 2 + 2.38125 s + 1.95 e t .
Remark 7.
The upper control computer output a control command every 0.02 seconds through the timer. The servo system took about 1 second from receiving the instruction to actually start to perform rotation. So a 1 s delay was added to the servo model.

5.2. Verification

In this section, the identified models were simulated by the turning test and the zig-zag test, and then the results were compared with the actual data, to verify the feasibility and correctness of the modeling and identification results.

5.2.1. Three DOF Model

First, we carried out the zig-zag test simulation experiment for model (37), and the results of the comparison are shown in Figure 7.
Figure 7a depicts the actual trajectory and the simulation trajectory. It can be seen from Figure 7a that the maximum lateral error between the actual data and simulated data is about 1.5 m, which accounts for 12.5% of the total lateral distance. The maximum longitudinal error is about 10 m, accounting for 4.5% of the total longitudinal distance. Due to the existence of external disturbance, the USV will have some drift, so the proportion of lateral error is greater than that of longitudinal, which is also reasonable. Figure 7b shows the comparison of u, v and r in the actual and simulated data. It can be clearly seen from the picture that the results of u and r basically coincided, and the simulation trend, v, was also the same. It is obvious from Figure 7c that the maximum difference between the actual and simulated course angles is about 3°. In addition, in the later stage of simulation, compared with the actual data, the phase difference between the simulated course and propulsion angle is caused by external disturbance.
Remark 8.
First of all, external disturbances include the actual wind, waves, and currents, and the currents generated by the USV’s own motion. To be more precise, in a stable sea state, the adverse effects of the currents generated by USV’s own motion may be greater than weather. Secondly, the sensors’ own factors can also cause measurement noise. Due to the above two points, the measured values of the yaw rate had a jagged appearance.
Remark 9.
Although the experimental environment is selected in a relatively calm sea state, the disturbance of wind, waves and currents still exists. Compared with large merchant ships, the USV has a smaller volume and is more sensitive to the external environment, especially lateral disturbance. Therefore, in the result of model validation, the error of v is larger than that of u and r.
Remark 10.
In the process of model identification, a certain modeling error is allowed. The main reasons are as follows:
(1) 
In terms of modeling theory, the derivation of the model is based on various assumptions and simplification. That is to say, it is difficult to fully reflect the characteristics of a USV with the mathematical model.
(2) 
During the voyage, due to the influences of sea condition, operation and various factors, the structure and parameters of a USV will change. In addition, when designing various USV controllers, designers take the uncertainty of the model parameters or structure into account [37].
(3) 
For the real ship, even though the related field experiments are carried out in a relatively calm sea area, the external interference is inevitable.
The comparison between simulation turning test and actual turning test is shown in Figure 8.
As can be seen from Figure 8a, the difference between the simulated radius of rotation and the actual radius was approximately 1.5 m, which accounts for 3.8% of the actual radius. The comparison results of u, v and r are provided in Figure 8b, and the differences between the simulation value and actual value were within a reasonable range. Figure 8c shows the response curves of the propulsion angle. The target propulsion angle was 10°, and the actual target propulsion angle was about 11°, which is due to external disturbance.

5.2.2. Response Model

In this section, we verify the correctness of the identification of the response model through the simulation and actual comparison. The comparison results of the zig-zag test are provided in Figure 9, and the comparison results of the turning test are plotted in Figure 10.
Figure 9 shows that as with the verification of the three DOF model, there was a certain phase difference in the course. Meanwhile, the maximum course error was approximately 3°, accounting for 7.5% of the total course angle. Figure 10 displays that the difference between the simulation trajectory and the actual trajectory was very small. The above two comparison results are sufficient to prove the correctness of the identification results of the final response model.
Remark 11.
The validation of the thruster thrust model and servo model were completed in the verification of three DOF model, so there was no need to verify them separately.

6. Course Keeping Field Experiment

On the basis of the obtained thruster servo model and response model, the numerical simulation and field experiment of course keeping are given to further verify the results of the theoretical research.

6.1. Numerical Simulation

In this paper, the proportional-derivative (PD) course keeping controller was used for the simulation and field experiments. At the same time, the PD control was compared with the proportional-integral-derivative (PID) control to verify the role of the integral term in the course control. The expression of the PID control is as follows [38]:
δ = K p ( ψ ψ d ) K i ( ψ ψ d ) K d ( ψ ψ d ) ,
where K p , K i and K d are the three positive control parameters, and ψ d is the target course. The parameters of the PID control are K p = 0.7 , K i = 0.00095 and K d = 1 . The parameters of the PD control are K p = 0.7 and K d = 1 . The gains of the PD controller were manually adjusted by the manipulator. The specific steps were as follows: (1) the control gain, K p , was gradually adjusted until the course appeared to have equal amplitude oscillation; and (2) K d was slowly adjusted from zero to optimize the control effect (this is also a Ziegler–Nichols tuning method. Of course, this is done manually). Meanwhile, the integral of time-weighted absolute error (ITAE) index was used to quantify the control accuracy of the two algorithms:
I T A E = 0 t t ψ ψ d d t .
Remark 12.
In actual ship motion control, PD control is often used instead of PID control. The reason for this is that the marine environment is often very bad, and external disturbance is inevitable. Persistent disturbance often causes saturation of the integral, so the integral term should be used with caution.
The white noise is used to drive the transfer function, 0.42 s s 2 + 036 s + 0.37 , to describe the external disturbance. The disturbance curve is depicted in Figure 11.
The initial course is 0°, and the target course is 100°. The simulation results are plotted in Figure 12, and the ITAE comparison results are shown in Table 3.
Figure 13 shows that the control effects of PID and PD were almost the same, and both were able to maintain the course of a USV near the target course. In addition to that, the ITAE indexes of PID and PD were 37.98 and 38.6 respectively. The introduction of an integral item in the same case of K p and K d played a positive role, but the value of this integral term was very small. Therefore, the integral item can be introduced in actual course control, but there is need to be cautious.

6.2. Field Experiment

The final purpose of modeling and identification is to carry out practical applications. In this subsection, we conduct course keeping in a field experiment to compare with the numerical simulation.
(1)
Experimental sea area: nearby waters (longitude:121.5548, latitude: 38.8612)
(2)
Sea state: three-level to four-level marine conditions.
(3)
Weather: the sea breeze was about three-levels north wind.
(4)
The driving speed of USV was kept at around 10 knots.
The initial course of the USV was 0°, and the target course was 100°. The results of the field experiment are shown in Figure 13.
The actual performances were very similar to the simulation results of the numerical simulation. The course angle was maintained near the target value, and the propulsion angle was constantly changing to resist the external disturbance. In addition, due to the course keeping, the field experiment was carried out in three- to four-level marine conditions, the initial propulsion angle was not 0°, in order to maintain stability. So, there appears to be a jump at the beginning of propulsion angle curve.

7. Conclusions

In this paper, a complete set of schemes has been proposed for vector propulsion of a USV, from data acquisition, model establishment, parameter identification, result verification, and finally, a real ship field experiment. Linking theory with practice is the greatest advantage of this paper. First, based on the hypothesis and simplification, the thrust model, thruster servo model, three DOF model and response model re established respectively. It was proven that vector propulsion of a USV belongs is an underactuated system, and its response model is in accordance with the Norrbin model structure. Then, the parameters of three DOF model, response model and thruster servo model were identified and verified. Finally, based on the thruster servo model and the response model, numerical simulation and real ship field experiment for course keeping were carried out. In future research, the trajectory tracking controller will be designed.

Author Contributions

The work presented here was performed in collaboration among all authors. D.M designed, analyzed, and wrote the paper. G.W guided the full text. Y.F conceived idea. X.S and B.Q analyzed the data. All authors have contributed to and approved the manuscript.

Funding

This research was funded by the Nature Science Foundation of China, grant number 51609033, the Nature Science Foundation of Liaoning Province of China, grant number 2015020022, and the Fundamental Research Funds for the Central Universities, grant numbers 3132014321, 3132016312 and 3132017133.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
USVunmanned surface vehicle
DOFdegree of freedom
MMGmanoeuvring mathematical model group
EKFextended Kalman filter
GPSglobal position system
IICinter integrated circuit
NMEAnational marine electronics association
CANcontroller area network
MEMSmicroelectro mechanical systems
MFCmicrosoft foundation classes
PIDproportional integral derivative
PDproportional derivative
ITAEthe integral of time-weighted absolute error

References

  1. Zheng, J.; Meng, F.; Li, Y. Design and Experimental Testing of a Free-Running Ship Motion Control Platform. IEEE Access 2018, 6, 4690–4696. [Google Scholar] [CrossRef]
  2. Hou, X.R.; Zou, Z.J.; Liu, C. Nonparametric Identification of Nonlinear Ship Roll Motion by Using the Motion Response in Irregular Waves. Appl. Ocean Res. 2018, 73, 88–99. [Google Scholar] [CrossRef]
  3. Jiao, J.; Ren, H.; Adenya, C.A.; Chen, C. Development of a Shipboard Remote Control and Telemetry Experimental System for Large-Scale Model’s Motions and Loads Measurement in Realistic Sea Waves. Sensors 2017, 17, 2485. [Google Scholar] [CrossRef] [PubMed]
  4. Mu, D.; Wang, G.; Fan, Y.; Qiu, B.; Sun, X. Adaptive Trajectory Tracking Control for Underactuated Unmanned Surface Vehicle Subject to Unknown Dynamics and Time-Varing Disturbances. Appl. Sci. 2018, 8, 547. [Google Scholar] [CrossRef]
  5. Mou, X.; Wang, H. Wide-Baseline Stereo-Based Obstacle Mapping for Unmanned Surface Vehicles. Sensors 2018, 18, 1085. [Google Scholar] [CrossRef] [PubMed]
  6. Wei, Y.; Zhang, Y. Effective Waterline Detection of Unmanned Surface Vehicles Based on Optical Images. Sensors 2016, 16, 1590. [Google Scholar] [CrossRef] [PubMed]
  7. Peng, Z.; Wang, D.; Wang, W.; Liu, L. Neural Adaptive Steering of an Unmanned Surface Vehicle with Measurement Noises. Neurocomputing 2016, 186, 228–234. [Google Scholar] [CrossRef]
  8. Yao, W.; Zhang, J.; Liu, Y.; Zhou, M.; Sun, M.; Zhang, G. Improved Vector Control for Marine Podded Propulsion Control System Based on Wavelet Analysis. J. Coast. Res. 2015, 73, 54–58. [Google Scholar] [CrossRef]
  9. Trancossi, M. Design of Acheon Thrust and Vector Propulsion System. In Proceedings of the SAE 2015 AeroTech Congress Exhibition, Seattle, WA, USA, 21–23 September 2015. [Google Scholar]
  10. Chen, C.W.; Cai, Q.W.; Yu, L.; Jiang, Y.; Kang, D.D.; Leng, J.X. Numerical Study of Hydrodynamic Performance of a High Maneuverability Ocean Engineering Ship with POD Propeller Design. In Proceedings of the OCEANS, Anchorage, AK, USA, 18–22 September 2017; pp. 1–5. [Google Scholar]
  11. Abkowitz, M.A. Lectures on Ship Hydrodynamics–Steering and Manoeuvrability; Hydro-and Aerodynamic Laboratory: Lyngby, Denmark, 1964. [Google Scholar]
  12. Abkowitz, M.A. Measurement of Hydrodynamic Characteristics from Ship Maneuvering Trials by System Identification. SNAME Trans. 1980, 88, 283–318. [Google Scholar]
  13. Asmara, I.P.S.; Kobayashi, E.; Pitana, T. Simulation of Collision Avoidance in Lamong Bay Port by Considering the PAW of Target Ship Using MMG Model and AIS Data. Mar. Eng. Front. 2014, 2, 31–38. [Google Scholar]
  14. Yasukawa, H.; Yoshimura, Y. Introduction of MMG Standard Method for Ship Maneuvering Predictions. J. Mar. Sci. Technol. 2015, 20, 37–52. [Google Scholar] [CrossRef]
  15. Jia, X.L.; Yang, Y.S. Mathematical Model of Ship Motion-Mechanism Modeling and Identification Modeling; Dalian Maritime University Press: Dalian, China, 1999. [Google Scholar]
  16. Mu, D.; Wang, G.; Fan, Y. Design of Adaptive Neural Tracking Controller for Pod Propulsion Unmanned Vessel Subject to Unknown Dynamics. J. Electr. Eng. Technol. 2017, 12, 2365–2377. [Google Scholar]
  17. Sun, X.; Wang, G.; Fan, Y.; Mu, D.; Qiu, B. An Automatic Navigation System for Unmanned Surface Vehicles in Realistic Sea Environments. Appl. Sci. 2018, 8, 193. [Google Scholar]
  18. Mu, D.; Wang, G.; Fan, Y.; Zhao, Y. Modeling and Identification of Podded Propulsion Unmanned Surface Vehicle and Its Course Control Research. Math. Probl. Eng. 2017, 2, 1–13. [Google Scholar] [CrossRef]
  19. Wirtensohn, S.; Wenzl, H.; Tietz, T.; Reuter, J. Parameter Identification and Validation Analysis for a Small USV. In Proceedings of the 20th International Conference on Methods and MODELS in Automation and Robotics, Miedzyzdroje, Poland, 24–27 August 2015. [Google Scholar]
  20. Sathyaseelan, D.; Hariharan, G.; Kannan, K. Parameter Identification for Nonlinear Damping Coefficient from Large-Amplitude Ship Roll Motion Using Wavelets. Beni-Suef Univ. J. Basic Appl. Sci. 2017, 6, 138–144. [Google Scholar] [CrossRef]
  21. Luo, W.; Li, X. Measures to Diminish the Parameter Drift in the Modeling of Ship Manoeuvring Using System Identification. Appl. Ocean Res. 2017, 67, 9–20. [Google Scholar] [CrossRef]
  22. Skjetne, R.; Smogeli, Y.N.; Fossen, T.I. A Nonlinear Ship Manoeuvering Model: Identification and Adaptive Control with Experiments for a Model Ship. Model. Identif. Control 2004, 25, 3–27. [Google Scholar] [CrossRef]
  23. Araki, M.; Hosseini, H.S.; Sanada, Y.; Tanimoto, K.; Umeda, N. Estimating Maneuvering Coefficients Using System Identification Methods with Experimental, System-Based, and CFD Free-Running Trial Data. Ocean Eng. 2012, 51, 63–84. [Google Scholar] [CrossRef] [Green Version]
  24. Herrero, E.R.; González, F.J.V. Two-step Identification of Non-Linear Manoeuvring Models of Marine Vessels. Ocean Eng. 2012, 53, 72–82. [Google Scholar] [CrossRef]
  25. Sonnenburg, C.R.; Woolsey, C.A. Modeling, Identification, and Control of an Unmanned Surface Vehicle. J. Field Robot. 2013, 30, 371–398. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, G.; Zhang, X.; Pang, H. Multi-Innovation Auto-Constructed Least Squares Identification for 4 DOF Ship Manoeuvring Modelling with Full-Scale Trial Data. Isa Trans. 2015, 58, 186–195. [Google Scholar] [CrossRef] [PubMed]
  27. Guo, H.P.; Zou, Z.J. System-Based Investigation on 4-DOF Ship Maneuvering with Hydrodynamic Derivatives Determined by RANS Simulation of Captive Model Tests. Appl. Ocean Res. 2017, 68, 11–25. [Google Scholar] [CrossRef]
  28. Du, J.L.; Guo, C.; Zhang, X.K. Simulation studies on adaptive nonlinear control of ship motion course. J. Syst. Simul. 2005, 17, 1445–1448. [Google Scholar]
  29. Meziou, M.T.; Ghommam, J.; Derbel, N. Adaptive Backstepping Neural Network Approach to Ship Course Control. In Proceedings of the 8th International Multi-Conference on Systems, Signals and Devices, Sousse, Tunisia, 22–25 March 2011; pp. 1–7. [Google Scholar]
  30. Du, J.; Guo, C.; Yang, C. Adaptive Nonlinear Design of Autopilot for Ship Course Tracking. J. Appl. Sci. 2006, 24, 83–88. [Google Scholar]
  31. Zhang, G.; Zhang, X. A Novel DVS Guidance Principle and Robust Adaptive Path-Following Control for Underactuated Ships Using Low Frequency Gain-Learning. Isa Trans. 2015, 56, 75–85. [Google Scholar] [CrossRef] [PubMed]
  32. Shojaei, K. Neural Adaptive Robust Control of Underactuated Marine Surface Vehicles with Input Saturation. Appl. Ocean Res. 2015, 53, 267–278. [Google Scholar] [CrossRef]
  33. Du, J.; Hu, X.; Sun, Y. Adaptive Robust Nonlinear Control Design for Course Tracking of Ships Subject to External Disturbances and Input Saturation. IEEE Trans. Syst. Man Cybern. Syst. 2017, 1–10. [Google Scholar] [CrossRef]
  34. Norrbin, N.H. Theory and Observations on the Use of a Mathematical Model for Ship Manoeuvring in Deep and Confined Waters; Statens Skeppsprovningsanstalt: Gothenburg, Sweden, 1971. [Google Scholar]
  35. Zhang, X.; Zhang, G. Design of Ship Course-Keeping Autopilot Using a Sine Function-Based Nonlinear Feedback Technique. J. Navig. 2016, 69, 246–256. [Google Scholar] [CrossRef]
  36. Nomoto, K.; Taguchi, K.; Honda, K.; Hirano, S. On the Steering Qualities of Ships. Int. Shipbuild. Prog. 1957, 4, 354–370. [Google Scholar] [CrossRef]
  37. Zhou, Z.; Zhang, B.; Mao, D. Robust Sliding Mode Control of PMSM Based on Rapid Nonlinear Tracking Differentiator and Disturbance Observer. Sensors 2018, 18, 1031. [Google Scholar] [CrossRef] [PubMed]
  38. Lopezfranco, C.; Gomezavila, J.; Alanis, A.; AranaDaniel, N.; Villaseñor, C. Visual Servoing for an Autonomous Hexarotor Using a Neural Network Based PID Controller. Sensors 2017, 17, 1865. [Google Scholar]
Figure 1. Lanxin USV.
Figure 1. Lanxin USV.
Sensors 18 01889 g001
Figure 2. Vector propulsion system.
Figure 2. Vector propulsion system.
Sensors 18 01889 g002
Figure 3. Multi-sensor structure.
Figure 3. Multi-sensor structure.
Sensors 18 01889 g003
Figure 4. The sea state of the field experiment.
Figure 4. The sea state of the field experiment.
Sensors 18 01889 g004
Figure 5. Schematic diagram of plane motion.
Figure 5. Schematic diagram of plane motion.
Sensors 18 01889 g005
Figure 6. A schematic diagram of the vector thrust distribution.
Figure 6. A schematic diagram of the vector thrust distribution.
Sensors 18 01889 g006
Figure 7. The comparison results of the zig-zag test.
Figure 7. The comparison results of the zig-zag test.
Sensors 18 01889 g007
Figure 8. The comparison results of the turning test.
Figure 8. The comparison results of the turning test.
Sensors 18 01889 g008
Figure 9. The comparison results of the zig-zag test.
Figure 9. The comparison results of the zig-zag test.
Sensors 18 01889 g009
Figure 10. The comparison results of the turning test.
Figure 10. The comparison results of the turning test.
Sensors 18 01889 g010
Figure 11. External disturbance curve.
Figure 11. External disturbance curve.
Sensors 18 01889 g011
Figure 12. The numerical simulation of course keeping.
Figure 12. The numerical simulation of course keeping.
Sensors 18 01889 g012
Figure 13. The field experiment for course keeping.
Figure 13. The field experiment for course keeping.
Sensors 18 01889 g013
Table 1. Specific Parameters of the Lanxin Unmanned Surface Vehicle (USV).
Table 1. Specific Parameters of the Lanxin Unmanned Surface Vehicle (USV).
ItemValue
Length between perpendiculars7.02 m
Breadth2.60 m
USV speed (max)35 kn
Draft (full load)0.32
Block coefficient0.6976
Displacement (full load)2.73 m3
Rudder area0.2091 m2
Propulsion angle (max)35 degrees
Distance Between gravity and center0.35 m
Pitch ratio0.3
Disk surface ratios0.516
Diameter of the propeller0.46 m
Table 2. The coefficients of the propeller thrust coefficient expression.
Table 2. The coefficients of the propeller thrust coefficient expression.
Item θ = 0.5 θ = 0.8 θ = 1.1
A 00 −0.1677−0.0517−0.2191
A 01 0.1747−0.03150.3013
A 02 −0.6720−0.5822−0.7309
A 10 0.80420.5853−0.8502
A 20 −0.1437−0.1026−0.1080
A 30 000
A 11 −0.8853−0.3381−1.0738
A 12 0.91300.66540.9908
A 21 0.34220.14170.3481
A 22 −0.3276−0.2215−0.3322
Table 3. ITAE index.
Table 3. ITAE index.
ItemValue
PID37.98
PD38.6

Share and Cite

MDPI and ACS Style

Mu, D.; Wang, G.; Fan, Y.; Sun, X.; Qiu, B. Modeling and Identification for Vector Propulsion of an Unmanned Surface Vehicle: Three Degrees of Freedom Model and Response Model. Sensors 2018, 18, 1889. https://doi.org/10.3390/s18061889

AMA Style

Mu D, Wang G, Fan Y, Sun X, Qiu B. Modeling and Identification for Vector Propulsion of an Unmanned Surface Vehicle: Three Degrees of Freedom Model and Response Model. Sensors. 2018; 18(6):1889. https://doi.org/10.3390/s18061889

Chicago/Turabian Style

Mu, Dongdong, Guofeng Wang, Yunsheng Fan, Xiaojie Sun, and Bingbing Qiu. 2018. "Modeling and Identification for Vector Propulsion of an Unmanned Surface Vehicle: Three Degrees of Freedom Model and Response Model" Sensors 18, no. 6: 1889. https://doi.org/10.3390/s18061889

APA Style

Mu, D., Wang, G., Fan, Y., Sun, X., & Qiu, B. (2018). Modeling and Identification for Vector Propulsion of an Unmanned Surface Vehicle: Three Degrees of Freedom Model and Response Model. Sensors, 18(6), 1889. https://doi.org/10.3390/s18061889

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

Article Metrics

Back to TopTop