2. 中国石油物资公司
2. CNPC Materials Company
0 引 言
跨接管是水下生产系统的重要组成部分之一,起着连接水下采油树、管汇及处理中心的作用[1]。根据刚柔性可以分为刚性跨接管和柔性跨接管[2],根据安装方式可以分为垂直式跨接管和水平式跨接管[3]。常见的垂直式刚性跨接管有U形和M形2种。
在水下油气生产过程中,海水流过跨接管表面时,会在尾流区产生周期性变化的脉动涡漩。若涡漩的脱落频率与跨接管的固有频率相同或相近,就可能产生共振现象,这就是涡激振动[4]。涡激振动可能会导致管道的疲劳损伤甚至失效。
近年来,国内外对海洋立管和隔水管等长跨度管线的涡激振动分析较多[5, 6, 7],但对跨接管的研究鲜有报道。笔者以垂直式M形跨接管为研究对象,基于ABAQUS和Fluent软件,对水下跨接管线发生涡激振动的可能性进行分析,以期为水下油气安全生产及输送提供指导。
1 几何模型及模态分布M形跨接管几何模型见图1。管道内径125 mm,外径D=172 mm,特征长度11 m,高度3 m。
假设跨接管安装在水下1 000 m处,其对应参数为[8]:壁厚23.5 mm,弹性模量210 GPa,密度7 850 kg/m3,海水密度ρ1=1 025 kg/m3、动力黏度1.2 mPa·s、流速0.1~0.4 m/s。当海水流速为0.4 m/s时,雷诺数为:
Re小于3×105,处于典型的亚临界雷诺数区域。
用Solidworks软件建立跨接管的三维模型,并将模型导入ABAQUS软件,定义其材料属性。由于跨接管两端与接头相连,故在两端施加固定约束。选用C3D8R六面体缩减积分单元为模型划分网格,结果如图2所示。整个分析过程基于线性摄动分析步,选用Lanczos方法进行迭代[9]。
基于ABAQUS软件的跨接管前4阶振型(变形放大1 000倍)如图3所示。
提取跨接管前4阶固有频率,如表1所示。
为得到海水流过跨接管表面时产生的涡脱频率,使用Fluent软件对三维跨接管进行绕流分析,在计算过程中选用RNG κ-ε湍流模型,对近壁区域选用增强型壁面处理[10, 11]。
2.1 计算域模型流体域的计算模型见图4,模型长15 m、宽5 m、高7 m。该计算域的网格划分结果见图5。为获取边界层区域的流动特征,对近壁区域的网格进行加密,网格数约820万。计算使用的边界条件为速度入口和压力出口,计算域上、下、左、右4个面均设置为对称边界,跨接管与海水接触的外表面粗糙度设置为0的无滑移壁面。由图6可见,海流在跨接管迎流面出现驻点及低速区;在跨接管上、下两侧速度升高、压力减小;在跨接管背流面速度升高、压力增大,发生了明显的流动分离现象。在尾流区域,漩涡呈周期性的上、下交替脱落。在达到管段50倍直径距离的区域,漩涡才逐渐消失,符合亚临界雷诺数下圆柱绕流的流态变化规律[9]。
2.2 计算结果计算得到不同海水来流速度下跨接管中心截面对应的绕流速度场如图6所示。
在计算过程中,监控跨接管表面的升力系数(Cl)和阻力系数(Cd)随时间的变化。以海流速度0.4 m/s工况为例,得到升力系数和阻力系数变化曲线,如图7和图8所示。
对升力系数曲线做快速傅里叶变换,得到能量谱密度与频率的关系曲线,如图9所示。
由图9可知,当频率为0.615 Hz时对应的能量谱密度值最大,即在海流速度0.4 m/s的情况下,涡脱频率为0.615 Hz。采用同样的方法,对其他流速下的绕流分析进行快速傅里叶变换,得到对应的涡脱频率,如表2所示。
对表2中数据进行处理,得到海水流速和涡脱频率关系曲线,见图10。从图可见,海水流速在0.1~0.4 m/s范围内,涡脱频率随着流速的升高而增大,最大频率达到0.615 Hz。结合ABAQUS软件计算出的跨接管固有频率,其1阶振型所对应的固有频率为3.045 4 Hz,其他振型对应的固有频率均高于此。显然,跨接管的固有频率超出了涡脱频率的范围,水下跨接管理论上不会因涡脱而产生共振现象。
3 结 论(1)通过ABAQUS有限元软件得到了水下跨接管前4阶模态下的振型及固有频率,其中1阶模态所对应的固有频率为3.045 4 Hz。
(2)当海水流速在0.1~0.4 m/s范围内时,跨接管尾流区的涡脱频率与海水流速呈线性关系,其固有频率超出了涡脱频率的范围,不会因涡脱而产生共振现象。
[1] | 王玮,孙丽萍,白勇.水下油气生产系统[J].中国海洋平台, 2009,24(6):41-45. |
[2] | 甄兴伟,黄一,张崎,等.基于新概念水中张力生产平台的柔性跨接管对刚性立管系统整体特性的参数影响研究[J].船舶力学,2014(6):711-723. |
[3] | 何同,李婷婷,段梦兰,等.深水刚性跨接管设计的主要影响因素分析[J].中国海洋平台,2012,27(4):50-56. |
[4] | 徐枫,欧进萍.方柱非定常绕流与涡激振动的数值模拟[J].东南大学学报:自然科学版,2005,35(增刊1):35-39. |
[5] | 蔡杰.海洋立管流固耦合特性数值模拟研究[D].上海:上海交通大学,2010. |
[6] | 赵鹏良,王嘉松,蒋世全,等.海洋立管涡激振动的流固耦合模拟计算[J].海洋技术,2010,29(3):73-77. |
[7] | 蔡杰,尤云祥,李魏,等.均匀来流中大长径比深海立管涡激振动特性[J].水动力学研究与进展,2010,25(1):50-58. |
[8] | 任大朋,黄一,刘刚.一种基于ANSYS和FLUENT的海洋立管的涡激响应分析方法[J].中国海洋平台,2007,22(4):32-36. |
[9] | 石亦平,周玉蓉.ABAQUS有限元分析实例详解[M].北京:机械工业出版社,2006:283-294. |
[10] | 詹昊,李万平,方秦汉,等.不同雷诺数下圆柱绕流仿真计算[J].武汉理工大学学报,2008,30(12):129-132. |
[11] | 王亚玲,刘应中,缪国平.圆柱绕流的三维数值模拟[J].上海交通大学学报,2001,35(10):1464-1469.,35(增刊1):35-39 |