首页/文章/ 详情

基于多尺度时间记忆Transformer框架的工业机器人退化感知剩余寿命预测

11月前浏览1952
     本期关注工业机器人性能退化预测的新进展:剩余使用寿命(remaining useful life, RUL)预测对于保障工业机器人的稳定运行至关重要。然而,现有深度学习方法在感知由健康状态向退化状态的转变方面存在局限,同时难以有效建模跨尺度的时序依赖关系,限制了预测精度。为破解这一难题,本文提出了一种多尺度时序记忆Transformer框架(multiscale temporal memory Transformer, MTMT):通过引入记忆自编码器Transformer网络实现状态转变点的识别,并结合多尺度时序Transformer网络(multiscale temporal Transformer network, MTTNet)挖掘局部与全局特征,从而实现跨阶段的高精度RUL预测。基于自主搭建的工业机器人加速寿命实验平台的验证结果表明,该方法不仅能够提前定位状态转变点,还在预测精度方面显著优于多种先进方法,为工业机器人健康管理与智能运维提供了一条全新的解决方案。  

    论文基本信息

    论文题目:

    Degradation-Aware Remaining Useful Life Prediction of Industrial Robot  

    via Multiscale Temporal Memory Transformer Framework  

    论文期刊:Reliability Engineering and System Safety

    论文日期:2025年

    论文链接:

    https://doi.org/10.1016/j.ress.2025.111176

    作者:Zhan Gao (a), Chengjie Wang (a), Jun Wu (a*),Yuanhang Wang (b), Weixiong Jiang (a), Tianjiao Dai (a)

    机构:

    a: School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan, China;  

    b: Sino-German College of Intelligent Manufacturing, Shenzhen Technology University, Shenzhen, China

    团队带头人简介:吴军老师现任华中科技大学船舶与海洋工程学院教授、博士生导师,副院长。作为IEEE高级会员、斯坦福大学访问学者,他长期致力于复杂机电装备的智能运维与健康管理研究,重点涵盖健康监测、异常检测、状态评估、故障诊断、寿命预测及维修决策等方向,并积极推动人工智能、大数据、物联网、数字孪生等新兴技术与装备工程的深度融合。在科研方面,吴教授主持国家自然科学基金项目3项、国家重点研发计划课题1项及多项省部级与企业委托课题,发表学术论文160余篇(其中JCR Q1区论文70余篇,ESI高被引论文8篇),授权专利与软件著作权60余件,参与制定国家标准5项,多项成果获省部级科技进步奖。其学术影响力连续入选斯坦福大学全球前2%顶尖科学家“年度影响力”榜单。(来源http://faculty.hust.edu.cn/Junwu/

    目录

    1 摘要

    2 引言

    理论基础

    3.1 逆离散傅里叶变换

    3.2 多尺度分块操作

    多尺度时序记忆Transformer

    4.1 特征提取与选择

    4.2 基于MATNet的状态变化识别

    4.3 基于MTTNet的RUL预测

    实验研究

    5.1 实验设置

    5.2 实验结果

    5.3 剩余寿命预测结果与分析

    6 讨论

    6.1 消融实验

    6.2 参数分析  

    6.3 与其他先进方法的对比

    7 结论

    1 摘要

    剩余使用寿命(RUL)的预测对于保障工业机器人(IR)的稳定运行至关重要。近年来,深度学习方法已在RUL预测任务中展现出良好性能,但其在感知IR由健康状态向退化状态转变的能力有限,且难以有效建模RUL预测所需的多重时序依赖关系。为解决这一问题,本文提出了一种多尺度时序记忆Transformer框架,将寿命预测与状态转变识别相融合。该框架由记忆自编码器Transformer网络与多尺度时序Transformer网络组成:前者用于捕捉时序信号中的潜在变化并检测状态转变点,后者则在此基础上挖掘跨尺度的时序依赖以实现精准的寿命预测。本文进一步构建了自主研发的IR实验平台对所提方法进行验证。实验结果表明,相较于其他先进方法,该框架不仅能够提前识别状态转变点,而且显著提升了IR寿命预测的精度。

    关键词:剩余使用寿命,状态变化识别,Transformer模型,多尺度时序特征,工业机器人

    2 引言

    工业机器人(IR)因技术进步和高产能需求在过去二十年中得到广泛应用。然而,复杂的工作环境导致其不可避免地出现损伤,因此需引入预测与健康管理(PHM)技术进行状态监测和健康维护。其中,剩余使用寿命(RUL)预测可有效防止设备退化并提升可靠性。

    随着传感技术的发展,传感器可获取大量健康状态数据,推动了数据驱动的 RUL 预测研究。深度学习(DL)方法因强大的非线性建模能力在工业机器人的 RUL 预测中表现优异,但仍难以捕获全局时间依赖,导致预测存在不确定性与滞后。Transformer 网络通过多头注意力机制可高效捕获全局依赖,已被用于增强 RUL 预测。然而,现有方法仍存在两大问题:一是 IR 在初期运行稳定,退化在不确定状态变化点开始,传统 Transformer 假设退化从初始服役开始,难以提前感知状态变化;二是 IR 数据为长寿命周期序列,包含多期(长、中、短期)潜在退化模式,现有方法难以捕获多期特征,造成时间信息丢失与预测精度下降

    针对上述问题,本文提出多尺度时序记忆 Transformer 框架实现退化感知 RUL 预测。框架包括两部分:首先,记忆自编码器 Transformer 用于信号重建并识别状态变化点;其次,多尺度时序 Transformer 捕获多期依赖以增强状态变化点后的 RUL 预测能力。为验证方法有效性,搭建加速寿命测试 IR 平台。主要贡献包括:

    • 提出退化感知 RUL 预测方案,通过多尺度时序记忆 Transformer 实现 IR 的高精度 RUL 预测,并可及时识别健康与退化状态。    

    • 设计记忆自编码器 Transformer 网络,挖掘动态时序依赖,实现状态变化点的提前检测。    

    • 构建多尺度时序 Transformer 网络,可选择预测路径并捕获多期依赖,实现无完全失效数据下的未来 RUL 预测。

    3 理论基础

    3.1 逆离散傅里叶变换
       

    离散傅里叶变换(DFT)及其逆变换(IDFT)共同组成一组可逆的变换对。通过 DFT,可将时域信号            映射为频域信号           ,其数学表达式为:

             

    其中,           表示样本总数,                     表示频率索引,           表示时间索引。符号            与            分别表示            的实部与虚部。其幅值与相位信息表示为:

             

    其中,           表示幅值,           表示相位。IDFT 能够有效地对输入序列进行频谱分析,从而揭示各频率所对应的幅值与相位信息,并将频域信号            转换回时域序列           。其表达式为:

             

    IDFT能够通过频率分量实现信号重构,从而支持数据压缩、噪声抑制,并保留关键频率以供分析。        

    3.2 多尺度切片操作        

    多尺度切片操作用于将输入数据或图像划分为固定大小的片段(patch),并将这些片段重构为向量序列,以作为 Transformer 网络的输入。切片操作的具体细节如图 1 所示。        

         

     1 多尺度切片操作的细节        

    对于输入序列           ,其中            表示时间维度,           表示特征维度,该序列可被划分为            个大小为            的片段,进而得到向量                      与            之间的关系表示为:

             

    为了获得多尺度切片集 合,首先基于            个不同的切片尺寸值定义切片尺寸集 合           。对于切片尺寸           ,可得到对应的切片           。进一步地,通过多尺度切片操作可得到切片集 合           。该切片操作能够从多个视角表征数据,从而在不同时间分辨率下增强时间特征表示。

    多尺度时序记忆Transformer

    如图 2 所示,提出了一种多尺度时间记忆 Transformer 框架,用于实现面向退化感知的 IR 剩余寿命(RUL)预测任务。该框架包含三个部分:特征提取与选择、记忆自编码 Transformer 网络(MATNet),以及多尺度时间 Transformer 网络(MTTNet)

     

     2 多尺度时间记忆 Transformer 框架流程图

    4.1 特征提取与选择    

    特征提取与选择模块旨在从测量信号中提取退化相关特征。提取的特征包括 7 维时域特征与维频域特征,并通过欧氏距离对特征进行选择。维时域特征为:均值幅值(tf1)、均方根(tf2)、峰值(tf3)、平方均根(tf4)、偏度(tf5)、峰度(tf6)以及标准差(tf7);5维频域特征为:平均频率(ff1)、频率方差(ff2)、频率中心(ff3)、一阶矩(ff4)及频率偏度(ff5)。由于部分特征对IR 退化现象的表征能力有限,同时高维特征增加了计算开销,因此采用欧氏距离方法对特征进行筛选。

    4.2 基于 MATNet 的状态变化识别    

    在正常状态向退化状态过渡过程中,可能出现微弱且受噪声干扰的异常信号。为准确捕捉这些细微变化并定位状态变化点,本文采用 MATNet。该网络由三部分组成:Transformer 编码器、增强记忆模块和密集解码器,其整体结构如图 3 所示。

     

     3 MATNet 结构示意图    

    4.2.1 Transformer 编码器    

    Transformer 编码器用于将测量信号嵌入到潜在空间中。设        为测量信号,其中        表示输入维度,       表示时间步数。Transformer 编码器采用标准 Transformer 架构,将测量信号        编码为长期状态。在编码器中,通过位置编码对测量信号        添加位置信息。这些信号被转换为查询向量       、键向量        和值向量       ,并输入多头注意力机制以捕捉长期时间依赖关系。

    多头注意力(Multi-head Attention)的表达式为:

         

    其中,       为变换维度,       表示可训练矩阵,       为输出维度。多头注意力的输出随后被输入前馈层,以完成维度变换并获得隐藏状态       

         

    4.2.2 Dense 解码器    

    Dense解码器用于根据隐藏状态进行信号重构。它采用两层线性层将潜在空间映射回输入空间:

         

    其中,       和        为可学习参数,       与        为偏置项。Transformer 编码器与密集解码器之间的重构损失表示为:

           

    4.2.3 增强记忆模块

    传统自编码器架构通常采用较小的重构损失来训练模型,但由于过度泛化及时间信息丢失,其重构能力有限。为解决该问题,引入了增强记忆模块,通过整合记忆机制来提升性能。该机制能够基于存储在记忆模块中的典型正常模式提取隐藏表示,同时限制编码器捕捉异常特征的能力。记忆机制包含两个阶段,即记忆更新阶段和查询更新阶段。

    在记忆更新阶段,记忆项          (        ) 被训练以存储输入信号中的典型正常模式,其中          表示记忆项编号。采用增量更新的方法对记忆项进行更新,这些记忆项包含在时间步          的正常模式下查询          的原型。该更新通过构建基于查询条件的记忆注意力          (        ) 实现:

           

    其中,         为可训练的温度参数,用于调节分布的集中程度。更新门          用于训练正常模式下的记忆项。它可控制通过查询获得的新模式在多大程度上被整合到当前记忆项的典型正常状态中。更新过程表示为:

           

           

    其中,         和          为可训练的变换矩阵,         与          分别表示 Sigmoid 函数和乘法运算。记忆更新阶段仅在训练阶段进行。

    在查询更新阶段,生成更新后的查询并将其输入密集解码器。采用 Softmax 函数对查询与记忆项进行计算,从而获得基于记忆的查询注意力:

             

    此外,通过将记忆项            与基于记忆的查询注意力相结合,可进一步获得检索后的记忆查询           

               

    进一步地,将              与              进行拼接,可得到更新后的查询             ,作为密集解码器的新输入。同时,潜在空间与输入空间被用于构建基于偏差的状态变化识别标准。潜在空间偏差(DLS)定义为查询与其在时间步              下最近记忆项              的距离。由于记忆项存储了正常状态的原型,因此退化状态下的 DLS 通常大于正常状态下的 DLS。输入空间偏差(DIS)定义为时间步              下输入信号与重构信号的距离。DIS 与 DLS 相乘,用于放大正常状态与退化状态之间的差异。状态评分              可表示为:

               

               

               

    退化状态下的状态评分高于正常状态下的评分。随后,通过设定阈值              来判定样本点的当前状态:

               

    当              个连续样本点的状态评分均超过阈值时,第一个样本点对应的时间点被识别为状态变化点。

    4.3 基于 MTTNet 的剩余寿命预测          

    为在确定的状态变化点之后实现 RUL 预测任务,MTTNet 被用于揭示IR 的退化行为。如图 4 所示,MTTNet 包含切片构建模块、多尺度注意力模块以及路径融合模块。

           

     4 MTTNet 结构示意图          

    4.3.1 切片构建模块           

    输入序列              首先被划分为              个不同切片大小              的片段,从而得到切片集 合             。由于序列具有特定的时间动态,不同序列可能适合不同尺度,而某些尺度可能引入冗余信息。为了捕捉多尺度时间特征,构建了切片构建模块。该模块用于构建路径集 合,并为多尺度注意力模块选择最优切片大小。

    在切片构建模块中,采用时间解缠(temporal disentanglement)方法从输入序列中获取周期模式和趋势模式。时间解缠包括周期分解和趋势分解,如图 5 所示。

           

     5 时间解缠结构示意图          

    周期分解利用离散傅里叶变换(DFT)将输入序列              分解为傅里叶基             

               

               

               

    其中,             和              分别表示幅值与相位,                         与              分别表示实部与虚部。为了保持频域的稀疏性并去除冗余信息,仅保留傅里叶基              的前              个频率,其中              设置为 3。随后,采用逆离散傅里叶变换(IDFT)挖掘频域中隐藏的周期特征,其构建基于              个傅里叶频率:

                 

    趋势分解用于区分周期性与非周期性成分,从而获得周期成分中潜在的趋势模式:

                 

    为提升趋势特征的表征能力,对                采用多种核尺寸的平均池化操作,以提取不同的趋势信息。利用 Softmax 激活函数对各趋势特征进行加权,得到高阶趋势特征               ,表示为:

                 

    其中,               表示采用第                个核的池化操作,               表示核的数量。此外,路径函数                用于生成路径权重,可用于确定应选择的切片大小。为避免重复选择某些尺度而忽略其他尺度,引入高斯噪声项以保持生成路径权重的稳定性:

                 

    其中,               和                表示可训练参数,               和                分别表示维度与切片数量。               表示 Mish 激活函数。为了保持路径的稀疏性,采用 top-               操作,仅保留前                条路径,其余路径权重置为 0。

    4.3.2 多尺度注意力模块             

    多尺度注意力模块用于从切片集 合中捕捉全局与局部时间依赖特性。该模块采用两种视角,即切片内视角(in-patch view)与切片间视角(patch-wise view),以同时挖掘全局和局部特征。在切片间视角中,建立切片内部不同时间步之间的关系,以提取局部特征。切片集 合中的每个切片                被嵌入为切片间表示               ,其中                表示嵌入维度。对切片间表示进行线性变换,以得到键向量               和值向量               

                   

    其中,                 和                  表示线性变换。定义初始化查询向量                  用于计算切片间注意力:

                   

    将所有切片的注意力结果整合,以形成切片间视角(patch-wise view)的最终输出:              

                   

    切片内视角用于建立切片之间的关系,从而提取全局特征。首先采用嵌入操作对每个切片进行映射,并重新排列以形成切片内表示                 ,其中                 。随后,对切片间表示进行线性变换,以获取键向量                  和值向量                 

                   

    其中,                 和                  表示线性变换,切片内注意力(in-patch attention)可以计算为:

                   

    其中,                 表示可训练向量。将所有切片的注意力输出整合,以形成切片内视角的最终输出:

                     

    双注意力机制的输出                    表示为:

                     

    其中,                   用于对                    进行维度变换,使其与                    的维度保持一致。

    4.3.3 路径融合模块                  

    路径融合模块用于整合多尺度补丁的特征。                     的维度与多尺度注意力模块的补丁大小相关。若 ,表示采用了对应补丁大小                      的多尺度注意力模块;若                     ,则表示忽略了补丁大小                      的多尺度注意力模块。假设在多尺度注意力模块中补丁大小                      的输出定义为                     ,则通过加权聚合对这些输出进行整合。不同补丁集的输出在融合后用于 RUL 预测:

    其中, 表示指示函数。当  时,指示函数的输出为 1;而当                      时,指示函数的输出为 0。

    5 实验研究

    5.1 实验设置

    为验证所提方法,构建了一套自研工业机器人(IR)平台。该平台由 6 个关节轴组成,用于驱动电机、减速器及传动系统,额定负载为 20 kg。在加速寿命测试中,IR 执行重复抓取铁块的任务,运行负载设定为额定负载的 1.5 倍,动作包括垂直抓取、320° 横向平移及垂直放置。每个轴配备振动、速度、位置和力矩传感器采集信号,每个采样周期包含 3000 个采样点,数据每 3 小时记录一次。测试过程中,第五轴与第六轴减速器在 1480 h 后因机械磨损出现噪声及润滑脂泄漏。平台结构示意如图 6 所示。

     6 自建 IR 平台的结构细节

    各传感器提取的 RMS 值如图 所示。速度、振动及位置信号未显示明显退化,主要因各关节采用闭环控制。第三轴与第四轴的力矩信号呈轻微退化,第五轴与第六轴表现明显退化。其中,第六轴减速器在 2325 h 发生故障,1480 h被定义为实际状态变化点。本研究以第六轴为 RUL 预测对象。  

    图 不同传感器提取的 RMS   

    5.2 实验结果  

    采集数据经归一化与移动平均平滑预处理后,对测量信号进行特征提取与选择。每个关节轴提取 7 维时域特征及维频域特征(图 8)。欧几里得距离阈值设为 0.1,最终选取转矩传感器的九个特征(tf3-tff4-ttf1-ttf6-ttf5-tff2-ttf1-tff3-ttf2-t),作为状态变化识别及IR 剩余寿命预测的输入。

    图 基于欧几里得距离的特征选择  

    为了实现 RUL 预测,首先训练MATNet 进行状态变化检测,其次训练 MTTNet 用于IR  RUL 预测。所提出 Transformer 框架的超参数列于表1。这些参数通过网格搜索确定,并采用 Adam 优化器,学习率分别为0.0005MATNet MTTNet)。训练过程中使用早停策略以防止过拟合。

    表 1 Transformer 框架的超参数设置

    为评价检测性能,引入伪劣化(FD)与伪正常(FN)两类指标。FD定义为预测的状态变化点之前的错误劣化标签,FN 定义为预测的状态变化点之后的错误正常标签。前 200 个样本用于训练 MATNet,其余样本用于测试。样本点的 nf 值设定为 3。状态变化检测结果如图 所示。

    图 9 MATNet 的检测结果

     N 分别表示劣化状态与正常状态,红色与绿色 区域分别表示正常阶段与劣化阶段。检测出的状态变化点为 488FD数量为 1FN 数量为0。实验结果表明,本文提出的 MATNet 能够有效检测状态变化,并能够提前预警劣化状态。为定量评估方法的预测性能,采用三种评价指标:归一化均方根误差(NRMSE)、得分(Score)与准确率(AC)。

       

       

    其中,     表示预测的 RUL 值,     表示预测 RUL 值与实际 RUL 值之间的误差,     为样本点总数。

    在已知样本点范围          内对 MTTNet 的 RUL 预测性能进行验证。已知样本点用于训练,其余样本点用于测试。图 10 表明,预测性能与已知样本点数量密切相关。随着已知样本点数量的增加,预测能力逐渐提升。当已知样本点数为      时,预测准确率达到 95.30%。

    图 10 已知样本点下的 RUL 预测性能  

    为进一步验证 MTTNet 的预测性能,分别选取退化状态下 90%60% 30% 的样本点进行短期、中期及长期预测,同时选取 20% 样本用于超长期预测。不同时间尺度下的不确定性量化 RUL 预测结果如图 11 所示,其中预测值取50% 分位数(QL),10% QL 90% QL 之间的区间以灰色标记。结果表明,预测 RUL 在初期及中期退化阶段存在波动,而末期退化阶段逐渐收敛;长期与超长期预测的不确定性区间较宽,这是训练样本有限所致。该结果表明 MTTNet 能有效刻画 IR 的退化趋势,为维护计划设计提供依据

    图 11 不同时间尺度下的不确定性量化 RUL 预测结果  

    6 讨论

    本节基于消融分析、参数分析以及与其他先进方法的对比,对所提出方法的性能进行了系统评估。在参数分析中,重点考察了补丁大小、Top-K 操作以及阈值      对模型性能的影响。此外,通过与现有先进方法的对比,进一步验证了所提方法的优越性。

    6.1 消融实验

    针对 MTTNet RUL 预测性能,前 200 个样本用于训练,其余样本用于测试。五种方法用于消融对比:M1 为完整MTTNetM2 去除状态变化检测;M3 去除补丁构建模块;M4去除 Patch-wise 视图;M5去除 In-patch 视图。结果表明(图12),完整方法能够实现最优 RUL 预测性能,说明状态变化检测、多尺度补丁构建及双视图注意力机制对预测精度均有显著提升。

    图 12 不同方法消融实验的预测结果  

     M2 相比,M1 NRMSE 与 Score 值更低,说明状态变化检测对于识别初始退化点具有重要作用。结合状态变化检测,NTNNet RUL 预测性能得到了显著提升。M3 的结果表明,多尺度操作能够使模型挖掘多时间尺度的特征,从而增强预测性能。当去除 Patch-wise 视图时,预测能力有所下降,这是因为 Patch-wise 视图提供了退化过程的局部时间信息。而去除 In-patch 视图则会导致全局时间信息的丢失。M4 的性能低于 M5,这可能是由于全局时间信息包含比局部时间信息更丰富的时间特征。综上所述,结合状态变化检测可以显著提升 RUL 预测性能;多尺度划分与多尺度注意力模块能够有效挖掘时间信息,从而增强预测能力。

    6.2 参数分析  

    6.2.1 补丁大小  

    为研究补丁大小及补丁数量对模型性能的影响,构建了五种不同的补丁集进行对比实验(表 2)。实验中,前200 个样本点用于训练所提出的 MATNet,结果取十次运行的NRMSE 平均值。需要注意的是,各补丁大小需能够整除时间窗口值。随着补丁大小的增加,补丁集的尺度层次变得更加丰富,从而为多尺度特征提取提供更多选择。

    表 不同补丁集的构建方案

     

    如图 13 所示,随着补丁数量的增加,NRMSE 呈下降趋势。然而,当补丁大小进一步增大时,预测性能提升幅度逐渐减小。在超长期预测中,P3 NRMSE 小于 P4,而在其他时间尺度下,P4 NRMSE 又高于 P3。此外,P5 NRMSE 超过 P4,可能因补丁数量过多导致模型过拟合。因此,选择补丁大小为 {1, 2, 4, 8} 的补丁集可获得最佳预测性能

     

    图 13 不同补丁集的 NRMSE 结果箱线图    

    6.2.2 Top-K 操作    

    Top-K 操作用于调控预测过程中补丁集的稀疏性。图 14 显示了不同值下的 NRMSE Score 结果。随着 K 值增加,NRMSE 改善幅度逐渐减小,而Score 并未呈现相同趋势。当 K 值超过 4 时,Score反而下降,这可能因引入过多多尺度时间信息而导致模型过拟合。综合考虑稀疏性与预测性能,K = 4 可实现最优的RUL 预测效果    

     

    图 14 基于不同 Top-K 参数的NRMSE  Score 指标对比

    6.2.3 阈值δ    

    为了研究不同阈值对检测性能的影响,选择了八个不同的阈值以确定最优阈值 δ。利用FNFD 以及检测位置指标评估 MATNet 的检测性能。从图 15 可以看出,随着阈值的增大,FN  FD 的数量逐渐减少。同时,随着阈值增大,检测到的状态变化位置呈现延迟趋势。综合不同指标的检测结果,当阈值 δ 设置为0.9 时,MATNet 的检测性能表现最佳

     

     15 不同阈值下的检测性能    

    6.3 与其他先进方法的对比    

    为进一步验证所提出方法的优越性,本文对 MATNet 的检测性能与MTTNet 的预测性能分别与多种先进方法进行了对比。为评估 MATNet 的状态变化检测能力,选取了七种先进方法进行比较:Autoencoder (AE):编码器基于Transformer,解码器为全连接层;RVAE:变分自编码器(VAE),编码器与解码器均基于TransformerSTOC:改进自编码器,由 Transformer 编码器与卷积解码器组成;FCVAE:能够提取局部与全局时间信息进行异常检测的 VAE 方法;SSPCL:基于对比学习的状态变化检测方法;RSTRN:基于关系网络,通过评估健康状态与故障状态的相似性实现检测;CWTKB:考虑零故障样本的智能故障检测方法。    

    检测性能以状态变化检测位置(Location)、FD、FN以及平均训练时间(ATT, Average Training Time为指标进行评估。状态变化识别结果如表 3 所示。

     3 与其他先进检测方法的比较    

     

    与其他先进检测方法相比,MATNet 能够更早地识别状态变化,其 FD 与 FN 总数最少。相较于 AE 方法,MATNet 借助增强记忆模块具备较强的时间依赖捕捉能力,可揭示时间信息的动态变化以实现精确检测。此外,MATNet 在检测性能上优于其他改进型自编码器方法,显示出相较于现有先进方法的显著优势。然而,由于增强记忆模块结构复杂,MATNet 平均训练时间(ATT高于部分其他方法。

    为进一步验证 MTTNet 的预测性能,选择了八种先进预测方法进行对比,具体如下:BiGRU-TSAM:基于双向GRU 与自注意力机制的 RUL 预测方法。CoT:结合多尺度卷积与Transformer 解码器的改进型 Transformer 网络。MTST:多分辨率Transformer,用于长期时间序列预测。TFT:基于多头稀疏注意力的改进型 Transformer,用于 RUL 预测。ATCN:基于自注意力的TCN,用于 RUL 预测。SCTCN:自校准模块增强的TCN,用于 RUL 预测。MR-LSTM:多尺度LSTM 网络,用于 RUL 预测。FF:增强频率信息的Transformer,用于时间序列预测。

    选取退化状态下前 150 个样本点用于训练MTTNet,其余样本用于测试训练模型。图 16 展示了预测结果。

     

    16 各方法RUL 预测结果对比

    结果表明,所提方法的预测分布更集中,离散度低于其他方法。尽管平均训练时间(ATT)显示MTTNet 在时间成本上相较于其他基于 Transformer 的方法并无优势,但在工业应用中,RUL 预测精度优先于训练效率。因此,本方法能够高效捕获时间依赖,实现高精度 RUL 预测。    

    基于 NRMSE 对预测结果进行了配对检验,如图 17 所示。星号 (*) 表示两方法结果存在显著性差异(*p<0.05**p<0.01***p<0.001****p<0.001)。结果表明,本方法在约5% 统计显著性水平下优于其他八种方法。在短期 RUL 预测任务中,各方法预测性能差异较小;随着预测期延长,NRMSE 显著高于本方法的其他方法数量逐渐增加。

     
    图 17 不同方法多期 RUL 预测结果的统计分析  

    6 结论

    本文提出了一种用于滚动轴承剩余寿命预测的渐进式混合超图注意力网络(PHHGATN)。该方法通过双分支混合超图结构,同时建模高阶与低阶相关性,并实现自适应特征融合。在网络设计上,采用渐进式结构构建与通道信息融合策略,以突出关键节点与超边的重要性。基于两个全寿命数据集的系统性实验表明,该方法能够有效提升预测精度,并具备一定的可解释性;在变工况数据集上的案例研究进一步验证了其良好的泛化能力与计算效率。此外,本文还探讨了基于双通道信息的健康预测特性,结果显示该融合策略不仅有效增强了预测结果的因果性,也展现了双通道信息在设备健康预测中的应用潜力


    编辑:Jin
    校核:李正平、陈凯歌、赵栓栓、曹希铭、赵学功、白亮、任超、海洋、Tina、陈宇航、陈莹洁、肖鑫鑫赵诚
    该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除


    来源:故障诊断与python学习
    ACTSystem振动复合材料非线性船舶海洋理论电机材料传动机器人多尺度数字孪生控制装配
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2025-09-07
    最近编辑:11月前
    故障诊断与python学习
    硕士 签名征集中
    获赞 87粉丝 152文章 328课程 0
    点赞
    收藏
    作者推荐

    SCI一区论文学习 | 基于年龄和状态依赖的跳变扩散过程的退化系统的剩余寿命预测

    本期给大家推荐一篇SCI论文基于年龄和状态依赖的跳变扩散过程的退化系统的剩余使用寿命预测。剩余寿命(Remaining Useful Life, RUL)的预测在预测和健康管理(Prognostics and Health Management, PHM)领域具有重要意义,可以提供有关系统的寿命信息。有效估计RUL的基础是为系统构建一个适用的退化模型。然而,现有的大多数退化模型只考虑了年龄依赖性问题,而忽略了状态依赖性。此外,系统运行环境的变化将导致退化状态跳变,这将影响退化路径和RUL估计的精度。文章提出了一种用于RUL估计的广义年龄状态依赖跳变扩散(Age-State Dependent Jump-Diffusion, ASDJD)模型。基于首达时间(First Hitting Time, FHT)推导了RUL分布的近似解析表达式。期望条件最大化(Expectation Conditional Maximization, ECM)和极大似然估计(Maximum Likelihood Estimate, MLE)用于估计模型的未知参数。通过模拟数据集和西安交通大学轴承数据集验证了所提出的模型。论文链接:通过点击本文左下角的阅读原文进行在线阅读及下载。论文基本信息论文题目: Remaining useful life prognostic for degrading systems with age- and state-dependent jump-diffusion processes论文期刊: ISA Transactions论文时间:2025年论文链接:https://doi.org/10.1016/j.isatra.2024.12.019作者:Bincheng Wen (a), Mingqing Xiao (a), Xilang Tang (b), Yawei Ge (c), Xin Zhao (a), Haizhen Zhu (a)机构:a: ATS Lab, Air Force Engineering University, 710038 Xi’an, China;b: Equipment Management and UAV Engineering College, Air Force Engineering University, 710038 Xi’an, China;c:Strategic Evaluation and Consultation Center, Academy of Military Sciences, 100000 Beijing, China.通讯作者邮箱:1255292204@qq.com 目录1 摘要 2 引言 3 ASDJD模型 4 基于跳变扩散模型的RUL估计 4.1 理论推导 4.2 ASDJD模型说明 4.3 在役单元的RUL预测 5 参数估计 6 案例研究 5.1 数值模拟 5.2 实验数据 7 讨论和结论 摘要剩余使用寿命(Remaining Useful Life, RUL)的预测在预测和健康管理(Prognostics and Health Management, PHM)领域具有重要意义,可以提供有关系统的寿命信息。有效估计RUL的基础是为系统构建一个适用的退化模型。然而,现有的大多数退化模型只考虑了年龄依赖性问题,而忽略了状态依赖性。此外,系统运行环境的变化将导致退化状态跳变,这将影响退化路径和RUL估计的精度。然而,现有的研究只考虑了一些影响因素,忽略了同时考虑年龄依赖、状态依赖和跳变。为了同时考虑年龄状态依赖(Age-State Dependent, ASD)和跳变,本文提出了一种用于RUL估计的广义年龄状态依赖跳变扩散(Age-State Dependent Jump-Diffusion, ASDJD)模型。基于首达时间(First Hitting Time, FHT)推导了RUL分布的近似解析表达式。期望条件最大化(Expectation Conditional Maximization, ECM)和极大似然估计(Maximum Likelihood Estimate, MLE)用于估计模型的未知参数。通过模拟数据集和西安交通大学轴承数据集验证了所提出的模型,结果表明在RUL估计过程中应考虑状态依赖性和跳变。关键词:剩余寿命,跳变扩散,年龄-状态相关,预测1 引言对于关键组件和系统,随着工艺的恶化,预测方法力求获得准确的剩余使用寿命估计。基于可靠的预测结果,管理者可以确定设备维修的最佳时间,并制定相关的维护计划,旨在提高运行可靠性,降低风险和成本。目前,已采用多种方法来估算剩余使用寿命。这些方法大致可以分为基于物理模型的方法、机器学习(Machine Learning, ML)技术以及统计数据分析方法[1]。基于物理模型的方法通过分析系统的失效原理并开发相应的数学模型来表征系统的退化。这种方法常用于金属材料中裂纹扩展和疲劳退化的研究。例如,Eshwar等人[2]提出了使用扩展有限元方法对疲劳裂纹扩展进行建模,并利用卡尔曼滤波框架进行在线寿命预测。Wang等人[3]选择了帕里斯定律作为机身面板裂纹扩展模型,并使用扩展卡尔曼滤波估计模型参数。Deng等人[4]基于动态工作条件下的Archard磨损模型,提出了一种用于预测滚珠丝杠寿命的集成数据驱动粒子滤波模型。然而,各种降解系统具有独特的降解机制,因此开发一个通用的降解模型具有挑战性。此外,由于系统退化原理的复杂性,很难获取准确的物理退化模型。与基于物理模型的方法不同,机器学习不需要考虑复杂的故障机制。目前,人工神经网络(Artificial Neural Network, ANN)[5]、相关向量机(Relevance Vector Machine, RVM)[6]、递归神经网络(Recurrent Neural Network, RNN)[7]、门控递归单元(Gated Recurrent Unit, GRU)[8]和长短期记忆(Long Short-Term Memory, LSTM)[9]等机器学习方法在PHM领域得到了广泛的应用。此外,为了应对预测不同操作环境的挑战,迁移学习在估计RUL领域的应用越来越普遍。赵等人[10]开发了一种基于Bi-LSTM和深度迁移学习的预测框架。Chen等人[11]提出了一种用于跨领域健康因素构建的迁移学习框架。ML的优点之一是它能够自动执行特征提取和融合。另一方面,多层网络增强了处理非线性映射关系的能力。然而,这种类型的模型,也称为黑盒模型,无法解释机器学习的结果。此外,数据的数量和质量对模型的训练有重大影响。在实际情况下,很难获得足够数量的退化数据。与上述方法相比,统计数据驱动方法不需要完全了解降解机制或大量历史降解数据,并且由于RUL估计的有效性而被广泛采用[12]。Si等人[13]系统地回顾了常见的统计数据驱动方法,包括基于回归的模型[14,15]、基于马尔可夫的模型[16,17]、维纳过程(Wiener Process, WP)[18,19]、伽玛过程[20,21]。在这些方法中,WP因其能够模拟非单调退化而得到了广泛的应用。Zhang等人[22]对基于WP的退化模型及其在PHM中的应用进行了深入分析。Cai等人[23]使用线性WP研究了海底采油树系统的退化机制。Si等人[24]介绍了一种基于时空转换的非线性随机退化模型。Yu等人[12]提出了一种考虑三种不确定性来源的非线性WP模型。Gao等人[25]认为,外部环境会影响系统失效阈值,并提出了可变失效阈值的维纳退化模型。Zhang等人[26]提出了一种混合非线性函数WP,用于拟合轴承退化数据中的温度变化曲线,采用RVM进行参数学习和更新。扩展了当前的非线性情况。Chen等人[27]将非线性WP与GRU相结合,并在预测过程中使用滤波算法来更新模型参数。尽管维纳退化模型在预测RUL方面表现出色,但它仍有两个问题需要进一步研究:(1)ASD退化过程,(2)随机跳变问题。许多工业部件的退化不仅与年龄有关,还与状态有关。例如,滚珠轴承[28]、疲劳裂纹[29]、金属合金[30]和气缸套[31]都表现出这种特性。现有文献[32]提出了一个新的表达式,该表达式假设增量与当前退化水平存在隐式依赖关系。通过著名的巴黎-埃尔多安定律来描述疲劳应力下应力强度因子范围与亚临界裂纹扩展之间的关系。同样,文献[29]提出了一种取决于裂纹状态和监测时间间隔的退化模型。因此,疲劳裂纹扩展是一个更典型的ASD退化过程。这些ASD退化现象在实践中经常遇到。然而,关于ASD的研究是有限的。Giorgio等人[30,31]提出了一个包含ASD的马尔可夫模型,该模型已应用于气缸套和金属合金的磨损过程。Orchard等人[33]考虑了退化过程中的状态依赖性,并构建了一个基于粒子滤波的预测框架来预测电池的RUL。同样,An等人[34]使用粒子滤波构建了一个状态相关的退化模型。Zhang等人[28]提出了一个包含ASD的扩散模型,用于估计轴承的RUL。Li等人[35,36]介绍了一种使用ASD扩散模型模拟退化过程的方法。该方法专门用于分析疲劳裂纹扩展和涡扇发动机。Pang等人[37]引入了ASD的退化模型,该模型考虑了单个单元之间的固有变异性。然而,在上述研究中,只有[28,35,37]与ASD扩散过程模型有关,没有发现更多相关研究。此外,马尔可夫模型[30,31]和扩散模型[28,35]都没有解释的跳变对RUL估计的影响。在工业应用中,随着系统退化的加速,系统的稳定性逐渐恶化,系统退化状态发生跳变的概率增加[38]。这种现象通常被定义为退化-冲击依赖性。例如,文献[38]区分了轴承磨损和MEMS系统退化的不同阶段。随着系统随运行时间的累积退化,其抗冲击性降低,冲击对系统的影响将在不同阶段发生变化。此外,系统运行环境的突然变化可能会导致其退化状态的跳变。Wen等人[39]、Kong等人[40]和Zhang等人[41]通过根据跳变的位置将系统划分为不同的阶段,考虑了具有跳变的多相退化模型。然而,由于跳变次数和位置的随机性,上述模型不能准确反映实际的退化过程,导致RUL估计结果与实际情况之间存在显著差异,特别是在开始阶段。Pang等人[42]提出了一种具有跳变的单相WP,其中该模型中的跳变被描述为复合非均匀泊松过程(Compound non-Homogeneous Poisson Process, NHPP)。Zhang等人[43]提出了单相跳变-扩散模型,得到了RUL分布的解析解。Sun等人[44]得出结论,跳变对退化增量、退化速率和退化过程有影响,他们使用非线性WP构建了跳变-退化模型。时变copula函数构建了不同阶段失效过程的依赖关系。Wang等人[38]将系统分为四种状态:正常、缓慢退化、快速退化和失效退化;系统对跳变的抵抗力因州而异。然而,多相和单相模型在考虑冲击时都忽略了状态依赖性对退化过程的影响。基于上述扩散过程分析退化模型表明,虽然已经对年龄-状态依赖的扩散模型和年龄依赖的冲击-扩散模型进行了研究,但它只考虑了影响系统退化过程的少数因素。年龄-状态依赖的扩散模型和年龄依赖的冲击扩散模型都没有同时考虑年龄依赖、状态依赖和跳变。ASD扩散模型中跳变的考虑仍然是一个很大程度上未被探索的领域。为了准确表示系统的退化状态,将随机跳变整合到ASD模型中至关重要。本文提出了一个综合考虑跳变的ASD扩散模型。使用ECM和MLE进行参数估计。本文的重要贡献总结如下:(1)同时考虑退化中的年龄相关、状态相关和跳变问题,并构建了ASDJD模型(2) 通过It&#39;o公式和Larmpertier变换,得到了ASDJD模型寿命分布的解析解。此外,基于两步ECM和MLE对模型参数进行了估计和更新。(3)利用数值模拟和西安交通大学的轴承数据对模型进行了验证。实验结果证明了ASDJD模型的有效性。本文的其余部分结构如下:第2节介绍了。在第三节中,基于It&#39;o公式和Larmperti变换,从理论上推导了RUL分布,并给出了一个例子。第4节介绍了一个使用ECM和MLE估计参数的框架。第5节包含一个数值模拟和一个真实案例,以验证所提出的模型。第6节总结了整篇文章。2 ASDJD模型 设 表示系统的退化过程,即具有跳跃的ASD扩散过程,变量 表示系统在时间 的状态值。系统的退化状态可以用以下形式表示: 其中 是漂移项, 是布朗随机积分,而 表示复合泊松过程(Composite Poisson Process, CPP)。 表示标准布朗运动, 是漂移系数函数,而 是扩散系数函数。 和 是系数函数的参数。 是NHPP,即强度为 的齐次泊松过程(Homogeneous Poisson Process, HPP)。 表示跳变幅度,是一个独立同分布的随机变量,有 。通常假设 。备注1.(1) 如果 , 则表示为参考文献[42]中的年龄依赖跳变扩散模型;(2) 如果 , 则表示为参考文献[24]中的年龄依赖扩散模型;(3) 如果 , 则表示为参考文献[45]中的状态依赖扩散模型。因此,所提出的模型扩展了现有方法。由于 是由布朗运动驱动的WP,其寿命被定义为FHT。 其中寿命的概率密度函数(Probability Density Function, PDF)表示为 RUL的定义如下: 其中 为RUL,其PDF为 。3 基于跳变扩散模型的RUL估计3.1 理论推导为了获得ASDJD模型的RUL分布,有必要将退化模型(1)转换为包含常数扩散系数的退化模型。因此,采用了Larmperti变换[47]和Ito公式[28,46],其中前者用于将随机过程 转换为随机过程 。本文定义Larmperti变换如下: 为了进行推导的下一步,首先给出Ito公式中的引理1。引理1.([48]):考虑随机过程 ,它代表一个具有跳跃的ASD扩散过程: 然后,对于任意 函数 ,过程 可以表示如下: 在本文中, 被定义为(4)。我们结合(1)、(4)和引理1。经过许多数学操作,(1)被转换为一个常数扩散系数。其中CPP 和漂移系数函数 由引理1获得。 其中, 描述了第 次冲击时刻退化状态的左极限,而 是拉姆佩蒂变换后退化状态的值。显然, 是连续的。根据定积分的第一中值定理, 可以简化为: 备注2. 通过利用引理1和Larmperti变换,广义ASDJD模型(1)可以等价地转换为一个具有常数扩散系数的随机过程(7)。在不牺牲一般性的前提下,接下来的部分将集中讨论具有常数扩散系数的ASDJD模型,其定义如下: ASD模型(11)中跳变的存在性给RUL的解析解带来了挑战,首先提供引理2来解决这个问题。引理2.对于泊松过程,均值定义为 值得注意的是, 是泊松过程的参数。如果函数 等于零,则退化过程的跳变不存在。如果函数 为常数,泊松过程可以归类为HPP。如果 表示强度函数,则泊松过程可以归类为NHPP。为简化表达式,下文中将, 记作 。随机变量 服从高斯分布且相互独立同分布。泊松过程 的均值和方差均为 。因此,基于高斯过程([42,43])的条件概率过程(CPP),可近似替换为高斯过程。 表达式可表示为 其中 , 。可将等式(11)所代表的模型的退化重新表述如下: 的寿命概率密度函数可以通过一定的数学表达式进行估算,如[24,35]的研究结果所示。 其中 有 且 和 可以表示为: 在实际应用中,给定时间点 处的退化条件可表示为 。RUL的失效分布函数(FDF)被定义为退化过程首次超过预定阈值的时刻。时间点 处的退化过程 可描述为: 类似地,RUL的PDF可以用数学方式表示为: 其中, 有 3.2 ASDJD模型说明在本文中,ASDJD模型被表示为M1,其中 , 。根据3.1的结果,我们可以推导出相应的变换结果。当系统起始时刻为0, 时,可以得到 , , , , 的解析表达式。基于(15),可以计算出寿命概率密度函数。 其中, 和 服从高斯分布。根据全概率定律,可以得到寿命的概率分布 。为了避免计算过程中的积分,这里首先引入以下公式。 当 且 根据(22),寿命PDF可以基于(21)得出。 其中, 需要注意的是,(23)式中仍包含期望值。由于其计算过程较为复杂,已无法给出解析表达式。不过,(23)式仅包含单变量积分,因此可以使用Matlab进行近似计算。 当 ,表示系统的退化状态为 。 时。可得到对应的函数 。RUL的概率密度函数也可用类似方法获得。 根据(22),可以得到以下表达式。 其中, 图1 所提出模型的流程图 4 参数估计要实现寿命和剩余使用时间(RUL)的估计,确定所提模型的参数至关重要。我们将模型(14)的参数向量记为 。 是漂移系数的参数,其中 。 是泊松过程的参数。假设在时间区间 内存在 次跳跃事件,这些事件无法直接观测,可视为系统的隐变量。在此特定情况下,最大似然估计法不适用。由于隐变量 具有时变特性,传统的混合高斯分布期望最大化求解方法失效。因此采用隐马尔可夫模型(ECM)进行参数估计。根据泊松过程的定义,若满足条件 ,则可推导出以下公式: 备注3. 根据公式(28),当 足够小时,若满足条件 ,则概率 。因此,当满足条件时,我们可以定义概率为 ,其中参数值 由BIC或AIC [42,43]确定。对于退化系统,设 表示退化数据, 表示退化状态增量。 服从混合高斯分布,其均值为 ,方差为 。为简化计算,令 。时间增量表示为 ,其中 。通常传感器的工作频率较高且固定,因此实际采样间隔足够小以保持恒定。在后续章节中, 。 表示隐变量 的权重。第一步: ECM参数估计E步: 我们需要计算完整似然的期望值,其表示为Q函数如下所示。 其中 表示对应的第 个参数向量。 是潜在变量。CM步骤: CM-1:固定 。令 ,并对其求最大值,即可得到如BoxI所示的 由于无法推导出 的解析解,因此我们采用粒子群优化技术来求解方程(32)。 CM-2: 固定 ,我们最大化 , 可以表示为: CM-3:固定 与CM-1类似, 的解析解无法推导,因此采用粒子群优化算法求解(34)。 两步法: 通过ECM方法,我们已获得参数 。现在只需确定参数 ,利用最大似然估计法来估计该参数。 其中 表示为: 通过以上讨论可知,利用ECM和MLE可以实现对所有未知参数的估计,为了便于理解,图2中总结了参数估计的步骤。 5 案例研究本节通过数值模拟和fXJTU-SY的实际案例两个实例,验证了所提出的模型的有效性。5.1 数值模拟本小节验证了参数估计和RUL估计模型的准确性。首先,基于离散化退化模型(37),我们使用欧拉近似生成D条退化路径[24,42]。 模拟模型的预设参数见表1。下面给出了泊松过程的强度函数 : 表1 模拟的参数设置 需要说明的是,本文采用的模拟冲击模型仅为示例模型,并不能完全反映工业退化的全貌。实际应用中,必须根据具体退化特征和运行环境构建对应的冲击模型强度函数。我们基于初始设定参数推导出退化路径:将历史数据集(即指定起始点前的数据)称为训练集,起始点后的数据则称为测试集。如图3(a)所示为退化路径示意图,图3(b)展示了数据的增量退化过程。根据图3呈现的数据,跳变频率呈现出逐步上升的趋势。 为评估所提模型的效能,我们采用蒙特卡洛(MC)模拟生成了15,000条降解轨迹,并从中获取实际概率密度函数分布。本文提出的模型标记为M1。通过两种对比方法验证模型有效性:(1)无跳跃的ASD扩散模型称为M2 [28];(2)与年龄相关的跳跃-扩散模型则标记为M3 [43]。根据第5节提出的ECM和MLE方法,我们成功获取了表2所列的ASDJD模型参数。初始参数通过训练数据学习获得,并在输入新测试数据时进行更新。从表2可以看出,退化参数能够通过ECM和MLE方法准确学习。表2 多个模型的参数估计结果图4展示了M1、M2和M3三种模型在不同初始时间点的RUL概率密度函数(PDF)结果。通过蒙特卡洛模拟获得的精确系统PDF表明,M1模型的RUL PDF与实际蒙特卡洛模拟结果高度吻合。实证研究表明,采用M1模型参数估算的解析PDF与蒙特卡洛方法得出的PDF具有极强的拟合度。值得注意的是,在 秒时,M2模型的估计误差明显大于M1和M3。这种误差值的差异揭示了状态跳跃对RUL估算的影响。随着设备老化时间增加,三种模型的估计误差均呈现下降趋势,但M1模型仍保持优于M2和M3的性能优势。需要特别说明的是,M2和M3的估算结果低于实际RUL值,这可能源于未考虑状态跳跃或状态依赖性的影响,从而导致系统维护延迟并增加制造成本。 5.2 实验数据示例验证中使用的模型的漂移和扩散系数函数如第3.2节所示,冲击的密度函数如(38)所示。由于模型函数已经过验证,所使用的退化数据应包含相应的退化特征,即指数退化,且冲击概率逐渐增加。5.2.1数据描述和评价指标本小节[49]采用的是XJTU-SY轴承退化数据集。该数据集包含15个滚动轴承样本,采样频率设定为25.6千赫兹(kHz),采样周期为一分钟。在评估剩余使用寿命(RUL)预测性能时,我们采用了三个关键指标:均方根误差(RMSE)、绝对误差(AE)以及累积相对精度(CRA)。 变量 代表RUL的真值,而 表示预测的RUL。 在[49]系统中,水平和垂直振动信号的最大振幅(MA)被用作系统退化特征,即每个采样点处采样振动信号的最大值。设定的退化特征失效阈值为20 g。采用轴承1-3和2-2来评估所提模型的有效性。完整信息详见表3。表3 XJTU-SY轴承数据集5.2.2.剩余寿命预测实验结果轴承退化增量如图5所示。显然,轴承1-3的增量变化趋势与仿真数据相似。预测的轴承1-3和轴承2-2的起始点分别为 秒和 秒。为更好地展示模型的应用效果,我们将在以下步骤中进行具体实施。模型流程如图1所示。首先,根据采集的振动信号,提取出模态分析(MA)作为退化数据 。由此可获得模态分析的增量参数 。 第二,通过所提出的参数估计方法,可获取1-3号轴承和2-2号轴承的模型参数。系统将利用历史数据及预测起始点前的数据来估算初始参数。随着新数据点的获取,参数将通过ECM和MLE方法进行更新。 第三,基于这些参数估计结果,我们能够预测设备剩余使用寿命(RUL)并计算其对应概率密度。 图6展示了轴承1-3和轴承2-2的RUL概率密度函数(PDF)。从图6(a)和(b)可以看出,随着退化数据的累积,RUL概率分布的标准差逐渐减小,数据点也变得更加集中。在轴承1-3和2-2中,M3的概率密度分布比M1更集中,这表明考虑状态依赖性可能会影响模型的不确定性。不过M1的预测精度仍优于M3,这一优势在图7和图8中均有体现。特别是对于轴承2-2,M1的表现明显优于M2和M3。这些结果表明,具有跳跃特性的年龄-状态依赖扩散模型能更好地拟合退化数据。通过对比M1、M2和M3可以看出,M1的泛化能力比M2和M3更强。对于轴承1-3,M1和M2的预测结果呈现相似趋势;但针对轴承2-2,M1的表现则显著优于M2和M3。 图7展示了RUL预测结果。黑色实线表示实际RUL值,蓝色阴影区域标示了误差的20%置信区间。在轴承1-3和轴承2-2工况下,M1和M3的预测表现优于M2,这表明考虑跳跃效应的模型能更好地拟合退化数据。对于轴承1-3工况,M1和M2的RUL估计值呈现相同趋势,但M1的预测精度更高。在轴承2-2工况中,M1的估计结果显著优于M3,这说明考虑状态依赖性能提升模型预测精度。由此可见,在轴承1-3和轴承2-2工况下,M1在多数情况下均优于M2和M3。为进一步验证结果的可靠性,我们采用均方根误差(RMSE)和平均绝对误差(AE)对M1、M2和M3模型进行评估。如图8所示,M1模型的平均绝对误差明显优于M2和M3,这表明所提出的M1模型在资源利用率(RUL)估计方面表现更优。值得注意的是,M2模型在两个数据集中的估计误差均最为显著,反映出其估计结果存在较大的波动性。表4展示了各模型的RMSE值:M1、M2和M3在基准方向1-3上的RMSE分别为134.02、1226.22和206.40。通过对比可以发现,M1模型的RMSE值显著低于M2和M3,这说明M1模型不仅稳定性更强,且方差更小。图8 RUL预测的AE表4 不同模型的均方根误差图9为三个模型的CRA [12,35],CRA是量化估计模型准确性的指标,CRA值越大,估计结果的精度越高。提出方法的结果表明,与另外两个模型相比,1-3号轴承和2-2号轴承的CRA值更高,这说明所提出的模型可以得到更精确的结果。图9 三种模型CRA的比较上述对比研究表明,所提出的方法能够显著提升剩余寿命(RUL)估计的准确性。实验结果表明,在建模退化过程时,必须充分考虑跳跃效应和状态依赖性的影响。但需注意的是,轴承2-2的预测精度低于轴承1-3。通过观察图5可以发现,轴承2-2的退化增量与轴承1-3存在偏差,这说明本文提出的冲击密度函数无法很好地拟合轴承2-2的特性。因此,轴承2-2的预测结果偏差明显大于轴承1-3。6 讨论和结论本文提出了一种具有随机跳跃特性的ASD扩散模型构建方法,旨在实现退化系统的剩余寿命(RUL)预测。现有模型(如文献[24,37,43]所述)可视为该方法的典型实例。通过引理1和Larmperti变换,我们将退化模型转化为具有恒定扩散系统的ASDJD模型,并采用NHPP分布描述时变随机跳跃特性。借助时空变换技术,我们获得了RUL概率密度函数(PDF)的解析表达式。未知参数估计采用两步法,结合ECM算法与最大似然估计(MLE)技术完成。为验证模型有效性,我们使用西安交通大学提供的仿真数据集和方位角数据进行对比分析。将本模型与未包含跳跃特性的ASD模型、年龄依赖型跳跃模型进行比较后发现,本模型能更精准地估算RUL。现有研究仍存在若干亟待深入探讨的局限性:(1)实际应用中,由于测量对象间的个体差异和系统性误差,数据变异性普遍存在。因此,如何将这种变异性纳入模型并推导解析解仍需进一步研究。(2)工业应用中的跳跃过程并非固定高斯随机变量,可能遵循时变高斯分布。这些特性将在后续研究中重点探讨。(3)系统退化过程并不总是符合指数漂移函数,冲击过程也不总遵循非负极点过程(NHPP)。如何提升模型的泛化能力,或根据数据退化特征自动选择适用的退化模型,仍是当前亟待攻克的全新难题。编辑:Tian校核:李正平、陈凯歌、赵栓栓、曹希铭、赵学功、白亮、任超、陈宇航、海洋该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

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