CN112198052B - 岩石在围压条件下的拉伸强度计算方法及应用 - Google Patents

岩石在围压条件下的拉伸强度计算方法及应用 Download PDF

Info

Publication number
CN112198052B
CN112198052B CN202011046352.2A CN202011046352A CN112198052B CN 112198052 B CN112198052 B CN 112198052B CN 202011046352 A CN202011046352 A CN 202011046352A CN 112198052 B CN112198052 B CN 112198052B
Authority
CN
China
Prior art keywords
rock
tensile strength
confining pressure
tensile
rock sample
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.)
Active
Application number
CN202011046352.2A
Other languages
English (en)
Other versions
CN112198052A (zh
Inventor
黄达
刘洋
岑夺丰
宫凤强
钟助
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hebei University of Technology
Original Assignee
Hebei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN202011046352.2A priority Critical patent/CN112198052B/zh
Publication of CN112198052A publication Critical patent/CN112198052A/zh
Application granted granted Critical
Publication of CN112198052B publication Critical patent/CN112198052B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001Type of application of the stress
    • G01N2203/0003Steady
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0017Tensile
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/006Crack, flaws, fracture or rupture
    • G01N2203/0067Fracture or rupture
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/025Geometry of the test
    • G01N2203/0256Triaxial, i.e. the forces being applied along three normal axes of the specimen
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property
    • G01N2203/0676Force, weight, load, energy, speed or acceleration

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)

Abstract

本发明公开了一种岩石在围压条件下的拉伸强度计算方法,该方法根据岩石围压与拉伸强度的关系准则计算岩石的拉伸强度。该方法包括步骤S1、进行岩石试样在不同围压条件下的拉伸试验,获取岩石试样在不同围压作用下的拉伸强度;步骤S2、根据不同围压及其对应的拉伸强度进行数据拟合,确定当前岩石材料属性相关的常数值,再根据确定的岩石围压与拉伸强度的关系准则,计算出当前岩石在不同围压条件下的拉伸强度。该方法提出的岩石围压与拉伸强度的关系准则为岩石破裂识别提供一种依据,为通过拉伸强度判定岩石破裂提供一个对比基准,同时为工程岩体的稳定性控制和失稳预警提供理论依据。

Description

