24小时热门版块排行榜    

查看: 703  |  回复: 2

bwj0319

新虫 (小有名气)

[交流] fluen中的UDF程序问题 已有1人参与

模拟往复式压缩机出口管道需要添加一个脉动的压力,脉动形式如图所示,表达式也在下面,周期是0.3s,针对这个问题,由于刚上手fluent不久,对UDF功能也不是很熟,外加C语言也很菜,于是仿照别人的程序编了一段,大家看能实现图示和表达式所需要的脉动压力吗?我的fluent开始算了,看监测点的曲线 好像开始就有脉动了,估计有问题啊 想听听大牛们的意见啊
程序:

/* unsteady.c */
/* UDF for specifying a transient pressure profile boundary condition */
#include "udf.h"
DEFINE_PROFILE(unsteady_pressure,thread,position)
{
face_t f;
begin_f_loop(f,thread)
{
real t=RP_Get_Real("flow-time";

int k;
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)
F_PROFILE(f,thread,position)=310000000;
else
F_PROFILE(f,thread,position)=310000000*(1+0.06*sin(2*(20.93*t-4.71)));
}
}
end_f_loop(f,thread)
}

模拟过程中设置监测点,结果监测点曲线如下附图,监测数据:"Convergence history of Static Pressure on point-1 (in SI units)"
"Flow Time" "Vertex Average Static Pressure"
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0020000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0030000001 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0040000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.099111e+08
0.0050000002 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991114e+08
0.0060000003 3.0991117e+08
0.0060000003 3.0991117e+08
0.0060000003 3.099112e+08
0.0060000003 3.099112e+08
0.0060000003 3.0991123e+08
0.0060000003 3.099113e+08
0.0060000003 3.0991142e+08
0.0060000003 3.0991155e+08
0.0060000003 3.0991168e+08
0.0060000003 3.0991187e+08
0.0070000003 3.0991187e+08
0.0070000003 3.0991216e+08
0.0070000003 3.0991248e+08
0.0070000003 3.0991293e+08
0.0070000003 3.099135e+08
0.0070000003 3.0991421e+08
0.0070000003 3.0991498e+08
0.0070000003 3.0991587e+08
0.0070000003 3.099169e+08
0.0070000003 3.0991805e+08
0.0070000003 3.0991936e+08
0.0070000003 3.099208e+08
0.0070000003 3.0992243e+08
0.0070000003 3.0992416e+08
0.0070000003 3.0992598e+08
0.0070000003 3.0992797e+08
0.0070000003 3.0993011e+08
0.0070000003 3.0993235e+08
0.0070000003 3.0993466e+08
0.0070000003 3.0993709e+08
0.0070000003 3.0993962e+08
0.0080000004 3.0993962e+08
0.0080000004 3.0994237e+08
0.0080000004 3.0994541e+08
0.0080000004 3.099487e+08
0.0080000004 3.0995213e+08
0.0080000004 3.0995571e+08
0.0080000004 3.0995946e+08
0.0080000004 3.0996333e+08
0.0080000004 3.0996723e+08
0.0080000004 3.099712e+08
0.0080000004 3.099752e+08
0.0080000004 3.0997923e+08
0.0080000004 3.0998323e+08
0.0080000004 3.099873e+08
0.0080000004 3.099913e+08
0.0080000004 3.099953e+08
0.0080000004 3.0999933e+08
0.0080000004 3.1000326e+08
0.0080000004 3.1000717e+08
0.0080000004 3.1001104e+08
0.0080000004 3.1001491e+08
0.0090000004 3.1001491e+08
0.0090000004 3.1001882e+08
0.0090000004 3.1002288e+08
0.0090000004 3.1002694e+08
0.0090000004 3.1003114e+08
0.0090000004 3.1003533e+08
0.0090000004 3.1003958e+08
0.0090000004 3.1004384e+08
0.0090000004 3.1004813e+08
0.0090000004 3.1005238e+08
0.0090000004 3.1005664e+08
0.0090000004 3.1006086e+08
0.0090000004 3.1006506e+08
0.0090000004 3.1006918e+08
0.0090000004 3.1007328e+08
0.0090000004 3.1007731e+08
0.0090000004 3.1008128e+08
0.0090000004 3.1008525e+08
0.0090000004 3.1008909e+08
0.0090000004 3.100929e+08
0.0090000004 3.1009664e+08
0.01 3.1009664e+08
0.01 3.1010042e+08
0.01 3.1010413e+08
0.01 3.1010784e+08
0.01 3.1011146e+08
0.01 3.101151e+08
0.01 3.1011869e+08
0.01 3.1012224e+08
0.01 3.1012573e+08
0.01 3.1012918e+08
0.01 3.1013264e+08
0.01 3.1013606e+08
0.01 3.1013946e+08
0.01 3.1014278e+08
0.01 3.1014614e+08
0.01 3.1014944e+08
0.01 3.101527e+08
0.01 3.1015594e+08
0.01 3.1015917e+08
0.01 3.1016234e+08
0.01 3.101655e+08
0.011000001 3.101655e+08
0.011000001 3.1016867e+08
。。。。。。
------------------------------------------------------
结果应该不对,没能实现需要的压力曲线啊 ,有懂UDF程序的大牛们 指导下俺吧 ,愁死了。。。。


压力曲线



压力表达式

监测点压力曲线



[ Last edited by bwj0319 on 2012-7-29 at 21:46 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2012-07-30 21:57:53
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)

k=0,if里面t就是从0开始了
2楼2012-07-29 19:26:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bwj0319

新虫 (小有名气)

引用回帖:
2楼: Originally posted by kylafree at 2012-07-29 19:26:27
for(k=0;k<=9;k=k+1)
{
if(t>=0.3*k&&t<0.225+0.3*k)

k=0,if里面t就是从0开始了

是啊  我就是要从0开始啊,这里是个周期函数,周期为0.3s,每周期里又分为两段了
3楼2012-07-29 22:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bwj0319 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见