›› 2012, Vol. 33 ›› Issue (6): 1679-1688.

• 基础理论与实验研究 • 上一篇    下一篇

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

袁小平,刘红岩,王志乔   

  1. 中国地质大学(北京) 工程技术学院,北京 100083
  • 收稿日期:2010-12-14 出版日期:2012-06-11 发布日期:2012-06-14
  • 通讯作者: 刘红岩,男,副教授,主要从事岩土力学数值方面的研究工作。E-mail: lhyan1204@126.com E-mail:袁小平,男,1987年生,主要从事岩土力学数值方面的研究工作
  • 作者简介:袁小平,男,1987年生,主要从事岩土力学数值方面的研究工作
  • 基金资助:

    国家自然科学基金资助(No. 41002113,No. 41162009 );中央高校基本科研业务费专项资金资助(No. 2011PY0190,No. 2010ZY45);教育部科学技术研究重点项目(No. 2011175)。

Crack tip plastic zone and elastoplastic fracture model for nonpenetrating jointed rock mass under uniaxial compression

YUAN Xiao-ping,LIU Hong-yan,WANG Zhi-qiao   

  1. School of Engineering & Technology, China University of Geosciences (Beijing), Beijing 100083, China
  • Received:2010-12-14 Online:2012-06-11 Published:2012-06-14

摘要: 基于Drucker-Prager(下简称D-P)准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。

关键词: 节理岩石, 应变能密度, 塑性软化, 断裂软化因子, 隐式积分算法, D-P准则, 尖端塑性区

Abstract: An elastoplastic fracture model for jointed rock mass under compressive loading is presented. In order to account for the mixed mode plastic zone model of wing crack tip of jointed rock under small scale yielding, the dimensionless function of the plastic zone diameter based on the Drucker-Prager (D-P) criterion is derived. The plastic zone models of mode I, mode II and mixed-mode (I+II) are compared with plastic zones obtained by Mises criterion and the influences of different Poisson's ratios for mixed-mode plastic zones of wing crack tip are studied. It is shown that the mode I and mixed-mode plastic zones of D-P criterion are larger than that of Mises criterion; and its mode II and mixed-mode plastic zones are not continuous between the upper and lower surfaces of wing crack. Moreover, plastic yield criterion is employed simultaneously with the fracture softening to describe the behaviors of jointed rock based on D-P criterion. The conception of fracture softening factor is introduced to characterize the softening feature of crack propagation. Regarding the strain energy density of wing cracks at joint-tip is responsible for the evolution of fracture softening of jointed rock. Thus, the behavior of fracture softening can be developed by proposed function of strain energy density. Borja’s strain hardening/softening function is employed as the plastic yield function, regarding the plastic inner variables and stress as two impact factors of the hardening function. The elastoplastic fracture model with its numerical algorithm is proposed; and the code of elastoplastic fracture model is implemented by using return mapping implicit integration algorithm. The influences of friction coefficients, initial joints orientation and stress intensity factor on the constitutive response under uniaxial compression are studied. The results are generally in good agreement with experimental data of uniaxial compression tests of jointed marble.

Key words: jointed rock mass, strain energy density, plastic softening, fracture softening factor, implicit integration algorithm, D-P criterion, crack tip plastic zone

中图分类号: 

  • TU 45
[1] 杜海民,马 巍,张淑娟,周志伟. 围压与含水率对冻结砂土破坏应变能密度影响特性研究[J]. , 2017, 38(7): 1943-1950.
[2] 谢 璨,李树忱,平 洋,李景龙,李术才,. 峰后裂隙岩石非线性损伤特性与数值模拟研究[J]. , 2017, 38(7): 2128-2136.
[3] 王仁超,曹婷婷,刘严如. 亚塑性模型不同积分算法的数值实现[J]. , 2017, 38(5): 1510-1516.
[4] 周永强,盛 谦,朱泽奇,付晓东. 基于修正Drucker-Prager准则的岩石次加载面模型[J]. , 2017, 38(2): 400-408.
[5] 李 舰,韦昌富,刘 艳,. 非饱和黏土状态相关本构模型的数值实现[J]. , 2017, 38(10): 2799-2808.
[6] 李晓锋 ,李海波 ,夏 祥 ,刘 博 ,冯海鹏,. 类节理岩石直剪试验力学特性的数值模拟研究[J]. , 2016, 37(2): 583-591.
[7] 崔 岚,郑俊杰,章荣军,张 威. 弹塑性软化模型下隧洞围岩变形与支护压力分析[J]. , 2014, 35(3): 717-722.
[8] 刘 艳 ,韦昌富 ,房 倩 ,陈 盼,. 非饱和土水-力本构模型及其隐式积分算法[J]. , 2014, 299(2): 365-370.
[9] 施明明,张友良,谭 飞. 修正应变能密度因子准则及岩石裂纹扩展研究[J]. , 2013, 34(5): 1313-1318.
[10] 袁小平,刘红岩,王志乔. 基于Drucker-Prager准则的岩石弹塑性损伤本构模型研究[J]. , 2012, 33(4): 1103-1108.
[11] 徐 磊,任青文. 岩体结构面新型本构模型研究[J]. , 2011, 32(S1): 217-224.
[12] 孙 倩,李树忱,冯现大,李文婷,袁 超. 基于应变能密度理论的岩石破裂数值模拟方法研究[J]. , 2011, 32(5): 1575-1582.
[13] 周 雷,张洪武. 非饱和土化学-塑性耦合本构行为的数值模拟[J]. , 2009, 30(7): 2133-2140.
[14] 姚国圣,李镜培,谷拴成. 考虑岩体扩容和塑性软化的软岩巷道变形解析[J]. , 2009, 30(2): 463-467.
[15] 张 凯,周 辉,冯夏庭,房敬年,张元刚. Biot固结理论中连续性方程形式的讨论[J]. , 2009, 30(11): 3273-3277.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 崔皓东,朱岳明. 二滩高拱坝坝基渗流场的反演分析[J]. , 2009, 30(10): 3194 -3199 .
[2] 何思明,吴 永,李新坡. 嵌岩抗拔桩作用机制研究[J]. , 2009, 30(2): 333 -337 .
[3] 刘清秉,项 伟,张伟锋,崔德山. 离子土壤固化剂改性膨胀土的试验研究[J]. , 2009, 30(8): 2286 -2290 .
[4] 况雨春,伍开松,杨迎新,马德坤. 三牙轮钻头破岩过程计算机仿真模型[J]. , 2009, 30(S1): 235 -238 .
[5] 杜文琪,王 刚. 土工结构地震滑动位移统计分析[J]. , 2011, 32(S1): 520 -0525 .
[6] 鄢治华,刘志伟,刘厚健. 黄河阶地上某电厂高边坡参数选取及其工程治理[J]. , 2009, 30(S2): 465 -468 .
[7] 许振浩 ,李术才 ,李利平 ,侯建刚 ,隋 斌 ,石少帅. 基于层次分析法的岩溶隧道突水突泥风险评估[J]. , 2011, 32(6): 1757 -1766 .
[8] 温世清 ,刘汉龙 ,陈育民. 浆固碎石桩单桩荷载传递特性研究[J]. , 2011, 32(12): 3637 -3641 .
[9] 李顺群 ,高凌霞 ,柴寿喜. 冻土力学性质影响因素的显著性和交互作用研究[J]. , 2012, 33(4): 1173 -1177 .
[10] 钟 声 ,王川婴 ,吴立新 ,唐新建 ,王清远. 点状不良地质体钻孔雷达响应特征 ——围岩及充填效应正演分析[J]. , 2012, 33(4): 1191 -1195 .