| 查看: 2321 | 回复: 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坐标就不变.求大神指点啊！！小女子不胜感激！ |
回复此楼» 猜你喜欢
298求调剂
已经有5人回复
-
272求调剂
已经有4人回复
-
285求调剂
已经有6人回复
-
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
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
4楼2018-03-07 19:44:54
