Next Article in Journal
Experimental and Numerical Studies on the Ultimate Bending Strength of Welded Plated Grillage with Combined Openings
Previous Article in Journal
Experiments on Cavitation Control around a Cylinder Using Biomimetic Riblets
Previous Article in Special Issue
Cooperative Maritime Search of Multi-Ship Based on Improved Robust Line-of-Sight Guidance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Model Predictive Control for a Dynamic Positioning Ship Based on the Laguerre Function

1
Zhilong (Dalian) Marine Technology Co., Ltd., 11th Floor, No. 523 Huangpu Road, Dalian 116020, China
2
Key Laboratory of Marine Dynamic Simulation and Control, Department of Navigation, Dalian Maritime University, Dalian 116026, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(2), 294; https://doi.org/10.3390/jmse12020294
Submission received: 28 December 2023 / Revised: 21 January 2024 / Accepted: 4 February 2024 / Published: 7 February 2024
(This article belongs to the Special Issue Optimal Maneuvering and Control of Ships)

Abstract

:
In this paper, we present a novel nonlinear model predictive control (NMPC) algorithm based on the Laguerre function for dynamic positioning ships to solve the problems of input saturation, unknown time-varying disturbances, and heavy computation. The nonlinear model of a dynamic positioning ship is presented as a linear model, transformed from a standard affine nonlinear state-space model by precise feedback linearization. The environmental disturbance is overcome using an integrator. The time cost of the proposed nonlinear control algorithm is decreased by inducing the Laguerre function to describe the feedback-linearization system input increments. The Laguerre function reduces the matrix dimensions of the nonlinear optimization problem. The simulation results for a DP supply vessel showed that the novel algorithm maintained the effective control performance of the original nonlinear model predictive control algorithm and had a reduced computation load to satisfy the requirements of real-time operation.

1. Introduction

According to definitions provided by the International Maritime Organization (IMO) and Det Norske Veritas (DNV), a dynamic positioning system (DPS) comprises all the necessary equipment for a dynamic positioning ship, including the power, propulsion, and dynamic positioning control systems [1,2]. Dynamic positioning systems have been widely used in marine engineering fields such as offshore oil and gas resource development and submarine pipeline laying. Shatto [3] created the first dynamic positioning system in 1961, equipped with four steering propellers. The first DPS was employed by the ship “Cuss 1” in the California Sea. In the same year, the drillship “Eureka” [3], belonging to the Royal Dutch Shell company, was equipped with the first analog-signal DP system. The first digital-signal DP ship, “Gloma Challenger” [3], was built in 1968 and traveled almost every ocean on Earth, providing a wealth of favorable evidence for geological discoveries, especially for the theory of plate and shell structures. Subsequently, the application of microwave position reference, underwater acoustic position reference, satellite positioning, and other position reference systems improved the accuracy of dynamic positioning systems.
The dynamic position control algorithm is the key component of a DPS. The first generation of dynamic positioning systems used conventional proportional–integral–derivative (PID) control laws. To avoid responding to high-frequency movements, low-pass or notch filters were employed to eliminate high-frequency components from the deviation signal [4].
In the mid-1970s, Balchen and others proposed the use of optimal control and Kalman filtering theory combined with dynamic positioning control, giving rise to the second generation of dynamic positioning systems, which were also widely applied [5,6,7]. Currently, the third generation of dynamic positioning systems focuses on advanced control algorithms such as nonlinear control algorithms and intelligent control theories [8,9,10,11]. Tannuri [12] designed a control algorithm for dynamic positioning systems based on sliding mode control and employed nonlinear multi-variable mathematical models, providing robustness to variations in displacement and environmental conditions. Do [13] developed a globally robust and adaptive output feedback controller for the dynamic positioning of surface ships under environmental disturbances based on Lyapunov’s direct method, forcing the ship’s position and orientation to globally asymptotically converge to the desired values. Hassani [14] proposed a new strategy for the design of robust DP controllers for marine vessels under different sea conditions using mixed- μ synthesis. Yang [15] proposed a robust adaptive NN-based output feedback control scheme for a dynamic positioning ship with uncertainties and unknown external disturbances. Zhang [16] developed a dynamic event-triggered mechanism for dynamic positioning vehicles with input saturation suited to marine applications due to its concision and flexibility. Cho [17] presented a sliding mode control algorithm as a robust dynamic positioning control technique applicable to various tasks in the marine industry.
Due to the influence of the mechanical performance of the thrusters, the control force and torque generated by a dynamic positioning system are subject to amplitude and increment constraints, i.e., the system inputs and their rates of change present saturation issues. The model predictive control (MPC) algorithm is the only advanced technology that can handle the online optimization control of multi-variable constrained systems in a systematic and intuitive manner [18], achieving favorable performance and robustness. Wang Yuanhui et al. [19], drawing inspiration from GreenDP, designed a dynamic positioning controller combining the Kalman filter and model predictive control. Subsequently [20], they used correlation-based non-switching analytical model predictive control theory to design a nonlinear model predictive controller for dynamic positioning, addressing the nonlinear control problem of dynamic positioning ship motion. Veksler [21] adopted the MPC algorithm to combine positioning control and thrust allocation into a single algorithm that could theoretically yield a near-optimal controller output. Miller [22] used the MPC regulator to control the LNG carrier training service ship “Dorchester Lady” and proved that a predictive controller could be built to steer an SS during an UNREP maneuver. However, its application was difficult due to the requirement of a ship dynamics’ linear incremental model. Based on this, Miller [23] combined the MPC algorithm with the line-of-sight (LOS) approach to control Maritime Autonomous Surface Ships (MASSs) through the ship trajectory tracking system. This was combined with a variable maneuver path advance, leading to effective trajectory tracking on turns and built-in integral action reference correlation.However, the above methods only considered the thrust amplitude constraints generated by the thrusters, ignoring the constraints on the thrust amplitude increment within a certain time and non-zero-mean disturbance issues.
Model predictive control (MPC) systematically and intuitively determines the current optimal control actions based on given constraints and performance requirements. However, it faces the challenge of substantial online computations. Liuping Wang [24] proposed a discrete model predictive control method using the Laguerre function to describe the control increment signal, thereby reducing the computational load of the algorithm. Kong Xiaobing et al. [25] introduced a nonlinear model predictive control algorithm based on input–output feedback linearization, employing an approximate optimization method to reduce the online computation load of solving nonlinear constrained problems in sequential quadratic programming during the rolling optimization process.
To address the nonlinear nature of dynamic positioning ship motion, the saturation of input amplitudes and their rates of change, and the susceptibility to disturbances from the marine environment, a novel nonlinear model predictive control algorithm for the dynamic positioning of ships is proposed based on the Laguerre function. Compared to the original algorithm, this new approach effectively reduces the online computational load while retaining control performance. It allows a ship to rapidly reach and maintain the desired position in real time.
This paper is organized as follows: Section 2 describes the problem and provides background knowledge. Section 3 presents the design process of the dynamic positioning controller and an analysis of the computation load. In Section 4, we provide simulation results demonstrating the effective performance of the proposed DP controller. Section 5 concludes the paper.

