CN113569903A - 数控机床刀具磨损预测方法、系统、设备、介质、终端 - Google Patents
数控机床刀具磨损预测方法、系统、设备、介质、终端 Download PDFInfo
- Publication number
- CN113569903A CN113569903A CN202110642404.0A CN202110642404A CN113569903A CN 113569903 A CN113569903 A CN 113569903A CN 202110642404 A CN202110642404 A CN 202110642404A CN 113569903 A CN113569903 A CN 113569903A
- Authority
- CN
- China
- Prior art keywords
- signal
- data set
- numerical control
- signals
- control machine
- 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
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000005299 abrasion Methods 0.000 title claims description 23
- 230000036541 health Effects 0.000 claims abstract description 78
- 238000012545 processing Methods 0.000 claims abstract description 38
- 238000011156 evaluation Methods 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 42
- 230000009467 reduction Effects 0.000 claims description 19
- 238000012216 screening Methods 0.000 claims description 16
- 238000000605 extraction Methods 0.000 claims description 11
- 230000002159 abnormal effect Effects 0.000 claims description 10
- 238000005520 cutting process Methods 0.000 claims description 10
- 238000007781 pre-processing Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 8
- 238000000513 principal component analysis Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 2
- 238000013210 evaluation model Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000003754 machining Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 102000008297 Nuclear Matrix-Associated Proteins Human genes 0.000 description 2
- 108010035916 Nuclear Matrix-Associated Proteins Proteins 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 210000000299 nuclear matrix Anatomy 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 235000010724 Wisteria floribunda Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000003801 milling Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- 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
- G06N3/088—Non-supervised learning, e.g. competitive learning
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Numerical Control (AREA)
Abstract
本发明属于机械技术领域，公开了一种数控机床刀具磨损预测方法、系统、设备、介质、终端，所述数控机床刀具磨损预测方法包括：进行信号采集与处理；进行信号特征提取；进行信号特征选择；进行健康指标数据集构建；生成健康指标序列数据集；进行刀具磨损预测。本发明通过采集数控机床工作过程中的控制器信号和传感器信号，利用多源信号建立刀具健康指标模型，考虑不同类型的信号反映的刀具磨损情况，提高刀具健康状态评估模型的泛化能力；通过动态时间规整算法构建刀具健康指标，利用不同时刻特征数据之间的差异来衡量刀具的磨损情况，在此基础上进行刀具磨损趋势预测，不需要依赖刀具磨损标签值即可开展刀具磨损预测，提高在实际工业场景中的应用范围与便捷性。
Description
技术领域
本发明属于机械技术领域，尤其涉及一种数控机床刀具磨损预测方法、系统、设备、介质、终端，可用于对数控机床刀具磨损进行预测。
背景技术
目前，作为数控加工中十分重要的部件，刀具磨损带来的问题是数控加工过程中面临的主要难题之一。在数控机床的铣削过程中，刀具的磨损退化是不可避免的。而一旦发生刀具失效，工件的表面质量会达不到要求，从而导致加工效率低下，磨损严重时甚至会造成机床的损坏。因此，有效地预测刀具的磨损情况对提高数控机床的生产效率具有十分重要的意义。目前，采用数据驱动结合机器学习手段是刀具磨损预测技术领域主流的方法与技术，但也存在一定的局限性。
同济大学在其申请的专利文献“基于编码器-解码器阶段注意力机制的刀具磨损预测方法”(专利申请号：202010777546.3，公开号：CN112070208A)中提出了一种刀具磨损预测方法。该方法的步骤是：首先对刀具运行数据进行预处理并提取统计特征，再根据最大信息系数和相关性筛选特征，然后将特征输入到编码器中进行编码，最后将编码特征输入到解码器中，输出刀具的磨损预测值。该方法虽然能实现刀具磨损值的准确预测。但是，该方法仍然存在的不足之处是，模型建立是基于有监督学习的算法，刀具的磨损值是已知的，而实际生产中刀具的磨损值不一定能实时获取，不利于实际生产中的推广使用。
通过上述分析，现有技术存在的问题及缺陷为：现有采用数据驱动结合机器学习进行刀具磨损预测的技术中，模型建立是基于有监督学习的算法，刀具的磨损值是已知的，而实际生产中刀具的磨损值不一定能实时获取，不利于实际生产中的推广使用。
解决以上问题及缺陷的难度为：
如何摆脱无法实时获取刀具磨损值的情况，使得模型能从控制器信号和传感器信号中学习到刀具的健康状态。
解决以上问题及缺陷的意义为：
现有技术针对刀具的磨损预测方法多是基于有监督学习的算法，然而在实际生产中刀具的磨损情况难以实时获取，通过利用不同时刻数据之间的差异来衡量刀具的磨损情况，从而提高模型在实际工业场景中的应用范围与便捷性。本发明旨在提出一种不需要依赖刀具磨损值的刀具磨损预测方法。
发明内容
针对现有技术存在的问题，本发明提供了一种数控机床刀具磨损预测方法、系统、设备、介质、终端，尤其涉及一种数控机床刀具磨损预测技术领域中的基于无监督学习的数控机床刀具磨损预测方法、系统、设备、介质、终端，旨在解决现有技术中利用有监督学习算法建立预测模型的应用场景局限问题。
本发明是这样实现的，一种数控机床刀具磨损预测方法，所述数控机床刀具磨损预测方法包括以下步骤：
步骤一，进行信号采集与处理；
步骤二，进行信号特征提取；
步骤三，进行信号特征选择；
步骤四，进行健康指标数据集构建；
步骤五，生成健康指标序列数据集；
步骤六，进行刀具磨损预测。
进一步，步骤一中，所述信号采集与处理，包括：
(1)采集数控机床工作过程中的信号，即控制器信号和传感器信号，控制器信号主要包括主轴负载和三个方向的机械坐标即x轴机械坐标、y轴机械坐标以及z轴机械坐标，传感器信号主要包括电流信号和三个方向的振动信号即x轴方向振动信号、y轴方向振动信号、z轴方向振动信号以及电流信号；
(2)将采集的信号进行预处理，首先根据三个方向的机械坐标和主轴负载剔除未接触加工对象时采集的信号，然后进行缺失值、异常值处理，利用最小二乘法将信号中的异常趋势项进行去除。
其中，所述利用最小二乘法将信号中的异常趋势项进行去除，包括：
首先利用高阶多项式拟合信号，再根据最小二乘法原理选取合适的多项式系数，使得信号与其拟合信号的误差平方和最小，然后将选取的多项式系数代入高阶多项式中，得到信号的趋势项，最后将信号减去趋势项，得到去除趋势项后的信号。
进一步，步骤二中，所述信号特征提取，包括：
从时域、频域以及时频域三个方面，将步骤一得到的信号进行特征提取。
进一步，步骤三中，所述信号特征选择，包括：
(1)利用单调性和趋势性这两个评价标准，将步骤二得到的信号特征进行筛选，包括：
1)按照下式，计算步骤二得到的每个信号特征单调性值和趋势性值：
其中，Smon为单个信号特征的单调性值，T为单个信号特征的样本长度，dH表示单个特征中每个样本与前一个样本的差值；
其中，Stred为单个信号特征的趋势性值，T为单个信号特征的样本长度，xi为单个信号特征的第i个样本，
为单个信号特征的均值，ti为xi对应的累计工作时间，
为累计工作时间序列的均值；
2)将每个信号特征单调性值和趋势性值的均值作为综合评价值，所有信号特征的综合评价值组成综合评价矩阵；
3)将综合评价矩阵进行归一化，选取综合评价值大于0.5的信号特征组成矩阵，得到筛选后的信号特征。
(2)利用核主成分分析算法对筛选得到的信号特征进行降维处理，包括：
1)将筛选得到的信号特征进行标准化处理；
2)选定核函数，计算核矩阵；
3)将核矩阵进行中心化，得到中心化核矩阵；
4)计算中心化核矩阵的特征值和特征向量，将特征值按降幂排序；
5)设定累计贡献率的阈值为90％，确定主元个数p，选取前p个特征值对应的特征向量组成矩阵，得到降维后的信号特征。
进一步，步骤四中，所述健康指标数据集构建，包括：
将步骤三得到的信号特征进行归一化处理，利用动态时间规整算法构建刀具健康指标数据集，包括：
(1)以一分钟为单位，对降维后的特征数据集进行划分；
(2)以第1个一分钟的信号特征作为刀具的健康基准矩阵，利用动态时间规整算法，计算每个一分钟的特征矩阵与健康基准矩阵的相似度量结果，得到距离向量；
(3)对距离向量进行归一化处理，利用五点三次平滑法对归一化后的距离向量进行平滑处理，得到刀具的健康指标数据集。
进一步，步骤五中，所述生成健康指标序列数据集，包括：
(1)将滑动窗口长度设置为m，窗口的移动步长设置为1，对步骤四得到的健康指标数据集进行滑窗，将每一次滑动的窗口内的所有数据作为的一个序列，将所有序列组成一个二维的序列数据集；其中，m表示任意正整数；
(2)将健康指标数据集中第m+1个数据至最后一个数据组成标签数据集；
(3)将序列数据集和标签数据集组成健康指标序列数据集。
进一步，步骤六中，所述刀具磨损预测，包括：
利用步骤五得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点和失效阈值，进行刀具磨损预测。
本发明的另一目的在于提供一种应用所述数控机床刀具磨损预测方法的数控机床刀具磨损预测系统，所述数控机床刀具磨损预测系统包括：
信号采集与处理模块，用于采集数控机床工作过程中的控制器信号和传感器信号，并将采集的信号进行预处理；
信号特征提取模块，用于从时域、频域以及时频域三个方面，将步信号采集与处理模块得到的信号进行特征提取；
信号特征选择模块，用于利用单调性和趋势性这两个评价标准，将信号特征提取模块得到的信号特征进行筛选，再利用核主成分分析算法对筛选得到的信号特征进行降维处理；
健康指标数据集构建模块，用于将信号特征选择模块得到的信号特征进行归一化处理，利用动态时间规整算法，构建刀具健康指标数据集；
健康指标序列数据集生成模块，用于生成健康指标序列数据集；
刀具磨损预测模块，用于利用健康指标序列数据集生成模块得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点和失效阈值，进行刀具磨损预测。
本发明的另一目的在于提供一种计算机设备，所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如下步骤：
采集数控机床工作过程中的控制器信号和传感器信号，并对信号进行预处理、特征提取以及特征选择，充分挖掘各种信号中与刀具磨损相关的信息，利用动态时间规整算法构建刀具健康指标，利用长短时记忆网络建立刀具磨损预测模型，实现数控机床刀具的磨损预测。
本发明的另一目的在于提供一种信息数据处理终端，所述信息数据处理终端用于实现所述数控机床刀具磨损预测系统。
结合上述的所有技术方案，本发明所具备的优点及积极效果为：本发明提供的基于无监督学习的数控机床刀具磨损预测方法，通过采集数控机床工作过程中的控制器信号和传感器信号，利用动态时间规整算法构建能表征刀具磨损退化状态的健康指标，然后基于长短时记忆网络构建刀具磨损趋势预测模型，有效地克服了无法实时获取磨损值情况下的刀具磨损预测问题，使得本发明提高了刀具磨损预测模型在实际工业场景中的应用范围和便捷性。
本发明通过采集数控机床工作过程中的控制器信号和传感器信号，利用多源信号建立刀具健康指标模型，充分考虑了不同类型的信号反映的刀具磨损情况，使得本发明提高了刀具健康状态评估模型的泛化能力。
本发明通过动态时间规整算法构建刀具健康指标，利用不同时刻特征数据之间的差异来衡量刀具的磨损情况，在此基础上基于长短时记忆网络模型进行刀具磨损趋势预测，使得本发明不需要依赖刀具磨损标签值即可开展刀具磨损预测，提高了在实际工业场景中的应用范围与便捷性。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例中所需要使用的附图做简单的介绍，显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的数控机床刀具磨损预测方法流程图。
图2是本发明实施例提供的数控机床刀具磨损预测方法原理图。
图3是本发明实施例提供的数控机床刀具磨损预测系统结构框图；
图中：1、信号采集与处理模块；2、信号特征提取模块；3、信号特征选择模块；4、指标数据集构建模块；5、指标序列数据集生成模块；6、刀具磨损预测模块。
图4是本发明实施例提供的01号刀具的预测结果图。
图5是本发明实施例提供的02号刀具的预测结果图。
图6是本发明实施例提供的03号刀具的预测结果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
针对现有技术存在的问题，本发明提供了一种数控机床刀具磨损预测方法、系统、设备、介质、终端，下面结合附图对本发明作详细的描述。
如图1所示，本发明实施例提供的数控机床刀具磨损预测方法包括以下步骤：
S101，进行信号采集与处理；
S102，进行信号特征提取；
S103，进行信号特征选择；
S104，进行健康指标数据集构建；
S105，生成健康指标序列数据集；
S106，进行刀具磨损预测。
本发明实施例提供的数控机床刀具磨损预测方法原理图如图2所示。
如图3所示，本发明实施例提供的数控机床刀具磨损预测系统包括：
信号采集与处理模块1，用于采集数控机床工作过程中的控制器信号和传感器信号，并将采集的信号进行预处理；
信号特征提取模块2，用于从时域、频域以及时频域三个方面，将步信号采集与处理模块得到的信号进行特征提取；
信号特征选择模块3，用于利用单调性和趋势性这两个评价标准，将信号特征提取模块得到的信号特征进行筛选，再利用核主成分分析算法对筛选得到的信号特征进行降维处理；
健康指标数据集构建模块4，用于将信号特征选择模块得到的信号特征进行归一化处理，利用动态时间规整算法，构建刀具健康指标数据集；
健康指标序列数据集生成模块5，用于生成健康指标序列数据集；
刀具磨损预测模块6，用于利用指标序列数据集生成模块得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点和失效阈值，进行刀具磨损预测。
下面结合实施例对本发明的技术方案作进一步描述。
实施例1
参照图2，对本发明实现的步骤作进一步的详细说明。
步骤1，信号采集与处理。
第1步，采集数控机床工作过程中的信号，即控制器信号和传感器信号，控制器信号主要包括主轴负载和三个方向的机械坐标即x轴机械坐标、y轴机械坐标以及z轴机械坐标，传感器信号主要包括电流信号和三个方向的振动信号即x轴方向振动信号、y轴方向振动信号、z轴方向振动信号以及电流信号。
第2步，将采集的信号进行预处理，首先根据三个方向的机械坐标和主轴负载剔除未接触加工对象时采集的信号，然后进行缺失值、异常值处理，利用最小二乘法将信号中的异常趋势项进行去除。
所述的利用最小二乘法将信号中的趋势项进行去除的步骤为：首先利用高阶多项式拟合信号，再根据最小二乘法原理选取合适的多项式系数，使得信号与其拟合信号的误差平方和最小，然后将选取的多项式系数代入高阶多项式中，得到信号的趋势项，最后将信号减去趋势项，得到去除趋势项后的信号。
步骤2，信号特征提取。
从时域、频域以及时频域三个方面，将步骤1得到的信号进行特征提取。
步骤3，信号特征选择。
利用单调性和趋势性这两个评价标准，将步骤2得到的信号特征进行筛选。
第1步，按照下式，计算步骤2得到的每个信号特征单调性值和趋势性值。
其中，Smon为单个信号特征的单调性值，T为单个信号特征的样本长度，dH表示单个特征中每个样本与前一个样本的差值。
其中，Stred为单个信号特征的趋势性值，T为单个信号特征的样本长度，xi为单个信号特征的第i个样本，x为单个信号特征的均值，ti为xi对应的累计工作时间，
为累计工作时间序列的均值。
第2步，将每个信号特征单调性值和趋势性值的均值作为综合评价值，所有信号特征的综合评价值组成综合评价矩阵。
第3步，将综合评价矩阵进行归一化，选取综合评价值大于0.5的信号特征组成矩阵，得到筛选后的信号特征。
再利用核主成分分析算法对筛选得到的信号特征进行降维处理。
第1步，将筛选得到的信号特征进行标准化处理。
第2步，选定核函数，计算核矩阵。
第3步，将核矩阵进行中心化，得到中心化核矩阵。
第4步，计算中心化核矩阵的特征值和特征向量，将特征值按降幂排序。
第5步，设定累计贡献率的阈值为90％，确定主元个数p，选取前p个特征值对应的特征向量组成矩阵，得到降维后的信号特征。
步骤4，健康指标数据集构建。
将步骤3得到的信号特征进行归一化处理，利用动态时间规整算法，构建刀具健康指标数据集。
所述的利用动态时间规整算法构建刀具健康指标数据集的步骤如下：
第1步，以一分钟为单位，对降维后的特征数据集进行划分。
第2步，以第1个一分钟的信号特征作为刀具的健康基准矩阵，利用动态时间规整算法，计算每个一分钟的特征矩阵与健康基准矩阵的相似度量结果，得到距离向量。
第3步，对距离向量进行归一化处理，利用五点三次平滑法对归一化后的距离向量进行平滑处理，得到刀具的健康指标数据集。
步骤5，生成健康指标序列数据集。
第1步，将滑动窗口长度设置为m，窗口的移动步长设置为1，对步骤(4)得到的健康指标数据集进行滑窗，将每一次滑动的窗口内的所有数据作为的一个序列，将所有序列组成一个二维的序列数据集，其中，m表示任意正整数。
第2步，将健康指标数据集中第m+1个数据至最后一个数据组成标签数据集。
第3步，将序列数据集和标签数据集组成健康指标序列数据集。
步骤6，刀具磨损预测。
利用步骤5得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点为70％，失效阈值为1，进行刀具磨损预测。
实施例2
本发明实施例使用的数据集是第二届工业大数据创新竞赛中的富士康机加刀具数据，数据采集自数控机床真实加工过程中，一把全新的刀具从开始正常加工直至刀具寿命终止时停止采集。
1)信号采集与处理。
1.1)采集数控机床工作过程中的信号，即控制器信号和传感器信号，控制器信号主要包括主轴负载和三个方向的机械坐标即x轴机械坐标、y轴机械坐标以及z轴机械坐标，传感器信号主要包括电流信号和三个方向的振动信号即x轴方向振动信号、y轴方向振动信号、z轴方向振动信号以及电流信号。
总共采集了3把刀具的信号，分别为01号刀具、02号刀具以及03号刀具，其中，控制器信号的采样频率为33Hz，传感器信号的采样频率为25600Hz。
1.2)将采集的信号进行预处理，首先根据三个方向的机械坐标和主轴负载剔除未接触加工对象时采集的信号，然后进行缺失值、异常值处理，利用最小二乘法将信号中的异常趋势项进行去除。
所述的利用最小二乘法将信号中的趋势项进行去除的步骤为：首先利用高阶多项式拟合信号，再根据最小二乘法原理选取合适的多项式系数，使得信号与其拟合信号的误差平方和最小，然后将选取的多项式系数代入高阶多项式中，得到信号的趋势项，最后将信号减去趋势项，得到去除趋势项后的信号。
2)信号特征提取。
从时域、频域以及时频域三个方面，将步骤1)得到的信号进行特征提取。
3)信号特征选择。
3.1)利用单调性和趋势性这两个评价标准，将步骤2)得到的信号特征进行筛选。
3.1.1)按照下式，计算步骤2)得到的每个信号特征单调性值和趋势性值。
其中，Smon为单个信号特征的单调性值，T为单个信号特征的样本长度，dH表示单个特征中每个样本与前一个样本的差值。
其中，Stred为单个信号特征的趋势性值，T为单个信号特征的样本长度，xi为单个信号特征的第i个样本，
为单个信号特征的均值，ti为xi对应的累计工作时间，
为累计工作时间序列的均值。
3.1.2)将每个信号特征单调性值和趋势性值的均值作为综合评价值，所有信号特征的综合评价值组成综合评价矩阵。
3.1.3)将综合评价矩阵进行归一化，选取综合评价值大于0.5的信号特征组成矩阵，得到筛选后的信号特征。
3.2)利用核主成分分析算法对筛选得到的信号特征进行降维处理。
3.2.1)将筛选得到的信号特征进行标准化处理。
3.2.2)选定核函数，计算核矩阵。
3.2.3)将核矩阵进行中心化，得到中心化核矩阵。
3.2.4)计算中心化核矩阵的特征值和特征向量，将特征值按降幂排序。
3.2.5)设定累计贡献率的阈值为90％，确定主元个数p，选取前p个特征值对应的特征向量组成矩阵，得到降维后的信号特征。
4)健康指标数据集构建。
将步骤3)得到的信号特征进行归一化处理，利用动态时间规整算法，构建刀具健康指标数据集。
所述的利用动态时间规整算法构建刀具健康指标数据集的步骤如下：
第1步，以一分钟为单位，对降维后的特征数据集进行划分。
第2步，以第1个一分钟的信号特征作为刀具的健康基准矩阵，利用动态时间规整算法，计算每个一分钟的特征矩阵与健康基准矩阵的相似度量结果，得到距离向量。
第3步，对距离向量进行归一化处理，利用五点三次平滑法对归一化后的距离向量进行平滑处理，得到刀具的健康指标数据集。
5)生成健康指标序列数据集。
5.1)将滑动窗口长度设置为m，窗口的移动步长设置为1，对步骤(4)得到的健康指标数据集进行滑窗，将每一次滑动的窗口内的所有数据作为的一个序列，将所有序列组成一个二维的序列数据集，其中，m表示任意正整数。
5.2)将健康指标数据集中第m+1个数据至最后一个数据组成标签数据集。
5.3)将序列数据集和标签数据集组成健康指标序列数据集。
6)刀具磨损预测。
利用步骤5)得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点为70％，失效阈值为1，进行刀具磨损预测。各刀具的预测结果如图4、图5和图6所示。
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现，所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘SolidState Disk(SSD))等。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种数控机床刀具磨损预测方法，其特征在于，所述数控机床刀具磨损预测方法包括以下步骤：
步骤一，进行信号采集与处理；
步骤二，进行信号特征提取；
步骤三，进行信号特征选择；
步骤四，进行健康指标数据集构建；
步骤五，生成健康指标序列数据集；
步骤六，进行刀具磨损预测。
2.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤一中，所述信号采集与处理，包括：
(1)采集数控机床工作过程中的信号，即控制器信号和传感器信号，控制器信号主要包括主轴负载和三个方向的机械坐标即x轴机械坐标、y轴机械坐标以及z轴机械坐标，传感器信号主要包括电流信号和三个方向的振动信号即x轴方向振动信号、y轴方向振动信号、z轴方向振动信号以及电流信号；
(2)将采集的信号进行预处理，首先根据三个方向的机械坐标和主轴负载剔除未接触加工对象时采集的信号，然后进行缺失值、异常值处理，利用最小二乘法将信号中的异常趋势项进行去除；
其中，所述利用最小二乘法将信号中的异常趋势项进行去除，包括：
首先利用高阶多项式拟合信号，再根据最小二乘法原理选取合适的多项式系数，使得信号与其拟合信号的误差平方和最小，然后将选取的多项式系数代入高阶多项式中，得到信号的趋势项，最后将信号减去趋势项，得到去除趋势项后的信号。
3.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤二中，所述信号特征提取，包括：
从时域、频域以及时频域三个方面，将步骤一得到的信号进行特征提取。
4.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤三中，所述信号特征选择，包括：
(1)利用单调性和趋势性这两个评价标准，将步骤二得到的信号特征进行筛选，包括：
1)按照下式，计算步骤二得到的每个信号特征单调性值和趋势性值：
其中，Smon为单个信号特征的单调性值，T为单个信号特征的样本长度，dH表示单个特征中每个样本与前一个样本的差值；
其中，Stred为单个信号特征的趋势性值，T为单个信号特征的样本长度，xi为单个信号特征的第i个样本，
为单个信号特征的均值，ti为xi对应的累计工作时间，
为累计工作时间序列的均值；
2)将每个信号特征单调性值和趋势性值的均值作为综合评价值，所有信号特征的综合评价值组成综合评价矩阵；
3)将综合评价矩阵进行归一化，选取综合评价值大于0.5的信号特征组成矩阵，得到筛选后的信号特征；
(2)利用核主成分分析算法对筛选得到的信号特征进行降维处理，包括：
1)将筛选得到的信号特征进行标准化处理；
2)选定核函数，计算核矩阵；
3)将核矩阵进行中心化，得到中心化核矩阵；
4)计算中心化核矩阵的特征值和特征向量，将特征值按降幂排序；
5)设定累计贡献率的阈值为90％，确定主元个数p，选取前p个特征值对应的特征向量组成矩阵，得到降维后的信号特征。
5.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤四中，所述健康指标数据集构建，包括：
将步骤三得到的信号特征进行归一化处理，利用动态时间规整算法构建刀具健康指标数据集，包括：
(1)以一分钟为单位，对降维后的特征数据集进行划分；
(2)以第1个一分钟的信号特征作为刀具的健康基准矩阵，利用动态时间规整算法，计算每个一分钟的特征矩阵与健康基准矩阵的相似度量结果，得到距离向量；
(3)对距离向量进行归一化处理，利用五点三次平滑法对归一化后的距离向量进行平滑处理，得到刀具的健康指标数据集。
6.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤五中，所述生成健康指标序列数据集，包括：
(1)将滑动窗口长度设置为m，窗口的移动步长设置为1，对步骤四得到的健康指标数据集进行滑窗，将每一次滑动的窗口内的所有数据作为的一个序列，将所有序列组成一个二维的序列数据集；其中，m表示任意正整数；
(2)将健康指标数据集中第m+1个数据至最后一个数据组成标签数据集；
(3)将序列数据集和标签数据集组成健康指标序列数据集。
7.如权利要求1所述数控机床刀具磨损预测方法，其特征在于，步骤六中，所述刀具磨损预测，包括：利用步骤五得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点和失效阈值，进行刀具磨损预测。
8.一种实施权利要求1～7任意一项所述数控机床刀具磨损预测方法的数控机床刀具磨损预测系统，其特征在于，所述数控机床刀具磨损预测系统包括：
信号采集与处理模块，用于采集数控机床工作过程中的控制器信号和传感器信号，并将采集的信号进行预处理；
信号特征提取模块，用于从时域、频域以及时频域三个方面，将步信号采集与处理模块得到的信号进行特征提取；
信号特征选择模块，用于利用单调性和趋势性这两个评价标准，将信号特征提取模块得到的信号特征进行筛选，再利用核主成分分析算法对筛选得到的信号特征进行降维处理；
健康指标数据集构建模块，用于将信号特征选择模块得到的信号特征进行归一化处理，利用动态时间规整算法，构建刀具健康指标数据集；
健康指标序列数据集生成模块，用于生成健康指标序列数据集；
刀具磨损预测模块，用于利用指标序列数据集生成模块得到的健康指标序列数据集，构造训练集，建立长短时记忆网络预测模型，设定预测起点和失效阈值，进行刀具磨损预测。
9.一种计算机设备，其特征在于，所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如下步骤：
采集数控机床工作过程中的控制器信号和传感器信号，并对信号进行预处理、特征提取以及特征选择，充分挖掘各种信号中与刀具磨损相关的信息，利用动态时间规整算法构建刀具健康指标，利用长短时记忆网络建立刀具磨损预测模型，实现数控机床刀具的磨损预测。
10.一种信息数据处理终端，其特征在于，所述信息数据处理终端用于实现如权利要求8所述数控机床刀具磨损预测系统。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110642404.0A CN113569903B (zh) | 2021-06-09 | 2021-06-09 | 数控机床刀具磨损预测方法、系统、设备、介质、终端 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110642404.0A CN113569903B (zh) | 2021-06-09 | 2021-06-09 | 数控机床刀具磨损预测方法、系统、设备、介质、终端 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113569903A true CN113569903A (zh) | 2021-10-29 |
| CN113569903B CN113569903B (zh) | 2024-04-09 |
Family
ID=78161829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110642404.0A Active CN113569903B (zh) | 2021-06-09 | 2021-06-09 | 数控机床刀具磨损预测方法、系统、设备、介质、终端 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113569903B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114091544A (zh) * | 2021-11-29 | 2022-02-25 | 西南交通大学 | 基于振动信号与神经网络的tbm盘形滚刀磨损识别系统 |
| CN114611633A (zh) * | 2022-05-10 | 2022-06-10 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种电磁阀的健康监测方法 |
| CN114952422A (zh) * | 2022-06-07 | 2022-08-30 | 中科航迈数控软件(深圳)有限公司 | 一种数控机床加工刀具状态实时预测方法 |
| CN115034504A (zh) * | 2022-06-28 | 2022-09-09 | 重庆大学 | 基于云边协同训练的刀具磨损状态预测系统及方法 |
| CN115122154A (zh) * | 2022-08-30 | 2022-09-30 | 南通市海门区海纳机械有限公司 | 一种金属构件加工数控铣床运行检测及控制方法 |
| CN115415851A (zh) * | 2022-10-08 | 2022-12-02 | 清华大学 | 一种基于函数型数据主成分分析的刀具健康监测方法 |
| CN115922443A (zh) * | 2023-01-09 | 2023-04-07 | 湖南大学 | 一种基于滑动窗口的刀具剩余使用寿命实时预测方法 |
| CN116070527A (zh) * | 2023-03-07 | 2023-05-05 | 南京航空航天大学 | 基于退化模型的铣削刀具剩余寿命预测方法 |
| CN117235678A (zh) * | 2023-11-14 | 2023-12-15 | 山东交通学院 | 一种木材切割过程中刀片磨损程度的检测方法 |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103674511A (zh) * | 2013-03-18 | 2014-03-26 | 北京航空航天大学 | 一种基于emd-svd与mts的机械磨损件性能评估与预测方法 |
| US20160189045A1 (en) * | 2014-12-31 | 2016-06-30 | Facebook, Inc. | Predicting computer model accuracy |
| CN107094093A (zh) * | 2017-04-13 | 2017-08-25 | 北京网瑞达科技有限公司 | 网络系统中的网络设备周期性的性能指标数据预测方法 |
| CN107138540A (zh) * | 2017-04-06 | 2017-09-08 | 首钢总公司 | 一种带钢断面板廓形状的拟合方法及评价方法 |
| JP2018086712A (ja) * | 2016-11-30 | 2018-06-07 | 株式会社日立製作所 | 工具摩耗予測装置およびその方法 |
| US20180272491A1 (en) * | 2017-03-24 | 2018-09-27 | National Cheng Kung University | Tool wear monitoring and predicting method |
| CN109141881A (zh) * | 2018-07-06 | 2019-01-04 | 东南大学 | 一种深度自编码网络的旋转机械健康评估方法 |
| CN110119551A (zh) * | 2019-04-29 | 2019-08-13 | 西安电子科技大学 | 基于机器学习的盾构机刀具磨损退化关联特征分析方法 |
| CN110355608A (zh) * | 2019-07-18 | 2019-10-22 | 浙江大学 | 基于自注意力机制与深度学习的刀具磨损量预测方法 |
| CN110533065A (zh) * | 2019-07-18 | 2019-12-03 | 西安电子科技大学 | 基于自编码特征和深度学习回归模型的盾构姿态预测方法 |
| CN110561191A (zh) * | 2019-07-30 | 2019-12-13 | 西安电子科技大学 | 基于pca与自编码器的数控机床刀具磨损数据处理方法 |
| CN111126255A (zh) * | 2019-12-23 | 2020-05-08 | 华中科技大学 | 基于深度学习回归算法的数控机床刀具磨损值预测方法 |
| US20200209831A1 (en) * | 2019-01-02 | 2020-07-02 | Industrial Technology Research Institute | Tool-life prediction system and method thereof |
| JP2020138254A (ja) * | 2019-02-27 | 2020-09-03 | ファナック株式会社 | 工作機械の加工不良発生予測システム |
| CN111687689A (zh) * | 2020-06-23 | 2020-09-22 | 重庆大学 | 一种基于lstm和cnn的刀具磨损状态预测方法及装置 |
| KR102249849B1 (ko) * | 2020-12-31 | 2021-05-10 | 주식회사 한국가스기술공사 | 기계적 결함에 대한 예측 모델을 처리하는 방법 |
-
2021
- 2021-06-09 CN CN202110642404.0A patent/CN113569903B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103674511A (zh) * | 2013-03-18 | 2014-03-26 | 北京航空航天大学 | 一种基于emd-svd与mts的机械磨损件性能评估与预测方法 |
| US20160189045A1 (en) * | 2014-12-31 | 2016-06-30 | Facebook, Inc. | Predicting computer model accuracy |
| JP2018086712A (ja) * | 2016-11-30 | 2018-06-07 | 株式会社日立製作所 | 工具摩耗予測装置およびその方法 |
| US20180272491A1 (en) * | 2017-03-24 | 2018-09-27 | National Cheng Kung University | Tool wear monitoring and predicting method |
| CN107138540A (zh) * | 2017-04-06 | 2017-09-08 | 首钢总公司 | 一种带钢断面板廓形状的拟合方法及评价方法 |
| CN107094093A (zh) * | 2017-04-13 | 2017-08-25 | 北京网瑞达科技有限公司 | 网络系统中的网络设备周期性的性能指标数据预测方法 |
| CN109141881A (zh) * | 2018-07-06 | 2019-01-04 | 东南大学 | 一种深度自编码网络的旋转机械健康评估方法 |
| US20200209831A1 (en) * | 2019-01-02 | 2020-07-02 | Industrial Technology Research Institute | Tool-life prediction system and method thereof |
| JP2020138254A (ja) * | 2019-02-27 | 2020-09-03 | ファナック株式会社 | 工作機械の加工不良発生予測システム |
| CN110119551A (zh) * | 2019-04-29 | 2019-08-13 | 西安电子科技大学 | 基于机器学习的盾构机刀具磨损退化关联特征分析方法 |
| CN110355608A (zh) * | 2019-07-18 | 2019-10-22 | 浙江大学 | 基于自注意力机制与深度学习的刀具磨损量预测方法 |
| CN110533065A (zh) * | 2019-07-18 | 2019-12-03 | 西安电子科技大学 | 基于自编码特征和深度学习回归模型的盾构姿态预测方法 |
| CN110561191A (zh) * | 2019-07-30 | 2019-12-13 | 西安电子科技大学 | 基于pca与自编码器的数控机床刀具磨损数据处理方法 |
| CN111126255A (zh) * | 2019-12-23 | 2020-05-08 | 华中科技大学 | 基于深度学习回归算法的数控机床刀具磨损值预测方法 |
| CN111687689A (zh) * | 2020-06-23 | 2020-09-22 | 重庆大学 | 一种基于lstm和cnn的刀具磨损状态预测方法及装置 |
| KR102249849B1 (ko) * | 2020-12-31 | 2021-05-10 | 주식회사 한국가스기술공사 | 기계적 결함에 대한 예측 모델을 처리하는 방법 |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114091544A (zh) * | 2021-11-29 | 2022-02-25 | 西南交通大学 | 基于振动信号与神经网络的tbm盘形滚刀磨损识别系统 |
| CN114611633A (zh) * | 2022-05-10 | 2022-06-10 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种电磁阀的健康监测方法 |
| CN114952422A (zh) * | 2022-06-07 | 2022-08-30 | 中科航迈数控软件(深圳)有限公司 | 一种数控机床加工刀具状态实时预测方法 |
| CN114952422B (zh) * | 2022-06-07 | 2023-10-27 | 中科航迈数控软件(深圳)有限公司 | 一种数控机床加工刀具状态实时预测方法 |
| CN115034504A (zh) * | 2022-06-28 | 2022-09-09 | 重庆大学 | 基于云边协同训练的刀具磨损状态预测系统及方法 |
| CN115034504B (zh) * | 2022-06-28 | 2024-05-28 | 重庆大学 | 基于云边协同训练的刀具磨损状态预测系统及方法 |
| CN115122154B (zh) * | 2022-08-30 | 2023-08-18 | 广东昭明电子集团股份有限公司 | 一种金属构件加工数控铣床运行检测及控制方法 |
| CN115122154A (zh) * | 2022-08-30 | 2022-09-30 | 南通市海门区海纳机械有限公司 | 一种金属构件加工数控铣床运行检测及控制方法 |
| CN115415851B (zh) * | 2022-10-08 | 2023-09-26 | 清华大学 | 一种基于函数型数据主成分分析的刀具健康监测方法 |
| CN115415851A (zh) * | 2022-10-08 | 2022-12-02 | 清华大学 | 一种基于函数型数据主成分分析的刀具健康监测方法 |
| CN115922443A (zh) * | 2023-01-09 | 2023-04-07 | 湖南大学 | 一种基于滑动窗口的刀具剩余使用寿命实时预测方法 |
| CN116070527A (zh) * | 2023-03-07 | 2023-05-05 | 南京航空航天大学 | 基于退化模型的铣削刀具剩余寿命预测方法 |
| CN117235678A (zh) * | 2023-11-14 | 2023-12-15 | 山东交通学院 | 一种木材切割过程中刀片磨损程度的检测方法 |
| CN117235678B (zh) * | 2023-11-14 | 2024-01-26 | 山东交通学院 | 一种木材切割过程中刀片磨损程度的检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113569903B (zh) | 2024-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113569903A (zh) | 数控机床刀具磨损预测方法、系统、设备、介质、终端 | |
| CN113560955B (zh) | 一种数控机床刀具剩余使用寿命预测方法、系统及应用 | |
| CN111459778B (zh) | 运维系统异常指标检测模型优化方法、装置及存储介质 | |
| CN116300691B (zh) | 一种多轴联动数控加工的状态监测方法及系统 | |
| CN112766342A (zh) | 一种电气设备的异常检测方法 | |
| CN111538311B (zh) | 一种基于数据挖掘的机械设备柔性多状态自适应预警方法及装置 | |
| CN112070208B (zh) | 基于编码器-解码器阶段注意力机制的刀具磨损预测方法 | |
| CN113579851B (zh) | 一种基于自适应分段pca的非平稳钻削过程监测方法 | |
| US11644812B2 (en) | Machine tool management method, machine tool management system and medium | |
| CN113798920A (zh) | 一种基于变分自动编码器与极限学习机的刀具磨损状态监测方法 | |
| US20220179393A1 (en) | Machine tool evaluation method, machine tool evaluation system and medium | |
| Dong | Combining unsupervised and supervised learning for asset class failure prediction in power systems | |
| CN114326593B (zh) | 一种刀具寿命预测系统及方法 | |
| CN113505531B (zh) | 基于结合传统诊断与机器学习方法且针对往复机械气门间隙异常故障的诊断模型构建方法 | |
| CN112016800B (zh) | 一种基于有效性指标的特征选择方法与系统 | |
| CN116383645A (zh) | 一种基于异常检测的系统健康度智能监测评估方法 | |
| CN116126807A (zh) | 一种日志分析方法及相关装置 | |
| CN115081551A (zh) | 基于K-Means聚类和优化RVM线损模型建立方法及系统 | |
| CN113535527A (zh) | 一种面向实时流数据预测性分析的降载方法及系统 | |
| CN111639664A (zh) | 一种基于多策略模式的线条痕迹批量比对系统 | |
| CN112632466A (zh) | 一种基于主成分分析和深度双向长短时记忆网络的轴承故障预测方法 | |
| CN116861204B (zh) | 基于数字孪生的智能制造设备数据管理系统 | |
| Di et al. | Tool Fault Diagnosis Method Based on Multiscale-Channel Attention Network | |
| CN116170200A (zh) | 电力监控系统时间序列异常检测方法、系统、设备及存储介质 | |
| CN117131425A (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 |