岩土力学 ›› 2020, Vol. 41 ›› Issue (7): 2485-2493.doi: 10.16285/j.rsm.2019.1343

• 数值分析 • 上一篇    下一篇

静水压力下内摩擦角对圆形洞室围岩局部破裂 影响的连续-非连续模拟

白雪元1，王学滨1, 2，舒芹1   

  1. 1. 辽宁工程技术大学 力学与工程学院，辽宁 阜新 123000；2. 辽宁工程技术大学 计算力学研究所，辽宁 阜新 123000
  • 收稿日期:2019-08-02 修回日期:2020-01-06 出版日期:2020-07-10 发布日期:2020-09-20
  • 通讯作者: 王学滨，男，1975年生，博士，教授，主要从事工程材料变形、破坏及稳定性方面的研究工作。E-mail: wxbbb@263.net E-mail: xueyuan_bai@126.com
  • 作者简介:白雪元，男，1991年生，博士研究生，主要从事计算固体力学方面的研究工作。
  • 基金资助:
    国家自然科学基金（No. 51874162）。

Continuum-discontinuum simulation of effects of internal friction angle on local fracture of circular cavern surrounding rock under hydrostatic pressure

BAI Xue-yuan1, WANG Xue-bin1, 2, SHU Qin1   

  1. 1. College of Mechanics and Engineering, Liaoning Technical University, Fuxin, Liaoning 123000, China; 2. Institute of Computational Mechanics, Liaoning Technical University, Fuxin, Liaoning 123000, China
  • Received:2019-08-02 Revised:2020-01-06 Online:2020-07-10 Published:2020-09-20
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (51874162).

PDF (Chinese)

164

PDF (English)

7

摘要： 为了更准确地模拟岩石由连续向非连续转化，降低网格依赖性，在自主开发的拉格朗日元与离散元耦合的连续?非连续方法的基础上，发展了一种内部开裂方法。模型被离散成精度较高的四边形单元，在开裂后，裂缝可沿四边形单元对角线和单元边界扩展。模拟了巴西圆盘岩样、单轴压缩岩样和静水压力下洞室围岩的变形?开裂过程。结果表明：对于巴西圆盘岩样，拉裂缝由圆盘中心向上下加载板扩展，直至贯穿岩样，拉裂缝较为平滑；对于单轴压缩岩样，剪裂缝较为平直，贯穿岩样的主剪裂缝方向与摩尔?库仑剪切面方向较为一致。由静水压力和内摩擦角对洞室围岩变形?开裂过程的影响规律可以发现：静水压力下，V形坑形成于洞室表面附近；剪裂缝向围岩内部扩展，形成狭长剪裂缝，其分布类似于剪切滑移线。静水压力越大，剪裂缝的范围越大。内摩擦角越大，剪裂缝的范围越小，狭长剪裂缝的扩展方向与环向的夹角越小。由天生桥二级水电站某洞室围岩的模拟可以发现，围岩呈现4个V形坑，这与现场观测结果基本吻合。

关键词: 网格依赖性, 洞室围岩, 变形?开裂过程, 静水压力, 内摩擦角, 剪切滑移线

Abstract: In order to simulate the transition of rock from continua to discontinua more accurately and reduce the mesh dependence, an element splitting method is developed based on the continuum-discontinuum method where the Lagrangian element method is coupled with discrete element method. The computational domain is discretized into high-precision quadrilateral elements; after cracking, the cracks can propagate along the diagonal lines of quadrilateral elements and the edges of triangular or quadrilateral elements. Deformation-cracking processes of the Brazilian disc rock specimen, uniaxial compression rock specimen and surrounding rocks under hydrostatic pressure are simulated. The following results are found. (1) For the Brazilian disc rock specimen, tensile cracks propagate from the center of the disc to the top and bottom until they penetrate through the disc, and the tensile cracks are relatively smooth; for the uniaxial compression rock specimen, shear cracks are relatively straight, and the direction of the main shear crack penetrating through the rock specimen is consistent with the shear plane direction obtained according to the Mohr-Coulomb criterion. (2) From the effects of hydrostatic pressure and internal friction angle on deformation-cracking processes, the following results are found. Under hydrostatic pressure, firstly, V-shaped notches appear near the cavern surface due to shear crack propagation; then, long and curved shear cracks appear due to further shear crack propagation whose distribution is similar to shear slip-lines. With an increase of hydrostatic pressure, the ranges of shear cracks increase, whereas with an increase of internal friction angle, the ranges of shear cracks decrease, and the angle between the propagation direction of the long shear crack and the annular direction decreases. (3) According to the simulation of the surrounding rock in phase-Ⅱ Tianshengqiao hydropower station, four V-shaped notches are observed, which is in a good agreement with field observation.

