CN109508010B - 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法 - Google Patents

一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法 Download PDF

Info

Publication number
CN109508010B
CN109508010B CN201811546402.6A CN201811546402A CN109508010B CN 109508010 B CN109508010 B CN 109508010B CN 201811546402 A CN201811546402 A CN 201811546402A CN 109508010 B CN109508010 B CN 109508010B
Authority
CN
China
Prior art keywords
grid
mobile
robot
map
grid points
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
CN201811546402.6A
Other languages
English (en)
Other versions
CN109508010A (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.)
Yancheng Institute of Technology
Original Assignee
Yancheng Institute 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 Yancheng Institute of Technology filed Critical Yancheng Institute of Technology
Priority to CN201811546402.6A priority Critical patent/CN109508010B/zh
Publication of CN109508010A publication Critical patent/CN109508010A/zh
Application granted granted Critical
Publication of CN109508010B publication Critical patent/CN109508010B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0287Control of position or course in two dimensions specially adapted to land vehicles involving a plurality of land vehicles, e.g. fleet or convoy travelling
    • G05D1/0289Control of position or course in two dimensions specially adapted to land vehicles involving a plurality of land vehicles, e.g. fleet or convoy travelling with means for avoiding collisions between vehicles
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0221Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving a learning process

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提出一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,该方法在分配栅格点时,每个栅格点至多只分配给一台移动机器人,能够有效避免移动机器人间发生碰撞。每台机器人每轮最多只能获得预先设置的NC个栅格点,避免造成系统中大量机器人的长时间等待现象,实现了多机器人共享栅格地图的目的,有效提高了多机器人系统整体效率。在分配过程中采用前瞻一步的策略,实时构建系统前瞻一步后的阻塞图模型,应用相应的环路死锁检测与避免方法,有效避免了环路死锁现象。

Description

