CN111268182A - 空间双臂机器人约束柔顺稳定控制方法 - Google Patents

空间双臂机器人约束柔顺稳定控制方法 Download PDF

Info

Publication number
CN111268182A
CN111268182A CN202010087291.8A CN202010087291A CN111268182A CN 111268182 A CN111268182 A CN 111268182A CN 202010087291 A CN202010087291 A CN 202010087291A CN 111268182 A CN111268182 A CN 111268182A
Authority
CN
China
Prior art keywords
robot
control
compliance
space
arm
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
CN202010087291.8A
Other languages
English (en)
Other versions
CN111268182B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202010087291.8A priority Critical patent/CN111268182B/zh
Publication of CN111268182A publication Critical patent/CN111268182A/zh
Application granted granted Critical
Publication of CN111268182B publication Critical patent/CN111268182B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G4/00Tools specially adapted for use in space
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G4/00Tools specially adapted for use in space
    • B64G2004/005Robotic manipulator systems for use in space

Abstract

本发明公开了一种空间双臂机器人的约束柔顺控制方法,推导了空间机器人的动力学模型;设计了调整期望运动的双层柔顺结构,协调了对目标操作以及机械臂之间挤压的柔顺;通过设计障碍李雅普诺夫函数,实现了对期望轨迹的约束跟踪;最后通过实例验证了本发明提出的方法的有效性。本发明通过生成柔顺期望轨迹实现了目标操作以及机械臂挤压之间的柔顺,并且约束了关节跟踪期望运动的控制性能,对抓捕后的非合作目标实现了可靠的双臂柔顺稳定控制。

Description

