这里是一个 同构造剥蚀 计算实例。
在终端或 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 |