首页/文章/ 详情

MSSP论文学习|一种基于物理知识嵌入与通道注意力的工业机器人寿命预测方法

9月前浏览836
     本期聚焦工业机器人预测性维护新范式:在柔性制造系统中,工业机器人长期运行于复杂多变工况,其核心部件的性能退化直接威胁生产精度与系统可靠性,而传统数据驱动方法严重依赖大量标注数据,在实际部署中成本高昂。为此,本文突破性地提出一种物理知识嵌入的通道注意力神经网络(PICAIU),通过将机器人动力学知识与伺服控制器信息深度融入数据驱动模型,有效引导网络学习具备物理意义的退化特征。该方法凭借知识-数据融合驱动,将训练所需轮次大幅下降;同时,模型通过通道注意力机制融合多物理特征,实现了对电机电流的精确估计,使其剩余使用寿命预测误差较最佳基线模型降低26%;在跨机器人验证中,面对从“数据丰富”到“数据贫乏”的苛刻条件,预测误差仍稳定于6.4%的低位,展现出在样本稀缺工业场景下的强大泛化能力。该研究为工业机器人状态评估与预测性维护提供了高精度、高效率与强泛化的全新解决方案。  

    论文基本信息

    论文题目:

    A degradation assess ment method for industrial robot in flexible manufacturing systems based on multiple physics-informed neural network  

    论文期刊:Mechanical Systems and Signal Processing

    论文日期:2025年

    论文链接:

    https://doi.org/10.1016/j.ymssp.2025.112793

    作者:Shijie Wang (a), Jianfent Tao (a*), Qinchent Jiang (a), Wei Chen (b), Chengliang Liu (a)

    机构:

    a. State Key Laboratory of Mechanical Systems and Vibration, Shang Hai Jiao Tong University, Shanghai 201100, China;

    b. STEP Robot Co., Ltd., Shanghai 201800, China

    第一作者简介:陶建峰老师现任上海交通大学机械与动力工程学院机电控制与物流装备研究所教授。陶老师长期致力于机电系统智能控制与故障诊断领域的研究,主要研究方向包括智能液压系统控制及基础元件设计、复杂机电系统及基础元件故障诊断等,推动智能运维技术在高端装备领域的创新应用。在科研项目方面,他主持及参与了多项国家级重要课题,包括国家重点研发计划课题“射流伺服阀性能退化机理及健康管理”,以及国家自然科学基金面上项目“类矩形盾构环臂式电液伺服管片拼装机振动机理及其预测控制方法”。相关研究成果已在《Mechanical Systems and Signal Processing》、《Journal of Intelligent Manufacturing》等国际高水平期刊上发表多篇学术论文。来源: https://me.sjtu.edu.cn/teacher_directory1/taojianfeng.html

    目录

    1 摘要

    2 引言

    3 基本概念

    3.1 物理信息神经网络

    3.2 机器人动力学牛顿-欧拉方程

    3.3 电机控制系统响应函数

    3.4 通道注意力机制

    3.5 改进的KL散度

    4 方法论

    4.1 物理信息神经网络理论

    4.2 通道注意力-Inception-Unet框架

    4.3 基于改进的KL散度的退化指标

    5 实验与分析

    5.1 实验设置

    5.2 时间序列建模实验

    5.3 减速器加速老化数据集的剩余寿命预测

    5.4 PICAIU鲁棒性测试

    5.5 多工况退化实验

    6 结论

    1 摘要

    随着工业机器人应用范围的扩大,性能退化监测对制造系统具有重要意义。本文提出了一种名为物理信息通道注意力Inception-Unet(Physical-Informed Channel-Attention-Inception-Unet, PICAIU)的框架,用于预测工业机器人的剩余使用寿命(RUL),该框架可应用于评估减速器性能退化。为提高训练效率和预测精度,该方法将机器人动力学知识和伺服控制器信息嵌入数据驱动神经网络,并在Inception-Unet中引入通道注意力机制以融合运动特征与功率特征,从而获得精确的电机电流估计值。基于机器人谐波齿轮传动的老化实验数据,建立了可适用于其他机器人的相关寿命曲线。通过系列实验验证了所提方法的鲁棒性与泛化能力:消融实验表明,物理信息约束的引入提升了训练效率与收敛精度,有效缓解了通道注意力机制对样本量的需求;RUL预测误差较最佳基线模型降低26%;在同批次条件下,PICAIU仅需基线模型42%的训练轮次即可使拟合误差低于10%,显著节约训练时间与样本成本。通过双机器人交叉验证,在从丰富数据集泛化至贫乏数据集时,预测误差仅为6.4%,这在样本量受限的实际工业场景中具有应用潜力。

    关键词:性能退化评估;柔性工况;知识驱动;工业机器人

    2 引言

    当前智能制造发展背景下,远程健康监测与故障诊断技术对保障复杂实时系统的稳定运行具有重要意义。工程师通过传感器采集现场数据,依据专业知识评估设备状态、分析损伤机理并制定维护策略。作为生产线核心装备,工业机器人的健康状况直接关系到生产安全与效率。健康指标通过单一量化参数综合反映关键部件健康水平,是确定维护时机的关键依据。例如,轴承变形会加剧偏心旋转导致磨损,减速器损坏则易引发齿轮异响。因此,利用不同损伤状态下部件的动态特性差异,实现对机器人内部健康状态的实时监测与准确评估,对保障其安全运行、延长服役寿命至关重要。

    在退化评估方面,轴承等重要部件已形成较成熟的技术路线。机器人领域虽有过基于扭矩偏差的评估方法探索,但因缺乏兼顾精度与效率的动力学仿真方案,相关研究进展有限。早期健康评估依赖人工特征(如有效值、峰值),在发动机等设备中效果显著。该类设备数据完整、要求严苛,而工业机器人作为柔性制造装备,其工况多变、性能波动常被接受,为健康指标构建带来挑战。目前主流方法通过基准数据对比或电流信号拟合工作循环,但仍受限于工况适应性

    工业机器人本质上是一种柔性复杂的机电系统。如图1所示,对于轴承等在稳定工况下运行的部件,其退化过程可借助历史序列数据进行推断。然而,在柔性制造环境中,机器人常需应对前所未有的新型工况,导致难以获得可资借鉴的历史数据。传统故障诊断方法大多基于开源数据集训练时间序列分类器以区分故障模式,这类研究在工业机器人场景中适用性有限。随着研究的深入,学者们逐渐认识到退化过程与控制理论中系统观测器的内在关联。

    图1. 不同工况的退化特性

    近期数据驱动研究侧重从退化数据中挖掘分类或回归模型,例如通过迁移学习与卷积网络识别故障,但过度依赖数据导致模型可解释性差。当前共识是,应结合领域知识推动人工智能发展。现有剩余寿命预测方法多针对固定工况设备,而工业机器人需摆脱对运动轨迹的依赖,建立适用于任意程序的统一健康评估框架。对于工业机器人,企业期望获得不依赖运动轨迹的健康评估方法,为任意运动程序建立统一的健康指标框架。基于上述研究现状,本文提出基于知识的通道注意力机制与多尺度卷积神经网络,以学习机器人关节中多维物理信号间的潜在动态关系。通过真实数据训练网络,节点计算模型输出与实际数据间的差异以判断系统是否存在退化。本文主要贡献如下:

    • 提出面向柔性制造系统工业机器人的健康评估框架,通过估计电流与实测电流的改进KL散度生成退化指标,实现跨轨迹的通用化健康评估;

    • 将牛顿-欧拉动力学与伺服电机控制关系嵌入损失项,通过机器人知识增强电机电流仿真神经网络的泛化能力与可靠性;

    • 设计基于通道注意力的多序列关系建模方法,描述多轴机械臂不同关节间的相互作用,提升电流估计精度;

    • 通过减速器加速寿命试验,系统验证所提退化预测方法相较传统方法的泛化优势,展示其对负载与轨迹变化的适应能力。

    3 基本概念

    3.1 物理信息神经网络

    物理信息神经网络(PINN)通过在损失函数中引入物理约束,该技术能够利用较少数据构建出泛化能力更强的模型。物理约束通过引入人类先验知识来指导神经网络学习,这一机制不仅能提升训练效率、降低数据采集成本,对于工业机器人的退化预测研究更具重要价值。此外,知识嵌入还有助于引导神经网络聚焦于关键特征的学习,从而避免被无意义噪声干扰,最终提高模型训练精度。受此启发,本文尝试将工业机器人的多层次知识体系引入深度学习任务,其损失函数如下式所示:

             

    损失函数由两部分构成:基于机器人领域知识推导的物理方程计算得出 物理约束损失           ;通过样本真实值与网络输出值之间的差异计算得出 数据拟合损失          。本研究的核心思路正是通过构建物理约束损失,有效提升神经网络的数据利用效率。

    3.2 机器人动力学牛顿-欧拉方程          

    机器人动力学方程如下式所示。等式左侧表征机器人的运动学与构型参数,右侧则对应关节扭矩等动力系统信号。其中,矩阵 D 为惯性矩阵,包含转动惯量与惯性积等参数:

               

    其中,C 为阻尼矩阵,涵盖粘滞力系数、阻尼比等系统参数;M 为质量矩阵,反映机器人位姿的初始配置及与质量分布相关的参数;τ 与 σ 则对应于摩擦力、电磁阻尼等二次效应。q 表示广义位移,本文中主要指电机的旋转位移。N 表示扭矩。此外,牛顿-欧拉方程包含一套从末端负载力到基座扭矩的迭代方程组。

               

    其中,             表示从关节              到关节              的旋转矩阵。             与              分别为关节              所承受的负载扭矩与驱动力矩。             与              分别表示关节              的质心位置及机器人所受合外力。在执行搬运任务时,合外力主要指机器人本体与负载对象的重力之和。鉴于伺服电机扭矩可通过扭矩系数与电流值计算获得,其扭矩满足以下关系式:

                 

    其中,系数                可直接从电机铭牌获取,在大多数运行工况下可视为稳定常数。机器人六个关节间的电流相互作用关系可表述为:

                 

    其中,                为运动微分方程对应的雅可比矩阵,               为相应的常数项。该公式描述了动力学约束关系,该约束是通过牛顿-欧拉方法推导得到的电流相互作用关系。

    3.3 电机控制系统响应函数          

    本文引入的另一项知识约束源于工业机器人采用的伺服电机。图 2 所示为伺服驱动下机器人关节系统的拉普拉斯模型:目标位置以给定电压            的形式通过电位器输入,经升压模块转换为            ,通过电机电枢电感与电阻转换为电流            ,再借助前述系数            在磁场中转化为给定电磁扭矩            。该扭矩与实际扭矩            反馈比较后,经过惯性环节与摩擦环节作用形成实际转速,最终通过减速器传动比            转换为机器人关节的实际角度            

             

    图2. 伺服电机系统模型

    电流信号直接反映电磁扭矩,故              可表征电流信号的幅值。该系统以电压为输入、角度为输出,其传递关系可表述为:

                 

    若将电流作为观测指标置于等号左侧,则该方程可转换为以下形式:

                 

    表 1 汇总了本文用于退化研究的机器人关节电机参数。由于电机退化程度相较于减速器或同步带而言相对较弱,其参数可作为实际分析中的现场校准参数。

    表1 电机系统方程参数表

             

    3.4 通道注意力机制            

    本文采用的通道注意力机制为挤压激励模块,该模块能够动态调整特征的重要性权重,突出更具价值的信息并抑制次要信息,从而提升模型对关键特征的关注度与整体性能。其实现过程包含三个步骤,具体如下:挤压操作:对于输入特征图             (其中              为通道数,             为特征图高度,             为宽度),首先进行全局平均池化处理,将每个通道的二维特征压缩为单个实数,得到             ,其计算公式为:

                 

    激励操作:随后对                执行两次全连接层操作。第一个全连接层将通道数降至                              为缩减率),并采用 ReLU 激活函数;第二个全连接层将通道数恢复至               。最终获得输出               ,其计算公式为:

                 

    其中                为 Sigmoid 函数,               为 ReLU 激活函数,               与                分别为两个全连接层的权重参数。

    缩放操作:最终将获得的通道权重                  与输入特征图                  逐通道相乘,得到输出特征图                 ,其计算公式为:

                   

    通过上述步骤,通道注意力机制实现了对不同通道的自适应加权,从而突出重要特征通道并抑制非重要通道。

    3.5 改进的KL散度              

    概率分布间的差异常采用KL散度进行度量。然而,由于其特定计算形式,KL散度的计算过程通常较为复杂。对于两个概率分布                  与                 ,其KL散度计算公式为:

                     

    对于连续概率分布                    与                   ,其计算公式为:

                     

    尽管KL散度具有重要应用价值,但其计算复杂度较高。为缓解此问题,本文提出改进KL散度以优化传统KL散度的计算效率。对于两个概率分布                    与                   ,MKL散度定义为:

                       

    对于连续概率分布                      与                     ,其计算公式可表示为:

                       

    该改进形式通过引入平方根运算替代原始对数分量,有效降低了计算复杂度。在保持对两分布差异捕捉能力的同时,显著提升了计算效率。

    方法论

    4.1 物理信息神经网络理论

    传统研究方法通常通过对比当前与历史时序数据,采用时域或频域指标来描述系统特定特征的退化过程。然而对于柔性场景中的工业机器人,轨迹与负载的改变同样会引起关节电流形态的变化,这一直是工业机器人剩余寿命预测面临的挑战。

    本文从机器人动力学仿真中获得启发,尝试以仿真数据与采集数据间的系统性偏差作为不同任务的统一退化评估框架。如图3所示,这一引入机器人仿真知识的框架称为物理信息神经网络(PINN),包含数据采集框架与深度学习转换框架。

     

    3. 基于物理信息神经网络(PINN)的退化指标评估框架

    该框架工作原理为:从PLC获取DH参数、关节角度、电源电压等机器人信息,将所有数据输入基于通道注意力机制的神经网络。本文以电源电压、运动角度和负载质量作为输入,以电机电流作为输出,训练这个引入知识约束的神经网络作为机器人动力学仿真模型。在物理约束引导和真实数据支撑下,该仿真器能够针对任何真实输入高精度稳定地输出对应电流。随后,利用估计电流与实际电流的改进KL散度作为参考,获得随时间变化的退化指标。通过机器人耐久性测试数据,可根据当前退化程度预测未来退化趋势,并制定相应的故障应对策略。  

    4.2 通道注意力-Inception-Unet框架      

    本文设计的跨工况评估框架利用采集的健康状态数据训练代理模型,将实际运行信号输入代理模型后获得理论健康电流,其与实际电流的偏差被视为退化参考指标。本文设计了一种基于通道注意力机制的神经网络,该网络包含通道注意力模块和改进的UNet结构,与当前主流的自编码器框架或生成对抗网络框架存在显著差异。

    该框架的主要特点之一是采用并行卷积结构实现图4中所示的多尺度卷积运算。该卷积模块被定义为引入模块,可改善传统卷积核不擅长提取多尺度特征的缺陷,从而更有效地学习时间序列的潜在信息。

       

    图4通道注意力模块

    多关节间的耦合效应会严重影响动力学求解精度,而如何从实时数据中识别当前关节耦合效应一直是机器人领域的难点。通道注意力机制能够增强多样本中的关键信息并抑制无效信息,因此本文将其自然应用于评估多关节多信号通道间的权重对应关系。针对本文研究的六关节机器人,输入信号包含六轴的关节角度、关节速度与关节电机功率。待处理的各序列长度为2048,故通道注意力模块的输入尺寸为2048。由于后续的初始网络将进行精细特征提取,本文未压缩通道数量,因此设定参数H=2048W=1C=18。涉及的部分神经网络计算函数如下所示

           

           

           

           

    基于SE注意力机制,实现了5所示的通道注意力模块。损失项d通过实际电流与估计电流的偏差计算得出。

           
    图5. 通道注意力- Inception-Unet 框架      

    基于双序列的损失函数计算公式如下:      

           

    神经网络输出的电流序列记为         ,实际采集的电流序列记为         ,基于牛顿-欧拉方程计算得到的电流记为         ,通过系统方程得到的电流记为         。三种损失项的计算形式如下:

           

           

           

            通过下式计算得出:

           

           

           

            通过下式计算得出:

           

    其中          表示拉普拉斯逆变换,         为实数,用于确保积分路径位于          的收敛域内,         为虚数单位。物理约束损失与样本约束损失对应关系如下:

             

             

    包含物理信息约束与数据标签约束的损失函数可表示为:

             

    4.3 基于改进的KL散度的退化指标        

    图 6 展示了利用KL散度计算机器人关节的剩余使用寿命。经过训练的网络等效于真实机器人的个性化数字孪生体,该孪生体记录了机器人处于完全健康状态时的潜在动力学关系。通过向孪生机器人输入实际运动与振动信号,可获得设备在完全健康状态下的理论扭矩输出。在此框架下,实际扭矩信号与理论扭矩信号构成两组不同的概率分布:当前状态下的实际电流记为          ,估计电流记为          

         

    图6. 基于改进KL散度的退化指标

    初始的实际电流为           ,估计电流为           

               

    改进KL散度能够有效度量先验分布            与后验分布            之间的差异,从而帮助我们理解观测数据对先验分布的影响程度。通过退化测试获得的模型将在每个实际应用场景中进行针对性调整,并根据现场任务差异持续优化。通过提出的PICAIU(物理信息通道注意力-Inception-Unet)架构,本文尝试引导神经网络从工业机器人关键状态序列中学习本质相关特征。

    5 实验与分析

    5.1 实验设置

    实验系统采用 JAKA Zu-12 和 JAKA Zu-16 两款机器人平台。两款机器人额定负载相差4kg且配备不同型号减速器,但末端法兰规格完全一致,可执行相同的加工任务。其中Zu-12老化速率较快,故以其数据作为主要实验样本,Zu-16 数据仅用于验证模型泛化性能。系统配备负载敏感型电主轴与激光扫描仪,可通过扫描工作台工件获取表面点云数据,并为焊缝抛光任务规划实时运动轨迹。实验系统采用JAKA Zu-12与JAKA Zu-16两台机器人平台开展验证。如图 7 所示,待加工件为高碳钢焊接构件,需对4条厚度约4cm的焊缝进行抛光处理。

           

    7. 机器人实验制造系统          

    由于工件更换耗时较长,实验选择0.1–0.15mm范围的随机进给量进行每次抛光加工。为加速部件磨损,实验前将减速器内部润滑脂全部清除,迫使齿轮以干摩擦形式啮合。该实验针对高碳钢焊缝抛光场景,机器人末端关节承受的负载力在600–1000N区间波动,经数日连续加工后可观察到减速器的显著损伤。

    图8展示了某测试工况下采集的实际数据,包含关节角度、电机驱动功率与电机电流。这些数据通过PLC网络端口上传至云端数据库。实验采用的协作机器人各电机功率不超过1.5kW,驱动电流小于2A,数据采集协议采用TCP/IP,采样频率为250Hz。

           

    图8. 实验采集的多源机器人信号

    采集信号包含负载运动与空载运动两类工况。负载运动时机器人通过线激光自动规划轨迹执行焊缝抛光任务,空载运动时则随机执行各类测试动作。作业期间,执行机构通过高压气体与控制阀确保工件额定压力,通过调节压力阀开度可将磨削过程中的反作用力控制在恒定值。工件安装有三个振动传感器,用于精确判定切削过程的起止时间点。图8所示曲线为根据振动幅度剔除换刀、退刀等无效动作后的有效数据。

    图9展示了更详细的硬件配置。实验采用的柔性制造系统通过可变位金属结构与差异化焊缝形态生成动态机器人作业轨迹。

           

    图9. 机器人柔性制造系统详图

    红色管路连接的高压气体为电主轴提供压力反馈与负载采集,抛光过程中的反作用力可直接通过蓝色压力表及转换系数获取。图9左下角点云图像呈现了加工过程中扫描获得的表面特征。由于焊接工艺的一致性难以保证,即便是同一工件,每次焊接后所需的磨削路径也会存在差异。

    5.2 时间序列建模实验        

    如表2所示,本实验根据动作类型、运行速度与负载条件划分为40组(4×3×16=192种工况)。每种工况均具有不同的轨迹特征,具体体现在运动形态与电流信号上。负载重量的变化通过安装在机器人末端的气动伺服磨削主轴实现。

    表2 工况类型划分

         

    其中,运动类型根据编程特性划分,本文采用的JAKA工业机器人基本运动功能包括直线、圆弧、插补和末端位姿四种类型,通过组合这些基本运动可实现任意复杂度的末端轨迹。在抛光图9所示工件时,四种运动类型均会涉及。为便于后续深入分析,根据运动功能调用对样本进行分类。图10展示了基于机器人学知识的电流预测精度。蓝色曲线代表实际电流,橙色曲线代表利用动力学递归关系与电机伺服系统传递函数估计的电流

         

    图10. 传统仿真方法电流估计结果

    可以看出:虽然体现机器人动力学知识的数据1估计电流与实际电流整体趋势一致,但无法准确估计高频电磁振荡;代表伺服系统知识的数据2估计电流虽不能准确反映宏观趋势,但保留了与电源和编码器相关的高频信息。由此可见,传统知识驱动的信号仿真方法只能在有限领域有效还原物理系统信息,但由于知识本身的前提假设限制,最终结果与真实信号间始终存在偏差。图11呈现了本文通道注意力模块学习得到的六路输入信号与六路输出电流间的通道权重矩阵。

         

    图11. 通道注意力模块生成的评分矩阵

    从热力图中可清晰观察到,高权重区域集中分布于主对角线附近及边界区域,且距离主对角线越远权重评分越高。当某个关节受到冲击时,距离会放大冲击产生的扭矩效应,因此沿轴向距离越远的关节承受的冲击影响越大。这表明通道注意力机制等效于学习了多关节机器人的串联特性。与预期相反的是,通道注意力获得的评分矩阵并非对称矩阵,这可能是由于机器人关节耦合效应会随传动方向产生不同形式的负载影响。

    图12展示了引入动力学模型与电机伺服特性物理约束后的电流估计样本。其电流估计性能优于图10所示结果,表明该模型能准确观测健康机器人的系统特性。通过输入特定任务场景参数并估计健康电流后,可基于其与实际电流的差异来表征机器人关节的剩余使用寿命。

         

    图12. 基于PICAIU框架的电流估计

    从192种工况采集的数据按不同工况变异因素划分为三个实验组。每组实验控制某一特定工况因素的变化,记录该工况变异实验中六个关节数据的拟合指数与波动方差。图13显示,纯动作变化工况获得的拟合指数显著低于其他工况组合,表明动作变化引起的信号估计误差大于其他工况因素。

         

    图13. 六关节电流估计效果

    5.3 减速器加速老化数据集的剩余寿命预测

    减速器损伤会导致关节阻力增大,并引起电流幅值上升。图14展示了150天加速寿命试验的全过程数据,主要包括关节3、4、5的电流变化。末端关节6及前两个关节1、2在加速实验期间处于固定状态。图中显示的负值表示电流方向与额定方向相反。

         

    14. 减速器失效过程的电流变化

    随时间推移,电流整体下移趋势表明关节阻力绝对值正在增大。在干摩擦条件下,减速器若持续摆动特定角度将导致失效,因此以当日关节摆动角度作为剩余寿命的参考基准。本实验完整采集了机器人减速器从初始安装到完全失效的全过程数据,并以当日关节摆动角度作为退化参考指标。        

    退化预测实验旨在检验所设计指标与真实工作周期数的契合程度。图15横轴表示采样数据组编号,编号越大表示采样时间越接近完全失效时刻。每完成一次加工轨迹,减速器运行周期数将增加整数倍。为便于分析,按完整轨迹执行次数对数据进行分组,因此实际剩余寿命在图1516中呈现为线性模型。经归一化处理后得到图15所示效果,各方法得到的RUL曲线越接近该直线,预测精度越高。

         

    15. 完整实验过程的退化曲线拟合效果        

    通过将加速寿命试验中的运动数据输入PICAIU网络,可获得健康状态下的扭矩曲线。利用改进KL散度计算实际电流与理论电流的偏差,并将其与其它常见特征量共同呈现在图15中,可发现本文提出的健康指标与退化实验中关节转动角度的累积增量具有更好的一致性。        

    16通过绘制四种变工况下的真实退化曲线与预测曲线,对比了本文提出的退化指标与其他方法的性能。可以看出,常见时域特征难以在柔性工况下实现与退化曲线的准确映射。从图16可见,当轨迹发生变化时,本文方法所提RUL的预测精度下降幅度最小。        

         

    16. 控制变量后的变工况实验结果        

    5.4 PICAIU鲁棒性测试        

    为评估运动姿态对预测精度的影响,本研究提取了不同运动类型的退化过程数据进行剩余寿命拟合测试。如表3所示,根据执行机构运动方向与重力方向的关系将测试数据分为三类。结果表明,当执行机构沿重力方向下移时拟合精度最高,基本达到表4所示水平。

    不同姿态下的预测精度

         

    剩余寿命预测的基线模型对比

         

    4将本文方法与近期剩余寿命预测研究进行了对比。在本研究采集的机器人减速器数据集上,表现最佳的基线模型为CNN-BiLSTMPICAIU的百分比预测误差RMSE相对该模型降低26%,绝对误差范围减少4.63,表明本方法的预测性能更加稳定。此外,使误差低于10%所需的训练周期仅为PF-BiGRU-TSAM模型的42%,证明本方法具有更优的样本利用效率。

    跨设备泛化能力测试结果

         

    为验证方法的泛化精度,采用JAKA Zu-12Zu-16两款机器人进行交叉验证。分别使用两款机器人的老化数据训练减速器RUL预测模型后,将其交叉应用于另一数据集计算预测误差。由于Zu-16额定负载更高且其减速器抗退化特性更强,本研究在Zu-12上采集的数据量更大。泛化实验中的百分比误差RMSE分别为12.66.4,表明从Zu-12Zu-16的泛化效果显著优于反向操作,其他指标也呈现相同趋势。

    5.5 多工况退化实验        

    本节通过消融实验详细检验物理约束的各项特性通过控制物理约束的作用范围,系统评估PICAIU所采用各项改进措施对前文所示剩余寿命预测精度的贡献程度,同时揭示性能局限并分析潜在原因。        

    实验对本文提出的健康评估方法与现有可靠性方法进行了量化对比。鉴于工业机器人系统的高度复杂性,研究领域目前缺乏针对关节退化程度的直接参考指标。因此,实验将本文提出的评估指标与航空发动机、轴承等领域的相关方法进行跨学科对比。对于这些跨领域健康评估方法,实验将其直接应用于本文关注的关节电流退化过程,最终结果如表6所示。

    不同序列建模方法的误差对比        

         

    从结果可以看出:在重复工况下,本文方法与其他健康评估方法相比未显现显著优势;但在柔性工况下,基于动力学观测神经网络的本文方法展现出明显的稳定性与精度提升。实验还发现一个重要现象:尽管多种机器学习回归方法能准确拟合现有数据,但其泛化能力仍显不足。即使神经网络完美拟合了输入特征与工作周期曲线的关系,仍无法准确预测新工况下的设备退化。        

    为验证物理约束对深度学习框架的改进效果,实验采用多种监督回归模型测试了物理损失项的增强作用。表 结果显示:对于常见RUL预测模型,添加物理约束仅小幅提升精度,但显著减少所需训练轮次并降低泛化误差。        

    剩余寿命预测的消融实验        

         

    表 展示了消融实验结果。实验进一步比较了两类知识约束对关节电流逼近性能的影响:移除所有物理约束后,神经网络迭代速度减缓且最终精度下降;添加完整物理约束后,网络训练速度加快且最终精度显著提升;仅保留动力学约束时,网络收敛速度明显加快但最终收敛精度不及双约束条件;仅保留电机伺服约束时,网络收敛速度略有提升但最终精度改善有限。实验表明,来自不同领域的知识共同提升了数据驱动方法的效率与拟合精度,其中动力学约束起主导作用,电机伺服约束起辅助作用。        

    两种物理损失项的消融实验        

         

    最后,通过表 的消融实验验证了退化曲线拟合误差。保留两类知识约束获得的扭矩偏差与退化曲线达到最佳拟合效果,而仅保留单一约束的扭矩偏差拟合效果略有下降,未使用知识约束的扭矩偏差则呈现不理想的拟合效果。实验证明,纯数据驱动方法难以从有限健康运动数据中学习真正的关节退化敏感特征。即使获得良好的电流估计效果,面对新退化数据时仍无法实现有效泛化。这印证了本文采用的PINN方法的核心优势:机器人学知识确保神经网络能从有限真实数据中学习系统本质信息,并有效外推至新的状态,从而避免过拟合陷阱。        

    面向概率与物理约束的消融实验        

         

    5 结论

    本文通过物理知识嵌入,开发了一种用于机器人损伤评估的理论数字孪生框架。通过可解释的模式识别方法,可实现工业机器人的在线系统辨识。这种基于模型退化程度的可解释健康评估方法融合了数据驱动与知识驱动方法的优势,取得了传统上需要大规模数据集才能实现的效果。

    该方法的优点在于借助精细化的机器人学知识降低了数据驱动模型的盲目性;其局限性在于人类经验难以推广至系统参数不透明的在役设备。对于参数标定困难或设备信息未开放的硬件,本文方法尚无法部署。下一步将深入研究机器人加工过程中的激励力传递模型,完善物理约束细节的实现,识别不同加工工艺间的差异,并探索跨个体差异的衰减过程建模理论方法。

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

    来源:故障诊断与python学习
    ACTMechanicalSystem振动复合材料旋转机械电源通用航空海洋焊接理论电机化机材料多尺度数字孪生物流
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2025-11-09
    最近编辑:9月前
    故障诊断与python学习
    硕士 签名征集中
    获赞 87粉丝 152文章 328课程 0
    点赞
    收藏
    作者推荐

    转速效应消除与多项式归一化|风机轴承早期故障检测与寿命预测新指标

    本期聚焦风力发电机轴承健康监测的新突破:在变速工况下实现早期故障检测(incipient fault detection, IFD)与剩余使用寿命(remaining useful life, RUL)预测。现有健康指标(health indicators, HIs)大多依赖稳态条件,其结果难以准确刻画退化轨迹,因而在非稳态下的应用受到严重限制。为解决这一难题,本文提出了一种多项式转速归一化健康指标(polynomial speed normalized health indicator, PSN-HI):首先利用基于能量比的转速效应消除方法弱化转速干扰,随后通过遗传算法(GA)优化的多项式归一化函数进一步校正振动信号幅值,从而在本质上去除转速影响。在此基础上,构建结合时域与频域特征的健康指标,实现了退化过程的定量刻画。实验结果表明,PSN-HI 相较传统与最新 HI,不仅能更准确反映退化趋势,还在 IFD 与 RUL 预测精度上实现显著提升,为变速风机轴承的健康管理与可靠运行提供了全新的解决方案。 论文链接:通过点击本文左下角的阅读原文进行在线阅读及下载。论文基本信息论文题目:A polynomial speed normalized health indicator for both incipient fault detection and prognosis of variable-speed wind turbine bearings 论文期刊:Advanced Engineering Informatics论文日期:2025年论文链接:https://doi.org/10.1016/j.aei.2025.103455作者:Dingliang Chen (a), Yi Wang (a), Yi Chai (b), Yuejian Chen (c), Yi Qin (a*)机构:a: State Key Laboratory of Mechanical Trans mission, Chongqing University, Chongqing 400044, China;b: School of Automation, Chongqing University, Chongqing 400044, China;c: Department of Mechanical Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada团队带头人简介:秦毅老师,现任重庆大学机械与运载工程学院教授、博士生导师,高端装备机械传动全国重点实验室成员。国家级青年人才,首批“重庆英才·青年拔尖人才”,重庆市学术技术带头人,连续入选全球前2%顶尖科学家“终身科学影响力排行榜”。秦教授长期致力于智能故障诊断与运维、动力学建模与数字孪生、智能无人系统及智能制造等方向的研究,重点关注复杂装备状态感知、健康评估、故障诊断与寿命预测等技术在工程实践中的应用。在科研方面,秦教授主持国家自然科学基金、国家重点研发计划、重庆市自然科学基金及多项企业合作项目;在国内外高水平期刊和会议发表论文230余篇(包括中科院一区及TOP期刊论文98篇,被引用7700余次),出版学术专著1部,参编教材1部。(来源:https://faculty.cqu.edu.cn/YiQin/zh_CN/index.htm)目录1 摘要2 引言3 基于转速归一化的健康指标构建方法3.1 基于能量比的转速效应消除方法3.2 多项式转速归一化3.3 风机轴承健康指标构建3.4 风机轴承的早期故障检测与剩余寿命预测4 实验验证与分析4.1 风力发电机轴承数据集描述4.2 基于所提方法的健康指标构建4.3 基于 PSN-HI 的风力机轴承早期故障检测4.4 基于 PSN-HI 的风力机轴承剩余寿命预测5 结论1 摘要现有的健康指标( health indicators, HIs)大多建立在稳态条件下,其结果往往难以真实反映退化轨迹,因而无法有效用于非稳态工况下的早期故障检测(IFD)和剩余寿命预测(RUL)。为此,本文提出了一种新的多项式转速归一化健康指标( polynomial speed normalized HI, PSN-HI),用于实现对变速风力机(WT)轴承的 IFD 与 RUL 预测。首先,提出了一种基于分解分量能量比与转速相关性的转速效应消除方法,以部分消除转速对监测数据的影响。随后,采用遗传算法与预设目标函数构建多项式转速归一化方法,对重构数据进行归一化,从而进一步消除转速对振动数据的干扰。在此基础上,通过两种混合模型生成时域与频域 HI,并通过线性加权构建 WT 轴承健康指标。最后,对比实验结果表明,PSN-HI 能够更为定量地刻画退化过程,并在变速条件下较传统与先进 HI 更好地实现 IFD 与 RUL 预测。关键词:故障检测,剩余使用寿命预测,健康指标,变速工况,风力发电机2 引言作为一种关键的机械装备,风力发电机(WT)能够将风能转化为电能,已广泛应用于陆上和海上风电场。然而,风机部件的故障率逐渐增加,导致运行与维护成本上升,同时降低了风电场的经济效益。因此,必须开展风机部件的早期故障检测(IFD)和剩余使用寿命(RUL)预测。健康指标(HI)的构建是实现 IFD 和RUL 预测任务的关键步骤,其目标在于定量评估风机部件的退化轨迹。合适的 HI 能够有效检测早期故障,并提高 RUL 预测的精度。近年来,研究人员在 HI 构建方法方面开展了大量探索。作为常用的基于统计参数的 HI,均方根值(RMS)和峭度(kurtosis)可分别从能量和冲击特性角度评估机械部件的退化过程。尽管基于统计参数的 HI 具有一定的物理意义,但它们往往仅从单一统计特征出发,难以全面描述退化过程。为克服上述不足,研究人员广泛采用了数据降维算法(如主成分分析 PCA、核主成分分析 KPCA、等距映射 ISOMAP)来融合多个统计特征,进而生成综合性的 HI 用于退化评估。在这些方法中,通常先从监测数据中提取多个统计特征,再通过特征融合构建单一 HI,以反映复杂的退化过程。尽管这类方法能够较好地刻画退化过程,但性能仍可能受冗余特征影响,同时还需依赖先验知识进行特征筛选。近年来,借助深度学习(DL)强大的特征挖掘能力,基于 DL 的 HI 构建方法受到广泛关注。根据构建方式,这类方法可分为有监督与无监督两类。然而,大多数基于 DL 的 HI 方法均假设工作条件恒定或近似恒定。在工业场景中,机电设备(如风力发电机、工业机器人、高速列车)常处于变工况下运行。此时,由于运行速度与载荷对监测数据的耦合作用,现有方法可能失效。通常情况下,运行速度可以测量,而载荷难以直接获取。由于载荷往往与速度相关,因此本研究主要关注变速条件下的 HI 构建,并在消除速度影响的同时抑制载荷的干扰。为减轻速度变化对监测信号的影响,已有研究开展了相关探索,速度归一化方法能够有效消除速度变化对信号幅值的影响。然而,现有方法多假设信号幅值与速度平方相关,但在实际工业场景中,这种关系可能并非严格成立,甚至可能与速度的高阶项相关。此外,速度变化还会导致信号能量比值发生变化。现有 HI 在变速条件下仍难以实现有效的 IFD 与 RUL 预测。针对上述问题,本文提出了一种适用于变速条件下风机轴承退化评估的 HI 构建方法。首先,提出基于能量比的速度效应消除方法,利用变分模态分解(VMD)将全寿命周期信号分解为多个分量,并通过能量比与速度的相关性分析,重构弱相关分量的信号。其次,提出基于多项式的速度归一化方法,以消除速度对监测信号幅值的影响,速度多项式的系数由遗传算法(GA)自适应确定。然后,采用高斯混合模型(GMM)、指数混合模型(EMM)和分布接触率度量(DCRM)在时域和频域中衡量监测数据与基准分布的差异,构建风机轴承 HI。所构建的 HI 能够在变速条件下有效表征退化趋势,并实现风机轴承的 IFD 与 RUL 预测。最后,通过风机发电机轴承数据集验证了所提出 HI 的有效性和优越性。本文的主要贡献如下:提出一种基于能量比的速度效应消除方法,利用 VMD 分解寿命周期信号,并通过能量比与速度相关性分析筛选弱相关分量以重构信号。 提出一种基于多项式的速度归一化方法,以消除速度对监测信号幅值的影响。速度多项式系数由 GA 自适应确定,归一化信号通过监测幅值与该速度多项式的比值计算获得。 提出一种基于分布差异的风机轴承 HI 构建方法,实验验证了其有效性和优越性。与其他 HI 相比,该方法在IFD 与 RUL 预测任务中展现出更强的故障检测与预测能力。3 基于转速归一化的健康指标构建方法在所提出的 HI 构建方法中,首先开发了一种基于能量比的转速效应消除方法,用以部分消除转速对监测数据的影响。随后,提出多项式转速归一化方法,以进一步去除转速对监测信号的影响。最终,通过 GMM、EMM 和 DCRM 在时域与频域中构建 HI。风机轴承的最终 HI 则通过线性加权时域与频域 HI 来实现健康监测。该方法的整体流程如图 1 所示。 图 1 所提 HI 构建方法整体流程示意图 3.1 基于能量比的转速效应消除方法 作为一种完全非递归的模态分解方法,VMD 可以克服经验模态分解方法中存在的端点效应与模态分量混叠问题,并且具备一定的数学理论基础。因此,本研究采用 VMD 对全寿命监测信号进行分解,以获得多个本征模态函数(IMFs)。风力发电机组轴承的全寿命数据假设为 。其中, 表示全寿命数据的数量, 表示单次采样数据的长度。VMD 的分解阶数 由能量差分法确定。寿命周期信号中的数据 通过 VMD 分解为模态分量的线性组合: 其中, 表示第 个 IMF 分量。通过 VMD 方法获得的多个 IMF 可以表征监测信号的不同组成部分。这些 IMF 的频率由低到高依次递增,因此每个 IMF 都能捕捉监测信号在特定频率范围内的振动模式或成分。为了分析不同分解信号与转速之间的关系,通过相关性分析选取与转速强相关的分解信号。每个 IMF 分量的能量计算公式如下: 其中, 表示第 个 IMF 分量的能量。随后,通过计算该分量在所有 IMF 分量中的能量占比来量化其相对重要性 。其公式如下: 接下来,通过全寿命周期数据 得到各 IMF 分量的能量比 , 。随后,利用皮尔逊相关系数分析每个分解信号的能量比 与寿命周期内的转速 之间的相关性: 将第 个 IMF 分量中相关性最高的分量作为与转速强相关的分解信号。其余分解信号则进行叠加,用以重构监测信号 : 最终,得到的寿命周期数据已在一定程度上削弱了转速的干扰,其表示为: 3.2 多项式转速归一化 基于能量变化的转速效应消除方法可以在一定程度上去除转速对监测数据的影响,但不能完全消除其作用。因此,本文提出一种新的多项式转速归一化方法,以进一步消除转速对监测信号的影响。该方法利用转速信号对原始监测信号进行除法处理,从而消除转速变化对监测数据幅值的影响。多项式各阶对应的系数通过遗传算法结合预设的目标函数求得,该目标函数定义如下:其中, 表示多项式转速归一化后的监测信号均方根(RMS)向量。速度多项式的具体阶数系数通过遗传算法(GA)最小化目标函数获得。监测信号通过速度多项式进行除法处理,从而得到归一化后的监测信号: 其中, 为速度多项式的系数, 为速度多项式的阶数, 为用于控制监测信号尺度的调节因子。最终,经多项式转速归一化处理后的监测数据表示为 3.3 风机轴承健康指标构建 首先,基于多项式转速归一化后的监测数据建立高斯混合模型(Gaussian mixture model, GMM),用于生成时域健康指标。GMM 中的高斯分量数量 通过贝叶斯信息准则(BIC)确定。GMM 参数通过 K-means 算法初始化,并使用期望最大化算法迭代更新。前 个健康数据被用作 GMM 的输入,用以估计基线分布 。监测数据 的数据分布 。通过 GMM 方法分别估计。随后,利用 与 计算分布接触率(DCRM)值 ,公式如下: 其中, 表示 与 的重叠区域面积,其通过积分公式计算如下: 通过计算得到的 DCRM 值生成时域健康指标向量 。时域 HI 可以在一定程度上刻画风机轴承的退化过程,但无法完全反映轴承的退化趋势。因此,引入频域数据以生成轴承的频域 HI。采用快速傅里叶变换方法对监测数据 进行频域转换。随后,构建指数混合模型(EMM),其指数分量数量 通过 BIC 确定。EMM 参数同样通过期望最大化算法迭代更新。前 个健康 FFT 数据用作 EMM 的输入,以估计基线分布 ,而整个 FFT 数据 的数据分布 通过构建的 EMM 分别估计。随后,通过以下公式计算 DCRM 值 : 其中, 表示基线分布 与 的重叠区域面积,其计算公式如下: 通过计算得到的 DCRM 值生成频域健康指标向量 。随后,通过对时域 HI 与频域 HI 进行线性加权,构建轴承的多项式转速归一化健康指标(PSN-HI)。该 PSN-HI 可结合时域与频域信息,用于描述风机轴承的退化过程。3.4 风机轴承初期故障检测与剩余使用寿命预测 采用 区间来确定报警阈值,从而准确检测风机轴承的初期故障点(IFP)。在确定 IFP 后,连续三个点中第一个超过报警阈值的点被视为真实 IFP;相反,没有连续超过三个点超出阈值的情况则被视为伪 IFP。作为长短期记忆网络的一种变体,MLMA-Net 模型已成功应用于 RUL 预测。MLMA-Net 通过记忆增强函数和分层机制提升 LSTM 的记忆能力。在本研究中,MLMA-Net 模型用于验证 PSN-HI 在 RUL 预测中的有效性。为了评估 RUL 预测的不确定性,在测试阶段引入了 dropout 操作以量化 MLMA-Net 的不确定性。MLMA-Net 的单元结构如图 2 所示。 图 2 MLMA-Net 的单元结构示意图遗忘门 、输入门 、输出门 以及候选状态 的计算公式如下: 其中, 表示当前时刻的输入数据; 表示上一时刻的隐藏状态; 表示 Sigmoid 函数; 表示双曲正切函数。基于遗忘门和输入门的输出值,注意力分布 的计算公式如下: 其中, 表示记忆增强函数; 表示门控的注意力比例; 用于判断门控输出值的正负性; 与 为主注意力门; 表示中层注意力层次; 表示高层注意力层次; 表示低层注意力层次。MLMA-Net 的输出由隐藏状态 决定,计算公式如下所示: 其中, 为 MLMA-Net 模型的输出; 与 分别为模型的权重和偏置参数。在训练 MLMA-Net 时,采用均匀初始化方法对权重参数进行初始化,并使用随机梯度下降(SGD)方法更新这些参数。4 实验验证与分析4.1 风力发电机轴承数据集 风机发电机轴承数据集采自中国北方平原的一个大型风电场。该数据集的现场传感器布置如图 3 所示。加速度传感器安装在发电机的主动端和非主动端。监测数据由多通道数据采集单元收集并存储。 图 3 风机发电机轴承数据集的现场传感器布置采样频率为 25.6 kHz,采样持续时间为1.0 s,覆盖整个寿命周期。在本研究中,采用 #2 和#10 风机发电机主动端的监测数据,以验证所提出的 HI 构建方法在变工况下的性能。风机发电机轴承全寿命周期数据的时域波形如图 4 所示,对应的转速曲线如图 5 所示。轴承寿命周期数据的样本长度分别为 #2 风机 204 个样本,#10风机 225 个样本。在维护过程中,发现轴承滚道存在严重磨损及元件滑移痕迹,磨损碎屑污染了润滑脂。风机发电机的缺陷轴承如图 6 所示。 图 4 风机发电机轴承全寿命周期数据的时域波形 图 5 风机发电机转速曲线 图 6 风机发电机的缺陷轴承 4.2 基于所提方法的健康指标构建 本文选用两台风机的发电机主动端轴承数据集,以验证所提 HI 构建方法的有效性。#2风机和 #10 风机的 VMD 分解数分别设为3 和 5。在确定分解数后,VMD 对寿命周期信号进行分解,并分析分解信号能量比与转速之间的相关性。随后,采用与转速相关性较弱的分量重构寿命周期信号,从而部分消除转速影响。基于重构后的监测信号,使用 GA 优化速度多项式系数。#2风机和 #10 风机所对应的GA种群规模与最大迭代次数分别设为 50、20、100和 100。确定速度多项式系数后,通过将监测信号幅值除以速度多项式获得归一化后的监测信号,从而进一步消除转速对信号幅值的影响。基于归一化监测数据,利用 GMM 和EMM 构建发电机主动端轴承 HI。选取健康状态下的前 50 个正常样本估计基线分布。构建的 HI 如图 7 所示。 图 7 构建的风机轴承健康指标根据图7,HI 值限制在 [0, 1] 范围内,可对变速工况下的风机轴承动态退化过程进行定量刻画。在正常阶段,HI 值较低,表明当前监测数据的分布与基线分布相似,轴承处于无缺陷状态。在退化阶段,HI 值明显高于正常阶段,显示当前监测数据分布与基线分布差异较大,表明轴承存在严重缺陷。由此可见,所构建的 HI 可有效、定量地描述风机轴承的退化过程。 为了对比所提出的 PSN-HI,还应用了多种传统与融合型健康指标。此外,还采用了基于 GMM 的HI(GMM-HI)、基于EMM 的 HI(EMM-HI)、以及基于置信值的GMM HI(GCV-HI)进行进一步对比。为了验证多项式速度归一化方法的优越性,还使用了未进行速度归一化的 HI(WSN-HI)以及平方速度归一化的HI(SSN-HI)进行对比分析。在WSN-HI 和 SSN-HI 中,两类HI 也是通过时间域与频率域 HI 的线性加权构建的。传统 HI 包括RMS、峭度(kurtosis)、Gini指数(GI)、负熵(NE)、近似熵(ApEn)、L2/L1范数和 Hoyer 指标(HM)。融合型HI 包括基于主成分分析的 HI(PCA-HI)、基于核主成分分析的HI(KPCA-HI)以及基于等距映射的HI(ISOMAP-HI),这些HI 通过 PCA、KPCA 和ISOMAP 对十个时间域与频率域特征进行构建。采用四个评价指标——单调性(Mon)、相关性(Cor)、鲁棒性(Rob)和综合指数(CI)——来评估不同方法构建的HI 性能。此外,还利用所提出的 PSN-HI 相对于其他 HI 的提升率(PR)来展示其优势。PR 定义如下: 不同 HI 的平均评价结果(基于WT 发电机轴承数据集)列于表 1。表 1 风电机组的不同健康指标(HI)平均评价结果 所提出的 PSN-HI 的综合指数(CI)最高,为0.531。此外,PSN-HI 的CI 分别比 RMS、峭度、GI、NE、ApEn、L2/L1范数、HM、PCA-HI、KPCA-HI、ISOMAP-HI、GMM-HI、EMM-HI、GCV-HI、WSN-HI及 SSN-HI 高17.0%、45.5%、60.4%、46.7%、74.7%、51.7%、50.9%、21.0%、23.2%、16.5%、45.5%、8.4%、52.6%、7.1%和 14.2%。PSN-HI的 Mon 值高于所有对比 HI。此外,Cor和 Rob 值也高于其他对比 HI,仅 EMM-HI 与 WSN-HI 除外。因此,PSN-HI 具有最佳的综合性能,能够定量描述 WT 发电机轴承在变速条件下的动态退化过程。这表明,基于能量比的速度效应消除方法和多项式速度归一化方法能够有效消除速度变化的影响,从而构建 WT 发电机轴承的HI。4.3 基于PSN-HI 的风力机轴承早期故障检测 在正常状态下,取前 50 个 HI 值用于计算 3σ 区间的均值和标准差 。采用所提出的 HI 对 #2 和 #10 风机进行 IFP 检测,结果如图 8 所示。IFP 能够被准确识别,其中 #2 和 #10 风机轴承的 IFP 分别位于位置 123 和 192。 图 8 基于所提出 HI 的初始故障点(IFP)检测结果 当轴的基频为 1 Hz 时,计算风机轴承的故障特征频率(FCFs)以验证检测结果,如表 2 所示。 表 2 WT 轴承的故障特征频率IFP 处的监测信号波形和频谱见图 9 和图10。在获得的频谱中明显观测到轴承的 FCF,验证了在检测位置确实存在初始故障,进一步证明了基于所提出 HI 的 IFD 精度。该验证表明,通过多项式速度归一化构建的 HI 在变速条件下用于风机轴承初始故障检测是可行且有效的。 图 9 #2 WT 第 123 位置监测数据的波形与频谱 图 10 #10 WT 第 192 位置监测数据的波形与频谱 对于 #2 风电机组,基于比较性HIs 的检测结果如图 11 所示。从图中可以直观地看到,基于 GI、NE、ApEn、L2/L1 Norm、HM、GMM、GCV和 SSN 构建的 HIs 在健康状态阶段均存在明显的数值波动现象。这种波动可能在健康状态监测过程中触发误报。基于 GI、NE、ApEn、L2/L1 Norm 和 HM 构建的 HIs 均未能成功检测出轴承的IFP。RMS、Kurtosis、PCA-HI、KPCA-HI、ISOMAP-HI、GMM-HI和 GCV-HI 检测出的 IFP 存在延迟现象。上述 HIs 对初始故障不敏感,这主要是由于速度变化引起的幅值效应造成的。SSN-HI 检测到的IFP 也存在延迟,这主要是因为平方归一化方法仅考虑速度平方,而未考虑速度的高阶影响。尽管 EMM-HI 和WSN-HI 检测到的 IFP 与 PSN-HI 相同,但在IFP 之后均表现出较大 波动,这不利于 RUL 预测。 图 11 #2 风电机组健康指标的 IFD 对比结果 对于 #10 风电机组,基于比较性HIs 的检测结果如图 12 所示。从图中可以观察到,大多数比较性 HIs 无法有效表征WT 轴承的退化过程。RMS、Kurtosis、GI、NE、ApEn、L2/L1 Norm、HM、PCA-HI、KPCA-HI、GMM-HI和 GCV-HI 构建的HIs 均未能检测出轴承 IFP。ISOMAP-HI、EMM-HI、WSN-HI和 SSN-HI 检测到的 IDP 分别为 223、192、192和 200,均发生在 PSN-HI 检测结果之后或与之相同。综上可知,所提出的 PSN-HI 拥有更强的 IFD 能力,并且比传统方法和现有先进 HIs 具有更高的 IFD 精度。 图 12 #10 风电机组健康指标的 IFD 对比结果 4.4. 基于 PSN-HI 的风电机组轴承RUL 预测 为 #2 和#10 风电机组构建的 HI 被用于通过所建立的 MLMA-Net 模型估算风电机组轴承的RUL。为了保证 RUL 预测结果的公平性和可靠性,进行了 100 次并行预测实验,并取这些预测结果的平均值作为风电机组轴承的最终 RUL。在预测阶段,启用 Monte-Carlo dropout 操作以量化预测模型的不确定性。#2 和 #10 风电机组的实际 RUL、预测 RUL、95%置信区间及 RUL 的概率分布如图 13 所示。 图 13 RUL 预测结果及不确定性分析。(a) #2 WT,(b) #10 WT 预测结果表明,MLMA-Net 模型预测的RUL 与实际 RUL 非常接近,并且预测值均落在 95% 置信区间内。根据RUL 的概率分布,可见概率分布的离散程度较小,说明预测不确定性较低。因此可以得出结论,PSN-HI 适用于RUL 预测,可基于 MLMA-Net 模型获得准确的RUL 预测结果。为了验证所提出的 PSN-HI 在RUL 预测中的优越性,还采用其他对比 HI 进行对比分析。各 HI 在RUL 预测中的评估结果列于表 4 和表 5。表 4 #2 WT 不同健康指标的预测性能评估 表 5 #10 WT 不同健康指标的预测性能评估 表 4 和表 5 显示,所提出的健康指标(HI)具有最低的预测误差和最高的评分(Score)值。对于#2 风电机组,PSN-HI 的预测性能优于传统HI,例如 RMS、Kurtosis、GI、NE、L2/L1范数和 HM。同时,PSN-HI 的预测效果也优于PCA-HI、KPCA-HI、ISOMAP-HI、GMM-HI、EMM-HI、GCV-HI、WSN-HI及 SSN-HI。对于#10 风电机组,所提出的 PSN-HI 优于所有对比方法。上述结果表明,PSN-HI 更适用于风电机组轴承的 RUL 预测,并能够提高RUL 预测的准确性。可以得出结论:基于能量比的速度效应消除方法与多项式速度归一化方法能够有效消除监测数据中速度的影响,从而构建性能优异的 HI。此外,所构建的PSN-HI 在预测性能上优于其他传统和先进的 HI 方法。5 结论本文提出了一种新的多项式速度归一化健康指标(PSN-HI),用于实现风电机组轴承在可变转速条件下的早期故障检测(IFD)与 RUL 预测。首先,提出了基于能量比的速度效应消除方法,该方法利用 VMD 将原始生命周期信号分解为多个模态分量,并分析各分量能量比与转速的相关性,以重构监测数据,部分消除速度对振动信号的影响。基于重构的生命周期信号,提出了多项式速度归一化方法,通过遗传算法(GA)优化与速度相关的目标函数获取多项式系数,并通过将重构信号除以速度多项式得到归一化的振动信号。随后,利用高斯混合模型(GMM)和指数混合模型(EMM)分别生成时域与频域健康指标,并通过线性加权构建最终的风电机组轴承健康指标。最后,通过风电机组发电机轴承数据集验证了所提出 PSN-HI 的有效性与优越性。同时,PSN-HI 被用于 IFD 与 RUL 预测。实验结果表明,PSN-HI 能够刻画非线性退化过程,其性能优于传统和现有先进健康指标,并且在可变转速条件下适用于 IFD 与 RUL 预测。与传统HI指标相比,PSN-HI 可获得更精确的故障检测与寿命预测结果。编辑:Jin校核:李正平、陈凯歌、赵栓栓、曹希铭、赵学功、白亮、任超、海洋、Tina、陈宇航、陈莹洁、肖鑫鑫、赵诚该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

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