哈工大机械设计电算实验(matlab2010).doc

机械设计电算实验一：普通V带传动设计内容和任务1、普通V带传动设计内容给定原始数据：传递的功率P，小带轮转速n，传动比i及工作条件。设计内容：带型号，基准长度Ld，根数Z，传动中心距a，带轮基准直径dd1、dd2，带轮轮缘宽度B，初拉力F0，和压轴力FQ。2、CAD任务： （1）编制V带传动设计程序框图。（2）编制V带传动设计原程序。（3）按习题或作业中数据运行程序，要求对每一组数据各按三种V带型号计算，对每一种带型号选三种小带轮直径进行计算并输出所有结果。二、变量标识符序号名称代号程序符号单位1名义功率Ppkw2设计功率Pdpdkw3工作情况系数KAkka4小带轮转速n1nr/min5工作机载荷特性载荷类型6原动机类型原动机类型7工作年限a年b班8带型号V带带型9小带轮基准直径dd1dd1mm10大带论基准直径dd2dd2mm11带速vvm/s12最大带速vmaxvmaxm/s13预估中心距a0a0mm14基准长度Ldldmm15实际中心距aaamm16小带轮包角1r度17V带根数Zz18基本额定功率P0pokw19功率增量P0pp0kw20传动比系数Kiki21长度系数KLkl22包角系数Kka23弯曲影响系数Kbkb24初拉力F0f0N25压轴力QfqN26轮缘宽度Bbbmm 三、程序框图开始输入原始数据P,n,i并选择电机计算设计功率由带型号确定最大带速，小带轮直径系，带的相关参数m，f，e 相关系数Kb确定带型号选择小带轮直径验算带速vvmax?计算大带轮直径给定范围，预估中心距AA计算Ld并取标准值，计算实际中心距并圆整计算小带轮包角计算选择相关系数：K， Ki 。插值计算P0 ，计算P0计算带根数ZZ=600r/min的内燃机)disp(II类原动机包括交流电动机(双笼型、滑环式、单相、大转差率)，);disp( 直流电动机（复激、串激），单缸发动机，n=600r/min的内燃机)d1=input(请选择原动机的类型，输入1或2。请输入:);disp(载荷类型：);disp(1、载荷平稳；2、载荷变动小；3、载荷变动较大；4、载荷变动很大。);e1=input(请选择载荷类型 请输入:);c=3 710 38 80;2.2 940 28 60;3 960 38 80;4 960 38 80;if p=3&n=710 d=38; e=80; disp(根据您的数据为您选择的电动机为 Y132M-8);elseif p=2.2&n=940 d=28; e=60; disp(根据您的数据为您选择的电动机为 Y112M-6);elseif p=3&n=960 d=38; e=80; disp(根据您的数据为您选择的电动机为 Y132S-6);elseif p=4&n=960 d=38; e=80; disp(根据您的数据为您选择的电动机为 Y132M1-8);else disp(欢迎继续使用本程序);end ka2=1.0 1.1 1.2;1.1 1.2 1.3;1.2 1.3 1.4;1.3 1.4 1.5;ka3=1.1 1.2 1.3;1.2 1.3 1.4;1.4 1.5 1.6;1.5 1.6 1.8;if d1=1 kka=1.1*k(e1,b);else kka=k(e1,b);endpd=kka*p;t2=1;while t2=1z=10;while z=10disp(V带带型：1、Z型:dd1min=50mm；2、A型:dd1min=75mm；3、B型:dd1min=125mm);x=input(请输入（1、2、3）来选择带型。 请选择：);x1=0.246 0.449 0.794;7.44 19.02 50.6;0.441 0.765 1.31;0.2925 0.7725 1.9875;k1=x1(1,x);k2=x1(2,x);k3=(10(-4)*x1(3,x);kb=10(-3)*x1(4,x);vmax=25;disp(优选带轮直径系列： 50 56 63 71 80 90 100)disp( 112 125 140 150 160 180 200)disp( 224 250 280 315 355 400 425)disp( 450 500 560 600 630 710 800)dd1=input(请参考dd1min以及优选系列，输入小带轮基准直径(mm) dd1=);v=pi*dd1*n/60/1000;while v=vmax disp(请减小小带轮基准直径dd1); dd1=input(请参考dd1min,重新输入小带轮基准直径(mm) dd1=); v=pi*dd1*n/60/1000;enddd2=round(i*dd1);a1=0.7*(dd1+dd2);a2=2*(dd1+dd2);disp(a1=);disp(a1);disp(a2=);disp(a2);a0=input(已知中心距a1=a0=a2,请输入a0的值（尽量为整数）。 a0=);while a0a2 disp(wrong!您输入的数值不在a1-a2之间，请重新输入); a0=input(请重新输入a0的值（尽量为整数）。 a0=);endld=2*a0+pi/2*(dd1+dd2)+(dd1-dd2)2/4/a0;kk=400 0.87 0 0;450 0.89 0 0;500 0.91 0 0;560 0.94 0 0;630 0.96. 0.81 0;710 0.9 0.83 0;800 1.00 0.85 0.82;900 1.03 0.87 0.84;. 1000 1.06 0.89 0.86;1120 1.08 0.91 0.88;1250 1.10 0.93 0.90;. 1400 1.14 0.96 0.92;1600 1.16 0.99 0.95;1800 1.18 1.01 0.98;. 2000 0 1.03 1.00;2240 0 1.06 1.03;2500 0 1.09 1.05;2800 0 . 1.11 1.07;3150 0 1.13 1.09;3550 0 1.17 1.13;4000 0 1.19 1.15;. 4500 0 0 1.18;for g=1:22 if ld=kk(g,1)&ld=(kk(g+1,1)-ld) ld1=kk(g+1,1); else ld1=kk(g,1); end break endendkl=kk(g,x+1);aa=round(a0+(ld-ld1)/2);r=180-(dd2-dd1)/aa*57.3;rr=220 1.20;210 1.15;200 1.10;190 1.05;180 1.00;170 0.98;160 0.95;. 150 0.92;140 0.89;130 0.86;120 0.82;110 0.78;100 0.73;90 0.68;for ii=1:13 if r=rr(ii+1,1) ka=rr(ii+1,2)+(r-rr(ii+1,1)/10*(rr(ii,2)-rr(ii+1,2); break endendp0=(k1*v(-0.09)-k2/dd1-k3*v2)*v;if i=1.8 ki=1.1202;elseif i=2 ki=1.1372;endm=0.06 0.1 0.17;pp0=kb*n*(1-1/ki);z=ceil(pd/(p0+pp0)/kl/ka);if z=vmax disp(请减小小带轮基准直径dd1); dd1=input(请参考dd1min,重新输入小带轮基准直径(mm) dd1=); v=pi*dd1*n/60/1000;enddd2=round(i*dd1);a1=0.7*(dd1+dd2);a2=2*(dd1+dd2);disp(a1=);disp(a1);disp(a2=);disp(a2);a0=input(已知中心距a1=a0=a2,请输入a0的值（尽量为整数）。 a0=);while a0a2 disp(wrong!您输入的数值不在a1-a2之间，请重新输入); a0=input(请重新输入a0的值（尽量为整数）。 a0=);endld=2*a0+pi/2*(dd1+dd2)+(dd1-dd2)2/4/a0;kk=400 0.87 0 0;450 0.89 0 0;500 0.91 0 0;560 0.94 0 0;630 0.96. 0.81 0;710 0.9 0.83 0;800 1.00 0.85 0.82;900 1.03 0.87 0.84;. 1000 1.06 0.89 0.86;1120 1.08 0.91 0.88;1250 1.10 0.93 0.90;. 1400 1.14 0.96 0.92;1600 1.16 0.99 0.95;1800 1.18 1.01 0.98;. 2000 0 1.03 1.00;2240 0 1.06 1.03;2500 0 1.09 1.05;2800 0 . 1.11 1.07;3150 0 1.13 1.09;3550 0 1.17 1.13;4000 0 1.19 1.15;. 4500 0 0 1.18;for g=1:22 if ld=kk(g,1)&ld=(kk(g+1,1)-ld) ld1=kk(g+1,1); else ld1=kk(g,1); end break endendkl=kk(g,x+1);aa=round(a0+(ld-ld1)/2);r=180-(dd2-dd1)/aa*57.3;rr=220 1.20;210 1.15;200 1.10;190 1.05;180 1.00;170 0.98;160 0.95;. 150 0.92;140 0.89;130 0.86;120 0.82;110 0.78;100 0.73;90 0.68;for ii=1:13 if r=rr(ii+1,1) ka=rr(ii+1,2)+(r-rr(ii+1,1)/10*(rr(ii,2)-rr(ii+1,2); break endendp0=(k1*v(-0.09)-k2/dd1-k3*v2)*v;if i=1.8 ki=1.1202;elseif i=2 ki=1.1372;endm=0.06 0.1 0.17;pp0=kb*n*(1-1/ki);z=ceil(pd/(p0+pp0)/kl/ka);if z10 f0=500*pd/z/v*(2.5-ka)/ka+m(x)*v2; fq=2*z*f0*sin(r/2/180*pi); bb=round(z-1)*15+20); disp(所选原动机种类：) if d1=1 disp(I类原动机); else disp(II类原动机); end disp(所选带型为：); if x=1 disp(Z型V带); elseif x=2 disp(A型V带); else disp(B型V带); end disp(小带轮直径(mm) dd1=) disp(dd1); disp(大带轮直径(mm) dd2=) disp(dd2); disp(实际中心距(mm) a0=); disp(aa); disp(小带轮包角（角度） =); disp(r); disp(V带根数 Z=); disp(z); disp(V带的基准长度(mm) Ld=); disp(ld); disp(初拉力(N) F0=); disp(f0); disp(作用在轴上的压力(N) FQ=); disp(fq); disp(带轮轮缘宽度(mm) B=); disp(bb);endt1=input(是否计算下一种带轮直径？ 1、yes；2、no。 请选择：);endt2=input(是否计算下一型号带？ 1、yes；2、no。 请选择：);if t2=2 disp(欢迎下次使用本程序)endendclear all;五、运算示例5.1.5方案方案电动机工作功率Pd/kW电动机满载转速nm/(r/min)工作机的转速nw/(r/min)第一级传动比轴承座中心高H/mm最短工作年限工作环境5.1.537108021703年3班室内、清洁运算结果截图：若选择