Next Article in Journal
Mapping Relation between Contour Error Components of Crankshaft Pin Journal and Axis Position Control Error of Oscillating Grinding Machine
Previous Article in Journal
The Assessment of the Condition of Knee Joint Surfaces with Acoustic Emission Analysis
Previous Article in Special Issue
Deep Learning for Efficient and Optimal Motion Planning for AUVs with Disturbances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LQG Control for Dynamic Positioning of Floating Caissons Based on the Kalman Filter

by
Jose Joaquin Sainz
1,
Elías Revestido Herrero
2,*,
Jose Ramon Llata
1,
Esther Gonzalez-Sarabia
1,
Francisco J. Velasco
2,
Alvaro Rodriguez-Luis
3,
Sergio Fernandez-Ruano
3 and
Raul Guanche
3
1
Department of Electronic Technology, Systems Engineering and Automatic Control, E.T.S. de Ingenieros Industriales y de Telecomunicacion, University of Cantabria, Av. de los Castros s/n, 39005 Santander, Spain
2
Department of Electronic Technology, Systems Engineering and Automatic Control, E.T.S. de Náutica, University of Cantabria, C/Gamazo 1, 39004 Santander, Spain
3
Environmental Hidraulics Institute, IH Cantabria, Universidad de Cantabria, c/Isabel Torres n 15, Parque Cientifico y Tecnológico de Cantabria, 39011 Santander, Spain
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(19), 6496; https://doi.org/10.3390/s21196496
Submission received: 8 September 2021 / Revised: 21 September 2021 / Accepted: 24 September 2021 / Published: 29 September 2021
(This article belongs to the Special Issue Marine Wireless Sensor Networks: Applications and New Challenges)

Abstract

:
This paper presents the application of an linear quadratic gaussian (LQG) control strategy for concrete caisson deployment for marine structures. Currently these maneuvers are carried out manually with the risk that this entails. Control systems for these operations with classical regulators have begun to be implemented. They try to reduce risks, but they still need to be optimized due to the complexity of the dynamics involved during the sinking process and the contact with the sea bed. A linear approximation of the dynamic model of the caisson is obtained and an LQG control strategy is implemented based on the Kalman filter (KF). The results of the proposed LQG control strategy are compared to the ones given by a classic controller. It is noted that the proposed system is positioned with greater precision and accuracy, as shown in the different simulations and in the Monte Carlo study. Furthermore, the control efforts are less than with classical regulators. For all the reasons cited above, it is concluded that there is a clear improvement in performance with the control system proposed.

1. Introduction

Nowadays, precast concrete (floating) caissons are used in the construction of maritime infrastructures all over the world. These types of constructions have been highly developed in the last few years [1]. Various operations are carried out to complete the construction of a maritime infrastructure, among which the most important are the following: sea transportation from the construction yard to the target site, positioning at the target site and sinking to its final position. Within the previously-mentioned operations, the most important ones are sinking and positioning. The achievement of these operations entails a level of risk. In addition, a small margin of deviation is allowed, less than 10 cm of deviation from theoretical values. This value is provided by the SAFE Project stakeholder, FCC, and is based in [2].
One of the standard procedures is the fine positioning of the caisson before proceeding with the sinking operations. Due to the complexity of the various operations involved during the whole process, the personnel involved in them must be adequately trained for the correct execution of the whole process. Furthermore, the complex operating environment in which the operations are carried out may provoke inaccuracies in the positioning of the caisson or collision with previous caissons. This may cause material damage to the equipment and even damage to the personnel who carry out the operations. Moreover, apart from accidents this maneuver may lead into inaccurate positioning of the caisson, something which may lead to economic losses. The sinking processes, which take place just after the fine positioning, is also a highly complex maneuver where the sinking velocity needs to be adjusted in order to avoid violent seabed contacts without increasing the positioning potential inaccuracies.
Therefore, and considering all of the above, in this work, we proposed a dynamic positioning (DP) system for anchored caissons, which contributes to the minimization of the operation risks in a way that safeguards the integrity of equipment and people. In addition, and simultaneously, it will make the operations economically less costly due to the reduction of the necessary personnel on board in the course of the operations.
In the state of art relative to vessels or floating structures [3,4], there are several DP approaches, which apply a wide variety of control methods. Among all these methods, the most outstanding are the following: PID control [4,5], acceleration feedback and Kalman filtering [4,6,7,8,9], robust control [10,11,12], Fuzzy control [13], backstepping method [14,15] and others [16,17,18]. Although with some similarities, all the cited references related to floating structures and ships differ from the dynamic behavior of the floating caissons presented in this work. Apart from that, the allocations system used in all the cited contributions is based on thrusters, which is completely different from the one proposed in this work due to being based on mooring lines. Besides, it is found a work where it is used a classical controller for a DP system of caissons [19]. The approach of the present paper differs from this one, since we propose the application of a linear quadratic gaussian (LQG) control strategy in a DP system, which normally demands a compound of staff and equipment when it is performed in the non-automated way.
As far as the LQG control strategy is concerned, some relevant contributions are found in the literature, related to the motion control of underwater vehicles [20,21,22] or their thrusters [23], to cite some of the most representative. In the same way, regarding to stabilization and motion control of ships, there is a variety of references [24,25,26,27]. All of these contributions include thrusters in the allocation system, see [28] for different types of thruster configurations, or [8] for marine vehicles and [29] for a general approach of control allocation algorithms. It is possible to check in the contributions published by [30,31,32,33,34] that DP control for moored vessels is also performed with thrusters, where the mooring system is modelled by means of finite elements [35].
The novelty of the present paper lies in the first-time application of the LQG control method to DP of floating caissons with an allocation system based on mooring lines. This approach differs from those previously reported in the literature, since it is performed without thrusters. We propose the application of an LQG control strategy based on the Kalman filter (KF) to a DP system of anchored floating caissons, see Figure 1. In this proposal, we implement a control allocation system without thrusters but based on teleoperated mooring lines. The relative position between the caisson and the anchor points changes the direction of the forces induced on the caisson by the different mooring lines. This highlights instability to the system, which constitutes a difficulty that the LQG control strategy must overcome.
Subsequently, the results between the classical controller applied to caissons [19] and the controller proposed in this paper will be compared. It will be observed that the proposed controller provides better results in the actuation signals and a reduction in the oscillations under the same perturbations conditions.

