|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
请各位高手看看错在哪里,我现在在做的课题需要用到动网格,比较急,非常感谢!希望你们能帮忙看看!QQ:364244643 Email: llj191@sina.com
#include "udf.h"
#include "dynamesh_tools.h"
#define L (90)
#define A (5)
#define FRE (10)
#define OMEGA (2.0*M_PI*FRE)
DEFINE_GRID_MOTION (tube, domain, dt, time, dtime)
{
Thread *tf = DT_THREAD (( Dynamic_Thread * ) dt );
face_t f;
Node *v;
real disp;
real x;
int n;
SET_DEFORMING_THREAD_FLAG(tf->t0);
disp = OMEGA*A/L*cos(OMEGA*time)* dtime;
begin_f_loop (f, tf )
{
f_node_loop (f, tf, n)
{
v=F_NODE (f, tf, n);
if(NODE_MARK(v)==1)
{
NODE_MARK(v) = 0;
x=NODE_X(v);
NODE_Y(v) += disp*x*(0.6-x);
}
}
Update_Face_Metrics(f,tf);
}
end_f_loop(f,tf);
} |
|