CN113359478B - 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法 - Google Patents

单自由度直线运动平台非均匀导轨摩擦参数的辨识方法 Download PDF

Info

Publication number
CN113359478B
CN113359478B CN202110802636.8A CN202110802636A CN113359478B CN 113359478 B CN113359478 B CN 113359478B CN 202110802636 A CN202110802636 A CN 202110802636A CN 113359478 B CN113359478 B CN 113359478B
Authority
CN
China
Prior art keywords
motion platform
guide rail
degree
formula
linear motion
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
CN202110802636.8A
Other languages
English (en)
Other versions
CN113359478A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202110802636.8A priority Critical patent/CN113359478B/zh
Publication of CN113359478A publication Critical patent/CN113359478A/zh
Application granted granted Critical
Publication of CN113359478B publication Critical patent/CN113359478B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position Or Direction (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了单自由度直线运动平台非均匀导轨摩擦参数的辨识方法，包括步骤：S1、先对单自由度直线运动平台进行加装，安装力传感器和悬挂系统；S2、构建单自由度直线运动平台的等效模型，并驱动运动平台本体从机械导轨的一端运动到另一端，对运动平台本体的位移和速度进行测量；S3、在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量。本发明具有辨识容易、快捷、精度高等优点。

Description

单自由度直线运动平台非均匀导轨摩擦参数的辨识方法
技术领域
本发明涉及参数辨识的技术领域，尤其涉及到单自由度直线运动平台非均匀导轨摩擦参数的辨识方法。
背景技术
微电子制造产业是当前全球经济发展的高速增长点，也是关系国家利益和国防安全的基础性和战略性产业，其作为基础与先导的作用已为世人共识，并一直被列为历次国家科技发展规划中的重点领域之一。
由于半导体制造工艺的进步和市场对微小芯片需求的急速增长，芯片I/O密度越来越高，芯片尺寸、芯片引线间距和焊盘直径持续减小。同时，为提高生产效率，封装速度逐年递增，因而对封装设备的定位精度和运行速度、加速度提出了极高的要求。而高速运动平台作为微电子制造装备的核心组成部分，其工作精度和工作效率是提升电子制造产业生产效率和产品质量的关键。因为，对于高精度的光学仪器设备、运动平台等等高精密设备来说，机械导轨的不均匀摩擦力等非线性干扰因素(空气阻力，电磁力等粘性阻力)是影响运动平台定位精度的主要原因，这就需要更加精确的摩擦建模方法，因此发明一种可以有效辨识出运动平台非均匀导轨摩擦参数的方法至关重要。
发明内容
本发明的目的在于克服现有技术的不足，提供单自由度直线运动平台非均匀导轨摩擦参数的辨识方法。
为实现上述目的，本发明所提供的技术方案为：
单自由度直线运动平台非均匀导轨摩擦参数的辨识方法，所述单自由度直线运动平台包括运动平台本体、机械导轨、光栅尺；平台本体设置在机械导轨上，沿着机械导轨进行直线运动，并由光栅尺测量位移量，辨识方法包括以下步骤：
S1、构建单自由度直线运动平台的等效模型，并获取运动平台本体运动时的驱动力以及运动过程中的位移和速度信息；
S2、在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量。
进一步地，所述步骤S1中，构建单自由度直线运动平台的等效模型前，先对所述单自由度直线运动平台进行加装，安装力传感器和悬挂系统，通过悬挂系统使得驱动单自由度直线运动平台的驱动力不变，通过力传感器获得驱动力在单自由度直线运动平台运动过程中力的动态信息。
进一步地，所述步骤S1中，构建的单自由度直线运动平台的等效模型为：
式(1)中，ft(t)为沿导轨非均匀分布的摩擦力，c为粘性阻力系数，为速度，ü为加速度，p(t)为驱动力，m为运动平台本体确切的总质量。
进一步地，所述步骤S2中，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数fx(x)、粘性阻力系数c以及运动平台本体确切的总质量m的具体过程如下：
采用驱动力p(t)驱动运动平台本体从机械导轨的一端运动到另一端，对运动平台本体的位移和速度进行测量，分别记为u(ti)，v(ti)，i＝0,1,...,n；
得到单自由度直线运动平台的等效模型后，使用分段线性插值将p(t)和ft(t)表示为：
p(t)＝[N0(t) N1(t) … Nn(t)][p(t0) p(t1) … p)tn)]T＝N(t)P (2)
ft(t)＝[N0(t) N1(t) … Nn(t)][ft(t0) ft(t1) … ft(tn)]T＝N(t)Ft (3)
将式(2)和式(3)代入等效模型后，整理得：
在区间[t0，t]上，对式(4)两端求积分得：
在式(5)分别令t＝tj，j＝0，1，...，n，可得方程组：
mΔV+cΔU+NIFt＝NIP (6)
其中，
式(7)中，Δt为步长，在x轴上取若干离散点xj，j＝1，2，...,q，设fj，fj‘分别为xj处的摩擦力与导数，则fx(x)可表示为：
fx(x)＝S(x)Fx (8)
S(x)＝[α1(x) β1(x) α2(x) β2(x) … αq(x) βq(x)]
Fx＝[f1 f1′ f2 f′2 … fq f′q]T (9)
在式(9)令x＝u(ti)，i＝0,1,...,n，可得：
ft(ti)＝S(u(ti))Fx (10)
记
A＝[S(u(t0))T S(u(t1))T … S(u(tn))T]T (11)
则
Ft＝AFx (12)
将式(12)代入式(6)可得：
mΔV+cΔU+NIAFx＝NIP (13)
记
则式(13)变为：
Bx＝y (15)
在E次实验条件下，式(15)可扩展为：
式(17)中，B与y的下标为实验序号，式(16)为矛盾方程组，为求得最小二乘解，在方程两端同时乘以系数矩阵的转置，即
求解式(18)可得待反求参数x，其中x的第1与第2个元素分别为运动平台本体确切的总质量m和c，剩余部分为导轨离散点处的摩擦与导数组成的向量Fx；
将驱动力p(t)，运动平台本体运动过程中的位移u(ti)和速度v(ti)反馈代入到式(18)即可得到导轨摩擦函数fx(x)，粘性阻力系数c和运动平台本体确切的总质量m。
与现有技术相比，本方案原理及优点如下：
1、由于运动平台本体运动过程中的阻力除摩擦力外，还有空气阻力和电磁阻力等粘性阻力，为排除电磁阻力的干扰，本方案通过外加一个力传感器跟悬挂系统来实现运动平台本体运动跟反馈，即外加悬挂系统使得驱动平台的驱动力不变，和用力传感器获得驱动力在平台运动本体过程中力的动态信息，并利用光栅尺反馈运动平台本体在运动过程中的位移和速度信息(知道时间)。
2、本方案需要反馈得到的信息较少，只需要得到驱动力、位移以及速度即可，容易、快捷。
3、本方案在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，即可求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量，而且辨识精度高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的服务作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明单自由度直线运动平台非均匀导轨摩擦参数的辨识方法的原理流程图；
图2为本发明中涉及到的单自由度直线运动平台的等效模型示意图。
具体实施方式
下面结合具体实施例对本发明作进一步说明：
本实施例所述的单自由度直线运动平台包括运动平台本体、机械导轨、光栅尺；平台本体设置在机械导轨上，沿着机械导轨进行直线运动，并由光栅尺测量位移量。
如图1所示，本实施例所述的单自由度直线运动平台非均匀导轨摩擦参数的辨识方法，包括以下步骤：
S1、先对所述单自由度直线运动平台进行加装，安装力传感器和悬挂系统，通过悬挂系统使得驱动单自由度直线运动平台的驱动力不变，通过力传感器获得驱动力在单自由度直线运动平台运动过程中力的动态信息。
S2、构建如图2所示的单自由度直线运动平台的等效模型：
式(1)中，ft(t)为沿导轨非均匀分布的摩擦力，c为粘性阻力系数，为速度，ü为加速度，p(t)为驱动力，m为运动平台本体确切的总质量；
并采用驱动力p(t)驱动运动平台本体从机械导轨的一端运动到另一端，对运动平台本体的位移和速度进行测量，分别记为u(ti)，v(ti)，i＝0,1,...,n。
图2中，摩擦力沿导轨非均匀分布，f为坐标x的函数fx(x)，还有的是，正常情况下，运动平台本体运动过程中的阻力除摩擦力外，还有空气阻力和电磁阻力等粘性阻力，由于步骤S1的加装已经省去了电磁阻力的影响，即省去电磁阻力后的粘性阻力r＝c*v，c为粘性阻力系数，v为移动的速度。由于运动平台本体上有附加部件及负载，且有制造误差，因此运动平台本体确切总质量m未知。
S3、在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量。
本步骤具体过程如下：
得到单自由度直线运动平台的等效模型后，使用分段线性插值将p(t)和ft(t)表示为：
p(t)＝[N0(t) N1(t) … Nn(t)][p(t0) p(t1) … p(tn)]T＝N(t)P (2)
ft(t)＝[N0(t) N1(t) … Nn(t)][ft(t0) ft(t1) … ft(tn)]T＝N(t)Ft (3)
将式(2)和式(3)代入等效模型后，整理得：
在区间[t0,t]上，对式(4)两端求积分得：
在式(5)分别令t＝tj，j＝0,1,...,n，可得方程组：
mΔV+cΔU+NIFt＝NIP (6)
其中，
式(7)中，Δt为步长，在x轴上取若干离散点xj,j＝1,2,...,q，设fj，fj‘分别为xj处的摩擦力与导数，则fx(x)可表示为：
fx(x)＝S(x)Fx (8)
S(x)＝[α1(x) β1(x) α2(x) β2(x) … αq(x) βq(x)]
Fx＝[f1 f1′ f2 f′2 … fq f′q]T (9)
在式(9)令x＝u(ti)，i＝0,1,...,n，可得：
ft(ti)＝S(u(ti))Fx (10)
记
A＝[S(u(t0))T S(u(t1))T … S(u(tn))T]T (11)
则
Ft＝AFx (12)
将式(12)代入式(6)可得：
mΔV+cΔU+NIAFx＝NIP (13)
记
则式(13)变为：
Bx＝y (15)
在E次实验条件下，式(15)可扩展为：
式(17)中，B与y的下标为实验序号，式(16)为矛盾方程组，为求得最小二乘解，在方程两端同时乘以系数矩阵的转置，即
求解式(18)可得待反求参数x，其中x的第1与第2个元素分别为运动平台本体确切的总质量m和c，剩余部分为导轨离散点处的摩擦与导数组成的向量Fx；
将驱动力p(t)，运动平台本体运动过程中的位移u(ti)和速度v(ti)反馈代入到式(18)即可得到导轨摩擦函数fx(x)，粘性阻力系数c和运动平台本体确切的总质量m。
由于运动平台本体运动过程中的阻力除摩擦力外，还有空气阻力和电磁阻力等粘性阻力，为排除电磁阻力的干扰，本实施例通过外加一个力传感器跟悬挂系统来实现运动平台本体运动跟反馈，即外加悬挂系统使得驱动平台的驱动力不变，和用力传感器获得驱动力在平台运动本体过程中力的动态信息，并利用光栅尺反馈运动平台本体在运动过程中的位移和速度信息(知道时间)。
还有的是，本实施例需要反馈得到的信息较少，只需要得到驱动力、位移以及速度即可，容易、快捷。
本实施例在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，即可求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量，而且辨识精度高。
以上所述之实施例子只为本发明之较佳实施例，并非以此限制本发明的实施范围，故凡依本发明之形状、原理所作的变化，均应涵盖在本发明的保护范围内。

