Dynamical Sliding Mode Control for Robust Dynamic Positioning Systems of FPSO Vessels

: The conventional proportional derivative (PD) control algorithm with appropriate gain scheduling is generally applied to a dynamic positioning (DP) system. However, ﬁnding appropriate gains through gain scheduling makes the DP system more complicated. A sliding-mode control algorithm controls an arbitrary point, such as the turret system on a ﬂoating production storage and ofﬂoading (FPSO) vessel. This algorithm was developed for DP and can be applied to FPSO vessels considering the uncertainty of the vessel dynamics, unknown time-varying environmental disturbances, and transient performance. To control an arbitrary point on the FPSO vessel using a DP controller, the Jacobian matrix in the kinematic equation is modiﬁed to present the arbitrary point in the control. The Lyapunov stability theory is applied in the design of the SM control algorithm to provide robustness to the control system. A time-domain simulation tool was developed to verify the effectiveness of the proposed SM control algorithm. The performance of the control algorithm was evaluated numerically to address its efﬁcacy. The results were compared with those obtained using the conventional PD control algorithm.


Introduction
Recently, the demand for advanced marine mechatronics systems has been increasing with the current situation where the region of human production activities is expanding to the deep ocean and polar regions.The dynamic positioning (DP) system is the representative marine mechatronics system.This system adjusts the positions and forces of a number of propulsion systems installed in the hulls and operates along a set route during vessel navigation.Vessels equipped with DP systems are operated in the offshore oil and gas industry because they have economic advantages when compared with traditional methods for maintaining the location of floating structures such as mooring systems or jacket structures.Specifically, DP systems have been applied in marine drilling, dredging, lifting, pipelaying, subsea installation and maintenance, and for transportation of workers between marine structures.These systems have become a major part of the marine industry [1][2][3][4].
The optimal performance of the control system is essential for the effective operation of DP systems used for various purposes in the ocean.When the DP system was first developed in the early 1960s, a single-input single-output (SISO) proportional-integraldifferential (PID) controller with a low-pass filter [5] was applied to independently control the surge, sway, and yaw motions of vessels.Balchen et al. [6] and Grimble et al. [7] applied multi-variable optimal control and Kalman filter to DP systems, which enabled reliability and stability of such systems to be used in a wide range of areas.Based on the nonlinear control theory, Lee et al. [8], Fossen et al. [1] and Kim et al. [9] developed weathervaning control to automatically set the heading angle to the direction of the wave loads acting on the ship to reduce fuel consumption and CO 2 emissions.Kim et al. [10] demonstrated the possibility of using weathervaning control applied to the process of loading and offloading for floating production storage and offloading (FPSO) with tandem arrangement and shuttle tanker.Sørensen [11] introduced DP control techniques for the application of marine mechatronics systems.Yu et al. [12] applied the H ∞ control technique to DP control and compared its performance with linear-quadratic-Gaussian control (LQG) through numerical simulations.Lee et al. [13] applied the theory of linear quadratic regulator (LQR) and LQR with an error integrator (LQI).Kim et al. [14] also developed a PID controller that included an anti-windup controller that restrained the propensity of divergence of the integral controller in the DP system of a shuttle tanker.Jeon et al. [15] studied the gain scheduling of the PID controller using a fuzzy system.Aleksander et al. [16] presented a dynamic positioning scheme based on a model-predictive control algorithm that combines positioning control and thrust allocation to theoretically yield a near-optimal controller output.Borkowski [17] suggested designing an expert system that includes a function to automatically stabilize ship's course with the inference engine.Additionally, Peng et al. [18] introduced a cooperative control scheme for the dynamic positioning of multiple offshore vessels based on a dynamic surface control technique with the input-to-state stability of a closed-loop network system.
Robust DP control techniques applicable to a variety of vessels and work scenarios are essential because of the recent demands in marine operations with respect to increasing types and frequencies.DP systems are based on conventional controllers; proportional derivative (PD) controller is a widely used linear controller.Conventional PD controllers typically use gain scheduling based on each characteristic of marine operation scenarios and ocean environmental changes to ensure DP performance.However, there are some problems related to the application of a linear PD controller to the DP system.Gain scheduling should be adjusted for each vessel model and ocean environment while operating the DP system.If the gain values are not appropriate in a specific vessel model and ocean environment, the conventional PD controller sometimes results in control difficulty with unsatisfactory performance in nonlinear plant controls, especially in FPSO vessels with turret mooring systems.
A sliding mode (SM) controller was applied to the DP system to overcome the limitation of the conventional linear controller.The SM controller approach enables the realization of control performance and stability requirements for the desired position.Because the SM controller uses vessel model information, an estimated mass is used as a parameter.Furthermore, the same performance is ensured for different loading conditions without requiring a gain-scheduling approach, and the SM controller can be easily adjusted with simple equations.Therefore, the SM controller can be applied in a wide range of environmental conditions without performance degradation.Tannuri et al. [19,20] applied an SM controller to the application of the motion control of the turret-moored offshore platform along the direction of the horizontal plane, such as surge, sway, and yaw.Kim et al. [21] also applied an SM controller to the depth adjustment and maneuverability control of an unmanned submarine and demonstrated its effectiveness using simulation.Liang et al. [22] presented a finite-time observer-based adaptive sliding mode output feedback control for the DP of ships considering input saturation, unmeasured states, and unknown environment disturbances.Zhao et al. [23] proposed a robust adaptive terminal sliding mode controller for the DP of a semi-submersible offshore platform using a state feedback controller and a robust adaptive terminal sliding mode compensator.Chen et al. [24] developed a novel iterative sliding mode-based output feedback controller for the DP of a ship subjected to ocean environmental disturbance with finite-time state observer.Agostinho et al. [25] proposed the application of a controller with sliding mode control technique to the DP of a floating vessel using numerical simulations and experimental validations.
In this study, a DP system based on the SM controller was developed for the application of FPSO carrying out turret connection work in the sea.The FPSO is equipped with six azimuth thrusters and a turret placed at the bow position of the vessel to carry out access work.A kinematic model is constructed that transforms the center of motion of the vessels because the motion of the vessel should be controlled around the turret installation point located in the bow position.The SM controller, a nonlinear control technique, is used to control the DP system to ensure control performance by selecting control parameters intuitively considering the response characteristics of vessels.In addition, a MATLAB/Simulink-based simulator, which can quantitatively simulate the motion of vessels, was used to evaluate the validity of the proposed control algorithms.The SM controller was designed by applying the Lyapunov stability theory to ensure the robustness of the controller.The effectiveness of the designed SM controllers was evaluated using time-domain simulations and compared with the results of a typical proportional-derivative (PD) controller.
Section 2 describes the configuration of simulators for the performance evaluation of the DP system and the theory of vessel motion utilized therein.Section 3 presents the design of the SM controllers that enables DP control around specific vessel locations.Section 4 presents the results of the DP control simulation in marine environment conditions and Section 5 summarizes the study and presents the conclusions.

Methodologies for the DP Control Application
This section describes the methodologies for the DP system.The DP system consists of sensors, controllers and filtering algorithms, and propellers.The sensors are used to measure the positions and angles of the vessel, while the control and thrust allocation algorithms compute the desired forces to be delivered by each propeller to counteract environmental loads such as wind, waves, and current.Figure 1 shows a block diagram of a DP system illustrating the connections between the controller, thruster allocation algorithm, propeller, vessel, sensors and wave filter.A wave filter suppresses the highfrequency motion from the measurements, and a thrust allocation algorithm distributes the desired forces to the propellers equipped in the vessel.To evaluate the motion and control response of vessels equipped with DP systems, a vessel operation simulator with three degrees of freedom (DOF) motions including surge, sway, and yaw was developed based on MATLAB/Simulink.This section describes a mathematical model that calculates the motion response of a marine surface vessel, considering the effects of sea loads such as from waves, currents, and wind loads.Vessels operating in the sea engage in six DOF motions, which include translational motions such as surge, sway, heave, and rotational motions including roll, pitch, and yaw.Specifically, the DP system in this study considers three DOF horizontal plain models, namely surge, sway, and yaw.
Equation (1) describes the surface vessel motion with a state variable vector representing the position (x, y) and angle (ψ) of the vessel, as defined in the global reference frame {e} shown in Figure 2. The vessel motion consists of a low-frequency motion η L , which is caused by the effects of wave drift force, current force, wind force, etc., and a high-frequency motion η H , which is caused by the wave exiting force as follows: (1) Equations ( 2) and ( 3) describe the 3 DOF equations of the motion of the vessel and represent the kinematic transformation between the global reference frame and the bodyfixed coordinate system for numerical simulations as follows: .
where ν is the velocity vector of the vessel defined in the fixed coordinate system such as the velocity of the surge, sway, and yaw directions.R(ψ L ) is the Euler angle transformation matrix that converts the velocity of the body coordinate system into the velocity of the global reference coordinate system.In Equation ( 3), M rb is the rigid body mass matrix, and M a is the added mass of the fluid around the hull at a low frequency of the vessel.C rb (ν) and C a are the matrices that express the Coriolis and centripetal forces, respectively, expressing the effect of the hull mass and added mass.Furthermore, ν r is the relative velocity defined by ν r = ν − ν c , and ν c is the current velocity.D 0 is a linear damping coefficient matrix at low frequencies, characterized by D 0 > 0. dNL(ν r ) is a nonlinear damping force vector that includes the current force and cross-flow damping force.τ wave and τ wind are the wave and wind loads, respectively, at low frequency; they are external forces of the marine environment.τ is the control force generated through the operation of the controller and the propeller.All parameters of matrices are described in Appendix A.
The vessel exposed above the waterline is subjected to sea wind loads according to its shape, as expressed in Equation ( 4) below, where ρ a is the air density and U w is the sea wind speed.A f w and A sw are the projection areas along the surge and sway directions of the vessel, respectively.C Xw , C Yw , and C Nw represent the wind force coefficients in the surge, sway, and yaw directions, respectively.γ w is the incident angle of sea wind.
The wave drift force acting on the vessel is described by the quadratic transfer function (QTF) in Equation (5).A j and A k refer to the amplitude of the j th and k th wave frequencies in the components of the N wave frequencies, which are calculated when generating the irregular wave from the given wave spectrum.The nonlinear transfer function is a function of complex numbers representing the wave drift force calculated through a numerical analysis program using the high-order boundary element method.QTF i R and QTF i I represent the real and imaginary parts of the quadratic transfer function, respectively [26].ω j and ω k are the j th and k th wave frequencies, respectively.ω d and d are the differences of wave frequencies and phases defined by ω d = ω i − ω k and d = i − k , respectivley.Finally, t denotes the time when the ocean waves are generated.
The high-frequency (or wave frequency) motion can be expressed by the response amplitude operator (RAO) as shown in Equation ( 6): where RAO is also expressed as a harmonic function representing the real and imaginary parts of the RAO.The noise generated using the sensor modules in the simulators for evaluating DP control responses, as shown in Figure 1, using the primary Gauss-Markov process is given as follows: .
where y i is the measured physical value, w is the Gaussian white noise, and µ is the inverse of the time constant that produces a long-term error in the physical measurement sensor.
Fossen [27] describes the configuration of the dynamical term.An extended Kalman filter (EKF) is applied to the simulator to separate the lowfrequency motion of the vessel model, including the low-frequency motion of the ves-sel, wave frequency, and sensor noise.The estimation algorithm for uncertain external forces was designed in this study by adapting the EKF with a disturbance model of the uncertain external force induced by wave, current, and wind.The EKF is a nonlinear ad-aptation of the Kalman filter, which observes a dynamic system in which zero-mean, Gaussian white noise is assumed in both the dynamics and measurement.The model as a continuous form that is observed by the extended Kalman filter is described by .
where x(t) is the state vector of the system described by x(t) = ξ(t), .
ξ(t) is the wave perturbation vector, .ξ(t) is the wave frequency vector, η(t) is vessel's motion vector, .η(t) is the vessel's velocity vector, and F ext is the uncertain external force vector induced by ocean environmental loads.u(t) is the input vector described by u(t) = τ, and z(t) is the measurement vector measured by motions of the vessel described by z(t) = η m (t).f is the nonlinear state dynamics, h is the nonlinear output function, G(t) is a vector that describes the noise of each state, and w(t) and v(t) are the process and measurement noises, respectively, described by zero-mean Gaussian white noise.The Equation ( 8) can be transformed into Equation ( 9) as discretized form with Euler discretization as follows: The EKF iterates through the prediction and update steps for each discrete time to estimate the states of the system as shown in Figure 3.The prediction step calculates a predicted state and covariance.The update step calculates the Kalman gain to adjust the covariance of each state and the estimated value based on the difference between the estimated and measured states.Finally, the covariance and state estimates of the next time step can be updated utilizing a Jacobian to linearize the conditions.In this procedure, the algorithm calculates the estimated state vector xk described by xk = ξ(k),  The propulsion system is modeled using a linear differential equation, and the maximum speed rate of each propeller characteristic is limited.In addition, the thrust distribution algorithm that distributes the control forces calculated from the controller to six azimuth propulsion systems is applied based on the Lagrange multiplier theory.Fossen [27] described the details of the extended Kalman filter, propulsion system, and thrust distribution algorithm.

Design of a Sliding Mode Controller
This section describes the design of the sliding mode controller for the dynamic motion control of a vessel model performing turret access operations in the sea.As shown in Figure 4, it is assumed that the turret (p {b} ) used to carry out the access work was located at the bow position of the vessel model.Equation (10) describes the kinematic equations for converting the motion of the vessel to ensure that the control algorithm places a specific point z on the target point z d .Equation (11) describes the kinematic transformation matrix T ψ, p {b} according to the position of the turret P {b} = x p y p T as defined in the bodyfixed coordinate system {b}.Additionally, Equation ( 12) describes the dynamic equation of the motion of the vessel with unknown external forces for the design of the SM controllers.
where M is the mass of a vessel including the low-frequency added mass described by M = M rb + M a , u is the control vector, and w is the uncertainty variable vector including the dynamic model and environmental external force.J ψ, p {b} is the transformation matrix that converts the body-fixed coordinate to the global reference coordinate system.Figure 5 describes the overall procedure of sliding model control.The control objective in this study was to set the turret position at a specific point and to ensure that the velocity of the vessel is zero; this can be expressed by Equation (13) as To achieve this control objective, the system error e can be defined using Equation ( 14) and the sliding surface for SM control can be defined using Equation ( 15) where I λ is the gain for feedback control when the system state is on a sliding plane, and all elements of I λ should have positive values.Therefore, the system error e automatically converges to zero when the s = 0 condition is satisfied.
For the design of the robust controller, the Lyapunov function of the control system with the sliding plane variable can be used as shown in Equations ( 16) and (17), which describes the derivative of the Lyapunov function to determine the stability of the control system.
The input u of the SM controller consists of the feedback controller u f and the switching controller u s .Equation ( 18) describes the feedback controller u f and Equation (19) describes the relation between the derivative of the Lyapunov function and u s .Equation ( 20) describes the switching controller u s .Here are the Equations ( 18)-( 21), as follows: where w is the estimated uncertainty of the external force, β 0 controls the time when the system reaches the sliding plane, and is a boundary layer that controls the tolerance in the sliding plane, which is the main issue of the sliding mode controller. is the vector multiplication operator.In practical applications of sliding mode controllers, the undesirable phenomenon of oscillations called "chattering" having finite frequency and amplitude may occur in the control inputs.Chattering leads to low control accuracy, high wear rate of mechanical systems, and high heat in electrical systems.To reduce the chattering phenomenon, the boundary layer , as shown in Figure 6, is set as smoothing forms of control inputs.As the system satisfies the condition .V < 0 from Equation (21), the system is globally uniformly asymptotically stable according to the LaSalle and Yoshizawa theorem [27][28][29].Khalil [30,31] described the details of the sliding mode controller.Figure 7 shows the block diagram of the sliding mode controller, illustrating its overall process.

Design of a Sliding Mode Controller
This section describes the contents of the DP simulation with the SM controller and their results.It also discusses the characteristics and validity of the SM controllers from the simulation results.

Configuration of Simulations
Implementing an SM controller requires a dynamic system model, uncertainty estimation, and external forces.In this study, an EKF was used to estimate uncertain external forces and apply them to SM controllers, as illustrated in Figure 8. Figure 9 shows a comparison between actual external disturbance F ext and estimated external disturbance Fext (blue and red lines, respectively).F i, ext (i = 1, 2, and 6) denotes external disturbances along the surge, sway, and yaw directions.The external disturbance error was calculated using Equation (22), and the values were 0.0622, 0.0487, and 0.0453 along the surge, sway, and yaw directions, respectively.The results show that the error is relatively small.It was assumed that the vessel model was built as a FPSO vessel.The FPSO performs turret connection work in the sea, and the SM controller adjusts the position of the turret and validates its performance.The length of the FPSO model was 300 m, and its main specifications and parameters are listed in Table 1.It was also assumed that six azimuth thrusters with a maximum thrust capacity of 500 kN were installed for DP control of the FPSO, and the information of the thruster is described in Table 2. Table 3 describes the ocean environmental conditions selected for the simulation.The control variables for the SM controller include the feedback control gain I λ , switching control gain F, and sliding plane boundary variable for chattering reduction.The control variables were set according to Table 4, considering the SM control variable proposed by Tannuri et al. [20].The reaching time t reach was set to 30 s when the sliding variable was located outside the sliding plane boundary.In addition, the maximum allowable error e max of surge, sway, and yaw were set as 2 m, 2 m, and 1 • , respectively.

