大跨度镂空网格结构风荷载研究
崔家春 周健 潘钧俊. 大跨度镂空网格结构风荷载研究[J]. 建筑结构,2020,50(18):16-21.
CUI Jiachun ZHOU Jian PAN Junjun. Research on wind load of large span hollow grid structure[J]. Building Structure,2020,50(18):16-21.
0 引言
镂空网格结构,是指没有覆盖屋面板的空间网格结构,在实际工程中作为太阳能板的支承结构,或者仅仅作为一种力与美统一融合的建筑表现形式。近几年镂空网格结构在我国的应用越来越多,已经建成的代表性项目有常州花博会主场馆、虹桥能源中心、南京牛首山佛顶宫大穹顶
与带封闭围护系统的屋面相比,镂空网格结构的风荷载具有以下几个特点:1)风可以穿过网格结构,每根杆件均被风场包围,相邻杆件形成相互的风荷载干扰体,所以屋面的风场更为复杂。2)风荷载分布更离散,需要获得每根杆件的风荷载,而无法像传统建筑那样获得屋面板的均布风荷载即可。3)由于需要获得杆件级别的风荷载,所以传统实体风洞试验几乎无法适应。对于大体量的镂空结构而言,建立缩尺风洞试验模型时,镂空结构的构件制作存在困难。因为风洞试验模型需要同时考虑几何缩尺以及风压测点与管路的排布,而在测压管与管路的几何尺度一定的情况下,两者的合理协调是不可回避的问题。4)在结构设计层面,需要把风荷载转化成杆件的线荷载或者节点荷载,才能在结构受力分析中予以准确的应用。因此,镂空结构由于自身的复杂几何外形特征,使得其所处的风环境以及风荷载特性非常复杂,《建筑结构荷载规范》(GB 50009—2012)
本文采用实体风洞试验和数值风洞模拟技术探索镂空网格结构平均风荷载的计算方法,并以上海G60科创云廊一期项目为背景,研究大跨度镂空网格结构的风荷载特性,为同类结构的设计提供参考。
1 数值风洞模拟可行性研究
由于镂空网格结构需要获得每根构件的风荷载,所以对于大型项目,采用缩尺实体模型进行风洞试验在模型制作、测压点布置等方面存在较大的实施困难。鉴于此,采用计算流体动力学 (computational fluid dynamics,CFD)方法进行数值风洞模拟,以获得镂空网格结构的风荷载分布。
为了研究数值风洞在镂空网格结构风荷载模拟分析中的可行性,设计如图1所示的模型。分别采用实体风洞试验、数值风洞模拟两种方法获得模型在X,Y,Z三个方向的风荷载总值,对比两种方法的差异。其中实体风洞试验(图2)采用测力天平方法获取风荷载总值
实体风洞试验和数值风洞模拟时,取B类地貌,基本风压为0.6kN/m2。模型中屋面分为4块,如图1所示。
模型的风向角定义见图3,其中0°,90°风向角方向分别与坐标系X,Y轴方向一致。
屋面的块1~3的计算结果如图4~6所示。数值风洞模拟与实体风洞试验在X,Y,Z三个方向风荷载总值Fx,Fy,Fz随风向角变化趋势一致。
数值模拟与风洞试验所得各块在X,Y,Z方向风荷载总值的最大值及其误差如表1所示。数值模拟结果总体上偏保守,其中,对于Fx,max,除块2的误差(27.0%)偏大外,其他2个分块的误差均在5%以下; 对于Fy,max,误差在33.3%~47.2%之间; 对于Fz,max,除块1的误差为20.4%外,其他分块的误差均在6.3%以下。各分块风荷载总值的最大值多出现在整体坐标Y方向对应的风向角,即90°,270°,个别出现在其附近风向角。
各块风荷载总值的最大值及误差对比 表1
分块编号 |
块1 | 块2 | 块3 | |
Fx,max |
试验结果/kN |
802.0 | 604.0 | 955.0 |
模拟结果/kN |
807.0 | 767.0 | 997.0 | |
误差 |
0.6% | 27.0% | 4.4% | |
Fy,max |
试验结果/kN |
1 083.0 | 1 033.0 | 1 470.0 |
模拟结果/kN |
1 498.0 | 1 521.0 | 1 960.0 | |
误差 |
38.3% | 47.2% | 33.3% | |
Fz,max |
试验结果/kN |
840.0 | 1 437.0 | 1 351.0 |
模拟结果/kN |
1 011.0 | 1 347.0 | 1 299.0 | |
误差 |
20.4% | -6.3% | -3.8% |
通过表1可以发现,数值模拟与风洞试验在不同风向角的风荷载总值大小总体趋势相同; Y方向的风荷载总值均大于X,Z两个方向; 大部分风向角下风荷载总值的误差可以控制在30%以内,个别风向角下风荷载总值的误差接近50%,数值模拟方法得到的结果更大。基于此,认为数值模拟结果满足建筑工程应用的精度要求。
2 截面类型对计算结果的影响分析
在工程应用中,杆件截面有工字形和箱形两种,其中工字形多为窄翼缘类型
为考察构件截面形式对风荷载计算结果的影响,构建了图7所示的分析模型,并对箱形截面、工字形截面、平板壳截面3种模型进行分析对比。计算模型的节点编号顺序与计算的风向定义如图7所示。计算输入风速为20m/s,湍动能、耗散能参数均设置为1,计算模型选取k-e模型。3种模型在顺风向的挡风投影面示意如图8所示。
图9、图10分别给出了3种模型模拟得到的0°,180°风向角作用下迎风向与竖向节点风荷载F。其中每个节点的风荷载取与其相连的所有杆件表面风荷载总和的1/2。
箱形截面、工字形截面模型的节点迎风向与竖向风荷载计算结果均较接近。平板壳截面模型的节点荷载在迎风向均偏小,这与该模型在顺风向的挡风投影面积失真有关(图8(a)),而在竖向,平板壳截面模型的节点荷载偏大。
因此,对于窄翼缘工字形截面构件镂空网格结构,可以采用箱形截面进行建模,平板壳截面过于简化,风荷载结果失真明显,在实际工程中不建议采用。
3 大跨度镂空网格结构数值风洞模拟
上海G60科创云廊一期项目平面尺寸宽126m、长734m,包括11栋80m高的钢筋混凝土超高层建筑及大跨度镂空屋面。镂空屋面采用铝合金单层网壳结构,塔楼顶部上方区域范围内的屋面覆盖有太阳能板,大约占屋面总面积的18.2%,其余区域均为镂空。镂空屋面采用树形柱支承在下部塔楼的顶部,树形柱顶部附近的部分屋面杆件为了获得更大的刚度和强度采用了钢构件。铝合金构件包括工字形和箱形两种截面形式,其中工字形构件采用侧封板装饰性封闭,所以在风荷载分析时可以全部采用箱形截面进行模拟,构件截面外包尺寸为220mm×550mm。
因为建筑体量较大,为了把计算规模控制在计算机运算能力以内,所以在CFD模型中,以镂空屋面结构水平与竖向受风面积保持不变为前提,对屋面结构的网格和构件尺寸进行了归并处理,对相邻平行向的构件进行合并,并使归并后的构件尺寸放大1倍,太阳能板的数量和尺寸则不简化。图11给出了构件归并前后网格示意图,图12给出了构件归并后的屋面几何模型。归并后典型网格尺寸由3~4m变为6.0~8.5m,构件截面外包尺寸为440mm×1 100mm。
在CFD模型中,考虑下部11栋超高层建筑、镂空屋面结构、太阳能板,按1∶1比例足尺建模,并通过对称边界条件减少数值模型的网格数量。本次数值风洞长3 500m、宽1 500m、高250m。整体几何模型如图13所示,数值风洞的CFD模型网格如图14所示,模型网格总计约2 800万个。为了研究节点等效风荷载 ,选择典型节点进行分析,图15给出了屋面宽度、长度方向中心节点编号。
计算选用k-e模型,数值风洞的入口、出口、侧面和地面边界条件如下:1)入口定义来流的平均风速和湍流强度。根据本工程周边条件,选为B类地面粗糙度类别,基本风压取0.6kN/m2。计算采用入口风剖面自保持技术,以进一步提高数值模拟的精度。2)出口采用压力出口。3)顶面和侧面采用对称边界条件。4)地面和建筑物采用无滑移壁面。
每隔22.5°为一个风向角,共计16个风向角,如图16所示。其中0°,90°风向角方向分别与整体坐标X,Y轴方向一致。
为了使设计人员方便使用风荷载数据,所以在数值风洞计算完成后,将杆件、太阳能板风荷载全部转换为对应的节点荷载,具体见图12。
图17、图18分别给出了0°,90°风向角对应的屋面风压分布云图。可以看到,屋面风压范围为-1.8~1.2kPa,风压变化总体连续、平滑。
图19给出了镂空屋面网格结构在X,Y,Z方向风荷载总值随风向角的变化曲线。可以看到,Y方向风荷载总值的极大值大于Z方向,表明该镂空屋面的风荷载除屋面法线方向的风压外,其余方向由于杆件腹板形成了迎风面,在面内的风荷载远大于传统的屋面结构,甚至是结构设计的控制荷载工况。 屋面风荷载总值见表2。
屋面风荷载总值 表2
风向角/° |
Fx/kN | Fy/kN | Fz/kN |
0 |
3 680.0 | -35.5 | 7.9 |
22.5 |
4 790.0 | 3 540.0 | 3 060.0 |
45.0 |
6 430.0 | 11 400.0 | 7 860.0 |
67.5 |
4 780.0 | 19 500.0 | 13 700.0 |
90.0 |
192.0 | 23 700.0 | 16 100.0 |
112.5 |
-4 470.0 | 19 800.0 | 14 200.0 |
135.0 |
-6 490.0 | 12 100.0 | 8 630.0 |
157.5 |
-4 460.0 | 3 420.0 | 3 290.0 |
180.0 |
-3 970.0 | -85.8 | 233.0 |
202.5 |
-4 680.0 | -3 680.0 | 2 920.0 |
225.0 |
-5 100.0 | -9 360.0 | 5 780.0 |
247.5 |
-3 540.0 | -14 000.0 | 8 000.0 |
270.0 |
-303.0 | -16 500.0 | 9 050.0 |
292.5 |
3 190.0 | -14 400.0 | 8 090.0 |
315.0 |
5 530.0 | -10 400.0 | 5 930.0 |
337.5 |
4 750.0 | -3 430.0 | 2 150.0 |
图20、图21分别给出了节点在不同剖面位置下X,Y,Z三个方向的最大等效节点风荷载Fxm,Fym,Fzm(节点所有风向角对应的等效风荷载绝对值的包络值,节点位置见图15)。可以看到,绝大部分节点在Y方向的最大等效风荷载Fym大于X,Z方向的,且Fym/Fxm,Fym/Fzm在1~5之间; X,Z方向的最大节点等效风荷载基本相当,且随节点位置不同变化幅度有限。沿结构宽度方向,Y方向最大等效节点风荷载在端部大、中间小的趋势明显。
4 结论
基于实体风洞试验和数值风洞模拟,对大跨度镂空网格结构的风荷载进行了研究,基本结论如下:
(1)镂空网格结构由于没有封闭的屋面围护系统,所以要获得其风荷载,则需要获得所有构件的表面风压,无法通过传统的屋面分块均布风压的形式获取。
(2)通过对测力天平实体风洞试验结果和数值风洞模拟结果进行对比,验证了可以采用数值风洞模拟方法获得镂空网格结构的平均风荷载。
(3)镂空屋面窄翼缘工字形截面构件与箱形截面构件在水平向与竖向风荷载差异较小,对于实际工程应用中高宽比较大的铝合金构件,为了减少数值风洞网格的复杂性,可以采用箱形截面近似模拟工字形截面。
(4)由于杆件腹板成为迎风面,镂空网格结构的水平向风荷载相对较大,甚至是结构设计的控制工况,这与传统带封闭屋面板的结构明显不同,需引起结构设计人员的重视。
本文研究了镂空网格结构的平均风荷载模拟方法、基本特征,限于技术条件,未涉及风振响应问题,有待进一步深入研究。
[2] 张雪峰,崔家春,张映洲.南京牛首山文化旅游区佛顶宫大穹顶大跨空间自由曲面屋盖风荷载特性研究[J].建筑结构,2018,48(14):14-18.
[3] 建筑结构荷载规范:GB 50009—2012[S].北京:中国建筑工业出版社,2012.
[4] 罗敏锋,姚曙.双曲结构内表皮系统在深圳机场T3 航站楼中的应用[J].施工技术,2013,42(4):68-71.
[5] 陈凯,魏庆鼎,朱凤荣,等.镂空围屏对建筑物表面风压的影响[J].北京大学学报(自然科学版),2010,46(1):1-6.
[6] 沈国辉,钱涛,杨晓强,等.设有外镂空装饰结构的扭转体型高层建筑风荷载研究[J].建筑结构学报,2013,34(6):68-74.
[7] 陈凯,魏庆鼎,朱凤荣,等.镂空材料对建筑物表面风荷载的影响[C]//第十三届全国结构风工程学术会议论文集.大连,2007:142-146.
[8] 左太辉,李庆祥,黄啟明,等.镂空双层幕墙建筑风压分布及非高斯特性研究[J].建筑结构,2016,46(16):108-112.
[9] 上海松江漕河泾新建生产及辅助用房高频天平测力风洞试验报告[R].上海:同济大学,2006.
[10] 杨联萍,邱枕戈.铝合金结构在上海地区的应用[J].建筑钢结构进展,2008,10(1):53-57.
[11] 姚念亮.铝合金格构结构的发展与应用[C]//第九届空间结构学术会议论文集.萧山,2000:28-35.
[12] 杨联萍,韦申,张其林.铝合金空间网格结构研究现状及关键问题[J].建筑结构学报,2013,34(2):1-19,60.