凸轮机构大作业西工大机械原理.doc

大作业（） （题号：） （一） 题目及原始数据···············（二） 推杆运动规律及凸轮廓线方程········· （）········· （） 计算程序·················（） 程序计算结果及分析·············（） 凸轮机构图·················（） 心得体会··················（） 参考书··················· 题目及原始数据 试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计 （1）推程运动规律为五次多项式运动规律，回程运动规律为运动规律； （2）打印出原始数据； （3）打印出理论轮廓和实际轮廓的坐标值； （4）打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角； （5）打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角； （）打印最后所确定的凸轮的基圆半径。 R0/mm 偏距 E/mm 滚子半径Rr/mm 推杆行程 h/mm 许用压力角 许用最小曲率半径[ρamin] [α1] [α2] 4-A 15 5 10 28 30° 70? 0.3Rr 计算点数：N=90 q1=60; 近休止角δ1 =180; 推程运动角δ2 =90; 远休止角δ3 =90; 回程运动角δ4 二 推杆运动规律及凸轮廓线方程 四 计算程序 1. #includestdio.h #includemath.h void main(){ double r0,or,rr,h,e,q1,q2,q3,q4,a,a11,a22,Q,pi,pa,paa,QQ,A1,A2,B1,B2,C1,C2; /*定义变量*/ double xz[90],yz[90],sz[90],x1z[90],y1z[90],Q1,Q2; double s0,s,x,y,y1,x1,dx,dxx,dy,dyy,ds,dss,sino,coso,p; int N,i,j; r0=19;e=5;h=28;rr=10;q1=60;q2=120;q3=90;q4=90;a11=30;a22=70;or=1;pi=3.141592653;pa=3; /*给已知量赋值*/ N=90;A1=0;B1=0;C1=1000; for(; ;){ Q=0; C1=1000; QQ=180/pi; r0=r0+or; s0=sqrt(r0*r0-e*e); for(i=1,j=0;i=N;i++,j++){ if(Q60){ /*近休阶段*/ s=0; ds=0; dss=0; a=atan(e/sqrt(r0*r0-e*e)); /*求压力角*/ if(aa11/QQ){ break; } else{ if(aA1) A1=a; A2=Q; } } else if(Q=60Q180){ /*五次多项式运动*/ Q1=Q-60; s=10*h*Q1*Q1*Q1/(q2*q2*q2)-15*h*Q1*Q1*Q1*Q1/(q2*q2*q2*q2)+6*h*Q1*Q1*Q1*Q1*Q1/(q2*q2*q2*q2*q2); ds=30*h*Q1*Q1*QQ/(q2*q2*q2)-60*h*Q1*Q1*Q1*QQ/(q2*q2*q2*q2)+30*h*Q1*Q1*Q1*Q1*QQ/(q2*q2*q2*q2*q2); dss=60*h*Q1*QQ*QQ/(q2*q2*q2)-180*h*Q1*Q1*QQ*QQ/((q2*q2*q2*q2))+120*h*Q1*Q1*Q1*QQ*QQ/((q2*q2*q2*q2*q2)); a=atan(fabs(ds-e)/(sqrt(r0*r0-e*e)+s)); if(aa11/QQ){ break; } else{ /*远休阶段*/ if(aA1) A1=a; A2=Q; } } else if(Q=180Q270){ s=28; ds=0;dss=0; a=atan(fabs(ds-e)/(sqrt(r0*r0-e*e)+s)); if(aa22/QQ){ break; } else{ if(aB1) B1=a; B2=Q;