深水浮式平台垂荡运动与水下柔性立管涡激振动的动力耦合

基于有限元数值模拟,进行了"平台垂荡-顶张力立管涡激振动"整体系统的动响应数值模拟。动响应模型考虑了立管尾迹流场的水动力与结构动力的耦合和垂荡引起的立管结构刚度的时变特性;分析了平台垂荡运动的频率、模态阶数等因素对水下顶张力立管涡激振动的影响。数值结果表明:与不考虑平台运动相比,立管的动响应位移会增大;立管响应幅值随着模态阶数的降低而增大;在响应过程中,尤其对于低阶模态,会出现响应的模态转换现象。鉴于在平台垂荡和涡激振动的共同作用下,立管的动响应会大于涡激振动、参数激励分别单独作用的响应,建议在立管实际工程设计中应该考虑平台运动和涡激振动耦合激励作用下的结构动响应。     

张力腿平台涡激运动模型试验与数值分析

为揭示张力腿平台涡激运动响应规律,采用数值模拟与模型试验相结合的方法研究了全水深系泊张力腿平台涡激运动响应。根据张力腿平台主尺度参数,按照几何相似制作了水池试验模型,在满足运动相似和动力相似的条件下开展了均匀流、剖面流模型试验,并将模型试验结果与数值模拟结果进行了对比,验证了数值模拟与模型试验结果的一致性。分析结果表明在均匀流作用下张力腿平台涡激运动的锁定区间在5.5<U_r<8.5,来流角对涡激运动影响较大。剖面流作用下平台涡激运动规律与均匀流作用基本一致,但剖面流造流引起的能量分散,使平台在XY平面的运动轨迹变得不规律,系泊系统提供的回复力对涡激运动幅值有抑制作用,来流角和流速对张力腿或立管模态影响明显。论文得到结论对于张力腿平台的工程设计有一定借鉴意义。     

波、流联合作用下海洋平台张力腿的涡激非线性振动

本文研究深水张力腿平台(以下简称TLP)的张力腿在波、流联合作用下涡街激起的振动.在所建立的数学模型中既考虑平方流体阻尼项,同时又计及中心激振影响,文中假定微幅波和线性分布的不均匀流是沿同一方向的.事实上,此方法可扩展到波、流呈任意角度的情况下,沿流向或垂直流向动力响应的研究.采用了多项伽僚金法和综合数值解法,以Conoco's Hutton张力腿平台为例,计算了当涡激频率与张力某固有频率相近发生谐振时,张力腿涡激非线性动力响应,给出一些重要结果.     

多频参激—涡激联合作用下顶张力立管分岔分析

复杂海洋环境中,顶张力立管顶部平台受到不同分量叠加的波浪作用,导致其轴向张力与多个频率相关,系统产生多频参激与涡激的联合振动,其动力学行为变得更为复杂,尤其是分岔和混沌等非线性振动特性有待进一步分析。首先,基于欧拉—伯努利梁理论,引入范德波尔尾流振子,建立和推导受多频参激—涡激联合作用下立管的动力学模型及其运动微分方程;接着,利用伽辽金法得到离散后的常微分方程组,采用多尺度法得到系统共振响应的调谐方程;最后,通过数值算例探究系统在多源和多频载荷联合作用下的共振响应。结果表明：对于多频参数激励,其中一个参数激励幅值在立管振动中起主导作用,而另一个会引起分岔的位置、数量和性质发生改变;随质量—阻尼参数增大,系统共振响应幅值整体先增大后减小,各分岔点出现的位置整体先向泄涡频率较大处移动,随后向泄涡频率较小处移动;质量—阻尼参数的变化还可能诱使系统发生倍周期分岔和混沌现象,危害结构安全。此外,采用直接数值积分求解微分方程组,验证近似分析结果,两者吻合较好。     

传统Spar平台和型组合共振运动响应研究

