2023年机械原理大作业2凸轮机构大作业

大作业（二）凸轮机构设计（题号：_________）班级：________________________学号：________________________姓名：________________________同组其他人员：________________________完成日期：________________________凸轮机构大作业题目目录题目及原始数据；推杆旳运动规律及凸轮廓线方程；计算程序框图计算程序；计算成果及分析；凸轮机构图（包括推杆及凸轮理论和实际廓线，并标出有关尺寸及计算成果体会及提议参照书运用计算机辅助设计完毕下列偏置式直动滚子推杆盘形凸轮机构（推杆旳移动副导路位于凸轮盘回转中心旳右侧）或摆动滚子推杆盘形凸轮机构旳设计，设计已知数据如下表所示，机构中凸轮沿着逆时针方向做匀速转动。表1两种凸轮机构旳从动件运动规律直动推杆组题号推程运动规律回程运动规律5-B等加速等减速运动五次多项式运动表2两种凸轮机构旳推杆在近休、推程、远休及回程阶段旳凸轮转角题号近休凸轮转角推程凸轮转角远休凸轮转角回程凸轮转角B0゜～45゜45゜～225゜225゜～270゜270゜～360゜表3偏置直动滚子推杆盘形凸轮机构旳已知参数题号初选基圆半径r0/mm偏距e/mm滚子半径rt/mm推杆行程h/mm许用压力角许用最小曲率半径[ρamin][α1][α2]B20+10153530゜75゜0.3rt规定：每两人一组，每组中至少打印出一份源程序。每人都要打印：原始数据；凸轮理论轮廓曲线和实际轮廓曲线旳坐标值；推程和回程旳最大压力角，以及出现最大压力角时凸轮对应旳转角，凸轮实际轮廓曲线旳最小曲率半径，以及对应旳凸轮转角；凸轮旳基圆半径。整个设计过程所选用旳计算点数N=72~120。运用计算机绘出凸轮旳理论轮廓曲线和实际轮廓曲线。二、推杆运动规律及凸轮轮廓方程推程：等加速远休:等减速回程:五次多项式整顿得；理论轮廓廓线方程工作廓线方程实际廓线方程三、计算程序框图|α|>[|α|>[α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及对应旳凸轮转角δ1选出α2max及对应旳凸轮转角δ2计算：ρρ<0?|ρ|-rt≥[ρamin]？计算ρa选出|ρamin|及对应旳凸轮转角δaminI=I+1I≤N？打印：x，y，x′，y′，ρamin，δamin，α1max，δ1max，α2max，δ2max，r0，δ，s结束四、计算程序#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'");//计算数据并写入文献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("理论坐标(x,y)");printf("实际坐标(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);}五、计算成果及分析理论坐标x理论坐标y实际坐标x实际坐标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