滚子摆动从动件凸轮设计matlab程序.doc

disp * 滚子摆动从动件凸轮设计 * disp 已知条件: disp 凸轮作顺时针方向转动,从动件做摆动 disp 从动件在推程作等加速/等减速运动,在回程作等加速等减速运动 rb =52;rt = 10;qm=15;ft = 60;fs = 10;fh = 60;alp = 35;a=140;l=122;q0=asin(rb/a)*180/pi; fprintf (1, 基圆半径 rb = %3.4f mm n,rb) fprintf (1, 滚子半径 rt = %3.4f mm n,rt) fprintf (1, 起始角度 q0= %3.4f mm n,q0) fprintf (1, 最大摆动角度 qm = %3.4f mm n,qm) fprintf (1, 推程运动角 ft = %3.4f 度 n,ft) fprintf (1, 远休止角 fs = %3.4f 度 n,fs) fprintf (1, 回程运动角 fh = %3.4f 度 n,fh) fprintf (1, 推程许用压力角 alp = %3.4f 度 n,alp) hd= pi / 180;du = 180 / pi; %角度弧度互换 d1 = ft + fs;d2 = ft + fs + fh; disp disp 计算过程和输出结果: disp 1- 计算凸轮理论轮廓的压力角和曲率半径 disp 1-1 推程(等加速/等减速运动) s = zeros(ft);ds = zeros(ft);d2s = zeros(ft);vt=zeros(ft);st1=zeros(ft);at=zeros(ft); at = zeros(ft);atd = zeros(ft);pt = zeros(ft); for f = 1 : ft if f atm atm = atd(f); end end fprintf (1, 最大压力角 atm = %3.4f 度n,atm) for f = 1 : ft if abs(atd(f) - atm) alp fprintf (1, * 凸轮推程压力角超过许用值,需要增大基圆!n) end ptn = rb ;ftn=0; for f = 1 : ft if pt(f) ptn ptn = pt(f); end end fprintf (1, 轮廓最小曲率半径 ptn = %3.4f mmn,ptn) for f = 1 : ft if abs(pt(f) - ptn) 0.1 ftn = f;break end end fprintf (1, 对应的位置角 ftn = %3.4f 度n,ftn) if ptn rt + 5 fprintf (1, * 凸轮推程轮廓曲率半径小于许用值,需要增大基圆或减小滚子!n) end disp 1-2 回程(等加速等减速运动) s = zeros(fh);ds = zeros(fh);d2s = zeros(fh); ah = zeros(fh);ahd = zeros(fh);ph = zeros(fh); for f = d1 : d2 k = f - d1; if k atm atm = atd(f); end end fprintf (1, 最大压力角 atm = %3.4f 度n,atm) for f = 1 : ft if abs(atd(f) - atm) alp fprintf (1, * 凸轮推程压力角超过许用值,需要增大基圆!n) end ptn = rb ;ftn=0; for f = 1 : ft if pt(f) ptn ptn = pt(f); end end fprintf (1, 轮廓最小曲率半径 ptn = %3.4f mmn,ptn) for f = 1 : ft if abs(pt(f) - ptn) 0.1 ftn = f;break end end fprintf (1, 对应的位置角 ftn = %3.4f 度n,ftn) if ptn rt + 5 fprintf (1, * 凸轮推程轮廓曲率半径小于许用值,需要增大基圆或减小滚子!n) end disp 2- 计算凸轮理论廓线与实际廓线的直角坐标 n = 360; s = zeros(n);ds = zeros(n);r = zeros(n);rp = zeros(n); x = zeros(n);y = zeros(n);dx = zeros(n);dy = zeros(n); xx = zeros(n);yy = zeros(n);xp = zeros(n);yp = zeros(n); xxp = zeros(n);yyp = zeros(n); for f = 1 : n if f ft/2 & f ft & f d1 & f (d2-fh/2)&f d2 & f = n s = 0;ds = 0; end xx(f) = a*cos(f*hd)-l*cos(s+q0-f)*hd);x=xx(f); yy(f) = a*sin(f*hd)+l*sin(s+q0-f)*hd);y=yy(f); dx(f) = -a*sin(f*hd)+l*sin(s+q0-f)*hd)*(ds-1); dx = dx(f); dy(f) = a*cos(f*hd)+l*cos(s+q0-f)*hd)*(ds-1); dy = dy(f); xp(f) = x-rt*(dy/sqrt(dx2+dy2);xxp=xp(f); yp(f) = y+rt*(dx/sqrt(dx2+dy2);yyp = yp(f); r(f) = sqrt (x 2 + y 2 ); rp(f) = sqrt (xxp 2 + yyp 2 ); end disp 2-1 推程(等加速/等减速运动) disp 凸轮转角 理论x 理论y 实际x 实际y for f = 10 : 10 :ft nu = f xx(f) yy(f) xp(f) yp(f); disp(nu) end disp 2-2 回程(等加速/等减速运动) disp 凸轮转角 理论x 理论y 实际x 实际y for f = d1 : 10 : d2 nu = f xx(f) yy(f) xp(f) yp(f); disp(nu) end disp 2-3 凸轮轮廓向径 disp 凸轮转角 理论r 实际r for f = 10 : 10 : n nu = f r(f) rp(f); disp(nu) end disp 绘制凸轮的理论轮廓和实际轮廓: plot(xx,yy,r-.) % 理论轮廓(红色,点划线) axis (-(150) (150) -(150) (150) % 横轴和纵轴的下限和上限 axis equal % 横轴和纵轴的尺度比例相同 text(50,0,X) % 标注横轴 text(0,50,Y) % 标注纵轴 text(-5,5,O) % 标注直角坐标系原点 title(摆动从动件盘形凸轮设计) % 标注图形标题 hold on; % 保持图形 plot(-(rb