深水钻井隔水管柱自激振动响应分析.pdfVIP

51卷 增刊 2 t{J 困 造 船 Vo1．51 Special2 2Ol0年 12月 SHJPBUJLDINGOFCHJNA Dec．2O】0 文章编号：1000—4882(2010)S2．028—06 深水钻井隔水管柱 自激振动响应分析 孙友义 1,2 陈国明 ，鞠少栋 (1．中国石油大学 海洋油气装备与安全技术研兖中心，东营 257061； 2．中海油研究总院技术研发中心 博士后工作站，北京 100027) 捅 姜 基于RNGk一￡湍流模型与动网格技术，针对隔水管在海洋环境中的实际雷诺数范围，对管柱涡激振动进 行数值模拟。数值模拟再现了圆柱响应的三个分支：初始分支、高幅分支与低幅分支，与经典实验吻合 良 好，同时清晰地捕捉到了漩涡泄放的2s模式与2P模式。研究发现在解锁 区内存在一个振幅显著低于低幅分 支且大小稳定的响应分支，其响应频率仍锁定在固有频率上，称之为 超“低幅分支”。在低幅分支与超低 幅分支中，发现升力系数中存在高阶频率；升力系数中高阶频率的出现是导致圆柱响应 由高幅分支向低幅 分支转变的主要原因。 关 键 词：隔水管；涡激振动；数值模拟；动网格；响应分支 0 引 言 关于圆柱体涡激振动的数值模拟，国内外学者已做 了较多的研究工作。这些工作多局限于低霄诺 数，同时对受迫振动的关注远多于 自激振动。实际的海洋细长结构物所处雷诺数范围在lxl04~lx10 之问，且其振动形式是 自激的。 管受迫振动的预测结果对于 自激振动预测具有借鉴意义，但二者从 本质_卜仍有很大不同，如受迫振动中观察到的某些漩涡泄放形式(如P+S模式)不会出现在 自激振动中。 由于VIV涉及到许多复杂的物理现象如湍流、流 耦合，加上数值计算内在的不稳定性及精度 问题， 使得VIV的数值模拟异常困难，少有计算结果能够精确再现经典的实验结果。目前模型实验难以模拟 高雷诺数环境，对高雷诺数下的圆柱VIV进行数值模拟就变得尤为重要。本文对 维圆柱的自激振动 进行数值模拟，针对隔水管在海洋环境 中的真实雷诺数火小进行仿真，较好地再现 了Khalak、 Williamson以及Govardhan等[1-5]关于低质量阻尼圆柱体涡激振动的实验结果。 1 分析方法与分析模型 1．1流体动力学控制方程 对于二维平面势流问题 ，忽略流体的呵压缩性，以及流动过程中的热量传递与气泡产生，则 质繁守恒方程(连续方程)与动量守恒方程(N—S方程)可以表示如下： V ·v=0 (1) 8t+(V。)‘v=VLVp十 (2) 。 金项 目：嘲家高技术研究发展计划(“863”汁划)资助项 目(No．2006AA09A106—4) 孙友义，等：深水钻 t隔水管柱 Irj激振动响应分析 29 式中， 为速度矢量：v表示散度；，为体力矢量； 为流体密度；P为流体微元体上的压力； 为 流体动力粘性系数。连续方程和N．S方程是流场控制基本方程，结合具体边界条件和初始条件使得方 程封闭，可以求解得到流场分布。 实际海洋环境下，隔水管所在海水流动往往处于湍流状态。湍流状态下，流动微粒随机性较强， 流动呈混乱状态，这时即使边界条件保持不变，流动也不是稳定的。一般认为，无论湍流运动多么复 杂，非稳态的连续方程StlN．S方程对于湍流的瞬时运动仍然是适用的。工程 中广为采刖的方法是对瞬 态N—S方程做时间平均处理，同时补充反映湍流特性的其他方程，如湍动能方程和湍动耗散率方程等。 Reynolds平均法是 目前使用最为广泛的湍流数值模拟方法。笔者采用RNG k一 模型模拟湍流宏观效 应，结合非平衡壁面函数能够很好地处