基于流固耦合的水翼涡振特性分析

2022-05-09 02:42刘新月马如相杜雷雨

兵器装备工程学报 2022年4期

刘新月，洪亮，马如相，南栩，杜雷雨

(南京理工大学能源与动力工程学院，南京 210094)

1 引言

近年来，随着船舶行业的蓬勃发展，舰艇的航速大大提升，舰艇水下结构的水弹性振动得到了造船界的重视。若舰艇的航速大，舵翼的结构强度不高，舵翼在水下运动时有可能会产生颤振，造成结构破坏，威胁航行安全。而在舰艇航速较低时，会产生一种容易被忽视的振动现象“涡振”，涡振即结构在流场中产生涡的频率与结构固有频率一致时发生的共振现象。涡振虽然是限幅运动，但会造成结构疲劳，长此以往，造成结构破坏。对于潜艇来说，发生涡振时会产生显著的噪声，破坏隐身性能。

针对水翼的振动现象，国内外的众多学者展开了研究。Chen等[1]采用SSTk-ω湍流模型对平行槽水翼叶片进行数值分析，研究了水翼位置和间距对涡产生的抑制效果和流动特性的影响。结果表明，叶片间距越小，越容易产生涡，并根据分析结果选取翼型。Huang等[2]选用NACA0012翼型进行了试验，着重探究了翼梢对涡激振动的影响，发现翼梢处容易产生涡量大的涡。Huerre等[3]以及Oertel[4]通过试验分别给出了流动不稳定性时涡发放特点和钝体尾流形态的概述。Chae[5-6]通过对二自由度下二维水翼的研究，讨论了水动力阻尼、非线性流固耦合响应等对水翼振动特性的影响。余志兴等[7-8]采用N-S方程以弹簧系统代替实际固体变形产生恢复力，发现要构成翼型涡振，攻角要足够大。侯磊等[9]采用离散涡方法，分析高雷诺数下来流速度及攻角对水翼涡发放频率影响规律。刘胡涛等[10]对二维水翼结构进行流固耦合运动分析，探讨了初始攻角、刚心位置以及来流速度对水翼振动的影响。郭春雨等[11]采用DEM模型对三维水翼梢涡处流场和气核运动特性的问题进行数值分析，为梢涡空化初生的机理研究建立基础。贾文隽等[12]分别采用传统CFD方法和格子玻尔兹曼流固耦合方法对平直机翼和仿生机翼进行数值模拟，并通过气结构响应对比找出了仿生机翼的优点及不足。赵宇等[13]采用数值方法对绕水翼的非定常空化流动进行了流动计算分析，结果表明空化现象对水气两相转换起到了促进作用，空化现象影响涡结构空间分布。陈衍茂等[14-15]采用同伦分析法、增量谐波法等研究含有立方非线性的二自由度机翼的极限环颤振，并对极限环的幅值和频率进行求解。

综上所述，对于水翼涡振特性的研究主要集中在单一翼型涡振机理等方面的研究，对于相似翼型的涡振特性差异方面存在空白。本文首先将所研究的翼型设置为刚体，探究其泄涡特性，然后通过双向流固耦合的方法，研究其涡振特性，并对几种相似翼型进行对比分析，为舰艇舵翼选型提供参考。

2 控制方程及流固耦合求解方法

2.1 流体控制方程

由于舵翼为流线型，后缘产生的涡多为细涡，雷诺平均湍流模型难以捕捉，故采用大涡模拟。大涡模拟综合采用了直接模拟和雷诺平均模拟的思想与方法。将湍流中的脉动划分尺度，大尺度湍流直接采用数值求解。而对小于计算网格尺度的湍流模型建立模型求解，其控制方程为：

(1)

(2)

2.2 结构振动控制方程

在本研究中，将三维水翼简化为二维截面，并引入弹簧阻尼系统，将水翼的振动简化为竖向和绕旋转轴扭转的两自由度的运动。

(3)

(4)

2.3 双向流固耦合求解方法

双向流固耦合即流场计算结果压力载荷通过流固耦合界面传递给结构场，同时结构场的位移变形结果通过流固耦合界面传递给流场，双向流固耦合分析流固耦合交界面处的数据传递是双向的。流场和结构场之间是相互影响的，考虑结构的变形对流场的影响。

本文在计算水翼涡激振动响应的每一个时间步长内，先求解流体控制方程，得到作用在水翼上的水动力FL(t)和力矩Mα(t)，将得到的水动力和力矩代入结构振动方程式。采用Newmark-β法，求解结构的振动响应，利用动网格技术，将得到的速度赋予结构及结构壁面附近的网格，使其随结构一起做刚体运动，以实现水翼与流场的双向耦合。

3 物理模型及网格

计算选取了对称翼型NACA0009、NACA0012和NACA0015，其中重点介绍了NACA0012翼型。NACA翼型末尾两位数字代表了相对厚度的大小，NACA0015翼型最厚。水翼弦长为0.1 m长，旋转轴设置在距前缘1/3弦长处，来流攻角为0°。

