双向剪切流作用下柔性立管平均阻力特性研究
Mean Drag Force of Flexible Riser Under Bidirectionally Sheared Flow
通讯作者: 付世晓,教授,博士生导师,电话(Tel.):021-34207053;E-mail:shixiao.fu@sjtu.edu.cn.
责任编辑: 王一凡
收稿日期: 2023-04-7 修回日期: 2023-05-18 接受日期: 2023-05-29
基金资助: |
|
Received: 2023-04-7 Revised: 2023-05-18 Accepted: 2023-05-29
作者简介 About authors
付雪鹏(1996—),博士生,从事涡激振动研究.
中国南海存在由内孤立波诱发形成的特殊双向剪切流场.开展了双向剪切流场下的涡激振动状态柔性立管的阻力特性实验研究,基于模态分析法以及张力梁弯曲理论处理实验测量应变时历得到了立管模型的平均阻力,分析了双向剪切流下的柔性立管初始位移、平均阻力以及平均阻力系数特性,确证了双向剪切流具有与线性剪切流相当的阻力系数放大效果.同时发现了双向剪切流下柔性立管独特的中部剪力极值特性,提出剪力系数并拟合得到了双向剪切流下柔性立管中部剪力极值经验公式,该公式将为南海立管设计提供强度校核外载荷输入.
关键词:
There exists a special bidirectionally sheared flow field due to solitons in the South China Sea. An experimental study of the mean drag force of a flexible riser undergoing vortex-induced vibration under bidirectionally sheared flow is conducted. Modal analysis and the beam theory are applied to processing the strain signal measured in the experiment to obtain the mean drag force. The initial displacement, mean drag force, and mean drag coefficient are investigated, and the amplification of the mean drag coefficient is confirmed under bidirectionally sheared flow with comparable amplification to linear shear flow. Meanwhile, a unique phenomenon that the shear force reaches an extreme value in the center of the riser under bidirectionally sheared flow is found in the experiment. The shear force coefficient is proposed and fitted to obtain an empirical formula for the extreme value of the central shear force of the flexible riser under bidirectionally sheared flow, which will provide external load input for the design of the riser in the South China Sea.
Keywords:
本文引用格式
付雪鹏, 付世晓, 张萌萌, 许玉旺, 任浩杰, 孙童晓.
FU Xuepeng, FU Shixiao, ZHANG Mengmeng, XU Yuwang, REN Haojie, SUN Tongxiao.
在背景洋流作用下,柔性立管会在顺流向(in-line,IL)以及横流向(cross flow,CF)发生涡激振动现象(vortex-induced vibration, VIV)[1⇓⇓⇓-5].对于横流向,立管所受水动力通常分解为均值为定值的平均阻力以及均值为零的脉动阻力.顺流向的涡激振动位移由平均阻力导致的初始位移以及脉动阻力导致的振动位移耦合组成.在立管设计中,平均阻力通常基于莫里森公式计算,其中最主要的参数为阻力系数.不同于刚性圆柱绕流时所受平均阻力,涡激振动状态下柔性立管通常伴随着明显的阻力放大系数[6-7],因此确定柔性立管涡激振动状态下的阻力系数是开展柔性立管的强度设计与校核工作的关键基础[8].许多研究学者已经针对不同流场下的柔性立管涡激振动特性开展了研究.Tognarelli等[9]开展了均匀流及剪切流下的柔性立管涡激振动响应特性实验,获取了不同雷诺数工况下顺流和横流向的涡激振动响应特征.Trim等[10]基于挪威深水计划开展了长细比达1 400的柔性立管涡激振动响应特性实验.针对柔性立管的阻力特性研究,Song等[11]开展了剪切流下柔性立管涡激振动响应特性试验,分析了阻力系数沿管长方向和雷诺数的分布特性以及涡激振动对阻力系数的放大效应,同时提出了当估算柔性立管发生涡激振动时阻力系数的经验公式[12-13].然而,在我国南海存在一种由内孤立波诱发的双向剪切流场,这种特殊的流场以流速分离点(O)为界呈现流速方向相反、幅值相等的流速分布特性.
当前针对双向剪切流场下的柔性立管阻力特性研究仍为空白,严重影响了南海油田的相关立管设计工作.本文设计了双向剪切流实验室模拟平台,开展了双向剪切流作用涡激振动状态时柔性立管阻力特性,基于光纤光栅应变传感器测得柔性立管模型顺流向应变时历,采用张力作用下梁弯曲理论计算柔性立管模型的平均阻力,并根据莫里森公式计算实验模型在不同流速工况下的阻力系数,研究了不同雷诺数下阻力系数的分布特性以及涡激振动状态对阻力系数的影响.实验中首次发现了双向剪切流流场诱发的剪力极值现象,对该现象进行了分析,总结提出了双向剪切流场下特有的立管模型剪力系数.这将为南海立管设计提供重要的参考依据.
1 实验装置、实验模型和实验方案
图1
图2
表1 实验立管模型物理参数
Tab.1
参数 | 取值 |
---|---|
模型长度,L/m | 7.64 |
外径,D/mm | 28.41 |
单位长度质量, | 1.24 |
弯曲刚度,EI/(N·m2) | 58.60 |
拉伸刚度,EA/kN | 940 |
预张力,Tpre/N | 980 |
阻尼比,ζ/% | 2.58 |
图3
为了测量立管模型的顺流向涡激振动应变响应,在实验模型的两个方向布置4组光纤应变片,于顺流向分别布置14个应变测点,如图4所示.实验立管应变数据采样频率为250 Hz.
图4
约化速度定义如下:
式中:f1为实验测得立管一阶固有频率.实验流速工况以边缘最大流速定义,实验流速工况由0.30~1.39 m/s共58个工况,对应的约化速度范围为6.59~30.79.
2 涡激振动数据处理方法
2.1 应变信号预处理
涡激振动状态下柔性立管顺流向应变信号通常可以分解为
式中:εCF-T(z, t)是张力导致的拉伸应变;εmb(z)是顺流向的初始弯曲应变;εviv-IL(z, t)为顺流向的涡激振动应变.涡激振动以及初始弯曲应变有如下关系:
2.2 模态分析法
以顺流向为例,对模态分析原理进行说明.模态分析法基于线性模态叠加,可以将模型表面测得的应变信号结合结构的模态振型计算得到测点的位移.一般地,模型表面的位移可以表示为
式中:pi(t)表示第i阶模态的位移权重;φi(z)表示模型的第i阶位移振型.
根据位移与曲率的关系,可以得到曲率κ(z,t)的表达式为
本实验中立管模型两端边界条件为简支端,其第i阶位移振型可以用正弦三角函数表示为
则应变可以表示为
式中:R为立管模型半径.
2.3 平均阻力逆向识别方法
根据直梁横向弯曲振动理论,带张力的梁弯曲的运动控制方程为
式中:w(z, t)为模型中和轴上的位移;T为顶端瞬时张力;c为阻尼系数;f(z, t)为模型所受到的水动力分布载荷.
基于模态分析法,可以获得位移模态权重pi(t),利用模态函数φi(x)可以得到位移对空间变量z的各阶导数:
进一步应用中心差分法,可获得模型各节点处的速度和加速度:
代入式(7)即可得到沿立管模型轴线方向上的水动力载荷分布.实验立管模型在顺流向所受平均拖曳力[18]为
式中:
对于本实验,水动力外径为D的柔性立管模型,其在背景流作用下沿管长所受拖曳力可表示为
式中:CD,m(z)为节点z处的平均阻力系数;ρ为流体密度;U为节点背景流速.节点处的平均阻力系数可以表示为
3 结果分析
3.1 顺流向初始位移
图5
图5
0.52 m/s流速工况下初始位移和背景流速沿管长分布情况
Fig.5
Distribution of initial displacement and flow velocity at Umax = 0.52 m/s
图6
图6
0.90 m/s流速工况下初始位移和背景流速沿管长分布情况
Fig.6
Distribution of initial displacement and flow velocity at Umax = 0.90 m/s
图7
图7
1.09 m/s流速工况下初始位移和背景流速沿管长分布情况
Fig.7
Distribution of initial displacement and flow velocity at Umax = 1.09 m/s
图8展示了不同流速工况下沿管长最大拖曳无量纲位移.可以看出最大拖曳无量纲位移随着约化速度增加而增加,实验工况下最大拖曳位移达到了0.55D.
图8
3.2 平均阻力系数特性
图9
图10
图10
0.52 m/s流速工况下平均阻力系数沿管长分布情况
Fig.10
Distribution of mean drag coefficient at Umax = 0.52 m/s
图11
图11
1.09 m/s流速工况下平均阻力沿管长分布情况
Fig.11
Distribution of mean drag at Umax=1.09 m/s
图12
图12
1.09 m/s流速工况下平均阻力系数沿管长分布情况
Fig.12
Distribution of mean drag coefficient at Umax = 1.09 m/s
图13展示了平均阻力系数随着雷诺数变化的趋势.图中:CD,m,a为总平均阻力系数,为立管模型沿管长平均阻力系数的平均值(除0.4<z/L<0.6).实验雷诺数定义为
图13
3.3 剪力分布特性
图14
图14
0.52 m/s流速工况下剪力沿管长分布情况
Fig.14
Distribution of shear force at Umax=0.52 m/s
图15
图15
0.90 m/s流速工况下剪力沿管长分布情况
Fig.15
Distribution of shear force at Umax=0.90 m/s
图16
图16
1.09 m/s流速工况下剪力沿管长分布情况
Fig.16
Distribution of shear force at Umax=1.09 m/s
为了定量分析该剪力极值,本文引入剪力系数CSF衡量双向剪切流诱导中部剪力极值的幅值,剪力系数与中部剪力极值的关系为
式中:FSF,max为中部剪力极值.
图17展示了剪力系数的拟合结果.采用最小二乘法进行系数拟合,得到剪力系数值为0.061 88.最终获得双向剪切流下中部剪力极值经验公式为
该剪力系数将为南海立管结构强度设计提供外载荷参数输入.
图17
4 结论
针对南海独有的双向剪切流场,本文开展了双向剪切流下柔性立管模型平均阻力特性实验研究,对双向剪切流下立管模型的初始位移、平均阻力、平均阻力系数以及剪力分布特性进行了研究,主要结论有:
(1) 不同于以往的均匀流以及剪切流等流场,双向剪切流下柔性立管会产生反对称的初始位移以及平均阻力分布.平均阻力极值点相较于初始位移极值点更偏向立管端部.
(2) 确证了双向剪切流下涡激振动状态柔性立管仍会产生阻力放大效应,实验结果表明平均阻力系数将会放大至1.47.这与剪切流下涡激振动柔性立管阻力放大效应相当.
(3) 首次发现了双向剪切流致使的柔性立管中部剪力极值现象,引入了剪力系数CSF并提出了中部剪力极值的经验公式FSF,max=0.5CSFρDL
参考文献
Vortex-induced vibrations
[J].
A critical review of the intrinsic nature of vortex-induced vibrations
[J].
内输高温高压流体海底悬跨管道的非线性涡激振动响应分析
[J].
Nonlinear dynamic response of suspended span pipe conveying high temperature and high pressure flow
[J].
State-of-the-art review of vortex-induced motions of floating offshore wind turbine structures
[J].
不同表面粗糙度下圆柱体涡激振动响应特性数值研究
[J].
DOI:10.16183/j.cnki.jsjtu.2018.04.006
[本文引用: 1]
使用数值模拟方法,研究了不同粗糙度下圆柱体的涡激振动响应特性,对圆柱体的涡激振动位移响应幅值、响应频率、涡激力、漩涡泄放模式以及涡激振动位移响应与涡激力之间的相位角等参数进行了分析.结果表明,圆柱体涡激振动响应幅值随着粗糙度的上升呈下降趋势.当圆柱体表面光滑或表面粗糙度较小时,涡激振动响应可分为初始分支、上分支以及低分支.在初始分支以及上分支区间,圆柱体尾部漩涡泄放模式呈2S模式;在低分支区间,圆柱体尾部漩涡泄放模式呈2P模式.当圆柱体表面粗糙度较大时,涡激振动响应仅存在初始分支和低分支.初始分支区间,尾部漩涡泄放模式呈2S模式;低分支区间,尾部漩涡泄放模式呈2P模式.
Numerical study of the vortex induced vibrations of a circular cylinder with different degrees of surface roughness
[J].
Novel geometries of serrated helical strakes to suppress vortex-induced vibrations and reduce drag
[J].
Vortex-induced vibrations of flexible cylinders predicted by wake oscillator model with random components of mean drag coefficient and lift coefficient
[J].
Drag coefficients of double unequal-diameter flexible cylinders in tandem undergoing vortex/wake-induced vibrations
[J].
VIV response of a long flexible cylinder in uniform and linearly sheared currents
[C]
Experimental investigation of vortex-induced vibration of long marine risers
[J].
An investigation into the hydrodynamics of a flexible riser undergoing vortex-induced vibration
[J].
Distribution of drag force coefficient along a flexible riser undergoing VIV in sheared flow
[J].
Hydrodynamics of flexible pipe with staggered buoyancy elements undergoing vortex-induced vibrations
[J].
南海流花海域内孤立波特征分析
[J].
Characteristics of internal solitons at Liuhua Area in the South China Sea
[J].
Experimental investigation of vortex-induced vibration of a flexible pipe in bidirectionally sheared flow
[J].
Comprehensive riser VIV model tests in uniform and sheared flow
[C]
Modal analysis of measurements from a large-scale VIV model test of a riser in linearly sheared flow
[J].
Hydrodynamic forces and coefficients on flexible risers undergoing vortex-induced vibrations in uniform flow
[J].
Laboratory measurements of vortex-induced vibrations of a vertical tension riser in a stepped current
[J].
剪切流下发生涡激振动的柔性立管阻力特性研究
[J].
DOI:10.6052/0459-1879-15-309
[本文引用: 1]
利用缩尺模型试验的方法研究了线性剪切流下涡激振动发生时柔性立管的阻力特性.文中基于光纤光栅应变传感器测得的模型应变信息,采用梁复杂弯曲理论计算了立管的平均阻力,继而分析了阻力系数沿管长方向和雷诺数的分布特性以及涡激振动对阻力系数的放大效应,并提出了用于估算柔性立管发生涡激振动时阻力系数的经验公式.结果表明:涡激振动对阻力系数有放大效应,使得立管局部阻力系数高达3.2;平均阻力系数在1.0×10<sup>4</sup>到1.2×10<sup>5</sup>的雷诺数区间内的值为1.3~2.0,并随雷诺数的增大而减小.本文提出的经验公式可准确估算高雷诺数下涡激振动发生时柔性立管的阻力系数,此经验公式考虑了流速、涡激振动主导模态以及主导频率对阻力系数的影响.
Investigation of drag forces for flexible risers undergoing vortex-induced vibration in sheared flow
[J].
DOI:10.6052/0459-1879-15-309
[本文引用: 1]
<p>The drag forces of a flexible riser undergoing vortex-induced vibration(VIV) in sheared flow were investigated using a scale model test.The mean drag forces were determined based on beam theory using the strains measured by the fiber bragg grating strain sensors of the riser.The distributions of the drag coefficients along the riser and Reynolds numbers(<em>Re</em>), and the VIV amplification of the drag coefficient, were studied, and a new empirical model to estimate the drag coefficient on a flexible riser undergoing VIV was proposed.The results show that VIV can amplify the drag coefficient, resulting in the local drag coefficient of the riser up to 3.2.For <em>Re</em> values from 1.0×10<sup>4</sup>~1.2×10<sup>5</sup>, the mean drag coefficient value was between 1.3 and 2.0, and decreased as <em>Re</em> increased.Furthermore, the proposed empirical prediction model, which accounts for the effects of current, the VIV dominant modal number and the frequency in the cross-flow direction, can predict riser drag coefficients under VIV accurately at high <em>Re</em> values.</p>
南海内波流对深水钻井的影响及对策
[J].
Impacts of internal waves in the South China Sea on deepwater drilling safety and corresponding countermeasures
[J].
/
〈 |
|
〉 |
