一种基于改进D-S证据理论的设备故障诊断方法与流程

本发明涉及一种故障诊断方法，尤其是涉及一种基于改进d-s证据理论的设备故障诊断方法。
背景技术：
：在设备运行过程中，其振动信号中包含着系统丰富的运行信息，因此振动信号常被作为设备故障监测和诊断的主要依据。由于一些设备结构复杂，其振动信号的采集需要多个传感器采集，信息量大，再加上各个传感器不同，测点不同，因此会导致诊断结果的不稳定性。要实现对设备故障的准确判断，需要利用信息融合处理不同信息源传递过来的数据，按照一定规则组合这些数据，然后做出全面、高效、准确合理的判断，增强整个故障诊断系统的可靠性。相比于其他不确定推理方法，d-s证据理论不仅可以处理因模糊性和随机性产生的不确定性，还可以把主客观因素都考虑进去，将其应用于设备的故障诊断可以减少因运行环境等不确定因素造成的误判情况的发生。然而，经典d-s证据理论只在证据冲突较小的时候效果好，当证据间存在高度冲突时，其合成结果常会产生与直觉相悖的结果，出现一票否决的现象，缺乏鲁棒性。针对这些问题，许多国内外学者提出了一些改进方法，主要包括修改组合规则和修正证据源两个方面。在修改组合规则方面，大多数的改进方法(《李弼程,王波,魏俊,钱曾波,黄玉琪.一种有效的证据理论合成公式[j].数据采集与处理,2002,(01):33-36.》、《孙全,叶秀清,顾伟康.一种新的基于证据理论的合成公式[j].电子学报,2000,(08):117-119.》、《邓勇,施文康.一种改进的证据推理组合规则[j].上海交通大学学报,2003,(08):1275-1278.》)都在yager公式基础上进行改进，提出新合成公式，但是由于去掉了原公式里的归一化因子，将绝大部分的支持率赋给了不确定项，提高了决策难度，仍会出现“一票否决”的情况。在修正证据源方面，一些改进方法(《杜峰,施文康,邓勇.证据特征提取及其在证据理论改进中的应用[j].上海交通大学学报,2004,(s1):164-168.》、《高锦秋，刘景林.ds证据理论的改进及在航空交流发电机故障诊断中的应用[j/ol].微特电机，2016,44(06):37-40.》、《费翔,周健.一种处理冲突证据的d-s证据权重计算方法[j].计算机工程,2016,42(02):142-145.》)利用证据间的相似度、冲突度或矛盾系数等来确定修改证据源不同的权重因子，对证据进行贴近度处理，这些方法都有一定的改进效果，但是没有利用好证据间的关联性，权重的确定存在一定的主观性。技术实现要素：本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于改进d-s证据理论的设备故障诊断方法。本发明的目的可以通过以下技术方案来实现：一种基于改进d-s证据理论的设备故障诊断方法，包括：步骤s1：获取设备运行过程中釆集的振动信号，并提取的时域分析特征参量和时频分析特征参量；步骤s2：基于预构建的rbf、bp和elman神经网络模型，将提取的时域分析特征参量和时频分析特征参量作为三个神经网络模型的输入量，进行初步诊断；步骤s3：将三个神经网络模型的输出结果进行融合；步骤s4：根据融合后的结果诊断故障。所述时域分析特征参量包括：高频波顶因数、波顶因数、峭度、高频峰值、峰值、加速度有效值包络线、加速度有效值、高频有效值、速度有效值和偏斜度；所述时频分析特征参量包括：速度有效值包络线和iso10816。所述步骤s2具体包括：步骤s21：将时域分析特征参量作为rbf神经网络模型的输入，得到第一初步诊断结果；步骤s22：将速度有效值包络线通过三层小波包分解得到8个频段信号并作为bp神经网络模型的输入，得到第二初步诊断结果；步骤s23：将iso10816通过三层小波包分解得到8个频段信号并作为elman神经网络模型的输入，得到第三初步诊断结果。所述步骤s22和步骤s23中的频段信号均通过归一化处理后作为神经网络模型的输入。所述步骤s3具体包括：步骤s31：对第一初步诊断结果、第二初步诊断结果和第三初步诊断结果进行归一化处理后作为原证据体；步骤s32：计算三个原证据体之间的jousselme距离和距离相似度；步骤s33：基于得到的jousselme距离和相似度，对原证据体进行融合得到新证据体。所述步骤s33中得到的新证据体为：其中：m(a)为新证据体，a为所有可能的故障种类的集合，k为新证据体的冲突系数，dsi为新证据相似度占新证据和原证据证据相似度之和的比重，f(a)为保留的一组原证据体，m1(ai)为第一个原证据体中发生第i种故障类型的概率，m2(aj)为第二个原证据体中发生第j种故障类型的概率，m3(al)为第三个原证据体中发生第l种故障类型的概率。所述新证据体的冲突系数具体为：所述步骤s4具体包括：步骤s41：提取所有故障类型中概率最大的一个，并判断其是否超过设定阈值，若为是，则执行步骤s42；步骤s42：判断该故障的概率是否超过相比其他故障类型的概率是否超过设定差值，若为是，则执行步骤s43；步骤s43：输出概率最大的故障类型作为诊断结果。所述设定差值为0.2。所述设定阈值为0.5。与现有技术相比，本发明具有以下有益效果：1)计算权重修正证据源，在保留原组合规则优点基础上对组合规则进行改进的算法，通过相关数据和信息，充分利用计算机的高速运算和多源信息的互补性来提高信息的质量。2)本发明的算法能够充分利用各个信息，避免了传统d-s证据理论悖论现象的出现，对提高设备诊断结果的精确性和可靠性有重要作用。附图说明图1为本发明的结构示意图。具体实施方式下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。一种基于改进d-s证据理论的设备故障诊断方法，如图1所示，包括：步骤s1：获取设备运行过程中釆集的振动信号，并提取的时域分析特征参量和时频分析特征参量，其中，时域分析特征参量包括：高频波顶因数、波顶因数、峭度、高频峰值、峰值、加速度有效值包络线、加速度有效值、高频有效值、速度有效值和偏斜度；时频分析特征参量包括：速度有效值包络线和iso10816，其中，iso10816振动监测评估标准。在设备运行过程中釆集振动信号时，由于运转设备存在多个振动源，而且振动信号在传输过程中会受到传输通道特性的影响，所以采集到的时域信号会呈现出杂乱无章的特点。此外，在釆集设备运行过程中的振动信号时，通常会包含很多的噪声，因此振动信号的频域内容就包含很多复杂的频率成分。为了能够有效利用采集到的时域信号，识别出与设备故障相关的频率及其特征，本发明采用基于时域和是时频域相结合的提取方法对故障信号进行特征提取。本发明用到的时域分析的特征参量为：高频波顶因数、波顶因数、峭度、高频峰值、峰值、加速度有效值包络线、加速度有效值、高频有效值、速度有效值、偏斜度；本发明用到的时频分析特征参量为：速度有效值包络线，iso10816。步骤s2：基于预构建的rbf、bp和elman神经网络模型，将提取的时域分析特征参量和时频分析特征参量作为三个神经网络模型的输入量，进行初步诊断，具体包括：步骤s21：将时域分析特征参量作为rbf神经网络模型的输入，得到第一初步诊断结果；步骤s22：将速度有效值包络线通过三层小波包分解得到8个频段信号并作为bp神经网络模型的输入，得到第二初步诊断结果；步骤s23：将iso10816通过三层小波包分解得到8个频段信号并作为elman神经网络模型的输入，得到第三初步诊断结果。本发明采用三种神经网络对设备振动特征进行初步诊断。这里将分别采用时域特征和时频特征对设备振动故障进行诊断。由于rbf神经网络诊断稳定，且精确度高，本发明采用此神经网络作为时域特征诊断网络。本发明对速度包络线有效值和iso10816两种振动特征采用小波包分解的时频分析方法得到信号的时频特征，再分别采用bp和elman神经网络分别对这两种信号的时频特征进行初步诊断。以上三种神经网络的网络结构建立完毕后，利用该三种不同的网络对设备振动故障进行故障初步诊断。将设备振动故障特征向量的训练样本通过训练得到合适的rbf、bp和elman神经网络模型，之后将测试样本分别经三种神经网络进行故障诊断。步骤s3：将三个神经网络模型的输出结果进行融合；由于bp神经网络和rbf、elman神经网络在对设备进行故障诊断时会出现结果不一致现象，因此本发明利用改进d-s证据理论将对这三种网络进行融合诊断，使之最终得到一个确定的结果。首先构建证据体，利用神经网络初步诊断的结果经归一化处理得到基本概率分配函数，在确定设备故障的基本概率分配函数后，根据以下步骤进行新证据构造和融合。其融合过程如下：a.计算证据体之间的jousselme距离和距离相似度。其中，jousselme距离由jousselem等人定义，具体为：其中，d表示2n×2n矩阵；n表示辨识框架的元素个数。d中的元素d定义为：d(b,c)＝|b∩c|/|b∪c|。b.计算新证据体矩阵：将计算出来的证据距离构成相似度矩阵，表示所有证据之间的相似度，反映了证据的冲突度，以矩阵形式作为修正对应焦元的权重系数，可以更全面的利用证据间的关联度，使新证据更加稳定可靠。得到的新证据体矩阵里的每一个元素经过归一化后得到每个新的证据焦元。c.计算新证据体与原证据体之间的平方根距离：用平方根距离可以简单直观地表示出新证据修正的变化及新老证据之间的偏差。由于将原来的证据进行修正后，会将证据平均化，原生证据无论有多大冲突，多组证据中必定有一组是正确的，如果保留原有的一组正确的证据，会进一步提高融合效果。本发明基于这个思想，根据修正后的证据与原生证据之间平方根距离选择一组原生优秀证据，距离越小，说明此证据越合理，保留最合理的证据给f(a)。d.改进组合规则。新组合规则表示为：其中，f(a)表示保留的一组优秀原生证据，为新证据体的冲突系数，dsi表示新证据的权重系数，新证据相似度占新证据和原证据证据相似度之和的比重:其中，ds(a,b)表示证据a与证据b之间的距离相似度。步骤s4：根据融合后的结果诊断故障。步骤s22和步骤s23中的频段信号均通过归一化处理后作为神经网络模型的输入。步骤s4具体包括：步骤s41：提取所有故障类型中概率最大的一个，并判断其是否超过设定阈值，若为是，优选的，设定阈值为0.5，则执行步骤s42；步骤s42：判断该故障的概率是否超过相比其他故障类型的概率是否超过设定差值，优选的，设定差值为0.2，若为是，则执行步骤s43；步骤s43：输出概率最大的故障类型作为诊断结果。下面结合具体应用对本申请进行说明，1.信号特征提取过程通过振动测试分析系统采集到的部分数据如表1所示。表1采集的振动特征数据前十种特征参量将直接用于时域分析，后两种特征参量(速度有效值包络线，iso10816)分别采用小波包分解用于时频分析，用其分解信号的均方根值表示信号能量的大小。本文选用“db1”小波作为基函数，并进行三层小波包分解。对设备采集到的每组振动信号进行小波包分解后可以得到8个频段信号，其频带分别为：(0～93.75hz)、(93.75～187.5hz)、(187.5～281.25hz)、(281.25～375hz)、(375～468.75hz)、(468.75～562.5hz)、(562.5～656.25hz)、(656.25～750hz)。如表2，表3分别为速度有效值包络线和iso10816两种振动信号的能量特征向量，其中，故障1表示轴承损坏，故障2表示机械不平衡，故障3表示周期性冲击。表2速度有效值包络线信号能量特征向量工况状态e1e2e3e4e5e6e7e8故障12.1104-0.00360.0261-0.1019-0.00010.0099-0.00990.0001故障23.0906-0.1044-0.1384-0.04440.06290.0114-0.0556-0.0041故障32.2299-0.1159-0.0524-0.0056-0.0269-0.03110.00030.0206表3iso10816信号能量特征向量2.设备初步诊断部分在本申请的rbf网络模型中，其输入层的10个输入量分别对应10个时域特征参量归一化后的特征向量；在bp神经网络模型中，输入层的8个输入量分别对应速度包络线有效值经小波变换后归一化后的特征向量；在elman反馈型神经网络中，其输入层的8个输入量分别对应iso10816经小波变换后归一化后的特征向量。三种神经网络的网络结构建立完毕后，利用该三种不同的网络对设备振动故障进行故障初步诊断。其中，故障类型1用(1，0，0)表示轴承损坏，故障类型2用(0，1，0)表示机械不平衡，故障类型3用(0，0，1)表示周期性冲击。上述表1中数据1的诊断结果如表4所示。表4神经网络诊断结果3.基于改进d-s证据理论的信息融合先构建证据体，假设设备的故障状态识别框架包括故障一(a1)、故障二(a2)、故障三(a3)，上小节的神经网络结果经归一化处理，得到基本概率分配函数m1、m2和m3,如表5所示。表5故障三状态下的证据体mk(ai)a1a2a3m10.04830.88300.0687m20.14700.12750.7255m30.13630.16500.6988其中ai(i＝1,2,3)分别代表设备的故障类型一、二、三；代表经bp神经网络诊断的结果，代表elman诊断结果，代表rbf网络诊断的结果；故障类型1用(1,0,0)表示，故障类型2用(0,1,0)表示，故障类型3用(0,0,1)表示。将本文改进算法和传统d-s证据理论算法作了诊断效果比较，如表6所示。表6本文改进算法的融合结果对比m123(a1)m123(a2)m123(a3)m123(u)d-s证据理论0.01780.34170.64050.0000本文改进方法0.07280.16210.76510.0000采用本文提出的改进算法，虽然证据一高度支持目标a2，但是后面两条证据均支持目标a3，从最终融合结果来看比较合理，符合实际情况，诊断效果比传统方法提高了很多。4.设备故障的决策诊断从表6的融合结果来看，传统证据理论融合和本文改进方法的结果的输出值最大分别为0.6405,0.7651，都高度支持目标a3，满足决策规则的三个条件，因此，诊断出设备的故障状态为故障三。本发明针对传统d-s证据理论对设备进行振动故障诊断时会出现一些悖论问题，提出了d-s理论改进算法。首先从时域、频域方面对振动信号进行特征提取，在频域中用到的小波包分解的方法使信号有效避免了频率混叠，便于频率特征向量的提取。然后通过rbf、bp和elman神经网络对设备振动特征信号进行初步诊断，利用初步诊断结果构造基本概率函数。之后利用相似度矩阵确定基本概率分配函数的权重系数，构造新证据体。考虑到新证据体可能会偏保守，丢失了原证据的优点，本发明又根据新证据与原证据的证据间相似度确定新证据与原证据的权重分配，使得新证据源与原证据以不同的权重一起通过改进的融合规则进行融合。新的改进方法既保留了传统证据理论的优点，又解决了冲突问题，有广泛的适用范围。通过实验表明此方法不仅正确诊断出设备运行系统的全部故障，而且诊断的精确度高于传统d-s和其他改进d-s理论的诊断方法，实现了对设备振动故障的精确诊断。当前第1页12