›› 2011, Vol. 32 ›› Issue (5): 1575-1582.

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

基于应变能密度理论的岩石破裂数值模拟方法研究

孙 倩,李树忱,冯现大,李文婷,袁 超   

  1. 山东大学 岩土与结构工程研究中心,济南 250061
  • 收稿日期:2010-06-07 出版日期:2011-05-10 发布日期:2011-09-23
  • 通讯作者: 李树忱,教授,博士生导师。E-mail: shuchenli@sdu.edu.cn E-mail:ningyezi1227@yahoo.com.cn
  • 作者简介:孙倩,女,1985年生,硕士研究生,主要从事岩石破坏机理以及数值计算等方面研究
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2010CB732002);山东省自然科学基金资助项目(No. 2009ZRB02285);山东大学自由探索(No. 2010TS038)

Study of numerical simulation method of rock fracture based on strain energy density theory

SUN Qian, LI Shu-chen, FENG Xian-da, LI Wen-ting, YUAN Chao   

  1. Research Center of Geotechnical & Structural Engineering, Shandong University, Ji?nan 250061, China
  • Received:2010-06-07 Online:2011-05-10 Published:2011-09-23

摘要: 采用双线性应变软化本构模型结合能量耗散原理建立了损伤本构方程,并通过应变能密度理论建立了细观单元岩石破坏的能量判别准则。当某一单元所存储的应变能超过某固定值时,单元进入损伤状态,同时单元的损伤程度随着能量耗散的增加而增加,损伤单元的材料属性也随之改变,直到变为具有一定残余强度的单元。随着荷载增加,单元损伤的程度变大,当单元储存的应变能超过所建立的能量判别准则时,定义单元破坏,随着破坏单元的数目不断增多,破坏单元相互连通形成宏观裂纹,实现了利用线性计算完成非线性计算的过程,避免了数值计算在单元断裂时的奇异性,模拟了岩石的峰后破裂行为。对上述算法利用FLAC中的FISH语言开发了岩石破裂化计算程序,并将该程序成功应用于巴西劈裂和中间裂隙拉伸试验的破裂模拟过程中,其模拟结果与相应的理论和试验结果吻合较好,说明该方法对于模拟岩石破裂过程的正确性和可行性。

关键词: 应变能密度, 能量耗散, 岩石破裂, FLAC2D

Abstract: For the nonlinear fracture phenomena of rock, many numerical methods are carried out for study, but the numerical simulation methods for cracked rock are still being studied. Therefore, this paper combines the bilinear strain softening constitutive model with energy dissipation principle to establish damage constitutive equations; and through the strain energy density theory, energy failure criterion of mesoscopic element is established. When the strain energy of an element stored exceeds a fixed value, the element comes into damage situation and the damage degree increases with the increasing of energy dissipation. Simultaneously, the material properties of the damaged element change until it becomes the element with residual strength. As load increases, the element’s damage degree increases. When the strain energy of the element stored exceeds the established value of the energy criterion, the element gets into fracture. With the increasing of fractured elements, macro-crack generates. This method implements linear calculation instead of the process of nonlinear computation, avoids the singularity of the numerical calculation in element fracture, and simulates the post-peak fracture behaviour of rock. The algorithm can be used to develop a rock fracture process computer program by the FISH language in the fast Lagrangian analysis of continual (FLAC). This program is successfully applied to the simulation of Brazilian splitting and tensile test with cracks. The simulation results agree well with the experimental results. It indicates that the simulation method of rock failure process is correct and feasible.

Key words: strain energy density, energy dissipation, rock fracture, fast Lagrangian analysis of continua in 2D (FLAC2D)

中图分类号: 

  • O 241
