CN109085203B - 一种材料动态剪切模量的测量方法及其测量装置 - Google Patents

一种材料动态剪切模量的测量方法及其测量装置 Download PDF

Info

Publication number
CN109085203B
CN109085203B CN201810742705.9A CN201810742705A CN109085203B CN 109085203 B CN109085203 B CN 109085203B CN 201810742705 A CN201810742705 A CN 201810742705A CN 109085203 B CN109085203 B CN 109085203B
Authority
CN
China
Prior art keywords
test piece
piezoelectric transducer
shear modulus
cylindrical test
measuring
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
CN201810742705.9A
Other languages
English (en)
Other versions
CN109085203A (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.)
Peking University
Original Assignee
Peking University
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 Peking University filed Critical Peking University
Priority to CN201810742705.9A priority Critical patent/CN109085203B/zh
Publication of CN109085203A publication Critical patent/CN109085203A/zh
Application granted granted Critical
Publication of CN109085203B publication Critical patent/CN109085203B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance

Landscapes

  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种材料动态剪切模量的测量方法及其测量装置。本发明采用扭转型压电换能器粘贴在被测圆柱试件的底面,测量得到导纳的一阶谐振峰对应的频率,根据公式剪切模量计算公式得到被测圆柱试件的剪切模量;本发明能够精确测出材料的动态剪切模量,且对材料本身没有任何损伤;本发明装置极其简单,易于制作,测量速度十分快捷;本发明适用于绝大多数材料剪切模量的测量,是一种通用的材料动态剪切模量的测量方法,在保证测量的无损性的同时还具有良好的精度,重复测量误差在千分之一左右;本发明具有很强的实用价值,有望进一步推动材料参数测量技术的发展。

Description

