CN115922443B - 一种基于滑动窗口的刀具剩余使用寿命实时预测方法 - Google Patents

一种基于滑动窗口的刀具剩余使用寿命实时预测方法 Download PDF

Info

Publication number
CN115922443B
CN115922443B CN202310029401.9A CN202310029401A CN115922443B CN 115922443 B CN115922443 B CN 115922443B CN 202310029401 A CN202310029401 A CN 202310029401A CN 115922443 B CN115922443 B CN 115922443B
Authority
CN
China
Prior art keywords
value
matrix
sliding window
time
wear
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
Application number
CN202310029401.9A
Other languages
English (en)
Other versions
CN115922443A (zh
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.)
Hunan University
Original Assignee
Hunan 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 Hunan University filed Critical Hunan University
Priority to CN202310029401.9A priority Critical patent/CN115922443B/zh
Publication of CN115922443A publication Critical patent/CN115922443A/zh
Application granted granted Critical
Publication of CN115922443B publication Critical patent/CN115922443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Machine Tool Sensing Apparatuses (AREA)

Abstract

本发明公开了一种基于滑动窗口的刀具剩余使用寿命实时预测方法，包括采集刀具切削过程中的切削力信号、振动信号及声发射信号并经采集卡转换为数字信号，每完成一次工件表面铣削，用显微镜离线测量刀具磨损值；对上述数字信号分别提取时域、频域和时频域特征，形成特征矩阵并进行规范化；设计滑动窗口提取全寿命已知的磨损值，与规范化后的特征矩阵形成输入矩阵对长短时记忆网络模型进行训练；按时序迭代输入前一时刻预测出的磨损值更新所提取的磨损值滑动窗，与规范化后的特征矩阵形成输入矩阵对训练好的长短时记忆网络模型进行验证，得到对应时刻刀具磨损量的预测值。只需要切削过程中的上述信号和前期少量测量磨损值，就能预测剩余使用寿命。

Description

