发明名称 基于数据驱动和随机直觉模糊策略的轴承故障诊断方法
摘要 本发明公开了基于数据驱动和随机直觉模糊策略的轴承故障诊断方法。该方法包括:1)以大数据思想构建模糊专家系统,实现数据驱动;2)待检轴承故障特征隶属度的建立;3)随机模糊集匹配与直觉模糊集决策融合。本发明以大量的实验数据为基础,进行频谱分析,通过幅值出现位置找特征频率。并将模糊证据融合转变为随机模糊集匹配和直觉模糊集决策。不确定度信息也是频谱信息的一部分能量,证据理论直接忽略掉有失客观性。直觉模糊决策充分了考虑了能量耗损。
申请公布号 CN104361242B 申请公布日期 2017.04.05
申请号 CN201410669226.0 申请日期 2014.11.20
申请人 青岛理工大学 发明人 孙显彬;谭继文;文妍
分类号 G06F19/00(2011.01)I;G01M13/04(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 段秋玲
主权项 基于数据驱动和随机直觉模糊策略的轴承故障诊断方法,包括以下步骤:一、以实验手段构建模糊专家系统利用实验手段模拟轴承各种典型故障运行模式,对于某一种故障在相同的时间间隔<img file="FDA0001216018630000016.GIF" wi="67" he="42" />内连续观测n次,作为一组,重复m组,代表m个专家;对于第k组数据,经频谱变换后,模糊认为频谱图中第i次出现振动幅值点处的频率为i倍故障特征频率,i=1,2,…,p;(1)计算k组数据在i倍特征频率处的平均值M<sub>i,k</sub>=(x<sub>i,k,1</sub>+ x<sub>i,k,2</sub>+ x<sub>i,k,3·····</sub>+ x<sub>i,k,n</sub>)/n(2)计算k组数据在i倍特征频率处的标准差:<img file="FDA0001216018630000011.GIF" wi="1814" he="231" />(3)利用M<sub>i,k</sub>,σ<sub>i,k</sub>构造高斯隶属度函数<img file="FDA0001216018630000012.GIF" wi="733" he="127" />(4)构造i倍频率处有m个专家组成模糊专家系统故障特征的隶属度函数<img file="FDA0001216018630000013.GIF" wi="1026" he="336" /><img file="FDA0001216018630000015.GIF" wi="59" he="59" />表示传感器获得的测量值;M<sub>i,a</sub>=min<sub>k=1,2,…,m</sub>(M<sub>i,k</sub>),M<sub>i,b</sub>=max<sub>k=1,2,…,m</sub>(M<sub>i,k</sub>);因此,取某一种故障处的i倍故障特征频率,i=1,2…,p,可以建立p个高斯隶属度函数;那么如果轴承有q种故障,则需要建立p*q个隶属度函数;二、待检轴承故障特征隶属度的建立在相同的时间间隔<img file="FDA0001216018630000017.GIF" wi="73" he="49" />内连续观测n次,作为一组待检样本,求出在i倍特征频率处的待检隶属度函数,则可以建立p个高斯隶属度函数;在i倍特征频率处的隶属度曲线<img file="FDA0001216018630000021.GIF" wi="902" he="127" />三、随机模糊集匹配与直觉模糊集融合将待检模式和模糊专家系统故障模式的隶属度函数进行匹配,就可以得到从采集信息中提取的特征对故障档案库中各个故障的支持程度,利用模糊随机集表示随机集似然测度,表示各个故障假设集全集的BPA;以第i次出现幅值处认为是i倍特征频率为例,待检模式的隶属度曲线与模糊专家系统中q种故障的隶属度曲线匹配相交,各自交点的纵坐标的最大值即为似然测度值;(5)得到BPA向量:μ<sub>i</sub>=(μ<sub>i1</sub>,μ<sub>i2</sub>,μ<sub>i3</sub>,….μ<sub>iq</sub>),i=1,2,…,p然后进行不确定度计算和向量元素归一化,Θ=1‑max(μ<sub>i1</sub>,μ<sub>i2</sub>,μ<sub>i3</sub>,…,μ<sub>iq</sub>);得到包含不确定度的归一化向量:v<sub>i</sub>=(v<sub>i1</sub>,v<sub>i2</sub>,v<sub>i3</sub>,….v<sub>iq</sub>,Θ<sub>i</sub>),其中v<sub>i1</sub>+v<sub>i2</sub>+v<sub>i3</sub>+…+v<sub>iq</sub>+Θ<sub>i</sub>=1此时得到单一传感器在i倍特征频率处的p个证据向量,传统方法是直接忽略不确定度Θ<sub>i</sub>,根据模糊证据理论融合;(6)直觉模糊集融合;充分考虑到不确定度也代表一定的信息,v<sub>i,j</sub>写成(m<sub>i,j</sub>,n<sub>i,j</sub>)其中i=1,2,…,p;j=1,2,…,q;其中m<sub>i,j</sub>=v<sub>i,j</sub>,代表支持j故障的隶属度,n<sub>i,j</sub>=(1‑m<sub>i,j</sub>‑Θ<sub>i</sub>)代表肯定不是j故障的隶属度;这样由随机模糊集的单一的似然测度值转变为二元数对来表示,随机模糊集转变为随机直觉模糊集,更接近客观现实;对于单一传感器的i个证据变为多元决策融合:<img file="FDA0001216018630000022.GIF" wi="1285" he="87" />其中,j=1,2,…,q代表q种故障;即得到:f(v)=f(&lt;m<sub>1</sub>,n<sub>1</sub>&gt;,&lt;m<sub>2</sub>,n<sub>2</sub>&gt;,&lt;m<sub>3</sub>,n<sub>3</sub>&gt;,……&lt;m<sub>q</sub>,n<sub>q</sub>&gt;)(7)多元决策判定规则判定故障类型应有数值最大的m<sub>i</sub>且应大于某一门限值θ;当m<sub>i</sub>=m<sub>j</sub>时, 比较n<sub>i</sub>和n<sub>j</sub>,首选n<sub>j</sub>中数值小的故障类型;不确定度应小于某一门限值,即1‑(m<sub>i</sub>+n<sub>i</sub>)&lt;ζ;数值最大的m<sub>i</sub>与数值第二大的m<sub>j</sub>之间应该有一定差,即m<sub>i</sub>‑m<sub>j</sub>&gt;ξ;基于以上规则,判定出是第j种故障;(8)N个传感器情况:如果有N个传感器,以上单一传感器的决策级融合则转变为N个传感器的特征级融合,N个传感器分布在不同的位置,应该有不同的权重系数;ω(ω<sub>1</sub>,ω<sub>2</sub>,….ω<sub>N</sub>),且<img file="FDA0001216018630000031.GIF" wi="299" he="79" />利用公式<img file="FDA0001216018630000032.GIF" wi="1118" he="78" />对N个传感器进行第二次融合;仍然利用以上融合判定规则进行判定属于第j种故障。
地址 266000 山东省青岛市市北区抚顺路11号