一种材料动态剪切模量的测量方法及其测量装置
技术领域
本发明涉及材料参数测量技术,具体涉及一种材料动态剪切模量的测量方法及其测量装置。
背景技术
材料的弹性模量是其基本的力学性质,对其进行准确测量具有重要而广泛的意义。对各向同性材料来说,杨氏模量、剪切模量和泊松比中只有两个是独立的。在应用中主要采用杨氏模量和剪切模量。目前已经有多种方法可以用来测量杨氏模量和剪切模量。静态模量是采用拉伸(压缩)和扭转试验机测量,测量的误差一般比较大,通常在5%左右甚至更大。动态弹性模量的测量方法相对比较准确。美国ASTM标准主要采用自由梁共振(ASTME1875-08),脉冲振动激励的方法(ASTM E1876-01),以及测量声速的方法(ASTM E 494-15)。这些方法理论上均能测量杨氏模量和剪切模量。但是它们在测量剪切模量时,均采用矩形条状样品。前两种方法在测量剪切模量时,需要根据样品的长厚比进行修正。而声速法测量剪切模量时,由于样品中传播的实际上是导波而不是体波,准确计算其传播时间相对比较困难。因此,上述这些方法在测量剪切模量时,通常会产生较大的误差,一般在百分之几的量级。
事实上,在早期的ASTM E1875标准中就指出,采用圆柱状样品基于扭转共振来测量剪切模量更为准确,而且计算公式也相对简单。然而,目前最大的困难是无法有效地激励样品的扭转共振。
发明内容
针对以上现有技术中存在的问题,本发明提出了一种材料动态剪切模量的测量方法及其测量装置。
本发明的一个目的在于提出一种材料动态剪切模量的测量装置。
本发明的材料动态剪切模量的测量装置包括:扭转型压电换能器和阻抗分析仪;其中,扭转型压电换能器包括两个尺寸完全相同的半圆环,外径为D,内径为d,厚度为h;两个半圆环均沿厚度方向极化;两个半圆环的面积为(D-d)h/2的侧面为电极面,两个半圆环的电极面相对,并且极化方向相反,固定粘接在一起组成圆环,从而形成扭转型压电换能器;被测圆柱试件的直径与扭转型压电换能器的外径一致;扭转型压电换能器固定粘接在被测圆柱试件的一个底面;扭转型压电换能器的两个电极面分别通过导线连接至阻抗分析仪;阻抗分析仪施加电压并测量被测圆柱试件的导纳随频率变化的曲线,得出导纳的一阶谐振峰对应的频率f;根据剪切模量计算公式
Figure GDA0001787002110000021
得到被测圆柱试件的剪切模量G;其中,ρM是被测圆柱试件的密度,ρP是扭转型压电换能器的密度,L是被测圆柱试件的长度。
扭转型压电换能器的材料采用压电材料;如压电陶瓷。
两个半圆环的电极面采用导电胶粘接在一起。
扭转型压电换能器与被测圆柱试件采用胶水粘接固定。
进一步,在两个半圆环的电极面上分别设置电极,两个半圆环的电极通过导电胶粘接在一起。
本发明的另一个目的在于提供一种材料动态剪切模量的测量方法。
本发明的材料动态剪切模量的测量方法,包括以下步骤:
1)将外径为D,内径为d,厚度为h的压电圆环沿厚度方向极化;
2)将极化后的圆环沿直径平均切割为两个半圆环,两个半圆环的面积为(D-d)h/2的侧面为电极面;
3)将两个半圆环的电极面相对,并且极化方向相反,固定粘接在一起组成圆环,从而形成扭转型压电换能器;
4)将扭转型压电换能器固定粘接在被测圆柱试件的一个底面;
5)扭转型压电换能器的两个电极面分别通过导线连接至阻抗分析仪;
6)阻抗分析仪施加电压,阻抗分析仪测量被测圆柱试件的导纳随频率变化的曲线,得出导纳的一阶谐振峰对应的频率f;
7)根据剪切模量计算公式:
Figure GDA0001787002110000022
得到被测圆柱试件的剪切模量G;其中,ρM是被测圆柱试件的密度,ρP是扭转型压电换能器的密度,L是被测圆柱试件的长度。
在步骤3)中,采用导电胶将两个半圆环的电极面相对粘接在一起。
在步骤4)中,采用胶水将扭转型压电换能器与被测圆柱试件粘接。
进一步包括,在两个半圆环的电极面上分别镀上电极,两个半圆环的电极通过导电胶粘接在一起。
本发明的优点:
本发明采用扭转型压电换能器粘贴在被测圆柱试件的底面,测量得到导纳的一阶谐振峰对应的频率,根据公式剪切模量计算公式得到被测圆柱试件的剪切模量;本发明能够精确测出材料的动态剪切模量,且对材料本身没有任何损伤;本发明装置极其简单,易于制作,测量速度十分快捷;本发明适用于绝大多数材料剪切模量的测量,是一种通用的材料动态剪切模量的测量方法,在保证测量的无损性的同时还具有良好的精度,重复测量误差在千分之一左右;本发明具有很强的实用价值,有望进一步推动材料参数测量技术的发展。
附图说明
图1为本发明的材料动态剪切模量的测量装置的一个实施例的示意图;
图2为本发明的材料动态剪切模量的测量装置的一个实施例的扭转型压电换能器的制备流程图;
图3为根据本发明的材料动态剪切模量的测量方法得到的被测圆柱试件的导纳随频率变化的曲线;
图4为本发明的扭转型压电换能器与被测圆柱试件粘贴在一起后的机电等效图。
具体实施方式
下面结合附图,通过具体实施例,进一步阐述本发明。
如图1所示,本实施例的材料动态剪切模量的测量方法,包括以下步骤:
1)在圆环1的上下表面分别制作电极,将外径为D,内径为d,厚度为h的压电陶瓷圆环沿厚度方向极化,然后除去电极,如图2(a)所示,图中的箭头方向为极化方向;
2)将极化后的圆环沿直径平均切割为两个半圆环,两个半圆环的面积为(D-d)h/2的侧面为电极面2,并在两个半圆环的电极面上喷制电极,如图2(b)所示;
3)将两个半圆环的电极面相对,并且极化方向相反,采用导电银胶固定粘接在一起组成圆环,从而形成扭转型压电换能器,如图2(c)所示;
4)将扭转型压电换能器采用502胶水固定粘接在被测圆柱试件的一个底面,被测圆柱试件的长度为L,如图2(d)所示;
5)扭转型压电换能器的两个电极面分别通过导线连接至阻抗分析仪6,如图1所示;
6)阻抗分析仪测量扭转型压电换能器与圆柱试件粘接后的导纳随频率变化的曲线,得出导纳的一阶谐振峰对应的频率f,如图3所示;
7)根据公式剪切模量计算公式:
Figure GDA0001787002110000041
得到被测圆柱试件的剪切模量G;其中,ρM是被测圆柱试件的密度,ρP是扭转型压电换能器的密度,L是被测圆柱试件的长度。
为了进一步验证本方法与装置的实用性,利用该方法分别对45号钢、铝、以及石英玻璃的圆柱试件进行测试。其中,扭转型压电换能器的材料为压电陶瓷PZT-5H,尺寸为:内径d=5mm±0.05mm;外径D=12mm±0.01mm;厚度h=1.8mm±0.01mm。密度ρP=7500kg/m3±7.5kg/m3,被测试件的直径均为D=12mm±0.01mm,其他相关参数已经列入下表。阻抗分析仪施加的电压幅值为5V;频率从0~60kHz扫描。测量的导纳曲线如图3所示,结果表明一阶谐振峰值相对于高阶谐振峰值较小,但通过对其进行局部扫频也可清晰的辨别出一阶谐振峰对应的谐振频率。最后将此谐振频率带入剪切模量计算公式:
Figure GDA0001787002110000042
即可获得被测材料的剪切模量。将三种材料测量的结果与相关资料给出的结果对比如下表:
Figure GDA0001787002110000043
可以看到,本发明的方法与装置能够精确的测量出材料的剪切模量。
最后介绍本发明的原理,即上述剪切模量计算公式的推导:
参见图1中的坐标系,首先根据压电半圆环的振动模式,认为其只有环向位移uθ,且环向位移uθ与环向电场Eθ的形式为:
Figure GDA0001787002110000051
其中ΘP是环向角位移,Vθ是环向施加的电压。
接着根据第一类压电方程:
Figure GDA0001787002110000052
以及运动方程:
Figure GDA0001787002110000053
即可得到振动方程:
Figure GDA0001787002110000054
其中γθz是剪应变;σθz是剪应力;
Figure GDA0001787002110000055
是压电陶瓷电学短路时的弹性柔顺系数;d15是压电系数;
Figure GDA0001787002110000056
是自由介电常数;Dθ是环向电位移;ρP是压电陶瓷的密度;cP是压电陶瓷电学短路时的横波波速。
再根据速度边界条件:
Figure GDA0001787002110000057
与力学边界条件:
Figure GDA00017870021100000511
即可获得压电半圆环的传输方程:
Figure GDA0001787002110000059
其中
Figure GDA00017870021100000510
代表环向角位移的导数,即环向角速度,U1与U2分别是压电半圆环在z=0与z=h处的角速度,F′1与F′2分别是压电半环在z=0与z=h处受的外力矩,dA是面积几分符号。j代表虚数单位,Z′P是对表达式的简化,kP是波数,N′是压电半圆环的机电转换系数,C′0是截止电容,I′是电流,ωP是角速度,
Figure GDA00017870021100000512
Figure GDA00017870021100000513
分别是压电陶瓷夹持介电常数与自由介电常数,其具体表达式如下:
Figure GDA0001787002110000061
扭转型压电换能器是将两个极化方向相反的半圆环共用电极面后得到的,此时圆环输出的力F1与F2以及电流I为原来F′1、F′2以及I′的两倍,而速度U1与U2以及电压Vθ不变,因此认为:当两个压电半圆环组成一个新的扭转型压电换能器时,其传输方程的形式不会改变,只是圆环相应的各参数,机电转换系数N,截止电容C0以及简化表达式ZP变为半圆环的两倍,即:N=2N′;ZP=2Z′P;C0=2C′0
参见图1中的坐标系,对被测圆柱试件做同样的分析可得到其传输方程:
Figure GDA0001787002110000062
其中U2与U3分别是被测圆柱试件在z′=0与z′=L处的角速度,F2与F3分别是被测圆柱试件在z′=0与z′=L处受的外力矩,ZM是简化表达式,kM是圆柱试件的波数,cM是圆柱试件的波速,ρM是被测圆柱试件的密度,IPM是被测圆柱试件的极惯性矩,ωM是圆柱试件的角速度。其具体表达式如下:
Figure GDA0001787002110000063
最后考虑到粘接处位移、速度、与力的连续性,可以得到:ωP=ωM=ω。即用ω统一表示,再根据扭转型压电换能器的传输方程与被测圆柱试件的传输方程即可画出它们粘接在一起后的机电等效图如图4所示。在图4中,六个等效阻抗的阻值依次分别为:
Figure GDA0001787002110000071
Figure GDA0001787002110000072
Figure GDA0001787002110000073
阻值相应地标明在对应的等效阻抗旁边。考虑到两端自由的边界条件后,即F1与F3等于0,可以得到:
Figure GDA0001787002110000074
其中
Figure GDA0001787002110000075
其中,I是电流,C0是电容,ω是角频率,N是扭转型压电换能器的机电转换系数,Z是对电路图中电阻简化的结果。
当Z=0时导纳取极大值,也即:
Figure GDA0001787002110000076
将ZP与ZM的表达式带入后得:
Figure GDA0001787002110000077
这是一个超越方程,在实验获得f后通过数值求解就能得到精确的剪切模量G。为了实际使用的方便,在试件达到一阶谐振频率时,对于绝大多数材料在L>>h的情况下:
Figure GDA0001787002110000078
由此可以得到:
Figure GDA0001787002110000079
最后需要注意的是,公布实施例的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。

