尧图建网站 尧图建网站 YAOTU WEB BUILD 免费咨询
ARTICLE DETAIL

资讯详情

深耕网站建设与建站编程的一线实战洞察。

【轴承故障诊断】用于轴承故障诊断的集中时频分析研究(Matlab代码实现)

【轴承故障诊断】用于轴承故障诊断的集中时频分析研究(Matlab代码实现) 欢迎来到本博客❤️❤️博主优势博客内容尽量做到思维缜密逻辑清晰为了方便读者。完整资源、论文复现、期刊合作、论文辅导及科研仿真定制事宜点击本文完整资源下载⛳️座右铭行百里者半于九十。⛳️赠与读者‍做科研涉及到一个深在的思想系统需要科研者逻辑缜密踏实认真但是不能只是努力很多时候借力比努力更重要然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览免得骤然跌入幽暗的迷宫找不到来时的路它不足为你揭示全部问题的答案但若能解答你胸中升起的一朵朵疑云也未尝不会酿成晚霞斑斓的别一番景致万一它给你带来了一场精神世界的苦雨那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。或许雨过云收神驰的天地更清朗.......第一部分——内容介绍用于轴承故障诊断的集中时频分析研究摘要滚动轴承是旋转机械设备的核心承载部件其运行状态直接决定整套装备的工作可靠性。轴承发生损伤时会产生周期性冲击振动该类振动信号具备显著的非平稳、非线性特征传统时域、频域分析手段难以同时刻画故障冲击发生的时刻与对应的频率信息。时频分析可以将一维振动信号映射至时间‑频率二维平面同步揭示信号随时间演变的频率组分是轴承故障特征提取的重要技术路径。传统时频方法普遍存在能量发散、时频聚集效果差、噪声抑制能力不足等缺陷微弱故障冲击容易被背景噪声掩盖造成早期故障漏检。集中时频分析以提升时频能量聚集度为核心目标通过对时频系数进行重分配、提取、聚合处理压缩信号能量的扩散范围强化故障瞬态冲击在时频平面上的特征表达能够在强噪声、变转速工况下凸显轴承故障特征。本文围绕集中时频分析在轴承故障诊断领域开展系统性研究梳理集中时频分析技术的发展脉络对比不同集中时频方法的技术特点分析该类方法在轴承故障识别中的优势与现存局限并探讨技术未来发展方向为滚动轴承故障特征提取与状态监测提供理论参考。关键词滚动轴承故障诊断非平稳振动时频分析时频聚集特征提取1 引言现代工业装备向着高速、重载、连续化运行方向发展滚动轴承广泛应用于风机、数控机床、轨道交通、大型发电机组等各类旋转机械当中轴承的局部损伤会逐步扩展最终引发整机停机事故带来经济损失甚至安全隐患。开展轴承故障诊断尽早识别早期微弱损伤对设备运维、故障预警具有重要工程价值。轴承滚道、滚动体出现点蚀、剥落等局部故障时滚动体反复撞击缺陷位置会生成周期性冲击脉冲。实际工业现场采集得到的轴承振动信号往往混杂齿轮啮合、电机电磁干扰以及环境噪声同时设备转速波动、负载变化进一步加剧信号的非平稳特性。单纯时域分析仅能观测冲击幅值变化无法定位故障对应的频率傅里叶变换类频域分析只能够获取全局频率构成丢失时间维度信息难以捕捉瞬态故障冲击的发生时刻面对变转速工况诊断效果大幅下降。时频分析技术把信号分解为时间与频率的联合分布实现时间、频率信息同步表征成为处理轴承非平稳振动信号的主流手段。短时傅里叶变换、小波变换等经典线性时频方法受海森堡测不准原理约束时间分辨率与频率分辨率无法同时达到最优故障冲击对应的能量会在时频平面发生弥散故障特征被噪声淹没。二次型时频分布存在交叉项干扰问题多分量信号下时频图谱可读性差不利于故障解读。集中时频分析是一类以提升时频能量聚集性能为核心的改进时频技术主要思路是基于原始时频结果结合信号瞬时特征对时频系数实施重排、筛选与聚合把原本发散分布的能量收拢到真实信号分量所在的时频脊线上弱化噪声带来的弥散伪影获得轮廓锐利、特征清晰的时频表达结果。近年来同步压缩变换、同步提取变换等一系列集中时频方法不断被提出逐步应用于轴承内圈、外圈、滚动体故障的特征提取在低信噪比、转速波动场景表现出优于传统时频手段的潜力。但集中时频分析在工程落地过程中依旧存在诸多现实问题面对强噪声、复合故障、转速剧烈变化的真实工况部分方法会出现特征丢失、时频畸变同时算法参数选择、计算效率、和后续智能诊断模型的适配性都有待进一步完善。基于此本文对面向轴承故障诊断的集中时频分析开展研究剖析轴承故障信号特点归纳集中时频分析技术体系对比各类方法的适用场景总结现存技术瓶颈展望后续研究方向。2 轴承故障振动信号与时频表征需求滚动轴承不同位置的局部故障会产生具有特定周期的冲击振动。当滚动体滚过损伤区域接触位置发生刚度突变激励起设备系统的高频共振故障冲击信息被调制在高频共振频段之上。早期故障冲击幅值微弱在原始振动波形中不易分辨随着故障不断发展冲击能量逐步提升。实际工况下轴承振动信号属于典型多分量非平稳信号包含转轴转动分量、齿轮啮合分量、各类谐波干扰再叠加现场环境噪声故障冲击特征被大量无关成分掩盖。在转速发生变化时故障冲击的时间间隔随之改变故障相关频率会随时间连续漂移传统依靠固定频谱的解调手段很难稳定提取故障特征。轴承故障诊断对时频分析结果提出明确需求。首先需要优异的时频聚集性能故障瞬态冲击应当在时频平面形成紧凑、边界清晰的能量团能量不能过度扩散其次要具备良好的抗噪声能力在低信噪比条件下依然可以保留故障对应的时频特征抑制噪声生成的伪时频分量再者面对转速波动能够跟踪时‑频脊线随时间的演变适配频率时变的故障信号同时应当尽可能减少交叉项、伪影干扰保证时频图谱物理可解释方便人工解读也便于作为图像输入给后续智能识别模型。传统时频方法很难同时满足上述全部要求。短时傅里叶变换受固定窗口约束时间、频率分辨率相互制约小波变换虽然实现多尺度分析但对于冲击信号依然存在能量扩散现象二次时频分布交叉项严重复杂信号下图谱杂乱。集中时频分析正是针对这些痛点发展而来通过后处理重分配机制在原有时频基底之上进一步压缩能量扩散专门适配轴承故障瞬态冲击的表征需求ResearchGa...。3 集中时频分析技术体系与技术特点集中时频分析本质上属于时频后处理增强技术并不完全抛弃传统时频变换框架大多以短时傅里叶变换、小波变换、S 变换作为基础依靠瞬时频率、群延迟等信号局部特征对分布在时频平面上的系数进行重新归置将分散的能量迁移到真实信号分量对应的时频位置以此实现能量集中得到高聚集度的时频表示。按照技术实现思路可以大致分为重分配类、同步压缩类、同步提取类三大技术路线。重分配变换是较早发展的集中时频思想该方法计算每一个时频点对应的局部瞬时参数把时频系数从计算得到的位置迁移到信号真实时频中心位置有效压缩能量扩散。重分配变换可以显著改善时频聚集效果但是该方法不具备信号逆变换能力完成时频增强之后无法重构回原始时域信号更多用于可视化特征观察难以和信号降噪、分量分离流程结合使用。同步压缩变换是目前故障诊断领域应用较多的集中时频方案该方法以小波变换或者短时傅里叶变换为基础依据估计得到的瞬时频率对时频系数进行沿频率方向的挤压聚合。不同于简单重分配同步压缩变换保留信号重构能力完成时频聚集之后可以逆变换得到各个信号分量既可以输出高聚集度时频图谱也能够实现多分量信号分解。在轴承故障场景中同步压缩能够收拢故障冲击对应的时频能量凸显随转速变化的故障时频脊线。当信号噪声过大、分量变化剧烈时同步压缩会出现脊线断裂部分故障能量会发生丢失影响微弱故障识别效果。同步提取变换在同步压缩的基础上进一步演化不再对全部时频系数做迁移挤压而是筛选保留与真实瞬时特征高度匹配的时频系数舍弃偏离信号脊线的弥散系数进一步强化时频聚集效果得到轮廓更加锐利的时频结果。该方法对瞬态冲击信号适配性较好适合提取轴承短时故障脉冲但是在噪声干扰下筛选阈值设置不当容易直接滤除部分微弱故障分量造成早期故障特征损失。除此之外学界还衍生出多阶同步压缩、多同步压缩、参数化重分配、瞬态提取变换等一系列改进集中时频方法针对调幅调频信号、强噪声、瞬态脉冲等不同信号特点做优化。部分改进方案同时在时间维度与频率维度开展系数重分配双向提升聚集性能更好适配轴承冲击这类时间高度局部化的信号特征。不同集中时频方法各有优劣不存在普适最优方案。重分配类聚集效果好但不可逆同步压缩支持信号重构适合多分量分解同步提取时频图视觉效果锐利但对参数与信噪比更加敏感。处理轴承故障信号时需要结合噪声水平、转速稳定程度、故障发展阶段选择适配的集中时频方案。4 集中时频分析在轴承故障诊断中的应用4.1 单故障特征提取轴承单故障场景指内圈、外圈或者滚动体单独发生损伤是实验室研究最常见的工况。集中时频分析接收轴承振动信号输出高聚集的二维时频图谱故障周期性冲击会在图谱中呈现周期性分布的高能时频斑块对应故障特征及其倍频成分。相较于传统时频方法集中时频处理之后故障冲击的时间位置更加精准频率脊线更加紧凑噪声生成的散乱背景被显著抑制。在恒定转速工况下故障相关频率近似固定集中时频可以把原本铺散开的故障能量收拢为清晰的横向时频条带方便人工观察故障特征的有无当设备转速存在波动故障特征频率随时间发生漂移集中时频方法能够追踪时变脊线把时变故障分量的能量集中在弯曲的脊线上克服传统频谱方法在变转速下失效的短板。借助时频图上聚集后的能量分布既可以直接人工判别故障也可以从时频分布当中提取二维特征参数结合峭度、熵等指标量化评估轴承损伤程度实现故障的初步定量分析。4.2 强噪声与早期微弱故障识别工程现场采集的轴承数据往往信噪比偏低早期故障冲击幅值很小故障能量远低于噪声能量传统时频处理后故障特征被噪声完全覆盖难以识别。集中时频分析通过系数筛选与重分配将分散的故障能量聚合同时分散噪声能量拉大故障分量与噪声之间的特征差异提升微弱冲击的可辨识程度。但该场景下也存在固有局限。噪声会干扰瞬时频率估计结果错误的瞬时参数会造成时频重分配发生偏差带来伪影或者特征断裂。部分集中时频方法为追求聚集效果会过滤掉幅值较低的分量有可能直接抹除早期微弱故障信息反而造成漏诊。因此面向早期故障诊断不能一味追求时频图谱视觉锐利需要兼顾特征保留与聚集程度合理控制算法筛选阈值。4.3 复合故障特征分离轴承会出现多处同时损伤的复合故障振动信号中叠加多种故障对应的冲击成分不同故障的时频特征互相交叠特征分离难度大。具备重构能力的集中时频分析不仅输出时频图像还可以将不同时频脊线对应的信号分量逆向分解出来实现不同故障源的分量分离。将分解得到的单分量进一步做解调分析可以分别提取不同故障对应的冲击特征区分内圈、外圈同时损伤等复合故障模式。复合故障工况下各类故障冲击互相调制信号模式复杂不同故障的时频脊线距离很近集中时频处理容易出现分量混叠分量分离效果下降是当前实际应用中的难点。4.4 与智能故障诊断的结合当下大量智能故障诊断研究将时频图像作为深度学习网络输入完成故障分类、故障程度识别。集中时频分析输出的高聚集度时频图谱故障模式轮廓清晰背景噪声干扰更少相比于普通时频图谱能够降低网络的特征学习负担提升模型分类准确率尤其在噪声污染样本、变工况样本上可以获得更好的泛化性能。同时也存在现实问题集中时频变换会改变信号原始特征分布过度的时频增强有可能人为引入伪特征造成网络过拟合不同集中时频算法计算开销差异较大部分复杂重分配策略运算耗时高难以满足设备在线实时监测的轻量化需求。5 现存关键问题尽管集中时频分析在轴承故障诊断中展现出良好潜力从实验室仿真与台架试验取得不错效果但面向复杂工业现场依旧存在多项亟待解决的问题。第一算法参数依赖性较强。多数集中时频方法的处理效果受基础变换窗口、迭代次数、筛选阈值等参数影响显著参数选取不合适会造成时频聚集不足或是故障特征被过度剔除。目前缺少完备的参数自适应选取机制大多依靠人工反复调参现场使用门槛高。第二强噪声下瞬时特征估计精度下降。集中时频分析的核心依赖瞬时频率、群延迟的估计结果噪声会造成估计偏差进而引发时频重分配错位产生伪脊线、脊线断裂、特征畸变严重损害故障提取效果低信噪比条件下性能衰减明显。第三复合故障与剧烈变工况适应性不足。当转速快速变化、多故障耦合时信号组分复杂多变集中时频容易出现分量混叠不同故障特征无法有效拆分难以实现多损伤的同时识别。第四计算效率与工程落地之间存在矛盾。高阶、多迭代的集中时频改进算法可以获得更好的聚集效果但计算量随之上升。工业在线监测场景对处理时延有严格约束高精度集中时频算法难以直接部署在嵌入式监测硬件上算法轻量化研究不足。第五时频结果的评价体系不够完善。目前大多依靠图谱视觉效果评价聚集性能缺少贴合轴承故障任务的量化评价标准。部分方法追求图像视觉锐利却丢失真实故障瞬态信息图谱好看但诊断失效如何平衡聚集程度、特征保真、噪声抑制三者之间的关系还需要进一步研究。6 发展展望针对当前集中时频分析技术存在的短板结合轴承故障诊断工程需求后续可以从如下方向开展研究。1参数自适应集中时频方法研究。结合信号本身的峭度、信噪比、冲击特征实现窗口、阈值、迭代次数等关键参数的自适应选择减少人工干预提升算法对不同故障信号的普适能力降低现场使用门槛。2噪声鲁棒的瞬时特征估计策略。在集中时频框架内部嵌入降噪机制优化瞬时参数估计流程抑制噪声带来的估计误差在强干扰环境下维持时频脊线完整提升早期微弱故障提取能力。3面向复合故障、大波动变转速工况的算法改进。针对转速剧烈变化、多故障耦合的真实场景结合阶次分析思想优化集中时频重分配逻辑减少多分量混叠实现多故障特征的有效分离。4算法轻量化与工程化实现。对集中时频核心流程做裁剪、优化在保证故障特征提取能力前提下降低计算开销适配嵌入式监测终端推动该类技术从实验室台架走向工业设备在线监测。5集中时频与人工智能深度融合。探索时频聚集与深度学习的协同机制一方面利用集中时频为网络提供高质量二维特征另一方面借助人工智能反向辅助时频参数优化、脊线识别、伪特征甄别减少人为调参构建端到端的故障诊断链路。同时建立面向轴承故障任务的时频质量量化评价指标不再单纯依靠图像视觉效果评判方法优劣。6多技术融合研究。将集中时频分析与解卷积、盲源分离、共振解调等传统故障特征手段结合发挥不同技术的优势先对原始振动信号做预处理再开展集中时频增强进一步拓展复杂工况下的故障诊断能力。7 结论滚动轴承故障振动属于典型非平稳瞬态冲击信号传统时频分析存在能量弥散问题不利于早期故障识别。集中时频分析以时频能量聚集为核心目标通过对时频系数重分配、筛选、挤压压缩能量扩散范围获得特征锐利的时频表达能够有效凸显轴承故障冲击信息在恒定转速、变转速、噪声干扰场景下均具备良好的故障特征提取潜力。不同技术路线的集中时频方法各有优缺点重分配变换聚集效果突出但不支持信号重构同步压缩变换兼顾聚集与重构能力同步提取变换时频视觉效果优异但对噪声与参数更加敏感。现有集中时频分析依旧面临参数敏感、强噪声性能退化、复合故障适应性差、计算开销大等现实瓶颈。未来需要向着参数自适应、噪声鲁棒、轻量化、多技术融合方向发展推动集中时频分析更好服务于工业滚动轴承的状态监测与故障预警工作。第二部分——运行结果2.1 算例12.2 算例22.3 算例32.4 算例4部分代码clearload(vib_data2.mat)fs 12000; N 1200; % sampling frequency and pointstime (1:N)/fs; % time sequencefre (fs/2)/(N/2):(fs/2)/(N/2):(fs/2); % frequency sequencefre2fliplr(fre);datadata(1:N);[tfr Te]TET_Y(data,100);tfr2WT2(data,fs,600);TsSST_Y(data,100);x10.079; x20.083;y12.1; y24;%...................Fig. 12..................................figuresuptitle(Fig. 12);subplot(511)plot(time,data);axis([0 0.1 -1.25 1.25]);%xlabel(Time / s);ylabel(Amp / g);rectangle(Position,[x1 -1.25 x2-x1 2.5],EdgeColor,red,Linewidth,1);axes(position,[0.92,0.815,0.07,0.15]);plot(time,data);axis off;xlim([x1,x2]);ylim([-1.25,2.5]);第三部分——参考文献文章中一些内容引自网络会注明出处或引用为参考文献难免有未尽之处如有不妥请随时联系删除。(文章内容仅供参考具体效果以运行结果为准)​​​​​​第四部分——本文完整资源下载资料获取更多粉丝福利MATLAB|Simulink|Python|数据|文档等完整资源获取本文完整资源下载
返回列表