CN108920812A - 一种机械加工表面粗糙度预测方法 - Google Patents

Abstract

本发明提供一种机械加工表面粗糙度的预测方法，包括以下步骤：选取足够数量的检测参数和对应的表面粗糙度作为样本数据；对样本数据进行归一化处理；采用试凑法确定网络结构；初始化Copula EDA的群种，通过计算预测值和实测值的误差选择S个优秀个体；通过Copula分布估计算法来优化网络参数，将最优的权值和阈值赋给神经网络；最后基于BP算法进一步修正网络的参数，直至达到期望误差或者迭代次数，保存最终的网络来计算表面粗糙度的预测值。本方法将Copula EDA和BP算法相结合来优化神经网络，既解决了单一BP算法容易陷入局部极小的情况，又解决了单一Copula EDA算法收敛速度慢的缺点，通过在线监测运行状态下的工艺系统，能够实现机械表面加工质量的预测预报。

Description

一种机械加工表面粗糙度预测方法

技术领域

本发明属于机械加工技术领域，尤其是涉及一种机械加工表面粗糙度预测方法。

背景技术

智能加工作为智能制造的关键技术，需要对机械加工过程实时在线监测，采集切削力、振动、声发射、温度等检测信号，提取相关特征，对切削过程的状态进行判定，并能对加工质量、刀具耐用度和工艺系统的稳定性进行预测，进而调整切削参数来实现智能控制。在机械加工中，表面粗糙度是衡量产品表面质量的重要指标之一。因此，如何有效地预测工件的表面粗糙度成为了机械制造领域的一个关键问题。

目前，国内外许多学者对表面粗糙度预测技术进行了深入研究。目前关于表面粗糙度预测研究方法主要分为三种：基于切削理论模型、基于切削参数的经验模型以及基于人工智能的预测方法。基于切削理论模型和基于切削参数的经验模型的预测方法，由于误差较大，更多的是采用人工智能的预测方法。BP神经网络具有理论上逼近任何非线性连续函数的能力，因而其在研究中得到广泛的应用，但是BP算法容易陷入局部极小值，无法达到全局最优点，而Copula分布估计算法(Copula Estimation of Distribution Algorithm，Copula EDA)则是基于全局来寻找最优解，因此很好的弥补了BP算法的劣势。

发明内容

针对智能制造中通过在线监测切削加工过程来判定系统的状态及预测，提出了基于Copula EDA和BP神经网络混合算法的表面粗糙的预测新方法，以监测信号作为变量参数，以表面粗糙度作为输出，选取足够数量的变量参数和对应的表面粗糙度作为样本数据；然后，对样本数据进行归一化处理；进而，基于公式，采用试凑法确定网络结构；然后，初始化Copula EDA的群种，通过计算预测值和实测值的误差选择S个优秀个体；通过Copula分布估计算法来优化网络参数，将最优的权值和阈值赋给神经网络；最后基于BP算法进一步修正网络的参数，直至达到期望误差或者迭代次数，保存最终的网络来计算表面粗糙度的预测值。

本发明通过以下技术方案实现：一种机械加工表面粗糙度的预测方法，包括以下步骤：

步骤S1、样本数据的获取：选取多组检测参数和对应的表面粗糙度作为样本数据；

步骤S2、样本数据预处理：将步骤S1中的检测参数和表面粗糙度数据进行归一化，并将归一化后的参数和表面粗糙度分别作为神经网络的输入和输出样本；

步骤S3、确定神经网络结构：利用试凑法确定确定网络结构；

步骤S4、初始化群种：随机产生P个神经网络的权值和阈值；

步骤S5、选择S个优秀个体：将S4产生的P个神经网络的权值和阈值逐一赋给神经网络，计算表面粗糙度R_a预测值和实测值的误差，选择S个优秀个体；

步骤S6、基于Copula分布估计算法优化网络参数：利用Copula分布估计算法优化网络参数，得到最优的权值和阈值，并将最优的权值和阈值赋给神经网络；

