汽車動力性計(jì)算matlab程序

1、%汽車動力性計(jì)算(自己編的動力性計(jì)算程序,供大家計(jì)算動力性時(shí)參考,具體參數(shù)大家根據(jù)所給程序?qū)?yīng)輸入,并對坐標(biāo)軸數(shù)值按需要進(jìn)行修改)clc;clear;closeall;%艮據(jù)所給發(fā)動機(jī)數(shù)據(jù)擬合外特性曲線(發(fā)動機(jī)數(shù)據(jù)根據(jù)你所得到的數(shù)據(jù)進(jìn)行輸入)n_test=5006007008009001000110012001300140015001600170018001900200021002200;T_test=97511081298149615461620217017851974197419701889182917481669170015241105;figure(1)plot(n_test,T_te

2、st,'g');holdongridon%p=polyfit(n_test,T_test,7);p=polyfit(n_test,T_test,2);n=450:1:2200;Ttq=polyval(p,n);plot(n,Ttq,'k');xlabel('發(fā)動機(jī)轉(zhuǎn)速n(r/min)');ylabel('發(fā)動機(jī)轉(zhuǎn)矩Ttq(N*m)');title('發(fā)動機(jī)轉(zhuǎn)矩曲線)；legend('測試曲線,擬合曲線)；%幀給車型動力總成相關(guān)參數(shù)ig=3.072.161.481.00.82;i0=4.0;eta=0.78;r=0.

