CN102566505A - 一种数控机床的智能故障诊断方法 - Google Patents

Abstract

本发明提供了一种数控机床的智能故障诊断方法，包括如下步骤：（1）确定故障特征参数集、故障类型集和样本集；（2）根据故障特征参数集和故障类型集，建立支持向量机模型；（3）利用样本集，训练并优化支持向量机；（4）检测数控机床的多个故障特征参数，将多个故障特征参数作为故障诊断的数据源；（5）将所述数据源输入到优化后的支持向量机，得出故障诊断结果。本发明根据数控机床部件故障的发生机理及特点，当数控机床部件发生故障时，通过检测数控机床，对检测数据进行处理后输入最小二乘支持向量机，得到诊断结果。本发明实现了对数控机床快速、有效的故障诊断，同时其故障诊断能力强、确诊率高、误判率低。

Description

一种数控机床的智能故障诊断方法

技术领域

本发明属于故障诊断领域，具体涉及一种数控机床的智能故障诊断方法。

背景技术

目前，数控机床在现代制造业中的应用日趋普及，但是数控机床在使用中经常出现故障，使连续生产系统无法正常工作。若设备出现故障而未能及时发现和排除，结果不仅导致机床本身损坏，甚至可能造成机毁人亡的严重后果。因此，在实际生产中用较短的时间和较低的维修成本，对数控机床所发生的故障及时做出准确判断，查明故障部位，找出故障原因和排除方法，能大大减少维修的盲目性，它对提高生产经济性和安全性具有重要的意义。

数控机床是机电一体化紧密结合的典范，是一个庞大的系统，涉及机、电、液、气、电子、光等各项技术，其结构复杂，部件多样。目前故障诊断基本上还是采用如下几种传统的方法： 

（1）直观检查法。直观检查法是不借助任何仪器仪表就能诊断出故障的一种简便有效方法。问：问现场人员故障产生的过程、故障现象及后果，询问故障是渐发性的还是突发性的等内容；看：看机床各部分是否处于正常工作状态，电控装置有无报警信息，熔丝是否熔断，元器件烟熏烧焦，电容器膨胀、开裂，电子元件有无断脚、虚焊等；看工件的表面粗糙度、颜色、伤痕等；听：听异常响声和机床运转声；闻：闻电气元件焦糊味及其它异味等；摸：在整机断电条件下用手感来判别机床的故障，如温升、振动、伤痕、波纹、爬行、松紧等。

（2）综合诊断法。仪器检查法：利用万用表等检测仪器仪表，对故障疑点进行电流、电压等测量，将测量值与正常值进行比较，从中找寻故障所在位置；接口状态检查法：现代数控系统多将PLC内置其中，而CNC与PLC之间以接口信号相互通讯。有些故障与接口信号错误或丢失相关，有的接口信号可在相应的接口板和 输入/输出板上有指示灯显示，有的可通过简单操作在CRT上显示，或用PLC编程器观察接口信号状态来诊断故障；参数检查法：数控系统、PLC、伺服驱动系统都设置了许多可修改的参数，以适应不同机床、不同工作状态的要求，这些参数是保证机床正常运行的前提条件。一旦因电池电量不足或受到外界干扰，可能导致部分参数丢失或变化，使机床无法正常工作。通过核对、调整参数，有时可以迅速排除故障。

（3）部件替换法。随着现代技术的发展，电路的集成度越来越高，按常规方法很难查找故障，现代诊断数控机床故障方法中，越来越多的采用部件替换法。部件替换法就是在故障范围大致确认，外部条件完全正确的情况下，利用同样的印制电路板、集成电路芯片或元器件替换有疑点部分，把故障范围缩小到印刷电路板或芯片一级的方法。替换部件可以是系统的备件，也可用机床上相同部件。

（4）交叉换位法。当发现了故障板或疑是故障板，又没备件时，可将系统中两个相同或相兼容的板互换，从中判断出故障板或故障部位。注意硬件接线的正确交换和相应的参数交换。

（5）原理分析法。从系统的工作原理上分析各点的电平和参数，用万用表等仪器仪表对其进行测量和比较，进而对故障进行系统检查的一种方法。

（6）敲击法。当系统故障表现为时而正常时而不正常时，基本可以断定为元器件接触不良或焊点虚焊，用敲击法敲击到虚焊或接触不良部位时，故障就会重现。

（7）信号跟踪法。按照控制系统框图从前往后或从后往前逐一跟踪有关信号的有无、大小及不同运行方式下的状态，与正常情况进行比较。

针对传统方法故障诊断成本高，速度慢的缺点，目前有很多专家开展了一系列故障诊断方法和技术的研究。典型的研究成果有：

