CN109581967A - 一种基于机电耦合作用的高速进给系统速度避振规划方法 - Google Patents

一种基于机电耦合作用的高速进给系统速度避振规划方法 Download PDF

Info

Publication number
CN109581967A
CN109581967A CN201811643600.4A CN201811643600A CN109581967A CN 109581967 A CN109581967 A CN 109581967A CN 201811643600 A CN201811643600 A CN 201811643600A CN 109581967 A CN109581967 A CN 109581967A
Authority
CN
China
Prior art keywords
speed
servo
mechanical
vibration
drive
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
CN201811643600.4A
Other languages
English (en)
Other versions
CN109581967B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201811643600.4A priority Critical patent/CN109581967B/zh
Publication of CN109581967A publication Critical patent/CN109581967A/zh
Application granted granted Critical
Publication of CN109581967B publication Critical patent/CN109581967B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/406Numerical 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 monitoring or safety
    • G05B19/4062Monitoring servoloop, e.g. overload of servomotor, loss of feedback or reference
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/33Director till display
    • G05B2219/33327Self diagnostic of control system, servo system
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/34Director, elements to supervisory
    • G05B2219/34013Servocontroller
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37634By measuring vibration

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Apparatuses For Generation Of Mechanical Vibrations (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供的一种基于机电耦合作用的高速进给系统速度避振规划方法,包括以下步骤:步骤1,计算伺服驱动输出的谐波特性,确定伺服驱动输出得到演变频率区间;步骤2,确定驱动机械系统的主要振动模态分布与振型;步骤3,建立不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系;步骤4,提取高速进给系统共振速度敏感参数点集;步骤5,将步骤4中得到的进给系统共振速度敏感参数点从运动速度指令中进行剔除,得到机电耦合作用关系的高速进给系统进给速度,完成速度避振规划;本发明能够有效的避免由于伺服驱动谐波激起机械共振而导致的精度恶化的问题,对于保证零件的加工精度和表面质量,实现复杂零件的智能制造具有重要的价值和意义。

Description

一种基于机电耦合作用的高速进给系统速度避振规划方法
技术领域
本发明属于高速高精数控机床智能制造领域,具体涉及一种考虑机电耦合作用的高速进给系统速度避振规划方法,适用于高速高精数控机床等场合。
背景技术
进给速度是数控机床切削加工中的重要参数,其大小对于零件加工精度和加工效率具有重要的影响。目前进给速度主要根据零件的加工精度以及表面粗糙度要求和刀具、工件的材料性质进行选取。在加工不同的零件时,会根据不同加工特点进行调整。最大进给速度受机床刚度以及进给系统的整体性能限制。
然而,在数控机床实际加工过程中,受驱动电路和电机结构非线性的影响,伺服驱动输出并非理想的恒定输出,而是夹杂了众多的谐波干扰成分,而且这些电机输出干扰的频率会随着进给速度的不同而发生演变。数控机床的机械系统受各个动静结合部的影响,具有多阶振动模态,而且在不同的运动位置,机械各阶模态也会发生变化。在运动过程中伺服驱动和机械系统之间存在着紧密的耦合作用关系。尤其面向复杂零件的智能制造过程,运动参数和机械位置在不断的发生变化,在某一速度下,伺服驱动输出的某一阶谐波成分频率接近机械系统某一阶固有频率,引起机械共振,会严重恶化系统的运动精度,影响零件最终的加工精度。目前的进给速度选择方法没有考虑伺服驱动输出的谐波特性和机械系统的振动特性,虽然在一定条件下能够满足使用需求,但是会出现在某些进给速度下,系统运动精度和零件加工精度急剧恶化的现象。尤其对于智能制造所需要的高速高精数控机床,伺服输出谐波频宽范围大,机械系统振动模态众多,如果在速度规划过程中不考虑伺服驱动和机械系统之间的机电耦合作用,很难满足实际的工程需求。
发明内容
本发明的目的在于提供一种基于机电耦合作用的高速进给系统速度避振规划方法,解决了现有的进给速度选择方法没有考虑伺服驱动输出的谐波特性和机械系统的振动特性,存在系统运动精度和零件加工精度急剧恶化的现象。
为了达到上述目的,本发明采用的技术方案是:
本发明提供的一种基于机电耦合作用的高速进给系统速度避振规划方法,包括以下步骤:
步骤1,根据进给系统的驱动方式和特点、驱动电路和电机结构非线性,计算伺服驱动输出的谐波特性,接着建立伺服驱动输出各项谐波成分与进给速度之间的演变关系,得到伺服驱动输出的演变频率区间;
步骤2,根据步骤1得到的伺服驱动输出的演变频率区间,得到驱动机械系统的主要振动模态分布与振型;
步骤3,根据步骤1得到的伺服驱动谐波特性和步骤2得到的机械振动模态分布,建立不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系;
步骤4,根据步骤3的耦合作用关系,得到不同速度下伺服驱动谐波与机械振动模态之间的相关性图,从中提取高速进给系统共振速度敏感参数点集;
步骤5,将步骤4中得到的进给系统共振速度敏感参数点从运动速度指令中进行剔除,得到机电耦合作用关系的高速进给系统进给速度,完成速度避振规划。
优选地,步骤1中,伺服驱动输出的谐波特性的表达式:
其中,F是伺服驱动输出,v是进给速度,为偏导运算符,t为时间,La、Lb、Lc分别为三相线圈的自感,Mab、Mac、Mbc分别为三相线圈间的互感,ia、ib、ic分别为三相电流,N为线圈匝数,l为线圈宽度,τ为电机极距,xa为电机气隙在进给方向上的不同位置,x为积分变量,Ba、Bb、Bc分别为永磁体产生的气隙三相磁场强度,F0是电机理想的输出,为伺服驱动输出的谐波,Ai、ωi分别为伺服驱动输出各项谐波的幅值、频率和相位。
优选地,步骤1中,伺服驱动输出各项谐波成分与进给速度之间的演变关系:
再根据实际工艺所需求的进给速度范围,得到伺服驱动输出的演变区间。
优选地,步骤1中,进给系统的驱动方式与特点是指高速进给系统采用旋转电机加滚珠丝杠驱动方式或直线电机的零传动驱动方式的驱动方式。
优选地,步骤1中,驱动电路和电机结构非线性分别是指驱动电路三相输出电流中的谐波成分和电机永磁体磁场中的谐波成分。
优选地,步骤2中,驱动机械系统的主要振动模态分布与振型具体是:
(A-λE)x=0
其中,A为机械模态矩阵,λ为机械固有频率,E为单位矩阵,x为振型向量。
优选地,步骤3中,不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系:
其中,Fvi为不同速度下对应的伺服输出谐波,M为机械振动模态。
优选地,步骤4中,提取高速进给系统共振速度敏感参数点集的具体方法是:将相关性图上伺服驱动谐波与机械振动模态之间的交点作为敏感点,进而提取到高速进给系统共振速度敏感参数点集。
与现有技术相比,本发明的有益效果是:
本发明提供的一种基于机电耦合作用的高速进给系统速度避振规划方法,综合考虑伺服驱动谐波特性以及其随速度的演变规律、数控机床机械振动特性以及二者之间的耦合作用关系,在传统进给速度选择的基础上,剔除共振速度敏感点,能够有效的避免由于伺服驱动谐波激起机械共振而导致的精度恶化的问题,对于保证零件的加工精度和表面质量,实现复杂零件的智能制造具有重要的价值和意义。
附图说明
附图1考虑机电耦合的速度避振规划流程图;
附图2伺服驱动输出谐波随进给速度的演变规律;
附图3机械系统主要振动模态及振型分布;
附图4伺服驱动和机械系统间的机电耦合作用关系。
具体实施方式
下面结合附图,对本发明进一步详细说明。
如图1所示,本发明提供的一种基于机电耦合作用的高速进给系统速度避振规划方法,包括以下步骤:
步骤一,根据进给系统的驱动方式和特点、驱动电路和电机结构非线性,计算伺服驱动输出的谐波特性,即
其中,F是伺服驱动输出,v是进给速度,为偏导运算符,t为时间,La、Lb、Lc分别为三相线圈的自感,Mab、Mac、Mbc分别为三相线圈间的互感,ia、ib、ic分别为三相电流,N为线圈匝数,l为线圈宽度,τ为电机极距,xa为电机气隙在进给方向上的不同位置,x为积分变量,Ba、Bb、Bc分别为永磁体产生的气隙三相磁场强度,F0是电机理想的输出,为伺服驱动输出的谐波,Ai、ωi分别为伺服驱动输出各项谐波的幅值、频率和相位。
进给系统的驱动方式与特点是指高速进给系统采用旋转电机加滚珠丝杠驱动方式或直线电机的零传动驱动方式的驱动方式;
驱动电路和电机结构非线性是指驱动电路三相输出电流中的谐波成分和电机永磁体磁场中的谐波成分。三相线圈的自感和互感可通过电机厂商或实际参数辨识得到,三相电流和气隙三相磁场强度可以通过传统理论方法计算得到也可以通过测试得到。
建立伺服驱动输出各项谐波成分随进给速度的演变关系,即
根据实际工艺所需求的进给速度范围,得到伺服驱动输出的演变区间。
步骤二,根据步骤一得到的伺服驱动输出的演变频率区间,得到驱动机械系统的主要振动模态分布与振型,即
(A-λE)x=0
其中,A为机械模态矩阵,λ为机械固有频率,E为单位矩阵,x为振型向量;
步骤三,根据步骤一得到的伺服驱动谐波特性和步骤二得到的机械振动模态分布,建立不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系:
其中,Fvi为不同速度下对应的伺服输出谐波,M为机械振动模态。
步骤四,根据步骤三的耦合作用关系,得到不同速度下伺服驱动谐波与机械振动模态之间的相关性图,将该相关性图上伺服驱动谐波与机械振动模态之间的交点作为敏感点进行提取,得到高速进给系统共振速度敏感参数点集,即
v=[v1v2…vi]
将敏感速度集从运动速度指令中进行剔除,得到考虑机电耦合作用关系的高速进给系统进给速度,完成速度避振规划。
选择某台配有直线电机进给系统的高速五轴立式加工中心中的x轴作为实施案例,具体步骤如下:
步骤一,根据进给系统的驱动方式和特点(案例中采用的是直线电机直接驱动方式),以及驱动电路和电机结构的非线性(主要考虑伺服驱动电流谐波以及电机端部效应和齿槽效应造成气隙磁场谐波),计算伺服驱动输出的谐波特性,即
其中,F是伺服驱动输出,v是进给速度,为偏导运算符,t为时间,La、Lb、Lc分别为三相线圈的自感,Mab、Mac、Mbc分别为三相线圈间的互感,ia、ib、ic分别为三相电流,N为线圈匝数,l为线圈宽度,τ为电机极距,xa为电机气隙在进给方向上的不同位置,x为积分变量,Ba、Bb、Bc分别为永磁体产生的气隙三相磁场强度,F0是电机理想的输出,为伺服驱动输出的谐波,Ai、ωi分别为伺服驱动输出各项谐波的幅值、频率和相位。
建立伺服驱动输出各项谐波成分随进给速度的演变关系,如附图2所示,即
根据实际工艺所需求的进给速度范围,得到伺服驱动输出的演变区间为0-100Hz。
步骤二,根据步骤一得到的伺服驱动输出演变主要频率区间,得到驱动机械系统的主要振动模态分布与振型,即
(A-λE)x=0
其中,A为机械模态矩阵,λ为机械固有频率,E为单位矩阵,x为振型向量;
机械系统在0-100Hz范围内主要有6阶模态,如附图3所示,由附图3提取机械前6阶模态分布,结果如表1所示。
表1机械系统前6阶模态分布表
步骤三,根据步骤一得到的伺服驱动谐波特性和步骤二得到的机械振动模态分布,建立不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系:
其中,Fvi为不同速度下对应的伺服输出谐波,M为机械振动模态。
不同速度下伺服驱动输出与机械动态性耦合关系如附图4所示。
步骤四,根据附图4所得分析结果,得到进给系统共振速度敏感参数点集,即
v(m/min)=[1931354658638690]m/min
根据传统方法初步选择满足要求的进给速度,将敏感速度集从运动速度指令中进行剔除,得到考虑机电耦合作用关系的高速进给系统进给速度,完成速度避振规划。

Claims (8)

1.一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,包括以下步骤:
步骤1,根据进给系统的驱动方式和特点、驱动电路和电机结构非线性,计算伺服驱动输出的谐波特性,接着建立伺服驱动输出各项谐波成分与进给速度之间的演变关系,得到伺服驱动输出的演变频率区间;
步骤2,根据步骤1得到的伺服驱动输出的演变频率区间,得到驱动机械系统的主要振动模态分布与振型;
步骤3,根据步骤1得到的伺服驱动谐波特性和步骤2得到的机械振动模态分布,建立不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系;
步骤4,根据步骤3的耦合作用关系,得到不同速度下伺服驱动谐波与机械振动模态之间的相关性图,从中提取高速进给系统共振速度敏感参数点集;
步骤5,将步骤4中得到的进给系统共振速度敏感参数点从运动速度指令中进行剔除,得到机电耦合作用关系的高速进给系统进给速度,完成速度避振规划。
2.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤1中,伺服驱动输出的谐波特性的表达式:
其中,F是伺服驱动输出,v是进给速度,为偏导运算符,t为时间,La、Lb、Lc分别为三相线圈的自感,Mab、Mac、Mbc分别为三相线圈间的互感,ia、ib、ic分别为三相电流,N为线圈匝数,l为线圈宽度,τ为电机极距,xa为电机气隙在进给方向上的不同位置,x为积分变量,Ba、Bb、Bc分别为永磁体产生的气隙三相磁场强度,F0是电机理想的输出,为伺服驱动输出的谐波,Ai、ωi分别为伺服驱动输出各项谐波的幅值、频率和相位。
3.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤1中,伺服驱动输出各项谐波成分与进给速度之间的演变关系:
再根据实际工艺所需求的进给速度范围,得到伺服驱动输出的演变区间。
4.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤1中,进给系统的驱动方式与特点是指高速进给系统采用旋转电机加滚珠丝杠驱动方式或直线电机的零传动驱动方式的驱动方式。
5.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤1中,驱动电路和电机结构非线性分别是指驱动电路三相输出电流中的谐波成分和电机永磁体磁场中的谐波成分。
6.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤2中,驱动机械系统的主要振动模态分布与振型具体是:
(A-λE)x=0
其中,A为机械模态矩阵,λ为机械固有频率,E为单位矩阵,x为振型向量。
7.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤3中,不同进给速度下伺服驱动系统和机械系统之间的机电耦合作用关系:
其中,Fvi为不同速度下对应的伺服输出谐波,M为机械振动模态。
8.根据权利要求1所述的一种基于机电耦合作用的高速进给系统速度避振规划方法,其特征在于,步骤4中,提取高速进给系统共振速度敏感参数点集的具体方法是:将相关性图上伺服驱动谐波与机械振动模态之间的交点作为敏感点,进而提取到高速进给系统共振速度敏感参数点集。
CN201811643600.4A 2018-12-29 2018-12-29 一种基于机电耦合作用的高速进给系统速度避振规划方法 Expired - Fee Related CN109581967B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811643600.4A CN109581967B (zh) 2018-12-29 2018-12-29 一种基于机电耦合作用的高速进给系统速度避振规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811643600.4A CN109581967B (zh) 2018-12-29 2018-12-29 一种基于机电耦合作用的高速进给系统速度避振规划方法

Publications (2)

Publication Number Publication Date
CN109581967A true CN109581967A (zh) 2019-04-05
CN109581967B CN109581967B (zh) 2020-03-17

Family

ID=65932779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811643600.4A Expired - Fee Related CN109581967B (zh) 2018-12-29 2018-12-29 一种基于机电耦合作用的高速进给系统速度避振规划方法

Country Status (1)

Country Link
CN (1) CN109581967B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112099424A (zh) * 2019-06-17 2020-12-18 施耐德电器工业公司 用于监控机器的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1905982A (zh) * 2004-10-21 2007-01-31 松下电器产业株式会社 焊接系统及自耗电极式焊接方法
US7733051B2 (en) * 2006-11-21 2010-06-08 Asm Assembly Automation Ltd. Model-based active electronic damping for stepper motors
CN108021039A (zh) * 2017-12-23 2018-05-11 西安交通大学 一种直线电机进给系统机电集成建模方法
CN108107734A (zh) * 2017-12-23 2018-06-01 西安交通大学 一种永磁同步直线电机进给系统机电耦合建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1905982A (zh) * 2004-10-21 2007-01-31 松下电器产业株式会社 焊接系统及自耗电极式焊接方法
US7733051B2 (en) * 2006-11-21 2010-06-08 Asm Assembly Automation Ltd. Model-based active electronic damping for stepper motors
CN108021039A (zh) * 2017-12-23 2018-05-11 西安交通大学 一种直线电机进给系统机电集成建模方法
CN108107734A (zh) * 2017-12-23 2018-06-01 西安交通大学 一种永磁同步直线电机进给系统机电耦合建模方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
卢秉恒: "高速高加速度下的进给系统机电耦合", 《机械工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112099424A (zh) * 2019-06-17 2020-12-18 施耐德电器工业公司 用于监控机器的方法
CN112099424B (zh) * 2019-06-17 2024-05-07 施耐德电器工业公司 用于监控机器的方法

Also Published As

Publication number Publication date
CN109581967B (zh) 2020-03-17

Similar Documents

Publication Publication Date Title
CN101692118A (zh) 直线电机的测试装置及测试方法
CN108107734B (zh) 一种永磁同步直线电机进给系统机电耦合建模方法
CN107017817B (zh) 一种高速内嵌式永磁同步电机电流解耦控制方法
CN108111088B (zh) 一种考虑气隙波动的永磁同步直线电机推力精确预测方法
CN104898568A (zh) 基于刚度辨识的数控机床进给系统控制参数优化方法
Yang et al. Analysis on the multi-dimensional spectrum of the thrust force for the linear motor feed drive system in machine tools
CN109581967A (zh) 一种基于机电耦合作用的高速进给系统速度避振规划方法
Bu et al. Rotor position tracking control for low speed operation of direct-drive PMSM servo system
Zhan et al. A review on parameters identification methods for asynchronous motor
CN103869245A (zh) 基于磁共能修正的开关磁阻电机磁链曲线测试方法及装置
Zhang et al. Influence of nonlinearity of servo system electrical characteristics on motion smoothness of precision CNC machine tools
CN106919143B (zh) 数值控制装置
Seok Frequency-spectrum-based antiwindup compensator for PI-controlled systems
CN107544435A (zh) 一种数控系统中的珩磨往复运动控制方法
CN110516300A (zh) 一种基于有限元分析的直线永磁电机涡流损耗计算方法
Zhang et al. Numerical iron loss calculation of a new axial flux machine with segmented-armature-torus topology
CN109709891B (zh) 一种直驱高速进给系统伺服参数多目标优化方法
CN103365241B (zh) R型变压器铁芯开料曲线动态补偿控制系统及方法
CN202388313U (zh) 镗床的控制装置
林献坤 et al. Servo Parameter Tuning Method for Direct Feed Drive Based on Dynamic Stiffness Evaluation
Li et al. Research on quartic polynomial velocity planning algorithm based on filtering
Wang et al. Experimental research on the linear motor micro-feed device with high-frequency response, long travel and high accuracy
CN111123834B (zh) 一种直线电机进给系统机电耦合强度评价方法
Li et al. Dynamic balance optimization of the cutting head for flexible materials based on solidworks
Liang et al. An investigation of spectral characteristic of displacement fluctuation and its effect on surface quality in precision machining

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200317

CF01 Termination of patent right due to non-payment of annual fee