CN116026593A - Cross-working-condition rolling bearing fault targeted migration diagnosis method and system - Google Patents

Cross-working-condition rolling bearing fault targeted migration diagnosis method and system Download PDF

Info

Publication number
CN116026593A
CN116026593A CN202211631915.3A CN202211631915A CN116026593A CN 116026593 A CN116026593 A CN 116026593A CN 202211631915 A CN202211631915 A CN 202211631915A CN 116026593 A CN116026593 A CN 116026593A
Authority
CN
China
Prior art keywords
rolling bearing
feature
loss
cross
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.)
Pending
Application number
CN202211631915.3A
Other languages
Chinese (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN202211631915.3A priority Critical patent/CN116026593A/en
Publication of CN116026593A publication Critical patent/CN116026593A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明提供了一种跨工况的滚动轴承故障靶向迁移诊断方法及系统,解决了传统滚动轴承故障诊断算法难以提取源域和目标域中网络深层特征信息、无法实现有效的跨域故障诊断的问题。本发明利用特征编码器从输入的滚动轴承振动信号中精确地提取信号的高维映射特征;将该特征进一步输入到图构建层,挖掘数据的深层特征,并利用多通道核图卷积网络对实例图建模;利用基于差异和对抗的训练来最小化源域和目标域分布之间的距离,分类器则使用提取的域不变特征来完成跨域故障识别。本发明与其他方法相比较,在滚动轴承跨工况条件下,可以更好的提取深层特征用于跨域传递,极大的提高了诊断准确率。

Figure 202211631915

The present invention provides a rolling bearing fault targeted migration diagnosis method and system across working conditions, which solves the problem that the traditional rolling bearing fault diagnosis algorithm is difficult to extract deep network feature information in the source domain and the target domain, and cannot realize effective cross-domain fault diagnosis. . The present invention uses a feature encoder to accurately extract the high-dimensional mapping feature of the signal from the input rolling bearing vibration signal; this feature is further input into the graph construction layer, the deep features of the data are mined, and the multi-channel kernel image convolution network is used to analyze Graph modeling; using difference-based and adversarial training to minimize the distance between the source and target domain distributions, the classifier uses the extracted domain-invariant features to complete cross-domain fault identification. Compared with other methods, the present invention can better extract deep features for cross-domain transmission under the condition of cross-working conditions of rolling bearings, and greatly improve the diagnosis accuracy.

Figure 202211631915

Description

跨工况的滚动轴承故障靶向迁移诊断方法及系统Rolling bearing fault targeted migration diagnosis method and system across working conditions

技术领域Technical Field

本发明涉及轴承故障诊断技术领域,特别涉及一种跨工况的滚动轴承故障靶向迁移诊断方法及系统。The present invention relates to the technical field of bearing fault diagnosis, and in particular to a rolling bearing fault targeted migration diagnosis method and system across working conditions.

背景技术Background Art

本部分的陈述仅仅是提供了与本发明相关的背景技术,并不必然构成现有技术。The statements in this section merely provide background art related to the present invention and do not necessarily constitute prior art.

滚动轴承作为大多数旋转机械设备的关键部位,由于工作环境恶劣、工况发生变化、超负荷运转时间长等特点,经常发生故障,进而带来整机机械故障并造成一定的经济损失。因此,对跨工况下滚动轴承的运行状态进行监测具有重要意义。Rolling bearings are the key parts of most rotating mechanical equipment. Due to the characteristics of harsh working environment, changing working conditions, long overload operation time, etc., they often fail, which leads to mechanical failure of the whole machine and causes certain economic losses. Therefore, it is of great significance to monitor the operating status of rolling bearings under different working conditions.

近年来,伴随着海量数据的出现,基于深度学习的数据驱动方法备受关注。它不同于传统的浅层网络架构,而是通过堆叠多层非线性处理单元来实现的,并且提供了一种端到端的解决方式,众多学者也进行了相应的研究。Jia等人提出了一种基于SAE的深度神经网络DNN,用于识别电机和齿轮箱中的故障。Ding等人提了使用卷积神经网络从小波包能量图像中挖掘能量波动的多尺度特征,用于主轴轴承的故障诊断。Chen等人提出了一种利用原始振动信号作为输入的自动特征学习神经网络,并使用两个具有不同核大小的CNN从原始数据中自动提取不同频率信号特征,然后根据学习到的特征,使用LSTM来识别故障类型。上述的深度学习算法已经在恒定工况下滚动轴承的故障诊断领域取得了良好的效果,然而在复杂的跨工况条件下却很难提取域间差异特性,使得模型的泛化性能有所下降。In recent years, with the emergence of massive data, data-driven methods based on deep learning have attracted much attention. It is different from the traditional shallow network architecture. It is realized by stacking multiple layers of nonlinear processing units and provides an end-to-end solution. Many scholars have also conducted corresponding research. Jia et al. proposed a deep neural network DNN based on SAE to identify faults in motors and gearboxes. Ding et al. proposed the use of convolutional neural networks to mine multi-scale features of energy fluctuations from wavelet packet energy images for fault diagnosis of spindle bearings. Chen et al. proposed an automatic feature learning neural network using raw vibration signals as input, and used two CNNs with different kernel sizes to automatically extract different frequency signal features from raw data, and then used LSTM to identify fault types based on the learned features. The above deep learning algorithms have achieved good results in the field of rolling bearing fault diagnosis under constant working conditions. However, it is difficult to extract inter-domain difference characteristics under complex cross-working conditions, which reduces the generalization performance of the model.

迁移学习作为一种减少跨域特征分布差异的方法,为建立从源域标记数据到目标域未标记数据之间的知识迁移提供了新的思路,有效的解决了跨工况下域间分布差异的问题。在智能故障诊断领域,已经提出了许多迁移学习方法来解决跨域诊断问题,它们基本可以分为基于实例的,基于模型的和基于特征的。Xiao等人使用TrAdaBoost来调整每个训练样本的权重因子,以此来增强故障分类器的诊断性能。Wang等人提出了一种基于估计伪标签的条件MMD,以缩短轴承故障诊断的分布距离。通过最小化MMD损失,在多个层中同时对齐边缘分布和条件分布。Han等人提出了一种深度对抗卷积神经网络DACNN,利用基于对抗的损失函数来缩小域间差异,用于提高齿轮箱和电机故障诊断的泛化性能。Li等人将具有拓扑结构的图数据作为输入,采用对数据关系挖掘更有效,特征表示更强大的图卷积网络GCN建模用于机械故障诊断,取得了优异的性能。可以看出,类别标签、领域标签和数据结构的信息在减少源域和目标域差异之间起着重要的作用,并且它们应该是互相完善,互相增强的。As a method to reduce the difference in cross-domain feature distribution, transfer learning provides a new idea for establishing knowledge transfer from source domain labeled data to target domain unlabeled data, and effectively solves the problem of inter-domain distribution differences under cross-working conditions. In the field of intelligent fault diagnosis, many transfer learning methods have been proposed to solve the cross-domain diagnosis problem, which can be basically divided into instance-based, model-based and feature-based. Xiao et al. used TrAdaBoost to adjust the weight factor of each training sample to enhance the diagnostic performance of the fault classifier. Wang et al. proposed a conditional MMD based on estimated pseudo-labels to shorten the distribution distance of bearing fault diagnosis. By minimizing the MMD loss, the marginal distribution and conditional distribution are aligned simultaneously in multiple layers. Han et al. proposed a deep adversarial convolutional neural network (DACNN) that uses adversarial loss functions to reduce inter-domain differences and improve the generalization performance of gearbox and motor fault diagnosis. Li et al. used graph data with topological structure as input and adopted graph convolutional network (GCN) modeling, which is more effective in data relationship mining and more powerful in feature representation, for mechanical fault diagnosis, and achieved excellent performance. It can be seen that the information of category labels, domain labels, and data structures plays an important role in reducing the differences between the source domain and the target domain, and they should complement and enhance each other.

然而,发明人发现,现有的方法仅仅考虑了源域和目标域的两种信息,未把数据结构融入深度神经网络模型之中,不同工况下数据分布差异较大,导致故障识别准确率低和泛化性能不足。However, the inventors found that the existing methods only considered two types of information, the source domain and the target domain, and did not integrate the data structure into the deep neural network model. The data distribution under different working conditions was quite different, resulting in low fault identification accuracy and insufficient generalization performance.

发明内容Summary of the invention

为了解决现有技术的不足,本发明提供了一种跨工况的滚动轴承故障靶向迁移诊断方法及系统,在滚动轴承跨工况条件下,可以更好的提取深层特征用于跨域传递,极大的提高了诊断准确率。In order to address the deficiencies in the prior art, the present invention provides a rolling bearing fault targeted migration diagnosis method and system across operating conditions. Under the rolling bearing cross-operating conditions, deep features can be better extracted for cross-domain transfer, greatly improving the diagnostic accuracy.

为了实现上述目的,本发明采用如下技术方案:In order to achieve the above object, the present invention adopts the following technical solution:

本发明第一方面提供了一种跨工况的滚动轴承故障靶向迁移诊断方法。A first aspect of the present invention provides a rolling bearing fault targeted migration diagnosis method across operating conditions.

一种跨工况的滚动轴承故障靶向迁移诊断方法,包括以下过程:A rolling bearing fault targeted migration diagnosis method across working conditions includes the following processes:

获取滚动轴承的振动信号;Obtain vibration signals of rolling bearings;

通过特征编码器、图构建层和多通道核图卷积网络提取振动信号中的深层高维特征;Extract deep high-dimensional features from vibration signals through feature encoders, graph construction layers, and multi-channel kernel graph convolutional networks;

根据提取的深层高维特征和分类器,得到滚动轴承故障诊断结果;According to the extracted deep high-dimensional features and classifiers, the rolling bearing fault diagnosis results are obtained;

其中,根据源域训练集结合分类器得到的深层高维特征得到分类损失,根据源域训练集和目标域训练集得到的深层高维特征得到结构差异损失,根据源域训练集和目标域训练集得到的深层高维特征以及对抗网络,得到对抗损失;Among them, the classification loss is obtained according to the deep high-dimensional features obtained by combining the source domain training set with the classifier, the structural difference loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set, and the adversarial loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set and the adversarial network;

以分类损失、第一参数与结构差异损失的乘积、第二参数与对抗损失的乘积三者的加和为总体损失函数,根据总体目标函数,通过反向传播算法优化特征提取器、分类器和判别器的参数。The overall loss function is the sum of the classification loss, the product of the first parameter and the structural difference loss, and the product of the second parameter and the adversarial loss. According to the overall objective function, the parameters of the feature extractor, classifier, and discriminator are optimized through the back propagation algorithm.

作为本发明第一方面可选的一种实现方式,图构建层用于获取邻接矩阵,包括:As an optional implementation of the first aspect of the present invention, the graph construction layer is used to obtain an adjacency matrix, including:

根据特征编码器网络,从样本数据中获得高维特征映射,即X=G(x);According to the feature encoder network, a high-dimensional feature map is obtained from the sample data, that is, X = G(x);

提取的高维特征映射,输入到线性层中,经过Softmax后表示为

Figure BDA0004006128250000031
The extracted high-dimensional feature map is input into the linear layer and expressed as
Figure BDA0004006128250000031

通过线性层的特征与其转置之间进行矩阵相乘计算,获得邻接矩阵A,

Figure BDA0004006128250000032
The adjacency matrix A is obtained by performing matrix multiplication between the features of the linear layer and its transpose.
Figure BDA0004006128250000032

通过KNN算法,构造边关系:

Figure BDA0004006128250000033
Through the KNN algorithm, construct edge relationships:
Figure BDA0004006128250000033

作为本发明第一方面进一步的限定,多通道核图卷积网络,包括:As a further limitation of the first aspect of the present invention, the multi-channel kernel graph convolutional network includes:

Figure BDA0004006128250000034
Figure BDA0004006128250000034

Figure BDA0004006128250000041
Figure BDA0004006128250000041

其中,X代表输入,A代表邻接矩阵,

Figure BDA0004006128250000042
代表可训练权值矩阵,G代表多通道核图卷积操作,
Figure BDA0004006128250000043
代表第ki个通道在第L层的高维特征表示,[·]代表特征拼接,H表示经过多通道核图卷积网络后的输出特征。Among them, X represents the input, A represents the adjacency matrix,
Figure BDA0004006128250000042
represents the trainable weight matrix, G represents the multi-channel kernel graph convolution operation,
Figure BDA0004006128250000043
represents the high-dimensional feature representation of the k i -th channel at the L layer, [ ] represents feature concatenation, and H represents the output feature after passing through the multi-channel kernel graph convolutional network.

作为本发明第一方面可选的一种实现方式,交叉熵损失LC,包括:As an optional implementation of the first aspect of the present invention, the cross entropy loss LC includes:

Figure BDA0004006128250000044
Figure BDA0004006128250000044

其中,

Figure BDA0004006128250000045
表示分类器的预测结果,E表示数学期望值,
Figure BDA0004006128250000046
为源域样本,
Figure BDA0004006128250000047
为其标签。in,
Figure BDA0004006128250000045
represents the prediction result of the classifier, E represents the mathematical expectation,
Figure BDA0004006128250000046
is the source domain sample,
Figure BDA0004006128250000047
Label it.

作为本发明第一方面可选的一种实现方式,结构差异损失Ls,包括:As an optional implementation of the first aspect of the present invention, the structural difference loss L s includes:

Figure BDA0004006128250000048
Figure BDA0004006128250000048

其中,

Figure BDA0004006128250000049
Figure BDA00040061282500000410
分别表示第i个源域样本和第j个目标域样本经过特征提取器后的特征映射,φ表示非线性特征映射,Ωk表示嵌入提取的特征到再生核希尔伯特空间RKHS中的距离度量,采用m个核的凸组合ku来对映射进行有效地估计:
Figure BDA00040061282500000411
其中,αu是不同核的加权参数,且
Figure BDA00040061282500000412
E表示数学期望值,
Figure BDA00040061282500000413
为源域样本,
Figure BDA00040061282500000414
为目标域样本。in,
Figure BDA0004006128250000049
and
Figure BDA00040061282500000410
They represent the feature maps of the i-th source domain sample and the j-th target domain sample after the feature extractor, φ represents the nonlinear feature map, Ω k represents the distance metric embedded in the extracted features to the reproducing kernel Hilbert space RKHS, and the convex combination of m kernels k u is used to effectively estimate the mapping:
Figure BDA00040061282500000411
where αu is the weighting parameter of different kernels, and
Figure BDA00040061282500000412
E represents the mathematical expectation,
Figure BDA00040061282500000413
is the source domain sample,
Figure BDA00040061282500000414
is the target domain sample.

作为本发明第一方面可选的一种实现方式,对抗损失LAD,包括:As an optional implementation of the first aspect of the present invention, counteracting the loss L AD includes:

Figure BDA00040061282500000415
Figure BDA00040061282500000415

其中,D(·)为经过判别器后的特征输出,

Figure BDA00040061282500000416
Figure BDA00040061282500000417
分别表示第i个源域样本和第j个目标域样本经过特征提取器后的特征映射,E表示数学期望值,
Figure BDA00040061282500000418
为源域样本,
Figure BDA00040061282500000419
为目标域样本。Among them, D(·) is the feature output after the discriminator,
Figure BDA00040061282500000416
and
Figure BDA00040061282500000417
They represent the feature maps of the i-th source domain sample and the j-th target domain sample after passing through the feature extractor, E represents the mathematical expectation value,
Figure BDA00040061282500000418
is the source domain sample,
Figure BDA00040061282500000419
is the target domain sample.

作为本发明第一方面可选的一种实现方式,通过反向传播算法优化特征提取器、分类器和判别器的参数,包括:As an optional implementation of the first aspect of the present invention, the parameters of the feature extractor, the classifier and the discriminator are optimized by a back propagation algorithm, including:

Figure BDA0004006128250000051
Figure BDA0004006128250000051

其中,

Figure BDA0004006128250000052
代表偏微分算子,η代表学习率,θF代表特征提取器的参数,θC代表分类器的参数,θD代表判别器的参数,LC为交叉熵损失,Ls为结构差异损失,LAD为对抗损失。in,
Figure BDA0004006128250000052
represents the partial differential operator, η represents the learning rate, θ F represents the parameters of the feature extractor, θ C represents the parameters of the classifier, θ D represents the parameters of the discriminator, LC is the cross entropy loss, Ls is the structural difference loss, and LAD is the adversarial loss.

本发明第二方面提供了一种跨工况的滚动轴承故障靶向迁移诊断系统。A second aspect of the present invention provides a rolling bearing fault targeted migration diagnosis system across operating conditions.

一种跨工况的滚动轴承故障靶向迁移诊断系统,包括:A rolling bearing fault targeted migration diagnosis system across working conditions, comprising:

数据获取模块,被配置为:获取滚动轴承的振动信号;The data acquisition module is configured to: acquire a vibration signal of the rolling bearing;

特征提取模块,被配置为:通过特征编码器、图构建层和多通道核图卷积网络提取振动信号中的深层高维特征;A feature extraction module, configured to: extract deep high-dimensional features from the vibration signal through a feature encoder, a graph construction layer, and a multi-channel kernel graph convolutional network;

故障诊断模块,被配置为:根据提取的深层高维特征和分类器,得到滚动轴承故障诊断结果;The fault diagnosis module is configured to: obtain a rolling bearing fault diagnosis result based on the extracted deep high-dimensional features and the classifier;

其中,根据源域训练集结合分类器得到的深层高维特征得到分类损失,根据源域训练集和目标域训练集得到的深层高维特征得到结构差异损失,根据源域训练集和目标域训练集得到的深层高维特征以及对抗网络,得到对抗损失;Among them, the classification loss is obtained according to the deep high-dimensional features obtained by combining the source domain training set with the classifier, the structural difference loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set, and the adversarial loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set and the adversarial network;

以分类损失、第一参数与结构差异损失的乘积、第二参数与对抗损失的乘积三者的加和为总体损失函数,根据总体目标函数,通过反向传播算法优化特征提取器、分类器和判别器的参数。The overall loss function is the sum of the classification loss, the product of the first parameter and the structural difference loss, and the product of the second parameter and the adversarial loss. According to the overall objective function, the parameters of the feature extractor, classifier, and discriminator are optimized through the back propagation algorithm.

本发明第三方面一种计算机可读存储介质,其上存储有程序,该程序被处理器执行时实现如本发明第一方面所述的跨工况的滚动轴承故障靶向迁移诊断方法中的步骤。A third aspect of the present invention is a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps in the method for targeted migration diagnosis of rolling bearing faults across operating conditions as described in the first aspect of the present invention.

本发明第四方面提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如本发明第一方面所述跨工况的滚动轴承故障靶向迁移诊断方法中的步骤。The fourth aspect of the present invention provides an electronic device, comprising a memory, a processor, and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the method for targeted migration diagnosis of rolling bearing faults across working conditions as described in the first aspect of the present invention are implemented.

与现有技术相比,本发明的有益效果是:Compared with the prior art, the present invention has the following beneficial effects:

1、本发明采用特征编码器从输入信号中自适应地提取信号特征,并且利用图构建层来获取特征编码器捕获特征中的数据结构,进而构造实例图,并且应用多通道核图卷积网络对其建模,进一步挖掘信号的高维特征,解决了跨工况下源域和目标域深层特征提取困难的问题。1. The present invention adopts a feature encoder to adaptively extract signal features from the input signal, and uses a graph construction layer to obtain the data structure in the features captured by the feature encoder, thereby constructing an instance graph and applying a multi-channel kernel graph convolutional network to model it, further mining the high-dimensional features of the signal, and solving the problem of difficulty in extracting deep features in the source domain and the target domain under cross-working conditions.

2、本发明针对跨工况下滚动轴承振动信号源域和目标域间数据差异大的问题,采用了基于数据结构差异的损失函数和基于对抗的域间对其损失函数来联合缩小域间差异,同时,分类器则使用提取的域不变特征来完成跨域故障识别。2. To address the problem of large data differences between the source domain and target domain of rolling bearing vibration signals under different working conditions, the present invention adopts a loss function based on data structure differences and a loss function based on adversarial domain-to-domain comparison to jointly reduce the inter-domain differences. At the same time, the classifier uses the extracted domain-invariant features to complete cross-domain fault identification.

3、本发明解决了工业场景中标记数据难以获得、不同工况下数据分布差异较大、导致故障识别准确率低和泛化性能不足的问题,本发明不仅可以提取深层高维特征用以缩小域间差异,而且可以获得较优的诊断性能。3. The present invention solves the problems that labeled data is difficult to obtain in industrial scenarios, data distribution varies greatly under different working conditions, resulting in low fault identification accuracy and insufficient generalization performance. The present invention can not only extract deep high-dimensional features to narrow the differences between domains, but also obtain better diagnostic performance.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

图1为本发明实施例1提供的跨工况下滚动轴承故障靶向迁移诊断方法的流程图;FIG1 is a flow chart of a rolling bearing fault targeted migration diagnosis method under cross-operating conditions provided by Embodiment 1 of the present invention;

图2为本发明实施例1提供的基于样本数据特征的KNN图构建流程;FIG2 is a KNN graph construction process based on sample data features provided in Example 1 of the present invention;

图3为本发明实施例1提供的HFZZ旋转机械故障模拟平台示意图;FIG3 is a schematic diagram of the HFZZ rotating machinery fault simulation platform provided in Example 1 of the present invention;

图4为本发明实施例1提供的不同跨域诊断任务下的实验结果。FIG. 4 shows the experimental results of different cross-domain diagnosis tasks provided by Example 1 of the present invention.

具体实施方式DETAILED DESCRIPTION

下面结合附图与实施例对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

应该指出,以下详细说明都是例示性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are all illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and/or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and/or combinations thereof.

在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。In the absence of conflict, the embodiments of the present invention and the features of the embodiments may be combined with each other.

实施例1:Embodiment 1:

如图1所示,本发明实施例1提供了一种跨工况的滚动轴承故障靶向迁移诊断方法,包括以下过程:As shown in FIG1 , Embodiment 1 of the present invention provides a rolling bearing fault targeted migration diagnosis method across working conditions, including the following process:

S1:信号采集和数据集划分S1: Signal acquisition and data set division

通过机械故障测试平台收集不同工况下的数据,包含带标签的源域数据集

Figure BDA0004006128250000071
和不带标签的目标域数据集
Figure BDA0004006128250000072
并且样本数据被进一步划分为训练数据和测试数据;Data under different working conditions are collected through the mechanical failure test platform, including labeled source domain datasets
Figure BDA0004006128250000071
and unlabeled target domain dataset
Figure BDA0004006128250000072
And the sample data is further divided into training data and test data;

S2:网络模型构建和高维特征提取S2: Network model construction and high-dimensional feature extraction

通过特征编码器、图构建层、多通道核图卷积网络构建模型,提取样本数据的深层高维特征;The model is constructed through feature encoder, graph construction layer, and multi-channel kernel graph convolutional network to extract deep high-dimensional features of sample data;

S3:总体目标函数和模型优化S3: Overall objective function and model optimization

通过分类损失函数,缩小分类误差,用迁移学习方法,基于数据结构差异的损失函数和基于对抗的域间对其损失函数来联合缩小域间差异,并通过反向传播算法来优化总体目标函数中的参数;Through the classification loss function, the classification error is reduced, and the loss function based on data structure difference and the loss function based on adversarial domain alignment are used to jointly reduce the domain difference, and the back propagation algorithm is used to optimize the parameters in the overall objective function;

S4:模型测试和诊断结果输出S4: Model testing and diagnostic results output

测试阶段,目标域测试数据可以采用S2中优化的网络进行深层特征提取,分类器可以直接用于故障分类。During the testing phase, the target domain test data can use the optimized network in S2 for deep feature extraction, and the classifier can be directly used for fault classification.

S2中,网络模型构建和高维特征提取的构建过程,包括:In S2, the construction process of network model construction and high-dimensional feature extraction includes:

(1)特征编码器(1) Feature Encoder

主要包含生成器G和分类器C,生成器G用于对输入数据进行编码以获得高维特征表示,分类器C将对源和目标任务进行最终分类;It mainly includes generator G and classifier C. Generator G is used to encode the input data to obtain high-dimensional feature representation, and classifier C will perform the final classification of source and target tasks;

本实施例中,构建了一个一维卷积神经网络CNN用于特征提取和故障分类,CNN的隐藏层对应于G用来进行非线性特征映射以获得深度特征表示,CNN输出层的Softmax用作分类器C来获得每个故障类型的概率输出;In this embodiment, a one-dimensional convolutional neural network CNN is constructed for feature extraction and fault classification. The hidden layer of CNN corresponds to G and is used to perform nonlinear feature mapping to obtain deep feature representation. The Softmax of the CNN output layer is used as a classifier C to obtain the probability output of each fault type.

CNN中的高维特性表示,可以表示为:The high-dimensional feature representation in CNN can be expressed as:

G(x)=GL(GL-1(…G2(G1(x,w1)))) (1)G(x)=G L (G L-1 (…G 2 (G 1 (x, w 1 )))) (1)

其中,x是输入,G是CNN的特征映射,wi是CNN中第i层的学习权重,L是CNN的层数。Where x is the input, G is the feature map of CNN, wi is the learned weight of the i-th layer in CNN, and L is the number of layers of CNN.

(2)图构建层(2) Graph Construction Layer

样本在经过特征编码器之后,可以利用KNN算法来生成图构建层,界定所有输出之间的邻接关系,反映样本间的局部特性,构建流程如图2所示。After the samples pass through the feature encoder, the KNN algorithm can be used to generate a graph construction layer to define the adjacency relationship between all outputs and reflect the local characteristics between samples. The construction process is shown in Figure 2.

边的构造可以由下式表示:The construction of the edge can be expressed as follows:

Aij=KNN(k,Liji),Aij∈A (2)A ij =KNN(k,L iji ),A ij ∈A (2)

其中,Ωi={Li1,Li2,…,Lin}表示节点hi与其他所有节点的距离集合,k是超参数。Wherein, Ω i ={L i1 ,L i2 ,…,L in } represents the distance set between node h i and all other nodes, and k is a hyperparameter.

图构建层用来获得邻接矩阵A,并从最小批输入矩阵中获得实例图,主要分为如下几步:The graph construction layer is used to obtain the adjacency matrix A and obtain the instance graph from the minimum batch input matrix. It is mainly divided into the following steps:

(A)利用(1)中的特征编码器网络,从样本数据中获得高维特征映射,即X=G(x);(A) Using the feature encoder network in (1), we obtain a high-dimensional feature map from the sample data, i.e., X = G(x);

(B)提取的高维特征映射,输入到线性层中,并经过Softmax后可表示为:

Figure BDA0004006128250000091
(B) The extracted high-dimensional feature map is input into the linear layer and can be expressed as follows after Softmax:
Figure BDA0004006128250000091

(C)通过线性层的特征与其转置之间进行矩阵相乘计算,获得邻接矩阵A,即

Figure BDA0004006128250000092
(C) By performing matrix multiplication between the features of the linear layer and its transpose, the adjacency matrix A is obtained, that is,
Figure BDA0004006128250000092

(D)通过KNN算法,构造边关系,即

Figure BDA0004006128250000093
因此,邻接矩阵的构建总结如下:(D) Through the KNN algorithm, the edge relationship is constructed, that is,
Figure BDA0004006128250000093
Therefore, the construction of the adjacency matrix can be summarized as follows:

Figure BDA0004006128250000094
Figure BDA0004006128250000094

其中,A是构造的邻接矩阵,

Figure BDA0004006128250000095
是高维特征经过线性层和Softmax后的输出,normlize(·)表示归一化函数,
Figure BDA0004006128250000096
是稀疏邻接矩阵,KNN(·)返回邻接矩阵A在行方向上的前k个最大值的索引。Among them, A is the constructed adjacency matrix,
Figure BDA0004006128250000095
is the output of the high-dimensional feature after the linear layer and Softmax, normlize(·) represents the normalization function,
Figure BDA0004006128250000096
is a sparse adjacency matrix, KNN(·) returns the indices of the first k largest values in the row direction of the adjacency matrix A.

(3)多通道核图卷积网络(3) Multi-channel kernel graph convolutional network

图卷积网络是一种可以由数据的几何形状和结构进行表示的网络,可以提供更多的信息表示,并且基于节点间的连接关系来进行网络的学习。一般可简化为G=(A,X),其中A是邻接矩阵,可以反映出节点之间的连接关系,X是节点的特征。L=IN-D-1/2AD-1/2是Laplacian矩阵,其中D可以由邻接矩阵获得,即

Figure BDA0004006128250000097
IN是单位矩阵,图卷积使用滤波器gθ=diag(θ)来平滑输入信号,可以表示为:Graph convolutional network is a network that can be represented by the geometric shape and structure of data, which can provide more information representation and learn the network based on the connection relationship between nodes. It can be generally simplified to G = (A, X), where A is the adjacency matrix, which can reflect the connection relationship between nodes, and X is the feature of the node. L = I N -D -1/2 AD -1/2 is the Laplacian matrix, where D can be obtained from the adjacency matrix, that is,
Figure BDA0004006128250000097
I N is the identity matrix. Graph convolution uses a filter g θ = diag(θ) to smooth the input signal, which can be expressed as:

gθ*G x=UgθUTx (4)g θ * G x=Ug θ U T x (4)

其中,θ是可学习的参数,*G是图卷积操作,U是Laplacian矩阵的特征向量,UTx代表信号在图上的傅里叶变换。Among them, θ is a learnable parameter, * G is the graph convolution operation, U is the eigenvector of the Laplacian matrix, and U T x represents the Fourier transform of the signal on the graph.

公式(4)中定义的图卷积操作并不是局部化的,并且具有较高的计算量。采用公式(5)将卷积核限制为多项式的展开:The graph convolution operation defined in formula (4) is not localized and has a high computational cost. Formula (5) is used to restrict the convolution kernel to a polynomial expansion:

Figure BDA0004006128250000101
Figure BDA0004006128250000101

其中,K是多项式的阶数,

Figure BDA0004006128250000102
λ表示Laplacian矩阵的特征值。Where K is the order of the polynomial,
Figure BDA0004006128250000102
λ represents the eigenvalue of the Laplacian matrix.

采用一种多核图卷积网络来获取更宽广感受野的特征表示,它的卷积操作可以定义为:A multi-core graph convolutional network is used to obtain feature representation with a wider receptive field. Its convolution operation can be defined as:

Figure BDA0004006128250000103
Figure BDA0004006128250000103

其中,

Figure BDA0004006128250000104
是可学习的参数,
Figure BDA0004006128250000105
代表第ki个卷积核的高维特征表示。in,
Figure BDA0004006128250000104
is a learnable parameter,
Figure BDA0004006128250000105
Represents the high-dimensional feature representation of the k i -th convolution kernel.

假设,网络中多通道卷积核的个数为ki,也即通道数,网络层数为L=(1,2,…,l),也即网络深度,使用具有不同感受野的多核图卷积网络在不同通道上分别进行图卷积,构造多通道核图卷积网络:Assume that the number of multi-channel convolution kernels in the network is k i , that is, the number of channels, and the number of network layers is L = (1, 2, …, l), that is, the network depth. Use a multi-core graph convolution network with different receptive fields to perform graph convolution on different channels and construct a multi-channel kernel graph convolution network:

Figure BDA0004006128250000106
Figure BDA0004006128250000106

Figure BDA0004006128250000107
Figure BDA0004006128250000107

其中,X代表输入,A代表邻接矩阵,

Figure BDA0004006128250000108
代表可训练权值矩阵,G代表多通道核图卷积操作,
Figure BDA0004006128250000109
代表第ki个通道在第L层的高维特征表示,[·]代表特征拼接,H表示经过多通道核图卷积网络后的输出特征。Among them, X represents the input, A represents the adjacency matrix,
Figure BDA0004006128250000108
represents the trainable weight matrix, G represents the multi-channel kernel graph convolution operation,
Figure BDA0004006128250000109
represents the high-dimensional feature representation of the k i -th channel at the L layer, [ ] represents feature concatenation, and H represents the output feature after passing through the multi-channel kernel graph convolutional network.

S3中,总体目标函数和模型优化的构建过程,包括:In S3, the overall objective function and model optimization construction process includes:

为了充分利用数据的特性和深度网络结构,来缩小源域和目标域数据映射后的特征差异,采用了三部分损失函数,分别是分类损失、结构差异损失、对抗损失来构成总体的目标函数。In order to make full use of the characteristics of the data and the deep network structure to narrow the feature differences between the source domain and the target domain after data mapping, a three-part loss function is adopted, namely classification loss, structural difference loss, and adversarial loss to constitute the overall objective function.

(1)分类损失(1) Classification loss

通过交叉熵损失来估计真实标签和预测标签之间的分类损失,交叉熵损失可以定义为:The classification loss between the true label and the predicted label is estimated by the cross entropy loss, which can be defined as:

Figure BDA0004006128250000111
Figure BDA0004006128250000111

其中,

Figure BDA0004006128250000112
表示标签分类器的预测结果,LC表示交叉熵损失函数。E表示数学期望值。in,
Figure BDA0004006128250000112
represents the prediction result of the label classifier, L C represents the cross entropy loss function, and E represents the mathematical expectation.

(2)结构差异损失(2) Structural difference loss

采用下式的度量方式来缩小源域和目标域之间的结构差异损失,定义为:The following metric is used to reduce the structural difference loss between the source domain and the target domain, which is defined as:

Figure BDA0004006128250000113
Figure BDA0004006128250000113

其中,

Figure BDA0004006128250000114
Figure BDA0004006128250000115
分别表示第i个源域样本和第j个目标域样本经过特征提取器(本发明中指特征编码器和多通道核图卷积网络)后的特征映射。φ表示非线性特征映射,Ωk表示嵌入提取的特征到再生核希尔伯特空间RKHS中的距离度量,本实施例中,采用m个核的凸组合ku来对映射进行有效地估计:in,
Figure BDA0004006128250000114
and
Figure BDA0004006128250000115
Respectively represent the feature mapping of the i-th source domain sample and the j-th target domain sample after passing through the feature extractor (the feature encoder and the multi-channel kernel graph convolution network in the present invention). φ represents a nonlinear feature mapping, Ω k represents the distance metric embedded in the extracted features into the reproducing kernel Hilbert space RKHS. In this embodiment, a convex combination of m kernels k u is used to effectively estimate the mapping:

Figure BDA0004006128250000116
Figure BDA0004006128250000116

其中,αu是不同核的加权参数,且

Figure BDA0004006128250000117
本发明中对其在[0,1]中进行了平均分配。where αu is the weighting parameter of different kernels, and
Figure BDA0004006128250000117
In the present invention, it is evenly distributed in [0,1].

(3)对抗损失(3) Fighting Losses

采用对抗训练的方式,来处理域斜变和域移位的问题,通过域判别器来判断提取的特征是来自源域还是目标域,并且训练特征提取器来欺骗判别器。当两者达到minmax博弈平衡时,可以获得域不变特性,采用公式(12)中的损失函数作为对抗损失:The adversarial training method is used to deal with the problems of domain skew and domain shift. The domain discriminator is used to determine whether the extracted features are from the source domain or the target domain, and the feature extractor is trained to deceive the discriminator. When the two reach the minmax game equilibrium, the domain invariance feature can be obtained, and the loss function in formula (12) is used as the adversarial loss:

Figure BDA0004006128250000121
Figure BDA0004006128250000121

其中,D(·)为经过判别器后的特征输出,取值介于0和1之间,可以区分来自源域还是目标域。Among them, D(·) is the feature output after the discriminator, and its value is between 0 and 1, which can distinguish whether it comes from the source domain or the target domain.

(4)总体目标函数(4) Overall objective function

总体目标函数可以定义为:The overall objective function can be defined as:

Figure BDA0004006128250000122
Figure BDA0004006128250000122

其中,τ和

Figure BDA0004006128250000123
是可调参数。Among them, τ and
Figure BDA0004006128250000123
It is an adjustable parameter.

(5)参数优化(5) Parameter Optimization

针对公式(13)中的总体目标函数,通过反向传播算法来优化各个部分的参数如下:For the overall objective function in formula (13), the parameters of each part are optimized by back propagation algorithm as follows:

Figure BDA0004006128250000124
Figure BDA0004006128250000124

Figure BDA0004006128250000125
Figure BDA0004006128250000125

Figure BDA0004006128250000126
Figure BDA0004006128250000126

其中,

Figure BDA0004006128250000127
代表偏微分算子,η代表学习率,θF代表特征提取器的参数,θC代表分类器的参数,θD代表判别器的参数。in,
Figure BDA0004006128250000127
represents the partial differential operator, η represents the learning rate, θ F represents the parameters of the feature extractor, θ C represents the parameters of the classifier, and θ D represents the parameters of the discriminator.

本实施例所提出的跨工况下滚动轴承故障靶向迁移诊断方法的详细网络结构如表1所示,其中C代表故障种类:The detailed network structure of the rolling bearing fault targeted migration diagnosis method under different working conditions proposed in this embodiment is shown in Table 1, where C represents the fault type:

表1:跨工况下滚动轴承故障靶向迁移诊断方法的详细网络结构表Table 1: Detailed network structure of rolling bearing fault targeted migration diagnosis method under cross-operating conditions

Figure BDA0004006128250000131
Figure BDA0004006128250000131

本实施例提供如下具体案例:This embodiment provides the following specific cases:

搭建HFZZ旋转机械故障模拟平台,如图3所示,该平台由电动机、控制系统、径向加载装置、加速度传感器等组成。原始数据通过加速度传感器采集,采集频率为12.8kHz。采用3种不同的转速,来模拟不同工况下的轴承运行状态,分别是工况A(1750rpm),工况B(2000rpm)和工况C(2250rmp)。轴承故障通过加工,共生成了包含正常(NM)、内圈故障(IR,IRU)、外圈故障(OR)、滚动体故障(BA)以及多种混合故障在内的9种健康状况,其中包含单故障和复合故障,其详细信息如表2所示。在实验中,在每个运行转速下,从每个健康状况中采集100组数据,每个数据包含1024个数据点。每种工况下,总共获得了900个样本(100*9种健康状况),则数据集一共包含2700个样本(900*3种不同转速),具体说明见表3,并按照6:4的比例划分训练集与测试集。The HFZZ rotating machinery fault simulation platform was built, as shown in Figure 3. The platform consists of a motor, a control system, a radial loading device, an accelerometer, etc. The original data was collected by an accelerometer with a collection frequency of 12.8kHz. Three different speeds were used to simulate the operating conditions of the bearings under different working conditions, namely, working condition A (1750rpm), working condition B (2000rpm) and working condition C (2250rmp). The bearing fault was processed to generate 9 health conditions including normal (NM), inner ring fault (IR, IRU), outer ring fault (OR), rolling element fault (BA) and multiple mixed faults, including single faults and compound faults. The detailed information is shown in Table 2. In the experiment, 100 sets of data were collected from each health condition at each operating speed, and each data contained 1024 data points. Under each working condition, a total of 900 samples (100*9 health conditions) were obtained, so the data set contains a total of 2700 samples (900*3 different speeds). The specific description is shown in Table 3, and the training set and test set are divided according to the ratio of 6:4.

表2:滚动轴承健康状况详细描述表Table 2: Detailed description of rolling bearing health status

Figure BDA0004006128250000141
Figure BDA0004006128250000141

表3:滚动轴承迁移任务详细描述表Table 3: Detailed description of rolling bearing migration tasks

Figure BDA0004006128250000142
Figure BDA0004006128250000142

为了检验所提出的跨工况下滚动轴承故障靶向迁移诊断方法的优越性,使用几种最先进的深度神经网络算法,例如CNN(Baseline),MKMMD,JAN,DANN,CDAN,对跨工况下跨域诊断任务进行比较。在故障诊断实验中进行了十次试验,取平均值作为实验结果如表4和图4所示。从结果可以看出,本实施例所提方法取得了优越的性能,总平均准确度方面有所提升。当源域和目标域差异变得更大,跨域任务变得更困难时,所提算法仍然可以获得性能上的提升。In order to verify the superiority of the proposed method for targeted migration diagnosis of rolling bearing faults under cross-working conditions, several state-of-the-art deep neural network algorithms, such as CNN (Baseline), MKMMD, JAN, DANN, and CDAN, were used to compare cross-domain diagnosis tasks under cross-working conditions. Ten tests were conducted in the fault diagnosis experiment, and the average values were taken as the experimental results as shown in Table 4 and Figure 4. It can be seen from the results that the method proposed in this embodiment has achieved superior performance and has improved in terms of overall average accuracy. When the difference between the source domain and the target domain becomes larger and the cross-domain task becomes more difficult, the proposed algorithm can still achieve performance improvement.

表4:滚动轴承数据集的实验结果表Table 4: Experimental results of rolling bearing dataset

Figure BDA0004006128250000143
Figure BDA0004006128250000143

Figure BDA0004006128250000151
Figure BDA0004006128250000151

实施例2:Embodiment 2:

本发明实施例2提供了一种跨工况的滚动轴承故障靶向迁移诊断系统,包括:Embodiment 2 of the present invention provides a rolling bearing fault targeted migration diagnosis system across working conditions, including:

数据获取模块,被配置为:获取滚动轴承的振动信号;The data acquisition module is configured to: acquire a vibration signal of the rolling bearing;

特征提取模块,被配置为:通过特征编码器、图构建层和多通道核图卷积网络提取振动信号中的深层高维特征;A feature extraction module, configured to: extract deep high-dimensional features from the vibration signal through a feature encoder, a graph construction layer, and a multi-channel kernel graph convolutional network;

故障诊断模块,被配置为:根据提取的深层高维特征和分类器,得到滚动轴承故障诊断结果;The fault diagnosis module is configured to: obtain a rolling bearing fault diagnosis result based on the extracted deep high-dimensional features and the classifier;

其中,根据源域训练集结合分类器得到的深层高维特征得到分类损失,根据源域训练集和目标域训练集得到的深层高维特征得到结构差异损失,根据源域训练集和目标域训练集得到的深层高维特征以及对抗网络,得到对抗损失;Among them, the classification loss is obtained according to the deep high-dimensional features obtained by combining the source domain training set with the classifier, the structural difference loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set, and the adversarial loss is obtained according to the deep high-dimensional features obtained by the source domain training set and the target domain training set and the adversarial network;

以分类损失、第一参数与结构差异损失的乘积、第二参数与对抗损失的乘积三者的加和为总体损失函数,根据总体目标函数,通过反向传播算法优化特征提取器、分类器和判别器的参数。The overall loss function is the sum of the classification loss, the product of the first parameter and the structural difference loss, and the product of the second parameter and the adversarial loss. According to the overall objective function, the parameters of the feature extractor, classifier, and discriminator are optimized through the back propagation algorithm.

所述系统的工作方法与实施例1提供的跨工况的滚动轴承故障靶向迁移诊断方法相同,这里不再赘述。The working method of the system is the same as the rolling bearing fault targeted migration diagnosis method across working conditions provided in Example 1, and will not be repeated here.

实施例3:Embodiment 3:

本发明实施例3提供了一种计算机可读存储介质,其上存储有程序,该程序被处理器执行时实现如本发明实施例1所述的跨工况的滚动轴承故障靶向迁移诊断方法中的步骤。Embodiment 3 of the present invention provides a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps in the method for targeted migration diagnosis of rolling bearing faults across operating conditions as described in Embodiment 1 of the present invention.

实施例4:Embodiment 4:

本发明实施例4提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如本发明实施例1所述跨工况的滚动轴承故障靶向迁移诊断方法中的步骤。Embodiment 4 of the present invention provides an electronic device, comprising a memory, a processor, and a program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for targeted migration diagnosis of rolling bearing faults across working conditions as described in Embodiment 1 of the present invention are implemented.

本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of hardware embodiments, software embodiments, or embodiments combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) containing computer-usable program codes.

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to the flowchart and/or block diagram of the method, device (system), and computer program product according to the embodiment of the present invention. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the process and/or box in the flowchart and/or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。A person skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims (10)

