首页/文章/ 详情

中科院一区Top论文学习|基于两阶段更新数字孪生的轴承剩余寿命预测

10月前浏览972

本期给大家推荐基于两阶段更新数字孪生的轴承剩余寿命预测在工业生产中,轴承作为核心部件直接影响系统运行稳定性,精准预测其剩余寿命(Remaining useful life, RUL)是保障生产的关键。当前主流RUL预测方法普遍忽略轴承实时运行健康状态,导致预测结果偏差较大。针对这一问题,本文提出融合两阶段更新数字孪生与双关联动态图卷积网络的轴承RUL预测方案:构建具有外圈缺陷扩展特性的轴承缺陷演化模型,建立孪生模型与真实轴承的实时交互机制,并校准缺陷优化全生命周期缺陷曲线,借助双向长短期记忆网络关联真实轴承振动特征与孪生缺陷,完成实时映射;最后将映射缺陷纳入RUL预测特征空间,提取物理与数字空间特征关联,实现最终预测。

论文链接:通过点击最左下角的阅读原文进行在线阅读及下载

论文基本信息

论文题目:Prediction of bearing remaining useful life based on a two-stage updated digital twin  

论文期刊:Advanced Engineering Informatics

Doihttps://doi.org/10.1016/j.aei.2025.103123

作者:Deqiang He (a), Jiayang Zhao (a), Zhenzhen Jin* (a,d), Chenggeng Huang (b), Fan Zhang (c), Jinxin Wu (a)
论文时间: 2025年
机构: 

a: Guangxi Key Laboratory of Manufacturing System & Advanced Manufacturing Technology, School of Mechanical Engineering of Guangxi University, Nanning 530004,  China 

b: School of Automation Engineering, University of Electronic Science and Technology, China 

c: Southwest Jiaotong University School of Design, China 

d: Guangxi Key Laboratory of Precision Navigation Technology and Application, Guilin University of Electronic Technology, Guilin 541004, China

作者简介:靳震震近五年围绕载运装备故障诊断与智能运维开展研究,主持和参与国家自然科学基金联合基金重点项目、面上项目、广西科技计划项目等项目10余项,发表高水平学术论文50余篇,其中第一作者或通讯作者30篇,ESI热点论文3篇,ESI高被引论文7篇,《光明日报》刊发文章1篇,申请发明专利24件,谷歌学术引用次数共计超1600次,H指数20,I10指数29。参编地方标准1部。获2024年广西大学优秀博士学位论文、2025年广西科技进步一等奖(排名第八)、2025年中国物流与采购联合会科学技术奖一等奖(排名第三)、2024年中国机械行业产教融合教育教学创新大赛全国二等奖(排名第六)、2025年广西大学本科教学成果一等奖(排名第三)、2025年广西大学研究生教学成果特等奖(排名第四)、2025年中国机械行业产教融合教育教学创新大赛全国一等奖(排名第四)各1项,入选全球前2%顶尖科学家榜单、广西青年科技人才托举工程。(来源于广西大学教师信息网)
通讯作者邮箱:sdkjdxjz@163.com

摘要

轴承作为工业生产中的关键部件,对系统的平稳运行至关重要。准确预测轴承剩余寿命(Remaining useful life, RUL)十分必要。然而,目前的RUL预测方法未考虑轴承运行的实时健康状态,导致RUL预测精度较低。为解决上述问题,本文提出一种基于两阶段更新数字孪生和双关联动态图卷积网络(Dual correlation dynamic graph convolutional network, DC-DGCN)的轴承RUL预测方法。首先,建立具有外圈缺陷扩展特征的轴承缺陷演化模型,通过多目标优化在第一阶段得到初始缺陷扩展曲线,实现孪生模型与真实轴承的实时交互。随后,利用第二阶段的校准缺陷进一步更新全生命周期缺陷曲线。采用双向长短期记忆网络(Bi-directional long short-term memory, Bi-LSTM)将真实轴承的振动特征与孪生缺陷相关联,完成实时映射。最后,将映射得到的缺陷融入用于RUL预测的特征空间,通过所提出的DC-DGCN方法提取物理空间与数字空间特征之间的关联性,实现最终预测。实验结果表明,该方法有效提高了轴承RUL预测的真实性。

关键词轴承,剩余寿命,实时健康状态,数字孪生,图卷积网络

目录

1 引言

2 基础理论

2.1 数字孪生

2.2 动态建模

3 数字孪生模型的两阶段更新

3.1 孪生建模

3.2 两阶段动态更新

3.3 基于双向长短期记忆网络的实时映射

4 动态图神经网络

4.1 图卷积网络

4.2 动态图卷积网络

5 模型有效性验证

5.1 数据集描述

5.2 孪生模型有效性验证

5.3 两阶段动态模型更新验证

5.4 缺陷尺寸实时映射

5.5 基于数字孪生的剩余寿命预测

6 总结

1 引言

随着《中国制造2025》和工业4.0等发展战略的持续推进,智能运维已成为各国制造业发展的主要趋势[1-3]。轴承在工业设备中起着关键作用,广泛应用于轨道交通、航空航天、风力发电等行业[4-6]。同时,工业设备中的大多数轴承在恶劣工况下工作,更容易发生腐蚀、磨损等故障。轴承故障不仅会影响系统的正常运行,严重时还可能引发重大安全事故[7, 8]。准确高效的剩余寿命(Remaining useful life, RUL)预测能够指导工业设备维护工作,制定合理的维护计划,避免因轴承意外故障造成的经济损失和人员伤亡。因此,为提高轴承RUL预测的真实性,工业界和学术界正积极探索和研究更优的预测方法[9, 10]。

近年来,随着人工智能技术的进一步发展,轴承RUL预测技术逐渐成熟。总体而言,预测方法可分为以下四类:基于物理模型的方法、基于数据驱动的方法、基于专家知识库的方法以及混合预测方法。基于物理模型的方法 会考虑设备的运行环境和退化趋势,结合物理结构建立表征轴承退化状态的数学函数。这类方法往往具有较高的预测精度,但需要深入了解被监测对象的物理特性,且预测真实性很大程度上依赖于物理模型的准确建立[11]。当前机械系统日益精密和复杂,构建能够准确反映轴承退化机理的物理模型面临巨大挑战。基于数据驱动的方法可直接对传感器获取的数据进行处理,从与退化趋势相关的历史状态信号中提取特征信息。基于这些信息,在状态监测数据与轴承RUL之间建立模糊函数关系,用于RUL预测。由于无需详细解释和建模轴承的失效机理,该方法特别适用于复杂机械系统[12-15]。基于专家知识库的方法利用专家系统特有的领域知识,通过将状态监测数据与预先定义的寿命数据进行对比,得出轴承的实时RUL[16]。但由于对领域知识要求较高,这类方法的应用范围并不广泛。混合预测方法将数据驱动技术与物理模型相结合,既利用了物理模型预测的准确性,又兼具数据驱动方法的实用性,能够有效捕捉轴承退化过程中的不确定性[17,18]。然而,这类方法的模型复杂度通常较高,且可能受到物理模型要求的限制。

随着大数据的发展与普及,基于数据驱动的方法逐渐成为RUL预测的主流。图神经网络(Graph convolutional network, GNNs)凭借其独特的关联性提取能力,在该领域受到越来越多的关注。图神经网络由顶点、边和图构成:顶点代表每个节点的特征,边代表节点之间的关系或相似性,图则是所有顶点和边的综合特征表示。捕捉特征之间的关联性能够有效从监测数据中提取退化信息,从而提高网络模型的鲁棒性和准确性。然而,现有文献常将多个传感器视为节点,将传感器之间的位置关系视为边。例如,Li等[19]建立了多传感器图神经网络用于轴承RUL预测,采用余弦相似度计算节点间的边权值。在单传感器领域,现有研究往往只考虑特征通道数据单维度的关联性。例如,Yang等[20]提出了一种基于单传感器的图神经网络,以时间序列作为节点连接判断的依据。仅考虑单维度关联性的图神经网络会丢失部分关联信息,导致预测精度下降。

上述所有方法在轴承RUL预测方面均取得了一定成果,但它们都未关注轴承退化过程中的运动行为,且忽略了特征空间与轴承物理空间之间的实时交互。数字孪生(Digital twin, DT)为连接物理实体与状态监测数据提供了可行途径[21]。

DT作为一种多尺度、多学科的仿真模型,能够实时模拟真实物理实体的运动,是连接物理空间与数字空间的桥梁[22]。同时,DT包含丰富的运行数据,可从中获取真实工况下难以得到的健康状态信息,这些信息能够为机械设备的制造、改进和维护提供指导[23]。利用DT技术可大幅降低实验成本,减少资源浪费[24]。目前,已有大量研究将DT应用于轴承全生命周期管理。

Ma等[25]构建了一个融合领域自适应和注意力机制的增强元迁移学习模型,并将其与基于有限元的DT模型相结合,实现了基于少量标记目标域数据的轴承故障诊断。Qin等[26]针对轴承故障诊断中的样本不平衡问题,提出了一种基于DT的数据集生成方法。他们构建反向物理信息神经网络来调整模型参数,以优化孪生信号的匹配度,并设计边界损失辅助网络实现快速收敛。实验结果表明,该方法在不同工况和故障模式下均能提高轴承故障诊断的准确性。Li等[27]对轴承的故障演化模式进行了研究,将其分为裂纹、剥落和缺陷三类,提出了轴承全生命周期的多尺度DT模型,并通过实验验证了该模型具有良好的准确性。

