深入理解振动分析

可视化演示

展示系统处理流水线、实时性能指标曲线及多方案对比分析，包含数据流转、吞吐延迟监控与方案评估。

振动基础概念：机器的脉搏

王师傅把听诊器挪到减速机轴承座正上方，对实习生说：“你听——每转一圈，机器就有一次心跳。心跳的强弱、快慢、节律，全藏在这股振动里。” 任何旋转设备都存在振动，它源于转子质量不平衡、轴不对中、齿轮啮合冲击、轴承微观缺陷等。当零件状态变化时，“心跳”也会随之改变，这就是振动分析能诊断疾病的根本前提。

要解读振动，先认识它的三个基本要素。频率，即每秒振动的次数，单位赫兹（Hz），好比心跳的节奏快慢；振幅，指振动幅度的大小，好比心跳的力度；相位，描述振动到达峰值的时间点，好比心跳的节拍位置。三要素各管一个维度，就像医生同时关注心率、血压和心律。

工程上常用位移、速度、加速度三种量来度量振动，三者通过频率相互换算。低频时位移明显，中频时速度主导，高频时加速度突出。轴承早期点蚀产生的冲击是数千赫兹的高频振动，必须用加速度传感器才能捕捉——就像心跳中的细微杂音，要用高灵敏度的听诊器头才能分辨。

王师傅常说：“波形图上那条正弦线，就是机器的'心电图'。学会读它，才算摸到了听诊的门把手。”

信号采集规范：听诊器的放置艺术

医生听诊从不乱放：心脏的声音在左侧锁骨中线处最清晰，传感器也一样。王师傅要求徒弟们“传感器必须贴在轴承座正上方，用螺栓刚性固定，方向对准受力线”。放错了位置，再好的分析也是白搭——你听到的可能只是“隔壁床病人的心跳”。

信号采集的核心法则是采样定理：要让采样后的数字信号不丢信息，采样率必须至少是信号最高频率的两倍。工程上为了留出抗混叠滤波器的过渡带，通常取最高频率的 2.56 倍以上。比如要分析 10 kHz 以内的轴承振动，采样率至少设为 25.6 kS/s。此外，采样时长决定了频率分辨率：采得越久，频域上看得越细。

传感器选型也有讲究：加速度计体积小、频响宽，是振动监测的通用选择。安装方式从优到劣依次为：螺栓安装、胶粘、磁吸座、手持探针。磁吸座方便但会衰减高频，测轴承早期故障时高频成分至关重要，此时必须用螺栓或胶粘。同时要避开薄壁罩壳、悬臂等柔性位置，并做好电缆固定与接地，防止 50 Hz 工频干扰混入信号。

“听诊器摆不对，诊断全白费。” 王师傅把这句话挂在嘴边——采集端做扎实，后面所有花哨算法才有意义。

时域特征分析：波形的“心电图”

拿到振动信号，王师傅第一件事不是算频谱，而是直接看波形。健康的设备波形平顺如呼吸，轴承出现剥落时，波形上会猛地冒出一个个尖锐的“扎刺”。这份原始的波形图就是设备的“心电图”。

为了把波形变成可量化、可报警的数字，工程师提炼出一组时域特征。均值反映信号中的直流分量，通常趋近于零；均方根值（RMS）反映振动能量，是设备整体状态的通用指标；峰值反映最大瞬时冲击；峰值因子是峰值与 RMS 的比值，能敏锐察觉“偶发的重击”；峭度则刻画波形分布的尖陡程度，正常信号的峭度接近 3，一旦超过 4，说明信号中出现了大量冲击性成分。

在实际诊断中，RMS 适合跟踪缓慢发展的磨损，它像体温计一样稳定可靠；但早期轴承点蚀产生的冲击持续时间极短，RMS 来不及反应，峭度和峰值因子却能第一时间捕捉到“尖刺”。因此，完善的监测系统往往同时保存 RMS 和峭度两条趋势线，互为补充。

“心电图上有尖刺，下一步就要查是哪个器官在闹脾气。” 王师傅知道，时域只是粗筛，真正的定位要靠频域。

频域特征分析：把混音拆成独奏

王师傅带实习生听交响乐排练，然后问：“整个乐队一起奏，你能听出是哪把小提琴跑了调吗？” 当然不能。但如果把每个人的声部单独放出来，问题立刻暴露。FFT（快速傅里叶变换）就是振动分析里的“分轨师”——它把杂乱叠加的时域波形，拆解成一个个独立的频率分量。

