CN113569478A - 基于元素影响度和置信规则库的旋转机械故障诊断方法 - Google Patents
基于元素影响度和置信规则库的旋转机械故障诊断方法 Download PDFInfo
- Publication number
- CN113569478A CN113569478A CN202110847686.8A CN202110847686A CN113569478A CN 113569478 A CN113569478 A CN 113569478A CN 202110847686 A CN202110847686 A CN 202110847686A CN 113569478 A CN113569478 A CN 113569478A
- Authority
- CN
- China
- Prior art keywords
- fault
- confidence
- rule base
- degree
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000012216 screening Methods 0.000 claims description 29
- 230000008859 change Effects 0.000 claims description 19
- 238000005457 optimization Methods 0.000 claims description 17
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 14
- 230000004927 fusion Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 10
- 238000013528 artificial neural network Methods 0.000 claims description 8
- 230000004913 activation Effects 0.000 claims description 5
- 238000011478 gradient descent method Methods 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 13
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000002159 abnormal effect Effects 0.000 abstract description 3
- 238000013507 mapping Methods 0.000 abstract description 3
- 238000012423 maintenance Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000036541 health Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/211—Selection of the most significant subset of features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Evolutionary Biology (AREA)
- Acoustics & Sound (AREA)
- Medical Informatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提出基于元素影响度和置信规则库的旋转机械故障诊断方法,针对旋转机械异常振动特征与其发生故障种类之间呈现的复杂非线性映射关系,以及将所有特征数据用于建模时,大量对诊断结果影响度低的信息造成建模计算量大且实时性差等问题,故利用平均影响度(MIV)算法挑选并加权故障特征为模型输入,以旋转机械的故障类型为输出,构建置信规则库(BRB)模型,规则库中的多条或某条置信规则将被输入的振动特征数据不同程度的激活,结合证据推理算法融合被激活规则,从融合的后项置信结构中得出发生故障类型,同时给出潜在故障类型发生概率,并且模型参数物理意义明确,可解释性强,最后以旋转机械中基本的电机转子为例,验证所提方法的有效性。
Description
技术领域
本发明涉及机械故障诊断技术领域,具体为基于元素影响度和置信规则库的旋转机械故障诊断方法。
背景技术
随着旋转机械运行速度越来越高、负载越来越大、自动化程度不断加强,加之其复杂的工业环境和长时间持续运行,设备一旦发生故障或任何损伤,将会直接影响其在工业应用中的运行,造成严重的经济损失和事故发生。因此,准确诊断出旋转机械设备故障,是保障其持续、稳定、安全运行的必要手段。旋转机械故障诊断的本质是通过诊断旋转机械中的转子系统、齿轮箱等异常辨别工况,由于旋转机械在高负载、长时间运行环境下,引起其发生故障因素较多,加之获取故障信号的传感器易受本身精度差异和环境噪声干扰,使得对其故障诊断具有不确定性、强随机性等特点,从而导致故障信号与故障类型之间存在很强的非线性关系。
实际工程中旋转机械经常工作在高负载、长时间运行等环境下,易出现转子不平衡、转子不对中、基座松动、连接器松动、齿轮缺齿等故障,并引发旋转机械不同程度的非线性振动。旋转机械异常振动特征与其发生故障种类之间非线性映射关系复杂;将所有特征数据用于建模时,大量对诊断结果影响度低的信息造成建模计算量大且实时性差。
发明内容
本发明提出了基于元素影响度和置信规则库的旋转机械故障诊断方法,可以监测机器情况,指导维修和保证机器正常运转。
本发明的技术方案如下:
基于元素影响度和置信规则库的旋转机械故障诊断方法,首先进行设备安装,包括在电机柔性转子试验台的转自支撑部位的不同方向上,安装多个频域振动加速度传感器,该方法包括,
步骤一,建立筛选模型,通过所述频域震动加速度传感器得到故障特征因子,并对其进行加权分配比重;
步骤二,分类故障类型,得到根据步骤一中所得模型,设置参考点,建立置信规则库,采用独热编码对其进行描述;
步骤三,设置置信规则库,筛选出故障特征的样本数据,以及筛选故障特征前项参考点和后项参考点,基于专家知识和历史数据构建初始置信规则库;
步骤四,执行故障诊断,根据步骤三所得的置信规则库,得到诊断结果;
步骤五,模型优化,基于独热编码和欧式距离进行模型参数优化。
作为进一步优化,设定电机转子的故障类型包括齿轮缺齿、连接器松动、基座松动、转子不对中、转子不平衡和正常运行,所述设备安装的准备步骤包括安装频域震动加速度传感器的个数为M。
作为进一步优化,计算机获取1X-3X倍频的频域幅值作为模型输入,在每个电机转子的故障类型下获取3个倍频的3×M个故障特征(c1,c2,…,c3M),
所述步骤一包括,根据与每种故障类型的相关性,从3×M个故障特征(c1,c2,…,c3M)中挑选相应的故障因子,并对所挑选的故障银子进行加权分配比重,从而得到筛选模型。
作为进一步优化,所述步骤一中,建立如下模型
Q(ci,δi)=q(c1,c2,…,cn,F),
其中,Q(ci,δi)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。
作为进一步优化,所述步骤一包括,
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新,定义故障特征的损失函数如下:
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
其中L表示故障特征因子数,n表示样本组数;
神经网络的拟合输出如下所示:
IV=[IV1,IV2,…,IVn]T,
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
根据故障特征因子对故障的平均影响程度,得出其相对应的约简加权信息Q(ci,δi)。
作为进一步优化,所述步骤二包括,
根据故障特征对故障诊断的平均影响程度筛选出相应的故障特征cI,cII,cIII,并根据筛选出的故障特征平均影响度值对其进行加权,加权属性权重为δI,δII,δIII,根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
对以上数据采用独热编码对其进行描述。
作为进一步优化,所述步骤四包括,
S4-1,计算输入量的匹配度,
S4-2,计算被激活规则的权重,
S4-3,融合被激活的规则并输出诊断故障结果。
作为进一步优化,所述S4-1包括,
所述S4-2包括,
所述S4-3包括,
O(X)={(Fj,βj),j=1,2,…,6},
其中
作为进一步优化,所述步骤五包括,
S5-1,对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,将系统输出与实际输出的欧式距离为参数设计优化目标函数
S5-2,以
为约束条件得到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性权重δk)为最小值时的最优参数集合P。
作为进一步优化,所述步骤五还包括,采用随机梯度下降法来对模型进行训练,获得ξ(P)的最优参数集,从多功能电机柔性转子试验台的M个频域振动加速度传感器获取到的3M个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
本发明的工作原理及有益效果为:
通过对非线性振动数据建模,及时诊断旋转机械发生故障类型,对提高旋转机械维护的“状态修”水平和实现健康周期管理起到辅助决策的作用。其优势在于利用监测振动数据根据专家经验知识建立的模型对旋转机械故障进行诊断,在不进行停机维护的情况下判断旋转机械的健康状态,从而自适应调整修护周期,保障旋转机械系统稳定、安全运行。通过多源传感器对旋转机械运行情况进行监测,利用与故障类型具有高相关性的监测数据建模诊断故障类型,对指导旋转机械维修和保证旋转机械健康状态有着积极作用。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明方法的流程框图;
图2是本发明方法实施例中ZHS-2型多功能电机柔性转子试验台所采集3个倍频15个故障特征对应故障类型的平均影响度;
图3是本发明方法实施例中通过MIV值筛选出故障特征c3,c5,c13的样本数据及其对应的故障类型;
图4是本发明方法实施例中模型优化前后对应的均方根误差图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都涉及本发明保护的范围。
具体实施例1,如说明书附图1所示,本发明提出的基于元素影响度和置信规则库的旋转机械故障诊断方法,利用平均影响度(MIV)算法挑选并加权故障特征为模型输入,以旋转机械的故障类型为输出,构建置信规则库(BRB)模型。规则库中的多条或某条置信规则将被输入的振动特征数据不同程度的激活,结合证据推理算法融合被激活规则,从融合的后项置信结构中得出发生故障类型,同时给出潜在故障类型发生概率,并且模型参数物理意义明确,可解释性强。最后以旋转机械中基本的电机转子为例,验证所提方法的有效性。
具体为:
准备步骤:设定在ZHS-2型多功能电机柔性转子试验台转子支撑部位的不同方向上,安装M个频域振动加速度传感器感知电机转子故障状态。由于不同电机转子故障类型会引起多频率成分的振动幅值变化,故获取1X~3X倍频的频域幅值变化为模型输入。设置电机转速为1500r/m,则基频1X为25Hz,n倍频nX(n=1,2,3,…)为(n×25)Hz。在每个电机转子故障类型下,获取3个倍频的3×M个故障特征(c1,c2,…,c3M)。
步骤一:置信规则库输入信息维数不能过大,否则会造成组合规则“爆炸”,严重影响其故障诊断的实时性,但同时故障特征信息过少又不足以表征故障类型。为此利用平均影响度算法将准备步骤中M个频域振动加速度传感器获取到的3×M个故障特征(c1,c2,…,c3M)中挑选出与故障类型高相关性的故障特征因子,并对其进行加权分配比重,将其用于建模诊断故障,将提高模型的确诊率和实时性。建立筛选模型如下
Q(ci,δi)=q(c1,c2,…,cn,F), (1)
其中,Q(ci,δi)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。具体步骤如下:
S1-1利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新。定义故障特征的损失函数如下:
S1-2对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
其中L表示故障特征因子数,n表示样本组数。神经网络的拟合输出如下所示:
IV=[IV1,IV2,…,IVn]T, (8)
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
根据故障特征因子对故障的平均影响程度,可以得出其相对应的约简加权信息Q(ci,δi)。
步骤二:根据步骤一中利用平均影响度筛选出对故障诊断影响程度高的故障特征cI,cII,cIII,并根据筛选出的故障特征平均影响度值对其进行加权,加权属性权重为δI,δII,δIII。根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
对于置信规则库的后项属性,采用独热编码对其进行描述,旋转机械故障编码见表1。
表1.故障类型及其编码
步骤三:根据步骤二中对筛选故障特征前项参考点和后项属性的定义,置信规则就可以被表示为Rk:如果且且那么Y是{(D1,βk,1),…,(DN,βk,N)},k=1,2,…,每条规则权重θk的初始值都设为1,属性权重初始值根据MIV值设置为δI,δII,δIII。
部分置信规则如表2所示
表2.初始置信规则库的部分规则
步骤四:在步骤三中基于专家知识和历史数据构建好初始置信规则库后,具体证据推理故障诊断过程如下:
S4-1,计算输入量的匹配度。将筛选故障特征的样本数据作为输入,输入变量xi与前项参考值存在一定的匹配度 如果xi小于等于中的最小值,或者大于等于中的最大值时,则xi与的匹配程度是0或1;否则xi与的匹配程度则为
O(X)={(Fj,βj),j=1,2,…,6}, (14)
其中
步骤五:基于独热编码和欧式距离的模型参数优化,具体步骤如下:
S5-1,由于旋转机械设备内部机理复杂,专家难以确定故障诊断模型中所有参数的精确值,为此需要将初始参数微调,提高模型的精度。但置信规则库推理的融合结果为发生各类故障的可信度,并不能直接度量模型诊断效果和优化模型参数,故对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,旋转机械故障编码见表1。
将系统输出与实际输出的欧式距离为参数设计优化目标函数
S5-2,利用S5-1给出的优化目标函数,找到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性0≤θk≤1,0≤δk≤1,0≤βj,k≤1和权重δk)为最小值时的最优参数集合P,约束条件为=1。采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)的最优参数集,从多功能电机柔性转子试验台的M个频域振动加速度传感器获取到的3M个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
具体实施例2,如说明书附图1-4所示,
在使用时,本发明专利以旋转机械中基本的电机转子为例,采用ZHS-2型多功能电机柔性转子系统,来验证所提方法的有效性。在该试验台转子支撑部位的不同方向上,安装5个频域振动加速度传感器感知电机转子故障状态。由于不同电机转子故障类型会引起多频率成分的振动幅值变化,故获取1X~3X倍频的频域幅值变化为模型输入。设置电机转速为1500r/m,则基频1X为25Hz,n倍频nX(n=1,2,3,…)为(n×25)Hz。在每个电机转子故障类型下,获取3个倍频的15个故障特征(c1,c2,…,c15),采集210次(其中100次数据用于模型训练,110次数据用于模型测试)。
步骤一,利用平均影响度算法将5个频域振动加速度传感器获取到的15个故障特征(c1,c2,…,c15)中挑选出与故障类型高相关性的故障特征因子,并对其进行加权分配比重,将其用于建模诊断故障,将提高模型的确诊率和实时性。建立筛选模型如下:
Q(ci,δi)=q(c1,c2,…,cn,F), (1)
其中,Q(ci,δi)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。具体步骤如下:
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新。定义15个故障特征的损失函数如下:
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
其中L=15表示故障特征因子数,n表示样本组数。神经网络的拟合输出如下所示:
S1-3,根据15个故障特征因子对故障的平均影响程度,可以得出其相对应的约简加权信息Q(ci,δi)。为了方便理解,设置网络的最大迭代次数为2000,期望目标误差最小值为1.0*e-5,增量δ=0.1。基于此训练稳定的网络,分别采用MIV算法求取c1-c15故障特征的平均影响值,并根据MIV值从大到小顺序排列。
其排序情况和计算结果如图2所示,选取c3,c5,c13为影响度高的故障特征,MIV值分别为1.1143、0.6876、0.4902。并根据公式6,计算选取故障特征的相对权重为0.4057,0.2504,0.1785。
步骤二,根据步骤一中利用平均影响度筛选出对故障诊断影响程度高的故障特征。根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
为了便于理解前项参考值的选定,根据步骤(2)中筛选出的高影响度故障特征c3,c5,c13样本数据600组,选取参考点如下:
对于置信规则库的后项属性,采用独热编码对其进行描述,旋转机械故障编码见表1。
步骤三,根据步骤二中对筛选故障特征前项参考点和后项属性的定义,置信规则就可以被表示为Rk:如果且且那么Y是{(D1,βk,1),…,(DN,βk,N)},k=1,2,…,每条规则权重θk的初始值都设为1,属性权重初始值根据MIV值设置为δI,δII,δIII。
为了便于理解,这里沿用步骤一和步骤二中筛选出的高影响度故障特征c3,c5,c13样本数据600组,以及筛选故障特征前项参考点和后项属性,可得到部分置信规则如表3所示:
表3.实施例2中初始置信规则库的部分规则
步骤四,在步骤三中基于专家知识和历史数据构建好初始置信规则库后,具体证据推理故障诊断过程如下:
S4-1,当第118组特征数据被输入时,输入量为X=[0.0776、0.0855、0.1380],则X1匹配参考点小正(PS)和正中(PM)的程度为0.4478和0.5522,X2匹配参考点正中(PM)和正大(ML)的程度为0.8773和0.1227,X3匹配小正(PS)和正中(PM)的程度为0.5585和0.4415,对于其他参考点的匹配程度均为0。
S4-2,利用公式11和12计算被激活规则的权重wk,对于第118组特征数据,可以得到其对R37-R38、R42-R43、R62-R63和R67-R68的激活权重分别为w37=0.2195,w38=0.1734,w42=0.0307,w43=0.0243,w62=0.2706,w63=0.2139,w67=0.0378,w68=0.0299而其他规则的激活权重均为0,也就是激活了8条规则。
S4-3,利用证据推理算法得到融合后的输出置信结构,可以分别通过公式15和16计算得出。将S4-2中关于第118组特征数据的及带入公式14,可以计算得出O(X)={(F1,0),(F2,0),(F3,0),(F4,0.0073),(F5,0.7968),(F6,0.1959)}。根据融合结果中对F1-F6的置信度,置信度最大的后项属性对应的故障类型即为诊断结果。对于第118组特征数据,其对应的故障诊断结果为连接器松动故障,潜在齿轮缺齿故障发生概率为19.59%,并利用公式17计算其基于独热编码的欧氏距离误差Δ=0.2032。
步骤五,基于独热编码和欧式距离的模型参数优化,具体步骤如下:
S5-1,由于旋转机械设备内部机理复杂,专家难以确定故障诊断模型中所有参数的精确值,为此需要将初始参数微调,提高模型的精度。但置信规则库推理的融合结果为发生各类故障的可信度,并不能直接度量模型诊断效果和优化模型参数,故对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,旋转机械故障编码见表1。
沿用步骤(5)中故障特征c3,c5,c13的600组样本数据作为训练数据,可以得到优化目标函数
S5-2,利用S5-1给出的优化目标函数,找到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性权重δk)为最小值时的最优参数集合P,约束条件为0≤θk≤1,0≤δk≤1,0≤βj,k≤1和(k=1,2,…125)。采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)<0.0016时的最优参数集,从多功能电机柔性转子试验台的5个频域振动加速度传感器获取到的15个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
具体实施例3
本发明方法的流程图如图1所示,核心部分是:针对旋转机械异常振动特征与其发生故障类型之间呈现的复杂非线性映射关系,以及将所有特征数据用于建模造成的计算量大且实时性差等问题。首先利用平均影响度(MIV)算法约简并加权重要特征信息为模型输入,以旋转机械的故障类型为输出,构建置信规则库(BRB)模型;然后规则库中的某条或多条置信规则将被输入的振动特征数据不同程度的激活,结合证据推理算法融合被激活规则,从融合的后项置信结构中换算出故障类型,其诊断结果不仅能判断当前故障类型,还可给出潜在故障发生概率,并且模型参数物理意义明确,可解释性强。
1、故障特征数据采集
本发明专利以旋转机械中基本的电机转子为例,采用ZHS-2型多功能电机柔性转子系统,来验证所提方法的有效性。在该试验台转子支撑部位的不同方向上,安装5个频域振动加速度传感器感知电机转子故障状态。由于不同电机转子故障类型会引起多频率成分的振动幅值变化,故获取1X~3X倍频的频域幅值变化为模型输入。设置电机转速为1500r/m,则基频1X为25Hz,n倍频nX(n=1,2,3,…)为(n×25)Hz。在每个电机转子故障类型下,获取3个倍频的15个故障特征(c1,c2,…,c15),采集210次(其中100次数据用于模型训练,110次数据用于模型测试)。
2、基于MIV的故障特征约简加权
本发明专利的仿真实验结果是在以下仿真条件下测得:设置网络的最大迭代次数为2000,期望目标误差最小值为1.0*e-5,增量δ=0.1。基于此训练稳定的网络,分别采用MIV算法求取c1-c15故障特征的平均影响值,并根据MIV值从大到小顺序排列。其排序情况和计算结果如图2所示。选取c3,c5,c13为影响度高的故障特征,MIV值分别为1.1143、0.6876、0.4902,并根据公式6,计算选取故障特征的相对权重为0.4057,0.2504,0.1785。
3、置信规则库模型的构建
通过电机柔性转子试验台采集故障特征c3,c5,c13的600(100次×6)组数据如图3所示,可以明显看出通过MIV值选取的故障特征更能准确表征出6种故障类型发生趋势。
设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
对于置信规则库的后项属性,采用独热编码对其进行描述。置信规则就可以被表示为Rk:如果且且那么Y是{(D1,βk,1),…,(DN,βk,N)}其中N=6,k=1,2,…,125,每条规则权重θk的初始值都设为1,属性权重初始值根据MIV值设置为δ1=0.4057,δ2=0.2504,δ3=0.1785。部分置信规则如表2所示。
2、置信规则库故障诊断的过程
在基于专家知识和历史数据构建好初始置信规则库后,具体故障诊断过程如下:
4.1计算输入量的匹配度
将故障特征c3,c5,c13的600组样本数据作为输入,利用公式13和14计算每一组故障特征数据对于各自参考点的匹配度。例如,当第118组特征数据被输入时,输入量为X=[0.0776、0.0855、0.1380],则X1匹配参考点小正(PS)和正中(PM)的程度为0.4478和0.5522,X2匹配参考点正中(PM)和正大(ML)的程度为0.8773和0.1227,X3匹配小正(PS)和正中(PM)的程度为0.5585和0.4415,对于其他参考点的匹配程度均为0。
4.2计算被激活规则的权重
获得输入量对每个规则中参考点的匹配度之后,利用公式11和12计算被激活规则的权重wk。例如对于第118组特征数据,可以得到其对R37-R38、R42-R43、R62-R63和R67-R68的激活权重分别为w37=0.2195,w38=0.1734,w42=0.0307,w43=0.0243,w62=0.2706,w63=0.2139,w67=0.0378,w68=0.0299而其他规则的激活权重均为0,也就是激活了8条规则。
4.3融合被激活的规则
利用证据推理算法得到融合后的输出置信结构
O(X)={(Fj,βj),j=1,2,…,6}
其中,Fj和βj可以分别通过公式16和17计算得出。例如,将第二步关于第118组特征数据的wk及βj,k带入公式15,可以计算得出
O(X)={(F1,0),(F2,0),(F3,0),(F4,0.0073),(F5,0.7968),(F6,0.1959)}。
4.4输出电机转子故障类型
根据融合结果中对F1-F6的置信度,置信度最大的后项属性对应的故障类型即为诊断结果。例如,对于第118组特征数据,其对应的故障诊断结果为连接器松动故障,潜在齿轮缺齿故障发生概率为19.59%,并利用公式18计算其基于独热编码的欧氏距离误差Δ=0.2032。
将600组c3、c5、c13数据样本输入其初始置信规则库中,得到的诊断结果如表4所示,其中本文所提方法对6种故障模式的确诊率分别为78%、88%、80%、91%、86%、80%,整体确诊率为83.83%。同时为分析该方法诊断的离散程度,随机选取50组样本数据为1大组,共计12大组,计算12大组数据的诊断结果与真实故障之间的均方根误差(RMSE)为0.4572,说明该方法诊断的离散程度较小,诊断信度更集中,通过MIV选取的故障特征更能表征故障发生类型。
表4初始置信规则库的诊断结果
5、置信规则库模型参数的优化。
由于初始置信规则库并不十分准确,需根据步骤五给出的优化算法,利用大量训练数据优化模型。具体步骤如下:
5.1:计算样本实际故障类型的独热编码与初始置信规则库输出之间的均方根误差
5.2:利用步骤五给出的优化模型,找到使ξ(P)为最小值时的最优参数集合P,约束条件为:0≤θk≤1,0≤δk≤1,0≤βj,k≤1和(k=1,2,…125)。本文利用600组训练数据,采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)<0.0016时的最优参数集,部分规则如表5所示。
表5.优化后置信规则库的部分规则
优化后的置信规则库对电机转子故障诊断的混淆矩阵如表6所示,6种故障类型的确诊率分别达到了91%、95%、93%、98%、97%、90%,总体故障准确率可达94%,要明显高于优化前的确诊率。进一步对比分析在12大组数据中,优化前后的置信规则库方法故障诊断结果的均方根误差,如图4所示。优化后的所提方法RMSE可达0.2457,对比说明优化后方法可以有效降低故障诊断的离散程度,提高诊断准确性。
表6优化后置信规则库的诊断结果
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.基于元素影响度和置信规则库的旋转机械故障诊断方法,首先进行设备安装,包括在电机柔性转子试验台的转自支撑部位的不同方向上,安装多个频域振动加速度传感器,其特征在于,该方法包括,
步骤一,建立筛选模型,通过所述频域震动加速度传感器得到故障特征因子,并对其进行加权分配比重;
步骤二,分类故障类型,得到根据步骤一中所得模型,设置参考点,建立置信规则库,采用独热编码对其进行描述;
步骤三,设置置信规则库,筛选出故障特征的样本数据,以及筛选故障特征前项参考点和后项参考点,基于专家知识和历史数据构建初始置信规则库;
步骤四,执行故障诊断,根据步骤三所得的置信规则库,得到诊断结果;
步骤五,模型优化,基于独热编码和欧式距离进行模型参数优化。
2.根据权利要求1所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,设定电机转子的故障类型包括齿轮缺齿、连接器松动、基座松动、转子不对中、转子不平衡和正常运行,所述设备安装的准备步骤包括安装频域震动加速度传感器的个数为M。
3.根据权利要求2所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,计算机获取1X-3X倍频的频域幅值作为模型输入,在每个电机转子的故障类型下获取3个倍频的3×M个故障特征(c1,c2,…,c3M),
所述步骤一包括,根据与每种故障类型的相关性,从3×M个故障特征(c1,c2,…,c3M)中挑选相应的故障因子,并对所挑选的故障银子进行加权分配比重,从而得到筛选模型。
4.根据权利要求3所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤一中,建立如下模型
Q(ci,δi)=q(c1,c2,…,cn,F),
其中,Q(ci,δi)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。
5.根据权利要求4所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤一包括,
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新,定义故障特征的损失函数如下:
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
其中L表示故障特征因子数,n表示样本组数;
神经网络的拟合输出如下所示:
IV=[IV1,IV2,…,IVn]T,
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
根据故障特征因子对故障的平均影响程度,得出其相对应的约简加权信息Q(ci,δi)。
7.根据权利要求6所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤四包括,
S4-1,计算输入量的匹配度,
S4-2,计算被激活规则的权重,
S4-3,融合被激活的规则并输出诊断故障结果。
8.根据权利要求7所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述S4-1包括,
所述S4-2包括,
所述S4-3包括,
O(X)={(Fj,βj),j=1,2,…,6},
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110847686.8A CN113569478A (zh) | 2021-07-27 | 2021-07-27 | 基于元素影响度和置信规则库的旋转机械故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110847686.8A CN113569478A (zh) | 2021-07-27 | 2021-07-27 | 基于元素影响度和置信规则库的旋转机械故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113569478A true CN113569478A (zh) | 2021-10-29 |
Family
ID=78167717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110847686.8A Pending CN113569478A (zh) | 2021-07-27 | 2021-07-27 | 基于元素影响度和置信规则库的旋转机械故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113569478A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114291675A (zh) * | 2021-11-18 | 2022-04-08 | 杭州电子科技大学 | 一种基于复杂网络和置信规则推理的电梯故障诊断方法 |
CN115222299A (zh) * | 2022-09-20 | 2022-10-21 | 南京智欧智能技术研究院有限公司 | 一种基于旋转库的动态分配存储方法及系统 |
CN117629637A (zh) * | 2024-01-24 | 2024-03-01 | 哈尔滨师范大学 | 一种航空发动机轴承故障诊断方法及诊断系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106022366A (zh) * | 2016-07-04 | 2016-10-12 | 杭州电子科技大学 | 一种基于近邻证据融合的旋转机械设备故障诊断方法 |
CN110057581A (zh) * | 2019-04-29 | 2019-07-26 | 杭州电子科技大学 | 基于区间型信度规则推理的旋转机械故障诊断方法 |
CN110132603A (zh) * | 2019-05-16 | 2019-08-16 | 杭州电子科技大学 | 基于并集置信规则库和蚁群算法的船舶柴油机故障定位方法 |
CN110363344A (zh) * | 2019-07-11 | 2019-10-22 | 安徽理工大学 | 基于miv-gp算法优化bp神经网络的概率积分参数预测方法 |
CN112016471A (zh) * | 2020-08-27 | 2020-12-01 | 杭州电子科技大学 | 一种不完备样本条件下的滚动轴承故障诊断方法 |
-
2021
- 2021-07-27 CN CN202110847686.8A patent/CN113569478A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106022366A (zh) * | 2016-07-04 | 2016-10-12 | 杭州电子科技大学 | 一种基于近邻证据融合的旋转机械设备故障诊断方法 |
CN110057581A (zh) * | 2019-04-29 | 2019-07-26 | 杭州电子科技大学 | 基于区间型信度规则推理的旋转机械故障诊断方法 |
CN110132603A (zh) * | 2019-05-16 | 2019-08-16 | 杭州电子科技大学 | 基于并集置信规则库和蚁群算法的船舶柴油机故障定位方法 |
CN110363344A (zh) * | 2019-07-11 | 2019-10-22 | 安徽理工大学 | 基于miv-gp算法优化bp神经网络的概率积分参数预测方法 |
CN112016471A (zh) * | 2020-08-27 | 2020-12-01 | 杭州电子科技大学 | 一种不完备样本条件下的滚动轴承故障诊断方法 |
Non-Patent Citations (5)
Title |
---|
K ZHANG ET AL: "Review of multiple fault diagnosis methods", 《CONTROL THEORY & APPLICATIONS》 * |
吕延卓 等: "基于置信规则库的机械故障诊断方法", 《计算机测量与控制》 * |
姬思雨: "面向一类信息物理系统的数据预处理与故障诊断方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李正辉 等: "基于置信规则库推理的电机转子故障诊断方法", 《科学技术创新》 * |
杨隆浩 等: "基于关联系数标准差融合的置信规则库规则约简方法", 《信息与控制》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114291675A (zh) * | 2021-11-18 | 2022-04-08 | 杭州电子科技大学 | 一种基于复杂网络和置信规则推理的电梯故障诊断方法 |
CN114291675B (zh) * | 2021-11-18 | 2024-05-03 | 杭州电子科技大学 | 一种基于复杂网络和置信规则推理的电梯故障诊断方法 |
CN115222299A (zh) * | 2022-09-20 | 2022-10-21 | 南京智欧智能技术研究院有限公司 | 一种基于旋转库的动态分配存储方法及系统 |
CN117629637A (zh) * | 2024-01-24 | 2024-03-01 | 哈尔滨师范大学 | 一种航空发动机轴承故障诊断方法及诊断系统 |
CN117629637B (zh) * | 2024-01-24 | 2024-04-30 | 哈尔滨师范大学 | 一种航空发动机轴承故障诊断方法及诊断系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113569478A (zh) | 基于元素影响度和置信规则库的旋转机械故障诊断方法 | |
CN110135079B (zh) | 一种海洋石油井控装备的宏观弹性评价方法及系统 | |
CN107725283B (zh) | 一种基于深度信念网络模型的风机故障检测方法 | |
CN110441065B (zh) | 基于lstm的燃气轮机在线检测方法与装置 | |
JP3993825B2 (ja) | 器械を取り付けられた機器およびプロセスのための推論信号生成装置 | |
CN103592141B (zh) | 用于测试复杂系统的可靠性的方法 | |
Li et al. | Multi-sensor data-driven remaining useful life prediction of semi-observable systems | |
CN105973594A (zh) | 一种基于连续深度置信网络的滚动轴承故障预测方法 | |
CN110362886B (zh) | 一种基于不确定性分析的城镇砌体住宅安全评估方法 | |
CN103678881B (zh) | 一种基于人工免疫与证据理论相结合的复合故障诊断方法 | |
CN104408322B (zh) | 综合多源故障概率似然信度的旋转机械设备故障诊断方法 | |
CN106090626B (zh) | 一种供水管网异常侦测方法 | |
Park et al. | Model updating method for damage detection of building structures under ambient excitation using modal participation ratio | |
KR20170053692A (ko) | 커널 회귀 모델의 앙상블을 위한 장치 및 방법 | |
CN106127300A (zh) | 一种旋转机械健康状态预测方法 | |
CN112149953B (zh) | 基于多模联动和多阶段协作的机电设备运行安全评估方法 | |
CN105975797A (zh) | 一种基于模糊数据处理的产品早期故障根原因识别方法 | |
CN110261771A (zh) | 一种基于传感器互补性分析的故障诊断方法 | |
CN114266289A (zh) | 一种复杂装备健康状态评估方法 | |
CN115577587A (zh) | 一种历史建筑健康状态监测方法及系统 | |
Samhouri et al. | An intelligent machine condition monitoring system using time-based analysis: neuro-fuzzy versus neural network | |
Jiang et al. | Paired ensemble and group knowledge measurement for health evaluation of wind turbine gearbox under compound fault scenarios | |
CN117009828B (zh) | 基于故障匹配算法的水电设备故障诊断方法 | |
CN112096693B (zh) | 水下生产液压控制系统在线诊断与校验方法、存储介质和控制终端 | |
CN102788955B (zh) | 基于Kalman滤波的ESN的涡轮发电机的分类子模型的剩余寿命预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211029 |
|
RJ01 | Rejection of invention patent application after publication |