1、分享一個Python寫的克林根貝爾格錐齒輪設計程序這是一個Python語言寫的克林根貝爾格錐齒輪的設計與校核程序獲取初始參數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<memin:de1+=1Re=de1/(2*sin(delta_s)b=int(round(Re*0.333,0)mnm=int(0.1*b/0.5)*0.5z1

3、=int(round(de1-b*sin(delta_s)*cos(beta_m_temp)/mnm,0)me_min=26.3*pow(K*T1*Y_FS/(z1*z1*sqrt(u*u+1)*sigmal_FP),1/3)Ri=Re-b*cos(theta_k)gamma=asin(mnm*z0/(2*r0)Rm=Re-0.5*b*cos(theta_k)zc=z1/sin(delta_s)beta_m=acos(mnm*zc/(2*Rm)Md=sqrt(Rm*Rm+r0*r0-2*Rm*r0*sin(beta_m-gamma)rb=Md/(1+z0/zc)phi_i=acos(Ri*Ri

4、+Md*Md-r0*r0)/(2*Re*Md)beta_i=atan(Ri-rb*cos(phi_i)/(rb*sin(phi_i)mni=2*Ri*cos(beta_i)/zczv_i1=z1/cos(delta_s)*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("z1:",z1)print("b:",b)print("mnm:",mnm)print("de1:",de

5、1)print("xg:",xg)初始參數獲取效果，(實際要慢慢修調齒寬系數和模數與齒寬的系數，會得到不同的效果)Python3.2(r32:88445,Feb202011,21:29:02)MSCv.150032bit(Intel)onWin32Type"copyright","credits"or"license。"formoreinformation.>>>=RESTART=>>>z1:42b:44mnm:4.0de1:244xg:-2.0193989028930326&g

6、t;>>詳細計算#-moduleimport-frommathimportpiaspifrommathimportsinassinfrommathimportcosascosfrommathimporttanastanfrommathimportatanasatanfrommathimportacosasacosfrommathimportasinasasinfrommathimportsqrtassqrtfrommathimportdegreesasdegreesfrommathimportradiansasradians#-selfdefinefunction-defanghe

7、x(x):x_result=x_int=int(x)x_minute=int(x-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# -char

8、actersblank=10# -basicelements-alpha_n=radians(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

9、_min:".rjust(blank),de1_min)print("me_min:".rjust(blank),me_min)#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("

10、;delta_b:".rjust(blank),degrees(delta_b),anghex(degrees(delta_b)print("delta_t:".rjust(blank),degrees(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(b

11、lank),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:".rjust(blank),b_try_h)#齒數mnm=4b=40beta_m_temp=radians(35)z1_try=(de1-b*sin(delta_s)*cos(beta_m_temp)/mnmprint("z1

12、_try:".rjust(blank),z1_try)#中點錐距theta_k=0Rm=Re-0.5*b*cos(theta_k)print("Rm:".rjust(blank),Rm)#假想平面齒輪齒數zc=z1/sin(delta_s)print("zc:".rjust(blank),zc)#實際中點螺旋角beta_m=acos(mnm*zc/(2*Rm)print("beta_m:".rjust(blank),degrees(beta_m),anghex(degrees(beta_m)#機床距檢驗r0=100z0=5M

13、d_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)ifMd>Md_minandMd<Md_max:print(“check:“.rjust(blank),“Mdcheckpass!")else:print("check:&

14、quot;.rjust(blank),"Mdcheckfail!")#基圓半徑rb=Md/(1+z0/zc)print("rb:".rjust(blank),rb)#Re處輔助角phi_e=acos(Re*Re+Md*Md-r0*r0)/(2*Re*Md)print("phi_e:".rjust(blank),degrees(phi_e),anghex(degrees(phi_e)#Ri處輔助角Ri=Re-b*cos(theta_k)phi_i=acos(Ri*Ri+Md*Md-r0*r0)/(2*Re*Md)print("

15、phi_i:".rjust(blank),degrees(phi_i),anghex(degrees(phi_i)#大端螺旋角beta_e=atan(Re-rb*cos(phi_e)/(rb*sin(phi_e)print("beta_e:".rjust(blank),degrees(beta_e),anghex(degrees(beta_e)#小端螺旋角beta_i=atan(Ri-rb*cos(phi_i)/(rb*sin(phi_i)print("beta_i:".rjust(blank),degrees(beta_i),anghex(d

16、egrees(beta_i)#大端法向模數mne=2*Re*cos(beta_e)/zcprint("mne:".rjust(blank),mne)#小端法向模數mni=2*Ri*cos(beta_i)/zcprint("mni:".rjust(blank),mni)#模數檢驗ifmne>mnmandmnm>mni:print("check:".rjust(blank),"modulecheckpass!")else:print("check:".rjust(blank),"

17、;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_t1*mnmEn_e1=

18、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)ifRi<RyandRy<Re:En_max=max(En_y1,En_y2)elifRy>Re:En_max=max(En_e1,En_e2)print(&quo

19、t;En_max:".rjust(blank),En_max)#法向最小齒槽寬En_i1=pi*mni/2-HwEn_i2=En_i1+4*x_t1*mnmEn_min=max(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_max<ke*En_min:<ke*En_mintcheckpass!&

20、quot;)print("check:".rjust(blank),"En_maxelse:print("check:".rjust(blank),"Pleaseredesignklingelnberggear!")# 小端二次切削檢驗ifEnmin>0.2*mnm:print("check:".rjust(blank),"En_min>0.2*mnmtcheckpass!")else:print("check:".rjust(blank),"

21、pleaseredesignklingelnberggear!")# 不產生根切的徑向變位系數zv_i1=z1/(cos(delta_s)*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)# 徑向變位系數(迭代運算)k=2*cos(beta_m)/(z2*sqrt(u*u+1)alpha_tm=atan(tan(alpha_n)/cos(beta_m)count=0x1=xg/2x

22、11=x1+0.02print("","-"*36)whileabs(x11-x1)>0,01:k11=sqrt(1+k*(ha_t-x1)*(1+k*(ha_t-x1)-cos(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*

23、k*(ha_t+x1)/(k22*k22*k22)f33=x1-f11/f22x22=x1x1=f33x11=x22count+=1print("iter:".rjust(blank),count,"x1:",x1)print("","-"*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

24、)*mnmprint("ha1:".rjust(blank),ha1)print("ha2:".rjust(blank),ha2)# 全齒高h=ha0+ha_t*mnmprint("h:".rjust(blank),h)# 當量齒輪齒數delta_n=alpha_nzv_n1=z1/(cos(delta_n)*pow(cos(beta_m),3)zv_n2=z2/(cos(delta_n)*pow(cos(beta_m),3)print("zv_n1:".rjust(blank),zv_n1)print("

25、;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=(ha0+x1*mnm-0.5*b*sin(theta_k)/t

26、an(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*sin(phi_e-lamda)YM=Md*cos(phi_e-

27、lamda)MA=sqrt(XA-XM)*(XA-XM)+(YA-YM)*(YA-YM)MB=sqrt(XB-XM)*(XB-XM)+(YB-YM)*(YB-YM)ifMA<r0+ha0*tan(alpha_n)andMB<r0+ha0*tan(alpha_n):print("check:".rjust(blank),"toolinterference,tcheckpass!")else:print("check:".rjust(blank),"toolinterference,tcheckfail!"

28、)# 小輪輪坯修正檢驗hap_t=1d_ani=mni*zv_i1+2*(hap_t+x1)*mnmpsi_ni=(pi*mni+4*mnm*(x_t1+x1*tan(alpha_n)/2*mni*zv_i1alpha_ani=acos(mni*zv_i1*cos(alpha_n)/d_ani)psi_ani=psi_ni+inv(alpha_n)-inv(alpha_ani)s_ani=psi_ani*d_aniifs_ani>=0.3*mnm:print("check:".rjust(blank),"s_ani>=0.3*mnm,tcheckpas

29、s!")else:print("check".rjust(blank),"pinionroughneedrevised!")# 中點分度圓直徑dm1=de1-b*cos(theta_k)*sin(delta_Es)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)

30、*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-2*b*cos(theta_k)*sin(delta_Es)dai2=dae2-2*b*cos(theta_k)*sin(delta_Eb)print("dai1:".rjust(blank),dai1)print(&qu

31、ot;dai2:".rjust(blank),dai2)# 實際齒寬b_ba=b*cos(theta_k)print("b_ba:".rjust(blank),b_ba)實際運行效果Python3.2(r32:88445,Feb202011,21:29:02)MSCv.150032bit(Intel)onWin32Type"copyright","credits"or"license。"formoreinformation.>>>=RESTART=>>>u:1.0d_m

32、in: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.54286461257294beta_m:35.4133536401132835,24,48gamma:5.7391704772667875,44,21Md:108.33492963933632check:Mdcheckpass!rb:97.82572235587