汽车理论课后题matlab程序.doc

1.3ｎ=600:１:4000; r＝0．３67; i0＝5.83； eff＝0．８5; f＝0．013；ｍ=3880； g=9.８； G=m＊g; CdＡ=2．77; a＝１.94７； hg=０．9; Ｌ=3.2；Iw1=１.798; Ｉw2=３.５9８; Iｗ=Ｉｗ1+Iw2； If＝0.218;Ttq＝-１9.31３+29５.２7＊ｎ/1000-1６5.44*(n／100０）.＾２+40.874*(n/10０0).^3-３.8４４5*(n／１００0)．^４;%驱动力行驶阻力平衡图ｆor ｉg＝[5.5６,2．7６9，1．64４,1.00,0.７93]　 Uａ=0.３７7*ｒ*n/ig／ｉ0； 　 Ft=Tｔq*ｉｇ*i0*ｅff／r; 　 plot（Ua,Ft); hoｌd　ｏn;endFf=G*ｆ；uａ=0:0.1:ｍａx(Ｕa)；Fw=CdA＊ua.^2/２1．15;plot(ｕa,(Ff＋Ｆw)); titlｅ('驱动力-行驶阻力平衡图＇)；　ｘlａbeｌ('Ｕa/(km/h)'）；ｙｌａbeｌ('Ｆｔ/N'）;ｇtext('Ft１'）,ｇｔｅｘｔ('Ｆt2'),gｔext（'Ft3'）,gtｅxt('Ft4＇),ｇtｅxt（'Fｔ5＇），gｔeｘt('Ｆf+Fw'）[ｘ,y］=ginput(1)； dｉｓp（'汽车的最高车速');dｉｓp(x);disp（'kｍ/h＇)；％最大爬坡度及最大爬坡度时的附着率Ua＝0．３77＊r＊n/5.56/i0;Ft=Ttｑ*5．5６*i0*eff/r；Ｆw=CdA＊Ｕａ．^２/２1.15；i＝taｎ(aｓiｎ（（Ｆt－（Ff+Fw))/G）);dｉsｐ('汽车的最大爬坡度')；dｉsp（max(i))；C=ｍａx（i)/(a/Ｌ+hｇ/L*max(i));disp('克服最大爬坡度时的附着率＇);disp(C);%加速度倒数曲线fｉgｕre；for ｉg=[5.56,2．7６９,1．644,１.00,0.79３] Ua＝0.３７7＊r*n/ｉg/i0; q＝1＋Iｗ/(m*ｒ^２）+If*ig＾2*i0^２*eff/(m*r＾2)； Ft=Ttq*ｉg*i0*efｆ/ｒ；　 Fw=CdA*Ｕa．^２/２1．15; 　 aｓ＝（Ft-(Fｆ+Fw）)/q/m;　　 　plot(Ua，１./as）； holｄ　on；eｎdａxｉs([0　98　0 10]);tｉtle（'行驶加速度倒数曲线')；xlabel('Uａ／(km/h）');ylabeｌ（'1/a')；ｇtｅxt（'1/a１')，ｇtｅｘｔ('１/a2'),gtext（'1／ａ３＇）,gtext('1/a４'),gtext(＇1/ａ5＇);%加速时间曲线u(1)＝0.３７７*r*6０0/i0/２.769;dｔ＝0．01;w=1;t(1)=0；ig＝2.７69;whilｅ(u（w)<９9.08) 　 n（w)=u（w)*ig*i0／0.3７７/ｒ；　　 Ttq（w)=-1９.313＋295.27*n（w）/1000-16５.44*(n(w)／１0０0）^２+40.8７4*(n(w)／1000)＾３-３.844５*(n（w）／１０0０)＾4;　 Ft(ｗ)＝Ｔtq（w)＊ig*i0＊eff／r; 　 Ff=G*ｆ; Fｗ(w)=CｄA＊u(w)^2/21.15; 　 q＝1+Iw/（ｍ*ｒ＾2）+Ｉｆ*ig^2*i0^2*eff／(m*ｒ^2);　 a（ｗ）＝3.6*(Ft(w）-(Ｆf+Fw(w)）)/q／ｍ; 　 u（ｗ+１）=u(ｗ)+a（w)*dｔ;　 　　ｉf(u(ｗ+１)＞=０.３77*４0００＊r／2.769/i0) 　 　iｇ＝１．644; end 　if(u(ｗ+1)>＝0.３７７*400０*ｒ/1.644/i０) 　 ig=1.0０; 　eｎd 　if(u（w+1)>=0.37７*４000*r/1.0０/ｉ0)　 　 　ig=0.793; ｅnｄ　 　t(ｗ+１)=（w+1)*dt； 　 ｗ=w+１;endfigure;ploｔ(t,u);axis([0 60 0 １０0］）;tｉtlｅ('２挡原地起步加速时间曲线＇);xlaｂel('时间t/s＇），yｌabel（＇车速u/（ｋｍ／h)')；[x１,ｙ1]=gｉnpｕt(1)；disp('加速到７０km／h的时间')；disp（x1);1.8Iw1＝2.２5; Iw２=2.2５;　Ｉw=Ｉw1+Ｉｗ2； If=0．２５；　m＝1６0０; ｇ=9．８;G=m＊g；Ttqmaｘ=１40; i1=3.85;　i０=4．08； efｆ=0.9； r=0．3； hｇ＝0．63; b=1．２5; L=２.7；q=1+Iw/(m＊r＾2)+If*i１^2*i0^2＊eff／(m*r^2）;Ｆｔ=Ttｑmａx*i1＊i0*efｆ／r;a=Ｆt／ｑ／m；ｑ=a/g；C1=q/(b/L-ｈg/L*q);if（C1>０.6) 　 dｉsｐ（'发动机转矩不能充分发挥'）;eｌｓe 　 disp(＇发动机转矩能够充分发挥')；endb=（q／0.6+hg/Ｌ*ｑ）*L；diｓｐ（'调整后的b值');disp(b)；disｐ（'调整后的前轴负荷率')，disp(ｂ/Ｌ)；1.9Ｉｗ=3.6; If=0.２2； ｍ＝2152; ｇ=9．８;　r=０.367; efｆ=0.9; i4=1;　i0＝４.55;amax=0.75； Ｕａ=５0/3．6; CdＡ=1.5；q=1＋Iw／(m＊r^2);f=q*０．59／3．6/ｇ;diｓｐ('汽车在该路面上的滚动阻力系数');disp（ｆ)；ｑ=１＋Iw/(m*r^2)＋If*i4^２*i0^2＊ｅff／(m*r^2);D4ｍax＝f＋q*aｍax／g;disｐ('最大动力因数＇）;dｉsｐ(D4maｘ）;Fw=CｄA*Uａ^２/21．１5;D1max=（Ｄ4max＊G+Fw)*i1／i４／G;iｍax＝tan(asiｎ((D1ｍax－ｆ＊ｓqrt(1-D１maｘ^2＋f^2))/(１+f^2))）;dｉｓp('最大爬坡度'),diｓp(imaｘ）;2．7n=６00：１：40０0；　r=０.367; i0＝5.８3; eｆｆ=０.85; f=０．０1３;CdA＝２.77;m＝3８８0; g＝9.8; G=m*g；Ｔtq=-19.３１３+295.２7*n/10００-16５．44*(ｎ/1000)．^２+４0.874*（n/100０).^3-3．84４5*(ｎ/１000).＾4;for iｇ=[５.56,２．７6９,1.６44，1.0０，0.79３] 　Uａ=0.３7７*r＊n/ig／i0;　　 Pｅ=Ttq．*n/９５５0; ploｔ（Ｕａ，Pe); ｈｏld on;ｅｎdUa=0:0.1：ｍax(Ua)；Pf=G*ｆ*Ua/36０0;Pw=ＣdA*Ua.^3/76１40；pｌｏｔ(Ua,（Ｐf＋Ｐw)／eff);titlｅ(＇汽车的功率平衡图'），xlabel（'Ua／(km/h)＇),ylａbel（'Ｐ/kw');％等速百公里燃油消耗n1=[８15,１207,１614,2012,2603，3006,３403，3８04］；b0=［1326．8,135４.7，１284.４,１122.9,１１４1.０,1０5１.２,１233．9，1１２９.7]；b１=［－４1６.46,-30３.98,-1８9．７5,－１２1．59,－９8．89３，－7３.7１4,－84．4７８,－45．29１］ ;b2=［７2．３７９,３6.657,14.524，７.0035,４.476３,2.85９3,２.9７88，0．7１113];b3＝[-５.８6２9,-2.05５3,-0.511８４，-０．１8517,－0．09１0７７，-0.05１38,-0.0４7４49，-０.00０７5215］;b4=[0．17７6８,0.0４３072,0.00681６４,0.0018555，0.０0０6８９06,０.0０03５03２,0.00０282３０，－0.０00０385６8];Ｂ０=sｐline（ｎ1，b０，n); B1=spline（n1,ｂ1,n)；　B２=spliｎe(ｎ1,ｂ2,n)； B3=spline(n1,b３，n)；　B4=ｓpline(ｎ1，b４,ｎ）;ｄeng=７；figuｒｅ;for ｉg=[1.０0,０.７９３] 　　 Uａ=0.377*r*n/ｉg/i0; Pｆ=Ｇ＊ｆ*Ua/3600； Pw=CdA*Ua.^３/76140; P=（Pｆ+Pw）/efｆ;　 b＝Ｂ0＋B１.*Ｐ+B２．*Ｐ.^2+B3.*P.^3＋B４．*P.＾4; 　 Q=P.*b.／(1．02*Uａ*deng); plot（Uａ,Q)； hｏld　ｏn;enｄaxis(［０ 100　１2 ２8]);title（'最高挡与次高挡等速百公里油耗曲线'),xlabel('Uａ/(km/h）'),ｙｌａbｅl('百公里油耗Qs/（Ｌ/100ｋm）');4.3ｂeta＝0.３８;　L＝3.9５； a1＝2.1; b1=L-a1;　hｇ1=0.84５;a2＝2．95； b2=L-a２; hg2＝1.17;ｇ=9.８;z=0:０.０1:1;ｆａiｆ1=bｅｔa*ｚ./（1/Ｌ*(b１+ｚ*hg１));fａir1＝(1-beta)*z.／(1/L*（a1-ｚ*hｇ1）);faｉf2=beｔa*z．/(1/L＊(ｂ２+z*ｈg2）);fair2=(1-betａ)＊z./(１/L*(a２-z*hg2))；pｌot(z,faｉf1，'r',z,fair1，'g＇，z,fａif2,＇ｂ',z,faｉr2，'k')；ｆａi=０:０.０1：１；Ef1＝(ｂ1/L)．/(ｂｅta-fai*hｇ1/L)；Eｒ1=(a1/Ｌ).／((1-betａ)+fai*hg１/L）;Ef２＝(b2／L)./(betａ-fａi＊hg２/Ｌ)；Ｅr2=(ａ2/L)./((1-ｂetａ)+fai*hg2／Ｌ);ｆigure;plot(faｉ,Ｅf１,＇r＇,fai,Ｅr１,'g'，ｆａi,Ｅf2,'b'，fai，Eｒ2，'k');aｘｉs（[0　1 0 １]);fai＝0.8；Er1=（a1／Ｌ)／((1－bｅｔa)+ｆai*ｈg1／L);Eｒ2=（ａ２/Ｌ)/((１－bｅta）+fai*hg2/L);ａｂmａx１＝fai*Eｒ1＊ｇ;ａbmax2=fai*Er2*g;t21=0.02; t22=０.2；　u0=３０;s1=1/3.6*(t21+t22/2）＊u0+u0^2/25．92/ａｂｍａx１;s2=1/3．6*(t2１+t22/２）*ｕ0＋u0^2/2５.９2/abmax2;dｉｓp(＇空载时的制动距离'),ｄiｓp(s1）;disp('满载时的制动距离'),ｄisp(s2);%前部管路损坏时制动距离z1=a1*faｉ／（L+ｆaｉ*hg1); z2=ａ２＊ｆai/（Ｌ+fai*hg2);ab。