CN111660141A - 一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 - Google Patents
一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 Download PDFInfo
- Publication number
- CN111660141A CN111660141A CN202010405649.7A CN202010405649A CN111660141A CN 111660141 A CN111660141 A CN 111660141A CN 202010405649 A CN202010405649 A CN 202010405649A CN 111660141 A CN111660141 A CN 111660141A
- Authority
- CN
- China
- Prior art keywords
- driving current
- milling cutter
- main shaft
- wear state
- signal
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0952—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining
- B23Q17/0957—Detection of tool breakage
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Automatic Control Of Machine Tools (AREA)
Abstract
本发明公开了一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,该方法基于刀具磨损和振动异常会导致主轴驱动电流信号出现不规则杂波成分的试验结果,利用傅里叶级数拟合将主轴驱动电流波形分解为反映电流有效值准静态变化的谐波成分和反映铣刀刃口和后刀面磨损状态以及振动异常的电流杂波信号,然后将电流杂波信号输入到卷积神经网络中进行铣刀状态特征提取和分类。通过上述方式,本发明能够实现复杂工况变切削条件下的铣刀状态的在线准确辨识,为复杂工况下预测铣刀剩余寿命和科学合理制定铣刀更换规则打下基础。
Description
技术领域
本发明涉及一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,属于数控机床状态检测技术领域。
背景技术
在金属铣削加工过程中,铣刀的后刀面及刀刃受切削应力与工件摩擦会逐渐累积磨损。同时,刀具与工件之间会产生相对振动,刀具磨损和振动异常不仅直接影响被加工表面质量和加工精度,使工件表面粗糙度增大,并且随着刀具磨损以及振动的加剧导致切削力不断的增大,严重的会造成刀具破损,致使工件报废、机床损坏,甚至影响工厂生产正常进行,给企业带来经济损失。目前国内外的刀具状态监测方法主要分为直接法和间接法两类。直接法采用测量刀刃外观、几何尺寸等,虽然对刀具状态识别精度高,但是大部分直接监测法需要停机检测,只适用于非加工过程,不能满足在线监测的要求,在实际生产加工过程中实用性不强。而间接法主要通过采集与刀具磨损相关的物理信号,例如,力信号、振动信号、电流信号和声发射信号等,然后运用信号处理技术提取与刀具磨损状态相关的特征,达到准确识别刀具磨损状态的目的。间接法主要优势在于可以实时监测刀具磨损状态,提高了生产效率。但是间接监测方法对切削力、振动信号的采集需要将传感器安装在机床加工区域内,干扰机床的正常加工,同时并不适用于复杂工况变切削条件下的铣刀磨损状态在线辨识。
发明内容
本发明的目的是为了解决现有技术不足的问题,设计了一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法。
为解决上述技术问题,本发明所采用技术的方案为一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法。该方法基于刀具磨损以及振动异常会导致主轴驱动电流信号出现不规则杂波成分的试验结果,利用傅里叶级数拟合提取与加工工艺参数相关的反应铣削力准静态变化的主轴驱动电流基波成分,然后剔除这部分电流成分对不同工况下刀具磨损状态在线辨识的影响,将剩余的与工况无关的反映铣刀磨损和振动状态的主轴驱动电流杂波信号输入到一维卷积神经网络中进行刀具磨损状态特征提取和分类。
铣削加工过程中的铣刀磨损状态分为正常磨损、急剧磨损、磨损振动异常和锋利振动异常四种状态。
数据采集是通过嵌入式数字采集系统与霍尔电流互感器以及三相加速度传感器建立通信获取金属铣削加工过程中的主轴驱动电流信号和主轴振动信号。
在保证拟合准确率和拟合效率的前提下,将获取到的主轴驱动电流信号经过四阶傅里叶级数拟合得到主轴驱动电流基波成分。
主轴驱动电流杂波信号是通过在主轴驱动电流原始信号中剔除电流基波成分后得到。
数据归一化处理将提取到的主轴驱动电流杂波信号限制在0到1范围内,保证程序运行时收敛加快。
深度学习神经网络采用的是一维卷积神经网络,网络结构由一层输入层,两层卷积层,两层池化层,一层扁平层,一层全连接层以及输出层组成。其中输入层是128*1的输入矩阵,两层卷积层分别使用的是24个9×1的卷积核和48个9×1的卷积核进行图像过滤。池化层使用的是6×1的池化核来对上一层数据图缩放,池化方式为最大池化。激活函数采用Relu函数,使用Softmax分类器对提取到的特征分类。最后根据实际分类与样本分类构建损失函数,利用逐层反向传播算法来调整网络的参数。
通过对比分析同磨损阶段不同工况下的主轴驱动电流基波成分,当加工工艺参数发生改变时,主轴驱动电流基波成分会随之改变。
对不同工况下的主轴驱动电流杂波成分进行对比分析得到,当刀具处于相同磨损阶段时,主轴驱动电流杂波信号在一定范围内波动;当刀具处于不同磨损阶段时,主轴驱动电流杂波成分会随着刀具磨损以及振动的加剧而变大。
本发明的显著效益是基于铣刀磨损和振动异常会导致主轴驱动电流信号出现不规则杂波成分的实验结果,通过四阶傅里叶级数拟合提取铣削加工过程中的电流杂波信号并输入到卷积神经网络中进行刀具磨损状态特征提取和分类,实现了数控机床复杂工况变切削条件下刀具磨损状态的准确在线辨识。一方面可以降低工件废品率和机床故障率,减少工厂成本,另一方面也为复杂工况下科学合理制定换刀规则以及铣刀剩余寿命的在线预测打下基础。
附图说明
图1为本发明的刀具磨损状态识别流程图。
图2为本发明的电流原始信号拟合图。
图3为本发明的不同状态下主轴驱动电流杂波信号。
图4为本发明采用的卷积神经网络结构图。
图5为本发明电流有效值曲线与刀具状态分类图。
具体实施方式
下面结合附图,对本发明进一步详细说明,以下实例仅用于解释本发明,并不构成对本发明保护范围的限定。
如图1所示为本发明的刀具磨损状态识别流程图。本发明提供的是一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,包括以下步骤:
a.采集金属铣削加工时主轴驱动电流的原始信号和主轴的振动信号;
b.对采集到的主轴驱动电流的原始信号通过傅里叶级数拟合,得到一种与工况无关的主轴驱动电流杂波信号;
c.将获取到的杂波信号经过数据归一化处理,使所有的信号幅值都处于0到1之间;
d.将处理后的数据输入到一维卷积神经网络中进行训练,根据实际分类与样本分类构建损失函数,利用逐层反向传播算法来调整网络的参数;
e.将测试数据输入到已训练完成的深度学习神经网络中进行铣刀磨损状态识别。
金属铣削加工过程中的数据采集首先在数控机床上安装三相加速度传感器和霍尔电流传感器,然后通过嵌入式数字采集系统与传感器之间建立通信采集铣削加工过程中的振动和电流的连续时域信号。
金属铣削加工过程中采集到的主轴驱动电流信号由电流基波和电流杂波两部分组成,其中电流基波信号反映的是铣削加工过程中的准静态铣削力分量,电流杂波信号反映的是铣削加工过程中的动态铣削力分量。准静态铣削力分量反映了由切削工艺参数、加工用量决定的随铣刀旋转和进给运动将金属材料从工件表面剥离下来所需的切削力;而动态铣削力分量一方面包含了铣削过程中刀具与工件之间相对振动导致的动态切削力,另一方面也包含了因刃口钝化和后刀面不规则磨损引起的铣削力的动态波动量。随着铣刀磨损和振动的加剧,动态铣削力变动量会逐渐变大,主轴转速出现波动,进一步引起电流信号出现波动。因此,随着铣刀磨损和振动加剧,动态铣削力增大会导致电流信号变化程度加大。
如图2所示为发明的电流原始信号拟合图,对铣削实验过程中采集到的主轴驱动电流信号使用不同阶次傅里叶级数进行拟合对比发现,采用四阶傅里叶级数拟合既能保证拟合准确率,最小化原始数据与电流趋势项的误差平方和,同时又能保证拟合效率,节省数据处理时间,所以原始信号拟合采用的是四阶傅里叶级数依据最小二乘法进行拟合。
如图3为本发明的在不同状态下主轴驱动电流杂波信号图,通过将不同状态下的电流杂波信号对比发现,在金属铣削平稳加工过程中,随着铣刀磨损程度加剧,电流杂波信号波动增大;而在非平稳加工过程中,刀具与工件之间相对振动加剧也会造成动态铣削力增大导致电流杂波波动增大。
如图4为本发明采用的一维卷积神经网络结构图。考虑到本次铣刀磨损状态分类样本有正常磨损、急剧磨损、磨损振动异常和锋利振动异常四类,所以对原有的LeNet-5模型进行了相应的修改。修改后的卷积神经网络由一个输入层、两个卷积层、两个池化层、一个扁平层、一个全连接层和一个输出层组成。其中两层卷积层分别使用的是24个9×1的卷积核和48个9×1的卷积核进行图像特征提取,池化层使用的是6×1的池化核来对上一层数据图缩放,激活函数由sigmoid函数改为Relu函数,池化层的采样方式改为最大池化,加快收敛速度,缓解了梯度消失现象,节省了训练时间。然后使用Softmax分类器根据提取到的刀具磨损状态特征进行分类,最后根据实际分类与样本分类的误差构建损失函数,并利用逐层反向传播算法来优化神经网络的参数,使整个网络达到最好的识别效果。
图5为本发明电流有效值曲线与刀具状态分类图。根据铣刀后刀面磨损量和工件加工表面质量等将刀具状态分为了正常磨损、急剧磨损、磨损振动异常以及锋利振动异常四种。正常磨损和急剧磨损是指在平稳加工过程中发生的后刀面以及刃口的磨损;磨损振动异常是指当刀具磨损严重时,后刀面和刃口出现蹦刃现象导致振动异常;而锋利振动异常指加工时由于切削工艺参数选择不合适导致的切削振动异常现象。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:包括以下步骤,
a.采集加工时主轴驱动电流的原始信号和主轴的振动信号;
b.对主轴驱动电流的原始信号进行数据预处理;
c.将处理后的数据输入到深度学习神经网络中进行训练;
d.将所需数据输入到训练已完成的深度学习神经网络中进行铣刀磨损状态识别;深度学习神经网络采用的式一维卷积神经网络;
深度学习神经网络训练包含网络训练,随机初始化网络参数之后开始逐层训练,得到训练结果;网络优化,根据训练结果与样本标签之间的误差进行误差反向传播调整各层参数,通过不断地训练和调整参数,使得整个网络识别效果达到最优。
2.根据权利要求1所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:铣削加工过程中的铣刀磨损状态分为正常磨损、急剧磨损、磨损振动异常和锋利振动异常四种状态。
3.根据权利要求1所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:主轴驱动电流的原始信号是由基波信号和杂波信号组成。
4.根据权利要求3中所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:基波信号反应的是切削力的准静态分量。
5.根据权利要求4中所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:影响基波信号幅值大小的因素是主轴转速、轴向切深、径向切深、进给速度等与工况相关的加工工艺参数。
6.根据权利要求3中所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:杂波信号反应的是切削力的动态分量,而影响切削力的动态分量的因素是由刀具刃口力、后刀面摩擦力以及切削振动异常共同作用产生的。
7.权利要求6中所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:影响杂波信号幅值大小的因素是不断变化的刃口力系数、后刀面摩擦力系数以及刀具与工件之间的相对振动。
8.根据权利要求1所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:所述的数据预处理至少需要以下步骤:
a.对主轴驱动电流的原始信号进行数据拟合,拟合的结果为主轴驱动电流的基波信号;
b.对主轴驱动电流的原始信号进行数据提取,提取的结果为主轴驱动电流的杂波信号;
c.对提取后的杂波信号进行数据归一化处理。
9.根据权利要求8所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:数据拟合的方法采用四阶傅里叶级数和最小二乘法。
10.根据权利要求8所述的一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法,其特征在于:数据提取的方法是从筛选后主轴驱动电流的原始信号中剔除数据拟合项,获得的数据是与工况无关的主轴驱动电流杂波信号;
数据归一化的方法是将获取到的主轴驱动电流杂波信号进行线性变化,将结果映射到0到1之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010405649.7A CN111660141B (zh) | 2020-05-14 | 2020-05-14 | 一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010405649.7A CN111660141B (zh) | 2020-05-14 | 2020-05-14 | 一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111660141A true CN111660141A (zh) | 2020-09-15 |
CN111660141B CN111660141B (zh) | 2022-02-15 |
Family
ID=72383592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010405649.7A Active CN111660141B (zh) | 2020-05-14 | 2020-05-14 | 一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111660141B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112621382A (zh) * | 2020-12-17 | 2021-04-09 | 北京工业大学 | 基于切削振动阻尼特性的在线刀具磨损监测方法 |
CN113352149A (zh) * | 2021-08-11 | 2021-09-07 | 西安博澜数字科技有限公司 | 一种刀具工况参数监控方法和装置 |
CN113414638A (zh) * | 2021-06-04 | 2021-09-21 | 西北工业大学 | 一种基于铣削力时序图深度学习的变工况铣刀磨损状态预测方法 |
CN113835396A (zh) * | 2021-11-26 | 2021-12-24 | 四川省机械研究设计院(集团)有限公司 | Cnc刀具监测方法和系统及调度管理方法和系统 |
CN114102260A (zh) * | 2021-11-22 | 2022-03-01 | 西安交通大学 | 机理-数据融合驱动的变工况刀具磨损状态监测方法 |
CN114184833A (zh) * | 2021-10-27 | 2022-03-15 | 中国科学院微电子研究所 | 自旋霍尔器件、霍尔电压的获取方法及最大池化的方法 |
CN114227381A (zh) * | 2021-12-24 | 2022-03-25 | 浙江大学台州研究院 | 机床刀具磨损检测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008254080A (ja) * | 2007-04-02 | 2008-10-23 | Mitsubishi Electric Corp | 加工工具寿命検出方法及び加工工具寿命検出装置 |
CN202726639U (zh) * | 2012-08-14 | 2013-02-13 | 沈阳机床(集团)设计研究院有限公司 | 基于数控铣床主轴伺服电机电流信号的切削状态监测系统 |
CN108356606A (zh) * | 2018-03-19 | 2018-08-03 | 西北工业大学 | 基于小波包分析与rbf神经网络的刀具磨损在线监测方法 |
CN108496123A (zh) * | 2016-03-31 | 2018-09-04 | 菲博罗有限公司 | 对机器进行振动诊断监测的方法 |
CN109909804A (zh) * | 2018-12-21 | 2019-06-21 | 北京工业大学 | 基于主轴驱动电流和工步的刀具磨损破损在线监测方法 |
CN111601564A (zh) * | 2018-03-08 | 2020-08-28 | 爱惜康有限责任公司 | 评估超声端部执行器的状态以及用于此的控制系统 |
-
2020
- 2020-05-14 CN CN202010405649.7A patent/CN111660141B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008254080A (ja) * | 2007-04-02 | 2008-10-23 | Mitsubishi Electric Corp | 加工工具寿命検出方法及び加工工具寿命検出装置 |
CN202726639U (zh) * | 2012-08-14 | 2013-02-13 | 沈阳机床(集团)设计研究院有限公司 | 基于数控铣床主轴伺服电机电流信号的切削状态监测系统 |
CN108496123A (zh) * | 2016-03-31 | 2018-09-04 | 菲博罗有限公司 | 对机器进行振动诊断监测的方法 |
CN111601564A (zh) * | 2018-03-08 | 2020-08-28 | 爱惜康有限责任公司 | 评估超声端部执行器的状态以及用于此的控制系统 |
CN108356606A (zh) * | 2018-03-19 | 2018-08-03 | 西北工业大学 | 基于小波包分析与rbf神经网络的刀具磨损在线监测方法 |
CN109909804A (zh) * | 2018-12-21 | 2019-06-21 | 北京工业大学 | 基于主轴驱动电流和工步的刀具磨损破损在线监测方法 |
Non-Patent Citations (1)
Title |
---|
张孟哲: "基于主轴电流的变工况铣刀磨损状态监测研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112621382A (zh) * | 2020-12-17 | 2021-04-09 | 北京工业大学 | 基于切削振动阻尼特性的在线刀具磨损监测方法 |
CN112621382B (zh) * | 2020-12-17 | 2022-05-24 | 北京工业大学 | 基于切削振动阻尼特性的在线刀具磨损监测方法 |
CN113414638A (zh) * | 2021-06-04 | 2021-09-21 | 西北工业大学 | 一种基于铣削力时序图深度学习的变工况铣刀磨损状态预测方法 |
CN113352149A (zh) * | 2021-08-11 | 2021-09-07 | 西安博澜数字科技有限公司 | 一种刀具工况参数监控方法和装置 |
CN113352149B (zh) * | 2021-08-11 | 2021-11-02 | 西安博澜数字科技有限公司 | 一种刀具工况参数监控方法和装置 |
CN114184833A (zh) * | 2021-10-27 | 2022-03-15 | 中国科学院微电子研究所 | 自旋霍尔器件、霍尔电压的获取方法及最大池化的方法 |
CN114102260A (zh) * | 2021-11-22 | 2022-03-01 | 西安交通大学 | 机理-数据融合驱动的变工况刀具磨损状态监测方法 |
CN114102260B (zh) * | 2021-11-22 | 2022-12-09 | 西安交通大学 | 机理-数据融合驱动的变工况刀具磨损状态监测方法 |
CN113835396A (zh) * | 2021-11-26 | 2021-12-24 | 四川省机械研究设计院(集团)有限公司 | Cnc刀具监测方法和系统及调度管理方法和系统 |
CN113835396B (zh) * | 2021-11-26 | 2022-03-04 | 四川省机械研究设计院(集团)有限公司 | Cnc刀具监测方法和系统及调度管理方法和系统 |
CN114227381A (zh) * | 2021-12-24 | 2022-03-25 | 浙江大学台州研究院 | 机床刀具磨损检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111660141B (zh) | 2022-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111660141B (zh) | 一种工况无关的基于主轴驱动电流的铣刀磨损状态识别方法 | |
CN110587377B (zh) | 一种在线监测铣削加工刀具缺损的方法 | |
Zhou et al. | Tool wear condition monitoring in milling process based on current sensors | |
CN108490880B (zh) | 一种数控机床切削刀具磨损状态实时监测方法 | |
CN110059442B (zh) | 一种基于零件表面粗糙度和功率信息的车削加工换刀方法 | |
CN110355608A (zh) | 基于自注意力机制与深度学习的刀具磨损量预测方法 | |
CN112192319A (zh) | 无监督模型的刀具磨损监测方法及系统 | |
Jemielniak | Contemporary challenges in tool condition monitoring | |
CN110909710A (zh) | 一种基于s试件的自适应主轴性能劣化识别方法 | |
CN111958321B (zh) | 基于深度神经网络的数控机床刀具磨损程度识别方法 | |
CN112207631A (zh) | 刀具检测模型的生成方法、检测方法、系统、设备及介质 | |
CN114273974A (zh) | 高速铣削加工中基于振动信号的刀具跳动参数在线估计法 | |
Wolszczak et al. | Monitoring of cutting conditions with the empirical mode decomposition | |
Yang et al. | Application of bispectrum diagonal slice feature analysis in tool wear states monitoring | |
Sun et al. | Tool breakage detection of milling cutter insert based on SVM | |
CN116690313B (zh) | 一种飞机结构件腹板加工刀具失效监控方法 | |
Watanabe et al. | Anomaly detection methods in turning based on motor data analysis | |
CN112946072A (zh) | 基于机器学习的砂带磨损状态监测方法 | |
CN110744359A (zh) | 一种数控车床刀具磨损监测系统及方法 | |
CN114888635A (zh) | 一种刀具状态监测方法 | |
CN114523338A (zh) | 一种基于噪声分析的切削刀具磨损状态监测方法 | |
CN114235043A (zh) | 外圆磨削颤振识别与尾架顶尖力在线监控测量装置及方法 | |
Li et al. | Tool fault diagnosis based on improved multiscale network and feature fusion | |
Reñones et al. | Industrial application of a multitooth tool breakage detection system using spindle motor electrical power consumption | |
CN115971970B (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 |