CN107421442B - 一种外部测量辅助的机器人定位误差在线补偿方法 - Google Patents
一种外部测量辅助的机器人定位误差在线补偿方法 Download PDFInfo
- Publication number
- CN107421442B CN107421442B CN201710364625.XA CN201710364625A CN107421442B CN 107421442 B CN107421442 B CN 107421442B CN 201710364625 A CN201710364625 A CN 201710364625A CN 107421442 B CN107421442 B CN 107421442B
- Authority
- CN
- China
- Prior art keywords
- robot
- coordinate system
- end effector
- laser tracker
- error
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Manipulator (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种外部测量辅助的机器人定位误差在线补偿方法,包括步骤:采用激光跟踪仪作为测量仪器,将激光跟踪仪的反射球固定在机器人末端执行器的刀具末端;分别建立机器人末端执行器坐标系及激光跟踪仪坐标系;初始化,利用机器人控制系统中任意设定的三个基准点的坐标计算机器人末端执行器坐标系到激光跟踪仪坐标系的变换矩阵;实时测量,检测任一位置时机器人末端执行器的位置误差并计算补偿向量;同步控制,在线补偿误差。采用本方法通过激光跟踪仪直接测量末端执行器的三维位置信息,进而对其定位信息进行补偿。
Description
技术领域
本发明涉及机器人,尤其涉及制造装备的定位误差补偿方法,利用外部测量仪器实现机器人末端执行器定位误差在线补偿。
背景技术
几何精度是工业机器人的重要性能指标,机器人出厂前或安装调试阶级,制造商可以提供标定服务保障机器人精度。然而,在使用过程中受到负载、重力、热效应、机械间隙、以及加工过程中的振颤等因素影响,机器人末端执行器精度逐渐降低,无法满足加工精度要求。将机器人返厂标定或请生产商现场标定,无疑会增加生产成本。因此通常采用以误差补偿的方式来提高机器人的末端精度。
当前主要有两类补偿方法:一类是离线方式,即在机器人使用之前通过标定或者建立误差补偿映射模型提高精度,这类方式建立的标定或补偿模型提前做好(通常是几个月之前),现场使用时无法改变;另一类是在线方式,即在机器人使用过程中监控误差,并实时对误差进行补偿。后一类在线对末端刀具位置进行补偿的方式显然比前一类方式的刀具定位精度高。
在线补偿的研究多数面向最小化(或补偿)导致末端刀具定位误差的各种内部机构参数和外部环境误差源。为了补偿所有可能的误差源,如几何误差、运动学误差、热误差、切削力误差等,不仅需要针对各种误差源建立对应的补偿模型,而且需要在机构内部或外部设置大量传感器监控实时状态,如温度传感器、压力传感器、激光干涉仪等。这些复杂的监控装置必然造成安装调试的耗时以及传感器之间的位置和信号干扰。
发明内容
本发明的目的在于克服已有技术的缺点,提供一种简单有效的、便于工业现场应用的机器人定位误差在线补偿方法,该方法通过激光跟踪仪直接测量末端执行器的三维位置信息,进而对其定位信息进行补偿。
本发明解决上述问题所采用的技术方案是:
本发明的一种外部测量辅助的机器人定位误差在线补偿方法,包括以下步骤:
(1)采用激光跟踪仪作为测量仪器,将激光跟踪仪的反射球固定在机器人末端执行器的刀具末端;
(2)分别建立机器人末端执行器坐标系及激光跟踪仪坐标系;
(3)初始化,利用机器人控制系统中任意设定的三个基准点的坐标计算机器人末端执行器坐标系到激光跟踪仪坐标系的变换矩阵;
(4)实时测量,检测任一位置时机器人末端执行器的位置误差并计算补偿向量,步骤如下:
(a)机器人控制器读取输入点位运动信息G代码并根据输入点位运动信息控制机器人末端运动到理论位置PC;
(b)激光跟踪仪实时测量机器人末端运动到输入点的位置信息为PL;
(c)在激光跟踪仪坐标系下计算机器人末端运动位置误差,具体过程为:利用初始化时得到的末端执行器坐标系和激光跟踪仪坐标系间变换关系,将在机器人末端执行器坐标系下的末端运动点的理论位置PC变换到激光跟踪仪坐标系下,并与该末端运动点的实际位置PL比较,计算得到机器人末端运动位置误差ΔL,如果误差大于提前设定的阈值,则记录误差值作为补偿向量,ΔL=PL-(R×PC+t),式中R为机器人末端执行器坐标系到跟踪仪坐标系变换的旋转矩阵,t为机器人末端执行器坐标系到激光跟踪仪坐标系变换的平移向量;
(d)在机器人末端执行器坐标系下补偿机器人末端运动位置误差,具体过程为:利用初始化所得的机器人末端执行器与激光跟踪仪之间的空间姿态变换关系的逆矩阵,将所得补偿向量变换回机器人末端执行器坐标系下,表示为ΔC,并补偿给机器人末端执行器的理论位置,得到该末端运动点补偿后的实际位置P′C,P′C=PC+ΔC,其中ΔC=R-1×ΔL。
(5)同步控制,在线补偿误差,控制流程如下:
(a)在机器人控制器内的实轴上,将末端运动点位的理想位置带入运动学逆解模型,求解实轴的转动参数;
(b)在机器人控制器内建立设定一个虚拟存储,作为与实轴同步的虚拟轴,将补偿向量ΔC带入运动学逆解模型,求解虚拟轴的转动参数;
(c)机器人控制器执行同步运动命令,然后将虚拟轴的运动叠加到实轴上,作为运动控制的实际输入变量,完成在线补偿。
(6)检测机器人在工作空间内的末端执行器位置误差,若不满足机器人精度需求,则重复步骤(1)至步骤(5),直至精度满足需求。
本发明的有益效果是简单有效、便于工业现场应用,该方法通过激光跟踪仪直接测量末端执行器的三维位置信息,进而对其定位信息进行补偿。因为末端执行器从测量系统得到的是绝对位置信息,所以末端执行器的定位精度是被一种可追溯的测量设备控制。
附图说明
图1是一台五自由度混联机器人的结构示意图;
图2是利用激光跟踪仪测量机器人末端位置误差示意图;
图3是本发明的在线补偿控制流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
下面结合公开号为CN104985596A的中国专利公开的一种含转动支架的五自由度混联机器人详细说明本发明的定位误差在线补偿方法。
本发明的一种外部测量辅助的机器人定位误差在线补偿方法,包括以下步骤:
(1)采用激光跟踪仪作为测量仪器,将激光跟踪仪的反射球固定在机器人末端执行器的刀具末端;
(2)分别建立机器人末端执行器坐标系及激光跟踪仪坐标系;
所述的激光跟踪仪坐标系为设备内置坐标系;
所述的机器人末端执行器坐标系建立包括以下步骤:
在机器人控制系统中任意设定三个基准点,三个基准点的坐标分别为:P1=(x1,0,0)T,P2=(0,0,0)T,P3=(0,y3,0)T,则末端执行器坐标系以P2为原点,向量Vx=P1-P2的方向为X轴方向,向量Vy=P3-P2的方向为Y轴方向,向量的方向为Z轴方向。
(3)初始化,利用机器人控制系统中设定的任意三个基准点的坐标计算机器人末端执行器坐标系到激光跟踪仪坐标系的变换矩阵,矩阵公式为:[R t],其中R为旋转矩阵,t为平移向量;其中矩阵公式中的平移向量t和旋转矩阵R的计算方法如下:
(a)控制机器人末端刀尖点分别运行到设定的三个基准点,并用激光跟踪仪测量每个基准点在跟踪仪坐标系下的位置,分别得到设定的三个基准点在激光跟踪仪坐标系下对应的点位信息为L1,L2,L3;
(b)计算平移向量:机器人末端执行器坐标系原点P2在激光跟踪仪坐标系下的坐标L2为机器人末端执行器坐标系到激光跟踪仪坐标系变换的平移向量t;
(c)计算旋转矩阵:首先,基于设定的三个基准点在激光跟踪仪坐标系下的点位信息计算三个坐标轴在激光跟踪仪坐标系下的方向向量;
其中,Nx表示X轴的方向向量,Ny表示Y轴的方向向量,Nz表示Z轴的方向向量。
然后,根据3个方向向量计算末端执行器坐标系和激光跟踪仪坐标系间坐标轴的欧拉角(θx,θy,θz),定义表示向量Nm的第i个分量,其中m={x,y,z},i={1,2,3},则
其中,Atan2表示反正切函数,sin-1表示反正弦函数;
最后,计算机器人末端执行器坐标系到跟踪仪坐标系变换的旋转矩阵R;
(4)实时测量,检测任一位置时机器人末端执行器的位置误差并计算补偿向量,步骤如下:
(a)机器人控制器读取输入点位运动信息G代码并根据输入点位运动信息控制机器人末端运动到理论位置PC;
(b)激光跟踪仪实时测量机器人末端运动到输入点的位置信息为PL;
(c)在激光跟踪仪坐标系下计算机器人末端运动位置误差,具体过程为:利用初始化时得到的末端执行器坐标系和激光跟踪仪坐标系间变换关系,将在机器人末端执行器坐标系下的末端运动点的理论位置PC变换到激光跟踪仪坐标系下,并与该末端运动点的实际位置PL比较,计算得到机器人末端运动位置误差ΔL,ΔL=PL-(R×PC+t),如果误差大于提前设定的阈值,则记录误差值作为补偿向量;
(d)在机器人末端执行器坐标系下补偿机器人末端运动位置误差,具体过程为:利用初始化所得的机器人末端执行器与激光跟踪仪之间的空间姿态变换关系的逆矩阵,将所得补偿向量变换回机器人末端执行器坐标系下,表示为ΔC,并补偿给机器人末端执行器的理论位置,得到该末端运动点补偿后的实际位置P′C,P′C=PC+ΔC,其中ΔC=R-1×ΔL。
(5)同步控制,在线补偿误差,控制流程如下:
(a)在机器人控制器内的实轴上,将末端运动点位的理想位置带入运动学逆解模型,求解实轴的转动参数;
(b)在机器人控制器内建立设定一个虚拟存储,作为与实轴同步的虚拟轴,将补偿向量ΔC带入运动学逆解模型,求解虚拟轴的转动参数;
(c)机器人控制器执行同步运动命令,然后将虚拟轴的运动叠加到实轴上,作为运动控制的实际输入变量,完成在线补偿。
(6)检测机器人在工作空间内的末端执行器位置误差,若不满足机器人精度需求,则重复步骤(1)至步骤(5),直至精度满足需求。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以作出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种外部测量辅助的机器人定位误差在线补偿方法,其特征在于包括以下步骤:
(1)采用激光跟踪仪作为测量仪器,将激光跟踪仪的反射球固定在机器人末端执行器的刀具末端;
(2)分别建立机器人末端执行器坐标系及激光跟踪仪坐标系;
(3)初始化,利用机器人控制系统中任意设定的三个基准点的坐标计算机器人末端执行器坐标系到激光跟踪仪坐标系的变换矩阵;
(4)实时测量,检测任一位置时机器人末端执行器的位置误差并计算补偿向量,步骤如下:
(a)机器人控制器读取输入点位运动信息G代码并根据输入点位运动信息控制机器人末端运动到理论位置PC;
(b)激光跟踪仪实时测量机器人末端运动到输入点的位置信息为PL;
(c)在激光跟踪仪坐标系下计算机器人末端运动位置误差,具体过程为:利用初始化时得到的末端执行器坐标系和激光跟踪仪坐标系间变换关系,将在机器人末端执行器坐标系下的末端运动点的理论位置PC变换到激光跟踪仪坐标系下,并与该末端运动点的实际位置PL比较,计算得到机器人末端运动位置误差ΔL,如果误差大于提前设定的阈值,则记录误差值作为补偿向量,ΔL=PL-(R×PC+t),式中R为机器人末端执行器坐标系到跟踪仪坐
标系变换的旋转矩阵,t为机器人末端执行器坐标系到激光跟踪仪坐标系变换的平移向量;
(d)在机器人末端执行器坐标系下补偿机器人末端运动位置误差,具体过程为:利用初始化所得的机器人末端执行器与激光跟踪仪之间的空间姿态变换关系的逆矩阵,将所得补偿向量变换回机器人末端执行器坐标系下,表示为ΔC,并补偿给机器人末端执行器的理论位置,得到该末端运动点补偿后的实际位置P′C,P′C=PC+ΔC,其中ΔC=R-1×ΔL。
(5)同步控制,在线补偿误差,控制流程如下:
(a)在机器人控制器内的实轴上,将末端运动点位的理想位置带入运动学逆解模型,求解实轴的转动参数;
(b)在机器人控制器内建立设定一个虚拟存储,作为与实轴同步的虚拟轴,将补偿向量ΔC带入运动学逆解模型,求解虚拟轴的转动参数;
(c)机器人控制器执行同步运动命令,然后将虚拟轴的运动叠加到实轴上,作为运动控制的实际输入变量,完成在线补偿;
(6)检测机器人在工作空间内的末端执行器位置误差,若不满足机器人精度需求,则重复步骤(1)至步骤(5),直至精度满足需求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710364625.XA CN107421442B (zh) | 2017-05-22 | 2017-05-22 | 一种外部测量辅助的机器人定位误差在线补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710364625.XA CN107421442B (zh) | 2017-05-22 | 2017-05-22 | 一种外部测量辅助的机器人定位误差在线补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107421442A CN107421442A (zh) | 2017-12-01 |
CN107421442B true CN107421442B (zh) | 2019-12-17 |
Family
ID=60428035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710364625.XA Active CN107421442B (zh) | 2017-05-22 | 2017-05-22 | 一种外部测量辅助的机器人定位误差在线补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107421442B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108015808B (zh) * | 2017-12-07 | 2021-04-20 | 天津大学 | 一种混联机器人的运动学标定方法 |
CN108527441B (zh) * | 2018-03-05 | 2024-02-13 | 中国计量大学 | 一种用于检测工业机器人轨迹误差的装置 |
CN108608462B (zh) * | 2018-05-07 | 2021-04-27 | 苏州尚贤医疗机器人技术股份有限公司 | 一种用于大型重载机械臂的四重精度补偿方法 |
CN108748150A (zh) * | 2018-06-04 | 2018-11-06 | 上海航天精密机械研究所 | 面向机器人加工的低成本实时补偿装置及方法 |
CN108705518A (zh) * | 2018-06-20 | 2018-10-26 | 珠海格力智能装备有限公司 | 工业机器人的处理方法、装置及系统 |
US20210291362A1 (en) * | 2018-07-16 | 2021-09-23 | Fastbrick Ip Pty Ltd | Active damping system |
KR20200091994A (ko) * | 2019-01-23 | 2020-08-03 | 삼성디스플레이 주식회사 | 레이저 가공 장치 및 레이저 가공 장치의 오차 보정 방법 |
CN110181509A (zh) * | 2019-05-14 | 2019-08-30 | 浙江树人学院(浙江树人大学) | 一种基于误差补偿的工业机器人运动控制方法 |
CN111390968B (zh) * | 2019-12-31 | 2022-05-10 | 宁波大正工业机器人技术有限公司 | 一种料盘及其点位补偿方法 |
CN111272166A (zh) * | 2020-02-27 | 2020-06-12 | 中国重汽集团济南动力有限公司 | 一种基于激光测距引导机器人的空间定位方法及系统 |
CN111366392B (zh) * | 2020-04-01 | 2022-02-08 | 伯朗特机器人股份有限公司 | 工业机器人最小定位时间的测量和计算方法 |
CN111546334B (zh) * | 2020-04-28 | 2023-02-03 | 中国科学院沈阳自动化研究所 | 一种减小轮廓误差的工业机器人在线位姿误差补偿方法 |
CN111650882A (zh) * | 2020-05-29 | 2020-09-11 | 天津大学 | 一种基于粗插补的混联机器人误差在线补偿系统及方法 |
CN112082482B (zh) * | 2020-09-09 | 2021-12-17 | 易思维(杭州)科技有限公司 | 对仅有棱边特征工件的视觉定位方法、应用及精度评估方法 |
CN112683163B (zh) * | 2020-12-03 | 2022-06-28 | 中国科学院光电技术研究所 | 一种适用于视觉测量相机的绝对位姿测量精度评估方法 |
CN113325430B (zh) * | 2021-05-10 | 2023-03-21 | 东风柳州汽车有限公司 | 末端执行器安装精度测量方法、装置、设备及存储介质 |
CN113400300B (zh) * | 2021-05-24 | 2024-05-03 | 陶建明 | 用于机器人末端的伺服系统及其控制方法 |
CN113510709B (zh) * | 2021-07-28 | 2022-08-19 | 北京航空航天大学 | 基于深度强化学习的工业机器人位姿精度在线补偿方法 |
CN113601508B (zh) * | 2021-08-16 | 2022-07-08 | 山东大学 | 一种机器人运动控制方法、系统及机器人 |
CN113778020B (zh) * | 2021-08-24 | 2023-05-26 | 哈尔滨岛田大鹏工业股份有限公司 | 一种基于向量法的毛刷补偿方法 |
CN114248270A (zh) * | 2021-12-28 | 2022-03-29 | 杭州亿恒科技有限公司 | 基于人工智能的工业机器人精度补偿方法 |
CN116673796B (zh) * | 2023-08-03 | 2023-11-10 | 成都飞机工业(集团)有限责任公司 | 一种用于机器人制孔系统的标定工具与标定方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101231749A (zh) * | 2007-12-20 | 2008-07-30 | 昆山华恒工程技术中心有限公司 | 一种工业机器人的标定方法 |
CN102062576A (zh) * | 2010-11-12 | 2011-05-18 | 浙江大学 | 基于激光跟踪测量的附加外轴机器人自动标定装置及方法 |
CN102087096A (zh) * | 2010-11-12 | 2011-06-08 | 浙江大学 | 一种基于激光跟踪测量的机器人工具坐标系自动标定装置及方法 |
CN102679925A (zh) * | 2012-05-24 | 2012-09-19 | 上海飞机制造有限公司 | 机器人定位误差测量方法 |
CN102873685A (zh) * | 2012-10-10 | 2013-01-16 | 浙江环球滤清器有限公司 | 机器人的机械误差补偿方法 |
CN103968761A (zh) * | 2014-05-28 | 2014-08-06 | 中科华赫(北京)科技有限责任公司 | 串联关节式机器人绝对定位误差校准方法及标定系统 |
CN106338990A (zh) * | 2016-08-12 | 2017-01-18 | 杭州亿恒科技有限公司 | 基于激光跟踪仪的工业机器人dh参数标定与零位标定方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8803055B2 (en) * | 2009-01-09 | 2014-08-12 | Automated Precision Inc. | Volumetric error compensation system with laser tracker and active target |
CN105739538B (zh) * | 2014-12-10 | 2019-07-23 | 上海机电工程研究所 | 基于机械臂运动装置的定位方法、定位测试装置及方法 |
-
2017
- 2017-05-22 CN CN201710364625.XA patent/CN107421442B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101231749A (zh) * | 2007-12-20 | 2008-07-30 | 昆山华恒工程技术中心有限公司 | 一种工业机器人的标定方法 |
CN102062576A (zh) * | 2010-11-12 | 2011-05-18 | 浙江大学 | 基于激光跟踪测量的附加外轴机器人自动标定装置及方法 |
CN102087096A (zh) * | 2010-11-12 | 2011-06-08 | 浙江大学 | 一种基于激光跟踪测量的机器人工具坐标系自动标定装置及方法 |
CN102679925A (zh) * | 2012-05-24 | 2012-09-19 | 上海飞机制造有限公司 | 机器人定位误差测量方法 |
CN102873685A (zh) * | 2012-10-10 | 2013-01-16 | 浙江环球滤清器有限公司 | 机器人的机械误差补偿方法 |
CN103968761A (zh) * | 2014-05-28 | 2014-08-06 | 中科华赫(北京)科技有限责任公司 | 串联关节式机器人绝对定位误差校准方法及标定系统 |
CN106338990A (zh) * | 2016-08-12 | 2017-01-18 | 杭州亿恒科技有限公司 | 基于激光跟踪仪的工业机器人dh参数标定与零位标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107421442A (zh) | 2017-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107421442B (zh) | 一种外部测量辅助的机器人定位误差在线补偿方法 | |
CN108297101B (zh) | 多关节臂串联机器人末端位姿误差检测和动态补偿方法 | |
CN108748159B (zh) | 一种机械臂工具坐标系自标定方法 | |
CN109655024B (zh) | 采用空间变换技术的位移传感器外部参数标定方法 | |
CN107020633B (zh) | 绝对机器人辅助定位方法 | |
Slamani et al. | Assessment of the positioning performance of an industrial robot | |
US7756608B2 (en) | System for calibration of an industrial robot and a method thereof | |
Wang et al. | Real-time error compensation of a three-axis machine tool using a laser tracker | |
CN109822574A (zh) | 一种工业机器人末端六维力传感器标定的方法 | |
CN111203861B (zh) | 一种机器人工具坐标系的标定方法及标定系统 | |
CN102087096A (zh) | 一种基于激光跟踪测量的机器人工具坐标系自动标定装置及方法 | |
CN111546334B (zh) | 一种减小轮廓误差的工业机器人在线位姿误差补偿方法 | |
WO2018196232A1 (zh) | 机器人和末端执行器的自动标定方法及系统 | |
CN110370271A (zh) | 工业串联机器人的关节传动比误差校准方法 | |
CN113636348A (zh) | 一种建筑安装使用的玻璃转运系统 | |
Guo et al. | A distance calibration method for kinematic parameters of serial industrial robots considering the accuracy of relative position | |
JPH0445841B2 (zh) | ||
CN113733155B (zh) | 六轴工业机器人标定装置和标定方法 | |
Wang et al. | Real-time laser tracker compensation of a 3-axis positioning system—dynamic accuracy characterization | |
JP2001038662A (ja) | 作業ロボットの校正方法 | |
Li et al. | Positional accuracy improvement of an industrial robot using feedforward compensation and feedback control | |
JP2011036956A (ja) | ロボットの精度調整方法およびロボット | |
JPS60128507A (ja) | ロボット座標系の較正方法 | |
Mayer et al. | Global kinematic calibration of a Stewart platform | |
CN116330287A (zh) | 在役精度智能诊断校准的工业机器人系统及诊断校准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |