屋盖静力风荷载的数值模拟研究.docxVIP

屋盖静力风荷载的数值模拟研究 1 大跨度复合屋盖结构风荷载模拟研究 上海旗忠森林体育城的网球场中心是一个旋转的开放钢屋顶。屋顶由八瓣大型开组合屋顶（或叶片）组成，像上海花和白玉兰。屋盖向内悬挑长度61.5m,屋盖最大平面直径123m,水平投影38000m2,檐口高度34.8m。“花瓣”单片结构自重200多吨,外形长71m,宽46m,高7m,坐落在重1800多吨的环梁上。环梁与每单片屋面之间设置三条导轨和机械传动装置来完成开合。此种“花瓣”旋转的开合屋盖形式国际上也不多见。 开合形式的屋盖是一种新型的建筑结构形式,近几年得到了很大的发展,现已成为现代体育建筑的主要发展趋势之一。对于大跨度开合屋盖,其最不利状态可能是完全开启或完全闭合状态,也有可能是中间的某个状态,所以工况较多,风荷载分析显得更为重要且比较复杂,我国对此类结构的风荷载研究较少,国内外相关的参考资料也很少。 到目前为止,国内外对大跨度屋盖的抗风研究,大多基于风洞试验结果的分析与研究,也有一些通过数值模拟的手段。Borri通过风洞试验研究了罗马奥林匹克体育场环状屋盖上不同测点的功率谱;Jiming Xie进行了一座大悬挑体育场的气动弹性模型风洞试验,通过测得的加速度响应获得了各个模态的惯性力;在开合屋盖研究方面,Irwin和Wardlaw对蒙特利尔奥林匹克运动场进行了弹性模型风洞试验,得到了平均风速、湍流密度以及功率谱,并对附加质量影响进行了研究;在数值风洞模拟方面,齐辉等对益阳体育场大悬挑屋盖平均风压风布进行了数值模拟,同时给出了附近风流场的规律性。 本文首先对上海旗忠体育城网球中心开合屋盖结构进行了平均风压的数值模拟,得到了屋盖上平均风压的数值模拟结果,并与风洞试验结果进行了比较。本文还在数值风洞模拟以及计算出静力风荷载和动力风荷载基础上,运用ANSYS有限元平台,重点分析了屋盖一个叶片导轨支撑点处在不同开合状态、不同风向下支反力变化的规律性,具有参考价值;所给出导轨上四个支撑点处在不同状态下的支反力和最不利工况,可为导轨抗风设计提供依据。 2 打开地板的平均风压和动力风负荷 2.1 湍流物理模型及数值模型 由于屋盖的开合过程有许多状态,经与设计单位商定,建模时考虑有代表性的全闭合(0°)、半开启(15°、30°)和完全开启(45°)共四种状态或工况。数值风洞模拟模型由国际领先的计算流体动力学软件CFX5完成,选用两阶的SST剪切应力输运湍流物理模型,按原型尺度建模,数值模型见图1(开启15°状态)。 数值风洞入口按B类地面粗糙度类别设置风剖面,平均风速10m以下为v0=29.67m/s,10m以上为v(z)=v0×(z/10)0.16;湍流强度参考国外有关资料输入,即5m以下为I0=0.23,5m以上为I(z)=0.094×(z/350)-0.21,这里z为离地高度;入口处的平均风剖面与风洞试验一致。出口设置相对静止压力为0.0(Pa)。建筑物壁面以及地面设置为无滑移的光滑壁面,数值风洞的侧面和顶面设置为自由滑移壁面。 2.2 风压值的平均模拟与风洞试验结果的比较 2.2.1 表吸力试验 由于篇幅的原因,本文仅给出屋盖在15°开启状态、0°风向时的结果,其上表面整体平均风压的等高线分布如图3所示,其中正号代表压力(指向表面),负号代表吸力(离开表面),单位为Pa,风向规定如图2所示。 从图3中的结果可以看出: (1)屋盖绝大部分区域为负压,少数部位出现正压。由于每一个屋盖叶片由六个小面组成,叶片大体中部高,周边界低,因此明显迎着风向的小面为正压。 (2)平均风压等高线比较密集的地方位于六个小面交界线附近,其压力最大值和最小值都在这些区域。其中正压(压力)最大值为981Pa,负压(吸力)的最小值达到了-1856Pa,但分布的区域很小,一般分布在屋盖叶片两顶点的脊背交接线附近。 2.2.2 叶片升力值对比 本开合屋盖结构的各个叶片是通过沿轨道滑行进行开合状态的变化,屋盖叶片在各个开合状态下的升力值是结构设计者关心的一个重要参数。通过数值风洞模拟,我们得到了屋盖在不同开合状态下叶片各个区以及整个叶片的总升力值,并将叶片总升力值与风洞试验相应结果进行比较。下面给出图2中布置有风洞试验测点的一个叶片在0°风向时的升力值,见表1。升力值为负号表示向上的吸力,叶片上表面A、B、C、D、E、F六个区的位置及编号亦见图2。 从表1中可以看出,叶片A区、B区、F区的升力值在全闭合状态下最大,E区的升力值在屋盖开启15°的状态下最大,C区、D区和屋盖下表面的最大升力值为完全开启(45°)状态。就整个叶片来讲,在0°风向下、屋盖完全开启(45°)时的升力值最大。 由表1中的数值模拟与风洞试验升力值比较看出,除屋盖开启45°(数值模拟方法大21%)外,其余三种状态的差别不超过10%,数