Processing math: 28%

平屋面摩擦速度的数值模拟研究

辛林桂, 周晅毅, 顾明

辛林桂, 周晅毅, 顾明. 平屋面摩擦速度的数值模拟研究[J]. 工程力学, 2024, 41(2): 151-159. DOI: 10.6052/j.issn.1000-4750.2022.03.0271
引用本文: 辛林桂, 周晅毅, 顾明. 平屋面摩擦速度的数值模拟研究[J]. 工程力学, 2024, 41(2): 151-159. DOI: 10.6052/j.issn.1000-4750.2022.03.0271
XIN Lin-gui, ZHOU Xuan-yi, GU Ming. NUMERICAL SIMULATION OF FRICTION VELOCITY OF FLAT ROOFS[J]. Engineering Mechanics, 2024, 41(2): 151-159. DOI: 10.6052/j.issn.1000-4750.2022.03.0271
Citation: XIN Lin-gui, ZHOU Xuan-yi, GU Ming. NUMERICAL SIMULATION OF FRICTION VELOCITY OF FLAT ROOFS[J]. Engineering Mechanics, 2024, 41(2): 151-159. DOI: 10.6052/j.issn.1000-4750.2022.03.0271

平屋面摩擦速度的数值模拟研究

基金项目: 国家自然科学基金项目(52078380);科技部国家重点实验室基金项目(SLDRCE19-B-14)
详细信息
    作者简介:

    辛林桂(1994−),男,河南人,博士生,主要从事结构风雪工程研究(E-mail: 1910271@tongji.edu.cn)

    顾 明(1957−),男,江苏人,教授,博士,博导,主要从事结构风工程研究(E-mail: minggu@tongji.edu.cn)

    通讯作者:

    周晅毅(1975−),男,湖南人,教授,博士,博导,主要从事结构风雪工程研究(E-mail: zhouxytj@tongji.edu.cn)

  • 中图分类号: TU973+.213

