| 查看: 3984 | 回复: 2 | |||
[求助]
关于Lookup_Thread（domain，id）宏求教各位前辈已有1人参与
|
|
我做的是二维水池中，使用define cg motion 使物体在漂在水面上做振动，要读取在其运动的过程中受到水的作用力，使用Lookup_Thread宏搜寻物体的壁面指针，在id中写入物体壁面id之后一运行就会出错，求各位前辈帮帮忙，金币不多，但情况紧急，不胜感激。 #include "udf.h" #include "mem.h" #include "dynamesh_tools.h" DEFINE_CG_MOTION(tuoyuan,dt,vel,omega,time,dtime) { face_t f; cell_t c; real a0=0.01; real w0=6.2832; real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND]; FILE *fp; Thread *t=DT_THREAD(dt); Domain *domain=Get_Domain(1); Thread *tf1=Lookup_Thread(domain,5); /*id 5为运动物体的壁面id*/ int i; NV_S(vel,=,0.0); NV_S(omega,=,0.0); if(!Data_Valid_P()) return; for(i=0;i<=ND_ND;i++) { f_glob=0; m_glob=0; } for(i=0;i<ND_ND;i++) x_cg=DT_CG(dt); vel[0]=a0*w0*cos(w0*time); Compute_Force_And_Moment(domain,tf1,x_cg,f_glob,m_glob,TRUE); fp=fopen("F:\\\\FLUENT\\\\zd1.txt","a+"  ;fprintf(fp,"%f,%.10lf,%.10lf,%.1f,%.1f,%f,%f,%f\\n",time,x_cg[0],x_cg[1],f_glob[0],f_glob[1],m_glob[0],m_glob[1],m_glob[2]); fclose(fp); } 另外在画网格的时候，画这个受力的物体有没有什么特殊的要求，比如边界条件上等等。  救命啊~~~ |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
有没有好心人告诉个管理员账号或者明天帮忙查询一下啊
已经有15人回复- 请大侠们帮我分析一下我的评审结果是几A几B呀？谢谢啦
已经有12人回复
- 一个UDF并行报错
已经有17人回复
- udf 重力方向变化
已经有17人回复
- 求助--投过science of the total environment 的前辈帮忙看一下，谢谢
已经有7人回复
小原子1224
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 83.1
- 帖子: 13
- 在线: 3.1小时
- 虫号: 2605251
- 注册: 2013-08-19
- 性别: MM
- 专业: 工程热物理相关交叉领域
2楼2014-10-09 20:19:23
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
3楼2014-10-10 09:21:39