CN108445835B - 一种数控机床的零件加工尺寸预测方法 - Google Patents
一种数控机床的零件加工尺寸预测方法 Download PDFInfo
- Publication number
- CN108445835B CN108445835B CN201810369668.1A CN201810369668A CN108445835B CN 108445835 B CN108445835 B CN 108445835B CN 201810369668 A CN201810369668 A CN 201810369668A CN 108445835 B CN108445835 B CN 108445835B
- Authority
- CN
- China
- Prior art keywords
- data
- processing
- scale error
- machine tool
- size
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/401—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for measuring, e.g. calibration and initialisation, measuring workpiece for machining purposes
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34242—For measurement only
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
本发明属于数控加工技术领域,并公开了一种数控机床的零件加工尺寸预测方法,该预测方法包括:通过软件与数控系统通讯并采集加工运行中的机床大数据;将机床大数据与铣削加工条件这些数据划分为两部分;将第一部分数据进行图形特征拟合得出图形尺寸数据;将第二部分数据与加工后所用测量仪测量的尺寸误差进行标定并建立预测模型,结合实际将第一部分数据拟合的尺寸误差与第二部分数据建立的预测模型评测的尺寸误差分配权重进行加和,得到最终的在线尺寸误差预测模型。本发明可减少测量仪器的使用,节省人工成本和设备使用维护成本,而且能够方便快捷的在线给出每个加工零件的尺寸误差结果且达到较好的精度,提高生产效率。
Description
技术领域
本发明属于数控加工技术领域,更具体地,涉及一种数控机床的零件加工尺寸预测方法。
背景技术
在机械加工领域,尺寸检验是保证产品质量的重要环节。在生产中,检验尺寸的一般的方法是在加工一定零件件数时抽取部分零件通过量规仪器测量其尺寸,虽然测量精确度高,但检测效率低,且抽样评估整体尺寸的可信度依然存在不合理的问题。
普通的零件尺寸测量方式如公告号为206037926U的专利公开了一种新型尺寸误差检定测量台利用测量小车与工件触碰利用磁栅来测量工件尺寸,测量过程通过感应开关智能控制,有效读数时间长,但仍需要从机床取出工件移动到工作台,会引起二次误差,且装置成本高。而新型的在机检测方式的引进能实现当加工完零件直接调用刀库预置测头并依据机床或额外装置的三坐标定位测量零件尺寸,公告号为CN 206989910 U的专利公开了一种数控加工中心检测装置,其底座的侧面分别焊接有顶部安装有非接触式测量装置的第一连接架和顶部安装有定位误差测量装置的第二连接架,通过定位误差测量装置测量定位误差和补偿、利用非接触测量装置对工件尺寸进行测量,由于该方法是采用“加工后测量”策略,测量与加工不能同时进行、测量成本大,在生产企业很难普及使用。
公告号为CN103235553A综合考虑数控加工系统内各种因素对数控加工精度的影响,从已造成误差的结果出发,通过对加工误差的测量、误差分析、辨识系统阶数,采用适当的迭代学习律计算下一个工件加工前需要施加的误差补偿量,该方法没有对误差源进行检测、建模,测量过程仍然需要在外部进行,不利于测量效率的提高。
综合上述,目前,零件尺寸评测方法存在不能全面评估问题,离线测量效率不高且引起二次装夹误差问题,高精度测量设备成本过高问题,零件尺寸受加工后测量策略瓶颈问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种数控机床的零件加工尺寸预测方法,该预测方法将加工的参数和条件与加工时采集的并按指令域方式提取的过程数据结合起来预测零件的尺寸误差。
为实现上述目的,按照本发明,提供了一种数控机床的零件加工尺寸预测方法,其特征在于,包括以下步骤:
(a)加工前确定数控加工的参数和条件;
(b)加工时按时域采集机床的电控数据;
(c)通过指令域分析方法计算电控数据中的主轴负载电流的特征值,其中,所述指令域分析方法首先在数控加工程序中找出最后一次半精加工或精加工其中一关键尺寸的所在行行号,根据该行号从所述机床总数据中分离出加工关键尺寸的机床数据,然后对该机床数据进行计算得到数据特征值,该数据特征值采用均值或均方根;
(d)将加工过程中的机床总数据划分为两部分:第一部分数据为进给轴位移数据,第二部分数据为G代码任务、电控数据、加工参数和加工条件数据;
(e)将第一部分数据的进给轴位移数据按G代码任务所示的加工图形进行图形拟合,获得拟合图形,将拟合图形的尺寸与加工该尺寸的理论值进行比较,获得尺寸误差;然后将第二部分数据与加工后用测量仪测量获得的尺寸误差进行标定并建立预测模型;
(f)将通过第一部分数据获得的尺寸误差与通过第二部分数据建立的预测模型预测的尺寸误差根据测试实例分配权重进行加权,得到最终的在线尺寸预测模型。
优选地,数控加工的参数和条件包括所选刀具直径、铣削方式、切削参数如主轴转速、进给量、切宽余量、切深和/或零件加工尺寸的类别,其中,零件加工尺寸的类别是指零件的内圆和外圆的尺寸误差的变化趋势(即内圆和外圆是铣大了还是铣小了),铣削方式为顺铣或逆铣。
优选地,电控数据包括主轴负载电流、进给轴电流、进给轴跟随误差和进给轴位移数据。
优选地,所述待预测的零件的关键尺寸包括零件上平行线之间距离尺寸、零件上直径尺寸和/或半径尺寸。
优选地,得到尺寸误差预测模型的具体过程如下:
第一部分数据进行拟合的方式是将采集的指令行的实际坐标插补点按给定的图形进行拟合,得出图形的特征尺寸;
第二部分数据建立预测模型的方式是:加工前的参数和条件结合加工时采集的数控机床电控数据,按实验后用测量仪测得的尺寸误差值通过使用神经网络方法建立非线性映射;
进一步利用这两部分数据获得的多组测试数据,将这些测试数据的结果分配合理的权重进行加权,最后得到尺寸误差预测模型。
优选地,所述图形拟合包括圆拟合和直线特征拟合。
优选地,第二部分数据建立预测模型时,将所选因素与尺寸误差通过建立神经网络模型来进行预测实现。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
1)本发明可减少测量仪器的使用,节省人工成本和设备使用维护成本;
2)本发明能够方便快捷的在线给出每个加工零件的尺寸误差结果且达到较好的精度,提高生产效率;
3)本发明尺寸误差模型易集成到数控系统,方便操作者使用;
4)本发明预测尺寸误差为加工过程补偿提供了判断依据。
附图说明
图1是本发明中加工零件的尺寸评测过程的流程图;
图2是本发明中根据本发明实施例的根据采集的坐标轴插补点位置进行拟合图形并计算尺寸误差的流程图;
图3是本发明中根据第二部分数据建立.神经网络模型预测尺寸误差的结构图;
图4是本发明中机床大数据提取与指令域分析方式的图;
图5是本发明中机床大数据在按指令行显示的特征值;
图6是本发明预测过程的具体流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
图1所示为本发明优选的在数控铣床(或加工中心)加工零件的尺寸误差预测方法。本实例选取对象为宝鸡机床厂的VMC850LT加工中心,数控系统为华中数控HNC-818B。本发明实施包括以下步骤:
步骤一:图1所示,为了使尺寸误差模型更加具体和明确,需要使制造资源一定的情况下即确定的加工环境,夹具工装、机床刚度、工件材料、冷却条件、噪声环境计划保持恒定的前提下尺寸误差仍包含随机性误差和稳定性误差,将模型输出与工作任务结合起来分析。因此在此前提下,明确加工参数和条件:刀具直径、铣削方式、尺寸类别、转速、切深、进给量。
在本实例中,加工参数譬如刀具直径、转速、切深、进给量是以原值代入,而铣削方式按顺铣、逆铣分别编码为1、2,针对尺寸类别,如果是圆尺寸,则按加工外圆与内圆分别编码为1、2,如果是距离尺寸,尺寸类别编码为0。其中,尺寸类别是指零件的内圆和外圆的尺寸误差的变化趋势。
步骤二:当用户在数控机床开始加工按下循环启动时,立即开始通过采集装置采集数控系统的电控数据,该采集装置以1000HZ的采样频率采集数控系统内部电控数据,包括:加工程序当前运行行号、主轴电流、进给轴(X轴、Y轴)电流、进给轴(X轴、Y轴)跟随误差、进给轴(X轴、Y轴)实际位置。
进一步,由于尺寸误差的产生主要发生在最后一次成型加工环节(半精加工或精加工),提取的数据为成型加工尺寸的时间段所采集的数据;
在本实例中采集的数据预处理和分析阶段采用指令域的分析方式,如图4,由于该成型加工对应G代码的某几行,结合G代码行,将采集的过程数据分成与G代码行的对应序列,这样既能保证获取加工成型尺寸期间的过程数据点在时间上对齐;
在本实例中采集的电控数据处理后的是数据在加工某尺寸的G代码行处的特征值,如图5,其中主轴电流的特征值该片段电流的平均值,进给轴电流的特征值取的是均方根值。
在本实例中采集的进给轴(X轴、Y轴)数据处理后仍是原值数据。
步骤三:判断尺寸评测模型是否标定,如果没有标定,需要将多组零件的尺寸在测量仪如三坐标测量机测量,并将实验结果输入评测装置中,则评测装置通过将步骤1和步骤2的数据作为因素,尺寸误差作为因变量进行基本训练得到尺寸评测模型;如果已经标定,则将步骤一和步骤二的数据输入到评测装置,评测装置给出尺寸误差结果。
在本实例中,尺寸误差评测模型将机床的大数据与铣削加工参数和条件这些数据划分为两部分:一部分为G代码任务和进给轴位移数据,另一部分为电流、误差数据、切削参数如进给量、主轴转速、切宽余量、切深和刀具直径数据;
在本实例中,如图3,将第一部分数据即进给轴位移数据结合G代码所示的加工几何图形进行图形特征拟合得出图形尺寸数据;
优选地,如果是图形圆,尺寸为直径或半径,假设采集的数据序列为{Mi(xi,yi),i=1,2,...n},其中n表示采集的数据点个数。则可以通过最小二乘法来代入圆模型方程x2+y2+ax+by+c=0拟合,其中a、b、c为待定参数,将方程改写成Aα=B。
其中
其中,α=(ATA)-1ATB;
则圆直径
尺寸误差y1=d1-d
优选地,如果是距离尺寸,则以其中一条直线L1为基准,加工L1时采集共m个插补点的数据序列为{Pj(x1j,y1j)},j=1,2,...m,拟合给定直线L1,加工与直线L1平行的直线L2时,采集的s个插补点的序列{Qk(x2k,y2k),k=1,2,...s}。则设L1:ax+by+c=0。
采集加工轮廓直线L2上所有点Q1,Q2,Q3,...Qs,对基准直线L1的距离分别为d1,d2,d3,...ds,且平均距离dmean=mean{d1,d2,...ds},则尺寸误差y1=dmean-d。
在本实例中,如图3,将第二部分数据即加工时采集的机床数据如跟随误差、主轴电流等和加工条件如刀具直径、铣削方式(顺铣、逆铣)与加工后所用高精度测量仪测量的尺寸误差进行标定并建立神经网络预测模型。
假如实验选的影响因素为刀具直径、加工方式、切削速度、切削深度、切宽余量、进给量、跟随误差、尺寸类别等,输出为尺寸误差,则建立模型为:
y2=f(d,c,n,ap,aw,f,sc,xc,yc,k,……)
y2——用第二部分数据预测的尺寸误差(μm)
d——为刀具直径(mm)
c——为加工方式,顺铣,逆铣,取值分别为1,2
n——为主轴转速(r/min)
ap——为切削深度(mm)
aw——为切宽余量(mm)
f——为进给量(mm/min)
sc—主轴电流均值(mA)
xc,yc—进给轴(X轴,Y轴)电流均方根(mA)
k——尺寸类别,加工内圆与加工外圆分别取1,2,距离尺寸取为0
在本实例中,考虑实际刀具磨损、切削条件的变化,结合实际将第一部分数据拟合的尺寸误差与第二部分数据建立的预测模型评测的尺寸误差通过额外的测试实例,多个测试实例用三坐标测量的真实结果序列{y0r,r=1,2,…p},其中p表示共p个测试实例,将两部分预测结果序列第一部分数据预测的尺寸误差(y11,y12,y13,...y1p),第二部分数据预测尺寸误差(y21,y22,y23,...y2p),按实际测量值计算权重进行加和,得到最终的在线尺寸预测模型:其中,测试实例是指测试实例是指通过是指通过加工实验完成两种预测方式标定建模后进行再加工实验,实验包括确定加工的参数和条件、采集加工过程电控数据、测量仪测量尺寸。
y=w1·y1+w2·y2
w1,w2为待求变量,令则
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种数控机床的零件加工尺寸预测方法,其特征在于,包括以下步骤:
(a)加工前确定数控加工的参数和条件;
(b)加工时按时域采集机床的电控数据;
(c)通过指令域分析方法计算电控数据中的主轴负载电流的特征值,其中,所述指令域分析方法首先在数控加工程序中找出最后一次半精加工或精加工其中一关键尺寸的所在行行号,根据该行号从机床总数据中分离出加工关键尺寸的机床数据,然后对该机床数据进行计算得到数据特征值,该数据特征值采用均值或均方根;
(d)将加工过程中的机床总数据划分为两部分:第一部分数据为进给轴位移数据,第二部分数据为G代码任务、电控数据、加工参数和加工条件数据;
(e)将第一部分数据的进给轴位移数据按G代码任务所示的加工图形进行图形拟合,获得拟合图形,将拟合图形的尺寸与加工尺寸的理论值进行比较,获得尺寸误差;然后将第二部分数据与加工后用测量仪测量获得的尺寸误差进行标定并建立预测模型;
(f)将通过第一部分数据获得的尺寸误差与通过第二部分数据建立的预测模型预测的尺寸误差根据测试实例分配权重进行加权,得到最终的在线尺寸预测模型。
2.根据权利要求1所述的一种数控机床的零件加工尺寸预测方法,其特征在于,数控加工的参数和条件包括所选刀具直径、铣削方式、主轴转速、进给量、切宽余量、切深和/或零件加工尺寸的类别,其中,零件加工尺寸的类别是指零件的内圆和外圆的尺寸误差的变化趋势,铣削方式为顺铣或逆铣。
3.根据权利要求1所述的一种数控机床的零件加工尺寸预测方法,其特征在于,电控数据包括主轴负载电流、进给轴电流、进给轴跟随误差和进给轴位移数据。
4.根据权利要求1所述的一种数控机床的零件加工尺寸预测方法,其特征在于,所述待预测的零件的关键尺寸包括零件上平行线之间距离尺寸、零件上直径尺寸和/或半径尺寸。
5.根据权利要求1所述的一种数控机床的零件加工尺寸预测方法,其特征在于,得到尺寸误差预测模型的具体过程如下:
第一部分数据进行拟合的方式是将采集的指令行的实际坐标插补点按给定的图形进行拟合,得出图形的特征尺寸;
第二部分数据建立预测模型的方式是:加工前的参数和条件结合加工时采集的数控机床电控数据,按实验后用测量仪测得的尺寸误差值通过使用神经网络方法建立非线性映射;
进一步利用这两部分数据获得的多组测试数据,将这些测试数据的结果分配合理的权重进行加权,最后得到尺寸误差预测模型。
6.根据权利要求1所述的一种数控机床的零件加工尺寸预测方法,其特征在于,所述图形拟合包括圆拟合和直线特征拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810369668.1A CN108445835B (zh) | 2018-04-24 | 2018-04-24 | 一种数控机床的零件加工尺寸预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810369668.1A CN108445835B (zh) | 2018-04-24 | 2018-04-24 | 一种数控机床的零件加工尺寸预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108445835A CN108445835A (zh) | 2018-08-24 |
CN108445835B true CN108445835B (zh) | 2019-11-12 |
Family
ID=63200720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810369668.1A Active CN108445835B (zh) | 2018-04-24 | 2018-04-24 | 一种数控机床的零件加工尺寸预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108445835B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381562B (zh) * | 2018-12-29 | 2023-08-15 | 鸿富锦精密电子(成都)有限公司 | 圆弧拐角加工设备误差补正方法及系统 |
CN109783987A (zh) * | 2019-02-28 | 2019-05-21 | 华中科技大学 | 一种基于神经网络模型的机床误差数据拟合与预测方法 |
CN110716494B (zh) * | 2019-11-13 | 2020-11-17 | 中国航发动力股份有限公司 | 刀具参数识别法和基于刀具参数的摆线加工参数优化方法 |
CN112059726B (zh) * | 2020-08-31 | 2022-05-24 | 西安精雕精密机械工程有限公司 | 一种多基准特征工件的在机测量与补偿加工方法 |
CN113536464B (zh) * | 2021-07-21 | 2024-03-08 | 浙江吉利控股集团有限公司 | 一种基于人工智能的尺寸控制方法及尺寸控制系统 |
CN113591386A (zh) * | 2021-08-03 | 2021-11-02 | 湖州绿色智能制造产业技术研究院 | 一种基于多级结构网络的液压机装配尺寸偏差预测方法 |
CN114310483B (zh) * | 2021-12-13 | 2023-04-28 | 华中科技大学 | 一种数控加工尺寸误差预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777785A (zh) * | 2015-03-02 | 2015-07-15 | 华中科技大学 | 一种基于指令域分析的数控加工工艺参数动态优化方法 |
CN105643370A (zh) * | 2016-03-09 | 2016-06-08 | 镇江智丰自动化科技有限公司 | 一种基于指令域分析的数控机床建模误差补偿方法 |
CN105974886A (zh) * | 2016-06-28 | 2016-09-28 | 华中科技大学 | 一种数控机床的健康监测方法 |
CN106695454A (zh) * | 2016-11-15 | 2017-05-24 | 华中科技大学 | 一种基于指令域分析的数控机床刀库的性能评估方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5691909A (en) * | 1995-12-29 | 1997-11-25 | Western Atlas | Method of virtual machining to predict the accuracy of part to be made with machine tools |
TWI481978B (zh) * | 2012-11-05 | 2015-04-21 | Univ Nat Cheng Kung | 工具機之加工品質的預測方法 |
-
2018
- 2018-04-24 CN CN201810369668.1A patent/CN108445835B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777785A (zh) * | 2015-03-02 | 2015-07-15 | 华中科技大学 | 一种基于指令域分析的数控加工工艺参数动态优化方法 |
CN105643370A (zh) * | 2016-03-09 | 2016-06-08 | 镇江智丰自动化科技有限公司 | 一种基于指令域分析的数控机床建模误差补偿方法 |
CN105974886A (zh) * | 2016-06-28 | 2016-09-28 | 华中科技大学 | 一种数控机床的健康监测方法 |
CN106695454A (zh) * | 2016-11-15 | 2017-05-24 | 华中科技大学 | 一种基于指令域分析的数控机床刀库的性能评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108445835A (zh) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108445835B (zh) | 一种数控机床的零件加工尺寸预测方法 | |
CN109623656B (zh) | 基于厚度在线检测的移动式双机器人协同打磨装置及方法 | |
CN101342664B (zh) | 一种复杂空间型面在机质量检测系统 | |
US6400998B1 (en) | Generation of measurement program in NC machining and machining management based on the measurement program | |
RU2087032C1 (ru) | Способ контроля и оценки качества отливки | |
CN102854841B (zh) | 一种曲面零件的形位误差原位补偿加工方法 | |
CN103777570B (zh) | 基于nurbs曲面的加工误差快速检测补偿方法 | |
CN108655820A (zh) | 一种数控加工基坐标系的标定方法 | |
CN115562161B (zh) | 一种基于在线监测的刀具切削路径加工误差补偿方法 | |
CN114237155A (zh) | 一种多轴数控加工的误差预测及补偿方法、系统及介质 | |
CN107081638A (zh) | 一种机床工作空间热误差辨识模型建立方法 | |
CN104503367A (zh) | 一种基于数控加工装置的涡轮叶片加工及在线检测系统 | |
CN107063149A (zh) | 数控电火花加工设备自动检测系统及方法 | |
CN105388842A (zh) | 一种零件表面加工缺陷的定位方法 | |
US10852709B2 (en) | Machine tool certification for part specific working volume | |
CN110405259A (zh) | 基于多传感器集成测量的自由曲面类零件加工系统 | |
CN108332642B (zh) | 一种直角头精度检测方法 | |
CN107044825A (zh) | 一种关于零件加工尺寸标准认证的检测方法及检测系统 | |
CN108919746B (zh) | 一种转摆台的热误差测试与分析方法 | |
CN112114557B (zh) | 五轴联动数控机床动态精度检测方法、系统及储存介质 | |
CN204423095U (zh) | 一种基于数控加工装置的涡轮叶片加工及在线检测系统 | |
CN110362038A (zh) | 识别五轴联动数控机床在线检测能力的试件与检测方法 | |
CN110524306A (zh) | 一种建立零件和刀具切削用量的匹配关系的系统及方法 | |
CN108080898A (zh) | 一种压气机转子叶片叶身波纹度加工工艺及检测方法 | |
Piórkowski et al. | Statistical testing of milled objects on numerically controlled three-axis milling machines |
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 |