大作业——凸轮机构.docVIP

一、原始数据 题号 11---B 推程运动规律 回程运动规律 正弦加速度运动 余弦加速度运动 近休凸轮转角 推程凸轮转角 远休凸轮转角 回程凸轮转角 0゜～60゜ 60゜～210゜ 210゜～270゜ 270゜～360゜ 初选基圆半径r0/mm 偏距e/mm 滚子半径rt/mm 推杆行程h/mm 许用压力角 许用最小曲率半径[ρamin] [α1] [α2] 15 +10 12 30 30゜ 75゜ 0.3rt 动规律及凸轮廓线方程 推程 正弦加速度运动规律 回程 余弦加速度运动规律 理论轮廓廓线方程 x=(s+s)sinδ+ecosδ y= (s+s)cosδ-esinδ 工作廓线方程 x=x-cosδ y=y-sinδ 实际廓线方程： 原始程序 #includestdio.h #includestdlib.h #includemath.h #define Pi 3.141592653 double r0=15,e=10,rt=12,h=30,t1=60,t2=210,t3=270,t4=360,a1=30,a2=75,P=3.6,N=120; void main() { int i,j; double zz,t,m,n,a,b,pa,s0,ds,dx,dy,d2x,d2y,d2s,T1,T2; double a1max=0,a2max=0,t1max,t2max,pamin=1000,tamin,x[4][121],s[121]; zz=Pi; T1=((t2-t1)*zz)/180; T2=((t4-t3)*zz)/180; flag: s0=sqrt(r0*r0-e*e); for(i=0;i=N;i++) { t=(2*zz*i)/N; if(i=0i=((t1*N)/360)) /*近休阶段*/ { s[i]=0; ds=0; d2s=0; } else if(i((t1*N)/360)i=((t2*N)/360)) /*推程阶段*/ { t=t-(t1*zz)/180; s[i]= h*(t/T1-sin((2*zz*t)/T1)/(2*zz)); ds= h/T1-(h*cos((2*zz*t)/T1))/T1; d2s=(2*zz*h*sin((2*zz*t)/T1))/(pow(T1,2)); b=(ds-e)/(sqrt(r0*r0-e*e)+s[i]); /*计算压力角*/ a=fabs(atan(b)); if(fabs(a)=((a1*zz)/180)) /*判断压力角*/ { r0=r0+e; goto flag; } else if(a1maxa) { a1max=a; t1max=t+(t1*zz)/180; } } else if(i((t2*N)/360)i=((t3*N)/360)) /*远休阶段*/ { s[i]=h; ds=0; d2s=0; } else if(i((t3*N)/360)i=((t4*N)/360)) /*回程阶段*/ { t=t-(t3*zz)/180; s[i]= (h*(1+cos((zz*t)/(T2))))/2 ; ds=- (h*zz*sin((zz*t)/(T2)))/(2*(T2)); d2s=-(h*zz*zz*cos((zz*t)/(T2)))/(2*pow((T2),2)); b=(ds-e)/(sqrt(r0*r0-e*e)+s[i]); /*计算压力角*/ a=fabs(atan(b)); if(fabs(a)=((a2*zz)/180)) /*判断压力角*/ { r0=r0+e; goto flag; } else if(a2maxa) { a2max=a; t2max=t+(t3*zz)/180; } } else { prin