Seakeeping analysis of dead ship condition in fishing ships based on Artificial Neural Networks

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Seakeeping analysis of dead ship condition in fishing ships based on Artificial Neural Networks'''</big></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
P. Romero-Tello & JE. Gutierrez-Romero</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
''Universidad Politécnica de Cartagena (UPCT), Cartagena, Murcia, Spain''</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
B. Serván-Camas</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
''Centre Internacional de Mètodes Numèrics en Enginyeria (CIMNE), Barcelona, Spain''</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
JI. Piazzese</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
''Compass Ingeniería y Sistemas (CompassIS), Barcelona, Spain''</div>
-->
==Abstract==

In the operation of ships, assessing seakeeping performance is crucial. Historically, this has been done through experimentation in towing tank basins or numerical computations. However, with the rise of Artificial Intelligence (AI) and increased computational resources, there are many opportunities to use AI in predicting seakeeping performance. This research will utilize a pre-trained Artificial Neural Network (ANN) to evaluate the behaviour of fishing vessels in various operational scenarios. One of the key advantages of using these algorithms is the ability to predict a large number of scenarios quickly, compared to traditional methods. By analysing millions of variations in the principal dimensions of a fishing ship and different sea states, the study aims to identify the optimal seakeeping performance in challenging conditions, ultimately improving ship safety by examining principal form coefficients and dimensions. The research will also determine significant conclusions.

'''Keywords''': Fishing ships optimization, artificial intelligence, artificial neural network, seakeeping optimization

==1. Introduction==

Dead ship condition (DSC) is known as the condition of a ship when subject to the waves without neither propulsive nor steering capabilities. This condition is one of the most dangerous ones in terms of safety at sea, and the danger increase as the seastate becomes more energetic. Based on the Spanish Comission for Maritime Accidents (''Comisión de Accidentes e Incidentes Marítimos del Gobierno Español'', CIAIM), the number of maritime accidents between 2010 and 2019 has increased [1], and a large percentage of this accidents happened in DSC. And most of the accidents under DSC had severe consequences (<span id='cite-_tabla1'></span>[[#_tabla1|Table 1]]).

<div class="center" style="font-size: 75%;">'''Table 1'''.  Number of maritime accidents between 2010 and 2019 [1]</div>

<div id='_tabla1'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! '''Year''' !! '''N*''' !! '''DS*'''!! '''VS*''' !! '''S*''' !! '''M*''' !! '''I*'''
|-style="text-align:center"
|  style="text-align: center;"|2010
|  style="text-align: center;"|91
|  style="text-align: center;"|2
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|-style="text-align:center"
|  style="text-align: center;"|2011
|  style="text-align: center;"|101
|  style="text-align: center;"|3
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|  style="text-align: center;"|-
|-style="text-align:center"
|  style="text-align: center;"|2012
|  style="text-align: center;"|60
|  style="text-align: center;"|5
|  style="text-align: center;"|2
|  style="text-align: center;"|1
|  style="text-align: center;"|0
|  style="text-align: center;"|2
|-
|  style="text-align: center;"|2013
|  style="text-align: center;"|123
|  style="text-align: center;"|7
|  style="text-align: center;"|1
|  style="text-align: center;"|1
|  style="text-align: center;"|3
|  style="text-align: center;"|2
|-
|  style="text-align: center;"|2014
|  style="text-align: center;"|97
|  style="text-align: center;"|16
|  style="text-align: center;"|4
|  style="text-align: center;"|3
|  style="text-align: center;"|4
|  style="text-align: center;"|5
|-
|  style="text-align: center;"|2015
|  style="text-align: center;"|204
|  style="text-align: center;"|87
|  style="text-align: center;"|1
|  style="text-align: center;"|50
|  style="text-align: center;"|4
|  style="text-align: center;"|32
|-
|  style="text-align: center;"|2016
|  style="text-align: center;"|316
|  style="text-align: center;"|209
|  style="text-align: center;"|0
|  style="text-align: center;"|202
|  style="text-align: center;"|6
|  style="text-align: center;"|1
|-
|  style="text-align: center;"|2017
|  style="text-align: center;"|298
|  style="text-align: center;"|181
|  style="text-align: center;"|0
|  style="text-align: center;"|177
|  style="text-align: center;"|2
|  style="text-align: center;"|2
|-
|  style="text-align: center;"|2018
|  style="text-align: center;"|274
|  style="text-align: center;"|178
|  style="text-align: center;"|0
|  style="text-align: center;"|176
|  style="text-align: center;"|2
|  style="text-align: center;"|0
|-
|  style="text-align: center;"|2019
|  style="text-align: center;"|283
|  style="text-align: center;"|208
|  style="text-align: center;"|1
|  style="text-align: center;"|196
|  style="text-align: center;"|6
|  style="text-align: center;"|5
|-
|  colspan='7'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 80%;">(*) N = Number of accidents; DS= ''Dead Ship''; VS= Very severe; S= Severe; M= Moderate; I= Incident.</span>
|}


Early 2020, after a long and demanding process, the second generation of intact ship stability criteria (SGISc) of the International Maritime Organization (OMI) were ultimate. [2]. And within the framework of the SGISc measure that provide orientation and limitations during navigation have been implemented to contribute to the design of safer ships [3]. A number of works regarding the implementation of the SGISc can be found in the literature. Among them we find Mata-Álvarez-Santullanol and Pérez-Rojas [4], where a set of 10 small fishing ships are analysed in terms of stability under DSC. In this work it is concluded that the SGISc provides orientation towards selecting among different options during the ship design stage. In Choi et al. [5] the probability of capsizing in DSC under the action of waves and wind is analysed. This analysis is also combined with a first order method to evaluate the fidelity when defining the potential scenarios leading to capsizing. Kubo et al. [6] developed a numerical model with coupled sway, heave, roll and pitch motions and validated it against model testing under artificial conditions of irregular waves and unsteady wind.

Gu et al. [7] remark that within the SGISc the IMO has included some method for directly evaluating the stability under DSC. The authors presented model testing and computational results to validate the latest concluding that more accurate numerical models need to be developed to evaluate stability in DSC, including the effects of drift and sway on roll motion. Shin and Chung [8] recently analysed thousands of tanker ships to transport oil derived chemical products under the SGISc concluding that the SGISc can provide clues to predict in advance the stability in DSC in the event of a ship engine failure. And Hu et al. [9] analysed the seakeeping response of a damaged ship (DTMB 5415) under beam waves in DSC using the method of volume of fluids. It was found that the roll motion amplitude of the damaged ship is smaller than that of the intact ship. Given the coupling motion of the wave excitation and the heeling moment produced by flooding water, various differences are observed in the value of the roll amplitude, and the flooding water effect turns into the damping of ship rolling.

The irruption of Artificial Intelligence in the last years is opening up new possibilities to cope with the seakeeping problem. Several works can be already found in the literature using machine learning techniques to predict the seakeeping performance of ships. Liu et al. [10] developed a Reservoir Computing capable of predicting the movements of the six degrees of freedom for the KVLCC2 ship under irregular seawaves. Taghva et al. [11] predicted heave roll and pitch response amplitude operators (RAOs) and added resistance in waves obtained by the strip theory for the S175 hull using artificial neural networks (ANNs). In [12,13] Cepowski continued that line of research applied to ferries and passenger ships. And in [14] Cepowski trained ANNs to predict the added resistance in waves using basic ship form factors. The ship database contained 14 different ships.

The operational requirements of ships have been traditionally assessed either by seakeeping computational models or by experimental tests. Both methods are time consuming and require the exact geometry and mass particulars of the ship, so they are not applicable in early design stages where only few form coefficients are known.

One of the most recent works in this line of research has been published by the authors of the present work [15]. In [15], Romero-Tello et al. trained ANNs to predict the seakeeping hydrodynamic loads. A large number of different ships were used for an extensive training, resulting on a dataset of more than 20.000 ships under different wave frequencies and directions. As a result, obtained an ANN capable of predicting the seakeeping hydrodynamic loads for any type of conventional displacement monohull ship.

==2. Objectives==

It would be desirable to be able to assess the safety of ships in the early design stages to analyse, for instance, DSCs, and try to minimize the effects on the crew, passage, and avoid capsizing. And the main objective of this work is to provide an answer based on AI techniques.

It this work, the operational parameters of an existing fishing ship (''Ciudad de Cartagena'', see [[#_figura1|Figure 1]] and [[#_figura2|Figure 2]]) will be optimized using the DSC under different seastates. The ANNs developed in [15] will be used to predict the seakeeping loads based on the form coefficients of the new hulls obtained by parametric transformation of the original ship. Since the ANNs can predict thousand of ships per second, the optimization process can be carried out through an extensive analysis of the parametric space.

The objective is to find the ship with the best response to the DSC subject to predefined wave scenarios. This search will be carried out across millions of parametric transformations of the original ship. And the selection of the best ship will be based on the following key performance indicators (KPI): motion sickness incidence (MSI); vertical and horizontal root mean square (RMS) accelerations; and maximum roll and pitch angles (MRA and MPA).

This work is organized as follows: the case study is presented in section 3. Section 4 introduces the optimization methodology. Section 5 shows the results obtained from the optimization, as well as the ship selected as the best. Section 6 provides the conclusions of this work.

==3. Case study==

[[#_figura1|Figure 1]] shows the original ship (''Ciudad de Cartagena'') to be used as case study. This is a fishing ship currently used as a school ship by the ''Universidad Politécnica de Cartagena (UPCT).'' [[#_tabla2|Table 2]] provides the main particulars under fully loaded conditions and [[#_figura2|Figure 2]]  provides the body plan.

<div id='_figura1'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[File:Review_551151435449_5135_Figura_1.png]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 1'''. Ciudad de Cartagena general view
|}


