clear all clc N=1000; T=02piN2pi; r=50;基圆半径为50 h =5; 齿高 z1=6;齿数 z2=z1; d=100;中心距 ri_1=dz1z1+z2;节圆半径 h1=08;小齿齿高 m=6;小齿齿数 z_1=z1m; z_2=z_1; l=h1sinz1mT;极坐标系下的齿廓 u=hsinz1T+ r ; c=hsinz1T+r+h1sin z1 mT;

title('Conjugate Profiles of Gear Teeth'); xlabel('Angle (rad)'); ylabel('Distance from Center (mm)'); legend('Original Gear Tooth', 'Conjugate Gear Tooth');

% 计算齿轮传动比 i=z2/z1 Ri=d/z1; Ro=d/z2; i=Ri/Ro; disp(['Gear ratio is ',num2str(i)]);

% 计算齿轮的模数和压力角 m=d/z1; phi=atan((h1z1)/(dpi)); disp(['Module is ',num2str(m),' mm']); disp(['Pressure angle is ',num2str(phi*180/pi),' degrees']);

% 计算齿面接触应力 F=1000; % 输入的力（N） b=25; % 齿宽（mm） K=1.5; % 应力集中系数 sigma_H=F/(bm); sigma_F=F/(bm*K); disp(['Hertzian contact stress is ',num2str(sigma_H),' MPa']); disp(['Maximum surface stress is ',num2str(sigma_F),' MPa']);

% 计算齿轮的弯曲疲劳寿命 Sut=1000; % 抗拉强度（MPa） Se=0.5Sut; % 弯曲疲劳极限（MPa） Ka=1.25; % 系数 Kb=1; % 系数 Kc=1; % 系数 Kd=1; % 系数 Sa=SeKaKbKcKd; % 弯曲疲劳极限强度（MPa） Nf=Sa^3/(sigma_F^3m*b^2); % 弯曲疲劳寿命（次） disp(['Bending fatigue life is ',num2str(Nf),' cycles']);

% 计算齿轮的接触疲劳寿命 Kf=2.7; % 应力循环裕度系数 qv=10; % 齿面速度（m/s） Z=P/m; % 等效齿数 Kt=1.5; % 动载荷系数 Kv=1; % 应力集中系数 sigma_HP=sqrt(sigma_H^2+3sigma_F^2); % 压力应力幅值（MPa） sigma_FP=Kfsigma_HP/Kt/Kv; % 最大等效应力幅值（MPa） Nf_contact=10^((6.67-log10(sigma_FP))^3); % 接触疲劳寿命（次） disp(['Contact fatigue life is ',num2str(Nf_contact),' cycles'])