2. Problem Description and Preliminaries

This section describes the nonlinear control problem of a DP ship and introduces several key lemmas for developing the formation control strategy.

2.1. Problem Statement

A motion model of a dynamic positioning vessel at sea can be divided into low-frequency and high-frequency motions. The former describes the low-frequency motion generated by external environmental disturbances, such as low-frequency wind, second-order wave forces, and ocean currents. The latter is associated with the high-frequency motion of the vessel caused by high-frequency wind and first-order wave forces. To reduce fuel consumption and propulsion equipment wear, typically only the low-frequency motion of dynamic positioning vessels is considered, focusing on the control of horizontal motion in three degrees of freedom: surge, sway, and yaw.
In Figure 1, vector η = [ x , y , ψ ] T represents the position and heading of the vessel in the north–east coordinate system O 0 X 0 Y 0 , and vector ν = [ u , v , r ] T represents the horizontal velocity and yaw angular velocity of the vessel in the body-fixed coordinate system. Referring to Fossen, the mathematical model for the low-frequency motion of a dynamic positioning vessel on the sea surface is expressed as [2]
η ˙ = R ( ψ ) ν M ν ˙ + D ( ν ) ν = τ + τ E n v .
In this equation, M is the inertia matrix incorporating additional mass and is reversible and positive definite. D is the linear damping matrix, also positive definite. τ represents the actual thrust generated by the thrusters, and τ E n v is an unknown time-varying low-frequency disturbance term, representing the disturbing forces and moments from slowly varying environmental disturbances such as wind, waves, and currents. R is the rotation matrix [2], expressed as
M = m 11 0 0 0 m 22 m 23 0 m 32 m 33
R ( ψ ) = cos ψ sin ψ 0 sin ψ cos ψ 0 0 0 1
D = d 11 0 0 0 d 22 d 23 0 d 32 d 33 .
Suppose that x = [ x , y , ψ , u , v , r ] T , u = τ , d = τ E n v , and
α 11 = d 11 / m 11 α 22 = ( d 22 m 33 d 32 m 23 ) / ( m 22 m 33 m 23 m 32 ) α 23 = ( d 23 m 33 d 33 m 23 ) / ( m 22 m 33 m 23 m 32 ) α 32 = ( d 22 m 32 d 32 m 22 ) / ( m 22 m 33 m 23 m 32 ) α 33 = ( d 23 m 32 d 33 m 22 ) / ( m 22 m 33 m 23 m 32 )
β 11 = 2 / m 11 β 22 = m 33 / ( m 22 m 33 m 23 m 32 ) β 23 = m 23 / ( m 22 m 33 m 23 m 32 ) β 32 = m 32 / ( m 22 m 33 m 23 m 32 ) β 33 = m 22 / ( m 22 m 33 m 23 m 32 ) .
Thus, the standard nonlinear state-space model for dynamic positioning vessels is obtained as follows:
x ˙ = f ( x ) + g ( x ) ( u + d ) y = h ( x ) ,
where
f ( x ) = x 4 c o s ( x 3 ) x 5 s i n ( x 3 ) x 4 s i n ( x 3 ) + x 5 c o s ( x 3 ) x 6 α 11 x 4 α 22 x 5 + α 23 x 6 α 32 x 5 + α 33 x 6
g ( x ) = 0 0 0 0 0 0 0 0 0 β 11 0 0 0 β 22 β 23 0 β 32 β 33
h ( x ) = h 1 ( x ) h 2 ( x ) h 3 ( x ) = x 1 x 2 x 3 .
To address the nonlinear dynamic positioning ship motion model described by Equations (1) and (7), we designed a nonlinear model predictive controller (NMPC). The goal was to enable the dynamic positioning ship to overcome environmental disturbances and reach and maintain the desired position while satisfying the system constraints.

2.2. Differential Geometry

