LES Research on Wind Effect on Tapered Super-tall Buildings
-
摘要: 风荷载往往是超高层建筑结构安全性和舒适性的控制性荷载。本文通过大涡模拟技术并结合作者建议的一种新的湍流入口生成方法(NSRFG方法),对四种不同锥度的楔形建筑模型进行了风效应的数值模拟对比研究,以比较不同锥度对结构风荷载和风振响应的影响,及检验本文数值模拟方法的适用性。数值模拟结果表明,虽然不同锥度的建筑模型平均风压和脉动风压分布规律相似,但在结构响应方面,结构横风向峰值基底弯矩响应随模型锥度增大显著减少,表明采用适当锥度的体形可以有效地减小结构横风向气动荷载。本文数值模拟结果和风洞试验结果整体规律相符,表明结合改进的入流湍流生成技术的大涡模拟具有一定的精度,可以给超高层建筑气动外形优化研究提供有价值的参考。Abstract: Wind load was often the key load for the safety and comfort of super-tall building structures.Large eddy simulation(LES)technique combined with a new inflow turbulence generator(NSRFG)was employed to study the wind effects on tapered building models with different tapers, which aimed to compare the effect of different taper on wind loads and wind-induced responses of the structure and validate the applicability of the numerical method of this paper.Numerical simulation showed thatthe distributions of the average and fluctuating pressure on building surfaces were similar, whilefor the wind-induced response, the peak base bending moment(BBM)of different models in across wind direction decreased remarkably with the increasing of the taper.Therefore, a building with bottom-up gradually shrinking and an appropriated tapershape can effectively decrease the aerodynamic wind load in across-wind direction.In addition, results from the numerical simulation were coincided with that from the wind tunnel test in the whole, which indicated that LES technique combined with an appropriated inflow turbulence generator would havecertain accuracy and canprovide valuable reference to the aerodynamic shape optimization of super-tall buildings.
-
1. 前言
随着社会经济的发展,国内近年来已出现了越来越多的超高层建筑。由于高层建筑结构具有柔度大、质量轻、自振频率以及阻尼较低的特点,结构的风效应逐步成为控制高层建筑安全性、舒适性和经济性的最重要的因素之一。超高层建筑结构属风敏感性结构,对风荷载所造成的扰动十分敏感,在强风作用下容易受到损坏。强风作用下超高层建筑结构的风荷载是结构安全性和舒适性的控制荷载,因此对风荷载的研究和分析已成为超高层建筑结构设计和体型优化的重要环节。
建筑外形是决定作用于结构上气动荷载的关键因素,通过合理改变建筑物的外形,对建筑外形进行优化设计,可以显著改变作用于建筑结构的风荷载和风致响应。通过传统的增强结构的强度与刚度已经很难满足现今结构设计者对于抗风设计的要求,于是合理地优化作用于超高层建筑结构上的风荷载已成为结构风工程的一个重要研究方向,如台北的101大厦及阿联酋迪拜大厦的外形都进行了抗风优化研究。
随着计算风工程在国内外的发展,数值模拟方法已经成为结构风工程的主要研究方法之一,大涡模拟(Large Eddy Simulation,LES)在建筑风效应的模拟中得到越来越广泛的应用。随着数值模拟技术和计算机性能的提升,近年来研究人员尝试将大涡模拟运用到工程实践中。例如,李秋胜等人[1-2]采用LES分析了高层建筑的空气动力学效应,Dagnew,A.等人[3]和Daniels,S.J.等人[4]运用LES对高层建筑的风荷载进行了评估。
目前,采用LES数值模拟方法对超高层建筑的气动外形进行参数优化研究还较少,分析可能原因,由于超高层建筑风荷载效应对气动外形变化较敏感,这需要数值模拟具有很高的精度和分辨率,而通常数值风洞技术还很难做到和物理风洞类比。
本文旨在运用计算流体动力学(Computational fluid dynamics,CFD)数值模拟技术,结合作者建议的一类新的入流湍流生成技术,对不同楔形外形的超高层建筑模型进行大涡模拟研究,以评估不同锥度楔形外形建筑模型的气动荷载和风致响应影响规律,并探讨采用新的入流湍流生成方法的LES数值模拟技术研究超高层建筑气动外形优化的可行性。
2. 数值模拟概况
2.1 算例模型与结构动力特性
本文选取4个不同外形的建筑模型进行计算,模型高度均为910mm,模型底边宽度为100mm,缩尺比取为1:600,相当于实际建筑高度546m,建筑高宽比9.1。4个模型中包含1个正方形等截面模型(模型1),3个锥度(模型底部截面边长与顶部截面边长差值与模型高度的比值)分别为2.2%,4.4%,6.6%的楔形建筑模型(模型2,模型3,模型4) 各模型形状和尺寸如图 1所示。为方便叙述,锥度用ξ表示并定义:
$\xi = \frac{{{D_1} - {D_2}}}{H}$
(1) 式中,D1和D2分别为模型底部和顶部截面边长;H为模型高度;
建筑模型动力计算所采用的前两阶自振周期均为10.2s(自振频率0.098 Hz),假定振型沿高度方向线性分布,阻尼比取为3.5%,基本风压参照深圳地区,选取重现期100年、50年和10年的基本风压分别为0.90kPa、0.75kPa和0.45kPa。
2.2 湍流入口的生成
入口湍流的模拟一直以来是大涡模拟的关键问题之一,其意义等同于物理风洞试验中需首先准确模拟边界层风场特性,因此有许多专门研究针对入口湍流的准确模拟、以及入口湍流影响等进行深入探讨。本文采用作者建议的一种湍流入口生成方法——NSRFG方法[5]生成大涡模拟的湍流入口,其表达式如下式所示:
$u\left( {x,t} \right) = \sum\limits_{n = 1}^N {{p_n}} \sin ({{\rm{k}}_n} \cdot \tilde x + 2\pi {f_n}t + {\varphi _n})$
(2) 其中,u={u,v,w},
$\tilde x = \frac{x}{{{L_n}}}$ ,x={x,y,z};${{\rm{p}}_n} = \sqrt {2{\rm{S}}({f_n})\Delta f} $ ;S={su,sv,sw}为三个方向上的风速谱(卡曼谱),N为功率谱离散数目,${f_n} = \frac{{2n - 1}}{{2\Delta f}}$ ,Δf为带宽,${\varphi _n} - U\left( { - \pi ,\pi } \right)$ ,另外参数Ln取值由下式得到:${L_n} = \frac{{{U_{av}}}}{{\gamma C{f_n}}}$
(3) 其中,
$\gamma = \left\{ {\begin{array}{*{20}{l}} {3.7{\beta ^{ - 0.3}},\beta < 6.0}\\ {2.1,\beta \ge 6.0} \end{array}} \right.$ ,β=CD/Lu,C为衰减系数;Lu为u方向湍流积分尺度;D=0.5~1.0h,h为建筑高度;Uav为平均风速。通过编程将NSRFG方法写成UDF用户自定义程序并与商用软件包ANSYS FLUENT链接进行计算。
本文研究对象原型为500m以上的超高层建筑,由于此高度已经超过了我国现行荷载规范的梯度风高度,故没有采用《建筑结构荷载规范(GB50009-2001)》的风场参数取值,而是依据英国工程科学数据库(ESDU)的方法定义了数值模拟入口的平均风剖面和湍流度剖面,如图 2所示。
2.3 网格划分与参数设定
数值模拟在华南理工大学风洞实验室的数值风洞平台商用CFD软件包ANSYSFLUENT上进行。参照AIJ[6]和相关文献的规定,确定流域尺寸如图 3a所示,堵塞率为0.8%,满足小于3%的要求。边界条件方面,模型表面和地面设置为无滑移壁面,流域左右壁面和顶面设置为对称面,流域入口为速度入口,出口采用压力出口,如图 3a所示。由于建筑模型较为规整,本次数值模拟采用结构网格进行计算。考虑到计算精度和计算资源的平衡,最小网格尺寸取为0.06D(D为模型底面宽度),网格拉伸比例为1.05。近墙壁y+值约为30~75,适用于壁面函数。不同模型网格划分方式一致,网格数量约为320万,网格划分情况如图 3b所示(仅给出模型1的情况)。
本次数值模拟计算采用压力隐式分割算法(Pressure Implicit with Splitting of Operators,PISO),空间离散采用中心差分格式,时间离散采用二阶隐式格式。亚格子模型采用壁面自适应局部涡粘模型(Wall-Adapting Local Eddy-viscosity model,WALE)。计算在华南理工大学风洞实验室的64核AMAX工作站进行,采用16核并行计算。时间步长选为0.002s,数值模拟计算总时间10s,并取后8s的数据作为分析数据。参考高度采用模型高度0.91m,在参考高度处的风速为11m/s。参考高度处建筑模型雷诺数约为7.8×104。
3. 结果分析
3.1 表面风压分布和流场分布
经计算,得到四个模型的表面平均风压系数和脉动风压系数分布,如图 5所示。由图 5可得,拥有不同锥度的模型平均风压系数和脉动风压系数大小差别不大,分布相似。随着锥度的增大,风压分布图形趋于拉长。值得注意的是,在模型背风面脉动风压系数分布方面,随着锥度的增大,图形辐射中心由两个化为一个,且有上移的趋势。具体风压系数分布情况见图 5。
图 6给出模型1在10s时刻湍流场示意图。由图 6可见,来流风在建筑迎风两侧发生分离,并产生分流涡旋,在来流风的作用下交替向下游泄出,产生规则的漩涡脱落现象。模型的尾流流场出现有别于来流湍流的复杂特征湍流,呈现出典型的钝体分离流动特性。
3.2 横风向基底气动弯矩和结构响应
强风作用下风荷载将成为影响结构安全的控制荷载,超高层建筑的横风向风致振动响应通常会超过顺风向响应而成为结构设计中要考虑的重点[7]。
当结构的振型沿高度线性分布时,基底弯矩和结构的1阶广义力成正比,因此基底弯矩是衡量结构整体气动性能的重要指标。为研究不同锥度ξ对模型基底弯矩的影响,在数值模拟过程中记录每一个时间步的横风向基底弯矩系数CMx,其定义如式(4) 所示:
${C_{{\rm{Mx}}}} = \frac{{{M_{\rm{x}}}}}{{\frac{1}{2}\rho {V_h}^2D{H^2}}}$
(4) 其中,Mx为建筑模型的横风向基底弯矩,Vh为参考高度风速,D为建筑模型的宽度(此处取底部宽度),H为建筑模型高度。
图 6给出了不同锥度ξ的模型横风向基底弯矩系数功率谱密度的对比。由图 6可见,气动弯矩随模型截面的收缩而衰减,同时基底弯矩的主峰频率随锥度ξ的增大而增大。此结果与文献[7]的风洞试验结果基本一致,进一步验证了本文数值模拟方法的适用性和准确性。
随着ξ的增大,基底弯矩的主峰频率也增大,使得建筑的基底气动弯矩的主能量分布向高频端移动,这可能会对结构带来不利的影响。图 6中竖线对应的频率值分别为100年、10年重现期下原型结构1阶自振频率换算到模型的值,其换算关系为:
${f_m} = \frac{{{f_p}}}{{{\lambda _f}}} = \frac{{{f_p}{\lambda _L}}}{{{\lambda _V}}}$
(5) 其中,fm和fp为模型与原型频率;λf、λL、λV分别为原型和模型的频率比、缩尺比和风速比。
根据结构动力学随机振动理论,由式(4) 记录的基底弯矩时程和结构动力特性可计算得到考虑结构动力放大效应的横风向基底弯矩功率谱密度和标准差,计算公式如式(6)~(8) 所示:
${S_{{M_D}}}\left( f \right) = {\left| {H\left( f \right)} \right|^2}{S_M}\left( f \right)$
(6) ${\sigma _{{M_D}}} = \int_0^\infty {{S_{{M_D}}}\left( f \right)} df \approx \sqrt {\sigma _{^M}^2 + \pi {f_0}{S_M}({f_0})/\left( {4\zeta } \right)} $
(7) $H{\left( f \right)^2} = \frac{1}{{\{ {{[1 - {{(f - {f_0})}^2}]}^2} + 4\zeta _ \cdot ^2{{\left( {f/{f_0}} \right)}^2}\} }}$
(8) 式中,|H(f)|2为结构机械导纳函数;f0和ζ分别为结构1阶自振频率和阻尼比。式(7) 中的σM2和πf0SM(f0)/(4ζ)分别为弯矩响应的背景分量和共振分量。在得到基底弯矩标准差之后,模型的峰值基底弯矩响应由下式计算得到:
${M_{\max }} = \bar M + g \cdot {\sigma _{{M_D}}}$
(9) 式中,Mmax为峰值基底弯矩响应;M为基底弯矩响应平均值(横风向取为0);g为峰值因子,由下式可得:
${g_x} = \sqrt {2\ln ({\upsilon _x}T)} + 0.577/\sqrt {2\ln ({\upsilon _x}T)} $
(10) 其中,T为统计时间,通常取600s;υx是单位时间内平均位移的水平交叉率数,由下式求得:
${\upsilon _x} = {[\int\limits_0^\infty {{{\rm{f}}^2}} {S_{x1}}\left( f \right)dn/\int\limits_0^\infty {{S_{x1}}\left( f \right)} dn]^{1/2}}$
(11) 关于高层建筑结构峰值响应计算的具体方法详见文献[8]。本文计算响应用的结构自振周期是根据已有工程和规范建议估计的,但实际的建筑结构周期会在这个基础上发生一定范围内变化。为全面考虑建筑外形对结构气动性能的影响必须分析该结构风致响应对周期变化的敏感性。图 7给出了不同ξ值下模型峰值基底弯矩随周期比的变化曲线图。此处周期比的定义是实际建筑结构周期与本文假定的计算周期之比,周期比的变化范围为0.6~1.4,相当于结构的计算周期范围为6.12~14.28s,这一范围基本覆盖了《建筑结构荷载规范(GB50009-2001)》计算结构周期的经验公式范围,基本满足此类高度超高层建筑自振周期的要求。
从图 7可得,当周期比大于等于0.85时,峰值基底弯矩随着锥度ξ的值增大而减小,由此看来建筑模型的截面收缩可十分有效地降低结构峰值基底弯矩。当周期比小于0.7时,ξ=2.2%和ξ=4.4%两种外形模型的峰值基底弯矩反而出现超过等截面模型的情况,可见变化规律与图 6所示基底弯矩功率谱密度的变化规律一致,且在所关心的周期比变化范围内,ξ=6.6%模型的峰值基底弯矩响应最小。此计算结果与文献[8]风洞试验结果变化规律整体相符,证明结合改进的入流湍流生成技术的大涡模拟具有相当的精度,可以给超高层建筑气动外形的优化研究提供有价值的参考,对不同外形的建筑模型风效应进行定性的比较和判断。
但不可否认的是,数值模拟结果与试验对比仍存在一定差距,可能与大涡模拟的滤波效应、数值模拟参数设置等因素有关,将在后续的工作中进行详细比较分析,以进一步提高数值模拟结果的精度。
4. 结论
本文采用CFD大涡模拟技术并采用作者建议的一种湍流生成方法——NSRFG方法对拥有不同锥度的楔形外形超高层建筑模型进行了风效应的研究,得出以下结论:
(1) 不同模型的平均风压系数和脉动风压系数大小差别不大,分布规律相似;随着锥度的增大,风压分布图形趋于拉长。
(2) 结构横风向峰值基底弯矩响应随模型锥度增大显著减少,这说明超高层建筑采用由下到上逐渐收缩的外形可以有效地减小结构横风向气动荷载。一定周期比范围内,锥度越大,横风向气动弯矩的削减效果越好。锥度分别为2.2%,4.4%,6.6%的建筑模型在周期比为1的情况下,相对于等截面模型横风向基底弯矩响应减少了58.8%,84.5%和86.2%。同时锥度越大,建筑的漩涡脱落频率也会相应增大,这有可能造成在结构自振频率高或者低风速情况下,锥度高的建筑风致响应比等截面还大的情况。此类建筑的外形优化还需和角区优化、避难层开透风槽等局部气动优化措施综合考虑。
(3) 采用作者建议的新的LES入流湍流生成方法NSRFG方法,数值模拟结果与风洞试验客观上虽存在一定差距(可能与大涡模拟滤波效应和数值模拟参数设置有关),但仍具有相当的精度和分辨率,可以较为准确地预测模型漩涡脱落频率,给超高层建筑气动外形的优化研究提供有价值的参考,可以达到在初步设计阶段对不同外形的建筑模型风效应进行体型优选的目的。
-
[1] 李秋胜, 刘顺.基于大涡模拟的平屋盖锥形涡数值模拟研究.湖南大学学报(自然科学版), 2015, 42(11): 72-79. DOI: 10.3969/j.issn.1674-2974.2015.11.010 [2] 卢春玲, 李秋胜, 黄生洪, 等.深圳平安国际金融大厦风致响应大涡模拟分析.建筑结构学报, 2012, 33(11): 1-11. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJB201211002.htm [3] Dagnew, A., Bitsuamlak, G.T.Large eddy simulation for wind-induced response of tall buildings located in a city center.In:Proceedings of the 2012 Engineering Mechanics Institute and 11th ASCE Joint Specialty Conference on Probabilistic Mechanics and Structural Reliability.Notre Dame, 2012.
[4] Daniels, S.J., Castro, I.P., Xie, Z.T.Peak loading and surface and surface pressure fluctuations of a tall model building.Journal of Wind Engineering and Industrial Aerodynamics, 2013, 120: 19-28. DOI: 10.1016/j.jweia.2013.06.014
[5] Yuanlin Yu, Yi Yang, Zhuangning Xie.An improved inflow turbulence generator for LES evaluation of wind load on high-rise buildings[C].Processing of the 8th international colloquium of buff bodies aerodynamics and application, Boston, 2016.
[6] Yoshihide Tominaga, Akashi Mochida, Ryuichiro Yoshie, Hiroto Kataoka, Tsuyoshi Nozu, Masaru Yoshikawa, Taichi Shirasawa.AIJ guidelines for practical applications of CFD to pedestrian wind environment around buildings.Journal of Wind Engineering and Industrial Aerodynamics, 2008.1749-1761. http://www.sciencedirect.com/science/article/pii/S0167610508000445
[7] 谢壮宁, 李佳.强风作用下楔形外形超高层建筑横风效应试验研究.建筑结构学报, 2011, 32(12): 118-126. http://www.cnki.com.cn/Article/CJFDTOTAL-JZJB201112014.htm [8] 谢壮宁, 方小丹, 倪振华.超高层建筑的等效静风荷载:扩展的荷载响应相关方法[J].振动工程学报, 2008, 21(4): 398-403. http://www.cnki.com.cn/Article/CJFDTOTAL-ZDGC200804015.htm