研究传统Spar平台在长周期规则波浪下的和型组合共振现象.考虑垂荡与纵摇运动的相互影响以及瞬时波面对于纵摇运动的参数激励作用,建立平台两自由度耦合运动方程.采用多尺度法导出了当波浪频率接近垂荡与纵摇固有频率之和时方程的一阶摄动解.针对一座实际Spar平台,通过求解幅值分岔方程,得到平台在不同波高下的运动响应并进行了数值模拟.计算结果发现,当波浪幅值超过某一极限值时,平台出现较大幅值的垂荡与纵摇亚谐共振运动.研究结果揭示了Spar平台出现和型组合共振的条件和组合共振时的动力特性.     

Truss Spar平台垂荡响应频域分析

利用已有水动力试验和数值模拟结果及Morison方程,简化估计了Truss Spar平台垂荡方向水动力贡献,快速有效地预报了平台极端海况下的垂荡响应.在频域内分别采用数值迭代计算和粘滞阻尼线性化方法得到了平台在不同数量垂荡板配置下的垂荡响应幅值算子.在墨西哥湾、西非和中国南海等3种不同地区极端海况以及墨西哥湾工作状态海况下计算了平台的响应极值,分析发现平台波频垂荡响应对入射波浪周期非常敏感.采用JONSWAP波能谱得到了Truss Spar平台垂荡响应谱.通过比较和参数分析,验证了平台在设计上的合理性与实际海况下良好的垂荡运动性能.     

悬跨海底电缆涡激振动特性及疲劳损伤分析

悬跨海底电缆作为细长柔性结构,在静力平衡状态下具有一定的垂度,在水流作用下的涡激振动特性与海底管线和海洋立管等结构也有很大的不同,其振动模态受垂跨比影响很大。通过物理模型试验开展了不同垂跨比下悬跨海缆的涡激振动和疲劳损伤特性研究。试验模型按照水弹性相似准则设计,试验中测量了不同流速下海缆模型产生涡激振动时的应变历时数据,采用模态分析法获得了模型涡激振动时的振动模态和振幅。分析了不同流速下海缆模型的振动模态、应变和疲劳损伤的变化和分布特征。试验结果表明：垂跨比显著影响了海缆的涡激振动模态和应变幅值大小。在本试验流速范围内,对一定长度的悬跨海缆模型,当垂跨比较大时,随着流速的增大,模型涡激振动的主响应振动依次出现反对称1阶和对称1阶模态;当垂跨比较小时,模型涡激振动的主响应模态依次出现反对称1阶和对称2阶模态。当涡激振动主响应模态为反对称1阶时,疲劳损伤最大值达到0. 1~0.7。     

畸形波作用下四立柱张力腿平台动力响应研究

介绍了畸形波生成的方法并对各种线性合成畸形波的效率进行了对比。选用改进的随机波加瞬态波的方法模拟强非线性波———"新年波"的波形。采用微元法对张力腿平台的立柱和浮筒进行离散,编制程序对张力腿平台在强非线性波作用下的耦合动力响应进行了数值模拟,重点比较分析了张力腿平台在随机波及畸形波中所受波浪力、平台动力响应、系泊系统张力特性及浪向角对平台运动的影响。研究表明:在"新年波"作用下,0°浪向时,平台在x和z方向所受波浪力较随机波增长了约1/4,纵荡、垂荡及纵摇响应值分别增加了33%,38%和12%,张力腿张力幅值增加约20%,位于平台四角的张力腿张力有所差异,相邻两根张力腿的张力差别不大,浪向角对平台运动响应、张力腿张力的影响畸形波大于随机波的作用。为今后考虑强非线性波浪影响时张力腿平台的设计提供了借鉴和参考。     

立柱直径比对四立柱平台涡激运动性能影响的数值研究

与单立柱平台(Spar)类似,四立柱平台在洋流作用下也容易发生涡激运动,对锚泊和立管系统的疲劳寿命产生不利影响。一种非对称半潜式平台具有两种不同尺寸的立柱,立柱直径比对非对称平台涡激运动性能的影响尚需深入研究。采用计算流体力学(CFD)方法,对直径比分别为1.0、1.5和2.0的立柱组合在两种不同角度水流作用下的运动响应进行了二维数值模拟,分析发现:立柱直径比不为1时四立柱平台也具有明显的涡激运动,其横荡运动具有显著的锁定现象;在广义折合速度定义下,与相同直径的四立柱平台相比,立柱直径比不为1时其横荡运动锁定区间变宽且滞后,横荡响应幅值的峰值具有明显差异;在所研究的直径比范围内,大直径立柱在上游时涡激运动响应比其在下游时更显著,且这种差异随着直径比变大而变大。     

