| 查看: 836 | 回复: 2 | |||
[求助]
求各位大侠看下这条体积力UDF的问题。 已有1人参与
|
|
要求：在TSTART之前，按照颗粒位置得到一个电场E0，然后将这个E0值导入到TSTART后的bforce中去。 计算结果是，bforce好像就是0了，即E0=0.0。 测试了下，<=TSTART时的bforce用C1*E0*P_DIAM(p)*P_DIAM(p)*(C2*pow(R,C3)+C4)来计算的话，体积力还是作用了的。 所以TSTART之后用的E0值不是之前计算的值，而是定义变量E0的初始赋值，请问要如何修改才能实现要求。 #define C1 0.0000000000383 #define C2 2361.0 #define C3 -1.0341 #define C4 44884.0 #define TSTART 0.1 DEFINE_DPM_BODY_FORCE(particle_body_force, p, i) { real bforce; real E0=0.0; cell_t c=RP_CELL(&(p->cCell)); Thread *t=RP_THREAD(&(p->cCell)) if(P_TIME(p)<=TSTART) { E0=C2*pow(P_POS(p)[0],C3)+C4; bforce=0.0; } else { if(i==1) bforce=0.0; else if(i==0) bforce=C1*E0*P_DIAM(p)*P_DIAM(p)*(C2*pow(R,C3)+C4); } return(bforce/P_MASS(p)); } |
回复此楼» 猜你喜欢
318求调剂
已经有8人回复
-
求调剂
已经有7人回复
-
0703化学
已经有13人回复
-
307分材料专业求调剂
已经有4人回复
-
0855求调剂材料
已经有8人回复
-
332求调剂
已经有9人回复
-
甲醇/二氯 1:15过柱子
已经有3人回复
-
272求调剂
已经有4人回复
-
一志愿西北农林畜牧专硕336分求调剂
已经有3人回复
-
求调剂
已经有10人回复
2楼2016-11-03 15:02:54
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
月只蓝: 金币-100, 应助指数-2, 屏蔽内容, 违规存档, 多次用完全一样、无实质性帮助的内容，回复不同应助帖！ 2017-03-24 10:58:08
月只蓝: 金币-100, 应助指数-2, 屏蔽内容, 违规存档, 多次用完全一样、无实质性帮助的内容，回复不同应助帖！ 2017-03-24 10:58:08
|
本帖内容被屏蔽 |
3楼2016-11-03 17:27:02
10