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)     

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     

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     

Predicting parametric rolling in irregular head seas with added resistance taken into account

Parametric rolling in irregular head seas is investigated using numerical simulation in the time domain in this paper. The effect of the change in speed due to added resistance on parametric rolling is taken into account. Three methods of calculating the added resistance in irregular head seas are applied to this problem, and then Pinkster’s method, which considers the time-varying added resistance in irregular head seas, is extended for further investigations in this paper. Furthermore, the Froude–Krylov component of restoring variation is estimated with coupling from heave and pitch motions, and the hydrodynamic effects are extrapolated nonlinearly with regards to roll angle. Finally, the “practical non-ergodicity” of parametric roll is investigated with a model experiment at zero forward speed and using the simulation, and the effect of added resistance is investigated with the simulation.     

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     

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.     

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

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

Stochastic averaging for estimating a ship roll in random longitudinal or oblique waves

In this work, the C11 container ship is taken as an example to analyze its rolling performances in random longitudinal or oblique waves. Firstly, a dynamic model of C11 roll in random waves is improved, and it is verified by the model test and numerical simulation. Mathematically, this dynamic model is a one-dimensional stochastic differential equation with random parametric (and external) excitation. Secondly, an enhanced stochastic averaging method is proposed to solve this stochastic differential equation. The validity of the solutions was verified by Monte Carlo simulation. At last, the probabilistic characteristics of the extreme rolling response were investigated based on the calculated results using enhanced stochastic averaging method. According to the analysis, some advices for ship's manoeuvring can be put forward when ships are navigating in random waves.     

On the Parametric Resonance of Tumblehome Hullforms in a Longitudinal Seaway

This work presents a numerical study of the influence of topside shape on parametric roll in longitudinal seas. The ONR 3 Topsides hulls are used to compare the single degree-of-freedom rolling response of a destroyer-sized vessel with flared (model 5613), tumblehome (model 5613-1), and wallsided (model 5613-2) shapes above the waterline. Through this relative comparison, the key aspects of parametric resonance of a tumblehome hullform are highlighted, namely the reduction in head seas forward speed required to encounter parametric roll. A variant on the Mathieu equation is used to simulate roll motions.     

An investigation of different methods for the prevention of parametric rolling

The parametric rolling of modern containerships is emerging as a serious problem, to the extent that its effects warrant a study into its prevention. In light of this, two methods for reduction of parametric rolling are proposed and examined by physical model experiments. The first is a sponson attached to the side of a ship, the purpose being to decrease the rate of change of the rollrestoring moment. The second is an antirolling tank to increase roll damping. By conducting free-running model experiments for a 6600-TEU post-Panamax container ship with sponsons under typical parametric rolling conditions, it was found that the sponsons could decrease the magnitude of parametric rolling. The antirolling tank could prevent parametric rolling completely in certain conditions, even in severe head seas. Using the damping coefficients from experimentally derived data of a model ship with an antiroll tank, a numerical simulation was established. The numerical model was then compared with the free-running model experiments. The results indicated that the numerical model could qualitatively verify the experimental results. Finally, an attempt to optimise the size of an antirolling tank for preventing parametric rolling for the subject post-Panamax container ship in the North Pacific Ocean is presented.     

Parametric rolling prediction in irregular seas using combination of deterministic ship dynamics and probabilistic wave theory

In recent years there have been reports of serious accidents of parametric rolling for modern container ships and car carriers. For avoiding such accidents, a prediction method of parametric rolling in irregular seas is required. Since parametric rolling is practically non-ergodic, repetitions of numerical simulations or experiments could be not feasible to ascertain the behaviour. Therefore, in this paper, a method combining a stochastic approach with a deterministic approach in order to estimate the probabilistic index without such simple repetitions is developed. The ship's response in regular seas is estimated by solving an averaged system of the original 1-DoF roll model, and random waves necessary for occurrence of parametric rolling is achieved by using Longuet-Higgins’s or Kimura’s wave group theory. As a result, a fast and robust computation method of the probabilistic index is established. Finally, it is concluded that the proposed method is considered to be one of the useful tools for discussing the new IMO Intact Stability Code.     

纵向规则波中参数横摇的数值模拟

基于ITTC推荐的参数横摇运动方程,采用估算方法确定方程中的阻尼系数和复原力矩系数,数值模拟出船舶在规则波中的参数横摇情况,将上述方法具体应用于一艘滚装船船模的参数横摇预报,并与试验数据进行对比,验证了所采用方法的有效性。     

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

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