<div id='_figura2'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image2.png|600px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 2'''. Ciudad de Cartagena body plan
|}


<div class="center" style="font-size: 75%;">'''Table 2'''. Ciudad de Cartagena main particulars under loaded condition</div>

<div id='_tabla2'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
|  style="text-align: left;"|'''Overall length (m)'''
|  style="text-align: right;"|24.000
|-
|  style="text-align: left;"|'''Waterline length (m)'''
|  style="text-align: right;"|21.650
|-
|  style="text-align: left;"|'''Breadth (m)'''
|  style="text-align: right;"|6.661
|-
|  style="text-align: left;"|'''Draft (m)'''
|  style="text-align: right;"|2.980
|-
|  style="text-align: left;"|'''Block coefficient C<sub>B</sub>'''
|  style="text-align: right;vertical-align: bottom;"|0.472
|-
|  style="text-align: left;"|'''Waterlines area coefficient C<sub>F</sub>'''
|  style="text-align: right;vertical-align: bottom;"|0.852
|-
|  style="text-align: left;"|'''Midship area coefficient C<sub>M</sub>'''
|  style="text-align: right;vertical-align: bottom;"|0.750
|-
|  style="text-align: left;"|'''Prismatic coefficient C<sub>P</sub>'''
|  style="text-align: right;"|0.630
|-
|  style="text-align: left;"|'''Longitudinal position buoyancy centre XB (m)'''
|  style="text-align: right;"|10.975
|-
|  style="text-align: left;"|'''Vertical position buoyancy centre KB (m)'''
|  style="text-align: right;"|1.969
|}


During its operational life as a fishing ship, the fishing campaigns took place in the areas of the Mediterranean Sea and South Atlantic. And for this analysis, a period of one operational year will be used to determine the wave conditions in the areas. [[#_tabla3|Table 3]] shows the significant wave height (<math display="inline">{H}_{s}</math>), peak period (<math display="inline">{T}_{p}</math>) and wave spectrum (<math display="inline">{S}_{\xi \xi }( \omega )</math>) selected based on [16].

<div class="center" style="font-size: 75%;">'''Table 3'''. Significant wave height, peak period, and wave spectrum (<math display="inline">{S}_{\xi \xi }( \omega )</math>)</div>

<div id='_tabla3'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! style="text-align: left;"| '''Zona''' !! <math display="inline">{\bf H}_S</math> (m)''' !!<math display="inline">{\bf T}_p</math> (s)'''!!'''Wave spectrum'''
|-
|  style="text-align: left;"|Mediterranean Sea
|  style="text-align: center;"|2.347
|  style="text-align: center;"|6.0
|  style="text-align: center;"|JONSWAP
|-
|  style="text-align: left;"|South Atlantic
|  style="text-align: center;"|3.570
|  style="text-align: center;"|8.8
|  style="text-align: center;"|JONSWAP
|}

==4. Optimization methodology==

===4.1. Pre-trained ANNs===

In this work, the ANN developed by Romero-Tello et al. [15] will be used to predict the seakeeping loads of convectional monohull ships. By seakeeping loads we refer to the wave radiation loads (added mass and damping), as well as the excitation loads (Froude-Krylov and diffraction wave loads), in the frequency domain. The ANN developed in Romero-Tello et al. [15] can predict these loads based only on the hull form coefficients, with no need of the exact hull geometry ([[#_figura3|Figure 3]]). And the average error when compared to traditional 3D frequency-domain seakeeping codes based on potential flow and the boundary element method (BEM) is below 5%. The dataset of ships used for training this ANN was composed of 2.0·10<sup>4</sup> different hulls. This dataset was built out of 50 different conventional monohull ships (bulkcarriers, cruise ships, crude carriers, fishing ships, etc.), and then augmented by carrying out 20x20 parametric transformations of L/B and B/T for each. The main advantage of using the ANN developed in Romero-Tello et al. [15] is the computational speed, which allows to compute the KPIs of two cases per second while keeping a good precision.

<div id='_figura3'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image3.png|282px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 3'''. ANN architecture for added mass prediction
|}


To assess the viability of the ANN, a comparison of the ANN and a BEM developed in Delhommeau [16] results for the original ship of the case study is carried out. It has to be said that the Ciudad de Cartagena was not used in the training of the ANN, so that the results are inferred from the results from other ships. Heave and pitch RAOs under head, following and beam seas are compared.

For the BEM computations, a mesh of 6400 elements has been generated for the Ciudad de Cartagena ship ([[#_figura4|Figure 4]]). [[#_figura5|Figures 5]], [[#_figura6|6]] and [[#_figura7|7]] show a comparison between RAOs obtained by ANN and BEM computations for head waves, beam waves and following waves respectively. Heave (RAO<sub>33</sub>) and pitch  (RAO<sub>55</sub>) response amplitude operator curves define the vessel response, for one particular degree of freedom, to one particular wave direction, amplitude and period.

<div id='_figura4'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image4-c.png|480px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 4'''. BEM mesh for Ciudad de Cartagena hull
|}


<div id='_figura5'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|align="center" | 
{|style="margin: 0em auto 0.1em auto;width:auto;" 
|+
|-
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image5.png|400px]] 
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image6.png|400px]] 
|}
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 5'''. RAOs comparison between ANN and BEM for head waves
|}


<div id='_figura6'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|align="center" | 
{|style="margin: 0em auto 0.1em auto;width:auto;" 
|+
|-
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image7.png|400px]] 
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image8.png|400px]] 
|}
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 6'''. RAOs comparison between ANN and BEM for beam waves
|}


<div id='_figura7'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|align="center" | 
{|style="margin: 0em auto 0.1em auto;width:auto;" 
|+
|-
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image9.png|400px]] 
|style="text-align: center;padding:10px;"| [[Image:Review_551151435449-image10.png|400px]]
|}
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 7'''. RAOs comparison between ANN and BEM for following waves
|}


===4.2. Key performance indicators (KPIs)===

Seakeeping stability criteria limit the maximum roll and pitch angles, as well as the maximum horizontal and vertical root mean square (RMS) accelerations. There exist a number of works investigating those criteria such as Rumawas et al. [17], Pattison and Sheridan [18], Ghaemi and Olszewski [19], and Stevens and Parsons [20]. They analyzed the seakeeping criteria provided by the North Atlantic Treat Organization (NATO) STANAG 4154, the NORDFORSK (Nordic Council of Ministers that provides funding for and facilitates Nordic cooperation on research and research infrastructure), and the USCG (United States Coast Guards, Cutters certification plan), which are provided in [[#_tabla4|Tables 4]], [[#_tabla5|5]], and [[#_tabla6|6]].

<div class="center" style="font-size: 75%;">'''Table 4'''. NATO STANAG 4154 [21]</div>

<div id='_tabla4'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
!style="text-align: left;"| '''Recommended criterium''' !! '''Threshold''' !! '''Location'''
|-
|  style="text-align: left;"|Motion Sickness Incidence (MSI)
|  style="text-align: center;"|20%
|  style="text-align: center;"|Any
|-
|  style="text-align: left;"|Motion Induced Interruption (MII)
|  style="text-align: center;"|1/m
|  style="text-align: center;"|Any
|-
|  style="text-align: left;"|Roll
|  style="text-align: center;"|4º
|  style="text-align: center;"|-
|-
|  style="text-align: left;"|Pitch
|  style="text-align: center;"|1.5º
|  style="text-align: center;"|-
|-
|  style="text-align: left;"|Vertical acceleration (RMS)
|  style="text-align: center;"|0.2 g
|  style="text-align: center;"|Bridge
|-
|  style="text-align: left;"|Horizontal acceleration (RMS)
|  style="text-align: center;"|0.1 g
|  style="text-align: center;"|Bridge
|}


<div class="center" style="font-size: 75%;">'''Table 5'''. NORDFORSK 1987 [22]</div>

<div id='_tabla4'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
!style="text-align: left;"| '''General criteria for ship operability''' !! '''Merchant''' !! '''Naval''' !! '''Fast ship'''
|-
|  rowspan='2' style="text-align: left;"|Vertical acceleration RMS at bow
|  style="text-align: center;"|0.275 g (<math display="inline">L \le 100</math> m)
|  rowspan='2' style="text-align: center;"|0.275 g
|  rowspan='2' style="text-align: center;"|0.65 g
|-
|  style="text-align: center;"|0.050 g (<math display="inline">L \le 330</math> m)
|-
|  style="text-align: left;"|Vertical acceleration RMS at bridge
|  style="text-align: center;"|0.15 g
|  style="text-align: center;"|0.20 g
|  style="text-align: center;"|0.275 g
|-
|  style="text-align: left;"|Horizontal acceleration RMS in bridge
|  style="text-align: center;"|0.12 g
|  style="text-align: center;"|0.10 g
|  style="text-align: center;"|0.10 g
|-
|  style="text-align: left;"|Roll
|  style="text-align: center;"|6.0 deg
|  style="text-align: center;"|4.0 deg
|  style="text-align: center;"|4.0 deg
|-
|  rowspan='2' style="text-align: left;"|Slamming probability
|  style="text-align: center;"|0.03 (<math display="inline">L \le 100</math> m)
|  rowspan='2' style="text-align: center;"|0.03
|  rowspan='2' style="text-align: center;"|0.03
|-
|  style="text-align: center;"|0.01 (<math display="inline">L \le 330</math> m)
|-
|  style="text-align: left;"|Greenwater probability
|  style="text-align: center;"|0.05
|  style="text-align: center;"|0.05
|  style="text-align: center;"|0.05
|-
!  style="text-align: left;"|'''Acceleration and roll angle criteria''' !! '''Vertical acceleration''' !! '''Horizontal acceleration''' !! '''Roll angle'''
|-
|  style="text-align: left;"|Light manual work (RMS)
|  style="text-align: center;"|0.20 g
|  style="text-align: center;"|0.10 g
|  style="text-align: center;"|6.0 deg
|-
|  style="text-align: left;"|Heavy manual work (RMS)
|  style="text-align: center;"|0.15 g
|  style="text-align: center;"|0.07 g
|  style="text-align: center;"|4.0 deg
|-
|  style="text-align: left;"|Intellectual work (RMS)
|  style="text-align: center;"|0.05 g
|  style="text-align: center;"|0.05 g
|  style="text-align: center;"|3.0 deg
|-
|  style="text-align: left;"|Passengers (RMS)
|  style="text-align: center;"|0.05 g
|  style="text-align: center;"|0.04 g
|  style="text-align: center;"|2.5 deg
|-
|  style="text-align: left;"|Cruise line (RMS)
|  style="text-align: center;"|0.02 g
|  style="text-align: center;"|0.03 g
|  style="text-align: center;"|2.0 deg
|}


