汽車動力性計算matlab程序

1、精品好資料學習推薦% 汽車動力性計算（自己編的動力性計算程序，供大家計算動力性時參考，具體參數(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 1829

2、 1748 1669 1700 1524 1105; figure(1) plot(n_test,T_test,g); hold ongrid 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)矩曲線); legend(測試曲線,擬合曲線); %所給車型動力總成相關(guān)參數(shù) ig=3.07 2.16 1.48 1.0 0.82; i0

3、=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 ongrid onfor i=length(n); Pe=Ttq.*n/9550; endAX,H1,H2=plotyy(n,Ttq,n,Pe); xlabel(發(fā)動機轉(zhuǎn)速n(r/min); ylabel(發(fā)動機轉(zhuǎn)矩Ttq(N*m); ylabel(AX(2),發(fā)動機功率Pe(Kw); title(發(fā)動機外特性曲線); % 各擋位速度曲線

4、 %計算各擋位車速 for i=1:length(ig); ua(i,:)=0.377*r*n/ig(i)/i0; end%計算各檔位最高車速 uamax=ua(:,length(ua(1,:); figure(3) hold onfor i=1:length(ig); plot(n,ua(i,:),k); endhold ongrid onxlabel(轉(zhuǎn)速n(r/min); ylabel(各擋位車速(km/h); title(各擋位車速-轉(zhuǎn)速表); legend(1擋車速,2擋車速,3擋車速,4擋車速,5擋車速); % 驅(qū)動力和行駛阻力平衡圖 %計算滾動阻力系數(shù) for i=1:lengt

5、h(ig); f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100).4; end%計算滾動阻力 for i=1:length(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;endfigure(4) hold

6、 onfor i=1:length(ig); plot(ua(i,:), Ft(i,:),k); plot(ua(i,:), F(i,:),r); plot(ua(i,:), Ff(i,:),b); endhold ongrid onxlabel(車速(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);end%計算風阻阻力功率 for i=1

7、: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,:);endfigure(5) hold onfor i=1:length(ig); plot(ua(i,:), P(i,:),k); plot(ua(i,:), Pz(i,:),r); endhold ongrid onxlabel(車速(km/h

8、); ylabel(發(fā)動機功率、阻力功率(kW); legend(發(fā)動機功率P,阻力功率Pz,Location,NorthWest); title(功率平衡圖); % 動力特性圖（動力因數(shù)圖） for i=1:length(ig); D(i,:)= (Ft(i,:)- Fw(i,:)/M/g; endfigure(6) hold onfor i=1:length(ig); plot(ua(i,:), D(i,:),k); plot(ua(i,:), f(i,:),r); endhold ongrid onxlabel(車速(km/h); ylabel(動力因數(shù)D); legend(動力因數(shù)D,

9、滾動阻力系數(shù)f); title(動力特性圖); % 爬坡度曲線圖 for i=1:length(ig); I(i,:)= (tan(asin(Ft(i,:)-(Ff(i,:)+Fw(i,:)/(M*g)*100; endfigure(7) hold onfor i=1:length(ig); if i=1 plot(ua(i,:),I(i,:),r); elseplot(ua(i,:),I(i,:),k); endendhold ongrid onxlabel(車速(km/h); ylabel(最大爬坡度（%）); legend(擋,高速檔); title(爬坡度曲線圖); % 加速度曲線圖

10、deta=1+1/M*4*Iw/r2+1/M*If*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)=ua(i,length(n) flag(i)=j; break; endend ad1=ad1 b(i+1,j:length(n); ad2=ad2 ua(i+1,j:length(n); endfigure(10) hold onplot(ad2,ad1,k); hold ongrid onxlabel(車速(km

11、/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; endendfor i=1:length(ig); for j=1:k at(i,j)=sum(t(i,1:j)/3.6; endendtotalat=at(1,:); for i=1:(length(ig)-1); for j=flag(i):k totalat=totalat totalat(length(totalat)