Claims (8)

1.一种材料动态剪切模量的测量方法,其特征在于,所述测量方法包括以下步骤:
1)将外径为D,内径为d,厚度为h的压电圆环沿厚度方向极化;
2)将极化后的圆环沿直径平均切割为两个半圆环,两个半圆环的面积为(D-d)h/2的侧面为电极面;
3)将两个半圆环的电极面相对,并且极化方向相反,固定粘接在一起组成圆环,从而形成扭转型压电换能器;
4)将扭转型压电换能器固定粘接在被测圆柱试件的一个底面;
5)扭转型压电换能器的两个电极面分别通过导线连接至阻抗分析仪;
6)阻抗分析仪施加电压,阻抗分析仪测量被测圆柱试件的导纳随频率变化的曲线,得出导纳的一阶谐振峰对应的频率f;
7)根据剪切模量计算公式:
Figure FDA0002398165630000011
得到被测圆柱试件的剪切模量G;其中,ρM是被测圆柱试件的密度,ρP是扭转型压电换能器的密度,L是被测圆柱试件的长度。
2.如权利要求1所述测量方法,其特征在于,在步骤3)中,采用导电胶将两个半圆环的电极面相对粘接在一起。
3.如权利要求1所述测量方法,其特征在于,在步骤4)中,采用胶水将扭转型压电换能器与被测圆柱试件粘接。
4.如权利要求1所述测量方法,其特征在于,进一步包括,在两个半圆环的电极面上分别镀上电极,两个半圆环的电极通过导电胶粘接在一起。
5.一种材料动态剪切模量的测量装置,其特征在于,所述测量装置包括:扭转型压电换能器和阻抗分析仪;其中,扭转型压电换能器包括两个尺寸完全相同的半圆环,外径为D,内径为d,厚度为h;两个半圆环均沿厚度方向极化;两个半圆环的面积为(D-d)h/2的侧面为电极面,两个半圆环的电极面相对,并且极化方向相反,固定粘接在一起组成圆环,从而形成扭转型压电换能器;被测圆柱试件的直径与扭转型压电换能器的外径一致;扭转型压电换能器固定粘接在被测圆柱试件的一个底面;扭转型压电换能器的两个电极面分别通过导线连接至阻抗分析仪;阻抗分析仪施加电压并测量被测圆柱试件的导纳随频率变化的曲线,得出导纳的一阶谐振峰对应的频率f;根据剪切模量计算公式
Figure FDA0002398165630000021
得到被测圆柱试件的剪切模量G;其中,ρM是被测圆柱试件的密度,ρP是扭转型压电换能器的密度,L是被测圆柱试件的长度。
6.如权利要求5所述的测量装置,其特征在于,所述扭转型压电换能器的材料采用压电材料。
7.如权利要求5所述的测量装置,其特征在于,所述两个半圆环的电极面采用导电胶粘接在一起。
8.如权利要求5所述的测量装置,其特征在于,进一步,在两个半圆环的电极面上分别设置电极,两个半圆环的电极通过导电胶粘接在一起。
CN201810742705.9A 2018-07-09 2018-07-09 一种材料动态剪切模量的测量方法及其测量装置 Active CN109085203B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810742705.9A CN109085203B (zh) 2018-07-09 2018-07-09 一种材料动态剪切模量的测量方法及其测量装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810742705.9A CN109085203B (zh) 2018-07-09 2018-07-09 一种材料动态剪切模量的测量方法及其测量装置