3、57;M=25000;g=9.8;c=1.5;f0=0.01;f1=0.0002;f4=0.0005;CD=1;A=8;Iw=3.6;If=0.04;%發(fā)動機(jī)外特性曲線圖figure(2)holdongridonfori=length(n);Pe=Ttq.*n/9550;endAX,H1,H2=plotyy(n,Ttq,n,Pe);xlabel('發(fā)動機(jī)轉(zhuǎn)速n(r/min)');ylabel('發(fā)動機(jī)轉(zhuǎn)矩Ttq(N*m)');ylabel(AX(2),'發(fā)動機(jī)功率Pe(Kw)');title('發(fā)動機(jī)外特性曲線)；%各擋位速度曲線%計(jì)算各

4、擋位車速fori=1:length(ig);ua(i,:)=0.377Tn/ig(i)/i0;end%計(jì)算各檔位最高車速uamax=ua(:,length(ua(1,:);figure(3)holdonfori=1:length(ig);plot(n,ua(i,:),'k');endholdongridonxlabel('轉(zhuǎn)速n(r/min)');ylabel('各擋位車速(km/h)');title('各擋位車速-轉(zhuǎn)速表');legend('1擋車速,2擋車速','3擋車速','4擋車速,

5、5擋車速);%驅(qū)動力和行駛阻力平衡圖%計(jì)算滾動阻力系數(shù)fori=1:length(ig);f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100)A4;end%計(jì)算滾動阻力fori=1:length(ig);Ff(i,:)=c*M*g*f(i,:);end%計(jì)算空氣阻力fori=1:length(ig);Fw(i,:)=CD*A*(ua(i,:).A2)/21.15;end%計(jì)算行駛阻力fori=1:length(ig);F(i,:)=Ff(i,:)+Fw(i,:);end%計(jì)算汽車驅(qū)動力fori=1:length(ig);Ft(i,:)=Ttq*ig(i)*i0

6、*eta/r;endfigure(4)holdonfori=1:length(ig);plot(ua(i,:),Ft(i,:),plot(ua(i,:),F(i,:),'r');plot(ua(i,:),Ff(i,:),'b');endholdongridonxlabel('車速(km/h)');ylabel('驅(qū)動力、行駛阻力(N)');legend('驅(qū)動力Ft','行駛阻力Ff+Fw','滾動阻力Ff');title('驅(qū)動力-行駛阻力平衡圖)；%汽車功率平衡圖%計(jì)算各

7、檔位功率fori=1:length(ig);P(i,:尸F(xiàn)t(i,:).*ua(i,:)/(3600*eta);end%計(jì)算風(fēng)阻阻力功率fori=1:length(ig);Pw(i,:)=CD*A*ua(i,:).A3/(76140*eta);end%計(jì)算滾動阻力功率fori=1:length(ig);Pf(i,:)=M*g*f(i,:).*ua(i,:)/(3600*eta);end%計(jì)算總阻力功率fori=1:length(ig);Pz(i,:尸Pw(i,:)+Pf(i,:);endfigure(5)holdonfori=1:length(ig);plot(ua(i,:),P(i,:),&

8、#39;k');plot(ua(i,:),Pz(i,:),'r');endholdongridonxlabel('車速(km/h)');ylabel('發(fā)動機(jī)功率、阻力功率(kW)');legend('發(fā)動機(jī)功率P','阻力功率Pz','Location','Northwest');title('功率平衡圖)；%動力特性圖(動力因數(shù)圖)fori=1:length(ig);D(i,：)=(Ft(i,:)-Fw(i,:)/M/g;endfigure(6)holdonfor

9、i=1:length(ig);plot(ua(i,:),D(i,:),'k');plot(ua(i,:),f(i,:),'r');endholdongridonxlabel('車速(km/h)');ylabel('動力因數(shù)D');legend('動力因數(shù)D','滾動阻力系數(shù)f');title('動力特性圖)；%爬坡度曲線圖fori=1:length(ig);I(i,:)=(tan(asin(Ft(i,:)-(Ff(i,:)+Fw(i,:)/(M*g)*100;endfigure(7)hold

10、onfori=1:length(ig);ifi=1plot(ua(i,:),I(i,:),'r');elseplot(ua(i,:),I(i,:),'k');endendholdongridonxlabel('車速(km/h)');ylabel('最大爬坡度(%)');legend('I擋,高速檔)；title('爬坡度曲線圖)；%加速度曲線圖deta=1+1/M*4*Iw/rA2+1/M*If*ig.A2*i0A2*eta/rA2;fori=1:length(ig);a(i,:)=(Ft(i,:)-Ff(i,:)

11、-Fw(i,:)./deta(i)/M;ifi=5forj=1:length(n)ifa(i,j)<0a(i,j)=0;elseendendendendfigure(8)holdonfori=1:length(ig);ifi=1plot(ua(i,:),a(i,:),'r');elseplot(ua(i,:),a(i,:),'k');endendholdongridonxlabel('車速(km/h)');ylabel('加速度a(m/sA2)');legend(I檔,高速檔)；title('加速度曲線圖)；axis

12、(012001,5)%加速度倒數(shù)曲線fori=1:length(ig);forj=1:length(n)b(i,j)=1,/a(i,j);endendfigure(9)holdonfori=1:length(ig)plot(ua(i,:),b(i,:),'k');endholdongridonxlabel('車速(km/h)');ylabel('各檔加速度倒數(shù)1/a');legend('各檔加速度倒數(shù)1/a曲線,Location','Northwest');title('各檔加速度倒數(shù)曲線圖)；axis(0

13、120010)ad1=b(1,:);ad2=ua(1,:);fori=1:(length(ig)-1);forj=1:length(n)ifua(i+1,j)>=ua(i,length(n)flag(i)=j；break;endendad1=ad1b(i+1,j:length(n);ad2=ad2ua(i+1,j:length(n);endfigure(10)holdonplot(ad2,ad1,'k');holdongridonxlabel('車速(km/h)');,'NorthWest');ylabel('加速度倒數(shù)1/a

14、9;);legend('加速度倒數(shù)1/a曲線','Location'title('加速度倒數(shù)曲線圖)；axis(0120010)%加速時(shí)間曲線k=length(n);fori=1:length(ig);t(i,1)=0;forj=2:kt(i,j)=abs(ua(i,j)-ua(i,j-1)*(b(i,j)+b(i,j-1)/2;endendfori=1:length(ig);forj=1:kat(i,j)=sum(t(i,1:j)/3.6;endendtotalat=at(1,:);fori=1:(length(ig)-1);forj=flag(i):ktotalat=totalattotalat(length(totalat)+t