柔性圆柱涡激振动流体力系数识别及其特性

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前,海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据,存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用,建立了柔性圆柱涡激振动流体力模型,运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数,设计并开展了拖曳水池模型实验,实验用柔性圆柱模型的质量比为1.82,长径比为195.5.通过与刚性圆柱流体力系数对比,深入分析了柔性圆柱流体力系数的特性.结果表明:柔性圆柱在一阶模态控制区,流体力系数随约化速度变化趋势与刚性圆柱大致相似;二阶模态控制区,升力系数和脉动阻力系数显著增大;附加质量系数在响应频率较低时与振动位移的相关性增强;当响应频率较低时,振动位移较大区域为能量耗散区,当响应频率较高时,振动位移较大区域为能量输入区.     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

正弦振荡来流下柔性立管涡激振动发展过程

在风浪流的作用下,海洋浮式结构物将带动悬链线立管在水中作周期性往复运动,从而在立管运动方向上产生相对振荡来流,这种振荡来流将激励立管悬垂段发生“间歇性” 的涡激振动. 在海洋工程水池中对一个4m 长的立管微段进行模型试验研究,以探索相对振荡来流作用下立管涡激振动产生的机理及其发展的物理过程. 试验通过振荡装置带动模型作正弦运动来模拟不同最大约化速度U_Rmax、不同KC（Keulegan-Carpenternumber）的相对振荡来流,利用光纤应变片测量立管涡激振动响应. 结合模态分析方法处理试验数据得到位移响应时历,继而提出相对振荡来流下柔性立管涡激振动发展的3 个阶段:建立阶段、锁定阶段以及衰减阶段. 并进一步总结了最大约化速度U_Rmax,KC 对涡激振动发展过程的影响规律. 最终获得不同最大约化速度U_Rmax下,涡激振动各发展阶段随KC 所占时间分布比例图.      

海洋热塑性增强管(RTP)涡激振动数值计算

基于Van der Pol尾流振子模型和多体系统传递矩阵法(transfer matrix method for multibody systems, MSTMM), 建立了可以快速预测海洋热塑性增强管(reinforced thermoplastic pipe, RTP)振动特性和涡激振动响应的动力学模型. 仿真结果与ANSYS软件仿真结果以及文献实验数据对比, 验证了本文模型的准确性. 研究了考虑RTP立管刚性接头, 不同顶张力, 不同来流分布等情况对RTP立管涡激振动响应的影响. 计算结果表明: 流速越大, 立管涡激振动激发出的模态越高; 立管涡激振动主要受低阶模态控制; 立管的刚性接头对立管的湿模态影响较小, 但是对较高阶模态为主所激发出的涡激振动振幅分布影响较大; 剪切流对沿立管轴向的涡激振动振幅分布影响较大, 低流速能量小所引起的涡激振动幅值较小, 但是当剪切流流速达到能激发出较高阶模态时, 相比同等流速的均匀流所引起的涡激振动振幅要大.      

基于能量传递规律的海洋立管涡激振动抑制研究

涡激振动是造成海洋立管疲劳损伤的重要因素, 抑制振动能够保障结构安全, 延长使用寿命. 多数涡激振动抑制方法基于干扰流场的方式, 但在复杂环境条件下, 仅通过干扰流场对振动的抑制效果有限. 因此, 从结构层面考虑开展了海洋立管涡激振动抑制研究. 基于能量传递的理论, 阐述了立管涡激振动过程中的能量传递规律. 振动能量以行波形式由能量输入区传播至能量耗散区, 主要在能量耗散区被消耗. 通过局部增大能量耗散区的阻尼, 增加振动能量在传播过程中的消耗, 实现涡激振动抑制. 为了求解立管涡激振动响应, 构建了尾流振子预报模型, 并根据实验结果验证了理论模型的可靠性. 基于理论计算得到的能量系数, 判定立管涡激振动的能量输入区和能量耗散区. 通过对比立管增大阻尼前后的响应, 分析了涡激振动抑制效果. 研究结果表明: 在能量输入区增大阻尼对涡激振动的抑制效果并不显著; 在能量耗散区增大阻尼使能量衰减系数达到临界值之后, 能够显著降低立管上部和底部的涡激振动位移; 当能量衰减系数超过临界值后, 继续增大耗散区阻尼对涡激振动抑制效果的提升不明显.      

深海采矿系统中悬臂式立管涡激振动分析

