CN110277941B - 永磁电机、转子位置控制方法、装置 - Google Patents
永磁电机、转子位置控制方法、装置 Download PDFInfo
- Publication number
- CN110277941B CN110277941B CN201910571898.0A CN201910571898A CN110277941B CN 110277941 B CN110277941 B CN 110277941B CN 201910571898 A CN201910571898 A CN 201910571898A CN 110277941 B CN110277941 B CN 110277941B
- Authority
- CN
- China
- Prior art keywords
- rotor
- permanent magnet
- magnet motor
- time
- position control
- 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 26
- 230000003044 adaptive effect Effects 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 238000012423 maintenance Methods 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
Images
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/0017—Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation
-
- 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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- 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
- H02P2203/00—Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
- H02P2203/03—Determination of the rotor position, e.g. initial rotor position, during standstill or low speed operation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electric Motors In General (AREA)
- Control Of Ac Motors In General (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开了一种永磁电机、转子位置控制方法、装置,所述方法:建立永磁电机的位置控制模型;提取所述位置控制模型中的时变特征参数;采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。采用上述方案,控制模型中的时变特征参数参与自适应控制器,使得时变特征参数最终会趋于恒定值,实现精确性较高且稳定性、鲁棒性较强的控制效果。
Description
技术领域
本发明涉及自动化控制领域,尤其涉及一种永磁电机、转子位置控制方法、装置。
背景技术
在高精度控制领域,良好的精密驱动平台是保证高质量控制效果的必要条件。
现有技术中,永磁电机作为驱动平台,而永磁电机是通过其的转子驱动其他设备。在永磁电机中往往采用比例积分控制器进行转子的驱动控制,然而比例积分控制器对系统的依赖性高,当存在参数的变化和外界扰动,尤其是在被驱动设备收到多方面因素的干扰的情况下,其控制效果的精确性较低。
发明内容
发明目的:针对现有技术缺陷,本发明旨在提供一种永磁电机、转子位置控制方法、装置,以实现精确性较高的控制效果。
技术方案:本发明实施例中提供一种永磁电机的转子位置控制方法,包括:建立永磁电机的位置控制模型;提取所述位置控制模型中的时变特征参数;采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
具体的,控制所述永磁电机的d轴电流id=0。
具体的,以永磁电机的位置环控制器、转速环控制器和电流环控制器建立三阶离散差分形式的特征模型:
θ(k)=f1(k)θ(k-1)+f2(k)θ(k-2)+f3(k)θ(k-3)+g0(k)Uq(k-1),
其中,θ为转子位置,Uq为永磁电机q轴输入电压,k为时序,f1,f2,f3和g0为时变特征参数。
具体的,以比例控制器建立位置环控制器;以永磁电机的转速环控制器和电流环控制器建立二阶离散差分形式的特征模型:
ω(k)=f1(k)ω(k-1)+f2(k)ω(k-2)+g0(k)+Uq(k-1),
其中,ω为转子角速度。
具体的,采用如下公式实时辨识时变特征参数:
φ(k)=[ω(k-1) ω(k-2) Uq(k-1)],
具体的,采用维持跟踪控制器:
其中,ωr *(k)为给定转子角转速,ω(k)为实际转子角转速,kw为正常参数;
采用黄金分割自适应控制器:
其中,l1=0.382,l2=0.618,e(k)为转速误差,η1,η2,μ,kL为正常参数,其中Uq=u0+u。
具体的,修正Uq因干扰而产生的误差,其中:
采用逻辑积分控制器:
uI(k)=uI(k-1)+kIe(k)
其中,kI为积分系数,k1,k2为正常参数,k2>k1;
采用逻辑微分控制器:
其中,kd为微分系数,kd1,kd2为正常参数,ε为给定误差带,y为给定转子角转速与实际转子角速度之差,△t为采样时间,Uq=u0+u+uI+ud。
具体的,采用混合粒子群算法对正常参数η1,η2,μ,kL和kw进行离线优化。
本发明实施例中还提供一种永磁电机的转子位置控制装置,包括:模型建立单元,提取单元和自适应单元,其中:所述模型建立单元,用于建立永磁电机的位置控制模型;所述提取单元,用于提取所述位置控制模型中的时变特征参数;所述自适应单元,用于采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
本发明实施例中还提供一种永磁电机,采用根据权利要求1~8任意一项所述的永磁电机的转子位置控制方法。
有益效果:与现有技术相比,本发明具有如下显著优点:由于永磁电机系统没有使用精确建模,无需传递函数推导,因此减少了大量的数值计算过程。采用控制模型中的时变特征参数参与自适应控制器,使得时变特征参数最终会趋于恒定值,实现精确性较高且稳定性、鲁棒性较强的控制效果。
附图说明
图1为本发明实施例中永磁电机的转子位置控制方法的流程示意图;
图2为本发明实施例中混合粒子群算法优化正常参数的过程示意图;
图3为本发明实施例中永磁电机的转子位置控制装置的结构示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
参阅图1,其为本发明实施例中永磁电机的转子位置控制方法的流程示意图,其中包括具体步骤,以下结合具体步骤进行详细说明。
步骤S101,建立永磁电机的位置控制模型。
本发明实施例中,控制所述永磁电机的d轴电流id=0。
本发明实施例中,以永磁电机的位置环控制器、转速环控制器和电流环控制器建立三阶离散差分形式的特征模型:
θ(k)=f1(k)θ(k-1)+f2(k)θ(k-2)+f3(k)θ(k-3)+g0(k)Uq(k-1),
其中,θ为转子位置,Uq为永磁电机q轴输入电压,k为时序,f1,f2,f3和g0为时变特征参数。
本发明实施例中,以比例控制器建立位置环控制器;以永磁电机的转速环控制器和电流环控制器建立二阶离散差分形式的特征模型:
ω(k)=f1(k)ω(k-1)+f2(k)ω(k-2)+g0(k)+Uq(k-1),
其中,ω为转子角速度。
在具体实施中,永磁电机的控制结果可以由永磁电机的转子位置体现。
在具体实施中,在永磁电机控制中,为了能够得到类似直流电机的控制特性,在电机转子上可以建立了坐标系,此坐标系与转子同步转动,取转子磁场方向为d轴,垂直于转子磁场方向为q轴,将电机的数学模型转换到此坐标系下,可实现d轴和q轴的解耦。控制所述永磁电机的d轴电流id=0,可以避免来自d轴的扰动,进一步提升控制结果的精确性。
在具体实施中,通过机械运动方程、q轴电压方程和电磁转矩方程,可以推导出
其中,θ为转子位置,J为电机转动惯量,B为阻尼系数,Te为电磁转矩,Pn为极对数,ψf为永磁体磁链,Rs为定子电阻,L为定子电感,iq为定子电流的q轴分量,uq为定子电压q轴分量。则可以建立以转子位置θ为输出,以q轴等效电压Uq为输入的三阶离散差分形式的特征模型为:
θ(k)=f1(k)θ(k-1)+f2(k)θ(k-2)+f3(k)θ(k-3)+g0(k)Uq(k-1),
其中,θ为转子位置,Uq为永磁电机q轴输入电压,k为时序,f1,f2,f3和g0为时变特征参数。
在具体实施中,k是离散的计数下的时序状态,θ(k)表示在k时序状态下的转子位置,其他相同的表示形式同理。
在具体实施中,为增强永磁电机系统的抗干扰能力,将三阶离散差分形式的特征模型可以分解为以比例控制器建立的位置环控制器,和永磁电机的转速环控制器和电流环控制器建立二阶离散差分形式的特征模型:
ω(k)=f1(k)ω(k-1)+f2(k)ω(k-2)+g0(k)+Uq(k-1),
其中,ω为转子角速度。
在具体实施中,f1,f2,f3和g0为时变特征参数,即随时序状态k而变化的特征参数,可以体现为是外界对永磁电机的影响因素和电机运行中自身变化的影响因素的体现。
在具体实施中,通过对上述输入二阶离散差分形式的特征模型的Uq进行补偿、抗干扰、自适应调整,使输出的转子角速度达到比较稳定的状态,进而转子角速度决定转子位置不受外界和自身变化的干扰影响,最终实现控制驱动效果的稳定性和精确性。
步骤S102,提取所述位置控制模型中的时变特征参数。
本发明实施例中,采用如下公式实时辨识时变特征参数:
φ(k)=[ω(k-1)ω(k-2)Uq(k-1)],
在具体实施中,由于时变特征参数是随时序变化的,因此可以实时的获取时变特征参数。
步骤S103,采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
在具体实施中,可以通过自适应控制器对输入电压进行控制Uq,实现快速跟踪给定目标的目的,进而使得电机系统最终的驱动控制效果稳定、精确,给定目标则可以是理想状态下的转子位置、角速度。
本发明实施例中,自适应控制器包括:维持跟踪控制器和黄金分割自适应控制器,其中:
维持跟踪控制器:
其中,ωr *(k)为给定转子角转速,ω(k)为实际转子角转速,kw为正常参数;
黄金分割自适应控制器:
其中,l1=0.382,l2=0.618,e(k)为转速误差,η1,η2,μ,kL为正常参数,其中Uq=u0+u。
在具体实施中,e(k)为转速误差,即给定转子角转速与实际转子角速度之差,给定转子角速度可以是指上述给定目标,kw、kL为可以防止表达式分母为零的正常参数,上标弧线表示该参数处于变化中。
在具体实施中,转子输入电压Uq=u0+u,表示经控制器控制后理想状态下的输入电压。
本发明实施例中,修正Uq因干扰而产生的误差,其中:
采用逻辑积分控制器:
uI(k)=uI(k-1)+kIe(k)
其中,kI为积分系数,k1,k2为正常参数,k2>k1;
采用逻辑微分控制器:
其中,kd为微分系数,kd1,kd2为正常参数,ε为给定误差带,y为给定转子角转速与实际转子角速度之差,△t为采样时间,Uq=u0+u+uI+ud。
在具体实施中,y(k)和e(k)可以表示相同的概念,ε为给定误差带,具体可以由多次试验情况获得。
在具体实施中,uI+ud表示经控制器对干扰的修正补偿,因此Uq=u0+u+uI+ud表示抗干扰状态下理想的转子q轴输入电压,转子角速度等因k变化而变化的变量在变化时,也可以导致Uq的变化。
本发明实施例中,可以采用混合粒子群算法对正常参数η1,η2,μ,kL和kw进行离线优化。
参阅图2,其为本发明实施例中混合粒子群算法优化正常参数的过程示意图。
在具体实施中,HPSO算法(HPSO,Hybrid Particle Swarm Optimization)与Simulink模型之间进行连接,Simulink模型可以将混合粒子群算法中产生的粒子群赋值给η1,η2,μ,kL和kw,然而运行控制系统模型,输出性能指标,混合粒子群算法根据性能指标判断η1,η2,μ,kL和kw是否符合优化标准,若不符合,则进行更新粒子,个体最优交叉,群体最优交叉,粒子变异,产生粒子群的循环操作,直到η1,η2,μ,kL和kw符合优化标准时停止。
在具体实施中,混合粒子群算法与Simulink之间连接的桥梁是粒子(即自适应控制器参数)和该粒子对应的适应度值(即控制系统的性能指标)。
在具体实施中,本发明提供永磁电机的转子位置控制方法由于设计的自适应抗扰控制可以自适应调整控制系统参数,因此抗扰性与鲁棒性均较好,适用于大口径天文望远镜等类似高精度且受较多外界干扰的设备。其次,由于电机系统不需要精确建模,无需传递函数推导,因此可以减少大量的数值计算过程,控制器设计简便,不需要现场反复调试。由于参数辨识最终会趋于恒定值,因此电机控制系统稳定可靠。
在具体实施中,由于望远镜口径的增大,且大口径天文望远镜伺服系统除了来自电机自身因素的干扰,还存在风阻力矩、摩擦力矩等其它干扰力矩的影响,因此在方位轴和高度轴的轴系驱动上的精确度和稳定度要求越来越高,而本发明实施例中提供的方法应用于永磁电机可以实现望远镜的俯仰轴的精准控制。
参阅图3,其为本发明实施例中永磁电机的转子位置控制装置30的结构示意图,具体包括:模型建立单元301,提取单元302和自适应单元303,其中:所述模型建立单元301,用于建立永磁电机的位置控制模型;所述提取单元302,用于提取所述位置控制模型中的时变特征参数;所述自适应单元303,用于采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
本发明还提供一种永磁电机,采用本发明实施例中任意一种的永磁电机的转子位置控制方法。
Claims (9)
1.一种永磁电机的转子位置控制方法,其特征在于,包括:
建立永磁电机的位置控制模型,以永磁电机的位置环控制器、转速环控制器和电流环控制器建立三阶离散差分形式的特征模型:
θ(k)=f1(k)θ(k-1)+f2(k)θ(k-2)+f3(k)θ(k-3)+g0(k)Uq(k-1),
其中,θ为转子位置,Uq为永磁电机q轴输入电压,k为时序,f1,f2,f3和g0为时变特征参数;
提取所述位置控制模型中的时变特征参数;
采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
2.根据权利要求1所述的永磁电机的转子位置控制方法,其特征在于,在所述建立永磁电机的位置控制模型之前,包括:
控制所述永磁电机的d轴电流id=0。
3.根据权利要求2所述的永磁电机的转子位置控制方法,其特征在于,所述建立永磁电机的位置控制模型,包括:
以比例控制器建立位置环控制器;以永磁电机的转速环控制器和电流环控制器建立二阶离散差分形式的特征模型:
ω(k)=f1(k)ω(k-1)+f2(k)ω(k-2)+g0(k)+Uq(k-1),
其中,ω为转子角速度。
7.根据权利要求6所述的永磁电机的转子位置控制方法,其特征在于,所述采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,包括:
采用混合粒子群算法对正常参数η1,η2,μ,kL和kw进行离线优化。
8.一种永磁电机的转子位置控制装置,其特征在于,包括:模型建立单元,提取单元和自适应单元,其中:
所述模型建立单元,用于建立永磁电机的位置控制模型,以永磁电机的位置环控制器、转速环控制器和电流环控制器建立三阶离散差分形式的特征模型:
θ(k)=f1(k)θ(k-1)+f2(k)θ(k-2)+f3(k)θ(k-3)+g0(k)Uq(k-1),
其中,θ为转子位置,Uq为永磁电机q轴输入电压,k为时序,f1,f2,f3和g0为时变特征参数;
所述提取单元,用于提取所述位置控制模型中的时变特征参数;
所述自适应单元,用于采用基于所述时变特征参数的自适应控制器对转子输入电压进行控制,由控制后的转子输入电压决定所述转子位置。
9.一种永磁电机,其特征在于,采用根据权利要求1~7任意一项所述的永磁电机的转子位置控制方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910571898.0A CN110277941B (zh) | 2019-06-28 | 2019-06-28 | 永磁电机、转子位置控制方法、装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910571898.0A CN110277941B (zh) | 2019-06-28 | 2019-06-28 | 永磁电机、转子位置控制方法、装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110277941A CN110277941A (zh) | 2019-09-24 |
CN110277941B true CN110277941B (zh) | 2021-01-05 |
Family
ID=67963118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910571898.0A Expired - Fee Related CN110277941B (zh) | 2019-06-28 | 2019-06-28 | 永磁电机、转子位置控制方法、装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110277941B (zh) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3417688A1 (de) * | 1984-05-12 | 1985-11-14 | Mergenthaler Linotype Gmbh, 6236 Eschborn | Verfahren und einrichtung zur gleichlaufregelung eines elektromotors, insbesondere zum antrieb eines drehbaren reflektierenden elements eines optischen abtastsystems |
JP2007252138A (ja) * | 2006-03-17 | 2007-09-27 | Nissan Motor Co Ltd | 電動機制御装置 |
CN101977009B (zh) * | 2010-10-09 | 2013-01-02 | 江苏中容电气有限公司 | 高精度数控机床进给驱动用正弦波直线电机的控制方法 |
CN102158164B (zh) * | 2011-03-02 | 2014-10-01 | 中国电力科学研究院 | 一种磁控式并联电抗器的梯形变权回归控制方法 |
CN102751935B (zh) * | 2012-07-03 | 2014-11-05 | 中国东方电气集团有限公司 | 一种在电动车转速开环下抑制异步电机干扰的方法 |
CN103701396B (zh) * | 2013-12-13 | 2016-01-20 | 天津大学 | 一种基于自适应模糊神经网络的电机转速跟踪控制方法 |
CH712829A1 (de) * | 2016-08-22 | 2018-02-28 | Lakeview Innovation Ltd | Verfahren zur sensorlosen Steuerung eines PMSM-Motors. |
CN108880358A (zh) * | 2018-07-02 | 2018-11-23 | 哈尔滨理工大学 | 基于角位移无时滞观测器的永磁同步电机控制方法及装置 |
CN109379003A (zh) * | 2018-11-27 | 2019-02-22 | 北京轩宇空间科技有限公司 | 高速永磁同步电机控制的方法、装置及计算机存储介质 |
CN109560736B (zh) * | 2018-12-18 | 2020-03-31 | 东南大学 | 基于二阶终端滑模的永磁同步电机控制方法 |
-
2019
- 2019-06-28 CN CN201910571898.0A patent/CN110277941B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110277941A (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110518846B (zh) | 基于惯量辨识的多电机伺服系统自抗扰滑模速度控制方法 | |
CN107612445B (zh) | 具有负载加速度反馈的随动调速系统控制方法 | |
Mao et al. | Design and implementation of continuous finite-time sliding mode control for 2-DOF inertially stabilized platform subject to multiple disturbances | |
CN111504360B (zh) | 基于时间坐标磁电编码器角度精分方法及装置 | |
CN107070341A (zh) | 基于鲁棒迭代学习控制的永磁同步电机转矩脉动抑制方法 | |
CN107577149B (zh) | 一种采用分数阶快速终端滑模控制的随动控制方法 | |
CN110212819B (zh) | 一种用于高速无刷直流电机的换相误差补偿方法 | |
CN108365787A (zh) | 一种基于内模控制的永磁同步电机调速系统及其设计方法 | |
CN112838797A (zh) | 一种基于改进指数趋近律的永磁同步电机模糊滑模控制方法 | |
CN112217436A (zh) | 一种抑制电流测量误差引起的永磁同步电机转速脉动方法 | |
CN111290263A (zh) | 一种改进的基于rbfnn及bas的pid最优化控制算法 | |
CN110569767B (zh) | 运动参数识别方法、系统、设备及计算机可读存储介质 | |
CN113890424B (zh) | 基于参数辨识的永磁同步电机速度环pi控制器调谐方法 | |
CN106817054B (zh) | 一种基于参数辨识的机械弹性储能用pmsg控制方法 | |
CN114955011A (zh) | 一种dgvscmg飞轮模式下框架系统定角控制方法 | |
CN110492814A (zh) | 粒子群算法优化滑膜变结构永磁同步电机控制参数的方法 | |
CN110277941B (zh) | 永磁电机、转子位置控制方法、装置 | |
CN107395080A (zh) | 基于级联非奇异终端滑模观测器的无速度传感器转矩控制系统及方法 | |
CN115566954B (zh) | 嵌入式电机调速控制补偿方法及系统 | |
CN116015133A (zh) | 一种基于有限时间扰动观测器的永磁同步电机滑模控制方法 | |
CN106849803B (zh) | 基于均匀分布边缘粒子滤波永磁同步电动机转速估计方法 | |
CN113067506B (zh) | 内模型等价输入干扰的永磁同步电机周期扰动抑制方法 | |
CN110048653B (zh) | 无轴承异步电机的滑模变结构mras转速辨识方法 | |
CN104734591B (zh) | 汽车电动转向电机磁场定向控制的串级系统稳定调速方法 | |
CN113890451A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210105 |
|
CF01 | Termination of patent right due to non-payment of annual fee |