高层建筑绕流场气动力效应的数值模拟

发布时间：2020-08-13 21:16

【摘要】：对高层建筑抗风的研究主要依靠风洞试验技术,近年来随着计算机技术和计算风工程理论的发展,高层建筑在风荷载作用下的数值模拟逐渐取得了一定的进展。计算流体动力学已经发展成为风工程应用中传统方法(风洞实验和现场测量)的补充工具。为了研究高层建筑三维绕流场的风荷载,本文采用雷诺平均和大涡模拟的方法对刚性结构表面的风荷载进行模拟,主要工作如下:(1)对CFD的基本问题——控制方程进行研究。基于有限体积法采用FORTRAN编制SIMPLE算法程序求解不可压缩流体的N-S方程;引入半交错网格,改进算法,采用两种不同的算法来计算N-S方程,对经典案例——二维稳态方腔流进行计算,验证了两种算法的可行性和程序的正确性;采用保正格式的有限体积法,针对四种高度扭曲的网格计算扩散方程,经过与基于有限体积法的SIMPLE算法,以及FLUENT进行对比发现,此方法计算结果稳定且能够保正;基于以上对控制方程的研究,分别确定使用雷诺平均方法(RANS)和大涡模拟方法(LES)时,两种方法的离散方程的方法,扩散项、对流项、时间项的离散格式,计算参数(时间步长、各变量松弛因子等),并为FLUENT自带UDF文件的编制提供基础。(2)数值模拟得到可靠结果需要稳定的入流边界条件(即满足特定条件的大气边界层风速场)和合理的网格。本文对目标边界层风洞的B、D类地貌的平均风剖面,湍流强度,湍流积分尺度的数据进行拟合,确定生成风场的基本参数;MATLAB自编程序实现CDRFG方法,对B类和D类地貌,有效生成满足Von-Karman谱的风速场;对CDRFG方法进行改进,显著提高计算效率,节约计算资源。网格的合理性对数值计算结果影响非常大,本文采用ANSYS ICEM划分网格,研究网格的第一层网格尺寸y+,网格延伸率等对计算结果的影响,经过反复调试,试算,最后生成了符合计算要求的网格,通过监测出流面的速度场验证网格的自保持性。(3)采用雷诺平均k-ε模型对高宽比为9的超高层建筑的表面平均风压进行模拟,内容包括几何建模,网格生成,设置边界条件,计算参数设定,通过FLUENT自带UDF程序将生成的风速文件加载到入流的速度边界,完成初始条件设定、边界条件设定、逐步控制计算过程并保存计算结果,将雷诺平均数值模拟的计算结果与风洞试验结果进行对比;为进一步验证雷诺平均方法的有效性,对置于高层建筑顶部的光伏板组件表面平均风压系数进行数值模拟,研究建筑顶层女儿墙和光伏板距对光伏板表面风荷载的影响。(4)采用大涡模拟Samagorinsky湍流模型针对高宽比为9的超高层建筑模型,分别在B类地貌和D类地貌下的风压进行数值模拟,并与雷诺平均和风洞试验的结果进行对比,得到B、D类地貌下超高层建筑表面风荷载的特性;分析计算结果得到结论,虽然需要更加精细的网格和更多的计算资源,采用大涡模拟技术研究超高层建筑的风荷载特性具有很好的稳定性和优越性,相比较风洞试验,数值模拟能够更加详细地了解结构表面的流动情况;将本文中高宽比为9的高层建筑与原有文献资料中相同截面尺寸,高宽比分别为4、5、6的高层建筑背风面风压系数进行对比,最终结果显示高宽比为9的高层建筑三维绕流效应是最弱的。
【学位授予单位】：武汉大学
【学位级别】：博士
【学位授予年份】：2018
【分类号】：TU973.2
【图文】：

斯托克斯方程,流动控制方程,物理问题,N-S方程

和技术的发展，计算流体动力学（Computational展成为风工程应用中传统方法（风洞实验和现场被人们接受应用于空气污染物扩散和复杂地形上住房规划，和结构设计等（Yan B. and Li Q.,20程（ComputationalWindEngineering，简称 CW维斯托克斯方程（Navier-Stokes Equations，简连续问题用离散的格式来近似的计算方法。为了近似值，首先要以给定微分方程或者基本物理定比如网格节点）未知量近似值之间的代数方程，散，然后解方程。目前的研究中数值计算方法、想都可以归结为如图 1-1。把原来在时间、空间场、温度场等），用有限个离散点上的值来代替离散得到关于这些值的代数方程，求解之，最后铨，2001）。

示意图,动量方程,物理意义,示意图

图 2-1 动量方程物理意义示意图) ( ) ( ) ( )( )yxxx zxx uu uv uw pfx y z x x y z ) ( ) ( ) ( )( )xy yy zyy vu vv vw pfx y z y x y z ) ( ) ( ) ( )( )yzxz zz wu wv ww px y z z x y z 体，切应力与速度的梯度成正比，代入上式可得：2) 1( ) ( )3x puV V u fx x 2) 1( ) ( )3y pvV V u fy y 2) 1( ) ( )3zw pwV V u f z z 缩流体，其连续性方程可进行以下推导：0 0u v w

示意图,能量方程,物理意义,示意图

超高层建筑绕流场气动力效应的数值模拟根据能量定律，物质体内总能量增长率等于作用在物质体上的体积力做功的表面力做功、物质体内的生成热和边界面上因热传导输入的热量之和（热第一定律）。( )( ) ( ) ( )Ediv uE f u div p u q div Tt (2-7 E 为总能（内能加动能），q 为流体微团单位质量单位时间的生成热，T 为温k 为热传导系数。

【参考文献】