齿轮故障诊断的实验数据集及python处理
本文首發(fā)于“52phm”公眾號,未經(jīng)授權(quán),嚴禁轉(zhuǎn)載!
1.摘要
該試驗臺在不同小齒輪條件下進行測試,并通過加速度計進行振動信號采集,加速度計采樣率為10KHz、采樣時長為10s,采樣數(shù)據(jù)共3包,每一包數(shù)據(jù)對應(yīng)著不同故障類型,分別是健康狀態(tài)、齒輪斷齒、齒輪磨損狀態(tài)下的數(shù)據(jù)集。該數(shù)據(jù)集被授權(quán)于用于任何學(xué)術(shù)和研究目的。
2.試驗臺結(jié)構(gòu)
試驗臺裝置整體傳動系統(tǒng)主要由電機作為驅(qū)動輸入,通過皮帶帶動齒輪箱,而齒輪箱輸出端通過皮帶帶動剎車系統(tǒng)。具體試驗臺裝置、原理圖及加速度計安裝示意圖如下圖所示。
試驗臺設(shè)置
測試原理圖
加速度計的方向及位置
3.數(shù)據(jù)集解讀
數(shù)據(jù)集有3包數(shù)據(jù),每包數(shù)據(jù)采樣率為10kHz,采樣時長為10s;
電機轉(zhuǎn)速1420rpm/min,小齒輪15齒數(shù)、大齒輪110齒數(shù);
嚙合頻率=(1420/60)*15=355Hz,在進行頻譜分析中來看,實際嚙合頻率是365Hz左右;
數(shù)據(jù)與故障對應(yīng)關(guān)系
-
斷齒?
Gearbox_a_chipped_tooth_full_load_03_December_2009_10kHz_pos1
-
正常 ?
Gearbox_no_fault_full_load_01_December_2009_10kHz_pos1
-
3個齒磨損 ?
Gearbox_three_worn_teeth_full_load_13_December_2009_10kHz_pos1
齒輪磨損和脫落
4.振動分析
這里使用python寫了一個數(shù)據(jù)處理的程序,將mat數(shù)據(jù)保存到csv文件、以及給出振動原始波形、頻譜圖的計算方法。簡單提供一個baseline,方便大家使用數(shù)據(jù)集和做一些分析學(xué)習(xí)。
數(shù)據(jù)處理
""" @日期:2022-08-06 @作者:xiaozhi @網(wǎng)站:http://www.52phm.cn @數(shù)據(jù)來源:https://www.researchgate.net/publication/303792317_EXPERIMENTAL_DATASET_FOR_GEAR_FAULT_DIAGNOSIS """ import pandas as pd from scipy import io ? ? if __name__ == '__main__':df = pd.DataFrame()data1 = io.loadmat("./Gearbox_a_chipped_tooth_full_load_03_December_2009_10kHz_pos1.mat")data1_acc = data1['acc'].flatten() * 9.8605df['Gearbox_a_chipped_tooth_full_load_03_December_2009_10kHz_pos1'] = data1_acc ?data1 = io.loadmat("./Gearbox_three_worn_teeth_full_load_13_December_2009_10kHz_pos1.mat")data1_acc = data1['acc'].flatten() * 9.8605df['Gearbox_three_worn_teeth_full_load_13_December_2009_10kHz_pos1'] = data1_acc ?data1 = io.loadmat("./Gearbox_no_fault_full_load_01_December_2009_10kHz_pos1.mat")data1_acc = data1['acc'].flatten() * 9.8605df['Gearbox_no_fault_full_load_01_December_2009_10kHz_pos1'] = data1_acc ?df.to_csv("齒輪箱數(shù)據(jù)集轉(zhuǎn)速1420主動軸15齒從動輪110齒嚙合頻率比實際355Hz有偏差在365Hz左右.csv", index=False) ? ? ?振動分析
這里涉及到時域波形和頻譜分析,本次baseline提供快速傅里葉變換的頻譜分析,具體實現(xiàn)如下代碼所示。
""" @日期:2022-08-06 @作者:xiaozhi @網(wǎng)站:http://www.52phm.cn @數(shù)據(jù)來源:https://www.researchgate.net/publication/303792317_EXPERIMENTAL_DATASET_FOR_GEAR_FAULT_DIAGNOSIS """ import numpy as np import pandas as pd import matplotlib.pyplot as plt plt.rcParams['font.family'] = 'SimHei' plt.rcParams['axes.unicode_minus'] = False ? ? def get_wav(data, fs, title):plt.clf()plt.title(title)t = np.arange(len(data)) * (1/fs)plt.plot(t, data)plt.xlabel("時間[s]")plt.ylabel("幅值[m/s^2]")plt.savefig(title + '_wav.png')# plt.show() ? ? def get_spec(data, fs, title):plt.clf()plt.title(title)f = np.fft.rfftfreq(len(data), 1/fs)amp = abs(np.fft.rfft(data)) * 2 / len(data)plt.plot(f, amp)plt.xlabel("頻率[Hz]")plt.ylabel("幅值[m/s^2]")plt.savefig(title + '_spec.png')plt.show() ? ? if __name__ == '__main__':df = pd.read_csv("齒輪箱數(shù)據(jù)集轉(zhuǎn)速1420主動軸15齒從動輪110齒嚙合頻率比實際355Hz有偏差在365Hz左右.csv")for col in df.columns.tolist():data1 = df[col].values# get_wav(data1, fs=10000, title=col)get_spec(data1, fs=10000, title=col) ? ? ?-
正常齒輪
?
-
齒輪斷齒
-
齒輪磨損
5.數(shù)據(jù)和代碼獲取
關(guān)注公眾號“52phm”,回復(fù)“齒輪箱”或“齒輪箱數(shù)據(jù)集”即可獲取本文數(shù)據(jù)集和代碼。
6.參考資料
https://www.researchgate.net/publication/303792317_EXPERIMENTAL_DATASET_FOR_GEAR_FAULT_DIAGNOSIS
52phm社區(qū),http://www.52phm.cn,與工業(yè)互聯(lián)網(wǎng)人一起成長,社區(qū)能夠提供學(xué)習(xí)到真實工業(yè)場景技術(shù)知識以及落地案例,其中涵蓋工業(yè)數(shù)據(jù)集、工業(yè)標準庫、機理模型、設(shè)備知識、機器學(xué)習(xí)、 深度學(xué)習(xí)、特征工程、振動分析、工業(yè)視覺、邊緣硬件及傳感器等技術(shù)知識!
?
?
關(guān)注公眾號“52phm”,回復(fù)“齒輪箱”或“齒輪箱數(shù)據(jù)集”即可獲取本文數(shù)據(jù)集和代碼。
總結(jié)
以上是生活随笔為你收集整理的齿轮故障诊断的实验数据集及python处理的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 百度云 x 中国联通 | 立标杆,中国联
- 下一篇: Python后端工程师学习路线——大学版