CN104935235A - 电机的电流谐波抑制方法 - Google Patents

电机的电流谐波抑制方法 Download PDF

Info

Publication number
CN104935235A
CN104935235A CN201510349263.8A CN201510349263A CN104935235A CN 104935235 A CN104935235 A CN 104935235A CN 201510349263 A CN201510349263 A CN 201510349263A CN 104935235 A CN104935235 A CN 104935235A
Authority
CN
China
Prior art keywords
current
harmonic
axle
qhe
dhe
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
CN201510349263.8A
Other languages
English (en)
Other versions
CN104935235B (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.)
Shanghai Step Electric Corp
Shanghai Sigriner Step Electric Co Ltd
Original Assignee
Shanghai Step Electric Corp
Shanghai Sigriner Step Electric Co Ltd
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 Shanghai Step Electric Corp, Shanghai Sigriner Step Electric Co Ltd filed Critical Shanghai Step Electric Corp
Priority to CN201510349263.8A priority Critical patent/CN104935235B/zh
Publication of CN104935235A publication Critical patent/CN104935235A/zh
Application granted granted Critical
Publication of CN104935235B publication Critical patent/CN104935235B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开了一种电机的电流谐波抑制方法,包括以下步骤:根据接收的d轴给定电流id *和d轴反馈电流id求差后的信号idh以及q轴给定电流iq *和q轴反馈电流iq求差后的信号iqh,估计出id的谐波电流idhe和iq的谐波电流iqhe;根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco,将d轴补偿电压udco和q轴补偿电压uqco分别加到d轴电流环控制器输出的d轴电压ud和q轴电流环控制器输出的q轴电压uq上。采用上述技术方案后,能够在电机矢量控制的过程中有效地消除电机的电流谐波。

Description

