块状物品推送机 机械原理课程设计

A=I*PI/180; B=ang1*PI/180; s=h-2*h*pow((B-A)/B,2); } else if(I<=ang2) { A=(I-ang1)*PI/180; B=(ang2-ang1)*PI/180; s=h-(10*h*pow(A/B,3)-15*h*pow(A/B,4)+6*h*pow(A/B,5)); } else s=0; return(s); }

double ds(int Q) { double A,B,C; if(Q<=ang1/2) { A=Q*PI/180; B=ang1*PI/180; C=4*h*A/(B*B); } else if((Q>ang1/2)&&(Q<=ang1)) { A=Q*PI/180; B=ang1*PI/180; C=4*h*(B-A)/(B*B); } else if(Q<=ang2) { A=(Q-ang1)*PI/180; B=(ang2-ang1)*PI/180; C=-30*h*A*A/pow(B,3)+60*h*pow(A,3)/pow(B,4)-30*h*pow(A,4)/pow(B,5); } else C=0; return C; }

double dss(int B3) { double A,B,C; if(B3<=ang1/2) { A=B3*PI/180; C=ang1*PI/180;

B=4*h/(C*C); } else if(B3>ang1/2&&B3<=ang1) { A=B3*PI/180; C=ang1*PI/180; B=-4*h/(C*C); } else if(B3<=ang2) { A=(B3-ang1)*PI/180; C=(ang2-ang1)*PI/180; B=-60*h* A/pow(C,3)+180*h*A*A/pow(C,4)-120*h*pow(A,3)/pow(C,5); }

else B=0; return(B); }

void xy(int ang) { double A,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; }

double a(int B1)/*****求解压力角****/ { double A,B; A=sqrt((ds(B1)-e)*(ds(B1)-e)); B=S(B1); return atan(A/(B+So)); }

double p(int B2) { double dx,dy,dxx,dyy; double A,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); }

//计算数据并写入文件 void main() { FILE *fp; int i; int k,h,l; double angle1max=0,angle2max=0,pmin=1000; if((fp=fopen(\ { printf(\ exit(0); }

fprintf(fp,\ The Kinematic Parameters of Point 4\\n\ fprintf(fp,\ x y x' 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)A2||p(i)

y' \ if(p(i)p(i)) { pmin=p(i); h=i; } } for(i=ang1;i<=ang2;i=i+5) { if(angle2maxp(i)) { pmin=p(i); h=i; } } for(i=0;i<72;i++) { fprintf(fp,\ {

fprintf(fp,\ } } fclose(fp);