CN113052018A - 一种钣金拉伸过程声发射信号识别方法 - Google Patents

一种钣金拉伸过程声发射信号识别方法 Download PDF

Info

Publication number
CN113052018A
CN113052018A CN202110258631.3A CN202110258631A CN113052018A CN 113052018 A CN113052018 A CN 113052018A CN 202110258631 A CN202110258631 A CN 202110258631A CN 113052018 A CN113052018 A CN 113052018A
Authority
CN
China
Prior art keywords
acoustic emission
signals
stretching process
tensile
sheet metal
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.)
Granted
Application number
CN202110258631.3A
Other languages
English (en)
Other versions
CN113052018B (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.)
Beijing University of Technology
Original Assignee
Beijing University of 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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202110258631.3A priority Critical patent/CN113052018B/zh
Publication of CN113052018A publication Critical patent/CN113052018A/zh
Application granted granted Critical
Publication of CN113052018B publication Critical patent/CN113052018B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种钣金拉伸过程声发射信号识别方法,利用声发射检测系统采集钣金拉伸过程中的声发射信号,对声发射信号的傅里叶变换或短时傅里叶变换进行非负矩阵分解,提取其在低维子空间映射的特征系数,用于构造训练字典和测试样本,并利用稀疏表示方法实现对拉伸过程不同阶段的表征与识别。该方法通过对声发射信号的非负矩阵分解和稀疏表示分析,可以实现声发射对应拉伸过程的自动识别,且识别率高。

Description

