首先声明,对于振动分析,本人是菜鸟。所以本文中可能有错误,或者不够专业和严谨。
最简单的振动检测指标是所谓通频值(Overall vibration value)。它是采集信号的均方根:
实例: 产生一个正弦波,并计算rms。该值越大,反映振动越大。
import numpy as np
from scipy import signal
import matplotlib.pyplot as plt
n = 512
dt = 0.01
t = np.linspace(1, n, n)*dt-dt
y = np.sin(2*np.pi*t)
plt.plot(t, y)
plt.show()
rms = np.sqrt(np.mean(y**2))
print(rms/N)
结果为 0.7 .正弦波的rms 是0.7。
RMS 的计算公式
波峰因数(Crest Factor)
Crest Factor (peak-to-RMS ratio) = (peak value)/(RMS value).
通频值和波峰因素振动仪表的方法的好处是产生一个单一的值来反映机器的健康。
CF+ 严重性
0 to 5 好 Good
6 to 10 满意 Satisfactory
11 to 15 不满意 Unsatisfactory above
15 不能接受 Unaccepatable
波峰因数+(Crest Factor +)
当轴承急剧损坏时,RMS 的值加大,与峰值接近,Crest Factor 反而会下降。
CF =oc Peak+oc RMS+oc (peak value)/(RMS value).
oc 代表正比于。我觉得就是乘一个比例系数就可以了。
PeakVue 方法
波峰因数计算量比较小,但是有许多的故障无法检测出来。实际上,轴承的许多早期故障是金属和金属摩擦产生的。表现为微弱的高频信号。如果和低频振动信号混合在一起。高频信号难以辨别。
PeakVue 方法是艾默生公司的专利技术,它的实现算法主要为:
1 使用高通滤波器(highpass Filter) 过滤掉低频信号,截止频率大约为1000Hz
2 使用希尔伯特变换 计算出包络线
3 对包络线 做FFT 变换。
我的项目
振动数据IEPE采集器
使用STM32F429 +ads127l01 实现IEPE 压电振动传感器ADC 采集,通过UDP 将数据传送到PC机。PC机使用Python 做后续数据处理。
测试床
python 程序
import socket
import numpy as np
from scipy.fftpack import fft
from scipy import signal
import struct
from matplotlib import pyplot as plt
def bytesToFloat(h1,h2,h3,h4):
ba = bytearray()
ba.append(h1)
ba.append(h2)
ba.append(h3)
ba.append(h4)
return struct.unpack("!f",ba)[0]
HOST = '192.168.31.98'
PORT = 3800
BUFSIZ = 512
ADDR = (HOST,PORT)
udpSerSock = socket.socket(socket.AF_INET,socket.SOCK_DGRAM)
udpSerSock.bind(ADDR)
t = np.linspace(0, 127, 128, True)
y =[]
for i in range(128):
y.append(0)
fig, (ax1, ax2) = plt.subplots(2, 1, sharex=False)
plt.show()
while True:
data, addr = udpSerSock.recvfrom(BUFSIZ)
for i in range(128):
y[i]=bytesToFloat(data[i*4+3],data[i*4+2],data[i*4+1],data[i*4+0])
b,a = signal.butter(16, 0.128, 'highpass')
filtered = signal.filtfilt(b, a, y)*20
yh = np.abs(signal.hilbert(filtered))
yy= np.abs(fft(yh))
yy[0] =0
ax1.clear()
ax1.set_title("Peakvue")
ax1.plot(t,yy)
ax1.axis([0, 128, 0,40])
ax2.clear()
ax2.plot(t,yh)
#ax2.axis([0, 128, -20, 20])
ax2.set_title("highpass")
plt.pause(0.1)
过几天慢慢解释