利用旋涡发生器减小大跨方形平屋盖风吸力的实验研究 Experimental study on reducing wind suction of large-span sguare flat roof using passive vortex generators 1 first-author 1 2 corresp 周蕾，博士，副研究员，从事结构及桥梁抗风研究；E-mail：lzhouau@connect.ust.hk 周蕾，博士，副研究员，从事结构及桥梁抗风研究；E-mail：lzhouau@connect.ust.hk lzhouau@connect.ust.hk 东北林业大学 土木工程学院，黑龙江 哈尔滨，150040 东北林业大学 土木工程学院，黑龙江 哈尔滨，150040 School of Civil Engineering, Northeast Forestry University, Harbin 150040, China School of Civil Engineering, Northeast Forestry University, Harbin 150040, China 香港科技大学 土木与环境工程学院，中国 香港，999077 香港科技大学 土木与环境工程学院，中国 香港，999077 Department of Civil and Environmental Engineering, Hong Kong University of Science and Technology, Hong Kong 999077, China Department of Civil and Environmental Engineering, Hong Kong University of Science and Technology, Hong Kong 999077, China 杨幼平 在强风作用下，屋盖作为风敏感部位往往最先发生破坏，如何抑制屋盖上的极值负风压是保障大跨建筑抗风安全的关键。基于此，提出一种简单有效的新型屋盖抗风气动措施。在平屋顶屋檐处安装旋涡发生器(PVG)，通过旋涡发生器的扰流作用来减小屋盖极值风压。采用风洞测压试验研究在不同风向角下PVG对屋面极值风压的气动控制效果，分析PVG对屋顶流场结构的影响并探讨PVG的工作机理。研究结果表明：PVG可以有效减小平屋盖屋面极值风压，尤其在倾斜风向角下的减压效果显著。在最不利风向角下，屋盖极值风压最大减小幅度可达73.3%。相对于长度来说，PVG高度对屋面风压的影响更明显，较密的安装间距和适宜的安装角度可以更好地降低屋盖风吸力。在安装PVG后，屋盖锥形旋涡的涡核位置和再附位置发生改变，屋顶风压低频部分能量占比大幅度减少，脉动能量减弱，大尺度旋涡减少，而小尺度的湍流造成的脉动能量成分增加。 The roof, as a wind-sensitive part, is often the first to suffer damages during strong wind events. How to mitigate the extreme negative wind pressure on the roof is the key to ensure the safety of the buildings. Consequently, a simple-effective new device, passive vortex generators(PVG), to mitigate extreme wind pressure was proposed. By the disturbing effect of the PVGs, they were installed at the eaves of flat roofs to reduce the roof extreme wind pressure. The aerodynamic control effects of PVGs on roof extreme wind pressure were presented at different wind angles. The effects of PVGs on flow field as well as their control mechanism were analyzed and discussed. The results show that the PVGs can effectively reduce the wind pressure on the flat roof, particularly at nonzero wind angles. At the most unfavorable wind angle, the maximum reduction in roof extreme wind pressure can be up to 73.3%. Compared with the length, the PVG height has a more obvious effect on mitigating the roof suction. A denser installed spacing and a suitable installed angle can better reduce the roof wind pressure. After installing the PVGs, the position of the vortex core and the reattachment position of the conical vortex are highly modified. The low-frequency portion of the wind pressure fluctuation energy is significantly reduced, and the total fluctuation energy is weakened. The large-scale vortices are then reduced, while the small-scale turbulence component is increased. 平屋盖 大跨建筑 风洞试验 极值风压 旋涡发生器 flat roof large-span building wind tunnel test extreme wind pressure passive vortex generator http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705825&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705811&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705826&type= 14.66796875 14.66793823 风灾一直是自然界最常见的灾害之一，极端风天气每年都会给经济造成巨大的破坏。其中，对住宅、商业和工业结构等大跨建筑的破坏造成的损失占风灾总损失的较大部分 [ 1 1 ] 。研究大跨建筑的抗风是十分必要的。然而对于大跨建筑来说，屋盖是建筑发生破坏的起始部位，风敏感性较强，总是受到很强的风吸力。 由于屋盖边缘的几何突变，流动分离总是发生在屋顶前缘附近，并伴随着各种流动结构的形成，其中最著名的是柱状涡和锥状涡 [ 2 2 ] 。这两种流动结构在高风速下会导致屋盖上的强风吸力，使得屋盖遭受破坏 [ 3 3 - 5 5 3-5 ] ，且相对柱状涡而言，锥状涡的作用形态更为复杂，破坏性更强 [ 6 6 ] 。一方面，屋盖破坏形式主要是负风压吸力导致的掀翻破坏，进而使得建筑的整体结构损坏，所以，降低相应部位的极值负风压是提高大跨屋盖结构抗风能力的关键。另一方面，风向角对风荷载的影响较大，在不同风向角下，屋盖的分离气泡和锥形旋涡的作用特性不同 [ 7 7 - 8 8 7-8 ] ，屋盖常在最不利风向角下遭受破坏。 气动措施是控制流动分离和漩涡脱落的有效方法。气动控制分为主动气动控制和被动气动控制两种。主动气动控制是借助外界辅助能量，在建筑物附近局部流场中施加适当的扰动形式来控制边界层流动，进而改善建筑物全局流场。马冬 [ 9 9 ] 通过风洞实验和数值模拟将等离子体流动控制应用于低矮房屋上，发现等离子体射流能对屋面的平均风压系数产生明显的影响。被动气动控制往往是通过改变建筑结构的气动外形来实现气动控制的目的，过程中不需要任何辅助能量的输入。例如，在屋盖处添加扰流板，多孔护栏及修建特殊形状的女儿墙等被动控制方式改变屋盖风场。被动控制技术较主动控制技术的优点是结构简单，无需添加额外的能量消耗装置，节省能源。KOPP等 [ 10 10 ] 研究了在屋顶安装各种护栏来减少角部面积平均荷载方面的效果。其中，扰流板和多孔连续护栏表现最好，在屋顶拐角处的峰值压力系数分别减少了44%和56%。BANKS [ 11 11 ] 在风洞中对模型采用了几种改造缓解技术，包括在屋顶安装扰流板、多孔网状护栏、圆孔多孔护栏、板栅多孔护栏、和实心护栏等，发现这些装置通过抑制气流分离，可以减缓屋顶上的旋涡作用，其中扰流板产生的效果较好。由此可见，相较于传统的被动装置，扰流板的减压效果更好，且安装简便，是优先选择的被动装置。周显鹏 [ 12 12 ] 通过采用风洞实验及CFD数值模拟的方法研究了扰流板高度和宽度对典型低矮建筑双坡屋面风压的影响。甘石等 [ 13 13 ] 利用风洞试验研究了扰流板宽度，高度和角度对双坡房屋屋面风压的影响，发现扰流板能有效地降低屋面的平均风压和峰值负压。WU [ 14 14 ] 通过现场实测试验和风洞试验，研究了扰流板对与平坡屋顶房屋的影响。然而，目前国内外对扰流板的研究多以改进传统的建筑构件为主，难免造成材料的浪费。考虑到部分节能型建筑的平屋顶上安装光伏板等集能装置，在屋盖上安装传统的扰流板对太阳能的收集有很大影响，需要开发一种实用的新型扰流装置来满足现代住房的需求。 旋涡发生器(PVG)是一种高效的边界层分离控制装置，具有结构简单、安装方便等优点。在来流风的作用下，PVG可生成一对反向旋转的流向涡，该类涡结构可有效扰乱或破坏大尺度涡 [ 15 15 ] ，因此PVG具有削弱引发极值风压的锥形涡和柱状涡的潜质，进而可减小平屋盖风吸力。MANOLESOS等 [ 16 16 ] 将旋涡发生器安装在风力涡轮机叶片机翼上进行风洞试验，发现旋涡发生器可以有效地延迟或抑制机翼上的分离现象。高峰等 [ 17 17 ] 利用数值模拟的方法研究了安装在超燃冲压发动机隔离段的多种微型旋涡发生器的流动控制性能，发现随着微型涡流发生器后缘宽度增加，隔离段的抗反压能力有较大提升，流场质量有所改善。XIN等 [ 18 18 ] 通过风洞试验的方法将旋涡发生器安装在桥梁下表面来抑制桥梁的涡激振动，发现旋涡发生器在合适的尺寸和展向安装间距下，可以完全抑制桥梁的竖向涡激振动和扭转振动。鉴于此，本文作者通过风洞试验方法，将被动式旋涡发生器(PVG)安装在大跨平屋盖建筑的屋檐上，研究旋涡发生器(PVG)对屋盖负风压的抑制效果。 1 风洞试验 1 1.1 　试验条件 2 风洞试验在东北林业大学风洞实验室进行，风洞试验段截面宽×高×长为0.8 m×1.0 m×5.0 m。风洞试验采用尖劈、地毯、粗糙元来模拟大气边界层流场。大气边界层流场参考GB 50009—2012《建筑结构荷载规范》 [ 19 19 ] 规定的B类地貌，地面粗糙度指数 α =0.16，模型顶部的紊流度 I 约为11.2%，风洞试验的参考速度 U 0 为8.6 m/s。 图1 图1 所示为风洞试验布置图。试验中大气边界层模拟试验测得的平均风速曲线和湍流强度曲线与中国规范规定的参考曲线拟合良好，如 图2 图2 所示。 1.2 　试验模型 2 在强风作用下，迎风前缘屋檐处由于流动分离产生柱状涡和锥形涡，使得屋面受到巨大的负压，导致屋盖被掀翻。本文将旋涡发生器放置在发生流动分离的屋檐位置。试验中大跨建筑长×宽×高为36 m×36 m×36 m，如 图3 图3 所示。将屋盖处的高度作为试验的参考高度，屋盖类型为平屋盖。本试验中建筑模型的几何比例为1꞉244，同时满足风洞阻塞率。屋盖的测点编号布置及来流方向如 图4 图4 所示，在屋面上共设置96个测点，在迎风前缘和角部区域适度加密。试验中风向角的间隔取15°，从0°~90°分别测量7个不同风向角下屋盖表面的风压。风压的采样频率为500 Hz，采样时间为30 s。 图5 图5 所示为旋涡发生器的几何参数，如旋涡发生器长度 L ，旋涡发生器高度 H ，旋涡发生器安装间距 R 以及旋涡发生器的安装角度 β ，旋涡发生器工作原理图见 图6 图6 。 数据单位：mm。 2 数据处理与验证 1 2.1 　数据处理方法 2 模型试验中符号约定以压力向内(压)为正，向外(吸)为负。屋盖表面第 i 个测压点所在位置的量纲一压力系数 <math id="M1"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705720&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705729&type= 8.97466660 4.48733330 为 <math display="block" id="M2"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mi>P</mi><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></mrow><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705730&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705706&type= 27.85533524 8.21266651 式中： P ( t i )为在该位置上测得的表面风压； P 0 和 P ∞ 分别为参考点处测得的平均总压和平均静压。 定义平均压力系数 <math id="M3"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705708&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705722&type= 8.55133343 3.72533321 为 <math display="block" id="M4"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mrow><mrow><mi>N</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705733&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705732&type= 26.41600037 12.61533356 定义脉动压力系数 <math id="M5"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705709&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705724&type= 7.11199999 3.72533321 为 <math display="block" id="M6"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfrac><mroot><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mo stretchy="false">[</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>t</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub><msup><mrow><mo stretchy="false">]</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705725&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705734&type= 52.23933029 10.41399956 对于极值风压的估算方法，“Cook-Mayne系数”法 [ 20 20 ] 是一种经典方法，因此，本文采用该方法估计屋盖模型极值风压。定义极值压力系数 <math id="M7"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">p</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705710&type= 7.95866632 3.72533321 为： <math display="block" id="M8"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">p</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">k</mi></mrow></msub><mo>=</mo><mi>m</mi><msub><mrow><mi>c</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705737&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705736&type= 17.86466599 3.72533321 <math display="block" id="M9"><msub><mrow><mi>c</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">d</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">1</mn><mo>+</mo><mn mathvariant="normal">0.636</mn><msub><mrow><mi>V</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705739&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705743&type= 24.12999916 3.72533321 式中： c ad 为调整系数； V x 为变量系数， V x = σ / m ； σ 为脉动压力； m 为时均压力。 选取迎风前缘区域的测压点计算迎风前缘区域的面积平均极值风压系数 <math id="M10"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">g</mi><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi mathvariant="normal">p</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">k</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705741&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705740&type= 19.47333336 3.72533321 <math display="block" id="M11"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">g</mi><mo>-</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi mathvariant="normal">p</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">k</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mrow></msub><mo>=</mo><mfrac><mrow><mfrac><mrow><mo>∑</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi mathvariant="normal">p</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">k</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>⋅</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mo>∑</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac></mrow><mrow><mn mathvariant="normal">0.5</mn><mi>ρ</mi><msup><mrow><mi>U</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705758&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705742&type= 45.97399902 16.00200081 式中： P i, peak ( t )为测压点 i 在 t 时刻的极值风压； A i 为测点 i 的所属面积。 2.2 　试验验证 2 为了验证试验结果的正确性，将不加旋涡发生器(PVG)屋盖试验数据与东京工艺大学(TPU)测压试验数据进行对比 [ 21 21 ] 。TPU试验模型与本试验模型相似，为边长16 cm的立方体平屋盖模型，缩尺比为1꞉100。TPU试验模拟日本规范中Ⅲ类地貌条件，粗糙度指数为0.2。本试验模型屋盖平均风压系数和极值压力系数与TPU试验进行对比如 图7 图7 所示。 (a) 0°风向本次试验；(b) 0°风向TPU试验 由 图7 图7 可见：0°风向角下模型屋盖平均风压分布与TPU试验结果相似，全屋盖受负压作用，其中迎风前缘所受负压最大，由迎风前缘向屋盖内侧递减，这符合平屋盖风压分布规律 [ 22 22 ] 。0°风向角下模型屋盖迎风前缘所受负压略小于TPU模型试验数据，这一方面是由于本试验模型略小于TPU试验模型，另一方面本试验模拟地貌曲线也与日本规范略有差异。总体来说，本次风洞试验结果与相关研究学者的风洞试验结果基本一致，证明本次风洞试验结果是可靠的。 3 结果与讨论 1 3.1 　旋涡发生器( PVG )的气动控制效果 2 由于屋顶转角区域所受的负压高于其内部的压力，所以划分出屋盖一侧角部区域，如 图4 图4 所示。为了研究PVG在不同风向角下的气动控制效果，将PVG安装在平屋顶模型屋檐处，其中PVG长 L =1.5 cm，高 H =1.0 cm，间距 R =2 cm，在屋盖4个转角附近对称布置，安装角度 β =60°和120°，共安装24个PVG，观察屋顶的极值风压变化情况。 图8 图8 所示为安装PVG和未安装PVG情况下的极值压力系数和角部面积平均极值风压系数的变化，风向角范围为0°~90°(间隔15°)。从 图8 图8 可以看出：未安装PVG时，屋盖大部分极值风压发生在转角区域。屋顶在0°风向和90°风向下受到的风吸力小于其他风向角下受到的风吸力，这说明屋顶对倾斜风向更为敏感。极值压力系数曲线呈对称分布，在45°风向角下测点11位置处的极值压力系数最小，此时极值压力系数为-4.72，45°风向角为最不利风向角。角部面积平均极值风压系数在30°~60°风向角下所受的风吸力较大，均超过-2.78。 在安装PVG后，屋盖在倾斜风向角下的极值风压有明显的降低，极值压力系数控制在-2.00以下，最大减小幅度可达59%，角部面积平均极值风压系数在倾斜风向角下也明显降低，最大减小幅度达到57%。这是由于在30°~60°倾斜风向角下，屋盖迎风角部区域受锥形旋涡的影响较大，使得此处产生了较大的风吸力，而安装在角部区域的PVG产生的顺流向涡干扰了锥形旋涡的形成，掺混了迎风前缘边界层内的气流，减小了该区域的风压。可见旋涡发生器(PVG)可以有效减小平屋盖的极值风压。 3.2 　旋涡发生器( PVG )几何因素的影响 2 为了研究PVG的几何因素对平屋盖风压的影响，测量在最不利风向角下(45°)，不同参数的PVG安装在试验模型上时屋盖的风压，试验具体工况见 表1 表1 。 工况 高度 H /cm 长度 L /cm 间距 R /cm 安装角度 β /(°) 工况1 0.5，0.7，1.0 1.5 2.0 60 工况2 1.0 1.1，1.3，1.5 2.0 60 工况3 1.0 1.5 2.0，3.0，4.0 60 工况4 1.0 1.5 2.0 30，60，90 3.2.1　高度 H 的影响 3 PVG高度 H 对屋盖风压系数影响的工况见 表1 表1 中工况1，其结果如 图9 图9 所示。由 图9 图9 可见：加入PVG后，迎风前缘测点行Row 1内风吸力急剧降低，最大降低幅度达59%(测点2)，相应地风压脉动也急剧减弱。测点列Column 2内存在很大的风吸力突变(测点48位置平均风压系数达到-2.0)，加入PVG后风压突变得到抑制。这是由于角部的气流遇PVG后在其附近产生大量漩涡，这些漩涡对分离气流产生强烈干扰，同时抑制了前缘锥形漩涡的形成。平均风压系数对PVG高度比较敏感，随PVG高度增加，其减压效果越来越好。 (a) Row 1测点，平均压力系数；(b) Row 1测点，脉动压力系数；(c) Column 2测点，平均压力系数； (d) Column 2测点，脉动压力系数 3.2.2　长度 L 的影响 3 PVG长度 L 对屋盖风压系数影响的工况见 表1 表1 中工况2，其结果如 图10 图10 所示。由图 10 10 (a)和 10 10 (b)可见：改变PVG长度对迎风前缘(Row 1)气流的干扰效果与其高度相比不那么明显。相对于另一侧迎风前缘而言(Column 2)， L =1.5 cm的减压效果与 L =1.3 cm的相比差别不大，然而当 L < 1.3 cm时，PVG对平均风压和脉动风压的抑制效果减弱。由此可见，在进行PVG设计时，为达到近似的减压效果，PVG的长度可以考虑在一定范围内有所调整。 (a) Row 1测点，平均压力系数；(b) Row 1测点，脉动压力系数；(c) Column 2测点，平均压力系数； (d) Column 2测点，脉动压力系数 3.2.3　间距 R 的影响 3 PVG安装间距 R 对屋盖风压系数影响的工况见 表1 表1 中工况3，其结果如 图11 图11 所示。由 图11 图11 可见：安装间距对屋盖前缘的风压影响显著，特别是对于脉动风压而言。在接近角部附近， R 的影响作用更加明显，即随安装间距增加，PVG对负风压即风压脉动的抑制作用逐渐降低。这是由于在迎风转角区域是锥形涡的形成区，较密的布置更有利于干扰分离流，抑制锥形涡的形成。 3.2.4　安装角度 β 的影响 3 PVG安装角度 β 对屋盖风压系数影响的工况见 表1 表1 中工况4，其结果如 图12 图12 所示。由 图12 图12 可见：3种安装角度的PVG均能对分离流产生干扰，减小屋盖迎风前缘的风吸力。然而，对比不同安装角度可以发现， β =90°时，PVG对前缘两侧靠近角部位置的减压效果不如PVG倾斜时的效果显著，其原因是：PVG倾斜时，不仅能够对气流产生干扰，而且倾斜的角度对撞向钝体边缘的气流起到疏导的作用，使得分离流沿着PVG倾斜方向偏转，从而弱化了来流风的分离作用，进而抑制了锥形涡的形成。 3.3 　最不利风向角下的气动控制 2 由于在最不利风向角下，屋盖所受的极值风压最大，所以控制屋盖在最不利风向角下的极值风压是提高屋盖抗风能力的关键。由前面分析可得PVG的最优参数，选取 L =1.5 cm， H =1.0 cm， R =2.0 cm， β =60°的PVG，将其安装在模型屋檐处，探究在最不利风向角下(45°)，PVG的气动控制效果。 图13 图13 所示为45°风向角下屋盖风压系数的等值线图。从 图13 图13 可以看出：屋顶表面风荷载以风吸力为主，在靠近迎风前缘位置风吸力较大，靠近内部风吸力逐渐减弱。加入PVG后，屋顶风压分布发生明显的变化。与未安装PVG的情况相比，安装PVG时屋顶明显表现出较低的极值压力系数且位于边缘和角部区域的压力降幅高于内部区域。角部的极值风压系数减少效果最为显著，最大减小幅度可达73.3%(从-4.5降低到-1.2)，平均压力系数降低了64.0%(从-2.5降低到-0.9)，脉动压力系数降低了77.8%(从0.9降低到0.2)。与未安装PVG相比，屋顶转角和边缘的压力梯度明显减小，锥形旋涡影响区域减小，尤其是迎风边缘两侧的压力分布基本平稳均匀，其他区域并没有出现风压增大等突变现象。可见，在最不利风向角下，旋涡发生器(PVG)可以有效减小平屋盖的极值风压。 (a) Row 1测点，平均压力系数；(b) Row 1测点，脉动压力系数；(c) Column 2测点，平均压力系数； (d) Column 2测点，脉动压力系数 (a) Row 1测点，平均压力系数；(b) Row 1测点，脉动压力系数；(c) Column 2测点，平均压力系数； (d) Column 2测点，脉动压力系数 (a) 未安装PVG的平均压力系数 C p mean ；(b) 安装PVG后的平均压力系数 C p mean ；(c) 未安装PVG的脉动压力系数 C p rms ；(d) 安装PVG后的脉动压力系数 C p rms ；(e) 未安装PVG的极值压力系数 C p peak ；(d) 安装PVG后的极值压力系数 C p peak 3.4 　脉动风压功率谱 2 图14 图14 所示为屋盖极值风吸力较大的风向角下最不利测点在安装PVG前后的脉动风功率谱变化。其中， f∙H / U 为量纲一频率， f∙S ( f )/ σ 2 为量纲一自功率谱， S ( f )为测点风压自谱， σ 2 为脉动风压的方差。风压功率谱密度选用对数刻度，同时采取适当的平均处理平滑谱估计的结果，这样可以从整体上评价PVG对屋盖动态脉动风压的减缓效果。对于45°风向角下测点1( 图14 图14 (a))，加入PVG后，来流风低频部分( f·H / U ≈0.1)能量占比大幅度减少，说明分离流中大尺度漩涡的占比减小，这对锥形旋涡的产生和运动起着重要的作用。对于45°风向角下测点2( 图14 图14 (b))，加入PVG后频谱整体能量降低，并且中低频部分能量降低较多，这说明大尺度漩涡携带的能量降低，脉动强度减弱。由 图 14 图 14 (c)可见：量纲一频率在10 -2 ~10 -1 之间时，未安装PVG的曲线为凸曲线，安装PVG后，曲线呈凹曲线，这说明低频部分整体脉动能量的降低。在高频区，屋盖脉动风能量增加，这说明PVG产生的小尺度漩涡携带高频脉动，与屋顶迎风前缘部位的来流风相互作用，消耗了一部分能量，同时扰乱和破坏了大尺度旋涡(见 图14 图14 (d))。 3.5 　旋涡发生器（ PVG ）对涡核的影响 2 在倾斜风向角下屋盖的极值风吸力主要是由迎风前缘的锥形旋涡引起的。因此，研究PVG对锥形旋涡的涡核位置和再附区域的影响是十分必要的。KAWAI [ 23 23 ] 通过垂直于迎风面前缘各截面的最大平均风压和最大脉动风压的位置来估计涡旋的轴线和锥形涡旋的再附范围。本文采用文献[ 23 23 ]中的方法来确定45°风向角下锥形旋涡的位置参数。 图15 图15 所示为45°风向角下锥形旋涡的变化范围。由 图15 图15 可见：找到屋盖每个横截面(垂直于 Y 轴)和纵截面(垂直于 X 轴)的最大平均风压系数和最大脉动风压系数所在的测点坐标( X ， Y )，对这些测点坐标进行曲线拟合(有95%的置信区间)，得到锥形涡旋的涡核曲线和再附曲线。从 图15 图15 可以看出：安装PVG后，前缘两侧锥形涡旋的再附区域明显减小，左侧涡核位置也向边缘移动。说明PVG产生的反向涡流破坏了锥形旋涡的结构。再附曲线和前缘之间的范围逐渐减小，这意味着锥形旋涡的影响范围逐渐减小，因此，PVG具有削弱引发极值风压的锥形旋涡的作用。 (a) 45°风向角下测点1 PSD；(b) 45°风向角下测点2 PSD；(c) 30°风向角下测点11 PSD；(d) 60°风向角下测点2 PSD (a) 未安装PVG；(b) 安装PVG后 4 结论 1 1) 在低矮平屋盖屋檐位置安装PVG，可以有效减少屋盖的极值风压，特别对于倾斜风向的气动控制效果最好。在最不利风向角下，PVG对于极值风压的控制程度最大可达73.3%。 2) 旋涡发生器的参数对屋盖风压的影响较大，与长度( L )相比，高度( H )的增加对锥形漩涡的干扰作用更强，为了更好地降低角部区域极值风压，PVG宜紧密布置。在适宜倾斜角下，PVG不仅能够干扰分离流，而且对迎风前缘的气流起到疏导的作用。 3) 安装旋涡发生器(PVG)后，来流风中低频部分能量占比大幅度减少，大尺度漩涡的脉动能量减弱，原因是来流中由PVG激发的小尺度的湍流成分大大削弱了大尺度湍流的能量。 4) 在45°风向角下，旋涡发生器(PVG)破坏了锥形旋涡的结构，改变了锥形旋涡的形成位置，从而减小的锥形旋涡的影响范围。 参考文献 HOLMES J D , PATON C , KERWIN R . Wind loading of structures [M]. London : Taylor & Francis , 2007 : 1 - 392 . Wind loading of structures TIELEMAN H W . Wind tunnel simulation of wind loading on low-rise structures: a review [J]. Journal of Wind Engineering and Industrial Aerodynamics , 2003 , 91 ( 12/13/14/15 ): 1627 - 1649 . Wind tunnel simulation of wind loading on low-rise structures: a review KRISHNA P . Wind loads on low rise buildings: a review [J]. Journal of Wind Engineering and Industrial Aerodynamics , 1995 , 54/55 : 383 - 396 . Wind loads on low rise buildings: a review LIN J X , SURRY D . The variation of peak loads with tributary area near corners on flat low building roofs [J]. Journal of Wind Engineering and Industrial Aerodynamics , 1998 , 77/78 : 185 - 196 . The variation of peak loads with tributary area near corners on flat low building roofs 陈学锐 , 顾志福 , 李燕 . 锥形涡诱导下建筑物顶面风荷载 [J]. 力学学报 , 2007 , 39 ( 5 ): 655 - 660 . 锥形涡诱导下建筑物顶面风荷载 CHEN Xuerui , GU Zhifu , LI Yan . Conical vortex induced wind loading on the roof of a building [J]. Chinese Journal of Theoretical and Applied Mechanics , 2007 , 39 ( 5 ): 655 - 660 Conical vortex induced wind loading on the roof of a building 孙瑛 , 苏宁 , 武岳 . 锥状涡作用下大跨度平屋盖表面脉动风压谱模型研究 [J]. 土木工程学报 , 2014 , 47 ( 1 ): 88 - 98 . 锥状涡作用下大跨度平屋盖表面脉动风压谱模型研究 SUN Ying , SU Ning , WU Yue . Modeling of conical vortex induced fluctuating wind pressure spectra on large-span flat roofs [J]. China Civil Engineering Journal , 2014 , 47 ( 1 ): 88 - 98 . Modeling of conical vortex induced fluctuating wind pressure spectra on large-span flat roofs 董欣 , 叶继红 , 丁洁民 . 低矮建筑表面破坏性旋涡及其诱导的风压特性研究综述 [J]. 建筑结构学报 , 2018 , 39 ( 1 ): 1 - 10 . 低矮建筑表面破坏性旋涡及其诱导的风压特性研究综述 DONG Xin , YE Jihong , DING Jiemin . Review of destructive vortices and vortex-induced wind pressures on low-rise buildings [J]. Journal of Building Structures , 2018 , 39 ( 1 ): 1 - 10 . Review of destructive vortices and vortex-induced wind pressures on low-rise buildings 秦乐 , 田玉基 . 某大跨平屋盖表面风荷载特性研究 [J]. 建筑结构 , 2012 , 42 ( 2 ): 160 - 164 . 某大跨平屋盖表面风荷载特性研究 QIN Le , TIAN Yuji . Research on wind load on a large-span plane roof [J]. Building Structure , 2012 , 42 ( 2 ): 160 - 164 . Research on wind load on a large-span plane roof 马冬 . 低矮结构的风压控制及抗风可靠度研究 [D]. 长沙 : 湖南大学 , 2017 : 1 - 79 . 低矮结构的风压控制及抗风可靠度研究 MA Dong . Study on the control of wind pressure and the wind-resistant reliability for low structure [D]. Changsha : Hunan University , 2017 : 1 - 79 . Study on the control of wind pressure and the wind-resistant reliability for low structure KOPP G A , MANS C , SURRY D . Wind effects of parapets on low buildings: part 4. Mitigation of corner loads with alternative geometries [J]. Journal of Wind Engineering and Industrial Aerodynamics , 2005 , 93 ( 11 ): 873 - 888 . Wind effects of parapets on low buildings: part 4. Mitigation of corner loads with alternative geometries BANKS D . The suction induced by conical vortices on low-rise buildings with flat roofs [D]. Fort Collins, CO, USA : Colorado State University , 2000 : 1 - 24 . The suction induced by conical vortices on low-rise buildings with flat roofs 周显鹏 . 水平悬挑女儿墙对低矮双坡屋面风压的影响 [D]. 泉州 : 华侨大学 , 2008 : 1 - 104 . 水平悬挑女儿墙对低矮双坡屋面风压的影响 ZHOU Xianpeng . Influence to wind pressures on the roof of low-rise gable-roofed buildings with cantilevered-parapets [D]. Quanzhou : Huaqiao University , 2008 : 1 - 104 . Influence to wind pressures on the roof of low-rise gable-roofed buildings with cantilevered-parapets 甘石 , 李钢 , 李宏男 . 扰流板减小低矮房屋屋面风压试验研究 [J]. 土木工程学报 , 2018 , 51 ( 6 ): 91 - 102 . 扰流板减小低矮房屋屋面风压试验研究 GAN Shi , LI Gang , LI Hongnan . Study on wind pressure on the roof of low-rise building alleviated by spoiler [J]. China Civil Engineering Journal , 2018 , 51 ( 6 ): 91 - 102 . Study on wind pressure on the roof of low-rise building alleviated by spoiler WU Fuqiang , SARKAR P P , MEHTA K C . Full-scale study of conical vortices and roof corners pressures [J]. Wind & Structures , 2001 , 4 ( 2 ): 131 - 146 . Full-scale study of conical vortices and roof corners pressures ZHANG Hongfu , XIN Dabo , OU Jinping . Wake control using spanwise-varying vortex generators on bridge decks: a computational study [J]. Journal of Wind Engineering and Industrial Aerodynamics , 2019 , 184 : 185 - 197 . Wake control using spanwise-varying vortex generators on bridge decks: a computational study MANOLESOS M , VOUTSINAS S G . Experimental investigation of the flow past passive vortex generators on an airfoil experiencing three-dimensional separation [J]. Journal of Wind Engineering and Industrial Aerodynamics , 2015 , 142 : 130 - 148 . Experimental investigation of the flow past passive vortex generators on an airfoil experiencing three-dimensional separation 高峰 , 王旭东 , 张涵 . 非对称来流下微型涡流发生器边界层流动控制 [J]. 空军工程大学学报(自然科学版) , 2016 , 17 ( 6 ): 6 - 11 . 非对称来流下微型涡流发生器边界层流动控制 GAO Feng , WANG Xudong , ZHANG Han . A research on boundary layer control of asymmetric supersonic flow past micro-vortex generator [J]. Journal of Air Force Engineering University(Natural Science Edition) , 2016 , 17 ( 6 ): 6 - 11 . A research on boundary layer control of asymmetric supersonic flow past micro-vortex generator XIN Dabo , ZHANG Hongfu , OU Jinping . Experimental study on mitigating vortex-induced vibration of a bridge by using passive vortex generators [J]. Journal of Wind Engineering and Industrial Aerodynamics , 2018 , 175 : 100 - 110 . Experimental study on mitigating vortex-induced vibration of a bridge by using passive vortex generators GB 50009—2012 . 建筑结构荷载规范 [S]. GB 50009—2012 . Load code for the design of building structures [S]. COOK N J , MAYNE J R . A refined working approach to the assessment of wind loads for equivalent static design [J]. Journal of Wind Engineering and Industrial Aerodynamics , 1980 , 6 ( 1 ): 125 - 137 . A refined working approach to the assessment of wind loads for equivalent static design Tokyo Polytechnic University . Aerodynamic database for isolated low-rise building [EB/OL]. [ 2007-01-01 ]. http: //www.wind.arch.t-kougei.ac.jp/system/eng/contents/code/tpu http://www.wind.arch.t-kougei.ac.jp/system/eng/contents/code/tpu . Aerodynamic database for isolated low-rise building 董欣 , 丁洁民 . 分离泡诱导下平屋盖表面风压特性研究 [J]. 建筑结构学报 , 2019 , 40 ( 5 ): 41 - 49 . 分离泡诱导下平屋盖表面风压特性研究 DONG Xin , DING Jiemin . Wind pressure characteristics on a flat roof induced by separation bubble [J]. Journal of Building Structures , 2019 , 40 ( 5 ): 41 - 49 . Wind pressure characteristics on a flat roof induced by separation bubble KAWAI H . Structure of conical vortices related with suction fluctuation on a flat roof in oblique smooth and turbulent flows [J]. Journal of Wind Engineering and Industrial Aerodynamics , 1997 , 69/ 70 / 71 : 579 - 588 . Structure of conical vortices related with suction fluctuation on a flat roof in oblique smooth and turbulent flows 10.11817/j.issn.1672-7207.2023.05.018 TU352.2 A 1672-7207（2023）05-1839-12 张洪福, 刘婷婷, 周蕾. 利用旋涡发生器减小大跨方形平屋盖风吸力的实验研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1839-1850. 张洪福, 刘婷婷, 周蕾. 利用旋涡发生器减小大跨方形平屋盖风吸力的实验研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1839-1850. ZHANG Hongfu, LIU Tingting, ZHOU Lei. Experimental study on reducing wind suction of large-span sguare flat roof using passive vortex generators[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1839-1850. ZHANG Hongfu, LIU Tingting, ZHOU Lei. Experimental study on reducing wind suction of large-span sguare flat roof using passive vortex generators[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1839-1850. 国家自然科学基金资助项目(51908107) Project(51908107) supported by the National Natural Science Foundation of China 2022-08-14 2022-10-18 2023-05-26 2023-06-20 38910307 中南大学学报(自然科学版) 5 54