1. A cross-working-condition rolling bearing fault targeted migration diagnosis method is characterized by comprising the following steps:
acquiring a vibration signal of a rolling bearing;
extracting deep high-dimensional features in the vibration signals through a feature encoder, a graph construction layer and a multi-channel kernel graph convolution network;
obtaining a rolling bearing fault diagnosis result according to the extracted deep high-dimensional characteristics and the classifier;
the method comprises the steps of obtaining classification loss according to deep high-dimensional characteristics obtained by combining a classifier with a source domain training set, obtaining structural difference loss according to the deep high-dimensional characteristics obtained by combining the source domain training set and a target domain training set, and obtaining countermeasures loss according to the deep high-dimensional characteristics obtained by combining the source domain training set and the target domain training set and a countermeasures network;
and optimizing parameters of the feature extractor, the classifier and the discriminator through a back propagation algorithm according to the overall objective function by taking the sum of the classification loss, the product of the first parameter and the structural difference loss and the product of the second parameter and the counterattack loss as an overall loss function.
2. The method for diagnosing the fault-targeted migration of the rolling bearing under the cross-working condition of claim 1, wherein the method comprises the following steps of,
the graph construction layer is used for acquiring an adjacency matrix and comprises the following steps:
obtaining a high-dimensional feature map, i.e., x=g (X), from the sample data according to the feature encoder network;
the extracted high-dimensional characteristic map is input into a linear layer and expressed as after Softmax
Figure FDA0004006128240000011
By performing matrix multiplication calculation between the features of the linear layer and its transpose, an adjacency matrix a is obtained,
Figure FDA0004006128240000012
by KNN algorithm, edge relations are constructed, i.e
Figure FDA0004006128240000013
3. A cross-operating mode rolling bearing fault targeted migration diagnostic method as claimed in claim 2, wherein,
a multi-channel kernel graph rolling network comprising:
Figure FDA0004006128240000021
Figure FDA0004006128240000022
wherein X represents the input, A represents the adjacency matrix,
Figure FDA0004006128240000023
represents a trainable weight matrix, G represents a multi-channel kernel graph convolution operation, ++>
Figure FDA0004006128240000024
Represents the kth i High-dimensional characterization of individual channels at layer L, [. Cndot.]Representing feature stitching, H represents the output features after passing through the multi-channel kernel graph convolution network.
4. The method for diagnosing the fault-targeted migration of the rolling bearing under the cross-working condition of claim 1, wherein the method comprises the following steps of,
cross entropy loss L C Comprising:
Figure FDA0004006128240000025
wherein ,
Figure FDA0004006128240000026
representing the prediction result of the classifier, E representing the mathematical expectation,/->
Figure FDA0004006128240000027
For the source domain sample, ++>
Figure FDA0004006128240000028
For its tag.
5. The method for diagnosing the fault-targeted migration of the rolling bearing under the cross-working condition of claim 1, wherein the method comprises the following steps of,
structural differential loss L s Comprising:
Figure FDA0004006128240000029
wherein ,
Figure FDA00040061282400000210
and
Figure FDA00040061282400000211
Respectively representing the feature mapping of the ith source domain sample and the jth target domain sample after the ith source domain sample passes through a feature extractor, phi represents the nonlinear feature mapping, omega k Representing distance metrics of embedded extracted features into regenerated kernel Hilbert space RKHS, employing convex combinations k of m kernels u To effectively estimate the mapping:
Figure FDA00040061282400000212
wherein ,αu Is a weighting parameter of different kernels and +.>
Figure FDA00040061282400000213
E represents a mathematical expectation value, +.>
Figure FDA0004006128240000031
For the source domain sample, ++>
Figure FDA0004006128240000032
Is a target domain sample.
6. The method for diagnosing the fault-targeted migration of the rolling bearing under the cross-working condition of claim 1, wherein the method comprises the following steps of,
countering loss L AD Comprising:
Figure FDA0004006128240000033
wherein D (·) is the feature output after passing through the discriminator,
Figure FDA0004006128240000034
and
Figure FDA0004006128240000035
Respectively representing the feature mapping of the ith source domain sample and the jth target domain sample after the feature extractor, E represents the mathematical expectation value,>
Figure FDA0004006128240000036
for the source domain sample, ++>
Figure FDA0004006128240000037
Is a target domain sample.
7. The method for diagnosing the fault-targeted migration of the rolling bearing under the cross-working condition of claim 1, wherein the method comprises the following steps of,
optimizing parameters of the feature extractor, classifier and arbiter by a back propagation algorithm, comprising:
Figure FDA0004006128240000038
wherein ,
Figure FDA0004006128240000039
represents partial differential operator, eta represents learning rate, theta F Parameters representative of feature extractor, θ C Representing the parameters of the classifier, θ D Representing parameters of the arbiter, L C For cross entropy loss, L s For structural differential loss, L AD To combat losses.
8. A cross-condition rolling bearing fault targeted migration diagnostic system, comprising:
a data acquisition module configured to: acquiring a vibration signal of a rolling bearing;
a feature extraction module configured to: extracting deep high-dimensional features in the vibration signals through a feature encoder, a graph construction layer and a multi-channel kernel graph convolution network;
a fault diagnosis module configured to: obtaining a rolling bearing fault diagnosis result according to the extracted deep high-dimensional characteristics and the classifier;
the method comprises the steps of obtaining classification loss according to deep high-dimensional characteristics obtained by combining a classifier with a source domain training set, obtaining structural difference loss according to the deep high-dimensional characteristics obtained by combining the source domain training set and a target domain training set, and obtaining countermeasures loss according to the deep high-dimensional characteristics obtained by combining the source domain training set and the target domain training set and a countermeasures network;
and optimizing parameters of the feature extractor, the classifier and the discriminator through a back propagation algorithm according to the overall objective function by taking the sum of the classification loss, the product of the first parameter and the structural difference loss and the product of the second parameter and the counterattack loss as an overall loss function.
9. A computer readable storage medium having stored thereon a program, which when executed by a processor, implements the steps of the cross-regime rolling bearing fault targeted migration diagnostic method according to any one of claims 1 to 7.
10. An electronic device comprising a memory, a processor and a program stored on the memory and executable on the processor, wherein the processor performs the steps in the cross-regime rolling bearing fault targeted migration diagnostic method of any one of claims 1-7 when the program is executed.
CN202211631915.3A 2022-12-19 2022-12-19 Cross-working-condition rolling bearing fault targeted migration diagnosis method and system Pending CN116026593A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211631915.3A CN116026593A (en) 2022-12-19 2022-12-19 Cross-working-condition rolling bearing fault targeted migration diagnosis method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211631915.3A CN116026593A (en) 2022-12-19 2022-12-19 Cross-working-condition rolling bearing fault targeted migration diagnosis method and system

