分享一個Python寫的克林根貝爾格錐齒輪設(shè)計程序

1、分享一個Python寫的克林根貝爾格錐齒輪設(shè)計程序這是一個Python語言寫的克林根貝爾格錐齒輪的設(shè)計與校核程序獲取初始參數(shù)frommathimportpiaspifrommathimportsinassinfrommathimportcosascosfrommathimporttanastanfrommathimportatanasatanfrommathimportacosasacosfrommathimportasinasasinfrommathimportsqrtassqrtfrommathimportdegreesasdegreesfrommathimportradiansasrad

2、ianstheta_t=radians(135)alpha_n=radians(20)thetak=0de1=150beta_m_temp=radians(35)u=1delta_s=atan(sin(theta_t)/(u+cos(theta_t)mnm=0.5memin=1z1=1b=1hat=1K=1.2SF=1.4YFS=4.35sigmal_FP=1000/S_FT1=1200/ur0=135z0=5whilemnm=RESTART=z1:42b:44mnm:4.0de1:244xg:-2.0193989028930326詳細(xì)計算#-moduleimport-frommathimpo

3、rtpiaspifrommathimportsinassinfrommathimportcosascosfrommathimporttanastanfrommathimportatanasatanfrommathimportacosasacosfrommathimportasinasasinfrommathimportsqrtassqrtfrommathimportdegreesasdegreesfrommathimportradiansasradians#-selfdefinefunction-defanghex(x):x_result=x_int=int(x)x_minute=int(x-

4、x_int)*60)x_second=int(round(x-x_int-x_minute/60.0)*3600,0)ifxsecond60:xminute=xminute+1xsecond=0ifxminute=60:xint=xint+1xminute=0x_result.append(x_int)x_result.append(x_minute)x_result.append(x_second)returnxresultdefinv(x):result=tan(x)-xreturnresult# -charactersblank=10# -basicelements-alpha_n=ra

5、dians(20)z1=43z2=z1u=z2/z1T1=1200/u#K=1.2S_H=1.1S_F=1.4Y_FS=4.35sigmal_HP=1500/S_Hsigmal_FP=1000/S_Fde1_min=1636*pow(K*T1/(u*sigmal_HP*sigmal_HP),1/3)me_min=26.3*pow(K*T1*Y_FS/(z1*z1*sqrt(u*u+1)*sigmal_FP),1/3)print(u:.rjust(blank),u)print(d_min:.rjust(blank),de1_min)print(me_min:.rjust(blank),me_mi

6、n)#theta_t=radians(135)delta_s=atan(sin(theta_t)/(u+cos(theta_t)delta_b=atan(u*sin(theta_t)/(1+u*cos(theta_t)delta_t=delta_s+delta_bprint(delta_s:.rjust(blank),degrees(delta_s),anghex(degrees(delta_s)print(delta_b:.rjust(blank),degrees(delta_b),anghex(degrees(delta_b)print(delta_t:.rjust(blank),degr

7、ees(delta_t),anghex(degrees(delta_t)# 外錐距de1=248de2=u*de1Re=de1/(2*sin(delta_s)print(Re:.rjust(blank),Re)print(de1:.rjust(blank),de1)print(de2:.rjust(blank),de2)# 齒寬b_try_q=Re*0.2b_try_l=Re*0.286b_try_h=Re*0.333print(b_try_q:.rjust(blank),b_try_q)print(b_try_l:.rjust(blank),b_try_l)print(b_try_h:.rj

8、ust(blank),b_try_h)#齒數(shù)mnm=4b=40beta_m_temp=radians(35)z1_try=(de1-b*sin(delta_s)*cos(beta_m_temp)/mnmprint(z1_try:.rjust(blank),z1_try)#中點錐距theta_k=0Rm=Re-0.5*b*cos(theta_k)print(Rm:.rjust(blank),Rm)#假想平面齒輪齒數(shù)zc=z1/sin(delta_s)print(zc:.rjust(blank),zc)#實際中點螺旋角beta_m=acos(mnm*zc/(2*Rm)print(beta_m:.r

9、just(blank),degrees(beta_m),anghex(degrees(beta_m)#機(jī)床距檢驗r0=100z0=5Md_min=0Md_max=400gamma=asin(mnm*z0/(2*r0)Md=sqrt(Rm*Rm+r0*r0-2*Rm*r0*sin(beta_m-gamma)print(gamma:.rjust(blank),degrees(gamma),anghex(degrees(gamma)print(Md:.rjust(blank),Md)ifMdMd_minandMdmnmandmnmmni:print(check:.rjust(blank),modul

10、echeckpass!)else:print(check:.rjust(blank),modulecheckfail!)#法向齒槽最大處的錐距Ry=sqrt(pow(zc-z0)/(zc+z0),2)*Md*Md+r0*r0)print(Ry:.rjust(blank),Ry)#輔助值x_t1=0x_t2=-x_t1ha_t=1.0ha0=1.25*ha_t*mnmHw=2*(x_t1*mnm+ha0*tan(alpha_n)print(Hw:.rjust(blank),Hw)#法向最大齒槽寬zp=z2/sin(delta_b)En_y1=pi*rb/zp-HwEn_y2=En_y1+4*x_

11、t1*mnmEn_e1=pi*mne/2-HwEn_e2=En_e1+4*x_t1*mnmprint(En_y1:.rjust(blank),En_y1)print(En_y2:.rjust(blank),En_y2)print(En_e1:.rjust(blank),En_e1)print(En_e2:.rjust(blank),En_e2)ifRiRyandRyRe:En_max=max(En_e1,En_e2)print(En_max:.rjust(blank),En_max)#法向最小齒槽寬En_i1=pi*mni/2-HwEn_i2=En_i1+4*x_t1*mnmEn_min=ma

12、x(En_i1,En_i2)print(En_i1:.rjust(blank),En_i1)print(En_i2:.rjust(blank),En_i2)print(En_min:.rjust(blank),En_min)#大端槽底檢驗ke=3.8ifEn_maxke*En_min:0.2*mnm:print(check:.rjust(blank),En_min0.2*mnmtcheckpass!)else:print(check:.rjust(blank),pleaseredesignklingelnberggear!)# 不產(chǎn)生根切的徑向變位系數(shù)zv_i1=z1/(cos(delta_s

