CN113496312A - 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质 - Google Patents

基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质 Download PDF

Info

Publication number
CN113496312A
CN113496312A CN202110755673.8A CN202110755673A CN113496312A CN 113496312 A CN113496312 A CN 113496312A CN 202110755673 A CN202110755673 A CN 202110755673A CN 113496312 A CN113496312 A CN 113496312A
Authority
CN
China
Prior art keywords
data
dimensional
cutter
value
network model
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.)
Pending
Application number
CN202110755673.8A
Other languages
English (en)
Inventor
袁东风
庞蓓
李东阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN202110755673.8A priority Critical patent/CN113496312A/zh
Publication of CN113496312A publication Critical patent/CN113496312A/zh
Pending legal-status Critical Current

Links

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • 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/044—Recurrent networks, e.g. Hopfield networks
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Human Resources & Organizations (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • General Business, Economics & Management (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质，包括：(1)采集多维刀具监测数据；(2)数据预处理；(3)特征提取；(4)构造训练集、测试集；(5)搭建基于LSTM算法网络模型；(6)训练基于LSTM算法网络模型；(7)通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测。本发明将长短期记忆网络运用到工业场景，对数控机床刀具的剩余寿命进行预测。本发明在不同特征维度和尺度提取了不同传感器采集到的多通道刀具监测数据的有效信号，充分反映了刀具的状态信息。

Description

基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预 测方法、设备及存储介质
技术领域
本发明涉及智能制造刀具故障诊断与健康管理领域，特别涉及一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质。
背景技术
随着制造业与新一代信息技术的深度融合，生产加工流程变得更加智能化、自动化，设备的故障诊断及预测性维护成为保障企业生产效率，降低生产成本的关键问题。刀具作为数控机床的“牙齿”，其实时的运行状态直接影响着机床加工效率和产品质量。精准预测刀具剩余寿命能够有效避免由于刀具状态异常而导致的工件质量问题，从而提高刀具的使用效率。
对生产过程中的刀具状态进行监测的研究最早可追溯到上世纪，甚至早于数控加工技术的大量应用。随着切削工艺的不断完善、自动化领域的不断发展、计算机技术的普及，刀具状态监测技术也在不断改进。在工业制造中，很多因素都影响着刀具的剩余寿命，例如：数控机床设置的切削参数、要加工的工件所使用的毛坯材料、用于切削的刀具的材料、切削液的温度等等，这些因素都或多或少地影响着刀具的剩余寿命，传统的基于累计加工时间预测和基于模型预测的刀具剩余寿命预测方法都过于受到实际加工工况的制约，导致预测的精度不高。因此，针对不同的工况和生产场景及时获取刀具状态变得尤为关键。
刀具剩余寿命预测与刀具状态监测密不可分，从监测方式来看主要分为离线与在线两种方式。离线监测的数据顾名思义来自切削的中止时刻，根据此时得到的刀具各项参数数据评估其剩余工作时长。但由于离线监测的必要条件是机床停机，在很大程度上妨碍了生产加工流程，并未被广泛使用。而在线监测则是在切削过程中不断地获取刀具各项参数，从而对其及时调整，也有利于更精准地把握刀具剩余寿命，因此，在线监测是目前主流的技术手段。
在线监测的最大优势就是在于其获取数据的实时性，因此，传感器在其中的作用不可忽视。
如今机械行业正与计算机技术深度交叉融合，当下热门的人工智能关键性技术深度学习也开始应用于刀具的寿命预测。Wen-Yang Chang等研究了反向传播神经网络算法的迭代梯度收敛性，该模型很好地预测了刀具刀片的磨损情况。Denis Boing等基于普通最小二乘法计算刀具磨损率，以刀具去除的材料体积作为刀具寿命结束的判据，估算刀具寿命所建立的刀具寿命模型能够在较高的切削速度下预测误差在4％以下的刀具寿命。黄鹤翔提出了一种基于3-KMMBS的刀具磨损检测方法，并搭建了更准确、更稳定、更有效的基于深度学习的刀具磨损检测系统。但是现有的这些工作忽略了多维特征融合以提高剩余工具使用寿命预测的准确性。由于在实际的工业生产中，描述刀具状况的数据并不是单一的信号，多传感器测量到的多通道数据才能尽可能完整地描述刀具的状态。刀具剩余寿命预测的关键就在于不同传感器测得的多维信号。传感器把切削力、振动信号等物理信号转换成电信号，这些数据就反映了刀具的实时状态。因此，使用从数据出发，根据这些传感器测得的多维信号进行分析才是研究的关键。
发明内容
针对刀具监测数据维度高难以有效处理的问题，本发明的目的在于提出一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法。将提取到的刀具监测信号特征输入到长短期记忆网络中，利用这些有效信号特征预测刀具剩余寿命。
本发明采用的刀具数据集的每份样本都有七个通道的数据，采用时频域分析和小波变换等刀具监测信号特征提取方法，充分挖掘监测数据中刀具磨损相关信息，结合特征筛选方法提取有效信号特征，尽可能完整地描述了刀具的状态。基于特征融合思路，将提取到的特征进行拼接，构成多维特征矩阵。
本发明采用长短期记忆网络(Long-Short Term Memory，LSTM)完成刀具磨损量预测任务，LSTM是一种可以学习长期依赖信息的特殊的循环神经网络，通过刻意设计避免了长期依赖问题。
本发明还提供了一种计算机设备及存储介质。
术语解释：
1、Adam(Adaptive Moment Estimation)算法利用梯度的一阶矩估计和二阶矩估计动态调整每个参数的学习率，能够在经过偏置校正后使每一次迭代学习率都有一个确定的范围，使得参数比较平稳。
2、LSTM，Long-Short Term Memory，长短期记忆网络，由Hochreiter&Schmidhuber提出，解决了循环神经网络在学习近期信息时，由于网络经过多层传播之后发生的梯度消失或梯度爆炸问题而导致网络参数更新缓慢，无法完成长期记忆的问题。
3、绝对均值，对轴承冲击振动为经绝对值处理后的平均值，检测值比峰值稳定。
4、峰值，反映某时刻振幅的最大值，适用于具有顺时冲击的故障诊断。
5、均方根值，对时间平均，用于反映信号的能量大小，适用于振幅值随时间缓慢变化的故障诊断。
6、方根幅值，算术平方根的平均值的平方。
7、歪度值，一般用来判断磨损情况。
8、峭度值，一般用来测量局部缺陷。
9、波形因子，均方根值与绝对均值的比值，即脉冲因子与峰值因子的比值。
10、脉冲因子，信号峰值与绝对均值的比值，用于检测信号中是否存在冲击。
11、歪度因子，三阶中心矩和标准差的三次方的比值。
12、峰值因子，用于诊断局部剥落、擦伤、刻痕、凹痕等离散型的缺陷，这类缺陷产生的脉冲波形总能量不大，但是波形的尖峰度明显，是信号峰值与均方根值的比值，用来检测信号中是否存在冲击的统计指标。
13、裕度因子，信号峰值与方根幅值的比值，用于检测机械设备的磨损情况。
14、峭度因子，表示波形平缓程度的，用于描述变量的分布。正态分布的峭度等于3，峭度小于3时分布的曲线会较“平”，大于3时分布的曲线较“陡”。
15、重心频率，表示功率谱曲线下面积重心。
16、均方频率，是与随机过程有关的一个量。
17、均方根频率，均方频率经过开方后的值。
18、频率方差，受功率谱和重心频率影响的统计量。
本发明的技术方案为：
一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，包括步骤如下：
(1)获取多维刀具监测数据；
多维刀具监测数据包括测有刀具的磨损值的样本数据及未测有刀具的磨损值的样本数据；
每个样本数据均包括X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号，X维力、Y维力、Z维力是指在切削过程中x轴、y轴、z轴这三个方向上刀具受到的切削力，X维振动、Y维振动、Z维振动是指在切削过程中x轴、y轴、z轴这三个方向上的振动信号(加速度)，x轴、y轴、z轴所处的坐标系为机床主轴坐标系，是根据右手笛卡尔直角坐标系建立的，主轴方向为z轴，垂直于主轴方向的平面上两个相互垂直的方向为x轴、y轴；声发射信号是金属在加工过程中由于内部分子晶格发生畸变、金属裂纹加剧以及其在塑性形变时所释放出来的一种超高频应力波脉冲信号；
(2)数据预处理；
(3)特征提取与融合；
(4)构造训练集、测试集；
(5)搭建基于LSTM算法网络模型；
(6)训练基于LSTM算法网络模型；
(7)通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测。
根据本发明优选的，步骤(2)中，对多维刀具监测数据数据预处理的具体实现过程包括：
读取测有刀具的磨损值的样本数据中的X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号这七通道数据并将其转换为DataFrame数据帧，转换后将这七个字段分别命名为Fx、Fy、Fz、Ax、Ay、Az、AE-rms；方便后续的特征提取工作。
将测有刀具的磨损值的样本数据中的刀具的磨损值在x轴、y轴、z轴这三个方向上的磨损量，求均值之后作为标签，得到若干标签矩阵；
将若干标签矩阵以.npy文件形式存储。
根据本发明优选的，步骤(3)中，对步骤(2)数据预处理后的数据进行特征提取的具体实现过程包括：
读取步骤(2)数据预处理后的数据帧，
在时域上，计算绝对均值、峰值、均方根值、方根幅值、歪度值、峭度值、波形因子、脉冲因子、歪度因子、峰值因子、裕度因子、峭度因子12个特征；
绝对均值(Absolute mean)的计算公式为：
Figure BDA0003147223200000041
峰值(Max)的计算公式为：max(zi)；
均方根值(Root mean square)的计算公式为：
Figure BDA0003147223200000042
方根幅值(Square root amplitude)的计算公式为：
Figure BDA0003147223200000043
歪度值(Skewness)的计算公式为：
Figure BDA0003147223200000044
峭度值(Kurtosis)的计算公式为：
Figure BDA0003147223200000045
波形因子(Shape factor)的计算公式为：
Figure BDA0003147223200000046
脉冲因子(Pulse factor)的计算公式为：
Figure BDA0003147223200000047
歪度因子(Skewness factor)的计算公式为：
Figure BDA0003147223200000051
峰值因子(Crest factor)的计算公式为：
Figure BDA0003147223200000052
裕度因子(Clearance factor)的计算公式为
Figure BDA0003147223200000053
峭度因子(Kurtosis factor)的计算公式为：
Figure BDA0003147223200000054
每一条数据长度为n，zij表示第i条数据第j个数据点，zi表示第i条数据，S(f)表示功率谱，max表示峰值；
在频域上，计算重心频率、均方频率、均方根频率、频率方差4个特征；
重心频率(FC)的计算公式为：
Figure BDA0003147223200000055
均方频率(MSF)的计算公式为：
Figure BDA0003147223200000056
均方根频率(RMSF)的计算公式为：
Figure BDA0003147223200000057
频率方差(VF)的计算公式为：
Figure BDA0003147223200000058
FC表示重心频率，S(f)为功率谱，f表示对输入的数据进行快速傅里叶变换之后得到的对应频率点，采样频率为1/50000；
在小波域上，用db3小波变换提取8个特征。
经过特征提取后，共得到24维特征，将这24维特征融合为特征矩阵，这样可以很好地反映刀具的磨损状态。均以.npy文件形式存储。
根据本发明优选的，步骤(4)中，构造训练集、测试集的具体实现过程包括：
读取标签矩阵和特征矩阵，对其进行数据标准差标准化处理后，得到符合标准正态分布即均值为0、方差为1的数据；
把每个标签矩阵(因变量)和对应的特征矩阵(自变量)相匹配作为样本数据集；
将样本数据集的部分数据作为训练集，剩余部分作为测试集。
进一步优选的，将样本数据集的80％作为训练集，剩余部分作为测试集。
根据本发明优选的，步骤(5)中，搭建基于LSTM算法网络模型，具体是指：
基于LSTM算法网络模型为长短期记忆网络，长短期记忆网络包括输入层、第一隐藏层、第二隐藏层、全连接层及输出层；
输入层的输入特征维数为6，第一隐藏层、第二隐藏层的节点个数均为64；
将训练集输入第一隐藏层，第二隐藏层的输入为第一隐藏层的计算结果，在第二隐藏层设置dropout率为0.05的dropout层，避免过拟合问题，全连接层输入64个神经元，输出10个神经元，输出层采用线性激活函数，完成第二隐藏层到输出层的线性变换，使输出层输出刀具的磨损值。
根据本发明优选的，步骤(6)中，训练基于LSTM算法网络模型，具体是指：
将训练集输入基于LSTM算法网络模型，将标签作为基于LSTM算法网络模型的输出，对基于LSTM算法网络模型进行训练；训练的同时记录每一次训练周期的损失函数；
训练参数为：Epoch＝500，BatchSize＝128，Learningrate＝0.0001；
将验证集输入长短期记忆网络模型进行验证，在基于LSTM算法网络模型训练过程中用Adam算法进行优化，优化更新学习率Learningrate，验证时记录每一次训练周期的损失函数，得到训练好的基于LSTM算法网络模型。
根据本发明优选的，步骤(7)中，通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测，具体是指：
将测试集输入训练好的基于LSTM算法网络模型，得到预测磨损值，画出刀具预测磨损值曲线，刀具预测磨损值曲线中，横坐标为切削次数，纵坐标为磨损值。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。
本发明的有益效果为：
1.本发明将长短期记忆网络运用到工业场景，对数控机床刀具的剩余寿命进行预测。
2.本发明在不同特征维度和尺度提取了不同传感器采集到的多通道刀具监测数据的有效信号，充分反映了刀具的状态信息。
3.本发明采用长短期记忆网络，能够充分利用刀具监测数据的多维特征，精准预测刀具剩余寿命，指导设备维护计划更新与优化，从而有效提高企业生产效率、降低生产成本、保证产品质量。
4.本发明搭建的网络利用了两层长短期记忆网络，第二层的长短期记忆网络接收第一层长短期记忆网络计算的结果，提升模型性能。
附图说明
图1是本发明基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的流程示意图；
图2是本发明基于LSTM算法网络模型的结构示意图；
图3(a)是基于LSTM算法网络模型训练集训练过程的损失函数曲线图；
图3(b)是基于LSTM算法网络模型验证集验证过程的损失函数曲线图；
图4是根据基于LSTM算法网络模型预测的刀具磨损值得到的刀具预测磨损值曲线示意图。
具体实施方式
下面结合说明书附图和实施例对本发明作进一步限定，但不限于此，
实施例1
一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，如图1所示，包括步骤如下：
(1)获取多维刀具监测数据；
多维刀具监测数据包括测有刀具的磨损值的样本数据及未测有刀具的磨损值的样本数据；
刀具的磨损值来自于2010年PHM协会(Prognostics Health ManagementSociety)高速数控机床刀具健康预测竞赛的开源数据集。
数据采集条件如表1所示。
表1
Figure BDA0003147223200000071
Figure BDA0003147223200000081
每个样本数据均包括X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号，X维力、Y维力、Z维力是指在切削过程中x轴、y轴、z轴这三个方向上刀具受到的切削力，由Kistler9265B三向测力仪获得，X维振动、Y维振动、Z维振动是指在切削过程中x轴、y轴、z轴这三个方向上的振动信号(加速度)，x轴、y轴、z轴所处的坐标系为机床主轴坐标系，是根据右手笛卡尔直角坐标系建立的，主轴方向为z轴，垂直于主轴方向的平面上两个相互垂直的方向为x轴、y轴；声发射信号是金属在加工过程中由于内部分子晶格发生畸变、金属裂纹加剧以及其在塑性形变时所释放出来的一种超高频应力波脉冲信号。
基于表1的切削条件，重复进行6次全寿命周期试验。端面铣削材料是正方形，每次走刀端面铣的长度为108mm。控制每次走刀时间长度相同，每次走刀后测量的刀具磨损量来自后刀面。采集到的数据集共包含了C1-C6共6份数据，C1、C4、C6这3组数据中测量了铣刀的磨损量，C2、C3、C5没有测量，作为该次比赛的测试集。每份数据共315个样本，所有数据均通过一台高速数控机床在铣削作业下采集获得，每条样本由7个通道数据组成，如表2所示，表中的声发射信号是金属在加工过程中由于内部分子晶格发生畸变、金属裂纹加剧以及其在塑性形变时所释放出来的一种超高频应力波脉冲信号。
表2
Figure BDA0003147223200000082
使用的数据为测量了刀具磨损量的C1、C4、C6三组数据，数据原始文件为.csv文件，在这三组数据的磨损值文件中，记录了刀具在X、Y、Z三个方向的磨损。
(2)数据预处理；具体实现过程包括：读取测有刀具的磨损值的样本数据中的X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号这七通道数据并将其转换为DataFrame数据帧，转换后将这七个字段分别命名为Fx、Fy、Fz、Ax、Ay、Az、AE-rms；方便后续的特征提取工作。
将测有刀具的磨损值的样本数据中的刀具的磨损值在x轴、y轴、z轴这三个方向上的磨损量，每次铣削结束后停机，用LEICA MZ12显微镜测量该刀具的磨损值。求均值之后作为标签，得到若干标签矩阵；
将若干标签矩阵以.npy文件形式存储。以便后续使用，3组标签矩阵的具体形状均为[315,]。
(3)特征提取与融合；具体实现过程包括：
读取步骤(2)数据预处理后的数据帧，在时域上，计算绝对均值、峰值、均方根值、方根幅值、歪度值、峭度值、波形因子、脉冲因子、歪度因子、峰值因子、裕度因子、峭度因子12个特征；
绝对均值(Absolute mean)的计算公式为：
Figure BDA0003147223200000091
峰值(Max)的计算公式为：max(zi)；
均方根值(Root mean square)的计算公式为：
Figure BDA0003147223200000092
方根幅值(Square root amplitude)的计算公式为：
Figure BDA0003147223200000093
歪度值(Skewness)的计算公式为：
Figure BDA0003147223200000094
峭度值(Kurtosis)的计算公式为：
Figure BDA0003147223200000095
波形因子(Shape factor)的计算公式为：
Figure BDA0003147223200000096
脉冲因子(Pulse factor)的计算公式为：
Figure BDA0003147223200000101
歪度因子(Skewness factor)的计算公式为：
Figure BDA0003147223200000102
峰值因子(Crest factor)的计算公式为：
Figure BDA0003147223200000103
裕度因子(Clearance factor)的计算公式为：
Figure BDA0003147223200000104
峭度因子(Kurtosis factor)的计算公式为：
Figure BDA0003147223200000105
每一条数据长度为n，zij表示第i条数据第j个数据点，zi表示第i条数据，S(f)表示功率谱，max表示峰值；
在频域上，计算重心频率、均方频率、均方根频率、频率方差4个特征；
重心频率(FC)的计算公式为：
Figure BDA0003147223200000106
均方频率(MSF)的计算公式为：
Figure BDA0003147223200000107
均方根频率(RMSF)的计算公式为：
Figure BDA0003147223200000108
频率方差(VF)的计算公式为：
Figure BDA0003147223200000109
FC表示重心频率，S(f)为功率谱，f表示对输入的数据进行快速傅里叶变换之后得到的对应频率点，采样频率为1/50000；
在小波域上，用db3(多贝西极限相位小波)小波变换提取8个特征。分析低频和高频部分的数据，采用的小波树深度为3，重建小波包变换来分析不同频段的特征。8个节点系数被参数化为时频域中提取的8个特征。
经过特征提取后，共得到24维特征，将这24维特征融合为特征矩阵，这样可以很好地反映刀具的磨损状态。均以.npy文件形式存储。以便后续使用，每组数据特征矩阵形状的具体形状均为[315,6，24]。
(4)构造训练集、测试集；具体实现过程包括：
读取标签矩阵和特征矩阵，对其进行数据标准差标准化处理后，得到符合标准正态分布即均值为0、方差为1的数据；
把每个标签矩阵(因变量)和对应的特征矩阵(自变量)相匹配作为样本数据集；
将样本数据集的80％作为训练集，剩余部分作为测试集。
(5)搭建基于LSTM算法网络模型；具体是指：
如图2所示，基于LSTM算法网络模型为长短期记忆网络，长短期记忆网络包括输入层、第一隐藏层(隐含层1)、第二隐藏层(隐含层2)、全连接层及输出层；
输入层的输入特征维数为6，第一隐藏层、第二隐藏层的节点个数均为64；
将训练集输入第一隐藏层，第二隐藏层的输入为第一隐藏层的计算结果，在第二隐藏层设置dropout率为0.05的dropout层，避免过拟合问题，全连接层输入64个神经元，输出10个神经元，输出层采用线性激活函数，完成第二隐藏层到输出层的线性变换，使输出层输出刀具的磨损值。
(6)训练基于LSTM算法网络模型；具体是指：
将训练集输入基于LSTM算法网络模型，将标签作为基于LSTM算法网络模型的输出，对基于LSTM算法网络模型进行训练；训练的同时记录每一次训练周期的损失函数；
训练参数为：Epoch＝500，BatchSize＝128，Learningrate＝0.0001；
将验证集输入长短期记忆网络模型进行验证，在基于LSTM算法网络模型训练过程中用Adam算法进行优化，优化更新学习率Learningrate，验证时记录每一次训练周期的损失函数，得到训练好的基于LSTM算法网络模型。
将验证样本集输入基于LSTM算法网络模型进行验证，更新其参数，得到训练好的基于LSTM算法网络模型。
验证结束后，将所属测试样本集输入训练好的基于LSTM算法网络模型，得到刀具剩余寿命预测曲线并计算解释回归模型的方差得分(explained_variance_score,ev)、平均绝对误差(mean_absolute_error,mae)、均方误差(mean_squared_error，mse)和决定系数R方(r2_score，r2)评估模型优劣。，计算得到的回归指标如表3所示。
表3
ev mae mse r2
LSTM 0.770603 0.141811 0.026706 0.406032
计算解释训练好的基于LSTM算法网络模型的方差得分(explained_variance_score,ev)、平均绝对误差(mean_absolute_error,mae)、均方误差(mean_squared_error,mse)和决定系数R方(r2_score,r2)评估基于LSTM算法网络模型优劣。
图3(a)是基于LSTM算法网络模型训练集训练过程的损失函数曲线图；图3(b)是基于LSTM算法网络模型验证集验证过程的损失函数曲线图；图3(a)中，横坐标为训练周期epoch，纵坐标为训练集损失函数train loss；图3(b)中，横坐标为训练周期epoch，纵坐标为验证集损失函数valid loss。
(7)通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测。具体是指：
将测试集输入训练好的基于LSTM算法网络模型，得到预测磨损值，画出刀具预测磨损值曲线，刀具预测磨损值曲线中，横坐标为切削次数，纵坐标为磨损值。预测的刀具剩余寿命曲线如图4所示。图4中，横坐标为切削次数，纵坐标为刀具磨损值。
实施例3
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现实施例1或2所述基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。
实施例4
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现实施例1或2所述基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。

Claims (10)

1.一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，包括步骤如下：
(1)获取多维刀具监测数据；
多维刀具监测数据包括测有刀具的磨损值的样本数据及未测有刀具的磨损值的样本数据；
每个样本数据均包括X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号，X维力、Y维力、Z维力是指在切削过程中x轴、y轴、z轴这三个方向上刀具受到的切削力，X维振动、Y维振动、Z维振动是指在切削过程中x轴、y轴、z轴这三个方向上的振动信号，x轴、y轴、z轴所处的坐标系为机床主轴坐标系，是根据右手笛卡尔直角坐标系建立的，主轴方向为z轴，垂直于主轴方向的平面上两个相互垂直的方向为x轴、y轴；声发射信号是金属在加工过程中由于内部分子晶格发生畸变、金属裂纹加剧以及其在塑性形变时所释放出来的一种超高频应力波脉冲信号；
(2)数据预处理；
(3)特征提取与融合；
(4)构造训练集、测试集；
(5)搭建基于LSTM算法网络模型；
(6)训练基于LSTM算法网络模型；
(7)通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测。
2.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(2)中，对多维刀具监测数据数据预处理的具体实现过程包括：
读取测有刀具的磨损值的样本数据中的X维力、Y维力、Z维力、X维振动、Y维振动、Z维振动、声发射信号这七通道数据并将其转换为DataFrame数据帧，转换后将这七个字段分别命名为Fx、Fy、Fz、Ax、Ay、Az、AE-rms；
将测有刀具的磨损值的样本数据中的刀具的磨损值在x轴、y轴、z轴这三个方向上的磨损量，求均值之后作为标签，得到若干标签矩阵；
将若干标签矩阵以.npy文件形式存储。
3.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(3)中，对步骤(2)数据预处理后的数据进行特征提取的具体实现过程包括：
读取步骤(2)数据预处理后的数据帧，
在时域上，计算绝对均值、峰值、均方根值、方根幅值、歪度值、峭度值、波形因子、脉冲因子、歪度因子、峰值因子、裕度因子、峭度因子12个特征；
绝对均值的计算公式为：
Figure FDA0003147223190000021
峰值的计算公式为：max(zi)；
均方根值的计算公式为：
Figure FDA0003147223190000022
方根幅值的计算公式为：
Figure FDA0003147223190000023
歪度值的计算公式为：
Figure FDA0003147223190000024
峭度值的计算公式为：
Figure FDA0003147223190000025
波形因子的计算公式为：
Figure FDA0003147223190000026
脉冲因子的计算公式为：
Figure FDA0003147223190000027
歪度因子的计算公式为：
Figure FDA0003147223190000028
峰值因子的计算公式为：
Figure FDA0003147223190000029
裕度因子的计算公式为：
Figure FDA00031472231900000210
峭度因子的计算公式为：
Figure FDA0003147223190000031
每一条数据长度为n，zij表示第i条数据第j个数据点，zi表示第i条数据，S(f)表示功率谱，max表示峰值；
在频域上，计算重心频率、均方频率、均方根频率、频率方差4个特征；
重心频率的计算公式为：
Figure FDA0003147223190000032
均方频率的计算公式为：
Figure FDA0003147223190000033
均方根频率的计算公式为：
Figure FDA0003147223190000034
频率方差的计算公式为：
Figure FDA0003147223190000035
FC表示重心频率，S(f)为功率谱，f表示对输入的数据进行快速傅里叶变换之后得到的对应频率点，采样频率为1/50000；
在小波域上，用db3小波变换提取8个特征；
经过特征提取后，共得到24维特征，将这24维特征融合为特征矩阵，均以.npy文件形式存储。
4.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(4)中，构造训练集、测试集的具体实现过程包括：
读取标签矩阵和特征矩阵，对其进行数据标准差标准化处理后，得到符合标准正态分布即均值为0、方差为1的数据；
把每个标签矩阵和对应的特征矩阵相匹配作为样本数据集；
将样本数据集的部分数据作为训练集，剩余部分作为测试集。
5.根据权利要求4所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，将样本数据集的80％作为训练集，剩余部分作为测试集。
6.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(5)中，搭建基于LSTM算法网络模型，具体是指：
基于LSTM算法网络模型为长短期记忆网络，长短期记忆网络包括输入层、第一隐藏层、第二隐藏层、全连接层及输出层；
输入层的输入特征维数为6，第一隐藏层、第二隐藏层的节点个数均为64；
将训练集输入第一隐藏层，第二隐藏层的输入为第一隐藏层的计算结果，在第二隐藏层设置dropout率为0.05的dropout层，避免过拟合问题，全连接层输入64个神经元，输出10个神经元，输出层采用线性激活函数，完成第二隐藏层到输出层的线性变换，使输出层输出刀具的磨损值。
7.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(6)中，训练基于LSTM算法网络模型，具体是指：
将训练集输入基于LSTM算法网络模型，将标签作为基于LSTM算法网络模型的输出，对基于LSTM算法网络模型进行训练；训练的同时记录每一次训练周期的损失函数；
训练参数为：Epoch＝500，BatchSize＝128，Learningrate＝0.0001；
将验证集输入长短期记忆网络模型进行验证，在基于LSTM算法网络模型训练过程中用Adam算法进行优化，优化更新学习率Learningrate，验证时记录每一次训练周期的损失函数，得到训练好的基于LSTM算法网络模型。
8.根据权利要求1所述的一种基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法，其特征在于，步骤(7)中，通过训练好的基于LSTM算法网络模型进行刀具剩余寿命预测，具体是指：
将测试集输入训练好的基于LSTM算法网络模型，得到预测磨损值，画出刀具预测磨损值曲线，刀具预测磨损值曲线中，横坐标为切削次数，纵坐标为磨损值。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-7任一所述基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-7任一所述基于多维特征提取与长短期记忆网络的刀具剩余寿命预测方法的步骤。
CN202110755673.8A 2021-07-05 2021-07-05 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质 Pending CN113496312A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110755673.8A CN113496312A (zh) 2021-07-05 2021-07-05 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110755673.8A CN113496312A (zh) 2021-07-05 2021-07-05 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质

Publications (1)

Publication Number Publication Date
CN113496312A true CN113496312A (zh) 2021-10-12

Family

ID=77998187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110755673.8A Pending CN113496312A (zh) 2021-07-05 2021-07-05 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质

Country Status (1)

Country Link
CN (1) CN113496312A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113919396A (zh) * 2021-10-13 2022-01-11 山东大学 一种基于半监督学习的振动信号和图像特征机床刀具磨损状态监测方法
CN114091544A (zh) * 2021-11-29 2022-02-25 西南交通大学 基于振动信号与神经网络的tbm盘形滚刀磨损识别系统
CN114178905A (zh) * 2021-12-01 2022-03-15 广州先进技术研究所 一种铣刀磨损量预测方法和系统
CN114595733A (zh) * 2022-05-10 2022-06-07 山东大学 基于长短期记忆网络的桥梁拉索断丝信号识别方法及系统
CN114749996A (zh) * 2022-05-25 2022-07-15 哈尔滨工业大学 基于深度学习和时序回归模型的刀具剩余寿命预测方法
CN115922443A (zh) * 2023-01-09 2023-04-07 湖南大学 一种基于滑动窗口的刀具剩余使用寿命实时预测方法
CN116911469A (zh) * 2023-09-12 2023-10-20 成都飞机工业（集团）有限责任公司 一种数控机床加工时间预测方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108520320A (zh) * 2018-03-30 2018-09-11 华中科技大学 一种基于多重长短期记忆网络和经验贝叶斯的设备寿命预测方法
CN109753923A (zh) * 2018-12-29 2019-05-14 晋西车轴股份有限公司 刀具磨损量的监测方法、系统、设备及计算机可读存储介质
CN111475921A (zh) * 2020-03-13 2020-07-31 重庆邮电大学 一种基于边缘计算和lstm网络的刀具剩余寿命预测方法
CN111687689A (zh) * 2020-06-23 2020-09-22 重庆大学 一种基于lstm和cnn的刀具磨损状态预测方法及装置
CN111767521A (zh) * 2020-06-22 2020-10-13 中国石油化工股份有限公司 一种基于卷积神经网络与长短期记忆网络的输油泵滚动轴承状态评估方法
CN112613646A (zh) * 2020-12-08 2021-04-06 上海交通大学烟台信息技术研究院 一种基于多维数据融合的设备状态预测方法及系统
US20210117796A1 (en) * 2020-12-01 2021-04-22 Harbin Engineering University Ship Motion Prediction Method Based on Long Short-Term Memory Network and Gaussian Process Regression

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108520320A (zh) * 2018-03-30 2018-09-11 华中科技大学 一种基于多重长短期记忆网络和经验贝叶斯的设备寿命预测方法
CN109753923A (zh) * 2018-12-29 2019-05-14 晋西车轴股份有限公司 刀具磨损量的监测方法、系统、设备及计算机可读存储介质
CN111475921A (zh) * 2020-03-13 2020-07-31 重庆邮电大学 一种基于边缘计算和lstm网络的刀具剩余寿命预测方法
CN111767521A (zh) * 2020-06-22 2020-10-13 中国石油化工股份有限公司 一种基于卷积神经网络与长短期记忆网络的输油泵滚动轴承状态评估方法
CN111687689A (zh) * 2020-06-23 2020-09-22 重庆大学 一种基于lstm和cnn的刀具磨损状态预测方法及装置
US20210117796A1 (en) * 2020-12-01 2021-04-22 Harbin Engineering University Ship Motion Prediction Method Based on Long Short-Term Memory Network and Gaussian Process Regression
CN112613646A (zh) * 2020-12-08 2021-04-06 上海交通大学烟台信息技术研究院 一种基于多维数据融合的设备状态预测方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
薛嫣;朱静;翟怡萌;邓艾东;: "基于LSTM的风机滚动轴承剩余寿命预测研究", 工业控制计算机, no. 06 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113919396A (zh) * 2021-10-13 2022-01-11 山东大学 一种基于半监督学习的振动信号和图像特征机床刀具磨损状态监测方法
CN114091544A (zh) * 2021-11-29 2022-02-25 西南交通大学 基于振动信号与神经网络的tbm盘形滚刀磨损识别系统
CN114178905A (zh) * 2021-12-01 2022-03-15 广州先进技术研究所 一种铣刀磨损量预测方法和系统
CN114595733A (zh) * 2022-05-10 2022-06-07 山东大学 基于长短期记忆网络的桥梁拉索断丝信号识别方法及系统
CN114595733B (zh) * 2022-05-10 2024-04-02 山东大学 基于长短期记忆网络的桥梁拉索断丝信号识别方法及系统
CN114749996A (zh) * 2022-05-25 2022-07-15 哈尔滨工业大学 基于深度学习和时序回归模型的刀具剩余寿命预测方法
CN115922443A (zh) * 2023-01-09 2023-04-07 湖南大学 一种基于滑动窗口的刀具剩余使用寿命实时预测方法
CN116911469A (zh) * 2023-09-12 2023-10-20 成都飞机工业（集团）有限责任公司 一种数控机床加工时间预测方法
CN116911469B (zh) * 2023-09-12 2024-01-12 成都飞机工业（集团）有限责任公司 一种数控机床加工时间预测方法

Similar Documents

Publication Publication Date Title
CN113496312A (zh) 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质
US11761930B2 (en) Prediction method of part surface roughness and tool wear based on multi-task learning
CN111366123B (zh) 一种基于多任务学习的零件表面粗糙度和刀具磨损预测方法
CN111085898A (zh) 一种工况自适应的高速铣削加工过程刀具监测方法及系统
CN114619292B (zh) 一种基于小波降噪和注意力机制融合gru网络的铣削加工刀具磨损监测方法
CN113255848B (zh) 基于大数据学习的水轮机空化声信号辨识方法
CN114048600A (zh) 一种数字孪生驱动的多模型融合工业系统异常检测方法
CN110561191B (zh) 基于pca与自编码器的数控机床刀具磨损数据处理方法
CN112257530B (zh) 基于盲信号分离和支持向量机的滚动轴承故障诊断方法
CN111832624A (zh) 一种基于对抗迁移学习的刀具剩余寿命预测方法
CN112518425B (zh) 基于多源样本迁移强化学习的智能加工刀具磨损预测方法
CN116380445B (zh) 基于振动波形的设备状态诊断方法及相关装置
CN112207631B (zh) 刀具检测模型的生成方法、检测方法、系统、设备及介质
CN114871850B (zh) 一种基于振动信号和bp神经网络的刀具磨损状态评估方法
CN115741235A (zh) 基于五轴加工中心刀具的磨损预测与健康管理方法
CN116871978A (zh) 基于多传感器融合的钻削刀具状态监测方法
Wang et al. Self-adaptive fusion of local-temporal features for tool condition monitoring: A human experience free model
CN114354184A (zh) 一种基于深度学习的大型回转装备主轴健康预警模型建立方法和装置
Wang et al. Research on tool remaining life prediction method based on CNN-LSTM-PSO
CN116372665A (zh) 一种基于opc ua的航空数控机床刀具健康监测方法
Wei et al. Research on tool wear state identification method driven by multi-source information fusion and multi-dimension attention mechanism
Zhou et al. Milling cutter wear prediction based on bidirectional long short-term memory neural networks
CN114764538B (zh) 一种设备声音信号模式识别方法
Xu et al. Total process of fault diagnosis for wind turbine gearbox, from the perspective of combination with feature extraction and machine learning: A review
Song et al. Cnc machine tool wear monitoring based on densely connected convolutional networks

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