Publications (1)

Publication Number Publication Date
CN116026593A true CN116026593A (en) 2023-04-28

Family

ID=86073278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211631915.3A Pending CN116026593A (en) 2022-12-19 2022-12-19 Cross-working-condition rolling bearing fault targeted migration diagnosis method and system

Country Status (1)

Country Link
CN (1) CN116026593A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116952554A (en) * 2023-07-05 2023-10-27 北京科技大学 Multi-sensor mechanical equipment fault diagnosis method and device based on graph rolling network
CN116977708A (en) * 2023-06-14 2023-10-31 北京建筑大学 Bearing intelligent diagnosis method and system based on self-adaptive aggregation visual view
CN117194983A (en) * 2023-09-08 2023-12-08 北京理工大学 A bearing fault diagnosis method based on progressive conditional domain adversarial network
CN117232577A (en) * 2023-09-18 2023-12-15 杭州奥克光电设备有限公司 Optical cable distributing box bearing interior monitoring method and system and optical cable distributing box

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116977708A (en) * 2023-06-14 2023-10-31 北京建筑大学 Bearing intelligent diagnosis method and system based on self-adaptive aggregation visual view
CN116977708B (en) * 2023-06-14 2024-04-12 北京建筑大学 Bearing intelligent diagnosis method and system based on self-adaptive aggregation visual view
CN116952554A (en) * 2023-07-05 2023-10-27 北京科技大学 Multi-sensor mechanical equipment fault diagnosis method and device based on graph rolling network
CN117194983A (en) * 2023-09-08 2023-12-08 北京理工大学 A bearing fault diagnosis method based on progressive conditional domain adversarial network
CN117194983B (en) * 2023-09-08 2024-04-19 北京理工大学 A bearing fault diagnosis method based on progressive conditional domain adversarial network
CN117232577A (en) * 2023-09-18 2023-12-15 杭州奥克光电设备有限公司 Optical cable distributing box bearing interior monitoring method and system and optical cable distributing box
CN117232577B (en) * 2023-09-18 2024-04-05 杭州奥克光电设备有限公司 Optical cable distributing box bearing interior monitoring method and system and optical cable distributing box

