download:fall.py

其中,fall.py 中内容如下:

######################################
# title: 自由落体
# date: 2017-04-28
# authors: 李长圣
# E-mail: sheng0619@163.com
# more info, see www.geovbox.com
#######################################
START ; 开始计算
BALL  ID=0  X=2.0  Y=0.5  RAD=0.5  COLOR=blue ;生成一个蓝色颗粒,圆心(2.0,0.5),半径0.5
BALL  ID=1  X=2.0  Y=2.5  RAD=0.5  COLOR=red
PROP  DEN=1e3  FRIC=0.0  KN=10e6  KS=10e6  DAMP=0.1;密度DEM,摩擦系数FRIC,法向刚度KN,切向刚度KS,阻尼DAMP
FIX   x y spin      RANGE id=(0, 0) ;固定id=0的颗粒的x和y方向的运动,并且固定其角速度spin
SET   DT=2.5e-04  GRAVTIY=(0.0, -9.81) ; 时间步DT,GRAVITY重力加速度
SET   vtk 100 ; 设置每100步保存一次vtk格式(使用Paraview打开)的计算结果
CYC   20000 ; 计算20000步
STOP ; 结束计算

程序必须以START、RESTORE和LOAD开始一个计算,上面的实例是以start开始的,生成了两个颗粒,固定蓝色颗粒,红色颗粒作自由落体运动。

  • 动态图如下:

















软件中所涉及的参数的单位均采用国际制标准单位(即千克kg,米m,秒s,牛N,帕Pa等)。

  • YADE手册中有一个类似的实例