CN109623827B - 一种高性能冗余度机械臂重复运动规划方法与装置 - Google Patents

一种高性能冗余度机械臂重复运动规划方法与装置 Download PDF

Info

Publication number
CN109623827B
CN109623827B CN201910055288.5A CN201910055288A CN109623827B CN 109623827 B CN109623827 B CN 109623827B CN 201910055288 A CN201910055288 A CN 201910055288A CN 109623827 B CN109623827 B CN 109623827B
Authority
CN
China
Prior art keywords
redundant manipulator
joint
repetitive motion
quadratic
representing
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
CN201910055288.5A
Other languages
English (en)
Other versions
CN109623827A (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.)
Lanzhou University
Original Assignee
Lanzhou 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 Lanzhou University filed Critical Lanzhou University
Priority to CN201910055288.5A priority Critical patent/CN109623827B/zh
Publication of CN109623827A publication Critical patent/CN109623827A/zh
Application granted granted Critical
Publication of CN109623827B publication Critical patent/CN109623827B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1602Programme controls characterised by the control system, structure, architecture
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)

Abstract

本发明提供了一种高性能冗余度机械臂重复运动规划方法与装置,包括如下步骤:1)设计二次型的重复运动性能指标,生成二次型优化冗余度解析方案;2)将步骤1)的二次型优化冗余度解析方案转化为二次规划的形式;3)将步骤2)的二次规划运用二次规划求解器进行求解;4)将步骤3)的求解结果传递给下位机控制器驱动冗余度机械臂运动。本发明克服了现有的重复运动方案中当冗余度机械臂的关节角度误差减小时,末端执行器的位置误差会增大的缺点,并首次在理论上消除了理论位置误差,显著提高了冗余度机械臂执行任务的精确度,对精密仪器的控制有着很大的意义,这在工程应用领域具有良好前景。

Description

