1、文章编号:1000-4750(2023)03-0036-08不同护栏透风率下桥面雪飘移的数值模拟张天歌,周晅毅,顾明(同济大学土木工程防灾国家重点实验室,上海200092)摘要:采用计算流体动力学(ComputationalFluidDynamics,CFD)方法,对典型桥梁断面上的雪飘移进行了数值模拟,得到桥面上风致积雪的重分布。为验证该文数值模拟方法的正确性,以平屋面风吹雪为案例,将该文数值模拟方法得到的结果与风洞试验结果进行了对比。在桥面雪飘移的数值模拟过程中,考虑了桥梁护栏的影响,对比分析了不同护栏透风率下桥面风致积雪重分布形式。研究发现:当护栏透风率大于 50%时,桥面上不会出现显著
2、的积雪沉积;当护栏透风率小于 50%时,桥面护栏附近出现了较显著的积雪沉积,且在迎风端护栏的背风侧沉积最大。为减小桥面风致积雪堆积对交通的不利影响,建议在桥梁设计时采用高透风率的护栏。关键词:计算流体动力学(CFD);雪飘移;积雪重分布;桥梁结构;护栏中图分类号:TU312+.1文献标志码:Adoi:10.6052/j.issn.1000-4750.2021.09.0733CFDSIMULATIONOFSNOWDRIFTINGAROUNDBRIDGEDECKWITHDIFFERENTBARRIERPOROSITIESZHANGTian-ge,ZHOUXuan-yi,GUMing(StateKe
3、yLaboratoryofDisasterReductioninCivilEngineering,TongjiUniversity,Shanghai200092,China)Abstract:Snowdriftingaroundatypicalbridgedeckwassimulatedusingcomputationalfluiddynamics(CFD)methodtoattainthesnowredistributiononthebridgedeck.Toverifytheaccuracyofthenumericalmodel,snowdriftingaroundaflatroofwas
4、firstlysimulated.Theresultofthenumericalsimulationwascomparedwiththatofwindtunneltests.Inthesimulationofsnowdriftingaroundthebridgedeck,thebarriersonabridgedeckwereconsideredandsnowredistributionsonthebridgedeckwithdifferentbarrierporositieswereanalyzed.Itisshownthatwhenthebarrierporosityisabove50%,
5、snowcouldnotbeaccumulatedheavilyonthebridgedeck.Whenthebarrierporosityisbelow50%,snowdriftwouldbeaccumulatednexttothebarriers,andsnowisaccumulatedsignificantlybehindthewindwardbarrier.Therefore,tomitigatetheeffectofwind-inducedsnowaccumulationonabridgedeckontraffic,itissuggestedthatthebarrierwithhig
6、hporosityshouldbeusedinthedesignofthebridge.Keywords:computationalfluiddynamics(CFD);snowdrifting;snowredistribution;bridge;barrier在我国北部的多雪地区,桥面积雪会导致车辆难以通行,致使道路受阻;由风引起的桥面不均匀分布雪荷载可能会加剧该交通问题。为减小桥面风致积雪对交通的不利影响,需要对桥面雪飘移的机理进行探究,合理地预测桥面积雪的重分布形式,进而采取相应的措施。由于桥梁附近的流场受护栏的影响较大,在探究桥面雪飘移时需要考虑护栏的影响。雪飘移问题的研究方法主要有现
7、场实测1、风洞试验2和数值模拟3。近年来,随着计算机技术的发展,基于计算流体动力学(CFD)的数值模拟方法逐渐兴起,成为工程领域中研究雪飘移问题的有效手段4,并广泛用于地面或建筑屋面雪飘移的模拟58。TOMINAGA 等7对建筑周边雪飘移的 CFD 模型进行了概述,并提出了一个模拟雪飘移的新模型。ZHOU 等对平屋面8、双坡屋面9上收稿日期:2021-09-23;修改日期:2021-11-12基金项目:国家自然科学基金项目(52078380);科技部国家重点实验室基金项目(SLDRCE19-B-14)通讯作者:周晅毅(1975),男,湖南人,教授,博士,博导,主要从事结构风雪荷载及风环境研究(
8、E-mail:).作者简介:张天歌(1995),男,山东人,博士生,主要从事结构风工程研究(E-mail:);顾明(1957),男,江苏人,教授,博士,博导,主要从事结构风工程研究(E-mail:).第40卷第3期Vol.40No.3工程力学2023 年3月Mar.2023ENGINEERINGMECHANICS36的雪飘移进行了模拟,预测了屋面上的风致积雪重分布。ZHU 等10对雪飘移的准定常模拟方法进行了改进,并将其用于高低屋面雪飘移的研究。周晅毅等11采用拉格朗日方法,对立方体周边的雪飘移进行了模拟,并考虑了颗粒运动后的质量浓度对雪层表面摩擦速度的修正。MA 等12结合现场实测和数值模拟
9、,对路堤周围风致积雪堆积进行了预测。然而,当前针对桥面雪飘移问题的研究依然较少。DUAN 等13采用数值模拟方法对桥面风致积雪重分布进行了预测,但其主要关注桥面的动力特性,没有对不同护栏下桥面雪飘移的特点进行探讨。本文基于 CFD 方法,探究不同护栏透风率下桥面雪飘移的特点。首先通过一个平屋面案例,对本文数值模拟方法的正确性进行了验证。随后采用该方法进行桥面雪飘移的数值模拟。桥面上共设置 4 道防撞栏,在模拟的过程中改变护栏的透风率(0%100%),对比不同护栏透风率下桥面附近流场形式、摩擦速度、桥面雪浓度和桥面风致积雪重分布的异同,基于模拟结果对实际工程中桥梁护栏的设计给出建议。1数值模拟方
10、法1.1数值模拟方案风致积雪运动属于两相流运动,本文采用定常 的 欧 拉-欧 拉 方 法 进 行 模 拟。空 气 相 采 用Reynolds 时间平均方法,湍流模型采用 realizablek-模型;雪相采用欧拉方法,通过附加一个雪相输运方程来模拟。模拟时假设空气相和雪相之间为单向耦合的,即雪在风的作用下发生飘移,而雪的飘移对风场不产生影响14。1.2雪相控制方程根据雪颗粒运动高度的不同,可以将雪的运动形式划分为:蠕移(小于0.01m)、跃移(0.01m0.1m)和悬移(0.1m100m),在模拟时通常将蠕移包含在跃移中15。在悬移运动时,雪颗粒会因重力作用产生向下的沉降;而在跃移运动时,由于
11、雪颗粒会在雪表面发生回弹,重力沉降的作用会被减弱。雪相的控制方程为16:t+(uj)xj=xj(tSctxj)+x3(wf)(1)式中:为雪浓度;t为湍流运动粘度;Sct为湍流施密特数;wf为沉降速度,在跃移层(wf,sal)和悬移层(wf,sus)取值不同。等式左侧分别为瞬态项和对流项,等式右侧分别为扩散项和考虑重力沉降的附加对流项。由于本文进行定常的数值模拟,因此不考虑雪随时间的瞬态变化(即瞬态项为0),只考虑雪的对流和扩散。1.3雪的侵蚀和沉积雪表面的摩擦速度是衡量雪侵蚀和沉积的关键因素,其表示雪表面所受剪切力的大小,定义为:u=0/a(2)式中:0为风作用下雪表面的剪切力;a为空气的密
12、度。雪的侵蚀和沉积通量可以采用 NAAIM 等17提出的经验公式进行计算:qero=Aero(u2u2t),u ut(3)qacc=wfu2tu2u2t,u ut(4)ut式中:Aero=0.0007 为经验系数;为阈值摩擦速度;下标 ero 表示侵蚀;acc 表示沉积。可以看出,当雪表面的摩擦速度大于阈值摩擦速度时,雪表面会发生侵蚀,反之则产生沉积。基于雪侵蚀和沉积通量的计算结果,可以得到单位时间积雪高度的变化量:h=(qero+qacc)/b(5)式中,b为雪的堆积密度。2数值模型2.1研究对象本文采用一个典型的分体箱梁断面作为研究对象18,桥面宽度(L)为 49.0m,断面高度(H)为4
13、.0m,桥面初始积雪厚度为 0.3m,桥面两侧和中部共设置 4 道防撞栏(以下简称“护栏”),护栏的高度为 1.6m(高出雪面部分),根据护栏位置将桥面划分为 5 个区域,如图 1 所示。护栏的透风率从 100%(无护栏)到 0%(不透风护栏)变化,变化间隔为 12.5%,根据护栏透风率的大小可以将其分为:无护栏、高透风率护栏和低透风率护栏三种类型,见表 1。2.2数值模拟计算域和网格划分由于桥梁轴向长度大,断面形状沿轴向基本保持不变,因此对桥梁断面建立二维模型。利用ANSYS/FLUENT 软件对桥面附近的风场进行模拟;将雪相的输运方程作为额外的方程对软件进工程力学37行二次开发,从而实现雪
14、飘移的模拟。图 2 为桥面数值模拟的计算域。网格方案采用渐变的结构化网格,网格总数为 7.0 万,在桥面附近对网格进行了加密处理。图 3 为各种护栏透风率工况下桥面附近的网格,桥面首层网格的高度为 0.05m,能够对跃移层和悬移层进行区分。经检验,在该网格方案下,计算结果具有较好的网格无关性。入口出口40 H8 L10 L加密区域(3 L10 H)图2数值模拟计算域Fig.2Computationaldomainofnumericalsimulation首层网格 度:0.05 m初始雪深:0.3 m透风率:0%12.5%25%37.5%50%62.5%75.5%87.5%100%透风率:护栏
15、度:1.6 m图3桥面附近网格Fig.3Gridaroundbridgedeck2.3数值模拟参数在数值模拟的过程中,需要对模拟参数进行合理的选取9。我国北方天气干燥,雪颗粒易在风的作用下发生飘移,且在桥面上发生飘移的雪大多是新降雪,于是本文考虑干燥、新降雪的情况,雪的堆积密度取 70kg/m31921,阈值摩擦速度取0.07m/s2223。湍流施密特数取 1.015。悬移层受重力沉降的影响较大,其沉降速度取 0.5m/s15,而跃移层受重力沉降的影响相对较小,其沉降速度取 0.15m/s。初始积雪深度取 30cm24。数值模拟的参数汇总见表 2。表2雪飘移的数值模拟参数汇总Table2Par
16、ametersinnumericalsimulationofsnowdrifting参数数值雪的堆积密度/(kg/m3)70雪表面的阈值摩擦速度/(m/s)0.07湍流施密特数1.0雪的沉降速度(悬移层)/(m/s)0.5雪的沉降速度(跃移层)/(m/s)0.15初始积雪深度/cm302.4边界条件设置来流风速在桥梁高度方向上的变化不显著,因此在速度入口采用均匀的风剖面。根据气象资料,对我国北部某地区冬季风速数据进行广义极值分布(GEV)拟合,如图 4 所示,得到具有 95%保证率的估计风速:8.5m/s。于是,在本文的数值模拟中,取入口处的风速为 8.5m/s,入口处的湍动能 k 和湍流耗散率 用下式计算25:k=32(uI)2(6)=aCk2(t)1(7)式中:I 为湍流强度,本文取 0.12;C=0.09 为常数;为湍流动力粘度;t/为湍流粘度比。速度出口边界设置为自由出流边界,上下边界设置为对称边界,雪表面采用标准壁面函数。1251258958351690835895275护栏1护栏2护栏3护栏4区域1区域2区域3区域4积雪表面区域530160H=400桥面宽度L=490012