机械原理大作业 凸轮设计

...wd......wd......wd...目录TOC\o"1-1"\h\u20436一、题目及原始数据32206二、推杆运动规律及凸轮廓线方程414500三、计算程序55267四、计算结果及分析122867五、凸轮机构图1223738六、体会及建议152527七、参考书15题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计表1凸轮构造的推杆运动规律直动推杆组题号推程运动规律回程运动规律10-A等加速等减速运动五次多项式运动表2凸轮构造的推杆在近休、推程、远休及回程段的凸轮转角题号近休凸轮转角推程凸轮转角远休凸轮转角回程凸轮转角10-A0°~60°60°~180°180°~270°270°~360°表3偏置直动滚子推杆盘形凸轮机构的参数题号初选的基圆半径r0/mm偏距e/mm滚子半径Rr/mm推杆行程h/mm许用压力角许用最小曲率半径[ρamin][α1][α2]10-A15+5102830700.3rr〔1〕打印出原始数据；〔2〕打印出理论轮廓和实际轮廓的坐标值；〔3〕打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；〔4〕打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；〔5〕打印最后所确定的凸轮的基圆半径。计算点数：100二、推杆运动规律及凸轮廓线方程1、推杆运动规律近休阶段：0°≤<60°推程阶段：60°≤<180°等加速段运动方程s=2h远休阶段：210°≤<280°回程阶段：280°≤<360°凸轮廓线方程理论廓线方程工作廓线方程x′=y′=计算程序#include<stdio.h>#include<math.h>voidmain(){//freopen("xxx.txt","w",stdout); doubler0,dr,rr,h,e,q1,q2,q3,q4,a,a11,a22,Q,pi,pa,paa,QQ,A1,A2,B1,B2,C1,C2;/*定义变量*/ doublexz[90],yz[90],sz[90],x1z[90],y1z[90],Q1,Q2; doubles0,s,x,y,y1,x1,dx,dxx,dy,dyy,ds,dss,sino,coso,p; intN,i,j; r0=22;e=14;h=35;rr=18;q1=30;q2=180;q3=70;q4=80;a11=35;a22=65;dr=1;pi=3.141592653;pa=6.3;/*给量赋值*/ N=90;A1=0;B1=0;C1=1000; for(;;){ Q=0; C1=1000; QQ=180/pi; r0=r0+dr; s0=sqrt(r0*r0-e*e); for(i=1,j=0;i<=N;i++,j++){ if(Q<30){/*近休阶段*/ s=0; ds=0; dss=0; a=atan(e/sqrt(r0*r0-e*e));/*求压力角*/ if(a>a11/QQ){ break; } else{ if(a>A1) A1=a; A2=Q; } } elseif(Q>=30&&Q<210){/*推程阶段*/Q1=Q-30;s=h*(1-cos(pi*Q1/q2))/2;ds=pi*h*sin(pi*Q1/q2)/(2*q2);dss=pi*pi*h*cos(pi*Q1/q2)/(2*q2*q2);a=atan(fabs(ds-e)/(sqrt(r0*r0-e*e)+s));if(a>a11/QQ){break;}else{/*远休阶段*/if(a>A1)A1=a;A2=Q;} } elseif(Q>=210&&Q<280){s=35;ds=0;dss=0;a=atan(fabs(ds-e)/(sqrt(r0*r0-e*e)+s));if(a>a22/QQ){break;}else{if(a>B1)B1=a;B2=Q;}}elseif(Q>=280&&Q<360){/*回程阶段*/Q2=Q-280;s=h*(1-(Q2/q4)+sin(2*pi*Q2/q4)/(2*pi));ds=h*(cos(2*pi*Q2/q4)-1)/q4;dss=-2*pi*h*sin(2*pi*Q2/q4)/(q4*q4);a=atan(fabs(ds-e)/(sqrt(r0*r0-e*e)+s));if(a>a22/QQ){break;}else{if(a>B1)B1=a;B2=Q;}}dx=(ds-e)*sin(Q/QQ)+(s0+s)*cos(Q/QQ);dy=(ds-e)*cos(Q/QQ)-(s0+s)*sin(Q/QQ);dxx=dss*sin(Q/QQ)+(ds-e)*cos(Q/QQ)+ds*cos(Q/QQ)-(s0+s)*sin(Q/QQ);dyy=dss*cos(Q/QQ)-(ds-e)*sin(Q/QQ)-ds*sin(Q/QQ)-(s0+s)*cos(Q/QQ);sino=dx/(sqrt(dx*dx+dy*dy));coso=-dy/(sqrt(dx*dx+dy*dy));x=(s0+s)*sin(Q/QQ)+e*cos(Q/QQ);y=(s0+s)*cos(Q/QQ)-e*sin(Q/QQ);x1=x-rr*coso;y1=y-rr*sino;sz[j]=s;yz[j]=y;xz[j]=x;x1z[j]=x1;y1z[j]=y1;p=pow(dx*dx+dy*dy,1.5)/(dx*dyy-dy*dxx);/*求理论轮廓曲率半径*/if(p<0){paa=(fabs(p)-rr);if(paa<pa){break;}else{if(paa<C1)C1=paa;C2=Q;}}Q=Q+4;}if(i==91){break;} } for(j=0;j<90;j++){printf("第%d组数据",j+1);/*输出数据*/printf("s=%f",sz[j]);printf("x=%f,y=%f;",xz[j],yz[j]);printf("x1=%f,y1=%f\n",x1z[j],y1z[j]);}printf("r0=%f\n",r0);printf("推程最大压力角(弧度)=%f,相应凸轮转角=%f\n",A1,A2-4);printf("回程最大压力角(弧度)=%f,相应凸轮转角=%f\n",B1,B2-4);printf("最小曲率半径=%f,相应凸轮转角=%f\n",C1,C2-4);}计算结果及分析凸轮机构图Matlab绘图x=[14.00000015.41071516.74635018.00039919.16675220.23972621.21409422.08510922.85417723.55698024.20913624.81984125.39659025.94493626.46830026.96783627.44235227.88829428.29978728.66873228.98497029.23649029.40969229.48969729.46068729.30628829.00996928.55545727.92716627.11061326.09283624.86278223.41167821.73335519.82453917.68509115.31818912.7304549.9320146.9365043.7609950.425869-3.045381-6.626403-10.288302-14.000000-17.728645-21.440063-25.099229-28.670776-32.119500-35.410884-38.511609-41.395702-44.073115-46.535808-48.771783-50.770147-52.521164-54.016303-55.248280-56.211093-56.900051-57.311798-57.444328-57.296994-56.870516-56.166970-55.189784-53.943719-52.434845-50.642721-48.444586-45.721768-42.406762-38.494686-34.046072-29.180734-24.063677-18.885153-13.837781-9.093986-4.786983-0.9978982.2492484.9841887.2793829.23316510.94973214.000000];y=[20.71231519.68527018.56232117.34893816.05103214.67492813.22732811.71528710.1552108.6051797.0711915.5479064.0283342.5041510.966043-0.595920-2.191898-3.831884-5.525313-7.280673-9.105140-11.004237-12.981525-15.038337-17.173559-19.383462-21.661588-23.998702-26.382794-28.799157-31.230513-33.657210-36.057462-38.407658-40.682702-42.856404-44.901901-46.792104-48.500165-49.999950-51.266511-52.276556-53.008891-53.444842-53.568639-53.367760-52.833222-51.959822-50.746310-49.195504-47.314339-45.113841-42.609040-39.827853-36.843216-33.679083-30.350869-26.874788-23.267776-19.547405-15.731802-11.839555-7.889627-3.9012610.1061114.1129668.09978412.04714015.93580319.74682923.46165027.05523930.46107133.57228036.25958738.39072539.85141740.56474340.50589039.70994038.27135836.33502034.07989331.69753629.36841227.23931725.40515023.89765222.68275820.712315];x1=[3.9200004.3150004.6889785.0401125.3666905.6671235.9399466.1838316.4075656.6562456.9399577.2629337.6276238.0347388.4832688.9705119.49207310.04189210.61225111.19382211.77573712.34568912.89008113.39421413.84252214.21884314.50673314.68980014.75207014.67835414.45462514.06838413.50901212.76808811.83968410.7206109.4106097.9125036.2322794.3791102.3653190.206274-2.079773-4.471913-6.946854-9.479256-12.042104-14.607122-17.145218-19.626947-22.022998-24.304672-26.444369-28.424511-30.262966-31.953983-33.489324-34.861508-36.063850-37.090493-37.936434-38.597553-39.070629-39.353357-39.444359-39.343192-39.050348-38.567256-37.896268-37.040653-36.004579-34.769401-33.220147-31.245851-28.787411-25.848578-22.498389-18.864538-15.118342-11.453008-8.057698-5.090569-2.654686-0.7817620.5707331.5052552.1546942.6474973.0799863.920000];y1=[5.7994485.5118765.1974504.8577034.4942894.1089803.7036523.2802802.8402312.4112191.9942371.5786531.1536480.7083800.232122-0.285610-0.854923-1.485412-2.185991-2.964711-3.828560-4.783237-5.832928-6.980083-8.225211-9.566690-11.000625-12.520742-14.118339-15.782286-17.499088-19.253010-21.026255-22.799200-24.550682-26.258332-27.898945-29.448882-30.884491-32.182550-33.320701-34.277896-35.034810-35.574247-35.881505-35.944704-35.755071-35.307167-34.599059-33.632433-32.412639-30.948670-29.253078-27.347942-25.298530-23.125866-20.840535-18.453672-15.976904-13.422298-10.802300-8.129674-5.417441-2.6788150.07