一种高性能冗余度机械臂重复运动规划方法与装置
技术领域
本发明涉及冗余度机械臂运动规划及控制领域,具体涉及一种冗余度机械臂的逆运动学求解方法。
背景技术
冗余度机械臂具有自由度大于任务空间所需最少自由度的特性,广泛应用于装备制造、产品加工、机器作业等国民经济生产活动中。冗余度机械臂的逆运动学问题是指已知冗余度机械臂末端位置,确定冗余度机械臂的关节角问题。当冗余度机械臂末端任务为一个封闭曲线时,其各个关节可能回不到初始位置,这种现象叫做关节角偏差,或称非重复运动问题。如果冗余度机械臂不能实现重复运动,则在运行过程中可能会产生不可预料的情况,甚至出现损害冗余度机械臂以及危害周围人员人身安全的问题;而重复运动规划方案就是要设计适当的性能参数,使得冗余度机械臂末端执行完封闭曲线任务时,各个关节角都能够回到其初始位置。
现有重复运动解析方法全都存在理论位置误差,而理论位置误差直接影响到了轨迹追踪任务的完成精度。具体而言,现有冗余度机械臂重复运动中的理论位置误差和关节角度误差存在耦合关系,由此会产生关节角度误差和理论位置误差之间的矛盾。例如:现有方法《一种冗余度机械臂重复运动规划的方法》(申请号:201010264141.6)中,关节角度反馈的增加会造成位置误差的增大,从而降低的精度。为了克服现有重复运动方案中出现理论位置误差偏大,从而对冗余度机械臂执行任务精度产生影响这一缺陷,本发明提出一种高性能冗余度机械臂重复运动规划方法与装置,有效地消除了执行任务时产生的理论误差,并能够更精确地完成冗余度机械臂的重复运动任务,这一点对于执行精密任务是非常有意义的。综上所述,本发明专利具有新颖性与实用性。
发明内容
本发明所要解决的技术问题是提供一种高性能冗余度机械臂重复运动规划方法与装置。
为解决上述技术问题,本发明通过以下技术方案予以实现:
1、一种高性能冗余度机械臂重复运动规划方法与装置包括如下步骤:
1)通过上位机采用二次型优化方案对冗余度机械臂的重复运动进行进行冗余度解析,设计的最小化性能指标为重复运动,受约束于雅可比等式,关节角度极限,关节角速度极限;
2)将步骤1)的二次型优化方案转化为二次规划;
3)将步骤2)的二次规划用二次规划求解器进行求解;
4)将步骤3)的求解结果传递给下位机控制器驱动冗余度机械臂运动。
2、一种高性能冗余度机械臂重复运动规划方法与装置设计为:最小化
Figure BDA0001952232630000021
受约束于雅可比等式
Figure BDA0001952232630000022
关节角度极限0-≤θ≤θ+,关节角速度极
Figure BDA0001952232630000023
上标T表示矩阵和向量的转置,上标
Figure BDA0001952232630000026
表示矩阵的伪逆,θ表示冗余度机械臂的关节角度,θ0表示冗余度机械臂的初始关节角度,
Figure BDA0001952232630000024
表示冗余度机械臂的关节角速度,κ>0表示关于重复运动的反馈系数,rd表示设定的末端执行器的期望轨迹,
Figure BDA0001952232630000025
表示设定的末端执行器期望速度,J表示冗余度机械臂的雅可比矩阵,I表示单位矩阵,
Figure BDA0001952232630000031
表示本发明中不同于现存方法的重要特征,起到促使位置误差与关节角度误差解耦的作用,θ±表示冗余度机械臂的关节角度极限,
Figure BDA0001952232630000032
表示冗余度机械臂的关节角速度的上下限;
3、上述重复运动的二次型优化方案可以转化为一个二次规划,即最小化
Figure BDA0001952232630000033
受约束于
Figure BDA0001952232630000034
其中
Figure BDA0001952232630000035
Figure BDA0001952232630000036
κ>0表示关于重复运动的反馈系数,β±为冗余度机械臂的关节速度上下限,它们的第i个元素分别定义为
Figure BDA0001952232630000037
i表示关节序号,i=1,2,3...n,正常数kα用来调节和保证关节速度有足够大的可行域。
4、此二次规划可以用二次规划求解器求解。
5、二次规划求解器求解结果转化为电机驱动所需要的控制信号,从而驱动各关节电机使冗余度机械臂重复运动。
6、通过控制装置的各模块之间的信息传递完成从求解器到冗余度机械臂重复运动的控制;
所述控制装置包括:
信息传递模块,用于获取二次规划求解器求解结果,并将求解结果传递给下位机控制冗余度机械臂的重复运动生成;
轨迹信息获取模块,用于获取所述目标冗余度机械臂设定的轨迹信息,并将轨迹信息传输给二次规划求解器;
等式构建模块,用于根据所述位置和速度信息,以及所述参考点轨迹信息,按照预设规则构建所述目标冗余度机械臂对应的雅可比矩阵等式;
控制信号确定模块,用于在所述雅可比矩阵等式、关节角度极限和关节角速度极限的约束下,根据二次型优化和标准二次规划方法确定所述目标冗余度机械臂的控制信号;
冗余度机械臂控制模块,用于根据所述控制信号控制目标冗余度机械臂,使目标冗余度机械臂实现重复运动。
与现有技术相比,本发明有如下优点:
现有的冗余度机械臂的重复运动方案都存在理论位置误差,理论位置误差与关节角度误差是成反比的关系,即关节角度误差越小,位置误差越大。本发明首次使末端执行器的位置误差与冗余度机械臂的关节角度误差解耦,并成功地消除了理论位置误差,令冗余度机械臂的重复运动能够达到更加精确的效果,这一点对于精密仪器是非常有意义的,这在实际应用中有重要作用和广阔前景。
附图说明
图1为本发明的流程图;
图2为实现本发明的冗余度机械臂三维模型图;
图3为实现本发明应用下的冗余度机械臂重复运动示意图;
图4为实现本发明应用下的冗余度机械臂重复运动过程图;
图5为现有重复运动方案应用下的冗余度机械臂重复运动末端执行器位置误差图和运动关节角度变化图;
图6为实现本发明应用下的冗余度机械臂重复运动末端执行器位置误差图和运动关节角度变化图;
具体实施方式
下面结合附图对本发明做进一步的说明。
图1所示的一种高性能冗余度机械臂重复运动规划方法与装置首先提出重复运动性能指标与约束;然后将其转化为相应的二次规划;接着应用二次规划求解器求解;最后下位机控制器根据求解结果控制冗余度机械臂的运动。
图2所示实现本发明的冗余度机械臂为一个空间六自由度的PUMA560机械臂。该冗余度机械臂由六个连杆组成,通过关节1、关节2、关节3、关节4、关节5和关节6组成。
本发明设计的重复运动规划方案为最小化:
Figure BDA0001952232630000051
约束条件:
Figure BDA0001952232630000052
θ-≤θ≤θ+   (3)
Figure BDA0001952232630000053
其中,
Figure BDA0001952232630000054
为重复运动的性能指标,上标T表示矩阵和向量的转置,上标
Figure BDA0001952232630000055
表示矩阵的伪逆,θ表示冗余度机械臂的关节角度,θ0表示冗余度机械臂的初始关节角度,
Figure BDA0001952232630000056
表示冗余度机械臂的关节角速度,κ>0表示关于重复运动的反馈系数,
Figure BDA0001952232630000057
表示设定的期望轨迹,J表示冗余度机械臂的雅可比矩阵,I表示单位矩阵,
Figure BDA0001952232630000058
表示冗余度机械臂的关节角速度的上下限,θ±表示冗余度机械臂的关节角度的上下限。
带物理约束的冗余度机械臂的二次型优化方案(1)-(4)便可描述为如下的二次规划方案:
最小化:
Figure BDA0001952232630000061
约束条件:
Figure BDA0001952232630000062
Figure BDA0001952232630000063
其中,
Figure BDA0001952232630000064
κ>0表示关于重复运动的反馈系数,
Figure BDA0001952232630000065
Figure BDA0001952232630000066
β±为冗余度机械臂的关节速度上下限,它们的第i个元素分别定义为
Figure BDA0001952232630000067
Figure BDA0001952232630000068
i表示关节序号,i=1,2,3...n,正常数kα用来调节和保证关节速度足够大的可行域。下面给出一个可行的求解器的例子(本发明不限于该求解器例子),通过使用梯度下降法加上速度补偿的方法,上述问题又可以等价为分段线性投影方程组:
Figure BDA0001952232630000069
Figure BDA00019522326300000610
其中,γ>0表示末端执行器位置误差反馈系数,r表示末端执行器实际轨迹,rd-r表示末端执行器的位置误差,w表示连接(8)和(9)的辅助变量,JTw表示速度补偿,v>0表示式子(9)的系数,
Figure BDA00019522326300000611
本发明的理论误差等于零,并且与关节误差解耦。然而对于现有的方案如《一种冗余度机械臂重复运动规划的方法》(申请号:201010264141.6)位置误差为
Figure BDA00019522326300000612
通过这个位置误差表达式,可以看出对于现有的重复运动方案当增加了关节角度误差反馈系数k时,关节角度误差会减小,可是位置误差会增大,这种矛盾的关系是现存冗余度机械臂重复运动方案的不足之处。由此可见,现存的冗余度机械臂重复运动规划方案理论上存在与关节角度误差有关的位置误差,然而所提出的一种高性能冗余度机械臂重复运动规划方法与装置不存在理论上的位置误差。对本发明进行求解器的求解,从而得到高性能冗余度机械臂重复运动规划的精确解。
通过上位机的二次规划求解器得到该二次规划的解后,再将结果转化为电机驱动所需要的控制信号,并通过该冗余度机械臂的控制器来驱动冗余度机械臂各关节电机运动,从而驱动各关节电机使冗余度机械臂重复运动。
现结合一个具体的实例操作对本发明的工作流程进行如下说明。
利用MATLAB软件对本发明方法进行六自由度的PUMA560冗余度机械臂重复运动轨迹追踪的实验仿真,并在相同参数设置与现存的方法进行结果对比。具体参数设置如下:γ=2000,k=10,v=0.001,冗余度机械臂的关节速度极限为β+=-β-=[0.9,0.9,0.9,0.9,0.9,0.9]弧度/秒,计算机仿真时间设为10秒,执行任务为三尖瓣图形,初始的角度设置为q0=[0,0,0.5,0,0,0]弧度。通过二次规划求解器求解,将计算得到的结果再传送给冗余度机械臂控制器从而控制冗余度机械臂的运动。结果相比之下,一种高性能冗余度机械臂重复运动规划方法与装置在完成冗余度机械臂轨迹任务时位置误差更小,能更加精确地完成任务。
图3为实现本发明应用下的冗余度机械臂重复运动示意图。从图3可以看出来所给定的三尖瓣轨迹追踪任务顺利被完成。
图4为实现本发明应用下的冗余度机械臂重复运动过程图。子图1为初始冗余度机械臂的状态,子图2,3,4为中间过程中冗余度机械臂的状态,子图5为一次任务循环结束时冗余度机械臂的状态,子图6为一个新的重复运动循环的冗余度机械臂的状态。图4体现了在一个任务循环中,结束时冗余度机械臂的状态和初始冗余度机械臂状态完全相同,最终完成了重复运动任务。
图5为现有重复运动方案应用下的冗余度机械臂重复运动末端执行器位置误差图和运动关节角度变化图。从图5看出,当角度反馈系数较大时,现有重复运动方案应用《一种冗余度机械臂重复运动规划的方法》(申请号:201010264141.6)中的冗余度重复运动末端执行器在整个运动过程中位置误差比较大,最大误差甚至接近0.01米,会影响任务完成的质量。另外,冗余度机械臂重复运动关节角度结束状态和初始状态相同,满足重复运动的要求。
图6为实现本发明应用下的冗余度机械臂重复运动末端执行器位置误差图和运动关节角度变化图。从图6看出,当角度反馈系数较大时,本发明应用下的冗余度机械臂重复运动末端执行器位置误差比较小,并整体保持在较小的状态下,最大位置误差小于2×10-4米。另外,冗余度机械臂重复运动关节角度结束状态和初始状态相同,满足重复运动的要求。

