CN107101887A - 一种声发射与数值计算相结合的岩石破裂过程分析方法 - Google Patents
一种声发射与数值计算相结合的岩石破裂过程分析方法 Download PDFInfo
- Publication number
- CN107101887A CN107101887A CN201710319573.4A CN201710319573A CN107101887A CN 107101887 A CN107101887 A CN 107101887A CN 201710319573 A CN201710319573 A CN 201710319573A CN 107101887 A CN107101887 A CN 107101887A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- rock
- damage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000011835 investigation Methods 0.000 title claims abstract description 6
- 230000006378 damage Effects 0.000 claims abstract description 73
- 238000011161 development Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 208000027418 Wounds and injury Diseases 0.000 claims description 4
- 238000012669 compression test Methods 0.000 claims description 4
- 208000014674 injury Diseases 0.000 claims description 4
- 230000001066 destructive effect Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 239000004575 stone Substances 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000003902 lesion Effects 0.000 abstract description 2
- 239000000463 material Substances 0.000 description 6
- 238000005065 mining Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 4
- 208000037656 Respiratory Sounds Diseases 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0067—Fracture or rupture
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明专利提出了一种声发射与数值计算相结合的岩石破裂过程分析方法,将声发射\微震监测与数值计算结合起来对岩石损伤与破裂过程进行分析。首先,利用声发射\微震信息反演破裂源损伤尺寸和程度,对岩石的弹性模量和强度进行相应地弱化;然后将这些参数输入数值模型之中,再进行数值计算分析。利用该方法,可以获取真实应力场并有效预测岩石破裂模式。本发明为研究岩石损伤与破裂提供了一种精确可靠的新方法。
Description
技术领域
本发明属于采矿工程与岩石工程领域,介绍一种利用声发射\微震监测与数值计算相结合的岩石破裂过程分析的方法。
背景技术
随着地表资源日趋减少,深部开采已成为未来采矿发展的趋势。在深部开采条件下,岩爆、塌方、冒顶、片帮等灾害日趋凸显,已经成为制约我国乃至世界深部矿产资源开采的一大瓶颈。因此,深入开展岩石损伤与破裂及其致灾机理的研究,寻求有关失稳的前兆规律,对于确保矿山安全和稳定生产具有重要意义。
数值计算是研究岩石损伤与破裂常用的技术手段之一。但岩石是一种不同于其他固体的特殊材料,具有非均质、非连续、非弹性、非线性和各向异性。岩体性质的准确表征一直是数值计算模型的难点问题。目前,通常采用概率随机分布描述岩石的非均匀性,然而,很少有文献能说明这种方法能在多大程度上真实地体现岩石的非均质性。因此,需要发展一种新的岩体性质表征方法和思路。
岩石的损伤破裂不仅与其内部的结构和成分有关,而且更为重要的是与受载后的损伤演化过程直接相关。岩石的早期损伤分布多取决于初始的细观结构,而后期的损伤则由前期损伤所引起的应力重分布决定。因此,忽略岩石真实细观结构,重点关注岩石的损伤演化可能是一个解决岩石力学问题的新思路。
岩石类脆性材料在受力或变形时将释放声波和次声波,它是脆性材料在变形破坏过程中伴随着局部微破裂的出现而产生的弹性波。材料学家将这种现象称之为声发射;在工程岩体中,采矿工程师称之为微震。岩石的声发射活动反映了岩石内部的损伤演化状况,但大部分研究仅限于定性评价。如何更有效地解读和利用声发射\微震数据以量化岩石的损伤并预测破坏依然是一个挑战。
为此,本发明专利提出了一种声发射\微震数据驱动的岩石损伤与破裂计算方法,将声发射\微震监测与数值计算结合起来对岩石损伤与破裂过程进行分析。首先,利用声发射\微震数据对岩体损伤进行定量表征,对岩石的弹性模量和强度进行相应地弱化;然后将这些参数输入数值模型之中,再进行数值计算分析。利用该方法,可以获取真实应力场并有效预测岩石破裂模式。
发明内容
本发明的目的是在于提出一种声发射\微震与数值计算相结合的岩石破裂过程分析方法。具体技术方案如下:
一种声发射与数值计算相结合的岩石破裂过程分析方法,步骤如下:
(1)选用一种岩石试件进行单轴压缩试验;
(2)采用声发射技术监测其破坏过程,以岩石试件底部圆心作为坐标原点,确定8个传感器位置为四个一组,每组设置在同一水平面上且均匀分布;
(3)将岩石试件放置于压力机的承压板间,先开始采集数据,再对上下承压板施加压力,直至试件破坏;
(4)基于声发射数据的岩石损伤定量表征方法:
震源半径方程为:
损伤变量D:
式中:σn表示垂直裂纹的应力,a表示裂纹长度的一半,E表示岩石的弹性模量,v表示岩石的泊松比,γ表示单位面积表面能,Ek表示动能,Ei表示内能,Ea表示表面能,D损伤变量;
(5)基于数值计算的损伤模型:
单元的损伤变量D按照如下关系给出:
εt0和εc0分别是当单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变,ε1和ε3分别是第一、第三主应变,n是表示单元损伤演化的系数;
单元的弹性模量按照如下关系式给出:E=(1-D)E0
E0和E分别为损伤前和损伤后的弹性模量;
单元的强度也按如上关系给出:σs=(1-D)σs0
其中σs0和σs分别为损伤前和损伤后的单元强度。
(6)按基于声发射数据的岩石损伤定量表征方法分析监测数据,反演岩石的损伤半径和损伤程度;利用式
E=(1-D)E0和σs=(1-D)σs0对岩石的弹性模量和强度进行相应地弱化;
(7)将声发射反演损伤作为岩石初始条件,按基于数值计算的损伤模型进行数值计算,即可得到岩石的最终破裂模式。
本发明改变了传统采用地震学理论反演室内试验和工程尺度破裂源尺寸的不恰当方法,提出了新的小尺度破裂源损伤范围和程度反演方法;突破了过去传统分别单独采用声发射\微震监测和数值计算分析岩石损伤与破裂的研究方法,将声发射\微震监测与数值计算结合起来对岩石损伤与破裂过程进行分析。本发明为研究岩石损伤与破裂提供了一种精确可靠的新方法。
附图说明
图1为传感器布置示意图,图中标号说明:1-1号传感器;2-2号传感器;3-3号传感器;4-4号传感器;5-5号传感器;6-6号传感器;7-7号传感器;8-8号传感器;x-x坐标;y-y坐标;z-z坐标。
图2-1为岩石损伤结果分布图,图2-2为基于声发射数据的反演损伤分布图;
图3-1为岩石最终破裂模式室内试验结果示意图,图3-2为岩石最终破裂模式数值模拟结果示意图。
具体实施方式
以下将结合图和具体实施过程对本发明做进一步详细说明:
本发明的原理如下:
(1)基于声发射监测的损伤模型
在声发射数据驱动模型中,假设每一个声发射事件对应一个裂纹。裂纹的长度等于声发射破裂源尺寸,即损伤区范围。一般认为声发射源破裂尺寸与其能量密切相关。破裂尺寸越大,释放的能量越大。微震和声发射数据处理通常借鉴地震学理论。地震学中的一个基本假设就是,岩石断裂是由剪切错动引起的。一个常用的源模型就是Brune模型,它假设震源为一个半径为r0的圆形断层面。震源半径可以通过以下公式计算
其中为cs介质横波波速,fc为横波拐角频率。但是,根据剪切模型得到的源尺寸非常大,与现场观测结果严重不符。造成这种结果的原因有两点:一是因为室内和现场声发射事件与天然地震的产生机理不同。天然地震是由断层剪切错动引起的,而室内和现场的岩石破裂以拉伸裂纹为主,且裂纹扩展方向与最大压应力方向平行。这个观点已经被许多研究者证实。另一个原因是剪切模型忽略了材料断裂表面能,不符合能量守恒定律。事实上,断裂表面能是岩石破裂过程中能量消耗的重要组成部分,其量级可能远远大于声发射事件能量。
因此,将地震学和断裂力学相结合,蔡明等提出了一种基于拉伸破坏和能量守恒的声发射源尺寸估计新模型。通过与现场观测进行对比,该方法得到的源尺寸更加合理。该模型假定拉伸断裂为岩石的主要失效机制,并且测得的声发射能量为总能量的一部分,它等于弹性应变能变化与表面能之差。忽略温度影响,热力学第一定律可以表示为:
Ei=Ea+Ed+Ek (2)
式中Ei表示内能,Ea表示表面能,Ed表示耗散能,Ek表示动能。内能为弹性应变能,一般认为是驱动岩石破裂的能量来源;当新的表面产生时,表面能将发生改变,它是裂纹扩展的最重要阻力;耗散能形式是多样的,但主要是塑性变形和滑动摩擦;动能是由于材料速度产生的。唐春安认为岩石细观非均匀性是造成其宏观非线性的根本原因,岩石性质在细观单元上表现为弹脆性。而且,一般为声发射是弹性应变能的释放造成的,这意味着声发射源应该是弹脆性的。因此,耗散能可以忽略
Ei=Ea+Ek (3)
岩石储存的内能,可根据下式计算:
其中,σn表示垂直裂纹的应力,a表示裂纹长度的一半,E表示岩石的弹性模量,v表示岩石的泊松比。表面能表达式为
其中γ表示单位面积表面能,KIC表示I型裂纹断裂韧度。因此,震源半径方程为
假定岩石的损伤程度与其能量变化密切相关。断裂表面能是阻碍岩石断裂的主要因素,其大小表征了岩石损伤程度的强弱。因此,定义了损伤变量D
(2)基于数值计算的损伤模型
假设岩体是理想弹性体,不会产生塑性变形,则满足如下运动微分方程:
式中,G表示剪切模量(Pa);v表示泊松比;ui(i=x,y,z)表示位移(m);t表示时间(s);Fi表示i方向体力分量(N/m3);ρ表示岩体密度。
采用最大拉应力准则和莫尔库伦准则作为破坏判断准则,可分别表达为:
则当介质的应力状态满足了最大拉应力准则和莫尔库伦准则时,其分别发生拉伸损伤和剪切损伤,即:
F1=0或F2=0 (10)
其中,ft0和fc0分别是单元的单轴抗拉和单轴抗压强度,σ1和σ3分别是第一、第三主应力,φ为内摩擦角。F1和F2是两个表示应力状态的函数,其数值大于零表示介质分别发生了拉伸和剪切损伤。这里需要说明的是,在任何受力条件下,拉伸损伤是优先判断的。当(10)式满足时,单元的损伤变量D按照如下关系给出:
这里εt0和εc0分别是当单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变,ε1和ε3分别是第一、第三主应变,n是表示单元损伤演化的系数。
(3)损伤对岩石力学参数的影响
按照弹性损伤理论,单元的弹性模量按照如下关系式给出:
E=(1-D)E0 (12)
其中E0和E分别为损伤前和损伤后的弹性模量。这里假定损伤及其演化都是各向同性的,故E、E0和D都是标量。
同理,单元的强度(单轴抗压强度和抗拉强度)也按如上关系给出:
σs=(1-D)σs0 (13)
其中σs0和σs分别为损伤前和损伤后的单元强度。这里假定损伤及其演化都是各向同性的,故E、E0、σs、σs0和D都是标量。
具体操作步骤如下:
(1)选用一种岩石进行单轴压缩试验。单轴压缩试件直径50mm,高度100mm,严格按照国际岩石力学学会(ISRM)建议方法进行加工。
(2)采用声发射技术监测其破坏过程,以岩石试件底部圆心作为坐标原点,确定8个传感器位置(见图1)。
(3)将岩石试件放置于压力机的承压板间。先开始采集数据,再对上下承压板施加压力,直至试件破坏。
(4)按基于声发射数据的岩石损伤定量表征方法分析监测数据,反演岩石的损伤半径和损伤程度;利用式
E=(1-D)E0 (12)
和
σs=(1-D)σs0 (13)
对岩石的弹性模量和强度进行相应地弱化(见图2-1、图2-2)。
(5)将声发射反演损伤作为岩石初始条件,按基于数值计算的损伤模型中的公式(8-13)进行数值计算,即可得到岩石的最终破裂模式(见图3-1、图3-2)。
Claims (1)
1.一种声发射与数值计算相结合的岩石破裂过程分析方法,其特征在于步骤如下:
(1)选用一种岩石试件进行单轴压缩试验;
(2)采用声发射技术监测其破坏过程,以岩石试件底部圆心作为坐标原点,确定8个传感器位置,四个一组,每组设置在同一水平面上且均匀分布;
(3)将岩石试件放置于压力机的承压板间,先开始采集数据,再对上下承压板施加压力,直至试件破坏;
(4)基于声发射数据的岩石损伤定量表征方法;
震源半径方程为;
<mrow>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<msubsup>
<mi>&pi;&sigma;</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mrow>
<mi>E</mi>
</mfrac>
<mo>-</mo>
<mn>4</mn>
<mi>&gamma;</mi>
<mi>a</mi>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
损伤变量D;
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>E</mi>
<mi>a</mi>
</msub>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
式中:σn表示垂直裂纹的应力,a表示裂纹长度的一半,E表示岩石的弹性模量,v表示岩石的泊松比,γ表示单位面积表面能,Ek表示动能,Ei表示内能,Ea表示表面能,D损伤变量;
(5)基于数值计算的损伤模型;
单元的损伤变量D按照如下关系给出:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mo>|</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>t</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
</mfrac>
<msup>
<mo>|</mo>
<mi>n</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>dF</mi>
<mn>1</mn>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mo>|</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>c</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>3</mn>
</msub>
</mfrac>
<msup>
<mo>|</mo>
<mi>n</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>dF</mi>
<mn>2</mn>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中εt0和εc0分别是当单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变,ε1和ε3分别是第一、第三主应变,n是表示单元损伤演化的系数;
单元的弹性模量按照如下关系式给出:E=(1-D)E0
E0和E分别为损伤前和损伤后的弹性模量;
单元的强度也按如上关系给出:σs=(1-D)σs0
其中σs0和σs分别为损伤前和损伤后的单元强度
(6)按基于声发射数据的岩石损伤定量表征方法分析监测数据,反演岩石的损伤半径和损伤程度;利用式E=(1-D)E0和σs=(1-D)σs0对岩石的弹性模量和强度进行相应地弱化;E0和E分别为损伤前和损伤后的弹性模量;σs0和σs分别为损伤前和损伤后的单元强度;
(7)将声发射反演损伤作为岩石初始条件,按基于数值计算的损伤模型进行数值计算,即可得到岩石的最终破裂模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710319573.4A CN107101887B (zh) | 2017-05-09 | 2017-05-09 | 一种声发射与数值计算相结合的岩石破裂过程分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710319573.4A CN107101887B (zh) | 2017-05-09 | 2017-05-09 | 一种声发射与数值计算相结合的岩石破裂过程分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107101887A true CN107101887A (zh) | 2017-08-29 |
CN107101887B CN107101887B (zh) | 2019-06-07 |
Family
ID=59668944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710319573.4A Expired - Fee Related CN107101887B (zh) | 2017-05-09 | 2017-05-09 | 一种声发射与数值计算相结合的岩石破裂过程分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107101887B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115888A (zh) * | 2018-09-03 | 2019-01-01 | 太原理工大学 | 一种声发射探头类型的选择方法 |
CN109187185A (zh) * | 2018-09-21 | 2019-01-11 | 中国石油大学(华东) | 一种聚合物损伤演化测试方法及系统 |
CN111007231A (zh) * | 2019-12-23 | 2020-04-14 | 鞍钢集团矿业有限公司 | 一种采动岩体内部新生剪切裂纹尺度的量化方法 |
CN111272883A (zh) * | 2020-03-04 | 2020-06-12 | 成都理工大学 | 一种基于声发射模型的岩石破裂模式智能探测识别方法 |
CN111322974A (zh) * | 2020-03-05 | 2020-06-23 | 北京理工大学 | 一种金属药型罩材料聚能射流平均直径的预测方法及其应用 |
CN111855410A (zh) * | 2020-08-01 | 2020-10-30 | 江西理工大学 | 一种尾砂充填组合体的弹性模量计算及破坏特性分析方法 |
CN112444448A (zh) * | 2020-11-17 | 2021-03-05 | 大连理工大学 | 基于聚类分析和信源熵原理的岩石受压微裂纹扩展方向性评价方法 |
CN112765781A (zh) * | 2020-12-30 | 2021-05-07 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN113203627A (zh) * | 2021-05-10 | 2021-08-03 | 北京科技大学 | 一种基于声发射的岩石力学非均质度的测试评价方法 |
CN113281148A (zh) * | 2021-05-19 | 2021-08-20 | 江西理工大学 | 一种岩石蠕变破坏微裂纹贯通阶段起始时间辨识方法 |
CN113283173A (zh) * | 2021-05-31 | 2021-08-20 | 沈阳工业大学 | 地下工程能量与参数综合反分析系统及其方法 |
CN113899879A (zh) * | 2021-10-08 | 2022-01-07 | 中国矿业大学 | 一种深部开挖扰动下岩石破碎后的动能确定方法 |
CN114113335A (zh) * | 2021-12-10 | 2022-03-01 | 东北大学 | 基于声发射/微震监测的岩石耗散能时空分布量化方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183411A (zh) * | 2011-02-15 | 2011-09-14 | 中国水电顾问集团华东勘测设计研究院 | 长期稳定荷载作用下脆性岩石破裂模式的预测方法 |
CN102998180A (zh) * | 2012-11-30 | 2013-03-27 | 北京交通大学 | 一种采用残余强度构建损伤岩石本构关系的方法 |
CN103592175A (zh) * | 2013-11-13 | 2014-02-19 | 北京航空航天大学 | 一种基于连续损伤力学退化模型的复合材料结构失效分析方法 |
CN103760024A (zh) * | 2014-01-29 | 2014-04-30 | 核工业北京地质研究院 | 基于累计声发射撞击数客观确定岩石启裂强度的方法 |
CN104198586A (zh) * | 2014-08-08 | 2014-12-10 | 西北矿冶研究院 | 轴应力下基于波速的岩石损伤变量的确定方法 |
CN104792621A (zh) * | 2015-04-21 | 2015-07-22 | 山东大学 | 一种用于监测岩石试样破坏多元前兆信息的系统及方法 |
CN105547825A (zh) * | 2016-01-21 | 2016-05-04 | 山东大学 | 单轴压缩过程中岩样损伤的监测装置及监测方法 |
CN106199700A (zh) * | 2016-06-30 | 2016-12-07 | 马克 | 一种地下水封储油洞库微震监测方法及系统 |
-
2017
- 2017-05-09 CN CN201710319573.4A patent/CN107101887B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183411A (zh) * | 2011-02-15 | 2011-09-14 | 中国水电顾问集团华东勘测设计研究院 | 长期稳定荷载作用下脆性岩石破裂模式的预测方法 |
CN102998180A (zh) * | 2012-11-30 | 2013-03-27 | 北京交通大学 | 一种采用残余强度构建损伤岩石本构关系的方法 |
CN103592175A (zh) * | 2013-11-13 | 2014-02-19 | 北京航空航天大学 | 一种基于连续损伤力学退化模型的复合材料结构失效分析方法 |
CN103760024A (zh) * | 2014-01-29 | 2014-04-30 | 核工业北京地质研究院 | 基于累计声发射撞击数客观确定岩石启裂强度的方法 |
CN104198586A (zh) * | 2014-08-08 | 2014-12-10 | 西北矿冶研究院 | 轴应力下基于波速的岩石损伤变量的确定方法 |
CN104792621A (zh) * | 2015-04-21 | 2015-07-22 | 山东大学 | 一种用于监测岩石试样破坏多元前兆信息的系统及方法 |
CN105547825A (zh) * | 2016-01-21 | 2016-05-04 | 山东大学 | 单轴压缩过程中岩样损伤的监测装置及监测方法 |
CN106199700A (zh) * | 2016-06-30 | 2016-12-07 | 马克 | 一种地下水封储油洞库微震监测方法及系统 |
Non-Patent Citations (2)
Title |
---|
MARTIN ENTACHER等: "《Rock Failure and Crack Propagation Beneath Disc Cutters》", 《ROCK MECH ROCK ENG》 * |
盖迪等: "《基于损伤力学理论的脉冲压裂模型及模拟研究》", 《采矿与安全工程学报》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115888A (zh) * | 2018-09-03 | 2019-01-01 | 太原理工大学 | 一种声发射探头类型的选择方法 |
CN109187185A (zh) * | 2018-09-21 | 2019-01-11 | 中国石油大学(华东) | 一种聚合物损伤演化测试方法及系统 |
CN111007231A (zh) * | 2019-12-23 | 2020-04-14 | 鞍钢集团矿业有限公司 | 一种采动岩体内部新生剪切裂纹尺度的量化方法 |
CN111007231B (zh) * | 2019-12-23 | 2022-04-15 | 鞍钢集团矿业有限公司 | 一种采动岩体内部新生剪切裂纹尺度的量化方法 |
CN111272883A (zh) * | 2020-03-04 | 2020-06-12 | 成都理工大学 | 一种基于声发射模型的岩石破裂模式智能探测识别方法 |
CN111322974B (zh) * | 2020-03-05 | 2021-03-26 | 北京理工大学 | 一种金属药型罩材料聚能射流平均直径的预测方法及其应用 |
CN111322974A (zh) * | 2020-03-05 | 2020-06-23 | 北京理工大学 | 一种金属药型罩材料聚能射流平均直径的预测方法及其应用 |
CN111855410A (zh) * | 2020-08-01 | 2020-10-30 | 江西理工大学 | 一种尾砂充填组合体的弹性模量计算及破坏特性分析方法 |
CN111855410B (zh) * | 2020-08-01 | 2022-05-03 | 江西理工大学 | 一种尾砂充填组合体的弹性模量计算及破坏特性分析方法 |
CN112444448B (zh) * | 2020-11-17 | 2022-02-18 | 大连理工大学 | 基于聚类分析和信源熵原理的岩石受压微裂纹扩展方向性评价方法 |
CN112444448A (zh) * | 2020-11-17 | 2021-03-05 | 大连理工大学 | 基于聚类分析和信源熵原理的岩石受压微裂纹扩展方向性评价方法 |
CN112765781A (zh) * | 2020-12-30 | 2021-05-07 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN112765781B (zh) * | 2020-12-30 | 2022-07-29 | 中铁二院工程集团有限责任公司 | 层状岩体破坏准则及确定方法 |
CN113203627A (zh) * | 2021-05-10 | 2021-08-03 | 北京科技大学 | 一种基于声发射的岩石力学非均质度的测试评价方法 |
CN113281148A (zh) * | 2021-05-19 | 2021-08-20 | 江西理工大学 | 一种岩石蠕变破坏微裂纹贯通阶段起始时间辨识方法 |
CN113283173A (zh) * | 2021-05-31 | 2021-08-20 | 沈阳工业大学 | 地下工程能量与参数综合反分析系统及其方法 |
CN113283173B (zh) * | 2021-05-31 | 2023-11-24 | 沈阳工业大学 | 地下工程能量与参数综合反分析系统及其方法 |
CN113899879A (zh) * | 2021-10-08 | 2022-01-07 | 中国矿业大学 | 一种深部开挖扰动下岩石破碎后的动能确定方法 |
CN114113335A (zh) * | 2021-12-10 | 2022-03-01 | 东北大学 | 基于声发射/微震监测的岩石耗散能时空分布量化方法 |
CN114113335B (zh) * | 2021-12-10 | 2023-12-22 | 东北大学 | 基于声发射/微震监测的岩石耗散能时空分布量化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107101887B (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107101887A (zh) | 一种声发射与数值计算相结合的岩石破裂过程分析方法 | |
Reches et al. | Gouge formation by dynamic pulverization during earthquake rupture | |
Kong et al. | An experimental investigation on stress-induced cracking mechanisms of a volcanic rock | |
Yoon et al. | Numerical investigation on optimized stimulation of intact and naturally fractured deep geothermal reservoirs using hydro-mechanical coupled discrete particles joints model | |
He et al. | A novel experimental technique to simulate pillar burst in laboratory | |
Li et al. | Dynamic behavior of aluminum honeycomb sandwich panels under air blast: Experiment and numerical analysis | |
Mason et al. | Seismic soil–foundation–structure interaction observed in geotechnical centrifuge experiments | |
Kusui et al. | Mechanical behaviour of scaled-down unsupported tunnel walls in hard rock under high stress | |
Hu et al. | Rock damage and energy balance of strainbursts induced by low frequency seismic disturbance at high static stress | |
Skrzypkowski et al. | Adjustment of the yielding system of mechanical rock bolts for room and pillar mining method in stratified rock mass | |
Hu et al. | 3D continuum-discrete coupled modelling of triaxial Hopkinson bar tests on rock under multiaxial static-dynamic loads | |
Wang | Identification of early-warning key point for rockmass instability using acoustic emission/microseismic activity monitoring | |
Liu et al. | Numerical studies on the failure process and associated microseismicity in rock under triaxial compression | |
Sielicki et al. | Masonry wall behaviour under explosive loading | |
You et al. | Experimental and numerical investigation on the mechanical responses and cracking mechanism of 3D confined single-flawed rocks under dynamic loading | |
Zhang et al. | Fractal characteristics of acoustic emission events based on single-link cluster method during uniaxial loading of rock | |
CN106124311A (zh) | 基于应变测试判别裂纹扩展演化过程的方法 | |
Wu et al. | A reformulated rate-dependent visco-elastic model for dynamic deformation and fracture of PMMA with peridynamics | |
Xiao et al. | A stability study of goaf based on mechanical properties degradation of rock caused by rheological and disturbing loads | |
Panteki et al. | Finite element analysis of dynamic concrete-to-rebar bond experiments in the push-in configuration | |
Liu et al. | Investigation of impact rock burst induced by energy released from hard rock fractures | |
Dong et al. | Acoustic emission location accuracy and spatial evolution characteristics of granite fracture in complex stress conditions | |
KR20150091635A (ko) | 암반의 취성 파괴 가능성 평가 방법 | |
Zhang et al. | Damage evolution analysis of coal samples under cyclic loading based on single-link cluster method | |
CN101603285A (zh) | 一种跑道拦阻系统及其拦阻性能的预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 |
|
CF01 | Termination of patent right due to non-payment of annual fee |