【正文】 ne i 120e6 define E 150e9 int main() {int angle,d,D。 double temps,tempc,Fy,Fz,Mx,e,l,l1,l2,l3,l4,F,FAy,FAz,FBy,FBz,W,Mx1,My1,Mz1,Mx2,My2,Mz2,DD,g,f,dd,h,b,A,s1,s2,s3,t2,t3,a,sr2,sr3,r,s11,s22,s33,t22,t33,sr22,sr33,tmax,tmin,R,tm,ta,nt,T,Kt,Et,B,pt,Q1,Q11,Q12,Q13,Q14,Q15,Q16,EI1,EI2,EI3,Q2,Q21,Q22,Q23,Q24,Q25,GI1,B1,u,G。 printf(enter the data:F,W,l1,l2,l3,e,angle\n)。 scanf(%lf,%lf,%lf,%lf,%lf,%lf,%d,amp。F,amp。W,amp。l1,amp。l2,amp。l3,amp。e,amp。angle)。 temps=sin(angle*Pi/180)。 tempc=cos(angle*Pi/180)。 Fy=F*temps。 Fz=F*tempc。 Mx=Fz*e。 printf(Fy=%lf,Fz=%lf,Mx=%lf\n,Fy,Fz,Mx)。 FBz=Fz*l1/(l1+l2)。 FBy=(Fy*l1W*(l1+l2+l3))/(l1+l2)。 FAz=l2*Fz/(l1+l2)。 FAy=(Fy*l2+W*l3)/(l1+l2)。 printf(FAy=%lf,FAz=%lf,FBy=%lf,FBz=%lf\n,FAy,FAz,FBy,FBz)。 l=*e。 l4=*l。 Mx1=FAz*e+(FzFAz)*e。 My1=FAz*(l1+l4)Fz*l4。 Mz1=FAy*(l1+l4)Fy*l4。 f=sqrt(Mx1*Mx1+My1*My1+Mz1*Mz1)。 g=32*f/(Pi*i)。 DD=pow(g,)。 My2=l1*FAz。 Mx2=e*FAz。 Mz2=FAy*l1。 f=sqrt(Mx2*Mx2+My2*My2+Mz2*Mz2)。 g=32*f/(Pi*i)。 dd=pow(g,)。 printf(D=%.4lfm,d=%.4lfm\n,DD,dd)。 D=(int)(1000*DD)。 24 d=(int)(1000*dd)。 if(d%2==0) d=d+2。 else d=d+1。 if(D%2==0) D=D+2。 else D=D+1。 printf(故 D 取 %dmm,d 取 %dmm \n,D,d)。 a=。 r=。 h=*(double)D/1000。 b=2*h/3。 A=h*b。 s1=FAy/A+FAz*e*6/(h*h*b)+FAy*6*(l1l4)/(b*b*h)。 printf(s1=%.2lfPa,s1)。 s2=FAy/A+FAy*(l1l4)*6/(b*b*h)。 t2=FAz*(l1l4)/(a*h*b*b)。 sr2=sqrt(s2*s2+4*t2*t2)。 printf(sr2=%lfPa,sr2)。 t3=r*t2。 s3=FAy/A+FAz*e*6/(h*h*b)。 sr3=sqrt(s3*s3+4*t3*t3)。 printf(sr3=%lfPa\n,sr3)。 if(s1=iamp。amp。sr2=iamp。amp。sr3=i) printf(左曲軸臂安全 \n)。 s11=(W+FBy)/A+Mx1*6/(h*h*b)+(FAy*(l1+l4)Fy*l4)*6/(h*b*b)。 printf(s11=%lfPa,s11)。 s22=(W+FBy)/A+(FAy*(l1+l4)Fy*l4)*6/(h*b*b)。 t22=(FAz*(l1+l4)Fz*l4)/(a*h*b*b)。 sr22=sqrt(s22*s22+4*t22*t22)。 printf(sr22=%lfPa,sr22)。 s33=(W+FBy)/A+Mx1*6/(h*h*b)。 t33=r*t22。 sr33=sqrt(s33*s33+4*t33*t33)。 printf(sr33=%lfPa\n,sr33)。 if(s11=iamp。amp。sr22=iamp。amp。sr33=i) printf(右曲軸臂安全 \n)。 Kt=。 Et=。 B=。 pt=。 T=160e6。 tmax=Mx*16/(Pi*((double)D/1000)*((double)D/1000)*((double)D/1000))。 tmin=0。 25 R=tmin/tmax。 tm=(tmax+tmin)/2。 ta=(tmaxtmin)/2。 nt=T/(Kt/(Et*B)*ta+pt*tm)。 printf(nt=%lf\n,nt)。 if(nt=n) printf(安全 \n)。 EI1=E*Pi*((double)D/1000)*((double)D/1000)*((double)D/1000)*((double)D/1000)/64。 EI2=E*h*b*b*b/12。 EI3=E*Pi*((double)d/1000)*((double)d/1000)*((double)d/1000)*((double)d/1000)/64。 Q11=(l1l4)*FAy*(l1l4)/2*(l1+2*l4+3*l2)/(3*(l1+l2))/EI1。 Q12=FAy*(l1l4)*e*(l2+l4)/(l1+l2)/EI2。 Q13=(FAy*l4*l4/2*(3*l2+l4)/(3*(l1+l2))+FAy*(l1l4)*l4*(2*l2+l4)/(2*(l1+l2))+(FyFAy)*l4*l4/2*(3*l2l4)/(3*(l1+l2))+(FAy*(l1+l4)Fy*l4)*l4*(2*l2l4)/(2*(l1+l2)))/EI3。 Q14=(l2l4)/(l1+l2)*e*(FAy*(l1+l4)Fy*l4)/EI2。 Q15=(l2l4)/(l1+l2)*(l2l4)/2*(FAy*(l1+l4+(l2l4)/3)Fy*(l4+(l2l4)/3))/EI1。 Q16=e/((l1+l2)*E*h*b)*(W+FByFAy)。 Q1=Q11+Q12+Q13+Q14+Q15+Q16。 printf(A 點處轉角θ z=%lf\n,Q1)。 u=。 G=E/(2*(1+u))。 B1=。 GI1=G*B1*h*b*b*b。 Q21=FAz*(l1l4)*(l1l4)/2*(l1+3*l2+2*l4)/(3*(l1+l2))/EI1。 Q22=FAz*(l1l4)*e*(l2+l4)/(l1+l2)/GI1。 Q23=(FAz*l4*l4/2*(3*l2+l4)/(3*(l1+l2))+FAz*(l1l4)*l4*(2*l2+l4)/(2*(l1+l2))+(FAz+Fz)*l4*l4/2*(3*l2l4)/(3*(l1+l2))+(FAz*(l1+l4)l4*Fz)*l4*(2*l2l4)/(2*(l1+l2)))/EI3。 Q24=(l2l4)/(l1+l2)*e*(FAz*(l1+l4)l4*Fz)/GI1。 Q25=(l2l4)/(l1+l2)*(l2l4)*(FAz*(l1+l4)l4*Fz)/(3*EI1)。 Q2=Q21+Q22+Q23+Q24+Q25。 printf(A 點處轉角θ y=%lf\n,Q2)。 return 0。 } 26 三、計算輸出結果