Next Article in Journal
A Stereolithographic Model-Based Dense Body Plan Generation Method to Construct a Ship Hydrodynamic Coefficients Database
Previous Article in Journal
Optimization of the Emissions Profile of a Marine Propulsion System Using a Shaft Generator with Optimum Tracking-Based Control Scheme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model Structures and Identification for Fully Embedded Thrusters: 360-Degrees-Steerable Steering-Grid and Four-Channel Thrusters

by
Gerben Peeters
1,2,*,
Muhammad Raheel Afzal
1,
Maarten Vanierschot
1,
René Boonen
1 and
Peter Slaets
1,*
1
KU Leuven, Mechanical Engineering Technology Cluster TC, Campus Group T, Andreas Vesaliusstraat 13, 3000 Leuven, Belgium
2
Research Foundation-Flanders (FWO), Egmontstraat 5, 1000 Brussels, Belgium
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(3), 220; https://doi.org/10.3390/jmse8030220
Submission received: 19 January 2020 / Revised: 10 March 2020 / Accepted: 14 March 2020 / Published: 20 March 2020
(This article belongs to the Section Ocean Engineering)

Abstract

:
The European Watertruck + project introduced a new fleet of self-propelled inland cargo barges to the European waters, in order to induce more sustainable freight transport in the European hinterland. An augmentation of the automation level of this fleet could further advance their competitiveness and potentially pave the way for unmanned inland cargo vessels. The motion control of such a vessel forms a key component in this envisaged automation chain and benefits from the knowledge of the capabilities of the propulsion system, which here envelops a 360-degrees-steerable steering-grid thruster in conjunction with a 360-degrees-steerable four-channel thruster. Therefore, this study details the mechanical design of both thrusters and lists their experimental towing-tank data. Furthermore, two different modelling methods are offered, one theoretically based and one using a multilayer neural network. A model structure comparison, based on a bias-variance trade-off, verifies the adequacy of the theoretical model which got expended with an angle-dependent thrust deduction coefficient. In addition, several multilayer feedforward neural network architectures exemplify their inherent capability to model the complex, nonlinear, flow phenomena inside the thrusters. These identified model structures can additionally improve thrust allocation algorithms and offer better plant models to study more advanced control strategies.

Graphical Abstract

1. Introduction

On average, current inland cargo vessels generate lower external cargo transportation costs compared to other modes of freight transport, where the authors in [1,2,3] define these costs as accidents, air pollution, climate, noise and congestion. Nevertheless, road-based cargo transport presently dominates the European hinterland freight transport [4]. Therefore, the European Commission targets inducing a cargo transport modal shift from road to rail and waterborne transport [5]. Subsequently, the European Watertruck + [6] project aims to facilitate this modal shift by constructing modular push vessels and barges which are either passive or self-propelled, in order to decouple sailing and transshipment time. The authors believe that an increase in the automation levels of these Watertruck + vessels, and inland cargo vessels in general, could further help the envisaged cargo shift. Furthermore, they believe that this automation augmentation could potentially pave the way for future unmanned or autonomous inland cargo vessels.
Therefore, in [7,8], they constructed a scale model of such a European class type I [9] self-propelled inland cargo barge. These self-propelled barges navigate with a non-conventional, over-actuated, fully-embedded actuation system configuration consisting of a 360-degrees-steerable steering-grid thruster in the bow in conjunction with a 360-degrees-steerable four-channel thruster in the stern. Conventionally, inland vessels tend to have one or more propellers, regularly ducted to protect them and to increase their performance, which are often placed in conjunction with multiple rudders to boost their manoeuvrability. Occasionally, the addition of an azimuth or bow thruster further improves their manoeuvrability [10,11]. In comparison, conventional marine vessels have propeller(s)–rudder(s) configurations positioned at the stern [12], and some marine vessels do carry more exotic propulsion systems such as transversal thrusters, azimuth thrusters, podded propellers, contra-rotating propellers, or even water jets for high speed vessels [13]. In addition, be aware that most small unmanned surface vehicles (USVs) also make use of propeller–rudder configurations [14]. However, some unmanned underwater vessels additionally utilise tunnel thrusters to increase their steering behaviour [15].
Within this aforementioned accumulation of more conventional actuation systems, the tunnel and azimuth thrusters seem to show the largest similarities with the thrusters of this study. The general design and performance of the tunnel thrusters are discussed in [16,17]. Whereas the authors in [18] review the different azimuth thruster designs, architectures, and working principles, the authors in [19] additionally review complementary thruster phenomena such as thruster–thruster or thruster–hull interactions, and the author in [20] published experimental data for an azimuth thruster during both static and dynamic operations. Unfortunately, the modelling literature for both systems is quite scarce. The existing research mostly uses the open-water propeller characteristics, as derived by [21], for its starting point on which additional flow phenomena can be added. For example, the authors in [22] built a mathematical manoeuvring model containing azimuth thrusters including thruster–thruster interactions. Similarly, the authors in [15] suggested a tunnel thruster model for an underwater vehicle exhibiting small yaw or pitch angles, and noted the scarcity of available models, probably due to the complexity of the flow phenomena at hand. Albeit certain parallels can be drawn between the actuation system in this study and these more conventional tunnel and azimuth systems, the usability and adequacy of these conventional modelling approaches need to be investigated. Therefore, this study aims to:
(i)
Detail the mechanical design of both the 360-degrees-steerable steering-grid thruster and the 360-degrees-steerable four-channel thruster nested inside the new fleet of Watertruck + barges.
(ii)
Discuss, list, and show experimental towing tank data from both thrusters at different azimuth angles and propeller speeds, at zero advance speed.
(iii)
Investigate the adequacy of the open-water propeller characteristics as the kernel of the theoretical thrust modelling approach.
(iv)
Incorporate azimuth-angle-dependent internal and external hull thrust deduction losses by means of an extension to the theoretical model from (iii)
(v)
Offer an additional artificial neural network modelling approach that might be useful to grip the current, and perhaps future, inherent complex flow phenomena occurring within both thrusters.
The successful completion of these aims offers the future researcher three different thrust characteristic descriptions (i.e., a data table (ii), an extended physically-based model (iii)+(iv), and an artificial neural network (v)) for both actuation systems (i). In addition to the insights that these data descriptions provide, their information can also be leveraged by advanced motion controllers, thrust allocation algorithms, and plant models. This paper continues as follows: first, Section 2 shows the constructed vessel and details its non-conventional over-actuated propulsion system. Afterwards, Section 3 expands the theoretical thrust model with angle-dependent thrust deductions and explains the multilayer neural network model structures. Subsequently, Section 4 lists the experimental towing-tank data and feeds these data to both modelling methods. Finally, Section 5 provides a comparison and discussion of the modelling results, and Section 6 concludes this study.

2. Material

This section starts with a brief description of the overall design philosophy of the new fleet of self-propelled inland barges in Section 2.1. Afterwards, it details the mechanical design of both the steering-grid and the four-channel thruster in Section 2.2 and Section 2.3, respectively. Finally, it concludes with a detailed summary of the main differences between these thrusters and the more conventional actuation systems in Section 2.4.

2.1. Self-Propelled Inland Barges

During its first phase, the Watertruck + project will construct 12 unpropelled barges, 16 self-propelled barges, and 3 push vessels of European Class type I–II. Their final deliverable aims towards building 500 vessels in the European Class range I–IV [6]. In order to gain insight in the automation potential of these vessels, and in order to study the technological feasibility, and current or future limitations, of unmanned inland cargo vessels, the authors constructed a scale model of a self-propelled barge of European Class type I. The real-size counterpart has a length of 38.5 m, a beam of 5.05 m, and a maximal draught of 2.8 m, whereas the KU Leuven scale model has a similar geometry, scaled down with a factor of 8. Figure 1a shows four real-size Watertruck + barges of European Class type I, and one push vessel (slightly visible on the right), and (b) the KU Leuven scale model deployed as an unmanned, autonomous, inland cargo vessel. Figure 2a depicts a three-dimensional design drawing of this scale model where the embedded four-channel thruster (b), and steering-grid thruster (c) can be seen on the left-hand side and right-hand side of the vessel, respectively. This figure shows the vessel with open hatches in (a), which are the grey rectangular extensions on the top of the image, and sails with closed hatches in Figure 1b. As shown on Figure 2, the axes of the body-fixed reference frames which are used for all components throughout this paper point as follows: the x-axis towards the bow, y-axis towards starboard, and the z-axis towards the bottom. Both thrusters have a Kaplan Ka-series propeller without nozzle with a blade area ratio 0.65, and a pitch-diameter ratio of 0.95. The diameter of the steering-grid propeller measures 100 mm, whereas the diameter of the four-channel propeller measures 150 mm.

2.2. Mechanical Design Steering-Grid Thruster

The 360-degrees-steerable steering-grid bow thruster consists of two main parts: the propeller providing the thrust force, and the steering grid which orients the outflow of the water stream and thus provides the steering capabilities. Figure 3a shows a longitudinal cut of this thruster, i.e., along the x z -plane through the origin of Figure 2c. The propeller and its tilted shaft can be seen on the right side of the cut, where an anti-debris grid covers the inlet hole. The left-hand side shows the axis of the steering grid, but the grid itself is not shown to increase the image readability.
Figure 3b shows an abstract top view of the bottom section of this thruster which illustrates the chosen angle convention of the internal control angle of the steering grid, θ i g . The drawn blue arrows show the orientation of the exiting water flow, which is opposite in direction compared to the thrust force, T g . It is assumed that internal control angle, θ i g , is equal to the output angle, θ o g , of T g , i.e., the flow exits the grid in alignment with the grid position in the x y -plane. Therefore, θ i g immediately gives the orientation of T g . On top of that, the steering grid has a static angle, α , of approximately 28° (curved surface), relative to the x-axis in the x z -plane, visible in Figure 2c. To illustrate this angle, the drawn blue arrows in Figure 3a show the exiting water flow in the x z -plane for θ i g = θ o g = 0 ° which can also be seen in in the x y -plane in Figure 3b on the left.

2.3. Mechanical Design Four-Channel Thruster

The 360-degrees-steerable four-channel stern thruster consists of two main parts: its propeller to induce thrust forces, and a 360-degrees-rotatable steering mechanism consisting of half a sphere with an opening of approximately 85° to orient its exiting water flow. Figure 4a shows a longitudinal section of this four-channel thruster, in the x z -plane in the origin of its body fixed reference frame. This section shows the propeller and its shaft together with the steering mechanism (of which only the bottom section is shown) which points towards the stern. This orientation would guide the exiting water flow towards the outlet channel on the left which points to the stern. The other shown outlet channel points to the bow, and the other two transversal outlet channels cannot be seen in this cut. Notice that the two shown outlets have a downwards angle of approximately 15° relative to the x-axis in the x z -plane, measured from the outlet of the internal steering mechanism—whereas the two transversal outlets have no downwards angle but are positioned a bit higher, directly in line of the internal outlet of the steering mechanism. This height–position difference of the outlets can be seen in Figure 2b.
Figure 4b depicts the angular convention of the internal control angle of the four-channel thruster θ i c which orients the theoretical outflow of the water stream, denoted by the blue arrows. Due to the geometric differences of the thruster channels, one cannot assume that θ i c will equal the output angle, θ o c , of the resulting thrust force, T c . The main geometric design effects that cause this discrepancy between θ i c and θ o c seem to be: (i) the fact that the water flow can only exit through one (e.g., Figure 4d at θ i c = 180 °), or a superposition of two (e.g., Figure 4d at θ i c = 45 °) channels simultaneously, (ii) the different channel lengths and shapes (e.g., downwards bends), (iii) internal deflections of the water stream in the x y -plane when θ i c 0°, 90°, 180°, or −90°, and (iv) a potentially remaining angular velocity of the water flow as there is no grid to align the outflow of the water. Therefore, in this study, the resulting four-channel thrust, T c , will be orthogonally decomposed in its longitudinal, T x c , and transversal, T y c , components for different propeller speeds, n c , according to:
T c ( n c , θ i c ) = T x c ( n c , θ i c ) 2 + T y c ( n c , θ i c ) 2 , θ o c = arctan T y c T x c .

2.4. Differences with Conventional Actuation Systems

In summary, the main design differences of both thrusters in this study compared to the more conventional thrusters appear to be the simultaneous occurrence of: (i) having the rotational axes of the propellers perpendicular (or almost perpendicular in case of the steering-grid) to the calm water plane, hence the propellers rotate in a plane parallel to the x y -plane, (ii) being completely nested inside the hull and thus not positioned directly in the ambient flow field, (iii) drawing in water from underneath the vessel hull, (iv) having both an inlet and an outlet positioned in the hull, not unlike jet systems, (v) having an internal deflection or rotation of the water stream from inlet to outlet, and (vi) having a 360-degrees-steerable outflow of the accelerated water provided by a steering mechanism, where the propeller itself does not change its position or orientation.

3. Modelling Structures and Identification Procedures

This study offers two different modelling approaches: one based on the open-water propeller characteristics, handled by Section 3.1, and one using artificial neural networks, explained by Section 3.2. Afterwards, Section 3.3 details the identification procedures for both modelling approaches.

3.1. Propeller Characteristics Thrust Model

This section uses the open-water propeller characteristics model of [21] as a starting point to construct new model structures. This open-water model takes a first order linearised lift force approximation of the propeller blades to calculate their theoretical thrust force. Given the fact that this lift force is a physical characteristic of the blade geometry, it is deemed suitable as a starting point to model both conventional and non-conventional actuation systems that use a propeller (for example, see its use for the tunnel and azimuth thrusters discussed in Section 1). Afterwards, this model will be expanded by adding the more complex flow phenomena that occur in the propulsion systems of this study. The model of [21] found the theoretical thrust force, T t h , to be quadratically dependent on the revolution speed of the propeller shaft, n, and bi-linearly dependent on the shaft speed and axial inflow velocity, V A :
T t h = T n n n 2 + T n V A n V A .
It is common to introduce the non-dimensional thrust coefficient, K T , and non-dimensional advance ratio, J, to express the ship propeller performance, by introducing:
K T = T t h ρ n 2 D p 4 ,
J = V A n D p ,
where ρ is the water density and D p the propeller diameter. Moreover, for the conventional open-water propeller configurations, K T is often adequately approximated by a linear expression in J:
K T = K T 0 + K T J J .
Combining Labels (2), (3), (4), and (5) gives:
T t h = K T 0 ρ D p 4 n 2 + K T J ρ D p 3 n V A ,
which expresses the same information as Label (2) with:
T n n = K T 0 ρ D p 4 , T n v = K T J ρ D p 3 .
Furthermore, V A can be written as a fraction of the total speed of the vessel, U, by introducing the wake factor, w, which accounts for the speed reduction of the flow field at the inlet of the propeller:
V A = ( 1 w ) U .
In addition, the available thrust force for the vessel, T, is a fraction of the generated T t h due to propeller–hull–interaction losses, accounted for by the thrust deduction number, t, which is assumed to be speed-independent:
T = ( 1 t ) T t h .
In conclusion, by accumulating Labels (6), (8), and (9), T becomes:
T = ( 1 t ) K T 0 ρ D 4 n 2 + K T J ρ D 3 n ( 1 w ) U = ( 1 t ) T n n n 2 + T n V A n ( 1 w ) U .
Note that (10) does not model the propeller dynamics, which can have a noticeable effect when manoeuvring at lower speeds [23]. Consequently, for some applications, it might be useful to extend T with these dynamics, as demonstrated in [24]. However, the experimental data of this study do not include transient dynamic measurements, and therefore their effects are omitted henceforth. Furthermore, as described in Section 4.1, the experiments were performed with zero advance speed, hence, K T = K T 0 . Combining this information, (10) tailors further down to:
T = ( 1 t ) T n n n 2
for this study. The model complexity of the thrusters in this study, which is summarised by Section 2.4, lies mainly in the angle-dependent thrust deduction, i.e., making t = f ( θ i ) . Moreover, both internal and external thrust deductions can occur as the water flows over the external hull, into the thruster inlet, and then through the nested thruster, and thus internal hull, itself, resulting in:
t = t i + t e = f ( θ i ) = t i ( θ i ) + t e ( θ i ) = t ( θ i )
where t i and t e respectively indicate the internal and external hull-interaction thrust deductions. Hence, (11) refines to:
T = [ 1 t ( θ i ) ] T n n n 2 = 1 [ t i ( θ i ) + t e ( θ i ) ] T n n n 2 .
Finally, in order to avoid a model structure bias in the identification results, the theoretical thrust force will be investigated and modelled as an unknown polynomial function in n, T t h = f ( n ) = T m ( n ) . This approach does not force T t h to equal T n n n 2 when U = 0 . Accordingly, this method can thus validate whether the originally derived quadratic propeller-speed dependent thrust model indeed provides the best data fit. Ultimately, the thrust forces T to be identified according to the theoretic model in this study turn into:
T = [ 1 t ( θ i ) ] T m ( n ) .

3.2. Feedforward Multilayer Network Thrust Model

This section describes the architecture of feedforward multilayer neural networks that offer an alternative modelling approach to describe the complex flow phenomena discussed in this study. In general, artificial neural networks have one input layer, a number of hidden layers, and one output layer. However, a multilayer feedforward network with a non-polynomial activation function and only one hidden layer can already serve as a universal approximation for any nonlinear function [25]. This statement means that, regardless of the complexity of the flow phenomena that occur in our system, there will always exist network structures that can approximate the data sets. The capacity of these neural networks to model nonlinear characteristics drove the decision to add this alternative modelling approach.
In order to describe the general layout of a multilayer feedforward network, imagine the following example. If one takes an input vector, x R m , and an output vector, y R l , separated by one hidden layer with n h hidden neurons, and connects the hidden layer to y via W R l × n h and to x via V R n h × m with the addition of a bias vector, β R n h , to the latter and an enveloping activation function, σ ( · ) , the following matrix–vector expression will describe the model structure:
y = W σ ( V x + β ) ,
which can also be expressed in an element-wise form, for i = 1, …, l, according to:
y i = r = 1 n h w i r σ ( j = 1 m v r j x j + β r )
Figure 5 illustrates such a network hierarchy with the following model structure: two inputs (m = 2), one hidden layer of three neurons ( n h = 3 ) of which the bias terms are not explicitly shown, two outputs (l = 2), and interconnection weights w i j , v i j .
Be aware that these descriptions, (15) and (16), implicitly use a linear activation function for the output layer, although this is not necessary. Furthermore, the σ ( · ) can have any shape, but typically a nonlinear function, such as sigmoid, tangent, hyperbolic tangent, etc., is used in order to exploit their characteristics when one trains the network with experimental data. Note that, similarly to (15), a second hidden layer could be inserted, resulting in:
y = W σ 2 ( V 2 σ 1 ( V 1 x + β 1 ) + β 2 ) ,
with n h 1 and n h 2 neurons in the hidden layers, connection-weight matrices, W R l × n h 2 , V 2 R n h 2 × n h 1 , V 1 R n h 1 × m , bias vectors, β 2 R n h 2 , β 1 R n h 1 , and activation functions, σ 1 ( · ) and σ 2 ( · ) . This hidden-layer insertion can be repeated to add multiple layers in an equivalent manner. The amount of hidden layers, their number of neurons, and the activation function types can all be chosen by the user. These design aspects offer great flexibility on the one hand, but an infinite amount of possible network architectures on the other hand. Accordingly, a few network designs will be shown in the results section which offers the possibility to compare their performance. These developed multilayer feedforward networks were trained to generate an output vector y (based on their input vector x), which aims to represent the desired output vector, i.e., the measured towing-tank data sets. More details of this network training procedure can be found in Section 3.3.2.

