Importance of several nonlinear factors on broaching prediction

In order to realize a more quantitative prediction of broaching and capsizing in following and quartering seas, existing mathematical modelling techniques should be upgraded. Therefore, it is necessary to systematically examine all factors relevant to capsizing in following and quartering seas. To this end, we first attempted to examine the prediction accuracy of wave-induced forces by comparing calculations with captive model experiments. As a result, we found that a wave-induced surge force has a certain nonlinearitiy with respect to wave steepness. The nonlinearity of sway–roll coupling with respect to sway velocity was also found, and our numerical result with these nonlinearities improves the prediction accuracy of the critical ship speed of capsizing in following and quartering seas. The importance of the wave effect on propeller thrust was also examined. We found that this effect is not negligibly small and could improve capsizing predictions in heavy following and quartering seas. Finally, we attempted to investigate the importance of nonlinear heel-induced hydrodynamic forces on numerical predictions of capsizing due to broaching. Here, we propose a new procedure for captive model experiments to obtain hydrodynamic forces with various heel angles up to 90°, and data on heel-induced hydrodynamic forces with respect to heel angle in calm water are provided. We then compare the numerical simulations with the nonlinear heel-induced hydrodynamic forces and without them. These time series comparisons show that the effect of nonlinear heel-induced hydrodynamic forces in calm water is not negligibly small for the case of ship capsizing due to broaching.     

Parametrical studies of a new numerical model for controlled ship motions in extreme astern seas

Parametrical studies based on numerical simulations were carried out for very steep regular waves to assess possible improvements in the state-of-the art numerical modelling of the control and capsizing behaviour of ships in following and quartering seas. A nonlinear 6-DOF numerical model has been developed with the inclusion of frequency-dependent terms, the so called memory effects, and a flexible axis system that allows straightforward combination of seakeeping and manoeuvring models while accounting for extreme motions. The previously undertaken validation analyses using extensive model test data provided qualitatively good agreement, whereas the comparison with numerical models without coupling of the vertical motions and frequency-dependent hydrodynamic terms embodied in radiation forces identified improvements in the accuracy. However, to broaden the assessment of the numerical model, further parametrical numerical analyses were carried out using two ships, which had previously been tested in the validation analyses, for various operational and environmental conditions. These parameters were changed in accordance with the recommendations from international organisations and experience from model tests to realise and avoid dangerous conditions that often result in capsizing, such as broaching associated with surf riding and low-cycle resonance. As a result of the parametric analysis, we discuss the sensitivity of the improvements in the numerical model for various critical operational and design parameters and its possible use to provide a link between the ship's behavior and these parameters.     

Nonlinear dynamics of ship capsizing due to broaching in following and quartering seas

To provide a theoretical methodology to predict the critical condition for capsizing due to broaching, a nonlinear dynamical system approach was applied to the surge–sway–yaw–roll motion of a ship running with an autopilot in following and quartering seas. Fixed points of a mathematical model for the ship motion and unstable manifolds of the fixed point near the wave crest were systematically investigated. As a result, the existence of heteroclinic bifurcation was identified. With numerical experiments, it was confirmed that this heteroclinic bifurcation reasonably well represents the critical condition for capsizing due to broaching. Thus the nonlinear dynamical approach can be substituted for tedious numerical experiments. Received for publication on Nov. 20, 1998; accepted on March 16, 1999     

Stability assessment for intact ships in the light of model experiments

A systematic method for assessing intact ship stability with a free-running model in a seakeeping and maneuvering basin is proposed in this paper. Model experiments were carried out in extremely steep regular waves for a model drifting, running in head seas, and quartering seas. This method was applied to two purse seiners, and efficiently identified thresholds in metacentric heights for capsizing of these ships. These capsizing thresholds are compared with requirements of the IMO Code on Intact Stability. This series of model experiments also confirms that capsizing at the threshold occurs only in quartering seas, and shows that capsizing is caused by broaching, loss of stability on a wave crest, or bow diving. Received for publication on Jan. 20, 1999; accepted on July 6, 1999     

Numerical prediction of the surf-riding threshold of a ship in stern quartering waves in the light of bifurcation theory

In order to develop design and operational criteria to be used at the International Maritime Organization (IMO), critical conditions for broaching are explored in the light of bifurcation analysis. Since surf-riding, which is a prerequisite to broaching, can be regarded as a heteroclinic bifurcation, one of global bifurcations, of a surge-sway-yaw-roll model in quartering waves, the relevant bifurcation condition is formulated with a rigorous mathematical background. Then an efficient numerical solution procedure suitable for tracing the surf-riding threshold hypersurface is presented with successful examples. This deals with all state and control variables in parallel, and excludes backward time integration and an orthogonal condition in the iteration process. The bifurcation conditions identified were compared with the results from a direct numerical simulation in the time domain. As a result, it was confirmed that the heteroclinic bifurcation provides a boundary between motions periodically overtaken by waves and nonperiodic motions such as surf-riding and broaching.     