DT的发展时间相对较短,目前尚无公认的最优轴承数值模型。此外,许多技术尚未成熟,如高精度模型构建、虚实数据传输、数据预处理等。在轴承RUL预测领域,DT的应用也处于逐步完善阶段。

Desai等[28]提出了一种基于摩擦学感知的DT模型,并将其与具有因果膨胀卷积的卷积神经网络相结合用于RUL预测,实现了无需磨损测量的自主维护。Liu等[29]基于振动现象构建了轴承DT模型,并采用领域对抗神经网络减小数字域与真实域之间的差异,该方法无需先验标记即可实现RUL预测。Zhang等[30]考虑到轴承RUL预测中未来数据缺失的问题,构建了基于敏感特征的DT框架,利用联合DT的集成学习Catboost方法对数据集进行补全,并通过实验验证了完整数据集的有效性。Zhao等[31]提出了一种通过DT模型生成样本数据的小样本RUL预测方法,采用多源回归对抗领域自适应方法对模拟数据进行领域转换。

对现有文献的分析表明,目前用于轴承的DT更倾向于建立仿真模型来模拟轴承的运行过程。然而,上述仿真模型并未实现DT的核心目标——实时交互与数据反馈。这导致现有RUL预测研究无法获取轴承的实时运行状态,且忽略了物理空间与数字空间数据之间的关联性。同时,大多数现有DT模型未考虑DT与物理实体之间的实时差异,无法准确反映轴承的退化过程。

为解决上述问题,本文提出一种基于DT模型两阶段更新的RUL预测方法。在孪生模型初始化后,通过多目标粒子群优化捕捉数字空间与物理空间之间的实时关联,更新孪生模型的动态响应,使其匹配最优的轴承全寿命缺陷尺寸。在第二阶段,通过校准缺陷尺寸进一步更新缺陷曲线。随后引入Bi-LSTM,将轴承的振动特征有效映射到孪生缺陷特征,以提高RUL预测的真实性。本文的主要贡献如下:

(1)提出了两阶段更新DT模型,与传统动态建模不同。该模型通过多目标优化捕捉两类空间之间的实时关联,更新模型参数以实现孪生模型与真实轴承的交互,并在第二阶段利用高保真校准缺陷进一步更新缺陷曲线。

(2)构建了基于Bi-LSTM的实时映射网络,并以DT模型两阶段更新得到的缺陷尺寸作为数据集对其进行训练。

(3)提出了基于双关联的动态图神经网络,能够更好地提取物理空间与数字空间之间的特征关联性。    

2 基础理论

DT为轴承RUL预测提供了新的可能,其核心是动态建模。构建能够充分反映轴承真实运行状态的动态模型,有助于DT更好地发挥作用,从而更准确地指导物理空间的RUL预测。    

2.1 数字孪生

Grieves在产品生命周期管理课程中首次提出数DT概念,引入了镜像空间模型[32]。此后,DT受到全球广泛关注,并在产品故障预测与健康管理领域不断发展。DT模型是物理对象的虚拟表示,通过与物理实体的实时交互,使孪生模型不断接近物理实体的真实运行状态。当达到一定相似度阈值后,可通过多种方法从孪生模型中获取数据,进而对物理实体进行优化,指导维护决策[33]。DT具有提高可靠性、降低成本等优势,已在多个领域成功应用[34, 35]。基于DT的轴承全生命周期管理一般流程如图1所示。

 

图1 DT

2.2 动态建模

本文引入赫兹接触理论建立轴承动力学模型,该模型能更准确地反映轴承局部缺陷的扩展过程,所建立的动力学模型由微分方程推导得出:  
     

其中,M为系统质量,C为阻尼系数,K为刚度,F为激励载荷。

为方便研究轴承动力学特性,假设滚动体与滚道之间的接触变形为非线性弹簧-阻尼单元。当滚动体经过缺陷位置时,会产生时变位移激励,进而引发动力学响应。本文建立模型时做出如下必要假设:

(1)轴承运行过程中温度恒定;

(2)滚动体等距排列,质量可忽略不计,且相互之间无作用力;

(3)仅考虑轴承径向方向的作用力,忽略轴向振动;

(4)轴承经过缺陷时发生局部接触变形,但轴承整体尺寸不变。

3 数字孪生模型的两阶段更新

本节构建了具有参数更新机制的两阶段DT模型,整体RUL预测框架如图2所示。首先,根据轴承的物理参数模拟其全生命周期的动力学响应;通过两阶段更新与物理实体交互,实时更新轴承的缺陷尺寸,更准确地捕捉轴承全生命周期的缺陷状态;随后,为提高DT交互的实时性,采用Bi-LSTM建立轴承物理空间特征与孪生缺陷之间的实时映射;最后,将实时缺陷融入特征空间,以提高轴承RUL预测的准确性。

 

图2 基于DTRUL预测框架

3.1 孪生建模

DT系统的动力学模型如图3所示。应用牛顿第二定律并结合式(1),可推导出考虑阻尼和恒定外载荷的滚动轴承二自由度动力学方程,其关系如式(2)所示:  
     

其中,      为轴承的等效质量,      为轴承的等效阻尼,            分别为xy方向的振动速度,            分别为xy方向的振动加速度,            分别为外圈与滚动体在xy方向的接触力分量,       、      分别为轴承的载荷分量。    

 

图3 孪生模型

将接触力视为非线性接触力,根据赫兹接触理论计算,可推导出式(3):  
     
对于圆柱滚子轴承,n=10/9;对于球轴承,n=3/2。将其代入式(2),可得到滚动轴承孪生模型的最终缺陷扩展模型,如式(4)所示:  
     
其中,      为轴承的总接触刚度,      为滚动体数量,      为滚动体      与外圈是否接触的判断因子,可由式(5)表示:  
     

   为滚道与滚动体i的接触变形,可由式(6)表示:

     

其中,    为时变位移激励,        时刻滚动体    的角位置,可通过式(7)-(9)计算:首先将轴的转速转换为弧度制,然后计算轴承保持架的转速,最后通过式(9)计算滚动体的位置。

             

其中,    为旋转系统中轴的转速,    为轴的角速度,    为滚动体直径,    为轴承的节圆直径,    为轴承保持架的转速,    为滚动体总数。 

对于时变位移激励    ,为模拟真实轴承缺陷演化产生的位移激励,通常将轴承全生命周期分为两个阶段进行计算。为简化动力学模型,忽略缺陷边缘与滚动体之间的接触变形,将缺陷视为长度为    、深度为    的矩形。本文所考虑的轴承局部缺陷扩展过程如图4所示,据此建立两个阶段的位移激励函数表达式。

 

图4 缺陷扩展过程

在初始缺陷生成阶段,由图4(a)可知,由于此时缺陷长度较小,滚动体无法完全接触到缺陷底部,其最大位移激励    始终小于缺陷深度    ,可由式(10)计算:    

此时,时变位移激励函数可表示为式(11)

其中,    为时刻        个滚动体的剩余角位置,        为缺陷中心对应的角度,    为缺陷对应的跨度角。图4(b)为缺陷扩展阶段:随着缺陷尺寸增大,最大位移激励逐渐增加,直至等于缺陷深度    。此时滚动体将完全陷入缺陷中,这种现象通常在故障中后期出现。时变位移激励函数定义如式(12)所示:

其中,    为滚动体接触缺陷底部时球心对应的跨度角。

3.2 两阶段动态更新    

基于3.1节的孪生建模,可调整各动态参数(如阻尼比、外载荷等)以匹配滚动轴承的真实动态响应。同时,为更准确地建模轴承全生命周期的缺陷尺寸,提出两阶段动态更新方法:采用多目标优化算法对孪生模型的缺陷参数进行持续迭代,得到帕累托最优前沿,使孪生模型生成的动态响应与真实响应高度匹配;在准确确定轴承全生命周期缺陷尺寸后,采用所提计算方法获取特定时刻的高保真缺陷尺寸,进而优化全生命周期缺陷曲线,最终得到精度更高的结果。    

3.2.1 多目标参数优化    

为使真实轴承与孪生轴承的动态响应差异最小化,首先获取初步的全生命周期缺陷尺寸。本文采用多目标优化对孪生模型参数进行迭代,优化目标同时考虑频域和时域相似度,以确保获得最匹配的缺陷参数。多目标优化的数学表达式如下:

     

其中,    为优化后的缺陷参数(包括缺陷长度    、缺陷深度    ),        分别为孪生轴承和真实轴承的动态响应,    表示快速傅里叶变换,    表示动态时间规整距离,    表示皮尔逊相关系数。 

多目标粒子群优化算法由Carlos A. Coello Coello等人提出,是将传统粒子群算法扩展到多目标优化问题的算法。粒子群算法作为一种群智能进化算法,其核心原理是将优化问题的最优解类比为食物位置,将粒子的运动类比为鸟类的飞行方向和速度。在种群更新过程中,每只“鸟”(粒子)通过同时考虑自身对食物位置的记忆以及群体记忆的影响来调整飞行方向,通过迭代更新种群位置,引导种群向最优解收敛,最终获得全局最优解。其具体更新过程可表示为式(14):

     

其中,    为第    次迭代后粒子i的速度向量,    为第    次迭代后粒子    的位置向量,    为惯性系数,        为[0,1]范围内的随机数,    为个体学习因子,    为社会学习因子,    为第    次迭代中粒子    的最优位置,    为第    次迭代中种群的最优位置。计算得到第    次迭代的速度后,采用式(15)更新每个粒子的位置:

     

