| 查看: 2415 | 回复: 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坐标就不变.求大神指点啊！！小女子不胜感激！ |
回复此楼» 猜你喜欢
国自然评审意见
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有12人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
4楼2018-03-07 19:44:54
18428346586
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 122.3
- 散金: 20
- 帖子: 29
- 在线: 25.3小时
- 虫号: 3226490
- 注册: 2014-05-22
- 专业: 多相流热物理学
2楼2015-12-02 20:49:09
3楼2016-03-16 17:56:25
30