Claims (3)

1.一种高性能冗余度机械臂重复运动规划方法,其特征在于,设计过程包括如下步骤:
1)通过上位机采用二次型优化方案对冗余度机械臂的重复运动进行冗余度解析,设计的最小化性能指标为重复运动,受约束于雅可比等式、关节角度极限、关节角速度极限,即为最小化
Figure FDA0004086230200000011
受约束于雅可比等式
Figure FDA0004086230200000012
关节角度极限θ-≤θ≤θ+,关节角速度极限
Figure FDA0004086230200000013
上标T表示矩阵和向量的转置,上标
Figure FDA00040862302000000113
表示矩阵的伪逆,θ表示冗余度机械臂的关节角度,θ0表示冗余度机械臂的初始关节角度,
Figure FDA0004086230200000014
表示冗余度机械臂的关节角速度,κ>0表示关于关节角度误差的反馈系数,rd表示设定的末端执行器的期望轨迹,
Figure FDA0004086230200000015
表示设定的末端执行器的期望速度,J表示冗余度机械臂的雅可比矩阵,I表示单位矩阵,
Figure FDA0004086230200000016
起到使末端执行器的位置误差与冗余度机械臂的关节角度误差解耦的作用,θ±表示关节角度的上下限,
Figure FDA0004086230200000017
表示冗余度机械臂的关节角速度的上下限;
2)将步骤1)的重复运动的二次型优化方案转化为二次规划,设计其性能指标为最小化
Figure FDA0004086230200000018
受约束于
Figure FDA0004086230200000019
其中
Figure FDA00040862302000000110
q=θ,q0=θ0
Figure FDA00040862302000000111
κ>0表示关于关节角度误差的反馈系数,β±为冗余度机械臂的关节速度上下限,它们的第i个元素分别定义为
Figure FDA00040862302000000112
i表示关节序号,i=1,2,3...n,正常数kα用来调节和保证关节速度有足够大的可行域;
3)将步骤2)的二次规划用二次规划求解器进行求解;
4)将步骤3)的求解结果传递给下位机控制器驱动冗余度机械臂运动。
2.根据权利要求1所述的一种高性能冗余度机械臂重复运动规划方法,其特征在于,本运动规划方法是由所述步骤3)的二次规划求解器求解。
3.根据权利要求2所述的一种高性能冗余度机械臂重复运动规划方法,其特征在于,可将所述步骤3)的二次规划求解器求解结果转化为电机驱动所需要的控制信号,从而在消除理论位置误差的情况下,驱动各关节电机使冗余度机械臂重复运动。
CN201910055288.5A 2019-01-21 2019-01-21 一种高性能冗余度机械臂重复运动规划方法与装置 Active CN109623827B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910055288.5A CN109623827B (zh) 2019-01-21 2019-01-21 一种高性能冗余度机械臂重复运动规划方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910055288.5A CN109623827B (zh) 2019-01-21 2019-01-21 一种高性能冗余度机械臂重复运动规划方法与装置

