简介：这份《电机噪音及振动分析》PDF 面向电机维修、检测与设备运维人员，聚焦电机运行中噪声与振动偏大这一常见故障，梳理从成因判别到整改落地的完整思路。内容按机械、电磁、通风三条线索展开：机械方面涉及转子固定键松动、风扇静平衡不足、转轴弯曲、铁心变形与装配不同心等；电磁方面覆盖三相绕组不平衡、绕组短路断路、断笼或端环开裂、改极时槽数配合不当、电源质量不佳等；通风方面关注风扇缺陷、风罩间隙与风路堵塞。文中还给出三类噪声的简易鉴别法，如去风扇或堵风口试运转、突然断电法、改变电压法、对拖法，并配套紧固螺栓、选用精密轴承、优化绕组型式与跨距、调整风扇外径等降噪减振措施。资源为单个 PDF，压缩包约 23KB，条目清晰，便于现场对照排查。已有 489 人学习下载，适合作为电机检修与故障诊断的速查参考。
1. 一台 3000 rpm 电机的啸叫，为什么要噪声和振动一起看
一台 15 kW 永磁同步电机在 2300 rpm 以上持续啸叫，售后先换了两套轴承，问题照旧。做电机噪音及振动分析时把噪声谱和机壳振动谱画到同一张阶次图上才看清：主峰在 12 阶，正好撞上定子二阶模态在 2400 rpm 附近的固有频率，是径向电磁力波激发的结构辐射噪声，和轴承没有半点关系。这类问题在伺服电机、电主轴、新能源驱动电机上都很常见，麻烦在于噪声和振动是两条测量通道、两套指标、两种传感器，只测声压级会丢掉激励源的位置，只测振动会漏掉辐射效率这一环。
这个方向要做的事情，是把"听上去很吵"还原成可量化的频率、阶次和幅值，再顺着激励源、传递路径、辐射面三段往下拆。适合做电控、结构、测试和故障诊断的工程师，尤其是需要给出可复现结论、而不是"换件试试"的人。下文按源头辨识、测试链路、数据解析、降噪验证四步展开，每一步都给到能直接跑的参数和代码。
2. 电机噪音与振动的三类源头：电磁力波、机械激励与空气动力
2.1 径向电磁力波：低阶力波为什么主导电磁噪声
气隙里的磁场不是理想正弦，定子绕组磁动势和转子磁场都含空间谐波，次数分别记为 ν 和 μ（以极对数为单位的空间次数）。两者在气隙中相互作用产生径向力波，阶次 r 由次数差或次数和决定，最常见的低阶力波来自定、转子一阶齿谐波：ν = Z1/p ± 1，μ = Z2/p ± 1，其中 Z1、Z2 为定转子槽数，p 为极对数。
定子铁心对这个力波的响应是强非线性的。在 r 阶力波作用下，铁心径向变形量近似与 1/(r²−1)² 成正比，r = 2 时的变形量比 r = 8 时高出一个数量级以上。极槽配合的噪声差异，本质就在这里——8 极 48 槽和 10 极 12 槽的最低阶力波阶次不同，辐射效率差出好几倍，跟绕组工艺关系不大。
力波频率方面，恒频供电的异步电机力波频率是电源频率的整数倍组合，并且随转差率略有偏移，在精细阶次谱上表现为稍微偏离整数阶的谱线。这一点可以用来和轴承特征频率区分：轴承频率与转频严格成固定比例，电磁力波频率跟电源频率绑定，跟转频只是近似比例。
变频驱动时，载波频率 f_c 及其边带 f_c ± k·f1 会在电流里产生高频谐波，对应力波落在载波频率附近。载波频率一旦落在 2~8 kHz 这段人耳最敏感的区间，啸叫就会格外刺耳。把 f_c 从 4 kHz 提到 8~10 kHz 是最直接的缓解手段，代价是开关损耗上升，驱动器温升和 EMC 表现需要重新评估。
2.2 机械与空气动力噪声的特征频率
机械类噪声的谱线有强规律性。设滚动体个数 n、滚动体直径 d、节圆直径 D、接触角 α、转频 fr = rpm/60（Hz），四个轴承特征频率可以直接算出来：
import math def bearing_fault_freqs(rpm, n, d, D, alpha_deg=0.0): """计算轴承四个特征频率，单位 Hz。 rpm: 轴转速 r/min；n: 滚动体个数；d: 滚动体直径 mm D: 节圆直径 mm；alpha_deg: 接触角（度） """ fr = rpm / 60.0 c = math.cos(math.radians(alpha_deg)) bpfo = n / 2.0 * fr * (1 - d / D * c) # 外圈通过频率 bpfi = n / 2.0 * fr * (1 + d / D * c) # 内圈通过频率 bsf = D / (2.0 * d) * fr * (1 - (d / D * c) ** 2) # 滚动体自转频率 ftf = 0.5 * fr * (1 - d / D * c) # 保持架频率 return {"BPFO": bpfo, "BPFI": bpfi, "BSF": bsf, "FTF": ftf}参数里最容易填错的是 d 和 D，务必查轴承手册的实测值而不是估算；接触角 α 对深沟球轴承一般为 0，角接触轴承必须按实际值代入。这四个频率在阶次谱上表现为固定的非整数阶，例如 BPFO 常见的 4.x、6.x 阶，是机械噪声里最容易识别的一类。内圈旋转时 BPFI 会带 ±fr 的边带，边带间隔就是转频，这也是判断损伤在内圈还是外圈的常用依据。
不平衡和转子偏心以 1X 转频为主，偏心还可能带出 2X，都属于与转频严格成比例的整阶。风扇的叶片通过频率 BPF = Z_b × fr，其中 Z_b 为叶片数，轴流风扇的 BPF 及其二、三次倍频是空气动力噪声的主要离散分量，其余是宽带涡流噪声，谱上没有明显尖峰，通常靠拆掉风扇前后对比总声压级来确认。
2.3 三类噪声源的辨识顺序：先断电，再扫转速
不要一上来就做精细谱。先用"断电惰转 + 变转速扫描"两个粗操作把三类源分开，再在对应频段用高分辨率谱锁细节，效率高得多。
具体步骤：
- 断电惰转记录。让电机从额定转速断电自由减速，同步录声压与振动。啸叫在断电瞬间消失，就是电磁源；断电后依然存在、并随转速扫过一个峰值，多半是机械激励或结构共振。
- 转速扫描。在 30%~100% 额定转速之间取 10~20 个稳态点，每点采 2~5 s。把所有谱按阶次而不是按频率对齐，看峰是随转速线性移动的阶次线，还是停在某频率不动的固定线。前者是电磁或机械激励，后者是结构模态或外部传递。
- 改变电气参数。调载波频率、调调制比、加减负载。电磁噪声对电流高频成分极度敏感，负载电流一变幅值就跟着变。
- 拆装对比。拆风扇、拆端盖分别复测，注意拆风扇会同时改变冷却条件，测点温度要同步记录。
| 特征 | 电磁噪声 | 机械噪声 | 空气动力噪声 |
|---|---|---|---|
| 谱线形态 | 与基频成整数倍的阶次线 | 转频倍频、轴承特征频率 | 叶片通过频率及倍频叠加宽带 |
| 随转速变化 | 阶次位置固定，幅值随电流变 | 阶次位置固定，幅值随转速上升 | 与转频成正比，幅值随风量变 |
| 断电惰转 | 瞬间消失 | 继续存在 | 继续存在 |
| 加载影响 | 幅值随负载电流明显变化 | 变化较小 | 变化较小 |
| 首选验证手段 | 电流谱与声压谱对照 | 振动测点加包络谱 | 拆风扇前后对比 |
3. 电机振动噪声测试链路：传感器、采集参数与转速跟踪
3.1 加速度计、麦克风与键相的安装要点
加速度计优先选 IEPE 型，通用灵敏度 100 mV/g（约 10 mV/(m/s²)），量程 ±50 g 能覆盖绝大多数中小型电机。低于 10 Hz 的低频段改用 500 mV/g 的高灵敏型号，代价是量程下降。安装方式直接决定可用频带：钢螺栓安装上限可到 10 kHz 以上，胶粘约 5 kHz，磁座只有 2 kHz 左右。测电磁噪声的高频力波时用磁座，数据基本不可信。
麦克风选 1/2 英寸自由场电容型，加防风罩，距电机表面 1 m、高度与轴心齐平，圆周均布 4~6 个方位点取能量平均。现场本底噪声至少要比目标声压级低 10 dB，否则测到的是环境不是电机。有条件就在半消声室或至少是硬地面、无强反射的空间做。
转速基准是阶次分析的前提。光电开关贴反光片每转一个脉冲，或用编码器 Z 相，脉冲沿的抖动要小于采样周期的 1%。没有可靠的键相信号，后面算出来的阶次分辨率全是虚的。
测点布置上，定子机壳圆周取水平径向、垂直径向、轴向三个方向，端盖和底座各加一个。轴向测点主要反映轴承和偏心，径向测点反映定子变形，两者分工不同，不要只测一个方向就下结论。
3.2 采样率、频率分辨率与窗函数的取值规则
采样率按分析带宽的 2.56 倍取，这是给抗混叠滤波器的过渡带留余量，不是随便定的经验数。频率分辨率由采样时长决定，Δf = fs/N = 1/T，想要 0.5 Hz 就必须采够 2 s。
| 参数 | 取值规则 | 3000 rpm 电机常用值 | 说明 |
|---|---|---|---|
| 采样率 fs | ≥ 2.56 × 分析带宽 | 25.6 kHz（带宽 10 kHz） | 留抗混叠过渡带 |
| 采样时长 T | 由阶次分辨率反推 | 2~5 s | T 越长 Δf 越小 |
| 频率分辨率 Δf | Δf = fs / N = 1 / T | 0.5 Hz（T = 2 s） | 稳态分析 |
| 阶次分辨率 ΔO | ΔO = 60 × Δf / rpm | 0.01 阶 @3000 rpm | 区分相邻边带 |
| 窗函数 | 稳态 Hanning，冲击 Flat-top | Hanning | 幅值精度优先选 Flat-top |
| 重叠率 | 50%~75% | 75% | 影响瀑布图密度 |
注意，阶次分辨率不是固定值，它随转速下降而变差。同样 0.5 Hz 的分辨率，在 600 rpm 时对应 0.05 阶，在 3000 rpm 时只有 0.01 阶。低速段的边带分析要相应加长采样时间。
3.3 用 Python 从时域数据算出阶次谱
拿到一段恒转速的加速度数据，第一步是把它从频率轴换到阶次轴：
import numpy as np from scipy.fft import rfft, rfftfreq def order_spectrum(acc, fs, rpm, f_max=3000): """恒转速段信号换算为阶次谱。 acc: 加速度时域序列，m/s^2 fs: 采样率，Hz rpm: 该段平均转速，r/min f_max: 分析上限频率，Hz """ x = np.asarray(acc, dtype=float) x = x - x.mean() # 去直流，避免 0 阶虚假大峰 n = x.size win = np.hanning(n) # Hann 窗抑制泄漏 X = np.abs(rfft(x * win)) * 2 / win.sum() # 单边谱幅值修正 f = rfftfreq(n, 1.0 / fs) order = f * 60.0 / rpm # 阶次 = 频率 × 60 / 转速 keep = f <= f_max return order[keep], X[keep]win.sum()做归一化，保证 1 m/s² 的正弦在谱上仍然读到 1 左右；乘以 2 是把双边谱折成单边谱。f_max建议控制在采样率的 0.4 倍以内，别踩进抗混叠滤波器的滚降区，那里幅值已经不准了。rpm 建议用键相信号在窗内取平均，而不是用变频器的设定转速——设定值和实际转速之间通常差几个百分点，在 3000 rpm 时就是几十转，足以让阶次线偏移一整个分辨率格。
4. 电机噪音振动数据的解析路径：时域指标、包络谱与阶次瀑布图
4.1 时域指标与包络谱判断轴承与偏心
时域指标是快速筛问题的手段，不同指标对应不同故障类型，选错了就看不出东西。
| 指标 | 计算方式 | 适用频段 | 典型用途 |
|---|---|---|---|
| 速度有效值 RMS | 10~1000 Hz 带通后算 RMS | 低频 | 整机振动烈度分级 |
| 加速度 RMS | 全带宽 RMS | 宽频 | 轴承、齿轮高频能量 |
| 峰值因子 | 峰值 / RMS | 宽频 | 早期冲击类缺陷 |
| 峭度 | 四阶矩归一化 | 宽频 | 轴承点蚀 |
| 包络谱峰值 | 带通加希尔伯特解调 | 高频共振带 | 轴承故障频率定位 |
轴承早期损伤的能量集中在高频共振带，直接在原始谱上看，故障特征频率往往被转频和电磁峰淹没。正确做法是先带通到共振带，再做包络解调：
import numpy as np from scipy.signal import butter, filtfilt, hilbert, welch def envelope_spectrum(x, fs, band=(2000, 8000), nperseg=8192): """带通 + 希尔伯特包络 + 包络谱，用于轴承早期损伤定位。""" b, a = butter(4, np.array(band) / (fs / 2.0), btype='band') xb = filtfilt(b, a, x) # 带通到轴承共振带，滤掉转频等低频成分 env = np.abs(hilbert(xb)) # 希尔伯特变换取包络 env = env - env.mean() f, Pxx = welch(env, fs, nperseg=nperseg) return f, Pxxband不能凭感觉填，要先看原始谱里 2~10 kHz 区间哪一段有明显的共振峰隆起，把带通设在那段上。nperseg决定包络谱的分辨率，8192 点在 25.6 kHz 采样下对应约 3 Hz，足以把 BPFO 和它的边带分开。解调后把峰值频率除以转频换成阶次，再和 2.2 节算出的四个特征频率比对，命中率高下立判。
4.2 阶次切片与瀑布图定位共振转速带
扫转速采到的数据最有价值的呈现方式是阶次切片：固定几个关心阶次，看它的幅值随转速怎么变，共振带会表现为一个明显的包络峰。
import numpy as np from scipy.fft import rfft, rfftfreq def order_slice(acc, fs, rpm_series, orders, win_sec=0.2, overlap=0.75): """按时间滑窗提取指定阶次的幅值随转速变化。 orders: 关心的阶次列表，如 [1, 2, 6, 12] win_sec: 分析窗长，决定阶次分辨率 overlap: 窗重叠率 """ n = int(win_sec * fs) step = max(1, int(n * (1 - overlap))) win = np.hanning(n) result = {o: [] for o in orders} rpm_axis = [] for start in range(0, len(acc) - n, step): seg = acc[start:start + n] rpm = float(np.mean(rpm_series[start:start + n])) if rpm <= 0: continue seg = seg - seg.mean() X = np.abs(rfft(seg * win)) * 2 / win.sum() f = rfftfreq(n, 1.0 / fs) order = f * 60.0 / rpm for o in orders: idx = int(np.argmin(np.abs(order - o))) # 最近谱线近似目标阶次 result[o].append(X[idx]) rpm_axis.append(rpm) return np.array(rpm_axis), {k: np.array(v) for k, v in result.items()}win_sec是关键参数：0.2 s 的窗在 3000 rpm 下给出约 0.01 阶的分辨率，够用；如果低转速阶次线密集，就要加长到 0.5~1 s，但时间分辨率会下降，快速升速工况下峰值可能被抹平。overlap取 0.75 让瀑布图更连续，代价是计算量翻四倍。得到结果后画成横轴转速、纵轴幅值的曲线，包络峰对应的转速就是共振点，再回头用恒转速工况确认一次。
4.3 结构模态与运行模态的交叉验证
阶次切片只能告诉你"在哪个转速共振"，不能告诉你"哪个结构在共振"。锤击法用安装了力传感器的力锤敲机壳，同时用加速度计测响应，算频响函数能得到固有频率、阻尼比和振型。注意加速度计要避开模态节点，否则那一阶模态在数据里根本不出现。
运行模态分析（OMA）不需要人工激励，直接用整机运行时的多通道振动做峰值拾取或随机子空间辨识，得到的是带装配预紧、带温度、带负载的真实模态。两者一对比，常常发现锤击测到的固有频率比运行模态高 5%~10%，原因是转子装配后刚度变化。
判据很直接：任何激励频率（基频及其倍频、齿谐波力波频率、载波频率及其边带）与固有频率的偏离小于 ±15%，就要考虑改结构或改激励。改哪个先看成本：变频参数可以现场改，极槽配合要动模具。
5. 降噪减振的验证闭环：参数调整、结构避让与报告固化
落地时的调整顺序一般是先电磁后结构，因为改载波频率和调制方式的代价最小。
电磁侧的首选动作是上调载波频率，把人耳最敏感频段的离散峰挪走；如果驱动器开关损耗受限，可以退而用随机 PWM 或周期频率调制，把载波能量摊成宽带，代价是高频段总能量略升。再往下才动极槽配合和斜极，目标都是提高最低阶径向力波阶次，斜极角一般取一个定子齿距对应的机械角（360/Z1 度），能有效压齿谐波力波，但轴向电磁力会增大、平均转矩略降，需要重新做电磁校核。
结构侧的三板斧是加厚定子轭部、机壳加筋、约束层阻尼贴片。前两者抬升模态频率或把它移出激励带，阻尼贴片则是不改频率、只压共振峰的 Q 值。隔振垫选型有个硬约束：隔振系统的固有频率不高于最低激励频率的 1/3，取 1/3~1/5 区间，否则不仅不隔振，还会放大。
| 手段 | 作用对象 | 关键参数 | 验证指标 |
|---|---|---|---|
| 载波频率上调 | 变频电磁噪声 | f_c 从 4 kHz 调到 10 kHz | 载波频段声压级 |
| 斜极/斜槽 | 齿谐波力波 | 斜极角约一个定子齿距 | 主阶次切片幅值 |
| 定子模态避让 | 结构共振 | 固有频率偏离激励 ±15% 以上 | 运行模态峰值 |
| 隔振垫 | 结构传递路径 | 系统固有频率 ≤ 激励的 1/3 | 底座振动烈度 |
验证必须用同一套采集脚本、同一组测点、同一批转速点做 A/B 对比，否则温度漂移和装配差异会盖过改进量。对比看三个数：总声压级 dB(A)、主阶次切片幅值、机壳振动速度有效值。经验上总声压级降 3 dB(A) 相当于声能量减半，但如果阵列数据里总声压级降了 3 dB(A)、主阶次却只降了 1 dB，通常说明能量转移到了相邻阶次或更高频段，这时候要回头看瀑布图，而不是盯着单点数值继续调——真正需要改的可能是另一个结构模态，而不是手上这个参数。
本文还有配套的精品资源，点击获取