不同于传统的海洋立管, 深海采矿系统中的垂直提升管道可以被视为一个底部无约束的柔性悬臂式立管, 工作过程中同样面临涡激振动和柔性变形问题. 本文采用一种无网格离散涡方法和有限元耦合的准三维时域求解数值模型, 系统性地研究了不同流速下悬臂式立管的涡激振动问题. 结果表明: 悬臂式立管的横向振动模态阶数随折合速度增加而增大, 在一定折合速度范围内主导振动模态保持不变; 当主导模态转变时, 对应的横向振幅会发生突降, 但是当新的高阶模态被激发后, 立管振幅随来流速度增加而再次逐渐增大; 在相同的振动模态下, 立管底部位移均方根值随折合速度线性增加, 主导振动频率在模态转变时会出现跳跃现象; 特别地, 本文讨论了三阶主导模态下悬臂式立管的振动响应, 无约束的立管底部呈现出较大的振动能量, 且振幅的驻波特征随折合速度增加而逐渐增强; 本文比较了两端铰支立管与悬臂式立管的涡激振动响应特征, 两者在振幅和主导振动频率两方面均表现出了相同的变化趋势.      

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法.鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究.先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解.对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析.分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成.当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导.当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导.随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.     

剪切流下发生涡激振动的柔性立管阻力特性研究

利用缩尺模型试验的方法研究了线性剪切流下涡激振动发生时柔性立管的阻力特性.文中基于光纤光栅应变传感器测得的模型应变信息,采用梁复杂弯曲理论计算了立管的平均阻力,继而分析了阻力系数沿管长方向和雷诺数的分布特性以及涡激振动对阻力系数的放大效应,并提出了用于估算柔性立管发生涡激振动时阻力系数的经验公式.结果表明:涡激振动对阻力系数有放大效应,使得立管局部阻力系数高达3.2;平均阻力系数在1.0×104到1.2×105的雷诺数区间内的值为1.3~2.0,并随雷诺数的增大而减小.本文提出的经验公式可准确估算高雷诺数下涡激振动发生时柔性立管的阻力系数,此经验公式考虑了流速、涡激振动主导模态以及主导频率对阻力系数的影响.      

海洋柔性结构涡激振动的流固耦合机理和响应

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

Numerical scheme for riser motion calculation during 3-D VIV simulation

This paper presents a numerical scheme for riser motion calculation and its application to riser VIV simulations. The discretisation of the governing differential equation is studied first. The top tensioned risers are simplified as tensioned beams. A centered space and forward time finite difference scheme is derived from the governing equations of motion. Then an implicit method is adopted for better numerical stability. The method meets von Neumann criteria and is shown to be unconditionally stable. The discretized linear algebraic equations are solved using a LU decomposition method. This approach is then applied to a series of benchmark cases with known solutions. The comparisons show good agreement. Finally the method is applied to practical riser VIV simulations. The studied cases cover a wide range of riser VIV problems, i.e. different riser outer diameter, length, tensioning conditions, and current profiles. Reasonable agreement is obtained between the numerical simulations and experimental data on riser motions and cross-flow VIV a/D. These validations and comparisons confirm that the present numerical scheme for riser motion calculation is valid and effective for long riser VIV simulation.     

Monitoring VIV fatigue damage on marine risers

Long flexible cylinders (e.g., risers, tendons and mooring lines) exposed to the marine environment encounter ocean currents leading to vortex-induced vibration (VIV). These oscillations, often driven at high frequencies over extended periods of time, may result in structural failure of the member due to fatigue damage accumulation. Recent developments in instrumentation and installation of data acquisition systems on board marine risers have made accurate measurement of riser responses possible. This paper aims at using the data from these data acquisition devices (typically strain gages and accelerometers) in order to understand the evolution of the riser VIV, with the final aim of estimating the fatigue damage. For this purpose we employ systematic techniques to reconstruct riser VIV response using the data from the available sensors. The reconstructed riser response allows estimation of the dynamic axial stresses due to bending and consequently the estimates of the fatigue damage along the entire riser. The above methods can take into account the fatigue damage arising from complicated riser motions involving the presence of traveling waves even with the use of very few sensors. An alternate approach using a Van der Pol wake oscillator model is also explored to obtain fatigue life estimates caused by riser VIV.     

Optimal lift force coefficient databases from riser experiments