Publications (2)

Publication Number Publication Date
CN109085203A CN109085203A (zh) 2018-12-25
CN109085203B true CN109085203B (zh) 2020-06-26

Family

ID=64837217

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810742705.9A Active CN109085203B (zh) 2018-07-09 2018-07-09 一种材料动态剪切模量的测量方法及其测量装置

Country Status (1)

Country Link
CN (1) CN109085203B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987595A (zh) * 2019-12-19 2020-04-10 北京大学 一种高低温环境下材料弹性模量与内耗的测量方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101046457A (zh) * 2007-05-08 2007-10-03 中国科学院上海微系统与信息技术研究所 一种应用于生物微质量检测的硅基压电薄膜传感器及制作方法
CN105865923A (zh) * 2016-05-31 2016-08-17 北京大学 一种软物质力学性能测量方法及其测量系统
CN106766580A (zh) * 2016-12-30 2017-05-31 广州芯康医疗科技有限公司 除霜装置、超低温存储箱及其除霜方法
CN106940346A (zh) * 2017-03-20 2017-07-11 北京大学 全指向性激励与接收水平剪切波的压电换能器及制备方法
CN107462634A (zh) * 2017-09-04 2017-12-12 中国计量大学 基于压电效应的土体抗剪强度原位测量方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120166095A1 (en) * 2010-12-23 2012-06-28 General Electric Company Highly selective chemical and biological sensors

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101046457A (zh) * 2007-05-08 2007-10-03 中国科学院上海微系统与信息技术研究所 一种应用于生物微质量检测的硅基压电薄膜传感器及制作方法
CN105865923A (zh) * 2016-05-31 2016-08-17 北京大学 一种软物质力学性能测量方法及其测量系统
CN106766580A (zh) * 2016-12-30 2017-05-31 广州芯康医疗科技有限公司 除霜装置、超低温存储箱及其除霜方法
CN106940346A (zh) * 2017-03-20 2017-07-11 北京大学 全指向性激励与接收水平剪切波的压电换能器及制备方法
CN107462634A (zh) * 2017-09-04 2017-12-12 中国计量大学 基于压电效应的土体抗剪强度原位测量方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Investigation on the dynamic shear modulus and damping ratio of steel slag sand mixtures;Wei Li 等;《Construction and Building Materials》;20171215;第162卷;第170-180页 *
大范围检测和监测的全向型SH波压电换能器;宦强 等;《2017远东无损检测新技术论坛》;20171230;第231-236页 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110987595A (zh) * 2019-12-19 2020-04-10 北京大学 一种高低温环境下材料弹性模量与内耗的测量方法及装置