引入自适应网格算法和帕累托支配实现基于粒子群优化的多目标优化:多目标粒子群优化算法采用自适应网格算法和轮盘赌选择法,从每次迭代得到的一系列非支配关系粒子中选择历史最优解,在外部存档中存储分布更均匀的帕累托前沿。

3.2.2 校准缺陷获取

通过多目标优化求解后,可得到轴承全生命周期缺陷尺寸。为进一步减小其与真实轴承缺陷尺寸的差异,需采用尽可能接近真实的缺陷尺寸对3.2.1节得到的全生命周期缺陷曲线再次调整。然而,在轴承运行过程中拆解轴承测量缺陷尺寸既不合理又成本高昂。为解决这一问题,本文引入文献[36]中的方法,通过分析轴承时域振动信号确定特定时间段内更准确的缺陷尺寸。    

滚动体经过局部缺陷时的物理过程如图5所示:当滚动体到达A点时,开始进入缺陷区域;此时由于接触面积减小,动态响应降低,直至到达最低局部点B(B点位置无需用于缺陷尺寸计算,故图中未标出);随后动态响应逐渐增大,并在C点达到峰值,此时滚动体即将脱离滚道;经过C点后,滚动体穿过缺陷到达D点,与缺陷边缘碰撞并产生高频响应。图中部分变量的几何关系如式(16)所示:

     

其中,        为A点、C点相对于缺陷入口E的角度,    为轴承节圆半径,    为保持架旋转频率,    为轴承从C点旋转到D点的时间,    为滚动体半径。由上式可知,从振动信号中提取各点位置后,可求得缺陷的跨度角    ,再通过式(17)得到精度更高的轴承缺陷长度:

     

获取部分时刻的校准缺陷尺寸后,采用式(18)对第一阶段生成的轴承全生命周期缺陷曲线进行修正,得到最终的两阶段更新缺陷尺寸:

     

其中,    为与第一个校准缺陷对应的第一阶段优化缺陷尺寸,    为轴承运行至失效的总时间,    为获取的校准缺陷数量,    为第    个时刻的校准缺陷尺寸。 

 

图5 滚动体通过缺陷区域

需说明的是,本文仅选择特定时刻计算校准缺陷,原因如下:1.上述校准缺陷计算过程基于轴承的双脉冲现象,而在轴承故障早期,该现象通常难以捕捉,因此该方法仅在缺陷扩展到一定尺寸后才具有较高精度;2.计算标定缺陷的成本较大,对轴承每一时刻计算标定缺陷是不切实际的,也是不合理的;3.计算每个轴承时刻的标定缺陷是不现实的,也是不合理的。

3.3 基于双向长短期记忆网络的实时映射    

为增强物理实体与孪生实体之间的交互性,加快信息交互速度,考虑到缺陷尺寸的增长是一个连续过程,且某一时刻的缺陷尺寸与其前后状态相关联,本文提出采用双向长短期记忆网络建立虚实空间之间的关联。该关联基于神经网络优异的学习能力,可接收轴承实时振动信号,通过训练后的神经网络输出孪生模型的实时缺陷尺寸,并反馈给真实轴承,为其RUL预测提供指导。    

由于缺陷尺寸也可视为一维特征,为实现最优关联效果,需选择合适的特征进行映射。均方根值受噪声影响较小,且能有效描述动态响应的大小,因此本文选择均方根值作为缺陷映射的特征。具体步骤如下:(1)通过3.2节的两阶段模型更新,获取轴承全生命周期缺陷尺寸;(2)计算轴承全生命周期的均方根值;(3)以均方根值为输入、全生命周期缺陷尺寸为输出,训练双向长短期记忆网络;(4)将轴承实时均方根值直接输入模型,得到实时映射的缺陷尺寸。    

4 动态图神经网络

为将孪生模型的退化特征更好地融入真实轴承的RUL预测,考虑到图卷积网络(Graph Convolutional Networks, GCNs)能够捕捉特征间的关联关系,且聚合该关联关系可提高轴承RUL预测的准确性,本文提出一种考虑双关联的动态图卷积网络(Dual correlation dynamic graph convolutional network, DC-DGCN)。基于DGCN的轴承RUL整体流程如图6所示,该网络主要包括三部分:真实轴承特征提取与孪生轴承特征映射、动态空间关联提取、时间关联提取。首先,通过与真实轴承交互得到最优全生命周期孪生模型,并将真实轴承特征与孪生模型缺陷特征融合;然后,划分时间步数据,构建动态空间关联图以提取融合特征通道的动态空间关联,并将其输入带有门控循环单元(Gated Recurrent Unit, GRU)的DGCN模块,提取双关联特征;最后,通过全连接层输出轴承RUL。

图6 DC-DGCN模型的框架

4.1 GCN

GCN可聚合具有关联信息的相邻节点,以更好地表示节点特征。图的结构可定义为    ,其中    表示节点的特征信息,    为节点数量(即特征数量),    为特征维度(即时间步长度),    表示节点之间的连接关系,    表示节点间关联的邻接矩阵。本文采用余弦相似度定义节点间的相似度矩阵,余弦相似度矩阵    中的元素    可通过式(19)计算,其中    为第    个轴承第    次采样数据第    行的特征,    为向量的二范数:

     

对余弦相似度矩阵    的元素应用式(20)定义的函数    ,可得到邻接矩阵    的元素    ,其中    为人工定义的阈值。    表示在时间关联图中两个时间节点无边线连接,即无时间关联性;    表示在时间关联图中两个时间节点有边线连接,即存在时间关联性。

     

计算空间关联邻接矩阵,并将其与输入特征矩阵结合。GCN模型将在傅里叶域构建图滤波器,该滤波器可作用于图节点,聚合一阶邻域特征,从而捕捉节点间的关联信息。归一化图拉普拉斯矩阵定义如式(21)所示:

     

其中,    为度矩阵,        为单位矩阵。GCN模型的信息传递过程可表示为式(22):

     

其中,     为带有自环的邻接矩阵(即在原始邻接矩阵中加入自连接),    为第    层的可训练权重矩阵,    为ReLU激活函数,    为第    层的输入(当    时,        为输入特征矩阵)

4.2 DGCN

由4.1节可知,GCN通过聚合相邻节点能够很好地捕捉特征间的空间关联,但仅依靠GCN无法获取特征间的时间依赖关系(即时间关联)。因此,本文设计DGCN模块,以同时获取融合特征通道的时空关联,并将其用于RUL预测。    

图7展示了动态图卷积网络模块的结构,其中      为时刻      的输入,      为隐藏层状态,                  分别为细胞状态、更新门和重置门。    

 
图7 DGCN模块的构建  

该模块的具体计算过程可表示为式(23),其中    为Sigmoid激活函数,    为图卷积网络处理过程,        分别为可训练权重矩阵和偏置:

     

综上,该模块在保留有效历史信息的同时,能够捕捉当前时空域的关联关系,因此其隐藏层同时包含融合特征通道的时间关联和空间关联信息。

5 模型有效性验证

为验证本文所提模型的有效性,本节采用公开数据集进行实验。所有实验均在搭载I7-11800H@2.30 GHz处理器和NVIDIA GeForce RTX 3090显卡的服务器上运行,RUL预测算法采用Pytorch框架实现。

5.1 数据集描述

采用XJTU-SY数据集的轴承退化数据验证所提两阶段动态更新DT模型的有效性,实验平台如图8所示。该数据集提供了15个轴承在三种不同工况下的全生命周期数据,表1列出了三组轴承的详细信息,采样间隔为1分钟,每个样本采样时长为1.28秒,采样频率为25.6kHz。

 
图8 XJTU-SY轴承实验台  

XJTU-SY数据集轴承工况

 

表2列出了实验轴承的实际参数,为最小化孪生模型与真实轴承的差异,将孪生模型的动力学参数设置为实际值。    

实验轴承主要参数    

 

XJTU-SY数据集每个样本包含32768个采样点,为便于观察和分析,每次仅使用总采样点中的10000个数据点优化孪生模型的动态参数。动态仿真参数如表3所示。采用上述动力学参数在Simulink中进行仿真,并在Matlab中使用ode45求解器执行后续多目标优化。    

3动态仿真主要参数    

 

5.孪生模型有效性验证    

孪生模型的动态响应如图9所示,其与真实信号的包络分析对比如图10所示。可以看出,由于局部缺陷的存在,孪生模型的动态响应具有明显的周期性。由于外圈缺陷位置固定,产生的冲击信号间隔均匀,滚动体经过缺陷产生冲击的周期        为外圈故障特征频率)。结合实验轴承的实际参数,通过式(24)可计算出其外圈故障特征频率为108Hz。由图10可知,孪生信号与真实信号的包络信号基本一致。由于真实轴承运行环境中存在噪声,其动态响应表现出不规则性,且振幅略大于孪生信号;相比之下,孪生信号更为规则。仿真数据的包络谱清晰地显示出轴承故障特征频率及其倍频,与真实轴承的包络谱差异极小。因此,可认为该孪生模型能够在一定程度上反映真实轴承的运行状态。

     

其中,    为滚动体数量,    为轴承旋转频率,    为滚动体直径,    为轴承节圆直径,    为轴承接触角。

 

时域孪生信号    

 

图10 实际信号与孪生信号的包络分析对比:(a)包络信号,(b)包络谱

