CN110658067B - 一种基于全场应变测量判别岩石裂纹应力门槛值的方法 - Google Patents

一种基于全场应变测量判别岩石裂纹应力门槛值的方法 Download PDF

Info

Publication number
CN110658067B
CN110658067B CN201910971125.1A CN201910971125A CN110658067B CN 110658067 B CN110658067 B CN 110658067B CN 201910971125 A CN201910971125 A CN 201910971125A CN 110658067 B CN110658067 B CN 110658067B
Authority
CN
China
Prior art keywords
rock sample
strain
crack
axial
stress
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
CN201910971125.1A
Other languages
English (en)
Other versions
CN110658067A (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.)
Wuhan Institute of Rock and Soil Mechanics of CAS
Original Assignee
Wuhan Institute of Rock and Soil Mechanics of CAS
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 Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Wuhan Institute of Rock and Soil Mechanics of CAS
Priority to CN201910971125.1A priority Critical patent/CN110658067B/zh
Publication of CN110658067A publication Critical patent/CN110658067A/zh
Application granted granted Critical
Publication of CN110658067B publication Critical patent/CN110658067B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • 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
    • G01N1/30—Staining; Impregnating ; Fixation; Dehydration; Multistep processes for preparing samples of tissue, cell or nucleic acid material and the like for analysis
    • 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/02—Details
    • G01N3/06—Special adaptations of indicating or recording means
    • G01N3/068—Special adaptations of indicating or recording means with optical indicating or recording means
    • 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
    • G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • G01N2001/2866—Grinding or homogeneising
    • 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
    • G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • G01N2001/2873—Cutting or cleaving
    • 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/006—Crack, flaws, fracture or rupture
    • G01N2203/0062—Crack or flaws
    • 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/0062—Crack or flaws
    • G01N2203/0066—Propagation of crack
    • 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
    • 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/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0676—Force, weight, load, energy, speed or acceleration
    • 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/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0682—Spatial dimension, e.g. length, area, angle

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于全场应变测量判别岩石裂纹应力门槛值的方法，包含如下步骤：制作岩石试样；岩石试样表面散斑处理；试验加载，采集岩石试样散斑面在试验加载过程中的图像；计算岩石试样散斑面全场位移；计算岩石试样散斑面全场应变；计算岩石试样体积应变、弹性体积应变及裂纹体积应变；获取岩石试样体积应变‑轴向应变曲线、裂纹体积应变‑轴向应变曲线及轴向应力‑轴向应变曲线，并绘制在同一副图中；判断岩石试样体积应变及裂纹体积应变的变化规律，确定岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种裂纹应力门槛值。本发明的优点是：对试验环境要求低，应用更方便，结果更可靠，更具有代表性和普适性。

Description

