汽車振動學(xué)實驗報告

1、汽車振動學(xué)實驗報告姓名:舒鵬學(xué)號:20110802班級:2011級交通運輸2班指導(dǎo)教師：蔣淑霞實驗主題:基于Matlab的汽車振動分析其中，固有頻率：其中，固有頻率：實驗一：單自由度阻尼振動分析：包括欠阻尼、臨界阻尼和過阻尼一、實驗內(nèi)容（1）將老師給的源程序讀懂（包含臨界阻尼和欠阻尼）。（2）參考上述程序，改變阻尼比=0丄0.2,1和畫出不同阻尼比下的響應(yīng)曲線共10條，并比較各響應(yīng)曲線特性。（3）同時考慮過阻尼曲線。二、實驗前期準備與內(nèi)容分析阻尼自由振動：前面講的自由振動系統(tǒng)沒有考慮阻力的影響，實際由于阻力的存在系統(tǒng)機械能不可能守恒，振動中這種阻力稱為阻尼，如摩擦阻尼、電磁阻尼、介質(zhì)阻尼和結(jié)構(gòu)

2、阻尼。對于實際系統(tǒng)中的阻尼很難確定，工程中最常見的一種阻尼力學(xué)模型是粘性阻尼，如在流體中運動、在潤滑表面的滑動。粘性阻尼力與相對速度成正比：卩廠Ic:粘性阻尼系數(shù)，或阻尼系數(shù)，單位：Ns/m建立平衡位置坐標系，受力分析：阻尼力、彈性力、重力.得到動力學(xué)方程：實驗代碼與結(jié)果第一種情況：零界阻尼ksai=1第一種情況：零界阻尼ksai=1%零界阻尼，ksai表示為阻尼比k=103;m=4;w0=sqrt(k/m);%固有頻率v0=10;%初始速度x0=5;%初始位移t=(0：300)/100;x1=exp(-w0*t).*(x0+v0*t)axes(handles.axes1)plot（t，x1）

3、%在圖形axes1中輸出圖像xlabel（時間（s）;%x坐標表示時間ylabel（位移（mm）;%y坐標表示位移第二種情況：欠阻尼ksai=0.2%阻尼比為0.2k=103;m=4;x0=5%初始位移v0=10%初始速度w0=sqrt(k/m);%固有頻率wd=sqrt(l-ksai“2)*wO;%阻尼固有頻率t=(0:300)/100;x2=exp(-ksai*w0*t).*(x0*cos(wd*t)+(v0+ksai*w0*x0)/wd)*sin(wd*t)y=x0*exp(-ksai*w0*t)%漸近線axes(handles.axes2)plot(t,x2,k,t,y,g-);在ax