For the control of nonlinear systems, researchers often neglect higher-order nonlinear terms by expanding the system model to a first-order linear term through Taylor expansion. However, precise feedback linearization technology uses differential geometry methods for linearization without ignoring any nonlinear terms, and the linearization of this method is accurate. To apply precise feedback linearization technology to nonlinear model predictive control, the following definition needs to be introduced.
Definition 1. 
If f is a smooth vector field on U and h is a smooth scalar function on U, then f ( h ) is a smooth function on U, defined as
L f h ( x ) = i = 1 n f i ( x ) ( h ( x ) x i )
The function L f h is the Lie derivative of function h along the vector field f.
Definition 2. 
For a multiple input multiple output (MIMO) affine nonlinear system,
x ˙ = f ( x ) + g ( x ) u y = h ( x ) ,
where x = R n is the n order state variables of the system. f , g , and h are the vectors field that are sufficiently smooth on the domain D R n . Mapping f: D R n and g: D R n are the vector fields in D.
For the MIMO system, there exists the integers ( r 1 , r 2 , , r m ) in the neighborhood of x 0 , which satisfies that
L g i L f k h ( x ) = 0 1 j m 1 i m k r i 1 ,
α ( x ) = L g 1 L f r 1 1 h 1 ( x ) L g 2 L f r 1 1 h 1 ( x ) L g m L f r 1 1 h 1 ( x ) L g 1 L f r 2 1 h 2 ( x ) L g 2 L f r 2 1 h 2 ( x ) L g m L f r 2 1 h 2 ( x ) L g 1 L f r m 1 h m ( x ) L g m L f r m 1 h m ( x ) L g m L f r m 1 h m ( x ) .
If α ( x ) is the nonsingular matrix when x = x 0 , the integers ( r 1 , r 2 , , r m ) are the relative degree of the output variable y i . The total relative degree is the sum, r = r 1 + r 2 + + r m . If the total relative degree r = n , the system can achieve the precise feedback linearization.

2.3. Laguerre Function

The Laguerre functions are a class of orthogonal functions primarily used for system identification. The discrete Laguerre network is shown in Figure 2 and expressed in Equation (15) [24].
Γ k ( z ) = Γ k 1 ( z ) z 1 a 1 a z 1 , Γ 1 ( z ) = 1 a 2 1 a z 1 ,
where a represents the poles of the discrete-time Laguerre network; the stability region is 0 a 1 ; and N is the order of the Laguerre network. By inverse-z-transforming each term in Equation (15), we obtain the discrete-time Laguerre functions, represented in vector form as
L ( k ) = l 1 ( k ) l 2 ( k ) l 3 ( k ) l N ( k ) .
The relationship between adjacent discrete-time Laguerre functions at consecutive instants can be derived from Equations (15) and (16):
L ( k + 1 ) = A l L ( k ) ,
where the matrix A l is a function of parameters a and β , with initial conditions
L ( 0 ) = β 1 a ( a ) N 1 T
A l = a 0 0 0 β a 0 0 a β β a 0 ( a ) N 2 β ( a ) N 3 β ( a ) N 4 β a .
Additionally, an important property of discrete-time Laguerre functions is orthogonality [24]:
k = 0 l i ( k ) l j ( k ) = 0 , i j k = 0 l i ( k ) l j ( k ) = 1 , i = j .
Based on the above properties, the impulse response of a stable system can be represented by the following discrete-time Laguerre functions:
H ( k ) = L ( k ) T η ¯ = L ( k ) T c 1 c 2 c 3 c N T ,
where η ¯ = L ( k ) T [ c 1 , c 2 , c 3 , , c N ] are coefficients determined by the system data, and leveraging the orthogonality of the Laguerre functions allows one to obtain the values of each coefficient.
c i = k = 0 H ( k ) l i ( k )

3. Steps for Dynamic Positioning Controller Design

3.1. Precise Feedback Linearization

The mathematical model of dynamic positioning vessel motion in Equation (7) is a standard affine nonlinear expression. Feedback linearization is required before applying the model predictive control algorithm. First, the Lie derivatives of the system output variables in Equation (7) are calculated:
L g 1 L f 0 h 1 L g 2 L f 0 h 1 L g 3 L f 0 h 1 L g 1 L f 0 h 2 L g 2 L f 0 h 2 L g 3 L f 0 h 2 L g 1 L f 0 h 3 L g 2 L f 0 h 3 L g 3 L f 0 h 3 = 0 0 0 0 0 0 0 0 0
a ( x ) = L g 1 L f 0 h 1 L g 2 L f 0 h 1 L g 3 L f 0 h 1 L g 1 L f 0 h 2 L g 2 L f 0 h 2 L g 3 L f 0 h 2 L g 1 L f 0 h 3 L g 2 L f 0 h 3 L g 3 L f 0 h 3 = β 11 cos x 3 β 22 sin x 3 β 23 sin x 3 β 11 sin x 3 β 22 cos x 3 β 23 cos x 3 0 β 32 β 33 .
According to Equations (23) and (24), the relative orders of the system are r 1 = 2 , r 2 = 2 , and r 3 = 2 , and the total relative order r = r 1 + r 2 + r 3 is equal to the number of system state variables. Therefore, the nonlinear model of dynamic positioning vessel motion can undergo exact state feedback linearization, resulting in new state variables:
x m 1 = L f 0 h 1 x = x 1 x m 2 = L f 0 h 2 x = x 2 x m 3 = L f 0 h 3 x = x 3 x m 4 = L f 1 h 1 x = x 4 cos ( x 3 ) x 5 sin ( x 3 ) x m 5 = L f 1 h 2 x = x 4 sin ( x 3 ) + x 5 cos ( x 3 ) x m 6 = L f 1 h 3 x = x 6 .
The nonlinear feedback control rate is given by
u = a ( x ) 1 ( b ( x ) + v ) d ,
in which
b ( x ) = L f 2 h 1 L f 2 h 2 L f 2 h 3 = cos x 3 ( α 11 x 4 x 5 x 6 ) sin x 3 ( α 22 x 5 + α 23 x 6 + x 4 x 6 ) sin x 3 ( α 11 x 4 x 5 x 6 ) + cos x 3 ( α 22 x 5 + α 23 x 6 + x 4 x 6 ) α 32 x 5 + α 33 x 6 .
The linear model obtained through feedback linearization is
x ˙ m = A c x m + B c v + E c d y = C c x m .
The coefficients of each matrix in this equation are A c = 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , B c = 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 1   C c = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 ,   E c = B c a ( x ) .

