disp（‘*********直齒輪的設(shè)計(jì)******‘）;
disp（‘********初步計(jì)算********‘）;
disp（‘輸入已知條件‘）;
PI=pi/180;%角度變換成弧度的系數(shù)；
alpha=20*PI;%齒輪壓力角；
P1=input（‘齒輪所傳遞的功率P1=‘）;
n1=input（‘小齒輪的轉(zhuǎn)速n1=‘）;
n2=input（‘大齒輪的轉(zhuǎn)速n2=‘）;
t=input（‘輸入壽命t=‘）;
CLCL=input（‘選擇大小齒輪的材料‘）;
disp（‘查圖2.1-16至2.1-27取得相關(guān)的參數(shù)‘）;
sigma_Hlim1=input（‘小齒輪的試驗(yàn)齒輪的接觸疲勞極限sigma_Hlim1=‘）;
sigma_Hlim2=input（‘大齒輪的試驗(yàn)齒輪的接觸疲勞極限sigma_Hlim2=‘）;
sigma_Flim1=input（‘小齒輪的試驗(yàn)齒輪的彎曲疲勞極限sigma_Flim1=‘）;
sigma_Flim2=input（‘大齒輪的試驗(yàn)齒輪的彎曲疲勞極限sigma_Flim2=‘）;
T1=9550*P1/n1;%計(jì)算小齒輪的轉(zhuǎn)矩；
K=input（‘初選載荷系數(shù)K=‘）;
u=input（‘輸入齒數(shù)比u=I;u=‘）;
if?sigma_Hlim1>=sigma_Hlim2
????????sigma_Hlim=sigma_Hlim2;
????else
????????sigma_Hlim=sigma_Hlim1;
end
????SHmin=input（‘查表達(dá)式2.1-45初選接觸應(yīng)力安全系數(shù)SHmin=‘）;
????sigma_HP=sigma_Hlim/SHmin;%sigma_HP為許用接觸應(yīng)力；
?????if?sigma_Flim1>=sigma_Flim2
????????sigma_Flim=sigma_Flim2;o=2;
????else
????????sigma_Flim=sigma_Flim1;o=1;
?????end
?????SFmin=input（‘查表達(dá)式2.1-45初選彎曲應(yīng)力安全系數(shù)SFmin=‘）;
????sigma_FP=sigma_Flim/SFmin;%sigma_FP為彎曲應(yīng)力；
????SJGS=input（‘輸入所使用的設(shè)計(jì)公式：‘JCYL’（接觸應(yīng)力），‘WQYL’（彎曲應(yīng)力）‘）;
if?SJGS==‘JCYL‘
????Ja=input（‘齒面接觸強(qiáng)度計(jì)算系數(shù)Ja=‘）;
????phi_a=input（‘初選齒寬系數(shù)（參考148頁(yè)）；phi_a=‘）;
????if?sigma_Hlim1>=sigma_Hlim2
?????o=2;
????else
?????o=1;
????end
????sigma_HP=sigma_Hlim/SHmin;%sigma_HP為許用接觸應(yīng)力；
????disp（‘***按齒面接觸強(qiáng)度計(jì)算****‘）;
????LH=input（‘輸入嚙合情況：（外嚙合）‘+’；（內(nèi)嚙合）‘-’；LH=‘）;
????if?LH==‘+‘
????????amin=Ja*（u+1）*（K*T1/（phi_a*u*sigma_HP^2））^（1/3）;%amin為最小中心距；
????????a=input（‘選取中心距a>=amin;a=‘）;
????elseif?LH==‘-‘
????????amin=Ja*（u-1）*（K*T1/（phi_a*u*sigma_HP^2））^（1/3）;%amin為最小中心距；
????????a=input（‘選取中心距a>=amin;a=‘）;????????
????end
????disp（‘根據(jù)經(jīng)驗(yàn)公式：平穩(wěn)載荷m=（.007~.01）a;中等沖擊m=（.01~.015）a;‘）;
????????disp（‘較大沖擊m=（.015~.02）a;‘）;
????????m=input（‘確定齒輪模數(shù)（查表2.1-487頁(yè)）m=‘）;
elseif?SJGS==‘WQYL‘
????phi_m=input（‘初選齒寬系數(shù)（參數(shù)141頁(yè)）；phi_m=‘）;
????Z1=input（‘初選小齒輪齒數(shù)（查表2.1-39147頁(yè)）；Z1=‘）;
????YFS=input（‘初選大小齒輪復(fù)合齒形系數(shù)（查2.1-342.1-35）;YFS=‘）;
????if?sigma_Flim1>=sigma_Flim2
?????o=2;
????else
???????o=1;
????end
????SFmin=input（‘查表達(dá)式2.1-45初選接觸應(yīng)力安全系數(shù)SFmin=‘）;
????sigma_FP=sigma_lim/SFmin;%sigma_FP為彎曲應(yīng)力；
????disp（‘****按齒根彎曲強(qiáng)度計(jì)算****‘）;
????mmin=12.5*（K*T1*YFS/（phi_m*Z1^2*sigma_FP））^（1/3）;%mmin為最小模數(shù)；
????m=input（‘選取模數(shù)m（參考2.1-487頁(yè)）；m=‘）;
????disp（‘根據(jù)經(jīng)驗(yàn)公式：平穩(wěn)載荷m=（.007~.01）a;中等沖擊m=（.01~.015）a;‘）;
????disp（‘較大沖擊m=（.015~.02）a;‘）;
????a=input（‘確定中心距a=‘）;
end
Z1=2*a/（m*（u+1））;
fprintf（1‘小齒輪齒數(shù)Z1約等于%3.3f\n‘Z1）;
Z1=input（‘選取小齒輪齒數(shù)Z1=‘）;
Z2=round（u*Z1）;%Z2為大齒輪的齒數(shù)；
d1=m*Z1;%d1為小齒輪分度圓直徑；
if?SJGS==‘JCYL‘
????b=phi_a*a;%b?為齒輪寬度；
elseif?SJGS==‘WQYL‘
????b=phi_m*m;
end
disp（‘****校核齒面接觸疲勞強(qiáng)度****‘）;
disp（‘***修正載荷系數(shù)計(jì)算***‘）;
KA=input（‘使用系數(shù)（查表2.1-46161頁(yè)）;KA=‘）;
K1=input（‘計(jì)算Kv中的K1（查表2.1-48161頁(yè)）；K1=‘）;
K2=input（‘計(jì)算Kv中的K2（查表2.1-48161頁(yè)）；K2=‘）;
Ft=2000*9550*P1/（d1*n1）;%Ft為分度圓切向力；
Kv=1+（K1*b/（KA*Ft）+K2）*Z1*pi*n1*d1/1.2e7*（u^2/（u^2+1））^.5;%Kv為動(dòng)載系數(shù)；
s=input（‘s為齒輪與軸中心（齒輪在支點(diǎn)中間）或軸支點(diǎn)（齒輪在支點(diǎn)外）的距離；s=‘）;
l=input（‘l為軸支點(diǎn)間的距離；l=‘）;
if?s/l<.1
????disp（‘齒輪的結(jié)構(gòu)布局為對(duì)稱(chēng)支承‘）;
elseif?s/l>.1&&s/l<.3
????disp（‘齒輪的結(jié)構(gòu)布局為非對(duì)稱(chēng)支承‘）;
elseif?s