水流对浮体作用的SPH方法模拟

水流与浮体相互作用时,运动情况十分复杂,目前很难准确有效地模拟。而SPH作为一种纯拉格朗日方法,无需构建网格,用核函数近似粒子进行离散,能较好地解决一些自由面大变形问题。文章利用SPH法对溃坝时引起的高速水流冲击浮体以及水体晃动时破损浮体的运动过程进行模拟。模拟结果表明,SPH法能有效地进行水流对浮体作用的研究。     

Green water loading on a floating structure with degree of freedom effects

The aim of the present work is to investigate whether the degree of freedom (DOF) of a floating body has a notable effect on the maximum impact pressure due to green water on deck. The analysis is carried out for a box-shaped floating structure with a deckhouse, using experimental and numerical means to model the green water load. Green water on deck and impact on the deckhouse is generated by the impingement of a focusing wave group on a floating structure. Computations are performed using a two-dimensional constrained interpolation profile-based model solving the Navier–Stokes (N–S) equations with free surface boundary condition to deal with nonlinear water–structure interactions. The free surface is captured by a volume of fluid (VOF)-type tangent of hyperbola for interface capturing/slope weighting (THINC/SW), which is more accurate than the original THINC scheme. The verifications of the simulation through a series of model-to-model comparisons are performed in a two-dimensional glass-wall wave tank. Experimental water surface elevations, body motions and impact pressure are compared satisfactorily with the computed results for different DOFs cases. As a result, the peak impact pressure due to green water decreases rapidly with the increasing DOF.     

Hydrodynamics of a body floating in a two-layer fluid of finite depth. Part 2. Diffraction problem and wave-induced motions

A linearized two-dimensional diffraction problem in a two-layer fluid of finite depth was solved for a general floating body and relevant wave-induced motions were studied. In a two-layer fluid, for a prescribed frequency, incident waves propagate with two different wave modes. Thus the wave-exciting forces and resulting motions must be computed separately for each mode of the incident wave. The boundary integral equation method developed by the authors in the Part-1 article was applied to directly obtain the diffraction potential (pressure) on the body surface. With the computed results, an investigation was carried out on the effects of the fluid density ratio and the interface position on the wave-exciting forces on the body and the motions of the body, including the case in which the body intersects the interface. By a systematic derivation using Green's theorem, all the possible reciprocity relations were derived theoretically in explicit forms for a system of finite depth; these relations were confirmed to be satisfied numerically with very good accuracy. Experiments were also carried out using water and isoparaffin oil as the two fluids and a Lewis-form body. Measured results for the sway- and heave-exciting forces and the heave motion were compared with the computed results, and a favorable agreement was found.     

Finite-difference simulation of green water impact on fixed and moving bodies

A numerical method was developed to predict the water impact pressure caused by green water phenomena. The density function method was employed in the framework of a locally refined overlapping grid system. The simple problem of a rectangular body placed in regular waves was considered and the simulation results were compared with tank experiments. Good agreement with experiment was shown for wave–body interactions and for the pressure values at three different positions on the body. The case of a rectangular body with a vertical wall on the deck was also considered and the comparison showed satisfactory agreement. It was demonstrated that this method could be extended to a moving body problem, in which the body was free to undergo heave, pitch, and surge motions.     

浮体升沉和纵摇运动对甲板上浪载荷影响数值研究（英文）

In this paper, the influence of heave and pitch motions on green water impact on the deck is numerically investigated. The vessel motions are determined using a potential theory based method and provided as input to finite volume based CFD computations of green water phenomenon. A dynamic mesh approach is adopted to determine instantaneous body positioning in the fluid domain.Detailed validation studies with published experimental results for 2D and 3D fixed vessel cases are initially performed to validate the present numerical approach before studying the moving vessel problem. The results show that inclusion of heave and pitch motion changes the disturbed wave field near the bow which influences the free surface as well as the impact loading due to green water. The effect of wave steepness on green water impact is also investigated and it is seen that the present numerical method is capable of capturing green water load. It is observed that the effects of vessel motions on green water load are not negligible and one should consider this effect too. The incorporation of vessel motions in the vertical plane affects the green water loading on the deck.     

物体在浮体上的纵向运动引起的浮体垂向运动

基于切片法建立物体在静水中浮体上纵向运动时浮体的水动力特征及运动微分方程。物体对浮体的垂向作用力考虑为瞬时压力，采用4阶龙格库塔法求解浮体运动微分方程，并将计算值与试验数据进行比较。     

船体湿表面脉动压力预报的三维非线性方法

