摘要:以广州国际会展中心为工程案例,进行了刚性模型风洞试验和有限元模态分析,在此基础上计算了屋盖的风致位移响应,结果表明,对于屋盖风振响应影响最大的因素是结构的基阶振型,其次才是风荷载;随着阻尼比的增加,位移响应谱的共振峰得到了有效平抑,峰值位移响应也随之减小。基于现场实测,提出功率谱点积识别大跨屋盖结构竖向整体振动固有频率的方法,并与有限元计算结果进行了对比分析,实测得到的前4阶固有频率与振型均能与有限元模型较好吻合,其中基阶固有频率的相对误差仅为0.3%,证明了该方法的有效性。
关键词:大跨屋盖 风洞试验 现场实测 模态识别
Wind tunnel and full-scale study of wind effects on a large roof structure
Jiyang Fu1, 2, Ruohong Zhao2, An Xu2, Jiurong Wu2
(1. School of Civil Engineering, Harbin Institute of Technology, Harbin 150001, China)
(2. Guangzhou University - Tamkang University Joint Research Center for Engineering Structure Disaster Prevention and Control, Guangzhou University, Guangzhou 510006, China)
Abstract: In this study, a wind tunnel test was conducted for the Guangzhou International Exhibition Centre. Simultaneous pressure measurements were made on its rigid roof model in a boundary layer wind tunnel. Meanwhile, a three-dimensional finite element model of the roof structure of Guangzhou International Exhibition Centre was established based on the structural design drawings; and the dynamic characteristics were analyzed based on the established finite element model. Furthermore, wind-induced responses of the large roof were computed; in which the wind loading inputs were derived from the simultaneous pressure measurements on the rigid roof model. The results showed that the influences of the structural fundamental modes on wind-induced vibration of the roof were generally more significant than wind loads. It was also observed from the wind tunnel measurements that there was a tendency for the peak displacement response to decrease with the increase of the damping ratio. On the other hand, full-scale measurements of wind effects on the large roof were conducted under strong wind action. Based on the field measurement results, a new method to identify the first several natural frequencies of the roof in vertical direction, using the dot matrix of power spectrum density approach, was presented. The frequency and model shape results determined by the proposed method and those obtained from the finite element model analysis are in good agreement with each other, and the relative error of fundamental frequencies obtained by the above two method is only 0.3%, which illustrates the applicability and accuracy of the proposed method.
Key words: wind tunnel test, full-scale measurements, wind effects, large roof
1 引言
大跨度屋盖结构具有质量轻、柔性大、阻尼小、自振频率低等特点,因而风荷载成为其结构设计的控制荷载。这类结构往往比较低矮,在大气边界层中处于风速变化大、湍流度高的区域,且屋顶形状多不规则,其绕流和空气动力作用十分复杂,所以这种大跨屋面对风荷载十分敏感,尤其是风荷载作用下的动态响应值得研究[1][2]。
本文以广州国际会展中心这一超大跨屋盖结构作为工程案例,进行了刚性模型风洞测压试验,结合有限元模态分析的结果计算了屋盖结构的风振响应并进行了分析;采用拾震器现场实测了屋盖结构在风场作用下的竖向振动速度时程,在此基础上提出功率谱点积法,识别了屋盖结构竖向整体振动的固有频率,并结合自互谱法识别了屋盖结构前四阶振型。上述分析结果与有限元模型的模态分析结果能够较好吻合,证明了有限元建模的合理性。
2大跨屋盖风振响应的风洞实验研究
2.1 大跨屋盖风振响应的计算原理
大跨屋盖的风振响应计算实质上是多自由度体系在时变外力作用下的响应问题:
(1)
其中, 、 、 为结构的质量阵、阻尼阵和刚度阵, 为屋盖有限元模型的自由度数, 为风洞试验中的测得的风压,此处假定其已转换为无量纲的脉动风压系数, 为测压点个数。
由完全二次型相关法(complete quadric correlation, CQC)可得结构位移响应谱的计算公式为[3]:
(2)
其中, 为振型矩阵, 为频率响应函数矩阵。上式考虑了所有振型交叉项及振型间的耦合,是计算大跨屋盖风振响应的精确算式[4]。针对CQC法计算量巨大的问题,本文采用了文献[5]提出的一种简化的算法——谐波激励法,其计算结果与上式完全等价。
2.2 广州国际会展中心风洞试验简介
广州国际会展中心位于广州东南部的琶洲岛,总用地面积70万m2,一期工程用地面积41.5万m2,建筑面积39.5万m2,建筑物总长度为458.5 m2,总宽度318.2 m2,包括六大功能区:展厅(含南面展厅和北面展厅)、室外展场、珠江散步道、车道、东入口车道及其他用房。广州国际会展中心风洞试验在汕头大学大气边界层风洞试验室进行,地貌类型按照国家《建筑结构荷载规范》(GB50009-2001)中的B类地貌考虑,地貌粗糙度指数 。刚性模型缩尺比为1:300,如图1所示。测点布置如图2所示,虚线所包含的范围为结构的E跨,该区域结构独立,本文将其作为风振响应的研究对象,因此在该区域上布置了更为密集的测压点,以获得更为全面的风荷载信息。风洞试验采样频率为312.5Hz,采样长度为20480帧。从0˚到360˚逆时针旋转每10˚为一个测试风向
2.3 风振响应分析
广州国际会展中心E跨屋盖采用ANSYS软件进行有限元建模和模态分析。其有限元模型如图3所示。图4为其前4阶固有频率及其对应的振型。从图中可知,结构第一阶振型主要表现为屋盖整体以半波长正弦波形式在竖向平面内振动,第二阶模态为屋盖的左、右半部分别以奇偶正弦波交替出现振动。各阶振型以竖向位移为主,因此本文的实测研究也针对屋盖的竖向振动。需要说明的是本文的风振响应计算中选取了前25阶振型,此处仅列出前4阶。
根据式(2)及上图给出的有限元模态分析结果按照1%阻尼比计算屋盖表面各节点的风致位移响应,全风向下峰值位移和均方根位移响应等值线分布如下图所示
由图5可见,跨中右侧为位移响应最大的区域,这是由于跨中区域在基阶振型中处于位移最大的位置,同时跨中右侧在0°风向角附近处于迎风面,所受的风荷载较大。图5显示在全风向下E跨屋盖角部风压绝对值较大,但图4表明屋盖角部的风致位移响应却较小,这表明对于这种外形相对规则的桁架梁结构屋盖,其风致振动是由其一阶振型所控制的,对风致位移响应影响最大的因素是结构的基阶振型,其次才是风荷载。为了考察阻尼比对计算结果的影响,将阻尼比取值分别设为0.005、0.01、0.02及0.05,分别计算在0o风向角下6榀桁架梁跨中节点(图3所示的1~6号节点)的位移响应。不同阻尼比对均方根位移响应计算结果的影响如表1所示。
由表1可见,阻尼比取值对均方根位移响应的影响比较大。大跨屋盖节点位移响应标准差随着阻尼比的增大而显著减小。其原因在于阻尼比的增大对于位移响应谱的共振峰有着明显的抑制作用,以1号节点为例,图6显式了不同阻尼比取值工况下的位移响应谱:
3 大跨屋盖结构模态参数的识别
针对大跨屋盖结构的风致效应实测研究目前多集中于风荷载特性方面,通过风致响应实测数据识别大跨屋盖结构的模态参数目前仍较少开展。这一方面是由于大跨屋盖在台风作用下的振幅通常远小于超高层建筑,导致其风致响应的实测数据信噪比难以达到较高的水平;另一方面大跨屋盖风致振动实测研究需要多个传感器同步采样,线路布设的施工通常有一定困难。
以广州国际会展中心E跨屋盖作为研究对象,进行了风致振动的现场实测,在此基础上识别了结构动力参数[3]。该跨由6榀跨度为126.6米的单向张弦梁平行布置并通过刚度很大的檩条及水平支撑构成。为了获得结构完整的振型,需要布置足够多的测试点,但由于传感器和信号采集仪通道数的不足,难以做到同时记录所有测试点上的振动信号,故而不得不采取分组分批测试的策略,这也导致了数据处理量和复杂性的增加。为了使分组测试的各组振动信号具有可比性,必须在所有测试点中选择一个作为参考点,参考点被包含在每一组中,重复测试。由于测试点的振幅越大,其信噪比就越高,因此,参考点的选择应注意避开各阶振型的驻点,也即使参考点在各阶模态振动下均有较大的振幅[6]。
本现场实测采用中国地震局生产的941B竖向拾振器,采样数据格式为竖向拾震器测得的屋盖竖向振动速度时程,单位为m/s。由于现场条件以及仪器条件的限制,拾振器不能直接安装在桁架上端,而只能安在两榀桁架之间的主檩条中间。由于屋面主檩条的竖向刚度较大可以假设主檩条能与两端的桁架共同振动。鉴于上述原因,在本文的测试中,首先选用一个拾振器作为参考测试点安装在隔壁马道的主檩条与马道连接处(即图7所示的32号测试点处),拾振器与主檩条固定在一起,以保证二者的振动的一致性。另外三个拾振器作为移动拾振器分别放在事先安排好的测试点处,与参考点的拾震器组成一个测试组。全部测试点及分组安排如图7所示。整跨屋盖共采用20组测试,共60个测试点。采用优泰软件的数据采集系统对每组进行同步数据采集,采样频率为25.6Hz,每组采集时间大约为6min。
已有的研究表明,第一阶固有频率和振型对于规则体型大跨屋盖结构的风致振动起重要作用。因此,实测数据对于屋盖结构第一阶固有频率的识别结果能否与有限元模型的计算结果较好吻合尤为重要。自互谱法[7]是识别结构固有频率和振型的常用方法,并考虑到跨中位置的信号信噪比较高,因此,本文首先考查跨中测点的竖向振动速度的自功率谱。以图7所示的8、22号测点为例,其实测竖向振动速度功率谱如下图所示:
由图8可见,各测试点的速度谱的首个峰值出现的频率位置基本相同,而通过有限元模态分析已知屋盖的第1阶模态振动为屋盖整体沿桁架纵向的半波长正弦波振动。这表明,各测点基本能够识别屋盖结构的第1阶固有频率,但由图8可见,2个测点的功率谱图形包含了部分的噪声信号,以至于图形中出现了一些并非结构固有频率所对应的峰值。文献[8]提出了一种基于相关函数的识别方法,但由于该方法涉及到矩阵的奇异值分解,仍较为复杂。
本文通过求功率谱点积的方法,可以简单快捷的判断屋盖竖向整体振动的前几阶固有频率。设第i个测试点的功率谱为 ,由于实际工程的采样信号是离散信号,因此其功率谱也是对应于频率点的离散值,其序列记为向量,同理第k测试点的功率谱向量记为 ,则上述两个向量的点积即为这两个测试点的功率谱点积。前述5个位于跨中位置的测试点的功率谱点积如图9(a)所示。功率谱点积在1.005Hz处出现了大于其余位置若干个数量级的峰值,非常清晰的显示了结构整体振动的第一阶固有频率。之所以出现这种现象,是由于功率谱的点积对于各榀桁架整体振动的谱峰位置有放大效应,而非整体振动的谱峰以及噪声信号所产生的“假峰”则被“抹平”了。因此通过功率谱点积来辨识结构整体振动固有频率,其精确度和可信度均高于通过观察一个或几个测点功率谱的方法。采用点积方法辨识得到的屋盖结构第一阶固有频率为1.005Hz,而有限元模型模态分析的结果为1.002Hz,两者相差不超过1%。在识别结构的第1阶固有频率后,其余各阶固有频率可通过令第1阶固有频率谱峰值为零再绘制功率谱图形的方法来识别,如图9(b),采用这种方法依次可得到屋盖结构的第2、3、4阶整体振动固有频率为1.375Hz、3.225Hz、4.025Hz。有限元模型模态分析的结果与实测结果的比较如表2所示。
结果表明,有限元模型模态分析和实测得到的结构第1阶固有频率误差仅为0.3%,这既说明了有限元建模的合理性,也说明本文提出的功率谱点积法是可行的。两种方法得到的第2、3阶固有频率的误差均不超过10%,第4阶固有频率的误差较大,其原因有待进一步研究,但基阶模态振动对屋盖总体风振响应的贡献约占70%,而第4阶模态对于屋盖振动的贡献很小。综上所述,有限元模型的模态分析结果于实测吻合较好,采用有限元模型结合风洞测压试验计算得到的屋盖响应是可信的。进一步根据自互谱识别结构的前4阶振型如下图所示:
由上述振型识别结果可知,第1、2、3、4阶振型分别为屋盖沿整跨0.5波长、1波长、1.5波长和2波张的正弦波竖向整体振动,这与有限元模态分析的结果是一致的。
4 结语
本文以广州国际会展中心E跨屋盖为案例,对其进行了基于风洞试验的风振响应分析和竖向风致振动的实测研究,得出以下基本结论:(1)对于类似于广州国际会展中心E跨屋盖的张弦梁结构,其风致振动通常为基阶振型所控制,因此跨中位置的峰值位移响应通常大于其他位置;(2)阻尼比取值对于结构均方根位移响应的计算结果有着较大影响,阻尼比的增大使得位移响应谱的共振峰被平抑,从而降低了结构的均方根位移响应;(3)功率谱点积法有效避免了自互谱法识别固有频率过程中所遇到的“假峰”或“毛刺”问题可以方便快捷的识别大跨屋盖结构竖向整体振动的固有频率。进一步采用自互谱法识别了结构的振型,结果与有限元模型的计算结果基本一致,证明了该方法的有效性。
参 考 文 献
[1]顾明, 周晅毅, 黄鹏. 大跨屋盖结构风致抖振响应研究[J]. 土木工程学报. 2006,27(2):37-42.
GU Ming, Zhou Xuan-yi, Huang Peng. A study on the wind-induced buffeting responses of large-span roof structures[J]. China Civil Engineering Journal. 2006,27(2):37-42.
[2]Uematsu Y, Watanabe K, Sasaki A, et al. wind-induced dynamic response and resultant load estimation of a circular flat roof [J]. J. Wind Eng. Ind. Aerodyn., 1999,83:251-261
[3]谢壮宁, 倪振华, 石碧青. 大跨度屋盖结构的等效静风荷载[J]. 建筑结构学报. 2007,28(1):113-118.
XIE Zhuang-ning, Ni Zhen-hua, Shi Bi-qing. Equivalent static wind loads on large span roof structures[J]. Journal of Building Structures. 2007,28(1):113-118.
[4]李方慧, 倪振华, 沈士钊. 大跨屋盖结构等效静风荷载研究[J]. 工程力学.2007,27(7):104-109.
LI Fang-hui,NI zhen-hua,SHEN Shi-zhao. Equivalent static wind loads on long span roofs[J]. Engineering mechanics. 2007,27(7):104-109.
[5]谢壮宁. 风致复杂结构随机振动分析的一种快速算法——谐波激励法[J].应用力学学报. 2007, 24(2): 263-267.
XIE Zhuang-ning. New rapid algorithm for wind-induced random vibration of complex structures[J]. Chinese j ournal of appl ied mechanics. 2007, 24(2): 263-267.
[6]朱乐东. 桥梁固有模态的识别[J].同济大学学报. 1999,27(2):179-183
ZHU Le-dong. Modal identification of bridges[J]. Journal of Tongji University. 1999,27(2):179-183.
[7]李国强, 陈素文, 李杰等. 上海金茂大厦结构动力特性测试[J].土木工程学报. 2000,33(2):35-39.
LI Guo-qiang, CHEN Su-wen, LI Jie, et,al. Measurements of dynamic properties of Shanghai Jinmao Building[J]. China Civil Engineering Journal. 2000,33(2):35-39.
[8]李中付, 宋汉文, 华宏星等. 一种白噪声激励下模态参数辨识方法[J].振动工程学报. 2002,15(1):52-56.
LIZhong-f u, SONG Han-wen, Hua Hong- xing, et,al. A method of modal parameter identification under white noise excitation[J]. Journal of Vibration Engineering. 2002,15(1):52-56.