船舶随浪运动稳性仿真计算

本文利用Ｌｉａｐｕｎｏｖ理论，研究了船舶在规则波浪运动的稳性；利用摄动理论，求解出船舶运动响应；并讨论了船舶横摇与垂荡运动频率、最大横摇角和波浪要素对稳性曲线ＧＺ的影响，以及流体动压力对稳性曲线的修正，从而给计算船舶在随浪中的稳性提供了一种方法。     

Estimation of survival probability for a ship in beam seas using the safe basin

The purpose of this paper is to analyze the nonlinear ship roll motion equation and the main parameters that induce ship capsizing in beam seas, estimate the survival probability of a ferry in random seas and to find out a risk assessment method for the ship’s intact stability. A single degree of freedom (1-DOF) dynamic system of ship rolling in beam seas is investigated and the nonlinear differential equation is solved in the time domain by the fourth order Runge-Kutta algorithm. The survival probability of a ferry in beam seas is investigated using the theory of “safe basin”. The survival probability is calculated by estimating erosion of “safe basin” during ship rolling motion by Monte Carlo simulations. From the results it can be concluded that the survival probability of a ship in beam sea condition can be predicted by combining Monte Carlo simulations and the theory of “safe basin”.     

A method for breach assessment onboard a damaged passenger ship

Reliable analysis of stability and safety level in a flooding emergency onboard a damaged passenger ship is extremely important for making correct decisions on evacuation and abandonment. Thus there is demand for an automated system that detects flooding and analyzes the severity of the situation. This procedure requires estimation of the actual breach in the hull, and calculation of possible progressive flooding to undamaged compartments. A new approach to the breach assessment, based on measurement data from the flood level sensors, is presented. The developed method is complemented by time-domain flooding simulations in order to separate progressive flooding from direct inflow through the breaches in the hull of the damaged ship. The developed approach is tested and demonstrated with a large passenger ship design for various damage scenarios. The results show that the size and location of the breach can be evaluated with reasonable accuracy on the basis of the level sensor data, provided that there are enough well-placed, working level sensors.     

Stability and capsizing analysis of nonlinear ship rolling in wind and stochastic beam seas

Considering the actual seaway condition, stability and capsizing of nonlinear ship rolling system in stochastic beam seas is of significant importance for voyage safety. Safe zone are defined in the phase space plan of the unperturbed Hamilton system to qualitatively distinguish ship motions as capsize and noncapsize. Capsize events are defined by solutions passing out of the safe zone. The probability of such an occurrence is studied by virtue of the random Melnikov function and the concept of phase space flux. In this paper, besides conventional wave excitation, the effect of wind load is also taken into account. The introduction of wind load will lead to asymmetry, in other words, it transforms the symmetric heteroclinic orbits into asymmetric homoclinic orbits. For asymmetric dynamical system, the orbital analytic solutions and its power spectrum are not readily available, and the technique of discrete time Fourier transformation (DTFT) is used. In the end, as verification of theoretical critical significant wave height, capsizing probability contour diagram is generated by means of numerical simulation. The contour diagram shows that these analytical methods provide reliable and predictive results about the likelihood of a vessel capsizing in a given seaway condition.     

基于AIS的船舶会遇局面紧迫度量化分析

近年来船舶大型化、高速化的发展以及数量的不断递增,导致船舶碰撞概率大大增加。为更好地保障船舶航行安全,提出船舶会遇局面紧迫度这一概念。该概念表征船舶在当前航行状态下由于一定范围内所有其他船舶的存在而处于的会遇局面的紧迫程度,其旨在帮助船舶驾驶员实时感知并掌握船舶在海上航行时所处会遇局面的状况,从而据此加以判断并及时采取措施规避碰撞风险。建立了船舶会遇局面紧迫度的计算模型,该模型以考虑到本船速度的动态圆形领域为空间约束条件,综合考虑其他船舶相对于本船的方位分布和距离这两个指标的影响,基于信息量量化紧迫度。利用2016年1月1日天津港AIS数据进行紧迫度的计算与分析,计算结果与实际相符,该计算模型能够准确反映船舶在航行过程中所处会遇局面的紧迫程度。利用双线性内插法建立天津港区域紧迫度二维彩色平面分布图,其表达效果优于传统热度图。     

Parametric roll due to hull instantaneous volumetric changes and speed variations

Parametric roll of a containership in head sea condition has been studied in the paper. A time domain routine for GZ righting arm calculation based on exact underwater hull geometry has been implemented into a two-degree-of-freedom procedure for roll response calculation. The speed variation due to e.g. added resistance has been accounted for in the model by the surge velocity. The ship roll motion due to a regular wave critical for parametric roll occurrence has been simulated, as well as the ship roll response in a severe stochastic sea. The present method has been compared with other existing methods for parametric roll prediction.     

