CN112084605A - 掘进机刀具紧固状态实时监测方法 - Google Patents

掘进机刀具紧固状态实时监测方法 Download PDF

Info

Publication number
CN112084605A
CN112084605A CN202010996302.4A CN202010996302A CN112084605A CN 112084605 A CN112084605 A CN 112084605A CN 202010996302 A CN202010996302 A CN 202010996302A CN 112084605 A CN112084605 A CN 112084605A
Authority
CN
China
Prior art keywords
bolt
cutter
sigma
fastening state
strain sensor
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.)
Granted
Application number
CN202010996302.4A
Other languages
English (en)
Other versions
CN112084605B (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.)
Chian Railway 14th Bureau Group Corp Tunnel Engineering Co ltd
Central South University
China Railway Construction Heavy Industry Group Co Ltd
Original Assignee
Chian Railway 14th Bureau Group Corp Tunnel Engineering Co ltd
Central South University
China Railway Construction Heavy Industry Group Co Ltd
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 Chian Railway 14th Bureau Group Corp Tunnel Engineering Co ltd, Central South University, China Railway Construction Heavy Industry Group Co Ltd filed Critical Chian Railway 14th Bureau Group Corp Tunnel Engineering Co ltd
Priority to CN202010996302.4A priority Critical patent/CN112084605B/zh
Publication of CN112084605A publication Critical patent/CN112084605A/zh
Application granted granted Critical
Publication of CN112084605B publication Critical patent/CN112084605B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/10Making by using boring or cutting machines
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Details Of Spanners, Wrenches, And Screw Drivers And Accessories (AREA)

Abstract

本发明提供了一种掘进机刀具紧固状态实时监测方法,包括:步骤1,将应变传感器安装在刀具联接螺栓上,通过应变传感器监测刀具预紧过程和运行过程中的螺栓应变量;步骤2,将应变传感器监测到的螺栓应变量传输并储存在计算机中;步骤3,计算机通过螺栓轴向拉力计算模型和螺栓紧固状态判断模型,对刀具紧固状态进行实时监测。本发明根据预紧过程和运行过程的应变及轴向拉力,提出了紧固状态判断模型,最终输出实时测量的轴向拉力和连接状态的判断;本发明能够准确、定量、科学的评价联接螺栓的轴向拉力变化,解决了掘进机运行过程刀具紧固状态实时监测问题。

Description