2. The Model

2.1. Hydrodynamic Model of the Caisson

The dynamic model of the caisson, object of study of this paper is represented by the following Cummins equation [36,37]:
( M + A ( ) ) + z ¨ ( t ) + 0 t K ¯ ( t τ ) z ˙ ( τ ) d τ + C z ( t )   =   F a 0 ( t ) + F w 0 ( t )
The function of delay and fluid memory effects is:
K ¯ ( t )   =   2 π 0 B ( w ) c o s ( w t ) d w
where z ( t )   =   [ x 0 , y 0 , z 0 , ϕ 0 , θ 0 , ψ 0 ] T are the position and the Euler angles of the caisson. M is the mass of the caisson, A ( ) is the added mass at infinite frequency, K is the function of delay and fluid memory effects,   B ( w ) is the damping coefficient for every frequency, C is the hydrostatic restoration coefficient, F a 0   =   [ X a 0 , Y a 0 , Z a 0 , K a 0 , M a 0 , N a 0 ] are the control forces and moments, F w 0   =   [ X w 0 , Y w 0 , Z w 0 , K w 0 , M w 0 , N w 0 ] are the forces and moments induced by disturbances.

2.2. Linear Approximation

The hydrodynamic model previously presented is nonlinear, see Equation (1). It is well-known that a linear quadratic regulator (LQR) provides a good performance if the dynamic behavior of the system is linear. That is why, in this section, we obtained a linear approximation of the model in Equation (1). Taking into consideration the operational conditions of this paper, which refer to dynamic positioning of caissons at low speed, a linear approximation can properly represent the dynamics of the caisson.
The structure of the linear approximation (Figure 2 and Table 1) is determined based on previous knowledge of the system. To do so, the step response was evaluated by observing fundamental aspects, such as the appropriate sampling period, the natural modes of the system, coupling among the different inputs and outputs of the system, delays, etc. Subsequently, input signals are generated for the model in Equations (1) and (2) based on square wave signals and pseudorandom binary sequence (PRBS) signals that are applied sequentially in each of the inputs of the system, obtaining the outputs response: x 0 ,   y 0 ,   z 0 ,   ϕ 0 ,   θ 0   and   ψ 0 .
Subsequently, by means of least squares [38,39,40,41] a family of parametric models is obtained in the form of transfer functions, which best fit to the data, see Figure 2 and Table 1. Then, in order to properly apply the LQG control strategy, the family models are converted into state space representation of Equations (3) and (4).
The cited linear approximation in state space form is:
x ˙ L ( t )   =   A L x L ( t )   +   B L u L ( t )
y L ( t )   =   C L x L ( t )   +   D L u L ( t )
A L   =   [ 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 0 0 0 0 0 0 0 0 0 1 ]
B L   =   [ b 11 0 0 0 b 22 0 0 0 0 0 0 0 0 0 0 0 0 b 63 ]
C L   =   e y e ( 6 )
D L   =   z e r o s ( 6 , 3 )
where:
  • A L 6 x 6 is the state matrix.
  • x L ( t )   =   [ x L , y L , z L , ϕ L , θ L , ψ L ] T is the state vector.
  • B L 6 x 3 is the input matrix.
  • u L ( t )   =   [ X L , Y L , N L ] T is the input vector.
  • C L 6 x 6 is the output matrix.
  • y L ( t )   =   [ x L , y L , z L , ϕ L , θ L , ψ L ] T is the output vector.
  • D L 6 x 3 is the feed through matrix.
Figure 3, exhibits the fit of the obtained linear approximation for the surge degree of freedom with respect to the data generated by the nonlinear model of Equation (1). This model validation is performed, as usual in the application of parameter estimation theory, by using different data than those used in the estimation. As can be seen, the linear approximation fits the response of the nonlinear model perfectly. There are no discrepancies between the two responses. This indicates that, for the present operating conditions, the linear approximation system can be used.
As previously discussed, due to the low speed at which the operations are performed, the linear approximation correctly represents the system dynamics. The proposed controller works correctly for linear systems. In case the operations were performed at a higher speed, it is possible that this linear approximation would not be valid and, therefore, another type of control would have to be implemented. However, floating caissons perform their operations in a slow motion because their own structure and characteristics prevent the operations to be performed at a high speed. As a result of that, the proposed controller can be considered to be fully valid and applicable to the system for which it has been designed.

2.3. Model of Wave Disturbances

According to [3] the effects of the waves can be divided in two effects:
  • The effects of the first-order waves ( w F ). The effects of first-order waves are small oscillations of zero mean.
  • The effects of the second-order waves (LF). The effect of second-order wave are typically represented by slow drift motions.
