CN109100142B - 一种基于图论的轴承半监督故障诊断方法 - Google Patents
一种基于图论的轴承半监督故障诊断方法 Download PDFInfo
- Publication number
- CN109100142B CN109100142B CN201810672129.5A CN201810672129A CN109100142B CN 109100142 B CN109100142 B CN 109100142B CN 201810672129 A CN201810672129 A CN 201810672129A CN 109100142 B CN109100142 B CN 109100142B
- Authority
- CN
- China
- Prior art keywords
- matrix
- label
- semi
- bearing
- degree
- 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.)
- Active
Links
Classifications
-
- 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/04—Bearings
- G01M13/045—Acoustic or vibration analysis
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种基于图论的轴承半监督故障诊断方法,该方法首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络;再计算复杂网络的结构参数,提取度分布的均值和标准差及网络复杂指数;最后利用基于图的半监督学***衡下情况下的轴承故障诊断,故障识别准确率高,具有显著的使用价值。
Description
技术领域
本发明涉及一种基于图论的轴承半监督故障诊断方法,属于机械零部件故障诊断技术领域。
背景技术
滚动轴承是使用最频繁的机械零部件之一。据统计,在所有机械故障中,超过40%的故障由轴承引起。在工业现场中,标签样本的收集是一件艰难的工作,尤其是稀缺的故障样本,但却存在着大量的无标签样本。针对这种只有少量标签数据且类别不平衡的情况,传统的故障诊断算法往往表现很差。此外,轴承一般运行在不确定的变工况下,这也会对故障诊断造成干扰。因此,通过利用无标签数据来提高故障诊断准确率有着十分重要的意义。
可视图算法是一种将时序数据转化成复杂网络的方法,该方法将样本中的每一个点当成图中的节点,利用物理可视的原理构建复杂网络。提取网络结构参数作为轴承故障特征,可视图特征对变工况具有天然的不敏感性。
半监督学***滑假设目标方程的最优解将少量标签传播到所有无标签数据上,实现无标签样本的状态辨识。
发明内容
本发明提出了一种基于图论的轴承半监督故障诊断方法,首先将信号通过可视图算法转换成复杂网络,提取复杂网络的结构特征作为轴承故障特征,作为半监督分类器的输入,然后利用基于图的半监督学习算法将标签传播到无标签数据上实现故障诊断。
本发明提出了一种基于图论的轴承半监督故障诊断方法,具体包括以下步骤:
(一)首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络;
(二)计算复杂网络的结构参数,提取度分布的均值和标准差及网络复杂指数。
(三)最后利用基于图的半监督学习处理无标签样本,实现轴承故障诊断。
优选地,步骤(一)采用可视图算法具体为:
针对一个具有N个数据点的样本,首先通过如下方法将振动加速度信号转化成可视图:
信号的任意两点(ta,xa),(tb,xb),对应可视图中的两个节点,设这两点之间任意一点为(tc,xc),两个节点相连接的条件满足
得到形容每个节点的连接情况邻接矩阵WD,矩阵WD为N*N维的对称矩阵,wij=1代表节点i与节点j相连;
对于任意一点,与其连接点的个数称为该点的度则度分布为DV=[d1,…,dn]。
优选地,步骤(二)复杂网络的结构特征参数具体计算方法如下:
(1)度分布均值:
(2)度分布标准差:
(3)图复杂指数GIC定义如下:
C=4c(1-c)
其中,
λmax表示可视图邻接矩阵的最大特征值。
优选地,步骤(三)实现轴承故障诊断的具体方法如下:
(1)图的构建
首先,通过核函数计算所有样本之间的相似度,计算得到点xi和xj之间邻接矩阵K:K∈Rn×n,Kij=k(xi,xj),
通过将邻接矩阵乘以二值矩阵B∈Bn×n和距离矩阵H∈Rn×n完成稀疏化:
稀疏化邻接矩阵后,通过高斯核加权算法计算重新计算两点之间的权重矩阵W,点xi和xj之间的权重计算公式如下:
其中,d(xi,xj)代表xi与xj之间的欧式距离,δ表示数据点的标准差;
(2)标签传播
采用平滑项和标签匹配项作为最终优化目标:
s.t.yij∈{0,1},
yij=1,for label(xi)=j,j=1,…,c.
其中,基于度的权重矩阵Λ=diag([λ1,…,λn]),Y={yij}∈Βn×c为数据的标签矩阵,ΛY为规范化标签变量;
其中,d为某样本的度,p是类别数据的先验分布且
A=PTLP+μ(PT-I)(P-I)=PTLP+μ(P-I)2
其中,P为传播矩阵,L为归一化的拉普拉斯矩阵;
通过贪婪梯度Max-Cut算法求解上式,最终得到无标签样本的标签,实现故障诊断。
优选地,所述贪婪梯度Max-Cut算法流程为:
(1)输入:图GA={X,A}标签样本Xl和标签Y;
(2)初始化:
通过标签样本Xl构建初始割集{Sj}:
Sj={xi|yij=1},j=1,2,…,c;
无标签样本集:
Xu=X\Xl;
(3)以下步骤循环,直至满足条件
对于z=0至|Xu|做如下计算:
计算连接性:
将无标签样本加上最优标签,并将其放入标签割集{Sj};
将xi放入标签样本Xl:Xl←Xl+xi;
从无标签样本中删除该样本Xu:Xu←Xu-xi;
(4)输出
Sj,j=1,2,...,c。
其中样本xi与标签样本子集S的连接性定义如下:
Sj={xi|yij=1},i=1,2,…,n;j=1,2,…,c。
本发明的优点与积极效果在于:
(1)可视图算法是一种将时序数列转换成复杂网络的算法,该方法将轴承振动信号映射到图中,图的结构特征对信号的水平和垂直方向的尺度变换(对应轴承的转速和负载)具有天然的不敏感性。
(2)利用改进规范化标签变量的能平衡类别标签影响的基于图的半监督轴承故障诊断算法处理类别不平衡无标签数据。为了解决数据类别不平衡带来的影响,利用双变量的目标函数,其中规范化标签变量用来平衡类别的影响。
(3)相比于其他故障诊断方法,半监督故障诊断算法只需要少量的标签样本即能获得高的故障诊断准确率。
(4)本发明方法属于数据驱动方法,无需建立模型,即可实现故障的检测和识别,降低了专业要求,增加了工程应用性;
附图说明
图1是本发明故障诊断方法的整体步骤流程图。
图2是可视图算法示意图。
图3是某故障样本振动信号图。
图4是对应故障样本度分布。
图5是轻微故障不同类别比率下的分类准确率曲线图。
图6是常规故障不同类别比率下的分类准确率曲线图。
图7是严重故障不同类别比率下的分类准确率曲线图。
具体实施方式
下面将结合附图和实例对本发明作进一步的详细说明。
本发明提出了一种可视图算法将信号转换成复杂网络,通过提取可视图特征作为半监督学习算法的输入,从而实现故障诊断的方法。如图1所示,具体步骤如下:
步骤一、将轴承振动信号通过可视图算法转换成复杂网络。
针对一个具有N个数据点的样本,首先通过如下方法将信号转化成可视图。如图2为转换示意图。
数学计算方式如下:
信号的任意两点(ta,xa),(tb,xb),对应可视图中的两个节点,两个节点相连接的条件如下:
设这两点之间任意一点为(tc,xc),满足
由此可得到形容每个节点的连接情况邻接矩阵WD,矩阵为N*N维的对称矩阵,wij=1代表节点i与节点j相连。
对于任意一点,与其连接点的个数称为该点的度则度分布为DV=[d1,…,dn]。
步骤二、计算复杂网络的结构特征参数
复杂网络的基本特征包括度分布均值及标准差和图复杂指数。计算方法如下:
(1)度分布均值:
(2)度分布标准差:
(3)图复杂指数GIC定义如下:
C=4c(1-c)
其中,
λmax表示可视图邻接矩阵的最大特征值。
步骤三、用基于图的半监督学习算法,实现故障识别与诊断
基于图的半监督学习算法原理包括图的构建和标签传播两部分如下:
数学符号定义:{(x1,y1),…,(xl,yl)}和{xl+1,…,xl+u}分别代表标签样本和无标签样本,Xl={x1,…,xl}和Xu={xl+1,…,xl+u}分别代表有标签和无标签输入数据。l和u表示标签样本和无标签样本的数量。Y={yij}∈Βn×c为数据的标签矩阵,其中yij=1代表样本xi的标签为j,j∈{1,2,…,c},每一个样本有且只有一个标签,即半监督学习的目标是识别无标签数据的标签{yl+1,…,yl+u},其中l<<n(l+u=n)。
(3)图的构建
首先,通过核函数计算所有样本之间的相似度,计算得到邻接矩阵K:K∈Rn×n,Kij=k(xi,xj)。然后,分别稀疏化邻接矩阵和重置样本之间的权重。
系数化:通过将邻接矩阵乘以一个二值矩阵B∈Bn×n和距离矩阵H∈Rn×n完成的。
二值矩阵是通过k最近邻算法计算得到的,数学表达如下:
最终得到二值矩阵B:
重置权重:
稀疏化邻接矩阵后,通过高斯核加权算法计算重新计算两点之间的权重矩阵W。点xi和xj之间的权重计算公式如下:
其中,d(xi,xj)代表xi与xj之间的欧式距离。
(4)标签传播
基于图的半监督学***滑假设的基础上,本发明采用平滑项和标签匹配项作为最终优化目标。
s.t.yij∈{0,1},
yij=1,for label(xi)=j,j=1,…,c.
上式中,第一项为平滑性的度量,第二项为标签匹配的度量。其中,L为归一化的拉普拉斯矩阵
L=D-1/2(D-W)D-1/2=I-D-1/2WD-1/2
为规范化标签变量
其中,基于度的权重矩阵Λ=diag([λ1,…,λn])
d为某样本的度,p是类别数据的先验分布且本文设pj=1/c。
求解步骤分为两步:
(1)F变量的优化:
固定变量Y,求导取零值点
其中,P=(L/μ+I)-1定义为传播矩阵。
(2)Y变量的优化:
将上一步骤求解出的最优值F*代入原目标函数,得
其中,A=PTLP+μ(PT-I)(P-I)=PTLP+μ(P-I)2
最终的目标函数可写成如下形式:
s.t.yij∈{0,1},
yij=1,for label(xi)=j,j=1,…,c.
通过贪婪梯度Max-Cut算法求解上式:
算法流程:
1、输入:图GA={X,A}标签样本Xl和标签Y;
2、初始化:
通过标签样本Xl构建初始割集{Sj}:Sj={xi|yij=1},j=1,2,…,c
无标签样本集Xu=X\Xl;
3、以下步骤循环,直至满足条件
对于z=0至|Xu|做如下计算:
计算连接性:xi∈Xu,j=1,…,c
计算
将无标签样本加上最优标签,并将其放入标签割集{Sj};
将xi放入标签样本Xl:Xl←Xl+xi;
从无标签样本中删除该样本Xu:Xu←Xu-xi;
4、输出:Sj,j=1,2,...,c。
其中样本xi与标签样本子集S的连接性定义如下:
Sj={xi|yij=1},i=1,2,…,n;j=1,2,…,c
最终得到无标签样本的标签,实现故障诊断。
本实例采取凯斯西储大学轴承数据中心提供的滚动轴承故障信号进行验证。分别使用正常、内环故障、外环故障和滚动体故障四种状态下的样本信号对本发明基于图论的轴承半监督故障诊断方法进行检测验证,具体步骤如下:
步骤一、对轴承振动信号进行以固定点数进行分割。
四种状态下的信号样本数如表1所示。
表1四种状态下样本数(故障直径:0.007英寸)
针对每一个样本,将其转换成节点数为N的复杂网络,图3和图4分别为某个轴承故障样本的时域图及其对应的度分布图。
步骤二、提取可视图特征
表1中对应的故障样本的可视图中xyz轴分别代表的分布的均值、标准差和图复杂指数,代表正常样本、内圈故障、外圈故障和滚动体故障样本在特征空间中的位置。
步骤三、利用基于图的半监督学习算法识别无标签数据的状态,实现故障识别与诊断
针对不同的类别不平衡程度,设样本比率为r,选取带标签每类故障样本3个,则正常样本个数为3*r。r的取值范围为1~20。图5、6、7所示为分别代表不同故障程度(故障点半径不同)下的分类准确率,图中横坐标为r,即图中曲线代表不同类别样本比例下的故障诊断准确率。可知,该算法相对经典的半监督学习算法(LGC和HFGF)具有更好的分类准确率。
为了更准确的说明该方法的有效性,统计了r=16情况下各方法针对每个类别的准确率。试验结果如表2所示。
表2试验对比结果
根据试验结果可以看出基于图的半监督轴承故障诊断方法优于LGC和GFHF等经典半监督学习方法。可见本发明的方法能够实现对轴承的故障检测与分类,故障准确率高,且正常样本与故障样本之间的准确率为100%,具有明显的实际应用价值。
Claims (2)
1.一种基于图论的轴承半监督故障诊断方法,其特征在于,具体包括以下步骤:
(一)首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络;
(二)计算复杂网络的结构参数,提取度分布的均值和标准差及网络复杂指数;
(三)最后利用基于图的半监督学习处理无标签样本,实现轴承故障诊断;
步骤(一)采用可视图算法具体为:
针对一个具有N个数据点的样本,首先通过如下方法将振动加速度信号转化成可视图:
信号的任意两点(ta,xa),(tb,xb),对应可视图中的两个节点,设这两点之间任意一点为(tc,xc),两个节点相连接的条件满足
得到形容每个节点的连接情况邻接矩阵WD,矩阵WD为N*N维的对称矩阵,wij=1代表节点i与节点j相连;
对于任意一点,与其连接点的个数称为该点的度则度分布为DV=[d1,…,dn];
步骤(二)复杂网络的结构特征参数具体计算方法如下:
(1)度分布均值:
(2)度分布标准差:
(3)图复杂指数GIC定义如下:
C=4c(1-c)
其中,
λmax表示可视图邻接矩阵的最大特征值;
步骤(三)实现轴承故障诊断的具体方法如下:
(1)图的构建
首先,通过核函数计算所有样本之间的相似度,计算得到点xi和xj之间邻接矩阵K:K∈Rn×n,Kij=k(xi,xj),
通过将邻接矩阵乘以二值矩阵B∈Bn×n和距离矩阵H∈Rn×n完成稀疏化:
稀疏化邻接矩阵后,通过高斯核加权算法计算重新计算两点之间的权重矩阵W,点xi和xj之间的权重计算公式如下:
其中,d(xi,xj)代表xi与xj之间的欧式距离,δ表示数据点的标准差;
(2)标签传播
采用平滑项和标签匹配项作为最终优化目标:
s.t.yij∈{0,1},
yij=1,for label(xi)=j,j=1,…,c.
其中,基于度的权重矩阵Λ=diag([λ1,…,λn]),Y={yij}∈Bn×c为数据的标签矩阵,ΛY为规范化标签变量;
其中,d为某样本的度,p是类别数据的先验分布且
A=PTLP+μ(PT-I)(P-I)=PTLP+μ(P-I)2
其中,P为传播矩阵,L为归一化的拉普拉斯矩阵;
通过贪婪梯度Max-Cut算法求解上式,最终得到无标签样本的标签,实现故障诊断。
2.如权利要求1所述的方法,其特征在于,
所述贪婪梯度Max-Cut算法流程为:
(1)输入:图GA={X,A}标签样本Xl和标签Y;
(2)初始化:
通过标签样本Xl构建初始割集{Sj}:
Sj={xi|yij=1},j=1,2,…,c;
无标签样本集:
Xu=X\Xl;
(3)以下步骤循环,直至满足条件
对于z=0至|Xu|做如下计算:
计算连接性:
将无标签样本加上最优标签,并将其放入标签割集{Sj};
将xi放入标签样本Xl:Xl←Xl+xi;
从无标签样本中删除该样本Xu:Xu←Xu-xi;
(4)输出
Sj,j=1,2,...,c;
其中样本xi与标签样本子集S的连接性定义如下:
Sj={xi|yij=1},i=1,2,…,n;j=1,2,…,c。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810672129.5A CN109100142B (zh) | 2018-06-26 | 2018-06-26 | 一种基于图论的轴承半监督故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810672129.5A CN109100142B (zh) | 2018-06-26 | 2018-06-26 | 一种基于图论的轴承半监督故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100142A CN109100142A (zh) | 2018-12-28 |
CN109100142B true CN109100142B (zh) | 2019-09-13 |
Family
ID=64845068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810672129.5A Active CN109100142B (zh) | 2018-06-26 | 2018-06-26 | 一种基于图论的轴承半监督故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100142B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110057588B (zh) * | 2019-05-09 | 2020-07-03 | 山东大学 | 基于奇异值与图论特征融合的轴承早期故障检测与诊断方法及*** |
CN110265002B (zh) * | 2019-06-04 | 2021-07-23 | 北京清微智能科技有限公司 | 语音识别方法、装置、计算机设备及计算机可读存储介质 |
CN110823574B (zh) * | 2019-09-30 | 2021-04-06 | 安徽富煌科技股份有限公司 | 一种基于半监督学习深度对抗网络的故障诊断方法 |
CN111721534B (zh) * | 2020-06-18 | 2021-09-24 | 山东大学 | 一种滚动轴承健康状态在线评估方法及*** |
CN113503960A (zh) * | 2021-07-21 | 2021-10-15 | 美核电气(济南)股份有限公司 | 基于可视图的机械***故障诊断方法及*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861481A (zh) * | 2017-10-31 | 2018-03-30 | 东北大学 | 基于工业大数据图半监督新类挖掘的故障诊断方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106510619B (zh) * | 2016-10-11 | 2019-04-05 | 天津大学 | 基于复杂网络的心电信号分析方法及在智能穿戴上的应用 |
CN106473736B (zh) * | 2016-10-11 | 2019-05-21 | 天津大学 | 基于复杂网络的脑电信号分析方法及应用 |
CN106503800B (zh) * | 2016-10-11 | 2018-11-30 | 天津大学 | 基于复杂网络的深度学习模型及在测量信号分析中的应用 |
CN106533742B (zh) * | 2016-10-31 | 2019-05-14 | 天津大学 | 基于时间序列模式表征的加权有向复杂网络建网方法 |
CN106596900A (zh) * | 2016-12-13 | 2017-04-26 | 贵州电网有限责任公司电力科学研究院 | 一种基于改进图的半监督分类的变压器故障诊断方法 |
CN107065828B (zh) * | 2017-03-24 | 2019-05-21 | 东北大学 | 基于知识数据的半监督核判别分析故障诊断方法 |
CN107976992B (zh) * | 2017-11-29 | 2020-01-21 | 东北大学 | 基于图半监督支持向量机的工业过程大数据故障监测方法 |
-
2018
- 2018-06-26 CN CN201810672129.5A patent/CN109100142B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861481A (zh) * | 2017-10-31 | 2018-03-30 | 东北大学 | 基于工业大数据图半监督新类挖掘的故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109100142A (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100142B (zh) | 一种基于图论的轴承半监督故障诊断方法 | |
WO2021135630A1 (zh) | 基于grcmse与流形学习的滚动轴承故障诊断方法 | |
CN110135459B (zh) | 一种基于双三元组深度度量学习网络的零样本分类方法 | |
Gómez et al. | Determining the accuracy in image supervised classification problems | |
CN114970605A (zh) | 一种多模态特征融合神经网络的制冷设备故障诊断方法 | |
CN105678343A (zh) | 基于自适应加权组稀疏表达的水电机组噪声异常诊断方法 | |
CN114004252A (zh) | 一种轴承故障诊断的方法、装置以及设备 | |
CN111189638B (zh) | 基于hmm和qpso优化算法的轴承故障程度辨识方法 | |
Liu et al. | Fault diagnosis of rod pumping wells based on support vector machine optimized by improved chicken swarm optimization | |
CN109298225B (zh) | 一种电压量测数据异常状态自动识别模型***及方法 | |
CN105606914A (zh) | 一种基于iwo-elm的航空功率变换器故障诊断方法 | |
CN114429152A (zh) | 基于动态指数对抗性自适应的滚动轴承故障诊断方法 | |
CN103278326A (zh) | 风力发电机组齿轮箱的故障诊断方法 | |
CN115187832A (zh) | 一种基于深度学习与格拉姆角场图像的能源***故障诊断方法 | |
CN111680875A (zh) | 基于概率基线模型的无人机状态风险模糊综合评价方法 | |
CN114091504A (zh) | 一种基于生成对抗网络的旋转机械小样本故障诊断方法 | |
CN107944761A (zh) | 基于人工智能消费维权指数企业投诉预警监测分析方法 | |
WO2019178930A1 (zh) | 一种机械设备故障诊断方法 | |
Zheng et al. | Research on feature extracted method for flutter test based on EMD and CNN | |
Zhao et al. | Parallel adversarial feature learning and enhancement of feature discriminability for fault diagnosis of a planetary gearbox under time-varying speed conditions | |
CN111985158A (zh) | 一种基于迁移学习的电液伺服阀故障诊断方法及*** | |
CN113610148B (zh) | 一种基于偏置加权AdaBoost的故障诊断方法 | |
CN113033683B (zh) | 一种基于静态与动态联合分析的工业***工况监测方法和*** | |
Wang et al. | A method of modal parameter identification for wind turbine blade based on binocular dynamic photogrammetry | |
Wang et al. | Rolling bearing diagnosis method based on improved standardized variable distance fusion hierarchical state space correlation entropy |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |