一种基于振动谱分析的电主轴轴承结构参数推断方法与流程

1.本发明属于轴承结构参数推断的技术领域，具体涉及一种基于振动谱分析的电主轴轴承结构参数推断方法。


背景技术：

2.机床电主轴通常工作在高转速、低负载的场景下，具有很高的加工和装配精度。轴承作为连接电主轴转子和定子的关键部件，同时又是电主轴系统的薄弱刚性环节，在机床长期的生产服役过程中可能出现由于润滑不良、动态磨损、冲击碰撞等引起的轴承故障，对生产安全和零件质量造成威胁。轴承故障信息通过声发射、热辐射、振动等形式与环境形成耦合反馈，现代故障诊断技术通过收集反馈信息，采用先进的信号处理手段剥离故障信息，可以实现轴承的状态评估和故障预警。其中振动信号是最常用的形式，具有感受直观、易于采集等特点，在制造行业中得到大量应用。
3.但故障诊断技术普遍需要轴承的结构参数作为支撑条件，而由于生产厂商技术保密和电主轴模块化封装等原因，很难直接获取轴承结构参数，故障诊断技术的应用受到极大限制。得益于电主轴超高的制造和装配精度，即使在高速运行状态下其振动信号频谱与理论仍表现出高度的一致性，可以通过谱分析手段对轴承结构参数进行可靠推断。
4.接触角是轴承参数推断的重要基础之一，电主轴与机械主轴相似，工作过程中同时承受径向和轴向负荷，但由于其高速的应用场景，轴向和径向负载均较机械主轴大幅降低，且其结构不涉及复杂的机械传动，因此广泛选配接触角为15
°
或25
°
的角接触球轴承。此外，轴承高度的iso标准化，同样为参数推断提供了便利。


技术实现要素：

