摘要
以中国东南沿海某超高层三塔连体建筑为研究对象,以风‑雨双向耦合算法为核心,基于计算流体动力学技术采用连续相和离散相模型进行风场和雨场的迭代模拟。首先基于9种风雨组合工况进行三塔连体建筑非定常脉动风场模拟,探讨超高层三塔连体建筑平均风压分布、表面速度流线和流场干扰机理。然后对比研究不同风雨组合工况下主塔表面雨滴附着数量、雨滴冲击力和雨压系数的分布规律,揭示风‑雨耦合场中结构表面速度流线、雨滴运行轨迹和最终速度的作用机理。最后提炼出超高层三塔连体建筑最不利风‑雨组合工况,并给出对应的雨压系数取值建议。研究表明:风‑雨耦合环境下超高层连体建筑迎风面雨荷载作用最为显著,此时雨荷载与风荷载最大比值可达23.81%,局部测点最大雨压系数达到0.301,100 a重现期风速和强大暴雨组合为风‑雨耦合作用的最不利组合工况。
相较于单体超高层建筑,超高层多塔连体建筑的相互干扰效应导致其表面风荷载分布及流场作用机理更加复
超高层建筑风致干扰效应属于结构风工程的传统研究方向,已有研究主要集中在多个单体建筑之间的静力风荷载干
关于风‑雨共同作用的研究主要集中于风单向驱动雨产生的效应,文献[13‑14]首次提出基于三维稳态风场下的雨滴轨迹数值计算方法,并得到了三维风场中建筑物迎风面风驱雨的分布情况;文献[15‑16]对斜拉桥拉索风雨激振机理及相关参数影响进行了系统研究;文献[17‑18]采用CFD数值模拟的手段研究了风驱雨作用下低矮房屋及特定群体建筑的迎风面雨压分布特性,相关研究表明降雨产生的雨压荷载在极端条件下可达到纯风压的30%以上,设计工程结构尤其是覆面结构时,应考虑极端气候下雨压所引起的荷载作用。已有研究主要集中于风驱雨作用下斜拉桥拉索、低矮建筑和输电塔等结构,对于日益兴建的超高层结构尤其是超高层多塔连体建筑,鲜有系统讨论风‑雨共同作用下主塔表面雨荷载分布规律及其与子塔之间的干扰效应。
鉴于此,基于风‑雨双向耦合算法并采用连续相和离散相模型进行风场和雨场的迭代模拟,以中国东南沿海地区某超高层三塔连体建筑为例,系统探讨不同风‑雨组合工况下主塔表面雨滴附着数量、雨滴冲击力和雨压系数的分布规律,揭示风‑雨耦合场中结构表面速度流线、雨滴运行轨迹和最终速度的作用机理。在此基础上,提炼出超高层三塔连体建筑最不利风‑雨组合工况,并给出对应的雨压系数取值建议。研究结论可为极端环境下超高层三塔连体建筑雨荷载和抗风设计提供科学依据。
降雨强度是指单位时间内的降雨量,与24 h平均降雨量相比,小时降雨量更能直观反映出极端天气条件下瞬时雨强对建筑结构安全性能的影响,
雨滴谱近似服从负指数分布,常用Marshall‑Palmer
(1) |
式中:Dp为雨滴直径;n(Dp)为不同直径雨滴个数浓度谱;N0为浓度,取常数值8 000;λ为尺度参数,其表达式为
(2) |
式中R为降雨强度。
降落高度大于20 m时几乎所有粒径雨滴都能达到最终末速
(3) |
式中:v(D)表示直径为D的雨滴垂直降落末速度。
在大暴雨时,雨滴占空气体积分数远小于10%,故可采用DPM(Discrete phase model)模型进行雨滴模拟,并在风场计算稳定后作为第二相插入连续相中进行风‑雨双向耦
(4) |
式中: up为离散相颗粒速度; u为连续相流体速度;FD(u-up)为单元颗粒质量的拖曳力;ρp和ρ分别为颗粒和流体密度;F为由离散相和连续相间的相互作用力。其中
(5) |
式中:μ为流体黏性系数;Dp为颗粒直径;Re为相对雷诺数,可表示为
(6) |
考虑雨滴离散相影响后,风连续相基本控制方程可表示为
(7) |
(8) |
式中:Sm为加入连续相的第二离散相质量,p为压力,为应力张量,为重力。其中应力张量可表示为
(9) |
式中:I为单位张量,等式右边第二项为体积膨胀作用。
雨滴冲击到塔筒壁面过程服从动量守恒定律,求解冲击力的关键在于碰撞时间。计算中忽略雨滴在冲击过程中可能发生的蒸发、飞溅、破裂等现象,认为雨滴与结构间相互作用遵循牛顿第二定律。由动量定理
(10) |
式中:f(t)为单个雨滴冲击力;v为雨滴速度。
雨滴在单位时间内对结构的冲击力F(τ)为
(11) |
将下落时雨滴近似看作球体,则
(12) |
由于雨滴直径一般在6 mm以下,且撞击前水平末速度相对较大,故为简化计算,将碰撞时间τ取为
(13) |
则雨滴对结构的冲击力可简化为
(14) |
超高层三塔连体建筑位于东南沿海地带,主体工程有3栋连体建筑,其中主塔为65层塔楼外加塔冠,总高度300 m,子塔B、C塔楼总高度158.7 m,底部为5层裙房,总高23.6 m,如

