一种CFRP铣削工艺参数的优化方法及系统

一种cfrp铣削工艺参数的优化方法及系统
技术领域
1.本发明涉及机械加工领域，特别是涉及一种cfrp铣削工艺参数的优化方法及系统。


背景技术：

2.碳纤维增强复合材料(carbon fiber reinforced polymer，cfrp)以其优异的综合性能，现已经被普遍的应用在航空航天、风电、运动器材、汽车等多个领域。cfrp构件一般采用一次整体成型的制造工艺，在成型后，为了保证设计要求的制造精度和装配精度等，通常需要对复合材料零件进行二次加工。由于cfrp的层合结构、非均匀分布、各向异性等特点，铣削过程中易产生毛刺、撕裂、分层等加工缺陷，加工质量难以控制。因此cfrp铣削工艺参数优化十分重要，这可以在一定程度上降低生产成本和废品率。


技术实现要素：

3.有鉴于此，本发明的目的在于提供一种cfrp铣削工艺参数的优化方法及系统，实现快速有效的cfrp铣削工艺参数优化，有效降低生产成本和废品率。
4.为实现上述目的，本发明采用如下技术方案：
5.一种碳纤维增强复合材料(carbon fiber reinforced polymer，cfrp)铣削工艺参数的优化方法，其特征在于，包括以下步骤：
6.步骤s1:根据铣削碳纤维材料的表面质量为目标，并基于预设优化变量及约束，构建实验方案；
7.步骤s2：根据实验结果获取训练集和测试集，并利用灰狼算法(grey wolf optimization，gwo)改进支持向量回归机(support vector machine，svr)svr算法进行拟合目标，得到gsvr模型；
8.步骤s3:基于gwo算法建立cfrp铣削工艺参数优化模型，将gsvr模型作为gwo算法的适应度函数，得到gwo-gsvr模型，并基于gwo-gsvr模型对cfrp铣削工艺参数寻优，得到最优表面质量下的工艺参数。
9.进一步的，所述优化变量包括切主轴转速nz、进给速度vj、轴向切深a
p
和径向切深ae；其中轴向切深a
p
和径向切深ae设为定值，设计主轴转速nz和进给速度vj不同水平做全因子实验。
10.进一步的，所述约束具体为：
11.主轴转速n
min
≤nz≤n
max
12.进给速度v
min
≤vj≤v
max
13.其中，n
min
、n
max
是主轴转速nz的最小值与最大值；v
min
、v
max
是进给速度vj的最小值与最大值。
14.进一步的，所述步骤s2具体为：
15.步骤s21:根据实验结果整理数据集为d＝{(x1,y1),(x2,y2),...(xi,yi)...,(xm,
ym)}，按照预设的比例进行划分数据集，划分为训练集d
train
和测试集d
test
；
16.步骤s22:对数据进行归一化处理；
17.步骤s23:利用gwo算法对svr算法进行参数寻优，得到gsvr模型。
18.进一步的，所述步骤s22具体为：原自变量为x＝{x1,x2,...xi...,xm},归一化方法其中变换后为d＝{(x
g1
,y1),(x
g2
,y2),...(x
gi
,yi)...,(x
gm
,ym)}。
19.其中，xi为变量初始值；x
gi
为归一化之后的变量；max(x)为原变量中的最大值；min(x)为原变量中的最小值。
20.进一步的，所述步骤s23具体为：
21.(1)构建svr算法目标形式：f(x)＝a
t
·
x b。
22.其中，f(x)是回归值，是cfrp三维表面粗糙度的表示；x是输入数据向量，是表示为主轴转速nz、进给速度vj、轴向切深a
p
、径向切深ae组合的向量；a是权重向量，b为截距向量；
23.(2)引入松弛因子，将问题转化为求解最小化目标函数：
[0024][0025]
约束为：
[0026][0027]
其中，ξi,为松弛变量；c为惩罚函数；ε为损失界限；
[0028]
(3)引入拉格朗日乘子法，利用对偶变换和非线性变换将目标函数进行转换；由于模型问题属于非线性问题，将通过核函数将自变量x映射到高维特征空间，从而把非线性回归问题转化为线性回归问题；
[0029]
(4)对svr算法的拟合性能影响最大的是惩罚系数c、核系数γ和损失界限ε，并利用灰狼算法进行智能寻优，利用最有参数进行拟合模型。
[0030]
进一步的，所述根据灰狼算法进行智能寻优，具体为：
[0031]
(1)在三维空间面进行初始化种群，设有n个灰狼个体，并且设置迭代次数为t，初始化个体灰狼为令第一个维度为惩罚系数c的大小，第二个维度为核函数系数γ大小，第三个维度为损失界限ε大小，且种群生成方法为佳点集生成法：
[0032]
pi＝[{d1*i},{d2*i},{d3*i}],i＝1,2,...,n
[0033][0034][0035][0036]
其中pi表示编号为i的灰狼对应的佳点集合；{dj*i}表示灰狼i在第j维度上对应的佳点元素，且在数值上取dj*i小数部分；p表示满足(p-3)/2≥3约束内的最小素数；表示
第i只灰狼在第j个维度上的数值；表示佳点系数，对应于佳点集合内第j维的元素；n表示灰狼种群的个体数；t表示需要迭代的最大次数；lb1,ub1表示惩罚系数c约束上的最小值和最大值；lb2,ub2表示核函数系数γ约束上的最小值和最大值；lb3,ub3表示损失界限ε约束上的最小值和最大值；
[0037]
(2)将初始化的灰狼种群个体带入svr算法进行拟合，得出拟合预测值y
p
＝{y
p1
,y
p2
,...y
pn
}，将预测试与实际值进行计算且灰狼的适应度值为ys＝mse；
[0038]
其中，y
p
代表灰狼群的预测值；y
pi
代表灰狼的预测值；yi代表灰狼的真实值；ys代表灰狼的适应度值；
[0039]
(3)对灰狼个体进行适应度排序，其中适应度最小为α灰狼，适应度第二小的为β灰狼，适应度第三小为δ灰狼，剩下灰狼称为的称为ω灰狼；
[0040]
(4)狼群进行包围猎物，灰狼靠近猎物；
[0041]
d＝|b
·
x
p
(t)-x(t)|
[0042]
x(t 1)＝x
p
(t)-a
·d[0043]
a＝2a
·r1-a
[0044]
b＝2
·
r2[0045][0046]
其中，d代表灰狼个体与猎物之间的距离；t表示目前的迭代次数；a和b表示系数向量；x
p
表示猎物位置；x表示当前灰狼的位置；a式收敛因子，随着迭代次数的增加由2非线性减小到0；r1和r2是随机向量，且它的模是[0,1]之间的随机数；
[0047]
(5)狩猎，灰狼进行识别猎物的位置，进行包围猎物，并且假设α、β和δ更加了解猎物的潜在位置，我们在迭代的过程中，保留这三个最优解，进行更新其他灰狼的个体位置，其中灰狼个体追踪猎物的模型如下：
[0048]dα
＝|b1·
x
α-x|
[0049]dβ
＝|b2·
x
β-x|
[0050]dδ
＝|b3·
x
δ-x|
[0051]
其中，d
α
、d
β
和d
δ
是表示α、β和δ灰狼与其他灰狼之间的距离；x
α
、x
β
和x
δ
表示α、β和δ灰狼的当前位置；b1、b2和b3是随机向量；
[0052]
(6)灰狼进行狩猎，ω灰狼群向猎物靠近
[0053]
x1＝x
α-a
·dα
[0054]
x2＝x
β-a
·dβ
[0055]
x3＝x
δ-a
·dδ
[0056][0057]
其中，x1、x2、x3表示ω灰狼朝向α、β、δ灰狼前进的步长和方向；x(t 1)代表的时ω灰狼更新后的位置；
[0058]
(7)进行迭代，当不满足迭代条件的时候，按照(3)-(6)进行迭代处理；当满足迭代条件时候，输出α灰狼，得到优化目标最优解；其中α灰狼的位置坐标表示着svr的惩罚系数c、核系数γ和损失界限ε；
[0059]
(8)将优化的惩罚系数c、核系数γ和损失界限ε代入svr算法，得到gsvr拟合模型，将训练集d
train
和测试集d
test
代入gsvr进行交叉训练，训练结束保留拟合模型。
[0060]
进一步的，所述核函数类型为高斯径向基核函数；具体回归形式转换为：
[0061][0062]
k(xi,x)＝exp(-γ||x
i-x||2)
[0063][0064]
其中，γ为核系数，且δ表示高斯核的带宽；αi,为拉格朗日乘子；k(xi,x)为高斯径向基核函数。
[0065]
进一步的，所述步骤s3具体为：
[0066]
步骤s31:灰狼种群初始化，设有no个灰狼个体，并且设置迭代次数为to，令第一个维度为主轴转速nz的大小，第二个维度为进给速度vj的大小，第三个维度为轴向切深a
p
的大小，第四个维度为径向切深ae的大小，初始化个体灰狼为且轴向切深a
p
和径向切深ae是定值，在二维平面利用为佳点集生成法生成初始化种群：
[0067]
p
oi
＝[{d
o1
*i},{d
o2
*i}],i＝1,2,...,no[0068][0069][0070][0071]
其中p
oi
表示编号为i的灰狼对应的佳点集合；{d
oj
*i}表示灰狼i在第j维度上对应的佳点元素，且在数值上取d
oj
*i小数部分；po表示满足(p
o-2)/2≥2约束内的最小素数；表示第i只灰狼在第j个维度上的数值；表示佳点系数，对应于佳点集合内第j维的元素；no表示灰狼种群的个体数；to表示需要迭代的最大次数；表示主轴转速nz约束下的最小值和最大值；表示进给速度vj约束下的最小值和最大值；表示轴向切深a
p
和径向切深ae的大小，是提前预设的定值。
[0072]
步骤s32:将初始化的灰狼种群个体gsvr模型中，得出拟合预测值
[0073]
其中，y
p
代表灰狼种群的预测值；y
pi
代表灰狼的预测值，此值为三维表面粗糙度，且作为灰狼的适应度值大小。
[0074]
步骤s33:进行灰狼种群等级分层，其中适应度最小为αo灰狼，适应度第二小的为βo灰狼，适应度第三小为δo灰狼，剩下灰狼称为的称为ωo灰狼；
[0075]
步骤s34:狼群进行包围猎物，灰狼靠近猎物；
[0076]
步骤s35:灰狼种群进行狩猎，并依据αo、βo和δo灰狼个体的位置更新ωo灰狼个体的位置。
[0077]
步骤s36:进行迭代，在满足迭代条件时，输出迭代结果αo灰狼，得到优化目标最优解；其中αo灰狼的位置坐标表示着cfrp铣削表面粗糙度最好时的工艺参数主轴转速nz、进给速度vj、轴向切深a
p
、径向切深ae。
[0078]
进一步的，所述灰狼社会的等级分类，具体为：
[0079]
第一层为：灰狼狼群的头狼记为α灰狼，α灰狼可以支配其他所有灰狼；
[0080]
第二层为：灰狼狼群中的β灰狼，其受α灰狼的支配，但可以支配狼群社会其他阶层的灰狼；
[0081]
第三层为：灰狼狼群中的δ灰狼，其受α灰狼和β灰狼的支配，但可以支配狼群社会其他阶层的灰狼；
[0082]
第四层为：灰狼狼群中的ω灰狼，其受α灰狼、β灰狼和δ灰狼的支配。
[0083]
一种cfrp铣削工艺参数的优化系统,其特征在于，包括处理器、存储器以及存储在所述存储器上的计算机程序，所述处理器执行所述计算机程序时，具体执行如上所述的铣削工艺参数的优化方法中的步骤。
[0084]
本发明与现有技术相比具有以下有益效果：
[0085]
本发明实现快速有效的cfrp铣削工艺参数优化，有效降低生产成本和废品率。
附图说明
[0086]
图1是本发明方法流程图。
具体实施方式
[0087]
下面结合附图及实施例对本发明做进一步说明。
[0088]
参考图1，在本实施例中，提供一种cfrp铣削工艺参数的优化方法，包括以下步骤：
[0089]
步骤s1:根据铣削碳纤维材料的表面质量为目标，并基于预设优化变量及约束，构建实验方案；
[0090]
步骤s2：根据实验结果获取训练集和测试集，并利用gwo算法改进svr算法进行拟合目标，得到gsvr模型；
[0091]
步骤s3:基于gwo算法建立cfrp铣削工艺参数优化模型，将gsvr模型作为gwo算法的适应度函数，得到gwo-gsvr模型，并基于gwo-gsvr模型对cfrp铣削工艺参数寻优，得到最优表面质量下的工艺参数。
[0092]
在本实施例中，步骤s1具体为：
[0093]
1.1 选择实验因素与各实验因素水平，设计实验。考虑到cfrp铣削加工中主要是精加工，轴向切深a
p
和径向切深ae在实际加工往往已经给定，所以在实验中将轴向切深和径向切深设为定值，设计主轴转速nz和进给速度vj不同水平做全因子实验。
[0094]
设置相应约束：
[0095]
主轴转速n
min
≤nz≤n
max
[0096]
进给速度v
min
≤vj≤v
max
[0097]
其中，n
min
、n
max
是主轴转速nz的最小值与最大值；v
min
、v
max
是进给速度vj的最小值与最大值。
[0098]
1.2 测量铣削实验后的cfrp板材的表面质量，评估方式用三维表面粗糙度sa进行评价。
[0099]
在本实施例中，步骤s2，具体为：
[0100]
2.1 根据实验结果整理数据集为d＝{(x1,y1),(x2,y2),...(xi,yi)...,(xm,ym)}，按照7：3的比例进行划分数据集，划分为训练集d
train
和测试集d
test
。
[0101]
2.2 对数据进行归一化处理，原数据集的自变量进行归一化处理，原自变量为x＝{x1,x2,...xi...,xm},归一化方法其中变换后为d＝{(x
g1
,y1),(x
g2
,y2),...(x
gi
,yi)...,(x
gm
,ym)}。
[0102]
其中，
[0103]
xi为变量初始值；
[0104]
x
gi
为归一化之后的变量；
[0105]
max(x)为原变量中的最大值；
[0106]
min(x)为原变量中的最小值。
[0107]
2.3 cfrp铣削表面质量拟合模型
[0108]
利用gwo算法对svr算法进行参数寻优，得到gsvr模型，利用gsvr模型进行训练cfrp表面质量的拟合模型，具体步骤方法为：
[0109]
2.3.1 gsvr模型是利用gwo算法进行优化svr算法，并利用高斯径向基核函数进行映射高维空间，将非线性问题转化为线性问题；高斯径向基核函数是具有偏差小，学习性能好的优势，还可以简化计算程序，提高模型的泛化能力；gwo算法可以进行智能寻优，选择最优的一组参数。
[0110]
2.3.2 gsvr模型具体优化步骤
[0111]
2.3.2.1 构建svr算法目标形式：f(x)＝a
t
·
x b。
[0112]
其中，
[0113]
f(x)是回归值，是cfrp三维表面粗糙度的表示；
[0114]
x是输入数据向量，是表示为主轴转速nz、进给速度vj、轴向切深a
p
、径向切深ae组合的向量；
[0115]
a是权重向量，b为截距向量。
[0116]
2.3.2.2 为了保证svr算法具有良好的性能，引入松弛因子，将问题转化为求解最小化目标函数：
[0117][0118]
约束为：
[0119][0120]
其中，
[0121]
ξi,为松弛变量；
[0122]
c为惩罚函数；
[0123]
ε为损失界限。
[0124]
2.3.2.3 引入拉格朗日乘子法，利用对偶变换和非线性变换将目标函数进行转换。由于模型问题属于非线性问题，将通过核函数将自变量x映射到高维特征空间，从而把非线性回归问题转化为线性回归问题，并且选择核函数类型为高斯径向基核函数；具体回归形式转换为：
[0125][0126]
k(xi,x)＝exp(-γ||x
i-x||2)
[0127][0128]
其中，
[0129]
γ为核系数，且δ表示高斯核的带宽；
[0130]
αi,为拉格朗日乘子；
[0131]
k(xi,x)为高斯径向基核函数。
[0132]
2.3.2.4 对svr算法的拟合性能影响最大的是惩罚系数c、核系数γ和损失界限ε；传统方法是人工调参进行优化拟合模型，本发明利用灰狼算法进行智能寻优，利用最有参数进行拟合模型；因为灰狼算法具有结构简单、需要调节的参数少的优点，适合作为svr算法的参数寻优。
[0133]
2.3.2.5 在三维空间进行初始化种群，设有n个灰狼个体，并且设置迭代次数为t，初始化个体灰狼为令第一个维度为惩罚系数c的大小，第二个维度为核函数系数γ大小，第三个维度为损失界限ε大小，且种群生成方法为佳点集生成法，可以使初始种群由更好的遍历性，从而更好的达到全局寻优的目的。
[0134]
pi＝[{d1*i},{d2*i},{d3*i}],i＝1,2,...,n
[0135][0136][0137][0138]
其中
[0139]
pi表示编号为i的灰狼对应的佳点集合；
[0140]
{dj*i}表示第i只灰狼在第j个维度上对应的佳点元素，且在数值上取dj*i小数部分；
[0141]
p表示满足(p-3)/2≥3约束内的最小素数；
[0142]
表示第i只灰狼在第j个维度上的数值；
[0143]
表示佳点系数，对应于佳点集合内第j维的元素；
[0144]
n表示灰狼种群的个体数；t表示需要迭代的最大次数；
[0145]
lb1,ub1表示惩罚系数c约束上的最小值和最大值；
[0146]
lb2,ub2表示核函数系数γ约束上的最小值和最大值；
[0147]
lb3,ub3表示损失界限ε约束上的最小值和最大值。
[0148]
2.3.2.6 将初始化的灰狼种群个体带入svr算法进行拟合，得出拟合预测值y
p
＝{y
p1
,y
p2
,...y
pn
}，将预测试与实际值进行计算且灰狼的适应度值为ys＝mse。
[0149]
其中，
[0150]yp
代表灰狼群的预测值；
[0151]ypi
代表灰狼的预测值；
[0152]
yi代表灰狼的真实值；
[0153]ys
代表灰狼的适应度值。
[0154]
2.3.2.7 对灰狼个体进行适应度排序，其中适应度最小为α灰狼，适应度第二小的为β灰狼，适应度第三小为δ灰狼，剩下灰狼称为的称为ω灰狼。
[0155]
灰狼社会的等级分类：
[0156]
第一层为：灰狼狼群的头狼记为α灰狼，α灰狼可以支配其他所有灰狼；
[0157]
第二层为：灰狼狼群中的β灰狼，其受α灰狼的支配，但可以支配狼群社会其他阶层的灰狼。
[0158]
第三层为：灰狼狼群中的δ灰狼，其受α灰狼和β灰狼的支配，但可以支配狼群社会其他阶层的灰狼。
[0159]
第四层为：灰狼狼群中的ω灰狼，其受α灰狼、β灰狼和δ灰狼的支配。
[0160]
2.3.2.8 狼群进行包围猎物，灰狼靠近猎物。
[0161]
d＝|b
·
x
p
(t)-x(t)|
[0162]
x(t 1)＝x
p
(t)-a
·d[0163]
a＝2a
·r1-a
[0164]
b＝2
·
r2[0165][0166]
其中，
[0167]
d代表灰狼个体与猎物之间的距离；
[0168]
t表示目前的迭代次数；
[0169]
a和b表示系数向量；
[0170]
x
p
表示猎物位置；
[0171]
x表示当前灰狼的位置；
[0172]
a式收敛因子，随着迭代次数的增加由2非线性减小到0；
[0173]
r1和r2是随机向量，且它的模是[0,1]之间的随机数。
[0174]
2.3.2.9 狩猎，灰狼进行识别猎物的位置，进行包围猎物，并且假设α、β和δ更加了解猎物的潜在位置，我们在迭代的过程中，保留这三个最优解，进行更新其他灰狼的个体位
置。其中灰狼个体追踪猎物的模型如下：
[0175]dα
＝|b1·
x
α-x|
[0176]dβ
＝|b2·
x
β-x|
[0177]dδ
＝|b3·
x
δ-x|
[0178]
其中，
[0179]dα
、d
β
和d
δ
是表示α、β和δ灰狼与其他灰狼之间的距离；
[0180]
x
α
、x
β
和x
δ
表示α、β和δ灰狼的当前位置；
[0181]
b1、b2和b3是随机向量。
[0182]
2.3.2.10 灰狼进行狩猎，ω灰狼群向猎物靠近
[0183]
x1＝x
α-a
·dα
[0184]
x2＝x
β-a
·dβ
[0185]
x3＝x
δ-a
·dδ
[0186][0187]
其中，
[0188]
x1、x2、x3表示ω灰狼朝向α、β、δ灰狼前进的步长和方向；
[0189]
x(t 1)代表的时ω灰狼更新后的位置。
[0190]
2.3.2.11 进行迭代，当不满足迭代条件的时候，按照2.3.2.7-2.3.2.10进行迭代处理；当满足迭代条件时候，输出α灰狼，得到优化目标最优解。其中α灰狼的位置坐标表示着svr的惩罚系数c、核系数γ和损失界限ε。
[0191]
2.3.2.12 将优化的惩罚系数c、核系数γ和损失界限ε代入svr算法，得到gsvr拟合模型，将训练集d
train
和测试集d
test
代入gsvr进行交叉训练，训练结束保留拟合模型。
[0192]
在本实施例中，步骤s3，具体为：
[0193]
3.1 建立cfrp铣削工艺参数优化模型，优化算法也是利用gwo算法，gwo算法拥有者快的收敛速度和求解精度。
[0194]
3.2 重复步骤二中的改进gwo算法，主要包括：
[0195]
3.2.1 灰狼种群初始化，设有no个灰狼个体，并且设置迭代次数为to，令第一个维度为主轴转速nz的大小，第二个维度为进给速度vj的大小，第三个维度为轴向切深a
p
的大小，第四个维度为径向切深ae的大小，初始化个体灰狼为且轴向切深a
p
和径向切深ae是定值，在二维平面利用为佳点集生成法生成初始化种群，可以使初始种群由更好的遍历性，从而更好的达到全局寻优的目的。
[0196]
p
oi
＝[{d
o1
*i},{d
o2
*i},{d
o3
*i}],i＝1,2,...,no[0197][0198][0199][0200]
其中，
[0201]
p
oi
表示编号为i的灰狼对应的佳点集合；
[0202]
{d
oj
*i}表示灰狼i在第j维度上对应的佳点元素，且在数值上取d
oj
*i小数部分；
[0203]
po表示满足(p
o-2)/2≥2约束内的最小素数；
[0204]
表示第i只灰狼在第j个维度上的数值；
[0205]
表示佳点系数，对应于佳点集合内第j维的元素；
[0206]
no表示灰狼种群的个体数；to表示需要迭代的最大次数；
[0207]
表示主轴转速nz约束下的最小值和最大值；
[0208]
表示进给速度vj约束下的最小值和最大值；
[0209]
表示轴向切深a
p
和径向切深ae的大小，是提前预设的定值。
[0210]
3.2.2 将初始化的灰狼种群个体带入2.3.2.12训练的gsvr模型中，得出拟合预测值
[0211]
其中，
[0212]yp
代表灰狼种群的预测值；
[0213]ypi
代表灰狼的预测值，此值为三维表面粗糙度，且作为灰狼的适应度值大小。
[0214]
3.2.3 与2.3.2.7步骤一致，进行灰狼种群等级分层，其中适应度最小为αo灰狼，适应度第二小的为βo灰狼，适应度第三小为δo灰狼，剩下灰狼称为的称为ωo灰狼。
[0215]
3.2.4 与2.3.2.8步骤一致，狼群包围猎物，靠近猎物。
[0216]
3.2.5 与2.3.2.9-2.3.2.10步骤一致，灰狼种群进行狩猎，并依据αo、βo和δo灰狼个体的位置更新ωo灰狼个体的位置。
[0217]
3.2.6 进行迭代，按照3.2.3-3.2.5进行迭代处理，在满足迭代条件时，输出迭代结果αo灰狼，得到优化目标最优解。其中αo灰狼的位置坐标表示着cfrp铣削表面粗糙度最好时的工艺参数主轴转速nz、进给速度vj、轴向切深a
p
、径向切深ae。
[0218]
实施例1：
[0219]
本实施例是在vmc850e型立式加工中心上对长度200mm，宽度为100mm，厚度为5mm的t300复合层cfrp板，对cfrp进行侧铣。考虑到cfrp铣削加工为精加工，轴向切深和径向切深在工序图上设定好的，所以在本次实验中将轴向切深和径向切深设为定值，取轴向切深为固定值5mm，径向切深为固定值1.5mm，将主轴转速和进给速度取4个水平，设计全因子实验表，全因子如实验表如表1所示。
[0220]
表1 全因子实验设计因素及水平表
[0221]
[0222]
表2 全因子实验设计表
[0223][0224][0225]
根据表2实验安排，对长度200mm，宽度为100mm，厚度为5mm的t300复合层cfrp板进行轮廓铣削。本次实验需要利用三维形貌仪器对铣削后的表面进行测量，其评价指标为三维表面粗糙度，表3给出铣削长度200mm，宽度为100mm，厚度为5mm的t300复合层cfrp板后所测量的三维表面粗糙度，测量方法为：基于一组工艺参数铣削后，测量三次表面粗糙度并取其平均值。表3给出cfrp轮廓铣削实验后的数据。
[0226]
表3 cfrp轮廓铣削实验数据
[0227][0228][0229]
1 gsvr模型训练
[0230]
在gsvr模型训练之前，需要指定算法的参数。首先对svr算法参数进行设置：选择svr核函数类型，由于铣削工艺参数与三维表面粗糙度是线性不可分的，所以svr的核函数选择高斯核函数(radialbasisfunction，rbf)；设置svr惩罚系数c、核函数系数γ、损失界限ε的大小，常规方法是更具手动调参，寻找最合适的参数，但手动调参时间长，需要多次尝试，并且调试效果不佳。针对svr惩罚系数c、核函数系数γ、损失界限ε大小的选取，本方法采用gwo算法智能寻优，设置gwo的种群规模n＝50，并用佳点集法初始化种群；迭代次数t＝100；并且设置灰狼维度为3，灰狼位置坐标为x＝{x1,x2,x3}，第一个维度为惩罚系数c的大小，第二个维度为核函数系数γ大小，第三个维度为损失界限ε大小；根据经验和参数推荐之设置三个参数的上下限为：0＜c＜100,0＜γ＜10,0＜ε＜1；此gwo算法的适应度函数为mse的大小。
[0231]
整理表3数据集，将数据进行标准化，代入svr模型，利用gwo智能优化算法得到svr最佳参数，获得最佳拟合效果。表4是svr的测试集准确率和用gwo优化后的gsvr测试集准确率。
[0232]
表4 模型验证机准确率表
[0233][0234]
由表4可知，gsvr的准确率更高，因此选取gsvr模型作为gwo算法的适应度函数。
[0235]
2 gwo算法优化结果
[0236]
在gwo算法优化前，设置gwo算法参数。设置gwo的种群规模n0＝50，并用佳点集法初始化种群；迭代次数t0＝100；并且设置灰狼维度为4，灰狼位置坐标为第一个维度为主轴转速nz的大小，第二个维度为进给速度vj的大小，第三个维度为轴向切深a
p
的大小、第四个维度为径向切深ae的大小；根据立式加工中心和cfrp铣削经验设置工艺参数的数值范围：5000(r/min)＜nz＜8000(r/min)，100(mm/min)＜vj＜100(mm/min)，a
p
＝5mm，ae＝1.5mm。
[0237]
将gsvr训练模型作为gwo算法适应度函数，进行迭代输出。表5是gwo-gsvr算法优化得到的工艺参数值(主轴转速、进给速度、轴向切深、径向切深)，基于此参数得到的cfrp轮廓铣削的三维表面粗糙度如表5所示。
[0238]
表5 gwo-gsvr算法优化结果表
[0239][0240]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0241]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0242]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0243]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计
算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0244]
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。