Simulation Results and Discussions
Multiple scenarios were considered for the performance evaluation of the SM controller designed to set the position of the FPSO performing turret access operations as shown in Table 5.It was assumed in all the scenarios that the turret installed on the FPSO vessel reached the set position to perform the approach operation, initially maintaining the zero target angle (angle on the bow) that the hull holds parallel to the angle of incident waves and currents.Then, the controller changes the target angle.To help understand the control scenario, Figure 10 illustrates the control scenarios for the dynamic positioning of a vessel with environmental disturbance in the clockwise and counter-clockwise DP control.Simulations were conducted to evaluate the effects of different target angles, for 5000 s.
Table 5. Scenarios for dynamic positioning control of a FPSO.

Scenario Initial Target Final Target
Case To evaluate the performance of the SM controllers, a proportional-integral (PD) controller was designed with appropriate gains.Then, the performance of the SM controller was compared with that of the PD controller in each control scenario described in Table 5.The PD controller gains used in the simulation are presented in Table 6. Figure 11 shows the simulation results of the dynamic motions of the FPSO vessel and their relative errors concerning surge, sway, and yaw, with a comparison between the SM and PD controllers.The first comparison between the controllers is of clockwise control against external loads to a target angle of −20 • in case 1.At the simulation time of 2500 s, a control command was entered to change the target angle from 0 to −20 • .The transient region occurs according to this step input.The control output from the SM controller, which rapidly increases in the transition region and then smoothly converges to the control objective is shown in Figure 11a.The SM controller exhibited outstanding performance as the SM controller stabilized the motion of the FPSO vessel in steady-state conditions at approximately 3000 s without overshooting.This means that the SM controller has the ability to eliminate steady-state errors as shown in Figure 11b.By contrast, the result of the PD controller has a different pattern in that it has overshoot, deviations, and offset values and is significantly different from the SM controller.The relative positioning error of the PD controller could not converge to zero and kept outputting deviations and offset values from the control objectives.Figure 11b shows that the PD controller has relatively large errors in comparison with the SM controller.Figure 12 shows the results of the surge, sway, and yaw motions of the FPSO vessel and their relative errors with the SM and PD controllers according to the case of counterclockwise control against external loads at a target angle of 20 • in case 2. At a simulation time of 2500 s, a control command was entered to change the target angle from 0 to 20 • .The variation of motions and relative errors in the PD controller are significant, and the characteristics of the PD controller are clearly evident.The motions of the FPSO vessel and relative errors from the PD controller become unstable and diverse with large oscillations.By contrast, the SM controller can react while maintaining the desired positions and target angles.It seems that the SM controller can effectively eliminate steady-state errors.This is because the SM controller is a multi-input-multi-output (MIMO) system that performs target angle control and simultaneously minimizes interference with other motions such as surge, sway, and yaw.The PD controller is a single-input-single-output (SISO) system that considers all motions as independent control elements.Therefore, the SM controller performs better than the PD controller with ocean environmental disturbances from winds and waves at different target angles.
The control forces and moments along the x, y, and z directions for cases 1 and 2 are shown in Figure 13 to compare the results of the SM and PD controllers.The control force required during the operation of the FPSO turret connection is expressed as the instantaneous control power produced by the controller.In the control forces from the SM controller in case 1, it can be seen that there are no significant variations in the control forces in the results of the SM controller.This means that the SM controller produces a similar control force indicating that less environmental external forces act in the surge when the FPSO controls the target angle against the wave, current, and wind loads.By contrast, there are huge variations in the control forces from 2500 to 2800 s in the transient region from the PD controller in case 1.The system becomes unstable and diverse with a large oscillation in case 2, and the control forces of the PD controller became more necessary.Regarding the sliding plane, Figure 14 shows the sliding plane response of the SM controller, as defined in Equation (15).s 1 , s 2 , and s 6 represent the sliding variables in the surge, sway, and yaw directions, respectively, and the dotted line represents the values of the sliding plane boundary variables to mitigate the chattering of the SM controller.The sliding variables in surge and sway motions are located inside the sliding plane boundary at the target angles of −20 and 20 • (cases 1 and 2), which means that the sliding plane response is stable.The sliding variable in yaw moves out of the sliding plane boundary when the step input is applied to the system and returns within a certain period (approximately 23 s).The sharp variations appearing on the sliding surface level were due to the abrupt change in the set point.This result indicated that the sliding plane boundary reached time t reach in less than 30 s, as shown in Table 4.   Figure 16 shows the trajectories of the FPSO vessel with DP systems in cases 1 and 2 in comparison with the SM and PD controllers.As shown in Figure 16b, the PD controller could not maintain the angle and position against the wind and wave loads in case 2. These results show that the PD controller is highly sensitive to winds and waves.By contrast, the system output from the SM controller is much more stable, which ensures the robustness of the system to winds and waves, as shown in Figure 16a,b.Figures 17 and 18 show the root-mean-square (RMS) values of the relative errors of motions and control forces of the FPSO vessel in steady-state conditions between −35 and 25 • of the target angles.The RMS values were calculated for the responses to FPSO after steady-state conditions for 30 min.The results show that the PD controller reacts as expected when the target angle is from −20 to 14.2 • , but when the target angle is beyond a certain limit above 14.3 • and below −20.4 • , the RMS values begin substantially increase.With substantially increasing RMS values, the FPSO vessel loses stability and becomes unstable to converge the control objective.By contrast, the RMS values of the errors of motions with the SM controller were close to zero.It is observed that most of the RMS of the vessel motion errors and control forces from the SM controller are generally smaller than the values of the PD controllers.Therefore, satisfactory results were obtained for every target angle between −32 and 20 • as shown in Figure 19, in the response of the FPSO with the SM controller.On the other hand, the PD controller can adjust the target angle of the vessel between −20 and 14 • .This is another advantage of the SM controller, in addition to its robust characteristic and its applicability in complex marine operations and systems.