本文给出了一种预报船体湿表面脉动压力响应的三维非线性时历模拟方法.计及瞬时物面变化引起的各种非线性因素,在时域内建立船体非线性运动方程.在每一时刻,认为船体以其瞬时平均湿表面为平衡位置做简谐振荡运动,利用三维线性频域解来简化运算.在此基础上,进一步将绕射和辐射问题在静平衡位置求解,仅考虑占主导地位的入射波和静恢复力变化的非线性,作为一种兼顾工程需要的计算精度和计算效率的简化方法.以此方法为基础编制了计算机程序,对实船进行了计算,并与模型试验进行了对比.分析比较表明,该方法能够给出令人满意的船体湿表面脉动压力响应预报结果.     

溃坝问题的数值仿真和实验(英文)

In this paper, two novel numerical computation methods are introduced which have been recently developed at Research Institute for Applied Mechanics ( RIAM ), Kyushu University, for strongly nonlinear wave-body interaction problems, such as ship motions in rough seas and resulting green-water impact on deck. The first method is the CIP-based Cartesian grid method, in which the free surface flow is treated as a multi-phase flow which is solved using a Cartesian grid. The second method is the MPS method, which is a so-called particle method and hence no grid is used. The features and calculation procedures of these numerical methods are described. One validation computation against a newly conducted experiment on a dam break problem, which is also described in this paper, is presented.     

Two-dimensional numerical simulation and experiment on strongly nonlinear wave–body interactions

A constrained interpolation profile (CIP)-based Cartesian grid method for strongly nonlinear wave–body interaction problems is presented and validated by a newly designed experiment, which is performed in a two-dimensional wave channel. In the experiment, a floating body that has a rectangular section shape is used. A superstructure is installed on the deck and a small floating-body freeboard is adopted in order to easily obtain water-on-deck phenomena. A forced oscillation test in heave and a wave–body interaction test are carried out. The numerical simulation is performed by the CIP-based Cartesian grid method, which is described in this paper. The CIP scheme is applied in the Cartesian grid-based flow solver. New improvements of the method include an interface-capturing method that applies the tangent of hyperbola for interface capturing (THINC) scheme and a virtual particle method for the floating body. The efficiency of the THINC scheme is shown by a dam-breaking computation. Numerical simulations on the experimental problem for both the forced oscillation test and the wave–body interaction test are carried out, and the results are compared to the measurements. All of the comparisons are reasonably good. It is shown, based on the numerical examples, that the present CIP-based Cartesian grid method is an accurate and efficient method for predicting strongly nonlinear wave–body interactions.     

深海系泊浮体物面非线性时域耦合动力分析

文章基于三维时域势流理论和弹性细长杆理论,研究并提出了深海系泊浮体物面非线性时域耦合动力分析方法。该方法采用时域物面非线性理论方法在瞬态位置直接时域模拟系泊浮体所需水动力,结合有限元方法计算系泊缆索的动力响应,利用异步耦合方法实现浮体和系泊缆索的时域耦合动力求解。既满足系泊浮体时域水动力耦合,又满足系泊浮体和系泊缆索动力耦合。通过对二阶非线性不规则波作用下深海系泊半潜式平台的时域耦合响应特性进行研究,将不同海况下物面非线性时域耦合静力响应和动力响应与间接时域耦合动力响应的三种方法计算结果进行比较。研究结果表明,系泊缆索动力响应明显,平台瞬态空间位置对垂荡低频运动影响较大,有必要在平台瞬时湿表面采用动力响应方法进行深海系泊浮体时域耦合响应分析。     

系泊多浮体系统流固耦合和浮体间耦合动力分析

利用多刚体力学运动学方法进行系统铁数学描述。以系统摇荡时各浮体质心位移和浮体角位移表达广义速率。进而表达出有关偏速度、偏角速度、加速度和角加速度。单个系泊浮体在规则波作用下，所受的流体动压力用势流理论计算。多浮体系统摇荡时除入射规则波外，每一浮体的摇荡辐射波对其他浮体也构成了入射波并作用以流体动压力。表达出系统的广义主动力和广义惯性力后代入Kane方程，补充以浮体间的运动约束方程，消去时间因子后得到求解各浮体质心位移和浮体角位移的复数线性代数方程组。     

浮式直立圆柱体辐射问题研究

采用特征函数展开及边界匹配的方法求解流场各子区域速度势的展开系数,得到了浮式直立圆柱体在静水面作纵荡,垂荡和纵摇运动时引起的水动力系数,分析了所取截断项数对计算结果的影响,所需最少截断项数取决于水深半径比和吃水半径比。讨论了水深对水动力系数的影响,得到了可近似表示深水情况的水深半径比。     

17万吨浮船坞两段建造法水上合拢对接工艺设计

