[求助] 新手求教UDF问题~

3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

3楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{
real time= RP_Get_Real("flow-time";
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
if (Q_tot<= Heat[i] )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_c_loop(c,t)
}

楼: Originally posted by xiegangmai at 2012-01-04 22:46:29:
请使用插入代码的功能

楼: Originally posted by wuming524 at 2012-01-04 22:06:19:
楼主太抬举俺了，出现这种问题你需要检查你的语句，大多出现调用无效语句出现的错误。本人也不是太懂。待高人解答吧~好运。

楼: Originally posted by kylafree at 2012-01-04 22:06:24:
real Q_tot=0;
real Heat[3]={278.5, 280.7, 278.5};

if (Q_tot<= Heat )
这里是一个数组和一个数在比大小吗？

楼: Originally posted by gyctju at 2012-01-08 21:11:19:
可能你想法错了。边界条件部分不应该跟累加传热量的循环放在一起，这样的实际效果跟你想的是不一样的。
还有边界条件设定时应该用begin-f-loop
建议修改成如下形式试试
#include "udf.h"
real Q_t ...

楼: Originally posted by gyctju at 2012-01-09 10:45:24:
我上边回复的那个修改大括号那块出错了，把begine_f_loop前面的大括号挪到最后去

16楼: Originally posted by gyctju at 2012-01-09 19:16:54:
again，边界条件用begin_f_loop
估计你都没看我修改后的语句。

#include "udf.h"
#define cp 800.
#define density 2180.
FILE *fp;

DEFINE_PROFILE(unsteady_temp,t,index)
{
float tempn,templ,volume,Q,Q_tot,time,a,b,tempad;
int i;
float Heat[5]={2780000.,        280.7,        278.5,        280.5,        285.};  
Domain *d;
Thread *ct,*t0;
cell_t c,c0;
face_t f;
time = RP_Get_Real("flow-time");
a = time/3600;
b = a/24;
d = Get_Domain(1);
i=(int)(((int)a)/24);

if (b==(int)b)   
Q_tot=0;
else
{
fp=fopen("Q_tot.txt","r");
fscanf(fp,"%f",Q_tot);
fclose(fp);
}
thread_loop_c(ct,d)  
{
        begin_c_loop(c,ct) /*Loop over all cells循环所有THREAD上所有CELL */
       {
volume=C_VOLUME(c,ct); /* get cell volume */
templ=C_T_M1(c,ct); /*Get cell tempertuare of previous step有问题！！！！！！！！*/
tempn=C_T(c,ct);   /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
       }
     end_c_loop(c,ct)
}
fp=fopen("Q_tot.txt","w");
  fprintf(fp,"%f",Q_tot);    /*写入累加后Q_tot值 */
fclose(fp);
if (Q_tot<= Heat[i])
  begin_f_loop(f,t)
{   
    F_PROFILE(f,t,index)=333;  /*BC:Temp=333K*/
}
end_f_loop(f,t)
else
begin_f_loop(f,t)
{   
c0 = F_C0(f,t);
t0 = THREAD_T0(t);
tempad=C_T(c0,t0); /*temperature of adjacent cell*/
F_PROFILE(f,t,index)=tempad;
}
end_f_loop(f,t)
}