步骤S7、基于BP算法优化网络参数：将S6得到的最优的权值和阈值赋给神经网络，直至达到期望误差或者迭代次数，得到最终的最优权值和阈值来计算表面粗糙度的预测值。

上述方案中，所述步骤S3具体为：根据输入层和输出层单元格数来确定隐含层的单元个数，基于公式其中m为输入层单元个数，n为输出层单元个数，利用试凑法确定误差最小时所对应的节点数。

上述方案中，所述步骤S4具体为：随机产生P个神经网络的权值和阈值{x_i1,x_i2,…x_ij}，i＝1,2,…,P，j是权值和阈值按照编码方式排列成长度为k(m+n)+n+k的一维数组的字符数，表示权值和阈值初始群种的维数，其中m为输入层单元个数，n为输出层单元个数，k为隐含层单元个数。

上述方案中，所述步骤S5具体为：将步骤S4产生的P个权值和阈值逐一赋给神经网络，并用归一化的切削力F作为输入得到的表面粗糙度R_a预测值，以实测值和预测值误差大小来评价这P个权值和阈值优劣，从其中选择误差最小的r个个体以及随机选择t个个体作为Copula EDA的S个优秀个体，S＝r+t。

上述方案中，所述步骤S6具体为：建立所述S个优秀个体{x_ij,i＝1,2,…,S}的概率分布模型F_j，利用Copula函数进行采样，产生R个服从联合分布函数C的变量u₁ ^k,u₂ ^k,…u_j ^k，表示第k组j维变量，k＝1,2,…,R,根据公式x_j ^k＝F_j ^-1(u_j ^k)，得到R个新个体x₁ ^k,x₂ ^k,…x_j ^k，同时从S个优秀个体中选择d个个体，并随机产生P-R-d个个体，将这三组组合起来作为下一代初始群种，返回步骤S5，直到达到期望误差或者迭代次数，停止优化，选出误差最小所对应的个体解码得到最优的权值和阈值。

上述方案中，所述步骤S7具体为：将步骤S6得到的最优的权值和阈值赋给神经网络，并在训练样本数据的过程中利用BP算法不断的修正权值和阈值，达到期望误差或者迭代次数停止训练，得到最终的最优权值和阈值。

上述方案中，还包括步骤S8；所述步骤S8为评价网络，具体为：将步骤S7得到的最终权值和阈值赋给网络，以验证实验的切削力F为输入，并计算预测的表面粗糙度R_a和实测值的误差，以此来评价网络结构的准确性和可靠性。

上述方案中，所述优秀个体S的个数不少于P/2个

上述方案中，选用正态分布作为优秀个体的边缘分布，选用Clayton Copula函数进行采样。

上述方案中，所述步骤6和步骤7在优化神经网络参数过程中，以表面粗糙度误差率最小为作为选择最优权值和阈值的评价标准。

与现有技术相比，本发明的有益效果是：本方法将Copula EDA和BP算法相结合来优化神经网络，既解决了单一BP算法容易陷入局部极小的情况，又解决了单一Copula EDA算法收敛速度慢的缺点，为自动化加工和机械加工过程在线监测提供了一个新的视角。通过在线监测运行状态下的工艺系统，能够实现机械表面加工质量的预测预报，对于实现智能制造和工业4.0具有重要的科学意义和较大的实用价值。

附图说明

图1是本发明一实施方式的Copula EDA-BP混合算法流程图；

图2是本发明一实施方式的1-4-1的神经网络模型。

具体实施方式

下面结合附图1以及具体实施例对本发明作进一步的说明，但本发明的保护范围并不限于此。

如图1所示，本发明所述一种机械加工表面粗糙度预测方法：包括以下步骤：

步骤S1、样本数据的获取；

步骤S2、样本数据预处理：

将实验结果中的检测参数和表面粗糙度数据进行归一化，并将归一化后的参数和表面粗糙度分别作为神经网络的输入和输出样本；

步骤S3、确定网络结构：

