五种传动轴静强度变形计算.doc

五种传动轴静强度 变形计算 设计题目: 传动轴的材料为优质碳素钢（牌号45），许用应力[σ]=80MPa， 经高频淬火处理。轴的表面，键槽均为端铣加工，E=210GPa。 已知数据 P/KW P1/KW n/rpm D/mm D1/mm D2/mm a/mm α/。 29.4 8.1 300 800 500 200 500 0 传动轴力学简图 传动轴零件图 %输入已知数据 sigma=80; E=210000; P=input(请输入大带轮传动的功率 P= （KW）); P1=input(请输入小带轮传动的功率 P1= （KW）); n=input(请输入小带轮的转速 n= (rpm)); D=input(请输入大带轮直径 D= （mm）); D1=input(请输入小带轮直径 D1= （mm）); D2=input(请输入齿轮直径 D2= （mm）); a=input(请输入 a= (mm)); alfa=input(请输入 α= ); %计算各轮受力并输出 F2=2*9.549*10^6*P/n/D; fprintf(大带轮D上作用的水平力：3*F2=%3.3f(N)\n,3*F2); m=9.549*10^6*P/n; fprintf(大带轮D上作用的力偶：m=%3.3f(Nmm)\n,m); F1=2*9.549*10^6*P1/n/D1; fprintf(小带轮D1上作用的铅垂力：3*F1=%3.3f(N)\n,3*F1); m1=9.549*10^6*P1/n; fprintf(小带轮D1上作用的力偶：m1=%3.3f(Nmm)\n,m1); F=2*(m-m1)/D2; fprintf(齿轮D2上作用的水平力：F*sinα=%3.3f(N)\n,F*sin(alfa)); fprintf(齿轮D2上作用的铅垂力：F*cosα=%3.3f(N)\n,F*cos(alfa)); fprintf(齿轮D2上作用的力偶：m-m1=%3.3f(Nmm)\n,m-m1); %根据各轮受力确定水平面、铅垂面、外力偶矩的载荷矩阵 x=[0 a 3*a 4*a 5*a];%分段数组（不同结构数据不同） MPQ1=[2 1 a 0];%齿轮处加单位力的载荷矩阵（不同结构数据不同） MPQH=input(请输入水平面内梁的载荷矩阵 MPQH= ); MPQV=input(请输入铅垂面内梁的载荷矩阵 MPQV= ); MT=input(请输入外力偶矩载荷矩阵MT= ); %根据梁的类型选择相应程序计算水平面、铅垂面玩具、单位力弯矩 string=input(请输入梁的类型（可选择J,YW,LDW）= ,s); if strcmp(string,LDW)==1 L1=input(请输入梁左端外伸段的长度L1= ); L2=input(请输入梁右端外伸段的长度L2= ); XQMH=QMLDW(x,L1,L2,MPQH); XQMV=QMLDW(x,L1,L2,MPQV); XQM1=QMLDW(x,L1,L2,MPQ1); else if strcmp(string,YW)==1 L2=input(请输入梁右端外伸段的长度L2= ) XQMH=QMYW(x,L2,MPQH); XQMV=QMYW(x,L2,MPQV); XQM1=QMYW(x,L2,MPQ1); else strcmp(string,J)==1 XQMH=QMDJ(x,MPQH); XQMV=QMDJ(x,MPQV); XQM1=QMDJ(x,MPQ1); end end %计算各截面扭矩 XT=TTT(x,MT); %水平面、铅垂面弯矩、单位力弯矩及各截面扭矩取值 x1=XQMH(:,1); mh=XQMH(:,3); mv=XQMV(:,3); m1=XQM1(:,3); TT=XT(:,2); %计算合成弯矩 [ix,dummy]=size(XQMH); for i=1:ix M(i)=sqrt(mh(i)^2+mv(i)^2); end %画出水平面弯矩图、铅垂面弯矩图、合成弯矩图、扭矩图 subplot(5,1,1) plot(x1,mh); grid; title(水平面内的弯矩图) subplot(5,1,2) plot(x1,mv);