骑浪/横甩薄弱性衡准和模型试验研究

目前国际海事组织(IMO)正在制定的第二代完整稳性衡准,其中就包括骑浪/横甩薄弱性衡准。文章首先介绍了最新骑浪/横甩薄弱性衡准方法,应用自编的衡准软件进行样船计算,分析了当前衡准的适用性。其次,开展了内倾船型在随浪和尾斜浪中的骑浪/横甩试验,试验中获得了四种与骑浪/横甩相关的运动特性:周期运动、稳定骑浪、横甩和横甩导致的倾覆,而且在某波浪条件下观察到船舶连续发生三次横甩的现象。最后,将内倾船型的骑浪/横甩薄弱性衡准计算结果与试验结果进行对比,验证了衡准方法对于内倾船型的适用性。     

Nonlinear analysis of parametric rolling in longitudinal and quartering seas with realistic modeling of roll-restoring moment

Parametric rolling of a containership in longitudinal and quartering seas is examined by applying nonlinear dynamics to a 1DOF mathematical model with realistic modeling of the wave effect on roll-restoring moment. In our previous work, we confirmed that a mathematical model with a roll-restoring moment in waves calculated with the Froude–Krylov assumption could considerably overestimate the danger of capsizing associated with parametric rolling. Therefore, in the present work, all numerical calculations based on nonlinear analysis were carried out with the direct aid of a measured roll-restoring moment in waves. For this purpose, captive model experiments were conducted for various sets of wavelengths in longitudinal seas. This experiment demonstrates that the Froude–Krylov prediction could not explain the wavelength effect on restoring moment as the wave-steepness effect. Using the numerical model with the aid of this measured roll-restoring moment, the Poincaré mapping technique was applied to identify bifurcation structures of roll motions not only in longitudinal seas, but also in quartering seas. As a result, it was confirmed that capsizing associated with parametric rolling is more likely to occur in following seas than in quartering seas. However, period-doubling and chaos appeared in quartering seas. Finally, an averaging method assuming a period-2 orbit was applied to the same model with the same conditions as the Poincaré map. Reasonably good agreement was obtained between the numerical results with a Poincaré map and those with the averaging method in longitudinal seas, but the averaging method has limited capability in quartering seas.     

URANS analysis of a broaching event in irregular quartering seas

Ship motions in a high sea state can have adverse effects on controllability, cause loss of stability, and ultimately compromise the survivability of the ship. In a broaching event, the ship losses control, naturally turning broadside to the waves, causing a dangerous situation and possibly capsizing. Classical approaches to study broaching rely on costly experimental programs and/or time-domain potential or system-based simulation codes. In this paper the ability of Reynolds averaged Navier–Stokes (RANS) to simulate a broaching event in irregular waves is demonstrated, and the extensive information available is used to analyze the broaching process. The demonstration nature of this paper is stressed, as opposed to a validated study. Unsteady RANS (URANS) provides a model based on first principles to capture phenomena such as coupling between sway, yaw, and roll, roll damping, effects of complex waves on righting arm, rudders partially out of the water, etc. The computational fluid dynamics (CFD) method uses a single-phase level-set approach to model the free surface, and dynamic overset grids to resolve large-amplitude motions. Before evaluating irregular seas two regular wave cases are demonstrated, one causing broaching and one causing stable surf riding. A sea state 8 is imposed following an irregular Bretschneider spectrum, and an autopilot was implemented to control heading and speed with two different gains for the heading controller. It is concluded that the autopilot causes the ship to be in an adverse dynamic condition at the beginning of the broaching process, and thus is partially responsible for the occurrence of the broaching event.     