（1）FTA分析法：它是对复杂系统进行安全性、可靠性分析的一种好方法，用它来预测和诊断故障，在数控机床的运行和维修中能收到良好的效果。

（2）CBR推理法：CBR 技术已经发展成为一套解决实际问题的方法论。CBR 故障诊断的基本思想是当寻找诊断方案时， 在过去类似诊断方案成功案例基础之上进行推理， 通过类比和修改来完成当前故障诊断任务。CBR 技术的特点使其在数控机床故障诊断领域有着独特的优势。 

（3）进化Monte Carlo优化法：是一种基于总体的蒙特卡洛方法， 它将遗传算法中的变异和交叉运算引入到模拟回火的框架，是一种具有“学习”能力的马尔可夫蒙特卡洛方法。

但是，以上几种故障诊断方法仍然不能从根本上解决数控机床故障诊断的快速性和准确性。

发明内容

本发明针对上述现有技术的不足，提供了一种具有快速性和准确性的数控机床智能故障诊断方法。

本发明是通过如下技术方案实现的：

一种数控机床的智能故障诊断方法，包括如下步骤：

（1）确定故障特征参数集、故障类型集和样本集，所述故障特征参数集中包含多个故障特征参数，所述故障类型集中包含多个故障类型；

（2）根据故障特征参数集和故障类型集，建立最小二乘支持向量机模型；

（3）利用样本集，训练并优化最小二乘支持向量机；

（4）检测数控机床的多个故障特征参数，将多个故障特征参数作为故障诊断的数据源；

（5）将所述数据源输入到优化后的最小二乘支持向量机，得出故障诊断结果。

进一步的，所述故障特征参数集包括数控机床的温度、检测元件的输出量、高低压电流、振动频率、噪声强度和伺服控制的输出电压；所述故障类型集包括伺服单元，加减速时间设定，检测元件，可控硅，旋转编码器，负载量，传动链润滑，可编程控制器或微处理器，检测控制系统或主轴驱动系统，继电器，熔断器，电源变压器。

进一步的，所述最小二乘支持向量机模型的输出向量为故障类型集，输入向量为故障特征参数集中每个故障特征参数的主元特征值。

进一步的，步骤（3）中采用不完全交叉验证网格搜索法在线优化最小二乘支持向量机。

进一步的，步骤（5）包括以下子步骤：

（5.1）将数据源进行小波去噪处理；

（5.2）计算数据源中每个故障特征参数的主元特征值；

（5.3）将步骤（5.2）中计算出的主元特征值输入到优化后的最小二乘支持向量机，得出故障诊断结果。

进一步的，所述主元特征值的计算方法为： 

（a）计算故障特征参数                                                

的信息熵

；其中，

为一个随机过程中

出现的概率，

，；

（b）归一化处理得到

；其中

；

（c）建立协方差矩阵

；其中

；

（d）计算协方差矩阵

的特征值

和特征向量

； 

（e）计算主元特征值

。

本发明根据数控机床部件故障的发生机理及特点，当数控机床部件发生故障时，通过检测数控机床的高低压电流、温度、振动频率、噪声强度和伺服控制的输出电压得到故障信息，对故障信息数据进行处理后输入最小二乘支持向量机，得到诊断结果。本发明将最小二乘支持向量机、熵带分析和主元特征提取相结合，实现了对数控机床快速、有效的故障诊断，同时其故障诊断能力强、确诊率高、误判率低。

附图说明

图1为数控机床的智能故障诊断方法的流程图；

图2为数控机床故障快速诊断方法优选实施例的流程图；

图3为数控机床故障快速诊断过程示意图。

具体实施方式

下面结合附图和具体实施方式对本发明做进一步详细的说明。

如图1-3所示，本发明提供了一种数控机床的智能故障诊断方法，包括如下步骤：

（1）确定故障特征参数集、故障类型集和样本集；

根据数控机床部件故障的发生机理及特点，本发明将数控机床的温度、检测元件的输出量，高低压电流、振动频率、噪声强度和伺服控制输出电压这六个参数作为故障特征参数。本发明采用向量R表示故障特征参数集，则，其中为六个故障特征参数。

根据数控机床常见的故障特征，本发明将数控机床的伺服单元、加减速时间设定、检测元件、可控硅、旋转编码器、负载量、传动链润滑、可编程控制器或微处理器、检测控制系统或主轴驱动系统、继电器、熔断器和电源变压器这12个参数作为故障类型。本发明采用向量F表示故障类型集，则

，其中

为12个故障类型。

收集数控机床部件运行状态的历史数据，作为样本集。为了清晰表达故障特征参数与故障类型参数之间的关系，可以按照表1中的形式建立故障数据库。  

表１

