数据下载:fall.py
- 视频教程:4自由落体和Paraview颜色配置
其中,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手册中有一个类似的实例