Publications (2)

Publication Number Publication Date
CN109623827A CN109623827A (zh) 2019-04-16
CN109623827B true CN109623827B (zh) 2023-04-07

Family

ID=66062294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910055288.5A Active CN109623827B (zh) 2019-01-21 2019-01-21 一种高性能冗余度机械臂重复运动规划方法与装置

Country Status (1)

Country Link
CN (1) CN109623827B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110103225B (zh) * 2019-06-04 2023-04-11 兰州大学 一种数据驱动的机械臂重复运动控制方法与装置
CN111890363B (zh) * 2020-07-27 2022-12-30 四川大学 一种基于快速自适应梯度神经网络算法的机械臂运动规划方法
CN112621761B (zh) * 2020-12-24 2022-06-24 中国科学院重庆绿色智能技术研究院 一种面向通信时滞的机械臂***多级优化协调控制方法
CN112828894B (zh) * 2021-01-21 2022-09-02 中国科学院重庆绿色智能技术研究院 一种冗余度机械臂的位置与力混合控制方法
CN113985738A (zh) * 2021-11-02 2022-01-28 长春工业大学 非凸约束的全向四轮移动机械臂重复运动的梯度神经网络协同控制

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108326844A (zh) * 2017-01-20 2018-07-27 香港理工大学深圳研究院 冗余度机械臂的可操作度优化的运动规划方法和装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8065060B2 (en) * 2006-01-18 2011-11-22 The Board Of Regents Of The University And Community College System On Behalf Of The University Of Nevada Coordinated joint motion control system with position error correction
JP5398588B2 (ja) * 2010-03-01 2014-01-29 本田技研工業株式会社 脚式移動ロボットの歩容生成装置
CN101927495B (zh) * 2010-08-25 2013-04-17 中山大学 一种冗余度机械臂重复运动规划方法
CN106426164B (zh) * 2016-09-27 2019-04-09 华南理工大学 一种冗余度双机械臂的多指标协调运动规划方法
CN106737670B (zh) * 2016-12-15 2019-01-25 华侨大学 一种具有抗噪特性的冗余度机械臂重复运动规划方法
CN106945041B (zh) * 2017-03-27 2019-08-20 华南理工大学 一种冗余度机械臂重复运动规划方法
CN108422424A (zh) * 2018-05-28 2018-08-21 兰州大学 一种具有饱和特性的抗扰动机械臂重复运动规划方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108326844A (zh) * 2017-01-20 2018-07-27 香港理工大学深圳研究院 冗余度机械臂的可操作度优化的运动规划方法和装置

