高地应力下地下洞室群开挖过程的数值模拟.pdfVIP

第26卷第4期 华 中 科 技 大 学 学 报 (城市科学版) Vol_26No．4 2009年 12月 J．ofHUST．(UrbanScienceEdition) Dec．2009 高地应力下地下洞室群开挖过程的数值模拟 陈景涛 ，朱进明，苏国韶。 (1．武汉理工大学 理学院，湖北 武汉 430070；2．中电武汉电子广场，湖北 武汉 430022； 3．广西大学 土木建筑工程学院，广西 南宁 530004) 摘 要：采用适用于高地应力下硬岩的弹一脆一塑性本构模型，对某水电站地下洞室群的开挖过程进行了数值模 拟 ，得到了高地应力下地下洞室群开挖过程中位移场 、应力场和塑性区的演化过程。比较发现，数值模拟结果 和实测位移及松动圈较为～致 ，表明数值计算的方法和结果正确可信，为高地应力条件下硬岩地下工程的计 算、设计、优化和围岩稳定性分析提供了科学依据。 关键词 ：高地应力； 地下洞室； 开挖过程 ； 数值模拟 中图分类号 ：TU452 文献标识码 ：A 文章编号：1672—7037(2009)40-0005-05 随着人类活动 向地下空间的延伸 ，在高地应 房洞室群为工程实例，采用适用于高地应力下硬 力区修建的地下工程越来越多。例如，在建中的 岩的弹一脆 一塑性本构模型 加。 ，以获得较为符 拉西瓦水电站地下厂房洞群最大地应力为29．6 合实际的破坏的范围、深度、应力和变形，为高地 MPa；设计中的锦屏一级水电站地下厂房区最大 应力下硬岩地下工程的分析与计算奠定基础。 主应力值达37．5lPa；锦屏二级水 电站的引水隧 洞群在长探洞 中实测最大主地应力值高达 42 1 数值计算模型 MPa；规划中的南水北调西线工程部分洞段开挖 时可能会遇到 50MPa左右的高水平挤压应力。 1．1 工程概况 根据相关研究 ¨ ，上述工程都属于高地应力下 地下洞室群主要 由主厂房、主变室、调压室、 硬岩地下工程。高地应力下大型地下工程的施工 尾闸室、尾水管、母线洞等构成，其中主厂房开挖 建设 ，往往需要投入大量的人力、物力和财力。如 尺寸(长 X宽 ×高)为311．75mx30．0mX73．84 果洞室群布置、施工顺序和支护措施合理，不仅可 m，主变室开挖尺寸 232．6m×29．0m×53．0m， 以保证地下工程在施工过程的安全性和使用期问 地下洞室群布置见图 1。地下洞室群位于微风 的长期稳定性，还能节省大量投资。开展高地应 化一新鲜块状结构花 岗岩体 内，主要 由Ⅱ类围岩组 力条件下大型地下工程的数值模拟研究，对确保 成，各类围岩的力学参数见表1，结构面分布见图2。 工程顺利进行、有效合理地开发地下空间具有十 分重要的意义，其成果可为我国水电工程、矿山工 程、能源储存工程等的安全建设与运营提供科学 依据 ，具备 良好的社会经济效益与应用前景 J。 传统的本构模型在模拟高地应力下硬岩破坏 的范围和深度方面不理想。比如采用传统的弹塑 性本构模型进行高地应力下水电站地下厂房洞群 的稳定性分析，由于高地应力地区地应力量值较 大，岩体弹模高，计算结果中常发现，洞室群围岩 出现贯穿性的大面积屈服，这与电站长期安全运 行的实际情况不符，还需耗用大量资金进行支护 图1 地下洞室群布置 处理 。本文以高地应力区的某水电站地下厂 收稿 日期 ：2009-04—16 作者简介：陈景涛 (1976一)，男，河南人，副教授，博士，研究方向为工程材料的力学行为，jtehen2006@163，eom。 基金项 目：国家 自然科学基金 。 · 6 · 华 中 科 技 大 学 学 报 (城市科学版)