4、信号采集与处理:采样率、抗混叠滤波、FFT频谱分析、时域与频域特征提取

振动信号采集与处理是设备异常检测联动PLC系统的感知前端,其质量直接决定了后续报警判断与停机保护动作的可靠性。工业现场振动信号往往混杂着工频干扰、电磁噪声、机械冲击脉冲等多种成分,若采样率设置不当或未进行抗混叠滤波,FFT频谱分析将产生虚假频率成分,导致PLC误报或漏报。本章围绕“采样率设定→抗混叠滤波→FFT频谱分析→时域/频域特征提取”这一完整链路展开,重点解决工程实践中信号采集参数如何确定、特征如何量化、以及如何将处理结果可靠地传递给PLC逻辑。

振动信号采集与处理链路 振动传感器 IEPE/压电式 灵敏度: 100mV/g 信号调理 电荷放大/恒流源 增益: 1~100 抗混叠低通滤波 截止频率 fc = fs/2.56 衰减斜率 ≥ -80dB/dec ADC采样 fs = 2.56 × fc 分辨率: 16~24bit FFT频谱分析 窗函数: Hanning/FlatTop 谱线数: 400~6400线 频率分辨率 Δf = fs/N 时域特征提取 峰值/峰峰值/均值 有效值(RMS)/峭度 波形因子/峰值因子 频域特征提取 1X/2X/3X幅值 边频带能量/谐波能量 频带能量占比 PLC联动决策 阈值比较 → 报警位置位 → 停机保护输出 关键参数: fs ≥ 2.56×fmax | 采样时长 T ≥ 1/Δf | 抗混叠滤波阶数 ≥ 8阶

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。

工程经验:对于PLC联动系统,采样率并非越高越好。过高的采样率会带来数据量激增、PLC扫描周期变长、存储压力增大等问题。建议根据设备最高分析频率(fmax)确定采样率:
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 良好 强干扰环境
实操建议:在PLC联动系统中,建议固定使用Hanning窗,并设置谱线数为1600线(N=4096)。这样在fs=25.6kHz时,Δf=6.25Hz,足以分辨大多数旋转机械的故障特征频率。

4.3 时域特征提取方法与工程应用

时域特征直接反映振动信号的幅值统计特性,计算简单、实时性好,适合PLC快速判断。常用特征包括:

  1. 峰值(Peak)与峰峰值(Peak-to-Peak):反映振动瞬时最大值,对冲击类故障敏感,但易受噪声干扰。
  2. 有效值(RMS):振动能量指标,ISO 10816标准采用RMS值评估设备状态,适合稳态振动监测。
  3. 峭度(Kurtosis):反映信号分布的尖锐程度,正常振动信号峭度≈3,轴承早期故障时峭度可升至5~10。
  4. 峰值因子(Crest Factor):峰值/RMS,用于区分磨损类故障(峰值因子下降)与冲击类故障(峰值因子上升)。
注意:时域特征虽然计算简单,但无法区分故障的具体频率成分。例如,不平衡与不对中都会导致RMS增大,但仅凭时域特征无法区分。因此,PLC联动系统应同时采集时域特征和频域特征,综合判断。

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报警阈值。

实际操作步骤:

  1. 将PLC模拟量输入模块采样率配置为25.6kS/s,启用内置抗混叠滤波器(截止频率10kHz)。
  2. 在PLC程序中增加FFT计算功能块,每2秒执行一次4096点FFT。
  3. 提取155Hz±5Hz频带幅值,设定报警阈值为5m/s²(正常值<1m/s²)。
  4. 增加峭度特征计算,设定报警阈值为5.0。
  5. 当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
总结:信号采集与处理是设备异常检测联动PLC系统的“眼睛”。采样率、抗混叠滤波、FFT参数、特征提取方法四者环环相扣,任何一环设置不当都会导致诊断结论失真。工程实践中应遵循“先确定分析频率范围→再定采样率→配置抗混叠滤波→选择FFT参数→提取多维特征→联动PLC逻辑”的流程,并针对每台设备建立特征频率档案,定期校准传感器和采集通道。