24小时热门版块排行榜

>论坛更新日志 (3315)
>虫友互识 (349)
>导师招生 (201)
>文献求助 (197)
>硕博家园 (107)
>公派出国 (90)
>考博 (85)
>找工作 (71)
>论文投稿 (69)
>博后之家 (68)
>基金申请 (64)
>论文道贺祈福 (60)
>考研 (60)
>休闲灌水 (60)
>教师之家 (51)
>招聘信息布告栏 (38)

返回列表

feir122

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1
帖子: 31
在线: 36.2小时
虫号: 65862
注册: 2005-04-30
性别: MM
专业: 工程地质学

[求助] 动网格边界面结点为什么不动呢，附程序，急求助！

#include "udf.h"
#include"dynamesh_tools.h"
DEFINE_GRID_MOTION(moving_wall, domain, dt, time, dtime)
{
  Thread *tf = DT_THREAD (dt);
  face_t f;
  Node *v;
  int n;
  cell_t c0;
  Thread *t0;
  real x,y,z;

  SET_DEFORMING_THREAD_FLAG (THREAD_T0 (tf));


  begin_f_loop (f, tf)
{

c0 = F_C0(f,tf);
t0 = THREAD_T0(tf);
C_UDMI(c0,t0,5)=C_UDSI_G(c0,t0,0)[0]+C_UDSI_G(c0,t0,1)[1];
F_UDMI(f,tf,5)=C_UDMI(c0,t0,5);
   f_node_loop (f, tf, n)
      {
      v = F_NODE (f, tf, n);
if (NODE_POS_NEED_UPDATE (v))
         {
{NODE_POS_UPDATED (v);

x=NODE_X(v);
y=NODE_Y(v);
z=NODE_Z(v);
z=z-C_UDMI(c0,t0,5)*dtime;
                     }
            }
      }
}
  end_f_loop (f, tf);
}

我是想控制底面边界上结点的z坐标下移一定的距离，下移的大小每个结点处都不太相同，取决于F——udmi(f,tf,5)中对应结点处值的大小，所以我就用了z=z-C_UDMI(c0,t0,5)*dtime或者z=z-F_UDMI(f,t,5)*dtime来定义，因为我觉得F_UDMI(f,t,5)中包含面上每个结点的值了。但是，动网格运行好久都不见变化，网格不动。是不是F_UDMI(f,t,5)不能直接这么使用给z坐标啊？在后处理中看UDM5都是一直在变化，一直计算的很正常，但是就是结点z坐标就不变.求大神指点啊！！小女子不胜感激！

回复此楼