轴承滚动体局部剥落故障诊断系统

资 源 简 介

该项目旨在通过MATLAB平台实现对滚动轴承中滚动体局部剥落故障的精准提取与诊断。系统核心功能包含对原始振动信号的读取、预处理、特征增强以及自动化的故障频率识别。针对滚动体局部剥落产生的周期性冲击信号，系统利用共振解调分析技术对高频细微冲击进行捕捉。通过希尔伯特变换（Hilbert Transform）获取信号的包络线，并进行包络谱分析，从而有效滤除噪声干扰并突出受故障调制的成分。系统将实时计算信号频率成分，并与理论计算得到的218.98Hz工作频率及其倍频进行精准比对。该程序支持自定义采样参数，能够处理

详 情 说 明

轴承滚动体局部剥落故障诊断系统

项目介绍

本系统是一个基于 MATLAB 平台的旋转机械健康监测工具，专门用于检测滚动轴承滚动体的局部剥落故障。通过模拟和分析振动信号，系统能够从强噪声背景中提取交变冲击特征，自动识别滚动体故障频率（BSF）及其倍频。该系统集成了信号模拟、时域统计、共振解调、希尔伯特变换包络分析以及自动诊断逻辑，为轴承的预测性维护提供技术支持。

功能特性

• 信号模拟仿真：能够生成包含周期性指数衰减冲击特征的滚动轴承故障信号，并模拟真实环境下的高斯白噪声干扰。
• 多维指标评估：计算有效值（RMS）、峭度值（Kurtosis）、峰值因子等时域特征，初步判断信号的冲击特性。
• 共振解调分析：利用带通滤波器提取故障冲击激起的结构高频共振成分，实现高低频解耦。
• 包络谱特征提取：通过希尔伯特变换（Hilbert Transform）提取信号包络，通过快速傅里叶变换（FFT）定位故障调制频率。
• 自动逻辑诊断：系统自动对比理论故障频率与实测频谱峰值，结合峭度准则，给出明确的故障诊断结论。
• 多维可视化展示：实时生成原始时域波形、包络线、包络频谱及局部精细频谱对比图。

使用方法

1. 环境准备：启动 MATLAB 软件（建议安装有信号处理工具箱 Signal Processing Toolbox）。
2. 执行程序：运行主脚本。程序将自动执行从信号生成到结果输出的所有流程。
3. 查看分析结果：

- 图形窗口：观察原始信号与包络线的重合情况，确认包络谱中红色虚线（理论故障频率）是否与蓝色谱带峰值对齐。 - 命令行终端：查看诊断报告，包括时域统计指标、捕捉到的倍频点以及最终的故障诊断结论。

1. 参数调整：根据实际工业需求，可在代码起始位置修改采样频率（fs）、故障特征频率（f_fault）或信噪比（snr）等参数。

系统要求

• 开发环境：MATLAB R2016b 或更高版本。
• 必备工具箱：Signal Processing Toolbox（用于执行希尔伯特变换、滤波器设计及寻找峰值）。

核心功能实现逻辑

系统流程严格遵循以下逻辑步骤：

1. 参数初始化

设定采样频率为 20000Hz，信号长度为 1 秒。定义滚动体故障特征频率为 218.98Hz，结构共振频率设为 3000Hz。

1. 故障信号合成

基于冲击力学模型，利用循环叠加算法生成一系列指数衰减脉冲。具体公式结合了阻尼系数和正弦振荡，模拟滚动体经过剥落坑时激发的瞬态振动。随后加入 -5dB 的强高斯白噪声以模拟复杂工况。

1. 时域特征计算

计算信号的均方根值（RMS）以反映能量大小，计算峭度值（Kurtosis）以衡量信号的冲击尖锐程度。针对轴承故障，重点监控峭度值是否超过正常范围（3.0）。

1. 共振解调预处理

设计一个 4 阶巴特沃斯（Butterworth）带通滤波器，通频带范围设定在共振频率 3000Hz 左右各 500Hz。该步骤旨在滤除低频振动干扰，保留包含故障信息的共振频段。

1. 希尔伯特包络分析

对滤波后的信号进行希尔伯特变换，解析其实时包络。对包络信号进行去直流处理（Detrend），并执行 FFT 变换得到包络谱。

1. 特征比对与自动化诊断

在 0-1000Hz 的范围内搜索包络谱峰值。程序自动检索理论故障频率（218.98Hz）的前 1-3 倍频。如果实测频率与理论频率的误差在 5Hz 以内，或信号峭度大于 4，系统则触发报警，判定存在剥落损伤。

关键算法说明

• 指数衰减冲击模型：模拟剥落坑引起的周期性宽频冲击，反映了由于滚动体通过载荷区表面缺陷时引发的物理过程。
• Butterworth 滤波器：采用 4 阶设计，在保证阻带衰减的同时具有较平坦的通带特性，适合提取特定频段的共振响应。
• 希尔伯特变换 (Hilbert Transform)：数学上构造信号的解析形式，是实现共振解调的核心，能够解调出被高频载波调制的低频冲击特征。
• 峰值匹配算法：通过设定幅值阈值（最大值的 30%）自动筛选显著谱峰，并利用最小误差原则执行实测频率与理论频率的关联匹配。