CROSS-FLOW VORTEX-INDUCED VIBRATION OF A FLEXIBLE RISER WITH INTERNAL FLOW IN SHEAR CURRENT
-
摘要: 立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉−伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动；立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显；立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.Abstract: As an important component transporting resources such as oil and mineral ores mixture from the seabed to the surface in ocean engineering, vortex-induced vibration (VIV) of flexible risers can be encountered when the risers are subjected to the external environmental conditions. As VIV can lead to structural fatigue for the riser system, which threatens to the facility safety during deepsea resource exploitation, it is of great significance to investigate VIV mechanism and dynamics. Therefore, VIV dynamics of a flexible fluid-conveying riser undergoing external shear current is studied based on the combination of the Euler-Bernoulli beam theory and the semi-empirical hydrodynamic model. The finite element method and Newmark-β method are adopted to discretize and solve the governing equation. The model is firstly validated by comparing with the experimental data in order to examine the accuracy of the present model. Subsequently, cross-flow (CF) VIV response of the fluid-conveying riser is mainly examined and analyzed while various internal flow velocity and fluid density are considered and changed. The results show that when the flexible riser is subjected to both internal flow and shear current, there appears multi-frequency response for CF VIV. And the CF vibrating frequency and the CF root mean square (RMS) displacement are evidently influenced by the internal flow velocity and fluid density. With the increase of the internal flow velocity and fluid density, the CF vibrating frequency decreases while the RMS displacement shows an increasing trend in CF direction. Furthermore, in addition to the variation of the CF vibrating frequency and RMS displacement, the change of internal flow densities can cause notable mode and frequency transitions.
-
引 言
海洋立管作为深海资源开发中输送海底油气或其他矿产混合物的重要部件, 一旦在海洋波浪和洋流的作用下发生涡激振动, 则极易引起立管的疲劳破坏. 因此海洋立管的涡激振动响应得到国内外学者的持续关注, 对海洋立管在外部环境荷载作用下产生的涡激振动的机理和响应已有大量的研究, 特别是针对大长细比的柔性立管[1-3]. 这些研究主要运用实验和数值模拟方法, 针对均匀流和剪切流中的柔性立管, 分析了涡激振动的模态、频率以及驻波和行波响应等, 为深入理解涡激振动机理和响应规律奠定了基础[1-3]. 在这些研究中, 一些学者关注了外部条件对立管涡激振动响应的影响规律, 如来流角度和边界条件等, 其中, Chaplin等[4]研究了不同来流流向对悬链线立管水动力系数以及振动响应的影响机理; Xu等[5]、高云等[6]和徐万海等[7]则通过实验和数值的方法研究分析了不同来流角度和不同剪切流中柔性圆柱体流体力特性及其系数的变化规律, 并且对结构的能量输入区域和输出区域进行了划分; Chen等[8]关注了三角形刚性柱体在不同来流流场中的振动响应和流场特性, 分析了不同来流角度下柱体在不同约化速度下的振动响应, 并且着重关注了随着约化速度增加产生的驰振现象; 李非凡等[9]和及春宁等[10-11]采用浸入边界法对细长柔性圆柱在均匀流、剪切流以及倾斜流条件下的涡激振动进行三维数值模拟, 主要研究了涡激振动的频率、阻力系数、能量传递以及主播和行波响应等; 同时, Seyed-Aghazadeh等[12]研究分析了不同边界条件对涡激振动响应的影响机理和规律. 除此之外, Song等[13]通过实验发现立管顺流向和横向的多模态响应并且立管振动频率与涡泻脱落频率一致, Hu等[14]进一步研究了不同质量比和阻尼比下的串列立管涡激振动响应, 并且发现质量比小的情况下立管涡激振动迟滞现象明显. 在柔性立管涡激振动数值模拟研究中, Bao等[15]运用厚切片理论预报了涡激振动的结构响应、水动力系数以及涡泻脱落形态等; Qu等[16-17]则改进了尾流振子模型, 数值预报了柔性圆柱体的锁频现象和相位角变化; Yuan等[18]和Zhang等[19]基于半经验时域模型研究了立管在稳定流和振荡流中的涡激振动响应, 主要分析其振动、频率和模态随不同参数的变化规律, Lin等[20]运用SDVM-FEM耦合方法研究了柔性立管产生的驻波和行波响应, 并且分析了其相应的流场特性, 这些研究为海洋工程界深入理解和认识立管涡激振动机理和响应规律奠定了坚实的基础. 然而, 这些研究主要考虑外流的作用, 海洋立管输送海底油气或其他矿产资源时, 内流对立管水动力响应的影响不可忽视, 因此, 研究含内流立管的涡激振动具有重要意义.
当介质在立管内部流动时, 立管受到流动介质产生的内部作用力, 可导致海洋立管复杂的涡激振动响应. 近几年, 越来越多的学者开始对考虑内外流共同影响下海洋立管的涡激振动进行探究. 研究证明, 在一定的内流速度下, 管道结构发生失稳现象[21]. 同时, 一些学者通过实验和数值模拟方法研究了含内流海洋立管涡激振动响应, Guo等[22]、李朋[23]和吴天昊[24]通过模型试验, 研究了含内流顶张立管和悬臂立管的动力响应, 主要关注立管振动频率和位移变化规律. Jiang等[25]、Duan等[26-27]、Thorsen等[28]、Wang等[29]和李艳等[30]则运用数值手段研究了含稳定内流海洋立管的涡激振动响应, 分析了内流速度和密度对立管涡激振动模态、频率和位移的影响机理. 然而在实际海洋工程中, 海洋立管输送的物质存在多相流动, 因此, Zhu等[31]、Meng等[32]和Li等[33]关注了含气液两相内流立管振动响应, 主要研究了不同的气液两相内流作用下立管对外流的动力响应. 这些研究证明, 在内流的影响下, 海洋立管涡激振动特性更加复杂, 涡激振动的频率和主导模态等都会发生变化.
虽然含内流海洋立管涡激振动特性研究已经引起国内外学者的广泛关注, 但是目前对含内流海洋立管涡激振动响应的研究刚刚起步, 相关的文献报道并不多见, 特别是剪切流作用的情况更是鲜有报道. 鉴于此, 本文旨在研究含内流海洋立管在剪切外流作用下的横向涡激振动响应, 分析内流对立管在剪切流中涡激振动响应机理和特性的影响.
1. 立管涡激振动数学模型
本研究主要针对两端简支立管(如图1所示), 研究其横向涡激振动响应, 立管被施以预张力
$T$ , 并且含有速度为${U_{\rm{i}}}$ 的内流. 含内流立管横向涡激振动控制方程可表示为[21]$$\begin{split} & \left( {{m_{\rm{r}}} + {m_{\rm{f}}}} \right)\frac{{{\partial ^2}y({\textit{z}},t)}}{{\partial {t^2}}} + c\frac{{\partial y({\textit{z}},t)}}{{\partial t}} + 2{m_{\rm{f}}}{U_{\rm{i}}}\frac{{{\partial ^2}y({\textit{z}},t)}}{{\partial {\textit{z}}\partial t}} +\\ & \qquad \left( {{m_{\rm{f}}}U_{\rm{i}}^2 - T} \right)\frac{{{\partial ^2}y({\textit{z}},t)}}{{\partial {{\textit{z}}^2}}} + EI\frac{{{\partial ^4}y({\textit{z}},t)}}{{\partial {{\textit{z}}^4}}} = {F_{{\rm{CF}}}}({\textit{z}},t) \\[-8pt] \end{split} $$ (1) ![]()
图 1 含内流立管横向涡激振动示意图Figure 1. Schematic of vortex-induced vibration (VIV) for a riser considering internal flow. CF stands for cross-flow其中,
${m_{\rm{r}}}$ 和${m_{\rm{f}}}$ 为单位长度立管和内流的质量,$c$ 为立管的阻尼系数,${U_{\rm{i}}}$ 为内流速度,$EI$ 为弯曲刚度,${F_{{\rm{CF}}}}$ 为外部流体力. 方程中第3项和第4项表示立管内流对立管涡激振动响应的影响, 其中第3项为内流对立管所施加的科氏力, 而第4项则表示内流流动所产生的离心力对立管涡激振动的影响.${F_{{\rm{CF}}}}$ 包含附加质量力和涡激力$${F_{{\rm{CF}}}} = - {m_{{\rm{a}}}}\ddot y + {F_{{\rm{ex,CF}}}}$$ (2) 其中, 附加质量
${m_{{\rm{a}}}} = {C_{{\rm{a}}}}{m_{\rm{e}}}$ , 附加质量系数${m_{{\rm{a}}}}$ 取为1.0,${m_{\rm{e}}}$ 为立管排开水单位长度质量.为了确定外部涡激力, 需要首先对立管在横向上的激励区和阻尼区进行判定, 可通过确定立管上无因次激励频率来判断激励区和阻尼区(如图2).
![]()
图 2 含内流立管横向激励区和阻尼区示意图Figure 2. Identification of excitation and damping regions of a fluid-conveying riser由模态分析可得到立管的各阶固有频率, 再结合流场并基于方程(2)获得立管各节点在每阶固有频率下的无因次频率
$$f_{{\rm{CF}},i}^*({\textit{z}}) = \frac{{{\omega _{{\rm{CF}},i}}}}{{2{\text{π}} }}\frac{{{D_{{\rm{e}}}}}}{{{U_{{\rm{e}}}}({\textit{z}})}}$$ (3) 其中,
$f_{{\rm{CF}},i}^*({\textit{z}})$ 表示第$ i $ 阶固有频率下节点${\textit{z}}$ 的无因次频率,${D_{{\rm{e}}}}$ 为立管的直径,${U_{{\rm{e}}}}({\textit{z}})$ 为外部剪切流在节点${\textit{z}}$ 处流速,${\omega _{{\rm{CF}},i}}$ 表示第$i$ 阶固有圆频率. 由于模型中的斯托哈尔数并不是恒定的, 根据实际雷诺数对节点${\textit{z}}$ 的斯托哈尔数进行修正[34].剪切流中, 沿立管轴向都存在激励区和阻尼区. 基于涡激振动的半经验频域预报模型VIVANA, 可根据以下标准来判定立管的激励区和阻尼区
$$0.125 \leqslant f_{{\rm{CF}},i}^*({\textit{z}}) \leqslant 0.2$$ (4) 当立管在节点
${\textit{z}}$ 处第$i$ 阶无因次激励频率在0.125 ~ 0.2区间内时, 则该节点处于第$ i $ 阶激励频率下的激励区内, 相反该节点则位于该阶激励频率的阻尼区内. 然而, 在多频响应下, 直接通过该标准得到的各激励频率的激励区可能发生重叠. 针对发生重叠的激励区, 本文采用等长度退让的方法处理多频响应下激励区发生重叠的问题[34].确定立管横向上的激励区后, 激励区内节点
${\textit{z}}$ 处的第$i$ 阶激励力可以表示为$$\begin{split} & {F_{{\rm{{\rm{ex,CF}}}}}}({\textit{z}}) = \\ & \qquad \frac{1}{2}{\rho _{\rm{e}}}{D_{{\rm{e}}}}U({\textit{z}})_{\rm{e}}^2{C_{{\rm{L}},{\rm{CF}},i}}({\textit{z}})\cos ({\omega _{{\rm{CF}},i}} \cdot t + {\varphi _{{\rm{CF}},i}}) \\[-6pt] \end{split} $$ (5) 其中,
${\rho _{\rm{e}}}$ 为外部流体密度,${C_{{\rm{L}},{\rm{CF}},i}}({\textit{z}})$ 为节点${\textit{z}}$ 处激励力系数,${\varphi _{{\rm{CF}},i}}$ 为相位角. 基于Gopalkrishnan[35]刚性圆柱体强迫振荡试验结果(图3), 根据同一无因次频率下的激励力系数${C_{\rm{L}}}$ 与无因次幅值之间的关系曲线(图4)可确定${C_{{\rm{L}},{\rm{CF}},i}}({\textit{z}})$ .![]()
相位角
${\varphi _{{\rm{CF}},i}}$ 则是根据Zhang等[19]提出的方法确定, 依不同的振型, 可以分别设为0和${\text{π}} $ .在阻尼区节点
${\textit{z}}$ 处的阻尼力则采用VIVANA中水动力阻力模型[35]$${F_{{\rm{damp}} ,{\rm{CF}}}}({\textit{z}}) = - {C_{{\rm{damp}} }}({\textit{z}})\frac{{\partial y({\textit{z}},t)}}{{\partial t}}$$ (6) 其中,
${F_{{\rm{damp}} ,{\rm{CF}}}}({\textit{z}})$ 为节点${\textit{z}}$ 处第$i$ 阶激励频率的水动力阻力,${C_{{\rm{damp}}}}({\textit{z}})$ 为阻尼区内阻尼系数.采用有限单元法对立管横向振动控制方程进行离散, 运用Newmark-β逐步积分方法求解方程, 进而研究含内流海洋立管横向涡激振动特性.
2. 模型验证
运用有限元方法将模型进行离散, 将模型离散为100个单元, 如图5所示, 之后运用Newmark-β法求解离散方程.
![]()
图 5 立管离散单元示意图Figure 5. Schematic of discretized elements of the fluid-conveying riser首先, 通过将数值计算结果与实验结果进行对比, 验证模型的准确性, 本文中选用的对比数据为宋磊建的模型实验数据[36], 该模型实验的具体参数如表1.
Parameter Value length L/m 6.75 outer diameter De/m 0.031 inner diameter di/m 0.027 bending stiffness EI/N·m2 1476.76 pretension T/N 3000 mass per unit mr/kg 1.768 damping ratio c/% 0.3 如图6所示, 在外部剪切流流速为1.5 m/s (线性剪切流中最大流速为1.5 m/s, 最小流速为0)情况下, 不含内流的立管在横向上的位移均方根与实验结果吻合较好, 数值模拟和实验中立管涡激振动的主导振动模态都为二阶模态. 同时, 数值模拟中立管节点25在横向上的主导振动频率为6.1 Hz (图7), 而实验中立管在横向上的主导振动频率为5.4 Hz. 值得注意的是立管在剪切流作用下会产生多频率涡激振动, 从图7中可以看出, 立管在横向上存在另一个振动频率2.7 Hz, 这是由于外部剪切流流速沿立管轴向变化, 导致立管不同位置产生不同的模态响应, 从而引起立管整体发生多模态多频率涡激振动响应. 然而, 由于数值模拟中所选取的附加质量系数和激励系数等原因, 数值模拟结果与实验结果存在一定的误差, 但是误差都较小.
![]()
图 7 横向振动时历曲线和振动频率Figure 7. Time history of vibration and vibrating frequency in CF direction3. 含内流海洋立管横向涡激振动特性
3.1 内流对立管固有频率的影响
立管结构的固有频率可以依据下式求得[37]
$$\left| {{\boldsymbol{K}} - {\omega ^2}{\boldsymbol{M}}} \right| = 0$$ (7) 式中,
$\omega $ 表示立管各阶圆频率,${\boldsymbol{K}}$ 和${\boldsymbol{M}}$ 分别为立管的刚度矩阵和质量矩阵, 其中立管的刚度矩阵受内流速度和质量的影响. 含内流立管的固有频率随着内流速度和密度变化规律如图8和图9所示.![]()
图 8 立管固有频率随着内流速度变化规律Figure 8. Variation of natural frequency with the increase of the internal flow velocity![]()
图 9 立管固有频率随着内流密度变化规律Figure 9. Variation of natural frequency with the increase of the internal fluid density从图中可以看出内流的速度和密度影响立管固有频率, 随着内流密度和速度增大, 立管的前3阶固有频率降低. 若内流密度设为2000 kg/m3, 当内流速度从0 m/s增大到50 m/s时, 立管的前3阶固有频率分别从2.33 Hz, 5.23 Hz和9.12 Hz减小到1.06 Hz, 3.2 Hz和6.68 Hz; 若内流速度设为10 m/s, 当内流密度从1000 kg/m3增大到4000 kg/m3, 立管的前3阶固有频率分别从2.51 Hz, 5.66 Hz和9.88 Hz减小到1.96 Hz, 4.45 Hz和7.81 Hz. 从式(7)看出, 立管内流的速度和密度增大可导致立管刚度减小, 从而减小立管固有频率.
3.2 横向振动频率
3.2.1 内流速度的影响
将立管内流密度设为2000 kg/m3, 立管外部剪切流速度为1.5 m/s, 图10展示了随着内流速度的增大, 立管节点50处横向振动频率的变化规律. 由于剪切流的作用, 立管在横向上出现多频率多模态涡激振动现象.
![]()
图 10 横向振动频率随着内流速度的变化趋势Figure 10. Variation of CF vibrating frequency with the increase of the internal flow velocity从图10可以看出, 立管横向振动的主导频率随着内流速度的增大而减小. 当立管内流速度从0增大到10 m/s时, 横向主导振动频率减小的幅度较小, 大约从9.15 Hz减小到9.08 Hz; 而当内流速度增大到20 m/s和30 m/s时, 横向主导振动频率随着内流速度的增大明显降低, 大约分别降低到8.8 Hz和8.2 Hz. 出现这种现象的原因是, 当存在内流时, 立管的固有频率随着内流速度的增大而降低, 立管涡泻频率在特定的内流速度下锁定到新的较低的固有频率, 从而引起在此频率下的涡激振动响应.
3.2.2 内流密度的影响
由于立管系统在深海资源开发中输送油气混合物以及固液矿产混合物, 所以需要研究不同内流密度对立管涡激振动的影响规律. 同样在立管外部剪切流速度为1.5 m/s, 将立管内流速度设为10 m/s, 分别研究内流密度1000 kg/m3, 2000 kg/m3和3000 kg/m3对立管横着涡激振动参数的影响规律, 其中密度1000 kg/m3接近油气混合物密度, 而密度2000 kg/m3和3000 kg/m3则接近固液矿产混合物密度. 内流密度的变化对立管横向涡激振动频率影响规律如图11所示. 从图11可以看出, 当立管内流密度为1000 kg/m3时, 立管横向振动的主导频率约为5.6 Hz, 当内流密度增大到2000 kg/m3时, 其主导频率转换到9.08 Hz. 出现这种现象的原因与新的主导模态响应激发以及固有频率的降低有关, 当内流速度从1000 kg/m3增大到2000 kg/m3时, 立管横向涡激振动的主导模态从二阶响应变换到3阶模态响应, 从而导致立管在横向上的主导振动频率从5.6 Hz跳跃到9.08 Hz; 而随着内流密度的继续增大, 立管固有频率降低并且没有发生模态响应转换, 从而导致立管横向涡激振动响应与新的固有频率锁定, 因此当内流密度增大到3000 kg/m3时, 立管横向上的主导振动频率降为8.4 Hz. 值得注意的是, 立管涡激振动的频率转换与模态转换同时发生.
![]()
图 11 横向振动频率随着内流密度的变化趋势Figure 11. Variation of CF vibrating frequency with the increase of the internal fluid density3.3 横向均方根位移
3.3.1 内流速度的影响
立管横向均方根位移随着内流速度的变化规律如图12所示. 从图12中可以看出, 随着内流速度的增大, 立管在横向上的最大均方根位移是不断增大的, 值得注意的是, 均方根位移的增大在内流速度较大的时候比较明显, 而内流速度较小时, 增大较慢. 当内流速度从0增大到10 m/s的时候, 立管横向的最大均方根位移几乎不变, 大约为0.24; 而当内流速度增大到20 m/s和30 m/s时, 立管横向的最大均方根位移显著增大, 大约为0.27和0.33. 立管横向最大均方根位移的增大是由于内流对立管涡激振动产生影响.
![]()
图 12 横向均方根位移随着内流速度变化趋势Figure 12. Variation of CF RMS with the increase of the internal flow velocity约化速度是研究立管涡激振动特性的一个重要参数, 有如下表达形式
$${U_{\rm{r}}} = \frac{{{U_{\rm{e}}}}}{{{f_n}{D_{\rm{e}}}}}$$ (8) 其中,
${f_n}$ 是立管的第n阶固有频率. 立管横向涡激振动振幅随着约化速度${U_{\rm{r}}}$ 先增大后减小, 并且一般在约化速度为6时达到最大值[38]. 当立管内流速度从0增大到10 m/s, 约化速度大约从4.3增大到4.45, 进而导致立管横着涡激振动振幅增大, 随着内流速度继续增大到20 m/s和30 m/s, 约化速度大约增大到4.6和4.8, 相应地, 立管在横向上的涡激振动振幅增大. 由此可见, 内流流速通过影响约化速度, 进而改变立管横向涡激振动振幅. 值得注意的是, 剪切流引起的立管横向涡激振动响应沿轴向一般应该是不对称的, 但本研究中立管的长度较短, 这种不对称特征未能充分体现.3.3.2 内流密度的影响
而立管横向均方根位移随内流密度的变化如图13所示. 从图13中可以看出, 当不发生模态转换时, 立管在横向上的最大均方根位移随内流密度的增大而增大, 当内流密度为2000 kg/m3时, 立管在横向上的最大均方根位移大约为0.24, 当内流密度增大到3000 kg/m3时, 其值约为0.35. 值得注意的是, 当立管在横向上发现模态转换时, 横向最大均方根位移值发生跳跃现象, 如当立管由二阶模态响应变为3阶模态响应, 其横向均方根位移最大值约从0.3 (1000 kg/m3)跳跃到0.24 (2000 kg/m3). 横向最大均方根位移随内流密度的增大而增大同样是由于密度变化引起约化速度变化, 从而导致立管横向涡激振动振幅改变. 当内流密度从1000 kg/m3增大到2000 kg/m3时, 其约化速度大约从4.8减小到4.45, 进而导致立管在横向上的涡激振动振幅减小, 而当内流密度继续增大到3000 kg/m3, 约化速度跳跃到接近5, 立管振动振幅相应增大.
![]()
图 13 横向均方根位移随着内流密度变化趋势Figure 13. Variation of CF RMS with the increase of the internal fluid density值得注意的是, 本研究中, 立管的主导振动模态在内流速度的影响下并没有改变, Duan等[26-27]证明了随着内流速度的增大, 立管在横向上的主导模态是增大的, 但是这种模态转换只是在内流速度较大时才发生. 显然, 本研究中内流速度还没有引起立管在横向上振动的模态转换. 而内流密度的增大对立管横向振动的主导模态影响明显, 如图13所示. 当内流密度为1000 kg/m3时, 立管横向振动的主导模态为二阶振动响应; 当内流密度增大到2000 kg/m3和3000 kg/m3时, 立管横向振动发生模态转化, 变为3阶模态主导响应. 发生模态转换的原因同样是由于内流密度增大降低立管固有频率, 从而激发新的涡激振动模态响应, 进而发生模态转换.
4. 结论
本文基于欧拉−伯努利梁理论和半经验时域水动力模型, 研究了内外流耦合作用下海洋立管在剪切流场中的横向涡激振动响应, 分析和讨论了立管在横向上的均方根位移、主导振动模态和频率等随内流速度和密度的变化规律.
海洋立管在剪切流作用下发生多模态多频率涡激振动响应, 深水立管在剪切流作用下的横向涡激振动响应沿立管轴向一般是不对称的, 立管较短时, 这种不对称特征不明显. 内流速度和密度的增大, 导致立管固有频率降低, 因此涡激力频率锁定在新的较低固有频率上, 激发新的横向涡激振动模态响应, 进而改变立管横向涡激振动振动频率, 导致立管在横向上的主导振动频率随内流速度和密度的增大而减小. 并且, 内流密度的增大可引起立管发生模态转换和频率转换. 同时, 立管固有频率的变化影响约化速度, 进而改变立管横向涡激振动均方根位移, 因此, 立管横向最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度和密度较大时, 横向最大均方根位移增大明显.
-
![]()
图 1 含内流立管横向涡激振动示意图
Figure 1. Schematic of vortex-induced vibration (VIV) for a riser considering internal flow. CF stands for cross-flow
![]()
图 2 含内流立管横向激励区和阻尼区示意图
Figure 2. Identification of excitation and damping regions of a fluid-conveying riser
![]()
![]()
图 5 立管离散单元示意图
Figure 5. Schematic of discretized elements of the fluid-conveying riser
![]()
图 7 横向振动时历曲线和振动频率
Figure 7. Time history of vibration and vibrating frequency in CF direction
![]()
图 8 立管固有频率随着内流速度变化规律
Figure 8. Variation of natural frequency with the increase of the internal flow velocity
![]()
图 9 立管固有频率随着内流密度变化规律
Figure 9. Variation of natural frequency with the increase of the internal fluid density
![]()
图 10 横向振动频率随着内流速度的变化趋势
Figure 10. Variation of CF vibrating frequency with the increase of the internal flow velocity
![]()
图 11 横向振动频率随着内流密度的变化趋势
Figure 11. Variation of CF vibrating frequency with the increase of the internal fluid density
![]()
图 12 横向均方根位移随着内流速度变化趋势
Figure 12. Variation of CF RMS with the increase of the internal flow velocity
![]()
图 13 横向均方根位移随着内流密度变化趋势
Figure 13. Variation of CF RMS with the increase of the internal fluid density
-
[1] Sarpkaya T. A critical review of the intrinsic nature of vortex-induced vibrations. Journal of Fluids and Structures, 2004, 19: 389-447
[2] Williamson CHK, Govardhan R. A brief review of recent results in vortex-induced vibrations. Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96: 713-735
[3] Bearman PW. Circular cylinder wakes and vortex-induced vibrations. Journal of Fluids and Structures, 2011, 27: 648-658
[4] Chaplin JR, King R. Laboratory measurements of the vortex-induced vibrations of an untensioned catenary riser with high curvature. Journal of Fluids and Structures, 2018, 79: 26-38 doi: 10.1016/j.jfluidstructs.2018.01.008
[5] Xu WH, Ma YX, Ji CN, et al. Laboratory measurements of vortex-induced vibrations of a yawed flexible cylinder at different yaw angles. Ocean Engineering, 2018, 154: 27-42 doi: 10.1016/j.oceaneng.2018.01.113
[6] 高云, 邹丽, 宗智. 两端铰接的细长柔性圆柱体涡激振动响应特性数值研究. 力学学报, 2018, 50(1): 9-20 (Gao Yun, Zou Li, Zong Zhi. Numerical study of response performance of vortex-induced vibration on a flexible cylinder with pinned-pinned boundary condition. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(1): 9-20 (in Chinese) doi: 10.6052/0459-1879-17-340 [7] 徐万海, 马烨璇, 罗浩等. 柔性圆柱涡激振动流体力系数识别及其特性. 力学学报, 2017, 49(4): 818-827 (Xu Wanhai, Ma Yexuan, Luo Hao, et al. Identification and characteristics of hydrodynamic coefficients for a flexible cylinder undergoing vortex-induced vibration. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(4): 818-827 (in Chinese) [8] Chen WL, Ji CN, Xu D, et al. Flow-induced vibrations of an equilateral triangular prism at various angles of attack. Journal of Fluids and Structures, 2020, 97: 103099 doi: 10.1016/j.jfluidstructs.2020.103099
[9] 李非凡, 及春宁, 许栋等. 细长柔性立管涡激振动的数值模拟研究. 工程力学, 2017, 34(4): 221-230 (Li Feifan, Ji Chunning, Xu Dong, et al. Numerical investigation on the vortex-induced vibration of a flexible slender riser. Engineering Mechanics, 2017, 34(4): 221-230 (in Chinese) [10] 及春宁, 花阳, 许栋等. 不同剪切率来流作用下柔性圆柱涡激振动数值模拟. 力学学报, 2018, 50(1): 21-31 (Ji Chunning, Hua Yang, Xu Dong, et al. Numerical simulation of vortex-induced vibration of a flexible cylinder exposed to shear flow at different shear rates. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(1): 21-31 (in Chinese) doi: 10.6052/0459-1879-17-222 [11] 及春宁, 邢国源, 张力等. 倾斜流作用下柔性立管涡激振动的数值模拟. 哈尔滨工程大学学报, 2018, 39(2): 324-331 (Ji Chunning, Xing Guoyuan, Zhang Li, et al. Numerical simulations of vortex-induced vibration of flexible riser subjected to inclined flow. Journal of Harbin Engineering University, 2018, 39(2): 324-331 (in Chinese) [12] Seyed-Aghazadeh B, Edraki M, Modarres-Sadeghi Y. Effects of boundary conditions on vortex-induced vibration of a fully submerged flexible cylinder. Experiments in Fluids, 2019, 60: 38
[13] Song JN, Lu L, Teng B, et al. Laboratory tests of vortex-induced vibrations of a long flexible riser pipe subjected to uniform flow. Ocean Engineering, 2011, 38: 1308-1322
[14] Hu Z, Wang J, Sun Y. Flow-induced vibration of one-fixed-one-free tandem arrangement cylinders with different mass-damping ratios using wind tunnel experiment. Journal of Fluids and Structures, 2020, 96: 103019 doi: 10.1016/j.jfluidstructs.2020.103019
[15] Bao Y, Zhu HB, Huan P, et al. Numerical prediction of vortex-induced vibration of flexible riser with thick strip method. Journal of Fluids and Structures, 2019, 89: 166-173 doi: 10.1016/j.jfluidstructs.2019.02.010
[16] Qu Y, Metrikine AV. A single van der pol wake oscillator model for coupled cross-flow and in-line vortex-induced vibrations. Ocean Engineering, 2020, 196: 106732 doi: 10.1016/j.oceaneng.2019.106732
[17] Qu Y, Metrikine AV. A wake oscillator model with nonlinear coupling for the vortex-induced vibration of a rigid cylinder constrained to vibrate in the cross-flow direction. Journal of Sound and Vibration, 2020, 469: 115161 doi: 10.1016/j.jsv.2019.115161
[18] Yuan YC, Xue HX, Tang WY. Numerical analysis of vortex-induced vibration for flexible risers under steady and oscillatory flows. Ocean Engineering, 2018, 148: 548-562 doi: 10.1016/j.oceaneng.2017.11.047
[19] Zhang MM, Fu SX, Song LJ, et al. A time domain prediction method for the vortex-induced vibrations of a flexible riser. Marine Structures, 2018, 59: 458-481
[20] Lin K, Wang J. Numerical simulation of vortex-induced vibration of long flexible risers using a SDVM-FEM coupled method. Ocean Engineering, 2019, 172: 468-486 doi: 10.1016/j.oceaneng.2018.12.006
[21] Païdoussis MP. Fluid-structure Interactions, Slender Structures and Axial Flow. volume 1. second ed.. California: Academic Press, 2014
[22] Guo HY, Lou M. Effect on internal flow on vortex induced vibration of risers. Journal of Fluids and Structures, 2008, 24(4): 496-504 doi: 10.1016/j.jfluidstructs.2007.10.002
[23] 李朋. 深水海洋输流立管涡激振动及干涉试验研究. [博士论文] . 青岛: 中国海洋大学, 2012 (Li Peng. Experimental study on vortex-induced vibration and interference of deepwater marine risers transporting fluid. [PhD Thesis]. Qingdao: Ocean University of China, 2012 (in Chinese))
[24] 吴天昊. 悬垂立管涡激振动响应特性研究. [硕士论文]. 上海: 上海交通大学, 2018 (Wu Tianhao. The study on VIV responses of a free-hang riser. [Master Thesis]. Shanghai: Shanghai Jiao Tong University, 2018 (in Chinese))
[25] Jiang TL, Dai HL, Wang L. Three-dimensional dynamics of fluid-conveying pipe simultaneously subjected to external axial flow. Ocean Engineering, 2020, 217: 107970 doi: 10.1016/j.oceaneng.2020.107970
[26] Duan JL, Chen K, You YX, et al. Numerical investigation of vortex-induced vibration of a riser with internal flow. Applied Ocean Research, 2018, 72: 110-121 doi: 10.1016/j.apor.2018.01.003
[27] Duan JL, Zhou JF, You YX, et al. Time-domain analysis of vortex-induced vibration of a flexible mining riser transporting flow with various velocities and densities. Ocean Engineering, 2021, 220: 108427 doi: 10.1016/j.oceaneng.2020.108427
[28] Thorsen MJ, Challabotla NR, Sævik S, et al. A numerical study on vortex-induced vibrations and the effect of slurry density variations on fatigue of ocean mining risers. Ocean Engineering, 2019, 174: 1-13 doi: 10.1016/j.oceaneng.2019.01.041
[29] Wang L, Jiang TL, Dai HL, et al. Three-dimensional vortex-induced vibrations of supported pipes conveying fluid based on wake oscillator models. Journal of Sound and Vibration, 2018, 422: 590-612 doi: 10.1016/j.jsv.2018.02.032
[30] 李艳, 成赟, 郭轶可等. 内外流共同作用下输流软管力学特性研究. 水资源与水工程学报, 2018, 29(3): 182-193 (Li Yan, Cheng Yun, Guo Yike, et al. Mechanical characteristics research on flexible hose under the influence of internal and external fluid. Journal of Water Resources and Water Engineering, 2018, 29(3): 182-193 (in Chinese) doi: 10.11705/j.issn.1672-643X.2018.03.31 [31] Zhu HJ, Gao Y, Zhao H. Experimental investigation on the flow-induced vibration of a free-hanging flexible riser by internal unstable hydrodynamic slug flow. Ocean Engineering, 2018, 164: 488-507 doi: 10.1016/j.oceaneng.2018.06.071
[32] Meng S, Chen Y, Che C. Slug flow's intermittent feature affects VIV responses of flexible marine risers. Ocean Engineering, 2020, 205: 106883 doi: 10.1016/j.oceaneng.2019.106883
[33] Li FQ, An C, Duan ML, et al. Combined damping model for dynamics and stability of a pipe conveying two-phase flow. Ocean Engineering, 2020, 195: 106683 doi: 10.1016/j.oceaneng.2019.106683
[34] Passano E, Larsen CM, Lie H, et al. VIVANA−Theory Manual Version 4.8. Trondheim, Norway, 2016
[35] Gopalkrishnan R. Vortex-induced forces on oscillating bluff cylinders. [PhD Thesis]. Cambridge, MA: Massachusetts Institute of Technology, 1993
[36] 宋磊建. 涡激振动状态下的柔性立管水动力特性研究. [博士论文]. 上海: 上海交通大学, 2016 (Song Leijian. Investigation on the hydrodynamics of a flexible riser under vortex-induced vibration. [PhD Thesis]. Shanghai: Shanghai Jiao Tong University, 2016 (in Chinese))
[37] 唐友刚. 高等结构动力学. 天津: 天津大学出版社, 2000 (Tang Yougang. Advanced Structural Dynamics. Tianjin: Tianjin University Press, 2000 (in Chinese))
[38] Blevins RD, Coughran CS. Experimental investigation of vortex-induced vibration in one and two dimensions with variable mass, damping, and Reynolds number. Journal of Fluids Engineering, 2009, 131(10): 101202
下载:
计量
- 文章访问数: 997
- HTML全文浏览量: 480
- PDF下载量: 134