3.3. Model Identification Procedures

The theoretically derived thrust model of Section 3.1 and the neural networks of Section 3.2 provide a θ i -dependent and n-dependent model structure for the modelled thrust forces, T ( n , θ i ) , to which their respective experimental data sets, D ( n , θ i ) , can be fitted. Equation (18) shows the chosen cost function which calculates the least squares error, E, between the models and their data, i.e., between T ( n , θ i ) and D ( n , θ i ) , for P data points. Section 3.3.1 and Section 3.3.2 explain the procedure to minimize (18), and thus maximize the model fit, for each afore-mentioned model structure of which the maximal model complexity was iteratively determined based on its impact on the final residuals:
E = p = 1 P T ( n p , θ i , p ) D ( n p , θ i , p ) 2

3.3.1. Identification Propeller Characteristics Model

The introduction of (14) into (18) produces the generic cost function for the model structures of Section 3.1:
E = p = 1 P [ 1 t ( θ i , p ) ] T m ( n p ) D ( n p , θ i , p ) 2 .
Within this study, the maximum polynomial orders for t ( θ i ) and T m ( n ) were respectively chosen to be quintic and cubic; hence, in this situation, E would equal:
E = p = 1 P [ 1 ( t 5 θ i , p 5 + t 4 θ i , p 4 + t 3 θ i , p 3 + t 2 θ i , p 2 + t 1 θ i , p + t 0 ) ] ( T n n n n p 3 + T n n n p 2 + T n n p ) D ( n p , θ i , p ) 2 .
Furthermore, Section 4 also identifies the model structures of T m ( n ) = T n n n 2 or T n n n n 3 , given that the former aligns with the theoretical model of [21] (see Section 3.1), whereas for t ( θ i ) all coefficients of its polynomial model will be identified. A Python (version 3.6) script performed the minimization of (20) by using the least_squares() function from the open source SciPy packages [26]. This function determined its step size of the optimization problem, i.e., changing the values of the coefficients of t ( θ i ) and T m ( n ) in order to minimize (20), based on a trust region reflective method, based on the algorithm from [27].

3.3.2. Identification Feedforward Multilayer Network Model

The term feedforward refers to the calculation of y based on x, as these data run through the network of Figure 5 from left to right. Afterwards, the cost function (21) can be determined by the insertion of y from, for example (15) or (17), into (18):
E = p = 1 P y p D ( n p , θ i , p ) 2 ,
where y p represents the network output of the p-th data point, which would be generated by the p-th input point, x p . Within this study, the deepest developed network structure embodied two hidden layers; subsequently, its cost function can be defined by:
E = p = 1 P W σ 2 ( V 2 σ 1 ( V 1 x p + β 1 ) + β 2 ) D ( n p , θ i , p ) 2 .
Note that the p inputs for this network, x p , are the same inputs as for the data points D ( n p , θ i , p ) ; hence, (22) refines to:
E = p = 1 P W σ 2 ( V 2 σ 1 ( V 1 ( n p , θ i , p ) + β 1 ) + β 2 ) D ( n p , θ i , p ) 2 .
Analogously to Section 3.3.1, this cost function (23) needs to decrease in order for the model fitness to increase. A common way to alter the distance between the data and the network is to use the back-propagation method as first suggested by [28]. This method derives an analytical expression of the gradient of the cost function for each neuron in the network by means of an iterative calculation from the right to the left in Figure 5, hence the name backpropagation, and uses this gradient to update the weights. Furthermore, one can augment this method by using the Levenberg–Marquardt algorithm to determine the updates of the connection weights [29]. This augmented back-propagation method was used in this study. Finally, in order to construct a network that generalises well to new inputs, one should avoid the over-fitting of data. Several regularisation methods exist to achieve this generalisation of the network. This study shows two regularisation methods: the early-stopping approach [30], and the Bayesian regularisation procedure of [31,32].
The early-stopping regulation divides x in a training, x t , and a validation, x v , set. Afterwards, it trains the network with the training data and tests the updated network with the validation data. When the error on the validation set, E ( x v ) , starts to rise, although the error of the training data, E ( x t ) , is still decreasing, the iterative algorithm stops, whereas, like the name suggests, Bayesian regulation uses the principles of the Bayesian probability theory as its regularisation principle. This regularisation is achieved by adding a term that penalises the size of weights, e.g., for a total of j weights w:
E w = j w j 2 ,
to the cost function (23). In addition, within this approach, the Bayesian probability theory automatically embeds Ocam’s razor principle by punishing too flexible and too complex module structures based on the Bayesian evidence principle [31]. Hence, this approach offers an automatic regularisation method, which boosts the generalisation of the trained networks. A collection of Matlab scripts constructed all the networks within this study. Consequently, the two afore-mentioned regularisation methods were applied by respectively calling the ‘trainlm’ [29], and ‘trainbr’ [33] function calls within the Matlab neural network environment. For the former, a data division of 70% training and 30% validation data was used. Note that this data separation is not necessary for the latter as it has an automatic regularisation approach as mentioned above.

4. Results

This section starts with listing and illustrating the raw experimental towing tank data from both thrusters in Section 4.1. Afterwards, Section 4.2 takes these data sets to identify the models described by (14) from Section 3.3.1, and similarly Section 4.3 trains the networks given by (15) and (17) from Section 3.3.2. All these results will be briefly discussed in each section, whereas an overall discussion and comparison between the results can be found in Section 5.

4.1. Experimental Towing-Tank Data

The experimental data for both thrusters were measured in a towing tank facility in deep water at zero advance speed, i.e., as a bollard pull test, and can be found in Section 4.1.1 for the steering-grid, and in Section 4.1.2 for the four-channel thruster.

4.1.1. Steering-Grid Thruster Data

The steering-grid thruster performed its bollard pull tests as a stand-alone device without enveloping ship hull hence only the internal angle-dependent hull losses occurred, i.e., t = t i = t i ( θ i g ) . In addition, given the assumption that θ i g = θ o g , as mentioned in Section 2.2, the resulting thrust force data from the steering-grid thruster, D g , were measured in alignment with the angle of the steering grid, using a planar beam load cell of the type LCPB series of OMEGA with a maximal error of 1%. Table A1 in Appendix A summarises these different D g measurements at different θ i g and different rotational speeds of the propeller, n g , i.e., D g = D g ( n g , θ i g ) . Subsequently, Figure 6 outlines these data points on a polar plot, assuming symmetry over the longitudinal x-axes. The decrease in thrust force for θ i g [ 150 , 180 ] ° might be explained by the potential occurrence of a recirculation zone of water flow, as in this orientation the exiting water stream is pointed towards the inlet of the steering-grid thruster.

4.1.2. Four-Channel Thruster Data

The bollard pull tests of the four-channel thruster were performed inside half a ship hull, which was split in the transversal direction, i.e., cut in the y z -plane at midship. Consequently, under the assumption that the most crucial external thrust deduction effects, t e , will be captured by using half a ship hull, t ( θ i c ) becomes the sum of t i ( θ i c ) and t e ( θ i c ) . Given the information that for the four-channel thruster, θ o c θ i c , as indicated in Section 2.3, the thrust forces were measured on both the longitudinal x-axis, D x c ( n c , θ i c ) , and transversal y-axis, D y c ( n c , θ i c ) , of the four-channel thruster for different θ i c and propeller speeds, n c . Therefore, similar to T c , the total resulting measured thrust force of the four-channel thruster, D c , can be derived from:
D c ( n c , θ i c ) = D x c ( n c , θ i c ) 2 + D y c ( n c , θ i c ) 2 , θ o c = arctan D y c D x c .
Table A2 in Appendix B lists all the four-channel thruster data points, measured by a UDW3 force/torque sensor from AMTI, with a maximal error of 3% per component (mainly due to the possibility of cross-talk). Thereupon, Figure 7 illustrates these measurements in two different perspectives. Firstly, plot Figure 7a,b show the decomposed measured thrust forces for θ i c [ 0 , 180 ] ° at different n c . Note that the small fluctuations in n c , reported in Table A2, are accumulated around their mean value in order to plot these data sets with one line, resulting in for example: n c = 340 ± 25 rpm. Secondly, plot Figure 7c,d shows exactly the same data as plot Figure 7a,b, but transformed into θ o c and D c , calculated according to (25).
Figure 7c confirms the nonlinear mapping of θ i c to θ o c hypothesis, which appeared to stay consistent over different n c -sets. This n c -independent consistency in the angle mapping seems to hint at a structural, geometrical origin of this nonlinearity. The relatively constant θ o c for θ i c [−60°, −120°] seems to originate from the incapability of the thruster to produce a significant longitudinal force, i.e., D x c , within this angular domain allowing the transversal force D y c to dominate the thrust generation, which results in an almost transversal θ o c . The inverse seems to be true for θ o c when θ i c [−150°, −180°]. Given the currently available data sets, it remains hard to judge the explicit cause of these nonlinearities. However, Section 2.3 ended with a list of possible contributors to this nonlinear mapping.

