資源描述:
《汽車動力性matlab仿真源程序.doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、clearallclcn=[1500:500:5500];%轉(zhuǎn)速范圍T=[78.5983.0485.0186.6387.0985.8784.6782.5080.54];%對應(yīng)各轉(zhuǎn)矩dt=polyfit(n,T,3);%對發(fā)動機輸出轉(zhuǎn)矩特性進行多項式擬合,階數(shù)取4n1=1000:100:5500;%????t=polyval(dt,n1);figure(1)title('發(fā)動機外特性')plot(n1,t,n,T,'o'),gridon%圖示發(fā)動機輸出轉(zhuǎn)矩特性%汽車驅(qū)動力計算G=input('整車重力/N,G=');%輸入970*9.8ig=[3.4161.8941.280.
2、9140.757];%變速器速比k=1:5;%5個前進檔r=0.272;i0=4.388;eta=0.9;ngk=[800800800800800];ngm=[55005500550055005500];ugk=0.377.*r.*ngk(k)./(ig(k).*i0);%計算每一檔發(fā)動機800rpm時的最低行駛速度ugm=0.377.*r.*ngm(k)./(ig(k).*i0);%計算每一檔發(fā)動機5400rpm最高行駛速度fork=1:5%依次計算5個檔的驅(qū)動力u=ugk(k):ugm(k);n=ig(k)*i0.*u./r/0.377;t=54.8179+2.2441
3、.*(n./100)-4.8003.*(n./1000).^2+2.815e-10.*n.^3Ft=t.*ig(k).*i0*eta/r;figure(2)plot(u,Ft)holdon,gridon%保證K的每次循環(huán)的圖形都保留顯示end%行駛阻力計算f0=0.009;f1=0.002;f4=0.0003;%三者都是轎車滾動阻力系數(shù)%disp'空氣阻力系數(shù)Cd=0.3--0.41,迎風(fēng)面積A=1.7--2.1'Cd=input('空氣阻力系數(shù)Cd=');%輸入0.3A=input('迎風(fēng)面積/m2,A=');%輸入2.3u=0:10:180;f=f0+f1.*(u./1
4、00)+f4.*(u./100).^4;Ff=G*f;%計算滾動阻力Fw=Cd*A.*u.^2./21.15;%計算空氣阻力F=Ff+Fw;%滾動阻力、空氣阻力之和title('驅(qū)動力-阻力圖(五檔速比為3.4161.8941.280.9140.757)')plot(u,F,'mo-');gridonholdon%圖解法求最高車速foru=50:180;k=5;n=ig(k)*i0.*u./r/0.377;t=54.8179+2.2441.*(n./100)-4.8003.*(n./1000).^2+2.815e-10.*n.^3;Ft=t.*ig(k).*i0*eta/r
5、;f=f0+f1.*(u./100)+f4.*(u./100).^4;Ff=G*f;Fw=Cd*A.*u.^2./21.15;F=Ff+Fw;ifabs(Ft-F)<1;umax=u;breakendenddisp('=====汽車動力性能仿真計算結(jié)果=====')disp('驅(qū)動力-行駛阻力平衡圖及最高車速')fprintf('汽車?yán)碚撟罡哕囁伲?qū)動力與行駛阻力曲線交點)Vmax=%3.3fkm/h',umax)