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

zc1987530 · 发表于 2009-4-24 14:32

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

zhanglw · 发表于 2010-8-11 08:49

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

seekyong · 发表于 2011-1-13 21:13

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

yejet · 发表于 2011-1-14 09:20

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

jsczjp2000 · 发表于 2011-6-22 16:11

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 · 发表于 2011-11-24 14:38

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

mwglikai · 发表于 2012-5-11 15:57

回复 5 # jsczjp2000 的帖子

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

S11021023 · 发表于 2012-7-22 21:27

回复 7 # mwglikai 的帖子

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

紫苏9296 · 发表于 2014-11-17 15:11

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

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

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

lihaotian · 发表于 2015-11-16 16:17

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

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

wjt1713573225 · 发表于 2016-3-10 15:01

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 · 发表于 2016-3-21 14:27

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

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

账号		自动登录	找回密码
密码			我要加入

[转子动力学] 谁有差分法求解雷诺方程的MATLAB程序？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

同楼上

点评