CN110333443B - 感应电机定子绕组温升测试方法 - Google Patents

感应电机定子绕组温升测试方法 Download PDF

Info

Publication number
CN110333443B
CN110333443B CN201910508381.7A CN201910508381A CN110333443B CN 110333443 B CN110333443 B CN 110333443B CN 201910508381 A CN201910508381 A CN 201910508381A CN 110333443 B CN110333443 B CN 110333443B
Authority
CN
China
Prior art keywords
stator winding
temperature
temperature rise
motor
thermal
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
CN201910508381.7A
Other languages
English (en)
Other versions
CN110333443A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201910508381.7A priority Critical patent/CN110333443B/zh
Publication of CN110333443A publication Critical patent/CN110333443A/zh
Application granted granted Critical
Publication of CN110333443B publication Critical patent/CN110333443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines
    • G01R31/343Testing dynamo-electric machines in operation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Abstract

本发明公开了一种感应电定子绕组温升测试方法。该方法基于电机瞬态热路模型,结合参数估计方法,利用定子绕组瞬态温升实测数据得到其温升曲线,从而获得定子绕组的温升限值。本发明具有测试时间明显减少的优势,且可获得定子绕组稳态温度与载流的关系曲线,仅需数个实验即可获得电机定子绕组的温升曲线。同时,根据曲线可得到绕组的温升限值及长时间可靠运行的最大电流值。

Description

