润滑轴承油膜压力分布matlab程序

clc;clear;clf;
D=0.52;
R=D/2;
c=0.00045;
D2=0.5209;
L=1.06;
l=0.003;
u=0.096;
namda=1.2;
omiga=89;
E=100000;
v=0.49;
v0=sqrt(2*v^2/(1-v));

e=0.85;
x=0;
y=-c*e;
vx=0;
vy=0;
pq=1;
p0=0;
m=60;
n=60;

deltL=2*(L/2)/m;
deltsita=2*pi/n;

ERR=1.0e-5;
PK=zeros(n+1,m+1); aa=zeros(n+1,m+1); bb=zeros(n+1,m+1); cc=zeros(n+1,m+1);dd=zeros(n+1,m+1); ee=zeros(n+1,m+1); ff=zeros(n+1,m+1);
k=1;
while k>0
P=PK;
for i=1:n+1
for j=1:m+1
H(i,j)=1+e*(cos((i-1)deltsita))+6l*(1-v0.^{2)*u*omiga*R.}2./(Ec.^3).P(i,j);
end
end
for i=2:n
for j=2:m
aa(i,j)=H(i,j)^{2*(H(i,j)-3/4*(H(i+1,j)-H(i-1,j)))/(deltsita}2);
bb(i,j)=H(i,j)^{2*(H(i,j)+3/4*(H(i+1,j)-H(i-1,j)))/(deltsita}2);
cc(i,j)=-2H(i,j)^{3*(1/(deltsita}2)+1/(deltL^2));
dd(i,j)=H(i,j)^{2*(H(i,j)-3/4*(H(i,j+1)-H(i,j-1)))/(deltL}2);
ee(i,j)=H(i,j)^{2*(H(i,j)+3/4*(H(i,j+1)-H(i,j-1)))/(deltL}2);
%f(j,i)=(H(j,i+1)-H(j,i-1))/(2deltsita)+2*(-VYcos(deltsita(i-1))+VXsin(deltsita(i-1)));
ff(i,j)=(H(i+1,j)-H(i-1,j))/(2*deltsita);
end
end
for i=2:n
for j=2:m
% PK(i,j)=(1-namda)P(i,j)+namda(ff(i,j)-(aa(i,j)*P(i-1,j)+bb(i,j)*P(i+1,j)+dd(i,j)*P(i,j-1)+ee(i,j)*P(i,j+1)))/cc(i,j);
PK(i,j)=(ff(i,j)-(aa(i,j)*P(i-1,j)+bb(i,j)*P(i+1,j)+dd(i,j)*P(i,j-1)+ee(i,j)*P(i,j+1)))/cc(i,j);
if PK(i,j)<0
PK(i,j)=0;
break;
end
end
end
sum1=0;
sum2=0;
for s=1:n+1
for t=1:m+1
sum1=sum1+abs((PK(s,t)-P(s,t)));
sum2=sum2+abs(P(s,t));
end
end
sum=sum1/sum2;
if sum<=ERR
break;
end
k=k+1;
end
VV=6.*u.*omiga.*2.*pi./60.*R.^2./c.2;
P=P.*VV;
sitas=(0:n)deltsita;
Ls=(0:m)(L/m);
[SITA,LL]=meshgrid(sitas,Ls);
mesh(SITA,LL,P’)
Fx=0;
Fy=0;
for i=1:n+1
for j=1:m+1
Fx=Fx+P(i,j)*sin((i-1)deltsita)deltsitadeltL;
Fy=Fy-P(i,j)cos((i-1)deltsita)deltsitadeltL;
end
end
fx=6.u.omigaR.^4./c.2Fx;
fy=6.u.omigaR.^4./c.2Fy;
f=sqrt(fx^2+fy2);
theta=180/piatan(fy/fx);
Mx=0;
My=0;
for i=1:n+1
for j=1:m+1
Mx=Mx-P(i,j)(deltL(m/2+1-j))cos((i-1)deltsita)deltsitadeltL;
My=My+P(i,j)(deltL(j-m/2-1))*sin((i-1)deltsita)deltsitadeltL;
end
end
mx=6.u.omigaR.^5./c.2Mx;
my=6.u.omigaR.^5./c.2My;
m=sqrt(mx^2+my2);
theta=180/piatan(my/mx);