岩石在围压条件下的拉伸强度计算方法及应用
技术领域
本发明涉及岩石力学研究技术领域,具体是一种岩石在围压条件下的拉伸强度计算方法及应用。
背景技术
在绝大多数岩石工程的建设施工中经常会遇到岩石拉伸破坏,如隧道开挖、地下洞室施工、高陡边坡卸荷以及采矿等。岩石破裂过程中往往伴随着较大的地应力,也就是岩石会受到围压,而在复杂地应力条件下的拉伸破坏往往更加剧烈,一旦处理不当,很可能造成施工进度延缓以及对人员设备等造成损害。目前普遍将岩石破裂过程等效为单纯的单向拉伸破坏,而忽略岩石在实际环境中会受到围压作用,因此有必要对围压条件下的岩石拉伸强度进行深入研究。
目前对于岩石在不同围压条件下拉伸强度的研究还比较少,更没有一个通用性的围压与拉伸强度之间的关系准则,无法对在围压作用下的岩石拉伸破坏强度变化进行预测。
发明内容
针对现有技术的不足,本发明拟解决的技术问题是,提供一种岩石在围压条件下的拉伸强度计算方法以应用。
本发明解决所述技术问题采用的技术方案是:
一种岩石在围压条件下的拉伸强度计算方法,其特征在于,该方法根据公式(1)的岩石围压与拉伸强度的关系准则计算岩石的拉伸强度;
Figure BDA0002708125540000011
其中,σ3为拉伸强度,σ1为围压值;m、n、k、q和a为常数。
该方法包括以下步骤:
步骤S1、进行岩石试样在不同围压条件下的拉伸试验,获取岩石试样在不同围压作用下的拉伸强度;
步骤S2、根据不同围压及其对应的拉伸强度进行数据拟合,确定公式(1)中的当前岩石材料属性相关的常数值,再根据确定的岩石围压与拉伸强度的关系准则,计算出当前岩石在不同围压条件下的拉伸强度。
步骤S1中拉伸试验的过程是:
1-1制备岩石试样,岩石试样的两端分别粘结带螺纹孔的金属端头;
1-2将岩石试样竖直安装在拉伸-压缩转换装置上,将岩石试样两端的金属端头与拉伸-压缩转换装置上、下部的连接端固连;然后将安装有岩石试样的拉伸-压缩转换装置整体密封在三轴压缩试验机的围压腔中,三轴压缩试验机的两个加载头分别与拉伸-压缩转换装置的加载端固连,向围压腔内充入硅油;
1-3对岩石试样施加围压,待围压值稳定后对岩石试样逐渐施加轴向压力;当岩石试样发生断裂,停止拉伸。
岩石试样采用花岗岩制成,岩石试样为高100mm,直径为50mm的圆柱体。
本发明还提供一种岩石在围压条件下的拉伸强度计算方法的应用,其特征在于,需要预先进行岩石试样的拉伸试验,确定不同岩石材料对应的m、n、k、q和a的取值;测量某处岩石的围压值,根据公式(1)的关系准则计算得到该处岩石破裂时拉伸强度的理论值;
Figure BDA0002708125540000021
实时监测该处岩石的实际强度,当岩石的实际强度小于拉伸强度的理论值时表示岩石未发生破裂,可以继续开挖;当岩石的实际强度大于或等于拉伸强度的理论值时表示岩石发生破裂,需要停止开挖操作。
与现有技术相比,本发明具有以下有益效果:
1.本发明提出的岩石围压与拉伸强度的关系准则为岩石破裂识别提供一种依据,为通过拉伸强度判定岩石破裂提供一个对比基准,同时为工程岩体的稳定性控制和失稳预警提供理论依据;在实际工程中,将拉伸强度的理论值与测量值进行对比,即可判断岩石是否发生破裂,同时可以对岩石的拉伸强度进行预测,方法简单可行。
2.本发明提出的一种岩石在围压条件下的拉伸强度计算方法,是一种全新的岩石拉伸强度确定方法,其采用分段函数的方式对岩石的拉伸强度进行确定,很好地揭示了岩石在不同围压下拉伸强度的演化规律;通过此方法能够更加直观的了解岩石的拉伸力学行为,认识岩石的拉伸破坏特征,为地下工程、边坡开挖和采矿等引起的岩石拉伸破坏失稳提供预警和防控建议。
附图说明
图1为本发明的整体流程示意图;
图2为本发明的实施例中花岗岩所受围压与拉伸强度的关系曲线图。
具体实施方式
为使本领域技术人员更好地理解本发明的技术方案,下面结合附图和具体实施例对本发明作进一步的详细说明。
本发明提供了一种岩石在围压条件下的拉伸强度计算方法(简称方法,参见图1-2),包括以下步骤:
步骤S1、进行岩石试样在不同围压条件下的拉伸试验,获取岩石试样在不同围压作用下的拉伸强度;
步骤S2、根据不同围压及其对应的拉伸强度进行数据拟合,建立如公式(1)的岩石围压与拉伸强度的关系准则,根据准则计算岩石的拉伸强度;
Figure BDA0002708125540000022
其中,σ3为拉伸强度,σ1为围压值;m、n、k、q和a为与岩石材料属性相关的常数,不同岩石材料这些常数的取值不同。
实施例
本实施例提供一种岩石在围压条件下的拉伸强度计算方法,包括以下步骤:
步骤S1、进行岩石试样在不同围压条件下的拉伸试验,获取岩石试样在不同围压作用下的拉伸强度;
1-1采用花岗岩制备岩石试样,将岩石试样加工为高100mm,直径为50mm的圆柱体;试样加工精度和表面平整度符合国标;岩石试样的两端通过强力胶(9900丙烯酸结构胶)分别粘结带螺纹孔的金属端头,并放置48小时,以达到粘结强度;
1-2将岩石试样竖直安装在拉伸-压缩转换装置上,分别通过螺栓将岩石试样两端的金属端头与拉伸-压缩转换装置上、下部的连接端固连;然后将安装有岩石试样的拉伸-压缩转换装置整体密封在三轴压缩试验机的围压腔中,三轴压缩试验机的两个加载头分别与拉伸-压缩转换装置的加载端固连,向围压腔内充入硅油,通过挤压硅油对岩石试样施加围压;
拉伸-压缩转换装置的结构及原理见专利CN201611055843.7,其作用是将三轴压缩试验机的压应力转换为拉伸应力;三轴压缩试验机可采用TFD-2000L微机伺服控制岩石三轴流变试验机,该试验机轴向最大试验力可达2000kN、最大侧压100Mpa、温度控制范围为室温至200℃;
1-3对岩石试样施加围压,待围压值稳定后对岩石试样逐渐施加轴向压力,拉伸-压缩转换装置将轴向压力转换为拉伸应力,实现岩石试样的拉伸试验;当岩石试样发生断裂,停止拉伸,三轴压缩试验机测量并记录此时的拉伸强度;改变围压值,再对岩石试样进行拉伸试验,重复试验多次,三轴压缩试验机测量并记录多组围压值及其对应的拉伸强度,结果如表1所示;
表1不同围压下的拉伸强度
围压(MPa) 拉伸强度(MPa)
0 -7.68
3 -8.71
6 -8.99
9 -8.40
12 -6.79
15 -7.41
18 -6.34
20 -4.34
21 -4.20
24 -2.31
28 -1.21
步骤S2、根据不同围压及其对应的拉伸强度进行数据拟合,建立如公式(1)的岩石围压与拉伸强度的关系准则;根据围压与拉伸强度的关系准则计算岩石的拉伸强度;
将表1中的试验数据进行拟合,得到花岗岩的拉伸强度与围压的关系曲线(如图2所示),通过分析发现,花岗岩的拉伸强度与围压的关系曲线由非常明显的两部分组成:即围压水平较小时,花岗岩的拉伸强度随围压增大呈现出二次函数曲线变化趋势;当围压水平较大时,花岗岩的拉伸强度随围压增大呈现出一次函数线性变化趋势,两段曲线在连接点处光滑连续且可导;
Figure BDA0002708125540000041
根据图2中的关系曲线推导出如公式(1)的围压与拉伸强度的关系准则,σ3为拉伸强度,σ1为围压值,此时m取0.0213,n取-0.2398,k取-7.938,q取-13.936,a取16.77。
应用岩石围压与拉伸强度的关系准则时,需要预先根据步骤S1和S2的操作进行岩石试样的拉伸试验,确定不同岩石材料对应的m、n、k、q和a的取值;该关系准则可以用于防止岩石垮塌,比如隧道开挖工程中,可以通过测量得到某处岩石的围压值,根据公式(1)的关系准则计算得到该处岩石破裂时拉伸强度的理论值;通过水压致裂法等现有方法实时监测该处岩石的实际强度,当岩石的实际强度小于拉伸强度的理论值时表示岩石未发生破裂,可以继续开挖;当岩石的实际强度大于或等于拉伸强度的理论值时表示岩石发生破裂,需要停止开挖操作。该关系准则还能用于实验室中对微小围压或者特别大围压条件下的岩石拉伸强度进行预测,由于微小围压不易控制,特大围压实施特别危险,均不容易实施,利用该关系准则就可以不需要进行试验直接获得这两种情况下的岩石拉伸强度。
本申请分别对砂岩、页岩和大理岩进行相同的操作,得出这三种岩石的拉伸强度与围压的关系曲线均满足:围压水平较小时,岩石的拉伸强度随围压增大呈现出二次函数曲线变化趋势;当围压水平较大时,岩石的拉伸强度随围压增大呈现出一次函数线性变化趋势,两段曲线在连接点处光滑连续且可导,也就是不同材料的岩石均存在这种关系准则。
本发明未述及之处适用于现有技术。

