不同切削力预测建模方法的比较研究

104377257

1 引言
# V" M8 p( r; E* ~在切削技术研究及实际切削加工中，有关切削力的数据是计算切削功率、设计和使用机床、刀具和夹具、开发切削数据库、实现加工中切削力控制等的重要依据。在实际生产中，为了在粗加工时充分利用机床功率，在精加工时有效保证工件质量，均需合理选择切削条件，并对选定切削条件下的切削力进行预测。
6 y+ g* @5 x3 @7 {4 y预测切削力的经验模型主要建立在最小二乘回归法的基础上，近年来，人工神经网络法和灰色理论建模法的应用也越来越多。这些建模方法具有不同的特点及使用条件，并各有利弊。本文结合实例对人工神经网络法和灰色理论建模法的建模特点及其优劣进行了较深入的分析，并与常用的最小二乘回归法进行比较，旨在为合理选择建模方法提供参考依据。
# a% c8 O% w* B+ y% u4 Q* c* O2 基于径向基神经网络的切削力预测建模
基于Kolmogorov定理的三层BP神经网络可较精确地拟合任意连续函数，当输入节点数为n时，隐层节点数为(2n+1)且常选择Sigmoid型传递函数。在实际应用中，往往需要大量的BP隐层节点，通过增加隐层数可减少各隐层上的节点数，但迄今尚无选取BP网络隐层数及其节点数的统一方法。此外，标准BP以及各种改进型BP算法均存在局部极小和收敛速度的问题。
径向基神经网络(RBF)精确拟合任意连续(或不连续)目标函数的能力及学习速度均优于BP网络。RBF的隐层节点采用径向基传递函数，其节点数不像BP网络那样需预先设定，而是在学习过程中不断增加直到满足误差指标为止。
/ {) q6 U3 l6 X+ W$ s. k根据切削力及其影响因素的特点，设计如下图所示的RBF网络。由图可知，RBF网络包括输入层、一个RBF隐层和输出层。输出层包含一个用于输出预测切削力的线性节点。隐层包含S1个RBF节点且S1值在学习过程中动态增加。输入层的R×Q阶输入矢量阵P表示有R个输入节点，在每个节点处输入Q个样本(Q等于试验组数m)。每个输入节点代表切削力的一个影响因子，且切削力的所有可量化影响因子均可抽象为一个输入节点。若考虑切削深度、进给量、切削速度、工件材料的剪切屈服应力、刀具材料、刀具的负倒棱宽度、主偏角、刃倾角、刀尖圆弧半径、刀具磨损、切削液等各种影响因素，则可有多个输入节点。根据实际建模经验，可主要考虑切削深度和进给量的影响，此时输入节点数R=2。

切削力预测的径向基神经网络结构图
  n  e$ t$ l1 ^(S2=1，S1动态确定)
5 X. i) I; I" }" o. N2 h% h设在m组切削条件下测得的试验数据为[P，T]，目标输出为T=[Fz1，…，Fzm]，R×Q阶输入矢量阵P表示在Q=m组试验中每组考虑R个切削力影响因子。在选定有关设计控制参数(如期望的网络输出误差平方和指标等)后，根据试验数据[P，T]，采用RBF设计算法可在较短时间内确定RBF网络隐层及输出层上的权矩阵W、偏差矩阵b及隐层上的节点数S1，完成切削力预测的神经建模。
8 h! `+ @1 p8 e3 切削力预测的灰色建模
4 B' \1 y0 a0 }4 h5 v* k灰色集合理论擅长于处理具有“部分信息已知、部分信息未知的小样本、贫信息”特点的不确定性对象，它通过对“部分已知信息”的生成与开发，从中提取有用信息，最终实现对研究对象内在规律的有效描述。切削加工实践表明，由于各种因素的影响，切削力通常表现出不确定性特征。GM(1，1)灰色建模原理简介如下：
设实测的原始数据序列为
) }& ]2 T8 f0 d$ p, U/ jY0=｛Y0(1)，Y0(2)，…，Y0(j)，…，Y0(n)｝
. ^! V+ y, D6 {& A: _$ [它的一次累加生成序列定义为
! b, z/ k4 j7 I1 f8 [- N8 LY1=｛Y1(1)，Y1(2)，…，Y1(j)，…，Y1(n)｝
其中，Y1(j)为
Y1(j)=
9 z6 e2 S9 ~; [! j1 t( pj
Y0(i)
Σ
; j5 k9 ^/ W! ]9 C& n8 g' y% Xi=1
对于序列Y1，其相邻平均生成定义为
Z1={Z1(2)，…，Z1(i)，…，Z1(n)｝
# L4 M$ `; ]# z$ v3 v; w! A% z! c0 \) E其中，Z1(i)可表达为
Z1(i)=0.5Y1(i)+0.5Y1(i-1)
假设列向量Y=Y0(2)，Y0(3)，…，Y0(n)]T，且矩阵B定义为
B=[-Z1(2)，1；-Z1(3)，1；…；-Z1(n)，1]
在灰色微分方程dY1/dt+aY1(t)=b中，参数a和b的估计值确定为
' X- |, e4 O( X[a，b]T=(BTB)Z-1BTY根据式(1)，有Y1(0)=Y0(1)。灰色微分方程的解为
) c3 y+ v& {6 gŶ1(1)=Y0(1)
- c. }  w0 b, R4 {! h4 {( I* o  L# V(2)
" B4 X6 W1 B8 O2 f) W* S9 R: TŶ1(i)=[Y0(1)-(b/a)]exp[-a(i-1)]+(b/a)
(3)
其中：i=2，3，…，n。式(2)、(3)可用于求Y1的模拟序列Y1。因此，Y0的模拟序列Y0可确定为
1 c0 x2 q' ]  h2 K) s- mŶ0(1)=Y0(1)
(4)
& D6 H1 f; Z  P& a/ CŶ0(i)=Ŷ1(i)-Ŷ1(i-1)
8 [6 E! `3 e9 G5 m& U! T3 B; s( I(5)
7 C1 q7 e0 _' ]其中：i=2，3，…，n。利用式(4)、(5)，i≤n用于原始序列Y0的模拟；i>n用于切削力不确定性预测。
; S7 V/ P. ~9 {! M$ Q* u4 模型的验证与分析
% W& T# Z6 S6 k% q/ O为了验证建模方法的有效性和准确性，需要获取建模数据和评价数据。表1为车削试验得出的切削力数据。试验条件：工件材料45钢(正火，HB=187)，工件直径81mm；YT15外圆车刀(416A)，前角15°，后角68°，副后角4°～6°，主偏角75°，副偏角10°～12°，刃倾角0°，刀尖圆弧半径R0.2mm，负倒棱宽度为0；主轴转速n=380r/min，切削速度v=96m/min。
) @  [8 [$ o# Y3 _2 b9 J; N9 y# b表1 切削力测量数据
& r# T1 O0 E! YNo.
1
2
3
4
5
* P9 ]8 {" c' C$ _' Z  c2 V! U6
7
# C1 W0 m, n; r) V8
* i% }4 o( ]1 p" r* d- s( Q9
ap(mm)
2
  |9 M1 J7 `, q' z% ^% [% v2
  z. M2 |/ y, E2 {9 ?( N+ O# ^% z3 x2
6 }6 g6 c8 K" l9 Y/ }* O2
2
1 Q. t) [! i" Y4 _( ?3
3 p  `7 y3 \9 P4 l! C  Z% T0 c3
( e* L3 ]3 F' I/ N) Y4 W0 T3 |* c3
& {! {; E/ }, J3
f(mm/r)
0 ]2 W  B0 ^! T" |+ H0.2
" A9 W& |1 Z, f  O8 ^* H0.3
9 E* S* W9 H, h) V4 c0.4
. h+ m& M: g  C, T$ N* g3 }- x9 k  O8 m" `0.5
0.1
# `% h& H2 I- N+ B% T6 ?: k- Q% S0.2
- v5 V* J' a2 g  |0.3
6 f# {- \& a9 B2 H4 K9 e, X) G7 l0.4
0.5
( @! t3 a2 \# s) |( ~主切削力Fz(N)
878
1129
1443
1 m4 k$ @) G# t1 ]# r) G% G+ g  X1756
& W' f" E; g4 D0 a) P- ~( J627
1255
1756
+ F5 Y4 Q1 e* n2 N2195
2760
神经网络模型的验证与分析
基于图1所示模型和表1数据，选择ap和f作为输入层节点，实测切削力Fz作为目标输出。选择表1中的数据样本No.1～8用于建模，数据样本No.9用于模型评价。采用径向基学习算法设计图1所示具体模型时，学习控制参数如下：网络输出误差平方和期望值e=0.01，径向基散布值sp=1.0，隐层最大节点数nr=1000，显示频率df=25。通过编程计算，得到表2所示切削力神经网络具体模型参数。基于表2模型和与表1切削条件对应的切削力神经网络计算结果列于表3。为便于比较，表3还列出了切削力实测值以及采用普通最小二乘多元线性回归模型的计算结果。最小二乘回归建模是基于表1中的数据样本No.1～8，数据样本No.9用于模型评价。回归模型的线性形式为
Yp= 5.135186 + 0.9719143[ln(ap)]p+ 0.862146[ln(f)]p
4 `9 Y- {- V+ _( K- K5 @(6a)
0 V7 ^! {- Q: P& f' n5 F+ zFz= exp(Yp)
(6b)
7 O0 Q7 ^7 O3 [) g! o6 |  ^: }7 q1 e表2 切削力神经网络具体模型参数
/ t1 s2 }1 m8 L6 {+ C隐层上权Wh
隐层偏差bh
输出层Wo(×106)
W11=0.4，W12=3.0
1 y% x+ D$ J( ?& U8 V" A! i; B! Cb1=0.8326
W11=0.4487
W21=0.3，W22=3.0
' V3 K* O' k' `b2=0.8326
- ~' R3 U$ f) m  U' QW12=-0.3428
' d" i4 l2 o3 M" P& Q' AW31=0.5，W32=2.0
b3=0.8326
W13=0.6923
: h- }. @6 J7 o* f1 R- SW41=0.4，W42=2.0
b4=0.8326
0 ?3 a% V% J) u6 J& C2 \/ ^& c; FW14=-1.0032
2 t8 {2 j( |5 f) n, ~% i# p  D$ HW51=0.1，W52=3.0
b5=0.8326
% h+ x' m2 }. u9 x# uW15=0.5320
W61=0.2，W62=2.0
b6=0.8326
+ v0 @' q: H( I! MW16=0.3526
$ z: k5 ^+ k/ m2 _W71=0.2，W72=3.0
b7=0.8326
. r! [) n) O( ^6 s( ~' X3 X7 k- t; lW17=-0.5968
- a! `% I! \1 |! [% l6 J; s其它结构参数：两个输入节点为ap，f；隐层节点数S1=7(通过学习确定)；隐层节点传递函数为radbas；输出层节点数S1=1(预选)；输出节点传递函数为线性函数；输出层偏差b0=-4.0992×104
% w( ]/ |, F- y表3 切削力Fz的神经网络预测值、最小二乘估计值与实测值的比较
+ W& J/ h/ p8 v  r) C) YNo.
% v& W7 c" s2 ^$ T切削力Fz(N)
$ z, e$ d$ ]/ t3 f# R2 Q6 L  o实测值
神经网络预测值
/ F" L$ k  X" j2 |% c相对误差B%
最小二乘估计值
0 a/ H1 E! J% d, I相对误差B%
1
" ~+ O3 o: C: h' u878
! D$ S4 ~# q4 q- _& t3 Q4 R, @878
0
815
& C4 [6 k- c* b3 l-7.2
$ Q7 G- P3 o6 H" p/ r  n2
3 S; C3 Z* @5 P: ]! y3 t1129
4 b4 n! p* w7 s$ k- U1129
3 ~9 ]% @' V2 Q6 g# `/ i0
1157
1 C0 m" {- W$ l( ]3 ]& \+2.5
; v6 K2 l' o* j4 k  J$ s; e, h3
1443
% G! |/ C- h2 p, X+ u2 J1443
0
) a. u- H2 d9 s% j4 r1482
3 Y& J  O7 N" A0 o  s+2.7
4
6 m/ H; |2 K" W8 J2 H1756
7 _4 c: r4 ]) B$ D. V1756
. G# f5 a; l: v$ L; i( s0
$ b& [4 j$ Z7 Y  m( D1797
; y2 P8 Y$ W6 n- c$ \8 H2 ~+2.3
& z; |$ z0 r  g. q9 I7 @5
* }4 U' ]3 q4 d( ]  K; a  @; \0 k627
% ~! ~# g6 S5 O  A' z# b627
, E" e$ F2 B$ a) [0
7 f6 |; |3 w0 F" H! E665
+6.1
6
1255
" Q$ i- Y8 l, J5 i7 m1255
6 m8 f- ?0 h0 c% I1 a6 n0
/ C  L$ A6 ]6 U( M  W" O$ h1209
-3.7
7
2 m6 e0 K! T! G' V3 ]9 P, q1756
) J+ ^/ }' V# i" u1756
9 w$ p; |6 b" W/ I1 F, H0
1715
-2.3
, k+ e0 S" L  i8
2195
2195
0
) E; K4 \% ~* b, O9 A2198
+0.14
9
9 Z8 I% U+ D6 c8 w1 h  t) N2760
2415
& E% P" M$ K3 l  }-12.5
3 l; u5 A7 U& ]! C/ [# V7 b' R2665
-3.4
表3中，相对误差定义为：B%=[(预测值-实测值)/实测值]×100%。
由表3可知，径向基神经网络(RBF)建模方法具有如下特点：①可精确拟合任意连续或非连续函数(如数据样本No.1～8)，其拟合精度高于常用的最小二乘回归法；由于RBF隐层上节点采用径向基传递函数，故节点数不需预先设定，而是在学习过程中不断增加直至满足误差指标为止。②当所选切削条件在建模试验样本的切削条件上限或下限之外时，径向基神经网络的预测效果较差，拟合精度低于最小二乘回归法(如数据样本No.9)。③为使神经网络法的预测范围较宽、预测结果较准确，选取建模用的各个试验样本之间切削条件的差异不应太大，且应采集尽可能多的试验样本，但此时RBF网络的隐层节点数将增多。
灰色预测模型的验证与分析
根据前述GM(1，1)灰色建模原理，选取表1中的数据样本No.1～8用于建模，数据样本No.9用于模型评价。可得到表4所示的切削力Fz的灰色模型预测值以及与最小二乘估计值和实测值的比较结果。由表4可知，灰色模型的拟合和预测精度低于常用的最小二乘回归法(即相对误差较大)，这是因为灰色模型的建模数据分布不能较好服从e指数规律分布所致。
& E4 B4 |0 ]( i8 e1 }& s从表1去掉数据样本No.4～6，将数据样本No.1～3和No.7～8用于建模，No.9仍用于模型评价，从而使建模数据更接近e指数规律分布，得到的灰色模型预测结果列于表5。显然，相对于表4而言，表5的灰色模型拟合和预测精度有了明显提高，且建模数据样本更少，这正是灰色建模方法的显著特点之一。
灰色模型预测法适合“少样本、贫信息”建模，且可获得较高的模型拟合和预测精度，但其先决条件是建模数据样本必须较好服从e指数分布规律。为进一步验证其建模特点，表6列出了另一计算实例，即采用灰色模型预测外圆磨削加工中的法向磨削力Fn。其中，数据样本No.1～5用于建模，数据样本No.6用于模型评价，法向磨削力Fn的实测值引自参考文献。从表6可知，只要建模数据样本较好服从e指数分布规律，则灰色模型预测法的拟合和预测精度优于最小二乘回归法。