大型汽车滚装船参数横摇研究

对于航行在纵向波浪中的船舶,参数横摇是横稳性中典型的不利情形。文章对参数横摇作了介绍,并对某大型汽车滚装船进行了实例分析。首先,对船舶在纵向波浪中的稳性变化进行了计算;然后,根据衡准进行验证;最后,采用数值仿真方法,研究了稳性变化、横摇阻尼、航速、航向等参数对参数横摇的影响。     

Nonlinear effects from wave-induced maximum vertical bending moment on a flexible ultra-large containership model in severe head and oblique seas

Vertical bending moment (VBM) is of crucial importance in ensuring the survival of vessels in rough seas. With regard to conventional vessels, wave-induced maximum VBM is normally considered to be experienced in head seas. It is conservative to determine the extreme VBM based on either numerical simulations or model tests in long-crested head seas. Extensive model tests have been conducted in head seas with focus on the nonlinear vertical responses in severe seas, and the measured results were compared with numerical calculations for validation. Unexpected phenomena, however, were observed during the model tests of an ultra-large containership. The maximum sagging and hogging VBMs were encountered in oblique seas. Furthermore, the significant wave height used in oblique seas was even smaller than that used in head seas. The nonlinear vertical load effects in oblique seas require further investigations for this particular vessel. Limited experimental results in oblique seas have been reported, in which the lateral responses were always more concerned than the vertical responses. Up to now, rare systematic comparisons of the nonlinear vertical responses between head and oblique seas have been published, especially when the hydroelastic effects are also accounted for. A 13000-TEU ultra-large containership model, which was designed by Hyundai Heavy Industries (HHI), has been tested in the towing tank and the ocean basin at the Marintek center in Trondheim. The experimental results in regular waves are first compared between head and oblique seas. The statistical characteristics of the VBM amidships under nineteen irregular wave conditions are then investigated. Next, the extreme hogging and sagging VBMs are compared under different wave conditions with focus on the extreme hogging VBMs. At the end of the paper, the uncertainties in the experiments are discussed.     

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

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

Theoretical study on the effect of parametric rolling on added resistance in regular head seas

A formula based on Maruo’s theory is presented for the added resistance in regular head seas when parametric rolling occurs. Firstly, the velocity potential of the radiation waves due to parametric rolling, sway and yaw motions should be superposed on the conventional velocity potential. Secondly, the added resistance is averaged within the time duration that is double the encounter period. Thirdly, the stationary-phase method which is used in Maruo’s method is also used to obtain the formula of added resistance in waves with parametric rolling taken into account as well as viscous roll damping. Fourthly, source distribution based on Maruo and Ishii's and Maruo and Iwase’s works is used to calculate the added resistance. Finally, the effect of parametric rolling on added resistance in regular head seas mainly results from viscous roll damping and it becomes larger as the amplitude of parametric rolling becomes larger.     

规则迎浪中船舶参数横摇的三维时域预报方法研究

参数横摇是船舶因复原特性改变而引起的典型非线性现象。文章采用三维时域方法预报规则迎浪中船舶的参数横摇运动。该方法引入匹配面将流域分为内域和外域,内域中采用Rankine源来满足物面条件和线性自由面条件,而外域中应用时域格林函数来满足线性自由面条件和远场辐射条件。数值方法中,Froude-Krylov力和恢复力是通过对船舶瞬时湿表面积积分获得,同时考虑了横摇、垂荡和纵摇三自由度之间的耦合作用,以及非线性横摇阻尼的影响。数值结果与试验结果吻合很好,说明该方法可以有效地预报参数横摇。     

船舶参数激励非线性随机横摇运动分析

研究了随机纵浪中船舶的参数激励横摇运动。将船舶初稳性高作为时变的随机参数激励,应用随机平均法和FPK方程,分析了船舶的参数激励横摇运动微分方程。考虑航速、船舶瞬时湿表面积及随机波浪,建立了船舶横摇幅值概率密度、概率分布函数及其均值的计算方法。以某渔政船为例,计算了不同航速时船舶随机横摇幅值的概率密度、概率分布函数及其均值。结果表明,参数激励谱峰频率为横摇固有频率的2倍时,发生主参数共振,船舶横摇严重,船舶倾覆发生的概率增加。