一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁 动态分配方法
技术领域
本发明涉及移动机器人领域,尤其是一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法。该方法主要应用于基于栅格地图的多移动机器人系统,主要解决当多移动机器人共用同一栅格地图时的路径冲突,多移动机器人间发生碰撞及系统死锁的技术问题,提高多移动机器人系统的整体效率和智能化水平。
背景技术
当多移动机器人共享同一栅格地图时,移动机器人间必然存在路径冲突,为了避免各移动机器人间发生碰撞,目前普遍采用的方法是在机器人路径规划时避开其他机器人已规划路径,机器人在运行到终点前,将一直独占整条路径栅格点,当系统中移动机器人密度较高时易造成大量机器人的长时间等待现象,降低了系统效率。
发明内容
发明目的:为解决上述技术问题,本发明提出一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,该方法主要应用于基于栅格地图的多移动机器人系统,主要解决当多移动机器人共用同一栅格地图时的路径冲突,避免多移动机器人间发生碰撞及系统死锁,提高多移动机器人系统的整体效率和智能化水平。
技术方案:为实现上述技术效果,本发明提出的技术方案为:
一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,所述多个移动机器人共用一个栅格地图;所述分配方法包括步骤;
(1)获取栅格地图中栅格点的坐标、栅格点的状态、栅格地图中的障碍物位置、各移动机器人的位置、各移动机器人已规划的路径;根据获取的信息标识所有栅格点及其状态:记所有栅格点的集合为V={vj|j∈N且1≤j≤NV},其中,vj表示第j个栅格点,NV为栅格点总数;vj的状态标识为Occupy(vj),Occupy(vj)=0表示栅格点vj未被占用,Occupy(vj)=1表示栅格点vj已被占用,占用栅格点vj的机器人记为R(vj);
(2)记第i个移动机器人为Ri,i∈N,且1≤i≤NR,N表示正整数,NR表示移动移动机器人的总数;记Ri已规划的路径为Path(Ri),将系统根据Path(Ri)已分配给Ri的栅格点中最靠近Path(Ri)终点的栅格点作为Ri的临时停靠点,记为FinOcc(Ri);记从FinOcc(Ri)开始沿着路径Path(Ri)前向第k个栅格点为
Figure BDA0001908245730000021
(3)逐一为移动机器人分配栅格点,其中,为移动机器人Ri分配栅格点的步骤包括:
(3-1)初始化k=1;
(3-2)根据
Figure BDA0001908245730000022
的状态标识判断栅格点
Figure BDA0001908245730000023
是否已被其他移动机器人占据,若是,则直接进入步骤(3-6);否则,执行步骤(3-3);
(3-3)假定将栅格点
Figure BDA0001908245730000026
分配给机器人Ri,根据假定的分配结果构建阻塞图GQ(VQ,EQ),GQ(VQ,EQ)为有向图,VQ为顶点集,EQ为边集,
VQ={vj|vj∈VT}
Figure BDA0001908245730000024
其中,vj表示第j个栅格点,VT表示所有移动机器人临时停靠点的集合,ejk表示从栅格点vj到vk的单向路径;
(3-4)根据构建好的阻塞图,检测阻塞图中是否存在环路死锁,若不存在,则执行步骤(3-5),否则,执行步骤(3-6);
(3-5)将栅格点
Figure BDA0001908245730000025
分配给机器人Ri;计算k=k+1,判断是否满足k>Nc,若满足,则执行步骤(3-6 ),否则,返回步骤(3-2);其中,Nc为预设的每次为移动机器人分配栅格数的最大值;
(3-6)结束为移动机器人Ri分配栅格点。
(4)当为所有移动机器人分配过栅格点后,等待所有移动机器人按照各自路径顺序通过已分配栅格点,并停留在临时停靠栅格点后,返回步骤(3),重新对所有移动机器人进行栅格点分配。
进一步的,判断阻塞图中是否存在环路死锁的方法为:当阻塞图中存在有向回路时表示存在环路死锁,一个有向回路表示一个发生环路死锁的阻塞路径。
进一步的,检索阻塞图中的环路死锁的方法为:
1)假定移动机器人Ri进入阻塞图GQ(VQ,EQ),初始化:
Figure BDA0001908245730000031
L={vcu},其中,vcu表示阻塞图中当前正在搜索的栅格点,L为发生环路死锁的阻塞路径;
2)判断阻塞图中是否存在vcu的前向邻节点Adj+(vcu);若不存在,则不会发生环路死锁;若存在,则更新:vcu=Adj+(vcu),L=L+{Adj+(vcu)};
3)判断是否满足
Figure BDA0001908245730000032
若满足,则表示发生环路死锁,并输出L;若不满足,则返回步骤2)。
有益效果:与现有技术相比,本发明具有以下优势:
本发明在分配栅格点的过程中采用前瞻一步的策略,实时构建系统前瞻一步后的阻塞图模型,应用相应的环路死锁检测与避免方法,有效避免了环路死锁现象。避免造成系统中大量机器人的长时间等待现象,实现了多机器人共享栅格地图的目的,有效提高了多机器人系统整体效率。
附图说明
图1为本发明的整体流程图;
图2为环路死锁的示意图;
图3为实施例中给出的栅格地图状态图;
图4为图1相应的阻塞图;
图5为环路死锁的阻塞路径搜索流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
图1所示为本发明的整体流程图
基于栅格地图的多移动机器人系统由栅格地图及若干移动机器人组成。NR为系统中的移动机器人数量。栅格地图可用图G(V,E)来描述,其中V={vj|j∈N且1≤j≤NV} 表示栅格点的集合,NV为栅格地图中的栅格点数,N表示正整数,E表示栅格边的集合,E={eij|eij=(vi,vj),其中vi,vj∈V且vi≠vj},NE为栅格边的条数,vi称为边eij的始点,vj称为边eij的终点。移动机器人可沿栅格边的方向双向运行,但只能停靠于栅格点,由于多台移动机器人共享同一栅格地图,因此机器人之间的相互冲突在所难免,为了避免碰撞,必须针对机器人间可能出现的干涉,提出对应的避碰方法。
在基于栅格地图的多移动机器人系统中,容易出现某台移动机器人被另一移动机器人阻塞而暂时不能运动的情况,如果若干移动机器人形成循环阻塞,就造成了如图2所示的环路死锁。
为解决这一问题,本发明提出一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,下面基于图3所示的栅格地图状态图,对本发明的具体步骤进行阐述,本发明的具体流程如图1所示,包括步骤:
(1)获取栅格地图中栅格点的坐标、栅格点的状态、栅格地图中的障碍物位置、各移动机器人的位置、各移动机器人已规划的路径;根据获取的信息标识所有栅格点及其状态:记所有栅格点的集合为V={vj|j∈N且1≤j≤NV},其中,vj表示第j个栅格点,NV为栅格点总数;vj的状态标识为Occupy(vj),Occupy(vj)=0表示栅格点vj未被占用,Occupy(vj)=1表示栅格点vj已被占用,占用栅格点vj的机器人记为R(vj);
(2)记第i个移动机器人为Ri,i∈N,且1≤i≤NR,N表示正整数,NR表示移动移动机器人的总数;记Ri已规划的路径为Path(Ri),将系统根据Path(Ri)已分配给Ri的栅格点中最靠近Path(Ri)终点的栅格点作为Ri的临时停靠点,记为FinOcc(Ri);记从FinOcc(Ri)开始沿着路径Path(Ri)前向第k个栅格点为
Figure BDA0001908245730000041
如:Path(Ri)={v1,v2,v7,v8,v13},其中{v1,v2,v7}已分配给机器人Ri,则FinOcc(Ri)=v7
Figure BDA0001908245730000042
(3)逐一为移动机器人分配栅格点,分配时的分配规则为:
规则1:为了避免移动机器人间发生碰撞,每个栅格点至多只允许分配给一台移动机器人。
规则2:为了避免某机器人一次分配的栅格点过多,造成系统中大量机器人的长时间等待现象,每个机器人每轮至多只允许分配NV个栅格点,提高多机器人系统的整体效率,实现多机器人共享栅格地图的目的。
规则3:为了避免出现环路死锁,在分配过程中采用前瞻一步的策略,实时构建系统阻塞图模型,进而预测并避免环路死锁。
基于上述规则,开始对每个移动机器人分配栅格点,其中为移动机器人Ri分配栅格点的步骤包括:
(3-1)初始化k=1;
(3-2)根据
Figure BDA0001908245730000051
的状态标识判断栅格点
Figure BDA0001908245730000052
是否已被其他移动机器人占据,若是,则直接进入步骤(3-6);否则,执行步骤(3-3);
(3-3)假定将栅格点
Figure BDA0001908245730000053
分配给机器人Ri,根据假定的分配结果构建阻塞图GQ(VQ,EQ),GQ(VQ,EQ)为有向图,VQ为顶点集,EQ为边集,
VQ={vj|vj∈VT}
Figure BDA0001908245730000054
其中,vj表示第j个栅格点,VT表示所有移动机器人临时停靠点的集合,ejk表示从栅格点vj到vk的单向路径;显然,由于每个栅格点至多只分配给一台移动机器人,所以阻塞图中的顶点数等于系统中的移动机器人数,边数小于移动机器人数;图1相应的阻塞图如图4所示。
(3-4)根据构建好的阻塞图,检测阻塞图中是否存在环路死锁,若不存在,则执行步骤(3-5),否则,执行步骤(3-6);
(3-5)将栅格点
Figure BDA0001908245730000055
分配给机器人Ri;计算k=k+1,判断是否满足k>Nc,若满足,则执行步骤(3-6 ),否则,返回步骤(3-2);其中,Nc为预设的每次为移动机器人分配栅格数的最大值;
(3-6)结束为移动机器人Ri分配栅格点。
(4)当为所有移动机器人分配过栅格点后,等待所有移动机器人按照各自路径顺序通过已分配栅格点,并停留在临时停靠栅格点后,返回步骤(3),重新对所有移动机器人进行栅格点分配。
上述过程中,环路死锁的搜索方法如图5所示:
定义环路死锁:阻塞图中的一个有向回路即表示一个环路死锁L,因此L可采用阻塞图的有序顶点集L{vi,vj,vk…vn}表示,对环路中任意两相邻顶点有(ej,ek)∈EQ
以移动机器人Ri为例,环路死锁的搜索流程为:
1)假定移动机器人Ri进入阻塞图GQ(VQ,EQ),初始化:
Figure BDA0001908245730000062
L={vcu},其中,vcu表示阻塞图中当前正在搜索的栅格点,L为发生环路死锁的阻塞路径;
2)判断阻塞图中是否存在vcu的前向邻节点Adj+(vcu);若不存在,则不会发生环路死锁;若存在,则更新:vcu=Adj+(vcu),L=L+{Adj+(vcu)};
3)判断是否满足
Figure BDA0001908245730000061
若满足,则表示发生环路死锁,并输出L;若不满足,则返回步骤2)。
本发明公开了基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,每个栅格点至多只分配给一台移动机器人,能够有效避免移动机器人间发生碰撞。每台机器人每轮最多只能获得预先设置好的NC个栅格点,避免造成系统中大量机器人的长时间等待现象,实现了多机器人共享栅格地图的目的,有效提高了多机器人系统整体效率。在分配过程中采用前瞻一步的策略,实时构建系统前瞻一步后的阻塞图模型,应用相应的环路死锁检测与避免方法,有效避免了环路死锁现象。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,其特征在于,所述多个移动机器人共用一个栅格地图;所述分配方法包括步骤;
(1)获取栅格地图中栅格点的坐标、栅格点的状态、栅格地图中的障碍物位置、各移动机器人的位置、各移动机器人已规划的路径;根据获取的信息标识所有栅格点及其状态:记所有栅格点的集合为V={vj|j∈N且1≤j≤NV},其中,vj表示第j个栅格点,NV为栅格点总数;vj的状态标识为Occupy(vj),Occupy(vj)=0表示栅格点vj未被占用,Occupy(vj)=1表示栅格点vj已被占用,占用栅格点vj的机器人记为R(vj);
(2)记第i个移动机器人为Ri,i∈N,且1≤i≤NR,N表示正整数,NR表示移动机器人的总数;记Ri已规划的路径为Path(Ri),将系统根据Path(Ri)已分配给Ri的栅格点中最靠近Path(Ri)终点的栅格点作为Ri的临时停靠点,记为FinOcc(Ri);记从FinOcc(Ri)开始沿着路径Path(Ri)前向第k个栅格点为
Figure FDA0002726488770000011
(3)逐一为移动机器人分配栅格点,其中,为移动机器人Ri分配栅格点的步骤包括:
(3-1)初始化k=1;
(3-2)根据
Figure FDA0002726488770000012
的状态标识判断栅格点
Figure FDA0002726488770000013
是否已被其他移动机器人占据,若是,则直接进入步骤(3-6);否则,执行步骤(3-3);
(3-3)假定将栅格点
Figure FDA0002726488770000014
分配给机器人Ri,根据假定的分配结果构建阻塞图GQ(VQ,EQ),GQ(VQ,EQ)为有向图,VQ为顶点集,EQ为边集,
VQ={vj|vj∈VT}
Figure FDA0002726488770000015
其中,vj表示第j个栅格点,VT表示所有移动机器人临时停靠点的集合,ejk表示从栅格点vj到vk的单向路径;
(3-4)根据构建好的阻塞图,检测阻塞图中是否存在环路死锁,若不存在,则执行步骤(3-5),否则,执行步骤(3-6);
(3-5)将栅格点
Figure FDA0002726488770000021
分配给机器人Ri;计算k=k+1,判断是否满足k>Nc,若满足,则执行步骤(3-6),否则,返回步骤(3-2);其中,Nc为预设的每次为移动机器人分配栅格数的最大值;
(3-6)结束为移动机器人Ri分配栅格点;
(4)当为所有移动机器人分配过栅格点后,等待所有移动机器人按照各自路径顺序通过已分配栅格点,并停留在临时停靠栅格点后,返回步骤(3),重新对所有移动机器人进行栅格点分配。
2.根据权利要求1所述的一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,其特征在于,判断阻塞图中是否存在环路死锁的方法为:当阻塞图中存在有向回路时表示存在环路死锁,一个有向回路表示一个发生环路死锁的阻塞路径。
3.根据权利要求2所述的一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法,其特征在于,检索阻塞图中的环路死锁的方法为:
1)假定移动机器人Ri进入阻塞图GQ(VQ,EQ),初始化:
Figure FDA0002726488770000022
L={vcu},其中,vcu表示阻塞图中当前正在搜索的栅格点,L为发生环路死锁的阻塞路径;
2)判断阻塞图中是否存在vcu的前向邻节点Adj+(vcu);若不存在,则不会发生环路死锁;若存在,则更新:vcu=Adj+(vcu),L=L+{Adj+(vcu)};
3)判断是否满足
Figure FDA0002726488770000023
若满足,则表示发生环路死锁,并输出L;若不满足,则返回步骤2)。
CN201811546402.6A 2018-12-17 2018-12-17 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法 Active CN109508010B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811546402.6A CN109508010B (zh) 2018-12-17 2018-12-17 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811546402.6A CN109508010B (zh) 2018-12-17 2018-12-17 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法