电机的电流谐波抑制方法
技术领域
本发明涉及电机的电流谐波抑制方法。
背景技术
由于电机齿槽效应等本身结构的原因,使得电机的反电势等并不是理想的形状,这会导致电机的电流不为正弦波,而含有谐波。谐波电流会引起转矩脉动,使得电机性能降低。对电流谐波进行抑制,可以提高电机的性能。
以三相异步电机举例来说,电流包含5次、7次、11次、13次、17次等谐波,当谐波频率较低时,电流环控制器可以较好地抑制谐波电流,但当频率升高时,由于电流环控制器带宽的限制,对谐波抑制的作用降低。要抑制谐波,可以对谐波电流进行单独控制。矢量控制由于其优良性能在异步电机控制中得到了广泛应用,矢量控制将三相电流转换到两相旋转坐标系下,实现了转矩和磁链的解耦。三相电流经过旋转坐标变换后,所含电流谐波变成6次、12次、18次等,对这些阶次的谐波进行补偿就可以起到谐波抑制的作用。
发明内容
本发明所要解决的技术问题在于提供一种电机的电流谐波抑制方法,其能够在电机矢量控制的过程中有效地消除电机的电流谐波。
为解决上述技术问题,本发明所采取的技术方案是:
电机的电流谐波抑制方法,其特点在于,包括以下步骤:
根据接收的d轴给定电流id *和d轴反馈电流id求差后的信号idh以及q轴给定电流iq *和q轴反馈电流iq求差后的信号iqh,估计出id的谐波电流idhe和iq的谐波电流iqhe
根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco,将d轴补偿电压udco和q轴补偿电压uqco分别加到d轴电流环控制器输出的d轴电压ud和q轴电流环控制器输出的q轴电压uq上。
本发明至少具有以下技术效果:
1、本发明根据估算出的电流谐波计算出d轴补偿电压udco和q轴补偿电压uqco,并将d轴补偿电压udco和q轴补偿电压uqco分别加到d轴电流环控制器输出的d轴电压ud和q轴电流环控制器输出的q轴电压uq上,实现了电机矢量控制过程中的谐波抑制,消除谐波后的电机相电流更加接近正弦波,从而提高了异步电机的性能;
2、本发明是通过估算出的电流谐波计算出d轴补偿电压udco和q轴补偿电压uqco,可以对采样和控制器的延时进行补偿,从而达到更好的谐波消除效果。
附图说明
图1是现有的电机矢量控制的控制框图。
图2示出了采用了根据本发明一实施例的电流谐波抑制方法的电机矢量控制的控制框图。
图3是根据本发明一实施例的电流谐波抑制模块的原理示意图。
图4是根据本发明一应用实例的d轴反馈电流id及其电流谐波实际值idh的波形示意图。
图5是根据本发明一应用实例的q轴反馈电流iq及其电流谐波实际值iqh的波形示意图。
图6是根据本发明一应用实例的d轴电流谐波实际值idh和估计值idhe的波形示意图。
图7是根据本发明一应用实例的q轴电流谐波实际值iqh和估计值iqhe的波形示意图。
图8是未采用根据本发明一应用实例的谐波抑制方法的电流波形示意图。
图9是采用了根据本发明一应用实例的谐波抑制方法的电流波形示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本发明是在电机矢量控制的基础上实现的。现有矢量控制的控制框图如图1所示。图1中,ω*和ω分别为给定转速和反馈转速,id *、id分别为d轴给定电流和反馈电流,iq *、iq分别为q轴给定电流和反馈电流,ia、ib分别为电机的A相和B相的相电流,iα、iβ为静止两相坐标系下的电流分量,ud、uq分别为d轴电压和q轴电压,uα,uβ为静止两相坐标系下的电压分量,ωe为电机同步角速度,速度控制器和电流环控制器均为PI调节器。
将电流转换到同步旋转坐标系后,d轴、q轴电流中含有6次、12次、18次、24次等谐波,根据傅立叶运算公式,d轴、q轴电流id、iq的表达式可以由下式来表示:
上式中,d0为id直流分量,da6n为id的6n次谐波正弦分量系数,db6n为id的6n次谐波余弦分量系数,q0为iq直流分量,qa6n为iq的6n次谐波正弦分量系数,qb6n为iq的6n次谐波余弦分量系数。
当电流环控制器可以很好地实现电流控制时,给定与反馈的差值为零。可以认为,电流环控制器反馈与给定的差值就是电流的谐波,那么,谐波就可以由下式来表示:
图2示出了采用了根据本发明一实施例的电流谐波抑制方法的电机矢量控制的控制框图。图2中的虚线框表示采用了根据本发明一实施例的电流谐波抑制方法的谐波抑制模块,K为谐波控制器增益,idhe和iqhe分别为估计出的d轴和q轴电流的谐波电流。
基于前述的分析,并结合图2所示,根据本发明一实施例的电机的电流谐波抑制方法,包括以下步骤:
步骤S11、根据接收的d轴给定电流id *和d轴反馈电流id求差后的信号idh以及q轴给定电流iq *和q轴反馈电流iq求差后的信号iqh,估计出id的谐波电流idhe和iq的谐波电流iqhe
上述的根据idh和iqh估计出id的谐波电流idhe和iq的谐波电流iqhe具体包括以下步骤:
a、选择要消除的谐波次数,然后根据以下公式计算出要消除的各次谐波的系数:
其中,da6n为id的6n次谐波正弦分量系数,db6n为id的6n次谐波余弦分量系数,qa6n为iq的6n次谐波正弦分量系数,qb6n为iq的6n次谐波余弦分量系数,ωe为电机同步角速度,T为电流周期,n=1,2,3......。
为使计算变得简便,在一种优选的实施方式中,可以采用低通滤波器来实现对上述各次谐波系数的计算。具体地,id的谐波电流idhe的各次谐波的系数通过以下方式计算获得:将idh*sin(6nωet)和idh*cos(6nωet)分别输入到一个低通滤波器,并将各低通滤波器的输出乘以2分别得到da6n和db6n;iq的谐波电流iqhe的各次谐波的系数通过以下方式计算获得:将iqh*sin(6nωet)和iqh*cos(6nωet)分别输入到一个低通滤波器,并将各低通滤波器的输出乘以2分别得到qa6n和qb6n
b、根据以下公式估计出id和iq的谐波电流idhe和iqhe
如要消除5次和7次,就相当于消除d轴和q轴电流中的6次谐波,只要计算6次谐波的系数,如果要消除所有谐波,就要计算6次、12次、18次、24次……,结合实际需要和控制器的运算能力进行选择。计算出各次谐波系数后,就可以求出谐波估计值idhe和iqhe
S12、根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco,将d轴补偿电压udco和q轴补偿电压uqco分别加到d轴电流环控制器输出的d轴电压ud和q轴电流环控制器输出的q轴电压uq上,从而实现谐波抑制。
请结合图3所示。在图3所示的实施例中,谐波控制器采用了比例调节器。前述的根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco具体包括以下步骤:
将0与idhe求差后的信号输入一比例调节器,该比例调节器计算出d轴补偿电压udco
将0与iqhe求差后的信号输入另一比例调节器,该另一比例调节器计算出q轴补偿电压uqco
为使本发明的技术方案更加清楚,以下结合一具体的应用实例进行说明。
本应用实例中使用了电压型变频器对异步电机进行控制来实现本发明的电流谐波抑制方法。变频器额定功率为11kw,额定电流为27A,异步电机为额定功率为3.7kw,额定电流为9.5A,系统载波频率为8k,未进行死区补偿。
首先,根据接收的d轴给定电流id *和d轴反馈电流id求差后的信号idh以及q轴给定电流iq *和q轴反馈电流iq求差后的信号iqh,估计出id的谐波电流idhe和iq的谐波电流iqhe
将电流环控制器反馈与给定的差值作为电流的实际谐波,如图4和图5所示,其中图4为id及其谐波idh,图5为iq及其谐波iqh
目标为消除d轴、q轴电流的6次、12次、18次谐波,只需将n设为1、2和3,先分别求出6次、12次、18次谐波的谐波系数da6、db6、da12、db12、da18、db18、qa6、qb6、qa12、qb12、qa18和qb18,得到id、iq电流谐波的表达式,根据表达式可以估计出电流谐波的大小。图6示出了id的电流谐波实际值idh与电流谐波估计值idhe,图7示出了iq的电流谐波的实际值iqh与电流谐波估计值iqhe
然后,计算补偿电压,将补偿电压加到d、q轴电流环控制器输出上,实现谐波抑制。
这样,就可以消除相电流中的5次、7次、11次、13次、17次和19次谐波。其中,谐波控制器增益设为2.5。请参阅图8和图9,图8为未进行谐波抑制的电流波形,图9为进行谐波抑制后的电流波形,通过比较可见,进行谐波抑制后的电机相电流更加接近正弦波,iq、id更加接近直流,这说明本方法可以很好地实现谐波抑制。
本发明不仅适用于异步电机,也适用于永磁同步电机。

