CN113569478A - 基于元素影响度和置信规则库的旋转机械故障诊断方法 - Google Patents

基于元素影响度和置信规则库的旋转机械故障诊断方法 Download PDF

Info

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
Application number
CN202110847686.8A
Other languages
English (en)
Inventor
李正辉
张娜
杨媛媛
侯园园
李晓丽
王智伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhengzhou Railway Vocational and Technical College
Original Assignee
Zhengzhou Railway Vocational and Technical College
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhengzhou Railway Vocational and Technical College filed Critical Zhengzhou Railway Vocational and Technical College
Priority to CN202110847686.8A priority Critical patent/CN113569478A/zh
Publication of CN113569478A publication Critical patent/CN113569478A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/211Selection of the most significant subset of features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability 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(cii)=q(c1,c2,…,cn,F),
其中,Q(cii)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。
作为进一步优化,所述步骤一包括,
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新,定义故障特征的损失函数如下:
Figure BDA0003181288200000021
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
Figure BDA0003181288200000022
Figure BDA0003181288200000023
其中L表示故障特征因子数,n表示样本组数;
神经网络的拟合输出如下所示:
Figure BDA0003181288200000031
Figure BDA0003181288200000032
S1-3令
Figure BDA0003181288200000033
则有
Figure BDA0003181288200000034
综上,取0.1≤δ≤0.3,i=1,2,…,n。
Figure BDA0003181288200000035
Figure BDA0003181288200000036
分别表示样本集
Figure BDA0003181288200000037
Figure BDA0003181288200000038
的网络输出结果,样本中故障特征变量对故障类型的影响程度表示为
IV=[IV1,IV2,…,IVn]T,
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
Figure BDA0003181288200000039
根据故障特征因子对故障的平均影响程度,得出其相对应的约简加权信息Q(cii)。
作为进一步优化,所述步骤二包括,
根据故障特征对故障诊断的平均影响程度筛选出相应的故障特征cI,cII,cIII,并根据筛选出的故障特征平均影响度值对其进行加权,加权属性权重为δI,δIIIII,根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
Figure BDA00031812882000000310
Figure BDA00031812882000000311
Figure BDA00031812882000000312
对以上数据采用独热编码对其进行描述。
作为进一步优化,所述步骤四包括,
S4-1,计算输入量的匹配度,
S4-2,计算被激活规则的权重,
S4-3,融合被激活的规则并输出诊断故障结果。
作为进一步优化,所述S4-1包括,
将筛选故障特征的样本数据作为输入,输入变量xi与前项参考值
Figure BDA0003181288200000041
存在匹配度
Figure BDA0003181288200000042
Figure BDA0003181288200000043
如果xi小于等于
Figure BDA0003181288200000044
中的最小值,xi
Figure BDA0003181288200000045
的匹配程度
Figure BDA0003181288200000046
是0;如果大于等于
Figure BDA0003181288200000047
中的最大值时,则xi
Figure BDA0003181288200000048
的匹配程度
Figure BDA0003181288200000049
是1;否则xi
Figure BDA00031812882000000410
的匹配程度
Figure BDA00031812882000000411
则为
Figure BDA00031812882000000412
Figure BDA00031812882000000413
所述S4-2包括,
获得输入量对每个规则中参考点的匹配度
Figure BDA00031812882000000414
之后,置信规则库中的规则将被输入故障特征数据不同程度的激活,激活权重wk
Figure BDA00031812882000000415
其中
Figure BDA00031812882000000416
为被激活规则的相对属性权重
Figure BDA00031812882000000417
所述S4-3包括,
根据wk
Figure BDA00031812882000000418
对不同程度激活的置信规则进行折扣,进而利用证据推理理论融合折扣后置信规则的后项结构,融合结果为
O(X)={(Fjj),j=1,2,…,6},
其中
Figure BDA00031812882000000419
Figure BDA00031812882000000420
式中,
Figure BDA00031812882000000421
且0≤βj≤1,βj表示融合结果对故障模式Fj的支持程度,若MAX(β1,β2,…,βj,…β6)=βj,则表示旋转机械发生故障Fj的可能性最高,诊断结果为
Figure BDA00031812882000000422
作为进一步优化,所述步骤五包括,
S5-1,对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,将系统输出与实际输出的欧式距离为参数设计优化目标函数
Figure BDA0003181288200000051
S5-2,以
Figure BDA0003181288200000052
为约束条件得到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性权重δk)为最小值时的最优参数集合P。
作为进一步优化,所述步骤五还包括,采用随机梯度下降法来对模型进行训练,获得ξ(P)的最优参数集,从多功能电机柔性转子试验台的M个频域振动加速度传感器获取到的3M个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
Figure BDA0003181288200000053
本发明的工作原理及有益效果为:
通过对非线性振动数据建模,及时诊断旋转机械发生故障类型,对提高旋转机械维护的“状态修”水平和实现健康周期管理起到辅助决策的作用。其优势在于利用监测振动数据根据专家经验知识建立的模型对旋转机械故障进行诊断,在不进行停机维护的情况下判断旋转机械的健康状态,从而自适应调整修护周期,保障旋转机械系统稳定、安全运行。通过多源传感器对旋转机械运行情况进行监测,利用与故障类型具有高相关性的监测数据建模诊断故障类型,对指导旋转机械维修和保证旋转机械健康状态有着积极作用。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图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(cii)=q(c1,c2,…,cn,F), (1)
其中,Q(cii)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。具体步骤如下:
S1-1利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新。定义故障特征的损失函数如下:
Figure BDA0003181288200000061
S1-2对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
Figure BDA0003181288200000071
Figure BDA0003181288200000072
其中L表示故障特征因子数,n表示样本组数。神经网络的拟合输出如下所示:
Figure BDA0003181288200000073
Figure BDA0003181288200000074
S1-3令
Figure BDA0003181288200000075
则有
Figure BDA0003181288200000076
综上,取0.1≤δ≤0.3,i=1,2,…,n。
Figure BDA0003181288200000077
Figure BDA0003181288200000078
分别表示样本集
Figure BDA0003181288200000079
Figure BDA00031812882000000710
的网络输出结果。样本中故障特征变量对故障类型的影响程度表示为
IV=[IV1,IV2,…,IVn]T, (8)
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
Figure BDA00031812882000000711
根据故障特征因子对故障的平均影响程度,可以得出其相对应的约简加权信息Q(cii)。
步骤二:根据步骤一中利用平均影响度筛选出对故障诊断影响程度高的故障特征cI,cII,cIII,并根据筛选出的故障特征平均影响度值对其进行加权,加权属性权重为δI,δIIIII。根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
Figure BDA00031812882000000712
Figure BDA00031812882000000713
Figure BDA00031812882000000714
对于置信规则库的后项属性,采用独热编码对其进行描述,旋转机械故障编码见表1。
表1.故障类型及其编码
Figure BDA0003181288200000081
步骤三:根据步骤二中对筛选故障特征前项参考点和后项属性的定义,置信规则就可以被表示为Rk:如果
Figure BDA0003181288200000082
Figure BDA0003181288200000083
Figure BDA0003181288200000084
那么Y是{(D1k,1),…,(DNk,N)},
Figure BDA0003181288200000085
k=1,2,…,每条规则权重θk的初始值都设为1,属性权重初始值根据MIV值设置为δI,δII,δIII
部分置信规则如表2所示
表2.初始置信规则库的部分规则
Figure BDA0003181288200000086
步骤四:在步骤三中基于专家知识和历史数据构建好初始置信规则库后,具体证据推理故障诊断过程如下:
S4-1,计算输入量的匹配度。将筛选故障特征的样本数据作为输入,输入变量xi与前项参考值
Figure BDA0003181288200000087
存在一定的匹配度
Figure BDA0003181288200000088
Figure BDA0003181288200000089
如果xi小于等于
Figure BDA00031812882000000810
中的最小值,或者大于等于
Figure BDA00031812882000000811
中的最大值时,则xi
Figure BDA00031812882000000812
的匹配程度
Figure BDA00031812882000000813
是0或1;否则xi
Figure BDA00031812882000000814
的匹配程度
Figure BDA00031812882000000815
则为
Figure BDA0003181288200000091
Figure BDA0003181288200000092
S4-2,计算被激活规则的权重。获得输入量对每个规则中参考点的匹配度
Figure BDA0003181288200000093
之后,置信规则库中的规则将被输入故障特征数据不同程度的激活,激活权重wk
Figure BDA0003181288200000094
其中
Figure BDA0003181288200000095
为被激活规则的相对属性权重
Figure BDA0003181288200000096
S4-3,融合被激活的规则并输出诊断故障结果。根据wk
Figure BDA0003181288200000097
对不同程度激活的置信规则进行折扣,进而利用证据推理理论融合折扣后置信规则的后项结构,融合结果为
O(X)={(Fjj),j=1,2,…,6}, (14)
其中
Figure BDA0003181288200000098
Figure BDA0003181288200000099
式中,
Figure BDA00031812882000000910
且0≤βj≤1,βj表示融合结果对故障模式Fj的支持程度,若MAX(β12,…,βj,…β6)=βj,则表示旋转机械发生故障Fj的可能性最高,诊断结果为
Figure BDA00031812882000000911
步骤五:基于独热编码和欧式距离的模型参数优化,具体步骤如下:
S5-1,由于旋转机械设备内部机理复杂,专家难以确定故障诊断模型中所有参数的精确值,为此需要将初始参数微调,提高模型的精度。但置信规则库推理的融合结果为发生各类故障的可信度,并不能直接度量模型诊断效果和优化模型参数,故对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,旋转机械故障编码见表1。
将系统输出与实际输出的欧式距离为参数设计优化目标函数
Figure BDA00031812882000000912
S5-2,利用S5-1给出的优化目标函数,找到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性0≤θk≤1,0≤δk≤1,0≤βj,k≤1和
Figure BDA0003181288200000101
权重δk)为最小值时的最优参数集合P,约束条件为=1。采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)的最优参数集,从多功能电机柔性转子试验台的M个频域振动加速度传感器获取到的3M个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
Figure BDA0003181288200000103
具体实施例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(cii)=q(c1,c2,…,cn,F), (1)
其中,Q(cii)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。具体步骤如下:
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新。定义15个故障特征的损失函数如下:
Figure BDA0003181288200000102
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
Figure BDA0003181288200000111
Figure BDA0003181288200000112
其中L=15表示故障特征因子数,n表示样本组数。神经网络的拟合输出如下所示:
Figure BDA0003181288200000113
Figure BDA0003181288200000114
S1-3,根据15个故障特征因子对故障的平均影响程度,可以得出其相对应的约简加权信息Q(cii)。为了方便理解,设置网络的最大迭代次数为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组,选取参考点如下:
Figure BDA0003181288200000115
Figure BDA0003181288200000116
Figure BDA0003181288200000117
对于置信规则库的后项属性,采用独热编码对其进行描述,旋转机械故障编码见表1。
步骤三,根据步骤二中对筛选故障特征前项参考点和后项属性的定义,置信规则就可以被表示为Rk:如果
Figure BDA0003181288200000119
Figure BDA00031812882000001110
Figure BDA00031812882000001111
那么Y是{(D1k,1),…,(DNk,N)},
Figure BDA0003181288200000118
k=1,2,…,每条规则权重θk的初始值都设为1,属性权重初始值根据MIV值设置为δI,δII,δIII
为了便于理解,这里沿用步骤一和步骤二中筛选出的高影响度故障特征c3,c5,c13样本数据600组,以及筛选故障特征前项参考点和后项属性,可得到部分置信规则如表3所示:
表3.实施例2中初始置信规则库的部分规则
Figure BDA0003181288200000121
步骤四,在步骤三中基于专家知识和历史数据构建好初始置信规则库后,具体证据推理故障诊断过程如下:
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组样本数据作为训练数据,可以得到优化目标函数
Figure BDA0003181288200000131
S5-2,利用S5-1给出的优化目标函数,找到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性权重δk)为最小值时的最优参数集合P,约束条件为0≤θk≤1,0≤δk≤1,0≤βj,k≤1和
Figure BDA0003181288200000132
(k=1,2,…125)。采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)<0.0016时的最优参数集,从多功能电机柔性转子试验台的5个频域振动加速度传感器获取到的15个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
Figure BDA0003181288200000133
具体实施例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),其值如下:
Figure BDA0003181288200000141
Figure BDA0003181288200000142
Figure BDA0003181288200000143
对于置信规则库的后项属性,采用独热编码对其进行描述。置信规则就可以被表示为Rk:如果
Figure BDA0003181288200000144
Figure BDA0003181288200000145
Figure BDA0003181288200000146
那么Y是{(D1k,1),…,(DNk,N)}其中N=6,
Figure BDA0003181288200000147
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计算每一组故障特征数据对于各自参考点
Figure BDA0003181288200000148
的匹配度。例如,当第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计算被激活规则的权重
获得输入量对每个规则中参考点的匹配度
Figure BDA0003181288200000149
之后,利用公式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)={(Fjj),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初始置信规则库的诊断结果
Figure BDA0003181288200000151
5、置信规则库模型参数的优化。
由于初始置信规则库并不十分准确,需根据步骤五给出的优化算法,利用大量训练数据优化模型。具体步骤如下:
5.1:计算样本实际故障类型的独热编码与初始置信规则库输出之间的均方根误差
Figure BDA0003181288200000152
5.2:利用步骤五给出的优化模型,找到使ξ(P)为最小值时的最优参数集合P,约束条件为:0≤θk≤1,0≤δk≤1,0≤βj,k≤1和
Figure BDA0003181288200000161
(k=1,2,…125)。本文利用600组训练数据,采用随机梯度下降法来对模型进行训练,则可以获得ξ(P)<0.0016时的最优参数集,部分规则如表5所示。
表5.优化后置信规则库的部分规则
Figure BDA0003181288200000162
优化后的置信规则库对电机转子故障诊断的混淆矩阵如表6所示,6种故障类型的确诊率分别达到了91%、95%、93%、98%、97%、90%,总体故障准确率可达94%,要明显高于优化前的确诊率。进一步对比分析在12大组数据中,优化前后的置信规则库方法故障诊断结果的均方根误差,如图4所示。优化后的所提方法RMSE可达0.2457,对比说明优化后方法可以有效降低故障诊断的离散程度,提高诊断准确性。
表6优化后置信规则库的诊断结果
Figure BDA0003181288200000171
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.基于元素影响度和置信规则库的旋转机械故障诊断方法,首先进行设备安装,包括在电机柔性转子试验台的转自支撑部位的不同方向上,安装多个频域振动加速度传感器,其特征在于,该方法包括,
步骤一,建立筛选模型,通过所述频域震动加速度传感器得到故障特征因子,并对其进行加权分配比重;
步骤二,分类故障类型,得到根据步骤一中所得模型,设置参考点,建立置信规则库,采用独热编码对其进行描述;
步骤三,设置置信规则库,筛选出故障特征的样本数据,以及筛选故障特征前项参考点和后项参考点,基于专家知识和历史数据构建初始置信规则库;
步骤四,执行故障诊断,根据步骤三所得的置信规则库,得到诊断结果;
步骤五,模型优化,基于独热编码和欧式距离进行模型参数优化。
2.根据权利要求1所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,设定电机转子的故障类型包括齿轮缺齿、连接器松动、基座松动、转子不对中、转子不平衡和正常运行,所述设备安装的准备步骤包括安装频域震动加速度传感器的个数为M。
3.根据权利要求2所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,计算机获取1X-3X倍频的频域幅值作为模型输入,在每个电机转子的故障类型下获取3个倍频的3×M个故障特征(c1,c2,…,c3M),
所述步骤一包括,根据与每种故障类型的相关性,从3×M个故障特征(c1,c2,…,c3M)中挑选相应的故障因子,并对所挑选的故障银子进行加权分配比重,从而得到筛选模型。
4.根据权利要求3所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤一中,建立如下模型
Q(cii)=q(c1,c2,…,cn,F),
其中,Q(cii)为约简加权信息,δ为属性权重,q(c1,c2,…,cn,F)为输入信息,c表示被约简加权特征信息,F为故障类型。
5.根据权利要求4所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤一包括,
S1-1,利用BP神经网络的输入层到隐含层之间的连接权重出现一个很小变化Δωij,这个变化将会传递到隐含层的输出Sj,并使其发生改变Δsj,进而产生网络输出的变化ek,通过反向传递将权重ωij和ωjk进行更新,定义故障特征的损失函数如下:
Figure FDA0003181288190000021
S1-2,对于故障特征数据集X=[x(1),x(2),…,x(L)],分别对各个样本数据中的某个特征变量进行自增自减操作,如下所示
Figure FDA0003181288190000022
Figure FDA0003181288190000023
其中L表示故障特征因子数,n表示样本组数;
神经网络的拟合输出如下所示:
Figure FDA0003181288190000024
Figure FDA0003181288190000025
S1-3令
Figure FDA0003181288190000026
则有
Figure FDA0003181288190000027
综上,取0.1≤δ≤0.3,i=1,2,…,n。
Figure FDA0003181288190000028
Figure FDA0003181288190000029
分别表示样本集
Figure FDA00031812881900000210
Figure FDA00031812881900000211
的网络输出结果,样本中故障特征变量对故障类型的影响程度表示为
IV=[IV1,IV2,…,IVn]T,
进而对按照观测列数对IV取平均值可计算得出故障特征对最终输出故障类型的平均影响程度,表示为:
Figure FDA0003181288190000031
根据故障特征因子对故障的平均影响程度,得出其相对应的约简加权信息Q(cii)。
6.根据权利要求5所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤二包括,
根据故障特征对故障诊断的平均影响程度筛选出相应的故障特征cI,cII,cIII,并根据筛选出的故障特征平均影响度值对其进行加权,加权属性权重为δI,δIIII,根据筛选加权的高影响度故障特征数据,设置5个前项属性的参考点,即正极小(VS),小正(PS),正中(PM),正大(ML)和极大(VL),其值如下:
Figure FDA0003181288190000032
Figure FDA0003181288190000033
Figure FDA0003181288190000034
对以上数据采用独热编码对其进行描述。
7.根据权利要求6所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤四包括,
S4-1,计算输入量的匹配度,
S4-2,计算被激活规则的权重,
S4-3,融合被激活的规则并输出诊断故障结果。
8.根据权利要求7所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述S4-1包括,
将筛选故障特征的样本数据作为输入,输入变量xi与前项参考值
Figure FDA0003181288190000035
存在匹配度
Figure FDA0003181288190000036
Figure FDA0003181288190000037
如果xi小于等于
Figure FDA0003181288190000038
中的最小值,xi
Figure FDA0003181288190000039
的匹配程度
Figure FDA00031812881900000310
是0;如果大于等于
Figure FDA00031812881900000311
中的最大值时,则xi
Figure FDA00031812881900000312
的匹配程度
Figure FDA00031812881900000313
是1;否则xi
Figure FDA00031812881900000314
的匹配程度
Figure FDA00031812881900000315
则为
Figure FDA0003181288190000041
Figure FDA0003181288190000042
所述S4-2包括,
获得输入量对每个规则中参考点的匹配度
Figure FDA0003181288190000043
之后,置信规则库中的规则将被输入故障特征数据不同程度的激活,激活权重wk
Figure FDA0003181288190000044
其中
Figure FDA0003181288190000045
为被激活规则的相对属性权重
Figure FDA0003181288190000046
所述S4-3包括,
根据wk
Figure FDA0003181288190000047
对不同程度激活的置信规则进行折扣,进而利用证据推理理论融合折扣后置信规则的后项结构,融合结果为
O(X)={(Fjj),j=1,2,…,6},
其中
Figure FDA0003181288190000048
Figure FDA0003181288190000049
式中,
Figure FDA00031812881900000410
且0≤βj≤1,βj表示融合结果对故障模式Fj的支持程度,若MAX(β1,β2,…,βj,…,β6)=βj,则表示旋转机械发生故障Fj的可能性最高,诊断结果为
Figure FDA00031812881900000411
9.根据权利要求8所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤五包括,
S5-1,对各类故障类型进行独热编码,将原本离散的属性特征扩展至欧式空间,将系统输出与实际输出的欧式距离为参数设计优化目标函数
Figure FDA00031812881900000412
S5-2,以
0≤θk≤1,0≤δk≤1,0≤βj,k≤1和
Figure FDA0003181288190000051
为约束条件得到使ξ(P)(规则权重θk、后项输出置信度βj,k、属性权重δk)为最小值时的最优参数集合P。
10.根据权利要求9所述的基于元素影响度和置信规则库的旋转机械故障诊断方法,其特征在于,所述步骤五还包括,采用随机梯度下降法来对模型进行训练,获得ξ(P)的最优参数集,从多功能电机柔性转子试验台的M个频域振动加速度传感器获取到的3M个故障特征及其数据,经过步骤一的平均影响度筛选加权处理后,重复步骤二到步骤四即可获得更为准确的旋转机械故障的诊断结果
Figure FDA0003181288190000052
CN202110847686.8A 2021-07-27 2021-07-27 基于元素影响度和置信规则库的旋转机械故障诊断方法 Pending CN113569478A (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 杭州电子科技大学 一种不完备样本条件下的滚动轴承故障诊断方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
K ZHANG ET AL: "Review of multiple fault diagnosis methods", 《CONTROL THEORY & APPLICATIONS》 *
吕延卓 等: "基于置信规则库的机械故障诊断方法", 《计算机测量与控制》 *
姬思雨: "面向一类信息物理系统的数据预处理与故障诊断方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
李正辉 等: "基于置信规则库推理的电机转子故障诊断方法", 《科学技术创新》 *
杨隆浩 等: "基于关联系数标准差融合的置信规则库规则约简方法", 《信息与控制》 *

Cited By (5)

* Cited by examiner, † Cited by third party
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