Broaching prediction in the light of an enhanced mathematical model,with higher-order terms taken into account

 We have attempted to develop a more consistent mathematical model for capsizing associated with surf-riding in following and quartering waves by taking most of the second-order terms of the waves into account. The wave effects on the hull maneuvring coefficients were estimated, together with the hydrodynamic lift due to wave fluid velocity, and the change in added mass due to relative wave elevations. The wave effects on the hydrodynamic derivatives with respect to rudder angles were estimated by using the Mathematical Modelling Group (MMG) model. Then captive ship model experiments were conducted, and these showed reasonably good agreements between the experiments and the calculations for the wave effects on the hull and the rudder maneuvring forces. It was also found that the wave effects on restoring moments are much smaller than the Froude–Krylov prediction, and the minimum restoring arm appears on a wave downslope but not on a wave crest amidship. Thus, an experimental formula of the lift force due to the heel angle of the ship is provided for numerical modelling. Numerical simulations were then carried out with these second-order terms of waves, and the results were compared with the results of free-running model experiments. An improved prediction accuracy for ship motions in following and quartering seas was demonstrated. Although the boundaries of the ship motion modes were also obtained with both the original model and the present one, the second-order terms for waves are not so crucial for predicting the capsizing boundaries themselves. Received: June 20, 2002 / Accepted: October 10, 2002 Acknowledgments. This research was supported by a Grant-in-Aid for Scientific Research of the Ministry of Education, Culture, Sports, Science and Technology of Japan (No. 13555270). The authors thank Prof. N. Rakhmanin of the Krylov Ship Research Institute for providing the Russian literature, as well as Mr. H. Murata of NHK (Japan Broadcasting Corporation) for translating it into Japanese. Address correspondence to: N. Umeda (e-mail: umeda@naoe.eng.osaka-u.ac.jp)     

Broaching prediction of a wave-piercing tumblehome vessel with twin screws and twin rudders

The new intact stability criteria which are under development at the International Maritime Organization (IMO) are required to cover a broaching phenomenon, well known as a great threat to high-speed vessels which can lead to capsizing. Some reports exist which demonstrate that their numerical models can predict a highly nonlinear phenomenon of broaching. However, additional validation studies are needed for unconventional vessels, in addition to conventional ones, to develop direct stability assessment methods for the new intact stability criteria. In this research, we selected as the subject ship a wave-piercing tumblehome vessel with twin screws and twin rudders, a design expected to be one of a new generation of high-speed monohull ships. Firstly, a series of captive model tests were conducted to measure the resistance, the manoeuvring forces, the wave-exciting forces, the heel-induced hydrodynamic forces, and the roll restoring variation for the unconventional tumblehome vessel. Secondly, the existing mathematical model which had been developed for broaching prediction of conventional vessels with a single propeller and a single rudder was extended to unconventional vessels with twin propellers and twin rudders. Finally, comparisons between numerical simulations and the existing free running model experiments were conducted. As a result, it was demonstrated that fair quantitative prediction of broaching is realised when the rudder force variation, the roll restoring variation and the heel-induced hydrodynamic force for large heel angles are taken into account.     

Melnikov integral formula for beam sea roll motion utilizing a non-Hamiltonian exact heteroclinic orbit: analytic extension and numerical validation

In the research field of nonlinear dynamical system theory, it is well known that a homoclinic/heteroclinic point leads to unpredictable motions, such as chaos. Melnikov’s method enables us to judge whether the system has a homoclinic/heteroclinic orbit. Therefore, in order to assess a vessel's safety with respect to capsizing, Melnikov’s method has been applied for investigations of the chaos that appears in beam sea rolling. This is because chaos is closely related to capsizing incidents. In a previous paper (Maki et al. in J Mar Sci Technol 15:102–106, 2010), a formula to predict the capsizing boundary by applying Melnikov’s method to analytically obtain the non-Hamiltonian heteroclinic orbit was proposed. However, in that paper, only limited numerical investigation was carried out. Therefore, further comparative research between the analytical and numerical results is conducted, with the result being that the formula is validated.     

船舶在波浪中航行时的安全操纵

船舶在波浪中航行,纯稳性丧失、参数激振和横甩是造成船舶倾覆的主要原因。针对波长与船长、波高与波长、波与船的波舷角三者对船舶稳性的影响进行讨论,揭示了船舶在波浪中航行时的稳性变化规律,提出了应合理地选择船舶的航行姿态、谨慎用舵等操船建议,保证船舶的航行安全。     

船舶横浪倾覆试验及其数值模拟

本文介绍船模倾覆试验的结果及其数值模拟。试验表明，在横浪情况下船舶因装载造成的横倾角，会大大增加倾覆的危险性。对这一结论进行的数值模拟与实测结果符合良好。     

随机横浪中甲板上浪船舶的倾覆研究

应用相空间转移率,定量研究了随机海浪中甲板上浪船舶的倾覆,给出了甲板上浪对船舶抗倾覆能力的影响.综合考虑非线性阻尼、非线性恢复力矩和随机横浪激励,建立了无甲板上浪和有甲板上浪时船舶随机非线性横摇运动的一般方程.以一艘倾覆的拖网船为例,分别求解了无甲板上浪和有甲板上浪时,不同海况激励下船舶横摇的相对相空间转移率.以相空间转移率作为船舶稳性损失的度量,定量比较了两种情况下船舶的抗倾覆能力.研究表明,甲板上浪后,船舶在较低的海况下会产生较大的相对相空间转移率,甲板上浪严重降低船舶的抗倾覆能力,从理论上进一步揭示了甲板上浪船舶的倾覆机理.     

