谁有差分法求解雷诺方程的MATLAB程序？

zc1987530

谁有差分法求解雷诺方程的MATLAB程序？很需要，能否发给我一份，zhangcheng1987530@foxmail.com
:handshake 十分感谢！

zhanglw

帮忙给我也发一份,急需要，感激不尽     zhanglw1987@163.com

seekyong

有maple的，但是结果和ISO标准有一定偏差。坛子里也没人研究解释

yejet

国内这个东西不是没人做，而且做得人还是不少
不过国内一般人是不会把这个东西放出来的
可以到国外一些学术网站找找

jsczjp2000

function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)
Fx=0;
Fy=0;
epsilon=0.6;
L=0.03;
d=0.06;
lambda=L/d;
m=40;
n=30;
phi1=0;
phi2=2*pi;
delta_phi=(phi2-phi1)/m;
delta_lambda=2/n;
k=1;
P0=zeros(m+1,n+1);
A=zeros(m+1,n+1);
B=zeros(m+1,n+1);
C=zeros(m+1,n+1);
D=zeros(m+1,n+1);
E=zeros(m+1,n+1);
F=zeros(m+1,n+1);
H=zeros(m+1,n+1);
for i=1:1:m+1
%      theta(i)=(i-1)*delta_phi;
        for j=1:1:n+1
            H(i,j)=1+epsilon*cos((i-1)*delta_phi);
        end
end
S=0;
T=0;
ERR=1e-3;
GAP=1;
while GAP>ERR
    k=k+1
    for i=1:1:m+1
        for j=1:1:n+1
            if (i==1)|(j==1)|(i==m+1)|(j==n+1)
                P(i,j)=0;
            else
                A(i,j)=(1+epsilon*cos((i+1/2-1)*delta_phi))^3;
                B(i,j)=(1+epsilon*cos((i-1/2-1)*delta_phi))^3;
                C(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;
                D(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;
                E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j);
                F(i,j)=6*delta_phi*((1+epsilon*cos((i+1/2-1)*delta_phi))-(1+epsilon*cos((i-1/2-1)*delta_phi)))+12*(delta_phi)^2*(X_dot*cos(phi1+(i-1)*delta_phi)+Y_dot*sin(phi1+(i-1)*delta_phi));
                P(i,j)=(A(i,j)*P0(i+1,j)+B(i,j)*P0(i-1,j)+C(i,j)*P0(i,j+1)+D(i,j)*P0(i,j-1)-F(i,j))/E(i,j);
                if P(i,j)<0
                    P(i,j)=0;
                else
                end
            end
        end
    end
%     for i=2:1:m
%         for j=2:1:n
            S=sum(sum(abs(P-P0)));
            T=sum(sum(abs(P)));
%         end
%     end
    GAP=S/T;
    P0=P;
end

for i=1:1:m+1
    for j=1:1:n+1
        Fx=Fx+P(i,j)*cos(phi1+(i-1)*delta_phi)*delta_phi*delta_lambda;
        Fy=Fy+P(i,j)*sin(phi1+(i-1)*delta_phi)*delta_phi*delta_lambda;
    end
end

0810064

请问楼上给出的是求解雷诺方程的MATLAB程序吗？

mwglikai

回复 5 # jsczjp2000 的帖子

下了,不错.不知能否发个动态雷诺方程的matlab程序给我.用差分法求解时变量t 不会处理.谢谢!mwglikai@163.com

S11021023

回复 7 # mwglikai 的帖子

师姐呀！不知你是哪位老师的爱徒

紫苏9296

mwglikai 发表于 2012-5-11 15:57
回复 5 # jsczjp2000 的帖子

下了,不错.不知能否发个动态雷诺方程的matlab程序给我.用差分法求解时变量t  ...

你好 我现在也用有限差分法解动态雷诺方程，里面的时间项不知道怎么写，您能发过来个程序我学习下吗？？？？不胜感谢！！！！！！！！

lihaotian

jsczjp2000 发表于 2011-6-22 16:11
function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)
Fx=0;
Fy=0;

请问下，您这个是求解稳态油膜力的吗？

wjt1713573225

jsczjp2000 发表于 2011-6-22 16:11
function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)
Fx=0;
Fy=0;

请问这个程序中的x,y ,x_dot y-dot 什么意思

NASA

wjt1713573225 发表于 2016-3-10 15:01
请问这个程序中的x,y ,x_dot y-dot 什么意思

x,y应该表示的是位移
x_dot y-dot表示的应该是速度