串列三圆柱的尾流特性及涡激振动研究.pdf
《串列三圆柱的尾流特性及涡激振动研究.pdf》由会员分享,可在线阅读,更多相关《串列三圆柱的尾流特性及涡激振动研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、第36卷第4期2023年8月Vol.36 No.4Aug.2023四川轻化工大学学报(自然科学版)Journal of Sichuan University of Science&Engineering(Natural Science Edition)收稿日期:2023-02-14基金项目:国家自然科学基金项目(11772273;12072298)通信作者:鲁 丽(1978-),女,副教授,博士,研究方向为流固耦合动力学,(E-mail)文章编号:20967543(2023)04001907DOI:10.11863/j.suse.2023.04.03串列三圆柱的尾流特性及涡激振动研究张 政1,
2、2,鲁 丽1,2(1.西南交通大学力学与航空航天学院,成都 611756;2.应用力学与结构安全四川省重点实验室,成都 611756)摘 要:采用重叠网格技术,研究了低雷诺数下串列三圆柱的尾流特性及中游圆柱横向振动的涡激振动响应。结果表明,对于静止三圆柱在不同节径比(2.06.0)下的流场,随着节径比的增大,尾流模态将从“2S”模式变为“2P”模式;其上游圆柱与下游圆柱的涡脱临界节径比在3.03.5之间,中游圆柱的涡脱临界节径比在2.53.0之间;对于上下游圆柱固定-中游圆柱可横向振动的串列三圆柱,其中游圆柱产生横向振动的原因是发生了结构模态失稳。虽然中游圆柱在涡脱临界节径比前后的振动响应趋势
3、相同,但是在临界节径比之后的情况将以更高的折减速度进入锁频状态,且锁频区间更宽。关键词:串列三圆柱;临界节径比;涡激振动;锁频区间中图分类号:O327 文献标志码:A引 言随着科技的进步,管束结构已经被广泛应用于各个工程领域当中,如高层建筑结构、海底隧道管线、海洋开采平台、大型机械结构中的热交换器、核电站中的蒸汽发生器等。目前针对管束结构的流致振动研究主要趋于以下两个目的:1)采取相应的措施减少共振现象的发生,以避免管束变形及疲劳造成的损害;2)充分利用管束的强烈振动,利用海洋资源进行能源的转换。因此对管束结构的绕流特性以及响应的研究十分有必要,可为实际工程提供参考。目前,对于多圆柱来说,诸多
4、学者已经针对双圆柱的流致振动问题开展了大量的研究1-4。而当圆柱个数增加为3个时,由于圆柱之间的相互影响,其 圆 柱 响 应 和 尾 流 特 性 都 将 变 得 更 加 复 杂。Harichandan等5模拟了雷诺数Re为100和200,节径比L/D为2.0和5.0的串列三圆柱流场发现,由于上游圆柱的存在,中下游圆柱受到很大且不稳定的流体力作用;在L/D=2.0、Re=100时(L为圆柱中心点距离,D为圆柱直径),三圆柱均无涡脱;在L/D=5.0时,两雷诺数下均存在明显涡脱。陈威霖等6-7系统地2023年8月四川轻化工大学学报(自然科学版)分析了节径比对三圆柱单自由度横向响应的影响,发现节径比
5、为1.2时将发生驰振现象,同时揭示了流体与结构之间相互作用的机理,并对节径比为1.2时串列三圆柱发生驰振现象的原因进行了分析。Shaaban等8着重分析了中上游静止圆柱节径比对于下游圆柱的流场及响应的影响,发现间隙流与柱体之间的相互作用使得升力系数出现高阶次频率,进而使得下游圆柱在较大的折减速度下也会大幅振动。涂佳黄等9-10对亚临界雷诺数下静止串列三圆柱和涡激振动下三圆柱不同的频率比进行了数值分析,发现涡脱频率会随着节径比的增加而显著增大,圆柱流体力系数能量谱密度(PSD)曲线的波动随频率比的增大而增强。张志猛等11主要针对上游固定条件下的三圆柱单自由度涡激振动响应和尾流特性进行分析,发现与
6、均振动三圆柱相比,此约束形式下未出现驰振现象。由于顺流向振幅比横向振幅小得多,因此学者们在研究中都以圆柱在横流方向上的振动为主,并忽略了顺流方向对横流方向振动的影响12-13。现有文献对于串列三圆柱的研究主要集中在3根圆柱均可以振动的情况,但是对于其他可能遇到的振动组合形式仍然缺乏深入的研究。因此本文将对层流条件下,三圆柱均固定以及上下游圆柱固定-中游圆柱可横向振动的两种串列三圆柱模型进行数值计算,以探究其流场特性及圆柱振动响应特性。1 计算方法与模型验证1.1 计算模型计算域尺寸为40D,70D,考虑到边界可能对圆柱产生影响,因此第一根圆柱的中心与前端和上下两端边界之间的距离均为20D,所预
7、留的尾涡长度至少为 40D,流场模型如图 1 所示。在低雷诺数(Re=150)下流场的计算采用 Laminar层流模型,左端为速度入口,均匀来流,右端为压力出口,上下采用对称边界,圆柱表面采用无滑移壁面。为了简化名称,后续将上游圆柱、中游圆柱及下游圆柱分别简称为 UC(Upstream Cylinder)、MC(Midstream Cylinder)和DC(Down-stream Cylinder)。20DL速度入口压力出口L50D20D20DUUCMCDC图1 串列三圆柱流场模型本文基于重叠网格技术,分别绘制前景网格(圆柱组件网格)和背景网格,重叠部分的数据通过插值实现数据传递。因此背景网格
8、需要对圆柱所在位置处进行加密,并且为了减小插值所带来的误差,重叠区域内各子区域网格的密度以及网格大小要尽量接近。流场计算模型如图2所示,图中前景网格与背景网格均采用结构化网格。前景网格采用直径为3D 的圆形网格,并且对圆柱表面边界层网格加密(图2(a));在背景网格上以第二根圆柱的中心为基准进行网格加密,加密区间为8D,12D(图2(b))。(a)前景网格 (b)流体域网格 图2 流场网格相比弹性网格和网格重构,重叠网格既可以保证结构在大变形下的计算精度,又可以避免因网格重构而使得计算效率降低。同时,对于本文模型,使用重叠网格法同样能避免使用网格重构法时产生网格畸变,从而避免了计算发散。1.2
9、 计算方法本文将二维圆柱简化为弹簧-质量-阻尼系统。以单圆柱为例,根据牛顿第二定律,具有两个自由20第36卷第4期张 政,等:串列三圆柱的尾流特性及涡激振动研究度的弹性支撑的圆柱控制方程可以写为:mx+cx+kx=Fd(tn)my+cy+ky=Fl(tn)(1)其中,m为圆柱体质量;k为弹簧刚度;c为结构阻尼;Fd(tn)为tn时刻圆柱在顺流向受到的阻力;Fl(tn)为tn时刻圆柱在横向受到的升力。式(1)还可以表示为:x+0 x+20 x=Fd(tn)/my+0y+20y=Fl(tn)/m(2)其中,为阻尼比,0为圆柱固有圆频率。连续性方程和粘性流体的非定常流动的 N-S方程14如下:uix
10、i=0(3)uit+ujuixj=-1pxi+v2ui(4)其中,xi、xj分别为i方向和j方向的笛卡尔坐标分量;ui、uj分别为i、j方向上的瞬时速度分量(m s-1);为不可压缩流体密度;p为流体压强(Pa);v为流体的运动粘度(m2 s-1);t为流动时间(s)。本文基于重叠网格技术,采用Fluent求解器对流场进行求解,获得作用在圆柱上的流体力,将其代入式(2)中,然后使用 C 语言编写的关于四阶Runge-Kutta法的UDF(User-Defined Function)程序对圆柱的控制方程进行求解,从而得到圆柱的振动响应。1.3 时间步长与模型验证在均匀来流作用下对节径比L/D=4
11、.0,Re=150的串列静止三圆柱模型进行时间步长的验证。表1列出了在节径比L/D=4.0、Re=150时,不同时间步长下的升、阻力系数以及其变化率,表中CMaxl为升力系数最大值,CMeand为阻力系数平均值。由表1可以看出,时间步长对本算法的影响不大。因此,将t=0.005 s作为后续的计算时间步长。表1 不同时间步长下的升、阻力系数圆柱类型UCMCDC升、阻力系数CMeandCMaxlCMeandCMaxlCMeandCMaxlt=0.001 s1.1800.3250.2601.1840.3910.365t=0.005 s1.1840.3300.2621.1870.3950.368变化率
12、/%0.331.530.771.631.020.82为保证计算准确性,对节径比L/D=4.0,Re=100的流场,计算 3 根圆柱上的升、阻力系数,并与文献15相比较,结果见表2。由表2可见,两者之间最大差值为1.44%。表2 静态验证圆柱类型UCMCDC升、阻力系数CMeandCMaxlCMeandCMaxlCMeandCMaxl本文1.2260.4430.4611.3600.2920.889文献151.2200.4450.4581.3620.2940.908差值/%0.490.450.650.150.681.442 串列三圆柱静止状态绕流规律为了探究串列三圆柱在静止状态下的绕流规律,给定R
13、e=150,L/D=2.06.0,计算并获得了静止串列三圆柱上受到的升、阻力系数和尾流模态随节径比变化的分布情况。图3(a)、图3(b)分别给出了升力系数最大值CMaxl以及阻力系数平均值CMeand在不同节径比下的分布规律。从图 3(a)可见,3根圆柱受到的升力分别在节径比L/D=2.02.5之间以及L/D=3.03.5之间发生较大变化。当L/D 3.0时,3根圆柱上CMaxl随着节径比的增大而增大,同时 DC上的CMaxl比 MC 大;而当L/D3.0 时,MC 上的CMaxl始终大于另外两根圆柱。从图 3(b)可见,UC 受到的CMeand始终大于另外两根圆柱。3 根圆柱受到的CMean
14、d同样在L/D=3.03.5之间发生较大改变,随着节径比的增大,MC上的阻力系数平均值由CMeand0;且当L/212023年8月四川轻化工大学学报(自然科学版)D3.0时,MC上的CMeand大于DC。234560.40.81.21.62.00CMaxlL/D UC MC DC(a)CMaxl23456-0.20.20.40.60.81.01.21.40CMeandL/DUCMCDC(b)CMeand图3 三圆柱静止下的CMaxl与CMeand分布图4所示为串列三圆柱静止绕流在不同节径比下的尾涡分布。(a)L/D=2.0(b)L/D=2.5(c)L/D=3.0(d)L/D=3.5(e)L/D
15、=4.0(f)L/D=5.0图4 三圆柱绕流在不同节径比下的尾涡分布从图4(a)中可以看出,当节径比L/D=2.0时,中下游圆柱的存在将抑制上游圆柱涡脱的产生。流体流经上游圆柱后形成自由剪切层,并包裹MC与DC,在3根圆柱上均未见明显的涡脱产生,而图4(a)后部产生的涡旋是由尾流的卷起效应引起的,与低雷诺数下的单圆柱情况类似。当2.5 L/D 3.0时,由图 4(b)、图 4(c)可以看出,每个周期内,由 MC 产生的涡撞击在 DC 上,与DC上产生的涡融合。融合后的涡分别在DC上下表面交替脱落,形成一排涡阶。上表面的涡为顺时针方向(蓝色),下表面的涡为逆时针方向(红色),尾流模态为“2S”模
16、式。与小节径比工况不同的是,当L/D3.5时,从图4(d)(f)可以看出,UC 后方的间隙中也开始形成涡脱,并撞击在MC上,之后被MC分为两个涡旋,分别从MC的两侧流过。若上游涡与MC上生成的涡旋向相同,则相互吸收;旋向相反的涡相互抵消。上游涡很快消散,由MC上产生的涡将顺着MC侧面流向后方。该涡流经DC两侧时,在DC两侧掠过并吸收从 DC 上产生的涡,从而形成两排稳定的涡街。此时的尾流模式与L/D=3.0时不同,MC上产生的涡不会直接撞击到DC上。因此,MC、DC的CMeand在该工况下发生变化,尾流模态为“2P”模式。对比图 4(a)(f),MC 产生涡脱的临界节径比在L/D=3.03.5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 串列 圆柱 特性 振动 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。