易贡滑坡数值模拟

1.案例介绍

本计算实例采用Massflow软件模拟易贡滑坡的运动。计算模型见下图1。所有边界条件为默认设置，即墙边界条件；基底模型采用库伦摩擦模型，相应库伦摩擦系数设为0.215。土压力系数设为常数1。滑坡地形区域如下图2所示。

2.操作步骤

第一步：新建一个文件夹，命名为caseX，将程序执行文件Massflow.exe拷贝到该文件夹目录下。
第二步：进行计算前的一些设置准备。包括基本参数设置、模型导入以及分析步设置等。本计算案例是地形高程信息是采用文件读取的方式进行导入的。因此将地形高程信息文件z.txt拷贝到工作文件目录下。针对计算案例的初始流体高度，用户可在dll_custom文件夹下自行修改相应的代码，通过Microsoft Visual Studio（VS）编译的动态链接库dll_custom.dll文件拷到当前工作文件夹目录下。通过用户自定义构建模型需要在create_h_custom.f90进行的设置如下：
do i=1-mbc,nx+mbc
do j=1-mbc,ny+mbc
if(x(i,j)>=691500.and.x(i,j)<=692000.and.y(i,j)>=3.346e+06.and.y(i,j)<=3.3465e+06)then
h(i,j)=150
else
h(i,j)=h_dry
endif
enddo
enddo
该案例的边界采用默认，即墙边界。其他模型参数直接在命令流文件input.txt中输入。
第三步：所有参数设置完成后，双击执行文件Massflow.exe运行。

3.本计算实例的命令流文件—input.txt

命令流1：通过读取文件的方式设置地形和流体高度
*type,1
*erosion,0
*density,0
*nregions,1
*ndim,2
*mbc,0
*remesh,0
*parallel,0
*create,grid,file,dem,all,z.txt
*create,topo,file,dem,all,z.txt
*create,h,custom,all
*create,bound,custom,all
*material,1,density,1000
*material,1,friction,1,0.215
*material_assign,1,0,1
*time,1,20,0.05,0.25
*output,1,tecplot,1,all,1,4