Publications (2)

Publication Number Publication Date
CN109508010A CN109508010A (zh) 2019-03-22
CN109508010B true CN109508010B (zh) 2020-11-27

Family

ID=65752869

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811546402.6A Active CN109508010B (zh) 2018-12-17 2018-12-17 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法

Country Status (1)

Country Link
CN (1) CN109508010B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112083722B (zh) * 2020-08-27 2023-08-04 广州赛特智能科技有限公司 一种移动机器人多机多楼层调度系统及调度方法
CN112372640A (zh) * 2020-11-30 2021-02-19 倍仪昇智能科技(苏州)有限公司 多机械臂防冲突运行方法
CN112859847B (zh) * 2021-01-06 2022-04-01 大连理工大学 一种通行方向限制下的多机器人协同路径规划方法
CN114003011B (zh) * 2021-11-03 2023-08-15 盐城工学院 一种多载量agvs防死锁任务调度方法
CN114035522B (zh) * 2021-11-09 2024-07-23 南京航空航天大学 一种oht搬运系统防死锁调度方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005225664A (ja) * 2004-02-16 2005-08-25 Mitsubishi Heavy Ind Ltd 搬送車監視制御装置、搬送システム、及び搬送車監視制御方法
US9707681B2 (en) * 2015-07-27 2017-07-18 Siemens Industry Software Ltd. Anti-collision management of overlapping robotic movements
CN106774305B (zh) * 2016-11-30 2019-11-08 上海振华重工电气有限公司 自动化集装箱码头多自动引导车辆路径冲突消解方法

