声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4623|回复: 11

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

[复制链接]
发表于 2009-4-24 14:32 | 显示全部楼层 |阅读模式

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

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

x
谁有差分法求解雷诺方程的MATLAB程序?很需要,能否发给我一份,zhangcheng1987530@foxmail.com
:handshake 十分感谢!
回复
分享到:

使用道具 举报

发表于 2010-8-11 08:49 | 显示全部楼层

同楼上

帮忙给我也发一份,急需要,感激不尽     zhanglw1987@163.com
发表于 2011-1-13 21:13 | 显示全部楼层
有maple的,但是结果和ISO标准有一定偏差。坛子里也没人研究解释
发表于 2011-1-14 09:20 | 显示全部楼层
国内这个东西不是没人做,而且做得人还是不少
不过国内一般人是不会把这个东西放出来的
可以到国外一些学术网站找找
发表于 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
发表于 2011-11-24 14:38 | 显示全部楼层
请问楼上给出的是求解雷诺方程的MATLAB程序吗?
发表于 2012-5-11 15:57 | 显示全部楼层
回复 5 # jsczjp2000 的帖子

下了,不错.不知能否发个动态雷诺方程的matlab程序给我.用差分法求解时变量t 不会处理.谢谢!mwglikai@163.com
发表于 2012-7-22 21:27 | 显示全部楼层
回复 7 # mwglikai 的帖子

师姐呀!不知你是哪位老师的爱徒
发表于 2014-11-17 15:11 | 显示全部楼层
mwglikai 发表于 2012-5-11 15:57
回复 5 # jsczjp2000 的帖子

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

你好 我现在也用有限差分法解动态雷诺方程,里面的时间项不知道怎么写,您能发过来个程序我学习下吗????不胜感谢!!!!!!!!
发表于 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;

请问下,您这个是求解稳态油膜力的吗?
发表于 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 什么意思

点评

x,y应该表示的是位移 x_dot y-dot表示的应该是速度  详情 回复 发表于 2016-3-21 14:27
发表于 2016-3-21 14:27 | 显示全部楼层
wjt1713573225 发表于 2016-3-10 15:01
请问这个程序中的x,y ,x_dot y-dot 什么意思

x,y应该表示的是位移
x_dot y-dot表示的应该是速度
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-3-29 09:44 , Processed in 0.088482 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表