一种钣金拉伸过程声发射信号识别方法
技术领域
本发明涉及一种基于非负矩阵分解和稀疏表示的板材拉伸过程声发射信号 识别方法,该方法适用于金属板材拉伸过程的声发射信号识别,属于无损检测领 域。
背景技术
钣金由于具有重量轻、强度高、成本低、大规模量产性能好等特点,在电子 电器、通信、汽车工业、医疗器械等领域广泛应用。钣金成形涉及冲压、弯曲、 拉伸等多种加工手段。同时,在钣金成形过程中,材料经历复杂的应力-应变状 态。为此,这些加工手段的工艺参数对钣金成形质量及稳定性有很大的影响。因 此,研究钣金成形过程中材料的应力-应变状态对于提高钣金成形质量及稳定性 具有重要意义。
声发射是指材料因局部能量的快速释放而发出瞬态弹性波的现象。作为一种 动态监测技术,声发射技术已经被应用于多种材料加工过程的监测[1-3]。例如, 张颖[4]对四种不同晶粒度的钢试件单轴拉伸过程进行了声发射检测,利用声发 射信号的幅值、能量和撞击数等特征参数描述了晶粒度对材料拉伸过程中声发射 特征的影响。丁利伟[5]利用声发射技术对复合材料拉伸过程进行了监测,通过 分析拉伸过程中声发射信号的能量、振铃技术、幅值等特征参数,将拉伸过程分 为了三个阶段。诸如以上文献所述,利用声发射技术进行拉伸过程监测的研究工 作很多,但目前已有文献仅停留在研究不同拉伸过程中声发射信号特征参数的变 化规律,鲜见有利用声发射技术进行拉伸过程状态表征的研究。
本发明涉及一种钣金拉伸过程声发射信号识别方法。通过对钣金拉伸过程声 发射信号进行非负矩阵分解,提取其在低维子空间映射的特征系数,用于构造训 练字典和测试样本,并利用稀疏表示方法实现对拉伸过程不同阶段的表征与识 别。
发明内容
本发明的目的在于提供一种基于非负矩阵分解和稀疏表示的拉伸过程声发 射信号识别方法,通过该方法可以自动实现声发射信号对应的拉伸过程的识别。 本方法通过对钣金拉伸过程声发射信号进行非负矩阵分解,提取其在低维子空间 映射的特征系数,用于构造训练字典和测试样本,并利用稀疏表示方法实现对拉 伸过程不同阶段的表征与识别。
本发明提出的一种基于非负矩阵分解和稀疏表示的钣金拉伸过程识别方法, 其基本原理如下:
在拉伸试件两端放置声发射传感器,并置于拉伸试验机上,用于接收拉伸过 程的声发射信号。根据拉伸试验机获取的应力应变曲线,判定拉伸过程不同力学 阶段的时间范围,进而提取不同阶段的声发射时域信号。对得到的声发射信号进 行傅里叶变换(FastFourier Transform,FFT)或短时傅里叶变换(Short Time Fourier Transformation,STFT),在此基础上对其进行非负矩阵分解(Nonnegative Matrix Factorization,NMF),构造出非负训练字典V。本发明以声发射信号的短 时傅里叶变换后的数据为例,用以说明基于非负矩阵分解和稀疏表示的钣金拉伸 过程声发射信号识别方法。
训练字典V∈Rm×n,其中m是数据长度,n是样本数目。V可以被两个低秩 非负矩阵W∈Rm×r和H∈Rr×n近似表示,V=WH,其中W称作基矩阵,H称作系 数矩阵,其中秩r的取值满足(m+n)*r<n<m*n。
根据Lee和Seung的乘法更新规则,同时以欧氏距离为目标函数求解非负矩 阵,得到低秩矩阵W和H。
Figure BDA0002968619120000021
基于欧式距离的目标函数的求解问题可以通过最优化算法进行求解,而梯度 下降法是用来求解极小值问题的有效方法,沿着负梯度方向,反复迭代求解,最 后就能解出局部最小值。应用梯度下降法迭代求解式(1):
Figure BDA0002968619120000022
在计算过程中,随机选取非负正态分布的初始化矩阵W与H,在定义适当 的迭代次数之后,迭代计算即可得到最优的W和H。迭代求解出的基矩阵W是V的基空间,根据穆尔-彭罗斯广义逆矩阵,求解其子空间矩阵W+[9]。因此,可 以得到矩阵V在子空间W+上的投影Y,Y=W+V,Y即是新的训练字典。
稀疏表示的核心思想是把同类训练样本视为一组基,将测试样本在不同类 基上线性表示,然后通过求解一个l1-范数最优化问题,计算出测试样本在不同 类别训练子字典上的最稀疏表示系数,最后根据测试样本稀疏表示后的重构误差 大小判定测试样本所属类别。信号稀疏表示的两个主要内容是字典的构造和稀疏 系数的求解,字典的构造是由声发射信号在非负矩阵子空间映射完成的,稀疏系 数的求解由l1范数求解。
假设有c类样本信号,其中每类有n个训练样本,则第i类训练子字典Yi
Yi=[yi1,yi2,···,yin]∈Rm×n (3)
其中,m是样本信号的数据长度。
假设y为第i类测试样本,y∈Rm,测试样本y可表示为所有训练样本的线 性组合:
y=Ya,a=(0,…0,ai,n1,…ai,ni,0,…0) (4)
其中,Y=[Y1,…Yi,…Yc]∈Rm×nc,Y是由n*c个训练样本构造的训练字典。a∈Rn,称为稀疏系数。
在对测试样本y进行稀疏表示分类时,需要对稀疏系数a进行求解,该求解 过程应该满足:
Figure BDA0002968619120000031
满足||y-Ya||2≤ε (5)
其中,ε是测试样本重构误差的容忍极限,a'是稀疏系数a的最稀疏解。
而当稀疏系数a足够稀疏时,l0范数等效为l1范数求解,并用拉格朗日定理 做正则化处理,故式(5)可以写为:
Figure BDA0002968619120000032
其中λ是正则化参数,其取值范围为(0,1)。
则测试样本y在不同类别子字典下的重构信号y′i
y′i=Yia′i (7)
a′i是测试样本y在第i类子字典上对应的稀疏系数,最后计算测试样本y与 重构y′i的误差值errori
errori=||y-y′i||2 (8)
根据误差值的最小原则判断当前测试样本y所属的类别:
Figure BDA0002968619120000033
本发明的技术方案如下:
本发明所采用的装置参见图1,包括拉伸试验机1、两块铝板试件2、声发 射传感器3、前置放大器4、全信息声发射分析仪5和计算机(配套分析软件)6。 设定拉伸试验机1的运行参数并拉伸铝板试件2,并由两个声发射传感器3同时 采集声发射信号,然后通过前置放大器4进行放大,并由全信息声发射分析仪5 记录双通道信号,最后在计算机6上进行信号分析、处理。
本发明提出的一种钣金拉伸过程声发射信号识别方法是通过以下步骤实现 的:
步骤一:按照图1所示的检测装置系统搭建拉伸声发射检测系统,系统包括 拉伸试验机1、铝板试件2、声发射传感器3、前置放大器4、全信息声发射分析 仪5和计算机6(配套分析软件)。根据拉伸试验机获得的应力应变曲线,判定 拉伸过程各个力学阶段的时间范围,随后在拉伸过程各阶段分别截取n组(共 5n组)信号。基于MATLAB软件对所有声发射信号进行截断,获得数据长度为 m的信号。在此基础上,对声发射信号进行短时傅里叶变换,得到5n个时频图, 其中每个时频图尺寸为w*l。
步骤二:分别随机抽取拉伸过程各阶段n组信号中的n1组(共5n1组)声发 射信号的时频图,并转换为列向量,每一列代表一个信号,组成各阶段训练子字 典vi(i=1,2,3,4,5),
Figure BDA0002968619120000041
剩余各阶段n2组(共5n2组)信号组成测试样本 集v0
Figure BDA0002968619120000042
步骤三:设定非负矩阵分解维度为N,迭代次数M,并使用随机非负正态 分布函数初始化W和H,通过式(2)对所有的训练子字典vi(i=1,2,3,4,5)进行 非负矩阵分解,得到基空间矩阵W,W∈Rwl×N
步骤四:基于穆尔-彭罗斯广义逆矩阵,求解基矩阵W的子空间W+, W+∈RN×wl。将所有训练子字典vi(i=1,2,3,4,5)分别向子空间W+进行投影,得 到每个训练子字典的特征系数矩阵yi(i=1,2,3,4,5),
Figure BDA0002968619120000043
所有训练子字典 的特征系数矩阵共同形成降维后的训练字典Vtrain
Figure BDA0002968619120000044
步骤五:同理,测试样本集v0向子空间W+进行投影,由所有测试样本的特 征系数共同形成降维后的测试样本集y0
Figure BDA0002968619120000045
步骤六:通过公式(6)求解测试样本集y0在训练字典上Vtrain的稀疏系数, 并通过公式(7)求取所有测试样本在不同阶段训练子字典上的重构信号。
步骤七:比较所有测试样本在不同阶段上重构信号的误差,并根据误差值的 最小原则判断所有测试样本的类别。最后,重复步骤二、三、四、五、六和七的 实验,得到十次平均后的基于非负矩阵分解和稀疏表示的声发射信号拉伸过程的 识别率。
本发明具有以下优点:(1)可以自动实现声发射信号对应拉伸过程的识别; (2)具有更高的识别准确率;(3)该算法运算速度快。
附图说明
图1拉伸声发射检测系统图。
图2铝板拉伸过程的应力应变曲线。
图3基于非负矩阵分解和稀疏表示的声发射信号识别率图。
图4本方法的实施流程图。
具体实施方式
下面结合具体实验对本发明作进一步说明:
本实验实施过程包括以下步骤:
步骤一:按照图1所示的检测装置系统搭建拉伸声发射检测系统,系统包括 拉伸试验机1、铝板试件2、声发射传感器3、前置放大器4、全信息声发射分析 仪5和计算机(配套分析软件)6。根据拉伸试验机获得的应力应变曲线如图2 所示,判定拉伸过程各个力学阶段的时间范围如表1所示,随后在拉伸过程各阶 段分别截取110组(共550组)信号。基于MATLAB软件对所有声发射信号进 行截断,获得频率3MHz、数据长度为2048的信号。最后,对每个声发射信号 进行短时傅里叶变换,得到550个时频图,其中每个时频图尺寸为129*18。
步骤二:分别随机抽取拉伸过程各阶段110组信号中的100组(共500组) 声发射信号的时频图,并转换为列向量,每一列代表一个信号,组成各阶段训练 子字典vi(i=1,2,3,4,5),vi∈R2322×100,剩余各阶段10组(共50组)信号组成测 试样本集v0,v0∈R2322×50
步骤三:设定非负矩阵分解维度为240,迭代次数500,并使用随机非负正 态分布函数初始化W和H。通过式(2)对所有的训练子字典vi(i=1,2,3,4,5)进 行非负矩阵分解,得到基空间矩阵W,W∈R2322×240
步骤四:基于穆尔-彭罗斯广义逆矩阵,求解基矩阵W的子空间W+, W+∈R240×2322。将所有训练子字典vi(i=1,2,3,4,5)分别向子空间W+进行投影, 得到每个训练子字典的特征系数矩阵yi(i=1,2,3,4,5),yi∈R240×100,所有训练子 字典的特征系数矩阵共同形成降维后的训练字典Vtrain,Vtrain∈R240×500
步骤五:同理,测试样本集v0向子空间W+进行投影,由所有测试样本的特 征系数共同形成降维后的测试样本集y0,y0∈R240×50
步骤六:通过公式(6)求解测试样本集y0在训练字典上Vtrain的稀疏系数, 并通过公式(7)求取所有测试样本在不同阶段训练子字典上的重构信号。
步骤七:比较所有测试样本在不同阶段上重构信号的误差,并根据误差值的 最小原则判断所有测试样本的类别。最后,重复步骤二、三、四、五、六和七的 实验,得到十次平均后的非负矩阵分解的稀疏表示分类算法识别率如图3所示。
与现有分类方法相比,本发明具有以下优点:(1)具有更高的识别准确率, 适应更复杂种类的样本;(2)适用于高维数据,非负分解后的训练字典具有低维 和稀疏性的特点;(3)具有更快的运算速度,较少程序运行时间。
本方法的实施流程图如图4所示。
表1拉伸过程各阶段时间范围
Figure RE-GDA0003063162650000061
以上是本发明的一个典型应用,本发明的应用不限于此。
参考文献
[1]Behrens B A,Bouguecha A,Buse C.Potentials of in situ monitoring ofaluminum alloy forging by acoustic emission[J].Archives of Civil andMechanical Engineering,2016,16(4):724-733.
[2]Hao S,Ramalingam S,Klamecki B E.Acoustic emission monitoring ofsheet metal forming:characterization of the transducer,the work material andthe process[J].Journal of Materials Processing Technology,2000,101(1–3):124-136.
[3]张进,柴孟瑜,项靖海,段权,李丽婵.316LN不锈钢断裂过程的声 发射特性[J].材料研究学报,2018,32(6):415-422.
[4]张颖,吴昊,高唅,王兵,张维.20#钢不同晶粒度试件拉伸损伤试验的声 发射特性研究,应用声学,2017,36(3):228-232.
[5]丁利伟,沈玉娣.复合材料拉伸过程的声发射特性研究.无损检 测,2009,31(10):781-784.

