不同尺寸的圆形隧洞剪切应变局部化过程模拟.pdfVIP

第 20卷 第 4期 中国地质灾害与防治学报 V0】．2O No．4 2009年 12月 TheChineseJournalofGeologicalHazardandControl Dec．2009 不同尺寸的圆形隧洞剪切应变局部化过程模拟 王学滨，潘一山，陶 帅 (辽宁工程技术大学力学与工程学 院，辽宁 阜新 123000) 摘要 ：利用 FLAC模拟 了两个不同直径 圆形隧洞的剪切应变局部化过程。为 了模拟隧洞开挖 ，利用编写的FISH 函数 删除隧洞 内部的单元。岩石服从莫尔库仑剪破坏与拉破坏复合 的破坏准则，破坏之后呈现应变软化 ～理想塑性行 为。本文的模拟分为 3步：首先，将静水压力施加在模型上 ，直至达到静力平衡状态 ；然后 ，利用编写的FISH函数 ，开 挖隧洞；最后 ，计算重新开始 ，直至达到静力平衡状态 (对于小孔隧洞)或者塑性流动状态 (对于大孔 隧洞)。模拟结果 表明，多个 “狗耳”形或 v形坑在小孔隧洞周边 附近产 生，最终，围岩处 于平衡状态。这 一结果与陆家佑和王 昌明 (1994)的实验结果及许多现场观察结果一致。对于大孔隧洞，由于在围岩中出现了多条剪切带，因而隧洞的整个断面 均遭到 了破坏 。这一现象-9现场观察到的猛烈破坏现象类似。隧洞 的剪切应变局部化受隧洞尺寸的影响。 关键词 ：圆形隧洞；剪切应变局部化；剪切带；V形坑 ；静水压力；直径 ；尺寸效应 文章编号 ：1003—8035(2OO9)O4—0101—08 中图分类号 ：P554 文献标识码 ：A 供的圆形巷道网格 (辐射状网格，即离巷道表面越近 ， 1 引言 网格越密)不能较好地模拟出巷道围岩的局部破坏现 对于圆形巷道 ，在模型实验及现场观测中经常会 象 ，仅能模拟出均匀的破坏现象 (圆环形的塑性区)或 发现 ，巷道发生了局部剪切破坏 ，即在 围岩中产生一 者模拟出剪切带离巷道表面越远越粗的不准确现象。 条或多条剪切带或剪切裂纹 ，或者形成 1个或多个 关于剪切带宽度的实验结果表明，其宽度主要依赖于 “狗耳”形 的 v形坑 。对于巷道 的非均匀 (或变 材料属性 ，对于同一种材料 ，其宽度大致为一个定值。 形局部化)破坏现象，开展理论方面的研究是非常困 为此 ，本文作者将尝试另外的处理方法。 难 的。理论研究方面的这一缺陷可以在一定程度上 本文的工作是采用 FLAC模拟不 同直径 圆形隧 被数值试验避免。在过去 ，针对矩形试样 的应变局部 洞的剪切应变局部化过程。先将模型剖分成大小相 化的数值分析较为多见，而关于圆形巷道 的应变局部 等的矩形单元 ，加载至平衡 ，然后再通过 FLAC的内 化 的数值分析非常少见。 嵌编程语言编写 FISH函数模拟开挖。研究结果表 材料的变形局部化与尺寸效应密切相关。就岩 明，这种处理方法可以得到比较符合实际的结果。在 石力学与工程而言，尺寸效应或尺度律的重要性是众 岩石破坏之后 ，本文采用了峰后应变软化一塑性的本 所周知的。Ba~ant和 Chen在其评述文章 的开头就指 构关系 ，而不是峰后脆一塑性的本构关系。后者 出，尺度律是一切物理理论 中最重要的方面…J。如 在应力陡然降低过程中没有塑性应变发生 ，并不完全 果不清楚尺度律，那么物理理论本身也是难于理解 与岩石的变形实质相符。 的 。因此 ，在很多物理和工程 问题中，尺度律 问题 2 计算模型、本构关系、计算步骤及 占据其中心位置就不足为奇了n 。最近 ，文献[12—