特大型冷却塔单塔内表面风荷载三维效应及其设计取值

邹云峰等

摘要：通过刚性模型测压风洞试验获得某核电站220 m高特大型冷却塔的内表面风荷载，分析冷却塔内表面风荷载的三维效应，采用有限元计算方法对内压设计取值简化进行探讨，并对结果产生的原因进行分析.结果表明，冷却塔内表面风压并非完全沿环向、高度均匀分布，这种不均匀性以风筒施工期间的空塔尤为严重；尽管内压的大小及分布特征会对响应产生一定的影响，但由于冷却塔的风致响应中以外压占主导地位，内压对响应的贡献较小，将具有“三维效应”的内压简化为沿高度、环向不变的常数进行设计可满足安全性要求，大小可取为-0.50.研究结果可为200 m级特大型冷却塔设计内压取值和相关规范修订提供参考与依据.

关键词：冷却塔；内压；三维效应；风洞试验；设计取值

中图分类号：TU311.3；V211.7文献标识码：A

为实现循环水的冷却，冷却塔风筒顶部敞开，底部由人字柱支撑而形成风通道，故内表面也受到风荷载作用.相对外表面风荷载而言，有关冷却塔内表面风荷载的研究成果较少，中国规范也未对其取值作出相关规定，但明确指出在进行塔筒局部弹性稳定性验算时必须要考虑内压产生的应力[1-2].研究表明[3]，考虑内压效应后，环向压应力增大，结构的整体稳定性与局部弹性稳定性降低，这一趋势可能随着冷却塔高度的增加而更加显著.

孙天风等[4]通过对茂名冷却塔的实测研究发现强风作用下的内压并非沿环向均匀分布；Kasperski等[5]通过风洞试验发现内压沿环向和高度均匀分布，压力系数值接近-0.50；李鹏飞等[6]的风洞试验结果表明内压基本均匀分布，但大小与填料层透风率密切相关；张陈胜[7]和沈国辉等[8]通过CFD方法对内压分布进行了研究，结果表明内压沿高度和环向变化明显；鲍侃袁[9]通过CFD数值模拟发现塔底尾流区内侧风压急剧减少；余关鹏[10]和沈国辉等[11]通过风洞试验发现内压在底部180°圆周角急剧减少.此外，一些学者认为内压沿环向、高度不变，根据经验假定它为某一数值，例如，Diver[12]认为内压压力系数为-0.40～-0.50，Sollenberger等[13]认为取值为-0.40，Scanlan等[14]在由内外压差实测数据获得外表面风压系数时取内压为-0.40，Kawarabata等[15]认为实际设计中内压可取为-0.45.由以上综述可以看出，虽然一些学者认为内压沿环向、高度不变，但也有研究结果（包括实测、数值模拟和风洞试验）表明内压沿高度变化，沿环向分布也并不均匀，内压的分布特征及其取值尚无统一认识.此外，已有研究大多在20世纪70，80年代进行，研究对象高度大多在100 m左右[16-17].

随着中国电力事业的发展，中国冷却塔高度即将突破200 m的世界纪录，此类特大型冷却塔设计的风荷载取值既无规范指导，又无实际工程经验借鉴，因此亟需开展内压相关研究.本文以中国某核电站拟建的220 m高双曲冷却塔为例，通过风洞试验的同步测压技术获得其内表面的风压系数，对该塔的内表面风荷载三维效应进行了分析，最后采用有限元方法对其内表面风压设计取值的简化方法进行了分析.研究成果可为此类巨型冷却塔内表面风荷载取值与中国相关规范的修订提供参考和依据.

1风洞试验概况

1.1试验模型