FFT 的输入是 N 个时域采样点，输出是 N 个频域复数值，其中第 k 个值 X(k) 代表着频率为 \(k \Delta f\) 的那个“纯音”有多大。把频率作为横轴、幅值作为纵轴画出来，就是频谱图。频谱上的每一个峰，都对应一种周期性振动源：1 倍转频的峰提示转子不平衡，2 倍转频的峰提示轴不对中，啮合频率的峰提示齿轮问题，特定故障频率的峰则能直接锁定轴承缺陷。

解读频谱有三条经验：一看峰值对应的频率是否与特征频率表吻合；二看谐波家族是否齐整（基频、2 倍、3 倍逐级排列）；三看峰周围是否有边带（以转频为间隔的“卫星峰”）。谐波家族常表示松动或摩擦，边带则往往意味调制——某个零件的缺陷在“打拍子”。

“从合奏里拎出跑调的音，频谱是最锋利的耳朵。” 但王师傅还有一句悄悄话：有些极早期的故障音太轻，直接听频谱听不见，得用上“隔墙听诊”的绝活。

包络分析：隔墙听诊的顺风耳

早期轴承缺陷像隔壁房间的窃窃私语，声音很轻，还被厚厚的墙壁——也就是正常振动——挡着。直接对原始信号做 FFT，只能看到一大片高频“轰响”，听不出里面藏着的节奏。王师傅的绝活是拿一个“听诊杯”扣在墙上：先把高频共振区域挑出来，再把它的轮廓描出来，最后听这个轮廓的起伏节奏。

这套绝活在技术上叫包络分析。缺陷滚动体每撞击一次损伤点，就会激起轴承结构数千赫兹的高频共振，而撞击本身按故障特征频率发生。因此高频共振信号的强弱，是被故障频率“调制”的。分析流程是：先用带通滤波器切出高频共振区，再用希尔伯特变换提取信号的包络，最后对包络做 FFT，得到包络谱。包络谱中突出的低频峰，正是轴承的故障特征频率。

包络分析的价值在于“剥离噪声、放大节奏”。外圈故障特征频率 89.6 Hz 对应的振动幅值原本埋没在 3000 Hz 的共振噪声里，包络谱却能把它清晰地托出来。正因如此，工程界将包络谱视为轴承早期诊断的“金标准”——往往能比频谱法提前数月发现缺陷。

“隔墙听诊”让王师傅总能在设备明显坏掉之前，就听见病魔的脚步声。但这个方法的适用面也有边界，下一招则是为“转速不听话”的设备准备的。

阶次分析：变转速下的定位罗盘

救护车从身边呼啸而过时，警笛声会由尖变沉——频率随着相对运动而改变。设备转速变化时也有类似效应：频谱上的峰会左右“漂移”，同一个故障频率，900 转和 1500 转时出现在不同位置，让人难以比对。

王师傅应对这种情况的武器是阶次分析。思路很朴素：既然频率随转频变化，那就把每个频率都除以当时的转频，得到一个无量纲的比值——阶次。工频 \(f\) 除以转频 \(f_r\) 得到 1 阶，两倍频得到 2 阶。阶次是“转速无关”的坐标：无论轴转得多快，不平衡永远出现在 1 阶位置，不对中永远在 2 阶位置，某个轴承故障永远在它的特征阶次位置，像罗盘上的固定方位角。

实现阶次分析需要两路信号：一路是振动信号，另一路是键相传感器测得的转速脉冲。系统根据脉冲实时跟踪转速，将振动信号重采样到“等角度域”——每转一圈固定采多少个点，再对重采样信号做 FFT，就得到阶次谱。在风电齿轮箱、汽车变速器这类转速大范围波动的场合，阶次谱是唯一能稳定定位故障的手段。

有了这张“罗盘”，即使设备转速忽快忽慢，王师傅也能稳住心神，按图索骥找到故障的“经纬度”。

轴承故障特征：关节的四种“咔哒”

王师傅让徒弟联想自己的膝关节：跑步时，关节里的滚动体沿着内外圈“轨道”飞驰。膝关节出问题会发出特定节奏的咔哒声，轴承也一样——不同位置的损伤，会踩出不同的节拍。滚动轴承有四种典型缺陷位置：外圈、内圈、滚动体和保持架，它们各自的“咔哒”频率都可由几何参数和转速计算。

