CN114692695B - 基于多指标融合的故障特征选择及分类方法研究 - Google Patents
基于多指标融合的故障特征选择及分类方法研究 Download PDFInfo
- Publication number
- CN114692695B CN114692695B CN202210379191.1A CN202210379191A CN114692695B CN 114692695 B CN114692695 B CN 114692695B CN 202210379191 A CN202210379191 A CN 202210379191A CN 114692695 B CN114692695 B CN 114692695B
- Authority
- CN
- China
- Prior art keywords
- component
- decomposition
- fault
- decomposed
- index
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000004927 fusion Effects 0.000 title claims abstract description 11
- 238000011160 research Methods 0.000 title description 5
- 238000005096 rolling process Methods 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 33
- 238000012549 training Methods 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 6
- 230000000295 complement effect Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 description 4
- 238000007637 random forest analysis Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 101100001677 Emericella variicolor andL gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
基于多指标融合的故障特征选择及分类方法研究,首先将原始信号通过改进CELMD方法分解为多个PF分量和单个残余分量;然后融合峭度、相关系数和K‑L散度三个指标提取最适合的PF分量进行信号重构;最后提取重构信号的多域特征和熵值特征,输入LightGBM分类器进行分类,实现滚动轴承故障的智能分类。本发明改善了端点效应,缓解了模态混叠的现象,有效减少了虚假分量对故障诊断的影响,进一步提高了故障诊断的分类准确度,并最终有效对不同故障类型进行分类。
Description
技术领域
本发明属于旋转机械装置中关键部件滚动轴承的故障诊断技术领域,具体涉及基于多指标融合的故障特征选择及分类方法研究。
背景技术
滚动轴承是旋转机械设备的关键组成部分,其性能的好与坏直接影响到整个设备的运行状况从而造成重大的经济损失,更有甚者会使操作人员的生命遭受威胁。据不完全统计,整个机械故障中有40%以上是由轴承故障问题造成的;在旋转机械中有约30%的故障问题是源于滚动轴承故障;在齿轮箱发生的故障中,轴承的故障约占到了19%;电机设备的故障当中轴承故障占比高达80%。因此,对滚动轴承及时进行检测运行状态和进行故障诊断意义重大。滚动轴承发生故障时会产生非平稳的、非高斯的信号,现有的信号分解技术可以对信号进行解调,将非平稳信号分解为多个敏感分量和单个残余分量之和,大部分信号分解技术仍然存在端点效应、模态混叠等问题,甚至出现过拟合、欠拟合等问题;此外,孤立的单一指标如相关系数、峭度或者K-L散度虽然能够作为选择某些敏感分量选择的依据,但是无法准确选择最佳的敏感分量,更有可能会丢弃某些包含特征的敏感分量,严重影响了重构信号的特征提取和故障诊断的准确性。综上所述,如何有效分解信号,缓解故障信号的端点效应、改善模态混叠的问题、避免过拟合、欠拟合等问题,并且筛选出分解之后的最佳敏感分量进行重构提取有效故障振动特征是滚动轴承故障诊断面临的一大难题。
发明内容
为了克服现有技术存在的缺点,本发明提供了基于多指标融合的故障特征选择及分类方法,利用改进的补充总体局部均值分解(Complementary Ensemble local meandecomposition,CELMD)算法将原始信号分解为多个具有物理意义的乘积函数(Productfunction,PF)和单个残余分量,其次提出多指标选择(Multi-Indicator Selection,MIS)算法选取最佳PF分量进行重构获得训练信号并提取时域、频域和时频域等多域特征和熵值特征,最后将重构信号的故障特征作为训练样本利用轻量级梯度提升机(Light GradientBoosting Machine,LightGBM)建模对轴承故障类别进行分类,高效、准确地实现滚动轴承的准确故障诊断。
为了达到上述目的,本发明所采用的技术方案为:
基于多指标融合的故障特征选择及分类方法研究,包括以下步骤:
A.获取滚动轴承在不同工况下的声发射信号,对于不同系统状态产生的信号(包括正常状态、不同故障直径下的内圈故障、滚动体故障、外圈故障),均使用本发明改进的CELMD算法进行分解,具体过程是:
a.对每种类型的原始信号x(t)成对加入所需的白噪声ni(t),即
其中,x(t)为原始信号,ni(t)为第i次加入的噪声且i∈[1,2,3…N],hi为噪声的幅值,为添加了正白噪声之后的信号,/>为添加了负白噪声之后的信号;
b.对加入正负白噪声的信号进行局部均值分解(Local Mean Decomposition,LMD)分解,分解成多个具有物理意义的乘积函数(Product function,PF),同时,修改同幅值的高斯白噪声ni(t),再次进行CELMD分解,重复N次后,将PF分量相加取平均值,则第j个PF分量为:
其中,为/>通过LMD分解并平均后的第i个PF分量,/>为/>通过LMD分解并平均后的第i个PF分量;
此时,x(t)被分解为M个PF分量和一个残余分量r(t),即
c.修改噪声幅值hi并重复以上步骤,重复N′次,则第k次分解的结果为
其中,为第k次CELMD分解的PF分量,r(k)(t)为第k次分解的残余分量;
将N’次分解的PF分量相加取平均值,即
其中,为最终分解后的第j个PF分量;
将x(t)分解为M个PF分量和残余分量,即
其中,r′(t)为分解后的残余分量;
B.采用本发明提出的多指标选择(Multi-Indicator Selection,MIS)算法对上述步骤所获的PF分量进行筛选,具体过程是:
a.首先,构建指标集L={L1,L2,L3},其中,L1为峭度,L2为相关系数,L3为Wassterstein距离的倒数;
b.将改进CELMD分解后得到的第i个PF分量分别计算指标集L中的值,从而得到mi1、mi2、mi3,要保证每一个对应的mij值符号相等,故将所有mij值进行非负处理,即
其中,j∈[1,2,3]且
c.计算每个在指标Lk的占比,并将同属于一个PF分量的占比相加得到/>即
其中,为每个PF分量中三个指标总和的占比;
d.设定固定阈值A(A为常数且满足0≤A≤1),将从大到小排序,将排序之后的从大到小进行选择,当选择的/>之和到达阈值A时,所选择的/>值对应的PF分量即为所选分量;
C.将上述MIS算法选择得到的PF分量进行重构,并将重构信号划分为训练集和测试集,再分别提取时域、频域特征和熵值特征,同时将不同故障类型打上标签作为该信号特征对应的系统状态,利用训练集训练轻量级梯度提升机(Light Gradient BoostingMachine,LightGBM)分类器,完成训练后,利用测试集评估分类效果,从而实现故障智能分类。
本发明的有益效果为:
本发明综合考虑了传统的信号分解技术对振动信号存在的缺点,改进了传统的CELMD分解技术,提出多次加入不同幅值的高斯白噪声取平均值进而降低了故障信号中噪声所带来的干扰,改善了端点效应,缓解了模态混叠的现象。综合对比了峭度、相关系数和K-L散度的优缺点,提出了MIS算法将多个指标进行融合,利用Wassterstein距离代替了K-L散度,分析了每个分量的融合指标占比,有效减少了虚假分量对故障诊断的影响,进一步提高了故障诊断的分类准确度,对重构信号提取了时域、频域等多域特征和熵域特征,构建LightGBM分类器进行故障诊断,最终有效对不同故障类型进行分类。
附图说明
图1为本发明的流程图;
图2为分类结果对比图。
具体实施方式
下面结合附图对本发明做进一步的详细描述。
参照图1,基于多指标融合的故障特征选择及分类方法研究,包括以下步骤:
A.采集滚动轴承在不同工况下的声发射信号,对于不同系统状态产生的信号(包括正常状态、不同故障直径下的内圈故障、滚动体故障、外圈故障),均使用本发明改进的CELMD算法进行分解,具体过程是:
a.对每种类型的原始信号x(t)成对加入所需的白噪声ni(t),即
其中,x(t)对应不同故障类型的原始信号,ni(t)对应第i次加入的高斯白噪声且i∈[1,2,3…N],N=10,hi对应噪声的幅值,hi依次取[0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50],为添加了正白噪声之后的信号,/>为添加了负白噪声之后的信号
b.对加入正负白噪声的信号进行LMD分解,分解成多个PF分量,同时,修改同幅值的高斯白噪声ni(t),再次进行CELMD分解,重复N次后,将PF分量相加取平均值,则第j个PF分量为:
其中,对应/>通过LMD分解并平均后的第i个PF分量,/>对应通过LMD分解并平均后的第i个PF分量;
此时,x(t)被分解为M个PF分量和一个残余分量r(t):
c.修改噪声幅值hi并重复以上步骤,重复N’次,则第k次分解的结果为
其中,为第k次CELMD分解的PF分量,r(k)(t)为第k次分解的残余分量;
将N’次分解的PF分量相加取平均值,即
其中,N’对应hi的个数对应第l次分解的第j个PF分量;
将x(t)分解为M个PF分量和残余分量,即
其中,r′(t)为残余分量,M对应最终分解之后PF的个数;
B.采用本发明提出的多指标选择(Multi-Indicator Selection,MIS)算法对上述步骤所获的PF分量进行筛选,具体过程是:
a.首先,构建指标集L={L1,L2,L3},其中,L1为峭度,L2为相关系数,L3为Wassterstein距离的倒数;
b.将改进CELMD分解后得到的第i个PF分量分别计算指标集L中的值,从而得到mi1、mi2、mi3,要保证每一个对应的mij值符号相等,故将所有mij值进行非负处理,即
其中,j∈[1,2,3]且Lj对应指标集中第j个分量;
c.计算每个在指标Lj的占比,并将同属于一个PF分量的占比相加得到/>即
其中,为每个PF分量中三个指标总和的占比;
d.设定固定阈值A(A为常数且满足0≤A≤1)A=0.85,将从大到小排序,将排序之后的/>从大到小进行选择,当选择的/>之和到达阈值A时,所选择的/>值所对应的PF分量即为所选分量;
C.将上述MIS算法选择得到的PF分量进行重构,并将重构信号划分为训练集和测试集,再分别提取时域、频域特征和熵值特征,同时将不同故障类型打上标签作为该信号特征对应的系统状态,利用训练集训练轻量级梯度提升机(Light Gradient BoostingMachine,LightGBM)分类器,完成训练后,利用测试集评估分类效果,从而实现故障智能分类;
本发明的高维混合特征向量由20个特征组成原始特征集F=(f1,f2,…,f20),共包括10个时域特征、5个频域特征和5个熵值特征。时域统计特征有均方根值(f1)、峰值(f2)、波峰系数(f3)、方根幅值(f4)、绝对平均值(f5)、偏度值(f6)、偏度因子值(f7)、峭度(f8)、形状因子(f9)、k因子(f10);频域统计特征有重心频率(f11)、均方频率(f12)、均方根频率(f13)、频率方差(f14)、频率标准差(f15);熵值特征有模糊熵(f16)、相位熵(f17)、增量熵(f18)、样本熵(f19)、余弦相似性熵(f20)
D.为了证明本发明方法的有效性,将其与不同的选择PF分量方法和不同分类方法进行比较,作进一步描述。滚动轴承数据集共有两种系统状态为正常状态和故障状态,而数据集包含3种故障直径,分别为0.1778mm、0.3556mm和0.5334mm,每种故障直径对应3种系统状态,分别为内圈故障、外圈故障和滚动体故障,所以共有10种系统状态,将每个数据集按照512个数据点进行切片进而划分训练集和数据集,使用本发明方法对对滚动轴承数据集进行故障分类,针对数据集,首先通过改进CELMD算法对原始信号进行分解,再利用MIS算法筛选出最优PF分量进行重构,最后提取特征划分数据集进行分类和诊断。对比分析了本发明方法和原始信号+RandomForest、原始信号+XGBoost、原始信号+LightGBM、correlation+RandomForest、correlation+XGBoost、correlation+LightGBM、K-L+RandomForest、K-L+XGBoost、K-L+LightGBM、MIS+RandomForest和MIS+XGBoost分类准确度,使用本发明方法可以达到99.24%的分类准确度,分别比上述方法高出24.78%、13.68%、10.24%、14.98%、9.16%、5.06%、13.47%、7.87%、4.85%、6.57%、4.53%,说明本发明方法对PF分量的选取具有稳健性,能够高效地对故障进行诊断,十二种方法分类结果对比如图2所示。
Claims (3)
1.基于多指标融合的故障特征选择及分类方法,特征在于包括以下步骤:
A.获取滚动轴承在不同工况下的声发射信号,对于不同系统状态产生的信号,包括正常状态、不同故障直径下的内圈故障、滚动体故障、外圈故障,均使用改进的补充总体局部均值分解(Complementary Ensemble local mean decomposition,CELMD)算法进行分解,具体过程是:
a.对每种类型的原始信号x(t)成对加入所需的白噪声ni(t),即
其中,x(t)为原始信号,ni(t)为第i次加入的噪声且i∈[1,2,3…N],hi为噪声的幅值,为添加了正白噪声之后的信号,/>为添加了负白噪声之后的信号;
b.对加入正负白噪声的信号进行局部均值分解(Local Mean Decomposition,LMD)分解,分解成多个具有物理意义的乘积函数(Product function,PF),同时,修改同幅值的高斯白噪声ni(t),再次进行CELMD分解,重复N次后,将PF分量相加取平均值,则第j个PF分量为:
其中,为/>通过LMD分解并平均后的第i个PF分量,/>为/>通过LMD分解并平均后的第i个PF分量;
此时,x(t)被分解为M个PF分量和一个残余分量r(t),即
c.修改噪声幅值hi并重复以上步骤,重复N’次,则第k次分解的结果为:
其中,为第k次CELMD分解的PF分量,r(k)(k)为第k次分解的残余分量;
将N’次分解的PF分量相加取平均值,即
其中,为最终分解后的第j个PF分量;
将x(t)分解为M个PF分量和残余分量,即
其中,r'(t)为分解后的残余分量;
B.采用提出的多指标选择(Multi-Indicator Selection,MIS)算法对上述步骤所获的PF分量进行筛选,具体过程是:
a.首先,构建指标集L={L1,L2,L3},其中,L1为峭度,L2为相关系数,L3为Wassterstein距离的倒数;
b.将改进CELMD分解后得到的第i个PF分量分别计算指标集L中的值,从而得到mi1、mi2、mi3,要保证每一个对应的mij值符号相等,故将所有mij值进行非负处理,即
其中,j∈[1,2,3]且
c.计算每个在指标Lj的占比,并将同属于一个PF分量的占比相加得到/>即
其中,为每个PF分量中三个指标总和的占比;
d.设定固定阈值A,A为常数且满足0≤A≤1,将从大到小排序,将排序之后的/>从大到小进行选择,当选择的/>之和到达阈值A时,所选择的/>值对应的PF分量即为所选分量;
C.将上述MIS算法选择得到的PF分量进行重构,并将重构信号划分为训练集和测试集,再分别提取时域、频域特征和熵值特征,同时将不同故障类型打上标签作为该信号特征对应的系统状态,利用训练集训练轻量级梯度提升机(Light Gradient Boosting Machine,LightGBM)分类器,完成训练后,利用测试集评估分类效果,从而实现故障智能分类。
2.如权利要求1所述的基于多指标融合的故障特征选择及分类方法,其特征在于:所述的加入同幅值高斯白噪声的次数N=10,加入不同幅值高斯白噪声的次数N’=10,噪声的幅值hi依次取值0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50。
3.根据权利要求1所述的基于多指标融合的故障特征选择及分类方法,其特征在于:设置的阈值系数A=0.85。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210379191.1A CN114692695B (zh) | 2022-04-12 | 2022-04-12 | 基于多指标融合的故障特征选择及分类方法研究 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210379191.1A CN114692695B (zh) | 2022-04-12 | 2022-04-12 | 基于多指标融合的故障特征选择及分类方法研究 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114692695A CN114692695A (zh) | 2022-07-01 |
CN114692695B true CN114692695B (zh) | 2024-05-14 |
Family
ID=82143699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210379191.1A Active CN114692695B (zh) | 2022-04-12 | 2022-04-12 | 基于多指标融合的故障特征选择及分类方法研究 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114692695B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118296470B (zh) * | 2024-06-03 | 2024-08-13 | 中国人民解放军海军工程大学 | 一种多分类器融合的智能分类方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677472A (zh) * | 2017-08-25 | 2018-02-09 | 中国矿业大学 | 面向网络变量筛选与特征熵融合的轴承状态噪声诊断算法 |
CN109612732A (zh) * | 2019-01-24 | 2019-04-12 | 北华大学 | 一种基于解卷积和包络谱的滚动轴承故障诊断方法 |
CN113092112A (zh) * | 2021-03-30 | 2021-07-09 | 北京工业大学 | 一种基于eemd多特征融合的轴承复合故障诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933881A (zh) * | 2019-03-06 | 2019-06-25 | 武汉大学 | 一种基于优化深度信念网络的电力电子电路故障诊断方法 |
-
2022
- 2022-04-12 CN CN202210379191.1A patent/CN114692695B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677472A (zh) * | 2017-08-25 | 2018-02-09 | 中国矿业大学 | 面向网络变量筛选与特征熵融合的轴承状态噪声诊断算法 |
CN109612732A (zh) * | 2019-01-24 | 2019-04-12 | 北华大学 | 一种基于解卷积和包络谱的滚动轴承故障诊断方法 |
CN113092112A (zh) * | 2021-03-30 | 2021-07-09 | 北京工业大学 | 一种基于eemd多特征融合的轴承复合故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114692695A (zh) | 2022-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111037365B (zh) | 基于生成式对抗网络的刀具状态监测数据集增强方法 | |
CN110110768B (zh) | 基于并行特征学习和多分类器的滚动轴承故障诊断方法 | |
CN105275833B (zh) | 一种基于CEEMD-STFT时频信息熵和multi-SVM的离心泵故障诊断方法 | |
CN112257530B (zh) | 基于盲信号分离和支持向量机的滚动轴承故障诊断方法 | |
CN108760327B (zh) | 一种航空发动机转子故障的诊断方法 | |
CN109102005A (zh) | 基于浅层模型知识迁移的小样本深度学习方法 | |
CN112084974A (zh) | 一种基于元学习的多标签滚动轴承故障诊断方法 | |
CN109186973A (zh) | 一种无监督深度学习网络的机械故障诊断方法 | |
CN110333077B (zh) | 一种基于近似支持向量的滚动轴承故障诊断方法 | |
CN114358123B (zh) | 一种基于深度对抗迁移网络的广义开放集故障诊断方法 | |
CN110044623A (zh) | 经验模态分解剩余信号特征的滚动轴承故障智能识别方法 | |
CN112364706A (zh) | 一种基于类不平衡的小样本轴承故障诊断方法 | |
CN114692695B (zh) | 基于多指标融合的故障特征选择及分类方法研究 | |
CN109932179A (zh) | 一种基于ds自适应谱重构的滚动轴承故障检测方法 | |
CN109034076A (zh) | 一种机械故障信号的自动聚类方法及自动聚类系统 | |
CN108509701A (zh) | 一种基于振动信号的旋转机械故障直接智能诊断方法 | |
CN109858345A (zh) | 一种适用于胀管设备的智能故障诊断方法 | |
CN107451760A (zh) | 基于时窗滑移受限玻尔兹曼机的滚动轴承故障诊断方法 | |
CN105865784A (zh) | 基于局部均值分解和灰色关联的滚动轴承检测方法 | |
Deng et al. | Rolling bearing fault diagnosis based on Deep Boltzmann machines | |
CN111707472B (zh) | 基于fft和全连接层-svm的滚动轴承故障诊断方法 | |
Zhao et al. | A novel deep fuzzy clustering neural network model and its application in rolling bearing fault recognition | |
CN115587290A (zh) | 基于变分自编码生成对抗网络的航空发动机故障诊断方法 | |
Yu et al. | Rolling bearing fault feature extraction and diagnosis method based on MODWPT and DBN | |
CN112146880B (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 |