《凸轮机构设计大作业》由会员分享,可在线阅读,更多相关《凸轮机构设计大作业(20页珍藏版)》请在金锄头文库上搜索。
1、大作业(二)凸轮机构设计(题号:8)班 级: 姓 名 、学 号: 成 绩: 完 成 日 期: 目录1. 凸轮机构大作业题目22. 推杆运动规律及凸轮廓线方程33. 程序流程图34. 源程序 55. 计算结果146. 凸轮机构图167. 体会及建议198. 参考资料20一、凸轮机构大作业题目试用计算机辅助设计完成下列摆动滚子推杆盘形凸轮机构的设计,已知数据如下表所示,凸轮沿着逆时针方向做匀速转动。表1 凸轮机构的推杆运动规律推程运动规律正弦加速度运动回程运动规律等加速等减速运动表2 两种凸轮机构的推杆在近休、推程、远休及回程阶段的凸轮转角题号近休凸轮转角推程凸轮转角远休凸轮转角回程凸轮转角A06
2、060180180270270360B04545210210260260360C03030210210280280360表3 摆动滚子推杆盘形凸轮机构的已知参数题号初选基圆半径r0/mm机架长度lOA/mm摆杆长度lAB/mm滚子半径rr/mm推杆摆角许用压力角许用最小曲率半径amin12A156055102435700.3rrB207065142640700.3rrC22727568182845650.35rr要求:每组(每三人为一组,每人一题)至少打印出一份源程序,每人打印出原始数据;凸轮理论轮廓和实际轮廓的坐标值;推程和回程的最大压力角,以及出现最大压力角时凸轮的相应转角;凸轮实际轮廓曲
3、线的最小曲率半径,以及相应的凸轮转角;和最后说确定的基圆半径。计算点数N=72120。绘出凸轮的理论轮廓和实际轮廓(可用计算机绘图)。二、推杆运动规律及凸轮廓线方程:推程(正弦加速度):s=h(/0)-sin(2/0)/(2) 回程(等加速段):s=h-2h2/02回程(等减速段):s=2h(0-)2/02凸轮理论廓线方程:x=lOAsin-lABsin(+0) y=lOAcos-lABcos(+0)式中,0为推杆的初始位置角,其值为: 0=arccos三、程序流程图|1?开始读入:r0,r0,rt,h或(),e或(lAB、lOA)1,2,3,4,1,2, amin,N计算:s0I=1计算:s
4、,x,y,ds/d,dx/d,dy/d,x,y计算:r0= r0+r0r0= r0=r0是回程?|2?选出1max及相应的凸轮转角1选出2max及相应的凸轮转角2计算:0?|-rtamin?计算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;%摆杆长度rr=18;%滚子半径fai=28*pi/180;%推杆摆角PI=3.141592653;alpha1=45;%许用压力角1alpha2=65
5、;%许用压力角阿尔法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=ones(1,121);Y=ones(1,121);XP=ones(1,121);YP=ones(1,121);FAI=ones(1,121);I=0;%初始
6、化while(I=N) F0=acos(a*a+L*L-r0*r0)/(2*a*L);delta=3*I*pi/180; if(delta=0) %正弦推程 F0=acos(a*a+L*L-r0*r0)/(2*a*L); 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*delta/delta1)/delta1); ddF=fai*2*pi*sin(
7、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; stheta=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; endendif(delta=180*pi/180&delta=250*pi/180&deltaalphamax2) alphamax2=alpha; deltamax2=delta*180/pi; end end i