掘进机刀具紧固状态实时监测方法
技术领域
本发明涉及隧道掘进机技术领域,特别涉及一种掘进机刀具紧固状态实时监测方法。
背景技术
伴随着现代化建设的不断深入,我国城市化进程进入快速发展模式,城市地下空间的利用和开发正成为现代化建设和可持续发展的重中之重;其中隧道掘进机广泛应用于地铁隧道施工,而刀盘刀具更是掘进机重要掘进部位,掘进机运行过程中,刀具联接是否紧固对掘进机的工作和效率影响很大。掘进机在运行过程中,会产生很大的变载荷,由于刀具与刀座通过螺栓紧固联接,剧烈的振动将对刀具联接螺栓的紧固状态产生很大影响。
目前国内外对掘进机刀具紧固状态实时监测的方法很少,普遍采用停机检查的方式来检查刀具联接的紧固状态,这种方法严重影响了施工进度和效率。
发明内容
本发明提供了一种掘进机刀具紧固状态实时监测方法,其目的是为了利用应变传感器及无线传输获得刀具联接螺栓应变实时监测数据,通过刀具联接螺栓应力和轴向拉力计算模型、螺栓紧固状态判断模型,实现对刀具紧固状态的实时监测。
为了达到上述目的,本发明的实施例提供了一种掘进机刀具紧固状态实时监测方法,包括:
步骤1,将应变传感器安装在刀具联接螺栓上,通过应变传感器监测刀具预紧过程和运行过程中的螺栓应变量;
步骤2,将应变传感器监测到的螺栓应变量传输并储存在计算机中;
步骤3,计算机通过螺栓轴向拉力计算模型和螺栓紧固状态判断模型,对刀具紧固状态进行实时监测。
其中,所述应变传感器包括一个与轴线平行的轴线应变传感器和一个与轴线垂直的垂直应变传感器,所述轴线应变传感器用于测量轴向应变值ε1,所述垂直应变传感器用于测量轴向垂直应变值ε2
其中,步骤1中,所述刀具连接螺栓通过预定的力矩扳手预紧,通过所述力矩扳手得到螺栓的计算轴向拉力
Figure BDA0002692715560000021
式中:T-扳手拧紧力矩,λ-螺纹升角,ρv-当量摩擦角,d2-螺纹中径,r-螺母和垫圈间接触面积的平均面积,μs-螺母和垫圈间的摩擦系数。
其中,步骤3中所述螺栓轴向拉力计算模型包括:
螺栓轴向应力值σx=E*(ε1+με2)/(1-μ2),测量轴向拉力FN=σx*πd1 2/4;
式中:E-螺栓的弹性模量,μ-螺栓的泊松比,d1-螺栓小径。
其中,根据螺栓的所述计算轴向拉力和测量轴向拉力得到修正系数K=F/FN
其中,步骤3中所述螺栓紧固状态判断模型包括:
运行过程中的螺栓应力σ=K*σx*[Cb/(Cb+Cm)],所述螺栓应力σ同时满足小于许用应力[σ]和许用应力幅[σa],要求σ≤[σ],[σ]=σs/S;
式中:Cb-螺栓的刚度,Cm-被连接件的刚度;σs-螺栓材料的屈服极限,S-安全系数;
运行过程中的螺栓应力幅σa=(K*σx/2)*[Cb/(Cb+Cm)]≤[σa];
所述许用应力幅[σa]=(A*Kt*Ku-1t)/(Kσ*Sa);
式中:A-尺寸系数,Ku-受力不均匀系数,Kt-螺纹制造工艺系数,Kσ-应力集中系数,Sa-安全系数,σ-1t-抗疲劳极限。
其中,步骤3中所述螺栓紧固状态判断模型还包括:
Fmin=[F-(2*[σa]*πd1 2)/4]/K;
当所述测量轴向拉力符合Fmin≤FN≤F,则判断螺栓紧固,当所述测量轴向拉力符合0≤FN≤Fmin,则判断螺栓松动。
本发明的上述方案有如下的有益效果:
本发明所述的掘进机刀具紧固状态实时监测方法根据预紧过程和运行过程的应变及轴向拉力,提出了紧固状态判断模型,最终输出实时测量的轴向拉力和连接状态的判断;本发明能够准确、定量、科学的评价联接螺栓的轴向拉力变化,解决了掘进机运行过程刀具紧固状态实时监测问题。
附图说明
图1为本发明的掘进机刀具紧固状态实时监测方法的流程图;
图2为本发明的掘进机刀具紧固状态实时监测方法的刀具结构示意图;
图3为本发明的掘进机刀具紧固状态实时监测方法的联接螺栓应变传感器安装示意图。
【附图标记说明】
1a-滚刀;1b-刀座楔块;1c-滚刀刀座;1d-刀座压块;1e-联接螺栓;1f-刀座C形块;2a-轴线应变传感器;2b-垂直应变传感器。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明针对现有的停机检查刀具紧固状态影响施工进度和效率的问题,提供了一种掘进机刀具紧固状态实时监测方法。
如图1所示,本发明的实施例提供了一种掘进机刀具紧固状态实时监测方法,包括:步骤1,将应变传感器安装在刀具联接螺栓上,通过应变传感器监测刀具预紧过程和运行过程中的螺栓应变量;步骤2,将应变传感器监测到的螺栓应变量传输并储存在计算机中;步骤3,计算机通过螺栓轴向拉力计算模型和螺栓紧固状态判断模型,对刀具紧固状态进行实时监测。
如图2所示,掘进机的滚刀刀具包括滚刀1a、刀座楔块1b、滚刀刀座1c、刀座压块1d、联接螺栓1e和刀座C形块1f,其中所述刀座楔块1b和刀座压块1d通过所述联接螺栓1e设置在所述滚刀刀座1c内,所述刀座C形块1f同样设置在所述滚刀刀座1c内,所述刀座C形块1f与刀座楔块1b配合将所述滚刀1c的转轴固定在所述滚刀刀座1c内;掘进机在工作过程中不停的掘进,所述滚刀刀具承受的变载荷较大,滚刀承受的载荷会传递到所述联接螺栓,对联接状态、刀具紧固状态产生影响。
如图3所示,所述应变传感器包括一个与轴线平行的轴线应变传感器2a和一个与轴线垂直的垂直应变传感器2b,所述轴线应变传感器2a用于测量轴向应变值ε1,所述垂直应变传感器2b用于测量轴向垂直应变值ε2
其中,步骤1中,所述刀具连接螺栓通过预定的力矩扳手预紧,通过所述力矩扳手得到螺栓的计算轴向拉力
Figure BDA0002692715560000041
式中:T-扳手拧紧力矩,λ-螺纹升角,ρv-当量摩擦角,d2-螺纹中径,r-螺母和垫圈间接触面积的平均面积,μs-螺母和垫圈间的摩擦系数。
其中,步骤3中所述螺栓轴向拉力计算模型包括:
螺栓轴向应力值σx=E*(ε1+με2)/(1-μ2),测量轴向拉力FN=σx*πd1 2/4;
式中:E-螺栓的弹性模量,μ-螺栓的泊松比,d1-螺栓小径。
其中,根据螺栓的所述计算轴向拉力和测量轴向拉力得到修正系数K=F/FN
本发明上述实施例所述的掘进机刀具紧固状态实时监测方法,由于在预紧过程中利用力矩扳手通过扭矩得到了所述联接螺栓的计算轴向拉力F,同时在掘进机运行过程中通过应变传感器检测到了所述联接螺栓的应变值进而获得了所述联接螺栓的实时测量轴向拉力F,从而进一步得出计算值和测量值的修正系数K=F/FN
其中,步骤3中所述螺栓紧固状态判断模型包括:
运行过程中的螺栓应力σ=K*σx*[Cb/(Cb+Cm)],所述螺栓应力σ同时满足小于许用应力[σ]和许用应力幅[σa],要求σ≤[σ],[σ]=σs/S;
式中:Cb-螺栓的刚度,Cm-被连接件的刚度;σs-螺栓材料的屈服极限,S-安全系数;
运行过程中的螺栓应力幅σa=(K*σx/2)*[Cb/(Cb+Cm)]≤[σa];
所述许用应力幅[σa]=(A*Kt*Ku-1t)/(Kσ*Sa);
式中:A-尺寸系数,Ku-受力不均匀系数,Kt-螺纹制造工艺系数,Kσ-应力集中系数,Sa-安全系数,σ-1t-抗疲劳极限。
本发明上述实施例所述的掘进机刀具紧固状态实时监测方法,掘进机在运行过程中,刀具联接螺栓所受轴向载荷为变载荷,因此联接螺栓需要同时满足小于许用应力[σ]和许用应力幅[σa]。
其中,步骤3中所述螺栓紧固状态判断模型还包括:
Fmin=[F-(2*[σa]*πd1 2)/4]/K;
当所述测量轴向拉力符合Fmin≤FN≤F,则判断螺栓紧固,当所述测量轴向拉力符合0≤FN≤Fmin,则判断螺栓松动。
本发明上述实施例所述的掘进机刀具紧固状态实时监测方法,由于掘进机刀具运行过程只受到压力,故取预紧时通过扭矩得到的计算轴向拉力F为最大拉应力,同时有Fmin=[F-(2*[σa]*πd1 2)/4]/K,即当所述测量轴向拉力符合Fmin≤FN≤F,则判断螺栓紧固,当所述测量轴向拉力符合0≤FN≤Fmin,则判断螺栓松动。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (7)

