多重网格法求解雷诺方程

张抛砖

哪位大神用多重网格法求解雷诺方程过，求资料，求matlab程序，求赐教，，邮箱 88984718@qq.com

huaijuliu

MATLAB程序执行效率很慢，请参考VENNER & LUBRECHT <MULTILEVEL METHOD IN LUBRICATION>一书，后面附录有c++语言源程序。

scott198510

1. clc;clear all;close all;
   
2. epison=0.4;
   
3. nx=180;
   
4. ny=100;
   
5. deltax=2*180/nx;
   
6. deltay=2/ny;
   
7. d=75;
   
8. L=50;
   
9. lamta=d/L;
   
10. lam=(lamta*deltax*pi/180*deltay)^2;
    
11. P=1.0*ones(nx+1,ny+1);
    
12. P(1,:)=1;
    
13. P(nx+1,:)=1;
    
14. P(:,1)=1;
    
15. P(:,ny+1)=1;
    
16. H=ones(nx+1,ny+1);
    
17. A=ones(nx+1,ny+1);
    
18. B=ones(nx+1,ny+1);
    
19. C=ones(nx+1,ny+1);
    
20. D=ones(nx+1,ny+1);
    
21. E=ones(nx+1,ny+1);
    
22. F=ones(nx+1,ny+1);
    
23. P2=ones(nx+1,ny+1);
    
24. P3=ones(nx+1,ny+1);
    
25. %计算间隙初始值
    
26. for i=1:1:nx+1
    
27. theta=(i-1)*deltax*pi/180;
    
28. for j=1:1:ny+1
    
29. H(i,j)=1+epison*cos(theta);
    
30. end
    
31. end
    
32. H
    
33. figure(1);%初始间隙分布'
    
34. set(gcf,'color','w');
    
35. [x,y]=meshgrid((0:deltay:2),(0:deltax:2*180));
    
36. mesh(x,y,H)
    
37. axis([0 2 0 2*180 0 1.5])
    
38. S=0;T=0;
    
39. wucha=1;
    
40. error=10^(-3);
    
41. count=0;
    
42. w=1.75;
    
43. while(wucha>=error)
    
44. PP=P;
    
45. for i=2:nx
    
46. for j=2:ny
    
47. P2(i,j)=P(i,j);
    
48. A(i,j)=H(i+1,j)^3;
    
49. B(i,j)=H(i-1,j)^3;
    
50. C(i,j)=lam*H(i,j+1)^3;
    
51. D(i,j)=lam*H(i,j-1)^3;
    
52. E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j);
    
53. F(i,j)=3*(deltax*pi/180)*(H(i+1,j)-H(i-1,j));
    
54. P3(i,j)=(A(i,j)*P(i+1,j)+B(i,j)*P(i-1,j)+C(i,j)*P(i,j+1)+D(i,j)*P(i,j-1)-F(i,j))/E(i,j);
    
55. P(i,j)=(1-w)*P2(i,j)+w*P3(i,j);
    
56. if P(i,j)<1
    
57. P(i,j)=1;
    
58. end
    
59. end
    
60. end
    
61. for i=2:nx
    
62. for j=2:ny
    
63. S=S+P(i,j)-PP(i,j);
    
64. T=T+P(i,j);
    
65. end
    
66. end
    
67. wucha=S/T;
    
68. count=count+1
    
69. end
    
70. P
    
71. figure(2);%压力分布图
    
72. [x,y]=meshgrid((0:deltay:2),(0:deltax:2*180));
    
73. mesh(x,y,P)
    
74. axis([0 2 0 2*180 0 3])
    
75. set(gcf,'color','w');

复制代码

多重网格解雷诺方程倒是没玩过，附个简单的采用一般差分形式解雷诺方程的例子，希望能有所启发。

hutshine

这个程序好像有点问题啊，压力分布是从0开始的，不是从1开始的