根据输入层和输出层单元格数来确定隐含层的单元个数，本文输入和输出各有一个单元，根据其中m为输入层单元个数，n为输出层单元个数，利用试凑法确定误差最小时所对应的节点数为最优；

步骤S4、初始化群种：

随机产生P个神经网络的权值和阈值{x_i1,x_i2,…x_ij}，i＝1,2,…,P，j是权值和阈值按照编码方式排列成长度为k(m+n)+n+k的一维数组的字符数，表示权值和阈值初始群种的维数，其中m为输入层单元个数，n为输出层单元个数，k为隐含层单元个数；

步骤S5、选择S个优秀个体：

将步骤S4产生的P个权值和阈值逐一赋给神经网络，并用归一化的切削力F作为输入得到的表面粗糙度R_a预测值，以实测值和预测值误差大小来评价这P个权值和阈值优劣，从其中选择误差最小的r个个体以及随机选择t个个体作为Copula EDA的S个优秀个体，S＝r+t；

步骤S6、基于Copula分布估计算法优化网络参数：

建立这S个优秀个体{x_ij,i＝1,2,…,S}的概率分布模型F_j，利用Copula函数进行采样，产生R个服从联合分布函数C的变量(u₁ ^k,u₂ ^k,…u_j ^k)，表示第k组j维变量，k＝1,2,…,R,根据公式x_j ^k＝F_j ^-1(u_j ^k)，得到R个新个体(x₁ ^k,x₂ ^k,…x_j ^k)，同时从S个优秀个体中选择d个个体，并随机产生P-R-d个个体，将这三组组合起来作为下一代初始群种，返回步骤S5。直到达到期望误差或者迭代次数，停止优化，选出误差最小所对应的个体解码得到权值和阈值；

步骤S7、基于BP算法优化网络参数：

将步骤S6得到的最优的权值和阈值赋给神经网络，并在训练样本数据的过程中利用BP算法不断的修正权值和阈值，达到期望误差或者迭代次数停止训练，最终得到最优的权值和阈值；

步骤S8、评价网络：

将步骤S7得到的最终权值和阈值赋给网络，以验证实验的切削力F为输入，并计算预测的表面粗糙度R_a和实测值的误差，以此来评价网络结构的准确性和可靠性。

具体实施例：

在铣削45钢试验中，本实例选用主切削力为检测参数。

S1、样本数据的获取

为了较好地描述切削力与表面粗糙度两变量的相关性，选取不同切削条件下所对应的切削力和表面粗糙度值作为样本数据，如表1所示。

其中，v代表切削速度，a_f代表切削速度每齿进给量，a_p代表切削深度，F代表切削力，R_a代表二维表面粗糙度。

表1铣削试验结果

S2、样本数据预处理

本实例采用的数据预处理方法是对数据进行归一化，它能够提高网络的收敛速度，同时还能避免数值问题，而且在BP算法中常采用sigmoid函数作为转移函数，归一化能够防止净输入绝对值过大而引起的神经元输入饱和现象。将S1中的检测参数和表面粗糙度数据进行归一化，并将归一化后的参数和表面粗糙度分别作为神经网络的输入和输出样本。本实施例使用的归一化公式为：

式一中x'是经过归一化公式变换后的数据，x是原始数据，x_max和x_min分别为原始数据中的最大值和最小值。

S3、确定神经网络结构

神经网络中各层节点数的选择对网络的性能影响很大，所以，需要对各层节点数进行优选。输入层和输出层的节点数目取决于实际问题的数据源，隐含层神经元个数的确定是一个十分复杂的问题，与具体的问题、输入和输出层神经元个数有着重要的关系，大大影响了神经网络的预测精度。如果隐含层神经元个数选取太少时，网络不能充分学习，局部极小的状况就多，学习过程可能不收敛；如果隐含层神经元个数选择太多时，造成学习时间太长，误差也比较大，网络容易出现训练过度的现象。所以隐含层神经元个数应该存在一个最佳的数目，但是目前并没有很好的确定方法，根据输入层和输出层单元格数来确定隐含层的单元个数，本实施例中输入和输出各有一个单元，本实例根据其中m为输入层单元个数，n为输出层单元个数，即隐含层神经元个数取值范围为[3,11]，使用相同的样本集进行训练，采用试凑法来计算最合适的节点数，使得网络误差最小的神经网络就是所需的隐含层节点个数，最终确定隐含层单元个数为4。

