Numerical analysis on ship maneuvering coupled with ship motion in waves

This paper considers a numerical analysis of ship maneuvering performance in the presence of incident waves and resultant ship motion responses. To this end, a time-domain ship motion program is developed to solve the wave–body interaction problem with the ship slip speed and rotation, and it is coupled with a modular-type 4-DOF maneuvering problem. In this coupled problem, the second-order mean drift force, which can play an important role in the ship maneuvering trajectory, is estimated by using a direct pressure integration method. The developed method is validated by observing the second-order mean drift force, and planar trajectories in maneuvering tests with and without the presence of incident waves. The comparisons are made for two ship models, Series 60 with block coefficient 0.7 and the S-175 containership, with existing experimental data. The maneuvering tests observed in this study include a zig-zag test in calm water, and turning tests in calm water and in regular waves. The present results show a fair agreement of overall tendency in maneuvering trajectories.     

A time-domain prediction method of ship motions

A time-domain analysis is used to predict wave loading and motion responses for a ship traveling at a constant speed in regular oblique waves. Considered as a distribution of normal velocities on the wetted hull surface, the combined diffraction and radiation perturbations caused by the forward moving ship and her motions are determined simultaneously. This way, the ship-hull boundary condition is exactly fulfilled. The 3-D time domain Green's function is used to express the combined diffraction/radiation potential in terms of impulsive and memory potentials. Application of the Bernoulli equation yields the pressure distribution and accordingly, the necessary hydrodynamic forces. The equations of motion of the ship are then developed and solved in the time domain.Forces and motions at forward speed are predicted for a Wigley ship-hull in head waves and for a catamaran-ferry in oblique waves. Comparison is made with published theoretical and experimental results for the Wigley ship-hull, and the agreement is good. For the catamaran, a self-propelled model is built and tested both in a large towing tank and in a seakeeping basin in order to measure the six-degrees-of-freedom forces, moments and motions at forward speed in regular waves of different directions. For the longitudinal motions, the agreement between measurements and predictions is generally good. For the transverse motions, however, acceptable discrepancy exists. The discrepancy is thought to be mainly due to the exclusion from the analysis of the rudder forces and viscous damping. The inclusion of such nonlinear effects in the time domain simulation involves complex analysis and this problem is left to a future research.     

Investigation of ship hydroelasticity

The importance of hydroelastic analysis of large and flexible container ships of today is pointed out. A methodology for investigation of this challenging phenomenon is drawn up and a mathematical model is worked out. It includes definition of ship geometry, mass distribution, structure stiffness, and combines ship hydrostatics, hydrodynamics, wave load, ship motion and vibrations. Based on the presented theory, a computer program is developed and applied for hydroelastic analysis of a flexible segmented barge for which model test results of motion and distortion in waves have been available. A correlation analysis of numerical simulation and measured response shows quite good agreement of the transfer functions for heave, pitch, roll, vertical and horizontal bending and torsion. The tool checked in such a way can be further used for reliable hydroelastic analysis of ship-like structures.     

Prediction of design wave loads of the ocean structure by equivalent irregular wave approach

The present study is employing the equivalent irregular wave approach to predict the wave loads for a ship encountering the worst sea state with respect to the critical dynamic loading parameter. Two different hydrodynamic numerical models, i.e. 3D pulsating source technique and 3D translating pulsating source technique, are applied to calculate the corresponding RAO of the ship moving in waves. Incorporating the RAO of the related physical properties, we can calculate the extreme value for the corresponding ship loading factor, which can be regarded as the worst sea state in the service lifetime of the ship. With the time and period of the occurrence of the corresponding extreme value, we can simulate the time history of the wave load in this period, which is so-called equivalent irregular wave approach. Comparing with the results calculated by the traditional equivalent regular wave approach, we find that the equivalent irregular wave approach can simulate the corresponding wave load more realistic, especially for dynamic pressure. Using the equivalent irregular wave approach can offer the effective and practical base for the ship structural analysis.     

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.     

Experimental investigation of hull girder vibrations of a flexible backbone model in bending and torsion

