CN108718165B - 一种基于误差补偿的感应电机零频稳定控制方法 - Google Patents

一种基于误差补偿的感应电机零频稳定控制方法 Download PDF

Info

Publication number
CN108718165B
CN108718165B CN201810614110.5A CN201810614110A CN108718165B CN 108718165 B CN108718165 B CN 108718165B CN 201810614110 A CN201810614110 A CN 201810614110A CN 108718165 B CN108718165 B CN 108718165B
Authority
CN
China
Prior art keywords
rotor
axis
motor
flux linkage
observer
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
CN201810614110.5A
Other languages
English (en)
Other versions
CN108718165A (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.)
Harbin Institute of Technology
Huazhong University of Science and Technology
Original Assignee
Harbin Institute of Technology
Huazhong University of Science and 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 Harbin Institute of Technology, Huazhong University of Science and Technology filed Critical Harbin Institute of Technology
Priority to CN201810614110.5A priority Critical patent/CN108718165B/zh
Publication of CN108718165A publication Critical patent/CN108718165A/zh
Application granted granted Critical
Publication of CN108718165B publication Critical patent/CN108718165B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/13Observer control, e.g. using Luenberger observers or Kalman filters
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/24Vector control not involving the use of rotor position or rotor speed sensors
    • H02P21/26Rotor flux based control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/09Motor speed determination based on the current and/or voltage without using a tachogenerator or a physical encoder
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/01Asynchronous machines

Abstract

本发明公开了一种基于误差补偿的感应电机零频稳定控制方法,该方法基于现有的感应电机无速度传感器驱动系统,通过在αβ坐标系下电机的定子电压指令输入值
Figure DDA0001696276370000011
和转子电压补偿模块的输入
Figure DDA0001696276370000012
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,或者,通过在dq坐标系下电机的定子电压指令输入值usd、usq和转子电压补偿模块的输入
Figure DDA0001696276370000013
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,实现在保证感应电机无速度传感器驱动系统零同步转速以及低同步转速时的稳定性的同时,解决采用虚拟电压注入法后感应电机无速度传感器驱动系统的转速观测精度下降的问题。

Description

