凸轮机构大作业 机械原理 课设.doc

大作业（二）凸轮机构设计（题号： 7-ABC ）班 级_学 号_姓 名_成 绩_同 组 者_完成日期 2010 年 12 月 日目录凸轮机构设计21.题目及原始数据22推杆运动规律及图轮廓线方程23计算程序34计算结果及分析65凸轮廓线图76体会及建议87参考书89凸轮机构设计1.题目及原始数据使用计算机辅助设计完成题目要求的摆动滚子推杆判刑凸轮机构的设计，已知数据如下表，凸轮沿逆时针方向运动。表1转角转角进休转角推程转角远休转角回城转角xxA200950180900xxxB2001020180900zsC200980180900表2已知参数数据（mm）初选基圆半径机架长度摆杆长度滚子半径推杆摆角许用压力角需用最小曲率半径xxA156055102435700.3xxxB207065142640700.3zsC227268182845650.352推杆运动规律及图轮廓线方程表3运动规律题号推程运动规律回程运动规律7A,B,C等加速等减速运动余弦加速度运动摆杆转角方程理论廓线方程工作廓线方程3计算程序%本段代码属621原创，版权所有，如有雷同纯属巧合，禁止复制粘贴。%此代码最终解释权归621所有clear %初始化数据q=1; %选择组别(1,2,3-A,B,C)r0=15,20,22; loa=60,70,72;lab=55,65,68; rr=10,14,18;h=24,26,28; alpha1=5,40,45;alpha2=70,70,65; u=0.3,0.3,0.35;j1=60,45,30; j2=180,210,210;j3=270,260,280;r0=r0(q);loa=loa(q);lab=lab(q);rr=rr(q);h=h(q);alpha1=alpha1(q);alpha2=alpha2(q);u=u(q);j1=j1(q);j2=j2(q);j3=j3(q);j=0:1:360; for i=1:100 s0=acosd(loa2+lab2-r02)/(2*loa*lab); alpha=0; for p=1:361 if j(p)j1 s(p)=0; dxdj(p)=loa*cosd(j(p)-lab*cosd(j(p)+s0); dydj(p)=-loa*sind(j(p)+lab*sind(j(p)+s0); ddxdj(p)=-loa*sind(j(p)+lab*sind(j(p)+s0); ddydj(p)=-loa*cosd(j(p)+lab*cosd(j(p)+s0); alpha(p)=atand(lab*0-(loa*cosd(s0+s(p)-lab)/loa*sind(s0+s(p); else if j(p)0.5*(j2+j1) s(p)=2*h*(j(p)-j1)2/(j2-j1)2); dxdj(p)=loa*cosd(j(p)-lab*cosd(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*(1+4*h*(j(p)-j1)/(j2-j1)2); dydj(p)=-loa*sind(j(p)+lab*sind(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*(1+4*h*(j(p)-j1)/(j2-j1)2); ddxdj(p)=-loa*sind(j(p)+lab*sind(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*(1+4*h*(j(p)-j1)/(j2-j1)2)2-4*lab*cosd(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*h/(j2-j1)2; ddydj(p)=-loa*cosd(j(p)+lab*cosd(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*(1+4*h*(j(p)-j1)/(j2-j1)2)2+4*lab*sind(j(p)+2*h*(j(p)-j1)2/(j2-j1)2+s0)*h/(j2-j1)2; alpha(p)=atand(lab*pi/180*4*h*(j(p)-j1)/(j2-j1)2-(loa*cosd(s0+s(p)-lab)/loa*sind(s0+s(p); else if j(p)j2 s(p)=h-2*h*(j2-j(p)2/(j2-j1)2); dxdj(p)=loa*cosd(j(p)+lab*cosd(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*(-1-4*h*(j2-j(p)/(j2-j1)2); dydj(p)=-loa*sind(j(p)+lab*sind(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*(-1-4*h*(j2-j(p)/(j2-j1)2); ddxdj(p)=-loa*sind(j(p)-lab*sind(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*(-1-4*h*(j2-j(p)/(j2-j1)2)2+4*lab*cosd(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*h/(j2-j1)2; ddydj(p)=-loa*cosd(j(p)+lab*cosd(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*(-1-4*h*(j2-j(p)/(j2-j1)2)2+4*lab*sind(-j(p)-h+2*h*(j2-j(p)2/(j2-j1)2-s0)*h/(j2-j1)2; alpha(p)=atand(lab*(-4*h*(j2-j(p)/(j2-j1)2)-(loa*cosd(s0+s(p)-lab)/loa*sind(s0+s(p); else if j(p)0 row(i)=10000; else row(i)=-row(i); end end alpha=abs(alpha); alphat=abs(alphat); if (max(alphat)alpha1)|(max(alpha)alpha2)|(min(row)-rr)u*rr) r0=r0+1; else break; end breakendmax_t,j_t=max(alphat);a=max_t,j_t;max_h,j_h=max(alpha);a=a,max_h,j_h;min_row,j_row=min(row);min_row=min_row-rr;a=a,min_row,j_row,r0sinj=dxdj./sqrt(dxdj.2+dydj.2);cosj=-dydj./sqrt(dxdj.2+dydj.2);x=loa*sind(j)-lab.*sind(j+s+s0);y=loa*cosd(j)-lab.*cosd(j+s+s0);x1=x-rr*cosj;y1=y-rr*sinj; % 绘图hold ongrid onaxis equalplot(x,y,LineWidth,2,LineStyle,-,color,b)plot(x1,y1,LineWidth,2,LineStyle,-,color,r)plot(0,0,.,MarkerSize,20,color,k)ylabel(Y)xlabel(X)title(理论，实际轮廓线)legend(理论轮廓线,实际轮廓线)4计算结果及分析表4结果数据（mm）推程最大压力角对应角度回程最大压力角对应角度最小曲率半径对应角度基圆半径xxA8.709712043.18843224.92608016xxxB