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

《滚子摆动从动件凸轮设计matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、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.4fmm',rb)fprintf(1,'滚子半径rt=%3.4fmm',rt)fprintf(1,'

2、起始角度q0=%3.4fmm',q0)fprintf(1,'最大摆动角度qm=%3.4fmm',qm)fprintf(1,'推程运动角ft=%3.4f度',ft)fprintf(1,'远休止角fs=%3.4f度',fs)fprintf(1,'回程运动角fh=%3.4f度',fh)fprintf(1,'推程许用压力角alp=%3.4f度',alp)hd=pi/180;du=180/pi;%角度弧度互换d1=ft+fs;d2=ft+fs+fh;disp''disp'计算过程和输出结

3、果:'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);forf=1:ftiff<=ft/2s(f)=2*(qm/ft^2)*f^2;st1(f)=s(f);s=s(f);%推程加速方程式ds(f)=(qm/ft^2)*f;

4、vt(f)=ds(f);ds=ds(f);d2s(f)=4*qm/ft;at(f)=d2s(f);d2s=d2s(f);elses(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);endat(f)=atan((-l*(1-ds))/(a*sin((s+q0)*hd))-(-1)

5、*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)*

6、hd)+l*ds*cos((s+q0-f)*hd);pt(f)=(p1^2+p2^2)^1.5/(p1*p4-p2*p3);p=pt(f);endatm=0;forf=1:ftifatd(f)>atmatm=atd(f);endendfprintf(1,'最大压力角atm=%3.4f度',atm)forf=1:ftifabs(atd(f)-atm)<0.1ftm=f;breakendendfprintf(1,'对应的位置角ftm=%3.4f度',ftm)ifatm>alpfprintf(1,'

7、*凸轮推程压力角超过许用值,需要增大基圆!')endptn=rb;ftn=0;forf=1:ftifpt(f)

8、disp'1-2回程(等加速等减速运动)'s=zeros(fh);ds=zeros(fh);d2s=zeros(fh);ah=zeros(fh);ahd=zeros(fh);ph=zeros(fh);forf=d1:d2k=f-d1;ifk<=fh/2s(f)=qm-2*qm*(k)^2/fh^2;st1(f)=s(f);s=s(f);ds(f)=-4*qm*k/fh^2;ds=ds(f);d2s(f)=-4*qm/fh^2;d2s=d2s(f);els