break.rar

clear m0=4080;m=9290; hg0=0.845;hg=1.170; L=3.95; a0=2.1;a=2.95; Beta=0.38; g=9.8; b0=L-a0;b=L-a; G=m*g; %理想的前、后制动器制动力分配曲线 %subplot(2,2,1); figure for phi=0.1:0.1:1; Fu1=0:10:G; Fu2_1=phi*G-Fu1; plot(Fu1,Fu2_1,'r'); hold on Fu2_2=Fu1/((b+phi*hg)/(a-phi*hg)); ii=find(Fu2_1>Fu2_2); Fu2_2(ii)=NaN; plot(Fu1,Fu2_2,'b'); end Fu2=0.5*(G*sqrt(b^2+4*hg*L*Fu1/G)/hg-(G*b/hg+2*Fu1)); hold on plot(Fu1,Fu2,'g') axis([0,G+10000,0,G+10000]); xlabel('Fu1');ylabel('Fu2'); title('理想的前、后制动器制动力分配曲线') hold off %Beta线的绘制 G2=G/1000;G02=m0*g/1000; Fu12=0:1:G2; Fu22=Fu12*(1-Beta)/Beta; %subplot(2,2,2); figure plot(Fu12,Fu22); hold on Fu222=0.5*(G2*sqrt(b^2+4*hg*L*Fu12/G2)/hg-(G2*b/hg+2*Fu12)); plot(Fu12,Fu222); Fu223=0.5*(G02*sqrt(b0^2+4*hg0*L*Fu12/G02)/hg0-(G02*b0/hg0+2*Fu12)); plot(Fu12,Fu223); axis([0,20,0,30]); xlabel('Fu1/(kN)');ylabel('Fu2/(kN)'); title('β线与I曲线'); text(10,5,'I曲线（空载）');text(1,20,'I曲线（满载）');text(14,20,'B');text(10,15,'β线'); hold off %f线组和r线组 %subplot(2,2,3); figure FXb1=0:10:G*a/hg+10000; for phi3=0.2:0.2:1; FXb2r=-FXb1*phi3*hg/(L+phi3*hg)+phi3*G*a/(L+phi3*hg); plot(FXb1,FXb2r); hold on FXb2f=FXb1*(L-phi3*hg)/(phi3*hg)-G*b/hg; k=find(FXb2f>FXb2r); FXb2f(k)=NaN; plot(FXb1,FXb2f,'r'); end hold on FXb2=0.5*(G*(sqrt(b^2+4*hg*L*FXb1/G))/hg-(G*b/hg+2*FXb1)); plot(FXb1,FXb2,'g'); hold on y=0; plot(FXb1,y,'k'); axis([0,G*a/hg+10000,-G*b/hg-10000,100000]); xlabel('FXb1');ylabel('FXb2'); title('f线组与r线组') text(160000,-10000,'(G*a/hg,0)');text(10000,-G*b/hg,'(0,G*b/hg)'); %不同phi值路面上汽车制动过程分析图的绘制 %r线组和f线组 %subplot(2,2,4); figure for phi=0.1:0.1:0.7; Fu1=0:10:G2; Fu2_1=phi*G2-Fu1; plot(Fu1,Fu2_1,'r'); hold on FXb1=0:0.1:80; FXb2r=-FXb1*phi*hg/(L+phi*hg)+phi*G2*a/(L+phi*hg); FXb2f=FXb1*(L-phi*hg)/(phi*hg)-G2*b/hg; m=find(FXb2f>FXb2r);n=find(FXb2f<0); %裁剪不要的曲线 FXb2f(m)=NaN;FXb2f(n)=NaN;FXb2r(m)=NaN; plot(FXb1,FXb2f,'r'); hold on plot(FXb1,FXb2r); end %I曲线 Fu2=0.5*(G2*sqrt(b^2+4*hg*L*Fu1/G2)/hg-(G2*b/hg+2*Fu1)); hold on plot(Fu1,Fu2,'g') Fu223=0.5*(G02*sqrt(b0^2+4*hg0*L*Fu12/G02)/hg0-(G02*b0/hg0+2*Fu12)); hold on plot(Fu12,Fu223); %β线 Fu12=0:1:G2; Fu22=Fu12*(1-Beta)/Beta; hold on plot(Fu12,Fu22); axis([0,G2-10,0,G2-10]); xlabel('f线组 FXb1/kN Fu1/kN');ylabel('r线组 FXb2/kN Fu2/kN'); title('不同phi值路面上汽车制动过程分析') text(60,50,'I曲线(满载)');text(60,15,'I曲线（空载）');text(30,60,'β线'); %利用附着系数与制动强度的关系 figure z=0:0.01:1; phif0=Beta*z*L./(b0+z*hg0); phif=Beta*z*L./(b+z*hg); phir0=(1-Beta)*z*L./(a0-z*hg0); phir=(1-Beta)*z*L./(a-z*hg); ph=z; plot(z,phif0,'r*',z,phif,'r',z,phir0,'k*',z,phir,'k',z,ph,'b') axis([0,1,0,1]); text(0.8,0.4,'ψf（空车）'); text(0.8,0.7,'ψf（满载）'); text(0.4,0.9,'ψr（空车）'); text(0.8,0.95,'ψr（满载）'); text(0.8,0.8,'ψ=z'); xlabel('制动强度z/g');ylabel('利用附着系数ψ') title('利用附着系数与制动强度的关系'); %前、后制动效率 figure Ef=100*z./phif; Er0=100*z./phir0; Er=100*z./phir; kk=find(Ef>Er); jj=find(Ef<Er); Ef(kk)=NaN; Er(jj)=NaN; plot(phif,Ef,'r',phir0,Er0,'b',phir,Er,'k'); axis([0,1,0,105]); xlabel('附着系数ψ');ylabel('制动效率（%）'); title('前、后制动效率曲线'); text(0.6,60,'空载'); text(0.6,97,'满载'); text(0.5,77,'Er'); text(0.3,95,'Ef'); text(0.8,90,'Er');