| 查看: 4210 | 回复: 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); } 另外在画网格的时候，画这个受力的物体有没有什么特殊的要求，比如边界条件上等等。  救命啊~~~ |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有没有好心人告诉个管理员账号或者明天帮忙查询一下啊
已经有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