1.一种掘进机刀具紧固状态实时监测方法,其特征在于,包括以下步骤:
步骤1,将应变传感器安装在刀具联接螺栓上,通过应变传感器监测刀具预紧过程和运行过程中的螺栓应变量;
步骤2,将应变传感器监测到的螺栓应变量传输并储存在计算机中;
步骤3,计算机通过螺栓轴向拉力计算模型和螺栓紧固状态判断模型,对刀具紧固状态进行实时监测。
2.根据权利要求1所述的掘进机刀具紧固状态实时监测方法,其特征在于,所述应变传感器包括一个与轴线平行的轴线应变传感器和一个与轴线垂直的垂直应变传感器,所述轴线应变传感器用于测量轴向应变值ε1,所述垂直应变传感器用于测量轴向垂直应变值ε2
3.根据权利要求1所述的掘进机刀具紧固状态实时监测方法,其特征在于,步骤1中,所述刀具连接螺栓通过预定的力矩扳手预紧,通过所述力矩扳手得到螺栓的计算轴向拉力
Figure FDA0002692715550000011
式中:T-扳手拧紧力矩,λ-螺纹升角,ρv-当量摩擦角,d2-螺纹中径,r-螺母和垫圈间接触面积的平均面积,μs-螺母和垫圈间的摩擦系数。
4.根据权利要求2所述的掘进机刀具紧固状态实时监测方法,其特征在于,步骤3中所述螺栓轴向拉力计算模型包括:
螺栓轴向应力值σx=E*(ε1+με2)/(1-μ2),测量轴向拉力FN=σx*πd1 2/4;
式中:E-螺栓的弹性模量,μ-螺栓的泊松比,d1-螺栓小径。
5.根据权利要求3或4中任一项所述的掘进机刀具紧固状态实时监测方法,其特征在于,根据螺栓的所述计算轴向拉力和测量轴向拉力得到修正系数K=F/FN
6.根据权利要求5所述的掘进机刀具紧固状态实时监测方法,其特征在于,步骤3中所述螺栓紧固状态判断模型包括:
运行过程中的螺栓应力σ=K*σx*[Cb/(Cb+Cm)],所述螺栓应力σ同时满足小于许用应力[σ]和许用应力幅[σa],要求σ≤[σ],[σ]=σs/S;
式中:Cb-螺栓的刚度,Cm-被连接件的刚度;σs-螺栓材料的屈服极限,S-安全系数;
运行过程中的螺栓应力幅σa=(K*σx/2)*[Cb/(Cb+Cm)]≤[σa];
所述许用应力幅[σa]=(A*Kt*Ku-1t)/(Kσ*Sa);
式中:A-尺寸系数,Ku-受力不均匀系数,Kt-螺纹制造工艺系数,Kσ-应力集中系数,Sa-安全系数,σ-1t-抗疲劳极限。
7.根据权利要求6所述的掘进机刀具紧固状态实时监测方法,其特征在于,步骤3中所述螺栓紧固状态判断模型还包括:
Fmin=[F-(2*[σa]*πd1 2)/4]/K;
当所述测量轴向拉力符合Fmin≤FN≤F,则判断螺栓紧固,当所述测量轴向拉力符合0≤FN≤Fmin,则判断螺栓松动。
CN202010996302.4A 2020-09-21 2020-09-21 掘进机刀具紧固状态实时监测方法 Active CN112084605B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010996302.4A CN112084605B (zh) 2020-09-21 2020-09-21 掘进机刀具紧固状态实时监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010996302.4A CN112084605B (zh) 2020-09-21 2020-09-21 掘进机刀具紧固状态实时监测方法