5.3 两阶段动态更新模型验证    

DT的关键过程是与物理空间的实时交互。为使孪生模型更准确地反映轴承的实时健康状态,采用两阶段孪生缺陷尺寸更新来跟踪真实缺陷尺寸。    

5.3.1 多目标优化求解    

在验证孪生模型有效性后,利用真实轴承的动态响应对其进行两阶段动态更新,以获取最优全生命周期缺陷尺寸。本节选择XJTU-SY数据集的Bearing1_1和Bearing1_3作为测试对象,通过两阶段更新得到两个轴承的全生命周期缺陷曲线,用于后续实时映射。    

多目标优化后得到的第一阶段轴承缺陷曲线如图11所示,每个时刻对应的不同缺陷尺寸为此时多目标优化求解得到的帕累托最优解。    

 

图11 多目标优化结果:aBearing1_1,bBearing1_3    

对每个时刻得到的缺陷尺寸取平均值,得到基于多目标优化的轴承第一阶段全生命周期缺陷尺寸曲线,如图12所示。

 

12第一阶段的缺陷尺寸(a)Bearing1_1,(b)Bearing1_3    

5.3.2 校准缺陷尺寸计算    

获取第一阶段轴承缺陷尺寸后,计算校准缺陷尺寸并进行第二阶段更新至关重要。文献研究表明,仅当轴承动态响应表现出明显的双脉冲现象时,该计算方法才具有较高精度。从第一阶段轴承缺陷尺寸曲线可以看出,Bearing1_1在78分钟时动态响应发生较大变化,可认为此时轴承缺陷尺寸已满足产生双脉冲现象的条件。因此,计算78分钟后的多个缺陷尺寸,对第一阶段轴承缺陷曲线进行更新;Bearing1_3采用类似方法进行推断和计算。表4列出了Bearing1_1的计算校准缺陷尺寸及其对应时刻的第一阶段缺陷尺寸。    

Bearing1_1高精度校准缺陷    

 

采用式(18)对第一阶段缺陷曲线进行修正,得到第二阶段更新缺陷尺寸,如图13所示。观察图13可知,第一阶段得到的轴承缺陷曲线能够大致跟踪整体趋势,但在轴承故障后期,多目标优化算法得到的缺陷尺寸曲线与两阶段更新后的曲线存在一定差异。原因是在缺陷较大的时间段内,轴承动态响应的不稳定性更强,导致缺陷匹配出现误差。两阶段校准缺陷可指导第一阶段缺陷曲线进行再次更新,从而得到高保真的轴承全生命周期缺陷曲线。

 

1二阶段更新过程(a)Bearing1_1,(b)Bearing1_3    

5.3.3 全生命周期缺陷尺寸验证    

为进一步验证两阶段更新得到的全生命周期缺陷尺寸的有效性和准确性,将其与最新研究成果进行对比,包括Li等[27]提出的多尺度故障演化DT模型、Shi等[32]提出的新型DT模型以及Qin等[38]提出的基于CycleGAN的DT模型。需说明的是,为便于对比,对部分相关文献的缺陷尺寸曲线进行了处理,但保留了其主要趋势和特征。对比结果如图14所示,可以看出,两阶段更新后的全生命周期缺陷尺寸曲线与其他相关文献的曲线趋势基本一致,验证了本文所提DT模型的有效性。由于本文采用多目标优化拟合第一阶段曲线,缺陷尺寸曲线的性能更稳定,有利于后续RUL预测。

 

1缺陷尺寸的对比    

5.4 缺陷尺寸实时映射    

本节选择Bearing1_1数据的80%作为双向长短期记忆网络的训练集,剩余20%的Bearing1_1数据及Bearing1_3数据作为测试集,以验证缺陷曲线的准确性和模型将振动特征实时映射到缺陷尺寸的有效性。    

本文设计的双向长短期记忆网络在输出时仅分别获取前向和后向部分的最后一步结果。通过调整隐藏层维度确定最优配置,实验表明隐藏层维度对预测结果的影响较小,因此将隐藏层维度设置为10,在保证预测性能的同时平衡参数数量。此外,采用Adam优化器对模型进行优化。    

模型的映射效果如图15所示,可见映射曲线与两阶段更新得到的缺陷尺寸曲线差异极小,验证了映射模型的有效性,该模型能够很好地将真实轴承的动态响应映射到孪生模型的缺陷尺寸。这种实时映射能够以优异的实时性和准确性揭示轴承振动信号与运行缺陷尺寸之间的关系。随后,采用训练后的模型获取轴承的实时运行状态,并将该信息融入RUL预测。

 

图15 基于Bi-LSTM的缺陷尺寸映射:(a)Bearing1_1,(b)Bearing1_3    

5.5 基于数字孪生的剩余寿命预测    

上述小节从时域、包络谱等方面验证了DT模型的有效性,并通过两阶段更新校准了轴承的全生命周期缺陷尺寸,训练得到性能较优的映射模型。为实现整体DT框架的实时交互目标,本部分将模型映射得到的缺陷尺寸作为特征加入RUL预测,并通过与其他先进预测方法对比,验证所提方法的优越性。    

5.5.1 超参数设置    

在本节的RUL预测实验中,共提取8个时域特征、11个频域特征和8个时频域特征。其中,8个时域特征包括方差、峰度、熵、最大最小值、峰峰值、偏度、均方根和标准差;频域特征通过对原始信号进行快速傅里叶变换提取,包括中心频率、平均频率、最大振幅、最大功率谱密度、平均功率谱密度、功率谱密度方差、偏度、峰度、频率标准差、频率分布不对称系数和频率分布峰度;时频域特征通过采用Daubechies小波基函数进行三层小波包分解得到,8个特征分别对应不同子带的能量占比。研究表明,这些特征在轴承RUL预测中具有有效性[37]。    

所提模型的超参数包括隐藏层维度和时间步长:时间步设置过小时,门控模块难以有效捕捉时间信息;时间步长设置过大则会导致模型参数过多,延长训练过程,因此将时间步长推荐范围设置为20-35。隐藏层维度决定模型的信息记忆能力,设置范围为2-10。通过网格搜索法优化,综合考虑预测性能后,最终确定时间步长为30,隐藏层维度为6。    

5.5.2 实验设计和结果    

基于XJTU-SY数据集设计三组RUL预测任务,以验证所提方法在不同工况下的有效性,具体配置如表5所示。    

预测任务的具体设置    

 

本文采用故障预测与健康管理领域广泛使用的均方根误差(Root mean square error, RMSE)和平均绝对误差(Mean square error, MAE)作为评价指标,衡量模型预测水平。这两个指标反映真实值与预测值之间的平均偏差,数值越小表示预测准确性越高。其中,均方根误差对异常值更敏感。此外,设计Score指标评价模型预测精度,该指标对预测值与真实值的差异高度敏感,且当真实值大于预测值时Score会显著增大,这一特性更符合实际需求。指标计算公式如下:

                 

其中,    为测试样本数量,    为第    个样本的真实值,    为第    个样本的预测值。

为验证本文获取的映射缺陷尺寸的有效性,将实验分为两种场景:S2表示仅提取上述特征;S1表示在提取上述特征的基础上加入全生命周期缺陷尺寸。    

在对比模型选择方面:首先选择仅提取原始特征序列时间特征的时间卷积网络(Temporal Convolutional Network, TCN);其次,为验证图神经网络特征提取的优越性,选择两种具有空间提取能力的模型——卷积门控循环单元(Convolutional GRU, ConvGRU)和卷积长短期记忆网络(Convolutional LSTM, C-LSTM),这两种模型既能提取时间特征,又能利用卷积神经网络提取空间特征;最后,选择图卷积网络(GCN)进行对比,以验证所提模型在双关联特征提取方面的优越性。此外,为更好地体现所提模型的优越性,选择文献[20]中的ChebGCN-LSTM模型进行对比,该模型将特定时刻的特征序列作为节点特征,以时间顺序作为节点间连接的判断依据,并采用长短期记忆网络获取时间依赖关系。    

训练过程采用Adam优化器,批大小设置为32,学习率设置为0.001,损失函数选择均方误差。    

各模型在两种场景下的预测结果如表6和表7所示。场景S2的对比数据表明,所提方法的预测准确性总体优于其他模型:其平均均方根误差比表现最优的对比模型ConvGRU低0.011,这得益于图神经网络提取的关联关系能够更有效地揭示轴承退化过程,最终提高RUL预测准确性;此外,所提模型的均方根误差比标准图卷积网络和ChebGCN-LSTM分别降低0.026和0.013,表明双关联提取考虑了更广泛的特征信息,验证了所提模型的优越性。同时,同时提取时间和空间特征的模型比仅提取时间特征的模型表现出更优的预测性能。    

综合场景S1来看,时间卷积网络在加入全周期缺陷尺寸作为特征后,预测性能未得到显著提升,表明仅提取时间特征的模型难以将实时缺陷尺寸很好地融入RUL预测;而ConvGRU、ChebGCN-LSTM和C-LSTM在加入该特征后,预测性能均有不同程度提升,总平均均方根误差分别降低0.11、0.10和0.12,验证了所提两阶段更新DT模型的有效性——通过映射得到的实时缺陷尺寸可提高模型的RUL预测能力。所提方法在S1场景下的平均均方根误差比S2场景降低0.020,降幅显著,表明与采用卷积神经网络进行空间特征提取相比,所提方法能更有效地提取实时缺陷特征与其他特征之间的关联关系。值得注意的是,所提模型的均方根误差降幅大于ChebGCN-LSTM,这得益于动态图的构建,使模型能够更有效地跟踪轴承退化趋势,获取适应性更强的关联信息。    