4.2. Identification of the Propeller-Characteristics-Based Thrust Models

Plugging the data sets from Section 4.1 into the cost function of (20) allows for the identification of the modelled thrust forces according to (14). By varying the order of the polynomials of t ( θ ) and T m ( n ) up to their respective maximums of quintic and cubic, one can compare the resulting model complexities based on the residuals of their cost functions. This comparison offers a bias-variance trade-off during the eventual model selection [34]. The final model selection depends on the requirements of its end user. Section 4.2.1 calculates T g ( n g , θ i g ) by using D g ( n g , θ i g ) , whereas Section 4.2.2 calculates the orthogonal decomposition of T c ( n c , θ i c ) into T x c ( n c , θ i c ) and T y c ( n c , θ i c ) by respectively using D x c ( n c , θ i c ) and D y c ( n c , θ i c ) .

4.2.1. Steering-Grid Thruster Identification Propeller-Characteristics-Based Thrust Models

Due to the absence of the hull in the towing-tank experiments of Section 4.1.1, no external hull thrust deductions occurred (except for the small external surface of the thruster itself which will be neglected here), i.e., t ( θ i g ) = t i ( θ i g ) . This paves the way to identify t e ( θ i ) by the future increase in thrust deduction when the actuation sits inside its hull by using: t ( θ i ) = t i ( θ i ) + t e ( θ i ) . Table 1 summarises the residuals of the cost functions after performing the afore-mentioned model fits. The absolute value of these residuals has no meaning; however, their relative size does as it indicates the relative model fitness of the underlying model structure. The columns denote the different model structures of T m ( n g ) , whereas the rows show the different polynomial orders of t ( θ i g ) . For example, a cubic-order model for t ( θ i g ) in combination with T m ( n g ) = T n n n g 2 can be identified by:
E = p = 1 P [ 1 ( t 3 θ i , p g 3 + t 2 θ i , p g 2 + t 1 θ i , p g + t 0 ) ] ( T n n n p g 2 ) D ( n p g , θ i , p g ) 2 .
Here, as visible in Table 1, the final cost of E would equal 5.38. To further clarify this representation of the results, Table A3 in Appendix C lists all the identified model coefficients for all the model structures of Table 1. This information sheds more light on the cost residuals; for example, Table A3e shows negative values for T n n n which would mean a thrust deduction with increasing n g , hence here the data seem to be over-fitted and thus losing physical meaning. Table 1 also seems to validate that taking T m ( n ) = T n n n 2 generates a simple but robust model which also aligns with the theoretically derived propeller-characteristics model of Section 3.1, as adding a linear propeller-speed dependency, T n , to T m ( n ) results in marginally lower residuals. Nevertheless, this linear term could be added if desired.
Given the adequacy of T m ( n ) = T n n n 2 , Figure 8a plots the identified quadratic, linear, and constant models for t ( θ i g ) with T m ( n ) = T n n n 2 , together with the data points D g ( n g , θ i g ) to illustrate the model fits. In addition, given the rather constant thrust decrease for θ i g [ 150 , 180 ] °, one could also perform the above-mentioned modelling procedure of T g only for the θ i g [ 0 , 150 ] ° domain, and assume T g ( n g , θ i g ) = T g ( n g , 150 °) to be a constant value for θ i g [ 150 , 180 ] °. Figure 8b illustrates this second methodology for T g , using a quadratic fit for t ( θ i g ) . Figure 9 provides another way to visualise these models fits by plotting T g over its n g domain for θ i g = 0°, 90°, and 180° by using a quadratic model fit for t ( θ i g ) and T m ( n ) = T n n n 2 for the aforementioned full θ i g -domain fit in (a), and for the fit incorporating the constant thrust zone in (b). The latter plot, Figure 9b, shows the benefit of adding the constant thrust zone as it provides a better fit for θ i g = 180 ° compared to Figure 9a.

4.2.2. Four-Channel Thruster Identification Propeller-Characteristics-Based Thrust Models

The four-channel thruster generated its data sets embedded inside half a hull, making it impossible to diversify between t i ( θ i c ) and t e ( θ i c ) as they both depend on the same angle. Thus, this section identifies the thrust deduction coefficients as a superposition of both external and internal hull losses. Note that future towing-tank experiments without a hull, as have been done for the steering-grid thruster, could identify t i separately and thus also give t e . Table 2 lists the different residuals of the cost functions of the model fitting algorithm for both T x c in (a), and T y c in (b). This table also gives the impression that T m ( n ) = T n n n 2 forms a simple and robust model structure to represent the theoretical thrust force, which again aligns with the theoretical derivation of Section 3.1. Here, too, adding T n (or T n and T n n n ) only generates marginal improvements compared to the quadratic n c -dependency. Furthermore, this table indicates that the polynomial order of t ( θ i c ) is better modelled by an odd function for T x c and by an even function for T y c .
Figure 10 illustrates the model fits of T x c in (a), and T y c in (b) for different polynomial orders of t ( θ i c ) with T m ( n ) = T n n n 2 . Moreover, this figure also confirms the better fitness of odd t ( θ i c ) functions for T x c and even functions for T y c .
Given the above-mentioned rather constant thrust output of D x c for θ i c [−60°, −120°], and of D y c for θ i c [−150°, −180°], one could model these regions by a constant value. Subsequently, Figure 11a shows a piecewise linear model consisting of a linear fit of T x c for θ i c [0°, −60°], a constant thrust zone where T x c = T x c ( 60 °) for θ i c [−60°, −120°], and finally again a linear fit for θ i c [−120°, −180°], whereas, similarly to Figure 10b, Figure 11b depicts a quartic and quadratic fit of t ( θ i c ) for this angle in the [0°, −150°] domain followed by a constant thrust zone where T y c = T y c ( 150 °). As can be seen, this constant thrust zone offers a better representation of the data compared to the non-physical polynomial over fitting of Figure 10b on the same angular domain.

4.3. Identification of the Feedforward Multilayer Network Thrust Models

Similarly to Section 4.2, this section plugged the data from Section 4.1 into different feedforward multilayer network designs. Within the subsequent figures and tables, the abbreviations ‘br’, and ‘es’ respectively denote a Bayesian regulation and early stopping regulation scheme. These abbreviations will be followed by two numbers which indicate the amount of neurons in each hidden layer (e.g., br5-0, will only have one hidden layer with five neurons, whereas br3-3 would have two hidden layers with three neurons each). Furthermore, all networks used a tan–sigmoid transfer function, t a n s i g ( · ) , for the activation functions of the hidden layers, i.e., σ ( · ) = t a n s i g ( · ) . Finally, note that all the training procedures of the artificial neural networks within this study were initialised with random weights for the network coefficients. This random initialisation can have an impact on the final results. Therefore, all networks were trained five times and the subsequent tables will show the average, minimum, and maximum final residual cost for each network structure over all five of these training cycles. Analogously to Section 4.2, the model structures can be selected based on a bias-variance trade-off which depends on the requirements of the end user.

4.3.1. Steering-Grid Thruster Identification Multilayer Networks

The networks for the steering-grid thruster used n g and θ i g for the input vector, x, as these inputs also generated the experimental data, i.e., D g ( n g , θ i g ) . The output vector of these networks, y, had only one element which constitutes the modelled thrust force, i.e., y = T g . To exemplify the cost function calculation and consequently the identification of its neural network, (27) calculates the residual for a network which models the steering grid with one hidden layer of three neurons by plugging (16) into (21), and thus by minimising:
E = p = 1 P r = 1 3 w r σ ( j = 1 2 v r j x j , p + β r ) D ( n p g , θ i , p g ) 2 ,
with n h = 3 , m = 2 , and l = 1 and thus its index is omitted. Accordingly, in matrix notation W R 1 × 3 , V R 3 × 2 , and β R 3 × 1 . The two input vector elements, i.e., x 1 , p and x 2 , p are n p g and θ i , p g , note that the subscript i from the latter refers to the internal actuation angle and not to an element notation. Depending on the selected training method, (27) can calculate the residual for es-3-0 or br-3-0, shown in Table 3 which lists the residuals of the final cost functions for all the identified network structures. Figure 12a plots the identification results of three different network topologies which only used one hidden layer which had two, three, or five neurons. Moreover, this figure compares the early stopping and Bayesian regulation approaches. Based on this plot, the Bayesian regulation outperformed the early stopping regulation. This performance difference can also be seen in Table 3a. In addition, Figure 12b shows similar identification results for two different network architectures which nested two hidden layers. Here too, one can observe the Bayes regulation method outperforming the early stopping approach, on which Table 3b provides further detail.

4.3.2. Four-Channel Thruster Identification Multilayer Networks