Although the coupled horizontal–torsional vibrations of open ships have been investigated numerically for decades, the available experimental data in oblique seas seem rare. Model tests, considering natural frequencies of bending and torsional modes, have been conducted by the Centre for Ships and Ocean Structures (CeSOS) in the towing tank and ocean basin. A flexible backbone model was designed with five cut-outs on the top side of the aluminum beam to realize approximately torsional stiffness as well as vertical and horizontal bending stiffness. This paper mainly deals with measured bending and torsional vibrations in regular and irregular waves. The damping ratios, mode shapes, and modal moments were derived from experimental data, and a numerical model based on modal superimposition is established according to the measured hydroelastic properties. Some time-domain simulations are carried out considering the structural characteristics, and compared with measured results. The test data in regular waves are presented to show the possible factors of influence on the vibrations. The effect of bending and torsional vibrations on the extreme response values in irregular waves is estimated. The uncertainties in the experiments are discussed and conclusions are presented at the end of this paper.     

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.     

Coupled analysis of nonlinear sloshing and ship motions

A coupled numerical model considering nonlinear sloshing flows and the linear ship motions has been developed based on a boundary element method. Hydrodynamic performances of a tank containing internal fluid under regular wave excitations in sway are investigated by the present time-domain simulation model and comparative model tests. The numerical model features well the hydrodynamic performance of a tank and its internal sloshing flows obtained from the experiments. In particular, the numerical simulations of the strong nonlinear sloshing flows at the natural frequency have been validated. The influence of the excitation wave height and wave frequency on ship motions and internal sloshing has been investigated. The magnitude of the internal sloshing increases nonlinearly as the wave excitation increases. It is observed that the asymmetry of the internal sloshing relative to still water surface becomes more pronounced at higher wave excitation. The internal sloshing-induced wave elevation is found to be amplitude-modulated. The frequency of the amplitude modulation envelope is determined by the difference between the incident wave frequency and the natural frequency of the internal sloshing. Furthermore, the coupling mechanism between ship motions and internal sloshing is discussed.     

Verification of A Numerical Harbour Wave Model

A numerical model for wave propagation in a harbour is verified by use of physical models.The extended time-dependent mild slope equation is employed as the governing equation,and the model is solved by use of ADI method containing the relaxation factor.Firstly,the reflection coefficient of waves in front of rubble-mound breakwaters under oblique incident waves is determined through physical model tests,and it is regarded as the basis for simulating partial reflection boundaries of the numerical model.Then model tests on refraction,diffraction and reflection of waves in a harbour are performed to measure wave height distribution.Comparative results between physical and numerical model tests show that the present numerical model can satisfactorily simulate the propagation of regular and irregular waves in a harbour with complex topography and boundary conditions.     

港域波浪数学模型的改进与验证

通过物理模型对改进的港内波浪传播变形数学模型进行验证。该数学模型以推广的时变缓坡方程为控制方程,采用含松弛因子的ADI法求解,并对波浪反射和透射边界模拟方法进行改进。先通过物理模型试验确定斜向浪入射条件下抛石防波堤前的波浪反射系数,作为数学模型中部分反射边界模拟的依据。然后进行了一个典型港口内波浪折射、绕射和反射的模型试验,测量港内波浪分布。对比模型试验和数学模型计算的结果表明,数学模型可较好地模拟港内复杂地形和边界条件下规则波和不规则波的传播变形。     

Sensitivity of extreme hydroelastic load effects to changes in ship hull stiffness and structural damping

A new hybrid method for the time-domain nonlinear simulation of the hydroelastic load effects and the peak over threshold (POT) method for the calculation of the short-term extreme responses are briefly described and applied to a flexible SL-7 class containership and a flexible liquefied natural gas (LNG) ship. Three stiffness levels, three stiffness distributions and three modal damping ratios are used to study the influence of the hull flexibility and structural damping on the short-term prediction of extreme vertical hydroelastic load effects. The results give justification for some simplified treatment of the first vertical flexible mode in early design stage when structural details are not available.     

Dynamic response of a monopile wind turbine in waves: Experimental uncertainty analysis for validation of numerical tools

The responses of a monopile offshore wind turbine subjected to irregular wave loads are investigated numerically and experimentally, considering a range of sea states. An extensive experimental campaign was carried out on a fully flexible model, representative of a 5 MW offshore wind turbine, at 1:40 scale. An assessment of the experimental results for the response amplitude operator for regular waves and the 90th percentile seabed bending moment in long-crested irregular waves is carried out using two models (analytical and numerical) for uncertainty propagation, suggesting that bias errors in the model properties and in the wave elevation contribute the most to the total uncertainty. The experimental results are also compared to a numerical model using beam elements and Morison-type wave loads with second order wave kinematics. The numerical model does not capture all of the responses within the level of uncertainty of the experiments, and possible reasons for the discrepancies are discussed.     