Dynamic stability in following seas: predictive and experimental approaches

This article presents work based on the development of a performance-based stability assessment method. It describes a numerical method used to determine the survival limit for a dynamic intact stability assessment procedure. The numerical method utilises a time-domain vessel motion program to assess the limit for a range of vertical centres of gravity (KG). The appropriateness of the numerical predictions was examined through comparison with model experiment results. Free-running model tests were conducted in regular following waves at discrete KGs. A comparison between the survival limits determined through the numerical and experimental methods is presented. The current International Maritime Organisation (IMO) stability criteria are also evaluated against the numerical and experimental dynamic performance-based stability assessment methods.     

Broaching probability for a ship in irregular stern-quartering waves: theoretical prediction and experimental validation

To avoid stability failure due to the broaching associated with surf riding, the International Maritime Organization (IMO) has begun to develop multilayered intact stability criteria. A theoretical model using deterministic ship dynamics and stochastic wave theory is a candidate for the highest layer of this scheme. To complete the project, experimental validation of the theoretical method for estimating broaching probability in irregular waves is indispensable. We therefore conducted free-running model experiments using a typical twin-propeller and twin-rudder ship in irregular waves. A simulation model of coupled surge–sway–yaw–roll motion was simultaneously refined. The broaching probability calculated by the theoretical method was within the 95 % confidence interval of that obtained from the experimental data. This could be an example of experimental validation of the theoretical method for estimating the broaching probability when a ship meets a wave.     

船舶低速横甩机理研究

本文利用非线性振动理论研究了参数激励频率和波浪的强迫激励频率相同且等于艏摇固有频率2倍时的船舶1/2亚谐共振运动,揭示了船舶在低速航行、遭遇频率相对较高时,积累艏摇运动导致横甩发生的机理.     

船舶大舵角转向时艏摇与横摇的耦合仿真研究

摘要：为进一步研究船舶在大幅度转向时艏摇和横摇的耦合机制,在非线性船舶运动数学模型的基础上,进行了不同情形下的操纵仿真试验。试验结果显示在大舵角转向时,船舶的艏摇和横摇运动存在较强的耦合作用,横摇幅度和艏摇幅度存在正相关性。指出大幅度的横摇使艏摇出现相位滞后和偏离基准航向的现象;在横摇过大的情况下,大幅度动舵和加速操纵将导致稳性迅速消失和航向打横。     

Prediction of capsizing probability for a ship with trapped water on deck

The authors have already examined a method for evaluating the capsizing probability of a ship in the dead ship condition based on a piecewise linear approximation of the restoring arm. Here, this method is extended to ships with trapped water on deck. This is because the stability of ships having a relatively high bulwark, such as fishing vessels, could substantially deteriorate due to trapped water on deck. First, the mean amount of water trapped on deck was estimated as a function of the significant wave height and the mean wave period using a model experiment in irregular beam seas. Second, the restoring arm curve with trapped water on deck was calculated hydrostatically and then approximated with a piecewise linear curve. Third, the roll angle was estimated using a nonlinear and uncoupled equation of absolute roll angle under stochastic wave and wind exciting moments. The short-term and long-term capsizing probabilities were calculated for a fishing vessel operating off Kyushu. Numerical results quantitatively demonstrated that the effect on capsizing probability of trapped water on deck cannot be ignored when accurately evaluating the stability of fishing vessels.     

On the time dependence of survivability of ROPAX ships

The time dependence of survivability of ROPAX vessels, when sustaining side collision damage in waves, is investigated herein by use of numerical simulations of ship motion and flooding. Conducted research confirms that ROPAX ships characteristically capsize fast, when sustaining damage leading to capsizing. A probabilistic analysis of the survive time after collision damage reveals that even for the most generic damage conditions assumed, the survival time in the case of capsizing remains short, which is characteristic of this type of ship design, disposing the typically large undivided deck to be subject to flooding in higher waves. In a case study, the unconditional survivability in waves, corresponding to survival s-factor of SOLAS regulations, is alternatively assessed with numerical simulations. The estimated survivability proves to be time independent in terms of practical implications. Observed deviations between current SOLAS formulation and simulations, suggest the employment of comprehensive simulation methods when more reliable estimations and assessments of survivability are required.