高速滑行艇在规则波中运动的RANSE数值模拟(英文)

This paper presents a study on the numerical simulation of planing crafts sailing in regular waves. This allows an accurate estimate of the seas keeping performance of the high speed craft. The simulation set in six-degree of freedom motions is based on the Reynolds averaged Navier Stokes equations volume of fluid (RANSE VOF) solver. The trimming mesh technique and integral dynamic mesh method are used to guarantee the good accuracy of the hydrodynamic force and high efficiency of the numerical simulation. Incident head waves, oblique waves and beam waves are generated in the simulation with three different velocities (Fn =1.0, 1.5, 2.0). The motions and sea keeping performance of the planing craft with waves coming from different directions are indicated in the flow solver. The ship designer placed an emphasis on the effects of waves on sailing amplitude and pressure distribution of planing craft in the configuration of building high speed crafts.     

Nonlinear dynamics of a submerged floating moored structure by incremental harmonic balance method with FFT

A submerged floating moored structure has a great potential in ocean engineering applications. The nonlinear dynamics of a submerged floating moored structure subjected to vertical excitation with possible slackness in the mooing system are investigated by incremental harmonic balance (IHB) method. Heaviside step function is introduced to describe the nonlinearity in axial stiffness arising from loss of mooring tension. The dimensionless governing equation is derived, and three parameters, frequency ratio η, damping ratio ζ and dimensionless net buoyancy W, are found to be independent. Due to the fact that the restoring force term is function of the unknown displacement and could barely be expressed in an explicit form of time, a fast Fourier transformation (FFT) is implemented in IHB method to simplify the Galerkin average procedure. Both stable and unstable solutions and both period-1 and bifurcated solutions are obtained by IHB method. The stability of the periodic solutions is investigated by Floquet theory. Parameter study is carried out. Results indicate that the system nonlinearity becomes stronger as dimensionless the net buoyancy W and damping ratio ζ decrease. A path to chaotic motions though a series of period doubling bifurcations is found. Multiple solutions are observed, and the domains of attraction are investigated by interpolated cell mapping (ICM) technique.     

Validation of a dynamic mooring model coupled with a RANS solver

Standard design procedures and simulation tools for marine structures are aimed primarily for use by the offshore oil and gas. Mooring system restoring forces acting on floating offshore structures are obtained from a quasi-static mooring model alone or from a coupled analysis based on potential flow solvers that do not always consider nonlinear mooring-induced restoring forces, fluid structure interactions, and associated hydrodynamic damping effects. This paper presents the validation of a dynamic mooring system analysis technique that couples the dynamic mooring model with a Reynolds-averaged Navier-Stokes (RANS) equations solver. We coupled a dynamic mooring model with a RANS equations solver, and analyzed a moored floating buoy in calm water, regular and irregular waves and validated our motion and mooring force predictions against experimental measurements. The mooring system consisted of three catenary chains. The analyzed response comprised decaying oscillating buoy motions, linear and quadratic damping characteristics, and tensile forces in mooring lines. The generally favorable comparison of predicted buoy motions and mooring forces to experimental data confirmed the reliability of our implemented coupling technique to predict system response. Additional comparative results from a potential flow solver demonstrated the benefits of the coupled dynamic mooring model with RANS equations. The successful validated tool of coupling the dynamic mooring model with the RANS solver is available as open source, and it shows the potential of the coupled methodology to be used for analyzing the moored offshore structures.     

Classification criterion of narrow/wide ice-resistant conical structures based on direct measurements

Ice-induced structural vibration generally decreases with an increase in structural width at the waterline. Definitions of wide/narrow ice-resistant conical structures, according to ice-induced vibration, are directly related to structure width, sea ice parameters, and clearing modes of broken ice. This paper proposes three clearing modes for broken ice against conical structures: complete clearing, temporary ice pile up, and ice pile up. In this paper, sea ice clearing modes and the formation requirements of dynamic ice force are analyzed to explore criteria determining wide/narrow ice-resistant conical structures. According to the direct measurement data of typical prototype structures, primary quantitative criterion of the ratio of a cone width at waterline (D) to sea ice thickness (h) is proposed. If the ratio is less than 30 (narrow conical structure), broken ice is completely cleared and a dynamic ice force is produced; however, if the ratio is larger than 50 (wide conical structure), the front stacking of broken ice or dynamic ice force will not occur.     

