汽車動力性計算matlab程序

1、精選優(yōu)質(zhì)文檔-傾情為你奉上% 汽車動力性計算（自己編的動力性計算程序，供大家計算動力性時參考，具體參數(shù)大家根據(jù)所給程序?qū)斎?，并對坐標軸數(shù)值按需要進行修改） clc; clear; close all; %根據(jù)所給發(fā)動機數(shù)據(jù)擬合外特性曲線（發(fā)動機數(shù)據(jù)按照你所得到的數(shù)據(jù)進行輸入） n_test=500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200; T_test=975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889

2、1829 1748 1669 1700 1524 1105; figure(1) plot(n_test,T_test,'g'); hold on grid on %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ā)動機轉(zhuǎn)速n(r/min)'); ylabel('發(fā)動機轉(zhuǎn)矩Ttq(N*m)'); title('發(fā)動機轉(zhuǎn)矩曲線'); legen

3、d('測試曲線','擬合曲線'); %所給車型動力總成相關(guān)參數(shù) ig=3.07 2.16 1.48 1.0 0.82; i0=4.0; eta=0.78; r=0.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ā)動機外特性曲線圖 figure(2) hold on grid on for i=length(n); Pe=Ttq.*n/9550; end AX,H1,H2=plotyy(n,Ttq,n,Pe); xlabel('發(fā)動機

4、轉(zhuǎn)速n(r/min)'); ylabel('發(fā)動機轉(zhuǎn)矩Ttq(N*m)'); ylabel(AX(2),'發(fā)動機功率Pe(Kw)'); title('發(fā)動機外特性曲線'); % 各擋位速度曲線 %計算各擋位車速 for i=1:length(ig); ua(i,:)=0.377*r*n/ig(i)/i0; end%計算各檔位最高車速 uamax=ua(:,length(ua(1,:); figure(3) hold on for i=1:length(ig); plot(n,ua(i,:),'k'); end hold o

5、n grid on xlabel('轉(zhuǎn)速n(r/min)'); ylabel('各擋位車速(km/h)'); title('各擋位車速-轉(zhuǎn)速表'); legend('1擋車速','2擋車速','3擋車速','4擋車速','5擋車速'); % 驅(qū)動力和行駛阻力平衡圖 %計算滾動阻力系數(shù) for i=1:length(ig); f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100).4; end %計算滾動阻力 for i=1:length

6、(ig); Ff(i,:)=c*M*g*f(i,:); end %計算空氣阻力 for i=1:length(ig); Fw(i,:)=CD*A*(ua(i,:).2)/21.15; end %計算行駛阻力 for i=1:length(ig); F(i,:)=Ff(i,:)+Fw(i,:);end %計算汽車驅(qū)動力 for i=1:length(ig); Ft(i,:)=Ttq*ig(i)*i0*eta/r;end figure(4) hold on for i=1:length(ig); plot(ua(i,:), Ft(i,:),'k'); plot(ua(i,:), F

7、(i,:),'r'); plot(ua(i,:), Ff(i,:),'b'); end hold on grid on xlabel('車速(km/h)'); ylabel('驅(qū)動力、行駛阻力(N)'); legend('驅(qū)動力Ft','行駛阻力Ff+Fw','滾動阻力Ff'); title('驅(qū)動力-行駛阻力平衡圖'); % 汽車功率平衡圖 %計算各檔位功率 for i=1:length(ig); P(i,:)=Ft(i,:).*ua(i,:)/(3600*eta)

8、;end %計算風阻阻力功率 for i=1:length(ig); Pw(i,:)=CD*A*ua(i,:).3/(76140*eta); end %計算滾動阻力功率 for i=1:length(ig); Pf(i,:)=M*g*f(i,:).*ua(i,:)/(3600*eta); end %計算總阻力功率 for i=1:length(ig); Pz(i,:)=Pw(i,:)+Pf(i,:);end figure(5) hold on for i=1:length(ig); plot(ua(i,:), P(i,:),'k'); plot(ua(i,:), Pz(i,:)

9、,'r'); end hold on grid on xlabel('車速(km/h)'); ylabel('發(fā)動機功率、阻力功率(kW)'); legend('發(fā)動機功率P','阻力功率Pz','Location','NorthWest'); title('功率平衡圖'); % 動力特性圖（動力因數(shù)圖） for i=1:length(ig); D(i,:)= (Ft(i,:)- Fw(i,:)/M/g; end figure(6) hold on for i=1:l

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

11、00; end figure(7) hold on for i=1:length(ig); if i=1 plot(ua(i,:),I(i,:),'r'); else plot(ua(i,:),I(i,:),'k'); end end hold on grid on xlabel('車速(km/h)'); ylabel('最大爬坡度（%）'); legend('擋','高速檔'); title('爬坡度曲線圖'); % 加速度曲線圖 deta=1+1/M*4*Iw/r2+1/M*If

12、*ig.2*i02*eta/r2; for i=1:length(ig); a(i,:)=(Ft(i,:)-Ff(i,:)-Fw(i,:)./deta(i)/M; if i=5 for j=1:length(n)if a(i,j)<0 a(i,j)=0; else end end end end figure(8) hold on for i=1:length(ig); if i=1 plot(ua(i,:),a(i,:),'r'); else plot(ua(i,:),a(i,:),'k'); end end hold on grid onxlabel(

13、'車速(km/h)'); ylabel('加速度a（m/s2）'); legend('檔','高速檔'); title('加速度曲線圖'); axis(0 120 0 1.5)% 加速度倒數(shù)曲線 for i=1:length(ig); for j=1:length(n) b(i,j)=1./a(i,j); end end figure(9) hold on for i=1:length(ig) plot(ua(i,:),b(i,:),'k'); end hold on grid on xlabel(

14、'車速(km/h)'); ylabel('各檔加速度倒數(shù)1/a'); legend('各檔加速度倒數(shù)1/a曲線','Location','NorthWest'); title('各檔加速度倒數(shù)曲線圖'); axis(0 120 0 10)ad1=b(1,:); ad2=ua(1,:); for i=1:(length(ig)-1);for j=1:length(n) if ua(i+1,j)>=ua(i,length(n) flag(i)=j; break; end end ad1=ad1 b

15、(i+1,j:length(n); ad2=ad2 ua(i+1,j:length(n); end figure(10) hold on plot(ad2,ad1,'k'); hold on grid onxlabel('車速(km/h)'); ylabel('加速度倒數(shù)1/a'); legend('加速度倒數(shù)1/a曲線','Location','NorthWest'); title('加速度倒數(shù)曲線圖'); axis(0 120 0 10)% 加速時間曲線 k=length(n); for i=1:length(ig); t(i,1)=0; for j=2:k t(i,j)=abs(ua(i,j)-ua(i,j-1)*(b(i,j)+b(i,j-1)/2; endend for i=1:length(ig); for j=1:k at(i,j)=sum(t(i,1:j)/3.6; end endtotalat=at(1,:); for i=1:(length(ig)-1); for j=flag(i):k totalat=totalat totalat(length(totalat)+t(i+