求助：PFC3D中如何循环删除指定区域的小球

civilchenxi

求助：PFC3D中，在重力作用下小球不断运动到某一区域，如何循环删除这一区域的小球，使得上面的小球不断下落？

rock_cumt

试一下这个：
;y as vertical direction
set y_bottom=..
def del_b
   bp = ball_head
  loop while bp # null
    bpnext = b_next(bp)
    y_pos = b_rad(bp)+b_y(bp)
    if y_pos<y_bottom then
      oo = b_delete(bp)
    end_if
    bp = bpnext
  end_loop
end

set fishcall  1 del_b


Good luck!  

civilchenxi

rock_cumt 发表于 2012-9-11 22:47
试一下这个：
;y as vertical direction
set y_bottom=..

我x与z方向也要限制，就是当小球自由下落到某长方体区域时，循环删除这一区域的小球

rock_cumt

定义区域
range name void_area x= (x1, x2) y= (y1, y2) 三维可以加上 z= (z1,z2)
循环删除
  bp = ball_head
  loop while bp # null
    bpnext = b_next(bp)
    if inrange( 'void_area', bp ) = 1 then
      oo = b_delete(bp)
    end_if
    bp = bpnext
  end_loop

civilchenxi

rock_cumt 发表于 2012-9-12 09:23
定义区域
range name void_area x= (x1, x2) y= (y1, y2) 三维可以加上 z= (z1,z2)
循环删除

谢谢无私的帮助，我试试

civilchenxi

rock_cumt 发表于 2012-9-12 09:23
定义区域
range name void_area x= (x1, x2) y= (y1, y2) 三维可以加上 z= (z1,z2)
循环删除

res dan.save
group rock range x 0 8 y 0 15 z 30 50
group ore range x 0 8 y 0 15 z 0 30
delet wall 13
delet wall 17
def del_b
command
  cyc 500
  delete ball range hole2
  delete ball range hole3
endcommand
end
del_b
删除wall 13 17 ,循环cyc500后，小球在自重下落入hole2与3，给定判断条件(hole2与hole3中rock与ore的个数比满足一定条件)若不满足删除hole2与3中的小球后再回到 cyc500，若满足则停止运行，不知如何实现？

civilchenxi

请教，求助

rock_cumt

def _cyc
    nstp = nstp + 1
    if nstp = 500 then
       判断条件
          满足{
                 command
                     delete ball range hole2
                     delete ball range hole3
                     pau
                 command
                }
        不满足{
              nstp = 0
        }
    endif
endif
end
这个过程的思路是每500判断一次，如果满足条件，运算暂定，否则进入下一个500步的循环
                 

civilchenxi

rock_cumt 发表于 2012-9-14 20:59
def _cyc
    nstp = nstp + 1
    if nstp = 500 then

您好，谢谢您的解答。一个区域（hole2或3）中如何统计两种不同小球的个数呢。group rock
group ore

rock_cumt

set ball extra 1
def set_extra
  bp = ball_head
  loop while bp # null
    bz = b_z(bp)
    if bz >=30 then
      b_extra(bp,1) =1 ;rock
   else
     b_extra(bp,1)=0; ore
   endif
  bp=b_next(bp)
endloop
end

def count
  bp = ball_head
  loop while bp # null
     ;if ball in hole2
        if (b_extra(bp,1) = 1 then
            h21 = h21 +1  ;for rock
       else h22 = h22 +1 ;for ore
   if ball in hole3
     ...........
  endloop
end

note： set fishcall should set after "res dan.sav"
then followed by @set_extra
function @counter can be invoke in @_cyc
For details, you should fix it yourself
good luck

civilchenxi

rock_cumt 发表于 2012-9-14 21:30
set ball extra 1
def set_extra
  bp = ball_head

有点思路了，谢谢了

stlyt

请问还在研究pfc吗