Temporal and spectral characteristics of seismic ground motions: Offshore versus onshore

To distinguish offshore and onshore seismic ground motions, conventional analyses in terms of peak ground acceleration (PGA) and earthquake response spectrum (ERS) have been carried out in a recent work by authors and other papers in literature. In the present study, distinct temporal and spectral characteristics between onshore and offshore earthquake ground motions are further investigated in time-domain and frequency-domain. The data used is 69 pairs of concurrent onshore and offshore ground motions collected from the Kyoshin Network (K-NET). Each pair of data are of approximately identical epicenter distances. Comparisons are made on zero-up-crossing period (T_z), peak-to-trough acceleration range (A_pt) and period (T_pt), duration of ground motion (T_d), predominant frequency (f_p) and the spectral bandwidth parameter (ε). The results indicate that for offshore horizontal and vertical seismic signals, statistics of T_z, T_pt, T_d and predominant period T_p tend to be larger than the onshore counterparts. Meanwhile, ε of the offshore vertical ground motions is also greater. Through the proposed energy ratio (ER) analysis, the spectral energy of offshore ground motion is found to shift to moderate and low frequency bands. The time-frequency analysis conducted by Hilbert-Huang transform (HHT) shows that the Hilbert spectra of offshore accelerations contain larger spectral energy than the onshore counterpart but the corresponding instantaneous frequencies at peak energy are smaller, especially for horizontal recordings. Therefore, larger dynamic response of offshore structures is prone to be induced by the offshore ground excitation. This is further validated through the dynamic analysis of a marine pipeline in case study.     

升船机提升系统中船厢、水体和钢缆相互耦合作用分析

本文在作者[1]先前采用水弹性理论建立的对整个提升系统的动力特性进行数值研究的理论模型的基础上,首先对船厢变形状态下的缆力表达式作了更为详细的推导,然后以中国船舶科学研究中心的1∶30的缩比模型为对象,从整体角度,对其流固耦合问题,主要是提升过程中船厢内水体的晃荡、提升缆的受力和船厢动态响应,作了数值研究.在提升系统受到船厢内一非对称水波扰动的假定下,分别对船厢简化为刚性和考虑实际弹性的两种情况进行了计算、分析和对比,得出了一些有益的结论.计算结果表明,对升船机这样一个复杂系统,文献[1]中所建立的理论模型及本文的数值计算方法是可行的.     

3D numerical modeling of wave forces on tandem fixed cylinders using the BEM

In this paper a 3D numerical model was developed to study the complicated interaction between waves and a set of tandem fixed cylinders.The fluid was considered to be inviscid and irrotational.Therefore,the Helmholtz equation was used as a governing equation.The boundary element method（BEM） was adopted to discretize the relevant equations.Open boundaries were used in far fields of the study domain.Linear waves were generated and propagated towards tandem fixed cylinders to estimate the forces applied on them.Special attention was paid to consideration of the effect on varying non-dimensional cylinder radius and distance between cylinders,ka and kd on forces and trapped modes.The middle cylinder wave forces and trapped modes in a set of nine tandem cylinders were validated utilizing analytical data.The comparisons confirm the accuracy of the model.The results of the inline wave force estimation on n tandem cylinders show that the critical cylinder in the row is the middle one for odd numbers of cylinders.Furthermore the results show that the critical trapped mode effect occurs for normalized cylinder radiuses close to 0.5 and 1.0.Finally the force estimation for n tandem cylinders confirms that force amplitude of the middle cylinder versus normalized separation distance fluctuates about that of a single cylinder.     

粘性流体晃荡与弹性结构的相互耦合作用

提出了处理粘性流体晃荡与弹性结构相互耦合作用的理论及相应数值计算方法.流体运动采用N-S方程描述,控制方程采用有限差分法离散,并由超松弛迭代法求解.液体自由表面通过流体体积法进行重构.为了考虑液舱结构变形对液体晃荡的影响,应用了FAVOR技术.由此建立了描述液体晃荡与结构相互耦合作用的水弹性理论,并由相应的计算方法进行了系列运算,以考察二维液舱不同结构刚度对液体晃荡的影响.计算结果有助于进一步理解此类问题的物理现象.     

Dynamic motions of planing vessels in head seas