Concluding Remarks
This paper dealt with a sliding mode control algorithm, which is considered as a robust dynamic positioning control technique applicable to various tasks in the marine industry; its validity was evaluated through a number of simulations.MATLAB/Simulink was used to configure the simulator for the performance evaluation of the DP system with a sliding mode (SM) controller.The simulator consists of a controller module, vessel motion module, sensor module, and filter module.The SM controller was designed to dynamically maintain the turret position by performing the connection work of the FPSO in the ocean.The Lyapunov stability theory was applied at the design stage to ensure the stability of the control system.
Conventional PD controllers were configured and compared to evaluate the validity of the SM-controller-based DP system.SM controllers showed good performance, with only small errors, low required control force, and stability.This is because the SM controller can correct the steady-state error, adjust the control parameters considering the physical characteristics of the vessel, and perform control by considering the interference between surge, sway, and yaw.In particular, the proposed DP system includes a function to adjust the control performance through the selection of control parameters, which intuitively considers the response characteristics of vessels and is highly desirable in diverse and complex marine operations.As an extension of the study, the simulation results from the SM-controller-based DP system should be experimentally validated.In addition, there is scope for future work in updating the proposed DP system with the capabilities to predict incident sea loads from future occurrences of swell, gust, etc.; a DP system with such capabilities can improve the robustness against harsh ocean environments.

