CN112949146B - 一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法 - Google Patents
一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法 Download PDFInfo
- Publication number
- CN112949146B CN112949146B CN202110357541.XA CN202110357541A CN112949146B CN 112949146 B CN112949146 B CN 112949146B CN 202110357541 A CN202110357541 A CN 202110357541A CN 112949146 B CN112949146 B CN 112949146B
- Authority
- CN
- China
- Prior art keywords
- air gap
- stator
- rotor
- tooth
- magnetic
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
Abstract
本发明涉及一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,包括以下步骤:建立定子永磁型磁通切换电机的解析模型,获得线性磁路无槽情况下气隙磁通密度;根据转子位置角,确定出定子齿、转子齿重叠结构模式及其对应的定子齿、转子齿之间的重叠角,采用有限元法分析电机铁芯饱和时的磁密分布,根据定、转子的重合角重叠角,估算出定子齿饱和磁阻、转子齿饱和磁阻;将定、转子齿饱和部分的磁阻等效为对应齿前气隙磁阻的增加,同时重新构建出气隙长度等效后的相对磁导函数;根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度。本发明可以快速分析电机的饱和情况下的磁场分布,以解决磁场解析法在分析电机饱和磁场时存在的缺陷。
Description
技术领域
本发明涉及电机设计领域,具体涉及一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法。
背景技术
定子永磁型磁通切换电机(包括机械轴承支撑磁通切换电机和无轴承磁通切换电机)定、转子显著的双凸极性特点,使得该种电机在齿尖处具有聚磁效应,可提高永磁材料的利用率。永磁体采用切向充磁方式且放置在定子上,不仅便于永磁体的安装,而且还有利于永磁体的散热,降低了永磁退磁的风险。在电机的绕组结构上,因具有绕组的一致性和互补性,使得电机的反电动势波形具有较高的正弦性。
当电机带重载运行时,为维持电机的正常运行,需要通入更大的绕组电流,由于铁芯材料B-H曲线具有非线性,当电机铁芯中的磁密超过膝点时,B-H磁化曲线不再线性变化,电机铁芯会发生磁饱和现象。此时,若继续通过增加定子绕组电流则电机的转矩并不再继续增大,相反带来电机的效率下降。为准确分析电机的磁场特点,必须充分考虑磁场的饱和现象。
传统的分析电机磁饱和问题,主要采用有限元法。在有限元模型中可以通过设置铁芯材料的磁化曲线数值,分析饱和情况下的磁场分布,但有限元法存在计算耗时过长等特点。在工程领域中,通常都是基于大量的数据进行类比估计,依赖一定的经验性。而目前电机磁饱和问题分析的方法十分匮乏。
发明内容
有鉴于此,本发明的目的在于提供一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,有效提高电机在传统线性磁路条件下的计算精度,拓宽了解析方法分析电机磁场的应用范围,能够对电机的饱和磁场进行分析。
为实现上述目的,本发明采用如下技术方案:
一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,包括以下步骤:
在线性磁路下,采用解析法,建立定子永磁型磁通切换电机的解析模型,获得考虑定子双边开槽情况下的气隙磁密Bslot;同时,选取磁场中的标量磁位为求解变量,计算定、转子双边开槽情况下的相对磁导函数,并进一步获得线性磁路无槽情况下气隙磁通密度;
根据转子位置角θr,确定出定子齿、转子齿重叠结构模式1~k及其对应的定子齿、转子齿之间的重叠角γi,i=1~k,采用有限元法分析电机铁芯饱和时的磁密分布,根据定、转子的重合角重叠角γi,估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri;
将定、转子齿饱和部分的磁阻等效为对应齿前气隙磁阻的增加,同时重新构建出气隙长度等效后的相对磁导函数;
根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度Bslot_sati;当转子位置角大于36°时,计算下一个周期的非线性磁路下的气隙磁通密度分布。
进一步的,所述估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri,具体如下:
设定、转子齿各自的饱和深度均呈线性变化,定子齿饱和深度gsi、转子齿饱和深度gri与重叠角γi,i=1~k,具有下列关系:
其中,hs、hr为定、转子齿的长度;αs、αr为定、转子齿的宽度;
根据磁导率为μi,计算出第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′,转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′分别为:
其中,Asi为定子齿的截面积、Ari为转子齿截面积。
进一步的,所述根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度Bslot_sati,具体如下:
设第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′相等;转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′相等,计算定子齿饱和段等效气隙长度Δgsi、转子齿饱和段等效气隙长度Δgri:
计算第i个定子齿、转子齿重叠结构模式中总的等效气隙长度gi:
gi=g+Δgsi+Δgri (9)
根据等效前后磁动势守恒原则及第i个定子齿、转子齿重叠结构模式中线性磁路无槽气隙磁密Bslotlessi,第i个定子齿、转子齿重叠结构模式中非线性磁路饱和无槽气隙磁密Bslotless_sati:
根据定子开槽时的气隙磁密可以由无槽气隙磁密与相对磁导函数相乘后获得原理,计算定子双边开槽时第i个定子齿、转子齿重叠结构模式中的磁密为:
Bslot_sati=Bslotless_sati·Λslot_sati (11)
其中,Bslot_sati为非线性磁路开槽下的磁密,Λslot_sati为气隙长度等效后的相对磁导函数。
进一步的,所述气隙长度等效后的相对磁导函数构建,具体如下:
根据获得的等效气隙总长度,重新构建带有饱和补偿的线性电机结构;
基于带有饱和补偿的线性电机结构,利用解析模型,计算第i个定子齿、转子齿重叠结构模式中相对磁导函数Λslot_sati。
进一步的,所述解析模型采用子域模型法或相对磁导函数法。
本发明与现有技术相比具有以下有益效果:
1、本发明有效提高分析计算的效率,缩短电机模型的计算时间。在电机模型尺寸改变时,避开了有限元法反复建模的缺点。
2、本发明有效提高电机在传统线性磁路条件下的计算精度,拓宽了解析方法分析电机磁场的应用范围,能够对电机的饱和磁场进行分析。
附图说明
图1是本发明方法流程图;
图2是本发明一实施例中12/10定子永磁型磁通切换电机横截面图;
图3是本发明一实施例中典型的3相绕组连接示意图;
图4是本发明一实施例中典型6相绕组连接示意图;
图5是本发明一实施例中典型双绕组连接示意图;
图6是本发明一实施例中定、转子重合角相对位置;
图7是本发明一实施例中定、转子齿间的等效气隙长度;
图8是本发明一实施例中饱和补偿以后的定子永磁型磁通切换电机横截面图;
图9是本发明一实施例中B-H非线性条件下的径向气隙磁密分布。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
本实施例中,磁通切换电机典型结构为12/10型,该电机定、转子呈双凸极结构,定子没有绕组时电机横截面结构如图2所示。定子部分由12个U型铁芯组成,并且相邻的两个U型铁芯之间夹着一个沿切向交替充磁的永磁体;转子部分由硅钢片堆叠而成,共有10个齿。电机共有12线圈,可以根据需要,将12个线圈连接成3相绕组(A相~C相)如图3所示,用于流过对称电流产生平稳的电磁转矩;也可以连接成对称6相绕组(A相~F相)如图4所示,同时流过转矩电流分量和悬浮电流分量,在产生平稳电磁转矩输出同时产生转子磁悬浮所需要的悬浮力;还可以在永磁体槽中同时放置三相悬浮绕组(a相~c相),与三相转矩绕组(A相~C相)共同构成无轴承磁通切换电机的双绕组结构如图5所示。
请参照图1,本发明提供一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,包括以下步骤:
在线性磁路下,采用解析法,建立定子永磁型磁通切换电机的解析模型,获得考虑定子双边开槽情况下的气隙磁密Bslot;同时,选取磁场中的标量磁位为求解变量,计算定、转子双边开槽情况下的相对磁导函数,并进一步获得线性磁路中无槽情况下气隙磁通密度;
根据转子位置角θr,确定出定子齿、转子齿重叠结构模式1~k及其对应的定子齿、转子齿之间的重叠角γi,i=1~k,采用有限元法分析电机铁芯饱和时的磁密分布,根据定、转子的重合角重叠角γi,估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri;
优选的,估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri,具体如下:
设定、转子齿各自的饱和深度均呈线性变化,定子齿饱和深度gsi、转子齿饱和深度gri与重叠角γi,i=1~k,具有下列关系:
其中,hs、hr为定、转子齿的长度,αs、αr为定、转子齿的宽度。
根据磁导率为μi,计算出第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′,转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′分别为:
其中,Asi为定子齿的截面积、Ari为转子齿截面积。
将定、转子齿饱和部分的磁阻等效为对应齿前气隙磁阻的增加,同时重新构建出气隙长度等效后的相对磁导函数;根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度Bslot_sati;
优选的,具体如下:设第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′相等;转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′相等,计算定子齿饱和段等效气隙长度Δgsi、转子齿饱和段等效气隙长度Δgri:
计算第i个定子齿、转子齿重叠结构模式中总的等效气隙长度gi:
gi=g+Δgsi+Δgri (9)
根据等效前后磁动势守恒原则及第i个定子齿、转子齿重叠结构模式中线性磁路无槽气隙磁密Bslotlessi,第i个定子齿、转子齿重叠结构模式中非线性磁路饱和无槽气隙磁密Bslotless_sati:
根据定子开槽时的气隙磁密可以由无槽气隙磁密与相对磁导函数相乘后获得原理,计算定子双边开槽时第i个定子齿、转子齿重叠结构模式中的磁密为:
Bslot_sati=Bslotless_sati·Λslot_sati (11)
其中,Bslot_sati为非线性磁路开槽下的磁密,Λslot_sati为气隙长度等效后的相对磁导函数。
当转子位置角大于36°时,计算下一个周期的非线性磁路下的气隙磁通密度分布。
优选的,在本实施例中,磁导率为μi获取具体如下,利用解析模型,例如子域模型法或相对磁导函数法,计算第i个定子齿、转子齿重叠结构模式中线性磁路下的气隙磁通密度B。根据B及铁芯材料B-μi曲线关系,查得磁导率为μi。
优选的,在本实施例中,根据获得的等效气隙总长度,重新构建带有饱和补偿的线性电机结构;基于带有饱和补偿的线性电机结构,利用解析模型,计算第i个定子齿、转子齿重叠结构模式中相对磁导函数Λslot_sati。
优选的,在本实施例中,解析模型采用子域模型法或相对磁导函数法。
通过本发明所提等效气隙长度变化的磁场饱和补偿方法,与有限元法获得的径向气隙磁密分布如图9所示。
从图9可以看出,针对磁路的非线性特点,本方法在磁密饱和程度严重时,补偿的效果明显;本发明补偿方法能够反应出与有限元法相接近的趋势,具有一定的磁饱和补偿效果,验证了本发明磁密饱和分析方法的可行性。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (4)
1.一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,其特征在于,包括以下步骤:
在线性磁路下,采用解析法,建立定子永磁型磁通切换电机的解析模型,获得考虑定子双边开槽情况下的气隙磁密Bslot;同时,选取磁场中的标量磁位为求解变量,计算定、转子双边开槽情况下的相对磁导函数,并获得线性磁路中无槽情况下气隙磁通密度;
根据转子位置角θr,确定出定子齿、转子齿重叠结构模式1~k及其对应的定子齿、转子齿之间的重叠角γi,i=1~k,采用有限元法分析电机铁芯饱和时的磁密分布,根据定、转子的重叠角γi,估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri;
将定、转子齿饱和部分的磁阻等效为对应齿前气隙磁阻的增加,同时重新构建出气隙长度等效后的相对磁导函数;
根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度Bslot_sati;当转子位置角大于36°时,计算下一个周期的非线性磁路下的气隙磁通密度分布;
所述估算出定子齿饱和磁阻Rsi、转子齿饱和磁阻Rri,具体如下:
设定、转子齿各自的饱和深度均呈线性变化,定子齿饱和深度gsi、转子齿饱和深度gri与重叠角γi,i=1~k,具有下列关系:
其中,hs、hr为定、转子齿的长度;αs、αr为定、转子齿的宽度;
根据磁导率为μi,计算出第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′,转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′分别为:
其中,Asi为定子齿的截面积、Ari为转子齿截面积。
2.根据权利要求1所述的一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,其特征在于,所述根据磁动势守恒原则,计算非线性磁路中补偿后的气隙磁通密度Bslot_sati,具体如下:
设第i个定子齿、转子齿重叠结构模式中定子齿饱和磁阻Rsi、定子齿饱和段等效气隙磁阻Rsi′相等;转子齿饱和磁阻Rri、转子齿饱和段等效气隙磁阻Rri′相等,计算定子齿饱和段等效气隙长度Δgsi、转子齿饱和段等效气隙长度Δgri:
计算第i个定子齿、转子齿重叠结构模式中总的等效气隙长度gi:
gi=g+Δgsi+Δgri (9)
根据等效前后磁动势守恒原则及第i个定子齿、转子齿重叠结构模式中线性磁路无槽气隙磁密Bslotlessi,第i个定子齿、转子齿重叠结构模式中非线性磁路饱和无槽气隙磁密Bslotless_sati:
根据定子开槽时的气隙磁密可以由无槽气隙磁密与相对磁导函数相乘后获得原理,计算定子双边开槽时第i个定子齿、转子齿重叠结构模式中的磁密为:
Bslot_sati=Bslotless_sati·Λslot_sati (11)
其中,Bslot_sati为非线性磁路开槽下的磁密,Λslot_sati为气隙长度等效后的相对磁导函数。
3.根据权利要求1所述的一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,其特征在于,所述气隙长度等效后的相对磁导函数构建,具体如下:
根据获得的等效气隙总长度,重新构建带有饱和补偿的线性电机结构;
基于带有饱和补偿的线性电机结构,利用解析模型,计算第i个定子齿、转子齿重叠结构模式中相对磁导函数Λslot_sati。
4.根据权利要求3所述的一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法,其特征在于,所述解析模型采用子域模型法或相对磁导函数法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110357541.XA CN112949146B (zh) | 2021-04-01 | 2021-04-01 | 一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110357541.XA CN112949146B (zh) | 2021-04-01 | 2021-04-01 | 一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112949146A CN112949146A (zh) | 2021-06-11 |
CN112949146B true CN112949146B (zh) | 2022-06-14 |
Family
ID=76232085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110357541.XA Active CN112949146B (zh) | 2021-04-01 | 2021-04-01 | 一种等效气隙长度变化的磁通切换电机磁场饱和补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112949146B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992663A (zh) * | 2017-11-27 | 2018-05-04 | 中国矿业大学 | 一种考虑铁耗的开关磁阻电机动态电路建模方法 |
CN108563912A (zh) * | 2018-05-11 | 2018-09-21 | 北京理工大学 | 一种表贴式永磁同步电机气隙磁场的解析方法 |
CN109617466A (zh) * | 2018-11-09 | 2019-04-12 | 中国矿业大学 | 一种开关磁阻电机复杂电力电子系统的非线性建模方法 |
CN111651914A (zh) * | 2020-05-14 | 2020-09-11 | 福州大学 | 一种无轴承磁通切换电机的磁场解析计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10944302B2 (en) * | 2018-04-09 | 2021-03-09 | Williams International Co., L.L.C. | Permanent-magnet generator incorporating a variable-reluctance stator system |
-
2021
- 2021-04-01 CN CN202110357541.XA patent/CN112949146B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992663A (zh) * | 2017-11-27 | 2018-05-04 | 中国矿业大学 | 一种考虑铁耗的开关磁阻电机动态电路建模方法 |
CN108563912A (zh) * | 2018-05-11 | 2018-09-21 | 北京理工大学 | 一种表贴式永磁同步电机气隙磁场的解析方法 |
CN109617466A (zh) * | 2018-11-09 | 2019-04-12 | 中国矿业大学 | 一种开关磁阻电机复杂电力电子系统的非线性建模方法 |
CN111651914A (zh) * | 2020-05-14 | 2020-09-11 | 福州大学 | 一种无轴承磁通切换电机的磁场解析计算方法 |
Non-Patent Citations (1)
Title |
---|
崔征山 等.基于改进子域模型的双绕组无轴承磁通切换电机磁场解析计算.《电机与控制学报》.2022, * |
Also Published As
Publication number | Publication date |
---|---|
CN112949146A (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhao et al. | Material-efficient permanent-magnet shape for torque pulsation minimization in SPM motors for automotive applications | |
Bonthu et al. | Comparisons of rare-earth and rare-earth-free external rotor permanent magnet assisted synchronous reluctance motors | |
Zhu et al. | Comparative design and analysis of new type of flux-intensifying interior permanent magnet motors with different q-axis rotor flux barriers | |
Bianchi et al. | An analytical approach to design the PM in PMAREL motors robust toward the demagnetization | |
Zhu et al. | Electromagnetic performance comparison of multi-layered interior permanent magnet machines for EV traction applications | |
Sebastian et al. | Analysis of induced EMF waveforms and torque ripple in a brushless permanent magnet machine | |
Hua et al. | Design of flux-switching permanent magnet machine considering the limitation of inverter and flux-weakening capability | |
Li et al. | A novel spoke-type PM machine employing asymmetric modular consequent-pole rotor | |
CN111654165B (zh) | 基于转子磁势模型的内置式永磁同步电机磁场计算方法 | |
De Donato et al. | No-load performance of axial flux permanent magnet machines mounting magnetic wedges | |
Ji et al. | A novel linear permanent-magnet vernier machine with improved force performance | |
Li et al. | Flux-focusing permanent magnet machines with modular consequent-pole rotor | |
Chen et al. | Magnet-frozen-permeability FEA and DC-biased measurement for machine inductance: Application on a variable-flux PM machine | |
Jara et al. | Rotor eddy-current losses reduction in an axial flux permanent-magnet machine | |
Li et al. | Analytical determination of optimal PM-arc ratio of consequent-pole permanent magnet machines | |
Wang et al. | Analysis and optimization of hybrid excitation permanent magnet synchronous generator for stand-alone power system | |
Du et al. | Analysis and design of a fault-tolerant permanent magnet vernier machine with improved power factor | |
Hao et al. | Analysis of cogging torque reduction techniques in axial-field flux-switching permanent-magnet machine | |
CN103560634B (zh) | 电动汽车用内置式永磁同步电机 | |
Liu et al. | Analysis and research on no-load air gap magnetic field and system multiobjective optimization of interior PM motor | |
Zhao et al. | Optimal design and experimental test of a SPM motor with cost-effective magnet utilization to suppress torque pulsations | |
Mahmoudi et al. | An analytical complementary FEA tool for optimizing of axial-flux permanent-magnet machines | |
Wang et al. | A novel spoke-type IPM rotor with hybrid radial and axial flux concentration for reduction of interpolar leakage flux | |
Niu et al. | Eddy-current analysis of double-stator inset-type permanent magnet brushless machines | |
Barcaro et al. | Air-gap flux density distortion and iron losses in anisotropic synchronous motors |
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 |