S4、初始化群种

随机产生P＝50个权值和阈值{x_i1,x_i2,…x_ij}，i＝1,2,…50，j是权值和阈值按照实数编码方式排列成长度为k(m+n)+n+k的一维数组的字符数，表示权值和阈值初始群种的维数，其中m为输入层单元个数，n为输出层单元个数，k为隐含层单元个数。本实施例的编码长度j＝4×(1+1)+1+4＝13。

图2是本实例1-4-1的神经网络模型其中1-4-1代表输入层有一个单元数，隐含层有4个单元数，输出层有一个单元数。隐含层到输入层的权值矩阵为阈值矩阵为输出层到隐含层的权值矩阵为阈值矩阵为B₂＝[b₆]，则神经网络权值和阈值编码后的形式为：

character_i＝[w₂₁w₃₁w₄₁w₅₁w₆₂w₆₃w₆₄w₆₅b₂b₃b₄b₅b₆] 式二

式二中i表示第i组权值和阈值经编码后的字符串。

S5、选择S个优秀个体

将随机产生的P＝50个权值和阈值一一赋给神经网络，将归一化处理后的64组样本数据中的切削力值作为输入，得到64个表面粗糙度的预测值，反归一化处理后计算与样本数据中的表面粗糙度值的误差，选出误差最小的r＝15个字符串，再随机选择t＝10个字符串作为Copula EDA的S(S＝r+t＝25)个优秀个体。

S6、基于Copula分布估计算法优化网络参数

建立这S个优秀个体{x_ij,i＝1,2,…,25,j＝1,2,…,13}的概率分布模型F_j，选用正态分布作为其边缘分布，即X_i～(mean(x_i)，std(x_i))，利用Copula函数，本实例选用Clayton Copula进行采样，产生R＝25个服从联合分布函数C的变量(u₁ ^k,u₂ ^k,…u_j ^k)，k＝1,2,…,25,根据公式x_j ^k＝F_j ^-1(u_j ^k)，得到25个新个体(x₁ ^k,x₂ ^k,…x_j ^k)，同时从S＝25个优秀个体中选择d＝10个个体，并随机产生P-R-d＝50-25-10＝15个个体，将这三组组合起来作为下一代初始群种，返回S5。直到达到期望误差或者迭代次数，停止优化，选出误差最小所对应的个体解码得到权值和阈值。

S7、基于BP算法优化网络参数

设置网络的迭代次数(net.trainParam.epochs＝10000)，网络的迭代步长(net.trainParam.show＝50)，网络的训练目标(net.trainParam.goal＝0.001)，相网络的学习系数(net.trainParam.lr＝0.05)，网络的动量因子(net.trainParam.mc＝0.9)，将S6得到的最优的权值和阈值赋给神经网络，并在训练样本数据的过程中利用BP算法不断的修正权值和阈值，达到期望误差或者迭代次数停止训练，最终得到最优的权值和阈值。输入层到隐含层的权值为[5.6000；-5.6000；5.6000；5.6000]，阈值[-5.6000；1.8667；1.8667；5.6000]，隐含层到输出层的权值为[0.3575；0.5155；0.4863；-0.2155]，阈值为0.3110。

S8、评价网络

保存经过Copula EDA-BP混合算法优化后的神经网络，将验证实验中的12组切削力F归一化处理后作为输入，得到12组表面粗糙度的预测值Ra_pre，并计算其与实测值的误差，结果如表2所示。

表2表面粗糙度预测模型的分析结果

从表2中可以看出，该预测模型的平均准确率是91.98％，所有预测数据的准确率均高于80％，低误差率说明了表面粗糙度预测值Ra_pre与测量值Ra的拟合情况良好。

所述实施例为本发明的优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。

