| 查看: 969 | 回复: 2 | ||
[求助]
计算弯管的磨损量，利用fluent帮助里的程序，出错了，各位帮忙看看怎么回事 已有1人参与
|
|
程序代码是： #include "udf.h" #define MIN_IMPACT_VELO -1000 Domain *domain; enum { NUM_OF_HITS, AVG_DIAMETER, AVG_RADI_VELO, NUM_OF_USED_UDM }; int UDM_checked = 0; void reset_UDM_s(void); int check_for_UDM(void) { Thread *t; if (UDM_checked) return UDM_checked; thread_loop_c(t,domain) { if (FLUID_THREAD_P(t)) if (NULLP(THREAD_STORAGE(t,SV_UDM_I))) return 0; } UDM_checked = 1; reset_UDM_s(); return UDM_checked; } void reset_UDM_s(void) { Thread *t; cell_t c; face_t f; int i; if (!check_for_UDM()) return; Message("Resetting User Defined Memory...\n"  ;thread_loop_f(t, domain) { if (NNULLP(THREAD_STORAGE(t,SV_UDM_I))) { begin_f_loop(f,t) { for (i = 0; i < NUM_OF_USED_UDM; i++) F_UDMI(f,t,i) = 0; } end_f_loop(f, t) } else { Message("Skipping FACE thread no. %d..\n", THREAD_ID(t)); } } thread_loop_c(t,domain) { if (NNULLP(THREAD_STORAGE(t,SV_UDM_I))) { begin_c_loop(c,t) { for (i = 0; i < NUM_OF_USED_UDM; i++) C_UDMI(c,t,i) = 0; } end_c_loop(c,t) } else { Message(" Skipping CELL thread no. %d..\n", THREAD_ID(t)); } } Message(" --- Done.\n" ;} DEFINE_DPM_EROSION(dpm_accr, p, t, f, normal, alpha, Vmag, Mdot) { real A[ND_ND], area; int num_in_data; Thread *t0; cell_t c0; real imp_vel[3], vel_ortho; #if RP_2D if (rp_axi) { real radi_pos[3], radius; radi_pos[0] = P_POS(p)[1]; radi_pos[1] = P_POS(p)[2]; radius = NV_MAG(radi_pos); NV_VS(radi_pos, =, radi_pos, /, radius); imp_vel[0] = P_VEL(p)[0]; imp_vel[1] = NVD_DOT(radi_pos, P_VEL(p)[1], P_VEL(p)[2], 0); } else #endif NV_V(imp_vel, =, P_VEL(p)); vel_ortho = NV_DOT(imp_vel, normal); if (vel_ortho < MIN_IMPACT_VELO) return; if (!UDM_checked) if (!check_for_UDM()) return; c0 = F_C0(f,t); t0 = THREAD_T0(t); F_AREA(A,f,t); area = NV_MAG(A); F_STORAGE_R(f,t,SV_DPMS_ACCRETION) += Mdot / area; MARK_PARTICLE(p, P_FL_REMOVED); if (THREAD_TYPE(t) == THREAD_F_JUMP) return; num_in_data = F_UDMI(f,t,NUM_OF_HITS); F_UDMI(f,t,AVG_DIAMETER) = (P_DIAM(p) + num_in_data * F_UDMI(f,t,AVG_DIAMETER)) / (num_in_data + 1); C_UDMI(c0,t0,AVG_DIAMETER) = F_UDMI(f,t,AVG_DIAMETER); F_UDMI(f,t,AVG_RADI_VELO) = (vel_ortho + num_in_data * F_UDMI(f,t,AVG_RADI_VELO)) / (num_in_data + 1); C_UDMI(c0,t0,AVG_RADI_VELO) = F_UDMI(f,t,AVG_RADI_VELO); F_UDMI(f, t, NUM_OF_HITS) = num_in_data + 1; C_UDMI(c0,t0,NUM_OF_HITS) = num_in_data + 1; } DEFINE_ON_DEMAND(reset_UDM) { domain = Get_Domain(1); reset_UDM_s(); } 各位看看哪错了，应该怎样解决  错误.png |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在fluent中初始化的时候出错，请高手看看我的编程哪里有错
已经有8人回复
- fluent运用UDF三维弯管模拟， 如何给定进口圆周速度（初学者）
已经有6人回复
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
moons2
金虫 (正式写手)
- 应助: 144 (高中生)
- 金币: 1440.5
- 红花: 8
- 帖子: 431
- 在线: 239.4小时
- 虫号: 3784607
- 注册: 2015-04-03
- 性别: GG
- 专业: 流体力学
2楼2015-04-23 16:03:31
3楼2015-04-24 15:45:25
