第一章作业。
车辆0802班孙铭辰学号：0120807250216
function fun()
m1=2000;m2=1800;m3=3880;g=9.81;a=1.947;l=3.2;hg=0.9;
r=0.367;nt=0.85;f=0.013;cda=2.77;
i0=5.83;if=0.218;iw1=1.798;iw2=3.589;
ig=[5.56 2.769 1.644 1.00 0.793];
iw=iw1*2+iw2*4输入已知参数。
for i=1:137
n(i)=(i+23)*25;
tq(i)=-19.313+295.27*(n(i)/1000)-165.
44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;
end从600~4000rpm油拟合公式计算发动机转距。
for j=1:5
for i=1:137
ft(i,j)=tq(i)*ig(j)*i0*nt/r;
ua(i,j)=0.377*r*n(i)/ig(j)/i0;
fz(i,j)=cda*ua(i,j)^2/21.15+f*m3*g; %计算各档位下的驱动力和行驶阻力。
endend
plot(ua,ft,ua,fz);
title('汽车驱动力—行驶助力图');
xlabel('ua(km/h)')
ylabel('ft(n作汽车驱动力与行驶阻力平衡图。
text(15,13000,'ft1');
text(26,6500,'ft2');
text(44,3800,'ft3');
text(75,23000,'ft4');
text(115,1550,'ft5');
text(110,2510,'ff+fw');
for i=1:200
n(i)=3300+i*0.1;
tq(i)=-19.313+295.27*(n(i)/1000)-165.
44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;
ft(i)=tq(i)*ig(5)*i0*nt/r;
ua(i)=0.377*r*n(i)/ig(5)/i0;
fz(i)=cda*ua(i)^2/21.15+f*m3*g;
e(i)=abs((ft(i)-fz(i)))
end从3300rpm开始以0.1rpm为步长搜索直到找到。
for i=1:200驱动力与行驶阻力之差绝对值最小点(ⅴ档处)
if(e(i)==min(e通过其他检验程序可估计出该点大致位置(3310rpm附近)
disp('最高车速输出结果。
disp(ua(i));
disp('km/h');
endend
for i=1:250
n(i)=1800+i;
tq(i)=-19.313+295.27*(n(i)/1000)-165.
44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;
ft(i)=tq(i)*ig(1)*i0*nt/r;
ua(i)=0.377*r*n(i)/ig(1)/i0;
fz(i)=cda*ua(i)^2/21.15+f*m3*g;
b(i)=asin((ft(i)-fz(i))/m3*g));
k(i)=tan(b(i));
co(i)=k(i)/(a/l+hg*k(i)/l); 爬坡度对应的附着率。
end从1800rpm开始以1rpm为步长收索直到找到。
for i=1:250
if(b(i)==max(b最大爬坡度点。
disp('最大爬坡度输出结果。
disp(k(i)*100
disp('%
disp('最大爬坡度相应的附着率=')
disp(co(i));
endend
for i=1:137
n(i)=(i+23)*25;
tq(i)=-19.313+295.27*(n(i)/1000)-165.
44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;
end从600~4000rpm油拟合公式计算发动机转距。
for j=1:5
for i=1:137
ft(i,j)=tq(i)*ig(j)*i0*nt/r;
ua(i,j)=0.377*r*n(i)/ig(j)/i0;
fz(i,j)=cda*ua(i,j)^2/21.15+f*m3*g; %计算各档位下的驱动力和行驶阻力。
q(j)=1+iw/(m3*r^2)+if*ig(j)^2*i0^2*nt/(m3*r^2);
a0(i,j)= ft(i,j)- fz(i,j))/q(j)*m3);
a0d(i,j)=1/a0(i,j);
endend
figure
for k=1:4由于第5档是超速档，一般不会用它去爬坡，另外在第五档会出现加速度很小的点，造成倒数很大，故只画出前四档。
plot(ua(:,k),a0d(:,k
hold on
endtitle('加速度倒数曲线');
xlabel('ua(km/h)')
ylabel('1/a0');
text(3.842,0.5,'1/a1');
text(28,0.7,'1/a2');
text(52,1.5,'1/a3');
text(77,4,'1/a4');
for i=1:137计算出2档换为3档，3档换为4档，以及到达70km/h速度时对应的转速。
m=ua(i,3)-ua(137,2);
if(m>0&m<1)
j=i; end
endfor i=1:137
n=ua(i,4)-ua(137,3);
if(n>0&n<1)
k=i;end
endfor i=1:137
if(ua(i,4)<=70)
m=i;end
endt=ua(1,2)*a0d(1,2汽车以2档起步加速到70km/h所用的时间。
for i=2:137
t1=(ua(i,2)-ua(i-1,2))*a0d(i,2)+a0d(i-1,2))/2;
t=t+t1;
endfor i=j:137
t2=(ua(i,3)-ua(i-1,3))*a0d(i,3)+a0d(i-1,3))/2;
t=t+t2;
endfor i=k:m
t3=(ua(i,4)-ua(i-1,4))*a0d(i,4)+a0d(i-1,4))/2;
t=t+t3;
endt=t+(ua(j,3)-ua(137,2))*a0d(137,2)+(ua(k,4)-ua(137,3))*a0d(137,3)+(70-ua(m,4))*a0d(m,4); 加上换挡处未加上的时间。
ts=t/3.6换为秒，因为1m/s=1km/h*3.6
disp('加速时间=')
disp(ts);
disp('s');
运行后：出现的图形。
运行结果为：最高车速=
km/h最大爬坡度=
最大爬坡度相应的附着率=
加速时间=
s第二章作业。
function fun()
m1=2000;m2=1800;m3=3880;g=9.81;a=1.947;l=3.2;hg=0.9;
r=0.367;nt=0.85;f=0.013;cda=2.77;
i0=5.83;if=0.218;iw1=1.798;iw2=3.589;
ig=[5.56 2.769 1.644 1.00 0.793];
iw=iw1*2+iw2*4;
n0=[815 1207 1614 2012 2603 3006 3403 3804];
for i=1:137
n(i)=(i+23)*25;
tq(i)=-19.313+295.27*(n(i)/1000)-165.
44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;
end从600~4000rpm油拟合公式计算发动机转距。
for j=1:5
for i=1:137
汽车理论作业
第一章。1 什么是汽车的动力性？其评价指标是什么？2 何谓汽车驱动力 行驶阻力平衡图，试以此图分析汽车动力性。3 什么是动力因数d，何谓汽车的动力特性图，试以此图分析汽车动力性。4 何谓汽车功率平衡方程式，请用结构参数写出该式。5 画出五挡汽车功率平衡图并利用该图分析汽车的动力性？计算ua6 何谓汽...
汽车理论作业
姓名 陈洪平。学号 021100701 班级 车辆一班。1.3解 1 绘制驱动力与行驶阻力平衡图的matlab程序 i0 5.83 r 0.367 nt 0.85 m 3880 g 9.8 f 0.013 cda 2.27 ig 6.09 3.09 1.71 1 n 600 1 4000 ua 0 ...
汽车理论作业
武汉理工大学。课程名称 汽车理论。开课学院 汽车工程学院。指导老师 邓亚东。学生姓名 王彬城。专业班级 车辆0902班。学生学号 0120907250220 第一章汽车的动力性。本题选用五档变速器。clcm1 2000 m2 1800 m3 3880 g 9.81 r 0.367 nt 0.85 f...