CN113092044B - 一种基于加权水平可见图的旋转机械故障诊断方法 - Google Patents

一种基于加权水平可见图的旋转机械故障诊断方法 Download PDF

Info

Publication number
CN113092044B
CN113092044B CN202110347411.8A CN202110347411A CN113092044B CN 113092044 B CN113092044 B CN 113092044B CN 202110347411 A CN202110347411 A CN 202110347411A CN 113092044 B CN113092044 B CN 113092044B
Authority
CN
China
Prior art keywords
graph
neural network
weighted
training
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.)
Active
Application number
CN202110347411.8A
Other languages
English (en)
Other versions
CN113092044A (zh
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN202110347411.8A priority Critical patent/CN113092044B/zh
Publication of CN113092044A publication Critical patent/CN113092044A/zh
Application granted granted Critical
Publication of CN113092044B publication Critical patent/CN113092044B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/022Vibration control arrangements, e.g. for generating random vibrations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • G01M7/025Measuring arrangements

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明公开了一种基于加权水平可见图的旋转机械故障诊断方法,本方法将传感器时序信号转换为水平可见图,可以在保持时间序列非线性动态特性的同时,引入图的拓扑结构特性,并利用采样序号的差值对边进行加权,以减弱时序上相距过远但仍有连接的影响。进一步地,使用图神经网络对转换得到的加权水平可见图进行学习,所使用的图神经网络聚合层使用不同线性变换函数对中心节点和邻接节点加权聚合,有效提高模型表达能力。本发明将传感器时序信号转换为加权水平可见图,结合图神经网络,可以挖掘更多有助于故障诊断的特征,有效提高故障诊断精度。

Description

一种基于加权水平可见图的旋转机械故障诊断方法
技术领域
本发明属于机械设备故障诊断技术,具体涉及一种基于加权水平可见图的旋转机械故障诊断方法。
背景技术
随着工业4.0的到来,传统工业正向自动化、数字化、智能化转型,科学高效的机械设备健康管理需求日益迫切。以卷积神经网络、自动编码器、生成对抗网络等为代表的深度学习模型在故障诊断领域得到了广泛应用,解决了传统机器学习模型如多层感知器、支持向量机等方法需要手工提取特征的问题,但是输入特征仍主要集中在时域、频域或时频域,侧重于数据数值本身,缺乏对数据中隐藏关系和结构的挖掘与探索。
近年来,常用的将时间序列转换为图的方法有递归网络、转移网络和水平可见图,基于图的时间序列分析在气象、经济、生物医学等领域都取得相关应用。同时,图神经网络作为新兴的人工智能技术,可以处理非欧式空间的图数据,实现数据结构关系的学习与推理,是学习图数据的强有力的工具。
因此,针对传感器采集的机械设备运行状态的时序信号,构建从原始时域到图结构数据的映射模型,利用图神经网络在学习与推理图数据上的优势,挖掘机械故障信号的潜在结构上的关系,提取更多有益于识别诊断的特征,将图神经网络理论应用于机械故障诊断领域,可以为机械系统的故障诊断提供新的方法和手段。
发明内容
为解决上述问题,本发明公开了一种基于加权水平可见图的旋转机械故障诊断方法,能够挖掘更多有助于故障诊断的特征,有效提高故障诊断精度。
为达到上述目的,本发明的技术方案如下:
一种基于加权水平可见图的旋转机械故障诊断方法,包括以下步骤:
(1)根据传感信号建立相应的加权水平可见图,
(2)构建图神经网络模型并对其进行训练,
(3)将待诊断传感信号输入模型进行分类识别。
其中,步骤(1)所述的根据传感信号建立相应的加权水平可见图包括:
步骤A1,将采集的传感信号划分为长度为N的样本序列,即每个样本中包含N个采样点,构成数据集X={x1,x2,…,x[M/N]},其中M为信号总长度,[·]为取整运算;
步骤A2,将等比例的健康与不同类型故障的数据随机划分为训练集,验证集和测试集以供模型的训练与测试;
步骤A3,使用Z-score标准化方法对原始传感信号样本x={x1,x2,…,xN}进行预处理得到标准化样本xznor
Figure BDA0003001217150000021
其中mean(·)和std(·)分别为求均值操作符和求标准差操作符;
步骤A4,使用加权水平可见图算法将预处理后的传感信号转换为有权图G,节点为采样点,节点初始特征为标准化的采样值;
步骤(2)和步骤(3)所述的构建图神经网络模型及其训练,识别待诊断样本包括:
步骤B1,构建图神经网络模型整体框架,包括聚合层,激活函数,归一化层,全局池化层,全连接分类模块;
步骤B2,设置图神经网络参数,包括批处理大小,初始学习率,学习率衰减速率与倍数,优化器,聚合层层数K,隐层特征维数,早停耐心轮次p;
步骤B3,将训练集图数据输入模型,进行训练与优化;
步骤B4,以验证集识别精度为指标监督模型训练状态,若连续p个训练轮次都不再提升,则提前停止模型训练,以验证集识别精度最高时的模型为最优图神经网络模型;
步骤B5,将测试集图数据输入步骤B4所述最优图神经网络模型进行测试。
优选的,所述步骤A4中,使用加权水平可见图算法将标准化后的样本
Figure BDA0003001217150000022
转换为有权图G的具体步骤顺序如下:
步骤A4-1,若
Figure BDA0003001217150000023
其中1≤i<n<j≤N,即
Figure BDA0003001217150000024
Figure BDA0003001217150000025
间的存在未被的打断的水平连接线,则两采样点间存在边;
步骤A4-2,计算
Figure BDA0003001217150000026
Figure BDA0003001217150000027
间连接边的权重wij,若连接边存在,则:
Figure BDA0003001217150000028
其中i,j为采样点所对应采样序号;否则wij=0;
步骤A4-3,若已完成所有采样点遍历,跳转至步骤A4-4;否则,继续遍历采样点,跳转至步骤A4-1;
步骤A4-4,完成时序信号至所述有权图G的转换,所得图为加权水平可见图。
优选的,所述步骤B1中,所述图神经网络模型的第k层聚合层为:
Figure BDA0003001217150000029
其中
Figure BDA0003001217150000031
为节点v的第k层的隐层特征,节点u属于节点v的邻接节点集合N(v),
Figure BDA0003001217150000032
Figure BDA0003001217150000033
为两个不同的线性变换函数,wvu是节点v与节点u边evu的权重。
优选的,所述步骤B1中,所述图神经网络的聚合层的非线性激活函数采用线性整流函数(Rectified Linear Unit,ReLU),归一化层采用批量归一化(BatchNormalization,BN)。
优选的,所述步骤B1中,所述全局池化层为求和池化,即图整体特征为图中所有节点隐层特征之和:
Figure BDA0003001217150000034
优选的,所述图神经网络损失函数采用负对数似然函数(Negative Log-likelihood Loss Function):
L(yG)=-log(yG)(5)
其中yG为诊断故障类型的标签。
本发明的有益效果为:
本发明提出一种新型加权水平可见图,使用采样序号间隔的倒数作为相连边的权重,削弱了在采样时序上过远的连接节点的干扰,同时,在节点聚合时,使用不同线性变换函数对中心节点和邻接节点加权聚合,提高模型表达能力;加权水平可见图与图神经网络相结合可以有效挖掘故障信号的潜在结构关系,为图分类提供更多有用特征,实现高精度的机械设备故障诊断。
附图说明
图1为本发明的方法流程示意图;
图2为示例采样序列;
图3为示例采样序列的加权水平可见图;
图4为凯斯西储大学轴承故障数据集试验平台;
图5为凯斯西储大学轴承故障数据集样本的标准化结果;
图6为凯斯西储大学轴承故障数据集样本的加权水平可见图;
图7为本方法与其他方法在凯斯西储大学轴承故障数据集上的诊断效果对比。
具体实施方式
下面结合附图和具体实施方式,进一步阐明本发明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
图1是本发明所述的一种基于加权水平可见图的旋转机械故障诊断方法的流程图,首先根据传感信号建立相对应的图信号,再构建图神经网络模型并对其进行训练,最后将待诊断传感信号输入模型进行分类识别。
图2为一自定义时间序列x={x1,x2,…,x6},对于未标准化时间序列,构建其加权水平可见图时,首先判断任意俩节点间是否存在连接:若xi,xj>xn,其中1≤i<n<j≤6,即xi和xj间的存在未被的打断的水平连接线,则两采样点间存在边;计算xi和xj间连接边的权重wij,若连接边存在,则:
Figure BDA0003001217150000041
其中i,j为采样点所对应采样序号;否则wij=0。图2采样序列对应的加权水平可见图如图3所示。
对于实际机械设备传感信号,本实施例使用凯斯西储大学轴承故障数据集作为实验数据集(数据集网站链接:https://csegroups.case.edu/bearingdatacenter/home),图4为其试验平台,包括一个2马力的电动机,一个扭矩传感器和一个功率测试计,使用安装在驱动端的振动传感器的信号进行识别,采样率为12kHz,电机转速为1750rpm。用电火花分别在内圈,外圈,滚动体三个部位加工单点损伤,各有三种直径:0.1778mm(0.007英寸)、0.3556mm(0.014英寸)、0.5334mm(0.021英寸),与正常工况一起,共有10类数据,每类100个样本,每个样本包含1024个采样点。建立振动传感信号相对应的图信号的步骤包括:
步骤A1,将采集的传感信号划分为长度为1024的样本序列,即每个样本中包含1024个采样点,构成数据集X={x1,x2,…,x[M/N]},其中M为信号总长度,[·]为取整运算;
步骤A2,将等比例的健康与不同类型故障的数据随机划分为训练集,验证集和测试集以供模型的训练与测试,训练集,验证集,测试集的比例为8:1:1;
步骤A3,使用Z-score标准化方法对原始传感信号样本x={x1,x2,…,x1024}进行预处理得到标准化样本xznor
Figure BDA0003001217150000042
其中mean(·)和std(·)分别为求均值操作符和求标准差操作符,某样本预处理后的时序图如图5所示;
步骤A4,使用加权水平可见图算法将预处理后的传感信号转换为有权图G,节点为采样点,节点初始特征为标准化的采样值,图G结构如图6所示;
其中,使用加权水平可见图算法将标准化后的样本
Figure BDA0003001217150000043
转换为有权图G的具体步骤顺序如下:
步骤A4-1,若
Figure BDA0003001217150000051
其中1≤i<n<j≤1024,即
Figure BDA0003001217150000052
Figure BDA0003001217150000053
间的存在未被的打断的水平连接线,则两采样点间存在边;
步骤A4-2,计算
Figure BDA0003001217150000054
Figure BDA0003001217150000055
间连接边的权重wij,若连接边存在,则:
Figure BDA0003001217150000056
其中i,j为采样点所对应采样序号;否则wij=0;
步骤A4-3,若已完成所有采样点遍历,跳转至步骤A4-4;否则,继续遍历采样点,跳转至步骤A4-1;
步骤A4-4,完成时序信号至所述有权图G的转换,所得图为加权水平可见图;
在时序数据集转为图结构数据后,构建图神经网络模型,并对其进行训练以识别待诊断样本,具体步骤包括:
步骤B1,构建图神经网络模型整体框架,包括4层聚合层,ReLU激活函数,批归一化层(Batch Normalization),全局池化层(Global Sum Pooling),全连接分类模块中激活函数选择Softmax函数;
其中,图神经网络模型的第k层聚合层表达式为:
Figure BDA0003001217150000057
其中
Figure BDA0003001217150000058
为节点v的第k层的隐层特征,节点u属于节点v的邻接节点集合N(v),
Figure BDA0003001217150000059
Figure BDA00030012171500000510
为两个不同的线性变换函数,wvu是节点v与节点u边evu的权重;
全局求和池化层为所有节点隐层特征之和,即图整体特征为:
Figure BDA00030012171500000511
损失函数采用负对数似然函数(negative log-likelihood loss function):
L(yG)=-log(yG)
其中yG为诊断故障类型的标签。
步骤B2,设置图神经网络参数,包括批处理大小为16,初始学习率0.01,每训练5份轮次学习率衰减1倍,优化器为Adam,聚合层层数K为4,隐层特征维数32,早停耐心轮次p为10;
步骤B3,将训练集图数据输入模型,进行训练与优化;
步骤B4,以验证集识别精度为指标监督模型训练状态,若连续10个训练轮次都不再提升,则提前停止模型训练,以验证集识别精度最高时的模型为最优图神经网络模型;
步骤B5,将测试集图数据输入步骤B4所述最优图神经网络模型进行测试。
图7对比了本方法与其他方法在凯斯西储大学轴承故障数据集上的诊断效果,取10次数据随机划分的平均精度,和标准差作为对比标准,节点嵌入维度、不同图神经网络聚合层数保持一致。相比其他方法,本方法取得最优诊断效果。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。

Claims (6)

1.一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:包括以下步骤:
(1)根据传感信号建立相应的加权水平可见图;
(2)构建图神经网络模型并对其进行训练;
(3)将待诊断传感信号输入模型进行分类识别;
其中,步骤(1)所述的根据传感信号建立相应的加权水平可见图包括:
步骤A1,将采集的传感信号划分为长度为N的样本序列,即每个样本中包含N个采样点,构成数据集X={x1,x2,…,x[M/N]},其中M为信号总长度,[·]为取整运算;
步骤A2,将等比例的健康与不同类型故障的数据随机划分为训练集,验证集和测试集以供模型的训练与测试;
步骤A3,使用Z-score标准化方法对原始传感信号样本x={x1,x2,…,xN}进行预处理得到标准化样本xznor
Figure FDA0003001217140000011
其中mean(·)和std(·)分别为求均值操作符和求标准差操作符;
步骤A4,使用加权水平可见图算法将预处理后的传感信号转换为有权图G,节点为采样点,节点初始特征为标准化的采样值;
步骤(2)和步骤(3)所述的构建图神经网络模型及其训练,识别待诊断样本包括:
步骤B1,构建图神经网络模型整体框架,包括聚合层,激活函数,归一化层,全局池化层,全连接分类模块;
步骤B2,设置图神经网络参数,包括批处理大小,初始学习率,学习率衰减速率与倍数,优化器,聚合层层数K,隐层特征维数,早停耐心轮次p;
步骤B3,将训练集图数据输入模型,进行训练与优化;
步骤B4,以验证集识别精度为指标监督模型训练状态,若连续p个训练轮次都不再提升,则提前停止模型训练,以验证集识别精度最高时的模型为最优图神经网络模型;
步骤B5,将测试集图数据输入步骤B4所述最优图神经网络模型进行测试。
2.根据权利要求1所述的一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:所述步骤A4中,使用加权水平可见图算法将标准化后的样本
Figure FDA0003001217140000012
转换为有权图G的具体步骤顺序如下:
步骤A4-1,若
Figure FDA0003001217140000013
其中1≤i<n<j≤N,即
Figure FDA0003001217140000014
Figure FDA0003001217140000015
间的存在未被的打断的水平连接线,则两采样点间存在边;
步骤A4-2,计算
Figure FDA0003001217140000021
Figure FDA0003001217140000022
间连接边的权重wij,若连接边存在,则:
Figure FDA0003001217140000023
其中i,j为采样点所对应采样序号;否则wij=0;
步骤A4-3,若已完成所有采样点遍历,跳转至步骤A4-4;否则,继续遍历采样点,跳转至步骤A4-1;
步骤A4-4,完成时序信号至所述有权图G的转换,所得图为加权水平可见图。
3.根据权利要求1所述的一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:所述步骤B1中,所述图神经网络模型的第k层聚合层为:
Figure FDA0003001217140000024
其中
Figure FDA0003001217140000025
为节点v的第k层的隐层特征,节点u属于节点v的邻接节点集合N(v),
Figure FDA0003001217140000026
Figure FDA0003001217140000027
为两个不同的线性变换函数,wvu是节点v与节点u边evu的权重。
4.根据权利要求1所述的一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:所述步骤B1中,所述图神经网络的聚合层的非线性激活函数采用线性整流函数,归一化层采用批量归一化。
5.根据权利要求1所述的一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:所述步骤B1中,所述全局池化层为求和池化,即图整体特征为图中所有节点隐层特征之和:
Figure FDA0003001217140000028
6.根据权利要求1所述的一种基于加权水平可见图的旋转机械故障诊断方法,其特征在于:所述图神经网络损失函数采用负对数似然函数:
L(yG)=-log(yG) (5)
其中yG为诊断故障类型的标签。
CN202110347411.8A 2021-03-31 2021-03-31 一种基于加权水平可见图的旋转机械故障诊断方法 Active CN113092044B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110347411.8A CN113092044B (zh) 2021-03-31 2021-03-31 一种基于加权水平可见图的旋转机械故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110347411.8A CN113092044B (zh) 2021-03-31 2021-03-31 一种基于加权水平可见图的旋转机械故障诊断方法

Publications (2)

Publication Number Publication Date
CN113092044A CN113092044A (zh) 2021-07-09
CN113092044B true CN113092044B (zh) 2022-03-18

Family

ID=76671682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110347411.8A Active CN113092044B (zh) 2021-03-31 2021-03-31 一种基于加权水平可见图的旋转机械故障诊断方法

Country Status (1)

Country Link
CN (1) CN113092044B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116226779A (zh) * 2021-11-29 2023-06-06 新智我来网络科技有限公司 一种联合学习参数聚合方法、装置及系统
CN117153188B (zh) * 2023-08-29 2024-05-28 安徽理工大学 一种基于声音频率的变功率设备故障分析方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102854015B (zh) * 2012-10-15 2014-10-29 哈尔滨理工大学 一种滚动轴承故障位置及性能退化程度诊断方法
EP3705953B1 (en) * 2019-03-06 2023-08-30 Robert Bosch GmbH Control of a physical system based on inferred state
CN112541233A (zh) * 2019-09-20 2021-03-23 宫文峰 基于改进卷积神经网络的旋转机械故障诊断方法
CN111966076B (zh) * 2020-08-11 2023-06-09 广东工业大学 基于有限状态机和图神经网络的故障定位方法
CN112269901B (zh) * 2020-09-14 2021-11-05 合肥中科类脑智能技术有限公司 一种基于知识图谱的故障判别推理方法

Also Published As

Publication number Publication date
CN113092044A (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN111914883B (zh) 一种基于深度融合网络的主轴轴承状态评估方法及装置
CN110162018B (zh) 基于知识蒸馏与隐含层共享的增量式设备故障诊断方法
CN106124212B (zh) 基于稀疏编码器和支持向量机的滚动轴承故障诊断方法
CN111046916A (zh) 一种基于空洞卷积胶囊网络的电机故障诊断方法及系统
CN113092044B (zh) 一种基于加权水平可见图的旋转机械故障诊断方法
CN109946080B (zh) 一种基于嵌入式循环网络的机械设备健康状态识别方法
CN114358123B (zh) 一种基于深度对抗迁移网络的广义开放集故障诊断方法
CN113505655A (zh) 面向数字孪生系统的轴承故障智能诊断方法
CN113723489A (zh) 一种改进关系网络的滚动轴承故障识别方法
CN110018322B (zh) 一种基于深度学习的转速检测方法及系统
CN111553297A (zh) 一种基于2d-cnn和dbn的涤纶长丝生产故障诊断方法和系统
CN116358871B (zh) 基于图卷积网络的滚动轴承微弱信号复合故障诊断方法
CN114548199A (zh) 一种基于深度迁移网络的多传感器数据融合方法
CN115290326A (zh) 一种滚动轴承故障智能诊断方法
CN115859077A (zh) 一种变工况下多特征融合的电机小样本故障诊断方法
CN113076920A (zh) 一种基于非对称域对抗自适应模型的智能故障诊断方法
CN112305388A (zh) 一种发电机定子绕组绝缘局部放电故障在线监测诊断方法
CN117056865B (zh) 一种基于特征融合的机泵设备运行故障诊断方法及装置
CN113758709A (zh) 结合边缘计算和深度学习的滚动轴承故障诊断方法及系统
CN116977708B (zh) 一种基于自适应聚合可视图的轴承智能诊断方法及系统
CN116821762A (zh) 一种基于多尺度图注意力融合网络的机械故障诊断方法
CN116738367A (zh) 基于多尺度图Transformer的滚动轴承故障诊断方法
CN116484258A (zh) 电梯曳引机轴承故障诊断方法
CN115876467A (zh) 一种基于伪标签传递式两阶段领域自适应的滚动轴承故障诊断方法
CN113807431A (zh) 一种基于多源信息融合的智能主轴状态评估方法和系统

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