UDF的求助

一个初始形状为规则矩形的管道，Y向高度为0.2，X向长度为0.5，Z垂直纸向面里，长度为1。

在水流过程中左右两侧的管道壁面会发生类似图中绿线所示的变形，在水流过程中，两侧壁会做随Y向位置变化的运动（上面位移变化小，下面位移变化大），其中黑色表示的位移变化。此外，在不同的时间中沟壁也会发生收缩或扩张现象，譬如说在t<1秒时，沟壁成收缩状态，变化至绿线所示位置；在1<t<2秒时，从收缩状态有恢复到初始红线位置；在2<t以后，管道壁面不再发生位移变化，保持不变。

下面是我自己凑写的UDF，仅仅作用在左侧壁面。右侧壁面想着是修改下速度符号即可。

还请大家帮忙指点下，谢谢了。

#include"udf.h"
DEFINE_CG_MOTION(value, dt, cg_vel, omega, time,dtime); x7 F* C8 Z1 {+ q. V( P
{
Thread *t;$ R& j: Q/ w1 S
real x[ND,ND,ND];

real y;
, face_t f;

begin_f_loop(f, t)
{
F_CENTROID (x, f, t);) H% x2 ^3 C/ i6 u$ X5 k e* p
y=x[0];

if(time>=0.00&time<=0.2)

cg_vel(f,t,position)=2*(0.2-y);

elseif (time>=0.2%time<=0.4)

cg_vel(f,t,position)=-2*(0.2-y);

if(time>=0.4)

cg_vel(f,t,position)=0;

}2 ]2 W, X: s5@
end_f_loop (f, t)! a"

}