首页/文章/ 详情

开源代码推荐 | 先验知识嵌入的首层可解释范式:面向轨道交通车辆人机协同故障监测

5月前浏览888

    近日,北京交通大学史红梅教授研究团队与香港理工大学和兰州理工大学团队在中科院计算机科学一区期刊《Journal of Industrial Information Integration》(IF=11.6)上发表题为《Prior knowledge-embedded first-layer interpretable paradigm for rail transit vehicle human-computer collaboration fault monitoring》的论文,提出一种基于可解释学习与多视图学习的故障诊断新范式——首层可解释范式。该研究通过嵌入信号处理先验知识,构建轻量化且可解释的胶囊网络PIFCapsule,并通过信号处理正则化损失优化整个网络结构,显著提升轨道交通车辆(高速列车、重载货车轮对、地铁列车)行走部轴箱(轴端)轴承的故障识别能力。


    代码链接:https://github.com/liguge/PIFCapsule

    若上述链接进不去,文末有下载方法。

    论文基本信息

    论文题目:Prior knowledge-embedded first-layer interpretable paradigm for rail transit vehicle human–computer collaboration fault monitoring

    论文期刊:Journal of Industrial Information Integration

    论文日期:2026

    论文链接https://doi.org/10.1016/j.jii.2026.101068  

    作者:Chao He[a,b], Hongmei Shi [a,b ,∗], Jing-Xiao Liao [c], Bin Liu [d], Qiuhai Liu [a,b], Jianbo Li [a,b],Zujun Yu [a,b]

    机构:

    a State Key Laboratory of Advanced Rail Autonomous Operation, Beijing Jiaotong University, Beijing 100044, China  

    b School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China  

    c Department of Industrial and Systems Engineering, The Hong Kong Polytechnic University, Hung Hom 999077, Hong Kong Special Administrative Region  

    d School of Mechanical and Electronical Engineering, Lanzhou University of Technology, Lanzhou 730050, China

    作者简介:

    何超,1995年出生于淄博市。目前正在北京交通大学攻读载运工具与运用工程博士学位。主要研究方向为迁移学习和工业大数据。目前已撰写并发表SCI和EI论文8篇,其中ESI高被引论文4篇。TIM/MST杰出审稿人,个人主页:https://liguge.github.io/  

    目录

    摘要

    1 引言

    2 相关工作

        2.1 基于物理的多视图信息的定义

        2.2 故障诊断中的人机协同智能范式

        2.3 由定制化制造引导的首层可解释瓶颈

        2.4 故障诊断中的胶囊网络

        2.5 基于物理的损失函数

        2.6 基于物理的特征融合

    3 提出的方法:PIFCapsule

        3.1 PIFCapsule概述

        3.2 具有信号处理先验知识的多重基于物理的可解释层

        3.3 注意力融合路由机制

        3.4 噪声阈值幅度比

        3.5 损失函数

    4 在轨道交通车辆轴承中的应用

        4.1 数据集和任务描述

        4.2 案例分析

    5 PIFCapsule的可解释性分析

        5.1 多视图首层瓶颈的可解释性

        5.2 路由系数Z的可解释性

    6 结论

    摘要

    轨道交通车辆承受着大载荷、高速度和恶劣环境的影响,易导致部件故障。首层可解释范式(FLIP)将人类先验知识嵌入智能设备,是受定制化制造和具身智能指导的智能范式之一。它由首层可解释模块、骨干网络和损失度量组成。然而,现有研究存在依赖单一来源信息、缺乏可解释骨干网络、无法进行特征融合等问题,因此难以处理多激励、耦合信号。为填补这一空白,本文提出了一种基于FLIP的单阶段多视图胶囊融合网络(PIFCapsule)。首先,设计了一个具有信号处理先验的第一层可解释模块,以实现参数自动优化,并突出不同信号处理算法产生的多视图特征之间的互补性。其次,采用可解释胶囊网络作为骨干网络。为解决信息融合效率低下和不足的问题,本文提出了一种高效的注意力融合路由(AFR),与传统的基于胶囊的网络相比,该路由将参数减少了约5.72倍,复杂度降低了约2.93倍。针对训练过程中缺乏基于物理的约束这一问题,本文提出了噪声阈值振幅比(NTAR)作为正则化项,通过抑制学习到的噪声来增强微弱的周期性瞬态脉冲。通过三个真实的轨道交通车辆数据集验证了该方法的有效性和可靠性:在仅使用10个样本的情况下,PIFCapsule的准确率比最先进的方法高出6.77%。鉴于其轻量化特性,该网络在智能边缘设备上的部署具有巨大潜力。代码可在https://github.com/liguge/PIFCapsule获取。  


     

    关键词:具身智能,多视图学习,首层可解释范式,注意力融合路由,噪声阈值幅度比,轨道交通车辆

    1 引言

    铁路运输安全与日常工作和生产活动密切相关。本研究关注高速、重载货运和地铁列车的监测与维护技术。车辆的走行部是一个至关重要的子系统。轴承作为走行部的关键部件,直接影响安全和可靠性。然而,它们面临高载荷、高速度、强噪声和复杂运行条件等挑战,导致点蚀、剥落和断裂等失效模式。

    具身智能(EI) 是实现轨道交通车辆智能运维的指导原则。EI将感知、认知和行动集成到制造系统中,实现设备上的"感知-认知-执行-反馈"闭环[1]。预测与健康管理(PHM) 平台也成为当前轨道交通车辆维护的标准设备。在PHM系统中,关键部件的故障预测和诊断是至关重要的环节。

    在故障预测和诊断算法中,一些适用于普通旋转机械的算法未能达到令人满意的性能。尽管面临这些挑战,智能故障诊断(IFD) 利用先进的信号处理和数据驱动算法,能够实现精确的故障识别和定位,从而提高铁路运输的安全性和可靠性[2-5]。然而,现有算法缺乏明确的可解释性,未能融入人类先验知识。这对实现定制化制造和具身智能(CMEI)[1]提出了重大挑战。

    具有人类先验知识的可解释IFD引起了广泛关注[13]。虽然计算机可以学习数据中的统计模式,但这是不可控的,可能无法与人类期望保持一致。内在可解释性旨在设计基于领域特定先验的透明且可理解的映射。因此,将人类先验知识嵌入系统,并通过可解释理论和物理定律约束优化过程[14]。

    首层可解释范式(FLIP) 代表了智能设备的典型智能范式。它将信号处理先验知识嵌入深度模型,以解决复杂的工业问题。如图1所示,FLIP是一种简单而有效的方法。FLIP被定义为使用先验知识修改或替换首层权重,或在骨干网络前部添加可微分信号处理层以模拟信号处理操作的框架。

    图1 首层可解释范式示意图  

    然而,根据时频不确定性原理(Gabor限制),时频分辨率难以充分表达弱故障特征。这是因为这些算法仅依赖于单一视图,这反映了信号的局部状态,容易受到信号失真和模糊性的影响。受多视图学习[18]启发,振动信号的特征可以从多种不同的信号处理技术进行描述,这些技术构成了原始信号的多个视图。显然,基于物理的多视图融合可以捕捉不同的特征,并为后续故障决策提供更可靠的依据[19,20]。  

    本文的主要贡献包括:   

    • 将多种散见于可解释故障诊断的多种方法(如基于小波、STFT、盲解卷积等的先验嵌入技术)系统性地总结并升华为一个统一的理论框架——首层可解释范式(First-Layer Interpretable Paradigm, FLIP)。该范式不仅清晰地刻画了以往研究的共同内核,更通过定义核心组件、设计原则和优化路径,为未来研究提供了可扩展的指导蓝图。

    • 在人机协同指导下,提出了一种单阶段基于物理的多视图融合胶囊网络(PIFCapsule),涵盖三个信号处理赋能的首层设计、一个注意力融合路由机制、基于物理的正则化项。    

    • 将注意力路由融合机制引入CapsNet,显著减少了原始动态路由的参数数量,同时提高了泛化能力。    

    • 引入噪声阈值幅度比作为基于物理的正则化项,可以使用可微分软阈值模块抑制噪声,而不依赖于故障先验周期。    

    • 通过高速列车、重载货运列车和地铁列车走行部数据集的综合结果证明了其卓越性能。

    2 相关工作

    2.1 基于物理的多视图信息的定义

    多视图特征融合是一种通过利用来自多个视图或特征的互补信息来提高性能的方法[18,32]。对于基于物理的多视图信息的定义,不同的研究者可能有不同的看法。这些视图可以是同一数据集的不同特征。在本文中,不同的信号处理技术被视为不同的视图。    

    尽管这三种信号处理方法都针对振动信号,但它们强调不同的物理特性:小波层使用小波权重和时域卷积从时域视角获取时频谱;STFT层使用窗权重和频域乘法获取频域表示;BD层不需要预设滤波器参数,仅使用NTRA优化的权重获取时域表示。    

    尽管由先验知识引导的自适应滤波权重可以实现时间和频率分辨率之间的相对平衡,但仍无法解决Gabor限制。为进一步缓解这一限制,不同的信号处理方法关注不同的物理方面,然后通过胶囊融合网络实现全局和局部特征融合。融合的特征是三个不同的物理视图,超出了原始振动信号单独所能提供的内容。    

    2.2 故障诊断中的人机协同智能范式

    在工业诊断领域,人机协作正在兴起,值得在循环具身智能制造的总体框架内重新审视。其中,智能设备的智能范式——包括预测与健康管理和人机协作等关键能力——作为定制化制造和具身智能的核心组成部分,为故障诊断提供系统支持。文献[33,34]主要在决策过程中添加人为干预。如果算法给出的故障置信度较低,人类需要基于其他指标和经验做出最终决策。Wen等人[35]引入了故障因果知识图谱,并使用大型语言模型实现人机协同决策。Li等人[36]也借鉴了这一思路,将可疑结果发送给人类进行二次判断。    

    上述方案仅在决策阶段引入人类知识以做出最终决策[19],而忽视了人类知识对计算机决策的指导作用。正如Kim所说[37],人类擅长从经验中发展先验知识,而机器擅长计算数据。因此,本文从不同角度解释人机协同诊断,将嵌入信号处理先验知识的可解释首层范式归类为典型的人机协作范式之一,并解决了FLIP的局限性。    

    2.3 由定制化制造引导的首层可解释瓶颈

    提出的首层可解释范式属于更广泛的信号处理赋能技术类别。它使用信号处理修改和优化深度模型,从而增强可解释性。受此启发,WaveletKernelNet[15]将连续小波变换集成到卷积中。此后,关于小波核的研究引起了广泛关注[7,38,39]。SincNet[40]强调了首层模块对提取特征的重要性。基于物理的卷积神经网络(PCNN)[41]基于轴承速度和故障特征频率生成卷积核。然而,PCNN只能识别单一故障。M-PINet[42]结合多个PCNN分支进行多故障识别,但计算复杂度高。M-IPISincNet[43]利用带通滤波器和逆傅里叶变换设计卷积以提取多尺度信息,但仅依赖于单一的基于物理的信息源。PICNN[44]基于不同特征频率对应不同故障设计物理特征加权层。它为某个故障的特征频率分配更高的权重,并减少其他频率的权重以抵抗干扰。Liu等人[45]将有限脉冲响应滤波器视为核,将中心频率和带宽视为频率的多项式,开发了首层可解释模型[46]。乘法卷积网络(MCN)[47]由一系列乘法滤波核组成,从频谱样本中提取差分模式。时频网络[48]利用STFT、Chirplet和Morlet小波模拟时频变换,获取丰富的特征[49]。尽管随机卷积层[50]不强调可解释性,但它随机改变核大小以生成新数据并保留全局信息,同时提高样本多样性和模型泛化能力。包络谱神经网络[51]根据集成经验模态分解指导首层模块设计,提取物理特征。尽管上述方法取得了有希望的结果,但它们本质上受到不可靠的单信息源的限制,容易受到不确定性的影响。这种限制导致分析不完整,进而降低了模型的鲁棒性。    

    2.4 故障诊断中的胶囊网络

    CapsNet在IFD中有广泛应用[52]。一些报告[53]使用谱图作为输入,并从事后解释中描绘热图。然而,这种方法计算成本高且参数量大。相比之下,一维信号具有某些优势。CapsNet还具有处理小样本诊断的能力。例如,CapsFormer[54]将CapsNet与自注意力机制集成。CapsNet用于特征提取,而自注意力机制关注谱图的空间和顺序方面。Wang等人[55]通过知识通知的卷积初始化胶囊,并提出谱模板方法建立胶囊与故障类型之间的映射。总体而言,大多数研究仍关注CapsNet在特征提取中的关键作用,例如使用CapsNet输出的特征进行视图自适应以计算源域和目标域之间的相似性度量[56-58],或将其与其他神经网络结合[59-62]。然而,高参数动态路由机制缺乏多源特征融合能力。此外,忽视基于物理的信息不仅降低了可解释性,还增加了对数据质量和数量的依赖。依赖单一数据模态或视图通常会产生模糊或不可靠的诊断决策,最终削弱泛化性能。    

    2.5 基于物理的损失函数

    基于物理的损失是基于物理的机器学习中的关键组成部分[63-65]。pyDSN[9]引入了平衡谱质量度量来评估STFT谱图质量。Qin等人[66]使用边界损失约束动态模型参数,从而最小化真实世界和模拟数据之间的频率差异。ClassBD[67]利用峰度和l2/l4范数优化盲解卷积滤波器。PGNN[68]将转子的不平衡振幅和相位角视为损失函数。Jia等人[69]提取状态特征频率和信号谱能量以推导物理伪标签并制定基于物理的损失。总之,基于物理的度量通常针对特定的信号处理流程或设备动力学进行定制。因此,基于物理的损失应针对不同的系统和算法进行定制。    

    2.6 基于物理的特征融合

    多源特征融合在IFD中得到了广泛研究[31,70]。然而,多源基于物理的特征融合仍然是一个新兴领域,值得注意的工作对此进行了研究。MPINet[42]结合多个基于物理的块(PIBs),并通过多分支特征连接实现融合。然而,这种方法需要训练多个PIBs。类似地,M-IPISincNet[43]集成电流和振动信号,使用逆傅里叶变换和带通滤波器设计基于物理的卷积层,并通过多尺度卷积提取多尺度特征。然而,MPINet和M-IPISincNet中的融合方法相对简单,缺乏互补信息视图之间的充分交互,仅考虑单一的基于物理的信息。Sun等人[71]引入了一种基于物理的多模态特征融合网络(PMFN),结合声学或振动信号与红外图像。PMFN集成了两种注意力机制:频域注意力,强调特定故障下的调制关系;区域注意力,突出谱图中更显著的特征,这些特征通过自注意力融合。然而,二维编码依赖于信号处理,其有效性严重依赖于手动参数选择。此外,通过两个自注意力分支的融合过程计算复杂度高,PMFN也忽视了基于物理的损失。Ying等人[72]提出了一种基于Dirichlet分布的可信融合模块,以确保动态融合的可靠性。然而,TMFEFN缺乏物理信息,需要大量训练数据。PFCG-Transformer[73]提出了一种集成声学和振动信号的物理框架,通过比较物理模型的数据分布与实际运行信号构建一致性和不确定性量化损失。然而,声振融合需要大量手动经验并涉及复杂计算,收集声学信号具有挑战性,图像信号处理高度复杂。    

    在实际轴承诊断中,振动信号对故障最敏感,并能充分反映故障特征[27,74]。PIFCapsule将信号处理赋能的权重视为可训练的,并通过数据驱动算法找到它们,而不是依赖于手动参数。此外,PIFCapsule通过注意力融合路由机制实现基于物理的多视图特征融合,使用范数激活操作关注信息最丰富的特征,并使用基于物理的正则化约束优化方向。

    3 提出的方法:PIFcapsule

    3.1 PIFcapsule概述

    图2 PIFCapsule的总体情况

    先前的首层可解释范式仅限于单视图信息,导致信息交互效率低下,并忽视了基于物理正则化的好处。因此,本文重点研究在首层可解释范式框架内基于CapsNet的基于物理的多视图特征融合。如图2所示,我们提出了一种注意力融合路由机制,增强了跨视图信息交互效率并保留了判别性特征。同时,我们实现胶囊路由,AFR促进了从初级胶囊到数字胶囊的高效信息流。此外,将基于物理的正则化纳入优化过程,以增强融合表示的物理保真度,产生更具泛化性的故障特征嵌入。    

     

    系统诊断流程    

    具体框架如图3所示。诊断工作流程包括以下步骤:    

    ● 收集三种类型的轨道交通车辆损伤数据集:高速列车牵引电机轴承数据集、重载货运列车轮对轴承数据集和铁路列车转向架数据集。收集的数据集分为训练集、验证集和测试集,每个类别的训练集大小不超过50个样本,符合小样本问题的定义。    

    ● 收集的训练集数据输入PIFCapsule,计算所有基于物理的损失以指导模型训练。训练后保存最优模型权重。    

    ● 在测试阶段,加载最优权重并评估性能。

    3.2 具有信号处理先验知识的多重基于物理的可解释层

    单源信息可靠性差且决策模糊,而多源特征融合可以捕捉具有抗干扰能力和鲁棒性的特征。因此,多源特征融合是一个优秀的解决方案。然而,多源特征融合仍然受到可解释性差的限制。因此,我们提出了一种基于物理的多视图特征融合框架。该方案集成了STFT、WT和BD信息视图:BD有效提取轴承故障的瞬态脉冲,而STFT和WT在这方面的能力较弱。STFT由于固定窗口而适应性和分辨率有限。WT提供更高分辨率并支持多尺度分析,但可能引入混叠和失真。

     3.2.1. 双阻尼拉普拉斯小波权重初始化 

    在小波源信息背景下,采用双阻尼拉普拉斯小波,这在小波权重初始化领域尚未被探索。实际振动信号中观察到的故障脉冲响应经常表现为双边不对称衰减波形,这一特性与双阻尼拉普拉斯小波的内在属性完全吻合[75]。

    其中         是采样频率,    是时间,    是双阻尼拉普拉斯小波的系数,    是时间常数,     和     是阻尼比。 因此,最终的小波视图层可以表示为:

       

    其中     是尺度因子,    是平移因子。     表示卷积操作。

    3.2.2. 窗权重初始化

    先前的研究表明,权重初始化方法在性能上优于信号处理层的设计[7]。信号处理层侧重于将关键参数转换为可微分形式。然而,这种方法存在大量计算开销,因为滑动窗口操作未被卷积替代。受权重初始化的启发,我们将权重的通道视为时域轴,最后一个维度视为频域轴。STFT核用作CNN的初始化权重,以有效模拟STFT。引入窗权重初始化:

    假设输入信号为     ,采样频率为     ,滤波器的输出通道数为     ,滤波器长度为     。定义的带宽          。对于每个通道     ,在范围      内设置滤波器      ,每个滤波器      使用窗函数(如Blackman窗)设计。通过归一化截止频率并考虑滤波器长度,计算理想脉冲响应     ,FIR滤波器响应可以表示为     与窗函数    的乘积。为确保初始权重更接近随机初始化,采用标准化来限制滤波器的范围。

       

    该操作的原理是每个滤波器捕捉原始信号的频率分量,将通道视为沿时间轴的滑动窗口捕捉这些频率分量的时间变化以模拟STFT:

       

    在STFT中,对加窗时间段执行傅里叶变换。为模拟此过程,我们使用傅里叶卷积而不是普通卷积:

       

    3.2.3. 盲卷积模块

    在BD领域,时域分量通过二次卷积神经网络(QCNN)  执行,而频域分量通过FFT处理。此过程通过集成峰度损失和     范数进一步增强。相比之下,我们设计噪声阈值幅度比,用于量化BD中的基本信息。 相比之下,QCNN产生大量参数且难以优化[67]。因此,我们选择标准CNN。类似地,我们通过减少时域中的周期性噪声并消除多余的频率分量来识别共振带。

       

    特别地,     的滤波器也在峰度损失和       范数度量的约束下学习。

    3.3 注意力融合路由机制

     

    图4 注意力融合路由机制 

    图4所示,第    层的胶囊由第     层的所有胶囊     预测,这是通过权重矩阵    和路由系数      实现的。 鉴于三个信息源具有不同的物理特性,使用单一的权重矩阵    不足以提取高度互补的特征。因此,利用三个线性层完成此变换。在这三个信息源中,通道表示时间轴,最后一个维度表示频率轴,它们共同映射到新空间:

       

    其中FC指全连接层。 随后,通过范数线性注意力机制实现特征融合。 自注意力机制通常使用Softmax激活函数,复杂度为     ,而线性自注意力机制复杂度为     。后者被认为是降低计算复杂度的有效方法,用设计的核函数近似相似性函数。在特征融合过程中,计算相似性有助于调整融合权重。给定具有特征表示     的核:

       

    其中     表示相似性函数。 根据线性注意力机制的定义,范数线性注意力可以表示如下:

       

    基于矩阵结合性质的性质,计算顺序可以从     改为     ,即:

       

    Softmax注意力机制提供非线性加权机制,从而容易关注最重要的特征。然而,线性注意力机制的输出更接近平均值,难以聚焦。因此,采用范数激活函数作为核函数来缓解这一缺点,即:

       

    其中     且      。    是高斯正态分布的累积分布函数。    表示范数,     用于控制聚焦程度,通常     。 与ReLU相比,GELU提供平滑性,有助于梯度优化,并对负输入保持非负性,从而防止神经元失活。然而,为满足分母的要求,输出必须保持非负,需要使用绝对值。最后,通过特征范数映射调节调整方向。 实际上,    通过自注意力机制获取三个信息源的全局融合特征,但也需要局部特征。为此,采用深度可分离卷积,可以理解为卷积注意力机制。它关注空间视图中的相邻局部信息,从而从每个局部区域获得不同的输出并增强特征多样性。如公式(12)所示:

       

    其中,    表示深度可分离卷积操作。 下一层胶囊的表示可以通过注意力权重获得:

       

    AFR有两个优点。首先,它在CapsNet中实现了路由机制,实现了非迭代、高速并行路由操作,减少了模型参数并提高了算法效率。其次,与传统注意力不同,AFR支持跨视图交互,并在模态内建立局部和全局依赖关系(长距离依赖),具有高效的信息适应性。它根据任务动态分配特征并关注关键信息,提高融合效率。 

    算法1 注意力融合路由机制

     

    3.4 噪声阈值幅度比

    盲卷积通过最大化信号的统计度量搜索滤波器来放大重复冲击。Fang等人提出了周期性噪声幅度比(PNAR) ,定义为噪声点的平均幅度与均方根的比值。如公式(14)所示:

       

    其中      表示欧几里得范数,     是信号长度,    是噪声长度。     可以解释为噪声位置的分布。 显然,      是一个关键参数,其值决定了PNAR的有效性。在PNAR中,噪声分布        是连续的,噪声长度由噪声比     和周期     决定,使得     。显然,连续噪声分布的假设过于理想化。实际上,由于环境干扰,噪声可能存在于整个周期,仅强度不同。类似地,预定义的噪声长度也不切实际。 为纠正这一限制,提出了噪声阈值幅度比(NTAR)。我们引入了一个新视角,将噪声位置视为噪声所占的权重——噪声系数。由于噪声分布在整个周期,噪声长度与信号长度相同(    )。假设噪声存在于整个周期,我们的目标是确定噪声强度与总信号强度的比例。因此,公式(14)可以重新表述为:

       

    NTAR主要通过噪声系数获取噪声,促进训练过程中噪声逐渐减少趋近于零。受小波系数启发,软阈值处理[76]也可用于减少噪声系数,其中阈值 至关重要。如图5所示,我们采用基于注意力的架构自动搜索      并保留噪声系数     ,如公式(16)所示:

           

    其中     是符号函数,    是阈值。保留噪声系数意味着保留所有低于阈值的值。 通过NTRA减少噪声比例以提高信噪比,从而提取更鲁棒的故障特征。

     

    图5 抗噪软阈值模块

    3.5 损失函数

    类别     的边界损失     如公式(18)所示。    是指示变量,如果类别     存在则     ,否则     。下界    ,上界     ,    表示输出的范数,     是降权因子,通常设置为0.5。边界损失旨在为每个类别输出真实的标签分布,同时抑制不存在类别的激活。

       

    其次,原始CapsNet包含重构损失以衡量生成信号与原始信号之间的相似性。然而,我们的目标是提高信噪比(SNR),而不是重构低SNR的原始信号。因此,采用    范数用于评估稀疏性,这有助于构建高SNR信号,有效抑制噪声,从而增强鲁棒性。如公式(19)所示:

       

    其中     ,    ,且     。 此外,小波和STFT都会生成谱图。在我们之前的工作中,提出了平衡谱质量(BSQ) [9]来平衡时频视图的信息并增强能量集中度。如第1节所述,小波和STFT信息是互补的,对信号的不同分量具有不同的敏感性,从而产生更敏感的能量集中区域。因此,使用BSQ来评估谱图的质量,如公式(20)所示:

       

    其中       和      分别是频率和时间视图的质量系数。 随后,采用    确定BD滤波器的最优参数。此过程旨在放大重复特征同时减轻噪声。最终,综合损失函数如公式(21)所示:

       

    确定     至关重要。为解决这一挑战,我们利用不确定性感知加权损失[67]来自动平衡每个损失分量对学习问题的重要性,并自动搜索特定的参数配置。

    4 在轨道交通车辆轴承中的应用

    4.1 数据集和任务描述

    图6 设备: (a)高速列车牵引电机; (b)超重货物列车轮对; (c)地铁车辆转向架齿轮箱

    本研究涵盖三个轨道交通车辆数据集,即高速列车牵引电机数据集、重载货运列车轮对轴承数据集和地铁列车转向架齿轮箱数据集。这些数据集的简要概述(表1)如下:

    轨道交通车辆的数据描述

     

    BJTU1[9]:数据集来自NTN牵引电机轴承试验台,这是一种具有国际先进标准的专用设备,用于进行高速列车牵引电机轴承实验,如图6(a)所示。该设备涵盖列车驱动系统的动态数据和信息收集与分析、地面平台实验验证以及全生命周期。它是实现高速列车牵引电机轴承预测与健康管理的特定实验平台。测试平台由电气控制柜、加速度计、多个测试轴承和信号传输设备组成。通过传感器和信号传输系统收集轴承的振动信号。为与实际应用场景保持一致,实验使用的轴承型号为HRB NU214 EM 32214H,其尺寸与实际应用中使用的轴承相同。包括四种机械状态:内圈故障(IF)、外圈故障(OF)、滚动体故障(BF)和正常(NC),采样频率为100 kHz,速度分别为200、250、300和350 km/h。    

    BJTU2[77]:数据集来自中车青岛四方开发的重载货运列车轮对轴承平台,如图6(b)所示。该试验台设计用于各种轴承型号的故障模拟、性能测试和剩余使用寿命评估。它可以模拟车辆轴承在实际运行期间的载荷、速度和相应的环境条件,从而最大限度地与真实铁路场景保持一致。具体而言,垂直载荷用于模拟轴重;横向载荷用于模拟不同的轨道条件,例如转弯、通过道岔和上下坡行驶时轨道上的侧向力;风扇用于模拟运行期间的横风,风速设置为8 m/s至10 m/s。实验样本是通过拆卸检查确认的故障轴承——这些轴承最初要么由在役铁路货车在线监测系统的警报识别为异常,要么通过人工检查识别。具体而言,标准轴承安装在轮对的一侧,实验轴承安装在另一侧。在轴箱外安装振动加速度传感器,采样频率为16 kHz。测试轴承为352226X2-2RZ。数据集包含七种故障类型:内圈剥落(IRP)、外圈剥落(ORP)、滚动体剥落(REP)、滚动体裂纹(REC)、保持架断裂(CF)、保持架半裂纹(CHC)和复合故障(外圈剥落+滚动体剥落)(CFs)。此外,还包括健康轴承(H)的数据。轴承转速设置在60 km/h至180 km/h范围内。

    4.1.2 地铁列车场景的公共数据描述

    BJTU3[78]:如图6(c)所示,实验平台基于地铁列车的转向架设计,与实际转向架的比例为1:2。此外,实验台上的传感器部署高度再现了实际测量点。采集的信号包含21个通道,采样频率为64 kHz,考虑了六种工况,转速设置为20 Hz、40 Hz和60 Hz,横向载荷为0 kN和10 kN。不同的速度用于模拟各种列车速度,而不同的横向载荷用于模拟直线行驶和转弯条件。评估列车传动系统的齿轮箱在九种不同健康状态下的表现:正常、齿轮裂纹齿(GCT)、齿轮磨损齿(GWT)、齿轮缺失齿(GMT)、齿轮崩齿(GCPT)、轴承内圈故障(BIR)、轴承外圈故障(BOR)、轴承保持架故障(BC)、轴承滚动体故障(BFE)。

    4.1.3 任务描述

    轨道交通车辆中的旋转机械设备面临几个挑战,包括重载荷、高速度、强噪声和复杂运行条件,再加上小样本。这种稀缺性通常阻碍训练出足够泛化的模型。为解决此问题,利用多源基于物理的特征融合CapsNet,特别是在有限样本下。本研究一方面涵盖三种不同类型的轨道交通车辆;另一方面,处理各种故障类型和设备,包括单点故障、复合故障以及齿轮箱和轴箱。    

    使用滑动窗口将信号分割成2048个样本的片段。滑动窗口之间没有重叠,以确保不会发生数据泄漏。根据先前的定义,如果每个类别的样本少于50个,则该问题被归类为小样本问题[6,79]。在小样本下,所提出的方法足以训练出泛化模型,从而在测试集上实现有效的故障识别。相比之下,其他主流方法在这方面表现相对较弱。训练样本大小为5、10、20、30、40和50。不同数据集和实验的训练样本设置各不相同,但都符合小样本的定义。测试集统一设置为每个故障1000个样本。    

    实验程序在PyTorch 2.5上执行,利用两个GTX 4090 24 GB GPU。每个实验重复五次以获得均值和标准差。参数配置详见表2    

    表2 配置参数    

     
    4.2 案例分析  

    4.2.1 在不同工况下不同数据集的PIFCapsule性能

     

    图7 在不同工作条件下使用小样本的性能能力    

    图7所示,评估了PIFCapsule模型在不同数据集和运行条件下的鲁棒性。具体而言,Data H的每个类别样本数设置为5,Data F和G设置为10。对于车辆,转速越慢,准确率越低。与其他条件相比,在低速和高载荷场景下,信噪比较低,噪声干扰在小样本下的影响更显著。(后续评估以这种最具挑战性的条件为基准。)即使在这样条件下,借助基于物理的信息和注意力融合路由机制等模块,PIFCapsule实现了超过90%的准确率,表现出鲁棒的性能.

    4.2.2 不同样本大小下PIFCapsule的性能评估

     

    图8 PIFCapsule在小样本下的性能能力

    因此,我们评估了PIFCapsule在解决小样本挑战方面的能力。如图8所示,性能随着样本大小的增加而稳步提高。具体而言,对于高速列车牵引电机的数据集,当样本大小达到10时,PIFCapsule达到100%。这种卓越性能可能归因于高SNR。相比之下,对于货运列车轮对轴承或地铁转向架,SNR相对较低,导致信号质量较差。由于性能取决于数据的质量和数量,因此必须设计在小样本下增强数据/特征质量的方法。为此,PIFCapsule从多个视图提取故障特征频率,从而进一步减少原始信号中的噪声并识别主要频带。单一信息视图可能无法完全捕捉特征频率。因此,小波视图、盲解卷积视图和STFT视图相互补充和加强,分别提取多尺度高分辨率特征、周期性脉冲分量以及不易混叠且保真度高的特征。随后,这些特征通过注意力路由特征融合和基于物理的损失函数充分整合,使模型能够学习更具判别性和高质量的特征。因此,即使只有10个样本,PIFCapsule也分别达到91.54%和95.62%。    

    4.2.3 不同路由机制的性能

     

    图9 三个模型的性能    

    表3 三种模型的性能、参数和计算复杂度对比

     

    PIFCapsule是一种基于胶囊的架构,其性能与原始CapsNet和EfficientCapsule进行了比较。EfficientCapsule使用自注意力机制替代原始动态路由进行视觉处理。如表3图9所示,PIFCapsule提高了27.35%和24.21%。值得注意的是,对于数据集F,这两个模型低于50%。然而,通过利用基于物理的多源信息和注意力融合路由机制,PIFCapsule有效减轻了噪声并提取了有价值的信息,性能超过90%。此外,PIFCapsule显著减少了参数数量和计算复杂度。它将参数减少了82%,FLOPs减少了66%,分别实现了0.4824 MB的参数大小和22.73 MFLOPs的计算成本。    

    注意力融合路由用于计算路由系数,从而减少迭代过程中的重复参数更新。此外,使用范数激活函数优化注意力权重,减轻了Softmax的高参数数量和复杂度。多视图全局特征融合直接在注意力融合路由模块中实现,无需复杂的多分支卷积融合结构。在注意力融合路由框架内实现深度可分离卷积以实现局部特征融合,降低了传统卷积操作固有的计算复杂度。    

     

    图10 混淆矩阵和t-sne    

    此外,通过混淆矩阵和T-SNE(如图10所示),我们分析了三个模型之间的性能差异。EfficientCapsule对于断齿仅达到68%,与轴承保持架故障混淆。原始CapsNet未能完全识别正常和异常,对于齿面磨损和轴承保持架故障的准确率低于62%。相反,PIFCapsule整合了基于物理的多视图信息,并通过AFR成功融合了这些特征。为精确提取鲁棒的物理特征,采用基于物理的损失来约束优化方向。

    4.2.4 与其他最先进(SOTA)模型的比较

    表4 与其他方法的性能对比(%)    

     

    表4所示,我们进行了与SOTA模型的性能比较的全面评估。模型涵盖三类:

    ① 基础模型:ResNet、WDCNN、RCL和MSResNet。

    ② 首层可解释范式模型:EWSNet、WaveletKernelNet、SincNet、FCC、MCNWFK、TFN、ClassBD、DFAWNet[81]、GTFENet[82]、wave_convNext[83]和RaVEL[84]。

    ③ 小样本下表现良好的模型:MCSwinT[85]、Convformer_NSE[86]、CLFormer[87]、Liconvformer[88]。    

    ResNet、WDCNN、RCL和MSResNet是数据驱动的。它们受样本大小的显著影响。在小样本下,不可能获得足够泛化的模型。因此,它们的平均值都低于82%。尽管ResNet表现出相对较好的性能,但它有大约380万个参数和较高的计算复杂度。另一方面,RCL对输入原始数据进行随机卷积数据增强。与WDCNN相比,虽然增加了计算复杂度,但提高了7.62%。总体而言,这些模型受到有限样本的限制,未能提取判别性特征,导致准确率较差。    

    此外,首层可解释模型基于信号处理赋能的算法,如小波、STFT、BD和时频变换,以减少对数据质量和数量的依赖。与数据驱动模型相比,这些算法取得了改进。其中,EWSNet和RaVEL表现最佳,分别为86.51%和86.79%。然而,与PIFCapsule相比,它们的性能落后7.05%和6.77%。EWSNet主要采用小波但忽视了其他基于物理的信息,并且采用多尺度结构,提取多尺度小波特征。RaVEL可视为将小波注入RCL。它在数据预处理阶段提高了SNR,增强了数据集的质量。因此,RaVEL对数据的要求较低,计算复杂度极低,为0.70 MB。与上述两个模型相比,

    PIFCapsule包含三种基于物理的信息,并使用AFR成功融合它们,达到93.56%。    

    测试了几种用于小样本诊断的轻量级Transformer模型。尽管MCSwinT达到了83.46%的最高准确率,但其参数数量和计算复杂度分别是PIFCapsule的4倍和17倍。CLFormer足够轻量,但平均准确率低于50%。PIFCapsule相比MCSwinT提高了10.1%。这表明即使面对最流行的Transformer架构,所提出的基于CapsNet的模型也具有竞争力。    

    4.2.5 与其他最先进多源融合方法的比较

    表5 在十个样本下与先进多源特征融合算法的对比    

     

    表5所示,我们对几种SOTA多源特征融合模型进行了比较分析,以评估AFR的优势。CSST_Net和MSIFT与PIFCapsule有相似之处,因为它们都采用自注意力结构。CSST_Net主要通过连接实现数据融合,而MSIFT利用交叉注意力机制实现特征融合。    

    表6 在数据集G下PIFCapsule中每个模块的功能对比    

     


    在本文中,CSST_Net通过连接小波、STFT和BD视图的三个特征进行数据融合。MSIFT输入小波和BD视图。如表6的消融实验所示,这两个视图的影响最显著。    

    结果表明,PIFCapsule在不同数据集上分别提高了0.98%和7.58%。尽管MSIFT具有竞争力,但其参数数量是PIFCapsule的2.46倍,计算复杂度是PIFCapsule的45.11倍。因此,其效率远低于PIFCapsule。此外,CSST_Net和MSIFT都使用相对复杂的Transformer架构,没有轻量化。

    4.2.6 消实验

       

    为评估不同模块的影响,如表6所示,STFT视图、小波视图和BD视图带来的改进分别为1.28%、43.03%和6.25%。显然,小波视图影响最显著,表明多尺度和高分辨率时频特征可以捕捉绝大多数故障特征,从而显著提高性能。STFT视图和BD视图起补充作用,BD视图特别提取周期性特征。    

    此外,我们研究了不同度量的影响。BSQ的主要目的是增强时频表示的能量集中度,而NTRA主要用于减少噪声并提取周期性脉冲特征。采用BSQ时,性能提高了1.94%。BSQ增强了提取的时频表示的能量集中度,并突出了特征频带。使用NTRA时,性能提高了3.68%。NTRA主要通过可微分软阈值模块学习周期性信号中的噪声,从而放大周期性脉冲并突出故障特征频率。总之,本研究中的两个主要损失函数相互补充,提高了原始信号的质量。

    5 PIFcapsule的可解释性分析

    5.1 多视图首层瓶颈的可解释性

    为验证首层瓶颈的可解释性,我们以货运轮对轴承为例,其故障频率主要低于1000 Hz。基于表6的消融实验,小波视图信息是最关键的因素。因此,将重点分析小波权重视图。首先,计算小波权重的幅频响应h(f):

       

    其中     表示每个通道的权重,     表示通道数。

     

    图11 小波权重的频率响应

     

    图12 8种状态FFT分析    

    图11所示,小波权重的幅频响应在351.56 Hz处出现最高峰值,而原始卷积的幅频响应近似一条直线,无法反映信号的物理特性。此外,我们使用快速傅里叶变换(FFT)对八种不同状态的原始信号进行了分析,如图12所示。对应于小波权重,各种状态在300 Hz附近的分布相对密集。这一观察表明,小波权重可以在300 Hz附近区分不同状态,并有效提取故障频率特征。    

     

    图13 BD视图输出结果的包络频谱分析    

    此外,如图13所示,对BD层滤波信号的包络谱幅度分析证实了这一发现。显然,幅度在300 Hz附近达到最大值,从而进一步证实小波核的平均频率响应主要集中在300 Hz附近。同时,观察到一个有趣的现象:如图(a)、(e)、(d)、(h)所示,最大或第二幅度点位于相同的频率304.96 Hz(289.06 Hz)和664.06 Hz,这可能导致混淆和误分类。因此,需要额外的视图信息——如小波或STFT视图——来进一步区分故障类型。    

     

    图14 不同通道中窗口权重的频率分布    

    最后,如图14所示,我们分析了不同通道上窗权重的频率分布。训练后,初始索引通道的频率主要分布在300 Hz附近,而这些通道在训练前分布在0 Hz附近,其他通道未观察到显著变化。这揭示了三个关键现象:

    ① 窗权重证实频率分布集中在300 Hz附近,可以提取判别性特征。

    ② 窗权重初始化近似最优窗权重,无需大量调整。

    ③ 在初始索引通道中,滤波器以300 Hz为中心,不同通道在不同频率提取原始信号的主要特征,主要位于低频范围内。

    5.2 路由系数Z的可解释性

    通过使用数据集F的案例研究,在PIFCapsule中,低级胶囊由从三个信息视图(小波、STFT和BD)提取的特征构建而成。使用了16组低级胶囊,这些胶囊随后通过AFR机制聚合成高级胶囊,产生8组。为了量化低级胶囊对高级胶囊的贡献,进行了以下分析。 如图4所示,低级胶囊的表示可以表示为 m

           

    对于任何给定的低级胶囊        ,原始贡献可以表示为:

           

    而实际贡献可以表示为:

           

    其中          表示与胶囊         相关的路由系数或注意力权重。PIFCapsule避免池化操作。在促进低级和高级胶囊之间的路由通信方面至关重要,从而实现三个基于物理的低级胶囊的融合。提供了对PIFCapsule优势的洞察。 因此,最终热力图计算为:        

    其中         是一个小常数。


     

    图15 热力图:从低级胶囊映射到高级胶囊    

    图15所示,热图直观地表示了低级胶囊对高级胶囊的贡献程度。尽管不同低级胶囊的贡献各不相同,但它们都是正的,表明整合了三个基于物理信息源的低级胶囊对故障定位做出了积极贡献。

    6 结论

    轨道交通车辆诊断面临高转速、重载荷和有限样本的挑战。FLIP——一个植根于定制化制造和具身智能的智能框架——系统地将人类先验知识嵌入智能诊断设备。基于此范式,提出了一种多视图信号处理先验嵌入的PIFCapsule。首先,提出了一种基于物理启发的多视图方法,结合小波变换、短时傅里叶变换和盲解卷积,提取不同的物理特性,并最大限度地利用跨时间、频率和时频的互补信息。其次,提出了一种新的非迭代且计算高效的注意力融合路由机制。AFR不仅能应对减少的胶囊数量,提高知识的内在转换能力,还能增强泛化能力。关键的是,在路由过程中,AFR充分整合了互补信息的变化,拓宽了信息空间维度,并成功捕捉了隐式的视图间关联。最后,提出了噪声阈值幅度比度量作为基于物理的正则化项,通过削弱学习到的噪声来放大与冲击相关的故障特征。

         

    然而,有几个方向值得进一步研究。首先,与一些单源基于物理的信息模型相比,PIFCapsule的参数数量和复杂度仍然相对较高。需要进一步研究在保持性能的同时实现更轻量级的模型。其次,PIFCapsule的特征具有更多的物理可解释性和泛化能力,这在不同机器间的迁移学习中具有潜力。第三,首层范式可以涵盖更多类型的信号处理,并开发更多信号处理赋能的卷积或基于Transformer的特征融合模块。最后,通过大型语言模型(LLMs)[92]增强现有算法的性能是一个有趣的研究方向。通过利用LLM代理或助手,研究人员旨在进一步优化现有算法的

    轻量级特性和性能。未来,可以探索将具身智能的感知-行动循环与LLMs的推理能力深度融合,从而构建能够主动交互、持续学习和自主决策的智能代理。这将实现对复杂工业系统中故障的动态理解、因果解释和自适应诊断。

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

    来源:故障诊断与python学习
    MechanicalSystem振动断裂非线性旋转机械轨道交通UG通信海洋UMLMS声学裂纹理论电机传动机器人数字孪生
    著作权归作者所有,欢迎分享,未经许可,不得转载
    首次发布时间:2026-03-25
    最近编辑:5月前
    故障诊断与python学习
    硕士 签名征集中
    获赞 87粉丝 152文章 327课程 0
    点赞
    收藏
    作者推荐

    “数字孪生/故障诊断”专题||数字孪生+视觉语言大模型:用于多场景工业故障诊断(下)

    工业设备的故障诊断,能不能像 ChatGPT 一样:看懂时频图、说清故障原因、还能给维修建议?这篇一区顶刊论文给出的答案是:可以,数字孪生 × 多模态大模型:工业故障诊断进入“可对话时代”🤖。本期推荐的这篇是西安交通大学许权宁的文章,本文系统综述了工业故障诊断从“数据驱动深度模型”走向“大模型时代”的完整技术脉络:首先回顾了传统规则/模型/数据三类范式的演进瓶颈,继而梳理了大型语言模型、视觉模型与多模态视觉-语言模型(LVLMs:Large Vision-Language Models)在故障预测与健康管理(PHM:Prognostics and Health Management) 领域的最新进展,指出通用大模型因缺乏工业知识、难以应对跨工况分布漂移与高质量样本稀缺等痛点;在此基础上,重点剖析了数字孪生驱动的深度合成数据、参数高效微调(LoRA:Low-Rank Adaptation)、跨模态对齐与对比学习等关键技术如何协同解决样本不足、域泛化与可解释性三大核心问题,最终提出面向多场景、可对话、可迁移的下一代智能故障诊断大模型研究框架与未来挑战。 由于文章篇幅过长,小编将分两次为大家翻译介绍这篇论文。第一篇推文系统梳理了大型视觉-语言模型在工业故障诊断中的瓶颈与突破:围绕“高质量跨工况数据难获取、模型通用知识难迁移”两大痛点,作者先回顾数据驱动、模型驱动与知识驱动三类智能诊断范式的演进,指出传统深度方法在跨域泛化上的先天缺陷;随后提出“深度数字孪生”思路,用六自由度轴承动力学+CycleGAN 虚拟-实测闭环映射,低成本生成高保真、多健康状态振动样本,为后续大型视觉-语言模型训练提供足量且物理一致的时频图谱数据,奠定多场景诊断基础。 本篇推文聚焦实验验证与落地价值。在CWRU、HUST、PU、XJTU四数据集上,FDGLM单工况准确率99.5%,跨工况仍达91.5%,跨数据集95.9%,均显著优于CNN、ViT、VisualGLM等基线;消融试验证实两阶段微调与双损失缺一不可,LoRA秩8、每类3000样本即达最优,数字孪生数据把真实信号相似度从0.5提升到0.92,为工业现场少样本、强噪声环境下的可迁移、可解释智能维护提供了即用型解决方案。 希望对大家的学习有所帮助,文章质量很高同时希望大家可以多多引用。如果有故障诊断相关方向研究人员希望宣传自己研究成果,欢迎大家在公众 号后台与小编联系投稿,大家一起交流学习。 论文链接:通过点击最左下角的阅读原文进行在线阅读及下载。 论文基本信息论文题目: Deep digital twin-powered large vision-language model for multi-scenario industrial fault diagnosis论文期刊:Advanced Engineering InformaticsDoi:https://doi.org/10.1016/j.aei.2025.103997作者: Quanning Xu(a), Guangrui Wen(b), Zihao Lei(c),Shulong Gu(d),Yu Su(e),Zhifen Zhang(f),Xuefeng Chen(g) 论文时间: 2025年 机构: a National Key Lab of Aerospace Power System and Plas ma Technology, Xi’an Jiaotong University, Xi’an 710049, Chinab State Key Laboratory for Manufacturing System Engineering, Xi’an Jiaotong University, Xi’an 710049, Chinac School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China作者简介:许权宁,西安交通大学机械工程学院博士研究生,国家航空航天动力系统与等离子体技术重点实验室、机械制造系统工程国家重点实验室核心成员。长期聚焦工业人工智能、数字孪生与多模态大模型融合研究,近三年在《Advanced Engineering Informatics》等顶刊以一作/通讯身份发表多篇论文,主导开发了面向多场景工业故障诊断的开源框架 FDGLM。(来源:ResearchGate) 摘要数据驱动的深度学习技术已被广泛采纳于工业健康监测与维护。然而,传统诊断方法常受场景特异性限制且智能能力不足。新近发展的大规模视觉-语言模型为智能故障诊断提供了先进的人机交互范式。本研究通过提出故障诊断通用语言模型(FDGLM:Fault Diagnosis General Language Model)框架,将LVLMs拓展至工业故障诊断,实现跨多种运行工况与数据集的智能化诊断。首先,依托先进数字孪生技术生成高质量振动样本以支撑模型训练;这些样本经短时傅里叶变换(STFT:Short-Time Fourier Transform)处理,生成增强的时频谱图,满足故障诊断通用语言模型(FDGLM:Fault Diagnosis General Language Model)训练对大规模数据的需求。随后,采用低秩自适应策略对预训练视觉模型进行微调,联合交叉熵损失与圆损失以同步提升对工业故障时频特征的判别能力与鲁棒性。接着,利用领域专用文本指令对语言模型进行微调,实现视觉与文本模态的深度对齐,使因果故障分析与维护决策成为可能。最后,在四个独立数据集上开展跨工况与跨数据集实验,结果显示:同工况诊断精度达0.995,变条件为0.890,跨数据集为0.947,同时呈现专业级对话诊断能力。实验结果证实,FDGLM以极小微调代价即可提供适用于工业场景的可靠诊断,性能优于现有框架,为下一代工业装备健康管理提供解决方案。 关键词:大语言模型;大规模视觉-语言模型;故障诊断;多模态大语言模型;数字孪生 目录1 引言2 相关研究 2.1 智能故障诊断 2.2 大规模模型3 深度数字孪生模型 3.1 滚动轴承的虚拟表示 3.2 虚拟-物理数据流融合4 所提出的大规模视觉-语言模型 4.1 模型架构 4.2 预训练视觉模型的微调 4.3 查询变换器 4.4 预训练语言模型的微调 4.5 实现流程5 实验验证 5.1 数据集配置与预处理 5.2 实验设计与实现细节 5.3 数字孪生模型验证 5.4 单工况验证 5.5 跨工况验证 5.6 跨数据集验证 5.7 超参数分析 5.8 消融实验6 结论注:小编能力有限,如有翻译不恰之处,请多多指正~ 若想进一步拜读完整版,请下载原论文进行细读。 5 实验验证为全面评估所提出的数字孪生驱动 LVLM 在故障诊断中的性能,本文在四种轴承数据集上开展大量实验,覆盖三大场景:单工况、跨工况与跨数据集泛化。单工况场景仅含一台机器且参数固定,主要考核模型对标准故障模式的识别能力;跨工况场景让同一台机器在不同负载或转速下运行,用以检验模型对工况变化的鲁棒性;跨数据集场景则引入不同机器或实验平台,测试模型对异构数据源的适应性。该评估方案系统验证了所提框架在多样应用场景下的精度与鲁棒性。5.1 数据集配置与预处理实验共使用三个公开数据集与一个自采数据集,详情如下:(1) CWRU 轴承数据集:该广泛使用的基准数据集包含健康、内圈故障、外圈故障及滚动体故障四种状态;每种故障按缺陷尺寸再分轻、中、重三级,共在四种负载-转速组合下采集,采样频率 12 kHz,轴承型号 SKF 6205。(2) HUST 轴承数据集:涵盖健康状态及滚动体、内圈、外圈的轻微与严重故障,涉及多转速工况;采用三轴加速度计记录振动信号,轴承型号 ER-16 K,采样频率 25.6 kHz。(3) PU 数据集:包含内圈、外圈及复合故障的多级严重程度;故障通过人工损伤与加速寿命试验引入,采样频率 64 kHz,轴承型号 SKF 6203。(4) XJTU 数据集(研究团队自采):利用轴承故障模拟平台采集,含健康、内圈故障、外圈故障、滚动体故障四种状态,共十二种运行工况,采样频率 12 kHz,轴承型号 SKF 6207。基于上述数据集,本文构建了对应的多模态图像-文本数据集,用于模型训练与评估。具体而言,以 STFT 将一维振动信号转为二维时频表征:截取 1024 点信号段并变换为时频图像,作为视觉输入。针对数据稀缺问题,利用所提出的深度数字孪生模型预先大规模生成训练样本。每张时频图像配对一条硬提示,如 “What is the fault with the bearing?”,引导语言模型;对应的文本标签依据专家知识预设模板自动生成,内容涵盖诊断结果、健康状态描述及维护建议,亦可借助 ChatGPT 等大语言模型辅助生成。5.2 实验设计与实现细节鉴于真实工业场景数据采集困难且大模型天然面临数据稀缺,本研究采用深度数字孪生方法,合成覆盖多种健康状态的大规模振动信号,用于扩充数据集并支撑所提模型的充分训练。为验证框架有效性,本文与传统机器学习算法、深度学习算法及先进计算机视觉分类模型开展对比实验,具体基线方法包括多层感知机(MLP:Multilayer Perceptron)、支持向量机(SVM:Support Vector Machine)、卷积神经网络(CNN:Convolutional Neural Network),以及 VGG16、ResNet18、DenseNet、ViT、Swin Transformer、ConvNeXt、ChatGLM2-6B-chat [10] 和 VisualGLM [32] 等深度模型。为保证评估严谨全面,采用四项广泛认可的指标:召回率(REC:Recall)、F1 分数(F1:F1-score)、准确率(ACC:Accuracy)与精确率(PRE:Precision) [50],每项指标均在多次独立实验上取平均,以降低随机性并增强统计鲁棒性。所提出的 FDGLM 框架以预训练 ViT为视觉编码器,ChatGLM-6B 为文本编码器。为提升微调性能,LoRA 矩阵被均匀且分层地注入两个编码器的选定层 [7,32,51]:具体而言,LoRA 矩阵嵌入 ViT 的第 0、13、26、38 层以及 ChatGLM-6B 的第 0、13、27 层,并设为可训练参数。训练过程分两阶段,分别进行 10 000 次与 5 000 次迭代,学习率 0.001,批次大小 8。所有训练均在显存 86 GB 的 NVIDIA A40 GPU 上完成。值得注意的是,每种健康状态类别通过第 3 章所述数字孪生框架扩充至 3 000 个样本,确保所有故障类别数据均衡。后续消融实验进一步探究样本量与诊断性能之间的关系。5.3 数字孪生模型验证为评估所提出数字孪生模型的有效性,本节对比合成孪生信号与 XJTU 数据集真实测量的相似度与保真度。虚拟-物理映射网络使用 Adam 优化器训练,批次大小为 10,共 500 个 epoch,每类包含 150 个训练样本;初始学习率设为 0.01,每 50 个 epoch 衰减一次。所有实验在 Windows 11 系统、PyTorch 2.1 环境下完成,计算平台为 Intel Core i7-12650H CPU、GeForce RTX 4060 GPU 及 16 GB RAM;更多实现细节见文献 [41]。轴承动态参数汇总于表 1,图 10 给出外圈故障工况下的信号对比。通过定义几何与运行参数,虚拟实体生成各健康状态的仿真振动响应;这些响应经预训练的虚拟-物理映射网络映射,得到对应孪生信号。以外圈故障为例,图 10 在时域、频域与包络域分别对比仿真信号、孪生信号与真实信号;为公平比较,所有信号幅值均归一化至 [-1, 1] 区间。表1 采用SKF - 6207中的轴承参数 图10 对比分析了XJTU数据集的外圈故障响应:( a )、( d )和( g )说明了测量数据的时间轮廓、频谱特征和解调频谱;( b )、( e )和( h )与模拟数据相对应;而( c ),( f )和( i )描述了虚拟孪生数据。在时域中,三种信号均呈现重复性冲击模式且变化趋势一致;仿真信号因虚拟实体理想化构建而波形清晰、冲击尖锐,不含测量噪声、制造公差或安装误差,真实信号则受环境噪声畸化。包络谱中,所有信号均清晰显示故障特征频率及其谐波;值得注意的是,仿真信号对内圈与滚动体缺陷相关频率成分予以强化,而真实信号中这些成分因多源振动与噪声被削弱,给直接模型更新与校准带来困难。包络谱幅值方面,仿真与真实信号在特征频率位置一致但幅度不同,孪生信号有效补偿该失配:在虚拟-物理映射过程中,引入环境与故障相关信息,同时保留测试观测到的高频共振响应,使仿真信号转化为孪生信号后故障频率及谐波得到凸显,表明该方法在吸收环境影响的同时保持故障信息。该性能得益于映射网络的多层级架构:Self-ONN 编码器实现层内高效特征提取,注意力增强跳跃连接强化层间表征,多深度编码器保证多尺度学习,而采用复合损失的对抗训练机制促进跨域信息融合与分布对齐 [41]。为进一步评估信号相似度,对归一化信号进行分布拟合分析。图 11 给出滚动体、内圈、外圈三类故障在单位时间内的振动响应高斯拟合结果,表明映射网络有效对齐跨域细节与空间特征,融合环境及振动信息,将仿真信号转化为近似真实测量的孪生信号,显著缩小分布差异。此外,采用皮尔逊相关系数量化信号相似度(表 2):仿真-真实信号相似度约 0.5,孪生-真实信号相似度超过 0.92,大幅提升其用于诊断模型训练的价值。 图11 数据分布拟合于XJTU数据集。( a )滚动体故障,( b )内圈故障,( c )外圈故障表2 不同状态下信号相似度的比较 5.4 单工况验证 表 3 与表 4 分别给出 XJTU 与 HUST 轴承数据集在单一工况下的诊断结果,展示各方法性能。经 STFT 处理后,振动信号的时频表示显著提升了轴承故障状态的可表征性 [52]。在所有轴承状态中,所提 FDGLM 始终优于其他方法,总体诊断准确率超过 99.5%。传统机器学习算法(包括 SVM 与 MLP)表现较差;在训练数据充足的前提下,基于 Transformer 的网络较传统卷积神经网络优势明显,因此 DCNN、ResNet18、DenseNet121 与 VGG16 性能均低于 ViT、Swin-Transformer 及所提 FDGLM。值得注意的是,ResNet18 与 DenseNet121 利用多通道特征复用提升分类性能,在 HUST 数据集上分别取得 95.678% 与 97.282% 的准确率。采用分层移位窗口注意力的 Swin-Transformer 位列第二,在 XJTU 与 HUST 数据集上分别达到 99.300% 与 99.100% 的准确率。ChatGLM2-6B-chat 作为代表性大语言模型,依赖 24 个手工提取的时域与频域特征输入 [10],在两个基准数据集上平均诊断准确率为 98.5% 与 98.2%;然而,其单一文本模态的依赖性限制了诊断能力的进一步提升。相比之下,VisualGLM 采用两阶段协同微调策略进行领域适配,输入为时频图像,在同一基准数据集上分类准确率达到 99.1% 与 98.9%。基于 Transformer 架构的 FDGLM 采用两阶段微调策略优化时频图像特征表示,并利用 Query Transformer 对齐视觉与文本模态。图 12 给出表现最优模型的混淆矩阵,DenseNet121、VGG16、ViT、Swin Transformer、ChatGLM2-6B-chat、VisualGLM 与 FDGLM 均实现接近完美的故障分类。图 12 中,坐标轴标签 1–7 分别对应以下轴承状态:健康状态、轻微与严重滚动体故障、轻微与严重内圈故障、轻微与严重外圈故障。健康状态易于区分,而同类型不同严重程度的故障分类难度较大,导致基线方法出现多数误分类。 表3 在XJTU数据集上的对比实验结果 表4 在HUST数据集上的对比实验结果 图12 对 HUST 数据集上各种方法的混淆矩阵进行比较。(a) 支持向量机 (SVM),(b) 多层感知器 (MLP),© 深度卷积神经网络 (DCNN),(d) ResNet18,(e) DenseNet121,(f) VGG16,(g) ViT,(h) Swin-Transformer,(i) ConvNeXt,(j) ChatGLM2-6B-chat,(k) VisualGLM,(l) FDGLM 由表 3 可见,除表现最差的 SVM 外,FDGLM 在微调阶段的可训练参数最少(0.86 M)。尽管总参数量达 72 亿,FDGLM 借助 LoRA 技术,仅在视觉编码器第 0、13、26、38 层与语言编码器第 0、13、27 层插入轻量可训练矩阵,其余预训练参数全部冻结,既保留基座模型表达能力,又实现工业故障诊断的高效领域适配,凸显以极低微调开销完成工业数据处理的巨大潜力。训练平台上,我们对各方法进行单样本推理延迟量化分析:如表 3 所示,传统机器学习与深度学习模型推理时间远短于大规模模型,差异主要源于参数量与结构复杂度。值得注意的是,ChatGLM2-6B-chat 的推理延迟高出一个数量级,根源在于其生成机制与架构设计:相比仅需一次前向的传统判别模型,LLM 通常采用自回归策略生成文本序列,输出长度 L 需 L 次前向,且每次均包含计算开销远高于卷积的自注意力机制;LVLM 额外引入视觉编码器与跨模态融合模块,推理时间最长。未来可通过量化(将权重与激活由高精度转低精度以减少内存、加速计算)、系统级推理优化与模型编译、架构级剪枝与知识蒸馏等手段提升大规模模型推理效率。我们进一步在搭载 i7-12650H 处理器、GeForce RTX 4060 GPU 与 16 GB 内存的笔记本上进行推理实验,延迟虽适度增加,模型仍稳定运行,验证了边缘部署的可行性。此外,FDGLM 通过人机交互实现智能故障诊断:其友好界面即使对非专业用户也能输出精准健康评估。如图 13 所示,界面分为左右两栏——左侧展示时频谱图与输入区,右侧呈现诊断结果与维护建议。与传统方法仅输出标签不同,FDGLM 依托知识增强语料,可生成详细健康评估、推断潜在故障根因并提供维护指导;同时支持文本与视觉多模态输入,弥合图像诊断与自然语言处理的鸿沟。图 14 与图 15 给出更多健康状态的定性诊断示例,可见模型在不同文本提示下均保持强鲁棒性,持续实现准确故障识别、因果分析与建议生成。作为 PHM(Prognostics and Health Management)领域的重大进展,该视觉-语言框架以智能、交互的解决方案超越传统方法,将精准诊断与易用部署相结合,显著提升工业设备的智能维护水平。 图13 轴承健康状态评估人机交互界面的演示 图14 轴承健康状态评估人机交互界面的演示在HUST数据集上使用FDGLM进行故障诊断的定性案例研究。( a )轻微滚动体故障,( b )严重滚动体故障,( c )轻微内圈故障 图15 不同文本提示的故障诊断定性案例研究。( a )严重的内圈故障,( b )轻微的外圈故障,( c )严重的外圈故障 为评估通用 LVLM 在工业时频图像识别中的有效性,本研究对比三款代表性模型:Qwen 3、DeepSeek-V3.1 与 GPT-5。图 16 以 HUST 数据集中“轻微内圈故障”时频图为输入,统一查询指令为 “What is the fault with the bearing?”。Qwen 3 仅检测到谱图中与机械缺陷相关的离散异常频率成分,未能指明具体故障类别,亦未提供失效机理解读或维护建议。DeepSeek-V3.1 起初未识别输入为轴承振动信号时频图;在补充提示 “This is a time–frequency diagram of the bearing vibration signal” 后,仍无法判定故障类别,仅能分析潜在原因并给出维护建议。GPT-5 表现最佳,直接识别出内圈故障并构建连贯推理链,但仍缺乏对根本失效机理的解释、决策支持洞察及故障严重度量化。实验揭示通用 LVLM 在工业诊断中的两大关键局限:(i) 专用时频特征辨识能力不足,导致诊断精度下降;(ii) 训练语料工业知识匮乏,制约其生成符合 PHM 需求的失效分析与维护策略。这些缺陷凸显了本文提出 FDGLM 框架在工业故障诊断中的必要性与技术优势。 图16 在HUST数据集上使用其他LVLM进行故障诊断的定性案例研究。( a ) Qwen 3,( b ) Deepseek-V3 . 1,( c ) GPT - 5 5.5 跨工况验证为评估所提 FDGLM 模型在未见运行工况下的诊断泛化能力,我们利用 CWRU 与 HUST 轴承数据集设计跨工况迁移诊断实验。如表 5 所示,每个数据集包含四种不同的转速-负载组合工况;据此,为每个数据集构建四项跨工况诊断任务(表 6),用以检验模型在不同运行工况下的诊断迁移性能。 表5 CWRU和HUST数据集的运行工况 表6 跨工况诊断任务配置与诊断准确性 实验结果表明,尽管不同运行条件导致分布偏移,FDGLM 在所有跨工况任务中仍保持高诊断准确率:在 CWRU 数据集上,面对转速与负载变化,准确率均超过 93.0%;在转速变化显著的 HUST 数据集上,任务 3 准确率最低(89.0%),四项任务平均准确率达 91.5%。这些结果证实 FDGLM 在单数据集跨工况诊断中具备强适应性与鲁棒性;尽管性能较单工况略有下降,但在未改动模型架构或损失函数的前提下取得该表现,凸显了框架固有的泛化能力。为深入评估跨工况诊断性能,本文聚焦 HUST 数据集任务 4 开展对比实验,基线方法与第 5.4 节保持一致,所有结果取多次独立运行平均,并给出最优运行的混淆矩阵。如表 7 所示,相较于单工况基准,各方法在跨工况场景下均出现性能下降;FDGLM 仍以 93.3% 的分类准确率领先。Swin-Transformer 与 ConvNeXt 降幅较小(约 1–1.5%),性能逼近 FDGLM;传统 CNN 方法鲁棒性不足,准确率普遍低于 90%。ChatGLM2-6B-chat 仅依赖交叉熵损失进行故障诊断,且受限于单一文本模态,跨域诊断效能受限;VisualGLM 虽采用结合交叉熵与对比损失的两阶段协同优化策略,但在工况差异显著的跨域工况下表现仍不理想,表明传统对比损失仍有较大改进空间。表7 在HUST数据集上针对任务4的对比实验结果 图 17 的坐标轴定义与图 12 保持一致。图 17 显示,与第 5.4 节中同一故障类别内的严重程度误分不同,跨工况验证导致更复杂且分散的错误模式:工况偏移使多数基线模型不仅混淆同一故障的不同严重程度,还错分不同故障类型,充分暴露这些模型在未见工况下面对分布偏移时的泛化局限。FDGLM 在跨工况任务中的鲁棒性部分源于视觉编码器微调阶段采用的 Circle 对比损失——该损失通过增大正样本对相似度、减小负样本对相似度,优化特征空间,促使模型学到在工况变化下仍保持稳定的判别性故障特征,从而提升分布偏移下的泛化能力。总体而言,结果表明 FDGLM 能够捕捉轴承信号在不同工况下的共性故障模式,并将其泛化至未见场景进行诊断,显著增强了模型在多变工况内的诊断泛化能力。 图17 在HUST数据集上对任务4的各种方法的混淆矩阵进行比较。( a ) SVM,( b ) MLP,( c ) Dcnn,( D ) Resnet18,( E ) Densenet121,( f ) Vgg16,( G ) VIT,( H ) Swin-Transform,( I ) Convnext,( J ) Chatglm2-6B-Chat,( K ) Visualglm,( L ) Fdglm 5.6 跨数据集验证为评估所提 FDGLM 模型的跨域诊断能力,本文利用 XJTU、CWRU、HUST 与 PU 四个独立轴承数据集开展全面跨数据集迁移诊断实验。实验设计确保各数据集健康状态类别严格一致,仅聚焦三类代表性状态:健康、内圈故障与外圈故障。各数据集在轴承型号(如 SKF6205、ER-16 K、SKF6207)、运行工况(转速与负载变化)及环境噪声等方面存在显著异质性,给故障诊断的领域适配带来巨大挑战。表 8 汇总了实验设置与诊断性能。 表8 跨数据集诊断任务配置与诊断准确率 FDGLM 在所有跨数据集任务中的平均分类准确率达 95.9%,彰显其出色的领域泛化能力;值得注意的是,跨数据集任务的平均诊断精度高于跨工况场景,这可归因于跨数据集设置中健康类别更少、诊断目标更一致。任务 1 获得最高诊断准确率,表明模型对 XJTU 数据集的适应性最佳。HUST 数据集采用 ER-16 K 轴承,与其他三个数据集使用的 SKF 系列差异显著,导致任务 3 分类准确率最低(94.7%),说明诊断模型对轴承型号变化存在一定敏感性——该敏感性源于轴承型号与负载-转速条件交互带来的显著分布偏移。结果证实,即使在完全不同数据集上训练与测试,FDGLM 仍保持强劲诊断性能,体现出其提取域不变故障特征并将知识有效迁移至未见领域的能力,显著提升了模型的跨域泛化水平。 5.7 超参数分析 5.7.1 LoRA矩阵的秩本节研究 LoRA 矩阵秩对所提模型微调性能的影响。LoRA 秩直接决定微调期间引入的可训练参数数量,并显著影响模型对轴承故障诊断的适应能力。尽管更高秩在理论上提供更大参数容量,实际微调性能仍受计算资源、样本规模等实验条件约束。为隔离这些因素,我们在受控实验设置下分析 LoRA 秩与诊断准确率的关系;不同层中的所有 LoRA 矩阵共享同一秩配置。如表 9 所示,秩按 2 的幂次系统设置为 4、8、16、32 与 64 [32],并在此配置下评估预训练 ViT 与 ChatGLM-6B 的微调性能。提高秩可扩大可训练参数规模、增强模型表征能力,但也增加数据需求与计算开销。实验结果表明,在计算与数据资源受限条件下,诊断准确率并非随 LoRA 秩单调上升,而是呈现非单调趋势:当秩设为 8 时,可训练参数总量为 0.86 M,FDGLM 取得最佳性能;秩超过 16 后,准确率较峰值略有下降,原因在于额外参数在有限数据与硬件条件下引入更高的过拟合风险与优化难度。因此,秩 = 8 在当前配置下达到最优,并被用于本研究所有基于 LoRA 的微调。 表9 不同LoRA矩阵秩次下的诊断性能 5.7.2 训练数据量如前所述,真实轴承振动数据不足以支撑大规模视觉-语言模型训练:例如,CWRU 与 HUST 数据集采样频率分别为 12 kHz 与 25.6 kHz,但每种健康状态仅提供约 10 s 有效数据;采用 1024 点无重叠分段,每类仅生成约 117 与 250 个样本。即便借助重采样等数据增广,样本总量仍远低于有效训练需求,且增广可能加剧过拟合风险。工业现场故障呈偶发、昂贵特性,进一步导致样本不平衡等复杂挑战。为缓解数据稀缺,本文采用深度数字孪生框架合成覆盖不同健康状态的多样化振动信号;本节考察合成数据集规模与诊断性能之间的关系。为评估训练数据量的影响并确定达到高准确率所需的最小样本规模,本文开展控制实验:在固定全部超参数与架构的前提下,测试每类 1 000、2 000、3 000、4 000 及 5 000 张时频图像五种样本量,这些图像均由合成振动信号经 STFT 生成;模型在每种配置下训练,并在独立测试集上评估。如图 18 所示,实验结果显示样本量与诊断准确率呈强相关:当每类样本为 1 000 与 2 000 时,模型准确率分别为 0.911 与 0.933,未达理论潜力,表明有限数据无法覆盖故障变异,削弱特征鲁棒性;当每类样本增至 3 000 时,准确率大幅提升至 0.995;继续增至 4 000 与 5 000 仅带来边际提升,无统计学显著改善。因此,在当前数字孪生配置、模型架构与计算约束下,每类 3 000 样本足以使 FDGLM 性能接近最优,结果证实深度数字孪生在扩充训练数据、缓解真实数据局限方面的有效性。 图18 不同样本量下的诊断性能 5.8 消融实验为严格评估 FDGLM 框架各核心组件的有效性与适应性,本节开展全面消融实验。整体 FDGLM 微调流程分两大阶段:视觉模型微调与语言模型微调。具体而言,在预训练视觉主干适配过程中引入两类损失函数进行联合优化:(1)交叉熵损失,用于增强故障类别可判别性;(2)Circle 对比损失,旨在提升所学特征表示的鲁棒性与泛化性。据此,本研究系统探究两微调阶段的单独及联合效应,以及上述两种视觉损失函数的各自贡献与协同作用。为保证公平与可比性,所有消融实验仅变动被测组件,其余配置完全一致;评估在 HUST 数据集单一工况下进行,结果汇总于表 10。若直接将预训练语言与视觉模型用于工业数据而不经 LoRA 微调,诊断性能急剧下降,表明预训练模型仅与通用视觉-文本输入对齐,难以从工业谱图中提取具有判别力的类别相关特征。仅微调语言模型时,诊断准确率显著下降且不稳定,说明视觉编码器若未适配,模型无法捕捉显著时频模式,诊断能力受限。单独使用交叉熵或 Circle 损失,诊断准确率分别为 0.977 与 0.973;二者在视觉模型微调中联合应用时,诊断性能进一步提升,凸显互补性:交叉熵损失强化决策边界学习,Circle 损失促进更具判别力且鲁棒的特征嵌入,协同提升编码器提取有效故障特征的能力。借助优化后的视觉编码器,再对语言模型进行微调,可将更丰富的视觉语义与文本表征融合,带来显著性能增益。综上,消融研究证实:视觉模型微调是 FDGLM 适配工业故障诊断的基础;交叉熵与 Circle 损失联合使用优于单独应用;所引入的两阶段微调方案通过联合优化视觉与语言组件,充分释放模型潜力,实现准确且鲁棒的故障诊断。 表10 不同成分的消融研究结果 6 结论本文提出一种新颖的数字孪生驱动大规模视觉-语言模型 FDGLM,为工业装备故障诊断提供智能、交互式解决方案。针对高质量多模态训练数据获取难题,采用深度数字孪生技术合成物理一致的振动信号,构建大规模训练语料;并引入两阶段领域适配策略,弥合通用预训练模型与特定工业应用之间的鸿沟。具体而言,视觉编码器以 LoRA 微调,联合交叉熵与 Circle 对比损失,增强时频图像的模式识别能力;语言编码器则借助专家诊断语料微调,实现跨模态语义对齐、故障推理与维护决策支持。大量实验表明,FDGLM 在多种诊断场景中表现优异:单工况准确率 99.6%,变工况 91.5%,跨数据集迁移 95.9%。超参数调优确定了最佳 LoRA 秩与最少样本量;消融研究验证了两阶段适配策略及集成视觉损失设计的有效性。本研究为工业系统智能运维开辟了新范式。未来工作将聚焦于解决真实工业环境中突发故障带来的小样本与零样本挑战;此外,编辑:赵栓栓校核:李正平、陈凯歌、曹希铭、赵学功、白亮、任超、海洋、陈宇航、Kira、Tina、王金、赵诚、肖鑫鑫、张优该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除来源:故障诊断与python学习

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