3.2. Linear Control Structure

Discretizing the state-space model (28) with a step size of h yields the discretized linear model:
x m ( k + 1 ) = A d x m ( k ) + B d v ( k ) + E d d ( k ) y ( k ) = C d x k .
Here, x m ( k ) represents the state of the linear system at time k; v ( k ) is the input to the linear system; y ( k ) is the output of the linear system at time k; and A d , B d , C d , and E d are the constant matrices after discretization.
Suppose that x m ( k + 1 ) = x m ( k ) + Δ x m ( k + 1 ) and v m ( k + 1 ) = v m ( k ) + Δ v m ( k + 1 ) ; then, based on Equation (29), we obtain
Δ x m ( k + 1 ) = A d Δ x m ( k ) + B d Δ v ( k ) y ( k + 1 ) = y ( k ) + C d ( A d Δ x m ( k ) + B d Δ v ( k ) ) .
Suppose that x u ( k ) = [ Δ x m ( k ) T , y ( k ) T ] T ; introducing an integral component, we obtain an augmented state-space model that eliminates the slowly varying disturbance term d, expressed as follows:
x u ( k + 1 ) = A u ( k ) + B Δ v ( k ) y ( k ) = C d x u ( k ) .
Here, A , B , and C represent the constant augmented matrix coefficients, which are as follows: A = A d 0 6 × 3 C d A d I 6 × 3 , B = B d C d B d , C = 0 3 × 6 I 3 × 3 .
According to Equations (29) and (30), the future N c system outputs N p in the control time domain can be predicted as follows ( N c N p ) :
Y = F x u ( k ) + Φ Δ V .
In the equation above, Y = y ( k + 1 | k ) y ( k + 1 | k ) y ( k + 1 | k ) ,   Δ V = Δ v ( k ) Δ v ( k + 1 ) Δ v ( k + N c 1 )   F = C A C A 2 C A N p , Φ = C B 0 0 C A B C B 0 C A N p 1 B C A N p 2 B C A N p N c B .
Additionally, based on Equations (26), (29), and (30), the relationship between the entire predicted time domain’s linear system state variables and the input v can be obtained:
X m k = Γ + Π Δ V ,
where X m ( k ) = x m ( k ) x m ( k + 1 ) x m ( k + N c 1 ) , Γ = x m ( k ) x m ( k ) + i = 1 1 A d i Δ x m ( k ) x m ( k ) + i = 1 N c 1 A d i Δ x m ( k )   Π = 0 0 0 i = 1 1 ( A m i 1 ) B m 0 0 i = 1 N p 1 ( A m i 1 ) B m i = 1 N p 1 ( A m i 1 ) B m 0 .
The optimization objective function is formulated as a quadratic function of Δ V :
J = ( N d Y ) T ( N d Y ) + Δ V T R w Δ V ,
where R w represents the input weight. Substituting Equation (32) into the previous expression yields the performance index J m i n with respect to Δ V :
J m i n = 1 2 Δ V T H Δ V + f T Δ V ,
where the Hessian matrix H and the vector f are defined as follows:
H = 2 ( Φ T Φ + R w )
f = 2 Φ ( N d F x ( k ) ) .
Under unconstrained conditions, the optimal solution for Δ V is given by
Δ V = ( Φ T Φ + R w ) 1 Φ T ( N d F x u ( k ) ) .
Due to the physical limitations of the propulsion system on a dynamic positioning vessel, the total control force and moment generated by the thrusters are constrained. Additionally, the amplitude of the total control force and moment changes within a certain time is also constrained, indicating the presence of input and input rate saturation in the system:
u m i n ( k ) u ( k ) u m a x ( k )
Δ u m i n ( k ) Δ u ( k ) Δ u m a x ( k ) .
The input u of the nonlinear dynamic positioning vessel motion model under the mapping of the nonlinear feedback control rate (26) results in the input v for the new linear system. The linear inequality constraints (39) and (40) with respect to u are transformed into nonlinear inequality constraints for v . However, the upper and lower bounds of V throughout the control time domain, V m i n and V m a x , cannot be directly determined. Their values are dependent on the system’s state at each sampling moment and need to be determined. According to Equation (26), the relationship between U and V throughout the entire control time domain can be obtained:
u ( k ) = a ( x ( k ) ) 1 ( b ( x ( k ) ) + v ( k ) ) u ( k + 1 ) = a ( x ( k + 1 ) ) 1 ( b ( x ( k + 1 ) ) + v ( k + 1 ) ) u ( k + N c 1 ) = a ( x ( k + N c 1 ) ) 1 ( b ( x ( k + N c 1 ) ) + v ( k + N c 1 ) ) .
That is,
U ( k ) = f x v ( X ( k ) , V ( k ) ) .
The original nonlinear system input u and the linear system input v can both be obtained by summing their control increments u and v ; that is,
u ( k + m ) = u ( k 1 ) + Δ u ( k ) + + Δ u ( k + m )
v ( k + m ) = v ( k 1 ) + Δ v ( k ) + + Δ v ( k + m ) .
Therefore, U and V over the entire control time domain are given by
U ( k ) = C 1 u ( k 1 ) + C 2 Δ U ( k )
V ( k ) = C 1 v ( k 1 ) + C 2 Δ V ( k ) ,
where C 1 = I I I , C 2 = I 0 0 I I 0 I I I .
Substituting the combined expressions (33) and (46) into (42) yields a nonlinear expression for U in terms of Δ V :
U = f Δ v ( Δ V ) .
The constraint conditions, transformed from the linear constraints on U in (39) and (40), are combined with (36) and (47) to form the following nonlinear programming problem:
m i n J = 1 2 Δ V T H Δ V + f T Δ V s u c h t h a t c ( Δ V ) 0 .
The left-hand side of the nonlinear constraint inequality in this equation is
c ( Δ V ) = f Δ v ( Δ V ) U m a x U m i n f Δ v ( Δ V ) C 2 1 ( f Δ v ( Δ V ) C 1 u ( k 1 ) ) Δ U m a x Δ U m i n C 2 1 ( f Δ v ( Δ V ) C 1 u ( k 1 ) ) .
The current optimal solution for the feedback-linearized system input at time k, denoted as Δ V o p t , is obtained by solving the nonlinear programming problem in Equation (48) using the sequential quadratic programming method. In this process, the initial point for iteration is set according to Equation (38), and the transformed values are derived to obtain the inputs for the original nonlinear system, representing the required thrust and moments for the dynamic positioning of the ship.
u ( k ) = a ( x ( k ) ) 1 [ a ( x ( k 1 ) ) u ( k 1 ) + b ( x ( k 1 ) ) b ( x ( k ) ) + Δ v ( k ) ]