空间双臂机器人约束柔顺稳定控制方法
技术领域
本发明涉及一种空间双臂机器人的柔顺控制方法,特别针对空间机器人抓捕非合作目标后的稳定控制任务。
背景技术
空间机器人抓取在轨目标是在轨服务任务所需的重要关键基础技术之一。在现有的研究中,研究人员通常将抓捕任务按序分为三个阶段,即抓捕前、抓捕中以及抓捕后。空间机器人在抓捕前、中两个阶段中操纵机械臂末端靠近抓捕点,并且实现机械臂末端与抓捕点的握紧。在抓捕后阶段,空间机器人的机械臂和基座协同操作,完成对目标的稳定或者重构任务。在整个抓捕任务中,抓捕前、中段的任务目标分别侧重于操作最优性的保证以及接触安全的保障。而抓捕后阶段与前两个阶段不同,抓捕后阶段任务要求空间机器人在保证机械臂末端与被抓捕目标接触安全的同时提高操作效率。相较于空间单臂机器人,空间双臂机器人能够满足大载荷操作的任务需求,适用于更多的空间非合作目标。因此,空间双臂机器人协调控制效率与操纵安全的控制方法是现阶段研究的重点。
对于惯性参数不确定的翻滚目标,根据估计参数设计的期望运动轨迹并不安全。在跟踪期望运动的过程中,末端接触可能对抓捕点造成损伤。针对机器人在抓取的过程中,机器人无法准确快速地跟踪上期望轨迹,进而导致因为跟踪的期望位移与力不统一而破坏被抓捕目标的问题,一部分学者将柔顺控制应用于抓捕后的消旋控制过程当中,以实现对翻滚目标的柔顺消旋。在实际的应用中,即使研究者们在规划优化的过程中考虑了接触安全的约束,在稳定翻滚的非合作目标时,由于目标的不确定性的影响,机械臂末端与目标的接触安全仍然无法保障。简而言之,机械臂准确地跟踪期望轨迹时,机械臂末端的接触力与力矩会因为目标的不确定性影响而违反规划过程中设计的接触约束。在已有研究成果中,研究者们通过设计基于位置的柔顺控制律,实现了空间单臂机器人对非合作目标的柔顺稳定,但对于空间双臂机器人,在对目标操作柔顺的基础上还需要考虑双臂之间挤压的柔顺,目前尚没有公开文献对这一问题进行研究。
发明内容
针对空间双臂机器人抓捕非合作目标后机械臂末端以及机械臂之间的挤压会造成目标损坏的问题,本发明提出一种针对空间双臂机器人的柔顺控制方法,根据双臂末端的接触对期望运动进行调整,采用两层柔顺结构,逐步调整操作过程中的期望轨迹,得到调整后的期望柔顺运动。本发明基于障碍李雅普诺夫函数,设计满足控制性能的跟踪控制器,通过约束关节跟踪期望轨迹的性能,实现空间双臂机器人对非合作目标的约束柔顺稳定,提高了控制方法的可靠性。
本发明的技术方案
所述一种空间双臂机器人约束柔顺稳定控制方法,其特征在于:包括以下步骤:
步骤1:建立空间双臂机器人动力学模型:
Figure BDA0002382507130000021
其中Hc为空间机器人的惯性矩阵,cc为空间机器人所受到的科式力,
Figure BDA0002382507130000022
为基座位置与机械臂关节角度组成的广义状态变量,维度为基座的6自由度加双机械臂的自由度2n,一共6+2n,n为机械臂关节个数,
Figure BDA0002382507130000023
为基座加速度与机械臂关节角加速度组成的广义变量,u为系统控制输入,维度为6+2n,Je为机械臂关节以及基座到机械臂末端的雅阁比矩阵,fe为机械臂末端测量的接触力;
步骤2:根据控制律
Figure BDA0002382507130000024
计算系统控制输入u并进行空间双臂机器人控制;其中α为中间控制量,是维度为6+2n的列向量,单个元素的表达式为:
Figure BDA0002382507130000025
xri为列向量xr的元素,xr为经过柔顺处理后的基座位置与机械臂关节角度的期望状态,维度为6+2n;esi为列向量es的元素,es=x-xr为跟踪误差;kcli=kci+xri,kcui=kci-xri,kci为列向量kc的元素,kc为设定的控制器性能约束参数,维度为6+2n,λ为设定的极小正常数,kei为列向量ke的元素,ke为设定的控制参数,维度为6+2n;KZ=diag(kz)为对角阵,kz也为设定的控制参数,维度为6+2n,
Figure BDA0002382507130000031
g为控制补偿量:
Figure BDA0002382507130000032
Figure BDA0002382507130000033
进一步的,系统控制输入u为由基座控制力及力矩以及机械臂关节控制力矩组成的列向量。
有益效果
本发明提出了一种空间双臂机器人约束柔顺稳定的控制方法,推导了空间双臂机器人的动力学模型,设计了内外双环结构对期望运动进行柔顺化处理,设计性能约束跟踪控制器,实现空间非合作目标的约束柔顺稳定,并且通过实例验证本发明能够实现利用约束跟踪控制器跟踪双层柔顺结构调整后的柔顺期望稳定轨迹。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1空间机器人系统示意图;
图2机械臂1末端接触力;
图3机械臂1末端接触力矩;
图4机械臂2末端接触力;
图5机械臂2末端接触力矩;
图6机械臂1关节跟踪误差;
图7机械臂2关节跟踪误差。
具体实施方式
下面详细描述本发明的实施例,所述实施例示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本实施例中给出本发明提出的空间双臂机器人约束柔顺稳定的控制方法的原理和稳定性证明,包括推导空间双臂机器人的动力学模型;设计内外双环结构对期望运动进行柔顺化处理;设计性能约束跟踪控制器,实现空间非合作目标的约束柔顺稳定。
建立空间双臂机器人动力学模型:
翻滚目标模型可以简化为刚体,在惯性参数确定已知的情况下,其动力学可由六自由度动力学方程表示。坐标系oxtytzt为定义在目标上的本体坐标系,其中,设o为目标质心,各坐标轴与通过质心的惯量主轴重合。定义刚体绕各轴的主惯量分别是Ιxyz。假设初始时刻本体坐标系和惯性坐标系重合,随后刚体以角速度ωo旋转,其在本体坐标系三个轴向的分量分别记为ωoxoyoz。控制力矩记为uo,目标的速度动力学方程可以表示为:
Figure BDA0002382507130000041
Figure BDA0002382507130000042
其中×为反对称矩阵符号,对于向量x=[x1x2x3],有以下的运算表示:
Figure BDA0002382507130000043
空间机器人的惯性参数准确已知,可以建立其多刚体的动力学模型。坐标系oxeyeze为定义在目标与机械臂末端抓捕点上的坐标系。坐标系oxbybzb为定义在空间机器人基座本体系上的坐标系。其中Hc为空间机器人的惯性矩阵,cc为空间机器人系统所受的科式力,机械臂以及基座的控制量为u,即系统控制输入。假设在稳定的过程中,末端的接触力可以通过末端的力传感器准确测量出,记作fe,Je为机械臂关节以及基座到机械臂末端的雅阁比矩阵,
Figure BDA0002382507130000044
为基座位置与机械臂关节角度组成的广义状态变量,维度为基座的6自由度加双机械臂的自由度2n,一共6+2n,n为机械臂关节个数;
则空间机器人的动力学方程可以写成以下的形式:
Figure BDA0002382507130000051
双层柔顺结构设计:
抓捕后阶段,空间双臂机器人末端测量得到的接触力可以分为两种:有效力与内力。有效力对目标的运动做出贡献,而内力则不影响目标的运动。根据抓捕矩阵,有效力FE与内力FI以及测量的接触力fe满足以下关系等式:
fe=FE+FI(4)
Figure BDA0002382507130000052
Figure BDA0002382507130000053
其中Jo为目标质心到抓捕点的雅各比矩阵。为了实现对目标稳定操作的柔顺,根据末端作用的有效力,可以设计对应的外环柔顺等式如下:
Figure BDA0002382507130000054
其中MO,DO,KO为柔顺等式参数,eO为质心期望运动与柔顺运动之差。利用外环柔顺等式,对期望运动进行调整可得到柔顺的目标稳定运动。根据抓捕点处的运动学以及动力学约束,可以根据内力对双臂末端的运动分别建立柔顺等式,对末端的期望运动进行调整,以实现机械臂之间挤压的柔顺。
Figure BDA0002382507130000055
其中MIj,DIj和KIj为柔顺等式参数,eIi为末端抓捕点处期望运动与柔顺运动的误差。根据内力建立的柔顺等式组成了内环柔顺结构。通过对内环的柔顺等式进行积分,可以得到各个机械臂末端的期望柔顺运动。在此基础上,以基座运动固定为前提,求解空间双臂的逆运动学得到期望关节运动。
约束跟踪控制器
为了提高柔顺控制的可靠性,本研究设计了基于障碍李雅普诺夫函数的约束跟踪控制器。空间双臂机器人跟踪期望轨迹的误差记为es=x-xr。根据中间控制量α,将状态量与中间控制量的差记为
Figure BDA0002382507130000061
选取障碍李雅普诺夫函数为:
Figure BDA0002382507130000062
其中
Figure BDA0002382507130000063
kcli=kci+xri,kcui=kci-xri
选取性能约束变量为kc。对障碍李雅普诺夫函数求导,可得:
Figure BDA0002382507130000064
利用退步法,可以设计抓捕后组合体系统的李雅普诺夫函数如下:
Figure BDA0002382507130000065
对李雅普诺夫函数进行求导,可得:
Figure BDA0002382507130000066
设计控制律如下:
Figure BDA0002382507130000067
Figure BDA0002382507130000068
Figure BDA0002382507130000069
其中xri为列向量xr的元素,xr为经过柔顺处理后的基座位置与机械臂关节角度的期望状态,维度为6+2n;esi为列向量es的元素,es=x-xr为跟踪误差;kci为列向量kc的元素,kc为设定的控制器性能约束参数,维度为6+2n,λ为设定的极小正常数,kei为列向量ke的元素,ke为设定的控制参数,维度为6+2n;KZ=diag(kz)为对角阵,kz也为设定的控制参数,维度为6+2n。
将控制律代入李雅普诺夫函数,可以得到李雅普诺夫函数的导数为负,本发明提出的控制律稳定性得以证明。
本实施例以带双机械臂的空间机器人为例,验证提出的空间机器人基座无扰控制方法,系统的运动学/动力学参数如表1所示。
表1双机械臂空间机器人的运动学/动力学参数
Figure BDA0002382507130000071
表1为实例中使用的空间机器人系统的运动学和动力学参数,双臂的动力学参数相同。
本实施例中,抓捕后阶段目标的初始翻滚速度为[-0.1745 0 0]Trad/s,抓捕点为[0 0.5 0]Tm和[0 -0.5 0]T。我们对惯性参数进行10%的拉偏,从而给出本文中的估计惯性参数,我们基于估计的惯性参数设计了消旋的轨迹,末端速度与角速度收敛至零,基座的位置与姿态保持不变。阻抗控制参数选择为:
MO=200E6,DO=640E6,KO=5E6
MIj=2000E6,DIj=1800diag(1,2,E4),KIj=10E6,j=1,2
控制参数选择为kei=1,kzi=3,kci=1.5,i=1,2,…,20。
根据利用约束跟踪控制器跟踪双层柔顺结构调整后的柔顺期望稳定轨迹,末端接触力与力矩如图2-5所示,关节的跟踪误差如图6-7所示。由图中的结果可以看出本方法成功地跟踪上了期望运动,同时保证了快速地调整了末端的接触力与力矩。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (2)