NUMERICAL SIMULATION OF FRICTION VELOCITY OF FLAT ROOFS

  • 摘要:

    建筑屋面的摩擦速度是决定其风致迁移雪荷载的重要因素。该文采用CFD (computational fluid dynamics)方法模拟了平屋面的摩擦速度,并结合其分布特征分析了屋面积雪传输率、侵蚀通量和暴露系数随跨度的变化趋势。为验证CFD方法的湍流模型对平屋面摩擦速度模拟的影响,将数值模拟结果与Irwin探头风洞试验进行了对比,结果表明Realizable k-ε湍流模型的模拟结果与风洞试验最为接近。根据4种不同跨度屋面的CFD模拟结果发现,当屋面跨度与高度之比为4时,屋面平均摩擦速度最小;并且随着跨度的增大,屋面摩擦速度最小值所在的位置和积雪沉积的区域均逐渐向迎风侧屋檐靠近。同时发现,屋面积雪的平均侵蚀通量随着跨度的增大而减小,屋面积雪传输率和雪荷载暴露系数则随着跨度的增大而增大。

    Abstract:

    The friction velocity of building roof is an important factor to determine the drifted snow load. In this paper, CFD (computational fluid dynamics) method is used to simulate the friction velocity of the flat roof. In addition, combined with its distribution characteristics, the variation trends of snow transport rate, snow erosion flux, and exposure factor with span are analyzed. To verify the effect of the turbulence model simulated by CFD on the friction velocity of the flat roof, the numerical simulation results are compared with the Irwin probe wind tunnel test. The results show that the simulation results of the Realizable k-ε turbulence model are most close to the friction velocity results obtained from the wind tunnel test. According to the CFD simulation results of four roofs with different spans, when the ratio of roof span to height is 4, the mean friction velocity of the roof is the smallest. And with the increase of roof span, the position of the minimum roof friction velocity and the area of snow deposition are both gradually closer to the windward eaves. Meanwhile, it is found that the mean snow erosion flux of the roof decreases with the increase of the roof span, while the snow transport rate and snow load exposure factor increase with the increase of the roof span.

  • 平屋面是一种典型的屋盖形式,其跨度对建筑周围的流场有着重要的影响[1-5]。风经过平屋面时会出现复杂的绕流、再附现象,在气流作用下雪颗粒将发生复杂的飘移运动,从而造成屋面和建筑周围积雪的不均匀分布[6-9]。大跨度的平屋面结构对于不均匀分布的雪荷载很敏感,雪荷载的不均匀分布形式往往是结构的最不利荷载分布之一。因此,准确预测不同跨度屋面上的积雪分布和传输特征有利于保障结构的安全。

    目前对平屋面积雪分布和传输的研究方法主要有现场观测[10-11]、风洞或水槽试验[12-13]和数值模拟[14-15]三种。现场观测记录了平屋面在自然环境中的积雪分布和传输规律,反映了最为真实积雪飘移的运动。THIIS和GJESSING[10]在挪威对平屋面进行了屋面积雪的实测研究,结果表明屋面尺寸对屋面积雪的飘移和分布有显著的影响;ZHANG等[11]采用现场观测研究了不同比例的平屋面和其他屋面模型的缩尺比对屋面积雪分布的影响,观测结果表明几何缩尺比对屋面积雪分布影响甚微。与现场测量相比,风洞或水槽试验可以提供一个稳定可控的试验平台来研究屋面的积雪飘移。O’ROURKE等[12]在水槽中利用核桃壳对平屋面积雪传输率进行了研究,结果表明,水槽模型与全比例模型的输运率吻合较好;QIANG等[13]在低温风洞中使用人工雪颗粒对积雪在平屋顶上的飘移进行了模拟,研究了降雪期间与无降雪期间积雪飘移的异同。不同于现场实测和风洞或水槽试验,对屋面积雪的数值模拟更偏向于机理研究。LIU等[14]利用拉格朗日方法建立了一个数值模型研究了平屋面的积雪飘移,探讨了击溅夹带对平屋顶积雪输运的作用及其贡献。QIANG等[15]介绍了一种积雪传输模型,用于计算平屋面的风致积雪传输量,在该模型中考虑了降雪量和阈值摩擦速度等变量的差异,模拟结果与前人的实测结果较为一致。上述的研究采用不同方法分析了平屋面积雪的分布特征,但忽略了屋面跨度可能对积雪传输造成的影响;同时,O’ROURKE等[16]通过研究发现饱和积雪传输率应按传输距离进行折减。由于目前对平屋面的积雪研究中,鲜见有学者考虑屋面跨度对屋面雪荷载的影响,因此结合屋面积雪传输率与传输距离分析跨度对平屋面雪荷载的影响是非常必要的。

    积雪野外实测的研究结果表明,传输率与摩擦速度之间存在指数关系[17-18]。并且一些研究结果表明,建筑屋面的摩擦速度是决定其风致迁移雪荷载的重要因素[19],如:康路阳等[20]采用分阶段的欧拉-欧拉方法模拟了有休止角的积雪表面的摩擦速度,并结合摩擦速度预测了平屋面的积雪分布;CAO等[21]通过风洞试验和数值模拟研究了低坡度屋面的积雪质量传输率与摩擦速度之间的关系。因此,准确获得平屋面的摩擦是研究屋面积雪传输的前提条件。

    本文首先采用CFD方法模拟了平屋面的摩擦速度,并且为了研究不同湍流模型的模拟精度,利用Irwin探头风洞试验与数值模拟结果进行了对比。接着,根据CFD模拟的结果分析了不同跨度平屋面流场和摩擦速度的分布特征。最后,基于数值模拟得到的摩擦速度,计算了平屋面积雪传输率,并讨论了不同来流风速下屋面跨度对积雪传输率和侵蚀通量的影响。

    为了获得不同跨度平屋面的摩擦速度,本文采用RANS方法对二维平屋面的周围的流场进行模拟。研究对象的跨度L分别为30 m、60 m、90 m和120 m,高度H为15 m,屋面跨度与高度之比L/H分别为2、4、6和8,如图1所示。

    图  1  研究对象
    Figure  1.  Research object

    CFD数值模拟采用ANSYS FLUENT软件进行,数值模型的几何缩尺比为1∶100,图2为计算域和网格划分示意图(未按比例画出),计算域的长度为16L,高度为20H。建筑迎风面到入口(velocity inlet)边界的距离为5L,出口(outflow)边界距离建筑背风面为10L, 如图2(a)所示。以模型L/H=4为例,高分辨率的结构化网格如图2(b)所示,最小网格尺寸约为0.03H,网格的横向和纵向的增长率均为1.05,总网格数量约为5.5万。

    图  2  计算域和网格划分示意图
    Figure  2.  Schematics of computational domain and computational grid

    计算域入口边界采用速度入口,参考《建筑结构荷载规范》对B类地貌的规定,并按照几何缩尺比(1∶100)缩放后在数值模拟中进行应用。入口平均风剖面U(y)如下所示[22]

    U(y)=U0.1(y0.1)α (1)

    式中:风剖面指数α取0.15;为了探究风速对屋面摩擦速度和积雪传输率的影响,本文取U0.1 (距地面0.1 m高度处风速)的风速大小分别为6 m/s、8 m/s和10 m/s。

    湍流强度剖面I(y)为[22]

    I(y)=0.14(y0.1)α (2)

    根据风速剖面U(y)和湍流强度剖面I(y),湍动能剖面k(y)和湍流耗散率ε(y)剖面为[23]

    k(y)=1.2[U(y)I(y)]2 (3)
    ε(y)=C1/2μk(y)U0.10.1α(y0.1)α1 (4)

    式中,Cμ 为0.09。

    下边界和建筑表面的壁面边界条件采用无滑移标准壁面函数。上边界设置为对称边界,出口边界设置为0梯度的出流边界。采用稳态RANS方法,对Standard k-ε、 RNG k-ε和Realizable k-ε三种湍流模型进行了比较分析。所有的传输方程均采用二阶迎风离散格式,压力速度耦合方法为SIMPLE算法,收敛准则为10−6

    采用三种尺寸的网格进行网格敏感性分析,分别为细网格、中网格和粗网格。细网格的最小网格尺寸为0.01H,屋面y+的范围为25~50;中网格的最小网格尺寸为0.03H,屋面y+的范围为31~70;粗网格的最小网格尺寸为0.05H,屋面y+的范围为38~95。图3为模型尺寸L/H=4时,采用Realizable k-ε模型和三种网格划分方案模拟的屋面摩擦速度对比。从图3可以看出,粗网格与细网格的差异较大,而中网格与细网格的差异较小,因此本文在下文中选择中网格的模拟结果进行分析。

    图  3  网格无关性分析
    Figure  3.  Grid-sensitivity analysis

    为了验证Standard k-ε、 RNG k-ε和Realizable k-ε三种湍流模型的模拟精度,并对比不同湍流模型对模拟结果的影响:首先,对Irwin探头进行了标定;接着,在不同来流风速下采用标定后Irwin探头测量了平屋面的摩擦速度;最后,将不同湍流模型的数值模拟和风洞试验进行了对比。风洞试验的风压测量采用美国Scanvival公司生产的ZOC33电子压力扫描模块,压力扫描阀的采样频率为312.5 Hz,采样时间为30 s。

    本文采用的Irwin探头为无方向探头[24],尺寸如图4所示,测压孔K1K2的垂直距离h为5 mm。

    图  4  Irwin探头尺寸 /mm
    Figure  4.  Dimensions of the Irwin sensor

    FERREIRA等[25]和FARIA等[26]认为壁面切应力与Irwin探头测得压差存在指数关系,即:

    τw=aΔPb (5)

    式中:ab为Irwin探头测量壁面切应力时的标定系数;ΔP = P0Ph,为Irwin探头测得的压差。由摩擦速度u = τw/ρaρa为空气密度,可得uΔP的关系为:

    u=mΔPn (6)

    式中,mn为Irwin探头测量摩擦速度时的标定系数。标定试验在同济大学TJ-2大气边界层风洞的空风场中,根据WANG等[27]的方法对Irwin探头测量摩擦速度进行了标定,如图5(a)所示。假设标定试验的近壁面风速符合对数分布:

    U(z)=uκln(zz0+1) (7)

    式中:U(z)为z/m高度处的风速,标定试验中不同风速下的风剖面U(z)分布如图5(b)所示;Von Karman常数κ为0.42;z0为标定平台的气动粗糙度高度。通过式(7)的关系对风剖面进行拟合后,可得到z0为1.173×10−5 m,五种不同风速下的u分别为0.17 m/s、0.23 m/s、0.32 m/s、0.41 m/s和0.50 m/s。于是再根据式(6)对高度h=5 mm的Irwin探头测得的ΔPu进行标定,通过重复性和方向性检验后,可得到h=5 mm的Irwin探头在测量摩擦速度时的标定系数m=0.183、n=0.338。

    图  5  Irwin探头测量摩擦速度的标定示意图
    Figure  5.  Calibration diagram of friction velocity measured by the Irwin sensor

    采用Irwin探头在同济大学TJ-2大气边界层风洞对平屋面的摩擦速度进行测量,风洞的试验段长15 m、宽3.0 m、高2.5 m,如图6所示。屋面模型的几何缩尺比为1∶100,跨度为0.6 m,高度为0.15 m,阻塞比小于5%。在进行风洞试验时仅模拟风场的二维效应,模型的两侧增加了侧挡板,并且为了减小挡板对风场的影响,挡板端部进行了削尖处理。Irwin探头在风洞试验模型上的布置如图7所示。考虑到Irwin探头相互遮挡可能会影响测量精度,故在屋面上对Irwin探头进行了交错布置,两排Irwin探头的间距为50 mm,Irwin探头距离迎风侧和背风侧屋檐的距离均为10 mm,其他的Irwin探头在屋面上均匀排布。

    图  6  风洞试验试验段示意图 /m
    Figure  6.  Schematic diagram of wind tunnel test section
    图  7  风洞试验模型测点布置 /mm
    Figure  7.  Layout of measuring points of wind tunnel test mode

    在测量平屋面的摩擦速度时,定义屋面高度H为参考高度,平屋面Irwin探头风洞试验采用了三种风速:参考风速U(H)分别为6 m/s、8 m/s和10 m/s,参考高度处的湍流强度I(H)约为16%,风速和湍流强度剖面如图8所示。参考了康路阳等[20]对风洞试验风速的拟合方法,风速剖面可由式(7)进行拟合,气动粗糙度高度z0为5.362×10−5 m,湍动能由式(2)进行计算,湍流耗散率的定义为:

    ε(z)=C1/2μk(z)U(z)z (8)

    式中,Cμ为经验常数,取0.09。

    图  8  风洞试验风速和湍流强度剖面
    Figure  8.  Wind speed and turbulence intensity profile of the wind tunnel test

    为了对比Standard k-ε、 RNG k-ε、Realizable k-ε三种湍流模型对数值模拟结果的影响,将CFD模拟获得的屋面摩擦速度u*与风洞试验结果进行了对比。以参考风速U(H)为8 m/s为例,数值模拟与风洞试验的对比如图9所示。

    图  9  数值模拟与风洞试验的摩擦速度对比(U(H)=8 m/s)
    Figure  9.  Comparison of friction velocity between numerical simulation and wind tunnel test (U(H)=8 m/s)

    在三种风速下,采用FAC2(the factor of 2)、线性相关系数R和平均百分比误差APD对数值模拟结果进行评价。FAC2、R和APD的定义如下所示[28-30]

    FAC2=1NiFi,Fi={1,12 (9)
    {{R}} = \frac{{\dfrac{1}{N}\displaystyle\sum\nolimits_i {( {{P_i} - \overline P} )} ( {{M_i} - \overline M} )}}{{{\sigma _P}{\sigma _M}}} (10)
    {{{\rm{APD}}}} = \frac{{{{\overline M}_{\rm abs}} - {{\overline P}_{\rm abs}}}}{{{{\overline M}_{\rm abs}}}} \times 100\text{%} (11)

    式中:PiMi (i=1, 2,…, N; N=48)分别为CFD模拟的模拟值和风洞试验测量值; {\sigma _P} {\sigma _M} 分别为CFD模拟和风洞试验结果的标准差;上划线代表平均值, {\overline P_{\rm abs}}{\text{ = }}\dfrac{1}{N}\displaystyle\sum\limits_{i = 1}^N {\left| {{P_i}} \right|} {\overline M_{\rm abs}}{\text{ = }}\dfrac{1}{N}\displaystyle\sum\limits_{i = 1}^N {\left| {{M_i}} \right|} ;FAC2、R理想值为1.0,APD理想值为0%。

    根据上述三种评价指标,不同湍流模型的数值模拟结果评价如表1所示。在FAC2评价中,Realizable k-ε模型的摩擦速度的评价结果与理想值最为接近,这表明Realizable k-ε模型的整体表现最好,很少发生明显过高或过低的预测。对比线性相关系数R可得,Realizable k-ε模型模拟的摩擦速度与风洞试验测量结果相比相关性最好。在平均百分比误差APD评价中,Realizable k-ε模型的摩擦速度评价结果与理想值0%最接近。综上所述,Realizable k-ε模型对摩擦速度的预测表现最好。因此,在后文的研究中选择Realizable k-ε模型的数值模拟结果进行分析。

    表  1  三种风速下不同湍流模型模拟结果对比
    Table  1.  Comparisons of simulation results of different turbulence models under three wind speeds
    评价指标摩擦速度u*
    FAC2RAPD/(%)
    理想值1.001.000.00
    湍流模型Standard k-ε0.710.58−26.36
    RNG k-ε0.740.63 −9.02
    Realizable k-ε0.780.74 −5.01
    下载: 导出CSV 
    | 显示表格

    U10=8 m/s为例,不同跨度平屋面的风速分布如图10所示。由图10可得,来流在屋盖迎风侧的前缘发生分离,迎风前缘最大风速约为来流风速U10(8 m/s)的1.5倍。来流在迎风前缘发生强烈的分离,形成旋涡后发生再附,并在屋面后方形成一个回流区。屋面跨度L=2H时(图10(a)),分离产生的旋涡几乎覆盖了整个屋面;随着屋面跨度的增大(图10(b)~图10(d)),旋涡的长度也在不断增大,但旋涡的长度与屋面跨度的比值在不断的缩小。同时,旋涡后方近屋面的风速也随着跨度的增大不断增大。

    图  10  不同跨度平屋面周围的风速分布
    Figure  10.  Wind speed distribution around flat roofs with different spans

    摩擦速度可以确定雪颗粒运动状态,在数值模拟中屋面摩擦速度 {u_*} 可由屋面剪切应力 {\tau _{\rm w}} 计算获得,即 {u_*}{\text{ = }}\sqrt {{\tau _{\rm w}}/{\rho _{\rm a}}} 。受近屋面风速的影响,不同跨度屋面的摩擦速度分布如图11所示。由于来流在屋盖前缘分离的影响,摩擦速度在屋盖前缘处较大,并且不同跨度屋面的摩擦速度在屋盖前缘处大小分布基本相等;而后经过分离产生的旋涡覆盖区域,摩擦速度先增大后减小;最后旋涡覆盖区域后方风速平稳发展,摩擦速度也逐渐增大。参考孙晓颖等[31]对雪颗粒阈值摩擦速度的取值,本文取雪颗粒的阈值摩擦速度u_{*t} 为0.15 m/s。当屋面跨度与高度之比L/H=2时,屋面摩擦在背风侧屋檐小于阈值,屋面其他位置均大于阈值,屋面整体呈现侵蚀状态。随着屋面跨度的增大,屋面摩擦速度的最小值所在的位置逐渐向迎风侧屋檐靠近,即屋面上产生积雪沉积的位置也逐渐靠近迎风侧屋檐。

    图  11  不同跨度屋面的摩擦速度分布
    Figure  11.  Friction velocity distribution of roofs with different spans

    将数值模拟获得的平屋面表面的平均摩擦速度定义为 {\overline u_*} ,来流风速分别为6 m/s、8 m/s和10 m/s时,不同跨度的平均摩擦速度如图12所示。由图12可得,屋面跨度相同时,屋面平均摩擦速度 {\overline u_*} 与来流风速{U_{10}}之比会随着来流风速的增大而减小,但总体差别很小,在5%左右。同时,在相同风速的来流作用下,屋面平均摩擦速度随着屋面跨度L与高度H之比的增大先减小后增大,其中L/H=4时屋面的平均摩擦速度在四种跨度的屋面中最小。

    图  12  不同跨度的平屋面表面的平均摩擦速度
    Figure  12.  Mean friction velocity of roofs with different spans

    本节采用CFD模拟获得的平屋面摩擦速度的分布讨论屋面跨度对屋面积雪传输的影响。屋面积雪传输率{Q_{{{\overline u}_*}}}与数值模拟获得的平均摩擦速度{\overline u_{\text{*}}}的关系可由OWEN[17]和KIND[18]的研究结果获得:

    {Q_{{{\overline u}_*}}} = a\frac{{{\rho _\rm a}}}{g}{\overline u_{*}^3}\left( {1 - \frac{{u_{*t}^2}}{{{{\overline u}_{*}^2}}}} \right) (12)

    式中: {\overline u_{\text{*}}}/(m/s)为CFD模拟获得的屋面平均摩擦速度; {u_{*t}} 为屋面雪颗粒起动的阈值摩擦速度,取0.15 m/s[31]a为无量纲系数,a{\text{ = }}0.25 + w/3{\overline u_*}[17]w为雪颗粒的沉降速度,取0.5 m/s [32-33]{\rho _{\rm a}}为空气密度;g为重力加速度。考虑到屋面跨度对积雪传输率的影响,采用O’ROURKE等[16]提出的折减公式\sqrt {L/210} 对积雪传输率进行折减,L为屋面跨度。折减后屋面积雪传输率为{Q_{\rm r}}为:

    {Q_{\rm r}} = {Q_{{{\overline u}_*}}}\sqrt {L/210} (13)

    根据屋面积雪传输率{Q_{\rm r}},屋面积雪沿跨度的平均侵蚀通量{q_{\rm r}}可以表示为{Q_{\rm r}}L比值:

    {q_{\rm r}} = {Q_{\rm r}}/L (14)

    图13为不同风速下跨度对平屋面的积雪传输率和侵蚀通量的影响。由于来流风速为屋面积雪的传输和侵蚀提供了动力,对比图13(a)图13(b)可得,屋面跨度相同时,风速越大屋面积雪的侵蚀就越剧烈,屋面积雪传输率和侵蚀通量也就越大。图13(a)显示屋面的积雪传输率随着屋面跨度的增大而增大。同时,从图13(b)可见,屋面积雪的沿跨度平均侵蚀通量随着屋面跨度的增大而减小。

    图  13  不同跨度平屋面的积雪传输率和侵蚀通量
    Figure  13.  Snow transport rate and erosion flux of flat roofs with different spans

    为了进一步分析屋面跨度对屋面雪荷载的影响,本文参考ASCE荷载规范[34]对平屋面雪荷载 {P_{\rm f}} 的定义:

    {P_{\rm f}}{\text{ = 0}}{\text{.7}}{C_{\rm e}}{C_{\rm t}}{I_{\rm s}}{P_{\rm g}} (15)

    式中: {C_{\rm e}} 为暴露系数; {C_{\rm t}} 为热力系数; {I_{\rm s}} 为重要系数; {P_{\rm g}} 为地面雪荷载,kg/m2。将 {C_{\rm t}} {I_{\rm s}} 取1.0后, {C_{\rm e}} 可以表示为:

    {C_{\rm e}} = \frac{{{P_{\rm f}}}}{{0.7{P_{\rm g}}}} = \frac{{(0.7{G_{\rm g}} - {q_{\rm r}})T}}{{(0.7{G_{\rm g}})T}}{\text{ = 1}} - \frac{{{q_{\rm r}}}}{{0.7{G_{\rm g}}}} (16)

    式中: {G_{\rm g}} /(kg/(m2·s))为地面积雪的沉积通量,即单位时间内在单位面积上的降雪质量,新降积雪密度为100 kg/m3时,这里将 {G_{\rm{g}}} 取为1.5×10−4 kg/(m2·s),对应的气象条件为暴雪天气[35] T /s为屋面风吹雪的持续时间; {q_{\rm r}} /(kg/(m2·s))为屋面积雪的侵蚀通量。

    图14为不同风速下屋面跨度对雪荷载暴露系数的影响。由图可得,屋面跨度相同时,雪荷载暴露系数随着来流风速的增大而减小;并且在同样的气候条件(风速)下,跨度越大的屋盖结构,来流风速对其屋面雪荷载的折减就越小。

    图  14  屋面跨度对雪荷载暴露系数的影响
    Figure  14.  Influence of roof span on snow load exposure factor

    本文采用CFD方法模拟了平屋面的摩擦速度,并结合其分布特征分析了屋面积雪传输率、侵蚀通量和暴露系数随跨度的变化趋势,主要结论如下:

    (1) 本文根据风洞试验对Irwin探头测量摩擦速度进行了标定,结合Irwin探头风洞试验对比不同湍流模型的CFD模拟结果发现,Realizable k-ε模型的模拟结果与风洞试验最为接近。

    (2) 对比不同跨度平屋面的CFD模拟结果可得,随着屋面跨度的增大,屋面摩擦速度的最小值所在的位置和积雪沉积的区域均逐渐向迎风侧屋檐靠近。不同跨度平屋面的平均摩擦速度随着风速的增大而增大,并且屋面跨度与高度之比L/H=4时,屋面的平均摩擦速度在本文选取的4种跨度的屋面中最小。

    (3) 屋面积雪的平均侵蚀通量随着跨度的增大而减小。并且在同样的气候条件(风速)下,跨度越大的屋盖结构,来流风速对其屋面雪荷载的折减就越小。

  • 图  1   研究对象

    Figure  1.   Research object

    图  2   计算域和网格划分示意图

    Figure  2.   Schematics of computational domain and computational grid

    图  3   网格无关性分析

    Figure  3.   Grid-sensitivity analysis

    图  4   Irwin探头尺寸 /mm

    Figure  4.   Dimensions of the Irwin sensor

    图  5   Irwin探头测量摩擦速度的标定示意图

    Figure  5.   Calibration diagram of friction velocity measured by the Irwin sensor

    图  6   风洞试验试验段示意图 /m

    Figure  6.   Schematic diagram of wind tunnel test section

    图  7   风洞试验模型测点布置 /mm

    Figure  7.   Layout of measuring points of wind tunnel test mode

    图  8   风洞试验风速和湍流强度剖面

    Figure  8.   Wind speed and turbulence intensity profile of the wind tunnel test

    图  9   数值模拟与风洞试验的摩擦速度对比(U(H)=8 m/s)

    Figure  9.   Comparison of friction velocity between numerical simulation and wind tunnel test (U(H)=8 m/s)

    图  10   不同跨度平屋面周围的风速分布

    Figure  10.   Wind speed distribution around flat roofs with different spans

    图  11   不同跨度屋面的摩擦速度分布

    Figure  11.   Friction velocity distribution of roofs with different spans

    图  12   不同跨度的平屋面表面的平均摩擦速度

    Figure  12.   Mean friction velocity of roofs with different spans

    图  13   不同跨度平屋面的积雪传输率和侵蚀通量

    Figure  13.   Snow transport rate and erosion flux of flat roofs with different spans

    图  14   屋面跨度对雪荷载暴露系数的影响

    Figure  14.   Influence of roof span on snow load exposure factor

    表  1   三种风速下不同湍流模型模拟结果对比

    Table  1   Comparisons of simulation results of different turbulence models under three wind speeds

    评价指标摩擦速度u*
    FAC2RAPD/(%)
    理想值1.001.000.00
    湍流模型Standard k-ε0.710.58−26.36
    RNG k-ε0.740.63 −9.02
    Realizable k-ε0.780.74 −5.01
    下载: 导出CSV
  • [1] 孙虎跃, 叶继红. 基于PIV技术的平屋盖表面分离泡流动结构研究[J]. 工程力学, 2016, 33(11): 121 − 131. doi: 10.6052/j.issn.1000-4750.2015.03.0245

    SUN Huyue, YE Jihong. 3D characteristics of separation bubbles around flat roofs by PIV technique [J]. Engineering Mechanics, 2016, 33(11): 121 − 131. (in Chinese) doi: 10.6052/j.issn.1000-4750.2015.03.0245

    [2] 陈波, 黄正源, 尹泽鹏, 等. 地貌类型和长宽比对平屋面建筑风荷载的影响分析[J]. 建筑结构学报, 2019, 40(7): 54 − 61. doi: 10.14006/j.jzjgxb.2017.0574

    Chen Bo, Huang Zhengyuan, Yin Zepeng, et al. Effects of terrain and length-width ratio on wind loads of flat-roof building [J]. Journal of Building Structures, 2019, 40(7): 54 − 61. (in Chinese) doi: 10.14006/j.jzjgxb.2017.0574

    [3] 周晅毅, 胡学富, 顾明, 等. 考虑摩擦速度修正的风致积雪重分布模拟方法研究[J]. 工程力学, 2017, 34(10): 19 − 25. doi: 10.6052/j.issn.1000-4750.2016.06.0432

    ZHOU Xuanyi, HU Xuefu, GU Ming, et al. Study on wind-induced snow redistribution considering modification of friction velocity [J]. Engineering Mechanics, 2017, 34(10): 19 − 25. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.06.0432

    [4]

    ZHOU X, ZU G, GU M, et al. LES and wind tunnel test on friction velocity on roof surfaces [J]. Cold Regions Science and Technology, 2018, 151: 109 − 118. doi: 10.1016/j.coldregions.2018.03.005

    [5] 崔子晗, 马文勇, 李飞强, 等. 长宽比对长方体周边风致积雪分布影响[J]. 工程力学, 2021, 38(增刊): 133 − 137. doi: 10.6052/j.issn.1000-4750.2020.05.S023

    CUI Zihan, MA Wenyong, LI Feiqiang, et al. The influence of aspect ratio on snow distribution around a cuboid model [J]. Engineering Mechanics, 2021, 38(Suppl): 133 − 137. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.05.S023

    [6] 马文勇, 崔子晗, 柴晓兵. 风速和持续时间对立方体周边风致积雪演化影响研究[J]. 工程力学, 2022, 39(1): 108 − 117. doi: 10.6052/j.issn.1000-4750.2020.12.0884

    MA Wenyong, CUI Zihan, CHAI Xiaobing. Effect of wind speed and duration on process of snow distribution around a cube [J]. Engineering Mechanics, 2022, 39(1): 108 − 117. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.12.0884

    [7]

    ZHOU X, KANG L, YUAN X, et al. Wind tunnel test of snow redistribution on flat roofs [J]. Cold Regions Science and Technology, 2016, 127: 49 − 56. doi: 10.1016/j.coldregions.2016.04.006

    [8]

    ZHOU X, KANG L, GU M, et al. Numerical simulation and wind tunnel test for redistribution of snow on a flat roof [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2016, 153: 92 − 105. doi: 10.1016/j.jweia.2016.03.008

    [9]

    O'ROURKE M, STIEFEL U. Roof snow loads for structural design [J]. Journal of Structural Engineering, 1983, 109(7): 1527 − 1537. doi: 10.1061/(asce)0733-9445(1983)109:7(1527)

    [10]

    THIIS T K, GJESSING Y. Large-scale measurements of snowdrifts around flat-roofed and single-pitch-roofed buildings [J]. Cold Regions Science & Technology, 1999, 30(1): 175 − 181.

    [11]

    ZHANG G, ZHANG Q, FAN F, et al. Field measurements of snowdrift characteristics on reduced scale building roofs based on the size effect study [J]. Structures, 2021, 32(1): 2020 − 2031.

    [12]

    O’ROURKE M, DEGAETANO A, TOKARCZYK J D. Snow drifting transport rates from water flume simulation [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2004, 92(14/15): 1245 − 1264. doi: 10.1016/j.jweia.2004.08.002

    [13]

    QIANG S, ZHOU X, KOSUGI K, et al. A study of snow drifting on a flat roof during snowfall based on simulations in a cryogenic wind tunnel [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2019, 188: 269 − 279. doi: 10.1016/j.jweia.2019.02.022

    [14]

    LIU Z, ZHOU X, GU M. A numerical model for snow drifting simulations on flat roofs using Lagrangian approach [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2022, 222: 104922. doi: 10.1016/j.jweia.2022.104922

    [15]

    QIANG S, ZHOU X, GU M, et al. A novel snow transport model for analytically investigating effects of wind exposure on flat roof snow load due to saltation [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 210(2): 104505.

    [16]

    O’ROURKE M, DEGAETANO A, TOKARCZYK J D. Analytical simulation of snow drift loading [J]. Journal of Structural Engineering, 2005, 131(4): 660 − 667. doi: 10.1061/(ASCE)0733-9445(2005)131:4(660)

    [17]

    OWEN P R. Saltation of uniform grains in air [J]. Journal of Fluid Mechanics, 1964, 20(2): 225 − 242. doi: 10.1017/S0022112064001173

    [18]

    KIND R J. A critical examination of the requirements for model simulation of wind-induced erosion/deposition phenomena such as snow drifting [J]. Atmospheric Environment , 1976, 10(3): 219 − 227. doi: 10.1016/0004-6981(76)90094-9

    [19]

    XIN L, ZHOU X, GU M. Wind tunnel test and CFD simulation of the near-roof wind speed and friction velocity on gable roofs [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2022, 225: 105009. doi: 10.1016/j.jweia.2022.105009

    [20] 康路阳, 周晅毅, 顾明. 考虑积雪休止角的屋面积雪漂移数值模拟方法[J]. 同济大学学报 (自然科学版), 2016, 44(1): 11 − 15.

    KANG Luyang, ZHOU Xuanyi, GU Ming. Numerical simulation method about snowdrift on roofs with consideration of repose angle of snow [J]. Journal of Tongji University (Natural Science), 2016, 44(1): 11 − 15. (in Chinese)

    [21]

    CAO R, YU Z, LIU Z, et al. An investigation of snow drifting on low-sloped gable roofs: Wind-tunnel tests and CFD simulations [J]. Cold Regions Science and Technology, 2021, 183(3): 103227.

    [22] GB 50009−2012, 建筑结构荷载规范[S]. 北京: 中国建筑工业出版社, 2012.

    GB 50009−2012, Load code for the design of building structures [S]. Beijing: China Building Industry Press, 2012. (in Chinese)

    [23]

    TOMINAGA Y, MOCHIDA A, YOSHIE R, et al. AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96(10/11): 1749 − 1761. doi: 10.1016/j.jweia.2008.02.058

    [24]

    IRWIN H. A simple omnidirectional sensor for wind-tunnel studies of pedestrian-level winds [J]. Journal of Wind Engineering and Industrial Aerodynamics, 1981, 7(3): 219 − 239. doi: 10.1016/0167-6105(81)90051-9

    [25]

    FERREIRA A D, THIIS T, FREIRE N A, et al. A wind tunnel and numerical study on the surface friction distribution on a flat roof with solar panels [J]. Environmental Fluid Mechanics, 2019, 19(3): 601 − 617. doi: 10.1007/s10652-018-9641-5

    [26]

    FARIA R, FERREIRA A D, LOPES A M G, et al. On the calibration of Irwin probes for flow in rectangular ducts with different aspect ratios [J]. Journal of Fluids Engineering, 2017, 139(10): 101401. doi: 10.1115/1.4036668

    [27]

    WANG W, DUN H, HE W, et al. Wind tunnel measurements of surface shear stress on an isolated dune downwind a bridge [J]. Applied Sciences, 2020, 10(11): 10114022. doi: 10.3390/app10114022

    [28]

    TOMINAGA Y, AKABAYASHI S, KITAHARA T, et al. Air flow around isolated gable-roof buildings with different roof pitches: Wind tunnel experiments and CFD simulations [J]. Building and Environment, 2015, 84: 204 − 213. doi: 10.1016/j.buildenv.2014.11.012

    [29]

    LIU Z, YU Z, CHEN X, et al. An investigation on external airflow around low-rise building with various roof types: PIV measurements and LES simulations [J]. Building and Environment, 2020, 169: 106583-1 − 106583-20. doi: 10.1016/j.buildenv.2019.106583

    [30]

    RICCI A, KALKMAN I, BLOCKEN B, et al. Impact of turbulence models and roughness height in 3D steady RANS simulations of wind flow in an urban environment [J]. Building and Environment, 2020, 171: 106617. doi: 10.1016/j.buildenv.2019.106617

    [31] 孙晓颖, 洪财滨, 武岳, 等. 建筑物周边风致雪漂移的数值模拟研究[J]. 工程力学, 2014, 31(4): 141 − 146. doi: 10.6052/j.issn.1000-4750.2012.11.0839

    SUN Xiaoying, HONG Caibin, WU Yue, et al. Numerical simulation of snow drifting around building model [J]. Engineering Mechanics, 2014, 31(4): 141 − 146. (in Chinese) doi: 10.6052/j.issn.1000-4750.2012.11.0839

    [32] 张天歌, 周晅毅, 顾明. 不同护栏透风率下桥面雪飘移的数值模拟[J]. 工程力学, 2023, 40(3): 36 − 43. doi: 10.6052/j.issn.1000-4750.2021.09.0733

    ZHANG Tiange, ZHOU Xuanyi, GU Ming. CFD simulation of snow drifting around bridge deck with different barrier porosity [J]. Engineering Mechanics, 2023, 40(3): 36 − 43. (in Chinese) doi: 10.6052/j.issn.1000-4750.2021.09.0733

    [33] 肖艳, 杨易. 拱形屋盖结构不平衡雪荷载的模拟研究[J]. 工程力学, 2018, 35(10): 152 − 161. doi: 10.6052/j.issn.1000-4750.2017.07.0551

    XIAO Yan, YANG Yi. A simulation study on uneven snow loads on an arch roof [J]. Engineering Mechanics, 2018, 35(10): 152 − 161. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.07.0551

    [34]

    ASCE/SEI 7-10, Minimum design loads for buildings and other structures [S]. Reston, VA: ASCE, 2011.

    [35] GB/T 28592−2012, 降水量等级[S]. 北京: 中国标准出版社, 2012.

    GB/T 28592−2012, Grade of precipitation [S]. Beijing: Standards Press of China, 2012. (in Chinese)

图(14)  /  表(1)
计量
  • 文章访问数:  256
  • HTML全文浏览量:  70
  • PDF下载量:  34
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-03-29
  • 修回日期:  2022-06-27
  • 录用日期:  2022-07-13
  • 网络出版日期:  2022-07-13
  • 刊出日期:  2024-01-31

目录

/

返回文章
返回