1、五種傳動軸靜強度變形計算設計題目：傳動軸地 材料為優質碳素鋼 V牌號45),許用應力o=80MPa,經高頻淬火處理軸地表面，鍵槽均為端銑加工,E=210GPa已知數據P/KWP1/KWn/rpmD/mmD1/mmD2fmma/mma /.29.48.13008005002005000傳動軸力學簡圖1傳動軸零件圖aaaa%輸入已知數據sigma=80。E=210000。P=input('請輸入大帶輪傳動地功率P= <KW)'>。P1=input('請輸入小帶輪傳動地功率P1= <KW)'>。n=input('請輸入小帶輪地轉速n=

2、(rpm>'>。D=input('請輸入大帶輪直徑D= <mm)'>。D仁input('請輸入小帶輪直徑 D1= <mm)'>。D2=input('請輸入齒輪直徑 D2= <mm)'>。a=input('請輸入 a= (mm>'>。alfa=input('請輸入 a ='>。%計算各輪受力并輸出F2=2*9.549*10A6*P/n/D。fprintf('大帶輪 D 上作用地水平力：3*F2=%3.3f(N>n',3*F

3、2>。m=9.549*10A6*P/n 。fprintf('大帶輪 D 上作用地力偶：m=%3.3f(Nmm>n',m>。F1=2*9.549*10A6*P1 /n/D1。fprintf('小帶輪 D1 上作用地鉛垂力：3*F%3.3f(N>n',3*F1>。m1=9.549*10A6*P1 /n。fprintf('小帶輪 D1 上作用地力偶：m1=%3.3f(Nmm>n',m1>。F=2*(m-m1>/D2。fprintf('齒輪 D2 上作用地水平力： F*sin a =%3.3f(N&

4、gt;n',F*sin(alfa>>。 fprintf('齒輪 D2 上作用地鉛垂力： F*cos a =%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 。 %分段數組 <不同結構數據不同)MPQ仁2 1 a 0。齒輪處加單位力地載荷矩陣<不同結構數據不同)MPQH=input(' 請輸入水平面內梁地載荷矩陣MP

5、QH= '。MPQV=input(' 請輸入鉛垂面內梁地載荷矩陣MPQV= '。MT=input(' 請輸入外力偶矩載荷矩陣 MT= '。 %根據梁地類型選擇相應程序計算水平面、鉛垂面玩具、單位力彎矩 string=input(' 請輸入梁地類型 可選擇 J,YW,LDW)= ','s'。 if strcmp(string,'LDW'=1L1=input(' 請輸入梁左端外伸段地長度L1= '。L2=input(' 請輸入梁右端外伸段地長度L2= '。XQMH=QMLDW(x

6、,L1,L2,MPQH。 XQMV=QMLDW(x,L1,L2,MPQV。 XQM1=QMLDW(x,L1,L2,MPQ1。else if strcmp(string,'YW'=1L2=i nput('請輸入梁右端外伸段地長度L2='XQMH=QMYW(x,L2,MPQH。XQMV=QMYW(x,L2,MPQV。XQM1=QMYW(x,L2,MPQ1。else strcmp(string,'J'=1XQMH=QMDJ(x,MPQH。XQMV=QMDJ(x,MPQV。XQM1=QMDJ(x,MPQ1。endend%計算各截面扭矩XT=TTT(x,M

7、T。%水平面、鉛垂面彎矩、單位力彎矩及各截面扭矩取值 x1=XQMH(:,1。mh=XQMH(:,3。mv=XQMV(:,3。m1=XQM1(:,3。TT=XT(:,2。%計算合成彎矩ix,dummy=size(XQMH。for i=1:ixM(i=sqrt(mh(iA2+mv(iA2。end %畫出水平面彎矩圖、鉛垂面彎矩圖、合成彎矩圖、扭矩圖 subplot(5,1,1plot(x1,mh> 。grid。title(' 水平面內地彎矩圖 '>subplot(5,1,2>plot(x1,mv>。grid。title(' 鉛垂面內地彎矩圖 

8、9;>。subplot(5,1,3>plot(x1,M> 。grid。title(' 合成彎矩圖 '>。subplot(5,1,4> 。plot(x1,TT>grid。title('扭矩圖'>。subplot(5,1,5>。plot(x1,m1>grid。title('單位力作用彎矩圖>。%利用第四強度理論設計軸徑for i=1:ixMM(i>=sqrt(M(i>A2+0.75*TT(i>A2>。 endM41=max(MM(1:25>> 。M3=max(MM

9、(25:(50+12>>> 。M2=max(MM(50+12>:(50+50+25>>> 。M1=max(MM(50+50+25>:(50+50+50+25>>> 。M42=max(MM(50+50+50+25>:(50+50+50+50>>> 。M4=max(M41,M42> 。d4=(32*M4/(pi*sigma>>A(1 /3>。d3=(32*M3/(pi*sigma>>A(1 /3>。d2=(32*M2/(pi*sigma>>A(1 /3&g

10、t;。d1=(32*M1/(pi*sigma>>A(1 /3>。d4 = %3.3fmmn',d4>。 d3 仝 %3.3fmmn',d3>。 d2= %3.3fmmn',d2> 。 di = %3.3fmmn',d1> 。fprintf(' 有第四強度理論：軸徑 fprintf(' 有第四強度理論：軸徑 fprintf(' 有第四強度理論：軸徑 fprintf(' 有第四強度理論：軸徑%計算齒輪處軸地撓度I=pi*d1A4/64。fh=0。fv=0。dummy,m=size(x>。

11、for i=1:m-1dx=(x(i+1>-x(i>>/50。for j=1:50 fh=fh+mh(i-1>*50+j>*m1(i-1>*50+j>*dx 。 fv=fv+mv(i-1>*50+j>*m1(i-1>*50+j>*dx 。endendfv=fv/(E*I> 。 fh=fh/(E*I> 。fprintf(' 齒輪處軸地豎直撓度 fh=%3.3fmmn',fh> 。 fprintf(' 齒輪處軸地豎直撓度 fv=%3.3fmmn',fv> 。rintf('

12、;齒輪處軸地合成撓度 f=%3.3fmm門'網11仲人2+22>>。運行結果 :請輸入大帶輪傳動地功率P= <KW)29.4請輸入小帶輪傳動地功率P1= <KW) 8 . 1請輸入小帶輪地轉速 n= (rpm>300 請輸入大帶輪直徑 D= <mm) 800 請輸入小帶輪直徑 D1= <mm) 500 請輸入齒輪直徑 D2= <mm)200 請輸入 a= (mm>500 請輸入a = 0大帶輪 D 上作用地水平力： 3*F2=7018.515(N> 大帶輪 D 上作用地力偶： m=935802.000(Nmm> 小帶輪

13、D1 上作用地鉛垂力： 3*F1=3093.876(N> 小帶輪 D1 上作用地力偶： m1=257823.000(Nmm> 齒輪 D2 上作用地水平力： F*sina =0.000(N> 齒輪 D2 上作用地鉛垂力： F*cosa =6779.790(N> 齒輪 D2上作用地力偶： m-m仁677979.000(Nmm> 請輸入水平面內梁地載荷矩陣 MPQH= 2 7018.515 4*a 0 請輸入鉛垂面內梁地載荷矩陣 MPQV= 2 3093.876 3*a 0 2 6779.79 a 0請輸入外力偶矩載荷矩陣MT= 677979 a257823 3*a935802 5*a請輸入梁地類