3.3. Introduction to Laguerre Functions

A nonlinear model predictive control algorithm typically selects a relatively large and appropriate control horizon ( N c ) to ensure a more favorable dynamic response and stability. However, increasing the control horizon also increases the solution time for the nonlinear programming problem, preventing the computer from issuing control commands in real time and affecting the positioning of the ship. Therefore, we introduce Laguerre functions to describe the control increments of the linearized system after feedback, proposing a new low-computational-cost nonlinear model predictive control algorithm for ship dynamic positioning.
We assume that the future m time steps of the linear system control input increments at the current time step k are represented by the following Laguerre function:
Δ v ( k + m ) = L ( m ) T η ¯ k ,
where L ( m ) T = L 1 ( m ) T 0 1 × N 0 1 × N 0 1 × N L 2 ( m ) T 0 1 × N 0 1 × N 0 1 × N L 3 ( m ) T , η ¯ k = η ¯ k 1 η ¯ k 2 η ¯ k 3 .
Substituting Equation (51) into the augmented state-space model (31), we obtain a state-space model incorporating the Lagrange function:
x u ( k + 1 ) = A x u ( k ) + B L ( 0 ) T η ¯ k y ( k ) = C x ( k ) .
Thus, based on the current state and output sampled at time k, we can predict the system’s state and output at future time k + m :
x u ( k + m | k ) = A m x u ( k ) + ϕ ( m ) T η ¯ k
y ( k + m | k ) = C A m x u ( k ) + ϕ ( m ) 6 T η ¯ k ,
where ϕ ( m ) T = i = 0 m 1 A m i 1 B L ( i ) T . Then, we substitute the Laguerre function into the optimization objective function (Equation (34)) to obtain
J = m = 1 N p [ η d + y ( k + m ) ] T [ η d + y ( k + m ) ] + n = 0 N p 1 [ L ( n ) T η ¯ k ] T [ L ( n ) T η ¯ k ] .
Due to the sufficiently large prediction horizon N p , the orthogonality of the Lagrange functions can be exploited to simplify the second term of Equation (55), resulting in
J = m = 1 N p [ η d + y ( k + m ) ] T [ η d + y ( k + m ) ] + η ¯ k T R L η ¯ k .
Introducing the variable x f ( k ) = [ Δ m ( k ) T , e ( k ) T ] T and substituting Equations (53) and (54) into Equation (56) yields the new performance index function:
J = 1 2 η ¯ T H ¯ η ¯ + f ¯ T η ¯ ,
where the Hessian matrix H ¯ and the vector f ¯ are represented as
H ¯ = 2 ( m = 1 N p ϕ ( m ) C T C ϕ ( m ) T + R L )
f ¯ = 2 m = 1 N p ( ϕ ( m ) C T C A m ) x f ( k ) .
Under unconstrained conditions, the optimal solution for η ¯ is
η ¯ o p t = H ¯ f ¯ .
To handle the nonlinear constraint conditions in the original nonlinear programming problem (48), the Laguerre function must be sequentially substituted into (33), (42) and (46) to obtain the expression of U in terms of η ¯ :
U = f η ( η ¯ ) .
By combining Equation (57) and Equation (61), we rewrite the original nonlinear programming problem (48) in terms of Δ V , forming a new nonlinear programming problem with respect to η ¯ :
m i n J = 1 2 η ¯ T H ¯ η ¯ + f ¯ η ¯ s u c h t h a t c ( η ¯ ) 0 ,
where
c ( η ¯ ) = f η ¯ ( η ¯ ) U m a x U m i n f η ¯ ( η ¯ ) C 2 1 ( f η ¯ ( η ¯ ) C 1 u ( k 1 ) ) Δ U m a x Δ U m i n C 2 1 ( f η ¯ ( Δ V ) C 1 u ( k 1 ) ) .
Finally, employing the sequential quadratic programming method, the optimal solution to Equation (62) is determined, with Equation (60) used as the initial point. This provides the required thrust and torque for the current time step k in the dynamic positioning of the vessel.
u ( k ) = a ( x ( k ) ) 1 [ a ( x ( k 1 ) ) u ( k 1 ) + b ( x ( k 1 ) ) b ( x ( k ) ) + L ( 0 ) T η ¯ k o p t ]