This section constructs different networks with two hidden layers and a Bayesian regulation approach for the four-channel thruster data. These specific network architectures were chosen based on their performance in the results of the steering-grid thruster, handled in Section 4.3.1. Although all the developed networks used n c and θ i c as their input values, three different approaches for the output, y, were selected to illustrate the network-architecture flexibility:
(i)
Fitting the orthogonally decomposed thrust forces separately, i.e., constructing two networks with one output each: one for deriving T x c based on D x c , and analogously one for finding T y c based on D y c .
(ii)
Fitting the orthogonally decomposed thrust forces simultaneously, i.e., constructing one network with two outputs, T x c and T y c .
(iii)
Fitting the resulting thrust force and its orientation simultaneously, i.e., construction one network with two outputs, T c and θ o c .
Table 4 compares the final cost function residuals of different network topologies for these three different cases. Firstly, regarding case (i), (a) summarises the residuals for T x c and (b) for T y c . Secondly, concerning case (ii), (c) shows the residuals of the chosen different network architectures to model this network which has two outputs. Thirdly, for case (iii), (d) shows a variety of model structures that represented this dual-output network. In addition, Figure 13 illustrates the fitness of several identified models for cases (i) to (iii). To plot these identified network structures, one of the five runs of the respective network designs was chosen at random. When one compares case (i) to case (ii), hence comparing the plots (a), (b) with (c), (d), the former seem to have less variance compared to the latter. However, caution should be taken when one compares these plots as they represent different network topologies, have different amounts of neurons and they are a random selection of one of the five network training cycles performed on each design.

5. Discussion

The pure experimental data of Section 4.1 already expose interesting results for both thruster designs. On the one hand, they show the potential occurrence of a recirculation zone of flowing water for certain positions of the steering grid. On the other hand, they highlight the rather biased output angles of the resulting thrust forces of the four-channel thruster. Nevertheless, additional experiments (e.g., setups with different channel geometries and lengths, propeller diameters, steering-mechanism designs, etc.) could provide complementary insights in the complex flow phenomena arising in these thrusters. In addition, experiments with U 0 could identify the wake factor which might depend on both the vessel velocity and the thruster location, i.e., w = f ( U , X t , Y t , Z t ) . Similarly, these experiments could further investigate the adequacy of (5) for these non-conventional thrusters. Note that the ambient flow, when U 0 , does not decrease the thruster performance, but the complex interactions between the water outlet and the surrounding flow can decrease the effective thrust [35]. Furthermore, experiments with dynamic propeller speeds could verify the adequacy of the existing dynamic thrust models, e.g., [24]. Here, one should note that the ambient flow field tends to not have a large impact on the propulsion system dynamics [36].
Given the complexity of the flow phenomena at hand, this study offered two model structure methodologies. Firstly, Section 4.2 validated the theoretical propeller-characteristics thrust model structures of Section 3.1 for both thruster designs by means of a bias-variance trade-off. This trade-off indicates that the quadratic propeller speed dependency indeed forms a simple but generic description of the resulting thrust forces for both embedded thruster designs at zero advance speed, although more complex model structures could be selected if desired by the end user. Moreover, Section 4.2 further illustrated the successful modelling approach to capture the internal control-angle dependent thrust deductions that exist in both propulsion systems. Evidently, the aforementioned additional experimental setups could also help to further refine the propeller-characteristics based thrust model with the above-discussed possible extensions to this model. Secondly, Section 4.3 demonstrated the possibility to use the multilayer feedforward networks from Section 3.2 to grasp the complex water flows. Bearing in mind that these networks can serve as universal approximations of nonlinear functions, they offer the capability to model more complex flow phenomena that might occur during future experiments with dynamic propeller speeds or when U 0 .
Figure 14 compares both model structure methodologies for the thrust forces generated by the four-channel thruster. Figure 14a depicts the longitudinal modelled forces, T x s , whereas (b) shows the transversal modelled forces, T y s . For T x s , a quintic model describes t ( θ i s ) in conjunction with T m ( n ) = T n n n 2 for the propeller characteristics-based model, whereas a Bayesian regulated network of two hidden layers, with three neurons each, forms the neural network structure. A similar approach describes the structures of T y s , only here t ( θ i s ) has a quartic order. In both cases, i.e., (a) and (b), the neural networks seem to outperform the theoretically based models. Particularly in (b), the theoretical model generates a non-physical fit for θ i c [−150°, −180°] whereas the neural network does not. In (a), for θ i c [−60°, −120°] at n c = 1530 rpm, the theoretical model also offers an incorrect higher thrust force compared to the neural network which is closer to the data points. These differences could origin from the tendency of the higher order polynomials of t ( θ i c ) to over fit to the data set, whereas the Ocam’s razor principle from the Bayesian regulation helps to avoid this situation.
Finally, the current and future identified model structures can provide inputs to several applications. For instance, they can illustrate the achieveable forces for a certain thruster over its whole control domain as exemplified by Figure 15. Furthermore, these models can be used to augment certain links in the automation chain such as thrust allocation methods and motion controllers. For example, an optimal thrust allocation algorithm calculates the optimal thrust magnitudes and orientations for all onboard propulsion systems in order to achieve the desired force and moment output which the motion control calculated [37]. Several implementations of these algorithms exist, such as a quadratic programming approach [38,39], or by using a linear pseudo-inverse model [40]. However, these approaches tend to model the thrust forces as angle-independent, whereas the models of this study could be used to provide more accurate, control-angle dependent, representations of these forces. In addition, motion controllers could take the aforementioned and modelled constant thrust zones into account, similarly to incorporating dead-zones [41,42], further improving their performance. Furthermore, model predictive control architectures could also benefit from the suggested physically-based models as they can intrinsically integrate propulsion system limitations and states into their motion control design [14,43]. Moreover, these predictive controllers could also embed the developed neural networks [44].

6. Conclusions

In conclusion, this study offers the mechanical design and experimental bollard-pull data of both non-conventional thrusters on the one hand, and it provides two different modelling methods to represent these data on the other hand. The design and data themselves deliver more insight into the operational capabilities of the thrusters, and help to better understand the occurring flow phenomena inside these thrusters. The suggested modelling methods can be used to augment the automation level of the novel fleet of self-propelled inland barges. The theoretical propeller-characteristics based thrust model provides an adequate data representation which got confirmed by means of a bias-variance trade-off and can be extended with an angle-dependent thrust deduction function for both propulsion systems. In addition, a multilayer feedforward neural network with Bayesian regulation exemplifies a second modelling approach which is inherently capable to model complex nonlinear functions and could offer more protection against the non-physical over-fitting of data.

Author Contributions

Conceptualization, G.P. and P.S.; methodology, G.P., M.V., R.B., and P.S.; software, G.P.; validation, G.P.; formal analysis, G.P.; investigation, G.P.; resources, G.P.; data curation, G.P.; writing—original draft preparation, G.P.; writing—review and editing, G.P., M.R.A., and M.V.; visualization, G.P.; supervision, M.V., R.B., and P.S.; project administration, P.S.; funding acquisition, P.S. All authors have read and agreed to the published version of the manuscript.

Funding

The Flanders Research Foundation (FWO) funds the doctoral scolarship of Gerben Peeters, 1S12517N. The scale model vessel and thrusters were funded by the EFRO-Flanders project “Autonoom Varen in de Westhoek”.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

( · ) c superscript for four-channel thruster variables
( · ) g superscript for steering-grid thruster variables
( · ) x subscript longitudinal decomposition body-fixed frame
( · ) y subscript transversal decomposition body-fixed frame
α static downwards angle steering-grid
β bias vector
σ ( · ) activation function
θ i internal angle steering mechanism
θ o output angle resulting thrust force
npropeller speed
n h amount of neurons in a hidden layer
tthrust deduction number
t i , t e thrust deduction number due to the internal, external hull
t ( θ i ) angle-dependent thrust deduction modelled as polynomial in θ i
t 0 , , t 5 coefficients quintic polynomial of t
wwake factor
xinput vector artificial neural network
youtput vector artificial neural network
Dexperimental data
Eresult cost function
Jadvance ratio
K T dimensionless thrust coefficient
K T 0 , K T J constant and linear term of linear approxmiation K T
Tmodelled thrust force
T m ( n ) theoretical thrust force modelled as a polynomial in n
T n , T n n , T n n n coefficients cubic polynomial of T t h
T t h theoretical thrust force
Utotal velocity vessel
W, Vconnection-weight matrices
X t , Y t , Z t three-dimensional position of thruster

Appendix A

Table A1. Data steering-grid thruster D g ( n g , θ i ) [N].
Table A1. Data steering-grid thruster D g ( n g , θ i ) [N].
θ i [°]500 rpm1000 rpm1500 rpm
01.857.4713.70
302.308.0815.16
601.757.6415.30 1
901.646.6814.52
1201.796.5713.11
1500.813.417.79
1800.793.368.12
1 data point calculated based on Table A3b with a quadratic fit for t ( θ ) .

Appendix B