4、es2輸出圖像holdon;ksai=0.4%阻尼比0.4k=103;m=4;x0=5%初始位移v0=10%初始速度w0=sqrt(k/m);%固有頻率wd=sqrt(l-ksai“2)*w0;%阻尼固有頻率t=(0:300)/100;x2=exp(-ksai*w0*t).*(x0*cos(wd*t)+v0*sin(wd*t)axes(handles.axes2)plot(t,x2,y)holdon;ksai=0.6%阻尼比0.6k=103;m=4;x0=5%初始位移v0=10%初始速度w0=sqrt(k/m);%固有頻率wd=sqrt(l-ksai“2)*w0;%阻尼固有頻率t=(0:300

5、)/100;x2=exp(-ksai*w0*t).*(x0*cos(wd*t)+v0*sin(wd*t)axes(handles.axes2)plot(t,x2,m)holdon;legend(欠阻尼0.2,零界線，欠阻尼0.4,欠阻尼0.6)%對曲線進行命名xlabel(時間(s);ylabel(位移(mm);第三種情況：過阻尼ksai=2%阻尼比為2k=103;m=4;x0=5%初始位移v0=10%初始速度w0=sqrt(k/m);%固有頻率w1=w0*sqrt(ksai“2-l);t=(0:300)/100shw1t=(exp(w1*t)-exp(-w1*t)/2chw1t=(exp(w

6、1*t)+exp(-w1*t)/2x3=exp(-ksai*w0*t).*(x0*chw1t+(v0+ksai*w0*x0)/w1)*shw1t)axes(handles.axes3);plot(t,x3)xlabel(時間(s);ylabel(位移(mm);在gui中輸出結(jié)果：5320.52310500.523620.523EEEEEE00-o1.5時間(a5320.52310500.523620.523EEEEEE00-o1.5時間(a1.5時間（51.5時間(s欠阻尼0丄霧界線矢阻尼04-欠阻尼0一6過阻尼霹界阻尼實驗二：單自由度簡諧振動幅頻相頻特性曲線一、實驗內(nèi)容1.觀察單自由度簡諧振

7、動的幅頻相頻特性曲線。2掌握放大因子B隨著頻率比變化的規(guī)律及曲線隨著阻尼比z變化的規(guī)律。二、實驗前期準備與內(nèi)容分析穩(wěn)態(tài)響應(yīng)特性以s為橫坐標，畫出幅頻特性曲線：&32&J75&32&J75三、實驗代碼與結(jié)果幅頻曲線：s=(0:300)/100;%頻率比ksai=0.1%阻尼比0.1bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1);%輸出圖像holdon;ksai=02%阻尼比0.2bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,k-)holdon;ksai=0.3%阻尼比0.3bl=l./sqrt(l-s.2).2+(2

8、*ksai*s).2)plot(s,b1,b)holdon;ksai=0.4%阻尼比0.4bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,y)holdon;ksai=0.5%阻尼比0.5bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,m)holdon;b1=(0:60)/10plot(1,b1)legend(ksai=0.1,ksai=0.2,ksai=0.3,ksai=0.4,ksai=0.5)xlabel(間歇力頻率與固有頻率之比)；ylabel(響應(yīng)振幅與靜位移之比)；相頻曲線：s=(0:500)/100;%頻率

9、比ksai=0.1;%阻尼比0.3A二atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,k);holdon;ksai=0.2;%阻尼比0.2A二atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,m);holdon;ksai=0.3;%阻尼比0.3A二atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;

10、A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,b);holdon;ksai=0.4;%阻尼比0.4A=atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,y);holdon;ksai=0.5%阻尼比0.5A二atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi);holdon;legend(阻尼比0.1,阻尼比0.2,阻

11、尼比0.3,阻尼比0.4,阻尼比0.5)xlabel(頻率比);ylabel(相位角);幅頻曲線圖蘭蘭N豫也鑑専曙黑也雪相頻曲線圖：實驗三基于MatLab的偏心質(zhì)量引起的強迫振動分析、實驗內(nèi)容基于MatLab的偏心質(zhì)量引起的強迫振動分析。（1）參考書上的內(nèi)容p6-p8，特別是公式1.2.16和1.2.17，求動態(tài)放大因子隨頻率比和阻尼比=的變化曲線。w2fkmn（2）求響應(yīng)相位角滯后于激振力的相位角隨頻率比和阻尼比=的變化曲線。w2jkmn二、實驗前期準備與內(nèi)容分析旋轉(zhuǎn)機械如電動機、水泵、汽車發(fā)動機等，由于偏心質(zhì)量而引起強迫振動是很普遍的現(xiàn)象。偏心質(zhì)量引起的強迫振動分析mew2me0基于Mat

12、Lab的偏心質(zhì)量引起的強迫振動分析四、實驗代碼與結(jié)果實驗二程序與實驗一程序相同，只是只是代表物理含義名稱不同s=(0:300)/100;%頻率比ksai=0.1%阻尼比0.1bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1);%輸出圖像holdon;ksai=02%阻尼比0.2bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,k-)holdon;ksai=0.3%阻尼比0.3bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,b)holdon;ksai=0.4%阻尼比0.4bl=l./sq