一种基于误差补偿的感应电机零频稳定控制方法
技术领域
本发明属于电机无速度传感器控制领域,更具体地,涉及一种基于误差补偿的感应电机零频稳定控制方法。
背景技术
感应电机(Induction Motor,简称IM,包括单相、三相和多相感应电机)及其驱动变频器(包括低压、中压和高压变频器)作为装备制造业的重要组成部分,被广泛应用于国民经济的各行各业,总体来说,无速度传感器驱动系统的优点包括:驱动系统成本低,系统集成配件少,尺寸小,可靠性高以及维护方便等。其缺点主要是电机低速运行带载能力弱,转速精度差以及低频运行不稳定性。
近年来,国内外学者对无速度传感器驱动系统中磁链观测和转速观测问题进行了深入全面的研究。根据每个算法特性,可以将以上观测方法分为三类:(1)基于电机模型构建的观测器:该方法已经在工业领域得到了广泛的应用,谐波含量、观测精度、动态响应等性能也已经得到了大幅度提升,但该类观测器在零同步转速时的转子转速不可观测问题一直没有得到解决,在实际应用中导致了感应电机无速度传感器驱动系统在电机低同步转速和零同步转速运行的不稳定问题。(2)基于电机各向异性构建的观测器:该方法利用感应电机的转子槽谐波、转子凸极效应、转子漏感等电机异向性,从采样得到的电流电压信号中提取出转子位置信息。该方法能够保证感应电机无速度传感器控制系统在电机同步转速等于零或者接近零时的长时间带载稳定运行,但是这类方法利用了感应电机的各向异性,依赖电机设计,而且存在转矩波动、噪声等问题,因此这类转速和磁链观测方法不适合广泛的工业应用。(3)虚拟电压注入法:该方法能够同时解决以上两类方法存在的问题。通过在观测器中注入虚拟电压信号,就能够解决感应电机无速度传感器驱动系统在电机低同步转速和零同步转速运行的不稳定问题,并且该方法仅在观测器中注入信号,不依赖电机各向异性,因此能够适合各种类型的感应电机,但是该方法存在会导致转速观测精度下降的问题。
发明内容
针对现有技术的缺陷,本发明的目的在于解决现有采用虚拟电压注入法后感应电机无速度传感器驱动系统的转速观测精度下降的技术问题。
为实现上述目的,第一方面,本发明提供了一种基于误差补偿的感应电机零频稳定控制方法,该方法基于现有的感应电机无速度传感器驱动系统,在αβ坐标系下电机的定子电压指令输入值
Figure BDA0001696276350000021
和转子电压补偿模块的输入
Figure BDA0001696276350000022
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系,包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure BDA0001696276350000024
S3.通过PARK反变换将
Figure BDA0001696276350000026
由dq坐标系变换到αβ坐标系,得到αβ坐标系下的转子补偿电压
Figure BDA0001696276350000027
Figure BDA0001696276350000028
并结合αβ坐标系下虚拟电压注入模块的输出值
Figure BDA0001696276350000029
来构成转子电压补偿模块的输出值
Figure BDA00016962763500000210
Figure BDA00016962763500000211
S4.基于
Figure BDA00016962763500000212
构建磁链观测器的动态数学模型;
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA0001696276350000031
使用磁链观测器观测转子磁链的旋转角度
S6.观测出的转子转速
Figure BDA0001696276350000033
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure BDA0001696276350000034
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制。
具体地,所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure BDA0001696276350000036
Figure BDA0001696276350000037
Figure BDA0001696276350000038
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000039
分别为d轴、q轴下的转子电流观测值。
具体地,转子补偿电压
Figure BDA00016962763500000311
的计算公式如下:
Figure BDA0001696276350000041
Figure BDA0001696276350000042
Figure BDA0001696276350000043
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000044
分别为d轴、q轴下的定子磁链观测值;
Figure BDA0001696276350000045
分别为d轴、q轴下的转子磁链观测值。
具体地,基于
Figure BDA0001696276350000046
构建磁链观测器的动态数学模型表示如下:
Figure BDA0001696276350000047
其中,
Figure BDA0001696276350000048
Figure BDA0001696276350000049
Figure BDA00016962763500000410
Figure BDA00016962763500000411
其中,
Figure BDA00016962763500000412
分别为α轴、β轴下的定子磁链观测值;
Figure BDA00016962763500000413
分别为α轴、β轴下的转子磁链观测值;
Figure BDA00016962763500000414
分别为α轴、β轴下的定子电流观测值;为转速观测器观测感应电机转子转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
具体地,所述基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA0001696276350000051
使用磁链观测器观测转子磁链的旋转角度
Figure BDA0001696276350000052
计算公式如下:
Figure BDA0001696276350000053
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;i、i分别为α轴、β轴下定子电流实际测量值;
Figure BDA0001696276350000055
分别为α轴、β轴下定子电流观测值;
Figure BDA0001696276350000056
分别为α轴、β轴下的转子磁链观测值;S1分别为
Figure BDA0001696276350000057
的时间积分值。
第二方面,本发明实施例提供了一种基于误差补偿的感应电机零频稳定控制方法,该方法基于现有的感应电机无速度传感器驱动系统,通过在dq坐标系下电机的定子电压指令输入值usd、usq和转子电压补偿模块的输入
Figure BDA0001696276350000058
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系,包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure BDA00016962763500000510
S3.利用
Figure BDA00016962763500000511
和dq坐标系下虚拟电压注入模块的输出值
Figure BDA00016962763500000512
构成转子电压补偿模块的输出值
Figure BDA00016962763500000513
Figure BDA00016962763500000514
S4.基于构建磁链观测器的动态数学模型;
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA0001696276350000061
使用观测出的电机转子转速计算转子磁链的旋转角度
Figure BDA0001696276350000062
S6.观测出的转子转速
Figure BDA0001696276350000063
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure BDA0001696276350000064
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制。
具体地,所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure BDA0001696276350000065
Figure BDA0001696276350000066
Figure BDA0001696276350000067
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000069
分别为d轴、q轴下的转子观测电流。
具体地,转子补偿电压
Figure BDA00016962763500000610
Figure BDA00016962763500000611
的计算公式如下:
Figure BDA0001696276350000071
Figure BDA0001696276350000072
Figure BDA0001696276350000073
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000074
分别为d轴、q轴下的定子磁链观测值;分别为d轴、q轴下的转子电流观测值。
具体地,所述基于构建磁链观测器的动态数学模型如下:
Figure BDA0001696276350000077
其中,
Figure BDA0001696276350000078
Figure BDA0001696276350000079
Figure BDA00016962763500000710
Figure BDA00016962763500000711
其中,分别为d轴、q轴下的定子磁链观测值;
Figure BDA00016962763500000713
分别为d轴、q轴下的转子磁链观测值;分别为d轴、q轴下的定子电流观测值;
Figure BDA00016962763500000715
为转速观测器观测感应电机转子转速,ωe为同步转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
具体地,所述基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA0001696276350000081
使用观测出的电机转子转速计算转子磁链的旋转角度
Figure BDA0001696276350000082
计算公式如下:
Figure BDA0001696276350000083
Figure BDA0001696276350000084
Figure BDA0001696276350000085
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;isd、isq分别为d轴、q轴下定子电流实际测量值;
Figure BDA0001696276350000086
分别为d轴、q轴下定子电流观测值;
Figure BDA0001696276350000087
分别为d轴、q轴下的转子磁链观测值;S2分别为
Figure BDA0001696276350000088
的时间积分值,ωs为滑差转速,Rr为电机转子电阻,Lr、Lm分别为电机转子侧电感和电机互感,S3代表对
Figure BDA0001696276350000089
的时间积分。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,通过在虚拟电压注入模块和磁链转速观测器之间增加一转子电压补偿模块,实现在保证感应电机无速度传感器驱动系统零同步转速以及低同步转速时的稳定性的同时,解决采用虚拟电压注入法后感应电机无速度传感器驱动系统的转速观测精度下降的问题。具有以下有益效果:
(1)在不对电机本体进行信号注入的情况下,实现无速度传感器驱动系统控制的感应电机在零同步转速或低同步转速下均能够输出200%电机额定转矩,并且在低同步转速下能够提升转速观测精度,在零同步转速时刻的转速估计误差不可避免,但能稳定运行。
(2)在不对电机本体进行信号注入的情况下,实现无速度传感器驱动系统控制的感应电机在0%电机额定转矩和零转子转速下长时间稳定运行,且转速观测误差为0,并能够在该状态长时间运行后正常加速启动。
(3)在不对电机本体进行信号注入的情况下,实现无速度传感器驱动系统控制的感应电机在负载为200%电机额定转矩不变的情况下,以任意加减速时间进行电机转速的正反转切换。
(4)转速估计误差的存在范围为同步转速在-0.5Hz-0.5Hz之间。其余范围内(大于0.5Hz和小于-0.5Hz时)均无转速估计误差。
附图说明
图1为本发明实施例一提供的基于虚拟电压注入的感应电机无速度传感器驱动控制系统结构示意图。
图2为本发明实施例一提供的离线标定法的试验平台结构示意图。
图3为本发明实施例二提供的基于虚拟电压注入的感应电机无速度传感器驱动控制系统结构示意图。
图4为本发明实施例提供的100%额定转矩时定子电流、观测转子转速、实际转子转速随时间变化的曲线图。
图5为本发明实施例提供的200%额定转矩时定子电流、观测转子转速、实际转子转速随时间变化的曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1为本发明实施例一提供的基于虚拟电压注入的感应电机无速度传感器驱动控制系统结构示意图。
感应电机无速度传感器驱动系统的硬件部分包括:三相电压源型逆变器和感应电机。三相交流电源经过不控整流得到直流母线电压UDC,供给电压源型逆变器,并利用逆变器来控制感应电机进行转矩和转速的控制。三相电压源型逆变器中包括电压和电流传感器。
感应电机无速度传感器驱动系统的软件部分包括:3相静止坐标/2相静止坐标变换模块,2相同步旋转坐标/2相静止坐标变换模块,电压空间矢量脉宽调制模块,电流PI调节模块,转速PI调节模块,磁链电流指令给定模块,转速指令给定模块,磁链和转速观测器模块、虚拟电压注入模块和转子电压补偿模块。
实施例一中控制系统在αβ坐标系下电机的定子电压指令输入值和转子电压补偿模块的输入
Figure BDA0001696276350000102
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块。下面介绍整个系统的控制方法,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系。包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure BDA0001696276350000103
Figure BDA0001696276350000111
Figure BDA0001696276350000112
Figure BDA0001696276350000113
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000114
分别为d轴、q轴下的转子电流观测值。
所述离线标定法计算等效补偿转子电阻Rrc,具体包括以下步骤:
(2-1)构建试验平台,图2为本发明实施例一提供的离线标定法的试验平台结构示意图。如图2所示,试验平台包括:变频器*2、电流传感器组*2、转速传感器*1、测功机*1、转矩传感器*1、被测电机*1。变频器一个起控制被测电机的作用,另一个起控制测功机的作用;电流传感器四个,其中两个用于测量被测电机的定子三相电流,另外两个用于测量测功机的定子三相电流;转速传感器用于测量电机实际转速;转矩传感器用于测量测功机的输出转矩;测功机用于输出固定负载转矩,运行在转矩控制模式;被测电机由基于虚拟电压注入的感应电机无速度传感器驱动系统控制,并运行在转速控制模式;两个电机进行同轴刚性连接。
(2-2)控制测功机输出2TL,其中,TL为被测电机额定负载转矩,被测电机转速由N1rpm逐渐降为-N1rpm,下降步长为n rpm。在每个转速下,调整Rrc,使得电机转速指令与电机实际转速相同,并记录下此时d轴定子电流值、q轴定子电流值和此时的Rrc
(2-3)测功机的输出依次减小0.1TL直至0,重复测量被测电机转速由N1rpm以步长n降低至-N1rpm过程的d轴定子电流值、q轴定子电流值和Rrc
(2-4)将上述数据整理为表格。
(2-5)控制系统在线运行时,根据当前d轴电流指令和q轴电流指令查询表格来确定Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure BDA0001696276350000121
Figure BDA0001696276350000122
Figure BDA0001696276350000123
Figure BDA0001696276350000124
Figure BDA0001696276350000125
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000126
分别为d轴、q轴下的定子磁链观测值;
Figure BDA0001696276350000127
分别为d轴、q轴下的转子磁链观测值。
S3.利用
Figure BDA0001696276350000129
通过PARK反变换由dq坐标系变换到αβ坐标系,得到αβ坐标系下的转子补偿电压
Figure BDA00016962763500001210
Figure BDA00016962763500001211
并结合αβ坐标系下虚拟电压注入模块的输出值
Figure BDA00016962763500001212
来构成转子电压补偿模块的输出值
Figure BDA00016962763500001213
Figure BDA00016962763500001214
S4.基于
Figure BDA00016962763500001215
构建磁链观测器的动态数学模型;
Figure BDA00016962763500001216
其中:
Figure BDA00016962763500001217
Figure BDA0001696276350000131
Figure BDA0001696276350000132
其中,
Figure BDA0001696276350000133
分别为α轴、β轴下的定子磁链观测值;
Figure BDA0001696276350000134
分别为α轴、β轴下的转子磁链观测值;
Figure BDA0001696276350000135
分别为α轴、β轴下的定子电流观测值;
Figure BDA0001696276350000136
为转速观测器观测感应电机转子转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA0001696276350000137
使用磁链观测器观测转子磁链的旋转角度
Figure BDA0001696276350000138
Figure BDA0001696276350000139
Figure BDA00016962763500001310
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;i、i分别为α轴、β轴下定子电流实际测量值;
Figure BDA00016962763500001311
分别为α轴、β轴下定子电流观测值;
Figure BDA00016962763500001312
分别为α轴、β轴下的转子磁链观测值;S1分别为
Figure BDA00016962763500001313
的时间积分值。
S6.观测出的转子转速
Figure BDA00016962763500001314
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure BDA00016962763500001315
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制。
步骤S6包括以下步骤:
(1)观测出的感应电机转子转速
Figure BDA00016962763500001316
作为转速PI调节模块的反馈量,与对应的转速指令
Figure BDA00016962763500001317
做差后进行转速PI控制;
(2)观测出的磁链旋转角度
Figure BDA0001696276350000141
被用于2相同步旋转坐标/2相静止坐标变换模块中的坐标变换计算;
(3)转速PI调节模块的输出
Figure BDA0001696276350000142
作为q轴电流PI调节模块的指令,磁链电流指令给定模块的输出
Figure BDA0001696276350000143
作为d轴电流PI调节模块的指令;由电流传感器采样得到的感应电机两相电流iSU、iSV,输入3相静止坐标/2相静止坐标变换模块,然后输出
Figure BDA0001696276350000144
至2相同步旋转坐标/2相静止坐标变换模块,最后得到2相同步旋转坐标系下的d轴电流isd和q轴电流isq,并将得到的d轴电流和q轴电流分别作为d轴电流PI调节器和q轴电流PI调节器的反馈量,与对应的磁链电流指令
Figure BDA0001696276350000145
Figure BDA0001696276350000146
做差后进行电流PI控制;
(4)d轴和q轴电流PI调节模块的输出usd和usq至2相同步旋转坐标/2相静止坐标变换模块,其将2相同步旋转坐标系下的电机输入电压指令转换为2相静止坐标系下的电机输入电压指令
Figure BDA0001696276350000147
(5)将
Figure BDA0001696276350000148
输出至电压空间矢量脉宽调制模块,产生能够控制开关器件SA、SB、SC的开关信号,进而达到控制感应电机转速和转矩的目的。
图3为本发明实施例二提供的基于虚拟电压注入的感应电机无速度传感器驱动控制系统结构示意图。
实施例二中控制系统在dq坐标系下电机的定子电压指令输入值usd、usq和转子电压补偿模块的输入
Figure BDA0001696276350000149
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块。下面介绍整个系统的控制方法,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系。包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure BDA0001696276350000151
Figure BDA0001696276350000152
Figure BDA0001696276350000153
Figure BDA0001696276350000154
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000155
分别为d轴、q轴下的转子电流观测值。
所述离线标定法计算等效补偿转子电阻Rrc,具体包括以下步骤:
(2-1)构建试验平台,试验平台与实施例一相同,包括:变频器*2、电流传感器组*2、转速传感器*1、测功机*1、转矩传感器*1、被测电机*1。变频器一个起控制被测电机的作用,另一个起控制测功机的作用;电流传感器两个用于测量被测电机的定子三相电流,另外两个用于测量测功机的定子三相电流;转速传感器用于测量电机实际转速;转矩传感器用于测量测功机的输出转矩;测功机用于输出固定负载转矩,运行在转矩控制模式;被测电机由基于虚拟电压注入的感应电机无速度传感器驱动系统控制,并运行在转速控制模式;两个电机进行同轴刚性连接。
(2-2)控制测功机输出2TL,其中,TL为被测电机额定负载转矩,被测电机转速由N1rpm逐渐降为-N1rpm,下降步长为n rpm。在每个转速下,调整Rrc,使得电机转速指令与电机实际转速相同,并记录下此时d轴定子电流值、q轴定子电流值和此时的Rrc
(2-3)测功机的输出依次减小0.1TL直至0,重复测量被测电机转速由N1rpm以步长n降低至-N1rpm过程的d轴定子电流值、q轴定子电流值和Rrc
(2-4)将上述数据整理为表格。
(2-5)控制系统在线运行时,根据当前d轴电流指令和q轴电流指令查询表格来确定Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure BDA0001696276350000161
Figure BDA0001696276350000162
Figure BDA0001696276350000163
Figure BDA0001696276350000164
Figure BDA0001696276350000165
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure BDA0001696276350000166
分别为d轴、q轴下的定子磁链观测值;
Figure BDA0001696276350000167
分别为d轴、q轴下的转子磁链观测值。
S3.利用
Figure BDA0001696276350000168
和dq坐标系下虚拟电压注入模块的输出值
Figure BDA0001696276350000169
构成转子电压补偿模块的输出值
Figure BDA00016962763500001611
S4.基于
Figure BDA00016962763500001612
构建磁链观测器的动态数学模型;
Figure BDA0001696276350000171
Figure BDA0001696276350000172
Figure BDA0001696276350000173
Figure BDA0001696276350000175
其中,
Figure BDA0001696276350000176
分别为d轴、q轴下的定子磁链观测值;分别为d轴、q轴下的转子磁链观测值;分别为d轴、q轴下的定子电流观测值;为转速观测器观测感应电机转子转速,ωe为同步转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure BDA00016962763500001710
使用观测出的电机转子转速计算转子磁链的旋转角度
Figure BDA00016962763500001711
Figure BDA00016962763500001712
Figure BDA00016962763500001713
Figure BDA00016962763500001714
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;isd、isq分别为d轴、q轴下定子电流实际测量值;分别为d轴、q轴下定子电流观测值;
Figure BDA0001696276350000181
分别为d轴、q轴下的转子磁链观测值;S2分别为
Figure BDA0001696276350000182
的时间积分值,ωs为滑差转速,Rr为电机转子电阻,Lr、Lm分别为电机转子侧电感和电机互感,S3代表对
Figure BDA0001696276350000183
的时间积分。
S6.观测出的转子转速
Figure BDA0001696276350000184
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure BDA0001696276350000185
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制。
步骤S6的具体实施过程与实施例一相同。
图4为本发明实施例提供的100%额定转矩时定子电流、观测转子转速、实际转子转速随时间变化的曲线图。iSU曲线为定子电流曲线,
Figure BDA0001696276350000186
为观测转子转速曲线,ωr为实际转子转速曲线。负载转矩为被测电机的100%额定负载转矩。转子转速指令由900rpm下降为-900rpm,当转子转速为负数时,电机运行在电动状态,当电机转子转速为正数时,电机运行在发电状态。从图中可以看出:(1)定子电流频率为0Hz时,系统能够稳定运行;(2)在定子电流频率大于0.5Hz和小于-0.5Hz时,
Figure BDA0001696276350000187
与ωr之间没有误差。
图5为本发明实施例提供的200%额定转矩时定子电流、观测转子转速、实际转子转速随时间变化的曲线图。iSU曲线为定子电流曲线,为观测转子转速曲线,ωr为实际转子转速曲线。负载转矩为被测电机的200%额定负载转矩。转子转速指令由120rpm下降为0rpm,再升到120rpm。从图中可以看出:定子电流频率为0Hz时,系统能够稳定运行。
以上,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。