Claims (3)

1.一种钣金拉伸过程声发射信号识别方法,其特征在于:该方法通过以下步骤实现:
步骤一:搭建拉伸声发射检测系统,根据拉伸试验机获得的应力应变曲线,判定拉伸过程各个力学阶段的时间范围,随后在拉伸过程各阶段分别截取n组信号,共5n组;基于MATLAB软件对所有声发射信号进行截断,获得数据长度为m的信号;最后,对所有声发射信号进行短时傅里叶变换,得到5n个时频图,其中每个时频图尺寸为w*l;
步骤二:分别随机抽取拉伸过程各阶段n组信号中的n1组声发射信号的时频图,共5n1组,并转换为列向量,每一列代表一个信号,组成各阶段训练子字典v,
Figure FDA0002968619110000011
剩余各阶段n2组信号组成测试样本集v0,共5n2组,
Figure FDA0002968619110000012
步骤三:设定非负矩阵分解维度为N,迭代次数M,并使用随机非负正态分布函数初始化W和H,对所有的训练子字典vi进行非负矩阵分解,得到基空间矩阵W,W∈Rwl×N,i=1,2,3,4,5;
步骤四:基于穆尔-彭罗斯广义逆矩阵,求解基矩阵W的子空间W+,W+∈RN×wl;将所有训练子字典vi分别向子空间W+进行投影,得到每个训练子字典的特征系数矩阵yi
Figure FDA0002968619110000013
所有训练子字典的特征系数矩阵共同形成降维后的训练字典Vtrain
Figure FDA0002968619110000014
步骤五:测试样本集v0向子空间W+进行投影,由所有测试样本的特征系数共同形成降维后的测试样本集y0
Figure FDA0002968619110000015
步骤六:求解测试样本集y0在训练字典上Vtrain的稀疏系数,并求取所有测试样本在不同阶段训练子字典上的重构信号;
步骤七:比较所有测试样本在不同阶段上重构信号的误差,并根据误差值的最小原则判断所有测试样本的类别;最后,重复步骤二、三、四、五、六和七的实验,得到十次平均后的基于非负矩阵分解和稀疏表示的声发射信号拉伸过程的识别率。
2.根据权利要求1所述的一种钣金拉伸过程声发射信号识别方法,其特征在于:拉伸声发射检测系统包括拉伸试验机(1)、两块铝板试件(2)、声发射传感器(3)、前置放大器(4)、全信息声发射分析仪(5)和计算机(6);设定拉伸试验机(1)的运行参数并拉伸铝板试件(2),并由两个声发射传感器(3)同时采集声发射信号,然后通过前置放大器(4)进行放大,并由全信息声发射分析仪(5)记录双通道信号,最后在计算机(6)上进行信号分析、处理。
3.根据权利要求1所述的一种钣金拉伸过程声发射信号识别方法,其特征在于:在钣金拉伸试件两端放置声发射传感器,并置于拉伸试验机上,用于接收拉伸过程的声发射信号;根据拉伸试验机获取的应力应变曲线,判定拉伸过程不同力学阶段的时间范围,进而提取不同阶段的声发射时域信号;对得到的声发射信号进行傅里叶变换FFT或短时傅里叶变换STFT,进行非负矩阵分解NMF,构造出非负训练字典V。
CN202110258631.3A 2021-03-09 2021-03-09 一种钣金拉伸过程声发射信号识别方法 Active CN113052018B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110258631.3A CN113052018B (zh) 2021-03-09 2021-03-09 一种钣金拉伸过程声发射信号识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110258631.3A CN113052018B (zh) 2021-03-09 2021-03-09 一种钣金拉伸过程声发射信号识别方法

Publications (2)

Publication Number Publication Date
CN113052018A true CN113052018A (zh) 2021-06-29
CN113052018B CN113052018B (zh) 2024-03-15

Family

ID=76511841

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110258631.3A Active CN113052018B (zh) 2021-03-09 2021-03-09 一种钣金拉伸过程声发射信号识别方法

Country Status (1)

Country Link
CN (1) CN113052018B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2112603A2 (en) * 2008-04-25 2009-10-28 Xerox Corporation Clustering using non-negative matrix factorization on sparse graphs
CN106326871A (zh) * 2016-08-26 2017-01-11 中国计量大学 一种基于字典分解和稀疏表示的鲁棒人脸识别方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2112603A2 (en) * 2008-04-25 2009-10-28 Xerox Corporation Clustering using non-negative matrix factorization on sparse graphs
CN106326871A (zh) * 2016-08-26 2017-01-11 中国计量大学 一种基于字典分解和稀疏表示的鲁棒人脸识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
PENGJING ZHAO 等: "Correlation between acoustic emission detection and microstructural characterization for damage evolution", ENGINEERING FRACTURE MECHANICS, 28 February 2020 (2020-02-28), pages 1 - 15 *
焦敬品 等: "基于超声局部谐振的板结构缺陷定量检测方法", 仪器仪表学报, vol. 40, no. 12, 31 December 2019 (2019-12-31), pages 1 - 8 *

