CN109318055B - 一种铣削刀具磨损状态特征提取多目标优化方法 - Google Patents

Abstract

本发明公开了一种铣削刀具磨损状态特征提取多目标优化方法，包括以下步骤：通过多种传感器通道采集多个物理场信号；计算每个传感通道信号的若干时域和频域统计特征参数、以及小波能量，构成特征参数候选集；以刀具状态预测精度和特征参数数量为优化目标，构造多目标优化模型；采用智能优化算法对优化模型进行全局寻优；将优化模型最优解所对应的传感特征参数集作为刀具状态监测所需的特征参数。本发明具有以下优点和效果：本发明从预测精度出发，以预测精度和特征参数个数为优化目标，考察各特征参数组合对应的刀具磨损量预测精度，避免了相关性强但预测精度低的现象。

Description

一种铣削刀具磨损状态特征提取多目标优化方法

技术领域

本发明涉及制造过程监测领域，特别涉及一种铣削刀具磨损状态特征提取多目标优化方法。

背景技术

随着市场竞争的日益激烈，制造型企业对生产过程的自动化需求不断增加，数控铣床的自动化是大多数制造过程自动化的重要组成部分。铣削刀具作为数控铣床最易损伤的部件，对其进行及时有效的状态监测与故障识别至关重要。

近年来，基于多传感特征融合的间接式铣削刀具状态监测(TCM)方法受到了国内外学者的广泛关注。基于多传感特征融合的TCM方法是通过多个物理场传感器获取切削加工过程的相关信号，借助信号处理和统计分析技术对刀具状态进行估计，以达到实时监测刀具运行状态的目的。目前，学者们开展了大量TCM研究，已提出了诸多比较有效的TCM的方法，这为高精度、高可靠性的TCM提供了一定的技术基础。然而，目前的大多数相关研究存在如下两个主要问题：

(1)多传感信号的组合分析固然可以提高TCM识别精度，但传感器并非越多越好，冗余信号过多反而会影响TCM的识别精度，亦会导致监测成本和维护成本的增加；

(2)目前众多研究所采用的特征选取算法，均是基于与刀具状态的线性/非线性相关度来考察。然而，相关性强的特征参数并不一定能获得良好的监测性能。

发明内容

本发明针对上述现有技术的不足，提供了一种铣削刀具状态特征提取多目标优化方法。

本发明的上述技术目的是通过以下技术方案得以实现的：一种铣削刀具磨损状态特征提取多目标优化方法，包括以下步骤：

S1、通过多种传感器通道采集多个物理场信号；

S2、计算每个传感通道信号的若干时域和频域统计特征参数、以及小波能量，构成特征参数候选集；

S3、以刀具状态预测精度和特征参数数量为优化目标，构造多目标优化模型；

S4、采用智能优化算法对优化模型进行全局寻优；

S5、将优化模型最优解所对应的传感特征参数集作为刀具状态监测所需的特征参数。

进一步设置是所述的步骤S1具体包括：

周期性采集刀具运行状态下的S个物理场传感通道的时域信号以及对应的刀具磨损量，共分别采集T次，组成训练样本集。

进一步设置是所述的步骤S2具体包括：

S21、计算每个传感通道x_it(i＝1,…,M)的F个统计特征参数，组成特征参数集f_it＝(f_it-1,f_it-2,...,f_it-F)；

S22、将所有S个传感通道的特征参数集合并，构造特征参数候选集F＝{f_it}。

进一步设置是所述的步骤S3具体包括：构造多目标优化模型

其中，Y_t′表示在某种特征参数组合F_δq下通过预测算法得到的X_t的预测值。

进一步设置是所述的预测算法包括如下步骤：

S301、找出{δ_q}(q＝1,…,Q)对应特征参数组合下的样本数据集，即选择{δ_q}中值为1的位置对应的特征参数；

S302、计算样本核矩阵

其中

为X_δq中两个样本点x_i和x_j的核函数值；

S303、按照下式计算X_t的预测值Y_t′：

其中，I表示单位矩阵，c表示正则化参数。

进一步设置是所述的步骤S4具体包括：

所述的智能优化算法包括粒子群算法、遗传算法和蚁群算法，并求得最优的特征参数集F*。

进一步设置是所述的步骤S5具体包括：

S51、计算F*对应的样本核矩阵

S52、在线监测，对于未知信号x，其对应的刀具磨损量预测值计算如下：

本发明的有益效果是：

(1)目前，对TCM的研究大多集中于多传感信号的组合分析。多传感信号固然可以提高TCM识别精度，但传感器并非越多越好，冗余信息过多反而会影响TCM的识别精度，亦会导致监测成本和维护成本的增加。本发明通过对“传感器—特征参数”的多目标组合优化分析，结合智能进化算法找出最佳的传感参数配置，如若某一传感通道对应的特征参数无一被选取，则该传感通道可以不予考虑，从而在保证预测精度的基础上获得最少传感数量。

(2)目前众多研究采用的特征选取算法有pearson相关分析、mRMR(最大相关最小冗余)等特征选取算法。但是，上述算法均是基于与刀具状态的线性/非线性相关度来考察，而相关性强的特征参数并不一定能获得良好的监测性能，特别是在切削加工这一复杂的过程中。本发明从预测精度出发，以预测精度和特征参数个数为优化目标，考察各特征参数组合对应的刀具磨损量预测精度，避免了相关性强但预测精度低的现象。

