汽车理论编程题.docx
- 文档编号:29409650
- 上传时间:2023-07-23
- 格式:DOCX
- 页数:14
- 大小:17.65KB
汽车理论编程题.docx
《汽车理论编程题.docx》由会员分享,可在线阅读,更多相关《汽车理论编程题.docx(14页珍藏版)》请在冰豆网上搜索。
汽车理论编程题
n=600:
1:
4000;r=;i0=;eff=;f=;m=3880;g=;G=m*g;CdA=;a=;hg=;L=;Iw1=;Iw2=;Iw=Iw1+Iw2;If=;
forig=[];
Ua=*r*n/ig/i0;
Ft=Ttq*ig*i0*eff/r;
Ff=m*g*f;
Fw=CdA*power(Ua,2)/;
Fs=Ff+Fw;
plot(Ua,Ft,Ua,Fs);
holdon;
end
text(15,13000,'F_t_1');
text(20,7000,'F_t_2');
text(40,4000,'F_t_3');
text(55,2500,'F_t_4');
text(70,1900,'F_t_5');
text(105,2100,'F_f+F_w');
xlabel('u_akm/h');
ylabel('F_tN');
title('驱动力-行驶阻力平衡图');
[x,y]=ginput
(1);disp('汽车的最高车速');disp(x);disp('km/h');
Ua=*r*n/ig/i0;
Ft=Ttq**i0*eff/r;
Ff=m*g*f;
Fw=CdA*power(Ua,2)/;
Fs=Ff+Fw;
i=atan((Ft-Fs)/G);
disp('汽车的最大爬坡度');disp(max(i));
C=max(i)/(a/L+hg/L*max(i));
disp('克服最大爬坡度时的附着率');disp(C);
figure;
forig=[];
Ua=*r*n/ig/i0;
Ft=Ttq*ig*i0*eff/r;
q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2);
Ff=m*g*f;
Fw=CdA*power(Ua,2)/;
Fs=Ff+Fw;
Ac=(Ft-Fs)/q/m;
plot(Ua,1./Ac);
holdon;
end
axis([098010]);
text(5,,'1');
text(10,1,'2');
text(20,,'3');
text(30,,'4');
text(40,4,'5');
title('行驶加速度倒数曲线');
xlabel('Ua/(km/h)');
ylabel('1/a');
figure;
u
(1)=*r*600/i0/;dt=;w=1;t
(1)=0;ig=;
while(u(w)<
n(w)=u(w)*ig*i0/r;
Ttq(w)=+*n(w)/*(n(w)/1000)^2+*(n(w)/1000)^*(n(w)/1000)^4;
Ft(w)=Ttq(w)*ig*i0*eff/r;
Ff=G*f;
Fw(w)=CdA*u(w)^2/;
q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2);
a(w)=*(Ft(w)-(Ff+Fw(w)))/q/m;
u(w+1)=u(w)+a(w)*dt;
if(u(w+1)>=*4000*r/i0)
ig=;
end
if(u(w+1)>=*4000*r/i0)
ig=;
end
if(u(w+1)>=*4000*r/i0)
ig=;
end
t(w+1)=(w+1)*dt;
w=w+1;
end
plot(t,u);
holdon
u=70;
plot(t,u);
axis([0600100]);
title('2挡原地起步加速时间曲线');
xlabel('时间t/s'),ylabel('车速u/(km/h)');
[x1,y1]=ginput
(1);disp('加速到70km/h的时间');disp(x1);
n=600:
1:
4000;r=;i0=;eff=;f=;m=3880;g=;G=m*g;CdA=;a=;hg=;L=;Iw1=;Iw2=;Iw=Iw1+Iw2;If=;
forig=[,,,,]
Ua=*r*n/ig/i0;
Pe=Ttq.*n/9550;
plot(Ua,Pe);
holdon;
end
Pmax=max(Pe);
plot(Ua,Pmax);
Ua=0:
:
max(Ua);
Pf=G*f*Ua/3600;
Pw=CdA*Ua.^3/76140;
plot(Ua,(Pf+Pw)/eff);
text(10,50,'1');
text(15,35,'2');
text(30,40,'3');
text(40,30,'4');
text(70,40,'5');
text(110,80,'(P_f+P_w)/η');
xlabel('u_akm/h');
ylabel('F_tN');
title('驱动力-行驶阻力平衡图');
title('汽车的功率平衡图'),xlabel('Ua/(km/h)'),ylabel('P/kw');
figure;
n1=[815,1207,1614,2012,2603,3006,3403,3804];
b0=[,,,,,,,];
b1=[,,,,,,,];
b2=[,,,,,,,];
b3=[,,,,,,,];
b4=[,,,,,,,];
B0=spline(n1,b0,n);
B1=spline(n1,b1,n);
B2=spline(n1,b2,n);
B3=spline(n1,b3,n);
B4=spline(n1,b4,n);
pg=7;
forig=[,]
Ua=*r*n/ig/i0;
Pf=G*f*Ua/3600;
Pw=CdA*Ua.^3/76140;
P=(Pf+Pw)/eff;
b=B0+B1.*P+B2.*P.^2+B3.*P.^3+B4.*P.^4;
Q=P.*b./*Ua*pg);
plot(Ua,Q);
holdon;
end
axis([01001228]);
text(40,14,'四挡');
text(80,16,'五档');
title('最高挡与次高挡等速百公里油耗曲线'),xlabel('Ua/(km/h)'),ylabel('百公里油耗Qs/(L/100km)');
3
主函数:
i0=[,,,,];
fori=1:
1:
5
y(i)=jiasushijian(i0(i));
end
y;
fori=1:
1:
5
b(i)=youhao(i0(i));
end
b;
plot(b,y,'+r');
holdon;
b1=linspace(b
(1),b(5),100);
y1=spline(b,y,b1);
plot(b1,y1);
text,,'i_0=');
text,,'i_0=');
text,,'i_0=');
text,,'i_0=');
text,,'i_0=');
title('燃油经济性—加速时间曲线');
xlabel('百公里油耗(L/100km)');
ylabel('加速时间s');
子函数1
functiony=jiasushijian(i0)
n1=linspace(0,5000);
nmax=4000;nmin=600;r=;yita=;CDA=;f=;G=(3880)*;ig=[,,,];
fori=1:
1:
4
uamax(i)=chesu(nmax,r,ig(i),i0);
uamin(i)=chesu(nmin,r,ig(i),i0);
ua(i,:
)=linspace(uamin(i),uamax(i),100);
n(i,:
)=zhuansu(ua(i,:
),r,ig(i),i0);
Ttq(i,:
)=zhuanju(n(i,:
));
Ft(i,:
)=qudongli(Ttq(i,:
),ig(i),i0,yita,r);
F(i,:
)=f*G+CDA*(ua(i,:
).^2)/;
delta(i,:
)=1+++*(ig(i)^2)*(i0^2)*yita)/(3880*r^2);
a(i,:
)=1./(delta(i,:
).*3880./(Ft(i,:
)-F(i,:
)));
F2(i,:
)=Ft(i,:
)-F(i,:
);
end
temp1(1,:
)=ua(2,:
)/;
temp1(2,:
)=1./a(2,:
);
n1=1;
forj1=1:
1:
100
ifua(3,j1)>max(ua(2,:
))&&ua(3,j1)<=70
temp2(1,n1)=ua(3,j1)/;
temp2(2,n1)=1./a(3,j1);
n1=n1+1;
end
end
n2=1;
forj1=1:
1:
100
ifua(4,j1)>max(ua(3,:
))&&ua(4,j1)<=70;
temp3(1,n2)=ua(4,j1)/;
temp3(2,n2)=1./a(4,j1);
n2=n2+1;
end
end
y=temp1(1,1)*temp1(2,1)+qiuji(temp1(1,:
),temp1(2,:
))+qiuji(temp2(1,:
),temp2(2,:
))+qiuji(temp3(1,:
),temp3(2,:
));
end
子函数2
function[ua]=chesu(n,r,ig,i0)
ua=*r*n/i0./ig;
end
子函数3
function[q]=jiasu(umax,umin,ig,a,ua0,i0)
globalfGCDAyitamrIfIw1Iw2pgB0B1B2B3B4n;
ua1=umin:
1:
umax;
delta=1+(Iw1+Iw2)/(m*r^2)+(If*ig^2*i0^2*yita)/(m*r^2);
P0=(G*f.*ua0./3600+CDA.*ua0.^3/76140+(delta*m.*ua0/3600)*a)/yita;
P=(G*f.*ua1/3600+CDA.*ua1.^3/76140+(delta*m.*ua1/3600)*a)/yita;
dt=1/*a)
fori=1:
1:
3401
b0(i)=B0(i)+B1(i)*P0(i)+B2(i)*(P0(i))^2+B3(i)*(P0(i))^3+B4(i)*(P0(i))^4;
end
b1=interp1(ua0,b0,ua1);
Qt=P.*b1./.*pg);
i1=size(Qt);
i=i1
(2);
Qt1=Qt(2:
i-1);
q=(Qt
(1)+Qt(i))*dt./2+sum(Qt1)*dt;
end
子函数4
function[y]=zhuanju(n)
y=+.*(n./1000).*(n./1000).^2+.*(n./1000).^.*(n./1000).^4;
end
子函数5
function[n]=zhuansu(ua,r,ig,i0)
n=ig*i0.*ua./*r);
end
子函数6
function[b]=youhao(i0)
globalfGCDAyitamrIfIw1Iw2pgB0B1B2B3B4n
ig=[,,,];r=;
yita=;CDA=;f=;%i0=;
G=(3880)*;If=;Iw1=;Iw2=;m=3880;
n0=[8151207161420122603300634033804];
B00=[];
B10=[];
B20=[];
B30=[];
B40=[];
n=600:
1:
4000;
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);
ua4=*r.*n./(i0*ig(4));
F4=f*G+CDA*(ua4.^2)/;
P_fw4=F4.*ua4./(yita**1000);
fori=1:
1:
3401
b4(i)=B0(i)+B1(i)*P_fw4(i)+B2(i)*(P_fw4(i))^2+B3(i)*(P_fw4(i))^3+B4(i)*(P_fw4(i))^4;
end
pg=;
ua4_m=[25,40,50];
s_m=[50,250,250];
b4_m=spline(ua4,b4,ua4_m);
F4_m=f*G+CDA*(ua4_m.^2)/;
P_fw4_m=F4_m.*ua4_m./(yita**1000);
Q4_m=P_fw4_m.*b4_m.*s_m./(102.*ua4_m.*pg);
Q4_a1=jiasu(40,25,ig(4),,ua4,i0);
Q4_a2=jiasu(50,40,ig(4),,ua4,i0);
Qid=;tid=;s=1075;
Q_i=Qid*tid;
Q4all=(sum(Q4_m)+Q4_a1+Q4_a2+Q_i)*100/s;
b=Q4all;
end
子函数7
function[p]=qiuji(x0,y0)
n0=size(x0);n=n0
(2);
x=linspace(x0
(1),x0(n),200);
y=spline(x0,y0,x);
plot(x,y);
p=trapz(x,y);
end
beta=;L=;a1=;b1=L-a1;hg1=;a2=;b2=L-a2;hg2=;g=;z=0:
:
1;
faif1=beta*z./(1/L*(b1+z*hg1));
fair1=(1-beta)*z./(1/L*(a1-z*hg1));
faif2=beta*z./(1/L*(b2+z*hg2));
fair2=(1-beta)*z./(1/L*(a2-z*hg2));
plot(z,faif1,'r',z,fair1,'g',z,faif2,'b',z,fair2,'k');
text,,'φ_r(空车)');
text,,'φ_r(满载)');
text,1,'φ=z');
text,,'φ_f(满载)');
text,,'φ_f(空车)');
xlabel('制动强度z/g');
ylabel('利用附着系数φ');
title('利用附着系数与制动强度的关系曲线');
fai=0:
:
1;
Ef1=(b1/L)./(beta-fai*hg1/L);
Er1=(a1/L)./((1-beta)+fai*hg1/L);
Ef2=(b2/L)./(beta-fai*hg2/L);
Er2=(a2/L)./((1-beta)+fai*hg2/L);
figure;
plot(fai,Ef1,'r',fai,Er1,'g',fai,Ef2,'b',fai,Er2,'k');
text,,'E_f');
text,,'满载');
text,,'E_r');
text,,'E_r');
text,,'空载');
xlabel('附着系数φ');
ylabel('制动效率');
title('前、后制动效率曲线');
axis([0101]);
fai=;
Er1=(a1/L)/((1-beta)+fai*hg1/L);
Er2=(a2/L)/((1-beta)+fai*hg2/L);
abmax1=fai*Er1*g;
abmax2=fai*Er2*g;
t21=;
t22=;
u0=30;
s1=1/*(t21+t22/2)*u0+u0^2/abmax1;
s2=1/*(t21+t22/2)*u0+u0^2/abmax2;
disp('空载时的制动距离'),disp(s1);disp('满载时的制动距离'),disp(s2);
z1=a1*fai/(L+fai*hg1);
z2=a2*fai/(L+fai*hg2);
abmax1=z1*g;
abmax2=z2*g;
s1=1/*(t21+t22/2)*u0+u0^2/abmax1;
s2=1/*(t21+t22/2)*u0+u0^2/abmax2;
disp('前部管路损坏时空载制动距离');disp(s1);disp('前部管路损坏时满载制动距离');disp(s2);
z1=b1*fai/(L-fai*hg1);z2=b2*fai/(L-fai*hg2);abmax1=z1*g;abmax2=z2*g;
s1=1/*(t21+t22/2)*u0+u0^2/abmax1;
s2=1/*(t21+t22/2)*u0+u0^2/abmax2;
disp('后部管路损坏时空载制动距离');disp(s1);disp('后部管路损坏时满载制动距离');disp(s2);
m=;Iz=3885;L=;a=;b=;k1=-62618;k2=-110185;i=20;g=;R0=15;u1=;K=m*(a/k2-b/k1)/L^2;
Uch=(1/K)^(1/2);
disp('稳定性因数(s^2/m^2)K=');disp(K);
disp('特征车速(m/s)Uch=');disp(Uch);u=0:
:
30;
S=u./(L*(1+K*u.^2));
plot(u,S);
title('汽车稳态横摆角速度增益曲线');xlabel('车速u(m/s)');
ylabel('稳态横摆角速度增益');
disp('u=s时,转向灵敏度为');disp(S(448));
SM=k2/(k1+k2)-a/L;ay=*g;A=K*ay*L;B=L/R0;R=L/(B-A);
C=R/R0;
disp('静态储备系数.=');disp(SM);
disp('侧向加速度为时前、后轮侧偏角绝对值之差(rad)?
a1-a2=');disp(A);
disp('侧向加速度为时转弯半径比值R/R0=');disp(C);
W0=L/u1*(k1*k2/(m*Iz)*(1+K*u1^2))^(1/2);
D=(-m*(k1*a^2+k2*b^2)-Iz*(k1+k2))/(2*L*(m*Iz*k1*k2*(1+K*u1^2))^(1/2));
t=atan((1-D^2)^(1/2)/(-m*u1*a*W0/(L*k2)-D))/(W0*(1-D^2)^(1/2));
E=atan((1-D^2)^(1/2)/D)/(W0*(1-D^2)^(1/2))+t;
disp('车速u=s时的瞬态响应参数分别为:
');disp('横摆角速度波动的固有(圆)频率(rad)为?
');disp(W0);
disp('阻尼比为');disp(D);
disp('反应时间(s)为');disp(t);
disp('峰值反应时间(s)为');disp(E);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 汽车 理论 编程