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)...