Similar Documents

Publication Publication Date Title
WO2024021246A1 (en) Cross-device incremental bearing fault diagnosis method based on continuous learning
Zhao et al. Deep multi-scale convolutional transfer learning network: A novel method for intelligent fault diagnosis of rolling bearings under variable working conditions and domains
Chen et al. Intelligent fault diagnosis for rotary machinery using transferable convolutional neural network
Zhang et al. Limited data rolling bearing fault diagnosis with few-shot learning
CN116026593A (en) Cross-working-condition rolling bearing fault targeted migration diagnosis method and system
US11042802B2 (en) System and method for hierarchically building predictive analytic models on a dataset
Li et al. Multiscale dynamic fusion prototypical cluster network for fault diagnosis of planetary gearbox under few labeled samples
CN112926641A (en) Three-stage feature fusion rotating machine fault diagnosis method based on multi-modal data
CN112101085A (en) Adaptive intelligent fault diagnosis method based on importance weighted domain impedance
CN113963165B (en) Small sample image classification method and system based on self-supervision learning
Graening et al. Shape mining: A holistic data mining approach for engineering design
CN117386344B (en) Drilling abnormal condition diagnosis method and system based on two-stage learning
Lv et al. Deep transfer network with multi-kernel dynamic distribution adaptation for cross-machine fault diagnosis
Liu et al. Structural discrepancy and domain adversarial fusion network for cross-domain fault diagnosis
Li et al. Unsupervised continual source-free network for fault diagnosis of machines under multiple diagnostic domains
CN114528547A (en) ICPS (information storage and protection System) unsupervised online attack detection method and device based on community feature selection
CN117272195A (en) Block chain abnormal node detection method and system based on graph convolution attention network
Wang et al. Bearing fault diagnosis algorithm based on granular computing
CN114757271B (en) Social network node classification method and system based on multichannel graph convolutional network
CN118410306B (en) Rolling bearing fault cross-equipment depth domain self-adaptive migration diagnosis method
CN117235490B (en) Fault self-adaptive diagnosis method integrating deep volume and self-attention network
Nemer et al. A review of research on industrial time series classification for machinery based on deep learning
CN109255433B (en) A Similarity-Based Approach to Community Detection
CN115019342B (en) Endangered animal target detection method based on class relation reasoning
Jiang et al. Fault Diagnosis of Gearbox in Multiple Conditions Based on Fine‐Grained Classification CNN Algorithm

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