关于PFC2D的contact生成的问题

lossxie

不知道大家有没有研究过PFC的contact的生成的问题，我用了一个很简单的模型去验证了一下，理论手册中所提的contact搜索及生成的方法，发现一个奇怪的问题。因为在手册中提过(包括cundall的论文)，生成contact首先同cell检索相邻ball，然后再计算两个ball之间的gap如果小于ctol就生成contact，但是我生成了三个ball，却让其中一个球(2号球)，从不同的角度接近另一个球（3号球），在生成contact瞬间通过fishcall函数计算两个ball之间的gap，发现2号球从不同角度接近3号球时，gap是不同的，这让我百思不得其解，而且还发现PFC2D程序根本不认函数cell_object，这是fish手册上的函数哦。请高手指点一下。
new
set random
set disk on
set max_balls=1
def Para
  bp1_num=3                    ;固定的ball的IDnumber
  bp2_num=2                    ;移动的ball的IDnumber
  xcor_value=15
  ycor_value=10
  vel_value=1e-5
  bp1=find_ball(bp1_num)
  bp2=find_ball(bp2_num)
  xvel_value=-1*vel_value*xcor_value/ycor_value
  yvel_value=vel_value*ycor_value/ycor_value  
end
Para
ball id=1 x=0.0 y=0.0 rad=0.5
ball id=2 x=xcor_value y=0.0 rad=0.5
ball id=3 x=0 y=ycor_value rad=0.5
plot ball Id=on
plot add cell Id=on
property dens=2000 kn=1e8 ks=1e8
ini xvel=0 yvel=0 spin=0
;ini xvel=-1e-4 yvel=0 spin=0 range ID=2
ini xvel=xvel_value yvel=yvel_value spin=0 range ID=2
fix x y spin
set dt dscale
set hist_rep=1
def num_contact
  num_contact=ccount
end
def Getball2_inf
  while_step
  ball2_xcor=b_x(bp2)
  ball2_ycor=b_y(bp2)
end
def GetCell_Tollinf
  Cellxlow=cell_xlow
  Cellylow=cell_ylow
  Cellxnum=cell_xnum
  Cellynum=cell_ynum
  Cellsize=cell_size
  command
    print Cellxlow
    print Cellylow
    print Cellxnum
    print Cellynum
    print Cellsize
  endcommand
end
def Creat_contact
  cp=fc_arg(0)
  bp2_xcor=b_x(bp2)
  bp2_ycor=b_y(bp2)
  dist=sqrt((b_x(bp2)-b_x(bp1))^2+(b_y(bp2)-b_y(bp1))^2)
  gap=dist-b_rad(bp1)-b_rad(bp2)
  ratio=2*gap/(b_rad(bp1)+b_rad(bp2))
  command
    print dist gap ratio
    pause
  endcommand
end
print num_contact
;set fishcall FC_CONT_CREATE Creat_contact   ;设定fishcall为contact生成
set fishcall 6 Creat_contact                 ;设定fishcall为contact生成
hist ball2_xcor
hist ball2_ycor
cycle 200000
请高手给点意见。

lossxie

还有就是当你使用print cell 时cell-size不知道是什么意思，根本不是cell的面积，而且还有两个量一个是actual cell 和 total cell 两个还不相等。这究竟是为什么呢？