一种基于滑动窗口的刀具剩余使用寿命实时预测方法
技术领域
本发明属于装备健康管理技术领域，特别是涉及一种基于滑动窗口的刀具剩余使用寿命实时预测方法。
背景技术
刀具作为加工系统中直接与工件接触的部分，对产品质量有着直接且重要的影响，刀具磨损退化会引起产品质量与精度的降低，甚至还有可能造成磨床故障甚至安全事故。因此在加工过程中对刀具使用状态进行实时监测并对其可用剩余寿命做出准确预测，在刀具失效前及时换刀，对保证加工质量和安全有重要意义。然而实际情况中，刀具退化过程存在一定差异性和随机性，传统基于失效机制或损伤第一原理建立数学模型来描述机械退化过程的方法严重依赖模型准确性，且实际应用中难以建立准确的数学模型导致这类方法难以满足个体刀具健康状态评估要求；目前基于数据驱动的寿命预测方法经常忽略传感器采集信号的时序特性，未充分利用时序信息导致预测精度低；此外，目前大部分方法需要已知刀具磨损全寿命周期数据，进行历史数据拟合，而未能充分利用刀具传感器实时数据进行剩余使用寿命实时预测。因此，亟需开发考虑传感器数据时序特征且实时的寿命预测方法。
发明内容
针对以上技术问题，本发明提供一种基于滑动窗口的刀具剩余使用寿命实时预测方法。
本发明解决其技术问题采用的技术方案是：
一种基于滑动窗口的刀具剩余使用寿命实时预测方法，方法包括以下步骤：
S100：测量并记录刀具切削过程中的
Figure 129338DEST_PATH_IMAGE001
通道切削力信号、振动信号及声发射信号并经采集卡转换为数字信号，每完成一次工件表面铣削，用显微镜离线测量刀具磨损值，将所有的刀具磨损值数据表达成预设矩阵形式；
S200：对切削力数字信号、振动数字信号及声发射数字信号分别提取时域、频域和时频域特征，形成特征矩阵，对特征矩阵进行规范化得到规范化后的特征矩阵；
S300：在模型训练阶段，设计滑动窗提取全寿命已知的磨损值，与规范化后的特征矩阵形成输入矩阵对长短时记忆网络模型进行训练，得到训练好的长短时记忆网络模型；
S400：在模型测试阶段，设计滑动窗提取磨损值，按时序迭代输入前一时刻预测出的磨损值更新所提取的磨损值滑动窗，与规范化后的特征矩阵形成输入矩阵对训练好的长短时记忆网络模型进行验证，得到对应时刻刀具磨损量的预测值。
优选地，S100中将所有的刀具磨损值数据表达成预设矩阵形式，具体为：
Figure 207015DEST_PATH_IMAGE002
其中，
Figure 849349DEST_PATH_IMAGE003
为切削次数，/>
Figure 594451DEST_PATH_IMAGE004
为切削/>
Figure 31249DEST_PATH_IMAGE003
次的磨损值向量。
优选地，时域特征包括绝对平均值、最大值、均方根、平方根振幅、偏斜、峰度、形状因子、脉冲因子、偏斜因子、波峰因子、间隙系数和峰度因子，频域特征包括重心频率、均方频率、均方根频率和频率方差。
优选地，S200包括：
S210：以一次切削过程所采集的信号为一个时间单位，每一次切削，每个通道提取时域特征、频域特征和时频域特征数据共
Figure 346955DEST_PATH_IMAGE005
个，横向拼接形成特征矩阵/>
Figure 527400DEST_PATH_IMAGE006
；
S220：根据特征矩阵计算得到每行特征数据的平均值和标准差，根据每行特征数据的平均值和标准差对特征矩阵进行规范化，具体为：
Figure 127009DEST_PATH_IMAGE007
其中，
Figure 734708DEST_PATH_IMAGE008
为规范化后的特征数据，/>
Figure 36245DEST_PATH_IMAGE009
为提取的时域特征、频域特征和时频域特征数据，/>
Figure 20381DEST_PATH_IMAGE010
为每行特征数据的平均值，/>
Figure 474497DEST_PATH_IMAGE011
为每行特征数据的标准差，/>
Figure 253097DEST_PATH_IMAGE012
，/>
Figure 58242DEST_PATH_IMAGE013
，/>
Figure 328293DEST_PATH_IMAGE014
；
将规范化后的特征矩阵的
Figure 636914DEST_PATH_IMAGE015
行横向拼接重塑成/>
Figure 586416DEST_PATH_IMAGE016
形式。
优选地，S300包括：
S310：以滑动窗形式提取全寿命已知的磨损值，对于第
Figure 613277DEST_PATH_IMAGE017
次切削，提取的磨损值窗为：
Figure 939217DEST_PATH_IMAGE018
其中，
Figure 351612DEST_PATH_IMAGE019
是预设滑动窗口长度，/>
Figure 472015DEST_PATH_IMAGE020
；
S320：将第
Figure 251752DEST_PATH_IMAGE017
次切削提取的磨损值滑动窗对应拼接到第/>
Figure 115803DEST_PATH_IMAGE017
次切削的特征向量后，对应特征矩阵形式变为/>
Figure 884170DEST_PATH_IMAGE021
，即为输入矩阵；
S330：提取从第
Figure 441053DEST_PATH_IMAGE019
次切削之后的磨损值/>
Figure 176928DEST_PATH_IMAGE022
为输出矩阵，与输入矩阵格式对应，具体为：
Figure 579091DEST_PATH_IMAGE023
步骤S340：使用输入矩阵和输出矩阵训练LSTM网络模型。
优选地，S400包括：
S410：滑动窗提取前
Figure 966078DEST_PATH_IMAGE019
个已知磨损值/>
Figure 693863DEST_PATH_IMAGE024
，横向拼接到第/>
Figure 651455DEST_PATH_IMAGE025
次切削的特征向量里，即第/>
Figure 388467DEST_PATH_IMAGE025
次的输入矩阵形式为/>
Figure 380693DEST_PATH_IMAGE026
，输入到训练好的LSTM模型，预测出第/>
Figure 761603DEST_PATH_IMAGE025
次的磨损值/>
Figure 472070DEST_PATH_IMAGE027
；
S420： 将
Figure 747193DEST_PATH_IMAGE027
按时序拼接到磨损值滑动窗末尾，同时剔除掉滑动窗的第一个值保持滑动窗长度为/>
Figure 593927DEST_PATH_IMAGE019
不变，将新的磨损值/>
Figure 912781DEST_PATH_IMAGE028
拼接到第/>
Figure 110545DEST_PATH_IMAGE025
次切削的特征向量里，输入到训练好的LSTM模型中，预测出第/>
Figure 658201DEST_PATH_IMAGE029
次的磨损值/>
Figure 890599DEST_PATH_IMAGE030
；
S430：每预测出一个新的磨损值就滑动更新磨损值滑动窗，拼接到下一时刻需要预测的磨损值的特征向量里，滚动迭代预测直至第
Figure 865508DEST_PATH_IMAGE031
次切削。
优选地，S100之后还包括对表达成预设矩阵形式的所有的刀具磨损值标签进行归一化，具体为：
Figure 301300DEST_PATH_IMAGE032
其中，
Figure 918226DEST_PATH_IMAGE033
为归一化后的磨损值标签数据，/>
Figure 739551DEST_PATH_IMAGE034
为测量的磨损值，/>
Figure 885362DEST_PATH_IMAGE035
为磨损值列向量中最小值，/>
Figure 41405DEST_PATH_IMAGE036
为磨损值列向量中最大值。
上述一种基于滑动窗口的刀具剩余使用寿命实时预测方法，带入之前时刻预测出的磨损值进行下一步预测，不需要测量全过程的磨损值，与其他用神经网络做预测的方法相比，本发明提出的数据输入考虑了数据的时序特征，提高了寿命预测精度，利用切削力、振动、声发射信号及前期少量测量磨损值，实现了三槽球头硬质合金铣刀的剩余使用寿命预测。
附图说明
图1为本发明一实施例提供的一种基于滑动窗口的刀具剩余使用寿命实时预测方法流程图；
图2为本发明一实施例提供的滑动窗思想提取磨损值示意图；
图3为本发明一实施例提供的LSTM网络原理结构图；
图4为本发明一实施例中针对只用时域、频域、时频域三种特征作为LSTM网络输入的预测结果示意图；
图5为本发明一实施例中针对用时域、频域、时频域三种特征加上滑动窗提取的磨损值作为LSTM网络输入的预测结果示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明的技术方案，下面结合附图对本发明作进一步的详细说明。
在一个实施例中，如图1所示，一种基于滑动窗口的刀具剩余使用寿命实时预测方法，方法包括以下步骤：
S100：测量并记录刀具切削过程中的
Figure 462023DEST_PATH_IMAGE001
通道切削力信号、振动信号及声发射信号并经采集卡转换为数字信号，每完成一次工件表面铣削，用显微镜离线测量刀具磨损值，将所有的刀具磨损值数据表达成预设矩阵形式；
S200：对切削力数字信号、振动数字信号及声发射数字信号分别提取时域、频域和时频域特征，形成特征矩阵，对特征矩阵进行规范化得到规范化后的特征矩阵；
S300：在模型训练阶段，设计滑动窗提取全寿命已知的磨损值，与规范化后的特征矩阵形成输入矩阵对长短时记忆网络模型进行训练，得到训练好的长短时记忆网络模型；
S400：在模型测试阶段，设计滑动窗提取磨损值，按时序迭代输入前一时刻预测出的磨损值更新所提取的磨损值滑动窗，与规范化后的特征矩阵形成输入矩阵对训练好的长短时记忆网络模型进行验证，得到对应时刻刀具磨损量的预测值。
与其他寿命预测方法相比较，本发明提出的寿命预测方法考虑了数据的时序特征，并且带入预测出的磨损值进行下一步预测。利用切削力信号、振动信号、声发射信号和前期少量磨损值，实现了磨床刀具的剩余使用寿命预测。
在一个实施例中，S100中将所有的刀具磨损值数据表达成预设矩阵形式，具体为：
Figure 137855DEST_PATH_IMAGE002
其中，
Figure 720146DEST_PATH_IMAGE003
为切削次数，/>
Figure 114218DEST_PATH_IMAGE004
为切削/>
Figure 844187DEST_PATH_IMAGE003
次的磨损值向量。
具体地，采集X、Y、Z三个轴向的切削力信号和振动信号，及铣削过程中的声发射信号，共计7个通道的信号。对于训练集，测量全过程315次磨损值，对于测试集，测量前10次工件表面铣削的刀具磨损值。
在一个实施例中，S200包括：
S210：以一次切削过程所采集的信号为一个时间单位，每一次切削，每个通道提取时域特征、频域特征和时频域特征数据共
Figure 640104DEST_PATH_IMAGE005
个，横向拼接形成特征矩阵/>
Figure 393297DEST_PATH_IMAGE006
；
S220：根据特征矩阵计算得到每行特征数据的平均值和标准差，根据每行特征数据的平均值和标准差对特征矩阵进行规范化，具体为：
Figure 9086DEST_PATH_IMAGE007
其中，
Figure 771506DEST_PATH_IMAGE008
为规范化后的特征数据，/>
Figure 671197DEST_PATH_IMAGE009
为提取的时域特征、频域特征和时频域特征数据，/>
Figure 595291DEST_PATH_IMAGE010
为每行特征数据的平均值，/>
Figure 963955DEST_PATH_IMAGE011
为每行特征数据的标准差，/>
Figure 264487DEST_PATH_IMAGE012
，/>
Figure 520150DEST_PATH_IMAGE013
，/>
Figure 615145DEST_PATH_IMAGE014
；
将规范化后的特征矩阵的
Figure 205526DEST_PATH_IMAGE015
行横向拼接重塑成/>
Figure 44169DEST_PATH_IMAGE016
形式。
在一个实施例中，时域特征包括绝对平均值、最大值、均方根、平方根振幅、偏斜、峰度、形状因子、脉冲因子、偏斜因子、波峰因子、间隙系数和峰度因子，频域特征包括重心频率、均方频率、均方根频率和频率方差。
具体地，对7个通道数据，每个通道提取时域特征（绝对平均值、最大值、均方根、平方根振幅、偏斜、峰度、形状因子、脉冲因子、偏斜因子、波峰因子、间隙系数、峰度因子）12个，频域特征（重心频率、均方频率、均方根频率、频率方差）4个，时频域特征8个，共24个特征，形成特征矩阵
Figure 934765DEST_PATH_IMAGE037
。将特征矩阵进行规范化后，重塑成/>
Figure 449929DEST_PATH_IMAGE038
的形式。
在一个实施例中，S300包括：
S310：以滑动窗形式提取全寿命已知的磨损值，对于第
Figure 793185DEST_PATH_IMAGE017
次切削，提取的磨损值窗为：
Figure 169940DEST_PATH_IMAGE018
其中，
Figure 649463DEST_PATH_IMAGE019
是预设滑动窗口长度，/>
Figure 834063DEST_PATH_IMAGE020
；
S320：将第
Figure 664616DEST_PATH_IMAGE017
次切削提取的磨损值滑动窗对应拼接到第/>
Figure 845062DEST_PATH_IMAGE017
次切削的特征向量后，对应特征矩阵形式变为/>
Figure 444670DEST_PATH_IMAGE021
，即为输入矩阵；
S330：提取从第
Figure 317948DEST_PATH_IMAGE019
次切削之后的磨损值/>
Figure 353906DEST_PATH_IMAGE022
为输出矩阵，与输入矩阵格式对应，具体为：
Figure 338043DEST_PATH_IMAGE023
步骤S340：使用输入矩阵和输出矩阵训练LSTM网络模型。
具体地，滑动窗形式提取磨损值的方法见图2所示，对于训练集，预设滑动窗口长度为10，第一个滑动窗为
Figure 792158DEST_PATH_IMAGE039
，拼接到第11次切削的特征向量里，迭代至第315次切削，故输入矩阵大小变为/>
Figure 836337DEST_PATH_IMAGE040
，输出矩阵也从第11次切削之后提前，大小变为
Figure 110324DEST_PATH_IMAGE041
。用/>
Figure 648884DEST_PATH_IMAGE042
作为输入、/>
Figure 957505DEST_PATH_IMAGE043
作为输出，训练LSTM网络，其中LSTM网络内部原理结构见图3所示。
在一个实施例中，S400包括：
S410：滑动窗提取前
Figure 907007DEST_PATH_IMAGE019
个已知磨损值/>
Figure 933868DEST_PATH_IMAGE024
，横向拼接到第/>
Figure 259808DEST_PATH_IMAGE025
次切削的特征向量里，即第/>
Figure 937782DEST_PATH_IMAGE025
次的输入矩阵形式为/>
Figure 58185DEST_PATH_IMAGE026
，输入到训练好的LSTM模型，预测出第/>
Figure 572343DEST_PATH_IMAGE025
次的磨损值/>
Figure 436394DEST_PATH_IMAGE027
；
S420： 将
Figure 467411DEST_PATH_IMAGE027
按时序拼接到磨损值滑动窗末尾，同时剔除掉滑动窗的第一个值保持滑动窗长度为/>
Figure 758715DEST_PATH_IMAGE019
不变，将新的磨损值/>
Figure 760169DEST_PATH_IMAGE028
拼接到第/>
Figure 427910DEST_PATH_IMAGE025
次切削的特征向量里，输入到训练好的LSTM模型中，预测出第/>
Figure 831210DEST_PATH_IMAGE029
次的磨损值/>
Figure 542683DEST_PATH_IMAGE030
；
S430：每预测出一个新的磨损值就滑动更新磨损值滑动窗，拼接到下一时刻需要预测的磨损值的特征向量里，滚动迭代预测直至第
Figure 765854DEST_PATH_IMAGE031
次切削。
具体地，对于测试集，第一个磨损值滑动窗为测量值
Figure 237286DEST_PATH_IMAGE044
，拼接到第11次切削的特征向量里，输入训练好的LSTM模型中预测出第11次磨损值/>
Figure 229513DEST_PATH_IMAGE045
，将/>
Figure 128199DEST_PATH_IMAGE045
带入磨损值滑动窗，并剔除掉/>
Figure 589398DEST_PATH_IMAGE046
，即滑动窗变为/>
Figure 598943DEST_PATH_IMAGE047
，拼接到第12次切削的特征向量里，输入LSTM网络模型预测出/>
Figure 711255DEST_PATH_IMAGE048
，将/>
Figure 780842DEST_PATH_IMAGE048
带入磨损值滑动窗预测
Figure 962294DEST_PATH_IMAGE049
，由此类推一直到第315次切削结束。
在一个实施例中，S100之后还包括对表达成预设矩阵形式的所有的刀具磨损值标签进行归一化，具体为：
Figure 41108DEST_PATH_IMAGE032
其中，
Figure 742348DEST_PATH_IMAGE033
为归一化后的磨损值标签数据，/>
Figure 982837DEST_PATH_IMAGE034
为测量的磨损值，/>
Figure 667896DEST_PATH_IMAGE035
为磨损值列向量中最小值，/>
Figure 767045DEST_PATH_IMAGE036
为磨损值列向量中最大值。
具体地，磨损值标签归一化后仍为
Figure 588371DEST_PATH_IMAGE050
的形式。
在一个详细的实施例中，本发明在高速数控机床Roders Tech RFM760上采用三槽球头硬质合金铣刀进行铣削试验，被切削工件材料为硬度HRC52不锈钢，主轴转速10400 r/min，进给率1555mm/min，径向切深0.2mm，轴向切深0.125mm。在切削过程中，以50 kHz采样速率采集切削力、振动及声发射信号共7个通道的信号。每完成一次工件表面铣削，用莱卡MZ12型号显微镜离线测量刀具磨损量，直至现场专业人员认为磨损量超出正常范围时停止试验，共完成315 次切削。图4为只输入时域、频域、时频域三种信号进行预测曲线图，图5为输入三种信号加磨损值滑动窗预测结果图，对比可以看到采用本发明中提出的考虑时序特征、加入迭代更新的预测出的磨损值进行下一步预测，预测曲线要更加平滑、贴近真实值、无首尾突变值，可以有效提升寿命预测准确率，表1为不同情形下预测结果的均方根误差，不考虑时序特征方法的输入为时域、频域、时频域特征，本方法考虑时序特征则是在普通方法基础上加上了磨损值滑动窗输入。
表1 不同输入下预测结果
Figure 999761DEST_PATH_IMAGE051
上述刀具剩余使用寿命预测方法，与其他寿命预测方法相比较，本发明提出的寿命预测方法带入之前时刻预测出的磨损值进行下一步预测，不需要测量全过程的磨损值。与其他用神经网络做预测的方法相比，本发明提出的数据输入考虑了数据的时序特征，提高了寿命预测精度。利用切削力、振动、声发射信号及前期少量测量磨损值，实现了三槽球头硬质合金铣刀的剩余使用寿命预测。
以上对本发明所提供的一种基于滑动窗口的刀具剩余使用寿命实时预测方法进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (3)