13、)*pow(cos(beta_i),3)xg=1.1*ha_t-(mni*zv_i1*sin(alpha_n)*sin(alpha_n)+b*sin(theta_k)/(2*mnm)print(xg:.rjust(blank),xg)# 徑向變位系數(shù)(迭代運(yùn)算)k=2*cos(beta_m)/(z2*sqrt(u*u+1)alpha_tm=atan(tan(alpha_n)/cos(beta_m)count=0x1=xg/2x11=x1+0.02print(,-*36)whileabs(x11-x1)0,01:k11=sqrt(1+k*(ha_t-x1)*(1+k*(ha_t-x1)-cos(

14、alpha_tm)*cos(alpha_tm)k22=sqrt(1+u*u*k*(ha_t+x1)*(1+u*u*k*(ha_t+x1)-cos(alpha_tm)*cos(alpha_tm)f11=u*u/k11-1/k22-(u*u-1)/sin(alpha_tm)f22=u*u*k*(1+k*(ha_t-x1)/(k11*k11*k11)+u*u*k*(1+u*u*k*(ha_t+x1)/(k22*k22*k22)f33=x1-f11/f22x22=x1x1=f33x11=x22count+=1print(iter:.rjust(blank),count,x1:,x1)print(,-*

15、36)x2=-x1print(x1:.rjust(blank),x1)print(x2:.rjust(blank),x2)# 齒頂高x1=round(x1,3)x2=-x1ha1=(ha_t+x1)*mnmha2=(ha_t+x2)*mnmprint(ha1:.rjust(blank),ha1)print(ha2:.rjust(blank),ha2)# 全齒高h(yuǎn)=ha0+ha_t*mnmprint(h:.rjust(blank),h)# 當(dāng)量齒輪齒數(shù)delta_n=alpha_nzv_n1=z1/(cos(delta_n)*pow(cos(beta_m),3)zv_n2=z2/(cos(del

16、ta_n)*pow(cos(beta_m),3)print(zv_n1:.rjust(blank),zv_n1)print(zv_n2:.rjust(blank),zv_n2)# 工藝分錐角delta_Es=delta_s-theta_kdelta_Eb=delta_b+theta_kprint(delta_Es:.rjust(blank),degrees(delta_Es),anghex(degrees(delta_Es)print(delta_Eb:.rjust(blank),degrees(delta_Eb),anghex(degrees(delta_Eb)# 刀盤干涉檢驗lamda=(

17、ha0+x1*mnm-0.5*b*sin(theta_k)/tan(alpha_n)+ha0*tan(alpha_n)/Redelta_h=Rm*tan(theta_k)*cos(delta_Eb)XA=sqrt(2*h*(Re+tan(delta_Eb)+ha2-delta_h)-(h/cos(delta_Eb)*(h/cos(delta_Eb)YA=Re-h*tan(delta_Eb)XB=sqrt(2*h*(Ri+tan(delta_Eb)+ha2-delta_h)-(h/cos(delta_Eb)*(h/cos(delta_Eb)YB=Ri-h*tan(delta_Eb)XM=Md*s

18、in(phi_e-lamda)YM=Md*cos(phi_e-lamda)MA=sqrt(XA-XM)*(XA-XM)+(YA-YM)*(YA-YM)MB=sqrt(XB-XM)*(XB-XM)+(YB-YM)*(YB-YM)ifMAr0+ha0*tan(alpha_n)andMB=0.3*mnm:print(check:.rjust(blank),s_ani=0.3*mnm,tcheckpass!)else:print(check.rjust(blank),pinionroughneedrevised!)# 中點分度圓直徑dm1=de1-b*cos(theta_k)*sin(delta_Es

19、)dm2=de2-b*cos(theta_k)*sin(delta_Eb)print(dm1:.rjust(blank),dm1)print(dm2:.rjust(blank),dm2)# 大端齒頂圓直徑dae1=de1+2*ha1*cos(delta_Es)-b*sin(theta_k)*cos(delta_Es)dae2=de2+2*ha2*cos(delta_Eb-pi/2)-b*sin(theta_k)*cos(delta_Eb)print(dae1:.rjust(blank),dae1)print(dae2:.rjust(blank),dae2)# 小端齒頂圓直徑dai1=dae1-

20、2*b*cos(theta_k)*sin(delta_Es)dai2=dae2-2*b*cos(theta_k)*sin(delta_Eb)print(dai1:.rjust(blank),dai1)print(dai2:.rjust(blank),dai2)# 實際齒寬b_ba=b*cos(theta_k)print(b_ba:.rjust(blank),b_ba)實際運(yùn)行效果Python3.2(r32:88445,Feb202011,21:29:02)MSCv.150032bit(Intel)onWin32Typecopyright,creditsorlicense。formoreinfo

21、rmation.=RESTART=u:1.0d_min:150.23521659011394me_min:3.9366944420085543delta_s:67.567,30,0delta_b:67.567,30,0delta_t:135.0135,0,0Re:134.21663283625685de1:248de2:248.0b_try_q:26.84332656725137b_try_l:38.38595699116946b_try_h:44.69413873447353z1_try:43.2194486685837Rm:114.21663283625685zc:46.542864612

22、57294beta_m:35.4133536401132835,24,48gamma:5.7391704772667875,44,21Md:108.33492963933632check:Mdcheckpass!rb:97.82572235587946phi_e:47.2216566589780447,13,18phi_i:68.5950080434168868,35,42beta_e:43.34789929041891643,20,52beta_i:32.7192686445736532,43,9mne:4.1940788707287355mni:3.4062016138360653check:modulecheckpass!Ry:132.7560688043775HW:3.6397023426620234En_y1:2.96342733669005En_y2:2.96342733669005En_e1:2.9