求助！压力出口udf，自己编了一个不对，哪位大侠给看看？

loober0932

我想编一个yz面上水压力随z方向改变却编程y方向的了，不知怎么回事，那位大侠救救我！
#include "udf.h"
DEFINE_PROFILE(outlet_z_pressure, thread, position)
{
  real y[ND_ND];
  real z;
  face_t f;
  begin_f_loop(f,thread)
  {
      F_CENTROID(y,f,thread);
      z= y[1];
  F_PROFILE(f,thread,position)=-9.8*998.2*z;
  }
  end_f_loop(f, thread);
}

[ 本帖最后由 loober0932 于 2009-5-29 16:16 编辑 ]

alicefhy

y[0]、y[1]、y[2]分别是x、y、z坐标

loober0932

尝试过了，没用唉。。。。出口的压力还是没变

yi_chao19791001

看看用x【0】，x【1】，x【2】；表示有什么问题。

fastbb

楼主检查一下gambit模型，是不是model中方向设错了

lxy732801

假设水位是z0，则压力应该为
9.8*998.2*（z0-z)