基于多参量状态信息融合的刀具磨损状态智能识别

HEATS

1 概述

刀具磨损是金属切削过程中必然发生的现象，现场统计证明，刀具失效、破损是引起NC类自动化加工机床利用率下降、生产成本提高的首要因素，是机床设备和切削加工安全事故的重要起因，由其引起的停机占故障停机时间的1/5～1/3。因此开展刀具磨损的在线检测在自动化加工中具有极其重要的现实意义。 4 Y* N. f& E; b) E7 v

刀具状态的检测方法主要有直接法与间接法。直接法是直接检测刀具磨损面的数值大小；间接法是通过检测与刀具磨损状态相关的物理信号来估计刀具磨损状态。当前的各种监控技术主要是集中于间接检测。其中最为常见也较为有效的有切削力检测、声发射检测和电机电流/功率检测等。切削力一直是表征切削过程的最重要的特征，切削过程的变化都是和切削力密切相关的，所以用切削力监测刀具状态有着灵敏度高、抗干扰性好等较为突出的优点。刀具的磨损量不仅与加工工况相关，而且与刀具使用时间的长短有着密切的关系。 7 o/ g! i- Y2 z: C/ F- d' k

本文主要从实用的角度出发，建立了基于切削力和基于相对切削时间的两种磨损模型。切削力模型是利用回归算法和模糊分类技术建立的，该模型描述了切削力、切削参数与刀具磨损状态之间的关系，通过检测切削力信号可在线识别刀具磨损状态。基于相对切削时间模型利用回归技术直接建立刀具磨损量与切削参数及时间的关系，可在较大的切削条件变化范围内实现对刀具磨损的识别。

2 两种磨损识别方法理论模型的建立

当加工过程中出现故障时(如刀具磨损或破损)，将导致加工状态信号的变化。因此传统的监控系统都是从检测的状态信号的特征参量的变化中获取故障特征，进而实现对加工状态故障的识别与决策。但是，在加工过程中由于切削参数的变化也将引起加工状态信号的变化，这种变化往往比由加工故障导致的加工状态信号的变化要大得多。传统的监控系统由于未能很好考虑切削条件变化对加工状态信号的影响而只能监视单一工况下的加工过程，而无法用于先进制造中切削条件多变的加工过程。就此，作者一方面利用模糊回归网络技术，根据磨损量建立了基于切削力信号的数学模型，利用切削力信号间接判断刀具磨损状态；另一方面，利用回归技术，基于相对切削时间建立了刀具磨损量与切削参数及相对切削时间的数学模型。 * v8 D8 u4 v" t7 G

$ S' r* T2 t$ ~3 @- u! j) b- E8 T2 w# B6 z7 P$ b, w% w+ N$ m* W* R; F. \7 b

图1 刀具磨损状态隶属度函数