Key words: mesh dependence, cavern surrounding rock, deformation-cracking process, hydrostatic pressure, internal friction angle, shear slip-line

中图分类号: 

  • TU 456
[1] 姜景山, 左永振, 程展林, 潘家军, . 不同密度粗粒料强度特性的大型真三轴试验[J]. 岩土力学, 2020, 41(8): 2601-2608.
[2] 杨福见, 胡大伟, 田振保, 周辉, 卢景景, 罗宇杰, 桂树强, . 高静水压力压实作用下疏松砂岩渗透 特性演化及其机制[J]. 岩土力学, 2020, 41(1): 67-77.
[3] 孙德安，张乾越，张 龙，朱赞成，. 高庙子膨润土强度时效性试验研究[J]. , 2018, 39(4): 1191-1196.
[4] 田 密, 张 帆, 李丽华, . 间接测量数据条件下岩土参数空间变异性 定量分析方法对比研究[J]. 岩土力学, 2018, 39(12): 4673-4680.
[5] 王国富，曹正龙，路林海，王 蓉，王 丹，韩 帅，. 黄河冲积层静止土压力系数原位测定与分析[J]. , 2018, 39(10): 3900-3906.
[6] 田 密，李典庆，曹子君，方国光，王 宇,. 基于贝叶斯理论的土性参数空间变异性量化方法[J]. , 2017, 38(11): 3355-3362.
[7] 王 艳，傅俊醒，唐 强，胡安详，刘佳鑫，. 表面活性剂改性粉土抗剪强度特性的研究[J]. , 2016, 37(S1): 329-333.
[8] 王仕传，孙本杰，邵 艳. 改进的主动土压力计算方法[J]. , 2015, 36(5): 1375-1379.
[9] 夏开宗 ，刘秀敏 ，陈从新 ，宋娅芬 ，欧 哲 ，龙 毅，. 考虑突变理论的顺层岩质边坡失稳研究[J]. , 2015, 36(2): 477-486.
[10] 周 辉 ，杨艳霜 ，刘海涛，. 岩石强度时效性演化模型[J]. , 2014, 35(6): 1521-1527.
[11] 李志平 ，彭振斌 ，肖尊群 ，张岳安 . 含水率影响下尾矿砂内摩擦角随标贯击数变化研究[J]. , 2013, 34(5): 1340-1344.
[12] 李丽华，肖衡林，唐辉明 ，马 强 ，陈 辉 ，孙 龙 . 轮胎碎片-砂混合土抗剪性能优化试验研究[J]. , 2013, 34(4): 1063-1067.
[13] 曹宇春，杨建辉. 基于有效固结应力法确定结构性黏性土不排水抗剪强度[J]. , 2013, 34(11): 3085-3090.
[14] 徐日庆 ，廖 斌 ，吴 渐 ，畅 帅 . 黏性土的非极限主动土压力计算方法研究[J]. , 2013, 34(1): 148-154.
[15] 蔡 建. 原状土的抗剪强度研究[J]. , 2012, 33(7): 1965-1971.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
[1] 姚仰平，侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 张力霆，齐清兰，魏静，霍倩，周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[4] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[5] 张明义，刘俊伟，于秀霞. 饱和软黏土地基静压管桩承载力时间效应试验研究[J]. , 2009, 30(10): 3005 -3008 .
[6] 刘 斌，李术才，李树忱，钟世航. 隧道含水构造直流电阻率法超前探测研究[J]. , 2009, 30(10): 3093 -3101 .
[7] 吴 亮，钟冬望，卢文波. 空气间隔装药爆炸冲击荷载作用下混凝土损伤分析[J]. , 2009, 30(10): 3109 -3114 .
[8] 吴世余，余金煌. 不透水地基上堤坝上游坡的附加渗径[J]. , 2009, 30(10): 3151 -3153 .
[9] 周晓杰，介玉新，李广信1. 基于渗流和管流耦合的管涌数值模拟[J]. , 2009, 30(10): 3154 -3158 .
[10] 吴昌瑜，张 伟，李思慎，朱国胜. 减压井机械淤堵机制与防治方法试验研究[J]. , 2009, 30(10): 3181 -3187 .
版权所有 © 《岩土力学》编辑部
地址：湖北武汉武昌小洪山中国科学院武汉岩土力学研究所(430071)　邮编：200042　电话：027-87198484, 87199252, 87199253, 87198752　E-mail: ytlx@whrsm.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发