涡激振动影响下细长海洋结构物的阻力特性

利用柔性圆柱体涡激振动水池模型试验对细长海洋结构物发生涡激振动时的阻力特性进行研究。利用试验测得的圆柱体模型顺流向平均弯曲应变,结合梁复杂弯曲理论,计算出模型各截面处的平均阻力,获得沿模型轴向分布的阻力系数,分析涡激振动影响下细长海洋结构物阻力系数的变化特性。研究结论表明:均匀流下细长海洋结构物的阻力为非均布载荷,涡激振动对阻力系数具有放大效应,导致阻力系数在1.4~2.5之间变化。     

柔性立管涡激振动响应轨迹特性研究

为了深入研究细长柔性立管的涡激振动响应特性,进行了柔性立管的拖曳水池试验。由拖车拖动立管产生相对来流,根据应变测试得到的应变数据,基于模态叠加法得到位移响应。试验分析前,通过数值方法先针对刚性立管的涡激振动响应轨迹特性进行了分析。紧接着,通过试验方法对柔性立管的单模态以及多模态涡激振动响应轨迹特性进行了深入的分析和讨论。通过分析发现:柔性立管在低速下具有与刚性立管类似的轨迹响应特性,均呈现经典的8字形状;柔性立管在高速下,其轨迹开始变得混乱,这主要是由位移的多模态响应特性所产生。     

基于OpenFOAM低质量比圆柱体双自由度涡激振动的数值研究

论文采用RANS方法对低质量比弹性支撑圆柱体的双自由度涡激振动进行了数值模拟,圆柱体的质量比为2.6,固有频率比为1,雷诺数范围为7300~55000,相应的约化速度范围2~15,包括了经典试验中出现的整个锁定范围。文中通过数值模拟得到了不同的响应分支,重点分析了各响应分支的振动特性,与相应试验的结果吻合较好;捕捉到对应于各分支的不同典型尾涡模型,同时还发现了2T到2P的过度态;伴随约化速度的改变,振动轨迹也相应地呈现出顺时针、逆时针以及偏向来流的上游或下游方向等多种不同的形式,在超上端分支范围内,振动轨迹由"8"字形逐渐向"新月"形转变。     

计及海底井口影响的隔水管涡激振动响应特性试验

在靠近海底井口位置,隔水管周围流场受井口系统的影响发生变化,流固耦合下涡激振动可能诱发隔水管在横流向(CF方向)更为剧烈的振动。为了研究受管土装置影响的隔水管涡激振动响应特性,使用自主设计的管土装置模拟海底井口,采用8 m柔性立管,进行了均匀流下单管和受管土装置影响的涡激振动对比试验。试验通过FBG光纤应变片测得应变信息,使用模态叠加法、FFT变换处理分析试验数据,对比分析2组隔水管的主导频率、应变时历与幅值谱、无因次振幅以及激励力系数等参数沿管径的分布情况。结果发现管土装置影响下,隔水管主导频率减小;CF方向涡激振动增大,振幅沿径向呈非对称特征;涡激振动振幅更大,隔水管受流场的激励与阻尼更为剧烈。     

低质量-阻尼因子圆柱体的涡激振动预报模型

本文考查了在均匀来流中作横向振荡的圆柱体与周围流体之间的能量转移,由此建立了基于受迫振荡实验数据的弹性支撑圆柱体在均匀流中的涡激振动响应预报模型.根据此模型,分析了低质量-阻尼因子圆柱体的涡激振动响应特性.就水中圆柱体涡激振动响应特性相关的几个关键性问题进行了深入的讨论,包括响应振幅的决定因素、附加质量对锁定范围及响应频率的影响.正确理解这些问题对于深水立管涡激振动响应的有效预报至关重要.     

立管横流-顺流双向涡激振动数值预报方法

海洋立管涡激振动属于典型的流固耦合问题,涡激振动在横流及顺流方向上会同时发生并存在耦合效应.在以往的相关研究中,常忽略顺流向振动或将横流及顺流向振动响应间相互作用进行弱化.本文以顶张式立管为研究对象,基于圆柱体双自由度受迫振动试验数据,使用有限元方法和能量平衡方程建立涡激振动频域预报模型.该数值模型能够同时得到立管在横...     

海洋立管涡激振动抑制装置模型试验研究

本文提出三种不同的涡激振动抑制装置设计方案：扰流纤维抑制装置、反向耦合单组双螺旋叶片式抑制装置以及反向耦合双组双螺旋叶片式抑制装置,分别对其进行试验研究,对比分析各模型的振动特性和抑制效果。实验结果表明：所设计的涡激振动抑制装置,均能在一定程度上降低立管的涡激振动响应。从抑制效果来看,对称来流情况下的双组双螺旋装置的效果最佳,其次为单组双螺旋装置和长扰流纤维装置,而短扰流纤维装置和非对称迎流情况下双组双螺旋装置的效果较差。对称来流和非对称迎流情况下的双组双螺旋装置是相同的,因此,采用双组双螺旋装置时,应根据海域的长期洋流方向,合理地布置该装置,使其最大限度地发挥抑制作用.     