表6 不同场景下各模型的RMSE值

 

表7 不同场景下各模型的MAE值

 

Score指标对比如图16和图17所示,可见所提模型的Score值最低,表明其预测精度最高,预测结果更有利于指导实际工作。尽管ConvGRU的均方根误差表现较好,但在Score指标上表现较差,说明其预测曲线难以很好地应用于实际场景。

 

1场景S2Score

 

1场景S1的Score值

为进一步验证所提方法的有效性,图18和图19分别展示了部分轴承在场景S2和S1下的RUL预测结果曲线。可以看出,大多数模型在整体趋势上能够反映轴承退化过程,但时间卷积网络的预测精度未能达到理想效果;与仅提取时间关联的模型相比,同时提取时空关联的模型表现出更高的准确性;且在加入全寿命缺陷尺寸作为特征后,各模型的预测曲线更接近实际寿命曲线。在两种场景下,动态图卷积网络均能比其他模型更好地跟踪轴承退化趋势,具有更优的RUL预测能力。    

两种场景的对比表明,动态图卷积网络在加入实时缺陷尺寸后,预测曲线的改善更为明显,说明所提方法通过融合缺陷尺寸特征,有效捕捉了轴承的整体退化特征,提高了轴承RUL预测的准确性,尤其在轴承寿命后期效果显著。这对RUL预测至关重要,因为轴承寿命前期的RUL通常无需过多关注,而准确的寿命后期RUL预测能够帮助提高设备可靠性,减少因突发故障造成的损耗。    

此外,需注意的是,所提方法在部分轴承的预测中,加入缺陷尺寸特征后预测效果未得到显著提升,推测原因是这些轴承在运行过程中未出现明显的外圈缺陷,而本文构建的孪生模型仅描述了外圈的缺陷扩展过程。    

为更清晰地说明所提方法的有效性和优越性,对比所提模型在两种场景下的预测结果,图20给出了部分轴承在两种场景下的预测曲线及误差大小。场景S1的模型预测精度更高,在预测后期的误差小于场景S2。这一优势源于动态图卷积网络能够有效利用全生命周期缺陷尺寸,且模型捕捉的关联关系显著提升了RUL预测性能。    

 

图18 场景S2下各模型的预测曲线:aBearing1-3,(b)Bearing1-4,(c)Bearing2-4,(d)Bearing3-3    

 

图19 场景S1下各模型的预测曲线:(a)Bearing1-3,(b)Bearing1-4,(c)Bearing2-4,(d)Bearing3-3    

 

2所提模型两种场景下的对比

6 总结

针对当前轴承RUL预测未考虑运行实时健康状态的问题,本文提出一种DC-DGCN轴承RUL预测方法。首先,通过构建含局部缺陷的孪生动力学模型来模拟真实轴承的运行过程,并采用第一阶段多目标优化方法实时更新孪生模型参数,使孪生模型生成的信号与真实轴承的动态响应更加吻合。在获取轴承全生命周期缺陷尺寸后,通过第二阶段校准进一步优化缺陷尺寸,并引入Bi-LSTM建立实时映射关系。最终,结合从DT中获取的缺陷尺寸特征,利用所提出的双关联动态图卷积网络实现轴承RUL预测。

实验结果表明,本文提出的动态孪生模型能够更准确地复现真实轴承的运行状态;在双关联动态图卷积网络的作用下,通过两阶段更新获取的缺陷特征显著提升了轴承RUL预测的准确性。此外,实验验证表明,本文所提方法的性能优于其他先进模型。

该模型将实时健康状态监测与先进数据分析技术相融合。随着人工智能技术的不断发展,该模型可进一步融入智能决策系统,实现设备的自动化维护与故障预警,有望为工业制造等领域的智能化转型提供重要支撑。

 

编辑:陈宇航

校核:李正平、陈凯歌曹希铭、赵学功、白亮、任超、海洋、赵栓栓、陈莹洁、Tina、王金、赵诚、肖鑫鑫、张优

该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除

来源:故障诊断与python学习


ACTMechanicalSystem振动碰撞非线性化学旋转机械航空航天轨道交通MATLABpython海洋裂纹参数优化理论电机化机数字孪生
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-10-26
最近编辑:10月前
故障诊断与python学习
硕士 签名征集中
获赞 87粉丝 152文章 328课程 0
点赞
收藏
作者推荐

域泛化开源代码分享| 平衡差异性与一致性:单域泛化故障诊断方法研究

本期给大家推荐的文章是《平衡差异性与一致性:单域泛化故障诊断方法研究》。本文以当前对于单域泛化故障诊断方法研究过程中的局限性为切入点,提出一种建模“差异性—一致性”对抗交互机制的单域泛化故障诊断框架。论文基本信息论文题目:Balancing Discrepancy and Consistency:Adversarial Single Domain Generalizationin Fault Diagnosis论文期刊:IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS论文日期:2025年作者:Guowei Zhang , Xianguang Kong , Qibin Wang , Jingli Du , and Hongbo Ma机构: College of Mechanical and Electrical Engineering, Xidian University, Xi’an 710071, China; State Key Laboratory of Electromechanical Integrated Manufacturing of High-performance Electronic Equipments, Xidian University,Xi’an 710071, China 第一作者简介: 张国伟,西安电子科技大学机电工程学院博士研究生,新加坡国立大学联合培养博士。研究聚焦于旋转机械智能故障诊断与跨域泛化建模,致力于提升诊断模型在复杂工业环境中的鲁棒性与实用性。以第一作者或通讯作者在IEEE TII、MSSP等国际权威期刊发表SCI论文10篇。入选首届中国科协青年人才托举工程博士生专项计划,曾获国家奖学金等多项荣誉。目录摘要 1 引言 2 问题定义 2.1 基于SDG的故障诊断问题建模 2.2 训练目标与理论动机 2.2.1 基于分布扰动的鲁棒学习 2.2.2 领域自适应理论的泛化保证 2.2.3 对抗结构设计的稳定性分析 3 所提出方法 3.1 基于域差异感知的可学习统计扰动 3.2 一致性引导的协方差与语义对齐 3.3 训练与推理流程 4 实验验证 4.1 数据集介绍 4.2 实验设置 4.2.1 参数配置 4.2.2 对比方法 4.2.3 评估指标 4.3 结果与讨论 4.3.1 主要结果分析 4.3.2 消融实验分析 4.3.3 计算开销评估 5 结论 摘要单域泛化(Single-domain Generalization, SDG)故障诊断旨在仅利用单一工况下的数据训练模型,使其具备在未知工况下的推广能力。尽管该策略具有良好的应用前景,但仍面临诸多挑战。现有研究主要通过生成新域特征以扩展数据分布,并提取域不变特征,然而在实践中存在三方面问题:(1)域扩展与域不变特征提取相互独立,限制了两者的协同作用与泛化能力;(2)域扩展过程常采用随机或方向不一致的增强策略,易引入噪声或产生不可靠特征;(3)域不变特征提取方法多侧重于全局或类别层面对齐,忽略了特征间关系的优化。为此,本文提出一种基于差异性与一致性对抗交互的SDG故障诊断框架,通过反复对抗训练动态平衡域间差异性与一致性,避免模型泛化受单向偏差影响。具体而言,设计了一个可学习的统计特征扰动模块以感知域间差异,引导模型生成能有效表征域差异的统计扰动,从而提升域扩展的可靠性。此外,提出协方差对齐与语义对齐策略,在原始特征与扰动特征之间保持全局及类别层面的一致性。协方差对齐用于减缓扰动特征引入的风格变化,并保留关键诊断信息,显式语义约束则进一步提升模型预测结果的一致性。最后,本文在五个公开数据集上开展大量实验,验证所提方法的有效性。1 引言随着工业领域加速迈向智能制造,设备智能化程度和运行标准不断提升。旋转机械作为关键组成部分,其运行需满足更高的可靠性与安全性要求。在此背景下,人工智能,尤其是深度学习,成为智能设备健康管理的重要技术手段,广泛应用于状态监测、故障诊断和剩余寿命预测等任务中,具备自动特征提取、模式识别鲁棒性强、预测建模精度高等优势。早期的智能故障诊断研究主要依赖于在固定工况下的监督学习,假设训练与测试数据分布一致。然而,在实际工业场景中,设备运行常伴随负载、转速及环境条件变化,导致数据分布发生偏移,严重影响模型的诊断性能。 为应对上述问题,领域自适应(Domain Adaptation, DA)方法被提出,用于将标注源域的诊断知识迁移至无标签但相关的目标域,通过对齐特征分布实现跨域泛化。尽管DA方法在跨工况故障诊断中取得一定进展,但通常依赖于训练阶段访问目标域数据,这一假设限制了其在工业场景中的应用价值。图1 不同领域泛化故障诊断方法流程对比图为摆脱对目标域数据的依赖,领域泛化(Domain Generalization, DG)故障诊断方法提出在训练阶段引入多个源域数据,使模型具备迁移至未知目标域的能力,如图 1(a)所示。此类方法主要分为两类:一类通过数据增强提升样本多样性,另一类从多域数据中提取域不变特征。总体来看,传统DG方法多依赖于来自多个源域的故障数据,但在实际工业中,多源故障数据往往难以获取,难以满足多源域泛化方法的训练需求。 近年来,单域泛化(Single-domain Generalization, SDG)故障诊断方法受到广泛关注,其核心在于无需多源训练数据,仅基于单一源域完成对多个未知目标域的故障识别。这一策略大幅降低了模型训练与部署所需数据量,缓解了工业场景中普遍存在的数据匮乏问题。因此,工业用户可在无需大量数据采集的前提下,快速部署智能诊断系统,降低应用成本并简化系统结构。目前,SDG 故障诊断研究主要通过数据或特征增强扩展源域分布,并结合域不变特征提取策略。Wang 等人提出风格补全模块,通过生成多样化分布样本提升模型泛化能力,缓解了SDG场景中目标样本分布不足的问题。Zhao 等人提出基于对抗互信息的SDG网络,通过生成与源域差异较大的伪目标域,提升模型泛化能力,有效缓解未知域偏差问题。Wang 等人提出多尺度风格生成与对抗对比网络,通过域生成模块构造多样化扩展域样本,并结合多尺度卷积与风格学习策略提取域不变特征,实现单一工况下的泛化诊断。Zhang等人构建伪目标域以模拟域偏差,推理域偏差原因并执行自适应迁移,解决了现有方法对源域与目标域偏差估计不足的问题。Li 等人提出基于域扩展融合的SDG框架,通过域扩展模块、双分支特征融合网络、多分类器状态识别模块以及对抗对比训练策略,在未知工况下提升故障诊断的鲁棒性与准确性。Jiang 等人提出基于不确定性自举的对抗增强域网络,通过生成模块与跨域辅助模块构建增强域,结合对抗训练与不确定性估计,解决了SDG场景下未知故障模式识别与跨域诊断的问题。如图1(b)所示,为应对更具挑战性的SDG任务,当前方法普遍采用数据或特征增强技术提升样本多样性,并结合域不变特征提取以提升泛化能力。尽管相关研究取得显著进展,但当前SDG方法仍面临三项关键挑战:现有方法普遍将域扩展与域不变特征提取视为两个相互独立的过程,实质上将SDG问题转化为多源域泛化问题。尽管此类方法在一定程度上提升了模型的泛化能力,但两阶段之间联系薄弱,限制了整体框架的协同效能。在实际工业环境中,设备运行条件(如负载、转速及环境因素)往往呈现不可预测的动态变化,现有方法未能充分挖掘潜在域迁移带来的多样性,导致模型整体泛化性能受限。多数方法在进行域扩展时,通常采用随机扰动方式,或刻意增强与源域分布的差异性。然而,这并不能确保生成特征与目标域内在分布保持一致。此外,增强方向缺乏可控性,难以实现特征学习的有效优化。传统域不变特征提取方法多基于距离度量或对比学习,但往往忽略了全局分布对齐与类别层面特异性之间的交互关系。尤其是全局对齐过程可能引入无关或误导性特征,导致与任务相关的关键信息丢失,影响诊断准确性。 为应对上述限制,源于对抗学习思想的对抗训练策略,已成为无监督领域自适应中提升模型泛化能力的重要手段。对抗式领域自适应方法通过极小极大优化过程对齐不同域间的特征分布,实现域不变表示的学习。同时,对抗式数据增强策略利用对抗扰动构建虚拟目标域分布,以提升模型在域偏移情形下的鲁棒性。受此启发,本文将SDG问题建模为“差异性–一致性”之间的对抗优化过程,如图1(c)所示。首先最大化原始特征与扰动特征分布之间的差异性以模拟多样化的域偏差,随后引入一致性引导下的协方差与语义对齐机制,增强模型预测的稳定性。通过差异性与一致性之间的对抗训练,逐步提升模型的泛化能力。上述机制协同作用,有效增强模型在工业应用中的可靠性与鲁棒性。本文的主要贡献如下: 提出一种建模“差异性–一致性”对抗交互机制的SDG故障诊断框架。该框架借助对抗训练思想,在迭代优化过程中动态平衡域间差异性与一致性,有效缓解单向泛化偏差问题。模型能够逐步学习更加稳健且具判别性的特征表示,保障在复杂多变工况下的诊断精度与稳定性,为工业场景中运行条件频繁变化的智能诊断系统部署提供切实可行的解决方案。设计一个可学习的统计扰动模块(Learnable Statistical Perturbation, LSP),用于建模域间差异感知。该模块引导模型学习能反映工况变化的统计扰动特征,有效捕捉工业设备在不同工作环境下的域特异性属性。与传统固定增强策略相比,所提模块基于域特异性距离度量生成更具针对性和目标导向的扰动,从而提升特征增强的可控性与相关性。提出一种一致性引导下的协方差与语义对齐策略,用于有效解耦域间变化与语义特征之间的耦合关系。该策略在表示层面对原始样本与增强样本的全局特征分布进行对齐,既保留关键诊断信息,又缓解特征扰动带来的风格变化。同时引入显式语义约束,进一步提升模型预测结果的一致性与稳定性。2 问题定义2.1 基于SDG的故障诊断问题建模在故障诊断场景下,单域泛化任务可表述为从输入数据 到预测标签 的映射学习,其中 来自单一源域 。设 为模型集 合,训练数据 服从源域分布 ,目标是从中学习到一个模型 ,使其能够在未知目标域分布 上具备良好的泛化能力。因此,SDG 可形式化为如下优化问题: 其中, 表示交叉熵损失函数, 表示在目标域分布上的期望风险。由于SDG故障诊断的主要挑战在于无法获取目标域分布 ,因此经验风险最小化(Empirical Risk Minimization, ERM)成为求解式1的基本策略,其通过最小化源域数据上的平均误差来训练模型: 其中, 表示训练集的经验分布。然而,已有研究表明,基于ERM的训练方法在泛化场景下表现出较差的鲁棒性,领域自适应理论指出,其根本原因在于源域训练风险与源目标域之间的分布差异。2.2 训练目标与理论动机为克服前述问题,受对抗式数据增强研究启发,本文将 SDG 问题建模为“差异性–一致性”之间的对抗优化过程。 其中, 表示由扰动函数 对源域潜在特征生成的虚拟目标域分布, 与 分别为源域特征与虚拟目标域特征, 表示两组特征的对齐损失, 为扰动强度的超参数。然而,从实践角度看,如何选择合适的域偏移幅度 具有挑战性,因此我们引入拉格朗日松弛形式,将其转化为如下优化目标: 其中, 表示概率分布空间中的距离度量, 为惩罚因子,用于控制扰动的幅度。从式4可见,该优化目标旨在最大化源域与虚拟目标域间的分布差异,同时最小化模型在目标分布上的风险。式4所定义的训练目标建立在两个经典理论基础之上:分布鲁棒优化理论与领域自适应泛化界理论。下文从三个方面对该目标函数的理论基础与结构设计进行解释。2.2.1 基于分布扰动的鲁棒学习为克服式2中 ERM 策略的局限性,本文引入分布鲁棒优化思想。参考Sinha等人之前研究,我们考虑所有位于以源域分布 为中心、半径为 的Wasserstein球体内的扰动分布集 合 ,其鲁棒优化目标可表示为: 其中, 表示在距离度量 下,从 出发半径为 的可接受分布集 合。为实现可行的优化,本文采用扰动函数 对源域进行变换以构造虚拟分布 ,并通过拉格朗日松弛策略引入正则项 ,从而得到式4的优化形式。该策略通过最小化潜在最坏情况下的预测风险,使模型在不同形式的域偏移下具备更强的鲁棒性与泛化能力。2.2.2 领域自适应理论的泛化保证从领域自适应理论角度出发,本文目标函数亦符合中提出的泛化上界理论。该理论指出目标域误差可由以下上界描述:其中, 与 分别为假设 在目标域与源域上的分类误差; 为源域与目标域在假设空间 下的 -距离,衡量两域分布差异; 表示源-目标共同最优假设下的不可约误差。该上界表明,若欲降低目标域误差,不仅需减少源域经验误差 ,亦需控制源目标域间的分布差异。由于 SDG 场景下无法访问真实目标域 ,本文采用虚拟目标域 近似代替,并在式4中引入特征对齐项 ,以缓解源域与扰动域之间的语义偏移。2.2.3 对抗结构设计的稳定性分析尽管式4中采用了极小极大的对抗优化结构,其并非典型生成对抗网络(GAN)式动态博弈模型。本文采用固定结构分支引入扰动,仅对浅层特征进行扰动变换,而深层特征则保持与原始流一致,并实施特征对齐。该结构实现了扰动生成与优化过程的解耦,有效规避了对抗优化过程中常见的不稳定性问题,从而实现了具有结构约束的稳定对抗优化目标。综上所述,式4所构建的训练目标在分布鲁棒性与领域泛化理论的指导下,统一引入分布扰动与特征一致性约束,使模型在单域条件下具备对未知目标域的良好泛化能力。下节将进一步分析式4中各组成部分的具体损失函数设计。首先,为增强源域与虚拟目标域之间的分布差异性,提出一种基于域差异感知的统计扰动机制。该机制不同于传统的样本级生成策略,而是在潜在空间中基于可学习的统计特征引入扰动。其次,为确保源域特征与其扰动版本之间的语义一致性,进一步引入一致性引导的协方差与语义对齐策略。该策略在缓解风格差异的同时,通过显式语义约束提升模型预测结果的一致性。由于源域与虚拟目标域的构造需满足式4所定义的结构,本文采用如图2所示的双流网络结构实现,该结构将在后文详细介绍。图2 通过差异性与一致性对抗机制实现的智能故障诊断 SDG 流程图图3 基于域差异感知的可学习统计扰动模块示意图3 所提方法3.1 基于域差异感知的可学习统计扰动在对抗式 SDG 框架(式4)基础上,本节首先设计用于域差异感知的可学习统计扰动模块(Learnable Statistical Perturbation, LSP)。其具体实现过程如图3所示。已有先前研究表明,通过扰动网络浅层中实例归一化后的仿射变换的均值与方差,可有效模拟域偏移,同时保留类别判别信息。设浅层特征表示为 ,其中 与 分别表示样本通道维度上的均值与标准差: 由于特征统计量编码了与域相关的信息,AdaIN(Adaptive Instance Normalization)方法通过将源域特征 的统计量替换为目标域特征 的统计量,以实现域风格迁移: 然而在 SDG 场景下,目标域特征统计量无法获取,成为建模过程中的关键难点。部分研究尝试通过线性插值或不确定性建模方式增强统计多样性,虽在一定程度上提升了泛化能力,但往往存在特征扰动空间探索不足、结构不明确等问题。为此,本文提出基于可学习统计偏移的 LSP 模块,生成统计可靠的伪目标域特征: 其中, 与 分别为可学习的均值与标准差偏移项。该模块克服了传统扰动策略的局限性,能够更有效地扩展原始特征与扰动特征之间的分布差异。鉴于对抗式 SDG 的关键目标在于最大化源域与虚拟目标域之间的分布差异,需引入合适的域感知距离度量引导模型学习合理的统计特征。先前研究表明,网络浅层特征主要捕捉风格信息,虽具有域相关性但与故障类别弱相关。受之前研究启发,深层协方差对齐(CORAL)可通过对齐源/目标域的二阶统计量降低域偏移。因此,本文利用协方差刻画网络浅层特征中的域差异性。整体网络由一系列顺序特征提取器 和一个分类器 构成,其结构可表示为: 其中 表示特征提取器的总层数。我们从浅层特征提取器中计算出一组协方差矩阵 ,并定义用于域差异感知的统计扰动损失为: 其中, 为特征向量维度, 表示矩阵的 Frobenius 范数平方, 为用于计算扰动损失的网络层数。 与 分别表示原始分支与扰动分支特征提取器, 为协方差矩阵计算操作,具体定义如下: 其中, 表示元素为1的列向量, 表示转置操作, 为样本数。3.2 一致性引导的协方差与语义对齐图4 一致性引导的协方差与语义对齐模块示意图延续对抗式 SDG 的建模思想(式4),另一项关键优化目标是在特征与语义层面对源域与虚拟目标域进行对齐。其具体实现过程如图4所示。为增强模型的一致性,我们从域特征与语义预测两个维度出发,对两个子网络进行对齐设计。在域级特征对齐方面,通过对浅层特征施加可学习统计扰动以模拟域偏移。然而,常用的对抗式学习或基于距离的对齐方法常存在训练不稳定或关键信息丢失的问题 。受实例白化(Instance Whitening)方法启发,去除对域变化敏感的协方差分量有助于提取域不变特征。因此,本文采用协方差对齐(Covariance Matching, CM)策略,在保留关键任务信息的同时,降低域偏移影响。具体地,分别计算原始特征与扰动特征的非中心协方差矩阵: 其中, 和 分别表示原始与扰动特征的非中心协方差矩阵,用于表征其二阶统计特性。通过最小化二者之间的 距离,实现协方差匹配以对齐两子网络: 尽管 CM 能在一定程度上对齐原始与扰动特征的整体分布,但未显式建模两者之间的直接关联。为此,本文进一步引入成对一致性损失(Pairwise Consistency, PC),通过交叉协方差矩阵刻画原始与扰动特征之间的一致性: 鉴于原始与扰动特征应保留相同的任务相关信息,理想情况下其交叉协方差矩阵应接近单位阵。因此,PC 损失通过约束该矩阵的对角元素接近1来保证特征表达的一致性: 其中, 表示提取交叉协方差矩阵的对角元素组成的列向量, 表示元素全为1的向量。在语义对齐方面,期望模型在嵌入空间中能够对具有相同语义标签的样本进行聚类,从而确保语义一致性。为此,本文采用最大化原始与扰动特征输出之间互信息(Mutual Information, MI)的方法,以保持语义表达一致性。互信息 定义为: 其下界可表示为: 其中, 为判别函数(critic function)。尽管直接最大化该下界有助于增强语义一致性,但可能抑制同类样本间的共享信息。为此,引入监督对比损失(Supervised Contrastive Loss, SCL),以显式增强同标签样本之间的互信息: 其中, , 表示与第 个样本类别一致的所有原始与扰动特征的输出集 合, 为温度系数。 此外,为进一步提升两子网络之间的语义一致性,最小化两分支网络的交叉熵损失: 3.3 训练与推理流程本文所提方法采用双分支网络结构,并通过交替训练的方式进行优化,训练过程分为两个阶段。在第一个阶段,执行一致性引导的协方差与语义对齐,以在源域与虚拟目标域特征之间保持全局与类别层面的一致性,从而在降低域间差异的同时保留任务相关信息。此外,引入显式语义约束以提升模型预测的稳定性与准确性。该阶段的损失函数定义如下: 其中, 、 和 分别为对应损失项的超参数权重。在第二个阶段,引导 LSP 模块在潜在特征空间中建模域间差异,生成更加多样且可靠的特征分布,以增强模型对域偏移的感知能力。其损失函数如下: 其中, 为域差异建模损失的超参数。最终,SDG 故障诊断通过上述两个阶段的交替对抗训练实现模型优化。在推理阶段,仅使用基于统计扰动特征训练得到的子网络进行预测,因为扰动特征有助于缓解域偏移对模型性能的影响。4 实验验证4.1 数据集介绍为全面评估所提 SDG 框架在故障诊断任务中的有效性,本文在五个具有代表性的旋转机械基准数据集上开展实验,覆盖多种运行工况与健康状态。数据集详情如下所示:帕德博恩大学(PU)数据集:该数据集包含 2100 条轴承振动样本,采集于三种运行工况(C1: 1500 RPM, 0.7 Nm, 400 N;C2: 1500 RPM, 0.1 Nm, 1000 N;C3: 900 RPM, 0.7 Nm, 1000 N)下的七种健康状态,每条样本包含 2048 个时域数据点,采样频率为 64 kHz。华中科技大学(HUST)数据集:该数据集包含 2100 条样本,覆盖三种运行工况(70 Hz、75 Hz、80 Hz)和七种健康状态(健康、内圈中/重度故障,外圈中/重度故障,滚动体中/重度故障),每条样本含有 2048 个时域数据点,采样频率为 25.6 kHz。北京工业大学(BJUT)数据集:该数据集采集自行星齿轮箱,在三种转速(1800、2100、2400 RPM)下记录了 3000 条样本,涵盖损齿、齿面磨损、齿根断裂和缺齿等四种故障模式,每种状态包含 200 条样本,采样频率为 48 kHz。山东科技大学(SDUST)数据集:该数据集包含 6000 条样本,覆盖三类故障位置、三个故障程度,共 10 种健康状态,采样于 2000、2500 和 3000 RPM 三种转速下。每条样本均包含 2048 个时域数据点,采样频率为 12.8 kHz。真实工厂轴承(RFB)数据集:该数据集来自工业生产线,包含 2400 条真实故障轴承样本,涉及四种状态(健康、滚动体缺失、纱卡、滚针偏移)与三种转速(100、200、300 RPM)。采样频率为 20.48 kHz,每条样本包含 2048 个时域数据点。4.2 实验设置4.2.1 参数配置所用特征提取器为五层卷积神经网络(CNN),后接两层全连接分类器。模型的学习率设为 0.0005,迭代轮数为 300,批次大小为 40。输入信号经快速傅里叶变换(FFT)转化为频域表示,未额外进行去噪或归一化处理。LSP 模块部署于 CNN 的前两层,其超参数设置为 -1; 和 部署于第三层,对应超参数分别为 10 和 1; 部署于分类器输出端,超参数设为 0.1。4.2.2 对比方法本文方法与多个现有主流 SDG 方法进行了性能对比,具体包括:基线模型:使用本文所构建的 CNN 结构,但不引入任何泛化策略。 先进 SDG 方法:包括 MEADA 、L2D、AMINet、MSG-ACN、ACL和 DEFSDG。4.2.3 评估指标在 SDG 故障诊断任务中,从一个运行工况中选取数据作为源域用于训练,其余工况数据作为目标域用于测试。为全面评估模型性能,采用分类准确率(Accuracy)、宏平均召回率(Macro Recall)和宏平均 F1 分数(Macro F1 Score)作为评估指标。准确率定义如下: 其中,正确预测样本数为模型预测结果与真实标签一致的样本数量,分母为测试集中样本总数。宏平均召回率表示对各类召回率求平均,具体定义为: 其中, 表示类别总数, 与 分别为第 类的真正例与假负例数。宏平均 F1 分数结合了精度与召回率,定义如下: 其中, 和 分别表示第 类的精度与召回率。4.3 结果与讨论4.3.1 主要结果分析表1至表5展示了在五个基准数据集上的 15 个领域泛化故障诊断任务中,本文方法与现有典型 SDG 方法的诊断性能对比结果。表头中的“S”表示每个泛化任务所使用的源域。所有任务均重复运行五次,并报告主要性能指标的平均值与标准差。 首先,可以观察到,基线 CNN 模型在所有任务中性能最差,表现出较弱的泛化能力。其次,与基线模型相比,多个先进的 SDG 方法在诊断性能上均取得显著提升,这主要得益于它们在训练中引入了数据增强、特征多样化策略,并融合了域泛化机制,从而增强了模型的特征表示能力。最后,在大多数实验设置下,本文方法在准确率、宏平均召回率和宏平均 F1 分数三项指标上均显著优于现有方法,且表现出更小的标准差,充分验证了所提方法在工业诊断任务中的鲁棒性与泛化能力。 为直观验证本文方法在缓解域偏差方面的有效性,我们在 HUST 数据集上设计了一个泛化任务(源域:70 Hz,目标域:75 Hz)。将不同方法学习到的源域与目标域特征表示通过 t-SNE 降维至一维后,采用核密度估计(KDE)进行可视化,所得密度曲线可用于评估两个域特征分布的重叠程度,从而判断其相似性。结果如图5所示。 由图可见,基线模型所学习的源域与目标域特征之间存在显著分布差异,主要由于其未能有效缓解源目标域之间的数据分布不一致问题。其他典型 SDG 方法(如 L2D、AMINet、MSG-ACN、DEFSDG)通过特征生成和域不变特征学习在一定程度上缓解了分布偏差。而本文方法通过交替执行扰动特征生成与域不变特征学习,更有效地优化了域对齐过程,最终在分布重叠程度上优于其他方法,表现出更强的特征泛化能力。 表1 不同方法在PU 数据集泛化任务中的结果对比表2 不同方法在HUST 数据集泛化任务中的结果对比 图5 不同方法在源域与目标域中的特征分布可视化:(a) CNN,(b) L2D, (c) AMINet,(d) MSG-ACN,(e) DEFSDG,(f) 本文方法 表3 不同方法在BJUT 数据集泛化任务中的结果对比表4 不同方法在SDUST 数据集泛化任务中的结果对比表5 不同方法在RFB 数据集泛化任务中的结果对比 4.3.2 消融实验分析为分析所提方法中各个模块的作用,表6给出了在五个数据集上不同消融策略下的诊断准确率(%)对比结果。表中 Method 1 至 Method 4 分别对应具体的模块组合策略。 从表中可观察到:与基线模型相比,Method 1 引入 LSP 模块后,模型能够探索更广泛的潜在域空间,显著提升了其泛化能力。Method 1、2 和 4 的对比表明,无论是增强协方差对齐,还是语义对齐,都能明显提升泛化准确率。此外,Method 3 与 Method 4 的对比结果表明,将协方差匹配与成对一致性联合使用,可有效抑制风格特征干扰,进一步改善模型的泛化性能。 最终,本文方法在整合全部模块的基础上,在五个数据集上均取得最佳性能,说明各模块协同优化对于提升故障诊断准确性具有重要作用。 表6 不同方法在各类数据集泛化任务中的结果对比为进一步验证各模块策略在 SDG 故障诊断中的有效性,本文在 HUST 数据集上进行泛化任务,并采用 t-SNE 将不同消融方法所得源/目标域特征降维为二维后进行可视化,结果如图6所示。图中使用不同颜色表示不同健康状态,不同形状表示不同域,以便于理解。图6 在HUST 数据集上不同消融设置下的特征表示二维可视化:(a) 基础 模型,(b) M1,(c) M2,(d) 本文方法图7 在不同协方差对齐策略下双流网络学习的特征协方差矩阵对比 首先,如图6(a) 所示,在基线模型下,源域与目标域特征高度重叠,类别边界模糊。引入 LSP 策略后(图6(b)),源/目标域特征间分布距离明显扩大,表明该策略能够有效引入域差异性。进一步引入协方差对齐(图6(c))后,源/目标域特征间距离缩小,域偏差得到缓解。当加入语义对齐策略后(图6(d)),同一类别的特征在源域与目标域中形成更加紧密的聚类。 为验证所提协方差对齐策略的有效性,图7展示了在不同对齐策略下,双流网络学习到的特征协方差矩阵对比。基线网络中的协方差矩阵差异明显,对角特征几乎无法分辨;引入协方差匹配策略后,矩阵整体更趋一致,但对角信息仍不明显,表明虽然实现了原始与扰动特征的整体对齐,但存在信息丢失风险。引入成对一致性后,对角特征增强,保留了部分任务信息,但矩阵仍存在不一致;联合两种策略后,不仅实现了矩阵结构的一致性,还保留了清晰的对角线特征,实现了原始与扰动特征在保留关键信息基础上的有效对齐。图8 模型第一层中原始与可学习扰动特征统计量在不同训练迭代下的对比图9 模型第一层中原始与扰动特征的可视化对比:(a)10 次迭代,(b)30 次迭代,(c)60 次迭代图10 基线模型与本文方法在不同数据集上的混淆矩阵对比:(a) HUST,(b) SDUST,(c) PU,(d) RFB,(e) BJUT 为直观展示模型如何学习统计特征的变化,图8展示了第一层卷积中 16 个通道在不同训练阶段的统计量变化。图中蓝色条表示原始特征的均值与方差,黄色条表示扰动特征。在第 10 次迭代时,两者差异较小;到第 30 次迭代,部分通道出现明显偏移;第 60 次迭代时,多个通道在均值与方差上差异显著,说明扰动特征逐步扩大了统计空间的分布差异。 进一步,图9利用 t-SNE 对上述原始与扰动特征进行可视化。图中不同颜色表示不同类别,不同形状区分原始与扰动特征。从结果可见,初期(10 次迭代)原始与扰动特征基本重合;随着训练推进,它们之间的分布距离逐渐扩大,体现出模型构造的域差异性增强。 图10展示了基线模型与本文方法在五个数据集上的混淆矩阵对比。可以看出,本文方法在所有数据集上均显著降低了误分类情况,提升了故障可分性。其中,HUST 与 SDUST 数据集中的分类混淆显著减少;而在 PU 数据集中,第四类故障的识别仍具挑战性,但本文方法在该类上也优于基线模型;在 RFB 与 BJUT 数据集上,本文方法相较基线模型展现出更强的泛化能力与诊断鲁棒性。4.3.3 计算开销评估表7 不同方法在模型复杂度与推理效率方面的对比表7对比分析了不同代表性方法在模型复杂度与推理效率方面的差异。评估指标包括参数量(Params)、浮点运算次数(FLOPs)、推理时间、内存占用与单样本能耗,旨在反映模型在实际部署中处理每个样本的资源开销。所有指标均在 Jetson Nano 平台上进行测量,推理批次大小设为 1,以模拟典型边缘设备的应用场景。内存使用表示模型在单样本推理过程中的运行时内存占用,单样本功耗为处理一个样本所需的平均能量。 从表中可见,本文方法在保持准确性的同时具备良好的计算效率,参数量为 16.34M,FLOPs 为 21.18M,平均推理时间仅为 31.27 ± 0.72 毫秒,内存占用为 73.24 MB,单位样本功耗为 0.15 W,展现出极低的部署开销。与 ACL 与 DEFSDG 等多分支设计方法相比,本文方法在推理成本方面具有显著优势。 值得注意的是,MEADA、L2D、AMINet 及 MSG-ACN 等方法虽在训练阶段引入样本生成策略,但其推理阶段结构与本文方法保持一致。因此,本文在表7中省略了这些方法,以避免冗余对比。结论本文针对工业设备的故障诊断任务,提出了一种基于单域泛化(SDG)的诊断方法,利用差异性与一致性的对抗建模机制提升模型在复杂工况下的泛化能力。该方法在浅层特征空间引入可学习统计扰动(LSP),主动模拟由负载、转速等运行条件变化引起的域偏移,从而增强模型对环境变化的鲁棒性与适应性。此外,结合协方差对齐与语义一致性对齐策略,从全局与类别层面对原始特征与扰动特征进行对齐,有效提升了模型在多变工业场景下的预测稳定性与一致性。为验证所提方法的实际效果,本文在五个典型旋转机械数据集上开展了系统实验,主要结论如下:所提方法在准确率、宏平均召回率和宏平均 F1 分数等指标上显著优于现有主流 SDG 方法,验证了其在工业相关泛化场景下的诊断可靠性;LSP 模块在潜在特征空间中显式模拟统计扰动,使模型能够更有效地捕捉复杂工况下可能出现的域偏移。与传统的随机扰动方式相比,LSP 提供了更具可控性和泛化性的分布建模能力,展现出良好的工业应用潜力;协方差对齐机制有效抑制风格特征干扰,语义对齐策略进一步增强了故障类别识别的稳定性,使该方法尤为适合部署于高动态变化的工业现场环境。需要指出的是,本文方法假设训练域与目标域在潜在空间中存在共享结构,在面对严重分布偏移时仍存在泛化失效的风险。为缓解该问题,可结合数字孪生技术与高保真仿真生成物理先验增强数据,从而提升训练数据多样性与模型适应性。此外,本文方法主要基于浅层网络中的统计特征建模,对于复杂域偏移的表达能力仍有限。未来可进一步探索高阶统计量、协方差结构或深度生成模型,以提升模型对复杂领域间差异的感知能力和跨域泛化性能。编辑:曹希铭校核:李正平、陈凯歌、赵栓栓、赵学功、白亮、任超、Tina、陈宇航、海洋、陈莹洁、王金、赵诚、肖鑫鑫该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