Stochastic inversion for the roll gyradius second moment mass property in ships at full-scale and model-scale

In current Naval Architecture practice, employing static considerations is an important and necessary step in assessing ship stability and seakeeping properties (e.g. inclining experiments, load line regulations, range of stability calculations). However, damaged vessels and vessels operating in heavy weather or in conditions where topside icing is a concern may require an additional assessment of stability that considers dynamic effects. Within such contexts, the actual (i.e. current) second moment properties of the vessel mass become very important in the associated equations of motion for a given ship. One such critical second mass moment property is the roll gyradius, as it is closely related to the occurrence of capsizing. The present paper furnishes a means for reckoning the actual roll gyradius of a given ship operating within a seaway. The approach hinges on the formulation and solution of a stochastic inverse problem that leverages existing seakeeping software against the shipboard inertial measurement unit (IMU) telemetry. The method is demonstrated at full-scale and validated at model scale.     

A Time-Domain Coupled Model for Nonlinear Wave Forces on A Fixed Box-Shaped Ship in A Harbor

A 2-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear wave forces on a fixed box-shaped ship in a harbor.The domain is divided into an inner domain and an outer domain.The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations.The other area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions.Along the interface boundaries between the inner domain and the outer domain,the volume flux is assumed to be continuous and the wave pressures are equal.Relevant physical experiment is conducted to validate the present model.It is shown that the numerical results agree with the experimental data.Compared with the coupled model with the flow in the inner domain governed by the Laplace equation,the present coupled model is more efficient and its solution procedure is more simple,which is particularly useful for the study on the effect of the nonlinear wave forces on a fixed box-shaped ship in a large harbor.     

瘫船稳性第二层薄弱性衡准研究

基于IMO第二代完整稳性草案中瘫船稳性第二层薄弱性衡准的有关研究,分析草案提出的简化计算方法,采用一个单自由度船舶横摇运动微分方程,利用例船CEHIPAR2792数据,编译相应程序,求得不规则横风横浪作用下瘫船横摇运动相应特征,进而计算特定海况下的倾覆概率及特定海区的平均倾覆指数,并与草案值进行对比。分析不同拟合扶正力臂曲线方法和剩余稳性高对计算结果的影响,并对瘫船稳性第二层薄弱性衡准方法进行修改。     

Theoretical and experimental study on dynamic behavior of a damaged ship in waves

Most of the large scaled casualties are caused by loss of structural strength and stability due to the progressive flooding and the effect of waves and wind. To prevent foundering and structural failure, it is necessary to predict the motion of the damaged ship in waves.This paper describes the motion of damaged ship in waves resulting from a theoretical and experimental study. A time domain theoretical model, which can be applied to any type of ship or arrangement, for the prediction of damaged ship motion and accidental flooding has been developed considering the effects of flooding of compartments. To evaluate the accuracy of the model, model tests are carried out in ship motion basin for three different damaged conditions: engine room bottom damage, side shell damage and bow visor damage of Ro–Ro ship in regular and irregular waves with different wave heights and directions.     

Numerical Treatments and Applications of the 3D Transient Green Function

How to evaluate time-domain Green function and its gradients efficiently is the key problem to analyze ship hydrodynamics in time domain.Based on the Bessel function,an Ordinary Differential Equation(ODE)was derived for time-domain Green function and its gradients in this paper.A new efficient calculation method based on solving ODE is proposed.It has been demonstrated by the numerical calculation that this method can improve the precision of the time-domain Green function.Numerical research indicates that it is efficient to solve the hydrodynamic problems.     

Domain decomposition and matching for time-domain analysis of motions of ships advancing in head sea

A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare.     

Total (practical) stability of ships

In this paper the “total (practical) stability” concept is introduced to nonlinear forced rolling motion of a ship. This is achieved by employing “boundedness” and “Lyapunov's function” approach. In this respect two new theorems are proved and conditions and domain of Practical Stability are evaluated. The Paper also contains a critical review of the present status of international intact ship stability regulations. A qualitative discussion of oscillatory rolling motion and the capsizing phenomena is presented.     

Time-domain simulation of wave–structure interaction based on multi-transmitting formula coupled with damping zone method for radiation boundary condition

Based on the Rankine source, this paper proposed a time-domain method for analyzing the three-dimensional wave–structure interaction problem in irregular wave. A stable integral form of the free-surface boundary condition (IFBC) is employed to update the velocity potential on the free surface. A multi-transmitting formula, with an artificial wave speed, is used to eliminate the wave reflection for radiation condition on the artificial boundary. An effective multi-transmitting formula, coupled with damping zone method, is further used to analyze the irregular wave diffraction at the artificial boundary. We investigate hydrodynamic forces on floating structure and compare our solution to the frequency-domain solution. It is shown that long time simulation can be done with high stability and the numerical results agree well with the solution obtained under the frequency domain. The efficiency of the proposed multi-transmitting formula and the coupled methods for radiation boundary make them promising candidates in studying the irregular water wave problem in time domain.     