波流作用下悬浮隧道运动响应纵向截断模型试验研究

利用波流水槽构建悬浮隧道纵向截断模型试验,分析了水流及波浪作用对悬浮隧道运动响应的影响,获得了管体结构的垂向及横向运动响应特性。研究结果表明:水流流速较小时,模型振幅随约化流速Vr增大而增大,而在约化流速为5.8附近,管体振幅达到最大并在这一临界值之后形成较为稳定的振幅值,产生明显的结构共振和涡激振动“锁定”现象。振幅稳定后垂向幅值约为管径的2倍,横向幅值与管径相当。对于波浪单独作用,在波浪较弱时,横向幅值与垂向几乎相等,随着特征KC_n数的增大,响应振幅先增后减。在特征KC_n数小于0.51范围内,垂向振幅大于横向振幅,其余范围则相反。垂向振幅和横向振幅分别在特征KC_n数为0.42和0.55时达到峰值,且两个方向的最大振幅值均与管体的管径相当。相对响应频率随着特征KC_n数的增大而减少,且垂向大于横向。     

Nonlinear dynamics of a fluid-conveying pipe under the combined action of cross-flow and top-end excitations

We report a theoretical investigation of an elastic and slender fluid-conveying pipe with a top-end excitation subjected to uniform cross flows. Considering the mean drag force and the time varying vortex-induced lift force which is modeled using a nonlinear van der Pol oscillator, the nonlinear partial differential equations of the motion of coupled fluid-structure system are constructed and simplified to a reduced-order model through the Galerkin-type discretization. By virtue of quasi-static displacement conditions, the characteristics of vortex-induced vibration of the pipe are evaluated for the first two lock-in modes. The results show that the top-end excitation can increase the vibration amplitude of the pipe when the cross-flow speed is out of the lock-in regions. When the cross-flow speed is within the lock-in region, however, the top-end oscillation causes a transition between quasi-periodic and periodic in the responses of the pipe, significantly reducing or increasing the vibration amplitudes depending on the excitation acceleration and frequency. This finding has an important guidance in suppressing vortex-induced vibrations by balancing the internal fluid velocity and the top-end excitation.     

Stability analysis of TLP tethers

Compliant offshore TLPs are essentially meant for deep oil/gas exploration and are usually constructed on the seashore and then towed down to the particular location for anchorage. They are connected to the sea bed by means of pretension cables. The increased use of TLPs in deep waters and necessity of reduction of usually high value of pretension make the effect of variable tension in the tether dynamics more significant. This work presents the dynamic analysis of tethers and TLPs considering the linearly varying tension along the tether length. The modal analysis considers a linear cable equation for tether modeling subjected to tension which varies along its length. A Mathieu stability analysis is then performed for TLPs of different shapes and different water depth vis-à-vis of 527.8, 872, and 1200 m respectively to obtain the amplitudes of tether vibrations. The unstable modes of vibration are also verified. The resultant modal forms for the tether's dynamic model are then obtained in form of Bessel's function. From the numerical studies conducted it is seen that increased tether tension not only leads to a stable platform but also improves the stability due to increased hydrodynamic loading contributing to added mass. From the studies conducted it is also seen that the triangular configuration TLPs with increased initial pretension are more stable compared to four leg TLP in the first mode of vibration.     

Nonlinear Dynamic Responses of Tension Leg Platform with Slack-Taut Tether

The slack-taut state of tether is a particular adverse circumstance,which may influence the normal operation state of tension leg platform (TLP).The dynamic responses of TLP with slack-taut tenor are studied with consideration of several nonlinear factors introduced by large amphliude motions.The time histories of stresses of tethers of a typical TLP in slacktaut state are given.In addition,the sensitivities of slack to stiffness and mass are investigated by varying the stiffness of tether and mass of TLP.It is found that slack is sensitive to the mass of TLP.The critical curved surfaces (over which indicates the slack) for the increase of mass are obtained.     

