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