Figure 1 .
Figure 1.Block diagram of a dynamic positioning system of marine vessels.

Figure 2 .
Figure 2. Coordinate systems: global reference frame {e} and the body-fixed coordinate {b}.
. ξ(k), η(k), .η (k), Fext (k) T .The EKF calculates estimated external disturbance Fext induced by the ocean environment, and it can be applied to the SM controller to compensate the ocean environmental loads.Additionally, Q k and R k are the variance matrices of the process and measurement noises, respectively.

Figure 5 .
Figure 5. Flowchart of the sliding mode control algorithm.

Figure 6 .
Figure 6.The boundary layer for chattering reduction.

Figure 7 .
Figure 7. Block diagram of the sliding mode controller.

Figure 8 .
Figure 8. Structure of the sliding mode control with extended Kalman filter.

Figure 9 .
Figure 9.Comparison between actual external disturbance and estimated external disturbance calculated by EKF (blue: actual, red: estimated).

Figure 10 .
Figure 10.Control scenarios for dynamic positioning of a vessel with environmental disturbance: (a) clockwise and (b) counter-clockwise DP control.

Figure 11 .
Figure 11.Simulation results of the dynamic motions of FPSO vessel and relative errors with a comparison of the SM and PD controllers according to case 1 according to: (a) dynamic motions and (b) relative errors.

Figure 12 .
Figure 12.Simulation results of the dynamic motions of FPSO vessel and relative errors with a comparison of the SM and PD controllers according to case 2 according to: (a) dynamic motions and (b) relative errors.

Figure 13 .
Figure 13.Simulation results of the control forces of FPSO vessel with comparison of the SM and PD controllers according to: (a) case 1 and (b) case 2.

Figure 14 .
Figure 14.Responses of the sliding surface of the SM controller depending on the control scenario: (a) case 1 and (b) case 2.

Figure 15
Figure15represents the yaw response of the phase plane that simultaneously shows the error e comprising the sliding variable s and the response of the derivative of the error .e in cases 1 and 2, respectively.The process of convergence to zero is shown as errors, and their derivatives are simultaneously reduced from the results.The blue dotted line represents the sliding boundary variable for mitigating the chattering of the SM controller.The trajectories in phase planes converged to the sliding surface, and the controller maintained the trajectories inside the boundary layer.These figures indicate that the system is stable because the initial and final sliding variables in the yaw motion are located inside the sliding plane boundaries.

Figure 15 .
Figure 15.Yaw responses of the phase plane of the SM controller: (a) case 1 and (b) case 2 (x: initial sliding variable, o: final sliding variable, blue dotted line: boundary layer).

Figure 16 .
Figure 16.Trajectories of a vessel with dynamic positioning systems with environmental disturbance in comparison of sliding mode (SM) and proportional-integral (PD) controllers: (a) case 1 and (b) case 2.

Figure 17 .
Figure 17.RMS values of errors of the motions of the FPSO vessel in steady state condition: (a) Target angle between −35 to 25 • and (b) Target angle between −20 to 14 • .

Figure 18 .
Figure 18.RMS values of the control forces of the FPSO vessel in steady-state condition: (a) Target angle between −35 to 25 • and (b) Target angle between −20 to 14 • .

Figure 19 .
Figure 19.Possible controllable range of DP systems: (blue) SM and (red) PD controllers.

Table 1 .
Vessel Specifications and parameters of FPSO.

Table 2 .
Specifications of thruster systems for FPSO.

Table 3 .
Environmental condition in the simulation.

Table 4 .
Control parameters for the sliding mode control.