CN103401504B - 一种永磁同步电机转子初始位置的修正方法 - Google Patents

一种永磁同步电机转子初始位置的修正方法 Download PDF

Info

Publication number
CN103401504B
CN103401504B CN201310339651.9A CN201310339651A CN103401504B CN 103401504 B CN103401504 B CN 103401504B CN 201310339651 A CN201310339651 A CN 201310339651A CN 103401504 B CN103401504 B CN 103401504B
Authority
CN
China
Prior art keywords
vibra
pos
current
rotor
current phasor
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
CN201310339651.9A
Other languages
English (en)
Other versions
CN103401504A (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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
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 Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN201310339651.9A priority Critical patent/CN103401504B/zh
Publication of CN103401504A publication Critical patent/CN103401504A/zh
Application granted granted Critical
Publication of CN103401504B publication Critical patent/CN103401504B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明提供一种永磁同步电机转子初始位置的修正方法,通过施加正负脉冲的电流矢量,使转子高频抖动,通过抖动量的大小,可以判断转子初始位置的偏差大小,通过改变电流矢量的方向,逐步减小偏差值。通过正负脉冲的作用,避免了转子向一个方向运动范围过大,或时间过长,从而不会对下一个电流矢量作用时转子的运动造成影响,避免了判断失误的问题。电流矢量的幅值和正负脉冲的作用时间可调,保证了转子振动的频率和幅值可调,增强了检测过程的灵活性和抗扰性。

Description

一种永磁同步电机转子初始位置的修正方法
技术领域
本发明属于电机控制领域,涉及永磁同步电机转子初始位置的修正方法。
背景技术
永磁同步电机由于功耗低、效率高、调速范围宽、坚固耐用、转矩纹波小等优点,在很多应用领域,尤其是高性能应用场合,如机器人和高精度数控机床以及伺服转台中,得到了越来越广泛的应用。
永磁同步电机采用电子换向,所以转子初始位置的偏差直接影响着电机输出力矩的大小,甚至会造成系统反转,影响系统的稳定运行。
在伺服转台等高精度控制的场合一般需要安装光电编码器、旋转变压器等位置传感器来获得电机转子的准确位置。对于绝对式位置传感器,初次安装时需要准确确定转子初始位置,而对于相对式位置传感器,每次上电时都需要确定转子初始位置。
目前常用的方法有预定位法,即在电机定子中加入固定方向的电压或电流矢量,电机转子将在电磁力的作用下旋转到某一固定位置,从而确定出转子初始位置。但采用此种方法电机会产生较大角度的偏转,这在一些应用场合是不允许的,可能会损坏仪器。并且如果电机带负载,或者电机的摩擦转矩较大,该方法将造成较大的估计误差。
另一类是如中国专利:永磁同步电机转子初始位置角的确定方法(授权公告号:CN101594114B)和中国专利:交流伺服系统中电机转子初始位置的检测方法及装置(公开号CN1148168A)提出了一类检测转子初始位置的方法,通过直接控制定子电流矢量的幅值和相位,以产生不同方向的定子磁场,控制转子的转动,使定子电流产生的磁场方向逐步逼近转子磁场方向。这类方法所施加的电流矢量在每个固定的方向上均为正值或均为负值,然后通过转子的逆时针旋转或顺时针旋转方向来判断下一步电流矢量的方向。此类方法存在如下问题:由于摩擦力的存在,尤其是静摩擦力往往要大于滑动摩擦力,这会造成如果所加电流矢量产生的电磁转矩太小,转子无法运动,或转动方向不明显,易受外界扰动的影响;产生的电磁转矩太大,又会使转子运动幅度过大,运动时间过长,这在一些场合是不安全的,另外也会影响下一个电流矢量作用时转子的运动方向,造成判断失误。
发明内容
本发明的技术解决方案是:克服现有技术的不足,提供一种方便、安全、可靠性高的永磁同步电机转子初始位置的修正方法。
本发明的技术解决方案:在电流矢量的每个方向上,施加具有正负幅值的电流值,正负脉冲的比例为1:1。在此正负脉冲作用下,可以降低摩擦的影响;另外,电机转子不会向一个方向做大幅运动,而是在初始位置附近做小幅度的振动,通过振动的幅度来判断转子的初始位置。
本发明的方法包括以下步骤:一种永磁同步电机转子初始位置的修正方法,包括以下步骤:
步骤一、初始化以下变量值:永磁同步电机d轴电流id_Ref=0,q轴电流的幅值iq_Amp=A0,电流矢量的方向Pos0(k)=0, 其中k=0;
步骤二、依次施加Pos0(k)、Posp(k)和Posn(k)方向的电流矢量,每个电流矢量包括n0个正脉冲和n0个负脉冲,正负脉冲交替施加,每个脉冲的宽度为T0,计算每个电流矢量对应的位置波动:Vibra0(k),Vibrap(k),Vibran(k);
步骤三、判断电流矢量的幅值iq_Amp是否小于限定值iq_Max,并且位置波动Vibra0(k),Vibrap(k)和Vibran(k)是否均小于某个限定值VibraMin,即判断(iq_Amp<iq_Max)&(Vibra0(k)<VibraMin)&(Vibrap(k)<VibraMin)&(Vibran(k)<VibraMin)是否为真;如果是,转向步骤七;如果否,转向步骤四;
步骤四、更新k=k+1,即将k加1,找出Vibra0(k-1),Vibrap(k-1)和Vibran(k-1)中的最小值,将该最小值赋值给Vibra0(k),同时将其对应的电流矢量的方向,即Pos0(k-1)、Posp(k-1)或Posn(k-1)赋值给Pos0(k);同时得到 Pos p ( k ) = Pos 0 ( k - 1 ) + π 2 k + 1 , Pos n ( k ) = Pos 0 ( k - 1 ) - π 2 k + 1 ; 施加Posp(k)和Posn(k)方向的电流矢量,每个电流矢量的作用时间为2n0·T0,并计算相应的位置波动量Vibrap(k)和Vibran(k);
步骤五、判断位置精度是否满足要求,如果是,转步骤六,如果否,转步骤三;
步骤六、找出步骤四中Vibra0(k),Vibrap(k)和Vibran(k)中的最小值,将该最小值对应的电流矢量的方向,即Pos0(k)、Posp(k)或Posn(k)赋值给PosCorrect0,并令iq_Ref=ΔA1,ΔA1≤ΔA0,转至步骤八;
步骤七、将电流幅值iq_Amp增加ΔA0,即iq_Amp=iq_Amp+ΔA0,其中ΔA0为每次增加的电流幅值,转至步骤二;
步骤八、施加幅值为iq_Amp,修正角度值脉冲宽度为T1的脉冲电流矢量,T1≤T0
步骤九、判断电机是否顺时针旋转,如果是,转向步骤十;如果否,转向步骤十一;
步骤十、获取转子初始位置修正值流程结束;
步骤十一、判断电机是否逆时针旋转,如果是,转向步骤十二;如果否,转向步骤十三;
步骤十二、获取转子初始位置修正值流程结束;
步骤十三、更新电流值iq_Amp=iq_Amp+ΔA1,转向步骤八。
本发明与现有技术相比的优点在于:
(1)本发明中施加正负脉冲的电流矢量,使转子高频抖动,可以克服静摩擦力的影响,提高位置检测精度。
(2)本发明通过正负脉冲的作用,避免了转子向一个方向运动范围过大,或时间过长,从而不会对下一个电流矢量作用时转子的运动造成影响,避免了判断失误的问题。
(3)本发明电流矢量的幅值和正负脉冲的作用时间可调,在保证转子振动的频率和幅值可调,增强了检测过程的灵活性和抗扰性。
附图说明
图1为本发明的流程图;
图2为驱动系统原理图;
图3是修正过程中电机转子的位置波动曲线;
图4是修正过程中转子电角度修正值的变化曲线。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
首先结合图2介绍一下永磁同步电机驱动控制系统的结构:主要包括永磁同步电机PMSM1、位置传感器2、电流传感器3、智能功率模块IPM4、数字控制器5。数字控制器5主要包括d轴电流控制器501、q轴电流控制器502、dq2abc变换模块503、abc2dq变换模块504、角度转换单元505。
其中永磁同步电机PMSM1的额定力矩为460Nm,额定电流为20A,极对数为12;位置传感器2可以是绝对式位置传感器或相对式位置传感器;电流传感器3采用高精度Hall电流传感器;智能功率模块IPM4为电压型逆变器,完成电压放大功能。数字控制器5主要完成基于磁场定向原理的电流闭环功能,并在此基础上完成图1所示的转子初始位置修正流程。
dq2abc变换模块503完成以下功能:将两维旋转坐标系转换到三维静止坐标系。具体如公式(1)所示:
C dq 2 abc = 2 3 cos θ e - sin θ e cos ( θ e - 2 π / 3 ) - sin ( θ e - 2 π / 3 ) cos ( θ e + 2 π / 3 ) - sin ( θ e + 2 π / 3 ) - - - ( 1 )
其中,θe表示电角度。
abc2dq变换模块504完成以下功能:将三维静止坐标系中的电流ia,ib,ic转换到两维旋转坐标系中的电流id,iq。具体如公式(2)所示:
C abc 2 dq = 2 3 cos θ e cos ( θ e - 2 π / 3 ) cos ( θ e - 4 π / 3 ) sin θ e sin ( θ e - 2 π / 3 ) sin ( θ e - 4 π / 3 ) - - - ( 2 )
其中,θe表示电角度。
d轴电流控制器501、q轴电流控制器502均采用PI算法,设计电流环带宽约500Hz。
角度转换单元505完成机械角度到电角度的转换。具体如公式(3)所示:
θe=Pair·θm+PosCorrect(3)
其中Pair为电机极对数,θm表示机械角度,θe表示电角度,PosCorrect表示电机转子初始位置的修正值(电角度)。
每次施加电流矢量的方向改变时,均将新的相位信息赋值给PosCorrect,即PosCorrect=Posx(k),其中x=0,p或n;Posx(k)表示为修正转子初始位置所加的电流矢量的方向;k表示方向改变的下标。
下面结合图1详细介绍转子初始位置的修正流程:
步骤一、初始化以下变量值:永磁同步电机d轴电流id_Ref=0,q轴电流的幅值iq_Amp=A0,电流矢量的方向Pos0(k)=0, (其中k=0);
步骤二、依次施加Pos0(k)、Posp(k)和Posn(k)方向的电流矢量,每个电流矢量的作用时间为2n0·T0(正负脉冲交替施加,各n0次;每个脉冲的宽度为T0;一般n0取值范围为2~5,此处取值为3,T0取值根据电机的响应时间选取,由于电机负载较大,此处选为100ms)。每次改变电流矢量的方向时,均将该方向赋值给PosCorrect,即PosCorrect=Posx(k),其中x=0,p或n。
计算每个电流矢量对应的位置波动:Vibra0(k),Vibrap(k),Vibran(k)。位置波动量采用以下算法:
Vibra x ( k ) = Σ i = 1 N | θ m ( i ) - 1 N Σ i = 1 N θ m ( i ) | - - - ( 4 )
其中,x=0,p或n;N表示在2n0·T0时间内位置采样点数;θm(i)表示电机转子的机械角度。
步骤三、判断电流矢量的幅值iq_Amp是否小于限定值iq_Max(一般可以设置为电机额定电流值的5%-50%),并且位置波动Vibra0(k),Vibrap(k)和Vibran(k)是否均小于某个限定值VibraMin(VibraMin根据位置传感器的精度和位置采样点数N设定),即判断(iq_Amp<iq_Max)&(Vibra0(k)<VibraMin)&(Vibrap(k)<VibraMin)&(Vibran(k)<VibraMin)是否为真。如果是,转向步骤七;如果否,转向步骤四。
步骤四、更新k=k+1,找出Vibra0(k-1),Vibrap(k-1)和Vibran(k-1)中的最小值,将该最小值赋值给Vibra0(k),将其对应的电流矢量的方向(即Pos0(k-1)、Posp(k-1)或Posn(k-1))赋值给Pos0(k);同时得到 施加Posp(k)和Posn(k)方向的电流矢量,每个电流矢量的作用时间为2n0·T0,并计算相应的位置波动量Vibrap(k)和Vibran(k)。
步骤五、判断位置精度是否满足要求(此处精度设为如果是,转步骤六,如果否,转步骤三。
步骤六、找出步骤四中Vibra0(k),Vibrap(k)和Vibran(k)中的最小值,将其对应的电流矢量的方向(即Pos0(k)、Posp(k)或Posn(k))赋值给PosCorrect0,并令iq_Ref=ΔA1(一般ΔA1≤ΔA0),转至步骤八。
步骤七、更新电流幅值iq_Amp=iq_Amp+ΔA0,其中ΔA0为每次增加的电流幅值,转至步骤二。
步骤八、施加幅值为iq_Amp,修正角度值脉冲宽度为T1的脉冲电流矢量(一般T1≤T0)。
步骤九、判断电机是否顺时针旋转。如果是,转向步骤十;如果否,转向步骤十一。
步骤十、获取转子初始位置修正值
步骤十一、判断电机是否逆时针旋转。如果是,转向步骤十二;如果否,转向步骤十三。
步骤十二、获取转子初始位置修正值结束。
步骤十三、更新电流值iq_Amp=iq_Amp+ΔA1(一般ΔA1≤ΔA0),转向步骤八。
经过上述过程,得到修正过程的主要曲线如图3和图4所示。图3是修正过程中电机转子的位置波动情况,可以看出在修正过程中,转子最大波动范围为235.5826°-235.2228°=0.36°(机械角度)。而图4是修正过程中转子电角度修正值PosCorrect的变化情况,在步骤六完成后,PosCorrect收敛到-27.78°(电角度),此时电机转子在幅值为5A,脉冲宽度为100ms的电流下始终保持静止,可以认为此时定子磁场和转子磁场方向平行,精度为(电角度),如果希望更高的精度,可以通过提高电流幅值和步骤五中收敛精度的判断值来达到。步骤八到步骤十三中,ΔA1=0.2A,T1=20ms,经过0.1s后,电机转子逆时针转动,所以执行步骤十二。至此,可以得到公式(3)中最终的修正值PosCorrect=-27.78°-90°=-117.78°(电角度),流程结束。
总之,本发明通过施加正负脉冲的电流矢量,使转子高频抖动,通过抖动量的大小,可以判断转子初始位置的偏差大小,通过改变电流矢量的方向,逐步减小偏差值。通过正负脉冲的作用,避免了转子向一个方向运动范围过大,或时间过长,从而不会对下一个电流矢量作用时转子的运动造成影响,避免了判断失误的问题。电流矢量的幅值和正负脉冲的作用时间可调,保证了转子振动的频率和幅值可调,增强了检测过程的灵活性和抗扰性。
本发明未详细阐述部分属于本领域公知技术。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭示的技术范围内,可想到的变换或替换,都应涵盖在本发明的包含范围之内。

Claims (1)

1.一种永磁同步电机转子初始位置的修正方法,其特征在于该方法包括以下步骤:
步骤一、初始化以下变量值:永磁同步电机d轴电流id_Ref=0,q轴电流的幅值iq_Amp=A0,电流矢量的方向Pos0(k)=0, 其中k=0;
步骤二、依次施加Pos0(k)、Posp(k)和Posn(k)方向的电流矢量,每个电流矢量包括n0个正脉冲和n0个负脉冲,正负脉冲交替施加,每个脉冲的宽度为T0,计算每个电流矢量对应的位置波动:Vibra0(k),Vibrap(k),Vibran(k);
步骤三、判断电流矢量的幅值iq_Amp是否小于限定值iq_Max,并且位置波动Vibra0(k),Vibrap(k)和Vibran(k)是否均小于某个限定值VibraMin,即判断(iq_Amp<iq_Max)&(Vibra0(k)<VibraMin)&(Vibrap(k)<VibraMin)&(Vibran(k)<VibraMin)是否为真;如果是,转向步骤七;如果否,转向步骤四;
步骤四、更新k=k+1,即将k加1,找出Vibra0(k-1),Vibrap(k-1)和Vibran(k-1)中的最小值,将该最小值赋值给Vibra0(k),同时将其对应的电流矢量的方向,即Pos0(k-1)、Posp(k-1)或Posn(k-1)赋值给Pos0(k);同时得到 施加Posp(k)和Posn(k)方向的电流矢量,每个电流矢量的作用时间为2n0·T0,并计算相应的位置波动量Vibrap(k)和Vibran(k);
步骤五、判断位置精度是否满足要求,如果是,转步骤六,如果否,转步骤三;
步骤六、找出步骤四中Vibra0(k),Vibrap(k)和Vibran(k)中的最小值,将该最小值对应的电流矢量的方向,即Pos0(k)、Posp(k)或Posn(k)赋值给PosCorrect0,并令iq_Amp=△A1,△A1≤△A0,转至步骤八;
步骤七、将电流幅值iq_Amp增加△A0,即iq_Amp=iq_Amp+△A0,其中△A0为每次增加的电流幅值,转至步骤二;
步骤八、施加幅值为iq_Amp,修正角度值脉冲宽度为T1的脉冲电流矢量,T1≤T0
步骤九、判断电机是否顺时针旋转,如果是,转向步骤十;如果否,转向步骤十一;
步骤十、获取转子初始位置修正值流程结束;
步骤十一、判断电机是否逆时针旋转,如果是,转向步骤十二;如果否,转向步骤十三;
步骤十二、获取转子初始位置修正值流程结束;步骤十三、更新电流值iq_Amp=iq_Amp+△A1,转向步骤八。
CN201310339651.9A 2013-08-06 2013-08-06 一种永磁同步电机转子初始位置的修正方法 Active CN103401504B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310339651.9A CN103401504B (zh) 2013-08-06 2013-08-06 一种永磁同步电机转子初始位置的修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310339651.9A CN103401504B (zh) 2013-08-06 2013-08-06 一种永磁同步电机转子初始位置的修正方法

Publications (2)

Publication Number Publication Date
CN103401504A CN103401504A (zh) 2013-11-20
CN103401504B true CN103401504B (zh) 2016-01-20

Family

ID=49565072

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310339651.9A Active CN103401504B (zh) 2013-08-06 2013-08-06 一种永磁同步电机转子初始位置的修正方法

Country Status (1)

Country Link
CN (1) CN103401504B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103916063B (zh) * 2014-03-31 2017-06-06 广东威灵电机制造有限公司 电机控制系统以及电机的转子的位置识别方法和装置
CN106125782B (zh) * 2016-08-31 2018-12-04 合肥金星机电科技发展有限公司 电机定位修正方法
CN106788072B (zh) * 2017-03-27 2019-04-16 安徽江淮汽车集团股份有限公司 永磁同步电机转子初始角度修正方法及修正系统
CN106772052B (zh) * 2017-03-27 2019-06-04 安徽江淮汽车集团股份有限公司 永磁同步电机转子初始角度修正方法及修正系统
CN107294459B (zh) * 2017-08-18 2019-08-06 安徽江淮汽车集团股份有限公司 永磁同步电机转子初始角度修正方法及修正系统
CN108700897A (zh) * 2017-10-09 2018-10-23 深圳市大疆灵眸科技有限公司 机械角度检测方法、云台及机器可读存储介质
CN108712128B (zh) * 2018-06-07 2021-10-01 南京信息职业技术学院 一种克服了摩擦力影响的交流伺服系统的对相方法
CN111224584A (zh) * 2020-03-02 2020-06-02 珠海格力电器股份有限公司 磁极位置的高效确定方法和装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369796A (zh) * 2007-08-17 2009-02-18 深圳市汇川技术股份有限公司 检测永磁同步电机转子磁极初始位置的方法和系统
CN101630938A (zh) * 2009-07-28 2010-01-20 哈尔滨工业大学 无位置传感器永磁同步电机转子初始位置辨识方法
CN102401626A (zh) * 2011-11-20 2012-04-04 中国科学院光电技术研究所 一种永磁同步电机转子位置传感器安装偏差估计方法
CN102414979A (zh) * 2009-05-27 2012-04-11 三菱电机株式会社 同步电动机的磁极位置估计装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013021843A (ja) * 2011-07-13 2013-01-31 Fuji Electric Co Ltd 永久磁石同期電動機の初期磁極位置調整装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101369796A (zh) * 2007-08-17 2009-02-18 深圳市汇川技术股份有限公司 检测永磁同步电机转子磁极初始位置的方法和系统
CN102414979A (zh) * 2009-05-27 2012-04-11 三菱电机株式会社 同步电动机的磁极位置估计装置
CN101630938A (zh) * 2009-07-28 2010-01-20 哈尔滨工业大学 无位置传感器永磁同步电机转子初始位置辨识方法
CN102401626A (zh) * 2011-11-20 2012-04-04 中国科学院光电技术研究所 一种永磁同步电机转子位置传感器安装偏差估计方法

Also Published As

Publication number Publication date
CN103401504A (zh) 2013-11-20

Similar Documents

Publication Publication Date Title
CN103401504B (zh) 一种永磁同步电机转子初始位置的修正方法
Li et al. Self-correction of commutation point for high-speed sensorless BLDC motor with low inductance and nonideal back EMF
CN103414427B (zh) 无刷直流电机控制方法
CN101594114B (zh) 永磁同步电机转子初始位置角的确定方法
CN103401502B (zh) 永磁同步电机初始角度辨识系统及方法
CN107465373B (zh) 基于线性霍尔传感器的直线电机自动门矢量控制的方法
CN103516284B (zh) 一种永磁同步电机电流增量预测算法
CN104779879B (zh) 永磁同步电机的电角度旋向和初始值辨识方法及系统
CN102938628A (zh) 一种永磁同步电机转子初始位置定位的方法
KR102588927B1 (ko) 모터 제어방법
CN103532465A (zh) 基于增量式模型参考自适应的永磁同步电机电感辨识算法
CN110212819B (zh) 一种用于高速无刷直流电机的换相误差补偿方法
CN104135197A (zh) 一种无传感器永磁同步电机调速控制策略
CN111371362B (zh) 高频注入法永磁直线电机动子位置估计的补偿方法
CN106130426A (zh) 基于ekf的无传感器超高速永磁同步电机转速控制方法
CN105915139A (zh) 一种矢量控制永磁同步伺服电机初始磁极位置搜索方法
CN104065319B (zh) 永磁同步电机零位初始角的标定方法
Jin et al. Closed-Loop compensation strategy of commutation error for sensorless brushless DC motors With nonideal asymmetric back-EMFs
CN105811828A (zh) 一种基于线性霍尔传感器的飞轮转速控制装置及方法
CN104167966A (zh) 一种带有霍尔传感器定位的永磁直线电机控制方法
CN109039198A (zh) 永磁电机霍尔位置的校正方法及装置
CN101282102A (zh) 兆瓦级永磁同步直驱风电机组无速度传感器控制方法
CN101814883B (zh) 基于位移传感器的永磁同步电机转子磁极位置检测方法
CN114157193B (zh) 一种寻优插值式同步电机转矩脉动抑制控制方法及系统
Huang et al. An iterative estimation algorithm of prepositioning focusing on the detent force in the permanent magnet linear synchronous motor system

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