1.一种基于滑动窗口的刀具剩余使用寿命实时预测方法，其特征在于，所述方法包括以下步骤：
S100：测量并记录刀具切削过程中的
Figure QLYQS_1
通道切削力信号、振动信号及声发射信号并经采集卡转换为数字信号，每完成一次工件表面铣削，用显微镜离线测量刀具磨损值，将所有的刀具磨损值数据表达成预设矩阵形式；
S100中将所有的刀具磨损值数据表达成预设矩阵形式，具体为：
Figure QLYQS_2
其中，/>
Figure QLYQS_3
为切削次数，/>
Figure QLYQS_4
为切削/>
Figure QLYQS_5
次的磨损值向量；
S200：对切削力数字信号、振动数字信号及声发射数字信号分别提取时域、频域和时频域特征，形成特征矩阵，对所述特征矩阵进行规范化得到规范化后的特征矩阵；
S200包括：
S210：以一次切削过程所采集的信号为一个时间单位，每一次切削，每个通道提取时域特征、频域特征和时频域特征数据共
Figure QLYQS_6
个，横向拼接形成特征矩阵/>
Figure QLYQS_7
；
S220：根据所述特征矩阵计算得到每行特征数据的平均值和标准差，根据所述每行特征数据的平均值和标准差对所述特征矩阵进行规范化，具体为：
Figure QLYQS_8
其中，/>
Figure QLYQS_12
为规范化后的特征数据，/>
Figure QLYQS_14
为提取的时域特征、频域特征和时频域特征数据，/>
Figure QLYQS_10
为每行特征数据的平均值，/>
Figure QLYQS_11
为每行特征数据的标准差，/>
Figure QLYQS_13
，/>
Figure QLYQS_15
，
Figure QLYQS_9
；
将规范化后的特征矩阵的
Figure QLYQS_16
行横向拼接重塑成/>
Figure QLYQS_17
形式；
S300：在模型训练阶段，设计滑动窗提取全寿命已知的磨损值，与所述规范化后的特征矩阵形成输入矩阵对长短时记忆网络模型进行训练，得到训练好的长短时记忆网络模型；
S300包括：
S310：以滑动窗形式提取全寿命已知的磨损值，对于第
Figure QLYQS_18
次切削，提取的磨损值窗为：
Figure QLYQS_19
其中，/>
Figure QLYQS_20
是预设滑动窗口长度，/>
Figure QLYQS_21
；
S320：将第
Figure QLYQS_22
次切削提取的磨损值滑动窗对应拼接到第/>
Figure QLYQS_23
次切削的特征向量后，对应特征矩阵形式变为/>
Figure QLYQS_24
，即为输入矩阵；
S330：提取从第
Figure QLYQS_25
次切削之后的磨损值/>
Figure QLYQS_26
为输出矩阵，与所述输入矩阵格式对应，具体为：
Figure QLYQS_27
步骤S340：使用所述输入矩阵和所述输出矩阵训练LSTM网络模型；
S400：在模型测试阶段，设计滑动窗提取磨损值，按时序迭代输入前一时刻预测出的磨损值更新所提取的磨损值滑动窗，与所述规范化后的特征矩阵形成输入矩阵对所述训练好的长短时记忆网络模型进行验证，得到对应时刻刀具磨损量的预测值；
S400包括：
S410：滑动窗提取前
Figure QLYQS_28
个已知磨损值/>
Figure QLYQS_29
，横向拼接到第/>
Figure QLYQS_30
次切削的特征向量里，即第/>
Figure QLYQS_31
次的输入矩阵形式为/>
Figure QLYQS_32
，输入到训练好的LSTM模型，预测出第/>
Figure QLYQS_33
次的磨损值/>
Figure QLYQS_34
；
S420： 将
Figure QLYQS_35
按时序拼接到磨损值滑动窗末尾，同时剔除掉滑动窗的第一个值保持滑动窗长度为/>
Figure QLYQS_36
不变，将新的磨损值/>
Figure QLYQS_37
拼接到第/>
Figure QLYQS_38
次切削的特征向量里，输入到训练好的LSTM模型中，预测出第/>
Figure QLYQS_39
次的磨损值/>
Figure QLYQS_40
；
S430：每预测出一个新的磨损值就滑动更新磨损值滑动窗，拼接到下一时刻需要预测的磨损值的特征向量里，滚动迭代预测直至第
Figure QLYQS_41
次切削。
2.根据权利要求1所述的方法，其特征在于，所述时域特征包括绝对平均值、最大值、均方根、平方根振幅、偏斜、峰度、形状因子、脉冲因子、偏斜因子、波峰因子、间隙系数和峰度因子，所述频域特征包括重心频率、均方频率、均方根频率和频率方差。
3.根据权利要求1所述的方法，其特征在于，S100之后还包括对表达成预设矩阵形式的所有的刀具磨损值标签进行归一化，具体为：
Figure QLYQS_42
其中，/>
Figure QLYQS_43
为归一化后的磨损值标签数据，/>
Figure QLYQS_44
为测量的磨损值，/>
Figure QLYQS_45
为磨损值列向量中最小值，/>
Figure QLYQS_46
为磨损值列向量中最大值。
CN202310029401.9A 2023-01-09 2023-01-09 一种基于滑动窗口的刀具剩余使用寿命实时预测方法 Active CN115922443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310029401.9A CN115922443B (zh) 2023-01-09 2023-01-09 一种基于滑动窗口的刀具剩余使用寿命实时预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310029401.9A CN115922443B (zh) 2023-01-09 2023-01-09 一种基于滑动窗口的刀具剩余使用寿命实时预测方法

