CN109060539B - 一种岩石微米尺度弹性模量及屈服强度获取方法 - Google Patents
一种岩石微米尺度弹性模量及屈服强度获取方法 Download PDFInfo
- Publication number
- CN109060539B CN109060539B CN201811099088.1A CN201811099088A CN109060539B CN 109060539 B CN109060539 B CN 109060539B CN 201811099088 A CN201811099088 A CN 201811099088A CN 109060539 B CN109060539 B CN 109060539B
- Authority
- CN
- China
- Prior art keywords
- rock
- elastic modulus
- indentation
- yield strength
- 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
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
- 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/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- 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
-
- 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/02—Details not specific for a particular testing method
- G01N2203/026—Specifications of the specimen
- G01N2203/0262—Shape of the specimen
- G01N2203/0278—Thin specimens
- G01N2203/0282—Two dimensional, e.g. tapes, webs, sheets, strips, disks or membranes
Abstract
本发明公开了一种微尺度岩石弹性模量及屈服强度参数的获取方法。所述方法包括:采用微米级压头开展岩石微米压痕实验,获取加载过程中的载荷—位移曲线,结合压痕实验公式可获取不同位移条件下岩石弹性模量;对岩石试样开展微CT扫描,建立压痕区域岩石骨架的有限元网格模型;以微米压痕实验得到的弹性模量作为输入参数模拟岩石单轴压缩过程,获得模型整体弹性模量并与岩芯单轴压缩实验对比,确定有效表征岩石微米弹性模量的压入深度RVE;然后,开展不同屈服强度条件下的岩样压痕实验数值模拟,并将模拟得到的加卸载载荷—位移曲线与压痕实验进行对比验证,从而确定岩石微米级的屈服强度。
Description
技术领域
本发明涉及岩土工程领域,特别涉及一种岩石微尺度力学参数的获取方法。
背景技术
岩石作为一种非均质多孔介质材料,其宏观力学性质与孔隙结构特征、矿物组分密切相关。大量理论、实验表明,岩石的变形及破裂源于微观。微尺度的岩石变形及破裂研究需要相应尺度的力学参数。然而,传统单轴、三轴压缩等试验的岩样尺寸均在厘米级以上,得到的岩石力学参数只能满足宏观岩石力学性质的研究,这些实验在岩样的制取、精度测量根本无法应用到微米级别岩石的参数的研究。对于微米尺度的材料力学性能测试也面临着众多难题,例如样品尺寸测量、应变—位移的测定以及样品的制备等。
针对目前岩石微米尺度力学参数难以获取的缺陷,本发明公开一种岩石微米尺度弹性模量及屈服强度参数的实验及数值模拟方法。
发明内容
本发明的目的在于提供一种可用于岩石微尺度弹性模量及屈服强度的获取方法,通过综合运用岩石微米压痕实验及数字岩心模拟技术,解决目前实验手段难以获取岩石微尺度力学参数的缺陷。
为了达到上述目的,本发明是这样实现的:
(一)将原始岩心制成厚度为5mm的切片,并进行抛光处理,确保切片上下两个底面平行、光滑。利用制备的岩样切片开展压痕实验,获取仪器压头压入样品过程中的载荷—位移曲线,利用压痕实验公式得出不同压入深度条件下的岩石弹性模量;对切片后岩样开展单轴压缩试验,获得宏观尺度下岩石弹性模量数据。
(二)对压痕实验中切片的压痕部分制样并开展微观CT扫描,获得岩石微观结构图像;基于数字岩心重建技术,建立岩石骨架的有限元网格模型;
(三)以岩石微米压痕实验得到的岩石弹性模量作为输入参数,开展单轴压缩数值模拟,并将模型整体表现出的弹性模量与单轴压缩实验结果进行对比验证,以确定可表征该岩样微米级弹性模量的压入深度RVE。
(四)基于数字岩心重建技术,建立岩石骨架及压头的有限元网格模型;假设岩石为理想Von-Mises各向同性强化弹塑性材料,压头和测试样品为刚—柔接触方式,遵循基本库仑—摩擦接触模型,计算过程考虑几何大变形;开展压痕实验过程的数值模拟,通过模拟不同屈服强度条件下的加卸载曲线并与压痕实验进行比对,从而确定该岩样微米尺度的屈服强度参数。
与传统技术相比,本发明的有益效果在于:
解决了传统技术在实验设备、方法、样品制取和测量精度等方面均无法应用至微观尺度的缺陷,为岩石微米级力学参数获取提供了新的研究方法。
附图说明
为了更清楚地说明本发明方法的技术方案,下面结合附图和具体实施方式对本申请实施例作进一步说明。
图1为本发明实施例所述的岩石微米级弹性模量及屈服强度获取方法流程图。
图2为本发明实施例提供的一个岩样切片及实验压痕分布图。
图3为本发明实施例提供的一个岩样压痕实验的加卸载过程位移—载荷分布曲线。
图4为本发明实施例提供的一个岩样切片压痕区的微CT图像及其骨架数字岩心模型。
图5为本发明实施例提供的岩样单轴压缩数值模拟得到的岩样应力及应变分布云图。
图6为本发明实施例提供的岩样压痕实验得到的弹性模量、数字岩心模型弹性模量与岩样单轴压缩弹性模量对比曲线。
图7为本发明实施例提供的压痕实验数值模拟采用的一个岩样骨架及压头模型图。
图8为本发明实施例提供的压痕实验数值模拟得到的岩样S7不同屈服强度下的位移—载荷曲线与压痕实验曲线对比图。
具体实施方式
为了使本发明实现的技术手段、达成目的和模型功效易于说明,下面结合附图及实施例,对本申请实施例中的技术方案进行详细描述。需要知悉的是,所描述的实施例仅是本申请一部分实施例,而不是全部实施例。基于本申请的实施例,本领域的其他技术人员在没有其他创新性劳动的前提下获得的所有其他实施例,都应当属于本申请保护范围。
图1是本申请一种岩石微米尺度弹性模量及屈服强度获取方法实施例的流程图,包括以下步骤。
S1:制备用于微米压痕实验的岩心切片,从原始岩样上制取2mm厚的岩石切片,并对其进行表面抛光处理,同时保证岩样切片上下表面的平行。本实施例共选用两类岩样作为研究对象,其中S7为砂岩,MS1为人造砂岩。
压痕实验过程将试样固定在样品台处,保证岩石切片下表面与样品台的完全接触,调整样品台位置以选取合适的测点位置;通过调整仪器右侧旋钮使实验仪压头至刚接触岩石样品的位置,点击实验仪配套软件界面上的下降按钮,对实验仪压头位置进行微调,直至压头完全接触;以50mN/s的加载速度加载至目标载荷,停留5s后开始卸载,记录加卸载力与位移。图2所示为本实施例所制备的岩样切片及实验压痕分布。
此时测试材料的简约弹性模量由下式计算。
式中,Ac是压头与测试样品接触区域沿加载方向的投影面积,是接触深度hc的函数。其中,S是卸载曲线的初始斜率,Pmax是压入载荷的最大值,可由如图3所示的本实施例一次实验得到的载荷—位移曲线求得。
对于圆锥形压头,
Ac=24.5(hc+0.011427)2 (3)
其中,α为中心线与锥面的夹角。
测试材料的弹性模量与简约弹性模量的关系为:
其中,E和ν为测试材料的弹性模量和泊松比(砂岩泊松比取0.31),E1和ν1为压头的弹性模量和泊松比。压头为金刚石,其弹性模量为1141GPa,泊松比为0.07。
S2:对压痕区域岩样开展微CT扫描,获得样品微结构特征图像;结合数字岩心技术,构建压痕处岩石骨架的有限元网格模型;图4所示为本实施例的一个岩样微CT图像及其岩石骨架有限元网格模型。
S3:将重建的岩石骨架数字岩心模型导入ANSYS软件,以不同压入深度条件下压痕实验得到的弹性模量作为输入参数,开展岩石单轴压缩数值模拟,模型沿z轴方向的上表面施加均布压力荷载,下表面施加固定约束条件;通过监测不同均布载荷条件下上表面的应变量,确定模型的弹性模量。图5所示为模型MS1-③在10MPa均布载荷下的Von-Mises应力及应变分布场图。当数值模拟得到的模型整体弹性模量与单轴压缩试验一致时,该压入深度为该岩样微米力学参数的RVE,此时压痕实验得到的为其微米尺度弹性模量。图6所示为本实施例的岩样压痕实验得到的弹性模量、数字岩心模型弹性模量与岩样单轴压缩弹性模量对比曲线。
S4:采用数字岩心建模技术,建立岩石骨架及压头的有限元网格模型。图7所示为本实施例采用的一个岩样骨架及压头的几何模型示意图,岩石模型尺寸为2mm×2mm×2mm,压头为120°圆锥体。鉴于实验过程压入深度在100μm以内,压头模型高100μm。由于该模型计算过程中涉及材料非线性与几何非线性的计算问题,网格质量及载荷步长显著影响数值模拟结果的收敛性。因此将圆锥形金刚石压头假定为刚体,有效减少了接触区域的网格数目,避免了因接触区域计算节点过多造成的收敛困难。假设岩石为理想Von-Mises各向同性强化弹塑性材料,压头和测试样品为刚—柔接触方式,遵循基本库仑—摩擦接触模型,计算过程考虑几何大变形。模型弹性模量数据选用压痕实验对应测点的弹性模量数据,采用多步加载载荷的方式,每个加载步位移收敛条件设置为0.1μm,接触面法向力收敛条件设置为1mN,应变能收敛条件设置为0.1mJ。
假设该模型微尺度屈服强度在某一区间上,通过模拟不同屈服强度条件下的加卸载曲线并与压痕实验进行对比,即可不断缩小真实屈服强度的所在区间,并最终确定其数值。以模型S7-①为例,数值模拟过程中所采用的屈服强度依次为50—70—80—83—85—100—120MPa,基于数值模拟结果绘制出相应的加卸载曲线,如图8所示,当屈服强度Y=83MPa时,数值模拟加卸载曲线与实验曲线较好吻合,因此微米尺度下岩样S7-①的屈服强度为83MPa,该数值大约为单轴压缩测得屈服强度(37MPa)的2.24倍。
以上所述仅为本发明的较佳实施例而已,用于描述本发明的基本原理、特征和主要优点,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种微尺度岩石弹性模量及屈服强度参数的获取方法,其特征在于,该方法包括:
利用制备的岩样切片开展压痕实验,获取压头压入样品过程中的载荷—位移曲线,利用压痕实验公式得出不同压入深度条件下的岩石微米级弹性模量;对切片后剩余岩样开展单轴压缩试验,获得宏观尺度下岩石弹性模量数据;
开展微观CT扫描,获得岩石微观结构图像;基于数字岩心重建技术,建立岩石骨架的有限元网格模型;
开展单轴压缩数值模拟,并将模型整体求得的弹性模量与单轴压缩实验结果进行对比验证,以确定可表征该岩样微米级弹性模量的压入深度;
基于数字岩心重建技术,建立岩石骨架及压头的有限元网格模型;开展压痕实验过程的数值模拟,通过模拟不同屈服强度条件下的加卸载过程载荷—位移曲线并与实验结果进行比对,从而确定该岩样微米尺度的屈服强度。
2.根据权利要求1所述的一种微尺度岩石弹性模量及屈服强度参数的获取方法,其特征在于,所述岩样切片厚度为5mm,样品表面进行抛光处理,确保切片上下两个底面平行、光滑。
3.根据权利要求1所述的一种微尺度岩石弹性模量及屈服强度参数的获取方法,其特征在于,对压痕实验中的岩样压痕部位钻取小岩心样品并利用微观CT对小岩样开展扫描成像。
4.根据权利要求1所述的一种微尺度岩石弹性模量及屈服强度参数的获取方法,其特征在于,采用ANSYS软件模拟岩石单轴压缩过程,以岩样微米压痕实验得到的弹性模量作为输入参数,模型沿z轴方向的上表面施加均布压力荷载,下表面施加固定约束条件;通过监测不同均布载荷条件下上表面的应变量,确定模型的弹性模量。
5.根据权利要求1所述的一种微尺度岩石弹性模量及屈服强度参数的获取方法,其特征在于,压痕实验数值模拟过程中,假设岩石为理想Von-Mises各向同性强化弹塑性材料,压头和测试样品为刚—柔接触方式,遵循基本库仑—摩擦接触模型,计算过程考虑几何大变形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811099088.1A CN109060539B (zh) | 2018-09-20 | 2018-09-20 | 一种岩石微米尺度弹性模量及屈服强度获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811099088.1A CN109060539B (zh) | 2018-09-20 | 2018-09-20 | 一种岩石微米尺度弹性模量及屈服强度获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109060539A CN109060539A (zh) | 2018-12-21 |
CN109060539B true CN109060539B (zh) | 2021-08-24 |
Family
ID=64762207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811099088.1A Active CN109060539B (zh) | 2018-09-20 | 2018-09-20 | 一种岩石微米尺度弹性模量及屈服强度获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109060539B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2789252C1 (ru) * | 2022-06-06 | 2023-01-31 | Общество с ограниченной ответственностью "Тюменский нефтяной научный центр" (ООО "ТННЦ") | Способ определения горизонтальных напряжений в массиве горных пород |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109612859B (zh) * | 2019-01-18 | 2024-02-27 | 四川大学 | 微观单元磨蚀性表达方法测量装置及检测方法 |
CN110765572B (zh) * | 2019-09-12 | 2021-04-09 | 中国科学院武汉岩土力学研究所 | 一种杏仁状玄武岩单三轴试验连续非连续数值模拟方法 |
CN110579400B (zh) * | 2019-09-25 | 2021-11-26 | 西南石油大学 | 一种脆性岩石微尺度强度及残余强度的测量计算方法 |
CN111398273B (zh) * | 2019-10-09 | 2021-07-06 | 天津大学 | 一种岩石几何-力学参数获取方法及全息扫描系统 |
CN111157338B (zh) * | 2020-01-08 | 2021-09-17 | 吉林大学 | 一种金属哑铃型试样大应变范围压缩硬化曲线的测量方法 |
CN111189701B (zh) * | 2020-01-08 | 2021-09-17 | 吉林大学 | 一种金属双曲线试样的大应变压缩硬化曲线的测量方法 |
CN111257536A (zh) * | 2020-01-20 | 2020-06-09 | 中国科学院武汉岩土力学研究所 | 一种岩石力学与储层工程参数评估方法 |
CN112179769A (zh) * | 2020-09-29 | 2021-01-05 | 西南石油大学 | 一种基于岩屑微纳米压痕实验的页岩弹性模量评价方法 |
CN112268825B (zh) * | 2020-10-16 | 2023-04-07 | 常州大学 | 一种基于纳米压痕测试的岩石强度获取方法 |
CN112525684A (zh) * | 2020-11-20 | 2021-03-19 | 西南石油大学 | 一种基于数字岩心技术的岩石粘聚力—内摩擦角计算方法 |
CN113138106B (zh) * | 2021-04-15 | 2022-08-30 | 东北石油大学 | 基于随钻岩屑录井资料的岩石弹性参数确定方法 |
CN113029746B (zh) * | 2021-04-16 | 2023-11-03 | 中国矿业大学(北京) | 一种通过细/微观力学测试来确定页岩宏观模量的方法 |
CN113125260B (zh) * | 2021-04-16 | 2023-08-18 | 西北工业大学 | 一种小尺寸非标准试验件弹性模量测试方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288499A (zh) * | 2011-08-31 | 2011-12-21 | 湖南大学 | 一种识别焊缝不同区域材料静态力学性能参数的检测方法 |
CN102313677A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊缝不同区域材料动态力学性能参数的检测方法 |
CN102313676A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊点不同区域材料静态力学性能参数的检测方法 |
CN102435514A (zh) * | 2011-08-31 | 2012-05-02 | 湖南大学 | 一种识别焊点不同区域材料动态力学性能参数的检测方法 |
CN104330298A (zh) * | 2014-10-30 | 2015-02-04 | 北京理工大学 | 表面梯度金属材料的细微观力学性能评估测量方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2989465B1 (fr) * | 2012-04-12 | 2014-11-21 | Total Sa | Procede de determination de parametres geomecaniques d'un echantillon de roche |
-
2018
- 2018-09-20 CN CN201811099088.1A patent/CN109060539B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288499A (zh) * | 2011-08-31 | 2011-12-21 | 湖南大学 | 一种识别焊缝不同区域材料静态力学性能参数的检测方法 |
CN102313677A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊缝不同区域材料动态力学性能参数的检测方法 |
CN102313676A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊点不同区域材料静态力学性能参数的检测方法 |
CN102435514A (zh) * | 2011-08-31 | 2012-05-02 | 湖南大学 | 一种识别焊点不同区域材料动态力学性能参数的检测方法 |
CN104330298A (zh) * | 2014-10-30 | 2015-02-04 | 北京理工大学 | 表面梯度金属材料的细微观力学性能评估测量方法及装置 |
Non-Patent Citations (3)
Title |
---|
"Nanoindentation Test Integrated with Numerical Simulation to Characterize Mechanical Properties of Rock Materials";Hoki Ban et al.;《Journal of Testing and Evaluation》;20140531;第42卷(第3期);第787-796页 * |
"基于微米压痕实验研究页岩力学特性";陈平 等;《石油勘探与开发》;20151031;第662-670页 * |
"多尺度复合材料力学研究进展";陈玉丽 等;《固体力学学报》;20180228;第39卷(第1期);第1-68页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2789252C1 (ru) * | 2022-06-06 | 2023-01-31 | Общество с ограниченной ответственностью "Тюменский нефтяной научный центр" (ООО "ТННЦ") | Способ определения горизонтальных напряжений в массиве горных пород |
Also Published As
Publication number | Publication date |
---|---|
CN109060539A (zh) | 2018-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109060539B (zh) | 一种岩石微米尺度弹性模量及屈服强度获取方法 | |
Zhu et al. | Replication of internal defects and investigation of mechanical and fracture behaviour of rock using 3D printing and 3D numerical methods in combination with X-ray computerized tomography | |
Jiang et al. | A simple and efficient approach to capturing bonding effect in naturally microstructured sands by discrete element method | |
Lu et al. | Meso-scale finite element model for FRP sheets/plates bonded to concrete | |
Zhao et al. | Study on failure characteristic of rock‐like materials with an open‐hole under uniaxial compression | |
CN103792143A (zh) | 一种单轴拉伸全程真应力应变曲线的快速获取方法 | |
CN106124316B (zh) | 散粒体岩土材料各向异性的室内真三轴试验方法 | |
Kim et al. | Microstructure characterization of cement paste from micro-CT and correlations with mechanical properties evaluated from virtual and real experiments | |
Xiao et al. | Experimental study on the mechanical properties of simulated columnar jointed rock masses | |
JPWO2010084840A1 (ja) | 押込試験方法および押込試験装置 | |
Zhou et al. | An experimental investigation of tensile fracturing behavior of natural and artificial rocks in static and dynamic Brazilian disc tests | |
CN111666699B (zh) | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 | |
Mohammed et al. | Modelling the microstructural evolution and fracture of a brittle confectionery wafer in compression | |
Han et al. | A relationship between tensile strength and loading stress governing the onset of mode I crack propagation obtained via numerical investigations using a bonded particle model | |
Lei et al. | Deformation localization and cracking processes of sandstone containing two flaws of different geometric arrangements | |
Tian et al. | Anisotropic creep behavior of soft-hard interbedded rock masses based on 3D printing and digital imaging correlation technology | |
CN104122205B (zh) | 一种利用压痕隆起量测量残余应力的方法 | |
Volz et al. | Comparison of three approaches to determine the projected area in contact from finite element Berkovich nanoindentation simulations in tungsten | |
Miguel et al. | Influence of the width of the loading strip in the Brazilian tensile test of concrete and other brittle materials | |
Li et al. | Three-dimensional micro flow-stress-damage (FSD) model and application in hydraulic fracturing in brittle and heterogeneous rocks | |
LeBaron et al. | Micromechanical analysis of energy release in snow fracture | |
CN111539140A (zh) | 一种有限元模拟单轴贯入强度试验的建模方法,设备及可读存储介质 | |
Huang et al. | 3D in-situ XCT image based meso-scale fracture modelling and validation of concrete using voxel hexahedron meshing and damage plasticity model | |
Wang et al. | Experimental study on the mechanical behavior and failure characteristics of rock analogs with filled internal fractures: A new method by sand powder 3D printing | |
Kuksenko et al. | Prediction of damage initiation and simulation of damage propagation in 3D woven composites during processing |
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 |