[原创]关于动网格DEFINE_GRID_MOTION的算例，请大家帮忙！！！

dabaiyang

[watermark]初始网格如第一个图：上下边界均为固体铜，左右为自由边界，原点在下面边界的中点，现准备让左边边界按照x=x+y/10000（y不变），右边边界按照x=x-y/10000(y不变)运动，但边界动一格后就出现负体积，由于我对DEFINE_GRID_MOTION理解还很粗浅，恳请大家帮忙。

dabaiyang

我的源程序是：
#include "udf.h"
DEFINE_GRID_MOTION(shiyan_l,domain,dt,time,dtime)
{
real x[ND_ND],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[ND_ND],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)
}

spacecxh

我也遇到类似问题，请高手帮忙解决下

gallongallon

你的动网格重新划分方法选的什么？参数是否合理？

dabaiyang

    我采用的是layering,那个collapse factor我从0.01到0.8都试过了，结果是一样的，而且运行的步数都是一样的。后来我又试过smoothing和remeshing，参数也试过好多组不同组合，但结果和以前的都一样，所以我就怀疑是不是那个宏用错了？高手帮忙呀！文字[color=#000000]文字

gallongallon

还有其他参数，比如时间步，最大最小单元SIZE，还有上面的边要设置为变形……
动网格说简单很简单说复杂也复杂

css

顶起

luckydog

友情顶

magician

将网格换成非结构网格是否会好些

fangdian

图呢

huhaijunhit

淘宝店铺  酷酷老哥 fluent 视频教程15天的 培训 视频。具体的看一下, 那里有手机号  ，店铺里有联系方式  发短信 联系
现在的联系方式是  qq1067620876，手机在qq资料里，网址是 http://sell.taobao.com/auction/goods/goods_on_sale.htm?tracelog=tszbb

打俩鸡蛋

你的单位确定调好了吗？你看一个时间步下网格坐标变化多少，你这用cg_motion都可以啊