3.4. Computational Analysis

The sequential quadratic programming algorithm involves approximating the nonlinear programming function at a certain point through a Taylor expansion into multiple quadratic programming problems and iteratively obtaining the optimal solution. Its computational complexity increases exponentially with the dimensionality of the variables.
The original nonlinear model predictive control method’s nonlinear optimization problem, expressed in Equation (48), and the new nonlinear optimization problem of the model predictive control algorithm based on the Laguerre function, expressed in Equation (62), have a similar form. However, the dimensions of their internal variables differ. In Equation (48), the variable Δ V has dimensions of 3 N C × 1 , the Hessian matrix H has dimensions of 3 N c × 3 N c , and the coefficient vector f has dimensions of 3 N c × 1 . Meanwhile, for the new algorithm’s optimization problem (62), the variable η ¯ has dimensions of 3 N × 1 , the Hessian matrix H ¯ has dimensions of 3 N × 3 N , and the coefficient vector F ¯ has dimensions of 3 N × 1 . In practical applications, the Laguerre function series N is much smaller than the control time domain N c . Therefore, in the process of sequential quadratic programming, with the same number of iterations, the computational load of the proposed ship dynamic positioning nonlinear model predictive control algorithm based on the Laguerre function is reduced compared to the original algorithm.

4. Simulation Verification and Data Analysis

To validate the effectiveness of the proposed ship dynamic positioning nonlinear model predictive control algorithm based on the Laguerre function, simulation verification was conducted using a dynamic positioning supply ship as the research object. The supply ship had a length of 76.1 m, a width of 18.8 m, a draft of 6.25 m, and a displacement of 4200 t. The main parameters are detailed in Table 1, and the non-dimensional inertia matrix and damping matrix were as follows [26]:
M = 1.1274 0.0000 0.0000 0.0000 1.8902 0.0744 0.0000 0.0744 0.1278
D = 0.0358 0.0000 0.0000 0.0000 0.1183 0.0124 0.0000 0.0041 0.0308 .
The simulation experiment was configured using a nonlinear model predictive controller with a prediction horizon ( N p ) of 10, a control horizon ( N c ) of 150, Laguerre function poles a = [ 0.46 , 0.46 , 0.46 ] T , and a series order of N = [ 4 , 4 , 4 ] T . The initial position of the ship was set to η 0 = [ 0 m , 0 m , 30 ] , and the desired position was η 0 = [ 20 m , 20 m , 40 ] . The low-frequency environmental disturbance was set to a certain value.
τ Env = R T ( ψ ) · 3 ( 1.2 + 1.5 sin ( 0.02 t ) + 1.4 sin ( 0.1 t ) ) kN 14 ( 0.9 + 2.0 sin ( 0.02 t ) + 1.3 sin ( 0.2 t ) ) kN 250 ( sin ( 0.09 t ) + 4 sin ( 0.01 t ) ) kN · m
Under simulated conditions equivalent to sea state four operations, the simulation time was set to 500 s with a time step of 0.5 s. The experiment was conducted on a Lenovo desktop computer equipped with an Intel Core(TM) i7-6700 CPU running at 3.4 GHz and with 8 GB of RAM. The results of the experiment are depicted in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, where the blue dashed line represents the original algorithm, and the red solid line represents the new NMPC algorithm based on the Laguerre function. Additionally, the results are summarized in Table 2.
Simulation experiment results analysis: Overall, the response curves of the proposed nonlinear model predictive control (NMPC) algorithm for ship dynamic positioning were almost identical to those of the original algorithm. Figure 3 and Figure 4 show the ship’s trajectory and time response under low-frequency environmental disturbance, indicating that the designed controller overcame ocean environmental disturbances to bring the ship to and maintain it at the desired position. Figure 5 represents the ship’s speed response curve; the speed components had relatively small values, and, under the influence of environmental disturbances, they oscillated near the zero point. Figure 6 illustrates the response curves of ship control forces and torques and their rate of change, showing that the control forces and torques were within the constraint range. Ultimately, they oscillated around a fixed value to counteract the impact of ocean environmental disturbances on the ship. Figure 7 represents the response curves for the rate of change in the control forces and torques, all within their constraint range, which oscillated around a zero mean. Figure 8 and Table 2 compare the computational load of the newly designed ship dynamic positioning NMPC algorithm based on the Laguerre function with that of the original algorithm. The graph indicates that the new algorithm’s single-computation time was significantly lower than that of the original algorithm. During the 1 s–17 s simulation period, the original algorithm’s single minimum computation time was 227.9 ms, with a maximum of 1847.0 ms and an average of 684.2 ms, preventing the real-time calculation of the control algorithm. During this period, the new algorithm’s single minimum computation time was 10.9 ms, the maximum was 227.9 ms, and the average was 134.6 ms, marking a decrease of 80.3% compared to the original algorithm. In the subsequent simulation period from 17 s to 500 s, as the sequential quadratic programming (SQP) iteration’s starting point was already the optimal point, the computational load was the same as when unconstrained. The new algorithm’s average computation time during this period was 13.6 ms, a 22.8% decrease compared to the original algorithm’s time of 10.5 ms. Overall, the average computation time of the new algorithm was 14.6 ms, representing a 59.1% decrease compared to the original algorithm’s time of 35.7 ms, ensuring the real-time performance of the algorithm solution.
The experimental results indicated that the ship dynamic positioning nonlinear model predictive control (NMPC) algorithm based on the Laguerre function designed in this work could overcome the impact of unknown time-varying ocean environmental disturbances under the conditions of input and input rate saturation. It enabled the ship to reach and maintain the desired position while retaining the control performance of the original nonlinear model predictive control algorithm. Additionally, the computational load significantly decreased, ensuring the real-time solution of the algorithm.