一种基于全场应变测量判别岩石裂纹应力门槛值的方法
技术领域
本发明涉及岩石力学特性测试领域，具体的说是涉及一种基于全场应变测量判别岩石裂纹应力门槛值的方法。
背景技术
岩石材料在受外力作用下，根据其应力-应变曲线关系可得到四种不同的裂纹应力门槛值，即：裂纹闭合应力σcc、裂纹起裂应力σci、裂纹损伤应力σcd及峰值应力σp。大部分学者认为，岩石裂纹起裂应力为岩石长期强度的下限值，而裂纹损伤应力则为岩石长期强度的上限值。在实际工程中，岩石的长期强度可被用来衡量工程的稳定性，故上述四种不同裂纹应力门槛值的确定及研究对实际工程设计具有较高的参考价值。
目前，测量岩石材料裂纹应力门槛值有如下几种方法，即：应变计算法、声发射法，但以上方法均具有一定的局限性。
其中，应变计算法是先通过在试样表面贴应变片或者安装环向应变规来对试样试验过程中所产生的轴向应变和环向应变进行记录并测量，然后再通过相应计算求得裂纹应力门槛值，此方法存在的缺点是：采取在试样上贴取的应变片数量是有限的，且测量的是试样局部应变变化数据，故所得结果可能不具有代表性；
其中，声发射法则是先通过在试样表面不同位置安装声发射探头来记录试验过程中声发射事件，然后再通过分析声发射事件变化规律来对裂纹应力门槛值进行确定，此方法存在的缺点是：受试验周边噪音影响较大，且岩石材料微观结构的非均质性，往往导致其测试结果在岩石力学测试中不理想。
发明内容
为消除背景技术中传统测量方法所带来的试验环境要求高及不具有代表性的缺点，本发明的目的在于提供一种对试验环境要求低，且测试结果更具有代表性和普适性的基于全场应变测量判别岩石裂纹应力门槛值的方法，其主要是采取的是先基于数字图像相关方法对岩石试样试验时其表面的全场应变进行测量，然后再通过计算得到岩石材料的裂纹应力门槛值。
为实现上述目的，本发明采取的技术方案为：一种基于全场应变测量判别岩石裂纹应力门槛值的方法，包含如下步骤：
S1、制作岩石试样；
S2、对岩石试样表面进行散斑处理；
S3、试验加载，采集岩石试样散斑面在试验加载过程中的图像；
S4、计算岩石试样散斑面的全场位移；
S5、计算岩石试样散斑面的全场应变；
S6、计算岩石试样体积应变、弹性体积应变及裂纹体积应变；
S7、获取岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，并绘制在同一副图中；
S8、判断岩石试样体积应变及裂纹体积应变相对岩石试样轴向应变的变化规律，确定岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值。
上述方法中，所述步骤S1的具体步骤为：先采用相应的切割设备将待测岩石材料切割成试验所需的试样形状及尺寸规格，然后再对其试样表面进行打磨，以便于后续的散斑处理。
上述方法中，所述步骤S2的具体步骤为：先采用相应的喷涂设备将白色哑光漆均匀地涂抹在步骤S1制备完成的岩石试样表面作为背景色，并经过20～30min的凝固后，再采用相应的喷涂设备将黑色哑光漆随机喷洒在已经凝固的白色哑光漆上面，形成黑白相间且可随时岩石试样变形的散斑面。
上述方法中，所述步骤S3的具体步骤为：
试验加载开始前，先将表面经过散斑处理的岩石试样放置在相应的岩石力学试验系统上，并将岩石试样散斑面朝向试验观测方向，然后再在试验观测方向设置一台CCD相机，并调整CCD相机的拍摄方向，使之垂直于岩石试样散斑面，与此同时设置好CCD相机的拍摄速率，最后再设置好岩石力学试验系统的轴向加载速率，以保证岩石试样在加载过程中处于准静态；
试验加载开始时，将岩石力学试验系统与CCD相机拍摄同步触发开始试验，直到岩石试样在轴向压力作用下发生破坏时，立即停止试验；
在整个试验加载过程中，通过CCD相机连续拍摄捕捉岩石试样散斑面的图像。
上述方法中，所述步骤S4的具体步骤为：先将步骤S3得到的岩石试样散斑面图像划分为若干个子区，然后通过下式(1)和(2)，计算出各个子区的x、y向位移，最后再将各子区位移场进行组合，即可得到岩石试样散斑面的全场位移分布；
Figure GDA0002518166400000031
Figure GDA0002518166400000032
式(1)和(2)中，
Figure GDA0002518166400000033
为参考子区位移的一阶偏导数，u0，v0分别为参考子区中心点沿x、y向的位移值；dx为岩石试样散斑面某一子区在x向的位移值，dy为岩石试样散斑面某一子区在y向的位移值；x向为与岩石试样轴向加载方向相垂直的试样侧向，y向为岩石试样的轴向加载方向；x0与y0为岩石试样散斑面某一子区变形前其中心点沿x与y方向的坐标，xi与yi为岩石试样散斑面某一子区变形前其内任一点沿x与y方向的坐标，x'i与y'i为岩石试样散斑面某一子区变形后其内任一点沿x与y方向的坐标；
所述步骤S5的具体步骤为：在通过上式(1)和(2)，得到各子区的x、y向位移后，先通过下式(3)和(4)计算出各子区的x、y向应变，然后将各子区应变进行组合，即可得到岩石试样散斑面的全场应变分布；
Figure GDA0002518166400000041
Figure GDA0002518166400000042
式(3)和(4)中，
Figure GDA0002518166400000043
为参考子区位移的一阶偏导数；εx为岩石试样散斑面某一子区的x向应变，εy为岩石试样散斑面某一子区的y向应变；
所述步骤S6的具体步骤为：在通过上式(3)和(4)得到各子区的x、y向应变后，先通过下式(5)计算出岩石试样的体积应变，然后通过下式(6)计算出岩石试样的弹性体积应变，最后再通过下式(7)求得试样的裂纹体积应变；
Figure GDA0002518166400000044
Figure GDA0002518166400000045
εcv＝εv-εe (7)
式(5)中，m为沿y向的子区总数，n为沿x向的子区总数，εx(i,j)为岩石试样散斑表面第i行，第j列相应子区的x向应变，εy(i,j)为岩石试样散斑表面第i行，第j列相应子区的y向应变，εv为岩石试样的体积应变；
式(6)中，v为岩石试样的泊松比，E为岩石试样的弹性模量，σ为岩石试样轴向应力，εe为岩石试样的弹性体积应变；
式(7)中，εcv为岩石试样的裂纹体积应变，εv为岩石试样的体积应变，εe为岩石试样的弹性体积应变；
所述步骤S7的具体步骤为：先将步骤S6中所得到的岩石试样体积应变εv与岩石试样裂纹体积应变εcv相对于岩石试样的轴向应变分别绘制出相应的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线；然后再将上述的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线与岩石试样的轴向应力-轴向应变曲线一起绘制在同一幅图中；其中，岩石试样的轴向应力为试验加载在试样轴向的轴向力，轴向应变为岩石试样的轴向应力与岩石试样的横截面积比值。
上述方法中，所述步骤S8的具体步骤为：根据步骤S7所绘制在同一副图中的岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，分析判断出体积应变和裂纹体积应变随轴向应变变化的规律，确定出岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值；
当岩石试样裂纹体积应变-轴向应变曲线开始趋于平直稳定段时，此时所对应的轴向应力即为岩石试样的裂纹闭合应力；
当岩石试样裂纹体积应变-轴向应变曲线从平直稳定段开始减小时，此时所对应的轴向应力即为岩石试样的裂纹起裂应力；
当岩石试样体积应变-轴向应变曲线达到其峰值时，其峰值所对应的轴向应力即为岩石试样的裂纹损伤应力；
当岩石试样轴向应力-轴向应变曲线趋于峰值时，其所对应的轴向应力即为岩石试样的峰值应力。
与现有技术相比，本发明的优点是：
与传统的应变测量方法相比，本发明方法所得测试结果更具有代表性及普适性；只需对试样表面进行散斑化处理即可且处理工序十分简单，试验过程中其对周围试验环境要求较低；
与传统的声发射方法相比，本发明方法具有较低的试验环境要求，在应用方面更方便且其结果更可靠。
综上所述，本发明提供的一种基于全场应变测量判别岩石裂纹应力门槛值的方法，在固体试验力学方面具有较好的适用性，且其所得结果相比传统方法更准确，同时对相关科研工作及测试方法的发展具有十分积极的作用。
附图说明
图1为本发明基于全场应变测量判别岩石裂纹应力门槛值方法的流程图；
图2为本发明一种试验实施例的试验示意图；
图3为本发明实施例所得裂纹体积应变，体积应变、轴向应力随轴向应变变化的曲线图；
图2中，“F”表示试验加载过程中，施加在试样上的轴向加载力；
图3中，左边纵坐标轴表示轴向应力，单位为MPa；右边坐标轴表示裂纹体积应变，单位为％；横坐标轴表示轴向应变，单位为％；曲线A表示试样轴向应力-轴向应变曲线，曲线B表示试样体积应变-轴向应变曲线，曲线C表示试样裂纹体积应变-轴向应变曲线。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合附图和具体实施方式，进一步阐述本发明是如何实施的。
参阅图1所示，本发明提供的一种基于全场应变测量判别岩石裂纹应力门槛值的方法，包含如下步骤：
S1、制作岩石试样；
S2、对岩石试样表面进行散斑处理；
S3、试验加载，采集岩石试样散斑面在试验加载过程中的图像；
S4、计算岩石试样散斑面的全场位移；
S5、计算岩石试样散斑面的全场应变；
S6、计算岩石试样体积应变、弹性体积应变及裂纹体积应变；
S7、获取岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，并绘制在同一副图中；
S8、判断岩石试样体积应变及裂纹体积应变相对岩石试样轴向应变的变化规律，确定岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值。
具体的说，步骤S1的具体步骤为：先采用相应的切割设备将待测岩石材料切割成试验所需的试样形状及尺寸规格，然后再对其试样表面进行打磨，以便于后续的散斑处理。
具体的说，步骤S2的具体步骤为：先采用相应的喷涂设备将白色哑光漆均匀地涂抹在步骤S1制备完成的岩石试样表面作为背景色，并经过20～30min的凝固后，再采用相应的喷涂设备将黑色哑光漆随机喷洒在已经凝固的白色哑光漆上面，形成黑白相间且可随时岩石试样变形的散斑面。
具体的说，步骤S3的具体步骤为：
试验加载开始前，先将表面经过散斑处理的岩石试样放置在相应的岩石力学试验系统上，并将岩石试样散斑面朝向试验观测方向，然后再在试验观测方向设置一台CCD相机，并调整CCD相机的拍摄方向，使之垂直于岩石试样散斑面，与此同时设置好CCD相机的拍摄速率，最后再设置好岩石力学试验系统的轴向加载速率，以保证岩石试样在加载过程中处于准静态；
试验加载开始时，将岩石力学试验系统与CCD相机拍摄同步触发开始试验，直到岩石试样在轴向压力作用下发生破坏时，立即停止试验；
在整个试验加载过程中，通过CCD相机连续拍摄捕捉岩石试样散斑面的图像。
具体的说，步骤S4的具体步骤为：先将步骤S3得到的岩石试样散斑面图像划分为若干个子区，然后通过下式(1)和(2)，计算出各个子区的x、y向位移，最后再将各子区位移场进行组合，即可得到岩石试样散斑面的全场位移分布；
Figure GDA0002518166400000081
Figure GDA0002518166400000082
式(1)和(2)中，
Figure GDA0002518166400000083
为参考子区位移的一阶偏导数，u0，v0分别为参考子区中心点沿x、y向的位移值；dx为岩石试样散斑面某一子区在x向的位移值，dy为岩石试样散斑面某一子区在y向的位移值；x向为与岩石试样轴向加载方向相垂直的试样侧向(参阅图2所示,XY坐标中的X轴)，y向为岩石试样的轴向加载方向(参阅图2所示，XY坐标中的Y轴)；x0与y0为岩石试样散斑面某一子区变形前其中心点沿x与y方向的坐标，xi与yi为岩石试样散斑面某一子区变形前其内任一点沿x与y方向的坐标，x'i与y'i为岩石试样散斑面某一子区变形后其内任一点沿x与y方向的坐标；
步骤S5的具体步骤为：在通过上式(1)和(2)，得到各子区的x、y向位移后，先通过下式(3)和(4)计算出各子区的x、y向应变，然后将各子区应变进行组合，即可得到岩石试样散斑面的全场应变分布；
Figure GDA0002518166400000091
Figure GDA0002518166400000092
式(3)和(4)中，
Figure GDA0002518166400000093
为参考子区位移的一阶偏导数；εx为岩石试样散斑面某一子区的x向应变，εy为岩石试样散斑面某一子区的y向应变；
步骤S6的具体步骤为：在通过上式(3)和(4)得到各子区的x、y向应变后，先通过下式(5)计算出岩石试样的体积应变，然后通过下式(6)计算出岩石试样的弹性体积应变，最后再通过下式(7)求得试样的裂纹体积应变；
Figure GDA0002518166400000094
Figure GDA0002518166400000095
εcv＝εv-εe (7)
式(5)中，m为沿y向的子区总数，n为沿x向的子区总数，εx(i,j)为岩石试样散斑表面第i行，第j列相应子区的x向应变，εy(i,j)为岩石试样散斑表面第i行，第j列相应子区的y向应变，εv为岩石试样的体积应变；
式(6)中，v为岩石试样的泊松比，E为岩石试样的弹性模量，σ为岩石试样轴向应力，εe为岩石试样的弹性体积应变；
式(7)中，εcv为岩石试样的裂纹体积应变，εv为岩石试样的体积应变，εe为岩石试样的弹性体积应变；
步骤S7的具体步骤为：先将步骤S6中所得到的岩石试样体积应变εv与岩石试样裂纹体积应变εcv相对于岩石试样的轴向应变分别绘制出相应的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线；然后再将上述的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线与岩石试样的轴向应力-轴向应变曲线一起绘制在同一幅图中；其中，岩石试样的轴向应力为试验加载在试样轴向的轴向力，轴向应变为岩石试样的轴向应力与岩石试样的横截面积比值。
步骤S8的具体步骤为：根据步骤S7所绘制在同一副图中的岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，分析判断出体积应变和裂纹体积应变随轴向应变变化的规律，确定出岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值；
当岩石试样裂纹体积应变-轴向应变曲线开始趋于平直稳定段时，此时所对应的轴向应力即为岩石试样的裂纹闭合应力；
当岩石试样裂纹体积应变-轴向应变曲线从平直稳定段开始减小时，此时所对应的轴向应力即为岩石试样的裂纹起裂应力；
当岩石试样体积应变-轴向应变曲线达到其峰值时，其峰值所对应的轴向应力即为岩石试样的裂纹损伤应力；
当岩石试样轴向应力-轴向应变曲线趋于峰值时，其所对应的轴向应力即为岩石试样的峰值应力。
下面以脆性花岗岩为例详细说明，其是如何通过本发明提供基于全场应变测量判别岩石裂纹应力门槛值的方法，来判别其裂纹应力门槛值的：
第1步，制作本实施例试验所需的花岗岩试样，其具体操作过程为：先采用圆盘切割机将本实施例提供的脆性花岗岩切割成如图2所示的尺寸为100mmX50mmX20mm(高度×宽度×厚度)的长方体试样100，然后采用打磨机将该长方体试样100表面打磨至平整；
第2步，对第1步制作而成的花岗岩试样表面进行散斑处理，得到如图2所示的散斑面101；
在本实施例中，制作该散斑面101的过程具体如下：
先采用喷涂机将白色哑光漆均匀涂抹在上述花岗岩试样表面进行背景色处理，并经过20-30min凝固后，再通过喷涂机将黑色哑光漆随机喷洒已经凝固的白色哑光漆上面，形成黑白相间且可以随上述花岗岩试样变形的散斑面101；
第3步，对第2步完成散斑处理的花岗岩试样进行试验加载处理，并采集本实施例花岗岩试样散斑面101在试验加载过程中的图像；
在本实施例中，加载上述花岗岩试样及采集该采集花岗岩试样散斑面101加载过程中的图像的具体操作过程如下：
试验加载开始前，先将花岗岩试样放置在MTS815.04岩石力学试验系统上，并将花岗岩试样散斑面101朝向试验观测方向，然后在试验观测方向设置一台CCD相机200及若干光源300，并调整CCD相机200的拍摄方向，使之垂直于花岗岩试样散斑面101，与此同时设置好CCD相机200的拍摄速率为50帧/s，最后再设置好岩石力学试验系统的轴向加载速率为0.0005mm/s，以保证花岗岩试样在加载过程中处于准静态；
试验加载开始时，将岩石力学试验系统与CCD相机200拍摄同步触发开始试验，通过岩石力学试验系统向花岗岩试样慢慢施加轴向加载力，直到当花岗岩试样在轴向压力作用下发生破坏时，立即停止试验；在整个试验加载过程，通过CCD相机200连续拍摄捕捉花岗岩试样散斑面101的图像。
第4步，基于第3步得到花岗岩试样散斑面图像，计算出本实施例花岗岩试样散斑面101全场位移；
在本实施例中，计算该散斑面101全场位移的具体操作过程如下：
先通过数字图像相关法将第3步得到的花岗岩试样散斑面图像划分为若干个子区，然后通过下式(1)和(2)，计算出花岗岩试样散斑面101各个子区的x、y向位移，最后再将求得的花岗岩试样散斑面101各子区位移场进行组合，即可得出花岗岩试样散斑面101的全场位移分布；
Figure GDA0002518166400000121
Figure GDA0002518166400000122
上式(1)和(2)中，
Figure GDA0002518166400000123
为参考子区位移的一阶偏导数，u0，v0分别为参考子区中心点沿x、y向的位移值；dx为花岗岩试样散斑面某一子区在x向的位移值，dy为花岗岩试样散斑面某一子区在y向的位移值；x向为与花岗岩试样轴向加载方向相垂直的试样侧向，y向为花岗岩试样的轴向加载方向；x0与y0为花岗岩试样散斑面某一子区变形前其中心点沿x与y方向的坐标，xi与yi为花岗岩试样散斑面某一子区变形前其内任一点沿x与y方向的坐标，x'i与y'i为花岗岩试样散斑面某一子区变形后其内任一点沿x与y方向的坐标；
第5步，计算本实施例花岗岩试样散斑面101全场应变；
在本实施例中，计算该散斑面101全场应变的具体过程如下：
先通过下式(3)和(4)计算出花岗岩试样散斑面101各子区的x、y向应变，然后再将求得的花岗岩试样散斑面101各子区应变进行组合，即可得到花岗岩试样散斑面101的全场应变分布；
Figure GDA0002518166400000131
Figure GDA0002518166400000132
上式(3)和(4)中，
Figure GDA0002518166400000133
为参考子区位移的一阶偏导数；εx为花岗岩试样散斑面某一子区的x向应变，εy为花岗岩试样散斑面某一子区的y向应变；
第6步，计算本实施例花岗岩试样的体积应变εv、弹性体积应变εe及裂纹体积应变εcv；
在本实施例中，计算上述体积应变εv、弹性体积应变εe及裂纹体积应变εcv的具体过程如下：
先通过下式(5)计算出花岗岩试样的体积应变，然后再通过下式(6)计算出花岗岩试样的弹性体积应变，最后再通过下式(7)求得花岗岩试样的裂纹体积应变；
Figure GDA0002518166400000134
Figure GDA0002518166400000135
εcv＝εv-εe (7)
上式(5)中，m为沿y向的子区总数，n为沿x向的子区总数，εx(i,j)为花岗岩试样散斑表面第i行，第j列相应子区的x向应变，εy(i,j)为花岗岩试样散斑表面第i行，第j列相应子区的y向应变，εv为花岗岩试样的体积应变；
上式(6)中，v为花岗岩试样的泊松比，E为花岗岩试样的弹性模量，σ为花岗岩试样轴向应力，εe为花岗岩试样的弹性体积应变；
上式(7)中，εcv为花岗岩试样的裂纹体积应变，εv为花岗岩试样的体积应变，εe为花岗岩试样的弹性体积应变；
第7步，获取本实施例花岗岩试样轴向应力-轴向应变曲线、体积应变-轴向应变曲线及裂纹体积应变-轴向应变曲线，并将他们绘制在同一副图中，参阅图3所示；
第8步，根据第7步绘制的曲线图，分析出本实施例花岗岩试样的体积应变及裂纹体积应变随轴向应变变换的规律，确定出本实施例花岗岩试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值。
图3为采用本发明基于全场应变测量判别岩石裂纹应力门槛值的方法以上述脆性花岗岩为试验实施例，所获得的轴向应力-轴向应变曲线、体积应变-轴向应变曲线及裂纹体积应变-轴向应变曲线；由该图结果分析可知，该脆性花岗岩的裂纹体积应变εcv和体积应变εv随轴向应变变化具有较为明显的特征：其中：其体积应变εv随轴向应变的变化是呈先增加后减小的趋势变化的，且具有一个峰值；而裂纹体积应变εcv随轴向应变的变化是呈先增加后平缓变化再下降的趋势变化，且存在两个较为明显的应变转折点。
如图3所示，在本实施例中，当花岗岩试样裂纹体积应变εcv刚进入平缓变化阶段时(见图3中a点)，其所对应的轴向应力值即为裂纹闭合应力σcc(见图3中b点)，其数值为7.6MPa，此时所对应的轴向应变为0.06％，裂纹闭合应力与峰值应力的比值为6.7％；
当花岗岩试样裂纹体积应变εcv从平缓变化阶段开始减小时(见图3中c点)，其所对应的轴向应力即为裂纹起裂应力σci(见图3中d点)，其数值为35.2MPa，此时所对应的轴向应变为0.116％，裂纹起裂应力与峰值应力的比值为31.2％；
当花岗岩试样体积应变εv达到峰值强度时(见图3中e点)，其所对应的轴向应力即为裂纹损伤应力σcd，其数值为82.3MPa，此时所对应的轴向应变为0.184％，裂纹损伤应力与峰值应力的比值为72.8％；
当花岗岩试样轴向应力-轴向应变曲线趋于峰值时(见图3中f点)，其所对应的轴向应力即为岩石试样的峰值应力σp。
另外说明，对比传统应变计算法或声发射方法获得的裂纹应力门槛值，由传统应变计算法或声发射方法获得的裂纹起裂应力与峰值应力的比值范围普遍为0.3-0.5，裂纹损伤应力与峰值应力的比值范围普遍为0.7-0.85。由上述实施例可知，本发明所确定的以上两个应力(即裂纹起裂应力和裂纹损伤应力)与峰值应力的比值均处于上述比值范围内，故而本发明提供的基于全场应变测量判别岩石裂纹应力门槛值方法所得到的结果能够符合大量数据的统计结果，且具有普适性和更精确的特点。
最后说明，以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (6)