Table A2. Towing tank measurements of the 360-degrees-steerable four-channel stern thruster.
Table A2. Towing tank measurements of the 360-degrees-steerable four-channel stern thruster.
θ i c [°] n c [rpm] D x c [N] D y c [N] D c [N] θ o c [°]
03546.980.106.980.82
−303482.49−3.093.97−51.14
−453281.71−3.563.95−64.34
−603300.60−5.625.65−83.91
−903250.37−6.866.87−86.91
−1203360.52−4.634.66−83.59
−135337−1.27−2.272.60−119.23
−150334−3.68−0.133.68−177.98
−180360−7.32−0.227.32−178.28
066023.780.6823.791.64
−306537.85−11.4613.89−55.59
−456286.20−12.7214.15−64.01
−606262.01−20.3320.43−84.35
−906172.48−26.1526.27−84.58
−1206302.22−16.6916.84−82.42
−135628−4.82−8.389.67−119.91
−150636−13.15−0.5213.16−177.74
−180665−24.91−0.1324.91−179.70
096051.061.2651.081.41
−3095216.60−24.5029.59−55.88
−4593013.23−27.5630.57−64.36
−609254.61−45.4545.68−84.21
−909254.83−58.3858.58−85.27
−1209385.39−36.3536.75−81.57
−135935−10.45−18.3021.07−119.73
−150934−29.63−1.4429.66−177.22
−180966−53.43−1.0653.44−178.86
0125587.641.9987.661.30
−30124828.34−42.0750.73−56.03
−45121923.00−47.1452.45−63.99
−6012167.72−79.0779.45−84.42
−9012158.45−102.88103.23−85.30
−12012269.23−64.2064.86−81.82
−1351226−18.17−32.1336.91−119.49
−1501226−51.58−2.5751.64−177.15
−1801260−92.01−2.2492.04−178.61
01545134.403.36134.441.43
−30153843.05−64.4477.50−56.25
−45151535.41−72.8380.98−64.07
−60151011.36−121.99122.52−84.68
−90151010.80−160.79161.15−86.16
−120151513.26−99.74100.62−82.43
−1351515−28.21−49.3356.83−119.76
−1501515−79.17−5.2179.34−176.23
−1801550−140.19−4.15140.25−178.30

Appendix C

Table A3. Identified thrust model coefficients of Table 1.
Table A3. Identified thrust model coefficients of Table 1.
(a)
Order t( θ )Cost t 5 t 4 t 3 t 2 t 1 t 0 T n
Quintic25.45 3 . 68 × 10 10 1 . 61 × 10 07 2 . 43 × 10 05 1 . 48 × 10 03 2 . 98 × 10 02 1 . 33 × 10 02 8 . 40 × 10 03
Quartic30.390 5 . 25 × 10 09 1 . 72 × 10 06 1 . 39 × 10 04 2 . 36 × 10 03 1 . 69 × 10 02 8 . 57 × 10 03
Cubic31.4400 1 . 67 × 10 07 7 . 23 × 10 05 4 . 87 × 10 03 7 . 75 × 10 03 8 . 18 × 10 03
Quadratic31.88000 2 . 61 × 10 05 1 . 72 × 10 03 1 . 94 × 10 02 8 . 29 × 10 03
Linear37.010000 2 . 58 × 10 03 4 . 81 × 10 03 9 . 36 × 10 03
Constant62.1200000 2 . 14 × 10 02 7 . 22 × 10 03
(b)
Order t( θ )Cost t 5 t 4 t 3 t 2 t 1 t 0 T nn
Quintic2.76 1 . 86 × 10 10 7 . 75 × 10 08 1 . 11 × 10 05 6 . 85 × 10 04 1 . 63 × 10 02 7 . 76 × 10 03 6 . 37 × 10 06
Quartic3.800 6 . 58 × 10 09 2 . 01 × 10 06 1 . 32 × 10 04 5 . 03 × 10 04 8 . 35 × 10 04 6 . 36 × 10 06
Cubic5.3800 3 . 52 × 10 07 1 . 30 × 10 04 9 . 39 × 10 03 5 . 85 × 10 03 6 . 23 × 10 06
Quadratic7.28000 3 . 06 × 10 05 2 . 55 × 10 03 8 . 08 × 10 03 6 . 57 × 10 06
Linear14.840000 2 . 60 × 10 03 6 . 83 × 10 03 7 . 24 × 10 06
Constant41.2400000 7 . 35 × 10 06 5 . 52 × 10 06
(c)
Order t( θ )Cost t 5 t 4 t 3 t 2 t 1 t 0 T nnn
Quintic27.15 8 . 28 × 10 11 3 . 04 × 10 08 3 . 75 × 10 06 2 . 34 × 10 04 8 . 61 × 10 03 1 . 52 × 10 02 4 . 40 × 10 09
Quartic27.320 7 . 20 × 10 09 2 . 15 × 10 06 1 . 31 × 10 04 1 . 74 × 10 03 9 . 32 × 10 04 4 . 34 × 10 09
Cubic29.1000 4 . 42 × 10 07 1 . 58 × 10 04 1 . 17 × 10 02 1 . 29 × 10 02 4 . 17 × 10 09
Quadratic31.72000 3 . 30 × 10 05 2 . 95 × 10 03 3 . 40 × 10 03 4 . 48 × 10 09
Linear39.850000 2 . 56 × 10 03 4 . 04 × 10 03 5 . 03 × 10 09
Constant64.7900000 2 . 24 × 10 08 3 . 83 × 10 09
(d)
Order t( θ )Cost t 5 t 4 t 3 t 2 t 1 t 0 T nn T n
Quintic1.50 2 . 00 × 10 10 8 . 44 × 10 08 1 . 23 × 10 05 7 . 56 × 10 04 1 . 71 × 10 02 9 . 49 × 10 02 5 . 74 × 10 06 1 . 70 × 10 03
Quartic2.960 6 . 65 × 10 09 2 . 05 × 10 06 1 . 38 × 10 04 1 . 56 × 10 04 3 . 61 × 10 02 5 . 28 × 10 06 1 . 19 × 10 03
Cubic4.4700 3 . 41 × 10 07 1 . 27 × 10 04 9 . 15 × 10 03 3 . 46 × 10 02 5 . 11 × 10 06 1 . 21 × 10 03
Quadratic6.12000 3 . 43 × 10 05 2 . 78 × 10 03 1 . 36 × 10 01 4 . 78 × 10 06 1 . 31 × 10 03
Linear13.350000 2 . 58 × 10 03 3 . 29 × 10 04 5 . 87 × 10 06 1 . 91 × 10 03
Constant39.8300000 7 . 91 × 10 02 4 . 11 × 10 06 1 . 35 × 10 03
Order t( θ )Cost t 5 t 4 t 3 t 2 t 1 t 0 T nnn T nn T n
Quintic0.99 2 . 43 × 10 10 1 . 03 × 10 07 1 . 51 × 10 05 9 . 27 × 10 04 2 . 06 × 10 02 4 . 62 × 10 02 2 . 34 × 10 09 1 . 01 × 10 05 1 . 05 × 10 03
Quartic2.620 9 . 61 × 10 09 2 . 97 × 10 06 2 . 02 × 10 04 2 . 27 × 10 05 5 . 09 × 10 01 1 . 31 × 10 09 6 . 50 × 10 06 5 . 91 × 10 04
Cubic4.1100 3 . 32 × 10 07 1 . 24 × 10 04 8 . 93 × 10 03 3 . 87 × 10 02 1 . 91 × 10 09 9 . 31 × 10 06 8 . 55 × 10 04
Quadratic5.67000 5 . 45 × 10 05 4 . 38 × 10 03 8 . 26 × 10 01 1 . 29 × 10 09 5 . 82 × 10 06 5 . 85 × 10 04
Linear12.770000 2 . 51 × 10 03 2 . 83 × 10 02 3 . 14 × 10 09 1 . 30 × 10 05 1 . 44 e 03
Constant39.2600000 1 . 44 × 10 01 2 . 75 × 10 09 1 . 13 × 10 05 1 . 29 × 10 03

