润滑理论的matlab求解程序包括弹流和刚流润滑.rar

【课题研究】润滑理论的matlab求解程序,包括弹流和刚流润滑.rar

function P=dispressure（emnratio）
%This function is to calculate the pressure distribution of journal bearing without taking the deformation into account
delta_m=2*pi/m;              %沿圆周方向均匀划分m格
delta_n=2/n;                  %沿轴向均匀划分n格
for j=1:n+1                   %对各节点赋压力初值0
for i=1:m+1
 P（ij）=0;
 Q（ij）=0;
end
end
S=0;
T=0;
for j=2:n                     %第一次计算各内部节点的压力值
for i=2:m
A=（1+e*cos（（i+1/2-1）*delta_m））^3;
B=（1+e*cos（（i-1/2-1）*delta_m））^3;
C=（（1/ratio）*delta_m/delta_n）^2*（1+e*cos（（i-1）*delta_m））^3;
D=C;
E=A+B+C+D;
F=delta_m*e*（cos（（i+1/2-1）*delta_m）-cos（（i-1/2-1）*delta_m））;
P（ij）=（A*P（i+1j）+B*P（i-1j）+C*P（ij+1）+D*P（ij-1）-F）/E;
if（P（ij）<=0）                     %引入雷诺边界条件
   for k=i:m
     P（kj）=0;
     S=S+abs（Q（kj））;
     T=T+abs（P（ij））;
     Q（kj）=P（kj）;
end
break;
end
S=S+abs（P（ij）-Q（ij））;
T=T+abs（P（ij））;
Q（ij）=P（ij）;
end	
end
while（S/T>0.001）         %循环计算各内部节点的压力值
S=0;
T=0;
for j=2:n
for i=2:m
  A=（1+e*cos（（i+1/2-1）*delta_m））^3;
  B=（1+e*cos（（i-1/2-1）*delta_m））^3;
  C=（（1/ratio）*delta_m/delta_n）^2*（1+e*cos（（i-1）*delta_m））^3;
  D=C;
  E=A+B+C+D;
  F=delta_m*e*（cos（（i+1/2-1）*delta_m）-cos（（i-1/2-1）*delta_m））;
  P（ij）=（A*P（i+1j）+B*P（i-1j）+C*P（ij+1）+D*P（ij-1）-F）/E;
  if（P（ij）<=0）
     for k=i:m
         P（kj）=0;
         S=S+abs（Q（kj））;
         T=T+abs（P（kj））;
         Q（kj）=P（kj）;
     end
      break;
end
  S=S+abs（P（ij）-Q（ij））;
  T=T+abs（P（ij））;
  Q（ij）=P（ij）;
end
end
end


%P;     %此部分是当本函数单独使用时，用于绘制离散压力分布
%for i=1:m+1
 % x（i）=（i-1）*delta_m;
 %end
%for j=1:n+1
 %y（j）=（j-1）*delta_n;
 %end
%figure
%[xy]=meshgrid（xy）
%mesh（P）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1758  2007-01-09 16:57  NewMatlab\dispressure.m

     文件        257  2007-01-09 16:57  NewMatlab\filmrupture.m

     文件        735  2007-01-09 16:57  NewMatlab\flow_quantity.m

     文件        816  2007-01-09 16:57  NewMatlab\friction.m

     文件        179  2007-01-09 16:57  NewMatlab\pressness.m

     文件        407  2007-01-09 16:56  NewMatlab\thickness.m

     文件        381  2007-01-09 16:56  NewMatlab\totalpressure.m

     文件       2249  2007-01-09 16:56  NewMatlab\Untitled55.m

     目录          0  2007-02-04 14:12  NewMatlab

----------- ---------  ---------- -----  ----

                 7000                    10