| 查看: 835 | 回复: 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)); } |
回复此楼» 猜你喜欢
材料与化工371求调剂
已经有5人回复
-
328分调剂
已经有4人回复
-
295求调剂
已经有7人回复
-
081700，311，求调剂
已经有13人回复
-
283分求调剂
已经有6人回复
-
0854求调剂
已经有3人回复
-
332求调剂
已经有4人回复
-
284求调剂
已经有5人回复
-
304求调剂（085602，过四级，一志愿985）
已经有11人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有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
