CN104091060B - 一种分段式Halbach阵列永磁电机磁场计算方法 - Google Patents
一种分段式Halbach阵列永磁电机磁场计算方法 Download PDFInfo
- Publication number
- CN104091060B CN104091060B CN201410309521.5A CN201410309521A CN104091060B CN 104091060 B CN104091060 B CN 104091060B CN 201410309521 A CN201410309521 A CN 201410309521A CN 104091060 B CN104091060 B CN 104091060B
- Authority
- CN
- China
- Prior art keywords
- theta
- permanent magnet
- alpha
- halbach array
- region
- 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
Landscapes
- Permanent Field Magnets Of Synchronous Machinery (AREA)
Abstract
本发明涉及一种分段式Halbach阵列永磁电机磁场计算方法,包括:确定求解区域;针对不同的求解区域,分别建立拉普拉斯方程或泊松方程;计及分段式Halbach阵列磁块之间的间隙,求解永磁体剩余磁化强度各次谐波下的径向及切向分量幅值;对所建立的拉普拉斯方程及泊松方程进行求解,得到求解区域内的标量磁位的表达式,进而得到各个区域磁密的径向分量及切向分量。本发明可以精确地求解具有任意每极磁块数及极对数的内外转子分段式Halbach阵列永磁电机的磁场。
Description
技术领域
本发明涉及一种分段式Halbach阵列永磁电机。
背景技术
Halbach阵列永磁电机的Halbach阵列一般有两种形式,即整环式Halbach阵列及分段式Halbach阵列。但由于整环式Halbach阵列需要特殊的充磁装置且加工较为困难,因此在工程上常用的是分段式Halbach阵列,即将预先充好磁的磁块按照特定的顺序进行组合排列得到所需的Halbach阵列,磁块之间的磁力会使得磁块之间不可避免的存在间隙,这会对气隙磁场造成影响,进而对电机的反电动势、电磁转矩等电磁性能造成影响,因此准确计算分段式Halbach阵列永磁电机的气隙磁场具有重要的意义。
对于气隙磁场的计算一般采用解析法和有限元法,其中有限元法计算结果准确,但是其所耗用的计算时间长,占用的计算资源多,这些因素对磁场分析造成了诸多不便。而相比于有限元法,解析法计算速度快,占用的计算资源少,且随着解析法的发展其计算精度也在逐步提高,同时利用解析法可以方便地分析电机参数对气隙磁场的影响。
近年来针对分段式Halbach阵列永磁电机磁场解析计算方面的方法越来越多,精确度也越来越高,这为分段式Halbach阵列永磁电机气隙磁场的准确计算奠定了基础。
发明内容
本发明的目的是提供一种能够精确计算分段式Halbach阵列永磁电机磁场的方法。本发明的技术方案如下:
一种分段式Halbach阵列永磁电机磁场计算方法,包括下列步骤:
第一步:确定求解区域:建立分段式Halbach阵列永磁电机的物理模型,在电机的物理模型中从里到外定义四个圆周位置,其半径分别为Ri、Rmi、Rmo及Ro,利用极坐标系下的半径坐标r来说明各个区域,其中,r<Ri及r>Ro的区域由铁磁材料组成,Ri<r<Rmi及Rmo<r<Ro的区域为空气,若电机为内转子电机,永磁体内径为Rmi,永磁体外径为Rmo,定子内径为Ro,此时Rmo<r<Ro的区域为有效气隙区域,Rmi<r<Rmo为Halbach阵列永磁体区域,定子铁心由铁磁材料组成,当Ri=Rmi时,转子铁心由铁磁材料组成,当Ri→0时,转子铁心由非铁磁材料组成;若电机为外转子电机,定子外径为Ri,永磁体内径为Rmi,永磁体外径为Rmo,此时Ri<r<Rmi的区域为有效气隙区域,Rmi<r<Rmo为Halbach阵列永磁体区域,定子铁心由铁磁材料组成,同时当Ro=Rmo时,转子铁心由铁磁材料组成,当Ro→∞时,转子铁心由非铁磁材料组成;
第二步:针对不同的求解区域,分别建立拉普拉斯方程或泊松方程:在有效气隙区域建立拉普拉斯方程在分段式Halbach阵列永磁体区域建立泊松方程其中,为有效气隙区域的标量磁位,为分段式Halbach阵列永磁体区域的标量磁位,M为剩余磁化强度矢量,μr为永磁体的相对磁导率,方程中的r代表极坐标系下的半径坐标,θ代表极坐标系下的角坐标;
第三步:计及分段式Halbach阵列磁块之间的间隙,求解永磁体剩余磁化强度各次谐波下的径向及切向分量幅值Mrn及Mθn,其中下标rn及θn分别代表n次谐波径向分量及n次谐波切向分量;
第四步:根据空气与铁磁材料和空气与永磁体交界面处的边界条件,即在空气与铁磁材料的交界面处,磁场强度H只有径向分量而切向分量为0;在空气与永磁体的交界面处,磁场强度H的切向分量相等;在空气与永磁体的交界面处,磁感应强度B的径向分量相等,结合第三步得到的剩余磁化强度n次谐波下的径向及切向分量的幅值Mrn及Mθn,对第二步所建立的拉普拉斯方程及泊松方程进行求解,得到求解区域内的标量磁位的表达式,进而得到各个区域磁密的径向分量及切向分量。
其中的第三步可采用下面的方法:
(1)计及分段式Halbach阵列磁块之间的间隙,将磁块的宽度d2与相邻磁块之间的距离d1的比值定义为相对极弧系数αpr=d2/d1,则相邻磁块之间的间隙为其中,p为极对数,l为每极磁块数,Rmi为永磁体内径,Rmo为永磁体外径;
(2)确定分段式Halbach阵列中各永磁体块的充磁方向,第i(i=1,2,…,2pl)块永磁体块的剩余磁化强度矢量与θ=0之间的夹角为θm,i=(1±p)θi,其中,θi为第i块永磁体中心线与θ=0之间的夹角,其值为
(3)永磁体剩余磁化强度的径向及切向分量是周期性的,其周期为2π/p,计及磁块间隙,写出一个周期内剩余磁化强度径向分量Mr及切向分量Mθ的表达式,其中M为剩余磁化强度矢量的幅值;
(4)对得到的剩余磁化强度径向及切向分量的周期性分段函数进行傅里叶分解得到表达式为
其中n为谐波阶次,Mrn及Mθn分别为剩余磁化强度各次谐波下的径向及切向分量的幅值,通过计算得到Mrn及Mθn的具体表达式如下
若np≠1,且n=1,3,5,7,…,则有
若np≠1,且n=2,4,6,8,…,则有
式中,n为谐波阶次;
“±”和“”中的上层符号表示外转子,下层符号表示内转子;h1n、h2n为待定系数,其表达式为
若np=1,对于外转子电机,则有
若np=1,对于内转子电机,则有
本发明具有如下的突出的有益效果:
1、本发明对气隙磁场进行计算时考虑了分段式Halbach阵列磁块之间的间隙,能够得到更为准确的气隙磁场分布,同时,可以对具有任意每极磁块数及极对数的内外转子分段式Halbach阵列永磁电机的磁场进行分析;
2、利用本发明的计算方法可以分析磁块间隙及计及磁块间隙时每极磁块数及极对数对气隙磁密基波幅值和波形畸变率的影响,为分段式Halbach阵列的组合拼装提供参考。
附图说明
图1 电机求解区域模型图。
图2 分段式Halbach阵列局部示意图。
图3 解析法与有限元法计算内转子Halbach阵列永磁电机气隙磁场结果对比图,其中(a)代表气隙磁密的径向分量对比图,(b)代表气隙磁密的切向分量对比图。
具体实施方式
下面结合附图和实施例对本发明进行说明。
第一步:确定求解区域;
由图1可以看到,求解区域被分为三个,区域I及区域III为空气域,若电机为内转子,则区域I为有效气隙区域,区域III根据转子铁心材料属性进行调节,当Ri=Rmi时,转子铁心由铁磁材料组成,当Ri→0时,转子铁心由非铁磁材料组成,若电机为外转子,则区域III为有效气隙区域,区域I根据转子铁心材料属性进行调节,当Ro=Rmo时,转子铁心由铁磁材料组成,当Ro→∞时,转子铁心由非铁磁材料组成;区域II为分段式Halbach阵列永磁体区域。
第二步:针对不同的求解区域,分别建立拉普拉斯方程或泊松方程;
在区域I及区域III建立拉普拉斯方程
在区域II建立泊松方程
式中:—区域I、II及III中的标量磁位;
M—永磁体剩余磁化强度矢量;
μr—永磁体的相对磁导率;
式中的r、θ代表极坐标系中的半径坐标和角坐标。
第三步:计及分段式Halbach阵列磁块之间的间隙,求解永磁体剩余磁化强度各次谐波下的径向及切向分量幅值Mrn及Mθn;
在分段式Halbach阵列永磁电机中,磁块之间的距离为d1,磁块的宽度为d2,如图2所示,将分段式Halbach阵列中磁块宽度与磁块之间距离的比值定义为相对极弧极数αpr,即
αpr=d2/d1 (4)
确定各块永磁体的充磁方向,分段式Halbach阵列永磁电机中每一块永磁体都有其特定的充磁方向,如图1所示。其中θm,i为第i(i=1,2,…,2pl)块永磁体块的剩余磁化强度矢量与θ=0之间的夹角。
θm,i=(1±p)θi (5)
式中,p—极对数;
l—每极磁块数;
+—外转子;
-—内转子;
θi—第i块永磁体中心线与θ=0之间的夹角。
(i=1,2,3,…,2pl) (6)
根据各块永磁体的充磁方向及各个位置相对于初始位置的夹角,得到一个周期内各个位置剩余磁化强度的径向及切向分量表达式
Mr与Mθ的周期为2π/p,在一个周期内剩余磁化强度径向及切向分量Mr与Mθ的表达式如下式(7)、(8)所示
式中,M—剩余磁化强度的幅值,M=Br/μ0,其中Br为永磁体的剩余磁密,μ0为空气的磁导率。
对式(7)及式(8)中得到的剩余磁化强度的径向及切向分量的表达式进行傅里叶分解,得到表达式如式(9)所示。通过计算得到剩余磁化强度各次谐波下的径向及切向分量幅值的表达式如式(10)~(13)所示
式中,Mrn及Mθn为剩余磁化强度n次谐波下的径向及切向分量的幅值。
通过计算
若np≠1,且n=1,3,5,7,…,则有
若np≠1,且n=2,4,6,8,…,则有
式中,n—谐波阶次;
“±”和“”中的上层符号表示外转子,下层符号表示内转子;h1n、h2n为待定系数,其表达式为
若np=1,对于外转子电机,则有
式中,
若np=1,对于内转子电机,则有
第四步:建立边界条件,对第二步所建立的拉普拉斯方程及泊松方程进行求解,得到三个求解区域内的标量磁位的表达式,进而得到各个区域中磁密的径向分量及切向分量;
边界条件的确立
在空气与铁磁材料的交界面处,磁场强度H只有径向分量而切向分量为0;
在空气与永磁体的交界面处,磁场强度H的切向分量相等;
在空气与永磁体的交界面处,磁感应强度B的径向分量相等。
各个区域磁密径向分量及切向分量表达式
当np≠1时
区域I
区域II
区域III
式中,b1n、g1n、g2n、b3n、Mn为待定系数,其表达式由下式列出
当np=1时
区域I
区域II
区域III
式中,b1、g1、g2、b3为待定系数,其表达式由下式列出
以一台内转子分段式Halbach阵列永磁同步电机为例来介绍提出的磁场计算方法,电机的参数如表1所示。
表1内转子电机参数
确定求解区域
电机为内转子结构,因此由图1可以看出,区域I为有效气隙区域,区域II为分段式Halbach永磁体阵列区域,因为Ri=Rmi,因此区域III为由铁磁材料构成的转子铁心区域。
针对不同求解区域,在极坐标系下,分别建立拉普拉斯方程或泊松方程
在区域I建立拉普拉斯方程
在区域II建立泊松方程
式中:—区域I、II中的标量磁位;
M—永磁体剩余磁化强度矢量;
μr—永磁体的相对磁导率。
永磁体剩余磁化强度的计算
确定各块永磁体的充磁方向
分段式Halbach阵列永磁电机中每一块永磁体都有其特定的充磁方向,其中第i(i=1,2,…,20)块永磁体块的剩余磁化强度矢量与与θ=0之间的夹角如式(28)所示。
θm,i=-θi (28)
θi—第i块永磁体中心线与θ=0之间的夹角,其表达式为
(i=1,2,3,…,20) (29)
Mr与Mθ的周期为π,在一个周期内Mr与Mθ的表达式如下式所示
式中,M—剩余磁化强度的幅值,M=Br/μ0。
对得到的剩余磁化强度的径向及切向分量进行傅里叶分解得到
式中,Mrn及Mθn为剩余磁化强度径向及切向分量n次谐波的幅值。
通过计算
若n=1,3,5,7,…,则有
若n=2,4,6,8,…,则有
式中,n—谐波阶次;
“±”和“”中的上层符号表示外转子,下层符号表示内转子;h1n、h2n为待定系数,其表达式为
边界条件的确立
在空气与铁磁材料的交界面处,磁场强度H只有径向分量而切向分量为0;
在空气与永磁体的交界面处,磁场强度H的切向分量相等;
在空气与永磁体的交界面处,磁感应强度B的径向分量相等。
各个区域磁密径向分量及切向分量表达式
区域I
区域II
式中,b1n、g1n、g2n、Mn为待定系数,其表达式由下式列出
本发明提出的磁场计算方法的正确性的验证
设定相对极弧系数αpr为0.95,建立内转子分段式Halbach阵列永磁同步电机的有限元模型,将由本发明的计算方法得到的结果与由有限元模型计算得到的结果进行对比,如图3所示。由图中可以看出,计算结果与有限元模型计算结果相一致,从而验证了本发明的计算方法的正确性。
Claims (1)
1.一种分段式Halbach阵列永磁电机磁场计算方法,包括下列步骤:
第一步:确定求解区域:建立分段式Halbach阵列永磁电机的物理模型,在电机的物理模型中从里到外定义四个圆周位置,其半径分别为Ri、Rmi、Rmo及Ro,利用极坐标系下的半径坐标r来说明各个区域,其中,r<Ri及r>Ro的区域由铁磁材料组成,Ri<r<Rmi及Rmo<r<Ro的区域为空气,若电机为内转子电机,永磁体内径为Rmi,永磁体外径为Rmo,定子内径为Ro,此时Rmo<r<Ro的区域为有效气隙区域,Rmi<r<Rmo为Halbach阵列永磁体区域,定子铁心由铁磁材料组成,当Ri=Rmi时,转子铁心由铁磁材料组成,当Ri→0时,转子铁心由非铁磁材料组成;若电机为外转子电机,定子外径为Ri,永磁体内径为Rmi,永磁体外径为Rmo,此时Ri<r<Rmi的区域为有效气隙区域,Rmi<r<Rmo为Halbach阵列永磁体区域,定子铁心由铁磁材料组成,同时当Ro=Rmo时,转子铁心由铁磁材料组成,当Ro→∞时,转子铁心由非铁磁材料组成;
第二步:针对不同的求解区域,分别建立拉普拉斯方程或泊松方程:在有效气隙区域建立拉普拉斯方程在分段式Halbach阵列永磁体区域建立泊松方程其中,为有效气隙区域的标量磁位,为分段式Halbach阵列永磁体区域的标量磁位,M为剩余磁化强度矢量,ur为永磁体的相对磁导率,方程中的r代表极坐标系下的半径坐标,θ代表极坐标系下的角坐标;
第三步:计及分段式Halbach阵列磁块之间的间隙,求解永磁体剩余磁化强度各次谐波下的径向及切向分量幅值Mrn及Mθn,其中下标rn及θn分别代表n次谐波径向分量及n次谐波切向分量,方法如下:
(1)计及分段式Halbach阵列磁块之间的间隙,将磁块的宽度d2与相邻磁块之间的距离d1的比值定义为相对极弧系数αpr=d2/d1,则相邻磁块之间的间隙为其中,p为极对数,l为每极磁块数,Rmi为永磁体内径,Rmo为永磁体外径;
(2)确定分段式Halbach阵列中各永磁体块的充磁方向,第i块永磁体块的剩余磁化强度矢量与θ=0之间的夹角为θm,i=(1±p)θi,其中,i=1,2,…,2pl,θi为第i块永磁体中心线与θ=0之间的夹角,其值为
(3)永磁体剩余磁化强度的径向及切向分量是周期性的,其周期为2π/p,计及磁块间隙,写出一个周期内剩余磁化强度径向分量Mr及切向分量Mθ的表达式,其中M为剩余磁化强度矢量的幅值;
(4)对得到的剩余磁化强度径向及切向分量的周期性分段函数进行傅里叶分解得到表达式为
其中n为谐波阶次,Mrn及Mθn分别为剩余磁化强度各次谐波下的径向及切向分量的幅值,通过计算得到Mrn及Mθn的具体表达式如下
若np≠1,且n=1,3,5,7,…,则有
若np≠1,且n=2,4,6,8,…,则有
式中,n为谐波阶次;
“±”和中的上层符号表示外转子,下层符号表示内转子;h1n、h2n为待定系数,其表达式为
若np=1,对于外转子电机,则有
若np=1,对于内转子电机,则有
第四步:根据空气与铁磁材料和空气与永磁体交界面处的边界条件,即在空气与铁磁材料的交界面处,磁场强度H只有径向分量而切向分量为0;在空气与永磁体的交界面处,磁场强度H的切向分量相等;在空气与永磁体的交界面处,磁感应强度B的径向分量相等,结合第三步得到的剩余磁化强度n次谐波下的径向及切向分量的幅值Mrn及Mθn,对第二步所建立的拉普拉斯方程及泊松方程进行求解,得到求解区域内的标量磁位的表达式,进而得到各个区域磁密的径向分量及切向分量的表达式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410309521.5A CN104091060B (zh) | 2014-06-30 | 2014-06-30 | 一种分段式Halbach阵列永磁电机磁场计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410309521.5A CN104091060B (zh) | 2014-06-30 | 2014-06-30 | 一种分段式Halbach阵列永磁电机磁场计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104091060A CN104091060A (zh) | 2014-10-08 |
CN104091060B true CN104091060B (zh) | 2017-01-18 |
Family
ID=51638776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410309521.5A Active CN104091060B (zh) | 2014-06-30 | 2014-06-30 | 一种分段式Halbach阵列永磁电机磁场计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104091060B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111627645A (zh) * | 2020-06-01 | 2020-09-04 | 北京卫星环境工程研究所 | 利用Halbach永磁阵列减小铁磁性材料磁性的方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106229106B (zh) * | 2016-09-12 | 2017-12-01 | 中国科学院苏州生物医学工程技术研究所 | 一种奇数层的Halbach阵列永磁体装置的设计与制作方法 |
DE102016219395A1 (de) * | 2016-10-06 | 2018-04-12 | Brose Fahrzeugteile GmbH & Co. Kommanditgesellschaft, Würzburg | Segmentmagnet und Permanentmagnetmotor mit Segmentmagneten |
CN108052764A (zh) * | 2017-12-26 | 2018-05-18 | 西安交通大学 | 一种环形Halbach永磁阵列的磁场建模方法 |
CN108363858B (zh) * | 2018-02-05 | 2021-04-27 | 东南大学 | 基于极坐标转换的圆筒直线振荡电机磁场确定方法 |
CN109307847B (zh) * | 2018-11-06 | 2021-01-05 | 哈尔滨工程大学 | 一种磁体的二阶欧拉反演法 |
CN109600006B (zh) * | 2018-11-30 | 2020-12-08 | 浙江大学 | 一种用于表贴式永磁电机电磁设计的求解方法 |
CN109948298B (zh) * | 2019-04-10 | 2023-07-18 | 天津工业大学 | 采用不均匀分段Halbach阵列的永磁电机磁场计算方法 |
CN110472336B (zh) * | 2019-08-15 | 2022-03-29 | 哈尔滨理工大学 | 内嵌式永磁同步电机的矩形永磁体等效方法及电磁性能解析方法 |
CN110492696B (zh) * | 2019-09-05 | 2021-01-19 | 合肥工业大学 | 优化Halbach阵列永磁电机裂比和极弧比的方法 |
-
2014
- 2014-06-30 CN CN201410309521.5A patent/CN104091060B/zh active Active
Non-Patent Citations (2)
Title |
---|
分段式Halbach阵列永磁同步电机磁场分析及稳健性优化设计;乔照威;《中国博士学位论文全文数据库》;20130531;第2013年卷(第5期);第二章 * |
基于Halbach阵列的永磁球形电动机磁场;夏长亮 等;《电工技术学报》;20070731;第22卷(第7期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111627645A (zh) * | 2020-06-01 | 2020-09-04 | 北京卫星环境工程研究所 | 利用Halbach永磁阵列减小铁磁性材料磁性的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104091060A (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104091060B (zh) | 一种分段式Halbach阵列永磁电机磁场计算方法 | |
CN109948298B (zh) | 采用不均匀分段Halbach阵列的永磁电机磁场计算方法 | |
CN104158458B (zh) | 一种具有斜槽结构的表贴式永磁电机磁场解析计算方法 | |
CN105005692A (zh) | 一种基于解析法的永磁电机磁场分析与转矩计算方法 | |
Jian et al. | Analytical calculation of magnetic field distribution in coaxial magnetic gears | |
Azzouzi et al. | Quasi-3-D analytical modeling of the magnetic field of an axial flux permanent-magnet synchronous machine | |
CN107612256A (zh) | 一种磁极分段型表贴式永磁同步电机的优化设计方法 | |
CN103853891A (zh) | 一种基于有限元分析的变参数永磁同步电动机建模方法 | |
CN103984864B (zh) | 一种含有不等厚磁极的表贴式永磁电机磁场优化方法 | |
CN108875168B (zh) | 一种考虑饱和的开关磁阻电机磁场解析计算方法 | |
CN112737173B (zh) | 一种分段斜极表贴式永磁电机磁场计算方法 | |
Jiang et al. | Modeling and analysis of spoke-type permanent magnet vernier machine based on equivalent magnetic network method | |
CN109508480A (zh) | 一种构造低频谐波电流计算电机高频电磁力的半解析方法 | |
Hafner et al. | Methods for computation and visualization of magnetic flux lines in 3-D | |
CN106021863A (zh) | 轴向磁通永磁涡流联轴器电磁转矩解析算法 | |
Zhang et al. | A complex-valued rotating magnetic property model and its application to iron core loss calculation | |
Poudel et al. | Deep learning based design methodology for electric machines: Data acquisition, training and optimization | |
Qiu et al. | Analytical analysis of sleeve permeability for output performance of high speed permanent magnet generators driven by micro gas turbines | |
CN105356804A (zh) | 提高同步磁阻电机凸极比的方法 | |
Zhang et al. | Analytical investigation of the magnetic-field distribution in an axial-flux magnetic-geared machine | |
Chen et al. | Investigation of the end leakage flux in fractional slot concentrated-winding surface permanent magnet machines | |
CN106897507A (zh) | 一种基于电感解析模型的内置式永磁同步电机设计方法 | |
Wang et al. | The calculation of eddy current losses density distribution in the permanent magnet of PMSM | |
Guo et al. | Analytical method for predicting air-gap magnetic field distribution of surface mounted permanent magnet machine with non-concentric magnetic poles | |
Kiriyakin et al. | Calculation of Eddy Current Losses in Permanent Magnets of Synchronous Turbine Generators by the Sweep Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |