齿轮接触疲劳强度分析程序_机械设计中使用_matlab

function [sH] = sH( ) %齿轮接触疲劳强度计算 T=73800*1000;%%%（N*mm）转矩 n1=20.83*98.47,%（r/min）小齿轮转速 mn=8;%%%（mm）法向模数 z1=36;%%%小齿轮齿数 z2=121;%%%大齿轮齿数 bb=11;%%%（度）螺旋角；直齿轮设为0 an=20;%%%(度）法面压力角,直齿即是压力角 b=160;%%%（mm）齿宽 aa=650;%%%实际中心距 ha=1;%%%齿顶高系数 x1=0.418;%%%小齿轮变位系数 x2=0.933;%%%大齿轮变位系数 %斜齿轮/直齿轮选择 if bb>=1 bbb=bb*3.14/180;%（rad）螺旋角弧度表示 else bbb=0;%（rad）螺旋角弧度表示 end %斜齿轮/直齿轮选择 d1=(mn/cos(bbb))*z1;%（mm）小齿轮分度圆直径 d2=(mn/cos(bbb))*z2;%（mm）大齿轮分度圆直径 Ft=2*T/d1;%（N）圆周力 kA=1.75*1.1;%%%使用系数 v=d1/2*n1*2*3.14/(60*1000);%%%（m/s）节圆线速度 kv=1.3;%%%动载系数 a=kA*Ft/b;%（N/mm）确定齿间载荷分配系数ka kb=1.2;%%%齿向载荷分配系数 ka=1.2;%%%齿间载荷分配系数 k=kA*kv*ka*kb;%载荷系数 %斜齿轮端面向重合度 if bb>=1 ann=an*3.14/180;%(rad)法面压力角 at=atan(tan(ann)/cos(bbb));%（rad)端面压力角 att=acos((d1+d2)*cos(at)/(2*aa));%(rad)实际端面压力角 rb1=d1*cos(at)/2;%小齿轮基圆半径 rb2=d2*cos(at)/2;%大齿轮基圆半径 yt=(z1+z2)/2*((cos(at)/cos(att))-1);%中心距变动系数 ytt=x1+x2-yt;%齿高变动系数 ra1=d1/2+(ha+x1-ytt)*mn/cos(bbb);%小齿顶圆半径 ra2=d2/2+(ha+x2-ytt)*mn/cos(bbb);%大齿顶圆半径 aa1=acos(rb1/ra1);%重合度中参数 aa2=acos(rb2/ra2);%重合度中参数 ea1=z1*(tan(aa1)-tan(att));%重合度中参数 ea2=z2*(tan(aa2)-tan(att));%重合度中参数 ea3=ea1+ea2;%重合度中参数 ea=ea3/(2*3.14);%重合度 else ea=1;%重合度 end %斜齿轮端面向重合度 ZE=189.8;%%%（根号下MPa）弹性影响系数 ZH=2.46;%%%区域系数 u=z1/z2;%小齿轮/大齿轮的传动比 %选择结构，给定内/外啮合判定公式 w=1;%外啮合判定条件，否则为0 if w>=1 uu=(u+1)/u; else uu=(u-1)/u; end %选择结构，给定内/外啮合判定公式 sH=sqrt((k*Ft*uu)/(b*d1*ea*u))*ZH*ZE;%（MPa）疲劳疲劳强度 end