（2）根据故障特征参数集和故障类型集，建立最小二乘支持向量机模型；

本发明选取Gauss 分布的径向基函数（RBF）作为最小二乘支持向量机的核函数，因为它能够很好地满足Mercer条件，能够尽量准确地反映训练样本数据的分布特征（对称分布），能够替代内积，通过一个非线性映射，在一个高维特征空间中给出一个最优分类超平面。RBF核函数的表达式如下：

，其中

是径向基核函数的参数。

本发明中最小二乘支持向量机模型的输出向量为故障类型集

。为了降低支持向量机模型训练的难度，本发明对故障特征参数集进行预处理后，作为支持向量机模型的输入向量。所述预处理是指对故障特征参数

的熵带进行样本归一化和主元特征提取。具体的讲，所述预处理包括以下步骤：

（a）熵计算：由于收集或采集到的数控机床故障数据受到各种环境因素的影响，有些数据呈现出不确定性。这种情况严重混乱了故障类别间的可分辨性，增大了故障分类与辨识的难度。本发明采用熵使故障的误判最小化。熵计算的过程如下：

　  设

为数控机床一个故障判别参数（

）。在一个随机过程中，如果

出现的概率

，且

，那么信息熵

为：

　　　　　　　　　

上式表明，若越大, 表明

中所蕴涵的信息量越丰富, 则

的不确定程度也越大。

（b）归一化处理：考虑到各种特征参数数值的差异性及分散性，为了降低它们之间由于量值差异过大造成的影响，对熵通过线性变换映射到[０，１]范围内的相对含量，以降低各特征参数数据之间的互斥性。归一化处理方法如下：

                   

（c）对通过线性变换的熵

进行如下变换（令：）：

　　　　　　　　　　　　　　　

式中，

为

维的正交矩阵，其中矩阵

中的第列

是协方差矩阵的第

个特征向量。

　　　　　　　　　　　　　　

即

　　　　　　　       

，

为协方差矩阵

的特征值，

为相应的特征向量，则可以得到主元

。

　　　　　　　　　　　　　　

，。

（3）利用样本集，训练并优化支持向量机；

　  本发明在样本集中随机抽取600组故障数据，分别求取熵，用最大的熵减去最小的熵得到熵带

。在熵带范围内取300组数控机床故障数据作为故障样本，对最小二乘支持向量机进行训练和优化。

训练支持向量机的过程如下：设给定N个训练样本

是输入数据，是输出数据。函数估计问题可以描述求解下面问题：

                   

                    

                  

其中，是核空间映射函数，权矢量，误差变量

，

是偏差量。损失函数

是SSE误差和规则化量之和，

是控制对错分样本惩罚的程度的可调参数，它控制对超出误差样本的惩罚的程度，实现在训练误差和模型复杂度之间的折衷。与标准SVM不同，这里只有等式约束，且优化目标的损失函数是误差

的二范数。核空间映射函数的目的是从原始空间中抽取特征，将原始空间中的样本映射为高维特征空间中的一个向量，以解决原始空间中线性不可分的问题。引入Lagrange函数，把有约束优化问题转化成无约束优化问题。

    

即：

  

其中，拉格朗日乘子

，

。对上式进行优化求解，根据KKT最优条件，把上式分别对

和

求偏微分并令它们等于0，即：

                     

                      

                        

                         

得到：

              

对于

消去

和

，得到如下方程：

                        

其中

   

即得到如下方程组：

              

        

优化支持向量机的过程如下：采用不完全交叉验证网格搜索法确定最小二乘支持向量机的最佳参数对(

)。具体的讲，包括如下子步骤：

(a1) 确定合适的正则化参数集和核参数集。实验发现，按照指数增长方式生成两种参数集是一种有效的方法，例如，

，

，网格搜索简单直接，因为每一个参数对(

)是独立的，可以并行地进行网格搜索。

(a2)应用网格搜索法在正则化参数集和核参数集中选择一个参数对(

)，用该参数对进行交叉验证。首先用一个步长为的(

) 组合, 得到学习精度最高的

 和

 的值。然后在这两个值旁的一定范围内进行一次更细致的网格搜索。在两个值的一定范围内把样本集D分为

组

，把任意的

组作为训练集，剩余的一组作为验证集。通过选择不同的验证集，可重复次。泛化性能可通过下式评价：

    

式中：

是第组验证集；

是验证集的样本；

是用

作为训练样本时得到的参数向量；

是LS-SVM的输出。

 (a3)循环选择参数对进行交叉验证，计算每个参数对的

，直到网格搜索停止。使得

最小的参数对(

)是最佳的，应用不完全交叉验证方法选择参数能够避免过拟合问题。