Claims (10)

1.一种机械加工表面粗糙度的预测方法，其特征在于，包括以下步骤：

步骤S1、样本数据的获取：选取多组检测参数和对应的表面粗糙度作为样本数据；

步骤S2、样本数据预处理：将步骤S1中的检测参数和表面粗糙度数据进行归一化，并将归一化后的参数和表面粗糙度分别作为神经网络的输入和输出样本；

步骤S3、确定神经网络结构：利用试凑法确定确定网络结构；

步骤S4、初始化群种：随机产生P个神经网络的权值和阈值；

步骤S5、选择S个优秀个体：将S4产生的P个神经网络的权值和阈值逐一赋给神经网络，计算表面粗糙度R_a预测值和实测值的误差，选择S个优秀个体；

步骤S6、基于Copula分布估计算法优化网络参数：利用Copula分布估计算法优化网络参数，得到最优的权值和阈值，并将最优的权值和阈值赋给神经网络；

步骤S7、基于BP算法优化网络参数：将S6得到的最优的权值和阈值赋给神经网络，直至达到期望误差或者迭代次数，得到最终的最优权值和阈值来计算表面粗糙度的预测值。

2.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，所述步骤S3具体为：根据输入层和输出层单元格数来确定隐含层的单元个数，基于公式其中m为输入层单元个数，n为输出层单元个数，利用试凑法确定误差最小时所对应的节点数。

3.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，所述步骤S4具体为：随机产生P个神经网络的权值和阈值{x_i1,x_i2,…x_ij}，i＝1,2,…,P，j是权值和阈值按照编码方式排列成长度为k(m+n)+n+k的一维数组的字符数，表示权值和阈值初始群种的维数，其中m为输入层单元个数，n为输出层单元个数，k为隐含层单元个数。

4.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，所述步骤S5具体为：将步骤S4产生的P个权值和阈值逐一赋给神经网络，并用归一化的切削力F作为输入得到的表面粗糙度R_a预测值，以实测值和预测值误差大小来评价这P个权值和阈值优劣，从其中选择误差最小的r个个体以及随机选择t个个体作为Copula EDA的S个优秀个体，S＝r+t。

5.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，所述步骤S6具体为：建立所述S个优秀个体{x_ij,i＝1,2,…,S}的概率分布模型F_j，利用Copula函数进行采样，产生R个服从联合分布函数C的变量u₁ ^k,u₂ ^k,…u_j ^k，表示第k组j维变量，k＝1,2,…,R,根据公式x_j ^k＝F_j ^-1(u_j ^k)，得到R个新个体x₁ ^k,x₂ ^k,…x_j ^k，同时从S个优秀个体中选择d个个体，并随机产生P-R-d个个体，将这三组组合起来作为下一代初始群种，返回步骤S5，直到达到期望误差或者迭代次数，停止优化，选出误差最小所对应的个体解码得到最优的权值和阈值。

6.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，所述步骤S7具体为：将步骤S6得到的最优的权值和阈值赋给神经网络，并在训练样本数据的过程中利用BP算法不断的修正权值和阈值，达到期望误差或者迭代次数停止训练，得到最终的最优权值和阈值。

7.根据权利要求1所述的机械加工表面粗糙度的预测方法，其特征在于，还包括步骤S8；所述步骤S8为评价网络，具体为：将步骤S7得到的最终权值和阈值赋给网络，以验证实验的切削力F为输入，并计算预测的表面粗糙度R_a和实测值的误差，以此来评价网络结构的准确性和可靠性。

8.根据权利要求1所述的表面粗糙度预测方法，其特征在于，所述优秀个体S的个数不少于P/2个。

9.根据权利要求1所述的表面粗糙度预测方法，其特征在于，选用正态分布作为优秀个体的边缘分布，选用Clayton Copula函数进行采样。

10.根据权利要求1所述的表面粗糙度预测方法，其特征在于，所述步骤6和步骤7在优化神经网络参数过程中，以表面粗糙度误差率最小为作为选择最优权值和阈值的评价标准。