1、远休凸轮转角回程凸轮转角A0 6060 180180 270270 360B0 4545 210210 260260 360C0 3030 210210 280280 360表 3 摆动滚子推杆盘形凸轮机构的已知参数初选基圆机架长度摆杆长度滚子半径推杆摆角许用压力角许用最小曲率半径 r 0/mml OA/mmlAB/mmr r/mm1 2半径 amin 156055102435700.3r r70652640C 22 7275 68 18 28 45 65 0.35r r要求:每组(每三人为一组,每人一题)至少打印出一份源程序,每人打印出原始数据; 凸轮理论轮廓和实际轮廓的坐标值; 推程和回程
2、的最大压力角, 以及出现最大压力角时凸轮的相应转角; 凸轮实际轮廓曲线的最小曲率半径, 以及相应的凸轮转角;和最后说确定的基圆半径。计算点数 N=72120。绘出凸轮的理论轮廓和实际轮廓(可用计算机绘图) 。二、推杆运动规律及凸轮廓线方程 :推程(正弦加速度):s=h(/)-sin(2/)/(2)0 02 2回程(等加速段):s=h-2h/ 回程(等减速段):s=2h(-) / 凸轮理论廓线方程: x=lOAsin-lABsin(+0)y=lOAcos -lABcos(+0)式中, 为推杆的初始位置角,其值为 :0=arccos (a2l 2r 0) / 2(al )三、程序流程图开始读入:
3、r0, r0,rt ,h 或( ),e 或( lAB、 lOA)1,2, 3,4, 1 ,2, amin,N计算: s0I=1 s,x,y,ds/d , dx/d ,dy/d , x, yr0= r0+ r0r0= r0= r0 | |是回程?2?选出及相应的凸轮转角 2max1max1 0?| |-r tamin?计算 a 选出 | amin| 及相应的凸轮转角 aminI=I+1IN?打印: x,y,x, y,amin,amin,1max,1max,2max,2max, r0, s四、源程序clear;r0=22;%初选的基圆半径dr0=0.05;a=72; %机架长度L=68;%摆杆长度
4、rr=18;%滚子半径fai=28*pi/180;% 推杆摆角PI=3.141592653;alpha1=45;%许用压力角 1alpha2=65;%许用压力角阿尔法 2lambda=6.3;%许用最小曲率半径N=120;%取用点的个数delta1=180*pi/180;% 推程凸轮最大转角delta2=70*pi/180; %远休凸轮最大转角delta3=80*pi/180;% 回程凸轮最大转角delta4=30*pi/180;% 近休凸轮最大转角alphamax1=0;% 推程最大压力角初值alphamax2=0; %回程最大压力角初值roumin=100; %凸轮最小曲率半径初值X=on
5、es(1,121);Y=ones(1,121);XP=ones(1,121);YP=ones(1,121);FAI=ones(1,121);I=0;%初始化while(I=N)F0=acos(a*a+L*L-r0*r0)/(2*a*L);delta=3*I*pi/180;if(delta=0) % 正弦推程F=fai* (delta/delta1)-sin(2*pi*delta/delta1)/(2*pi); x=a*sin(delta)-L*sin(delta+F+F0);y= a*cos(delta)-L*cos(delta+F+F0);dF= fai*(1/delta1-cos(2*pi
6、*delta/delta1)/delta1); ddF=fai*2*pi*sin(2*pi*delta/delta1)/(delta1*delta1);dx=a*cos(delta)-L*cos(delta+F+F0)*(1+dF);ddx=-a*sin (delta)+L* sin(delta+F+F0)*(1+dF)2-L*cos(delta+F+F0)*ddF;dy=-a*sin(delta)+L*sin(delta+F+F0)*(1+dF);ddy=-a*cos (delta)+L* cos(delta+F+F0)*(1+dF)2+ L*sin(delta+F+F0)*ddF;sthe
7、ta=dx/sqrt(dx*dx+dy*dy);ctheta=-dy/sqrt(dx*dx+dy*dy);xp=x-rr*ctheta;yp=y-rr*stheta;alpha=atan(L*abs(dF)+(a*cos(F+F0)-L)/(a*sin(F+F0)*180/pi;alpha=abs(alpha);if (alphaalphamax1)alphamax1=alpha;deltamax1=delta*180/pi;endif(delta=180*pi/180&deltaalpha1)|(alphamax2alpha2) %优化条件 r0=r0+dr0;continuerou=(dx
8、2+dy2)(3/2)/(dx*ddy-dy*ddx); % 计算曲率半径if (rou0)rou=-rou;if (rou-rr)(0.35*rr) % 优化条件r0=r0+dr0;if(rouroumin)roumin=rou;deltamin=delta *180/pi;roumina=roumin-rr;X(I+1)=x;Y(I+1)=y;XP(I+1)=xp;YP(I+1)=yp;I=I+1;figure(1);axis equalhold ont=0:3:360;X_1=r0*cosd(t); %画基圆Y_1=r0*sind(t);X_2=rr*cosd(t)+X(10); %画滚
9、子Y_2=rr*sind(t)+Y(10);plot(X_1,Y_1,m-,X,Y,:,XP,YP,k,X_2,Y_2,c-);legend(基圆 ,理论轮廓 实际轮廓 plot(0,0,ko %固定凸轮点plot(X(10),Y(10), %固定滚子点title( 凸轮轮廓曲线图 xlabel(X/mmylabel(Y/mmhold offdisp(推程最大压力角: alphamax1推程最大压力角相应转角:deltamax1回程最大压力角:alphamax2回程最大压力角相应转角:deltamax2凸轮最小曲率半径:roumin凸轮最小曲率半径相应转角:roumina最后确定的凸轮基圆半径
10、:r0五、计算结果A 组:alphamax1 =34.9492deltamax1 =45alphamax2 =46.7626deltamax2 = 255roumin =15.0000roumina =5.0000r0 =19.7000B 组:alphamax1 =24.2568deltamax1 = 177.0000alphamax2 = 51.9666deltamax2 = 291.0000roumin = 24.3000roumina =6.3000r0 =32.8000C组:六、凸轮机构图C 组:七、体会及建议本次凸轮机构的设计,我们熟悉了解析法在机构设计中的应用,锻炼了编程能力。这次
11、大作业用到了非常多的 MATLAB编程知识,我们查阅了大量资料,也请教了不少同学,最后编出来程序实属不易。通过这次亲手设计实践, 我们也认识到具体的操作远远比想象中的难,尤其是工程问题, 涉及到方方面面的知识, 这一点在以后的学习生活中也应该引起注意,注重细节才能做得更好。另外,大作业带给我们最大的好处就是我们更加熟悉凸轮的运动规律, 这比起上课老师的讲解印象更深,因为加入了自己的思考和亲手设计, 我们也能更加熟练地理解凸轮的设计实践,我想这是大作业带给我们最大的好处。这次凸轮设计的大作业也引起了我们的思考, 因为这实在是涉及到太多的编程知识, C 语言、 MATLAB等都可能用到,然而我们平常的教学过程中没有怎么教到这些知识,尤其是 MATLAB,我们想学校
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1