凸轮机构运动分析.docx

1、凸轮机构运动分析大作业（二）凸轮机构设计（题号： 6-ABC摆动滚子推杆盘型凸轮机构）班 级_05020905_成 绩_小组成员周海龙2009301399 王江 2009301464 王建辉2009301395_完成日期 2011 年 11 月 27日目录凸轮机构设计 31.题目及原始数据 32推杆运动规律及图轮廓线方程 33计算程序 44计算结果及分析 85凸轮廓线图 86体会及建议 107参考书 108小组互评 10凸轮机构设计1.题目及原始数据使用计算机辅助设计完成题目要求的摆动滚子推杆判刑凸轮机构的设计，已知数据如下表，凸轮沿逆时针方向运动。表1 转角转角进休转角推程转角远休转角回程转

2、角王江06060180180270270360周海龙04545210210260260360王建辉03030210210280280360表2 已知参数数据（mm）初选基圆半径机架长度摆杆长度滚子半径推杆摆角许用压力角需用最小曲率半径王江156055102435700.3周海龙207065142640700.3王建辉227268182845650.352推杆运动规律及图轮廓线方程表3 运动规律题号推程运动规律回程运动规律6A,B,C等加速等减速运动余弦加速度运动摆杆转角方程：理论廓线方程：工作廓线方程：3计算程序clearq=1; %选择组别(1,2,3-A,B,C)r0=15,20,22;

3、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-r0

4、2)/(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、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*(

6、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)/lo

7、a*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

8、-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

9、)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

10、_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,-,co