车削切削力不确定性的模糊—灰色预测

HEATS

1 引言

# ^$ r" L! Q& A9 u) G

加工过程中出现的不确定性现象已成为一个突出的问题。切削力作为过程特征量之一，其不确定性来自许多过程因素(包括可量化、难量化因素以及可预测、难预测因素)。若忽略难量化和难预测因素的影响，可近似认为切削力的不确定性主要源于两个方面：实际的力测量系统的影响(包括标定方法、漂移补偿方法、软件中的平均值和圆整方法、测力仪温度等的影响)和过程参数变化的影响(包括切削用量三要素的影响)。 # Y$ r. s5 _5 J. z8 v

常见的对系统不确定性的估计方法有概率分布法、模糊估计法和最小二乘法等，但至今所见关于加工过程不确定性研究的文章或报道为数不多。为了评价切削力的不确定性，考虑测量系统的标定误差和切削参数误差，根据最小二乘回归法和误差补偿理论，有研究者提出了TLSRM模型，此模型产生的不确定性预测结果是一个确定值而不是一个范围；考虑磨削参数变化和忽略测量系统的标定误差影响，基于模糊理论和线性规划理论，有关文献提出了一个模糊预测模型GFPM，而此模型产生的磨削力不确定性预测结果是一个范围而非一个确定值。

模糊预测模型GFPM使用区间形式来预测和评价磨削力的不确定性，从工程实践的角度看，要比预测一个确定点更有意义。但是模型GFPM不确定性预测的相对误差仍然较大。为了进一步研究这个问题，笔者注意到灰色理论擅长于估计任意系统的不确定性，因而试图建立一种模糊—灰色预测方法，并将其应用于车削力不确定性预测。为了验证所提出的方法，作者还建立了一个车削力检测系统。 - n" j, k: _/ C% _% S5 e

2 车削力不确定性的模糊—灰色预测原理

车削中，主切削力的经验数学模型可表达为 " K& ?: U5 b. U$ i" u- ]9 y' X) O( G( O3 S2 k) V& P5 a6 s5 \4 ] P# m: _ c6 ^. u2 N

Fz=Cz(ap)x(f)y

(1)

式中，a_p为切削深度(mm)，f 为进给量(mm/r)，x和y为指数，C_z为系数。

对式(1)取对数，可得 9 V3 c: G! r8 ]* `5 P4 ]* H, ]% {' F# `. L2 G9 n0 C9 b! o+ K n; [

log(Fz)=log(Cz)+xlog(ap)+ylog(f)

(2)

令：Y=log(F_z)，b₀=log(C_z)，b₁=x，X₁=log(a_p)，b₂=y，X₂=log(f)，代入式(2)得线性模型分明形式为 5 l- c% e" U% Q' c$ p* r7 Z

Y=b0+b1X1+b2X2

(3)

0 a* m) _; I0 |* t9 R

