| 查看: 915 | 回复: 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)); } |
回复此楼» 猜你喜欢
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
不应该看fileCode
已经有8人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
-
重要来源：本周末出结果
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
好奇怪的filecode
已经有6人回复
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