The following state space representation can be used to represent the first-order wave effects:
x ˙ w ( t )   =   [ 2 λ w 0 0 0 0 2 λ w 0 0 0 0 2 λ w 0 1 0 0 0 1 0 0 0 1 w 0 2 0 0 0 w 0 2 0 0 0 w 0 2 0 0 0 0 0 0 0 0 0 ] x w ( t ) + [ 2 λ ω 0 σ 2 λ ω 0 σ 2 λ ω 0 σ 0 0 0 ] u w ( t )
y w ( t )   =   [ I d e n t i t y 3 x 3 0 3 x 3 ] x w ( t ) ,
where λ is the damping, the gain K ω   is K ω = 2 λ ω 0 σ , σ   is   the   wave   intensity , ω 0 is the wave dominant frecuency and w is the Gaussian white noise, the state vector is x w ( t )   =   [ x w , y w , z w , φ w , θ w , ψ w ] T and the input vector is u w ( t )   =   [ w x , w y , w ψ , 0 , 0 , 0 ] T . These effects must be added to each of the outputs of the system to be contaminated.
Moreover, the second-order wave drift forces F w = [ X w , Y w , N w ] T are modelled as slow varying bias terms (wiener processes) [3,6]:
X ˙ w ( t )   =   w X
Y ˙ w ( t )   =   w Y
N ˙ w ( t )   =   w N ,
where w X , w Y and w N are sequences of white noise.
Then, the complete model used in this work including first and second-order wave effects is:
x ˙ L ( t )   =   A L x L ( t )   +   B L [ u L ( t ) + F w ( t ) ]   +   w ( t )
y L ( t )   =   C L x L ( t )   +   D L [ u L ( t ) + F w ( t ) ]
y t ( t )   =   y L ( t )   +   y w ( t ) + v ( t ) ,
where v ( t ) is the noise caused by the sensors and this effect is modeled with Gaussian white noise, w ( t ) is the process noise and is also modeled with Gaussian white noise and y t ( t ) is the output vector with the addition of first-order waves and sensor noise.

3. Kalman Filter

The LQG control strategy of this work is based on the KF. We use this filter to estimate the states and the wave effects previously defined in Equations (5) and (6). It must be noted that only the first-order wave effects are estimated with the KF. The second-order wave effects are compensated with a proportional integral (PI) controller.
In order to estimate and filter, the first-order wave Equations (5) and (6), it is necessary to increase the states of the model with which the KF works, as indicated in [7]. The model with augmented states corresponds to:
x ˙ f ( t )   =   [ A a n ] x a n ( t )   +   [ B a n ] u a n ( t )   +   w ( t )
y ˙ f ( t )   =   [ C a n ] x a n ( t )   +   v ( t )
with,
x f ( t )   =   [ x L , y L , z L , ϕ L , θ L , ψ L , x w , y w , z w , φ w , θ w , ψ w ] T
y f ( t )   =   [ x f , y f , ψ f ] T
u a n ( t )   =   [ u L ( t ) u w ( t ) ] ;
A a n   =   [ A L 6 x 6 0 6 x 6   0 6 x 6   A w 6 x 6   ] ;
B a n = [   B L 6 x 3 0 6 x 3     0 6 x 3   B w 6 x 3   ] ;
C a n = [ I d e n t i t y 6 x 6   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 0 0 0 0 0 0 0 0 0 1   ] ;
u a n (t) are the inputs of the caisson and u w ( t ) is the wave input vector, y f ( t ) is the KF output vector, v ( t ) is the noise caused by the sensors and this effect is modeled with Gaussian white noise and w ( t ) is the process noise and also is modeled with Gaussian white noise. This process noise is incorporated to take into account the discrepancy between the model approximation and real dynamic of the caisson.
The algorithm used in the implementation of a KF as indicated in [3,8] is:
  • The matrices design
    Q ( t )   =   Q T   >   0
    R ( t )   =   R T   >   0
    x ^ f ( 0 )   =   x f 0
  • Initial conditions
    P ( 0 )   =   E [ ( x f ( 0 ) x ^ f ( 0 ) ) ( x f ( 0 ) x ^ f ( 0 ) ) T ]   =   P 0
  • The Kalman profit matrix propagations
    K ( t )   =   P ( t ) H T ( t ) R 1 ( t )
    P is the solution of the Riccati algebraic equation
  • Propagation of the estimated state
    x ^ ˙ f ( t )   =   A a n ( t ) x ^ f ( t ) + B a n ( t ) u a n ( t ) + K ( t ) [ y f ( t ) H ( t ) x ^ f ( t ) ]
The tuning parameters are the covariance matrix of the noise R and the covariance matrix of the states Q. The tuning of the KF is done through the entries of the state and measurement noise [8]. If the model is believed to be uncertain, we need to increase the state covariance.
The noise from the sensors was considered to be uncorrelated. Therefore, R was implemented as the diagonal of the sensor noise covariances.
R   =   d i a g o n a l ( σ x , σ y , σ z , σ ϕ , σ θ , σ ψ ) .
The state covariance matrix was also established diagonally for the same considerations
Q   =   d i a g o n a l ( Q x , Q y , Q z , Q ϕ , Q θ , Q ψ ) .

4. LQG Control Strategy

