車輛工程基于MATLAB的動力性仿真分析及優化設計程序_第1頁
車輛工程基于MATLAB的動力性仿真分析及優化設計程序_第2頁
車輛工程基于MATLAB的動力性仿真分析及優化設計程序_第3頁
車輛工程基于MATLAB的動力性仿真分析及優化設計程序_第4頁
車輛工程基于MATLAB的動力性仿真分析及優化設計程序_第5頁
已閱讀5頁，還剩8頁未讀， 繼續免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、n=linspace(600,4000,100);%均分計算指令 ,600 最低轉速 ,4000 最高轉速 ,均分為 100 等分r=0.367;i0=5.83;nt=0.85;G=3880*9.8;f=0.013;CDA=2.77;If=0.218;Iw1=1.798;Iw2=3.598;m=3880;L=3.2;a=1.947;hg=0.9; ig=6.09,3.09,1.71,1.00;%輸入已知參數ua1=0.377*r*n/i0/ig(1);ua2=0.377*r*n/i0/ig(2);ua3=0.377*r*n/i0/ig(3); ua4=0.377*r*n/i0/ig(4);%各

2、轉速各擋位下的速度Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;%從 6004000rpm 油擬合公式計算發動機轉距Ft1=Tq*i0*ig(1)*nt/r;Ft2=Tq*i0*ig(2)*nt/r;Ft3=Tq*i0*ig(3)*nt/r;Ft4=Tq*i0*ig(4)*nt/r;%從 6004000rpm 各擋位的驅動力Ff=G*f; ua=linspace(0,200,100);Fw=CDA*ua.*ua/ 21.15; %空氣阻力 plot(ua1,Ft1,ua2,Ft2

3、,ua3,Ft3,ua4,Ft4,ua,Ff+Fw);%畫出各擋位的 Ua-Ft,及 Ua-Ff+Ftxlabel(ua/ km/h);ylabel(F/N);%標注橫縱軸title( 汽車驅動力 -行駛阻力平衡圖 ); %標注圖形題目 gtext(Ft1),gtext(Ft2),gtext(Ft3),gtext(Ft4),gtext(Ff+Fw);%給每根線條添加符號legend(Ft1,Ft2,Ft3,Ft4,Ff+Fw);%標注圖例umax=max(ua4); disp( 汽車最高車速 =);disp(umax);disp(km/h); imax=tan(asin(max(Ft1-(Ff

4、+Fw)/G); %最大爬坡度的公式 disp( 汽車最大爬坡度 =);disp(imax);%輸出最高車速 ,與最大爬坡度的結果n=600:1:4000;%600最低轉速 ,4000 最高轉速 ,相鄰數組間隔 1r=0.367;i0=5.83; eff=0.85; f=0.013;CdA=2.77;m=3880;g=9.8;%輸入已知參數G=m*g;Ttq=-19.313+295.27*n/1000-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4; %從 6004000rpm 油擬合公式計算發動機轉距for ig=6.09,3.09

5、,1.71,1.00Ua=0.377*r*n/ig/i0;%各轉速各擋位下的速度Pe=Ttq.*n/9550;%各轉速下的功率plot(Ua,Pe);hold on; %使當前軸及圖形保持而不被刷新 ，準備接受此后將繪制的圖形 ，多圖共存 endUa=0:0.1:max(Ua);Pf=G*f*Ua/ 3600;%滾動阻力Pw=CdA*Ua.3/76140; %空氣阻力 plot(Ua,(Pf+Pw)/eff);title( 汽車的功率平衡圖 ), xlabel(Ua/(km/h),ylabel(P/kw);%畫出汽車的功率平衡圖gtext(Ft1),gtext(Ft2),gtext(Ft3),

6、gtext(Ft4),gtext(Pf+Pw)/nt);legend( , , , ,Pf+Pw/nt);n=600:1:4000; %600最低轉速 ,4000最高轉速 ,相鄰數組間隔 r=0.367;i0=5.83;nt=0.85;f=0.013;CdA=2.77;m=3880;g=9.8; %輸入已知參數G=m*g;Ttq=-19.313+295.27*n/1000-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4; %從 6004000rpm 油擬合公式計算發動機轉距for ig=6.09,3.09,1.71,1.00 Ua=0

7、.377*r*n/ig/i0;Ft=Ttq*i0*ig*nt/r;Fw=CdA*Ua.2/21.15;D=(Ft-Fw)/G%汽車動力因子公式plot(Ua,D);%畫出汽車動力特性圖hold on; %使當前軸及圖形保持而不被刷新，準備接受此后將繪制的圖形，多圖共存 endf=0.0076+0.000056*Ua %滾動阻力與速度之間的關系 plot(Ua,f);%畫出速度與滾動阻力圖title( 汽車動力特性圖 ), %給圖加題目 xlabel(Ua/(km/h),ylabel(D);gtext(),gtext(),gtext(),gtext( ),gtext(f);legend( , ,

8、 , ,f);n=600:10:4000; %600最低轉速 ,4000最高轉速 ,相鄰數組間隔 10 m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=6,09 3.09 1.71 1.00;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;%輸入已知參數Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4; %從 6004

9、000rpm 油擬合公式計算發動機轉距Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r; %各轉速各擋位下的驅動力ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;%各擋位各轉速下的速度Fw1=CDA*ua1.2/21.15;Fw2=CDA*ua2.2/21.15;Fw3=CDA*ua3.2/21.15;Fw4=CDA*ua4.2/21.15;

10、%不同速度下的空氣阻力Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r2)+(If*ig(1)2*i02*nT)/(m*r2);deta2=1+(Iw1+Iw2)/(m*r2)+(If*ig(2)2*i02*nT)/(m*r2);deta3=1+(Iw1+Iw2)/(m*r2)+(If*ig(3)2*i02*nT)/(m*r2);%不同擋位下的汽車旋轉質量換deta4=1+(Iw1+Iw2)/(m*r2)+(If*ig(4)2*i02*nT)/(m*r2); 算系數a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m

11、);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4; %各擋位下的加速度plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4);title( 汽車的加速度倒數曲線 );xlabel(ua(km/h); ylabel(1 /a ); %作汽車加速度倒數曲線 gtext(1/a1),gtext(1 /a2),gtext(1 /a3),gtext(1 /a4); legend(1/a1,1/a2,1/a3,1/a4);n=600:10:4000;m=3880;g=9

12、.8;nmin=600;nmax=4000;G=m*g;ig=6.09 3.09 1.71 1.00;nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; %輸入已知參數Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4; %從 6004000rpm 油擬合公式計算發動機轉距Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;

13、Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;Fw1=CDA*ua1.2/21.15;Fw2=CDA*ua2.2/21.15;Fw3=CDA*ua3.2/21.15;Fw4=CDA*ua4.2/21.15;Ff=G*f;i1=asin(Ft1-Ff-Fw1)/G);i2=asin(Ft2-Ff-Fw2)/G);i3=asin(Ft3-Ff-Fw3)/G);i4=asin(Ft

14、4-Ff-Fw4)/G);%各轉速各擋位下的驅動力%各擋位各轉速下的速度%不同速度下的空氣阻力%不同檔位下的坡度plot(ua1,i1,ua2,i2,ua3,i3,ua4,i4);title( 汽車的爬坡度圖 );xlabel(ua/(km*h-1);ylabel(i/%);%作汽車的坡度圖gtext(),gtext(),gtext(),gtext( );m=3880;g=9.8;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;pg=7.1;%汽油的重度取 7.1N/Lig=6.09 3.09 1.71 1;n=600:1:4000;n0=815 1207 1

15、614 2012 2603 3006 3403 3804;B00=1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7;B10=-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291;B20=72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113;B30=-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215;B40

16、=0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568;%輸入已知參數B0=spline(n0,B00,n);B1=spline(n0,B10,n);B2=spline(n0,B20,n);B3=spline(n0,B30,n);B4=spline(n0,B40,n); %使用三次樣條插值，保證曲線的光滑連續 ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;%求出發動機轉速范圍內對應的 3、4 檔車速Pe3=(m*g*f*ua3 /3600+CdA*ua3.3/76140)/0.85;Pe4=(m*g*f*ua4 /3600+CdA*ua4.3/76140)/0.85; %發動機功率 for i=1:1:3401%用擬合公式求出各個燃油消耗率b3(i)=B0(i)+B1(i)*Pe3(i)+B2(i)*Pe3(i).2+B3(i)*Pe3(i).3+B4(i)*Pe3(i).4; b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*Pe4(i).2+B3(i)*Pe4(i).3+B4(i)*Pe4(i).4;%插值得出對

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論