Significant past effort has gone into understanding the complicated flow–structure interaction problem of vortex-induced vibration (VIV) of long flexible cylindrical structures (e.g., risers, mooring lines, tendons, conductors) in the ocean environment. However, major challenges persist with regard to riser VIV modeling and response prediction. The existing prediction schemes are based on a number of hypotheses, experimental facts and data like strip theory, energy balance, correlation length and, most importantly, the use of lift force coefficient databases. Recent advances in observing the VIV motions on experimental risers with high confidence shows that some of these assumptions may not be valid. One important source of the discrepancies between theoretical estimates and experimental observations arise from the use of experimentally obtained lift coefficient databases. These databases were obtained under the laboratory conditions of limited Reynolds number, and under the assumption that the cross-flow motions are not influenced by restraining the in-line motions. In this paper we develop a method to improve the modeling capability of riser VIV by extracting empirical lift coefficient databases from field riser VIV measurements. The existing laboratory-based lift coefficient databases are represented in a flexible parameterized form using a set of carefully chosen parameters. Extraction of the lift coefficient parameters is posed as an optimization problem, where the objective is to minimize the error between the prediction using a theoretical model and the experimental data. Application of the method to data from the Norwegian Deepwater Programme experiments shows that the new optimal databases significantly reduce the error in estimating the riser VIV cross-flow response.     

Effect of internal flow on vortex-induced vibration of risers

Although there are many studies dedicated to the problem of vortex-induced vibration (VIV) of marine risers, VIV experiments with internally flowing fluid have not been carried out before. In order to investigate this area, the present experiment with an internally flowing fluid and external current was designed. The riser was towed in the water flume with varying internal flow speeds. The tests in still water and in a current were conducted successfully. Various measurements were obtained including the frequency responses and the time-domain tracing of in-line and cross-flow responses. The experimental results exhibit several valuable features. First, with an increase in internal flow speed, the response amplitude increases while the vibration frequency decreases. Secondly, internally flowing fluid lessens the correlation of the vibration between different sections. In addition, by plotting both in-line strain and cross-flow strain simultaneously, a figure-of-eight for bending strain is also observed, and the trajectories in different cycles are more concordant with the increase of internal flow speed.     

Experimental study on dynamic response of smooth and straked risers in tandem arrangement coupling interference effect

An experiment was conducted in a combined wave–current water flume on two tandem risers subjected to uniform flow. The riser model has an effective length of 2.0 m. The aspect ratio is 111.11. The upstream riser is smooth and the downstream riser fitted with three-strand helical strakes with pitch 17.5D and height 0.25D. By varying the external flow velocities and spacing ratios, through comparisons with the dynamic response of isolated smooth and isolated straked riser, the paper observe how interference effect impacts the dynamic characteristics and dynamic response of two risers in tandem arrangement, reveal how the suppression efficiency of the three-strand helical strakes responds to spacing ratio and external flow velocity, and explore the wake excitation effect of inter-riser fluids on the downstream riser and their dynamic feedback to the upstream riser. The results show that the dominant frequency of the upstream smooth riser is sensitivity to the change of the spacing ratio is low, and the displacement response is offset or enhanced in different degree due to the difference of the interference efficiency. The downstream straked riser dominates frequency and displacement higher than the isolated straked riser. The wake vortex of the upstream smooth riser acts on the downstream riser, occupying a dominant position in the vibration of the downstream riser. It degrades the vibration suppression efficiency of the three-strand helical strakes: the suppression efficiency is the highest at spacing ratio of 8, being a merely 70.57%. With the increase of the spacing ratio, the CF displacement of the riser gradually decreases, and the IL displacement gradually increases. The interference efficiency partition and suppression efficiency at different spacing ratios reflects the dynamic feedback of the upstream smooth riser is much smaller than the interference effect of the downstream suppression riser.     

A systematic approach to riser VIV response reconstruction

Vortex-induced vibration (VIV) of long flexible cylindrical structures (e.g. risers, pipelines, tendons, mooring lines) enduring ocean currents is ubiquitous in the offshore industry. Though significant effort has gone into understanding this complicated fluid–structure interaction problem, major challenges remain in modeling and predicting the response of such structures (for example a riser). The work presented in this paper provides a systematic approach to estimate and analyze the vortex-induced motions of a marine riser. A systematic framework is developed, which allows reconstruction of the riser motion from a limited number of sensors placed along its length. A full reconstruction criterion is developed, which classifies when the measurements from the sensors contain all information pertinent to riser VIV response, and when they do not, in which case additional, analytical methods must be employed. Reconstruction methods for both scenarios are developed and applied to experimental data. Finally, a systematic study on the error during the reconstruction is also undertaken. The methods developed in this paper can be applied to: improve understanding of the vortex shedding mechanisms, including the presence of traveling waves and higher-harmonic forces; develop tools for in-situ estimation of fatigue damage on marine risers; and estimate the vortex-induced forces on marine risers.     

Fluid-solid coupling and dynamic response of vortex-induced vibration of slender ocean cylinders

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.