Claims (5)

1.电机的电流谐波抑制方法,其特征在于,包括以下步骤:
根据接收的d轴给定电流id *和d轴反馈电流id求差后的信号idh以及q轴给定电流iq *和q轴反馈电流iq求差后的信号iqh,估计出id的谐波电流idhe和iq的谐波电流iqhe
根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco,将d轴补偿电压udco和q轴补偿电压uqco分别加到d轴电流环控制器输出的d轴电压ud和q轴电流环控制器输出的q轴电压uq上。
2.根据权利要求1所述的电机的电流谐波抑制方法,其特征在于,所述的根据idh和iqh估计出id的谐波电流idhe和iq的谐波电流iqhe包括以下步骤:
选择要消除的谐波次数,然后根据以下公式计算出要消除的各次谐波的系数:
d a 6 n = 2 T ∫ 0 T i d h sin ( 6 nω e t ) d t d b 6 n = 2 T ∫ 0 T i d h cos ( 6 nω e t ) d t q a 6 n = 2 T ∫ 0 T i q h sin ( 6 nω e t ) d t q b 6 n = 2 T ∫ 0 T i q h cos ( 6 nω e t ) d t
其中,da6n为id的6n次谐波正弦分量系数,db6n为id的6n次谐波余弦分量系数,qa6n为iq的6n次谐波正弦分量系数,qb6n为iq的6n次谐波余弦分量系数,ωe为电机同步角速度,T为电流周期,n=1,2,3......;
根据以下公式估计出id和iq的谐波电流idhe和iqhe
i d h e = Σ n = 1 ∞ da 6 n s i n ( 6 nω e t ) + d b 6 n c o s ( 6 n ω e t ) i q h e = Σ n = 1 ∞ qa 6 n s i n ( 6 nω e t ) + q b 6 n c o s ( 6 n ω e t )
3.根据权利要求2所述的电机的电流谐波抑制方法,其特征在于,id的谐波电流idhe的各次谐波的系数通过以下方式计算获得:
将idh*sin(6nωet)和idh*cos(6nωet)分别输入到一个低通滤波器,并将各低通滤波器的输出乘以2分别得到da6n和db6n
iq的谐波电流iqhe的各次谐波的系数通过以下方式计算获得:
将iqh*sin(6nωet)和iqh*cos(6nωet)分别输入到一个低通滤波器,并将各低通滤波器的输出乘以2分别得到qa6n和qb6n
4.根据权利要求1所述的电机的电流谐波抑制方法,其特征在于,根据估计出的谐波电流idhe和iqhe计算d轴补偿电压udco和q轴补偿电压uqco具体包括以下步骤:
将0与idhe求差后的信号输入一比例调节器,该比例调节器计算出d轴补偿电压udco
将0与iqhe求差后的信号输入另一比例调节器,该另一比例调节器计算出q轴补偿电压uqco
5.根据权利要求1所述的电机的电流谐波抑制方法,其特征在于,所述的电机为异步电机或永磁同步电机。
CN201510349263.8A 2015-06-19 2015-06-19 电机的电流谐波抑制方法 Active CN104935235B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510349263.8A CN104935235B (zh) 2015-06-19 2015-06-19 电机的电流谐波抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510349263.8A CN104935235B (zh) 2015-06-19 2015-06-19 电机的电流谐波抑制方法