Claims (5)

1.一种岩石在围压条件下的拉伸强度计算方法,其特征在于,该方法根据公式(1)的岩石围压与拉伸强度的关系准则计算岩石的拉伸强度;
Figure FDA0002708125530000011
其中,σ3为拉伸强度,σ1为围压值;m、n、k、q和a为常数。
2.根据权利要求1所述的岩石在围压条件下的拉伸强度计算方法,其特征在于,该方法包括以下步骤:
步骤S1、进行岩石试样在不同围压条件下的拉伸试验,获取岩石试样在不同围压作用下的拉伸强度;
步骤S2、根据不同围压及其对应的拉伸强度进行数据拟合,确定公式(1)中的当前岩石材料属性相关的常数值,再根据确定的岩石围压与拉伸强度的关系准则,计算出当前岩石在不同围压条件下的拉伸强度。
3.根据权利要求2所述的岩石在围压条件下的拉伸强度计算方法,其特征在于,步骤S1中拉伸试验的过程是:
1-1制备岩石试样,岩石试样的两端分别粘结带螺纹孔的金属端头;
1-2将岩石试样竖直安装在拉伸-压缩转换装置上,将岩石试样两端的金属端头与拉伸-压缩转换装置上、下部的连接端固连;然后将安装有岩石试样的拉伸-压缩转换装置整体密封在三轴压缩试验机的围压腔中,三轴压缩试验机的两个加载头分别与拉伸-压缩转换装置的加载端固连,向围压腔内充入硅油;
1-3对岩石试样施加围压,待围压值稳定后对岩石试样逐渐施加轴向压力;当岩石试样发生断裂,停止拉伸。
4.根据权利要求3所述的岩石在围压条件下的拉伸强度计算方法,其特征在于,岩石试样采用花岗岩制成,岩石试样为高100mm,直径为50mm的圆柱体。
5.一种岩石在围压条件下的拉伸强度计算方法的应用,其特征在于,需要预先进行岩石试样的拉伸试验,确定不同岩石材料对应的m、n、k、q和a的取值;测量某处岩石的围压值,根据公式(1)的关系准则计算得到该处岩石破裂时拉伸强度的理论值;
Figure FDA0002708125530000012
实时监测该处岩石的实际强度,当岩石的实际强度小于拉伸强度的理论值时表示岩石未发生破裂,可以继续开挖;当岩石的实际强度大于或等于拉伸强度的理论值时表示岩石发生破裂,需要停止开挖操作。
CN202011046352.2A 2020-09-29 2020-09-29 岩石在围压条件下的拉伸强度计算方法及应用 Active CN112198052B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011046352.2A CN112198052B (zh) 2020-09-29 2020-09-29 岩石在围压条件下的拉伸强度计算方法及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011046352.2A CN112198052B (zh) 2020-09-29 2020-09-29 岩石在围压条件下的拉伸强度计算方法及应用

