随着海洋油气资源的开发逐步向深海推进,各种适应深海水域环境的采油平台纷纷涌现。不论采用何种形式的结构物开发海洋油气资源,都不可避免地要使用圆柱形结构,例如支撑立管、输油管线等。它们连接了海底矿藏与海面的作业平台,是海洋基础结构的关键组成部分,是进行深水石油天然气开采不可或缺的设备。一直以来,立管问题都是海洋工程的重要研究内容,通常将立管的运动简化为圆柱体的受迫振动和自由振动两种基本类型,二者都是用来研究涡激振动的简化模型。相比而言,受迫振动较自激振动能更多地反映出涡与柱体振动之间的相互作用以及能量转移的信息,因此本文主要针对圆柱体的受迫振动进行数值研究。本文在Navier-Stokes方程和k-ω湍流模型的基础上,利用流线迎风有限元方法结合ALE动网格技术,对稳定流中低雷诺数(Re=100,Re=200)和亚临界雷诺数(Re=10000)下的圆柱受迫振动问题开展了数值模拟研究。着重从升阻力历时曲线,升阻力随振动频率的变化关系,受迫振动锁定区间的划分以及尾流场涡流形态几个方面进行分析。计算中,首先在圆柱振幅一定的条件下,通过固定圆柱运动的振幅,得到了圆柱受力随振动频率的变化关系:研究表明,锁定发生前,升阻力值随振动频率的增大而增大;当振动频率接近自然泄涡频率时,升阻力值都将达到最大值,此时振动发生锁定;随着振动频率的继续增加,振动逐渐远离锁定区,升阻力值在经历短暂减小过程后再次增大。另外,结合升阻力变化关系,通过FFT频谱分析,划分了各个振幅下的锁定区间,并和已有实验和数值结果比较良好。通过对尾流场泄涡形式的分析,分别得到了2S,2P,P+S三种典型受迫振动尾流模式。总体来讲,本文采用的二维数值模型能够基本正确地反映出圆柱发生受迫振动时的涡激振动特性以及有关的受力变化趋势,这为今后进一步开展三维数值分析工作奠定了重要的基础。