CN113239585A - 测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 - Google Patents
测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 Download PDFInfo
- Publication number
- CN113239585A CN113239585A CN202110498290.7A CN202110498290A CN113239585A CN 113239585 A CN113239585 A CN 113239585A CN 202110498290 A CN202110498290 A CN 202110498290A CN 113239585 A CN113239585 A CN 113239585A
- Authority
- CN
- China
- Prior art keywords
- sample
- asymmetric
- fracture toughness
- rock
- grooving
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种测量岩石Ⅱ型断裂韧度的非对称切槽短棒方法,通过制作在岩石两侧引入平行于试样端面的半通切槽的压缩短芯试样。利用有限元分析软件研究静压力加载条件下试样的最佳几何形状。进而通过数值模拟确定试样的应力强度因子并由此推导出断裂韧度的公式。本发明优点是:简化了试样的制备并使用岩芯作为基础的试样,便于静态和动态试验的进行;并在数值模拟时考虑了次生裂纹的存在。
Description
技术领域
本发明涉及岩土工程岩石断裂技术领域,特别涉及一种用于测量岩石Ⅱ型断裂韧度的非对称切槽短棒的方法。
背景技术
土木工程和采矿结构的破坏对工程活动有重大影响。因此,分析岩石和类岩石材料的破坏机理至关重要,因为这有助于工程师提高安全性和防止采矿中的破坏。理解岩石破坏机理的理论方法之一是断裂力学,它规定了岩石在静态或动态载荷下的断裂条件。岩石断裂力学研究的是岩石材料中裂纹的萌生和扩展。根据加载类型,有三种基本断裂模式:Ⅰ型(张开模式)、Ⅱ型(滑动模式或面内剪切模式)、Ⅲ型(滑开模式或面外剪切模式)。
作为断裂参数之一,断裂韧度即临界应力强度因子(SIF),代表材料抵抗裂纹扩展的能力。断裂韧度是材料最重要的力学性能之一,每种断裂模式都有一个断裂韧度。前人已提出多种Ⅱ型断裂韧度的测量方法,但是已有的方法需要复杂的试件夹具和大量的加工工作来制备试件,而且这些方法产生的裂缝轨迹偏离了原始的切槽面,导致人们对真正的裂缝起裂模式产生怀疑。在断裂力学中,断裂韧度是抗裂纹扩展的度量,计算断裂韧度需要一个初始裂纹。但目前的多种测量方法往往忽略了这个前提,因此,计算断裂韧度方法的可靠性不足。
发明内容
本发明针对现有技术的缺陷,提供了一种测量岩石Ⅱ型断裂韧度的非对称切槽短棒方法,解决了现有技术中存在的缺陷。
为了实现以上发明目的,本发明采取的技术方案如下:
一种测量岩石Ⅱ型断裂韧度的非对称切槽短棒方法,包括以下步骤:
步骤1:制备压缩短芯试样;
步骤2:利用有限元分析软件研究静压力加载条件下试样三维模型的应力场来确定最大剪应力所在平面;
步骤3:由上述平面建立2D模型来确定最佳C/H值;
步骤4:使用J积分法的有限元程序确定无量纲的应力强度因子;
步骤5:由步骤4的结果推导出断裂韧度的公式:
其中:其中:C为两个槽口之间的距离;
H为试样的长度;
P是施加在试件端面上的荷载;
D为岩芯直径;
a为切槽的深度,为D/2。
进一步地,所述压缩短芯试样是通过从相反的侧面引入两个平行于核心端部的半通槽口来制作的,切槽间隙不超过1mm;缺口和其最近的试样端部之间的距离相同,切槽前端平行,在穿过试样中心轴线的中间平面形成矩形岩桥。
进一步地,压缩短芯试样非对称切槽间距为试样厚度的1/10~1/5,C/H值,非对称切槽中心线与试样厚度中心线重合。
作为优选,所述测量岩石Ⅱ型静态断裂韧度的非对称切槽短棒形状为长径比为2:1的圆柱试样,例如:直径50mm、长度100mm。
作为优选,所述测量岩石Ⅱ型动态断裂韧度的非对称切槽短棒形状为长径比为1:1的圆柱试样,例如:直径38mm、长度38mm。
与现有技术相比,本发明的优点在于:
简化了试样的制备并使用岩芯作为基础的试样;并在数值模拟时考虑了次生裂纹的存在,而次生裂纹也是断裂力学理论计算应力强度因子的先决条件。
附图说明
图1是本发明实施例非对称切槽短棒方法流程图;
图2是本发明实施例压缩短芯试样结构示意图。
图3是本发明实施例实验数据对比图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
如图1所示,一种测量岩石Ⅱ型断裂韧度的非对称切槽短棒方法,包括以下步骤:
(1)如图2所示,制备半通切槽的压缩短芯试样(切槽间隙不超过1mm)。
(2)利用有限元分析软件研究静压力加载条件下试样三维模型的应力场来确定最大剪应力所在平面。具体过程如下:
使用有限元程序建立试样的三维(3D)模型,压缩载荷P固定。沿上缺口提取剪切应力的节点,可得到剪应力的最大值位于中心平面(垂直于破坏面)。
(3)建立中心平面的2D压缩短芯试样模型,以研究C/H值对剪应力分布的影响。对应力场进行分析后发现,对于压缩短芯试样,当H/D=1、2时,C/H值应小于0.3,以确保剪切破坏是沿断裂面进行的。
(4)通过在压缩短芯试件的每个缺口尖端引入次生裂纹。使用J积分法的有限元程序确定无量纲的应力强度因子。结果表明,次生裂纹深度hc在0.2mm至1.0mm范围内,应力强度因子作为hc的函数近似线性增加:
KⅡ=1.474hc+7.4954
(5)结合上式推导出断裂韧度的公式:
其中:C为两个槽口之间的距离;
H为试样的长度;
P是施加在试件端面上的荷载;
D为岩芯直径;
a为切槽的深度,为D/2。
该方法通过制作在岩石两侧引入平行于试样端面的半通切槽的压缩短芯试样。利用有限元分析软件研究静压力加载条件下试样的最佳几何形状。根据断裂力学理论计算应力强度因子的前提条件,考虑沿剪切断裂的次生裂纹的存在。进而通过数值模拟确定试样的应力强度因子并由此推导出断裂韧度的公式。
与PTS方法对比。
PTS方法已被接受为ISRM建议的测定岩石Ⅱ型静态模式断裂韧度的方法。
通过实验比对两者接近(如图3所示)。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (5)
2.根据权利要求书1所述的非对称切槽短棒方法,其特征在于:所述压缩短芯试样是通过从相反的侧面引入两个平行于核心端部的半通槽口来制作的,切槽间隙不超过1mm;缺口和其最近的试样端部之间的距离相同,切槽前端平行,在穿过试样中心轴线的中间平面形成矩形岩桥。
3.根据权利要求书2所述的非对称切槽短棒方法,其特征在于:压缩短芯试样非对称切槽间距为试样厚度的1/10~1/5,C/H值,非对称切槽中心线与试样厚度中心线重合。
4.根据权利要求书3所述的非对称切槽短棒方法,其特征在于:测量岩石Ⅱ型静态断裂韧度的非对称切槽短棒形状为长径比为2:1的圆柱试样。
5.根据权利要求书3所述的非对称切槽短棒方法,其特征在于:测量岩石Ⅱ型动态断裂韧度的非对称切槽短棒形状为长径比为1:1的圆柱试样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498290.7A CN113239585A (zh) | 2021-05-08 | 2021-05-08 | 测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498290.7A CN113239585A (zh) | 2021-05-08 | 2021-05-08 | 测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113239585A true CN113239585A (zh) | 2021-08-10 |
Family
ID=77132464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110498290.7A Withdrawn CN113239585A (zh) | 2021-05-08 | 2021-05-08 | 测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113239585A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4299120A (en) * | 1979-03-19 | 1981-11-10 | Terra Tek, Inc. | Method for determining plane strain fracture toughness of non-elastic fracture mechanics specimens |
CN104833590A (zh) * | 2014-07-08 | 2015-08-12 | 吴礼舟 | 岩石ii型断裂韧度测试的新方法 |
CN105043849A (zh) * | 2015-07-08 | 2015-11-11 | 四川大学 | 用于岩石i-ii复合断裂韧度测试的深梁试件及复合断裂韧度测试方法 |
CN109781548A (zh) * | 2019-01-25 | 2019-05-21 | 四川大学 | 基于ndb试样的非对称三点弯曲加载测试岩石复合断裂韧度的方法 |
US20210116341A1 (en) * | 2019-12-31 | 2021-04-22 | Chengdu University Of Technology | experimental test method for subcritical propagation rate of rock fractures based on triaxial stress - strain curve |
-
2021
- 2021-05-08 CN CN202110498290.7A patent/CN113239585A/zh not_active Withdrawn
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4299120A (en) * | 1979-03-19 | 1981-11-10 | Terra Tek, Inc. | Method for determining plane strain fracture toughness of non-elastic fracture mechanics specimens |
CN104833590A (zh) * | 2014-07-08 | 2015-08-12 | 吴礼舟 | 岩石ii型断裂韧度测试的新方法 |
CN105043849A (zh) * | 2015-07-08 | 2015-11-11 | 四川大学 | 用于岩石i-ii复合断裂韧度测试的深梁试件及复合断裂韧度测试方法 |
CN109781548A (zh) * | 2019-01-25 | 2019-05-21 | 四川大学 | 基于ndb试样的非对称三点弯曲加载测试岩石复合断裂韧度的方法 |
US20210116341A1 (en) * | 2019-12-31 | 2021-04-22 | Chengdu University Of Technology | experimental test method for subcritical propagation rate of rock fractures based on triaxial stress - strain curve |
Non-Patent Citations (1)
Title |
---|
YING XU.ETC: ""Evaluation of the short core in compression (SCC) method for measuring mode II fracture toughness of rocks"", 《ENGINEERING FRACTURE MECHANICS》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105043849B (zh) | 用于岩石i‑ii复合断裂韧度测试的深梁试件及复合断裂韧度测试方法 | |
CN106932253B (zh) | 测试岩石i-ii复合型动态断裂韧度的试件组件及测试方法 | |
Zhu et al. | Softening damage analysis of gypsum rock with water immersion time based on laboratory experiment | |
Manouchehrian et al. | Numerical analysis of confinement effect on crack propagation mechanism from a flaw in a pre-cracked rock under compression | |
Brook | The equivalent core diameter method of size and shape correction in point load testing | |
Yang | Crack coalescence behavior of brittle sandstone samples containing two coplanar fissures in the process of deformation failure | |
Coviello et al. | On the measurement of the tensile strength of soft rocks | |
CN109781548B (zh) | 基于ndb试样的非对称三点弯曲加载测试岩石复合断裂韧度的方法 | |
Haeri et al. | On the strength and crack propagation process of the pre-cracked rock-like specimens under uniaxial compression | |
CN105424435A (zh) | 一种制作三维节理裂隙圆柱型岩石试样的方法 | |
CN106546482A (zh) | 用于脆性材料i‑iii复合型断裂韧度测试的试件组件及测试方法 | |
CN104833590A (zh) | 岩石ii型断裂韧度测试的新方法 | |
Lin et al. | Determination of strength and deformation properties of columnar jointed rock mass using physical model tests | |
CN109916739A (zh) | 一种测定混凝土i-ii复合型裂缝荷载-位移全曲线的方法 | |
Xiaohu et al. | Quantification of geological strength index based on discontinuity volume density of rock masses | |
Haeri | Crack analysis of pre-cracked brittle specimens under biaxial compression | |
CN106546477A (zh) | 用于脆性材料iii型动态断裂韧度测试的试件组件及测试方法 | |
CN113239585A (zh) | 测量岩石ⅱ型断裂韧度的非对称切槽短棒方法 | |
Haeri et al. | Numerical simulation of compressive to tensile load conversion for determining the tensile strength of ultra-high performance concrete | |
Liu et al. | Loading rate effect on mixed-mode I+ II fracture of V-notched Brazilian disc rock specimen under impact loads | |
CN106769340B (zh) | 岩石试样内部预制裂纹的制作方法 | |
Liu et al. | The coupled influence of surface and internal crack propagation on rock breakages by indentations in biaxial states | |
CN113281148B (zh) | 一种岩石蠕变破坏微裂纹贯通阶段起始时间辨识方法 | |
Li et al. | Characterization of Marcellus Shale fracture properties through size effect tests and Computations | |
Andersen et al. | Evaluating the crack resistance of carbon anodes: Implementation of a measurement system for tensile strength and fracture toughness |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210810 |