笔者之前曾从事过风电机组的轴承故障诊断工作，对于其中轴承故障频率的计算总是一知半解，最近又重新翻阅了一些资料进行研究，以防忘记，通过此博客将其中的重点进行记录总结。
轴承结构与故障频率
常见的滚动轴承主要包括下面五个组成部分：
- 内圈：通常与轴配合
- 外圈：通常支撑滚动体
- 滚动体：在内圈与外圈间承受载荷
- 保持架：分离滚动体以减少摩擦
- 密封圈：延长使用寿命，保存润滑剂，防止异物侵入
滚动轴承可能由于润滑不良、载荷过大、材质不当、轴承内落入异物、锈蚀等原因，引起轴承工作表面上的剥落、裂纹、压痕、腐蚀凹坑等缺陷或局部损伤。当滚动轴承另一工作表面通过某个缺陷点时，就会产生一个微弱的冲击脉冲信号。随着转轴的旋转，工作表面不断与缺陷点接触冲击，从而产生一个周期性的冲击振动信号。缺陷点位于不同的元件工作表面，冲击振动信号的周期间隔即频率是不相同的，这个频率就称为滚动轴承的故障特征频率。一般来讲，故障特征频率包括：外圈损伤频率 f b p f o f_{bpfo} fbpfo、内圈损伤频率 f b p f i f_{bpfi} fbpfi、滚珠损伤频率 f b s f f_{bsf} fbsf以及保持架损伤频率 f f t f f_{ftf} fftf。
轴承基本参数
我们先看看上面图中的字母标记所代表的参数含义：
- v i v_i vi：内圈圆周速度
- v m v_m vm：滚动体的公转圆周速度
- v e v_e ve：外圈圆周速度
- w i w_i wi：内圈角速度
- w m w_m wm：滚动体的公转角速度
- w e w_e we：外圈角速度
- D i D_{i} Di：内圈直径
- D m D_m Dm：滚动体中心圆直径
- D e D_e De：外圈直径
- d：滚动体直径
- α \alpha α：接触角
这些参数满足下列关系：
- D i = D m − d c o s ( α ) D_{i}=D_{m}-dcos(\alpha) Di=Dm−dcos(α)
- D e = D m + d c o s ( α ) D_{e}=D_{m}+dcos(\alpha) De=Dm+dcos(α)
- v m = v i + v e 2 v_{m}=\frac{v_{i}+v_{e}}{2} vm=2vi+ve
- w = 2 π f w=2 \pi f w=2πf
- v = w r v=wr v=wr
故障特征频率计算
这里先介绍一些频率的意义：
- f m f_m fm：保持架旋转频率
- f i f_i fi：内圈旋转频率
- f e f_e fe：外圈旋转频率
保持架旋转频率 f m f_m fm
由于
v
m
=
v
i
+
v
e
2
v_m=\frac{v_i+v_e}{2}
vm=2vi+ve，我们可以得到保持架的公转频率
f
m
f_m
fm：
w
m
D
m
2
=
w
i
D
i
2
+
w
e
D
e
2
2
π
f
m
D
m
=
π
f
i
D
i
+
π
f
e
D
e
2
f
m
=
f
i
2
(
1
−
d
c
o
s
(
α
)
D
m
)
+
f
e
2
(
1
+
d
c
o
s
(
α
)
D
m
)
\begin{aligned} w_m\frac{D_m}{2}&=\frac{w_i\frac{D_i}{2}+w_e\frac{D_e}{2}}{2} \\ \pi f_m D_m&=\frac{\pi f_i D_i+ \pi f_e D_e}{2} \\ f_m= \frac{f_i}{2}&(1-\frac{dcos(\alpha)}{D_m})+\frac{f_e}{2}(1+\frac{dcos(\alpha)}{D_m}) \end{aligned}
wm2DmπfmDmfm=2fi=2wi2Di+we2De=2πfiDi+πfeDe(1−Dmdcos(α))+2fe(1+Dmdcos(α))
通常我们也将该频率称为FTF。
外圈损伤特征频率 f b p f o f_{bpfo} fbpfo
实质上，计算
f
b
p
f
o
f_{bpfo}
fbpfo就是在计算滚动体相对于外圈的转动频率，我们可以先计算出保持架相对于外圈的转动频率
f
r
o
f_{ro}
fro，由于保持架上有n个滚动体，所以
f
b
p
f
o
=
n
f
r
o
f_{bpfo}=nf_{ro}
fbpfo=nfro
f
r
o
=
f
m
−
f
e
f
r
o
=
f
i
−
f
e
2
(
1
−
d
c
o
s
(
α
)
D
m
)
\begin{aligned} f_{ro}&=f_{m}-f_{e} \\ f_{ro}&=\frac{f_i-f_e}{2}(1-\frac{dcos(\alpha)}{D_m}) \end{aligned}
frofro=fm−fe=2fi−fe(1−Dmdcos(α))
内圈损伤特征频率 f b p f i f_{bpfi} fbpfi
计算
f
b
p
f
i
f_{bpfi}
fbpfi的过程和计算
f
b
p
f
o
f_{bpfo}
fbpfo的过程是类似的，我们可以先计算出保持架相对于内圈的转动频率
f
r
i
f_{ri}
fri，自然可得
f
b
p
f
i
=
n
f
r
i
f_{bpfi}=nf_{ri}
fbpfi=nfri
f
r
i
=
f
m
−
f
i
f
r
i
=
f
e
−
f
i
2
(
1
+
d
c
o
s
(
α
)
D
m
)
\begin{aligned} f_{ri}&=f_m-f_i \\ f_{ri}&=\frac{f_e-f_i}{2}(1+\frac{dcos(\alpha)}{D_m}) \end{aligned}
frifri=fm−fi=2fe−fi(1+Dmdcos(α))
滚动体损伤特征频率 f b s f f_{bsf} fbsf
计算
f
b
s
f
f_{bsf}
fbsf时，实质上就是计算滚动体的自转相对于内圈或者外圈的转动频率。以内圈为例，由于滚动体与内圈的接触点线速度是一致的，设自转角速度为
w
b
s
f
w_{bsf}
wbsf，我们便可得到下式
w
b
s
f
d
=
(
w
m
−
w
i
)
D
i
f
b
s
f
=
D
i
d
(
f
m
−
f
i
)
=
D
m
−
d
c
o
s
(
α
)
d
f
e
−
f
i
2
(
1
+
d
c
o
s
(
α
)
D
m
)
=
f
e
−
f
i
2
D
m
d
(
1
−
d
c
o
s
(
α
)
D
m
)
(
1
+
d
c
o
s
(
α
)
D
m
)
=
f
e
−
f
i
2
D
m
d
(
1
−
(
d
c
o
s
(
α
)
D
m
)
2
)
\begin{aligned} w_{bsf}d&=(w_m-w_i)D_i \\ f_{bsf}&=\frac{D_i}{d}(f_m-f_i) \\ &=\frac{D_{m}-dcos(\alpha)}{d}\frac{f_e-f_i}{2}(1+\frac{dcos(\alpha)}{D_m}) \\ &=\frac{f_e-f_i}{2}\frac{D_m}{d}(1-\frac{dcos(\alpha)}{D_m})(1+\frac{dcos(\alpha)}{D_m}) \\ &=\frac{f_e-f_i}{2}\frac{D_m}{d}(1-(\frac{dcos(\alpha)}{D_m})^2) \end{aligned}
wbsfdfbsf=(wm−wi)Di=dDi(fm−fi)=dDm−dcos(α)2fe−fi(1+Dmdcos(α))=2fe−fidDm(1−Dmdcos(α))(1+Dmdcos(α))=2fe−fidDm(1−(Dmdcos(α))2)
至此，所有的轴承故障特征频率均计算完毕。
参考
滚动轴承故障频率推导 - 知乎 (zhihu.com)
滚珠轴承故障频率的计算 - 知乎 (zhihu.com)
(67条消息) 滚动轴承特征频率计算推导过程_雪花儿帅帅的博客-CSDN博客
3万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
