1、第 卷 第期 年月合 肥 工 业 大 学 学 报（自 然 科 学 版）（）收稿日期：；修回日期：基金项目：国家自然科学基金资助项目（）；安徽省自然科学基金杰出青年基金资助项目（）和中央高校基本科研业务费专项资金资助项目（）作者简介：郅伦海（），男，山东菏泽人，博士，合肥工业大学教授，博士生导师：大跨度体育馆风效应的大涡模拟及风洞试验的对比研究郅伦海，程磊，孟磊，刘毛方（合肥工业大学 土木与水利工程学院，安徽 合肥 ；安徽省金田建筑设计咨询有限责任公司，安徽 合肥 ）摘要：文章利用大涡模拟（，）湍流模型对合肥某大跨度体育馆的风效应进行研究，得到体育馆表面平均风压和脉动风压的分布规律，并与刚性模型
2、测压风洞试验数据进行对比，验证数值模拟的可靠性。分析结果表明，副馆及周边建筑群对主馆平均风压和脉动风压分布有较大影响，数值模拟得到的风压系数与风洞试验结果吻合较好。研究结果可为大跨度体育场馆的抗风设计提供技术参考。关键词：大涡模拟（）；风洞试验；风压系数；大跨度体育馆中图分类号：文献标志码：文章编号：（），（，；，）：（）：（）；因为大跨度屋盖结构的自重轻、阻尼小、柔度大，自振周期与风速的卓越周期接近，通常位于大气边界层中风速变化大、湍流度高的近地区域，主要受到气流的分离和再附流作用，所以大跨度屋盖结构对风荷载极其敏感。大跨度屋盖结构的风荷载特性不仅与结构造型有关，还与所处的建筑环境有关。周边
3、群体建筑分布愈加密集，大跨度屋盖结构受到的遮挡效应和狭缝效应会愈加显著。现场实测和风洞试验数据显示，在大跨度屋盖迎风前端，气流分离现象严重，脉动风压极值会达到平均风压的倍以上，其影响十分严重。因此有必要开展群体布局下大跨度屋盖结构平均风压及脉动风压的研究。本文以合肥某大跨度体育馆为工程背景，利用计算流体动力学软件 ，采用大涡模拟（，）湍流模型，基于随机湍流生成法得到入口脉动风速，对体育馆主副两馆进行风荷载模拟，计算并分析主馆表面的平均风压系数及脉动风压系数分布特性，并将计算结果与风洞动态测压试验进行对比。结果表明数值模拟较好地反映了群体布局下主馆表面风压分布特性，得到的风压系数与风洞试验数据有
4、较好的吻合。本文方法和研究结果对于确定群体布局下，复杂体型的大跨屋盖结构屋面的风荷载具有一定的工程意义和科学价值。工程概况该大跨体育馆由主馆和副馆组成，如图所示。副馆南北方向跨度为，东西方向跨度为 ，屋盖侧立面由四边形与三角形构成，屋面由倾斜的四边形屋面板拼接而成，相邻屋面板空隙处由三角形屋面板封闭，屋面最高点标高为，最低点标高为；主馆南北方向跨度为 ，东西方向跨度为 ，屋盖构造与副馆相同，屋面标高为。主馆与副馆之间水平最短距离为。体育馆东面无遮挡，北面、西面、南面存在较多高低层建筑，其中最高建筑高达。合肥市 重现期基本风速为，按 类地貌考虑。图建筑效果图大涡模拟数值风洞计算控制方程大涡模拟采
5、用空间过滤方法，分离湍流流场中的大尺度涡和小尺度涡；对大尺度涡，采用直接模拟直接计算；对小尺度涡，采用亚格子模型来模拟其对 流 场 的 影 响。过 滤 后 的 连 续 性 方 程 和 方程式如下：（?）（）（?）（?）?（）（）其中：为流体的密度，当流体是不可压缩流体时，；、为速度分量，经过单元平均后得到上方带有短线的变量；为压力；为运动黏度系数；?为亚格子应力，模拟小尺度涡对湍流流场的作用，为求解控制方程，需采用亚格子 模 型 来 模 拟，本 文 选 用 广 泛 应 用 的 模型。计算模型及网格划分数值风洞计算几何建模时按实际尺寸建模，考虑周围 内建筑群，建立体育馆群体模型，几何模型如图所示
6、。文献 表明，当计算域阻塞率小于时，可以忽略计算域边界对建筑表面风压分布的影响，阻塞率计算公式为：（）（）其中：为建筑物最大的迎风面积；为计算域界面的面积。图体育馆群体模型基于以上要求，并考虑数值模拟准确性和模拟计算成本，经多次试算后，确定计算域尺寸为 ，经过计算，最大阻塞率为，满足要求。为提高计算效率，将整个计算域分成个部分：一部分为外围区域，其形状规整；另一部分为核心区域，包含体育馆模型及周边建筑，如图所示（单位为）。模拟不同风向角来风工况时，将外围区域旋转相应角度，即可建立该工况下数值模型。针对不同风向角工况，分别进行网格划分。由于该大跨度体育馆造型复杂，特别是屋面凸起部分存在尖角，难以
7、在计算域中生成高质量六面体结构化网格。为正确模拟大气边界层的风场特性，获得气流撞击复杂造型屋盖后的流动现象；划分高质量网格，提高计算效率；故在 中，采用多面体网格对整个计算域进行网格划分。在网格划分前，先在建模软件 第期郅伦海，等：大跨度体育馆风效应的大涡模拟及风洞试验的对比研究 中将核心区域与外围区域的交界面、交界线共享拓扑，使个计算域的交界面、交界线在网格划分时共节点，提高数据传递效率和计算收敛性，交界面在 中自动识别为 。外围区域形状规整，采用较大尺度网格进行划分，面网格尺寸为；核心区域形状复杂，采用较小尺度网格进行划分，裙房面网格为，各风向下，体育馆表面网格最小尺寸为，最大尺寸为。总体
8、网格数量为 左右。在体育馆表面生成层贴体棱柱体网格，第层壁面网格高度取 ，网格尺寸增长率为，最大网格歪斜度为。文献 证明了高雷诺数瑞流下 壁面函数对近壁面处理的正确性。本文体育馆近壁面值为，在对数率区域内，使用壁面函数对近壁面进行处理，使近壁面满足 湍流模型的要求，保证模拟结果的准确性。图计算域示意图边界条件入口边界采用速度入口（）。根据文献，大气边界层内风速随高度变化规律采用指数律风剖面来描述，即（）（）其中：为任意高度处的平均风速；为标准高度处的平均风速；为建筑物计算位置离地面的高度；为标准高度，取；为地面粗糙度指数，在类地貌条件下，取值为.。通过湍动能和湍流耗散率来确定入口处湍流边界条件
9、，由等效边界层假设给出表达式，即（）.（）（）（）（）（）其中：为模型常数，取值为.；为卡门常数，取值为.；、分别为湍流强度和湍流积分尺度。参考文献，、的计算公式分别为：.，；.（）.，；.，（）（）.（）其中：；。入口湍流通过 方法生成，由于计算域出口距离建筑物足够远，可认为出口处尾流已充分发展，出口边界采用自由出流边界（），流场变量散度为；计算域的顶面和侧面采用对称边界（），在 中对称边界即无壁面剪应力的自由滑移壁面；地面和建筑物表面采用壁面条件（），为无滑移光滑壁面。求解参数设置在数值求解中，认为空气是不可压缩流体。采用 亚格子模型对 方程进行求解，计算离散方程组采用分离式解法中的 算法
10、来求解。为保证数值离散精度，动量方程中的对流项采用数值耗散低的二阶中心差分格式，时间项采用二阶隐式格式。为加快数值模拟计算收敛，在 进 行 非 定 常 大 涡 模 拟 计 算 前，先 采 用 模型进行定常计算，并把计算结果瞬态化处理，作为 的初始流场。计算的时间步长取，相当于 的采样频率，能正确反映实际脉动风压特性，共进行 步的 非定常计算。数值模拟结果分析大跨度体育馆为封闭结构，结构表面某点的风压系数计算公式为：（）（）（）结构表面点处的平均风压系数的计算公式为：?（）（）同时由（）式、（）式可得，结构表面点的脉动风压系数均方根为：?（）?（）其中：（）为试验模型上第个测压孔所在位合肥工业大
11、学学报（自然科学版）第 卷置的风压系数；（）为该位置上测得的表面风压值；、分别为参考点处测得的平均总压和平均静压，参考点高度取 处；为某一时刻；为时域点数。模拟时将正北方向定义为 风向，逆时针旋转，每隔 为一个试验工况，共 个试验工况，如图所示。限于篇幅，本文仅给出 风向角工况下数值模拟结果。图风向角示意图平均风压系数主馆表面平均风压系数分布如图所示。经仔细研究 风向角工况下主馆表面数值模拟结果，可以得到如下结论：（）主馆北立面（迎风面）由屋盖部分和主体部分组成。屋盖部分以正压为主，主要受压力影响。屋盖形状对风压影响较大，当来流遇到棱角部位，会发生流动分离，形成离散的漩涡，使局部的压力增大。由
12、于屋盖形状不对称，屋盖下部东侧倾斜率大于西侧，流动分离现象更严重，导致屋盖东侧平均压力系数大于西侧。屋盖部分表面最大平均风压系数为，位置出现在屋盖东侧上部。由于副馆的存在，对主馆起到遮挡作用，且气流经过副馆会发生流动分离，在主馆、副馆之间形成漩涡致使北立面主体部分平均风压系数呈上下对称分布，中部以负压为主，平均风压系数达到 。（）主馆整体造型大致对称，东西立面形状相同，平均风压系数分布规律相同，均为负值。气流在东西立面前端发生分离，平均风压系数沿立面方向递 减，最 大 负 压 系 数 发 生 在 立 面 前 端，值 为 。平均风压系数等值线基本垂直于气流流动方向，立面前端等值线密集，平均风压系
13、数变化梯度较大；沿立面方向，平均风压等值线逐渐稀疏，变化梯度逐渐减小。对比东西立面平均风压系数分布云图可以发现，由于主馆西侧存在建筑群，气流流经时会发生绕流，对主馆西面风压产生影响，致使西面平均风压系数绝对值略大于主馆东面。图主馆表面平均风压系数分布云图第期郅伦海，等：大跨度体育馆风效应的大涡模拟及风洞试验的对比研究（）主馆南立面（背风面）的风压以负压为主。虽然主馆造型大致对称，但由于西面建筑群的存在及气流的绕流作用，南立面风压系数并不是对称分布。平均风压系数绝对值由西至东逐渐减小，由负变正。最大负平均风压系数出现在凹进部分西侧，值为；最大正平均风压系数出现在凹进部分东侧，值为。（）主馆顶面的
14、风压以负压为主，整体上平均风压系数绝对值沿来流方向递减，等值线逐渐变稀疏，变化梯度逐渐减小。最大负平均风压系数出现在屋盖顶面前缘，值为。由于屋面上存在凸起部分，当气流流经时，会发生明显的流动分离及回流现象。流动分离会使凸起部分后方形成较大的负压区，最大负压系数为；回流会使屋面凸起部分前方压力系数由负变正，最大正压系数为。（）由规范类同体型的体型系数反推得到主馆表面平均风压系数，迎风面为 ，侧风面为 ，背风面为 ，屋面为 ；模拟结果整体上略小于规范值，在屋盖前端及凸起部位气流分离严重，数值模拟结果大于规范值。脉动风压系数脉动风压系数反映体育馆表面湍流特性，其大小及分布特征受到来流湍流特性和漩涡脱
15、落特性等因素影响。主馆表面的脉动风压系数分布云图如图所示。从图可以看出：（）主馆北、东、西立面的脉动风压系数分布规律与平均风压系数分布规律相似，但脉动风压系数等值线分布更加复杂，这是受到来流中的湍流成分影响造成的。图主馆表面的脉动风压系数分布云图（）来流在顶面前端迎风边缘发生分离、附着现象，因此整体上顶面前端脉动风压系数大于顶面后端，前端的脉动风压系数的变化梯度大于后端；在屋面凸起部分，气流流动分离现象严重，因此顶面凸起部位脉动风压系数急剧增大。（）因为主馆造型整体上对称，所以南立面脉动风压系数分布大致对称；因为体育馆西侧、西南侧建筑群的干扰作用，所以南立面西侧脉动风压系数略大于东侧。数值模拟
16、与风洞试验对比风洞试验该大跨度体育馆项目的风洞试验在湖南大学进行，采用 直流式矩形截面边界层风洞进行高速试验。体育馆模型是由 板制成的刚体模型，模型缩尺比为 ，高度为 ，如图所示。为了测得体育馆主馆表面的风压分布，在模型个立面及顶面布置了 个测点，如图所示。试验工况与数值模拟一致，以北面来风定义为 风向，逆时针旋转，每隔 为一个工况，共 个工况，风向角示意图见图。合肥工业大学学报（自然科学版）第 卷图风洞试验模型图各测点布置该体育馆周边建筑密集，为 类地貌，即。为模拟类地貌的大气边界层流场，在建筑模型前方设置尖塔阵和粗糙元。在测压模型风洞试验中，为保证试验的准确性，应准确模拟平均风速、风的湍流强度及其积分尺度等，风洞试验的风速剖面及湍流度剖面如图所示。图风洞试验的风速剖面及湍流度剖面第期郅伦海，等：大跨度体育馆风效应的大涡模拟及风洞试验的对比研究在大气边界层中，平均风速剖面符合指数分布律，其中参考高度 的风速 ，湍流强度自下而上逐渐衰减，在近地面达到最大值。平均风压系数对比在体育馆主馆迎风面、背风面各选取个典型测点，个侧风面各选取个典型测点，顶面选取个典型测点，共计 个测点进行分析，并