Also Published As

Publication number Publication date
CN109508010A (zh) 2019-03-22

Similar Documents

Publication Publication Date Title
CN109508010B (zh) 一种基于栅格地图的多移动机器人系统栅格点前瞻性防死锁动态分配方法
CN111638717B (zh) 一种分布式自主机器人交通协调机制的设计方法
CN108267149B (zh) 多移动机器人的冲突管理方法及系统
JP6671507B2 (ja) ロボット現場戻りのための方法及び装置{method and device for returning robots from site}
Silver Cooperative pathfinding
WO2019141222A1 (zh) 多移动机器人的冲突管理方法及系统
CN110471418B (zh) 智能停车场中的agv调度方法
WO2019141223A1 (zh) 移动机器人的路径规划方法及系统
CN110530369A (zh) 基于时间窗的agv任务调度方法
WO2019141219A1 (zh) 多移动机器人的调度方法及系统
CN111311959B (zh) 多路口协同控制方法、装置、电子设备及存储介质
WO2019141224A1 (zh) 多移动机器人的冲突管理方法及系统
WO2022032444A1 (zh) 一种多智能主体避障方法、系统和计算机可读存储介质
CN114815802A (zh) 一种基于改进蚁群算法的无人天车路径规划方法和系统
CN114489062B (zh) 面向车间物流的多自动导引车分布式动态路径规划方法
CN112835364A (zh) 一种基于冲突检测的多机器人路径规划方法
WO2023020213A1 (zh) 任务分配的方法、装置、设备、存储介质及程序产品
Xu et al. Dynamic spare point application based coordination strategy for multi-AGV systems in a WIP warehouse environment
CN113064436B (zh) 一种agv系统中动态路径规划和去中心化避障方法
Tai et al. A time-efficient approach to solve conflicts and deadlocks for scheduling AGVs in warehousing applications
CN117908475A (zh) 调度方法和调度装置
KR102397338B1 (ko) 자율 주차 로봇 시스템의 작업 계획 관리 장치 및 관리 방법
CN111553404B (zh) 一种通信网络节点自适应动态规划方法
CN111582579A (zh) 一种基于拓扑地图的预测性多agv任务分配方法
CN116700265A (zh) 一种基于改进cbs算法的多机器人路径规划方法及系统

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