1.一种空间双臂机器人约束柔顺稳定控制方法,其特征在于:包括以下步骤:
步骤1:建立空间双臂机器人动力学模型:
Figure FDA0002382507120000011
其中Hc为空间机器人的惯性矩阵,cc为空间机器人所受到的科式力,
Figure FDA0002382507120000012
为基座位置与机械臂关节角度组成的广义状态变量,维度为基座的6自由度加双机械臂的自由度2n,一共6+2n,n为机械臂关节个数,
Figure FDA0002382507120000013
为基座加速度与机械臂关节角加速度组成的广义变量,u为系统控制输入,维度为6+2n,Je为机械臂关节以及基座到机械臂末端的雅阁比矩阵,fe为机械臂末端测量的接触力;
步骤2:根据控制律
Figure FDA0002382507120000014
计算系统控制输入u并进行空间双臂机器人控制;其中α为中间控制量,是维度为6+2n的列向量,单个元素的表达式为:
Figure FDA0002382507120000015
xri为列向量xr的元素,xr为经过柔顺处理后的基座位置与机械臂关节角度的期望状态,维度为6+2n;esi为列向量es的元素,es=x-xr为跟踪误差;kcli=kci+xri,kcui=kci-xri,kci为列向量kc的元素,kc为设定的控制器性能约束参数,维度为6+2n,λ为设定的极小正常数,kei为列向量ke的元素,ke为设定的控制参数,维度为6+2n;KZ=diag(kz)为对角阵,kz也为设定的控制参数,维度为6+2n,
Figure FDA0002382507120000016
g为控制补偿量:
Figure FDA0002382507120000017
Figure FDA0002382507120000018
2.根据权利要求1所述一种空间双臂机器人约束柔顺稳定控制方法,其特征在于:统控制输入u为由基座控制力及力矩以及机械臂关节控制力矩组成的列向量。
CN202010087291.8A 2020-02-11 2020-02-11 空间双臂机器人约束柔顺稳定控制方法 Active CN111268182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010087291.8A CN111268182B (zh) 2020-02-11 2020-02-11 空间双臂机器人约束柔顺稳定控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010087291.8A CN111268182B (zh) 2020-02-11 2020-02-11 空间双臂机器人约束柔顺稳定控制方法

Publications (2)

Publication Number Publication Date
CN111268182A true CN111268182A (zh) 2020-06-12
CN111268182B CN111268182B (zh) 2022-08-23

Family

ID=70995007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010087291.8A Active CN111268182B (zh) 2020-02-11 2020-02-11 空间双臂机器人约束柔顺稳定控制方法

Country Status (1)

Country Link
CN (1) CN111268182B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112078836A (zh) * 2020-09-02 2020-12-15 西北工业大学 一种大型充气式柔性消旋装置及方法
CN112372637A (zh) * 2020-10-27 2021-02-19 东方红卫星移动通信有限公司 低轨卫星空间机械臂自适应阻抗柔顺控制方法、模块及系统
CN113664813A (zh) * 2021-08-25 2021-11-19 深圳市优必选科技股份有限公司 双臂机器人控制方法、装置、双臂机器人和可读存储介质
CN115625711A (zh) * 2022-11-08 2023-01-20 沈阳工业大学 一种考虑末端力的双臂机器人协同控制方法
CN117001675A (zh) * 2023-09-28 2023-11-07 江苏云幕智造科技有限公司 一种双臂协作操控非合作目标避障轨迹规划方法
CN117182929A (zh) * 2023-11-06 2023-12-08 北京控制工程研究所 一种双臂机器人在轨装配的柔顺控制方法及装置
CN115625711B (zh) * 2022-11-08 2024-04-30 沈阳工业大学 一种考虑末端力的双臂机器人协同控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101733746A (zh) * 2009-12-22 2010-06-16 哈尔滨工业大学 空间机器人非合作目标自主识别与捕获方法
CN108927803A (zh) * 2018-07-25 2018-12-04 西北工业大学 一种在连续碰撞条件下空间机器人目标抓捕稳定方法
CN109606753A (zh) * 2018-11-11 2019-04-12 上海宇航系统工程研究所 一种空间双臂机器人协同抓捕目标的控制方法
CN110125936A (zh) * 2019-05-15 2019-08-16 清华大学深圳研究生院 一种空间机器人的柔顺控制方法和地面实验验证系统
CN111290269A (zh) * 2020-02-11 2020-06-16 西北工业大学深圳研究院 一种空间机器人的自适应柔顺稳定控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101733746A (zh) * 2009-12-22 2010-06-16 哈尔滨工业大学 空间机器人非合作目标自主识别与捕获方法
CN108927803A (zh) * 2018-07-25 2018-12-04 西北工业大学 一种在连续碰撞条件下空间机器人目标抓捕稳定方法
CN109606753A (zh) * 2018-11-11 2019-04-12 上海宇航系统工程研究所 一种空间双臂机器人协同抓捕目标的控制方法
CN110125936A (zh) * 2019-05-15 2019-08-16 清华大学深圳研究生院 一种空间机器人的柔顺控制方法和地面实验验证系统
CN111290269A (zh) * 2020-02-11 2020-06-16 西北工业大学深圳研究院 一种空间机器人的自适应柔顺稳定控制方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112078836A (zh) * 2020-09-02 2020-12-15 西北工业大学 一种大型充气式柔性消旋装置及方法
CN112078836B (zh) * 2020-09-02 2022-10-28 西北工业大学 一种大型充气式柔性消旋装置及方法
CN112372637A (zh) * 2020-10-27 2021-02-19 东方红卫星移动通信有限公司 低轨卫星空间机械臂自适应阻抗柔顺控制方法、模块及系统
CN112372637B (zh) * 2020-10-27 2022-05-06 东方红卫星移动通信有限公司 低轨卫星空间机械臂自适应阻抗柔顺控制方法、模块及系统
CN113664813A (zh) * 2021-08-25 2021-11-19 深圳市优必选科技股份有限公司 双臂机器人控制方法、装置、双臂机器人和可读存储介质
CN113664813B (zh) * 2021-08-25 2022-07-29 深圳市优必选科技股份有限公司 双臂机器人控制方法、装置、双臂机器人和可读存储介质
CN115625711A (zh) * 2022-11-08 2023-01-20 沈阳工业大学 一种考虑末端力的双臂机器人协同控制方法
CN115625711B (zh) * 2022-11-08 2024-04-30 沈阳工业大学 一种考虑末端力的双臂机器人协同控制方法
CN117001675A (zh) * 2023-09-28 2023-11-07 江苏云幕智造科技有限公司 一种双臂协作操控非合作目标避障轨迹规划方法
CN117182929A (zh) * 2023-11-06 2023-12-08 北京控制工程研究所 一种双臂机器人在轨装配的柔顺控制方法及装置
CN117182929B (zh) * 2023-11-06 2024-01-23 北京控制工程研究所 一种双臂机器人在轨装配的柔顺控制方法及装置

