CN112561306A - 一种基于Hankel矩阵的滚动轴承健康状态评价方法 - Google Patents
一种基于Hankel矩阵的滚动轴承健康状态评价方法 Download PDFInfo
- Publication number
- CN112561306A CN112561306A CN202011457626.7A CN202011457626A CN112561306A CN 112561306 A CN112561306 A CN 112561306A CN 202011457626 A CN202011457626 A CN 202011457626A CN 112561306 A CN112561306 A CN 112561306A
- Authority
- CN
- China
- Prior art keywords
- matrix
- rolling bearing
- hankel matrix
- hankel
- time domain
- 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
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 122
- 238000005096 rolling process Methods 0.000 title claims abstract description 56
- 238000011156 evaluation Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 18
- 230000002159 abnormal effect Effects 0.000 claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 10
- 230000015556 catabolic process Effects 0.000 claims description 10
- 238000006731 degradation reaction Methods 0.000 claims description 10
- 230000003862 health status Effects 0.000 claims description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 6
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 2
- 238000010801 machine learning Methods 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- 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
- G06F2218/16—Classification; Matching by matching signal segments
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Strategic Management (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Optimization (AREA)
- Economics (AREA)
- Mathematical Analysis (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Computing Systems (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于Hankel矩阵的滚动轴承健康状态评价方法,对已有健康状态下振动时域信号进行分段并对获取的分段信号建立起对应的Hankel矩阵;其次建立每个分段信号的Hankel矩阵并计算平均Hankel矩阵的特征向量;利用特征向量对平均Hankel矩阵进行矩阵分解并将分解后的矩阵分解为对角矩阵和非对角矩阵;再计算非对角矩阵的1范数并设置10倍非对角矩阵的1范数为阈值;对于待检测的时域信号建立对应的Hankel矩阵,使用前述过程得到的特征向量对图连接矩阵进行矩阵分解,并计算异常值;若大于前述阈值则该段时域信号对应的滚动轴承为非健康状态。本发明克服了传统机器学习方法学习效率不好等问题,可以快速实现滚动轴承性能健康状态的快速辨识。
Description
技术领域
本发明涉及设备运行状态评价技术领域,特别涉及一种基于Hankel矩阵的滚动轴承健康状态评价方法。
背景技术
机械系统故障检测对于减少机械系统停机时间,预防灾难性故障具有重要意义。目前,针对故障特征提取提出了许多算法,但如何从包含大量干扰噪声的信号中监测机械系统的状态仍然是一个挑战。时域、频域、时频域的特征提取是实现滚动轴承健康状态评价的重要手段。为鲁棒评价滚动轴承健康状态,多种机器学习方法被提出并取得了较好的效果。然而,全寿命周期数据不完整导致机器学习方法的学习效率依然不高。
发明内容
本发明的目的在于,提供一种基于Hankel矩阵的滚动轴承健康状态评价方法。本发明可以快速实现滚动轴承性能健康状态的快速辨识。
本发明的技术方案:一种基于Hankel矩阵的滚动轴承健康状态评价方法,包括以下步骤:
S1:对已有健康状态下振动时域信号进行分段,分段信号长度为主轴旋转一圈时的时域信号;
S2:建立分段信号的Hankel矩阵;
S3:利用分段信号的Hankel矩阵构建分段信号的平均Hankel矩阵,并计算平均Hankel矩阵的特征向量;
S4:使用步骤S3中的特征向量对平均Hankel矩阵进行矩阵分解;
S5:提取矩阵非对角元素并计算其1范数;
S6:设置滚动轴承健康状态评价阈值;
S7:对于待检测的时域信号,建立待检测时域信号的Hankel矩阵,其中分段的信号长度与步骤S1中分段信号长度保持一致;
S8:使用步骤S3中的特征向量对待检测时域信号Hankel矩阵进行矩阵分解,并计算异常值;
S9:比较异常值与步骤S6中的阈值大小,若大于前述阈值则该段时域信号对应的滚动轴承为非健康状态。
上述的基于Hankel矩阵的滚动轴承健康状态评价方法,所述的步骤S1中,滚动轴承每转数据长度计算方法为:
其中,L为滚动轴承每转数据长度,fs为数据采样频率,r为滚动轴承旋转频率。
前述的基于Hankel矩阵的滚动轴承健康状态评价方法,所述的步骤S2中,Hankel矩阵的建立算法为:
其中,X为所建立的Hankel矩阵,x=[x1,x2…,xN]为原始采集信号序列,n为嵌入维数,N-n+1为各行子序列的长度;
平均Hankel矩阵的建立算法为:
前述的基于Hankel矩阵的滚动轴承健康状态评价方法,所述的步骤S4中,使用步骤S3中的特征向量对平均Hankel矩阵进行矩阵分解的算法为:
前述的基于Hankel矩阵的滚动轴承健康状态评价方法,提取矩阵非对角元素并计算其1范数,算法为;
其中,Mnon-diag为矩阵非对角元素序列,|| ||1为1范数计算符号,Mnon-diag为矩阵非对角元素序列第i个元素,n为矩阵非对角元素序列长度。
前述的基于Hankel矩阵的滚动轴承健康状态评价方法,所述的步骤S8中,计算的异常值为滚动轴承性能退化指标,算法为:
At=||non-diag(Yt)||1-||Mnon-diag||1
其中,At为滚动轴承性能退化指标,|| ||1为1范数计算符号,non-diag()为非对角元素提取符号,Yt为待测时域信号Hankel矩阵分解后的矩阵信息,||Mnon-diag||1为步骤S5中所得到非对角元素1范数。
前述的基于Hankel矩阵的滚动轴承健康状态评价方法,所述的步骤S9中,比较异常值与步骤S6中的阈值的算法为:
其中,Ft为待测时域信号的运行状态判断结果,t为周期数据序号,At为待测时域信号对应的滚动轴承性能退化指标,Fth为步骤S6中设置的阈值。
与现有技术相比,本发明具有以下有益效果:
本发明首先对已有健康状态下振动时域信号进行分段并对获取的分段信号建立起对应的Hankel矩阵;其次建立分段信号的Hankel矩阵并计算平均Hankel矩阵的特征向量;随后利用特征向量对平均Hankel矩阵进行矩阵分解;再提取矩阵非对角元素并计算其1范数,并设置10倍矩阵非对角元素的1范数为阈值;对于待检测的时域信号分段并建立对应的Hankel矩阵,使用前述过程得到的特征向量对待检测时域信号Hankel矩阵进行矩阵分解,并计算异常值;若大于前述阈值则该段时域信号对应的滚动轴承为非健康状态。本发明克服了传统机器学习方法学习效率不好等问题,可以快速实现滚动轴承性能健康状态的快速辨识,具有准确性高的优点。
附图说明
图1是本发明的流程图;
图2是本发明的实施例中振动时域信号分割示意图;
图3是本发明步骤S3生成的平均Hankel矩阵;
图4是本发明步骤S3计算的特征向量;
图5是本发明中实施例中滚动轴承性能健康的指标,其中其阈值设置为1000。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明,但并不作为对本发明限制的依据。
实施例:一种基于Hankel矩阵的滚动轴承健康状态评价方法,下面结合具体滚动轴承健康状态评价对本发明进一步说明,其流程图如图1所示,包括以下步骤,
S1:如附图2所示,对已有健康状态下振动时域信号进行分段,分段信号长度为主轴旋转一圈时的时域信号;
滚动轴承每转数据长度计算方法为:
其中,L为滚动轴承每转数据长度,fs为数据采样频率,r为滚动轴承旋转频率;
S2:建立每个分段信号的Hankel矩阵;
Hankel矩阵的建立算法为:
其中,X为所建立的Hankel矩阵,x=[x1,x2…,xN]为原始采集信号序列,n为嵌入维数,N-n+1为各行子序列的长度。
S3:利用每个分段信号的Hankel矩阵构建分段信号的平均Hankel矩阵,如附图3所示,随后,计算平均Hankel矩阵的特征向量,如图4所示;
平均Hankel矩阵的建立算法为:
S4:使用步骤S3中的特征向量对平均Hankel矩阵进行矩阵分解,算法为:
S5:提取矩阵非对角元素并计算其1范数,算法为;
其中,Mnon-diag为矩阵非对角元素序列,|| ||1为1范数计算符号,Mnon-diag为矩阵非对角元素序列第i个元素,n为矩阵非对角元素序列长度;
S6:设置滚动轴承健康状态评价阈值,10倍非对角矩阵的1范数为阈值,本实施例中的阈值设置为1000;
S7:对于待检测的时域信号,建立待检测时域信号的Hankel矩阵,其中分段的信号长度与步骤S1中分段信号长度保持一致;
S8:使用步骤S3中的特征向量对待检测时域信号Hankel矩阵进行矩阵分解,并计算异常值;如附图5所示(其中水平方向为时间序列,垂直方向为所计算的异常值),计算的异常值为滚动轴承性能退化指标,算法为:
At=||non-diag(Yt)||1-||Mnon-diag||1
其中,At为滚动轴承性能退化指标,|| ||1为1范数计算符号,non-diag()为非对角元素提取符号,Yt为当前周期数据Hankel矩阵分解后的矩阵信息,||Mnon-diag||1为步骤S5中所得到非对角元素1范数;
S9:比较异常值与步骤S6中的阈值大小,若大于前述阈值则该段时域信号对应的滚动轴承为非健康状,算法为:
其中,Ft为待测时域信号的运行状态判断结果,t为周期数据序号,At为待测时域信号对应的滚动轴承性能退化指标,Fth为步骤S7中设置的阈值。
从图5中看出,该段时域信号对应的滚动轴承性能退化指标在大部分情况下在阈值1000之上,说明该滚动轴承处于非健康状态,由此证明本发明可以快速实现滚动轴承性能健康状态的快速辨识。
Claims (7)
1.一种基于Hankel矩阵的滚动轴承健康状态评价方法,其特征在于:包括以下步骤:
S1:对已有健康状态下振动时域信号进行分段,分段信号长度为主轴旋转一圈时的时域信号;
S2:建立分段信号的Hankel矩阵;
S3:利用分段信号的Hankel矩阵构建分段信号的平均Hankel矩阵,并计算平均Hankel矩阵的特征向量;
S4:使用步骤S3中的特征向量对平均Hankel矩阵进行矩阵分解;
S5:提取矩阵非对角元素并计算其1范数;
S6:设置滚动轴承健康状态评价阈值;
S7:对于待检测的时域信号,建立待检测时域信号的Hankel矩阵,其中分段的信号长度与步骤S1中分段信号长度保持一致;
S8:使用步骤S3中的特征向量对待检测时域信号Hankel矩阵进行矩阵分解,并计算异常值;
S9:比较异常值与步骤S6中的阈值大小,若大于前述阈值则该段时域信号对应的滚动轴承为非健康状态。
6.根据权利要求1所述的基于Hankel矩阵的滚动轴承健康状态评价方法,其特征在于:所述的步骤S8中,计算的异常值为滚动轴承性能退化指标,算法为:
At=||non-diag(Yt)||1-||Mnon-diag||1
其中,At为滚动轴承性能退化指标,||||1为1范数计算符号,non-diag()为非对角元素提取符号,Yt为待测时域信号Hankel矩阵分解后的矩阵信息,||Mnon-diag||1为步骤S5中所得到非对角元素1范数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011457626.7A CN112561306B (zh) | 2020-12-11 | 2020-12-11 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011457626.7A CN112561306B (zh) | 2020-12-11 | 2020-12-11 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112561306A true CN112561306A (zh) | 2021-03-26 |
CN112561306B CN112561306B (zh) | 2023-12-08 |
Family
ID=75062159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011457626.7A Active CN112561306B (zh) | 2020-12-11 | 2020-12-11 | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112561306B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114486256A (zh) * | 2021-08-22 | 2022-05-13 | 北京燃气绿源达清洁燃料有限公司 | 一种cng压缩机滚动轴承故障特征提取方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100280772A1 (en) * | 2007-10-24 | 2010-11-04 | Abb Research Ltd. | Method for detection and automatic identification of damage to rolling bearings |
US20140039809A1 (en) * | 2012-07-31 | 2014-02-06 | Eurocopter | Method of detecting defects of a rolling bearing by vibration analysis |
CN104155108A (zh) * | 2014-07-21 | 2014-11-19 | 天津大学 | 一种基于振动时频分析的滚动轴承故障诊断方法 |
CN104198187A (zh) * | 2014-09-04 | 2014-12-10 | 昆明理工大学 | 一种机械振动故障特征时域盲提取方法 |
CN104807640A (zh) * | 2015-05-08 | 2015-07-29 | 重庆交通大学 | 基于振动敏感时频特征的空间滚动轴承寿命状态检测方法 |
CN105975749A (zh) * | 2016-04-28 | 2016-09-28 | 上海交通大学 | 一种轴承健康评估与预测的方法及系统 |
CN106198020A (zh) * | 2016-07-07 | 2016-12-07 | 华北电力大学(保定) | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 |
CN107121285A (zh) * | 2017-04-21 | 2017-09-01 | 南京理工大学 | 一种滚动轴承振动信号故障特征提取方法 |
CN108960328A (zh) * | 2018-07-04 | 2018-12-07 | 温州大学 | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 |
WO2019061006A1 (en) * | 2017-09-26 | 2019-04-04 | Schaeffler Technologies AG & Co. KG | METHOD AND DEVICE FOR DIAGNOSING BEARING FAULT, READABLE STORAGE MEDIUM, AND ELECTRONIC DEVICE |
CN109582917A (zh) * | 2018-10-24 | 2019-04-05 | 加康康健有限公司 | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 |
-
2020
- 2020-12-11 CN CN202011457626.7A patent/CN112561306B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100280772A1 (en) * | 2007-10-24 | 2010-11-04 | Abb Research Ltd. | Method for detection and automatic identification of damage to rolling bearings |
US20140039809A1 (en) * | 2012-07-31 | 2014-02-06 | Eurocopter | Method of detecting defects of a rolling bearing by vibration analysis |
CN104155108A (zh) * | 2014-07-21 | 2014-11-19 | 天津大学 | 一种基于振动时频分析的滚动轴承故障诊断方法 |
CN104198187A (zh) * | 2014-09-04 | 2014-12-10 | 昆明理工大学 | 一种机械振动故障特征时域盲提取方法 |
CN104807640A (zh) * | 2015-05-08 | 2015-07-29 | 重庆交通大学 | 基于振动敏感时频特征的空间滚动轴承寿命状态检测方法 |
CN105975749A (zh) * | 2016-04-28 | 2016-09-28 | 上海交通大学 | 一种轴承健康评估与预测的方法及系统 |
CN106198020A (zh) * | 2016-07-07 | 2016-12-07 | 华北电力大学(保定) | 基于子空间和模糊c均值聚类的风电机组轴承故障诊断法 |
CN107121285A (zh) * | 2017-04-21 | 2017-09-01 | 南京理工大学 | 一种滚动轴承振动信号故障特征提取方法 |
WO2019061006A1 (en) * | 2017-09-26 | 2019-04-04 | Schaeffler Technologies AG & Co. KG | METHOD AND DEVICE FOR DIAGNOSING BEARING FAULT, READABLE STORAGE MEDIUM, AND ELECTRONIC DEVICE |
CN108960328A (zh) * | 2018-07-04 | 2018-12-07 | 温州大学 | 一种基于汉克尔矩阵的诊断汽车轮毂轴承故障的方法 |
CN109582917A (zh) * | 2018-10-24 | 2019-04-05 | 加康康健有限公司 | 一种基于ssa的信号去噪方法、装置、终端设备及存储介质 |
Non-Patent Citations (3)
Title |
---|
SHUAI ZHANG 等: "Multi-Channel Hankel Matrix Completion through Nonconvex Optimization", 《IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING》, vol. 12, no. 4, pages 617 - 632 * |
石远程 等: "基于 MEEMD 的滚动轴承故障诊断方法", 《煤矿机械》, vol. 35, no. 6, pages 266 - 268 * |
苗乃树 等: "基于卷积稀疏滤波的轴承微弱故障检测方法", 《农业装备与车辆工程》, vol. 58, no. 9, pages 131 - 134 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114486256A (zh) * | 2021-08-22 | 2022-05-13 | 北京燃气绿源达清洁燃料有限公司 | 一种cng压缩机滚动轴承故障特征提取方法 |
CN114486256B (zh) * | 2021-08-22 | 2023-10-31 | 北京燃气绿源达清洁燃料有限公司 | 一种cng压缩机滚动轴承故障特征提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112561306B (zh) | 2023-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111562108A (zh) | 一种基于cnn和fcmc的滚动轴承智能故障诊断方法 | |
CN112766342A (zh) | 一种电气设备的异常检测方法 | |
CN109406949B (zh) | 基于支持向量机的配电网早期故障检测方法及装置 | |
CN111881594B (zh) | 一种核动力设备的非平稳信号状态监测方法及系统 | |
CN115481657A (zh) | 一种基于电信号的风力发电机组通信滑环故障诊断方法 | |
CN113155462A (zh) | 基于辛几何模式分解与图结构增强动态时间规整的轴承故障诊断方法 | |
CN112561306A (zh) | 一种基于Hankel矩阵的滚动轴承健康状态评价方法 | |
CN114492534B (zh) | 一种跨尺寸电机轴承故障诊断模型的构建方法及应用 | |
CN116168720A (zh) | 一种电机声音异常故障诊断方法、系统、及可存储介质 | |
CN111929489B (zh) | 故障电弧电流的检测方法及系统 | |
CN115683319A (zh) | 一种电力变压器绕组状态评估方法 | |
CN114861774A (zh) | 一种电网中虚假数据识别方法及系统 | |
CN112464209B (zh) | 一种电力终端指纹认证方法及装置 | |
CN107147143B (zh) | 一种风机连锁脱网故障早期预警模型建立方法 | |
CN117743120A (zh) | 基于大数据技术的网卡异常检测方法 | |
CN116337447B (zh) | 一种非平稳工况下轨道车辆轮对轴承故障诊断方法及设备 | |
CN114112390A (zh) | 一种非线性复杂系统早期故障诊断方法 | |
CN115561575A (zh) | 基于多维矩阵轮廓的海上风电场电气异常状态辨别方法 | |
CN112014821B (zh) | 一种基于雷达宽带特征的未知车辆目标识别方法 | |
CN110376436B (zh) | 一种多尺度噪声功率谱线谱检测方法 | |
CN107884188B (zh) | 一种基于模式转换的变转速/变负载滚动轴承故障诊断方法 | |
CN112595515A (zh) | 一种动力轴系轴承故障检测方法及系统 | |
CN112697268A (zh) | 一种基于t-sne的电机异常检测集成算法 | |
CN118514566B (zh) | 一种电瓶车违规充电的检测方法 | |
CN116643165B (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 |