千斤顶 福州大学.doc

六. 附表及计算机程序 计算机程序： #includestdio.h #includemath.h # define pi 3.141592 void main() { double A(int i,int j); double F,H1,Fp; double xcgm=(4.8e8)/4, xcgmb=60E6, xtao=40E6; /* 螺杆许用应力 ，螺母许用应力，许用剪应力 */ double cgm, cgmb, tao; /* 螺杆应力,螺母应力,剪应力 */ /* 螺杆参数 */ double d1,d,D,d2,P,fy,fyv; /*依次为：外螺纹小径,外螺纹大径,内螺纹大径,螺纹中径,螺距P,螺纹升角,当量摩擦角*/ double l, Fcr , efficiency , f=0.08 ; /*依次为:工作长度,临界载荷,传动总效率,摩擦系数 */ /* 螺母参数 */ double H,u,xp=20e6; /*依次为:螺母的高度,旋合圈数,许用压力*/ /* 手柄参数 */ double L,Lp,dp; /*依次为:手柄总长度，有效长度，直径*/ /* 中间变量 */ double y,h1=0,DD,DD1,h,b,I,T1,T2,D0,d0; int i; /* printf(please input F(KN),H1(mm),Fp(N)=); scanf(%lf,%lf,%lf,F,H1,Fp); printf(%lf,%lf,%lf\n,F,H1,Fp); F=F*1000; H1=H1*0.001; */ F=70000,H1=0.27,Fp=200; d1=1.25*sqrt(F/xcgm)*1000; i=0; do{ for(;d1=A(i,0);i++) ; d1=0.001*A(i,0);/*螺杆小径*/ d=0.001*A(i,1); D=0.001*A(i,2); d2=0.001*A(i,3); P=0.001*A(i,4); y=P/(pi*d2); fy=atan(y); /*螺纹升角*/ y=f/cos(15*pi/180); fyv=atan(y); /*当量摩擦角*/ if(i==17) printf(please design again !); h=0.5*P; H=F*P/(pi*d2*h*xp); H=1000*H; H=(int)H+1; H=H/1000; /*螺母高度*/ u=H/P; /*旋合圈数*/ b=0.634*P; tao=F/(pi*D*b*u);/*剪应力*/ cgmb=3*F*h/(pi*D*b*b*u); /*弯曲应力*/ DD=1.6*d; DD1=0.6*d; T1=F*d2*tan(fy+fyv)/2; D0=DD-1.5*0.001; d0=DD1+1.5*0.001; T2=0.12*F*(D0*D0*D0-d0*d0*d0)/(3*(D0*D0-d0*d0)); Lp=(T1+T2)/Fp; L=Lp+1.6*d/2; /*手柄总长度*/ dp=Fp*Lp/(0.1*150*1e6); dp=1000*pow(dp,1.0/3); dp=(int)dp+1; dp=dp/1000; /*手柄直径*/ d1=d1*1000; I=pi*d1*d1*d1*d1/64; d1=d1/1000; h1=1.9*dp; l=h1+H1+H/2; Fcr=pi*pi*2.06*1e11*I/(2*l*2*l)/1e12; /*临界应力*/ i++; } while(taoxtao||cgmbxcgmb||Fcr/F3.5||fyfyv||u10) ; efficiency=F*P/(2*pi*(T1+T2));/*传动总效率*/ printf(d1=%4.2lfmm , ,d1*1000); printf(H=%4.2lfmm , ,H*1000); printf(u=%4.3lf , \n,u); printf(L=%6.2lfmm, ,L*1000); printf