简介《电机噪音及振动分析》是一份面向电机维修、检测与运行维护人员的实用技术文档围绕电机修理后噪声与振动增大的典型问题展开。内容从机械、电磁、通风三条成因线索切入梳理转子不平衡、铁心松动、三相绕组不平衡、断笼、风扇缺陷等常见故障并给出通风噪声鉴别法、机械噪声鉴别法、电磁噪声鉴别法及突然断电法、改变电压法、对拖法等排查思路同时汇总紧固装配、选用专用轴承、优化绕组跨距、调整风扇间隙等降噪减振措施兼顾故障定位与检修落地。资源包共1个PDF文件约23KB篇幅精炼便于随时查阅所涉内容可与ISO相关标准对照。目前已有489人学习下载适合需要快速判断噪声来源、制定修理方案的电气检修人员与相关专业学习者参考。1. 一台噪音偏大的电机问题往往不在噪音本身生产线上一台 7.5 kW 异步电机空载听着还行一上负载就冒出 2 kHz 附近又尖又细的哨声换到另一台设备上噪音却是 100 Hz 附近沉闷的嗡嗡声手摸机脚能明显感到抖动。这两种声音背后的机理完全不同前者多半是气隙磁场力波激励定子铁芯后者往往指向转子动平衡或对中不良。电机噪音及振动分析这件事说到底就是把混在一起的声音按频率、按阶次、按传递路径拆开落到某个具体零件或某段电磁设计上。做设备点检的、做电机本体设计的、做 NVH 整改的都得靠这套拆解方法把听着不对变成这里不对。2. 电机噪音与振动的三类源头及其频谱指纹现场判断噪音来源第一步不是拆机而是看频谱长什么样。电机的噪声和振动基本可以归到三类电磁、机械含轴承、空气动力。三者的频谱形态差别很大电磁类集中在工频整数倍和齿槽相关的离散谱线上并且对电源频率、负载电流高度敏感机械类跟转速严格成正比升降速时谱线跟着平移空气动力类则是围绕叶片通过频率的宽带成分。先按这个框架归类再往下算具体特征频率排查范围能一次缩小到三分之一。2.1 电磁噪声气隙磁场力波与齿槽频率定转子开槽之后气隙磁导沿圆周周期性变化磁场能量也随之脉动产生随时间和空间分布的电磁力波。力波阶数越低、频率越接近定子铁芯的固有频率辐射出的噪声就越大。工程上最常查的几个频率是电源频率的 2 倍2f工频电机为 100 Hz来自磁致伸缩和径向力波、转子槽通过频率 R×frR 为转子槽数fr 为转频、定子槽通过频率 Z×fr以及它们的组合频率 |Z ± R|×fr。判断是不是电磁源有个很省事的现场试验把供电频率从 50 Hz 微调到 45 Hz 或 55 Hz变频器拖动时直接改给定如果那根可疑谱线跟着电源频率走而不跟着转速走基本可以定性为电磁噪声。反之谱线锁死在转速倍数上就往机械方向查。提示变频驱动的电机在 PWM 载波频率附近会出现一族边带载波频率 ± 整数倍调制频率这是驱动引入的不要误判成电机本体缺陷。把载波频率改一档再测一次谱线跟着搬家的就是它。2.2 轴承与机械噪声特征频率的解析计算轴承故障诊断的起点是算出四个特征频率外圈通过频率 BPFO、内圈通过频率 BPFI、滚动体自转频率 BSF、保持架频率 FTF。它们由转速和轴承几何尺寸唯一确定算出来之后在谱图上按这些频率找峰、找谐波、找以转频为间隔的边带就能定位到是哪个圈出了问题。import math def bearing_freqs(n_rpm, n_balls, ball_dia, pitch_dia, contact_angle_deg0.0): 计算滚动轴承四个特征频率(Hz) n_rpm : 轴转速 r/min n_balls : 滚动体个数 ball_dia : 滚动体直径 mm pitch_dia : 节圆直径 mm contact_angle_deg: 接触角深沟球轴承取 0 fr n_rpm / 60.0 # 轴旋转频率 r ball_dia / pitch_dia # 直径比 ca math.radians(contact_angle_deg) bpfo (n_balls / 2.0) * fr * (1 - r * math.cos(ca)) # 外圈故障特征频率 bpfi (n_balls / 2.0) * fr * (1 r * math.cos(ca)) # 内圈故障特征频率 bsf (pitch_dia / (2 * ball_dia)) * fr * (1 - (r * math.cos(ca)) ** 2) ftf 0.5 * fr * (1 - r * math.cos(ca)) # 保持架 return {FTF: ftf, BSF: bsf, BPFO: bpfo, BPFI: bpfi} print({k: round(v, 2) for k, v in bearing_freqs(1480, 8, 8.73, 39.0).items()})这段代码的核心是那个直径比 r 和接触角的余弦项它决定了内外圈特征频率的分裂方式。参数从轴承型号手册里查滚动体个数和直径、节圆直径、接触角四个量缺一个都算不准。有一点容易踩坑公式默认外圈静止、内圈随轴旋转如果装在皮带轮侧的轴承外圈跟着转比如托辊BPFO 和 BPFI 的表达式要对调。另外算出来的频率和实测往往差 1% 到 2%原因是实际转速不等于铭牌转速用键相或光电传感器直接测转速代入误差基本能压到 0.5% 以内。噪声/振动源频谱特征对转速的依赖现场验证手段电磁力波2f、Z×fr、R×fr 及组合频率部分跟随电源频率改变供电频率谱线是否跟走转子不平衡1×fr 及其小幅谐波严格正比升降速1×fr 幅值随转速平方增长对中不良2×fr 突出常伴 1×fr、3×fr严格正比联轴器两侧轴向振动明显增大机械松动多阶谐波底部噪声抬高正比但谱线发毛紧固后复测谐波族整体下降轴承局部缺陷BPFO/BPFI 及其谐波、转频边带正比包络解调后信噪比明显提升风扇/气动叶片数×fr宽带成分正比且随流量变化拆掉风扇罩试转对比2.3 空气动力噪声与结构共振如何区分风扇噪声的特征频率是叶片通过频率等于叶片数乘以转频通常在几百赫兹。它和机械类一样严格正比于转速区分点是拆掉风扇或改变风道阻力后这类成分会明显变化而轴承和不平衡的谱线几乎不动。结构共振则是另一回事——它不是声源而是放大器。当某个激励频率撞上机座、端盖或安装支架的某阶固有频率振动和噪声会被放大十几分贝。判断共振最直接的办法是敲击试验用力锤敲机脚采集响应看频响函数峰值在哪。如果这个峰值频率和频谱上那根突出的谱线重合问题就在结构而不是激励源。这种情况下只做动平衡或者换轴承往往收效有限得从加筋、改变支撑刚度或者错开共振区入手。做这类分析时频率分辨率要足够细Δf 最好小于 1 Hz否则峰值会糊成一片看不出到底是 96 Hz 还是 104 Hz。3. 振动信号采集与频谱分析的最小可复现流程很多人分析不出来问题不是算法不够高级而是采集环节就已经失真了。加速度计的量程、安装方式、采样率这三项定死之后后面再怎么处理都救不回来。这一章按先定参数、再写代码、最后调细节的顺序走一遍从一段原始加速度信号出发得到一根能用来下结论的谱线。3.1 采样率、量程与安装方式三个先决条件采样率按最低 2.56 倍最高分析频率来定工程上通常取 2.56 或 5 倍剩下的余量留给抗混叠滤波器。想看 5 kHz 以内的轴承高频共振采样率至少要 12.8 kHz取 25.6 kHz 更稳妥。量程要按预估的峰值来选加速度计满量程 50 g 结果实际冲击只有 2 g信噪比会差到没法看反过来量程不够会削顶削顶后的频谱会出现一大片虚假谐波比故障本身还显眼。安装方式是高频测量里最容易被忽视的一环。磁吸座在 1 kHz 以上就开始明显衰减而且会引入自己的共振峰蜂蜡和双面胶只适合临时低频测量。下面是几种常见安装方式的可用频率上限实际值随传感器和表面粗糙度变化仅作选型参考。安装方式可用频率上限参考适用场景钢制螺栓直接拧紧10 kHz 以上长期在线监测、轴承高频诊断粘接剂环氧/瞬干胶58 kHz短期精测、无法打孔的场合磁吸座12 kHz快速巡检、低频不平衡与对中手持探针500 Hz1 kHz粗筛不适合诊断轴承缺陷注意测点位置要固定并做标记同一台设备前后两次对比测点差 2 cm高频段的幅值就可能差一倍。水平和垂直方向都要测轴向振动是判断对中不良的关键很多现场只测一个方向结果把对中问题当成不平衡处理。3.2 用 Python 把一段加速度信号读成频谱拿到一段以 m/s² 为单位的加速度时域数据后先用 Welch 法做功率谱密度估计再转成 dB。相比直接做一次 FFTWelch 法分段加窗平均之后随机噪声被压下去离散谱线会凸显出来。import numpy as np from scipy import signal fs 25600 # 采样率 Hz x np.loadtxt(acc.csv) # 单列加速度单位 m/s^2 x x - x.mean() # 去直流否则 0 Hz 附近会拖尾 nper 4096 # 每段点数决定频率分辨率 nover nper // 2 # 50% 重叠提高平均次数利用率 win signal.windows.hann(nper, symFalse) # 汉宁窗抑制泄漏 f, Pxx signal.welch(x, fsfs, windowwin, npersegnper, noverlapnover, scalingdensity, detrendFalse) f, Pxx f[1:], Pxx[1:] # 丢掉 0 Hz 分量 Pxx_db 10 * np.log10(Pxx / 1e-6) # 参考 (1e-3 m/s^2)^2即 dB np.savetxt(spectrum.csv, np.c_[f, Pxx_db], delimiter,, headerf_hz,psd_db, comments) # 时域指标用于快速判断冲击性 rms np.sqrt(np.mean(x ** 2)) crest np.max(np.abs(x)) / rms print(fRMS{rms:.4f} m/s^2, 峰值因子{crest:.2f})逻辑上分三步先去直流防止 0 Hz 泄漏污染低频段再用汉宁窗加分段平均得到功率谱密度最后转 dB 并落盘。参数里最关键的是 nperseg频率分辨率 Δf fs / nperseg 25600 / 4096 6.25 Hz。这个分辨率用来分辨 BPFO 和 BPFI 是够的但如果要看电机转差频率50 Hz 电机满载滑差约 1%3%即 0.51.5 Hz6.25 Hz 完全不够必须把 nperseg 提到 32768 以上。峰值因子同时输出是为了给一个粗判正常运转的电机峰值因子多在 35超过 6 说明存在明显冲击这时候再去上包络解调。3.3 加窗、平均与重叠把谱线从噪声里拔出来窗函数的选择直接决定幅值准不准。汉宁窗适合连续随机信号幅值要修正平顶窗幅值精度最高适合做精确的幅值比较代价是频率分辨率变差相邻谱线容易被抹平矩形窗只在信号整周期采样时用比如阶次分析的角度域重采样之后。如果不确定就先跑一遍汉宁窗再跑一遍平顶窗看幅值差多少差值超过 10% 说明频率分辨率不足。平均次数建议 20 到 50 次。次数太少随机噪声压不住谱线毛刺多次数太多转速有波动的设备会出现谱线糊掉——因为每次采样的转速略不同同一根谱线在频率轴上位置轻微漂移平均之后就被抹宽了。重叠率取 50% 到 75%能在同样的数据长度下把平均次数翻两三倍算力允许的情况下这是性价比最高的一步。变频驱动的电机做变速分析时干脆放弃平均改用短时傅里叶变换做瀑布图看谱线随转速的走向会更直观。3.4 包络解调轴承早期缺陷的谱线长什么样轴承早期缺陷的冲击能量很弱直接看原始频谱那几根特征谱线常常被 1×fr 和背景噪声盖住。这时候用带通滤波把结构共振区取出来做希尔伯特包络再对包络求频谱信噪比能提升一截。from scipy.signal import butter, filtfilt, hilbert def envelope_spectrum(x, fs, band(2000, 6000), nperseg2048): b, a butter(4, [band[0] / (fs / 2), band[1] / (fs / 2)], btypeband) xb filtfilt(b, a, x) # 零相位带通避免冲击位置偏移 env np.abs(hilbert(xb)) # 希尔伯特包络 env env - env.mean() f, E signal.welch(env, fsfs, npersegnperseg) return f, E f, E envelope_spectrum(x, 25600, band(2000, 6000)) # 在 BPFO/BPFI 及其 2~5 阶谐波处找峰同时看有没有以 fr 为间隔的边带带通区间一般选在系统的高频共振带上26 kHz 是最常见的起点具体要看锤击试验测出的频响峰值在哪。filtfilt 做零相位滤波是必须的普通 filter 会引入相位延迟让冲击在时间轴上偏移做时间同步分析时对不上。包络谱上判断轴承缺陷看三件事特征频率处有没有峰、有几阶谐波、峰的两侧有没有以转频为间隔的边带。内圈故障因为随轴旋转受载方向周期性变化边带特别明显外圈故障位置固定边带相对弱但谐波更丰富。4. 变转速工况下的阶次分析与噪音源定位定转速下的频谱分析解决不了全部问题。风机、水泵、伺服驱动的电机大半时间在升降速用固定采样率的 FFT 去看谱线会随转速拖着尾巴根本没法对比。同时就算找到了那根可疑谱线也还要回答一个更麻烦的问题它是源头还是被结构放大的结果这一章处理这两件事。4.1 为什么升降速工况要看阶次而不是频率阶次是频率除以转频得到的无量纲量。转子不平衡永远出现在 1 阶对中不良在 2 阶风扇在叶片数阶。用阶次作横轴不同转速下的谱线就能对齐叠加一眼看出哪一阶在什么转速段冒头。这个特性叫阶次跟踪本质上要求采样在角度域上等间隔而不是在时间域上等间隔。实现路径有两条一个是硬件级的键相触发采样采集卡按外部脉冲触发另一个是软件重采样用等时间采的数据加一路键相信号事后插值到角度域。前者精度高后者灵活、成本低日常诊断用后者就够了。4.2 角度域重采样的实现做法是先用键相脉冲的上升沿时刻反推每个时刻对应的转角和转速再在角度轴上重新插值。下面的实现假设每转一个脉冲samples_per_rev 决定角度域的最高阶次取 256 意味着能分析到 128 阶对电机诊断绰绰有余。import numpy as np from scipy.interpolate import interp1d def order_track(x, t, key_t, samples_per_rev256): x,t : 等时间采样的信号与时间轴 key_t : 每转一个脉冲的上升沿时刻 返回角度轴(转)和角度域重采样信号 rev np.arange(len(key_t)) # 脉冲序号作为转数 ang np.arange(0, rev[-1], 1.0 / samples_per_rev) # 由 转数 - 时间 插值再取角度对应的时间 t_ang interp1d(rev, key_t, kindlinear, fill_valueextrapolate)(ang) xa interp1d(t, x, kindlinear, fill_valueextrapolate)(t_ang) return ang, xa ang, xa order_track(x, t, key_t, samples_per_rev256) X np.abs(np.fft.rfft(xa * np.hanning(len(xa)))) orders np.arange(len(X)) * (1.0 / (len(xa) / 256) if False else 1.0) / 256 orders np.arange(len(X)) / 256.0插值方式用线性就够追求精度可以换三次样条但要注意样条在高转速变化率下会产生过冲。角度域信号做完 FFT横轴直接就是阶次第 k 根谱线对应 k/256 阶。真正的坑在键相信号质量上脉冲有抖动或者一转多个脉冲没做分频处理重采样后的阶次谱会出现虚假的高阶成分看起来像啮合频率其实是采样本身的伪影。上机之前先画一张转速-时间曲线看看转速变化是否平滑突跳的地方基本就是脉冲丢失。4.3 锤击法模态与传递路径判断是源还是放大器确认了激励频率之后下一个问题是这个频率的振动为什么在某个测点特别大。锤击法模态测试能回答。用带力传感器的力锤敲击同时采力和响应算频响函数和相干函数。相干函数在共振区应该接近 1如果全频段都低于 0.8说明要么锤击能量不够要么响应信号被背景噪声污染了这次数据不能用。测点布置上判断源还是放大器的关键是比较源侧和响应侧的振动量级。源侧比如电机轴承座和响应侧比如机脚、基座、连接管路如果响应侧远大于源侧说明中间有放大环节光做动平衡没用。这时候做一次工作变形分析ODS看整机在某根谱线频率下的振型节点和腹点在哪加强筋该加在哪里就一目了然。环节典型测点判断依据处理方向激励源电机轴承座、端盖该点振动显著高于同机型基准动平衡、轴承更换、电磁整改传递路径机脚、底板、联轴器源与响应之比异常存在共振峰加筋、改支撑刚度、加隔振器响应体外壳、管路、护罩相干函数高、振型腹点集中局部加强、避开共振频率有一点要注意加隔振器不是万能的。隔振器本身有固有频率如果电机启停过程要穿越这个频率反而会在某段转速下振动陡增。选隔振器时系统固有频率要低于最低激励频率的 1/3 左右同时确认启停过程的穿越时间足够短。5. 从谱图到结论判据阈值与几个规避误判的技巧拿到一张干净的谱图之后最后一步是把它翻译成要不要停机的结论。宽带振动烈度是最常用的量化指标按 ISO 20816 系列原 ISO 10816分成 A、B、C、D 四个评价区。中小型电机常用的参考限值见下表实际判定务必核对标准原文和设备的安装方式刚性安装与弹性安装的限值不同测点取在轴承座还是机脚也会差一档。评价区含义中小型机组刚性安装参考限值RMS 速度A新交付机组≤ 1.4 mm/sB可长期运行≤ 2.8 mm/sC不宜长期运行需安排检修≤ 4.5 mm/sD可能造成损坏 4.5 mm/s单看宽带 RMS 有个明显缺陷它对早期轴承缺陷不敏感等到 RMS 涨上去往往已经磨到后期了。所以实际点检里更实用的是趋势 分频段的组合把频谱分成低频段101000 Hz看不平衡和对中、中频段13 kHz、高频段310 kHz看轴承和润滑分别记录趋势。任何一个频段在两周内涨幅超过 50%即使绝对值还在 B 区也值得停机检查。几个高频出现的误判。第一把 2 倍电源频率误判成机械问题判别方法是改一次供电频率谱线跟着电源走的就是电磁源这台电机的定子绕组或气隙均匀性有问题。第二把转速波动导致的谱线展宽当成故障判定前先看转速是否稳定波动超过 1% 就该改用阶次谱。第三把外壳或护罩的局部共振当成轴承故障护罩的固有频率往往就落在 24 kHz正是包络分析的带通区间测的时候顺手把护罩按住再采一次谱线消失就是它。第四忽略传感器自身的共振压电式加速度计在 2030 kHz 有安装共振如果分析带宽设得过高这个峰会被当成真实信号。一个比较实用的做法是把阶次谱和工作变形测点叠在一张图上横轴排阶次纵轴排测点位置用颜色深浅表示幅值哪一阶在哪个位置最亮源和路径一眼分得开。配合两周一次的固定测点复测同一根谱线的幅值变化率比绝对幅值更能说明问题很多故障都是先在变化率上冒头几周之后才反映到 RMS 上。本文还有配套的精品资源点击获取 电机噪音与振动分析：频谱、特征频率与阶次诊断实践
简介《电机噪音及振动分析》是一份面向电机维修、检测与运行维护人员的实用技术文档围绕电机修理后噪声与振动增大的典型问题展开。内容从机械、电磁、通风三条成因线索切入梳理转子不平衡、铁心松动、三相绕组不平衡、断笼、风扇缺陷等常见故障并给出通风噪声鉴别法、机械噪声鉴别法、电磁噪声鉴别法及突然断电法、改变电压法、对拖法等排查思路同时汇总紧固装配、选用专用轴承、优化绕组跨距、调整风扇间隙等降噪减振措施兼顾故障定位与检修落地。资源包共1个PDF文件约23KB篇幅精炼便于随时查阅所涉内容可与ISO相关标准对照。目前已有489人学习下载适合需要快速判断噪声来源、制定修理方案的电气检修人员与相关专业学习者参考。1. 一台噪音偏大的电机问题往往不在噪音本身生产线上一台 7.5 kW 异步电机空载听着还行一上负载就冒出 2 kHz 附近又尖又细的哨声换到另一台设备上噪音却是 100 Hz 附近沉闷的嗡嗡声手摸机脚能明显感到抖动。这两种声音背后的机理完全不同前者多半是气隙磁场力波激励定子铁芯后者往往指向转子动平衡或对中不良。电机噪音及振动分析这件事说到底就是把混在一起的声音按频率、按阶次、按传递路径拆开落到某个具体零件或某段电磁设计上。做设备点检的、做电机本体设计的、做 NVH 整改的都得靠这套拆解方法把听着不对变成这里不对。2. 电机噪音与振动的三类源头及其频谱指纹现场判断噪音来源第一步不是拆机而是看频谱长什么样。电机的噪声和振动基本可以归到三类电磁、机械含轴承、空气动力。三者的频谱形态差别很大电磁类集中在工频整数倍和齿槽相关的离散谱线上并且对电源频率、负载电流高度敏感机械类跟转速严格成正比升降速时谱线跟着平移空气动力类则是围绕叶片通过频率的宽带成分。先按这个框架归类再往下算具体特征频率排查范围能一次缩小到三分之一。2.1 电磁噪声气隙磁场力波与齿槽频率定转子开槽之后气隙磁导沿圆周周期性变化磁场能量也随之脉动产生随时间和空间分布的电磁力波。力波阶数越低、频率越接近定子铁芯的固有频率辐射出的噪声就越大。工程上最常查的几个频率是电源频率的 2 倍2f工频电机为 100 Hz来自磁致伸缩和径向力波、转子槽通过频率 R×frR 为转子槽数fr 为转频、定子槽通过频率 Z×fr以及它们的组合频率 |Z ± R|×fr。判断是不是电磁源有个很省事的现场试验把供电频率从 50 Hz 微调到 45 Hz 或 55 Hz变频器拖动时直接改给定如果那根可疑谱线跟着电源频率走而不跟着转速走基本可以定性为电磁噪声。反之谱线锁死在转速倍数上就往机械方向查。提示变频驱动的电机在 PWM 载波频率附近会出现一族边带载波频率 ± 整数倍调制频率这是驱动引入的不要误判成电机本体缺陷。把载波频率改一档再测一次谱线跟着搬家的就是它。2.2 轴承与机械噪声特征频率的解析计算轴承故障诊断的起点是算出四个特征频率外圈通过频率 BPFO、内圈通过频率 BPFI、滚动体自转频率 BSF、保持架频率 FTF。它们由转速和轴承几何尺寸唯一确定算出来之后在谱图上按这些频率找峰、找谐波、找以转频为间隔的边带就能定位到是哪个圈出了问题。import math def bearing_freqs(n_rpm, n_balls, ball_dia, pitch_dia, contact_angle_deg0.0): 计算滚动轴承四个特征频率(Hz) n_rpm : 轴转速 r/min n_balls : 滚动体个数 ball_dia : 滚动体直径 mm pitch_dia : 节圆直径 mm contact_angle_deg: 接触角深沟球轴承取 0 fr n_rpm / 60.0 # 轴旋转频率 r ball_dia / pitch_dia # 直径比 ca math.radians(contact_angle_deg) bpfo (n_balls / 2.0) * fr * (1 - r * math.cos(ca)) # 外圈故障特征频率 bpfi (n_balls / 2.0) * fr * (1 r * math.cos(ca)) # 内圈故障特征频率 bsf (pitch_dia / (2 * ball_dia)) * fr * (1 - (r * math.cos(ca)) ** 2) ftf 0.5 * fr * (1 - r * math.cos(ca)) # 保持架 return {FTF: ftf, BSF: bsf, BPFO: bpfo, BPFI: bpfi} print({k: round(v, 2) for k, v in bearing_freqs(1480, 8, 8.73, 39.0).items()})这段代码的核心是那个直径比 r 和接触角的余弦项它决定了内外圈特征频率的分裂方式。参数从轴承型号手册里查滚动体个数和直径、节圆直径、接触角四个量缺一个都算不准。有一点容易踩坑公式默认外圈静止、内圈随轴旋转如果装在皮带轮侧的轴承外圈跟着转比如托辊BPFO 和 BPFI 的表达式要对调。另外算出来的频率和实测往往差 1% 到 2%原因是实际转速不等于铭牌转速用键相或光电传感器直接测转速代入误差基本能压到 0.5% 以内。噪声/振动源频谱特征对转速的依赖现场验证手段电磁力波2f、Z×fr、R×fr 及组合频率部分跟随电源频率改变供电频率谱线是否跟走转子不平衡1×fr 及其小幅谐波严格正比升降速1×fr 幅值随转速平方增长对中不良2×fr 突出常伴 1×fr、3×fr严格正比联轴器两侧轴向振动明显增大机械松动多阶谐波底部噪声抬高正比但谱线发毛紧固后复测谐波族整体下降轴承局部缺陷BPFO/BPFI 及其谐波、转频边带正比包络解调后信噪比明显提升风扇/气动叶片数×fr宽带成分正比且随流量变化拆掉风扇罩试转对比2.3 空气动力噪声与结构共振如何区分风扇噪声的特征频率是叶片通过频率等于叶片数乘以转频通常在几百赫兹。它和机械类一样严格正比于转速区分点是拆掉风扇或改变风道阻力后这类成分会明显变化而轴承和不平衡的谱线几乎不动。结构共振则是另一回事——它不是声源而是放大器。当某个激励频率撞上机座、端盖或安装支架的某阶固有频率振动和噪声会被放大十几分贝。判断共振最直接的办法是敲击试验用力锤敲机脚采集响应看频响函数峰值在哪。如果这个峰值频率和频谱上那根突出的谱线重合问题就在结构而不是激励源。这种情况下只做动平衡或者换轴承往往收效有限得从加筋、改变支撑刚度或者错开共振区入手。做这类分析时频率分辨率要足够细Δf 最好小于 1 Hz否则峰值会糊成一片看不出到底是 96 Hz 还是 104 Hz。3. 振动信号采集与频谱分析的最小可复现流程很多人分析不出来问题不是算法不够高级而是采集环节就已经失真了。加速度计的量程、安装方式、采样率这三项定死之后后面再怎么处理都救不回来。这一章按先定参数、再写代码、最后调细节的顺序走一遍从一段原始加速度信号出发得到一根能用来下结论的谱线。3.1 采样率、量程与安装方式三个先决条件采样率按最低 2.56 倍最高分析频率来定工程上通常取 2.56 或 5 倍剩下的余量留给抗混叠滤波器。想看 5 kHz 以内的轴承高频共振采样率至少要 12.8 kHz取 25.6 kHz 更稳妥。量程要按预估的峰值来选加速度计满量程 50 g 结果实际冲击只有 2 g信噪比会差到没法看反过来量程不够会削顶削顶后的频谱会出现一大片虚假谐波比故障本身还显眼。安装方式是高频测量里最容易被忽视的一环。磁吸座在 1 kHz 以上就开始明显衰减而且会引入自己的共振峰蜂蜡和双面胶只适合临时低频测量。下面是几种常见安装方式的可用频率上限实际值随传感器和表面粗糙度变化仅作选型参考。安装方式可用频率上限参考适用场景钢制螺栓直接拧紧10 kHz 以上长期在线监测、轴承高频诊断粘接剂环氧/瞬干胶58 kHz短期精测、无法打孔的场合磁吸座12 kHz快速巡检、低频不平衡与对中手持探针500 Hz1 kHz粗筛不适合诊断轴承缺陷注意测点位置要固定并做标记同一台设备前后两次对比测点差 2 cm高频段的幅值就可能差一倍。水平和垂直方向都要测轴向振动是判断对中不良的关键很多现场只测一个方向结果把对中问题当成不平衡处理。3.2 用 Python 把一段加速度信号读成频谱拿到一段以 m/s² 为单位的加速度时域数据后先用 Welch 法做功率谱密度估计再转成 dB。相比直接做一次 FFTWelch 法分段加窗平均之后随机噪声被压下去离散谱线会凸显出来。import numpy as np from scipy import signal fs 25600 # 采样率 Hz x np.loadtxt(acc.csv) # 单列加速度单位 m/s^2 x x - x.mean() # 去直流否则 0 Hz 附近会拖尾 nper 4096 # 每段点数决定频率分辨率 nover nper // 2 # 50% 重叠提高平均次数利用率 win signal.windows.hann(nper, symFalse) # 汉宁窗抑制泄漏 f, Pxx signal.welch(x, fsfs, windowwin, npersegnper, noverlapnover, scalingdensity, detrendFalse) f, Pxx f[1:], Pxx[1:] # 丢掉 0 Hz 分量 Pxx_db 10 * np.log10(Pxx / 1e-6) # 参考 (1e-3 m/s^2)^2即 dB np.savetxt(spectrum.csv, np.c_[f, Pxx_db], delimiter,, headerf_hz,psd_db, comments) # 时域指标用于快速判断冲击性 rms np.sqrt(np.mean(x ** 2)) crest np.max(np.abs(x)) / rms print(fRMS{rms:.4f} m/s^2, 峰值因子{crest:.2f})逻辑上分三步先去直流防止 0 Hz 泄漏污染低频段再用汉宁窗加分段平均得到功率谱密度最后转 dB 并落盘。参数里最关键的是 nperseg频率分辨率 Δf fs / nperseg 25600 / 4096 6.25 Hz。这个分辨率用来分辨 BPFO 和 BPFI 是够的但如果要看电机转差频率50 Hz 电机满载滑差约 1%3%即 0.51.5 Hz6.25 Hz 完全不够必须把 nperseg 提到 32768 以上。峰值因子同时输出是为了给一个粗判正常运转的电机峰值因子多在 35超过 6 说明存在明显冲击这时候再去上包络解调。3.3 加窗、平均与重叠把谱线从噪声里拔出来窗函数的选择直接决定幅值准不准。汉宁窗适合连续随机信号幅值要修正平顶窗幅值精度最高适合做精确的幅值比较代价是频率分辨率变差相邻谱线容易被抹平矩形窗只在信号整周期采样时用比如阶次分析的角度域重采样之后。如果不确定就先跑一遍汉宁窗再跑一遍平顶窗看幅值差多少差值超过 10% 说明频率分辨率不足。平均次数建议 20 到 50 次。次数太少随机噪声压不住谱线毛刺多次数太多转速有波动的设备会出现谱线糊掉——因为每次采样的转速略不同同一根谱线在频率轴上位置轻微漂移平均之后就被抹宽了。重叠率取 50% 到 75%能在同样的数据长度下把平均次数翻两三倍算力允许的情况下这是性价比最高的一步。变频驱动的电机做变速分析时干脆放弃平均改用短时傅里叶变换做瀑布图看谱线随转速的走向会更直观。3.4 包络解调轴承早期缺陷的谱线长什么样轴承早期缺陷的冲击能量很弱直接看原始频谱那几根特征谱线常常被 1×fr 和背景噪声盖住。这时候用带通滤波把结构共振区取出来做希尔伯特包络再对包络求频谱信噪比能提升一截。from scipy.signal import butter, filtfilt, hilbert def envelope_spectrum(x, fs, band(2000, 6000), nperseg2048): b, a butter(4, [band[0] / (fs / 2), band[1] / (fs / 2)], btypeband) xb filtfilt(b, a, x) # 零相位带通避免冲击位置偏移 env np.abs(hilbert(xb)) # 希尔伯特包络 env env - env.mean() f, E signal.welch(env, fsfs, npersegnperseg) return f, E f, E envelope_spectrum(x, 25600, band(2000, 6000)) # 在 BPFO/BPFI 及其 2~5 阶谐波处找峰同时看有没有以 fr 为间隔的边带带通区间一般选在系统的高频共振带上26 kHz 是最常见的起点具体要看锤击试验测出的频响峰值在哪。filtfilt 做零相位滤波是必须的普通 filter 会引入相位延迟让冲击在时间轴上偏移做时间同步分析时对不上。包络谱上判断轴承缺陷看三件事特征频率处有没有峰、有几阶谐波、峰的两侧有没有以转频为间隔的边带。内圈故障因为随轴旋转受载方向周期性变化边带特别明显外圈故障位置固定边带相对弱但谐波更丰富。4. 变转速工况下的阶次分析与噪音源定位定转速下的频谱分析解决不了全部问题。风机、水泵、伺服驱动的电机大半时间在升降速用固定采样率的 FFT 去看谱线会随转速拖着尾巴根本没法对比。同时就算找到了那根可疑谱线也还要回答一个更麻烦的问题它是源头还是被结构放大的结果这一章处理这两件事。4.1 为什么升降速工况要看阶次而不是频率阶次是频率除以转频得到的无量纲量。转子不平衡永远出现在 1 阶对中不良在 2 阶风扇在叶片数阶。用阶次作横轴不同转速下的谱线就能对齐叠加一眼看出哪一阶在什么转速段冒头。这个特性叫阶次跟踪本质上要求采样在角度域上等间隔而不是在时间域上等间隔。实现路径有两条一个是硬件级的键相触发采样采集卡按外部脉冲触发另一个是软件重采样用等时间采的数据加一路键相信号事后插值到角度域。前者精度高后者灵活、成本低日常诊断用后者就够了。4.2 角度域重采样的实现做法是先用键相脉冲的上升沿时刻反推每个时刻对应的转角和转速再在角度轴上重新插值。下面的实现假设每转一个脉冲samples_per_rev 决定角度域的最高阶次取 256 意味着能分析到 128 阶对电机诊断绰绰有余。import numpy as np from scipy.interpolate import interp1d def order_track(x, t, key_t, samples_per_rev256): x,t : 等时间采样的信号与时间轴 key_t : 每转一个脉冲的上升沿时刻 返回角度轴(转)和角度域重采样信号 rev np.arange(len(key_t)) # 脉冲序号作为转数 ang np.arange(0, rev[-1], 1.0 / samples_per_rev) # 由 转数 - 时间 插值再取角度对应的时间 t_ang interp1d(rev, key_t, kindlinear, fill_valueextrapolate)(ang) xa interp1d(t, x, kindlinear, fill_valueextrapolate)(t_ang) return ang, xa ang, xa order_track(x, t, key_t, samples_per_rev256) X np.abs(np.fft.rfft(xa * np.hanning(len(xa)))) orders np.arange(len(X)) * (1.0 / (len(xa) / 256) if False else 1.0) / 256 orders np.arange(len(X)) / 256.0插值方式用线性就够追求精度可以换三次样条但要注意样条在高转速变化率下会产生过冲。角度域信号做完 FFT横轴直接就是阶次第 k 根谱线对应 k/256 阶。真正的坑在键相信号质量上脉冲有抖动或者一转多个脉冲没做分频处理重采样后的阶次谱会出现虚假的高阶成分看起来像啮合频率其实是采样本身的伪影。上机之前先画一张转速-时间曲线看看转速变化是否平滑突跳的地方基本就是脉冲丢失。4.3 锤击法模态与传递路径判断是源还是放大器确认了激励频率之后下一个问题是这个频率的振动为什么在某个测点特别大。锤击法模态测试能回答。用带力传感器的力锤敲击同时采力和响应算频响函数和相干函数。相干函数在共振区应该接近 1如果全频段都低于 0.8说明要么锤击能量不够要么响应信号被背景噪声污染了这次数据不能用。测点布置上判断源还是放大器的关键是比较源侧和响应侧的振动量级。源侧比如电机轴承座和响应侧比如机脚、基座、连接管路如果响应侧远大于源侧说明中间有放大环节光做动平衡没用。这时候做一次工作变形分析ODS看整机在某根谱线频率下的振型节点和腹点在哪加强筋该加在哪里就一目了然。环节典型测点判断依据处理方向激励源电机轴承座、端盖该点振动显著高于同机型基准动平衡、轴承更换、电磁整改传递路径机脚、底板、联轴器源与响应之比异常存在共振峰加筋、改支撑刚度、加隔振器响应体外壳、管路、护罩相干函数高、振型腹点集中局部加强、避开共振频率有一点要注意加隔振器不是万能的。隔振器本身有固有频率如果电机启停过程要穿越这个频率反而会在某段转速下振动陡增。选隔振器时系统固有频率要低于最低激励频率的 1/3 左右同时确认启停过程的穿越时间足够短。5. 从谱图到结论判据阈值与几个规避误判的技巧拿到一张干净的谱图之后最后一步是把它翻译成要不要停机的结论。宽带振动烈度是最常用的量化指标按 ISO 20816 系列原 ISO 10816分成 A、B、C、D 四个评价区。中小型电机常用的参考限值见下表实际判定务必核对标准原文和设备的安装方式刚性安装与弹性安装的限值不同测点取在轴承座还是机脚也会差一档。评价区含义中小型机组刚性安装参考限值RMS 速度A新交付机组≤ 1.4 mm/sB可长期运行≤ 2.8 mm/sC不宜长期运行需安排检修≤ 4.5 mm/sD可能造成损坏 4.5 mm/s单看宽带 RMS 有个明显缺陷它对早期轴承缺陷不敏感等到 RMS 涨上去往往已经磨到后期了。所以实际点检里更实用的是趋势 分频段的组合把频谱分成低频段101000 Hz看不平衡和对中、中频段13 kHz、高频段310 kHz看轴承和润滑分别记录趋势。任何一个频段在两周内涨幅超过 50%即使绝对值还在 B 区也值得停机检查。几个高频出现的误判。第一把 2 倍电源频率误判成机械问题判别方法是改一次供电频率谱线跟着电源走的就是电磁源这台电机的定子绕组或气隙均匀性有问题。第二把转速波动导致的谱线展宽当成故障判定前先看转速是否稳定波动超过 1% 就该改用阶次谱。第三把外壳或护罩的局部共振当成轴承故障护罩的固有频率往往就落在 24 kHz正是包络分析的带通区间测的时候顺手把护罩按住再采一次谱线消失就是它。第四忽略传感器自身的共振压电式加速度计在 2030 kHz 有安装共振如果分析带宽设得过高这个峰会被当成真实信号。一个比较实用的做法是把阶次谱和工作变形测点叠在一张图上横轴排阶次纵轴排测点位置用颜色深浅表示幅值哪一阶在哪个位置最亮源和路径一眼分得开。配合两周一次的固定测点复测同一根谱线的幅值变化率比绝对幅值更能说明问题很多故障都是先在变化率上冒头几周之后才反映到 RMS 上。本文还有配套的精品资源点击获取