柔性拖缆涡激振动响应特性数值模拟

为分析水下拖缆在一定流速下的振动响应,使用ANSYS Workbench平台并基于双向流固耦合计算方法构建三维水下拖缆涡激振动的数值仿真框架,并对1 kn均匀流速下的拖缆涡激振动响应进行数值模拟分析。仿真得到了拖曳力、升力,振幅以及位移轨迹变化结果,为进一步研究拖缆振动响应对水下探测系统的影响提供了一定研究基础。     

两自由度不同截面形式柱体涡激振动的CFD数值模拟

文章采用计算流体力学（CFD）方法,结合SST k-棕湍流模型,对低质量比柱体进行两自由度涡激振动数值模拟,得到了柱体升力、曳力系数的时程曲线,并观察了柱体进入锁振状态的幅值变化,研究了不同截面形式柱体在外流速处于0.1-1.0 m/s范围内的振动响应。将圆柱体在不同流速下两向振动的CFD数值模拟与实验数据进行比较,得到了较为满意的结果。通过分析不同截面柱体在不同外流速下的振动幅值发现,带有抑振装置的柱体截面形式能够有效地减小涡激振动,其中,板状截面柱体抑振效果较好。     

基于一种固体区域迭代算法的圆柱涡激振动数值计算

利用Fluent平台的用户自定义程序（UDF）以及动网格模型,实现了圆柱运动方程的一种迭代求解算法,分别对层流、湍流状态下,弹性支承圆柱体在一定约化速度下的涡激响应进行了数值模拟,探讨了不同阻尼比对涡激响应的影响。结果表明：采用该迭代求解算法对弹性支承圆柱涡激振动的预测结果较为合理;随着阻尼比的逐渐增加,初始支振幅、升阻力系数时程曲线将由多频率拍振,最终变为单一频率主导的振动,且涡激振幅逐渐减小;除了质量-阻尼比联合参数m*ζ外,阻尼比ζ本身也应作为一个重要的涡激影响参数单独进行考量。     

Experimental study on dynamic response of a floating offshore wind turbine under various freak wave profiles

This study performed experimental investigation on the dynamic response of an in-place floating offshore wind turbine (FOWT) under freak wave actions. Based on the method of wave profile modulation, various freak wave profiles embedded in unidirectional Gaussian seas were generated in wave basin and the action of these waves on the FOWT was measured and analyzed, which has not been done before. The motions of FOWT were analyzed in time domain as well as time-frequency domain. The effect of freak wave parameters on FOWT motions was addressed, i.e., freak wave height, freak wave period, large crest, and deep trough. The dynamic response of FOWT was observed as a spike at the occurrence of freak wave in a conventional random wave, where the impact of freak wave can last for 17 spectral peak periods of wave. Data analysis shows that the motions of FOWT increased linearly with the freak wave height. In addition, the occurrence of freak wave induced the coupled effect on surge and pith, which was strengthen with the increase of freak wave height and wave period. Compared to a large crest, a deep trough of freak wave led to stronger motions and was supposed to be a key concern on the safety of the FOWT. The novel findings in this study provided a reference for the design of survival load on a FOWT and benchmarks for validating numerical models.     

Numerical simulation of wave-induced nonlinear motions of a two-dimensional floating body by the moving particle semi-implicit method

The moving particle semi-implicit (MPS) method was applied to compute nonlinear motions of a floating body influenced by the water on deck. To compute the motions of a rigid body, the fluid pressure at the position of each particle on the body surface was directly integrated in space and the equations of translational and rotational motions were integrated in time to determine the correct position of the rigid-body surface at each time step of the time-domain calculation. The performance of this method was validated through a comparison with measured results in an experiment that was newly conducted using a model of a box-shaped floating body with a small freeboard. Although the overall agreement was good, some discrepancies were observed for a shorter wave period, especially for the drift motion in sway. The effect of numerical resolution on the results was checked by changing the number of particles. With a higher number of particles, no obvious improvement was seen in the global body motions, but the resolution of the local free-surface profile, including the water on deck, was improved.     

通过水下储能系统的蓄电池单元的流:三个接触的似花结构的球状物(英文)

An LES simulation of flow over an accumulator unit of an underwater compressed air energy storage facility was conducted. The accumulator unit consists of three touching underwater balloons arranged in a floral configuration. The structure of the flow was examined via three dimensional iso surfaces of the Q criterion. Vortical cores were observed on the leeward surface of the balloons. The swirling tube flows generated by these vortical cores were depicted through three dimensional path lines. The flow dynamics were visualized via time series snapshots of two dimensional vorticity contours perpendicular to the flow direction; revealing the turbulent swinging motions of the aforementioned shedding-swirling tube flows. The time history of the hydrodynamic loading was presented in terms of lift and drag coefficients. Drag coefficient of each individual balloon in the floral configuration was smaller than that of a single balloon. It was found that the total drag coefficient of the floral unit of three touching balloons, i.e. summation of the drag coefficients of the balloons, is not too much larger than that of a single balloon whereas it provides three times the storage capacity. In addition to its practical significance in designing appropriate foundation and supports, the instantaneous hydrodynamic loading was used to determine the frequency of the turbulent swirling-swinging motions of the shedding vortex tubes; the Strouhal number was found to be larger than that of a single sphere at the same Reynolds number.     

