2. 风工程四川省重点实验室 (西南交通大学),成都 610031
2. Wind Engineering Key Laboratory of Sichuan Province (Southwest Jiaotong University), Chengdu 610031, China
在高层建筑顺风向等效静力风荷载及风振响应计算时,为描述脉动风荷载空间分布特性,需要明确脉动风荷载空间相关性和气动导纳[1].就脉动风荷载相关性而言,国内外相关规范沿用了Davenport[2]提出的“脉动风压与脉动风速具有一致的空间相关性”这一假定.已有研究表明[3-4]:脉动风压的相关性要显著高于脉动风.气动导纳函数对脉动风荷载具有空间相关性折减作用[1],对于采用测压方法直接获得建筑结构表面脉动风压的研究,可不用考虑.
Davenport[2, 5]最早提出脉动风速相干函数可采用频率及两点间距的简单指数函数表示.ESDU[4]及Krenk[6]指出该模型的衰减因子存在较大不确定性,但由于缺乏高层建筑脉动风荷载空间相关性数据,该不确定性已成为结构响应计算不可忽略的误差之一.中国现行《建筑结构荷载规范》[7]采用了Shiotani[8]提出的与两点间距有关的相干函数模型,这对简化计算有利,但该函数与频率无关,在整个结构频率区段内取为常数可能导致计算误差.顾明等[9]针对不同截面形式的高层建筑进行了大量研究,给出一种便于工程应用的脉动风荷载相干函数经验公式.黄冬梅等[10]通过某一特殊外形超高层建筑测压试验,分析了超高层建筑脉动风荷载的空间相关性,并得到了相应的相干函数表达式.
在以上研究中,一般未考虑紊流风场空间结构与脉动风荷载相关性之间的关系.研究表明,当紊流积分尺度大于或者接近结构尺寸时,脉动风压相关性会得到增强[11].此外,在钝体 (流线型箱梁主梁断面、矩形断面) 相干函数研究中,认为紊流度对脉动风荷载空间相关性影响较小[12-13],如文献[11, 14-15]提出的相干函数模型中,均主要考虑紊流积分尺度影响,该结论在高层建筑中是否一致仍需讨论.
为进一步探讨高层建筑顺风向脉动风荷载及脉动风速的空间相关性,减小因相干函数不确定性导致的误差,本文在三类大尺度紊流风场中,通过刚性模型同步测压试验,对矩形模型顺风向脉动风荷载空间结构进行研究.以典型钝体结构——矩形断面作为研究对象,考虑到矩形高层建筑宽厚比较少有超过3的情况,试验选取模型宽厚比为2:1.根据试验结果,对脉动风速、顺风向脉动风荷载、迎风面和背风面脉动风压的空间相关性进行对比分析.实验结果验证了高层建筑脉动风荷载相关性强于脉动风速这一结论.在传统相干函数模型基础上,综合考虑紊流风场、间距及结构特征尺寸对脉动风荷载相干函数的影响,提出可计及流场参数,竖向间距和结构特征尺寸的矩形建筑顺风向脉动风荷载相干函数经验公式,可为矩形高层建筑顺风向脉动风荷载精细化分析提供参考.
1 试验概况 1.1 测压试验简介本试验在西南交通大学XNJD-3工业风洞进行.试验段截面为22.5 m (宽)×4.5 m (高),风速范围为1.0~16.5 m/s.模型高2.3 m,矩形截面长20 cm,宽10 cm,见图 1.为增加高度方向 (竖向) 测点层组合个数,该模型高宽比较大,考虑到本文主要从基础研究的角度出发,分析矩形高层建筑脉动风荷载空间相关性,模型高宽比可适当放宽.在高度方向上共布置11层测点,离地高度分别为19、24、44、67、90、110、125、130、133、151、181 cm.每层测点均为52个,测点布置见图 2.试验采用电子扫描阀 (型号DSM-3400) 测量风向角为0°和90°时模型的表面风压,并定义来流垂直长边时风向角为0°(图 2).
紊流风场由尖塔、挡板和粗糙元模拟产生,得到建筑结构荷载规范[7]中的B、C和D三类风场,见图 3.模拟的平均风速剖面与规范[7]要求的B、C和D类地表一致,紊流度剖面比规范值稍大,但符合试验要求.图 3中,zg为梯度风高度,z为竖向高度,Ug为梯度高度风速,U/Ug为平均风速与梯度风速比值.
本文采用两个Cobra Probe风速探头测量空间不同位置的脉动风速时程,测量间距为0.03~1.62 m.图 4以10#断面高度为例,对C类风场顺风向脉动风速谱进行拟合 (f为风速频率).可见模拟风场的无量纲风速谱与von-Karman谱拟合较好,并通过拟合Karman谱获取顺风向紊流积分尺度.试验风场对应的实际紊流积分尺度在110~160 m范围,与实际情况基本相符.
频域内,脉动风的竖向相关性可由无量纲标准交叉谱来表达,一般可近似取标准交叉谱的实部简化表示:
$ {\rm{coh}}\left( {\Delta z,f} \right) = \frac{{\left| {{S_{{\rm{uu}}}}\left( {{z_1},{z_2},f} \right)} \right|}}{{\sqrt {{S_{\rm{u}}}\left( {{z_1},f} \right){S_{\rm{u}}}\left( {{z_2},f} \right)} }}, $ | (1) |
式中:Suu为互谱实部,Su为各位置点谱,f为风速频率,Δz=z2-z1.
$ {\rm{coh}}\left( {\Delta z,f} \right) = \exp \left( { - c\Delta zf/U} \right), $ | (2) |
式中:coh (Δz, f) 为脉动风速相干函数,c为指数衰减系数,一般竖向建议c=10[5],文中将作为待定系数拟合.
式 (2) 结构简单,应用广泛,但存在一定的缺陷:在低频区,即使间距很大,相干函数仍可接近1,这与大量实测结果不符.Krenk[6]通过引入修正频率,提出了修正指数公式来解决上述问题,该修正是基于标准化的Karman谱模型,
$ \begin{array}{*{20}{c}} {{\rm{coh}}\left( {{k_1},\gamma } \right) = }\\ {\frac{2}{{\Gamma \left( \gamma \right)}}\left[ {{{\left( {\frac{{{k_1}\gamma }}{2}} \right)}^\gamma }{{\rm{K}}_\gamma }\left( {{\kappa _1}\gamma } \right) - {{\left( {\frac{{{k_1}\gamma }}{2}} \right)}^{\gamma + 1}}{{\rm{K}}_{1 - \gamma }}\left( {{k_1}\gamma } \right)} \right],} \end{array} $ | (3) |
式中:Γ为伽玛函数,Kγ和K1-γ为第二类修正Bessel函数,k1=2πf/U为波数.当γ=1/2时,上式可简化为
$ {\rm{coh}}\left( {\Delta z,f} \right) = \left( {1 - \frac{{c\Delta z{f_x}}}{{2U}}} \right)\exp \left( { - \frac{{c\Delta z{f_x}}}{U}} \right), $ | (4) |
其中fx为修正频率
$ {f_x} = f \cdot \sqrt {1 + {{\left( {1/{k_1}L} \right)}^2}} . $ | (5) |
长度模量L=1.34Lu,Lu为脉动风顺风向积分尺度.高频和空间距离较小的标准交叉谱可用上述公式得到准确表示.但在低频和空间距离较大时,还需要发展更为精确的湍流理论,如快速畸变理论等.
采用式 (1) 对试验数据拟合,可得不同间距下的相干函数.采用非线性最小二乘法,根据式 (2)~(4),对每类风场55个间距组合的风速相干函数进行拟合.结果表明:Davenport模型的衰减系数可取为10.4,Krenk模型衰减系数取为6.同时图 5还给出了风速相干函数在不同风场中的分布趋势.结果表明,三类风场中的脉动风速相干函数分布规律基本一致,可认为受风场类型影响较小.
在频域内,一般采用相干函数[9, 10, 16-17]描述空间任意两点脉动风荷载的相关度,而在时域内则采用相关系数来描述.在0°风向角下,矩形高层建筑的顺风向脉动风荷载相关系数随间距变化如图 6所示.图 6结果再次证实高层建筑脉动风荷载比脉动风速具有更强的相关性,与Dyrbye[1]和Kareem[3]结论一致.同时,B类风场中顺风向脉动风荷载相关性要略强于C、D类风场,其原因可能是B类风场的紊流积分尺度略大于C、D类风场.图 6中,横坐标为约化间距Δz/B,B为矩形断面迎风面宽度.
由于缺乏顺风向脉动风荷载在迎风面和背风面的相关性数据,一般假定迎风面与背风面脉动风压全相关[1],但该假定未得到有效验证.在此基础上,认为该假定可能产生的误差与采用脉动风速相关性代替脉动风荷载相关性产生的误差能在一定程度上相互抵消[3].以上假定为高层建筑风振响应计算带来了不确定性,可能导致较大的计算误差.
图 7给出了B类风场0°风向角下,迎风面和背面相关系数沿竖向变化规律.发现迎风面脉动风压相关系数与矩形断面整体荷载趋势一致,但相关性会略小.随着间距增大,迎风面脉动压力相关系数会逐渐低于断面总体力.该结果表明顺风向脉动风荷载主要受迎风面脉动风压控制,但仅考虑迎风面脉动风压相关性则可能低估顺风向脉动风荷载.对于背风面,当竖向间距较小时,脉动压力相关性并不能完全忽略.当Δz/B>1,脉动风压相关系数快速衰减,明显小于迎风面和断面整体阻力.该结果表明若忽略迎风面和背风面脉动风压的弱相关性,即认为顺风向风荷载相关性仅与迎风面脉动风压有关,可能会引起高层建筑顺风向脉动风荷载取值偏不安全.采用已经反映了脉动风及其气动贡献、实际结构的风荷载相干函数,是获得更精确结构风荷载的有效方法.
通过测压试验,同步获得各测点的脉动风压时程,然后将各截面 (或层) 压力时程积分可得到顺风向脉动风荷载 (脉动阻力).这里讨论的各层间脉动阻力相干特性已经考虑了脉动风速、迎风面和背风面脉动风压的综合影响.
图 8给出了0°和90°风向角下,顺风向风荷载相干函数在高度方向 (竖向) 的变化情况.相干函数具有如下特性:1) 相干函数与频率、间距成反比,相干函数衰减量大小取决于旋涡尺寸;2) 间距较大时,相干函数在低频区取值小于1,并随间距增大不断减小;3) 同一风场中,当间距组合的Δz/U接近时,如:5# vs 6#和8# vs 10#,(Δz/U≈0.03 B类和D类,Δz/U≈0.04 C类),其相干函数具有基本一致的分布规律,Δz/U的值越大,衰减速率越大,表明顺风向脉动风荷载相干函数同时受竖向间距和紊流参数空间变化控制;4) 宽边迎风时相关性强于窄边迎风,宽厚比对顺风向脉动风荷载相关性存在一定影响;5) 风场类别对相干函数有一定影响,表现出高紊流度风场顺风向脉动风荷载相干函数衰减量更大,这可能是由于B类风场的积分尺度大于C、D类风场引起的.
文献[10-11, 18]在传统脉动风速相干函数的基础上进行了一系列改进,应用较广泛的顺风向脉动风荷载相干函数模型为
$ {\rm{coh}}\left( {\Delta z,f} \right) = \delta \exp \left( { - \frac{{c\Delta zf}}{U}} \right), $ | (6) |
式中δ、c为待拟合参数.上式中,考虑了间距和频率对脉动风荷载相干函数的影响,可通过δ对低频值进行修正,δ在0~1内取值.指数衰减系数c通过试验结果拟合,比较符合实际情况.式 (6) 相对于传统脉动风速相干函数模型已有较大改进,但仍需进一步完善:式中未考虑紊流风场空间变化对相干函数的影响,并在低频和间距较大或较小时未考虑频率修正,与实际结果存在一定偏差.
根据顺风向脉动风荷载相干函数分布特性,特别是反映脉动风荷载相关性高于脉动风速相关性这一事实,参考已有文献[6, 10-11, 18-19],提出矩形高层建筑顺风向脉动风荷载相干函数表达式:
$ {\rm{coh}}\left( {\Delta z,s,f} \right) = \left( {1 - \frac{{\eta c\left( s \right)f'\Delta z}}{U}} \right)\exp \left( { - \frac{{c\left( s \right)f'\Delta z}}{U}} \right), $ | (7) |
$ s = \Delta z{B^{\alpha /2}}/{H^{1 + \alpha /2}}, $ | (8) |
$ f' = \sqrt {{f^2} + {{\left( {U/2\pi L} \right)}^2}} . $ | (9) |
式中:无量纲参数s反映了迎风面宽度、风场类别和间距的综合影响,α为风剖面指数,H为模型高度,η反映高频影响,f′为修正频率,c(s) 为指数衰减系数,L=1.34Lu,系数1.34为Karman谱修正系数理论值[11].
采用非线性最小二乘法,图 9给出了B类风场0°风向角下的部分拟合结果.将式 (6) 和Davenport模型所代表的相干函数模型作为对比,同时给出了拟合结果.对比结果表明,本文给出的顺风向脉动风荷载相干函数公式具有更高的拟合精度.
通过拟合,η取值为常数,η=1/3.c可用无量纲系数s的二次项表达式较好拟合,结果见图 10.
$ c\left( s \right) = - 66.26{s^2} + 27.39s + 2.354. $ | (10) |
图 11以B类风场中0°风向角的部分组合为例,给出顺风向脉动风荷载相干函数随无量纲参数s的变化情况.本文结果符合实际相干函数空间分布规律,在低频,随着s增大,顺风向脉动风荷载相关性逐渐降低.此外,在考虑了紊流风场和结构特征尺寸修正后,本文模型在一定程度上克服了传统脉动风速相干函数的不足.
1) 矩形高层建筑顺风向脉动风荷载相关性强于脉动风速的相关性.
2) 矩形模型顺风向脉动风荷载相关性强于迎风面脉动风压的相关性,背风面脉动风压相关性较弱,但在小间距时不能完全忽略.
3) 风场类别对脉动风速和顺风向脉动风荷载相关性有一定影响,相干函数在高紊流度风场中的衰减速率要大于低紊流度风场.
4) 矩形高层建筑顺风向脉动风荷载相关性与频率、空间间距、紊流风场特性密切相关.其与频率和间距成反比,与紊流积分尺度成正比,同类风场中,当相干函数的Δz/U相同时,其相干函数分布基本一致.
5) 本文给出的矩形高层建筑顺风向脉动风荷载相干函数经验公式对高层建筑脉动风荷载的精细化分析具有一定参考价值.
[1] | DYRBYE C, HANSEN S O. Wind loads on structures[M]. New York: Wiley, 1997. |
[2] | DAVENPORT A G. The application of statistical concepts to the loading of structures[J]. ICE Proceedings, 1961, 19(4): 449-472. DOI: 10.1680/iicep.1961.11304 |
[3] | KAREEM A. Synthesis of fluctuating along wind loads on buildings[J]. Journal of Engineering Mechanics, 1986, 112(1): 121-125. DOI: 10.1061/(ASCE)0733-9399(1986)112:1(121) |
[4] | ES DU. Characteristic of atmospheric turbulence near the ground. part Ⅲ: variations in space and time for strong winds: neutral atmosphere[M]. London: Engineering Sciences Date Unit, 1986. |
[5] | DAVENPORT A G. Gust loading factors[J]. Journal of Structural Division, ASCE, 1967, 93: 11-34. |
[6] | KRENK S. Wind field coherence and dynamic wind forces [C]// IUTAM symposium on advances in nonlinear stochastic mechanics. Netherlands: Springer, 1996: 269-278. |
[7] |
建筑结构荷载规范: GB 50009-2012[S]. 北京: 中国建筑工业出版社, 2012.
Load code of the design of building structures: GB 50009-2012[S]. Beijing: China Architecture & Building Press, 2012. |
[8] | SHIOTANI M, AVAI H. Lateral structure of gusts in high winds [C]// International Conference on the Wind Effect on Buildings and Structures. Cambridge:Cambridge University Press, 1967:535-555. |
[9] |
顾明, 张建国. 高层建筑顺风向脉动荷载相干性研究[J].
土木工程学报, 2008, 41(11): 18-22.
GU Ming, ZHANG Jianguo. Coherence analysis of along-wind fluctuating loads on high-rise building[J]. China Civil Engineering Journal, 2008, 41(11): 18-22. DOI: 10.3321/j.issn:1000-131X.2008.11.004 |
[10] |
黄东梅, 朱乐东. 超高层建筑风力空间相干函数研究-分析归纳法[J].
土木工程学报, 2010, 43(9): 32-39.
HUANG Dongmei, ZHU Ledong. Study on spatial correlation functions of wind loads on a super-tall building-analysis and induction method[J]. China Civil Engineering Journal, 2010, 43(9): 32-39. |
[11] | LI S P, LI M S, LIAO H L. The lift on an aerofoil in grid-generated turbulence[J]. Journal of Fluid Mechanics, 2015, 771: 16-35. DOI: 10.1017/jfm.2015.162 |
[12] | HANN F L. The effects of turbulence on the aerodynamics of long-span bridge [D]. South Bend Indiana: University of Notre Dame, 2000. |
[13] | LI Q S, MELBOURNE W H. The effect of large-scale turbulence on pressure fluctuations in separated and reattaching flows[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1999, 83: 159-169. DOI: 10.1016/S0167-6105(99)00069-0 |
[14] | KIMURA K, FUJINO Y, NAKATO H, et al. Characteristics of buffeting forces on flat cylinders[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1997, 69/70/71: 365-374. |
[15] | LAROSE G L, MANN J. Gust loading on streamlined bridge decks[J]. Journal of Fluid and Structures, 1998, 12(5): 511-536. DOI: 10.1006/jfls.1998.0161 |
[16] |
顾明, 叶丰. 典型超高层建筑风荷载频域特性研究[J].
建筑结构学报, 2006, 27(1): 30-36.
GU Ming, YE Feng. Frequency domain characteristics of wind loads on typical super-tall buildings[J]. Journal of Building Structure, 2006, 27(1): 30-36. |
[17] | SIMIU E, SCANLAN R H. Wind effects on structures: An introduction to wind engineering[M]. New York: Wiley, 1986. |
[18] | VICKERY B J, KAO K H. Drag or along-wind response of slender structures[J]. Journal of the Structural Division, 1972, 98(1): 21-36. |
[19] | LIANG S, LIU S, LI Q S, et al. Mathematical model of acrosswind dynamic loads on rectangular tall buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(12): 1757-1770. |