某核电站拟建冷却塔淋水面积为20 000 m2，塔顶标高220.0 m，喉部标高169.4 m，进风口标高13.45 m，人字柱底面标高0.0 m，塔顶直径109.0 m，喉部直径103.5 m，底部直径169.9 m，风筒采用分段等厚，最小厚度在喉部断面，壁厚0.23 m，最大厚度在下环梁位置，壁厚1.4 m，由均匀分布的56对1.4 m人字柱支撑.刚性模型测压风洞试验在湖南大学HD2风洞的高速试验段进行，试验段长17.0 m，宽3.0 m，高2.5 m.试验模型采用6 mm厚的有机玻璃制作，保证模型具有足够的强度和刚度，在试验风速下不发生变形及不出现明显的振动.模型内表面在外形上与实际结构保持几何相似，几何缩尺比为1/200，外表面几何相似则由于结构壁厚太薄而无法满足，模型底部由严格几何缩尺的人字柱支撑，保证人字柱之间的空隙使得空气可以自由出入，确保真实模拟冷却塔内部空气流动特征，试验模型照片如图1所示.在模型表面共布置14层测点，每层沿环向等间距布置36个测点，共计504个测点，测点布置及圆周角定义如图2所示.

1.2风场模拟

项目厂址周边地形与GB 50009—2001《建筑结构荷载规范》规定的B类地貌类似，在湖南大学HD2风洞高速试验段模拟了B类地貌风场，转盘中心处的模拟结果如图3所示.从图3（a）可以看出，风洞中模拟的平均风速剖面与GB 50009—2001《建筑结构荷载规范》规定的B类风场基本一致，湍流度剖面也与实际大气中的情况基本一致；图3（c）给出了转盘中心50 cm高处的顺风向脉动风谱，可以看出，模拟的顺风向脉动风谱与常用的von Karman，Kaimal和Davenport等理论谱基本一致.

1.3符号定义

任意测点i处的风压系数CPi表示为：

2风荷载三维效应

冷却塔内表面风压与塔底填料层透风率大小及有无十字挡风板密切相关，考虑到风筒在施工期间塔底尚未安装填料层，此时透风率为100%（无十字挡板），而冷却塔实际运行状态下填料层透风率一般为30%（有十字挡板），因此以这两种工况的测试结果为例，对内压三维效应及其设计取值进行分析.

图4为各测层平均风压沿环向分布情况.由图4可知，当透风率为100%（无十字挡板）时（图4（a）），塔底风压系数在180°圆周角附近突然增大（最大值达-0.24），这一现象与文献[4，8]等的研究结果一致，其解释为从塔底迎风面进来的气流撞击在尾流区内壁，使得该区域的风压增大；其他高度的平均风压系数沿环向基本不变，但不同高度处的平均风压系数值略有不同，约为-0.50～-0.60.当透风率为30%（有十字挡板）时（图4（b）），塔底内压沿环向分布更为均匀，没有在180°圆周角附近发生突变，各高度平均风压系数沿环向基本不变，这是因为填料层具有“整流”作用，使塔内气流分布较空塔更为均匀，但不同高度处的平均风压系数值略有不同，约在-0.45～-0.55内微小波动.图5为各测层压力系数平均值沿高度变化曲线，从图5可以更为清晰地看出，内压并非完全沿高度均匀分布.图6为各高度截面阻力系数沿高度变化曲线，阻力系数绝对值与0偏差越远，表明该高度风压沿环向分布越不均匀，由图6可知，当透风率为100%（无十字挡板）时，这种不均匀性在塔底尤为显著.

3响应计算与结果分析

3.1有限元模型及荷载取值

采用大型通用有限元分析软件ANSYS对原型结构进行有限元分析.建模时，冷却塔筒体采用shell63壳单元模拟，人字柱采用beam188 Timoshenko梁单元模拟，支柱上端节点与风筒末节圆的有关节点位置保持一致，边界条件为人字柱底端固结.划分网格时，子午向根据模板节数划分，环向等分为人字柱对数的适当倍数，保证适当的网格密度以确保计算结果的准确性.有限元模型及1阶模态分析结果如图7所示，结构基频为0.738 Hz.

考虑到冷却塔风致响应中外表面风荷载的贡献占主导地位，因此各计算工况均考虑了外表面风荷载的作用，选取中国有肋曲线（该塔的设计风压曲线）为外表面风压曲线.

3.2静力响应

