基于FLAC^3D的隧洞衬砌与围岩接触性状分析.pdfVIP

第32卷第8期 水 电 能 源 科 学 Vo1．32No．8 2014年 8月 W aterResourcesandPower Aug．2014 文章编号 ：1000～7709(2014)08—0112—04 基于 FLAC3D的隧洞衬砌与围岩接触性状分析 张 凌 (中水电海外投资有限公司，北京 100048) 摘要：为了更好地了解围岩与衬砌之间的相互作用，利用三维拉格朗日有限差分法，在商业软件FLAC 中通 过设定接触面来研究不 同工况下围岩与衬砌 间的接触状态及 围岩与衬砌 的应力、应变规律 。通过数值模拟破 坏性试验揭示高应力下待开挖区域变形规律，通过接触面来反映隧洞围岩和衬砌嵌入形态。模拟结果表明， 岩性和开挖支护方式选择在实际地下隧洞施工中起着重要作用。 关键词 ：隧洞 ；衬砌 ；接 触性状 ；FLAC ；数值模 拟 中图分类号 ：TU311 文献标志码 ：A 1 引言 2 数值计算模型的建立 水工隧洞围岩的稳定性和支护结构的可靠性 2．1 计算模型 是实现工程安全、经济、顺利建设 的重要保障，然 使用高地应力下的隧洞模型 (图1)来反映接 而对深埋引水隧洞特别是高地应力区域 的隧洞 ， 触面产生的作用 ，对于衬砌使用 Brick的实体单 无论是设计计算理论还是施工建设技术 目前 尚处 元进行模拟，通过改变岩体及衬砌的力学参数来 于研究探索阶段 。地下岩体开挖前，岩体处于平 反映围岩与衬砌之间的作用。建立 3OmX30m× 衡状态 ，隧洞开挖时破坏 了原应力平衡状态 ，使岩 10m三维模拟模型，共产生 8640个实体单元 ， 体内部的应力重新分布 。从力学的观点来看，地 9515个节点，550个接触面单元。围岩物理力学 下洞室的围岩失稳是由于围岩的应力水平达到或 参数为：弹性模量 89．5～8950．0MPa，泊松 比 超过岩体 的强度范围，形成一个连续贯通的塑性 0．26，粘聚力 1．6MPa，内摩擦角 5O。，侧压力系数 区或滑动面 ，产生较大位移最终导致失稳的问题 。 K1．4，重度 1900kN／m。，饱和抗压强度 50MPa， 隧洞工程的围岩承受着相对较大的岩体 自重和初 剪切模量 61．6～6160．0MPa，埋深 710m。衬 始地应力 的作用 ，若 围岩压力较大 ，容易出现塌 砌结构参数为：弹性模量 l1．7×10。MPa，剪切模 落 、岩爆等现象_l】]，因此隧}同支护在整个隧洞建设 量 5．4×10。MPa，泊松 比0．167，重度 2500kN／m3， 过程中尤为重要。FLAC∞是一种三维显函数有 衬砌厚度r。0．5m。 限差分程序 ，其算法的基础是三维快速拉格朗 日 算法，将计算区域划分为若干四面体单元，每个单 元在给定的边界条件下遵循指定 的线性或非线性 本构关系，若单元应力使得材料屈服或产生塑性 流动，则单元网格可以随着材料的变形而变形，特 别适用于岩石力学的非线性大变形或不稳定(如 滑动或分离)问题，其计算模型网格能以大应变模 图 1 隧洞网格 图 式变 形 并 随 材 料 流 动。鉴 此 ，本 文 通 过 在 Fig．1 Tunnelgrid FLAC如中设定接触面来分析围岩与衬砌的相互 2．2 衬砌与围岩 间接触面的模拟分析理论 作用及作用状态，为隧洞工程支护设计、施工方式 地下洞室开挖衬砌后 ，围岩变形受到衬砌抑 提供参考 。