DEFINE_ON_DEMAND(print_force)
{
   real A[ND_ND];
   face_t f;
   real dfx=0,dfy=0;
   int zone_ID=3;
   Thread*thread_name=Lookup_Thread(domain,zone_ID);
   begin_f_loop(f,thread)
   {
         F_AREA(A,f,thread)
         dfx+=F_P(f,thread)*A[1];
         dfy+=F_P(f,thread)*A[0];
     }
     end_f_loop(f,thread)
     printf("x-force=%f y-force=%f,dfx,dfy)

2楼: Originally posted by gxw13131 at 2015-05-11 12:58:09
不明白你的意思，但你的udf里那个thread_name得改成thread,其实就是一个遍历用的指针。单元面积分量就是通过F_AREA()这个宏得到的，放在A[]数组里面，A为x方向，A为y方向。三维的话以此类推。
...

3楼: Originally posted by liubojia at 2015-05-11 14:53:15
现在的问题就是如何在F_AREA()中定义A[]这个你会不帮我下出来行不
然后我把金币给你  
F_AREA(A,f,thread)
dfx+=F_P(f,thread)*F_AREA(A,f,thread);
  dfy+=F_P(f,thread)*F_AREA(A,f,thread)...

4楼: Originally posted by gxw13131 at 2015-05-11 19:05:23
A[]你不是在前面已经定义了吗？F_AREA(A,f,thread)这个语句实际上已经把单元各个面积分量输出到A[]里去了。
你这个udf不是想计算圆柱的升力和阻力吗？完全没必要用udf。在monitors里面，残差下面那个create，就可 ...

4楼: Originally posted by gxw13131 at 2015-05-11 19:05:23
A[]你不是在前面已经定义了吗？F_AREA(A,f,thread)这个语句实际上已经把单元各个面积分量输出到A[]里去了。
你这个udf不是想计算圆柱的升力和阻力吗？完全没必要用udf。在monitors里面，残差下面那个create，就可 ...