The dynamic response of planing vessels in regular head seas is investigated numerically. Nonlinear time domain simulations were performed using a 2D + t theory (two-dimensional plus time dependent theory). A prismatic hull form was assumed. We employed a two-dimensional (2D) boundary element method to solve the initial boundary value problems in 2D cross planes, in which nonlinear free-surface conditions and exact body boundary conditions were satisfied. At each time step, the total force and moment on the hull could be obtained by using the sectional forces calculated in those 2D planes. Heave and pitch motions were then acquired by solving the equations for those motions. The calculated heave and pitch responses were compared with the experiments by Fridsma (A systematic study of the rough-water performance of planing boats. Davidson Laboratory Report R-1275, 1969) for two different Froude numbers. Three-dimensional (3D) corrections at the transom stern were applied to show the influence of the 3D effect at the stern on the numerical results. Ship motions were affected by the 3D corrections, especially near the resonance frequency, while the phase angles were slightly affected and the acceleration peaks at the bow near the resonance frequency were sensitive to the 3D corrections. Other error sources in the theoretical results are also mentioned.     

高速三体船波浪中运动与增阻CFD计算研究

基于数值波浪水池技术,对波浪中高速三体船运动及增阻进行CFD计算研究。控制方程—RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波。首先对规则波顶浪中单个主船体的运动和增阻进行了计算,并与高速细长体理论计算结果进行了比较;随后进行了三体船运动和增阻计算,分析了侧片体对主船体阻力增加的影响。为高速三体船耐波性研究提供了数值工具。     

On the slack phenomena and snap force in tethers of submerged floating tunnels under wave conditions

Under severe sea wave conditions, the mooring tethers of submerged floating tunnel (SFT) might go slack. It may cause the structure failure during the service lifetime of SFT. The paper investigated SFT dynamics when going through tether slacking and the related snap force under wave conditions. Besides the nonlinearity of fluid drag and of structural geometry for a relative large structure displacement, the problem is characterized by the nonlinearity due to the discontinuity in axial stiffness of the tethers. To include these nonlinearities, the method of Lagrange energy is used to build the governing equations of SFT motion, and a bilinear oscillator is introduced to simulate the mooring tether operating in an alternating slack-taut state. The sensitivities of the occurrence of tether slacking to wave height and wave period are investigated. Results show that at a large wave height SFT tether will go slack and snap force occurs. SFT responses are categorized into three types of state according to the dynamic response characteristics of tether tension. Effects of two fundamental structure parameters, buoyancy-weight ratio (BWR) and inclined mooring angle (IMA), on the dynamic responses of SFT are analyzed. A slack-taut map of SFT tethers is built. It intuitively describes the occurrences of slack and snap force with different combinations of the two parameters. An analytical approach for slack prediction by deriving the slack criterion is provided to reveal the mechanism of the presented slack-taut map. By present research, the authors tried to make their effort to provide an alternative philosophy for SFT structural design on concerning preventing the occurrence of tether slacking and snap force.     

Riser-soil interaction model effects on the dynamic behavior of a steel catenary riser

A mathematical model employed to analyze the global dynamics of a Steel Catenary Riser (SCR) taking into account the interaction with the seafloor and the effect of the soil reaction forces is established. The choice of soil model plays an important role for the dynamic behavior of SCRs. The riser has been modeled using flexible beam with large curvature and elastic foundation beam to describe the riser-soil interaction by means of realistic nonlinear load-deflection (P–y) curves. The study is made to improve an existing finite element numerical code for dynamic analysis of mooring lines and risers, known as CABLE 3D, which is based on a slender rod assumption. Effects of nonlinear seabed model on the dynamic behavior of SCRs under vessel cyclic perturbation have further been investigated and discussed using a realistic P–y curve to simulate soil deformation and resistance forces. The interaction model depends on several factors, such as soil strength, penetration depth and riser characteristics. The dynamic responses of the riser Touchdown Point (TDP) excited by vessel periodic heave motion are studied and the results are compared with those from the linear spring model. SCR has been perturbed by 10 regular sinusoidal cycles and the responses calculated by improved code show a number of features such as suction force mobilization, gradual increasing penetration depth, and gradual reduction of soil resistance at maximum penetration. The riser behavior at the touchdown zone (TDZ) depends on the riser top motion amplitude, nonlinear soil stiffness and suction force. The impact of the riser-soil interaction model on the dynamic behavior in the TDZ has been thoroughly studied in this paper.     

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)     