<div class="center" style="font-size: 75%;">'''Table 6'''. USCG [19]</div>

<div id='_tabla6'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
!'''Criterium''' !! '''Threshold''' !! '''Location'''
|-
|  style="text-align: center;"|Vertical acceleration 
|  style="text-align: center;"|0.4 g
|  style="text-align: center;"|Bow
|-
|  style="text-align: center;"|Vertical acceleration
|  style="text-align: center;"|0.2 g
|  style="text-align: center;"|Bridge
|-
|  style="text-align: center;"|Horizontal acceleration
|  style="text-align: center;"|0.2 g
|  style="text-align: center;"|Bridge
|-
|  style="text-align: center;"|Motion Sickness Incidence (MSI)
|  style="text-align: center;"|5 % under 30 min. exposition
|  style="text-align: center;"|Any
|-
|  style="text-align: center;"|Motion Induced Interruption (MII)
|  style="text-align: center;"|2.1 stops per minute
|  style="text-align: center;"|Any
|-
|  style="text-align: center;"|Roll angle
|  style="text-align: center;"|8.0 º
|  style="text-align: center;"|-
|-
|  style="text-align: center;"|Pitch angle
|  style="text-align: center;"|3.0 º
|  style="text-align: center;"|-
|}


In this work for the optimization analysis considering DSC, the NORDFORSK 1987 criteria for heavy manual work (recommended for fishing ships) have been adopted because those are the most demanding. Along with the NORDFORSK criteria, we have included the STANAG maximum RMS pitch angle in order to also reduce the pitch movement. The KPIs are then computed at three different locations: point 1, bow; point 2, stern-starboard; and point 3, navigation bridge ([[#_figura8|Figure 8]]). [[#_tabla7|Table 7]] shows the dimensionless coordinates of the three locations to be analysed.

<div class="center" style="font-size: 75%;">'''Table 7'''. Dimensionless coordinates of analysis points</div>

<div id='_tabla7'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! !! <math display="inline">{\bf X}/L_w</math> !! <math display="inline">{\bf Y}/B </math>  !! '''Z/T'''
|-
|  style="text-align: center;"|'''Bow'''
|  style="text-align: center;"|0.923
|  style="text-align: right;"|0.000
|  style="text-align: right;"|2.769
|-
|  style="text-align: center;"|'''Stern-board'''
|  style="text-align: center;"|0.024
|  style="text-align: right;"|0.460
|  style="text-align: right;"|2.469
|-
|  style="text-align: center;"|'''Navigation bridge'''
|  style="text-align: center;"|0.701
|  style="text-align: right;"|0.000
|  style="text-align: right;"|2.783
|}


<span id='_figura8'></span><div id="_figura8" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Review_551151435449-image11.png|408px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Figure 8. 3D view of the locations to compute the KPIs on the Ciudad de Cartagena.''</span></div>

Next the KPI are described. The MSI KPI measures the percentage of sickness after two hours navigation. This is evaluated as:

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>MSI\, \left( \%\right) =100\times \Phi \left\{ \frac{\mathrm{log}\,\left| \frac{{\ddot{\eta }}_{33}}{g}\right| -{\mu }_{MSI}}{0.4}\right\} ,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

where:

:* <math display="inline">\Phi</math>  is the cumulative normal distribution,

:* <math display="inline">\ddot{{\eta }_{33}}</math> is the heave acceleration obtained from 0.784 <math display="inline">\sqrt{{m}_{4}}</math>,

:* <math display="inline">{m}_{i}</math> is the i-th spectral moment,

:* g is the gravity acceleration,

:* <math display="inline">{\mu }_{MSI}=0.654+3.697\, (\mathrm{log}\,\left( {f}_{e}\right) +</math><math>2.32\, {\left( \mathrm{log}\,\left( {f}_{e}\right) \right) }^{2}</math> is a factor given by O’Hanlon and McCauley in 1974 [24] and depends on:

::* <math display="inline">{f}_{e}={\omega }_{e}/(2\pi )</math>,

::* <math display="inline">{w}_{e}=\, \sqrt{\frac{{m}_{4}}{{m}_{2}}}</math>.

The vertical and horizontal RMS accelerations are given by:

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{a}_{v}^{RMS}=\sqrt{{m}_{{4}_{33}}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}
{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{a}_{h}^{RMS}=\sqrt{{m}_{{4}_{22}}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

And maximum and significant roll and pitch angles are given by Eqs. 4 and 5 respectively.

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{X}_{max}=\sqrt{2\mathrm{ln}\,\left( N\right) }\, \sqrt{{m}_{0}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}
{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{X}_{S}=2\sqrt{{m}_{0}}, </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}
{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>N=\frac{t}{2\pi \, \sqrt{\frac{{m}_{0}}{{m}_{2}}}},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

where m<sub>0</sub>, m<sub>2</sub>, m<sub>4</sub> are the zero, second, and fourth spectral moments.

==4.3. Analysis==

Before optimizing the operational parameters, it will be verified whether the selected criteria are fulfilled by the Ciudad de Cartagena. Then, geometrical transformation of the original hull will be performed to improve the seakeeping performance. [[#_tabla8|Table 8]] provide the ranges of variation of the form coefficients for the geometrical transformations (<math display="inline">\pm</math> 10%, of initial values), where L is the ship length, B is the breadth, T is the draft, Vol is the displacement, A<sub>F</sub> is the flotation area, A<sub>M</sub> is the midship area, A<sub>C</sub> is the amidship area, and X<sub>b</sub> and Z<sub>b</sub> are the horizontal and vertical position of the buoyancy centre respectively. Z<sub>b</sub> is measured from free surface.

For each resulting combination, dimensionless form coefficients will be obtained, discarding unrealistic combinations such as those with block coefficient larger than one. These dimensionless hull form coefficients are the inputs for the ANN developed in [15]. Longitudinal and transversal distance from the gravity centre to the metacentric centre (<math display="inline">G{M}_{L}</math> and <math display="inline">G{M}_{T}</math>) are assumed to remain unchanged respect to the ''Ciudad de Cartagena''.

<span id='_tabla8'></span><div id="_tabla8" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Table 8. Ranges for geometric transformations.''</span></div>

{| style="width: 90%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"| 
|  style="border: 1pt solid black;text-align: center;"|'''L (m)'''
|  style="border: 1pt solid black;text-align: center;"|'''B (m)'''
|  style="border: 1pt solid black;text-align: center;"|'''T (m)'''
|  style="border: 1pt solid black;text-align: center;"|'''Vol (m<sup>3</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|'''A<sub>F</sub> (m<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|'''A<sub>M</sub> (m<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|'''AC (m<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|'''X<sub>B</sub> (m)'''
|  style="border: 1pt solid black;text-align: center;"|'''Z<sub>B</sub> (m)'''
|-
|  style="border: 1pt solid black;text-align: center;"|'''Initial value'''
|  style="border: 1pt solid black;text-align: center;"|21.650
|  style="border: 1pt solid black;text-align: center;"|6.661
|  style="border: 1pt solid black;text-align: center;"|2.980
|  style="border: 1pt solid black;text-align: center;"|203.024
|  style="border: 1pt solid black;text-align: center;"|122.922
|  style="border: 1pt solid black;text-align: center;"|14.879
|  style="border: 1pt solid black;text-align: center;"|56.239
|  style="border: 1pt solid black;text-align: center;"|10.975
|  style="border: 1pt solid black;text-align: center;"|-1.011
|-
|  style="border: 1pt solid black;text-align: center;"|'''Minimum value'''
|  style="border: 1pt solid black;text-align: center;"|19.485
|  style="border: 1pt solid black;text-align: center;"|5.995
|  style="border: 1pt solid black;text-align: center;"|2.682
|  style="border: 1pt solid black;text-align: center;"|182.722
|  style="border: 1pt solid black;text-align: center;"|110.630
|  style="border: 1pt solid black;text-align: center;"|13.391
|  style="border: 1pt solid black;text-align: center;"|50.615
|  style="border: 1pt solid black;text-align: center;"|9.878
|  style="border: 1pt solid black;text-align: center;"|-0.910
|-
|  style="border: 1pt solid black;text-align: center;"|'''Maximum value'''
|  style="border: 1pt solid black;text-align: center;"|23.815
|  style="border: 1pt solid black;text-align: center;"|7.327
|  style="border: 1pt solid black;text-align: center;"|3.278
|  style="border: 1pt solid black;text-align: center;"|223.326
|  style="border: 1pt solid black;text-align: center;"|135.215
|  style="border: 1pt solid black;text-align: center;"|16.367
|  style="border: 1pt solid black;text-align: center;"|61.862
|  style="border: 1pt solid black;text-align: center;"|12.073
|  style="border: 1pt solid black;text-align: center;"|-1.112
|}


For each geometrical factors, equidistant values are obtained within the range of variation. And for each of them the KPI is obtained. Then the KPI evaluation process consists of the following steps:

:* Step 1: Generate analysis cases within the range of parameter variations given in [[#_tabla8|Table 8]].

:* Step 2: Obtain RAOs using the ANNs developed in [15].

:* Step 3: Obtain spectral movements for each new ship: got the operational areas defined in [[#_tabla3|Table 3]], and for 7 wave directions covering from following to head waves.

:* Step 4: Obtain the maximum value of each KPI evaluated at each analysis point given in [[#_tabla7|Table 7]], considering the two operational areas and wave directions every 30 degrees, Eqs 7 - 11. Where Med. and Atl. refer to the Mediterranean Sea and the South Atlantic, respectively, and pt. refers to each of the three analysis points on the vessel.

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>MS{I}_{max}=\mathrm{max}\,\, \left\{ {\left( MS{I}_{pt1},\, MS{I}_{pt2},\, MS{I}_{pt3}\right) }_{Med},{\left( MS{I}_{pt1},\, MS{I}_{pt2},\, MS{I}_{pt3}\right) }_{Atl}\, \right\},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{{a}_{v}^{RMS}}_{max}=\mathrm{max}\,\, \left\{ {\left( {{a}_{v}^{RMS}}_{pt1},\, {{a}_{v}^{RMS}}_{pt2},\, {{a}_{v}^{RMS}}_{pt3}\right) }_{Med},{\left( {{a}_{v}^{RMS}}_{pt1},\, {{a}_{v}^{RMS}}_{pt2},\, {{a}_{v}^{RMS}}_{pt3}\right) }_{Atl}\, \right\},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{{a}_{h}^{RMS}}_{max}=\mathrm{max}\,\, \left\{ {\left( {{a}_{h}^{RMS}}_{pt1},\, {{a}_{h}^{RMS}}_{pt2},\, {{a}_{h}^{RMS}}_{pt3}\right) }_{Med},{\left( {{a}_{h}^{RMS}}_{pt1},\, {{a}_{h}^{RMS}}_{pt2},\, {{a}_{h}^{RMS}}_{pt3}\right) }_{Atl}\, \right\}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{Roll}_{max}=\mathrm{max}\,\, \left\{ {\left( {Roll}_{pt1},\, {Roll}_{pt2},\, {Roll}_{pt3}\right) }_{Med},{\left( {Roll}_{pt1},\, {Roll}_{pt2},\, {Roll}_{pt3}\right) }_{Atl}\, \right\},</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;"
|-
|
{| style="text-align: left; margin:auto;width: 100%"
|-
| style="text-align: center;" | <math>{Pitch}_{max}=\mathrm{max}\,\, \left\{ {\left( {Pitch}_{pt1},\, {Pitch}_{pt2},\, {Pitch}_{pt3}\right) }_{Med},{\left( {Pitch}_{pt1},\, {Pitch}_{pt2},\, {Pitch}_{pt3}\right) }_{Atl}\, \right\}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

The optimal geometry is obtained after a ship selection process. This is as follows:

:1. Select all ships fulfilling all stablished criteria.

:2. In case there is none, select all ships improving the maximum value of each KPI compared to the initial vessel.

:3. If none, select all ships improving the higher number of KPIs respect to the initial ship.

:4. After the previous process, all selected ships will be ranked based on each KPI, from lowest to highest.

:5. For each ship, the worst ranked position out of all KPI ranks will be taken, the candidate ships are ordered based on their position in the previously selected ranking.

:6. Finally, the ship with the lowest rank is selected.

=5. RESULTS=

First, it has been verified the fulfilment of the NORDFORSK criteria for the original ship under the wave conditions provided in [[#_tabla3|Table 3]]. [[#_tabla9|Table 9]] provide the maximum values at the three analysis points of the ship. It is observed how the initial ship does not fulfil the selected requirements.

<span id='_tabla9'></span><div id="_tabla9" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Table 9. Maximum KPIs for Ciudad de Cartagena.''</span></div>

{| style="width: 78%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-right: 1pt solid black;text-align: center;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|'''Max.'''
|  rowspan='2' style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|'''Max.'''
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|'''Bow'''
|  style="border: 1pt solid black;text-align: center;"|'''Stern-board'''
|  style="border: 1pt solid black;text-align: center;"|'''Nav. bridge'''
|-
|  style="border: 1pt solid black;text-align: center;"|'''MSI'''
|  style="border: 1pt solid black;text-align: center;"|34.863
|  style="border: 1pt solid black;text-align: center;"|28.146
|  style="border: 1pt solid black;text-align: center;"|25.981
|  style="border: 1pt solid black;text-align: center;"|34.863
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Horizontal acceleration(m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|7.419
|  style="border: 1pt solid black;text-align: center;"|5.740
|  style="border: 1pt solid black;text-align: center;"|7.243
|  style="border: 1pt solid black;text-align: center;"|7.419
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Vertical acceleration (m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|1.433
|  style="border: 1pt solid black;text-align: center;"|1.121
|  style="border: 1pt solid black;text-align: center;"|1.134
|  style="border: 1pt solid black;text-align: center;"|1.433
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Roll angle (º)'''
|  colspan='4'  style="border: 1pt solid black;text-align: center;"|54.768
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Pitch angle (º)'''
|  colspan='4'  style="border: 1pt solid black;text-align: center;"|15.288
|}


==5.1. KPI evaluation and ship selection process==

For each geometrical factors, four equidistant values are obtained within the range of variation. Then, 262,144 geometrical variations are considered. After discarding unrealistic combinations, the resulting dataset to be analysed is reduced to 190,080 combinations. Then the KPI are evaluated at the three locations, and for the seastates and wave directions defined for the operational areas. This results in 1,140,480 case scenarios simulated.

After obtaining all the KPIs for the 190,080 candidates, none of them fulfil all the four selected requirements. Then, it is found that 3,858 ships improved all maximum KPIs respect to the initial ship. Following the ranking method described for the ship selection process, it is found the best ship. [[#_tabla10|Table 10]] compares the geometrical particulars of the selected ship respect to the ''Ciudad de Cartagena''. And [[#_tabla11|Table 11]] provides the KPIs for the selected ship.

<span id='_tabla10'></span><div id="_tabla10" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Table 10. Initial and best ships particulars.''</span></div>

{| style="width: 60%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|+
| style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;" |
| style="border: 1pt solid black;text-align: center;" |'''Ciudad de Cartagena'''
| style="border: 1pt solid black;text-align: center;" |'''Best ship'''
|-
| style="border: 1pt solid black;text-align: center;" |'''LW (m)'''
| style="border: 1pt solid black;text-align: center;" |21.650
| style="border: 1pt solid black;text-align: center;" |20.928
|-
| style="border: 1pt solid black;text-align: center;" |'''B (m)'''
| style="border: 1pt solid black;text-align: center;" |6.661
| style="border: 1pt solid black;text-align: center;" |5.995
|-
| style="border: 1pt solid black;text-align: center;" |'''T (m)'''
| style="border: 1pt solid black;text-align: center;" |2.980
| style="border: 1pt solid black;text-align: center;" |3.079
|-
| style="border: 1pt solid black;text-align: center;" |'''Displacement (t)'''
| style="border: 1pt solid black;text-align: center;" |208
| style="border: 1pt solid black;text-align: center;" |229
|-
| style="border: 1pt solid black;text-align: center;" |'''CB'''
| style="border: 1pt solid black;text-align: center;" |0.472
| style="border: 1pt solid black;text-align: center;" |0.578
|-
| style="border: 1pt solid black;text-align: center;" |'''CF'''
| style="border: 1pt solid black;text-align: center;" |0.852
| style="border: 1pt solid black;text-align: center;" |0.947
|-
| style="border: 1pt solid black;text-align: center;" |'''CM'''
| style="border: 1pt solid black;text-align: center;" |0.750
| style="border: 1pt solid black;text-align: center;" |0.725
|-
| style="border: 1pt solid black;text-align: center;" |'''CP'''
| style="border: 1pt solid black;text-align: center;" |0.630
| style="border: 1pt solid black;text-align: center;" |0.797
|-
| style="border: 1pt solid black;text-align: center;" |'''XB / LW'''
| style="border: 1pt solid black;text-align: center;" |0.507
| style="border: 1pt solid black;text-align: center;" |0.507
|-
| style="border: 1pt solid black;text-align: center;" |'''KB / T'''
| style="border: 1pt solid black;text-align: center;" |0.661
| style="border: 1pt solid black;text-align: center;" |0.639
|}


<span id='_tabla11'></span><div id="_tabla11" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Table 11. KPI values for selected best ship''</span></div>

{| style="width: 78%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-right: 1pt solid black;text-align: center;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|'''Max.'''
|  rowspan='2' style="border-top: 1pt solid black;border-left: 1pt solid black;border-right: 1pt solid black;text-align: center;"|'''Max.'''
|-
|  style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;"|'''Bow'''
|  style="border: 1pt solid black;text-align: center;"|'''Stern-board'''
|  style="border: 1pt solid black;text-align: center;"|'''Nav. bridge'''
|-
|  style="border: 1pt solid black;text-align: center;"|'''MSI'''
|  style="border: 1pt solid black;text-align: center;"|29.464
|  style="border: 1pt solid black;text-align: center;"|28.987
|  style="border: 1pt solid black;text-align: center;"|22.164
|  style="border: 1pt solid black;text-align: center;"|29.464
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Horizontal acceleration(m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|3.902
|  style="border: 1pt solid black;text-align: center;"|3.469
|  style="border: 1pt solid black;text-align: center;"|3.915
|  style="border: 1pt solid black;text-align: center;"|3.915
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Vertical acceleration (m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|1.198
|  style="border: 1pt solid black;text-align: center;"|1.141
|  style="border: 1pt solid black;text-align: center;"|1.005
|  style="border: 1pt solid black;text-align: center;"|1.198
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Roll angle (º)'''
|  colspan='4'  style="border: 1pt solid black;text-align: center;"|28.352
|-
|  style="border: 1pt solid black;text-align: center;"|'''Max. Pitch angle (º)'''
|  colspan='4'  style="border: 1pt solid black;text-align: center;"|14.901
|}

==5.2. Discussion of results==

[[#_tabla12|Table 12]] compares the maximum KPI values for the Ciudad de Cartagena and the best selected ship. It is observed that a considerable improvement of the seakeeping performance is achieved.

<span id='_tabla12'></span><div id="_tabla12" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">''Table 12. KPI values comparison for selected best ship and the Ciudad de Cartagena''</span></div>

{| style="width: 70%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  colspan="2"; style="border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"| 
|  style="border: 1pt solid black;text-align: center;"|'''Ciudad de Cartagena'''
|  style="border: 1pt solid black;text-align: center;"|'''Best ship'''
|-
|  colspan="2"; style="border: 1pt solid black;text-align: center;"|'''MSI'''
|  style="border: 1pt solid black;text-align: center;"|34.863
|  style="border: 1pt solid black;text-align: center;"|29.464
|-
|  colspan="2"; style="border: 1pt solid black;text-align: center;"|'''Max. Horizontal acceleration(m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|7.419
|  style="border: 1pt solid black;text-align: center;"|3.915
|-
|  colspan="2"; style="border: 1pt solid black;text-align: center;"|'''Max. Vertical acceleration (m/s<sup>2</sup>)'''
|  style="border: 1pt solid black;text-align: center;"|1.433
|  style="border: 1pt solid black;text-align: center;"|1.198
|-
|  colspan="2"; style="border: 1pt solid black;text-align: center;"|'''Max. Roll angle (º)'''
|  style="border: 1pt solid black;text-align: center;"|54.768
|  style="border: 1pt solid black;text-align: center;"|28.352
|-
|  colspan="2"; style="border: 1pt solid black;text-align: center;"|'''Max. Pitch angle (º)'''
|  style="border: 1pt solid black;text-align: center;"|15.288
|  style="border: 1pt solid black;text-align: center;"|14.901
|}


The selected recommendations are appropriate to assess seakeeping performance since they stablished thresholds based on the work to be carried out on board. In the case study of this work, a fishing ship, those of heavy manual work are taken as appropriated.

=6. CONCLUSIONS=

A wise selection of the form coefficients in the early design stage of ships is key to ensure the seakeeping performance and safety of the ship, for instance, in dead ship conditions. In this context, the ANN developed in [15] can be used for fast evaluation of KPIs and hull form optimization, providing ship designer with a powerful tool to this end.

It has been shown how the ANN from [15] provides the capability of testing a very large number of hulls in a negligible time. In this work, approximately 200,000 ships have been evaluated under two operative conditions and three points of analyses (1.14 x 10<sup>6</sup> cases) in 172 hours.

Out of the analyses carried out, it is observed that the most critical KPIs are the maximum roll angle and RMS horizontal acceleration. And the ship selected out of the optimization process improved the KPIs respect to the initial ship.

=7. REFERENCES=

[1] CIAIM (2019). Informe de 2019 de la Comisión de Investigación de Accidentes e Incidentes Marítimos, ministerio de transportes movilidad y agenda urbana. Gobierno de España.

[2] IMO SDC 1/INF.8. (2013). Information collected by the correspondence group on intact stability regarding the second generation intact stability criteria development. Submitted by Japan. 13 November 2013.

[3] Petacco, N. and Gualeni, P. (2020). IMO Second Generation Intact Stability Criteria: General Overview and Focus on Operational Measures. Journal of Marine Science and Engineering 8(7):494. [https://doi.org/10.3390/jmse8070494 https://doi.org/10.3390/jmse8070494]

[4] Mata-Álvarez-Santullanol, F. and Pérez-Rojas, L. (2015). Application of IMO Second Generation Intact Stability Criteria for Dead Ship Condition to Small Fishing Vessels. Proceedings of the 12th International Conference on the Stability of Ships and Ocean Vehicles, 14-19 June 2015, Glasgow.

[5] Choi, J.; Jensen, J.; Kristensen, H.; Nielsen, U. and Erichsen, H. (2017). Intact Stability Analysis of Dead Ship Conditions using FORM. J. Ship. Res. 61:167–176. [https://doi.org/10.5957/JOSR.170005. https://doi.org/10.5957/JOSR.170005.]

[6] Kubo, T. Umeda, N. Izawa, S. and Matsuda, A. (2019). Total Stability Failure Probability of a Ship in Beam Wind and Waves: Model Experiment and Numerical Simulation. In: Belenky, V., Spyrou, K., van Walree, F., Almeida Santos Neves, M., Umeda, N. (eds) Contemporary Ideas on Ship Stability. Fluid Mechanics and Its Applications, vol 119. Springer, Cham 2019. [https://doi.org/10.1007/978-3-030-00516-0_35 https://doi.org/10.1007/978-3-030-00516-0_35]

[7] Gu, M. Lu, J. and Wang, Th. (2015). Stability of a tumblehome hull under the dead ship condition. J Hydrodyn. Ser. B. 27(3), 452–457. [https://doi.org/10.1016/S1001-6058(15)60503-0. https://doi.org/10.1016/S1001-6058(15)60503-0.]

[8] Shin, D.M. and Chung, J. (2021) Application of dead ship condition based on IMO second-generation intact stability criteria for 13K oil chemical tanker. Ocean Engineering. 238:109776. [https://doi.org/10.1016/j.oceaneng.2021.109776. https://doi.org/10.1016/j.oceaneng.2021.109776.]

[9] Hu, L. Wu, H. Yuan, Z. Li, and W. Wang, X. (2021). Roll motion response analysis of damaged ships in beam waves. Ocean Engineering. 227:108558. [https://doi.org/10.1016/j.oceaneng.2020.108558 https://doi.org/10.1016/j.oceaneng.2020.108558]

[10] Liu, L.; Yang, Y. and Peng, T. (2022). Machine learning prediction of 6-DOF motions of KVLCC2 ship based on RC model. J. Ocean Eng. Sci. [https://doi.org/10.1016/j.joes.2022.08.004. https://doi.org/10.1016/j.joes.2022.08.004.]

[11] Taghva, H.R. Ghassemi, H. and Nowruzi, H. (2018). Seakeeping performance estimation of the container ship under irregular wave condition using artificial neural network. Am. J. Civil Eng. Archit. 6:147-153. [https://doi.org/10.12691/ajcea-6-4-3. https://doi.org/10.12691/ajcea-6-4-3.]

[12] Cepowski, T. (2008). Determination of optimum hull form for passenger car ferry with regard to its sea-keeping qualities and additional resistance in waves. Pol. Marit. Res. 15:3-11. [https://doi.org/10.2478/v10012-007-0057-2. https://doi.org/10.2478/v10012-007-0057-2.]

[13] Cepowski, T. (2009). On the modeling of car passenger ferry ship design parameters with respect to selected sea-keeping qualities and additional resistance in waves Pol. Marit. Res. 16: 3-10, [https://doi.org/10.2478/v10012-008-0026-4. https://doi.org/10.2478/v10012-008-0026-4.]

[14] Cepowski, T. (2020). The prediction of ship added resistance at the preliminary design stage by the use of an artificial neural network. Ocean Eng. 195:106657. [https://doi.org/10.1016/j.oceaneng.2019.106657. https://doi.org/10.1016/j.oceaneng.2019.106657.]

[15] Romero-Tello, P. Gutiérrez-Romero, J.E. and Serván-Camas, B. (2022). Prediction of seakeeping in the early stage of conventional monohull vessels design using artificial neural network. Journal of Ocean Engineering and Science. [https://doi.org/10.1016/j.joes.2022.06.033. https://doi.org/10.1016/j.joes.2022.06.033.]

[16] Delhommeau, G. (1993). Numerical Simulation of Hydrodynamics: Ships and Offshore Structures. Seakeeping codes AQUADYN and AQUAPLUS. Ecole Centrale de Nantes.

[17] Rumawas, V.; Asbjørnslett, B. E. and Klöckner, C. A. (2017). Human Factors Evaluation in Ship Design: A Case Study on Offshore Supply Vessels in the Norwegian Sea, Part II: Multivariate Analyses and Structural Modelling. Naval Engineers Journal 129 (1):45.

[18] Pattison, J. H. and Sheridan, D. J. (2004). Human Performance Factors and Measures in Hull Form Selection. NATO unclassified/unlimited, RTO-MP-AVT-110, Meeting Proceedings-Paper 34 date 10/04/2004

[19] Ghaemi, M.H. and Olszewski, H. (2017). Total Ship Operability-Review, Concept And Criteria. Polish Maritime Research. Special Issue 2017 S1 (93). 24:74-81. [https://doi.org/10.1515/pomr-2017-0024. https://doi.org/10.1515/pomr-2017-0024.]

[20] Stevens, S.C. and Parsons, M.G. (2002). Effects of Motion at Sea on Crew Performance: A Survey. Marine Technology. 39(1):29–47.

[21] NATO. (2000). Standardization Agreement (STANAG): Subject: Common Procedures for Seakeeping in the Ship Design Process, NATO, Military Agency for Standardization.

[22] NORDFORSK. (1987). Assessment of Ship Performance in a Seaway. Copenhagen, Nordforsk.

[23] O'Hanlon, J. and McCauley, M. (1974). Motion Sickness Incidence as a Function of Vertical Sinusoidal Motion. Aerospace medicine. Vol. 45; pp. 366-369