Application of the conditional optimization method in the problem of vessel stormy sailing Serhii Zinchenko1, Oleh Tovstokoryi1, Vadym Mateichuk 1, Pavlo Nosov1, Ihor Popovych2 and Kostiantyn Kyrychenko1 1 Kherson state maritime academy, 20, Ushakova ave., Kherson, 73000, Ukraine 2 Kherson state university, 27, Universytetska str., Kherson, 73003, Ukraine Abstract Stormy sailing conditions are among the most extreme and dangerous. The rocking of the ship in the roll and trim channels greatly exhausts the crew, which leads to an increase in the number of errors when making management decisions. The situation worsens also due to the fact that during a storm such dangerous phenomena as harmonic and parametric resonances, a decrease in stability on passing waves, broaching, impacts of group waves in the stern, which can lead to the , hull destruction, overloads of the power plant and occurrence of blackouts. The most radical way to combat such dangers is the automation of control processes. The authors have developed a method of automatic and optimal control of the ship's movement in storm conditions, which allows to avoid the occurrence and development of dangerous phenomena that lead to the capsize of the ship. The obtained result is explained by: the use of an on-board computer in the ship's motion control system; constant measurement of ship movement parameters and waves; finding, at each step of the on-board computer, safe and optimal values of the ship's course and speed by solving the optimization problem; taking into account when solving the optimization problem linear and non-linear constraints of the type of inequalities that define dangerous areas; maintaining safe and optimal movement parameters using the automatic control system. The theoretical significance of the obtained result lies in the development of a method of automatic and optimal control of the ship's movement in conditions of stormy sailing. The practical significance of the obtained results consists in: verification of the developed method by mathematical modeling; the possibility of using the method in the automatic control system, which allows to automate and optimize the processes of controlling a ship in a storm, reduce the influence of the human factor on control processes, crew fatigue, risks of losing the ship and cargo, and generally increase the safety of shipping. Keywords Navigation safety, human factor, intelligent transport systems, automated system, resonance zone, stormy sailing 1 1. Introduction Steering in a storm is an important stage in sailing a ship. capsizing caused by harmonic resonance, parametric resonance of the first and second type, loss of stability in following seas, broaching and impacts of group waves in the stern of the ship. Also, during a storm, the ship is subjected to large loads, which can significantly increase with improper control and lead to the destruction of the hull. To facilitate the task of steering a ship in a storm, a number of scientists have suggested using storm diagrams. The universal Remes diagram allows to determine dangerous areas of harmonic resonance and choose safe parameters of storming (course and speed of the vessel). At various times, the International Maritime Organization has also developed guidelines for safe sailing in storms. This is "Resolution A.562(14), adopted on 20 November 1985, Recommendation on a severe wind and rolling criterion (weather criterion) for the intact stability of passenger and cargo ships of 24 meters in length and over, MSC/Circ.707 from 19 October 1995. Guidance to the master for avoiding dangerous situations in following and quartering seas (MSC/Circ.707 from 19 October 1995), MSC.1/Circ.1228 from 11 January 2007. Revised guidance to the master for avoiding dangerous situations in adverse weather and sea conditions. Storm diagrams and IMO guidelines involve visual ICST-2024: Information Control Systems & Technologies, September 23-25, 2023, Odesa, Ukraine. srz56@ukr.net (S. Zinchenko); otovstokory@gmail.com (O. Tovstokoryi); mateichykv@gmail.com (V. Mateichuk); pason@ukr.net (P. Nosov); Ihorpopovych999@gmail.com (I. Popovych); kvklecturer@gmail.com (K. Kyrychenko); 0000-0001-5012-5029 (S. Zinchenko); 0000-0003-3048-0028 (O. Tovstokoryi); 0000-0001-9328-0651 (V. Mateichuk); 0000-0002-5067-9766 (P. Nosov); 0000-0002-1663-111X (I. Popovych); 0000-0002-0974-6904 (K. Kyrychenko) © 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR ceur-ws.org Workshop ISSN 1613-0073 Proceedings determination of wave parameters and subsequent manual calculations, using formulas and graphs, which gives low accuracy and contains an element of human factors. To reduce sway, stabilizers of sway built on various physical principles, vane thrusters, special forms of the ship's hull, etc. are also used. Recently, automated navigation control systems have been developed: HULLMOS navigation parameters monitoring system, OCTOPUS-DSS decision support system, IMDSS integrated marine decision support system, VOSS navigation optimization assessment system, fiber optic system "SENSFIB" and others, which automatically measure the parameters of the own vessel's movement and disturbances, perform the necessary calculations for assessing the seaworthiness of the vessel. At the same time, automated navigation control systems involve the presence of a person in the control circuit and associated decision-making delays and the human factor. The authors of this article believe that the best solution in this case is the use of automatic storming modules in automated systems. An example of automatic module in automated systems is autopilot, which was introduced in the last century and is successfully used on modern vessels. The captain only needs to make a decision to use the automatic module, and then observe its operation. The use of automatic control modules to solve various problems was considered by us earlier, in particular: ot point [18], automatic reset of kinetic energy in case of an inevitable collision [19], automatic control of redundant structures of executive devices [20], a wave [21], devising an approach to safety management of vessel control through the identification of navigator's state [22] and other. Based on the above, a very important issue is the development of the automatic control methods, which would ensure safe sailing in stormy conditions, the implementation of these methods in the automatic control modules of the ship's control system. 2. Related works Many works are devoted to the study of the processes of controlling a ship in a storm. In particular, work [1] considered the optimization of route planning in order to reduce the probability of the ship encountering severe weather conditions and to minimize the accumulation of fatigue in the ship's structure. Storm profiles were defined as a set of waves, the parameters of which were obtained on the basis of full-scale measurements for two years of flights. The advantages of using route optimization are considered on the example of a container ship. The work [2] describes a method that allows you to prepare polar maps for various assumed sea states and load conditions in a relatively short time, without long-term modeling of irregular wave conditions. If ships are equipped with a wave radar and have access to a weather service, the polar charts will correspond to the actual sea state, allowing them to be used in practice to form a route taking into account the weather, including the dangers of parametric rolling. The work [3] investigated the dynamics, controllability, and stability of oscillations in a roll channel on a regular wave. The vessel was modeled as a classical Helmholtz-Duffing oscillator with highly nonlinear asymmetric recovery moment characteristics. The method of incremental harmonic balance was used to obtain frequency characteristics. The effect of gain and time delay in the control circuit of the fundamental and subharmonic responses is investigated. Steady, periodic, and stationary solutions obtained by the incremental harmonic balance method are verified by numerical simulations. The solutions are supplemented with phase portrait, Poincaré map, time chronology and Fourier spectrum for better understanding. In work [4] the issue of calculating the damping moment in the roll channel was studied. An accurate calculation of the damping moment improves the prediction of oscillations. There are several methods for calculating the damping moment, but most of them do not take into account the effects of viscosity. In the study, the authors used the method of computational hydrodynamics to determine the characteristics of oscillation and damping coefficients in the roll channel for various conditions. The influence of scale effects on vibration characteristics and damping coefficients was studied. At various times, the International Maritime Organization has also developed guidelines for safe sailing in storms. In the latest document MSC.1/Circ.1228, the requirements for avoiding the dangers associated with capsizing of the vessel from: synchronous resonance, parametric resonance, reduction of stability in following seas, broaching and encounter wave grouping phenomena. In addition to these dangers, there are also dangers associated with the destruction of the vessel's hull. Works [5-7] are devoted to the study of issues of hull strength. Storm diagrams and IMO guidelines involve visual determination of wave parameters and subsequent manual calculations, using formulas and graphs, which gives low accuracy and contains an element of human factors. To reduce wobble, sway stabilizers [8], hydrodynamic stabilizers [9], gravity stabilizers [10, 11], gyroscopic stabilizers, vane thrusters, special forms of the vessel's hull [12], etc. are also used. Work [13] gives recommendations to shipmasters when preparing a ship for stormy sailing in terms of preparing the steering system, engine room, water tight doors and other important systems. In the work [14], chapter 7 "Ship Handling in Rough Sea", the author describes the dangers of sailing in counter and following seas. In counter seas is this propeller racing; speed reduction and the torque rich effect on the engine; shipping seas; slamming phenomenon; countermeasures for rough weather in countering seas. In following seas it are encounter wave grouping phenomena; parametric rolling phenomena; reduction of stability; surf-riding (broaching-to) phenomena; countermeasures for rough weather in following seas. The given material is accompanied by illustrations, tables and graphs. In doctoral thesis [15], in the storm sailing section, the author investigates the forces and moments acting on the ship's hull, conducts mathematical modeling of the processes of the ship's angular motion in the roll channel and the stability of the ship during waves. The author suggested using different models, depending on the navigation task to be solved. The linear 6-degree model of the ship's angular motion in the roll channel, which can be used to calculate the amplitude-frequency characteristics, is advisable to use for evaluating the roll amplitude, optimizing the ship's route in channels and narrows, and planning operations at shallow depths. It is advisable to use the nonlinear 6-degree model of the ship in studies of stability, parametric resonance, broaching, etc. The article [16] discusses the issue of parametric resonance in a roll channel on a regular wave. The research was carried out experimentally on a typical Norwegian fishing vessel with a blunt hull and a small ratio of length to width, as well as numerical simulation. Nonlinearities in Froude-Krylov loads and recoverable loads were taken into account by integrating the pressure over the instantaneous wetted surface of the hull. It was found that near the instability limit of the Mathieu diagram, experimental and numerical predictions regarding the appearance of parametric roll differ from each other. The instability limits in the experiment also differed from the instability limits of the Mathieu diagram. The region of instability for the 6-DOF experiments and simulations spans a wider range of frequency ratios, and the amplitude threshold of the metacentric height is found to be lower than that of the Mathieu diagram. The results also showed that the region of instability, in the presence of speed, shifts towards smaller ratios of natural and forced oscillation frequencies. The most common and most dangerous phenomenon on container ships is parametric rocking in countering and following seas, which occurs as a result of dynamic instability of the vessel. The causes of parametric roll are periodic (with a period that is twice the period of the wave) changes in the metacentric height during the movement of the ship across the waves. Parametric roll is also facilitated by special hull shapes (use of extended sides in the front of the hull, etc.). The authors of the article [17], by means of mathematical modeling in the MATLAB environment, investigated the influence of wave height, wavelength and ship speed on the parametric resonance of a container ship. The enveloping stability curve for countering and following excitation is determined. Damping effects were taken into account when estimating the range of metacentric height changes. As can be seen from the above review of literary sources, the issue of automatically controlling the movement of a ship in a storm has not been considered before, so their solution remains an urgent scientific and technical task. In this article, the authors developed a storming method that allows to automatically find safe and optimal parameters of the ship's movement (speed and course) by solving at each step of the on- board computer an optimization problem with linear and nonlinear constraints of the inequalities 3. Materials and Methods The object of research is the processes automatic stormy sailing. The subject of research is methods and models of automatic stormy sailing. The purpose of research is to develop -board computer of automatic control system, which would allow automating processes of storm sailing, reduce influence of the human factor on the control processes, reduce crew fatigue, reduce the risks of losing the ship and cargo, and generally increase navigation safety. To solve the problem, authors used numerical methods of conditional optimization to find safe and optimal parameters of the ship's movement in a storm at each step of the on-board computer; recommendations of the IMO, theoretical work and practical experience of other authors, for the formation of inequalities that define the dangerous areas of storm sailing, and are used as constraints when solving the problem of conditional optimization; methods of automatic control theory, mathematical modeling. 4. Results During stormy sailing, there is a danger of the vessel capsizing caused by synchronous resonance, parametric resonance, reduction of stability in following seas, broaching, impacts of encounter group waves in the stern of the vessel, and the danger of destruction of the vessel's hull due to exceeding the maximum permissible forces and moments. In fig. 1 shows areas of synchronous resonance GR , parametric resonance of the first type  PR1 and parametric resonance of the second type  PR 2 for the wavelength  = 90m and the period of the vessel's own oscillations TС = 12,5 s . Figure 1: Areas of synchronous and parametric resonance The task of storm sailing is to choose safe parameters of the vessel's movement - speed V * and course * , which would not belong to dangerous zones. Taking into account that in addition to the indicated areas, there are other areas, the impact of which can lead to the capsizing of the vessel or the destruction of the hull, it becomes clear that it is quite difficult to take into account all the existing dangers when manually controlling the vessel, and even in the presence of difficult sailing conditions. The authors believe the most radical way to solve this problem is the use of automatic storming modules in automated systems, the algorithms of which can find not only safe, but also optimal storming parameters. The task of the automatic storming module is constantly, with on- taking into account the permissible range of storming speeds and all areas dangerous for capsizing the vessel and destroying the hull, supporting the automatic control system of safe and optimal movement parameters. In the most general form, the task of finding safe and optimal storm parameters can be formulated as follows F (V * , q* , , TC ) = min F (V , q, , TC ), (1) V ,q  f1(V , q, , TC )  0   f 2 (V , q, , TC )  0  , (2) ...  f n (V , q, , TC )  0 where F (V , q, , TC ) is the aim function to be optimized, f j (V , q, , TC )  0, j = 1..n is the linear and non-linear constraints of inequalities type. The objective function (1) and constraints (2) depend on V and the heading angle q of the wave, as well as on the wavelength  and the period of the vessel's own oscillations TC . Of the named parameters, only two are available for control - V and the heading angle of the wave q =  − KW + 180dg (due to the course of the vessel), which we will use for conditional optimization. Other parameters on which we have no influence (wave parameters, maximum and minimum vessel speeds, period of natural oscillations of the vessel, etc.) are used as external data when solving the optimization problem. As the first limitation, we will accept a linear limitation on the vessel's speed. To maintain controllability, the vessel's speed must be greater than minimum speed Vmin and less than maximum st in storm conditions, i.e. the condition must be fulfilled speed Vmax Vmin  V  Vmax st (3) As the second limitation, we will accept the condition of finding the vessel's movement parameters (speed and course) outside the limits of synchronous resonance. It is known [1] that synchronous resonance occurs under the condition when the ratio of the vessel's own oscillations period TC to the imaginary period of the waves  rocking the vessel is within T 0,7  C  1,3 , (4)  where the imaginary period  of the waves is determined by the formula  = . (5) 1,25  + V cos q From the inequalities (4), we determine the area of the storm that is safe from synchronous resonance  TC    1,3  (6)  TC  0,7   As the third and fourth constraints, we accept the condition of finding the vessel's motion parameters outside the parametric resonance of the 1st and 2nd types. It is known that the most dangerous is parametric resonance of the 2nd type, when the period of the natural oscillations of the vessel is twice as long as the imaginary period of the wave. TC −   T , (7) 2 where T is the difference between the imaginary periods of the waves and the half-period of the vessel's own oscillations. Also, parametric resonance can occur when the period of the vessel's own oscillations coincides with the period of the waves. TC −   T (8) As the fifth limitation, we accept the condition of finding the parameters of the vessel's movement outside the region of stability reduction in following seas. The vessel loses its stability in following seas under the condition [5]  = L  V = C , (9) q = 180dg  where L , C = 1,25  is the wave speed. From the system (9), we find the region safe from the reduction of stability in following seas   − L  LLS . or . V − C  VLS     q , (10)  q − 180 LS  where LLS , VLS , q LS are the dimensions of the region of loss of stability in following seas. As the sixth limitation, we accept the condition of finding the parameters of the vessel's motion outside the broaching limits. Broaching often results in loss of control during accelerated sliding along the steep leading edge of a high wave. Broaching occurs under the condition [5] V cos(180 − q )   1,4  L (11)    135  q  225 From system (11), we find the area of storming that is safe from broaching V cos q  −1,4 L  (12) q  225dg. or. q  135dg As the seventh limitation, we accept the condition of finding the vessel's movement parameters outside the area of dangerous impacts of encounter group waves in the stern. As the seventh limitation, we accept the condition of finding the vessel's movement parameters outside the area of dangerous impacts of encounter group waves in the stern. In the accompanying disturbance, when the speed of the waves slightly exceeds the speed of the vessel, then the handling of the vessel becomes uncontrollable. Also, the action of such waves can damage the hull and the steering device. The dangerous area is determined by the V/T diagram [5], based on which we write down the system of inequalities that define the boundaries of this dangerous area. The dangerous area is determined by the V/T diagram [5], based on which we write down the system of inequalities that define the boundaries of this dangerous area.  V cos(180 − q) 1,3  2  TW (13) 135dg  q  225dg  From system (13), we find the area safe from the impact of group waves in the stern of the vessel. V cos q  −2TW . or.V cos q  −1,3TW  , (14) q  135dg . or. q  225dg where TW is the wave period. Restrictions can be extended further to determine areas of navigational hazards, damage to the vessel's hull, etc. As aim function we can select a function F = ( * −SET ) 2 → min , (15) which minimizes the deviation of the safe course  * from the given  SET , or function F = (V * −VSET ) 2 → min , (16) which minimizes the deviation of the safe speed V * from the given VSET , or any other target function that will provide the desired control quality. To find the parameters V * ,  * , which optimize the objective function (15), (16), or another and are within the safe regions (3), (6), (7), (8), (10), (12), ( 14), we apply the procedure of nonlinear optimization of the type f min con (•) MATLAB f min con (@ fun, x0, A, b, Aeq , beq , lb, ub, @ nonlcon ) (17) where @ fun ( ) is a link to the target function file, x0 = V (0), q (0) is the initial vector of parameters to be optimized, A is a matrix of the linear inequalities system, absent, b is a vector of the right part linear inequalities system, absent, Aeq is a matrix of the linear equalities system, absent, beq is a vector of the right part linear equalities system, absent, lb = − Vmin ,−  is a vector of the lower limit of the parameter area, ub = Vmax ,  is a vector of the upper limit of the parameter area, @ nonlcon is a link to a nonlinear constraint system file. The safe parameters of the vessel's movement V * , * , which are found at each step of the on- board controller by solving a conditional optimization problem with linear and nonlinear constraints of the inequalities type, are fed to the inputs of the units for determining the control of the power plant and the rudder.  V*  =  max  Vmax , (18)  * ( *  = k ( −  ) + k z + k    −  dt) where  is a deviation angle of the power plant telegraph,  max is a maximum deviation angle of the power plant telegraph, K is a current course,  z is a yaw angular rate, k  , k , k  are the gain coefficients of the PID - regulator. 4.1. Experiment To confirm the efficiency and effectiveness of the developed methods of automatic stormy sailing, mathematical modeling of stormy sailing processes in the MATLAB environment was carried out. 4.1.1. Simulation of automatic control with synchronous resonance Simulation of automatic control with synchronous resonance. The initial parameters of the vessel's movement and waves correspond to point A (synchronous resonance area), Fig. 1: the specified course of the vessel is (0) = 180dg , the specified speed of the vessel is V (0) = 4,1 m s ; wave course is KW = 0dg wave speed is C = 13,7 m s , wave length is  = 120m . Fig. 2 shows graphs of changes in time of longitudinal speed V x [ m / s ] , longitudinal displacement X g [m] , roll angular rate  x [dg / s ] , roll angle [dg ] , yaw angular rate  z [dg / s ] and yaw angle [dg ] . Figure 2: Longitudinal and angular movement parameters change in roll and yaw channel in experiment with synchronous resonance Fig. 3 shows graphs of changes in time of lateral speed V y [ m / s ] , lateral movement Yg [m] , trim rate  y [dg / s ] , trim angle [dg ] , telegraph deflection teta[dg ] and stern deflection delta[dg ] . The control system supports the set course of the vessel (t ) = 0dg for up to t = 50 s . In this area, a significant amplitude of oscillations (t ) = 50dg is observed in the roll channel. For t = 50 s , the automatic storming module is turned on, which calculates the safe course and speed of the vessel. To change the course and speed, the rudder delta (t ) and the telegraph teta (t ) are deflected. After that, a change in course and speed is observed. As can be seen from the graph (t ) , the course of the vessel begins to change from (t ) = 0dg to (t ) = −77dg , the speed of the vessel increases from V x (t ) = 4.1 m / s to V x (t ) = 5.2 m / s . Starting with t = 80 s amplitude of oscillations in the roll channel sharply decreases to (t ) = 5dg , which corresponds to the amplitude of forced oscillations. The change in lateral speed V y (t ) and lateral displacement Yg (t ) is caused by a change in the ship's course (t ) . The trim angular rate  y (t ) and the trim angle  (t ) change with the frequency of the forced oscillations and are within the limits  y (t ) = 0.5dg and  (t ) = 0.2dg . Figure 3: with synchronous resonance 4.1.2. Simulation of automatic control with parametric resonance The initial parameters of the vessel's movement and the wave correspond to point B, which is in the area of parametric resonance  PR 2 , fig. 1: the vessel's course is (0) = 0dg , the vessel's speed is V (0) = 2,5 m s ; wave course is KW = 180dg , wave speed is C = 11,9 m s . Fig. 4 shows graphs of changes in time of longitudinal speed V x [ m / s ] , longitudinal displacement X g [m] , roll angular rate  x [dg / s ] , roll angle [dg ] , yaw angular rate  z [dg / s ] and yaw angle [dg ] . Fig. 5 shows graphs of changes in time of lateral speed V y [ m / s ] , lateral movement Yg [m] , trim rate  y [dg / s] , trim angle [dg ] , telegraph deflection teta[dg ] and stern deflection delta[dg ] . The automatic storming module continuously calculates and maintains the optimal and safe parameters of the vessel's movement. For optimization, the aim function (15) was used, which ensured the minimum deviation of safe course, calculated at each step on-board controller, from set course. The calculated optimal and safe storming parameters are V * = 2.1 m / s, * = 0dg . It can be seen from the graphs V (t ), (t ) that the control system maintains a safe course with precision (t )  2dg and at the same time reduces the vessel's speed to a safe one. After reaching a safe speed for about t = 100 s , the amplitude of oscillations in the roll channel sharply decreases to the value of forced oscillations. As for the channel of lateral movement, it can be seen from the graphs that the lateral speed during the entire experiment did not exceed V y (t )  0.2m / s , and the lateral deviation did not exceed Yg (t )  6.0m . In the trim channel, the angular speed does not exceed  y (t )  0.5dg / s , and the trim angle does not exceed  (t )  0.2dg . Figure 4: Longitudinal and angular movement parameters change in roll and yaw channel in experiment with parametric resonance Figure 5: with parametric resonance The analysis of the simulation results shows that the automatic storming module allows to effectively calculate the optimal and safe parameters of the vessel's movement and maintain them during stormy sailing. Further research may be related to the extension of the conditional optimization method to other hazardous areas to avoid hull damage, bypass navigational hazards, etc. 5. Discussion A which allows to avoid the occurrence and development of dangerous phenomena, such as harmonic and parametric resonances, loss of stability in the following seas, broaching, impacts of the group The obtained result is explained by: the use of on-board controller in the ship's motion control system; constant measurement of ship movement and wave parameters; finding at each step of the on-board controller safe and optimal values of the ship's course and speed by solving the optimization problem; taking into account, when solving the optimization problem, linear and non- linear constraints of the inequalities type, which define dangerous areas; maintaining safe and optimal movement parameters by means of the automatic control system. The obtained results differ from the known solutions in that they allow to automate and optimize the control processes a ship in a storm, reduce the influence of the human factor on the, reduce the risks of capsizing the ship, save cargo and human life. The method of determining the safe and optimal parameters of the ship's movement in a storm is designed only for use in the on-board computer of the automated ship movement control system and cannot be used for manual control. The obtained results are reproducible and can be used in the development of automated systems with automatic storming modules. Further research may be related to expanding the range of possible applications of the developed method. 6. Conclusion A method of automatic stormy sailing has been developed, which allows to automatically find safe and optimal stormy sailing parameters and maintain them. The obtained results are explained by the solution at each step of the on-board controller of the optimization problem with constraints of the inequalities type, which describe the dangerous areas. The scientific novelty of the obtained results lies in the fact that, for the first time, a method of automatically avoiding the dangers of capsizing a ship during stormy sailing has been developed, unlike existing solutions, it allows to automate and optimize the processes of ship control in a storm. The practical significance of the obtained results is that: the workability and effectiveness of the method are verified by mathematical modeling in the MATLAB environment. The method can be used for the development of automatic modules for controlling the movement of a ship in stormy conditions, which will allow automating the control processes in a storm, reducing the influence of the human factor on the control processes, reducing the fatigue of the crew, reducing the risks of losing the ship and cargo, and generally increasing the safety of shipping. References [1] L. Gracia, H. Wang, W. Mao, N. Osawa, A study on the difference in wave statistics based on storm model for the ship structural assessment, in proc. 39th International Conference on Ocean, Offshore and Arctic Engineering (ASME 2020), ISOPE I 17 583, 2020. doi: 10.1115/OMAE2020- 19048. [2] K. Takeda, M. Akagi, K. Ishibashi, Introduction of Guidelines on Preventive Measures against Parametric Rolling, ClassNK Technical Journal 7 (2023) 13 21. [3] K. Ranjan, M. Ranjan, Controlling period-doubling route to chaos phenomena of roll oscillations of a biased ship in regular sea waves. Nonlinear Dynamics 111 4 (2023) 13889 13918. doi: 10.21203/rs.3.rs-2506716/v1. [4] S. Kianejad, H. Enshaei, J. Duffy, N. Ansarifard, D. Ranmuthugala, Ship Roll Damping Coefficient Prediction Using CFD, Journal of Ship Research, 63 (2019) 108 122. doi: 10.5957/JOSR.09180061. [5] V. Marasanov, D. Stepanchikov, A. Sharko, O. Sharko, Operator of the Dynamic Process of the Appearance of Acoustic Emission Signals during Deforming the Structure of Materials, in: 2020 IEEE 40th International Conference on Electronics and Nanotechnology, ELNANO 2020 Proceedings, 2020, pp. 646 650. doi: 10.1109/ELNANO50318.2020.9088893. [6] V. Marasanov, D. Stepanchikov, A. Sharko, O. Sharko, Technology for determining the residual life of metal structures under conditions of combined loading according to acoustic emission measurements, Communications in Computer and In-formation Science 1158 (2020) 202 217. doi:10.1007/978-3-030-61656-4_13. [7] M. Sharko, O. Liubchuk, G. Krapivina, N. Petrushenko, O. Gonchar, K. Vorobyova, N. Vasylenko, Using Markov Chains, Lecture Notes on Data Engineering and Communications Technologies 149 (2023) 197 213. doi:10.1007/978-3-031-16203-9?page=1#toc. [8] M. V. Gorp, J. Goyens, M. Alfaro, S. Wassenbergh, Keels of boxfish carapaces strongly improve stabilization against roll, Journal of The Royal Society Interface 19 189 (2022). doi: 19. 10.1098/rsif.2021.0942. [9] J. Lian, H. Wang, Y. Guo, W. Pengwen, Effects of damping plate on the motion response of transport ships under waves, Applied Ocean Research 134 (2023). doi: 10.1016/j.apor.2023.103507. [10] H. Wu, R. Rao, H. Guo, D. Zhang, X. Li, Li Zhao, Z. Li, Y. Peng, Research on performance of solid-liquid triboelectric nanogenerators based on anti rolling tank, Applied Energy 353 122153 (2024). doi: 10.1016/j.apenergy.2023.122153. [11] G. Bulian, J. Cercos Pita, Co simulation of ship motions and sloshing in tanks, Ocean Engineering, 152 (2018) 353 376. doi: 10.1016/j.oceaneng.2018.01.028. [12] V. Khramushin, Stormy Seakeeping and Navigation Safety Researches for Hull Form Design, in: 2nd Asia Pacific Workshop on Marine Hydrodynamics APHydro 2004, at: Sangnam International House, Pusan National University, Busan, Korea, 2004. [13] A. Wankhede, What to do when a ship encounters rough weather? Marine sight, 2021. https://www.marineinsight.com/marine-safety/what to do when ship encounters rough weather/. [14] Capt. Takuzo Okada, Marine Weather Ship Handling in Rough Seas, Japan P&I Club. P&I Loss Prevention Bulletin 45 (2019) 108. [15] O. D. Pipchenko, Development of the theory and practice of risk management in solving complex navigation problems, The Doctoral Thesis, National University Odesa Maritime Academy, 2021. [16] I. Ghamari, M. Greco, M. Faltinsen, C. Lugni, Numerical and experimental study on the parametric roll resonance for a fishing vessel with and without forward speed, Applied Ocean Research, 101 (2020). doi: 10.1016/j.apor.2020.102272. [17] A. Manoj, J. Jose, S. Kumar, Measures to minimize parametric rolling in container ships, in: Conference: 35th National Convention of Marine Engineers at: Mylapore, Chennai, Oktober 2023. [18] S. Zinchenko, V. Kobets, O. Tovstokoryi, K. Kyrychenko, P. Nosov, I. Popovych, Control of the pivot point position of a conventional single screw vessel, CEUR WS.org (ICST 2023), Vol. 3513, 2023, pp. 130 140. https://ceur ws.org/Vol 3513/paper11.pdf [19] S. Zinchenko, K. Kyrychenko, O. Grosheva, P. Nosov, I. Popovych, P. Mamenko, Automatic reset of kinetic energy in case of inevitable collision of ships, in: IEEE Xplore, 13th International 2023, pp. 496 500. doi: 10.1109/ACIT58437.2023.10275545. [20] S. Zinchenko, V. Kobets, O. Tovstokoryi, P. Nosov, I. Popovych, Intelligent System Control of the Vessel Executive Devices Redundant Structure, CEUR Workshop Proceedings, vol. 3403, 2023, pp. 582 594. https://ceur ws.org/Vol 3403/. [21] S. Zinchenko, O. Tovstokoryi, V. Mateichuk, P. Nosov, I. Popovych, V. Perederyi, Automatic 3668, 2024, pp. 235 246. https://ceur ws.org/Vol 3668/paper16.pdf. [22] P. Nosov, O. Koretsky, S. Zinchenko, Yu. Prokopchuk, I. Gritsuk, I. Sokol, K. Kyrychenko, Devising an approach to safety management of vessel control through the identification of navigator's state, Eastern European Journal of Enterprise Technologies, 4 3 (2023) 19 32. doi: 10.15587/1720-4061.2023.286156. .