1. 工程概况
钦州北部湾中心项目位于广西壮族自治区钦州市主城区,项目规划设计为超高层主楼、高层附楼、裙楼的城市综合体。超高层主塔楼平面设计为43.6 m × 43.6 m方形布局,高220 m,采用框架核心筒结构。高层附塔楼平面设计为48.4 m × 21.2 m长方形布局,高83 m。图1为项目的效果图。由于钦州属于强/台风地区,为了保证该项目在使用阶段的抗风安全性,对主塔楼进行了刚性模型高频压力积分试验。基于风洞试验结果,本文对该超高层建筑物表面极值风压分布特性、塔楼主体风荷载分布以及顶部楼层舒适度进行了详细讨论,并与商业结构分析软件提供的结果及荷载规范建议值进行了对比分析。
2. 风洞试验
2.1. 试验设备及测量系统
试验在湖南科技大学风工程试验研究中心的大气边界层风洞试验室中进行。该风洞为直流单试验段风洞,截面为4.0 m × 3.0 m (宽 × 高)的矩形,试验段的风速在1.0~30.0 m/s内连续可调。试验过程中,风速测量采用澳大利亚TFI公司的眼镜蛇三维脉动风速测量仪,压力测量采用美国Scanivalve公司的Hyscan-1000电子扫描阀测压系统。
2.2. 风场模拟
试验开始之前,采用二元尖塔和粗糙元来模拟大气边界层平均风速和湍流强度分布。根据该高层建筑物周围的地貌情况,并结合我国《建筑结构荷载规范》(GB 50009-2012) [1] 中对地面粗糙度的分类规定,在本次试验中模拟了相应的B类地面粗糙度风场,平均风速剖面指数α = 0.15。图2为流场调试后得到的B类地面粗糙度的平均风速剖面和湍流强度分布剖面。
2.3. 试验模型及测点布置
测压试验模型采用ABS板制作,几何外形与实际建筑保持相似。主体模型具有足够的强度和刚度,在试验中不会发生明显的变形和振动现象,可以保证压力测量的精度。风洞试验的模型缩尺比均为1:300,

Figure 2. Mean wind speed and turbulence intensity profiles
图2. 平均风速剖面及紊流度剖面
满足堵塞度 < 5%的要求 [2] 。0˚风向角来流时,测压试验模型整体布置如图3所示。
为获取建筑物的表面风压,在主模型表面上共布置了471个测点,沿高度共设置15层(图4)。建筑物顶部为外露的玻璃挡板,为获取外、内表面压力差,在该处布置了28对双测点,分别同步测试外、内表面的风压;2~14测点层每层布置28个单测点(图5)。每个测点埋设外径为1.2 mm、内径为1.0 mm的

Figure 3. Models in wind tunnel test
图3. 风洞试验模型布置图
不锈钢管,测压管垂直于建筑物外表面,并保证与模型表面平齐。测压孔与电子扫描阀之间采用PVC管连接,PVC管的长度控制在1 m以下,以保证测量压力的可靠性 [3] 。
2.4. 试验工况
为考虑周边建筑物对主塔楼表面风效应的干扰作用,刚性模型测压试验分两种情况进行:第一种情况为单体(不考虑周边建筑物的影响);第二种情况为所有建筑(考虑周边建筑物的影响)。不同情况下,均以每个风向角为一个工况,风向角从0˚~360˚变化,顺时针变化,变化间隔为10˚,同时增加了45˚、135˚、225˚、315˚的试验结果,试验一共有80种工况。试验时取参考点高度为70 cm,对应大气边界层风场原型的高度为210.0 m。试验控制风速为12 m/s。试验过程中,不同角度下各测点风压时程采样频率为332 Hz,采样时间为30.12 s,即每个测点在每个风向角下记录了10,000个风压时域信号数据。

