这里是一个 同构造剥蚀 计算实例。

在终端或 xshell 中运行 vboxdaily syn_erosion.py

syn_erosion.py.py 中内容如下:

######################################
# title: 同构造沉积
# date: 2019-01-19
# authors: 李长圣
# E-mail: sheng0619@163.com
# more info, see www.geovbox.com
#######################################
#从 0000046000.sav <版本1.3> 计算节点恢复,由 一个实例学会VBOX 生成
RES 0000046000.sav
# 每次 100 步更新一次进度条
SET stepbar 100
#设置墙的挤压速度 x方向速度为2.0
WALL id 1 xv 2.0
#设置墙的挤压量x方向推进1000.0,每挤压1000.0保存一次计算结果
IMPLE wall id 1 xmove 1000.0 save 1000.0 print 1000.0 ps 1000.0

##################################### 剥蚀 #####################################
#删除4000米以上的颗粒
DEL RANGE y 4000.0 999000.0
################################################################################

#设置墙的挤压速度 x方向速度为2.0
WALL id 1 xv 2.0
#设置墙的挤压量x方向推进5000.0,每挤压2000.0保存一次计算结果
IMPLE wall id 1 xmove 5000.0 save 2000.0 print 1000.0 ps 1000.0

#计算停止
STOP

得到的演化过程如下:

演化过程

单位 (km)

构造形态及平均应力、最大剪切应力

微观参数表

颗粒直径         
d (m)
颗粒密度         
ρ (kg∙m-3)
颗粒摩擦系数         
缩短速率         
υ (m∙s-1)
时间步         
s
60,80 2500 0.3 2.0 0.05