| 查看: 1069 | 回复: 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 |
回复此楼» 猜你喜欢
2027广东省杰青
已经有5人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在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
