CN114689448B - 张开裂隙压密闭合及时效扩展的损伤断裂分析方法 - Google Patents
张开裂隙压密闭合及时效扩展的损伤断裂分析方法 Download PDFInfo
- Publication number
- CN114689448B CN114689448B CN202011605157.9A CN202011605157A CN114689448B CN 114689448 B CN114689448 B CN 114689448B CN 202011605157 A CN202011605157 A CN 202011605157A CN 114689448 B CN114689448 B CN 114689448B
- Authority
- CN
- China
- Prior art keywords
- crack
- fracture
- length
- time
- wing
- 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/28—Investigating ductility, e.g. suitability of sheet metal for deep-drawing or spinning
-
- 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/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0252—Monoaxial, i.e. the forces being applied along a single axis of the specimen
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
本发明属于土木工程领域,具体地,涉及基于断裂力学与损伤力学概念,提出张开裂隙压密闭合及时效扩展的损伤断裂分析方法,步骤如下:(1)、将含多个张开裂隙岩体蠕变过程分为两个阶段:张开裂隙闭合阶段和新生裂纹蠕变扩展破坏阶段;(2)、根据裂隙压密阶段的变形和受力特性推导应变和时间表达式;(3)、推导应力强度因子与翼裂纹长度关系;(4)、推导岩石轴向应变与翼裂纹长度的关系;(5)、推导翼裂纹长度与时间的关系;(6)、推导裂纹扩展阶段轴向应变与时间的关系;(7)、进行单轴分级压缩蠕变试验,将试验结果与理论结果进行对比分析。
Description
技术领域
本发明属于土木工程领域,具体地,涉及基于断裂力学与损伤力学概念,提出一种综合考虑张开裂隙压密闭合及时效扩展的损伤断裂分析方法。
背景技术
裂隙岩体的蠕变效应已经成为影响工程稳定性的重要因素,尤其对于脆性岩石而言,其破坏时塑性变形很小,因此更容易受蠕变行为影响导致失效破坏。另外,裂隙岩体在破坏过程中往往伴随着原生裂隙尖端翼裂纹产生与扩张以及新生裂纹的萌生与发展,随着工程开挖造成应力重分布或地应力发生变化,裂隙会迅速进入不稳定扩展阶段进而发生贯通破坏。
目前,对含闭合型裂隙的岩体展开研究较多,闭合裂隙在荷载作用下,裂隙尖端会产生翼裂纹,随着荷载进一步增加,翼裂纹扩展并伴随着新的微裂纹生成和扩展,直至裂纹连通,岩体因无法继续承受外部荷载而破坏,渐进破坏过程中不存在裂隙压密的过程。但是,对于具有多个张开型裂隙的岩体,在压缩过程中裂隙之间的孔隙可以被压密实。我们要想完全掌握工程中裂隙岩体的蠕变破坏及裂纹扩展机理,研究张开型裂隙的闭合阶段是非常必要的。
为此,本发明基于损伤力学、断裂力学和流变力学的基本理论,将含多个张开裂隙岩体的蠕变破坏过程分为两个阶段:考虑时间效应的张开裂隙闭合阶段和新生裂纹蠕变扩展破坏阶段,然后分别建立可以分析两个阶段蠕变机理的理论模型。
目前国内外相关的裂隙压密与裂隙扩展理论研究现状如下:
1、《Investigation ofMacroscopic Brittle Creep Failure Caused byMicrocrack Growth Under Step Loading and Unloading in Rocks》一文,通过理论与试验结合,在损伤力学和断裂力学的基础上,推导出脆性岩石在三轴压缩情况下应变和裂纹长度、时间的关系式,并验证了理论模型的准确性。(参见《Rock Mechanics and RockEngineering》,2016,49:2581–2593.Xiao-zhao Li,Zhu-shan Shao.);
2、《A new method to model the non-linear crack closure behavior ofrocks underuniaxial compression》一文通过一种新的DEM方法用于模拟岩石裂纹中的闭合行为,并讨论了裂纹强度、裂纹间隙和裂纹方向三个因素对裂纹闭合行为的影响。(参见《International Journal ofRock Mechanics and Mining Sciences》,2018,112:171-183.Pei-Qi Ji,Xiao-Ping Zhang,Qi Zhang.)。
发明内容
为克服现有技术的不足,本发明基于断裂力学与损伤力学概念,提出一种综合考虑张开裂隙压密闭合及时效扩展的损伤断裂分析方法。为实现上述目的,本方法采用下述技术方案:
基于断裂力学与损伤力学概念,对张开裂隙压密闭合及时效扩展阶段的研究步骤如下:
(1)、将含多个张开裂隙的岩体蠕变过程分为两个阶段:张开裂隙闭合阶段和新生裂纹蠕变扩展破坏阶段;
(2)、根据裂隙压密阶段的变形和受力特性推导应变和时间表达式;
(3)、推导裂隙扩展阶段应力强度因子与翼裂纹长度关系;
(4)、推导岩石轴向应变与翼裂纹长度的关系;
(5)、推导翼裂纹长度与时间的关系;
(6)、推导裂纹扩展阶段轴向应变与时间的关系;
(7)、进行单轴分级压缩蠕变试验,将试验结果与理论结果进行对比验证。
基于上述,其具体过程如下:
1、如图1所示,含有张开裂隙的岩体受到外部荷载产生压缩变形时,其渐进破坏所经历的过程分为四个阶段:原生裂缝闭合过程、新生微裂纹稳定扩展过程、不稳定扩展过程和破坏过程。
本发明按照岩石所处的应力水平将裂隙的发展状态分为两个主要阶段,分别为裂隙闭合阶段和裂隙扩展阶段。
2、根据裂隙压密阶段的变形和受力特性推导应变和时间表达式。
理论模型中,设裂隙长度为L,宽度为b,长方体岩石试样中的裂隙数量为n,试样的长度为W,宽度为I,高度为H,体积为V,轴向荷载σ1均匀地施加在试样上下表面,则裂隙闭合阶段的轴向应变(体积应变)表达式如下:
式中:ε1表示岩石的轴向应变,ξ为裂隙闭合应变率。
随着加载的进行,具体的应变率与时间的表达式如下:
3、推导裂隙扩展阶段应力强度因子与翼裂纹长度关系。
本发明研究的裂纹类型为Ι型裂纹,Ι型裂纹的应力强度因子在单轴压缩情况中表达式如下:
μ表示裂隙平面间的摩擦系数,β表示材料常数,σ1表示轴向应力,a表示翼裂纹长度。
当不考虑单轴压缩过程中裂纹的亚临界扩展阶段时,轴向应力σ1与翼裂纹长度a的关系表达式如下:
式中:KΙC表示裂纹断裂韧度。
4、推导岩石的轴向应变与翼裂纹长度变化的关系。
本发明研究的模型中含有多个不连续裂隙,基于岩土体中原生裂隙的尺寸及密度,由损伤力学的理论,岩石轴向应变与翼裂纹的长度关系表达式如下:
式中:Nv表示单位体积内原生裂隙的数量,描述裂隙的密度对损伤变量的贡献,本发明中Nv是裂隙数量与岩石试样体积的比值,m和ε0表示与材料相关的常量参数。
5、推导翼裂纹长度与时间的关系。
岩体受力变形过程中,当裂纹尖端强度因子KΙ小于裂纹断裂韧性度KΙC时,岩体中的裂纹会随外部载荷加大而稳定扩展,此过程称为亚临界裂纹扩展阶段。考虑亚临界裂纹扩展阶段的作用,翼裂纹长度与时间的关系表达式如下:
式中:υ表示的是翼裂纹扩展速度,λ表示应力腐蚀经验指数,这两个参数通过亚临界裂纹扩展试验获取。
6、推导岩石轴向应变与时间的关系。
对翼裂纹扩展速率积分求得翼裂纹长度随时间变化的表达式,综合岩石翼裂纹长度与时间的关系式和岩石轴向应变与翼裂纹长度关系式得到岩石轴向应变与时间的关系:
7、对岩石进行单轴分级压缩蠕变试验,将试验结果与理论结果对比。
本发明建立的理论模型研究裂隙角度和裂隙长度对蠕变行为的影响。采用伺服控制的单轴加载装置,以控制变量的方法,对岩样进行分级压缩蠕变试验。得到的应变时间关系曲线,通过本发明中理论模型拟合试验结果,验证本发明中模型的正确性与实用性。
相对于现有技术,本发明具有如下有益效果:
1、本发明通过基于损伤力学和断裂力学的理论计算和推导,引入裂隙闭合应变率的概念,建立了翼裂纹长度与应力强度的关系及翼裂纹长度与时间关系,推导了轴向应变与时间的关系表达式,研究了张开型裂隙闭合阶段以及新生裂纹扩展阶段的理论模型,为预防工程灾害和对裂隙岩体中裂缝扩展的预测提供了理论指导。
2、该发明方法可广泛应用于土木、水电、采矿、能源等领域的裂隙岩体力学特性的研究,应用领域广泛。
附图说明
图1为多裂隙岩体压缩过程的典型应力-应变曲线示意图;
图2为单轴压缩过程张开型裂隙岩体受力示意图;
图3为裂隙压密过程示意图;
图4为岩石试样理论结果和试验结果对比示意图。
具体实施方式
如图1所示,含有张开型裂隙的岩体受到外部荷载产生压缩变形时,其渐进破坏所经历的过程分为四个阶段:原生裂缝闭合过程、新生微裂纹稳定扩展过程、不稳定扩展过程和破坏过程。本发明按照岩石所处的应力水平将裂隙的发展状态分为两个主要阶段,分别为裂隙闭合阶段和裂隙扩展阶段,步骤如下:
1、对于裂隙压密阶段的理论分析,具体方法如下:
如图2所示,单轴压缩情况中,裂隙岩体及内部裂隙在轴向应力σ1作用下,裂隙表面的法向应力为σα和切向应力τα根据下式求得:
式中:α表示裂隙方向与水平方向的夹角。
在理论模型中,设裂隙长度为L,宽度为b,长方体岩石试样中的裂隙数量为n,试样的长度为W,宽度为I,高度为H,体积为V,轴向荷载σ1均匀地施加在试样上下表面。假设裂隙闭合阶段的应变是由裂隙状态的变化造成,裂隙闭合阶段的轴向应变表达为:
式中:ε1表示岩石的轴向应变,ξ为裂隙闭合应变率。
随着加载的进行,具体的应变率根据试验结果表达为如下形式:
式中:ξ表示裂隙闭合应变率,E1表示粘弹性模量,E0表示岩石材料的弹性模量,ω表示裂隙闭合应变率影响系数,η表示粘滞系数。
2、推导裂隙扩展阶段应力强度因子与翼裂纹长度关系,具体方法如下:
本发明研究的裂纹类型为Ι型裂纹,Ι型裂纹的应力强度因子在单轴压缩情况中表达式如下:
μ表示裂隙平面间的摩擦系数,β表示材料常数,σ1表示轴向应力,a表示翼裂纹长度。
当不考虑单轴压缩过程中裂纹的亚临界扩展阶段时,轴向应力σ1与翼裂纹长度a的关系表达式如下:
式中:KΙC表示裂纹断裂韧度。
3、推导岩石的轴向应变与翼裂纹长度变化的关系,具体方法如下:
本发明研究的模型中含有多个不连续裂隙,基于岩土体中原生裂隙的尺寸及密度,由损伤力学的理论,岩石轴向应变与翼裂纹的长度关系表达式如下:
式中:Nv表示单位体积内原生裂隙的数量,可以描述裂隙的密度对损伤变量的贡献,本发明中Nv是裂隙数量与岩石试样体积的比值,m和ε0表示与材料相关的常量参数。
4、推导翼裂纹长度与时间的关系,具体方法如下:
岩体受力变形过程中,当裂纹尖端强度因子KΙ小于裂纹断裂韧性度KΙC时,岩体中的裂纹随外部载荷加大而稳定扩展,此过程称为亚临界裂纹扩展阶段。考虑亚临界裂纹扩展阶段的作用,翼裂纹长度与时间的关系表达式如下:
式中:υ表示的是翼裂纹扩展速度,λ表示应力腐蚀经验指数,这两个参数通过亚临界裂纹扩展试验获取。
本发明研究的是分级加载下裂隙扩展阶段的变形规律,即σ1是时间t的函数,c1和c2是裂纹长度a的函数,将式(8)数值积分得到翼裂纹长度随时间的变化关系。当KΙ=KΙC时,翼裂纹开始扩展,将此时的轴向应力σ1代入式(8)得到初始翼裂纹长度a0,然后将初始值代入对整个裂隙扩展阶段的积分过程得到每个时间点的翼裂纹长度值。
5、推导轴向应变与时间的关系,具体方法如下:
对翼裂纹扩展速率积分求得翼裂纹长度随时间变化的表达式,综合岩石翼裂纹长度与时间的关系式和岩石轴向应变与翼裂纹长度关系式得到岩石轴向应变与时间的关系:
6、对岩样进行单轴分级压缩蠕变试验,将试验结果与理论结果对比,具体方法如下:
本发明建立的理论模型可以研究裂隙角度和裂隙长度对岩石蠕变行为的影响。采用伺服控制的单轴加载装置,以控制变量的方法,对长150mm,宽150mm,高300mm的长方体岩样进行单轴分级压缩蠕变试验。在试验过程中,每级施加的荷载是试样单轴抗压强度的10%~20%。首先,以0.12mm/min的加载速率对试样进行第一级加载,加载过程中计算机自动记录应力应变随时间变化的数据。当荷载达到第一级的规定应力时,应力保持恒定64小时。以同样方式对试样进行第二级加载。在试样破坏前,荷载逐级增加。其中岩石试样中的裂隙采用具有一定厚度钢片制成,其裂隙角度为30°,裂隙长度为30mm。表1所示为理论模型中岩石试样参量,图4所示为岩石试样理论结果和试验结果对比图。左侧为试样通过室内试验得到的轴向应变—时间曲线与理论拟合的曲线对比图,右侧为试样轴向应力—轴向应变试验曲线与理论曲线对比图。通过应变时间关系曲线看出,由本发明中的理论模型式(9)拟合得到的应变变化规律与室内试验结果基本一致,验证了本发明中模型的正确性与实用性。
表1理论模型中的参数值
Claims (6)
1.张开裂隙压密闭合及时效扩展的损伤断裂分析方法,具体步骤如下:
(1)、将含多个张开裂隙的岩体蠕变过程分为两个阶段:张开裂隙闭合阶段和新生裂纹蠕变扩展破坏阶段;
(2)、根据裂隙压密阶段的变形和受力特性推导应变和时间表达式;
(3)、推导应力强度因子与翼裂纹长度关系;
(4)、推导轴向应变与翼裂纹长度的关系;
(5)、分析翼裂纹长度与时间的关系;
(6)、推导裂纹扩展阶段轴向应变与时间的关系;
(7)、进行单轴分级压缩蠕变试验,将试验结果与理论结果进行对比验证;
所述的张开裂隙压密闭合及时效扩展的损伤断裂分析方法,其特征在于,根据裂隙压密阶段的变形和受力特性推导应变和时间表达式,方法如下:
在理论模型中,设裂隙长度为L,宽度为b,长方体岩石试样中的裂隙数量为n,试样的长度为W,宽度为I,高度为H,体积为V,轴向荷载σ1均匀地施加在试样上下表面;假设裂隙闭合阶段的应变是由裂隙状态的变化造成,裂隙闭合阶段的轴向应变表达式如下:
式中:ε1表示岩石的轴向应变,ξ为裂隙闭合应变率;
随着加载的进行,具体的应变率根据试验结果表达如下:
4.根据权利要求1所述的张开裂隙压密闭合及时效扩展的损伤断裂分析方法,其特征在于,分析翼裂纹长度与时间的关系,方法如下:
岩体受力变形过程中,当裂纹尖端强度因子KI小于裂纹断裂韧性KIC时,岩体中的裂纹随外部载荷加大而稳定扩展,此过程称为亚临界裂纹扩展阶段;考虑亚临界裂纹扩展阶段的作用,得到翼裂纹长度与时间的关系式:
式中:υ表示的是翼裂纹扩展速度,λ表示应力腐蚀经验指数,这两个参数通过亚临界裂纹扩展试验获取;
权利要求1研究的是分级加载下裂隙扩展阶段的变形规律,即σ1是时间t的函数,c1和c2是裂纹长度a的函数,将式(6)进行数值积分得到翼裂纹长度随时间的变化关系;当KI=KIC时,翼裂纹开始扩展,将此时的轴向应力σ1代入式(6),得到初始翼裂纹长度a0,然后将初始值代入对整个裂隙扩展阶段的积分过程得到每个时间点的翼裂纹长度值。
6.根据权利要求1所述的张开型裂隙压密闭合及时效扩展的损伤断裂分析方法,其特征在于,进行单轴分级压缩蠕变试验,将试验结果与理论结果进行对比验证,方法如下:
权利要求1建立的理论模型研究裂隙角度和裂隙长度对蠕变行为的影响;采用伺服控制的单轴加载装置,以控制变量的方法,对含裂隙的长方体岩样进行分级压缩蠕变试验,得到应变时间关系曲线,通过权利要求1建立的理论模型拟合试验结果,验证权利要求1建立的模型的正确性与实用性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011605157.9A CN114689448B (zh) | 2020-12-30 | 2020-12-30 | 张开裂隙压密闭合及时效扩展的损伤断裂分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011605157.9A CN114689448B (zh) | 2020-12-30 | 2020-12-30 | 张开裂隙压密闭合及时效扩展的损伤断裂分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114689448A CN114689448A (zh) | 2022-07-01 |
CN114689448B true CN114689448B (zh) | 2023-03-24 |
Family
ID=82132293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011605157.9A Active CN114689448B (zh) | 2020-12-30 | 2020-12-30 | 张开裂隙压密闭合及时效扩展的损伤断裂分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114689448B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9305121B2 (en) * | 2010-06-28 | 2016-04-05 | Exxonmobil Upstream Research Company | Method and system for modeling fractures in ductile rock |
CN102928512B (zh) * | 2012-11-05 | 2015-01-21 | 中国科学院武汉岩土力学研究所 | 一种岩石时效劣化过程的测试方法 |
CN106599476B (zh) * | 2016-12-15 | 2018-03-09 | 长江水利委员会长江科学院 | 岩体二维细观时效破裂幂函数型模型方法 |
RU2667316C1 (ru) * | 2017-10-25 | 2018-09-18 | Российская Федерация в лице Министерства промышленности и торговли Российской Федерации (Минпромторг России) | Способ определения коэффициентов интенсивности напряжений для трещин |
CN108509711B (zh) * | 2018-03-29 | 2020-03-17 | 中国地质大学(北京) | 平面压剪应力下非贯通裂隙岩体损伤本构模型的构建方法 |
CN111159794A (zh) * | 2018-11-08 | 2020-05-15 | 中国石油大学(华东) | 多裂隙类岩石试样力学性质的几何损伤流变分析方法 |
CN109902444B (zh) * | 2019-03-29 | 2020-08-18 | 中南大学 | 一种预测翼裂纹扩展轨迹的方法 |
CN111220482A (zh) * | 2020-03-05 | 2020-06-02 | 中国石油大学(北京) | 岩石微观裂缝扩展模拟方法、装置、设备及存储介质 |
-
2020
- 2020-12-30 CN CN202011605157.9A patent/CN114689448B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114689448A (zh) | 2022-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cheng et al. | Multi-peak deformation behavior of jointed rock mass under uniaxial compression: Insight from particle flow modeling | |
Huang et al. | Deformation response of roof in solid backfilling coal mining based on viscoelastic properties of waste gangue | |
Dai et al. | Loading-rate-dependent progressive fracturing of cracked chevron-notched Brazilian disc specimens in split Hopkinson pressure bar tests | |
Ma et al. | Experimental investigation of the mechanical properties of rock salt under triaxial cyclic loading | |
Shan et al. | Study on the triaxial unloading creep mechanical properties and damage constitutive model of red sandstone containing a single ice-filled flaw | |
Yang et al. | Numerical study on failure behavior of brittle rock specimen containing pre-existing combined flaws under different confining pressure | |
Wang et al. | Modelling of crack propagation in rocks under SHPB impacts using a damage method | |
Meng et al. | Effects of size and strain rate on the mechanical behaviors of rock specimens under uniaxial compression | |
Huang et al. | Experimental and numerical study on loading rate effects of rock-like material specimens containing two unparallel fissures | |
Wang et al. | Creep properties and constitutive model for salt rock subjected to uniaxial trapezoidal cyclic loading | |
Xia et al. | Simulation of rock deformation and mechanical characteristics using clump parallel-bond models | |
Zhao et al. | Mathematical methods to unloading creep constitutive model of rock mass under high stress and hydraulic pressure | |
Zhao et al. | Deformation and damage evolution of rock salt under multilevel cyclic loading with constant stress intervals | |
Yang et al. | New crack initiation model for open-flawed rock masses under compression–shear stress | |
Zhang et al. | A uniaxial compressive damage constitutive model for rock mass with persistent joints considering the joint shear strength | |
CN114689448B (zh) | 张开裂隙压密闭合及时效扩展的损伤断裂分析方法 | |
Zhou et al. | Micromechanics damage modeling of brittle rock failure processes under compression | |
Lin et al. | A data assimilation method for blast load prediction | |
CN113899879A (zh) | 一种深部开挖扰动下岩石破碎后的动能确定方法 | |
Han et al. | Study on hysteretic deformation characteristics and energy dissipation of granite under cyclic loading | |
Bajpai et al. | Theoretical assessment of railway ballast degradation under cyclic loading | |
Wang et al. | Dynamic shear failure of freeze-thawed Tibet hornfels subjected to multilevel cyclic shear (MLCS) loads: insights into structural dependent failure characteristics | |
Chu et al. | Acoustic emission parameters and energy dissipation law of cyclic load and unload damage of dihydrate gypsum under different loading rate | |
Gao et al. | Crack closure and initiation stresses of coal subjected to thermo-gas-mechanical coupling | |
Zhang et al. | Dynamic Shear Fracture Behaviour of Granite Under Axial Static Pre-force by 3D High-Speed Digital Image Correlation |
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 |