一种大型数控铣床刀具磨损程度识别方法
【专利摘要】本发明提供了一种大型数控铣床刀具磨损程度识别方法，该识别方法通过大型数控铣床运行状态下振动时域信号的获取，采用快速傅里叶变换得到振动信号的频域分布，选取时域与频域中的多个统计特征参数进行基于扩散映射法的刀具磨损特征参数的降维，采用留一交叉验证法与优化搜索算法确定扩散映射的尺度参数选取问题；结合Nystrom扩展与核回归算法对未知待测刀具进行磨损程度的识别。本发明可以有效克服大型数控铣床刀具磨损样本缺失之缺点，提高大型数控铣床刀具磨损程度的识别精度，减少因刀具磨损识别不及时造成的维修成本与时间。
【专利说明】一种大型数控铣床刀具磨损程度识别方法
【技术领域】
[0001] 本发明属于大型数控铣床领域，具体涉及一种大型数控铣床刀具磨损程度的识别 方法。
【背景技术】
[0002] 数控铣床除能铣削平面、沟槽、轮齿、螺纹和花键轴外，还能加工比较复杂的型面， 生产效率较高，在机械制造行业中被广泛应用。特别是大型数控铣床（如龙门铣床），因其 加工精度和生产效率均较高，常应用于大型工件的成批生产。
[0003]刀具作为大型数控铣床最易损伤的部件，对其进行及时有效的状态识别与监测尤 为重要。据统计，刀具磨损是引起铣床故障的首要因素，由此引起的停机时间占铣床总停 机时间的20% -30%。而且，在铣削加工中，一旦刀具发生损伤故障而没有及时发现，会直 接影响产品加工质量、加工精度和生产效率，严重的还将导致铣床功能失效及停机、工件报 废、甚至危害人员安全。因此，如何有效地识别大型数控铣床刀具的磨损程度，监测刀具运 行状态，已成为大型数控铣床智能化发展急需解决的问题。
[0004] 近年来，国内外学者在铣床刀具磨损程度的识别上做了大量的研究工作，提出了 诸多有效的高精度、高可靠性诊断方法，如时间序列分析、频谱分析、小波分析、神经网络、 支持向量机、混合智能等，这为大型数控铣床刀具磨损程度的识别提供了一定的技术基础。 然而，对于大型数控铣床而言，其刀具磨损程度的识别还面临着如下几个问题：（1)大型数 控铣床的加工对象一般较大，使得不同刀具磨损程度下的数据采集比较困难，成本较高，训 练样本数据少；(2)大多数研究方法需要人为确定关键参数的取值（如支持向量机的惩罚 函数和尺度因子），主观性较强，在没有较多经验可供参考的大型数控铣床刀具磨损程度的 识别上误诊率较高。
【发明内容】
[0005] 本发明针对上述现有技术的不足，提供了一种准确性更高的大型数控铣床刀具磨 损程度识别方法。
[0006] 本发明是通过如下技术方案实现的：
[0007] -种大型数控铣床刀具磨损程度识别方法，包括以下步骤：
[0008] (1)采集大型数控铣床在Μ种刀具磨损状态下的振动时域信号；从每种磨损状态 的振动时域信号中截取连续的采样数为t的不重叠的S组信号；并利用快速傅里叶变换将 每组时域信号的波形转换成频域分布；其中M、t和S均为大于1的正整数；
[0009] (2)分别计算每种磨损状态下S组信号数据的8个时频域无量纲统计特征参数， 组合成样本数据集{(C^nii)} (i = 1，2，…，Ν,Ν = MXS)，其中 Q = {cn, ci2, ···，ci8}为第 i 个样本的特征参数集，m,为Q对应的刀具后刀面最大磨损量；
[0010] 其中，所述8个时频域无量纲统计特征参数，包括4个时域无量纲统计参数 (Cn、Ci2、Ci3、Ci4)和4个频域无量纲统计参数（C i5、Ci6、Ci7、Ci8);设第i个样本信号数据 为：Xi = {xn, xi2,. . .，xin},经过FFT变换后的频域信号数据为fi = {fn, fi2,..., fj， 广 _ ^i-rms /-f 一 k=\ ~^^ ， Li2= ， Ci3 - "TT ? ， x, X. S n - \ i ι-rms i " f(A'-/；-,")· A Σ(Χ,,-Χ,)4 Cii= -~^-， ci4=-^--3，Ci5 ，C =-^, 1 Yfik S： n-\ N^, A, I h nil 广 _ k;ni4 · 一 n/4 ， 免二 1
[0011] xik为样本点Xi的第k个信号，n为样本点Xi的数据量，七= -hl·， n k=l ' =J丄?χ2，xi-m = max ixik I k = 1，…，n}，Ss = Χ(χα -X,·)2，pik 为频率 fik nil ( 、 nil ( 、 的功率谱，Ph = max {Pik I k = 1，-,n}，Nt_m = --，Nt_0 = -， 1Σ{λ2·^-} I ?Α- I k 二 ? 丨 Jt=l Ρ = Ι^νΣ^2;
[0012] (3)结合留一交叉验证建立尺度参数ε的数学优化模型，根据样本数据集Ω，采 用优化搜索算法对尺度参数ε进行寻优，找出使得目标函数取值最小的ε取值；
[0013] (4)根据扩散映射（DM)方法的降维原理，构造样本数据集的邻接矩阵W = {w.j} nxn :wij - exP (_ I I Ci-Cj 11 / ε ); 2 6
[0014] 其中|C;-C;| =Σ(&)2为&和(^的欧氏距离；然后，对W按行进行标准化处 ￥=1 Ν 理，令 w' = {w1 iji ,w'i^wu /(Σ^); i=i
[0015] (5)求解W'的特征值和特征向量:W' φ,= λ,φ,，XkSW'的第k个最大特 征值，且有ι= Κ> λι> λ2>…，小，为特征值λ,对应的单位特征向量；根据预先设 定的降维维数Κ选取特征值：Λ = { λ 1； λ 2, ···, λ κ}，对应的特征向量构成映射矩阵Ω = { Φ Φ 2，···，小 iJ NXK ;
[0016] (6)计算各样本点在映射矩阵Ω下的映射坐标： I TV 0
[0017] φ (Q) = {￠,(^),k = 1，2，...，K}， Α(Φ=^7^ΣεχΡ(-|C/-C7.|「/f) n * Λ j=\ Φ--?Α)为Q在映射Ω下的第k个分量，Φ jk为<K的第j个分量；
[0018] (7)每隔固定时间间隔采集一次大型数控铣床在运行状态下的振动时域信号，构 成待诊断信号X ;并将时域波形转换成频域分布；然后，计算待诊断信号数据的8个时频域 统计特征参数C(X) = {CX1，CX2，…，CX8}；
[0019] (8)对X进行Nystrom扩展，计算X在映射矩阵Ω下的映射坐标： γ Ν 0
[0020] φ ⑴=. ^exp(-||C(X)-Cy|| I sy<f>jk 丄、.凡k J=l
[0021] (9)对X进行核回归分析，得到对应的磨损程度值，即后刀面最大磨损量估计值 为：
[0022] mx = X[exp(-||C(X) - Ct f !ε)· ] j[exp(-|C(X) - Cq『/ ε)。 /=1 / q=\
[0023] 优选的，步骤（3)中采用禁忌搜索算法对尺度参数ε进行寻优，具体包括如下步 骤：
[0024] (3. 1)确定包含尺度参数ε的核函数为：
[0025] Κε(€η€\) = exp(-|c,. -q|2 Ιε)/^，其中 Q 和 C」为样本点；
[0026] (3. 2)采用留一交叉验证法，计算各样本点的磨损程度估计值：
[0027] ,咖)=hoy tnc") /=1. }Φ1 f ρ=\^ρφ?
[0028] (3. 3)建立关于尺度参数ε的数学优化模型： 1 Ν
[0029] Μιη Ζ(6·) =--ih^a·))2 , ε e (0,1) Ν i=l
[0030] (3. 4)采用禁忌搜索对上述数学模型进行优化求解，得出使得预测误差最小的尺 度参数ε取值。
[0031] 本发明具有如下有益效果：
[0032] (1)目前对数控铣床刀具的磨损研究，集中于对中小型铣床的刀具磨损研究，因为 中小型铣床刀具磨损样本数据比较容易收集，而对于大型数控铣床刀具的磨损识别研究很 少。同时，现有的大多数故障诊断方法是在样本数据量大的前提下开展的，在小样本情形 下，这些方法的训练效果很差，对刀具的磨损识别无能为力。本发明的提出可以克服上述弊 端，本发明能够小样本情形下有效识别大型数控铣床刀具的磨损程度。
[0033] (2)目前，对刀具磨损程度的识别研究大多只考虑对三种磨损状态（初期磨损、 中度磨损和严重磨损）的分类上，没有研究刀具磨损程度的渐进式非线性演变规律，本发 明通过建立刀具磨损程度的回归模型，可有效地揭示大型数控铣床刀具磨损程度的演变规 律。
[0034] (3)本发明可以有效克服大型数控铣床刀具磨损样本缺失之缺点，提高大型数控 铣床刀具磨损程度的识别精度，减少因刀具磨损识别不及时造成的维修成本与时间。
【具体实施方式】
[0035] 本发明提供了一种大型数控铣床刀具磨损程度识别方法，包括以下步骤：
[0036] (1)采集大型数控铣床（一般情况下多为龙门铣床）在Μ种刀具磨损状态下的振 动时域信号；
[0037] 其中，Μ根据刀具后刀面最大磨损量来确定。本实施例中取Μ = 5，其对应的5种 刀具磨损状态，根据后刀面最大磨损量划分的不同分为正常状态、轻微磨损、中度磨损、较 大磨损、急剧磨损，如表1所示。
[0038] 表1后刀面最大磨损量与道具磨损状态对应表
[0039]
【权利要求】
1. 一种大型数控铣床刀具磨损程度识别方法，其特征在于，包括以下步骤： (1)采集大型数控铣床在Μ种刀具磨损状态下的振动时域信号；从每种磨损状态的振 动时域信号中截取连续的采样数为t的不重叠的S组信号；并利用快速傅里叶变换将每组 时域信号的波形转换成频域分布；其中M、t和S均为大于1的正整数； (2)分别计算每种磨损状态下S组信号数据的8个时频域无量纲统计特征参数，组合成 样本数据集{(Ci; mj } (i = 1，2，…，Ν, Ν = MX S),其中 Q = {cn, ci2,…，ci8}为第 i 个样 本的特征参数集，m,为Q对应的刀具后刀面最大磨损量； 其中，所述8个时频域无量纲统计特征参数,包括4个时域无量纲统计参数（Cn、 Ci2、Ci3、Ci4)和4个频域无量纲统计参数（Ci5、C i6、Ci7、Ci8);设第i个样本信号数据 为：Xi = {xn, xi2,. . .，xin},经过FFT变换后的频域信号数据为fi = {fn, fi2,. . .，fj， 1 Σ(-ν/?- --^/)5 ，。手，,_? ， 叉 i i-rms ^ i n Συ?-ιι,η)·ρ? 1 Σ(χ?- _X<)4N p Cn = - ^72 ， ^ m Pi-rms I nil -是二《/4 . - "7m ; 是二 1 xik为样本点\的第k个信号，n为样本点\的数据量，巧=丄文?，U , n \n k=l Xi-m = max {xik I k = 1，…，n}，St = ?-^-Υ(x;1 -xj2，Pik 为频率 fik 的功率谱，= Vn-lfc=i nil { . nil ( 、 Ek2·及丨 _ max{Pik|k= 1，…，η}，#,- ,,: ^-- , ^,-0= -， 1 ^p2 I p =\yJ-Pik ; (3)结合留一交叉验证建立尺度参数ε的数学优化模型，根据样本数据集Ω，采用优 化搜索算法对尺度参数ε进行寻优，找出使得目标函数取值最小的ε取值； (4)根据扩散映射（DM)方法的降维原理，构造样本数据集的邻接矩阵W= {?^·}ΝΧΝ =exp (-1 I Ci-Cj I2/ ε ); . 6 其中卜―c I =Σ(~-?)2为&和(^的欧氏距离；然后，对w按行进行标准化处理， ^=1 Ν 令 r = {w1 ij}, Κ- = wtJ /(Σ wu) ?' i=i (5)求解W'的特征值和特征向量Φ, = λ,φ,，的第k个最大特征 值，且有1= λ0 > > λ2>…，（^为特征值λ ,对应的单位特征向量；根据预先设定 的降维维数Κ选取特征值：Λ = { λ 1； λ 2, ···, λ κ}，对应的特征向量构成映射矩阵Ω = { Φ Φ 2，···，小 iJ NXK ; (6)计算各样本点在映射矩阵Ω下的映射坐标： j Ν ? Φ (Ci) = k = 1，2，…，κ}，llq-Cij7 幻為々 7V . Ak J=] Φ--?Α)为Q在映射Ω下的第k个分量，ΦΛ为的第j个分量； (7)每隔固定时间间隔采集一次大型数控铣床在运行状态下的振动时域信号，构成待 诊断信号X ;并将时域波形转换成频域分布；然后，计算待诊断信号数据的8个时频域统计 特征参数 C(X) = {CX1，CX2，…，CX8}； (8)对X进行Nystrom扩展,计算X在映射矩阵Ω下的映射坐标： 1 7V 0 Φ ⑴={ Φ k (X), k = 1,2, ···, Κ}, A (^) - , ^exp(-||C(Z)-Cy I /ε) φβ Λ/ · Ak y=1 (9)对X进行核回归分析，得到对应的磨损程度值，即后刀面最大磨损量估计值为： mx = X[exp(-||C(Z) - C, f Iε)·nh]/|>p(-|c(JQ - Cqf /ε)。 /=1 / q=\
2.根据权利要求1所述的大型数控铣床刀具磨损程度识别方法，其特征在于，步骤（3) 中采用禁忌搜索算法对尺度参数ε进行寻优，具体包括如下步骤： (3. 1)确定包含尺度参数ε的核函数为： Κε(€(,(^) = exp(-||c; -/^)/72^，其中 Q 和 C」为样本点； (3. 2)采用留一交叉验证法，计算各样本点的磨损程度估计值： /?,⑷=；|χ.(?,?;)·/?]/ χ^(?,?,) j=lj^i j p=Up 淖 i (3. 3)建立关于尺度参数ε的数学优化模型： 1 /V 、 Min Ζ{ε) = -V (mi -ηι^ε))' ， ε e (0,1) Nti (3. 4)采用禁忌搜索对上述数学模型进行优化求解，得出使得预测误差最小的尺度参 数ε取值。
【文档编号】G06F17/50GK104050340SQ201410321361
【公开日】2014年9月17日 申请日期:2014年7月7日 优先权日:2014年7月7日
【发明者】周余庆, 李峰平, 梁薇薇, 郑静 申请人:温州大学