Re数5000时圆柱受迫振动的数值的研究.pdfVIP

50卷增刊 中 国 造 船 V01．50Special 2009f#11月 SHIPBUILDlNGOFCHINA Nov．2009 文章编号：1000-4882(2009)S．0186．06 Re数5 000时圆柱受迫振动的数值研究 赵静1，董国海2，吕林2 (I．中国石油集团工程技术研究院。天津300451； 2．大连理工大学海岸和近海工程国家重点实验室，大连116024) 摘要 在Navier-Stokes方程和Ic啪湍流模型的基础上，利用流线迎风有限元方法对Re一5000时的圆柱受迫振动 问题开展了数值模拟研究．从计算的结果来看，本文的数值模型成功模拟了Re=5000条件下圆柱发生受迫振 动时尾迹区内的2S。2P和P+S尾漉模式，同时升力及拖曳力的计算结果与已有文献符合较好，为今后进一步 开展更高雷诺数条件下的受迫振动研究工作奠定了基础． 关 键词：受迫振动；湍流模型；流线迎风有限元；N．S方程 0引言 现代海洋开发中，海洋石油的开发最为引人注目。不论采用何种形式的结构物开发海洋油田，都 不可避免地要使用到立管，它们是海洋基础结构的关键组成部分。因此，立管等细长柔性结构的涡激 Induced 振动(VortexVibration，VlV)问题受到了越来越多的重视。 通常将圆柱体在稳定流中的运动分为受迫振荡和自激振荡两种基本类型，二者都是用来研究涡激 振动的简化模型。对于深海隔水管而言，如果从几千米长的总体上考虑，这些柔性立管系统的振动形 式可视为自由振动，但就某一局部微段而言，则是在端部位移约束作用下的受迫振动。相比而言，受 迫振荡较自激振荡能更多地反映出涡与柱体振动之间的相互作用以及能量转移的信息，因此本文主要 针对圆柱体的受迫振动问题进行数值研究。 受迫振动的数值模拟研究中，较常用是离散涡法和有限元法。由于和离散涡方法相比，有限元法 的适用性更强，因此更多学者选择采用有限元法11-51开展研究工作。相对而言，国内对受迫振动的数 值研究较少，潘志远【6】比较了圆柱体涡激振动实验和受追振动实验各自的特点和二者的关系，并从实 验现象和数值模拟的结果分析入手，剖析了两类机理实验研究的异同。梁亮文和万德成17J利用 FLUENT软件对Re,=200时均匀流中圆柱的受迫振动问题进行了数值模拟。 本文将建立和发展基于流线迎风有限元的稳定流固耦合运动数值分析模型。通过求解二维雷诺平 动时升力、拖曳力和尾流形态的变化特征，通过与数值结果的对比来验证模型的可靠性，为今后开展 更高雷诺数下的圆柱受迫振动数值研究奠定基础。 1 数值模型 下的无量纲形式， 尝=o (1) a譬， 50卷增刊 赵静。等：Re数5000时嘲柱受迫振动的数值研究 187 c2， --鲁+(／41--U?，考=一苦(p+；七)丢(去考+2_岛] 的流体运动速度分量，Ⅳ7表示．，方向的网格结点运动速度，P为压力，k为湍动能，，为时间， =《锄，／教，+锄，／魏)／2为平均应变张量，h为湍流粘性系数，Re=uod／v表示雷诺数，坳为均匀来流 流速，d为圆柱直径，v为流体的运动学粘性系数。 k-to湍流模型对基本控制方程进行封闭， 豢+(吩一u，．，、qak=考[(去+d‘_]鸶]+B一∥’破 (3) 害+(吩一哆)考=考[(去+吖)考]+警只一肋2 (4) 式(3)、(4)中，只=2咋(钆，％)，表示湍动能k的产生率，其余各模型常数的取值如表l所示。