Claims (1)

1.单自由度直线运动平台非均匀导轨摩擦参数的辨识方法，所述单自由度直线运动平台包括运动平台本体、机械导轨、光栅尺；运动平台本体设置在机械导轨上，沿着机械导轨进行直线运动，并由光栅尺测量位移量，其特征在于，辨识方法包括以下步骤：
S1、构建单自由度直线运动平台的等效模型，并获取运动平台本体运动时的驱动力以及运动过程中的位移和速度信息；
S2、在得到单自由度直线运动平台的等效模型以及驱动力、运动过程中的位移和速度信息的情况下，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数、粘性阻力系数以及运动平台本体确切的总质量；
所述S1中，构建单自由度直线运动平台的等效模型前，先对所述单自由度直线运动平台进行加装，安装力传感器和悬挂系统，通过悬挂系统使得驱动单自由度直线运动平台的驱动力不变，通过力传感器获得驱动力在单自由度直线运动平台运动过程中力的动态信息；
所述S1中，构建的单自由度直线运动平台的等效模型为：
式(1)中，ft(t)为沿导轨非均匀分布的摩擦力，c为粘性阻力系数，为速度，/>为加速度，p(t)为驱动力，m为运动平台本体确切的总质量；
所述S2中，通过非均匀导轨摩擦参数求取算法，求得非均匀导轨摩擦函数fx(x)、粘性阻力系数c以及运动平台本体确切的总质量m的具体过程如下：
采用驱动力p(t)驱动运动平台本体从机械导轨的一端运动到另一端，对运动平台本体的位移和速度进行测量，分别记为u(ti)，v(ti)，i＝0,1,...,n；
得到单自由度直线运动平台的等效模型后，使用分段线性插值将p(t)和ft(t)表示为：
p(t)＝[N0(t) N1(t) … Nn(t)[p(t0) p(t1) … p(tn)]T＝N(t)P (2)
ft(t)＝[N0(t) N1(t) … Nn(t)][ft(t0) ft(t1) … ft(tn)]T＝N(t)Ft (3)
将式(2)和式(3)代入等效模型后，整理得：
在区间[t0，t]上，对式(4)两端求积分得：
在式(5)分别令t＝tj，j＝0，1，...，n，可得方程组：
mΔV+cΔU+NIFt＝NIP (6)
其中，
式(7)中，Δt为步长，在x轴上取若干离散点xj，j＝1，2，...，q，设fj，fj‘分别为xj处的摩擦力与导数，则fx(x)可表示为：
fx(x)＝S(x)Fx (8)
在式(9)令x＝u(ti)，i＝0，1，...，n，可得：
ft(ti)＝S(u(ti))Fx (10)
记
A＝[S(u(t0))T S(u(t1))T … S(u(tn))T]T (11)
则
Ft＝AFx (12)
将式(12)代入式(6)可得：
mΔV+cΔU+NIAFx＝NIP (13)
记
B＝[ΔV ΔU NIA]
y＝NIP (14)
则式(13)变为：
Bx＝y (15)
在E次实验条件下，式(15)可扩展为：
式(17)中，B与y的下标为实验序号，式(16)为矛盾方程组，为求得最小二乘解，在方程两端同时乘以系数矩阵的转置，即
求解式(18)可得待反求参数x，其中x的第1与第2个元素分别为运动平台本体确切的总质量m和c，剩余部分为导轨离散点处的摩擦与导数组成的向量Fx；
将驱动力p(t)，运动平台本体运动过程中的位移u(ti)和速度v(ti)反馈代入到式(18)即可得到导轨摩擦函数fx(x)，粘性阻力系数c和运动平台本体确切的总质量m。
CN202110802636.8A 2021-07-15 2021-07-15 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法 Active CN113359478B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110802636.8A CN113359478B (zh) 2021-07-15 2021-07-15 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110802636.8A CN113359478B (zh) 2021-07-15 2021-07-15 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法

Publications (2)

Publication Number Publication Date
CN113359478A CN113359478A (zh) 2021-09-07
CN113359478B true CN113359478B (zh) 2023-07-25

Family

ID=77539587

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110802636.8A Active CN113359478B (zh) 2021-07-15 2021-07-15 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法

Country Status (1)

Country Link
CN (1) CN113359478B (zh)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101226068B (zh) * 2008-02-01 2010-04-07 西安电子科技大学 动态摩擦参数测试系统及测试方法
AT509857B1 (de) * 2010-08-05 2011-12-15 Trumpf Maschinen Austria Gmbh Biegepresse mit einer werkteil-positionier- vorrichtung sowie ein verfahren zum betrieb
JP5996127B2 (ja) * 2013-12-06 2016-09-21 三菱電機株式会社 摩擦同定方法および摩擦同定装置
CN104915498B (zh) * 2015-06-08 2018-02-23 广东工业大学 基于模型识别与等效简化的高速平台运动参数自整定方法
CN205004938U (zh) * 2015-10-14 2016-01-27 南京智御恒品智能科技有限公司 一种直线电机
CN106002312B (zh) * 2016-06-29 2018-01-23 广东工业大学 一种单驱动刚柔耦合精密运动平台及其实现方法及应用
CN106230339A (zh) * 2016-08-24 2016-12-14 北京信息科技大学 一种直线电机大行程位移精度精密控制方法
CN108015564A (zh) * 2017-12-19 2018-05-11 广东工业大学 一种超声辅助克服摩擦直线运动平台
CN112462611B (zh) * 2020-11-30 2022-05-20 华中科技大学 一种精密机电系统滑动摩擦建模方法
CN112526887B (zh) * 2020-12-18 2022-06-07 北京工业大学 一种电动助力制动系统自适应摩擦补偿控制方法
CN112737451B (zh) * 2020-12-29 2022-04-01 华中科技大学 一种面向控制的永磁同步直线电机系统辨识方法

Also Published As

Publication number Publication date
CN113359478A (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
US11484982B2 (en) Online CNC machine tool geometric/thermal error measurement and compensation system
CN102735386B (zh) 考虑弯曲刚度的斜拉索索力数值计算方法
CN106994644B (zh) 半导体晶片的厚度分布测定系统及方法、研磨系统及研磨方法、厚度余量分布测定方法
CN101556138A (zh) 一种时栅直线位移传感器
US20110061441A1 (en) Gantry stage orthogonality error measurement method and error compensation method for position processing
CN109613889B (zh) 基于微分变换的数控机床在机测量系统综合误差补偿方法
CN113359478B (zh) 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法
CN106644253A (zh) 用于恒力磨削的三维力传感器解耦标定及滤波方法与装置
TW201440413A (zh) 磁對準系統及其對準方法
CN106695458A (zh) 一种数控机床进给系统的摩擦力建模方法
CN110470263A (zh) 一种基于梯度提升树的回转体测量系统误差补偿方法
CN106289058B (zh) 一种光栅尺定位分配精度补偿的方法
CN112880567B (zh) 一种附面层厚度测量方法
CN112003501B (zh) 干扰磁场下电机正弦误差的出力补偿方法及装置
CN114119501B (zh) 一种微铣削加工未变形切削厚度测量方法及系统
CN105700470A (zh) 一种用于减小机床伺服进给系统跟踪误差的方法
CN106989661B (zh) 一种测试机床静压导轨表面形状误差的方法
CN205748217U (zh) 一种直齿锥齿轮安装距检测装置
CN110243462B (zh) 运输测振仪振动检测方法
Aksonov et al. The motion differential characteristics estimation using incremental encoders in the CNC feedback loop
CN109241675A (zh) 一种用于并联构型主轴头的动力学标定方法
CN107471838A (zh) 一种宏微运动控制装置及控制方法
CN105674895B (zh) 非接触测量拉索非线性动应变的计算方法
CN221404239U (zh) 一种高精度线位移传感器全自动校准装置
CN113358066B (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