4、信号采集与处理:采样率、抗混叠滤波、FFT频谱分析、时域与频域特征提取
振动信号采集与处理是设备异常检测联动PLC系统的感知前端,其质量直接决定了后续报警判断与停机保护动作的可靠性。工业现场振动信号往往混杂着工频干扰、电磁噪声、机械冲击脉冲等多种成分,若采样率设置不当或未进行抗混叠滤波,FFT频谱分析将产生虚假频率成分,导致PLC误报或漏报。本章围绕“采样率设定→抗混叠滤波→FFT频谱分析→时域/频域特征提取”这一完整链路展开,重点解决工程实践中信号采集参数如何确定、特征如何量化、以及如何将处理结果可靠地传递给PLC逻辑。
4.1 采样率设定与抗混叠滤波原理
核心原理:采样率(Sampling Rate, fs)决定了振动信号能够被准确还原的最高频率范围。根据奈奎斯特采样定理,理论上fs必须大于信号最高频率成分的2倍,工程上通常取2.56倍以确保足够的频谱保真度。例如,若设备转速为3000rpm(50Hz),需要分析至10倍频(500Hz),则采样率应设为1280Hz以上。但实际工业振动监测中,往往需要覆盖到轴承故障特征频率(BPFO/BPFI)的数倍频,因此采样率通常设定在设备转频的100~200倍。
抗混叠滤波的必要性:当信号中存在高于fs/2的频率成分时,这些成分会“折叠”到低频区域,形成虚假的频谱峰值,严重干扰故障诊断。抗混叠滤波器(低通滤波器)必须在ADC采样之前将高于fs/2的频率成分衰减到可忽略水平。工程上要求滤波器在截止频率fc处的衰减不低于-80dB,且过渡带尽量陡峭。典型设计为8阶巴特沃斯或切比雪夫滤波器,截止频率设定为fc = fs/2.56。
fs = 2.56 × fmax
例如:分析至10kHz,则fs = 25.6kHz,取整为25.6kS/s或26kS/s。
4.2 FFT频谱分析的关键参数与窗函数选择
FFT分析流程:将时域信号截取为N点数据块,加窗函数后执行快速傅里叶变换,得到频谱幅值。关键参数包括:
- 频率分辨率(Δf):Δf = fs / N,决定了频谱中能区分的最小频率间隔。对于旋转机械,Δf应小于转频的1/4,否则相邻谐波无法分辨。
- 谱线数(Lines):常用400、800、1600、3200线,对应N=1024、2048、4096、8192点。
- 窗函数选择:Hanning窗适用于一般振动信号,幅值精度较好;FlatTop窗适用于需要精确幅值测量的场合(如振动烈度评估);矩形窗仅适用于瞬态冲击信号。
| 窗函数 | 主瓣宽度 | 旁瓣衰减 | 幅值精度 | 适用场景 |
|---|---|---|---|---|
| 矩形窗 | 最窄 | -13dB | 差 | 瞬态冲击、校准 |
| Hanning窗 | 较窄 | -31dB | 良好 | 一般振动监测 |
| FlatTop窗 | 较宽 | -44dB | 极佳 | 幅值精确测量 |
| Blackman窗 | 宽 | -58dB | 良好 | 强干扰环境 |
4.3 时域特征提取方法与工程应用
时域特征直接反映振动信号的幅值统计特性,计算简单、实时性好,适合PLC快速判断。常用特征包括:
- 峰值(Peak)与峰峰值(Peak-to-Peak):反映振动瞬时最大值,对冲击类故障敏感,但易受噪声干扰。
- 有效值(RMS):振动能量指标,ISO 10816标准采用RMS值评估设备状态,适合稳态振动监测。
- 峭度(Kurtosis):反映信号分布的尖锐程度,正常振动信号峭度≈3,轴承早期故障时峭度可升至5~10。
- 峰值因子(Crest Factor):峰值/RMS,用于区分磨损类故障(峰值因子下降)与冲击类故障(峰值因子上升)。
4.4 频域特征提取与故障频率对应关系
频域特征通过FFT频谱提取特定频率成分的幅值,与设备结构参数和运行状态直接关联。典型频域特征包括:
- 1X幅值(转频幅值):反映转子不平衡程度,通常以振动速度mm/s表示。
- 2X幅值:反映不对中、轴弯曲等故障。
- 谐波能量(2X+3X+...):反映松动、摩擦等非线性故障。
- 边频带能量:围绕啮合频率或固有频率两侧的边带,反映齿轮磨损、轴承故障。
- 高频频带能量(1kHz~10kHz):反映轴承早期点蚀、剥落等冲击性故障。
4.5 案例拆解:离心泵振动异常联动PLC停机保护
背景与约束:某化工厂一台离心泵(转速2970rpm,4极电机驱动,轴承型号6312)频繁出现振动超标报警,但PLC系统未触发停机,导致轴承最终抱死,造成非计划停车。事后分析发现,原系统仅采集振动速度RMS值,阈值设定为4.5mm/s(ISO 10816区域C上限),但RMS值在故障初期仅缓慢上升,未及时触发报警。
初始现象与目标:泵自由端轴承座振动加速度频谱中出现明显的BPFO(外圈故障特征频率)及其谐波,但速度频谱中1X幅值正常。目标是通过优化信号采集与处理策略,提前识别轴承早期故障并联动PLC停机。
关键参数与数据:
- 轴承6312参数:滚动体数Z=8,节径D=85mm,滚动体直径d=18.5mm,接触角α=0°
- BPFO = (Z/2) × (1 - d/D × cosα) × fr = 4 × (1 - 0.2176) × 49.5 = 155.2Hz
- 采样率设定:fs = 25.6kHz,谱线数1600,Δf = 6.25Hz
分析判断过程:采集加速度信号后,FFT频谱在155Hz、310Hz、465Hz处出现明显峰值,且边频带间隔为49.5Hz(转频),确认外圈故障。时域峭度值从3.2升至6.8,峰值因子从4.1升至7.5。而RMS速度值仅从3.2mm/s缓慢升至3.8mm/s,未达到4.5mm/s报警阈值。
实际操作步骤:
- 将PLC模拟量输入模块采样率配置为25.6kS/s,启用内置抗混叠滤波器(截止频率10kHz)。
- 在PLC程序中增加FFT计算功能块,每2秒执行一次4096点FFT。
- 提取155Hz±5Hz频带幅值,设定报警阈值为5m/s²(正常值<1m/s²)。
- 增加峭度特征计算,设定报警阈值为5.0。
- 当BPFO幅值超过阈值且峭度>5时,触发停机保护输出。
结果对比:优化后系统在轴承故障早期(RMS仅3.5mm/s时)即触发报警,提前72小时预警,避免了轴承抱死事故。PLC停机保护动作时间从原来的“RMS超标后停机”提前至“BPFO特征超标后停机”,设备维修成本降低约60%。
风险与代价:增加FFT计算导致PLC扫描周期从10ms延长至25ms,但该泵控制回路对实时性要求不高,可接受。需注意FFT计算对PLC CPU负载的影响,建议使用独立信号处理模块或高速计数器配合DSP处理。
复盘经验:单一RMS阈值无法捕捉早期轴承故障,必须结合频域特征。采样率、谱线数、窗函数的选择直接影响故障频率分辨率,工程中应针对每台设备建立特征频率库,并定期更新阈值。
4.6 常见误区与总结复盘
| 误区 | 后果 | 正确做法 |
|---|---|---|
| 采样率设置过低 | 高频故障成分混叠,频谱出现虚假峰值 | fs ≥ 2.56 × fmax,并验证抗混叠滤波效果 |
| 忽略抗混叠滤波 | 噪声折叠到低频,误报率升高 | 必须使用硬件低通滤波器,截止频率=fs/2.56 |
| 窗函数选择不当 | 幅值误差大,泄漏严重 | 一般用Hanning,精确幅值用FlatTop |
| 仅依赖时域RMS | 早期故障漏报 | 结合峭度、峰值因子、频域特征 |
| FFT点数不足 | 频率分辨率差,无法分辨相邻谐波 | Δf ≤ 转频/4,必要时增加N |