网站首页 公文大全 个人文档 实用范文 讲话致辞 实用工具 心得体会 哲学范文 总结范文 范文大全 报告 合同 文书 信函 实用
  • 汇报体会
  • 节日庆典
  • 礼仪
  • 毕业论文
  • 评语寄语
  • 导游词
  • 口号大全
  • 其他范文
  • 百花范文网 > 实用范文 > 其他范文 > 基于K-epsilon模型的某大跨度渡槽桥梁风流数值模拟

    基于K-epsilon模型的某大跨度渡槽桥梁风流数值模拟

    时间:2023-02-27 09:20:06来源:百花范文网本文已影响

    石梓钰,张 翔,章雯雯,孙 文,邓 颖

    (安徽建筑大学 土木工程学院,安徽 合肥 230601)

    渡槽作为一种水利设施,在引江济淮这项引水工程中发挥重大作用。大跨度的渡槽桥梁多存在于高山峡谷之间,受到较大的风荷载,需要进行抗风分析[1-14]。而目前,国内较为缺乏对大跨度渡槽桥梁的抗风的模拟以及分析。

    计算流体力学[8](CFD)技术,即通过计算机的数值计算以及相关图像处理,对流体运动等物理现象进行分析模拟。CFD技术与传统的理论分析方法、实际测量方法(风洞试验)组成了流体流动分析的完整体系。李正农[1]通过定义结构风振反应的峰值因子和峰值均方差因子,对渡槽的抗风可靠度进行了分析;
    李正农[3]通过风洞试验对U形渡槽槽体的分析,得到了渡槽风压分布和风载体形系数与渡槽槽体的高宽比和流场状况有关;
    龙俊廷[12]对大跨度悬索桥桥面局部风场实测和CFD技术模拟相对比作了研究。

    为拓宽数值风洞模拟的模拟方式,现利用ANSYS Fluent流体力学计算软件进行桥梁风流模拟。通过使用Launder和Spalding等[4-6]提出及修正的三种湍流方程进行数值风洞模拟,将大跨度桥梁最常用的湍流模型Standard k-epsilon模型作为标准分析模型[10],通过使用对其修正后的两种湍流模型,来验证这两种湍流模型针对大跨度渡槽风流模拟的适用性,并探究不同分析方式下更适用的湍流模型。

    1.1 湍流模型

    湍流模型中的两方程模型是目前流体计算中常用到的模型,通过引入湍动能k和耗散率ε来形成相应的控制方程。三种湍流模型均是针对已经充分发展的湍流有效。

    1.1.1 S tandard k-eps ilon模型

    Standard k-epsilon模型是桥梁风流计算模型较为常用的两方程模型,在科学研究以及工程实践中得到了广泛的应用,对于不可压流体计算有良好的适用性。其不可压流体的两方程模型如下[4,8]:

    式中:Gk表示由平均速度梯度引起湍动能的产生项,C1ε、C2ε、Cμ为经验常数,σk、σε分别为湍动能和耗散率对应的Prandtl数,其数值如表1所列。

    表1 标准k-e ps ilon模型相关参数表

    1.1.2 RNG k-eps ilon模型

    RNG k-epsilon模型在标准k-epsilon模型的基础上考虑了平均流动的旋转,以及旋转流动的情况,对于带有弯曲壁面的流动带有更好的适应性。其不可压流体的两方程模型如下[5,8]:

    式中:Eij为时均应变率;
    其余参数如表2所列。

    表2 RNG k-e ps ilon模型相关参数表

    1.1.3 Realizable k-eps ilon模型

    Realizable k-epsilon模型在标准k-epsilon模型的基础上对湍流黏度公式进行了改变,引入了曲率和旋转的相关内容,同时对耗散率方程进行了较大改变,能够更为合理的表示能量传输。其不可压流体的两方程模型如下[6,8]:

    表3 Re a liza ble k-e ps ilon模型相关参数表

    利用以上三种湍流模型进行数值模拟与实际风洞试验结果存在一定误差[13],需要继续对湍流模型的选择与修正做进一步的研究分析。

    1.2 边界条件

    边界条件是指求解流场边界所需进行求解的变量,通过给定边界条件来驱动流场内流体的流动。为通过计算得到流场的解,必须给定合理的流场边界。

    1.2.1 流动进口边界

    定义流动进口处的流体流动速度,现设置流场进口类型为速度进口[8](velocity-inlet),流体类型为不可压缩流体。

    1.2.2 流动出口边界

    现定义计算流场出口类型为压力出口[8](pressureoutlet),且在出口边界设置相对压力,其余的流动属性由内部判断。

    1.2.3 壁面边界

    本文将桥梁内壁以及流场壁面设为光滑无滑移内壁[8](wall),对后续近壁面区的计算处理也有良好效果。

    1.3 近壁面处理方法

    由本文1.1节所介绍的三种模型均适用于已发展完全的湍流模型,这些模型均为高Re数的湍流模型。然而在近壁面区,流体运动常常会受到壁面的影响,导致湍流发展并不充分,Re数较低。尤其在黏性底层,流动几乎为层流。因此这个区域不能采用k-epsilon模型进行计算。

    壁面函数法[8](wall functions),是将壁面上的物理量与湍流核心区域待求解的物理量联系起来,直接采用半经验公式来求得近壁面区流体的流动状态。

    ANSYS Fluent软件[8]建议使用壁面函数法处理湍流近壁面问题时,距离参数处于(30≤y+≤300)。流体此时处于过渡层和对数率层,即此时速度呈对数率分布。

    式中,对于光滑壁面,κ=0.4,E=9.793。直接对过渡层和对数率层采用半经验公式求解,在不改变计算模型的前提下对湍流问题进行求解。

    现采用可缩放壁面函数[7](Scalable Wall Functions),其适合于高Re数流动,且避免了标准壁面函数在y+<11时,结果恶化[8]。该壁面函数在应对划分较粗的网格时,也能拥有良好的计算效果,以便于提高计算的准确性。

    2.1 大跨度渡槽桥梁基本模型

    某大跨度渡槽桥梁由斜拉结构和钢桁架拱组合体系组成。主跨长度为240 m,索塔高52 m,通过中央分隔带将渡槽分为双向航道,单个航道宽24 m。如图1于主跨钢桁架梁上布置50个顺桥向计算测点(A1~A50),每个测点间距4.8 m。通过流体计算来测定测点瞬时风速、风压。

    图1 顺桥向测点分布图

    横桥向计算测点布置于跨中截面处,数量为50个(B1~B50),布置方式与顺桥向布置方式类似,测点间距为1.2 m(见图2)。

    图2 横桥向测点分布图

    2.2 计算模拟工况

    2.2.1 计算模型建立

    现采用流体力学计算软件ANSYS Fluent建立大跨度桥梁计算模型。并生成相应的计算网格如图3所示。

    图3 渡槽桥梁风场计算网格图

    2.2.2 计算工况

    现将Standard k-epsilon模型作为标准分析模型,通过测定各个测点的风压、流速来比较RNG k-epsilon模 型、Realizable k-epsilon模 型 与 标 准k-epsilon模型的吻合程度,并通过数据分析得到各个模型的适用性。

    现共设置六个计算工况,各工况速度进口的速度量级(Velocity Magnitude)均设为26 m/s,其余变量如表4所列。

    表4 各计算工况设置一览表

    3.1 平均风压系数定义

    为了研究需要和便于比较,通过对计算机数值模拟的风压进行处理,将无量纲的风压系数作为研究结构表面风压分布特性的基础参数,其定义为[9-10]:

    式中:ξi表示平均风压系数;
    Pi表示观测点瞬时风压值;
    Prs表示静压力值;
    ρ表示空气密度(取1.225 kg/m3);
    Vr表示参考点风速(现选取渡槽桥面上方20 m处风速为参考风速)。

    3.2 数值模拟结果及分析

    3.2.1 数值结果

    通过ANSYS Fluent软件对本文2.2.2节设置的六个工况进行分析,得到了如图4的风压分布图。由于各湍流模型风压分布图相似,仅列出Standardk-ε模型的风压分布图。

    图4 S ta nda rd k-ε模型各工况风压分布图

    在0°风向角的风流作用下,分别对工况一、二、三的顺桥向、横桥向测点的风压值进行提取,根据公式(9)得到不同湍流模型下各测点的平均风压系数如图5、图6所示。

    图5 0°风向角顺桥向测点平均风压系数曲线图

    图6 0°风向角横桥向测点平均风压系数曲线图

    在0°风向角工况中,根据图5可知,两种湍流模型的模拟结果与Standard k-epsilon模型数值较为接近,Realizable k-epsilon模型的模拟效果比RNG k-epsilon模型的模拟效果较为良好;
    根据图6可知,两种湍流模型的数值模拟结果与Standard k-epsilon模型均存在一定偏差,但RNG k-epsilon模型模拟的总体效果较为良好。

    在90°风向角的风流作用下,分别对工况一、二、三的顺桥向、横桥向测点的风压值进行提取。根据公式(10)得到各测点的平均风压系数如图7、图8所示。

    在90°风向角工况中,根据图8可知,两种湍流模型的数值模拟结果与Standard k-epsilon模型数值较为接近。根据图7可知,两种湍流模型的数值模拟结果与Standard k-epsilon模型均存在一定偏差,但Realizable k-epsilon模型模拟的总体效果较为良好。

    图7 90°风向角顺桥向测点平均风压系数曲线图

    图8 90°风向角横桥向测点平均风压系数曲线图

    3.2.2 结果分析

    由本文3.2.1节的分析可知RNG k-epsilon模型、Realizable k-epsilon模型对于模拟大跨度渡槽桥梁风流均具有一定的适用性。由于渡槽桥梁存在较多的弯曲壁面,且受到水流作用的影响,更适合采用RNG k-epsilon模型、Realizable k-epsilon模型进行分析模拟。为更简便反应各工况下风流模拟的拟合程度,将所有工况的数值模拟结果整合分为两种情况:

    当桥梁上的测点分布方向与风流的运动方向一致时,如图5和图8所示,RNG k-epsilon模型与的数值计算效果与Standard k-epsilon模型的计算效果较为吻合,较为接近风洞试验的试验结果。两种模型均适用于大跨度渡槽桥梁的风流模拟。

    当桥梁上的测点分布方向与风流的运动方向呈90°时,如图6和图7所示,RNG k-epsilon模型与Realizable k-epsilon模型的数值计算效果与Standard k-epsilon模型的数值计算效果存在一些偏差。为显示两种模型针对各工况的适应性,继续定义了平均风压系数差值ξj,以及测点平均风压系数偏离率χj,如下所示:

    式中:ξiS表示Standard k-epsilon模型平均风压系数;
    ξSmax表示Standard k-epsilon模型平均风压系数最大值;
    ξiRNG表示RNG k-epsilon模型平均风压系数;
    ξiRea表示平均风压系数。

    由图9可知,在0°风向角工况下,为测定横桥向布置的测点,RNG k-epsilon模型的平均风压系数差值较为稳定,偏差率χj2<25%,数值计算效果吻合程度优于Realizable k-epsilon模型的数值计算效果。由图10可知,在90°风向角工况下,为测定顺桥向布置的测点,Realizable k-epsilon模型的平均风压系数差值较小,偏差率χj2<15%,数值计算效果优于RNG k-epsilon模型的数值计算效果。由此可见,在0°风向角工况下测定横桥向测点更适合选用RNG k-epsilon模型进行数值计算模拟;
    在90°风向角工况下测定顺桥向测点更适合选用Realizable k-epsilon模型进行数值计算模拟。

    图9 0°风向角横桥向测点平均风压系数差值曲线图

    图10 90°风向角顺桥测点平均风压系数差值曲线图

    (1)采用ANSYS Fluent软件进行大跨度渡槽风流模拟具有良好的适用性,但相较于风洞试验存在测量精度的限制,可以通过对湍流模型的选取、修正,以及对壁面函数的选取进行进一步的精确处理。

    (2)Standard k-epsilon模型作为桥梁最常用的数值模拟模型,为RNG k-epsilon模型、Realizable k-epsilon模型提供了良好的参考标准,且后两种模型在Standard k-epsilon模型基础上进行了一部分改进,对于模拟大跨度渡槽桥梁风流有着更良好的适用性。且针对不同风向角以及不同的测量工况下选取不同的分析模型可以提高数值模拟的精确度。

    猜你喜欢 渡槽风压湍流 浅谈水利工程渡槽基础承台施工建材发展导向(2022年24期)2022-12-22天山煤电公司106 煤矿自然风压的规律研究与应用煤炭与化工(2022年5期)2022-06-17论工况环境温度对风压传感器精度的影响科学与财富(2021年33期)2021-05-10“湍流结构研究”专栏简介空气动力学学报(2020年1期)2020-11-29傍山道路泥石流渡槽设计中华建设(2020年5期)2020-07-24——以东风渠天宫桥渡槽和黎家沟渡槽为例">渡槽工程结构安全监测方案应用
    ——以东风渠天宫桥渡槽和黎家沟渡槽为例水利与建筑工程学报(2020年2期)2020-05-31高架渡槽不同支撑结构型式动力特性研究水利规划与设计(2020年1期)2020-05-25均匀来流下方柱表面风压非高斯特性的流场机理同济大学学报(自然科学版)(2020年2期)2020-04-07重气瞬时泄漏扩散的湍流模型验证中国特种设备安全(2018年10期)2018-12-18湍流十章国外科技新书评介(2014年5期)2014-12-17

    相关热词搜索:渡槽 数值 桥梁

    • 范文大全
    • 说说大全
    • 学习资料
    • 语录
    • 生肖
    • 解梦
    • 十二星座