【机械设计】基于matlab二级齿轮减速器斜齿齿轮传动（齿面接触疲劳强度判定）【含Matlab源码 8054期】.zip

% 机械设计----斜齿齿轮传动 clear;clc P=input('请输入P='); n_1=input('请输入n_1='); z_1=input('请输入z_1='); i=input('请输入i='); z_2=ceil(i.*z_1); disp('z_2的值为：');disp(z_2); % 按齿根弯曲疲劳强度设计 T_1=9.55.*10.^6.*(P./n_1); disp('T_1的值为：');disp(T_1); K_t=input('请输入K_t='); beta=input('请输入beta='); fai_d=input('请输入fai_d='); Sigma_FE1=input('请输入Sigma_FE1='); Sigma_FE2=input('请输入Sigma_FE2='); S_Fmin=input('请输入S_Fmin='); j=input('请输入j=');L_h=input('请输入L_h='); N_1=60.*n_1.*j.*L_h; disp('N_1的值为：'); disp(N_1); N_2=N_1./i; disp('N_2的值为：'); disp(N_2); K_FN1=input('请输入K_FN1='); K_FN2=input('请输入K_FN2='); Sigma_F1_xy=K_FN1.*Sigma_FE1./S_Fmin; disp('Sigma_F1_xy的值为：'); disp(Sigma_F1_xy); Sigma_F2_xy=K_FN2.*Sigma_FE2./S_Fmin; disp('Sigma_F2_xy的值为：'); disp(Sigma_F2_xy); % 计算当量齿数 z_v1=z_1./((cos(beta)).^3); z_v2=z_2./((cos(beta)).^3); disp('z_v1的值为：'); disp(z_v1); disp('z_v2的值为：'); disp(z_v2); Y_Fa1=input('请输入Y_Fa1='); Y_Sa1=input('请输入Y_Sa1='); Y_Fa2=input('请输入Y_Fa2='); Y_Sa2=input('请输入Y_Sa2='); YYSigma1=Y_Fa1.*Y_Sa1./Sigma_F1_xy; YYSigma2=Y_Fa2.*Y_Sa2./Sigma_F2_xy; disp('YYSigma1的值为：');disp(YYSigma1); disp('YYSigma2值为：');disp(YYSigma2); Epsilon_a=input('请输入Epsilon_a='); Epsilon_beta=0.318.*fai_d.*z_1.*tan(beta); if Epsilon_beta > 1; Epsilon_beta = 1; else Epsilon_beta = Epsilon_beta; end disp('Epsilon_beta的值为：');disp(Epsilon_beta); a_n=input('请输入a_n='); a_t=atan(tan(a_n)./cos(beta)); disp('a_t的值为：');disp(a_t); beta_b=atan(tan(beta).*cos(a_t)); disp('beta_b的值为：');disp(beta_b); Y_Epsilon=0.25+(0.75.*(cos(beta_b)).^2)./Epsilon_a; Y_beta=1-Epsilon_beta.*beta./(120*pi/180); disp('Y_Epsilon的值为：');disp(Y_Epsilon); disp('Y_beta的值为：');disp(Y_beta); m_n=((2.*K_t.*T_1.*(cos(beta)).^2./(fai_d.*(z_1).^2)).*max(YYSigma1,YYSigma2).*Y_Epsilon.*Y_beta).^(1/3); disp('m_n的值为：');disp(m_n); m_n_bz=ceil(m_n); disp('m_n_bz的值为：');disp(m_n_bz); d_1=m_n_bz.*z_1./cos(beta); disp('d_1的值为：');disp(d_1); v=pi.*d_1.*n_1./60000; disp('v的值为：');disp(v); b=fai_d.*d_1; disp('b的值为：');disp(b); % 计算载荷系数K K_A=input('请输入K_A='); K_v=input('请输入K_v='); K_Ha=input('请输入K_Ha='); K_Fa=input('请输入K_Fa='); K_H_beta=1.05+0.26.*(1+0.6.*(fai_d).^2).*(fai_d).^2+0.16.*10.^(-3).*b;% 硬齿面齿轮中小齿轮相对支撑非对称布置，6级精度 disp('K_H_beta的值为：');disp(K_H_beta); K_F_beta=K_H_beta; disp('K_F_beta的值为：');disp(K_F_beta); K= K_A.*K_v.*K_Fa.*K_F_beta; disp('K的值为：');disp(K); % 几何尺寸计算 a=(z_1+z_2).*m_n_bz./(2.*cos(beta)); a_bz=ceil(a); disp('a_bz的值为：');disp(a_bz); beta_xz=acos((z_1+z_2).*m_n_bz./(2.*a_bz)); disp('beta_xz的值为：');disp(beta_xz); d_1_xz=m_n_bz.*z_1./cos(beta_xz); d_2_xz=m_n_bz.*z_2./cos(beta_xz); disp('d_1_xz的值为：');disp(d_1_xz); disp('d_2_xz的值为：');disp(d_2_xz); b_xz=fai_d.*d_1_xz; disp('b_xz的值为：');disp(b_xz); b_xz_bz_2=ceil(b_xz); disp('b_xz_bz_2的值为：');disp(b_xz_bz_2); b_xz_bz_1=b_xz_bz_2+5; disp('b_xz_bz_1的值为：');disp(b_xz_bz_1); % 校核齿面接触疲劳强度 Sigma_Hlim1=input('请输入Sigma_Hlim1='); Sigma_Hlim2=input('请输入Sigma_Hlim2='); S_Hlim=input('请输入S_Hlim='); K_HN1=input('请输入K_HN1='); K_HN2=input('请输入K_HN2='); Sigma_H1_xy=K_HN1.*Sigma_Hlim1./S_Hlim; Sigma_H2_xy=K_HN2.*Sigma_Hlim2./S_Hlim; Sigma_H_xy=(Sigma_H1_xy+Sigma_H2_xy)./2; disp('Sigma_H1_xy的值为：');disp(Sigma_H1_xy); disp('Sigma_H2_xy的值为：');disp(Sigma_H2_xy); disp('Sigma_H_xy的值为：');disp(Sigma_H_xy); Z_E=input('请输入Z_E='); Z_H=input('请输入Z_H='); Epsilon_a_jc=input('请输入Epsilon_a_jc='); Epsilon_beta_jc=input('请输入Epsilon_beta_jc='); Z_Epsilon=((4-Epsilon_a_jc).*(1-Epsilon_beta_jc)./3+Epsilon_beta_jc./Epsilon_a_jc).^(1/2); disp('Z_Epsilon的值为：');disp(Z_Epsilon); Z_beta=(cos(beta_xz)).^(1/2); disp('Z_beta的值为：');disp(Z_beta); % 外啮合 Sigma_H=Z_E.*Z_H.*Z_Epsilon.*Z_beta.*((2.*K.*T_1)./(fai_d.*(d_1_xz).^3).*((i+1)./i)).^(1/2); disp('Sigma_H的值为：');disp(Sigma_H); if Sigma_H <= Sigma_H1_xy; Sigma_H <= Sigma_H2_xy; disp('齿面接触疲劳强度合格'); else disp('齿面接触疲劳强度不合格'); end