1.一种基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：包含如下步骤：
S1、制作岩石试样；
S2、对岩石试样表面进行散斑处理；
S3、试验加载，采集岩石试样散斑面在试验加载过程中的图像；
S4、计算岩石试样散斑面的全场位移：即先将步骤S3得到的岩石试样散斑面图像划分为若干个子区，然后通过下式(1)和(2)，计算出各个子区的x、y向位移，最后再将各子区位移场进行组合，即可得到岩石试样散斑面的全场位移分布；
Figure FDA0002518166390000011
Figure FDA0002518166390000012
式(1)和(2)中，
Figure FDA0002518166390000013
为参考子区位移的一阶偏导数，u0，v0分别为参考子区中心点沿x、y向的位移值；dx为岩石试样散斑面某一子区在x向的位移值，dy为岩石试样散斑面某一子区在y向的位移值；x向为与岩石试样轴向加载方向相垂直的试样侧向，y向为岩石试样的轴向加载方向；x0与y0为岩石试样散斑面某一子区变形前其中心点沿x与y方向的坐标，xi与yi为岩石试样散斑面某一子区变形前其内任一点沿x与y方向的坐标，xi'与yi'为岩石试样散斑面某一子区变形后其内任一点沿x与y方向的坐标；
S5、计算岩石试样散斑面的全场应变：即在通过上式(1)和(2)，得到各子区的x、y向位移后，先通过下式(3)和(4)计算出各子区的x、y向应变，然后将各子区应变进行组合，即可得到岩石试样散斑面的全场应变分布；
Figure FDA0002518166390000021
Figure FDA0002518166390000022
式(3)和(4)中，
Figure FDA0002518166390000023
为参考子区位移的一阶偏导数；εx为岩石试样散斑面某一子区的x向应变，εy为岩石试样散斑面某一子区的y向应变；
S6、计算岩石试样体积应变、弹性体积应变及裂纹体积应变：即在通过上式(3)和(4)得到各子区的x、y向应变后，先通过下式(5)计算出岩石试样的体积应变，然后通过下式(6)计算出岩石试样的弹性体积应变，最后再通过下式(7)求得试样的裂纹体积应变；
Figure FDA0002518166390000024
Figure FDA0002518166390000025
εcv＝εv-εe (7)
式(5)中，m为沿y向的子区总数，n为沿x向的子区总数，εx(i,j)为岩石试样散斑表面第i行，第j列相应子区的x向应变，εy(i,j)为岩石试样散斑表面第i行，第j列相应子区的y向应变，εv为岩石试样的体积应变；
式(6)中，v为岩石试样的泊松比，E为岩石试样的弹性模量，σ为岩石试样轴向应力，εe为试样的弹性体积应变；
式(7)中，εcv为岩石试样的裂纹体积应变，εv为岩石试样的体积应变，εe为岩石试样的弹性体积应变；
S7、获取岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，并绘制在同一幅 图中；
S8、判断岩石试样体积应变及裂纹体积应变相对岩石试样轴向应变的变化规律，确定岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值。
2.根据权利要求1所述的基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：所述步骤S1的具体步骤为：先采用相应的切割设备将待测岩石材料切割成试验所需的试样形状及尺寸规格，然后再对其试样表面进行打磨，以便于后续的散斑处理。
3.根据权利要求1所述的基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：所述步骤S2的具体步骤为：先采用相应的喷涂设备将白色哑光漆均匀地涂抹在步骤S1制备完成的岩石试样表面作为背景色，并经过20～30min的凝固后，再采用相应的喷涂设备将黑色哑光漆随机喷洒在已经凝固的白色哑光漆上面，形成黑白相间且可随时岩石试样变形的散斑面。
4.根据权利要求1所述的基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：所述步骤S3的具体步骤为：
试验加载开始前，先将表面经过散斑处理的岩石试样放置在相应的岩石力学试验系统上，并将岩石试样散斑面朝向试验观测方向，然后再在试验观测方向设置一台CCD相机，并调整CCD相机的拍摄方向，使之垂直于岩石试样散斑面，与此同时设置好CCD相机的拍摄速率，最后再设置好岩石力学试验系统的轴向加载速率，以保证岩石试样在加载过程中处于准静态；
试验加载开始时，将岩石力学试验系统与CCD相机拍摄同步触发开始试验，直到岩石试样在轴向压力作用下发生破坏时，立即停止试验；
在整个试验加载过程中，通过CCD相机连续拍摄捕捉岩石试样散斑面的图像。
5.根据权利要求1所述的基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：所述步骤S7的具体步骤为：先将步骤S6中所得到的岩石试样体积应变εv与岩石试样裂纹体积应变εcv相对于岩石试样的轴向应变分别绘制出相应的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线；然后再将上述的体积应变-轴向应变曲线和裂纹体积应变-轴向应变曲线与岩石试样的轴向应力-轴向应变曲线一起绘制在同一幅图中；其中，岩石试样的轴向应力为试验加载在试样轴向的轴向力，轴向应变为岩石试样的轴向应力与岩石试样的横截面积比值。
6.根据权利要求1所述的基于全场应变测量判别岩石裂纹应力门槛值的方法，其特征在于：所述步骤S8的具体步骤为：根据步骤S7所绘制在同一幅 图中的岩石试样体积应变-轴向应变曲线、裂纹体积应变-轴向应变曲线及轴向应力-轴向应变曲线，分析判断出体积应变和裂纹体积应变随轴向应变变化的规律，确定出岩石试样包含裂纹闭合应力、裂纹起裂应力、裂纹损伤应力及峰值应力在内的这四种不同裂纹应力门槛值；
当岩石试样裂纹体积应变-轴向应变曲线开始趋于平直稳定段时，此时所对应的轴向应力即为岩石试样的裂纹闭合应力；
当岩石试样裂纹体积应变-轴向应变曲线从平直稳定段开始减小时，此时所对应的轴向应力即为岩石试样的裂纹起裂应力；
当岩石试样体积应变-轴向应变曲线达到其峰值时，其峰值所对应的轴向应力即为岩石试样的裂纹损伤应力；
当岩石试样轴向应力-轴向应变曲线趋于峰值时，其所对应的轴向应力即为岩石试样的峰值应力。
CN201910971125.1A 2019-10-14 2019-10-14 一种基于全场应变测量判别岩石裂纹应力门槛值的方法 Active CN110658067B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910971125.1A CN110658067B (zh) 2019-10-14 2019-10-14 一种基于全场应变测量判别岩石裂纹应力门槛值的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910971125.1A CN110658067B (zh) 2019-10-14 2019-10-14 一种基于全场应变测量判别岩石裂纹应力门槛值的方法