Figure 5. Definitions of pressure taps and wind loads
图5. 标准层测点布置图
2.5. 数据处理
模型试验中符号约定以压力向内(压)为正,向外(吸)为负。表面各点的风压系数由下列公式给出:
(1)
式中:
是试验模型上第i个测压孔所在位置的风压系数时程,
是该位置上测得的表面风压时程,
和
分别为参考点处测得的平均总压和平均静压。外露的玻璃挡板处的内外风压处理方法参考文献 [4] 。
对风压系数时程
进行数学统计即可得到平均风压系数
和脉动风压系数
。极大风压系数
和极小风压系数
可以分别用平均风压系数和脉动风压系数的代数和来表示:
(2)
(3)
其中:
是峰值因子。由于脉动风常近似作为高斯过程考虑 [5] ,所以本文假设脉动风压也是高斯分布,取峰值因子g = 3.5,符合实际工程需要。
通过对不同测点层的风压结果进行积分,便可得到各测定层的风荷载时程。在一阶线性振型的假定下 [6] ,根据随机振动和结构动力学理论,可得到结构的均方根位移响应:
(4)
相应的脉动风引起的各层等效静力风荷载:
(5)
式中
是峰值因子,本文取2.5,具有99.38%保证率 [7] 。建筑结构各楼层最大等效静力风荷载等于平均风荷载与脉动风引起的等效静力风荷载之和,最小等效静力风荷载等于平均风荷载与脉动风引起的等效静力风荷载之差,即:
(6)
(7)
结构楼顶高度z的加速度响应
的功率谱与相应的位移相应功率谱有如下关系:
(8)
因此加速度均方根响应可如下式计算:
(9)
结构楼顶加速度峰值为:
(10)
3. 结果分析
3.1. 极值风压分布特性
极值风压是幕墙设计关注的重点。为了讨论周边建筑物对主体建筑表面极值风压分布的影响,图6和图7分别绘出了全风向角下主体建筑在无、有周边建筑物干扰情况下的极大风压系数和极小风压系数分布的等值线图。
3.1.1. 全风向极大风压系数
无周边建筑物干扰时,从图6(a)可以看出,各立面的极大风压系数沿高度的变化规律相近,在建筑物的中上部及立面拐角区域较大,表现为较大的压力,该区域的极大风压系数均匀分布在1.4附近,与荷载规范建议值1.45~1.47 (150 m~200 m高度范围内)接近。中下部的压力相对较小,大多在1.1~1.2之间。在南立面的下部中间区域风压变化剧烈,数值也减小到0.9,这主要是由于受到高度较低(83 m)的附楼的遮挡作用。在顶部外露的玻璃挡板处,气流表现为强烈的三维绕流效应,极大风压系数最大达到1.93,约为规范建议值1.44的1.34倍,应该引起幕墙设计人员的重视。
有周边建筑物干扰时,从图6(b)可以看出主体建筑各立面的极大风压系数沿高度的变化规律与无周边建筑干扰时基本一致。局部区域的风压等值线变化加剧,这可能与扰动作用下气流的脉动增强有关。中上部的极大风压系数维持在1.5左右,比无周边建筑干扰时略有增加。
3.1.2. 全风向极小风压系数
无周边建筑物干扰时,从图7(a)可以看出,主体建筑各立面的极小风压系数沿高度方向没有明显的
注:从左到右依次为西、南、东、北立面,下同。
Figure 6. Contour line of the maximum wind pressure coefficients under all wind directions. (a) Without interfering buildings; (b) With interfering buildings
图6. 全风向角下极大风压系数分布等值线图。(a) 无周边建筑物;(b) 有周边建筑物