Investigation of Hydroelastic Ship Responses of An ULOC in Head Seas

Investigation of hydroelastic ship responses has been brought to the attention of the scientific and engineering world for several decades. There are two kinds of high-frequency vibrations in general ship responses to a large ocean-going ship in its shipping line, so-called springing and whipping, which are important for the determination of design wave load and fatigue damage as well. Because of the huge scale of an ultra large ore carrier (ULOC), it will suffer seldom slamming events in the ocean. The resonance vibration with high frequency is springing, which is caused by continuous wave excitation. In this paper, the wave-induced vibrations of the ULOC are addressed by experimental and numerical methods according to 2D and 3D hydroelasticity theories and an elastic model under full-load and ballast conditions. The influence of loading conditions on high-frequency vibration is studied both by numerical and experimental results. Wave-induced vibrations are higher under ballast condition including the wave frequency part, the multiple frequencies part, the 2-node and the 3-node vertical bending parts of the hydroelastic responses. The predicted results from the 2D method have less accuracy than the 3D method especially under ballast condition because of the slender-body assumption in the former method. The applicability of the 2D method and the further development of nonlinear effects to 3D method in the prediction of hydroelastic responses of the ULOC are discussed.     

Calculation of a capsizing rate of a ship in stochastic beam seas

This study introduces a method of calculating a capsizing rate of a ship. The phenomenon ‘capsizing’ is described as a jump of local equilibrium point from that near the upright position of a ship to what describes the upside-down attitude of the capsized ship; the rate of occurrence of such jumps was calculated. The potential function corresponding to the roll restoring moment have two potential wells located at the roll displacement angle 0 and 180°, respectively. A nonlinear Fokker–Planck equation for the joint probability density function of roll angle and velocity was solved. The excitation to the ship was assumed to be a combination of a regular harmonic wave and a white noise process.     

A Semi-Analytical Method for the PDFs of A Ship Rolling in Random Oblique Waves

The PDFs (probability density functions) and probability of a ship rolling under the random parametric and forced excitations were studied by a semi-analytical method. The rolling motion equation of the ship in random oblique waves was established. The righting arm obtained by the numerical simulation was approximately fitted by an analytical function. The irregular waves were decomposed into two Gauss stationary random processes, and the CARMA (2, 1) model was used to fit the spectral density function of parametric and forced excitations. The stochastic energy envelope averaging method was used to solve the PDFs and the probability. The validity of the semi-analytical method was verified by the Monte Carlo method. The C11 ship was taken as an example, and the influences of the system parameters on the PDFs and probability were analyzed. The results show that the probability of ship rolling is affected by the characteristic wave height, wave length, and the heading angle. In order to provide proper advice for the ship''s manoeuvring, the parametric excitations should be considered appropriately when the ship navigates in the oblique seas.     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.     

循环载荷作用下损伤船体梁极限强度研究

为了研究循环载荷作用下扶强材初始损伤对其极限强度的影响,进行了14组扶强材的循环加载试验和分析。构造了考虑材料累积损伤完整、断筋和大变形的扶强材单元极限承载力计算公式,提出了相应循环载荷作用下损伤扶强材单元的端缩曲线表达式和船体梁极限强度计算的简化逐步破坏法。编制了循环载荷作用下船体梁损伤极限强度计算程序,进行了船体梁极限强度计算,并与有限元结果进行对比。研究结果表明：改进的损伤扶强材模型可较为准确地描述扶强材材料损伤的完整、断筋和大变形的极限承载力退化情况,扶强材腹板断裂的损伤相较初始大变形及材料累积损伤形式承载力下降程度更明显;所提出的循环载荷作用下损伤船体梁极限强度计算的简化逐步迭代方法,能定量地计算扶强材在不同类型损伤下的极限承载力退化程度,具有较高精度,方便易行,可应用于工程设计。     

Study on Ship Automatic Berthing System with Mooring Lines

This study describes an automatic berthing system with mooring lines. It is designed to be berthed by using mooring device on the upper deck of a ship. It is to berth once maintaining parallel with the quay by controlling both forward and aft breast lines. Berthing method is used through length adjustment of mooring lines connected between ship and quay by controlling the angular velocity and the torque of hydraulic motor in mooring device. The study is conducted under three changing conditions of draft, such as even-keel, rise of the gravity center and trim to stern. Variables affecting berthing stability are determined based on the control performance of each condition. Bond graphs method is used to model the system. Controller is designed as PID control method of reference-model algorithm. The control program is composed of synchronous control system based on the equations derived with the numerical analysis. The tank test is conducted to verify the usefulness of the control program.