Publications (2)

Publication Number Publication Date
CN112198052A CN112198052A (zh) 2021-01-08
CN112198052B true CN112198052B (zh) 2022-07-12

Family

ID=74007035

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011046352.2A Active CN112198052B (zh) 2020-09-29 2020-09-29 岩石在围压条件下的拉伸强度计算方法及应用

Country Status (1)

Country Link
CN (1) CN112198052B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112525671B (zh) * 2021-01-21 2021-07-30 东北石油大学 一种真三轴围压条件下岩石直接拉伸试验装置
CN113008307B (zh) * 2021-03-09 2021-11-30 华侨大学 一种确定等间距破裂岩层裂隙间距的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106556536A (zh) * 2016-11-14 2017-04-05 重庆大学 一种可在三轴压缩试验机上使用的岩石三轴拉压试验装置
CN107463740A (zh) * 2017-07-27 2017-12-12 中南大学 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法
CN107991176A (zh) * 2018-01-30 2018-05-04 中南大学 一种岩石三轴的拉伸试验装置及其方法
CN109443943A (zh) * 2018-11-26 2019-03-08 长江大学 岩石样本强度准则材料参数的构建方法及设备

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201196615Y (zh) * 2007-09-26 2009-02-18 河南理工大学 围压下岩样直接拉伸试验装置
US8783091B2 (en) * 2009-10-28 2014-07-22 Halliburton Energy Services, Inc. Cement testing
US20150152724A1 (en) * 2013-06-19 2015-06-04 Conocophillips Company Core sample testing protocol
US9746403B2 (en) * 2014-10-06 2017-08-29 CNPC USA Corp. Method of testing a polycrystalline diamond compact cutter
CN104849134B (zh) * 2015-05-26 2016-01-20 中国石油大学(华东) 多级应力分级加载蠕变力学试验确定岩石长期强度的方法
US11578564B2 (en) * 2018-05-30 2023-02-14 Saudi Arabian Oil Company Systems and methods for predicting shear failure of a rock formation
CN109855973A (zh) * 2019-01-22 2019-06-07 东北大学 一种岩石三轴直接拉伸室内实验装置和方法
CN110031304B (zh) * 2019-04-25 2020-04-21 四川大学 一种岩石变形破坏模式预测方法
CN110470529A (zh) * 2019-09-19 2019-11-19 中国科学院武汉岩土力学研究所 三维应力状态下岩石间接拉伸强度测试方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106556536A (zh) * 2016-11-14 2017-04-05 重庆大学 一种可在三轴压缩试验机上使用的岩石三轴拉压试验装置
CN107463740A (zh) * 2017-07-27 2017-12-12 中南大学 考虑中间主应力效应的岩石类材料真三轴试验数值模拟方法
CN107991176A (zh) * 2018-01-30 2018-05-04 中南大学 一种岩石三轴的拉伸试验装置及其方法
CN109443943A (zh) * 2018-11-26 2019-03-08 长江大学 岩石样本强度准则材料参数的构建方法及设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A strength criterion based on critical state mechanics for intact rocks;Singh M;《Rock Mechanics and Rock Engineering》;20051230(第3期);全文 *
Practical estimates of rock mass strength;Hoek E;《International Journal of Rock》;19971230(第8期);全文 *
Strain Rate Dependency of Coarse Crystal Marble Under Uniaxial Compression: Strength, Deformation and Strain Energy;Yanrong Li;《Rock Mechanics and Rock Engineering 》;20141230(第4期);全文 *
高应力卸荷条件下大理岩破裂面细微观形态特征及其与卸荷岩体强度的相关性研究;黄达;《岩土力学》;20121230;全文 *

