gusucode.com > 润滑理论的求解程序包括弹流和刚流润滑不成功源码程序 > thickness.m

    function hmin=thickness(m,n,e,P)  %本函数的功能是求最小油膜厚度
delta_m=2*pi/m;     %沿圆周方向均匀划分m格，圆周方向的无量纲等分数
delta_n=2/n;        %沿轴向均匀划分n格，轴向的无量纲等分数
for i=1:m+1
for j=1:n+1
   H(i,j)=1+e*cos((i-1)*delta_m);
end
end      %形成油膜无量纲厚度矩阵

hmin=H(1,n/2+1);    %求最小油膜厚度
for i=2:m
    if(H(i,n/2+1)<hmin)
        hmin=H(i,n/2+1);
    end
end