Influence of wave approach angle on TLP's response

The current study focuses on the response analysis of triangular tension leg platform (TLP) for different wave approach angles varying from 0° through 90° and its influence on the coupled dynamic response of triangular TLPs. Hydrodynamic loading is modeled using Stokes fifth-order nonlinear wave theory along with various other nonlinearities arising caused by change in tether tension and change in buoyancy caused by set down effect. Low frequency surge oscillations and high frequency tension oscillations of tethers are ignored in the analysis. Results show that wave approach angle influences the coupled dynamic response of triangular TLP in all degrees of freedom except heave. Response in roll and sway degrees of freedom are activated which otherwise are not present in TLP's response to unidirectional waves. Pitch and roll responses are highly stochastic in nature indicating high degree of randomness. Variation in surge, sway and heave responses are nonlinear and are not proportional to change in wave height for the same period.     

水下悬浮隧道锚索稳定性分析

为研究水下悬浮隧道锚索的稳定性,建立锚索在涡街激励作用下振动的数学方程,采用伽辽金法对其进行化简,锚索振动的稳定性通过Lyapunov指数法进行判断,分析了锚索动静张力比的大小、悬浮隧道系统的阻尼比和锚索参数激励频率对锚索振动稳定性的影响。计算结果表明,锚索振动失稳的范围取决于锚索的动静张力比、锚索的阻尼比和频率比(锚索参数激励频率与锚索1阶模态固有频率的比值);当频率比为1和2左右时,随着锚索动静张力比的增大和阻尼比的减小,锚索逐渐从稳定状态转变为不稳定状态,且锚索的不稳定区域逐渐增大。     

Numerical Simulation of Self-Excited and Forced Vibrating Circular Cylinders in Current

Numerical simulations of a low-mass-damping circular cylinder which can oscillate freely at transverse and streamwise directions are presented in this work. The Navier-Stokes equations are solved with finite volume method, and large eddy simulation of vortex is also performed in the calculation. In order to implement dynamic mesh, overlapping grids are generated to lessen the computation for mesh field itself. Self-excited vibrations are firstly calculated to obtain the average amplitudes and frequencies of the target circular cylinder in the current flow situation, and then forced oscillations are implemented with parameters obtained in vortex-induced vibrations previously. With slight amplitude modulation, time series of displacements in vortex-induced vibrations are essentially harmonic. Regarding the fluid force, which are larger in forced oscillations than those in corresponding self-excited cases because the fluid subtracts energy from the forced cylinders. The phase angles between forces and displacements are 0o and 180o for self-excited case and forced case respectively. In vortex-induced vibrations, the interactions between fluid and structure produce some weakly energetic vortices which induce the modulations of amplitude and frequency.     

Numerical Simulation of Self-Excited and Forced Vibration of Circular Cylinders in Current

Numerical simulations of a low-mass-damping circular cylinder which can oscillate freely at transverse and streamwise directions are presented in this work.The Navier-Stokes equations are solved with finite volume method,and large eddy simulation of vortex is also performed in the calculation.In order to implement dynamic mesh,overlapping grids are generated to lessen the computation for mesh field itself.Self-excited vibrations are firstly calculated to obtain the average amplitudes and frequencies of the target circular cylinder in the current flow situation,and then forced oscillations are implemented with parameters obtained in vortex-induced vibrations previously.With slight amplitude modulation,time series of displacements in vortex-induced vibrations are essentially harmonic.Regarding the fluid force,which are larger in forced oscillations than those in corresponding self-excited cases because the fluid subtracts energy from the forced cylinders.The phase angles between forces and displacements are 0° and 180° for self-excited case and forced case respectively.In vortex-induced vibrations,the interactions between fluid and structure produce some weakly energetic vortices which induce the modulations of amplitude and frequency.     

Numerical Investigation on Vortex-Induced Rotations of A Triangular Cylinder Using An Immersed Boundary Method

A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).