fluent周期边界_FLUENT流固耦合柱体结构涡激振动仿真案例解析

作者 | 陈东阳博士 仿真秀科普作者

一、柱体结构涡激振动定义

对于海洋工程、风工程上普遍采用的圆柱形断面结构物，流体绕过柱体时会产生交替发放的泻涡，这种交替发放的泻涡又会在柱体上生成顺流向及横流向周期性变化的脉动压力。

如果此时柱体是弹性支撑的，或者柔性管体允许发生弹性变形，那么脉动流体力将引发柱体(管体)的周期性振动，这种规律性的柱状体振动反过来又会改变其尾流的泻涡发放形态。这种流体-结构物相互作用的问题被称作"涡激振动"(Vortex-Induced Vibration :VIV)，它是典型的一种流固耦合振动现象。

二、柱体结构涡激振动动力学模型及仿真方法简介

进行二维弹性支撑柱体的数值模拟是研究海洋工程和风工程中柱体结构涡激振动现象和机理的重要手段。

根据牛顿第二定律，2-DOF弹性支撑的柱体运动的控制方程可以写为：

(1)

式中， m为圆柱体的质量， c为结构阻尼系数， k为结构刚度系数。

式(1)又可以写为：

(2)

式(2)中，柱体固有频率

，

阻尼比

对于涡激振动的数值模拟大多以柱体的横向振动研究为主，同时考虑横向和来流向的耦合振动研究相对较少。其主要原因在于，对柱体结构涡激振动数值模拟，必须保证柱体周围网格质量非常好，才能有效预测涡激振动响应。

一般使用网格质量高的结构化网格，但是当柱体结构发生较大振动位移时，周围流场网格会发生畸变，甚至产生负网格，导致计算失败。如果再考虑流向耦合振动，计算难度将非常大，计算成功率也将明显降低。

如果采用非结构化网格，且采用网格重构技术，可以吸收柱体较大的振动位移，但是非结构化网格质量相比结构化网格质量较差，且采用非结构化网格必定需要大大增加网格量，从而大大增加了计算时间。

此外，如果对于涡激振动抑制装置设计研究，也就是在柱体结构表面形状变复杂的情况下，网格划分难度和计算难度也将大大增加。因此，寻求一种既可以保证网格质量，又能不大幅度增加网格数量，且可以避免网格畸变或者负网格问题的方法十分重要。

基于CFD商业软件FLUENT和结构动力学原理，通过用户自定义函数(UDF)及嵌套网格技术，建立了2-DOF弹性支撑柱体结构VIV数值模型。非定常不可压缩流体RANS方程为：

(4)

(5)

式(5)中，