中文摘要：

针对结构涡激振动数值模拟问题,以开源流体动力学计算软件OpenFOAM 3.0为平台建立结构涡激振动数值模拟方法,即将结构涡激振动问题简化为单自由度振子模型,采用Newmark-β法进行结构振动方程求解,应用OpenFOAM动网格求解器进行动网格计算与更新。以宽高比为5的矩形断面为例,首先采用3种不同的雷诺平均（RANS）湍流模型（即SST k-ω,k-ε和k-ω模型）进行风攻角为a=0°,2°,4°,6°时静止绕流计算,以检验不同湍流模型的计算精度;然后对矩形断面在0°风攻角下竖向涡振响应进行了数值模拟,并与已有文献和试验结果进行比较。研究表明：3种湍流模型的计算结果总体较为接近,且与试验结果较为吻合,其中SST k-ω湍流模型的计算结果与试验结果吻合相对最好;对于升力系数和升力矩系数,3种湍流模型的计算结果都与试验结果存在一定的差异;总体而言采用OpenFOAM的SST k-ω湍流模型所得静止矩形断面绕流计算结果与试验结果吻合较好。矩形断面涡激振动响应数值模拟结果与试验结果相比锁定风速区间有一定前移,且最大涡振振幅较试验结果略偏大。

英文摘要：

For the numerical simulation of VIV of structures,the numerical simulation method of VIV of structures based on the open source software OpenFOAM 3. 0 is established,i. e.,to simplify the VIV of structures into a single degree of freedom oscillator model,to solve the structure vibration equation by Newmark-beta method,and to calculate and update the dynamic grid with OpenFOAM dynamic mesh solver.Taking a rectangular section with aspect ratio of 5 for example,the flows around the stationary rectangular section at attack angles of 0°,2°,4° and 6° respectively are calculated by 3 Reynolds averaged（ RANS）turbulence models（ SST k-ω model,k-ε model and k-ω model） to check the simulation accuracy of different turbulence models. Then,the numerical simulation of vertical VIV response of the rectangular section with wind attack angles of 0° is conducted using dynamic mesh technology,and the numerical result is compared with the existing literature results and experimental results. The result shows that（ 1） the calculation results of the 3 turbulence models are generally close to each other and are consistent with the experimental result,the calculation result of SST k-ω turbulence model agrees well with the experimental result relatively;（ 2） for the lift and lift moment coefficients of the rectangular section,the calculation results of the 3 turbulence models are different from the experimental result;（ 3） generally,the numerical result of flow around the stationary rectangular section with SST k-ω turbulence model in OpenFOAM agrees well with the experimentalresult;（ 4） the general trend of numerical simulation result of VIV response of the rectangular section agrees well with the experimental result,however,the numerical simulation result of the VIV response of the rectangular section in lock-in wind velocity region is less than the experimental result,and the numerical simulation result of VIV maximal amplitude is slightly larger than the experimental result.