机械设计计算齿轮强度校核M文件

disp（‘*********直齿轮的设计******‘）;
disp（‘********初步计算********‘）;
disp（‘输入已知条件‘）;
PI=pi/180;%角度变换成弧度的系数；
alpha=20*PI;%齿轮压力角；
P1=input（‘齿轮所传递的功率P1=‘）;
n1=input（‘小齿轮的转速n1=‘）;
n2=input（‘大齿轮的转速n2=‘）;
t=input（‘输入寿命t=‘）;
CLCL=input（‘选择大小齿轮的材料‘）;
disp（‘查图2.1-16至2.1-27取得相关的参数‘）;
sigma_Hlim1=input（‘小齿轮的试验齿轮的接触疲劳极限sigma_Hlim1=‘）;
sigma_Hlim2=input（‘大齿轮的试验齿轮的接触疲劳极限sigma_Hlim2=‘）;
sigma_Flim1=input（‘小齿轮的试验齿轮的弯曲疲劳极限sigma_Flim1=‘）;
sigma_Flim2=input（‘大齿轮的试验齿轮的弯曲疲劳极限sigma_Flim2=‘）;
T1=9550*P1/n1;%计算小齿轮的转矩；
K=input（‘初选载荷系数K=‘）;
u=input（‘输入齿数比u=I;u=‘）;
if sigma_Hlim1>=sigma_Hlim2
        sigma_Hlim=sigma_Hlim2;
    else
        sigma_Hlim=sigma_Hlim1;
end
    SHmin=input（‘查表达式2.1-45初选接触应力安全系数SHmin=‘）;
    sigma_HP=sigma_Hlim/SHmin;%sigma_HP为许用接触应力；
     if sigma_Flim1>=sigma_Flim2
        sigma_Flim=sigma_Flim2;o=2;
    else
        sigma_Flim=sigma_Flim1;o=1;
     end
     SFmin=input（‘查表达式2.1-45初选弯曲应力安全系数SFmin=‘）;
    sigma_FP=sigma_Flim/SFmin;%sigma_FP为弯曲应力；
    SJGS=input（‘输入所使用的设计公式：‘JCYL’（接触应力），‘WQYL’（弯曲应力）‘）;
if SJGS==‘JCYL‘
    Ja=input（‘齿面接触强度计算系数Ja=‘）;
    phi_a=input（‘初选齿宽系数（参考148页）；phi_a=‘）;
    if sigma_Hlim1>=sigma_Hlim2
     o=2;
    else
     o=1;
    end
    sigma_HP=sigma_Hlim/SHmin;%sigma_HP为许用接触应力；
    disp（‘***按齿面接触强度计算****‘）;
    LH=input（‘输入啮合情况：（外啮合）‘+’；（内啮合）‘-’；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（‘根据经验公式：平稳载荷m=（.007~.01）a;中等冲击m=（.01~.015）a;‘）;
        disp（‘较大冲击m=（.015~.02）a;‘）;
        m=input（‘确定齿轮模数（查表2.1-487页）m=‘）;
elseif SJGS==‘WQYL‘
    phi_m=input（‘初选齿宽系数（参数141页）；phi_m=‘）;
    Z1=input（‘初选小齿轮齿数（查表2.1-39147页）；Z1=‘）;
    YFS=input（‘初选大小齿轮复合齿形系数（查2.1-342.1-35）;YFS=‘）;
    if sigma_Flim1>=sigma_Flim2
     o=2;
    else
       o=1;
    end
    SFmin=input（‘查表达式2.1-45初选接触应力安全系数SFmin=‘）;
    sigma_FP=sigma_lim/SFmin;%sigma_FP为弯曲应力；
    disp（‘****按齿根弯曲强度计算****‘）;
    mmin=12.5*（K*T1*YFS/（phi_m*Z1^2*sigma_FP））^（1/3）;%mmin为最小模数；
    m=input（‘选取模数m（参考2.1-487页）；m=‘）;
    disp（‘根据经验公式：平稳载荷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‘小齿轮齿数Z1约等于%3.3f\n‘Z1）;
Z1=input（‘选取小齿轮齿数Z1=‘）;
Z2=round（u*Z1）;%Z2为大齿轮的齿数；
d1=m*Z1;%d1为小齿轮分度圆直径；
if SJGS==‘JCYL‘
    b=phi_a*a;%b 为齿轮宽度；
elseif SJGS==‘WQYL‘
    b=phi_m*m;
end
disp（‘****校核齿面接触疲劳强度****‘）;
disp（‘***修正载荷系数计算***‘）;
KA=input（‘使用系数（查表2.1-46161页）;KA=‘）;
K1=input（‘计算Kv中的K1（查表2.1-48161页）；K1=‘）;
K2=input（‘计算Kv中的K2（查表2.1-48161页）；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为动载系数；
s=input（‘s为齿轮与轴中心（齿轮在支点中间）或轴支点（齿轮在支点外）的距离；s=‘）;
l=input（‘l为轴支点间的距离；l=‘）;
if s/l<.1
    disp（‘齿轮的结构布局为对称支承‘）;
elseif s/l>.1&&s/l<.3
    disp（‘齿轮的结构布局为非对称支承‘）;
elseif s