Also Published As

Publication number Publication date
CN112198052A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
CN112198052B (zh) 岩石在围压条件下的拉伸强度计算方法及应用
US10928282B2 (en) Method and device for determining elasticity of cement stone utilized in well cementing of oil-gas well
CN110595909B (zh) 模拟深部岩体不同温度影响下的真三轴试验系统及方法
Lang et al. Creep experiment and rheological model of deep saturated rock
CN102221504B (zh) 一种测定破裂后岩石三轴强度的加载路径方法
CN106124313B (zh) 混凝土及类似材料受压下综合变形性能的测试装置
CN111238931B (zh) 基于能量演化的页岩脆性指数评价方法
CN103163023B (zh) 一种石油管材环境断裂韧性测试方法及装置
CN205103090U (zh) 一种岩石力学三轴压力试验机的压力室结构
CN106769456A (zh) 一种持久荷载下全级配混凝土长期性能测试装置及方法
CN107907409A (zh) 一种确定岩石起裂应力的方法、设备及存储设备
CN101509852A (zh) 一种获取厚壁圆筒试样环形断裂的试验方法
CN108896415B (zh) 一种岩体蠕变冲击试验方法
CN116879408A (zh) 一种基于声发射b值特性指标的岩爆前兆获取方法
CN109030242B (zh) 一种电磁动力岩石直剪仪及操作方法
CN108519264B (zh) 一种裂隙试样制备方法
CN102230869A (zh) 现场岩体流变试验仪及其测试方法
Yu, X.*, Da Gama, CD**, Na, Y.***, Wang, Q**** & Xie Deformation behaviour of rocks under compression and direct tension
CN115931606A (zh) 基于循环加卸载试验作用下岩石损伤程度计算的方法
CN112198051B (zh) 侧压作用下基于能量演化的岩石拉伸破裂识别方法及应用
CN111855412B (zh) 一种基于应力能量比的岩爆倾向性等级判别方法
CN114965055A (zh) 利用圆环变形法测量岩石拉伸应变及其力学参数的方法
Brenne et al. Hydraulic Fracturing of a Devonian Slate under confining pressure–with emphasis on cleavage inclination
CN113125266A (zh) 一种岩石粘聚力和内摩擦角时效劣化演化方程获取方法
Gao et al. Crack closure and initiation stresses of coal subjected to thermo-gas-mechanical coupling

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