Publications (2)

Publication Number Publication Date
CN112084605A true CN112084605A (zh) 2020-12-15
CN112084605B CN112084605B (zh) 2022-06-10

Family

ID=73739321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010996302.4A Active CN112084605B (zh) 2020-09-21 2020-09-21 掘进机刀具紧固状态实时监测方法

Country Status (1)

Country Link
CN (1) CN112084605B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006053822A (ja) * 2004-08-13 2006-02-23 Komatsu Ltd 締結された構造物の解析方法
US20070078627A1 (en) * 2005-10-05 2007-04-05 Denso Corporation Screw tightening designing method and screw tightening designing apparatus
US20100076739A1 (en) * 2008-09-25 2010-03-25 Livermore Software Technology Corporation Method of initializing bolt pretension in a finite element analysis
CN103226151A (zh) * 2013-01-25 2013-07-31 中南大学 一种掘进机刀盘盘形滚刀群体运行状态监测系统和方法
CN107271097A (zh) * 2017-06-09 2017-10-20 北京工业大学 一种实时测量螺栓夹紧力衰减的装置及方法
CN109944594A (zh) * 2019-02-23 2019-06-28 鞍山师范学院 全断面岩石掘进机连接处紧固螺栓应变监测方法
CN110345019A (zh) * 2019-08-06 2019-10-18 浙江未来技术研究院(嘉兴) 一种风力发电机叶片紧固螺栓健康状态的检测方法及系统
CN111307357A (zh) * 2019-11-12 2020-06-19 北京理工大学 一种基于微分法的螺栓紧固轴力超声检测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006053822A (ja) * 2004-08-13 2006-02-23 Komatsu Ltd 締結された構造物の解析方法
US20070078627A1 (en) * 2005-10-05 2007-04-05 Denso Corporation Screw tightening designing method and screw tightening designing apparatus
US20100076739A1 (en) * 2008-09-25 2010-03-25 Livermore Software Technology Corporation Method of initializing bolt pretension in a finite element analysis
CN103226151A (zh) * 2013-01-25 2013-07-31 中南大学 一种掘进机刀盘盘形滚刀群体运行状态监测系统和方法
CN107271097A (zh) * 2017-06-09 2017-10-20 北京工业大学 一种实时测量螺栓夹紧力衰减的装置及方法
CN109944594A (zh) * 2019-02-23 2019-06-28 鞍山师范学院 全断面岩石掘进机连接处紧固螺栓应变监测方法
CN110345019A (zh) * 2019-08-06 2019-10-18 浙江未来技术研究院(嘉兴) 一种风力发电机叶片紧固螺栓健康状态的检测方法及系统
CN111307357A (zh) * 2019-11-12 2020-06-19 北京理工大学 一种基于微分法的螺栓紧固轴力超声检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李涛: "螺栓紧固扭矩计算方法的探讨", 《炼油与化工》, 31 December 2014 (2014-12-31), pages 19 - 22 *
王小杰 等: "基于应变的螺栓轴向紧固力测量及衰减规律的研究", 《新技术新工艺》, 31 December 2017 (2017-12-31), pages 68 - 72 *