Publications (2)

Publication Number Publication Date
CN115922443A CN115922443A (zh) 2023-04-07
CN115922443B true CN115922443B (zh) 2023-06-16

Family

ID=86552550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310029401.9A Active CN115922443B (zh) 2023-01-09 2023-01-09 一种基于滑动窗口的刀具剩余使用寿命实时预测方法

Country Status (1)

Country Link
CN (1) CN115922443B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116551467B (zh) * 2023-07-05 2023-09-01 南京维拓科技股份有限公司 基于物联网监控数据的刀具磨损状态识别方法
CN117961646A (zh) * 2024-03-29 2024-05-03 赣州海盛硬质合金有限公司 一种重载加工用数控刀片剩余使用寿命预测方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021174525A1 (zh) * 2020-03-06 2021-09-10 大连理工大学 一种基于多任务学习的零件表面粗糙度和刀具磨损预测方法
CN112247674B (zh) * 2020-10-10 2021-09-21 北京理工大学 一种刀具磨损预测方法
CN113569903B (zh) * 2021-06-09 2024-04-09 西安电子科技大学 数控机床刀具磨损预测方法、系统、设备、介质、终端
CN113305645B (zh) * 2021-06-22 2022-07-15 重庆邮电大学工业互联网研究院 一种基于混合神经模型的数控机床刀具剩余寿命预测方法
CN113496312A (zh) * 2021-07-05 2021-10-12 山东大学 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质
CN113779724B (zh) * 2021-09-13 2022-07-22 江苏仅一联合智造有限公司 一种充填包装机故障智能预测方法及其系统

