数字孪生 是当前故障诊断领域的研究热点之一。它们只需要很少的实际信号就可以产生故障样本,从而降低了实验成本。小编整理搜集了一些相关研究与大家进行分享。本期分享的是不平衡样本下基于逆物理信息的数字孪生轴承故障诊断神经网络 ,该论文发表于Knowledge-Based Systems ,适合具备一定有信号基础并打算在故障诊断领域深入挖掘的学习者,非常值得阅读 ! 论文链接 : 通过点击最左下角的阅读 原文 进行在线阅读及下载。 1 论文基本信息 论文题目 :Inverse physics–informed neural networks for digital twin–based bearing fault diagnosis under imbalanced samples
论文期刊 :Knowledge-Based Systems
Doi : https://doi.org/10.1016/j.knosys.2024.111641
论文时间 :2024年
作者 :Yi Qin (a), Hongyu Liu (a), Yi Wang (a), Yongfang Mao(b)
机构:
a: State Key Laboratory of Mechanical Transmission for Advanced equipment, Chongqing University, Chongqing 400044, PR China; b: School of Automation, Chongqing University, Chongqing 400044, PR China. 2 目录 注:本文只选中原论文部分进行分享,若想拜读,请下载原论文进行细读。 3 摘要 在实际工程中,各种故障类别的轴承样本不足,给滚动轴承的智能故障诊断带来了很大障碍。为了解决样本不平衡的问题,本工作探索了一种基于数字孪生技术的轴承故障数据生成方法。首先,通过将轴承动态模型嵌入到神经网络中,建立逆物理信息神经网络(Inverse Physics–informed Neural Network, PINN)来识别动态模型参数。在该网络中,设计了一个边界损失来快速确定加速网络收敛的参数的近似范围,并构造了一个真值损失来评估模拟数据和实际数据之间的频谱差异。然后,利用逆PINN、轴承故障动态模型和实际振动数据,提出了一种基于数字孪生的故障数据生成方法,用于在多种工况和故障模式下生成高质量的轴承故障样本。最后,将该方法应用于特定工况下的轴承故障振动样本生成。将生成样本用于训练诊断网络,解决了样本不平衡的问题。多个实验的对比结果表明,所提出的数据生成方法有效地提高了跨工况轴承故障诊断的精度,并超过了多种先进的诊断方法。
关键词 : 数字孪生,物理信息神经网络,轴承动态模型,故障诊断,样本不平衡
4 引言 已提出的轴承故障诊断方法大致可分为基于信号处理的方法[1]和基于人工智能的技术[2]。深度学习方法,包括卷积神经网络[3]、深度神经网络[4]和联邦学习[5],在智能轴承故障诊断中占据主导地位。然而,这些深度学习方法通常需要大量的故障样本。在许多工程应用中,正常的轴承振动信号样本非常丰富,而故障振动信号样本的可用性有限。此外,一些故障的样本可能比其他故障少得多。因此,轴承数据集表现出样本不平衡,这可能会显著影响故障诊断的效率。故障样本的稀缺性和数据集的不平衡性可能会阻碍故障诊断模型识别某些故障类型。为了解决这些问题,Zhao等人[6]提出了一种基于小波变换和CNN的机械故障诊断方法。[7]研究了一种能够进行多领域特征提取和选择的代价敏感学习方法,有效地降低了漏报率和平均代价。此外,Li 等人[8]在训练数据类型较少的情况下,引入持续学习来提高故障诊断的效果。另一种可用的方法是添加相关的故障样本。 基于实际信号的数据生成是一种可行的采样增强技术。Wang等人[9]建立了一种自适应、高效的亚像素网络DAESPN用于变速条件下的故障样本扩充和智能故障诊断。An等人[10]开发了一种主动构造的数据增强模型以增强轴承故障诊断。Pei等人[11]设计了一种改进的基于元学习的Wasserstein自动编码器以增加滚动轴承数据,提高轴承故障诊断能力。然而,这些数据驱动的方法往往需要大量的实际样本来训练生成式模型,而且生成的样本往往与训练样本相似,限制了多样性,阻碍了故障诊断准确率的提高。 另一种方法是利用故障动态模型来生成故障样本。Shao等人[12]建立了转子-轴承系统动力学模型,该模型产生了不同健康状态的轴承振动信号;这些信号随后被用于具有生成式对抗性网络(Generative Adversarial Network, GAN)的实际轴承的故障诊断。Qin等人提出了一种利用时变激励生成故障样本的故障动力学模型[13]。类似地,Liu 等人的研究[14]建立了用于故障样本产生的现实的轴承振动模型,并将其用于基于GAN的轴承故障诊断。Liu等人[15]提出了一种个性化故障诊断方法,通过有限元仿真训练智能传感器网络,提高了故障样本缺失时故障诊断的有效性。然而,基于动态模型的数据生成方法需要确定故障动态模型的关键参数,如刚度和阻尼比。这些参数通常是由经验设定的。估计参数和实际参数之间的差异导致了仿真数据和测量数据之间的差异,影响了仿真数据对诊断模型训练的效果。数字孪生模型只需要很少的实际信号就可以生成故障样本,从而节省了实验材料和时间。 基于数字孪生的样本扩充技术结合数据驱动和基于模型的方法,并引起了人们的极大兴趣。它们只需要很少的实际信号就可以产生故障样本,从而降低了实验成本。提出了几种提高轴承振动数据的方法。Piltan等人[16]建立了机械模型与自适应观测器相结合的轴承数字孪生框架,实现了轴承裂纹尺寸的精确检测。为了阐述退化过程中的轴承振动响应,Qin等人[17]建立了全生命周期的轴承数字孪生模型。值得注意的是,一种新的基于循环GAN的映射网络(CycleGAN)[18]最大限度地减少了模拟信号和实际信号之间的差异信号进行调整的。此外,动态模型没有完全集成到神经网络中。 为了将物理模型与神经网络相结合,人们提出了各种PINN,其中一些已成功地应用于滚动轴承的故障预测和诊断。例如,Chen等人[19]提出了一种物理信息-退化一致性递归神经网络来预测机器RUL,并将机器部件的退化能力集成到神经网络中。Yucesan等人[20]建立了基于物理信息的PINN模型,用于风电机组主轴承疲劳建模。然而,这些方法主要是利用物理信息进行RUL预测,而不是求解物理模型。已经提出了一些用于求解各种物理模型的PINN。例如,Ling等人[21]用PINN和Mao等人预测湍流和高速气动扰流。[22]用PINN对高速气动流动进行了模拟。不幸的是,现有的PINN无法计算二自由度轴承动力学模型的振动响应,因为两个微分方程组是相互耦合的。因此,需要对PINN进行改进,以便对故障滚动轴承进行动态建模。 刚度和阻尼比是轴承动力学模型中的两个重要参数,一般都是凭经验设定的,往往会导致仿真信号与实测信号存在较大的偏差。为了解决这一问题,提出了一种滚动轴承的逆PINN。与传统的PINN方法不同的是,直接利用滚动轴承的二自由度微分方程组而不是损失函数来计算振动响应,并将网络损耗视为模拟信号与采集的真实信号之间的频谱差异,以防止沿X和Y方向的动态微分方程组之间的耦合效应。同时,逆PINN 的优化对象是轴承的动态参数,而不是振动响应,因此可以精确地辨识轴承动力学模型的参数。利用逆PINN和轴承动力学模型,提出了一种新的故障数据生成方法,用于产生各种工况下的高精度轴承振动信号。最后,将生成的数据成功地应用于样本不平衡情况下的跨工况轴承故障诊断。本研究的贡献如下: 1. 为了辨识动态模型的参数,我们提出了一种逆PINN方法,将轴承动态模型集成到神经网络中。在该网络中,边界损失被设计为快速确定参数的近似范围,从而加速网络收敛,并构造了真值损失来评估模拟数据与实际数据之间的频谱差异。利用所提出的逆PINN方法,可以准确地辨识轴承的动态参数。 2. 利用逆PINN、轴承故障动态模型和采集的振动信号,提出了一种新的基于数字孪生的数据生成方法,可以在多个工作条件下生成高质量的故障样本。该方法从单一工况下采集的少量轴承样本中生成出不同工况下高质量的轴承故障振动样本。由此可见,该方法实现了故障滚动轴承的数字孪生。 3. 对于具有不平衡轴承振动样本的跨工况故障诊断,该方法被用来生成特定工况下的轴承故障振动样本。然后,将样本用于诊断网络的训练。实验结果表明,与现有的各种SOTA数据生成和不平衡故障诊断方法相比,该方法能够提高跨工况轴承故障诊断的精度。 5 准备工作 5.1 PINN
PINN可以结合经验数据和抽象的数学运算符[23],例如偏微分方程(Partial Differential Equations, PDE)。一旦将物理原理或约束嵌入到PINN中,网络就会表现出更强的可解释性和抗干扰性。一般来说,一个典型的Pinn的损失被定义为: 其中, 表示输出数据和实际数据之间的差距, 表示初始值和边界损失, 由微分方程式决定, , 和 用于平衡三个损失之间相互作用的影响的权重,可以由使用者定义或自动调整。PINN的参数通过各种基于梯度的优化器(例如ADAM和L-BFGS)来最小化损失来进行优化。然而,传统的PINN算法主要用于求解单个偏微分方程组。 5.2 数字孪生
作为一种新的设备运行和维护方法,数字孪生技术起源于2003年格里夫斯教授在密歇根大学的一次演讲[24]。在本讲座中,定义了数字孪生模型,它包含三个基本组成部分:虚拟空间、物理空间和数据流的互联性。2012年,NASA[25]将数字孪生视为集成多个物理域、尺度和概率的模拟过程。这一创新方法在各种设备的设计、制造和运行维护方面已被证明是成功的。Deon等人[26]建立了火电厂发电机组预知检修的数字孪生模型。如图1所示,数字孪生框架主要由物理空间、虚拟空间以及它们之间的数据流或映射关系组成。物理空间表示被分析的系统或部件;虚拟空间是所建立的分析模型,诸如有限元模型或集总参数动态模型。并且数据流或映射用于桥接物理空间和虚拟空间之间的关系。针对轴承故障诊断中样本不均衡的问题,利用数字孪生技术生成轴承故障样本。
6 提出的轴承故障振动数据生成方法 6.1 轴承故障动力学模型
首先建立了轴承故障动力学模型。通过将外圈上的轴承故障简化为立方体,降低了模型的复杂性(图2)。立方体的长度、宽度和高度分别由L、B和H表示。根据故障和滚动轴承参数,通过下式计算故障引起的最大位移激励 : 式中,D表示滚子的直径, 表示位移激励的理论最大值。 考虑半正弦激励,我们可以定义位移激励函数如下[27]: 其中 表示进入故障区域时滚动体与X轴的夹角, 表示离开故障区域时滚动体与X轴的夹角, 表示第j个滚动体的位置。 滚动体与滚子接触的法向变形包含了轴承故障位移激励函数的影响。接触变形表达式为: 其中, 和 分别表示滚动轴承在X和Y方向上的位移,ε表示轴承公差。我们建立了故障轴承的两自由度动力学方程: 其中m表示当量轴承质量;K表示力变形系数(即刚度);c表示当量阻尼比; 表示第j滚子在接触区的判断系数,下式给出: 刚度K和阻尼比c是轴承故障动力学模型中的两个关键参数,对模拟振动响应的质量有很大影响。 然而,它们通常是通过实际经验来确定的。 在下一节中,我们将讨论如何通过测量的轴承故障振动信号来准确识别K和c。
在这一部分中,将详细阐述所提出的逆PINN。
逆PINN原理 :与现有的PINN不同,本文提出的逆PINN的优化目标是获得两个动力学参数,而不是一个轴承振动响应。逆PINN的结构图如图3所示,逆PINN的输入是从实际物理空间测得的滚动轴承故障振动信号,其输出是动态参数。在网络训练期间,动态参数被更新。根据滚动轴承的最新动态参数和已知的滚子数、滚子直径、故障位置、初始角位置等结构参数,采用Runge-Kutta方法结合张量计算,求解了非线性轴承故障动力学方程,得到了故障轴承在X、Y方向的振动响应。然后,将测量信号和模拟信号的FFT之差计算为损失。此外,还设计了两个动态参数的边界损失,用于网络训练。
由于刚度K和阻尼比c之间的差高一个数量级,所以K不是由神经网络直接输出的,而是通过将网络输出值乘以七个数量级而得到的。这种方法可以确保网络输出的两个参数尽可能地处于相似的数量级,并有利于所提出的网络的PINN收敛。逆PINN的目标与原PINN不同,前者的结构比后者更简单。
逆PINN的损失函数 :逆PINN的损失函数由两部分组成。第一部分是边界损失,它包含所有未知参数的边界条件。传统PINN中的边界损失主要是用来约束某一时间点或空间点的响应值或物理量。例如,在基于流体力学的PINN中,某一位置的流体速度为零或某一值,可以利用这些信息来表示边界损失。然而,在逆PINN中,目标值不再是特定的响应,而是未知的物理参数。这些参数的近似值范围被用作边界条件。当动态参数在特定范围内更新时,可以加速网络训练。边界损失表示 为: 其中, 表示由网络模型输出的第j个参数; 和 分别是参数的下限和上限;以及n是参数的数目。本文考虑了两个动力参数:轴承刚度K和阻尼比c,并根据这两个参数的取值范围设定了具体的上下限。根据目标方位的经验知识,可以确定参数值的近似范围。例如,轴承的阻尼比通常在10到1000之间,而刚度通常在 到 之间。
损失函数的第二部分是模拟振动与测量信号之间的偏差。通过在频域中计算这部分的损失函数,抑制了被测信号样本中噪声的影响。通过将时域信号转换为频域信号来防止由相位差和其他因素造成的干扰。此外,动态参数,包括刚度和阻尼比,对频谱的影响比时域信号更显著。因此,基于频域信号的损失函数保证了这些参数的准确、快速的优化。将测量的振动信号和模拟的振动信号的频域变换结果之间的均方根误差(RMSE)作为损失,称为 。这一损失被定义为:
其中,X和Y是测量的振动信号沿X和Y方向的傅里叶变换;X和Y是对应于通过逆PINN更新的动态参数的模拟振动信号的傅里叶变换;real(·)和imag(·)分别表示数据的实部和虚部。真值损失有效地减小了模拟信号和测量信号之间的频谱差异。利用这一损失来训练逆PINN,可以得到更接近实际情况的动态参数。然后,将这些动态参数引入到轴承故障动态模型中,从而提高了仿真信号与实测信号之间的相似性。 总而言之,建议的逆PINN的损失公式如下: 6.3 轴承故障数据生成 在此基础上,提出了一种基于数字孪生模型的故障数据生成方法,并给出了具体的实现步骤。 第一步:输入滚动轴承的结构参数,将在某一工况下采集到的测量振动信号送入逆变器。根据专家经验设置两个动态参数的取值范围,然后分别在两个范围内初始化两个动态参数。 第二步:根据故障轴承的结构参数、工况参数和动力学参数,采用Runge-Kutta算法结合张量计算,求解故障轴承的二自由度动力学方程。然后,得到模拟的振动信号。 第三步:计算模拟振动信号和测量信号的傅里叶变换。用逆PINN法计算输出动态参数。 第四步:通过公式计算 。然后,用EQ计算边界损失 。 第五步:使用获得的损失,用ADAM优化器训练反PINN。通过损耗调整网络参数。 第六步:重复步骤2到5,直到逆PINN收敛。完成网络训练后,准确识别两个动态参数K和c。 第七步:将获得的动态参数提交到轴承故障动态模型中,计算出各种工况和故障模式(包括故障大小)下的模拟振动信号。 第八步:在轴承振动响应中加入高斯白噪声,以提高生成数据的多样性。 7 实验结果与分析 7.1 实验描述
提出的方法将通过两个实验进行验证。第一个实验基于Case Western Reserve University(CWRU)的轴承故障测试。两个测试SKF轴承(6205-2RS)位于轴承处。采用驱动端故障直径为0.54 mm(即L=B=H=0.54 mm)的故障滚动轴承的振动信号。在实验过程中,电机转速为1750rpm,信号采集频率为12 kHz。第二个实验是基于西安交通大学的轴承故障测试[28]。主要用于轴承退化试验。采用两台由PCB公司生产的加速度计,以25.6 kHz的采样频率沿X和Y方向测量轴承加速度信号。
表4 由逆PINN和专家经验在两个实验中获得的动力学参数
图7 CWRU实验中模拟信号与实测信号的比较:
(A)由优化参数得到的模拟结果;(B)由经验参数得到的模拟结果
在CWRU和XJTU实验中,故障滚动轴承的动态参数(表4)是通过逆PINN估计的。同时,表4也列出了其他文献[17,27]设定的经验参数。观察到这两类动力学参数之间有很大的差异。不同的动力参数会导致不同的振动响应。在CWRU实验中,利用这两种动力学参数计算了振动信号。模拟信号和测量信号同时显示在图7中。在时间域波形中,很难确定K优化的参数是否优于经验参数。
图8 CWRU实验中仿真频谱与实验频谱的比较:
(A)基于经验参数的仿真结果;(B)基于逆PINN优化参数的仿真结果
图9 XJTU实验中仿真频谱与实验频谱的比较:
(A)基于经验参数的仿真结果;(B)基于逆PINN优化参数的仿真结果
动态参数K和c不影响旋转频率和故障特征频率,但影响其他频率成分,如固有频率。为了弄清动力参数对振动响应的影响,我们计算了由两种动力参数得到的振动信号的频谱,并分别将它们与被测信号的频谱进行了比较。CWRU实验的比较结果如图8(A)和(B)所示。提出的逆PINN方法得到的频谱与实测信号的频谱接近,而基于经验参数的频谱的主成分与实际信号的主成分有很大的偏差。XJTU实验中的比较结果如图9(A)和(B)所示。也得到了类似的结论。
根据振动信号的频谱和文章中的时频频谱,根据最大峰值可以估计出轴承的固有频率。为了说明该参数优化方法相对于其他方法的优越性,计算了逆PINN法和专家经验估计的固有频率的误差和错误率。由所提出的逆PINN估计的固有频率具有很高的精度,这意味着所提出的逆PINN可以求解动态参数,并且提高了仿真模型的精度。
7.3 数据生成示例
图11以外圈故障为样本,在40赫兹的转速和10kN的载荷下,
模拟和测量了轴承振动信号的频谱。
在用逆PINN识别出轴承的动态参数后,利用所开发的基于数字孪生故障数据生成方法,可以模拟故障轴承在任何工况下的振动响应。以XJTU实验中的轴承故障数据生成为例。利用所获得的动力学参数和轴承故障动力学模型,生成了在40 Hz旋转频率和10kN载荷下发生外环故障的轴承振动数据。计算了产生的和测量的信号的频谱(图11)。这两个频谱具有相似的特征谱线,即有效地模拟和产生了轴承振动信号中的主要特征频率分量。
8 总结 本文研究了一种轴承数据生成方法,以提高对不平衡样本的跨工况轴承故障诊断能力。为了准确识别轴承的动态参数,将轴承动力学模型与神经网络相结合,建立了一种新的逆PINN模型。这种网络的原理与传统的PINN相反。在逆PINN中,建立边界损失和真值损失来加速网络收敛,并将其用于评估模拟频域数据与真实频域数据之间的差异。提出了一种基于逆PINN、故障轴承动力学模型和实际数据的故障数据生成方法,用于生成特定工况下补充样本的轴承故障信号。然后,用它们来训练故障诊断网络。实验结果表明,该方法能有效克服轴承故障诊断中的样本平衡问题,优于SOTA数据生成方法和不平衡故障诊断方法。未来,我们将从强背景噪声的测量信号中识别模型参数,进一步提高生成数据的质量和多样性。
编辑:曹希铭
校核:李正平 、陈凯歌、赵栓栓、董浩杰、 赵学功、 白亮、陈少华
该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除
点击左下角阅读原文 ,即可在线阅读论文。