Also Published As

Publication number Publication date
CN112084605B (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
CN1210599A (zh) 刀具状态的自动监测
CN201885697U (zh) 一种预应力筋自动张拉监测装置
CN206132290U (zh) 一种盾构机滚刀测力装置
CN112084605B (zh) 掘进机刀具紧固状态实时监测方法
CN105538040A (zh) 一种加工中心床身与立柱螺栓连接状态检测方法及装置
US20090044976A1 (en) Arrangement and method for controlling drilling parameters
CN112077344B (zh) 一种车床用自适应检测夹紧力卡盘
WO2023185429A1 (zh) 定子铁心压紧状态故障检测方法、装置、发电机组监测系统及计算机可读存储介质
CN1235570A (zh) 用于制造螺纹连接的方法
CN112857641B (zh) 用于铁钻工扭矩的测量组件和铁钻工扭矩的测量装置
CN116858415A (zh) 一种高强螺栓紧固时的轴力实时监测方法
CN116379060A (zh) 一种用于检测高强螺栓松动的垫片及制作方法
CN100513735C (zh) 大口径石油套管的拧接工艺
CN110985023A (zh) 一种刀盘偏载检测装置及其检测方法
CN211042553U (zh) 一种扭矩测试设备
CN212337853U (zh) 一种随钻工具用防脱扣螺栓
CN114029893A (zh) 一种采用电动扭矩扳手的螺栓拧紧方法
CN211121951U (zh) 全自动高强螺栓复合检测仪
CN212774300U (zh) 一种用于隧道掘进机的刀柄的紧固测力装置
CN114840949B (zh) 基于负载识别的电动工具转速自适应辅助调节方法
CN219084270U (zh) 一种能够监测水轮机顶盖螺栓的智能螺栓传感器
CN217156074U (zh) 混凝土桥梁预应力精轧螺纹钢筋张拉力扭矩检测系统
CN216954934U (zh) 一种测量机床主轴拉力的仪器
CN215762886U (zh) 砻谷胶辊安装结构
CN112502726B (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