The LQG control is derived from an LQR controller, the difference between them is that the state feedback in an LQG is obtained by the KF [19]. In this work, we use the KF for the states and waves estimation as indicated in the previous section. Then, we apply an LQG control for three degrees of freedom: x f , y f and ψ f , in the DP system of the Figure 1.
Figure 4 shows the structure of the LQG controller. It must be noted that it includes a PI controller in order to compensate the drift effects induced by the second-order wave effect previously described. The controller calculates the vector of forces and moments F c ( t )   =   [ X c , Y c , N c ] T by means of the current position of the caisson and the reference vector R e f ( t )   =   [ x r e f , y r e f , ψ r e f ] T .
The implementation of the LQR regulator was carried out with the model in state space of Equations (3) and (4).
The cost function is stablished by doing [3]:
J ( u )   =   0 [ x L T Q L Q R x L + u L T R L Q R u L ] d t ,
where Q L Q R R nxn is a negative undefined matrix and R L Q R R pxp is positive definite.
The linear feedback of the control state is calculated by minimizing cost function, Equation (23), being:
u L Q R   =   K L Q R   x f
where:
K L Q R   =   R L Q R 1 B L Q R T P L Q R .
To do this, it is necessary to calculate P L Q R , which is the positive and symmetric definite solution of the Riccati algebraic equation:
A L T P L Q R   +   P L Q R A L     P L Q R B L R L Q R 1 B L T P L Q R   +   Q L Q R   =   0 .

5. Control Allocation

The deployment procedure includes eight winches that are the actuators of the system. See Appendix A for more details relative to the winches and mooring lines configuration. These actuators act by means of cables attached to the caisson. Therefore, winches exert different tensions that result in different forces and moments on the caisson.
The control allocation calculates the tensions that each winch must exert on the caisson with the objective to obtain as a result of these actions the forces commanded by the dynamic control in the three control variables already mentioned.
T c b w   =   [ T c b w 1 , T c b w 2 , T c b w 3 , T c b w 4 , T c b w 5 , T c b w 6 , T c b w 7 , T c b w 8 ] .
The algorithm of the control allocation is the following (Algorithm 1) [19]
Algorithm 1 The algorithm of the control allocation
Initialize T m a x = 10 and n l = 8;
Load the control signals [ X c ,   Y c ,   N c ];
for k   ( 1 , , ) do
if X > 0 then
       f 12 = A c 1 ;
else
    if X < 0 then
        f 12 = A c 2
    end if
end if
if Y > 0 then
          f 34 = A c 3 ;
else
    if Y < 0 then
         f 34 = A c 4 ;
    end if
end if
if N > 0 then
           f 56 = A c 5 ;
else
    if N < 0 then
           f 56 = A c 6 ;
    end if
end if
Add the contributions of all of the axis   S = [ f 12 + f 34 + f 56 ] T ;
Adjust the gain  S = S 2   ;
Calculate the tensions of the n l lines T c = S * a b s ( [ y f x , y f y , y f N ] ) ;
for j   ( 1 , ... , n l )   do
    if T(j) > T m a x   then
       T(j) =   T m a x   ;
    end if
end
end
T max maximum value of the tension, n l number of anchor lines.
The matrices of the Algorithm are the following:
A c 1 = [ 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] A c 2 = [ 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] A c 3 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 0 0 0 ]
A c 4 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 ] A c 5 = [ 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1 0 0 0 0 1 0 0 1 1 ] A c 6 = [ 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 1 1 0 1 0 ]
The tensions T c b w are multiplied for gain of the winches and then, forces and moments are provoked in the center of gravity of the caissons, for details relative to this aspect see [19].

6. Simulation Results

The simulations performed in this section are done to verify the behavior of the LQG controller implemented in the DP system shown in Figure 1. All of them are referred to the caisson, whose main particulars are summarized in the Appendix A. The simulations were carried out in the Matlab-Simulink environment with a time step of 0.1 s. These simulations are divided into two groups: in the former, the DP system of Figure 1 with the LQG controller is simulated and it is compared to the results given by a classical controller, see [19] for more details relative to the controller. In the latter, a Monte Carlo study of 200 realizations is developed, where the performance of the LQG controller is compared to the one provided by the classical controller, both of them in the DP system of Figure 1.
In the simulation, the reference vector was set to Ref(t) = [ 5 m ,   4 m ,   0.175 r a d ] T . The KF matrices were tuned by doing:
  • Q = d i a g ( [ 0.001 , 1 , 0.001 , 1 , 0.001 , 1 , 1 × 10 12 , 1 , 1 × 10 12 , 1 , 1 × 10 15 , 1 × 10 12 ] ) .
  • R = d i a g ( [ 1 × 10 3 , 1 × 10 3 , 1 × 10 3 , 1 × 10 3 , 1 × 10 3 , 1 × 10 3 ] ) .
    The LQR matrices were tuned by doing:
  • Q L Q R = [ 4 × 10 6 , 0 , 0 , 0 , 0 , 0 ; 0 , 4 × 10 6 , 0 , 0 , 0 , 0 ; 0 , 0 , 4 × 10 6 , 0 , 0 , 0 ; 0 , 0 , 0 , 4 × 10 6 , 0 , 0 ;
    0 , 0 , 0 , 0 , 4 × 10 6 , 0 ; 0 , 0 , 0 , 0 , 0 , 4 × 10 6 ] .
  • R L Q R = [ 12 × 10 4 , 0 , 0 ; 0 , 12 × 10 4 , 0 ; 0 , 0 , 8 × 10 5 ] .
