低转速、大直径的回转支承类轴承，它的故障冲击频率往往落在几十赫兹，时域波形和常规频谱都被轴频与齿轮啮合盖住，直接看不出毛病。我们做在线监测时，第一步是把 BPFO、BPFI、BSF、FTF 这四个特征频率算准，再用包络谱把藏在共振带里的冲击尖峰抓出来。这里以 KMF 的 KVE 08 0350 钢丝轨道轴承走一遍完整过程，方法本身对所有球轴承通用。
作者：Chans｜杭州融奕科技有限公司 应用工程师
一、从样本读几何，把节圆直径和接触角钉死
算特征频率只需要四个量：节圆直径 D_pw、钢球直径 D_W、球数 Z、接触角 α。KMF 的 KVE 08 0350 把前两个直接印在样本里。钢丝轨道轴承没有整体的内外套圈，钢球就排在四根硬化钢丝围成的圆形轨道上，所以钢球的节圆直径正好等于样本的滚动圆直径 D_L。样本 P14 的尺寸表给出 D_L = 350 mm，同一行还列出 d=338、D=362、H=12，这三数是钢丝轨道所在的结构包络，用来校核安装空间而不是算频率。我们实际装机前就是用 d/D/H 去核对轴肩和座孔，频率计算则只读 D_L，两套尺寸各管各的事，别混。
钢球直径在 P12 的选型表里按系列给出，不是每个型号单独列。KVE 08 这一行的球径是 8.000 mm，KVE 05、KVE 06 各系列球径不同，我们选型号时得先按前缀确认球径，不能拿错一档。接触角来自四点接触的结构：样本 P6 写明该轴承压力角（两条接触方向之间的夹角）为 90°，换算到单条接触线的接触角 α = 45°。
唯一卡壳的是球数 Z。样本没有列每型号的球数，我们按 cages 的节距估算：Z ≈ π·D_L / (D_W + 0.3)，代入 350 和 8.0 得到 Z ≈ 132 颗。这一步是假设，不是样本数字，所以后面用 Z=130/132/134 三档做敏感性，看它对结果影响多大。
载荷与转速上限也一并记下：P14 表给出动载 C = 22.4 kN、静载 C0 = 119.0 kN，脂润滑极限转速 nG = 546 min⁻¹（油润滑 654 min⁻¹）。这张尺寸表同时把动载、静载和两种润滑下的极限转速列在一行，我们后面定采样上限就靠 nG 这一列。
这几个数后面既用来校核，也用来定采样的频率上限。〔来源：德国KMF钢丝轨道轴承 P6/P12/P14〕
二、四个公式代入 KVE 08 0350
标准轴承缺陷频率公式（ISO 13373-1《机器状态监测与诊断—振动状态监测》给出的滚动轴承特征频率计算式）如下，其中 f_r = n/60 是轴频。四个频率各自对应一种失效位置：BPFO 看外圈滚道、BPFI 看内圈滚道、BSF 看钢球自转、FTF 看保持架（即节圆）整体打滑。公式里 (1 − D_W/D_pw·cosα) 这一项，是外圈固定后钢球绕节圆滚动、相邻两次压过同一处缺陷的几何比例；cosα 把接触角的倾斜折进转速投影，四点接触取 45° 时这一项只有 0.016 上下，所以 BPFO 和 BPFI 相差很小、阶次只差约 2。
BPFO = Z/2 · (1 − D_W/D_pw · cosα) · f_r BPFI = Z/2 · (1 + D_W/D_pw · cosα) · f_r BSF = D_pw/(2·D_W) · (1 − (D_W/D_pw · cosα)²) · f_r FTF = 1/2 · (1 − D_W/D_pw · cosα) · f_r把几何代进去：D_pw=350、D_W=8.0、α=45°、Z=132。先算共同项 D_W/D_pw·cosα = 8/350×0.7071 = 0.01617，于是 (1−0.01617)=0.98383、(1+0.01617)=1.01617，Z/2 = 66。下面用 Python 跑出额定转速 n=546 min⁻¹（f_r=9.10 Hz）和一档低转速 n=20 min⁻¹，并打印阶次。
importmath DL,DW,alpha,nG=350.0,8.0,45.0,546.0Z=round(math.pi*DL/(DW+0.3))# 样本未给球数，按 cages 节距估算c=math.cos(math.radians(alpha))ratio=DW/DLdeffreqs(Z,n):fr=n/60.0BPFO=Z/2*(1-ratio*c)*fr BPFI=Z/2*(1+ratio*c)*fr BSF=DL/(2*DW)*(1-(ratio*c)**2)*fr FTF=0.5*(1-ratio*c)*frreturnfr,BPFO,BPFI,BSF,FTF fr,BPFO,BPFI,BSF,FTF=freqs(Z,nG)print("Z=%d n=546 min^-1 fr=%.4f Hz"%(Z,fr))print("BPFO=%.2f Hz (阶次 %.3f)"%(BPFO,BPFO/fr))print("BPFI=%.2f Hz (阶次 %.3f)"%(BPFI,BPFI/fr))print("BSF =%.2f Hz (阶次 %.3f)"%(BSF,BSF/fr))print("FTF =%.3f Hz (阶次 %.4f)"%(FTF,FTF/fr))fr2,b2,i2,s2,t2=freqs(Z,20.0)print("n=20 min^-1: BPFO=%.2f BPFI=%.2f BSF=%.2f FTF=%.3f Hz"%(b2,i2,s2,t2))输出结果：BPFO = 590.89 Hz（阶次 64.933），BPFI = 610.31 Hz（阶次 67.067），BSF = 199.01 Hz（阶次 21.869），FTF = 4.48 Hz（阶次 0.492）。阶次不随转速变，所以现场拿到转速先乘阶次就能得到实时频率。低转速 n=20 min⁻¹ 时 f_r 只有 0.333 Hz，BPFO 降到 21.6 Hz、FTF 低到 0.16 Hz——这正是必须上包络谱的理由：这些尖峰和轴频、不对中混在同一段低频里，普通频谱根本分不出来。
球数假设带来的波动也有限。我们把 Z 在 130/132/134 之间挪一档，对比一下：
| Z | BPFO@546 (Hz) | 相对 132 档偏差 |
|---|---|---|
| 130 | 581.94 | −1.51% |
| 132 | 590.89 | 基准 |
| 134 | 599.85 | +1.51% |
Z 每差 2 颗，BPFO 约差 1.5%，远小于现场转速波动带来的误差，所以球数按节距估算够用，不必为了一颗半颗反复纠结。工程上我们不追求把 Z 算到整数，差一两颗对诊断结论没有影响，关键是把 D_L、D_W、α 这三个样本硬数读对、读准。
三、用包络谱反查 BPFO（仿真信号）
我们用一段仿真信号验证公式对不对。做法是按 BPFO 的周期在共振中心 4000 Hz 上叠加一串衰减冲击，再压一层轴频正弦和噪声当背景，模拟现场难辨的情况；然后做正交混频（I = x·cos2πf_res t，Q = x·sin2πf_res t），滑动平均低通得到包络，再对包络扫频找主峰。这段信号是仿真，没有动真实采集设备，只为验证算法能定位尖峰。正交混频的本质是把 4000 Hz 的载波搬到低频：乘法之后，冲击包络项留在直流附近，载波项跑到 8000 Hz 被低通滤掉，剩下的就是一串按 BPFO 间隔排列的包络脉冲。
importmath,random random.seed(7)fs,N,f_res,fr=12000.0,6000,4000.0,546.0/60.0BPFO=590.89Timp=fs/BPFO x=[0.8*math.sin(2*math.pi*fr*n/fs)+0.15*(random.random()*2-1)forninrange(N)]imp=[]k=0whileTrue:tk=k*Timp+random.uniform(-1,1)iftk>=N:breakimp.append(tk);k+=1fortkinimp:ki=int(tk)forninrange(max(0,ki-30),min(N,ki+120)):dt=(n-ki)/fsifdt>=0:x[n]+=math.exp(-dt*1200.0)*math.sin(2*math.pi*f_res*dt)I=[x[n]*math.cos(2*math.pi*f_res*n/fs)forninrange(N)]Q=[x[n]*math.sin(2*math.pi*f_res*n/fs)forninrange(N)]deflp(s):out=[0.0]*Nforninrange(N):a=max(0,n-14);out[n]=sum(s[a:n+1])/(n-a+1)returnout Il,Ql=lp(I),lp(Q)env=[math.hypot(Il[n],Ql[n])forninrange(N)]best_f,best_m=0.0,0.0forfiinrange(300,2001):f=float(fi);re=im=0.0forninrange(0,N,3):ph=2*math.pi*f*n/fs;re+=env[n]*math.cos(ph);im+=env[n]*math.sin(ph)m=math.hypot(re,im)ifm>best_m:best_m,best_f=m,fprint("包络谱主峰 = %.1f Hz, 理论 BPFO = %.2f Hz, 误差 %.2f%%"%(best_f,BPFO,100*(best_f-BPFO)/BPFO))输出：包络谱主峰 591.0 Hz，与理论 BPFO 590.89 Hz 误差 0.02%；2×、3× 谐波分别在 1181.8 Hz、1772.7 Hz。实测冲击平均间隔 20.31 samples，对应 590.9 Hz，和理论对上。这说明公式算出来的频率，在包络谱上确实能定位到尖峰；反过来，现场拿到一个 591 Hz 的包络尖峰，就能锁定外圈缺陷。谐波成串出现，说明冲击是周期性、重复性的，而不是随机磕碰，这点对区分真故障和偶发干扰很关键。
四、把采样参数一次钉死
光知道频率还不够，采样设错了包络尖峰照样抓不到。我们按额定点定参数：
- 共振带取 3600–4400 Hz、中心 4000 Hz，所以采样率至少 2×4400 = 8800 Hz；取 fs = 12000 Hz，Nyquist 6000 Hz 留了余量，避免带内信号在采集中被折叠。
- 频率分辨率 Δf 想压到 1 Hz，采样时长 T ≥ 1/Δf = 1 s；实测 0.5 s（6000 点）已能把 591 Hz 峰值定位到整数赫兹内，只是相邻阶次若靠得太近时需要拉到 1 s 才分得开。
- 检查全速域：BPFO 超过带上限 4400 Hz 要到 n ≈ 4066 min⁻¹，而样本 nG 才 546 min⁻¹，所以这台轴承在额定转速内 BPFO 一直落在带内，包络解调全程有效。
结论很直接：fs=12000 Hz、T≥1 s、带通 3600–4400 Hz，这套参数对 KVE 08 0350 从 20 到 546 min⁻¹ 都吃得住，不用随转速换档。
五、用样本自带算例交叉校验
KMF 样本 P68 给了一个 KVE 10 0300 的寿命算例：DL=300、DW=10、Fa=2.50 kN、Fr=1.00 kN、Mk=0.38 kNm、n=20 min⁻¹，算得 P=5.75 kN、Lh=77000 h；P67 的静强度算例给 C0a=99 kN、Foa=7.5 kN，得到静安全系数 S0=2.36。我们把同一型号的几何代进特征频率公式：DL=300、DW=10、α=45°、Z≈92，n=20 min⁻¹ 时 BPFO ≈ 15.0 Hz。这一档比 KVE 08 0350 的 21.6 Hz 还低，因为 KVE 10 的节圆直径 300 更小、球数约 92 也更少；换型号时只要把 D_L、D_W、Z 三数换掉，同一段代码直接重算，不必改结构。
同一个型号，样本用 DIN ISO 281 的寿命指数 p=3 算出 77000 小时，我们用振动公式算出 15 Hz 的故障特征——两条线对上了：低速下寿命还很长，但包络谱已经能给出外圈特征频率，说明监测要趁早，不能等寿命快到才查。静安全系数 S0=2.36 落在样本图表认可的可用区间内，校核结论是这台轴承在给定载荷下静态不超限、动态可监测，可以放心纳入在线诊断清单。
〔来源：德国KMF钢丝轨道轴承 P67/P68〕
数据来源：德国KMF钢丝轨道轴承 P6/P12/P14/P67/P68；ISO 13373-1（振动状态监测）；DIN ISO 281
参考链接：机器人关节薄壁轴承
AI 声明：本文由 AI 协助整理，数据经人工核对，仅供参考