Claims (8)

1.一种基于误差补偿的感应电机零频稳定控制方法,在感应电机无速度传感器驱动系统αβ坐标系下电机的定子电压指令输入值
Figure FDA0002229855160000011
和转子电压补偿模块的输入
Figure FDA0002229855160000012
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系,其特征在于,包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure FDA0002229855160000013
Figure FDA0002229855160000014
S3.通过PARK反变换将
Figure FDA0002229855160000015
由dq坐标系变换到αβ坐标系,得到αβ坐标系下的转子补偿电压
Figure FDA0002229855160000017
Figure FDA0002229855160000018
并结合αβ坐标系下虚拟电压注入模块的输出值
Figure FDA0002229855160000019
来构成转子电压补偿模块的输出值
Figure FDA00022298551600000110
Figure FDA00022298551600000111
S4.基于
Figure FDA00022298551600000112
构建磁链观测器的动态数学模型;
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure FDA00022298551600000113
使用磁链观测器观测转子磁链的旋转角度
Figure FDA00022298551600000114
S6.观测出的转子转速
Figure FDA00022298551600000115
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure FDA00022298551600000116
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制;
所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure FDA0002229855160000021
Figure FDA0002229855160000022
Figure FDA0002229855160000023
Figure FDA0002229855160000024
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure FDA0002229855160000025
分别为d轴、q轴下的转子电流观测值。
2.如权利要求1所述的控制方法,其特征在于,转子补偿电压
Figure FDA0002229855160000026
Figure FDA0002229855160000027
的计算公式如下:
Figure FDA0002229855160000028
Figure FDA0002229855160000029
Figure FDA00022298551600000210
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure FDA0002229855160000031
分别为d轴、q轴下的定子磁链观测值;
Figure FDA0002229855160000032
分别为d轴、q轴下的转子磁链观测值。
3.如权利要求1所述的控制方法,其特征在于,基于
Figure FDA0002229855160000033
构建磁链观测器的动态数学模型表示如下:
Figure FDA0002229855160000034
其中,
Figure FDA0002229855160000035
Figure FDA0002229855160000036
Figure FDA0002229855160000037
Figure FDA0002229855160000038
其中,
Figure FDA0002229855160000039
分别为α轴、β轴下的定子磁链观测值;
Figure FDA00022298551600000310
分别为α轴、β轴下的转子磁链观测值;
Figure FDA00022298551600000311
分别为α轴、β轴下的定子电流观测值;
Figure FDA00022298551600000312
为转速观测器观测感应电机转子转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
4.如权利要求1所述的控制方法,其特征在于,所述基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速使用磁链观测器观测转子磁链的旋转角度
Figure FDA0002229855160000042
计算公式如下:
Figure FDA0002229855160000044
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;i、i分别为α轴、β轴下定子电流实际测量值;
Figure FDA0002229855160000045
分别为α轴、β轴下定子电流观测值;
Figure FDA0002229855160000046
分别为α轴、β轴下的转子磁链观测值;S1分别为
Figure FDA0002229855160000047
的时间积分值。
5.一种基于误差补偿的感应电机零频稳定控制方法,通过在感应电机无速度传感器驱动系统dq坐标系下电机的定子电压指令输入值usd、usq和转子电压补偿模块的输入
Figure FDA0002229855160000048
之间增加一虚拟电压注入模块,在虚拟电压注入模块、磁链转速观测器之间增加一转子电压补偿模块,αβ坐标系是2相静止坐标系,dq坐标系是2相同步旋转坐标系,其特征在于,包括以下步骤:
S1.基于虚拟电压注入值或者离线标定法计算等效补偿转子电阻Rrc
S2.基于等效补偿转子电阻Rrc计算转子补偿电压
Figure FDA0002229855160000049
S3.利用
Figure FDA00022298551600000411
和dq坐标系下虚拟电压注入模块的输出值
Figure FDA00022298551600000412
Figure FDA00022298551600000413
构成转子电压补偿模块的输出值
Figure FDA00022298551600000414
Figure FDA00022298551600000415
S4.基于
Figure FDA00022298551600000416
构建磁链观测器的动态数学模型;
S5.基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速
Figure FDA0002229855160000051
使用观测出的电机转子转速计算转子磁链的旋转角度
S6.观测出的转子转速
Figure FDA0002229855160000053
被用于转速PI调节模块和磁链观测器,观测出的转子磁链旋转角度
Figure FDA0002229855160000054
用于2相同步旋转坐标/2相静止坐标变换模块,从而实现感应电机无速度传感器转速和转矩的控制;
所述基于虚拟电压注入值计算等效补偿转子电阻Rrc,具体包括以下步骤:
(1-1)利用dq坐标系下的电机定子电压指令输入值usd、usq和虚拟电压注入模块中的比例系数k,计算d轴下的虚拟电压注入值usd_inj和q轴下的虚拟电压注入值usq_inj
usd_inj=(k-1)usd
usq_inj=(k-1)usq
(1-2)利用usd_inj和usq_inj计算等效补偿转子电阻Rrc
Figure FDA0002229855160000055
Figure FDA0002229855160000056
Figure FDA0002229855160000057
Figure FDA0002229855160000058
其中,ωe为同步转速;ωs为滑差转速;Rs、Rr、Rrc分别为电机定子电阻、转子电阻、等效补偿转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;分别为d轴、q轴下的转子电流观测值。
6.如权利要求5所述的控制方法,其特征在于,转子补偿电压
Figure FDA0002229855160000061
Figure FDA0002229855160000062
的计算公式如下:
Figure FDA0002229855160000063
Figure FDA0002229855160000064
其中,Rrc为等效补偿转子电阻,Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感;
Figure FDA0002229855160000066
分别为d轴、q轴下的定子磁链观测值;
Figure FDA0002229855160000067
分别为d轴、q轴下的转子磁链观测值。
7.如权利要求5所述的控制方法,其特征在于,所述基于
Figure FDA0002229855160000068
构建磁链观测器的动态数学模型如下:
Figure FDA0002229855160000069
其中,
Figure FDA00022298551600000610
Figure FDA00022298551600000611
Figure FDA00022298551600000613
其中,
Figure FDA00022298551600000614
分别为d轴、q轴下的定子磁链观测值;分别为d轴、q轴下的转子磁链观测值;
Figure FDA0002229855160000071
分别为d轴、q轴下的定子电流观测值;
Figure FDA0002229855160000072
为转速观测器观测感应电机转子转速,ωe为同步转速;Rs、Rr分别为电机定子电阻和转子电阻;Ls、Lr、Lm分别为电机定子侧电感、电机转子侧电感和电机互感。
8.如权利要求5所述的控制方法,其特征在于,所述基于磁链观测器的动态数学模型,使用转速观测器观测感应电机转子转速使用观测出的电机转子转速计算转子磁链的旋转角度计算公式如下:
Figure FDA0002229855160000075
Figure FDA0002229855160000076
Figure FDA0002229855160000077
其中,kp、ki分别为转速观测器的比例环节增益和积分环节增益;isd、isq分别为d轴、q轴下定子电流实际测量值;
Figure FDA0002229855160000078
分别为d轴、q轴下定子电流观测值;分别为d轴、q轴下的转子磁链观测值;S2分别为的时间积分值,ωs为滑差转速,Rr为电机转子电阻,Lr、Lm分别为电机转子侧电感和电机互感,S3代表对
Figure FDA00022298551600000711
的时间积分。
CN201810614110.5A 2018-06-14 2018-06-14 一种基于误差补偿的感应电机零频稳定控制方法 Active CN108718165B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810614110.5A CN108718165B (zh) 2018-06-14 2018-06-14 一种基于误差补偿的感应电机零频稳定控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810614110.5A CN108718165B (zh) 2018-06-14 2018-06-14 一种基于误差补偿的感应电机零频稳定控制方法