13、rt(l-s.2).2+(2*ksai*s).2)plot(s,b1,y)holdon;ksai=0.5%阻尼比0.5bl=l./sqrt(l-s.2).2+(2*ksai*s).2)plot(s,b1,m)holdon;b1=(0:60)/10plot(1,b1)legend(ksai=0.1,ksai=0.2,ksai=0.3,ksai=0.4,ksai=0.5)xlabel(頻率比)；ylabel(放大因子)；s=(0:500)/100;%頻率比ksai=0.1;%阻尼比0.3A=atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i

14、)=A(1,i)+pi;endendplot(s,A*180/pi,k);holdon;ksai=0.2;%阻尼比0.2A=atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,m);holdon;ksai=0.3;%阻尼比0.3A=atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,b);holdon;ksai=0.4;%阻尼比0.4A=atan(2

15、*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi,y);holdon;ksai=0.5%阻尼比0.5A二atan(2*ksai*s./(l-s2);fori=1:size(A,2);ifA(1,i)0;A(1,i)=A(1,i)+pi;endendplot(s,A*180/pi);holdon;4,阻尼比legend(阻尼比0.1,阻尼比0.2,阻尼比0.3,阻尼比0.4,阻尼比xlabel(頻率比);ylabel(相位角);幅頻曲線圖11.5頻率比ksai=0.1ksai=0.2

16、ksai=0.3ksai=0.4ksai=0.52.51.522.533.544.55頻率比四：基于Matlab的三自由度振動分析、實驗內(nèi)容參考源程序1motaiqiujie.m符號方程的求解，解特征方程得到三個振動固有頻率；源程序2chapter2_3freed.m（1）首先讀懂程序（2）系統(tǒng)M和K矩陣如下，且m=1000kg；k=10人6，求固有頻率和三個振型的實頻和虛頻曲線、幅頻和相頻曲線和Nyquist圖。理解主振動含義。m0013k-k廠00M0卜-*-2Jt-40000-1-340二、實驗前期準備與內(nèi)容分析Matlab求特征值和特征向量:三、實驗代碼與結(jié)果a=0.001;b=0.0

17、01;M=100,0,0;0,100,0;0,0,100;K1=2,-1,0;-1,2,-1;0,-1,2;K=1.0e+5*K1;C=a*M+b*K;%比例阻尼V,D=eig(K,M);%V,D=EIG(A,B)producesadiagonalmatrixDofgeneralized%eigenvaluesandafullmatrixVwhosecolumnsarethe%correspondingeigenvectorssothatA*V=B*V*D.forj=1:1:3omega(j)=sqrt(D(j,j);fork=1:1:3Y(k,j)=V(k,j)/V(3,j);endendo

18、megaY%頻率取值范圍和步長%動剛度矩陣%頻率取值范圍和步長%動剛度矩陣%頻響函數(shù)矩陣H=zeros(701,3,3);forj=1:1:701Z二K-W(j廠2*M+w(j)*C*i;HH=inv(Z);H(j,:,:)=HH;end;Hr=real(H);Hi=imag(H);Habs=abs(H);Hang=(180/pi)*angle(H);%下面分別是三個振型的實頻和虛頻曲線；figure(1)subplot(3,2,1),plot(w,Hr(:,1,1),gridontitle(實頻曲線l),xlabel(omega),ylabel(Hr)subplot(3,2,2),plot(

19、w,Hi(:,1,1),gridontitle(虛頻曲線1),xlabel(omega),ylabel(Hi)subplot(3,2,3),plot(w,Habs(:,1,1),gridontitle(幅頻曲線1),xlabel(omega),ylabel(|H|)subplot(3,2,4),plot(w,Hang(:,1,1),gridontitle(相頻曲線1),xlabel(omega),ylabel(psi)subplot(3,2,6);plot(Hr(:,1,1),Hi(:,1,1),gridontitle(乃奎斯特圓l),xlabel(Re(H),ylabel(Im(H)%畫跨點

20、的頻響figure(2)subplot(3,2,1);plot(w,Hr(:,2,1),gridontitle(實頻曲線2),xlabel(omega),ylabel(Hr)subplot(3,2,2);plot(w,Hi(:,2,1),gridontitle(虛頻曲線2),xlabel(omega),ylabel(Hi)subplot(3,2,3);plot(w,Habs(:,2,1),gridontitle(幅頻曲線2),xlabel(omega),ylabel(|H|)subplot(3,2,4);plot(w,Hang(:,2,1),gridontitle(相頻曲線2),xlabel(

21、omega),ylabel(psi)subplot(3,2,6);plot(Hr(:,2,1),Hi(:,2,1),gridontitle(乃奎斯特圓2),xlabel(Re(H),ylabel(Im(H)%畫跨點的頻響figure(3)subplot(3,2,1);plot(w,Hr(:,3,1),gridontitle(實頻曲線3),xlabel(omega),ylabel(Hr)subplot(3,2,2);plot(w,Hi(:,3,1),gridontitle(虛頻曲線3),xlabel(omega),ylabel(Hi)subplot(3,2,3);plot(w,Habs(:,3,

22、1),gridontitle(幅頻曲線3),xlabel(omega),ylabel(|H|)subplot(3,2,4);plot(w,Hang(:,3,1),gridontitle(相頻曲線3),xlabel(omega),ylabel(psi)subplot(3,2,6);plot(Hr(:,3,1),Hi(:,3,1),gridontitle(乃奎斯特圓3),xlabel(Re(H),ylabel(Im(H)w=sqrt(eig(K/M);未求和答案，請優(yōu)化程序并計算結(jié)果。實驗五：多自由度汽車振動分析未求和答案，請優(yōu)化程序并計算結(jié)果。一、實驗內(nèi)容給定發(fā)動機的參數(shù)書上的題3.5（p98）

23、二、實驗前期準備與內(nèi)容分析某4支撐發(fā)動機懸置系統(tǒng)，已知如下參數(shù)：質(zhì)量m轉(zhuǎn)動慣量J剛度矩陣K線位移矩陣X角位移矩陣0求系統(tǒng)的質(zhì)量矩陣、剛度矩陣、固有頻率和振型，并指出系統(tǒng)的耦合情況。三、實驗代碼與結(jié)果clearall;m=0.51;Jx=236;Jy=1087;Jz=957;Jxy=0;Jyz=0;Jxz=0;fori=1:4ou(i)=0;ov(i)=0;endku(1)=249;fori=1:4ou(i)=0;ov(i)=0;endku(1)=249;ku(2)=249;ku(3)=1849;ku(4)=1849;x(1)=53.17*0.01;x(2)=53.17*0.01;x(3)=-2

24、7.33*0.01;x(4)=-27.33*0.01;ov(3)=-3*pi/180;ou(2)=40*pi/180;os(i)=0;kv(1)=249;kv(2)=249;kv(3)=623;kv(4)=623;y(1)=26.18*0.01;y(2)=-27.42*0.01;y(3)=29.38*0.01;y(4)=-30.62*0.01;ov(4)=-3*pi/180;ou(3)=5*pi/180;ks(1)=878;ks(2)=878;ks(3)=971;ks(4)=971;z(1)=-14.6*0.01;z(2)=-14.6*0.01;z(3)=-11.6*0.01;z(4)=-11

25、.6*0.01;ou(1)=-40*pi/180;ou(4)=5*pi/180;%求系統(tǒng)的質(zhì)量矩陣、剛度矩陣、固有頻率和振型%求解過程與答案fori=1:6forj=1:6M(i,j)=0;endendM(1,1)=m;M(2,2)=m;M(3,3)=m;M(4,4)=Jx;M(5,5)=Jy;M(6,6)=Jz;M(4,5)=-Jxy;M(5,4)=-Jxy;M(4,6)=-Jxz;M(6,4)=-Jyz;M(5,6)=-Jyz;M(6,5)=-Jyz;Msum1=zeros(6,6)fork=1:4Cu(:,:,k)=100;0cos(ou(k)sin(ou(k);0-sin(ou(k)c

26、os(ou(k);Cv(:,:,k)=cos(ov(k)0-sin(ov(k);010;sin(ov(k)0cos(ov(k);Cs(:,:,k)=cos(os(k)sin(os(k)0;-sin(os(k)cos(os(k)0;001;C(:,:,k)=Cs(:,:,k)*Cu(:,:,k)*Cv(:,:,k);D(:,:,k)=ku(k)00;0kv(k)0;00ks(k)T(:,:,k)=1000z(k)-y(k);010-z(k)0 x(k);001y(k)-x(k)0;K(:,:,k)=(T(:,:,k)*(C(:,:,k)*D(:,:,k)*C(:,:,k)*T(:,:,k);en

27、dKK=sum(K,3)C=inv(M)*KK%A表示主振型；B%A表示主振型；B表示固%固有頻率%固有頻率C=K*inv(M)f1=w/2/pi1.0e+003*1.0e+003*0.00050.00050000.00050000.00050000000000000000000.23600001.08700000.95701.0e+0034.19120.1.0e+0034.19120.00320.0921-0.0002-0.47600.02510.00322.2691-0.06030.4601-0.01690.19920.0921-0.06033.1778-0.0267-0.13740.01

28、71-0.00020.4601-0.02670.3507-0.00110.1192-0.4760-0.0169-0.1374-0.00110.54900.00170.02510.19920.01710.11920.00170.7495w二90.743178.883266.64560.67110.83371.0771其中，質(zhì)量矩陣為M、剛度矩陣為KK、固有頻率W.實驗六：利用行列式值求矩陣特征根一、實驗內(nèi)容利用行列式值求矩陣特征根！二、實驗前期準備與內(nèi)容分析無阻尼三自由度系統(tǒng)運動微分方程：m、k、x分別為氣門、挺桿和搖臂等效質(zhì)量、等效剛度和垂直位移。F0-氣門推力-搖臂轉(zhuǎn)角系統(tǒng)運動微分方程求解：

29、尤1+xco%2+%3-m、k、x分別為氣門、挺桿和搖臂等效質(zhì)量、等效剛度和垂直位移。F0-氣門推力-搖臂轉(zhuǎn)角系統(tǒng)運動微分方程求解：尤1+xco%2+%3-nt兀2+(少；+謝孑)尤2123_閔強1如;歹(多)特征頻率方程：i少3兀20兀3+(況+w|)%3+殞刼給出：=0.175kg,m2=0.077kg,zn3=0.036kg,kD=34kN/miTi代=47kN/nmi，毎二總2Ql十kr27.8kN/mm求出固有頻率：01二0.72kHz。2=*7k缶-8.4kHz三、實驗代碼與結(jié)果%z,x,y方向的傳感器標定值分別為10.19，10.22，10.74Fs=1000;loadydata

30、(:,2)=y(2001:4000);data(:,3)=y(4001:6000);data(:,4)=y(6001:8000);x=data(:,2)/10.22;y=data(:,3)/10.74;z=data(:,4)/10.19;figure(1)subplot(3,1,1)plot(x)title(x)subplot(3,1,2)plot(y)title(y)subplot(3,1,3)plot(z)title(z)figure(2)subplot(3,1,1);pz,fz=filtbank(z,Fs);hz=bankdisp(pz,fz);title(座椅Z方向1/3倍頻程分析)；

31、wz=0.630.710.80.911110.80.630.50.40.3150.250.20.160.1250.1；subplot(3,1,2)；px,fx=filtbank(x,Fs)；hx=bankdisp(px,fx)；title(座椅X方向1/3倍頻程分析)；wx=110.80.630.500.400.3150.250.200.160.1250.100.080.0630.050.040.03150.025；subplot(3,1,3)；py,fy=filtbank(y,Fs)；hy=bankdisp(py,fy)；title(座椅Y方向1/3倍頻程分析)；wy=110.80.630.500.400.3150.25