Publications (2)

Publication Number Publication Date
CN110658067A CN110658067A (zh) 2020-01-07
CN110658067B true CN110658067B (zh) 2020-08-04

Family

ID=69040770

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910971125.1A Active CN110658067B (zh) 2019-10-14 2019-10-14 一种基于全场应变测量判别岩石裂纹应力门槛值的方法

Country Status (1)

Country Link
CN (1) CN110658067B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111413199A (zh) * 2020-04-10 2020-07-14 中国建筑第六工程局有限公司 一种确定岩石应变局部化启动应力水平的方法
CN112945724B (zh) * 2021-01-28 2023-09-15 中国科学院武汉岩土力学研究所 基于纳米压痕的裂纹尖端断裂过程区域确定方法及设备
CN112801982B (zh) * 2021-01-28 2022-06-17 中国科学院武汉岩土力学研究所 确定岩石混合压剪裂纹尖端断裂过程区长度的方法及设备
CN112903443B (zh) * 2021-01-28 2021-12-14 中国科学院武汉岩土力学研究所 确定岩石材料裂纹尖端断裂过程区传播速度的方法及设备
CN113155017B (zh) * 2021-04-23 2022-05-24 西南石油大学 一种山区管道应变在线监测系统
CN113567243B (zh) * 2021-07-15 2022-05-03 中南大学 一种基于轴向应力响应的裂纹闭合应力确定方法
CN113776942B (zh) * 2021-09-17 2022-05-20 东北大学 一种岩石三轴压缩下识别闭合应力和启裂应力的试验方法
CN114858551B (zh) * 2022-04-28 2023-07-14 煤炭科学研究总院有限公司 一种可实现监测层状试样分层界面处的变形场的方法
CN114858076B (zh) * 2022-04-28 2023-07-18 煤炭科学研究总院有限公司 一种可实现监测样品内部变形场的方法
CN116296820B (zh) * 2023-03-31 2024-04-05 四川大学 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品
CN118190631A (zh) * 2024-05-15 2024-06-14 四川大学 一种岩石应力阈值确定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0457788B1 (de) * 1989-02-13 1994-05-04 Fraunhofer-Gesellschaft Zur Förderung Der Angewandten Forschung E.V. Verfahren und vorrichtung zum überwachen des auftretens oder ausbreitens von brüchen oder rissen in materialien, werkstücken, bauteilen, bauwerken oder himmelskörpern, insbesondere in erdbebengefährdeten gebieten oder bei gebirgsbewegungen im tunnelbau
CN104502198A (zh) * 2015-01-05 2015-04-08 大连理工大学 一种基于数字图像相关的裂纹j积分测量方法
WO2015056790A1 (ja) * 2013-10-18 2015-04-23 国立大学法人佐賀大学 亀裂検出システム及び亀裂検出方法
CN106124311A (zh) * 2016-07-29 2016-11-16 中国科学院武汉岩土力学研究所 基于应变测试判别裂纹扩展演化过程的方法
CN106294984A (zh) * 2016-08-08 2017-01-04 中国科学院武汉岩土力学研究所 一种基于细观裂纹增长率确定岩体损伤门槛的方法
CN107907409A (zh) * 2017-11-10 2018-04-13 中国地质大学（武汉） 一种确定岩石起裂应力的方法、设备及存储设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0457788B1 (de) * 1989-02-13 1994-05-04 Fraunhofer-Gesellschaft Zur Förderung Der Angewandten Forschung E.V. Verfahren und vorrichtung zum überwachen des auftretens oder ausbreitens von brüchen oder rissen in materialien, werkstücken, bauteilen, bauwerken oder himmelskörpern, insbesondere in erdbebengefährdeten gebieten oder bei gebirgsbewegungen im tunnelbau
WO2015056790A1 (ja) * 2013-10-18 2015-04-23 国立大学法人佐賀大学 亀裂検出システム及び亀裂検出方法
CN104502198A (zh) * 2015-01-05 2015-04-08 大连理工大学 一种基于数字图像相关的裂纹j积分测量方法
CN106124311A (zh) * 2016-07-29 2016-11-16 中国科学院武汉岩土力学研究所 基于应变测试判别裂纹扩展演化过程的方法
CN106294984A (zh) * 2016-08-08 2017-01-04 中国科学院武汉岩土力学研究所 一种基于细观裂纹增长率确定岩体损伤门槛的方法
CN107907409A (zh) * 2017-11-10 2018-04-13 中国地质大学（武汉） 一种确定岩石起裂应力的方法、设备及存储设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于DIC类岩石材料裂纹扩展试验与数值模拟研究;徐丽海;《中国优秀硕士学位论文全文数据库 基础科学辑》;20190115(第1期);第49-74页 *

