2023年機械原理大作業(yè)2凸輪機構(gòu)大作業(yè)

大作業(yè)（二）凸輪機構(gòu)設(shè)計（題號：_________）班級：________________________學(xué)號：________________________姓名：________________________同組其他人員：________________________完成日期：________________________凸輪機構(gòu)大作業(yè)題目目錄題目及原始數(shù)據(jù)；推桿旳運動規(guī)律及凸輪廓線方程；計算程序框圖計算程序；計算成果及分析；凸輪機構(gòu)圖（包括推桿及凸輪理論和實際廓線，并標(biāo)出有關(guān)尺寸及計算成果體會及提議參照書運用計算機輔助設(shè)計完畢下列偏置式直動滾子推桿盤形凸輪機構(gòu)（推桿旳移動副導(dǎo)路位于凸輪盤回轉(zhuǎn)中心旳右側(cè)）或擺動滾子推桿盤形凸輪機構(gòu)旳設(shè)計，設(shè)計已知數(shù)據(jù)如下表所示，機構(gòu)中凸輪沿著逆時針方向做勻速轉(zhuǎn)動。表1兩種凸輪機構(gòu)旳從動件運動規(guī)律直動推桿組題號推程運動規(guī)律回程運動規(guī)律5-B等加速等減速運動五次多項式運動表2兩種凸輪機構(gòu)旳推桿在近休、推程、遠休及回程階段旳凸輪轉(zhuǎn)角題號近休凸輪轉(zhuǎn)角推程凸輪轉(zhuǎn)角遠休凸輪轉(zhuǎn)角回程凸輪轉(zhuǎn)角B0゜～45゜45゜～225゜225゜～270゜270゜～360゜表3偏置直動滾子推桿盤形凸輪機構(gòu)旳已知參數(shù)題號初選基圓半徑r0/mm偏距e/mm滾子半徑rt/mm推桿行程h/mm許用壓力角許用最小曲率半徑[ρamin][α1][α2]B20+10153530゜75゜0.3rt規(guī)定：每兩人一組，每組中至少打印出一份源程序。每人都要打?。涸紨?shù)據(jù)；凸輪理論輪廓曲線和實際輪廓曲線旳坐標(biāo)值；推程和回程旳最大壓力角，以及出現(xiàn)最大壓力角時凸輪對應(yīng)旳轉(zhuǎn)角，凸輪實際輪廓曲線旳最小曲率半徑，以及對應(yīng)旳凸輪轉(zhuǎn)角；凸輪旳基圓半徑。整個設(shè)計過程所選用旳計算點數(shù)N=72~120。運用計算機繪出凸輪旳理論輪廓曲線和實際輪廓曲線。二、推桿運動規(guī)律及凸輪輪廓方程推程：等加速遠休:等減速回程:五次多項式整頓得；理論輪廓廓線方程工作廓線方程實際廓線方程三、計算程序框圖|α|>[|α|>[α1]？開始讀入：r0，Δr0，rt，h或（φ），e或（lAB、lOA）δ1，δ2，δ3，δ4，[α1]，[α2]，[ρamin]，N計算：s0I=1計算：s，x，y，ds/dδ，dx/dδ，dy/dδ，x′，y′計算：αr0=r0+Δr0r0=r0=Δr0是回程？|α|>[α2]？選出α1max及對應(yīng)旳凸輪轉(zhuǎn)角δ1選出α2max及對應(yīng)旳凸輪轉(zhuǎn)角δ2計算：ρρ<0?|ρ|-rt≥[ρamin]？計算ρa選出|ρamin|及對應(yīng)旳凸輪轉(zhuǎn)角δaminI=I+1I≤N？打?。簒，y，x′，y′，ρamin，δamin，α1max，δ1max，α2max，δ2max，r0，δ，s結(jié)束四、計算程序#include<stdio.h>#include<math.h>#include<stdlib.h>#definePI3.doublefact[72][2];doubletheory[72][2];intang1=180,ang2=225,ang3=315;doubleh=35,rb=20,b=1;doubleA1=30*PI/180,A2=75*PI/180;doubleP=18.2,e=10;doubleSo,r=15;doubleS(intI){ doubles; doubleA; doubleB; if(I<=ang1/2) { A=I*PI/180; B=ang1*PI/180; s=2*h*pow(A/B,2); } elseif((I>ang1/2)&&(I<=ang1)) { A=I*PI/180; B=ang1*PI/180; s=h-2*h*pow((B-A)/B,2); } elseif(I<=ang2)s=h; elseif(I<=ang3) { A=(I-ang2)*PI/180; B=(ang3-ang2)*PI/180; s=h*(1-10*pow(A/B,3)+15*pow(A/B,4)-6*pow(A/B,5)); } elses=0; return(s);}doubleds(intQ){ doubleA,B,C; if(Q<=ang1/2) { A=Q*PI/180; B=ang1*PI/180; C=4*h*A/(B*B); } elseif((Q>ang1/2)&&(Q<=ang1)) { A=Q*PI/180; B=ang1*PI/180; C=4*h*(B-A)/(B*B); } elseif(Q<=ang2)C=0; elseif(Q<=ang3) { A=(Q-ang2)*PI/180; B=(ang3-ang2)*PI/180; C=h*(-30*A*A/pow(B,3)+60*pow(A,3)/pow(B,4)-30*pow(A,4)/pow(B,5)); } elseC=0; returnC;}doubledss(intB3){ doubleA,B,C; if(B3<=ang1/2) { A=B3*PI/180; C=ang1*PI/180; B=4*h/(C*C); } elseif(B3>ang1/2&&B3<=ang1) { A=B3*PI/180; C=ang1*PI/180; B=-4*h/(C*C); } elseif(B3<=ang2)B=0; elseif(B3<=ang3) { A=(B3-ang2)*PI/180; C=(ang3-ang2)*PI/180; B=h*(-60*A/pow(C,3)+240*A*A/pow(C,4)-120*A*A*A/pow(C,5)); }elseB=0; return(B);}voidxy(intang){ doubleA,B,C,E,F,dx,dy; A=ang*PI/180; B=S(ang); C=ds(ang);dx=(So+B)*cos(A)+sin(A)*C-e*sin(A);dy=-sin(A)*(So+B)+C*cos(A)-e*cos(A); E=r*dy/sqrt(dx*dx+dy*dy); F=r*dx/sqrt(dx*dx+dy*dy); theory[ang/5][0]=(So+B)*sin(A)+e*cos(A); theory[ang/5][1]=(So+B)*cos(A)-e*sin(A); fact[ang/5][0]=theory[ang/5][0]-E; fact[ang/5][1]=theory[ang/5][1]+F;}doublea(intB1)/*****求解壓力角****/{ doubleA,B; A=sqrt((ds(B1)-e)*(ds(B1)-e)); B=S(B1); returnatan(A/(B+So));}doublep(intB2){ doubledx,dy,dxx,dyy; doubleA,B,C,D,E; A=B2*PI/180; B=ds(B2); C=S(B2);D=dss(B2); dx=(So+C)*cos(A)+sin(A)*B-e*sin(A);dy=-sin(A)*(So+C)+B*cos(A)-e*cos(A); dxx=-(C+So)*sin(A)+cos(A)*B+D*sin(A)-e*cos(A); dyy=-cos(A)*(So+C)-B*sin(A)+D*cos(A)-sin(A)*B+e*sin(A); E=sqrt(pow(dx*dx+dy*dy,3))/sqrt(pow((dx*dyy-dxx*dy),2)); return(E);}voidmain(){FILE*fp; inti; intk,h,l; doubleangle1max=0,angle2max=0,pmin=1000;if((fp=fopen("sanying","w"))==NULL){ printf("Cann'topenthisfile.\n"); exit(0); } fprintf(fp,"\nTheKinematicParametersofPoint4\n"); fprintf(fp,"xyx'y'");//計算數(shù)據(jù)并寫入文獻for(;i!=360;){ rb=rb+b; So=sqrt(rb*rb-e*e); for(i=0;i<=ang1;i=i+5) { if(a(i)>A1||p(i)<P) break; } if(ang1+5-i)continue; for(i=ang1+5;i<=ang2;i=i+5) { if(p(i)<P)break; } if(ang2+5-i)continue; for(i=ang2+5;i<=ang3;i=i+5) { if(a(i)>A2||p(i)<P) break; } if(ang3+5-i)continue; for(i=ang3+5;i<360;i=i+5) { if(p(i)<P) break; }} for(i=0;i<360;i=i+5) { xy(i); } for(i=0;i<=ang1;i=i+5) { if(angle1max<a(i)) { angle1max=a(i); k=i; } if(pmin>p(i)) { pmin=p(i); h=i; } } for(i=ang2;i<=ang3;i=i+5) { if(angle2max<a(i)) { angle2max=a(i); l=i; } if(pmin>p(i)) { pmin=p(i); h=i; } }for(i=0;i<72;i++){ fprintf(fp,"\n"); {fprintf(fp,"%12.3f\t%12.3f\t%12.3f\t%12.3f\t",theory[i][0],theory[i][1],fact[i][0],fact[i][1]); }} fclose(fp); printf("理論坐標(biāo)(x,y)");printf("實際坐標(biāo)(x,y)");printf("\n"); for(i=0;i<72;i++) { printf("%f",theory[i][0]); printf(""); printf("%f",theory[i][1]);printf("");printf("%f",fact[i][0]); printf(""); printf("%f",fact[i][1]); printf("\n"); } printf("基圓半徑是:%f\n",rb); printf("推程最大壓力角是:%f\n",angle1max*180/PI); printf("此時角度是是:%d\n",k); printf("回程最大壓力角是:%f\n",angle2max*180/PI); printf("此時角度是是:%d\n",l); printf("最小曲率半徑是:%f\n",pmin); printf("此時角度是:%d\n",h);}五、計算成果及分析理論坐標(biāo)x理論坐標(biāo)y實際坐標(biāo)x實際坐標(biāo)y10.000000115.56816211.293103130.51232120.039083114.31063922.471855129.11204629.953795112.28870733.510983126.86081839.696314109.51162944.356597123.76932049.219138105.99042754.955228119.85034458.474959101.73793865.253713115.11883867.41655796.76890275.199098109.59202375.99672591.10010384.738564103.28945384.16820984.75052793.81938796.23333091.88369177.741555102.38893288.44855899.09581070.097185110.39469879.963042105.75722361.844259117.78439870.807912111.82071753.012716124.50610461.017774117.23936343.635834130.50843250.630961121.96673133.750482135.74079239.689778125.95714423.397353140.15368328.240731129.16599812.621184143.69904916.334732131.5501161.470946146.3306834.027279133.068162-10.000000148.004678-8.621416133.573478-21.724363148.573242-21.808621132.927116-33.592903147.847706-35.134332131.110804-45.483796145.811871-48.463499128.119609-57.273502142.463995-61.659683123.962059-68.838237137.816884-74.586611118.660100-80.055446131.897811-87.109735112.248890-90.805265124.748264-99.097770104.776442-100.971946116.423537-110.42420596.303108-110.445244106.992155-120.96875486.900923-119.12172596.535156-130.61874876.652799-126.90601085.145215-139.27043065.651605-133.71191572.925643-146.83016053.999112-139.46348559.989261-153.21550741.804845-144.09591246.457159-158.35620529.184835-147.55631932.457365-162.19498516.260293-149.80440818.123425-164.6882493.156226-150.8129553.592926-165.806597-10.000000-150.568162-10.994037-165.535189-23.084827-149.123647-25.379543-163.947084-35.993964-146.544211-39.571896-161.111243-48.629166-142.849486-53.463083-157.049248-60.894270-138.067589-66.947384-151.792023-72.695933-132.234915-79.922175-145.379550-83.944335-125.395853-92.288709-137.860660-94.553870-117.602453-103.952871-129.292568-104.443793-108.914028-114.825889-119.740481-113.538836-99.396700-124.825012-109.277097-121.727552-89.087472-133.938976-97.798450-128.743333-77.939306-141.897783-85.147665-134.320234-66.002685-148.296955-71.447391-138.263460-53.439531-152.846268-56.952612-140.469540-40.484954-155.392740-42.000903-140.930485-27.409448-155.923766-26.960506-139.731023-14.484083-154.553686-12.184372-137.035532-1.950857-151.4989942.025104-133.06816210.000000-147.04718715.439380-128.08806221.244452-141.52279127.915887-122.36209131.730004-135.25256139.400453-116.13746341.471701-128.53101449.921547-109.61679950.539030-121.58764459.577772-102.93796259.034539-114.55843768.519504-96.16081367.065477-107.46611876.923981-89.26278174.71021