| 查看: 930 | 回复: 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)); } |
回复此楼» 猜你喜欢
面上没中，邀请各位路过的虫友分析一下分数
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有17人回复
-
Ei源刊怎么投
已经有4人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
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