CN116578832A - 预测机械加工零件表面粗糙度的装置 - Google Patents
预测机械加工零件表面粗糙度的装置 Download PDFInfo
- Publication number
- CN116578832A CN116578832A CN202310854707.8A CN202310854707A CN116578832A CN 116578832 A CN116578832 A CN 116578832A CN 202310854707 A CN202310854707 A CN 202310854707A CN 116578832 A CN116578832 A CN 116578832A
- Authority
- CN
- China
- Prior art keywords
- data
- correlation
- surface roughness
- sparrow
- analysis
- 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
- 230000003746 surface roughness Effects 0.000 title claims abstract description 67
- 238000004458 analytical method Methods 0.000 claims abstract description 36
- 238000005520 cutting process Methods 0.000 claims abstract description 31
- 238000010606 normalization Methods 0.000 claims abstract description 21
- 238000003754 machining Methods 0.000 claims abstract description 20
- 238000003801 milling Methods 0.000 claims abstract description 18
- 241000287127 Passeridae Species 0.000 claims description 52
- 238000000034 method Methods 0.000 claims description 41
- 238000012545 processing Methods 0.000 claims description 32
- 238000013527 convolutional neural network Methods 0.000 claims description 30
- 230000008569 process Effects 0.000 claims description 25
- 238000012549 training Methods 0.000 claims description 24
- 238000013528 artificial neural network Methods 0.000 claims description 23
- 238000010219 correlation analysis Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 10
- 230000001360 synchronised effect Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 230000006978 adaptation Effects 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 4
- 210000002569 neuron Anatomy 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000006399 behavior Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000002431 foraging effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007418 data mining Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011176 pooling Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- GJNGXPDXRVXSEH-UHFFFAOYSA-N 4-chlorobenzonitrile Chemical compound ClC1=CC=C(C#N)C=C1 GJNGXPDXRVXSEH-UHFFFAOYSA-N 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004439 roughness measurement Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- 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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/20—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring workpiece characteristics, e.g. contour, dimension, hardness
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0464—Convolutional networks [CNN, ConvNet]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
预测机械加工零件表面粗糙度的装置,通过将振动、切削力、温度的动态变量输入数据转换器中,实现对零件加工过程中表面粗糙度的在线预测,并通过预先设置的粗糙度值优劣区间判断零件表面质量的优劣状况。可降低人工经验的干涉,提高机械加工零件表面粗糙度预测的精度。装置包括:采集模块、分析模块、归一化模块、网络模型、预测模块、铣削数据库平台。
Description
技术领域
本发明涉及机械加工中不规则表面的计量的技术领域,尤其涉及一种预测机械加工零件表面粗糙度的装置。
背景技术
表面粗糙度是衡量机械加工过程性能和评定零件表面质量的一个重要指标。主要影响零件的强度、抗腐蚀性、耐磨性和密封性等性能。随着高端制造领域对零件加工精度的要求越来越高,实际应用中对产品表面质量要求也变高,而表面粗糙度是提高其表面质量的一项非常重要的指标,所以对机械加工(尤其是高速铣削)状态下表面粗糙度的预测显得尤为重要。
表面粗糙度测量方法一般是利用粗糙度测量仪和光传感器,分别通过接触式触针测量和非接触式光学测量,电子显微镜与三维检测技术的应用极大提高了粗糙度预测的精度。但是以上方法都是在零件加工前后对表面粗糙度进行检测,而精准高效地对零件表面粗糙度的在线预测正逐渐成为机械加工领域的发展趋势之一。
随着机械加工数据量的增长,数据挖掘技术和深度神经网络学习逐渐应用到表面粗糙度预测这一领域,并提供了新的研究发展方向。在实际加工过程中,刀具参数、零件参数,切削参数等静态数据,以及力信号、振动信号、温度信号等动态数据都是影响零件表面粗糙度变化的因素。因其高维、非线性的数据特征,使得数据挖掘技术应用于探究异构数据对零件表面粗糙度的影响变化规律显得十分重要,再结合深度神经网络的精准预测,将大大提高加工过程中零件的表面质量。
国内外的研究者主要利用传统的理论方法对零件表面进行粗糙度几何形貌的建模,以及通过数据分析和理论建模方法预测表面粗糙度;利用传统神经网络方法直接对加工后的表面粗糙度进行建模,分析工艺参数对其影响规律,进而预测其表面粗糙度。而利用深度神经网络学习法建立高维的、非线性的动态数据与表面粗糙度之间的关系,相较于传统方法可提高预测的精准度。
发明内容
为克服现有技术的缺陷,本发明要解决的技术问题是提供了一种机械加工零件表面粗糙程度预测的装置,其能够解决传统方法中依赖人工经验提取数据信号特征和一般神经网络预测精度低、效率不高的问题,降低了人工经验的干涉,提高了机械加工零件表面粗糙度预测的精度。
本发明的技术方案是:一种预测机械加工零件表面粗糙度的装置,其包括:
采集模块,其对采集加工过程中的力、振动和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,并将处理后的动态信号数据上传云端;
分析模块,其对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度;
归一化模块,其进行数据的归一化处理;
网络模型,其进行数据归一化后由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型;
预测模块,其通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测;
铣削数据库平台,采集模块、分析模块、归一化模块、网络模型、预测模块的功能运行都依附于云平台,并通过云服务器和云备份为数据处理以及粗糙度预测提供算力与数据同步储存。
本发明通过采集模块对采集加工过程中的力、振动和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,分析模块对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度,进行数据的归一化处理,网络模型由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型,预测模块通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测,因此能够解决传统方法中依赖人工经验提取数据信号特征和一般神经网络预测精度低、效率不高的问题,降低了人工经验的干涉,提高了机械加工零件表面粗糙度预测的精度。并且全程通过连接云端服务器,以提高数据处理效率和神经网络运算速度,以及粗糙度实时在线预测响应速率,并在数据处理分析过程时刻将数据同步并上传云平台备份,以增强海量铣削数据的存储能力。
附图说明
图1是根据本发明的预测机械加工零件表面粗糙度的装置工作原理图。
图2是灰色关联度分析的流程图。
图3是相关性分析流程图。
图4是SSA-CNN算法的工作原理图。
图5是表面粗糙度的预测值与实际值对比图。
图6是表面粗糙度的误差直方图。
图7是表面粗糙度的回归拟合图。
图8是均方根误差收敛图。
具体实施方式
如图1所示,一种预测机械加工零件表面粗糙度的装置,其包括:
采集模块,其对采集加工过程中的力、振动和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,并将处理后的动态信号数据上传云端;
分析模块,其对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度;
归一化模块,其进行数据的归一化处理;
网络模型,其进行数据归一化后由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型;
预测模块,其通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测;
铣削数据库平台,采集模块、分析模块、归一化模块、网络模型、预测模块的功能运行都依附于云平台,并通过云服务器和云备份为数据处理以及粗糙度预测提供算力与数据同步储存。
本发明通过采集模块对采集加工过程中的振动、切削力和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,分析模块对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度,进行数据的归一化处理,网络模型由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型,预测模块通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测,因此能够解决传统方法中依赖人工经验提取数据信号特征和一般神经网络预测精度低、效率不高的问题,降低了人工经验的干涉,提高了机械加工零件表面粗糙度预测的精度。并且全程通过连接云端服务器,以提高数据处理效率和神经网络运算速度,以及粗糙度实时在线预测响应速率,并在数据处理分析过程时刻将数据同步并上传云平台备份,以增强海量铣削数据的存储能力。
优选地,所述采集模块中,将压电式三向测力仪固定于零件底部;选取三向加速度传感器分别安装在零件旁和机械加工刀具旁;将红外热像仪与笔记本相连接,并放在合适位置,聚焦于零件;切削力信号通过Kistler三向测力仪采集;振动信号通过LMS振动噪声测试分析仪来进行信号采集、整理和分析;温度信号通过Optris测温仪对温度信息进行采集与处理;最后通过多传感器采集系统对三种信号进行汇总,并实现同步采集。
优选地,所述分析模块中,对采集的振动、切削力、温度信号进行关联度分析为:
通过指标正向化处理得到矩阵X:
,
将母序列,系统的因变量Y,记为x0;将子序列,系统的自变量X,记为[x1,x2,…,xn],
母、子序列中各个指标的元素分别除以对应指标的均值,得到的结果记为zmn:
,
式中xmn为序列中的指标元素,为指标均值,
排列后得到标准化矩阵,记为Z:
,
通过以下公式计算子序列中各个指标与母序列的关联系数,
,
其中ρ为分辨系数,a为两极最小差,b为两极最大差,a,b=1,2,3,…,N; a、b取值如下式:
,
定义y(x0,xi) 为灰色关联度,并通过以下公式求取每列关联系数矩的均值,以此计算灰色关联度,确定关联程度,
,
比较灰色关联度,将子序列中影响程度大,关联度高的因素筛出。
优选地,所述分析模块中,对采集的振动、切削力、温度信号进行相关性分析为:
在R语言环境下使用统计软件SPSS,通过正态性检验来研究机械加工动态信号数据的各类影响因素对表面粗糙度的相关性显著程度,
构建两组数据变量矩阵散点图,判断两组数据变量之间是否有线性关系,对分析的两组数据进行一个正态性检验,判断是否服从正态分布,若两组连续数据服从正态分布且二者属于线性相关关系则利用Pearson系数,反之则使用Spearman系数;
选择Pearson系数进行相关性分析,计算公式如下:
,
其中ρ表示相关系数,cov表示协方差,E表示数学期望/均值;
选择Spearman系数进行相关性分析,计算公式如下:
,
其中d i 表示第i个数据对应的位次之差,n表示总的观测样本数;
相关性系数取值范围为[-1,1],0表示两者之间没有相关性,-1表示负相关性,1表示正相关性,越接近-1或1,则相关性越强;将相关性强的数据筛选出来。
优选地,所述归一化模块中,
将机械加工过程中的动态信号分组进行归一化处理,使数据在[0,1]区间,依据公式如下:
,
式中X'为此组数据归一化后的结果,X max 和X min 分别为此组数据的最大值和最小值,通过MATLAB中的随机函数randperm将归一化的数据打乱排序,再将数据以9:1的比例划分训练集和测试集。
优选地,所述网络模型中,
确定麻雀算法的种群数、最大迭代次数和初始参数的优化个数,并分别对其优化边界进行设置;
种群数和最大迭代次数设置为10;优化参数包括学习率e、迭代次数n、每批迭代数minibatch、第一层卷积层的核大小ks1和数量kn1、第2层卷积层的核大小ks2和数量kn2,以及两个全连接层的神经元数量N1、N2,共九个参数需要优化;参数的优化范围分别是, e∈[0.001-0.01]、n∈[10-100]、minibatch∈[100-512]、ks∈[1-3]、kn∈[1-20]、N∈[1-50];
确定发现者在全部种群中的占比,
使麻雀种群中的发现者占麻雀群总数的20%;
初始化种群并更新各群体的位置;
发现者位置更新公式:
,
式中,k为迭代次数,表示第i个成员在第j维所处的位置,i,j=1,2,3,…,N;α∈(0,1]且为均匀分布的随机数;M为最大迭代次数;Q为服从标准正态分布的随机数;L为1×d的单位矩阵;R 2 为预警值,范围在[0,1];ST为安全值,范围在[0.5,1];
当R 2 <ST,则麻雀群搜索范围附近没有危险,可继续扩大搜索范围,反之则有危险,需要移动位置;
追随者位置更新公式:
,
其中,表示发现者在k次迭代时所处的最差位置;/>表示发现者在k+1次迭代的最优位置;A为1×d的多维矩阵,其元素为1或-1的随机数,A + =A T (AA T )-1;
当 i>n/2,则表示麻雀群的适应度低,需要扩大搜索范围,继续移动觅食,反之则适应度值高,可在最优位置附近随机觅食;
预警者位置更新公式:
,
其中 ,表示发现者在k次迭代时所处的最优位置;β为步长参数,且服从N(0,1);/>为[-1,1]的随机数,体现麻雀群的移动方向,并控制步长;ε避免分母为零所设的极小值常数;f i 为第i个麻雀的适应度值;f w 为麻雀群的最差适应度值;f g 为麻雀群的最优适应度值;
通过觅食与反捕食行为,更新麻雀最优个体的位置,当符合RMSE小于设定的误差阈值或大于迭代次数时,输出最优个体所对应的最优参数;
将SSA优化后的初始参数输入CNN中,进行神经网络训练,当误差在允许内,则输出已训练好的模型,并对表面粗糙度Ra进行预测。
需进行说明的是,本文并不是直接预测工件表面粗糙度Ra,而是通过预测峰谷任一测点到基准线距离的绝对值后,再通过粗糙度计算公式计算出Ra,以达到预测粗糙度的目的。
Ra的公式如下:
,
式中n表示测点数;Z(x)表示被测轮廓线上各点至基准线距离的绝对值。
还提供了一种机械加工零件表面的粗糙程度计量预测的方法,其包括以下步骤:
(1)对采集加工过程中的力、振动和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,并将处理后的动态信号数据上传云端;
(2)对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度;
(3)数据经归一化后由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型;
(4)通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测;
(5)将数据同步并上传云平台,并进行备份。
优选地,所述步骤(1)中,将压电式三向测力仪固定于零件底部;选取三向加速度传感器分别安装在零件旁和机械加工刀具旁;将红外热像仪与笔记本相连接,并放在合适位置,聚焦于零件;切削力信号通过Kistler三向测力仪采集;振动信号通过LMS振动噪声测试分析仪来进行信号采集、整理和分析;温度信号通过Optris测温仪对温度信息进行采集与处理;最后通过多传感器采集系统对三种信号进行汇总,并实现同步采集。
优选地,如图2所示,所述步骤(3)中,对采集的振动、切削力、温度信号进行关联度分析为:
通过指标正向化处理得到矩阵X:
,
将母序列,系统的因变量Y,记为x0;将子序列,系统的自变量X,记为[x1,x2,…,xn],
母、子序列中各个指标的元素分别除以对应指标的均值,得到的结果记为zmn:
,
式中xmn为序列中的指标元素,为指标均值,
排列后得到标准化矩阵,记为Z:
,
通过以下公式计算子序列中各个指标与母序列的关联系数,
,
其中a为两极最小差,b为两极最大差,ρ为分辨系数,i,k=1,2,3,…,N; a、b取值如下式:
,
定义y(x0,xi) 为灰色关联度,并通过以下公式求取每列关联系数矩的均值,以此计算灰色关联度,确定关联程度,
,
比较灰色关联度,将子序列中影响程度大,关联度高的因素筛出。
优选地,如图3所示,所述步骤(3)中,对采集的振动、切削力、温度信号进行相关性分析为:
在R语言环境下使用统计软件SPSS,通过正态性检验来研究机械加工动态信号数据的各类影响因素对表面粗糙度的相关性显著程度,
构建两组数据变量矩阵散点图,判断两组数据变量之间是否有线性关系,对分析的两组数据进行一个正态性检验,判断是否服从正态分布,若两组连续数据服从正态分布且二者属于线性相关关系则利用Pearson系数,反之则使用Spearman系数;
选择Pearson系数进行相关性分析,计算公式如下:
,
其中ρ表示相关系数,cov表示协方差,E表示数学期望/均值;
选择Spearman系数进行相关性分析,计算公式如下:
,
其中d i 表示第i个数据对应的位次之差,n表示总的观测样本数;
可对生成的相关性系数矩阵进行可生成视化操作,并成热力图。
表1
如表1所示,对于相关性系数来说,0表示两者之间没有相关性,-1表示负相关性,1表示正相关性,越接近-1或1,则相关性越强。将相关性强的数据筛选出来。
优选地,所述步骤(4)中,
将机械加工过程中的动态信号分组进行归一化处理,使数据在[0,1]区间,依据公式如下:
,
式中X'为此组数据归一化后的结果,X max 和X min 分别为此组数据的最大值和最小值,通过MATLAB中的随机函数randperm将归一化的数据打乱排序,再将数据以9:1的比例划分训练集和测试集。
优选地,如图4所示,所述步骤(5)中,
确定麻雀算法的种群数、最大迭代次数和初始参数的优化个数,并分别对其优化边界进行设置;
种群数和最大迭代次数设置为10;优化参数包括学习率e、迭代次数n、每批迭代数minibatch、第一层卷积层的核大小ks1和数量kn1、第2层卷积层的核大小ks2和数量kn2,以及两个全连接层的神经元数量N1、N2,共九个参数需要优化;参数的优化范围分别是, e∈[0.001-0.01]、n∈[10-100]、minibatch∈[100-512]、ks∈[1-3]、kn∈[1-20]、N∈[1-50];
确定发现者在全部种群中的占比,
使麻雀种群中的发现者占麻雀群总数的20%;
初始化种群并更新各群体的位置;
发现者位置更新公式:
,
式中,k为迭代次数,表示第i个成员在第j维所处的位置,i,j=1,2,3,…,N;α∈(0,1]且为均匀分布的随机数;M为最大迭代次数;Q为服从标准正态分布的随机数;L为1×d的单位矩阵;R 2 为预警值,范围在[0,1];ST为安全值,范围在[0.5,1];
当R 2 <ST,则麻雀群搜索范围附近没有危险,可继续扩大搜索范围,反之则有危险,需要移动位置;
追随者位置更新公式:
,
其中,表示发现者在k次迭代时所处的最差位置;/>表示发现者在k+1次迭代的最优位置;A为1×d的多维矩阵,其元素为1或-1的随机数,A + =A T (AA T )-1;
当i>n/2,则表示麻雀群的适应度低,需要扩大搜索范围,继续移动觅食,反之则适应度值高,可在最优位置附近随机觅食;
预警者位置更新公式:
,
其中, 表示发现者在k次迭代时所处的最优位置;β为步长参数,且服从N(0,1);/>为[-1,1]的随机数,体现麻雀群的移动方向,并控制步长;ε避免分母为零所设的极小值常数;f i 为第i个麻雀的适应度值;f w 为麻雀群的最差适应度值;f g 为麻雀群的最优适应度值;
通过觅食与反捕食行为,更新麻雀最优个体的位置,当符合RMSE小于设定的误差阈值或大于迭代次数时,输出最优个体所对应的最优参数;
将SSA优化后的初始参数输入CNN中,进行神经网络训练,当RMSE小于设定的误差阈值或大于迭代次数时,如图8所示当误差在允许范围RMSE=0.1以下并趋于平稳时,则输出已训练好的模型,并对表面粗糙度Ra进行预测。
需进行说明的是,本文并不是直接预测工件表面粗糙度Ra,而是通过预测峰谷任一测点到基准线距离的绝对值后,再通过粗糙度计算公式计算出Ra,以达到预测粗糙度的目的。
Ra的公式如下:
,
式中n表示测点数;Z(x)表示被测轮廓线上各点至基准线距离的绝对值。
下面以一台三轴立式数控机床铣削加工为例,详细说明本发明的实施方式。
第一步,准备实验器材:
实验所用零件尺寸为100mm×60mm×25mm的TC4钛合金;实验刀具选用PCBN铣刀片,选用直径14cm、16cm和18cm的两刃硬质合金铣刀杆;机床类型为三轴立式数控铣削机床;加工方式采用顺铣干切。
第二步,铣削动态数据信号的采集:
将Kistler 9129AA压电式三向测力仪固定于零件底部;选取PCB的356A16 SN102355型号的三向加速度传感器安装在零件旁和铣削刀具旁,并连接到LMS振动测试仪;将Optris PI红外热像仪与笔记本相连接,并置于合适位置,调整镜头角度和焦距。铣削力信号通过Kistler 5080A电荷放大器,由Kistler 5697A数据采集器采集;振动信号通过软件LMS Test.Lab来进行信号采集、整理和分析;温度信号通过Optris PI Connect软件对温度信息进行处理。将采集频率统一设为4096Hz,采集时间为10s,最后通过多传感器采集系统对三种信号进行汇总,并实现同步采集。
第三步,对采集信号数据的筛选处理:
对采集的力、振动、温度信号进行关联度和相关性分析。
采用灰色关联度分析法分析各影响因素关联程度的大小,找出影响表面粗糙度变化的主要因素 。
对信号数据的相关性分析。
对关联程度高的影响因素所对应的数据进行相关性分析,在R语言环境下使用统计软件SPSS,通过正态性检验来研究铣削动态信号数据的各类影响因素对表面粗糙度的相关性显著程度。
第四步,数据的归一化处理:
将铣削过程中的动态信号分组进行归一化处理,使数据在[0,1]区间。将归一化的数据通过MATLAB中的随机函数randperm打乱排序,再将数据以9:1的比例划分训练集和测试集。其中,训练集样本数为1800,测试集样本数为200。
第五步,构建神经网络结构和训练模型:
确定麻雀算法的种群数、最大迭代次数和初始参数的优化个数,并分别对其优化边界进行设置。
种群数和最大迭代次数设置为10;优化参数包括学习率e、迭代次数n、每批迭代数minibatch、第一层卷积层的核大小ks1和数量kn1、第2层卷积层的核大小ks2和数量kn2,以及两个全连接层的神经元数量N1、N2,共九个参数需要优化;参数的优化范围分别是, e∈[0.001-0.01]、n∈[10-100]、minibatch∈[100-512]、ks∈[1-3]、kn∈[1-20]、N∈[1-50]。其中因池化层会在一定程度上过滤掉显著的有用信息,故在此不添加池化层。
使麻雀种群中的发现者占麻雀群总数的20%,并初始化种群并更新各群体的位置,进而输出最优位置所对应的超参数。优化后的各超参数分别是为,e=0.0062、n=30、minibitch=112、ks1=9、kn1=288、ks2=13、kn2=177、N1=42、N2=26。
通过卷积神经网络自动提取动态特征,并利用优化后的参数训练神经网络,经过两个卷积层,两个全连接层,辍学率为0.5进行拟合回归输出预测值。最终预测准确率达到97.34%,均方根误差为0.010384。
由粗糙度模型的预测结果可绘制出表面粗糙度预测值与实际值对比图,通过SSA-CNN神经网络预测误差图观察可知,预测值与真实值的误差范围控制在±0.04μm的较小范围内。同时表面粗糙度预测误差直方图(参见图6)以零误差为界限向两边呈现出良好的正态分布。并且当表面粗糙度回归拟合(参见图7)时,R=0.98672,且决定系数R2=0.973,此时均方根误差RMSE=0.0108786。
如图5所示,以CNN的预测值与实际值的均方误差最小化作为适应度函数,通过训练得到的CNN的误差能够最小化,且预测准确率在90%以上,则输出对应的神经网络模型。
第六步,零件铣削过程表面粗糙度的实时预测:
在零件高速铣削过程中,将采集到的动态参数经预处理和筛选后输入SSA-CNN模型中,并进行动态数据特征的自动提取。随后通过神经网络预测模型训练并输出此工况下零件表面粗糙度的实时预测值,进而反映零件表面质量的优劣状况。
第七步,上传实验数据和结果:
实验全程连接云端服务器,并在数据处理分析过程中,时刻将实验数据和预测结果同步并上传云平台备份,以提高运算处理速度、增强存储能力。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属本发明技术方案的保护范围。
Claims (6)
1.预测机械加工零件表面粗糙度的装置,其特征在于:其包括:
采集模块,其对采集加工过程中的力、振动和温度信号进行预处理,通过高低通滤波去除波动较大的异常频段,并将处理后的动态信号数据上传云端;
分析模块,其对处理后的振动、切削力、温度信号进行关联度和相关性分析,利用灰色关联度分析法对数据进行关联度分析,找出影响表面粗糙度变化的主要因素,之后利用Pearson和Spearman作为评估标准分析振动、切削力、温度和表面粗糙度之间的相关性显著程度;
归一化模块,其进行数据的归一化处理;
网络模型,其进行数据归一化后由麻雀群优化算法搭建卷积神经网络模型,包括神经网络结构设计、特征自动提取与迭代训练,训练过程中,当评估误差在规定范围内获取优良的卷积神经网络模型;
预测模块,其通过在线预测将振动、切削力、温度的动态数据特征值输入优良的卷积神经网络模型,实现对零件加工过程中表面粗糙度的在线预测;
铣削数据库平台,采集模块、分析模块、归一化模块、网络模型、预测模块的功能运行都依附于云平台,并通过云服务器和云备份为数据处理以及粗糙度预测提供算力与数据同步储存。
2.根据权利要求1所述的预测机械加工零件表面粗糙度的装置,其特征在于:所述采集模块中,将压电式三向测力仪固定于零件底部;选取三向加速度传感器分别安装在零件旁和机械加工刀具旁;将红外热像仪与笔记本相连接,并放在合适位置,聚焦于零件;切削力信号通过Kistler三向测力仪采集;振动信号通过LMS振动噪声测试分析仪来进行信号采集、整理和分析;温度信号通过Optris测温仪对温度信息进行采集与处理;最后通过多传感器采集系统对三种信号进行汇总,并实现同步采集。
3.根据权利要求2所述的预测机械加工零件表面粗糙度的装置,其特征在于:所述分析模块中,对采集的振动、切削力、温度信号进行关联度分析为:
通过指标正向化处理得到矩阵X:
,
将母序列,系统的因变量Y,记为x0;将子序列,系统的自变量X,记为[x1,x2,…,xn],
母、子序列中各个指标的元素分别除以对应指标的均值,得到的结果记为zmn:
,
式中xmn为序列中的指标元素,为指标均值,
排列后得到标准化矩阵,记为Z:
,
通过以下公式计算子序列中各个指标与母序列的关联系数,
,
其中ρ为分辨系数,a为两极最小差,b为两极最大差,a,b=1,2,3,…,N; a、b取值如下式:
,
定义y(x0,xi) 为灰色关联度,并通过以下公式求取每列关联系数矩的均值,以此计算灰色关联度,确定关联程度,
,
比较灰色关联度,将子序列中影响程度大,关联度高的因素筛出。
4.根据权利要求3所述的预测机械加工零件表面粗糙度的装置,其特征在于:所述分析模块中,对采集的振动、切削力、温度信号进行相关性分析为:
在R语言环境下使用统计软件SPSS,通过正态性检验来研究机械加工动态信号数据的各类影响因素对表面粗糙度的相关性显著程度,
构建两组数据变量矩阵散点图,判断两组数据变量之间是否有线性关系,对分析的两组数据进行一个正态性检验,判断是否服从正态分布,若两组连续数据服从正态分布且二者属于线性相关关系则利用Pearson系数,反之则使用Spearman系数;
选择Pearson系数进行相关性分析,计算公式如下:
,
其中ρ表示相关系数,cov表示协方差,E表示数学期望/均值;
选择Spearman系数进行相关性分析,计算公式如下:
,
其中d i 表示第i个数据对应的位次之差,n表示总的观测样本数;
相关性系数取值范围为[-1,1],0表示两者之间没有相关性,-1表示负相关性,1表示正相关性,越接近-1或1,则相关性越强;将相关性强的数据筛选出来。
5.根据权利要求4所述的预测机械加工零件表面粗糙度的装置,其特征在于:所述归一化模块中,
将机械加工过程中的动态信号分组进行归一化处理,使数据在[0,1]区间,依据公式如下:
,
式中X'为此组数据归一化后的结果,X max 和X min 分别为此组数据的最大值和最小值,通过MATLAB中的随机函数randperm将归一化的数据打乱排序,再将数据以9:1的比例划分训练集和测试集。
6.根据权利要求5所述的预测机械加工零件表面粗糙度的装置,其特征在于:所述网络模型中,
确定麻雀算法的种群数、最大迭代次数和初始参数的优化个数,并分别对其优化边界进行设置;
种群数和最大迭代次数设置为10;优化参数包括学习率e、迭代次数n、每批迭代数minibatch、第一层卷积层的核大小ks1和数量kn1、第2层卷积层的核大小ks2和数量kn2,以及两个全连接层的神经元数量N1、N2,共九个参数需要优化;参数的优化范围分别是, e∈[0.001-0.01]、n∈[10-100]、minibatch∈[100-512]、ks∈[1-3]、kn∈[1-20]、N∈[1-50];
确定发现者在全部种群中的占比,
使麻雀种群中的发现者占麻雀群总数的20%;
初始化种群并更新各群体的位置;
发现者位置更新公式:
,
式中,k为迭代次数,表示第i个成员在第j维所处的位置,i,j=1,2,3,…,N;α∈(0,1]且为均匀分布的随机数;M为最大迭代次数;Q为服从标准正态分布的随机数;L为1×d的单位矩阵;R 2 为预警值,范围在[0,1];ST为安全值,范围在[0.5,1];
当R 2 <ST,则麻雀群搜索范围附近没有危险,可继续扩大搜索范围,反之则有危险,需要移动位置;
追随者位置更新公式:
,
其中,表示发现者在k次迭代时所处的最差位置;/>表示发现者在k+1次迭代的最优位置;A为1×d的多维矩阵,其元素为1或-1的随机数,A + =A T (AA T )-1;
当 i>n/2,则表示麻雀群的适应度低,需要扩大搜索范围,继续移动觅食,反之则适应度值高,可在最优位置附近随机觅食;
预警者位置更新公式:
,
其中 ,表示发现者在k次迭代时所处的最优位置;β为步长参数,且服从N(0,1);/>为[-1,1]的随机数,体现麻雀群的移动方向,并控制步长;ε避免分母为零所设的极小值常数;f i 为第i个麻雀的适应度值;f w 为麻雀群的最差适应度值;f g 为麻雀群的最优适应度值;
通过觅食与反捕食行为,更新麻雀最优个体的位置,当符合RMSE小于设定的误差阈值或大于迭代次数时,输出最优个体所对应的最优参数;
将SSA优化后的初始参数输入CNN中,进行神经网络训练,当误差在允许内,则输出已训练好的模型,并对表面粗糙度Ra进行预测,通过预测峰谷任一测点到基准线距离的绝对值后,再通过粗糙度计算公式计算出Ra,以达到预测粗糙度的目的,公式如下:
,
式中n表示测点数;Z(x)表示被测轮廓线上各点至基准线距离的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310854707.8A CN116578832B (zh) | 2023-07-13 | 2023-07-13 | 预测机械加工零件表面粗糙度的装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310854707.8A CN116578832B (zh) | 2023-07-13 | 2023-07-13 | 预测机械加工零件表面粗糙度的装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116578832A true CN116578832A (zh) | 2023-08-11 |
CN116578832B CN116578832B (zh) | 2023-09-12 |
Family
ID=87534503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310854707.8A Active CN116578832B (zh) | 2023-07-13 | 2023-07-13 | 预测机械加工零件表面粗糙度的装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116578832B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117495211A (zh) * | 2024-01-03 | 2024-02-02 | 东北大学 | 基于自适应周期发现的工业母机加工工件质量预测方法 |
CN117520825A (zh) * | 2024-01-04 | 2024-02-06 | 东北大学 | 基于多尺度特征融合的工业母机加工工件质量预测方法 |
CN117520825B (zh) * | 2024-01-04 | 2024-05-10 | 东北大学 | 基于多尺度特征融合的工业母机加工工件质量预测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2012145217A (ru) * | 2012-10-25 | 2014-04-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Самарский государственный технический университет" | Способ оценки шероховатости поверхности детали и устройство для его реализации |
CN103761429A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 铣削加工工件表面粗糙度的预测方法 |
CN106407669A (zh) * | 2016-09-07 | 2017-02-15 | 江苏大学 | 一种切削加工表面粗糙度的预测方法 |
CN111177911A (zh) * | 2019-12-24 | 2020-05-19 | 大连理工大学 | 一种基于sdae-dbn算法的零件表面粗糙度在线预测方法 |
CN111366123A (zh) * | 2020-03-06 | 2020-07-03 | 大连理工大学 | 一种基于多任务学习的零件表面粗糙度和刀具磨损预测方法 |
CN113065207A (zh) * | 2021-03-25 | 2021-07-02 | 沈阳工业大学 | 一种砂带磨削螺杆曲面表面粗糙度的预测方法 |
US20210364482A1 (en) * | 2020-03-06 | 2021-11-25 | Dalian University Of Technology | Prediction method of part surface roughness and tool wear based on multi-task learning |
CN115655193A (zh) * | 2022-09-27 | 2023-01-31 | 沈阳工业大学 | 一种铣削工件表面粗糙度的预测方法 |
CN116117598A (zh) * | 2023-01-16 | 2023-05-16 | 江苏师范大学 | 湿铣削条件下的表面质量测量与预测系统 |
CN116186499A (zh) * | 2023-02-23 | 2023-05-30 | 武汉数字化设计与制造创新中心有限公司 | 一种数控加工表面粗糙度预测方法 |
-
2023
- 2023-07-13 CN CN202310854707.8A patent/CN116578832B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2012145217A (ru) * | 2012-10-25 | 2014-04-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Самарский государственный технический университет" | Способ оценки шероховатости поверхности детали и устройство для его реализации |
CN103761429A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 铣削加工工件表面粗糙度的预测方法 |
CN106407669A (zh) * | 2016-09-07 | 2017-02-15 | 江苏大学 | 一种切削加工表面粗糙度的预测方法 |
CN111177911A (zh) * | 2019-12-24 | 2020-05-19 | 大连理工大学 | 一种基于sdae-dbn算法的零件表面粗糙度在线预测方法 |
CN111366123A (zh) * | 2020-03-06 | 2020-07-03 | 大连理工大学 | 一种基于多任务学习的零件表面粗糙度和刀具磨损预测方法 |
US20210364482A1 (en) * | 2020-03-06 | 2021-11-25 | Dalian University Of Technology | Prediction method of part surface roughness and tool wear based on multi-task learning |
CN113065207A (zh) * | 2021-03-25 | 2021-07-02 | 沈阳工业大学 | 一种砂带磨削螺杆曲面表面粗糙度的预测方法 |
CN115655193A (zh) * | 2022-09-27 | 2023-01-31 | 沈阳工业大学 | 一种铣削工件表面粗糙度的预测方法 |
CN116117598A (zh) * | 2023-01-16 | 2023-05-16 | 江苏师范大学 | 湿铣削条件下的表面质量测量与预测系统 |
CN116186499A (zh) * | 2023-02-23 | 2023-05-30 | 武汉数字化设计与制造创新中心有限公司 | 一种数控加工表面粗糙度预测方法 |
Non-Patent Citations (3)
Title |
---|
张宝磊;熊艺文;王为庆;朱帅玲;: "高速铣削TC4表面粗糙度预测模型研究", 组合机床与自动化加工技术, no. 03 * |
王卓 等: "光学材料磨削加工亚表面损伤层深度测量及预测方法研究", 航空精密制造技术, vol. 43, no. 5 * |
闵莉;王哲;高龙飞;刘大任;: "基于BP神经网络的车削表面粗糙度检测", 沈阳建筑大学学报(自然科学版), no. 05 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117495211A (zh) * | 2024-01-03 | 2024-02-02 | 东北大学 | 基于自适应周期发现的工业母机加工工件质量预测方法 |
CN117495211B (zh) * | 2024-01-03 | 2024-03-19 | 东北大学 | 基于自适应周期发现的工业母机加工工件质量预测方法 |
CN117520825A (zh) * | 2024-01-04 | 2024-02-06 | 东北大学 | 基于多尺度特征融合的工业母机加工工件质量预测方法 |
CN117520825B (zh) * | 2024-01-04 | 2024-05-10 | 东北大学 | 基于多尺度特征融合的工业母机加工工件质量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116578832B (zh) | 2023-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108803486B (zh) | 基于并联深度学习网络的数控机床热误差预测与补偿方法 | |
CN110647943B (zh) | 基于演化数据聚类分析的切削刀具磨损监测方法 | |
CN116578832B (zh) | 预测机械加工零件表面粗糙度的装置 | |
CN106842922B (zh) | 一种数控加工误差优化方法 | |
CN112247674B (zh) | 一种刀具磨损预测方法 | |
CN113688534B (zh) | 一种基于多特征融合模型寻找最优铣削参数的研究方法 | |
CN113487533B (zh) | 一种基于机器学习的零件装配质量数字化检测系统及方法 | |
CN110889091A (zh) | 基于温度敏感区间分段的机床热误差预测方法与系统 | |
CN113609789A (zh) | 一种基于时空特征并行提取的刀具磨损状态预测方法 | |
CN111881860B (zh) | 滚刀磨损在位识别模型建模方法及滚刀磨损在位识别方法 | |
Li et al. | Roughness prediction model of milling noise-vibration-surface texture multi-dimensional feature fusion for N6 nickel metal | |
CN111650891B (zh) | 具有恒力自适应控制方法的五轴精密小龙门数控加工中心 | |
CN114161227A (zh) | 一种基于仿真特征和信号特征融合的刀具磨损量监控方法 | |
Kumar et al. | Performance evaluation for tool wear prediction based on Bi-directional, Encoder–Decoder and Hybrid Long Short-Term Memory models | |
CN112475410A (zh) | 一种铣削温度与多元影响因子的关联分析系统及方法 | |
CN113601261B (zh) | 一种面向刀具的在线快速优化模型的监测方法 | |
CN114186337A (zh) | 一种基于多源数据融合的压气机旋转失速预测方法 | |
CN113836662A (zh) | 凸轮曲线槽机构设计缺陷的动态辨识及去特征化修复方法 | |
CN115099266A (zh) | 一种基于梯度提升决策树的硬车表面白层预测方法 | |
Liu et al. | Selection of key temperature measuring points for thermal error modeling of CNC machine tools | |
CN117592223B (zh) | 一种面向航空航天材料的孔加工刀具智能设计方法 | |
Zou et al. | Development of a tool wear online monitoring system for dry gear hobbing machine based on new experimental approach and DAE-BPNN-integrated mathematic structure | |
CN111580463A (zh) | 机床热关键点选取方法 | |
Lan et al. | Study on the Predictive of Dynamic Milling Force of Milling Process Based on Data Mining | |
KR102558133B1 (ko) | 인공지능을 통한 가공 설비 최적 운전 조건 추적 방법 |
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 |