Publications (2)

Publication Number Publication Date
CN104935235A true CN104935235A (zh) 2015-09-23
CN104935235B CN104935235B (zh) 2018-05-08

Family

ID=54122244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510349263.8A Active CN104935235B (zh) 2015-06-19 2015-06-19 电机的电流谐波抑制方法

Country Status (1)

Country Link
CN (1) CN104935235B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105915136A (zh) * 2016-05-29 2016-08-31 南京理工大学 基于模糊神经网络的电机电流谐波抑制控制系统及方法
CN106788127A (zh) * 2017-01-17 2017-05-31 浙江大学 基于二维查表与插值法的逆变器非线性谐波补偿方法
CN109660116A (zh) * 2017-10-12 2019-04-19 美的集团股份有限公司 三相变流器及其控制方法、装置
CN110574280A (zh) * 2017-03-07 2019-12-13 罗伯特·博世有限公司 用于电机的调节设备、电驱动系统和用于调节电机的方法
CN111371363A (zh) * 2020-04-01 2020-07-03 上海新时达电气股份有限公司 基于信号注入的电机最大转矩电流比控制方法及其装置
CN111464085A (zh) * 2020-04-29 2020-07-28 华南理工大学 一种基于阶次提取的电机电流谐波及转矩脉动抑制方法
CN111551781A (zh) * 2020-05-18 2020-08-18 湖南大学 等效畸变电压直接检测方法及逆变器非线性因素补偿方法
CN111835258A (zh) * 2020-07-27 2020-10-27 广东希塔变频技术有限公司 电机驱动控制方法、控制装置和驻车空调控制器
CN112134495A (zh) * 2020-08-19 2020-12-25 湖州慧微电子科技有限公司 永磁同步电机电流的开环谐波补偿方法
CN112737447A (zh) * 2020-12-28 2021-04-30 广东美芝制冷设备有限公司 压缩机组件、制冷设备和可读存储介质
CN112787558A (zh) * 2020-12-28 2021-05-11 广东美芝制冷设备有限公司 压缩机组件的控制方法、压缩机组件和制冷设备
CN113507250A (zh) * 2021-07-14 2021-10-15 浙江大学先进电气装备创新中心 一种内置式永磁同步电机电流谐波抑制方法
CN114785220A (zh) * 2022-03-28 2022-07-22 吉林大学 一种基于车载永磁同步电机转矩波动抑制控制系统的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102201770A (zh) * 2011-05-30 2011-09-28 重庆大学 一种注入谐波电压抑制永磁同步电机谐波电流的控制方法
CN103718452A (zh) * 2011-06-03 2014-04-09 Trw有限公司 具有电压谐波整形的电动机控制
WO2014077596A1 (ko) * 2012-11-14 2014-05-22 포스코에너지 주식회사 인버터의 리플 및 옵셋 보상 장치 및 그 방법
CN103973191A (zh) * 2014-04-09 2014-08-06 东南大学 一种九相磁通切换永磁电机的缺相容错控制方法
CN104579042A (zh) * 2013-10-22 2015-04-29 广东美的制冷设备有限公司 永磁同步电机的控制系统及其力矩波动抑制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102201770A (zh) * 2011-05-30 2011-09-28 重庆大学 一种注入谐波电压抑制永磁同步电机谐波电流的控制方法
CN103718452A (zh) * 2011-06-03 2014-04-09 Trw有限公司 具有电压谐波整形的电动机控制
WO2014077596A1 (ko) * 2012-11-14 2014-05-22 포스코에너지 주식회사 인버터의 리플 및 옵셋 보상 장치 및 그 방법
CN104579042A (zh) * 2013-10-22 2015-04-29 广东美的制冷设备有限公司 永磁同步电机的控制系统及其力矩波动抑制方法
CN103973191A (zh) * 2014-04-09 2014-08-06 东南大学 一种九相磁通切换永磁电机的缺相容错控制方法

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105915136A (zh) * 2016-05-29 2016-08-31 南京理工大学 基于模糊神经网络的电机电流谐波抑制控制系统及方法
CN105915136B (zh) * 2016-05-29 2018-10-12 南京理工大学 基于模糊神经网络的电机电流谐波抑制控制系统及方法
CN106788127A (zh) * 2017-01-17 2017-05-31 浙江大学 基于二维查表与插值法的逆变器非线性谐波补偿方法
CN106788127B (zh) * 2017-01-17 2018-11-20 浙江大学 基于二维查表与插值法的逆变器非线性谐波补偿方法
CN110574280A (zh) * 2017-03-07 2019-12-13 罗伯特·博世有限公司 用于电机的调节设备、电驱动系统和用于调节电机的方法
CN110574280B (zh) * 2017-03-07 2023-02-21 罗伯特·博世有限公司 用于电机的调节设备、电驱动系统和用于调节电机的方法
CN109660116A (zh) * 2017-10-12 2019-04-19 美的集团股份有限公司 三相变流器及其控制方法、装置
CN111371363A (zh) * 2020-04-01 2020-07-03 上海新时达电气股份有限公司 基于信号注入的电机最大转矩电流比控制方法及其装置
CN111371363B (zh) * 2020-04-01 2023-08-15 上海新时达电气股份有限公司 基于信号注入的电机最大转矩电流比控制方法及其装置
CN111464085B (zh) * 2020-04-29 2022-03-29 华南理工大学 一种基于阶次提取的电机电流谐波及转矩脉动抑制方法
CN111464085A (zh) * 2020-04-29 2020-07-28 华南理工大学 一种基于阶次提取的电机电流谐波及转矩脉动抑制方法
CN111551781B (zh) * 2020-05-18 2021-06-11 湖南大学 等效畸变电压直接检测方法及逆变器非线性因素补偿方法
CN111551781A (zh) * 2020-05-18 2020-08-18 湖南大学 等效畸变电压直接检测方法及逆变器非线性因素补偿方法
CN111835258A (zh) * 2020-07-27 2020-10-27 广东希塔变频技术有限公司 电机驱动控制方法、控制装置和驻车空调控制器
CN112134495A (zh) * 2020-08-19 2020-12-25 湖州慧微电子科技有限公司 永磁同步电机电流的开环谐波补偿方法
CN112737447A (zh) * 2020-12-28 2021-04-30 广东美芝制冷设备有限公司 压缩机组件、制冷设备和可读存储介质
CN112787558A (zh) * 2020-12-28 2021-05-11 广东美芝制冷设备有限公司 压缩机组件的控制方法、压缩机组件和制冷设备
CN112737447B (zh) * 2020-12-28 2023-01-31 广东美芝制冷设备有限公司 压缩机组件、制冷设备和可读存储介质
CN113507250A (zh) * 2021-07-14 2021-10-15 浙江大学先进电气装备创新中心 一种内置式永磁同步电机电流谐波抑制方法
CN114785220A (zh) * 2022-03-28 2022-07-22 吉林大学 一种基于车载永磁同步电机转矩波动抑制控制系统的方法

