CN114136565A - 一种多振源系统的故障诊断系统及方法 - Google Patents

一种多振源系统的故障诊断系统及方法 Download PDF

Info

Publication number
CN114136565A
CN114136565A CN202111230641.2A CN202111230641A CN114136565A CN 114136565 A CN114136565 A CN 114136565A CN 202111230641 A CN202111230641 A CN 202111230641A CN 114136565 A CN114136565 A CN 114136565A
Authority
CN
China
Prior art keywords
vibration
signal
composite
layer
fault
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.)
Granted
Application number
CN202111230641.2A
Other languages
English (en)
Other versions
CN114136565B (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN202111230641.2A priority Critical patent/CN114136565B/zh
Publication of CN114136565A publication Critical patent/CN114136565A/zh
Application granted granted Critical
Publication of CN114136565B publication Critical patent/CN114136565B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures
    • G01M7/02—Vibration-testing by means of a shake table
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00—Testing of vehicles
    • G01M17/08—Railway vehicles
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种多振源系统的故障诊断系统及方法，涉及机械故障诊断技术领域。本发明针对多振源系统中的复合振动信号，使用盲源分离方法实现复合振动信号的分离，获取独立振源的振动信号，生成正常状态下分离信号数据集。使用该分离信号数据集训练对抗式生成网络，通过判别网络对的独立振源的振动信号进行辨识，判断该独立振源是否处于异常状态，从而实现无样本下故障源的定位和达到对故障的定位和异常预警功能。该方法适用性强，能够在缺少故障样本的情况下，充分利用不同独立振源的正常稳定状态工作振动信号训练模型，并在故障发生时进行及时的定位和异常预警。

Description

一种多振源系统的故障诊断系统及方法
技术领域
本发明涉及一种电力机车转向架构架装配设备，特别是一种多振源系统的故障诊断系统及方法。
背景技术
高速列车是我国提升运输能力和带从产业发展的重要工具，其转向架旋转部件的振动信号监测和故障诊断对保障安全性与可持续性具有重要工程价值。转向架结构复杂、振源数多、空间密集，使其传感器布置数量有限，传统针对单部件的监控策略难以得到全面的监测信息。同时由于庞大的列车保有量和保守的维修策略，监测数据具有正常运行数据量巨大、故障运行数据稀少的特点，使基于分类方法的人工智能的诊断模型难以适用。传感器数量少、故障数据少成为了目前高速列车旋转部件状态监测和故障诊断的难点。因此，对有限数目传感器信号进行分离和识别，充分利用运行数据，对转向架旋转部件进行状态监测和故障预警，可以有效提升高速列车的安全性和可持续性，具有重要的工程应用价值。
现有的转向架旋转部件的振动监测系统，其传感器布置方式主要针对单个部件，使用变换域方法对信号进行处理，通过设置频谱或能量谱的阈值进行识别和预警。而高速列车转向架结构紧凑，传感器接收到的往往是多个振源的复合信号，其信号成分复杂。因此导致现有系统诊断精度低、易误诊，只能对单一部件严重故障进行有效的预警。
盲源分离技术可在信号源及其复合方式不确定的条件下，利用信源统计独立性，从复合信号中自适应分离信源信息，有效提升检测目标信噪比。同时，基于深度学习的人工智能模型在特征识别和分类任务上展现了卓越的效果，因此在故障诊断领域有着巨大的发展潜力。然而大部分模型在训练时需要大量的、模式类别丰富的样本数据。如卷积神经网络(CNN)，循环神经网络(RNN)、自动编码器(AE)等模型在有充足样本进行训练后，达到了较高的诊断精度。但在转向架旋转部件的诊断问题中大量数据模式单一，且故障样本稀少、模式不全，给模型的训练带来了巨大的挑战。目前针对少样本的故障诊断问题，提出了生成样本补充、模型精简、元学习等解决措施，但仍旧存在泛化能力弱、诊断精度不高等问题。
发明内容
本发明所要解决的技术问题是，针对现有技术不足，提供一种多振源系统的故障诊断系统及方法，从多振源系统中分离出独立振源的振动信号，并判断该振动信是否为故障信号。
为解决上述技术问题，本发明所采用的技术方案是：一种多振源系统的故障诊断方法，所述多振源系统包括多个独立振源，其特征在于，所述多振源系统还包括n个数据采集装置，所述数据采集装置的采样信号长度为m，包括以下步骤：
A，获取分离信号数据集；搭建生成对抗神经网络模型；
分离信号数据集的具体实现方式包括：
S1，获取多振源系统的复合振动信号；
S2，基于所述复合振动信号获取复合信号矩阵，对所述复合信号矩阵依次进行零均值化与白化处理，获取白化处理后的复合信号矩阵；
S3，基于白化处理后的复合信号矩阵获取分离信号矩阵；
S4，将所述分离信号矩阵切分为
Figure BDA0003315779960000021
个一维数组，输出单通道的波形.WAV文件，对于每一个独立振源，制作分离信号数据集；l为含振动特征的样本长度，根据所述独立振源的转速和所述信号采集装置的采样率确定；
所述生成对抗神经网络模型包括生成网络G和判别网络D；
所述生成网络G的网络结构包括：
第一层为全连接Dense层；
第二层为Reshape层，且其后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第三层至第六层均为1-D转置卷积层，卷积核大小为25，步长为4，且每层后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第七层为为1-D转置卷积层，卷积核大小为25，步长为4，且接有一个批标准化单元和LeakyRelu激活单元；
所述判别网络D的网络结构包括：
第一层为至第五层均为1-D卷积层，卷积核大小为25，步长为4，且每层后接一个LeakyReLu激活单元和一个PhaseShuffle相位重置单元；
第六层为Reshape层；
第七层为全连接Dense层；
B，对于每一个独立振源，基于所述分离信号数据集和所述生成对抗神经网络模型进行迭代训练；
C，对于每一个独立振源，获取振动状态判别器。
具体地，所述步骤B的具体实现方式包括：
S601，获取损失函数：
Figure BDA0003315779960000031
Figure BDA0003315779960000032
为真实样本和生成样本的线性插值，其计算公式为：
Figure BDA0003315779960000033
xr为真实样本信号；k为线性插值的权重，取值范围为k∈[0,1]；xf为生成样本信号；pr为真实样本的概率分布；pg为生成样本的概率分布；D(xr)为以xr为输入且满足一阶Lipschitz约束的判别网络函数，D(xf)为以xf为输入且满足一阶Lipschitz约束的判别网络函数；
Figure BDA0003315779960000034
为以
Figure BDA0003315779960000035
为输入且满足一阶Lipschitz约束的判别网络函数；
Figure BDA0003315779960000036
为判别网络函数D(xr)的期望；
Figure BDA0003315779960000037
为判别网络函数D(xf)的期望；
Figure BDA0003315779960000038
为函数
Figure BDA0003315779960000039
的期望。
S602，解码样本WAV文件为张量数据，裁剪长度并归一化幅值，打乱顺序后进行采样打包；
S603，对步骤S602获取的数据进行迭代训练。
具体地，获取所述损失函数的方法包括：
S6011，获取Wasserstein距离项W(pr,pg)，计算方法为：
Figure BDA00033157799600000310
θ为生成网络参数；
S6012，获取惩罚项GP的计算公式为：
Figure BDA00033157799600000311
Wasserstein又称为推土机距离(Earth-Move Distance,EM距离)，其定义为：
Figure BDA00033157799600000312
其中，∏(p,q)为分布p，q组合起来的所有可能的联合分布的集合，γ为服从∏(p,q)分布的γ样本，(x,y)服从联合分布γ，
Figure BDA00033157799600000313
为以(x,y)为输入的‖x-y‖的期望值，inf{}表示集合的下确界。
由于在计算生成样本分布pg与真实数据分布pr的是不现实的，因此将W(pr,pg)的计算方法转换为：
Figure BDA0003315779960000041
梯度惩罚项是为了迫使判别网络函数满足一阶Lipschitz约束.
具体地，步骤S603中，迭代训练的具体实现方式包括：
S6031，设定收敛条件为最大迭代次数N，随机化初始化生成器网络参数θ，判别网络参数φ；
S6032，进行次数为m的判别训网络训练：为从所述分离信号数据集中采集一批真实样本xr，输入判别网络输出D(xr；φ)，再从N(0,1)分布中采样一批随机向量h，输入生成器输出xf＝G(h；θ)，随后计算损失函数，使用Adam优化器使，设置超参数α,β1,β2，进行参数φ优化并更新，使损失函数值最大化：
Figure BDA0003315779960000042
S6033，进行生成器训练：采样N(0,1)分布输入生成器，计算损失函数，使用Adam优化器，设置超参数α1,β1,β2，进行参数θ的优化并更新，使损失函数
Figure BDA0003315779960000047
的值最小化：
Figure BDA0003315779960000043
S6034，训练迭代直至达到最大迭代次数N。
具体地，基于步骤S1至步骤S3，对于每一个独立振源，周期性获取分离信号矩阵，输入至所述振动状态判别器，输出判别结果，当所述判别结果为负时，判别为故障，并进行报警。
具体地，步骤S1的具体实现方式包括：步骤S1的具体实现方式包括：每个所述信号采集装置采集一个复合振动信号，组成n×m尺寸的复合信号矩阵X，n与m均为自然数。
具体地，步骤S2的具体实现方式包括：
S201，对所述复合信号矩阵按行去中心化，基于公式
Figure BDA0003315779960000044
获取去中心化处理后的复合信号矩阵Xd；ej为除第j行元素为1，其余均为零的n行列向量；Xj为复合信号矩阵X的第j行数据；
Figure BDA0003315779960000045
为复合信号矩阵第j行第i列数值；
S202，基于公式
Figure BDA0003315779960000046
计算协方差矩阵C；
S203，基于公式CV＝VΛ进行特征分解；V＝[v1,v2,…,vn]是特征向量按列组成的矩阵；Λ＝diag(ζ1,ζ2,…,ζn)是特征值组成的对角矩阵；
S204，进行白化处理，基于公式
Figure BDA0003315779960000051
获取复合信号矩阵Xw。
具体地，步骤S3的具体实现方式包括：
S301，通过拟牛顿迭代运算，计算最优复合矩阵W；
S302，基于公式S＝WT·X获取分离信号矩阵S。
具体地，步骤S301的具体实现方式包括：
S3011，初始化最优复合矩阵W的第i列，i的初始值为1，初始化为长度为m，向量元素和为1的随机列向量，得到迭代初始值
Figure BDA0003315779960000052
设定迭代阈值为∈；
S3012，进行第j次迭代运算，j的初始值为1，基于公式
Figure BDA0003315779960000053
Figure BDA0003315779960000054
获取第i个复合向量wi的第j次迭代值；
Figure BDA0003315779960000055
a1＝1；
S3013，标准化迭代数值，令
Figure BDA0003315779960000056
S3014，计算
Figure BDA0003315779960000057
若
Figure BDA0003315779960000058
认为迭代收敛，进入步骤S3015，否则返回步骤S3012；
S3015，将向量wi保存至最优复合矩阵W的第i列，对于i+1，返回步骤S3012，当i＝n时，迭代结束。
基于同一个技术构思，本发明还提供了一种多振源系统的故障诊断系统，包括多个独立振源，振动信号采集装置和计算机设备；所述振动信号采集装置用于采集复合振动信号；所述计算机设备被配置或编程为用于执行所述多振源系统的故障诊断方法。
与现有技术相比，本发明所具有的有益效果为：本发明针对多振源系统中的复合振动信号，使用盲源分离方法实现复合振动信号的分离，获取独立振源的振动信号，生成正常状态下分离信号数据集。使用该分离信号数据集训练对抗式生成网络，通过判别网络对的独立振源的振动信号进行辨识，判断该独立振源是否处于异常状态，从而实现无样本下故障源的定位和达到对故障的定位和异常预警功能。该方法适用性强，能够在缺少故障样本的情况下，充分利用不同独立振源的正常稳定状态工作振动信号训练模型，并在故障发生时进行及时的定位和异常预警。
附图说明
图1为本发明一实施例的一种多振源系统中无样本故障诊断方法的流程示意图。
图2为本发明一实施例中网络模型训练流程示意图。
图3为本发明一实施例中仿真正常运行状态下源信号的时域图。
图4为本发明一实施例中仿真正常运行状态下源信号的频域图。
图5为本发明一实施例中仿真正常运行状态下分离信号的时域图。
图6为本发明一实施例中仿真正常运行状态下分离信号的频域图。
图7为本发明一实施例中仿真故障运行状态下源信号的时域图。
图8为本发明一实施例中仿真故障运行状态下源信号的频域图。
图9为本发明一实施例中仿真故障运行状态下分离信号的时域图。
图10为本发明一实施例中仿真故障运行状态下分离信号的频域图。
图11为本发明一实施例中正常信号样本与故障信号样本的判别结果。
其中，部件A，部件B与部件C均为多振源系统中的独立振源。
具体实施方式
本发明一实施例的多振源系统，包括多个独立振源，所述多振源系统的故障诊断系统包括多个振动信号采集装置和计算机设备，所述振动信号采集装置为传感器，多个传感器用于采集所述多振源系统的复合振动信号，所述计算机设备被配置或编程为用于执行多振源系统的故障诊断方法，如图1所示，该方法包括以下步骤：
S1，使用多个传感器获取多振源系统的复合振动信号；设定有n个传感器，每个所述传感器采集一个复合振动信号，一共n个复合信号，每个复合信号的长度为m，组成n×m尺寸的复合信号矩阵X，n与m均为自然数。
S2，对所述复合振动信号进行前处理，以去除输入信号中的冗余信息。首先对所述复合信号矩阵X进行零均值化处理，然后使用PCA算法进行白化处理，具体实现方法包括以下步骤：
S201，按行去中心化，通过公式
Figure BDA0003315779960000061
计算获取去中心化处理后的复合信号矩阵Xd；
其中，Xj为信号矩阵X的第j行数据；
Figure BDA0003315779960000062
为信号矩阵第j行第i列数值；
ej为除第j行元素为1，其余均为零的n行列向量；
S202，去相关：通过公式
Figure BDA0003315779960000071
计算协方差矩阵C；其中，
Figure BDA0003315779960000072
为矩阵Xd的转置矩阵；
S203，基于公式CV＝VΛ进行特征分解，得到矩阵V和矩阵Λ；
其中，V＝[v1,v2,…,vn]是特征向量按列组成的矩阵；
Λ＝diag(ζ1,ζ2,…,ζn)是特征值组成的对角矩阵；
S204，根据公式
Figure BDA0003315779960000073
得到白化处理后的复合信号矩阵Xw；
其中，VT为矩阵V的转置矩阵；
S3，使用FastICA算法进行信号分离：从白化处理后的复合信号矩阵Xw分离出包含n个独立振源的振动信号的分离信号矩阵S，包括以下步骤：
S301，通过拟牛顿迭代运算，计算最优复合矩阵W，计算方法为：
S3011，初始化最优复合矩阵W第i列，i的初始值为1，初始化为长度为m，向量元素和为1的随机列向量，得到迭代初始值
Figure BDA0003315779960000074
设定迭代阈值为∈；
S3012，进行第j次迭代运算，j的初始值为1，令
Figure BDA0003315779960000075
Figure BDA0003315779960000076
其中
Figure BDA0003315779960000077
a1＝1；
S3013，标准化迭代数值，令
Figure BDA0003315779960000078
S3014，计算
Figure BDA0003315779960000079
若
Figure BDA00033157799600000710
则认为迭代收敛，进入步骤S3015，否则返回步骤S3012。
S3015，将向量wi保存至最优复合矩阵W的第i列，对于i+1,返回步骤S3012，当i＝n时，迭代结束。
S302，基于公式S＝WT·X，计算得到分离信号矩阵S；
S4，使用分离信号矩阵S制作分离信号数据集。具体实现方式包括：根据所述独立振源的转速和所述传感器的采样率，确定能够包含振动特征的样本长度l，然后将分离信号矩阵S切分为
Figure BDA00033157799600000711
个一维数组，对于每一个独立振源，输出单通道的波形.WAV文件，按照所代表的独立振源进行标记分组，以供模型训练使用。
S5，搭建生成对抗神经网络模型：该网络模型包括生成网络G和判别网络D；
生成网络G的输入为[batch_size,100]的随机张量，其网络结构包括：
第一层为全连接Dense层，连接单元数为16384，输出尺寸为[batch_size,16384]；
第二层为Reshape层，输出尺寸为[batch_size,16,1024]；且其后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第三层为至第六层为1-D转置卷积层，卷积核大小为25，步长为4，输出尺寸分别为[batch_size,64,512]、[batch_size,256,256]、[batch_size,1024,128]、[batch_size,4096,64]，且每层后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第七层为为1-D转置卷积层，卷积核大小为25，步长为4，输出尺寸为[batch_size,16384]；且其后接有一个批标准化单元和LeakyRelu激活单元。
所述判别网络D的输入为[batch_size,16384,channel_size]，其网络结构包括：
第一层为至第五层为1-D卷积层，卷积核大小为25，步长为4，输出尺寸分别为[batch_size,4096,64]、[batch_size,1024,128]、[batch_size,256,256]、[batch_size,64,512]、[batch_size,16,1024]；且每层后接一个LeakyReLu激活单元和一个PhaseShuffle相位重置单元；
第六层为Reshape层，输出尺寸为[batch_size,16384]；
第七层为全连接Dense层，输出为判别概率，尺寸为[batch_size,1]。
S6，训练网络模型。
如图2所示，训练过程包括编写训练函数，所述训练函数包含损失函数，数据加载、优化迭代模块。
S601，损失函数包含Wasserstein距离项和梯度惩罚项，Wasserstein又称为推土机距离(Earth-Move Distance,EM距离)，其定义为：
Figure BDA0003315779960000081
其中，∏(p,q)为分布p，q组合起来的所有可能的联合分布的集合，γ为服从∏(p,q)分布的γ样本，(x,y)服从联合分布γ，
Figure BDA0003315779960000082
为以(x,y)为输入的‖x-y‖的期望值，inf{}表示集合的下确界。
由于在计算生成样本分布pg与真实数据分布pr的是不现实的，因此将W(pr,pg)的计算方法转换为：
Figure BDA0003315779960000083
其中，xr为真实样本信号；
xf为生成样本信号；
pr为真实样本的概率分布；
pg为生成样本的概率分布；
Figure BDA0003315779960000091
为以xr为输入的D(xr)的期望；
Figure BDA0003315779960000092
为以xf为输入的D(xf)的期望；
D(xf)为满足一阶Lipschitz约束的判别网络函数。
梯度惩罚项是为了迫使判别网络函数满足一阶Lipschitz约束，惩罚项GP定义为：
Figure BDA0003315779960000093
其中，
Figure BDA0003315779960000094
为真实样本和生成样本的线性插值，其计算公式为；
Figure BDA0003315779960000095
其中k为线性插值参数，其取值范围为k∈[0,1]
得到损失函数
Figure BDA0003315779960000096
为：
Figure BDA0003315779960000097
S602，数据加载，具体包括：对于每一个独立振源数据集，添加分离信号数据集的文件路径，按路径提取样本文件。解码样本.WAV文件为张量数据，裁剪长度并归一化幅值。将所有数据打乱顺序，然后按一定数量进行采样打包，供迭代器使用。
S603，，迭代训练，如图2所示，训练过程具体包括：
S6031，设定收敛条件为最对于每一个独立振源大迭代次数N，随机化初始化生成器网络参数θ，判别网络参数φ；
S6032，进行次数为m的判别训网络训练：为从所述分离信号数据集中采集一批真实样本xr，输入判别网络输出D(xr；φ)，再从N(0,1)分布中采样一批随机向量h，输入生成器输出xf＝G(h；θ)，随后计算损失函数，使用Adam优化器使，设置超参数α,β1,β2，进行参数φ优化并更新，使损失函数值
Figure BDA00033157799600000910
最大化：
Figure BDA0003315779960000098
S6033，进行生成器训练：采样N(0,1)分布输入生成器，计算损失函数，使用Adam优化器，设置超参数α1,β1,β2，进行参数θ的优化并更新，使损失函数
Figure BDA0003315779960000099
的值最小化：
Figure BDA0003315779960000101
训练迭代直至达到最大迭代次数N。
S7，提取收敛后判别网络的参数φ，构造输入已经搭好的网络模型中，对于每一个独立振源，形成振动状态判别器。
S8，基于步骤S1至步骤S3，对于每一个独立振源，周期性获取分离信号矩阵，输入至所述振动状态判别器，输出判别结果，当所述判别结果为负时，判别为故障，并进行报警。
以列车转向架列为多振源系统的案例，所述车转向架包括三个独立振源，分别为轴承、电机和齿轮箱，使用三个传感器采集所述转向架的复合振动信号，使用本发明的多振源系统的故障诊断方法进行仿真，具体包括：
电机的振动信号的主要成分有转子转速的一倍频率、二倍频率与少量高频信号，当定子线圈松动时会产生高幅值的线圈通过频率振动；齿轮箱的振动信号主要包含转速的一倍频率、二倍频率和齿轮啮合频率，当发生故障是会发生一定程度的频率调制和幅度调制现象；所述轴承的振动信号主要是其内外圈的转速倍频还有少量的零件共振高频成分，当其发生故障时会产生一定频率范围内的高频冲激信号。
根据所述轴承、电机和齿轮箱的振动特性建立振动仿真信号，得到各部件正常状态下信号函数如下：
x＝10sin(152×2π×t)+5sin(304×2π×t)+2sin(700×2π×t)+cos(800×2π×t)
y＝5sin(496×2π×t)+2sin(993×2π×t)+2cos(1440×2π×t)
z＝10sin(273×2π×t)+5sin(446×2π×t)+2sin(993×2π×t)+cos(1986×2π×t)
x为电机的正常信号函数，y为齿轮箱的正常信号函数，z为滚动轴承的正常信号函数；
选取轴承、电机和齿轮箱的部分故障特征添加入正常信号函数中形成故障信号：
x_fault＝10sin(152×2π×t)+5sin(304×2π×t)+2sin(700×2π×t)+cos(800×2π×t)+10sin(750×2π×t)
y_fault＝5sin(496×2π×t+5sin(50×2π×t))+2sin(993×2π×t)+2cos(1440×2π×t)
Figure BDA0003315779960000111
x_fault为电机的故障信号函数，y_fault为齿轮箱的故障信号函数，z_fault为滚动轴承的故障信号函数；
分别使用上述正常信号函数和故障信号函数生成六条具有8.192×106个采样点，采样率为16kHz，共计100s的仿真数据。
随后使用线性复合方式模拟信号的复合过程，将三个独立源信号复合成为三个复合信号，以模拟传感器的观测过程。其复合信号矩阵为：
Figure BDA0003315779960000112
本案例中独立振源的数量与传感器数量相同，即独立振源的数量与观测信号源数量相同，且数据本身是零均值化的，所以无需进行数据前处理，直接使用FastICA算法对进行分离。
如图3-图10所示，基于分离前后的正常信号的时域波形图与频谱图，对比源信号可以发现，特征频率在分离前后变化较小，分离效果较好。
将分离后的三条正常信号进行归一化处理，并分别等距切分500份，得到1500条正常样本数据，分别为轴承的正常样本数据a，电机的正常样本数据b和齿轮箱的正常样本数据c。
将分离后的三条故障信号进行归一化处理，并分别等距切分500份，得到1500条故障样本数据，分别为轴承的故障样本a_fault、电机的故障样本b_fault和齿轮箱的故障样本c_fault。
将上述六个据按比特编码输出为.wav文件保存，采样率为16kHz，比特位深为16。将正常样本数据输入搭建好的一维生成对抗神经网络模型中进行模型训练，在经过15010次循环后达到收敛。
提取收敛模型中的判别网络对振动信号进行状态判别。输入500个轴承的正常样本数据a，其输出值的范围在[8.1,19.1]；输入500个电机的正常样本数据b，其输出值的范围在[5.1,17.2]；输入500个齿轮箱的正常样本数据c，其输出值的范围在[21.8,26.1]；输入500个轴承的故障样本a_fault，其输出值的范围在[-34.4,-33.1]；输入500个电机的故障样本b_fault，其输出值的范围在[-7.5,-2.2]；输入500个齿轮箱的故障样本c_fault，其输出值的范围在[-17.3,-14.9]。
如图11所示，将输出值归一化后的分布中，故障与正常信号有着明显的界限。因此判别网络对故障信号有着强大的分辨能力，当判别网络输出为负时，认为该部件故障，并进行预警。

Claims (10)

1.一种多振源系统的故障诊断方法，所述多振源系统包括多个独立振源，其特征在于，所述多振源系统还包括n个数据采集装置，所述数据采集装置的采样信号长度为m，包括以下步骤：
A，获取分离信号数据集；搭建生成对抗神经网络模型；
分离信号数据集的具体实现方式包括：
S1，获取多振源系统的复合振动信号；
S2，基于所述复合振动信号获取复合信号矩阵，对所述复合信号矩阵依次进行零均值化与白化处理，获取白化处理后的复合信号矩阵；
S3，基于白化处理后的复合信号矩阵获取分离信号矩阵；
S4，将所述分离信号矩阵切分为
Figure FDA0003315779950000011
个一维数组，输出单通道的波形.WAV文件，对于每一个独立振源，制作分离信号数据集；l为含振动特征的样本长度，根据独立振源的转速和所述信号采集装置的采样率确定；
所述生成对抗神经网络模型包括生成网络G和判别网络D；
所述生成网络G的网络结构包括：
第一层为全连接Dense层；
第二层为Reshape层，且其后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第三层至第六层均为1-D转置卷积层，卷积核大小为25，步长为4，且每层后接有一个批标准化单元，一个LeakyRelu激活单元和一个tanh激活单元；
第七层为为1-D转置卷积层，卷积核大小为25，步长为4，且接有一个批标准化单元和LeakyRelu激活单元；
所述判别网络D的网络结构包括：
第一层为至第五层均为1-D卷积层，卷积核大小为25，步长为4，且每层后接一个LeakyReLu激活单元和一个PhaseShuffle相位重置单元；
第六层为Reshape层；
第七层为全连接Dense层；
B，对于每一个独立振源，基于所述分离信号数据集和所述生成对抗神经网络模型进行迭代训练；
C，对于每一个独立振源，获取振动状态判别器。
2.如权利要求1所述的多振源系统的故障诊断方法，其特征在于，步骤B的具体实现方式包括：
S601，获取损失函数
Figure FDA0003315779950000021
Figure FDA0003315779950000022
Figure FDA0003315779950000023
为真实样本和生成样本的线性插值，其计算公式为：
Figure FDA0003315779950000024
xr为真实样本信号；k为线性插值的权重，取值范围为k∈[0，1]；xf为生成样本信号；pr为真实样本的概率分布；pg为生成样本的概率分布；D(xr)为以真实样本信号xr为输入且满足一阶Lipschitz约束的判别网络函数，D(xf)为以生成样本信号xf为输入且满足一阶Lipschitz约束的判别网络函数；
Figure FDA0003315779950000025
为以真实样本和生成样本的线性插值
Figure FDA0003315779950000026
为输入且满足一阶Lipschitz约束的判别网络函数；
Figure FDA0003315779950000027
为判别网络函数D(xr)的期望；
Figure FDA0003315779950000028
为判别网络函数D(xf)的期望；
Figure FDA0003315779950000029
为函数
Figure FDA00033157799500000210
的期望；λ为梯度惩罚项系数；
S602，解码样本WAV文件为张量数据，裁剪长度并归一化幅值，打乱顺序后进行采样打包；
S603，对步骤S602获取的数据进行迭代训练。
3.如权利要求2所述的多振源系统的故障诊断方法，其特征在于，获取所述损失函数的方法包括：
S6011，获取Wasserstein距离项W(pr，pg)，计算方法为：
Figure FDA00033157799500000211
θ为生成网络参数；
S6012，获取惩罚项GP的计算公式为：
Figure FDA00033157799500000212
4.如权利要求3所述的用于多振源系统的故障诊断方法，其特征在于，步骤S603中，迭代训练的具体实现方式包括：
S6031，设定收敛条件为最大迭代次数N，随机化初始化生成器网络参数θ，判别网络参数φ；
S6032，进行次数为m的判别训网络训练：从所述分离信号数据集中采集一批真实样本xr，输入判别网络输出D(xr；φ)，再从N(0，1)分布中采样一批随机向量h，输入生成器输出xf＝G(h；θ)，随后计算损失函数，使用Adam优化器，设置超参数α，β1，β2，进行参数φ优化并更新，使损失函数值
Figure FDA0003315779950000031
最大化：
Figure FDA0003315779950000032
S6033，进行生成器训练：采样N(0，1)分布输入生成器，计算损失函数，使用Adam优化器，设置超参数α1，β1，β2，进行参数θ的优化并更新，使损失函数
Figure FDA0003315779950000033
的值最小化：
Figure FDA0003315779950000034
S6034，训练迭代直至达到最大迭代次数N。
5.如权利要求1～4任一所述的用于多振源系统的故障诊断方法，其特征在于，基于步骤S1至步骤S3，周期性获取分离信号矩阵，输入至所述振动状态判别器，输出判别结果，当所述判别结果为负时，判别为故障，并进行报警。
6.如权利要求5所述的用于多振源系统的故障诊断方法，其特征在于，步骤S1的具体实现方式包括：每个所述信号采集装置采集一个复合振动信号，组成n×m尺寸的复合信号矩阵X，n与m均为自然数。
7.如权利要求5所述的用于多振源系统的故障诊断方法，其特征在于，步骤S2的具体实现方式包括：
S201，对所述复合信号矩阵按行去中心化，基于公式
Figure FDA0003315779950000035
获取去中心化处理后的复合信号矩阵Xd；ej为除第j行元素为1，其余均为零的n行列向量；Xj为复合信号矩阵X的第j行数据；
Figure FDA0003315779950000036
为复合信号矩阵X第j行第i列数值；
S202，基于公式
Figure FDA0003315779950000037
计算协方差矩阵C；
S203，基于公式CV＝VA进行特征分解；V＝[v1，v2，...，vn]是特征向量按列组成的矩阵；Λ＝diag(ζ1，ζ2，...，ζn)是特征值组成的对角矩阵；
S204，进行白化处理，基于公式
Figure FDA0003315779950000038
获取白化处理后的复合信号矩阵Xw。
8.如权利要求5所述的用于多振源系统的故障诊断方法，其特征在于，步骤S3的具体实现方式包括：
S301，通过拟牛顿迭代运算，计算最优复合矩阵W；
S302，基于公式S＝WT·X获取分离信号矩阵S。
9.如权利要求8所述的用于多振源系统的故障诊断方法，其特征在于，步骤S301的具体实现方式包括：
S3011，初始化最优复合矩阵W的第i列，i的初始值为1，初始化为长度为m，向量元素和为1的随机列向量，得到迭代初始值
Figure FDA0003315779950000041
设定迭代阈值为∈；
S3012，进行第j次迭代运算，j的初始值为1，基于公式
Figure FDA0003315779950000042
Figure FDA0003315779950000043
获取第i个复合向量wi的第j次迭代值；
Figure FDA0003315779950000044
a1＝1；
S3013，标准化迭代数值，令
Figure FDA0003315779950000045
S3014，计算
Figure FDA0003315779950000046
若
Figure FDA0003315779950000047
认为迭代收敛，进入步骤S3015，否则返回步骤S3012；
S3015，将向量wi保存至最优复合矩阵W的第i列，对于i+1，返回步骤S3012，当i＝n时，迭代结束。
10.一种多振源系统的故障诊断系统，包括多个独立振源，其特征在于，还包括振动信号采集装置和计算机设备；所述振动信号采集装置用于采集复合振动信号；所述计算机设备被配置或编程为用于执行权利要求1～9任一所述的方法。
CN202111230641.2A 2021-10-22 2021-10-22 一种多振源系统的故障诊断系统及方法 Active CN114136565B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111230641.2A CN114136565B (zh) 2021-10-22 2021-10-22 一种多振源系统的故障诊断系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111230641.2A CN114136565B (zh) 2021-10-22 2021-10-22 一种多振源系统的故障诊断系统及方法

Publications (2)

Publication Number Publication Date
CN114136565A true CN114136565A (zh) 2022-03-04
CN114136565B CN114136565B (zh) 2022-08-09

Family

ID=80394496

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111230641.2A Active CN114136565B (zh) 2021-10-22 2021-10-22 一种多振源系统的故障诊断系统及方法

Country Status (1)

Country Link
CN (1) CN114136565B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140079248A1 (en) * 2012-05-04 2014-03-20 Kaonyx Labs LLC Systems and Methods for Source Signal Separation
CN110567720A (zh) * 2019-08-07 2019-12-13 东北电力大学 非平衡小样本场景下风机轴承故障深度对抗诊断方法
CN111898634A (zh) * 2020-06-22 2020-11-06 西安交通大学 一种基于深度对抗域自适应的智能故障诊断方法
CN112330538A (zh) * 2020-09-29 2021-02-05 电子科技大学 一种基于特征点优化提取的损伤温度重构图像拼接方法
CN113095402A (zh) * 2021-04-12 2021-07-09 江南大学 一种基于编码输入的生成对抗网络故障检测方法及系统
CN113375941A (zh) * 2021-06-08 2021-09-10 西安交通大学 高速动车组轴承的开集故障诊断方法
CN113378655A (zh) * 2021-05-24 2021-09-10 电子科技大学 一种基于深度神经网络的对抗性能量分解方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140079248A1 (en) * 2012-05-04 2014-03-20 Kaonyx Labs LLC Systems and Methods for Source Signal Separation
CN110567720A (zh) * 2019-08-07 2019-12-13 东北电力大学 非平衡小样本场景下风机轴承故障深度对抗诊断方法
CN111898634A (zh) * 2020-06-22 2020-11-06 西安交通大学 一种基于深度对抗域自适应的智能故障诊断方法
CN112330538A (zh) * 2020-09-29 2021-02-05 电子科技大学 一种基于特征点优化提取的损伤温度重构图像拼接方法
CN113095402A (zh) * 2021-04-12 2021-07-09 江南大学 一种基于编码输入的生成对抗网络故障检测方法及系统
CN113378655A (zh) * 2021-05-24 2021-09-10 电子科技大学 一种基于深度神经网络的对抗性能量分解方法
CN113375941A (zh) * 2021-06-08 2021-09-10 西安交通大学 高速动车组轴承的开集故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SHENGNAN TANG ET AL.: "Data Preprocessing Techniques in Convolutional Neural Network Based on Fault Diagnosis Towards Rotating Machinery", 《IEEE ACCESS》 *
徐金东 等: "注意机制生成对抗网络单通道盲图像分离", 《中国科技论文》 *

Also Published As

Publication number Publication date
CN114136565B (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
CN109580215B (zh) 一种基于深度生成对抗网络的风电传动系统故障诊断方法
CN112200244B (zh) 一种基于递阶对抗训练的航天发动机异常智能检测方法
CN109555566B (zh) 一种基于lstm的汽轮机转子故障诊断方法
CN107657250B (zh) 轴承故障检测及定位方法及检测定位模型实现系统和方法
Gao et al. ASM1D-GAN: An intelligent fault diagnosis method based on assembled 1D convolutional neural network and generative adversarial networks
CN111582551B (zh) 风电场短期风速预测方法、系统及电子设备
CN110657984A (zh) 一种基于强化胶囊网络的行星齿轮箱故障诊断方法
CN105973595A (zh) 一种滚动轴承故障的诊断方法
CN111855202A (zh) 一种齿轮箱故障诊断方法及系统
CN114548199A (zh) 一种基于深度迁移网络的多传感器数据融合方法
CN113255458A (zh) 基于多视图关联特征学习的轴承故障诊断方法
CN113642508A (zh) 基于参数自适应vmd与优化svm的轴承故障诊断方法
CN114897138A (zh) 基于注意力机制和深度残差网络的系统故障诊断方法
CN109615109A (zh) 基于小样本数据的深度学习风电告警信息分析方法
CN118296470A (zh) 一种多分类器融合的智能分类方法
CN117171544B (zh) 基于多通道融合卷积神经网络的电机振动故障诊断方法
CN113361782A (zh) 基于改进mkpls的光伏发电功率短期滚动预测方法
CN117828531A (zh) 一种基于多传感器多尺度特征融合的轴承故障诊断方法
CN114383845B (zh) 一种基于嵌入式零样本学习模型的轴承复合故障诊断方法
CN114136565B (zh) 一种多振源系统的故障诊断系统及方法
CN115356599B (zh) 一种多模态城市电网故障诊断方法及系统
Priyadarshini et al. Detection of islanding and non-islanding fault disturbances in microgrid using LMD and deep stacked RVFLN based auto-encoder
GB2623358A (en) Method and system for fault diagnosis of nuclear power circulating water pump based on optimized capsule network
CN115456013A (zh) 基于运行数据的风电机组变流器功率模块故障诊断方法
CN115130500A (zh) 一种基于sdae算法的电机故障诊断系统及方法

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