1、欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟欟收稿日期：基金项目：中铁二十局集团类科研课题（）；河北建筑工程学院科研课题（）作者简介：郝勇（），男，副教授，研究方向：装配式建筑结构。通信作者：王小东（），男，助教，研究方向：结构抗震。犱 狅 犻：大跨屋盖结构屋面板抗风揭研究郝勇周建彬韩天骄武航飞董华海周垠巴盼锋董文丽王小东（河北建筑工程学院，河北 张家口 ；中铁二十局集团第二工程有限公司，北京 ；天津城建大学，天津 ）摘要：以张家口奥体中心游泳馆金属屋面工程实例为背景，将风洞试验结果与 数值模拟结果进行对比，结果
2、表明风洞试验测点数据与数值模拟结果数据整体趋势吻合较好，规律性基本一致。在此基础上，设计出三种屋盖造型，利用 流体仿真软件计算出三种屋盖造型下 风向角屋盖的风压系数。为便于统计，将屋盖结构划分为个区域，计算出局部极值风压与极值风吸力，将其风压结果输入到有限元计算软件中绘制屋面板的节点位移和 应力曲线并分析屋面板的力学性能。结果表明：屋盖造型会影响风荷载的变化，而且屋盖的边缘区域为受力薄弱区域，在大风天气下更容易遭受风揭破坏，需要在实际施工过程中对屋盖边缘区域进行加固处理。关键词：抗风揭；风向角；金属屋面板；风洞试验；有限元分析中图分类号：犜 犝 文献标识码：犃 引言随着经济的飞速发展与建造技术
3、的日益革新，空间尺度更大、几何外形更新颖、屋盖表面曲率更复杂的大跨屋盖结构建筑被广泛应用在体育场馆、车站、航站楼等公共建筑中金属屋面作为大跨屋盖围护结构体系中最主要的部分，其抗风揭性能决定着整体建筑的抗风水平工程实践表明，屋面的风揭破坏已成为大跨屋盖结构建筑中最主要的破坏形式 在此背景下，研究金属屋面抗风揭性能将对我国大跨屋盖结构建筑的发展具有十分重要的意义金属屋面各构件之间多通过机械连接固定，在大风天气下会受到强大的风荷载吸力与不定常的脉动效应共同作用，引起屋面结构产生风振作用，造成构件之间连接失效，导致金属屋面发生风揭破坏目前，屋盖结构抗风揭研究多基于屋盖基本几何参数，、屋盖风向角、倾角变
4、化对风压分布影响 ，屋面板连接方式不同对抗风揭性能的影响 等，未考虑屋盖造型变化对屋面风揭的影响因此，本篇首先利用 软件对张家口奥体中心大跨屋盖结构游泳馆建立足尺模型，与风洞试验模型进行了 风向角风压的对比，之后考虑屋盖造型变化对结构的影响，计算出不同风向角下各测点的标准风压最后，将统计意义下的局部极值风压值与极值风吸力值的风压导入有限元分析软件，对比分析造型变化与受风向角度的改变对金属屋面板的力学性能影响 工程实例 工程概况游泳馆位于张家口奥林匹克体育中心东北角，屋盖外观整体为中心部分内凹的自由曲面屋盖桁架结构最大 长 度 ，最 大 宽 度 ，跨 中 厚 度 ，建 筑 最 高 点 ，水 平
5、投 影 面 积 约 整体效果如图所示，游泳馆整体结构如图所示第 卷 第期 年月河 北 建 筑 工 程 学 院 学 报犑 犗 犝 犚 犖 犃 犔犗 犉犎 犈 犅 犈 犐犐 犖 犛 犜 犐 犜 犝 犜 犈犗 犉犃 犚 犆 犎 犐 犜 犈 犆 犜 犝 犚 犈犃 犖 犇犆 犐 犞 犐 犔犈 犖 犌 犐 犖 犈 犈 犚 犐 犖 犌犞 狅 犾 犖 狅 犑 狌 狀 图奥体中心整体效果图图游泳馆整体结构图 风洞试验游泳馆风洞试验为刚性模型测压试验，几何缩尺比为：，其中屋面部分布置了 个测点模型定义来流从正北方向吹向正南方向为 风向角，并以 为间隔逆时针旋转，共测试了 风向角范围内的 个风向角工况风洞中的试验模
6、型见图，游泳馆屋面测点布置见图图试验模型图图游泳馆屋盖测点布置图张家口奥体中心游泳馆，地表粗糙度类型为类（），风洞试验中模拟得到的平均风速和顺风向紊流度（）剖面如图所示图（）中的实线和散点分别代表建筑结构荷载规范 平均风速的理论值和风洞试验实际的平均风速模拟值，可知风洞试验风速剖面与规范风速剖面整体吻合较好（）平均风速剖面（）紊流度剖面图平均风速剖面和紊流度剖面第期郝勇 等大跨屋盖结构屋面板抗风揭研究风洞试验采用无量纲风压系数来表述结构表面测点的风压：犆狆 犻，犘犻，犘狊犘狋犘狊犘犻，犘狊 犞（）式中：犆狆 犻，为屋盖结构测点犻在不同风向角下的风压系数；犘犻，为屋盖结构测点犻在风向角下总压；犘
7、狊为参考点处静压平均值；犘狋为参考高度犺处总压；为空气密度；犞为参考高度 处的平均风速根据风洞试验得到了张家口奥体中心游泳馆模型在 风向角下的表面压力分布情况风洞测压试验得到以下结论：（）张家口奥体中心游泳馆屋面的极值负压较大，且负压较大区域位于游泳馆屋面结构顶部西北侧外边缘悬挑处（）屋盖中间内凹区域的极值正、负风压较边缘区域比较小 数值模拟鉴于张家口奥林匹克体育中心建筑群巨大的体量和复杂的体型，需考虑游泳馆周边建筑干扰影响，不利于大跨屋盖结构抗风揭研究的推广因此在现有风洞试验基础上，对游泳馆采用流体仿真软件 进行 数值模拟，探究屋盖中心凹凸体型变化对金属屋面的风荷载影响计算模型采用足尺模型，
8、流域尺寸参考国内外相关文献 设置大小为 （）（）（），将游泳馆模型放置于迎风向轴三分之一处游泳馆阻塞率：（）（），流域满足要求，见图模型整体采用非结构化网格进行网格划分，游泳馆模型处进行加密处理，最终整体模型网格数量达到 万，游泳馆网格划分及局部加密见图图模型与整体流域关系图图游泳馆网格划分及局部加密图计算流域入口边界条件定义为速度入口（），采用 建筑结构荷载规范（）给定的风速沿高度呈指数率分布风速剖面，使用 进行编译表达式如下：犞犞犣犣（）（）犠犞犞 （）式中：犣、犞为某高度及对应的平均风速；犣为参考高度（取 ）；犞为参考高度处的风速，选取 年重现期基本分压 ，根据公式（）换算可得犞 ；为地
9、面粗糙度指数，类地区取 出口边界设置为完全发展出流边界（）；建筑表面和地面采用无滑移壁面（）；流域两侧与顶面采用自由滑移壁面（）风速入口湍流强度、湍动能、湍动能比耗散率参考日本规范，模拟结果验证大跨屋盖结构金属屋面抗风揭研究需要尽可能接近真实情况，为保证风压模拟的准确性，数值模拟的风载体型系数需要与风洞试验的数据进行比较，保证较高的重合度但风洞试验中周围建筑群对游泳馆存在干扰效应，故本篇仅选取干扰效应影响最小的风向角 时游泳馆屋盖风载体型系数的模拟结果为例，与风洞试验进行比较，验证数值模拟的准确性河 北 建 筑 工 程 学 院 学 报第 卷风载体型系数为实际压力与来流风压的比值，可由测点的风压
10、系数计算得到：狀 犻，犆狆 狀 犻，（犣犻犺）（）式中：狀 犻，为屋盖结构测点犻在不同风向角下的风载体型系数；犣犻为测点犻所处的高度；犺为参考点处高度；为地貌粗糙度指数图游泳馆 风向角体型系数对比图由图可以看出，风洞试验测点数据和数值模拟结果数据整体趋势吻合较好，规律基本一致，局部屋面凸起部分体型系数存在偏差分析原因：其一，风洞试验模型与数值模拟模型外部流场环境不同，风洞试验中游泳馆周围建筑群对游泳馆风压存在干扰影响；其二，风洞试验本身属于缩尺模型试验，受边界效应、支架干扰与粗糙度模拟等影响，自身存在一定的缺陷；其三，湍流模型方程采用的是与真实湍流具有较高的精度吻合度的 模型，雷诺数相似比不能
11、完全满足风洞试验与数值模拟的一致性，存在雷诺数效应影响 三种屋盖计算结果对比为了更为细致的研究大跨屋盖金属屋面抗风揭影响，在考虑张家口奥体中心游泳馆真实建筑长轴方向剖面为“元宝”型（后称模型一，高度为 ）的 风向角的基础上，进一步考虑改变屋盖中心部分造型对风压的影响结合工程实际与造型要求，增加“平屋面”（后称模型二，高度为 ）、“凸屋面”（后称模型三，高度为 ）屋盖造型，对另外两种屋面同样采用了流体仿真软件 进行 数值模拟三种屋面长轴方向剖面风迹流线示意图如图所示为了方便对不同的屋盖造型进行风荷载描述，根据模型一的造型与风荷载分布特点，将屋面分为九个部分，屋面分区如图 所示图三种屋面风迹流线示
12、意图图 屋面分区图 标准风压计算屋面属于围护结构，依据我国建筑结构荷载规范围护结构计算公式，将得到三种造型的屋面体型系数代入，得到屋面的风荷载标准值狑犽 犻，犓犜犵 狕 犻犆狆 犻，狑（）上式中，犠犽 犻，为测点犻在不同风向角下的风荷载标准值，犵 狕 犻为测点犻的阵风系数（可由建筑结构荷载规范查得），犓犜为风场转换系数，类地区犓犜为 ，犆狆 犻，为屋盖结构测点犻在不同风向角下的风压系数第期郝勇 等大跨屋盖结构屋面板抗风揭研究 计算结果分析统计三种屋盖形式 风向角下各测点的标准风压，汇总各分区的极值风压值、极值风吸值，见表（表中正值为风压，负值为风吸力）表三种模型各区域的极值风压、风吸力（）模型
13、一 模型二 模型三 从表可以看出，模型一的最大极值风压与最大极值风吸力发生在一区，模型二与模型三最大极值风压与最大极值风吸力发生在八区由三种模型正负风压数值范围得出：游泳馆屋盖中心部分凹凸造型的改变对屋盖风吸力的影响大于风压伯努力效应指出风速流动越快，其压力越小，故风吸力受屋盖造型的变化对风场绕流特性影响明显，结合图风迹流线图可知，屋盖造型的改变会导致风速产生分流效应，使得模型一、模型三较模型二拥有更大的风速综上可知，当屋盖中心区域部分高度增加时，屋盖表面风荷载并没有随着屋盖高度增加而增大，说明在大跨屋盖结构建筑中，屋盖表面风荷载主要与屋盖上方的风场绕流特性相关，即屋盖的造型较屋盖的高度更能影
14、响风荷载的变化 金属屋面板抗风揭数值分析选取金属屋面中最不利的风压与风吸力的极值区域采用静态风载分析其抗风揭性能首先，对 节 风向角三种屋盖形式各分区最强极值风压值与最强极值风吸力值进行筛选，找出各模型的最大的极值风压值与风吸力值对应的区域之后，根据屋面板的单板面积大小将该极值区域划分为若干部分，每一部分即为一块屋面板最后，对屋面板表面的标准风压基于微小单元均匀性假设理论，将屋面板上受到相同的标准风压导入有限元分析软件，对该屋面板进行力学性能模拟分析 模型建立屋面系统采用 型金属板，单块板宽 ，总长 ，板厚，截面形状如图 所示金属板抗拉强度为 ，屈服强度为 采用有限元软件进行建模，其中，屋面板
15、采用实体单元，板间接触设置为绑定，为了减小计算量简化模型，并未对支座进行建模，而是通过在两端设置位移约束来模拟边界条件，边界条件设置如图 所示图 金属屋面板截面示意图 三种屋盖形式有限元模拟结果对比风荷载施加方式根据 计算结果在金属板上施加均布压力，荷载方向始终垂直于屋面板荷载设置如图 所示河 北 建 筑 工 程 学 院 学 报第 卷图 两跨金属屋面板的边界条件设定图 两跨屋面板的荷载设置考虑到金属屋面板的破坏通常发生在较大的风压的情况下，便于计算，参照前文奥体中心游泳馆屋盖区域划分图和表所示的三种模型在 风向角下个区域的最大风压值和最大风吸力值，在有限元中对屋面板施加均布压力荷载，荷载大小为
16、风压值，荷载方向向上模拟风吸力值，向下模拟风压值由风洞试验结果可知，对于游泳馆最强的极值风吸力值发生在结构顶部外边缘处，即屋面板外边缘处由有限元软件对三种模型的模拟结果可知，最大节点位移均发生在屋面板 边缘处 节点位置处，如图 所示将三个模型各区域的 节点在最大风压值、最大风吸力值下的位移结果绘制成点线图，如图、所示；最大 应力均发生在板 的 单元位置处，如图 所示将三个模型各区域的 单元在最大风压值、最大风吸力值下的应力结果绘制成点线图，如图、所示图 板 节点位置示意图图 板 单元位置示意图 图 三种模型极值风压值下的节点位移结果图 三种模型极值风压值下的 应力结果第期郝勇 等大跨屋盖结构屋面板抗风揭研究 图 三种模型极值风吸力值下的节点位移结果图 三种模型极值风吸力值下的 应力结果 结果对比分析从上述点线图中可以看出，模型的节点位移曲线和应力结果曲线的走势大致相同，应力和节点位移成正比变化在极值风压的作用下，模型一的、区域变形与模型二、模型三相比较大；模型三的、区域变形与模型一、模型二相比较大由于屋面节点较高，模型三各个区域的屋面板均受到较大的风吸力值影响，故大部分区域的节点位移和