基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法

1.本发明涉及故障诊断技术领域，尤其涉及基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法。


背景技术：

2.随着现代科学技术的飞速发展与智能工业技术的迅速进步，当前制造业的机械设备日益呈现自动化、大型化、精密化和智能化，这不仅对机械设备关键部件的结构设计更加复杂，制造精度要求更高，还对机械设备中关键部件的安装和工作中的运行平稳性提出了更高的要求。滚动轴承是几乎所有类型的旋转机械的重要组成部分。滚动轴承作为机械设备中的关键部件，30％的机械故障是由轴承故障引起的。滚动轴承有助于机器平稳高速地旋转，减小摩擦，并能承载大量负载，是支撑旋转轴的关键部件。及时、准确地判断轴承的故障形式和故障发生位置，对提高机械设备安全性、可靠性和使用效率都有重要意义。
3.基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法，主要通过选择和考察合适的信号动态分析指标，通过分析振动信号的混合时域特征，结合图论的相关算法分析数据特征间的相互联系，以此来对轴承故障进行诊断。


技术实现要素：

4.本发明的目的是为了解决现有技术中存在现有分析数据特征方法对滚动轴承故障诊断结果精度不高的缺点，而提出的一种基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法。
5.一种基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法，包括以下步骤：
6.s1、滚动轴承的振动数据经过提取后，对滚动轴承振动信号的时域无量纲特征进行提取，包括峭度cq、峰值x
p
、脉冲因子cf、裕度因子ce，以及小波包分解的8维能量，其中，在对小波包能量提取时，采用小波包变换法通过对数据不断分解，同时选取分解后的低频部分和高频部分的能量特征。
7.s2、利用主成分分析法提取小波包分解能量中贡献度占95％的能量，经过分析，提取前5维的小波包能量，结合选取的4类时域无量纲多特征构成9维混合故障特征向量。
8.s3、分别计算每维混合故障特征向量中数据之间的欧氏距离，通过图论算法，使用计算出的数据间的欧式距离建立无向图模型，无向图模型通过邻接矩阵表示，利用奇异值分解法(svd)对每维特征的邻接矩阵进行分解，提取每个邻接矩阵的奇异值作为分类器输入向量。
9.s4、使用支持向量机(svm)作为故障分类器，建立故障分类规则，将经过奇异值分解处理的输入向量以一定比例分为训练集和测试集。利用训练集训练svm模型，并设置合适的svm分类器的惩罚因子和核参数，建立基于训练集的支持向量机(svm)轴承故障诊断模型。
10.s5、将测试集作为输入，使用训练后的svm模型进行故障分类，测试本算法训练的
模型分类准确度，并输出结果。
11.优选的，根据振动信号时域特性对提取的振动数据进行分析，分别提取信号的峭度cq、峰值x
p
、脉冲因子cf、裕度因子ce进行数据多方面的无量纲特性分析。
12.优选的，将数据的无量纲时域信号和数据经过pca算法提取的小波包能量组成为混合特征向量，实现对数据的多角度分析。
13.优选的，根据图论算法的理论核心，利用欧氏距离算法计算混合数据特征向量各数据之间的欧式距离，将数据作为点，将欧式距离作为边，无差别进行点对点连接，组合成为无向图，通过邻接矩阵的方式展现图形结构关系，从而对数据点之间的相关性进行分析。
14.优选的，使用奇异值算法计算无向图邻接矩阵的奇异值，作为分类器的输入数据特征向量。
15.优选的，使用svm算法建立故障诊断模型，建立故障分类规则，实现通过数据特征向量的输入对滚动轴承的故障信号进行有效地分析与诊断。
16.本发明的有益效果
17.1、在对数据进行特征提取时，选取了多个时域特性无量纲指标和小波包分解能量组成数据混合特征。无量纲数据指标对早期故障较为敏感，受工作条件变化的影响较小，具有更好的数据分析效果。
18.2、小波包能量通过对数据进行小波包分解得到各分解频带信号，将各分解频带信号所携能量与信号总能量之比作为特征向量以表征设备的运行状态。信号经小波包分解后所得分解频带相互独立且无冗余。混合数据特征不仅能够表征数据的各项信息，表示数据的故障变化，还能够从多角度对数据进行故障分析，提升了数据故障诊断的精准性。
19.3、为了分析故障数据的变化，利用欧式距离算法计算每个数据之间的距离，不仅可以表示数据的变化，还可以表征数据之间的相关关系。结合数据点与欧氏距离建立数据的无向图结构，更能代表数据的特点，使故障分解更具全面性，进一步提升数据特征的精度。
20.4、对数据进行上述处理后，利用奇异值算法提取图结构的奇异值作为故障分类器的输入向量，既保留了原始数据的特性，又转变了数据形式，更适合分类器的输入。使用支持向量机(svm)模型作为故障分类器，选择高斯核函数作为支持向量机的核函数，更适合当前的输入向量。利用选定好的不同类型故障数据对svm模型进行训练，提高故障诊断结果的准确度，使分类结果更接近真实故障类型，并减小误差的存在。
附图说明
21.图1为基于k-means分类原则的pls改进算法的流程图
22.图2为本发明运行过程流程图
具体实施方式
23.下面将结合本发明实施例中的附图，对本发明中涉及的方法更明白、更完整地描述。显然，下面描述中的实施例子，仅仅是本发明的一部分实施例子，而不是全部的实施例子。
24.实施例
25.针对现有利用数据特征进行滚动轴承故障诊断的问题，目前存在使用单一数据特征作为特征向量进行故障诊断精度不高，故障类型分类出现错分、漏分等缺点，基于此，提出一种基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法。
26.参照图1，一种基于混合特征的无向图邻接矩阵的滚动轴承故障诊断方法，基于目前由于直接采用滚动轴承的振动数据作为输入向量进行故障诊断时，数据的特征不够明显，使模型精度受到影响。并且，当数据中存在异常数据时，会严重影响模型的准确度，并且该异常数据的诊断结果误差很大，影响操作人员对滚动轴承故障发生的判断，造成严重的损失。因此，建立一种可以利用分析数据混合特征的改进算法：
27.首先对采集的滚动轴承振动数据进行特征提取。在实际采集中，振动数据包含了很多有用的信息，这些信息都会对故障的精准诊断产生影响，使用没有经过特征提取的原始振动数据很难准确诊断出故障的类型，为了避免这类问题的发生，采用分析数据时域信号的方法，其中有量纲特征指标对于早期故障较为敏感，但会因工作条件的变化而变化，不能稳定的反映轴承状态，而且有量纲时域诊断参数其值在实际应用中很难进行区分。而无量纲指标对轴承故障较敏感，不易受周围环境的影响，因此提取数据时域信号的无量纲特征参数。在选择无量纲特征参数上选取了峭度cq、峰值x
p
、脉冲因子cf、裕度因子ce这四个参数，峭度对冲击信号特别敏感，特别适用于表面损伤类故障、尤其是早期故障的诊断；轴承故障时的峰值明显大于正常状态时的峰值；脉冲因子是用来检测信号中有无冲击的指标；裕度因子常用来检测机械设备的磨损状况。计算公式如下：
[0028][0029]
x
p
＝e(max|xi|)
[0030][0031][0032]
其中脉冲因子cf为一个周期内信号最高值和最低值之间差的值，即峰值x
p
与一周期内信号的平均幅值的比值。裕度因子ce为峰值x
p
与一周期内信号的方根幅值xr的比值。
[0033]
它描述了信号值的变化范围的大小与平均幅值表示为了提取数据更多振动信号的更多数据特征，对振动数据使用haar小波包进行3层小波包分解，提取了数据的8维小波包能量，利用pca主成分分析算法对8维小波包能量进行降维，按照95％的贡献度提取了前4维的小波包能量作为使用。将上述步骤提取到的4类数据混合参数与降维后的小波包能量进行结合，成为8维数据混合特征向量。
[0034]
为了分析每一维数据之间的相关关系，使用欧式距离算法计算每一维中各个数据之间的欧氏距离，用来体现数据点之间的相似性，欧氏距离是最常见的用来衡量两组高维数据在空间上距离的方法，公式如下：
[0035][0036]
其中,x＝[x1,x2,
…
,xn]，y＝[y1,y2,
…
,yn]分别表示高维空间中的两组数据。
[0037]
参照图2，图谱理论的基本思想是在图和矩阵之间建立对应的关系，通过矩阵的相关属性来研究图相关的问题。基于图论的方法是把聚类转换为一个组合优化问题并利用图论和相关的启发式算法来解决该问题，本发明利用图论算法的思想核心，将数据点看作图的顶点，欧氏距离看作图的连接边，这样可以将数据之间的相关关系用无向图的结构表示出来，更有利于故障特征的分析。表示数据之间相关关系的数据利用邻接矩阵的形式表示。为了方便输入，利用奇异值算法将邻接矩阵转换成向量的形式。
[0038]
最后，经过上述步骤处理后的数据特征向量包含了大量的数据信息，用此去训练基于svm(支持向量机)的故障诊断模型，实现对待测数据的定性分析和定量计算，以达到准确诊断出轴承故障类型的目的
[0039][0040]
表1.故障诊断不同类型诊断率
[0041]
在上表1的测试结果中可以看出，利用本发明的方法诊断的轴承故障类型分为四类，其中，轴承的正常情况、内圈故障与外圈故障的诊断率均为100％，滚动体故障出现一个样本错分的情况，诊断准确率为99.54％。
[0042]
以上所述仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以同等替换或改变，都应涵盖在本发明的保护范围之内。