5.本发明的目的在于提供一种基于振动谱分析的电主轴轴承结构参数推断方法，旨在解决上述问题。
6.本发明主要通过以下技术方案实现：
7.一种基于振动谱分析的电主轴轴承结构参数推断方法，包括以下步骤：
8.步骤s100：采用数字滤波器对电主轴的振动信号x(t)进行带通滤波，消除低频和高频分量对频谱分析过程的干扰；
9.步骤s200：读取并验证电主轴转频f
r
，在指定的频谱区间进行最大值搜索，确定基频f
b
，并通过其谐波频率信息进行修正；
10.步骤s300：电主轴轴承接触角为15
°
或25
°
，以此为基础构建比例因子λ的候选集θ＝{θ
15
°
,θ
25
°
}，通过计算和逻辑判定确定比例因子λ；
11.步骤s400：定义调制系数ρ，指定中心调制频率f
c
和载波频率f
e
，在频谱上滑动计算形成调制系数序列p
(κ)
,κ＝1,2，绘制p
(κ)
茎叶图构建索引序号构成集合i，j，则z＝i∩j；
12.步骤s500：粗估轴径d
est
，以比例因子λ和滚子数量z为约束条件查阅轴承手册，确定轴承外径d和内径d。
13.为了更好地实现本发明，进一步地，在步骤s100中，通过传感器采集电主轴空运行、稳态的振动信号x(t)，获取转频f
r
；传感器的采集位置为电主轴前端轴承外壳。
14.为了更好地实现本发明，进一步地，采用数字滤波器对振动信号x(t)进行带通滤波的滤波通带为[10,20f
r
]。
[0015]
为了更好地实现本发明，进一步地，在步骤s200中，fft变换获取稳态振动信号x(t)的幅值谱x(f)，并根据谐波频率信息验证转频f
r
；在[0.4f
r
,0.47f
r
]区间内搜索局部最大值确定基频f
b
。
[0016]
为了更好地实现本发明，进一步地，在步骤s200中，利用谐波频率信息对f
b
进行匹配修正：计算nf
b
,n＝1,2,
…
,10，检查x(f)序列是否存在局部极值落在nf
b
的小邻域范围内([nf
b
‑
nf
δ
,nf
b
nf
δ
]，f
δ
为频谱分辨率，若超出范围，以该极值频率f
m
对f
b
进行修正f
b
＝f
m
/n。
[0017]
为了更好地实现本发明，进一步地，在步骤s300中，计算主轴轴承比例因子λ：
[0018]
构建比例因子λ候选集θ＝{θ
15
°
,θ
25
°
}，其中θ
τ
通过下式计算获得：
[0019][0020]
若|round(θ
15
°
)
‑
θ
15
°
|≤|round(θ
25
°
)
‑
θ
25
°
|，则比例因子λ＝round(θ
15
°
)，反之λ＝round(θ
25
°
)。
[0021]
为了更好地实现本发明，进一步地，在步骤s400中，计算主轴轴承滚子数量z：
[0022]
步骤s401：定义调制系数ρ：
[0023]
ρ＝sig<δ(f)
·
x(f)>
ꢀꢀꢀ
(2)
[0024]
其中，
[0025]
为有界非线性变换函数；
[0026][0027]
其中：f
c
为中心调制频率；f
e
为载波频率；
[0028]
步骤s402：分别取(f
c
,f
e
)＝(if
b
,f
r
)和(jf
r
‑
jf
b
,f
r
)，依次取i＝1,2,
…
,50和j＝1,2,
…
,35在x(f)上滑动计算调制系数序列
[0029]
步骤s403：绘制p
(1)
和p
(2)
序列的茎叶图并识别局部极值，以极值对应的索引序号构成集合i，j，则z＝i∩j；
[0030]
步骤s404：由几何关系推导轴承比例因子λ与滚子数上限z
max
之间的明确数量关系为z
max
＝floor(λπ)，计算并验证z<z
max
。
[0031]
本发明的有益效果：
[0032]
本发明基于电主轴振动信号频谱与理论的一致性开展振动谱分析，实现电主轴轴承结构参数的逆向推断。本发明依次进行信号采集、滤波预处理、fft变换、基频搜索、比例因子计算、调制系数序列构建等操作，最终获得的电主轴轴承结构参数包括：比例因子、滚子个数、轴承内外径尺寸，本发明方法的实用性和有效性在某航空结构件生产加工中心电主轴的实施例中得到验证，具有较好的实用性。
附图说明
[0033]
图1为本发明的流程图；
[0034]
图2为轴承比例因子λ的示意图；
[0035]
图3为电主轴时域振动信号及频域幅值谱；
[0036]
图4为p
(1)
序列茎叶图；
[0037]
图5为p
(2)
序列茎叶图。
具体实施方式
[0038]
实施例1：
[0039]
一种基于振动谱分析的电主轴轴承结构参数推断方法，包括以下步骤：
[0040]
步骤s100：采用数字滤波器对电主轴的振动信号x(t)进行带通滤波，消除低频和高频分量对频谱分析过程的干扰；
[0041]
步骤s200：读取并验证电主轴转频f
r
，在指定的频谱区间进行最大值搜索，确定基频f
b
，并通过其谐波频率信息进行修正；
[0042]
步骤s300：电主轴轴承接触角为15
°
或25
°
，以此为基础构建比例因子λ的候选集θ＝{θ
15
°
,θ
25
°
}，通过计算和逻辑判定确定比例因子λ；
[0043]
步骤s400：定义调制系数ρ，指定中心调制频率f
c
和载波频率f
e
，在频谱上滑动计算形成调制系数序列p
(κ)
,κ＝1,2，绘制p
(κ)
茎叶图构建索引序号构成集合i，j，则z＝i∩j；
[0044]
步骤s500：粗估轴径d
est
，以比例因子λ和滚子数量z为约束条件查阅轴承手册，确定轴承外径d和内径d。
[0045]
本发明基于电主轴振动信号频谱与理论的一致性开展振动谱分析，实现电主轴轴承结构参数的逆向推断。本发明依次进行信号采集、滤波预处理、fft变换、基频搜索、比例因子计算、调制系数序列构建等操作，最终获得的电主轴轴承结构参数包括：比例因子、滚子个数、轴承内外径尺寸，本发明方法的实用性和有效性在某航空结构件生产加工中心电主轴的实施例中得到验证，具有较好的实用性。
[0046]
实施例2：
[0047]
一种基于振动谱分析的电主轴轴承结构参数推断方法，包括以下步骤：
[0048]
1.布置传感器采集电主轴空运行稳态振动信号x(t)，获取转频f
r
；
[0049]
2.设计数字滤波器，对信号x(t)进行带通滤波，滤波通带为[10,20f
r
]；
[0050]
3.fft变换获取x(t)的幅值谱x(f)，并根据谐波频率信息验证转频f
r
；
[0051]
4.在[0.4f
r
,0.47f
r
]区间内搜索局部最大值确定基频f
b
；
[0052]
5.利用谐波频率信息对f
b
进行匹配修正：计算nf
b
,n＝1,2,
…
,10，检查x(f)序列是否存在局部极值落在nf
b
的小邻域范围内([nf
b
‑
nf
δ
,nf
b
nf
δ
]，f
δ
为频谱分辨率)，若超出范围，以该极值频率f
m
对f
b
进行修正f
b
＝f
m
/n；
[0053]
6.计算主轴轴承比例因子λ，其物理意义如图2所示；
[0054]
6.1.构建比例因子λ候选集θ＝{θ
15
°
,θ
25
°
}，其中θ
τ
通过下式计算获得：
[0055][0056]
6.2.若|round(θ
15
°
)
‑
θ
15
°
|≤|round(θ
25
°
)
‑
θ
25
°
|，则比例因子λ＝round(θ
15
°
)，反之
λ＝round(θ
25
°
)；
[0057]
7.计算主轴轴承滚子数量z；
[0058]
7.1.定义调制系数ρ：
[0059]
ρ＝sig<δ(f)
·
x(f)>
ꢀꢀꢀ
(2)
[0060]
式中，
[0061]
为有界非线性变换函数；
[0062]
f
c
为中心调制频率；f
e
为载波频率。
[0063]
7.2.分别取(f
c
,f
e
)＝(if
b
,f
r
)和(jf
r
‑
jf
b
,f
r
)，依次取i＝1,2,
…
,50和j＝1,2,
…
,35在x(f)上滑动计算调制系数序列
[0064]
7.3.绘制p
(1)
和p
(2)
序列的茎叶图并识别局部极值，以极值对应的索引序号构成集合i，j，则z＝i∩j；
[0065]
7.4.由几何关系推导轴承比例因子λ与滚子数上限z
max
之间的明确数量关系为z
max
＝floor(λπ)，计算并验证z<z
max
；
[0066]
8.根据主轴外形尺寸粗估轴径d
est
，以比例因子λ和滚子数量z为约束条件查阅轴承手册，确定轴承外径d和内径d。
[0067]
本发明基于电主轴振动信号频谱与理论的一致性开展振动谱分析，实现电主轴轴承结构参数的逆向推断。本发明依次进行信号采集、滤波预处理、fft变换、基频搜索、比例因子计算、调制系数序列构建等操作，最终获得的电主轴轴承结构参数包括：比例因子、滚子个数、轴承内外径尺寸，本发明方法的实用性和有效性在某航空结构件生产加工中心电主轴的实施例中得到验证，具有较好的实用性。
[0068]
实施例3：
[0069]
一种基于振动谱分析的电主轴轴承结构参数推断方法，某航空结构件生产加工中心呈卧式结构，其主轴为整体式电主轴，型号为：mfws
‑
2307/24。在有限的参数条件下，开展本发明所述的基于振动谱分析的轴承结构参数推断，可以可靠地获得该型号电主轴的轴承结构参数，为后续开展轴承状态评估和故障预警奠定基础。布置加速度传感器采集主轴空载、平稳运行状态下的振动信号如图3(a)所示，采集位置为电主轴前端轴承外壳，采样频率为25600hz，主轴转速为15000rpm，即250hz，时域信号rms为6.857m/s
‑2，滤波并做fft变换得到幅值谱如图3(b)。
[0070]
图3(b)中可以精准识出250hz、500hz、750hz、1000hz、1500hz谱线，与加工中心面板设置的主轴转速信息一致，因此判定f
r
＝250hz。f
b
搜索区间为[100,117.3]hz，区间内113.05hz处取得最大值，验证113.05hz的前10阶谐波处均有局部极值落在邻域范围内，即可判断f
b
＝113.05hz。比例因子λ候选集θ计算结果为{10.10,9.48}，根据步骤6.2.提供的准则容易判断λ＝10。
[0071]
选取(f
c
,f
e
)＝(i
×
113.05,250)，绘制p
(1)
如图4，由图可以构建索引集合i＝{25,42,45,50}。选取(f
c
,f
e
)＝(j
×
139.95,250)，绘制p
(2)
如图5，得索引集合j＝{25,28}，进而z＝i∩j＝25。z
max
＝floor(10π)＝31，z<z
max
。
[0072]
根据图纸提供的电主轴外形尺寸d
est
＝68mm，在λ＝10和z＝25的约束条件下查阅轴承参数手册，确定轴承外径d＝100mm和内径d＝70mm，skf品牌的角接触轴承参考型号为71914cd/p4a(混合陶瓷球轴承)。
[0073]
本发明基于电主轴振动信号频谱与理论的一致性开展振动谱分析，实现电主轴轴承结构参数的逆向推断。本发明依次进行信号采集、滤波预处理、fft变换、基频搜索、比例因子计算、调制系数序列构建等操作，最终获得的电主轴轴承结构参数包括：比例因子、滚子个数、轴承内外径尺寸，本发明方法的实用性和有效性在某航空结构件生产加工中心电主轴的实施例中得到验证，具有较好的实用性。
[0074]
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。