滾子擺動從動件凸輪設(shè)計matlab程序

disp'********滾子擺動從動件凸輪設(shè)計********'disp'已知條件:'disp'凸輪作順時針方向轉(zhuǎn)動,從動件做擺動'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=%mm\n',rb)fprintf(1,'滾子半徑rt=%mm\n',rt)fprintf(1,'起始角度q0=%mm\n',q0)fprintf(1,'最大擺動角度qm=%mm\n',qm)fprintf(1,'推程運動角ft=%度\n',ft)fprintf(1,'遠休止角fs=%度\n',fs)fprintf(1,'回程運動角fh=%度\n',fh)fprintf(1,'推程許用壓力角alp=%度\n',alp)hd=pi/180;du=180/pi;%角度弧度互換d1=ft+fs;d2=ft+fs+fh;disp''disp'計算過程和輸出結(jié)果:'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;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)*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)^(p1*p4-p2*p3);p=pt(f);endatm=0;forf=1:ftifatd(f)>atmatm=atd(f);endendfprintf(1,'最大壓力角atm=%度\n',atm)forf=1:ftifabs(atd(f)-atm)<ftm=f;breakendendfprintf(1,'對應(yīng)的位置角ftm=%度\n',ftm)ifatm>alpfprintf(1,'*凸輪推程壓力角超過許用值,需要增大基圓!\n')endptn=rb;ftn=0;forf=1:ftifpt(f)<ptnptn=pt(f);endendfprintf(1,'輪廓最小曲率半徑ptn=%mm\n',ptn)forf=1:ftifabs(pt(f)-ptn)<ftn=f;breakendendfprintf(1,'對應(yīng)的位置角ftn=%度\n',ftn)ifptn<rt+5fprintf(1,'*凸輪推程輪廓曲率半徑小于許用值,需要增大基圓或減小滾子!\n')enddisp'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);elses(f)=2*qm*(d2-f)^2/fh^2;st1(f)=s(f);s=s(f);ds(f)=-4*qm*(d2-f)/fh^2;ds=ds(f);d2s(f)=4*qm/fh^2;d2s=d2s(f);endat(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)^(p1*p4-p2*p3);p=pt(f);endatm=0;forf=1:ftifatd(f)>atmatm=atd(f);endendfprintf(1,'最大壓力角atm=%度\n',atm)forf=1:ftifabs(atd(f)-atm)<ftm=f;breakendendfprintf(1,'對應(yīng)的位置角ftm=%度\n',ftm)ifatm>alpfprintf(1,'*凸輪推程壓力角超過許用值,需要增大基圓!\n')endptn=rb;ftn=0;forf=1:ftifpt(f)<ptnptn=pt(f);endendfprintf(1,'輪廓最小曲率半徑ptn=%mm\n',ptn)forf=1:ftifabs(pt(f)-ptn)<ftn=f;breakendendfprintf(1,'對應(yīng)的位置角ftn=%度\n',ftn)ifptn<rt+5fprintf(1,'*凸輪推程輪廓曲率半徑小于許用值,需要增大基圓或減小滾子!\n')enddisp'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);forf=1:niff<=ft/2s(f)=2*(qm/ft^2)*f^2;s=s(f);ds(f)=(qm/ft^2)*f;ds=ds(f);elseiff>ft/2&f<=fts(f)=qm-2*qm*(ft-f)^2/ft^2;s=s(f);ds(f)=4*qm*(ft-f)/ft^2;ds=ds(f);elseiff>ft&f<=d1s=qm;ds=0;elseiff>d1&f<=(d2-fh/2)k=f-d1;s(f)=qm-2*qm*(k)^2/fh^2;s=s(f);ds(f)=-4*qm*k/fh^2;ds=ds(f);elseiff>(d2-fh/2)&f<=d2s(f)=2*qm*(d2-f)^2/fh^2;s=s(f);ds(f)=-4*qm*(d2-f)/fh^2;ds=ds(f);elseiff>d2&f<=ns=0;ds=0;endxx(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(dx^2+dy^2));xxp=xp(f);yp(f)=y+rt*(dx/sqrt(dx^2+dy^2));yyp=yp(f);r(f)=sqrt(x^2+y^2);rp(f)=sqrt(xxp^2+yyp^2);enddisp'2-1推程(等加速/等減速運動)'disp'凸輪轉(zhuǎn)角理論x理論y實際x實際y'forf=10:10:ftnu=[fxx(f)yy(f)xp(f)yp(f)];disp(nu)enddisp'2-2回程(等加速/等減速運動)'disp'凸輪轉(zhuǎn)角理論x理論y實際x實際y'forf=d1:10:d2nu=[fxx(f)yy(f)xp(f)yp(f)];disp(nu)enddisp'2-3凸輪輪廓向徑'disp'凸輪轉(zhuǎn)角理論r實際r'forf=10:10:nnu=[fr(f)rp(f)];disp(nu)enddisp'繪制凸輪的理論輪廓和實際輪廓:'plot(xx,yy,'r-.')%理論輪廓(紅色,點劃線)axis([-(150)(150)-(150)(150)])%橫軸和縱軸的下限和上限