CN109159122B - 采用椭圆型终态神经网络的冗余机器人重复运动规划方法 - Google Patents
采用椭圆型终态神经网络的冗余机器人重复运动规划方法 Download PDFInfo
- Publication number
- CN109159122B CN109159122B CN201811063399.2A CN201811063399A CN109159122B CN 109159122 B CN109159122 B CN 109159122B CN 201811063399 A CN201811063399 A CN 201811063399A CN 109159122 B CN109159122 B CN 109159122B
- Authority
- CN
- China
- Prior art keywords
- redundant
- neural network
- redundant robot
- joint
- robot
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
Abstract
一种基于椭圆型终态神经网络的冗余机器人重复运动规划方法,在笛卡尔空间中给定机器人末端执行器的期望轨迹rd(t),并给出各个关节的期望回拢角度θd(0);对于机器人的重复运动,采用渐近收敛性能指标,通过将冗余机器人轨迹规划的二次优化问题转化为时变矩阵方程求解问题,以椭圆型终态神经网络作为求解器。本发明在初始位置偏移的情况下,实现冗余机器人快速有限时间收敛的重复运动规划任务。本发明提供一种具有有限时间收敛、计算精度高、易于实现的基于椭圆型终态神经网络的冗余机器人运动规划方法。
Description
技术领域
本发明涉及工业机器人运动规划技术,具体地,提出一种可以有限时间收敛的、初始位置偏离期望轨迹情形下的、采用椭圆型终态神经网络的冗余机器人重复运动控制方法。
背景技术
工业机器人是面向工业领域的多关节机械手或多自由度的机器装置,它能自动执行工作,是靠自身动力和控制能力来实现各种功能的一种机器。工业机器人技术是近年来高新技术发展的重要领域之一。这一技术在工业、农业、国防、医疗卫生、办公自动化及生活服务等众多领域有着越来越多的应用。
机器人的分类方法很多,按照应用可以将机器人分为工业机器人、极限作业机器人和娱乐机器人;按照自由度的多少可以分为冗余和非冗余机器人。冗余机器人是指末端执行器在执行给定的任务时比其所必需的自由度有更多的自由度。冗余机器人具有更多的操作空间,动作更灵活。考虑n自由度冗余机器人各关节角度和末端执行器位移的关系
r(t)=f(θ(t)) (1)
其中,r(t)是末端执行器在笛卡尔坐标系下的位姿变量,θ(t)表示关节角度。末端笛卡尔空间和关节空间的微分之间的关系为
其中,J(θ(t))+=J(θ(t))T(J(θ(t))J(θ(t))T)-1是雅克比矩阵J(θ(t))的伪逆。
具有等式约束的最小速度范数性能指标作为运动规划的目标函数(D.E.Whitney,Resolved motion rate control of manipulators and human prostheses(操纵器和人工假肢的运动速率控制),IEEE Trans.Man-Machine Syst.1969,10(2):47-53):
式中,A(t)为正定加权矩阵。建立拉格朗日函数
解得
基于二次优化(Quadratic Optimization,QP)的冗余解析方案受到关注(F.T.Cheng,T.-H.Chen,and Y.Y.Sun,Resolving manipulator redundancy underinequality constraints(不等式约束条件下的冗余机器人轨迹规划方法),IEEE Trans.Robotics Automat1994,10(1):65-71;):
冗余机器人轨迹规划的性能直接关系到冗余机器人能否实现给定的末端任务。当末端执行器的运动轨迹是闭合的,在冗余机器人完成末端工作任务后,各个关节角变量在运动空间中的轨迹不一定封闭。这种非重复性问题可能产生不期望的关节位形,使得冗余机器人末端封闭轨迹的重复作业出现预料之外的情况,甚至会导致意外及危险情况的发生。应用最为广泛的伪逆控制方法不能保证运动的重复性。为了完成重复运动,通常采用自运动的方法进行弥补,而自运动进行调整往往效率不高(详见Klein C A and Huang C,Review of Pseudo Inverse Control for use with Kinematically RedundantManipulators(基于伪逆控制方法的冗余机器人运动规划).IEEETrans.Syst.Man.Cybern.1983,13(2):245-250;Tchon K, Janiak M.Repeatableapproximation of the Jacobian pseudo-inverse(雅可比伪逆的可重复逼近).Systemsand Control Letters,2009,58(12):849-856)。
为了执行重复运动任务,引入重复运动指标作为优化准则,形成重复运动规划(Repetitive motion planning,RMP)方案(Zhang Y,Wang J,Xia Y.A dual neuralnetwork for redundancy resolution of kinematically redundant manipulatorssubject to joint limits and joint velocity limits(基于关节角度和角速度限制的冗余机器人轨迹规划方法).IEEE Trans Neural Netw,2003,14(3):658-667)。常用的重复运动指标为如下渐近收敛性能指标AOC(Asympototically-Convengent OptimalityCriterion):
这里的二次规划(QP)采用递归神经网络(RNN)计算方法求解。通常的神经网络求解器具有渐近收敛性能,只要计算时间足够长,能够获得有效解。
基于二次规划描述冗余度解析问题,已发表文献往往采用递归神经网络求解。相比于具有渐近收敛动态特性的递归神经网络,基于椭圆型终态神经网络的动态特性具有有限时间收敛性,收敛速度快、计算精度高。值得指出的是,冗余机器人轨迹规划问题归结为时变计算问题,而求解时变问题的有效方法是采用具有有限时间收敛性能的递归神经网络计算方法。另外,已发表文献中的有限时间收敛的神经网络多采用线性激励函数,或具有无限值激励函数,实际实现时,由于能量有限,无限值激励函数神经网络实现时也存在本质困难。
发明内容
为了克服现有冗余机器人重复运动规划方法的不能在有限时间收敛、计算精度较低、不易实现的不足,本发明提供一种具有有限时间收敛、计算精度高、易于实现的基于椭圆型终态神经网络的冗余机器人运动规划方法。本发明采用渐近收敛性能指标,通过将冗余机器人轨迹规划的二次优化问题转化为时变矩阵方程求解问题,以基于椭圆型终态神经网络作为求解器,在初始位置偏移的情况下,实现冗余机器人快速有限时间收敛的重复运动规划任务。
为了实现这一目的,本发明给出如下的技术方案。
一种基于椭圆型终态神经网络的冗余机器人重复运动规划方法,包括以下步骤:
1)在笛卡尔空间中给定冗余机器人末端执行器的期望的运动轨迹rd(t),并给出各个关节的期望关节角度θd(0);
2)对于重复运动的冗余机器人,其初始关节角度为θ(t)=θ0,初始关节角度θ0不同于期望关节角度,即θ0≠θd(0);
3)将冗余机器人重复运动规划描述为二次规划问题,给出渐近收敛性能指标,形成重复运动规划方案:
其中,g(θ)=-βθ(θ(t)-θd(0)),θ(t)、分别表示冗余机器人的关节角度和关节角速度,βθ是一设计参数,用来形成关节位移的动态性能;θ(t)-θd(0)表示各个关节角与各个关节的期望关节角度之间的关节角位移偏差,rd(t)表示冗余机器人末端执行器期望的运动轨迹,表示冗余机器人末端执行器期望的速度向量。由于冗余机器人的初始位置可能不在期望的运动轨迹上,通过减小冗余机器人末端执行器期望的运动轨迹与实际运动轨迹位置间的误差rd(t)-f(θ(t)),改变冗余机器人末端执行器的运动方向,βr>0表示位置的参数增益,用来调节冗余机器人末端执行器运动到期望的运动轨迹的速率;J(θ(t))是冗余机器人雅可比矩阵,f(θ(t))是冗余机器人实际运动轨迹;
4)构建如下描述椭圆型终态神经网络的动态特性方程
其中,
Eij(t)为误差变量,ε>0为用于调整误差吸引速度的常值,σ>0为F(Eij(t),σ)分段函数分段边界;
式(10)是有限区间稳定的,其有限时间收敛性需分两种情况说明:
当|Eij(t)|≤σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
当|Eij(t)|>σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
5)定义拉格朗日函数
W(t)Y(t)=V(t) (14)
其中,
I为具有相应维数单位矩阵;
6)为求解步骤3)中的二次规划问题,由式(14)定义误差
E(t)=W(t)Y(t)-V(t) (15)
根据椭圆型终态神经网络动态特性方程式(10),给出椭圆型终态神经网络方程
通过神经网络计算过程,得到冗余机器人各个关节角度。
本发明的技术构思为:采用渐近收敛的性能指标
其中,
g(θ(t))=- βθ(θ(t)-θd(0))
βr>0,θ(t)-θd(0)表示各个关节角与各个关节的期望关节角度之间的关节角位移偏差,由于人为构造的收敛关系,这个误差会不断地缩小,直至为零。参数βr用来调节冗余机器人末端执行器运动到期望的运动轨迹的速率。
对冗余机器人重复运动规划构建一种椭圆型终态神经网络的动态特性方程
其中,
Eij(t)为误差变量,ε>0为用于调整误差吸引速度的常值,σ>0为F(Eij(t),σ)分段函数分段边界。
式(10)是有限区间稳定的,其有限时间收敛性需分两种情况说明:
当|Eij(t)|≤σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
当|Eij(t)|>σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
当Eij(t)有限时间收敛于零时,表示冗余机器人末端执行器沿着期望的运动轨迹在运动。
本发明的有益效果主要表现在:将一种椭圆型终态神经网络方法应用于冗余机器人重复运动规划中,在初始位置偏移情况下实现冗余机器人的有限时间收敛重复运动规划任务。相比于已有重复运动规划方法,该方法误差收敛快,且冗余机器人各个关节偏离期望关节角度的误差精度高。相比于具有渐近收敛特性的递归神经网络,椭圆型终态神经网络的有限时间收敛特性适用于时变问题求解(式 (14)为时变矩阵方程);该神经网络求解器采用的网络动态方程(10)是有限值函数,在工程应用中易于实现,且成本低,符合工程实际需要。
附图说明
图1为本发明提供的重复运动规划方案的流程图。
图2取不同σ值时的函数F(Eij(t),σ)。
图3为本发明重复运动规划方案的冗余机器人PA10。(摘自Y Zhang,Z Zhang.Repetitive motion planning and control of redundant robot manipulators(冗余机器人的重复运动规划与控制).Springer Berlin Heidelberg,2013:23-24)
图4为冗余机器人PA10末端执行器的运动轨迹。
图5为冗余机器人PA10各个关节的运动轨迹,其中,1为PA10机械臂基座,2为肘部轨迹,3为手腕轨迹,4为末端执行器运动轨迹。
图6为冗余机器人PA10各个关节角。
图7为冗余机器人PA10各个关节速度。
图8为冗余机器人PA10末端执行器的各位置误差。
图9为以椭圆型终态神经网络和渐近神经网络求解时的误差。
具体实施方式
下面结合附图对本发明做进一步描述。
参照图1~图9,一种基于椭圆型终态神经网络的冗余机器人重复运动规划方法,图1为冗余机器人重复运动规划方案的流程图,由以下三个步骤组成:1、确定冗余机器人末端执行器期望的运动轨迹和期望关节角度;2、采用渐近收敛性能指标,并形成冗余机器人重复运动二次规划方案;3、以基于椭圆型终态神经网络求解二次规划问题,获得各关节角轨迹,包括以下步骤:
1)确定期望的运动轨迹
设定冗余机器人PA10的期望关节角度确定圆轨迹的圆心坐标(x=0.2m,y=0m,z=0m),将圆的半径设为0.2m,其圆面与平面的夹角为末端执行器完成圆轨迹时间为T=10s。考虑到冗余机器人的初始位置不在期望的运动轨迹上,将冗余机器人的七个关节角度初值设为
2)建立冗余机器人重复运动的二次规划方案
将冗余机器人重复运动规划描述为二次规划问题,给出渐近收敛性能指标,形成重复运动规划方案:
其中,g(θ)=-βθ(θ(t)-θd(0)),θ(t)、分别表示冗余机器人的关节角度和关节角速度,βθ是一设计参数,用来形成关节位移的动态性能。θ(t)-θd(0)表示各个关节角与各个关节的期望关节角度之间的关节角位移偏差,rd(t)表示冗余机器人末端执行器期望的运动轨迹,表示冗余机器人末端执行器期望的速度向量。由于冗余机器人的初始位置可能不在期望的轨迹上,通过减小冗余机器人末端执行器期望的运动轨迹与实际运动轨迹位置间的误差rd(t)-f(θ(t)),改变冗余机器人末端执行器的运动方向,βr>0表示位置的参数增益,用来冗余机器人调节末端执行器运动到期望的运动轨迹的速率。J(θ(t))是冗余机器人雅可比矩阵,f(θ(t))是冗余机器人实际运动轨迹。
3)构造椭圆型终态神经网络求解上述二次规划问题
其中,
Eij(t)为误差变量,ε>0为用于调整误差吸引速度的常值,σ>0为F(Eij(t),σ)分段函数分段边界。
式(10)是有限区间稳定的,其有限时间收敛性需分两种情况说明:
当Eij(0)≤σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
当Eij(0)>σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
定义拉格朗日函数
W(t)Y(t)=V(t) (14)
其中,
I为具有相应维数单位矩阵。
为求解二次规划问题,由式(14)定义误差
E(t)=W(t)Y(t)-V(t) (15)
根据椭圆型终态神经网络动态特性方程式(10),给出椭圆型终态神经网络方程
通过神经网络计算过程,得到冗余机器人各个关节角度。
图2为取不同σ值时的F(Eij(t),σ)。从图中可以看出,随着σ减小,该函数导数趋于无穷大,从而使得神经网络终态收敛。
图3所示的冗余机器人PA10用于实现本发明的重复运动规划方案。该冗余机器人有七个自由度(四个旋转轴和三个枢轴),PA10机械手的主要规格是机械手的全长为1.345m,机械手的重量为343N,所有轴的最大组合速度为1.55m/s,有效载荷重量为98N。臂的长度为d3=0.45m,d5=0.5m,d7=0.08m。
图4为冗余机器人末端执行器在空间中的运动轨迹。图中给出目标圆轨迹及冗余机器人末端执行器运动轨迹。可以看出,冗余机器人末端执行器的初始位置不在期望的运动轨迹上。如图8所示,随着时间增加(T=10s),冗余机器人末端执行器的终值位置误差精度在XYZ轴三个方向上达到4×10-5,实际轨迹与期望的运动轨迹吻合。
图5所示为冗余机器人各个关节的运动轨迹。从图中可以看出,各个关节的轨迹在运行一个周期后是闭合的,实现了重复运动控制。
为了验证基于椭圆型终态神经网络方法在重复运动规划中的有效性,冗余机器人PA10末端执行器完成圆轨迹过程中得到的关节角瞬态轨迹和关节角速度瞬态轨迹如图6、7所示。可以看出,冗余机器人的各关节角最终收敛于期望的关节角位置,关节角速度收敛于0。
当T=10s时,采用椭圆型终态神经网络求解得到的冗余机器人运动前后关节角与其期望关节角之间的最大偏差为9.2608×10-5rad,用渐近神经网络求解得到的各个关节角的终值误差最大偏差为6.8740×10-4rad,详见表1:
表1 其中,渐近神经网络动态特性
椭圆型终态神经网络动态特性
其中,
图9为将基于椭圆型终态的神经网络与基于渐近神经网络的冗余机器人重复轨迹规划误差的F范数进行对比,从图中可以看出,当误差的F范数||E||F收敛至0.01时,基于椭圆型吸引律的神经网络需要时间T(||Eij(t)||F=0.01)=0.650s;基于渐近神经网络需要时间T(||Eij(t)||F=0.01)=6.932s。
Claims (1)
1.一种基于椭圆型终态神经网络的冗余机器人重复运动规划方法,其特征在于,所述方法包括以下步骤:
1)在笛卡尔空间中给定冗余机器人末端执行器的期望的运动轨迹rd(t),并给出各个关节的期望关节角度θd(0);
2)对于重复运动的冗余机器人,其初始关节角度为θ(0)=θ0,初始关节角度θ0不同于期望关节角度,即θ0≠θd(0);
3)将冗余机器人重复运动规划描述为二次规划问题,给出渐近收敛性能指标,形成重复运动规划方案:
其中,g(θ)=-βθ(θ(t)-θd(0)),θ(t)、分别表示冗余机器人的关节角度和关节角速度,βθ是一设计参数,用来形成关节位移的动态性能;θ(t)-θd(0)表示各个关节角与各个关节的期望关节角度之间的关节角位移偏差,rd(t)表示冗余机器人末端执行器期望的运动轨迹,表示冗余机器人末端执行器期望的速度向量;由于冗余机器人的初始位置可能不在期望的运动轨迹上,通过减小冗余机器人末端执行器期望的运动轨迹与实际运动轨迹位置间的误差rd(t)-f(θ(t)),改变冗余机器人末端执行器的运动方向,βr>0表示位置的参数增益,用来调节冗余机器人末端执行器运动到期望的运动轨迹的速率;J(θ(t))是冗余机器人雅可比矩阵,f(θ(t))是冗余机器人实际运动轨迹;
4)构建如下描述椭圆型终态神经网络的动态特性方程
其中,
Eij(t)为误差变量,ε>0为用于调整误差吸引速度的常值,σ>0为F(Eij(t),σ)分段函数分段边界;
式(10)是有限区间稳定的,其有限时间收敛性需分两种情况说明:
当|Eij(t)|≤σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
当|Eij(t)|>σ时,从Eij(0)到Eij(t)=0变化,收敛时间满足
5)定义拉格朗日函数
W(t)Y(t)=V(t) (14)
其中,
I为具有相应维数单位矩阵;
6)为求解步骤3)中的二次规划问题,定义误差
E(t)=W(t)Y(t)-V(t) (15)
根据椭圆型终态神经网络动态特性方程式(10),给出椭圆型终态神经网络方程
通过神经网络计算过程,得到冗余机器人各个关节角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811063399.2A CN109159122B (zh) | 2018-09-12 | 2018-09-12 | 采用椭圆型终态神经网络的冗余机器人重复运动规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811063399.2A CN109159122B (zh) | 2018-09-12 | 2018-09-12 | 采用椭圆型终态神经网络的冗余机器人重复运动规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109159122A CN109159122A (zh) | 2019-01-08 |
CN109159122B true CN109159122B (zh) | 2021-01-01 |
Family
ID=64894776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811063399.2A Active CN109159122B (zh) | 2018-09-12 | 2018-09-12 | 采用椭圆型终态神经网络的冗余机器人重复运动规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109159122B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110695994B (zh) * | 2019-10-08 | 2021-03-30 | 浙江科技学院 | 一种面向双臂机械手协同重复运动的有限时间规划方法 |
CN111230860B (zh) * | 2020-01-02 | 2022-03-01 | 腾讯科技(深圳)有限公司 | 机器人控制方法、装置、计算机设备及存储介质 |
CN111185909B (zh) * | 2020-01-14 | 2022-03-18 | 深圳众为兴技术股份有限公司 | 机器人运行工况获取方法、装置、机器人及存储介质 |
CN117506934B (zh) * | 2024-01-03 | 2024-03-12 | 创新奇智(北京)科技有限公司 | 工业机器人自动化控制系统及方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107595392B (zh) * | 2012-06-01 | 2020-11-27 | 直观外科手术操作公司 | 使用零空间回避操纵器臂与患者碰撞 |
CN106456265B (zh) * | 2014-03-17 | 2019-07-30 | 直观外科手术操作公司 | 使用逆运动学在关节极限下的遥控操作手术系统和控制方法 |
CN106985138B (zh) * | 2017-03-13 | 2019-05-31 | 浙江工业大学 | 基于终态吸引优化指标的冗余机械臂轨迹规划方法 |
CN107127754A (zh) * | 2017-05-09 | 2017-09-05 | 浙江工业大学 | 一种基于终态吸引优化指标的冗余机械臂重复运动规划方法 |
CN107972031B (zh) * | 2017-11-10 | 2020-07-07 | 浙江科技学院 | 一种冗余机械臂可重复运动的初始位置定位方法 |
CN107962566B (zh) * | 2017-11-10 | 2021-01-08 | 浙江科技学院 | 一种移动机械臂重复运动规划方法 |
CN108098777B (zh) * | 2018-01-12 | 2021-04-30 | 华侨大学 | 一种冗余度机械臂力矩层重复运动控制方法 |
-
2018
- 2018-09-12 CN CN201811063399.2A patent/CN109159122B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109159122A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109159122B (zh) | 采用椭圆型终态神经网络的冗余机器人重复运动规划方法 | |
CN107962566B (zh) | 一种移动机械臂重复运动规划方法 | |
CN107972030B (zh) | 一种冗余机械臂重复运动中的初始位置修正方法 | |
CN107891424B (zh) | 求解冗余机械臂逆运动学的有限时间神经网络优化方法 | |
CN107972031B (zh) | 一种冗余机械臂可重复运动的初始位置定位方法 | |
CN106985138B (zh) | 基于终态吸引优化指标的冗余机械臂轨迹规划方法 | |
Yin et al. | Direct adaptive robust tracking control for 6 DOF industrial robot with enhanced accuracy | |
CN108908347B (zh) | 一种面向冗余移动机械臂容错型重复运动规划方法 | |
CN109940615B (zh) | 一种面向双臂机械手同步重复运动规划的终态网络优化方法 | |
Vasu | Fuzzy PID based adaptive control on industrial robot system | |
CN109159124B (zh) | 采用快速双幂次终态神经网络的冗余机器人重复运动规划方法 | |
CN108908340B (zh) | 采用有限区间神经网络的冗余机器人重复运动规划方法 | |
CN110695994B (zh) | 一种面向双臂机械手协同重复运动的有限时间规划方法 | |
CN109159121B (zh) | 采用抛物线型终态神经网络的冗余机器人重复运动规划方法 | |
CN109015657B (zh) | 一种面向冗余机械臂重复运动规划的终态网络优化方法 | |
Liu et al. | Research on kinematic modeling and analysis methods of UR robot | |
CN108908341B (zh) | 二次根式终态吸引的冗余机器人重复运动规划方法 | |
Jun et al. | Continuous terminal sliding mode control of a 6-DOF wire-driven parallel robot | |
Li et al. | Fault-tolerant motion planning of redundant manipulator with initial position error | |
Si et al. | A complete solution to the inverse kinematics problem for 4-DOF manipulator robot | |
Ren et al. | Adaptive synchronization control of a planar parallel manipulator | |
Zhang et al. | Adaptive kinematic control of redundant manipulators | |
Xiang et al. | An optimal trajectory control strategy for underwater welding robot | |
Yao et al. | A motion simulation for dual-arm robot based on the state-space method | |
Zeng et al. | Time-optimal trajectory planning based on particle swarm optimization |
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 |