图1 超高层三塔连体建筑及测点布置示意图
Fig.1 Schematic diagram of model and measuring point layout of the three‑tower connected tall buildings
该超高层三塔连体建筑位于B类地貌,重现期10、50和100 a对应的最大风速分别为28.28、33.47和35.77 m/s,大暴雨气候条件下3种雨强分别为:弱大暴雨(64 mm/h)、中大暴雨(100 mm/h)和强大暴雨(200 mm/h)。

图2 对比工况组合示意图
Fig.2 Sketch map of contrast working conditions
三塔连体结构主体建筑外轮廓约为140 m×180 m×300 m(长×宽×高),为同时保证超高层建筑处于降雨区域和尾流充分发展,计算域尺寸取顺风向X=4 000 m,横风向Y=2 000 m,高度方向Z=600 m。为了兼顾计算效率与精度,同时考虑到建筑物外形复杂,网格划分采用混合网格离散形

图3 计算域及加密网格划分示意图
Fig.3 Computational domain and cipher gridding
数值模拟采用3D单精度分离式求解器,计算参数及边界条件设定如

图4 速度及湍流度剖面示意图
Fig.4 Velocity and turbulence profile
脉动风场求解稳定后加入离散相进行风雨场耦合迭代运算。采用直径为1.0~6.0 mm的雨滴来模拟连续直径分布的降雨(见
风雨耦合迭代完成后,可输出建筑物表面捕捉到的雨滴信息,据此计算雨滴对建筑结构的冲击作用力,计算雨滴冲击荷载时,假定0.01 s时
对该超高层连体建筑进行测压风洞试验,试验详细介绍见文献[

图5 超高层典型截面风场数值模拟结果与风洞试验曲线对比示意图
Fig.5 Comparison between numerical simulation results and wind tunnel test at typical interfaces of the three-tower connected tall buildings

图6 超高层建筑80 m和250 m截面处速度流线图
Fig.6 Speed streamline at 80 m and 250 m sections of the three‑tower connected tall buildings
(1)两个子塔对主塔的流场分布形成了显著的干扰效应。由于子塔的存在,主塔与子塔的“峡谷效应”使得三塔之间的流体出现明显的加速现象,同时子塔的干扰效应也改变了矩形流场的对称分布。
(2)在250 m高度截面处,结构的尾流漩涡区域出现明显的减少现象,但建筑矩形断面的流场仍未出现对称分布,说明子塔对主塔的干扰作用在子塔高度以上(Z≥157.8 m)仍然存在,超高层三塔连体建筑在静力干扰方面呈现明显的三维效应。

图7 雨滴数量与水平末速度分布曲线
Fig.7 Distribution curves of raindrop number and horizontal velocity
表4—6对比给出不同组合工况下建筑物迎风面不同高度范围雨荷载、风荷载及二者比值。由表4—6可知:
(1)迎风面各高度范围内均有雨滴附着,不同组合工况下建筑外表面雨荷载随着高度的增加而增大;不同高度范围雨荷载与风荷载的比值大多分布在3%~15%范围,但在工况3的223.6~261.3 m高度范围,雨荷载达到风荷载的23.81%。
(2)在同一重现期风速下,雨荷载随着雨强的增大而增大,固定雨强下风速的提高将导致建筑物迎风面的雨荷载增大,但雨荷载与风荷载的比值减小。

图8 迎风面典型测点雨压系数分布曲线
Fig.8 Distribution curves of typical rain pressure coefficients at the windward surface
为更清晰展示建筑物各侧面和高度雨滴位置、数量和对应的压力系数,


图9 各工况下超高层建筑外表面雨滴和雨压系数三维分布示意图
Fig.9 Three‑dimensional distribution of raindrop and rain pressure coefficient of three‑tower connected tall buildings
(1)各组合工况下雨滴撞击位置主要集中在建筑物迎风面区域,来流在迎风面的分流和绕流携带部分雨滴远离迎风面中心点。因此,迎风面两侧区域雨压系数较大,中间区域雨压系数较小,受尾流区域漩涡的驱动作用,背风面区域有少量雨滴附着。
(2)各工况下建筑外表面收集到的雨滴总数随着雨强的增大逐渐变多,而随风速的增加迅速减少,其中工况3下雨滴数量最多。
(3)各工况下雨压系数主要集中于迎风面,其余范围数值基本为0,雨压系数最大值为0.243,发生在工况3的0.73H~0.86H高度范围内。
基于最不利组合工况3和超高层建筑各表面雨滴分布特性,
本文系统探讨了风‑雨耦合环境下超高层三塔连体建筑气动力与雨压分布特性,主要涉及风‑雨双向耦合数值模拟、风荷载分布、干扰效应、流场作用机理、雨荷载分布、风‑雨等效压力系数等。得到如下主要结论:
(1)超高层三塔连体建筑的静力干扰呈现明显的三维效应,子塔的干扰效应减小了主塔迎风面和侧风面的风荷载。
(2)主塔表面雨滴捕捉数量随风速增大而减少,随雨强增大而增加,工况3下雨滴捕捉数量最多,并以5 mm直径雨滴占比最大,雨滴水平速度主要分布在10~30 m/s范围内。
(3)不同工况下雨滴撞击位置主要集中在迎风区域,且迎风面雨荷载随着高度的增加逐渐增大,迎风面雨荷载与风荷载最大比值为23.81%,局部测点最大雨压系数达到0.301,发生在工况3的塔顶范围内。
(4)风雨耦合的最不利组合工况为100 a重现期风速(35.77 m/s)和强大暴雨组合(200 mm/h),主塔迎风面雨压系数整体取值建议为0.12,侧风区和背风区整体雨压系数数量级在1
参考文献
余先锋, 谢壮宁, 顾明. 群体高层建筑风致干扰效应研究进展[J]. 建筑结构学报, 2015, 36(3): 1‑11. [百度学术]
YU Xianfeng, XIE Zhuangning, GU Ming. Research progress of wind‑induced interference effects on tall buildings[J]. Journal of Building Structures, 2015, 36(3): 1‑11. [百度学术]
Huang M F, Chan C M, Lou W J. Optimal performance‑based design of wind sensitive tall buildings considering uncertainties[J]. Computers & Structures, 2012, 98: 7‑16. [百度学术]
Blocken B, Carmeliet J. A review of wind‑driven rain research in building science[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2004, 92(13): 1079‑1130. [百度学术]
Kasperski M. Extreme wind load distributions for linear and nonlinear design[J]. Engineering Structures, 1992, 14(1): 27‑34. [百度学术]
Khanduri A C, Stathopoulos T, Bédard C. Wind‑induced interference effects on buildings—A review of the state‑of‑the‑art[J]. Engineering Structures, 1998, 20(7): 617‑630. [百度学术]
Lam K M, Leung M Y H, Zhao J G. Interference effects on wind loading of a row of closely spaced tall buildings[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2008, 96(5): 562‑583. [百度学术]
Yu X F, Xie Z N, Wang X, et al. Interference effects between two high‑rise buildings on wind‑induced torsion[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2016, 159: 123‑133. [百度学术]
Huang D, Zhu X, He S, et al. Characteristics of the aerodynamic interference between two high‑rise buildings of different height and identical square cross‑section[J]. Wind & Structures An International Journal, 2017, 24(5): 501‑528. [百度学术]
Zhang J W, Li Q S. Field measurements of wind pressures on a 600 m high skyscraper during a landfall typhoon and comparison with wind tunnel test[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2018, 175: 391‑407. [百度学术]
Xie J, Irwin P A. Wind‑induced response of a twin‑tower structure[J]. Wind & Structures An International Journal, 2001, 4(6): 495‑504. [百度学术]
Shen G, Sun B, Lou W, et al. Distribution of wind load on symmetrical double tower tall building[J]. Journal of Building Structures, 2003, 24(1): 64‑66. [百度学术]
Song J, Tse K T. Dynamic characteristics of wind‑excited linked twin buildings based on a 3‑dimensional analytical model[J]. Engineering Structures, 2014, 79(79): 169‑181. [百度学术]
Choi E C C. Parameters affecting the intensity of wind‑driven rain on the front face of a building[J]. Journal of Wind Engineering & Industrial Aerodynamics, 1994, 53(1/2): 1‑17. [百度学术]
Choi E C C. Numerical modelling of gust effect on wind‑driven rain[J]. Journal of Wind Engineering & Industrial Aerodynamics, 1997, 72(1): 107‑116. [百度学术]
Cosentino N, Flamand O, Ceccoli C. Rain‑wind induced vibration of inclined stay cables—Part I: Experimental investigation and physical explanation[J]. Wind & Structures An International Journal, 2003, 6(6): 471‑484. [百度学术]
Jing H, He X, Li H. Initial condition and damping effects on rain‑wind induced cable vibration[J]. Journal of Wind Engineering & Industrial Aerodynamics, 2018, 175: 376‑383. [百度学术]
王辉, 李新俊, 潘竹. 建筑(群)立面风驱雨压荷载的数值模拟研究[J]. 土木工程学报, 2014, 47(9): 94‑100. [百度学术]
WANG Hui, LI Xinjun, PAN Zhu. Numerical simulation study of WDR pressure load on facades of buildings[J]. China Civil Engineering Journal, 2014, 47(9): 94‑100. [百度学术]
Arch G, Hajra B, Moravej M, et al. An experimental study to assess the effect of soffit louvered vents on wind loads and wind driven rain intrusion on low rise buildings[J]. Sustainable Cities & Society, 2017, 34: 43‑55. [百度学术]
杨俊涛, 楼文娟. 风驱雨CFD模拟及平均雨荷载计算方法研究[J]. 空气动力学学报, 2011, 29(5): 600‑606. [百度学术]
Yang Juntao, Lou Wenjuan. Research on wind‑driven rain CFD simulation and method calculating mean rain load[J]. Acta Aerodynamica Sinica, 2011, 29(5): 600‑606. [百度学术]
Mcfarquhar G M, List R. The raindrop mean free path and collision rate dependence on rainrate for three‑peak equilibrium and Marshall‑Palmer distributions[J]. Journal of the Atmospheric Sciences, 2010, 48(48): 1999‑2004. [百度学术]
Gunn R, Kinzer G D. The terminal fall velocity for water droplets in stagnant air[J]. Journal of the Atmospheric Sciences, 1949, 6(4): 243‑248. [百度学术]
Rigby E C, Marshall J S, Hitschfeld W. The development of the size distribution of raindrops during their fall[J]. Journal of the Atmospheric Sciences, 2010, 11(5): 362‑372. [百度学术]
刘顺, 黄生洪, 李秋胜,等. 基于欧拉多相流模型的桥梁主梁三维风驱雨数值研究[J]. 工程力学, 2017, 34(4): 63‑71. [百度学术]
Liu Shun, Huang Shenghong, LI Qiusheng, et al. 3D numerical simulation of wind‑driven rain on bridge deck sections based on Eulerian multiphase model[J]. Engineering Mechanics, 2017, 34(4): 63‑71. [百度学术]
Douvi E, Margaris D. Aerodynamic performance investigation under the influence of heavy rain of a NACA 0012 airfoil for wind turbine applications[J]. International Review of Mechanical Engineering, 2012, 6(6):1-8. [百度学术]
Elshaer A, Bitsuamlak G, Damatty A E. Enhancing wind performance of tall buildings using corner aerodynamic optimization[J]. Engineering Structures, 2017, 136: 133‑148. [百度学术]
陈博文. 低矮房屋表面风雨压力CFD数值模拟[D]. 哈尔滨:哈尔滨工业大学, 2009. [百度学术]
CHEN Bowen. Low building surface wind pressure CFD numerical simulation[D]. Harbin: Harbin Industrial University, 2009. [百度学术]
Ke S, Wang H, Ge Y. Wind load effects and equivalent static wind loads of three‑tower connected tall buildings based on wind tunnel tests[J]. Structural Engineering & Mechanics, 2016, 58(6): 967‑988. [百度学术]