CN107769636B - 一种无位置传感器永磁同步电机转子位置检测方法 - Google Patents

一种无位置传感器永磁同步电机转子位置检测方法 Download PDF

Info

Publication number
CN107769636B
CN107769636B CN201710997952.9A CN201710997952A CN107769636B CN 107769636 B CN107769636 B CN 107769636B CN 201710997952 A CN201710997952 A CN 201710997952A CN 107769636 B CN107769636 B CN 107769636B
Authority
CN
China
Prior art keywords
rotor
motor
phase
permanent magnet
speed
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
CN201710997952.9A
Other languages
English (en)
Other versions
CN107769636A (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.)
Northwest University of Technology
Original Assignee
Northwest University of 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 Northwest University of Technology filed Critical Northwest University of Technology
Priority to CN201710997952.9A priority Critical patent/CN107769636B/zh
Publication of CN107769636A publication Critical patent/CN107769636A/zh
Application granted granted Critical
Publication of CN107769636B publication Critical patent/CN107769636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/14Electronic commutators
    • H02P6/16Circuit arrangements for detecting position
    • H02P6/18Circuit arrangements for detecting position without separate position detecting elements
    • H02P6/182Circuit arrangements for detecting position without separate position detecting elements using back-emf in windings
    • 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
    • H02P21/18Estimation of position or speed

Abstract

本发明涉及一种无位置传感器永磁同步电机转子位置检测方法,采用信号处理的方法对三相绕组反电势的变化过程进行数字处理,提取出电机旋转过程中反向电动势的某些相对稳定的特征位置,诸如零点与极值点,得到位置特征点,对处于特征点之间的转子位置信息采用扩展卡尔曼滤波方法估计,最终得到转子的准确位置。该方法数据采集方便,计算量小,实时性好,适用于多种电机控制器实现。

Description

一种无位置传感器永磁同步电机转子位置检测方法
技术领域
本发明涉及一种检测无位置传感器永磁同步电机转子位置检测方法,属于电机控制技术领域。
背景技术
永磁同步电机要稳定运行,必须实时检测到永磁转子位置,实现闭环控制。想让永磁同步电机处于闭环运行,驱动系统必须实时知道永磁转子的磁极位置。只有这样才能满足对速度和电流的精确控制,因此转子位置和转速信号的获取是整个驱动系统非常重要的环节。
高精度的永磁同步电机系统对速度控制和位置控制提出很高的要求,对传感器的要求也相应地提高。由于机械传感器(如编码器、霍尔感应器、旋转变压器等)的安装带来了系统成本增加、体积增大、受工作环境影响、可靠性降低等缺陷。无传感器永磁同步电机系统是在不安装位置传感器的前提下,利用检测到的电机电压、电流和数学模型估算出转子位置和转速,即利用电气特性来反映出电机的机械运动特性,由于无需改造电机,省去昂贵的机械传感器、降低维护费用和不怕粉尘与潮湿环境影响等优点,而且从根本上避免了由于加装机械传感器引起的电机轴抖动、机械惯量增加等硬件传感器不可避免的缺陷,因此无位置传感器的永磁同步电机系统获得广泛的应用。
永磁同步电机无位置传感器矢量控制的方法按照其理论基础可以划分为两类,一类是基于电机基波模型的无位置传感器控制方法,另一类是基于电机谐波模型的无位置传感器控制方法。基于电机基波模型方法的基本思想是利用电机反电势和转子位置的关系来估算转子位置。主要包括反电势直接计算法、模型参考自适应法、观测器法以及卡尔曼滤波法等。该方法本质上都是从电机的反电势中获取电机转子位置和转速,一般适用于中高速范围。基于电机谐波模型的无位置传感器控制方法基本思想是利用转子的“物理凸极”结构对定子电信号的调制作用来估算转子位置,主要包括电感测量法、电压脉冲法、高频信号注入法和载波频率法。该类方法可以在低速甚至零速时估算转子位置,但是当电机运行于高速区时,反电动势过大,电压方程中的旋转分量不可忽略,使得位置估算精度降低,稳定性变差。
发明内容
要解决的技术问题
常用的无位置传感器转子位置检测方法有:全降维观测器法、模型参考自适应观测器法、滑模观测器法等,然而上述方法均需要大量的软件计算,给实际应用带来一定的困难,从而影响矢量控制方法的精度。
为了避免上述现有技术的不足之处,本发明提出一种无位置传感器永磁同步电机转子位置检测方法。
技术方案
一种无位置传感器永磁同步电机转子位置检测方法,其特征在于步骤如下:
步骤1:当电机从速度开环、电流闭环的启动状态切换到速度、电流双闭环的运行状态,通过电流传感器采集电机三相绕组电流,计算电机的三相反向电动势U、V、W分别为:
式中,DA、DB、DC分别是控制三相桥臂的PWM信号的占空比,Ud是施加到绕组上驱动电压,iA、iB、iC是相电流,Ro是绕组的等效直流电阻;
步骤2:通过三相反电动势A、B、C、A-B、B-C、C-A与电角度轴的交点将电角度轴分割为30度等长的12个区间,在一个电周期内得到12个准确的转子位置的特征点,对处于特征点之间的转子位置信息采用扩展卡尔曼滤波方法估计;
步骤3:由电流传感器采样定子电流,一方面通过dq变换输入到电流控制器,另一方面输入到速度-位置观测器中由反电势计算得到转子位置和速度,再分别输入到坐标变换和速度控制器中,构成空间矢量双闭环控制系统;
永磁同步电机在α、β坐标系下的电压方程为:
式中,iα、iβ为电流在α、β轴上的分量,ωe是转子旋转的电角速度、θe是转子旋转的电角度;
状态方程为:
其中,Rs为电机相电阻,Ls为电机相电感,λr为电机转子永磁体磁链,Np是电机的极对数,ωr是电机转子的角速度;
步骤4:选取状态变量、输入矢量、输出矢量为:x=[iα iβ ωr θe]T,u=[uα uβ]T, y=[iα iβ]T
设系统的采样周期为Tc,对公式(2)和(3)进行线性化和离散化得到:
其中,vk是测量噪声,wk是过程噪声,状态转移矩阵和测量矩阵分别为:
步骤5:扩展卡尔曼滤波器的输入为经过Clark变换得到的iα、iβ和经过Park逆变换得到的uα、uβ,输出为永磁同步电机的转子位置角和转速。
有益效果
本发明提出的一种无位置传感器永磁同步电机转子位置检测方法,该方法无需位置传感器,降低了系统成本;计算量小,实时性好,采集数据量小,便于电机控制器中多种微处理器(DSP、FPGA、单片机等)实现;适用于空调变频调速、水泵、风机等电机控制的应用场合。
附图说明
图1无位置传感器永磁同步电机控制系统结构框图
图2 2个电周期内的分区解析示意图
图3转子位置检测流程
具体实施方式
永磁同步电机转子在旋转过程中,绕组会因切割磁力线而产生反向电动势,而反向电动势与转子相位有直接的对应关系,因此可以利用反电势信息来监测转子的相位。在电机运行过程中凭反向电动势的幅值来确定转子的各个相位实现方法复杂,一方面各相电流分身存在有高频成份,另一方面转速不稳定时,反向电动势本身也不稳定,所以不适合直接从反向电动势的幅值来换算转子相位,可采用信号处理的方法对三相绕组反电势的变化过程进行数字处理,提取出电机旋转过程中反向电动势的某些相对稳定的特征位置,诸如零点与极值点,得到位置特征点,对处于特征点之间的转子位置信息采用扩展卡尔曼滤波方法估计,最终得到转子的准确位置。
(1)电机绕组反向电动势检测及转子位置特征点获取方法
电机的每相绕组为可等效为电阻电感串联支路(RL支路),绕组两端的电压由电感压降和电阻压降组成。R作为电机绕组的RL支路在电机转子转动过程中会呈现感应电动势,此时电机的U、V、W三相反向电动势分别为:
式中,DA、DB、DC分别是控制三相桥臂的PWM信号的占空比,Ud是施加到绕组上驱动电压,iA、iB、iC是相电流,Ro是绕组的等效直流电阻。根据上式可在整个电周期范围内测量到各相绕组的反电势。
检测到三相电机绕组的反向电动势A、B、C,根据这三相反电势的过零点以及极大值点可以将一个电角度周期以30度为步长分为12个区间。每个区间分界处是根据三相反电势的过零点和正负极大值点确定的。由于三相绕组在空间是对称安置的,因此某相绕组反电势的正负极大值点与另外两相反电势之差的过零点位置相同,在具体实现过程中某相的正负极极大值点可通过另外两相的差得到。如附图2中显示的六条曲线A、C、C、C-B、B-A、A-C与电角度轴的交点就将电角度轴分割为30度等长的区间。因此在一个电周期内通过反电势可以得到12个准确的转子位置的特征点,这 12个位置特征点对应附图2中曲线与0轴的12个交点(方格表示)。
(2)基于扩展卡尔曼滤波的无位置传感器位置估计算法
在一个电周期内通过反电势可以得到12个准确的转子位置的特征点。这些特征点之间的位置信息采取扩展卡尔曼滤波估计的方法进行递推。
永磁同步电机在两相定子静止坐标系(αβ坐标系)下的电压方程为:
iα、iβ为电流在α、β轴上的分量。ωe是转子旋转的电角速度、θe是转子旋转的电角度。
状态方程为:
其中,Rs为电机相电阻,Ls为电机相电感,λr为电机转子永磁体磁链,Np是电机的极对数,ωr是电机转子的角速度。考虑到电机的惯性时间常数比系统的采样周期大很多,可以认为
选取状态变量、输入矢量、观测输出矢量为:x=[Iα Iβ ωe θe]T,u=[uα uβ]T, y=[iα iβ]T
设系统的采样周期为Tc,对公式(2)和(3)进行线性化和离散化得到:
vk是测量噪声,wk是过程噪声,状态转移矩阵和测量矩阵分别为:
由上述方程可进行扩展卡尔曼器滤波器设计,进行转子位置的估计。
扩展卡尔曼滤波器的输入为经过Clark变换得到的iα、iβ和经过Park逆变换得到的uα、uβ,输出为永磁同步电机的转子位置角和转速。扩展卡尔曼滤波器输出经过闭环反馈送到PI控制器,输出电机的控制电压,通过空间矢量控制三相逆变器驱动永磁同步电机运转。
可以看出,可以利用电机绕组的反电势在α、β轴的分量,采取卡尔曼滤波方法估算出转子位置θe和转速ωe。计算过程中由于受到AD采样精度和量化误差等因素影响,计算得到的位置θe和实际转子位置θ存在误差。如前所述,在电机转子旋转一个电周期时通过反电势可以得到12个准确的转子位置特征点,系统检测到一个特征位置,相当于转子从上一特征点到该特征位置点旋转了30°,这些特征点之间的转子角度信息采取扩展卡尔曼滤波估计的方法进行递推。每检测到一次位置特征点(相当于电机旋转30°),利用该角度信息对卡尔曼滤波估计的角度信息进行校正,防止卡尔曼滤波的发散,同时也保证了转子角度估计的准确。
现结合实施例、附图对本发明作进一步描述:
(1)电机启动
启动时采取速度开环,电流闭环的启动方式,首先将电机带到一定转速,当反电势算法中估计的转子位置和速度信息稳定时,切换到反电势算法的矢量控制模式中,使得电机由低转速到高转速平稳过渡。
切换过程的实质是给定电流矢量与转子d轴间夹角θ由一个较小的值变化到90°的过程。通过给定角度连续增加的电压空间矢量将电机拖入同步,实现从速度开环、电流闭环的启动状态切换到速度、电流双闭环的运行状态,电压空间矢量变化的快慢和方向决定了起动速度的快慢和方向。
(2)电机绕组反向电动势检测及转子位置特征点获取
电机转子在旋转过程中,绕组会因切割磁力线而产生反向电动势,而反向电动势与转子相位有直接的对应关系,因此可以利用反电势信息来监测转子的相位。
控制电路通过电流传感器采集电机三相绕组电流,由于同一时刻三相绕组电流之和为零,也可采集两相绕组电流,通过计算得到第三相电流。绕组的等效直流电阻与电机参数有关,是已知量,驱动电压是控制电路给出的已知量,因此根据公式(1)可在线计算得到三相绕组的反向电动势。
检测到三相电机绕组的反向电动势A、B、C,根据这三相反电势的过零点以及极大值点可以将一个电角度周期以30度为步长分为12个区间。每个区间分界处是根据三相反电势的过零点和正负极大值点确定的。但三相的极大值处是信号变化较为缓慢的过程,直接检测准确度相对较低。由于三相绕组在空间是对称安置的,因此三相反电势的正负极大值点与另外两相反电势之差的过零点位置相同,并且该另外两相反电势之差得到的过零过程比该相反电势本身的过零过程变化快,便于检测。因此在具体实现过程中某相的正负极极大值点可通过另外两相的差得到。最终在一个电周期内通过反电势可以得到12个准确的转子位置的特征点,特征点之间相差30度。
(3)基于扩展卡尔曼滤波的转子位置估计
当启动状态完成后,电机就切换到基于反电势无位置算法的速度—电流双闭环矢量控制状态。由电流传感器采样定子电流,一方面通过转子坐标系(dq坐标系)dq 变换输入到电流控制器,另一方面输入到速度-位置观测器中由反电势计算得到转子位置和速度,再分别输入到坐标变换和速度控制器中,构成空间矢量双闭环控制系统。
永磁同步电机在α、β坐标系下的电压方程为:
状态方程为:
其中,Rs为电机相电阻,Ls为电机相电感,λr为电机转子永磁体磁链,考虑到电机的惯性时间常数比系统的采样周期大很多,可以认为
选取状态变量、输入矢量、输出矢量为:x=[iα iβ ωr θe]T,u=[uα uβ]T, y=[iαiβ]T
设系统的采样周期为Tc,该周期可根据实际的电机控制系统载频选取,对上述公式数字化处理得到:
xk=f(xk-1,uk-1,wk-1)=Akxk-1+BTcuk-1+wk-1
yk(xk,vk)=Hkxk+vk
vk是测量噪声,wk是过程噪声。
状态转移矩阵和测量矩阵分别为:
由上述方程可得到扩展卡尔曼器滤波器。
扩展卡尔曼滤波器的输入为经过Clark变换得到的iα、iβ和经过Park逆变换得到的uα、uβ,输出为永磁同步电机的转子位置角和转速。扩展卡尔曼滤波器输出经过闭环反馈送到PI控制器,输出电机的控制电压,通过空间矢量控制三相逆变器驱动永磁同步电机运转。在电机运转过程中,系统每检测到一个特征位置,就利用该特征位置点对应的角度信息对卡尔曼滤波估计的角度信息进行校正,防止卡尔曼滤波的发散,同时也保证了转子角度估计的准确。

Claims (1)

1.一种无位置传感器永磁同步电机转子位置检测方法,其特征在于步骤如下:
步骤1:当电机从速度开环、电流闭环的启动状态切换到速度、电流双闭环的运行状态,通过电流传感器采集电机三相绕组电流,计算电机的三相反向电动势U、V、W分别为:
式中,DA、DB、DC分别是控制三相桥臂的PWM信号的占空比,Ud是施加到绕组上驱动电压,iA、iB、iC是相电流,Ro是绕组的等效直流电阻;
步骤2:通过三相反电动势A、B、C、A-B、B-C、C-A与电角度轴的交点将电角度轴分割为30度等长的12个区间,在一个电周期内得到12个准确的转子位置的特征点,对处于特征点之间的转子位置信息采用扩展卡尔曼滤波方法估计;
步骤3:由电流传感器采样定子电流,一方面通过dq变换输入到电流控制器,另一方面输入到速度-位置观测器中由反电势计算得到转子位置和速度,再分别输入到坐标变换和速度控制器中,构成空间矢量双闭环控制系统;
永磁同步电机在α、β坐标系下的电压方程为:
式中,iα、iβ为电流在α、β轴上的分量,ωe是转子旋转的电角速度、θe是转子旋转的电角度;
状态方程为:
其中,Rs为电机相电阻,Ls为电机相电感,λr为电机转子永磁体磁链,Np是电机的极对数,ωr是电机转子的角速度;
步骤4:选取状态变量、输入矢量、输出矢量为:x=[iα iβ ωr θe]T,u=[uα uβ]T,y=[iαiβ]T
设系统的采样周期为Tc,对公式(2)和(3)进行线性化和离散化得到:
其中,vk是测量噪声,wk是过程噪声,状态转移矩阵和测量矩阵分别为:
步骤5:扩展卡尔曼滤波器的输入为经过Clark变换得到的iα、iβ和经过Park逆变换得到的uα、uβ,输出为永磁同步电机的转子位置角和转速。
CN201710997952.9A 2017-10-24 2017-10-24 一种无位置传感器永磁同步电机转子位置检测方法 Active CN107769636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710997952.9A CN107769636B (zh) 2017-10-24 2017-10-24 一种无位置传感器永磁同步电机转子位置检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710997952.9A CN107769636B (zh) 2017-10-24 2017-10-24 一种无位置传感器永磁同步电机转子位置检测方法

Publications (2)

Publication Number Publication Date
CN107769636A CN107769636A (zh) 2018-03-06
CN107769636B true CN107769636B (zh) 2019-11-22

Family

ID=61269127

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710997952.9A Active CN107769636B (zh) 2017-10-24 2017-10-24 一种无位置传感器永磁同步电机转子位置检测方法

Country Status (1)

Country Link
CN (1) CN107769636B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108599655B (zh) * 2018-03-21 2020-10-27 泉州装备制造研究所 基于权重控制永磁同步电机无速度传感器的转速估计方法
CN110030916A (zh) * 2019-04-18 2019-07-19 安徽拓信电气科技有限公司 一种初级永磁式直线位置传感器
FR3097090B1 (fr) * 2019-06-04 2021-05-21 Renault Sas Procédé d’estimation du couple électromagnétique d’une machine électrique synchrone
DE102019211800B4 (de) 2019-08-06 2022-12-15 Conti Temic Microelectronic Gmbh Verfahren und Vorrichtung zum Ermitteln der Drehzahl und des Drehwinkels einer Motorwelle eines mechanisch kommutierten Gleichstrommotors
CN111009992B (zh) * 2019-12-05 2021-05-28 北京动力机械研究所 一种高、低电压兼容内置式无刷直流起发电系统
CN111638453A (zh) * 2020-06-15 2020-09-08 哈尔滨理工大学 一种伺服同步电机磁场旋转位置和速度的检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176653A (zh) * 2011-01-19 2011-09-07 哈尔滨工业大学 带指数渐消因子的卡尔曼滤波器的感应电机转速观测方法
CN102611381A (zh) * 2012-03-12 2012-07-25 浙江工业大学 永磁同步电机直接转矩控制系统
CN103684178A (zh) * 2013-12-17 2014-03-26 清华大学 一种永磁同步电机转速滤波装置和滤波方法
CN107276479A (zh) * 2017-07-28 2017-10-20 北京控制工程研究所 一种两相正交绕组永磁同步电机转速确定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7276877B2 (en) * 2003-07-10 2007-10-02 Honeywell International Inc. Sensorless control method and apparatus for a motor drive system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176653A (zh) * 2011-01-19 2011-09-07 哈尔滨工业大学 带指数渐消因子的卡尔曼滤波器的感应电机转速观测方法
CN102611381A (zh) * 2012-03-12 2012-07-25 浙江工业大学 永磁同步电机直接转矩控制系统
CN103684178A (zh) * 2013-12-17 2014-03-26 清华大学 一种永磁同步电机转速滤波装置和滤波方法
CN107276479A (zh) * 2017-07-28 2017-10-20 北京控制工程研究所 一种两相正交绕组永磁同步电机转速确定方法

Also Published As

Publication number Publication date
CN107769636A (zh) 2018-03-06

Similar Documents

Publication Publication Date Title
CN107769636B (zh) 一种无位置传感器永磁同步电机转子位置检测方法
CN102946227B (zh) 凸极式永磁无刷直流电机电磁转矩观测方法及装置
CN105391364A (zh) 一种无刷直流电机无位置传感器控制系统及控制方法
CN104158462B (zh) 一种无位置传感器的永磁同步电机初始位置检测方法
CN108418480A (zh) 一种无传感永磁同步电机的启动方法
CN103414427A (zh) 无刷直流电机控制方法
CN107302328B (zh) 基于低频电压注入法的永磁同步电机无位置传感器的控制方法
CN105680742A (zh) 一种无刷直流电机无位置传感器转子位置识别系统及方法
EP2719071A1 (en) Electric drive unit
CN107425774B (zh) 永磁同步电机直轴电感辨识方法及装置
CN108900127A (zh) 考虑交叉耦合效应的ipmsm低速段无位置传感器控制方法
CN108270377A (zh) 一种电机参数测量方法及其装置
CN104485868A (zh) 表贴式永磁同步电机电流预测控制方法
CN207780217U (zh) 一种旋转变压器的零位角测试装置
CN104919696B (zh) 用于确定多相电动机的转子位置的设备
CN104242748A (zh) 一种用于开关磁阻电机系统的转子位置估计方法及装置
CN105811843A (zh) 一种永磁同步电机任意转速快速启动方法
CN102170262B (zh) 一种直驱永磁同步风电机组无速度传感器控制方法
CN110212832A (zh) 用于异步电机的磁链观测方法、系统及介质
Song et al. Sensorless control of surface permanent magnet synchronous motor using a new method
CN104135198A (zh) 永磁同步马达的转子角度估测方法
CN104767451A (zh) 电梯门机无位置传感器电机转子初始位置的检测方法
CN109510525B (zh) 一种永磁同步电机初始状态检测方法
CN101977000A (zh) 电励磁同步电动机的转子位置与速度测量方法及控制装置
Stănică et al. A brief review of sensorless motors position control

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