滾動(dòng)軸承的振動(dòng)信號(hào)特征分析報(bào)告

1、南昌航空大學(xué)實(shí)驗(yàn)報(bào)告課程名稱(chēng)： 數(shù)字信號(hào)處理實(shí)驗(yàn)名稱(chēng)：滾動(dòng)軸承的振動(dòng)信號(hào)特征分析實(shí)驗(yàn)時(shí)間： 2013年5月14日班 級(jí) ： 100421學(xué) 號(hào) ： 10042134姓 名 ： 吳涌濤成 績(jī) ： 滾動(dòng)軸承的振動(dòng)信號(hào)特征分析一、 實(shí)驗(yàn)?zāi)康睦脭?shù)字信號(hào)處理課程中學(xué)習(xí)的序列運(yùn)算、周期信號(hào)知識(shí)、DFT知識(shí)，對(duì)給定的正常軸承數(shù)據(jù)、內(nèi)圈故障軸承數(shù)據(jù)、外圈故障軸承數(shù)據(jù)、滾珠故障軸承數(shù)據(jù)進(jìn)行時(shí)域特征或頻域特征提取和分析，找出能區(qū)分四種狀態(tài)（滾動(dòng)軸承的外圈故障、內(nèi)圈故障、滾珠故障和正常狀態(tài)）的特征。二、 實(shí)驗(yàn)原理振動(dòng)機(jī)理分析：機(jī)械在運(yùn)動(dòng)時(shí)，由于旋轉(zhuǎn)件的不平衡、負(fù)載的不均勻、結(jié)構(gòu)剛度的各向異性、間隙、潤(rùn)滑不良、支撐

2、松動(dòng)等因素，總是伴隨著各種振動(dòng)。振動(dòng)的幅值、頻率和相位是振動(dòng)的三個(gè)基本參數(shù)，稱(chēng)為振動(dòng)三要素。幅值：幅值是振動(dòng)強(qiáng)度的標(biāo)志，它可以用峰值、有效值、平均值等方法來(lái)表示。頻率：不同的頻率成分反映系統(tǒng)內(nèi)不同的振源。通過(guò)頻譜分析可以確定主要頻率成分及其幅值大小，從而尋找振源，采取相應(yīng)的措施。相位：振動(dòng)信號(hào)的相位信息十分重要，如利用相位關(guān)系確定共振點(diǎn)、測(cè)量振型、旋轉(zhuǎn)件動(dòng)平衡、有源振動(dòng)控制、降噪等。對(duì)于復(fù)雜振動(dòng)的波形分析，各諧波的相位關(guān)系是不可缺少的。 在振動(dòng)測(cè)量時(shí)，應(yīng)合理選擇測(cè)量參數(shù)，如振動(dòng)位移是研究強(qiáng)度和變形的重要依據(jù)；振動(dòng)加速度與作用力或載荷成正比，是研究動(dòng)力強(qiáng)度和疲勞的重要依據(jù)；振動(dòng)速度決定了噪聲的高

3、低，人對(duì)機(jī)械振動(dòng)的敏感程度在很大頻率范圍內(nèi)是由速度決定的。速度又與能量和功率有關(guān)，并決定動(dòng)量的大小。提取振動(dòng)信號(hào)的幅域、時(shí)域、頻域、時(shí)頻域特征，根據(jù)特征進(jìn)行故障有無(wú)、故障類(lèi)型和故障程度三個(gè)層次的判斷。三、 實(shí)驗(yàn)內(nèi)容Step1、使用importdata（）函數(shù)導(dǎo)入振動(dòng)數(shù)據(jù)。Step2、把大量數(shù)據(jù)分割成周期為單元的數(shù)據(jù)，分割方法為：設(shè)振動(dòng)信號(hào)為xk(k=1,2,3,n)采樣頻率為fs，傳動(dòng)軸的轉(zhuǎn)動(dòng)速率為Vr。采樣間隔為： （1）旋轉(zhuǎn)頻率為： （2）傳動(dòng)軸的轉(zhuǎn)動(dòng)周期為： （3）由式（1）和（3）可推出振動(dòng)信號(hào)一個(gè)周期內(nèi)采樣點(diǎn)數(shù)N： （4）由式（2）可得到傳動(dòng)軸的轉(zhuǎn)動(dòng)基頻fr=29.95Hz，再由式（