5. Conclusions

This paper proposes a novel nonlinear model predictive control (NMPC) algorithm for ship dynamic positioning based on the Laguerre function. The algorithm uses the Laguerre function to describe the control increment signal of the linear system after feedback linearization, reducing the dimensions of the coefficient matrices in the nonlinear constraint problems. This addresses the issue of the high computational load in the original nonlinear model predictive control algorithm. Finally, the algorithm was validated through simulation experiments on a supply ship. The results demonstrated that the improved ship dynamic positioning NMPC algorithm, under input and input rate saturation conditions, overcame the impact of unknown time-varying disturbances, such as wind, waves, and currents. It ensured that the ship reached and maintained the desired position and heading angle. The new algorithm retains the favorable control performance of the original NMPC algorithm while addressing the computational load issue, satisfying real-time computational requirements.
In future work, we will study the stability and reliability of the NMPC algorithm based on the Laguerre function and consider thrust allocation to solve the overdriving problem for dynamic positioning ships.

Author Contributions

Conceptualization, Y.Y.; methodology, Y.Y.; formal analysis, D.G. and X.Q.; software, X.Q. and H.S.; data curation, X.Q.; resources, X.Q.; writing—original draft preparation, X.Q.; writing—review and editing, X.Q. and D.G.; project administration, Y.Y.; funding acquisition, X.Q. and Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the following: (a) the National Key R&D Program of China (No. 2022YFB4300803, No. 2022YFB4301402) by Ministry of Science and Technology of the People’s Republic of China; (b) the 2023 Liaoning Provincial Science and Technology Plan (Key) Project (No. 2023JH1/10400047) by Liaoning Provincial Department of science and technology; and (c) the 2023 Dalian Science and Technology Talent Innovation Support Policy Project (No. 2023RY005) by Dalian Bureau of Science and Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

Author Xiaobin Qian was employed by the company Zhilong (Dalian) Marine Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Sørensen, A.J. A survey of dynamic positioning control systems. Annu. Rev. Control 2011, 35, 123–136. [Google Scholar] [CrossRef]
  2. Fossen, T.I. Handbook of Marine Craft Hydrodynamics and Motion Control; Wiley: Hoboken, NJ, USA, 2011. [Google Scholar]
  3. Bian, X.Q.; Fu, M.Y.; Wang, Y.H. Ship Dynamic Positioning; Science Press: Beijing, China, 2011. [Google Scholar]
  4. Shi, Y.; Shen, C.; Fang, H.; Li, H. Advanced control in marine mechatronic systems: A survey. IEEE/ASME Trans. Mechatron. 2017, 22, 1121–1131. [Google Scholar] [CrossRef]
  5. Perez, T.; Donaire, A. Constrained control design for dynamic positioning of marine vehicles with control allocation. Model. Identif. Control 2009, 30, 32–46. [Google Scholar] [CrossRef]
  6. Fossen, T.I.; Perez, T. Kalman filtering for positioning and heading control of ships and offshore rigs. IEEE Control Syst. Mag. 2009, 29, 32–46. [Google Scholar]
  7. Fung, P.T.-K.; Grimble, M. Dynamic ship positioning using self-tuning Kalman filter. IEEE Trans. Autom. Control 1983, 28, 339–349. [Google Scholar] [CrossRef]
  8. Zhang, Y.H.; Jiang, J.G.; Gao, D.K. Novel disturbance compensating dynamic positioning of dredgers based on adaptive backstepping. J. Southeast Univ. (Eng. Ed.) 2011, 6, 680–688. [Google Scholar]
  9. Muhammad, S.; Dòria-Cerezo, A. Passivity-based control applied to the dynamic positioning of ships. IET Control Theory Appl. 2012, 6, 680–688. [Google Scholar] [CrossRef]
  10. Fossen, T.I.; Grøvlen, A. Nonlinear output feedback control of dynamically positioned ships using vectorial observer backstepping. IEEE Trans. Control Syst. Technol. 1998, 6, 121–128. [Google Scholar] [CrossRef]
  11. Du, J.L.; Hu, X.; Krstić, M.; Sun, Y.Q. Robust dynamic positioning of ships with disturbances under input saturation. Automatica 2016, 73, 207–214. [Google Scholar] [CrossRef]
  12. Tannuri, E.A.; Agostinho, A.C.; Morishita, H.M.; Moratelli, L. Dynamic positioning systems: An experimental analysis of sliding mode control. Control Eng. Pract. 2010, 18, 1121–1132. [Google Scholar] [CrossRef]
  13. Do, K.D. Global robust and adaptive output feedback dynamic positioning of surface ships. J. Mar. Sci. Appl. 2011, 10, 325–332. [Google Scholar] [CrossRef]
  14. Hassani, V.; Sørensen, A.J.; Pascoal, A.M.; Athans, M. Robust Dynamic Positioning of offshore vessels using mixed-μ synthesis modeling, design, and practice. Ocean Eng. 2017, 129, 389–400. [Google Scholar] [CrossRef]
  15. Yang, Y.; Guo, C.; Du, J.L. Robust adaptive NN-based output feedback control for a dynamic positioning ship using DSC approach. Sci. China (Inform. Sci.) 2014, 57, 1–3. [Google Scholar] [CrossRef]
  16. Zhang, G.Q.; Lv, S.N.; Huang, C.F.; Zhang, X.K. Robust adaptive control for dynamic positioning vehicles in presence of adjustable threshold rule and input constraints. Ocean. Eng. 2023, 282, 114950. [Google Scholar] [CrossRef]
  17. Cho, S.; Shim, H.; Kim, Y.A. Dynamical Sliding Mode Control for Robust Dynamic Positioning Systems of FPSO Vessels. J. Mar. Sci. Eng. 2022, 10, 474. [Google Scholar] [CrossRef]
  18. Xi, Y.G.; Li, D.W.; Lin, S. Model Predictive Control—Status and Challenges. Acta Autom. Sin. 2013, 39, 222–236. [Google Scholar] [CrossRef]
  19. Wang, Y.H.; Shi, X.C.; Bian, X.Q. Restriction control of marine dynamic positioning system based on model predictive control. Ship Eng. 2007, 29, 22–25. [Google Scholar]
  20. Wang, Y.H.; Sui, Y.F.; Wu, J. Marine dynamic position system based on nonlinear model predictive control. J. Harbin Eng. Univ. 2013, 34, 110–115. [Google Scholar]
  21. Veksler, A.; Johansen, T.A.; Borrelli, F.; Realfsen, B. Dynamic Positioning with Model Predictive Control. IEEE Trans. Control Syst. Technol. 2016, 24, 1340–1353. [Google Scholar] [CrossRef]
  22. Gierusz, W.; Miller, A. Prediction Control Systems in Marine Applications. Int. J. Mar. Navig. Saf. Sea Transp. 2020, 14, 361–366. [Google Scholar] [CrossRef]
  23. Miller, A. Model predictive ship trajectory tracking system based on line of sight method. Bull. Pol. Acad.-Sci.-Tech. Sci. 2023, 71, e145763. [Google Scholar] [CrossRef]
  24. Wang, L.P. Model Predictive Control System Design and Implementation Using MATLAB; Springer-Verlag London Limited: London, UK, 2009. [Google Scholar]
  25. Kong, X.B.; Liu, X.J. Continuous-time nonlinear model predictive control with inputoutput linearization. Control Theory Appl. 2012, 29, 217–224. [Google Scholar]
  26. Fossen, T.I.; Sagatun, S.I.; Sørensen, A.J. Identification of dynamically positioned ships. Control Eng. Pract. 1996, 4, 369–376. [Google Scholar] [CrossRef]
