发布时间 : 星期三 文章内燃机课程设计6200柴油机曲轴设计动力计算更新完毕开始阅读ec13d356091c59eef8c75fbfc77da26925c59625
?r?0.923(rp/dp)=1.69；
(?0.2205?0.1015e/dp)?(7.8955?10.654(b/dp)?5.3482(b/dp)2?0.8570(b/dp)3)由计算结果可知，d=149.2mm＜150mm，故设计的曲轴可用。
. . .
附录 Matlab计算程序
>> %内燃机课程设计动力计算% a1 =0 : 5 : 720;%曲柄转角%
Pg1=[3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3.05,3.1,3.15,3.2,3.25,3.3,3.35,3.35,3.4,3.45,3.45,3.5,3.75,4,4.25,4.5,4.75,5,5.5,6.5,7,8,9,10.5,13,15,18,21.5,26,32,40,49,59,65,80,105,119,124,125,115,101,87.5,72.5,60,50,43,36.5,32,28,25,22,19.5,18,16.5,15.5,14,13,12.5,12,11,10.5,10.25,10,9.75,9.5,9,8.5,8,7.75,7,6.5,6,5.5,5,4,3.5,3,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,2.5,3];%气缸压力，kg/cm^2% a = 0 : 1 : 720;
Pg = interp1(a1,Pg1,a,'spline');
>> Pg = Pg/10.197;%气缸压力单位转化，Mpa% Ne = 440;%单位是kw%
mj = 35.76 + 0.347 * 34.76; %活塞组等效质量，kg% mb = 0.653 * 34.76; %连杆组算到大端的质量，kg% D = 200;%活塞直径，mm% L = 540;%连杆长度，mm% R = 135;%曲柄半径，mm% z = 6;%气缸数；
x = R/L;%曲柄连杆比%
B = asin(x*sin(a*pi/180));%连杆摆角% w = 600*pi/30;%转速，rad/s%
Pj = - mj * R * w^2 *(cos(a*pi/180) + x * cos(a*pi/90))/(pi * D^2/4 * 10^3);%往复惯性力，Mpa%
Pb = 0.1;%活塞底部气体压力，取为大气压力，Mpa% P = Pj - Pb + Pg;%总作用力，Mpa% figure(1);%打开新图版;
plot(a,Pg,a,Pj,'--',a,P,'-.');%蓝色的为气缸压力与曲轴转角的关系，黄色为往复惯性力与曲柄转角的关系,红色为总作用力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('(P,Pg,Pj)/Mpa');%加纵坐标% legend('Pg','Pj','P') grid on ;%添加网格%
Ph = P .* tan(B);%活塞侧推力，单位是Mpa% Pc = P./cos(B);%连杆力，单位是Mpa%
Pn = Pc .* cos(a*pi/180 + B);%法向力，单位是Mpa% Pt = Pc .* sin(a*pi/180 + B);%切向力，单位是Mpa% SumPt = Pt ;%为总切力,单位是Mpa;% for i=1:721
for j=1:5
m=i+720*j/z;
. . .
if m>721 m=m-720; end
SumPt(i)=SumPt(i)+Pt(m); end
end
avePt = mean(SumPt);%平均切向力,单位是Mpa% Rbh = Pt;%曲柄销负荷水平分量，单位是Mpa% Prb = mb * R * w^2/(pi * D^2/4 * 10^3);
Rbv = Pn - Prb; %曲柄销负荷垂直分量，单位是Mpa% Rb = (Rbh .* Rbh + Rbv .* Rbv).^0.5; %曲柄销总负荷% figure(2);%打开新图版%
plot(a,Ph);%画侧推力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('侧推力Ph/Mpa');%加纵坐标% grid on ;%添加网格% figure(3);%打开新图版%
plot(a,Pc);%画连杆力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('连杆力Pc/Mpa');%加纵坐标% grid on ;%添加网格%
figure(4);%打开新图版%
plot(a,Pn);%画法向力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('法向力Pn/Mpa');%加纵坐标% grid on ;%添加网格% figure(5);%打开新图版%
plot(a,Pt);%画切向力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('切向力Pt/Mpa');%加纵坐标% grid on ;%添加网格% figure(6);%打开新图版% A = 0 :1: 720;
plot(A,SumPt);%画总切向力与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('总切向力SumPt/Mpa');%加纵坐标% grid on ;%添加网格% figure(7);%打开新图版%
plot(a,Rb);%画曲柄销负荷与曲柄转角的关系% xlabel('曲柄转角a/deg');%加横坐标% ylabel('曲柄销负荷Rb/Mpa');%加纵坐标% grid on ;%添加网格%
Ni = avePt * pi * D^2 * R * w / (4 * 10^6)%由总切力计算指示功率% d = ( Ni - Ne/0.85 ) / Ni %计算误差%
. . .
. . .