车辆工程基于MATLAB的动力性仿真分析和优化设计程序文件.docVIP

word文档整理分享 参考资料 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); %各转速各挡位下的速度 Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; %从600~4000rpm油拟合公式计算发动机转距 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; %从600~4000rpm各挡位的驱动力 Ff=G*f; ua=linspace(0,200,100); Fw=CDA*ua.*ua/21.15; %空气阻力 plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua,Ff+Fw); %画出各挡位的Ua-Ft,及Ua-Ff+Ft xlabel(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+Fw))/G))); %最大爬坡度的公式 disp(汽车最大爬坡度=); disp(imax); %输出最高车速,与最大爬坡度的结果 n=600:1:4000; %600最低转速,4000最高转速,相邻数组间隔1 r=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; %从600~4000rpm油拟合公式计算发动机转距 for ig=[6.09,3.09,1.71,1.00] Ua=0.377*r*n/ig/i0; %各转速各挡位下的速度 Pe=Ttq.*n/9550; %各转速下的功率 plot(Ua,Pe); hold on; %使当前轴及图形保持而不被刷新，准备接受此后将绘制的图形，多图共存 end Ua=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),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; %从600~4000rpm油拟合公式计算发动机转距 for ig=[6.09,3.09,1.71,1.00] Ua=0.