图8所示为各工况下的静力响应比较，可以看出，位移响应对内压大小及其分布特征不敏感，各工况下的位移响应基本一致，当塔底风压分布较为不均匀时（工况6），位移甚至略微偏小；相对而言，内压的分布特征对子午向内（应）力影响较小，对环向内（应）力影响较大；内压沿子午向分布的不均匀性对内（应）力的影响较小，环向分布的不均匀性影响较大；内压分布越不均匀，环向压内（应）力越大，拉内（应）力越小，且塔底风压沿环向分布的不均匀性对环向内（应）力的影响随高度的增加而减小，当内压为均匀分布时，环向压内（应）力与内压系数绝对值成正比，拉内（应）力则与之成反比.

对而言，内压沿环向分布的不均匀性对屈曲稳定影响较大，沿高度的不均匀性影响较小；当内压均匀分布时，临界风速与内压绝对值成反比.

综上所述，静力响应与屈曲稳定计算结果均表明，将具有“三维效应”的内压简化为沿高度、环向不变的常数进行设计可保证该塔的安全性，且内压取值为-0.50可满足要求.

为进一步分析内压大小及分布特征对响应的影响机理，图10给出了不考虑外压时表1中各内压工况下的响应计算结果.由图10可以看出，尽管塔底风压沿环向不均匀分布时会使得响应增大，但它引起的响应相对外压较小，例如，最大位移约为0.25 cm，不到外压的10%.因此，尽管内压的分布特征会对响应产生一定的影响，但由于冷却塔的风致响应中以外压占主导地位，内压对响应的贡献较小，故将具有“三维效应”的内压简化为沿高度、环向不变的常数进行设计可满足安全性要求.

4结论

基于刚性模型测压风洞试验获得了某核电站220 m高超大型冷却塔沿高度、环向变化的三维内表面风压系数，分析了冷却塔内表面风荷载的三维效应，采用有限元计算方法对内压设计取值简化进行了讨论，并对结果产生原因进行了分析，主要结论如下：

1）冷却塔内表面风压并非完全沿高度、环向均匀分布，风压系数沿高度方向约在0.1范围内波动，环向不均匀性以风筒施工期间的空塔塔底风压尤为严重.

2）尽管内压的大小及分布特征会对响应产生一定的影响，但由于冷却塔的风致响应中以外压占主导地位，内压对响应的贡献较小，将具有“三维效应”的内压简化为沿高度、环向不变的常数进行设计可满足安全性要求，大小可取为-0.50.

参考文献

[1]GB/T 50102-2003工业循环水冷却设计规范[S]. 北京： 中国计划出版社， 2003：21-23.

GB/T 50102-2003Code for design of cooling for industrial recalculating water[S]. Beijing： China Planning Press， 2003：21-23. （In Chinese）

[2]NDGJ5-88火力发电厂水工设计技术规定[S]. 北京： 中国电力出版社， 2006：120-122.

NDGJ5-88Technical specification for hydraulic design of thermal power plant[S]. Beijing： China Electric Power Press， 2006：120-122. （In Chinese）

[3]赵林， 李鹏飞， 葛耀君. 等效静风荷载下超大型冷却塔受力性能分析[J]. 工程力学， 2008， 25（7）：79-86.

ZHAO Lin， LI Pengfei， GE Yaojun. Numerical investigation on performance of super large cooling towers under equivalent static wind load[J]. Engineering Mechanics，2008， 25（7）： 79-86. （In Chinese）

[4]孙天凤， 周良茂. 无肋双曲线冷却塔风压分布的全尺寸测量和风洞研究[J]. 空气动力学学报， 1983， 12（4）：68-76.

SUN Tianfeng， ZHOU Liangmao. A full scale and wind tunnel study of wind pressure distribution around a ribless hyperbolic cooling tower[J]. Acta Aerodynamica Sinica， 1983，12（4）：68-76. （In Chinese）

[5]KASPERSKI M， NIEMANN H J. On the correlation of dynamic wind loads and structural response of naturaldraught cooling towers[J]. Jouranl of Wind Engineering and Industrial Aerodynamics， 1988， 30（2）：67-75.

