附程序：
%copyright gejianyong
clc
clear
close all;
g=9.8
ma=9290%满载质量
m0=4080%空载质量
Ga=ma*g%满载重力
G0=m0*g%空载重力
hga=1.17%满载质心高度
hg0=0.6%空载质心高度
L=3.95%轴距
ba=1%满载质心至后轴距离
b0=1.85%空载质心至后轴距离
aa=2.95%满载质心至前轴距离
a0=2.1%空载质心至前轴距离
B=0.38%制动力分配系数
% f1前轮制动器制动力
% f2a满载后轮理想制动器制动力
%以下为满载时制动过程
f1=0:10:60000;
f2a=0.5*(Ga*((ba*ba+4*hga*L*f1/Ga).^0.5)/hga-(Ga*ba/hga+2*f1));%满载I曲线公式
%f2Ba满载后轮实际制动器制动力
f2Ba=f1*(1-B)/B;%满载B线
figure(1)
plot(f1/1000,f2a/1000,'k',f1/1000,f2Ba/1000,'k')%画出I曲线，B线
%P附着系数
for P=0.1:0.1:1
fxbfa=(L-P*hga)*f1/P/hga-Ga*ba/hga;%fxbfa满载f线
fxbfa1=fxbfa(fxbfa<=f2a);%取I曲线下方f线
f1f=f1(fxbfa<=f2a);
fxbra=-P*hga*f1/(L+P*hga)+P*Ga*aa/(L+P*hga);%fxbra满载r线
fxbra1=fxbra(fxbra>=f2a);%取I曲线上方r线
f1r=f1(fxbra>=f2a);
hold on
plot(f1f/1000,fxbfa1/1000,'k',f1r/1000,fxbra1/1000,'k')%画出f线
axis([0 60 0 60])
%axis square
end
%title('满载时不同φ值路面的制动过程分析')
xlabel('{\itf} 线组 {\itF}_{μ1}/kN,{\itF}_{Xb1}/kN')
ylabel('{\itr} 线组 {\itF}_{μ2}/kN,{\itF}_{Xb2}/kN')
%以下为空载时制动过程
f1=0:10:30000;
f20=0.5*(G0*((b0*b0+4*hg0*L*f1/G0).^0.5)/hg0-(G0*b0/hg0+2*f1));%空载I曲线公式
%f2B0空载后轮实际制动器制动力
f2B0=f1*(1-B)/B;%空载B线
figure(2)
plot(f1/1000,f20/1000,'k',f1/1000,f2B0/1000,'k')%画出I曲线，B线
%P附着系数
for P=0.1:0.1:1
fxbf0=(L-P*hg0)*f1/P/hg0-G0*b0/hg0;%fxbf0空载f线
fxbf01=fxbf0(fxbf0<=f20);%取I曲线下方f线
f1f=f1(fxbf0<=f20);
fxbr0=-P*hg0*f1/(L+P*hg0)+P*G0*a0/(L+P*hg0);%fxbr0空载r线
fxbr01=fxbr0(fxbr0>=f20);%取I曲线上方r线
f1r=f1(fxbr0>=f20);
hold on
plot(f1f/1000,fxbf01/1000,'k',f1r/1000,fxbr01/1000,'k')%画出f线
axis([0 30 0 30])
%axis square
end
%title('空载时不同φ值路面的制动过程分析')
xlabel('{\itf} 线组 {\itF}_{μ1}/kN,{\itF}_{Xb1}/kN')
ylabel('{\itr} 线组 {\itF}_{μ2}/kN,{\itF}_{Xb2}/kN')
%以下为利用附着系数与制动强度的关系
z=0.01:0.01:1;
%z=0.2:0.2:1%计算数据用
Pfa=B*z*L./(ba+z*hga);%满载前轴利用附着系数
Pra=(1-B)*z*L./(aa-z*hga);%满载后轴利用附着系数
Pf0=B*z*L./(b0+z*hg0);%空载前轴利用附着系数
Pr0=(1-B)*z*L./(a0-z*hg0);%空载后轴利用附着系数
Pz=z;%理想利用附着系数
Pl=(z+0.07)/0.85;%法规
Pll=Pl(0.2<=Pl&Pl<=0.8);
zl=z(0.2<=Pl&Pl<=0.8);
figure(3)
plot(z,Pfa,'k',z,Pra,'k',z,Pf0,'k--',z,Pr0,'k--',z,Pz,'k--','LineWidth',1.5)
hold on
plot(zl,Pll,'k')
fplot('[z-0.08,z+0.08]',[0.15,0.3],'k')
fplot('(z-0.02)/0.74',[0.3,1],'k')
axis([0 1 0 1])
%title('利用附着系数与制动强度的关系曲
线')
xlabel('制动强度{\itz/g}')
ylabel('利用附着系数{\itφ}')
%以下为制动效率与附着系数的关系曲线
P=0:0.01:1;
%P=0.2:0.2:1%计算数据用
Ef=ba./L./(B-P*hga./L);
Er=aa./L./((1-B)+P*hga./L);
Er0=a0./L./((1-B)+P*hg0./L);
figure(4)
plot(P,Ef*100,P,Er*100,P,Er0*100,'color',[0 0 0])
axis([0 1 0 100])
%title('前、后制动效率曲线')
xlabel('附着系数{\itφ}')
ylabel('制动效率（%）')
%以下为评价
P=0.8%同步附着系数为0.8
P0a=(L*B-ba)/hga%满载同步附着系数
P00=(L*B-b0)/hg0%空载同步附着系数
%计算知后轮先抱死
v=60%正常行驶国标制动初速度
sl=0.15*v+v*v/130%正常行驶国标制动距离
vb=50%失效行驶国标制动初速度
slba=0.15*vb+100*vb*vb/30/115%失效行驶满载国标制动距离
slb0=0.15*vb+100*vb*vb/25/115%失效行驶空载国标制动距离
za=P*aa/(L*(1-B)+P*hga)%满载制动强度
aamax=za*g%满载最大制动减速度
z0=P*a0/(L*(1-B)+P*hg0)%空载制动强度
a0max=z0*g%空载最大制动减速度
sa=(0.02+0.2/2)*v/3.6+v*v/25.92/aamax%满载正常行驶制动距离计算公式
s0=(0.02+0.2/2)*v/3.6+v*v/25.92/a0max%空载正常行驶制动距离计算公式
B=0%前管路损坏 后轮先抱死
zaf=P*aa/(L*(1-B)+P*hga)%满载制动强度
aafmax=zaf*g%满载最大制动减速度
z0f=P*a0/(L*(1-B)+P*hg0)%空载制动强度
a0fmax=z0f*g%空载最大制动减速度
saf=(0.02+0.2/2)*v/3.6+v*v/25.92/aafmax%满载失效行驶制动距离计算公式
s0f=(0.02+0.2/2)*v/3.6+v*v/25.92/a0fmax%空载失效行驶制动距离计算公式
B=1%后管路损坏 前轮先抱死
zar=P*ba/(L*B-P*hga)%满载制动强度
aarmax=zar*g%满载最大制动减速度
z0r=P*b0/(L*B-P*hg0)%空载制动强度
a0rmax=z0r*g%空载最大制动减速度
sar=(0.02+0.2/2)*v/3.6+v*v/25.92/aarmax%满载失效行驶制动距离计算公式
s0r=(0.02+0.2/2)*v/3.6+v*v/25.92/a0rmax%空载失效行驶制动距离计算公式