| 查看: 2275 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
[求助]
动网格边界面结点为什么不动呢,附程序,急求助!
|
||
|
#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坐标就不变.求大神指点啊!!小女子不胜感激! |
» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有9人回复
最失望的一年
已经有17人回复
为什么nbs上溴 没有产物点出现呢
已经有6人回复
求推荐博导
已经有4人回复
存款400万可以在学校里躺平吗
已经有34人回复
求助一下有机合成大神
已经有4人回复
求推荐英文EI期刊
已经有5人回复
26申博
已经有3人回复
基金委咋了?2026年的指南还没有出来?
已经有10人回复
疑惑?
已经有5人回复
3楼2016-03-16 17:56:25
18428346586
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 122.3
- 散金: 20
- 帖子: 29
- 在线: 25.3小时
- 虫号: 3226490
- 注册: 2014-05-22
- 专业: 多相流热物理学
2楼2015-12-02 20:49:09
4楼2018-03-07 19:44:54













回复此楼