Publications (2)

Publication Number Publication Date
CN108718165A CN108718165A (zh) 2018-10-30
CN108718165B true CN108718165B (zh) 2020-01-10

Family

ID=63912982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810614110.5A Active CN108718165B (zh) 2018-06-14 2018-06-14 一种基于误差补偿的感应电机零频稳定控制方法

Country Status (1)

Country Link
CN (1) CN108718165B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111175603B (zh) * 2020-01-06 2022-03-08 广州视源电子科技股份有限公司 变频器的老化方法、老化装置、存储介质和处理器
CN112615578B (zh) * 2020-12-17 2023-11-28 深圳市迈凯诺电气股份有限公司 一种异步电机开环矢量控制系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0496681A (ja) * 1990-08-13 1992-03-30 Meidensha Corp 誘導電動機のベクトル制御装置
CN1647359A (zh) * 2002-04-02 2005-07-27 株式会社安川电机 交流电动机的无传感器控制装置及控制方法
CN101478283A (zh) * 2009-01-22 2009-07-08 北京清能华福风电技术有限公司 不平衡电网电压下双馈异步风力发电机转子侧逆变器控制方法
JP2015204651A (ja) * 2014-04-11 2015-11-16 株式会社明電舎 誘導電動機の制御装置および制御方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0496681A (ja) * 1990-08-13 1992-03-30 Meidensha Corp 誘導電動機のベクトル制御装置
CN1647359A (zh) * 2002-04-02 2005-07-27 株式会社安川电机 交流电动机的无传感器控制装置及控制方法
CN101478283A (zh) * 2009-01-22 2009-07-08 北京清能华福风电技术有限公司 不平衡电网电压下双馈异步风力发电机转子侧逆变器控制方法
JP2015204651A (ja) * 2014-04-11 2015-11-16 株式会社明電舎 誘導電動機の制御装置および制御方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Inverter Nonlinear Error Compensation Using Feedback Gains and Self-Tuning Estimated Current Error in Adaptive Full-Order Observer;Wei Sun 等;《IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS》;20160229;第52卷(第1期);472-482 *
Zero Stator Current Frequency Operation of Speed-Sensorless Induction Motor Drives Using Stator Input Voltage Error for Speed Estimation;Wei Sun 等;《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》;20160331;第63卷(第3期);1490-1498 *
感应电机无速度传感器矢量控制低速运行技术研究;孙伟;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;中国学术期刊(光盘版)电子杂志社;20180131(第01期);全文 *

Also Published As

Publication number Publication date
CN108718165A (zh) 2018-10-30

Similar Documents

Publication Publication Date Title
CN108599651B (zh) 基于虚拟电压注入的感应电机无速度传感器驱动控制方法
CN108092567B (zh) 一种永磁同步电动机转速控制系统及方法
CN110022106B (zh) 一种基于高频信号注入的永磁同步电机无位置传感器控制方法
CN103825525B (zh) 一种改进的无传感器永磁同步电机速度估测方法
CN101931361A (zh) 一种感应电动机的矢量控制装置
CN108988725B (zh) 一种采用改进复矢量pi控制器的永磁同步电机电流谐波抑制系统及方法
CN111555680A (zh) 一种永磁同步电机无差拍预测电流控制方法
Tian et al. Rotor position estimation of sensorless PMSM based on extented Kalman filter
CN110798116B (zh) 一种基于电枢模型独立前馈补偿的电机矢量复合控制器
CN111342719B (zh) 一种无速度传感器驱动的异步电机控制方法
CN110429891B (zh) 一种无位置传感器永磁电机直驱发电控制方法
CN108718165B (zh) 一种基于误差补偿的感应电机零频稳定控制方法
CN113422550B (zh) 基于复矢量解耦及延时补偿的高速电机低载波比控制方法
CN113676088B (zh) 带有谐波抑制的永磁同步电机无速度传感器控制方法
Chi Position-sensorless control of permanent magnet synchronous machines over wide speed range
CN108540031B (zh) 无轴承同步磁阻电机的转速估计方法及转矩控制系统
CN114301361B (zh) 一种基于母线电流控制的无电解电容永磁同步电机驱动系统控制方法
CN116365937A (zh) 一种高速永磁同步电机无位置传感器控制方法
CN109560741A (zh) 一种基于测量误差补偿器的永磁同步电机系统及补偿方法
CN114337426A (zh) 一种d-q轴静止坐标系下永磁同步电机偏差解耦控制方法
CN112701970B (zh) 一种低载波比下的pmsm精确离散自抗扰控制方法
CN114465530A (zh) 永磁同步电机的速度控制方法及系统
Hu et al. Vector Control of Permanent Magnet Synchronous Motor Based on New MRAS
KR20180080730A (ko) 모터 제어 시스템 및 방법
CN113141139A (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