采用的网格策略为“刚性边界层运动区域+动网格变形区域+外流场”的划分策略，翼型表面附有贴体网格，计算Y+值，保证近壁面处网格的质量。翼型外圆形区域为随动网格，随翼型一起运动，以保证翼型计算的收敛性及精度。圆环区域为动网格区采用三角形网格，以保证具有较好的运动性。最外围为流场网格。图1为所使用网格和局部放大图。

边界条件设置：左侧弧形边界以及上下边界设置为速度入口(velocity-inlet)，右侧边界设置为压力出口(pressure-outlet)，翼型采用无滑移壁面边界(wall)。

图1 网格示意图

湍流模型选用大涡模拟。计算算法采用“Coupled”，该算法同时求解流体的连续方程、动量方程和能量方程，适合求解本文关于不可压缩及等温流动的耦合问题。采用二阶迎风格式，以提高计算精度。水翼由于其结构特点，泄涡频率较高，时间步的选取影响计算精度，本文经过时间步无关性验证，选取时间步长为0.000 5 s，限于篇幅不再赘述。

本文二维水翼的结构参数取值来自于工程实例。水翼弦长为0.1 m，展长为0.2 m，材料为POM(聚甲醛树脂)，采用一端固支的方法，基于ANSYS平台对其提取质量、模态等结构参数。竖弯刚度由式(5)计算得出:

(5)

式中:f为水翼一阶竖弯固有频率，m为水翼单位展长质量。

在对水翼进行模态求解时发现：水翼的扭转固有频率要远大于竖弯固有频率。竖弯方向的涡振在低速下更容易发生，本文只分析了竖弯方向单自由度的一阶涡振现象.本文使用的翼型结构参数取值如表1。

表1 水翼结构参数Table 1 Structural parameters of hydrofoil

4 刚性水翼泄涡特性分析

水翼设置为刚体，分别计算了NACA0012水翼在流速0.5 m/s、1.0 m/s、1.5 m/s、2.0 m/s下的升力系数。图2、图3分别为水翼在流速0.5 m/s、1.0 m/s下的升力系数随时间变化曲线，可以看出，升力系数变化的范围不大，在计算稳定后，升力系数曲线呈规律的正弦波形，反映了水翼表面受到了周期性变化的垂向载荷。随着流速的增加，水翼受到的垂向载荷频率变大。

图2 0.5 m/s时升力系数曲线

图3 1.0 m/s时升力系数曲线

对4种流速下的水翼的升力曲线进行快速傅里叶变换(FFT)，得到水翼受到的垂向载荷的频率，即水翼的泄涡频率，表2为提取到的各个流速下的泄涡频率。图4表示流速与泄涡频率的关系，可直观看出泄涡频率与流速呈现线性关系。

表2 泄涡频率Table 1 Vortex shedding frequency

图4 流速与泄涡频率关系曲线

以上为NACA0012翼型的分析，采用同样的方法，计算了NACA0009以及NACA0015翼型在同流速下的泄涡频率，如图5所示。

图5 各翼型泄涡频率曲线

从图5可以看出，本文选取的3种对称翼型的泄涡频率都随着流速的增加而增加。在同流速的比较下，水翼的厚度对泄涡频率产生影响，厚度越大的水翼，泄涡频率越低。

5 弹性水翼涡振特性分析

前文对刚性水翼型的泄涡特性进行了分析，只考虑了流场对水翼的作用，得到的规律较为简单。本节将水翼设置为弹性水翼，使用双向流固耦合的方法，对水翼的涡振特性进行分析。

在计算刚性水翼NACA0012时，水翼在1 m/s时泄涡频率达到了50.37 Hz，与该翼型的自振频率55.15 Hz十分接近，现在1 m/s左右的区间进行细化，提取水翼在不同流速下的竖向位移图。

图6、图7和图8分别给出了水翼在0.5 m/s、1.0 m/s、1.5 m/s时的竖向振动时域图。在低流速时水翼振动较为规律，0.5 m/s时水翼振幅达到最大数值后振幅略有缩小，而后以稳定振幅振动。1 m/s时水翼振幅呈现“喇叭形”，推测在该流速下发生涡振。1.5 m/s时水翼振幅不规则，反映流场在较高流速下湍流度变大，水翼受力情况更复杂。

图6 0.5 m/s竖向位移曲线

图7 1.0 m/s竖向位移曲线

图8 1.5 m/s竖向位移曲线

图9—图11为水翼的振动频谱曲线。在0.5 m/s的流速下，水翼的振动频率较低，与前文计算刚性水翼在该流速下的泄涡频率23.33 Hz相近。在1 m/s时，水翼振动频率与自振频率相近，符合发生涡振的特征。随着流速增加，频谱图出现双峰特征，左侧峰集中在水翼的自振频率附近，呈现明显宽谱特性，右侧峰为水翼的泄涡频率，较为集中。