CN102928200A - 干涉测量光学材料均匀性的方法 - Google Patents
干涉测量光学材料均匀性的方法 Download PDFInfo
- Publication number
- CN102928200A CN102928200A CN2012104060661A CN201210406066A CN102928200A CN 102928200 A CN102928200 A CN 102928200A CN 2012104060661 A CN2012104060661 A CN 2012104060661A CN 201210406066 A CN201210406066 A CN 201210406066A CN 102928200 A CN102928200 A CN 102928200A
- Authority
- CN
- China
- Prior art keywords
- sample
- interferometer
- theta
- cos
- reflection mirror
- 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
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种干涉测量光学材料均匀性的方法,将光学材料加工成前后表面具有一定的夹角α的样品,利用干涉仪通过五次波面测量,然后计算得到样品的光学均匀性,为绝对测量结果。特别适合于待测元件超过干涉仪测试口径的大口径元件测量。
Description
技术领域
本发明涉及光学材料均匀性,特别是一种干涉测量光学材料均匀性的方法,主要用于光学材料均匀性的测量,适用于一般数字干涉仪,特别适合于测量尺寸超过干涉仪测试口径的待测元件。
背景技术
光学均匀性的定义是同一光学材料中各点折射率的不一致性,也称为折射率非均匀性,通常以材料最大折射率和最小折射率的差值来表示。
在许多光学应用领域,采用几百毫米量级的大尺寸光学元件,对材料的要求较高,需要测试光学均匀性。在先技术包括定性测量和定量测量,定量测量又分为一般测量和绝对测量,一般测量是指测试结果包含样品抛光波面误差或贴置板波面误差及测试仪器标准镜波面误差等。绝对测量是指测试结果去除样品抛光波面误差或贴置板波面误差及测试仪器标准镜波面误差。本发明是一种绝对测量的定量测量方法。
在先技术[1](参见GBT 7962.2-1987无色光学玻璃测试方法光学均匀性平行光管测试方法)中所描述的测试方法采用一对平行光管装置,其一作为准直管,其二作为望远镜,用分辨率法和星点法确定玻璃的均匀性,是一种定性测量方法,不能给出定量结果。
在先技术[2](参见GBT 7962.2-2010无色光学玻璃测试方法光学均匀性斐索平面干涉法)中所描述的测试方法采用斐索平面干涉仪测量无色光学玻璃的光学均匀性,是一种定量测量方法。样品垂直地放置在干涉仪的测试光路中,通过样品的透射波面的波峰-波谷值除以样品厚度,即为光学均匀性。测试结果中包含样品前后表面波面误差和干涉仪的两个标准镜组成的空腔波面误差,是一般测量方法。
在先技术[3](参见GBT 7962.3-2010无色光学玻璃测试方法光学均匀性全息干涉法)中所描述的测试方法是利用全息差分干涉原理,将反射干涉条纹和透射干涉条纹两组条纹记录在一张全息图上,从全息图的再现波面,求得折射率变化和厚度 变化的数值。是一种绝对测量的定量测量方法。但是要求样品放置垂直放置在测试光路中,待测元件的口径不能超过干涉仪口径。
在先技术[4](参见Johannes Schwider,R.Burow,K.-E.Elssner,R.Spolaczyk,and J.Grzanna,“Homogeneity testing by phase sampling interferometry”,APPLIED OPTICS,Vol.24,No.18,p3059,1985)中所描述的测量方法是利用干涉仪测量材料的光学均匀性,将样品垂直放置在干涉仪的测试光路中,测量样品前表面波面误差,透过材料的后表面波前误差和透过材料的标准反射镜波面误差,然后取出样品测量干涉仪的两个标准镜组成的空腔波面误差。将四次测量的结果通过公式计算得到光学均匀性。是一种绝对测量的定量测量方法。但是要求样品放置垂直放置在测试光路中,待测元件的口径不能超过干涉仪口径。
发明内容
本发明要解决的技术问题在于克服上述已有技术困难,提供一种干涉测量光学材料均匀性的方法,通过五次波面测量,计算得到样品的光学均匀性。特别适合于待测元件超过干涉仪测试口径的大口径元件测量。
本发明的技术解决方案如下:
一种干涉测量光学材料均匀性的方法,其特点在于该方法包括下列步骤:
①将待测的光学材料加工成样品,该样品的前后表面具有一定的夹角α,样品的前后表面的面形应符合干涉测量要求,一般波面误差不大于3.3微米,满足如下条件,使前后表面的反射光分开:
其中δ为干涉仪的分辨角度,n为样品平均折射率,θ为入射角。
②调整干涉仪的标准透射镜和标准反射镜平行,启动干涉仪,测量干涉仪的标准反射镜的波面误差C(x,y);
③在所述的标准透射镜和标准反射镜之间放入所述的样品,干涉仪出射光束对样品的入射角为θ,定义样品表面法线转动到入射光束时逆时针为正,测量透过样品的标准反射镜的波面误差为T1(x,y);
④重新放置所述的样品,该样品的位置与第③步中样品的位置关于干涉仪光轴对称,即干涉仪出射光束对样品的入射角为-θ,再测量透过样品的标准反射镜的波面误差为T2(x,y);
⑤移动所述的标准反射镜,使所述的标准反射镜与所述的样品的前表面的反射光垂直,测量经样品的前表面反射的标准反射镜的波面误差为A(x,y);
⑥调整所述的标准反射镜,使所述的标准反射镜与透过样品的后表面的反射光垂直,测量透过样品经过样品后表面反射的标准反射镜的波面误差为B(x,y);
⑦数据处理:
计算样品前表面的面形分布Za(x,y):
计算Za(x+Δx,y):将Za(x,y)水平移位Δx=2t·tgθ′·cosθ;
计算样品后表面的面形分布Zb(x,y):
计算θ方向均匀性分布Δn+(x,y):
计算-θ方向均匀性分布Δn-(x,y):
Δn-(x,y)=Δn+(x,y)-T1(x,y)+T2(x,y)
计算折射率均匀性Δn:
所述的干涉仪为数字干涉仪。
本发明的技术效果如下:
在先技术[1]用平行光管测量光学材料均匀性,无法得到定量结果。在先技术[2]用一般干涉测量方法,所得结果中包含样品表面波面误差和干涉仪标准镜波面误差。在先技术[3]和[4]分别利用全息差分干涉方法和干涉仪测量方法,为绝对测量方法,但是样品放置要求与入射光路垂直,限制了样品测试口径不能超过干涉仪标准镜口径。本发明采用干涉测量光学均匀性,共五步测量,样品任意角度放置在光路中,可以测量超过干涉仪标准镜口径的样品。测试结果中去除了样品前后表面波面误差 和干涉仪两个标准镜波面误差,为绝对测量结果。
本发明特别适合于待测元件超过干涉仪测试口径的大口径元件的测量。
附图说明
图1为本发明的测量方法流程示意图
图2为本发明中使用坐标系示意图
具体实施方式
图1为本发明的测量方法流程示意图,图中1为干涉仪的标准透射镜,2为干涉仪的标准反射镜,3是待测样品。本发明干涉仪测量光学材料均匀性的方法,包括下列步骤:
①将待测的光学材料加工成样品,该样品的前后表面具有一定的夹角α,样品的前后表面的面形应符合干涉测量要求,一般波面误差不大于3.3微米,满足如下条件,使前后表面的反射光分开:
其中δ为干涉仪的分辨角度,n为样品平均折射率,θ为入射角。
②调整干涉仪的标准透射镜和标准反射镜平行,启动干涉仪,如如图1(a),测量干涉仪的标准反射镜的波面误差C(x,y);
③如图1(b),在所述的标准透射镜和标准反射镜之间放入所述的样品,干涉仪出射光束对样品的入射角为θ,定义样品表面法线转动到入射光束时逆时针为正,测量透过样品的标准反射镜的波面误差为T1(x,y);
④如图1(c),重新放置所述的样品,该样品的位置与第③步中样品的位置关于干涉仪光轴对称,即干涉仪出射光束对样品的入射角为-θ,再测量透过样品的标准反射镜的波面误差为T2(x,y);
⑤如图1(d),移动所述的标准反射镜,使所述的标准反射镜与所述的样品的前表面的反射光垂直,测量经样品的前表面反射的标准反射镜的波面误差为A(x,y);
⑥如图1(e),调整所述的标准反射镜,使所述的标准反射镜与透过样品的后表面的反射光垂直,测量透过样品经过样品后表面反射的标准反射镜的波面误差为 B(x,y);
⑦数据处理:
本发明数据处理的基本原理如下:
通过五次测量,得到波面误差分布C(x,y),T1(x,y),T2(x,y),A(x,y),B(x,y)。满足以下五个方程,其中n为样品平均折射率,θ为入射角,θ′为折射入样品中的折射角,由折射定律求解,见下式。t为样品厚度,以上在测量时为已知量。
五个未知量包括第一步干涉仪的标准反射镜的波面误差C(x,y),前表面面形Za(x,y),后表面面形Zb(x,y),样品第二步放置时,入射光路上折射率变化Δn+(x,y),样品第三步放置时入射光路上折射率变化Δn-(x,y)。其中C(x,y)由直接测量得到。坐标系如图2所示,以干涉仪的标准透射镜1在成像CCD中心的对应点为原点,测量准直光出射方向为z轴,水平方向为x轴,垂直方向为y轴。前表面面形Za(x,y)和后表面面形Zb(x,y)的数值均为垂直光学元件表面方向上的表面高度变化。
A(x,y)=C(x,y)-4Za(x,y)cosθ (3)
上式中Za(x+Δx,y)为样品前表面面形二维分布Za(x,y)在水平方向移动Δx的结果,Δx等于样品前后表面反射光,由于样品厚度引起的错位,沿样品前表面水平方向移动的距离。
Δx=2t·tgθ′·cosθ (5)
数据处理过程如下:
测量值为波面误差分布C(x,y),T1(x,y),T2(x,y),A(x,y),B(x,y)。
计算水平方向移动Δx:A(x,y)和B(x,y)的水平方向长度之差。或者由公式(5)计算。
计算样品前表面面形分布
计算Za(x+Δx,y):将Za(x,y)水平移位Δx。
计算样品后表面面形分布Zb(x,y):
计算θ方向均匀性分布Δn+(x,y):
计算-θ方向均匀性分布Δn-(x,y):Δn-(x,y)=Δn+(x,y)-T1(x,y)+T2(x,y)
计算折射率均匀性
上述方程求解出未知量Za(x,y),Zb(x,y),Δn+(x,y),Δn-(x,y)。其中Δn+(x,y)和Δn-(x,y)的峰谷值为(Δn+)PV和(Δn-)PV,均为样品的最大折射率与最小折射率的差值。利用求平均值得到折射率均匀性 测量结果中去除了前后表面的影响,为绝对测量结果。
样品放置角度θ改变,不影响测试结果。但是当样品口径大于干涉仪标准镜口径时,更适用于增大放置角度,使测量光束覆盖样品全口径。
样品的前后表面加工面形应符合干涉测量要求,一般波面误差不大于3.3微米。
样品的前后表面需要加工成一定的角度α,使前后表面的反射光分开。满足的条件如下:
其中干涉仪的分辨角度为δ,n为样品平均折射率,θ为入射角。
下面是具体测量实施例:
待测样品3的尺寸为100mm×80mm,厚度t为15mm,波长λ为633nm,平均折射率为1.5,入射角θ为43.5度,折射角θ′为27.3度。第一步测试得到标准反射镜波面误差分布C(x,y),第二步测试样品与光轴45度放置,透过样品的标准反射镜波面误差分布为T1(x,y),第三步测试样品与光轴-45度放置,透过样品的标准反射镜波面误差为T2(x,y),第四步测试通过样品前表面反射,标准反射镜的波面误差为 A(x,y),第五步测试通过样品前表面反射,标准反射镜的波面误差为B(x,y)。
Δn-(x,y)=Δn+(x,y)-T1(x,y)+T2(x,y)
光学均匀性
测试步骤结构请参阅图1,图2。
Claims (2)
1.一种干涉测量光学材料均匀性的方法,采用干涉仪测量光学材料均匀性,其特征在于该方法包括下列步骤:
①将待测的光学材料加工成样品(3),该样品的前后表面具有一定的夹角α,样品的前后表面的面形应符合干涉测量要求,一般波面误差不大于3.3微米,满足如下条件,使前后表面的反射光分开:
其中δ为干涉仪的分辨角度,n为样品平均折射率,θ为入射角。
②调整干涉仪的标准透射镜(1)和标准反射镜(2)平行,启动干涉仪,测量干涉仪的标准反射镜(2)的波面误差C(x,y);
③在所述的标准透射镜(1)和标准反射镜(2)之间放入所述的样品(3),干涉仪出射光束对样品(3)的入射角为θ,定义样品表面法线转动到入射光束时逆时针为正,再测量透过样品的标准反射镜(2)的波面误差为T1(x,y);
④重新放置所述的样品(3),该样品的位置与第③步中样品(3)的位置关于干涉仪光轴对称,即干涉仪出射光束对样品(3)的入射角为-θ,再测量透过样品的标准反射镜2的波面误差为T2(x,y);
⑤移动所述的标准反射镜(2),使所述的标准反射镜(2)与所述的样品(3)的前表面的反射光垂直,测量经样品的前表面反射的标准反射镜(2)的波面误差为A(x,y);
⑥调整所述的标准反射镜(2),使所述的标准反射镜(2)与透过样品的后表面的反射光垂直,测量透过样品经过样品后表面反射的标准反射镜(2)的波面误差为B(x,y);
⑦数据处理:
计算样品前表面的面形分布Za(x,y):
计算Za(x+Δx,y):将Za(x,y)水平移位Δx=2t·tgθ′·cosθ;
计算样品后表面的面形分布Zb(x,y):
计算θ方向均匀性分布Δn+(x,y):
计算-θ方向均匀性分布Δn-(x,y):
Δn-(x,y)=Δn+(x,y)-T1(x,y)+T2(x,y)
计算折射率均匀性Δn:
2.根据权利要求1所述的干涉测量光学材料均匀性的方法,其特征在于所述的干涉仪为数字干涉仪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210406066.1A CN102928200B (zh) | 2012-10-22 | 2012-10-22 | 干涉测量光学材料均匀性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210406066.1A CN102928200B (zh) | 2012-10-22 | 2012-10-22 | 干涉测量光学材料均匀性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102928200A true CN102928200A (zh) | 2013-02-13 |
CN102928200B CN102928200B (zh) | 2015-02-18 |
Family
ID=47643053
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210406066.1A Active CN102928200B (zh) | 2012-10-22 | 2012-10-22 | 干涉测量光学材料均匀性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102928200B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217274A (zh) * | 2013-04-16 | 2013-07-24 | 北极光电(深圳)有限公司 | 检测和判定变形薄膜滤光片光特性合格区域的方法 |
CN103335982A (zh) * | 2013-06-21 | 2013-10-02 | 中国科学院上海光学精密机械研究所 | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 |
CN109406107A (zh) * | 2018-10-19 | 2019-03-01 | 中国兵器工业标准化研究所 | 红外光学材料均匀性测试的试样面形误差的控制方法 |
CN109406106A (zh) * | 2018-10-19 | 2019-03-01 | 中国兵器工业标准化研究所 | 光学材料折射率均匀性的评价方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002090946A1 (en) * | 2001-05-07 | 2002-11-14 | Asml Us, Inc. | Interferometric determination of three-dimensional refractive index distribution |
DE102004057442A1 (de) * | 2004-11-25 | 2005-04-14 | Carl Zeiss Smt Ag | Verfahren zum Messen einer optischen Eigenschaft eines Substrates und Verfahren zum Herstellen eines optischen Elements aus einem Substrat |
JP2006162453A (ja) * | 2004-12-08 | 2006-06-22 | Canon Inc | 測定方法及び装置、露光装置、並びに、デバイス製造方法 |
JP2007093498A (ja) * | 2005-09-30 | 2007-04-12 | Canon Inc | 測定方法及び装置、露光装置、並びに、デバイス製造方法 |
CN101251497A (zh) * | 2008-03-19 | 2008-08-27 | 中国科学院上海光学精密机械研究所 | 光学玻璃均匀性测试装置及其测试方法 |
CN101587012A (zh) * | 2009-06-30 | 2009-11-25 | 成都光明光电股份有限公司 | 光学玻璃光学均匀性测试装置及其测试方法 |
US20110051127A1 (en) * | 2009-01-22 | 2011-03-03 | Panasonic Corporation | Refractive index measuring apparatus |
CN102062678A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院上海光学精密机械研究所 | 大口径光学元件透射率和反射率的测量装置和测量方法 |
-
2012
- 2012-10-22 CN CN201210406066.1A patent/CN102928200B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002090946A1 (en) * | 2001-05-07 | 2002-11-14 | Asml Us, Inc. | Interferometric determination of three-dimensional refractive index distribution |
DE102004057442A1 (de) * | 2004-11-25 | 2005-04-14 | Carl Zeiss Smt Ag | Verfahren zum Messen einer optischen Eigenschaft eines Substrates und Verfahren zum Herstellen eines optischen Elements aus einem Substrat |
JP2006162453A (ja) * | 2004-12-08 | 2006-06-22 | Canon Inc | 測定方法及び装置、露光装置、並びに、デバイス製造方法 |
JP2007093498A (ja) * | 2005-09-30 | 2007-04-12 | Canon Inc | 測定方法及び装置、露光装置、並びに、デバイス製造方法 |
CN101251497A (zh) * | 2008-03-19 | 2008-08-27 | 中国科学院上海光学精密机械研究所 | 光学玻璃均匀性测试装置及其测试方法 |
US20110051127A1 (en) * | 2009-01-22 | 2011-03-03 | Panasonic Corporation | Refractive index measuring apparatus |
CN101587012A (zh) * | 2009-06-30 | 2009-11-25 | 成都光明光电股份有限公司 | 光学玻璃光学均匀性测试装置及其测试方法 |
CN102062678A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院上海光学精密机械研究所 | 大口径光学元件透射率和反射率的测量装置和测量方法 |
Non-Patent Citations (6)
Title |
---|
CHIAYU AI等: "Measurement of the inhomogeneity of a window", 《OPTICAL ENGINEERING》, vol. 30, no. 9, 30 September 1991 (1991-09-30), pages 1399 - 1404, XP000230033 * |
JOHANNES SCHWIDER等: "Homogeneity testing by phase sampling interferometry", 《APPLIED OPTICS》, vol. 24, no. 18, 15 September 1985 (1985-09-15), pages 3059 - 3061, XP002464165, DOI: doi:10.1364/AO.24.003059 * |
V.S.DOLADUGINA: "Technique for estimating wavefront deformation when monitoring the homogeneity of large blanks of optical glasses", 《JOURNAL OF THE OPTICAL TECHNOLOGY》, vol. 70, no. 2, 28 February 2003 (2003-02-28), pages 120 - 124 * |
刘旭等: "干涉法测量光学材料光学非均匀性", 《激光技术》, vol. 35, no. 2, 31 March 2011 (2011-03-31), pages 189 - 192 * |
林娟: "大口径光学玻璃光学均匀性干涉绝对测量方法", 《应用光学》, vol. 29, no. 1, 31 January 2008 (2008-01-31), pages 120 - 123 * |
郭培基等: "光学玻璃光学均匀性的绝对测量技术", 《激光杂志》, vol. 24, no. 3, 30 June 2003 (2003-06-30), pages 26 - 27 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217274A (zh) * | 2013-04-16 | 2013-07-24 | 北极光电(深圳)有限公司 | 检测和判定变形薄膜滤光片光特性合格区域的方法 |
CN103335982A (zh) * | 2013-06-21 | 2013-10-02 | 中国科学院上海光学精密机械研究所 | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 |
CN103335982B (zh) * | 2013-06-21 | 2016-05-11 | 中国科学院上海光学精密机械研究所 | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 |
CN109406107A (zh) * | 2018-10-19 | 2019-03-01 | 中国兵器工业标准化研究所 | 红外光学材料均匀性测试的试样面形误差的控制方法 |
CN109406106A (zh) * | 2018-10-19 | 2019-03-01 | 中国兵器工业标准化研究所 | 光学材料折射率均匀性的评价方法 |
WO2020078212A1 (zh) * | 2018-10-19 | 2020-04-23 | 中国兵器工业标准化研究所 | 红外光学材料均匀性测试的试样面形误差的控制方法 |
CN109406106B (zh) * | 2018-10-19 | 2020-06-30 | 中国兵器工业标准化研究所 | 光学材料折射率均匀性的评价方法 |
US10955364B2 (en) * | 2018-10-19 | 2021-03-23 | China North Standardization Center | Method for controlling surface errors of sample in homogeneity measurement of infrared optical materials |
Also Published As
Publication number | Publication date |
---|---|
CN102928200B (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103335982B (zh) | 利用波长调谐移相干涉仪测量平行平板光学均匀性的方法 | |
CN108061639B (zh) | 一种结合自适应光学技术的大动态范围、高精度相位差法波前测量仪 | |
CN102788558B (zh) | 散斑相关和散斑干涉相结合的三维变形测量系统及方法 | |
CN103454249B (zh) | 基于白光干涉的光学玻璃均匀性检测方法及装置 | |
CN102928200B (zh) | 干涉测量光学材料均匀性的方法 | |
CN105352915B (zh) | 一种折射率二维分布的动态测量方法 | |
CN104165758B (zh) | 基于斐索干涉仪的透镜焦距测量装置及方法 | |
CN110057543B (zh) | 基于同轴干涉的波面测量装置 | |
CN105823563B (zh) | 动态高分辨率波前测量装置及方法 | |
CN103063156A (zh) | 一种高温环境下双波长剪切干涉测量物体表面曲率的方法 | |
CN103471561B (zh) | 一种三维小角度测量装置及方法 | |
CN101986097B (zh) | 在球面面形干涉检测中高精度消除离焦误差及倾斜误差的方法 | |
CN101614526A (zh) | 测量厚度和折射率的双共焦方法及测量装置 | |
CN105043612A (zh) | 一种光学材料应力测量系统 | |
CN205537546U (zh) | 基于psd和楔形平晶微分干涉法的晶圆表面检测装置 | |
Park et al. | Measurement of film thickness by double-slit experiment | |
CN104655029A (zh) | 一种位相增强型薄膜厚度测量方法和系统 | |
CN204612666U (zh) | 一种位相增强型薄膜厚度测量系统 | |
Shepard et al. | Measurement of internal stress in glass articles | |
Shao et al. | Analysis of The Surface Shape Effect on Optical Homogeneity Measurement of Large Calibre Optical Materials | |
CN116295984B (zh) | 一种平面钢化玻璃内应力分布的无损检测方法 | |
CN115684086B (zh) | 一种基于双劈尖干涉模型测量溶液折射率的装置及方法 | |
Sasaki et al. | Exact measurement of flat surface profiles by object shifts in a phase-conjugate Fizeau interferometer | |
Leite et al. | Measurement of the refractive index of glass by optical metrology | |
CN205067057U (zh) | 光学设备标准检测仪 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |