CN106487304B - 一种基于滑模反电动势观测器的永磁同步电机状态估计方法 - Google Patents
一种基于滑模反电动势观测器的永磁同步电机状态估计方法 Download PDFInfo
- Publication number
- CN106487304B CN106487304B CN201610953899.8A CN201610953899A CN106487304B CN 106487304 B CN106487304 B CN 106487304B CN 201610953899 A CN201610953899 A CN 201610953899A CN 106487304 B CN106487304 B CN 106487304B
- Authority
- CN
- China
- Prior art keywords
- electromotive force
- observer
- sliding mode
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 21
- 238000009415 formwork Methods 0.000 title abstract 4
- 230000003068 static effect Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 238000004804 winding Methods 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 1
- 238000001914 filtration Methods 0.000 abstract description 2
- 244000145845 chattering Species 0.000 abstract 1
- 238000001514 detection method Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P21/0007—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode control
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开了一种基于滑模反电动势观测器的永磁同步电机状态估计方法,根据实时检测到的定子电压和定子电流,采用sigmoid函数取代符号函数,构造了一种滑模反电动势观测器对反电动势进行估计。该类型滑模反电动势观测器的特点在于对开关噪声不敏感并且具有滤波特性,因而使得估计的反电动势无需通过低通滤波器进行滤波。对估计的反电动势进行归一化处理后通过一种相当于简化的扩展卡尔曼滤波器的PI‑I结构的锁相环类型位置观测器提取转子位置和速度信息。本发明的永磁同步电机状态估计方法即使在存在噪声和电机参数变化的情况下仍具有较好的估计性能,并且减小了传统滑模观测器中存在的抖振现象,提高了对转子位置和速度的动态追踪性能和估计精度。
Description
技术领域
本发明涉及永磁同步电机无位置传感器控制领域,具体涉及一种永磁同步电机状态估计方法。
背景技术
近几十年来,永磁同步电动机以其体积小、效率高、功率密度大、动态性能好等优点在高性能伺服系统中获得了广泛的应用。对于永磁同步电机的矢量控制而言,获得准确的位置和速度信息十分关键。转子的位置和速度可以通过在转子转轴上安装机械传感器例如光电编码器来检测,但这些传感器不仅增加了成本而且降低了系统的可靠性,同时限制了永磁同步电机的应用范围。因此,研究永磁同步电机无传感器技术具有重要的理论和实用价值。
根据永磁同步电机在不同转速下转子位置估算效果,把无位置传感器控制方法分为两大类:
(1)基于电机的凸极(饱和)效应的转子位置估计方法。这类方法主要通过在电机注入高频信号条件下跟踪转子的空间凸极效应以实现对转子位置的估计。由于该类方法利用的是电机自身的凸极性,与电机转速及反电动势大小无关,因此能够实现在低速及零速条件下对转子位置的估计。但这类方法的有效性易受到实际电机的非线性磁性行为的影响,同时在提取转子位置信息的解调过程中低通滤波器或带通滤波器的使用会带来相移及幅值畸变等问题。
(2)基于基波激励下电机数学模型的转子位置估计方法。这类方法不需要注入额外的信号,适用于在中高速运行时对转子位置的估计,具体方法有:利用定子端电压和电流的直接计算法,基于观测器的闭环算法,包括扩展卡尔曼滤波法,模型参考自适应法,滑模观测器法等。
基于基波激励下电机数学模型的转子位置估计方法易受到电机参数变化、负载扰动以及未建模动态的影响。与其他方法相比,在中高速条件下,滑模观测器法具有响应快速,对参数变化及扰动不灵敏,计算量较小,观测结果具有很好的鲁棒性的优点。由于滑模变结构控制本质上是不连续的开关控制,因此对一个实际的滑模观测器而言,由于开关时间和空间上的滞后,使得滑模观测器不可避免的存在抖振现象。在传统滑模观测器中估计的反电动势抖振现象较为严重,解决方法通常是采用低通滤波器进行滤波,但这会造成估计的转子位置存在相位滞后。
发明内容
为解决传统滑模观测器所存在的问题,本发明提出一种能够有效减小传统滑模观测器中存在的抖振,提高对永磁同步电机转子位置和速度的动态追踪性能和估计精度的方法。
为实现上述目的,本发明所采用的技术方案是一种基于滑模反电动势观测器的状态估计方法,该方法包括以下步骤:
步骤(1):检测电机的三相绕组电流ia、ib、ic,通过CLARK变换得到定子电流在两相静止坐标系下a、β轴的分量ia、iβ,根据定子电流在a、β轴的分量ia、iβ,定子电流估计值在a、β轴的分量以及定子电压在a、β轴的分量ua、uβ构造一种滑模反电动势观测器,利用所构造的滑模反电动势观测器对反电动势进行估计,估计的反电动势为
步骤(2):对滑模反电动势观测器估计得到的反电动势经归一化处理后输入到采用PI-I结构的锁相环类型位置观测器中提取转子位置和转速信息。
所述步骤(1)构造的滑模观测器为:
其中,ia、iβ分别为定子电流在两相静止坐标系下a、β轴的分量,ua、uβ分别为定子电压在两相静止坐标系下a、β轴的分量;
分别为α、β轴定子电流的估计值,分别为a、β轴反电动势的估计值;
H为sigmoid函数,k1~4是观测器增益;
RS、LS分别为定子电阻和电感;
sigmoid函数表达式为:
设
将滑模超平面定义为:
S=[Sa Sβ]T=[e1 e2]T
观测器增益取值范围为:
本发明相比现有技术具有以下优点及有益效果:
(1)本发明所提出永磁同步电机状态估计方法采用sigmoid函数代替符号函数,有效的抑制了传统滑模观测器所存在的抖振现象。
(2)通过滑模反电动势观测器直接对反电动势进行估计,由于没有传统滑模观测器估计所需的低通滤波器环节,从而简化了系统结构和消除了低通滤波器所造成的时延。
(3)采用一种比例积分(PI)加积分(I)结构的锁相环类型位置观测器从估计的反电动势中提取转子位置和速度信息,提高了对转子位置和速度的动态追踪性能和估计精度。
附图说明
图1是本发明的滑模反电动势观测器原理框图;
图2是本发明的PI-I结构的锁相环类型位置观测器原理框图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
图1和图2分别是本发明的滑模反电动势观测器原理框图和PI-I结构的锁相环类型位置观测器原理框图,具体实施步骤为:
步骤(1):构造滑模反电动势观测器,通过该观测器估计的反电动势为
在两相静止坐标系下定子电流状态方程可表示为:
其中,永磁同步电机的反电动势为:
其中,ia、iβ分别为定子电流在两相静止坐标系下a、β轴的分量,ua、uβ分别为定子电压在两相静止坐标系下a、β轴的分量,RS、LS分别为定子电阻和电感,ωe表示转子电角速度,ψf表示永磁磁链,θe表示转子位置。
假设反电动势的变化与其他变量相比足够慢,由两相静止坐标系下定子电流状态方程可得到如下的动态方程:
所构造的滑模反电动势观测器为:
其中,分别为a、β轴定子电流的估计值,k1~4是观测器增益;
H为sigmoid函数,sigmoid函数表达式为:
a表示倾斜系数。
滑模反电动势观测器状态估计动态误差方程为:
其中,
将滑模超平面定义为:
S=[Sa Sβ]T=[e1 e2]T
取李雅普诺夫函数为:
对李雅普诺夫函数求导:
为了满足使观测器增益k1、k2的取值范围为:
k1>|e3|max/Ls,k2>|e4|max/Ls
当系统处于滑动模态时:
在这种情况下,有:
此时,简化的反电动势动态方程为:
当时,跟踪误差e3、e4指数收敛至零。
因此,观测器增益k1~4的取值范围为:
步骤(2):对估计反电动势经归一化处理后输入到采用PI-I结构的锁相环类型位置观测器中提取出转子位置和速度具体如下:
步骤(21):对估计的反电动势进行归一化处理:
其中,是归一化的滑模反电动势观测器输出,作为锁相环类型位置观测器的输入,na、nβ分别代表中存在的谐波和噪声。
利用外差法计算转子位置估计误差:
式中,ε包括一个包含有转子位置估计误差的主要分量ε'以及两个经高频信号和调制的噪声相关项,表示转子位置估计值。
当转子位置估计误差时,上式可化简为:
步骤(22):采用PI-I结构的锁相环类型的位置观测器从ε'中提取出转子位置及速度信息并消除噪声相关项。
转子速度估计值表示为:
转子位置估计值表示为:
式中,K0表示积分系数,K1、K2表示比例系数;该表达式增加了一个与ε成比例的前馈信号以增强观测器对转子位置估计的追踪性能。
整个位置观测器的传递函数可以表示为:
上式的特征多项式为:
S3+K2S2+K1S+K0=(S+ωn)(S2+2ζnωnS+ωn 2)
式中,K0=ωn 3,K1=(2ζn+1)ωn 2,K2=(2ζn+1)ωn
可以看出观测器的性能由ωn和ζn决定,ωn表示固有频率,ζn表示阻尼系数。
注意到本发明位置观测器的传递函数与扩展卡尔曼滤波法类似,将传递函数写成如下形式:
其中,Ts为采样周期,w'为ε(k)的积分,ke1,ke2和ke3为卡尔曼滤波器增益,与位置观测器增益的关系如下:
由此可见,本发明的位置估计精度能够达到类似于扩展卡尔曼滤波器的位置估计效果。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,包括如下步骤:
步骤(1),构造滑模反电动势观测器,通过该观测器估计反电动势为包括:
检测电机的三相绕组电流ia、ib、ic,通过CLARK变换得到定子电流在两相静止坐标系下α、β轴的分量iα、iβ,根据定子电流在α、β轴的分量iα、iβ,定子电流估计值在α、β轴的分量以及定子电压在α、β轴的分量uα、uβ,构造一种滑模反电动势观测器,利用所构造的滑模反电动势观测器对反电动势进行估计,得到反电动势为
步骤(2),对滑模反电动势观测器估计得到的反电动势经归一化处理后输入到采用PI-I结构的锁相环类型位置观测器中提取出转子位置估计值和速度
2.根据权利要求1所述的一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,所述步骤(1)中构造的滑模反电动势观测器的设计方法包括:
建立滑模反电动势观测器的表达式:
其中,iα、iβ分别为定子电流在两相静止坐标系下α、β轴的分量,uα、uβ分别为定子电压在两相静止坐标系下α、β轴的分量;
分别为α、β轴定子电流的估计值,分别为α、β轴反电动势的估计值;
H(x)为sigmoid函数,k1~4是观测器增益;
RS、LS分别为定子电阻和电感;
建立滑模反电动势观测器状态估计动态误差方程:
其中,其中,eα、eβ为永磁同步电机反电动势在α、β轴的实际值分量;
定义滑模超平面:
S=[Sα Sβ]T=[e1 e2]T
定义李雅普诺夫函数为:
对李雅普诺夫函数求导:
为了满足使观测器增益k1、k2的取值范围为:
k1>|e3|max/Ls,k2>|e4|max/Ls,
3.根据权利要求2所述的一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,所述sigmoid函数表达式为:
其中a表示倾斜系数。
4.根据权利要求1所述的一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,所述步骤(2)具体步骤包括如下:
步骤(21):对估计的反电动势采用下式进行归一化处理:
其中,是归一化的滑模反电动势观测器输出,作为锁相环类型位置观测器的输入,nα、nβ分别代表中存在的谐波和噪声,θe表示转子位置实际值;
步骤(22):采用PI-I结构的锁相环类型位置观测器提取转子位置及速度信息并消除噪声相关项。
5.根据权利要求4所述的一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,所述PI-I结构的锁相环类型位置观测器的传递函数的表达式为:
S为滑模超平面。
6.根据权利要求5所述的一种基于滑模反电动势观测器的永磁同步电机状态估计方法,其特征在于,所述步骤(22)的实现包括:利用外差法计算转子位置估计误差ε,将ε输入所述位置观测器中,得到转子速度估计值为:
S为滑模超平面
转子位置估计值表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610953899.8A CN106487304B (zh) | 2016-10-27 | 2016-10-27 | 一种基于滑模反电动势观测器的永磁同步电机状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610953899.8A CN106487304B (zh) | 2016-10-27 | 2016-10-27 | 一种基于滑模反电动势观测器的永磁同步电机状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106487304A CN106487304A (zh) | 2017-03-08 |
CN106487304B true CN106487304B (zh) | 2018-10-09 |
Family
ID=58271561
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610953899.8A Expired - Fee Related CN106487304B (zh) | 2016-10-27 | 2016-10-27 | 一种基于滑模反电动势观测器的永磁同步电机状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106487304B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107154761B (zh) * | 2017-05-26 | 2019-12-06 | 北京工业大学 | 发电机全控整流系统和控制方法 |
CN107809191B (zh) * | 2017-10-24 | 2020-01-07 | 厦门金龙汽车新能源科技有限公司 | 一种同步电机无速度传感器的角度观测方法 |
CN107911057B (zh) * | 2017-11-13 | 2020-02-21 | 南京工程学院 | 一种用于飞轮储能系统的转子位置鲁棒观测方法 |
CN108390605A (zh) * | 2018-03-27 | 2018-08-10 | 江苏大学 | 基于模糊变结构模型参考自适应观测器的永磁同步电机状态估计方法 |
CN108681255A (zh) * | 2018-05-16 | 2018-10-19 | 江苏大学 | 一种基于滑模变结构控制的削弱磁悬浮飞轮抖振的方法 |
CN109194224B (zh) * | 2018-08-29 | 2022-02-18 | 南京理工大学 | 基于扩张状态观测器的永磁同步电机无传感器控制方法 |
CN108964546A (zh) * | 2018-09-07 | 2018-12-07 | 上海理工大学 | 基于自适应滑模观测器检测电机转子位置及转速的方法 |
CN109660168B (zh) * | 2018-12-29 | 2020-11-17 | 珠海格力电器股份有限公司 | 一种电机控制方法、系统及电机 |
CN110426957B (zh) * | 2019-07-31 | 2020-03-13 | 深圳信息职业技术学院 | 一种基于时延观测器的水厂投药系统自适应滑模控制方法 |
CN110677087A (zh) * | 2019-11-20 | 2020-01-10 | 四川虹美智能科技有限公司 | 变频电机控制方法 |
CN110995094B (zh) * | 2019-12-09 | 2021-04-13 | 哈尔滨工业大学 | 一种低载波比控制下的永磁同步电机位置和速度估算方法 |
CN110971167B (zh) * | 2019-12-24 | 2021-10-12 | 江苏大学 | 基于检波滤波器的变漏磁电机无位置传感器控制方法 |
CN111478637B (zh) * | 2020-04-16 | 2022-06-14 | 广东美的智能科技有限公司 | 电机控制方法及电机控制系统 |
CN111669093B (zh) * | 2020-06-23 | 2022-07-08 | 南京邮电大学 | 一种基于自适应扩展卡尔曼滤波的电机参数估计方法 |
CN111969921B (zh) * | 2020-07-20 | 2022-04-19 | 四川虹美智能科技有限公司 | 电机转速的确定方法及装置、电子设备 |
CN111987959A (zh) * | 2020-08-11 | 2020-11-24 | 哈尔滨理工大学 | 永磁同步电机无传感器控制方法 |
CN112671298B (zh) * | 2021-01-13 | 2023-03-21 | 湖南进芯电子科技有限公司 | 改进的用于永磁同步电机控制的pll无感控制算法 |
CN113098350B (zh) * | 2021-05-24 | 2022-05-31 | 无锡威孚高科技集团股份有限公司 | 基于滑模观测器的表贴式永磁同步电机转子位置计算方法和装置 |
CN113904605B (zh) * | 2021-09-29 | 2024-03-19 | 上海氢恒汽车电子有限公司 | 永磁同步电机无位置传感器转子位置的确定方法及装置 |
CN114598206B (zh) * | 2022-03-08 | 2024-04-09 | 天津农学院 | 永磁同步电机宽速域转子位置观测器设计方法 |
CN115459656A (zh) * | 2022-09-19 | 2022-12-09 | 江苏科技大学 | 用于agv搬运机器人的永磁同步电机无传感器控制方法 |
CN117040342A (zh) * | 2023-08-10 | 2023-11-10 | 西安电子科技大学 | 一种基于复数型归一化的电机位置和速度估计方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103414423A (zh) * | 2013-08-22 | 2013-11-27 | 东南大学 | 一种面贴式永磁同步电机无位置传感器直接转矩控制方法 |
CN105227010A (zh) * | 2015-10-23 | 2016-01-06 | 哈尔滨工业大学 | 一种永磁同步电机无位置传感器位置观测误差谐波脉冲消除方法 |
CN105490608A (zh) * | 2015-12-30 | 2016-04-13 | 上海贝思特门机有限公司 | 一种永磁电梯门机控制器及其控制方法 |
CN106026801A (zh) * | 2016-07-04 | 2016-10-12 | 珠海格力节能环保制冷技术研究中心有限公司 | 永磁同步电机转子位置的检测方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4480696B2 (ja) * | 2005-08-26 | 2010-06-16 | 三洋電機株式会社 | モータ制御装置 |
-
2016
- 2016-10-27 CN CN201610953899.8A patent/CN106487304B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103414423A (zh) * | 2013-08-22 | 2013-11-27 | 东南大学 | 一种面贴式永磁同步电机无位置传感器直接转矩控制方法 |
CN105227010A (zh) * | 2015-10-23 | 2016-01-06 | 哈尔滨工业大学 | 一种永磁同步电机无位置传感器位置观测误差谐波脉冲消除方法 |
CN105490608A (zh) * | 2015-12-30 | 2016-04-13 | 上海贝思特门机有限公司 | 一种永磁电梯门机控制器及其控制方法 |
CN106026801A (zh) * | 2016-07-04 | 2016-10-12 | 珠海格力节能环保制冷技术研究中心有限公司 | 永磁同步电机转子位置的检测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106487304A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106487304B (zh) | 一种基于滑模反电动势观测器的永磁同步电机状态估计方法 | |
An et al. | Frequency-adaptive complex-coefficient filter-based enhanced sliding mode observer for sensorless control of permanent magnet synchronous motor drives | |
CN107482977B (zh) | 一种永磁同步电机转子位置和转速检测方法 | |
CN103199779B (zh) | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 | |
CN110198150B (zh) | 一种永磁同步电机多参数在线辨识方法 | |
Feng et al. | Hybrid terminal sliding-mode observer design method for a permanent-magnet synchronous motor control system | |
US7098623B2 (en) | Method in connection with permanent magnet synchronous machines | |
Wang et al. | Review of sensorless control techniques for PMSM drives | |
CN105227010B (zh) | 一种永磁同步电机无位置传感器位置观测误差谐波脉冲消除方法 | |
CN103199788B (zh) | 永磁同步电机转子位置观测器 | |
CN107579690B (zh) | 一种基于滑模观测的超高速永磁同步电机转速估计方法 | |
CN103414423A (zh) | 一种面贴式永磁同步电机无位置传感器直接转矩控制方法 | |
CN108900127B (zh) | 考虑交叉耦合效应的ipmsm低速段无位置传感器控制方法 | |
CN114598206B (zh) | 永磁同步电机宽速域转子位置观测器设计方法 | |
CN105515488B (zh) | 一种基于自适应滤波器同步电机低速段无传感器控制方法 | |
CN115173774B (zh) | 一种永磁同步电机无位置传感器控制方法及系统 | |
CN106026831B (zh) | 一种无位置传感器控制下的电流解耦方法 | |
CN106411209A (zh) | 永磁同步电机无位置传感器控制方法 | |
Wang et al. | Rotor position estimation method for permanent magnet synchronous motor based on super-twisting sliding mode observer | |
CN108512475A (zh) | 一种基于脉振高频电压注入法的转子位置检测方法 | |
CN108322121A (zh) | 一种无人机电机转子初始位置检测方法、装置、设备及存储介质 | |
CN110661466B (zh) | 准比例谐振自适应观测器及永磁同步电机位置估算方法 | |
CN113037163A (zh) | 一种新型具有动态误差补偿的永磁同步电机位置观测方法 | |
CN110971167B (zh) | 基于检波滤波器的变漏磁电机无位置传感器控制方法 | |
CN108270373B (zh) | 一种永磁同步电机转子位置检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181009 Termination date: 20191027 |
|
CF01 | Termination of patent right due to non-payment of annual fee |