CN108958169B - 一种基于s曲线的pvt控制方法 - Google Patents

一种基于s曲线的pvt控制方法 Download PDF

Info

Publication number
CN108958169B
CN108958169B CN201810842983.1A CN201810842983A CN108958169B CN 108958169 B CN108958169 B CN 108958169B CN 201810842983 A CN201810842983 A CN 201810842983A CN 108958169 B CN108958169 B CN 108958169B
Authority
CN
China
Prior art keywords
curve
time
acceleration
period
movement
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
CN201810842983.1A
Other languages
English (en)
Other versions
CN108958169A (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN201810842983.1A priority Critical patent/CN108958169B/zh
Publication of CN108958169A publication Critical patent/CN108958169A/zh
Application granted granted Critical
Publication of CN108958169B publication Critical patent/CN108958169B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/41Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Feedback Control In General (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种基于S曲线的PVT控制方法,它解决了现有技术中加速度变化过大在运行过程中易发生抖动的问题,具有能够保证速度、加速度连续,速度在时间内变化较大时加速度不会发生过大抖动现象,使得数控***加工时更平稳,加工效率更高的效果;其技术方案为:通过给定的始末位置、始末速度、始末运动时间以及运动参数限制,利用S曲线加减速得到运动总时间T’;判断运动总时间T’与给定始末运动时间T的关系,在T’不满足给定始末运动时间T的情况下,将T’按照等周期S曲线圆整至时间T;当T’满足给定始末运动时间T时,计算插补周期下的插补位移。

Description

一种基于S曲线的PVT控制方法
技术领域
本发明涉及数控***运动控制领域,尤其涉及一种基于S曲线的PVT控制方法。
背景技术
PVT(位置-速度-时间)控制方式是输入首末位置P、速度V、时间T,同时设置运动参数限制,控制***自动规划出周期插补位移,在指定的时间T内,满足首末速度V的前提下达到指定位置P。
现有技术中提出了使用3次多项式拟合PVT曲线,已知首末位置和首末速度,构造四元一次方程组,求解得到以时间t做自变量的一元三次方程。在周期插补时,控制***输入插补时间求得插补位移,驱动执行器至指定位置。这种方式计算简单,容易规划,但是其加速度不连续,在始末位置处会有较大冲击。
现有技术中提出了使用5次多项式拟合PVT曲线,已知首末位置、始末速度,始末加速度均等于0,构造六元一次方程组求解一元六次方程。这种方式计算较复杂,其始末加速度均约束为0,不会发生冲击,但是当速度V在时间T内变化较大时加速度变化过大发生抖动现象。
综上所述,针对现有技术如何保证数控***PVT运动控制中加速度连续且不会发生抖动现象,尚缺乏有效的解决方案。
发明内容
为了克服现有技术的不足,本发明提供了一种基于S曲线的PVT控制方法,其具有能够保证速度、加速度连续,速度在时间内变化较大时加速度不会发生过大抖动现象,使得数控***加工时更平稳,加工效率更高的效果。
本发明采用下述技术方案:
一种基于S曲线的PVT控制方法,通过给定的始末位置、始末速度、始末运动时间以及运动参数限制,利用S曲线加减速得到运动总时间T’;判断运动总时间T’与给定始末运动时间T的关系,在T’不满足给定始末运动时间T的情况下,将T’按照等周期S曲线圆整至时间T;当T’满足给定始末运动时间T时,计算插补周期下的插补位移。
进一步的,构造以时间t为自变量的等周期S曲线函数:t=Φ(t');
其中,t’表示圆整时间。
进一步的,当T'>T时,构造等周期S曲线缩短时间T'至T;圆整之后运动时间变为T,比例函数为Φ'(t)。
进一步的,求解Φ'max和Φ”'max,其中,
Figure BDA0001745984870000021
根据Φ'max和Φ”max得到Φ(t)和Φ'(t)的函数方程。
进一步的,由S曲线规划出运行时间区间为[0,T'],经等周期S曲线圆整后的时间区间为[0,T],得到运动总位移:
Figure BDA0001745984870000022
进一步的,将插补周期代入运动总位移公式中得到插补位移。
进一步的,当T'<T时,构造等周期S曲线延长时间T'至T。
进一步的,求解运动总时间T'时,给定的始末位置P、始末速度V、始末加速度均为0。
与现有技术相比,本发明的有益效果是:
本发明基于等周期S曲线进行PVT控制,得到运动的总时间T',在不满足给定的始末运动时间T的情况下,将T'按照等周期S曲线变化圆整至T,保证在规定的时间、始末速度下,到达指定位移;在整个运行过程中加速度连续,当速度变化较大时,加速度不会发生过大的抖动现象,保证数控***的稳定性;使得数控***加工时更平稳,加工效率更高。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为本发明的流程图;
图2为S曲线加减速示意图;
图3为本发明的等周期S曲线缩短时间示意图;
图4为本发明的等周期S曲线图;
图5为本发明的等周期S曲线延长时间示意图;
图6为本发明的PVT大变速速度曲线图;
图7为本发明的PVT大变速加速度曲线图;
图8为本发明的PVT小变速速度曲线图;
图9为本发明的PVT小变速加速度曲线图;
图10为五次多项式PVT加速度曲线图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
正如背景技术所介绍的,现有技术中存在加速度变化过大在运行过程中易发生抖动的不足,为了解决如上的技术问题,本申请提出了一种基于S曲线的PVT控制方法。
本申请的一种典型的实施方式中,提供了一种基于S曲线的PVT控制方法,采用等周期S曲线圆整运行时间T,通过S曲线加减速度求得该位移下的总时间T',当T'=T时,计算插补周期下的插补位移,并进行插补输出;当T'>T或者T'<T时,将T'按照等周期S曲线圆整至时间T,保证在规定的时间、始末速度下,到达指定位移。
在整个的运行过程中加速度连续,当速度变化较大时,加速度不会发生过大的抖动现象,保证数控***的稳定性。
S曲线加减速全过程中的速度、加速度和加加速度曲线如图2所示,在PVT控制中,始末点速度V不相同,在整个过程中可分为三大段:加速段、匀速段和减速段,如果距离很短时可能不存在匀速段情况。
加速段又分为加加速段①、匀加速段②和减加速段③,减速段分为加减速段⑤、匀减速段⑥和减减速段⑦。
由于始末速度不同,则加减速的时间也不同;由于S型曲线对称,加加速和减加速的时间相同t1=t3,减加速和减减速的时间也相同t5=t7。在实际的运动过程中,S曲线可能不同时包含匀速段、匀加速段、匀减速段。
PVT控制中,通过给定的始末位置P、始末速度V、始末加速度均为0,以及运动参数限制,求解运动总时间T',比较和给定时间T的关系,进而构造以插补时间t'为自变量的等周期S曲线函数:
t=Φ(t') (1)
其中,t'表示圆整时间。
本申请基于S曲线的PVT控制方法具体如下:
(1)当T'>T时,构造等周期S曲线缩短时间T'至T。
如图3所示,横轴为时间t,纵轴为比例函数Φ'(t),未圆整之前运动时间为T',比例函数横为1;圆整之后,运动时间变为T,比例函数变为Φ'(t),阴影部分的面积与凸起的面积相等,圆整前后的面积不变。
①计算等周期S曲线方程:
S曲线为分段函数,求解S曲线方程先求解其中的参数:Φ'max、Φ”max,如图4所示,等周期S曲线中比例函数的最大值在
Figure BDA0001745984870000041
处,所以
Figure BDA0001745984870000042
等周期S曲线图如图4所示,S1=S2=S3=S4,S曲线的积分面积与虚线三角形的面积相等。
已知t1=t2=t3=t4=t5=t6,由等积法得:
Figure BDA0001745984870000043
Figure BDA0001745984870000051
Figure BDA0001745984870000052
通过公式(4)、(5)及S曲线特征方程,得到Φ(t)和Φ'(t)的函数方程。
②、PVT位移插补
由S曲线规划出运行时间区间为[0,T'],经等周期S曲线圆整后的时间区间为[0,T],插补周期发生变化,运动总位移与速度的关系表示为:
Figure BDA0001745984870000053
将插补周期t'代入Φ(t)和Φ'(t)的函数方程,得到当前插补周期的时间Φ(t')和比例函数Φ'(t'),通过公式(6)求得插补位移。
(2)、当T'<T时,构造等周期S曲线延长时间T'至T。
如图5所示,阴影部分的面积和凹下的面积相等,圆整前后的面积不变。
①、计算等周期S曲线方程:
Figure BDA0001745984870000054
时,等周期S曲线方程的求解方式与(1)中相同,此处不再赘述。
②、PVT位移插补
由S曲线规划出运行时间区间为[0,T'],经等周期S曲线延长后的时间区间为[0,T],代入插补周期由公式(6)求解插补位移。
对于公式(6),已知位移对时间t是连续可导的,对位移求一阶导数可得到速度:
v(t)=S'(t)=v(Φ(t))Φ'(t)
Φ(t)为等周期S曲线方程,在定义区间[0,T]上连续可导,因此,v(t)在区间内连续,对速度求一阶导数可得加速度:
a(t)=S”(t)=v'(t)=a(Φ(t))Φ'2(t)+v(Φ(t))Φ”(t)
由S曲线特性知Φ'(t)、Φ”(t)在定义区间上连续,a(t)在区间内也连续。
因此,本申请的方法是可行的。
本申请的另一种实施方式中,以数控机床单轴做PVT控制,其运动参数限制为:最大速度限制400°/s,最大加速度限制1000°/s2,最大加加速度限制80000°/s3
分别给定两组PVT,表1为相邻周期T速度变化较大,表2为相邻周期T速度变化较小。
图10为同条件下的五次多项式加速度曲线图,由图中可以看出加速度变化波动很大。
图6~图7为表1的周期速度、加速度曲线图,从图中可以看出,速度和加速度在时间区间内连续,且当速度变化较大时,加速度在±17000内浮动。
图8~图9为表2的周期速度、加速度曲线图,从图中可以看出,速度和加速度在时间区间内平稳连续,无波动,将图7~图9与图10对比,完整地证明了本申请的可行性。
表1PVT参数表
序号 位置 速度 周期
1 1258.739899 63.973827 21.95
2 1261.940364 64.009293 22.00
3 1265.142598 54.044679 22.05
4 1268.346597 44.079985 22.10
5 1271.552357 34.115211 22.15
6 1274.759875 24.150359 22.20
7 1277.969147 14.185428 22.25
8 1281.180168 4.220419 22.30
9 1284.392934 0.0 22.35
表2PVT参数表
序号 位置 速度 周期
1 180.057763 0.0 0.0
2 180.168636 2.217459 0.05
3 180.347390 3.575093 0.10
4 180.589727 4.846733 0.15
5 180.891812 6.041703 0.20
6 181.250210 7.167962 0.25
7 181.661828 8.232346 0.30
8 182.123866 9.240765 0.35
9 182.633784 10.198358 0.40
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (2)

1.一种基于S曲线的PVT控制方法,其特征在于,通过给定的始末位置、始末速度、始末运动时间以及运动参数限制,利用S曲线加减速得到运动总时间T’;判断运动总时间T’与给定始末运动时间T的关系,在T’不满足给定始末运动时间T的情况下,将T’按照等周期S曲线圆整至时间T;当T’满足给定始末运动时间T时,计算插补周期下的插补位移;
构造以时间t为自变量的等周期S曲线函数:t=Φ(t');
其中,t’表示圆整时间;
当T'>T时,构造等周期S曲线缩短时间T'至T;圆整之后运动时间变为T,比例函数为Φ'(t);
S曲线为分段函数,求解S曲线方程先求解其中的参数:Φ'max、Φ”max,等周期S曲线中比例函数的最大值在
Figure FDA0002753664640000011
处,所以
Figure FDA0002753664640000012
等周期S曲线图中,S1=S2=S3=S4,S曲线的积分面积与虚线三角形的面积相等;
已知t1=t2=t3=t4=t5=t6,由等积法得:
Figure FDA0002753664640000013
Figure FDA0002753664640000014
Figure FDA0002753664640000015
通过公式(4)、(5)及S曲线特征方程,得到Φ(t)和Φ'(t)的函数方程;
由S曲线规划出运行时间区间为[0,T'],经等周期S曲线圆整后的时间区间为[0,T],得到运动总位移:
Figure FDA0002753664640000016
将插补周期代入运动总位移公式中得到插补位移;
当T'<T时,构造等周期S曲线延长时间T'至T。
2.根据权利要求1所述的一种基于S曲线的PVT控制方法,其特征在于,求解运动总时间T'时,给定的始末位置P、始末速度V、始末加速度均为0。
CN201810842983.1A 2018-07-27 2018-07-27 一种基于s曲线的pvt控制方法 Active CN108958169B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810842983.1A CN108958169B (zh) 2018-07-27 2018-07-27 一种基于s曲线的pvt控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810842983.1A CN108958169B (zh) 2018-07-27 2018-07-27 一种基于s曲线的pvt控制方法

Publications (2)

Publication Number Publication Date
CN108958169A CN108958169A (zh) 2018-12-07
CN108958169B true CN108958169B (zh) 2021-04-06

Family

ID=64466007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810842983.1A Active CN108958169B (zh) 2018-07-27 2018-07-27 一种基于s曲线的pvt控制方法

Country Status (1)

Country Link
CN (1) CN108958169B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112532146B (zh) * 2020-11-16 2023-02-28 深圳市科曼医疗设备有限公司 加减速运动控制方法、装置、设备和介质
CN113253678B (zh) * 2021-04-29 2022-07-05 中国科学院微电子研究所 精密运动平台的运动轨迹规划方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0764622A (ja) * 1993-08-26 1995-03-10 Matsushita Electric Ind Co Ltd ロボットの軌道補間装置
EP1517205A2 (en) * 2003-08-05 2005-03-23 Fanuc Ltd Controller for machine
CN1971457A (zh) * 2005-11-25 2007-05-30 中国科学院沈阳计算技术研究所有限公司 一种用于数控机床的加减速控制方法
CN101493687A (zh) * 2009-03-02 2009-07-29 广西大学 实时前瞻全程加减速控制的nurbs曲线自适应分段插补方法
CN102298359A (zh) * 2011-06-18 2011-12-28 山东大学 一种数控***中加减速过渡平滑控制方法
CN102419570A (zh) * 2011-09-29 2012-04-18 上海大学 数控机床高速加工的加减速前瞻控制方法
CN102540978A (zh) * 2010-12-09 2012-07-04 中国科学院沈阳计算技术研究所有限公司 面向高速加工的表面质量优先的样条曲线实时插补方法
CN102799146A (zh) * 2012-08-08 2012-11-28 成都乐创自动化技术股份有限公司 用于数控设备控制***的s形加减速控制的速度规划方法
CN103076759A (zh) * 2012-12-31 2013-05-01 北京配天大富精密机械有限公司 一种获取数控***的s型曲线的方法及其装置
CN103801973A (zh) * 2012-11-09 2014-05-21 沈阳高精数控技术有限公司 基于伺服刀架控制器的伺服电机的优化控制方法
CN103970073A (zh) * 2013-01-24 2014-08-06 北京配天大富精密机械有限公司 一种用于数控***的加减速规划方法、装置及数控机床
CN104135212A (zh) * 2014-07-21 2014-11-05 华南理工大学 一种柔性化运动控制ip核及实现方法
CN106168790A (zh) * 2016-02-29 2016-11-30 华南理工大学 一种在线改变目标速度和位置的s形加减速控制方法
CN107765639A (zh) * 2017-09-30 2018-03-06 山东大学 一种s曲线加减速的圆整误差补偿方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004059966B3 (de) * 2004-12-13 2006-06-22 Siemens Ag Verfahren und Einrichtung zur Bewegungsführung eines bewegbaren Maschinenelements einer numerisch gesteurten Maschine

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0764622A (ja) * 1993-08-26 1995-03-10 Matsushita Electric Ind Co Ltd ロボットの軌道補間装置
EP1517205A2 (en) * 2003-08-05 2005-03-23 Fanuc Ltd Controller for machine
CN1971457A (zh) * 2005-11-25 2007-05-30 中国科学院沈阳计算技术研究所有限公司 一种用于数控机床的加减速控制方法
CN101493687A (zh) * 2009-03-02 2009-07-29 广西大学 实时前瞻全程加减速控制的nurbs曲线自适应分段插补方法
CN102540978A (zh) * 2010-12-09 2012-07-04 中国科学院沈阳计算技术研究所有限公司 面向高速加工的表面质量优先的样条曲线实时插补方法
CN102298359A (zh) * 2011-06-18 2011-12-28 山东大学 一种数控***中加减速过渡平滑控制方法
CN102419570A (zh) * 2011-09-29 2012-04-18 上海大学 数控机床高速加工的加减速前瞻控制方法
CN102799146A (zh) * 2012-08-08 2012-11-28 成都乐创自动化技术股份有限公司 用于数控设备控制***的s形加减速控制的速度规划方法
CN103801973A (zh) * 2012-11-09 2014-05-21 沈阳高精数控技术有限公司 基于伺服刀架控制器的伺服电机的优化控制方法
CN103076759A (zh) * 2012-12-31 2013-05-01 北京配天大富精密机械有限公司 一种获取数控***的s型曲线的方法及其装置
CN103970073A (zh) * 2013-01-24 2014-08-06 北京配天大富精密机械有限公司 一种用于数控***的加减速规划方法、装置及数控机床
CN104135212A (zh) * 2014-07-21 2014-11-05 华南理工大学 一种柔性化运动控制ip核及实现方法
CN106168790A (zh) * 2016-02-29 2016-11-30 华南理工大学 一种在线改变目标速度和位置的s形加减速控制方法
CN107765639A (zh) * 2017-09-30 2018-03-06 山东大学 一种s曲线加减速的圆整误差补偿方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
始末速度不为零的S型加减速时间规划算法研究;杨亮亮 等;《机械工程学报》;20161230;第52卷(第23期);第199-206页 *

Also Published As

Publication number Publication date
CN108958169A (zh) 2018-12-07

Similar Documents

Publication Publication Date Title
CN108549328B (zh) 自适应速度规划方法及***
CN108958169B (zh) 一种基于s曲线的pvt控制方法
US10768602B2 (en) Curve velocity planning method, device, and numerical control machining route data processing method thereof
CN107368639B (zh) 速度规划方法、装置、计算机设备和存储介质
CN107825424B (zh) 一种高速机械手减少残余振动的非对称s型轨迹规划方法
CN102419570B (zh) 数控机床高速加工的加减速前瞻控制方法
WO2019119268A1 (zh) 数控***、具有存储功能的装置及nurbs曲线插补方法
CN111897216A (zh) 一种多运动段速度规划和插补方法
EP2022608A2 (en) Improved blending algorithm for trajectory planning
CN107850887A (zh) S型曲线规划方法、装置及数控机床
CN112659126A (zh) 一种基于非均匀b样条曲线的工业机器人轨迹段过渡方法
CN114077226B (zh) S型曲线速度规划方法、控制终端及计算机可读存储介质
CN110874082B (zh) 工件拐角的切割路径确定方法、设备和存储介质
CN102298359B (zh) 一种数控***中加减速过渡平滑控制方法
JP6386511B2 (ja) 工具経路生成装置、工具経路生成方法及び工具経路生成プログラム
CN104181860A (zh) 数控机床s型加减速控制方法
CN111310106B (zh) 一种基于逐步逼近原始轮廓的切割轮廓拟合优化方法
CN101976308A (zh) 一种非圆曲线工艺设计及数控加工方法
CN110185523A (zh) 尿素喷射量控制方法及装置
CN105278462A (zh) 机器人控制***轨迹加减速插补算法
CN112327758B (zh) 基于b样条拟合及分段插补的小线段刀轨局部光顺方法
CN111381565B (zh) 三次多项式s曲线加减速的加加速度递归计算限制方法
CN106444636B (zh) 一种适用于连续混合曲线的速度规划参数选择方法
CN110703684A (zh) 一种端点速度不限的轨迹规划方法及装置
CN110175339B (zh) 一种s型加减速的最优时间分配方法及装置

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