[1] 赵国彦, 李振阳, 吴浩, 王恩杰, 刘雷磊. 含非贯通裂隙砂岩的动力破坏特性研究[J]. 岩土力学, 2019, 40(S1): 73-81.
[2] 陈炳瑞, 吴昊, 池秀文, 刘辉, 伍梦蝶, 晏俊伟, . 基于STA/LTA岩石破裂微震信号实时识 别算法及工程应用[J]. 岩土力学, 2019, 40(9): 3689-3696.
[3] 马德鹏, 周岩, 刘传孝, 商岩冬, . 不同卸围压速率下煤样卸荷破坏能量演化特征[J]. 岩土力学, 2019, 40(7): 2645-2652.
[4] 赵振华, 张晓君, 李晓程, . 含卸压孔硬岩应力松弛特性试验研究[J]. 岩土力学, 2019, 40(6): 2192-2199.
[5] 肖晓春, 樊玉峰, 吴迪, 丁鑫, 王磊, 赵宝友, . 组合煤岩破坏过程能量耗散特征及冲击危险评价[J]. 岩土力学, 2019, 40(11): 4203-4212.
[6] 汪轶群,洪 义,国 振,王立忠, . 南海钙质砂宏细观破碎力学特性[J]. , 2018, 39(1): 199-206.
[7] 姚旭龙,张艳博,刘祥鑫,梁 鹏,孙 林, . 岩石破裂声发射关键特征信号优选方法[J]. , 2018, 39(1): 375-384.
[8] 何明明,李 宁,陈蕴生,朱才辉. 不同循环加载条件下岩石阻尼比和阻尼系数研究[J]. , 2017, 38(9): 2531-2538.
[9] 张鹏海,杨天鸿,徐 涛,于庆磊,周靖人,赵永川. 蚀变花岗片麻岩破坏过程中声发射事件的演化规律[J]. , 2017, 38(8): 2189-2197.
[10] 杜海民,马 巍,张淑娟,周志伟. 围压与含水率对冻结砂土破坏应变能密度影响特性研究[J]. , 2017, 38(7): 1943-1950.
[11] 谢 璨,李树忱,平 洋,李景龙,李术才,. 峰后裂隙岩石非线性损伤特性与数值模拟研究[J]. , 2017, 38(7): 2128-2136.
[12] 杜海民,马 巍,张淑娟,周志伟. 三轴循环加卸载条件下高含冰冻结砂土变形特性试验研究[J]. , 2017, 38(6): 1675-1681.
[13] 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171.
[14] 周 磊,朱哲明,刘 邦. 裂纹对直墙拱形隧道围岩损伤破坏模式的影响规律研究[J]. , 2017, 38(12): 3688-3697.
[15] 肖晓春,丁 鑫,赵 鑫,潘一山,王爱文,王 磊,. 加载速率影响的煤体破裂过程声-电荷试验研究[J]. , 2017, 38(12): 3419-3426.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 张文杰,陈云敏. 垃圾填埋场抽水试验及降水方案设计[J]. , 2010, 31(1): 211 -215 .
[2] 宫伟力,安里千,赵海燕,毛灵涛. 基于图像描述的煤岩裂隙CT图像多尺度特征[J]. , 2010, 31(2): 371 -376 .
[3] 万 智,董 辉,刘宝琛. 基于正交设计下SVM滑坡变形时序回归预测的超参数选择[J]. , 2010, 31(2): 503 -508 .
[4] 王明年,郭 军,罗禄森,喻 渝,杨建民,谭忠盛. 高速铁路大断面黄土隧道深浅埋分界深度研究[J]. , 2010, 31(4): 1157 -1162 .
[5] 胡勇刚,罗 强,张 良,黄 晶,陈亚美. 基于离心模型试验的水泥土搅拌法加固斜坡软弱土地基变形特性分析[J]. , 2010, 31(7): 2207 -2213 .
[6] 谭峰屹,姜志全,李仲秋,颜惠和. 附加质量法在昆明新机场填料压实密度检测中的应用研究[J]. , 2010, 31(7): 2214 -2218 .
[7] 柴 波,殷坤龙,肖拥军. 巴东新城区库岸斜坡软弱带特征[J]. , 2010, 31(8): 2501 -2506 .
[8] 王维铭,孙 锐,曹振中,袁晓铭. 国内外地震液化场地特征对比研究[J]. , 2010, 31(12): 3913 -3918 .
[9] 杨召亮,孙冠华,郑 宏. 基于潘氏极大值原理的边坡稳定性的整体分析法[J]. , 2011, 32(2): 559 -563 .
[10] 王光进,杨春和,张 超,马洪岭,孔祥云,侯克鹏. 超高排土场的粒径分级及其边坡稳定性分析研究[J]. , 2011, 32(3): 905 -913 .