《matlab-汽車?yán)碚?作業(yè)(word文檔良心出品).doc》由會員分享，可在線閱讀，更多相關(guān)《matlab-汽車?yán)碚?作業(yè)(word文檔良心出品).doc（4頁珍藏版）》請?jiān)趨R文網(wǎng)上搜索。

1、%習(xí)題1.3%平衡圖clear;n=600:1:4000;T=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;r=0.367;e=0.85;io=5.83;f=0.013;CdA=2.77;G=38024;m=3880;If=0.218;Iw1=1.798;Iw2=3.598;Ig=5.56,2.769,1.644,1.00,0.793;ig1=5.56;ig2=2.769;ig5=0.793;for ig=Ig; Ft=T*ig*io*e/r; Ua=0.377*r*n/(io*ig)

2、; subplot(1,2,1); plot(Ua,Ft); axis square; title(驅(qū)動力與行駛阻力平衡圖); xlabel(Ua/(km/h) ylabel(F/N) hold on;endF=G*f+(CdA*Ua.2)/21.15;hold on; plot(Ua,F,r);Ft=T*ig5*io*e/r;k=abs(Ft-F);imin=find(min(k)=k);text(Ua(imin),Ft(imin), max=,num2str(Ua(imin);plot(Ua(imin),Ft(imin),r*);fprintf(最高車速為：);fprintf(%g,Ua(

3、imin);fprintf(km/hn);%爬坡度Ft1=T*ig1*io*e/r;Ua=0.377*r*n/(io*ig1);F=G*f+(CdA*Ua.2)/21.15;alpha=asin(Ft1-F)/G);i=tan(alpha);imax=find(max(i)=i);fprintf(最大爬坡度為：);fprintf(%g,i(imax);%附著率a1=1.947;L=3.2;hg=0.9;q=i+(1./cos(alpha)*(1./9.8);C=q/(a1/L+hg*q/L);fprintf(n附著率為：);fprintf(%g,C);fprintf(n);%加速度倒數(shù)曲線fo

4、r ig=Ig; Ft=T*ig*io*e/r; Ua=0.377*r*n/(io*ig); F=G*f+(CdA*Ua.2)/21.15; delta=1+(Iw1+Iw2)/(m*r2)+(If*ig2*io2*e)/(m*r2); a=(Ft-F)/(m*delta); v=1./a; subplot(1,2,2); plot(Ua,v); axis square; axis(0 99 0 10) title(汽車行駛加速度倒數(shù)曲線); xlabel(Ua/(km/h) ylabel(1/a) hold onend%加速時間yitaT=0.85;r=0.367;f=0.013;CdA=2

5、.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;ig=5.56 2.769 1.644 1.00 0.793;nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5 deta(i)=1+(Iw1+Iw2)/(m*r2)+(If*(ig(i)2*i02*yitaT)/(m*r2);enddeta_u=0.01;ua=6:deta_u:99;N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*

6、ua;delta=0*ua;Ff=G*f;Fw=CdA*ua.2/21.15;for i=1:N k=i; if ua(i)=u2(2) n=ua(i)*(ig(2)*i0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig(2)*i0*yitaT/r; inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i); delta(i)=deta_u*inv_a(i)/3.6; elseif ua(i)=u2(3) n=ua(i)*(ig(3)*i

7、0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig(3)*i0*yitaT/r; inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i); delta(i)=deta_u*inv_a(i)/3.6; elseif ua(i)=u2(4) n=ua(i)*(ig(4)*i0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/10

8、00)4; Ft=Tq*ig(4)*i0*yitaT/r; inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i); delta(i)=deta_u*inv_a(i)/3.6; else n=ua(i)*(ig(5)*i0/r)/0.377; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)2+40.874*(n/1000)3-3.8445*(n/1000)4; Ft=Tq*ig(5)*i0*yitaT/r; inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i); delta(i)=deta_u*inv_a(i)/3.6; end

9、a=delta(1:k); t(i)=sum(a);endfigure(2);plot(t,ua);axis(0 40 0 70);title(汽車加速時間曲線);xlabel(t/s);ylabel(ua/km/h);fprintf(由圖可知，加速到70KM/h大約需要25.72sn);%習(xí)題1.8clear;m=1600;a=1450;b=1250;L=2700;hg=630;Mmax=140;ig1=3.85;i0=4.08;e=0.9;r=0.3;If=0.25;Iw1=2.25;Iw2=2.25;delta=1+(Iw1+Iw2)/(m*r2)+(If*ig12*i02*e)/(m*

10、r2);Ft=Mmax*ig1*i0*e/r;Acc=Ft/(delta*m);Fz1=m*9.8*b/L-Acc*(m*hg/L+(Iw1+Iw2+If*i0*ig1)/(L*r);q=Acc/9.8;C1=q/(b/L-hg*q/L);fprintf(最大加速度時需要的附著率為:);fprintf(%g,C1);fprintf(n此時的附著率大于0.6，故不能產(chǎn)生最大的加速度。n);b1=q*L/0.6+hg*q;x=b1/L;fprintf(調(diào)整后b值為：);fprintf(%g,b1);fprintf(n此時前軸負(fù)荷率為：);fprintf(%gn,x);%習(xí)題1.9clear;m=2

11、152;i0=4.55;ig1=3.79;ig2=2.17;ig3=1.41;ig4=1;ig5=0.86;hg=0.57;CdA=1.5;L=2300;If=0.22;Iw=3.6;r=0.367;e=0.9;delta=1+(Iw)/(m*r2);Acc=0.59/3.6;f=Acc*delta/9.8;fprintf(滾動阻力系數(shù)為：);fprintf(%gn,f);delta=1+(Iw)/(m*r2)+(If*ig42*i02*e)/(m*r2);D=f+delta*0.75/9.8;fprintf(最大動力因數(shù)為：);fprintf(%gn,D);q=0.52/(1/0.6-hg/

12、L);i=ig1*D-f;fprintf(最大爬坡度為：);fprintf(%gn,min(i,q);%習(xí)題2.7%功率平衡圖clear;n0=600:1:4000;T=-19.313+295.27*(n0/1000)-165.44*(n0/1000).2+40.874*(n0/1000).3-3.8445*(n0/1000).4;r0=0.367;e=0.85;io=5.83;f0=0.013;CdA=2.77;G0=38024;m0=3880;If0=0.218;Iw10=1.798;Iw20=3.598;Ig=5.56,2.769,1.644,1.00,0.793;ig5=0.793;i

13、g4=1;for ig=Ig; Ua=0.377*r0*n0/(io*ig); Ft=T*ig*io*e/r0; Pe=Ua.*Ft/3600; subplot(1,2,1); plot(Ua,Pe); title(功率平衡圖); xlabel(Ua/(km/h) ylabel(Pe/KW) hold on;endP=(G0*f0*Ua/3600+CdA*Ua.3/76140)./e;plot(Ua,P,r);%第二問clear;global f G CDA nT m r If Iw1 Iw2 pg i0 B0 B1 B2 B3 B4 n ig=5.56,2.769,1.644,1.00,0.

14、793;r=0.367;pg=7;i0=5.83;nT=0.85;CDA=2.77;f=0.013;G=(3880)*9.8;If=0.218;Iw1=1.798;Iw2=3.598;m=3880; n0=815 1207 1614 2012 2603 3006 3403 3804;B00=1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7;B10=-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291;B20=72.379 36.657 14.524 7.003

15、5 4.4763 2.8593 2.9788 0.71113;B30=-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215;B40=0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568;n=600:1:4000;B0=spline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n); B3=spline(n0,B30,n);B4=spline

16、(n0,B40,n);ua3=0.377*r.*n./(i0*ig(3);ua4=0.377*r.*n./(i0*ig(5);F3=f*G+CDA*(ua3.2)/21.15;F4=f*G+CDA*(ua4.2)/21.15; Pe3=F3.*ua3./(nT*3.6*1000); Pe4=F4.*ua4./(nT*3.6*1000);for i=1:1:3401 b3(i)=B0(i)+B1(i)*Pe3(i)+B2(i)*(Pe3(i)2+B3(i)*(Pe3(i)3+B4(i)*(Pe3(i)4;b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*(Pe4(i)2+B3(i)*

17、(Pe4(i)3+B4(i)*(Pe4(i)4;endQ3=Pe3.*b3./(1.02.*ua3.*pg);Q4=Pe4.*b4./(1.02.*ua4.*pg);subplot(1,2,2);plot(ua3,Q3,r);axis(0 110 10 30)hold on;plot(ua4,Q4);title(等速百公里油耗曲線);xlabel(Ua(km/h);ylabel(百公里油耗(L/100km);text(42,21,四擋),text(62,17,五擋);%習(xí)題4.3%第一問clear;me=4080;hge=0.845;Le=3.950;ae=2.10;betae=0.38;be

18、=Le-ae;mf=9290;hgf=1.170;Lf=3.950;af=2.950;betaf=0.38;bf=Lf-af;z=0:0.01:1;faife=betae*z*Le./(be+z*hge);%空載時的ffaiff=betaf*z*Lf./(bf+z*hgf);%滿載時的ffaire=(1-betae)*z*Le./(ae-z*hge);%空載時的rfairf=(1-betaf)*z*Lf./(af-z*hgf);%滿載時的ffai=z;figure(1);plot(z,faife,k-.,z,faiff,b,z,faire,k-.,z,fairf,b,z,fai,k);axis

19、 square;text(0.6,1.3,r(空車),text(0.62,0.9,r(滿載),text(0.92,1.01,=z),text(0.85,0.8,f(滿載),text(0.85,0.38,f(空車);title(利用附著系數(shù)與制動強(qiáng)度的關(guān)系曲線);xlabel(制動強(qiáng)度z/g);ylabel(利用附著系數(shù));Efm=z./faiff*100;Erm=z./fairf*100;Erk=z./faire*100;Efk=z./faife*100;figure(2);plot(faiff,Efm,b,fairf,Erm,b,faire,Erk,k,faife,Efk,k);axis s

20、quare;axis(0 1 0 100);text(0.25,88,Ef),text(0.6,97,滿載 ),text(0.75,92,Er),text(0.65,2,Er),text(0.55,75,Er),text(0.55,65,空載);title(前后制動效率曲線);xlabel(附著系數(shù));ylabel(制動效率%);%第二問t2=0.02+0.5*0.2;Ua=30;aMaxe=0.68*0.8*9.8;aMaxf=0.87*0.8*9.8;se=t2*Ua/3.6+(Ua2)/(25.92*aMaxe);sf=t2*Ua/3.6+(Ua2)/(25.92*aMaxf);fpri

21、ntf(空載時車輪不抱死的制動距離為：);fprintf(%gn,se);fprintf(滿載時車輪不抱死的制動距離為：);fprintf(%gnn,sf);%第三問%前部損壞ze=ae*0.8/(Le+0.8*hge);zf=af*0.8/(Lf+0.8*hgf);amaxe=ze*9.8;amaxf=zf*9.8;se1=t2*Ua/3.6+(Ua2)/(25.92*amaxe);sf1=t2*Ua/3.6+(Ua2)/(25.92*amaxf);fprintf(前部管路損壞時，空載，車輪不抱死的制動距離為：);fprintf(%gn,se1);fprintf(前部管路損壞時，滿載，車輪不

22、抱死的制動距離為：);fprintf(%gn,sf1);%后部損壞ze=be*0.8/(Le-0.8*hge);zf=bf*0.8/(Lf-0.8*hgf);amaxe=ze*9.8;amaxf=zf*9.8;se2=t2*Ua/3.6+(Ua2)/(25.92*amaxe);sf2=t2*Ua/3.6+(Ua2)/(25.92*amaxf);fprintf(后部管路損壞時，空載，車輪不抱死的制動距離為：);fprintf(%gn,se2);fprintf(后部管路損壞時，滿載，車輪不抱死的制動距離為：);fprintf(%gn,sf2);%習(xí)題4.5clear;m=1600;a=1450;b

23、=1250;L=2700;hg=630;Mmax=140;ig1=3.85;i0=4.08;e=0.9;r=0.3;If=0.25;Iw1=2.25;Iw2=2.25;beta=0.65;%第一問fai0=(L*beta-b)/hg;fprintf(同步附著系數(shù)為：);fprintf(%gn,fai0);%第二問Ef=b/(L*beta-0.7*hg);fprintf(制動效率為：);fprintf(%gn,Ef);%第三問abmax=9.8*Ef*0.7;fprintf(最大制動減速度為：);fprintf(%g,abmax);fprintf( m/s2nn)%第四問fprintf(制動系統(tǒng)

24、增益 a b cn);fprintf(雙回路系統(tǒng) G G Gn);fprintf(1回路失效時 0.7G G Gn);fprintf(2回路失效時 1.3G G Gnn);%第五問%a系統(tǒng)前部回路失效時z1=a*0.7/(L+0.7*hg);abmaxa1=z1*9.8;Ef1=z1/0.7;%a系統(tǒng)后部回路失效時z2=b*0.7/(L-0.7*hg);abmaxa2=z2*9.8;Ef2=z2/0.7;%對于b c 因?yàn)閜hiphi0,前輪先抱死,當(dāng)車輪剛要抱死時zbc=0.7*b/(2*L*beta-0.7*hg);abmaxbc=zbc*9.8;Erbc=zbc/0.7;fprintf(

25、對于a系統(tǒng)前部回路失效時制動效率為);fprintf(%gn,Ef1);fprintf(對于a系統(tǒng)前部回路失效時最大制動減速度為);fprintf(%gn,abmaxa1);fprintf(對于a系統(tǒng)后部回路失效時制動效率為);fprintf(%gn,Ef2);fprintf(對于a系統(tǒng)后部回路失效時最大制動減速度為);fprintf(%gn,abmaxa2);fprintf(對于b、c系統(tǒng)某一個回路失效時制動效率為);fprintf(%gn,Erbc);fprintf(對于b、c系統(tǒng)某一個回路失效時最大制動減速度為);fprintf(%gn,abmaxbc);%第六問A：一個回路失效時，不易發(fā)生跑偏。B：左右制動力不易相等。C：同上。