复杂多等离子体大跨空间建筑风荷载的数值模拟研究.docxVIP

复杂多等离子体大跨空间建筑风荷载的数值模拟研究 随着经济和社会的发展，城市地标的形状越来越多样化，大型多环形建筑的结构也越来越注重其生动的审美特征。该结构表面的折叠不规则，曲线多变，并经常干扰内各单体结构之间的风场，这使周围的风场非常复杂。目前，我国的负荷规范中很难反映复杂不规则、多因素多链结构的抗风设计和安全评价所需的结构表面风压分布和风荷载值，这是计算和抗风灾的难题。 基于计算流体动力学(CFD)的数值模拟方法是获取复杂空间建筑表面风压和风荷载的有效手段之一.在这方面,卢旦等 本文以上海世博会博物馆景观建筑“欢庆之云”为对象,研究不规则多连体大跨空间建筑结构风压分布、阵风系数和等效静风荷载.运用Rhino 3D软件构建结构几何模型,运用ANSYS/FLUENT软件开展结构表面风压、风荷载和建筑周围风场的数值模拟分析.通过平均风模拟获得结构表面风压分布,计算结构表面的平均风压 1 为物体模型的构建和理论基础 1.1 博博物馆“狂欢之云”建筑的几何模型 以上海世博博物馆“欢庆之云”为背景对象,进行风压分布数值模拟和等效静风荷载分析.该工程为空间曲面大跨度网格结构,位于上海市世博会地区浦西文化博览区,建筑效果图如图1所示. 将“欢庆之云”建筑分为3个柱支撑的云状大跨空间建筑及其相互间的连接体建筑.运用Rhino软件,构建包含世博博物馆“欢庆之云”建筑(简称云结构)和周边建筑物在内的建筑群的几何模型,3个云结构编号依次为1 号云、2 号云和3 号云(见图2(a)).鉴于建筑形体的复杂性,在计算结构风压分布时,进行建筑结构表面的区域划分.“欢庆之云”建筑上表面称为0 号分域,具体包括建筑上表面的0-1区、0-2区、0-3区和0-4区,而0-5区是与0-4区相对的建筑下表面.1号云结构对应的立柱状网壳称为1号立柱网壳,1 号立柱网壳表面分为1-1,1-2,1-3,1-4,1-5和1-6等6个分区;1-2区是1-1区逆时针方向的相邻曲面,以此类推.同理划分2号立柱网壳、3号立柱网壳的分区,均划分为6个分区,它们的编号方式与1号柱相同.结构的分区如图2(b)所示. 根据“欢庆之云”建筑的尺度和体型确定风场计算域,即为950m×400m×200m.图3为风场计算域的网格划分,其中对建筑群和地面边界层位置的计算网格进行适当加密,网格数量约为2.5×10 1.2 基于rans的k-湍流模型 围绕上述复杂形体建筑结构的风压和风场问题,运用CFD方法进行模拟计算,即求解Navier-Stokes方程并引入基于RANS(Reynolds Average Navier-Stokes)的k-ε 湍流模型: 式中:ρ为流体密度;p为流体所受压力;μ为空气分子黏度;μ 根据Richards Hoxey假定,将平均流方程和上述k-ε 表达式写为 2 模拟分析平均风压值 2.1 风特性分析参数 本工程所处地区为B类地貌(粗糙度系数α=0.16),取上海地区50 年重现期基本风压为0.55kPa.入口边界采用指数风剖面,选用基于雷诺时均的标准k-ε2方程涡黏湍流模型,以湍流动能k、湍流耗散率ε2个参数来表征风的湍流特性: 式中:I为湍流强度;l为湍流积分尺度;v为风速. 出口边界条件采用完全发展的出流边界条件;数值风洞顶部和两侧墙面采用自由滑移的壁面条件;结构表面和地面采用无滑移的壁面条件. 气体选用不可压缩的常密度空气模型.计算中采用3D单精度隐式算法,分离式求解器,流项的离散格式选择精度较高的二阶迎风格式,速度-压力耦合选取SIMPLEC算法. 2.2 平均风压分布 以风向角θ= 0°为例进行建筑表面风压分布规律的分析,风压云图如图4所示.图4(a)所示为θ=0°时结构上表面平均风压p 图4(b)所示为θ=0°时结构下表面平均风压分布云图.由图可见,下表面为负风压且其变化较为平缓;对于立柱网壳,在面向来流方向的表面负风压较大,而其背面负风压较小;2号立柱网壳与3号立柱网壳之间迎风面气流分离明显,造成很大的负风压,而在背风面负风压很小;其他的来风风向角下结构风场风压分析方法相同,不再赘述.综上可见,由于该大跨度空间建筑结构的形体多变、竖向起伏且存在建筑物遮挡,导致结构表面的风场在竖向、水平向呈现出复杂的分离、再附现象,结构表面风压分布十分复杂. 2.3 分区表面负风压 平均风压系数C 对云结构的上表面和下表面(见图5(a)),气流在迎风面产生分离,其上、下表面均承受负风压,如0-5区处于云结构的下底面,曲面形状复杂,平均风压系数随风向角变化剧烈.对1号立柱网壳(见图5(b)),在θ=135°~180°下直接迎风,其几个分区表面承受正风压;在θ=0°~90°及θ=225°~315°下迎风面的上游有结构物遮挡,其6个分区表面承受负风压且风压分布相似.对2 号立柱