请教各位前辈: 关于球和球间接触力的问题

smalldds

请问各位前辈
我今天设一个试体在顶部以墙来设一个-y 方向的速度来向下压
我想看球和球之间的接触力时
我打print information contact force可以显示出所有球和球的接触力
但是数据太多的话command window会看不到跑到上面的数值
我试着把球和球的接触力抓到.txt or .dat中
假设我要看球(ID 184)和球(ID206)的接触力
hist id=1 ball id 206 cforce id 184  or  hist id=1 ball cforce id 206 id 184
但是都失败了不知道为什么会失败
还请各位前辈能否指导一下该怎么样才能把数值写出来
感谢

ahzhuguiqing

能用history这样直接提取两个球的接触力么？？？可以用较笨的方法得到指定球的接触力：遍历其中一个球（184）的所有接触，当与其接触的球的id为206时，将法向力与切向力分别赋值给变量，用history记录这两个变量的值即可。付上程序。

def nfor
  bp=find_ball(184)
  cp=b_clist(bp)
  nfor=0
  sfor=0
  loop while cp # null
    nfor0=c_nforce(cp)
    sfor0=c_sforce(cp)
    if c_ball1(cp)=bp then
      bp_other=c_ball2(cp)
      cp=c_b1clist(cp)
    else
      bp_other=c_ball1(cp)
      cp=c_b2clist(cp)
    end_if
    id_other=b_id(bp_other)
    if id_other = 206 then
      nfor=nfor0
      sfor=sfor0
    end_if
  end_loop
end
history id 121 nfor
history id 122 sfor