附图说明

图1为本发明的方法流程图；

图2为本发明的整体流程运行图。

具体实施方式

下面结合附图对本发明做进一步详细描述：

如图1和图2所示，为本发明实施例中，提出的一种铣削刀具磨损状态特征提取多目标优化方法，包括以下步骤：

S1、通过多种传感器通道采集多个物理场信号。其主要包含：

针对新切削刀具，周期性采集刀具运行状态下的S个物理场(如振动、电流、声发射、切削力、声音等)传感通道的时域信号以及对应的刀具磨损量，共分别采集T次，组成训练样本集(X，Y)＝{(X_t,Y_t)}，X_t∈R^N×T，Y_t∈R，t＝1,2,...,T，X_t＝(x_1t,…,x_St)表示第t次采集信号时所有S个传感通道的信号，Y_t表示第t次采集信号时对应的刀具磨损量，t表示第t次采集的信号，N表示每个传感通道每次采样的信号点数，其中M和N均为大于1的正整数。

S2、计算每个传感通道信号的若干时域和频域统计特征参数、以及小波能量，构成特征参数候选集。其主要包含：

计算每个传感通道x_it(i＝1,…,M)的F个统计特征参数(可以包含时域、频域、时频域小波能量等统计参数)，组成x_it的特征参数集f_it＝(f_it-1,f_it-2,...,f_it-F)，并将所有S个通道的特征参数集合并，构造特征参数候选集F＝{f_it}。为描述方便，将X_t对应的候选参数集F记为F_t＝{f_tk,t＝1,…,T,k＝1,…,Q}，Q为候选参数总个数。

S3、以刀具状态预测精度和特征参数数量为优化目标，构造多目标优化模型。其主要包含：构造多目标优化模型

式(1)、(2)和(3)中，Y′_t表示在某种特征参数组合F_δq下通过预测算法得到的X_t的预测值，预测算法如下：

a)找出{δ_q}(q＝1,…,Q)对应特征参数组合下的样本数据集，即选择{δ_q}中值为1的位置对应的特征参数。例如，{δ_q}＝{1,0,0,1,1,0......,}表示选择候选集F中第1、4、5个参数，进而提炼出样本数据集(X_δq,Y)＝{(f_δq-t,Y_t)}。

b)计算样本核矩阵

其中

为X_δq中两个样本点x_i和x_j的核函数值，核函数k(x_i,x_j)可以取线性核函数、高斯核函数、多项式核函数等常用的核函数。

c)按照下式计算X_t的预测值Y′_t：

式(4)中，I表示单位矩阵，c表示正则化参数(一般取2)。

S4、采用智能优化算法对优化模型进行全局寻优。其主要包含：

采用多目标智能进化算法,如粒子群算法、遗传算法、蚁群算法等,求解上一步骤中优化模型的全局最优解，求得最优的特征参数集F*。

S5、将优化模型最优解所对应的传感特征参数集作为刀具状态监测所需的特征参数。其主要包含：

a)计算F*对应的样本核矩阵

b)在线监测，对于未知信号x，其对应的刀具磨损量预测值计算如下：

以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (1)

1.一种铣削刀具磨损状态特征提取多目标优化方法，其特征在于，包括以下步骤：

S1、通过多种传感器通道采集多个物理场信号；

S2、计算每个传感通道信号的若干时域和频域统计特征参数、以及小波能量，构成特征参数候选集；

S3、以刀具状态预测精度和特征参数数量为优化目标，构造多目标优化模型；

S4、采用智能优化算法对优化模型进行全局寻优；

S5、将优化模型最优解所对应的传感特征参数集作为刀具状态监测所需的特征参数；

所述的步骤S1具体包括：

周期性采集刀具运行状态下的S个物理场传感通道的时域信号以及对应的刀具磨损量，共分别采集T次，组成训练样本集。

所述的步骤S2具体包括：

S21、计算每个传感通道x_it(i＝1,…,M)的F个统计特征参数，组成特征参数集f_it＝(f_it-1,f_it-2,...,f_it-F)；

S22、将所有S个传感通道的特征参数集合并，构造特征参数候选集F＝{f_it}；

所述的步骤S3具体包括：构造多目标优化模型

其中，Y_t′表示在某种特征参数组合F_δq下通过预测算法得到的X_t的预测值；

所述的预测算法包括如下步骤：

S301、找出{δ_q}(q＝1,…,Q)对应特征参数组合下的样本数据集，即选择{δ_q}中值为1的位置对应的特征参数；

S302、计算样本核矩阵

其中

为X_δq中两个样本点x_i和x_j的核函数值；

S303、按照下式计算X_t的预测值Y_t′：

其中，I表示单位矩阵，c表示正则化参数；

所述的步骤S4具体包括：

所述的智能优化算法包括粒子群算法、遗传算法和蚁群算法，并求得最优的特征参数集F*；

所述的步骤S5具体包括：

S51、计算F*对应的样本核矩阵

S52、在线监测，对于未知信号x，其对应的刀具磨损量预测值计算如下：

Images