References

  1. Al Enezy, O.; Van Hassel, E.; Sys, C.; Vanelslander, T. Developing a cost calculation model for inland navigation. Res. Transp. Bus. Manag. 2017, 23, 64–74. [Google Scholar] [CrossRef]
  2. Van Essen, H.; Schroten, A.; Otten, M.; Sutter, D.; Schreyer, C.; Zandonella, R.; Maibach, M.; Doll, C. External Costs of Transport in Europe—Update Study for 2008; Technical Report September; CE Delft, Infras, Fraunhofer ISI: Delft, The Netherlands, 2011. [Google Scholar]
  3. Van Essen, H. Sustainable Transport Infrastructure Charging and Internalisation of Transport Externalities; Technical Report December; CE Delft: Delft, The Netherlands, 2018. [Google Scholar]
  4. European Commission. Eurostat; Technical report; European Commission: Brussels, Belgium, 2019. [Google Scholar]
  5. Kallas, S. Transport 2050: Comission Outlines Ambitious Plan to Increase Mobility and Reduce Emmisions; Technical Report March; European Commission: Brussels, Belgium, 2011. [Google Scholar]
  6. The European Watertruck+ Project. Available online: www.watertruckplus.eu (accessed on 14 January 2020).
  7. Peeters, G.; Catoor, T.; Afzal, M.R.; Kotze, M.; Geenen, P.; Van Baelen, S.; Vanierschot, M.; Boonen, R.; Slaets, P. Design and build of a scale model unmanned inland cargo vessel: Actuation and control architecture. J. Phys. Conf. Ser. 2019. [Google Scholar] [CrossRef]
  8. Peeters, G.; Kotze, M.; Afzal, M.R.; Catoor, T.; Van Baelen, S.; Geenen, P.; Vanierschot, M.; Boonen, R.; Slaets, P. An unmanned inland cargo vessel: Design, build, and experiments. Ocean Eng. 2020, 201, 107056. [Google Scholar] [CrossRef]
  9. CEMT. Resolution No.92/2: On New Classification of Inland Waterways. Technical Report 92; European Conference of Ministers of Transport: 1992. Available online: https://www.itf-oecd.org/resolution-no-922-new-classification-inland-waterways (accessed on 9 March 2020).
  10. Liu, J. Impacts of Rudder Configurations on Inland Vessel Manoeuvrability. Ph.D. Thesis, TU Delft, Delft, The Netherlands, 2017. [Google Scholar] [CrossRef]
  11. Liu, J.; Hekkenberg, R.; Rotteveel, E.; Hopman, H. Literature review on evaluation and prediction methods of inland vessel manoeuvrability. Ocean Eng. 2015, 106, 458–471. [Google Scholar] [CrossRef]
  12. Papanikolaou, A. Ship Design; Springer: Berlin/Heidelberg, Germany, 2014; pp. 1–16. [Google Scholar] [CrossRef]
  13. Fossen, T.I. Handbook of Marine Craft Hydrodynamics and Motion Control; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  14. Liu, Z.; Zhang, Y.; Yu, X.; Yuan, C. Unmanned surface vehicles: An overview of developments and challenges. Annu. Rev. Control 2016, 41, 71–93. [Google Scholar] [CrossRef]
  15. Palmer, A.; Hearn, G.E.; Stevenson, P. Modelling tunnel thrusters for autonomous underwater vehicles. IFAC Proc. Vol. IFAC 2008, 41, 91–96. [Google Scholar] [CrossRef] [Green Version]
  16. Beveridge, J.L. Design and Performance of Bow Thrusters; Technical report; Naval Ship Research and Development Center: Bethesda, MD, USA, 1971. [Google Scholar]
  17. Baniela, S.I. The performance of a tunnel bow thruster with slow speed ahead: A revisited issue. J. Navig. 2009, 62, 631–642. [Google Scholar] [CrossRef]
  18. Desai-Patil, V.; Ayare, A.; Mahajan, B.; Bade, S. A review of azimuth thruster. Int. J. Mech. Eng. 2015, 2, 38–41. [Google Scholar] [CrossRef] [Green Version]
  19. Dang, J.; Laheij, H. Dynamic Positioning Conference; Technical Report June; Wärtsilä Propulsion Netherlands BV (WPNL): Rotterdam, The Netherlands, 2004. [Google Scholar]
  20. Triantafyllou, M.S. Steady and Unsteady Dynamics of an Azimuthing Podded Propulsor Related to Vehicle Maneuvering. Ph.D. Theis, MIT, Cambridge, MA, USA, 2004. [Google Scholar]
  21. Blanke, M. Ship Propulsion Losses Related to Automatic Steering and Prime Mover Control. Ph.D. Theis, Technical University of Denmark, Kongens Lyngby, Denmark, 1981. [Google Scholar]
  22. Delefortrie, G.; Tello Ruiz, M.; Vantorre, M. Manoeuvring model of an estuary container vessel with two interacting Z-drives. J. Mar. Sci. Technol. 2018, 23, 739–753. [Google Scholar] [CrossRef]
  23. Yoerger, D.R.; Cooke, J.G.; Slotine, J.J. The influence of thruster dynamics on underwater vehicle behavior and their incorporation into control system design. IEEE J. Ocean. Eng. 1990, 15, 167–178. [Google Scholar] [CrossRef] [Green Version]
  24. Blanke, M.; Lindegaard, K.; Fossen, T. Dynamic model for thrust generation of marine propellers. IFAC Proc. Vol. 2000, 33, 353–358. [Google Scholar] [CrossRef] [Green Version]
  25. Leshno, M.; Lin, V.Y.; Pinkus, A.; Schocken, S. Multilayer feedforward networks with a nonpolynomial activation function can approximate any function. Neural Netw. 1993, 6, 861–867. [Google Scholar] [CrossRef] [Green Version]
  26. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0—Fundamental algorithms for scientific computing in Python. arXiv 2019, arXiv:1907.10121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Branch, M.A.; Coleman, T.F.; Li, Y. A subspace, interior, and conjugate gradient method for large-scale bound-constrained minimization problems. SIAM J. Sci. Comput. 1999, 20, 1–23. [Google Scholar] [CrossRef]
  28. Rumelhart, D.E.; Hinton, G.E.; Williams, R.J. Learning representations by back-propagating errors. Nature 1986, 323, 533–536. [Google Scholar] [CrossRef]
  29. Hagan, M.T.; Menhaj, M.B. Brief papers. Brain Cogn. 1996, 32, 273–344. [Google Scholar] [CrossRef] [Green Version]
  30. Yao, Y.; Rosasco, L.; Caponnetto, A. On early stopping in gradient descent learning. Constr. Approx. 2007, 26, 289–315. [Google Scholar] [CrossRef]
  31. MacKay, D.J.C. A practical Bayesian framework for backpropagation networks. Neural Comput. 1992, 4, 448–472. [Google Scholar] [CrossRef]
  32. Mackay, D.J.C. Probable networks and plausible predictions—A review of practical Bayesian methods for supervised neural networks. Netw. Comput. Neural Syst. 1995, 6, 469–505. [Google Scholar] [CrossRef]
  33. Foresee, F.D.; Hagan, M.T. Gauss-Newton approximation to Bayesian learning. In Proceedings of the International Conference on Neural Networks (ICNN’97), Houston, TX, USA, 12 June 1997. [Google Scholar] [CrossRef]
  34. Pillonetto, G.; Dinuzzo, F.; Chen, T.; De Nicolao, G.; Ljung, L. Kernel methods in system identification, machine learning and function estimation: A survey. Automatica 2014, 50, 657–682. [Google Scholar] [CrossRef] [Green Version]
  35. English, J. The design and performance of lateral thrust units for ships: Hydrodynamic considerations. Quart. Trans. Roy. Inst. Nav. Arch. 1963, 105, 251–279. [Google Scholar]
  36. Saunders, A.; Nahon, M. The effect of forward vehicle velocity on through-body AUV tunnel thruster performance. Ocean. Conf. Rec. 2002, 1, 250–259. [Google Scholar]
  37. Sørdalen, O.J. Optimal thrust allocation for marine vessels. Control Eng. Pract. 1997, 5, 1223–1231. [Google Scholar] [CrossRef]
  38. Liang, C.C.; Cheng, W.H. The optimum control of thruster system for dynamically positioned vessels. Ocean Eng. 2004, 31, 97–110. [Google Scholar] [CrossRef]
  39. Wei, Y.; Fu, M.; Ning, J.; Sun, X. Quadratic programming thrust allocation and management for dynamic positioning ships. TELKOMNIKA Indones. J. Electr. Eng. 2013, 11, 1632–1638. [Google Scholar] [CrossRef]
  40. Ahani, A.; Ketabdari, M.J. Alternative approach for dynamic-positioning thrust allocation using linear pseudo-inverse model. Appl. Ocean Res. 2019, 90, 101854. [Google Scholar] [CrossRef]
  41. Rubio, J.D.J.; Zamudio, Z.; Pacheco, J.; Mújica Vargas, D. Proportional derivative control with inverse dead-zone for pendulum systems. Math. Probl. Eng. 2013, 2013. [Google Scholar] [CrossRef]
  42. Magyar, B.; Hós, C.; Stépán, G. Influence of control valve delay and dead zone on the stability of a simple hydraulic positioning system. Math. Probl. Eng. 2010, 2010. [Google Scholar] [CrossRef] [Green Version]
  43. Li, Z.; Sun, J. Disturbance compensating model predictive control with application to ship heading control. IEEE Trans. Control Syst. Technol. 2012, 20, 257–265. [Google Scholar] [CrossRef]
  44. Mohanty, S. Artificial neural network based system identification and model predictive control of a flotation column. J. Process Control 2009, 19, 991–999. [Google Scholar] [CrossRef]
