基于频率响应法的海底悬跨管道涡激振动预报.docxVIP

基于频率响应法的海底悬跨管道涡激振动预报 原油长距离输送是海洋油气资源开发的重要组成部分。由于管道支墩或海床表面凹凸不平, 或管道周围海床的冲刷等因素的影响, 平铺在海床上的管道容易出现悬空段。在一定条件下, 水流流经管线时, 管道的尾流区会产生旋涡。当旋涡的泄放频率接近于管道的自振频率时, 将出现旋涡泄放频率在一个较大的流速范围内接近于管道自振频率的 “频率锁定”现象, 并发生大振幅振动, 即“涡激共振”, 这种大振幅振动是造成悬空管道疲劳破坏的主要原因。 在过去的三十多年里, 很多国家的学者对海洋工程细长结构物的涡激振动现象进行了大量的研究, 提出了一系列涡激振动预报模型, 其中大部分是针对深水立管的涡激振动问题。从水动力的角度来看, 这些模型主要可以分为三类:一类是基于切片理论的CFD方法, 如Norsk Hydro、USP等, 这类模型大都将二维N-S方程的求解与梁模型结合在一起, 通过求解Navier-Stokes方程得到作用在圆柱上的力, 然后将力反馈到圆柱上求得结构的动力响应。从与实验结果的比较来看, 计算结果大部分低估了结构在锁定条件下的响应;第二类是基于圆柱受迫振动实验数据的经验模型, 这类模型的基本假定是涡激振动发生在一个或有限个离散频率上, 目前这类模型用的较多的主要有SHEAR7、VIVA、VIVANA等, Larsen对这一类模型进行了详细的讨论;还有一类是介于上述两者之间的尾流振子模型, 通过相互独立的方式分别建立圆柱振子运动方程和流体振子运动方程, 然后利用它们共同预报流体-弹性系统的动力响应, 这些模型通常不进行任何流场分析, 且它们的值最大限度只能解释和模拟一些实验结果。尾流振子模型一般只用于横流向振动。 海底悬跨管道作为一种常见的海洋细长结构物, 其涡激振动问题与深水立管有较大的不同。管线与海床之间的间隙比沿管线轴线的连续变化使海底悬跨管道涡激振动问题的求解变得相当复杂。Halse采用基于切片理论的CFD方法提出了在均匀流和剪切流中海底悬跨管道的涡激振动分析方法。Hansen et al分析了管沟附近海底悬跨管道的动力特性。Choi考察了轴力和边界条件对管跨自振频率、容许管跨长度的影响, 其响应计算方法采用了DNV规范G-14的方法。Furnes et al对海底悬跨管线建立了一个时域模型, 对顺流向 (IL) 和横流向 (CF) 两个方向振动耦合的力学机理进行了研究, 模型假设IL和CF挠曲是通过任意时刻的轴向张力耦合的, 张力沿空间均匀分布, 可通过计算位移引起的伸长得到;采用的水动力表达式足够光滑以便处理随时间和空间连续变化的任意水流。傅强和郭海燕在考虑内流影响的情况下, 采用复模态分析方法考察了流体阻尼对海底悬跨管道自振频率的影响。最新的DNV规范RP-F105提供了几个响应模型 (response models) 用于计算VIV响应;响应模型中, 结构的响应幅值是折减速度的函数;在计算均匀流中的VIV时 (不考虑涡流湍动和波的影响, 忽略海床近壁和管沟的影响) , IL和CF方向的响应曲线分别决定于稳定性参数和频率比。Drago et al对冲刷引起的海底悬跨管线在波浪和水流作用下的疲劳特性进行了分析。Maincon et al提出了一种基于实验数据用于计算线性结构涡激振动载荷和动力响应的逆有限元方法 (an inverse finite element method) 。 CFD方法对计算机的计算能力要求特别高, 尚不成熟, 目前实际工程仍然基于经验模型。基于由挪威海洋技术研究所 (MARINTEK) 和挪威科技大学 (NTNU) 共同研发的VIVANA模型, 编制了计算程序, 探讨了不同流速、不同悬跨长度和不同外径的海底悬跨管道的涡激振动特征。 1 动力响应频率 采用频率响应法计算海底悬跨管道涡激振动响应, 其动力平衡方程为 Μ¨r(t)+C˙r(t)+Κr(t)=R(t)(1)Mr¨(t)+Cr˙(t)+Kr(t)=R(t)(1) 假定结构各点的响应都是简谐的, 各点的响应之间存在一个相位差, 载荷和响应之间存在一个相位相关。这样结构的位移响应r(t) 和载荷R(t) 可以分别用一个随时间做简谐变化的复位移矢量和复载荷矢量来表示 r(t)=xeiωt(2)R(t)=Xeiωt(3)r(t)=xeiωt(2)R(t)=Xeiωt(3) 动力平衡方程式 (1) 可表示为 [-ω2(ΜS+ΜΗ)+iω(CS+CΗ)+Κ]x=X(4)[?ω2(MS+MH)+iω(CS+CH)+K]x=X(4) 式中:位移矢量x和激励力矢量X均为复矢量;MS为结构质量矩阵;MH为附加质量矩阵, 是ω和水流速度的函数;K为刚度矩阵, 包括材料刚度和几何刚度;CS为结构阻尼矩阵;CH为水动力阻尼矩阵,