简谐振动条件下平板表面附近速度场特性研究论文.pdfVIP

第十四届全国结构风工程学术会议论文集 简谐振动条件下平板表面附近速度场特性研究 黄林 (宁波工程学院，浙江宁波315016) 摘要基于二维不可压缩N．s方程结合RNGk-5湍流模型，对扭转和竖向简谐振动条件下的平板绕流进行了 数值模拟．应用位移厚度概念对简谐振动过程中平板表面附近的速度场进行的分析计算表明，简谐振动条件下平板 存在钝化现象；同时结合兰金涡核模型，通过旋涡对位移厚度的影响实现了速度场中对旋涡的定位及影响范围的确 定，发现平板扭转简谐振动和高频竖向简谐振动过程中上下表面附近靠迎风侧存在交替生成的旋涡．该旋涡在从迎 风侧向背风侧移动的过程中，存在从生成到变大到消失的全过程，未移动到尾部脱落． 关键词平板、位移厚度、简谐振动，运动钝化，兰金涡核 1． 引 言 旧Tacoma大桥的风毁引起了桥梁界对桥梁抗风的高度关注。关于Tacoma大桥的风毁，存在3 扭转颤振组合作用是旧Tacoma大桥风毁的原因。60多年来，世界各地的研究者们围绕桥梁的风振 响应机理展开了大量的研究。自英国Severn桥建成以来，流线形扁平箱梁在现代大跨度悬索桥中得 到了广泛应用。但目前世界上跨径最大的明石海峡大桥则采用了类似金门大桥的桁架加劲梁。这说 明，实际桥梁结构气动性能的好坏取决于气流和桥梁间的流固耦合作用。 河，采用离散涡方法计算平板的绕流场并提取气动导数。其后，不少研究者都对平板的流场及气动 导数进行过研究5~81。随着PIV(Particle Image 梁振动过程中的流场特性已成必然趋势。应用PIV技术已成功获得了分离双箱桥梁断面尾迹流场外、 导流板对桥梁底部旋涡的影响10)及H形断面振动过程中旋涡的生成变化规律¨)等流场细部信息，为 深入研究大跨度悬索桥的颤振失效奠定了基础，同时也为开展桥梁抗风的数值模拟提供了可靠的参 考依据。 本文采用改进的网格分区方法8)计算平板作简谐振动的流场，借用稳态流场中的边界层的位移 厚度概念12)，对平板上下表面的位移厚度进行了研究。进一步结合兰金涡核模型13)，对平板简谐振 动过程上下表面附近的旋涡进行了研究。 2．平板简谐振动数值模拟 计算模型如图1所示，平板宽B=0．7m，宽厚比200。计算区域如图2所示，入口距平板模型中 心6m，出口距平板模型中心14m，上下距平板模型中心4m。采用改进的网格分区方法8)将整个流 场网格划分为四边形刚性边界层网格、三角形动网格和四边形静止网格，平板表面贴体网格高度采 来流攻角为零度，以板宽计算雷诺数Re=7．2×105，来流湍流强度根据风洞试验采用O．5％；上下边 界和出口为压力边界条件，参考压力为零，湍流强度采用0．5％。计算时，首先固定平板模型作定常 计算，待流场计算比较稳定时将定常计算结果的流场作为平板强迫振动的初始流场，对平板分别施 弟十日届奎国结构AIg学术叁试论文集 2Afi(向上为正)作非定常计算， 加扭转及竖向简谐振动a0)=吼sin2研(逆时针为正)、h0】=丸sin 4 s。 湍流模型进行数值模拟计算。计算迭代时间步长山=3,0x10 目I平扳计算模型(单位mm 目2平桓计算目格 3．速度场特性研究 一 3l静止平板 在高雷诺数条件下，按边界层理论”’位移厚度很小，若要实现对边界层的分辨，需要将边界层 范国内的网格划分得非常密，这将导致简谐振动的选代时问过长。故为检验流场计算结果的台理性， 将来流风速改为匕=0Irrds(咀板宽计算雷诺数Re4792)．采用遗代时间步长山=0005s，其余参 数不变。由平均速度场计算得到的平板上下表面的位移厚度与Blasius的理论解”’的对比如图3所示． x为平板沿流向的坐标(原点为平板中心)。图中实线为计算结果，点线为理论解，上／r部曲线分别 对应平板『．门：表面(下部曲线的计算值反