螺旋列板立管受迫振动时的水动力学研究
上海交通大学 工程力学系; 水动力学教育部重点实验室, 上海 200240
Hydrodynamics Study of Riser with Helical Strakes Oscillating in Flow
Department of Engineering Mechanics; Key Laboratory of Hydrodynamics of the Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China
通讯作者: 孙 仁,男,教授,博士生导师,电话(Tel.):18121002767;E-mail:drrsun@sjtu.edu.cn.
责任编辑: 石易文
收稿日期: 2020-07-27
基金资助: |
|
Received: 2020-07-27
作者简介 About authors
李昂(1990-),男,河北省沧州市人,博士生,现主要从事流体力学研究
为了减小立管受到的振动损伤,本文对无列板立管以及三螺旋对称分布列板立管的水动力参数进行了实验评估.自然环境中的波浪与来流被抽象为振荡流与均匀流叠加,实验工况分为静止流场和相对运动流场.同时,立管受迫振动方向分别为流场纵向、横向以及纵向45° 夹角.结合Morison公式,由实验数据计算得到立管的附加质量系数Cm和托曳力系数Cd.实验结果表明,带螺旋列板立管的Cm与Keulegan-Carpenter(KC)数无关,Cm与受迫振动方向无关;Cd与流场速度和振动最大速度间的比率以及KC数的1/3次幂成反比关系.该发现与经典理论振荡流中平板结构Morison参数的形式相吻合,使经典理论在立管设计中获得新的验证与应用.在相同条件下,带螺旋列板立管的Cd值平均高出光滑立管273%.该发现表明在复杂工况下,螺旋列板能有效降低振动造成的影响.该结果可为海洋结构物的消振设计提供新的思路, 在海洋工程领域具有较高的应用价值.
关键词:
To reduce the vibration damage to the riser, hydrodynamic parameters of a smooth riser and a riser with triple symmetrically distributed helical strakes are evaluated in an experiment. The wave and current in natural environment are conceptualized to the experimental condition of the oscillatory flow and the uniform flow. Experiment cases are divided into the static flow and the uniform flow, while the risers oscillate along the in-line, the transverse, and the diagonal direction. The added mass coefficient Cm and the drag coefficient Cd are calculated from experimental data by using the Morison equation. The results of the riser with strakes indicate that Cm is independent of the Keulegan-Carpenter (KC) number and the oscillating direction. Cd is found to be in inverse proportion to the ratio of 1/3 power of the KC number and the maximum flow speed to oscillatory velocity. This finding is consistent with the format of Morison parameters of a flat plate in oscillatory flow, which endows the classical theory with novel validation and implication in the design of a riser. Under the same condition, the Cd of the helical strake riser is promoted over 273% than the bare riser, which indicates that the helical strakes efficiently reduce the influence of oscillation under complex conditions. The results provide a new solution for vibration reduction of offshore structures, which is of high value in ocean engineering.
Keywords:
本文引用格式
李昂, 孙仁.
LI Ang, SUN Ren.
海洋立管在复杂工况,如波浪与流共同作用下产生的振动会对其结构造成疲劳损害, 其来源包括圆柱体结构物尾流中形成交替泻涡导致的涡激振动[1],波浪以及海洋平台传导的振动等[2,3].如何更好地预报与控制海洋立管的振动是学术界和工业界的重要研究课题.目前,常用的抑制振动手段之一是在光滑圆管外侧加装螺旋列板[4],破坏漩涡沿立管轴向的泻涡三维关联性,从而达到减小振动的效果.Sui等[5]利用风洞试验测试了装有列板圆柱的受激振动,在雷诺数为2.7×103~5.79×104范围内,优化了列板设计参数.沙勇等[6]基于水池模型实验结果和工程设计经验,分析了螺旋列板抑制立管涡激振动的效果,对立管设计提供了相关参考.此外在数值计算方面,李艳潇等[7]采用大涡模拟方法, 对雷诺数为 3900 均匀来流下的串列双立管的涡激振动响应进行三维数值分析.针对不同的立管间距以及附加螺旋列板情况, 建模分析立管的水动力系数, 探究螺旋列板对双立管涡激振动的抑制效果.而对于柔性立管,近年来许多实验验证了列板亦能够有效抑制涡激振动[8,9].高云等[10]系统地研究了不同螺旋列板状态下立管应变、位移响应等参数, 发现流向响应与横向响应同样重要不可忽略,裸管响应特性与带螺旋列板的立管响应特性区别较大,响应特性与螺旋列板几何形状紧密相关.然而,加装螺旋列板会提高立管制造的难度以及相应成本,如何更高效地在立管轴向关键位置布置螺旋列板,以达到最大的抑制效果,对海洋油气开发设施的成本控制至关重要.因此,研究螺旋列板圆柱的水动力特性是海洋立管结构设计的重要前提条件[11].
尽管,至今已有大量涉及光滑与带列板的圆柱立管绕流研究,但是主要的工作集中在单纯均匀流引起的涡激振动的水动力响应.而在复杂工况,如均匀流与振荡流共同作用下造成振动的水动力参数的实验数据很少,亟待研究补充.针对波浪与流共同作用下的海洋立管振动问题,本文利用拖拽水池实验研究,把振荡流和均匀流叠加,建立了以Morison公式[20,21,22]为基础的螺旋列板圆柱在静水和均匀来流中受迫振动的模型,参数包含了拖曳力系数Cd、附加质量系数Cm以及Keulegan-Carpenter (KC)数与流场恒定速度、振动最大速度比率.实验结果与光滑圆柱实验数据进行了比较,并通过线性回归拟合发现了螺旋列板圆柱阻力系数基于KC数的两段区间表达式,为解决海洋立管受激振动的问题提供了重要参考.
1 圆柱立管两自由度运动的Morison参数
1.1 Morison公式
Morison公式通常被用于研究物体在振荡流中的水动力响应[23], 本文在此基础上将其拓展到振荡流与均匀流共同作用下的水动力分析.定义
式中:v为流体与物体的相对速度;
在纵向x方向,物体的位置、速度、加速度为
式中:
在横向y方向,物体的位置、速度、加速度为
式中:Ay为y方向振幅.
则物体速度v和加速度矢量
1.2 系数求解
在一个振动周期T中,力与速度的点积随时间积分可得:
同理,力与加速度积分可得:
测量的瞬态量可能存在随机误差,故本文统计时间平均意义上的物理量,把Morison公式代入式(6)和(7),得到Cd和Cm的表达式为
传感器测量立管瞬态受力F,并记录受迫振动参数v,根据式(8)和(9)即可得到时域平均的Cd、Cm.
1.3 主要参数
研究正弦振荡流体中的物体受力(或流体静止而物体振动)时,通常以无量纲KC数描述物体所受黏性力和相对惯性之间的关系,即:
式中:v0为振动速度的绝对值.
均匀来流速度v∞与v0之比,则有:
当
式中:ν为流体运动黏度.本实验的雷诺数范围为3.81×103~5.72×104.
2 实验装置与条件
2.1 实验装置
图1
2.2 实验条件
以两种形状的立管进行实验:光滑圆柱立管和螺旋列板圆柱立管.立管一端固定在拖曳臂,保持整体浸没水中.以KC数和v∞/v0无量纲数作为实验条件参数,实验工况参数如表1所示,其中N/A为空缺.实验中,分别测量立管两端受力(ATI Gamma 六分力传感器,ATI Industrial Automation, Inc., USA),同时记录下受迫振动的幅值和频率,共进行 2160 次实验,单次实验耗时约6 min.
表1 实验工况参数
Tab.1
KC | v∞/v0 | ||||||||
---|---|---|---|---|---|---|---|---|---|
工况1 | 工况2 | 工况3 | 工况4 | 工况5 | 工况6 | 工况7 | 工况8 | 工况9 | |
0.5 | 0.00 | N/A | N/A | N/A | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
1 | 0.00 | N/A | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
2 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
3 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
4 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
5 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
7 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
10 | 0.00 | 0.25 | 0.50 | 0.75 | 1.00 | 1.25 | 1.50 | 3.00 | 5.00 |
3 结果与讨论
3.1 光滑立管实验
首先对光滑圆柱在均匀流中横向简谐振动进行验证性实验
图2
图3
图4
图4
光滑立管在静止流场中受迫振动时的Morison参数
Fig.4
Morison parameters of a bare riser oscillating without current
当流体相对立管速度不为0时,本文测试了3种不同振动形式,分别为立管沿着流场横向振动,沿流场纵向振动以及与流向呈45° 夹角振动,得到了拖曳力系数、附加质量系数随振动系数v∞/v0变化的曲线,如图5~7所示.由图5~7可知,除KC=0.5外,Cd不受v∞/v0的影响,而Cm会随着v∞/v0的增加而下降.当KC数增加时,Cd也没有显著区分,而Cm则明显降低,同时Cm与振动方向的关联较大.Aronsen[25]研究了光滑圆柱在静止与均匀流场中的Cd,其雷诺数与本文大致相同,测得Cd值的范围为0.90~1.85,与本文实验结果接近.以上表明,对于光滑圆柱立管,在受迫振动情况下会持续受到显著的黏性拖曳力,同时附加质量系数也较低,这对海洋立管在流动中保持稳定是不利的.
图5
图5
光滑立管在匀速流场中沿流场横向振动时的Morison参数
Fig.5
Morison coefficients of a bare riser oscillating vertical with uniform current
图6
图6
光滑立管在匀速流场中沿流场纵向振动时的Morison参数
Fig.6
Morison coefficients of a bare riser oscillating parallel with uniform current
图7
图7
光滑立管在匀速流场中与流场呈45° 方向振动时的Morison参数
Fig.7
Morison coefficients of a bare riser oscillating along direction of 45° with uniform current
3.2 螺旋列板立管实验
用带螺旋列板的立管替换光滑圆柱形立管,其他实验条件与3.1节保持一致.当处于静止流场中时,圆柱立管的Cd与KC数的关系如图8(a)所示.Cd在KC=0.5时最大,最大值为12,并随着KC数的增加单调递减.KC数越低时,Cd越大,结构物附近的流体动能耗散越快,对振动产生的拖曳力的抵消效果越强.当KC=0.5时,螺旋列板立管Cd值比普通立管的Cd值高出8.75,百分比差距为269%;从整个KC数范围来看,螺旋列板立管Cd值比普通立管Cd值平均高出4.76.这说明特别是在KC数较低时,相较于光滑圆柱,螺旋列板结构更能够有效地抑制振动.另外,螺旋列板立管的附加质量系数Cm与KC数无明显关联(见图8(b)).
图8
图8
螺旋列板立管在静止流场中受迫振动时的Morison参数
Fig.8
Morison parameters of a helical strake riser oscillating without current
同样,当流体相对立管速度不为0时,实验工况分成3种不同振动形式,分别为立管沿着流场横向振动,沿流场纵向振动以及与流向呈45° 夹角振动,如图9~11所示.由图9~11可知,Cd 与振动方向关联较小,Cd 随着v∞/v0的增加而下降.而当v∞/v0相同时,KC数越小而Cd越大,意味着在高频振动、低流速时,螺旋列板的抑制振动效果更明显.Cm会随着v∞/v0的增加而略微降低,并且不受KC数影响,Cm与振动方向关联不大.与光滑立管数据比较可知,相同实验条件下螺旋列板结构的Cd明显大于光滑圆柱立管.Cd越高意味着流体的动能耗散更快,更有利于消除黏性拖曳力的影响,从而减小涡激振动对于立管结构的损害.特别当v∞/v0≤1时,该情况较多出现在深海(水流速度很低)或者波浪振动大于流致振动的海洋浅表区域,此时立管的Cd值远大于光滑圆柱,螺旋列板立管比普通立管的Cd平均值高出273%,这说明螺旋列板对于抑制波浪与流共同造成的振动有巨大的优势.
图9
图9
螺旋列板立管在匀速流场中沿流场横向振动时的Morison参数
Fig.9
Morison coefficients of a helical strake riser oscillating vertical with uniform current
图10
图10
螺旋列板立管在匀速流场中沿流场纵向振动时的Morison参数
Fig.10
Morison coefficients of a helical strake riser oscillating parallel with uniform current
图11
图11
螺旋列板立管在匀速流场中与流场呈45° 方向振动时的Morison参数
Fig.11
Morison coefficients of a helical strake riser oscillating along direction of 45° with uniform current
图12
本文拟合了
3 结语
本文利用实验方法测出光滑与三螺旋列板立管的水动力参数Cd、Cm,并对比了其差异.结果发现,螺旋列板立管的Cm不受KC数和受迫振动方向影响,说明它的振动固有频率比较稳定.螺旋列板立管的Cd随KC数的增加而下降,Cd随流速振动速度比率的增加而下降,Cd与受迫振动方向无关.相较于光滑圆柱立管,螺旋列板立管比普通立管的Cd平均提高273%,尤其在KC数较低时,螺旋列板立管的Cd相比光滑立管的绝对差值更大.Cd值高意味着流体动能耗散大,受黏性拖曳力影响小,是一种有效降低振动对结构影响的设计,最后本文拟合了螺旋列板立管的Cd与KC数的函数关系,并根据结构物在振荡流中的水动力学机理对其进行了解释.研究证明,螺旋列板立管可为海洋工程领域中包含波浪与洋流等因素导致的振动问题提供解决思路,有广阔的应用前景.
参考文献
A review on flow-induced vibration of offshore circular cylinders
[J]. ,DOI:10.1007/s42241-020-0032-2 URL [本文引用: 1]
随机波浪和涡流联合作用下海洋立管多频参数激励振动响应
[J]. ,
Multi-frequency parametrically excited vibrations of marine riser under simultaneous random waves and vortex
[J].
半主动调谐质量阻尼器控制管道振动实验研究
[J]. ,
Experiment study on semi-active pipeline vibration control with tuning mass dampers
[J].
The effectiveness of helical strakes in the suppression of high-mode-number VIV
[C]//Offshore Technology Conference. ,
VIV suppression for a large mass-damping cylinder attached with helical strakes
[J]. ,DOI:10.1016/j.jfluidstructs.2016.01.005 URL [本文引用: 1]
抑制涡激振动的螺旋列板设计参数研究
[J]. ,
Helical strakes for marine risers VIV suppression and its design parameters selection
[J].DOI:10.1016/S0029-8018(03)00107-0 URL [本文引用: 1]
串列双立管螺旋列板抑制涡激振动的数值模拟
[J]. ,
Three dimensional numerical simulation of VIV of tandem double risers by helical strakes
[J].
Investigation on the effectiveness of helical strakes in suppressing VIV of flexible riser
[J]. ,DOI:10.1016/j.apor.2013.11.006 URL [本文引用: 1]
VIV response of a long flexible riser fitted with strakes in uniform and linearly sheared currents
[J]. ,DOI:10.1016/j.apor.2015.05.006 URL [本文引用: 1]
柔性立管涡激振动抑制装置试验研究
[J]. ,
Experimental investigation on the suppression device of VIV of a flexible riser
[J].
Hydrodynamic coefficients for risers with strakes
[C]//Proceedings of ASME 2002 21st International Conference on Offshore Mechanics and Arctic Engineering. ,
A robotic Intelligent Towing Tank for learning complex fluid-structure dynamics
[J]. ,
Experimental investigation on vortex induced forces of oscillating cylinder at high Reynolds number
[J]. ,DOI:10.12989/ose.2013.3.3.167 URL [本文引用: 1]
均匀流中近壁面垂直流向振荡圆柱水动力特性研究
[J]. ,
Hydrodynamic characters of a near-wall circular cylinder oscillating in cross flow direction in steady current
[J].
Dynamic response and hydrodynamic coefficients of a cylinder oscillating in crossflow with an upstream wake interference
[J]. ,DOI:10.1016/j.oceaneng.2020.107520 URL [本文引用: 1]
Kill line model cross flow inline coupled vortex-induced vibration
[C]//Proceedings of ASME 2017 36th International Conference on Ocean, Offshore and Arctic Engineering. ,
Experimental investigation of cross flow-inline coupled Vortex-Induced Vibration on riser with finite length buoyancy module
[C]//OCEANS 2016 MTS/IEEE Monterey. ,
Mapping the properties of the vortex-induced vibrations of flexible cylinders in uniform oncoming flow
[J]. ,DOI:10.1017/jfm.2019.738 URL [本文引用: 1]
Vortex induced vibration of riser with low span to diameter ratio buoyancy modules
[C]//The 27th International Ocean and Polar Engineering Conference. ,
Forces on cylinders and plates in an oscillating fluid
[J]. ,DOI:10.6028/jres.060.043 URL [本文引用: 1]
Force on a circular cylinder in viscous oscillatory flow at low Keulegan-Carpenter numbers
[J]. ,DOI:10.1017/S0022112086002999 URL [本文引用: 1]
Drag coefficient enhancement of dual cylinders in oscillatory flow
[C]//The 27th International Ocean and Polar Engineering Conference. ,
Sinusoidal flow relative to circular cylinders
[J]. ,DOI:10.1017/S0022112085001756 URL [本文引用: 1]
An experimental investigation of in-line and combined in-line and cross-flow vortex induced vibrations
[D]. ,
The forces on sharp-edged cylinders in oscillatory flow at low Keulegan-Carpenter numbers
[J]. ,DOI:10.1017/S0022112080002595 URL [本文引用: 1]
/
〈 | 〉 |