matlab角接触球轴承和圆锥滚子轴承轴承寿命校核

**

matlab角接触球轴承和圆锥滚子轴承轴承寿命校核

**
自己修改轴承所受的力，以及受力方向等

Fr1=1000; %轴承1上的径向力，单位N%%%%%%%%%%%%%%%%%%%%%%%%%%
Fr2=2000; %轴承2上的径向力,单位N%%%%%%%%%%%%%%%%%%%%%%%%%%
Fa=1000;  %轴上的轴向力，单位N，规定向右为正%%%%%%%%%%%%%%%%%%%%%%%%%%
C=33200;  %基本额定动载荷，单位N%%%%%%%%%%%%%%%%%%%%%%%%%%
n=500;    %转速，单位r/min
type=7;   %轴承类型 3为圆锥滚子轴承 7为角接触球轴承 1为其他轴承%%%%%%%%%%%%%%%%%%%%%%%%%%
D=1;      %1为轴承正装，2为轴承反装%%%%%%%%%%%%%%%%%%%%%%%%%%
coe=(10^6)/(60*n);%寿命公式系数
T=10000;  %轴承要求寿命
if type==3 %判断轴承类型%%%%%%%%%%%%%%%%%%%%%%%%%%
    fprintf('轴承类型：圆锥滚子轴承。\n');
    e=0.80;%判断系数%%%%%%%%%%%%%%%%%%%%%%%%%%
    ep=10/3;%Epsilon
    if D==1 %判断正反装，，，，此为正装
         fprintf('轴承安装方式：正装。\n');
         %1求派生轴向力
         Y=e;
         Fd1=Fr1/(2*Y); %轴承1上的派生轴向力
         Fd2=-Fr2/(2*Y); %轴承2上的派生轴向力
         %2判断哪个轴承被压紧并求出每个轴承上所受的轴向力
         if Fd1+Fd2+Fa<=0 %轴承1被压紧
             fprintf('轴承1被压紧。\n');
             Fa1=abs(Fa+Fd2); %轴承1上的轴向力
             Fa2=Fd2;         %轴承2上的轴向力
         else  %轴承2被压紧
              fprintf('轴承2被压紧。\n');
              Fa1=Fd1;        %轴承1上的轴向力
              Fa2=abs(Fa+Fd1);%轴承2上的轴向力
         end
         %3轴承当量动载荷
         fd=1; %载荷系数%%%%%%%%%%%%%%%%%%%%%%%%%
         %轴承1的当量动载荷
         if Fa1/Fr1<=e
             X=1; %径向动载荷系数%%%%%%%%%%%%%%%%%%%%%%%%%%
             P1=X*Fr1;
         else
             X=0.41; %径向动载荷系数%%%%%%%%%%%%%%%%%%%%%%%%%%
             Y=0.87; %轴向动载荷系数%%%%%%%%%%%%%%%%%%%%%%%%%%
             P1=(X*Fr1+Y