请大侠帮我指正一下我编写的UDF存在什么问题？

如烟

我想做一个注射器液体流动仿真，推动力由一个弹簧提供，随着弹簧的放松，推理减小，我想编写一个UDF来设定入口出的压力。编写的UDF如下：
#include "udf.h"
DEFINE_PROFILE(unsteady_pressure,thread,position)
{
cell_t f;
real preasure;
real velocity;
begin_c_loop(f,thread)
{
real time_step=RP_Get_Real("physical-time-step");//获取计算的步长
preasure=C_P(f,thread);//读取压力值
velocity=C_U(f,thread);//读取当前速度
F_PROFILE(f, thread, position)=preasure;//想将度读取的入口压力复制回去，判断UDF能否正常工作
//F_PROFILE(c, thread, position)=preasure-(182.6*time_step-0.0002)*velocity/7.6;
}
end_c_loop(f,thread)
}
出现的问题：把此UDF作为出口条件，仿真后和定常流结果一样，当作为入口条件后，原来的入口压力被变小，与出口压力相当，这是为什么？应该也是定常流才对啊，请大侠帮帮忙，我的问题出在哪里？

cora

preasure改成入口边界条件的数值

如烟

本人想仿真一个类似注射器的非定常流体运动，遇到两个问题，一个是压力入口的压力随弹簧的放松而减小，需要编写UDF，另一个问题是，压力入口的边界是运动的，寻找高手帮我解决这两个问题，必有重谢！QQ:158445498

[ 本帖最后由 如烟 于 2008-10-9 21:47 编辑 ]