本发明通过样本集的学习，完成训练过程寻找到支持向量集，确定拉格朗日乘子

和偏移因子

，得到优化后的支持向量机。　

（4）检测所述故障特征参数集中的故障特征参数，作为故障诊断的数据源，即检测数控机床的温度、检测元件的输出量、高低压电流、振动频率、噪声强度和伺服控制的输出电压等故障特征参数作为数据源。

（5）将所述数据源输入到优化后的最小二乘支持向量机，得出故障诊断结果。具体的讲，包括以下子步骤：

（5.1）将数据源进行小波去噪处理；

数控机床特征参数在采集和传输过程中，不可避免会受到大量噪声信号的影响而产生杂波，必然会影响对故障判断的准确度，因此本发明采用小波去噪技术对数据源进行小波去噪处理。小波去噪处理的算法如下:

令检测特征量为：

式中

为有用信号；

为噪声序列。假设是零均值且服从高斯分布的随即序列，即服从分布。对上式两边做小波变换，有：

　　　　　　　　　　　

式中

表示小波基对

做时移为

，尺度因子为

的小波变换。经正交小波变换后，最大程度地去除了

的相关性，其能量集中在少数的小波系数上，而这些系数即是在各个尺度下的模极大值。

（5.2）计算数据源中每个故障特征参数的主元特征值，其中计算主元特征值的方法与步骤（3）中的方法相同，在此不再赘述。

（5.3）将步骤（5.2）中计算出的主元特征值输入到优化后的最小二乘支持向量机，得出故障诊断结果，即正常或有故障，以及故障所属的故障类型。

采用本发明所述的智能故障诊断方法对浙江乐清液压制造有限公司的数控车床进行故障诊断，结果如表2所示。

表2

故障类型	故障诊断的识别率
		伺服单元故障	98%
加减速时间设定太大或太小	98%
		检测元件故障	99%
可控硅烧毁	100%
		旋转编码器故障	100%
负载超量	99%
		传动链润滑不良	98%
可编程控制器或微处理器损坏	100%
		检测控制系统或主轴驱动系统故障	100%
继电器损坏	100%
		熔断器损坏	100%
电源变压器损坏	100%

可见，本发明能够实现对数控机床有效的故障诊断。

本发明不仅局限于上述具体实施方式，本领域一般技术人员根据本发明公开的内容，可以采用其它多种具体实施方式实施本发明，因此，凡是采用本发明的设计结构和思路，做一些简单的变化或更改的设计，都落入本发明保护的范围。

Claims (6)

1.一种数控机床的智能故障诊断方法，包括如下步骤：

（1）确定故障特征参数集、故障类型集和样本集，所述故障特征参数集中包含多个故障特征参数，所述故障类型集中包含多个故障类型；

（2）根据故障特征参数集和故障类型集，建立最小二乘支持向量机模型；

（3）利用样本集，训练并优化最小二乘支持向量机；

（4）检测数控机床的多个故障特征参数，将多个故障特征参数作为故障诊断的数据源；

（5）将所述数据源输入到优化后的最小二乘支持向量机，得出故障诊断结果。

2.根据权利要求1所述的智能故障诊断方法，其特征在于，所述故障特征参数集包括数控机床的温度、检测元件的输出量、高低压电流、振动频率、噪声强度和伺服控制的输出电压；所述故障类型集包括伺服单元、加减速时间设定、检测元件、可控硅、旋转编码器、负载量、传动链润滑、可编程控制器或微处理器、检测控制系统或主轴驱动系统、继电器、熔断器和电源变压器。

3.根据权利要求1所述的智能故障诊断方法，其特征在于，所述最小二乘支持向量机模型的输出向量为故障类型集，输入向量为故障特征参数集中每个故障特征参数的主元特征值。

4.根据权利要求1所述的智能故障诊断方法，其特征在于，步骤（3）中采用不完全交叉验证网格搜索法在线优化最小二乘支持向量机。

5.根据权利要求1所述的智能故障诊断方法，其特征在于，步骤（5）包括以下子步骤：

（5.1）将数据源进行小波去噪处理；

（5.2）计算数据源中每个故障特征参数的主元特征值；

（5.3）将步骤（2）中计算出的主元特征值输入到优化后的最小二乘支持向量机，得出故障诊断结果。

6.根据权利要求3或5所述的智能故障诊断方法，其特征在于，所述主元特征值的计算方法为： 

（a）计算故障特征参数                                               的信息熵

；其中，

为一个随机过程中出现的概率，

，

；

（b）归一化处理得到

；其中

；

（c）建立协方差矩阵；其中；

（d）计算协方差矩阵

的特征值

和特征向量

； 

（e）计算主元特征值。

Images