| 查看: 1995 | 回复: 4 | ||
hellowqy铁虫 (小有名气)
|
[求助]
求助：我编写的UDF在一个时间步内被调用两次
|
|
给位前辈 小弟我编写的UDF在一个时间步内被调用两次，导致输出的数据重复了一次（如图），恳请大家帮忙，献上200金币，万分感谢！！ # include "udf.h" # include "mem.h" # include "dynamesh_tools.h" # include "math.h" DEFINE_CG_MOTION(ts,dt,vel,omega,time,dtime) { FILE *fp = fopen("fluid_forces_and_moments.txt","a"  ;Domain *d = Get_Domain(3); Thread *tf1 = Lookup_Thread(d, 3); Thread *t; cell_t c; face_t f; real CG[3]; real FORCE[3]; real MOMENT[3]; real cx,cy; real wcx=0.0; real wcy=0.0; real va_liq=0.0; real z[ND_ND]; if (!Data_Valid_P ()) return; omega[0]=0.0; omega[1]=0.0; vel[0]=1*time; vel[1]=0.0; thread_loop_c(t,d) /*loops over all cell threads in domain*/ { begin_c_loop(c, t) /*loops over cells in a cell thread */ { if(0<C_VOF(c,t)&&C_VOF(c,t)<=1.0) { va_liq = va_liq + C_VOLUME(c,t); C_CENTROID(z,c,t); wcx = wcx + z[0]*C_VOLUME(c,t); wcy = wcy + z[1]*C_VOLUME(c,t); } } end_c_loop(c, t) } cx=wcx/va_liq; cy=wcy/va_liq; CG[0]=cx; CG[1]=cy; Compute_Force_And_Moment(d, tf1, CG, FORCE, MOMENT,FALSE); fprintf(fp,"%g,%g,%g,%g,%g\n",time,FORCE[0],FORCE[1],MOMENT[2],CG[0]); fclose(fp); }  00000000000.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
2楼2015-11-12 12:28:51
3楼2015-11-12 12:30:27
hellowqy
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 682
- 散金: 50
- 帖子: 64
- 在线: 61.1小时
- 虫号: 1433904
- 注册: 2011-10-09
- 性别: GG
- 专业: 机械动力学
4楼2015-11-17 10:57:48
kokal
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 1279.6
- 散金: 332
- 红花: 3
- 帖子: 335
- 在线: 244.3小时
- 虫号: 1058733
- 注册: 2010-07-16
- 专业: 流体力学
5楼2016-01-08 17:15:01
