[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
初始网格如第一个图:上下边界均为固体铜,左右为自由边界,原点在下面边界的中点,现准备让左边边界按照x=x+y/10000(y不变),右边边界按照x=x-y/10000(y不变)运动,但边界动一格后就出现负体积,由于我对DEFINE_GRID_MOTION理解还很粗浅,恳请大家帮忙。[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
我的源程序是:#include "udf.h"
DEFINE_GRID_MOTION(shiyan_l,domain,dt,time,dtime)
{
real x,p,q;
face_t f;
Node *v;
int n;
Thread *tf=DT_THREAD(dt);
Thread *t0=THREAD_T0(tf);
SET_DEFORMING_THREAD_FLAG(t0);
begin_f_loop(f,tf)
{
f_node_loop(f,tf,n)
{
v=F_NODE(f,tf,n);
if(NODE_POS_NEED_UPDATE(v))
{
NODE_POS_UPDATED(v);
p=NODE_X(v);
q=NODE_Y(v);
p=p+q/10000;
NV_D(NODE_COORD(v),=,p,q,0);
}
}
}end_f_loop(f,tf)
}
DEFINE_GRID_MOTION(shiyan_r,domain,dt,time,dtime)
{
real x,p,q;
face_t f;
Node *v;
int n;
Thread *tf=DT_THREAD(dt);
Thread *t0=THREAD_T0(tf);
SET_DEFORMING_THREAD_FLAG(t0);
begin_f_loop(f,tf)
{
f_node_loop(f,tf,n)
{
v=F_NODE(f,tf,n);
if(NODE_POS_NEED_UPDATE(v))
{
NODE_POS_UPDATED(v);
p=NODE_X(v);
q=NODE_Y(v);
p=p-q/10000;
NV_D(NODE_COORD(v),=,p,q,0);
}
}
}end_f_loop(f,tf)
}
[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
我也遇到类似问题,请高手帮忙解决下[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
你的动网格重新划分方法选的什么?参数是否合理?[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
我采用的是layering,那个collapse factor我从0.01到0.8都试过了,结果是一样的,而且运行的步数都是一样的。后来我又试过smoothing和remeshing,参数也试过好多组不同组合,但结果和以前的都一样,所以我就怀疑是不是那个宏用错了?高手帮忙呀!文字文字[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
还有其他参数,比如时间步,最大最小单元SIZE,还有上面的边要设置为变形……动网格说简单很简单说复杂也复杂
[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
顶起[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
友情顶[原创]关于动网格DEFINE_GRID_MOTION的算例,请大家帮忙!!!
将网格换成非结构网格是否会好些 图呢 淘宝店铺酷酷老哥 fluent 视频教程15天的 培训 视频。具体的看一下, 那里有手机号,店铺里有联系方式发短信 联系现在的联系方式是qq1067620876,手机在qq资料里,网址是 http://sell.taobao.com/auction/goods/goods_on_sale.htm?tracelog=tszbb 你的单位确定调好了吗?你看一个时间步下网格坐标变化多少,你这用cg_motion都可以啊
页:
[1]