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