留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

考虑材料非线性的RC双曲冷却塔风致破坏过程

张军锋 裴昊 朱冰 刘庆帅

张军锋, 裴昊, 朱冰, 刘庆帅. 考虑材料非线性的RC双曲冷却塔风致破坏过程[J]. 工程力学, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
引用本文: 张军锋, 裴昊, 朱冰, 刘庆帅. 考虑材料非线性的RC双曲冷却塔风致破坏过程[J]. 工程力学, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
Jun-feng ZHANG, Hao PEI, Bing ZHU, Qing-shuai LIU. WIND-INDUCED FAILURE PROCESS OF RC HYPERBOLIC COOLING TOWERS CONSIDERING THE MATERIAL NONLINEARITY[J]. Engineering Mechanics, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
Citation: Jun-feng ZHANG, Hao PEI, Bing ZHU, Qing-shuai LIU. WIND-INDUCED FAILURE PROCESS OF RC HYPERBOLIC COOLING TOWERS CONSIDERING THE MATERIAL NONLINEARITY[J]. Engineering Mechanics, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324

考虑材料非线性的RC双曲冷却塔风致破坏过程

doi: 10.6052/j.issn.1000-4750.2020.05.0324
基金项目: 国家自然科学基金项目(51508523)
详细信息
    作者简介:

    裴 昊(1995−),男,河南郑州人,硕士生,主要从事结构抗风研究(E-mail: 1315491065@qq.com)

    朱 冰(1995−),男,江西九江人,硕士生,主要从事结构抗风研究(E-mail: 535453402@qq.com)

    刘庆帅(1996−),男,河南新乡人,硕士生,主要从事结构抗风研究(E-mail: 727428364@qq.com)

    通讯作者: 张军锋(1983−),男,河南平顶山人,副教授,博士,主要从事结构和桥梁抗风研究(E-mail: brilliantshine@163.com)

