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.     

Qualitative aspects of nonlinear ship motions in following and quartering seas with high forward velocity

By utilizing a four-degrees-of-freedom numerical model with dense grids of control parameters and the sudden-change concept, the qualitative aspects of the nonlinear motions of a fishing vessel complying with the International Maritime Organization's intact stability criteria in following and quartering seas were intensively explored. As a result, capsizing due to broaching, capsizing without broaching, broaching without capsizing, stable surf-riding, and steady periodic motion were identified. The natures of the boundaries of these motions in the control parameter plane were investigated, and the effects of the initial conditions and the nonlinearity of calm-water maneuvering forces are also discussed. Furthermore, comparisons with a model experiment showed that the numerical model used here qualitatively explains capsizing phenomena, but quantitatively overestimates the danger of capsizing. Received: June 11, 2001 / Accepted: October 9, 2001     

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.     

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     

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)     

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.     

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.     

Analytical formulae for predicting the surf-riding threshold for a ship in following seas

Making use of Melnikov’s method, a generalized formula for predicting the surf-riding threshold is developed as an extension to the applications of Kan and Spyrou. A new analytical formula for calculating the surf-riding threshold of a ship in following seas is also proposed in light of nonlinear dynamical system theory. By applying a continuous piecewise linear approximation to the wave-induced surge force, a heteroclinic bifurcation point is obtained analytically with an uncoupled surge equation. Results calculated using these formulae are presented, and they show good agreement with those obtained utilizing numerical bifurcation analysis. Further, it was confirmed that the surf-riding threshold obtained using the proposed formulae agrees reasonably well with that obtained experimentally for an unconventional vessel.     

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.     

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

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

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.     

风帆助航船非线性横摇倾覆概率研究

风帆助航技术为船舶利用可再生能源提供了途径,但随着横摇的幅度加大,导致船舶倾覆的概率也增大。文章采用数值分析的方法对风帆助航船舶的非线性横摇方程进行分析求解,对有风帆助航与无风帆助航船在风浪作用下产生倾覆的概率作了比较,为采用风帆助航技术的远洋船舶在符合初稳性和气象衡准稳性条件下的横摇倾覆防范提供研究思路。     

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.     

Large-amplitude motion responses of a Ro–Ro ship to regular oblique waves in intact and damaged conditions

This article presents a nonlinear time-domain simulation method for the prediction of large-amplitude motions of a Ro–Ro ship in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried out to investigate the dynamic motion responses of Ro–Ro ship Dextra to various wave amplitudes at three different wave headings. The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements show good agreement except in the roll-resonant region. Nonlinear effects are significant in horizontal modes of motion, and resonant roll motion, and there is strong coupling between all modes of motion in the roll-resonant region for large-amplitude responses. On the other hand, the time-domain simulation technique suffers from numerical drift in horizontal modes of motion as wave amplitude increases. This is due to nonlinear equations of motion and the lack of a restoring force and moment in horizontal motion. Received: April 30, 2002 / Accepted: August 9, 2002 Acknowledgments. II Programme of the European Community Commission under contract No. BRPR-CT97-0513. Address correspondence to: H.S. Chan (hoi-sang.chan@ncl.ac.uk)     

船舶参数激励非线性横摇运动方程

讨论了船舶在海浪上运动时参数激励产生的机理，建立了船舶参数激励和波浪激励非线性横摇运动方程。本文工作为研究倾覆机理奠定了基础。     

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

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

随机横浪中船舶倾覆过程模拟与安全域计算

研讨了随机海浪中船舶非线性随机横摇运动倾覆过程和安全域问题.考虑非线性阻尼、非线性复原力矩和随机波浪,建立了随机横浪中船舶运动的随机非线性微分方程.考虑两参数海浪谱和不同的有义波高,计算随机波浪力矩,采用简谐加速度数值计算方法求解随机微分方程,进行了随机倾覆过程的数值模拟.考虑船舶运动的瞬时状态,发展了胞映射方法计算随机非线性横摇运动安全域的方法,并编制了计算机程序.模拟了一条渔船不同有义波高作用下横摇倾覆过程,考虑不同的初始条件及随机波浪参数,在初始值平面上采用胞映射法构造了渔船航行的安全域.波高较小时,安全域是一个连续的区域;随着有义波高增大,安全域面积减小并且不再连续,倾覆域和安全域相互包围.计算结果表明,船舶运动的瞬时状态对安全域具有重要影响,而且随机海浪中安全域具有随机特性.     

Study of damaged ship motions taking into account floodwater dynamics

This paper presents an application of shallow water theory to describe the motion of floodwater inside a rolling ship in damage condition. The time domain theoretical approach to the coupled problems of ship and water inside compartment motions is briefly described, including the method used to solve for the water motion characteristics and forces exerted on the ship. This approach is applied to the study of the behaviour of a passenger Ro–Ro ship in regular beam seas and numerical results are given for the intact and damaged conditions. Comparison is made with experimental results. For the damaged condition, the characteristics of the floodwater motion are studied in the time domain for a number of different wave frequencies. The shape of the free-surface and phase of water motion in relation to the ship roll motion are shown for several wave frequencies. The dynamic floodwater roll moment is also shown and compared with the static roll moment (flat horizontal free surface), allowing the conclusion that the dynamic roll moment is much larger than the static roll moment, for high wave frequencies, and is in phase opposition in relation to the roll motion.     

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

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

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.