›› 2005, Vol. 26 ›› Issue (S2): 189-195.

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

单轴拉伸岩样破坏过程及尺寸效应数值模拟

王学滨   

  1. 辽宁工程技术大学 力学与工程科学系,辽宁 阜新123000
  • 收稿日期:2005-05-11 发布日期:2005-12-16
  • 作者简介:王学滨,男, 1975年生,博士,讲师。主要从事非均质材料(混凝土、岩土、韧性金属等)变形、破坏、稳定性理论研究及变形局部化的数值计算工作。
  • 基金资助:
    国家自然科学青年基金项目(No.50309004).

Numerical simulation of failure process and size effect of rock specimen in uniaxial tension

WANG Xue-bin   

  1. Department of Mechanics and Engineering Sciences, Liaoning Technical University, Fuxin 123000, China
  • Received:2005-05-11 Published:2005-12-16

摘要: 由于从实验及理论角度研究岩样单轴拉伸条件下的破坏全过程及尺寸效应难度都很大。因此采用拉格朗日元法来研究这些问题。在峰值强度之前后,岩石材料的本构模型分别取为线弹性及拉破坏线性应变软化模型。为了使拉伸塑性区不出现在试样的端部,在试样的两侧面中部预制了2个凹槽。数值模拟结果表明,全程拉应力-拉应变曲线分为峰前和峰后阶段。在接近峰值的峰前阶段,由于两凹槽附近具有明显的拉应力集中现象,拉伸塑性区最先出现在两凹槽附近。随着轴向拉应变的增加,发生拉伸破坏的单元的数目增加,新发生拉伸破坏的单元越来越接近试样的中心,直到两块拉伸塑性区在应变软化阶段贯通。两凹槽连线上各单元拉应力的分布呈现3个阶段,“澡盆型”(“U型”)阶段,“双峰型”(“M型”)阶段及“单峰型”(“Π型”)阶段。“澡盆型”阶段对应于全程拉应力-拉应变曲线的弹性阶段。“双峰型”阶段及“单峰型”阶段对应于全程拉应力-拉应变曲线的非弹性阶段(包括峰值强度之前的一小段,即应变硬化阶段及峰后的应变软化阶段)。增加试样高度及降低试样宽度,拉应力-拉应变曲线的软化段变得越来越陡峭,因而试样越容易发生失稳破坏。由于试样宽度较大时,试样内部的单元并非处于单向拉应力状态,因此,增加试样宽度,全程拉应力-拉应变曲线的峰值强度增加。当试样宽度较小时,从出现塑性区,到塑性区贯通所需要的时间步较小,或应变范围较窄。这说明试样的脆性较强,前兆不明显。前兆不明显的脆性破坏对应常见的是洞室岩爆、冲击地压及地震等灾害。

关键词: 单轴拉伸, 应变软化, 塑性区, 尺寸效应, 破坏过程, 应力-应变曲线, 失稳

Abstract: Due to the fact that there is a great deal of experimental and analytical difficulty in investigating the failure process and the size effect of rock specimen in uniaxial tension, in the present paper these problems are modeled numerically by FLAC. In elastic and strain-softening stages, the adopted constitutive relations of rock are linear elastic and linear strain-softening due to tensile failure, respectively. To ensure that elements subjected to tensile failure are localized into a narrow band within the specimen, two null elements at two edges of specimen are prescribed to model two notches. It is found from numerical results that the complete tensile stress-tensile strain curve exhibits pre-peak and post-peak branches. Prior to the peak stress, tensile stress concentration appears near the two notches so that two elements close to the two notches due to tension yield. The number of yielded elements is increased with axial tensile strain and the new plastic elements are even closer to the center of specimen until a horizontal and longer plastic zone is formed. Tensile stress distribution undergoes three distinct stages: upward concave stage, stage of downward concavity with two peaks and stage of downward concavity with one peak. The first stage corresponds to the elastic stage of complete stress-strain curve; however, the latter two stages correspond to the inelastic stage including strain-hardening and strain-softening stages. Softening branch of complete stress-strain curve becomes steeper so that the instability of rock specimen easily occurs as the height of specimen is increased or the width of specimen is decreased. Peak stress of complete stress-strain curve increases with the width of specimen, which is due to the fact that the stress state within the specimen is not a simple uniaxial state of stress for wider specimen. For thinner specimen, the needed time steps form the initial formation of yielded elements to the full development of a plastic zone are less than those for wider specimen. This reflects that the precursor for thinner specimen is less apparent. Brittle failure without obvious precursor is concerned with seriously natural hazards, such as rockburst and earthquake.