感应电机定子绕组温升测试方法
技术领域
本发明属于电机控制技术,具体为一种感应电定子绕组温升测试方法。
背景技术
感应电机定子绕组运行温度是电机安全运行的重要参数,如果定子绕组工作温度在一定时间内超过其绝缘材料的极限温度,将严重影响电机性能,降低电机使用寿命,严重时会导致电机损毁。因此,实际中需对电机定子绕组进行温升测试,得到其稳态温度及温升限值,以验证电机绝缘耐热性能,并为后续电机设计和改进提供指导依据。
实验上电机绕组温升测量方法主要有温度计、热电偶和热电阻等,但存在测量精度易受外界环境的影响、热态电阻不易测量、试验时间长、无法获得稳态温度等缺点。
理论上电机绕组的温升测量方法主要有等效热路法和数值分析法两类。等效热路法是将温度场简化为带有集总参数的热路进行计算,其准确度很大程度上取决于热路拓扑结构和热路参数。数值分析法如有限元分析(FEA)和计算流体动力学(CFD),其精度受电机的结构参数以及所用材料的物性参数影响,且材料的物性参数较难确定,另外,数值分析法尤其是有限元法计算量庞大。
发明内容
本发明的目的在于提出了一种基于瞬态热路模型的感应电机定子绕组温升测量方法。
实现本发明目的的技术方案为:一种感应电机定子绕组温升测试方法,包括以下步骤:
步骤1、根据电机定子绕组瞬态热路模型,确定定子绕组温度表达式,将此热路模型等效为一个非线性多输出多输入系统,得到状态方程;
步骤2、建立参数向量,采用参数估计方法获得定子绕组瞬态热路模型的多组热参数估计值,将多组热参数带入离散状态方程计算得到多个定子绕组温度;
步骤3、以定子绕组温度计算结果与实测结果的最小误差为目标建立定子绕组瞬态热路模型热参数的优化函数,通过对比定子绕组温度计算值与实际测量定子绕组温度的差值确定满足优化函数设定值的最佳热参数;
步骤4、根据最佳热参数确定电机定子绕组瞬态热路模型具体形式,绘制电机定子绕组的温升曲线获得定子绕组的温升限值。
优选地,步骤1确定的定子绕组温度表达式为:
Figure BDA0002091061990000021
式中,Peq为电机热损耗,Rj表示第j个定子热阻,Ci表示第i个定子热容,Ti为第i 个节点温度,Te为环境温度,n为自然数。
优选地,所述定子绕组热损耗Peq具体为:
Figure BDA0002091061990000022
其中,P0为初始温度Te下的电机热损耗。
优选地,步骤1将定子绕组热模型等效为一个非线性多输入多输出系统,得到的线性状态方程:
Figure BDA0002091061990000023
其中,T=[T1 T2... Tn]T
Figure BDA0002091061990000024
U=[Peq Te]T
Figure BDA0002091061990000025
Figure BDA0002091061990000031
矩阵A中各系数如下:
Figure BDA0002091061990000032
Figure BDA0002091061990000033
优选地,步骤2建立的参数向量为:
X=[C1 C2 ...Cn R1 R2 ... Rn+1]T
优选地,步骤3建立的定子绕组瞬态热路模型热参数的优化函数为:
Figure BDA0002091061990000034
其中,T(k)为k时刻定子绕组的计算温度,
Figure BDA0002091061990000062
为k时刻感应电机定子绕组的实测温度。
本发明与现有技术相比,其显著优点为:(1)本发明无需等待绕组温升达到稳态,明显减少测试时间;(2)本发明仅需数个实验即可获得电机定子绕组的温升曲线;(3) 本发明根据温升曲线可得到绕组的温升限值及长时间可靠运行的最大电流值。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1是本发明流程图。
图2是感应电机定子绕组瞬态热路模型示意图。
图3是不同电流条件下定子绕组温升实验结果示意图。
图4是94A、108A和120A条件下定子绕组温升曲线图。
图5是定子绕组温升计算与实测值对比图。
具体实施方式
一种基于瞬态热路模型的感应电机定子绕组温升测量方法,通过选择合适的电机定子绕组瞬态热路模型,并基于控制理论,得到定子绕组温度的状态方程,通过遗传算法得到状态方程的各个参量拟合值,将参量带入离散状态方程求解定子绕组温度,以定子绕组温度计算值与实测值的误差最小值为目标优化函数,若该最小差值满足目标优化函数设定值,进而得到最佳热参数值,将具体的状态方程绘制温升特性曲线,根据温升特性曲线获得定子绕组的温升限值,具体步骤为:
步骤1、根据电机定子绕组瞬态热路模型,确定定子绕组温度表达式,将此热路模型等效为一个非线性多输出多输入系统,得到状态方程;
步骤2、建立参数向量,采用参数估计方法获得定子绕组瞬态热路模型的多组热参数估计值,将多组热参数带入离散状态方程计算得到多个定子绕组温度;
步骤3、以定子绕组温度计算结果与实测结果的最小误差为目标建立定子绕组瞬态热路模型热参数的优化函数,通过对比定子绕组温度计算值与实际测量定子绕组温度的差值确定满足优化函数设定值的最佳热参数;
步骤4、根据最佳热参数确定电机定子绕组瞬态热路模型具体形式,绘制电机定子绕组的温升曲线获得定子绕组的温升限值。
进一步的实施例中,步骤1确定的定子绕组温度表达式为:
Figure BDA0002091061990000041
式中,Peq为电机热损耗,Rj表示第j个定子热阻,Ci表示第i个定子热容,Ti为第i 个节点温度,Te为环境温度,n为自然数。
进一步的实施例中,所述定子绕组热损耗Peq具体为:
Figure BDA0002091061990000051
其中,P0为初始温度Te下的电机热损耗。
进一步的实施例中,步骤1将定子绕组热模型等效为一个非线性多输入多输出系统,得到的线性状态方程:
Figure BDA0002091061990000052
其中,T=[T1 T2... Tn]T
Figure BDA0002091061990000053
U=[Peq Te]T
Figure BDA0002091061990000054
Figure BDA0002091061990000055
矩阵A中各系数如下:
Figure BDA0002091061990000056
Figure BDA0002091061990000057
进一步的实施例中,步骤2建立的参数向量为:
X=[C1 C2 ... Cn R1 R2 ... Rn+1]T
进一步的实施例中,步骤3建立的定子绕组瞬态热路模型热参数的优化函数为:
Figure BDA0002091061990000061
其中,T(k)为k时刻定子绕组的计算温度,
Figure BDA0002091061990000062
为k时刻感应电机定子绕组的实测温度
实施例
如图1所示,本实施例针对四阶瞬态热路模型,一种感应电机定子绕组温升测试方法,包括以下步骤:
步骤1、采用图1所示瞬态热路模型,根据热电类比理论,得到绕组温度表达式:
Figure BDA0002091061990000063
其中,Peq为电机热损耗,R1为定子绕组绝缘层热阻,R2为定子槽绝缘层热阻,R3为忽略齿部的定子叠片热阻,R4为定子叠片与机壳的热阻,R5为气隙热阻,C1为定子绕组热容,C2为定子槽绝缘层热容,C3为定子齿热容,C4为定子轭热容,T1为电机定子绕组温度,T2为定子槽绝缘层温度,T3为定子齿温度,T4为定子轭温度,Te为环境温度。在转子堵转、电机通直流电的情况下,转子温度、机壳温度和环境温度近似相等,均为 Te,进一步地,所述定子绕组热损耗Peq具体为:
Figure BDA0002091061990000064
其中,P0为初始温度Te下的电机热损耗。
将温度T作为状态变量,Peq、Te作为输入量,则式(1)可转换为如下所示的线性状态方程:
Figure BDA0002091061990000071
其中,T=[T1 T2 T3 T4]T
Figure BDA0002091061990000072
u=[Peq Te]T
Figure BDA0002091061990000073
Figure BDA0002091061990000074
矩阵A中各系数如下:
Figure BDA0002091061990000075
Figure BDA0002091061990000076
步骤2、建立参数向量:
X=[C1 C2 C3 C4 R1 R2 R3 R4 R5]T
根据测量所得温度,通过遗传算法求解得到多组电机定子绕组温升模型中的未知参数向量X;
利用遗传算法对待辨识参数进行估计时,遗传算法的操作参数如表1所示。
Figure BDA0002091061990000077
Figure BDA0002091061990000081
表1
通过前期的温度测量,将时间作为自变量,温度测量值作为因变量,利用遗传算法得到温度随时间变化曲线方程的各项参数值,即获取多组不同的电机定子绕组温升模型中的未知参数向量X。
将热参数带入离散状态方程求得定子绕组温度,具体方法为:
k时刻的定子各部分温度为Ti (k),环境温度为Te (k),电机热损耗为Peq (k),其中i=1,2, 3,4;k∈N。
对(3)式按时间进行离散,令其中:T(k)=[T1 (k) T2 (k) T3 (k) T4 (k)]T; U(k)=[Peq (k) Te (k)]T
得到迭代方程:T(k)=CT(k-1)+DU(k-1) (6)
其中,C=ΔtA+E,D=ΔtB
其中,Δt为离散时间,E为单位矩阵。
步骤3、以热路模型的定子绕组温度计算结果与实测结果误差最小为目标,建立目标优化函数如下:
Figure BDA0002091061990000082
其中,X为待辨识参数向量表示为:
X=[C1 C2 C3 C4 R1 R2 R3 R4 R5]T
Figure BDA0002091061990000091
为k时刻感应电机定子绕组的实测温度。
根据目标优化函数式(7),求出多组定子绕组温度计算值与实际温度测量值的误差绝对值的最小值,当该最小误差温度满足给定的设定值时即可确定所求热参量,进而得到电机定子绕组的温升特性。若不满足,则重新收集温度数据进行步骤2获得新的热参数并重新进行步骤3。
步骤4、绘制电机定子绕组的温升曲线,根据温升曲线获得定子绕组的温升限值。
对本实施例进行验证,具体过程如下:
通过验证平台获得电机温升数据。验证平台由鼠笼式感应电机、直流电源、测温模块以及PC机组成。其中,感应电机参数如表2所示;测温模块采用NTC503J3950温度传感器,其B值精度为1%,阻值精度为1%;直流电源为电机提供72V直流电;温度、电流和电压数据通过232总线发送至PC机。
Figure BDA0002091061990000092
表2
将温度传感器探头置于电机机盖内部(所测温度为定子绕组端部温度,可近似为电机定子侧温度)。本验证在恒温恒湿、自然对流条件下进行,环境温度为30℃,测量时间间隔为0.3s。对电机定子绕组加载不同大小的直流电,测量定子绕组的温升。采用直流电进行温升测试,电机损耗仅在定子绕组内,因而无需考虑定子铁耗和转子损耗。
按电流大小不同共进行6组,分别为80A、94A、100A、108A、120A和134A。定子绕组的最高温度为110℃,未超过电机绝缘等级的要求。验证结果如图3所示。
进一步地,对估计参数进行实验验证与分析。
分别选取94A、108A和120A载流条件下30℃至80℃区间的数据进行参数估计。采用遗传算法操作参数,计算得到定子绕组热路模型热参数,如表3所示,可得到电机定子绕组温升的具体状态方程。
Figure BDA0002091061990000101
表3
图4给出了定子绕组在94A、108A和120A载流条件下从初始温度30℃上升到80℃的实测曲线和计算曲线。由图可见,定子绕组温度还未达稳态,计算值与实测值较为吻合,电机定子绕组通流后,其温度迅速升高,且加载电流越大,温升速率越高。
表4给出了计算值与实测值的计算误差,最大计算误差小于2%,平均计算误差小于 1%。
Figure BDA0002091061990000102
表4
进一步地,选取80℃至110℃区间定子绕组温度的数据来验证参数估计方法的有效性。图5为电机定子绕组温升的计算与实测值对比图,其中虚线为计算值,实线为实测值,两者基本吻合。
表5给出了从80℃至110℃区间的计算误差,最大计算误差在1%以内,平均计算误差在1%以内,表明本方法所用电机定子绕组温升模型参数估计方法是有效性的。
Figure BDA0002091061990000111
表5
因此,本发明可通过加载较大的负载电流,来减少测试时间。本发明还可获得传统测试方法由于电机绝缘等级的限制而无法得到的超温条件下定子绕组稳态温升。