Figure 1. Horizontal coordinate frames of ship motion.
Figure 1. Horizontal coordinate frames of ship motion.
Jmse 12 00294 g001
Figure 2. Discrete-time Laguerre network.
Figure 2. Discrete-time Laguerre network.
Jmse 12 00294 g002
Figure 3. Response trajectory curves of the DP vessel.
Figure 3. Response trajectory curves of the DP vessel.
Jmse 12 00294 g003
Figure 4. Response curves of positions and headings.
Figure 4. Response curves of positions and headings.
Jmse 12 00294 g004
Figure 5. Response curves of ship velocity and yaw rate.
Figure 5. Response curves of ship velocity and yaw rate.
Jmse 12 00294 g005
Figure 6. Response curves of control forces and moments.
Figure 6. Response curves of control forces and moments.
Jmse 12 00294 g006
Figure 7. Response curves of control forces and moment rates.
Figure 7. Response curves of control forces and moment rates.
Jmse 12 00294 g007
Figure 8. Comparison of computation load.
Figure 8. Comparison of computation load.
Jmse 12 00294 g008
Table 1. The main parameters of the DP supply vessel.
Table 1. The main parameters of the DP supply vessel.
ParameterValueParameterValue
L o a 76.1 m B m 18.8 m
d m 6.25 mDWT4200 t
X−1000∼1000 kNY−300∼300 kN
N−7620∼7620 kN·m Δ X −200∼200 kN/s
Δ Y −60∼60 kN/s Δ N −1524∼1524 kN·m/s
Table 2. Comparison of mean computation time.
Table 2. Comparison of mean computation time.
Time RangeMarkNew Algorithm (ms)Origin
Algorithm (ms)
Difference
1 s∼17 sMin10.9227.9
1 s∼17 sMax304.41847.0
1 s∼17 sMean134.6684.280.3%
17 s∼500 sMin10.113.1
17 s∼500 sMax12.315.9
17 s∼500 sMean10.513.622.8%
TotalMean14.635.759.1%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qian, X.; Shen, H.; Yin, Y.; Guo, D. Nonlinear Model Predictive Control for a Dynamic Positioning Ship Based on the Laguerre Function. J. Mar. Sci. Eng. 2024, 12, 294. https://doi.org/10.3390/jmse12020294

AMA Style

Qian X, Shen H, Yin Y, Guo D. Nonlinear Model Predictive Control for a Dynamic Positioning Ship Based on the Laguerre Function. Journal of Marine Science and Engineering. 2024; 12(2):294. https://doi.org/10.3390/jmse12020294

Chicago/Turabian Style

Qian, Xiaobin, Helong Shen, Yong Yin, and Dongdong Guo. 2024. "Nonlinear Model Predictive Control for a Dynamic Positioning Ship Based on the Laguerre Function" Journal of Marine Science and Engineering 12, no. 2: 294. https://doi.org/10.3390/jmse12020294

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