CN110808703A - 考虑铁损电阻的永磁同步电机转速及转子位置估计方法 - Google Patents

考虑铁损电阻的永磁同步电机转速及转子位置估计方法 Download PDF

Info

Publication number
CN110808703A
CN110808703A CN201911250599.3A CN201911250599A CN110808703A CN 110808703 A CN110808703 A CN 110808703A CN 201911250599 A CN201911250599 A CN 201911250599A CN 110808703 A CN110808703 A CN 110808703A
Authority
CN
China
Prior art keywords
current
electromotive force
coordinate system
back electromotive
permanent magnet
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
CN201911250599.3A
Other languages
English (en)
Other versions
CN110808703B (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201911250599.3A priority Critical patent/CN110808703B/zh
Publication of CN110808703A publication Critical patent/CN110808703A/zh
Application granted granted Critical
Publication of CN110808703B publication Critical patent/CN110808703B/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/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
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

本发明提出一种考虑铁损电阻的永磁同步电机转速及转子位置估计方法,包括以下步骤:首先,根据检测到的绕组状态得到两相静止坐标系下的电流iα、iβ和电压uα、uβ;然后,基于考虑铁损电阻的永磁同步电机等效模型在两相静止坐标系下设计滑模观测器,得到电流观测值
Figure DDA0002308929030000011
以电流真实值iα、iβ和观测值
Figure DDA0002308929030000012
的差值作为滑模面函数,观测出扩展反电动势
Figure DDA0002308929030000013
最后,利用稳态前馈+反馈的控制方法对反电动势进行二次估计,得到当前时刻电机电角速度和电角度的估计值
Figure DDA0002308929030000014
Figure DDA0002308929030000015
本发明在滑模观测器算法的设计中考虑了铁损电阻,提高了估计精度与观测器稳定性;采用稳态前馈+反馈的控制方法对扩展反电动势进行二次估计,减轻了估计结果中的高频抖阵现象。

Description

考虑铁损电阻的永磁同步电机转速及转子位置估计方法
技术领域
本发明属于电机控制技术领域,具体涉及一种考虑铁损电阻的永磁同步电机转速及转子位置估计方法。
背景技术
随着电动汽车的发展,永磁同步电机作为目前常用的车用驱动电机,对其控制算法的研究备受关注。通过永磁同步电机无位置传感器算法估计电机的运行状态,可省去机械传感器的安装,节约电机成本,降低控制系统风险。
目前针对永磁同步电机无位置传感器状态估计算法主要可分为两类:高频注入法和观测器法。其中,高频注入法根据电机的凸极效应,通过外加激励的方法实现对转子位置的估计,该方法需要多个滤波器,具体实现复杂且高频激励的注入会带来噪声。观测器法根据电机的电压信号和电流信号作为观测器的输入,利用控制算法估计转子位置和转速,主要控制算法包括:滑模观测器算法、模型参考自适应算法、扩展卡尔曼滤波器算法。如中国专利公布号CN108574440A,公布日2018-09-25,通过模型参考自适应算法对永磁同步电机的状态进行估计;中国专利公布号CN110138301A,公布日2019-08-16,通过扩展卡尔曼滤波算法对永磁同步电机运行状态进行预测。但模型参考自适应算法和扩展卡尔曼滤波算法在一定程度上都依赖于模型精度,且扩展卡尔曼滤波算法比较复杂,实际应用困难。
滑模观测器算法对模型精度要求不高,具有抗干扰能力强,鲁棒性好等优点。如中国专利公布号CN107579690A,公布日2018-01-12;中国专利公布号CN108599645A,公布日2018-09-28等,通过改进的滑模观测器算法对永磁同步电机状态进行估计,但这些发明在设计滑模观测器时均忽略了电机铁损电阻的影响,铁损电阻的存在会对滑模观测器的稳定性有较大影响。
发明内容
本发明的目的是提供一种考虑铁损电阻的永磁同步电机转速及转子位置估计方法,提高了电机电角速度和电角度的估计精度与观测器稳定性;同时,采用稳态前馈+反馈的控制方法对扩展反电动势进行二次估计,减轻了状态观测结果的高频抖阵现象。
实现本发明目的的技术方案为:一种考虑铁损电阻的永磁同步电机转速及转子位置估计方法,具体包括下述步骤:
步骤1、将检测到的三相电流ia、ib、ic和三相电压ua、ub、uc通过Clark变换得到两相静止坐标系下的电流iα、iβ和电压uα、uβ,根据考虑铁损电阻的永磁同步电机在d、q轴的等效电路,得到考虑铁损电阻的表贴式永磁同步电机在两相静止坐标系下的电压状态方程为:
Figure BDA0002308929010000021
式中,
Figure BDA0002308929010000022
其中,iα、iβ分别为两相静止坐标系下的定子电流分量,Eα、Eβ分别为两相静止坐标系下的扩展反电动势分量,Rm为等效电阻,Ri为定子电阻,Rf为铁损电阻,A为等效铁损电阻,ωe为电角速度,ψm为永磁体磁链,iq为三相电流变换到d-q轴两相旋转坐标系下的q轴电流分量,p为微分算子,Ls为d轴电感,表贴式永磁同步电机的d轴电感与q轴电感相同;
步骤2、根据电压状态方程设计电流状态观测方程,以uα、uβ做为输入,得到两相静止坐标系下的电流观测值
Figure BDA0002308929010000028
,以电流真实值iα、iβ和观测值
Figure BDA0002308929010000029
的差值作为滑模面函数,建立滑模观测器,对所建立的滑模观测器进行稳定性分析,具体步骤为:
(1)、两相静止坐标系下的电流状态观测方程为:
Figure BDA0002308929010000023
其中,为两相静止坐标系下的定子电流观测值,uα、uβ为两相静止坐标系下的电压分量,其作为滑模观测器的控制输入;
Figure BDA0002308929010000025
为两相静止坐标系下扩展反电动势观测值;
Figure BDA0002308929010000026
为滑模观测器的滑模面函数;K为观测器反馈增益。
(2)、定义Lynapunov函数为以下形式:
Figure BDA0002308929010000027
其中,为滑模面函数;
Figure BDA0002308929010000032
为两相静止坐标系下反电动势的观测值与真实值的差值,η为敏感性系数;
为保证滑模观测器稳定,即
Figure BDA0002308929010000033
得到观测器稳定性条件为:
Figure BDA0002308929010000034
步骤3、根据滑模观测器观测出扩展反电动势利用稳态前馈+反馈的控制方法对观测得到的反电动势进行二次估计,以消除反电动势观测结果中的抖动现象,假定在一个控制周期内反电动势变化缓慢,则稳态下
Figure BDA0002308929010000036
的表达式可表示为:
Figure BDA0002308929010000037
利用稳态前馈+反馈的控制方法所得到反电动势估计值的计算公式为:
其中,k、k为积分项系数;分别为采用稳态前馈+反馈控制后的反电动势估计值;
步骤4、根据扩展反电动势的估计值
Figure BDA00023089290100000310
计算求得当前时刻电机电角速度和电角度的估计值
Figure BDA00023089290100000311
Figure BDA00023089290100000312
为:
Figure BDA00023089290100000313
其中,id、iq分别为定子电流经Park变换得到d-q轴下的电流分量。
本发明的有益效果是:
1、根据考虑铁损电阻后的表贴式永磁同步电机电压状态方程,设计滑模观测器时考虑铁损电阻,提高了观测结果的精度与滑模观测器的稳定性;
2、设计电流状态观测方程时,采用连续函数代替开关函数作为反电动势观测结果的补偿,可减小观测结果的抖振现象;
3、采用稳态前馈+反馈的控制方法对反电动势观测值进行二次估计,省去了滤波器,在保证观测结果精度的基础上有效减小了电机电角速度的高频抖动现象。
附图说明
下面结合附图对本发明作进一步的说明:
图1为本发明提出的考虑铁损电阻的永磁同步电机转速及转子位置估计方法流程图;
图2为本发明提出的考虑铁损电阻的永磁同步电机d-q轴等效电路图;
图3为本发明提出的滑模观测器系统框图;
图4为本发明提出的基于稳态前馈+反馈的反电动势控制系统框图。
具体实施方式:
下面结合附图对本发明作进一步说明。
为了解决在永磁同步电机无传感器控制算法中,观测器稳定性易受铁损电阻的影响而变差和观测结果存在高频抖动的问题,本发明提出考虑铁损电阻的永磁同步电机转速及转子位置估计方法,流程图如图1所示,具体步骤如下:
步骤1、将检测到的三相电流ia、ib、ic和三相电压ua、ub、uc通过Clark变换得到两相静止坐标系下的电流iα、iβ和电压uα、uβ,根据图2所示的考虑铁损电阻的永磁同步电机d-q轴等效电路图,列出永磁同步电机在d-q轴下的电压状态方程为:
Figure BDA0002308929010000041
其中,id、iq分别为两相旋转坐标系下(d、q轴下)的定子电流分量,ud、uq分别为两相旋转坐标系下的定子电压分量,Ls为d轴电感,表贴式永磁同步电机的d轴电感与q轴电感相同;
将式(1)经过反Park变换得到两相静止坐标系下的电压方程为:
Figure BDA0002308929010000051
式(1)和式(2)中,
Figure BDA0002308929010000052
其中,iα、iβ分别为两相静止坐标系下的定子电流分量,Eα、Eβ分别为两相静止坐标系下的扩展反电动势分量,Rm为等效电阻,Ri为定子电阻,Rf为铁损电阻,A为等效铁损电阻,ωe为电角速度,ψm为永磁体磁链,p为微分算子。
步骤2、根据电压状态方程设计电流状态观测方程,以uα、uβ做为输入,得到两相静止坐标系下的电流观测值
Figure BDA0002308929010000053
如式(4)所示。
其中,为两相静止坐标系下的定子电流观测值,uα、uβ为两相静止坐标系下的电压分量,其作为滑模观测器的控制输入;
Figure BDA0002308929010000056
为两相静止坐标系下扩展反电动势观测值;为滑模观测器的滑模面函数;K为观测器反馈增益。
以电流真实值iα、iβ和观测值
Figure BDA0002308929010000058
的差值作为滑模面函数,建立滑模观测器如图3所示。滑模面函数的表达式为:
Figure BDA0002308929010000059
滑模面函数的微分值由式(2)和式(4)可得:
Figure BDA0002308929010000061
为使观测器稳定,同时获取反电动势的估计方程,定义李雅普诺夫函数为:
Figure BDA0002308929010000062
其中,
Figure BDA0002308929010000063
为滑模面函数;
Figure BDA0002308929010000064
为两相静止坐标系下反电动势的观测值与真实值的差值,η为敏感性系数。
为保证滑模观测器稳定,即
Figure BDA0002308929010000065
得到观测器稳定性条件为:
步骤3、根据滑模观测器观测出反电动势
Figure BDA0002308929010000067
利用稳态前馈+反馈的控制方法对反电动势进行二次估计,如图4所示。
假定在一个控制周期内反电动势变化缓慢,则稳态下
Figure BDA0002308929010000068
的表达式可表示为:
Figure BDA0002308929010000069
利用稳态前馈+反馈的控制方法所得到反电动势估计值的计算公式为:
Figure BDA00023089290100000610
其中,k、k为积分项系数;
Figure BDA00023089290100000611
分别为采用稳态前馈+反馈控制后的反电动势估计值。
步骤4、根据扩展反电动势值的估计值
Figure BDA0002308929010000071
计算求得当前时刻电机电角速度和电角度的估计值
Figure BDA0002308929010000072
Figure BDA0002308929010000073
计算式如式(11)所示。
Figure BDA0002308929010000074

Claims (1)

1.一种考虑铁损电阻的永磁同步电机转速及转子位置估计方法,其特征在于,包括以下步骤:
步骤1、将检测到的三相电流ia、ib、ic和三相电压ua、ub、uc通过Clark变换得到两相静止坐标系下的电流iα、iβ和电压uα、uβ,根据考虑铁损电阻的永磁同步电机在d、q轴的等效电路,得到考虑铁损电阻的表贴式永磁同步电机在两相静止坐标系下的电压状态方程为:
Figure FDA0002308927000000011
式中,
Figure FDA0002308927000000012
其中,iα、iβ分别为两相静止坐标系下的定子电流分量,Eα、Eβ分别为两相静止坐标系下的扩展反电动势分量,Rm为等效电阻,Ri为定子电阻,Rf为铁损电阻,A为等效铁损电阻,ωe为电角速度,ψm为永磁体磁链,iq为三相电流变换到d-q轴两相旋转坐标系下的q轴电流分量,p为微分算子,Ls为d轴电感,表贴式永磁同步电机的d轴电感与q轴电感相同;
步骤2、根据电压状态方程设计电流状态观测方程,以uα、uβ做为输入,得到两相静止坐标系下的电流观测值以电流真实值iα、iβ和观测值
Figure FDA0002308927000000014
的差值作为滑模面函数,建立滑模观测器,对所建立的滑模观测器进行稳定性分析,具体步骤为:
(1)、两相静止坐标系下的电流状态观测方程为:
Figure FDA0002308927000000015
其中,为两相静止坐标系下的定子电流观测值,uα、uβ为两相静止坐标系下的电压分量,其作为滑模观测器的控制输入;
Figure FDA0002308927000000017
为两相静止坐标系下扩展反电动势观测值;
Figure FDA0002308927000000022
为滑模观测器的滑模面函数;K为观测器反馈增益;
(2)、定义Lynapunov函数为以下形式:
Figure FDA0002308927000000023
其中,
Figure FDA0002308927000000024
为滑模面函数;
Figure FDA0002308927000000025
为两相静止坐标系下反电动势的观测值与真实值的差值,η为敏感性系数;
为保证滑模观测器稳定,即
Figure FDA0002308927000000026
得到观测器稳定性条件为:
Figure FDA0002308927000000027
步骤3、根据滑模观测器观测出扩展反电动势
Figure FDA0002308927000000028
采用稳态前馈+反馈的控制方法对观测得到的反电动势进行二次估计,以消除反电动势观测结果中的抖动现象,假定在一个控制周期内反电动势变化缓慢,则稳态下
Figure FDA0002308927000000029
的表达式可表示为:
Figure FDA00023089270000000210
利用稳态前馈+反馈的控制方法所得到反电动势估计值的计算公式为:
Figure FDA00023089270000000211
其中,k、k为积分项系数;分别为采用稳态前馈+反馈控制后的反电动势估计值;
步骤4、根据扩展反电动势值的估计值
Figure FDA00023089270000000213
计算求得当前时刻电机电角速度和电角度的估计值为:
Figure FDA0002308927000000031
其中,id、iq分别为定子电流经Park变换得到d-q轴下的电流分量。
CN201911250599.3A 2019-12-09 2019-12-09 考虑铁损电阻的永磁同步电机转速及转子位置估计方法 Active CN110808703B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911250599.3A CN110808703B (zh) 2019-12-09 2019-12-09 考虑铁损电阻的永磁同步电机转速及转子位置估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911250599.3A CN110808703B (zh) 2019-12-09 2019-12-09 考虑铁损电阻的永磁同步电机转速及转子位置估计方法

Publications (2)

Publication Number Publication Date
CN110808703A true CN110808703A (zh) 2020-02-18
CN110808703B CN110808703B (zh) 2021-01-22

Family

ID=69492762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911250599.3A Active CN110808703B (zh) 2019-12-09 2019-12-09 考虑铁损电阻的永磁同步电机转速及转子位置估计方法

Country Status (1)

Country Link
CN (1) CN110808703B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112713824A (zh) * 2021-01-04 2021-04-27 潍柴动力股份有限公司 一种永磁同步电机无传感器控制方法及系统
CN113938062A (zh) * 2021-09-30 2022-01-14 南京航空航天大学 一种三相永磁同步电机的扩展反电动势观测方法
CN115037212A (zh) * 2022-06-17 2022-09-09 合肥工业大学 一种基于位置观测的永磁同步电机转子位置估计方法
CN115580192A (zh) * 2022-10-25 2023-01-06 宁波奥克斯电气股份有限公司 电机控制方法、装置、控制器及计算机可读存储介质
CN116620042A (zh) * 2023-06-05 2023-08-22 小米汽车科技有限公司 电机控制参数标定方法、装置及车辆
RU2825800C1 (ru) * 2024-02-06 2024-08-29 Федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский Томский политехнический университет" Устройство получения сигнала оценки частоты вращения ротора и сигнала оценки момента сопротивления на валу синхронного электродвигателя с постоянными магнитами

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100237817A1 (en) * 2009-03-23 2010-09-23 Jingbo Liu Method and Apparatus for Estimating Rotor Position in a Sensorless Synchronous Motor
CN103997270A (zh) * 2014-06-09 2014-08-20 浙江理工大学 一种隐极式永磁同步电机的无传感器矢量控制装置及方法
CN105281631A (zh) * 2014-06-16 2016-01-27 现代自动车株式会社 电动机的无传感器控制方法和使用该方法的系统
CN107579690A (zh) * 2017-08-28 2018-01-12 南京理工大学 一种基于滑模观测的超高速永磁同步电机转速估计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100237817A1 (en) * 2009-03-23 2010-09-23 Jingbo Liu Method and Apparatus for Estimating Rotor Position in a Sensorless Synchronous Motor
CN103997270A (zh) * 2014-06-09 2014-08-20 浙江理工大学 一种隐极式永磁同步电机的无传感器矢量控制装置及方法
CN105281631A (zh) * 2014-06-16 2016-01-27 现代自动车株式会社 电动机的无传感器控制方法和使用该方法的系统
CN107579690A (zh) * 2017-08-28 2018-01-12 南京理工大学 一种基于滑模观测的超高速永磁同步电机转速估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
N. R. ABJADI ET AL: "A Nonlinear Adaptive Controller for Speed Sensorless PMSM Taking the Iron Loss Resistance into Account", 《2005 INTERNATIONAL CONFERENCE ON ELECTRICAL MACHINES AND SYSTEMS》 *
苑婷 等: "考虑铁损的超高速永磁同步电机无传感器控制", 《电力电子技术》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112713824A (zh) * 2021-01-04 2021-04-27 潍柴动力股份有限公司 一种永磁同步电机无传感器控制方法及系统
CN113938062A (zh) * 2021-09-30 2022-01-14 南京航空航天大学 一种三相永磁同步电机的扩展反电动势观测方法
CN113938062B (zh) * 2021-09-30 2023-12-19 南京航空航天大学 一种三相永磁同步电机的扩展反电动势观测方法
CN115037212A (zh) * 2022-06-17 2022-09-09 合肥工业大学 一种基于位置观测的永磁同步电机转子位置估计方法
CN115037212B (zh) * 2022-06-17 2024-03-08 合肥工业大学 一种基于位置观测的永磁同步电机转子位置估计方法
CN115580192A (zh) * 2022-10-25 2023-01-06 宁波奥克斯电气股份有限公司 电机控制方法、装置、控制器及计算机可读存储介质
CN116620042A (zh) * 2023-06-05 2023-08-22 小米汽车科技有限公司 电机控制参数标定方法、装置及车辆
CN116620042B (zh) * 2023-06-05 2024-02-06 小米汽车科技有限公司 电机控制参数标定方法、装置及车辆
RU2825800C1 (ru) * 2024-02-06 2024-08-29 Федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский Томский политехнический университет" Устройство получения сигнала оценки частоты вращения ротора и сигнала оценки момента сопротивления на валу синхронного электродвигателя с постоянными магнитами

Also Published As

Publication number Publication date
CN110808703B (zh) 2021-01-22

Similar Documents

Publication Publication Date Title
CN110808703B (zh) 考虑铁损电阻的永磁同步电机转速及转子位置估计方法
Sul et al. Sensorless control of IPMSM: Past, present, and future
CN104901600B (zh) 宽转速范围内永磁同步电机的无位置传感器控制方法
CN110912485B (zh) 一种考虑铁损电阻的永磁同步电机无位置传感器控制方法
CN105915142B (zh) 一种基于解耦自适应观测器的永磁同步电机转子位置和转速估算方法
CN103414423A (zh) 一种面贴式永磁同步电机无位置传感器直接转矩控制方法
CN110995095B (zh) 无位置传感器的永磁同步电机控制方法及汽车动力系统
CN108900131A (zh) 表贴式永磁同步电机低速无位置传感器矢量控制系统与方法
CN106026817A (zh) 一种基于采用卡尔曼滤波器的滑模观测器的无速度传感器控制系统
CN108347207B (zh) 基于复数pi控制器的永磁同步电机位置和速度估算方法
CN106208864A (zh) 一种基于smo的无速度传感器控制系统
CN114465543B (zh) 一种永磁同步电机无位置传感器控制方法
CN106059424A (zh) 一种基于改进卡尔曼观测器的无速度传感器控制方法
CN112671298B (zh) 改进的用于永磁同步电机控制的pll无感控制算法
CN106026831B (zh) 一种无位置传感器控制下的电流解耦方法
CN112886880A (zh) 三电平永磁同步电机无位置传感器模型预测电流控制方法
CN114598206A (zh) 永磁同步电机宽速域转子位置观测器设计方法
CN106026834A (zh) 一种永磁同步电机的无速度传感器控制方法
CN110649849A (zh) 一种基于新型滑模观测器的磁通切换型永磁直线电机无位置控制方法
CN115276487A (zh) 一种永磁同步电机控制方法
CN108880351A (zh) 永磁同步电机转子位置的估算方法及系统
CN113364375B (zh) 变结构电流调节器的pmsm驱动系统无传感器控制方法
CN104811117A (zh) 一种永磁同步电机转子转速估计的方法
CN118157526A (zh) 一种基于改进型线性超螺旋的永磁同步电机控制方法
CN113037163A (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