Claims (6)

1.一种感应电机定子绕组温升测试方法,其特征在于,包括以下步骤:
步骤1、根据电机定子绕组瞬态热路模型,确定定子绕组温度表达式,将此热路模型等效为一个非线性多输出多输入系统,得到状态方程;
步骤2、建立参数向量,采用参数估计方法获得定子绕组瞬态热路模型的多组热参数估计值,将多组热参数带入离散状态方程计算得到多个定子绕组温度;
步骤3、以定子绕组温度计算结果与实测结果的最小误差为目标建立定子绕组瞬态热路模型热参数的优化函数,通过对比定子绕组温度计算值与实际测量定子绕组温度的差值确定满足优化函数设定值的最佳热参数;
步骤4、根据最佳热参数确定电机定子绕组瞬态热路模型具体形式,绘制电机定子绕组的温升曲线获得定子绕组的温升限值。
2.根据权利要求1所述的感应电机定子绕组温升测试方法,其特征在于,步骤1确定的定子绕组温度表达式为:
Figure FDA0002091061980000011
式中,Peq为电机热损耗,Rj表示第j个定子热阻,Ci表示第i个定子热容,Ti为第i个节点温度,Te为环境温度,n为自然数。
3.根据权利要求2所述的感应电机定子绕组温升测试方法,其特征在于,所述定子绕组热损耗Peq具体为:
Figure FDA0002091061980000012
其中,P0为初始温度Te下的电机热损耗。
4.根据权利要求1所述的感应电机定子绕组温升测试方法,其特征在于,步骤1将定子绕组热模型等效为一个非线性多输入多输出系统,得到的线性状态方程:
Figure FDA0002091061980000013
其中,T=[T1 T2 … Tn]T
Figure FDA0002091061980000021
U=[Peq Te]T
Figure FDA0002091061980000022
Figure FDA0002091061980000023
矩阵A中各系数如下:
Figure FDA0002091061980000024
Figure FDA0002091061980000025
5.根据权利要求1所述的感应电机定子绕组温升测试方法,其特征在于,步骤2建立的参数向量为:
X=[C1 C2 … Cn R1 R2 … Rn+1]T
6.根据权利要求1所述的感应电机定子绕组温升测试方法,其特征在于,步骤3建立的定子绕组瞬态热路模型热参数的优化函数为:
Figure FDA0002091061980000026
其中,T(k)为k时刻定子绕组的计算温度,
Figure FDA0002091061980000027
为k时刻感应电机定子绕组的实测温度。
CN201910508381.7A 2019-06-12 2019-06-12 感应电机定子绕组温升测试方法 Active CN110333443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910508381.7A CN110333443B (zh) 2019-06-12 2019-06-12 感应电机定子绕组温升测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910508381.7A CN110333443B (zh) 2019-06-12 2019-06-12 感应电机定子绕组温升测试方法

Publications (2)

Publication Number Publication Date
CN110333443A CN110333443A (zh) 2019-10-15
CN110333443B true CN110333443B (zh) 2021-07-06

Family

ID=68140403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910508381.7A Active CN110333443B (zh) 2019-06-12 2019-06-12 感应电机定子绕组温升测试方法

Country Status (1)

Country Link
CN (1) CN110333443B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112213638B (zh) * 2020-10-28 2021-12-14 瑞声新能源发展(常州)有限公司科教城分公司 一种电机线圈的发热状态监测方法、相关设备和介质
CN113720495A (zh) * 2021-08-26 2021-11-30 南京理工大学 一种电机转子实时温度估计方法
CN114647263A (zh) * 2022-03-11 2022-06-21 杭州新剑机器人技术股份有限公司 温度控制方法、系统以及机器人的驱动系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1234640A (zh) * 1997-12-19 1999-11-10 施耐德电器公司 模拟电动机内温升的电子装置
CN101620013A (zh) * 2008-07-03 2010-01-06 艾默生电气公司 校准电动机控制电路以改进电动机温度测量的方法和系统
TW201209382A (en) * 2010-08-25 2012-03-01 Nat Univ Chung Cheng An error compensation apparatus for the built-in motor spindle
CN103383433A (zh) * 2013-07-03 2013-11-06 中国人民解放军海军工程大学 船舶发电机定子铁心的状态监测和故障预警方法
CN104978477A (zh) * 2014-04-14 2015-10-14 现代自动车株式会社 用于估计电动机的转子温度的系统和方法
CN105181173A (zh) * 2015-09-02 2015-12-23 南车株洲电力机车研究所有限公司 一种监测电机绕组温升的方法及装置
CN107966659A (zh) * 2017-12-07 2018-04-27 重庆车辆检测研究院有限公司 永磁同步电机中定子绕组温升在线检测方法及检测系统
CN108562854A (zh) * 2018-04-08 2018-09-21 华中科技大学 一种电机异常状态在线预警方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10298168B2 (en) * 2016-07-07 2019-05-21 Schweitzer Engineering Laboratories, Inc. Slip-dependent motor model

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1234640A (zh) * 1997-12-19 1999-11-10 施耐德电器公司 模拟电动机内温升的电子装置
CN101620013A (zh) * 2008-07-03 2010-01-06 艾默生电气公司 校准电动机控制电路以改进电动机温度测量的方法和系统
TW201209382A (en) * 2010-08-25 2012-03-01 Nat Univ Chung Cheng An error compensation apparatus for the built-in motor spindle
CN103383433A (zh) * 2013-07-03 2013-11-06 中国人民解放军海军工程大学 船舶发电机定子铁心的状态监测和故障预警方法
CN104978477A (zh) * 2014-04-14 2015-10-14 现代自动车株式会社 用于估计电动机的转子温度的系统和方法
CN105181173A (zh) * 2015-09-02 2015-12-23 南车株洲电力机车研究所有限公司 一种监测电机绕组温升的方法及装置
CN107966659A (zh) * 2017-12-07 2018-04-27 重庆车辆检测研究院有限公司 永磁同步电机中定子绕组温升在线检测方法及检测系统
CN108562854A (zh) * 2018-04-08 2018-09-21 华中科技大学 一种电机异常状态在线预警方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
stator-winding thermal models for short-time thermal transients:definition and validation;Aldo Boglietti et.;《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》;20160531;第63卷(第5期);正文第2713至2720页 *
基于效率及温升的轴向磁通永磁电机优化设计;陈晨 等;《中国电机工程学报》;20160320;第36卷(第6期);正文第1686至1693页 *

Also Published As

Publication number Publication date
CN110333443A (zh) 2019-10-15

Similar Documents

Publication Publication Date Title
CN110333443B (zh) 感应电机定子绕组温升测试方法
Huber et al. A low-order thermal model for monitoring critical temperatures in permanent magnet synchronous motors
CN105160202B (zh) 确定油浸式变压器绕组热点温度和顶层油温的方法
Huber et al. Monitoring critical temperatures in permanent magnet synchronous motors using low-order thermal models
CN107391822B (zh) 一种基于自适应时间步长的瞬变电磁-温度场耦合计算方法
CN109510544B (zh) 基于热阻网络模型的电机参数估计方法及系统
CN106849011A (zh) 一种伺服电机过热保护方法
CN110991123A (zh) 基于l-m算法优化的变压器绕组热点瞬态温度计算方法
CN101132129A (zh) 一种基于热积累的交流电机过载反时限保护方法
CN102880747A (zh) 一种光电复合海缆温度场建模计算分析方法
Zhang Online thermal monitoring models for induction machines
Shi et al. Electromagnetic-thermal coupled simulation under various fault conditions of a triple redundant 9-phase PMASynRM
CN110431392B (zh) 用于确定循环泵中的输送介质温度的方法以及循环泵
Zhu et al. Coupled electromagnetic-thermal-fluidic analysis of permanent magnet synchronous machines with a modified model
Josue et al. Transformer hot-spot temperature estimation for short-time dynamic loading
Lystianingrum et al. On estimating instantaneous temperature of a supercapacitor string using an observer based on experimentally validated lumped thermal model
CN102937670A (zh) 一种基于Elman神经网络的感应电机转子电阻参数辨识方法
CN108038268B (zh) 一种廊管内gil运行温度评估方法
Saheba et al. Virtual thermal sensing for electric machines
Radakovic et al. Loading of transformers in conditions of controlled cooling system
Wang et al. Coupled electromagnetic and thermal analysis of permanent magnet rectifier generator based on LPTN
Radakovic et al. Thermal model of oil power transformers with a tap changer
CN112234911A (zh) 一种永磁电机转子永磁体温度的实时监测方法及模型
CN111274669A (zh) 一种具备热参数校准功能的爪极发电机瞬态温度监测方法
CN113489423B (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