由于船台限制，17万吨浮船坞无法在船台整体合拢。经过技术论证，采取将坞体分为前后两段分别建造，将浮箱甲板以下端部密封后下水，再在水上进行浮态对接；利用对接工装将坞体对接缝处的局部区域与水隔离，再进行焊接。此工艺操作难度较大，特别是前后段坞体对中调整。针对对中难度大的问题，研发快速定位工装，并编制详细的合拢过程控制操作程序，确保水上对接合拢的顺利进行。同时为其它大型船坞或船舶两段建造、水上合拢积累经验，为其它船舶修理改装、水上定位作业施工提供一定的借鉴和指导作用。     

紧致插值曲线CIP方法及其应用

波浪与建筑物的相互作用过程会涉及到波浪破碎、水气掺混和结构物的动力响应等复杂过程,对数值算法提出了更高要求。文章基于紧致插值曲线CIP（Constrained Interpolation Profile）方法建立了可模拟波浪破碎、翻滚等自由面大变形流动问题的数学模型。模型以CIP方法为流场基本求解器,离散了纳维—斯托克斯（Navier-Stokes：N-S）方程,同时还以CIP方法捕捉了自由面,通过多相流理论描述了流—固—气之间的相互作用。对强非线性自由表面流动问题的典型算例溃坝问题开展了数值模拟,并通过与他人实验结果的比较验证了模型的有效性。最后开展了极端波浪对浮式结构冲击过程的模拟,准确地预测了甲板上浪和结构的动力响应等问题。     

全球FPSO船型开发领域核心技术的演化分析

以国内外专利数据为分析对象,采用专利计量分析的技术方法深入研究浮式生产储油卸油装置(FPSO)开发技术的演化态势,对FPSO船型开发领域的核心专利进行分析,总结出FPSO船型开发的发展路线.结果表明:当前FPSO船型开发领域的主要研究方向是旧船改造、新概念FPSO和通用FPSO.旧船改造的技术点主要聚焦于超大型船舶设计;新概念FPSO的技术点主要集中于浮体形状设计、上部甲板大型化和刚性立管作业问题;通用FPSO的技术点主要侧重于标准化船壳设计、模块化FPSO设计和FPSO一体化建造,未来企业可考虑在模块化FPSO的上甲板承载与布局、管线布置、生活区布置、减摇和抑制垂荡问题上进行专利布局研究.     

Three-dimensional numerical analysis of shipping water onto a moving ship using a particle method

The objective of this study was to develop a numerical analysis method based on the moving particle semi-implicit method for simulating shipping water on a moving ship. Towing tests of a very large crude carrier were numerically analyzed for three typical wavelengths. The ship was forced to move in order to express previously measured ship oscillations, and the calculated fluid behavior and the impact pressure on the deck were compared with the experimental results.     

Numerical simulation of the spreading dynamic responses of the multibody system with a floating base

To simulate the dynamic responses of the multibody system with a floating base when the upper parts spread with a certain sequence and relative speed, the homogeneous matrix method is employed to model and simulate a four-body system with a floating base and the motions are analyzed when the upper parts are spread sequentially or synchronously. The rolling, swaying and heaving temporal variations are obtained when the multibody system is under the conditions of the static water along with the wave loads and the mean wind loads or the single pulse wind loads, respectively. The moment variations of each joint under the single pulse wind load are also gained. The numerical results showed that the swaying of the floating base is almost not influenced by the spreading time or form when the upper parts spread sequentially or synchronously, while the rolling and the heaving mainly depend on the spreading time and forms. The swaying and heaving motions are influenced significantly by the mean wind loads. The single pulse wind load also has influences on the dynamic responses. The torque of joint 3 and joint 4 in the single pulse wind environment may be twice that in the windless environment when the system spreads with 60 s duration.     

二维甲板上浪模型试验研究

甲板上浪属于非常复杂和强非线性的波体相互作用过程。研究甲板上浪的基本方法包括理论分析、模型试验和数值计算。本文利用浪高仪测试方法对特定实验工况捕捉甲板上浪特征参数信息,分析甲板上浪的水位高度、上浪量与浮态之间的关系。利用PIV方法捕捉了波浪沿模型首端爬行、变形、破碎的过程_[1]。获得了上浪的统计规律,有助于自由液面的非线性流体力学模型试验技术的发展,推动相应学科的进步和发展。     

内河汽车滚装船活动滚装区域的水喷淋系统研究

升降式甲板作为一种滚装船舶专用的特种设备,适用于需要在装载小型车与大型车间切换的滚装船舶。配置升降式甲板的船舶可根据装载需求,灵活切换装车甲板状态,提高运输效率。然而,在升降式甲板底部设置的传统压力水雾管路与消防总管的固定连接会对甲板升降造成阻碍,需要拆除部分管路并在甲板位置调整后重新连接。这无论从安全性和便利性来看均极为不利,一方面无法保证管路每次拆装后的密性,同时每次改变升降式甲板状态都需要拆装大量管路。该文主要提出侧向压力水雾喷淋的设计方案,用于解决升降式甲板的压力水雾消防问题。     

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

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