WIND-INDUCED FAILURE PROCESS OF RC HYPERBOLIC COOLING TOWERS CONSIDERING THE MATERIAL NONLINEARITY

  • 摘要: 为明确RC双曲冷却塔的极限静风荷载及其风致破坏过程,以一座大型冷却塔为例,采用ABAQUS中分层壳单元模拟塔筒的混凝土和双向钢筋网,在线性计算和非线性计算验证的基础上,采用弥撒开裂模型和双线性模型模拟混凝土和钢筋的非线性受力特征,分析了自重+风荷载(工况Ⅰ)和自重+冬温+风荷载(工况Ⅱ)两种工况下的静风破坏过程,并结合荷载位移曲线、内力、应变和应力的发展和分布探究了其破坏机理。结果表明:对于工况Ⅰ,风荷载系数为λ=1.384时,首先在迎风子午线相对高度hs/Hs=0.37位置因子午向受拉出现环向裂缝且为贯通裂缝,之后整个迎风区的环向裂缝不断增加并沿环向扩展,迎风区内表面和侧风区外表面也因环向弯矩而出现子午向裂缝,开裂位置的钢筋应力迅速增加,荷载-位移曲线也由线性进入非线性并且迎风区和侧风区位移迅速增加,迎风区和侧风区内力也表现出明显的重分布特征,最后因混凝土持续开裂和喉部区域双向钢筋屈服而破坏,极限风荷载系数为λ=2.007;对于工况Ⅱ,温度效应的计入使λ=1.0时即在侧风区上部首先出现子午向裂缝,子午向裂缝开展也较工况Ⅰ略严重,但因温度的弯矩效应有限,其荷载位移曲线和工况Ⅰ基本一致,且仍以迎风区环向开裂最为显著,结构破坏依然由风荷载控制,最终的极限荷载荷载系数为λ=1.842。
  • 图  1  剪切保留模型

    Figure  1.  Shear retention model

    图  2  拉伸硬化模型

    Figure  2.  Tension stiffening model

    图  3  钢筋双线性模型

    Figure  3.  Bilinear model for reinforcement

    图  4  冷却塔特征尺度示意

    Figure  4.  Characteristic dimensions of cooling tower

    图  5  塔筒结构参数

    Figure  5.  Structural parameters of tower shell

    图  6  分层壳示意图

    Figure  6.  Layered shell element

    图  7  Port Gibson塔的荷载位移曲线

    Figure  7.  Load displacement curves of Port Gibson tower

    图  8  自重作用下塔筒内力分布

    Figure  8.  Internal forces of shell under gravity

    图  9  冬温荷载作用下塔筒内力分布

    Figure  9.  Internal forces of tower shell under gravity

    图  10  风荷载作用下塔筒内力分布

    Figure  10.  Internal forces of tower shell under wind load

    图  11  塔筒横向截面变形

    Figure  11.  Deformed shapes of tower shell: sliced views

    图  12  两种荷载工况所得荷载位-移曲线

    Figure  12.  Load-displacement curves of two load cases

    图  13  塔筒内力分布

    Figure  13.  Distribution of internal forces of tower shell

    图  14  塔筒混凝土最大主应变(仅给出εt=7.354×10−5~5.2×10−3范围内的拉应变)

    Figure  14.  Maximum principal strain of concrete in tower shell

    图  15  塔筒外侧混凝土和钢筋应力分布

    Figure  15.  Distribution of stress of exterior concrete and steel in tower shell

    表  1  两种荷载工况结果

    Table  1.   Results of two load cases

    特征结果G+λWG+T+λW计入温度效应后的结果偏差/(%)
    开裂荷载系数1.384(A)1.000(A')−28
    屈服荷载系数1.876(B)钢筋未屈服
    极限荷载系数2.007(C)1.842(C')−8
    极限位移0.836 m0.531 m−26
    注:括号中的字母如图12中所示。
    下载: 导出CSV
  • [1] 武际可. 大型冷却塔结构分析的回顾与展望[J]. 力学与实践, 1996, 18(6): 1 − 5.

    Wu Jike. Review and prospect about structural analysis of cooling towers [J]. Mechanics in Engineering, 1996, 18(6): 1 − 5. (in Chinese)
    [2] 张军锋, 葛耀君, 赵林. 渡桥电厂冷却塔风毁事故介绍和分析[J]. 工业建筑, 2013, 43(2): 153 − 158.

    Zhang Junfeng, Ge Yaojun, Zhao Lin. Introductions and analyses of collapse of cooling tower for Ferrybridge power station [J]. Industrial Construction, 2013, 43(2): 153 − 158. (in Chinese)
    [3] 张军锋, 葛耀君, 赵林. 双曲冷却塔结构特性新认识[J]. 工程力学, 2013, 30(6): 67 − 76.

    Zhang Junfeng, Ge Yaojun, Zhao Lin. New perceptions on the structure behavior of hyperboloidal cooling towers [J]. Engineering Mechanics, 2013, 30(6): 67 − 76. (in Chinese)
    [4] Medwadowski S J. Buckling of concrete shells: An overview [J]. Journal of the International Association for Shell and Spatial Structures, 2004, 45(144): 51 − 63.
    [5] Der T J, Fidler R. A model study of the buckling behaviour of hyperbolic shells [J]. ICE (Institution of Civil Engineers) Proceedings, London, England, 1968, 41(1): 105 − 118.
    [6] Mungan I. Buckling of reinforced concrete cooling tower shells: BSS approach [J]. ACI Structural Journal, 1982, 79(5): 387 − 391.
    [7] Abel J F, Billington D P, Nagy D A, Wiita D C. Buckling of cooling towers [J]. Journal of the Structural Division, ASCE, 1982, 108(10): 2162 − 2174.
    [8] 孙琪超. 自然通风冷却塔整体稳定性模型实验及有限元分析研究[D]. 杭州: 浙江大学, 2017.

    Sun Qichao. Model experiment and FEM analysis of global stability of natural draft cooling tower [D]. Hangzhou: Zhejiang University, 2017. (in Chinese)
    [9] Mang H A, Floegl H, Trappel F, et al. Wind-loaded reinforced-concrete cooling towers: buckling or ultimate load? [J]. Engineering Structures, 1983, 5(3): 163 − 180. doi:  10.1016/0141-0296(83)90014-7
    [10] Milford R V, Schnobrich W C. The effect of cracking on the ultimate strength of reinforced concrete cooling towers [J]. Engineering Structure, 1986, 8(1): 46 − 54. doi:  10.1016/0141-0296(86)90019-2
    [11] Mahmoud B E H, Gupta A K. Inelastic large displacement behavior and buckling of cooling tower [J]. Journal of Structural Engineering, 1995, 121(6): 981 − 985. doi:  10.1061/(ASCE)0733-9445(1995)121:6(981)
    [12] Noh H C. Nonlinear behavior and ultimate load bearing capacity of reinforced concrete natural draught cooling tower shell [J]. Engineering Structures, 2006, 28(3): 399 − 410. doi:  10.1016/j.engstruct.2005.08.016
    [13] Waszczyszyn Z, Pabisek E, Pamin J, Radwańska M. Nonlinear analysis of a RC cooling tower with geometrical imperfections and a technological cut-out [J]. Engineering Structures, 2000, 22(5): 480 − 489. doi:  10.1016/S0141-0296(98)00122-9
    [14] Meschke G, Mang H A. Finite element analyses of cracked cooling tower shell [J]. Journal of Structural Engineering, 1991, 117(9): 2620 − 2638. doi:  10.1061/(ASCE)0733-9445(1991)117:9(2620)
    [15] Noh S Y, Krätzig W B, Meskouris K. Numerical simulation of serviceability, damage evolution and failure of reinforced concrete shells [J]. Computers and Structures, 2003, 81(8): 843 − 857.
    [16] 朱冰. 考虑钢筋混凝土材料非线性的冷却塔极限风荷载分析[D]. 郑州: 郑州大学, 2020.

    Zhu Bing. Ultimate wind load analysis of cooling tower considering material nonlinearity of reinforced concrete [D]. Zhengzhou: Zhengzhou University, 2020. (in Chinese)
    [17] 石俊彪. 自然通风冷却塔非线性有限元分析[D]. 杭州: 浙江大学, 2008.

    Shi Junbiao. Nonlinear finite element analysis of natural draught cooling tower [D]. Hangzhou: Zhejiang University, 2008. (in Chinese)
    [18] 赵林, 李鹏飞, 葛耀君. 等效静风荷载下超大型冷却塔受力性能分析[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)
    [19] 赵海, 白国良, 徐亚洲, 等. 初始缺陷对超大型双曲冷却塔极限承载力的影响[J]. 西安建筑科技大学学报(自然科学版), 2011, 43(6): 777 − 783.

    Zhao Hai, Bai Guoliang, Xu Yazhou, et al. The influence of initial imperfections on the ultimate bearing capacity of super-large hyperbolic cooling towers [J]. Journal of Xi'an University of Architecture & Technology (Natural Science Edition), 2011, 43(6): 777 − 783. (in Chinese)
    [20] 刘家宝. 混凝土薄壳结构非线性稳定性分析研究[D]. 北京: 中国建筑科学研究院, 2017.

    Liu Jiabao. Nonlinear stability analysis of concrete shell structures [D]. Beijing: China Academy of Building Research, 2017. (in Chinese)
    [21] Yu Q Q, Gu X L, Li Y, Lin F. Collapse mechanism of reinforced concrete superlarge cooling towers subjected to strong winds [J]. Journal of Performance of Constructed Facilities, 2017, 31(6): 04017101.
    [22] 吴鸿鑫, 柯世堂, 王飞天. 基于CFD与LS-DYNA耦合技术的风致冷却塔倒塌破坏分析[C]// 第28届全国结构工程学术会议论文集(第I册), 中国, 南昌, 2019: 307 − 310.

    Wu Kaixin, Ke Shitang, Wang Feitian. Wind-indeued collapse analysis of cooling tower based on CFD and LS-DYNA [C]// Proceedings of the 28th National Academic Conference on Structual Engineering (Volume I), Nanchang China, 2019: 307 − 310. (in Chinese)
    [23] 王飞天, 柯世堂, 王晓海, 等. 风致特大型冷却塔结构连续性倒塌分析[C]. 第28届全国结构工程学术会议论文集(第I册), 中国, 南昌, 2019: 274 − 277.

    Wang Feitian, Ke Shitang, Wang Xiaohai, et al. Progressive collapse analysis of super large cooling tower under wind load [C]. 28th National Academic Conference on Structual Engineering, China Nanchang (Volume I), Nanchang, China, 2019: 274 − 277. (in Chinese)
    [24] 聂建国, 王宇航. ABAQUS中混凝土本构模型用于模拟结构静力行为的比较研究[J]. 工程力学, 2013, 30(4): 59 − 67, 82.

    Nie Jianguo, Wang Yuhang. Comparison study of constitutive model of concrete in ABAQUS for static analysis of structures [J]. Engineering Mechanics, 2013, 30(4): 59 − 67, 82. (in Chinese)
    [25] Lu X Z, Ye L P, Teng J G, et al. Meso-scale finite element model for FRP sheets/plates bonded to concrete [J]. Engineering Structures, 2005, 4(27): 564 − 575.
    [26] GB/T 50102−2014, 工业循环水冷却设计规范[S]. 北京: 中国计划出版社, 2014.

    GB/T 50102−2014, Code for design of cooling for industrial recalculating water [S]. Beijing: Chinese Planning Press, 2014. (in Chinese)
    [27] 张军锋, 朱冰, 丁玉玺. 冷却塔塔筒荷载效应和简化设计方法—各荷载单独作用分析[J]. 防灾减灾工程学报, 2018, 38(1): 161 − 167.

    Zhang Junfeng, Zhu Bing, Ding Yuxi. Load effects and simplified design procedure of the shell of hyperbolic cooling tower-Effects of every single load [J]. Journal of Disaster Prevention and Mitigation Engineering, 2018, 38(1): 161 − 167. (in Chinese)
    [28] 张军锋, 葛耀君. 双曲冷却塔温度效应分析[J]. 华中科技大学学报(自然科学版), 2012, 40(11): 63 − 67.

    Zhang Junfeng, Ge Yaojun. Study on the thermal effects of hyperboloidal cooling towers [J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2012, 40(11): 63 − 67. (in Chinese)
    [29] GB 50010−2018, 混凝土结构设计规范 [S]. 北京: 中国建筑工业出版社, 2018.

    GB 50010−2018, Code for design of concrete structures [S]. Beijing: China Architecture & Building Press, 2018. (in Chinese)
    [30] Zhang J F, Ge Y J, Zhao L, et al. Wind induced dynamic responses on hyperbolic cooling tower shells and the equivalent static wind load [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 169: 280 − 289. doi:  10.1016/j.jweia.2017.08.002
  • [1] 毕继红, 霍琳颖, 乔浩玥, 赵云.  单向受拉状态下的钢纤维混凝土本构模型 . 工程力学, 2020, 37(6): 155-164. doi: 10.6052/j.issn.1000-4750.2019.08.0443
    [2] 尹世平, 盛杰, 贾申, 叶桃.  TRC加固RC梁的弯曲疲劳破坏过程和应变发展的试验研究 . 工程力学, 2015, 32(增刊): 142-148. doi: 10.6052/j.issn.1000-4750.2014.04.S010
    [3] 顾祥林, 付武荣, 汪小林, 洪丽.  混凝土材料与结构破坏过程模拟分析 . 工程力学, 2015, 32(11): 9-17. doi: 10.6052/j.issn.1000-4750.2014.07.ST09
    [4] 张军锋, 葛耀君, 赵林, 柯世堂.  双曲冷却塔子午线型对结构受力性能的影响 . 工程力学, 2014, 31(3): 187-196. doi: 10.6052/j.issn.1000-4750.2012.10.0779
    [5] 张军锋, 葛耀君, 赵林, 柯世堂.  双曲冷却塔表面三维绕流特性及风压相关性研究 . 工程力学, 2013, 30(9): 234-242. doi: 10.6052/j.issn.1000-4750.2012.05.0382
    [6] 张军锋, 葛耀君, 赵林.  双曲冷却塔结构特性新认识 . 工程力学, 2013, 30(6): 67-76. doi: 10.6052/j.issn.1000-4750.2012.02.0124
    [7] 张军锋, 葛耀君, 赵林.  基于风洞试验的双曲冷却塔静风整体稳定研究 . 工程力学, 2012, 29(5): 68-77.
    [8] 刘智光, 陈健云.  考虑材料组成特性的混凝土轴拉破坏过程细观数值模拟 . 工程力学, 2012, 29(7): 136-146. doi: 10.6052/j.issn.1000-4750.2010.09.0682
    [9] 任青文, 吴 晶.  极限荷载的突变模型研究 . 工程力学, 2011, 28(7): 1-005.
    [10] 孟焕陵.  组合构件双重非线性分析模型研究与应用 . 工程力学, 2010, 27(7): 244-249.
    [11] 高令飞, 郑小平, 姚振汉.  用边界元行列子域法模拟非均匀脆性材料的弹性特性和破坏过程 . 工程力学, 2010, 27(5): 8-013,.
    [12] 熊铁华, 梁枢果, 邹良浩.  风荷载下输电铁塔的失效模式及其极限荷载 . 工程力学, 2009, 26(12): 100-104,.
    [13] 王晓峰, 杨庆山.  空间薄壁截面梁的材料非线性有限元模型 . 工程力学, 2009, 26(8): 138-142.
    [14] 杨菊英, 张娟霞, 王大国, 马天辉.  颗粒增强复合材料中界面破裂过程数值模拟研究 . 工程力学, 2009, 26(1): 18-024.
    [15] 于庆磊, 唐春安, 朱万成, 唐世斌.  基于数字图像的混凝土破坏过程的数值模拟 . 工程力学, 2008, 25(9): 0-078.
    [16] 桂国庆, 王玉娥.  铝合金单层球面网壳的非线性稳定分析 . 工程力学, 2006, 23(增刊Ⅱ): 0-035,.
    [17] 刘福林.  用Mises屈服条件求内边界固支环板的极限荷载 . 工程力学, 2003, 20(1): 162-165.
    [18] 程进, 江见鲸, 肖汝诚, 项海帆.  大跨度钢拱桥结构极限承载力分析 . 工程力学, 2003, 20(2): 7-10.
    [19] 陈永强, 林葱郁, 姚振汉, 郑小平.  非均匀材料破坏过程数值模拟的边界元法研究 . 工程力学, 2003, 20(3): 19-25.
    [20] 刘小强, 吴惠弼.  高层钢框架的非线性分析模型 . 工程力学, 1993, 10(4): 42-51.
  • 加载中
图(16) / 表 (1)
计量
  • 文章访问数:  51
  • HTML全文浏览量:  12
  • PDF下载量:  18
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-05-27
  • 修回日期:  2020-10-26
  • 网络出版日期:  2021-02-03
  • 刊出日期:  2021-03-25

考虑材料非线性的RC双曲冷却塔风致破坏过程

doi: 10.6052/j.issn.1000-4750.2020.05.0324
    基金项目:  国家自然科学基金项目(51508523)
    作者简介:

    裴 昊(1995−),男,河南郑州人,硕士生,主要从事结构抗风研究(E-mail: 1315491065@qq.com)

    朱 冰(1995−),男,江西九江人,硕士生,主要从事结构抗风研究(E-mail: 535453402@qq.com)

    刘庆帅(1996−),男,河南新乡人,硕士生,主要从事结构抗风研究(E-mail: 727428364@qq.com)

    通讯作者: 张军锋(1983−),男,河南平顶山人,副教授,博士,主要从事结构和桥梁抗风研究(E-mail: brilliantshine@163.com)

摘要: 为明确RC双曲冷却塔的极限静风荷载及其风致破坏过程,以一座大型冷却塔为例,采用ABAQUS中分层壳单元模拟塔筒的混凝土和双向钢筋网,在线性计算和非线性计算验证的基础上,采用弥撒开裂模型和双线性模型模拟混凝土和钢筋的非线性受力特征,分析了自重+风荷载(工况Ⅰ)和自重+冬温+风荷载(工况Ⅱ)两种工况下的静风破坏过程,并结合荷载位移曲线、内力、应变和应力的发展和分布探究了其破坏机理。结果表明:对于工况Ⅰ,风荷载系数为λ=1.384时,首先在迎风子午线相对高度hs/Hs=0.37位置因子午向受拉出现环向裂缝且为贯通裂缝,之后整个迎风区的环向裂缝不断增加并沿环向扩展,迎风区内表面和侧风区外表面也因环向弯矩而出现子午向裂缝,开裂位置的钢筋应力迅速增加,荷载-位移曲线也由线性进入非线性并且迎风区和侧风区位移迅速增加,迎风区和侧风区内力也表现出明显的重分布特征,最后因混凝土持续开裂和喉部区域双向钢筋屈服而破坏,极限风荷载系数为λ=2.007;对于工况Ⅱ,温度效应的计入使λ=1.0时即在侧风区上部首先出现子午向裂缝,子午向裂缝开展也较工况Ⅰ略严重,但因温度的弯矩效应有限,其荷载位移曲线和工况Ⅰ基本一致,且仍以迎风区环向开裂最为显著,结构破坏依然由风荷载控制,最终的极限荷载荷载系数为λ=1.842。

English Abstract

张军锋, 裴昊, 朱冰, 刘庆帅. 考虑材料非线性的RC双曲冷却塔风致破坏过程[J]. 工程力学, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
引用本文: 张军锋, 裴昊, 朱冰, 刘庆帅. 考虑材料非线性的RC双曲冷却塔风致破坏过程[J]. 工程力学, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
Jun-feng ZHANG, Hao PEI, Bing ZHU, Qing-shuai LIU. WIND-INDUCED FAILURE PROCESS OF RC HYPERBOLIC COOLING TOWERS CONSIDERING THE MATERIAL NONLINEARITY[J]. Engineering Mechanics, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
Citation: Jun-feng ZHANG, Hao PEI, Bing ZHU, Qing-shuai LIU. WIND-INDUCED FAILURE PROCESS OF RC HYPERBOLIC COOLING TOWERS CONSIDERING THE MATERIAL NONLINEARITY[J]. Engineering Mechanics, 2021, 38(3): 228-238. doi: 10.6052/j.issn.1000-4750.2020.05.0324
  • RC双曲冷却塔作为一种高耸空间薄壁结构,风荷载是其设计控制荷载[1]。渡桥电厂冷却塔风毁事故[2]和不断攀升的结构高度使结构的抗风安全性始终备受关注,既有研究主要关注结构的风荷载特性、群塔干扰效应、风振效应以及弹性稳定问题[3]。实际上,渡桥电厂事故被一些学者认为是风致失稳破坏,从而引发了近四十年的稳定研究热潮[4-8]。然而,Mang 等[9]于1983年采用有限元方法对Port Gibson冷却塔首次进行了考虑材料非线性的极限荷载分析,并认为其风致破坏应属于材料破坏而非弹性失稳。之后,冷却塔的静风极限荷载分析受到关注。既有研究基本均对塔筒采用分层壳单元,沿厚度划分为多个层次并用钢板层模拟钢筋。对于荷载主要考虑自重和风荷载[9-13],部分学者亦考虑了温度效应[14-15]。但由于问题本身的复杂性以及各学者所采用的软件、单元、材料本构、计算方法以及荷载取值等的差异,即使对Port Gibson同一冷却塔的计算结果也表现出较大的离散性[16]:所得开裂荷载系数和极限荷载系数分别分布在0.75~1.58和1.49~3.15,并且其中亦有文献[11]认为是失稳破坏。

    近年来,随着国内冷却塔建设规模的快速发展,国内学者亦注意到对冷却塔应考虑材料非线性进行分析。但初期研究往往比较粗略,未能对混凝土的材料非线性进行准确模拟[17-18]。文献[19]进行了考虑初始缺陷的极限风荷载分析,但未就关键的材料非线性模型和模拟方法进行介绍。另外,对塔筒破坏的原因也有不同的结论:文献[20]关注了动力极限风荷载,但仅对塔筒的受压损伤进行分析而对更关键的受拉损伤却没有说明,并认为迎风区混凝土的破坏是受压破坏;文献[21]模拟了冷却塔的倒塌过程,其结果显示冷却塔破坏源于混凝土受拉开裂和环向钢筋断裂引起的材料强度失效,这均与大多认为冷却塔源于迎风区子午向受拉破坏的结论[9-10, 12]不一致。鉴于塔筒RC材料非线性模拟的困难,文献[22-23]将塔筒视为均质材料并赋予弹塑性材料特性,但这种处理无法准确模拟混凝土开裂后的特性以及钢筋和混凝土之间的协同工作特点,且两者对同一冷却塔所取的弹塑性参数也差异显著。

    以一座大型冷却塔为例,采用ABAQUS软件对其进行极限静风荷载计算。采用分层壳单元模拟塔筒,采用弥撒开裂模型模拟混凝土的非线性受压特性、受拉开裂和拉伸硬化效应,并计入几何非线性。从荷载-位移曲线、变形模式、内力、应变和应力的发展对塔筒的破坏过程进行系统分析,并从中探究其破坏机理。

    • 钢筋混凝土的承载特性表现出高度非线性,合理地反映这些非线性因素是准确获得冷却塔破坏过程的必要前提。本文采用弥散开裂模型模拟混凝土非线性本构,该模型可以准确地模拟混凝土构件受弯时的非线性受力特性[24],也是既有冷却塔极限风荷载分析的主要方式[9-15]

      该模型通过4个方面的材料行为模拟混凝土的非线性特性:非线性受压特性,破坏面,剪切保留和拉伸硬化,详见文献[16]。

      剪切保留模型如式(1)[25]图1[13]所示。

      $$ G'{\rm{ = }}\beta (1 - \varepsilon _{{\rm{cr}}}/\varepsilon _{\rm{m}})G_0 $$ (1)

      式中:G0为闭合混凝土的剪切模量;G'为开裂混凝土的剪切模量;β为裂缝闭合后的剪切刚度系数;εtεcrεm分别为混凝土开裂应变、开裂后垂直于裂缝方向的应变和剪切刚度为0时的应变。参考文献[12-13],下文取β=0.25,εm=0.004。

      拉伸硬化效应[10-12]用以模拟混凝土开裂后的行为[24]。与文献[10-11]一致,下文同样假定混凝土开裂后应力线性卸载(图2),在应变达到极限应变εu=TSεt时应力降至0,TS即为拉伸硬化参数。既有研究[10-12]表明TS对结果影响较大,在板壳结构试验结果验证模型的基础上[16],下文取TS=15。文献[16]还分析了TS对结果的影响。

      图  1  剪切保留模型

      Figure 1.  Shear retention model

      图  2  拉伸硬化模型

      Figure 2.  Tension stiffening model

    • 尽管钢筋被模拟为平面应力钢板层,但其实际是单向受力,故采用单轴双线性应力-应变曲线模拟其非线性本构[9, 12, 15],如图3所示。其中:σyεy分别为屈服应力和屈服应变;EsEp分别为弹性模量和塑性模量,下文取Ep=0.05Es

      图  3  钢筋双线性模型

      Figure 3.  Bilinear model for reinforcement

    • 本文所选冷却塔如图4所示,塔高177 m。下文采用塔筒相对高度hs/Hs表示所关注的高度位置hs,其中Hs为塔筒总高度(图4)。塔筒为渐变厚度(图5(a)),最小厚度0.271 m,塔筒下缘最大厚度1.40 m,上缘厚度渐增为0.40 m,塔筒通过48对Φ1.30 m人字柱与基础连接。塔筒和顶端檐口均采用C40混凝土和HRB400级钢筋,塔筒环向和子午向配筋率如图5(b)所示,下支柱采用C45混凝土。

      图  4  冷却塔特征尺度示意

      Figure 4.  Characteristic dimensions of cooling tower

      图  5  塔筒结构参数

      Figure 5.  Structural parameters of tower shell

      依规范[26]取基本风速V0=31 m/s,场地类别为B类,对应风振系数β=1.9。温度和风荷载联合作用下的环向弯矩亦需重视[27],考虑到冬温和夏温效应接近[28],故仅考虑冬温作用:依规范[26]塔外设计气温取−9.6 ℃,塔内设计气温取(环梁有挡水设施+单元系统)模式。

      由于结构和荷载均对称,为降低计算代价,仅建立半结构模型[9-12, 14]。塔筒采用壳单元S4R模拟,经不同网格密度划分对比将塔筒划分为环向×子午向=96×72个单元,并沿厚度划分9个混凝土层加4个单向受力钢板层(图6)。

      图  6  分层壳示意图

      Figure 6.  Layered shell element

      其中,最外层(也即图6中的第1和第9层)厚30 mm用以模拟混凝土保护层,其它各层按1∶1∶1.5∶2∶1.5∶1∶1确定。不同网格和层次划分的对比分析详见文献[16]。钢板层厚度根据配筋面积相等进行确定,其位置均在相应钢筋的中心位置。下支柱和顶端檐口均采用梁单元B31模拟。

      由于渡桥电厂冷却塔的破坏以及文献[15, 21]的计算结果均显示冷却塔的风致破坏源于塔筒而非下支柱,故仅考虑塔筒和顶端檐口中的材料非线性。材料属性依规范取其标准值[16, 29],由此得混凝土的开裂应变为εt=ft/Ec=7.354×10−5,其中抗拉强度ft=2.39 MPa,弹性模量Ec=32.5 GPa。此外,对于混凝土的非线性受压特性,取初始受压屈服强度为0.4倍抗压强度,即$f_{\rm{c}}' $=0.4fc

      考虑两种荷载工况并采用牛顿迭代法进行计算:自重+风荷载(G+λW)和自重+冬温+风荷载(G+T+λW),其中G、T和W分别为规范[26]所列荷载标准值。对于荷载效应,轴力F以受拉为正,弯矩M以外表面受拉为正;轴力和弯矩均以X方向为环向,Y方向为子午向,如FYMX分别为子午向轴力和环向弯矩,具体可参考文献[30]。径向位移以内凹为正,外凸为负。环向角度θ以迎风点为0°,逆时针转动为正。

      需要说明,在对本算例进行分析前,还首先根据文献[10, 12]对Port Gibson冷却塔进行了重复计算[16]。尽管所用软件和单元不同,材料非线性模型和计算方法等也存在差异,但所得结果基本吻合(图7),从而验证了本文所用软件、单元以及材料非线性模型和计算方法的合理性。

      图  7  Port Gibson塔的荷载位移曲线

      Figure 7.  Load displacement curves of Port Gibson tower

    • 作为非线性分析的基础,首先应了解冷却塔在各荷载作用下的内力特征。自重作用下(图8)塔筒以双向受压为主,且子午向压力远大于环向压力,双向弯矩可以忽略。自重产生的轴压力将有助于抵抗风荷载产生的拉力。冬温荷载作用下(图9)塔筒以双向正弯矩为主,双向轴力可以忽略。

      图  8  自重作用下塔筒内力分布

      Figure 8.  Internal forces of shell under gravity

      图  9  冬温荷载作用下塔筒内力分布

      Figure 9.  Internal forces of tower shell under gravity

      风荷载作用下塔筒内力分布较为复杂(图10),除MY外,各荷载效应均较显著,并以子午向轴力最为显著。对各内力特征的详细介绍可见文献[27]。

      图  10  风荷载作用下塔筒内力分布

      Figure 10.  Internal forces of tower shell under wind load

    • 在两种工况下,施加风荷载前,结构变形可以忽略。施加风荷载后,各截面的变形规律基本一致,并且考虑到对称性仅给出0°≤θ≤180°范围的结果(图11):即迎风区和侧风区分别产生显著的内凹和外凸位移,且前者比后者增加更为显著,背风区的变形始终较小。在进入非线性阶段后,迎风区内凹位移最大值均始终在θ=0°位置,但随风荷载的增加而有上下移动;而侧风区外凸位移最大点则同时有环向和子午向移动:图11所示从下到上的4个断面的外凸位移最大点逐渐由θ=40°区域后移至θ=75°区域,而在子午向则逐渐由hs/Hs=0.65区域上移至塔筒顶部。

      图  11  塔筒横向截面变形

      Figure 11.  Deformed shapes of tower shell: sliced views

    • 为直观体现结构的非线性效应,选取Ⅰ点(θ=0°、hs/Hs=0.782,迎风点喉部)和Ⅱ点(θ=60°、hs/Hs=0.70)这两个内凹和外凸位移始终较为显著的位置,以及Ⅲ点(θ=120°、hs/Hs=0.70)作为对照位置,给出荷载-位移曲线如图12所示。可以看出,两种荷载工况下,冷却塔破坏过程中的非线性主要体现在迎风区和侧风区,其他区域的非线性极不明显。

      图  12  两种荷载工况所得荷载位-移曲线

      Figure 12.  Load-displacement curves of two load cases

      对于工况G+λW,以Ⅰ点为例,其荷载位移曲线可划分为3个阶段。第1阶段(OA段,λ<1.384)为线性响应阶段:位移随风荷载线性增加,整个塔筒未出现裂缝。第2阶段(AB段,1.384≤λ<1.876)为裂缝扩展阶段:由于风荷载在塔筒迎风区产生的子午向轴拉力逐渐达到并超过自重作用下的子午向轴压力,故λ=1.384时,迎风子午线上hs/Hs=0.37位置混凝土应力首先达到峰值并进入卸载阶段,也即该处产生裂缝,且该裂缝为贯穿塔筒壁厚的环向裂缝;随着风荷载的增加,环向裂缝在塔筒迎风区不断形成并扩展,塔筒局部刚度逐渐下降,结构响应进入明显的非线性阶段,迎风区和侧风区位移迅速增大,内力也出现重分布,详见下节。第3阶段(BC段,1.876≤λ≤2.007)为硬化阶段:λ=1.876时,迎风子午线上hs/Hs=0.70位置内侧子午向钢筋首先受拉屈服,结构进入硬化阶段,风荷载又有进一步增加,钢筋屈服区域沿环向逐渐扩展至θ=16°,沿子午向则扩展至hs/Hs=0.52~0.80。并且λ=1.529时,迎风子午线hs/Hs=0.62高度内侧和θ=75°、hs/Hs=0.95高度外侧出现子午向裂缝。最终,由于迎风区混凝土的持续开裂和钢筋屈服,λ=2.007时计算终止(表1)。这一荷载位移曲线形式与文献[9-15]所得结果较为类似。

      表 1  两种荷载工况结果

      Table 1.  Results of two load cases

      特征结果G+λWG+T+λW计入温度效应后的结果偏差/(%)
      开裂荷载系数1.384(A)1.000(A')−28
      屈服荷载系数1.876(B)钢筋未屈服
      极限荷载系数2.007(C)1.842(C')−8
      极限位移0.836 m0.531 m−26
      注:括号中的字母如图12中所示。

      对于工况G+T+λW,由于冬温荷载产生双向正弯矩,使塔筒外表面混凝土在风荷载作用前即处于双向受拉状态,而风荷载作用后侧风区的环向正弯矩效应进一步增加,迎风区子午向也将处于拉弯联合受力状态,所以塔筒提前开裂:λ=1.000时,塔筒侧风区θ=73°、hs/Hs=0.93位置外表面即首先出现子午向裂缝;λ=1.080时,塔筒迎风子午线上hs/Hs=0.08位置外表面开始产生环向裂缝。因此,该工况的荷载-位移曲线略提前进入非线性阶段,且整体刚度偏小。

      对比表1可知,冬温荷载显著减小了冷却塔的开裂荷载和极限位移,但对极限荷载影响不大。这是由于尽管冬温荷载产生的双向正弯矩导致塔筒提前开裂,但此弯矩效应较小且固定不变,且随着风荷载的持续增加,其与风荷载产生的荷载效应相比也越来越小,故冷却塔的最终破坏依然由风荷载控制。也正因此,两种荷载工况所得荷载-位移曲线几乎一致;二者变形和裂缝开展情况也相似(见4.2节和4.3节),只是在同一λ下,计入冬温效应的变形和裂缝开展更为严重。因此,对于本例,温度作用只是增加了环向和子午向弯曲损伤,这在文献[14-15]中亦有体现。

    • 荷载-位移曲线表现出的明显非线性特征直接源于混凝土开裂导致的结构刚度损伤,这也使结构刚度分布发生变化,进而导致内力重分布并加剧结构开裂。由于两种荷载工况所得荷载-位移曲线变形模式相差不大,为观察塔筒破坏过程中的内力重分布情况,以工况G+λW为例,给出hs/Hs=0.37和0.70这两个断面的内力分布如图13所示。

      图  13  塔筒内力分布

      Figure 13.  Distribution of internal forces of tower shell

      图13可知,在λ≤1.384时,整个塔筒所有内力均呈线性变化,并且由于计入了自重效应,FY在迎风区逐渐由轴压力表现为显著的轴拉力,在侧风区始终为显著的轴压力,在背风区亦为轴压力但幅值较小;自重效应的存在也使FX在整个环向均为轴压力;MX在整个环向幅值较小,在迎风区和侧风区分别为负弯矩和正弯矩。

      λ=1.384时,由于hs/Hs=0.37高度首先出现环向裂缝(见下文4.5节),θ=0~15°的各内力均开始表现出非线性并逐渐扩展到θ=45°位置,即出现内力重分布。并且环向裂缝的出现不仅影响子午向受力,也影响环向受力:与此同时,hs/Hs=0.70高度θ=0°的FXMX也开始出现内力重分布,而该高度FY依然保持线性变化,直至λ=1.48时开始非线性变化。而此时此处尚未开裂,这同样是内力重分布的结果。此后,随着塔筒开裂的继续发展,塔筒内力重分布显著区域也不断扩展至θ=90°位置并在两个高度表现出相似的规律。

      对于子午向的内力重分布,θ=0°位置开裂后,此处的FY略有下降,但其两侧的FY快速增加,用以补偿θ=0°子午线上FY的下降。迎风区大范围受拉开裂后,整个迎风区的FY不再线性增加,还造成侧风区的FY亦不再随荷载线性增加。与环向裂缝主要因子午向轴拉力引起不同,子午向裂缝主要由环向弯矩引起,内力重分布也主要表现为迎风区和侧风区MXFX幅值和环向波动的急剧增加。

    • 图14给出了两种荷载工况所得塔筒最外层和最内层混凝土最大主应变εpr分布。由于冬温和风荷载的弯矩效应有限,两种工况下,塔筒均主要依靠薄膜作用抵抗风荷载,所以塔筒内外层混凝土εpr随风荷载的变化规律基本相同,尤其是工况G+λW。对于工况G+λW,λ=1.50时,迎风区中下部已有较大面积混凝土产生了环向裂缝,故结构开始表现出明显的非线性特征。λ=1.80时,迎风区已产生大面积开裂,并存在4个明显的贯穿壁厚的环向裂缝带(亦可见图15(c));由于迎风区和侧风区分别存在较大的环向负弯矩和正弯矩,故上述区域内侧和外侧混凝土均产生了子午向裂缝,但其宽度相对较小,直到结构破坏也未贯穿壁厚。λ=2.00时,迎风区hs/Hs=0.08~0.80区域全部开裂,且开裂区域与λ=1.90时近乎一致,仅裂缝宽度急速增大。最终,冷却塔因大面积的混凝土开裂和钢筋屈服而破坏。计入冬温荷载后,冬温荷载产生的双向正弯矩导致塔筒提前开裂,且同一λ下裂缝开展情况更为严重,主要表现在侧风区外侧。

      图  14  塔筒混凝土最大主应变(仅给出εt=7.354×10−5~5.2×10−3范围内的拉应变)

      Figure 14.  Maximum principal strain of concrete in tower shell

      图  15  塔筒外侧混凝土和钢筋应力分布

      Figure 15.  Distribution of stress of exterior concrete and steel in tower shell

    • 由于迎风子午线上的子午向受拉开裂和hs/Hs=0.70断面的环向拉弯开裂最为显著,为进一步了解塔筒破坏过程中的应力分布情况和破坏过程,给出G+λW工况下迎风子午线上最外层混凝土和外侧钢筋的子午向应力分布和hs/Hs=0.70断面表层混凝土和钢筋的环向应力分布,同时给出上述两个位置处的应力荷载曲线,如图15所示。

      图15(a)可以看出,λ=0.50时,除塔筒顶部和底部边界外,塔筒整个高度范围均处于双向受压状态。λ=1.0时,风荷载作用下的子午向轴拉力已经超过自重的压力效应,迎风区大部分区域已进入子午向受拉状态。λ=1.50时,hs/Hs=0.10~0.50范围内混凝土子午向应力已达到抗拉强度,且hs/Hs=0.37处混凝土已进入卸载阶段。λ=1.80时,迎风区已有4个位置混凝土应力完全降至0,即迎风区产生了4个完全开裂的环向裂缝带,裂缝区荷载完全由钢筋承担导致钢筋应力激增(图15(b))。λ=2.00时,迎风区子午向开裂区域几乎不再随λ增加,hs/Hs=0.08~0.80范围混凝土完全退出工作,hs/Hs=0.52~0.80区域钢筋屈服(图15(b)图15(c))。

      另外,从图15(c)更可直观看出,混凝土开裂后其应力急剧减小至0,与此同时裂缝区钢筋的应力则急剧增加。4个主裂缝带产生的位置顺序依次为hs/Hs=0.37、0.52、0.17和0.70,对应的λ依次为1.384、1.451、1.525和1.631。尽管喉部位置混凝土最后开裂,但由于该处壁厚和配筋均较小,故该处钢筋首先屈服。这也说明,子午向受拉开裂和钢筋屈服是冷却塔在风荷载作用下破坏的直接原因。

      对于环向应力,由图15(d)可以看出,hs/Hs=0.70高度外侧混凝土环向应力的分布与图13所示MX的分布规律一致,这也表明冷却塔的环向应力状态受弯矩控制。λ≤1.50时,除迎风子午线局部区域,环向应力基本处于线性状态;进入非线性状态后,受环向正弯矩影响,θ=25°~65°区域外侧混凝土应力达到抗拉强度,但此后其应力并未快速卸载,而是维持在一个较高的水平,这是由于该区域的裂缝宽度在塔筒破坏过程中始终较小(图14),这也体现在该区域钢筋应力在混凝土开裂后增长不大(图15(e))。而由图15(f)可知,λ=1.529和1.696时,hs/Hs=0.70高度θ=0°内侧和θ=50°外侧混凝土环向应力先后达到抗拉强度并进入卸载阶段,这两处混凝土产生子午向裂缝,且钢筋的拉应力迅速增加。尽管这两处环向钢筋均未屈服,但λ=1.98时迎风子午线上hs/Hs=0.57处内侧环向钢筋屈服,并迅速扩展至hs/Hs=0.53~0.61范围。这也表明,尽管结构破坏源于子午向受拉破坏,但环向受力亦需关注。

    • 以一座大型RC冷却塔为例,通过数值模拟研究了冷却塔在两种荷载工况下的非线性静风响应,从荷载位移曲线、变形模式、内力、应变和应力分布对其破坏过程进行了系统阐述,并从中探究了其破坏机理。主要结论如下:

      (1)通过对Port Gibson冷却塔的计算和与既有结果对比,表明ABAQUS的分层壳单元配合弥散开裂模型可以有效进行冷却塔极限静风荷载分析。

      (2)在自重+风荷载作用下,在λ=1.384时,迎风点hs/Hs=0.37位置在子午向轴拉力作用下首先开裂并产生贯穿厚度的环向裂缝,之后裂缝随风荷载增加而不断形成并沿环向扩展至整个迎风区,沿子午向几乎扩展至整个高度,开裂位置的钢筋应力迅速增加;λ=1.529时,迎风子午线hs/Hs=0.62高度内侧和θ=75°、hs/Hs=0.95高度外侧混凝土因环向弯矩而开裂;λ=1.876时,hs/Hs=0.70处子午向钢筋首先受拉屈服,之后钢筋屈服区域沿环向逐渐扩展至θ=16°,沿子午向则扩展至hs/Hs=0.52~0.80;λ=1.98时迎风子午线上hs/Hs=0.57处内侧环向钢筋屈服,并迅速扩展至hs/Hs=0.53~0.61范围。λ=2.007时,因混凝土持续开裂和钢筋屈服,冷却塔达到抗拉极限状态而破坏。

      (3)在自重+风荷载作用下,塔筒开裂之前各响应均表现为线性特征,塔筒开裂之后,塔筒局部刚度下降,荷载-位移曲线也由线性进入非线性并且迎风区和侧风区位移迅速增加,同时迎风区和侧风区内力也表现出明显的重分布特征。结构破坏时,迎风区子午向受拉导致大面积的受拉贯通环向裂缝和钢筋屈服,而环向仅有塔筒中上部受弯开裂且裂缝宽度有限,钢筋屈服范围也较小,因此可认为迎风子午向受拉导致的塔筒开裂和钢筋屈服是结构破坏的首要原因,但环向受力亦需关注。

      (4)在自重+冬温荷载+风荷载作用下,由于冬温荷载产生的双向正弯矩效应增加了塔筒的弯曲损伤,导致塔筒侧风区在λ=1.0时提前产生子午向裂缝。但由于此弯矩效应固定不变且远小于风荷载作用下的子午向受拉效应,故结构的破坏依然由风荷载控制,极限荷载系数为λ=1.842。

      (5)两种荷载工况下结构均在λ=1.50时明显进入非线性段,且之后塔筒变形急速增加,此荷载系数这也与结构设计时的风荷载分项系数1.4相当,说明在设计风荷载下结构的冗余安全度较为有限。

参考文献 (30)

目录

    /

    返回文章
    返回