| 查看: 1951 | 回复: 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 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
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