4、3）可得到一個(gè)周期內(nèi)采樣點(diǎn)數(shù)N=400.67，取N=400。Step3、提取振動(dòng)信號(hào)的特征，分析方法包括：1、時(shí)域統(tǒng)計(jì)分析指標(biāo)（波形指標(biāo)(Shape Factor)、峰值指標(biāo)(Crest Factor)、脈沖指標(biāo)(Impulse Factor)、裕度指標(biāo)(Clearance Factor)、峭度指標(biāo)(KurtosisValue) ）等，相關(guān)計(jì)算公式如下：（1） 波形指標(biāo)： (5)其中，為峰值，為均值。計(jì)算公式如下：（2） 波峰指標(biāo)： (6)其中，為均方根值。均方根值計(jì)算如下： (7)（3） 峭度指標(biāo)： (8)其中，為峭度，計(jì)算式如下： (9)（4） 脈沖指標(biāo)： (10)（5） 偏度指標(biāo)： (1

5、1)其中，為偏度，計(jì)算式如下： (12)（6） 裕度指標(biāo)： (13)其中，為方根幅值，計(jì)算公式如下： 2、頻域分析，包括離散頻譜序列統(tǒng)計(jì)分析指標(biāo)，階次分析等。滾動(dòng)軸承的理論特征頻率計(jì)算公式如下：（1） 內(nèi)圈故障頻率： =161.73 Hz（2） 外圈故障頻率： =107.82 Hz（3） 滾珠故障頻率： =141.17 Hz其中，N 為滾珠數(shù)，d為滾珠直徑，a 為接觸角,D為軸承節(jié)徑。Step4（提高要求）、提取未知信號(hào)的特征，與已知信號(hào)的特征進(jìn)行比對(duì)分類(lèi)，診斷信號(hào)是否存在故障以及故障類(lèi)型。未知故障類(lèi)型信號(hào)：1.mat四、實(shí)驗(yàn)的程序代碼（1）、待測(cè)信號(hào)1function y1,y2,y3,y4