Figure 7. Contour line of the minimum wind pressure coefficients under all wind directions. (a) Without interfering buildings; (b) With interfering buildings
图7. 全风向角下极小风压系数分布等值线图。(a) 无周边建筑物;(b) 有周边建筑物
变化规律。各立面在靠近拐角区域的极小风压系数的绝对值较大,而在靠近立面中心区域则有逐渐减小的趋势,这主要与拐角区域气流的加速分离导致脉动风压增加有关。根据我国荷载规范的建议 [1] ,可估算出150 m~200 m高度范围内幕墙立面中心区域的极小风压系数在−1.45~−1.47之间,拐角区域的极小风压系数在−2.03~−2.06之间。将试验结果与规范建议值进行对比分析可知,我国荷载规范对立面中心区域的极小风压系数的取值偏于不安全,对拐角区域的取值则偏于保守。
有周边建筑物干扰时,从图7(b)可以看出各立面的极小风压系数的绝对值在下部区域略有减小,这可能是由于受到周边建筑物的遮挡作用。靠近中上部的极小风压系数相对无周边建筑物是变化不大,仅在局部区域的极小风压系数的绝对值有所增加,这可能跟上游建筑物的尾流干扰有关。总的来看,各立面的极小风压系数的绝对值均要大于相同位置的极大风压系数,这表明极小风压是幕墙抗风设计的控制压力。
3.2. 等效静力风荷载
等效静力风荷载是风工程研究人员和结构工程师之间的纽带 [8] 。对于超限超高层建筑,结构工程师通常需要利用风洞试验提供的等效静力风荷载对主体结构的强度和刚度进行验算,并形成超限审查报告提供给全国超限高层审查委员会进行审查。图8给出了有周边建筑干扰时50年一遇基本风压作用下钦州北部湾中心主楼基底弯矩等效静力风荷载随风向角的变化规律。由于荷载规范未能直接给出钦州地区的基本风压,因此风荷载的计算参照钦州气象局提供的气象资料 [9] 取50年一遇基本风压为0.54 kN/m2。结构的动力特性由设计部门提供,第一阶振型为X向平动,自振频率为0.375 Hz;第二阶振型为Y向平动,自振频率为0.426 Hz;第三阶振型为扭转,自振频率为0.490 Hz。考虑结构为钢筋混凝土核心筒结构,结构阻尼比取0.05。可以看出:1) 在X、Y以及Z三个方向上,最大和最小等效静力风荷载均与相应方向上的平均风荷载随风向角的变化规律一致,这说明总的等效静力风荷载主要由平均风荷载控制,脉动风荷载的贡献较小。2) Mx、My随风向角的变化规律与三角函数随角度的变化规律类似。由于建筑物的截面为方形,Mx的等效静力风荷载绝对值最大发生在270˚风向角,数值为2.41 × 109 N,My的等效静力风荷载绝对值最大发生在180˚风向角,数值为2.43 × 109 N,两者数值相近。3) Mz的数值大多为Mx或My的5%以内,但在个别角度由于脉动风的不对称性加强,Mz的等效静力风荷载绝对值也能达到相应角度下Mx的27%,需引起结构工程师的重视。
对于外形规则的高层建筑,结构设计师可根据荷载规范中对相关系数的规定,结合商业结构计算分析软件预先计算各楼层的风荷载。考虑到风洞试验提供的Fx(180˚)和Fy(270˚)最不利结果均发生在正交方向,与结构计算软件计算风荷载时一致。图9给出了50年一遇基本风压作用时最不利风向角下X和Y方向楼层风荷载试验结果与有限元结果的比较。可以看出,1) Fx的试验值随楼层高度的变化规律在大部分楼层均与PKPM值吻合较好,数值比PKPM结果略微偏大;由于顶部三维流和风洞底面的影响,顶部和底部试验结果与PKPM值差异较大。2) Fy的试验结果随楼层高度的变化规律与Fx基本一致,但相对Fx而言,Fy的试验结果比PKPM结果大,这与周边干扰有关。
3.3. 风致加速度
与结构的安全性相比,超限超高层建筑在风荷载作用下的使用性更加引人关注。我国《高层建筑混凝土结构技术规程》 [10] 规定,高度不小于150 m的高层建筑结构应具有良好的使用条件,要满足相应的舒适度要求。即在10年一遇的基本风压作用下计算的顺风向和横风向楼层顶部加速度:住宅和公寓不超过0.15 m/s2,办公楼和旅馆不超过0.25 m/s2。图10给出了10年一遇基本风压作用下由风洞试验得到的顶部楼层风致加速度随风向角变化曲线。可以看出:1) X向加速度最大值发生在90˚和270˚风向角,Y向加速度最大值发生在0˚和180˚风向角,均表现出明显的横风向效应。2) 扭转向加速度相对X和Y向较小,这与建筑的截面为方形对称截面有关。3) X向加速度最大值为
0.123 m
/s2,Y向加速度最大值为0.109 m/s2,均满足高规的限值要求。
(a)Mx (b)My
(c)Mz
Figure 8. Variations of equivalent static overturning moments with wind direction
图8. 基底弯矩等效静力风荷载随风向角的变化规律
(a) 180˚风向角下Fx (b)270˚风向角下Fy
Figure 9. Variations of the local wind forces with elevations under the unfavorable wind direction
图9. 最不利风向角下楼层风荷载随高度变化曲线图

Figure 10. Variations of wind-induced accelerations in the top floor with wind direction
图10. 顶部楼层风致加速度随风向角变化曲线图
4. 结论
风洞试验是评估超限超高层建筑风荷载的常有效方法。本文以钦州北部湾中心主楼为研究对象,详细研究了该方形截面超高层建筑表面极值风压系数、主塔楼风荷载及风致加速度的特性。可以得到以下结论:
1) 极值风压为幕墙设计的控制风压。荷载规范对立面中心区域的极小风压系数的取值偏于不安全,对拐角区域的取值则偏于保守。
2) 基底弯矩等效静力风荷载由平均风荷载控制。Mx、My随风向角的变化规律与三角函数随角度的变化规律类似。由于脉动风的不对称性,个别角度Mz的等效静力风荷载绝对值也能达到相应角度下Mx的27%,需引起结构工程师的重视。
3) 对于方形截面的超限超高层建筑,PKPM可以较好地估算楼层风荷载,但周边建筑物的存在也会导致PKPM计算的结果偏小。因此,有必要对通过风洞试验来评估周边建筑物的干扰作用。
4) 水平方向加速度表现明显的横风向效应,方形截面的扭转加速度相对水平方向较小。
基金项目
国家自然科学基金资助(51508183),湖南省教育厅开放基金资助(15K044)。
*通讯作者。