The wave parameters for all the simulations are σ   =   0.125 , ω 0 = 1.2 and λ = 0.1 , see Figure 5.
The DP system with the LQG controller counteract the second-order drift effects. This counteracting is due to the PI controller included in the first control loop shown in the Figure 4.
In addition, the implemented system follows the references trajectory (a line between the initial and final point, see Figure 6. It can be observed that the response is stable and has no steady state error. Furthermore, the first-order wave effects were filtered by KF implemented in the system. As a result of that, the controller receives the signals free of oscillations and the levels of noise corresponding to the standard instrumentation used on board the caisson. This can be seen in the zoomed part of Figure 7. Moreover, the actuator signals are not saturated. This means that oscillations are quite limited, see Figure 8. In the same way, the tension signal commanded to the winches does not present saturations or excessive oscillations, see Figure 9.
Taking into account the results, the system is able to perform the dynamic positioning of the caisson with high precision, compensating and filtering the effects of the environment in which the operations will take place.
Figure 10 expose the response of the DP system with the classical controller previously cited with a first-order network as a filter. By comparing Figure 6 and Figure 10, it can be seen a more elevated oscillation in the response of the later one. Moreover, a sharp behavior of controller, seen in Figure 11, represents a significant increment of the control efforts, even with saturations with respect to Figure 8. It can be concluded that the LQG control strategy brings out a controller performance closer to their optimal operating range. This will extend the service life of the winches and contribute to a safely maneuvers development.
Finally, a Monte Carlo study was carried out for dynamic position of the caisson at the location x = 0 and y = 0. It is concluded in Figure 12 that the performance of the DP system with the classical controller (Figure 12b) clearly provides a higher dispersion of the position of the caisson compared to the LQG controller (Figure 12a).

7. Conclusions

In this paper, based on a KF, an LQG control strategy has been applied to for dynamic positioning of floating moored concrete caissons. The simulation results evidenced a good compensation of the second-order wave drift motions thanks to the PI controller included the LQG control loop. Additionally, the low oscillation in the performance of the caisson’s motion also indicates a good compensation of the first-order wave induced motions.
Furthermore, the results given by the classical approaches to the ones provided by the LQG control strategy were compared. From this comparison it can be concluded that the LQG control strategy provides a reduction in the oscillations. Moreover, it also contributes to reducing in the control efforts, avoiding potential saturations. Therefore, this may lead to longer life cycles of winches and actuators, as well as to a safer performance of the manor lines. Finally, the Monte Carlo results shows that the LQG control strategy present a significant lower dispersion for station keeping maneuvers compared to the classical controller.
As it is well-known, the application of the LQG control method gives good results if the dynamic behavior of the system is linear. As previously explained, due to the low speed at which the operations are performed, the linear approximation can properly represent the dynamics of the caisson. However, if these operational conditions significantly change, the linear approximation might not correctly represent the behavior of the caisson and the performance of the controller may be degraded.

Author Contributions

Data curation, E.G.-S.; Investigation, J.J.S. and J.R.L.; Methodology, J.J.S., E.R.H. and J.R.L.; Software, A.R.-L. and S.F.-R.; Supervision, F.J.V. and R.G.; Writing–original draft, J.J.S.; Writing–review & editing, E.R.H. All authors have read and agreed to the published version of the manuscript.

Funding

The Spanish FEDER/Ministry of Science, Innovation and Universities—State Research Agency is greatly acknowledged for partially funding our research through the SAFE Project (Desarrollo de un Sistema Autónomo para el Fondeo de Estructuras para Obras Marítimas), GrantAgreement: RTC-2017-6603-4. The Regional Ministry of Universities, Equality, Culture and Sports of the Gov-ernment of Cantabria has supported this work through the ControlFond project (Control De Ve-hículos Subacuáticos No Tripulados Para Supervisión De Estructuras Para Obras Marítimas Fondeadas). The authors would like to thank FCC Construcción CO as a collaborator in the de-velopment of the SAFE Project, specially Victor Florez Casillas and Nuria Cotallo Aguado (Tech-nical Direction/Hydraulic and Maritime Works) and Alvaro de Toro Mingo (Machinery Direction). R. Guanche also acknowledges financial support from the Ramon y Cajal Program (RYC-2017-23260) of the Spanish Ministry of Science, Innovation and Universities.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The Spanish FEDER/Ministry of Science, Innovation and Universities—State Research Agency is greatly acknowledged for partially funding our research through the SAFE Project (Desarrollo de un Sistema Autónomo para el Fondeo de Estructuras para Obras Marítimas), GrantAgreement: RTC-2017-6603-4. The Regional Ministry of Universities, Equality, Culture and Sports of the Government of Cantabria has supported this work through the ControlFond project (Control De Vehículos Subacuáticos No Tripulados Para Supervisión De Estructuras Para Obras Marítimas Fondeadas). The authors would like to thank FCC Construcción CO as a collaborator in the development of the SAFE Project, specially Victor Florez Casillas and Nuria Cotallo Aguado (Technical Direction/Hydraulic and Maritime Works) and Alvaro de Toro Mingo (Machinery Direction). R. Guanche also acknowledges financial support from the Ramon y Cajal Program (RYC-2017-23260) of the Spanish Ministry of Science, Innovation and Universities. Sensors 21 06496 i001

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Caisson Main Features

The caisson and its lines’ distribution are shown in Figure A1. The main data of the caisson, including the moments of inertia, the center of the gravity and center of buoyancy are found in Table A1.
Figure A1. (a) Mooring lines’ distribution of the caisson; left plant view. (b) Mooring lines’ distribution of the caisson; right lateral view.
Figure A1. (a) Mooring lines’ distribution of the caisson; left plant view. (b) Mooring lines’ distribution of the caisson; right lateral view.
Sensors 21 06496 g0a1
Table A1. Weight distribution of the caisson.
Table A1. Weight distribution of the caisson.
N a m e ValueDescription
m (kg) 6,366,000Total mass
LenghtC (m) 33.69Length of the caisson.
WidthC (m) 19.6Width of the caisson.
HeightC (m) 17.3Height of the caisson.
Dl (m) 14Depth of the mooring point, see Figure A1
Ll (m) 60Length of the line, see Figure A1
Lc1 (m) HeightCD + DlHeight from the mooring point to the winche.
Lc2 (m) L l 2 L c 1 2 Distance from the mooring point to the winche.
Cog (m) −2.474Position of the center of gravity in the 𝑍 axis.
Cob (m) −4875Position of the center of buoyancy in the 𝑍 axis.
D9.75Draft.
Ixx (kg m2) 464,290,000Inertia roll
Iyy (kg m2) 865,050,000Inertia pitch
Izz (kg m2) 865,050,000Inertia yaw

References

  1. Cejuela, E.; Negro, V.; Del Campo, J.M.; MartAnt Esteban, M.D.; L-Gutiez, J.S. Recent history, types, and future of modern caisson technology: The way to more sustainable practices. Sustainability 2018, 10, 3839. [Google Scholar] [CrossRef] [Green Version]
  2. Burgos Teruel, Marcelo. Guía de Buenas Prácticas Para la Ejecución de Obras Marítimas; Puertos del Estado: Madrid, Spain, 2008; ISBN 978-84-88975-68-6. [Google Scholar]
  3. Fossen, T.I. Handbook of Marine Craft Hydrodynamics and Motion Control; John Wiley and Sons, Ltd.: Hoboken, NJ, USA, 2011; ISBN 9781119994138. [Google Scholar]
  4. Fossen, T.I. Marine Control Systems: Guidance, Navigation and Control of Ships, Rigs and Underwater Vehicles; Marine Cybernetics: Trondheim, Norway, 2002. [Google Scholar]
  5. Grinyak, V.M.; Pashin, S.S. Control of the Vessel Course using of PID-Regulator under Parametric Uncertainty. IOP Conf. Ser. Earth Environ. Sci. 2020, 459, 22011. [Google Scholar] [CrossRef]
  6. Pérez, T. Ship Motion Control: Course Keeping and Roll Stabilisation Using Rudder and Fins; Springer: Berlin/Heidelberg, Germany, 2005; ISBN 1849969787/1852339594/9781849969789/9781852339593. [Google Scholar]
  7. Grewal, S.M.; Andrews, P.A. Kalman Filtering: Theory and Practice Using MATLAB, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2001; ISBN 0471392545/9780471392545. [Google Scholar]
  8. 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]
  9. Fossen, T.I.; Perez, T. Marine Systems Simulator (MSS). Available online: https://github.com/cybergalactic/MSS (accessed on 23 September 2021).
  10. Zhang, G.; Huang, C.; Zhang, X.; Tian, B. Robust adaptive control for dynamic positioning ships in the presence of input constraints. J. Mar. Sci. Technol. 2019, 24, 1172–1182. [Google Scholar] [CrossRef]
  11. Huang, Y.; Liu, Z.; Huang, W.; Chen, S. Robust H control for nonlinear course system of unmanned surface vessel with polytopic uncertainty based on sum of squares. Trans. Inst. Meas. Control 2021, 43, 390–399. [Google Scholar] [CrossRef]
  12. Lin, X.; Liang, K.; Li, H.; Jiao, Y.; Nie, J. Robust Finite-Time H-Infinity Control with Transients for Dynamic Positioning Ship Subject to Input Delay. Math. Probl. Eng. 2018, 2018, 2838749. [Google Scholar] [CrossRef]
  13. Xu, S.; Wang, X.; Yang, J.; Wang, L. A fuzzy rule-based PID controller for dynamic positioning of vessels in variable environmental disturbances. J. Mar. Sci. Technol. 2020, 25, 914–924. [Google Scholar] [CrossRef]
  14. Chen, H.; Peng, Y.; Zhang, D.; Xie, S.; Yan, H. Dynamic positioning for underactuated surface vessel via L1 adaptive backstepping control. Trans. Inst. Meas. Control 2021, 43, 355–370. [Google Scholar] [CrossRef]
  15. Witkowska, A.; Śmierzchalski, R. Adaptive dynamic control allocation for dynamic positioning of marine vessel based on backstepping method and sequential quadratic programming. Ocean. Eng. 2018, 163, 570–582. [Google Scholar] [CrossRef]
  16. Tutturen, S.A.; Skjetne, R. Hybrid control to improve transient response of integral action in dynamic positioning of marine vessels. IFAC-Pap. OnLine 2015, 28, 166–171. [Google Scholar] [CrossRef] [Green Version]
  17. 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] [Green Version]
  18. Wu, B.-Q.; Zhao, T. Double Closed-loop Control with Input Saturation for DP Vessels based on Disturbance Observer. In Proceedings of the 2020 5th International Conference on Mechanical, Control and Computer Engineering (ICMCCE), Haerbin, China, 25 December 2020; pp. 646–652. [Google Scholar]
  19. Revestido Herrero, E.; Llata, J.R.; Gonzalez-Sarabia, E.; Velasco, F.J.; Sainz, J.J.; Rodriguez-Luis, A.; Guanche, R. Dynamic positioning of floating caissons based on the UKF filter under external perturbances induced by waves. Ocean. Eng. 2021, 235, 109055. [Google Scholar] [CrossRef]
  20. Carlucho, I.; Menna, B.; De Paula, M.; Acosta, G.G. Comparison of a PID controller versus a LQG controller for an autonomous underwater vehicle (Universidad Nacional del Centro de la Prov. de Buenos Aires, Grupo INTELYMEC, Argentina. In Proceedings of the 2016 3rd IEEE/OES South American International Symposium on Oceanic Engineering (SAISOE), Buenos Aires, Argentina, 15 June 2016; p. 6. [Google Scholar]
  21. Wadoo, S.A.; Sapkota, S.; Chagachagere, K. Optimal Control of an Autonomous Underwater Vehicle (New York Institute of Technology, Department of Electrical and Computer Engineering, Old Westbury, NY, United States). In Proceedings of the 2012 IEEE Long Island Systems, Applications and Technology Conference, Farmingdale, NY, USA, 4 May 2012; p. 6. [Google Scholar]
  22. Noh, M.M.; Arshad, M.R.; Mokhtar, R.M. Depth and pitch control of USM underwater glider: Performance comparison PID vs. LQR. Indian J. Mar. Sci. 2011, 40, 200–206. [Google Scholar]
  23. Ahad, M.A.; Ahmad, S.M. Investigation of a 2-DOF Active Magnetic Bearing Actuator for Unmanned Underwater Vehicle Thruster Application. Actuators 2021, 10, 79. [Google Scholar] [CrossRef]
  24. Li, Q.N.T.; Yu, R.; Shan, Q.; Wu, Y.; Ning, J. An LQG Optimal Linear Controller for Fin Stabilizer System of Marine Vessels. In Proceedings of the 2020 International Conference on System Science and Engineering (ICSSE), Kagawa, Japan, 31 August–3 September 2020. [Google Scholar]
  25. Lu, B. LQG Control for a Class of Complex Disturbance Systems Based on Kalman Filter. In Proceedings of the 2018 Chinese Automation Congress (CAC), Xian, China, 30 November 2018; pp. 1997–2000. [Google Scholar]
  26. Sun, H.-F.; Zhao, Z.-B.; Zhao, X. Study on slanting rudder of ship rolling stabilization based on LQG control method. J. Syst. Simul. 2009, 21, 503–506. [Google Scholar]
  27. Peng, X.-Y.; Zhao, X.-R.; Yin, Z.-F. Statistics analysis on LQG control robustness of a ship lateral motion. J. Syst. Simul. 2007, 19, 3–5. [Google Scholar]
  28. Fossen, T.I.; Johansen, T.A. A survey of control allocation methods for ships and underwater vehicles. In Proceedings of the 2006 14th Mediterranean Conference on Control and Automation, Piscataway, NJ, USA, 28 June 2006; p. 6. [Google Scholar]
  29. Johansen, T.A.; Fossen, T.I. Control allocation a survey. Automatica 2013, 49, 1087–1103. [Google Scholar] [CrossRef] [Green Version]
  30. Aamo, O.M.; Fossen, T.I. Controlling line tension in thruster assisted mooring systems. In Proceedings of the 1999 IEEE International Conference on Control Applications (Cat. No.99CH36328), Kohala Coast, HI, USA, 22 August 1999; pp. 1104–1109. [Google Scholar]
  31. Berntsen, P.I.B.; Aamo, O.M.; Leira, B.J.; Sorensen, A.J. Structural reliability-based control of moored interconnected structures. Control Eng. Pract. 2008, 16, 495–504. [Google Scholar] [CrossRef]
  32. Chen, M.; Ge, S.S.; How, B.V.E.; Choo, Y.S. Robust adaptive position mooring control for marine vessels. IEEE Trans. Control. Syst. Technol. 2013, 21, 395–409. [Google Scholar] [CrossRef]
  33. Jayasiri, A.; Nandan, A.; Imtiaz, S.; Spencer, D.; Islam, S.; Ahmed, S. Dynamic Positioning of Vessels Using a UKF-Based Observer and an NMPC-Based Controller. IEEE Trans. Autom. Sci. Eng. 2017, 14, 1778–1785. [Google Scholar] [CrossRef]
  34. Tiwari, K.; Krishnankutty, P. Comparison of PID and LQR Controllers for Dynamic Positioning of An Oceanographic Research Vessel. In Proceedings of the Fifth International Conference in Ocean Engineering (ICOE2019), Madras, India, 28 July–1 August 2019. [Google Scholar]
  35. Aamo, O.M.; Fossen, T.I. Finite element modelling of moored vessels. Math. Comput. Model. Dyn. Syst. 2001, 47–75. [Google Scholar] [CrossRef]
  36. Cummins, W.E. The impulse response function and ship motions. Schiffstechnik 1962, 9, 101–109. [Google Scholar]
  37. Armesto, J.A.; Guanche, R.; Del Jesus, F.; Iturrioz, A.; Losada, I.J. Comparative analysis of the methods to compute the radiation term in Cummins’ equation. J. Ocean. Eng. Mar. Energy 2015, 1, 377–393. [Google Scholar] [CrossRef]
  38. Lennart, L. System Identification: Theory for the User; Prentice-Hall: Hoboken, NJ, USA, 1987; ISBN 0138816409/9780138816407. [Google Scholar]
  39. Norton, J.P. An Introduction to Identification; Academic Press: London, UK, 1988; ISBN 9780125217309/0125217307. [Google Scholar]
  40. Torsten, S.; Petre, S. System Identification; Academic Press: New York, NY, USA, 1989; ISBN 9780131276062/0131276069. [Google Scholar]
  41. Revestido Herrero, E.; Velasco González, F.J.; López García, E.; Moyano Pérez, E. Diseño de Experimentos para la Estimación de Parámetros de Modelos de Maniobra Lineales de Buques. Rev. Iberoam. Autom. Inf. Ind. 2012, 9, 123. [Google Scholar] [CrossRef] [Green Version]
Figure 1. LQG control loop for dynamic positioning of caissons based on the KF.
Figure 1. LQG control loop for dynamic positioning of caissons based on the KF.
Sensors 21 06496 g001
Figure 2. Structure of the linear approximation for a caissons’ draft of 10.75 m.
Figure 2. Structure of the linear approximation for a caissons’ draft of 10.75 m.
Sensors 21 06496 g002
Figure 3. Validation of the linear approximation for the surge degree of freedom. Step input = 10,000 N.
Figure 3. Validation of the linear approximation for the surge degree of freedom. Step input = 10,000 N.
Sensors 21 06496 g003
Figure 4. LQG control strategy for three degrees of freedom x, y and ψ.
Figure 4. LQG control strategy for three degrees of freedom x, y and ψ.
Sensors 21 06496 g004
Figure 5. (a) Systems with LQG controller control results with waves. First-order waves forces and moments applied to the caisson. (b) Second-order waves forces and moments applied to the caisson.
Figure 5. (a) Systems with LQG controller control results with waves. First-order waves forces and moments applied to the caisson. (b) Second-order waves forces and moments applied to the caisson.
Sensors 21 06496 g005
Figure 6. System with LQG controller control results with waves, positions (x,y).
Figure 6. System with LQG controller control results with waves, positions (x,y).
Sensors 21 06496 g006
Figure 7. (a) System with LQG controller control results with waves, position x. (b) System with LQG controller control results with waves, position y and ψ .
Figure 7. (a) System with LQG controller control results with waves, position x. (b) System with LQG controller control results with waves, position y and ψ .
Sensors 21 06496 g007
Figure 8. System with LQG controller control results with waves. Control signals calculated by the double-loop controller.
Figure 8. System with LQG controller control results with waves. Control signals calculated by the double-loop controller.
Sensors 21 06496 g008
Figure 9. (a) System with LQG controller control results with waves, tensions applied in each of the winches 1–4. (b) System with LQG controller control results with waves, tensions applied in each of the winches 5–8.
Figure 9. (a) System with LQG controller control results with waves, tensions applied in each of the winches 1–4. (b) System with LQG controller control results with waves, tensions applied in each of the winches 5–8.
Sensors 21 06496 g009
Figure 10. (a) System with classical controller control results with waves, position x. (b) System with classical controller control results with waves, position y and ψ .
Figure 10. (a) System with classical controller control results with waves, position x. (b) System with classical controller control results with waves, position y and ψ .
Sensors 21 06496 g010
Figure 11. System with classical controller control result with waves. Force in X, Y and N.
Figure 11. System with classical controller control result with waves. Force in X, Y and N.
Sensors 21 06496 g011
Figure 12. (a) System with LQG controller control result with waves. Monte Carlo study of 200 realizations. (b) System with classical controller control result with waves. Monte Carlo study of 200 realizations.
Figure 12. (a) System with LQG controller control result with waves. Monte Carlo study of 200 realizations. (b) System with classical controller control result with waves. Monte Carlo study of 200 realizations.
Sensors 21 06496 g012
Table 1. Family of transfer functions of the linear approximation.
Table 1. Family of transfer functions of the linear approximation.
N a m e T r a n s f e r   F u n c t i o n   A p r o x i m a t i o n
G 11 a n 11 s 2
G 15 a n 15 s 2 + b n 15 s + c n 15 a d 15 s 3 + b d 15 s 2 + c d 15 s
G 22 a n 22 s 2
G 24 a n 24 s 2 + b n 24 s + c n 24 a d 24 s 3 + b d 24 s 2 + c d 24 s
G 26 a n 26 s 2
G 33 a n 33 a d 33 s 2 + b d 33 s + c d 33
G 42 a n 42 s 2 + b n 42 s + c n 42 a d 42 s 3 + b d 42 s 2 + c d 42 s
G 44 a 44 s 2 + b 44 s + c 44
G 55 a n 55 a d 55 s 2 + b d 55 s + c d 55
G 66 a n 66 s 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sainz, J.J.; Revestido Herrero, E.; Llata, J.R.; Gonzalez-Sarabia, E.; Velasco, F.J.; Rodriguez-Luis, A.; Fernandez-Ruano, S.; Guanche, R. LQG Control for Dynamic Positioning of Floating Caissons Based on the Kalman Filter. Sensors 2021, 21, 6496. https://doi.org/10.3390/s21196496

AMA Style

Sainz JJ, Revestido Herrero E, Llata JR, Gonzalez-Sarabia E, Velasco FJ, Rodriguez-Luis A, Fernandez-Ruano S, Guanche R. LQG Control for Dynamic Positioning of Floating Caissons Based on the Kalman Filter. Sensors. 2021; 21(19):6496. https://doi.org/10.3390/s21196496

Chicago/Turabian Style

Sainz, Jose Joaquin, Elías Revestido Herrero, Jose Ramon Llata, Esther Gonzalez-Sarabia, Francisco J. Velasco, Alvaro Rodriguez-Luis, Sergio Fernandez-Ruano, and Raul Guanche. 2021. "LQG Control for Dynamic Positioning of Floating Caissons Based on the Kalman Filter" Sensors 21, no. 19: 6496. https://doi.org/10.3390/s21196496

APA Style

Sainz, J. J., Revestido Herrero, E., Llata, J. R., Gonzalez-Sarabia, E., Velasco, F. J., Rodriguez-Luis, A., Fernandez-Ruano, S., & Guanche, R. (2021). LQG Control for Dynamic Positioning of Floating Caissons Based on the Kalman Filter. Sensors, 21(19), 6496. https://doi.org/10.3390/s21196496

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