伞形膜结构风荷载分布特性的数值模拟研究

2013-05-24杨洪智

中国建筑金属结构 2013年6期

杨洪智

（中国建筑西南设计研究院有限公司，四川成都610042 ）

0 引言

索膜结构作为一种新颖的结构形式，以其独特的建筑造型和卓越的轻质高强性能[2]，正得到日益广泛的应用。传统的风洞试验方法所需的时间长，费用高，而且由于某些相似条件在实验中很难满足，从而在一定程度上降低了实验精度。近年来，伴随着计算机硬件技术的迅速提升以及数值计算科学的发展[3]，在计算机上综合运用CFD数值模拟技术对结构周围绕流流场进行数值模拟已成为一种新的解决该问题的途径。与风洞试验方法相比，CFD数值模拟技术优点在于[4]：费用省、周期短、效率高，可以构造与实际结构尺寸相同的计算模型，方便地变化各种参数等。因此CFD数值模拟技术在索膜结构抗风设计中具有广阔的前景。

基于上述考虑，本文运用CFD数值模拟技术对处于大气边界层底层的伞形膜结构风荷载分布特性进行了系统研究，对其表面风压分布进行了数值模拟。数值模拟依托CFX软件平台，基于雷诺平均模拟（RANS）方法进行。

1 伞形膜结构绕流特性的数值模拟

1.1 计算模型

（1）结构建模

如图1所示，计算模型为正六边形平面伞形膜结构，底部开敞，其中结构水平投影的外接圆直径为10m，即L=10m，中心为一个半径0.5m的钢环。不设膜内索，形成水平环向曲率变化连续的膜曲面。膜厚度t=lmm，初始预张力T0=2.5kN/m，张拉刚度Et=255kN/m。索的弹性模量为160Gpa，泊松比为0.3。本文选取如下参数进行计算：

1）风向角（a）： 0°，90°；

2）矢跨比（f/L）：1/12，1/10，1/8；

图1 伞形膜结构计算示意图

膜结构有限元模型中，膜面采用shell 41单元，索采用link 10单元，假定索和膜之间没有相对滑动。本文膜结构找形采用了支座位移法和小弹性模量法进行初始找形。

（2）流域建模

根据结构尺寸，计算流域取为120m×80m×40m,膜结构置于整个流域沿流向的1/3位置处。流域设置满足阻塞率<3%的要求，以尽量消除计算域对所关心的模型附近流动状态的影响。对整个流域进行混合网格划分，即在靠近结构的区域使用较密的四面体网格，其余部分采用相对稀疏的结构化六面体网格，在四面体网格与六面体网格过渡处和结构表面附近使用棱柱体网格。流场和网格划分如图2、图3所示，其中网格节点数为6万，网格数为13万，其中非结构网格为8万，结构网格为5万。采用二阶迎风格式离散控制方程。本文的流体计算模型中流体质量密度（空气）均取为1.29kg/m3，粘度μ均取为1.74×10-5 kg/m·s。场地类型为C类，地面粗糙度指数a= 0.22。

图2 0°风向角伞形膜结构流场透视图

图3 0°风向角伞形膜结构流域网格划分

1.2 边界条件

（1）进口边界条件

式中 zb为标准参考高度，我国规范取作10m；z为任一离地高度(m)；为标准参考高度 zb对应的平均风速(m/s)；a为地面粗糙度指数。入口处采用的湍流强度 Iu、湍流动能k和湍流动能耗散率ε的表达式分别为：

（2）出口及壁面边界条件

出流面接近完全发展，采用完全发展出流边界条件(outflow)；流域顶部和两侧采用对称边界条件(symmetry)，等价于自由滑移的壁面；膜结构表面和地面采用无滑移的壁面条件(wall)，并在地面引入粗糙壁面修正。

2 伞形膜结构风荷载分布特性

2.1 风向角对膜面风荷载分布的影响

伞形膜结构矢跨比f/L =1/10，风速v=20m/s，矢跨比f/L=1/4，使用瞬态求解。来流风向分别为0°，90°，计算结构表面的风压系数分布。不同风向角作用下的膜面风压系数分布，详见图4～图5。

图4 0°风向角伞形膜表面风压系数分布图

图5 90°风向角伞形膜表面风压系数分布图

由于是开敞式膜结构，将膜面对应点的上下表面风压系数相减，得到膜面的净风压系数。计算模型和边界条件关于流向（X轴）对称，数值模拟的结果整体上也成对称分布，初步验证了数值模拟的合理性。根据不同风向角下膜表面风压系数分布图，可以看出，0°风向角时在膜面前缘有强烈的气流分离现象,膜面主要分为两个区域，迎风面为正压区，背风面大部分为负压区，并且在靠近顶部圆环处出现了一个小的正压区。90°风向角时膜面风压分布情况与0°风向角时相似，风压分布系数值稍有不同，最大负压的绝对值比0°风向角时稍大。