均匀流中印体振荡的频率%27锁定%27现象.pdfVIP

第十届全国计算流体力学会议论文集 均匀流中钝体振荡的频率“锁定现象 刘松符松 清华大学工程力学系，北京100084 摘要为研究流致振菇问题中柱体振动频率与近尾迹周期性的频率“锁定。现象，本文连取 具有代表性的腰柱绕流作为研究对象．用有限体积法对平行于均匀来流方向强迫振荡的圆柱 绕流问压进行了二维数值模拟．本文给出了频率锁定发生的杀件和类型。入流雷诺数选取 Re=200、855等几种亚临界冒诺数情况．计算条件采用了较小的折台振幅和较宽的折合频率 范围。对应于不同的涡脱落模态和频率锁定的发生．流场结构、升阻力系数与涡脱落SUov．ha】 数有强烈变化． 关键词：流致振动，圆柱绕流，涡脱落，频率锁定 0引言 1880年，Rayleigh在观察弦线在空气中的振动时，注意到弦的振动不是沿着风向， 而主要发生在与风垂直的方向，从此流致振荡问题为人们所关注并研究。对均匀不变的 来流，由于钝物体的尾迹流动伴随着复杂的分离和旋涡脱落现象，诱导的旋涡脱落可能 诱发交变的流场，导致物体大幅度的振动。这是流致或涡致动态激振的基本形式。 流致振动现象在工程上较为常见，高层建筑、烟囱、近海工程结构以及海底管线、 架空电缆中．经常会遇到旋涡脱落诱发的流体动力载荷和结构振动问题，促使人们对于 流固耦台问题、振动的机理及后果进行研究。由于横向振动一般要比流向振动大得多， 因此目前大量研究工作都集中于横向振动问题上；而且，考虑到问题的复杂性和相似性， 圆柱在均匀来流中的强迫振荡成为许多人研究的主要对象。 设F为反映物体振荡及流动的特征物理量， 对自激振荡而言，F=F(￡乙，D，v，‘，^，p，p呷柱)=F(Re，U，，＆h，P／O咖柱) 对强迫振荡． v，爿，正)=F(Re，彳／D，厂r／厶) F=F(U…D 体，当折合速度U，(=u。／，=rD)增加到一定范围之内时，涡脱落频率f将为物体振动频 率所俘获而与其保持一致。此时，物体振幅显著增加，相位角发生突变。这就是流致振 (1973)证实了流向振荡中也存在类似的频率“锁定”现象。 对横向振荡条件下q／D，．丘／厶)参数平面内锁定发生的区域及相应的尾涡脱落模 位于2正。附近约1．20--．2．50如的一个范围内，涡脱落与圆柱振荡几乎同步，且具体表现 为对称和反对称两种不同的涡脱落方式。Ongoren＆Rockwell”1对与均匀来流成各种不 同角度振荡的圆柱绕流问题进行了详细的实验观察，根据涡脱落的时空分布特征对涡脱 落模态进行了统一的细致分类，作者还发现，在圆柱尾流中存在着不同模态涡结构的相 互竞争。Yokoi＆Kamcmoto[zl的流向振动实验发现：涡脱落结构依赖于频率比t／t有三 种典型涡脱落方式，涡脱落频率趋向于为疋的n倍(n=o．5，1，2，3…)。有关该问题的数值计 算结果近年来也有所增加(chan窟嘲；zhaagE”)。但总体来说，对流向振荡问题的研究远没 一229 有完善。 本文集中于研究频率锁定发生的频率和幅值条件及相应的涡脱落模式。对于实验中 观察到的双子对称涡脱落模态．本文试图通过研究尾涡与圆柱的相对速度来说明其发生 机理。 1数值计算方法概述 本文的数值求解建立在同樯随时间移动的运动坐标系中，以流动速度和压力为求解 变量。时间项采用具有二阶精度的三层全隐式格式盯LFI)，对流项的离散采用QUICK格 式。 将直角坐标系下物理平面中化M0的方程通过数学关系变换到计算平面g，q，f)中 得到计算平面下的动量方程： 警+扣_砜一(V--y，h]a瞒u+专[-(…挑+(V-y。砘】等 u’ =井胁～筹，；)]+孝(害+剥 为在动网格坐标中求解绝对速度，计算网格划分为三个区域：紧贴圆柱的内层网格， 远离圆柱的外层网格，在上两个区域的中间部分为过渡层网格。内层网格点随圆柱刚体 运动，外围厢格点固定于惯性系中，而在过渡层控制体如弹性体收张变形。在三个区域 的交界处网格一阶导数光滑过渡，同时在圆柱表面和外边界处保持网格的良好正交性。 边界条