Also Published As

Publication number Publication date
CN113052018B (zh) 2024-03-15

Similar Documents

Publication Publication Date Title
Pan et al. Compressed sensing based on dictionary learning for reconstructing blade tip timing signals
CN107247687B (zh) 一种基于特征正交分解的非平稳随机过程快速模拟方法
CN111168569B (zh) 一种磨削材料去除量预测方法、装置、设备及存储介质
CN104200065A (zh) 基于随机共振与核主元分析相结合的机械振动信号特征提取方法
Bai et al. Phased array ultrasonic signal compressive detection in low-pressure turbine disc
CN111595541A (zh) 一种海量振动传递率数据卷积神经网络处理的多维结构损伤识别方法
Zhai et al. Pattern recognition approach to identify loose particle material based on modified MFCC and HMMs
CN110487499B (zh) 一种基于奇异值分解和支持向量机的螺栓预紧力辨识方法
CN114118174A (zh) 一种滚动轴承故障诊断方法及系统
CN113052018B (zh) 一种钣金拉伸过程声发射信号识别方法
CN116881712A (zh) 一种混凝土大坝活动裂纹电磁脉冲信号识别方法
Yang et al. Independent component analysis–based feature extraction technique for defect classification applied for pulsed eddy current NDE
CN111275003A (zh) 基于类最优高斯核多分类支持向量机的微震信号识别方法
CN111191502A (zh) 基于钻柱振动信号的粘滑和跳钻异常工况识别方法
CN115901954A (zh) 超声导波管道缺陷无损检测方法
CN115545202A (zh) 一种基于量子门获取待测系统本征值的方法及装置
CN112183297B (zh) 一种超声相控阵信号稀疏特征提取方法
Chen et al. Deep compressed sensing-based accurate blade vibration reconstruction using sub-sampled tip-timing measurements
Chen et al. Deep compressed sensing-based blade vibration reconstruction from sub-sampled tip-timing signals
Lu et al. Denoising method for capillary electrophoresis signal via learned tight frame
CN107977680A (zh) 一种信号稀疏重构的方法
Zheng et al. Research on composite impact localization by MUSIC based guided wave monitoring
Jaiswal et al. Wavelet transform of modal data of beam in damage detection exercise
CN109406084B (zh) 一种对复合材料产品实现全频域快速冲击监测的方法
Ravanfar et al. Damage detection optimization using wavelet multiresolution analysis and genetic algorithm

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