动力分析中，颗粒碰撞逃离问题

dragonclouder

做了一个简单的例子，三个小球，参看所附程序，当中间小球速度大于26时，小球不会与前面进行阻挡的小球产生碰撞，而是直接逃离到远端去了，这是什么原因呢？我把用于阻挡的小球换成wall，逃离的速度会更小，这是为什么呢？
new
set pint 20  ; update plots every this many cycles
ball id 1 rad 0.6 x 50 y 0.5 z 0.5
ball id 2 rad 0.6 x 70 y 0.5 z 0.5
ball id 3 rad 0.6 x 90 y 0.5 z 0.5
fix  x y z xspin yspin zspin range id 1
;fix  x y z xspin yspin zspin range id 3
ini xvel=-26 range id 2
plot create footing
plot add ball
plot show
prop kn 1e8 ks 1e8 fric 0.2 density 2e5
prop n_bond 1e5 s_bond 1e5
trace energy on                             ;monitor the energy curve
his energy body
his energy bond
his energy strain
his energy kinetic
hist diag muf
hist diag mcf
his ball xpos id=2
his ball zpos id=2
his ball yspin id=2
his ball xvel  id=2
plot create view_1
plot set back white
plot add his 1 2 3 4
plot set title text &
'History of measured energies Eb, Epb, Ec and Ek (file Energy_1.DAT)'

cycle 10000
;save hopp_3a.sav
plot show

sddzyc2

时步太大了 这个问题也曾困扰了我好久

scvictory

你的问题原因有二
(1)动力问题，你没有设置时间步。
(2)第二，你没设置阻尼，pfc默认阻尼是0.7，也就是下一步的速度会是无阻尼的0.3倍，所以你的颗粒速度快速降低了。
你可以用自由落体尝试，如果设置加速度为10，默认阻尼得到的加速度只有3.
我费了2周才找到这个问题所在。

程序设计

3# scvictory
请问设置时间步大小 依据什么啊  大神