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

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 = ft / 2 s(f)=2*(qm/ft^2)*f^2;st1(f)=s(f);s = s(f); %推程加速方程式 ds(f)=(qm/ft^2)*f;vt(f)=ds(f);ds = ds(f); d2s(f)=4*qm/ft;at(f)=d2s(f);d2s = d2s(f); else s(f)=qm-2*qm*(ft-f)^2/ft^2;st1(f)=s(f); s = s(f); %推程减速方程式 ds(f)=4*qm*(ft-f)/ft^2;vt(f)=ds(f);ds = ds(f); d2s(f)=-4 *qm/ft^2;at(f)=d2s(f);d2s = d2s(f); end at(f)= atan((-l*(1-ds))/(a*sin((s+q0)*hd))-(-1)*cos((s+q0)*hd)/sin((s+q0)*hd));atd(f) = at(f) * du; %推程压力角的角度和弧度表达式 p1= -a*sin(f*hd)+l*sin((s+q0-f)*hd)*(ds-1); p2= a*cos(f*hd)+l*cos((s+q0-f)*hd)*(ds-1); p3=-a*cos(f*hd)+l*(ds-1)^2*cos((s+q0-f)*hd)+l*d2s*sin((s+q0-f)*hd); p4=-a*sin(f*hd)-l*(ds-1)^2*sin((s+q0-f)*hd)+l*ds*cos((s+q0-f)*hd); pt(f)= (p1^2+p2^2)^1.5/(p1*p4-p2*p3) ;p = pt(f); end atm = 0; for f = 1 : ft if atd(f) atm atm = atd(f); end end fprintf (1,