機械原理大作業(yè) 凸輪設計

...wd......wd......wd...目錄TOC\o"1-1"\h\u20436一、題目及原始數(shù)據(jù)32206二、推桿運動規(guī)律及凸輪廓線方程414500三、計算程序55267四、計算結果及分析122867五、凸輪機構圖1223738六、體會及建議152527七、參考書15題目及原始數(shù)據(jù)試用計算機輔助設計完成偏置直動滾子推桿盤形凸輪機構的設計表1凸輪構造的推桿運動規(guī)律直動推桿組題號推程運動規(guī)律回程運動規(guī)律10-A等加速等減速運動五次多項式運動表2凸輪構造的推桿在近休、推程、遠休及回程段的凸輪轉(zhuǎn)角題號近休凸輪轉(zhuǎn)角推程凸輪轉(zhuǎn)角遠休凸輪轉(zhuǎn)角回程凸輪轉(zhuǎn)角10-A0°~60°60°~180°180°~270°270°~360°表3偏置直動滾子推桿盤形凸輪機構的參數(shù)題號初選的基圓半徑r0/mm偏距e/mm滾子半徑Rr/mm推桿行程h/mm許用壓力角許用最小曲率半徑[ρamin][α1][α2]10-A15+5102830700.3rr〔1〕打印出原始數(shù)據(jù)；〔2〕打印出理論輪廓和實際輪廓的坐標值；〔3〕打印出推程和回程的最大壓力角，以及出現(xiàn)最大壓力角時凸輪的相應轉(zhuǎn)角；〔4〕打印出凸輪實際輪廓曲線的最小曲率半徑，以及相應的凸輪轉(zhuǎn)角；〔5〕打印最后所確定的凸輪的基圓半徑。計算點數(shù)：100二、推桿運動規(guī)律及凸輪廓線方程1、推桿運動規(guī)律近休階段：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組數(shù)據(jù)",j+1);/*輸出數(shù)據(jù)*/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,相應凸輪轉(zhuǎn)角=%f\n",A1,A2-4);printf("回程最大壓力角(弧度)=%f,相應凸輪轉(zhuǎn)角=%f\n",B1,B2-4);printf("最小曲率半徑=%f,相應凸輪轉(zhuǎn)角=%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