《机械原理课程设计——牛头刨床 2 待》由会员分享，可在线阅读，更多相关《机械原理课程设计——牛头刨床 2 待（5页珍藏版）》请在金锄头文库上搜索。
1、【转】机械原理课程设计牛头刨床 2 待【转】机械原理课程设计牛头刨床(2)待续2011年05月23日2.5源程序及计算结果 1）程序 #includestdio.h #includestdlib.h #includemath.h const double PI=3.14159; double L10; double X10,Y10; double V10,U10; double A10,B10; double F10,W10,E10; double S10,C10; double Sgn(double Xin) double Resf; if(Xin=0) Resf=1.0; if(Xin1e
2、-10) Resf=atan(Yin/Xin); Resf=Resf-(Sgn(Xin)-1)*PI/2; else Resf=PI/2; Resf=Resf-(Sgn(Yin)-1)*Resf; return(Resf); void mcrank(int i,int j,int a,int b,double F9) Fj= Fj+F9; Si=Li*sin(Fj); Ci=Li*cos(Fj); Xb=Xa+Ci; Yb=Ya+Si; Vb=Va-Wj*Si; Ub=Ua+Wj*Ci; Ab=Aa-Wj*Wj*Ci-Ej*Si; Bb=Ba-Wj*Wj*Ci+Ej*Si; int mrpr(
3、int i,int j,int k,int b,int c,int d,int e, int m,double Res3) double A0,B0,C0,X1,Y1,F1,Ar,Ak; double G1,G4,G5,G6,s1,v1,a1; A0=Xb-Xd; B0=Yb-Yd; C0=Li+Lk; G1=A0*A0+B0*B0-C0*C0; if(G10) return(0); s1=sqrt(G1); X1=C0-B0; Y1=A0+m*s1; F1=Angle(X1,Y1); if(F1PI|F10) Fj=2*(F1+Sgn(X1)*PI); if(fabs(F1)0.001) F
4、j=2*PI; Si=Li*sin(Fj); Ci=Li*cos(Fj); Sk=Lk*sin(Fj); Ck=Lk*cos(Fj); Sj=Lj*sin(Fj); Cj=Lj*cos(Fj); Xc=Xb-Si; Yc=Yb+Ci; Xe=Xc+Cj-s1*cos(Fj); Ye=Yc+Sj-s1*sin(Fj); G6=(Xb-Xd)*cos(Fj)+(Yb-Yd)*sin(Fj); Wj=(Ub-Ud)*cos(Fj)-(Vb-Vd)*sin(Fj) /G6; v1=(Vb-Vd)*(Xb-Xd)+(Ub-Ud)*(Yb-Yd )/G6; Vc=Vb-Wj*Ci; Uc=Ub-Wj*Si
5、; Ve=Vd-Wj*(Sj-Ck); Ue=Ud+Wj*(Cj+Sk); G4=Ab-Ad+Wj*Wj*(Xb-Xd)+2*Wj*v1*sin(F j); G5=Bb-Bd+Wj*Wj*(Xb-Xd)-2*Wj*v1*cos(F j); Ej=(G5*cos(Fj)-G4*sin(Fj)/G6; a1=(G4*(Xb-Xd)+G5*(Yb-Yd)/G6; Ar=a1; Ak=2*Wj*v1; Ae=Ad-Ej*(Sj-Ck)-Wj*Wj*(Cj+Sk); Be=Bd+Ej*(Cj+Sk)-Wj*Wj*(Sj-Ck); Res0=s1; Res1=v1; Res2=a1; return(1);
6、 int mrrp(int i,int j,int b,int c,int r,int m) double B0,C0,Z1,S1,X1,Y1,F1; double Q1,Q2,Q3,Q4,Q5,A1,V1; B0=2*(Xr-Xb)*cos(Fj)+2*(Yr-Yb)*sin(Fj) ; Sj=Lj*sin(Fj); Cj=Lj*cos(Fj); C0=pow(Xr-Xb),2)+pow(Yr-Yb),2)+ pow(Lj,2)- pow(Li,2)-2*(Xr-Xb)*Sj+2*(Yr-Yb)*Cj; if(B0*B0-4*C00) return(0); Z1=sqrt(B0*B0-4*C
7、0); S1=(-B0+m*Z1)/2; Xc=Xr+S1*cos(Fj)-Sj; Yc=Yr+S1*sin(Fj)+Cj; X1=Xc-Xb; Y1=Yc-Yb; F1=Angle(X1,Y1); Fi=F1; Si=Li*sin(Fi); Ci=Li*cos(Fi); Q1=Vr-Vb-Wj*(S1*sin(Fj)+Cj); Q2=Ur-Ub+Wj*(S1*cos(Fj)-Sj); Q3=Si*sin(Fj)+Ci*cos(Fj); Wi=(-Q1*sin(Fj)+Q2*cos(Fj)/Q3; V1=-(Q1*Ci+Q2*Si)/Q3; Vc=Vb-Wi*Si; Uc=Ub+Wi*Ci;
8、Q4=Ar-Ab+Ci*pow(Wi,2)-Ej*(S1*sin(Fj)+C j)-pow(Wj,2)*(S1*cos(Fj)-Sj)-2*Wj*V1*sin( Fj); Q5=Br-Bb+Si*pow(Wi,2)+Ej*(S1*cos(Fj)-S j)-pow(Wj,2)*(S1*sin(Fj)+Cj)+2*Wj*V1*cos( Fj); A1=(-Q4*Ci-Q5*Si)/Q3; Ei=(-Q4*sin(Fj)+Q5*cos(Fj)/Q3; Ac=Ab-Ei*Si-Ci*(Wi,2); Bc=Bb+Ei*Ci-Si*(Wi,2); return(1); void main() int ii
9、,index,iFlagea,iFlageb; double p1,F9,Res3,N1,K,M,N,P,T,R; p1=PI/180; L1=90; L2=0; L3=580; L4=0;L5=174; L6=0; X1=0;Y1=350;N1=64;X4=0;Y4=0;X7=0; printf(L1=90; L2=0; L3=580; L4=0;L5=174; L6=0;n); printf(F1DEG X6mm Y6mm V6m/s A6m/s2n); T=sqrt(Y1*Y1-L1*L1); P=T*L3/Y1; R=(L3-P)/2; Y7=L3-R; W1=-N1*PI/30; M
10、=L1/Y1; K=asin(M); F9=0;F1=-PI+K; mcrank(1,1,1,2,F9); iFlagea=mrpr(2,3,4,2,3,4,5,1,Res); if(iFlagea=0) printf(Because of wrong data,the Caculation failedn); F6=0; iFlageb=mrrp(5,6,5,6,7,1); N=X6; X1=-N; Y1=-Y7+Y1; X4=-N; Y4=-Y7; X7=-N; Y7=0; for(ii=0;ii=12;ii+) F1=-PI+K+ii*(-30)*p1; F9=0; mcrank(1,1
11、,1,2,F9); iFlagea=mrpr(2,3,4,2,3,4,5,1,Res); if(iFlagea=0) printf(Because of wrong data,the Caculation failedn); F6=0; iFlageb=mrrp(5,6,5,6,7,1); if(iFlageb=1) printf(%8.2f,%8.2f,%8.2f,%8.2f,%8.2fn,-(F1-K+ PI)/p1,X6,Y6,V6/1000,A6/1000); else printf(Because of wrong data,the Caculation failed!n); get
12、ch(); 2）计算结果 L1=90; L2=0; L3=580; L4=0;L5=174; L6=0; F1DEG X6mm Y6mm V6m/s A6m/s2 0.00, 0.00, 0.00, -0.00, 7.69 30.00, 17.89, 0.00, 0.42, 5.41 60.00, 61.37, 0.00, 0.67, 2.96 90.00,118.51, 0.00, 0.78, 0.71 120.00,180.27, 0.00, 0.78, -1.17 150.00,237.76, 0.00, 0.67, -2.50 180.00,280.99, 0.00, 0.41, -3.50 210.00,298.28, 0.00, -0.00, -4.88 240.00,275.04, 0.00, -0.62, -6.73 270.00,200.09, 0.00, -1.25, -4.27 300.00, 97.78, 0.00, -1.23, 4.18 330.00, 23.38, 0.00, -0.63, 8.41 360.00, 0.0