6、=KuttosisFactor();s=importdata(E:原始信號(hào)內(nèi)圈故障105.mat);s_value=s.X105_DE_time;s_cutvalue=s_value(1:24000);tree=wpdec(s_cutvalue,3,sym6);det=wpcoef(tree,2);sigma=median(abs(det)/0.675;alpha=2;thr=wpbmpen(tree,sigma,alpha);keepapp=1;xd=wpdencmp(tree,s,nobest,thr,keepapp);T=400;Vaule1=zeros(1,50);for n=1:1:

7、50 B=sum(xd(1+(n-1)*T):(n*T).4)/T; Rms=sqrt(sum(xd(1+(n-1)*T):(n*T).2)/T); Vaule1(n)=B/(Rms4);end; s=importdata(E:原始信號(hào)正常97.mat);s_value1=s.X097_DE_time;s_cutvalue=s_value1(1:24000);tree=wpdec(s_cutvalue,3,sym6);det=wpcoef(tree,2);sigma=median(abs(det)/0.675;alpha=2;thr=wpbmpen(tree,sigma,alpha);keep

8、app=1;xd=wpdencmp(tree,s,nobest,thr,keepapp);Vaule2=zeros(1,50);for n=1:1:50 B=sum(xd(1+(n-1)*T):(n*T).4)/T; Rms=sqrt(sum(xd(1+(n-1)*T):(n*T).2)/T); Vaule2(n)=B/(Rms4);end; s=importdata(E:原始信號(hào)外圈故障130.mat);s_value1=s.X130_DE_time;s_cutvalue=s_value1(1:24000);tree=wpdec(s_cutvalue,3,sym6);det=wpcoef(t

9、ree,2);sigma=median(abs(det)/0.675;alpha=2;thr=wpbmpen(tree,sigma,alpha);keepapp=1;xd=wpdencmp(tree,s,nobest,thr,keepapp);Vaule3=zeros(1,50);for n=1:1:50 B=sum(xd(1+(n-1)*T):(n*T).4)/T; Rms=sqrt(sum(xd(1+(n-1)*T):(n*T).2)/T); Vaule3(n)=B/(Rms4);end; s=importdata(E:原始信號(hào)滾珠故障118.mat);s_value1=s.X118_DE

10、_time;s_cutvalue=s_value1(1:24000);tree=wpdec(s_cutvalue,3,sym6);det=wpcoef(tree,2);sigma=median(abs(det)/0.675;alpha=2;thr=wpbmpen(tree,sigma,alpha);keepapp=1;xd=wpdencmp(tree,s,nobest,thr,keepapp);Vaule4=zeros(1,50);for n=1:1:50 B=sum(xd(1+(n-1)*T):(n*T).4)/T; Rms=sqrt(sum(xd(1+(n-1)*T):(n*T).2)/T

11、); Vaule4(n)=B/(Rms4);end; s=importdata(E:待診斷信號(hào)1.mat);s_value1=s;s_cutvalue=s_value1(1:24000);tree=wpdec(s_cutvalue,3,sym6);det=wpcoef(tree,2);sigma=median(abs(det)/0.675;alpha=2;thr=wpbmpen(tree,sigma,alpha);keepapp=1;xd=wpdencmp(tree,s,nobest,thr,keepapp);T=400;Vaule5=zeros(1,50);for n=1:1:50 B=su

12、m(xd(1+(n-1)*T):(n*T).4)/T; Rms=sqrt(sum(xd(1+(n-1)*T):(n*T).2)/T); Vaule5(n)=B/(Rms4);end; y1(1)=(Vaule1(1)+Vaule1(2)+Vaule1(3)/3;y2(1)=(Vaule2(1)+Vaule2(2)+Vaule2(3)/3;y3(1)=(Vaule3(1)+Vaule3(2)+Vaule3(3)/3;y4(1)=(Vaule4(1)+Vaule4(2)+Vaule4(3)/3;y5(1)=(Vaule5(1)+Vaule5(2)+Vaule5(3)/3;y1(50)=(Vaule

13、1(50)+Vaule1(49)+Vaule1(48)/3;y2(50)=(Vaule2(50)+Vaule2(49)+Vaule2(48)/3;y3(50)=(Vaule3(50)+Vaule3(49)+Vaule3(48)/3;y4(50)=(Vaule4(50)+Vaule4(49)+Vaule4(48)/3;y5(50)=(Vaule5(50)+Vaule5(49)+Vaule5(48)/3;for n=2:1:49 y1(n)=(Vaule1(n-1)+Vaule1(n)+Vaule1(n+1)/3; y2(n)=(Vaule2(n-1)+Vaule2(n)+Vaule2(n+1)/3; y3(n)=(Vaule3(n-1)+Vaule3(n)+Vaule3(n+1)/3; y4(n)=(Vaule4(n-1)+Vaule4(n)+Vaule4(n+1)/3; y5(n)=(Vaule5(n-1)+Vaule5(n)+Vaule5(n+1)/3; end;%繪圖k=linspace(1,50,50);plot(k,y1(k),-.k,k,y2(k),-r,k,y3(k),-pb,k,y4(k),-*g,k,y5(k),m:d);legend(內(nèi)圈故障,正常外圈故障,故障,診斷信號(hào));小結(jié)：