[6]李鹏飞， 赵林， 葛耀君，等. 超大型冷却塔风荷载特性风洞试验研究[J]. 工程力学， 2008， 25（6）：60-67.

LI Pengfei， ZHAO Lin， GE Yaojun， et al. Investigation on wind load characteristics for super large cooling tower in wind tunnel[J]. Engineering Mechanics，2008， 25（6）： 60-67. （In Chinese）

[7]张陈胜. 大型双曲冷却塔风荷载的数值模拟研究[D]. 杭州： 浙江大学建筑工程学院， 2008：33-37.

ZHANG Chensheng. Numerical simulation study of wind load on large hyperbolic cooling tower[D].Hangzhou： College of Civil Engineering and Architecture， Zhejiang University， 2008：33-37. （In Chinese）

[8]沈国辉， 张陈胜， 孙炳楠， 等. 大型双曲冷却塔内表面风荷载的数值模拟[J]. 哈尔滨工业大学学报， 2011， 43（4）：104-108.

SHEN Guohui， ZHANG Chensheng， SUN Bingnan， et al. Numerical simulation of wind load on inner surface of large hyperbolic cooling tower[J]. Journal of Harbin Institute of Technology， 2011， 43（4）：104-108. （In Chinese）

[9]鲍侃袁， 沈国辉， 孙炳楠. 大型双曲冷却塔平均风荷载的数值模拟研究[J]. 空气动力学学报， 2009， 27（6）：650-655.

BAO Kanyuan， SHEN Guohui， SUN Bingnan. Numerical simulation of mean wind load on large hyperbolic cooling tower[J]. Acta Aerodynamic Sinica ， 2009， 27（6）：650-655. （In Chinese）

[10]余关鹏. 大型双曲冷却塔风荷载特性和风致干扰效应研究[D]. 杭州： 浙江大学建筑工程学院， 2010：70-80.

YU Guanpeng. Study on characteristics of wind load and windinduced interference for large hyperbolic cooling tower[D]. Hangzhou： College of Civil of Engineering and Architecture， Zhejiang University， 2010：70-80. （In Chinese）

[11]沈国辉， 余关鹏， 孙炳楠， 等. 考虑内外压共同作用的大型冷却塔风荷载分析[J].空气动力学学报， 2011， 29（4）：439-446.

SHEN Guohui， YU Guanpeng， SUN Bingnan，et al. Analysis of wind load on large hyperbolic cooling tower considering interaction between internal and external pressure[J].Acta Aerodynamic Sinica， 2011， 29（4）：439-446. （In Chinese）

[12]DIVER M. Large cooling towers the present trend[J]. Jouranl of Structural Engineer，1977， 55（10）：130-137.

[13]SOLLENBERGER N J， BILLINGTON D P. Wind loading and response of cooling towers[J]. Journal of the Structural Division， ASCE， 1980， 103（3）：601-621.

[14]SCANLAN R H， LEONARD J F. Turbulent winds and pressure effects around a rough cylinder at high reynolds number[J]. Jouranl of Wind Engineering and Industrial Aerodynamics， 1982， 9：207-236.

[15]KAWARABATA Y， NAKAE S， HARADA M. Some aspects of the wind design of cooling towers[J]. Journal of Wind Engineering and Industrial Aerodynamics， 1983， 14：167-180.

[16]邹云峰， 牛华伟， 陈政清. 基于完全气动弹性模型的冷却塔干扰效应风洞试验研究[J]. 湖南大学学报： 自然科学版， 2013， 40（12）：1-6.

ZOU Yunfeng， NIU Huawei， CHEN Zhengqing. Wind tunnel test on windinduced interference effect of cooling towers based on full aeroelastic model[J]. Journal of Hunan University： Natural Sciences， 2013， 40（12）：1-6. （In Chinese）

[17]柯世堂， 葛耀君， 赵林. 基于气弹试验大型冷却塔结构风致干扰特性分析[J]. 湖南大学学报： 自然科学版， 2010， 37（11）：18-23.