Also Published As

Publication number Publication date
CN109623827A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109623827B (zh) 一种高性能冗余度机械臂重复运动规划方法与装置
CN110103225B (zh) 一种数据驱动的机械臂重复运动控制方法与装置
Zhao et al. System identification of the nonlinear residual errors of an industrial robot using massive measurements
CN108015766B (zh) 一种非线性约束的原对偶神经网络机器人动作规划方法
CN107607918B (zh) 一种基于机器人的柱面近场测量馈源定位及偏焦方法
CN105773623A (zh) 基于预测型间接迭代学习的scara机器人轨迹跟踪控制方法
CN108098777B (zh) 一种冗余度机械臂力矩层重复运动控制方法
CN105598968B (zh) 一种并联机械臂的运动规划与控制方法
CN109129486B (zh) 一种抑制周期噪声的冗余度机械臂重复运动规划方法
CN104760041A (zh) 一种基于突加度的障碍物躲避运动规划方法
CN112417755A (zh) 一种主从式手术机器人轨迹预测控制方法
US20220305646A1 (en) Simulation-in-the-loop Tuning of Robot Parameters for System Modeling and Control
CN109129487B (zh) 在周期噪声下基于泰勒型离散周期节律神经网络的冗余度机械臂重复运动规划方法
CN109623812B (zh) 考虑航天器本体姿态运动的机械臂轨迹规划方法
CN113146600A (zh) 基于运动学迭代学习控制的柔性机器人轨迹规划方法及装置
CN111230882A (zh) 一种水果分拣并联机器人夹持机构的自适应变阻抗控制方法
CN115157262B (zh) 关节受限冗余机械臂的障碍物躲避与优化控制方法、***及机器人
CN110561441B (zh) 一种冗余度机械臂位姿控制的单94lvi迭代算法
CN114055467B (zh) 基于五自由度机器人的空间位姿在线仿真***
CN110434854B (zh) 一种基于数据驱动的冗余度机械臂视觉伺服控制方法与装置
CN114310914A (zh) 多自由度机械臂模糊自适应迭代轨迹跟踪控制方法及***
Lämmle et al. Simulation-based learning of the peg-in-hole process using robot-skills
JPH0683427A (ja) 柔軟マニピュレータの作業点位置制御システム
JP2003340755A (ja) サーボモータ制御用コントローラにおけるゲイン設定法、コントローラの有効性検証法およびロボット制御法
CN109648567B (zh) 一种具有容噪特性的冗余度机械臂规划方法

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