桁架式Spar平台与系泊/立管系统的全时域非线性耦合动态分析

开发了对浮式平台系统进行耦合动态分析的全时域程序。采用二阶时域方法计算水动力荷载,在此方法中,对物面边界条件和自由水面边界条件进行泰勒级数展开,利用Stokes摄动展开分别建立相应的一阶、二阶边值问题,而且此边值问题的计算域不随时间变化。采用高阶边界元方法计算每一时刻流场中的速度势,利用四阶预报校正法对二阶自由水面边界条件进行数值积分。在自由表面加入一个人工阻尼层来避免波浪的反射。对于系泊缆索/立管/张力腿的动力分析,在一个总体坐标系中对控制方程进行描述,采用基于细长杆理论的有限元方法进行求解。在耦合动态分析中,采用Newmark方法对平台和系泊缆索/立管/张力腿的运动方程同时进行求解。利用开发的耦合分析程序对一个桁架式Spar平台的运动响应进行了数值模拟,给出了平台的位移和系泊缆索/立管上端点的张力,并得到了一些重要结论。     

水下拖缆抛物构型分析研究(英文)

This article discusses the dynamic state analysis of underwater towed-cable when tow-ship changes its speed in a direction making parabolic profile path. A three-dimensional model of underwater towed system is studied. The established governing equations for the system have been solved using the central implicit finite-difference method. The obtained difference non-linear coupled equations are solved by Newton's method and satisfactory results were achieved. The solution of this problem has practical importance in the estimation of dynamic loading and motion, and hence it is directly applicable to the enhancement of safety and the effectiveness of the offshore activities.     

Analyzing parabolic profile path for underwater towed-cable

This article discusses the dynamic state analysis of underwater towed-cable when tow-ship changes its speed in a direction making parabolic profile path. A three-dimensional model of underwater towed system is studied. The established governing equations for the system have been solved using the central implicit finite-difference method. The obtained difference non-linear coupled equations are solved by Newton＇s method and satisfactory results were achieved. The solution of this problem has practical importance in the estimation of dynamic loading and motion, and hence it is directly applicable to the enhancement of safety and the effectiveness of the offshore activities.     

气浮结构的运动特性研究

文章通过理论分析和试验研究探讨了气浮结构的运动特性。按单自由度刚体弹簧体系建立了气浮结构的升沉和摇摆运动方程,进而求气浮结构运动的固有周期,并对气浮结构的运动特性进行了讨论分析。文中引进的气浮力折减系数ka体现了气浮体的恢复力刚度系数与普通浮体的差异,该参数对气浮体的浮态、稳性、运动特性都有很大的影响。理论分析和试验表明气体浮体的固有周期比相应实浮体的大。结合理论和试验,确定了气浮结构作升沉运动时附连水质量系数的取值。摇摆运动的附连水质量系数的确定较为复杂,还需作进一步的研究。     

斜壁结构的水动力收敛方法研究

在分析波浪中的船舶运动或者计算大型结构的水动力系数时,往往采用时域格林函数方法。时域格林函数的一个重要局限性在于它在计算具有倾斜壁面的水动力系数时,结果很容易发散。时域格林函数本身的奇异性以及高频振动特性显然是水动力发散的一个重要原因。但即使该奇异性通过增加粘性以及表面张力的方式使之消失,计算具有斜壁结构的水动力时,发散现象依然存在。因此,该文提出一种滤波方法,除去时域格林函数的高频部分,留下其低频部分,并定义滤波系数,从而使作用于斜壁结构上的水动力值收敛。通过文中结果与频域兴波格林函数法的比较来确定最优滤波系数。结果表明：最优滤波系数几乎与运动幅值无关,但会受到运动频率以及物面形状的影响。     

Numerical and experimental motion simulations of nonsymmetric ship collisions

A calculation model to simulate nonsymmetric ship collisions, implying an arbitrary impact location and collision angle, is described in the paper. The model that is introduced is based on the time integration of twelve equations of motion, six for each ship. The motions of the ships are linked together by a mutual contact force. The contact force is evaluated as an integral over the surface tractions at the contact interface. The calculation model provides full time histories of the ship motions and the acting forces. Physical understanding of the underlying phenomena was obtained by a series of model-scale experiments in which a striking ship collided with an initially motionless struck ship. In this paper, numerical simulations of four nonsymmetric collisions are presented and the calculations are validated with the results of the experiments.