基于cfd的房屋双坡房屋风压系数研究.docxVIP

基于cfd的房屋双坡房屋风压系数研究 低土地房屋是工业和民用建筑中常见的房屋形式。部分低矮房屋,如目前开发应用的冷弯薄壁型钢结构房屋及其屋面材料向着轻质高强的方向发展,对于该类结构体系,风荷载在结构设计中往往起着主要的控制作用。历次的台风灾害调查表明屋面破坏是低层双坡屋面房屋的主要破坏形式之一。1969年的“卡迈尔飓风”席卷了美国整个密西西比州之后,一美国海军基地对300幢低矮建筑物的调查后发现,有12%全部倒塌,78%需要加盖新屋面。影响低层房屋屋面风荷载的因素很多,如风向角、屋面坡角、房屋长宽比、房屋高宽比等,但是我国现行《建筑结构荷载规范》(GB50009-2001)中,仅给出了考虑屋面坡角的双坡屋面房屋各面的体型系数,对于其它影响因素均未提及,这已明显不能满足日异月新的工程应用需求。因此,较全面准确地掌握这类房屋的屋面风压分布规律十分必要。 大气边界层中钝体绕流是非常复杂的,对于结构的风荷载,很难从流体力学理论上直接进行分析。传统的研究方法有两种:现场实测和风洞试验。相对于高层建筑,低层房屋的全尺寸场地试验开展得较多。例如Hoxey等人[4―5]对12个不同山墙形状的双坡屋面房屋进行了实尺度模型风压测试,结果显示房屋的几何外形如高、宽、坡度等因素均对其受风作用有较明显的影响。在低矮建筑的原型试验中颇具影响的是美国德州理工大学风工程研究现场试验室(Texas Tech Wind Engineering Research Field Laboratory,WERFL)的TTU(Texas Tech University Building Model)建筑模型场地试验[6―8]。风洞试验是目前公认较为准确的确定结构风荷载的方法。但是,不论是风洞试验,还是现场实测都存在着试验经费大、周期长等问题。随着计算机软硬件水平的飞速发展和计算流体动力学(CFD,Computational Fluid Dynamics)技术的不断完善,出现了与试验相对应的数值模拟方法,并已逐步成为继风洞试验后预测建筑物表面风压、周围风速和湍流特性的一种新的有效方法。 王辉[10―11]采用标准k-ε和RNG k-ε湍流模型,对一幢低层双坡屋面房屋在不同风向角下的表面风压进行了数值模拟,并与风洞试验结果作了比较分析。但是,该文的风洞试验和数值模拟均忽略了大气边界层风速的影响,来流风取为均匀流,且仅对单一建筑进行了分析,并未深入研究影响其表面风压的影响因素。陈水福采用Realizable k-ε湍流模型对文献中的12个双坡屋面房屋模型的屋面风压进行了数值模拟,数值模拟结果与足尺试验的结果误差普遍在20%―30%。并在此基础上,对同类房屋在不同高宽比和不同坡度情况下的风压变化规律进行了参数分析,但是该文的数值分析也是建立在均匀流的基础上,且仅针对特定的12个足尺试验模型进行分析,部分模型可比性较差,规律不够明显。顾明对TTU标准模型进行了风洞试验研究,并采用Fluent的RSM(雷诺应力)模型和CFX的SST(剪切应力运输)模型研究了TTU实尺度模型的定常绕流场,并重点比较分析了数值模拟结果、实测结果及缩尺模型风洞试验之间的差异,并未深入研究低矮建筑的风压分布情况。文献[14―15]由标准k-ε模型得到的建筑表面风压及湍动能等时均值与大气边界层风洞试验值有一定的吻合,但在建筑物背风面及侧风面处计算值偏小,回流区域也明显偏小。 本文基于专业CFD软件平台Fluent 6.3,首先采用基于Reynolds时均(RANS,Reynolds Averaged Navier-Stokes)的标准k-ε等湍流模型计算了TTU标准低层建筑模型的静态绕流的风流场,建立了适用于低层建筑结构的数值风洞模型,并将各湍流模型的计算结果进行了比较分析。其次,采用RNG k-ε模型对低层双坡房屋进行了数值分析。较系统地研究了来流风向角、屋面坡度、挑檐长度、檐口高度和房屋长宽比对屋面风压系数以及建筑物各表面体型系数的影响。 1 ttu标准建筑模型的数值模拟 1.1 引入新的湍流模型方程 在计算风工程中,钝体绕流问题的控制方程是粘性不可压Navier-Stokes方程。当前应用最广的是基于RANS的N-S方程,并由湍流模型进行封闭,再作离散求解的数值模拟方法。湍流时均流动的控制方程如下[16―17]: 上面三式分别为时均连续方程、雷诺方程和标量φ的时均输运方程。要使方程组封闭,必须对Reynolds应力项作出某种假定,即引入新的湍流模型方程,把湍流的脉动值与时均值等联系起来。标准k-ε模型是涡粘模型中最具代表性的两方程模型,并在科学研究及工程实践中得到了最为广泛的检验和成功应用。在标准k-ε模型中,对于Reynolds应力的各个分量,假定粘性系数μt是各向同性的标量。由于建筑物背面、侧面处