吉大材力设计传动轴6-e

/30 材料力学课程设计{return(sqrt(pow(a,2)+pow(b,2)));}floatmax2(floata,floatb){return(a>b?a:b);}floatmax3(floata,floatb,floatc){return(max2(max2(a,b),c));}floatqd(floata)/*求直径的函数*/{floate,D,ks=0.08;e=32/Pi/ks*a;D=pow(e,1.0/3);return(D);}voidmain(){floatd1,d2,d3,d4;floatM,M1,F,F1,F2;floatFy1,Fz1,Fy2,Fz2;floatP,P1,n,D,D1,D2,G2,G1,a,Alpha;floatMxA,MxB,MxC,MxD,MxE,MxF,MxP,MxQ,MxU,MxV,MxW;floatMyA,MyB,MyC,MyD,MyE,MyF,MyP,MyQ,MyU,MyV,MyW;floatMzA,MzB,MzC,MzD,MzE,MzF,MzP,MzQ,MzU,MzV,MzW;floatMwA,MwB,MwC,MwD,MwE,MwF,MwP,MwQ,MwU,MwV,MwW;floatPa=80*pow(10,6),E=200*pow(10,9);floatMmaxA,MmaxB,MmaxC,MmaxD,MmaxE,MmaxF,mmaxP,MmaxQ,MmaxU,MmaxV,MmaxW;floatE1,I,t,fDy,fDz,f;longints;inti,j,i1,j1,i2;printf("inputP:19.1kW\n");scanf("%f",&P);printf("inputP1:11.0kW\n");scanf("%f",&P1);printf("inputn:700r/min\n");scanf("%f",&n);printf("inputD:600mm\n");scanf("%f",&D);printf("inputD1:300mm\n");scanf("%f",&D1);printf("inputD2:150mm\n");scanf("%f",&D2);printf("inputG2:600N\n");scanf("%f",&G2);printf("inputG1:200N\n");scanf("%f",&G1);printf("inputa:600mm\n");scanf("%f",&a);printf("inputAlpha:80degrees\n");scanf("%f",&Alpha);Alpha*=Pi/180; /*单位转换*/D/=1000;D1/=1000;D2/=1000;a/=1000;M=9549*P/n;/*求扭矩和力*/M1=9549*P1/n;F=2*(M-M1)/D2;F1=2*M1/D1;F2=2*M/D;printf("M=%0.3fNm,M1=%0.3fNm,F=%0.3fN,F1=%0.3fN,F2=%0.3fN\n",M,M1,F,F1,F2);Fy2=(5*(G1+3*F1)+G2+3*F*cos(Alpha))/4;/*求支座反力*/Fy1=G2+cos(Alpha)*F-Fy2+G1+3*F1;Fz2=(3*sin(Alpha)*F-3*F2)/4;Fz1=-Fz2-3*F2+sin(Alpha)*F;printf("Fy1=%0.3fN,Fz1=%0.3fN,Fy2=%0.3fN,Fz2=%0.3fN\n\n",Fy1,Fz1,Fy2,Fz2);MxB=-M;/*各轴的扭矩*/MxD=M-M1;MxE=M1;printf("MxB=%0.3fNm,MxD=%0.3fNm,MxE=%0.3fNm\n",MxB,MxD,MxE);MzB=Fy1*a;MzD=G2*a/4+cos(Alpha)*F*3*a/4-3.0/4*(G1+3*F1)*a;MzE=-(G1+3*F1)*a;MyB=Fz1*a;MyD=3.0/4*a*F*sin(Alpha)-3.0/4*F2*3*a/3;MyE=0;printf("MzB=%0.3fNm,MzD=%0.3fNm,MzE=%0.3fNm\n",MzB,MzD,MzE);printf("MyB=%0.3fNm,MyD=%0.3fNm,MyE=%0.3fNm\n\n",MyB,MyD,MyE);MzP=1.0/2*MzB;MzQ=3.0/4*MzB+1.0/4*MzD;MzC=1.0/2*MzB+1.0/2*MzD;MzU=1.0/4*MzB+3.0/4*MzD;MzV=1.0/2*MzD+1.0/2*MzE;MzW=1.0/2*MzE;MzF=0;MyP=1.0/2*MyB;MyQ=-1.0/4*MyD-3.0/4*MyB;MyC=-1.0/2*MyD-1.0/2*MyB;MyU=-3.0/4*MyD-1.0/4*MyB;MyV=1.0/2*MyD;MyW=0;MyF=0;MxP=0;MxF=MxE;MxQ=MxC=MxU=MxB;MxV=MxW=MxE;printf("MzP=%0.3fNm,MzQ=%0.3fNm,MzC=%0.3fNm,MzU=%0.3fNm,\nMzV=%0.3fNm,MzW=%0.3fNm,MzF=%0.3fNm\n\n",MzP,MzQ,MzC,MzU,MzV,MzW,MzF);printf("MyP=%0.3fNm,MyQ=%0.3fNm,MyC=%0.3fNm,MyU=%0.3fNm,\nMyV=%0.3fNm,MyW=%0.3fNm,MyF=%0.3fNm\n\n",MyP,MyQ,MyC,MyU,MyV,MyW,MyF);printf("MxP=%0.3fNm,MxQ=%0.3fNm,MxC=%0.3fNm,MxU=%0.3fNm,\nMxV=%0.3fNm,MxW=%0.3fNm,MxF=%0.3fNm\n\n",MxP,MxQ,MxC,MxU,MxV,MxW,MxF);MmaxB=Mmax(MxB,MyB,MzB); /*求危险截面上的力矩*/MmaxD=Mmax(MxD,MyD,MzD);MmaxE=Mmax(MxE,MyE,MzE);printf("MmaxB=%0.3fNm,MmaxD=%0.3fNm,MmaxE=%0.3fNm\n",MmaxB,MmaxD,MmaxE);d1=qd(MmaxE); /*设计直径*/d2=qd(MmaxD);d3=qd(MmaxB);printf("d1=%0.3fmm,d2=%0.3fmm,d3=%0.3fmm\n",d1,d2,d3);d1=max3(d1,1.1*d2,1.21*d3);printf("可以取的安全直径为%dmm\n",int(d1/2)*2+2);i=int(d1/2)*2+2;j=int(i/1.1/2)*2+2;i1=int(j/1.1/2)*2+2;j1=int(i1/1.1/2)*2+2;printf("D1=%d,D2=%d,D3=%d,D4=%d\n\n",i,j,i1,j1);s=int(d1/2)*2+2; /*求挠度*/t=a*a*a;E1=200;I=Pi*s*s*s*s/64;fDy=(19*G2*t/24);fDy=fDy+3*t*cos(Alpha)*F/4;fDy=fDy-31*t*(G1+3*F1)/32;fDy=fDy/E1/I*1000000;printf("fDy=%0.3fmm\n",fDy);fDz=(19*t*F2/8.0-3*t*F*sin(Alpha)/4)*1000/E1/I*1000;printf("fDz=%0.3fmm\n",fDz);f=Mmax(fDy,fDz,0);printf("f=%0.2fmm\n",f);floata1[10],b[10],c[10];intnum;charh;floatko,kt,bat,d,Eyz,Ex,NO;floatomax,tmax,no,nt,w,n1;floato1=300,t1=155,T=0.10;charj2;NO=0;a1[0]=MyP;/*赋值语句*/b[0]=MzP;c[0]=MxP;a1[1]=MyB;b[1]=MzB;c[1]=MxB;

a1[2]=MyQ;b[2]=MzQ;c[2]=MxQ;a1[3]=MyD;b[3]=MzD;c[3]=MxD;a1[4]=MyV;b[4]=MzV;c[4]=MxV;a1[5]=MyE;b[5]=MzE;c[5]=MxE;a1[6]=MyW;b[6]=MzW;c[6]=MxW;a1[7]=MyF;b[7]=MzF;c[7]=MxF;for(i2=0;i2<28;i2++){scanf("%c",&h);printf("截面\n");if(h=='P'||h=='p'){/*赋值/*赋值*/ko=1.73,kt=1.41,bat=2.4,Eyz=0.81,Ex=0.76;d=j1/10.0;}elseif(h=='B'||h=='b'){num=1;ko=1.81,kt=1.62,bat=2.4,Eyz=0.81,Ex=0.76;d=i1/10.0;}elseif(h=='Q'||h=='q'){num=2;ko=1.75,kt=1.44,bat=2.4,Eyz=0.81,Ex=0.76;d=i1/10.0;}elseif(h=='D'||h=='d')

{num=3;ko=1.81,kt=1.62,bat=2.4,Eyz=0.78,Ex=0.74;d=j/10.0;}elseif(h=='V'||h=='v'){num=4;ko=1.80,kt=1.46,bat=2.4,Eyz=0.78,Ex=0.74;d=j/10.0;}elseif(h=='W'||h=='w'){num=6;ko=2.17,kt=1.59,bat=2.4,Eyz=0.81,Ex=0.0.76;d=j1/10.0;}elseif(h=='F'||h=='f'){num=7;ko=1.81,kt=1.62,bat=2.4,Eyz=0.81,Ex=0.76;d=j1/10.0;}else{printf("请你输入要测定疲劳强度的截面\n");continue;}/*检测疲劳强度过程/*检测疲劳强度过程*/no=%0.3f，nt=%0.3f，if(tmax<NO)tmax=-tmax;no=o1/ko*bat*Eyz/omax;nt=t1/(kt*tmax/Ex/bat/2.0+T*0.5*tmax);n1=no*nt/Mw(no,nt);if(tmax!=NO){if(omax!=NO){printf("omax=%0.3f，tmax=%0.3f，n=%0.3f\n",omax,tmax,no,nt,n1);

if（n1>=2）printf（"n>=2在%c截面处安全\n",h）;elseprintf（"*2在％c截面处不安全\n",h）;}else{n1=nt;printf（"omax=%0.3f，tmax=%0.3f，nt=%0.3f\n",omax,tmax,nt,n1）;if（n1>=2）printf（"n>=2在%c截面处安