Key words: uniaxial tension, strain-softening, plastic zone, size effect, failure process, stress-strain curve, instability

中图分类号: 

  • TU452
[1] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[2] 徐 鹏, 蒋关鲁, 雷 涛, 刘 琪, 王智猛, 刘 勇, . 考虑填土强度的加筋土挡墙动位移计算[J]. 岩土力学, 2019, 40(5): 1841-1846.
[3] 朱仁杰, 车爱兰, 严 飞, 文 海, 葛修润, . 含贯通性结构面岩质边坡动力演化规律[J]. 岩土力学, 2019, 40(5): 1907-1915.
[4] 王凤云, 钱德玲, . 基于统一强度理论深埋圆形隧道围岩的剪胀分析[J]. 岩土力学, 2019, 40(5): 1966-1976.
[5] 王娟娟, 郝延周, 王铁行. 非饱和压实黄土结构特性试验研究[J]. 岩土力学, 2019, 40(4): 1351-1357.
[6] 蒋海明, 李 杰, 王明洋, . 块系岩体滑移失稳中低摩擦效应的理论与试验研究[J]. 岩土力学, 2019, 40(4): 1405-1412.
[7] 王 腾, 吴 瑞. 黏土中海底管线竖向贯入阻力研究[J]. 岩土力学, 2019, 40(3): 871-878.
[8] 宋洪强, 左建平, 陈 岩, 黎立云, 洪紫杰, . 基于岩石破坏全过程能量特征改进的能量跌落系数[J]. 岩土力学, 2019, 40(1): 91-98.
[9] 巨能攀, 邓天鑫, 李龙起, 蒋金阳, 张陈羊. 强震作用下陡倾顺层斜坡倾倒变形 机制离心振动台试验[J]. 岩土力学, 2019, 40(1): 99-108.
[10] 王文沛, 李 滨, 冯 振, 张博文, 高 杨, . 考虑场地效应的高陡岩质斜坡地震失稳机制[J]. 岩土力学, 2019, 40(1): 297-304.
[11] 秦庆词, 李克钢, 杨宝威, 王 庭, 张雪娅, 郭 文. 岩石全应力-应变过程关键特征点损伤特征分析[J]. 岩土力学, 2018, 39(S2): 14-24.
[12] 高桂云,王成虎,王春权,. 双圆环直接拉伸试验试样最优尺寸范围研究[J]. , 2018, 39(S1): 191-202.
[13] 李术才,潘东东,许振浩,李利平,林 鹏,袁永才,高成路,路 为, . 承压型隐伏溶洞突水灾变演化过程模型试验[J]. , 2018, 39(9): 3164-3173.
[14] 王凤云,钱德玲. 基于切向应变软化的深埋圆形隧道围岩弹塑性分析[J]. , 2018, 39(9): 3313-3320.
[15] 高 樯,温 智,王大雁,牛富俊,谢艳丽,苟廷韬,. 基于冻融交界面直剪试验的冻土斜坡失稳过程研究[J]. , 2018, 39(8): 2814-2822.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 徐远杰,潘家军,刘祖德. 混凝土面板堆石坝的一种坝坡修整算法[J]. , 2009, 30(10): 3139 -3144 .
[2] 杨建平,陈卫忠,田洪铭,于洪丹. 应力-温度对低渗透介质渗透率影响研究[J]. , 2009, 30(12): 3587 -3594 .
[3] 宫伟力,安里千,赵海燕,毛灵涛. 基于图像描述的煤岩裂隙CT图像多尺度特征[J]. , 2010, 31(2): 371 -376 .
[4] 贾明涛,王李管. 基于区域化变量及RMR评价体系的金川III矿区矿岩质量评价[J]. , 2010, 31(6): 1907 -1912 .
[5] 柴 波,殷坤龙,肖拥军. 巴东新城区库岸斜坡软弱带特征[J]. , 2010, 31(8): 2501 -2506 .
[6] 文松霖,徐文强. GRF基础承载机制试验研究[J]. , 2010, 31(S1): 35 -40 .
[7] 李建军,梁仁旺. 水泥土抗压强度和变形模量试验研究[J]. , 2009, 30(2): 473 -477 .
[8] 孙 勇. 利用旁压试验确定抗滑桩地基系数的研究[J]. , 2009, 30(3): 699 -702 .
[9] 苏国韶,张克实,吕海波. 位移反分析的粒子群优化-高斯过程协同优化方法[J]. , 2011, 32(2): 510 -515 .
[10] 牛文杰,叶为民,陈 宝. 高压实膨润土的非饱和渗透膨胀模型[J]. , 2009, 30(S2): 88 -92 .