[求助]请高手帮帮忙！

stan · 发表于 2005-8-4 01:08:45

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

您需要登录才可以下载或查看，没有账号？注册

x

[这个贴子最后由stan在 2005/09/13 10:02pm 第 2 次编辑]

哪位高手可以帮我把下面这个函数编个udf：

谢谢！

xiaofish · 发表于 2005-8-5 13:48:40

以前的帖子上copy的
改一改就可以了
#include "udf.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t = CURRENT_TIME;
real t3,t2,t1 = 1.6 e-05;
real k,b=0.006;
  k=(4 * b) / t1;
begin_f_loop(f, thread)
{
t2 = t / t1;
t3 = t - floor(t2) * t1;
if  (t3<=t1/4)
   F_PROFILE(f, thread, position) =0.021+k*t3;
else if  (t3>t1/4 && t3<=(3*t1)/4)
   F_PROFILE(f, thread, position) =0.021-k*t3+2*b;

else
F_PROFILE(f, thread, position) =0.021+k*t3-4*b;
}
end_f_loop(f, thread)
}

stan · 发表于 2005-8-5 15:27:47

谢谢！
但这个与我的函数差别还是很大的，我还是不会。

zhubin · 发表于 2005-8-7 14:03:26

我来写写看,高手看对不对,刚学了点udf,呵呵,如果有错误还请多多指导啊!
#include "udf.h"
#include "math.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t = CURRENT_TIME;
begin_f_loop(f, thread)
  {
if  (t>0&&t<0.24)
F_PROFILE(f, thread, position) =240*exp(-7.557t)*sin(13.09t);
else if  (t>0.24&&t<=0.2743)
F_PROFILE(f, thread, position) =-(34.3)*exp(-7.557)*sin(13.09*(t-0.24))

else
F_PROFILE(f, thread, position) =0;
  }
end_f_loop(f, thread)
}

stan · 发表于 2005-8-8 12:19:50

[这个贴子最后由stan在 2005/08/08 12:28pm 第 1 次编辑]

下面引用由zhubin在 2005/08/07 02:03pm 发表的内容：
我来写写看,高手看对不对,刚学了点udf,呵呵,如果有错误还请多多指导啊!
#include "udf.h"
#include "math.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
...

调试时出现
Error: D:\01\fluent\v.c: line 14: parse error.
Error: D:\01\fluent\v.c: line 20: parse error.
另外请问在哪里设置t<=0.8?
谢谢！

bluesharpe · 发表于 2005-12-2 10:54:11

我表示怀疑，但是不一定怀疑的就对。DEFINE_PROFILEMIN

bluesharpe · 发表于 2005-12-2 11:03:47

DEFINE_PROFILE命令要在迭代循环前就算好了，那我估计你的这个程序就会一直以第一种情况作为速度初值计算。在进入迭代后就不会去根据你的流体时间来判断你的速度初值的。

Mclaren · 发表于 2005-12-2 11:11:57

LBSALE[1]LBSALE[这个贴子最后由Mclaren在 2005/12/02 11:14am 第 1 次编辑]

我来试一下
#include "udf.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t = RP_Get_Real("flow-time");
begin_f_loop(f, thread)
{
if  (t>0&&t<0.24)
F_PROFILE(f,thread,position) =240*exp(-7.557t)*sin(13.09t);
  else if  (t>0.24&&t<=0.2743)
F_PROFILE(f, thread, position) =-(34.3)*exp(-7.557)*sin(13.09*(t-0.24));
  else
  F_PROFILE(f, thread, position) =0;
}
end_f_loop(f, thread)
}

Mclaren · 发表于 2005-12-2 11:17:15

我在上面的帖子里面把程序贴出来了
想看的自己点
一个币，应该的

cocowei · 发表于 2005-12-23 23:58:40

注意14行，应该为13.09*t

cocowei · 发表于 2005-12-24 00:00:12

可以用DEFINE_ADJUST来试试

sanmango · 发表于 2006-5-18 11:08:36

下面引用由stan在 2005/08/08 00:19pm 发表的内容：
调试时出现
Error line 14 D line 20: parse error.
另外请问在哪里设置t<=0.8?
谢谢！

renwoxing · 发表于 2007-8-20 10:58:35

下面引用由zhubin在 2005/08/07 02:03pm 发表的内容：
我来写写看,高手看对不对,刚学了点udf,呵呵,如果有错误还请多多指导啊!
#include "udf.h"
#include "math.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
...

你的十四行,少个*号
在十九行,少个;号
编译误把１９，２０当成一行．故有错误．

		自动登录	找回密码
密码			注册

[求助]请高手帮帮忙！

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