Figure 1. (a) Four real-size barges, with individual lengths of 38.50 m, and a push vessel from Watertruck + [6], and (b) the KU Leuven scale model with a length of 4.81 m from [7,8] with additional sensors for environment perception.
Figure 1. (a) Four real-size barges, with individual lengths of 38.50 m, and a push vessel from Watertruck + [6], and (b) the KU Leuven scale model with a length of 4.81 m from [7,8] with additional sensors for environment perception.
Jmse 08 00220 g001
Figure 2. (a) slightly tilted top view of the design drawing of the complete vessel, (b) tilted bottom view of the four-channel thruster, adapted from [7,8], and (c) tilted bottom view of the steering-grid thruster, adapted from [7,8].
Figure 2. (a) slightly tilted top view of the design drawing of the complete vessel, (b) tilted bottom view of the four-channel thruster, adapted from [7,8], and (c) tilted bottom view of the steering-grid thruster, adapted from [7,8].
Jmse 08 00220 g002
Figure 3. 360-degrees-steerable steering-grid bow thruster: (a) longitudinal cross section at its symmetry plane, and (b) abstract top view of its bottom section for θ i g = 0 ° and 90°, adapted from [7,8].
Figure 3. 360-degrees-steerable steering-grid bow thruster: (a) longitudinal cross section at its symmetry plane, and (b) abstract top view of its bottom section for θ i g = 0 ° and 90°, adapted from [7,8].
Jmse 08 00220 g003
Figure 4. 360-Degrees-Steerable four-channel stern thruster: (a) longitudinal cross section at its symmetry plane, and (b) abstract top view of its bottom section for θ i c = 45 ° and = 180 °, adapted from [7,8].
Figure 4. 360-Degrees-Steerable four-channel stern thruster: (a) longitudinal cross section at its symmetry plane, and (b) abstract top view of its bottom section for θ i c = 45 ° and = 180 °, adapted from [7,8].
Jmse 08 00220 g004
Figure 5. Multilayer feedforward network with one hidden layer.
Figure 5. Multilayer feedforward network with one hidden layer.
Jmse 08 00220 g005
Figure 6. Polar plot of D g ( n g , θ i g ) [N], adapted from [7,8].
Figure 6. Polar plot of D g ( n g , θ i g ) [N], adapted from [7,8].
Jmse 08 00220 g006
Figure 7. Data measurements four-channel thruster for different n c and θ i c [ 0 , 180 ] °: (a) longitudinal component of D c , (b) transversal component of D c , (c) angle mapping between θ i c and θ o c , and (d) D c , adapted from [7,8].
Figure 7. Data measurements four-channel thruster for different n c and θ i c [ 0 , 180 ] °: (a) longitudinal component of D c , (b) transversal component of D c , (c) angle mapping between θ i c and θ o c , and (d) D c , adapted from [7,8].
Jmse 08 00220 g007
Figure 8. Model fits for T g by: (a) varying the order of t ( θ i g ) with T m ( n ) = T n n n 2 , and (b) a quadratic fit of t ( θ i g ) for θ i g [ 0 , 150 ] ° and a constant thrust zone for θ i g [ 150 , 180 ] ° with T m ( n ) = T n n n 2 .
Figure 8. Model fits for T g by: (a) varying the order of t ( θ i g ) with T m ( n ) = T n n n 2 , and (b) a quadratic fit of t ( θ i g ) for θ i g [ 0 , 150 ] ° and a constant thrust zone for θ i g [ 150 , 180 ] ° with T m ( n ) = T n n n 2 .
Jmse 08 00220 g008
Figure 9. Model fits for T g with θ i g = 0°, 90°, and 180° and T m ( n ) = T n n n 2 , based on the identification methodology of Figure 8a in (a), and of Figure 8b in (b).
Figure 9. Model fits for T g with θ i g = 0°, 90°, and 180° and T m ( n ) = T n n n 2 , based on the identification methodology of Figure 8a in (a), and of Figure 8b in (b).
Jmse 08 00220 g009
Figure 10. Exemplary model fits for (a) T x c and (b) T y c .
Figure 10. Exemplary model fits for (a) T x c and (b) T y c .
Jmse 08 00220 g010
Figure 11. Alternative model fits: (a) piecewise linear model for T x c , and (b) constant thrust zone model for T y c .
Figure 11. Alternative model fits: (a) piecewise linear model for T x c , and (b) constant thrust zone model for T y c .
Jmse 08 00220 g011
Figure 12. Identification results and regulation comparison multilayer feedforward neural networks for the steering-grid thruster: (a) one hidden layer, (b) two hidden layers.
Figure 12. Identification results and regulation comparison multilayer feedforward neural networks for the steering-grid thruster: (a) one hidden layer, (b) two hidden layers.
Jmse 08 00220 g012
Figure 13. Identification results multilayer feedforward neural networks for the four-channel thruster. First, for case (i), the identified T x c in (a) and T y c in (b) with respective residuals 301 . 59 and 514 . 55 . Second, for case (ii), the identified T x c in (c) and T y c in (d) with a network residual of 26.42. Third, for case (iii), the identified T c in (e) and θ o c in (f) with a final residual of 412.87.
Figure 13. Identification results multilayer feedforward neural networks for the four-channel thruster. First, for case (i), the identified T x c in (a) and T y c in (b) with respective residuals 301 . 59 and 514 . 55 . Second, for case (ii), the identified T x c in (c) and T y c in (d) with a network residual of 26.42. Third, for case (iii), the identified T c in (e) and θ o c in (f) with a final residual of 412.87.
Jmse 08 00220 g013
Figure 14. Comparison neural network and propeller-characteristics-based model fits for (a) T x c and (b) T y c .
Figure 14. Comparison neural network and propeller-characteristics-based model fits for (a) T x c and (b) T y c .
Jmse 08 00220 g014
Figure 15. Visualisation of the possible magnitudes and orientations of T g modelled by the quadratic thrust deduction model, i.e., polynomial order of t ( θ i g ) equals two, and a pure quadratic n g -depedency, i.e., T m ( n ) = T n n n 2 , for a continuous area of ( θ i g , n g ) pairs. The outer circumference of this plot shows the same information as the quadratic polynomial at 1500 rpm of Figure 8a.
Figure 15. Visualisation of the possible magnitudes and orientations of T g modelled by the quadratic thrust deduction model, i.e., polynomial order of t ( θ i g ) equals two, and a pure quadratic n g -depedency, i.e., T m ( n ) = T n n n 2 , for a continuous area of ( θ i g , n g ) pairs. The outer circumference of this plot shows the same information as the quadratic polynomial at 1500 rpm of Figure 8a.
Jmse 08 00220 g015
Table 1. Cost function residuals for the different T g models.
Table 1. Cost function residuals for the different T g models.
Model t ( θ i g ) T nnn , T nn , T n T nn , T n T nnn T nn T n
Constant39.2639.8364.7941.2462.12
Linear12.7713.3539.8514.8437.01
Quadratic5.676.1231.727.2831.88
Cubic4.114.4729.105.3831.44
Quartic2.622.9627.323.8030.39
Quintic0.991.5027.152.7625.45
Table 2. Cost function residuals for the different T c models: (a) T x c , and (b) T y c .
Table 2. Cost function residuals for the different T c models: (a) T x c , and (b) T y c .
(a)
Model t ( θ i c ) T n T nn T nnn T nn , T n T nnn , T nn , T n
Constant38,795.7638,795.9538,796.0838,795.4338,795.34
Linear7292.645092.895719.025085.505074.75
Quadratic7066.954860.485500.834853.594843.96
Cubic3511.661267.952168.831267.481267.44
Quartic3035.18738.331630.70737.48737.46
Quintic2644.39313.481207.49312.45312.37
(b)
Model t ( θ i c ) T n T nn T nnn T nn , T n T nnn , T nn , T n
Constant22,722.7521,261.8422,199.0321,255.8321,223.65
Linear22,173.9020,707.7721,684.9720,700.1020,662.85
Quadratic7808.535129.446204.615128.915128.64
Cubic6888.014515.325245.354144.764144.50
Quartic3620.45562.711654.30559.84559.09
Quintic3414.50356.561462.72354.04353.52
Table 3. Residuals cost functions of the multilayer feedforward neural networks for Figure 12: (a) one hidden layer, (b) two hidden layers.
Table 3. Residuals cost functions of the multilayer feedforward neural networks for Figure 12: (a) one hidden layer, (b) two hidden layers.
(a)
ModelAverageMinimumMaximum
es5-037.9112.8671.26
es3-022.8212.8651.23
es2-022.668.3352.80
br5-07.637.627.63
br3-08.388.378.38
br2-010.6010.6010.60
(b)
ModelAverageMinimumMaximum
es3-258.544.16192.96
es3-180.164.16235.70
br3-20.420.420.42
br3-10.370.370.37
br2-18.448.448.44
Table 4. Comparison of the identification residuals: (a) for the networks of (i) that model T x c , (b) for the networks of (i) that model T y c , (c) for the networks of case (ii), and (d) for the networks of case (iii).
Table 4. Comparison of the identification residuals: (a) for the networks of (i) that model T x c , (b) for the networks of (i) that model T y c , (c) for the networks of case (ii), and (d) for the networks of case (iii).
(a)
ModelAverage CostMinimumMaximum
br2-217,390.562338.7077,598.00
br3-2962.83292.151130.50
br3-3253.10220.82301.59
br4-231,144.55211.3177,597.00
br4-3192.64192.64192.64
br4-489.8328.53181.78
(b)
ModelAverage CostMinimumMaximum
br2-213,528.721252.9062,631.00
br3-213,081.49693.3662,634.00
br3-3299.16114.91636.50
br4-2581.44521.76624.52
br4-3388.454.24518.54
br4-4194.910.42486.64
(c)
ModelAverage CostMinimumMaximum
br3-4703.2532.141279.90
br3-5299.0733.31695.43
br5-3587.82558.14657.58
br4-4336.9226.42575.19
br5-3587.82558.14657.58
br5-416.853.7367.21
(d)
ModelAverage CostMinimumMaximum
br3-4719.98585.921256.20
br3-5345.46116.16500.89
br5-327,349.43748.1534,017.00
br4-4408.6788.74513.91
br5-327,349.43748.1534,017.00
br5-433,841.0033,810.0033,863.00

Share and Cite

MDPI and ACS Style

Peeters, G.; Afzal, M.R.; Vanierschot, M.; Boonen, R.; Slaets, P. Model Structures and Identification for Fully Embedded Thrusters: 360-Degrees-Steerable Steering-Grid and Four-Channel Thrusters. J. Mar. Sci. Eng. 2020, 8, 220. https://doi.org/10.3390/jmse8030220

AMA Style

Peeters G, Afzal MR, Vanierschot M, Boonen R, Slaets P. Model Structures and Identification for Fully Embedded Thrusters: 360-Degrees-Steerable Steering-Grid and Four-Channel Thrusters. Journal of Marine Science and Engineering. 2020; 8(3):220. https://doi.org/10.3390/jmse8030220

Chicago/Turabian Style

Peeters, Gerben, Muhammad Raheel Afzal, Maarten Vanierschot, René Boonen, and Peter Slaets. 2020. "Model Structures and Identification for Fully Embedded Thrusters: 360-Degrees-Steerable Steering-Grid and Four-Channel Thrusters" Journal of Marine Science and Engineering 8, no. 3: 220. https://doi.org/10.3390/jmse8030220

APA Style

Peeters, G., Afzal, M. R., Vanierschot, M., Boonen, R., & Slaets, P. (2020). Model Structures and Identification for Fully Embedded Thrusters: 360-Degrees-Steerable Steering-Grid and Four-Channel Thrusters. Journal of Marine Science and Engineering, 8(3), 220. https://doi.org/10.3390/jmse8030220

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