Also Published As

Publication number Publication date
CN104935235B (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN104935235A (zh) 电机的电流谐波抑制方法
Andersson et al. Motion sensorless IPMSM control using linear moving horizon estimation with Luenberger observer state feedback
CN110429886B (zh) 一种永磁同步电机低速域转子位置辨识方法
EP2779414B1 (en) Motor control system having bandwidth compensation
CN103997272B (zh) 永磁同步电机的负载扰动补偿装置及方法
CN110518852B (zh) 基于谐波注入的永磁同步电机电流谐波抑制方法
CN105680756B (zh) 一种用于双三相异步电机的控制方法以及装置
US9048778B2 (en) Rotor phase/speed estimating device for an AC motor
US11316451B2 (en) Method for determining current-dependent and/or rotational angle position-dependent characteristic variables of an electrical machine, and frequency converter
US9998043B2 (en) Rotary machine controller
CN105024615A (zh) 一种永磁同步电机低速无传感器控制方法及装置
CN116232154A (zh) 基于复矢量离散滤波器的电机参数估计及控制方法和系统
EP2747273B1 (en) Method and arrangement for torque estimation of a synchronous machine
CN109495027B (zh) 带lc滤波器的永磁同步电机磁极位置估计方法及系统
CN108494307B (zh) 一种永磁同步电机的逆变器非线性效应补偿方法
CN111106775A (zh) 永磁同步电机角度滤波方法及系统
Pravica et al. Interior permanent magnet wind generator torque estimation considering low-pass filter phase shift compensation
Jin et al. Torque ripple reduction based on flux linkage harmonics observer for an interior PM synchronous motor including back EMF harmonics
CN112653360B (zh) 一种高速永磁同步电机无位置传感器控制方法
Zhang et al. Research on a sensorless SVM-DTC strategy for induction motors based on modified stator model
KR20180080730A (ko) 모터 제어 시스템 및 방법
Yuan et al. Cross-Cancellation Sliding Mode Observer-Based Sensorless Control of PMSLM
Sahoo et al. Implementation and simulation of direct torque control scheme with the use of FPGA circuit
CN108923701B (zh) 一种无辅助信号注入的永磁同步电机转子位置权重观测器
CN115001341A (zh) 考虑死区效应的无传感器控制方法、装置以及硬件设备

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant