CN113686646B - 一种确定岩石长期强度的方法 - Google Patents
一种确定岩石长期强度的方法 Download PDFInfo
- Publication number
- CN113686646B CN113686646B CN202110896551.0A CN202110896551A CN113686646B CN 113686646 B CN113686646 B CN 113686646B CN 202110896551 A CN202110896551 A CN 202110896551A CN 113686646 B CN113686646 B CN 113686646B
- Authority
- CN
- China
- Prior art keywords
- rock
- stress
- bulk modulus
- obtaining
- strain
- 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
Links
Images
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
-
- 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/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0071—Creep
-
- 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/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
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
本发明公开了一种确定岩石长期强度的方法,属于岩石力学性能技术领域。所述确定岩石长期强度的方法包括:获取岩石在不同应力条件下的体积应变;根据应力球张量和体积应变获取岩石的在不同应力条件下体积模量;取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图;根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式;根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值;取多个应力值的平均值,得到岩石长期强度。本发明确定岩石长期强度的方法加载施加短,可以确定岩石长期强度。
Description
技术领域
本发明涉及岩石力学性能技术领域,特别涉及一种确定岩石长期强度的方法。
背景技术
岩石在深部复杂地质环境下的流变性质是关乎深部围岩支护以及长期稳定性的关键力学性能,其中岩石的长期强度是所有深部地下工程都需要重视的力学性质之一。经过大量的试验和工程实例证明,地下岩体在长期荷载作用下变形破坏时的强度值是小于岩石自身的峰值强度值,且岩体的变形失稳也表现出了与时间有关的强度性能。同时,学者们还将岩石的长期强度定义为一种具有时效性的强度特性,将在地下工程试用期间不会发生变形破坏的最大作用荷载称为长期强度。因此,研究岩石的长期强度以及提出一种确定岩石长期强度的方法对实际工程具有重要的指导意义。
目前,工程中和学术研究中最常用确定岩石长期强度的方法一般为:
1,等时应力-应变曲线的拐点确定法,等时应力-应变曲线的拐点在实际中较难确定,该拐点的变化规律较为离散,导致对应的应力值无法明确确定;
2,过渡蠕变确定法,过渡蠕变确定法确定得到的岩石长期强度值是一个范围,没有确定的长期强度数值,该范围的大小也与制定蠕变方案的应力水平等级有关,故该方法具有一定的人为误差和局限性。
发明内容
本发明提供一种确定岩石长期强度的方法,解决了或部分解决了现有技术中岩石长期强度的应力值无法明确确定,存在误差的技术问题。
为解决上述技术问题,本发明提供了一种确定岩石长期强度的方法包括:获取岩石在不同力条件下的体积应变;假设岩石在加载瞬间产生的瞬时应变全部为弹性应变,根据应力球张量和体积应变获取岩石的在不同应力条件下体积模量;取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图;根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式;根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值;取多个应力值的平均值,得到岩石长期强度。
进一步地,所述获取岩石在不同应力条件下的体积应变包括:获取岩石在不同应力条件下的轴向蠕变历时曲线图,根据轴向蠕变历时曲线图获取岩石在不同应力条件下的轴向应变ε1;获取岩石在不同应力条件下的径向蠕变历时曲线图,根据径向蠕变历时曲线图获取岩石在不同应力条件下的径向应变ε3;根据εV=ε1+2ε3;式中:ε1为轴向应变;ε3为径向应变;εV为体积应变。
进一步地,根据岩石在不同应力条件下的体积应变,获取体积蠕变历时曲线。
进一步地,根据岩石的在不同应力条件下体积模量,获取岩石在不同应力作用下的体积模量变化规律图。
进一步地,所述获取不同时间的岩石的体积模量与应力之间的关系表达式包括:时间点取0h、4h、8h、12h、16h、20h、24h;则
1/K=-0.00101exp(σ1/21.12375)+0.01933,t=0h;
1/K=-0.001exp(σ1/20.67437)+0.01866,t=4h;
1/K=-0.00111exp(σ1/20.91978)+0.019,t=8h;
1/K=-0.00107exp(σ1/20.42317)+0.01874,t=12h;
1/K=-4.78497·10-4exp(σ1/17.02925)+0.01473,t=16h;
1/K=-0.00186exp(σ1/22.92788)+0.02298,t=20h;
1/K=-0.00185exp(σ1/22.77147)+0.0232,t=24h。
本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
由于采用了获取岩石在不同应力条件下的体积应变,假设岩石在加载瞬间产生的瞬时应变全部为弹性应变,根据应力球张量和体积应变获取岩石在不同应力条件下体积模量,取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图,根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式;根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值,取多个应力值的平均值,得到岩石长期强度,加载施加短,可以确定岩石长期强度,避免误差。
附图说明
图1为本发明实施例提供的确定岩石长期强度的方法的时间点为0h的倒数与应力之间的关系图;
图2为本发明实施例提供的确定岩石长期强度的方法的时间点为4h的倒数与应力之间的关系图
图3为本发明实施例提供的确定岩石长期强度的方法的时间点为8h的倒数与应力之间的关系图
图4为本发明实施例提供的确定岩石长期强度的方法的时间点为12h的倒数与应力之间的关系图
图5为本发明实施例提供的确定岩石长期强度的方法的时间点为16h的倒数与应力之间的关系图
图6为本发明实施例提供的确定岩石长期强度的方法的时间点为20h的倒数与应力之间的关系图
图7为本发明实施例提供的确定岩石长期强度的方法的时间点为24h的倒数与应力之间的关系图
图8为本发明实施例提供的确定岩石长期强度的方法的轴向蠕变历时曲线图
图9为本发明实施例提供的确定岩石长期强度的方法的径向蠕变历时曲线图
图10为本发明实施例提供的确定岩石长期强度的方法的体积蠕变历时曲线图;
图11为本发明实施例提供的确定岩石长期强度的方法的岩石50MPa、70MPa、80MPa及90MPa应力作用下的体积模量变化规律图;
图12为本发明实施例提供的确定岩石长期强度的方法的岩石60MPa应力作用下的体积模量变化规律图。
具体实施方式
参见图1-7,本发明实施例提供的一种确定岩石长期强度的方法包括:
步骤S1,获取岩石在不同应力条件下的体积应变。
步骤S2,假设岩石在加载瞬间产生的瞬时应变全部为弹性应变,根据应力球张量和体积应变获取岩石在不同应力条件下体积模量。
步骤S3,取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图。
步骤S4,根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式。
步骤S5,根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值。
步骤S6,取多个应力值的平均值,得到岩石长期强度。
本申请具体实施方式由于采用了获取岩石在不同应力条件下的体积应变,假设岩石在加载瞬间产生的瞬时应变全部为弹性应变,根据应力球张量和体积应变获取岩石的在不同应力条件下体积模量,取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图,根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式;根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值,取多个应力值的平均值,得到岩石长期强度,加载施加短,可以确定岩石长期强度,避免误差。
详细介绍步骤S1。
参见图8-9,所述获取岩石在不同应力条件下的体积应变包括:
步骤S11,获取岩石在不同应力条件下的轴向蠕变历时曲线图,根据轴向蠕变历时曲线图获取岩石在不同应力条件下的轴向应变ε1。
步骤S12,获取岩石在不同应力条件下的径向蠕变历时曲线图,根据径向蠕变历时曲线图获取岩石在不同应力条件下的径向应变ε3。
步骤S13,根据εV=ε1+2ε3;式中:ε1为轴向应变;ε3为径向应变;εV为体积应变。
参见图10,在本实施方式中,根据岩石在不同应力条件下的体积应变,获取体积蠕变历时曲线,可以直接由体积蠕变历时曲线获取不同应力条件下体积应变。
详细介绍步骤S2。
所述根据应力球张量和体积应变获取岩石的在不同应力条件下体积模量包括:
参见图11-12,在本实施方式中,根据岩石在不同应力条件下体积模量,获取岩石不同应力作用下的体积模量变化规律图。由图12可知,在应力为60MPa作用下体积模量开始由正值向负值转化,且岩石的体积蠕变应变也开始由正值向负值转化,即岩石体积由压缩状态向扩容状态转化。因此,可以认为岩石的长期强度值的范围值50-70MPa。
参见图1-7,详细介绍步骤S4,。
参见图1-7,所述获取不同时间的岩石的体积模量与应力之间的关系表达式包括:
时间点取0h、4h、8h、12h、16h、20h、24h;则
1/K=-0.00101exp(σ1/21.12375)+0.01933,t=0h;
1/K=-0.001exp(σ1/20.67437)+0.01866,t=4h;
1/K=-0.00111exp(σ1/20.91978)+0.019,t=8h;
1/K=-0.00107exp(σ1/20.42317)+0.01874,t=12h;
1/K=-4.78497·10-4exp(σ1/17.02925)+0.01473,t=16h;
1/K=-0.00186exp(σ1/22.92788)+0.02298,t=20h;
1/K=-0.00185exp(σ1/22.77147)+0.0232,t=24h。
令体积模量K=0,可以得到体积模量由正值向负值转化时对应的应力值。时间点取0、4、8、12、16、20和24h时,所对应的应力值分别为62.35、60.50、59.41、58.47、58.36、57.64、57.59MPa,将上述所有应力取平均值即可得到岩石的长期强度值为59.19MPa。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种确定岩石长期强度的方法,其特征在于,包括:
获取岩石在不同应力条件下的轴向蠕变历时曲线图,根据轴向蠕变历时曲线图获取岩石在不同应力条件下的轴向应变ε1;
获取岩石在不同应力条件下的径向蠕变历时曲线图,根据径向蠕变历时曲线图获取岩石在不同应力条件下的径向应变ε3;
根据εV=ε1+2ε3获取岩石在不同应力条件下的体积应变;式中:ε1为轴向应变;ε3为径向应变;εV为体积应变;
根据岩石在不同应力条件下的体积应变,获取体积蠕变历时曲线;
假设岩石在加载瞬间产生的瞬时应变全部为弹性应变,根据应力球张量和体积应变获取岩石的在不同应力条件下体积模量;
根据岩石的在不同应力条件下体积模量,获取岩石不同应力作用下的体积模量变化规律图;
取不同时间点的岩石的体积模量的倒数,获得体积模量的倒数与应力之间的关系图;
根据体积模量的倒数与应力之间的关系图,获取不同时间的岩石的体积模量与应力之间的关系表达式;
根据不同时间的岩石的体积模量与应力之间的关系表达式,且视岩石的体积模量为0,获取体积模量由正值向负值转化时对应的多个应力值;
取多个应力值的平均值,得到岩石长期强度。
4.根据权利要求2所述的确定岩石长期强度的方法,其特征在于,所述获取不同时间的岩石的体积模量与应力之间的关系表达式包括:
时间点取0h、4h、8h、12h、16h、20h、24h;则
1/K=-0.00101exp(σ1/21.12375)+0.01933,t=0h;
1/K=-0.001exp(σ1/20.67437)+0.01866,t=4h;
1/K=-0.00111exp(σ1/20.91978)+0.019,t=8h;
1/K=-0.00107exp(σ1/20.42317)+0.01874,t=12h;
1/K=-4.78497·10-4exp(σ1/17.02925)+0.01473,t=16h;
1/K=-0.00186exp(σ1/22.92788)+0.02298,t=20h;
1/K=-0.00185exp(σ1/22.77147)+0.0232,t=24h。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896551.0A CN113686646B (zh) | 2021-08-05 | 2021-08-05 | 一种确定岩石长期强度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896551.0A CN113686646B (zh) | 2021-08-05 | 2021-08-05 | 一种确定岩石长期强度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113686646A CN113686646A (zh) | 2021-11-23 |
CN113686646B true CN113686646B (zh) | 2022-07-12 |
Family
ID=78578919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110896551.0A Active CN113686646B (zh) | 2021-08-05 | 2021-08-05 | 一种确定岩石长期强度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113686646B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2339816C1 (ru) * | 2007-02-21 | 2008-11-27 | Государственное образовательное учреждение высшего профессионального образования Московский государственный горный университет (МГГУ) | Способ определения длительной прочности горных пород |
CN104849134A (zh) * | 2015-05-26 | 2015-08-19 | 中国石油大学(华东) | 多级应力分级加载蠕变力学试验确定岩石长期强度的方法 |
CN107977521A (zh) * | 2017-12-11 | 2018-05-01 | 安徽理工大学 | 一种延期突出煤的蠕变本构及其数值化实现技术 |
CN111678787A (zh) * | 2020-06-16 | 2020-09-18 | 中国矿业大学 | 一种基于数值拟合的岩石长期强度确定方法 |
CN111861107A (zh) * | 2020-06-08 | 2020-10-30 | 中国地质大学(武汉) | 一种基于强度折减法的滑坡动态稳定性评价方法 |
-
2021
- 2021-08-05 CN CN202110896551.0A patent/CN113686646B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2339816C1 (ru) * | 2007-02-21 | 2008-11-27 | Государственное образовательное учреждение высшего профессионального образования Московский государственный горный университет (МГГУ) | Способ определения длительной прочности горных пород |
CN104849134A (zh) * | 2015-05-26 | 2015-08-19 | 中国石油大学(华东) | 多级应力分级加载蠕变力学试验确定岩石长期强度的方法 |
CN107977521A (zh) * | 2017-12-11 | 2018-05-01 | 安徽理工大学 | 一种延期突出煤的蠕变本构及其数值化实现技术 |
CN111861107A (zh) * | 2020-06-08 | 2020-10-30 | 中国地质大学(武汉) | 一种基于强度折减法的滑坡动态稳定性评价方法 |
CN111678787A (zh) * | 2020-06-16 | 2020-09-18 | 中国矿业大学 | 一种基于数值拟合的岩石长期强度确定方法 |
Non-Patent Citations (4)
Title |
---|
基于应力和时间双重影响下岩石蠕变模型研究;刘文博 等;《中南大学学报(自然科学版)》;20200831;第51卷(第08期);第2256-2265页 * |
基于流变试验的向家坝砂岩长期强度评价;李良权 等;《工程力学》;20101130;第27卷(第11期);第127-136和143页 * |
基于西原模型的变参数蠕变方程;阎岩 等;《岩土力学》;20101031;第31卷(第10期);第3025-3035页 * |
岩石长期强度确定方法介绍;齐亚静 等;《水利科技与经济》;20160630;第22卷(第06期);第18-22页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113686646A (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110705165A (zh) | 一种构建岩石材料弹塑性-损伤耦合力学本构模型的方法 | |
CN110296892B (zh) | 基于能量分析的岩石损伤演化过程中特征应力确定方法 | |
Chow et al. | On evolution laws of anisotropic damage | |
Bacca et al. | A model for the Mullins effect in multinetwork elastomers | |
CN113686646B (zh) | 一种确定岩石长期强度的方法 | |
CN112434473A (zh) | 一种考虑损伤渗流应力耦合的数值模拟方法 | |
CN114894616A (zh) | 一种基于变形模量的岩石蠕变模型粘、弹性参数获取方法 | |
CN115795916A (zh) | 一种岩石脆性评价方法及装置 | |
Xia et al. | Strength criterion for rocks subjected to cyclic stress and temperature variations | |
CN113607547A (zh) | 一种基于岩石裂隙初始点的弹性应变能获取方法 | |
Wang et al. | Lateral response and installation recommendation of hybrid monopile foundation for offshore wind turbines under combined loadings | |
CN112051141A (zh) | 一种复合材料应变率相关压缩本构模型的构建方法 | |
CN112730057B (zh) | 一种压缩荷载作用脆性固体材料长期寿命评价方法 | |
Xie et al. | Numerical and experimental study on rod-fastened rotor dynamics using semi-analytical elastic-plastic model | |
CN116562055B (zh) | 一种加卸载环境中锚索支护锚固性能确定方法 | |
Ming-xing et al. | Transfer matrix solutions for lateral behavior of pile foundation considering the skin and end resistance effect | |
DU et al. | Experimental study on energy response of argillaceous sandstone under cyclic loading | |
Pei et al. | Particle bond-degradation model for cumulative damage analysis of concrete under cyclic loading | |
He et al. | Analytical and numerical solutions for shear mechanical behaviors of structural plane | |
Wang et al. | Experimental investigation and constitutive modelling of creep-damage behaviours in monzogranite | |
CN114398805A (zh) | 一种水岩耦合作用下裂隙岩石蠕变模型构建方法及系统 | |
Sabelkin et al. | Investigation into relative slip during fretting fatigue under partial slip contact condition | |
Shen et al. | A method of analyzing laterally loaded monopiles considering the influence of additional resistance | |
CN112834319A (zh) | 一种软土小应变硬化模型力学参数确定方法 | |
CN117213956A (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 |