24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

你好CPC

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 62
帖子: 25
在线: 29.9小时
虫号: 2608914
注册: 2013-08-21
性别: GG
专业: 工程热力学

[求助] 求各位大侠看下这条体积力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));
}

回复此楼

» 猜你喜欢

1楼 2016-03-23 20:50:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

小七工作室1

新虫 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
月只蓝: 金币-100, 应助指数-2, 屏蔽内容, 违规存档, 多次用完全一样、无实质性帮助的内容，回复不同应助帖！ 2017-03-24 10:58:08

本帖内容被屏蔽