Also Published As

Publication number Publication date
CN111268182B (zh) 2022-08-23

Similar Documents

Publication Publication Date Title
CN111268182B (zh) 空间双臂机器人约束柔顺稳定控制方法
CN110104216B (zh) 一种针对运动学冗余双臂空间机器人的协同路径规划方法
CN111290269B (zh) 一种空间机器人的自适应柔顺稳定控制方法
CN110202574B (zh) 基于环境刚度估计的机器人自适应混合阻抗/导纳控制方法
Erhart et al. An impedance-based control architecture for multi-robot cooperative dual-arm mobile manipulation
EP1728600B1 (en) Controlling the trajectory of an effector
Nguyen et al. Aerial tool operation system using quadrotors as rotating thrust generators
CN109606753B (zh) 一种空间双臂机器人协同抓捕目标的控制方法
CN110340898B (zh) 一种自由漂浮空间机械臂自适应容错控制方法
CN107505947B (zh) 一种空间机器人捕获目标后消旋及协调控制方法
CN109426147B (zh) 捕获卫星后组合航天器的自适应增益调整控制方法
CN109591018B (zh) 一种用于目标捕获后稳定的自由漂浮空间机械臂控制方法
CN109352656B (zh) 一种具有时变输出约束的多关节机械臂控制方法
CN105773617A (zh) 空间机器人的三指式抓持器碰撞预测方法
CN112809666B (zh) 一种基于神经网络的5-dof机械臂力位跟踪算法
Chwa et al. Online trajectory planning of robot arms for interception of fast maneuvering object under torque and velocity constraints
CN111975770B (zh) 一种采用空间双臂机器人的自旋目标抓捕方法
CN114851193B (zh) 一种空间机械臂与未知环境接触过程的智能柔顺操控方法
Chunting et al. Adaptive hybrid impedance control for a dual-arm robot manipulating an unknown object
Jin et al. A motion planning method based vision servo for free-flying space robot capturing a tumbling satellite
Yang et al. An adaptive force control method for 7-dof space manipulator repairing malfunctioning satellite
Nagendran et al. Dynamic capture of free-moving objects
Doulgeri et al. Force/position tracking of a robot in compliant contact with unknown stiffness and surface kinematics
CN115625711B (zh) 一种考虑末端力的双臂机器人协同控制方法
Bai et al. Dynamics model and adaptive backstepping control on the 7-DOF manipulators of a humanoid robot

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