




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領
文檔簡介
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+
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 教育科技融合啟航盛典見證新高度
- 全球鈾礦資源分布與核能產(chǎn)業(yè)技術(shù)創(chuàng)新與產(chǎn)業(yè)升級路徑研究報告
- 公交優(yōu)先戰(zhàn)略2025年城市交通擁堵治理的公共交通與城市社區(qū)建設協(xié)同報告
- ChEMBL22003-生命科學試劑-MCE
- 浙江工業(yè)職業(yè)技術(shù)學院《核醫(yī)學基礎》2023-2024學年第一學期期末試卷
- 甘肅省武威五中學2025屆化學九上期末檢測模擬試題含解析
- 廣東信息工程職業(yè)學院《行為矯正學》2023-2024學年第一學期期末試卷
- 2024-2025學年江蘇省南京市鼓樓區(qū)鼓樓實驗中學七年級數(shù)學第一學期期末經(jīng)典模擬試題含解析
- 江西師范高等專科學?!度诤峡臻g設計》2023-2024學年第一學期期末試卷
- 廣州華南商貿(mào)職業(yè)學院《概率論與數(shù)理統(tǒng)計》2023-2024學年第一學期期末試卷
- 重力式混凝土擋土墻施工方案
- 2024年盤扣式腳手架安裝服務合同一
- PCR實驗室(新冠核酸檢測實驗室)SOP文件 (一)
- 委托代辦購買水果合同范例
- 2024至2030年輕鋼隔墻龍骨項目投資價值分析報告
- 透析器首次使用綜合癥的應急預案
- 河南開放大學本科《公司法律實務(本)》形考任務1-3+終考試題及答案
- 舊房加裝電梯基礎施工方案
- 國開2024年《JavaScript程序設計》形考1-3答案
- 2024年中國沖擊波醫(yī)療器械市場調(diào)查研究報告
- 廣東省四校2023-2024學年高二下學期期末聯(lián)考+英語試卷(含答案)
評論
0/150
提交評論