Numerical simulation of irregular wave overtopping against a smooth sea dike

Based on the filtered Navier-Stokes equations and Smagorinsky turbulence model,a numerical wave flume is developed to investigate the overtopping process of irregular waves over smooth sea dikes.Simulations of fully nonlinear standing wave and regular wave’s run-up on a sea dike are carried out to validate the implementation of the numerical wave flume with wave generation and absorbing modules.To model stationary ergodic stochastic processes,several cases with different random seeds are computed for each specified irregular wave spectrum.It turns out that the statistical mean overtopping discharge shows good agreement with empirical formulas,other numerical results and experimental data.     

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.     

Oblique Wave-Induced Responses of A VLFS Edged with A Pair of Inclined Perforated Plates

This paper is concerned with the hydroelastic responses of a mat-like, rectangular very large floating structure (VLFS) edged with a pair of horizontal/inclined perforated anti-motion plates in the context of the direct coupling method. The updated Lagrangian formulae are applied to establish the equilibrium equations of the VLFS and the total potential formula is employed for fluids in the numerical model including the viscous effect of the perforated plates through the Darcy''s law. The hybrid finite element-boundary element (FE-BE) method is implemented to determine the response reduction of VLFS with attached perforated plates under various oblique incident waves. Also, the numerical solutions are validated against a series of experimental tests. The effectiveness of the attached perforated plates in reducing the deflections of the VLFS can be significantly improved by selecting the proper design parameters such as the porous parameter, submergence depth, plate width and inclination angle for the given sea conditions.     

Hydroelastic analysis of floating structures with liquid tanks and comparison with experimental tests

In this study, we develop a numerical method for a 3D linear hydroelastic analysis of floating structures with liquid tanks subjected to surface regular water waves and compare the numerical results with experimental tests. Considering direct couplings among structural motion, sloshing, and water waves, a mathematical formulation and a numerical method extended from a recent work [1] are developed. The finite element method is employed for the floating structure and internal fluid in tanks, and the boundary element method is used for the external fluid. The resulting formulation completely incorporates all the interaction terms including hydrostatic stiffness and the irregular frequency effect is removed by introducing the extended boundary integral equations. Through various numerical tests, we verify the proposed numerical method. We also performed 3D hydroelastic experimental tests of a floating production unit (FPU) model in an ocean basin. The measured dynamic motions are compared with the numerical results obtained using the proposed method.     

Dynamics of a moored barge under regular and random waves

Developments in the study of wave forces and construction techniques in deep water by the offshore oil industry have increased the use of marine terminals at deep water locations. A thorough understanding of moored ship dynamics when subjected to waves, wind and current combined with the use of flexible mooring lines would help to design berthing terminals for exposed areas. In this paper, the three dimensional problem of wave interactions with a barge moored to a single point is dealt with, based on the finite element method. The effect of flexibility of the mooring line and the point of mooring on the response of the barge as well as the mooring line tension is investigated. The paper compares the numerical results with model tests carried out on a barge moored to a fixed support under regular and random waves in head sea. The effect of stiffness of the mooring line on the barge response for different mooring points is discussed, which would be useful for the designers. The effect of viscous damping is also considered. The analytical results are in good agreement with the experimental results in both regular and random waves.     

Wave Forces Acting on Vertical Walls

Regular and irregular wave forces acting on vertical walls are studied by a previously developed numerical model. The computed wave forces are compared with the available experimental data to verify the numerical model, and satisfactory agreements are obtained. The variation of wave forces with incident angles and the shape of simultaneous pressure distribution are investigated, and the comparisons between numerical results and Goda＇ s predictions are also carried out. It is concluded that the maximum wave forces acting on the unit length of vertical wall is often induced by the obliquely incident waves instead of normally incident waves, while Goda＇ s formula may be inapplicable for oblique wave incidence. The shape of simultaneous pressure distribution is not significantly influenced by incident angles, and it can be favorably predicted by Goda＇ s formula. When regular wave heights are taken as the same as irregular wave height H1%, the irregular wave forces Ph. 1% are slightly larger than regular wave forces in most cases.     

动力定位船舶二阶低频慢漂力模型试验研究

对一艘动力定位船舶二阶低频慢漂力进行了模型试验,并将试验得到的纵向慢漂力谱、横向慢漂力谱与势流理论方法得到的理论值进行比较,结果表明,该模型试验方法与理论计算较为吻合。可为动力定位系统的设计和应用提供参考。