外圈固定在轴承座上，缺陷位置静止，滚动体每经过一次就撞击一次，脉冲规律稳定；内圈随轴旋转，缺陷位置不断改变，撞击信号要穿过滚动体和油膜才传到传感器，路径长短周期性变化，产生明显的调制现象；滚动体自身缺陷的脉冲频率与保持架转速有关，频响最高；保持架故障频率最低，通常只有转频的零点几倍。

例如某轴承 9 个滚动体、\(d/D \approx 0.15\)、转速 1500 r/min，算得 BPFO ≈ 89.6 Hz。王师傅在包络谱上看到 89.6 Hz 的尖峰，便果断写下诊断：“外圈剥落，建议两周内更换。” 听懂关节的节拍，设备就能免于“腿断停工”。

齿轮箱故障特征：牙齿的咬合节拍

齿轮箱像一嘴牙齿，齿与齿不断咬合传递动力。王师傅说，健康齿轮的“咬合节拍”均匀而安静；牙齿磨损、崩缺后，节拍变得刺耳、零乱。齿轮的每对齿啮合一次，就产生一次冲击，这个基础节拍叫啮合频率。

啮合频率等于齿数乘以所在轴的转频。正常齿轮的频谱以啮合频率及其谐波为主；当某个齿轮的齿面磨损时，啮合频率幅值升高，两侧会出现以该轴转频为间隔的边带，像主峰两边对称列队的小卫星。边带的间距直接告诉你：问题出在哪根轴上。齿面点蚀会产生宽带噪声，使频谱“地坪”抬高；断齿则每转一次产生一次巨大冲击，时域波形上出现周期性“大坑”。

王师傅总爱补一句：“边带是齿轮在喊'我在这根轴上'”。只要听懂牙齿的咬合节拍，就能在齿轮“掉牙”前安排它休假。

AI 应用场景

Python 演示

import numpy as np def hilbert_env(x):    """简易包络提取: 通过 Hilbert 变换取模"""    n = len(x)    X = np.fft.fft(x)    h = np.zeros(n)    if n % 2 == 0:        h[0] = h[n // 2] = 1        h[1:n // 2] = 2    else:        h[0] = 1        h[1:(n + 1) // 2] = 2    return np.abs(np.fft.ifft(X * h)) def simulate_bearing_signal(fs=25600, seconds=2, rpm=1500):    """模拟轴承外圈故障信号: 转频 + 故障冲击串 + 噪声"""    t = np.linspace(0, seconds, int(fs * seconds), endpoint=False)    fr = rpm / 60.0    bpfo = 3.585 * fr    shaft = 0.8 * np.sin(2 * np.pi * fr * t)    impulse = 0.35 * np.exp(-300 * (t % (1.0 / bpfo))) * np.sin(2 * np.pi * 3000 * t)    noise = 0.05 * np.random.randn(len(t))    return t, shaft + impulse + noise def detect_defect(signal, fs):    n = len(signal)    freqs = np.fft.rfftfreq(n, 1.0 / fs)    spec = np.abs(np.fft.rfft(signal)) / n    env = hilbert_env(signal)    env_spec = np.abs(np.fft.rfft(env)) / n    spec_peak = freqs[np.argmax(spec[5:]) + 5]    env_peak = freqs[np.argmax(env_spec[10:]) + 10]    return spec_peak, env_peak if __name__ == '__main__':    np.random.seed(42)    fs = 25600    _, sig = simulate_bearing_signal(fs)    spec_peak, env_peak = detect_defect(sig, fs)    theory = 3.585 * (1500 / 60.0)    print('FFT 频谱主峰: %.2f Hz (转频)' % spec_peak)    print('包络谱主峰: %.2f Hz' % env_peak)    print('理论外圈故障频率: %.2f Hz' % theory)

1. def hilbert_env(x):定义希尔伯特变换函数，用于提取信号包络。
2. h[1:n // 2] = 2在频域把正频率分量加倍，这是构造解析信号的关键步骤。
3. impulse = 0.35 * np.exp(-300 * (t % (1.0 / bpfo))) * np.sin(2 * np.pi * 3000 * t)模拟故障冲击串：每经过一个 BPFO 周期产生一次高频衰减冲击。
4. env = np.abs(np.fft.ifft(X * h))对解析信号取模，得到时域包络曲线。
5. env_spec = np.abs(np.fft.rfft(env)) / n对包络做 FFT，得到包络谱。
6. peak_idx = np.argmax(env_spec[10:]) + 10跳过低频区段，在包络谱中寻找故障特征峰。
7. print('包络谱主峰: %.2f Hz' % env_peak)输出检测到的故障频率，与理论值对比验证。

常见误区

面试题