Also Published As

Publication number Publication date
CN115922443A (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN115922443B (zh) 一种基于滑动窗口的刀具剩余使用寿命实时预测方法
TWI640390B (zh) 刀具磨耗監測與預測方法
CN110561193B (zh) 一种基于特征融合的刀具磨损评估及监控的方法与系统
US6845340B2 (en) System and method for machining data management
CN113569903B (zh) 数控机床刀具磨损预测方法、系统、设备、介质、终端
CN113378725B (zh) 一种基于多尺度-通道注意力网络的刀具故障诊断方法、设备及存储介质
JP6832327B2 (ja) 生産プロセス中に異常なワークピースを自動検出するためのデータ駆動方法
CN109333160B (zh) 高温合金钻削过程钻头磨损形式及磨损状态的在线监测方法
CN106514434A (zh) 一种基于数据的铣削刀具磨损监测方法
CN110928237B (zh) 一种基于振动信号的数控加工中心颤振在线辨识方法
CN113496312A (zh) 基于多维特征提取融合与长短期记忆网络的刀具剩余寿命预测方法、设备及存储介质
CN112171376B (zh) 一种基于电流信号分割的机床工件实时统计方法
CN112160761B (zh) 一种基于现场参数的硬岩tbm刀盘滚刀磨损实时评估方法
CN114326593B (zh) 一种刀具寿命预测系统及方法
CN113762182A (zh) 基于深度网络自适应的铣削刀具磨损状态预测方法及装置
CN113505531A (zh) 基于结合传统诊断与机器学习方法且针对往复机械气门间隙异常故障的诊断模型构建方法
CN111993158B (zh) 刀具检测模型的生成方法、检测方法、系统、设备及介质
CN115344951A (zh) 基于时间卷积网络和辅助学习的刀具磨损量预测方法
CN117577137B (zh) 切刀健康评估方法、装置、设备及存储介质
CN111428440B (zh) 一种基于条件概率的时序日志样本自动标注方法及装置
CN116187148A (zh) 设备寿命预测及模型训练方法、装置和电子设备
JP7370044B2 (ja) 監視条件決定装置、監視条件決定方法、及びプログラム
CN117943892A (zh) 基于时序数据分析的刀具监测数据预测优化方法
CN116805041A (zh) 一种刀具磨损量预测方法及系统
CN116900815A (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