定义式(3)对应的模糊形式为 5 d2 K W* {1 Y0 }+ u U3 G: r. q* M. e( U3 y

Y(Xp)=A0+A1X1p+A2X2p

(4)

这里选择对称三角形形式的模糊隶属度函数。设三角形底边的中心和宽度分别为C和W，则对称模糊数表达为Y(X_p)=[C(X_p)，W(X_p)]，A₀=[C₀，W₀]，A₁=[C₁，W₁]，A₂=[C₂，W₂]。分明变量为X_1p=log(a_p)_p，X_2p=log(f)_p(p表示试验代码，且p=1，2，3，…，n)。

为了使力的测量值落在式(4)所得预测区间内，并且使对称模糊数的宽度W_i之和最小，应当满足以下线性规划关系：

规划目标为 / H( x2 R2 w0 g& F' [8 c! I7 X0 r* K \3 I) B, m' P# C1 s9 k1 Y% R6 P, k

∑W(Xp)→min.

(5)

$ s3 l# u+ k7 @& [) ~. y; j

约束条件为 5 A% T3 ~; U% f6 m: A( B" y0 ~7 s6 `7 e3 x7 O5 s- k+ x+ T+ u" h: Q% ?! C4 ^' [6 o: C

log(Fz)p≤C(Xp)+(1-h)W(Xp)	(6)
log(Fz)p≥C(Xp)-(1-h)W(Xp)	(7)

式中，W(X_p)=W₀+W₁X_1p+W₂X_2p，C(X_p)=C₀+C₁X_1p+C₂X_2p，宽度W_i≥0，且中心Ci≥0(下标i=0，1，2)。模糊隶属度h 满足0z)_p为实测力的对数值。考虑到进给量01p=log(10^mf)_p且取m≥2。

此外，式(4)～(7)可以逻辑地推广到其它切削特征量(如切削温度等)的建模。

在求出式(4)中的模糊数Y(X_p)=[C(X_p)，W(X_p)]后，再取反对数，即可得到车削力不确定性的预测区间。依据给定的切削参数范围，安排一组测量试验，试验编号为p=1，2，…，n。对于第p 个测量试验，测量实际的切削力值；同时将相应的切削参数代入式(4)～(7)，以便得到切削力不确定性的模糊预测区间；再考察切削力的实测值是否落在预测范围内。 " w0 y' g; a( i' @0 U# f+ u+ o

已有的研究结果表明，以上描述的模糊预测模型具有如下特征：要得到较好的预测结果，建模的数据样本越多，则在建模数据范围以内模型的输出越满意，但在建模数据范围以外模型输出越不满意。而灰色集合理论擅长于预测不确定性，且能解决模糊预测模型存在的问题。

基于式(4)，假设模糊模型的输出序列为：EY={Y(X_p)}={[C_p，W_p]}={[Y_u(p)，Y_s(p)]}其中Y_u(p)=(C_p-W_p)，Y_s(p)=(C_p+W_p)，且p=1，2，…，n。选出q(q

对于序列Y_u0和Y_s0，它们的一次累加生成序列分别定义为 Y_u1={Y_u1(1)，…，Y_u1(j)，…，Y_u1(q)}
Y_s1={Y_s1(1)，…，Y_s1(j)，…，Y_s1(q)}其中j=1，2，…，q1 x2 q* v2 N7 C% v) k- y

元素Y_u1(j)和Y_s1(j)表达为 2 Y1 g8 d$ s% z+ x$ _8 Z/ i& d4 T( R r i$ T& u5 m% r! [' V, Q5 B' ]1 I+ ?. L. @6 {2 v1 q0 N8 Z4 ^8 A! ]/ e- Y/ i! H

Yk1(j)=∑Yk0(i)

(8)

其中j=1，2，…，q2 M5 R& S! | x u/ G* f

对于序列Y_u1和Y_s1，它们的相邻平均生成分别定义为 " V! {/ V* `" b5 l3 I' q( L! u) u+ p, P5 p4 Q+ }( p& O" U4 G8 O n3 Y

Zu1={Zu1(2)，…，Zu1(i)，…，Zu1(q)}	(i=2，…，q
Zs1={Zs1(2)，…，Zs1(i)，…，Zs1(q)}

其中，元素Z_u1(i)和Z_s1(i)可表达为 . G% t% V: o$ SZ_k1(i)=0.5Y_k1(i)+0.5Y_k1(i-1)(k=u 或s) 4 m. ^: v/ p! |/ ]4 D/ d3 @

假设列向量Y=[Y_k0(2)，Y_k0(3)，…，Y_k0(q)]T，且矩阵B定义为 7 d$ V, v+ {. E) c# m: w, _ F* MB=[-Z_k1(2)，1；-Z_k1(3)，1；…；-Z_k1(q)，1](k=u 或s) ! P+ S7 K; I5 F. O

灰色微分方程dY_k1/dt+aY_k1(t)=b 中的参数a和b的估计值确定为： . `& b$ ~* D+ C' K9 @ { Y- h/ X. ~8 Y

[a，b]T=(BTB)-1 BTY

(9)

) j9 t2 }. e0 @- a$ y- n) @) r

根据式(8)，有Y_k1(0)=Y_k0(1)。

灰色微分方程的解为 9 F6 |) t6 H1 R, o6 a9 ]5 d; D; m7 T! R, [" T9 }( Q9 t% U: |. B' e( r" y# [1 F5 q' F% Y8 @4 Y w- o. O; d( y9 S

Yk1(1)=Yk0(1)	(10a)
Yk1(i)=[Yk0(1)-b/a]exp[-a(i-1)]+b/a	(10b)

式中，i=2，3，…，q。

式(10)可用于求Y_k1的模拟序列Y_k1。因此，Y_k0的模拟序列Y_k0可确定为 : X9 M( ^) h* {- W9 O! }/ y1 Q A' y6 o1 k6 _5 @. W8 {0 K3 \* ?3 ~2 M- h: I/ U0 m3 R5 q1 @0 u- Y- [

Yk0(1)=Yk0(1)(11a)
Yk0(i)=Yk1(i)-Yk1(i-1)	(11b)

其中i=2，3，…，q；k=u 或s。i≤q 用于原始序列Y_u0和Y_s0的模拟；i>q 用于切削力不确定性预测。

3 切削力检测试验

. W+ Y d7 I2 r2 U2 a

6 I* w' c7 K- m& E* i3 `& W3 F: C. G- Q0 {& R' W- c* \) P, W G8 s3 h# U5 }2 k3 j, @7 ^; E0 q: g' T6 E* g m5 m: P# H

9 L6 M3 e" H1 \% I& E, M! B9 w, N, H' q: }( W: o+ w" k9 Y# \( i& N ^' m, k M( M2 ^4 F2 Q% `$ s3 ]9 c* Z9 b0 g, U1 I" V% o- O" }* D# C$ o2 b4 V5 c# z2 s* R) _* w z) o0 r2 I+ H! `3 n! s3 p& @& n: ~1 {" W3 n; B5 s) o' l$ q3 E4 b* S, A$ ] H" }9 ?0 C2 C5 s; J9 F5 Y. i4 x" P$ `* l0 c& P1 F* `+ ?) h3 @1 S9 a. H" g5 S5 c. Q2 c* s* N6 h2 z" N' p 表1 检测系统主要组成部分   组成部分 说明 1 测量传感器 电阻应变式 2 信号放大器 YD-4A电阻应变式电压放大 3 数字测量仪器 12位A/D 转换卡 4 检测软件 用C++和MASM 语言开发 5 采集参数 采样频率50Hz；样本500～1000 6 计算机 586微机

+ F8 S, U2 V) C8 t8 F1 s$ ~& Z8 D/ r! B3 ^3 Z/ A% f3 r# u2 v6 D0 ]! i( l! r0 {& q& L7 s. V# c% T7 a, K6 k% |5 J$ n( j! X& ?( M- z+ W$ y, P1 v; h" {, c0 I$ Q/ v9 l9 Z6 {8 s7 S( `) d7 [; ~; s6 j9 F4 X9 L7 b! z0 i9 O9 a0 f: t; _9 H& _& n8 ~7 g8 U/ `4 j" P- N: G/ l# F F2 F9 ?) C2 i9 o' O5 Z/ ^/ E% a( c! k$ y5 x- c; w3 `2 r# F" R, T 表2 车削力试验记录 No. 切削深度ap (mm) 进给量f (mm/r) 主切削力Fz (N) 1 2 0.1 439 2 2 0.2 878 3 2 0.3 1129 4 2 0.4 1443 5 2 0.5 1756

1. 切削力检测试验条件 4 Z2 M6 \* C# l5 c* p# }5 j+ G/ x
   为验证模糊—灰色预测方法的有效性，建立了切削力检测系统，以获取建模和评价数据。检测系统的主要组成部分见表1。其它试验条件包括：工件材料45钢(正火处理，硬度HB187)，工件直径81mm；CA6140车床，主轴转速380r/min；切削速度96m/min；YT15硬质合金刀具。
   1 ~& E: P. I- c0 }, l' K6 a1 \
2. 切削力试验与检测
   基于上述切削力检测试验条件和表2所列切削参数进行切削力试验，并将切削力检测结果记录在表2中。

4 预测和试验评价

5 w) Y% k; z" ]1 T3 _

1. 切削力不确定性的模糊回归模型
   根据式(4)～(7)，用表2中全部数据进行切削力不确定性模糊建模(考虑到在建模数据范围以外模糊估计较差)。为便于编程计算，令m=2且X_1p=log(10^mf)_p。表3 列出了切削力不确定性的模糊估计结果，对应的模糊回归模型为 , M* i6 `$ K1 J& D4 [& [Y(X<sub>p</sub>)=[logF<sub>z</sub>]<sub>p</sub>=A<sub>0</sub>+A<sub>1</sub> [log(a<sub>p</sub>)]<sub>p</sub>+A<sub>2</sub> [log(100f)]<sub>p</sub>式中，A<sub>0</sub>=[1.272603，0.1092808]，A<sub>1</sub>=[0.8845206，0]，A<sub>2</sub>=[0.856115，0]。模糊隶属度h=0.5。相对误差∆<sub>u</sub>(或∆<sub>s</sub>)定义为 7 ]$ D5 S: X2 ]. i p& M$ {∆=(估计值-测量值)/测量值 0 m4 z( g( v0 I9 @+ }- a4 n( g8 m0 @7 `7 }; ^2 a9 R8 |, ?( B2 v% y# G" _) J( c- s& v6 Y! B% g: j/ z7 w7 [7 m5 X8 D& K! A* q5 c+ M- i, L/ |, d( n2 S9 Z+ B% K$ L+ k0 ]6 a' N

   " ?0 V) R0 p, k& I: h: }6 [& W$ i! r( [5 s3 y% B5 X5 R+ d& O8 \! P8 ]' I' R1 S9 _8 N% L5 M# X A) I( l4 b$ u* H) x. i( B; z: r8 {$ X$ O) w& D, M: P: V1 P4 X% s% a8 s# R) d, M) `5 R% Q( J9 e+ `$ ?+ D m7 n6 B: S, l7 C0 ]% j2 ^! S/ K9 I* t2 l1 p, t$ Z' S, w8 h% g& l$ ]* i) P+ z' h# B9 |8 M6 y$ Z j3 \7 f) u ]9 B0 l) u3 F+ t h8 b0 D: f2 s: Q5 A- a/ t4 q0 `! o( [( `$ O: ]* X+ I( O& W5 ?$ o& {- G( s1 B; B" m% r' q! Q- K# I' N: p' S. K3 ^% W. G6 o 表3 车削力不确定性模糊估计结果 No. 实测力 模糊估计 下限∆u(%) 上限∆s(%) 1 439 [415,516] -5.5 +17.5 2 878 [749,932] -14.7 +6.20 3 1129 [1063,1324] -5.8 +17.3 4 1443 [1363,1696] -5.5 +17.5 5 1756 [1646,2047] -6.35 +16.6

   1 \0 v. R( Y/ W) K% r$ _2 ^9 e. P+ N) ^8 ?& g& L- x4 L( [5 y/ u1 r7 C5 {4 J: }& V) m8 z% F- ~# Z o4 I7 F! n' O: e2 h1 ^1 n/ q# ?0 q; n6 o$ V: @4 q$ C' b9 s& I' S1 _2 g6 F8 @$ @( u5 }; D1 c: h7 m' Z% }' O% P/ g3 _ S( Z4 B/ | b* p/ |- u$ u* F6 o. m# w) E+ k9 V! J3 K- `3 {% R0 S( s' B- Z4 x. `; e5 v, ], W0 O( f+ g# [# X4 ~1 N' g4 k' ?% o! {% P6 u6 }. Z; c, _* v7 s/ e" P1 O. ?0 A+ a6 [# R' m" u" M: N$ H- b 表4 车削力不确定性模糊—灰色预测 (含3个建模数据的部分信息模型) No. 实测力 模糊—灰色预测 下限∆u(%) 上限∆s(%) 1 439 [415.00,516.00] -5.5 +17.5 2 878 [740.10,920.80] -15.7 +4.9 3 1129 [1046.6,1303.4] -7.3 +15.4 4 1443 [1480.1,1845.1] +2.6 +27.9 5 1756 [2093.2,2611.8] +19.2 +48.7

   " _# `' _# k# K5 w" X# `: m# i( {/ k @4 v' O1 c+ W: a6 S R) x# A6 t4 T) u0 ^6 ^( g/ x z( U) J9 Y9 \1 ]6 {& y) p r5 y/ K* ]* V: R; n$ M& f6 m* S8 y1 \+ I C3 i# z' e" ^! {+ A8 g1 ]' B0 o/ k# G2 s, g1 t' V: @, @0 G9 b: c9 a/ r, v5 G% C/ ^' ?2 j: u# M, y& E0 Y/ H: b+ g# s( j- F: k d$ g' b- s: z) x. \9 {( R9 B6 g) E( C9 k6 ?: n, w4 j+ ]0 E* U6 v6 q- F I6 Z; {2 e) D0 C) E, x1 d3 F0 j# K* ?& N4 i1 y! z+ f: B/ D/ i4 w" n; R9 n+ A) ~( ]; f+ ^ R+ u1 \8 L7 v0 H 表5 车削力不确定性模糊—灰色预测 (含4个建模数据的新信息模型) No. 实测力 模糊—灰色预测 下限∆u(%) 上限∆s(%) 1 439 [415.00,516.00] -5.5 +17.5 2 878 [765.60,953.30] -12.8 +8.5 3 1129 [1020.6,1270.5] -9.6 +12.5 4 1443 [1360.4,1693.2] -5.7 +17.3 5 1756 [1813.5,2256.6] +3.3 +28.5

   2 p1 K1 L1 v3 X- B, V b2 D! i/ ~7 f& T# t/ B1 F1 y3 L$ ^* ?, ^% ^6 L5 |" [- Q1 [5 U3 i" p( ]1 e) s. I( v/ J& _5 \) I2 p0 {6 ~% ]# n/ F, N0 y; o2 r# U" V! R) D- ?1 F# K9 [9 h% d7 m. [5 D% _4 p- N2 j/ F8 D: ]2 U# n2 m p5 m+ E9 g' P: i2 V7 ?1 }% M3 C2 g9 T: Y3 p1 O+ d' e% E2 H0 X2 M7 M! l, E9 o E3 d% w* ]) n8 s" F) W) S8 @' y" W: o+ R$ s: S6 w K 表6 车削力不确定性模糊—灰色预测 (含3个建模数据的新陈代谢模型) No. 实测力 模糊—灰色预测 下限∆u(%) 上限∆s(%) 2 878 [749.00,932.00] -14.7 +6.2 3 1129 [1056.9,1316.4] -6.4 +16.6 4 1443 [1353.4,1684.2] -6.2 +16.7 5 1756 [1733.2,2154.6] -1.3 +22.7

2. 车削力不确定性的模糊—灰色预测 " s$ m( ^& i- I8 N& d# A
   1. 选择表3中No.1～3的数据用于模糊—灰色预测建模，则车削力不确定性的模糊—灰色预测结果如表4 所示，其中数据No.1～3对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.1～3)，而数据No.4～5表示模糊—灰色建模数据范围之外的预测值。此外，灰色微分方程中的参数确定为 [a_s,b_s]=[-0.3475,590.7376]
      [a_u,b_u]=[-0.3466,475.3764]
      相对误差定义为 , ^* ^, s) i/ D: V∆=(预测值-测量值)/测量值
   2. 选择表3中的数据No.1～4用于模糊—灰色预测建模，将车削力不确定性模糊—灰色预测结果列于表5，其中数据No.1～4对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.1～4)，而数据No.5 表示模糊—灰色建模数据范围之外的预测。灰色微分方程中的参数确定为 [a_s,b_s]=[-0.2872,674.6824]
      [a_u,b_u]=[-0.2874,541.5835] ; x- v" D, E& p M
   3. 选择表3中的数据No.2～4用于模糊—灰色预测建模，则车削力不确定性模糊—灰色预测结果列于表6，其中数据No.2～4对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.2～4)，而数据No.5 表示模糊—灰色建模数据范围之外的预测值。灰色微分方程中的参数确定为 1 L% {6 ?5 l% n( h- H[a_s,b_s]=[-0.2464,931.306]
      [a_u,b_u]=[-0.2473,746.3059]

! l$ \+ X" o# S3 V- p2 N

5 预测结果分析

t% @* [* I1 `! f# e+ @! B

由表4数据可知：数据No.1～3是原始数据的较好模拟(原始数据来自表3中No.1～3)，但预测区间稍微向左漂移；虽然数据No.4在建模数据范围之外，但它的预测较好且含较小的相对误差(+2.6%)，即实测力落在预测区间下限的左侧邻近；数据No.5的预测较差，即实测力落在预测区间下限的左侧且含较大的相对误差(+19.2%)。 4 c. Y; q$ r& {7 q

由表5数据可知：数据No.1～4是原始数据的较好模拟(原始数据来自表3中数据No.1～4)；虽然数据No.5在建模数据范围之外，但它的预测较好且含较小的相对误差(+3.3%)，即实测力落在预测区间下限的左侧邻近。

由表6数据可知：数据No.2～4是原始数据的较好模拟(原始数据来自表3中数据No.2～4)；虽然数据No.5是在建模数据范围之外，但它的预测较好且含较小的相对误差(-1.3%)；实测力落在预测区间下限的右侧邻近。

基于上述分析可知：模糊—灰色预测方法可以较好地模拟建模数据，同时可以在建模数据样本范围之外有效地预测切削力不确定性区间，不过真正的有效预测点仅局限于最后那个建模数据。此外还可发现：表6中的新陈代谢模型的预测结果优于表4中的部分信息模型，而表4的部分信息模型又优于表5中的新信息模型。与常规的最小二乘回归法、神经网络法和模糊回归法明显不同，模糊—灰色预测法仅需要少数试验数据样本来建模，并能有效地预测将来的不确定性区间。 K o S+ `# _8 `; S

6 结论

在建立的车削力检测系统中验证了所提出的切削力不确定性模糊—灰色预测方法的有效性。模糊—灰色预测结果与实测切削力数据的比较表明：所提出的预测新方法不仅能较好地模拟建模数据样本，而且还能以较小的相对误差预测建模数据范围以外的其它数据点(即预测未来数据点)。此外，与传统建模方法明显不同的是：模糊—灰色预测方法在建模中仅需要少量试验数据样本，可大大减少切削试验的成本和工作量。