CN106198020A - 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 - Google Patents
基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 Download PDFInfo
- Publication number
- CN106198020A CN106198020A CN201610530259.6A CN201610530259A CN106198020A CN 106198020 A CN106198020 A CN 106198020A CN 201610530259 A CN201610530259 A CN 201610530259A CN 106198020 A CN106198020 A CN 106198020A
- Authority
- CN
- China
- Prior art keywords
- matrix
- fault
- subspace
- formula
- wind turbines
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了基于子空间和模糊C均值聚类的风电机组轴承故障诊断法,利用随机子空间的线性状态空间模型提取出反映故障本质信息的特征向量;将得到的每个特征向量作为一个数据样本,所有的特征向量组成总的数据样本,然后利用FCM对数据样本进行分类,得到故障的分类。本发明的有益效果是能够有效地识别风电机组轴承不同类型故障。
Description
技术领域
本发明属于风电技术领域,涉及基于子空间和模糊C均值聚类的风电机组轴承故障诊断法。
背景技术
由于风场大部分位于环境复杂恶劣的地区,经常受极端天气的影响。随着机组累计运行时间的增加,机组部件不断老化,极易出现故障。风电机组上的主轴、偏航、变桨、发电机、齿轮箱等许多部位都装配有轴承,轴承故障在机组故障中占有很高的比例。为了减少风电机组的停机时间,降低机组的维修费用,对风电机组的重要轴承部件进行状态监测很有必要。目前的轴承故障特征提取方法主要有时域分析法和频域分析法。时域分析法通过计算有效值、裕度因数、峭度系数、峰值等时域参数,与正常情况下的参数指标进行对比来判断轴承是否发生故障;频域分析法一般计算振动信号的幅值谱、功率谱及阶次包络谱等,根据频谱图中幅值突出的频率成分提取轴承的故障特征频率,实现轴承的故障诊断。由于风电机组轴承的运行条件恶劣,受力复杂,并经常受到交变载荷、冲击载荷的影响,使得轴承成为风电机组常见的易发生故障的部件之一。目前振动信号分析技术是旋转机械设备状态监测与故障诊断的主要手段,但风电机组与其他大型旋转机械设备存在重大区别,风机在运行过程中转速是变化的,其振动信号往往表现出非平稳的特性,使得传统的故障诊断方法难以实现故障状态的精确判别。风力发电***越来越复杂,所包含的零部件件也越来越多,使得其振动信号频率成分十分复杂,很难检测出其故障特征信息。实际风场中的风电机组振动信号一般为非平稳信号,传统的时域分析方法和频域分析方法只能够处理平稳信号,但是对属于非平稳信号的轴承振动信号则几乎无能为力,无法有效地诊断实际风电机组发生的故障。
发明内容
本发明的目的在于提供基于子空间和模糊C均值聚类的风电机组轴承故障诊断法,解决了风机在运行过程中转速是变化的,其振动信号往往表现出非平稳的特性,使得传统的目前的振动信号分析技术诊断故障方法难以实现故障状态的精确判别的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:利用随机子空间的线性状态空间模型提取出反映故障本质信息的特征向量;
其中,X是状态量,Y是采集到的振动信号,w是***噪声,v是测量噪声,A是***矩阵,C是输出矩阵;
步骤2:将得到的每个特征向量作为一个数据样本,所有的特征向量组成总的数据样本,然后利用FCM对数据样本进行分类,得到故障的分类。
进一步,所述步骤1中提取特征向量的步骤如下:
1)正交投影
首先,定义由量测量组成的分块Hankel矩阵Y:
其中,{yk}k=1,…,N+i+j为***的(N+i+j)个测量输出,Yp和Yf均为分块Hankel矩阵,p表示历史测量值,f表示未来测量值;i和j+1分别为历史量测量和未来测量形成分块矩阵的行数,min{i,j+1}≥n,矩阵Yp的维数为i×N,矩阵Yf的维数为(j+1)×N,N>>max{i,j+1};
对Y阵进行重新分块,如式(2)所示
其中,为Yf的第一行移到了Yp的末行之后的矩阵,为Yf去掉了第一行后的矩阵;对于式(1),将Yf正交投影到Yp的空间上,则正交投影Pm定义为
其中,为求Moore-Penrose逆,为矩阵Yp的逆矩阵;同样对于式(2),将正交投影到的空间上,则有正交投影Pm-1为
2)奇异值分解
对投影矩阵Pm进行奇异值分解,则有
其中,U1=[u1…un]和V1=[v1…vn]为酋矩阵,S1=diag{σ1…σn}为对角矩阵,σi是S1的第i个奇异值,U0、V0和S0为零矩阵;根据式(5)所示的奇异值分解,计算得到状态量Xk,对投影矩阵Pm-1进行奇异值分解,计算得到状态量Xk+1;
3)将计算出的状态量Xk和Xk+1,代入随机状态空间模型:
利用最小二乘法,估计出***矩阵A,***矩阵A的特征值[x1,x2,…,xn]即为随机子空间方法提取出的故障特征向量。
进一步,所述步骤2中对特征向量组成总的数据样本进行分类的方法如下:
步骤1:确定聚类中心个数K和加权指数m,初始化分类矩阵B,迭代次数l=0;
步骤2:根据下式计算所有样本[x1,x2,…,xn]的聚类中心Z=[z1,z2,…,zk]:
式中,bij为第j个故障特征向量对第i个聚类中心的隶属度,m为加权指数,n为故障特征向量的总个数,xj为第j个故障特征向量;
步骤3:根据下式更新分类矩阵B:
式中,dij为第j个故障特征向量对第i个聚类中心的欧式距离,dej为第j个故障特征向量对第e个聚类中心的欧式距离;
步骤4:对于给定判别精度ε,如果两次迭代的分类矩阵之差小于ε,则停止迭代,否则设置l=l+1,返回步骤2,直至满足条件。
本发明的有益效果是能够有效地识别风电机组轴承不同类型故障。
附图说明
图1是本发明轴承故障诊断方法流程示意图;
图2是四种状态下的滚动轴承振动信号时域波形图;
图3是随机子空间方法提取的特征向量图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
1子空间算法提取故障特征
随机子空间的线性状态空间模型描述如下:
其中,X是状态量,一般是不可测量,没有实际的物理意义,只是为了便于描述这个***动态的数学对象;Y是观测量,在这里就是采集到的振动信号;w是***噪声,由于建模不精确和一些干扰造成的;v是测量噪声。这些量均为相应维数的列向量。A是***矩阵;C是输出矩阵。
一、利用随机子空间方法提取出反映故障本质信息的特征向量。
1、正交投影
首先,定义由量测量组成的分块Hankel矩阵Y:
其中,{yk}k=1,…,N+i+j为***的(N+i+j)个测量输出;Yp和Yf均为分块Hankel矩阵,p表示历史测量值,f表示未来测量值;i和j+1分别为历史量测量和未来测量形成分块矩阵的行数,一般不小于***模型最大阶数,min{i,j+1}≥n;矩阵Yp的维数为i×N,矩阵Yf的维数为(j+1)×N,N>>max{i,j+1}。
为了计算投影矩阵Pm-1,从而进一步求取状态量Xk+1,对Y阵进行重新分块,如式(2)所示
其中,为Yf的第一行移到了Yp的末行之后的矩阵,为Yf去掉了第一行后的矩阵。
对于式(1),将Yf正交投影到Yp的空间上,则正交投影Pm可以定义为
其中,为求Moore-Penrose逆,为矩阵Yp的逆矩阵。
同样对于式(2),将正交投影到的空间上,则有正交投影Pm-1为
2、奇异值分解
对投影矩阵Pm进行奇异值分解,则有
其中,U1=[u1…un]和V1=[v1…vn]为酋矩阵,分别包含了n个左右奇异值向量;S1=diag{σ1…σn}为对角矩阵,σi是S1的第i个奇异值;U0、V0和S0为零矩阵。
根据式(5)所示的奇异值分解,可计算得到状态量Xk,同理,对投影矩阵Pm-1进行奇异值分解,可计算得到状态量Xk+1。
3、估计***矩阵A,得到反映故障本质信息的故障特征向量。
将计算出的状态量Xk和Xk+1,代入随机状态空间模型:再利用最小二乘法,估计出***矩阵A。
***矩阵A的特征值[x1,x2,…,xn]即为随机子空间方法提取出的故障特征向量。
二、模糊C均值聚类识别故障类型
将得到的每个特征向量xj(j=1,2,…,n)作为一个数据样本,所有的特征向量[x1,x2,…,xn]组成总的数据样本,然后利用FCM对数据样本进行分类。
模糊C均值聚类(Fuzzy C Means Clustering,FCM)是一种基于目标函数的聚类算法,以极小化所有数据点与各聚类中心的欧式距离及模糊隶属度的加权和为目标,不断修正聚类中心和分类矩阵到符合终止准则,将具有类似特征的数据样本聚为一类,非常适用于风电机组轴承的故障类型识别,本发明将子空间方法提取的故障特征向量输入到FCM中进行模式识别,具体步骤如下:
步骤1:确定聚类中心个数K和加权指数m,初始化分类矩阵B,迭代次数l=0。
步骤2:根据下式计算所有样本[x1,x2,…,xn]的聚类中心Z=[z1,z2,…,zk]:
式中,bij为第j个故障特征向量对第i个聚类中心的隶属度,m为加权指数,n为故障特征向量的总个数,xj为第j个故障特征向量。
步骤3:根据下式更新分类矩阵B:
式中,dij为第j个故障特征向量对第i个聚类中心的欧式距离,dej为第j个故障特征向量对第e个聚类中心的欧式距离。
步骤4:对于给定判别精度ε,如果两次迭代的分类矩阵之差小于ε,优选ε为分类矩阵B范数值的1%,则停止迭代,否则设置l=l+1,返回步骤2,直至满足条件。至此,FCM实现了对轴承故障特征向量的分类。
二、算例分析验证
本发明在风电机组实验台上进行了轴承故障试验。该***主要由风轮、主轴轴承与发电机组成,主轴轴承为滚动轴承,主要承受径向载荷,同时也可承受由于风对风轮作用产生的部分轴向载荷,在轴承座上安装加速度传感器,用来采集振动加速度信号。传感器采样率fs为10KHz,采样点数为3000。考虑到滚动轴承的局部损伤可能发生在外圈、内圈和滚动体上,在轴承的内圈、外圈与滚动体上切割沟槽来模拟轴承的局部损伤,设置的故障程度属较轻微的程度。选取正常、内圈故障、滚动体故障和外圈故障4类状态下的振动信号各40组数据作为样本。
图2为在风力发电机轴承座上测取的滚动轴承4种振动加速度信号时域图,从上到下依次为正常、内圈故障、滚动体故障和外圈故障。图3为随机子空间方法提取得到的特征向量,从图中可见不同状态的特征具有不同的分布。最后,将提取得到的特征向量输入到FCM中进行分类,得到识别的轴承故障类型。结果表明,本发明所述方法能够有效地识别风电机组轴承不同类型故障。
本发明利用随机子空间方法对风电机组轴承的时域振动信号进行建模,估计出***矩阵A并计算A的特征值,将其作为轴承的故障特征向量。将子空间方法提取出的故障特征向量作为输入量,输入到模糊C均值聚类方法中进行模式识别,对不同状态的特征进行分类。
本发明的优点在于:将子空间方法应用于风机轴承的故障特征提取,采用的是振动数据,提高了判断准确率。本发明采用的振动信号分析法是一种有效的状态检测方法,特别适用于轴承等旋转机械设备。风电机组的轴承振动信号大都是非平稳信号,传统的分析方法无法取得好的效果。本发明的随机子空间方法是一种分析非平稳振动信号的有效方法,直接在时域里分析数据建立模型,可挖掘出反映故障本质信息的特征。模糊C均值聚类是应用最为广泛的一种分类方法,受噪声的影响小,能够准确稳定地识别轴承不同类型的故障。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (3)
1.基于子空间和模糊C均值聚类的风电机组轴承故障诊断法,其特征在于按照以下步骤进行:
步骤1:利用随机子空间的线性状态空间模型提取出反映故障本质信息的特征向量;
其中,X是状态量,Y是采集到的振动信号,w是***噪声,v是测量噪声,A是***矩阵,C是输出矩阵;
步骤2:将得到的每个特征向量作为一个数据样本,所有的特征向量组成总的数据样本,然后利用FCM对数据样本进行分类,得到故障的分类。
2.按照权利要求1所述基于子空间和模糊C均值聚类的风电机组轴承故障诊断法,其特征在于:所述步骤1中提取特征向量的步骤如下:
1)正交投影
首先,定义由量测量组成的分块Hankel矩阵Y:
其中,{yk}k=1,…,N+i+j为***的(N+i+j)个测量输出,Yp和Yf均为分块Hankel矩阵,p表示历史测量值,f表示未来测量值;i和j+1分别为历史量测量和未来测量形成分块矩阵的行数,min{i,j+1}≥n,矩阵Yp的维数为i×N,矩阵Yf的维数为(j+1)×N,N>>max{i,j+1};
对Y阵进行重新分块,如式(2)所示
其中,为Yf的第一行移到了Yp的末行之后的矩阵,为Yf去掉了第一行后的矩阵;对于式(1),将Yf正交投影到Yp的空间上,则正交投影Pm定义为
其中,为求Moore-Penrose逆,为矩阵Yp的逆矩阵;同样对于式(2),将正交投影到的空间上,则有正交投影Pm-1为
2)奇异值分解
对投影矩阵Pm进行奇异值分解,则有
其中,U1=[u1…un]和V1=[v1…vn]为酋矩阵,S1=diag{σ1…σn}为对角矩阵,σi是S1的第i个奇异值,U0、V0和S0为零矩阵;根据式(5)所示的奇异值分解,计算得到状态量Xk,对投影矩阵Pm-1进行奇异值分解,计算得到状态量Xk+1;
3)将计算出的状态量Xk和Xk+1,代入随机状态空间模型:
利用最小二乘法,估计出***矩阵A,***矩阵A的特征值[x1,x2,…,xn]即为随机子空间方法提取出的故障特征向量。
3.按照权利要求1所述基于子空间和模糊C均值聚类的风电机组轴承故障诊断法,其特征在于:所述步骤2中对特征向量组成总的数据样本进行分类的方法如下:
步骤1:确定聚类中心个数K和加权指数m,初始化分类矩阵B,迭代次数l=0;
步骤2:根据下式计算所有样本[x1,x2,…,xn]的聚类中心Z=[z1,z2,…,zk]:
式中,bij为第j个故障特征向量对第i个聚类中心的隶属度,m为加权指数,n为故障特征向量的总个数,xj为第j个故障特征向量;
步骤3:根据下式更新分类矩阵B:
式中,dij为第j个故障特征向量对第i个聚类中心的欧式距离,dej为第j个故障特征向量对第e个聚类中心的欧式距离;
步骤4:对于给定判别精度ε,如果两次迭代的分类矩阵之差小于ε,则停止迭代,否则设置l=l+1,返回步骤2,直至满足条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610530259.6A CN106198020A (zh) | 2016-07-07 | 2016-07-07 | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610530259.6A CN106198020A (zh) | 2016-07-07 | 2016-07-07 | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106198020A true CN106198020A (zh) | 2016-12-07 |
Family
ID=57472399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610530259.6A Pending CN106198020A (zh) | 2016-07-07 | 2016-07-07 | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106198020A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106885697A (zh) * | 2017-03-17 | 2017-06-23 | 华东交通大学 | 基于fcm‑hmm的滚动轴承的性能退化评估方法 |
CN108827458A (zh) * | 2018-06-15 | 2018-11-16 | 西安交通大学 | 一种风力发电机叶片固有频率在线识别方法 |
CN110057588A (zh) * | 2019-05-09 | 2019-07-26 | 山东大学 | 基于奇异值与图论特征融合的轴承早期故障检测与诊断方法及*** |
CN111458141A (zh) * | 2020-04-10 | 2020-07-28 | 中国工程物理研究院机械制造工艺研究所 | 基于经验模态分解和奇异值分解的振源数估计方法 |
CN111638449A (zh) * | 2020-06-03 | 2020-09-08 | 国网北京市电力公司 | 一种配电自动化开关故障诊断方法、设备及可读存储介质 |
CN111898690A (zh) * | 2020-08-05 | 2020-11-06 | 山东大学 | 一种电力变压器故障分类方法及*** |
CN112561306A (zh) * | 2020-12-11 | 2021-03-26 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
CN113447813A (zh) * | 2020-09-03 | 2021-09-28 | 鲁能集团有限公司 | 海上风力发电机组的故障诊断方法及设备 |
US11539317B2 (en) | 2021-04-05 | 2022-12-27 | General Electric Renovables Espana, S.L. | System and method for detecting degradation in wind turbine generator bearings |
-
2016
- 2016-07-07 CN CN201610530259.6A patent/CN106198020A/zh active Pending
Non-Patent Citations (6)
Title |
---|
刘长良 等: ""基于变分模态分解和模糊C均值聚类的滚动轴承故障诊断"", 《中国电机工程学报》 * |
张小田 等: ""基于状态监测的风电机组主轴承早期故障预测方法"", 《广东电力》 * |
王军辉 等: ""基于EEMD和模糊C均值聚类的风电机组齿轮箱故障诊断"", 《太阳能学报》 * |
赵洪山 等: ""基于子空间方法的风机齿轮箱故障预测算法"", 《电力自动化设备》 * |
赵洪山 等: ""基于最大信噪比的风电机组主轴承的故障特征提取"", 《可再生能源》 * |
赵洪山 等: ""基于确定性随机子空间方法的风机齿轮箱故障预警"", 《可再生能源》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106885697A (zh) * | 2017-03-17 | 2017-06-23 | 华东交通大学 | 基于fcm‑hmm的滚动轴承的性能退化评估方法 |
CN106885697B (zh) * | 2017-03-17 | 2019-09-20 | 华东交通大学 | 基于fcm-hmm的滚动轴承的性能退化评估方法 |
CN108827458A (zh) * | 2018-06-15 | 2018-11-16 | 西安交通大学 | 一种风力发电机叶片固有频率在线识别方法 |
CN110057588A (zh) * | 2019-05-09 | 2019-07-26 | 山东大学 | 基于奇异值与图论特征融合的轴承早期故障检测与诊断方法及*** |
CN111458141A (zh) * | 2020-04-10 | 2020-07-28 | 中国工程物理研究院机械制造工艺研究所 | 基于经验模态分解和奇异值分解的振源数估计方法 |
CN111638449B (zh) * | 2020-06-03 | 2022-08-05 | 国网北京市电力公司 | 一种配电自动化开关故障诊断方法、设备及可读存储介质 |
CN111638449A (zh) * | 2020-06-03 | 2020-09-08 | 国网北京市电力公司 | 一种配电自动化开关故障诊断方法、设备及可读存储介质 |
CN111898690A (zh) * | 2020-08-05 | 2020-11-06 | 山东大学 | 一种电力变压器故障分类方法及*** |
CN111898690B (zh) * | 2020-08-05 | 2022-11-18 | 山东大学 | 一种电力变压器故障分类方法及*** |
CN113447813A (zh) * | 2020-09-03 | 2021-09-28 | 鲁能集团有限公司 | 海上风力发电机组的故障诊断方法及设备 |
CN113447813B (zh) * | 2020-09-03 | 2022-09-13 | 中国绿发投资集团有限公司 | 海上风力发电机组的故障诊断方法及设备 |
CN112561306A (zh) * | 2020-12-11 | 2021-03-26 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
CN112561306B (zh) * | 2020-12-11 | 2023-12-08 | 领伟创新智能***(浙江)有限公司 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
US11539317B2 (en) | 2021-04-05 | 2022-12-27 | General Electric Renovables Espana, S.L. | System and method for detecting degradation in wind turbine generator bearings |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106198020A (zh) | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 | |
CN110410282B (zh) | 基于som-mqe和sfcm的风电机组健康状态在线监测及故障诊断方法 | |
WO2021135630A1 (zh) | 基于grcmse与流形学习的滚动轴承故障诊断方法 | |
CN107003663B (zh) | 具有活动部分的装置的监视 | |
CN110376522B (zh) | 一种数据融合的深度学习网络的电机故障诊断方法 | |
Xiong et al. | Low-speed rolling bearing fault diagnosis based on EMD denoising and parameter estimate with alpha stable distribution | |
CN111562108A (zh) | 一种基于cnn和fcmc的滚动轴承智能故障诊断方法 | |
CN112257530B (zh) | 基于盲信号分离和支持向量机的滚动轴承故障诊断方法 | |
CN106644484A (zh) | Eemd与邻域粗糙集结合的涡桨发动机转子***故障诊断方法 | |
CN102768115A (zh) | 一种风电机组齿轮箱健康状态实时动态监控方法 | |
CN108092623B (zh) | 一种光伏阵列多传感器故障检测定位方法 | |
CN106650122B (zh) | 一种设备变工况运行风险评估方法 | |
CN111680875B (zh) | 基于概率基线模型的无人机状态风险模糊综合评价方法 | |
CN102944416A (zh) | 基于多传感器信号融合技术的风电机组叶片故障诊断方法 | |
CN105891546A (zh) | 基于大数据的风电机组偏航***中风向标故障诊断的方法 | |
CN109406147B (zh) | 一种变速工况下的列车轴承轨边声学诊断方法 | |
CN111678699B (zh) | 一种面向滚动轴承早期故障监测与诊断方法及*** | |
CN111964909A (zh) | 滚动轴承运行状态检测方法、故障诊断方法及*** | |
CN112507479B (zh) | 一种基于流形学习和softmax的石油钻机健康状态评估方法 | |
CN106127184A (zh) | 一种风机齿轮箱故障诊断方法 | |
CN111076934A (zh) | 一种基于s变换诊断轴承潜在故障的方法 | |
CN109494716A (zh) | 基于Bootstrap的风电输出功率置信区间预测方法 | |
CN111474476B (zh) | 一种电机故障预测方法 | |
CN116337447A (zh) | 一种非平稳工况下轨道车辆轮对轴承故障诊断方法及设备 | |
CN116304551A (zh) | 一种基于bcb模型的电机轴承故障诊断及特征提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161207 |