Also Published As

Publication number Publication date
CN110658067A (zh) 2020-01-07

Similar Documents

Publication Publication Date Title
CN110658067B (zh) 一种基于全场应变测量判别岩石裂纹应力门槛值的方法
CN110686971B (zh) 一种精确判断裂隙岩体起裂应力及裂纹类型的方法
CN108548834B (zh) 一种基于sem原位成像系统的应变测量方法
CN110779797B (zh) 一种测量金属拉伸试验过程中塑性应变比的方法
KR101720845B1 (ko) 인장시험에서 디지털 이미지 해석 기법을 이용하여 고 변형 구간에서 재료의 진응력-진변형률 곡선을 추출하는 방법
JP5435352B2 (ja) 板状材料の破断ひずみ特定方法
CN111426563A (zh) 一种3d打印试样拉伸破坏性能分析方法
CN104777046B (zh) 基于小时间尺度的疲劳裂纹扩展机理测试方法
CN112629725B (zh) 一种活塞杆喷涂层残余应力测试方法
CN106018179A (zh) 一种基于图像处理的胶料表面黏性测量方法及装置
CN104122205B (zh) 一种利用压痕隆起量测量残余应力的方法
CN107991149B (zh) 一种利用单向拉伸获得金属板材拉剪复合变形的方法
CN109115602A (zh) 一种基于轴向压力确定柱体承压截面偏移量的方法和系统
CN106525584A (zh) 一种用于透射电子显微镜下原位力学加载工具的制备方法
CN108956665B (zh) 一种用于脆性材料微观结构研究的中子测量方法
KR20200045335A (ko) 소성 응력 확대 계수 산출 시스템
CN113670722A (zh) 一种金属抗拉伸强度性能的检测方法
CN111982946A (zh) 一种测定高碳钢盘条珠光体球团及领域尺寸的方法
CN210037436U (zh) 用于涂层材料原位拉伸观测系统的原位拉伸设备
CN112331574A (zh) 半导体器件高度分布的检测方法
CN112504149A (zh) 一种基于dic技术的板料扩孔率获取方法
CN110455711A (zh) 一种聚合物涂层材料界面结合强度测试装置
CN113218796B (zh) 材料应变-硬度关系检测方法
CN117929121A (zh) 一种围岩相似材料裂纹与位移演化测量方法
CN110029302B (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