1. 磨损状态模糊分类
   刀具磨损是相对缓慢的过程，通常磨损过程分为3个阶段，即初期磨损阶段、正常磨损阶段和急剧磨损阶段。三个阶段属于不同的磨损状态，各状态间的界限有一定的模糊性和重叠性。对于自动化加工而言，检测刀具磨损量不一定要精确的量，只要知道其在一定的磨损范围内即可。针对铣削，磨损状态很复杂，根据铣削过程的要求及实验数据，把刀具磨损状态分为A、B、C、D、E五类，各类的平均磨损量分别为：0.1mm、0.15mm、0.2mm、0.25mm、0.3mm。根据磨损状态的模糊性，建立如图1所示的磨损状态隶属度函数。之所以选择B类为梯形是因为正常加工时，在此区间其模糊程度相对要求低些。
   0 [5 g2 N; | M! J$ h5 ?
2. 基于切削力信号磨损模型 * M! k9 {; ?3 P( _
   切削力是表征切削过程最重要特征的物理量。在切削过程中所产生的多数变化都和切削力密切相关。大量研究证明，随着刀具磨损的增加将导致切削力增大，刀具破损将引起切削力幅值的骤变。切削速度v、进给速率f、切削深度a_p的变化都会影响切削力的大小。同时，切削力的大小也随工件材料、刀具材料等具体切削环境的不同而有所不同。根据金属切削原理，对于新刃刀具切削力与切削用量的关系有如下公式： 8 [% S! B1 R! B6 R! C, e8 {& p) N4 d: h6 _9 o6 j9 {' A. V# Y- U

   F=Cvxfyapz

   (1)

   式中：C——决定于刀具几何尺寸和材料性质的系数； 3 d$ c8 c" H( Z
   x,y,z——切削用量的指数。
   式(1)表示了加工过程中切削力同各切削参数的关系，该切削力模型属于静态非线形模型，经线形化后，模型变为： * z* ~9 Q* y. c0 s% P# h4 _% Z- C+ U% I0 b2 J' C* {! O( r# J: O' k' C3 N; s: @% r

   InF=InC+xInv+yInf+xInap % F( J( ~: N G M& R1 I2 B* u

   (2)

   9 C4 T0 u5 O. Z6 ~* S
   针对不同磨损状态刀具，在A、B、C、D、E五类中，相应的有如下表示式： " O' D/ B! X& i7 N. N8 b4 F6 M3 S% r: S& c7 c' A" b

   A类：Sa=InFa=a11+a12Inv+a13Inf+a14Inap & a$ K2 D7 C7 { B类：Sb=InFb=a21+a22Inv+a23Inf+a24Inap 9 o) U( l7 S- {7 t' \8 N" [ C类：Sc=InFc=a31+a32Inv+a33Inf+a34Inap D类：Sd=InFd=a41+a42Inv+a43Inf+a44Inap E类：Se=InFe=a51+a52Inv+a53Inf+a54Inap ) |" M! x1 @# b8 q- j

   (3)

   0 L; ^5 r4 [+ Y5 c6 C
   利用矩阵形式表达为： 3 i% H6 T/ v# U0 J: L$ b1 a

   S=W*X

   (4)

   
   ) @2 n5 b% S: t) `# C! b% z( d: _

   图2 切削力模型网络示意图

   0 ]9 u2 L6 T9 R2 M9 f. ?
   式(4)可以用图2的神经网络形式来表达。对该网络模型不采用传统的BP算法来调节权重值，而是直接利用回归技术来调节权重值，这是由于该网络输入输出之间的关系已经是线性的，而且基于回归的神经网络的学习速度快，可以在线学习识别。
   切削力模型是以切削力幅值为模型的特征值，在已知v、f、a_p的情况下，利用式(3)推得F_a、F_b、F_c、F_d、F_e作为图1中的模糊聚类的中心，将检测到的切削力值与聚类中心比较，从而来确定当前刀具磨损隶属度而达到检测目的。
   $ M; x$ ^' f- K2 N
3. 基于相对切削时间的磨损模型 - }& ~, r, W8 N0 E/ N+ B
   根据加工过程中刀具磨损规律的先验知识可知，刀具磨损量随加工时间的发展而变化。在不同的切削条件下，增长的快慢不同，研究发现，在同一切削条件下，磨损量的变化率近乎是一定值，但基于实际工况考虑，将磨损量的变化与时间的变化的关系定位为指数关系。再根据刀具寿命公式建立式(5)所示的基于相对切削时间的磨损检测模型。模型的输入量为切削速度v,进给速度f,切削深度a_p,初始的刀具磨损量VB₀切削时间间隔Dt,模型输出为刀具磨损量估测值VB(t)。 * r6 U% l) H2 U0 r5 \1 f+ g+ \5 A1 G+ V$ i

   DVB(t)=VB(t)-VB0=KwvxfyapzDtm

   (5)

   $ i! }' i( q% y6 k1 X+ w
   线性化为： 1 r4 x4 s) T. {( O4 Q* b# a7 ~6 T0 M0 `! k

   InDVB(t)=InKw+xInv+yInf+zInap+mInDt 9 Z7 N+ Z( K, l" f* Q n% b

   (6)

   式中，Dt—同一切削条件下磨损检测的切削时间间隔； # K! v0 O0 {6 u/ C2 u' S* w
   DVB(t)—同一切削条件下一定时间间隔内——磨损量的改变值； . C7 w& N: i! r3 F7 ` u% Z
   K_w—与刀具几何尺寸和材料性质有关的系数。

; Y x3 ~4 ^/ I" r7 z( S% k. q( d; U) W2 Q3 p

图3 铣削刀具监控实验系统

) x, c; n0 \( N, ^/ D- o( l. b

3 磨损模型系数确立及实验验证

% `1 C" Y7 O Q" ]8 R

5 U( T( S" F! X% {4 l& O
1. 实验装置及实验条件
   实验是在HURCO-BMC-20立式加工中心上进行的。一台Kistler9257A型三向测力仪，三台YE5850型电荷放大器，放大倍数为四十倍的工具显微镜，一台PC586微机。实验原理示意图如图3所示。传感器的信号经放大和低通滤波后将由微机控制进行采样和处理，因铣削过程有刀齿切削和刀具回转两个周期，对信号实行与刀具回转周期同步采样的方法来获得每转平均铣削力的大小。以两齿螺旋刃键槽铣刀为研究对象，以不同切削参数进行测试实验，根据实测数据识别模型系数。实验条件如表1所示。 / J* z! Q) l1 ` G7 M* b2 `, c* [4 }; I4 ?# I1 V9 Z5 R, i. r2 P4 N$ s8 q u6 v, L

   表1 铣削实验条件

   刀具 5 f. g! f2 G4 |) w t	材料 6 H3 ]0 I: u: {% p( w" w, o	铣削方式	工件材料	切削速度 m/min	进给速度 mm/min # a7 D/ [/ e) J	切削深度 mm	切削液
   键槽铣刀f10	高速钢 * J0 c* ~ C% h& |	顺铣 . E; D G$ y4 q1 N. `	不锈钢	21.98 # ^/ s0 [5 [) O5 `	35 / D* h) N( U- J	3.5 + _. f# i2 w( s+ S; W0 P% p	无
   				18.84	30 % ^! d- t$ f8 f5 J; Y* Z# A! J, s	3 ! S! B2 s7 O* ], A7 T
   				15.70	25 / T/ N ~6 @/ l1 N" y/ n$ a; \	2.5

2. 模型系数获取及其校验 7 M! I9 O% H! {
   为了验证所创建模型的合理性和实用性，进行了大量的实验。图4是利用磨损量分别为0.15、0.2、0.25、0.3mm的刀具，在不同组合的切削条件下进行实验所得到的铣削力模型训练用样本，共45组。图中的每条曲线是在同一组合的切削条件下进行的。利用回归算法得网络权值矩阵为： 1 u0 R1 I- `: h# n , N( Q" c9 K8 W
   上述权重值每行的相关系数R<sup>2</sup>分别为0.9871,0.9765,0.9654,0.9453,0.9645。因此可以看出所创建的铣削力模型反映出力信号与切削参数的密切关系。
   2 b3 C; n' d8 C; S4 d' Q. R3 E1 `. F; Q9 M& B0 T M2 B# s/ H' z+ ~

   图4 切削力模型训练用样本

   图5 时间模型训练用样本

   图5是在不同的切削条件下、不同的时间间隔内对不同磨损初值的刀具进行实验所获得的36组训练样本。图中的每条曲线也是在同一组合的切削条件下进行的，直接利用回归算法训练得式，6，回归系数依次为-2.1348,0.8607,0.1111,0.6143,0.8333,相关系数R²=0.9857表明磨损变化量与切削参数及相对切削时间有较密切的关系。 ; J0 I- b/ r3 K4 g1 f( @
   为了对上述模型进行验证，表2列出了具体的校验用数据样本。校验结果分别如表4及图6、图7所示。表3和图6表明了铣削力模型对刀具磨损类别有较高的识别准确率。图7表明刀具处于正常磨损阶段时，模型有较高的识别率；在初期磨损和急剧磨损阶段，磨损识别效果差。但因为实际加工过程中，我们关心的磨损状态主要是在正常磨损阶段，所以该模型可在工程上应用。但由于实际加工的不稳定因素的出现(如硬质点等)，时间模型不能给予体现，所以实际应用时，应将切削力模型和时间模型结合使用。

, S. ], |8 E/ f# d( U* H; m" e. V2 F3 Q* R& h: \! S) e8 e2 f& r. Y3 M: X

% j1 q! m4 a+ a8 J4 D& B5 Q6 @- N' M! D0 i$ _6 H3 ?* N. }. c. S/ z9 M4 R- Z' {* x+ S* s" m5 \+ p9 I8 Y0 c4 c3 e& H) H O" P c) ~7 v* n( C/ m1 _; u) }. B# v0 O _+ o9 ~; r6 j6 a8 Q, \) O5 b }- V) ]. e& l0 ]0 F1 |1 s( K, ~! c( ?" J' j- D9 v3 R& c0 K* u' Q+ L0 b! N5 v4 Q# _- O8 \8 ~: X, `! I7 ?6 d% }6 G1 b% o! F5 K' i- a( t; o+ \! y' d9 k0 H8 F1 K$ z* @: {- S& ` o. T! q% B |, }- z5 d) h, G1 }! v2 X" a4 x# R% G9 Z1 C8 u, m$ s ]0 _0 P! [' j8 ]' {9 I( M' I+ i- `, r" c# R8 y+ }5 x6 k& h5 ?" q! i* O# ?+ }$ Z3 c3 O6 @: T# o0 w% ~, m* \, m Y( \ d3 V1 E+ c9 Q$ h }' \* v' S0 I/ K0 {4 t4 S! [' I7 v5 }; ?8 [/ [: u( ^- j. p) T$ j- V( f; q* L2 ]+ R: U% Z2 P0 u' J" T$ I' V' _" _9 \1 v4 ?6 ? j; `! e# `' @0 C# V# J) z7 k6 ]! f' a. d8 I6 [6 y1 a& {& l' Q0 B4 a9 m0 W( R; M- W9 j( Y/ ]- t$ B$ h( ^: o! m/ X P K; v! I$ k5 U7 {, ]$ n* Y, n8 u% O8 S q$ ^. S# D0 W. q' o% P8 {( X% W4 h8 P" A( j7 ?4 R+ [; ?3 z4 v9 V! V' y; @& F3 A K* i& B9 v) t. Z% R R1 j" v, ]8 r; ^. d) k2 _, d d1 M1 S+ B5 ?2 d+ c: Y8 y7 Z3 U* B* I# {% J' u D# H0 \& ^3 V+ W5 p8 H" A- {' \" \7 A+ N. v9 C( c7 w V7 @. M+ ?4 A' a% D% Q2 S/ ]9 ^* I3 [' V# I. ^3 m! G* b9 U1 N/ a* A9 ]* O+ k, n: u! ~/ R+ N. o5 {, w. C; z0 ~0 k' ?9 B6 p& o. |/ u8 G, U$ z% b; M/ y) {4 P. o0 a6 ?, H9 J. h! g: z2 d: a, ^4 I& e$ W$ e. M+ D# L( G9 R0 R5 Z+ g1 r# |3 i* }4 u) S: K- A7 i2 }& m) |5 | c% s* `0 r& r. E8 i+ v0 h+ v8 W. c9 N6 V: \% P; T1 p% E0 F7 R' q1 K. }2 I4 c0 }8 ?, O4 g5 W! A+ R1 Z# Y' C: i# K 表2 校验用样本 序号 切削速度 m/min 进给速度 mm/min 切削深度 mm 切削时间 min 实测切削力 N 初始磨损量 mm 实际磨损量 mm 1 21.98 30 3.5 3 406.42 0 0.08 2 21.98 30 3.5 5 594.89 0.21 0.232 3 21.98 30 2.5 3 305.75 0.094 0.105 4 15.70 35 3 2.5 630.21 0.105 0.12 5 18.84 25 3 2.5 392.95 0.06 0.077 6 18.84 35 2.5 4.5 616.82 0.31 0.344 7 15.70 35 3 4.5 721.26 0.224 0.238 8 15.70 30 2.5 4.5 522.47 0.18 0.216 9 15.70 25 3.5 2 574.40 0 0.049 10 21.98 35 3 4.5 591.28 0.304 0.238 11 15.70 25 3 7.5 560 0.11 0.16 12 18.84 35 3.5 7 725.38 0.186 0.235

# G* a& T# L. \3 O/ ^4 u* X. N9 M' N. R% A9 x) h; R. K6 Z7 g* A; H! C, H: Z1 T1 O" v, k7 w5 l$ w- b" Q! T9 t, p- z1 X: ], D) W" R @5 k$ R8 L" }3 g$ c; \$ {% j( z" @+ ]- a, [+ R( c/ n5 [( t6 m8 S+ [& |/ e" G8 y6 P7 w- y$ S# G1 A& s# t5 X0 [4 k9 s2 I* `. O% X: k* o6 E H) l( j8 D e) x! `0 A8 l. j, w/ d% C7 z6 \$ n8 _8 _! X; `5 |- x3 k5 P. ?5 }2 F; w' }9 E# c* A) u# I, B1 ^0 f0 g, L1 a& S' |9 e* X9 t! B# G, T1 ~; j2 R3 ~9 F* {, R& W2 j Q& j. j& P ~3 a/ r$ G0 y% d' A! g! L2 w/ q ?) E! v9 c/ y9 ?9 d6 r, l: o6 L1 R+ L6 |/ J9 Y& {- J: o! R5 D4 L( B, ~0 Z. O8 v9 A0 q3 H0 w9 m8 e6 ~# ]# A* @& x* @( M: ~+ w: i" a% W: z' S/ @7 `) l( _/ m! h- i7 V8 J3 W8 K( D( T# U( ~% N- j! Y1 \9 c D9 h; b1 v: a5 o% t [- Z, C% y/ v% J4 |4 M1 o, w& k' {; `9 M+ Q, L# t! I; x9 M4 z# m/ i. N( }$ _9 T, X9 H8 L5 u; `0 |5 ~; t2 l7 e: y. {5 G+ N4 W; q* Z& ^% Y! z1 E2 Y D: F# n* c( i9 W4 a) {0 @: l% x% h5 i5 W! p& H% b5 v4 M+ m5 I8 l/ p O5 [- X; c! X, f/ h! |- r9 y; }; Z5 A& `; _* \9 k9 u; u- l* x1 ]; i: ~/ B3 r; T/ U x2 m; C4 s Q" @' J6 x! R9 f5 g o/ A! e& c1 N, D: X* `. U( u- p- w2 L+ {, R5 D! L0 X( R5 l5 J: h( P5 a% t% t0 u! M! B3 t$ B7 \7 ^0 u" W4 ]5 }5 R; t, }- j/ r- U- @1 b/ c4 q+ h, T# I8 x# G8 d( u# ~# Y. i4 A7 M( G# f0 |1 J5 B" E/ w! ]+ F; L0 _5 V$ p) D; H$ y2 d0 z, k8 B) ~, R" \# k* U, T, G6 @: {/ ]( x9 j# C p5 _- F/ G8 h9 W7 `7 @* a: m* Y9 a, Y/ I, s2 L/ j- H/ e1 ~5 z" A8 R- H7 `3 e4 U6 S2 E4 y' x* g" m. ^9 u9 C5 @2 V: k. s7 z) A D' Y8 C2 S6 y4 R- f0 G; ~2 k+ b7 J# D/ o* w/ y) t# H9 @5 x( `! B/ ]6 X( L6 z" \2 D% P3 _0 R# d# Q* w) I0 Y: _: P. g* i( z( k- ~3 v" {" T7 `' ?/ n/ E1 W7 L8 _- a. w8 c9 _# R$ g9 a2 o" y- k: L6 b1 m. g$ B P8 y- C1 B7 n r9 ]" a+ H8 {" v6 M8 S- n- E. E! z+ S+ e+ m8 O8 Q* W7 {& ]# f+ ?# U; y( o9 k1 ?5 U+ z" I0 Z* Y& m! Z2 @% {0 b0 S5 x, S' @% l9 J( T7 J& i* `! N( Q0 [- \* s% r! g* O8 a7 e* O' y: D/ c6 h, G0 @7 x: M; F* l8 Q, }, ? T1 w: B/ ^' `9 g0 R/ g# [ 表3 切削力模型的校验结果 序号 实测磨损量mm 实测磨损类别 隶属度 估测磨损类别 估测磨损量 mm µA µB µC µD µE 1 0.08 A 1 0 0 0 0 A 0.1 2 0.232 C 0 0.449 0.551 0 0 C 0.227 3 0.105 B 1 0 0 0 0 A 0.1 4 0.12 B 0.547 0.453 0 0 0 B 0.123 5 0.077 A 1 0 0 0 0 A 0.1 6 0.344 E 0 0 0 0 1 E 0.3 7 0.238 C 0 0 0.722 0.278 0 D 0.276 8 0.216 C 0 0.732 0.268 0 0 C 0.213 9 0.049 A 1 0 0 0 0 A 0.1 10 0.38 E 0 0 0 0 1 E 0.3 11 0.16 B 0 1 0 0 0 B 0.15 12 0.235 C 0 0.323 0.677 0 0 C 0.234 类别识别率 84.6%

* \6 l4 D8 b) n8 o$ q. [9 u* ^4 v0 U, H) [* b! q# M6 Q' m2 l

图6 切削力模型校验效果图 ) |! v+ L* K) P

图7 时间模型校验效果图

4 结论

本文主要从实用的角度出发，建立了基于切削力和基于相对切削时间的两种磨损模型。主要结论如下： : B2 a) q' `1 k. M

. r/ i2 v6 G" F- L% [% F6 k/ {# p
1. 利用模糊分类和回归分析所建立的铣削力模型，真实地反映了切削参数与切削力的关系，较为理想地对刀具磨损状态进行分类和识别。校验结果表明，模型对刀具状态的有效识别率为84.6%,在工程上有应用价值。 ) N3 j! {4 M. m0 @" p
2. 探讨了基于时间的磨损建模方法。利用回归分析，直接建立刀具磨损量与切削参数和时间参数之间的关系，校验结果表明，该模型在一定程度上可满足工程上的要求。

1 W% r& q0 i9 I3 a4 T