Design of an adaptive backstepping controller for autoberthing a cruise ship under wind loads
 DOI : 10.2478/IJNAOE20130184
 Author: Park JongYong, Kim Nakwan
 Organization: Park JongYong; Kim Nakwan
 Publish: International Journal of Naval Architecture and Ocean Engineering Volume 6, Issue2, p347~360, 30 June 2014

ABSTRACT
The autoberthing of a ship requires excellent control for safe accomplishment. Crabbing, which is the pure sway motion of a ship without surge velocity, can be used for this purpose. Crabbing is induced by a peculiar operation procedure known as the pushpull mode. When a ship is in the pushpull mode, an interacting force is induced by complex turbulent flow around the ship generated by the propellers and side thrusters. In this paper, three degrees of freedom equations of the motions of crabbing are derived. The equations are used to apply the adaptive backstepping control method to the autoberthing controller of a cruise ship. The controller is capable of handling the system nonlinearity and uncertainty of the berthing process. A control allocation algorithm for a ship equipped with two propellers and two side thrusters is also developed, the performance of which is validated by simulation of autoberthing .

KEYWORD
Autoberthing , Crabbing , Adaptive backstepping control , Control allocation

INTRODUCTION
Berthing is the process of positioning and mooring a ship beside a quay, jetty, or floating dock, usually for the purpose of loading or unloading. For large ships such as a container or a cruise ship, berthing is done with the aid of tug boats. When the ship approaches the berthing position, forward tug boats are used to hold the bow to prevent the ship from contacting the quay. Aft tug boats are then used to push the ship towards the quay. If the lateral speed of the ship is higher than the desired speed, the tug boats would be used to retard it. By careful operation of the propellers and rudder, the ship is positioned a few meters away from the quay, and thereafter brought nearer by means of tug boats and mooring ropes. The entire operation is actually very complex and time consuming.
Crabbing is the pure sway motion of a ship without surge velocity, and is induced by a peculiar operation method known as the pushpull mode. The pushpull mode is induced by the combined manipulation of the main propeller and side thrusters. The two propellers are made to generate the same amount of thrust while rotating in opposite directions, thereby exerting a yawing moment on the vessel without inducing longitudinal motion. By the simultaneous operation of the side thrusters, the pushpull mode is implemented, resulting in the generation a large lateral force. When a ship is in the pushpull mode, an interaction force is induced by complex turbulent flow around the ship generated by the propellers and side thrusters. Crabbing is a slow sway motion and can thus be applied to berthing if the ship is equipped with propellers and the thrusters are close to the berthing position. Some requirements for effective crabbing and berthing were presented by Quadvlieg (1998), namely the maintenance of a lateral speed of 0.25(
m/s ) against a Beaufort 6 wind. For cruise ships, the requirement for effective berthing is often the ability to berth and unberth in a wind speed corresponding to Beaufort 7 without the aid of tug boats. Without the use of tug boats, it is difficult to manually manipulate the propellers and side thrusters during the berthing operation, which has given rise to the need to apply a controller.There have been some recent studies on crabbing. Lee et al. (2000) simulated crabbing of a ship with twin rudders and twin skegs. However, whereas crabbing involves lowspeed maneuvering and pure lateral motion, their mathematical model was a conventional model used for simulating longitudinal maneuvers. According to Yoshimura (1988), the forces acting on a ship are generated by the cross flow dragwhich is proportional to the square of the lateral speedand not by the lift in low speed manoeuver. Experiments on the hydrodynamic forces of a ferry in crabbing motion involving pure sway and pure yaw without longitudinal velocity were conducted by Yoo et al. (2006). Based on experimental investigation, Quadvlieg and Toxopeus (1998) proposed simple techniques for calculating the interaction force, but the techniques were shown to be inadequate for estimating crabbing capability. Studies of autoberthing are mainly based on the optimal control theory, the neural network theory, or an expert system (Yamato et al., 1993; Hasegawa and Kitera, 1993; Im and Hasegawa, 2001). A system for planning the optimal berthing path of a ship was developed by Djouani and Hamam (1995). Yamato et al. (1993) and Hasegawa (1994) proposed autoberthing controllers based on the neural network theory for considering the nonlinear characteristic of lowspeed ship maneuvering. Im and Hasegawa (2001) proposed a parallel hidden layer neural network controller for improving autoberthing performance beyond that achieved by a conventional neural controller. They, however, focused on berthing with the aid of a rudder, which is not realistic for extremely lowspeed maneuvering, wherein the rudder is incapable of generating sufficient lateral force and yaw moment. There has hardly been any study that took into consideration the interaction force and severe environmental disturbances such as strong winds. Moreover, there have been few studies on autoberthing by crabbing.
An autoberthing controller should be able to handle the uncertainty in modeling the interaction force, as well as the nonlinearity of the force acting on the hull during crabbing. Adaptive control is a control method that adapts to a system with uncertainty, whereas backstepping control was developed by Kokotovic (1992) for designing the controller of nonlinear dynamic systems. Adaptive backstepping control is a combination of the two methods and is used for systems characterized by nonlinearity and uncertainty, such as those used for crabbing.
In this paper, we develop a mathematical model of crabbing for enhanced berthing simulation. The interaction force and wind force are incorporated in the model, and are assumed to be unknown disturbances during the design of an adaptive backstepping controller that can adapt to unknown disturbances. The Lyapunov stability theory is employed in the design of the controller, which is verified by simulations of crabbing, berthing, and unberthing under wind conditions. The paper is finally concluded by drawing some conclusions.
MATHEMATICAL MODEL
> Overview
The coordinate system used in this study is shown in Fig. 1. It consists of the bodyfixed coordinate
O_{b} −x_{b} y_{b} z_{b} and spacefixed coordinateO_{s} −x_{s} y_{s} z_{s} . The origin of the bodyfixed coordinate is located at the middle of the ship, on the centerline. The zaxis is positive downward and all angles are positive in the clockwise direction. The kinematic equations of the motion on the horizontal plane can be expressed aswhere
u andv are respectively the linear surge and sway velocities defined in the bodyfixed coordinate,r is the yaw rate, andβ is the drift angle. The orientation of the bodyfixed coordinate relative to the inertia coordinate is described by the Euler angles rotated by the yaw angleψ .The design of the autoberthing controller must be preceded by the mathematical modeling of the ship motion. The three degrees of freedom equations of the surge, sway, and yaw can be expressed as follows:
The terms on the righthand side of Eq. (2) represent the external forces acting on the ship. In this study, the following modular type mathematical model is used to express the external forces:
The modular type mathematical model is expressed in terms of the force acting on the hull, the propeller force, the rudder force, the side thruster force, the interaction force, and the wind force. The rudder is excluded from the control algorithm because it is incapable of generating sufficient force and moment when the longitudinal speed of the ship is almost zero. Here, the subscript
W represents wind. In the harbor, the effects of the current and wave are negligible; hence, only the wind force is considered in the environmental disturbance. When the speed of the ship is very low or nearly zero, a high lateral force is generated by the propellers and side thrusters. One main propeller astern is operated at maximum power and the other is operated ahead to balance the longitudinal force to the required value during berthing. The pushpull mode is thus implemented because one main propeller pushes the ship while the other pulls it. This pushpull mode is illustrated in Fig. 2.> Hull force
The hull force consists of the added mass term, hydrodynamic damping term, and restoring term. The added mass is the pressureinduced force due to the motion of the ship, and it is proportional to the acceleration. The restoring term is the hydrostatic force due to the weight and buoyancy of the ship, and does not have a horizontal component. The added mass is determined by the empirical formula of Motora (1959). The hydrodynamic damping is produced by the wave, drag, and vortex shedding effect. During the crabbing motion for berthing, the surge velocity is nearly zero and the sway velocity is increased by the side thrusters and propellers in the pushpull mode. The sway velocity is also very small, which makes the wavemaking damping effect negligible. The inclusion of these terms in a mathematical model for predicting the lateral force affords greater accuracy than the conventional model. A static sway test, dynamic sway test, static yaw test, and dynamic yaw test without forward speed were conducted by Yoo (2006) and the results were used to develop a mathematical model of the force acting on the hull during crabbing. The equations of the model are as follows:
> Propeller force
The thrust of the propeller is a function of its thrust coefficient, rotational speed, and diameter. The sign function is adopted to consider the direction of the thrust of the propeller:
where
D_{P} ,n_{P} are respectively the diameter and rotational speed (rpm ) of the propeller. The loss in the thrust during its transmission to the ship is given by the thrust deduction coefficientt_{P} .where
y_{p} is the lateral distance between the propeller and the center of gravity of the ship.> Side thruster force
The side thruster is used when the rudder is incapable of generating sufficient lateral force and yaw moment at low speeds. In the pushpull mode, the propellers produce a large yaw moment, and the side thruster is therefore used to compensate for the yaw moment. The thrust of the thruster is a function of its thrust coefficient, rotational speed, and diameter. The direction of the thrust of the side thruster during berthing is different from that during unberthing; hence, a sign function is adopted:
where
D_{ST} ,n_{ST} are respectively the diameter of the side thruster and its rotational speed (rpm ). The loss in the thrust of the side thruster during its transmission to the ship can be expressed as the thrust deduction coefficient of the thruster.where
y_{ST} is the longitudinal distance from the side thruster to the center of gravity of the ship.> Interacting force
In the pushpull mode, the fluid particles accelerated by the reverse rotation of the propeller impact the stern of the ship, whereas those accelerated by the side thruster change the flow field on the lateral side of the ship. Quadvlieg (2011) conducted a pushpull mode model test and observed a difference between the measured force and the sum of the forces of the propeller and the side thruster. The difference is defined as the interaction force, which is expressed as follows:
According to Eq. (9), the interaction force is a function of the thrust of the reverse propeller and the force generated by the side thrusters. The direction of the interaction force is the same as the crabbing direction.
a_{I} is the interaction force coefficient, which is dependent on the shape of the stern; andx_{I} is the longitudinal distance from the center of gravity of the ship to the point of action of the interaction force.> Wind force
The wind force is considered as an environmental disturbance. The wind force is considered in the berthing problem because it has been suggested that the requirements for effective berthing are dependent on the wind conditions Quadvlieg (1998). The equations of the wind force can be obtained as expressed by Eq. (10).
where
A_{F _ ab}, A_{S _ ab} are respectively the lateral and longitudinal projection areas of the superstructure, andU_{R _W} is the relative velocity vector of the wind and the ship. The relative velocity vector can be obtained by coordinate transform of the wind velocity vector, which is described in the spacefixed coordinate.The wind coefficients
C_{X_W}, C_{Y_W}, C_{N_W} are functions of the shape and area of the superstructure. Fujiwara (2001) integrated the results of wind tunnel tests conducted on 34 types of ships such as VLCC, PCC, and LNG carrier, and used the stepwise method of linear multiple regression analysis to propose a method for estimating the wind coefficient. The wind coefficient estimation equations can be expressed as follows:where
ψ_{R _W} is the wind direction in the bodyfixed coordinate. The wind direction in the bodyfixed coordinate system used in this paper can be defined as follows:The coefficients
X_{0}, X_{1}…,N_{3} can be obtained by Fujiwara’s empirical formula. There are some differences between the wind coefficients obtained by the empirical formula and those obtained from the wind tunnel experimental data for a cruise ship; hence, the coefficients are reproduced by the least square method. The wind coefficients obtained from the experimental data and the reproduced values are shown in Fig. 3.CONTROLLER DESIGN
A backstepping controller can deal with the system nonlinearity of the berthing maneuver. An adaptive controller can also be used to compensate for the unknown disturbance. Based on the work of Fossen (1994), an adaptive backstepping controller for the vector system is designed in this section.
> Adaptive backstepping controller
Three degrees of freedom nonlinear equations of the planar motion can be expressed as
where
M is the inertia matrix that includes the added mass and added mass moment of inertia,C (υ ) are the Coriolis and centripetal matrices of the rigid body,τ is the control force vector, andW is the disturbance caused by the interaction and the wind.υ = [u, v, r ]^{T} is the velocity vector in body fixed coordinate andη =[xs, ys, ψ ]^{T} is the position vector in space fixed coordinate. For crabbing, the matricesM ,C (υ ),D (υ ) can be constructed as follows:The spacefixed velocity vector and acceleration are given by
Using Eqs. (16) and (17), Eq. (14) can be expressed as
The matrices can be defined as
Eq. (18) can thus be reformulated as
The first backstepping variable is defined as the error between the position vector and reference model position vector for the ship to track:
The virtual control is defined as
where
α _{1} is the function used to stabilizez _{1}, and can be chosen by the designer.z _{2} is the second backstepping variable, which is the error between the virtual control and the stabilizing function. The stabilizing function is chosen as follows:Λ is the positive definite matrix. The
z _{2} system can be obtained by differentiating Eq. (23) with respect to time:A Lyapunov candidate function is
where
K_{p} and Γ are the positive definite matrices. The empirical formula of the interaction force has both modeling and parameter uncertainties, and it is difficult to use it to determine the wind speed and direction with high accuracy. In this study, it is assumed that the interaction and wind forces are unknown disturbanceW . Unknown disturbanceW is estimated by adaptive law. Estimated disturbance is defined asŴ . The parameter estimation error term is included in the Lyapunov candidate function and can be expressed asThe time differential of the Lyapunov candidate function is given by
By substituting Eq. (25) into Eq. (28), we obtain
The control law is chosen as follows:
K_{D} is the positive matrix. By substituting Eq. (30) into Eq. (29), we obtainThe adaptive law is defined as
The time differential of the Lyapunov candidate function can thus be made negative semidefinite as follows:
By LaSalle’s theorem, the backstepping variables 1 z and 2 z approach zero as the time tends to infinity. Fig. 5 is a block diagram of the adaptive backstepping control.
> Control allocation
The determined control force vector is allocated to the propellers and side thrusters. The rudder is incapable of generating sufficient lateral force and yaw moment during crabbing and is therefore excluded from the control allocation. The relationship between the control force vector and the thruster forces can be expressed as follows:
where the subscripts
P _P andP _S respectively represent the port side propeller and the starboard propeller, and the subscriptsST1 andST2 are used to distinguish the two bow side thrusters. The second entry in the control force vector represents the force in the lateral direction. The control force in the lateral direction is equally allocated to the side thrusters.Based on Eqs. (8), (9), and (35), the command for the rotational speed of the side thrusters can be expressed as
The propellers are used to compensate for the yaw moment generated by the side thrusters and balance the longitudinal force. Based on Eqs. (5), (6), and (34), the command for the rotational speed of the propellers can be expressed as
In this paper, it is assumed that the dynamics of the propellers and side thrusters are in accordance with first order actuator dynamics, and that the time constant
T is 1.> Reference model
A reference model is used to prevent wastage of control effort as a result of excessive position error in the initial phase. In this study, a second order dynamic system is established for use in referencing the model dynamics. The reference model dynamics can thus be written as
η_{d _ C} is the command position (i.e., the berthing position) vector andη_{d} is the reference model position vector for the ship to track. Initial reference model position vector is same with the initial position vector of a ship. When the command position vectorη_{d _ C} is set, the reference model position vectorη_{d} is calculated by the reference model dynamics in realtime. If the reference model dynamics is stable, the reference model position vectorη_{d} converges to command position vectorη_{d _ C} as time goes on. If the adaptive backstepping controller designed well, the position vector of a ship converges to reference model position vectorη_{d} , so the position vector of a shipη is converges to command position vectorη_{d _ C} as time goes on. The natural frequencyω_{n} and damping ratioζ of the reference model are respectively selected as 0.005 and 0.9.SIMULATIONS
A cruise ship equipped with twin propellers and two side thrusters is used for the simulation. The principal dimensions of the ship are given in Table 1.
The control gain matrices are diagonal and their components are given in Table 2. The values of control gain matrices are chosen by trial and error in order to obtain the satisfactory berthing performance.
Three scenarios were simulated, namely crabbing in wind, autoberthing in wind, and autounberthing in wind. The establishment of the wind condition is based on the study of Quadvlieg (1998). It is assumed that the behavior of wind velocity and direction can be specified as a normal distribution function.
> Performance criteria 1
The simulation of crabbing in wind is used to check whether the cruise ship satisfies performance criteria 1 given by Hooren (1985), namely that a lateral speed of 0.25
m/s should be maintained against a 22knots wind. The initial position vector of the ship is [0m 0m 0 deg]^{T} , the desired velocity vector isη_{d } = [0m /s 0.25m /s 0 deg /s ]^{T} , and the reference model position vector isη_{d } = [0m 0.25 ×tm 0 deg]^{T} . The mean wind speed is 22knots and standard deviation of wind speed is 1.9knots . The mean wind direction is 270°, which means that the wind force acts on the starboard side of the ship. The standard deviation of wind direction is 10°. The simulation time is 1,500s and the results are shown in Figs. 57.It can be observed from Fig. 5 that the ship trajectories are in good agreement with the reference model position vector. Fig. 6 shows the velocity vector of the ship. The lateral velocity is 0.25
m/s in the steady state, which indicates that the ship satisfies the crabbing criteria. Fig. 7 shows the control inputs, which means the rotational speed of propellers and side thrusters. The wind speed and direction has a noise, so the control inputs fluctuated with some noise.> Performance criteria 2
The simulation of autoberthing in wind is used to verify whether the ship satisfies performance criterion 2, namely the ability to berth and unberth under a wind speed corresponding to Beaufort 7 without the aid of tug boats. The initial position vector of the ship is [0
m 100m 0 deg]^{T} and the command position (i.e., the berthing position) vectorη_{d _ C} is [0m 0m 0 deg]^{T} . The mean wind direction and speed are 90° and 35knots , respectively. The standard deviation of wind direction is 10° and the standard deviation of wind speed is 1.9knots .Although a wind force acts on the ship, the simulation results indicate stable berthing. Moreover, the heading angle error and the longitudinal error is almost zero. The maximum lateral velocity is about −0.2
m/s . From the control input responses, it can be seen that the side thrusters generate port side thrust and negative moment. The port propeller rotates normally and the starboard propeller inversely to balance the longitudinal force and yaw moment. After completion of the berthing maneuver, there are still steady control inputs to maintain the reference model position against the wind force.> Autounberthing simulation
The simulation of autounberthing is used to verify the control allocation algorithm. If the control allocation is not designed properly, inappropriate results may be obtained depending on the direction of the lateral motion. The wind condition is the same as that of scenario 2. The initial position vector and initial reference model position vector are both [0
m 0m 0deg ]^{T} . The command position vectorη_{d_C} = [x _{d_C}, y _{d_C}, ψ _{d_C} ]^{T} is [0m 100m 0deg ]^{T} .It can be seen from Fig. 11 that the designed control algorithm produced stable unberthing. The rotational speed of the side thrusters is negative and the lateral wind force is positive, which indicates that the rotational direction of the side thrusters must be negative to compensate for the wind force as the ship approaches the reference model position.
CONCLUSIONS
The equations of the crabbing maneuver of a ship were derived in this paper. The hull force, propeller force, side thruster force, interaction force, and wind force were considered as external forces acting on the ship. Because crabbing is an extremely low speed maneuver, the rudder was excluded from the control algorithm. Errors may be present in the model of the interaction force and the wind force of actual berthing. The interaction and wind forces were thus defined as uncertainty terms in the design of the controller. The terms were estimated by the adaptive law and compensated for by a control input. The other terms were controlled by the backstepping control method. The results of autoberthing simulations that were performed to verify the derived control law confirmed the effectiveness of the designed controller.

[]

[Fig. 1] Coordinate systems.

[]

[]

[Fig. 2] Interacting force induced by pushpull mode.

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[Fig. 3] Wind coefficients.

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[Fig. 4] Block diagram of adaptive backstepping control.

[]

[]

[]

[]

[]

[]

[Table 1] Principal dimensions of cruise ship.

[Table 2] Components of control gain matrices.

[Fig. 5] Trajectory and position vector of ship (scenario 1).

[Fig. 6] Velocity vector of ship (scenario 1).

[Fig. 7] Rotational speed of propellers and side thrusters (scenario 1).

[Fig. 8] Trajectory and position vector of ship (scenario 2).

[Fig. 9] Velocity vector of ship (scenario 2).

[Fig. 10] Rotational speed of propellers and side thrusters (scenario 2).

[Fig. 11] Trajectory and position vector of ship (scenario 3).

[Fig. 12] Velocity vector of ship (scenario 3).

[Fig. 13] Rotational speed of propellers and side thrusters (scenario 3).