桁架式Spar平台在波浪中的运动预报与实验测量

基于细长体水动力公式和刚体全非线性运动方程,以及桁架式Spar平台主体与其锚泊系统相互耦合的力和位移边界条件,建立了桁架式Spar平台在波浪中的运动响应预报模型,对其在波浪中的运动和锚泊力进行了预报。通过运动光学测量系统和拉力传感器在哈尔滨工程大学水池对Spar平台模型在波浪中的运动和锚泊力进行了测量。数值预报结果和实验测量进行了比较,结果表明二者符合得很好。     

Numerical prediction of the surf-riding threshold of a ship in stern quartering waves in the light of bifurcation theory

In order to develop design and operational criteria to be used at the International Maritime Organization (IMO), critical conditions for broaching are explored in the light of bifurcation analysis. Since surf-riding, which is a prerequisite to broaching, can be regarded as a heteroclinic bifurcation, one of global bifurcations, of a surge-sway-yaw-roll model in quartering waves, the relevant bifurcation condition is formulated with a rigorous mathematical background. Then an efficient numerical solution procedure suitable for tracing the surf-riding threshold hypersurface is presented with successful examples. This deals with all state and control variables in parallel, and excludes backward time integration and an orthogonal condition in the iteration process. The bifurcation conditions identified were compared with the results from a direct numerical simulation in the time domain. As a result, it was confirmed that the heteroclinic bifurcation provides a boundary between motions periodically overtaken by waves and nonperiodic motions such as surf-riding and broaching.     

Efficient calculation of slamming pressures on ships in irregular seas

An efficient method for calculation of the slamming pressures on ship hulls in irregular waves is presented and validated for a 290-m cruise ship. Nonlinear strip theory was used to calculate the ship–wave relative motions. The relative vertical and roll velocities for a slamming event were input to the slamming calculation program, which used a two-dimensional boundary element method (BEM) based on the generalized 2D Wagner formulation presented by Zhao et al. To improve the calculation efficiency, the method was divided into two separate steps. In the first step, the velocity potentials were calculated for unit relative velocities between the section and the water. In the next step, these precalculated velocity potentials were used together with the real relative velocities experienced in a seaway to calculate the slamming pressure and total slamming force on the section. This saved considerable computer time for slamming calculations in irregular waves, without significant loss of accuracy. The calculated slamming pressures on the bow flare of the cruise ship agreed quite well with the measured values, at least for time windows in which the calculated and experimental ship motions agreed well. A simplified method for calculation of the instantaneous peak pressure on each ship section in irregular waves is also presented. The method was used to identify slamming events to be analyzed with the more refined 2D BEM method, but comparisons with measured values indicate that the method may also be used for a quick quantitative assessment of the maximum slamming pressures.     

基于重叠网络技术的有限体积法水动力学仿真

A finite volume algorithm was established in order to investigate two-dimensional hydrodynamic problems. These include viscous free surface flow interaction with free rigid bodies in the case of large and/or relative motions. Two-phase flow with complex deformations at the interface was simulated using a fractional step-volume of fluid algorithm. In addition, body motions were captured by an overlapping mesh system. Here, flow variables are transferred using a simple fully implicit non-conservative interpolation scheme which maintains the second-order accuracy of implemented spatial discretisation. Code was developed and an appropriate set of problems investigated. Results show good potential for development of a virtual hydrodynamics laboratory.     

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.     

船舶内共振动力学行为的研究

本文考虑船舶横摇与纵摇的非线性耦合运动，采用摄动分析及数值计算方法，研究了船舶存在内共振时的非线性动力学行为。研究表明，船舶非线性运动的重要特征是响应不对称、调幅调相及饱和等现象，为揭示船舶在高海情下的倾覆机理奠定了基础。     

Grounding experiments on soft bottoms

To verify a theoretical analysis procedure for calculation of the hull girder response of ships running aground, a series of large-scale ship grounding experiments was performed on an artificial island made of engineered fill. The tests were conducted by running a condemned fishing vessel up on selected beaches of the island with velocities ranging between 2m/s and 5.5m/s. During the tests, surge, heave, and pitch accelerations were measured and also the deformations of the beach and the ship bow. Based on these accelerations, rigid body velocities and motions were determined. The forces arising from the interaction between the bow of the vessel and the seabed were determined by solving the equations of motions. This article describes the analysis of the measured results and a comparison of the results from the full-scale measurements with results from a recently developed analysis procedure for grounding on soft bottoms. Presented at the International Conference on Technology for Marine Environment Preservation (MARIENV '95), Tokyo, Japan, September 24–29, 1995.