Also Published As

Publication number Publication date
CN109085203A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
Malocha et al. Recent measurements of material constants versus temperature for langatate, langanite and langasite
CN106895803A (zh) 一种利用两个声表面波谐振器分离温度影响的装置及方法
CN114778698B (zh) 基于复合压电薄膜体声波谐振的材料弹性模量测量方法
Xie et al. A modified piezoelectric ultrasonic composite oscillator technique for simultaneous measurement of elastic moduli and internal frictions at varied temperature
CN106840056A (zh) 一种双声表面波应变传感器及其设计方法
CN109085203B (zh) 一种材料动态剪切模量的测量方法及其测量装置
Wang et al. An online viscosity measurement method based on the electromechanical impedance of a piezoelectric torsional transducer
CN110987595B (zh) 一种高低温环境下材料弹性模量与内耗的测量方法及装置
Jiang et al. Determining the optimal pre-tightening force of a sandwich transducer by measuring resonance resistance
CN105115858B (zh) 一种基于机电阻抗法的粘度计及其检测方法
Castro et al. Apparent negative mass in QCM sensors due to punctual rigid loading
Wang et al. A U-shape shear horizontal waveguide sensor for on-line monitoring of liquid viscosity
JP5225284B2 (ja) 電気機械変換素子の電気機械特性検査方法
Purohit et al. Low force measurement based on impedance of a piezo-resonator
Lee et al. Improvements in electrical properties of piezoelectric microcantilever sensors by reducing parasitic effects
Sánchez et al. A piezoelectric minirheometer for measuring the viscosity of polymer microsamples
JP4734657B2 (ja) トランスデューサ及びこのトランスデューサを備えた計測装置
Levassort et al. Characterisation of thin layers of parylene at high frequency using PZT thick film resonators
RU2730127C1 (ru) Способ неразрушающего контроля пьезопакетов
Jeyaseelan et al. Smart piezo-ceramic coating for multifunctional applications
JP6694565B2 (ja) 粘弾性特性の測定方法およびこの方法を用いた測定装置
CN115468866B (zh) 压电材料的霍普金森一维动态压缩力电特性的试验方法
Xu et al. High-precision low-power quartz tuning fork temperature sensor with optimized resonance excitation
JP3550113B2 (ja) 硬さセンサー及びそれを用いた硬さ測定方法
Kolkoori et al. Modeling of 3D electro-mechanical wave propagation in 1-3 piezocomposite based ultrasonic sensors: simulation and experimental comparison

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